Analytic-numerical description of asymptotic solutions of a Cauchy problem in a neighborhood of singularities for a linearized system of shallow-water equations

S. Yu. Dobrokhotov, B. Tirozzi, S. Ya. Sekerzh-Zenkovich, A. I. Shafarevich, and their co-authors suggested new effective asymptotic formulas for solving a Cauchy problem with localized initial data for multidimensional linear hyperbolic equations with variable coefficients and, in particular, for a linearized system of shallow-water equations over an uneven bottom in their cycle of papers. The solutions are localized in a neighborhood of fronts on which focal points and self-intersection points (singular points) occur in the course of time, due to the variability of the coefficients. In the present paper, a numerical realization of asymptotic formulas in a neighborhood of singular points of fronts is presented in the case of the system of shallow-water equations, gluing problems for these formulas together with formulas for regular domains are discussed, and also a comparison of asymptotic solutions with solutions obtained by immediate numerical computations is carried out.     

Asymptotic behaviors of the stress fields in the vicinity of dislocations and dislocation segments

We obtain the singular asymptotic behavior of the stress field in the vicinity of a non-planar dislocation in three dimensions and the nearly singular behavior of the full self-force of the dislocation including both glide and climb forces, using asymptotic analysis. We also derive asymptotic formulas for the stress field in the vicinity of a curved dislocation segment. Numerical examples are presented to examine the asymptotic formulas. The obtained formulas can be used for qualitative understanding of the stress tensor associated with dislocations and efficient and accurate calculation of the stress tensor in dislocation dynamics simulations.     

Low-frequency scattering of acoustic waves by a bounded rough surface in a half-plane

The problem of the scattering of harmonic plane waves by a rough half-plane is studied here. The surface roughness is finite. The slope of the irregularity is taken as arbitrary. Two boundary conditions are considered, those of Dirichlet and Neumann. An asymptotic solution is obtained, when the wavelength lambda of the incident wave is much larger than the characteristic length of the roughness iota, by means of the method of matched asymptotic expansions in terms of the small parameter epsilon= 2piiota/lambda. For the Dirichlet problem, the solution of the near and far fields is obtained up to O(epsilon2). The far field solution is given in terms of a coefficient that have a simple explicit expression, which also appears in the corresponding solution to the Neumann problem, already solved. Also the scattering cross section is given by simple formulas to O(epsilon3). It is noted that, for the Dirichlet problem, the leading term is of order epsilon3 which, by contrast, is different from that of the circular cylinder in full space, that is, of order epsilon(-1) (log epsilon)(-2). Some examples display the simplicity of the general results based on conformal mapping, which involve arcs of circle, polygonal lines, surface cracks and the like.     

Cookbook asymptotics for spiral and scroll waves in excitable media

Algebraic formulas predicting the frequencies and shapes of waves in a reaction-diffusion model of excitable media are presented in the form of four recipes. The formulas themselves are based on a detailed asymptotic analysis (published elsewhere) of the model equations at leading order and first order in the asymptotic parameter. The importance of the first order contribution is stressed throughout, beginning with a discussion of the Fife limit, Fife scaling, and Fife regime. Recipes are given for spiral waves and detailed comparisons are presented between the asymptotic predictions and the solutions of the full reaction-diffusion equations. Recipes for twisted scroll waves with straight filaments are given and again comparisons are shown. The connection between the asymptotic results and filament dynamics is discussed, and one of the previously unknown coefficients in the theory of filament dynamics is evaluated in terms of its asymptotic expansion. (c) 2002 American Institute of Physics.     

Propagation of finite amplitude sound waves radiated from a pulsating sphere

The propagation of weakly non-linear acoustic waves is studied, which are radiated from a harmonically pulsating sphere in an inviscid perfect gas. A representation of the solution is presented for a far field equation of the first order, which is closely related to the solution obtained by the method of renormalization. The applicability of the method to the present problem is proved within the first order approximation. The asymptotic solution in the far field is obtained up to the second order by using the method of renormalization. A uniformly valid solution is, to the second approximation, constructed by matching the near field solution with the far field solution. Also investigated is the effect of dimensionless parameters and the second order correction on the acoustic shock formation distances and the non-linear distortion of waveforms.     

An improved variational method

In order to improve the unitarity of the S-matrix, an improved variational formulism is derived by proposing new generating functionals and adopting proper asymptotic boundary conditions for trial relative wave functions. The formulas with the weighted line-column balance for the single-channel and multi-channel scatterings, where the non-central interaction is implicitly considered, are presented. A numerical check is performed with a soluble model in a four coupled channel scattering problem. The result shows that the high accuracy and the unitarity of the S-matrix are reached.     

Analytical formulas for complex permittivity of periodic composites. Estimation of gain and loss enhancement in active and passive composites

ABSTRACT

The asymptotic homogenization method is applied to complex dielectric periodic composites. An equivalence to coupled dielectric problems with real coefficients is shown. This is similar to a piezoelectric problem: an out-plane mechanical displacement and an in-plane electric potential establishing a correspondence principle. Closed-form formulas for the complex dielectric effective tensor in the case of a square array of circular inclusions embedded in a matrix are given. These formulas are written in terms of a real and symmetric matrix which facilitates the implementation of the computational scheme. We also get similar formulas for multilayered complex dielectric composites. The real closed-form formulas are advantageous for estimating gain and loss enhancement properties of active and passive composites in certain volume fraction intervals. Numerical computations are performed and the results are compared with other approaches showing the usefulness of the obtained formulas. This may be of interest in the context of metamaterials.     

Transitive ensembles of random matrices related to orthogonal polynomials

Transitive correlations of eigenvalues for random matrix ensembles intermediate between real symmetric and hermitian, self-dual quaternion and hermitian, and antisymmetric and hermitian are studied. Expressions for exact n-point correlation functions are obtained for random matrix ensembles related to general orthogonal polynomials. The asymptotic formulas in the limit of large matrix dimension are evaluated at the spectrum edges for the ensembles related to the Legendre polynomials. The results interpolate known asymptotic formulas for random matrix eigenvalues.     

An improved variational method

In order to improve the unitarity of the S-matrix, an improved variational formulism is derived by proposing new generating functionals and adopting proper asymptotic boundary conditions for trial relative wave functions. The formulas with the weighted line-column balance for the single-channel and multi-channel scatterings, where the non-central interaction is implicitly considered, are presented. A numerical check is performed with a soluble model in a four coupled channel scattering problem. The result shows that the high accuracy and the unitarity of the S-matrix are reached.     

Numerical and analytical study of wave processes in periodic stratified media

We study the behaviour of the solutions of the Cauchy problem with discontinuous initial data for nonstandard linear partial differential equations modeling wave processes in periodic stratified media. Asymptotic formulas at large t are derived. The found asymptotic formulas are in a good agreement with the results of numerical experiments done by using the analytical computation system REDUCE 3.8.     

A splitting extrapolation for solving nonlinear elliptic equations with d-quadratic finite elements

Nonlinear elliptic partial differential equations are important to many large scale engineering and science problems. For this kind of equations, this article discusses a splitting extrapolation which possesses a high order of accuracy, a high degree of parallelism, less computational complexity and more flexibility than Richardson extrapolation. According to the problems, some domain decompositions are constructed and some independent mesh parameters are designed. Multi-parameter asymptotic expansions are proved for the errors of approximations. Based on the expansions, splitting extrapolation formulas are developed to compute approximations with high order of accuracy on a globally fine grid. Because these formulas only require us to solve a set of smaller discrete subproblems on different coarser grids in parallel instead of on the globally fine grid, a large scale multidimensional problem is turned into a set of smaller discrete subproblems. Additionally, this method is efficient for solving interface problems.     

Exponentially small splittings in Hamiltonian systems

Families of area preserving analytical maps, depending on a small parameter epsilon, are considered, with the case epsilon=0 corresponding to an integrable map. The asymptotic formulas for the splittings of separatrices are derived by the method of analytical continuation of the separatrices to the complex domain. The main terms of the asymptotics are exponentially small with respect to the size of the perturbation. As epsilon tends to zero, the intersection angle of the separatrices can oscillate. The exponent and the oscillatory multiplier of the asymptotic formulas are determined by the position of poles of the homoclinic (heteroclinic) orbit of the limiting flow. Pre-exponential coefficients in the asymptotic formulas contain a multiplier obtained by the numerical study of separatrices of "model" maps in the complex domain.     

Fast asymptotic solutions for sound fields above and below a rigid porous ground

The current study simultaneously addresses the problem of reflection and refraction of sound from a rigid porous ground surface. A more rigorous approach is used to derive more accurate asymptotic solutions that can be cast in a convenient form for ease of numerical implementations. The solutions provide means for rapid computations of the sound fields above and below the rigid porous ground. The improved asymptotic formulas for both situations agree well with numerical results obtained by other numerical schemes, which are more accurate but computationally more intensive. More importantly, the asymptotic solutions can be written in the well-known form of the Weyl-van der Pol formula, which provides a direct correlation between the reflected wave term for the sound field above the porous ground and the transmitted (refracted) wave term for the sound field below.     

Analytic approximations of statistical quantities and response of noisy oscillators

Nonlinear oscillators have been utilized in many contexts because they encompass a large class of phenomena. For a reduced phase oscillator model with weak noise forcing that is necessarily multiplicative, we provide analytic formulas for the stationary statistical quantities of the random period. This is an important quantity which we term ‘response’ (i.e., the spike times, instantaneous frequency in neuroscience, the cycle time in chemical reactions, etc.) that is often analytically intractable in noisy oscillator systems. The analytic formulas are accurate in the weak noise limit so that one does not have to numerically solve a time-varying Fokker-Planck equation. The steady-state and dynamic responses are also analyzed with deterministic forcing. A second order analytic formula is derived for the steady-state response, whereas the dynamic response with time-varying forcing is more complicated. We focus on the specific case where the forcing is sinusoidal and accurately capture the frequency response with an analytic approximation that is obtained with a rescaling of the equation. By utilizing various techniques in the weak noise regime, this work leads to a better understanding of how the random period of nonlinear oscillators are affected by multiplicative noise and external forcing. Comparisons of the asymptotic formulas with a full oscillator system confirm the qualitative accurateness of the theory.     

Behavior near the focal points of asymptotic solutions to the Cauchy problem for the linearized shallow water equations with initial localized perturbations

We study the behavior of the wave part of asymptotic solutions to the Cauchy problem for linearized shallow water equations with initial perturbations localized near the origin. The global representation for these solutions based on the generalized Maslov canonical operator was given earlier. The asymptotic solutions are also localized in the neighborhood of certain curves (fronts). The simplification of general formulas and the behavior of asymptotic solutions in a neighborhood of the regular part of fronts was also given earlier. Here the behavior of asymptotic solutions in a neighborhood of the focal point of the fronts is discussed in detail and the proof of formulas announced earlier for the wave equation is given. This paper can be regarded as a continuation of the paper in Russiian Journal of Mathematical Physics 15 (2), 192–221 (2008). In memoriam V.A. Borovikov     

Transmission capability of a sequential relativistic quantum communication channel with limited observation time

It is shown that the identity of particles must inevitably be taken into account, because states in quantum field theory are nonlocalizable. This circumstance, together with finite limiting velocity, is responsible for the asymptotic character of formulas for the transmission capability of nonrelativistic communication channels (they are formally valid only for infinite time delay between messages, when the identity of particles is negligible, and, correspondingly, for infinitely slow transmission in time—bits per second per message). The transmission capability of a sequential relativistic quantum communication channel is obtained in real time with allowance for the identity of particles.     

Asymptotic theory of linear water waves in a domain with nonuniform bottom with rapidly oscillating sections

A linear problem for propagation of gravity waves in the basin having the bottom of a form of a smooth background with added rapid oscillations is considered. The formulas derived below are asymptotic ones; they are quite formal, and we do not discuss the problem concerning their uniformness with respect to these parameters.     

An asymptotic approach to the connection problem for the first and the second painlevé equations

A new asymptotic method of attack on the connection problem around the point at infinity for Painlevé transcendents of the first and second kind is developed. Connection results between two general angles of approach to infinity, not on special sector boundaries, are obtained for generic four-real-parameter asymptotic behaviours of the Painlevé transcendent. Here the term generic means that two of the four free real parameters are constrained to be nonzero.     

Oscillatory and Asymptotic Behavior of a Second-Order Nonlinear Functional Differential Equations

This paper is concerned with oscillatory and asymptotic behavior of solutions of a class of second order nonlinear functional differential equations. By using the generalized Riccati transformation and the integral averaging technique, new oscillation criteria and asymptotic behavior are obtained for all solutions of the equation. Our results generalize and improve some known theorems.     

The asymptotic behavior of the eigenvalues and eigenfunctions of the Laplace operator on a compact Riemannian manifold

The asymptotic properties of the spectrum of a Laplace operator on a Riemannian manifold are studied. New asymptotic formulas are derived for spectrum series, which are associated with stable geodesics.Translated from Izvestiya Vysshikh Uchebnykh Zavedenii, Radiofizika, Vol. 15, No. 3, pp. 448–453, March, 1972.The authors thank S. I. Al'ber for the statement of the problem and his interest in the work.