The inverse problem for mathematical models of heart excitation

The inverse problem for mathematical models of heart excitation is stated; this problem is to determine the initial condition in the initial-boundary value problem for an evolutionary system of partial differential equations given the volume potential whose density is determined by the solution to the evolutionary system. It is proved that the solution of the inverse problem in the generic statement is not unique.     

Differential quadrature for multi-dimensional problems

The technique of differential quadrature for the solution of partial differential equations, introduced by Bellman et al., is extended and generalized to encompass partial differential equations involving multiple space variables. Approximation formulae for a variety of first and second order partial derivatives and typical weighting coefficients are presented. Application of these formulae is demonstrated on the solution of the convection-diffusion equation for the two- and three-dimensional space dependent cases and for both the transient and steady-state dispersion of inert, neutrally buoyant pollutants from continuous sources into an unbounded atmosphere.     

A solvable hyperbolic free boundary problem modelling tumour regrowth

Recently, Tian and Friedman et al. developed a mathematical model on brain tumour recurrence after resection [J.P. Tian, A. Friedman, J. Wang and E.A. Chiocca, Modeling the effects of resection, radiation and chemotherapy in glioblastoma, J. Neuro-Oncol. 91(3) (2009), pp. 287–293]. The model is a free boundary problem with a hyperbolic system of nonlinear partial differential equations. In this article, we conduct a rigorous analysis on this hyperbolic system and prove the local and global existence and uniqueness of the solution. It is well known that most nonlinear free boundary problems are impossible to solve in terms of explicit analytical solutions. In contrast, the free boundary problem in this study is solvable, and the explicit solution is found using the backward characteristic curve method. This explicit solution is then validated by numerical simulation results. An interesting finding in this study is that the problem can be treated as a hyperbolic system defined on an infinite domain where the initial condition has a first-type discontinuity.     

一类渗透率反演问题解的存在性与唯一性

本根据Thierry Bourbie et al建立的测定致密岩心的渗透率的装置，交换相应的数学模型中的边界条件和附加条件位置，得到了相应正问题的解析解．尔后，运用偏微分方程反问题中的系数反演方法，构造出了反演渗透率的关系式，在此基础上，运用不动点定理讨论了解析反演解的存在性与唯一性．反演的结果表明：只要在L端持续测量t1时间间隔，则所给的附加条件可以唯一确定渗透率．     

Time decay of the solution to the Cauchy problem for a three-dimensional model of nonsimple thermoelasticity

This paper is devoted to the investigation of the solution to the Cauchy problem for a system of partial differential equations describing thermoelasticity of nonsimple materials in a three-dimensional space. The model of linear dynamical thermoelasticity of nonsimple materials is considered as the system of partial differential equations of fourth order. In this paper, we proposed a convenient evolutionary method of approach to the system of equations of nonsimple thermoelasticity. We proved the L^p−L^q time decay estimates for the solution to the Cauchy problem for linear thermoelasticity of nonsimple materials.     

Parabolic First and Second Order Differential Equations

We are concerned with an inverse problem for a first-order linear evolution equation. Moreover, a complete second-order evolution equation will be considered, too. We indicate sufficient conditions for existence and uniqueness of a solution. All the results apply well to inverse problems for equations from mathematical physics. As a possible application of the abstract theorems, some examples of partial differential equations are given.     

Degenerate first-order inverse problems in Banach spaces

We are concerned with an inverse problem for a degenerate linear evolution equation of first-order. Both hyperbolic and parabolic cases will be considered. We indicate sufficient conditions for the existence and uniqueness of a solution. All the results can be applied to inverse problems for equations from mathematical physics. As a possible application of the abstract theorems, some examples of partial differential equations are given.     

Integrability of the Equations for Nonsingular Pairs of Compatible Flat Metrics

We solve the problem of describing all nonsingular pairs of compatible flat metrics (or, in other words, nonsingular flat pencils of metrics) in the general N-component case. This problem is equivalent to the problem of describing all compatible Dubrovin–Novikov brackets (compatible nondegenerate local Poisson brackets of hydrodynamic type) playing an important role in the theory of integrable systems of hydrodynamic type and also in modern differential geometry and field theory. We prove that all nonsingular pairs of compatible flat metrics are described by a system of nonlinear differential equations that is a special nonlinear differential reduction of the classical Lamé equations, and we present a scheme for integrating this system by the method of the inverse scattering problem. The integration procedure is based on using the Zakharov method for integrating the Lamé equations (a version of the inverse scattering method).     

Preconditioners Based on Fundamental Solutions

We consider a new preconditioning technique for the iterative solution of linear systems of equations that arise when discretizing partial differential equations. The method is applied to finite difference discretizations, but the ideas apply to other discretizations too. If E is a fundamental solution of a differential operator P, we have E*(Pu) = u. Inspired by this, we choose the preconditioner to be a discretization of an approximate inverse K, given by a convolution-like operator with E as a kernel. We present analysis showing that if P is a first order differential operator, KP is bounded, and numerical results show grid independent convergence for first order partial differential equations, using fixed point iterations. For the second order convection-diffusion equation convergence is no longer grid independent when using fixed point iterations, a result that is consistent with our theory. However, if the grid is chosen to give a fixed number of grid points within boundary layers, the number of iterations is independent of the physical viscosity parameter. AMS subject classification (2000) 65F10, 65N22     

Finite-difference methods for solving the reaction-diffusion equations of a simple isothermal chemical system

Numerical methods are proposed for the numerical solution of a system of reaction-diffusion equations, which model chemical wave propagation. The reaction terms in this system of partial differential equations contain nonlinear expressions. Nevertheless, it is seen that the numerical solution is obtained by solving a linear algebraic system at each time step, as opposed to solving a nonlinear algebraic system, which is often required when integrating nonlinear partial differential equations. The development of each numerical method is made in the light of experience gained in solving the system of ordinary differential equations, which model the well-stirred analogue of the chemical system. The first-order numerical methods proposed for the solution of this initialvalue problem are characterized to be implicit. However, in each case it is seen that the numerical solution is obtained explicitly. In a series of numerical experiments, in which the ordinary differential equations are solved first of all, it is seen that the proposed methods have superior stability properties to those of the well-known, first-order, Euler method to which they are compared. Incorporating the proposed methods into the numerical solution of the partial differential equations is seen to lead to two economical and reliable methods, one sequential and one parallel, for solving the travelling-wave problem. © 1994 John Wiley & Sons, Inc.     

Integrability and beyond

An algebraic approach to solving nonlinear functional equations in the Riemann theta functions is stated. By the inverse scattering method and some general methods of the theory of partial differential equations, the solution of the initial boundary value problem for the nonlinear Schrödinger equation is presented. Bibliography:17 titles.     

Backward SDEs and Cauchy problem for semilinear equations in divergence form

 We extend the definition of solutions of backward stochastic differential equations to the case where the driving process is a diffusion corresponding to symmetric uniformly elliptic divergence form operator. We show existence and uniqueness of solutions of such equations under natural assumptions on the data and show its connections with solutions of semilinear parabolic partial differential equations in Sobolev spaces. Received: 22 January 2002 / Revised version: 10 September 2002 / Published online: 19 December 2002 Research supported by KBN Grant 0253 P03 2000 19. Mathematics Subject Classification (2002): Primary 60H30; Secondary 35K55 Key words or phrases: Backward stochastic differential equation – Semilinear partial differential equation – Divergence form operator – Weak solution     

Thermo-chemical modelling of the formation of a composite material

We consider a composite material composed of carbon or glass fibres included in a resin which becomes solid when it is heated up (the reaction of reticulation).

A mathematical model of the cure process is given by a kinetic equation describing the evolution of the reaction of reticulation coupled with the heat equation. The geometry of the composite material is periodic, with a small period ? >0, thus we get a coupled system of nonlinear partial differential equations.

First we prove the existence and uniqueness of a solution by using a fixed point theorem and we obtain a priori estimates. Then we derive the homogenized problem which describes the macroscopic behaviour of the material. We prove the convergence of the solution of the problem to the solution of the homogenized problem when ? tends to zero as well as the estimates for the difference of the exact and the approximate solutions.     

Deterministic algorithms for matrix completion

The goal of the matrix completion problem is to retrieve an unknown real matrix from a small subset of its entries. This problem comes up in many application areas, and has received a great deal of attention in the context of the Netflix challenge. This setup usually represents our partial knowledge of some information domain. Unknown entries may be due to the unavailability of some relevant experimental data. One approach to this problem starts by selecting a complexity measure of matrices, such as rank or trace norm. The corresponding algorithm outputs a matrix of lowest possible complexity that agrees with the partially specified matrix. The performance of the above algorithm under the assumption that the revealed entries are sampled randomly has received considerable attention (e.g., Refs. Srebro et al., 2005; COLT, 2005; Foygel and Srebro, 2011; Candes and Tao, 2010; Recht, 2009; Keshavan et al., 2010; Koltchinskii et al., 2010). Here we ask what can be said if the observed entries are chosen deterministically. We prove generalization error bounds for such deterministic algorithms, that resemble the results of Refs. Srebro et al. (2005); COLT (2005); Foygel and Srebro (2011) for the randomized algorithms. We still do not understand which sets of entries in a given matrix can be used to properly reconstruct it. Our hope is that the present work sheds some light on this problem. © 2013 Wiley Periodicals, Inc. Random Struct. Alg., 45, 306–317, 2014     

On the resolution and optimization of a system of fuzzy relational equations with sup-T composition

This paper provides a thorough investigation on the resolution of a finite system of fuzzy relational equations with sup-T composition, where T is a continuous triangular norm. When such a system is consistent, although we know that the solution set can be characterized by a maximum solution and finitely many minimal solutions, it is still a challenging task to find all minimal solutions in an efficient manner. Using the representation theorem of continuous triangular norms, we show that the systems of sup-T equations can be divided into two categories depending on the involved triangular norm. When the triangular norm is Archimedean, the minimal solutions correspond one-to-one to the irredundant coverings of a set covering problem. When it is non-Archimedean, they only correspond to a subset of constrained irredundant coverings of a set covering problem. We then show that the problem of minimizing a linear objective function subject to a system of sup-T equations can be reduced into a 0–1 integer programming problem in polynomial time. This work generalizes most, if not all, known results and provides a unified framework to deal with the problem of resolution and optimization of a system of sup-T equations. Further generalizations and related issues are also included for discussion.     

Recovery of differential equations with nonlinear dependence on the spectral parameter

The inverse spectral problem of recovering pencils of second-order differential operators on the half-line is studied. We give a formulation of the inverse problem, prove the uniqueness theorem and provided a procedure for constructing the solution of the inverse problem. We also establishe connections with inverse problems for partial differential equations.     

Solution of inverse diffusion problems by operator-splitting methods

In this paper the inverse solution of the general (non-linear) diffusion problem or backward heat conduction problem. It is assumed that the direct solution can be satisfactorily modelled, for example by the finite difference method. The nature of the problem and typical approaches to its solution are briefly reviewed.An operator-splitting method is introduced as a means of solving the inverse diffusion problem. An error analysis of the method is given, particularly for the application of the method to the simple diffusion equation. The method is applied to a range of test problems to illustrate the points of the analysis and to demonstrate the properties and performance of the method.     

Continuation of Galerkin approximations for reaction-diffusion problems

This paper deals with the use of the Galerkin approximation for calculating branches of steady-state solutions. It is motivated by the analysis of a reaction-diffusion system modeled by a pair of nonlinear partial differential equations on a two-dimensional domain. The goal is to check the possibility of closed loops emerging from a trivial branch. This issue is of importance in recent theories on morphogenesis in embryos (Kauffman et al. [3]). Numerical methods for continuing Galerkin approximations of the steady states give arcs of stable or unstable solutions. The numerical results are in agreement with the predictions of Brezzi et al. [6–8]. In particular, bifurcations from the trivial steady-state or symmetry-breaking bifurcations remain bifurcations for the approximate problem.The whole connected set of solutions thus obtained gives new insight into the behavior of solutions to reaction-diffusion equations and strongly advocates Kauffman's theory.     

Existence and uniqueness results for a nonlinear differential system

We study the existence, uniqueness and asymptotic behavior of the strong and weak solutions of a nonlinear differential system with 2N equations in a real Hilbert space H, subject to a boundary condition and initial data. This problem is a discrete version with respect to spatial variable x of some partial differential problems (with H = ℝⁿ ), which have applications in integrated circuits modelling     

A nonlinear diffusion model for image restoration

In this paper, we study the ill posed Perona-Malik equation of image processing^[14] and the regularized P-M model i.e. C-model proposed by Catte et al.^[4]. The authors present the convex compound of these two models in the form of the system of partial differential equations. The weak solution for the equations is proved in detail. The additive operator splitting (AOS) algorithm for the proposed model is also given. Finally, we show some numeric experimental results on images.