• Journal of computational fluids engineering
• /
• v.14 no.4
• /
• pp.78-85
• /
• 2009

Cavitating flow simulation is of practical importance for many engineering systems, such as marine propellers, pump impellers, nozzles, injectors, torpedoes, etc. The present work has focused on the simulation of cavitating flow past cylinders with strong side flows. The governing equation is the Navier-Stokes equation based on the homogeneous mixture model. The momentum and energy equation is in the mixture phase while the continuity equation is solved liquid and vapor phase, separately. An implicit dual time and preconditioning method are employed for computational analysis. For the code validation, the results from the present solver have been compared with experiments and other numerical results. A fairly good agreement with the experimental data and other numerical results have been obtained. After the code validation, the strong side flow was applied to include the wake flow effects of the submarine or ocean tide.

• Journal of the Society of Naval Architects of Korea
• /
• v.44 no.2 s.152
• /
• pp.130-138
• /
• 2007

A coupled variational equation for fluid-structure interaction (FSI) problems is derived from a steady state Navier-Stokes equation for incompressible Newtonian fluid and an equilibrium equation for geometrically nonlinear structures. For a fully coupled FSI formulation, between fluid and structures, a traction continuity condition is considered at interfaces where a no-slip condition is imposed. Under total Lagrange formulation in the structural domain, finite rotations are well described by using the second Piola-Kirchhoff stress and Green-Lagrange strain tensors. An adjoint shape design sensitivity analysis (DSA) method based on material derivative approach is applied to the FSI problem to develop a shape design optimization method. Demonstrating some numerical examples, the accuracy and efficiency of the developed DSA method is verified in comparison with finite difference sensitivity. Also, for the FSI problems, a shape design optimization is performed to obtain a maximal stiffness structure satisfying an allowable volume constraint.

• Journal of computational fluids engineering
• /
• v.16 no.4
• /
• pp.7-13
• /
• 2011

A comparative analysis of thermal models in the lattice Boltzmann method(LBM) for the simulation of laminar natural convection in a square cavity is presented. A HYBRID method, in which the thermal equation is solved by the Navier-Stokes equation method while the mass and momentum conservation are resolved by the lattice Boltzmann method, is introduced and its merits are explained. All the governing equations are discretized on a cell-centered, non-uniform grid using the finite-volume method. The convection terms are treated by a second-order central-difference scheme with a deferred correction method to ensure stability of the solutions. The HYBRID method and the double-population method are applied to the simulation of natural convection in a square cavity and the predicted results are compared with the benchmark solutions given in the literatures. The predicted results are also compared with those by the conventional Navier-Stokes equation method. In general, the present HYBRID method is as accurate as the Navier-Stokes equation method and the double-population method. The HYBRID method shows better convergence and stability than the double-population method. These observations indicate that this HYBRID method is an efficient and economic method for the simulation of incompressible fluid flow and heat transfer problem with the LBM.

• Journal of applied mathematics & informatics
• /
• v.16 no.1_2
• /
• pp.383-405
• /
• 2004

We study the incomprssible Navier Stokes equations for the flow inside contraction geometry. The governing equations are expressed in the vorticity-stream function formulations. A rectangular computational domain is arised by elliptic grid generation technique. The numerical solution is based on a technique of automatic numerical generation of acurvilinear coordinate system by transforming the governing equation into computational plane. The transformed equations are approximated using central differences and solved simultaneously by successive over relaxation iteration. The time dependent of the vorticity equation solved by using explicit marching procedure. We will apply the technique on several irregular-shapes.

• Journal of Nuclear Fuel Cycle and Waste Technology(JNFCWT)
• /
• v.3 no.2
• /
• pp.77-84
• /
• 2005

A voloxidizer for a hot cell demonstration, that handles spent fuels of a high radiation level in a limited space should be small and spent fuel powders should not be dispersed out of the equipment involved. In this study a density rate equation as well as the Stokes'equation has been proposed in order to obtain the theoretical terminal velocity of powders. The terminal velocity of U$_{3}$O$_{8}$ has been predicted by using the terminal velocity of SiO$_{2}$, and then determination has been the optimum air flow rate which is able to prevent powders from scattering. An equation which has shown a relationship between theoretical terminal velocities of U$_{3}$O$_{8}$ and SiO$_{2}$ has been derived with the help of the Stokes'equation, and then an experimental verification made for the theoretical Stokes' equation of SiO$_{2}$ by means of an experimental device made of acryl. The theoretical terminal velocity based on the proposed density rate equation has been verified by detecting U$_{3}$O$_{8}$ powders in a filter installed in the mock-up voloxidizer. As the results, the optimum air flow rates seem to be 20 LPM by the Stokes'equation while they are 14.5 L/min by the density rate equation. At the experiments with the mock-up voloxidizer, a trace amount of U$_{3}$O$_{8}$ seems to be detectable at the air flow rate of 14.5 L/min by the density rate equation, but U$_{3}$O$_{8}$ powders of 7$\mu$m diameter seem detectable at the air flow rate of 20 L/min by the Stokes'equation. It is revealed that 14.5 L/min is the optimum air flowe rate which is capable of preventing U$_{3}$O$_{8}$ powders from scattering in the UO$_{2}$ voloxidizer and the proposed density rate equation is proper to calculate the terminal velocity of U$_{3}$O$_{8}$ powders.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.24 no.3
• /
• pp.9-16
• /
• 1987

In this paper, the fluid flow in the two-dimensional tank is analyzed by the Finite Difference Method. The Navier-Stokes equation is modified for the tank fixed coordinate system. For the treatment of the free surface, the Volume of Fluid Method by Hirt and Nichols is adopted. The continuity equation and the Poisson equation which is derived from the Navier-Stokes equation to find the pressure are solved by the Successive-Line-Overrelaxation Method. The comparison of the calculated results with experimental data show a favorable agreement. The fluid flow in the two-dimensional tank can be predicted reasonably before the free surface reaches breaking by this numerical method.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.11 no.2
• /
• pp.157-163
• /
• 2003

This work was undertaken for the numerical analysis of defrosting phenomena of automobile windshield. To analyze the defrost, the flow and temperature field of cabin interior, heat transfer through the windshield glass, and phase change of the frost should be analyzed simultaneously. The flow field was obtained by solving the 3-D unsteady Navier-Stokes equation and the temperature field was computed by energy equation. The phase-change process of Stefan problem was solved by enthalpy method. For code validation, the temperature field of the driven cavity was calculated. The result of calculation shows a good agreement with the other numerical results. Then, the present code was applied to the defrosting analysis of a real automobile and, also, a good agreement with experiment was obtained.

• Journal of Ocean Engineering and Technology
• /
• v.36 no.2
• /
• pp.101-107
• /
• 2022

Dean (1965) proposed the use of the root mean square error (RMSE) in the dynamic free surface boundary condition (DFSBC) and kinematic free-surface boundary condition (KFSBC) as an error evaluation criterion for wave theories. There are well known wave theories with RMSE more than 1%, such as Airy theory, Stokes theory, Dean's stream function theory, Fenton's theory, and trochodial theory for deep-water waves. However, none of them can be applied for deep-water breaking waves. The purpose of this study is to provide a closed-form solution for deep-water waves with RMSE less than 1% even for breaking waves. This study is based on a previous study (Shin, 2016), and all flow fields were simplified for deep-water waves. For a closed-form solution, all Fourier series coefficients and all related parameters are presented with Newton's polynomials, which were determined by curve fitting data (Shin, 2016). For verification, a wave in Miche's limit was calculated, and, the profiles, velocities, and the accelerations were compared with those of 5^th-order Stokes theory. The results give greater velocities and acceleration than 5^th-order Stokes theory, and the wavelength depends on the wave height. The results satisfy the Laplace equation, bottom boundary condition (BBC), and KFSBC, while Stokes theory satisfies only the Laplace equation and BBC. RMSE in DFSBC less than 7.25×10^-2% was obtained. The series order of the proposed method is three, but the series order of 5^th-order Stokes theory is five. Nevertheless, this study provides less RMSE than 5^th-order Stokes theory. As a result, the method is suitable for offshore structural design.

• 유체기계공업학회:학술대회논문집
• /
• 2002.12a
• /
• pp.209-214
• /
• 2002

The present study is to compare the performance of turbulence models in the analysis of the complex flowfield of an axial flow compressor. Baldwin-Lomax turbulence model and k-$\omega$ turbulence model were selected for the comparison. The thin-layer Wavier-Stokes equation was calculated by explicit, finite-difference numerical scheme. A spatially-varying time-step and an implicit residual smoothing were used to improve convergence. Experimental measurements for NASA rotor 37 were cited fer the comparison with numerical data. The compared two turbulence models gave similar performance over all except for total pressure.

• 한국전산유체공학회:학술대회논문집
• /
• 2000.05a
• /
• pp.91-98
• /
• 2000

The three-dimensional unsteady compressible Full Navier-Stokes equation solver with sliding multi-block method has been applied to analyze three dimensional characteristics of the flow field and compression wave around the high speed train which Is entering into a tunnel. The numerical scheme of AF + ADI was used to efficiently solve Navier-Stokes equations in the curvilinear coordinate system. The vortex formation around the nose region was found and the generation of compression wave due to the blockage effects was observed ahead of the train in the form of plane wave. The three dimensional characteristics of the flow field compared to the analytic results were discussed in detail. The variation of pressure of tunnel wall surface and velocity profile of the train are identified as the train enters into a tunnel. The changes in aerodynamic forces and streamlines of each specific sections are also discussed and presented.