The Incompressible Navier–Stokes for the Nonlinear Discrete Velocity Models

Abstract

We establish the incompressible Navier–Stokes limit for the discrete velocity model of the Boltzmann equation in any dimension of the physical space, for densities which remain in a suitable small neighborhood of the global Maxwellian. Appropriately scaled families solutions of discrete Boltzmann equation are shown to have fluctuations that locally in time converge strongly to a limit governed by a solution of Incompressible Navier–Stokes provided that the initial fluctuation is smooth, and converges to appropriate initial data. As applications of our results, we study the Carleman model and the one-dimensional Broadwell model.     

A deconvolution enhancement of the Navier–Stokes-αβ model

A deconvolution enhancement of the Navier–Stokes-αβ model for turbulent flow is introduced. The energy and energy-dissipation rate for the enhanced model are derived. It is also shown that the consistency error, relative to the Navier–Stokes equations, and the microscale of the enhanced model are less than those of the Navier–Stokes-αβ model. The proposed model is used to simulate the Taylor–Green vortex problem and results show a qualitatively improved representation of the mean-square vorticity when compared to the Navier–Stokes-αβ model. Numerical studies of the energy spectrum and the alignment between the vorticity and the eigenvectors of the stretching tensor for three-dimensional turbulent flows with Re = 200 are used to explore the utility of the model. A benchmark problem of a two-dimensional channel flow over a step for Re = 600 also indicates that this model can be applied to more general flows than those involving periodic boundary conditions.     

Applications of a Cole–Hopf transform to the 3D Navier–Stokes equations

The Navier–Stokes equations written in the vector potential can be recast as non-linear Schrödinger equations at imaginary times, i.e. heat equations with a potential term, using the Cole–Hopf transform. On this basis, we study two kinds of Navier–Stokes flows by means of direct numerical simulations. In an experiment on vortex reconnection, it is found that the potential term takes large negative values in regions where intensive reconnection takes place, whereas the signature of the non-linear term is more broadly spread. For decaying turbulence starting from a random initial condition, such a correspondence is also observed in the early stage when the flow is dominated by vorticity layers. At later times, when the flow features several tubular vortices, this correspondence becomes weaker. Finally, a similar set of transformations is presented for the magneto–hydrodynamic equations, which reduces them to a set of heat equations with suitable potential terms, thereby obtaining new criteria for the regularity of their solutions.     

Self-oscillations of electronic and nuclear polarization during optical excitation of an Si/CaF2 nanostructure in a transverse magnetic field

We examine theoretically low-frequency and high-frequency self-oscillations of electronic and nuclear polarization in an Si/CaF₂ nanostructure in a transverse magnetic field. We show that the low-frequency self-oscillations are stable in zero field, and the analogous high-frequency oscillations are stable beyond the region of the maximum on the Hanle curve. The frequency of the low-frequency oscillations is 0.001–0.500 of the reciprocal nuclear longitudinal relaxation time; the frequency of the high-frequency oscillations is 10⁸–10⁹ Hz, and their amplitude reaches 50% of the initial electronic spin polarization. __________ Translated from Zhurnal Prikladnoi Spektroskopii, Vol. 73, No. 3, pp. 363–369, May–June, 2006.     

Numerical comparisons of time–space iterative method and spatial iterative methods for the stationary Navier–Stokes equations

This paper makes some numerical comparisons of time–space iterative method and spatial iterative methods for solving the stationary Navier–Stokes equations. The time–space iterative method consists in solving the nonstationary Stokes equations based on the time–space discretization by the Euler implicit/explicit scheme under a weak uniqueness condition (A2). The spatial iterative methods consist in solving the stationary Stokes scheme, Newton scheme, Oseen scheme based on the spatial discretization under some strong uniqueness assumptions. We compare the stability and convergence conditions of the time–space iterative method and the spatial iterative methods. Moreover, the numerical tests show that the time–space iterative method is the more simple than the spatial iterative methods for solving the stationary Navier–Stokes problem. Furthermore, the time–space iterative method can solve the stationary Navier–Stokes equations with some small viscosity and the spatial iterative methods can only solve the stationary Navier–Stokes equations with some large viscosities.     

MC-Smooth: A mass-conserving,smooth vorticity–velocity formulation for multi-dimensional flows

A novel vorticity–velocity formulation of the Navier–Stokes equations – the Mass-Conserving, Smooth (MC-Smooth) vorticity–velocity formulation – is developed in this work. The governing equations of the MC-Smooth formulation include a new second-order Poisson-like elliptic velocity equation, along with the vorticity transport equation, the energy conservation equation, and N_spec species mass balance equations. In this study, the MC-Smooth formulation is compared to two pre-existing vorticity–velocity formulations by applying each formulation to confined and unconfined axisymmetric laminar diffusion flame problems. For both applications, very good to excellent agreement for the simulation results of the three formulations has been obtained. The MC-Smooth formulation requires the least CPU time and can overcome the limitations of the other two pre-existing vorticity–velocity formulations by ensuring mass conservation and solution smoothness over a broader range of flow conditions. In addition to these benefits, other important features of the MC-Smooth formulation include: (1) it does not require the use of a staggered grid, and (2) it does not require excessive grid refinement to ensure mass conservation. The MC-Smooth formulation is a computationally attractive approach that can effectively extend the applicability of the vorticity–velocity formulation.     

Laminar incompressible viscous flow past a cylinder performing rotary oscillations

The numerical simulation of the laminar viscous flow past a cylinder performing rotary oscillations around its axis is carried out. The Navier–Stokes equations are solved by finite volume method using the program package OpenFOAM. The values of the amplitude and frequency of forced oscillations are found, at which the maximum reduction of the drag coefficient of the cylinder is achieved.     

Numerical solution for viscous flow for two-dimensional domains using orthogonal coordinate systems

A method to generate two-dimensional orthogonal grids in simply and doubly connected domains is given. The method not only generates the grids but also finds the modulus of the domain simultaneously. Also, the Navier0Stokes equations are solved is some doubly connected domains. The insteady vorticity stream function approach is used. The stream function on one of the boundaries has to be updated at every time step. Comparisons are made between numerical and experimental results; quite good agreement can be achieved.     

A gradient stable scheme for a phase field model for the moving contact line problem

In this paper, an efficient numerical scheme is designed for a phase field model for the moving contact line problem, which consists of a coupled system of the Cahn–Hilliard and Navier–Stokes equations with the generalized Navier boundary condition [1], [2], [4]. The nonlinear version of the scheme is semi-implicit in time and is based on a convex splitting of the Cahn–Hilliard free energy (including the boundary energy) together with a projection method for the Navier–Stokes equations. We show, under certain conditions, the scheme has the total energy decaying property and is unconditionally stable. The linearized scheme is easy to implement and introduces only mild CFL time constraint. Numerical tests are carried out to verify the accuracy and stability of the scheme. The behavior of the solution near the contact line is examined. It is verified that, when the interface intersects with the boundary, the consistent splitting scheme [21], [22] for the Navier Stokes equations has the better accuracy for pressure.     

Fourier spectral and wavelet solvers for the incompressible Navier–Stokes equations with volume-penalization: Convergence of a dipole–wall collision

In this study, we use volume-penalization to mimic the presence of obstacles in a flow or a domain with no-slip boundaries. This allows in principle the use of fast Fourier spectral methods and coherent vortex simulation techniques (based on wavelet decomposition of the flow variables) to compute turbulent wall-bounded flow or flows around solid obstacles by simply adding one term in the equation. Convergence checks are reported using a recently revived, and unexpectedly difficult dipole–wall collision as a benchmark computation. Several quantities, like the vorticity isolines, truncation error, kinetic energy and enstrophy are inspected for a collision of a dipole with a no-slip wall and compared with available benchmark data obtained with a standard Chebyshev pseudospectral method. We quantify the possible deteriorating effects of the Gibbs phenomenon present in the Fourier based schemes due to continuity restrictions of the penalized Navier–Stokes equations on the wall. It is found that Gibbs oscillations have a negligible effect on the flow evolution allowing higher-order recovery of the accuracy on a Fourier basis by means of postprocessing. An advantage of coherent vortex simulations, on the other hand, is that the degrees of freedom of the flow computation can strongly be reduced. In this study, we quantify the possible reduction of degrees of freedom while keeping the accuracy. For an optimal convergence scenario the penalization parameter has to scale with the number of Fourier and wavelet modes. In addition, an implicit treatment of the Darcy drag term in the penalized Navier–Stokes equations is beneficial since this allows one to set the time step independent from the penalization parameter without additional computational or memory requirements.     

Modified periodic-shell boundary condition for dissipative particle dynamics simulations

We have developed the modified periodic-shell boundary condition (BC) for dissipative particle dynamics (DPD) simulations, which enables us to simulate an outer flow problem around an obstacle using a small simulation region. In order to clarify the validity of this BC, we have treated a uniform flow past a circular cylinder. The present BC has been compared with the ordinary BC such as the uniform flow condition. Also, the present results have been compared with those of the numerical results of the Navier–Stokes equation. The ordinary uniform BC is seen to give rise to significantly distorted flow fields and also to significant disappearance of dissipative particles from the simulation region. In contrast, for the present modified periodic-shell BC, the number density of dissipative particles is kept almost constant during a simulation run, and the flow field is in reasonable agreement with the result, which has been obtained by numerical simulations of the Navier–Stokes equation.     

星地下行孔径接收闪烁频谱的理论研究

 基于孔径接收的对数振幅时间相关函数,结合星地激光下行链路的实际情况,推导了适用于星地下行孔径接收的闪烁频谱表达式。进而分析了孔径尺度和湍流外尺度对闪烁频谱的影响。结果表明,下行孔径接收闪烁频谱分为低频段和高频段两个区间,低频段为常数,高频段呈幂指数规律变化,且幂律的绝对值为11/3。当接收孔径增加时,频谱的幅值和宽度逐渐减小,其形状保持不变。随着等效湍流外尺度的增加,低频谱的幅值逐渐增大,高频谱保持不变。     

Collective oscillations of the magnetic moments of a chain of spherical magnetic nanoparticles with uniaxial magnetic anisotropy

Magnetic particles moving freely in a fluid can organize dense phases (3D clusters or linear chains). We analyze the spectrum of magnetic oscillations of a chain of spherical magnetic particles taking into account the magnetic anisotropy of an individual particle for an arbitrary relation between the anisotropy energy and the energy of the dipole interaction of particles. For any relation between these energies, the spectrum contains three branches of collective oscillations: a high-frequency branch and a weakly split doublet of low-frequency branches. The frequency of the high-frequency branch is determined by a stronger interaction, while the frequencies of the low-frequency branches are determined by the weakest interaction. Accordingly, the dispersion is maximal for oscillations formed by the dipole-dipole interaction of particles, which have high frequencies in the case of a strong dipole interaction or low frequencies in the case of a strong anisotropy.     

Analogue of the Cole-Hopf transform for the incompressible Navier–Stokes equations and its application

We consider the Navier–Stokes equations written in the stream function in two dimensions and vector potentials in three dimensions, which are critical dependent variables. On this basis, we introduce an analogue of the Cole-Hopf transform, which exactly reduces the Navier–Stokes equations to the heat equations with a potential term (i.e. the nonlinear Schrödinger equation at imaginary times). The following results are obtained. (i) A regularity criterion immediately obtains as the boundedness of condition for the potential term when the equations are recast in a path-integral form by the Feynman-Kac formula. (ii) This in turn gives an additional characterisation of possible singularities for the Navier–Stokes equations. (iii) Some numerical results for the two-dimensional Navier–Stokes equations are presented to demonstrate how the potential term captures near-singular structures. Finally, we extend this formulation to higher dimensions, where the regularity issues are markedly open.     

Conservative integral form of the incompressible Navier–Stokes equations for a rapidly pitching airfoil

This study provides a simple moving-grid scheme which is based on a modified conservative form of the incompressible Navier–Stokes equations for flow around a moving rigid body. The modified integral form is conservative and seeks the solution of the absolute velocity. This approach is different from previous conservative differential forms [1], [2], [3] whose reference frame is not inertial. Keeping the reference frame being inertial results in simpler mathematical derivation to the governing equation which includes one dyadic product of velocity vectors in the convective term, whereas the previous [2], [3] needs to obtain the time derivative with respect to non-inertial frames causing an additional dyadic product in the convective term. The scheme is implemented in a second-order accurate Navier–Stokes solver and maintains the order of the accuracy. After this verification, the scheme is validated for a pitching airfoil with very high frequencies. The simulation results match very well with the experimental results [4], [5], including vorticity fields and a net thrust force. This airfoil simulation also provides detailed vortical structures near the trailing edge and time-evolving aerodynamic forces that are used to investigate the mechanism of the thrust force generation and the effects of the trailing edge shape. The developed moving-grid scheme demonstrates its validity for a rapid oscillating motion.     

Large-eddy simulation of impinging jets with embedded azimuthal vortices

We have carried out large-eddy simulations of an impinging jet with embedded azimuthal vortices, a model of the wake of a helicopter hovering in ground effect. The azimuthal vortices are generated by sinusoidal forcing of the velocity at the jet exit. They strengthen while they are advected towards the ground; when they are close to the solid surface, a layer of opposite-sign vorticity is formed at the wall, and lifted up to form a secondary vortex that interacts with the primary one. Regions of reversed flow are caused by the strong, localised, adverse pressure gradient. After this interaction, the primary vortices begin to decay, mostly due to the Reynolds shear stresses, which contribute to the turbulent diffusion of vorticity term in the budget of the phase-averaged azimuthal vorticity. This mechanism is extremely robust, and plays the most important role in the vortex decay even if no turbulence is initially present in the jet, or if the no-slip condition is removed. A three-dimensional instability also plays a role: removing it leads to slower decay. Our results also point out some challenges for turbulence models for the unsteady Reynolds-averaged Navier–Stokes equations.     

h-Multigrid for space-time discontinuous Galerkin discretizations of the compressible Navier–Stokes equations

Being implicit in time, the space-time discontinuous Galerkin discretization of the compressible Navier–Stokes equations requires the solution of a non-linear system of algebraic equations at each time-step. The overall performance, therefore, highly depends on the efficiency of the solver. In this article, we solve the system of algebraic equations with a h-multigrid method using explicit Runge–Kutta relaxation. Two-level Fourier analysis of this method for the scalar advection–diffusion equation shows convergence factors between 0.5 and 0.75. This motivates its application to the 3D compressible Navier–Stokes equations where numerical experiments show that the computational effort is significantly reduced, up to a factor 10 w.r.t. single-grid iterations.     

Approximation of the hydrostatic Navier–Stokes system for density stratified flows by a multilayer model: Kinetic interpretation and numerical solution

We present a multilayer Saint-Venant system for the numerical simulation of free surface density-stratified flows over variable topography. The proposed model formally approximates the hydrostatic Navier–Stokes equations with a density that varies depending on the spatial and temporal distribution of a transported quantity such as temperature or salinity. The derivation of the multilayer model is obtained by a Galerkin-type vertical discretization of the Navier–Stokes system with piecewise constant basis functions. In contrast with classical multilayer models in the literature that assume immiscible fluids, we allow here for mass exchange between layers. We show that the multilayer system admits a kinetic interpretation, and we use this result to formulate a robust finite volume scheme for its numerical approximation. Several numerical experiments are presented, including simulations of wind-driven stratified flows.     

An effective explicit pressure gradient scheme implemented in the two-level non-staggered grids for incompressible Navier–Stokes equations

In this paper, an improved two-level method is presented for effectively solving the incompressible Navier–Stokes equations. This proposed method solves a smaller system of nonlinear Navier–Stokes equations on the coarse mesh and needs to solve the Oseen-type linearized equations of motion only once on the fine mesh level. Within the proposed two-level framework, a prolongation operator, which is required to linearize the convective terms at the fine mesh level using the convergent Navier–Stokes solutions computed at the coarse mesh level, is rigorously derived to increase the prediction accuracy. This indispensable prolongation operator can properly communicate the flow velocities between the two mesh levels because it is locally analytic. Solution convergence can therefore be accelerated. For the sake of numerical accuracy, momentum equations are discretized by employing the general solution for the two-dimensional convection–diffusion–reaction model equation. The convective instability problem can be simultaneously eliminated thanks to the proper treatment of convective terms. The converged solution is, thus, very high in accuracy as well as in yielding a quadratic spatial rate of convergence. For the sake of programming simplicity and computational efficiency, pressure gradient terms are rigorously discretized within the explicit framework in the non-staggered grid system. The proposed analytical prolongation operator for the mapping of solutions from the coarse to fine meshes and the explicit pressure gradient discretization scheme, which accommodates the dispersion-relation-preserving property, have been both rigorously justified from the predicted Navier–Stokes solutions.     

On Numerical Simulation of Sound Damping Mechanisms in the Cell of a Sound-Absorbing Structure

Numerical simulation of acoustic processes with an interferometer at high acoustic pressure levels is one of the ways to noise reduction processes using samples of sound-absorbing structures (SAS). In the course of research, an SAS sample consisting of a single Helmholtz resonator of circular shape was used. Studies were conducted at sound pressure levels of 110, 130, 140, and 150 dB. Results obtained with the interferometer with normal wave incidence were taken as a basis. In the calculations, systems of linearized and complete Navier–Stokes equations were used. The results obtained with the linearized Navier–Stokes equations make it possible to determine the acoustic characteristics of the sample for linear modes of operation with sufficient accuracy. The complete system of Navier–Stokes equations, taking into account compressibility, allowed good qualitative and quantitative agreement with the experiment at high acoustic pressure levels.