• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2004.03a
• /
• pp.186-191
• /
• 2004

The reactive flowfield of the transverse injecting combustor has been studied using Euler-Lagrange method in order to develop an efficient solution procedure for the understanding of liquid spray combustion in the transverse injecting combustor which has been widely used in ramjets and turbojet afterburners. The unsteady two-dimensional gas-phase equations have been represented in Eulerian coordinates and the liquid-phase equations have been formulated in Lagrangian coordinates. The gas-phase equations based on the conservation of mass, momentum, and energy have been supplemented by combustion. The vaporization model takes into account the transient effects associated with the droplet heating and the liquid-phase internal circulation. The droplet trajectories have been determined by the integration of the Lagrangian equation in the flow field obtained from the separate calculation without considering the iterative effect between liquid and gas phases. The reported droplet trajectories had been found to deviate from the initial conical path toward the flow direction in the very end of its lifetime when the droplet size had become small due to evaporation. The integration scheme has been based on the TEACH algorithm for gas-phase equation, the second order Runge-Kutta method for liquid-phase equations and the linear interpolation between the two coordinate systems. The calculation results has shown that the characteristics of the droplet penetration and recirculation have been strongly influenced by the interaction between gas and liquid phases in such a way that most of the vaporization process has been confined to the wake region of the injector, thereby improving the flame stabilization properties of the flowfield.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.04a
• /
• pp.165-169
• /
• 2009

The effect of the Dirichlet boundary condition for the redistance equation of level set method on the solutionof sloshing problem is investigated by adopting four Dirichlet boundary conditions. For the solution of the incompressible Navier-Stokes equations, P1P1 four-step fractional finite element method is employed and a least-square finite element method is used for the solutions of the two hyperbolic type equations of level set method; advection and redistance equation. ALE (Arbitrary Lagrangian Eulerian) method is used to deal with a moving computational domain. It has been shown that the free surface motion in a sloshing tank is strongly dependent on the type of the Dirichlet boundary condition and the results of broken dam and sloshing problems using various Dirichlet boundary conditions are discussed and compared with the existing experimental results.

• 한국전산유체공학회:학술대회논문집
• /
• 2007.10a
• /
• pp.141-148
• /
• 2007

We present an immersed boundary (IB) method for 3D simulation of flappingflags in a uniform flow. The proposed formulation is manipulated on the basis of an efficient Navier-Stokes solver adopting the fractional step method and a staggered Cartesian grid system. A direct numerical method is developed to calculate the flag motion, with the elastic force treated implicitly. The fluid motion defined on an Eulerian grid and the flag motion defined on a Lagrangian grid are independently solved and the mass of flag is handled in a natural way. An additional momentum forcing is formulated from the flag motion equation in a way similar with the direct-forcing IB formulation and acts as the interaction force between the flag and ambient fluid. A series of numerical tests are performed and the present results are compared qualitatively and quantitatively with previous studies. The instantaneous flag motion is analyzed under different conditions and surrounding vortical structures are identified. The effects of physical parameters on the flapping frequency are studied.

• Journal of computational fluids engineering
• /
• v.12 no.2
• /
• pp.14-20
• /
• 2007

We present an improved immersed boundary method for computing incompressible viscous flow around an arbitrarily moving body on a fixed computational grid. The main idea is to incorporate feedback forcing scheme of virtual boundary method with Peskin's regularized delta function approach in order to use large CFL number and transfer quantities between Eulerian and Lagrangian domain effectively. From the analysis of stability limits and effects of feedback forcing gains, optimum regions of the feedback forcing are suggested.

• 한국전산유체공학회:학술대회논문집
• /
• 2011.05a
• /
• pp.483-487
• /
• 2011

In the present study, a CFD program based on a finite volume method was developed by using an unstructured polyhedral grid system for the accurate simulation with the complex geometry of computational domain. To simulate the transient flow induced by the moving solid object, the program used a fractional step method and a ALE (Algebric Lagrangian-Eulerian) method. The grid deformation for the moving of solid object were performed with a spring analogy based on the center coordinate of each computational grid. To verify the present program with these methodologies, the numerical results of the flow around the fixed and oscillating circular cylinder were compared with the previous numerical results.

• 한국전산유체공학회:학술대회논문집
• /
• 2007.04a
• /
• pp.23-29
• /
• 2007

We present an improved immersed boundary method for computing incompressible viscous flow around an arbitrarily moving body on a fixed computational grid The main idea is to incorporate feedback forcing scheme of virtual boundary method with Peskin's regularized delta function approach in order to use large CFL number and transfer quantities between Eulerian and Lagrangian domain effectively. From the analysis of stability limits and effects of feedback forcing gains, optimum regions of the feedback forcing are suggested.

• International Journal of Ocean System Engineering
• /
• v.1 no.1
• /
• pp.37-45
• /
• 2011

The nonlinear free-surface motions interacting with a floating body were investigated using the Moving Particle Semi-implicit (MPS) method proposed by Koshizuka and Oka [6] for incompressible flow. In the numerical method, more realistic Lagrangian moving particles were used for solving the flow field instead of the Eulerian approach with a grid system. Therefore, the convection terms and time derivatives in the Navier-Stokes equation can be calculated more directly, without any numerical diffusion, instabilities, or topological failure. The MPS method was applied to a numerical simulation of predicting the efficiency of floating-type breakwater interacting with waves.

• Korean Journal of Hydrosciences
• /
• v.5
• /
• pp.85-97
• /
• 1994

A hybrid finite difference method for the longitudinal dispersion equation, which is based on combining the Holly-Preissmann scheme with fifth-degree Hermite interpolating polynomial and the generalized Crank-Nicholson scheme, is described and comparatively evaluated with other characteristics-based numerical methods. Longitudinal dispersion of an instantaneously-loaded pollutant source is simulated, and computational results are compared with the exact solution. The present method is free from wiggles regardless of the Courant number, and exactly reproduces the location of the peak concentration. Overall accuracy of the computation increases for smaller value of the weighting factor, $\theta$of the model. Larger values of $\theta$ overestimates the peak concentration. Smaller Courant number yields better accuracy, in general, but the sensitivity is very low, especially when the value of $\theta$ is small. From comparisons with the hybrid method using cubic interpolating polynomial and with splitoperator methods, the present method shows the best performance in reproducing the exact solution as the advection becomes more dominant.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.25 no.12
• /
• pp.866-873
• /
• 2015

The transverse vibration of the deploying beam from rigid hub was analyzed by simulation and experiment. The linear governing equation of the deploying beam was obtained using the Euler-Bernoulli beam theory. To discretize the governing equation, the Galerkin method was used. After transforming the governing equation into the weak form, the weak form was discretized. The discretized equation was expressed by the matrix-vector form, and then the Newmark method was applied to simulate. To consider the damping effect of the beam, we conducted the modal test with various beam length. The mass proportional damping was selected by the relation of the first and second damping ratio. The proportional damping coefficient was calculated using the acquired natural frequency and damping ratio through the modal test. The experiment was set up to measure the transverse vibration of the deploying beam. The fixed beam at the carriage of the linear actuator was moved by moving the carriage. The transverse vibration of the deploying beam was observed by the Eulerian description near the hub. The deploying or retraction motion of the beam had the constant velocity and the velocity profile with acceleration and deceleration. We compared the transverse vibration results by the simulation and experiment. The observed response by the Eulerian description were analyzed.

• Proceedings of the KSME Conference
• /
• 2008.11a
• /
• pp.1487-1491
• /
• 2008

Detailed knowledge on the motion of blood cells flowing in micro-channels under simple shear flow and the influence of blood flow is essential to provide a better understanding on the blood rheological properties and blood cell aggregation. The microscopic behavior of red blood cell (RBCs) is numerically investigated using a fluid-structure interaction (FSI) method based on the Arbitrary-Lagrangian-Eulerian (ALE) approach and the dynamic mesh method (smoothing and remeshing) in FLUENT (ANSYS Inc., USA). The employed FSI method could be applied to the motions and deformations of a single blood cell and multiple blood cells, and the primary thrombogenesis caused by platelet aggregation. It is expected that, combined with a sophisticated large-scale computational technique, the simulation method will be useful for understanding the overall properties of blood flow from blood cellular level (microscopic) to the resulting rheological properties of blood as a mass (macroscopic).