• 스마트미디어저널
• /
• 제4권2호
• /
• pp.9-16
• /
• 2015

Fluids appear in innumerable phenomena; therefore, it is interesting to reproduce those phenomena by computer graphics techniques. However, this process is not trivial. We work with a fluid simulation that uses Navier-Stokes equations to model the fluid, a semi-Lagrangian approach to solve it and the level set method to track the surface of the fluid. Modified versions of the Navier-Stokes equations for computer graphics allow us to create a wide diversity of effects. In this paper, we propose a technique that allows us to integrate a force inspired by surface tension into the model. We describe which information we need and how to modify the model with this new approach. We end up with a modified simulation that has additional effects that might be suitable for computer graphics purposes. The effects that we are able to recreate are small waves and droplet-like formations close to the surface of the fluid. This model preserves the overall behavior governed by the Navier-Stokes equations.

• 대한기계학회논문집
• /
• 제17권9호
• /
• pp.2329-2338
• /
• 1993

An accurate analysis procedure to solve the flow about a flat plate at various incidences has been developed. The Navier-Stokes equations of stream function and vorticity form are solved in a sufficiently large computational domain, in which the grid lines are mutually orthogonal. The details of the flow near the singularity at the tip of the plate is well captured by the analytic solution which is asymptotically matched to the numerically generated outer solution. The solution for each region is obtained iteratively : the solution of one (inner or outer) region uses that of the other as the boundary condition after each cycle. The resulting procedure is accurate everywhere and also computationally efficient as the singularity has been removed. It is applied to the flat plate for a wide range of Re : the results agree very well with the existing computation and experiment.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2009년 추계학술대회논문집
• /
• pp.27-33
• /
• 2009

In this paper, the flow past a rotary oscillating circular cylinder is simulated. The high-order and high-resolution numerical schemes with the characteristic boundary conditions are used for the compressible Navier-Stokes equation. The frequencies of rotating oscillation are $0.19\;{\leq}\;S_f\;{\leq}\;0.25$ for the maximum angular $\theta_{max}=10^{\circ}$ and $17^{\circ}$. The flow conditions are Mach number of 0.3 and Reynolds number of 1000. At Lock-on and Non-lock-on region which are defined by the relation between the vortex shedding frequency and the oscillating frequency, the drag and lift coefficient are analyzed.

• 한국전산유체공학회지
• /
• 제4권3호
• /
• pp.12-20
• /
• 1999

The computations of the flowfield and pollutant dispersion over a flat plate and the Russian hills of various slopes are described. The Gaussian plume and the puff model have been used to calculate concentration of pollutant. The Reynolds-averaged unsteady incompressible Navier-Stokes equation with low Reynolds κ-ε model has been used to calculate the flowfield. The flow data of a flat plate and the Russian hills from Navier-Stokes equation solutions has been used as the input data for the puff model. The computational results of flowfield agree well with experimental results of both a flat plate and Russian hills. The concentration prediction by the Gaussian plume model and the Gaussian puff model also agrees flirty well with experiments.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제14권2호
• /
• pp.125-140
• /
• 2010

This paper develops a least-squares approach to the solution of the optimal control problem for the Navier-Stokes equations. We recast the optimality system as a first-order system by introducing velocity-flux variables and associated curl and trace equations. We show that a least-squares principle based on $L^2$ norms applied to this system yields optimal discretization error estimates in the $H^1$ norm in each variable.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제14권1호
• /
• pp.43-56
• /
• 2010

The representative numerical algorithms to solve the time dependent Navier-Stokes equations are projection type methods. Lots of projection schemes have been developed to find more accurate solutions. But most of projection methods [4, 11] suffer from inconsistency and requesting unknown datum. E and Liu in [5] constructed the gauge method which splits the velocity $u=a+{\nabla}{\phi}$ to make consistent and to replace requesting of the unknown values to known datum of non-physical variables a and ${\phi}$. The errors are evaluated in [9]. But gauge method is not still obvious to find out suitable combination of discrete finite element spaces and to compute boundary derivative of the gauge variable ${\phi}$. In this paper, we define 4 gauge algorithms via combining both 2 decomposition operators and 2 boundary conditions. And we derive variational derivative on boundary and analyze numerical results of 4 gauge algorithms in various discrete spaces combinations to search right discrete space relation.

• 대한수학회지
• /
• 제48권6호
• /
• pp.1153-1170
• /
• 2011

In this paper, spatial decay estimates for the time dependent compressible viscous isentropic flow in a semi-infinite three dimensional pipe are derived. An upper bound for the total energy in terms of the initial boundary data is obtained as well. The results established in this paper may be viewed as a version of Saint-Venant's principle in transient compressible Navier-Stokes flow.

• 대한수학회보
• /
• 제55권3호
• /
• pp.881-898
• /
• 2018

This article deals with implementation of a high-order finite difference scheme for numerical solution of the incompressible Navier-Stokes equations on curvilinear grids. The numerical scheme is based on pseudo-compressibility approach. A fifth-order upwind compact scheme is used to approximate the inviscid fluxes while the discretization of metric and viscous terms is accomplished using sixth-order central compact scheme. An implicit Euler method is used for discretization of the pseudo-time derivative to obtain the steady-state solution. The resulting block tridiagonal matrix system is solved by approximate factorization based alternating direction implicit scheme (AF-ADI) which consists of an alternate sweep in each direction for every pseudo-time step. The convergence and efficiency of the method are evaluated by solving some 2D benchmark problems. Finally, computed results are compared with numerical results in the literature and a good agreement is observed.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2001년도 춘계 학술대회논문집
• /
• pp.149-154
• /
• 2001

Parallel performance of a Myrinet based PC-cluster was tested and compared with a conventional Fast-Ethernet system. A preconditioned Navier-Stokes code was parallelized with domain decomposition technique, and used for the parallel performance test. Speed-up ratio was examined as a major performance parameter depending on the number of processor and the network topology. As was expected, Myrinet system shows a superior parallel performance to the Fast-Ethernet system even with a single network adpater for a dual processor SMP machine. A test for the dependency on problem size also shows that network communication speed is a crucial factor for parallelized computational fluid dynamics analysis and the Myrinet system is a plausible candidate for high performance parallel computing system.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2001년도 추계 학술대회논문집
• /
• pp.91-97
• /
• 2001

The thin-layer Navier-Stokes equations are solved for the hypersonic flow over blunt cone configurations with applications to laminar as well as turbulent flows. The equations are expressed in the forms of flux-vector splitting and explicit algorithm. The upwind schemes of Steger-Warming and van Leer are investigated in their ability to accurately predict the heating loads along the surface of the body. A comparison with the second order extensions of these schemes is made and a hybrid scheme incorporating a combination of central differencing and flux-vector-splitting is presented. This scheme is also investigated in its ability to accurately predict heat transfer distributions.