• 한국전산유체공학회:학술대회논문집
• /
• 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.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.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.03a
• /
• pp.251-255
• /
• 2008

Flow field around helicopter involves incompressible flow near the blade root and compressible flow at the blade tip. A problem occurs for low Mach number flow due to the stiffness of the governing equations. Time-derivative preconditioning techniques have been incorporated to reduce the stiffness that occurs at low speed region. The preconditioned form of the compressible Navier-Stokes and Euler equations is used. Computations are performed for the Caradonna-Tung's hovering and non-lifting forward flight case. Computational results are in good agreement with the experimental data.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.10a
• /
• pp.251-255
• /
• 2008

Flow field around helicopter involves incompressible flow near the blade root and compressible flow at the blade tip. A problem occurs for low Mach number flow due to the stiffness of the governing equations. Time-derivative preconditioning techniques have been incorporated to reduce the stiffness that occurs at low speed region. The preconditioned form of the compressible Navier-Stokes and Euler equations is used. Computations are performed for the Caradonna-Tung's hovering and non-lifting forward flight case. Computational results are in good agreement with the experimental data.

• Journal of computational fluids engineering
• /
• v.18 no.2
• /
• pp.56-65
• /
• 2013

Two compression ramp problems and an impinging shock problem are computed to investigate influence of turbulence models and eddy viscosity on the shock-wave / boundary layer interaction. A Navier-Stokes boundary layer generation code was applied to the generation of inflow boundary conditions. Computational results are validated well with the experimental data and effects of turbulence models are investigated. It is shown that the behavior of turbulence (eddy) viscosity directly affects both the extent of the separation and shock-wave positions over the separation.

• Journal of Ocean Engineering and Technology
• /
• v.4 no.2
• /
• pp.223-223
• /
• 1990

A circular cylinder is oscillated in th otherwise quiescent viscous fluid. Numerical analysis performed for this problem by using the fourth-order Runge-kutta method for the unsteady Navier-stokes equations. For K(Kelegan-Carpenter''s No.)=5, the flow developed symmetrically, while for K=10, it revealed random patterns. The coefficient of the rms force is overestimated by 20-30% compared with the experimental result.

• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 2006.11a
• /
• pp.113-120
• /
• 2006

In this paper, generation mechanisms of ocean freak waves are briefly introduced in the context of wave-current interaction phenomena. As an accurate and efficient numerical tool, the spectral element method is presented with general features and specific treatment for the wave-current interaction problem. The present model of the fluid motion is based on the Navier-Stokes equations incorporating a velocity-pressure formulation. In order to deal with the free surface motion, an Arbitrary Lagrangian-Eulerian (ALE) description is adopted. As an intermediate stage of development, solution procedure and characteristic aspects of the present modeling and numerical method features are addressed in detail, and numerical results for wave-current interaction is left as further study.

• Smart Structures and Systems
• /
• v.13 no.4
• /
• pp.637-657
• /
• 2014

This paper presents a dynamic analysis of three-dimensional beams. Structures made of functionally graded materials are considered. Several higher-order as well as classical theories are derived by means of a compact notation for the a-priori expansion order of the displacement field over the beam cross-section. The governing differential equations and boundary conditions are obtained in a condensed nuclear form that does not depend on the kinematic hypotheses. The problem is, then, exactly solved in space by means of a Navier-type solution, whereas time integration is performed by means of Newmark's solution scheme. Slender and short simply supported beams are investigated. Results are validated towards three-dimensional FEM results obtained via the commercial software ANSYS. Numerical investigations show that good accuracy can be obtained through the proposed formulation provided that the appropriate expansion order is considered.

• 한국전산유체공학회:학술대회논문집
• /
• 1995.04a
• /
• pp.5-29
• /
• 1995

A three-dimensional flow field induced by two trains passing by each other inside a tunnel is studied based on the numerical simulation of the three-dimensional compressible Euler/Navier-Stokes equations formulated in the finite difference approximation. Domain decomposition method with the FSA(fortified solution algorithm) interface scheme is used to treat this moving-body problem. The computed resluts show basic characteristic of the flow field created when two trains passing by each other. History of the pressure distributions and the aerodynamic forces acting on the trains are mailnly discussed. The results indicate that the phenomenon is complicated due to the interaction of the flow induced by two trains. Strong side force occurs between the two trains when the front portion of the opposite train passes by. It fluctuates rapidly and maximum suction force occurs when two trains are aligned side by side. The results also indicate the effectiveness of the present numerical method for moving boundary problems.

• 한국전산유체공학회:학술대회논문집
• /
• 2005.04a
• /
• pp.119-122
• /
• 2005

In this paper a method for analysing aerodynamic characteristics of a rotating parachute in steady descending motion is presented. Because of a complex geometric configuration of the parachute associated with side vents and discontinuous skirts, special procedure was adopted th handle the geometry in the analysis. A panel method was successfully applied to the present problem and yielded good results using above procedure. A CFD code using the full Navier-Stokes equations was also applied and provided good results. Parachute free drop and wind tunnel tests were performed to determine the fully developed canopy configuration and aerodynamic characteristics. The method can be used for optimizing the parachute size and side vent configurations in the design period.