A Systematic Approach for Correcting Nonlinear Instabilities. The Lax-Wendroff Scheme for Scalar Conservation Laws.
A three dimensional finite element method for biological active soft tissue formulation in cylindrical polar coordinates
A hyperelastic constitutive law, for use in anatomically accurate finite element models of living structures, is suggested for the passive and the active mechanical properties of incompressible biological tissues. This law considers the passive and active states as a same hyperelastic continuum medium, and uses an activation function in order to describe the whole contraction phase. The variational and the FE formulations are also presented, and the FE code has been validated and applied to describe...
A three dimensional finite element method for biological active soft tissue Formulation in cylindrical polar coordinates
A hyperelastic constitutive law, for use in anatomically accurate finite element models of living structures, is suggested for the passive and the active mechanical properties of incompressible biological tissues. This law considers the passive and active states as a same hyperelastic continuum medium, and uses an activation function in order to describe the whole contraction phase. The variational and the FE formulations are also presented, and the FE code has been validated and applied to describe...
A two parameter iterative method for solving algebraic systems of domain decomposition type
An iterative procedure containing two parameters for linear algebraic systems originating from the domain decomposition technique is proposed. The optimization of the parameters is investigated. A numeric example is given as an illustration.
A two-fluid hyperbolic model in a porous medium
The paper is devoted to the computation of two-phase flows in a porous medium when applying the two-fluid approach. The basic formulation is presented first, together with the main properties of the model. A few basic analytic solutions are then provided, some of them corresponding to solutions of the one-dimensional Riemann problem. Three distinct Finite-Volume schemes are then introduced. The first two schemes, which rely on the Rusanov scheme, are shown to give wrong approximations in some...
A uniformly controllable and implicit scheme for the 1-D wave equation
This paper studies the exact controllability of a finite dimensional system obtained by discretizing in space and time the linear 1-D wave system with a boundary control at one extreme. It is known that usual schemes obtained with finite difference or finite element methods are not uniformly controllable with respect to the discretization parameters and . We introduce an implicit finite difference scheme which differs from the usual centered one by additional terms of order and . Using a discrete...
A uniformly controllable and implicit scheme for the 1-D wave equation
This paper studies the exact controllability of a finite dimensional system obtained by discretizing in space and time the linear 1-D wave system with a boundary control at one extreme. It is known that usual schemes obtained with finite difference or finite element methods are not uniformly controllable with respect to the discretization parameters h and Δt. We introduce an implicit finite difference scheme which differs from the usual centered one by additional terms of order h2 and Δt2. Using...
A wavelet based space-time adaptive numerical method for partial differential equations
A wavelet Galerkin method applied to partial differential equations with variable coefficients.
Acceleration of a fixed point algorithm for fluid-structure interaction using transpiration conditions
In this work, we address the numerical solution of fluid-structure interaction problems. This issue is particularly difficulty to tackle when the fluid and the solid densities are of the same order, for instance as it happens in hemodynamic applications, since fully implicit coupling schemes are required to ensure stability of the resulting method. Thus, at each time step, we have to solve a highly non-linear coupled system, since the fluid domain depends on the unknown displacement of the structure....
Acceleration of a fixed point algorithm for fluid-structure interaction using transpiration conditions
In this work, we address the numerical solution of fluid-structure interaction problems. This issue is particularly difficulty to tackle when the fluid and the solid densities are of the same order, for instance as it happens in hemodynamic applications, since fully implicit coupling schemes are required to ensure stability of the resulting method. Thus, at each time step, we have to solve a highly non-linear coupled system, since the fluid domain depends on the unknown displacement of...
Accurate numerical discretizations of non-conservative hyperbolic systems
We present an alternative framework for designing efficient numerical schemes for non-conservative hyperbolic systems. This approach is based on the design of entropy conservative discretizations and suitable numerical diffusion operators that mimic the effect of underlying viscous mechanisms. This approach is illustrated by considering two model non-conservative systems: Lagrangian gas dynamics in non-conservative form and a form of isothermal Euler equations. Numerical experiments demonstrating...
Accurate numerical discretizations of non-conservative hyperbolic systems
We present an alternative framework for designing efficient numerical schemes for non-conservative hyperbolic systems. This approach is based on the design of entropy conservative discretizations and suitable numerical diffusion operators that mimic the effect of underlying viscous mechanisms. This approach is illustrated by considering two model non-conservative systems: Lagrangian gas dynamics in non-conservative form and a form of isothermal Euler equations. Numerical experiments demonstrating...
Adapting meshes and time-steps for phase change problems
We address the numerical approximation of the two-phase Stefan problem and discuss an adaptive finite element method based on rigorous a posteriori error estimation and refinement/coarsening. We also investigate how to restrict coarsening for the resulting method to be stable and convergent. We review implementation issues associated with bisection and conclude with simulations of a persistent corner singularity, for which adaptivity is an essential tool.
Adaptive fast marching and level set methods for propagating interfaces.
Adaptive finite element method for shape optimization
We examine shape optimization problems in the context of inexact sequential quadratic programming. Inexactness is a consequence of using adaptive finite element methods (AFEM) to approximate the state and adjoint equations (via the dual weighted residual method), update the boundary, and compute the geometric functional. We present a novel algorithm that equidistributes the errors due to shape optimization and discretization, thereby leading to coarse resolution in the early stages and fine resolution...
Adaptive finite element relaxation schemes for hyperbolic conservation laws
We propose and study semidiscrete and fully discrete finite element schemes based on appropriate relaxation models for systems of Hyperbolic Conservation Laws. These schemes are using piecewise polynomials of arbitrary degree and their consistency error is of high order. The methods are combined with an adaptive strategy that yields fine mesh in shock regions and coarser mesh in the smooth parts of the solution. The computational performance of these methods is demonstrated by considering scalar...
Adaptive Finite Element Relaxation Schemes for Hyperbolic Conservation Laws
We propose and study semidiscrete and fully discrete finite element schemes based on appropriate relaxation models for systems of Hyperbolic Conservation Laws. These schemes are using piecewise polynomials of arbitrary degree and their consistency error is of high order. The methods are combined with an adaptive strategy that yields fine mesh in shock regions and coarser mesh in the smooth parts of the solution. The computational performance of these methods is demonstrated by considering scalar...
Adaptive finite volume discretization of density driven flows in porous media.
Adaptive mesh refinement strategy for a non conservative transport problem
Long time simulations of transport equations raise computational challenges since they require both a large domain of calculation and sufficient accuracy. It is therefore advantageous, in terms of computational costs, to use a time varying adaptive mesh, with small cells in the region of interest and coarser cells where the solution is smooth. Biological models involving cell dynamics fall for instance within this framework and are often non conservative to account for cell division. In that case...