• The KSFM Journal of Fluid Machinery
• /
• v.18 no.2
• /
• pp.54-59
• /
• 2015

The performance of UV reactor which is used in water treatment is strongly affected by UV fluence rate and water flow in the UV reactor. Therefore, CFD tools are widely used in designing process of UV reactors. This paper describes the development of a computational fluid dynamics (CFD) methodology that can be used to calculate the performance of open channel type UV reactor used in wastewater treatment plant. All computations were performed using commercial CFD code, CFX, by considering three dimensional, steady, incompressible flow. The Eulerian-Eulerian multi-phase method were used to capture the water-air interface. The MSSS model, provided by UVCalc3D, was used to calculate the UV intensity field. The numerical predictions and calculated UV Dose were compared with experimental dataset to validate the CFD methodology. The reactor performance based on MS2 log reduction was well matched with measurements within 6%.

• Wind and Structures
• /
• v.1 no.1
• /
• pp.111-125
• /
• 1998

This paper deals with the development of variable-node element and its application to the adaptive h-version mesh refinement-recovery for the incompressible viscous flow analysis. The element which has variable mid-side nodes can be used in generating the transition zone between the refined and unrefined element and efficiently used for the construction of a refined mesh without generating distorted elements. A modified Guassian quadrature is needed to evaluate the element matrices due to the discontinuity of derivatives of the shape functions used for the element. The penalty function method which can reduce the number of the independent variables is adopted for the purpose of computational efficiency and the selective reduced integration is carried out for the convection and pressure terms to preserve the stability of solution. For the economical analysis of transient problems in which the locations to be refined are changed in accordance with the dynamic distribution of velocity gradient, not only the mesh refinement but also the mesh recovery is needed. The numerical examples show that the optimal mesh for the finite element analysis of a wind around the structures can be obtained automatically by the proposed scheme.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.11 no.1
• /
• pp.530-541
• /
• 2019

To improve our current understanding of tsunami-like solitary waves interacting with a row of vertical slotted piles on a sloping beach, a 3D numerical wave tank based on the CFD tool $OpenFOAM^{(R)}$ was developed in this study. The Navier-Stokes equations were employed to solve the two-phase incompressible flow, combining with an improved VOF method to track the free surface and a LES model to resolve the turbulence. The numerical model was firstly validated by our laboratory measurements of wave, flow and dynamic pressure around both a row of piles and a single pile on a slope subjected to solitary waves. Subsequently, a series of numerical experiments were conducted to analyze the breaking wave force in view of varying incident wave heights, offshore water depths, spaces between adjacent piles and beach slopes. Finally, a slamming coefficient was discussed to account for the breaking wave force impacting on the piles.

• Journal of the Korean Society of Visualization
• /
• v.20 no.2
• /
• pp.45-54
• /
• 2022

We analyzed the low wall-shear stress area in the intracranial aneurysm that occurred at an anterior communicating artery with a special emphasis on vortical structures close to the wall. We reconstructed the aneurysm model from patient CTA data. We assumed blood as an incompressible Newtonian fluid and treated the blood vessel as a solid wall. The pulsatile boundary condition was applied at the inlet of the anterior cerebral artery. From the instantaneous flow field, we computed the histogram of the wall-shear stress over the aneurysm wall and found the low wall-shear stress event (< 0.4 Pa). This extreme event was due to the low wall-shear stress area that occurred at the daughter sac. We found that the merging of two vortices induced the low wall-shear stress area; one arises from the morphological characteristics of the daughter sac, and the other is formed by a jet flow into the aneurysm sac. The latter approaches the daughter sac, which ultimately leads to the strong ejection event near the daughter sac.

• Proceedings of the KSME Conference
• /
• 2007.05b
• /
• pp.2951-2954
• /
• 2007

Tumor hemodynamics in vascular state is numerically simulated using pressure node solution. The tumor angiogenesis pattern in our previous study is used for the geometry of vessel networks. For tumor angiogenesis, the equation that governed angiogenesis comprises a tumor angiogenesis factor (TAF) conservation equation in time and space, which is solved numerically using the Galerkin finite element method. A stochastic process model is used to simulate vessel formation and vessel. In this study, we use a two-dimensional model with planar vessel structure. Hemodynamics in vessel is assumed as incompressible steady flow with Newtonian fluid properties. In parent vessel, arterial pressure is assigned as a boundary condition whereas a constant terminal pressure is specified in tumor inside. Kirchhoff's law is applied to each pressure node to simulate the pressure distribution in vessel networks. Transient pressure distribution along with angiogenesis pattern is presented to investigate the effect of tumor growth in tumor hemodynamics.

• Journal of computational fluids engineering
• /
• v.12 no.1
• /
• pp.35-42
• /
• 2007

This paper describes a recent development on the divergence free basis function based on a hermite stream function and verifies its validity by comparing results with those from a modified residual method known as one of stabilized finite element methods. It can be shown that a proper choice of degrees of freedom at a node with a proper arrangement of the hermite interpolation functions can yield solenoidal or divergent free interpolation functions for the velocities. The well-known cavity problem has been chosen for validity of the present algorithm. The comparisons from numerical results between the present and the modified residual showed the present method yields better results in both the velocity and the pressure within modest Reynolds numbers(Re = 1,000).

• Journal of the Korea Computer Graphics Society
• /
• v.14 no.4
• /
• pp.1-5
• /
• 2008

In this paper, we propose an efficient technique for improving the grid based fluid simulation by sub-grid visuals. The detailed turbulency generated efficiently by Vortex Particle Method are blended with the flow fields coming from the traditional incompressible Navier-Stokes solver. The algorithm enables large- and small- scale detail to be edited separately.

• Journal of computational fluids engineering
• /
• v.11 no.3 s.34
• /
• pp.59-63
• /
• 2006

Accurate prediction of a supersonic missile base drag continues to defy even well-rounded CFD codes. In an effort to address the accuracy and predictability of the base drags, the influence of grid system and competitive turbulence models on the base drag is analyzed. Characteristics of some turbulence models is reviewed through incompressible turbulent flow over a flat plate, and performance for the base drag prediction of several turbulence models such as Baldwin-Loman(B-L), Spalart-Allmaras(S-A), k-$\varepsilon$, k-$\omega$ model is assessed. When compressibility correction is injected into the S-A model, prediction accuracy of the base drag is enhanced. The NSWC wind tunnel test data are utilized for comparison of CFD and semi-empirical codes on the accuracy of base drag predictability: they are about equal, but CFD tends to perform better. It is also found that, as angle of attack of a missile with control fins increases, even the best CFD analysis tool we have lacks the accuracy needed for the base drag prediction.

• Journal of Advanced Marine Engineering and Technology
• /
• v.27 no.4
• /
• pp.494-502
• /
• 2003

The near-tip field of mode-I dynamic cracks steadily propagating in a strain softening material is investigated under plane strain conditions. The material is assumed to be incompressible and its deformation obeys the $J_2$ flow theory of plasticity. A power-law stress-strain relation with strain softening is adopted to account for the damage behavior of materials near the dynamic crack tip. By assuming that the stresses and strain have the same singularity at the crack tip. this paper obtains a fully continuous dynamic crack-tip field in the damage region. Results show that the stress and strain components the same logarithmic singularity of (In(R/r))$\delta$, and the angular variations of filed quantities are identical to those corresponding to the dynamic cracks in the elastic-perfectly plastic material.

• Journal of Ocean Engineering and Technology
• /
• v.15 no.4
• /
• pp.46-52
• /
• 2001

SMAC method, one of the numerical simulation techniques, is modified from the original MAC method for the time-dependent variation of fluid flows. The Navier-Stokes equations for incompressible time-dependent viscous flow are applied, and Also marker particles which present the visualization of fluid flows are used. In this study, two-dimensional numerical simulations of the overtopping are carried out by SMAC method, and the simulation results are visualized, In addition to, motion pictures are made for efficient visualization of the simulation results. This numerical simulation could also be applied to the design of coastal structures as dike and revetment.