Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userStability analysis for a new model of multispecies convection-diffusion-reaction in poroelastic tissue
84
0
0.0
(
0
)
Added by
Ricardo Ruiz Baier I
Publication date
2020
fields
Informatics Engineering
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
We perform the linear stability analysis for a new model for poromechanical processes with inertia (formulated in mixed form using the solid deformation, fluid pressure, and total pressure) interacting with diffusing and reacting solutes convected in the medium. We find parameter regions that lead to spatio-temporal instabilities of the coupled system. The mutual dependences between deformation and diffusive patterns are of substantial relevance in the study of morphoelastic changes in biomaterials. We provide a set of computational examples in 2D and 3D (related to brain mechanobiology) that can be used to form a better understanding on how, and up to which extent, the deformations of the porous structure dictate the generation and suppression of spatial patterning dynamics, also related to the onset of mechano-chemical waves.
rate research
Read More
We analyse a PDE system modelling poromechanical processes (formulated in mixed form using the solid deformation, fluid pressure, and total pressure) interacting with diffusing and reacting solutes in the medium. We investigate the well-posedness of the nonlinear set of equations using fixed-point theory, Fredholms alternative, a priori estimates, and compactness arguments. We also propose a mixed finite element method and rigorously demonstrate the stability of the scheme. Error estimates are derived in suitable norms, and numerical experiments are conducted to illustrate the mechano-chemical coupling and to verify the theoretical rates of convergence.
In this article, using the weighted discrete least-squares, we propose a patch reconstruction finite element space with only one degree of freedom per element. As the approximation space, it is applied to the discontinuous Galerkin methods with the upwind scheme for the steady-state convection-diffusion-reaction problems over polytopic meshes. The optimal error estimates are provided in both diffusion-dominated and convection-dominated regimes. Furthermore, several numerical experiments are presented to verify the theoretical error estimates, and to well approximate boundary layers and/or internal layers.
A study is presented on the convergence of the computation of coupled advection-diffusion-reaction equations. In the computation, the equations with different coefficients and even types are assigned in two subdomains, and Schwarz iteration is made between the equations when marching from a time level to the next one. The analysis starts with the linear systems resulting from the full discretization of the equations by explicit schemes. Conditions for convergence are derived, and its speedup and the effects of difference in the equations are discussed. Then, it proceeds to an implicit scheme, and a recursive expression for convergence speed is derived. An optimal interface condition for the Schwarz iteration is obtained, and it leads to perfect convergence, that is, convergence within two times of iteration. Furthermore, the methods and analyses are extended to the coupling of the viscous Burgers equations. Numerical experiments indicate that the conclusions, such as the perfect convergence, drawn in the linear situations may remain in the Burgers equations computation.
We present and analyze a new iterative solver for implicit discretizations of a simplified Boltzmann-Poisson system. The algorithm builds on recent work that incorporated a sweeping algorithm for the Vlasov-Poisson equations as part of nested inner-outer iterative solvers for the Boltzmann-Poisson equations. The new method eliminates the need for nesting and requires only one transport sweep per iteration. It arises as a new fixed-point formulation of the discretized system which we prove to be contractive for a given electric potential. We also derive an accelerator to improve the convergence rate for systems in the drift-diffusion regime. We numerically compare the efficiency of the new solver, with and without acceleration, with a recently developed nested iterative solver.
We present a Petrov-Gelerkin (PG) method for a class of nonlocal convection-dominated diffusion problems. There are two main ingredients in our approach. First, we define the norm on the test space as induced by the trial space norm, i.e., the optimal test norm, so that the inf-sup condition can be satisfied uniformly independent of the problem. We show the well-posedness of a class of nonlocal convection-dominated diffusion problems under the optimal test norm with general assumptions on the nonlocal diffusion and convection kernels. Second, following the framework of Cohen et al.~(2012), we embed the original nonlocal convection-dominated diffusion problem into a larger mixed problem so as to choose an enriched test space as a stabilization of the numerical algorithm. In the numerical experiments, we use an approximate optimal test norm which can be efficiently implemented in 1d, and study its performance against the energy norm on the test space. We conduct convergence studies for the nonlocal problem using uniform $h$- and $p$-refinements, and adaptive $h$-refinements on both smooth manufactured solutions and solutions with sharp gradient in a transition layer. In addition, we confirm that the PG method is asymptotically compatible.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
