• Journal of the Society of Naval Architects of Korea
• /
• v.45 no.1
• /
• pp.29-41
• /
• 2008

The present paper provides a CFD analysis of diffraction problem for a ship with forward speed using an unsteady RANS simulation method, a WAVIS code. The WAVIS viscous solver adopting a finite volume method has second order accuracy in time and field discretizaions for the RANS equations. A two phase level-set method and a realizable ${\kappa}-{\varepsilon}$ turbulence model are adopted to compute the free surface and to meet the turbulence closure, respectively. To validate the capability of the present numerical methods for the simulation of an unsteady progressive regular wave, computations are performed for three grid sets with refinement ratio of ${\sqrt{2}}$. The main simulation is performed for a DTMB5512 model with a forward speed in a regular head sea condition. Validation of the present numerical method is carried out by comparing the present CFD results with available unsteady experimental data published in the 2005 Tokyo CFD Workshop: resistance, heave force, pitch moment, unsteady free surface elevations and velocity fields.

• Journal of Ocean Engineering and Technology
• /
• v.11 no.2
• /
• pp.100-106
• /
• 1997

Viscous flow around an actual ship is calculated by an use of RANS(Reynolds-averaged Navier-Stokes) solver. Reynolds stress is modelled by using k-$\varepsilon$ turbulence model and the law of wall is applied near the body. Body fitted coordinates are introduced for the treatment of the complex boundary of the ship hull form. The transformed equations in the computational domain are numerically solved by an employment of FVM(Finite Volume Method). SIMPLE(Semi-Implcit Pressure Linked Equation) method is adopted in the calculation of pressure and the solution of the disssssssscretized equation is obtained by the line-by-line method with the use of TDMA(Tri-Diagonal Matrix Algorithme). The subject ship model of actual calculation is 4,410 TEU class container carrier. For 4 geosim models the calculated viscous resistancce values are compared with the model test results and analyzed on their componentss. The resistance performance of an actual ship is predicted very resonably, so this mothod may be utilized as a design tool of hull form.

• Journal of computational fluids engineering
• /
• v.11 no.4 s.35
• /
• pp.32-38
• /
• 2006

The objective of this study is to evaluate the prediction accuracy of development large eddy simulation(LES) program for turbulent flow behind a bluff-body. The LES solver was implemented on parallel computer consisting 16 processors. To verify the capability of LES code, the results were compared with those of Reynolds Averaged Navier-Stokes(RANS) using standard ${\kappa}-{\varepsilon}$ model as well as experimental data. The results showed that the LES and RANS qualitatively well predicted the experimental results, such as mean axial, radial velocities and turbulent kinetic energy. In the quantitative analysis, however, the LES showed a better prediction performance than RANS. Specially, the LES well described characteristics of the recirculation zones, such as air stagnation point and jet stagnation point. Finally, the unsteady phenomena on the Bluff-body, such as the transition of recirculation region and vorticity, was examined with LES methodology.

• Journal of Ocean Engineering and Technology
• /
• v.13 no.2 s.32
• /
• pp.138-146
• /
• 1999

The viscous flow around a ship hull is calculated by the use of RANS(Reynolds-averaged Navier-Stokes) solver. Reynolds stresses are midelled by using the k-${epsilon}$ turbulence model and the law is applied near the body. Body fitted corrdinates are introduced for the treatment of the complex boundary of the ship hull form and the governing equations in the physical domain transformed into ones in the computational domain. The transformed equations are numerically solved by an employment of FVM(Finite Volume Method). SIMPLE(Semi-Implicit Pressure Linked Equation) method is adopted in the calculation of pressure and the solution of the sidcretized equation is obtained by the line-by-line method with the use of TDMA(Tri-Diagonal Matrix Algorithme). To assure the proprietty of this computing method, HSVA tanker and Dyne hull are calculated ar both model and ship scale Reynolds number. Their reaults of pressure distributions on fore and aft body, axial velocity contours and transverse velocity velocity vectors and viscous resistance coefficients are compared with other's experiments and calculations.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.7 no.3
• /
• pp.595-609
• /
• 2015

The wake characteristics of Contra-Rotating Propeller (CRP) were investigated using numerical simulation and flow measurement. The numerical simulation was carried out with a commercial CFD code based on a Reynolds Averaged Navier-Stokes (RANS) equations solver, and the flow measurement was performed with Stereoscopic Particle Image Velocimetry (SPIV) system. The simulation results were validated through the comparison with the experiment results measured around the leading edge of rudder to investigate the effect of propeller operation under the conditions without propeller, with forward propeller alone, and with both forward and aft propellers. The evolution of CRP wake was analyzed through velocity and vorticity contours on three transverse planes and one longitudinal plane based on CFD results. The trajectories of propeller tip vortex core in the cases with and without aft propeller were also compared, and larger wake contraction with CRP was confirmed.

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

Deforming mesh should be used when bodies are deforming or moving relative to each other due to the presence of aerodynamic forces and moments. Also, the flow solver for such a flow problem should satisfy the geometric conservation law to ensure the accuracy of the solutions. In this paper, a RANS(Reynolds Averaged Navier-Stokes) solver including automatic mesh capability using TFI(Transfinite Interpolation) method and GCL is developed and applied to flows induced by oscillating wings with given frequencies. The computations are performed both on deforming meshes and on rigid meshes. The computational results are compared with experimental data, which shows a good agreement.

• Journal of computational fluids engineering
• /
• v.15 no.1
• /
• pp.37-45
• /
• 2010

Deforming mesh should be used when bodies are deforming or moving relative to each other due to the presence of aerodynamic forces and moments. Also, the flow solver for such a flow problem should satisfy the geometric conservation law to ensure the accuracy of the solutions. In this paper, a RANS(Reynolds Averaged Navier-Stokes) solver including automatic mesh capability using TFI(Transfinite Interpolation) method and GCL is developed and applied to flows induced by oscillating wings with given frequencies. The computations are performed both on deforming meshes and on rigid meshes. The computational results are compared with experimental data, which shows a good agreement.

• Journal of the Korea Institute of Military Science and Technology
• /
• v.26 no.2
• /
• pp.159-170
• /
• 2023

In the present study, the dynamic behavior characteristics of an amphibious assault vehicle during water entry were analyzed using STAR-CCM+, a commercial computational fluid dynamics(CFD) code. All computations were performed using an overset mesh system and a RANS based flow-solver coupled with dynamic fluid-body interaction(DFBI) solver for simulating three degrees of freedom motion. For numerical validation of the solver, a water entry simulation of inclined circular cylinder was conducted and it was compared between an existing experiment data and CFD results. The pitch angle variation and the trajectory of the circular cylinder during water entry shows good agreement with previous experimental and numerical studies. For the water entry simulations of the amphibious assault vehicle, the analysis of dynamic behaviors of the amphibious assault vehicle with different slope angles, submerged depths and initial velocities were conducted. It is confirmed that the steep slope angle increases the submerged volume of the amphibious assault vehicle, so the buoyancy acting on the vehicle is increased and the moved distance for the re-flotation is decreased. It is also revealed that the submerged volume is increased, bow-up phenomenon occur earlier.

• Nuclear Engineering and Technology
• /
• v.56 no.1
• /
• pp.9-18
• /
• 2024

In-vessel retention through external reactor vessel cooling (IVR-ERVC) is a severe accident management (SAM) strategy that has been adopted and used in many nuclear reactors such as AP1000, APR1400, and light water reactor etc. Some reactor accidents have raised concerns about nuclear reactors among residents, leading to a decrease in residents' acceptability and many studies on SAM are being conducted. Experiments on IVR-ERVC are almost impossible due to its specificity, so fluid characteristics are analyzed through BALI experiments with similar condition. In this study, computational fluid dynamics (CFD) via Reynolds-averaged Navier-Stokes (RANS) and large eddy simulation (LES) for BALI experiments were performed. Steady-state CFD analysis was performed on three turbulence models, and SST k-ω model was in good agreement with the experimental measurement temperature within the maximum error range of 1.9%. LES CFD analysis was performed based on the RANS analysis results and it was confirmed that the temperature and wall heat flux for depth was consistent within an error range of 1.0% with BALI experiment. The LES CFD analysis results were compared with those of the Lagrangian-based solver. LES matched the temperature distribution better than SOPHIA, but SOPHIA calculated the position of boundary between stratified layer and convective layer more accurately. On the other hand, Lagrangian-based solver predicted several small eddy behaviors of the convective layer and LES predicted large vortex behavior. The vibration characteristics near the cooling part of the BALI experimental device were confirmed through Fast Fourier Transform (FFT) investigation. It was found that the power spectral density for pressure at least 10 times higher near the side cooling than near the top cooling.

• Journal of the Society of Naval Architects of Korea
• /
• v.44 no.1 s.151
• /
• pp.16-25
• /
• 2007

Studies of unsteady-airfoil flows have been motivated mostly by efforts to avoid. or reduce such undesirable effects as flutter, noise and vibrations, dynamic stall. In this paper, we carry out a computational study of viscous flows around a two-dimensional oscillating airfoil to investigate unsteady effects in these important and challenging flows. A fully implicit incompressible RANS solver has been used for calculating unsteady viscous flows around an airfoil. The cell-centered End order finite volume method is utilized to discretize governing equations. in order to ease the flow computation for fluid region changing in time, improve the qualify of solution and simplify the grid generation for an oscillating airfoil flow, the computational method adopts a moving and deforming grid generation technique based on the multi-block grid topology. The numerical method is applied for calculating viscous flows of an oscillating NACA 0012 in uniform flow. The computational results are compared with available experimental data. Computed results are compared with experimental data and flow characteristics of the experiment are reproduced well In the computed results.