• 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.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.19 no.2
• /
• pp.197-215
• /
• 2015

Modeling water flow in variably saturated, porous media is important in many branches of science and engineering. Highly nonlinear relationships between water content and hydraulic conductivity and soil-water pressure result in very steep wetting fronts causing numerical problems. These include poor efficiency when modeling water infiltration into very dry porous media, and numerical oscillation near a steep wetting front. A one-dimensional finite element formulation is developed for the numerical simulation of variably saturated flow systems. First order backward Euler implicit and second order Crank-Nicolson time discretization schemes are adopted as a solution strategy in this formulation based on Picard and Newton iterative techniques. Five examples are used to investigate the numerical performance of two approaches and the different factors are highlighted that can affect their convergence and efficiency. The first test case deals with sharp moisture front that infiltrates into the soil column. It shows the capability of providing a mass-conservative behavior. Saturated conditions are not developed in the second test case. Involving of dry initial condition and steep wetting front are the main numerical complexity of the third test example. Fourth test case is a rapid infiltration of water from the surface, followed by a period of redistribution of the water due to the dynamic boundary condition. The last one-dimensional test case involves flow into a layered soil with variable initial conditions. The numerical results indicate that the Crank-Nicolson scheme is inefficient compared to fully implicit backward Euler scheme for the layered soil problem but offers same accuracy for the other homogeneous soil cases.

• 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.

• Journal of Korean Port Research
• /
• v.13 no.1
• /
• pp.167-174
• /
• 1999

Numerical analysis of two-fluid flows for both water and air is carried out. Free-Surface flows with an arbitrary deformation have been simulated around two dimensional submerged hydrofoil. The computation is performed using a finite volume method with unstructured meshes and an interface capturing scheme to determine the shape of the free surface. The method uses control volumes with an arbitrary number of faces and allows cell-wise local mesh refinement. the integration in space is of second order based on midpoint rule integration and linear interpolation. The method is fully implicit and uses quadratic interpolation in time through three time levels The linear equation systems are solved by conjugate gradient type solvers and the non-linearity of equations is accounted for through picard iterations. The solution method is of pressure-correction type and solves sequentially the linearized momentum equations the continuity equation the conservation equation of one species and the equations or two turbulence quantities.

• Proceedings of the KSME Conference
• /
• 2001.06e
• /
• pp.305-310
• /
• 2001

This paper presents the computational method for analyzing the compressible flow fields in a high voltage gas circuit breaker. There are many difficult problems in analyzing the gas flow in GCB due to complex geometry, moving boundary, shock wave and so on. In particular, the distortion problem of the grid due to the movement of moving parts can be worked out by the fixed grid technique. Numerical simulations are based on a fully implicit finite volume method of the compressible Reynolds-averaged Navier-Stokes equations to obtain the pressure, density, and velocity through the entire interruption process. The presented method is applied to the real circuit breaker model and the pressure in front of the piston is good agreement with the experimental one.

• Journal of Advanced Marine Engineering and Technology
• /
• v.32 no.8
• /
• pp.1215-1220
• /
• 2008

To obtain the shape of the free surface more accurately, computations are carried out by a finite volume method using unstructured meshes and an interface capturing method. Free-surface flow, which is very important in the fields of ship and marine engineering, is numerically simulated for flows of both water and air. Control volumes are used with an arbitrary number of faces and allows a local mesh refinement. The integration is of second order, with a midpoint rule integration and linear interpolation. The method is fully implicit and uses quadratic interpolation. The solution method of pressure-correction type solves sequentially equations of momentum, continuity, conservation, and two-equations turbulence model. Comparison are quantitatively made between the computation and experiment in order to confirm the solution method.

• 한국전산유체공학회:학술대회논문집
• /
• 2004.03a
• /
• pp.166-171
• /
• 2004

Direct numerical simulations are peformed to investigate the physics of a spatially developing turbulent boundary layer flow suddenly subjected to spanwise oscillating electro-magnetic forces in the near-wall region. The Reynolds number based on the inlet momentum thickness and free-stream velocity is $Re_\theta=300$. A fully-implicit fractional step method is employed to simulate the flow. The mean flow properties and the Reynolds stresses are obtained to analyze the near-wall turbulent structure. It is found that skin-friction and turbulent kinetic energy can be reduced by the electro-magnetic forces. Instantaneous flow visualization techniques are used to observe the response of streamwise vortices to spanwise oscillating forces. The near-wall vortical structures are clearly affected by spanwise oscillating electro-magnetic forces.

• Journal of computational fluids engineering
• /
• v.10 no.1
• /
• pp.16-23
• /
• 2005

The objective of this research is to estimate two dimensional ablating and charring of heat shield materials in severe aero-thermal heat transfer. This estimation requires an accurate and rapid technique for its serious heat transfer with a moving boundary. Aerodynamic heating is obtained by an explicit relation which is a function of Mach number and air condition, while a fully implicit method is used for heat transfer calculations. Moving boundary is captured by FLIAR method which is a subgroup of VOF. Thickness of ablating and charring of heat shield, temperature of the moving surface and rate of radiation heat are calculated and compared with references. The results are in good agreement with other calculations.