• 한국전산유체공학회:학술대회논문집
• /
• 2009.04a
• /
• pp.355-362
• /
• 2009

Most material of engineering interest undergoes solidification process from liquid to solid state. Identifying the underlying mechanism during solidification process is essential to determine the microstructure of material which governs the physical properties of final product. In this paper, we expand our previous two-dimensional numerical technique to three-dimensional simulation for computing dendritic solidification process with fluid convection. We used Level Contour Reconstruction Method to track the moving liquid-solid interface and Sharp Interface Technique to correctly implement phase changing boundary condition. Three-dimensional results showed clear difference compared to two-dimensional simulation on tip growth rate and velocity.

• Korea-Australia Rheology Journal
• /
• v.13 no.3
• /
• pp.131-139
• /
• 2001

A coupled approach is proposed for the prediction of sheet profile in sheet casting process, which combines one-dimensional analytical method on planar elongational flow region and three-dimensional numerical method on the other region. The strategy is constructed from the observations that the flow domain of sheet casting process can be separated into two parts based old the flow kinematics. The flow field in the central region of sheet, over which the planar elongational flow dominates, is possibly replaced by one-dimensional analytical solution. Then only a partial flow domain near the edge region of sheet, where the flow kinematics cannot be described by the planar elongational flow itself, requires three-dimensional numerical simulation. Good agreement is observed between the coupled approach developed in this study and the full three-dimensional numerical simulation previously developed and reported by the authors. This coupled approach may have provided flexibility with low costs to accommodate a wide range of die sizes in sheet casting process.

• Journal of Environmental Science International
• /
• v.2 no.2
• /
• pp.103-113
• /
• 1993

The land and sea breeze over the Pusan coastal area is studied by three dimensional mesoscale numerical model. According to the results of the simulation experiments, both Pusan areas and Kimhae areas, the sea breeze began at 0800LST and the strongest at 1500LST and then at 1800LST. After midnight, the sea breeze changed about the land breeze and become weaker than that of the sea breeze in the daytime. Comparisons between calculations and observations showed that the characteristics of diurnal variation and v-component of the wind velocity relatively is similar to the Pusan areas. On the Kimhae areas, however, observations showed time lag which compared to the results of simulation experiments in the velocity of sea breeze and diurnal variation. From the above results, comparisons between calculations and observations is much more similar to the coastal areas than on the inland area.

• Proceedings of the KSME Conference
• /
• 2008.11b
• /
• pp.2401-2405
• /
• 2008

In the present study, a three-dimensional least square/level set based two-phase flow code was developed for the simulation of three-dimensional sloshing problems using finite element discretization. The present method can be utilized for the analysis of a free surface flow problem in a complex geometry due to the feature of FEM. Since the finite element method is employed for the spatial discretization of governing equations, an unstructured mesh can be naturally adopted for the level set simulation of a free surface flow without an additional load for the code development except that solution methods of the hyperbolic type redistancing and advection equations of the level set function should be devised in order to give a bounded solution on the unstructured mesh. From the numerical experiments of the present study, it is shown that the proposed method is both robust and accurate for the simulation of three-dimensional sloshing problems.

• Journal of Mechanical Science and Technology
• /
• v.15 no.12
• /
• pp.1750-1756
• /
• 2001

Three-dimensional and unsteady flow analysis is a practical target of high performance computation. As recently advances of computers, a numerical prediction by the large eddy simulation (LES) are introduced and evaluated for various engineering problems. Its advanced methods for the complex turbulent flows are discussed by several examples applied for aerodynamic designs, analysis of fluid flow mechanisms and their interaction to complex phenomena. These results of time-dependent and three-dimensional phenomena are visualized by interactive graphics and animations.

• Journal of Ocean Engineering and Technology
• /
• v.4 no.1
• /
• pp.55-61
• /
• 1990

For the tracking of three dimensional free surface motions, a method referred to as the Volume of Fluid(VOF) algorithm is extended. In order to calculate the slope of three dimensional free surface which is the most important for the advection algorithm that decides the amount of fluid from cell to cell and for the application of free surface boundary condition, a simple method utilizing two dimensional slope informations is introduced. The extended algroithm is tested by demonstrating the simulation of a propagating sinusoidal wave through the channel whose width changes abruptly.

• Proceedings of the Korea Water Resources Association Conference
• /
• 1987.07a
• /
• pp.113-122
• /
• 1987

A numerical simulation technique of three-dimensional finite difference model is developed to study the groundwater flow system in Dcjima, an upland area which faces Kasumigaura Lake. For general perspectives of the groundwater flow system, a steady state three-dimentional model is simulated. For the sedimentary mud formations which are found in the representative formation, three situations of hydraulic conductivity are considered, representing an isotropic condition and situations where the horizontal permeability is equal to 10 times and 100times of the vertical one. The finite difference grid used in the simulation has 60x50x30=90,000 nodes. A converged solution with a tolerance of 0.001 meter of hydraulic head is set. Having determined the flow net by using a steady state three-dimensional model. the results for the three cases of hydraulic conductivity are compared with the results of tracer methods (Bae and Kayane 1987) With the aid of four representative vertical cross-sections, groundwater flow systems in the study area are assumed. Water balances for the three cases indicate very good agreement between total recharge and discharge in each case Analyses of groundwater flow system based on the tritium concentrations and water quality measurements (Bae and Kayane 1987) are confirmed by the numerical simulation and the results obtained by these two methods appeared to be in close agreement.

• Wind and Structures
• /
• v.1 no.2
• /
• pp.127-144
• /
• 1998

In this paper, the wake structures behind a square cylinder at the Reynolds number of 22,000 are simulated using the large eddy simulation, and the main features of the wake structure associated with unsteady vortex-shedding are investigated. The Smagorinsky model is used for parametrization of the subgrid scales. The finite element method with isoparametric linear elements is employed in the computations. Unsteady computations are performed using the explicit method with streamline upwind scheme for the advection term. The time integration incorporates a subcycling strategy. No-slip condition is enforced on the wall surface. A comparative study between two-and three-dimensional computations puts a stress on the three-dimensional effects in turbulent flow simulations. Simulated three-dimensional wake structures are compared with numerical and experimental results reported by other researchers. The results include time-averaged, phase-averaged flow fields and numerically visualized vortex-shedding pattern using streaklines. The results show that dynamics of the vortex-shedding phenomenon are numerically well reproduced using the present method of finite element implementation of large eddy simulation.

• Advances in aircraft and spacecraft science
• /
• v.2 no.1
• /
• pp.77-94
• /
• 2015

The paper focuses on the integration of a non-linear one-dimensional model of Synthetic Jet (SJ) actuator in a well-assessed numerical simulation method for turbulent compressible flows. The computational approach is intended to the implementation of a numerical tool suited for flow control simulations with affordable CPU resources. A strong compromise is sought between the use of boundary conditions or zero-dimensional models and the full simulation of the actuator cavity, in view of long-term simulation with multiple synthetic jet actuators. The model is integrated in a multi-domain numerical procedure where the controlled flow field is simulated by a standard CFD method for compressible RANS equations, while flow inside the actuator is reduced to a one-dimensional duct flow with a moving piston. The non-linear matching between the two systems, which ensures conservation of the mass, momentum and energy is explained. The numerical method is successfully tested against three typical test cases: the jet in quiescent air, the SJ in cross flow and the flow control on the NACA0015 airfoil.

• Proceedings of the KSME Conference
• /
• 2001.11b
• /
• pp.570-577
• /
• 2001

Three-dimensional time-dependent flow past a circular cylinder is numerically investigated using direct numerical simulation for Reynolds number 280 and 300. The higher-order finite difference scheme is employed for the spatial distributions along with the second order Adams-Bashforth and the first order backward-Euler time integration. In x-y plane, the convection term is applied by the 5th order upwind scheme and the pressure and viscosity terms are applied by the 4th order central difference. And in spanwise, Navier-Stokes equation is distributed using of Spectral Method. At Reynolds number 259 the two-dimensional wake becomes linearly unstable to a second branch of modes with wavelength about 1.0 diameters at onset (B-mode). Present results of three-dimensional effects of in wake of a circular cylinder is represented with spanwise and streamwise vorticity contours as Reynolds numbers.