非线性耦合Klein-Gordon方程组和耦合Schr\"{o}dinger-Boussinesq方程组的双行波解

本文提出了一种全新复合$(\frac{G''}{G})$展开方法,运用这种新方法并借助符号计算软件构造了非线性耦合Klein-Gordon方程组和耦合Schr\"{o}dinger-Boussinesq方程组的多种双行波解,包括双双曲正切函数解,双正切函数解,双有理函数解以及它们的混合解. 复合$(\frac{G''}{G})$展开方法不但直接有效地求出了两类非线性偏微分方程的双行波解,而且扩大了解的范围.这种新方法对于研究非线性偏微分方程具有广泛的应用意义.     

Accurate spectral solutions for the parabolic and elliptic partial differential equations by the ultraspherical tau method

We present a double ultraspherical spectral methods that allow the efficient approximate solution for the parabolic partial differential equations in a square subject to the most general inhomogeneous mixed boundary conditions. The differential equations with their boundary and initial conditions are reduced to systems of ordinary differential equations for the time-dependent expansion coefficients. These systems are greatly simplified by using tensor matrix algebra, and are solved by using the step-by-step method. Numerical applications of how to use these methods are described. Numerical results obtained compare favorably with those of the analytical solutions. Accurate double ultraspherical spectral approximations for Poisson's and Helmholtz's equations are also noted. Numerical experiments show that spectral approximation based on Chebyshev polynomials of the first kind is not always better than others based on ultraspherical polynomials.     

The Moore-Penrose and group inverse of square matrices and the Jordan canonical form

This paper is a continuation of the paper [3], where the first two Penrose's equations are solved for a square matrix, which is transformed into the Jordan canonical form. Our main aim is to solve the other two Penrose's equations for the same class of matrices. We also find a block representation of the group inverse and the form of the corresponding Jordan matrix, by solving the corresponding set of equations.     

A transformation of the equations of mechanics

The structure of the equations obtained from Hamilton's equations by a coordinate Legendre transformation of the Hamilton function is discussed. Examples are considered.     

Newton's second law in field theory

In this article we present a natural generalization of Newton's Second Law valid in field theory, i.e., when the parameterized curves are replaced by parameterized submanifolds of higher dimension. For it we introduce what we have called the geodesic k-vector field, analogous to the ordinary geodesic field and which describes the inertial motions (i.e., evolution in the absence of forces). From this generalized Newton's law, the corresponding Hamilton's canonical equations of field theory (Hamilton-De Donder-Weyl equations) are obtained by a simple procedure. It is shown that solutions of generalized Newton's equation also hold the canonical equations. However, unlike the ordinary case, Newton equations determined by different forces can define equal Hamilton's equations.     

Local Stability of Euler's and Kahan's Methods

Two discretization methods, the forward Euler's method and the Kahan's reflexive method, are compared by looking at the local stabilities of fixed points of a system of differential equations. We explain why forward Euler's method is not as good from the viewpoint of complex analysis. Conformal mappings are used to relate the eigenvalues of the Jacobian matrices of the differential equations system and the resulting difference equations system. The Euler's method will not preserve Hopf bifurcation. The Kahan's method preserves the local stability of the fixed points of the differential equations.     

Use of Halley's Method in the Nonlinear Finite Element Analysis

Halley's method is a higher order iteration method for the solution of nonlinear systems of equations. Unlike Newton's method, which converges quadratically in the vicinity of the solution, Halley's method can exhibit a cubic order of convergence. The equations of Halley's method for multiple dimensions are derived using Padé approximants and inverse one-point interpolation, as proposed by Cuyt. The investigation of the performance of Halley's method concentrates on eight-node volume elements for nonlinear deformations using Staint Venant-Kirchhoff's constitutive law, as well as a geometric linear theory of von Mises plasticity. The comparison with Newton's method reveals the sensibility of Halley's method, in view of the radius of attraction but also demonstrates the advantages of Halley's method considering simulation costs and the order of convergence. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Analysis of Operator Models Arising in Problems of Hereditary Mechanics

The correct solvability of initial-value problems for integrodifferential equations with unbounded operator coefficients in Hilbert spaces is determined. Numerous problems of hereditary mechanics and thermal physics have motivated the study of such equations. Models taking into account a Kelvin–Voight friction are considered. A spectral analysis of the operator functions, which are the symbols of integro-differential equations considered in this paper, is made.     

The Euler Method in the Field of Operators

Approximate solutions of a class of nonlinear differential equations are constructed in the field of Mikusiński operators by using a method similar to Euler's. Their character is analyzed and the error of approximation is estimated. The results are applied to a class of partial integro–differential equations.     

Splitting scheme for poroelasticity and thermoelasticity problems

Boundary value problems in thermoelasticity and poroelasticity (filtration consolidation) are solved numerically. The underlying system of equations consists of the Lamé stationary equations for displacements and nonstationary equations for temperature or pressure in the porous medium. The numerical algorithm is based on a finite-element approximation in space. Standard stability conditions are formulated for two-level schemes with weights. Such schemes are numerically implemented by solving a system of coupled equations for displacements and temperature (pressure). Splitting schemes with respect to physical processes are constructed, in which the transition to a new time level is associated with solving separate elliptic problems for the desired displacements and temperature (pressure). Unconditionally stable additive schemes are constructed by choosing a weight of a three-level scheme.     

On approximations of the Euler-Peano scheme for multivalued stochastic differential equations

It is known that a unique strong solution exists for multivalued stochastic differential equations under the Lipschitz continuity and linear growth conditions. In this paper we apply the Euler-Peano scheme to show that existence of weak solution and pathwise uniqueness still hold when the coefficients are random and satisfy one-sided locally Lipschitz continuous and an integral condition (i.e. Krylov's conditions put forward in On Kolmogorov's equations for finite-dimensional diffusions, Stochastic PDE's and Kolmogorov Equations in Infinite Dimensions (Cetraro, 1998), Lecture Notes in Math., 1715, Springer, Berlin, 1999, pp. 1–63). When the coefficients are nonrandom and possibly discontinuous but only satisfy some integral conditions, the sequence of solutions of the Euler-Peano scheme converges weakly, and the limit is a weak solution of the corresponding MSDE. As a particular case, we obtain a global semi-flow for stochastic differential equations reflected in closed, convex domains.     

A complete steady state model of solute and water transport in the kidney

The purpose of this paper is to incorporate a detailed model, along with an optimized set of parameters for the proximal tubule, into J. L. Stephenson's current central core model of the nephron. In this model a set of equations for the proximal tubule are combined with Stephenson's equations for the remaining four tubules and interstitium, to form a complete nonlinear system of 34 ordinary differential and algebraic equations governing fluid and solute flow in the kidney. These equations are then discretized by the Crank-Nicholson scheme to form an algebraic system of nonlinear equations for the unknown concentrations, flows, hydrostatic pressure, and potentials. The resulting system is solved via factored secant update with a finite-difference approximation to the Jacobian. Finally, numerical simulations performed on the model showed that the modeled behavior approximates, in a general way, the physiological mechanisms of solvent and solute flow in the kidney.     

Computer-based proof of the existence of superconvergence points in the finite element method; superconvergence of the derivatives in finite element solutions of Laplace's,Poisson's,and the elasticity equations

In this article we address the problem of the existence of superconvergence points for finite element solutions of systems of linear elliptic equations. Our approach is quite different from all other studies of superconvergence. We prove that the existence of superconvergence points can be guaranteed by a numerical algorithm, which employs a finite number of operations (provided that there is no roundoff-error). By employing this approach, we can reproduce all known results on superconvergence of finite element solutions for linear elliptic problems and we can obtain many new results. Here, in particular, we address the problem of the superconvergence points for the gradient of finite element solutions of Laplace's and Poisson's equations and we show that the sets of superconvergence points are very different for these two cases. We also study the superconvergence of the components of the gradient of the displacement, the strain and stress for finite element solutions of the equations of elasticity. For Laplace's and Poisson's equations (resp. the equations of elasticity), we consider meshes of triangular as well as square elements of degree p, 1 ? p ? 7 (resp. 1 ? p ? 4). For the meshes of triangular elements we investigate the effect of the geometry of the mesh by considering four mesh patterns that typically occur in practical meshes, while in the case of square elements, we study the effect of the element type (tensor-product, serendipity, or other). © 1996 John Wiley & Sons, Inc.     

Boundary value problems in the theory of thermoporoelasticity for materials with double porosity

In this paper the linear quasi-static theory of thermoelasticity for solids with double porosity is considered. The system of equations of this theory is based on the equilibrium equations for solids with double porosity, conservation of fluid mass, constitutive equations, Darcy's law for materials with double porosity and Fourier's law for heat conduction. The basic internal and external boundary value problems (BVPs) of steady vibrations are formulated. The uniqueness and existence theorems for classical solutions of the above mentioned BVPs are proved by means of the potential method (boundary integral equation method) and the theory of singular integral equations. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Some applications of a galerkin-collocation method for boundary integral equations of the first kind

Here we apply the boundary integral method to several plane interior and exterior boundary value problems from conformal mapping, elasticity and fluid dynamics. These are reduced to equivalent boundary integral equations on the boundary curve which are Fredholm integral equations of the first kind having kernels with logarithmic singularities and defining strongly elliptic pseudodifferential operators of order - 1 which provide certain coercivity properties. The boundary integral equations are approximated by Galerkin's method using B-splines on the boundary curve in connection with an appropriate numerical quadrature, which yields a modified collocation scheme. We present a complete asymptotic error analysis for the fully discretized numerical equations which is based on superapproximation results for Galerkin's method, on consistency estimates and stability properties in connection with the illposedness of the first kind equations in L₂. We also present computational results of several numerical experiments revealing accuracy, efficiency and an amazing asymptotical agreement of the numerical with the theoretical errors. The method is used for computations of conformal mappings, exterior Stokes flows and slow viscous flows past elliptic obstacles.     

Mathematical studies of the solution of Burgers' equations by Adomian decomposition method

In this paper, a novel Adomian decomposition method (ADM) is developed for the solution of Burgers' equation. While high level of this method for differential equations are found in the literature, this work covers most of the necessary details required to apply ADM for partial differential equations. The present ADM has the capability to produce three different types of solutions, namely, explicit exact solution, analytic solution, and semi-analytic solution. In the best cases, when a closed-form solution exists, ADM is able to capture this exact solution, while most of the numerical methods can only provide an approximation solution. The proposed ADM is validated using different test cases dealing with inviscid and viscous Burgers' equations. Satisfactory results are obtained for all test cases, and, particularly, results reported in this paper agree well with those reported by other researchers.     

Different types of spdes in the eyes of girsanov's theorem

We prove Girsanov's theorem for continuous orthogonal martingale measures. We then define space-time SDEs, and use Girsanov's theorem to establish a oneto- one correspondence between solutions of two space-time SDEs differing only by a drift coefficient. For such stochastic equations, we give necessary conditions under which the laws of their solutions are absolutely continuous with respect to each other. Using Girsanov's theorem again, we prove additional existence and uniqueness results for space-time SDEs. The same one-to-one correspondence and absolute continuity theorems are also proved for the stochastic heat and wave equations     

Existence and Numerical Solutions of Nonlinear Quadratic Integral Equations Defined on Unbounded Intervals

The first goal of this article is to discuss the existence of solutions of nonlinear quadratic integral equations. These equations are considered in the Banach space L ^p (?₊). The arguments used in the existence proofs are based on Schauder's and Darbo's fixed point theorems. In particular, to apply Schauder's fixed point theorem based method, a special care is devoted to the proof of the L ^p -compactness of the operators associated with our nonlinear quadratic integral equations. The second goal of this work is to study a numerical method for solving nonlinear Volterra integral equations of a fairly general type. Finally, we provide the reader with some examples that illustrate the different results of this work.