• Coupled systems mechanics
• /
• v.7 no.4
• /
• pp.421-440
• /
• 2018

In this study, a computational method for wall shear stress combined with an implicit direct-forcing immersed boundary method is presented. Near the immersed boundaries, the sub-grid stress is determined by a wall model in which the wall shear stress is directly calculated from the Lagrangian force on the immersed boundary. A coupling mathematical model of the transition process for a model Francis turbine comprising turbulent flow and rotating rigid guide vanes is established. The spatiotemporal distributions of pressure, velocity, vorticity and turbulent quantity are gained with the transient process; the drag and lift coefficients as well as other forces (moments) are also obtained as functions of the attack angle. At the same time, analysis is conducted of the characteristics of pressure pulsation, velocity stripes and vortex structure at some key parts of flowing passage. The coupling relations among the turbulent flow, the dynamical force (moment) response of blade and the rotating of guide vane are also obtained.

• 한국전산유체공학회:학술대회논문집
• /
• 2006.05a
• /
• pp.41-42
• /
• 2006

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

IB (immersed boundary) method is one of the prominent tool in computational fluid dynamics for the analysis of flows over complex geometries. The IB technique simplyfies the solution procedure by eliminating the requirement of complex body fitted grids and it is also superior in terms of memory requirement. In this study we have developed numerical code (FOTRAN) for the analysis of 2D flow over a cylinder using IB technique. The code is validated by comparing the wake lengths and separation angles given by Guo et. al. We employed fractional-step procedure for solving the Navier-Stokes equations governing the flow and discrete forcing IB technique for imposing boundary conditions. Also we have developed a 3D code for the backward-facing-step flow and flow over a sphere. The reattachment length in backward-facing-step flow was compared with the one given by Nie and Armaly, which has proven the validity of our code.

• Journal of the Society of Naval Architects of Korea
• /
• v.58 no.5
• /
• pp.271-280
• /
• 2021

The code developed using a pressure-based method for unified conservation laws of incompressible/compressible fluids is expanded to handle moving or deforming body boundaries using the hybrid Cartesian/immersed boundary method. An instantaneous pressure field is calculated from a pressure Poisson equation for the whole fluid domain, including the compressible gas region. The polytropic gas is assumed for the compressible fluid so that the energy equation is decoupled. Immersed boundary nodes are identified based on edges crossing body boundaries. The velocity vector is reconstructed at the immersed boundary node using an interpolation along the assigned local normal line. The developed code is validated by comparing the time histories of pressure and wave elevation for sloshing in a rectangular and a membrane-type tank. The validated code is applied to simulate air cushion effects in a rectangular tank under sway motion. Time variations of pressure fields are analyzed in detail as the air pocket pulsates. It is shown that the contraction and expansion of the air pocket dominate the pressure loads on the wall of the tank. The present results are in good agreement with other experimental and computational results for the amplitude and the decay of the pressure oscillations measured at the pressure gauges.

• Geomechanics and Engineering
• /
• v.18 no.6
• /
• pp.651-659
• /
• 2019

The IB-SEM numerical method combines the spectral/hp element method and the rigid immersed boundary method. This method avoids the problems of low computational efficiency and errors that are caused by the re-division of the grid when the solids move. Based on the Fourier transformation and the 3D immersed boundary method, the 3D IB-SEM system was established. Then, using the open MPI and the Hamilton HPC service, the computational efficiency was increased substantially. The flows around a cylinder and a sphere were simulated by the system. The surface of the cylinder generates vortices with alternating shedding, and these vortices result in a periodic force acting on the surface of the cylinder. When the shedding vortices enter the flow field behind the cylinder, a recirculation zone is formed. Finally, the three-dimensional pore flow was successfully investigated.

• 한국전산유체공학회:학술대회논문집
• /
• 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.14 no.4
• /
• pp.41-47
• /
• 2009

An immersed boundary method(IBM, Kim et al.(2001)) for simulating flows over complex geometries is applied to computation of three-dimensional Floquet stability of a periodic wake. Floquet stability analysis is employed to extract different modes of three-dimensional instability. To verify the present method, a fully-resolved Floquet stability calculation for flow past a circular cylinder is considered. There are two different instability modes with long(mode A) and short (mode B) spanwise wavelengths for the periodic wake of a circular cylinder. The critical Reynolds number and the most unstable spanwise wavelengths of modes A and B are computed using the present method, and compared with other authors' results currently available.

• Proceedings of the KSME Conference
• /
• 2007.05b
• /
• pp.2731-2736
• /
• 2007

In this paper, incompressible flow over a cylinder near a plane wall using the Immersed Boundary. Finite Difference Lattice Boltzmann Method (IB-FDLBM) is implemented. In this present method, FDLBM is mixed with IBM by using the equilibrium velocity. We introduce IBM so that we can easy to simulate bluff-bodies. With this numerical procedure, the flow past a circular cylinder near a wall is simulated. We calculated the flow patterns about various Reynolds numbers and gap ratios between a circular cylinder and plane wall. So these are enabled to observe for vortex shedding. The numerical results are found to be in good agreement with those of previous studies.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.29 no.8 s.239
• /
• pp.940-947
• /
• 2005

An immersed boundary method for simulation of density-stratified flows has been developed and applied to computation of viscous flows past three different types of obstacle under table density stratification, namely laminar flows past a vertical barrier, a cosine hill, and a sphere, respectively. Density forcing is introduced on the body surface or inside the body. Significant changes in flow characteristics are observed depending on Fr. The numerical results are in good agreement with other authors' experimental and numerical results currently available, and shed light on computation of density-stratified flows in complex geometries.

• Proceedings of the KSME Conference
• /
• 2004.04a
• /
• pp.1909-1914
• /
• 2004

An immersed boundary method for simulation of density-stratified flows is developed and applied to computation of viscous flows over two-dimensional obstacles in a bounded domain under stable density stratification. Density sources/sinks are introduced on the body surface. Two obstacle shapes are used, a vertical barrier and a smooth cosine-shaped hill; weak stratification, defined by $K=ND/{\pi}U{\leq}1$, where U, N, and D are the upstream velocity, buoyancy frequency, and domain height, respectively, is considered. The results are consistent with other authors' calculations, and shed light on computation of density-stratified flows in complex geometries.