Constructive approximation of solution for fourth‐order nonlinear boundary value problems

In this paper, we obtain a sequence of approximate solution converging uniformly to the exact solution of a class of fourth‐order nonlinear boundary value problems. Its exact solution is represented in the form of series in the reproducing kernel space. The n‐term approximation u_n(x) is proved to converge to the exact solution u(x). Moreover, the derivatives of u_n(x) are also convergent to the derivatives of u(x). Copyright © 2008 John Wiley & Sons, Ltd.     

A best approximation for the solution of one‐dimensional variable‐coefficient Burgers' equation

In this article, an iterative method for the approximate solution to one‐dimensional variable‐coefficient Burgers' equation is proposed in the reproducing kernel space W_(3,2). It is proved that the approximation w_n(x,t) converges to the exact solution u(x,t) for any initial function w₀(x,t) ε W_(3,2), and the approximate solution is the best approximation under a complete normal orthogonal system . Moreover the derivatives of w_n(x,t) are also uniformly convergent to the derivatives of u(x,t).© 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

The Instability of Nonmonotonic Wave Solutions of Parabolic Equations

We consider the class of equations u_t=f(u_xx, u_x, u) under the restriction that for all a,b,c. We first consider this equation over the unbounded domain ? ∞ < x < + ∞, and we show that very nearly every bounded nonmonotonic solution of the form u(t, x)=?(x?ct) is unstable to all nonnegative and all nonpositive perturbations. We then extend these results to nonmonotonic plane wave solutions u(t, x, y)=?(x?ct) of u_t = F(u_xx, u_xy, u_x, u_y, u). Finally, we consider the class of equations u_t=f(u_xx, u_x, u) over the bounded domain 0 < x < 1 with the boundary conditions u(t, x)=A at x=0 and u(t, x)=B at x=1, and we find the stability of all steady solutions u(t, x)=?(x).     

An inverse polynomial method for the identification of the leading coefficient in the Sturm–Liouville operator from boundary measurements

An inverse polynomial method of determining the unknown leading coefficient k=k(x) of the linear Sturm–Liouville operator Au=−(k(x)u^′(x))^′+q(x)u(x), x(0,1), is presented. As an additional condition only two measured data at the boundary (x=0,x=1) are used. In absence of a singular point (u^′(x)≠0,u^″(x)≠0,x[0,1]) the inverse problem is classified as a well-conditioned . If there exists at least one singular point, then the inverse problem is classified as moderately ill-conditioned (u^′(x₀)=0,x₀(0,1);u^′(x)≠0,x≠x₀;u^″(x)≠0,x[0,1]) and severely ill-conditioned (u^′(x₀)=u^″(x₀)=0,x₀(0,1);u^′(x)≠0,u^″(x)≠0,x≠x₀). For each of the cases direct problem solution is approximated by corresponding polynomials and the inverse problem is reformulated as a Cauchy problem for to the first order differential equation with respect the unknown function k=k(x). An approximate analytical solution of the each Cauchy problems are derived in explicit form. Numerical simulations all the above cases are given for noise free and noisy data. An accuracy of the presented approach is demonstrated on numerical test solutions.     

A Threshold Result for a Non-local Parabolic Equation

In this paper, we consider the Cauchy problem: (ECP) u_t−Δu+p(x)u=u(x,t)∫u²(y,t)/∣x−y∣dy; x∈ℝ³, t>0, u(x, 0)=u₀(x)⩾0 x∈ℝ³, (0.2) The stationary problem for (ECP) is the famous Choquard–Pekar problem, and it has a unique positive solution ū(x) as long as p(x) is radial, continuous in ℝ³, p(x)⩾ā>0, and lim_{∣x∣→∞}p(x)=p¯>0. In this paper, we prove that if the initial data 0⩽u₀(x)⩽(≢)ū(x), then the corresponding solution u(x, t) exists globally and it tends to the zero steady-state solution as t→∞, if u₀(x)⩾(≢)ū(x), then the solution u(x,t) blows up in finite time. © 1997 B. G. Teubner Stuttgart–John Wiley & Sons Ltd.     

Existence and Uniqueness of Strong Solutions for a Class of Quasi‐Linear Hyperbolic Equations with Order Degeneration

Let k(y) > 0, 𝓁(y) > 0 for y > 0, k(0) = 𝓁(0) = 0 and lim_{y → 0}k(y)/𝓁(y) exists; then the equation L(u) ≔ k(y)u_xx – ∂_y(𝓁(y)u_y) + a(x, y)u_x = f(x, y, u) is strictly hyperbolic for y > 0 and its order degenerates on the line y = 0. Consider the boundary value problem Lu = f(x, y, u) in G, u|_AC = 0, where G is a simply connected domain in ℝ² with piecewise smooth boundary ∂G = AB∪AC∪BC; AB = {(x, 0) : 0 ≤ x ≤ 1}, AC : x = F(y) = ∫^y₀(k(t)/𝓁(t))^1/2dt and BC : x = 1 – F(y) are characteristic curves. Existence of generalized solution is obtained by a finite element method, provided f(x, y, u) satisfies Carathéodory condition and |f(x, y, u)| ≤ Q(x, y) + b|u| with Q ∈ L²(G), b = const > 0. It is shown also that each generalized solution is a strong solution, and that fact is used to prove uniqueness under the additional assumption |f(x, y, u₁) – f(x, y, u₂| ≤ C|u₁ – u₂|, where C = const > 0.     

Multiple solutions of sublinear Lane-Emden elliptic equations

We study the sublinear elliptic equation, −Δ u = |u|^psgn u + f(x,u) in the bounded domain Ω under the zero Dirichlet boundary condition. We suppose that 0 < p < 1 and |f(x,u)| is small enough near u = 0 and do not suppose that f(x,u) is odd on u. Then we prove that this problem has infinitely many solutions. Supported in part by the Grant-in-Aid for Scientific Research (C) (No. 16540179), Ministry of Education, Culture, Sports, Science and Technology.     

The Generalized Riemann Problem for First Order Quasilinear Hyperbolic Systems of Conservation Laws II

In this paper, we are concerned with the existence and uniqueness of the local solution to the generalized Riemann problem for first order quasi-linear hyperbolic systems of conservation laws in the presence of the shock wave with large amplitude and the centered wave. Apart from some exceptions, we prove the problem admits a unique piecewise smooth solution u=u(t,x), and this solution has a structure similar to the similarity solution u=u(x/t) of the corresponding Riemann problem in the neighborhood of the origin, provided that the coefficients of the system and the initial conditions are sufficiently smooth. The application of our results in rich system is also given.     

A global Newton method for the zeros of cylinder functions

The zeros of cylinder functions C _u (x)=cos α, J _u (x) - sin α, Y _u(x) coincide with those of the ratios H _u (x)=C _u (x)/C _u-1 (x) except, perhaps, at x = 0. We show monotonicity properties of H u(x) and f _u (x) = x ^2v-1 H u(x) and their derivatives for x > 0. We then build a Newton-Raphson iterative method based on the monotonic function f u(x) which is shown to be convergent, for any real values of u and α and any starting value x ₀ > 0, to an sth positive root c ,s of C _u (x) = 0, s being such that c _,s and x0 belong to the same interval (c _u-1 ,s', c _{u -1} ,s'+1]. We also show applications of the method. In particular, taking advantage of the fact that the ratio H _u (x) for first kind Bessel functions J u(x) can be evaluated by using a continued fraction, a very simple algorithm is built; it becomes especially efficient for low values of u and s and it allows the evaluation of the real zeros for arbitrary orders u, positive or negative. This revised version was published online in August 2006 with corrections to the Cover Date.     

On the existence of shock fronts for a Burgers equation with a singular source term

In this paper we consider the Cauchy problem for the equation ∂u/∂t + u ∂u/∂x + u/x = 0 for x > 0, t ⩾ 0, with u(x, 0) = u⁰₋(x) for x < x₀, u(x, 0) = u⁰₊(x) for x > x₀, u⁰₋(x₀) > u⁰₊(x₀). Following the ideas of Majda, 1984 and Lax, 1973, we construct, for smooth u⁰₋ and u⁰₊, a global shock front weak solution u(x, t) = u₋(x, t) for x < ϕ(t), u(x, t) = u₊(x, t) for x > ϕ(t), where u₋ and u₊ are the strong solutions corresponding (respectively) to u⁰₋ and u⁰₊ and the curve t → ϕ(t) is defined by dϕ/dt (t) = 1/2[u₋(ϕ(t), t) + u₊(ϕ(t), t)], t ⩾ 0 and ϕ(0) = x₀. © 1998 B. G. Teubner Stuttgart—John Wiley & Sons, Ltd.     

Comparison results for nonlinear elliptic equations with lower–order terms

We consider a solution u of the homogeneous Dirichlet problem for a class of nonlinear elliptic equations in the form A(u) + g(x, u) = f, where the principal term is a Leray–Lions operator defined on and g(x, u) is a term having the same sign as u and satisfying suitable growth assumptions. We prove that the rearrangement of u can be estimated by the solution of a problem whose data are radially symmetric.     

Initial boundary value problem for a class of fourth-order wave equation with viscous damping term

In this article we study the initial boundary value problem for a class of fourth-order nonlinear wave equation with viscous damping term u _tt ???αu _xxt ?+?u _xxxx ?=?f(u _x ) _x . By argument related to the potential well-convexity method, we prove the global existence and nonexistence of the solution. Further, we give some sharp conditions for global existence and nonexistence of the solution. This generalizes the results obtained in Chen and Lu [G. Chen and B. Lu, The initial-boundary value problems for a class of nonlinear wave equations with damping term, J. Math. Anal. Appl. 351 (2009), pp. 1–15].     

The Ambrosetti–Prodi Problem for the p-Laplace Operator

In this work we study the existence of a solution for the problem ? Δ _p u = f(u) + tΦ(x) + h(x), with homogeneous Dirichlet boundary conditions. Here the nonlinear term f(u) is a so-called jumping nonlinearity. In the proofs we use topological arguments and the sub-supersolutions method, together with comparison principles for the p-Laplacian.     

Parameter dependence for a class of ordinary differential equations with measurable right‐hand side

The paper studies differential equations of the form u′(x) = f(x, u(x), λ(x)), u(x₀) = u₀, where the right‐hand side is merely measurable in x. In particular sufficient conditions for the continuous and the differentiable dependence of solution u on the data and on the parameter λ are stated. © 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim     

Solution of mixed problems with boundary elastic-force control for the telegraph equation

In terms of a finite-energy generalized solution of the telegraph equation, for any time interval T, we consider the problem on the boundary elastic-force control u _x (0, t) = μ(t) at the endpoint x = 0 for the process described by the Klein-Gordon-Fock equation under the condition that the other endpoint x = l is either fixed, or free, or is controlled by an elastic force. For any time interval T, we obtain the solution u(x, t) in closed form.     

Periodic and Homoclinic Solutions of Extended Fisher–Kolmogorov Equations

In this paper we study the existence of periodic solutions of the fourth-order equations u^iv − pu″ − a(x)u + b(x)u³ = 0 and u^iv − pu″ + a(x)u − b(x)u³ = 0, where p is a positive constant, and a(x) and b(x) are continuous positive 2L-periodic functions. The boundary value problems (P₁) and (P₂) for these equations are considered respectively with the boundary conditions u(0) = u(L) = u″(0) = u″(L) = 0. Existence of nontrivial solutions for (P₁) is proved using a minimization theorem and a multiplicity result using Clark's theorem. Existence of nontrivial solutions for (P₂) is proved using the symmetric mountain-pass theorem. We study also the homoclinic solutions for the fourth-order equation u^iv + pu″ + a(x)u − b(x)u² − c(x)u³ = 0, where p is a constant, and a(x), b(x), and c(x) are periodic functions. The mountain-pass theorem of Brezis and Nirenberg and concentration-compactness arguments are used.     

Using a Mixed Integer Quadratic Programming Solver for the Unconstrained Quadratic 0-1 Problem

In this paper, we consider problem (P) of minimizing a quadratic function q(x)=x ^t Qx+c ^t x of binary variables. Our main idea is to use the recent Mixed Integer Quadratic Programming (MIQP) solvers. But, for this, we have to first convexify the objective function q(x). A classical trick is to raise up the diagonal entries of Q by a vector u until (Q+diag(u)) is positive semidefinite. Then, using the fact that x _i ²=x _i, we can obtain an equivalent convex objective function, which can then be handled by an MIQP solver. Hence, computing a suitable vector u constitutes a preprocessing phase in this exact solution method. We devise two different preprocessing methods. The first one is straightforward and consists in computing the smallest eigenvalue of Q. In the second method, vector u is obtained once a classical SDP relaxation of (P) is solved. We carry out computational tests using the generator of (Pardalos and Rodgers, 1990) and we compare our two solution methods to several other exact solution methods. Furthermore, we report computational results for the max-cut problem.     

The radial basis functions method for identifying an unknown parameter in a parabolic equation with overspecified data

Parabolic partial differential equations with overspecified data play a crucial role in applied mathematics and engineering, as they appear in various engineering models. In this work, the radial basis functions method is used for finding an unknown parameter p(t) in the inverse linear parabolic partial differential equation u_t = u_xx + p(t)u + φ, in [0,1] × (0,T], where u is unknown while the initial condition and boundary conditions are given. Also an additional condition ∫₀¹k(x)u(x,t)dx = E(t), 0 ≤ t ≤ T, for known functions E(t), k(x), is given as the integral overspecification over the spatial domain. The main approach is using the radial basis functions method. In this technique the exact solution is found without any mesh generation on the domain of the problem. We also discuss on the case that the overspecified condition is in the form ∫₀^s(t) u(x,t)dx = E(t), 0 < t ≤ T, 0 < s(t) < 1, where s and E are known functions. Some illustrative examples are presented to show efficiency of the proposed method. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

Analysis for the identification of an unknown diffusion coefficient via semigroup approach

This paper presents a semigroup approach for the mathematical analysis of the inverse coefficient problems of identifying the unknown coefficient k(u_x) in the inhomogenenous quasi‐linear parabolic equation u_t(x, t)=(k(u_x)u_x(x, t))_x +F(u), with the Dirichlet boundary conditions u(0, t)=ψ₀, u(1, t)=ψ₁ and source function F(u). The main purpose of this paper is to investigate the distinguishability of the input–output mappings Φ[·]:??→C¹[0, T], Ψ[·]:??→C¹[0, T] via semigroup theory. Copyright © 2009 John Wiley & Sons, Ltd.