Characteristics with singularities and the boundary values of the asymptotic solution of the Cauchy problem for a degenerate wave equation

We consider the Cauchy problem with spatially localized initial data for the twodimensional wave equation degenerating on the boundary of the domain. This problem arises, in particular, in the theory of tsunami wave run-up on a shallow beach. Earlier, S. Yu. Dobrokhotov, V. E. Nazaikinskii, and B. Tirozzi developed a method for constructing asymptotic solutions of this problem. The method is based on a modified Maslov canonical operator and on characteristics (trajectories) unbounded in the momentum variables; such characteristics are nonstandard from the viewpoint of the theory of partial differential equations. In a neighborhood of the velocity degeneration line, which is a caustic of a special form, the canonical operator is defined via the Hankel transform, which arises when applying Fock’s quantization procedure to the canonical transformation regularizing the above-mentioned nonstandard characteristics in a neighborhood of the velocity degeneration line (the boundary of the domain). It is shown in the present paper that the restriction of the asymptotic solutions to the boundary is determined by the standard canonical operator, which simplifies the asymptotic formulas for the solution on the boundary dramatically; for initial perturbations of special form, the solutions can be expressed via simple algebraic functions.     

Nonstrictly hyperbolic systems of conservation laws of the conjugate type

We consider an inverse boundary problem for a general second order self-adjoint elliptic differential operator on a compact differential manifold with boundary. The inverse problem is that of the reconstruction of the manifold and operator via all but finite number of eigenvalues and traces on the boundary of the corresponding eigenfunctions of the operator. We prove that the data determine the manifold and the operator to within the group of the generalized gauge transformations. The proof is based upon a procedure of the reconstruction of a canonical object in the orbit of the group. This object, the canonical Schrödinger operator, is uniquely determined via its incomplete boundary spectral data.     

Traces, extensions and co-normal derivatives for elliptic systems on Lipschitz domains

For functions from the Sobolev space Hs(Ω), , definitions of non-unique generalized and unique canonical co-normal derivative are considered, which are related to possible extensions of a partial differential operator and its right-hand side from the domain Ω, where they are prescribed, to the domain boundary, where they are not. Revision of the boundary value problem settings, which makes them insensitive to the generalized co-normal derivative inherent non-uniqueness are given. It is shown, that the canonical co-normal derivatives, although defined on a more narrow function class than the generalized ones, are continuous extensions of the classical co-normal derivatives. Some new results about trace operator estimates and Sobolev spaces characterizations, are also presented.     

Optimal fishery with coastal catch

In many spatial resource models, it is assumed that an agent is able to harvest the resource over the complete spatial domain. However, agents frequently only have access to a resource at particular locations at which a moving biomass, such as fish or game, may be caught or hunted. Here, we analyze an infinite time‐horizon optimal control problem with boundary harvesting and (systems of) parabolic partial differential equations as state dynamics. We formally derive the associated canonical system, consisting of a forward–backward diffusion system with boundary controls, and numerically compute the canonical steady states and the optimal time‐dependent paths, and their dependence on parameters. We start with some one‐species fishing models, and then extend the analysis to a predator–prey model of the Lotka–Volterra type. The models are rather generic, and our methods are quite general, and thus should be applicable to large classes of structurally similar bioeconomic problems with boundary controls. Recommedations for Resource Managers

• Just like ordinary differential equation‐constrained (optimal) control problems and distributed partial differential equation (PDE) constrained control problems, boundary control problems with PDE state dynamics may be formally treated by the Pontryagin's maximum principle or canonical system formalism (state and adjoint PDEs).
• These problems may have multiple (locally) optimal solutions; a first overview of suitable choices can be obtained by identifying canonical steady states.
• The computation of canonical paths toward some optimal steady state yields temporal information about the optimal harvesting, possibly including waiting time behavior for the stock to recover from a low‐stock initial state, and nonmonotonic (in time) harvesting efforts.
• Multispecies fishery models may lead to asymmetric effects; for instance, it may be optimal to capture a predator species to protect the prey, even for high costs and low market values of the predators.

     

Uniform asymptotics of the boundary values of the solution in a linear problem on the run-up of waves on a shallow beach

We consider the Cauchy problem with spatially localized initial data for a two-dimensional wave equation with variable velocity in a domain Ω. The velocity is assumed to degenerate on the boundary ?Ω of the domain as the square root of the distance to ?Ω. In particular, this problems describes the run-up of tsunami waves on a shallow beach in the linear approximation. Further, the problem contains a natural small parameter (the typical source-to-basin size ratio) and hence admits analysis by asymptotic methods. It was shown in the paper “Characteristics with singularities and the boundary values of the asymptotic solution of the Cauchy problem for a degenerate wave equation” [1] that the boundary values of the asymptotic solution of this problem given by a modified Maslov canonical operator on the Lagrangian manifold formed by the nonstandard characteristics associatedwith the problemcan be expressed via the canonical operator on a Lagrangian submanifold of the cotangent bundle of the boundary. However, the problem as to how this restriction is related to the boundary values of the exact solution of the problem remained open. In the present paper, we show that if the initial perturbation is specified by a function rapidly decaying at infinity, then the restriction of such an asymptotic solution to the boundary gives the asymptotics of the boundary values of the exact solution in the uniform norm. To this end, we in particular prove a trace theorem for nonstandard Sobolev type spaces with degeneration at the boundary.     

On a variational problem with unknown boundaries and the determination of optimal shapes of elastic bodies : PMM vol. 39, n 6, 1975, pp. 1082–1092

The optimization problem is considered for a partial differential equation of elliptic type. The boundary of the domain in which the equation is given emerges as the control function and is to be determined from the condition of the extremum of the integral of the solution of the boundary value problem. Seeking the extremals is reduced to solving a va national problem without differential constraints. Necessary conditions for optimality are obtained, and shapes of elastic bars possessing the maximum stiffness under torsion are found with their aid.     

The quasilinearization method on an unbounded domain

We apply a method of quasilinearization to a boundary value problem for an ordinary differential equation on an unbounded domain. A uniquely determined Green's function, which is integrable and of fixed sign, is employed. The hypotheses to apply the quasilinearization method imply uniqueness of solutions. The quasilinearization method generates a bilateral iteration scheme in which the iterates converge monotonically and quadratically to the unique solution.

     

Universal boundary value problem for equations of mathematical physics

A new statement of a boundary value problem for partial differential equations is discussed. An arbitrary solution to a linear elliptic, hyperbolic, or parabolic second-order differential equation is considered in a given domain of Euclidean space without any constraints imposed on the boundary values of the solution or its derivatives. The following question is studied: What conditions should hold for the boundary values of a function and its normal derivative if this function is a solution to the linear differential equation under consideration? A linear integral equation is defined for the boundary values of a solution and its normal derivative; this equation is called a universal boundary value equation. A universal boundary value problem is a linear differential equation together with a universal boundary value equation. In this paper, the universal boundary value problem is studied for equations of mathematical physics such as the Laplace equation, wave equation, and heat equation. Applications of the analysis of the universal boundary value problem to problems of cosmology and quantum mechanics are pointed out.     

Analysis of united boundary‐domain integro‐differential and integral equations for a mixed BVP with variable coefficient

The mixed (Dirichlet–Neumann) boundary‐value problem for the ‘Laplace’ linear differential equation with variable coefficient is reduced to boundary‐domain integro‐differential or integral equations (BDIDEs or BDIEs) based on a specially constructed parametrix. The BDIDEs/BDIEs contain integral operators defined on the domain under consideration as well as potential‐type operators defined on open sub‐manifolds of the boundary and acting on the trace and/or co‐normal derivative of the unknown solution or on an auxiliary function. Some of the considered BDIDEs are to be supplemented by the original boundary conditions, thus constituting boundary‐domain integro‐differential problems (BDIDPs). Solvability, solution uniqueness, and equivalence of the BDIEs/BDIDEs/BDIDPs to the original BVP, as well as invertibility of the associated operators are investigated in appropriate Sobolev spaces. Copyright © 2005 John Wiley & Sons, Ltd.     

Asymptotic behavior for the filtration equation in domains with noncompact boundary

We consider the initial value boundary problem with zero Neumann data for an equation modeled after the porous media equation, with variable coefficients. The spatial domain is unbounded and shaped like a (general) paraboloid, and the solution u is integrable in space and nonnegative. We show that the asymptotic profile for large times of u is one dimensional and given by an explicit function, which can be regarded as the fundamental solution of a one-dimensional differential equation with weights. In the case when the domain is a cone or the whole space (Cauchy problem), we obtain a genuine multidimensional profile given by the well-known Barenblatt solution.     

Boundary value problem for a second-order elliptic equation in the exterior of an ellipse

We consider a boundary value problem for a second-order linear elliptic differential equation with constant coefficients in a domain that is the exterior of an ellipse. The boundary conditions of the problem contain the values of the function itself and its normal derivative. We give a constructive solution of the problem and find the number of solvability conditions for the inhomogeneous problem as well as the number of linearly independent solutions of the homogeneous problem. We prove the boundary uniqueness theorem for the solutions of this equation.     

A multiplicative method of solving boundary value problems of shell theory

The boundary conditions are transferred to an arbitrarily chosen point by multiplication of matrices (multiplicatively). The transfer matrices of the boundary conditions are an analytic solution of a system of first-order linear ordinary differential equations in canonical form of the mechanics of the deformation of shells in the form of values of Cauchy–Krylov functions. At an arbitrarily chosen point, the boundary conditions are combined in a system of algebraic equations in matrix form, columns of the unknown quantities of which are parameters of the required values of the problem. The effectiveness of the method – the simplicity with which it can be realized on a computer, the low costs of computer time and the RAM – is based on the multiplicative transfer of the boundary conditions into matrix form. The class of problems is limited by the possibilities of the Fourier method of separation of the variables in partial differential equations.     

The method of lines for Poisson's equation with nonlinear or free boundary conditions

Summary The method of lines is used to solve Poisson's equation on an irregular domain with nonlinear or free boundary conditions. The partial differential equation is approximated by a system of second order ordinary differential equations subject to multi-point boundary conditions. The system is solved with an SOR iteration which employs invariant imbedding for each one dimensional problem. An application of the method to a boundary control problem and to a free surface problem arising in electrochemical machining is described. Finally, some theoretical convergence results are presented for a model problem with radiative boundary conditions on fixed boundaries.This work was supported by the U.S. Army Research Office under Grant DA-AG29-76-G-0261     

Simple Asymptotics for a Generalized Wave Equation with Degenerating Velocity and Their Applications in the Linear Long Wave Run-Up Problem

Asymptotic solutions of the wave equation degenerating on the boundary of the domain (where the wave propagation velocity vanishes as the square root of the distance from the boundary) can be represented with the use of a modified canonical operator on a Lagrangian submanifold, invariant with respect to theHamiltonian vector field, of the nonstandard phase space constructed by the authors in earlier papers. The present paper provides simple expressions in a neighborhood of the boundary for functions represented by such a canonical operator and, in particular, for the solution of the Cauchy problem for the degenerate wave equation with initial data localized in a neighborhood of an interior point of the domain.     

A numerical study of boundary layer flow of viscous incompressible fluid past an inclined stretching sheet and heat transfer using nonpolynomial spline method

In the present paper, we study the boundary layer flow of viscous incompressible fluid over an inclined stretching sheet with body force and heat transfer. Considering the stream function, we convert the boundary layer equation into nonlinear third-order ordinary differential equation together with appropriate boundary conditions in an infinite domain. The nonlinear boundary value problem has been linearized by using the quasilinearization technique. Then, we develop a nonpolynomial spline method, which is used to solve the flow problem. The convergence analysis of the method is also discussed. We study the velocity function for different angles of inclination and Froude number with the help of various graphs and tables. Then using these in heat convection flow, we obtain the expression for temperature field. Skin friction is also calculated. The various results have been given in tables. At last, we calculated the Nusselt number.     

New closed-form Green function and integral formula for a thermoelastic quadrant

In this study a new Green’s function and a new Green-type integral formula for a boundary value problem (BVP) in thermoelastostatics for a quadrant are derived in closed form. On the boundary semi-straight-lines twice mixed homogeneous mechanical boundary conditions (one boundary semi-straight-line is free of loadings and normal displacements and tangential stresses are prescribed on the other one) are prescribed. The thermoelastic displacements are subject by a heat source applied in the inner points of the quadrant and by mixed non-homogeneous boundary heat conditions (on one boundary semi-straight-line the temperature is prescribed and the heat flux is given on the other one). When thermoelastic Green’s function is derived the thermoelastic displacements are generated by an inner unit point heat source, described by δ Dirac’s function. All results are obtained in elementary functions that are formulated in a special theorem. A closed-form solution for a particular BVP of thermoelastostatics for a quadrant also is included. Using the proposed approach it is possible to extend the obtained for quadrant results to any other canonical Cartesian domain.     

On variations in the solutions of integro-differential equations of mixed problems of elasticity theory for domain variations

The behaviour of the solution of the boundary value problem for a pseudodifferential equation (PDE), Green's function of this problem, and also some of their local and global characteristics, during variation of the domain is investigated. Formulas are proposed that enable the solution of a broad class of PDE in a domain to be expressed in terms of the solution in the near domain. Local characteristics of the solution are expressed in terms of the local characteristics of the solution in the near domain. A double asymptotic form of Green's function for both arguments tending to the domain boundary occurs in the variation formula. The variation of this double asymptotic form as the domain varies is expressed in terms of this same asymptotic form. The system of variation formulas obtained is closed. It enables the PDE solution in the domain to be reduced to the solution of an ordinary differential equation in functional space. The local characteristics of the solution can also be found by this method without calculating the solution itself. If there is sufficient symmetry in the initial operator, then conservation laws in the Noether sense are obtained for its Green's function and its asymptotic form. The behaviour of the quantities under investigation is studied under inversion.

The investigation of variations of the solutions of problems for the variation of the domain occurs in the paper by Hadamard /1/, who studied the variation in conformal mapping and obtained a formula similar to (1.4). The formula for the variation of the solution of the boundary value problem for an elliptic differential equation is obtained in /2/. Variation formulas for the case of the operator of the problem about a crack and a circular domain are obtained in /3, 4/. The Irwin formula /5/ is obtained from formulas (1.4) and (1.21) by substitution.     

Theoretical study of a Bénard-Marangoni problem

In this paper we prove the existence of strong solutions for the stationary Bénard-Marangoni problem in a finite domain flat on the top, bifurcating from the basic heat conductive state. The Bénard-Marangoni problem is a physical phenomenon of thermal convection in which the effects of buoyancy and surface tension are taken into account. This problem is modelled with a system of partial differential equations of the type Navier-Stokes and heat equation. The boundary conditions include crossed boundary conditions involving tangential derivatives of the temperature and normal derivatives of the velocity field. To define tangential derivatives at the boundary, intended in the trace sense, it is necessary order two derivatives in the interior of the domain and thus the boundary term contains as high derivatives as the interior term. We overcome this difficulty by considering the weak formulation, and transforming the boundary integral into an equivalent integral defined in the whole domain. This allows us to reformulate the weak problem with a temperature having only order one weak derivatives. Concerning regularity results, we obtain strong solutions for the stationary Bénard-Marangoni problem.     

Error estimations for some meshless boundary interpolation methods

Meshless methods have become quite popular in numerical treatment of partial differential equations because of their simplicity and the fact that they require neither domain nor boundary mesh. In general, however, they convert the original problem to a highly ill‐conditioned linear system of algebraic equations with a dense matrix. Recently, a special technique has been proposed which circumvents this computational difficulty. This method, called Direct Multi‐Elliptic Interpolation Method, is based on a scattered data interpolation which defines the interpolation function as a solution of a higher order multi‐elliptic equation. Here the boundary version of this meshless method which is based on a multi‐elliptic boundary interpolation is considered. Error estimations are derived justifying the interpolation function to be a good approximation of the solution of the original boundary value problem as well. At the same time, the problem of large, dense and ill‐conditioned matrices as well as the mesh generation are completely avoided. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Compact gradient tracking in shape optimization

In the present paper we consider the minimization of gradient tracking functionals defined on a compact and fixed subdomain of the domain of interest. The underlying state is assumed to satisfy a Poisson equation with Dirichlet boundary conditions. We proof that, in contrast to the situation of gradient tracking on the whole domain, the shape Hessian is not strictly H ^1/2-coercive at the optimal domain which implies ill-posedness of the shape problem under consideration. Shape functional and gradient require only knowledge of the Cauchy data of the state and its adjoint on the boundaries of the domain and the subdomain. These data can be computed by means of boundary integral equations when reformulating the underlying differential equations as transmission problems. Thanks to fast boundary element techniques, we derive an efficient algorithm to solve the problem under consideration.