• Geomechanics and Engineering
• /
• v.28 no.5
• /
• pp.437-454
• /
• 2022

The hydro-mechanical (HM) two-way sequential coupling in the TOUGH2-FLAC3D linking algorithm is validated completely and successfully in both M to H and H to M directions, which are initiated by mechanical surface loading for geomechanical validation and hydrological groundwater pumping for hydrogeological validation, respectively. For such complete and successful validation, a TOUGH2-FLAC3D linked numerical model is developed first by adopting the TOUGH2-FLAC3D linking algorithm, which uses the two-way (fixed-stress split) sequential coupling scheme and the implicit backward time stepping method. Two geomechanical and two hydrogeological validation problems are then simulated using the linked numerical model together with basic validation strategies and prerequisites. The second geomechanical and second hydrogeological validation problems are also associated with the Mandel effect and the Noordbergum and Rhade effects, respectively, which are three phenomenally well-known but numerically challenging HM effects. Finally, sequentially coupled numerical solutions are compared with either analytical solutions (verification) or fully coupled numerical solutions (benchmarking). In all the four validation problems, they show almost perfect to extremely or very good agreement. In addition, the second geomechanical validation problem clearly displays the Mandel effect and suggests a proper or minimum geometrical ratio of the height to the width for the rectangular domain to maximize agreement between the numerical and analytical solutions. In the meantime, the second hydrogeological validation problem clearly displays the Noordbergum and Rhade effects and implies that the HM two-way sequential coupling scheme used in the linked numerical model is as rigorous as the HM two-way full coupling scheme used in a fully coupled numerical model.

• Coupled systems mechanics
• /
• v.1 no.4
• /
• pp.361-379
• /
• 2012

This paper presents a fully coupled three-dimensional solver for the analysis of interaction between pulsatile flow and large deformation structure. A partitioned time marching algorithm is employed for the solution of the time dependent coupled discretised problem, enabling the use of highly developed, robust and well-tested solvers for each field. Conservative transfer of information at the fluid-structure interface is combined with an effective multi-predict-correct iterative scheme to enable implicit coupling of the interacting fields at each time increment. The three-dimensional unsteady incompressible fluid is solved using a powerful implicit time stepping technique and an ALE formulation for moving boundaries with second-order time accurate is used. A full spectrum of total variational diminishing (TVD) schemes in unstructured grids is allowed implementation for the advection terms and finite element shape functions are used to evaluate the solution and its variation within mesh elements. A finite element dynamic analysis of the highly deformable structure is carried out with a numerical strategy combining the implicit Newmark time integration algorithm with a Newton-Raphson second-order optimisation method. The proposed model is used to predict the wave flow fields of a particular flow-induced vibrational phenomenon, and comparison of the numerical results with available experimental data validates the methodology and assesses its accuracy. Another test case about three-dimensional biomedical model with pulsatile inflow is presented to benchmark the algorithm and to demonstrate the potential applications of this method.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.8
• /
• pp.1128-1138
• /
• 2000

Pumping action in ejector systems is generally achieved through the mixing of a high-velocity and high-energy stream with a lower-velocity and lower-energy stream within a duct. The design and performance evaluation of the ejector systems has developed as a combination of scale-model experiments, empiricism and theoretical analyses applicable only to very simplified configurations, because of the generic complexity of the flow phenomena. In order to predict the detailed performance characteristics of such systems, the flow phenomena throughout the operating regimes of the ejector system should be fully understood. This paper presents the computational results for the two-dimensional supersonic ejector system with a second throat. The numerical simulations are based on a fully implicit finite volume scheme of the compressible Reynolds-averaged Navier-Stokes equation in a domain that extends from the stagnation chamber to the diffuser exit. For a wide range of the operating pressure ratio the flow field inside the ejector system is investigated in detail. The results show that the supersonic ejector systems have an optimal throat area for the operating pressure ratio to be minimized.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.4 no.4
• /
• pp.70-77
• /
• 2000

Supersonic Dual Bell Nozzle (SDBN) is an altitude-adaptive propulsion nozzle achieved only by a nozzle wall inflection. In order to investigate the altitude adaptive capability and the effectiveness of this nozzle concept, the present study addresses a computational work of the flow through SDBN. Several types of the SDBNs are tested for a wide range of the pressure ratio which covers from an over-expended flow to a fully under-expended flow at the exit of the SDBN. Axisymmetric, compressible, Wavier-Stokes equations are numerically solved using a fully implicit finite volume differencing scheme. The present computational results reveal that the base nozzle length affects the shock wave system occurring inside SDBN. For a quit wide range of the pressure ratio the flow separation occurs at the nozzle inflection point. It is found that the maximum thrust coefficient is obtainable for the correct expansion state at the exit of SDBN.

• Proceedings of the KSME Conference
• /
• 2007.05b
• /
• pp.3191-3196
• /
• 2007

The present study addresses a method to operate a fuel-cell system effectively using a recirculation ejector which recycles wasted hydrogen gas. Configuration of a recirculation ejector is changed to investigate the flow behavior through it under varying operating conditions, and how such conditions affect the fuel-cell hydrogen cycle. The numerical simulations are based on a fully implicit finite volume scheme of the axisymmetric, compressible, Reynolds-Averaged, Navier-Stokes equations for hydrogen gas, and are compared with available experimental data for validation. The results show that a hydrogen recirculation ratio is effectively controlled by a configurational alteration within the operational region in which the recirculation passage doesn't plugged by a sonic line.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.10 no.5
• /
• pp.160-167
• /
• 2002

The transient thermoelastic analysis of automotive disk brakes with frictional contact is performed by using the finite element method. To analyze the thermoelastic behaviors occurring in disk brakes, the coupled heat conduction and elastic equations are solved. The fully implicit transient scheme is used to improve the computation accuracy at every time step. The numerical results of the thermoelastic behaviors are obtained during the repeated braking condition. The computational results show that the thermoelastic instability(TEI) phenomenon(the growth of non-uniformities in contact pressure) occurs in disk brakes. Also, the effect of material properties on the thermoelastic behaviors is investigated to facilitate the conceptual design of the brake system.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.13 no.3
• /
• pp.167-173
• /
• 2001

This paper represents the numerical analysis for microenvironments various temperature and humidity in bed mattress. He purpose of this study is for healthful bed mattress by controling a bacteria with a prediction of the vapor and temperature distributions in the bed mattress. The numerical model is one dimensional unsteady state and the governing equations were discretized by fully implicit scheme. The numerical results were compared with experimental data, and showed a good agreement with them. Specially, the excess-relative humidity shows a lower distribution near the surface of mattress, meaning that the optimum living condition for bacteria will be caused.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2003.10a
• /
• pp.95-98
• /
• 2003

The physics of the plume-induced shock and separation particulary at a high plume to exit pressure ratio and supersonic speeds up to Mach 3.0 with aid without a passive control method, porous extension, were studied using computational techniques. Mass-averaged Navier-Stokes equations with the RNG k-$\varepsilon$ turbulence model were solved using a fully implicit finite volume scheme and a 4-stage Runge-Kutta method. The courol methodology for plume-afterbody interactions is to use a perforated wall attached at either the nozzle exit or the edge of the missile base. The Effect of porous wall length on plume interference is also investigated. The computational results show the main effect of the porous extension on plume-afterbody interactions is to in the plume from strongly underexpanding during a change in flight conditions. With control, a change in porous extension length has no significant effect on plume interference.

• 유체기계공업학회:학술대회논문집
• /
• 2006.08a
• /
• pp.223-226
• /
• 2006

Computational Study of a sphere subjected to free stream flow and simultaneously subjected to spinning motion is carried out. Three dimensional compressible Navier-Stokes equations are solved using fully implicit finite volume scheme. SST(Shear Stress Transport) $k-{\omega}$ turbulence model is used. Aerodynamic characteristics being affected are studied. Validation of the numerical process is done for the no spin condition. Variation of drag coefficient and shock wave strength with increase in spinning rate is reported. Changes in the wake region of the sphere with respect to spinning speed are also observed.

• Proceedings of the KSME Conference
• /
• 2000.11b
• /
• pp.485-490
• /
• 2000

This paper dipicts the computational results for the axisymmetric subsonic/sonic ejector systems with a second throat. The numerical simulations are based on a fully implicit finite volume scheme of the compressible Reynolds-Averaged Navier-Stokes equation in a domain that extends form the stagnation chamber to the ejector diffuser exit. In order to obtain practical design factors for subsonic/sonic ejector systems, the ejector throat area, the mixing section configuration, and the ejector throat length were changed in computations. For the subsonic/sonic ejector systems operating in the range of low operation pressure ratio, the effects of the design factors on the flow are discussed.