• Journal of Mechanical Science and Technology
• /
• v.15 no.5
• /
• pp.665-670
• /
• 2001

The interaction of a travelling shock with the shear layer of a flat plate is studied computationally. The Euler and Navier-Stokes equations are solved numerically on quadrilateral unstructured adaptive grids. The flat plate is installed horizontally on the central axis of a shock tube. The shear layer is first created by two shock waves at different speeds splitted by a flat plate. A series of small vortices is developed as a consequence in the shear layer. The shock wave reflected at the end wall impinges the shear layer. The complicated shock dynamics in the evolution to the pseudo-steady state is represented with the morphological transformation of a planar shock into an oblique shock.

• Steel and Composite Structures
• /
• v.36 no.2
• /
• pp.179-186
• /
• 2020

In this paper, the forced resonance vibration of porous functionally graded (FG) curved nanobeam is examined. In order to capture the hardening and softening mechanisms of nanostructure, the nonlocal strain gradient theory is employed to build the size-dependent model. Using the Timoshenko beam theory together with the Hamilton principle, the equations of motion for the curved nanobeam are derived. Then, Navier series are used in order to obtain the dynamical deflections of the porous FG curved nanobeam with simply-supported ends. It is found that the resonance position of the nanobeam is very sensitive to the nonlocal and strain gradient parameters, material variation, porosity coefficient, as well as geometrical conditions. The results indicate that the resonance position is postponed by increasing the strain gradient parameter, while the nonlocal parameter has the opposite effect on the results. Furthermore, increasing the opening angle or length-to-thickness ratio will result in resonance position moves to lower-load frequency.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.11a
• /
• pp.121-126
• /
• 1998

A numerical optimization procedure is developed to find the inlet velocity profile that yields the most uniform epitaxial layer in a vertical MOCVD reactor. It involves the solution of fully elliptic equations of motion, temperature, and concentration; the finite volume method based on SIMPLE algorithm has been adopted to solve the Navier-Stokes equations. The overall optimization process is highly nonlinear and has been efficiently treated by the sequential linear programming technique that breaks the non-linear problem into a series of linear ones. The optimal profile approximated by a 6th-degree Chebyshev polynomial is very successful in reducing the spatial non-uniformity of the growth rate. The optimization is particularly effective to the high Reynolds number flow. It is also found that a properly constructed inlet velocity profile can suppress the buoyancy driven secondary flow and improve the growth-rate uniformity.

• 한국전산유체공학회:학술대회논문집
• /
• 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.11 no.4 s.35
• /
• pp.39-47
• /
• 2006

In this paper, the authors developed a multidimensional interpolation method inside a finite volume cell in the computation of high-order accurate numerical flux such as the fifth order WEND (weighted essentially non-oscillatory) scheme. This numerical method starts from a simple Taylor series expansion in a proper spatial order of accuracy, and the WEND filter is used for the reconstruction of sharp nonlinear waves like shocks in the compressible flow. Two kinds of interpolations are developed: one is for the cell-averaged values of conservative variables divided in one mother cell (Type 1), and the other is for the vertex values in the individual cells (Type 2). The result of the present study can be directly used to the cell refinement as well as the convective flux between finer and coarser cells in the Cartesian adaptive grid system (Type 1) and to the post-processing as well as the viscous flux in the Navier-Stokes equations on any types of structured and unstructured grids (Type 2).

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

In the present study, the numerical prediction of the oil amount leaked from the hole of a damaged tank is investigated using the improved MPS (Moving Particle Semi-implicit) method, which was originally proposed by Koshizuka and Oka (1996) for incompressible flow. The governing equations, which consist of the continuity and Navier-Stokes equations, are solved by Lagrangian moving particles, and all terms expressed by differential operators should be replaced by the particle interaction models based on a Kernel function. The simulation results are validated though the comparison with the analytic solution based on Torricelli's equilibrium relation. Furthermore, a series of numerical simulations under the various conditions are performed in order to estimate more accurately the initial amount of leaked oil.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.11
• /
• pp.1467-1475
• /
• 2001

In the present paper, approximate coordinate transformations for simulation of turbulent flows with active wall motions, leading to a significant reduction in the computational cost while maintaining the numerical accuracy, are presented: the Navier-Stokes equations are coordinate-transformed with an approximation of Taylor-series truncation and neglect of some less-significant terms. The performance of the proposed transformations is evaluated in simulation of the channel flow at Re$\sub$$\tau$/=140 with wall deformations of │η$\sub$m/$\^$+/ 5. The approximate transformations provide flow structures as wall as turbulence statistics in good agreement with those from a complete coordinate transformation [Phys. Fluids 12, 3301 (2000)] and allow 25-30% savings in the CPU time as compared to the complete one.

• 한국전산유체공학회:학술대회논문집
• /
• 1999.05a
• /
• pp.66-72
• /
• 1999

A numerical simulation method has been under development for solving turbulent flows around a ship model in maneuvering motion using the Reynolds Averaged Navier-Stokes equations. The method used second-order finite differences, collocated grids, pressure-Poisson equation and four-stage Runge-Kutta scheme as key components of the solution method. A modified Baldwin-Lomax model is used for the turbulence closure. This paper presents a preliminary result of the computational study on turbulent flows past a ship model in drift motion. Calculations are carried out for a Series 60 $C_B=0.6$ ship model, for which detailed experimental data are available. The results of the present calculations are compared with the experimental data for hydrodynamic forces acting on the model as well as velocity distributions at longitudinal sections. Only fair agreements has been achieved. The computational results show the complex asymmetrical shear flow patterns including three-dimensional separations followed by formation of bilge vortices both in bow and stern regions.

• 한국전산유체공학회:학술대회논문집
• /
• 2002.10a
• /
• pp.35-40
• /
• 2002

This paper deals with a thrust generation on flapping-airfoil by dynamic stall. Dynamic stall refers to a series of complicated aerodynamic phenomena accompanied by a stall delay in unsteady motion. In most cases, once it occurs, the dynamic stall may lead to an abrupt fluctuation of aerodynamic forces. An inverse $k\acute{a}rm\acute{a}n$ vortex has been considered as a main reason for a thrust generation. In this paper, however, we have found out that a thrust is closely related to reduced frequency and leading edge vortex in addition to inverse Karman vortex. In order to certify our opinion, picking and plunging motions were calculated with the parameter of amplitude and frequency by using the unsteady, incompressible Navier-Stokes flow solver with a two-equation turbulence model. For more efficient computation, it is parallelized by MPI programming method.

• Journal of Ship and Ocean Technology
• /
• v.6 no.2
• /
• pp.37-50
• /
• 2002

In the present numerical simulation of viscous free surface flow around a ship, two-fluids in-compressible Reynolds-averaged Navier-Stokes equations with the standard $\textsc{k}-\varepsilon$turbulence model are discretized on a regular grid by using a finite volume method. A local level-set method is introduced for capturing the free surface movement and the influence of the viscous layer and dynamic boundary condition of the free surface are implicitly considered. Partial differential equations in the level-set method are discretized with second order ENO scheme and explicit Euler scheme in the space and time integration, respectively. The computational results for the Series-60 model with $C_B=0.6$ show a good agreement with the experimental data, but more validation studies for commercial complicated hull forms are necessary.