• Title/Summary/Keyword: Navier solution

Search Result 408, Processing Time 0.024 seconds

Parallel Preconditioner for the Domain Decomposition Method of the Discretized Navier-Stokes Equation (이산화된 Navier-Stokes 방정식의 영역분할법을 위한 병렬 예조건화)

  • Choi, Hyoung-Gwon;Yoo, Jung-Yul;Kang, Sung-Woo
    • Transactions of the Korean Society of Mechanical Engineers B
    • /
    • v.27 no.6
    • /
    • pp.753-765
    • /
    • 2003
  • A finite element code for the numerical solution of the Navier-Stokes equation is parallelized by vertex-oriented domain decomposition. To accelerate the convergence of iterative solvers like conjugate gradient method, parallel block ILU, iterative block ILU, and distributed ILU methods are tested as parallel preconditioners. The effectiveness of the algorithms has been investigated when P1P1 finite element discretization is used for the parallel solution of the Navier-Stokes equation. Two-dimensional and three-dimensional Laplace equations are calculated to estimate the speedup of the preconditioners. Calculation domain is partitioned by one- and multi-dimensional partitioning methods in structured grid and by METIS library in unstructured grid. For the domain-decomposed parallel computation of the Navier-Stokes equation, we have solved three-dimensional lid-driven cavity and natural convection problems in a cube as benchmark problems using a parallelized fractional 4-step finite element method. The speedup for each parallel preconditioning method is to be compared using upto 64 processors.

Analysis of Laminar Flow and Heat Transfer in Asymmetric, Sudden Expansion Channel (비대칭급확대채널의 층류유동 및 열전달 해석)

  • Won, Seung-Ho;Maeng, Joo-Sung;Son, Byung-Jin
    • The Magazine of the Society of Air-Conditioning and Refrigerating Engineers of Korea
    • /
    • v.13 no.1
    • /
    • pp.5-13
    • /
    • 1984
  • This analysis of numerical procedure is prediction of laminar flow and heat transfer at two dimension and steady flow in asymmetric sudden expansion channel. At former study, to analyse the flows with separation, the full Navier-Stokes equation is used, but there are many difficulties to analyse, and although significant progress has been made in the development of efficient computational methods for the Navier-Stokes equations, very large computation times are still required. In case of reward-facing flow, boundary-layer equation is used instead of full Navier-Stokes equation to analyse velocity fields, and result of this numerical analysis is good agreement with the given experimental study. In this case, since the computer time required for the boundary-layer calculation is an order of magnitude less than required for the solution of the full Navier-Stokes equation, this boundary-layer model provides a good approximate solution.

  • PDF

ANALYSIS AND COMPUTATIONS OF OPTIMAL AND FEEDBACK CONTROL PROBLEMS FOR NAVIER-STOKES EQUATIONS

  • Lee, Hyung-Chun
    • Journal of the Korean Mathematical Society
    • /
    • v.34 no.4
    • /
    • pp.841-857
    • /
    • 1997
  • We present analysis and some computational methods for boundary optimal and feedback control problems for Navier-Stokes equations. We use one example to illustrate our methodology and ideas which are applicable to general control problems for Navier-Stokes equations. First, we discuss the existence of optimal solutions and derive an optimality system of equations from which an optimal solution may be computed. Then we present a gradient type iterative method. Finally, we present some numerical results.

  • PDF

AN IMPLICIT NUMERICAL SCHEME FOR SOLUTION OF INCOMPRESSIBLE NAVIER-STOKES EQUATIONS ON CURVILINEAR GRIDS

  • Fayyaz, Hassan;Shah, Abdullah
    • Bulletin of the Korean Mathematical Society
    • /
    • v.55 no.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.

A BOUNDARY CONTROL PROBLEM FOR THE TIME-DEPENDENT 2D NAVIER-STOKES EQUATIONS

  • Kim, Hongchul;Kim, Seon-Gyu
    • Korean Journal of Mathematics
    • /
    • v.16 no.1
    • /
    • pp.57-84
    • /
    • 2008
  • In this paper, a boundary control problem for a flow governed by the time-dependent two dimensional Navier-Stokes equations is considered. We derive a mathematical formulation and a relevant process for an appropriate control along the part of the boundary to minimize the drag due to the flow. After showing the existence of an optimal solution, the first order optimality conditions are derived. The strict differentiability of the state solution in regard to the control parameter shall be exposed rigorously, and the necessary conditions along with the system for the optimal solution shall be deduced in conjunction with the evaluation of the first order Gateaux derivative to the performance functional.

  • PDF

ON THE LINEARIZATION OF DEFECT-CORRECTION METHOD FOR THE STEADY NAVIER-STOKES EQUATIONS

  • Shang, Yueqiang;Kim, Do Wan;Jo, Tae-Chang
    • Journal of the Korean Mathematical Society
    • /
    • v.50 no.5
    • /
    • pp.1129-1163
    • /
    • 2013
  • Based on finite element discretization, two linearization approaches to the defect-correction method for the steady incompressible Navier-Stokes equations are discussed and investigated. By applying $m$ times of Newton and Picard iterations to solve an artificial viscosity stabilized nonlinear Navier-Stokes problem, respectively, and then correcting the solution by solving a linear problem, two linearized defect-correction algorithms are proposed and analyzed. Error estimates with respect to the mesh size $h$, the kinematic viscosity ${\nu}$, the stability factor ${\alpha}$ and the number of nonlinear iterations $m$ for the discrete solution are derived for the linearized one-step defect-correction algorithms. Efficient stopping criteria for the nonlinear iterations are derived. The influence of the linearizations on the accuracy of the approximate solutions are also investigated. Finally, numerical experiments on a problem with known analytical solution, the lid-driven cavity flow, and the flow over a backward-facing step are performed to verify the theoretical results and demonstrate the effectiveness of the proposed defect-correction algorithms.

Development of An Unsteady Navier-Stokes Solver using Implicit Dual Time Stepping Method and DADI Scheme (내재적 이중시간 전진기법과 DADI 기법을 이용한 비정상 Navier-Stokes 코드개발)

  • Lee, Eun-Seok
    • Journal of the Korean Society for Aeronautical & Space Sciences
    • /
    • v.33 no.9
    • /
    • pp.34-40
    • /
    • 2005
  • In present study, a two dimensional unsteady Navier-Stokes solver has been developed using the Diagonalized ADI (DADI) method and implicit dual time stepping method. The jacobian matrices in steady state Navier-Stokes equations are introduced from inviscid flux terms. The implicit treatment of artificial dissipation terms results in a block penta-diagonal matrix system and it becomes a scalar penta-diagonal matrix by diagonalization. In steady state equations about fictitious time, a new residual including a real time derivative term is introduced. From a converged solution about fictitious time, a real time unsteady solution can be obtained, which is called 'implicit dual time stepping method'. For code validation, an oscillating flat plate, a regular Karman vortices past a circular cylinder and shock buffeting around a bicircular airfoil problems are numerically solved. And they are compared with a theoretical solution, experiments and other researcher's computations.