Numerical studies of a nonlinear heat equation with square root reaction term

Interest in calculating numerical solutions of a highly nonlinear parabolic partial differential equation with fractional power diffusion and dissipative terms motivated our investigation of a heat equation having a square root nonlinear reaction term. The original equation occurs in the study of plasma behavior in fusion physics. We begin by examining the numerical behavior of the ordinary differential equation obtained by dropping the diffusion term. The results from this simpler case are then used to construct nonstandard finite difference schemes for the partial differential equation. A variety of numerical results are obtained and analyzed, along with a comparison to the numerics of both standard and several nonstandard schemes. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

A method for solving nonlinear spectral problems for a class of systems of differential algebraic equations

A method for calculating eigenvalues of a nonlinear spectral problem for one class of linear differential algebraic equations is proposed under the assumption of an analytical dependence on spectral parameter of the matrices appearing in the system of equations and the matrices determining boundary conditions.     

The King-Werner method for solving nonsymmetric algebraic Riccati equation

For the nonsymmetric algebraic Riccati equation arising from transport theory, we concern about solving its minimal positive solution. In [1], Lu transferred the equation into a vector form and pointed out that the minimal positive solution of the matrix equation could be obtained via computing that of the vector equation. In this paper, we use the King-Werner method to solve the minimal positive solution of the vector equation and give the convergence and error analysis of the method. Numerical tests show that the King-Werner method is feasible to determine the minimal positive solution of the vector equation.     

On the solution of a mixed nonlinear integral equation

In this paper, we consider a mixed nonlinear integral equation of the second kind in position and time. The existence of a unique solution of this equation is discussed and proved. A numerical method is used to obtain a system of Harmmerstein integral equations of the second kind in position. Then the modified Toeplitz matrix method, as a numerical method, is used to obtain a nonlinear algebraic system. Many important theorems related to the existence and uniqueness solution to the produced nonlinear algebraic system are derived. The rate of convergence of the total error is discussed. Finally, numerical examples when the kernel of position takes a logarithmic and Carleman forms, are presented and the error estimate, in each case, is calculated.     

Existence and multiplicity of positive solutions for an elastic beam equation

This paper investigates the boundary value problem for elastic beam equation of the form

Sufficient close-to-necessary conditions for the blowup of solutions to a strongly nonlinear generalized Boussinesq equation

An initial boundary value problem for the generalized Boussinesq equation with allowance for linear dissipation and free electron sources is considered. The strong generalized time-local solvability of the problem is proved. Sufficient conditions are obtained for the blowup of the solution and for time-global solvability. Two-sided estimates of the blowup time are derived.     

A necessary and sufficient algebraic condition for the controllability of a thermoelastic plate equation

In this paper we give a necessary and sufficient algebraic conditionfor the approximate controllability of the following thermoelasticplate equation with Dirichlet boundary conditions w_tt + ²w + = a₁(x)u₁ + ... + a_m(x)u_m, t 0, x , _t – ß – w_t = d₁(x)u₁ + ... + d_m(x)u_m,t 0, x , = w = w = 0, t 0, x , where 0, ß > 0, is a sufficiently regular boundeddomain in R^N, a_i, d_i, L² (; R), the control functions u_i L²(0, t₁; R); i = 1, 2, ..., m. This condition is easy to checkand is given by Rank [P_jBA_jP_jBA²_jP_jB ... A^3_j–1_jP_jB] = 3_j,BU=b₁U₁+...+b_mU_m,b_i=[0, a_i, d_i], A_j=[0, –²_j, 0, 1, 0, –_j, 0, _j, –ß_j]P_j, j1, where _j, S are the eigenvalues of – with Dirichlet boundarycondition and P_j, S are the projections on the correspondingeigenspace.     

Exact multiplicity of solutions to a diffusive logistic equation with harvesting

An Ambrosetti-Prodi type exact multiplicity result is proved for a diffusive logistic equation with harvesting. We show that a modified diffusive logistic mapping has exactly either zero, or one, or two pre-images depending on the harvesting rate. It implies that the original diffusive logistic equation with harvesting has at most two positive steady state solutions.     

Krylov subspace methods of Hessenberg based for algebraic Riccati equation

In this paper, we propose a class of special Krylov subspace methods to solve continuous algebraic Riccati equation (CARE), i.e., the Hessenberg-based methods. The presented approaches can obtain efficiently the solution of algebraic Riccati equation to some extent. The main idea is to apply Kleinman-Newton"s method to transform the process of solving algebraic Riccati equation into Lyapunov equation at every inner iteration. Further, the Hessenberg process of pivoting strategy combined with Petrov-Galerkin condition and minimal norm condition is discussed for solving the Lyapunov equation in detail, then we get two methods, namely global generalized Hessenberg (GHESS) and changing minimal residual methods based on the Hessenberg process (CMRH) for solving CARE, respectively. Numerical experiments illustrate the efficiency of the provided methods.     

Involutive characteristic sets of algebraic partial differential equation systems

This paper presents an algorithm to reduce a nonlinear algebraic partial differential equation system into the involutive characteristic set with respect to an involutive prolongation direction, which covers the existing algorithms based on Riquier method, Thomas method, and Pommaret method. It also provides new algorithms for computing involutive characteristic sets due to the existence of new involutive directions. Experiments show that these new algorithms may be used to significantly reduce the computational steps in Wu-Ritt's characteristic set method for algebraic partial differential equations.     

Construction of iteration functions for the simultaneous computation of the solutions of equations and algebraic systems

We construct iteration functions for the simultaneous computation of the solutions of a system of equations, with local quadratic convergence: they generalize to the multivariate case the well-known Weierstrass function for polynomials, which is expected to be globally convergent except on a zero-measured set of starting points. We clarify these functions using univariate interpolation. Both for polynomials and algebraic systems with real coefficients, we extend the conjecture of global convergence to the research of real roots or solutions.     

A note on the necessary conditions for the algebraic Riccati equation

We derive some useful and easily computable necessary conditionsfor the existence of a positive semi-definite solution to thealgebraic Ricatti equation (ARE). A motivating example is givento highlight the usefulness of the conditions for controllerand observer designs for nonlinear systems. Further, an upperbound on the trace of the solution to the ARE is also derived.     

Conditions for the existence of nonoscillatory solutions of a second-order nonlinear differential equation

For the nonlinear differential equation

Existence and multiplicity results for nonlinear Schrödinger-Poisson equation with general potential

This paper is concerned with the Schrödinger-Poisson equation-\mathrm{\Delta }u+V(x)u+\phi (x)u=f(x,u),x\in {ℝ}^{\mathrm{3}},-\mathrm{\Delta }\phi =u^{\mathrm{2}},\underset{\left|x\right|\to +\infty }{\lim }\phi (x)=\mathrm{0.}Under certain hypotheses on V and a general spectral assumption, the existence and multiplicity of solutions are obtained via variational methods.     

Energy decay of solutions of a degenerate Kirchhoff equation with a weak nonlinear dissipation

In this paper we study decay properties of the solutions to the degenerate Kirchhoff equation with a weak nonlinear dissipative term.     

An efficient multistep iteration scheme for systems of nonlinear algebraic equations associated with integral equations

The aim of this research is to present a new iterative procedure in approximating nonlinear system of algebraic equations with applications in integral equations as well as partial differential equations (PDEs). The presented scheme consists of several steps to reach a high rate of convergence and also an improved index of efficiency. The theoretical parts are furnished, and several computational tests mainly arising from practical problems are given to manifest its applicability.     

A modified Newton method for solving non-symmetric algebraic Riccati equations arising in transport theory

Liang Bao ^{The non-symmetric algebraic Riccati equation arising in transporttheory can be rewritten as a vector equation and the minimalpositive solution of the non-symmetric algebraic Riccati equationcan be obtained by solving the vector equation. In this paper,we apply the modified Newton method to solve the vector equation.Some convergence results are presented. Numerical tests showthat the modified Newton method is feasible and effective, andoutperforms the Newton method.}