• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.39 no.6
• /
• pp.487-497
• /
• 2011

Numerical analysis has been performed to investigate the rotor speed and pitch range variations when the airspeed is increased in autorotation. Transient Simulation Method(TSM) was used to obtain the steady states of autorotation. The rotor blade was analyzed by the two-dimensional compressible Navier-Stokes solver in order to adapt to the airspeed increase and the results were used in the transient simulation method. Meanwhile, the Pitt/Peters inflow theory was used to supply the induced velocity fields. For the prescribed torque equilibrium state, the combinations of velocity, shaft angle, and pitch angle were produced to investigate the rotor speeds and variable ranges. The rotor tip Mach number and rotor speed were correlated and the trim range of pitch angle was observed with respect to the shaft angle decrease.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.32 no.7
• /
• pp.13-18
• /
• 2004

This paper describes an efficient and robust optimization method for helicopter rotor airfoil design in forward flight. Navier-Stokes analysis was employed to compute the dynamic response of an airfoil, which simulates the unsteady rotor flow-field in forward flight. The optimization system consists of two categories; Response Surface Method to construct the response surface model based on D-optimal 3-level factorial design, and Genetic Algorithm to obtain the optimum solution of a defined objective function including penalty terms of constraints. The influence of design variables and their interactions on the aerodynamic performance was examined through the optimization process.

• Transactions of the Society of Information Storage Systems
• /
• v.5 no.1
• /
• pp.19-35
• /
• 2009

The present work is an experimental and analytical study on a flexible disk rotating close to a rigid rotating disk in open air. In the analytical study, the air flow in the gap between the flexible disk and the rigid disk is modeled using Navier-Stokes and continuity equations while the flexible disk is modeled using the linear plate theory. The flow equations are discretized using the cell centered finite volume method (FVM) and solved numerically with semi-implicit pressure-linked equations (SIMPLE algorithm). The spatial terms in the disk equation are discretized using the finite difference method (FDM) and the time integration is performed using fourth-order Runge-Kutta method. An experimental test-rig is designed to investigate the dynamics of the flexible disk when rotating close to a co-rotating, a counter-rotating and a fixed rigid disk, which works as a stabilizer. The effects of rotational speed, initial gap height and inlet-hole radius on the flexible disk displacement and its vibration amplitude are investigated experimentally for the different types of stabilizer. Finally, the analytical and experimental results are compared.

• Ocean Systems Engineering
• /
• v.1 no.1
• /
• pp.31-57
• /
• 2011

Impact pressure due to sloshing is of great concern for the ship owners, designers and builders of the LNG carriers regarding the safety of LNG containment system and hull structure. Sloshing of LNG in partially filled tank has been an active area of research with numerous experimental and numerical investigations over the past decade. In order to accurately predict the sloshing impact load, a new numerical method was developed for accurate resolution of violent sloshing flow inside a three-dimensional LNG tank including wave breaking, jet formation, gas entrapping and liquid-gas interaction. The sloshing flow inside a membrane-type LNG tank is simulated numerically using the Finite-Analytic Navier-Stokes (FANS) method. The governing equations for two-phase air and water flows are formulated in curvilinear coordinate system and discretized using the finite-analytic method on a non-staggered grid. Simulations were performed for LNG tank in transverse and longitudinal motions including horizontal, vertical, and rotational motions. The predicted impact pressures were compared with the corresponding experimental data. The validation results clearly illustrate the capability of the present two-phase FANS method for accurate prediction of impact pressure in sloshing LNG tank including violent free surface motion, three-dimensional instability and air trapping effects.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.7 no.4
• /
• pp.566-576
• /
• 1995

A Computational code has been developed in order to design axial fans by the numerical optimization techniques incorporated with flow analysis code solving three-dimensional Navier-Stokes equation. The steepest descent method and the conjugate gradient method are used to look for the search direction in the design space, and the golden section method is used for one-dimensional search. To solve the constrained optimization problem, sequential unconstrained minimization technique, SUMT, is used with imposed quadratic extended interior penalty functions. In the optimization of two-dimensional cascade design, the ratio of drag coefficient to lift coefficient is minimized by the design variables such as maximum thickness, maximum ordinate of camber and chord wise position of maximum ordinate. In the application of this numerical optimization technique to the design of an axial fan, the efficiency is maximized by the design variables related to the sweep angle distributed by quadratic function along the hub to tip of fan.

• Wind and Structures
• /
• v.17 no.3
• /
• pp.263-274
• /
• 2013

A numerical technique for fluid-structure interaction, which is based on the finite element method (FEM) and computational fluid dynamics (CFD), was developed for application to an industrial chimney equipped with a pendulum tuned mass damper (TMD). In order to solve the structural problem, a one-dimensional beam model (Navier-Bernoulli) was considered and, for the dynamical problem, the standard second-order Newmark method was used. Navier-Stokes equations for incompressible flow are solved in several horizontal planes to determine the pressure in the boundary of the corresponding cross-section of the chimney. Forces per unit length were obtained by integrating the pressure and are introduced in the structure using standard FEM interpolation techniques. For the fluid problem, a fractional step scheme based on a second order pressure splitting has been used. In each fluid plane, the displacements have been taken into account considering an Arbitrary Lagrangian Eulerian approach. The stabilization of convection and diffusion terms is achieved by means of quasi-static orthogonal subscales. For each period of time, the fluid problem was solved and the geometry of the mesh of each fluid plane is updated according to the structure displacements. Using this technique, along-wind and across-wind effects have been properly explained. The method was applied to an industrial chimney in three scenarios (with or without TMD and for different damping values) and for two wind speeds, showing different responses.

• Structural Engineering and Mechanics
• /
• v.81 no.4
• /
• pp.523-527
• /
• 2022

The objective of this study is to simulate the flow in pipes with various boundary conditions. Free-pressure fluid model, is used in the pipe based on Navier-Stokes equation. The models are solved by using the numerical method. A problem called "stability of pipes" is used in order to compare frequency and critical fluid velocity. When the initial conditions of problem satisfied the instability conditions, the free-pressure model could accurately predict discontinuities in the solution field. Employing nonlinear strains-displacements, stress-strain energy method the governing equations were derived using Hamilton's principal. Differential quadrature method (DQM) is used for obtaining the frequency and critical fluid velocity. The results of this paper are analyzed by hyperbolic numerical method. Results show that the level of numerical diffusion in the solution field and the range of well-posedness are two important criteria for selecting the two-fluid models. The solutions for predicting the flow variables is approximately equal to the two-pressure model 2. Therefore, the predicted pressure changes profile in the two-pressure model is more consistent with actual physics. Therefore, in numerical modeling of gas-liquid two-phase flows in the vertical pipe, the present model can be applied.

• The KIPS Transactions:PartA
• /
• v.16A no.6
• /
• pp.419-426
• /
• 2009

To realistically simulate fluid, the Navier-Stokes equations are generally used. Solving these Navier-Stokes equations on the Eulerian framework, the non-linear advection terms invoke heavy computation and thus Semi-Lagrangian methods are used as an approximated way of solving them. In the Semi-Lagrangian methods, the locations of advection sources are traced and the physical values at the traced locations are interpolated. In the case of Stam's method, there are relatively many chances of numerical losses, and thus there have been efforts to correct these numerical errors. In most cases, they have focused on the numerical interpolation processes, even simultaneously using particle-based methods. In this paper, we propose a new approach to reduce the numerical losses, through improving the tracing method during the advection calculations, without any modifications on the Eulerian framework itself. In our method, we trace the grids with the velocities which will let themselves to be moved to the current target position, differently from the previous approaches, where velocities of the current target positions are used. From the intuitive point of view, we adopted the simple physical observation: the physical quantities at a specific position will be moved to the new location due to the current velocity. Our method shows reasonable reduction on the numerical losses during the smoke simulations, finally to achieve real-time processing even with enhanced realities.

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

The moment-of-fluid (MOF) method is a new volume-tracking method that accurately treats evolving material interfaces. Based on the moment data (volume and centroid) for each material, the material interfaces are reconstructed with second-order spatial accuracy in a strictly conservative manner. The MOF method is coupled with a stabilized finite element incompressible Navier-Stokes solver for two fluids, namely water and air. The effectiveness of the MOF method is demonstrated with a free-surface dam-break problem.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.32 no.6B
• /
• pp.379-385
• /
• 2012

Recent development of computing power and theoretical advances in computational fluid dynamics have made possible numerical simulations of water waves with full Navier-Stokes equations. In this study, an internal wave maker using the mass source function approach was combined with the level set finite element method for generation of waves. The model is first applied to the two-dimensional linear wave generation and propagation. Then, it is applied to the three-dimensional simulation of the same problem. To effectively utilize computational resources and enhance the speed of execution, parallel algorithms are developed and applied for the three-dimensional problem. The results of numerical simulations are compared with theoretical values and good agreements are observed.