A numerical comparison of Chebyshev methods for solving fourth order semilinear initial boundary value problems

In solving semilinear initial boundary value problems with prescribed non-periodic boundary conditions using implicit-explicit and implicit time stepping schemes, both the function and derivatives of the function may need to be computed accurately at each time step. To determine the best Chebyshev collocation method to do this, the accuracy of the real space Chebyshev differentiation, spectral space preconditioned Chebyshev tau, real space Chebyshev integration and spectral space Chebyshev integration methods are compared in the L² and W^2,2 norms when solving linear fourth order boundary value problems; and in the L^∞([0,T];L²) and L^∞([0,T];W^2,2) norms when solving initial boundary value problems. We find that the best Chebyshev method to use for high resolution computations of solutions to initial boundary value problems is the spectral space Chebyshev integration method which uses sparse matrix operations and has a computational cost comparable to Fourier spectral discretization.     

On an exact control problem for a semilinear wave equation with an unknown source

The exact controllability of a semilinear wave equation, with Dirichlet boundary control on a part of the boundary and an unknown source, is shown. The nonlinear term has at most a linear growth, the initial and target spaces are L²(Ω)×H⁻¹(Ω).     

The mixed initial-boundary value problem for diagonalizable quasilinear hyperbolic systems in the first quadrant

In this paper, we investigate the mixed initial-boundary value problem for diagonalizable quasilinear hyperbolic systems with nonlinear boundary conditions on a half-unbounded domain . Under the assumptions that system is strictly hyperbolic and linearly degenerate, we obtain the global existence and uniqueness of C¹ solutions with the bounded L¹∩L^∞ norm of the initial data as well as their derivatives and appropriate boundary condition. Based on the existence results of global classical solutions, we also prove that when t tends to infinity, the solutions approach a combination of C¹ travelling wave solutions. Under the appropriate assumptions of initial and boundary data, the results can be applied to the equation of time-like extremal surface in Minkowski space R^1+(1+n).     

Another approach to the heat flow for Plateau problem

For the minimal surfaces in Rⁿ with Plateau boundary condition and establish the global existence and uniqueness of the flow as well as the continuous dependence of the initial datum.     

Global solvability for the Kirchhoff equations in exterior domains of dimension three

We consider the initial (boundary) value problem for the Kirchhoff equations in exterior domains or in the whole space of dimension three, and show that these problems admit time-global solutions, provided the norms of the initial data in the usual Sobolev spaces of appropriate order are sufficiently small. We obtain uniform estimates of the L¹(R) norms with respect to time variable at each point in the domain, of solutions of initial (boundary) value problem for the linear wave equations. We then show that the estimates above yield the unique global solvability for the Kirchhoff equations.     

Heat flow and Hardy inequality in complete Riemannian manifolds with singular initial conditions

Upper bounds are obtained for the heat content of an open set D with singular initial condition f on a complete Riemannian manifold, provided (i) the Dirichlet-Laplace-Beltrami operator satisfies a strong Hardy inequality, and (ii) f satisfies an integrability condition. Precise asymptotic results for the heat content are obtained for an open bounded and connected set D in Euclidean space with C² boundary, and with initial condition f(x)=δ(x)^−α,0<α<2, where δ(x) is the distance from x to the boundary of D.     

A nonlinear singular diffusion equation with source

In this paper, the existence, uniqueness and dependence on initial value of solution for a singular diffusion equation with nonlinear boundary condition are discussed. It is proved that there exists a unique global smooth solution which depends on initial data in L¹ continuously.     

Global existence and decay of energy to systems of wave equations with damping and supercritical sources

This paper is concerned with a system of nonlinear wave equations with supercritical interior and boundary sources and subject to interior and boundary damping terms. It is well-known that the presence of a nonlinear boundary source causes significant difficulties since the linear Neumann problem for the single wave equation is not, in general, well-posed in the finite-energy space H ¹(Ω) × L ²(?Ω) with boundary data from L ²(?Ω) (due to the failure of the uniform Lopatinskii condition). Additional challenges stem from the fact that the sources considered in this article are non-dissipative and are not locally Lipschitz from H ¹(Ω) into L ²(Ω) or L ²(?Ω). With some restrictions on the parameters in the system and with careful analysis involving the Nehari Manifold, we obtain global existence of a unique weak solution and establish (depending on the behavior of the dissipation in the system) exponential and algebraic uniform decay rates of energy. Moreover, we prove a blow-up result for weak solutions with nonnegative initial energy.     

On decay and blow-up of the solution for a viscoelastic wave equation with boundary damping and source terms

In this paper we consider the decay and blow-up properties of a viscoelastic wave equation with boundary damping and source terms. We first extend the decay result (for the case of linear damping) obtained by Lu et al. (On a viscoelastic equation with nonlinear boundary damping and source terms: Global existence and decay of the solution, Nonlinear Analysis: Real World Applications 12 (1) (2011), 295-303) to the nonlinear damping case under weaker assumption on the relaxation function g(t). Then, we give an exponential decay result without the relation between g^′(t) and g(t) for the linear damping case, provided that ‖g‖_L¹(0,∞) is small enough. Finally, we establish two blow-up results: one is for certain solutions with nonpositive initial energy as well as positive initial energy for both the linear and nonlinear damping cases, the other is for certain solutions with arbitrarily positive initial energy for the linear damping case.     

Smooth global Lagrangian flow for the 2D Euler and second-grade fluid equations

We present a very simple proof of the global existence of a C^∞ Lagrangian flow map for the 2D Euler and second-grade fluid equations (on a compact Riemannian manifold with boundary) which has C^∞ dependence on initial data u₀ in the class of H^s divergence-free vector fields for s > 2.     

Convergence rate of the solution toward boundary layer solution for initial-boundary value problem of the 2-D viscous conservation laws

In this paper, we study the large time behavior of the solution to the initial boundary value problem for 2-D viscous conservation laws in the space x ? bt. The global existence and the asymptotic stability of a stationary solution are proved by Kawashima et al. [1]. Here, we investigate the convergence rate of solution toward the boundary layer solution with the non-degenerate case where f′(u₊) − b < 0. Based on the estimate in the H² Sobolev space and via the weighted energy method, we draw the conclusion that the solution converges to the corresponding boundary layer solution with algebraic or exponential rate in time, under the assumption that the initial perturbation decays with algebraic or exponential in the spatial direction.     

Uniform decay rates for the energy of weakly damped defocusing semilinear Schrödinger equations with inhomogeneous Dirichlet boundary control

In this paper, we study the open loop stabilization as well as the existence and regularity of solutions of the weakly damped defocusing semilinear Schrödinger equation with an inhomogeneous Dirichlet boundary control. First of all, we prove the global existence of weak solutions at the H¹-energy level together with the stabilization in the same sense. It is then deduced that the decay rate of the boundary data controls the decay rate of the solutions up to an exponential rate. Secondly, we prove some regularity and stabilization results for the strong solutions in H²-sense. The proof uses the direct multiplier method combined with monotonicity and compactness techniques. The result for weak solutions is strong in the sense that it is independent of the dimension of the domain, the power of the nonlinearity, and the smallness of the initial data. However, the regularity and stabilization of strong solutions are obtained only in low dimensions with small initial and boundary data.     

Analysis and control of a scalar conservation law modeling a highly re-entrant manufacturing system

In this paper, we study a scalar conservation law that models a highly re-entrant manufacturing system as encountered in semi-conductor production. As a generalization of Coron et al. (2010) [14], the velocity function possesses both the local and nonlocal character. We prove the existence and uniqueness of the weak solution to the Cauchy problem with initial and boundary data in L^∞. We also obtain the stability (continuous dependence) of both the solution and the out-flux with respect to the initial and boundary data. Finally, we prove the existence of an optimal control that minimizes, in the Lp-sense with 1?p?∞, the difference between the actual out-flux and a forecast demand over a fixed time period.     

Periodic solutions of a nonlinear propagation equation via a fixed point argument

In this work the initial value problem for the equation $$u_t + \beta u_x + yf(u)_x - \delta u_{xxt} = g,\forall x \in R, \forall t \in [0,T],$$ with periodic boundary conditions is interpreted in the sense of periodic distributions and studied via fixed point arguments. Weak solutions exist iff∈C ⁰ (R) andg∈L ^∞(L ²(0,1)). Moreover, regularity inf, g and the initial data implies regularity of solutions.     

Initial boundary value problems for the compressible viscoelastic fluid

This paper is concerned with the initial boundary value problems of the system modeling the compressible viscoelastic fluid of Oldroyd type. The global in time solution is proved to exist uniquely near the equilibrium state in H².     

On smooth dependence on initial conditions for dissipative PDEs, an ODE-type approach

Using elementary ODE tools we show that for a class of dissipative PDEs, including Navier-Stokes equations on the plane with periodic boundary conditions, Galerkin projections are C¹-convergent, i.e. the partial derivatives with respect to initial conditions for Galerkin projections converge giving rise to smooth dependence on initial conditions for full PDE.     

Smoothness of the trajectories of ideal fluid particles with Yudovich vorticities in a planar bounded domain

We consider the incompressible Euler equations in a (possibly multiply connected) bounded domain Ω of R², for flows with bounded vorticity, for which Yudovich (1963) proved in [29] global existence and uniqueness of the solution. We prove that if the boundary ∂Ω of the domain is C^∞ (respectively Gevrey of order M?1) then the trajectories of the fluid particles are C^∞ (respectively Gevrey of order M+2). Our results also cover the case of “slightly unbounded” vorticities for which Yudovich (1995) extended his analysis in [30]. Moreover if in addition the initial vorticity is Hölder continuous on a part of Ω then this Hölder regularity propagates smoothly along the flow lines. Finally we observe that if the vorticity is constant in a neighborhood of the boundary, the smoothness of the boundary is not necessary for these results to hold.     

Mixed initial boundary value problem for hyperbolic geometric flow on Riemann surfaces

The hyperbolic geometric flow equations is introduced recently by Kong and Liu motivated by Einstein equation and Hamilton Ricci flow. In this paper, we consider the mixed initial boundary value problem for hyperbolic geometric flow, and prove the global existence of classical solutions. The results show that, for any given initial metric on R² in certain class of metric, one can always choose suitable initial velocity symmetric tensor such that the solutions exist, and the scalar curvature corresponding to the solution metric gij keeps bounded. If the initial velocity tensor does not satisfy the certain conditions, the solutions will blow up at a finite time. Some special explicit solutions to the reduced equation are given.     

On the heat flow of equation of surfaces of constant mean curvatures

We consider the initial and boundary value problem of heat flow of equation of surfaces of constant mean curvatures. We give sufficient conditions on the initial data such that the heat flow develops finite time singularity. We also provide a new set of initial data to guarantee the existence of global regular solutions to the heat flow that converges to zero in H ¹ exponentially as time goes to infinity.     

A fast boundary element algorithm for time-dependant potential problems

The numerical solution of time-dependant potential problems via the boundary element has been crippled by the high computational cost due to the inherent time history constraint in the integral representation. Using a boundary-only formulation, the time integrations, at any instant in time, have to be evaluated starting from the initial time. This time-history dependence becomes impractical and inadequate for problems where computations are to be performed for extended times. This also made the boundary element uncompetitive compared to the domain-mesh based methods, such as finite difference and finite element methods, for the solution of transient potential problems. Generally, the evaluation of the potential at N domain points using M boundary points at the Kth time step requires an amount of computer operations of the order O(KM²+KNM). This paper presents an algorithm which requires a computational cost of the order of only O(M²+NM), where the dependence from the past K-steps is removed. The algorithm combines the boundary element method and a scheme, which uses virtual collocation points and radial basis functions to approximate the domain integral.