This work is focused on the analysis of non-linear flows of slightly compressible fluids in porous media not adequately described by Darcy's law. We study a class of generalized nonlinear momentum equations which covers all three well-known Forchheimer equations, the so-called two-term, power, and three-term laws. The non-linear Forchheimer equation is inverted to a non-linear Darcy equation with implicit permeability tensor depending on the pressure gradient. This results in a degenerate parabolic equation for the pressure. Two classes of boundary conditions are considered, given pressure and given total flux. In both cases they are allowed to be unbounded in time. The uniqueness, Lyapunov and asymptotic stabilities, and other long-time dynamical features of the corresponding initial boundary value problems are analyzed. The results obtained in this paper have clear hydrodynamic interpretations and can be used for quantitative evaluation of engineering parameters. Some numerical simulations are also included.
The application of the new numerical approach for elastodynamics problems developed in our previous paper and based on the new solution strategy and the new time-integration methods is considered for 1D and 2D axisymmetric impact problems. It is not easy to solve these problems accurately because the exact solutions of the corresponding semi-discrete elastodynamics problems contain a large number of spurious high-frequency oscillations. We use the 1D impact problem for the calibration of a new analytical expression describing the minimum amount of numerical dissipation necessary for the new time-integration method used for filtering spurious oscillations. Then, we show that the new numerical approach for elastodynamics along with the new expression for numerical dissipation for the first time yield accurate and non-oscillatory solutions of the considered impact problems. The comparison of effectiveness of linear and quadratic elements as well as rectangular and
In this paper we study a monolithic Newton-Krylov solver with exact Jacobian for the solution of incompressible FSI problems. A main focus of this work is on the use of geometric multigrid preconditioners with modified Richardson smoothers preconditioned by an additive Schwarz algorithm. The definition of the subdomains in the Schwarz smoother is driven by the natural splitting between fluid and solid. The monolithic approach guarantees the automatic satisfaction of the stress balance and the kinematic conditions across the fluid-solid interface. The enforcement of the incompressibility conditions both for the fluid and for the solid parts is taken care of by using inf-sup stable finite element pairs without stabilization terms. A suitable Arbitrary Lagrangian Eulerian (ALE) operator is chosen in order to avoid mesh entanglement while solving for large displacements of the moving fluid domain. Numerical results of two and three-dimensional benchmark tests with Newtonian fluids and nonlinear hyperelastic solids show a robust performance of our fully incompressible solver especially for the more challenging direct-to-steady-state problems.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.