Correct initial boundary value problems for dispersive equations

This paper deals with correctness of initial boundary value problems for general dispersive equations of finite odd orders. For the Kawahara and KdV equations we prove existence, uniqueness and stability of strong global solutions in a bounded domain for different signs of a coefficient of the highest derivative as well as their asymptotics when the coefficient of the higher-order derivative in the Kawahara equation approaches zero.     

A problem with normal derivatives in boundary conditions for a system of differential equations

We seek for a solution to a system of differential equations, using linear relations connecting normal derivatives of the desired functions at the domain boundary.     

A modified sensitivity equation method for the Euler equations in presence of shocks

The continuous sensitivity equation method allows to quantify how changes in the input of a partial differential equation (PDE) model affect the outputs, by solving additional PDEs obtained by differentiating the model. However, this method cannot be used directly in the framework of hyperbolic PDE systems with discontinuous solution, because it yields Dirac delta functions in the sensitivity solution at the location of state discontinuities. This difficulty is well known from theoretical viewpoint, but only a few works can be found in the literature regarding the possible numerical treatment. Therefore, we investigate in this study how classical numerical schemes for compressible Euler equations can be modified to account for shocks when computing the sensitivity solution. In particular, we propose the introduction of a source term, that allows to remove the spikes associated to the Dirac delta functions in the numerical solution. Numerical studies exhibit a strong impact of the numerical diffusion on the accuracy of this strategy. Therefore, we propose an anti-diffusive numerical scheme coupled with the approximate Riemann solver of Roe for the state problem. For the sensitivity problem, two different numerical schemes are implemented and compared: one which takes into account the contact wave and another that neglects it. The effects of the numerical diffusion on the convergence of the schemes with respect to the grid are discussed. Finally, an application to uncertainty propagation is investigated and the different numerical schemes are compared.     

On the two-dimensional initial boundary value problem for the Navier-Stokes equations with discontinuous boundary data

We consider the initial boundary value problem for the Navier-Stokes equations with boundary conditions . We assume that may have jump discontinuities at finitely many points ξ1;. . .,ξm of the boundary ϖΩ of a bounded domain Ω ⊂ ℝ². We prove that this problem has a unique generalized solution in a finite time interval or for small initial and boundary data. The solution is found in a class of vector fields with infinite energy integral. The case of a moving boundary is also considered. Bibliography: 11 titles. Dedicated to O. A. Ladyzhenskaya on the occasion of her 70^th birthday. Translated fromZapiski Nauchnykh Seminarov POMI, Vol. 197, pp. 159–178, 1992. Translated by E. V. Frolova.     

A method for calculating invariant subspaces of symmetric hyperbolic equations

An algorithm is constructed for calculating invariant subspaces of symmetric hyperbolic systems arising in electromagnetic, acoustic, and elasticity problems. Discrete approximations are calculated for subspaces that correspond to minimal eigenvalues and smooth eigenfunctions. Difficulties related to the presence of an infinite-dimensional kernel in the differential operator are successfully handled. The efficiency of the algorithm is demonstrated using acoustics equations.     

Superconvergence analysis of the finite element method for nonlinear hyperbolic equations with nonlinear boundary condition

This paper discusses the semidiscrete finite element method for nonlinear hyperbolic equations with nonlinear boundary condition. The superclose property is derived through interpolation instead of the nonlinear H^1 projection and integral identity technique. Meanwhile, the global superconvergence is obtained based on the interpolated postprocessing techniques.     

Initial boundary value problems for quasilinear symmetric hyperbolic systems with characteristic boundary

This paper is devoted to initial boundary value problems for quasi-linear symmetric hyperbolic systems in a domain with characteristic boundary. It extends the theory on linear symmetric hyperbolic systems established by Friedrichs to the nonlinear case. The concept on regular characteristics and dissipative boundary conditions are given for quasilinear hyperbolic systems. Under some assumptions, an existence theorem for such initial boundary value problems is obtained. The theorem can also be applied to the Euler system of compressible flow. __________ Translated from Chinese Annals of Mathematics, Ser. A, 1982, 3(2): 223–232     

A wavelet operational method for solving fractional partial differential equations numerically

Fractional calculus is an extension of derivatives and integrals to non-integer orders, and a partial differential equation involving the fractional calculus operators is called the fractional PDE. They have many applications in science and engineering. However not only the analytical solution existed for a limited number of cases, but also the numerical methods are very complicated and difficult. In this paper, we newly establish the simulation method based on the operational matrices of the orthogonal functions. We formulate the operational matrix of integration in a unified framework. By using the operational matrix of integration, we propose a new numerical method for linear fractional partial differential equation solving. In the method, we (1) use the Haar wavelet; (2) establish a Lyapunov-type matrix equation; and (3) obtain the algebraic equations suitable for computer programming. Two examples are given to demonstrate the simplicity, clarity and powerfulness of the new method.     

On the numerical integration of multi-dimensional,initial boundary value problems for the Euler equations in quasi-linear form

A matricial formalism to solve multi-dimensional initial boundary values problems for hyperbolic equations written in quasi-linear based on the λ scheme approach is presented. The derivation is carried out for nonorthogonal, moving systems of curvilinear coordinates. A uniform treatment of the integration at the boundaries, when the boundary conditions can be expressed in terms of combinations of time or space derivatives of the primitive variables, is also presented. The methodology is validated against a two-dimensional test case, the supercritical flow through the Hobson cascade n.2, and in three-dimensional test cases such as the supersonic flow about a sphere and the flow through a plug nozzle. © 1998 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 14: 781–814, 1998     

Non-standard finite difference method applied to an initial boundary value problem describing hepatitis B virus infection

In this paper, two non-standard finite difference (NSFD) schemes are proposed for a mathematical model of hepatitis B virus (HBV) infection with spatial dependence. The dynamic properties of the obtained discretized systems are completely analyzed. Relying on the theory of M-matrix, we prove that the proposed NSFD schemes is unconditionally positive. Furthermore, we establish that the NSFD method used preserves all constant steady states of the corresponding continuous initial boundary value problem (IBVP) model. We prove that the conditions for those equilibria to be asymptotically stable are consistent with the continuous IBVP model independently of the numerical grid size. The global asymptotical properties of the HBV-free equilibrium of the proposed NSFD schemes are derived via the construction of a suitable discrete Lyapunov function, and coincides with the continuous system. This confirms that the discretized models are dynamically consistent since they maintain essential properties of the corresponding continuous IBVP model. Finally, numerical simulations are performed from which it is demonstrated that the proposed NSFD method is advantageous over the standard finite difference (SFD) method.     

The initial boundary value problems for the semilinear diffusion equations with data in L p spaces

In this paper, I have established conditions under which the existence and uniqueness of weak solutions of some Semilinear diffusion equations with initial and boundary data in fractional LL ^p spaces can be established in a bounded domain with smooth boundary The interest of this method relies on the fact that it is by successive approximations and hence amenable to numerical treatment. The paper also considers the semigroups theory on the existence of weak classical solutions.     

Elliptic equations with low regularity boundary data via the unified method

We use novel integral representations developed by the second author to prove certain rigorous results concerning elliptic boundary value problems in convex polygons. Central to this approach is the so-called global relation, which is a non-local equation in the Fourier space that relates the known boundary data to the unknown boundary values. Assuming that the global relation is satisfied in the weakest possible sense, i.e. in a distributional sense, we prove there exist solutions to Dirichlet, Neumann and Robin boundary value problems with distributional boundary data. We show that the analysis of the global relation characterises in a straightforward manner the possible existence of both integrable and non-integrable corner singularities.     

Global singularity structure of weak solutions to 3-D semilinear dispersive wave equations with discontinuous initial data

We study the global singularity structure of solutions to 3-D semilinear wave equations with discontinuous initial data. More precisely, using Strichartz’ inequality we show that the solutions stay conormal after nonlinear interaction if the Cauchy data are conormal along a circle.     

A class of nonlinear singularly perturbed initial boundary value problems for reaction diffusion equations with boundary perturbation

A class of nonlinear initial boundary value problems for reaction diffusion equations with boundary perturbation is considered. Under suitable conditions and using the theory of differential inequalities the asymptotic solution of the initial boundary value problems is studied.     

A general method for writing the dynamical equations of nonholonomic systems with ideal constraints

The basic notions of the dynamics of nonholonomic systems are revisited in order to give a general and simple method for writing the dynamical equations for linear as well as non-linear kinematical constraints. The method is based on the representation of the constraints by parametric equations, which are interpreted as dynamical equations, and leads to first-order differential equations in normal form, involving the Lagrangian coordinates and auxiliary variables (the use of Lagrangian multipliers is avoided). Various examples are illustrated.      

Method for solving hyperbolic systems with initial data on non-characteristic manifolds with applications to the shallow water wave equations

We are concerned with hyperbolic systems of order-one linear PDEs originated on non-characteristic manifolds. We put forward a simple but effective method of transforming such initial conditions to standard initial conditions (i.e. when the solution is specified at an initial moment of time). We then show how our method applies in fluid mechanics. More specifically, we present a complete solution to the problem of long waves run-up in inclined bays of arbitrary shape with nonzero initial velocity.     

A fast Fourier-collocation method for second boundary integral equations

In this paper we develop a fast collocation method for second boundary integral equations by the trigonometric polynomials. We propose a convenient way to compress the dense matrix representation of a compact integral operator with a smooth kernel under the Fourier basis and the corresponding collocation functionals. The compression leads to a sparse matrix with only O(nlog²n) number of nonzero entries, where 2n+1 denotes the order of the matrix. Thus we develop a fast Fourier-collocation method. We prove that the fast Fourier-collocation method gives the optimal convergence order up to a logarithmic factor. Moreover, we design a fast scheme for solving the corresponding truncated linear system. We establish that this algorithm preserves the quasi-optimal convergence of the approximate solution with requiring a number of O(nlog³n) multiplications.     

A Tau method with perturbed boundary conditions for certain ordinary differential equations

Ortiz' recursive formulation of the Lanczos Tau method (TM) is a powerful and efficient technique for producing polynomial approximations for initial or boundary value problems. The method consists in obtaining a polynomial which satisfies (i) aperturbed version of the given differential equation, and (ii) the imposed supplementary conditionsexactly. This paper introduces a new form of the TM, (denoted by PTM), for a restricted class of differential equations, in which the differential equations as well as the supplementary conditions areperturbed simultaneously. PTM is compared to the classical TM from the point of view of their errors: it is found that the PTM error is smaller and more oscillatory than that of the TM; we further find that approximations nearly as accurate as minimax polynomial approximations can be constructed by means of the PTM. Detailed formulae are derived for the polynomial approximations in TM and PTM, based on Canonical Polynomials. Moreover, various limiting properties of Tau coefficients are established and it is shown that the perturbation in PTM behaves asymptotically proprtional to a Chebyshev polynomial. Dedicated to Eduardo L. Ortiz on the occasion of his 70th birthday     

A free boundary problem for the Navier–Stokes equations with measure data

The existence of a weak solution of an initial boundary-value problem for the plane nonstationary Navier–Stokes equations with Radon measure data on the free boundary, is established. The problem may be considered as a model of the blood flow around the heart valves. An inverse problem is studied, it allows us to find the boundary forces acting on the valve from the observed values of the velocity of the fluid in a fixed subregion.