The extrapolation of Padé approximants in the closed-loop simulation of human thermoregulation

This paper presents a mathematical model of the human thermoregulatory system which has been developed to simulate the reaction of the human body to instantaneous changes in the temperature of the environment. The model combines two widely accepted approaches: the expression of heat flow within and from the human body in the form of partial differential equations, and the use of control theory. The closed-loop simulation is generated by introducing physical phenomena such as sweating, shivering, vasodilation and vasoconstriction, thus transforming the model equations into a system of nonlinear partial differential equations. The central difference operator is used to approximate spatial partial derivatives at a large number of mesh points, thus achieving high accuracy in the space dimension, and the problem is reduced to that of solving a nonlinear system of ordinary differential equations. A novel extrapolation algorithm, which copes well with the discontinuities between initial and boundary conditions caused by instantaneous environmental changes, is used to produce a high accuracy approximation in the time domain, and generate a solution to the problem. The model is used to simulate several important physical problems.     

Qualitative analysis of the invasion free boundary problem in biofilms

The work presents the qualitative analysis of the free boundary value problem related to the invasion model for multispecies biofilms. This model is based on the continuum approach for biofilm modeling and consists of a system of nonlinear hyperbolic partial differential equations for microbial species growth and spreading, a system of semilinear elliptic partial differential equations describing the substrate trends and a system of semilinear elliptic partial differential equations accounting for the diffusion and reaction of motile species within the biofilm. The free boundary evolution is regulated by a nonlinear ordinary differential equation. Overall, this leads to a free boundary value problem essentially hyperbolic. By using the method of characteristics, the partial differential equations constituting the invasion model are converted to Volterra integral equations. Then, the fixed point theorem is used for the uniqueness and existence result. The work is completed with numerical simulations describing the invasion of nitrite oxidizing bacteria in a biofilm initially constituted by ammonium oxidizing bacteria.     

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.     

一个药物作用肿瘤生长自由边界问题整体解的存在唯一性

该文研究一个描述药物作用下肿瘤生长的数学模型,这个肿瘤模型是对Jackson模型的一个改进,其数学形式是由一个二阶非线性抛物型方程与两个一阶非线性偏微分方程组耦合而成的自由边界问题.通过运用抛物型方程的L~P理论与一阶偏微分方程的特征方法,并利用Banach不动点定理,证明了该问题存在唯一的整体经典解.     

Hopf-Cole transformation to some systems of partial differential equations

In this paper we study a system of nonlinear partial differential equations which we write as a Burgers equation for matrix and use the Hopf-Cole transformation to linearize it. Using this method we solve initial value problem and initial boundary value problems for some systems of parabolic partial differential equations. Also we study an initial value problem for a system of nonlinear partial differential equations of first order which does not have solution in the standard distribution sense and construct an explicit solution in the algebra of generalized functions of Colombeau. Received November 1999     

Study of the solvability of a boundary value problem for the system of nonlinear differential equations of the theory of shallow shells of the Timoshenko type

We study the solvability of a boundary value problem for a system of nonlinear second-order partial differential equations under given boundary conditions, which describes the equilibrium of elastic shallow shells with hinged edges in the framework of the Timoshenko shear model. The study method implies the reduction of the original system of equations to a single nonlinear differential equation whose solvability is proved with the use of the contraction mapping principle.     

Free boundary problem for an initial cell layer in multispecies biofilm formation

The initial attached cell layer in multispecies biofilm growth is considered. The corresponding mathematical model leads to discuss a free boundary problem for a system of nonlinear hyperbolic partial differential equations, where the initial biofilm thickness is equal to zero. No assumptions on initial conditions for biomass concentrations and biofilm thickness are required. The data that the problem needs are the concentration of biomass in the bulk liquid and biomass flux from the bulk liquid. The method of characteristics is used to convert the differential system to Volterra integral equations for which an existence and uniqueness theorem is proved. Subsequently, we show that the free boundary is an increasing function of time and biomass concentrations are positive in agreement with the biological process.     

A Chebyshev spectral collocation method for solving Burgers’-type equations

In this paper, we elaborated a spectral collocation method based on differentiated Chebyshev polynomials to obtain numerical solutions for some different kinds of nonlinear partial differential equations. The problem is reduced to a system of ordinary differential equations that are solved by Runge–Kutta method of order four. Numerical results for the nonlinear evolution equations such as 1D Burgers’, KdV–Burgers’, coupled Burgers’, 2D Burgers’ and system of 2D Burgers’ equations are obtained. The numerical results are found to be in good agreement with the exact solutions. Numerical computations for a wide range of values of Reynolds’ number, show that the present method offers better accuracy in comparison with other previous methods. Moreover the method can be applied to a wide class of nonlinear partial differential equations.     

THE NONLOCAL INITIAL PROBLEMS OF A SEMILINEAR EVOLUTION EQUATION

The purpose of this paper is to investigate the existence of solutions to a nonlocal Cauchy problem for an evolution equation. The methods used here include the semigroup methods in proper spaces and Schauder's theorem. And the results are applied to a system of nonlinear partial differential equations with nonlinear boundary conditions.     

The method of lines for Poisson's equation with nonlinear or free boundary conditions

Summary The method of lines is used to solve Poisson's equation on an irregular domain with nonlinear or free boundary conditions. The partial differential equation is approximated by a system of second order ordinary differential equations subject to multi-point boundary conditions. The system is solved with an SOR iteration which employs invariant imbedding for each one dimensional problem. An application of the method to a boundary control problem and to a free surface problem arising in electrochemical machining is described. Finally, some theoretical convergence results are presented for a model problem with radiative boundary conditions on fixed boundaries.This work was supported by the U.S. Army Research Office under Grant DA-AG29-76-G-0261     

Similarity representation of wave propagation in a nonlinear viscoelastic rod on a group theoretic basis

Similarity representation of wave propagation in a nonlinear viscoelastic rod, subjected to a velocity impact is constructed. By the use of a multiparameter dimensional group of transformations, under which the system of basic equations and auxiliary conditions are invariant, similarity transformations are obtained which are used to construct the similarity representation in the form of an ordinary nonlinear boundary value problem of the original system, consisting of a set of simultaneous nonlinear partial differential equations and the association boundary and initial conditions.     

Schwarz problem for higher‐order complex partial differential equations in the upper half plane

Linear and nonlinear elliptic complex partial differential equations of higher‐order are considered under Schwarz conditions in the upper‐half plane. Firstly, using the integral representations for the solutions of the inhomogeneous polyanalytic equation with Schwarz conditions, a class of integral operators is introduced together with some of their properties. Then, these operators are used to transform the problem for linear equations into singular integral equations. In the case of nonlinear equations such a transformation yields a system of integro‐differential equations. Existence of the solutions of the relevant boundary value problems for linear and nonlinear equations are discussed via Fredholm theory and fixed point theorems, respectively.     

半导体问题的向后差分配置法及误差分析

Collocation method is put forward to solve the semiconductor problem with heat-conduction, whose mathematical model is described by an initial and boundary problem for a nonlinear partial differential equation system. One elliptic equation is for the electric potential, and three parabolic equations are for the electron concentration, hole concentration and heat-conduction. Using the prior estimate and technique of differential equations, we obtained almost optimal error estimates in L2.     

The problem of the nonlinear diffusive predator–prey model with the same biotic resource

In recent years, the research on the diffusive predator–prey model has attracted much attention. In these models, the carrying capacity is considered as a constant. In 2013, H. M. Safuan investigated the system of a predator and prey that shares the same biotic resource, where the carrying capacity is a function of the time. The spatial component of ecological interactions has been recognized as an important factor. So, we will discuss the problem of the nonlinear diffusive predator–prey model with the same biotic resource. This model is the system of the nonlinear partial differential equations with zero-flux boundary condition. The main objective of the present paper is to investigate the existence and uniqueness of the solution of this model. In this paper, we also obtain that there is a unique solution of the nonlinear partial differential equations with Dirichlet boundary condition.     

在纵向自由薄膜凝固期内温度分布的数值分析

考虑到薄膜表面的热通量主要是来自辐射，因而采用一个依赖时间的拟二维拟线性扩散方程的Ｓｔｅｆａｎ问题（混合初边值问题）作为该问题的数学模型。用一种具有Ｃｒａｎｋ－Ｎｉｃｈｏｌｓｏｎ格式的无条件稳定的有限差分析来求解抛物型偏微分方程的定解问题。用追赶法求解离散化的三对角方程组，然后用预估校正法求解拟线性扩散方程，从而避免了示解非线性差分方程组，琥珀亚硝酸酯在纵向自由薄膜凝固期内的温度分布数值计算结果和     

ON A SEMILINEAR MRABOLIC EQUATION SYSTEM WITH NONLOCAL AND COUPLED BOUNDARY CONDITIONS

1IntroductionInthestudyofquasi-statethermoelasticity,Deng[1-2]derivedamathemati-calmodelwhichinvolvesalinearparabolicequationwithanonlocalboundarycondition.Thismodelhasbeenextendedtomoregeneralsemilinearparabo-licequationsinhigh-dimensiondomainsbyFriedman[5]andKawohlI6],andmorerecentlybyDeng[3],Yin[13],Paol8-lo]andWang[11],andvariouscomparison,estimateandstabilityresultshavebeenobtained.InhtispaperweextendtheproblemofPao[9]tothefollowingproblemwithmoregeneralcoupledboundaryconditions(PE):w…     

The geometrically nonlinear dynamic responses of simply supported beams under moving loads

This paper presents a method for determining the nonlinear dynamic responses of structures under moving loads. The load is considered as a four degrees-of-freedom system with linear suspensions and tires flexibility, and the structure is modeled as an Euler–Bernoulli beam with simply supported at both ends. The nonlinear dynamic interaction of the load–structure system is discussed, and Kelvin−Voigt material model is employed for the beam. The nonlinear partial differential equations of the dynamic interaction are derived by using the von Kármán nonlinear theory and D'Alembert's principle. Based on the Galerkin method, the partial differential equations of the system are transformed into nonlinear ordinary equations, which can be solved by using the Newmark method and Newton−Raphson iteration method. To validate the approach proposed in this paper, the comparison are performed using a moving mass and a moving oscillator as the excitation sources, and the investigations demonstrate good reliability.     

An Age-structured Mathematical Model for the Within Host Dynamics of Malaria and the Immune System

In the paper, we use a mathematical model to study the population dyna mics of replicating malaria parasites and their interaction with the immune cells within a human host. The model is formulated as a system of age-structured partial differential equations that are then integrated over age to obtain a system of nonlinear delay differential equations. Our model incorporates an intracellular time delay between the infection of the red blood cells by the merozoites that grow and replicate within the infected cells to produce new merozoites. The infected red blood cells burst approximately every 48 h releasing daughter parasites to renew the cycle. The dynamical processes of the parasites within the human host are subjected to pressures exerted by the human immunological responses. The system is then solved using a first-order, finite difference method to give a discrete system. Numerical simulations carried out to illustrate stability of the system reveal that the populations undergo damped oscillations that stabilise to steady states.      

Integral equations related to the study of an inverse coefficient problem for a system of partial differential equations

We consider an inverse coefficient problem for a linear system of partial differential equations. The values of one solution component on a given curve are used as additional information for determining the unknown coefficient. The proof of the uniqueness of the solution of the inverse problem is based on the analysis of the unique solvability of a homogeneous integral equation of the first kind. The existence of a solution of the inverse problem is proved by reduction to a system of nonlinear integral equations.     

UPWIND FINITE VOLUME SCHEMES FOR SEMICONDUCTOR DEVICE

The mathematical model of semiconductor devices is described by the initial boundary value problem of a system of three nonlinear partial differential equations. One equation in elliptic form is for the electrostatic potential; two equations of convection-dominated diffusion type are for the electron and hole concentrations. Finite volume element procedure are put forward for the electrostatic potential, while upwind