A streaming approach for sparse matrix products and its application in Galerkin multigrid methods.
A study of Galerkin method for the heat convection equations
The paper investigates the Galerkin method for an initial boundary value problem for heat convection equations. New error estimates for the approximate solutions and their derivatives in strong norm are obtained.
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 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...
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 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 modeling for free-surface flows
This work represents a first step towards the simulation of the motion of water in a complex hydrodynamic configuration, such as a channel network or a river delta, by means of a suitable “combination” of different mathematical models. In this framework a wide spectrum of space and time scales is involved due to the presence of physical phenomena of different nature. Ideally, moving from a hierarchy of hydrodynamic models, one should solve throughout the whole domain the most complex model (with...
Adattività nel metodo degli elementi finiti
Additive Schwarz algorithms for parabolic convection-diffusion equations.
An a posteriori error analysis for dynamic viscoelastic problems
In this paper, a dynamic viscoelastic problem is numerically studied. The variational problem is written in terms of the velocity field and it leads to a parabolic linear variational equation. A fully discrete scheme is introduced by using the finite element method to approximate the spatial variable and an Euler scheme to discretize time derivatives. An a priori error estimates result is recalled, from which the linear convergence is derived under suitable regularity conditions. Then, an a posteriori...
An a posteriori error analysis for dynamic viscoelastic problems
In this paper, a dynamic viscoelastic problem is numerically studied. The variational problem is written in terms of the velocity field and it leads to a parabolic linear variational equation. A fully discrete scheme is introduced by using the finite element method to approximate the spatial variable and an Euler scheme to discretize time derivatives. An a priori error estimates result is recalled, from which the linear convergence is derived under suitable regularity conditions. Then, an a posteriori error...
An adaptive finite element method in reconstruction of coefficients in Maxwell's equations from limited observations
We propose an adaptive finite element method for the solution of a coefficient inverse problem of simultaneous reconstruction of the dielectric permittivity and magnetic permeability functions in the Maxwell's system using limited boundary observations of the electric field in 3D. We derive a posteriori error estimates in the Tikhonov functional to be minimized and in the regularized solution of this functional, as well as formulate the corresponding adaptive algorithm. Our numerical experiments...