• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.23 no.1
• /
• pp.104-113
• /
• 1999

A hybrid spatial differencing scheme for the discrete ordinates method is proposed to predict radiative heat transfer in two-dimensional rectangular enclosures. Since this scheme takes the advantages of the diamond scheme and step scheme and includes the characteristics of medium, more accurate and stable results can be obtained. In its development several spatial differencing schemes are examined to address the effect of numerical smearing (or false scattering). Predictions from the proposed hybrid scheme are compared to those of other schemes for transparent, purely absorbing, purely scattering, or absorbing-emitting-isotropically scattering media. It is found that the proposed scheme predicts stable and less smeared results than others.

• Journal of Advanced Research in Ocean Engineering
• /
• v.2 no.2
• /
• pp.48-60
• /
• 2016

Accurate simulation of free-surface wave flows around a ship is very important for better hull-form design. In this paper, a computational fluid dynamics (CFD) code which is based on the open source libraries, OpenFOAM, was developed to predict the wave patterns around a ship. Additional anti-diffusion source term for minimizing a numerical diffusion, which was caused by convection differencing scheme, was considered in the volume-fraction transport equation. The influence of the anti-diffusion source term was tested by applying it to free-surface wave flow around the Wigley and KCS model ships. In results, the wave patterns and hull wave profiles of the Wigley and KCS model ships for various anti-diffusion coefficients showed quite close patterns. While, the band width of the water volume-fraction values between 0.1 to 0.9 at the Wigley and KCS model hull surfaces was narrowed by considering the anti-diffusion term. From the results, anti-diffusion source term decreased free-surface smearing.

• 한국연소학회:학술대회논문집
• /
• 2000.05a
• /
• pp.189-196
• /
• 2000

A numerical simulation was conducted to analyze the interaction of flame and vortices. The characteristic scales of flame and vortices were limited in the thin laminar flamelet regime. Within this regime, flame is assumed as discontinuity surface and its motion in flow field was described by G-equation instead of full governing equations. Additional approximations include distribution of line volume sources on flame surface to simulate effect of volume expansion. Contrast to previous calculations, current study employed vortex transport equation to evaluate attenuation and smearing of vortices. Two extreme conditions of frozen vortex and frozen flame were considered to validate the current method. Comparison with direct numerical simulation resulted in satisfactory quantitative agreement with higher computational efficiency which warrants the usefulness of the present model in more complex situation.

• Journal of the Society of Naval Architects of Korea
• /
• v.41 no.4
• /
• pp.17-29
• /
• 2004

Internal waves are computed using a ghost fluid method on an unstructured grid. Discontinuities in density and dynamic pressure are captured in one cell without smearing or oscillations along a multimaterial interface. A time-accurate incompressible Navier-Stokes/Euler solver is developed based on a three-point backward difference formula for the physical time marching. Artificial compressibility is introduced with respect to pseudotime and an implicit method is used for the pseudotime iteration. To track evolution of an interface, a level set function is coupled with the governing equations. Roe's flux difference splitting method is used to calculate numerical fluxes of the coupled equations. To get higher order accuracy, dependent variables are reconstructed based on gradients which are calculated using Gauss theorem. For each edge crossing an interface, dynamic pressure is assigned for a ghost node to enforce the continuity of total pressure along the interface. Solitary internal waves are computed and the results are compared with other computational and experimental results.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.12
• /
• pp.1555-1569
• /
• 2000

Numerical simulation of fluid flow with moving free surface has been carried out. For the free surface flow, a VOF(Volume of Fluid)-based algorithm utilizing a fixed grid system has been investigated. In order to reduce numerical smearing at the free surface represented on a fixed grid system, a new free surface tracking algorithm based on the donor-acceptor scheme has been presented. Novel features of the proposed algorithm are characterized as two numerical tools; the orientation vector to represent the free surface orientation in each cell and the baby-cell to determine the fluid volume flux at each cell boundary. The proposed algorithm can be easily implemented in any irregular non-uniform grid systems that are usual in finite element method (FEM). Moreover, the proposed algorithm can be extended and applied to the 3-D free surface flow problem without additional efforts. For computation of unsteady incompressible flow, a finite element approximation based on the explicit fractional step method has been adopted. In addition, the SUPG(streamline upwind/Petrov-Galerkin) method has been implemented to deal with convection dominated flows. Combination of the proposed free surface tracking scheme and explicit fractional step formulation resulted in an efficient solution algorithm. Validity of the present solution algorithm was demonstrated from its application to the broken dam and the solitary wave propagation problems.