Internal gravity wave beams as invariant solutions of Boussinesq equations in geophysical fluid dynamics

It is shown that Lie group analysis of differential equations provides the exact solutions of two-dimensional stratified rotating Boussinesq equations which are a basic model in geophysical fluid dynamics. The exact solutions are obtained as group invariant solutions corresponding to the translation and dilation generators of the group of transformations admitted by the equations. The comparison with the previous analytic studies and experimental observations confirms that the anisotropic nature of the wave motion allows to associate these invariant solutions with uni-directional internal wave beams propagating through the medium. It is also shown that the direction of internal wave beam propagation is in the transverse direction to one of the invariants which corresponds to a linear combination of the translation symmetries. Furthermore, the amplitudes of a linear superposition of wave-like invariant solutions forming the internal gravity wave beams are arbitrary functions of that invariant. Analytic examples of the latitude-dependent invariant solutions associated with internal gravity wave beams that have different general profiles along the obtained invariant and propagating in the transverse direction are considered. The behavior of the invariant solutions near the critical latitude is illustrated.     

Dynamics of a rotating layer of an ideal electrically conducting incompressible fluid

A system of nonlinear partial differential equations is considered that models perturbations in a layer of an ideal electrically conducting rotating fluid bounded by spatially and temporally varying surfaces with allowance for inertial forces. The system is reduced to a scalar equation. The solvability of initial boundary value problems arising in the theory of waves in conducting rotating fluids can be established by analyzing this equation. Solutions to the scalar equation are constructed that describe small-amplitude wave propagation in an infinite horizontal layer and a long narrow channel.     

A new scheme for the solution of reaction diffusion and wave propagation problems

In this paper, a robust numerical scheme is presented for the reaction diffusion and wave propagation problems. The present method is rather simple and straightforward. The Houbolt method is applied so as to convert both types of partial differential equations into an equivalent system of modified Helmholtz equations. The method of fundamental solutions is then combined with the method of particular solution to solve these new systems of equations. Next, based on the exponential decay of the fundamental solution to the modified Helmholtz equation, the dense matrix is converted into an equivalent sparse matrix. Finally, verification studies on the sensitivity of the method’s accuracy on the degree of sparseness and on the time step magnitude of the Houbolt method are carried out for four benchmark problems.     

The Well-Posedness of Dynamical Equations of Magneto-electro-elasticity

A mathematical model of wave propagation in magneto-electro-elastic materials is obtained in the form of a symmetric hyperbolic system of the first-order partial differential equations. This model is a result of the qualitative analysis of the coupled time-dependent Maxwell’s equations and equations of elastodynamics which are considered together with constitutive relations in non-homogeneous anisotropic magneto-electro-elastic materials. Applying the theory and methods of symmetric hyperbolic systems, we have proved that the reported model of wave propagation in magneto-electro-elastic materials satisfies the Hadamards correctness requirements: solvability, uniqueness and stability with respect to perturbation of data.     

Reconstructing parameters of a medium from incomplete spectral information

The problem of the synthesis of a stratified medium with specified amplitude and phase properties is investigated. The wave propagation in the medium is described by a system of differential equations. The synthesis problem considered in the paper relates to inverse problems of spectral analysis with incomplete spectral information. Using the contour integral method we study properties of spectral characteristics and obtain algorithms for the solution of the synthesis problem for differential equations with singularities.     

3D dynamic responses of a multi-layered transversely isotropic saturated half-space under concentrated forces and pore pressure

Dynamic Green's function plays an important role in the study of various wave radiation, scattering and soil-structure interaction problems. However, little research has been done on the response of transversely isotropic saturated layered media. In this paper, the 3D dynamic responses of a multi-layered transversely isotropic saturated half-space subjected to concentrated forces and pore pressure are investigated. First, utilizing Fourier expansion in circumferential direction accompanied by Hankel integral transform in radial direction, the wave equations for transversely isotropic saturated medium in cylindrical coordinate system are solved. Next, with the aid of the exact dynamic stiffness matrix for in-plane and out-of-plane motions, the solutions for multi-layered transversely isotropic saturated half-space under concentrated forces and pore pressure are obtained by direct stiffness method. A FORTRAN computer code is developed to achieve numerical evaluation of the proposed method, and its accuracy is validated through comparison with existing solutions that are special cases of the more general problems addressed. In addition, selected numerical results for a homogeneous and a layered material model are performed to illustrate the effects of material anisotropy, load frequency, drainage condition and layering on the dynamic responses. The presented solutions form a complete set of Green's functions for concentrated forces (including horizontal load in x(y)-direction, vertical load in z-direction) as well as pore pressure, which lays the foundation for further exploring wave propagation of complex local site in a layered transversely isotropic saturated half-space by using the BEMs.     

Propagation of TM waves in a layer with arbitrary nonlinearity

A boundary value problem for Maxwell’s equations describing propagation of TM waves in a nonlinear dielectric layer with arbitrary nonlinearity is considered. The layer is located between two linear semi-infinite media. The problem is reduced to a nonlinear boundary eigenvalue problem for a system of second-order nonlinear ordinary differential equations. A dispersion equation for the eigenvalues of the problem (propagation constants) is derived. For a given nonlinearity function, the dispersion equation can be studied both analytically and numerically. A sufficient condition for the existence of at least one eigenvalue is formulated.     

Optimal control problems for wave equations

The problem on minimizing a quadratic functional on trajectories of the wave equation is considered. We assume that the density of external forces is a control function. A control problem for a partial differential equation is reduced to a control problem for a countable system of ordinary differential equations by use of the Fourier method. The controllability problem for this countable system is considered. Conditions of the noncontrollability for some wave equations were obtained.     

EXACT EXPLICIT SOLUTIONS OF THE NONLINEAR SCHRDINGER EQUATION COUPLED TO THE BOUSSINESQ EQUATION

A system comprised of the nonlinear Schrodinger equation coupled to the Boussinesq equation (S-B equations) which dealing with the stationary propagation of coupled non-linear upper-hybrid and magnetosonic waves in magnetized plasma is proposed. To examine its solitary wave solutions, a reduced set of ordinary differential equations are considered by a simple traveling wave transformation. It is then shown that several new     

Formalism of two potentials for the numerical solution of Maxwell’s equations

A new formulation of Maxwell’s equations based on the introduction of two vector and two scalar potentials is proposed. As a result, the electromagnetic field equations are written as a hyperbolic system that contains, in contrast to the original Maxwell system, only evolution equations and does not involve equations in the form of differential constraints. This makes the new equations especially convenient for the numerical simulation of electromagnetic processes. Specifically, they can be solved by applying powerful modern shock-capturing methods based on the approximation of spatial derivatives by upwind differences. The cases of an electromagnetic field in a vacuum and an inhomogeneous material are considered. Examples are given in which electromagnetic wave propagation is simulated by solving the formulated system of equations with the help of modern high-order accurate schemes.     

Travelling wave solutions for some models in phytopathology

In this work several models of fungal disease propagation are considered. They consist of reaction-diffusion systems coupled with ordinary differential equations with or without time delay as well as integro-differential system of equations. We derive some conditions that ensure the existence and uniqueness of travelling wave solutions for these various models. Our proof is based on a suitable re-formulation in the form of a nonlinear integral equation with measure kernel convolutions.     

Propagation of harmonic waves in prestressed fiber viscoelastic thick tubes filled with a viscous dusty fluid

In this work, propagation of harmonic waves in initially stressed cylindrical viscoelastic thick tubes filled with a Newtonian fluid is studied. The tube, subjected to a static inner pressure P_i and a positive axial stretch λ, will be considered as an incompressible viscoelastic and fibrous material. The fluid is assumed as an incompressible, viscous and dusty fluid. The field equations for the fluid are obtained in the cylindrical coordinates. The governing differential equations of the tube’s viscoelastic material are obtained also in the cylindrical coordinates utilizing the theory of small deformations superimposed on large initial static deformations. For the axially symmetric motion the field equations are solved by assuming harmonic wave solutions. A closed form solution can be obtained for equations governing the fluid body, but due to the variability of the coefficients of resulting differential equations of the solid body, such a closed form solution is not possible to obtain. For that reason, equations for the solid body and the boundary conditions are treated numerically by the finite-difference method to obtain the effects of the thickness of the tube on the wave characteristics. Dispersion relation is obtained using the long wave approximation and, the wave velocities and the transmission coefficients are computed.     

Volume integral equations of the scattering problem of poroelasticity and their properties

Scattering of monochromatic waves on an isolated inhomogeneity (inclusion) in an infinite poroelastic medium is considered. Wave propagation in the medium and the inclusion are described by Biot's equations of poroelasticity. The problem is reduced to 3D‐integro‐differential equations for displacement and pressure fields in the region occupied by the inclusion. Properties of the integral operators in these equations are studied. Discontinuities of the fields on the inclusion boundary are indicated. The case of a thin inclusion with low permeability is considered. The corresponding scattering problem is reduced to a 2D integral equation on the middle surface of the inclusion. The unknown function in this equation is the pressure jump in the transverse direction to the inclusion middle surface. An inclusion with a thin layer of low permeability on its interface is considered. The appropriate boundary conditions on the inclusion interface are pointed out. Methods of numerical solution of the volume integral equations of the scattering problems of poroelasticity are discussed.     

Maxwell’s equations in inhomogeneous bi-anisotropic materials: Existence,uniqueness and stability for the initial value problem

In the present paper inhomogeneous bi-anisotropic materials characterized by matrices of electric permittivity, magnetic permeability and magnetoelectric characteristics are considered. All elements of these matrices are functions of the position in three dimensional space. The time-dependent Maxwell’s equations describe the electromagnetic wave propagation in these materials. Maxwell’s equations together with zero initial data are analyzed and a statement of the initial value problem (IVP) is formulated. This IVP is reduced to the IVP for a symmetric hyperbolic system of partial differential equations of the first order. Applying the theory of a symmetric hyperbolic system, new existence, uniqueness and stability estimate theorems have been obtained for the IVP of Maxwell’s equations in inhomogeneous bi-anisotropic materials.     

Modelling of circumferential waves in cylindrical thermoelastic plates with voids

The propagation of thermoelastic waves along circumferential direction in homogeneous, isotropic, cylindrical curved solid plates with voids has been investigated in the context of linear generalized theory of thermoelasticity. The plate is subjected to stress free or rigidly fixed, thermally insulated or isothermal boundary conditions. Mathematical modeling of the problem for the considered cylindrical curved plate with voids leads to a system of coupled partial differential equations. The model has been simplified by using the Helmholtz decomposition technique and the resulting equations are solved by using the method of separation of variables. The formal solution obtained by using Bessel’s functions with complex arguments is utilized to derive the secular equations which govern the wave motion in the plate with voids. The longitudinal shear motion and axially symmetric shear vibration modes get decoupled from the rest of the motion in contrast to non-axially symmetric plane strain vibrations. These modes remain unaffected due to thermal variations and presence of voids. In order to illustrate theoretical developments, numerical solutions have been carried out for a stress free, thermally insulated or isothermal magnesium plate and are presented graphically. The obtained results are also compared with those available in the literature.     

A finite difference technique for solving optimization problems governed by linear functional differential equations

Aspects of the approximation and optimal control of systems governed by linear retarded nonautonomous functional differential equations (FDE) are considered. First, certain FDE are shown to be equivalent to corresponding abstract ordinary differential equations (ODE). Next, it is demonstrated that these abstract ODE may be approximated by difference equations in finite dimensional spaces. The optimal control problem for systems governed by FDE is then reduced to a sequence of mathematical programming problems. Finally, numerical results for two examples are presented and discussed.     

Numerical Solution of Systems of Loaded Ordinary Differential Equations with Multipoint Conditions

A system of loaded ordinary differential equations with multipoint conditions is considered. The problem under study is reduced to an equivalent boundary value problem for a system of ordinary differential equations with parameters. A system of linear algebraic equations for the parameters is constructed using the matrices of the loaded terms and the multipoint condition. The conditions for the unique solvability and well-posedness of the original problem are established in terms of the matrix made up of the coefficients of the system of linear algebraic equations. The coefficients and the righthand side of the constructed system are determined by solving Cauchy problems for linear ordinary differential equations. The solutions of the system are found in terms of the values of the desired function at the initial points of subintervals. The parametrization method is numerically implemented using the fourth-order accurate Runge–Kutta method as applied to the Cauchy problems for ordinary differential equations. The performance of the constructed numerical algorithms is illustrated by examples.     

Stabilization of the uniform Timoshenko beam by one locally distributed feedback

In this article, we study the energy decay rate for an elastic Timoshenko system. This system consists of two coupled wave equations. Only the equation about the rotation angle is damped by one locally distributed feedback at the neighbourhood of the boundary. The equation for the transverse displacement of the beam is only indirectly damped through the coupling. First, we establish an exponential energy decay rate in the case of the same speed of propagation. Next, when the wave speeds are different, a polynomial-type decay rate is obtained. These results are proved by verifying the frequency domain conditions.     

A mathematical model for the three-dimentional ocean sound propagation

A model is developed mathematically to represent sound propagation in a three-dimensional ocean. The complete development is based on characteristics of the physical environment, mathematical theory, and computational accuracy.While the two-dimentional underwater acoustic wave propagation problem is not yet solved completely for range-dependent environments,three-dimentional environmental effects, such as fronts and eddies, often cannot be neglected. To predict underwater sound propagation, one usually deals with the solution of the Helmholtz (reduced wave) equation. This elliptical equation, along with a set of boundary conditions including a wall condition at the maximum range, forms a well-posed problem, which is pure boundary-value problem. An existing approach to economically solve this three-dimensional range-dependent problem is by means of a two-dimensional parabolic partial differential equation. This parabolic approximation approach, within the limitation of mathematical and acoustical approximations, offers efficient solutions to a class of long-range propagation problems. The parabolic wave equation is much easier to solve than the elliptic equation; one major saving is the removal of the wall boundary condition at the maximum range. The application of the two-dimensional parabolic wave equation to a number of realistic problems has been successful.We discuss the extension of the parabolic equation approach to three-dimensional problems. This paper begins with general considerations of the three-dimensional elliptic wave equation and shows how to transform this equation into parabolic equations which are easier to solve. The development of this paper focuses on wide angle three-dimensional underwater acoustic propagation and accommodates as a special case prevoius developments by other authors. In the course of our development, the physical properties, mathematical validity, and computational accuracy are the primary factors considered. We describe how parabolic wave equations are derived and how wide angle propagation is taken into consideration. Then, a discussion of the limitations and the advantages of the parabolic equation approximation is highlighted. These provide the background for the mathematical formulation of three-dimensional underwater acoustic wave propagation models.Modelling the mathematical solution to three-dimensional underwater acoustic wave propagation involves difficulties both in describing the theoretical acoustics and in performing the large scale computations. We have used the mathematical and physical properties of the problem to simplify considerably. Simplications allow us to introduce a three-dimensional mathematical model for underwater acoustic propagation predictions. Our wide angle three-dimensional parabolic equation model is theoretically justifiable and computationally accurate. This model offers a variety of capabilities to handle a class of long-range propagation problems under acoustical environments with three-dimensional variations.