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.
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...
The Fourier problem on planar domains with time variable boundary is considered using integral equations. A simple numerical method for the integral equation is described and the convergence of the method is proved. It is shown how to approximate the solution of the Fourier problem and how to estimate the error. A numerical example is given.
We examine a heterogeneous alternating-direction method for the approximate solution of the FENE Fokker–Planck equation from polymer fluid dynamics and we use this method to solve a coupled (macro-micro) Navier–Stokes–Fokker–Planck system for dilute polymeric fluids. In this context the Fokker–Planck equation is posed on a high-dimensional domain and is therefore challenging from a computational point of view. The heterogeneous alternating-direction scheme combines a spectral Galerkin method for...
The newly developed unifying discontinuous formulation named the correction procedure via reconstruction (CPR) for conservation laws is extended to solve the Navier-Stokes equations for 3D mixed grids. In the current development, tetrahedrons and triangular prisms are considered. The CPR method can unify several popular high order methods including the discontinuous Galerkin and the spectral volume methods into a more efficient differential form....
In this paper we construct analytic-numerical solutions for initial-boundary value systems related to the equation , where is an arbitrary square complex matrix and ia s matrix such that the real part of the eigenvalues of the matrix is positive. Given an admissible error and a finite domain , and analytic-numerical solution whose error is uniformly upper bounded by in , is constructed.
In this paper we consider the problem of detecting pollution in some non linear parabolic systems using the sentinel method. For this purpose we develop and analyze a new approach to the discretization which pays careful attention to the stability of the solution. To illustrate convergence properties we give some numerical results that present good properties and show new ways for building discrete sentinels.
In this article we introduce an adaptive multi-level method in space and time for convection diffusion problems. The scheme is based on a multi-level spatial splitting and the use of different time-steps. The temporal discretization relies on the characteristics method. We derive an a posteriori error estimate and design a corresponding adaptive algorithm. The efficiency of the multi-level method is illustrated by numerical experiments, in particular for a convection-dominated problem.
We study error estimates and their convergence rates for approximate solutions of spectral Galerkin type for the equations for the motion of a viscous chemical active fluid in a bounded domain. We find error estimates that are uniform in time and also optimal in the L2-norm and H1-norm. New estimates in the H(-1)-norm are given.
This work is concerned with the asymptotic analysis of a time-splitting scheme for the Schrödinger equation with a random potential having weak amplitude, fast oscillations in time and space, and long-range correlations. Such a problem arises for instance in the simulation of waves propagating in random media in the paraxial approximation. The high-frequency limit of the Schrödinger equation leads to different regimes depending on the distance of propagation, the oscillation pattern of the initial...