• Journal of Korea Water Resources Association
• /
• v.56 no.spc1
• /
• pp.1049-1058
• /
• 2023

The flow passing through river-crossing structures such as weirs and low-fall dams is dominated by rapidly varied flow including hydraulic jump. The intense unsteadiness of flow velocity and free surface profile affects the stability of such hydraulic structures. In particular, the steady hydraulic jump generated at high Froude number conditions includes remarkably air entrainment, making the flow characteristics more complicated. In this study, a large-eddy simulation was performed for turbulence effect and the hybrid VoF technique to simulate the steady hydraulic jump at the Froude number of 7.3 and the Reynolds number of 15,700. The results of the numerical simulation showed that the instantaneous maximum pressure and time-average pressure distribution calculated on the bottom surface downstream of the structure could be reasonably well reproduced being in good agreement with the experimental values. However, the instantaneous minimum pressure distribution in the direct downstream of the structure shows the opposite pattern to the target experimental measurement value. However, the numerical simulation performed in this study is considered to reasonably predict the minimum pressure distributions observed in various experiments conducted at similar conditions. The vertical distributions of flow velocity and air concentration computed in the center of the hydraulic jump were found to be in good agreement with the experimental results measured under similar conditions, showing self-similarity. These results show that the large eddy simulation and hybrid VoF techniques applied in this study can reproduce the hydraulic jump with strong air entrainment and the resulting intense free surface and pressure fluctuations at high Froude number conditions.

• Fibers and Polymers
• /
• v.7 no.4
• /
• pp.389-397
• /
• 2006

The beneficial effects of graduated compression stockings (GCS) in prophylaxis and treatment of venous disorders of human lower extremity have been recognized. However, their pressure functional performances are variable and unstable in practical applications, and the exact mechanisms of action remain controversial. Direct surface pressure measurements and indirect material properties testing are not enough for fully understanding the interaction between stocking and leg. A three dimensional (3D) biomechanical mathematical model for numerically simulating the interaction between leg and GCS in dynamic wear was developed based on the actual geometry of the female leg obtained from 3D reconstruction of MR images and the real size and mechanical properties of the compression stocking prototype. The biomechanical solid leg model consists of bones and soft tissues, and an orthotropic shell model is built for the stocking hose. The dynamic putting-on process is simulated by defining the contact of finite relative sliding between the two objects. The surface pressure magnitude and distribution along the different height levels of the leg and stress profiles of stockings were simulated. As well, their dynamic alterations with time processing were quantitatively analyzed. Through validation, the simulated results showed a reasonable agreement with the experimental measurements, and the simulated pressure gradient distribution from the ankle to the thigh (100:67:30) accorded with the advised criterion by the European committee for standardization. The developed model can be used to predict and visualize the dynamic pressure and stress performances exerted by compression stocking in wear, and to optimize the material mechanical properties in stocking design, thus, helping us understand mechanisms of compression action and improving medical functions of GCS.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.37 no.4
• /
• pp.383-391
• /
• 2009

The detonation wave characteristics of JP-7 and oxygen mixture is investigated by one-step induction parameter model (IPM) obtained from a detailed chemistry mechanism. A general procedure of obtaining reliable one-step kinetics IPM for hydrocarbon mixture from the fully detailed chemistry is described in this study. The IPM is obtained by the reconstruction of the induction time database obtained from a detailed kinetics library. The IPM was confirmed by the comparison of the induction time calculations with that from detailed kinetics. The IPM is later implemented to a fluid dynamics code and applied for the numerical simulation of detonation wave propagation. The numerical results show the detailed characteristics of the detonation wave propagation in JP-7 and oxygen mixture at affordable computing time, which is not be possible by the direct application of the detailed chemical kinetics mechanism of hydrocarbon fuel combustion.

• Advances in nano research
• /
• v.6 no.3
• /
• pp.243-255
• /
• 2018

Polystyrene granules were combined with nanosilver to form a nanocomposite film. One-side migration was conducted to test into three food simulants (3% acetic acid, 10% ethanol and 95% ethanol) at $40^{\circ}C$ temperature on different period of time (2, 4, 6, 8 and 10 days). It was found that, among the simulants, the highest migration amount was obtained with 3% acetic acid, while the 95% ethanol revealed the least migration level. Diffusion coefficients of nanosilver particles into simulants were estimated by inverse simulation using experimental data of concentration variation in the simulants. The finite element method used to solve the mass transfer equation and the numerical results indicates the sameresponse with the experimental data. The numerical results confirmed that the highest diffusion coefficient for acetic acid 3% (1.82E-10 to $1.76E-9m^2\;s^{-1}$) and the lowest diffusion coefficient for ethanol 95% from 2 to 10 days were obtained, respectively. Also, results of diffusion coefficient - concentration relation showed, the diffusion coefficient had in direct correlation with time and concentration. The results indicated that, in the 3% acetic acid, due to the increasing of diffusion coefficient of silver nanoparticles, they are released faster and distributed uniformly.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.9
• /
• pp.1262-1272
• /
• 2003

A two-dimensional direct numerical simulation is performed to investigate the dynamic behaviors of a single vortex in counter reacting and non-reacting flow field. A predictor-corrector-type numerical scheme with a low Mach number approximation is used in this simulation. A 16-step augmented reduced mechanism is adopted to treat the chemical reaction. The budget of the vorticity transport equation is examined to reveal a mechanism leading to the formation, destruction and transport of a single vortex according to the direction of vortex generation in reacting and non-reacting flows. The results show that air-side vortex has more larger strength than that of fuel-side vortex in both non-reacting and reacting flows. In reacting flow, the vortex is more dissipated than that in non-reacting flow as the vortex approach the flame. The total circulation in reacting flow, however, is larger than that in non-reacting flow because the convection transport of vorticity becomes much large by the increased velocity near the flame region. It is also found that the stretching and the convection terms mainly generate vorticity in non-reacting and reacting flows. The baroclinic torque term generates vorticity, while the viscous and the volumetric expansion terms attenuate vorticity in reacting flow. Furthermore, the contribution of volumetric expansion term on total circulation for air-side vortex is much larger than that of fuel-side vortex. It is also estimated that the difference of total circulation near stagnation plane according to the direction of vortex generation mainly attributes to the convection term.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.19 no.1
• /
• pp.29-36
• /
• 2015

Base isolation is considered as a seismic protective system in the design of next generation Nuclear Power Plants (NPPs). If seismic isolation devices are installed in nuclear power plants then the safety under a seismic load of the power plant may be improved. However, with respect to some equipment, seismic risk may increase because displacement may become greater than before the installation of a seismic isolation device. Therefore, it is estimated to be necessary to select equipment in which the seismic risk increases due to an increase in the displacement by the installation of a seismic isolation device, and to perform research on the seismic performance of each piece of equipment. In this study, modified NRC-BNL benchmark models were used for seismic analysis. The numerical models include representations of isolation devices. In order to validate the numerical piping system model and to define the failure mode, a quasi-static loading test was conducted on the piping components before the analysis procedures. The fragility analysis was performed by using the results of the inelastic seismic response analysis. Inelastic seismic response analysis was carried out by using the shell finite element model of a piping system considering internal pressure. The implicit method was used for the direct integration time history analysis. In addition, the collapse load point was used for the failure mode for the fragility analysis.

• Computers and Concrete
• /
• v.21 no.6
• /
• pp.623-635
• /
• 2018

In order to research the influence of different recycled aggregate contents on the mechanical properties of pervious concrete, the experimental study and numerical simulation analysis of the mechanical properties of pervious concrete with five kinds of recycled aggregates contents (0%, 25%, 50%, 75% and 100%) are carried out in this paper. The experimental test were first performed on concrete specimens of different sizes in order to determine the influence of recycled aggregate on the compressive strength and splitting tensile strength, direct tension strength and bending strength. Then, the development of the internal cracks of pervious concrete under different working conditions is studied more intuitively by $PFC^{3D}$. The experimental results show that the concrete compressive strength, tensile strength and bending strength decrease with the increase of the recycled aggregate contents. This trend of reduction is not only related to the brittleness of recycled aggregate concrete, but also to the weak viscosity of recycled aggregate and cement paste. It is found that the fracture surface of pervious concrete with recycled aggregate is smoother than that of natural aggregate pervious concrete by $PFC^{3D}$, which means that the bridging effect is weakened in the stress transfer between the left and right sides of the crack. Through the analysis of the development of the internal cracks, the recycled aggregate concrete generated more cracks than the natural aggregate concrete, which means that the recycled aggregate concrete is easier to form a coalescence fracture surface and eventually break.

• Nuclear Engineering and Technology
• /
• v.56 no.4
• /
• pp.1310-1319
• /
• 2024

The reactor cavity cooling system (RCCS) is a passive reactor safety system commonly present in the designs of High-Temperature Gas-cooled Reactors (HTGR) that removes heat from the reactor pressure vessel by means of natural convection and radiation. It is one of the factors responsible for ensuring that the reactor does not melt down under any plausible accident scenario. For the simulation of accident scenarios, which are transient phenomena unfolding over a span of up to several days, intermediate fidelity methods and system codes must be employed to limit the models' execution time. These models can quantify radiation heat transfer well, but heat transfer caused by natural convection must be quantified with the use of correlations for the heat transfer coefficient. It is difficult to obtain reliable correlations for HTGR RCCS heat transfer coefficients experimentally due to such a system's size. They could, however, be obtained from high-fidelity steady-state simulations of RCCSs. The Rayleigh number in RCCSs is too high for using a Direct Numerical Simulation (DNS) technique; thus, a Reynolds-Averaged Navier-Stokes (RANS) approach must be employed. There are many RANS models, each performing best under different geometry and fluid flow conditions. To find the most suitable one for simulating an RCCS, the RANS models need to be validated. This work benchmarks various RANS models against three experiments performed on the HTTR RCCS Mockup by the Japanese Atomic Energy Agency (JAEA) in 1993. This facility is a 1/6 scale model of a vessel cooling system (VCS) for the High Temperature Engineering Test Reactor (HTTR), which is operated by JAEA. Multiple RANS models were evaluated on a simplified 2d-axisymmetric geometry. They were found to reproduce the experimental temperature profiles with errors of up to 22% for the lowest temperature benchmark and 15% for the higher temperature benchmarks. The results highlight that the pragmatic turbulence models need to be validated for high Rayleigh natural convection-driven flows and improved accordingly, more publicly available experimental data of RCCS resembling experiments is needed and indicate that a 2d-axisymmetric geometry approximation is likely insufficient to capture all the relevant phenomena in RCCS simulations.

• Journal of Advanced Marine Engineering and Technology
• /
• v.33 no.8
• /
• pp.1261-1267
• /
• 2009

The Direct Torque Control[DTC] controls torque and flux by restricting the flux and torque errors within respective hysteresis bands, and motor torque and flux are controlled by the stator voltage space vector using optimum inverter switching table. And the Current Error Compensation method is on the basis of compensating current difference between the induction motor and its numerical model, in which the identical stator voltage is supplied for both the actual motor and the model so that the gaps between stator currents of the two can be forced to decay to zero as time proceeds. Consequently, the rotor speed approaches to the model speed, namely, setting value and the system can control motor speed precisely. This paper proposes a new sensorless speed control of induction motor using DTC and Current Error Compensation, which requires neither shaft encoder, speed estimator nor PI controllers. And through computer simulation, confirm effectiveness of proposed method.

• Journal of the Korean Society for Precision Engineering
• /
• v.29 no.5
• /
• pp.563-571
• /
• 2012

In 3D integration package using TSV technology, bonding is the core technology for stacking and interconnecting the chips or wafers. During bonding process, however, warpage and high stress are introduced, and will lead to the misalignment problem between two chips being bonded and failure of the chips. In this paper, a finite element approach is used to predict the warpages and stresses during the bonding process. In particular, in-plane deformation which directly affects the bonding misalignment is closely analyzed. Three types of bonding technology, which are Sn-Ag solder bonding, Cu-Cu direct bonding and SiO2 direct bonding, are compared. Numerical analysis indicates that warpage and stress are accumulated and become larger for each bonding step. In-plane deformation is much larger than out-of-plane deformation during bonding process. Cu-Cu bonding shows the largest warpage, while SiO2 direct bonding shows the smallest warpage. For stress, Sn-Ag solder bonding shows the largest stress, while Cu-Cu bonding shows the smallest. The stress is mainly concentrated at the interface between the via hole and silicon chip or via hole and bonding area. Misalignment induced during Cu-Cu and Sn-Ag solder bonding is equal to or larger than the size of via diameter, therefore should be reduced by lowering bonding temperature and proper selection of package materials.