• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 한국소성가공학회 2004년도 추계학술대회논문집
• /
• pp.75-78
• /
• 2004

A mesh generation algorithm adapted to the mesh density map using the Delaunay mesh generation technique is developed. In the finite element analyses of the forging processes, the numerical error increases as the process goes on because of discrete property of the finite elements or severe distortion of elements. Especially, in the region where stresses and strains are concentrated, the numerical discretization error will be highly increased. However, it is too time consuming to use a uniformly fine mesh in the whole domain to reduce the expected numerical error. Therefore, it is necessary to construct locally refined mesh at the region where the error is concentrated such as at the die corner. In this study, the point insertion algorithm is used and the mesh size is controlled by moving nodes to optimized positions according to a mesh density map constructed with a posteriori error estimation. An optimization technique is adopted to obtain a good position of nodes. And optimized smoothing techniques are also adopted to have smooth distribution of the mesh and improve the mesh element quality.

• International Journal of Aeronautical and Space Sciences
• /
• 제12권4호
• /
• pp.318-330
• /
• 2011

This paper investigates the convergence characteristics of the modified artificial compressibility method proposed by Turkel. In particular, a focus is mode on the convergence characteristics due to variation of the preconditioning factor (${\alpha}_u$) and the artificial compressibility (${\beta}$) in conjunction with an upwind method. For the investigations, a code using the modified artificial compressibility is developed. The code solves the axisymmetric incompressible Reynolds averaged Navier-Stokes equations. The cell-centered finite volume method is used in conjunction with Roe's approximate Riemann solver for the inviscid flux, and the central difference discretization is used for the viscous flux. Time marching is accomplished by the approximated factorization-alternate direction implicit method. In addition, Menter's k-${\omega}$ shear stress transport turbulence model is adopted for analysis of turbulent flows. Inviscid, laminar, and turbulent flows are solved to investigate the accuracy of solutions and convergence behavior in the modified artificial compressibility method. The possible reason for loss of robustness of the modified artificial compressibility method with ${\alpha}_u$ >1.0 is given.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• 제38권8호
• /
• pp.790-797
• /
• 2010

Mode transition of the axi-symmetric screech tone in the low supersonic Mach number range from 1.0 to 1.20 is numerically analyzed. The axi-symmetric Navier-Stokes equations and the k-e turbulence model are solved in the cylindrical coordinate system. The dispersion-relation-preserving(DRP) scheme is applied for space discretization and the optimized four levels marching method are used for time integration. At low supersonic Mach numbers with an axi-symmetric A1 mode in the simulation, it is shown that acoustic propagation due to the nonlinear effects is seen in the lateral direction and the screech tone frequency is the same as the vortex passing frequency due to the generation of intense large-scale vortical motions.

• Proceedings of the Korean Society For Composite Materials Conference
• /
• 한국복합재료학회 2000년도 추계학술발표대회 논문집
• /
• pp.228-231
• /
• 2000

The active vibration control of laminated composite shell has been performed with the optimized sensor/actuator system. PVDF film is used fur the material of sensor/actuator. Finite element method is utilized to model the whole structure including the piezoelectric sensor/actuator system, The distributed selective modal sensor/actuator system is established to prevent the adverse effect of spillover. In the finite element discretization process, the nine-node shell element with five nodal degrees of freedoms is used. Electrode patterns and lamination angles of sensor/actuator are optimized using genetic algorithm. Sensor is designed to minimize the observation spillover, and actuator is designed to minimize the system energy of the control modes under a given initial condition. Modal sensor/actuator profiles are optimized for the first and the second modes suppression of singly curved cantilevered composite shell structure. Discrete LQG method is used as a control law. The real time vibration control with profile optimized sensor/actuator system has been performed. Experimental result shows successful performance of the integrated structure for the active vibration control.

• Wind and Structures
• /
• 제17권1호
• /
• pp.1-19
• /
• 2013

A multiscale finite element method is applied to the Spalart-Allmaras turbulence model based detached-eddy simulation (DES). The multiscale arises from a decomposition of the scalar field into coarse (resolved) and fine (unresolved) scales. It corrects the lack of stability of the standard Galerkin formulation by modeling the scales that cannot be resolved by a given spatial discretization. The stabilization terms appear naturally and the resulting formulation provides effective stabilization in turbulent computations, where reaction-dominated effects strongly influence near-wall predictions. The multiscale DES is applied in the context of high-Reynolds flow over the Commonwealth Advisory Aeronautical Council (CAARC) standard tall building model, for both uniform and turbulent inflows. Time-averaged pressure coefficients on the exterior walls are compared with experiments and it is demonstrated that DES is able to resolve the turbulent features of the flow and accurately predict the surface pressure distributions under atmospheric boundary layer flows.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• 제28권11호
• /
• pp.1302-1314
• /
• 2004

Turbulent heat transfer to $CO_2$ at supercritical pressure flowing in vertical tubes is investigated using direct numerical simulation (DNS). A conservative space-time discretization scheme for variable-density flows at low Mach numbers is adopted in the present study to treat steep variations of fluid properties at supercritical pressure just above the thermodynamic critical point. The fluid properties at these conditions are obtained using PROPATH and used in the form of tables in the simulations. The buoyancy influence induced by strong variation of density across the pseudo-critical temperature proved to play a major role in turbulent heat transfer at supercritical state. Depending on the degree of buoyancy influence, turbulent heat transfer may be enhanced or significantly deteriorated, resulting in local hot spots along the heated surface. Based on the results of the present DNS combined with theoretical considerations, the physical mechanism of this local heat transfer deterioration is elucidated.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 한국추진공학회 2004년도 제22회 춘계학술대회논문집
• /
• pp.426-432
• /
• 2004

Large-Eddy Simulation (LES) is applied for the simulation of compressible flat plate boundary with Reynolds number up to 5 X 10$^{5}$ . Numerical examples include shock/boundary layer interaction and boundary layer transition, aiming future application to the analysis of transonic fan/compressor cascades. The present LES code uses hybrid com-pact/WENO scheme for the spatial discretization and compact diagonalized implicit scheme for the time integration. The present code successfully predicted the bypass transition of subsonic boundary layer. As for supersonic turbulent boundary layer, mean and fluctuation velocity of the attached boundary, as well as the evolution of the friction coefficient and the displacement thickness both upstream and downstream of the separation region are all in good agreement with experiment. The separation point also agreed with the experiment. In the simulation of the shock/laminar boundary layer interaction, the dependence of the transition upon the shock strength is reproduced qualitatively, but the extent of the separation region is overpredicted. These numerical examples show that LES can predict the behavior of boundary layer including transition and shock interaction, which are hardly managed by the conventional Reynolds-averaged Navier-Stokes approach, although there needs to be more effort before achieving quantitative agreement.

• Journal of the Society of Naval Architects of Korea
• /
• 제36권4호
• /
• pp.37-47
• /
• 1999

A Finite Volume Method for the discretization of the two-dimensional incompressible Navier-Stokes equation is used to analyse water entry & exit problems in a generalized coordinate system. The free-surface deformations generated by the water entry or exit of a rigid body are simulated by the Level-Set scheme[11]. In the water entry problems for a wedged section and a flared-ship section, the calculation results of water impact force are compared with the experimental results[5] and the time varying free-surface deformations and flow characteristics of the water exit of a cylinder are investigated.

• Advances in aircraft and spacecraft science
• /
• 제5권3호
• /
• pp.349-362
• /
• 2018

A non-linear numerical simulation technique for predicting the unsteady performances of an airbreathing engine is developed. The study focuses on the simulation of integrated propulsion systems, where a closer coupling is needed between the airframe and the engine dynamics. In fact, the solution of the fully unsteady flow governing equations, rather than a lumped volume gas dynamics discretization, is essential for modeling the coupling between aero-servoelastic modes and engine dynamics in highly integrated propulsion systems. This consideration holds for any propulsion system when a full separation between the fluid dynamic time-scale and engine transient cannot be appreciated, as in the case of flow instabilities (e.g., rotating stall, surge, inlet unstart), or in case of sudden external perturbations (e.g., gas ingestion). Simulations of the coupling between external and internal flow are performed. The flow around the nacelle and inside the engine ducts (i.e., air intakes, nozzles) is solved by CFD computations, whereas the flow evolution through compressor and turbine bladings is simulated by actuator disks. Shaft work balance and rotor dynamics are deduced from the estimated torque on each turbine/compressor blade row.

• Journal of the Society of Naval Architects of Korea
• /
• 제38권3호
• /
• pp.1-13
• /
• 2001

RANS solver has been developed to solve the flows past a ship with a drift angle. The solver employs a finite volume method for the spatial discretization and Euler implicit method for the time integration. Turbulent flows are simulated by Spalart-Allmaras one-equation model. Developed solver is applied to analyze the hydrodynamic forces and flows of two tankers with a same forebody but different afterbodies. The computed flows and hydrodynamic forces are compared with the measured flows and captive model test data. The computed results show good agreements with experimental data and show clearly the effects of stern hull form on the hydrodynamic forces and the flows.