Numerical solutions with a priori error bounds for coupled self-adjoint time dependent partial differential systems

This paper is concerned with the construction of accurate continuous numerical solutions for partial self-adjoint differential systems of the type (P(t) u_t)_t = Q(t)u_xx, u(0, t) = u(d, t) = 0, u(x, 0) = f(x), u_t(x, 0) = g(x), 0 ≤ x ≤ d, t >- 0, where P(t), Q(t) are positive definite oR^r×r-valued functions such that P′(t) and Q′(t) are simultaneously semidefinite (positive or negative) for all t ≥ 0. First, an exact theoretical series solution of the problem is obtained using a separation of variables technique. After appropriate truncation strategy and the numerical solution of certain matrix differential initial value problems the following question is addressed. Given T > 0 and an admissible error ϵ > 0 how to construct a continuous numerical solution whose error with respect to the exact series solution is smaller than ϵ, uniformly in D(T) = {(x, t); 0 ≤ x ≤ d, 0 ≤ t ≤ T}. Uniqueness of solutions is also studied.     

Hitting time of a moving boundary for a diffusion

We obtain asymptotic estimates for the quantity r = log P[T_f[rang]t] as t → ∞ where T_f = inf\s{s : |X(s)|[rang]f(s)\s} and X is a real diffusion in natural scale with generator a(x) d²(·)/dx² and the ‘boundary’ f(s) is an increasing function. We impose regular variation on a and f and the result is expressed as r = ∫^t₀ λ₁ (f(s) ds(1 + o(1)) where λ₁(f) is the smallest eigenvalue for the process killed at ±f.     

On the behavior of the solution of the Burgers equation as the viscosity goes to zero

The solution of the initial boundary-value problem u_?′ ? ?D²u_? + u_?Du_? = f on (a, b) x(0, T), u_?(a, t) = u_?(b, t) = 0 and u_?(x, 0) = 0 on (a, b), is shown to converge to the solution of the limiting equation as the viscosity tends to zero. Estimates on the rate of convergence are given.     

On a nonconvolution Volterra resolvent

Under fairly weak assumptions, the solutions of the system of Volterra equations x(t) = ∝₀^ta(t, s) x(s) ds + f(t), t > 0, can be written in the form x(t) = f(t) + ∝₀^tr(t, s) f(s) ds, t > 0, where r is the resolvent of a, i.e., the solution of the equation r(t, s) = a(t, s) + ∝₀^ta(t, v) r(v, s)dv, 0 < s < t. Conditions on a are given which imply that the resolvent operator f ∝₀^tr(t, s) f(s) ds maps a weighted L¹ space continuously into another weighted L¹ space, and a weighted L^∞ space into another weighted L^∞ space. Our main theorem is used to study the asymptotic behavior of two differential delay equations.     

A direct numerical procedure for the Cauchy problem for the heat equation

For the Cauchy problem, u_t = u_xx, 0 < x < 1, 0 < t ? T, u(0, t) = f(t), 0 < t ? T, u_x(0, t) = g(t), 0 < t ? T, a direct numerical procedure involving the elementary solution of υ_t = υ_xx, 0 < x, 0 < t ? T, υ_x(0, t) = g(t), 0 < t ? T, υ(x, 0) = 0, 0 < x and a Taylor's series computed from f(t) ? υ(0, t) is studied. Continuous dependence better than any power of logarithmic is obtained. Some numerical results are presented.     

Free and forced vibrations of nonlinear wave equations in ball

First, we shall deal with the free vibrations of a nonlinear radially symmetric wave equation (∂_t²−△)u=f(r,u) in n-dimensional ball B_a with center at the origin and radius a, where f is smooth, monotone decreasing in u, and satisfies f(r,0)=0. f(r,u) has asymptotic properties . For n=1,3 we shall show the existence of infinitely many radially symmetric time-periodic solutions with different periods of irrational multiple of a. Second, we shall deal with BVP for a forced nonlinear wave equation (∂_t²−△)u=εg(r,t,u), where g is T-periodic in t and ε is a small parameter. Under some Diophantine condition on a/T we shall show the existence of time-periodic solutions of the BVP. For 1?n?5 we shall construct infinitely many T satisfying the above Diophantine inequality, using asymptotic expansions of the zero points of the Bessel functions.     

Weak solution of (r, (x, u(x)) u′(x))′ = (Fu)′(x) with r(0, u(0))u′(0) = ku(0), r(L, u(L))u′ (L) = hu(L) and k, h are suitable elements of [0, ∞]

We consider weak solutions to the nonlinear boundary value problem (r, (x, u(x)) u′(x))′ = (Fu)′(x) with r(0, u(0)) u′(0) = ku(0), r(L, u(L)) u′(L) = hu(L) and k, h are suitable elements of [0, ∞]. In addition to studying some new boundary conditions, we also relax the constraints on r(x, u) and (Fu)(x). r(x, u) > 0 may have a countable set of jump discontinuities in u and r(x, u)^?1?L_q((0, L) × (0, p)). F is an operator from a suitable set of functions to a subset of L_p(0, L) which have nonnegative values. F includes, among others, examples of the form (Fu)(x) = (1 ? H(x ? x₀)) u(x₀), (Fu)(x) = ∫_x^Lf(y, u(y)) dy where f(y, u) may have a countable set of jump discontinuities in u or F may be chosen so that (Fu)′(x) = ? g(x, u(x)) u′(x) ? q(x) u(x) ? f(x, u(x)) where q is a distributional derivative of an L₂(0, L) function.     

Potential theory associated with Uhlenbeck-Ornstein process

Some parallel results of Gross' paper (Potential theory on Hilbert space, J. Functional Analysis1 (1967), 123–181) are obtained for Uhlenbeck-Ornstein process U(t) in an abstract Wiener space (H, B, i). Generalized number operator $\text{N}$ is defined by $\text{N}$f(x) = ?lim_∈←0{E[f($\text{U}$(τ^∈_ξ))] ? f(x)}/E[τ^∈_ξ, where τ_x^? is the first exit time of U(t) starting at x from the ball of radius ? with center x. It is shown that $\text{N}$f(x) = ?trace D²f(x)+〈Df(x),x〉 for a large class of functions f. Let r_t(x, dy) be the transition probabilities of U(t). The λ-potential G_λf, λ > 0, and normalized potential Rf of f are defined by G_λf(X) = ∫₀^∞e^?λtr_tf(x) dt and Rf(x) = ∫₀^∞ [r_tf(x) ? r_tf(0)] dt. It is shown that if f is a bounded Lip-1 function then trace D²G_λf(x) ? 〈DG_λf(x), x〉 = ?f(x) + λG_λf(x) and trace D²Rf(x) ? 〈DRf(x), x〉 = ?f(x) + ∫_Bf(y)p₁(dy), where p₁ is the Wiener measure in B with parameter 1. Some approximation theorems are also proved.     

Single point blow-up for a semilinear initial value problem

The initial value problem on [?R, R] is considered: u_t(t, x) = u_xx(t, x) + u(t, x)^γu(t, ±R) = 0u(0, x) = ?(x), where ? ? 0 and γ is a fixed large number. It is known that for some initial values ? the solution u(t, x) exists only up to some finite time T, and that ∥u(t, ·)∥_∞ → ∞ as t → T. For the specific initial value ? = kψ, where ψ ? 0, ψ_xx + ψ^γ = 0, ψ(±R) = 0, k is sufficiently large, it is shown that if x ≠ 0, then lim_{t → T}u(t, x) and lim_{t → T}u_x(t, x) exist and are finite. In other words, blow-up occurs only at the point x = 0.     

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.     

Solvability of Mixed Problems for the Klein–Gordon–Fock Equation in the Class <Emphasis Type="Italic">L</Emphasis><Subscript><Emphasis Type="Italic">p</Emphasis></Subscript> for <Emphasis Type="Italic">p</Emphasis> ≥ 1

We prove that the mixed problem for the Klein–Gordon–Fock equation u _tt (x, t) ? u _xx (x, t) + au(x, t) = 0, where a ≥ 0, in the rectangle Q _T = [0 ≤ x ≤ l] × [0 ≤ t ≤ T] with zero initial conditions and with the boundary conditions u(0, t) = μ(t) ∈ L _p [0, T ], u(l, t) = 0, has a unique generalized solution u(x, t) in the class L _p (Q _T ) for p ≥ 1. We construct the solution in explicit analytic form.     

Asymptotic behavior for differential equations which cannot be locally linearized

For functions ? which are continuous and locally Lipschitz the authors define a multi-valued differential D_f and prove that if all solutions of the multi-valued differential equation u′ ? D_f(u) approach zero as t → ∞, then all solutions x(·) of x′ = ?(x) with small |x(0)| approach zero exponentially as t → ∞. If ? is continuously differentiable, then D_f coincides with the (single-valued) Frechet differential of ?. Other results on the asymptotic behavior of solutions of perturbed, multi-valued differential equations are presented.     

Existence of a solution to a nonlinear equation

Equation (−Δ+k²)u+f(u)=0 in D, u_|∂D=0, where k=const>0 and D⊂R³ is a bounded domain, has a solution if is a continuous function in the region |u|?a, piecewise-continuous in the region |u|?a, with finitely many discontinuity points uj such that f(uj±0) exist, and uf(y)?0 for |u|?a, where a?0 is an arbitrary fixed number.     

Local existence of solutions to some degenerate parabolic equation associated with the p-Laplacian

The existence of local (in time) solutions of the initial-boundary value problem for the following degenerate parabolic equation: ut(x,t)−Δpu(x,t)−|u|^q−2u(x,t)=f(x,t), (x,t)∈Ω×(0,T), where 2?p<q<+∞, Ω is a bounded domain in RN, is given and Δp denotes the so-called p-Laplacian defined by Δpu:=∇⋅(|∇u|^p−2∇u), with initial data u₀∈Lr(Ω) is proved under r>N(q−p)/p without imposing any smallness on u₀ and f. To this end, the above problem is reduced into the Cauchy problem for an evolution equation governed by the difference of two subdifferential operators in a reflexive Banach space, and the theory of subdifferential operators and potential well method are employed to establish energy estimates. Particularly, Lr-estimates of solutions play a crucial role to construct a time-local solution and reveal the dependence of the time interval [0,T₀] in which the problem admits a solution. More precisely, T₀ depends only on Lr|u₀| and f.     

Analytic-numerical solutions with a priori error bounds for a class of strongly coupled mixed partial differential systems

This paper deals with the construction of analytic-numerical solutions with a priori error bounds for systems of the type u_t = Au_xx, u(0,t) + u_x(0,t) = 0, Bu(1,t) + Cu_x(1,t) = 0, 0 < x < 1, t > 0, u(x,0) = f(x). Here A, B, C are matrices for which no diagonalizable hypothesis is assumed. First an exact series solution is obtained after solving appropriate vector Sturm-Liouville-type problems. Given an admissible error ε and a bounded subdomain D, after appropriate truncation an approximate solution constructed in terms of data and approximate eigenvalues is given so that the error is less than the prefixed accuracy ε, uniformly in D.     

Scattering length and perturbations of −Δ by positive potentials

This paper extends a result of Fujita [On the blowing up of solutions to the Cauchy problem for u_t = Δu + u^{1 + a}, J. Faculty Science, U. of Tokyo 13 (1966), 109–124] to show that solutions u = u(t, x) for t > 0 and x?R² to the equation u_t = Δu + u² with u(0, x) = a(x) must grow at a rate faster than exp(∥x∥²) at some finite time t, as long as a(x) is nonnegative and not almost everywhere zero.     

Inverse Problems for Time-Dependent Singular Heat Conductivities: Multi-Dimensional Case

We consider an inverse boundary value problem for the heat equation ? _t u = div (γ? _x u) in (0, T) × Ω, u = f on (0, T) × ?Ω, u| _t=0 = u ₀, in a bounded domain Ω ? ? ⁿ , n ≥ 2, where the heat conductivity γ(t, x) is piecewise constant and the surface of discontinuity depends on time: γ(t, x) = k ² (x ∈ D(t)), γ(t, x) = 1 (x ∈ Ω?D(t)). Fix a direction e* ∈ 𝕊 ^n?1 arbitrarily. Assuming that ?D(t) is strictly convex for 0 ≤ t ≤ T, we show that k and sup {e*·x; x ∈ D(t)} (0 ≤ t ≤ T), in particular D(t) itself, are determined from the Dirichlet-to-Neumann map : f → ?_ν u(t, x)|_{(0, T)×?Ω}. The knowledge of the initial data u ₀ is not used in the proof. If we know min_0≤t≤T (sup _x∈D(t) x·e*), we have the same conclusion from the local Dirichlet-to-Neumann map. Numerical examples of stationary and moving circles inside the unit disk are shown. The results have applications to nondestructive testing. Consider a physical body consisting of homogeneous material with constant heat conductivity except for a moving inclusion with different conductivity. Then the location and shape of the inclusion can be monitored from temperature and heat flux measurements performed at the boundary of the body. Such a situation appears for example in blast furnaces used in ironmaking.     

Periodic solutions for wave equations with variable coefficients with nonlinear localized damping

We discuss the existence of periodic solutions to the wave equation with variable coefficients utt−div(A(x)∇u)+ρ(x,ut)=f(x,t) with Dirichlet boundary condition. Here ρ(x,v) is a function like ρ(x,v)=a(x)g(v) with g^′(v)?0 where a(x) is nonnegative, being positive only in a neighborhood of a part of the domain.     

Cauchy problem for semilinear wave equation with time-dependent metrics

We establish the existence of a weak solution u of the semilinear wave equation where a(t,x) is equal to 1 outside a compact set with respect to x and a non-linear term f_k which satisfies |f_k(u)|≤C|u|^k. For some non-trapping time-periodic perturbations a(t,x), we obtain the long time existence of a solution from little initial data.