A spring-dashpot system for modelling lung tumour motion in radiotherapy.
We present in this paper a stability study concerning finite volume schemes applied to the two-dimensional Maxwell system, using rectangular or triangular meshes. A stability condition is proved for the first-order upwind scheme on a rectangular mesh. Stability comparisons between the Yee scheme and the finite volume formulation are proposed. We also compare the stability domains obtained when considering the Maxwell system and the convection equation.
It is well known that the classical local projection method as well as residual-based stabilization techniques, as for instance streamline upwind Petrov-Galerkin (SUPG), are optimal on isotropic meshes. Here we extend the local projection stabilization for the Navier-Stokes system to anisotropic quadrilateral meshes in two spatial dimensions. We describe the new method and prove an a priori error estimate. This method leads on anisotropic meshes to qualitatively better convergence behavior...
In this work we consider a stabilized Lagrange multiplier method in order to approximate the Coulomb frictional contact model in linear elastostatics. The particularity of the method is that no discrete inf-sup condition is needed. We study the existence and the uniqueness of solution of the discrete problem.
We outline a solution method for mixed finite element discretizations based on dissecting the problem into three separate steps. The first handles the inhomogeneous constraint, the second solves the flux variable from the homogeneous problem, whereas the third step, adjoint to the first, finally gives the Lagrangian multiplier. We concentrate on aspects involved in the first and third step mainly, and advertise a multi-level method that allows for a stable computation of the intermediate and final...
Method for numerical solution of Volterra integral equations, based on the O.I.M. methods, is suggested. It is known that the class of O.I.M. methods includes -stable methods of arbitrary high order of asymptotic accuracy. In part 5, it is proved that these methods generate methods for numerical solution of Volterra equations which are also -stable and of an arbitrarily high order. There is one advantage of the methods. Namely, they need no matrix inversion in the course of their numerical realization....
We propose a new reduced basis element-cum-component mode synthesis approach for parametrized elliptic coercive partial differential equations. In the Offline stage we construct a Library of interoperable parametrized reference components relevant to some family of problems; in the Online stage we instantiate and connect reference components (at ports) to rapidly form and query parametric systems. The method is based on static condensation at the interdomain level, a conforming eigenfunction “port”...
We propose a simple numerical method for capturing the steady state solution of hyperbolic systems with geometrical source terms. We use the interface value, rather than the cell-averages, for the source terms that balance the nonlinear convection at the cell interface, allowing the numerical capturing of the steady state with a formal high order accuracy. This method applies to Godunov or Roe type upwind methods but requires no modification of the Riemann solver. Numerical experiments on scalar...