• Environmental Engineering Research
• /
• v.19 no.2
• /
• pp.165-174
• /
• 2014

A model combining multi-dimensional discretized population balance equations with a computational fluid dynamics simulation (CFD-DPBE model) was developed and applied to simulate turbulent flocculation and sedimentation processes in sediment retention basins. Computation fluid dynamics and the discretized population balance equations were solved to generate steady state flow field data and simulate flocculation and sedimentation processes in a sequential manner. Up-to-date numerical algorithms, such as operator splitting and LeVeque flux-corrected upwind schemes, were applied to cope with the computational demands caused by complexity and nonlinearity of the population balance equations and the instability caused by advection-dominated transport. In a modeling and simulation study with a two-dimensional simplified pond system, applicability of the CFD-DPBE model was demonstrated by tracking mass balances and floc size evolutions and by examining particle/floc size and solid concentration distributions. Thus, the CFD-DPBE model may be used as a valuable simulation tool for natural and engineered flocculation and sedimentation systems as well as for flocculant-aided sediment retention ponds.

• Journal of computational fluids engineering
• /
• v.4 no.1
• /
• pp.53-61
• /
• 1999

Two-dimensional, unsteady, incompressible and compressible Navier-Stokes codes are developed for the computation of the viscous turbulent flow over high-lift airfoils. The compressible code involves a conventional upwind-differenced scheme for the convective terms and LU-SGS scheme for temporal integration. The incompressible code with pseudo-compressibility method also adopts the same schemes as the compressible code. Three two-equation turbulence models are evaluated by computing the flow over single and multi-element airfoils. The compressible and incompressible codes are validated by predicting the flow around the RAE 2822 transonic airfoil and the NACA 4412 airfoil, respectively. In addition, both the incompressible and compressible code are used to compute the flow over the NLR 7301 airfoil with flap to study the compressible effect near the high-loaded leading edge. The grid systems are efficiently generated using Chimera overlapping grid scheme. Overall, the κ-ω SST model shows closer agreement with experiment results, especially in the prediction of adverse pressure gradient region on the suction surfaces of high-lift airfoils.

• Coupled systems mechanics
• /
• v.8 no.2
• /
• pp.129-146
• /
• 2019

The main goal of this work is to develop a numerical simulator to study an isothermal single-phase two-component flow in a naturally fractured oil reservoir, taking into account advection and diffusion effects. We use the Peng-Robinson equation of state with a volume translation to evaluate the properties of the components, and the discretization of the governing partial differential equations is carried out using the Finite Difference Method, along with implicit and first-order upwind schemes. This process leads to a coupled non-linear algebraic system for the unknowns pressure and molar fractions. After a linearization and the use of an operator splitting, the Conjugate Gradient and Bi-conjugated Gradient Stabilized methods are then used to solve two algebraic subsystems, one for the pressure and another for the molar fraction. We studied the effects of fractures in both the flow field and mass transport, as well as in computing time, and the results show that the fractures affect, as expected, the flow creating a thin preferential path for the mass transport.

• Journal of computational fluids engineering
• /
• v.19 no.3
• /
• pp.91-97
• /
• 2014

A high resolution numerical method aimed at solving cavitating flow was proposed and applied to gas-liquid two-phase shock tube problem with arbitrary void fraction. The present method with compressibility effects employs a finite-difference 4th-order Runge-Kutta method and Roe's flux difference splitting approximation with the MUSCL TVD scheme. The Jacobian matrix from the inviscid flux of constitute equation is diagonalized analytically and the speed of sound for the two-phase media is derived by eigenvalues. So that the present method is appropriate for the extension of high order upwind schemes based on the characteristic theory. By this method, a Riemann problem for Euler equations of one dimensional shock tube was computed. Numerical results of high speed flow phenomena such as detailed observations of shock and expansion wave propagations through the gas-liquid two-phase media and some data related to computational efficiency are made. Comparisons of predicted results and solutions at isothermal condition are provided and discussed.

• Journal of the Society of Naval Architects of Korea
• /
• v.39 no.1
• /
• pp.8-15
• /
• 2002

A two-dimensional higher order panel method using B-splines has been developed to overcome the disadvantages of the low order panel method and to obtain more accurate solution. The sources and the normal dipoles are distributed on both the body and the free surface. Instead of applying the upwind finite difference schemes to satisfy the linearized free surface and the radiation condition, the derivatives of the basis functions of the B-splines are directly applied to the linearized free surface condition. Numerical damping in the Dawson's method are avoided in the Present computations. In order to validate the present method, numerical computations are carried out for a submerged cylinder and a two-dimensional hydrofoil steadily moving beneath a free surface. The numerical results show that fast convergence and better accuracies have been achieved by the present method.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.33 no.4
• /
• pp.1-6
• /
• 2005

The ability to predict radar return from flying vehicles becomes a critical technology issue in the development of stealth configurations. Toward developing a CEM code based on Maxwell's equations for analysis of RCS reduction schemes, an explicit upwind scheme suitable for multidisciplinary design is presented. The DFFT algorithm is utilized to convert the time-domain field values to the frequency-domain. A Green's function based on near field-to-far field transformation is also employed to calculate the bistatic RCS. To verify the numerical calculation the two-dimensional field around a perfectly conducting cylinder is considered. Finally results are obtained for the scattering electromagnetic field around an airfoil in order to illustrate the feasibility of applying CFD based methods to CEM.