Numerical solutions of 2D Navier-Stokes equations based on generalized harmonic polynomial cell method with non-uniform grid

2022 ◽  
pp. 1-15
Author(s):  
Xueying Yu ◽  
Yanlin Shao ◽  
David R. Fuhrman
Keyword(s):  
Poisson Equation ◽  
Numerical Solutions ◽  
Stokes Equations ◽  
Harmonic Polynomial ◽  
Second Order ◽  
Navier Stokes ◽  
Immersed Boundary ◽  
Navier Stokes Equations ◽  
Uniform Grids ◽  
Global Accuracy

Abstract It is essential for a Navier-Stokes equations solver based on a projection method to be able to solve the resulting Poisson equation accurately and efficiently. In this paper, we present numerical solutions of the 2D Navier-Stokes equations using the fourth-order generalized harmonic polynomial cell (GHPC) method as the Poisson equation solver. Particular focus is on the local and global accuracy of the GHPC method on non-uniform grids. Our study reveals that the GHPC method enables use of more stretched grids than the original HPC method. Compared with a second-order central finite difference method (FDM), global accuracy analysis also demonstrates the advantage of applying the GHPC method on stretched non-uniform grids. An immersed boundary method is used to deal with general geometries involving the fluid-structure-interaction problems. The Taylor-Green vortex and flow around a smooth circular cylinder and square are studied for the purpose of verification and validation. Good agreement with reference results in the literature confirms the accuracy and efficiency of the new 2D Navier-Stokes equation solver based on the present immersed-boundary GHPC method utilizing non-uniform grids. The present Navier-Stokes equations solver uses second-order FDM for the discretization of the diffusion and advection terms, which may be replaced by other higher-order schemes to further improve the accuracy.

A 2D Navier-Stokes Equations Solver Based on Generalized Harmonic Polynomial Cell Method With Non-Uniform Grid

2021 ◽  
Author(s):  
Xueying Yu ◽  
Yanlin Shao ◽  
David R. Fuhrman
Keyword(s):  
Poisson Equation ◽  
Stokes Equation ◽  
Laplace Equation ◽  
Stokes Equations ◽  
Harmonic Polynomial ◽  
Second Order ◽  
Navier Stokes ◽  
Navier Stokes Equation ◽  
Navier Stokes Equations ◽  
Uniform Grids

Abstract It is essential for a Navier-Stokes equations solver based on a projection method to be able to solve the resulting Poisson equation accurately and efficiently. In this paper, we present a new 2D Navier-Stokes equation solver based on a recently proposed fourth-order method, namely the generalized harmonic polynomial cell (GHPC) method, as the Poisson equation solver. The GHPC method is a generalization of the 2D HPC method originally developed for the Laplace equation. In the recent development of the HPC method, loss of accuracy on highly stretched or distorted grids has been reported when solving the Laplace equation, while the performance of the GHPC method on non-uniform grids is still not explored and discussed in the literature. Therefore, the local accuracy of the GHPC method is investigated in detail in the present study, which reveals that the GHPC method allows for the use of much larger grid aspect ratio than that for the original HPC method. Global accuracy of the GHPC method on stretched non-uniform girds is also thoroughly analyzed by considering cases with analytical solutions. Obvious advantages of using the GHPC method in terms of accuracy are demonstrated by comparing with a second-order central Finite Difference Method (FDM). The present Navier-Stokes equations solver uses second-order FDMs for the discretization of the diffusion and advection terms, which may be replaced by other higher-order schemes to further improve the accuracy. Meanwhile, an immersed boundary method [1] is used to study the fluid-structure-interaction problems. The Taylor-Green vortex and flow around a smooth circular cylinder are studied to confirm the accuracy and efficiency of the new 2D Navier-Stokes equation solver. The predictions show good agreements with the experimental and numerical results in the literature.

Steady flow through a channel with a symmetrical constriction in the form of a step

1980 ◽  
Vol 372 (1750) ◽  
pp. 393-414 ◽  
Keyword(s):  
Asymptotic Theory ◽  
Numerical Solutions ◽  
Stokes Equations ◽  
Volumetric Flow Rate ◽  
Reynolds Numbers ◽  
Navier Stokes ◽  
Navier Stokes Equations ◽  
Uniform Grids ◽  
Flow Through ◽  
Good Agreement

Numerical solutions of the Navier-Stokes equations are given for the steady, two-dimensional, laminar flow of an incompressible fluid through a channel with a symmetric constriction in the form of a semi-infinite step change in width. The flow proceeds from a steady Poiseuille velocity distribution far enough upstream of the step in the wider part of the channel to a corresponding distribution downstream in the narrower part and is assumed to remain symmetrical about the centre line of the channel. The numerical scheme involves an accurate and efficient centred difference treatment developed by Dennis & Hudson (1978) and results are obtained for Reynolds numbers, based on half the volumetric flow rate, up to 2000. For a step that halves the width of the channel it is found that very fine uniform grids, with 80 intervals spaced across half of the wider channel upstream, are necessary for resolution of the solution for the flow downstream of the onset of the step. Slightly less refined grids are adequate to resolve the flow upstream. The calculated flow ahead of the step exhibits very good agreement with the asymptotic theory of Smith (1979 b)for Reynolds numbers greater than about 100; indeed, comparisons of the upstream separation position and of the wall vorticity nearby are believed to yield the best agreement between numerical and asymptotic solutions yet found. Downstream there is also qualitative agreement regarding separation and reattachment as the grid size is refined in the numerical results.

Pontryagin maximum principle and second order optimality conditions for optimal control problems governed by 2D nonlocal Cahn–Hilliard–Navier–Stokes equations

Analysis ◽  
2020 ◽  
Vol 40 (3) ◽  
pp. 127-150
Author(s):  
Tania Biswas ◽  
Sheetal Dharmatti ◽  
Manil T. Mohan
Keyword(s):  
Optimal Control ◽  
Optimal Control Problem ◽  
Control Problem ◽  
Stokes Equations ◽  
Minimum Principle ◽  
Immiscible Fluids ◽  
Second Order ◽  
Navier Stokes ◽  
Distributed Optimal Control ◽  
Navier Stokes Equations

AbstractIn this paper, we formulate a distributed optimal control problem related to the evolution of two isothermal, incompressible, immiscible fluids in a two-dimensional bounded domain. The distributed optimal control problem is framed as the minimization of a suitable cost functional subject to the controlled nonlocal Cahn–Hilliard–Navier–Stokes equations. We describe the first order necessary conditions of optimality via the Pontryagin minimum principle and prove second order necessary and sufficient conditions of optimality for the problem.

A numerical study of the interaction between unsteay free-stream disturbances and localized variations in surface geometry

1989 ◽  
Vol 209 ◽  
pp. 285-308 ◽  
Author(s):  
R. J. Bodonyi ◽  
W. J. C. Welch ◽  
P. W. Duck ◽  
M. Tadjfar
Keyword(s):  
Numerical Solutions ◽  
Free Stream ◽  
Stokes Equations ◽  
Numerical Study ◽  
Navier Stokes ◽  
Surface Geometry ◽  
Navier Stokes Equations ◽  
S Waves ◽  
Tollmien Schlichting

A numerical study of the generation of Tollmien-Schlichting (T–S) waves due to the interaction between a small free-stream disturbance and a small localized variation of the surface geometry has been carried out using both finite–difference and spectral methods. The nonlinear steady flow is of the viscous–inviscid interactive type while the unsteady disturbed flow is assumed to be governed by the Navier–Stokes equations linearized about this flow. Numerical solutions illustrate the growth or decay of the T–S waves generated by the interaction between the free-stream disturbance and the surface distortion, depending on the value of the scaled Strouhal number. An important result of this receptivity problem is the numerical determination of the amplitude of the T–S waves.

NUMERICAL SOLUTIONS OF 2D AND 3D SLAMMING PROBLEMS

2021 ◽  
Vol 153 (A2) ◽  
Author(s):  
Q Yang ◽  
W Qiu
Keyword(s):  
Free Surface ◽  
Numerical Solutions ◽  
Stokes Equations ◽  
Type Equation ◽  
The Body ◽  
Navier Stokes ◽  
Time Step ◽  
Navier Stokes Equations ◽  
Highly Nonlinear ◽  
2D And 3D

Slamming forces on 2D and 3D bodies have been computed based on a CIP method. The highly nonlinear water entry problem governed by the Navier-Stokes equations was solved by a CIP based finite difference method on a fixed Cartesian grid. In the computation, a compact upwind scheme was employed for the advection calculations and a pressure-based algorithm was applied to treat the multiple phases. The free surface and the body boundaries were captured using density functions. For the pressure calculation, a Poisson-type equation was solved at each time step by the conjugate gradient iterative method. Validation studies were carried out for 2D wedges with various deadrise angles ranging from 0 to 60 degrees at constant vertical velocity. In the cases of wedges with small deadrise angles, the compressibility of air between the bottom of the wedge and the free surface was modelled. Studies were also extended to 3D bodies, such as a sphere, a cylinder and a catamaran, entering calm water. Computed pressures, free surface elevations and hydrodynamic forces were compared with experimental data and the numerical solutions by other methods.

Sign in / Sign up

Export Citation Format

Share Document