• Journal of the Society of Naval Architects of Korea
• /
• v.51 no.5
• /
• pp.409-418
• /
• 2014

Simulation of flow past a complex marine structure requires a fine resolution in the vicinity of the structure, whereas a coarse resolution is enough far away from it. Therefore, a lot of grid cells may be wasted, when a simple Cartesian grid system is used for an Immersed Boundary Method (IBM). To alleviate this problems while maintaining the Cartesian frame work, we adopted an Adaptive Mesh Refinement (AMR) scheme where the grid system dynamically and locally refines as needed. In this study, We implemented a moving IBM and an AMR technique in our basic 3D incompressible Navier-Stokes solver. A Volume Of Fluid (VOF) method was used to effectively treat the free surface, and a recently developed Lagrangian Dynamic Subgrid-scale Model (LDSM) was incorporated in the code for accurate turbulence modeling. To capture vortex induced vibration accurately, the equation for the structure movement and the governing equations for fluid flow were solved at the same time implicitly. Also, We have developed an interface by using AutoLISP, which can properly distribute marker particles for IBM, compute the geometrical information of the object, and transfer it to the solver for the main simulation. To verify our numerical methodology, our results were compared with other authors' numerical and experimental results for the benchmark problems, revealing excellent agreement. Using the verified code, we investigated the following cases. (1) simulating flow around a floating sphere. (2) simulating flow past a marine structure.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.37 no.2
• /
• pp.113-120
• /
• 2009

To numerically simulate aerodynamics of rotor-airframe interaction in a rigorous manner, we need to solve the Navier-Stokes system for a rotor-airframe combination in a single computational domain. This imposes a computational burden since rotating blades and a stationary body have to be simultaneously dealt with. An efficient alternative is a momentum source method in which the action of rotor is approximated as momentum source in a stationary mesh system built around the airframe. This makes the simulation much easier. The magnitude of the momentum source is usually evaluated by the blade element theory, which often results in a poor accuracy. In the present work, we evaluate the momentum source from the simulation data by using the Navier-Stokes equations only for a rotor system. Using this data, we simulated the time-averaged steady rotor-airfame interaction and developed the unsteady rotor-airframe interaction. Computations were carried out for the simplified rotor-airframe model (the Georgia Tech configuration) and the results were compared with experimental data. The results were in good agreement with experimental data, suggesting that the present approach is a usefull method for rotor-airframe interaction analysis.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.6 no.2
• /
• pp.219-235
• /
• 2014

With the development of the technology of underwater moving bodies, the need for developing the knowledge of surface effect interaction of free surface and underwater moving bodies is increased. Hence, the two-phase flow is a subject which is interesting for many researchers all around the world. In this paper, the non-linear free surface deformations which occur during the water-exit of a circular cylinder due to its buoyancy are solved using finite volume discretization based code, and using Volume of Fluid (VOF) scheme for solving two phase flow. Dynamic mesh model is used to simulate dynamic motion of the cylinder. In addition, the effect of cylinder mass in presence of an external force is studied. Moreover, the oblique exit and entry of a circular cylinder with two exit angles is simulated. At last, water-exit of a circular cylinder in six degrees of freedom is simulated in 3D using parallel processing. The simulation errors of present work (using VOF method) for maximum velocity and height of a circular cylinder are less than the corresponding errors of level set method reported by previous researchers. Oblique exit shows interesting results; formation of waves caused by exit of the cylinder, wave motion in horizontal direction and the air trapped between the waves are observable. In 3D simulation the visualization of water motion on the top surface of the cylinder and the free surface breaking on the front and back faces of the 3D cylinder at the exit phase are observed which cannot be seen in 2D simulation. Comparing the results, 3D simulation shows better agreement with experimental data, specially in the maximum height position of the cylinder.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.14 no.5
• /
• pp.32-44
• /
• 2010

A two-dimensional thermal response and ablation analysis code for predicting charring material ablation and shape change on solid rocket nozzle is presented. The thermogravimetric analysis (TGA) techniques have been used to characterize the thermal decomposition constants for Arrhenius parameters. Two heterogeneous reactions involving carbon and the oxidizing species of $H_2O$ and $CO_2$ are considered and determined by Zvyagin's ablation model and kinetic constants. The moving boundary problem and mesh moving are solved by remeshing-rezoning method in MSC-Marc-ATAS program. The difference between the calculated and experimental value of char and ablation thickness is up to 20%. For the performance prediction of thermal protection systems, this method will be integrated with a three-dimensional finite-element thermal and structure analysis code through the real time sensing of in-depth temperature and heat flux.

• International Journal of Highway Engineering
• /
• v.10 no.4
• /
• pp.213-224
• /
• 2008

Flexible pavement responses to vehicular loading, such as critical stresses and strains, in each pavement layer, could be predicted by the multilayered elastic analysis. However, multilayered elastic theory suffers from major drawbacks including spatial dimension of a numerical model, material properties considered in the analysis, boundary conditions, and ill-presentation of tire-pavement contact shape and stresses. To overcome these shortcomings, three-dimensional finite element (3D FE) models are developed and numerical analyses are conducted to calculate pavement responses to moving load in this study. This paper introduces a methodology for an effective 3D FE to simulate flexible pavement structure. It also discusses the mesh development and boundary condition analysis. Sensitivity analyses of flexible pavement response to loading are conducted. The infinite boundary conditions and time-dependent history of calculated pavement responses are considered in the analysis. This study found that the outcome of 3D FE implicit dynamic analysis of flexible pavement that utilizes appropriate boundary conditions, continuous moving load, viscoelastic hot-mix asphalt model is comparable to field measurements.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.25 no.2
• /
• pp.139-148
• /
• 2012

The MLS(Moving Least Squares) Difference Method is a numerical scheme that combines the MLS method of Meshfree method and Taylor expansion involving not numerical quadrature or mesh structure but only nodes. This paper presents an dynamic algorithm of MLS difference method for solving transient solid mechanics problems. The developed algorithm performs time integration by using Newmark method and directly discretizes strong forms. It is very convenient to increase the order of Taylor polynomial because derivative approximations are obtained by the Taylor series expanded by MLS method without real differentiation. The accuracy and efficiency of the dynamic algorithm are verified through numerical experiments. Numerical results converge very well to the closed-form solutions and show less oscillation and periodic error than FEM(Finite Element Method).

• Journal of Mechanical Science and Technology
• /
• v.17 no.6
• /
• pp.888-895
• /
• 2003

Combustion chamber crevices in SI engines are identified as the largest contributors to the engine-out hydrocarbon emissions. The largest crevice is the piston ring-pack crevice. A numerical simulation method was developed, which would allow to predict and understand the oxidation process of piston crevice hydrocarbons. A computational mesh with a moving grid to represent the piston motion was built and a 4-step oxidation model involving seven species was used. The sixteen coefficients in the rate expressions of 4-step oxidation model are optimized based on the results from a study on the detailed chemical kinetic mechanism of oxidation in the engine combustion chamber. Propane was used as the fuel in order to eliminate oil layer absorption and the liquid fuel effect. Initial conditions of the burned gas temperature and in-cylinder pressure were obtained from the 2-zone cycle simulation model. And the simulation was carried out from the end of combustion to the exhaust valve opening for various engine speeds, loads, equivalence ratios and crevice volumes. The total hydrocarbon (THC) oxidation in the crevice during the expansion stroke was 54.9% at 1500 rpm and 0.4 bar (warmed-up condition). The oxidation rate increased at high loads, high swirl ratios, and near stoichiometric conditions. As the crevice volume increased, the amount of unburned HC left at EVO (Exhaust Valve Opening) increased slightly.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.14 no.4
• /
• pp.39-48
• /
• 2006

In this research, in order to study the spray, combustion, and emission characteristics of the common rail DME engine, the target engine was disassembled, and 3D CAD file was constructed using a 3D measurement machine and a rapid prototyping machine. Using the obtained 3D geometry, fine moving meshes are generated, and three dimensional non-steady turbulence flow field and combustion phenomenon including spray were numerically analyzed. As a result, IMEP of DME and diesel in medium and high speed revolution showed similar performance. As the DME fuel start to burn in spray area, the vaporized fuel rapidly spreads squish area in low speed revolution. In the case of DME engine, CO and NOx are relatively consistent with experiment results. It was found that the break-up, evaporation, collision model of DME fuel need to be properly adjusted through matching the characteristics of fuel and injector for further improvement.

• The KSFM Journal of Fluid Machinery
• /
• v.16 no.1
• /
• pp.40-46
• /
• 2013

Recently, the duct system with a cross flow fan was used to improve the ventilation in various industrial fields. For the efficient ventilation, it is necessary to design the duct system based on the flow characteristics around the cross flow fan. In the present study, the flow characteristics around a cross flow fan in the ventilation duct were predicted by using the moving mesh and sliding interface techniques for the rotation of blades. To design the duct system with the high performance of ventilation, the CFD simulations were repeated with the revised duct model based on the DOE. With the numerical results of flow rate through the ventilation duct with various geometric parameters, the optimized geometry of ventilation duct to maximize the flow rate was obtained by using the Kriging approximation method. From the performance curves of cross flow fan in the original and optimized models of ventilation duct, it was observed that the flow rate through the optimized model is about 16 percent larger than that through the original model.

• Journal of the Korean Society for Precision Engineering
• /
• v.29 no.10
• /
• pp.1043-1049
• /
• 2012

The FEM analysis of machine tools is the general analysis process to evaluate machine performance in the industry for a long time. Despite advances in FEM software, because of difficult simplicity of CAD drawing, little experience of joints stiffness modeling and troublesome manual contact area divide for bindings, the industry designers think the FEM analysis is still an area of FEM analysis expert. In this paper, the automation of modeling process with simplicity of drawing, modeling of joints and contact area divide is aimed at easy FEM analysis to enlarge utilization of a virtual machine tools. In order to verify the effects of modeling automation, a slant bed type model with tilting table is analyzed. The results show FEM modeling automation method only needed 45 minutes to complete the whole modeling process, while manual modeling method requires almost one month with 8200 calculations for coordinate transformations and stiffness data input.