The discrete continuity equation in primitive variable solutions of incompressible flow

doi:10.1016/0021-9991(91)90260-R

Journal of Computational Physics

Volume 95, Issue 1, July 1991, Pages 212-227

https://doi.org/10.1016/0021-9991(91)90260-R Get rights and content

Abstract

The use of a non-staggered computational grid for the numerical solutions of the incompressible flow equations has many advantages over the use of a staggered grid. A penalty, however, is inherent in the finite-difference approximations of the governing equations on non-staggered grids. In the primitive-variable solutions, the penalty is that the discrete continuity equation does not converge to machine accuracy. Rather it converges to a source term which is proportional to the fourth-order derivative of the pressure, the time increment, and the square of the grid spacing. An approach which minimizes the error in the discrete continuity equation is developed. Numerical results obtained for the driven cavity problem confirm the analytical developments.

References (12)

A.J. Chorin
J. Comput. Phys.
(1967)
S. Abdallah
J. Comput. Phys.
(1987)
S. Abdallah
J. Comput. Phys.
(1987)
F. Sotiropoulos et al.
J. Comput. Phys.
(1990)
J.C. Strikwerda et al.
J. Comput. Phys.
(1988)
F.H. Harlow et al.
Phys. Fluids
(1965)

There are more references available in the full text version of this article.

Cited by (80)

Analysis of factors affecting the internal seawater flow performance of pressure-retaining drilling tool for seafloor drill during drilling process
2023, Geoenergy Science and Engineering
In order to solve the low core collection rate during the natural gas hydrate pressure-retaining exploration, a pressure coring technology for seafloor drill is introduced, and the reasons for the decline in the core collection rate are analyzed. The mathematical model for predicting the water pressure change in the core barrel is established, and the design criteria of drainage path parameters of the drilling tool is formed and verified by engineering application. The results show that: when the valve spool opening x reaches 3.0 mm, the half cone angle α reaches 30°, and the equivalent annular area A₁ reaches 30 mm², the water pressure ΔP in the core barrel tends to be stable. Under the condition that the drilling speed is 20mm/s-50 mm/s and the ratio of the effective flow area of the screen is 0.2–1, the change of water pressure ΔP is little. When the valve spool opening is 4.0 mm, the velocity and pressure change of the seawater at the cone of the valve spool are stable. The engineering application in sediment formation shows that when adopting pressure coring method, the core collection rate exceeds 70% within the drilling speed of 20–30 mm/s.
Optimisation of drainage performance of the thin-walled core barrel sealing technology for pressure preservation sampling
2022, Ocean Engineering
Citation Excerpt :
Fig. 8 shows the velocity contour variation with different openings. Assuming that the minimum sectional area perpendicular to fluid velocity is the cross-section, according to the continuity equation, the sectional area at the cross-section is the smallest, and the fluid velocity is the largest (Sotiropoulos and Abdallah, 1991; DeGroff et al., 1998). It can be seen from Fig. 8 that when the opening is 2.0 mm, the maximum velocity appears in the cone area between the valve spool and the valve cover.
Pressure preservation sampling technology is a crucial method to assess natural gas hydrates. There are limitations to directly applying conventional pressure preservation technology for seafloor drills, such as large and complex structures. The working and design principle of the thin-walled core barrel sealing technique is introduced. To solve the problem of low core recovery percent during sediment drilling, the flow field characteristics of the valve spool in different working conditions are discussed. The area of the cross-section and the calculation method under the typical position are obtained. The valve spool structural parameters are determined and verified by a laboratory drilling test. The results show that: the larger the cone angle of the spool or the smaller the percentage of valve spool cone height to valve spool full cone height (the following is called P), the smaller the opening value when the cross-section area reaches a stable level. When the cone angle is 30, and the P is 20%, the cross-sectional area reaches stable after the opening is 4.5 mm. When the opening exceeds 3.0 mm, the flow rate changes little. When drilling in the sediment, the average coring recovery percent exceeds 70%.
Unstructured Cartesian refinement with sharp interface immersed boundary method for 3D unsteady incompressible flows
2016, Journal of Computational Physics
A novel numerical method is developed for solving the 3D, unsteady, incompressible Navier–Stokes equations on locally refined fully unstructured Cartesian grids in domains with arbitrarily complex immersed boundaries. Owing to the utilization of the fractional step method on an unstructured Cartesian hybrid staggered/non-staggered grid layout, flux mismatch and pressure discontinuity issues are avoided and the divergence free constraint is inherently satisfied to machine zero. Auxiliary/hanging nodes are used to facilitate the discretization of the governing equations. The second-order accuracy of the solver is ensured by using multi-dimension Lagrange interpolation operators and appropriate differencing schemes at the interface of regions with different levels of refinement. The sharp interface immersed boundary method is augmented with local near-boundary refinement to handle arbitrarily complex boundaries. The discrete momentum equation is solved with the matrix free Newton–Krylov method and the Krylov-subspace method is employed to solve the Poisson equation. The second-order accuracy of the proposed method on unstructured Cartesian grids is demonstrated by solving the Poisson equation with a known analytical solution. A number of three-dimensional laminar flow simulations of increasing complexity illustrate the ability of the method to handle flows across a range of Reynolds numbers and flow regimes. Laminar steady and unsteady flows past a sphere and the oblique vortex shedding from a circular cylinder mounted between two end walls demonstrate the accuracy, the efficiency and the smooth transition of scales and coherent structures across refinement levels. Large-eddy simulation (LES) past a miniature wind turbine rotor, parameterized using the actuator line approach, indicates the ability of the fully unstructured solver to simulate complex turbulent flows. Finally, a geometry resolving LES of turbulent flow past a complete hydrokinetic turbine illustrates the potential of the method to simulate turbulent flows past geometrically complex bodies on locally refined meshes. In all the cases, the results are found to be in very good agreement with published data and savings in computational resources are achieved.
Effect of vortex shedding in unsteady aerodynamic forces for a low Reynolds number stationary wing at low angle of attack
2016, Journal of Fluids and Structures
This study elucidates the relation between wake vortex shedding and aerodynamic force fluctuations for a low Reynolds number wing from time resolved particle image velocimetry (TR-PIV) experimental measurements. The results reveal a periodic lift and drag variation within the shedding cycle and resolve the frequencies of those fluctuations from a proper orthogonal decomposition (POD) and power spectral density (PSD) analysis. To show the effect of vortex shedding on the body force fluctuations, the evolution of instantaneous aerodynamic forces is compared to the pressure field of the fluid flow and to the vortical structures in the wake of the airfoil. A six step model describing the vortex-force relation is proposed. It shows that changes in lift such as maximum lift and minimum lift are associated with the detachment of a vortex. It also shows that the minimum or local minimum drag value is obtained at the onset formation of a vortex on the airfoil wake. Similarly, the maximum or local maximum drag is obtained at the onset formation of the saddle on the airfoil wake. The model further explains the asymmetry observed in the unsteady drag force evolution. The model can be used to optimize flow control and fluid-structure interaction applications.
Explicit filtering and exact reconstruction of the sub-filter stresses in large eddy simulation
2016, Journal of Computational Physics
Citation Excerpt :
These exchanges are important because they affect the small-scale flow dynamics, while spurious pressure modes are negligible on unstructured meshes [46]. Stabilisation of the P1 CG–P1 CG element pair is achieved by including a fourth-order pressure derivative term into the continuity equation [47,48]. The disadvantage is that the continuity equation is no longer exactly satisfied, i.e. the velocity field is not divergence-free.
Explicit filtering has the effect of reducing numerical or aliasing errors near the grid scale in large eddy simulation (LES). We use a differential filter, namely the inverse Helmholtz operator, which is readily applied to unstructured meshes. The filter is invertible, which allows the sub-filter scale (SFS) stresses to be exactly reconstructed in terms of the filtered solution. Unlike eddy viscosity models, the method of filtering and reconstruction avoids making any physical assumptions and is therefore valid in any flow regime. The sub-grid scale (SGS) stresses are not recoverable by reconstruction, but the second-order finite element method used here is an adequate source of numerical dissipation in lieu of an SGS model. Results for incompressible turbulent channel flow at ${Re}_{τ} = 180$ are presented which show that explicit filtering and exact SFS reconstruction is a significant improvement over the standard LES approach of implicit filtering and eddy-viscosity SGS modelling.
An overset grid method for integration of fully 3D fluid dynamics and geophysics fluid dynamics models to simulate multiphysics coastal ocean flows
2014, Journal of Computational Physics
Citation Excerpt :
The framework will be based on a hybrid of an F3DFD and a GFD model via a DDM implemented by an overset grid method. The solver for incompressible flow on overset meshes (SIFOM) is adopted as the F3DFD model, which was developed by us and co-workers to simulate small-scale, fully 3D, unsteady, incompressible flows with high-fidelity (e.g., [14,33,43–46]). The unstructured grid finite volume coastal ocean model (FVCOM) is selected as the GFD model, which is a popular model in the coastal ocean community, for large-scale background flows [6,47].
It is now becoming important to develop our capabilities to simulate coastal ocean flows involved with distinct physical phenomena occurring at a vast range of spatial and temporal scales. This paper presents a hybrid modeling system for such simulation. The system consists of a fully three dimensional (3D) fluid dynamics model and a geophysical fluid dynamics model, which couple with each other in two-way and march in time simultaneously. Particularly, in the hybrid system, the solver for incompressible flow on overset meshes (SIFOM) resolves fully 3D small-scale local flow phenomena, while the unstructured grid finite volume coastal ocean model (FVCOM) captures large-scale background flows. The integration of the two models are realized via domain decomposition implemented with an overset grid method. Numerical experiments on performance of the system in resolving flow patterns and solution convergence rate show that the SIFOM-FVCOM system works as intended, and its solutions compare reasonably with data obtained with measurements and other computational approaches. Its unparalleled capabilities to predict multiphysics and multiscale phenomena with high-fidelity are demonstrated by three typical applications that are beyond the reach of other currently existing models. It is anticipated that the SIFOM-FVCOM system will serve as a new platform to study many emerging coastal ocean problems.

View all citing articles on Scopus

View full text

The discrete continuity equation in primitive variable solutions of incompressible flow

Abstract

J. Comput. Phys.

J. Comput. Phys.

J. Comput. Phys.

J. Comput. Phys.

J. Comput. Phys.

Phys. Fluids