Time-dependent Navier-Stokes Equations Research Articles

The 3-D time-dependent Navier–Stokes equations on multi-connected domains with inhomogeneous boundary conditions and spectral hyperviscosity

We consider the time-dependent 3-D Navier–Stokes equations (NSE) on a multi-connected bounded domain Ω⊂R3\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\Omega \\subset \\mathbb {R^{\ extrm{3}}}$$\\end{document} with inhomogeneous boundary data β∈H1/2(Γ)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\beta \\in H^{1/2}(\\Gamma )$$\\end{document} on ∂Ω=Γ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\partial \\Omega =\\Gamma $$\\end{document}, where Γ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\Gamma $$\\end{document} is a union of Lipschitz continuous surfaces Γ0,Γ1,⋯,Γl\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\Gamma _{0},\\Gamma _{1},\\dots ,\\Gamma _{l}$$\\end{document}. This assumption includes the particular case when the Γi\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\Gamma _{i}$$\\end{document} are disjoint, the stationary version of which is classically known as Leray’s problem. Existence results for Leray’s problem have either assumed flux conditions beyond the general flux condition necessitated by compatibility constraints, or required size restrictions on the data. Here we incorporate a spectral hyperviscosity term in the time-dependent case and obtain existence and foundational results, assuming only the general flux condition and without imposing size restrictions on the boundary data β\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\beta $$\\end{document}. For any such β∈H1/2(Γ)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\beta \\in H^{1/2}(\\Gamma )$$\\end{document} we establish global existence and uniqueness of mild solutions. Then on any interval [0, T] on which these solutions and the corresponding NSE solution share a common H1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$H^{1}$$\\end{document}-bound (as is present on local intervals of existence of strong solutions, in certain special cases, or as is commonly assumed in achieving strong convergence results in numerical studies) we show for slightly more regular β\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\beta $$\\end{document} that the spectrally-hyperviscous solutions converge strongly and uniformly in H1(Ω)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$H^{1}(\\Omega )$$\\end{document} to the NSE solution as the spectral hyperviscosity term vanishes in the limit of key parameters. To achieve this robust sense of approximation of the NSE system, an involved setup and specially-adapted semigroup techniques assume essential roles, and the exposition is new for the case β=0\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\beta =0$$\\end{document} as well. Our final results adapt the NSE reformulation in Liu et al. (J Comput Phys 229(9):3428–3453, 2010) to recast our approximating system in a form potentially more adaptable to computation.

Superconvergence error analysis of linearized semi-implicit bilinear-constant SAV finite element method for the time-dependent Navier–Stokes equations

In this paper, based on the scalar auxiliary variable (SAV) approach, the superconvergence error analysis is investigated for the time-dependent Navier–Stokes equations. In which, an equivalent system of the Navier–Stokes equations with three variables and a fully-discrete scheme is developed with semi-implicit Euler discretization for the temporal direction and low-order bilinear-constant finite element discretization for the spatial direction, respectively. With the help of the high-precision estimations of the bilinear-constant finite element pair on the rectangular meshes, the superclose error estimates for velocity in H1-norm and pressure in L2-norm are obtained by treating the trilinear term carefully and skillfully. The global superconvergence results are also derived in terms of a simple and efficient interpolation post-processing technique. Finally, some numerical results are provided to demonstrate the correctness of the theoretical analysis.

Full discretization of the time dependent Navier–Stokes equations with anisotropic slip boundary condition coupled with the convection–diffusion–reaction equation

Full discretization of the time dependent Navier–Stokes equations with anisotropic slip boundary condition coupled with the convection–diffusion–reaction equation

Semi-Lagrangian finite element exterior calculus for incompressible flows

We develop a semi-Lagrangian discretization of the time-dependent incompressible Navier-Stokes equations with free boundary conditions on arbitrary simplicial meshes. We recast the equations as a nonlinear transport problem for a momentum 1-form and discretize in space using methods from finite element exterior calculus. Numerical experiments show that the linearly implicit fully discrete version of the scheme enjoys excellent stability properties in the vanishing viscosity limit and is applicable to inviscid incompressible Euler flows. We obtain second-order convergence and conservation of energy is achieved through a Lagrange multiplier.

A Performance Simulation Methodology for a Whole Turboshaft Engine Based on Throughflow Modelling

To accurately predict the matching relationships between the various components and the engine performance in the whole aero-engine environment, this study introduces a two-dimensional throughflow simulation method for the whole aero-engine. This method is based on individual throughflow solvers for the turbo-machinery and the combustor. It establishes a throughflow simulation model for the whole engine by integrating with the compressor-turbine co-operating equations and boundary conditions. The turbo-machinery throughflow solver employs a circumferentially averaged form of the time-dependent Navier–Stokes equations (N-S) as the governing equation. The combustor solver uses the Reynolds Average Navier–Stokes (RANS) method to solve flow and chemical reaction processes by constructing turbulence, combustion, and radiation models. The accuracy of the component solver is validated using Pratt and Whitney’s three-stage axial compressor (P&W3S1) and General Electric’s high-pressure turbine (GE-EEE HPT), and the predicted results are consistent with the experimental data. Finally, the developed throughflow method is applied to simulate the throttling characteristics of the WZ-X turboshaft engine. The results predicted by the throughflow program are consistent with the GasTurb calculations, including the trends of shaft power delivered, specific fuel consumption (SFC), inlet airflow, and total pressure ratio of the compressor. The developed method to perform throughflow simulation of the whole aero-engine eliminates the dependence on a general component map. It can quickly obtain the meridian flow field parameters and overall engine characteristics, which is expected to guide the design and modification of the engine in the future.

A numerical study of interface dynamics in fluid materials

This paper deals with the approximation of the dynamics of two fluids having non-matching densities and viscosities. The modeling involves the coupling of the Allen-Cahn equation with the time-dependent Navier-Stokes equations. The Allen-Cahn equation describes the evolution of a scalar order parameter that assumes two distinct values in different spatial regions. Conversely, the Navier-Stokes equations govern the movement of a fluid subjected to various forces like pressure, gravity, and viscosity. When the Allen-Cahn equation is coupled with the Navier-Stokes equations, it is typically done through a surface tension term. The surface tension term accounts for the energy required to create an interface between the two phases, and it is proportional to the curvature of the interface. The Navier-Stokes equations are modified to include this term, which leads to the formation of a dynamic interface between the two phases. The resulting system of equations is known as the two-phase Navier-Stokes/Allen-Cahn equations. In this paper, the authors propose a mathematical model that combines the Allen-Cahn model and the Navier-Stokes equations to simulate multiple fluid flows. The Allen-Cahn model is utilized to represent the diffuse interface between different fluids, while the Navier-Stokes equations are employed to describe the fluid dynamics. The Allen-Cahn-Navier-Stokes model has been employed to simulate the generation of bubbles in a liquid subjected to an acoustic field. The model successfully predicted the size of the bubbles and the frequency at which they formed. The numerical outcomes were validated against experimental data, and a favorable agreement was observed.

A Mixed Finite Element Approximation for Time-Dependent Navier–Stokes Equations with a General Boundary Condition

In this paper, we present a numerical scheme for addressing the unsteady asymmetric flows governed by the incompressible Navier–Stokes equations under a general boundary condition. We utilized the Finite Element Method (FEM) for spatial discretization and the fully implicit Euler scheme for time discretization. In addition to the theoretical analysis of the error in our numerical scheme, we introduced two types of a posteriori error indicators: one for time discretization and another for spatial discretization, aimed at effectively controlling the error. We established the equivalence between these estimators and the actual error. Furthermore, we conducted numerical simulations in two dimensions to assess the accuracy and effectiveness of our scheme.

Numerical Modeling of Aircraft Fire: Postcrash Fire

The Fire Dynamics Simulator (FDS) model was utilized in this study to replicate a full-scale aircraft postcrash experiment conducted within the C-133 test facility by the Federal Aviation Administration. FDS is a computational fire field model that incorporates submodels for soot formation, pyrolysis, and thermal radiation transport. It solves three-dimensional time-dependent Navier–Stokes equations and is grounded in the large-eddy simulation approach and the eddy dissipation concept, serving as turbulence and combustion models. The obtained results, including the heat release rate and temperature, were validated against experimental data and compared with earlier prediction studies employing different turbulence and combustion models. The results from this simulation closely align with the experiment’s findings. The impact of fire-blocking layers and carry-on baggage on interior material was examined. Moreover, two boundary conditions were imposed on the fuselage structure: 1) the adiabatic wall, and 2) heat loss within the wall. Both the fire-blocking layers and the adiabatic boundary condition played a significant role in the flashover occurrence. The large-eddy simulation and eddy dissipation concept approaches have demonstrated a reliable ability to predict flashover and general fire properties to a considerable extent.

Numerical analysis of the SAV scheme for the EMAC formulation of the time-dependent Navier–Stokes equations

In this work, motivated by the scalar auxiliary variable (SAV) method, we construct the first-order implicit–explicit type scheme for the EMAC formulation of the time-dependent Navier–Stokes equations. This scheme is linear and only requires solving a series of Stokes type equations with constant coefficients at each time step. What is more, the stability is given out without any condition on the time step, where the Gronwall’s constant does not dependent on the Reynolds number. Especially, this scheme also conserves momentum and angular momentum in the fully discrete case with only weakly enforced divergence constraint. Furthermore, we give out the priori error estimates for the velocity. Finally, we carry out several numerical tests to confirm the theoretical results and the effectiveness of the proposed scheme.

Numerical Flowfield Visualization Over Heat Shield Of Satellite Launch Vehicle at Mach 0.8 - 3.0

The paper presents a numerical flowfield visualization over a typical payload shroud of a satellite launch vehicle at Mach number range of 0.80 -3.0 and Reynolds number range of 33 x 106 - 46 x 106 /m. The numerical simulation over the bulbous heat shield is carried out by solving time-dependent compressible axisymmetric turbulent Reynolds-averaged NavierStokes equations. The closure of these equations is obtained employing the Baldwin-Lomax turbulence model. A three-stage Runge-Kutta time-stepping scheme has been used in conjunction with finite-volume discretization of the computational domain. The flowfield features over the bulbous heat shield have been analyzed from the velocity vector and the density contour plots. Flow separation on the payload shroud due to the normal shock wave is observed at Mach number 0.80 and 0.90. The normal shock movement on the heat shield is simulated for various transonic Mach number. A recirculation zone of flowfield is formed in the boat-tail region of the heat shield. The vorticity formation ahead of the heat shield moves close to the heat shield with the increasing transonic Mach number, which is observed in the density contours plots of the flowfield. Shock stand-off distance at the supersonic speed is calculated and compared with asymptotic formula of Frank and Zierep. They are found in good agreement.

Flow Field Simulations Over Reentry Modules at High Speed

Reentry capsule configurations significantly differ from each other due to entry conditions and mission requirements. This paper describes numerical simulations of the viscous flow over the Beagle and the OREX (Orbital Reentry EXperiments) configurations for freestream Mach numbers in the range of 1.2 - 5.0. The flow fields over the reentry modules are obtained by solving time-dependent axisymmetric compressible laminar Navier-Stokes equations. The fluid mechanics equations are discretized in spatial coordinates employing a finite volume method, which reduces the governing equations to semi-discretized ordinary differential equations. Temporal integration is carried out using a two-stage Runge-Kutta time-stepping scheme. A local time-stepping is employed to get the steady state solution. The numerical simulation is performed on a single-block structured grid. The flow field features around the reentry capsules such as bow shock wave, sonic line, expansion fan and recirculation flow in the base-shell region are well captured by the present numerical computation. The effects of the geometrical parameters of the module, such as aspect ratio, frontal segment bluntness, fore body cone angle and shoulder rounding radius on the wake throat width, distance from the wake throat to the base of the model, flow departure angle and aerodynamic drag are analyzed using the numerically simulated flow field.

Aerothermodynamic Analysis of Space Recovery Experiment Capsule at Mach 6

The main focus of this paper is to compute the flowfield and assess the aerodynamic drag of Space Recovery Experiment (SRE) modules at freestream Mach 6 at zero angle of incidence. Time-dependent Compressible axisymmetric Navier-Stokes equations are numerically simulated using a finite volume discretization method with a three-stage Runge-Kutta time-stepping scheme. The computed shock stand-off distance, pressure ratio across normal shock and stagnation point heat flux are showing satisfactory results with the analytical solutions.Available wind tunnel data such as Schlieren images, oil flow pictures and fore-body aerodynamic drag are compared with numerical results.Influence of various ramp angles of the SRE modules on surface pressure, flow separation, skin friction coefficient, wall heat flux and fore and base aerodynamic drag coefficients are investigated using present numerical data.

Numerical analysis of projection methods for the time-dependent Navier-Stokes equations with modular grad-div stabilization

Numerical analysis of projection methods for the time-dependent Navier-Stokes equations with modular grad-div stabilization

Pressure robust SUPG-stabilized finite elements for the unsteady Navier–Stokes equation

Abstract In the present contribution, we propose a novel conforming finite element scheme for the time-dependent Navier–Stokes equation, which is proven to be both convection quasi-robust and pressure robust. The method is built combining a ‘divergence-free’ velocity/pressure couple (such as the Scott–Vogelius element), a discontinuous Galerkin in time approximation and a suitable streamline upwind Petrov–Galerkin-curl stabilization. A set of numerical tests, in accordance with the theoretical results, is included.

Unconditional Superconvergence Error Estimates of Semi-Implicit Low-Order Conforming Mixed Finite Element Method for Time-Dependent Navier–Stokes Equations

In this paper, the unconditional superconvergence error analysis of the semi-implicit Euler scheme with low-order conforming mixed finite element discretization is investigated for time-dependent Navier–Stokes equations. In terms of the high-accuracy error estimates of the low-order finite element pair on the rectangular mesh and the unconditional boundedness of the numerical solution in L∞-norm, the superclose error estimates for velocity in H1-norm and pressure in L2-norm are derived firstly by dealing with the trilinear term carefully and skillfully. Then, the global superconvergence results are obtained with the aid of the interpolation post-processing technique. Finally, some numerical experiments are carried out to support the theoretical findings.

Enhancing the viscosity-splitting method to solve the time-dependent Navier–Stokes equations

The viscosity-splitting (VS) scheme for solving the incompressible time-dependent Navier–Stokes equations consists to decouple the full problem into two easy subproblems to resolve. This method splits the nonlinearity (convection term) and the incompressibility into two different steps. As a result, every subproblem is easier to solve than the Navier–Stokes system. The semi-implicit approach is mostly used in the first step for the nonlinearity term and the intermediate velocity is calculated by dropping the pressure gradient. We then resolve a linear Stokes problem in the next step.We present a modified viscosity-splitting scheme (Incremental viscosity-splitting) (IVS) to improve this method. This scheme adds the gradient of the pressure (calculated at the previous time step) in the first step and modifies the Stokes equations in the second step to predict a better value of the intermediate velocity. A manufactured solution and two well-known benchmarks will be used to demonstrate the effectiveness of the suggested strategy.

A data-driven surrogate modeling approach for time-dependent incompressible Navier-Stokes equations with dynamic mode decomposition and manifold interpolation

This work introduces a novel approach for data-driven model reduction of time-dependent parametric partial differential equations. Using a multi-step procedure consisting of proper orthogonal decomposition, dynamic mode decomposition, and manifold interpolation, the proposed approach allows to accurately recover field solutions from a few large-scale simulations. Numerical experiments for the Rayleigh-Bénard cavity problem show the effectiveness of such multi-step procedure in two parametric regimes, i.e., medium and high Grashof number. The latter regime is particularly challenging as it nears the onset of turbulent and chaotic behavior. A major advantage of the proposed method in the context of time-periodic solutions is the ability to recover frequencies that are not present in the sampled data.

Modeling and analysis of dynamic load factor in a fail-safe crane mechanism equipped with a hydraulic equalizing cylinder

In the design and manufacturing of sensitive crane mechanisms, in particular, those utilized in nuclear facilities, safety is of utmost importance. A few designs are based on a hydraulic equalizing cylinder where the purpose is mainly to dampen the dynamic load factor generated by the sudden failure of the wire rope system. To the best of the authors' knowledge, design equations and engineering models for such an intricate dynamic mechanism are not reported anywhere in the literature. Toward achieving this goal, a model is proposed for the dynamic behavior of a hydraulic equalizing cylinder when coupled to a single-failure-proof cross-reeved crane hoist. The constitutive equations involve the time-dependent incompressible Navier-Stokes equations. These equations are coupled with the dynamics of the lower block and the elastic wire rope. A nonlinear second-order ordinary differential equation is obtained, which is solved analytically. This yields a closed-form relation for the dynamic load factor. The validity of the new model is investigated using experimental studies, which are presented in this work. The normalized root mean square error for the dynamic load factor obtained by the model compares favorably to the experimental results, where the maximum observed error is 3.9 percent. The results confirm that the current model can be used as a design equation by practitioners.

Convergence of the MAC scheme for the incompressible Navier-Stokes equations with variable density and viscosity

The present paper addresses the convergence of the implicit Marker-and-Cell scheme for time-dependent Navier–Stokes equations with variable density and density-dependent viscosity and forcing term. A priori estimates on the unknowns are obtained, and thanks to a topological degree argument, they lead to the existence of an approximate solution at each time step. Then, by compactness arguments relying on these same estimates, we obtain the convergence (up to the extraction of a subsequence), when the space and time steps tend to zero, of the numerical solutions to a limit; this latter is shown to be a weak solution to the continuous problem by passing to the limit in the scheme.

Optimal tubes for non-cylindrical Navier-Stokes flows with Navier boundary condition

We consider a shape optimization problem where state equations are the time-dependent Navier-Stokes equations with Navier boundary condition occupying a non-cylindrical domain. In the first part of our work, we prove existence of weak solutions for the state equations using the elliptic regularization method. In the second part, we show stability of the state problem solutions with respect to admissible shape variations and derive an existence result for the considered shape minimization problem.