On the traveling wave solutions for a parabolic–hyperbolic system in linear thermoelasticity

In this paper, the parabolic–hyperbolic system of linear thermoelasticity with variable coefficients is transformed into a system of two coupled equations. We discuss first the conditions which govern this separation in the case of a system of two coupled equations for which a general result on the separability is formulated. It is then shown that the explicit traveling wave solutions are obtained in the exact form.     

Conditional linearizability criteria for a system of third-order ordinary differential equations

We provide linearizability criteria for a class of systems of two third-order ordinary differential equations that is cubically nonlinear in the first derivative, by differentiating a system of second-order quadratically nonlinear ordinary differential equations and using the original system to replace the second derivatives. The procedure developed splits into two cases: those for which the coefficients are constant and those for which they are variables. Both cases are discussed and examples given.     

筛选法解非线性方程组

本文给出了一种新方法解非线性方程组，也是筛选法的一个推广方程组被分成两部分，一部分被当作约束条件，另一部分的最小二乘被当作目标函数．本质上，两种不同方法被用于解同一非线性方程组．     

Global Solution of an Initial-Value Problem for Two-Dimensional Compressible Euler Equations

This paper is concerned with the existence of global continuous solutions of the expansion of a wedge of gas into a vacuum for compressible Euler equations. By hodograph transformation, we first prove that the flow is governed by a partial differential equation of second order, which is further reduced to a system of two nonhomogeneous linearly degenerate equations in the phase space under an irrotationality condition. Then this conclusion is applied to solving the problem that a wedge of gas expands into a vacuum, which is actually a Goursat-type problem for these two equations in the supersonic domain.     

The inverse problem for Lagrangian systems with certain non-conservative forces

We discuss two generalizations of the inverse problem of the calculus of variations, one in which a given mechanical system can be brought into the form of Lagrangian equations with non-conservative forces of a generalized Rayleigh dissipation type, the other leading to Lagrangian equations with so-called gyroscopic forces. Our approach focusses primarily on obtaining coordinate-free conditions for the existence of a suitable non-singular multiplier matrix, which will lead to an equivalent representation of a given system of second-order equations as one of these Lagrangian systems with non-conservative forces.     

可行逐步二次规划方法解非线性方程组

We transform the system of nonlinear equations into a nonlinear programming problem, which is attacked by feasible sequential quadratic programming(FSQP) method.We do not employ standard least square approach.We divide the equations into two groups. One group, which contains the equations with zero residual,is treated as equality constraints. The square of other equations is regarded as objective function. Two groups are updated in every step.Therefore, the subproblem is updated at every step, which avoids the difficulty that it is required to lie in feasible region for FSQP.     

Well-posedness and exponential stability of Maxwell-like systems coupled with strongly absorbing layers

This work deals with a family of PDEs which are derived from the Maxwell system set into an anisotropic medium. The system consists of four coupled linear equations. The first two correspond to the Maxwell system perturbed by zero-order operators which are represented by diagonal tensors with compact support. The last two equations are ODEs. Each system is written in a bounded domain and its boundary is modelled by a Silver–Müller condition. Firstly we establish that each problem is well-posed, providing the tensors are positive with bounded terms in Ω. Secondly we address the question of the long-time stability of each model. We prove that there exists a functional of energy which is exponentially decreasing if the domain is strictly star-shaped.     

Blowup with vorticity control for a 2D model of the Boussinesq equations

We propose a system of equations with nonlocal flux in two space dimensions which is closely modeled after the 2D Boussinesq equations in a hyperbolic flow scenario. Our equations involve a vorticity stretching term and a non-local Biot-Savart law and provide insight into the underlying intrinsic mechanisms of singularity formation. We prove stable, controlled finite time blowup involving upper and lower bounds on the vorticity up to the time of blowup for a wide class of initial data.     

Hypergeometric functions of two variables

Summary The systematic investigation of contour integrals satisfying the system of partial differential equations associated with Appell's hypergeometric functionF ₁ leads to new solutions of that system. Fundamental sets of solutions are given for the vicinity of all singular points of the system of partial differential equations. The transformation theory of the solutions reveals connections between the system under consideration and other hypergeometric systems of partial differential equations. Presently it is discovered that any hypergeometric system of partial differential equations of the second order (with two independent variables) which has only three linearly independent solutions can be transformed into the system ofF ₁ or into a particular or limiting case of this system. There are also other hypergeometric systems (with four linearly independent solutions) the integration of which can be reduced to the integration of the system ofF ₁.     

An iterative splitting approach for linear integro-differential equations

The motivation for this paper is to solve a model based on the dynamics of electrons in a plasma using a simplified Boltzmann equation. Such problems have arisen in active plasma resonance spectroscopy, which is used for plasma diagnostic techniques; see Braithwaite and Franklin (2009) [1]. We propose a modified iterative splitting approach to solve the Boltzmann equations as a system of integro-differential equations. To enable solution by fast and iterative computations, we first transform the integro-differential equations into second order differential equations. Second, we split each second order differential equations into two first order differential equations via a splitting approach. We carry out an error analysis of the higher order iterative approach. Numerical experiments with a simplified Boltzmann equation will be discussed, along with the benefits of computing with this splitting approach.     

Analysis of a nutrient-phytoplankton model in the presence of viral infection

In this paper, a system of reaction-diffusion equations arising in a nutrient-phytoplankton populations is investigated. The equations model a situation in which phytoplankton population is divided into two groups, namely susceptible phytoplankton and infected phytoplankton. A number of existence and non-existence results about the non-constant steady states of a reaction diffusion system are given. If the diffusion coefficient of the infected phytoplankton is treated as bifurcation parameter, non-constant positive steady-state solutions may bifurcate from the constant steady-state solution under some conditions.     

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.     

Complex Lie symmetries for scalar second-order ordinary differential equations

A scalar complex ordinary differential equation can be considered as two coupled real partial differential equations, along with the constraint of the Cauchy–Riemann equations, which constitute a system of four equations for two unknown real functions of two real variables. It is shown that the resulting system possesses those real Lie symmetries that are obtained by splitting each complex Lie symmetry of the given complex ordinary differential equation. Further, if we restrict the complex function to be of a single real variable, then the complex ordinary differential equation yields a coupled system of two ordinary differential equations and their invariance can be obtained in a non-trivial way from the invariance of the restricted complex differential equation. Also, the use of a complex Lie symmetry reduces the order of the complex ordinary differential equation (restricted complex ordinary differential equation) by one, which in turn yields a reduction in the order by one of the system of partial differential equations (system of ordinary differential equations). In this paper, for simplicity, we investigate the case of scalar second-order ordinary differential equations. As a consequence, we obtain an extension of the Lie table for second-order equations with two symmetries.     

The use of splitting of the equations of motion in the numerical solution of problems of the dynamics of elastic systems

Linear elastic systems with a finite number of degrees of freedom, the initial equations of motion of which are constructed using the finite element method or other discretization methods, are considered. Since, in applied dynamics problems, the motions are usually investigated in a frequency range with an upper bound, the degrees of freedom of the initial system of equations are split into dynamic and quasi-dynamic degrees. Finally, the initial system of equations is split into a small number of differential equations for the dynamic degrees of freedom and into a system of algebraic equations for determining the quasi-static displacements, represented in the form of a matrix series. The number of terms of the series taken into account depends on the accuracy required.     

Investigation of variational problems by direct methods

A direct method is proposed for solving variational problems in which an extremal is represented by an infinite series in terms of a complete system of basis functions. Taking into account the boundary conditions gives all the necessary conditions of the classical calculus of variations, that is, the Euler-Lagrange equations, transversality conditions, Erdmann-Weierstrass conditions, etc. The penalty function method reduces conditional extremum problems to variational ones in which the isoperimetric conditions described by constraint equations are taken into account by Lagrangian multipliers. The direct method proposed is applied to functionals depending on functions of one or two variables.     

Shear and Flexural Buckling Modes of a Spherical Sandwich Shell in a Centrosymmetric Temperature Field Inhomogeneous across the Thickness

Problems on buckling modes (BMs) are considered for a spherical sandwich shell with thin isotropic external layers and a transversely soft core of arbitrary thickness in a centrosymmetric temperature field inhomogeneous across the shell thickness. For their statement, the two-dimensional equations of the theory of moderate bending of thin Kirchhoff–Love shells are used for the external layers, with regard for their interaction with the core; for the core, maximum simplified geometrically nonlinear equations of thermoelasticity theory, in which a minimum number of nonlinear summands is retained to correctly describe its pure shear BM, are utilized. An exact analytical solution to the problem on initial centrosymmetric deformation of the shell is found, assuming that the temperature increments in the external layers are constant across their thickness. It is shown that the three-dimensional equations for the core, linearized in the neighborhood of the solution, can be integrated along the radial coordinate and reduced to two two-dimensional differential equations, which supplement the six equations that describe the neutral equilibrium of the external layers. It is established that the system of eight differential equations of stability, upon introduction of new unknowns in the form of scalar and vortical potentials, splits into two uncoupled sets of equations. The first of them has two kinds of solutions, by which the pure shear BM is described at an identical value of the parameter of critical temperature. The second system describes a mixed flexural BM, whose realization, at definite combinations of determining parameters of the shell and over wide ranges of their variation, is possible for critical parameters of temperature by orders of magnitude exceeding the similar parameter of shear BM.     

Tracking within a time interval on the basis of data supplied by finite observers

We solve the tracking control problem, in which one should bring a trajectory of a system of linear ordinary differential equations into a neighborhood of a trajectory of another system within a given time interval. After getting into this neighborhood, one should keep the trajectory of the first subsystem in it for a time interval of given duration. For the control synthesis, we use incomplete and imprecise information on the online deviation of one trajectory from the other, which is obtained in real time from linear equations of observation. We consider distinct structures of observers, which substantially affect the solution of control problems for such systems. The equations of dynamics and admissible measurements contain uncertainty for which one knows only some hard pointwise constraints. To solve the main problem, we use an approach that can be reduced to the construction of auxiliary information sets and weakly invariant sets with a subsequent “aiming” of one set at a tube. We suggest an efficient method for an approximate solution on the basis of ellipsoidal calculus techniques. The results of the algorithm operation are illustrated by an example of the solution of a tracking control problem for two fourth-order subsystems.     

On the approximate solution of some two‐dimensional singular integral equations

An approximation method for a wide class of two‐dimensional integral equations is proposed. The method is based on using a special function system. Orthonormality and good interaction with fundamental integral operators arising in partial differential equations are remarkable properties of this system. In addition, all the basis elements can easily be calculated by recurrence relations. Taking into account these properties we construct a numerical algorithm which does not require additional effort (such as quadrature) to compute the values of the fundamental operators on the basis elements. Copyright © 2001 John Wiley & Sons, Ltd.     

Two‐dimensional shifted Legendre polynomial collocation method for electromagnetic waves in dielectric media via almost operational matrices

In this paper, a numerical solution of fractional partial differential equations (FPDEs) for electromagnetic waves in dielectric media will be discussed. For the solution of FPDEs, we developed a numerical collocation method using an algorithm based on two‐dimensional shifted Legendre polynomials approximation, which is proposed for electromagnetic waves in dielectric media. By implementing the partial Riemann–Liouville fractional derivative operators, two‐dimensional shifted Legendre polynomials approximation and its operational matrix along with collocation method are used to convert FPDEs first into weakly singular fractional partial integro‐differential equations and then converted weakly singular fractional partial integro‐differential equations into system of algebraic equation. Some results concerning the convergence analysis and error analysis are obtained. Illustrative examples are included to demonstrate the validity and applicability of the technique. Copyright © 2017 John Wiley & Sons, Ltd.     

A semi-analytical direct optimal control solution for strongly excited and dissipative Hamiltonian systems

A semi-analytical direct optimal control solution for strongly excited and dissipative Hamiltonian systems is proposed based on the extended Hamiltonian principle, the Hamilton-Jacobi-Bellman (HJB) equation and its variational integral equation, and the finite time element approximation. The differential extended Hamiltonian equations for structural vibration systems are replaced by the variational integral equation, which can preserve intrinsic system structure. The optimal control law dependent on the value function is determined by the HJB equation so as to satisfy the overall optimality principle. The partial differential equation for the value function is converted into the integral equation with variational weighting. Then the successive solution of optimal control with system state is designed. The two variational integral equations are applied to sequential time elements and transformed into the algebraic equations by using the finite time element approximation. The direct optimal control on each time element is obtained respectively by solving the algebraic equations, which is unconstrained by the system state observed. The proposed control algorithm is applicable to linear and nonlinear systems with the quadratic performance index, and takes into account the effects of external excitations measured on control. Numerical examples are given to illustrate the optimal control effectiveness.