• 콘크리트학회논문집
• /
• 제13권3호
• /
• pp.268-276
• /
• 2001

본 논문은 고온에서의 콘크리트 재료모델을 연구하였다. 콘크리트 응력-변형률 곡선은 온도가 증가함에 따라 그 형태가 변한다. 온도에 따른 콘크리트 재료거동의 변화를 나타내기 위하여 변형된 Saenz 제안식을 이용하여 응력-변형률 관계를 표시하였다. 고온에서의 급격한 변형률의 증가현상을 설명하기 위하여, 콘크리트의 변형률 성분을 순수 열팽창 변형률, 열적크리프 변형률, 과도 변형률 및 역학적 변형률로 구분하여 나타내었다. 열적크리프 변형률은 Baily-Norton의 장기크리프 곡선 식을 수정.제안하여 1축 실험 결과를 온도, 시간 및 응력의 함수로 표현하였고, 또한 유효응력 및 유효변형률 개념을 도입하여 다차원에서도 적용할 수 있는 모델을 제시하였다. 과도 변형률을 제안하여 다공탄성 거동을 가정한 콘크리트 내에 포함된 공극 및 수분의 작용을 역학적 거동의 영향을 분석하고자 하였다. 마지막으로, 본 논문에서 제시한 고온에서의 콘크리트 재료모델을 이용한 해석결과를 실제 화재실험자료와 비교하였다.

• 한국수자원학회논문집
• /
• 제32권4호
• /
• pp.417-426
• /
• 1999

횡유입이 없는 경우(sign곡선의 홍수파 유입)와 횡유입이 있는 경우(강우 유입)를 대상으로 선형과 비선형 Muskingum-Cunge법에 의한 단위폭 사면에서의 유출수문곡선의 특성을 비교·고찰하였다. 유출곡선은 선형법에서는 확산효과에 의해서 거의 대칭적으로 확산되나, 비선형법에서는 비선형효과와 확산효과의 상호작용을 의하여 상승부는 급하게 되고, 하강부는 완만하게 된다. 선형법은 유입된 질량을 정확히 보존하나, 비선형법은 질량의 증가나 손실을 초래한다. 비선형법에서 기저유량이 작고 유량변화의 비가 큰 경우, 완경사에서는 파형의 감쇠와 질량의 감소를 가져오나, 급경사에서는 kinematic shock의 생성과 함께 질량의 증가를 가져온다. 같은 특성을 갖는 사면의 경우, 유출곡선의 전파속도와 파형의 변화는 선형법에서는 단위폭당 기준유량에 좌우되나, 비선형법에서는 기저유량 및 유량 변화의 비(= 첨두유량/기저유량)에 좌우된다. 횡유입(강우)이 있는 비선형법의 경우 횡유입이 없는 경우와 마찬가지로, 사면경사가 완만하고 유량변화의 비가 클수록 첨두유량은 작아지고 파형은 확산되나, 질량의 손실은 거의 없다.

• 한국수자원학회:학술대회논문집
• /
• 한국수자원학회 2017년도 학술발표회
• /
• pp.281-281
• /
• 2017

하천에서 하도불안정(stream instability)으로 인하여 하상의 형태가 변화하고 하상파(sand wave)가 발생한다. 사련(ripple), 사구(dune) 등과 같은 하상파는 흐름저항을 유발하여 홍수시 수위를 증가시킨다. 수리실험 및 수치모형을 이용하여 사련 및 사구의 발달과정 그리고 이를 지나는 난류흐름에 대한 연구가 국외에서는 이루어지고 있지만 국내의 경우 거의 찾아보기 힘들다. 수치모형을 활용한 연구는 주로 횡방향으로 하상파가 일정하다는 가정하에 연직 2차원 수치모형을 적용하였으나 최근 컴퓨터 기술 및 수치기법의 고도화로 3차원 RANS(Reynolds averaged Navier-Stokes) 또는 LES(Large eddy simulation)를 이용한 수치모형이 개발되고 있다. 본 연구에서는 LES에 유사이송 및 하상변동 모형과 결합하여 사구발달에 대한 수치모의를 수행하였다. LES와 유사이송 및 하상변동 모형의 결합은 순간유속성분을 하상변동모형에 직접 적용되기 때문에 난류영향을 고려할 수 있는 것이 장점이다. 특히 사구의 발달에 따라 복잡한 흐름이 발생하며 3차원 와구조가 발생하므로 난류특성의 고려는 필수적이다. 수치모의는 Delft Hydraulics (Bakker et al., 1986)에서 수행한 수리실험 T39를 활용하였다. 수리실험은 길이 100 m, 폭 0.5 m 개수로에서 수행되었으며 평균유속은 0.611 m/s, 수심은 0.436 m이다. 하상파 실험에 사용된 유사입경은 0.78 mm 균일사를 사용하였다. 수치모의 조건은 수리실험과 동일하게 하였으나 계산시간의 효율을 고려하여 흐름방향의 계산영역은 4.0 m로 하고 주기경계조건(periodic boundary condition)을 부여하여 계산을 수행하였다. 수치모의 계산은 사구의 길이 및 파고가 평형상태에 이를 때까지 수행되었다. 수치모의 통해 사구발달에 따른 흐름 및 하상변동 특성을 분석하였다.

• Wind and Structures
• /
• 제28권5호
• /
• pp.285-298
• /
• 2019

For super long-span bridges, the aerodynamic forces induced by the flow passing the box girder should be considered carefully. And the Reynolds number sensitively of aerodynamic characteristics is one of considerable issue. In the study, a numerical study on the Reynolds number sensitivity of aerodynamic characteristic (flow pattern, pressure distribution and aerodynamic forces) of a twin-box girder were carried out using large eddy simulation (LES) with the dynamic Smagorinsky-Lilly subgrid model. The results show that the aerodynamic characteristics have strong correlation with the Reynolds number. At the leading edge, the flow experiences attachment, departure, and reattachment stages accompanying by the laminar transition into turbulence, causing pressure plateaus to form on the surface, and the pressure plateaus gradually shrinks. Around the gap, attributing that the flow experiences stages of laminar cavity flow, the wake with alternate shedding vortices, and turbulent cavity flow in sequence with an increase in the Reynolds number, the pressures around the gap vary greatly with the Reynold number. At the trailing edge, the pressure gradually recovers as the flow transits to turbulence (the flow undergoes wake instability, shear layer transition-reattachment station), In addition, at relative high Reynolds numbers, the drag force almost does not change, however, the lift force coefficient gradually decreases with an increase in Reynolds number.

• Nuclear Engineering and Technology
• /
• 제54권6호
• /
• pp.2276-2296
• /
• 2022

This study validates the applicability of the CUPID code for simulating subcooled wall boiling under high-pressure conditions against number of DEBORA tests. In addition, a new numerical technique in which the interfacial momentum non-drag forces are calculated at the cell faces rather than the center is presented. This method reduced the numerical instability often triggered by calculating these terms at the cell center. Simulation results showed good agreement against the experimental data except for the bubble sizes in the bulk. Thus, a new model to calculate the Sauter mean diameter is proposed. Next, the effect of the relationship between the bubble departure diameter (D_dep) and the nucleation site density (N) on the performance of the Wall Heat Flux Partitioning (WHFP) model is investigated. Three correlations for D_dep and two for N are grouped into six combinations. Results by the different combinations show that despite the significant difference in the calculated D_dep, most combinations reasonably predict vapor distribution and liquid temperature. Analysis of the axial propagations of wall boiling parameters shows that the N term stabilizes the inconsistences in D_dep values by following a behavior reflective of D_dep to keep the total energy balance. Moreover, ratio of the heat flux components vary widely along the flow depending on the combinations. These results suggest that separate validation of D_dep correlations may be insufficient since its performance relies on the accompanying N correlations.

• Smart Structures and Systems
• /
• 제28권6호
• /
• pp.745-760
• /
• 2021

Dynamic buckling of structure is one of the failure modes that needs to be considered since it may result in catastrophic failure of the structure in a short period of time. For a thin fiber-reinforced polymer (FRP) plate under compression, buckling is an inherent hazard which will be intensified by the existence of defects like holes, cracks, and delamination. On the other hand, the growth of the delamination is another prime concern for thin FRP plates. In the current paper, reinforcing the plates against buckling is realized by using SMA wires in the form of stitches. A numerical framework is proposed to simulate the dynamic instability emphasizing the effect of the SMA stitches in suppressing delamination growth. The suggested algorithm is more accurate than the other methods when considering the transformation point of the SMA wires and the modeling of the cohesive zone using simple and yet reliable technique. The computational design of the method by producing the line by line orders leads to a simple algorithm for simulating the super-elastic behavior. The Lagoudas constitutive model of the SMA material is implemented in the form of user material subroutines (VUMAT). The normal bilinear spring model is used to reproduce the cohesive zone behavior. The nonlinear finite element formulation is programmed into FORTRAN using the Newmark-beta numerical time-integration approach. The obtained results are compared with the results obtained by the finite element method using ABAQUS/Explicit solver. The obtained results by the proposed algorithm and those by ABAQUS are in good agreement.

• Geomechanics and Engineering
• /
• 제33권2호
• /
• pp.141-154
• /
• 2023

This study investigated the effect of water retention characteristics between aged and fresh Municipal Solid Waste (MSW) on the stability of the landfill. A series of transient numerical modeling for the slope of an MSW landfill was performed considering the variation of water retention characteristics due to leachate circulation. Four different scenarios were considered in this analysis depending on how to obtain hydraulic conductivity and the aging degree of materials. Unsaturated hydraulic properties of the MSW used for the modeling were evaluated through modified hanging column tests. Different water retention properties and various landfill conditions, such as subgrade stiffness, leachate injection frequency, and gas and leachate collection system, were considered to investigate the pore water distribution and slope stability. The stability analyses related to the factor of safety showed that unsaturated properties under those varied conditions significantly impacted the slope stability, where the factor of safety decreased, ranging between 9.4 and 22%. The aged materials resulted in a higher factor of safety than fresh materials; however, after 1000 days, the factor of safety decreased by around 10.6% due to pore pressure buildup. The analysis results indicated that using fresh materials yielded higher factor of safety values. The landfill subgrade was found to have a significant impact on the factor of safety, which resulted in an average of 34% lower factor of safety in soft subgrades. The results also revealed that a failed leachate collection system (e.g., clogging) could result in landfill failure (factor of safety < 1) after around 298 days, while the leachate recirculation frequency has no critical impact on stability. In addition, the accumulation of gas pressure within the waste body resulted in factor of safety reductions as high as 24%. It is essential to consider factors related to the unsaturated hydraulic properties in designing a landfill to prevent landfill instability.

• Nuclear Engineering and Technology
• /
• 제55권6호
• /
• pp.1974-1987
• /
• 2023

The load-following operation of small modular reactors (SMRs) requires accurate prediction of transient behaviors that can occur in the balance of plants (BOP) and the nuclear steam supply system (NSSS). However, 1-D thermal-hydraulics analysis codes developed for safety and performance analysis have conventionally excluded the BOP from the simulation by assuming ideal boundary conditions for the main steam and feed water (MS/FW) systems, i.e., an open loop. In this study, we introduced a lumped model of BOP fluid system and coupled it with NSSS without any ideal boundary conditions, i.e., in a closed loop. Various methods for coupling boundary conditions at MS/FW were tested to validate their combination in terms of minimizing numerical instability, which mainly arises from the coupled boundaries. The method exhibiting the best performance was selected and applied to a transient simulation of an integrated NSSS and BOP system of a SMART. For a transient event with core power change of 100-20-100%, the simulation exhibited numerical stability throughout the system without any significant perturbation of thermal-hydraulic parameters. Thus, the introduced boundary-condition coupling method and BOP fluid system model can expectedly be employed for the transient simulation and performance analysis of SMRs requiring daily load-following operations.

• Structural Engineering and Mechanics
• /
• 제89권6호
• /
• pp.557-570
• /
• 2024

Due to interfacial ageing, chemical action and interfacial damage, the interface debonding may appear in the interfaces of composite laminates. Particularly, the laminates display a side-dependent effect at small scale. In this work, a three-dimensional (3D) and anisotropic thick nanoplate model is proposed to investigate the effects of imperfect interface and nonlocal parameter on the bending deformation, vibrational response and buckling stability of one-dimensional (1D) hexagonal quasicrystal (QC) layered nanoplates. By combining the linear spring model with the transferring matrix method, exact solutions of phonon and phason displacements, phonon and phason stresses of bending deformation, the natural frequencies of vibration and the critical buckling loads of 1D hexagonal QC layered nanoplates are derived with imperfect interfaces and nonlocal effects. Numerical examples are illustrated to demonstrate the effects of the imperfect interface parameter, aspect ratio, thickness, nonlocal parameter, and stacking sequence on the bending deformation, the vibrational response and the critical buckling load of 1D hexagonal QC layered nanoplate. The results indicate that both the interface debonding and nonlocal effect can reduce the stiffness and stability of layered nanoplates. Increasing thickness of QC coatings can enhance the stability of sandwich nanoplates with the perfect interfaces, while it can reduce first and then enhance the stability of sandwich nanoplates with the imperfect interfaces. The biaxial compression easily results in an instability of the QC layered nanoplates compared to uniaxial compression. QC material is suitable for surface layers in layered structures. The mechanical behavior of QC layered nanoplates can be optimized by imposing imperfect interfaces and controlling the stacking sequence artificially. The present solutions are helpful for the various numerical methods, thin nanoplate theories and the optimal design of QC nano-composites in engineering practice with interfacial debonding.

• 대기
• /
• 제24권1호
• /
• pp.49-68
• /
• 2014

This study analyzed atmospheric conditions for the convergent cloud band (Cu-Cb line) in developing stage and its neighbouring convections formed over the East Sea on 1 February 2012, by using synoptic, satellites data, and WRF numerical simulation output of high resolution. In both satellite images and the WRF numerical simulation outputs, the Cu-Cb line that stretched out toward northwest-southeast was shown in the East Sea, and cloud lines of the L mode were aligned in accordance with the prevailing surface wind direction. However, those of the T mode were aligned in the direction of NE-SW, which was nearly perpendicular direction to the surface winds. The directions of the wind shear vectors connecting top winds and bottom winds of the moist layers of the L mode and the T mode were identical with those of the cloud lines of L mode and T mode, respectively. From the WRF simulation convection circulations with a convergence in the lower layer of atmosphere and a divergence above 1.5 km ASL (Above Sea Level) were identified in the Cu-Cb line. A series of small sized vortexes (maximum vortex: $320{\times}10^{-5}s^{-1}$) of meso-${\gamma}$-scale formed by convergences was found along the Cu-Cb lines, suggesting that Cu-Cb lines, consisting of numerous convective clouds, were closely associated with a series of the small vortexes. There was an absolute unstable layer (${\partial}{\theta}/{\partial}z$ < 0) between sfc and ~0.3 km ASL, and a stable layer (${\partial}{\theta}/{\partial}z$ > 0) above ~2 km ASL over the Cu-Cb line and cloud zones. Not only convectively unstable layers (${\partial}{\theta}_e/{\partial}z$ < 0) but also neutral layers (${\partial}{\theta}_e/{\partial}z{\approx}=0$) in the lower atmosphere (sfc~1.5 km ASL) were scattered around over the cloud zones. Particularly, for the Cu-Cb line there were convectively unstable layers in the surface layer, and neutral layers (${\partial}{\theta}_e/{\partial}z{\approx}=0$) between 0.2 and ~1.5 km ASL over near the center of the Cu-Cb line, and the neutralization of unstable layers came from the release of convective instability.