• Ocean Systems Engineering
• /
• v.12 no.1
• /
• pp.63-86
• /
• 2022

A refined projection-based purely Lagrangian meshfree method is presented towards reliable numerical analysis of fluid flow interactions with saturated/unsaturated porous media of uniform/spatially-varying porosities. The governing equations are reformulated on the basis of two-phase mixture theory with incorporation of volume fraction. These principal equations of mixture are discretized in the context of Incompressible SPH (Smoothed Particle Hydrodynamics) method. Associated with the consideration of governing equations of mixture, a new term arises in the source term of PPE (Poisson Pressure Equation), resulting in modified source term. The linear and nonlinear force terms are included in momentum equation to represent the resistance from porous media. Volume increase of fluid particles are taken into consideration on account of the presence of porous media, and hence multi-resolution ISPH framework is also incorporated. The stability and accuracy of the proposed method are thoroughly examined by reproducing several numerical examples including the interactions between fluid flow and saturated/unsaturated porous media of uniform/spatially-varying porosities. The method shows continuous pressure field, smooth variations of particle volumes and regular distributions of particles at the interface between fluid and porous media.

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

A similarity solution of the Navier-Stokes equation for the axisymmetric stagnation flow near a plane wall coated with a magnetic fluid of uniform thickness is constructed. The shape functions representing the flow in two (magnetic and normal) fluid layer are determined from a third order boundary value problem, which is solved by the Runge-Kutta method with two shooting parameters. Features of the flow including streamline pattern and interface velocity are investigated for the varying values of density ratio, viscosity ratio, and Reynolds number. The results for the interface and wall shear stress, boundary layer and displacement thickness are also presented.

• Journal of computational fluids engineering
• /
• v.14 no.2
• /
• pp.46-51
• /
• 2009

The slip effect from the molecular interaction between fluid particles and solid surface atoms plays a key role in microscale fluid transport and heat transfer since the relative importance of surface forces increases as the size of the system decreases to the microscale. There exist two models to describe the slip effect: the Maxwell slip model in which the slip correction is made on the basis of the degree of shear stress near the wall surface and the Langmuir slip model based on a theory of adsorption of gases on solids. In this study, as the first step towards developing a general purpose numerical code of the compressible Navier-Stokes equations for computational simulations of microscale fluid flow and heat transfer, two slip models are implemented into a finite element numerical code of a simplified equation. In addition, a pressure-driven gas flow in a microchannel is investigated by the numerical code in order to validate numerical results.

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

In this presentation, I talk about various fluid simulation methods that have been developed for computer graphics special effects since 1996. They are all based on CFD but sacrifice physical reality for visual plausability and time. But as the speed of computer increases rapidly and the capability of GPU (graphics processing unit) improves, methods for more physical realism have been tried. In this talk, I will focus on four aspects of fluid simulation methods for computer graphics: (1) particle level-set methods, (2) particle-based simulation, (3) methods for exact satisfaction of incompressibility constraint, and (4) GPU-based simulation. (1) Particle level-set methods evolve the surface of fluid by means of the zero-level set and a band of massless marker particles on both sides of it. The evolution of the zero-level set captures the surface in an approximate manner and the evolution of marker particles captures the fine details of the surface, and the zero-level set is modified based on the particle positions in each step of evolution. (2) Recently the particle-based Lagrangian approach to fluid simulation gains some popularity, because it automatically respects mass conservation and the difficulty of tracking the surface geometry has been somewhat addressed. (3) Until recently fluid simulation algorithm was dominated by approximate fractional step methods. They split the Navier-Stoke equation into two, so that the first one solves the equation without considering the incompressibility constraint and the second finds the pressure which satisfies the constraint. In this approach, the first step introduces error inevitably, producing numerical diffusion in solution. But recently exact fractional step methods without error have been developed by fluid mechanics scholars), and another method was introduced which satisfies the incompressibility constraint by formulating fluid in terms of vorticity field rather than velocity field (by computer graphics scholars). (4) Finally, I want to mention GPU implementation of fluid simulation, which takes advantage of the fact that discrete fluid equations can be solved in parallel.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2006.11a
• /
• pp.848-854
• /
• 2006

In this study, a computational analysis system has been developed in order to investigate flow-induced vibration(FIV) phenomenon for general stator-rotor cascade configurations. Relative movement of the rotor with respect to stator is reflected by modeling independent two computational domains. Fluid domains are modeled using the unstructured grid system with dynamic moving and local deforming methods. Unsteady, Reynolds-averaged Navier-Stokes equations with one equation Spalart-Allmaras and two-equation SST $k-\omega$ turbulence models are solved for unsteady flow problems. A fully implicit time marching scheme based on the Newmark direct integration method is used flow computing the coupled governing equations of the fluid-structure interaction problem. Detailed FIV responses for different flow conditions are presented with respect to time and vibration characteristics are also physically investigated in the time domain.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.16 no.10 s.115
• /
• pp.1082-1089
• /
• 2006

In this study, advanced computational analysis system has been developed in order to investigate flow-induced vibration(FIV) phenomenon for general stator-rotor cascade configurations. Relative movement of the rotor with respect to stator is reflected by modeling Independent two computational domains. Fluid domains are modeled using the unstructured grid system with dynamic moving and local deforming methods. Unsteady, Reynolds-averaged Wavier-stokes equations with one equation Spalart-Allmaras and two-equation SST ${\kappa}-{\varepsilon}$ turbulence models are solved for unsteady flow problems and also relative moving and vibration effects of the rotor cascade are fully considered. A coupled implicit time marching scheme based on the Newmark integration method is used for computing the governing equations of fluid-structure interaction problems. Detailed vibration responses for different flow conditions are presented and then vibration characteristics are physically investigated in the time domain as computational virtual tests.

• Bulletin of the Korean Chemical Society
• /
• v.18 no.8
• /
• pp.841-850
• /
• 1997

A simple molecular theory of mixtures is formulated based on the nonrandom two-fluid lattice-hole theory of fluids. The model is applicable to mixtures over a density range from zero to liquid density. Pure fluids can be completely characterized with only two molecular parameters and an additional binary interaction energy is required for a binary mixture. The thermodynamic properties of ternary and higher order mixtures are completely defined in terms of the pure fluid parameters and the binary interaction energies. The Quantitative prediction of vapor-liquid, and solid-vapor equilibria of various mixtures are demonstrated. The model is useful, in particular, for mixtures whose molecules differ greatly in size. For real mixtures, satisfactory agreements are resulted from experiment. Also, the equation of state (EOS) is characterized well, even the liquid-liquid equilibria behaviors of organic mixtures and polymer solutions with a temperature-dependent binary interaction energy parameter.

• Bulletin of the Korean Chemical Society
• /
• v.18 no.9
• /
• pp.965-972
• /
• 1997

An approximate molecular theory of classical fluids based on the nonrandom lattice statistical-mechanical theory is presented. To obtain configurational Helmholtz free energy and equation of state (EOS), the lattice-hole theory of the Guggenheim combinatorics is approximated by introducing the nonrandom two-fluid theory. The approximate nature in the derivation makes the model possible to unify the classical lattice-hole theory and to describe correctly the configurational properties of real fluids including macromolecules. The theory requires only two molecular parameters for a pure fluid. Results obtained to date have demonstrated that the model correlates quantitatively the first- and second-order thermodynamic properties of real fluids. The basic simplicity of the model can readily be generalized to multicomponent systems. The model is especially relevant to (multi) phase equilibria of systems containing molecularly complex species.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.3
• /
• pp.422-432
• /
• 2001

An analytical characterization is presented on the transient heat transfer by an oscillating flow through a porous slab. Based on a two-equation model, analytic solutions are obtained for both the fluid and solid temperature variations. Two parameters are identified as the Stanton number for the internal heat exchange and the ratio of the thermal capacities between the solid and fluid phases. The heat transfer characteristics are shown to be classified into four regimes according to the two parameters and physical interpretation is presented on the particular heat transfer processes within each regime. In addition, the condition for the local thermal equilibrium between the phases is examined and the relevant criterion is suggested.

• Journal of the Korean Society of Safety
• /
• v.18 no.1
• /
• pp.101-107
• /
• 2003

Dynamic behavior of flexible rectangular liquid storage structures is analysed by the developed method. The rectangular liquid storage structures are assumed to be fixed to the ground and a moving coordinate system is used. The irrotational motion of invicid and incompressible ideal fluid is represented by two analytic solutions. One is the solution of the fluid motion in the rigid rectangular liquid storage structure due to ground motions and the other is the solution of the fluid motion by the motion of the wall in the flexible rectangular liquid storage structure. The motion of structure is modeled by finite elements. The fluid-structure interaction effect is reflected into the coupled equation of motion as added fluid mass matrix. The free surface sloshing motion and hydrodynamic pressure acting on the wall in the flexible rectangular liquid storage structure due to the horizontal ground motion are obtained by the developed method and verified.