• Journal of the Korea Institute of Military Science and Technology
• /
• v.9 no.2 s.25
• /
• pp.152-160
• /
• 2006

Time marching methods are well-suited for high speed compressible flow computations. However, it is well known that the time marching methods suffer a slow down in convergence due to disparity in Eigenvalues. A local preconditioning method is one of numerical methods to enhance convergence characteristics of low mach number flows by modifying Eigenvalues of the governing equations. In this paper, the local preconditioning method of Weiss is applied to a 2 dimensional Navier-Stokes code and the efficiency of the preconditioning method is shown through a number of computational examples.

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

The local preconditioning method has both robust convergence and accurate solutions by using local flow properties for parameters in the preconditioning matrix. Preconditioning methods have been very effective to low speed inviscid flows. In the viscous and turbulent flows, deterioration of convergence should be overcame on the high aspect ratio grids to get better convergence and accuracy. In the present study, the local time stepping and min-CFL/max-VNN definitions are applied to compare the results and we propose the method that switches between two methods. The min-CFL definition is applied for inviscid flow problems and the min-CFL/max-VNN definition is implemented to viscous and turbulent flow problems.

• 한국전산유체공학회:학술대회논문집
• /
• 2004.10a
• /
• pp.7-14
• /
• 2004

It is well known that preconditioning methods are efficient for convergence acceleration at compressible low Mach number flows. In this study, the original Euler equations and three preconditioners nondimensionalized differently are implemented in two dimensional inviscid bump flows using the 3rd order MUSCL and DADI schemes as flux discretization and time integration respectively. The multigrid and local time stepping methods are also used to accelerate the convergence. The test case indicates that a properly modified local preconditioning technique involving concepts of a global preconditioning one produces Mach number independent convergence. Besides, an asymptotic analysis for properties of preconditioning methods is added.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.33 no.8
• /
• pp.8-17
• /
• 2005

It is well known that preconditioning methods are efficient for convergence acceleration in the compressible low Mach number flows. In this study, the original Euler equations and three differently nondimensionalized preconditioning methods are implemented in two dimensional inviscid bump flows using the 3rd order MUSCL and DADI schemes as numerical flux discretization and time integration, respectively. The multigrid and local time stepping methods are also used to accelerate the convergence. The test case indicates that a properly modified local preconditioning technique involving concepts of a global preconditioning allows Mach number independent convergence. Besides, an asymptotic analysis for properties of preconditioning methods is added.

• Journal of the Korean Mathematical Society
• /
• v.41 no.2
• /
• pp.345-368
• /
• 2004

A preconditioning algorithm is developed in this paper for the iterative solution of the linear system of equations resulting from the p-version boundary element approximation of the three dimensional integral equation with hypersingular operators. The preconditioner is derived by first making the nodal and side basis functions locally orthogonal to the element internal bases, and then by decoupling the nodal and side bases from the internal bases. Its implementation consists of solving a global problem on the wire-basket and a series of local problems defined on a single element. Moreover, the condition number of the preconditioned system is shown to be of order $O((1＋ln/p)^{7})$. This technique can be applied to discretization with triangular elements and with general basis functions.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.03b
• /
• pp.13-16
• /
• 2008

After the collapse of the Tacoma bay bridge at Tacoma Washington, the accurate prediction of aerodynamics became crucial to the sound design of bridges. CFD(Computational Fluid Dynamics) becomes important tool for the prediction on wind effects on the bridge due to the recent development of CFD. The usage of CFD is further prompted by the advantages in using CFD, such as low-cost and fast feed-back of design. In this paper, an unsteady compressible Reynolds averaged Navier-Stokes code is used for the computation of the flow over bridges. Coakley's ��q-${\omega}$ �� two-equation turbulence model is used for the turbulent eddy viscosity. For accurate and stable computations, the local preconditioning method is adapted to the code. Aerodynamic characteristics of a couple bridges are presented to show the validity and the accuracy of the method.

• Communications of the Korean Mathematical Society
• /
• v.16 no.1
• /
• pp.153-161
• /
• 2001

In the course of working on the preconditioning $C^1$ polynomial spline collocation method, one has to deal with the exponential decay of $C^1$ Lagrange polynomial splines. In this paper we show the exponential decay of $C^1$ Lagrange polynomial splines using the Chebyshev-Gauss points as the local data points.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.36 no.10
• /
• pp.939-946
• /
• 2008

Euler codes or Navier-Stokes codes for compressible flows suffer severe degradation in convergence as Mach number approaches zero. The convergence problem arose from the wide disparity in characteristic speeds can be solved using preconditioning methods without large modifications. In this paper, a preconditioned RANS(Reynolds Averaged Navier-Stokes) solver is developed for analysis of low Mach number flows. In order to validate the method, computational examples are chosen and the results are compared with the experimental data and the existing computed results showing a good accuracy and convergence characteristics for steady inviscid, laminar and turbulent flows at low Mach number.

• Nuclear Engineering and Technology
• /
• v.52 no.2
• /
• pp.258-270
• /
• 2020

A new task of using the Jacobian-Free-Newton-Krylov (JFNK) method for the PWR core transient simulations involving neutronics, thermal hydraulics and mechanics is conducted. For the transient scenario of PWR, normally the Picard iteration of the coupled coarse-mesh nodal equations and parallel channel TH equations is performed to get the transient solution. In order to solve the coupled equations faster and more stable, the Newton Krylov (NK) method based on the explicit matrix was studied. However, the NK method is hard to be extended to the cases with more physics phenomenon coupled, thus the JFNK based iteration scheme is developed for the nodal method and parallel-channel TH method. The local gap conductance is sensitive to the gap width and will influence the temperature distribution in the fuel rod significantly. To further consider the local gap conductance during the transient scenario, a 1D mechanics model is coupled into the JFNK scheme to account for the fuel thermal expansion effect. To improve the efficiency, the physics-based precondition and scaling technique are developed for the JFNK iteration. Numerical tests show good convergence behavior of the iterations and demonstrate the influence of the fuel thermal expansion effect during the rod ejection problems.