Algebraic and Geometric Properties of Matrix Solutions of Nonlinear Wave Equations

We improve the construction of exact matrix solutions for nonlinear wave equations by using unitary anti-Hermitian and anticommuting matrices. We prove the theorem that constructs the matrix functions u _n satisfying the nonlinear wave equation for a set of special potentials. In this case, the graph of complex solution u ₁ has a soliton-like form with a finite number of coils. Exponential representation of matrix solutions u _n is associated with continuous rotations that can be used for describing intrinsic rotations and state changes of elementary particles. We also prove the theorem on the decomposition of continuous rotation (described by solution u ₂) onto three simultaneous rotations about coordinate vectors. Each of the three constructed matrix solutions u ₃ is also decomposed into the triplet of elementary matrix solutions.     

Equations for the description of wave scattering by multi-valued random surfaces

We consider a scattering theory for multi-valued rough surfaces which cannot be described by the conventional equation of the type z = ζ(x,y). Both Dirichlet and Neumann problems are analyzed. Starting with Green's theorem we obtain a representation of the scattered field, the surface integral equation, and the extinction theorem for such surfaces. In contrast to conventional theory, these equations contain three random functions x = x(u ₁,u ₂), y = y(u ₁,u ₂), and z = z(u ₁,u ₂), where u ₁ and u ₂ are the parameters describing the surface. We introduce two scattering amplitudes S _± for describing the scattered wave above and below the surface. The extinction theorem, if formulated in terms of S _?, allows us to determine S _? for an arbitrary multi-valued surface and after this it becomes possible to derive a simple integral equation for surface sources. Knowledge of the surface sources allows us to find S ₊ by integration.     

Long-Time Asymptotics for Strong Solutions¶of the Thin Film Equation

In this paper we investigate the large-time behavior of strong solutions to the one-dimensional fourth order degenerate parabolic equation u _t =−(u u _xxx ) _x , modeling the evolution of the interface of a spreading droplet. For nonnegative initial values u ₀(x)∈H ¹(ℝ), both compactly supported or of finite second moment, we prove explicit and universal algebraic decay in the L ¹-norm of the strong solution u(x,t) towards the unique (among source type solutions) strong source type solution of the equation with the same mass. The method we use is based on the study of the time decay of the entropy introduced in [13] for the porous medium equation, and uses analogies between the thin film equation and the porous medium equation. Received: 2 February 2001 / Accepted: 7 October 2001     

Dynamical Localization of Quantum Walks in Random Environments

We study shock statistics in the scalar conservation law ∂ _t u+∂ _x f(u)=0, x∈ℝ, t>0, with a convex flux f and spatially random initial data. We show that the Markov property (in x) is preserved for a large class of random initial data (Markov processes with downward jumps and derivatives of Lévy processes with downward jumps). The kinetics of shock clustering is then described completely by an evolution equation for the generator of the Markov process u(x,t), x∈ℝ. We present four distinct derivations for this evolution equation, and show that it takes the form of a Lax pair. The Lax equation admits a spectral parameter as in Manakov (Funct. Anal. Appl. 10:328–329, 1976), and has remarkable exact solutions for Burgers equation (f(u)=u ²/2). This suggests the kinetic equations of shock clustering are completely integrable.     

Equations with Singular Diffusivity

Recently models of faceted crystal growth and of grain boundaries were proposed based on the gradient system with nondifferentiable energy. In this article, we study their most basic forms given by the equations u _t=(u _x/|u _x|) _x and u _t=(1/a)(au _x/|u _x|) _x , where both of the related energies include a |u _x| term of power one which is nondifferentiable at u _x=0. The first equation is spatially homogeneous, while the second one is spatially inhomogeneous when a depends on x. These equations naturally express nonlocal interactions through their singular diffusivities (infinitely large diffusion constant), which make the profiles of the solutions completely flat. The mathematical basis for justifying and analyzing these equations is explained, and theoretical and numerical approaches show how the solutions of the equations evolve.     

The lie algebra of invariant group of the KdV,MKdV, or Burgers equation

Let (E): u _t=H(u) denote the KdV, MKdV or Burgers equation, and U(s)=(D^j _s)/u _j, where D=d/dx, u _i=Dⁱ _u, s=s(u, u ₁, ..., u _n) is a polynomial of u _i with constant coefficients, be the generator of invariant group of equation (E). We prove in this paper that all such generators form a commutative Lie algebra, from which it follows that for any symmetry s(u, ..., u _n) of (E), the evolution equation u _t=s(u, ..., u _n) possesses an infinite number of symmetries (or conservation laws in the case of KdV and MKdV equations).     

Zeitliche Entwicklung des Reflexionskoeffizienten einer Plasmaschicht bei starker elektromagnetischer Einstrahlung

We have calculated numerically the temporal evolution of the nonlinear reflection coefficient R of an overdense plasma layer by solving the system of partial differential equations consisting of the wave equation for the slowly varying amplitude of the electric field and the hydrodynamic equations for the ion motion, including a ponderomotive force term. In dependence of the (normalized) amplitude of the incident wave u_A two regimes exist: Below a critical amplitude u_A^* ? 1 the reflection coefficient is approximately independent on the amplitude u_A and temporally constant. In the opposite case u_A>u_A^*, on the other hand, ‖R‖² decreases slowly with time down to a minimum value and after that it increases rapidly to the initial value. We think, that our results are important to interpret the anomalous reflectivity observed in some experiments when strong electromagnetic waves are incident on an overdense plasma.     

On Integration of the Nonlinear d’Alembert-Eikonal System and Conditional Symmetry of Nonlinear Wave Equations

Abstract

We study integrability of a system of nonlinear partial differential equations consisting of the nonlinear d’Alembert equation □u = F (u) and nonlinear eikonal equation u _xµ u _x µ = G(u) in the complex Minkowski space R(1, 3). A method suggested makes it possible to establish necessary and sufficient compatibility conditions and construct a general solution of the d’Alembert-eikonal system for all cases when it is compatible. The results obtained can be applied, in particular, to construct principally new (non-Lie, non-similarity) solutions of the non-linear d’Alembert, Dirac, and Yang-Mills equations. Solutions found in this way are shown to correspond to conditional symmetry of the equations enumerated above. Using the said approach, we study in detail conditional symmetry of the nonlinear wave equation □w = F ₀(w) in the four-dimensional Minkowski space. A number of new (non-Lie) reductions of the above equation are obtained giving rise to its new exact solutions which contain arbitrary functions.     

The field-theoretic description of dynamics of interfaces in disordered media

The time evolution of an interface in a disordered media is described by using the propagator method. The method enables one to represent the perturbation expansions of different quantities characterizing the interface by means of diagrams which are familiar from the field theory. By the analysis of the divergences in the vicinity of the critical dimension d_c = 4 we found that the regularization of the theory demands the renormalization of the mobility and all moments of the disorder correlator excepting the zero one. The renormalization group (RG) calculations of the average velocity of the interface, the roughness, and the functional RG equation of the disorder correlator are presented to order ? = 4 - d. The latter coincides with the result obtained by D. S. Fisher in the equilibrium case. The RG equations have a pole at the value of the driving force, which coincides with the value of the threshold below which the interface becomes pinned as predicted by Bruinsma and Aeppli. The behavior of the mobility in the vicinity of the pole is discussed.     

Variational Operators,Symplectic Operators,and the Cohomology of Scalar Evolution Equations

For a scalar evolution equation u_t = K(t, x, u, u_x, . . . , u_2m+1) with m ≥ 1, the cohomology space H^1,2() is shown to be isomorphic to the space of variational operators and an explicit isomorphism is given. The space of symplectic operators for u_t = K for which the equation is Hamiltonian is also shown to be isomorphic to the space H^1,2() and subsequently can be naturally identified with the space of variational operators. Third order scalar evolution equations admitting a first order symplectic (or variational) operator are characterized. The variational operator (or symplectic) nature of the potential form of a bi-Hamiltonian evolution equation is also presented in order to generate examples of interest.     

On the structure of exceptional evolution polynomials

Given an evolution equation u _t =P(u, u_x,..., u_M), P polynomial, a set of obstructions is derived whose fulfillment is necessary for the existence of (nontrivial polynomial) conserved densities p(u,..., u_N) for infinitely many values of N and some nontrivial infinitesimal symmetry.     

Renormalization-Group Method for Reduction of Evolution Equations; Invariant Manifolds and Envelopes

The renormalization group (RG) method as a powerful tool for reduction of evolution equations is formulated in terms of the notion of invariant manifolds. We start with derivation of an exact RG equation which is analogous to the Wilsonian RG equations in statistical physics and quantum field theory. It is clarified that the perturbative RG method constructs invariant manifolds successively as the initial value of evolution equations, thereby the meaning to set t₀=t is naturally understood where t₀ is the arbitrary initial time. We show that the integral constants in the unperturbative solution constitutes natural coordinates of the invariant manifold when the linear operator A in the evolution equation is semi-simple, i.e., diagonalizable; when A is not semi-simple and has a Jordan cell, a slight modification is necessary because the dimension of the invariant manifold is increased by the perturbation. The RG equation determines the slow motion of the would-be integral constants in the unperturbative solution on the invariant manifold. We present the mechanical procedure to construct the perturbative solutions hence the initial values with which the RG equation gives meaningful results. The underlying structure of the reduction by the RG method as formulated in the present work turns out to completely fit to the universal one elucidated by Kuramoto some years ago. We indicate that the reduction procedure of evolution equations has a good correspondence with the renormalization procedure in quantum field theory; the counter part of the universal structure of reduction elucidated by Kuramoto may be Polchinski's theorem for renormalizable field theories. We apply the method to interface dynamics such as kink–anti-kink and soliton–soliton interactions in the latter of which a linear operator having a Jordan-cell structure appears.     

一类非线性演化方程新的显式行波解

借助Mathematica软件,采用三角函数法和吴文俊消元法,获得了一类非线性演化方程u_tt+au_xx+bu+cu³=0的三组行波解,其中包括新的行波解、扭状孤波解和钟状孤波解.从而作为该方程的特例,如Duffing方,Klein-Gordon方程、Landau-Ginburg-Higgs方程和⁴方程等也都获得了相应的若干行波解.这种方法也适用于其他非线性方程. 关键词：     

On Properties of Hamiltonian Structures for a Class of Evolutionary PDEs

In a recent paper we proved that for certain class of perturbations of the hyperbolic equation u _t = f (u)u _x , there exist changes of coordinate, called quasi-Miura transformations, that reduce the perturbed equations to the unperturbed one. We prove in the present paper that if in addition the perturbed equations possess Hamiltonian structures of certain type, the same quasi-Miura transformations also reduce the Hamiltonian structures to their leading terms. By applying this result, we obtain a criterion of the existence of Hamiltonian structures for a class of scalar evolutionary PDEs and an algorithm to find out the Hamiltonian structures.      

Solitary waves in the Madelung's fluid: Connection between the nonlinear Schrödinger equation and the Korteweg-de Vries equation

An investigation to deepen the connection between the family of nonlinear Schr?dinger equations and the one of Korteweg-de Vries equations is carried out within the context of the Madelung's fluid picture. In particular, under suitable hypothesis for the current velocity, it is proven that the cubic nonlinear Schr?dinger equation, whose solution is a complex wave function, can be put in correspondence with the standard Korteweg-de Vries equation, is such a way that the soliton solutions of the latter are the squared modulus of the envelope soliton solution of the former. Under suitable physical hypothesis for the current velocity, this correspondence allows us to find envelope soliton solutions of the cubic nonlinear Schr?dinger equation, starting from the soliton solutions of the associated Korteweg-de Vries equation. In particular, in the case of constant current velocities, the solitary waves have the amplitude independent of the envelope velocity (which coincides with the constant current velocity). They are bright or dark envelope solitons and have a phase linearly depending both on space and on time coordinates. In the case of an arbitrarily large stationary-profile perturbation of the current velocity, envelope solitons are grey or dark and they relate the velocity u₀ with the amplitude; in fact, they exist for a limited range of velocities and have a phase nonlinearly depending on the combined variable x-u_{0 s} (s being a time-like variable). This novel method in solving the nonlinear Schr?dinger equation starting from the Korteweg-de Vries equation give new insights and represents an alternative key of reading of the dark/grey envelope solitons based on the fluid language. Moreover, a comparison between the solutions found in the present paper and the ones already known in literature is also presented. Received 20 February 2002 and Received in final form 22 April 2002 Published online 6 June 2002     

Symmetry groups and Gauss kernels of Schrödinger equations

The relationship between symmetries and Gauss kernels for the Schrödinger equation iu_t=u_xx+f(x)u is established. It is shown that if the Lie point symmetries of the equation are nontrivial, a classical integral transformations of the Gauss kernels can be obtained. Then the Gauss kernels of Schrödinger equations are derived by inverting the integral transformations. Furthermore, the relationship between Gauss kernels for two equations related by an equivalence transformation is identified.     

A higher order lattice BGK model for simulating some nonlinear partial differential equations

In this paper, we consider a one-dimensional nonlinear partial differential equation that has the form ut + αuux + βunux - γuxx + δuxxx = F(u). A higher order lattice Bhatnager-Gross-Krook (BGK) model with an amending-function is proposed. With the Chapman-Enskog expansion, different kinds of nonlinear partial differential equations are recovered correctly from the continuous Boltzmann equation. The numerical results show that this method is very effective.     

Nonlinear kinematic response of premixed flames to harmonic velocity disturbances

This paper analyzes the nonlinear dynamics of premixed flames responding to harmonic velocity disturbances. These nonlinear dynamics were studied by solving a constant flame speed front tracking equation for the flame’s response to harmonically oscillating velocity disturbances. The solution to these equations is used to quantify the transfer function relating the ratio of the normalized flame area to velocity fluctuations, G = (A′/A_o)/(u′/u_o), upon the amplitude of velocity oscillations, ε = u′/u_o. Due to nonlinearities, the amplitude of this transfer function relative to its linear value decreases with increasing amplitude of velocity oscillation, u′/u_o. In contrast, the transfer function phase exhibits almost no amplitude dependence. The velocity amplitude where transfer function nonlinearities become significant depends strongly upon three parameters: a Strouhal number, St = ωL_f/u_o (where L_f is the flame length), the ratio of the flame length to width, β = L_f/R, and the flame shape in the absence of perturbations (i.e., conical, inverted wedge, etc.). In the linear case, the transfer function, G, depends only upon an algebraic combination of the first two parameters, given by St₂ = St (1 + β²)/β². In general, however, G exhibits a distinct dependence upon both parameters St and β. In particular, we show that the nonlinear response of G is an intrinsically dynamic phenomenon; i.e., its quasi-steady response (St 1) is purely linear. As such, nonlinearity is enhanced with increasing Strouhal numbers. In contrast, nonlinearity is suppressed at large β values; as such, the response of a long flame remains quite similar to its linear value, even at large ε values where the flame front exhibits substantial corrugation and cusping. Finally, we show that the response of conical flames remains much more linear at comparable disturbance amplitudes than for “V” or wedge-shaped flames. These predictions are shown to be consistent with available experimental data.