Search

In a recent work (Castanon Quiroz & Di Pietro (2020) A hybrid high-order method for the incompressible Navier–Stokes problem robust for large irrotational body forces. Comput. Math. Appl., 79, 2655–2677), we have introduced a pressure-robust hybrid high-order method for the numerical solution of the incompressible Navier–Stokes equations on matching simplicial meshes. Pressure-robust methods are characterized by error estimates for the velocity that are fully independent of the pressure. A crucial question was left open in that work, namely whether the proposed construction could be extended to general polytopal meshes. In this paper, we provide a positive answer to this question. Specifically, we introduce a novel divergence-preserving velocity reconstruction that hinges on the solution inside each element of a mixed problem on a subtriangulation, then use it to design discretizations of the body force and convective terms that lead to pressure robustness. An in-depth theoretical study of the properties of this velocity reconstruction, and their reverberation on the scheme, is carried out for arbitrary polynomial degrees $k\geq 0$ and meshes composed of general polytopes. The theoretical convergence estimates and the pressure robustness of the method are confirmed by an extensive panel of numerical examples.

4D Flow Magnetic Resonance Imaging (MRI) is the state-of-the-art technique to comprehensively measure the complex spatio-temporal and multidirectional patterns of blood flow. However, it is subject to artifacts such as noise and aliasing, which due to the 3D and dynamic structure is difficult to detect in clinical practice. In this work, a new mathematical and computational model to determine the quality of 4D Flow MRI is presented. The model is derived by assuming the true velocity satisfies the incompressible Navier–Stokes equations and that can be decomposed by the measurements plus an extra field . Therefore, a non-linear problem with as unknown arises, which serves as a measure of data quality. A stabilized finite element formulation tailored to this problem is proposed and analyzed. Then, extensive numerical examples—using synthetic 4D Flow MRI data as well as real measurements on experimental phantom and subjects—illustrate the ability to use for assessing the quality of 4D Flow MRI measurements over space and time.

This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY); International audience; An Eulerian method to numerically solve incompressible bifluid problems with high density ratio is presented. This method can be considered as an improvement of the Ghost Fluid method, with the specificity of a sharp second-order numerical scheme for the spatial resolution of the discontinuous elliptic problem for the pressure. The Navier–Stokes equations are integrated in time with a fractional step method based on the Chorin scheme and discretized in space on a Cartesian mesh. The bifluid interface is implicitly represented using a level-set function. The advantage of this method is its simplicity to implement in a standard monofluid Navier–Stokes solver while being more accurate and conservative than other simple classical bifluid methods. The numerical tests highlight the improvements obtained with this sharp method compared to the reference standard first-order methods.

This paper studies fully discrete approximations to the evolutionary Navier{ Stokes equations by means of inf-sup stable H1-conforming mixed nite elements with a grad-div type stabilization and the Euler incremental projection method in time. We get error bounds where the constants do not depend on negative powers of the viscosity. We get the optimal rate of convergence in time of the projection method. For the spatial error we get a bound O(hk) for the L2 error of the velocity, k being the degree of the polynomials in the velocity approximation. We prove numerically that this bound is sharp for this method., MINECO grant MTM2016-78995-P (AEI), Junta de Castilla y León grant VA024P17, Junta de Castilla y León grant VA105G18, MINECO grant MTM2015-65608-P

We study the Cauchy problem of the 3D incompressible Navier–Stokes equations with nonlinear damping term $\alpha|\mathbf{u}|^{\beta-1}\mathbf{u}\ (\alpha>0\ \text{and}\ \beta\geq1)$. In [J. Math. Anal. Appl. 377(2011), 414–419], Zhang et al. obtained global strong solution for $\beta>3$ and the solution is unique provided that $33$.

We present a numerical method specifically designed for simulating three-dimensional fluid--structure interaction (FSI) problems based on the reference map technique (RMT). The RMT is a fully Eulerian FSI numerical method that allows fluids and large-deformation elastic solids to be represented on a single fixed computational grid. This eliminates the need for meshing complex geometries typical in other FSI approaches, and greatly simplifies the coupling between fluid and solids. We develop the first three-dimensional implementation of the RMT, parallelized using the distributed memory paradigm, to simulate incompressible FSI with neo-Hookean solids. As part of our new method, we develop a new field extrapolation scheme that works efficiently in parallel. Through representative examples, we demonstrate the method's accuracy and convergence, as well as its suitability in investigating many-body and active systems. The examples include settling of a mixture of heavy and buoyant soft ellipsoids, lid-driven cavity flow containing a soft sphere, and swimmers actuated via active stress., Comment: 12 pages, 5 figures (+ Supplemental Information: 26 pages, 10 figures, 3 tables, 6 movies, 1 dataset)

Lattice-Boltzmann methods are known for their simplicity, efficiency and ease of parallelization, usually relying on uniform Cartesian meshes with a strong bond between spatial and temporal discretization. This fact complicates the crucial issue of reducing the computational cost and the memory impact by automatically coarsening the grid where a fine mesh is unnecessary, still ensuring the overall quality of the numerical solution through error control. This work provides a possible answer to this interesting question, by connecting, for the first time, the field of lattice-Boltzmann Methods (LBM) to the adaptive multiresolution (MR) approach based on wavelets. To this end, we employ a MR multi-scale transform to adapt the mesh as the solution evolves in time according to its local regularity. The collision phase is not affected due to its inherent local nature and because we do not modify the speed of the sound, contrarily to most of the LBM/Adaptive Mesh Refinement (AMR) strategies proposed in the literature, thus preserving the original structure of any LBM scheme. Besides, an original use of the MR allows the scheme to resolve the proper physics by efficiently controlling the accuracy of the transport phase. We carefully test our method to conclude on its adaptability to a wide family of existing lattice Boltzmann schemes, treating both hyperbolic and parabolic systems of equations, thus being less problem-dependent than the AMR approaches, which have a hard time guaranteeing an effective control on the error. The ability of the method to yield a very efficient compression rate and thus a computational cost reduction for solutions involving localized structures with loss of regularity is also shown, while guaranteeing a precise control on the approximation error introduced by the spatial adaptation of the grid. The numerical strategy is implemented on a specific open-source platform called SAMURAI with a dedicated data-structure relying on set algebra.

We present a rigorous convergence result for smooth solutions to a singular semilinear hyperbolic approximation, called vector-BGK model, to the solutions to the incompressible Navier-Stokes equations in Sobolev spaces. Our proof deeply relies on the dissipative properties of the system and on the use of an energy which is provided by a symmetrizer, whose entries are weighted in a suitable way with respect to the singular perturbation parameter. This strategy allows us to perform uniform energy estimates and to prove the convergence by compactness.

We discretize the incompressible Navier–Stokes equations on a polytopal mesh by using mimetic reconstruction operators. The resulting method conserves discrete mass, momentum, and kinetic energy in the inviscid limit, and determines the vorticity such that the global vorticity is consistent with the boundary conditions. To do this we introduce a dual mesh and show how the dual mesh can be completed to a cell-complex. We present existing mimetic reconstruction operators in a new symmetric way applicable to arbitrary dimension, use these to interpolate between primal and dual mesh and derive properties of these operators. Finally, we test both 2- and 3-dimensional versions of the method on a variety of complicated meshes to show its wide applicability. We numerically test the convergence of the method and verify the derived conservation statements.

We present the spatial discretization of the nonlinear convective flux in the incompressible Navier–Stokes equations using local boundary-value problems. The finite-volume discretization of the incompressible Navier–Stokes equations over staggered grids requires the approximation of the cell-face velocity components. In the proposed method, the cell-face velocity components are computed from local boundary-value problems, with the pressure gradient and the gradient of the transverse flux as source terms. The approximation is expressed as the sum of the homogeneous part, corresponding to the convection–diffusion operator, and the inhomogeneous part, corresponding to the source terms. The homogeneous part of the approximation is shown to be a weighted average of the central and the upwind approximation and thus, is first-order convergent over coarse grids and second-order over finer grids. We show that the inhomogeneous part of the approximation effectively removes the numerical diffusion introduced by the homogeneous part. The inclusion of the source terms in the local boundary-value problem results in a more accurate approximation that shows uniform second-order convergence and does not introduce any spurious oscillations.

In this work we present a novel exact Riemann solver for an artificial Equation of State (EoS) based modification of the incompressible Euler equations. Differently from the well known artificial compressibility method, this new approach overcomes the lack of the pressure-velocity coupling by using a suitably designed artificial EoS. The modified set of equations fits into the framework of first order hyperbolic conservation laws. Accordingly, an exact Riemann solver is derived. This new solver has the advantage to avoid a wave pattern violation issue which can affect the exact Riemann problem solution obtained by the standard artificial compressibility approach. The new artificial EoS based Riemann solver can be used as a tool for the definition of the advective Godunov fluxes in a Finite Volume or a discontinuous Galerkin discretization of the incompressible Navier–Stokes (INS) equations. We assess and analyse the new exact Riemann solver on 1D Riemann problems. Its capability and effectiveness are then shown in the context of a high-order accurate discontinuous Galerkin discretization of the INS equations. In particular, we verify the convergence properties, both in space and time, considering two test cases in 2D, namely, the Kovasznay test case and a damped travelling waves test case. Finally, we perform an implicit large eddy simulation of the incompressible turbulent flow over periodic hills where the classic forcing term formulation is modified to deal with variable time steps. Solution comparisons with respect to numerical and experimental results from the literature are given.

The kinetic energy of a flow is proportional to the square of the L2(Ω) norm of the velocity. Given a sufficient regular velocity field and a velocity finite element space with polynomials of degree r, then the best approximation error in L2(Ω) is of order r+1. In this survey, the available finite element error analysis for the velocity error in L∞(0,T;L2(Ω)) is reviewed, where T is a final time. Since in practice the case of small viscosity coefficients or dominant convection is of particular interest, which may result in turbulent flows, robust error estimates are considered, i.e., estimates where the constant in the error bound does not depend on inverse powers of the viscosity coefficient. Methods for which robust estimates can be derived enable stable flow simulations for small viscosity coefficients on comparatively coarse grids, which is often the situation encountered in practice. To introduce stabilization techniques for the convection-dominated regime and tools used in the error analysis, evolutionary linear convection–diffusion equations are studied at the beginning. The main part of this survey considers robust finite element methods for the incompressible Navier–Stokes equations of order r−1, r, and r+1∕2 for the velocity error in L∞(0,T;L2(Ω)). All these methods are discussed in detail. In particular, a sketch of the proof for the error bound is given that explains the estimate of important terms which determine finally the order of convergence. Among them, there are methods for inf–sup stable pairs of finite element spaces as well as for pressure-stabilized discretizations. Numerical studies support the analytic results for several of these methods. In addition, methods are surveyed that behave in a robust way but for which only a non-robust error analysis is available. The conclusion of this survey is that the problem of whether or not there is a robust method with optimal convergence order for the kinetic energy is still open, The research of Bosco García-Archilla has been supported by Spanish MICIU under grant PGC2018-096265-B-I00 and MICI under grant PID2019-104141GB-I00. The Research of Julia Novo has been supported by Spanish MICI under grants PID2019-104141GB-I00 and VA169P20 (Junta de Castilla y Leon, ES) cofinanced by FEDER funds

A generalized approach to decompose the compressible Navier-Stokes equations into an equivalent set of coupled equations for flow and acoustics is introduced. As a significant extension to standard hydrodynamic/acoustic splitting methods, the approach provides the essential coupling terms, which account for the feedback from the acoustics to the flow. A unique simplified version of the split equation system with feedback is derived that conforms to the compressible Navier-Stokes equations in the subsonic flow regime, where the feedback reduces to one additional term in the flow momentum equation. Subsonic simulations are conducted for flow-acoustic feedback cases using a scale-resolving run-time coupled hierarchical Cartesian mesh solver, which operates with different explicit time step sizes for incompressible flow and acoustics. The first simulation case focuses on the tone of a generic flute. With the major flow-acoustic feedback term included, the simulation yields the tone characteristics in agreement with reference results from K\"uhnelt based on Lattice-Boltzmann simulation. On the contrary, the standard hybrid hydrodynamic/acoustic method with the feedback-term switched off lacks the proper tone. As the second simulation case, a thick plate in a duct is studied at various low Mach numbers around the Parker-beta-mode resonance. The simulations reveal the flow-acoustic feedback mechanism in very good agreement with experimental data of Welsh et al. Simulations and theoretical considerations reveal that the feedback term does not reduce the stable convective flow based time step size of the flow equations., Comment: Submitted to Journal of Computational Physics

This paper is devoted to the problem of solving the steady incompressible Navier–Stokes equations discretized by the Galerkin method on the spaces generated by the B-spline/NURBS basis functions, which is called isogeometric analysis. Two pressure-correction methods are presented for the solution of the incompressible flow in the benchmark backward facing step and also in Kaplan water turbine as conforming multipatch domains. The velocity and pressure under-relaxation is considered and the computational examples are compared with the coupled approach.

A novel reduced-order model (ROM) formulation for incompressible flows is presented with the key property that it exhibits non-linearly stability, independent of the mesh (of the full order model), the time step, the viscosity, and the number of modes. The two essential elements to non-linear stability are: (1) first discretise the full order model, and then project the discretised equations, and (2) use spatial and temporal discretisation schemes for the full order model that are globally energy-conserving (in the limit of vanishing viscosity). For this purpose, as full order model a staggered-grid finite volume method in conjunction with an implicit Runge-Kutta method is employed. In addition, a constrained singular value decomposition is employed which enforces global momentum conservation. The resulting ‘velocity-only’ ROM is thus globally conserving mass, momentum and kinetic energy. For non-homogeneous boundary conditions, a (one-time) Poisson equation is solved that accounts for the boundary contribution. The stability of the proposed ROM is demonstrated in several test cases. Furthermore, it is shown that explicit Runge-Kutta methods can be used as a practical alternative to implicit time integration at a slight loss in energy conservation.

SUMMARY An implicit fractional-step method for the numerical solution of the time-dependent incompressible Navier–Stokes equations in primitive variables is studied in this paper. The method, which is first-orderaccurate in the time step, is shown to converge to an exact solution of the equations. By adequately splitting the viscous term, it allows the enforcement of full Dirichlet boundary conditions on the velocity in all substeps of the scheme, unlike standard projection methods. The consideration of this method was actually motivated by the study of a well-known predictor–multicorrector algorithm, when this is applied to the incompressible Navier–Stokes equations. A new derivation of the algorithm in a general setting is provided, showing in what sense it can also be understood as a fractional-step method; this justifies, in particular, why the original boundary conditions of the problem can be enforced in this algorithm. Two different finite element interpolations are considered for the space discretization, and numerical results obtained with them for standard benchmark cases are presented. © 1998 John Wiley & Sons, Ltd.

In this work we present a mimetic spectral element discretization for the 2D incompressible Navier–Stokes equations that in the limit of vanishing dissipation exactly preserves mass, kinetic energy, enstrophy and total vorticity on unstructured triangular grids. The essential ingredients to achieve this are: (i) a velocity–vorticity formulation in rotational form, (ii) a sequence of function spaces capable of exactly satisfying the divergence free nature of the velocity field, and (iii) a conserving time integrator. Proofs for the exact discrete conservation properties are presented together with numerical test cases on highly irregular triangular grids.

In this work we consider the numerical solution of incompressible flows on two-dimensional manifolds. Whereas the compatibility demands of the velocity and the pressure spaces are known from the flat case one further has to deal with the approximation of a velocity field that lies only in the tangential space of the given geometry. Abandoning $H^1$-conformity allows us to construct finite elements which are -- due to an application of the Piola transformation -- exactly tangential. To reintroduce continuity (in a weak sense) we make use of (hybrid) discontinuous Galerkin techniques. To further improve this approach, $H(\operatorname{div}_{\Gamma})$-conforming finite elements can be used to obtain exactly divergence-free velocity solutions. We present several new finite element discretizations. On a number of numerical examples we examine and compare their qualitative properties and accuracy., Comment: 30 pages, 16 figures, 1 table

International audience; In this work we propose a novel Hybrid High-Order method for the incompressible Navier– Stokes equations based on a formulation of the convective term including Temam's device for stability. The proposed method has several advantageous features: it supports arbitrary approximation orders on general meshes including polyhedral elements and non-matching interfaces; it is inf-sup stable; it is locally conservative; it supports both the weak and strong enforcement of velocity boundary conditions; it is amenable to efficient computer implementations where a large subset of the unknowns is eliminated by solving local problems inside each element. Particular care is devoted to the design of the convective trilinear form, which mimicks at the discrete level the non-dissipation property of the continuous one. The possibility to add a convective stabilisation term is also contemplated, and a formulation covering various classical options is discussed. The proposed method is theoretically analysed, and an energy error estimate in $h^{k+1}$ (with $h$ denoting the meshsize) is proved under the usual data smallness assumption. A thorough numerical validation on two and three-dimensional test cases is provided both to confirm the theoretical convergence rates and to assess the method in more physical configurations (including, in particular, the well-known two-and three-dimensional lid-driven cavity problems).

This Accepted Manuscript will be available for reuse under a CC BY-NC-ND licence after 24 months of embargo period, In this paper we derive a posteriori error estimates for inf–sup stable mixed finite element approximations to the evolutionary Navier–Stokes equations. We reduce the problem of getting a posteriori error estimations of a non-linear evolutionary problem to that of getting a posteriori estimations of a linear steady problem. The main idea is based on the fact that the solutions of some Newton-type linearized Navier–Stokes equations around the plain Galerkin approximations approach the solution of the original Navier–Stokes equations with a bigger rate of convergence than the plain Galerkin method. As a consequence, the difference between the Galerkin approximations and the solution of the linearized problem is an estimator of the Galerkin error. Moreover, since the Galerkin approximation of the evolutionary Navier–Stokes equations is also the Galerkin approximation of the linearized equations, any a posteriori estimator of the error in the Newton-type linearized Navier–Stokes equations is also an a posteriori error estimator of the full Navier–Stokes equations, Research supported by Spanish MEC under grant MTM2016-78995-P (AEI/FEDER, UE) and VA024P17 (Junta de Castilla y Leon, ES)

A generic particle–mesh method using a hybridized discontinuous Galerkin (HDG) framework is presented and validated for the solution of the incompressible Navier–Stokes equations. Building upon particle-in-cell concepts, the method is formulated in terms of an operator splitting technique in which Lagrangian particles are used to discretize an advection operator, and an Eulerian mesh-based HDG method is employed for the constitutive modeling to account for the inter-particle interactions. Key to the method is the variational framework provided by the HDG method. This allows to formulate the projections between the Lagrangian particle space and the Eulerian finite element space in terms of local (i.e. cellwise) l 2 -projections efficiently. Furthermore, exploiting the HDG framework for solving the constitutive equations results in velocity fields which excellently approach the incompressibility constraint in a local sense. By advecting the particles through these velocity fields, the particle distribution remains uniform over time, obviating the need for additional quality control. The presented methodology allows for a straightforward extension to arbitrary-order spatial accuracy on general meshes. A range of numerical examples shows that optimal convergence rates are obtained in space and, given the particular time stepping strategy, second-order accuracy is obtained in time. The model capabilities are further demonstrated by presenting results for the flow over a backward facing step and for the flow around a cylinder.

A strategy for choosing the forcing terms in inexact Newton iterations is presented. The final goal is to obtain fast steady state marching strategies for the solution of PDEs. The new approach is analyzed and tested in the context of inexact Newton methods but is also well suited to be applied in the pseudo-transient continuation framework. To validate the strategy and assess its gains in terms of computational costs we seek approximate solutions of the incompressible Navier-Stokes equations at high-Reynolds numbers. In particular we consider the well known 2D lid-driven cavity flow and backward-facing step problem focusing on the efficiency of the time marching strategy. Residual history and computation time are monitored and compared with many fixed and adaptive forcing term choices of reference.

The aim of this paper is to prove the strong convergence of the solutions to a vector-BGK model under the diffusive scaling to the incompressible Navier-Stokes equations on the two-dimensional torus. This result holds in any interval of time [ 0 , T ] , with T > 0 . We also provide the global in time uniform boundedness of the solutions to the approximating system. Our argument is based on the use of local in time H s -estimates for the model, established in a previous work, combined with the L 2 -relative entropy estimate and the interpolation properties of the Sobolev spaces.

Grad-div stabilization is a classical remedy in conforming mixed finite element methods for incompressible flow problems, for mitigating velocity errors that are sometimes called poor mass conservation. Such errors arise due to the relaxation of the divergence constraint in classical mixed methods, and are excited whenever the spatial discretization has to deal with comparably large and complicated pressures. In this contribution, an analogue of grad-div stabilization for Discontinuous Galerkin methods is studied. Here, the key is the penalization of the jumps of the normal velocities over facets of the triangulation, which controls the measure-valued part of the distributional divergence of the discrete velocity solution. Our contribution is twofold: first, we characterize the limit for arbitrarily large penalization parameters, which shows that the stabilized nonconforming Discontinuous Galerkin methods remain robust and accurate in this limit; second, we extend these ideas to the case of non-simplicial meshes; here, broken grad-div stabilization must be used in addition to the normal velocity jump penalization, in order to get the desired pressure robustness effect. The analysis is performed for the Stokes equations, and more complex flows and Crouzeix-Raviart elements are considered in numerical examples that also show the relevance of the theory in practical settings., 24 pages, 3 figures, 5 tables

International audience; This paper presents the common mathematical features which are leading to convergence properties for a family of numerical schemes applied to the discretisation of the steady and transient incompressible Navier-Stokes equations with homogeneous Dirichlet’s boundary conditions. Thisfamily includes the Taylor-Hood scheme, the MAC scheme, the Crouzeix-Raviart scheme generalized into the Hybrid Mixed Mimetic scheme, which can be combined with a variety of discretisations for the nonlinear convection term, each of them being more efficient than the others in particular situa-tions. We provide tools for analyzing all the combined methods, and proving their convergence to a weak solution of the problem.

International audience; In this work a new higher-order (>2) accurate finite volume method for the resolution of the incompressible Navier-Stokes equations on unstructured grids is presented. The formulation is based on the use of Moving Least Squares (MLS) approximants. Third and fourth order accurate discretizations of the convective and viscous fluxes are obtained on collocated meshes. In addition, MLS is employed to design a new Momentum Interpolation Method that allows interpolations better than linear on any kind of mesh. The accuracy and performance of the proposed method is demonstrated by solving different steady and unsteady benchmark problems on unstructured grids.

The application of pseudo-symplectic Runge–Kutta methods to the incompressible Navier–Stokes equations is discussed in this work. In contrast to fully energy-conserving, implicit methods, these are explicit schemes of order p that preserve kinetic energy to order q, with q>p. Use of explicit methods with improved energy-conservation properties is appealing for convection-dominated problems, especially in case of direct and large-eddy simulation of turbulent flows. A number of pseudo-symplectic methods are constructed for application to the incompressible Navier–Stokes equations and compared in terms of accuracy and efficiency by means of numerical simulations.

We introduce a space–time discontinuous Galerkin (DG) finite element method for the incompressible Navier–Stokes equations. Our formulation can be made arbitrarily high order accurate in both space and time and can be directly applied to deforming domains. Different stabilizing approaches are discussed which ensure stability of the method. A numerical study is performed to compare the effect of the stabilizing approaches, to show the method’s robustness on deforming domains and to investigate the behavior of the convergence rates of the solution. Recently we introduced a space–time hybridizable DG (HDG) method for incompressible flows [S. Rhebergen, B. Cockburn, A space–time hybridizable discontinuous Galerkin method for incompressible flows on deforming domains, J. Comput. Phys. 231 (2012) 4185–4204]. We will compare numerical results of the space–time DG and space–time HDG methods. This constitutes the first comparison between DG and HDG methods.

The objective of the present work is to extend our FDS-based third-order upwind compact schemes by Shah et al. (2009) [8] to numerical solutions of the unsteady incompressible Navier–Stokes equations in curvilinear coordinates, which will save much computing time and memory allocation by clustering grids in regions of high velocity gradients. The dual-time stepping approach is used for obtaining a divergence-free flow field at each physical time step. We have focused on addressing the crucial issue of implementing upwind compact schemes for the convective terms and a central compact scheme for the viscous terms on curvilinear structured grids. The method is evaluated in solving several two-dimensional unsteady benchmark flow problems.

The Artificial Compressibility Method (ACM) for the incompressible Navier-Stokes equations is (link-wise) reformulated (referred to as LW-ACM) by a finite set of discrete directions (links) on a regular Cartesian mesh, in analogy with the Lattice Boltzmann Method (LBM). The main advantage is the possibility of exploiting well established technologies originally developed for LBM and classical computational fluid dynamics, with special emphasis on finite differences (at least in the present paper), at the cost of minor changes. For instance, wall boundaries not aligned with the background Cartesian mesh can be taken into account by tracing the intersections of each link with the wall (analogously to LBM technology). LW-ACM requires no high-order moments beyond hydrodynamics (often referred to as ghost moments) and no kinetic expansion. Like finite difference schemes, only standard Taylor expansion is needed for analyzing consistency. Preliminary efforts towards optimal implementations have shown that LW-ACM is capable of similar computational speed as optimized (BGK-) LBM. In addition, the memory demand is significantly smaller than (BGK-) LBM. Importantly, with an efficient implementation, this algorithm may be one of the few which is compute-bound and not memory-bound. Two- and three-dimensional benchmarks are investigated, and an extensive comparative study between the present approach and state of the art methods from the literature is carried out. Numerical evidences suggest that LW-ACM represents an excellent alternative in terms of simplicity, stability and accuracy., Comment: 62 pages, 20 figures

In this paper we prove some properties of the maximal solution of Navier–Stokes equations. If the maximum time T ν ∗ is finite, we establish that the growth of ‖ u ( t ) ‖ H ˙ s is at least of the order of ( T ν ∗ − t ) − s / 3 (see Eq. (1.4)), also we give some new blow-up results. Specific properties and standard techniques are used.

In the present work, a highly efficient incompressible flow solver with a semi-implicit time advancement on a fully staggered grid using a high-order compact difference scheme is developed firstly in the framework of approximate factorization. The fourth-order compact difference scheme is adopted for approximations of derivatives and interpolations in the incompressible Navier–Stokes equations. The pressure Poisson equation is efficiently solved by the fast Fourier transform (FFT). The framework of approximate factorization significantly simplifies the implementation of the semi-implicit time advancing with a high-order compact scheme. Benchmark tests demonstrate the high accuracy of the proposed numerical method. Secondly, by applying the proposed numerical method, we compute turbulent channel flows at low and moderate Reynolds numbers by direct numerical simulation (DNS) and large eddy simulation (LES). It is found that the predictions of turbulence statistics and especially energy spectra can be obviously improved by adopting the high-order scheme rather than the traditional second-order central difference scheme.

Finite element methods for problems given in complex domains are often based on tetrahedral meshes. This paper demonstrates that the so-called rational Large Eddy Simulation model and a projection-based Variational Multiscale method can be extended in a straightforward way to tetrahedral meshes. Numerical studies are performed with an inf–sup stable second order pair of finite elements with discontinuous pressure approximation.

The artificial compressibility method for the incompressible Navier-Stokes equations is revived as a high order accurate numerical method (fourth order in space and second order in time). Similar to the lattice Boltzmann method, the mesh spacing is linked to the Mach number. An accuracy higher than that of the lattice Boltzmann method is achieved by exploiting the asymptotic behavior of the solution of the artificial compressibility equations for small Mach numbers and the simple lattice structure. An easy method for accelerating the decay of acoustic waves, which deteriorate the quality of the numerical solution, and a simple cure for the checkerboard instability are proposed. The high performance of the scheme is demonstrated not only for the periodic boundary condition but also for the Dirichlet-type boundary condition.

The present paper proves that all limit points of sequences of renormalized solutions of the Boltzmann equation in the limit of small, asymptotically equivalent Mach and Knudsen numbers are governed by Leray solutions of the Navier-Stokes equations. This convergence result holds for hard cutoff potentials in the sense of H. Grad, and therefore completes earlier results by the same authors [Invent. Math. 155, 81-161(2004)] for Maxwell molecules., 56 pages, LaTeX, a few typos have been corrected, a few remarks added, one uncited reference removed

The electro-diffusion model, which arises in electrohydrodynamics, is a coupling between the Nernst-Planck-Poisson system and the incompressible Navier-Stokes equations. For the generally smooth doping profile, the quasineutral limit (zero-Debye-length limit) is justified rigorously in Sobolev norm uniformly in time. The proof is based on the elaborate energy analysis and the key point is to establish the uniform estimates with respect to the scaled Debye length., 20 pages

The quasineutral limit (zero-Debye-length limit) of viscous quantum hydrodynamic model for semiconductors is studied in this paper. By introducing new modulated energy functional and using refined energy analysis, it is shown that, for well-prepared initial data, the smooth solution of viscous quantum hydrodynamic model converges to the strong solution of incompressible Navier–Stokes equations as the Debye length goes to zero.

We present a well-posed stochastic Galerkin formulation of the incompressible Navier–Stokes equations with uncertainty in model parameters or the initial and boundary conditions. The stochastic Galerkin method involves representation of the solution through generalized polynomial chaos expansion and projection of the governing equations onto stochastic basis functions, resulting in an extended system of equations. A relatively low-order generalized polynomial chaos expansion is sufficient to capture the stochastic solution for the problem considered. We derive boundary conditions for the continuous form of the stochastic Galerkin formulation of the velocity and pressure equations. The resulting problem formulation leads to an energy estimate for the divergence. With suitable boundary data on the pressure and velocity, the energy estimate implies zero divergence of the velocity field. Based on the analysis of the continuous equations, we present a semi-discretized system where the spatial derivatives are approximated using finite difference operators with a summation-by-parts property. With a suitable choice of dissipative boundary conditions imposed weakly through penalty terms, the semi-discrete scheme is shown to be stable. Numerical experiments in the laminar flow regime corroborate the theoretical results and we obtain high-order accurate results for the solution variables and the velocity divergence converges to zero as the mesh is refined. Funding agencies: SUPRI-B at Stanford University; U.S. Department of Energy Office of Science, Office of Advanced Scientific Computing Research [DE-SC0006402]; Uni Research, Norway

The possibility to use a Lagrangian frame to solve problems with large time-steps was successfully explored previously by the authors for the solution of homogeneous incompressible fluids and also for solving multi-fluid problems (Idelsohn etal. 2012; 2014; 2013). The strategy used by the authors was named Particle Finite Element Method second generation (PFEM-2).The objective of this paper is to demonstrate in which circumstances the use of a Lagrangian frame with particles is more accurate than a classical Eulerian finite element method, and when large time-steps and/or coarse meshes may be used. Fil: Idelsohn, Sergio Rodolfo. Institució Catalana de Recerca i Estudis Avancats; España. International Center for Numerical Methods in Engineering; España. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Santa Fe. Centro de Investigaciones en Métodos Computacionales. Universidad Nacional del Litoral. Centro de Investigaciones en Métodos Computacionales; Argentina Fil: Oñate, Eugenio. International Center for Numerical Methods in Engineering; España Fil: Nigro, Norberto Marcelo. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Santa Fe. Centro de Investigaciones en Métodos Computacionales. Universidad Nacional del Litoral. Centro de Investigaciones en Métodos Computacionales; Argentina Fil: Becker, Pablo Javier. International Center for Numerical Methods in Engineering; España Fil: Gimenez, Juan Marcelo. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Santa Fe. Centro de Investigaciones en Métodos Computacionales. Universidad Nacional del Litoral. Centro de Investigaciones en Métodos Computacionales; Argentina

A fourth-order accurate solution method for the three-dimensional Helmholtz equations is described that is based on a compact finite-difference stencil for the Laplace operator. Similar discretization methods for the Poisson equation have been presented by various researchers for Dirichlet boundary conditions. Here, the complicated issue of imposing Neumann boundary conditions is described in detail. The method is then applied to model Helmholtz problems to verify the accuracy of the discretization method. The implementation of the solution method is also described. The Helmholtz solver is used as the basis for a fourth-order accurate solver for the incompressible Navier-Stokes equations. Numerical results obtained with this Navier-Stokes solver for the temporal evolution of a three-dimensional instability in a counter-rotating vortex pair are discussed. The time-accurate Navier-Stokes simulations show the resolving properties of the developed discretization method and the correct prediction of the initial growth rate of the three-dimensional instability in the vortex pair.

Weighted essentially non-oscillatory (WENO) finite difference schemes have been recommended in a competitive study of discretizations for scalar evolutionary convection-diffusion equations [20]. This paper explores the applicability of these schemes for the simulation of incompressible flows. To this end, WENO schemes are used in several non-incremental and incremental projection methods for the incompressible Navier-Stokes equations. Velocity and pressure are discretized on the same grid. A pressure stabilization Petrov-Galerkin (PSPG) type of stabilization is introduced in the incremental schemes to account for the violation of the discrete inf-sup condition. Algorithmic aspects of the proposed schemes are discussed. The schemes are studied on several examples with different features. It is shown that the WENO finite difference idea can be transferred to the simulation of incompressible flows. Some shortcomings of the methods, which are due to the splitting in projection schemes, become also obvious.

In this paper we present a model for incompressible chemically reacting flows where reactants enter the domain, react, and then leave the domain. In the interior of the domain we have the Navier–Stokes equations for the fluid flow coupled with reaction diffusion equations for the chemistry and temperature. On the boundary we use mixed Neumann and inhomogeneous Robin boundary conditions for the chemistry and temperature equations which fix the amount of chemicals and heat flowing into the system. For three dimensional domains we then show the existence of weak solutions which are physically reasonable, i.e., which satisfy certain maximum principle properties. These solutions are sufficient to provide the basis for a dynamical systems study of these reacting flows. We do not show that the solutions are unique; however, we are able to show the existence of global attractors for the system.

Harlow and Welch (Phys. Fluids 8:2182-2189, 1965) introduced a discretization method for the incompressible Navier-Stokes equations conserving the secondary quantities kinetic energy and vorticity, besides the primary quantities mass and momentum. This method was extended to fourth order accuracy by several researchers [25], [14] and [21]. In this paper we propose a new consistent boundary treatment for this method, which is such that continuous integration-by-parts identities (including boundary contributions) are mimicked in a discrete sense. In this way kinetic energy is exactly conserved even in case of non-zero tangential boundary conditions. We show that requiring energy conservation at the boundary conflicts with order of accuracy conditions, and that the global accuracy of the fourth order method is limited to second order in the presence of boundaries. We indicate how non-uniform grids can be employed to obtain full fourth order accuracy.

International audience; In this paper we study a three-dimensional fluid-structure interaction problem. The motion of the fluid is modeled by the Navier-Stokes equations and we consider for the elastic structure a finite dimensional approximation of the equation of linear elasticity. The time variation of the fluid domain is not known a priori, so we deal with a free boundary value problem. Our main result yields the local in time existence and uniqueness of strong solutions for this system.

International audience; In this work we introduce a new class of lowest order methods for diffusive problems on general meshes with only one unknown per element. The underlying idea is to construct an incomplete piecewise affine polynomial space with sufficient approximation properties starting from values at cell centers. To do so we borrow ideas from multi-point finite volume methods, although we use them in a rather different context. The incomplete polynomial space replaces classical complete polynomial spaces in discrete formulations inspired by discontinuous Galerkin methods. Two problems are studied in this work: a heterogeneous anisotropic diffusion problem, which is used to lay the pillars of the method, and the incompressible Navier-Stokes equations, which provide a more realistic application. An exhaustive theoretical study as well as a set of numerical examples featuring different difficulties are provided.

In this work we present a pressure-correction scheme for the incompressible Navier-Stokes equations combining a discontinuous Galerkin approximation for the velocity and a standard continuous Galerkin approximation for the pressure. The main interest of pressure-correction algorithms is the reduced computational cost compared to monolithic strategies. In this work we show how a proper discretization of the decoupled momentum equation can render this method suitable to simulate high Reynolds regimes. The proposed spatial velocity-pressure approximation is LBB stable for equal polynomial orders and it allows adaptive p-refinement for velocity and global p-refinement for pressure. The method is validated against a large set of classical two- and three-dimensional test cases covering a wide range of Reynolds numbers, in which it proves effective both in terms of accuracy and computational cost.

From the literature, it is known that the Least-Squares Spectral Element Method (LSSEM) for the stationary Stokes equations performs poorly with respect to mass conservation but compensates this lack by a superior conservation of momentum. Furthermore, it is known that the Least-Squares Spectral Collocation Method (LSSCM) leads to superior conservation of mass and momentum for the stationary Stokes equations. In the present paper, we consider mass and momentum conservation of the LSSCM for time-dependent Stokes and Navier–Stokes equations. We observe that the LSSCM leads to improved conservation of mass (and momentum) for these problems. Furthermore, the LSSCM leads to the well-known time-dependent profiles for the velocity and the pressure profiles. To obtain these results, we use only a few elements, each with high polynomial degree, avoid normal equations for solving the overdetermined linear systems of equations and introduce the Clenshaw–Curtis quadrature rule for imposing the average pressure to be zero. Furthermore, we combined the transformation of Gordon and Hall (transfinite mapping) with the least-squares spectral collocation scheme to discretize the internal flow problems.

This paper presents a computational vascular fluid-structure interaction (FSI) methodology and its application to patient-specific aneurysm models of the middle cerebral artery bifurcation. A fully coupled fluid-structural simulation approach is reviewed, and main aspects of mesh generation in support of patient-specific vascular FSI analyses are presented. Quantities of hemodynamic interest such as wall shear stress and wall tension are studied to examine the relevance of FSI modeling as compared to the rigid arterial wall assumption. We demonstrate the importance of including the flexible wall modeling in vascular blood flow simulations by performing a comparison study that involves four patient-specific models of cerebral aneurysms varying in shape and size.

We present a new coupled discretization approach for species transport in an incompressible fluid. The Navier-Stokes equations for the flow are discretized by the divergence-free Scott-Vogelius element on barycentrically refined meshes guaranteeing LBB stability. The convection-diffusion equation for species transport is discretized by the Voronoi finite volume method. In accordance to the continuous setting, due to the exact integration of the normal component of the flow through the Voronoi surfaces, the species concentration fulfills discrete global and local maximum principles. Besides of the the numerical scheme itself, we present important aspects of its implementation. Further, for the case of homogeneous Dirichlet boundary conditions, we give a convergence proof for the coupled scheme. We report results of the application of the scheme to the interpretation of limiting current measurements in an electrochemical flow cell with cylindrical shape.