A second-order backward difference time-stepping scheme for penalized Navier-Stokes equations modeling filtration through porous media

2015 ◽  
Vol 32 (2) ◽  
pp. 681-705 ◽  
Author(s):  
S. S. Ravindran
2019 ◽  
Vol 19 (3) ◽  
pp. 681-701 ◽  
Author(s):  
Max Gunzburger ◽  
Nan Jiang ◽  
Zhu Wang

AbstractWe consider settings for which one needs to perform multiple flow simulations based on the Navier–Stokes equations, each having different initial condition data, boundary condition data, forcing functions, and/or coefficients such as the viscosity. For such settings, we propose a second-order time accurate ensemble-based method that to simulate the whole set of solutions, requires, at each time step, the solution of only a single linear system with multiple right-hand-side vectors. Rigorous analyses are given proving the conditional stability and establishing error estimates for the proposed algorithm. Numerical experiments are provided that illustrate the analyses.


2021 ◽  
Vol 0 (0) ◽  
Author(s):  
Riccardo Milani ◽  
Jérôme Bonelle ◽  
Alexandre Ern

Abstract We investigate artificial compressibility (AC) techniques for the time discretization of the incompressible Navier–Stokes equations. The space discretization is based on a lowest-order face-based scheme supporting polytopal meshes, namely discrete velocities are attached to the mesh faces and cells, whereas discrete pressures are attached to the mesh cells. This face-based scheme can be embedded into the framework of hybrid mixed mimetic schemes and gradient schemes, and has close links to the lowest-order version of hybrid high-order methods devised for the steady incompressible Navier–Stokes equations. The AC time-stepping uncouples at each time step the velocity update from the pressure update. The performances of this approach are compared against those of the more traditional monolithic approach which maintains the velocity-pressure coupling at each time step. We consider both first-order and second-order time schemes and either an implicit or an explicit treatment of the nonlinear convection term. We investigate numerically the CFL stability restriction resulting from an explicit treatment, both on Cartesian and polytopal meshes. Finally, numerical tests on large 3D polytopal meshes highlight the efficiency of the AC approach and the benefits of using second-order schemes whenever accurate discrete solutions are to be attained.


Analysis ◽  
2020 ◽  
Vol 40 (3) ◽  
pp. 127-150
Author(s):  
Tania Biswas ◽  
Sheetal Dharmatti ◽  
Manil T. Mohan

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


Author(s):  
Hisham Elsafti ◽  
Hocine Oumeraci

In this study, the fully-coupled and fully-dynamic Biot governing equations in the open-source geotechFoam solver are extended to account for pore fluid viscous stresses. Additionally, turbulent pore fluid flow in deformable porous media is modeled by means of the conventional eddy viscosity concept without the need to resolve all turbulence scales. A new approach is presented to account for porous media resistance to flow (solid-to-fluid coupling) by means of an effective viscosity, which accounts for tortuosity, grain shape and local turbulences induced by flow through porous media. The new model is compared to an implemented extended Darcy-Forchheimer model in the Navier-Stokes equations, which accounts for laminar, transitional, turbulent and transient flow regimes. Further, to account for skeleton deformation, the porosity and other model parameters are updated with regard to strain of geomaterials. The presented model is calibrated by means of available results of physical experiments of unidirectional and oscillatory flows.


Sign in / Sign up

Export Citation Format

Share Document