Blowup of solutions for the “bad” Boussinesq-type equation

The paper studies the blowup of solutions to the initial boundary value problem for the “bad” Boussinesq-type equation u_tt−u_xx−bu_xxxx=σ(u)_xx, where b>0 is a real number and σ(s) is a given nonlinear function. By virtue of the energy method and the Fourier transform method, respectively, it proves that under certain assumptions on σ(s) and initial data, the generalized solutions of the above-mentioned problem blow up in finite time. And a few examples are shown, especially for the “bad” Boussinesq equation, two examples of blowup of solutions are obtained numerically.     

On the damped nonlinear evolution equation wtt = σ(w)xx − ywt

Initial boundary value problems for the damped nonlinear wave equation w_tt = σ(w)_xx ? yw_t arise in several areas of applied mathematics and, in particular, in studies of shearing flow in a nonlinear viscoelastic fluid; the problems of global existence and nonexistence of smooth solutions have been extensively studied in the strictly hyperbolic case σ′(δ) ? ε > 0, ?δ?R¹ as well as in the case where σ′(0) > 0 and the initial data are chosen so small that σ′(w) > 0 for as long as a smooth solution w(x, t) exists. In this paper the global nonexistence problem is studied for the cases σ′(0) = 0 and σ′(0) > 0 but σ′(δ) < 0 for ¦δ¦ sufficiently large and growth estimates which are valid on the maximal interval of existence of a sufficiently smooth solution are derived.     

A study on 3-D stress distributions in the bi-adhesively bonded T-joints

In this study, the normal (σ_x, σ_y, σ_z) and shear stress (τ_xy, τ_yz, τ_zx) distributions occurring in a bi-adhesively bonded T-joint with was investigated via a non-linear three dimensional finite element analysis. For this purpose, first of all, using 2024-T3 aluminum alloy as the adherend and the support, a two-part paste (DP 460) and a film type (SBT 9244) as adhesive, two different types of T-joint samples (single-adhesively bonded T-joint and bi-adhesively bonded T-joint) were produced for experimental studies. After experimental studies on the three different T-joint types were conducted, stress analyses in the T-joints were performed with a three-dimensional finite element analysis by considering the geometrical non-linearity and the material non-linearities of the adhesive (DP460 and SBT9244) and adherend (AA2024-T3). Finally, for a given adherend, the lower the stiffness of the adhesive used in the overlap, the lower the stress concentration, leading to potentially higher joint strength. The use of relatively low stiffness adhesives at the ends of the overlap in a bi-adhesive can decrease the stress concentration and, therefore, potentially lead to higher joint strength.     

Control parameter estimation in a semi-linear parabolic inverse problem using a high accurate method

Parabolic inverse problems have an important role in many branches of science and technology. The aim of this research work is to solve these classes of equations using a high order compact finite difference scheme. We consider the following inverse problem for finding u(x, t) and p(t) governed by u_t = u_xx + p(t)u + φ(x, t) with an over specified condition inside the domain. Spatial derivatives are approximated using central difference scheme. The time advancement of the simulation is performed using a “third order compact Runge-Kutta method”. The convergence orders for the approximation of both u and p are of o(k³ + h²) which improves the results obtained in the literature. An exact test case is used to evaluate the validity of our numerical analysis. We found that the accuracy of the results is better than that of previous works in the literature.     

Analysis of a semigroup approach in the inverse problem of identifying an unknown parameters

This article presents a semigroup approach to the mathematical analysis of the inverse parameter problems of identifying the unknown parameters p(t) and q in the linear parabolic equation u_t(x, t)  = u_xx + qu_x(x, t) + p(t)u(x, t), with Dirichlet boundary conditions u(0, t) = ψ₀, u(1, t) = ψ₁. The main purpose of this paper is to investigate the distinguishability of the input-output mapping Φ[·]:P→H^1,2[0,T], via semigroup theory. In this paper, it is shown that if the nullspace of the semigroup T(t) consists of only zero function, then the input-output mapping Φ[·] has the distinguishability property. It is also shown that the types of the boundary conditions and the region on which the problem is defined play an important role in the distinguishability property of the mapping. Moreover, under the light of the measured output data u_x(0, t) = f(t) the unknown parameter p(t) at (x, t) = (0, 0) and the unknown coefficient q are determined via the input data. Furthermore, it is shown that measured output data f(t) can be determined analytically by an integral representation. Hence the input-output mapping Φ[·]:P→H^1,2[0,T] is given explicitly interms of the semigroup.     

A finite difference solution to a two-dimensional parabolic inverse problem

In this paper, we consider the problem of finding u = u(x, y, t) and p = p(t) which satisfy u_t = u_xx + u_yy + p(t)u + ? in R × [0, T], u(x, y, 0) = f(x, y), (x, y) ∈ R = [0, 1] × [0, 1], u is known on the boundary of R and u(x^∗, y^∗, t) = E(t), 0 < t ? T, where E(t) is known and (x^∗, y^∗) is a given point of R. Through a function transformation, the nonlinear two-dimensional diffusion problem is transformed into a linear problem, and a backward Euler scheme is constructed. It is proved by the maximum principle that the scheme is uniquely solvable, unconditionally stable and convergent in L_∞ norm. The convergence orders of u and p are of O(τ + h²). The impact of initial data errors on the numerical solution is also considered. Numerical experiments are presented to illustrate the validity of the theoretical results.     

Numerical similarity reductions of the (1+3)-dimensional Burgers equation

We consider the (1+3)-dimensional Burgers equation u_t = u_xx + u_yy + u_zz + uu_x which has considerable interest in mathematical physics. Lie symmetries are used to reduce it to certain ordinary differential equations. We employ numerical methods to solve a number of these ordinary differential equations.     

Stable periodic traveling waves for a predator-prey model with non-constant death rate and delay

In this paper we will consider a predator-prey model with a non-constant death rate and distributed delay, described by a partial integro-differential system. The main goal of this work is to prove that the partial integro-differential system has periodic orbitally asymptotically stable solutions in the form of periodic traveling waves; i.e. N(x, t) = N(σt − μ · x), P(x, t) = P(σt − μ · x), where σ > 0 is the angular frequency and μ is the vector number of the plane wave, which propagates in the direction of the vector μ with speed c = σ/∥μ∥; and N(x, t) and P(x, t) are the spatial population densities of the prey and the predator species, respectively. In order to achieve our goal we will use singular perturbation’s techniques.     

Blowup of solutions for improved Boussinesq type equation

The paper studies the existence and uniqueness of local solutions and the blowup of solutions to the initial boundary value problem for improved Boussinesq type equation u_tt−u_xx−u_xxtt=σ(u)_xx. By a Galerkin approximation scheme combined with the continuation of solutions step by step and the Fourier transform method, it proves that under rather mild conditions on initial data, the above-mentioned problem admits a unique generalized solution u∈W^2,∞([0,T];H²(0,1)) as long as . In particular, when σ(s)=as^p, where a≠0 is a real number and p>1 is an integer, specially a<0 if p is an odd number, the solution blows up in finite time. Moreover, two examples of blowup are obtained numerically.     

Critical Exponents of Quasilinear Parabolic Equations

In this paper we study the critical exponents of the Cauchy problem in Rⁿ of the quasilinear singular parabolic equations: u_t = div(|∇u|^m − 1∇u) + t^s|x|^σu^p, with non-negative initial data. Here s ≥ 0, (n − 1)/(n + 1) < m < 1, p > 1 and σ > n(1 − m) − (1 + m + 2s). We prove that p_c ≡ m + (1 + m + 2s + σ)/n > 1 is the critical exponent. That is, if 1 < p ≤ p_c then every non-trivial solution blows up in finite time, but for p > p_c, a small positive global solution exists.     

Continuous spectrum, point spectrum and residual spectrum of operator matrices

Let M_C denote a 2 × 2 upper triangular operator matrix of the form , which is acting on the sum of Banach spaces X⊕Y or Hilbert spaces H⊕K. In this paper, the sets and ?_C∈B(K,H)σ_r(M_C) are, respectively, characterized completely, where σ_c(·) denotes the continuous spectrum, σ_p(·) denotes the point spectrum and σ_r(·) denotes the residual spectrum. Moreover, some corresponding counterexamples are given.     

A new off-step high order approximation for the solution of three-space dimensional nonlinear wave equations

In this paper, we propose a new high accuracy numerical method of O(k² + k²h² + h⁴) based on off-step discretization for the solution of 3-space dimensional non-linear wave equation of the form u_tt = A(x,y,z,t)u_xx + B(x,y,z,t)u_yy + C(x,y,z,t)u_zz + g(x,y,z,t,u,u_x,u_y,u_z,u_t), 0 < x,y,z < 1,t > 0 subject to given appropriate initial and Dirichlet boundary conditions, where k > 0 and h > 0 are mesh sizes in time and space directions respectively. We use only seven evaluations of the function g as compared to nine evaluations of the same function discussed in  and . We describe the derivation procedure in details of the algorithm. The proposed numerical algorithm is directly applicable to wave equation in polar coordinates and we do not require any fictitious points to discretize the differential equation. The proposed method when applied to a telegraphic equation is also shown to be unconditionally stable. Comparative numerical results are provided to justify the usefulness of the proposed method.     

Weak separation axioms and weak covering properties

We identify some remnants of normality and call them rudimentary normality, generalize the concept of submetacompact spaces to that of a weakly subparacompact space and that of a weakly^? subparacompact space, and make a simultaneous generalization of collectionwise normality and screenability with the introduction of what is to be called collectionwise σ-normality. With these weak properties, we show that,1) on weakly subparacompact spaces, countable compactness = compactness, ω₁-compactness = Lindelöfness;2) on weakly subparacompact Hausdorff spaces with rudimentary normality, regularity = normality = countable paracompactness; and3) on weakly subparacompact regular T₁-spaces with rudimentary normality, collectionwise σ-normality = screenability = collectionwise normality = paracompactness.The famous Normal Moore Space Conjecture is thus given an even more striking appearance and Worrell and Wicke?s factorization of paracompactness (over Hausdorff spaces) along with Krajewski?s are combined and strengthened. The methodology extends itself to the factorization of paracompactness on locally compact, locally connected spaces in the manner of Gruenhage and on locally compact spaces in that of Tall, and to the factorization of subparacompactness and metacompactness in the genre of Katuta, Chaber, Junnila and Price and Smith and that of Boone, improving all of them.     

Peakon, loop and solitary travelling wave solutions for the general modified CH-DP equation

Travelling wave solutions for the general modified CH-DP equation u_t − u_xxt + αu²u_x − βu_xu_xx = uu_xxx are developed. By using the dynamical system method, a peakon and a dark soliton are found to coexist for the same wave speed. Exact explicit blow-up solutions are given. By using numerical simulation, a loop solution for a special case is discussed.     

One-point pseudospectral collocation for the one-dimensional Bratu equation

The one-dimensional planar Bratu problem is u_xx + λ exp(u) = 0 subject to u(±1) = 0. Because there is an analytical solution, this problem has been widely used to test numerical and perturbative schemes. We show that over the entire lower branch, and most of the upper branch, the solution is well approximated by a parabola, u(x) ≈ u₀ (1 − x²) where u₀ is determined by collocation at a single point x = ξ. The collocation equation can be solved explicitly in terms of the Lambert W-function as u(0) ≈ −W(−λ(1 − ξ²)/2)/(1 − ξ²) where both real-valued branches of the W-function yield good approximations to the two branches of the Bratu function. We carefully analyze the consequences of the choice of ξ. We also analyze the rate of convergence of a series of even Chebyshev polynomials which extends the one-point approximation to arbitrary accuracy. The Bratu function is so smooth that it is actually poor for comparing methods because even a bad, inefficient algorithm is successful. It is, however, a solution so smooth that a numerical scheme (the collocation or pseudospectral method) yields an explicit, analytical approximation. We also fill some gaps in theory of the Bratu equation. We prove that the general solution can be written in terms of a single, parameter-free β(x) without knowledge of the explicit solution. The analytical solution can only be evaluated by solving a transcendental eigenrelation whose solution is not known explicitly. We give three overlapping perturbative approximations to the eigenrelation, allowing the analytical solution to be easily evaluated throughout the entire parameter space.     

Length problems for automorphisms of modules over local rings

Let R be a local ring and M a free module of a finite rank over R. An element τ ∈ Aut_RM is said to be simple if τ ≠ 1 fixes a hyperplane of M.We shall show that for any σ ∈ Aut_RM there exist a basis X for M and ρ ∈ Aut_RM such that ρ acts as a permutation on X and ρ⁻¹σ is a product of m or less than m simple elements in Aut_RM, where m is the order of the invariant factors of σ modulo the maximal ideal of R.Also we shall investigate the problem treated by E.W. Ellers and H. Ishibashi [Factorizations of transformations over a valuation ring, Linear Algebra Appl. 85 (1987) 17-27], in which they showed that σ is a product of simple elements and gave an upper bound of the smallest number of such factors of σ, whereas in the present paper we will give lower bounds of σ in case that R is a local domain. Moreover we will factorize θσ as a product of symmetries and transvections for some θ the matrix of which is diagonal.     

Construction of nonnegative symmetric matrices with given spectrum

Let σ = (λ₁, … , λ_n) be the spectrum of a nonnegative symmetric matrix A with the Perron eigenvalue λ₁, a diagonal entry c and let τ = (μ₁, … , μ_m) be the spectrum of a nonnegative symmetric matrix B with the Perron eigenvalue μ₁. We show how to construct a nonnegative symmetric matrix C with the spectrum

(λ₁+max{0,μ₁-c},λ₂,…,λ_n,μ₂,…,μ_m).     

Alexandrov’s inequality and conjectures on some Toeplitz matrices

We study determinant inequalities for certain Toeplitz-like matrices over C. For fixed n and N ? 1, let Q be the n × (n + N − 1) zero-one Toeplitz matrix with Q_ij = 1 for 0 ? j − i ? N − 1 and Q_ij = 0 otherwise. We prove that det(QQ^∗) is the minimum of det(RR^∗) over all complex matrices R with the same dimensions as Q satisfying ∣R_ij∣ ? 1 whenever Q_ij = 1 and R_ij = 0 otherwise. Although R has a Toeplitz-like band structure, it is not required to be actually Toeplitz. Our proof involves Alexandrov’s inequality for polarized determinants and its generalizations. This problem is motivated by Littlewood’s conjecture on the minimum 1-norm of N-term exponential sums on the unit circle. We also discuss polarized Bazin-Reiss-Picquet identities, some connections with k-tree enumeration, and analogous conjectured inequalities for the elementary symmetric functions of QQ^∗.