• Advances in aircraft and spacecraft science
• /
• v.5 no.3
• /
• pp.363-383
• /
• 2018

One-dimensional (1D) models of incompressible flows, can be of interest for many applications in which fast resolution times are demanded, such as fluid-structure interaction of flows in compliant pipes and hemodynamics. This work proposes a higher-order 1D theory for the flow-field analysis of incompressible, laminar, and viscous fluids in rigid pipes. This methodology is developed in the domain of the Carrera Unified Formulation (CUF), which was first employed in structural mechanics. In the framework of 1D modelling, CUF allows to express the primary variables (i.e., velocity and pressure fields in the case of incompressible flows) as arbitrary expansions of the generalized unknowns, which are functions of the 1D computational domain coordinate. As a consequence, the governing equations can be expressed in terms of fundamental nuclei, which are invariant of the theory approximation order. Several numerical examples are considered for validating this novel methodology, including simple Poiseuille flows in circular pipes and more complex velocity/pressure profiles of Stokes fluids into non-conventional computational domains. The attention is mainly focused on the use of hierarchical McLaurin polynomials as well as piece-wise nonlocal Lagrange expansions of the generalized unknowns across the pipe section. The preliminary results show the great advantages in terms of computational costs of the proposed method. Furthermore, they provide enough confidence for future extensions to more complex fluid-dynamics problems and fluid-structure interaction analysis.

• Journal of the Society of Naval Architects of Korea
• /
• v.40 no.5
• /
• pp.36-42
• /
• 2003

A CFD simulation technique has been developed to handle the unsteady body motion with large amplitude by use of overlapping multi-block grid system. The three-dimensional, viscous and incompressible flow around body is investigated by solving the Navier-Stokes equations, and the motion of body is represented by moving effect of the grid system. Composite grid system is employed in order to deal with both the body motion with large amplitude and the condition of numerical wave maker in convenience at the same time. The governing equations, Navier-Stokes (N-S) and continuity equations, are discretized by a finite volume method, in the framework of an O-H type boundary-fitted grid system (inner grid system including test model) and a rectangular grid system (outer grid system including simulation equipments for generation of wave environments). If this study, several flow configurations, such as an oscillating cylinder with large KC number, are studied in order to predict and evaluate the hydrodynamic forces. Furthermore, the motion simulation of a Series 60 model advancing in a uniform flow under the condition of enforced roll motion of angle 20$^{\circ}$ is performed in the developed numerical wave tank.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.10 no.5
• /
• pp.98-106
• /
• 2002

The airflow characteristics of an air-conditioning duct with multiple diffusers were investigated through one-dimensional analysis, CFD simulation and experimental measurement. One-dimensional program based on Bernoulli's equation and minor loss equations was developed in order to evaluate the air distribution rate at each diffuser. In CFD simulation, three-dimensional flow characteristics inside air-conditioning duct were computed for incompressible viscous flow, adopting the RNG k-$\xi$turbulence model. Also, in an effort to equalize the discharge flow rate at each outlet, the optimization procedure has been performed to obtain the optimum diffuser area. In this process, square of difference between maximum discharge rate and minimum discharge rate is used as an object function. Diffuser area and discharge velocity are established as constraints. After optimization process, determined design variables are applied again in CFD simulation and experiment to validate the optimized result by one-dimensional program. Comparison with the experimental data of airflow rate distribution showed that the developed program seems to be acceptable and can be useful design tool for an automotive air-conditioning duct in an initial design stage.

• 한국전산유체공학회:학술대회논문집
• /
• 1996.05a
• /
• pp.99-104
• /
• 1996

The jet pump is being used in many fields for several purposes because of its simple construction and easy operation. The characteristics of the geometrical variables, pressure gradient and velocity distribution of the jet pump are studied using the CFD technique. The flow calculations through a bended nozzle. a mixing chamber and a venturi are presented and phenomenological aspects are discussed. This study solve 3-D steady incompressible Navier-Stokes equations using the Iterative time marching scheme. The governing equations are differenced with 1st-order accurate backward difference scheme for the time derivatives and 3rd-order accurate QUICK scheme for the convective terms. The Mark-and-cell concept was applied efficiently to solve continuity equation, which is differenced 2nd-order accurate central differenced scheme. The 4th-order artificial damping is added to the continuity equation for numerical stability. A O-type of grid system is generated inside a nozzle and venturi of the jet pump. It has concluded that the results of present study properly agree with physical flow phenomena.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.22 no.11
• /
• pp.767-772
• /
• 2010

The incompressible viscous flow past a sphere under forced oscillation is numerically investigated at a Reynolds number of 300. The immersed boundary method is used to handle the sphere oscillating vertically to the streamwise direction. There are two important variables to characterize the oscillating state of a sphere. One is an oscillating amplitude normalized by the sphere diameter is set as a fixed number of 0.2. Another is the frequency ratio which is defined by $f_e/f_o$, where fe and fo are the excited frequency and the natural frequency of vortex shedding for the stationary sphere. In this study, three different frequency ratios of 0.8, 1.0 and 1.2 are considered. The results show a periodic flow with hairpin vortices shedding from upper and lower positions as well as vortical legs obliquely extended by oscillating motion of sphere. The enveloping vortical structure experience rupture twice in one period of oscillation. As the frequency of oscillation is increased, the vortical legs are getting shorter and eventually the hairpin vortices are much closer to the adjacent one.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.1 no.1
• /
• pp.1-12
• /
• 2009

Wing-in-Ground vehicles and aerodynamically assisted boats take advantage of increased lift and reduced drag of wing sections in the ground proximity. At relatively low speeds or heavy payloads of these craft, a flap at the wing trailing-edge can be applied to boost the aerodynamic lift. The influence of a flap on the two-dimensional NACA 4412 airfoil in viscous ground-effect flow is numerically investigated in this study. The computational method consists of a steady-state, incompressible, finite volume method utilizing the Spalart-Allmaras turbulence model. Grid generation and solution of the Navier-Stokes equations are completed using computer program Fluent. The code is validated against published experimental and numerical results of unbounded flow with a flap, as well as ground-effect motion without a flap. Aerodynamic forces are calculated, and the effects of angle of attack, Reynolds number, ground height, and flap deflection are presented for a split and plain flap. Changes in the flow introduced with the flap addition are also discussed. Overall, the use of a flap on wings with small attack angles is found to be beneficial for small flap deflections up to 5% of the chord, where the contribution of lift augmentation exceeds the drag increase, yielding an augmented lift-to-drag ratio.

• Journal of Ocean Engineering and Technology
• /
• v.15 no.4
• /
• pp.46-52
• /
• 2001

SMAC method, one of the numerical simulation techniques, is modified from the original MAC method for the time-dependent variation of fluid flows. The Navier-Stokes equations for incompressible time-dependent viscous flow are applied, and Also marker particles which present the visualization of fluid flows are used. In this study, two-dimensional numerical simulations of the overtopping are carried out by SMAC method, and the simulation results are visualized, In addition to, motion pictures are made for efficient visualization of the simulation results. This numerical simulation could also be applied to the design of coastal structures as dike and revetment.

• Journal of Advanced Research in Ocean Engineering
• /
• v.2 no.2
• /
• pp.61-66
• /
• 2016

This study presents the wave run-up height and depression depth around offshore cylindrical structures according to the wave period. The present study employs the volume of fluid method with the realizable turbulence model based on a commercial computational fluid dynamics software called the "STAR-CCM+" to simulate a 3D incompressible viscous two-phase turbulent flow. The present results for the wave run-up height and depression depth with regard to the wave period are compared with those of the relevant previous experimental and numerical studies.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1998.04a
• /
• pp.60-67
• /
• 1998

This paper deals with the development of a variable-node element and its application to the adaptive h-version mesh refinement-recovery for the incompressible viscous flow analysis. The element which has variable mid-side nodes can be used in generating the transition zone between the refined and unrefined elements and efficiently used for construction of a refined mesh without generating distorted elements. A modified Gaussian quadrature is needed to evaluate the element matrices due to the discontinuity of derivatives of the shape functions used for the element. The penalty function method which can reduce the number of independent variables is adopted for the purpose of computational efficiency and the selective reduced integration is carried out for the convection and pressure terms to preserve the stability of solution. For the economical analysis of transient problems, not only the mesh refinement but also the mesh recovery is needed. The numerical examples show that the optimal mesh for the finite element analysis of a wind around the structures can be obtained automatically by the proposed scheme.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 1994.10a
• /
• pp.55-63
• /
• 1994

The analysis of superplastic sheet forming process is studied by the use of the finite element method using a convected coordinate system and a skew boundary condition. In the formulation, the large inelastic behavior of the superplastic material is described as incompressible, nonlinear, viscous flow. The formulation is then approximated to the finite dimensional space with the use of membrane elements, which results in algebraic linear equations. In addition to the finite element formulation, a pressure cycle control algorithm is combined in the analysis for optimization of the forming time, which deals with the maximization of the strain rate sensitivity, the protection of the thickness reduction, the consistency of the desired strain rate and improvement of formability.