We present a new semi-discrete central scheme for one-dimensional shallow water flows along channels with non-uniform rectangular cross sections and bottom topography. The scheme preserves the positivity of the water height, and it is preserves steady-states of rest (i.e., it is well-balanced). Along with a detailed description of the scheme, numerous numerical examples are presented for unsteady and steady flows. Comparison with exact solutions illustrate the accuracy and robustness of the numerical...

We present a new semi-discrete central scheme for one-dimensional shallow water flows along channels with non-uniform rectangular cross sections and bottom topography. The scheme preserves the positivity of the water height, and it is preserves steady-states of rest (i.e., it is well-balanced). Along with a detailed description of the scheme, numerous numerical examples are presented for unsteady and steady flows. Comparison with exact solutions illustrate the accuracy and robustness of the numerical...

In this paper, a phase field system of Penrose–Fife type with non–conserved order parameter is considered. A class of time–discrete schemes for an initial–boundary value problem for this phase–field system is presented. In three space dimensions, convergence is proved and an error estimate linear with respect to the time–step size is derived.

The aim of this paper is to demonstrate how the variational equations from can be formulated and solved in some abstract Banach spaces without any a priori construction of special linearization schemes. This should be useful e.g. in the analysis of heat conduction problems and modelling of flow in porous media.

We analyze a nonlinear discrete scheme depending on second-order finite differences. This is the two-dimensional analog of a scheme which in one dimension approximates a free-discontinuity energy proposed by Blake and Zisserman as a higher-order correction of the Mumford and Shah functional. In two dimension we give a compactness result showing that the continuous problem approximating this difference scheme is still defined on special functions with bounded hessian, and we give an upper and a lower...

We analyze a nonlinear discrete scheme depending on second-order finite differences. This is the two-dimensional analog of a scheme which in one dimension approximates a free-discontinuity energy proposed by Blake and Zisserman as a higher-order correction of the Mumford and Shah functional. In two dimension we give a compactness result showing that the continuous problem approximating this difference scheme is still defined on special functions...

In this paper, we study a model for the magnetization in thin ferromagnetic films. It comes as a variational problem for $S^1$-valued maps $m^{\text{'}}$ (the magnetization) of two variables $x^{\text{'}}$: $E_{\epsilon }\left(m^{\text{'}}\right)=\epsilon \int {|{\nabla }^{\text{'}}·m^{\text{'}}|}^{2}d{x}^{\text{'}}+\frac{1}{2}\int {\left||{\nabla }^{\text{'}}{|}^{-1/2}{\nabla }^{\text{'}}·m^{\text{'}}\right|}^{2}d{x}^{\text{'}}$. We are interested in the behavior of minimizers as $\epsilon \to 0$. They are expected to be $S^1$-valued maps $m^{\text{'}}$ of vanishing distributional divergence ${\nabla }^{\text{'}}·m^{\text{'}}=0$, so that appropriate boundary conditions enforce line discontinuities. For finite $\epsilon >0$, these line discontinuities are approximated by smooth transition layers, the so-called Néel walls. Néel...

The aim of this work is to compare a new uncoupled solver for the cardiac Bidomain model with a usual coupled solver. The Bidomain model describes the bioelectric activity of the cardiac tissue and consists of a system of a non-linear parabolic reaction-diffusion partial differential equation (PDE) and an elliptic linear PDE. This system models at macroscopic level the evolution of the transmembrane and extracellular electric potentials of the anisotropic cardiac tissue. The evolution equation is...

We compare numerical experiments from the String Gradient Weighted Moving Finite Element method and a Parabolic Moving Mesh Partial Differential Equation method, applied to three benchmark problems based on two different partial differential equations. Both methods are described in detail and we highlight some strengths and weaknesses of each method via the numerical comparisons. The two equations used in the benchmark problems are the viscous Burgers’ equation and the porous medium equation, both...

The paper has been presented at the 12th International Conference on Applications of Computer Algebra, Varna, Bulgaria, June, 2006A MATHEMATICA package for finding Lie symmetries of partial differential equations is presented. The package is designed to create and solve the associated determining system of equations, the full set of solutions of which generates the widest permissible local Lie group of point symmetry transformations. Examples illustrating the functionality of the package's tools...

A method to approximate the Euler equations is presented. The method is a multi-domain approximation, and a variational form of the Euler equations is found by making use of the divergence theorem. The method is similar to that of the Discontinuous-Galerkin method of Cockburn and Shu, but the implementation is constructed through a spectral, multi-domain approach. The method is introduced and is shown to be a conservative scheme. A numerical example is given for the expanding flow around a point...

In this paper we study the finite element approximations to the Sobolev and viscoelasticity type equations and present a direct analysis for global superconvergence for these problems, without using Ritz projection or its modified forms.

The aim of this paper is to develop a crowd motion model designed to handle highly packed situations. The model we propose rests on two principles: we first define a spontaneous velocity which corresponds to the velocity each individual would like to have in the absence of other people. The actual velocity is then computed as the projection of the spontaneous velocity onto the set of admissible velocities (i.e. velocities which do not violate the non-overlapping constraint). We describe here the...

The aim of this paper is to develop a crowd motion model designed to handle highly packed situations. The model we propose rests on two principles: we first define a spontaneous velocity which corresponds to the velocity each individual would like to have in the absence of other people. The actual velocity is then computed as the projection of the spontaneous velocity onto the set of admissible velocities (i.e. velocities which do not violate the non-overlapping constraint). We describe here...