The scheme for the numerical solution of the incompressible Navier-Stokes equations coupled with the equation for temperature through the temperature dependent viscosity and thermal conductivity coefficients is presented. It is applied, together with the spectral element method, to the 2D calculations of flow around heated cylinder. High order polynomial approximation is combined with the decomposition of whole computational domain to only a few elements. Resulting data are compared with the experimental...

We consider in this paper a mathematical and numerical model to design an industrial software solution able to handle real complex furnaces configurations in terms of geometries, atmospheres, parts positioning, heat generators and physical thermal phenomena. A three dimensional algorithm based on stabilized finite element methods (SFEM) for solving the momentum, energy, turbulence and radiation equations is presented. An immersed volume method (IVM) for thermal coupling of fluids and solids is introduced...

In this paper, we introduce a set of methods for processing and analyzing long time series of 3D images representing embryo evolution. The images are obtained by in vivo scanning using a confocal microscope where one of the channels represents the cell nuclei and the other one the cell membranes. Our image processing chain consists of three steps: image filtering, object counting (center detection) and segmentation. The corresponding methods are based on numerical solution of nonlinear PDEs, namely...

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.