A single cell-based model of the ductal tumour microarchitecture.
A sixth-order finite volume method for the 1D biharmonic operator: Application to intramedullary nail simulation
A new very high-order finite volume method to solve problems with harmonic and biharmonic operators for onedimensional geometries is proposed. The main ingredient is polynomial reconstruction based on local interpolations of mean values providing accurate approximations of the solution up to the sixth-order accuracy. First developed with the harmonic operator, an extension for the biharmonic operator is obtained, which allows designing a very high-order finite volume scheme where the solution is...
A Slideing Mesh-Mortar Method for a two Dimensional Currents Model of Electric Engines
The paper deals with the application of a non-conforming domain decomposition method to the problem of the computation of induced currents in electric engines with moving conductors. The eddy currents model is considered as a quasi-static approximation of Maxwell equations and we study its two-dimensional formulation with either the modified magnetic vector potential or the magnetic field as primary variable. Two discretizations are proposed, the first one based on curved finite elements and the...
A sliding Mesh-Mortar method for a two dimensional Eddy currents model of electric engines
The paper deals with the application of a non-conforming domain decomposition method to the problem of the computation of induced currents in electric engines with moving conductors. The eddy currents model is considered as a quasi-static approximation of Maxwell equations and we study its two-dimensional formulation with either the modified magnetic vector potential or the magnetic field as primary variable. Two discretizations are proposed, the first one based on curved finite elements and the...
A smoothing Levenberg-Marquardt method for the complementarity problem over symmetric cone
In this paper, we propose a smoothing Levenberg-Marquardt method for the symmetric cone complementarity problem. Based on a smoothing function, we turn this problem into a system of nonlinear equations and then solve the equations by the method proposed. Under the condition of Lipschitz continuity of the Jacobian matrix and local error bound, the new method is proved to be globally convergent and locally superlinearly/quadratically convergent. Numerical experiments are also employed to show that...
A Sobolev gradient method for treating the steady-state incompressible Navier-Stokes equations
The velocity-vorticity-pressure formulation of the steady-state incompressible Navier-Stokes equations in two dimensions is cast as a nonlinear least squares problem in which the functional is a weighted sum of squared residuals. A finite element discretization of the functional is minimized by a trust-region method in which the trustregion radius is defined by a Sobolev norm and the trust-region subproblems are solved by a dogleg method. Numerical test results show the method to be effective.
A software package for the numerical integration of ODEs by means of high-order Taylor methods.
A solution of the continuous Lyapunov equation by means of power series
A Solver for Complex-Valued Parametric Linear Systems
This work reports on a new software for solving linear systems involving affine-linear dependencies between complex-valued interval parameters. We discuss the implementation of a parametric residual iteration for linear interval systems by advanced communication between the system Mathematica and the library C-XSC supporting rigorous complex interval arithmetic. An example of AC electrical circuit illustrates the use of the presented software.* This work was partly supported by the DFG grant GZ:...
A SOR Acceleration of Self-Adjoint and m-Accretive Splitting Iterative Solver for 2-D Neutron Transport Equation
We present an iterative method based on an infinite dimensional adaptation of the successive overrelaxation (SOR) algorithm for solving the 2-D neutron transport equation. In a wide range of application, the neutron transport operator admits a Self-Adjoint and m-Accretive Splitting (SAS). This splitting leads to an ADI-like iterative method which converges unconditionally and is equivalent to a fixed point problem where the operator is a 2 by 2 matrix...
A spatially sixth-order hybrid -CCD method for solving time fractional Schrödinger equations
We consider highly accurate schemes for nonlinear time fractional Schrödinger equations (NTFSEs). While an strategy is employed for approximating the Caputo fractional derivative in the temporal direction, compact CCD finite difference approaches are incorporated in the space. A highly effective hybrid -CCD method is implemented successfully. The accuracy of this linearized scheme is order six in space, and order in time, where is the order of the Caputo fractional derivative involved. It...
A special finite element method based on component mode synthesis
The goal of our paper is to introduce basis functions for the finite element discretization of a second order linear elliptic operator with rough or highly oscillating coefficients. The proposed basis functions are inspired by the classic idea of component mode synthesis and exploit an orthogonal decomposition of the trial subspace to minimize the energy. Numerical experiments illustrate the effectiveness of the proposed basis functions.
A special Gaussian rule for trigonometric polynomials.
A spectral characterization of the behavior of discrete time AR–representations over a finite time interval
In this paper we investigate the behavior of the discrete time AR (Auto Regressive) representations over a finite time interval, in terms of the finite and infinite spectral structure of the polynomial matrix involved in the AR-equation. A boundary mapping equation and a closed formula for the determination of the solution, in terms of the boundary conditions, are also gived.
A spectral Galerkin approximation of the Orr-Sommerfeld eigenvalue problem in a semi-finite domain.
A spectral method for the eigenvalue problem for elliptic equations.
A spectral regularization method for a Cauchy problem of the modified Helmholtz equation.
A spectral-Tau approximation for the Stokes and Navier-Stokes equations
A splitting method using discontinuous Galerkin for the transient incompressible Navier-Stokes equations
In this paper we solve the time-dependent incompressible Navier-Stokes equations by splitting the non-linearity and incompressibility, and using discontinuous or continuous finite element methods in space. We prove optimal error estimates for the velocity and suboptimal estimates for the pressure. We present some numerical experiments.
A splitting method using discontinuous Galerkin for the transient incompressible Navier-Stokes equations
In this paper we solve the time-dependent incompressible Navier-Stokes equations by splitting the non-linearity and incompressibility, and using discontinuous or continuous finite element methods in space. We prove optimal error estimates for the velocity and suboptimal estimates for the pressure. We present some numerical experiments.