Wave propagation in spatially modulated tubes

Energy Technology Data Exchange (ETDEWEB)

Ziepke, A., E-mail: ziepke@itp.tu-berlin.de; Martens, S.; Engel, H. [Institut für Theoretische Physik, Hardenbergstraße 36, EW 7-1, Technische Universität Berlin, 10623 Berlin (Germany)

2016-09-07

We investigate wave propagation in rotationally symmetric tubes with a periodic spatial modulation of cross section. Using an asymptotic perturbation analysis, the governing quasi-two-dimensional reaction-diffusion equation can be reduced into a one-dimensional reaction-diffusion-advection equation. Assuming a weak perturbation by the advection term and using projection method, in a second step, an equation of motion for traveling waves within such tubes can be derived. Both methods predict properly the nonlinear dependence of the propagation velocity on the ratio of the modulation period of the geometry to the intrinsic width of the front, or pulse. As a main feature, we observe finite intervals of propagation failure of waves induced by the tube’s modulation and derive an analytically tractable condition for their occurrence. For the highly diffusive limit, using the Fick-Jacobs approach, we show that wave velocities within modulated tubes are governed by an effective diffusion coefficient. Furthermore, we discuss the effects of a single bottleneck on the period of pulse trains. We observe period changes by integer fractions dependent on the bottleneck width and the period of the entering pulse train.

On the propagation of low-hybrid waves of finite amplitude

International Nuclear Information System (INIS)

Kozyrev, A.N.; Piliya, A.D.; Fedorov, V.I.

1979-01-01

Propagation of low-hybrid waves of a finite amplitude with allowance for variation in plasma density caused by HF field pressure is studied. Considered is wave ''overturning'' which takes place in the absence of space dispersion. With taking account of dispersion the wave propagation is described by the third-order nonlinear equation which differs in shape from the complex modified Korteweg-de-Vries (Hirota) equation. Solutions of this equation of the space solution type are found

Use of conformal mapping to describe MHD wave propagation

International Nuclear Information System (INIS)

Bulanov, S.V.; Pegoraro, F.

1993-01-01

A method is proposed for finding explicit exact solutions of the magnetohydrodynamic equations describing the propagation of magnetoacoustic waves in a plasma in a magnetic potential that depends on two spatial coordinates. This method is based on the use of conformal mappings to transform the wave equation into an equation describing the propagation of waves in a uniform magnetic field. The basic properties of magnetoacoustic and Alfven waves near the critical points, magnetic separatrices, and in configuration with magnetic islands are discussed. Expressions are found for the dimensionless parameters which determine the relative roles of the plasma pressure, nonlinearity, and dissipation near the critical points. 30 refs

A fast-multipole domain decomposition integral equation solver for characterizing electromagnetic wave propagation in mine environments

KAUST Repository

Yücel, Abdulkadir C.

2013-07-01

Reliable and effective wireless communication and tracking systems in mine environments are key to ensure miners\\' productivity and safety during routine operations and catastrophic events. The design of such systems greatly benefits from simulation tools capable of analyzing electromagnetic (EM) wave propagation in long mine tunnels and large mine galleries. Existing simulation tools for analyzing EM wave propagation in such environments employ modal decompositions (Emslie et. al., IEEE Trans. Antennas Propag., 23, 192-205, 1975), ray-tracing techniques (Zhang, IEEE Tran. Vehic. Tech., 5, 1308-1314, 2003), and full wave methods. Modal approaches and ray-tracing techniques cannot accurately account for the presence of miners and their equipments, as well as wall roughness (especially when the latter is comparable to the wavelength). Full-wave methods do not suffer from such restrictions but require prohibitively large computational resources. To partially alleviate this computational burden, a 2D integral equation-based domain decomposition technique has recently been proposed (Bakir et. al., in Proc. IEEE Int. Symp. APS, 1-2, 8-14 July 2012). © 2013 IEEE.

Periodic and solitary-wave solutions of the Degasperis-Procesi equation

International Nuclear Information System (INIS)

Vakhnenko, V.O.; Parkes, E.J.

2004-01-01

Travelling-wave solutions of the Degasperis-Procesi equation are investigated. The solutions are characterized by two parameters. For propagation in the positive x-direction, hump-like, inverted loop-like and coshoidal periodic-wave solutions are found; hump-like, inverted loop-like and peakon solitary-wave solutions are obtained as well. For propagation in the negative x-direction, there are solutions which are just the mirror image in the x-axis of the aforementioned solutions. A transformed version of the Degasperis-Procesi equation, which is a generalization of the Vakhnenko equation, is also considered. For propagation in the positive x-direction, hump-like, loop-like, inverted loop-like, bell-like and coshoidal periodic-wave solutions are found; loop-like, inverted loop-like and kink-like solitary-wave solutions are obtained as well. For propagation in the negative x-direction, well-like and inverted coshoidal periodic-wave solutions are found; well-like and inverted peakon solitary-wave solutions are obtained as well. In an appropriate limit, the previously known solutions of the Vakhnenko equation are recovered

Analysis of pulse thermography using similarities between wave and diffusion propagation

Science.gov (United States)

Gershenson, M.

2017-05-01

Pulse thermography or thermal wave imaging are commonly used as nondestructive evaluation (NDE) method. While the technical aspect has evolve with time, theoretical interpretation is lagging. Interpretation is still using curved fitting on a log log scale. A new approach based directly on the governing differential equation is introduced. By using relationships between wave propagation and the diffusive propagation of thermal excitation, it is shown that one can transform from solutions in one type of propagation to the other. The method is based on the similarities between the Laplace transforms of the diffusion equation and the wave equation. For diffusive propagation we have the Laplace variable s to the first power, while for the wave propagation similar equations occur with s2. For discrete time the transformation between the domains is performed by multiplying the temperature data vector by a matrix. The transform is local. The performance of the techniques is tested on synthetic data. The application of common back projection techniques used in the processing of wave data is also demonstrated. The combined use of the transform and back projection makes it possible to improve both depth and lateral resolution of transient thermography.

Electromagnetic Wave Propagation in Random Media

DEFF Research Database (Denmark)

Pécseli, Hans

1984-01-01

The propagation of a narrow frequency band beam of electromagnetic waves in a medium with randomly varying index of refraction is considered. A novel formulation of the governing equation is proposed. An equation for the average Green function (or transition probability) can then be derived...

Separate P‐ and SV‐wave equations for VTI media

KAUST Repository

Pestana, Reynam C.; Ursin, Bjø rn; Stoffa, Paul L.

2011-01-01

In isotropic media we use the scalar acoustic wave equation to perform reverse time migration RTM of the recorded pressure wavefleld data. In anisotropic media P- and SV-waves are coupled and the elastic wave equation should be used for RTM. However, an acoustic anisotropic wave equation is often used instead. This results in significant shear wave energy in both modeling and RTM. To avoid this undesired SV-wave energy, we propose a different approach to separate P- and SV-wave components for vertical transversely isotropic VTI media. We derive independent pseudo-differential wave equations for each mode. The derived equations for P- and SV-waves are stable and reduce to the isotropic case. The equations presented here can be effectively used to model and migrate seismic data in VTI media where ε - δ is small. The SV-wave equation we develop is now well-posed and triplications in the SV wavefront are removed resulting in stable wave propagation. We show modeling and RTM results using the derived pure P-wave mode in complex VTI media and use the rapid expansion method REM to propagate the waveflelds in time. © 2011 Society of Exploration Geophysicists.

Nonlinear propagation of short wavelength drift-Alfven waves

DEFF Research Database (Denmark)

Shukla, P. K.; Pecseli, H. L.; Juul Rasmussen, Jens

1986-01-01

Making use of a kinetic ion and a hydrodynamic electron description together with the Maxwell equation, the authors derive a set of nonlinear equations which governs the dynamics of short wavelength ion drift-Alfven waves. It is shown that the nonlinear drift-Alfven waves can propagate as two-dim...

Wave fields in real media wave propagation in anisotropic, anelastic, porous and electromagnetic media

CERN Document Server

Carcione, José M

2007-01-01

This book examines the differences between an ideal and a real description of wave propagation, where ideal means an elastic (lossless), isotropic and single-phase medium, and real means an anelastic, anisotropic and multi-phase medium. The analysis starts by introducing the relevant stress-strain relation. This relation and the equations of momentum conservation are combined to give the equation of motion. The differential formulation is written in terms of memory variables, and Biot's theory is used to describe wave propagation in porous media. For each rheology, a plane-wave analysis is performed in order to understand the physics of wave propagation. The book contains a review of the main direct numerical methods for solving the equation of motion in the time and space domains. The emphasis is on geophysical applications for seismic exploration, but researchers in the fields of earthquake seismology, rock acoustics, and material science - including many branches of acoustics of fluids and solids - may als...

Study of the Electromagnetic Waves Propagation over the Improved Fractal Sea Surface Based on Parabolic Equation Method

Directory of Open Access Journals (Sweden)

Wenwan Ding

2016-01-01

Full Text Available An improved fractal sea surface model, which can describe the capillary waves very well, is introduced to simulate the one-dimension rough sea surface. In this model, the propagation of electromagnetic waves (EWs is computed by the parabolic equation (PE method using the finite-difference (FD algorithm. The numerical simulation results of the introduced model are compared with those of the Miller-Brown model and the Elfouhaily spectrum inversion model. It has been shown that the effects of the fine structure of the sea surface on the EWs propagation in the introduced model are more apparent than those in the other two models.

Analysis of stress wave propagation in an elasto-viscoplastic plate

International Nuclear Information System (INIS)

Nakagawa, Noritoshi; Kawai, Ryoji; Urushi, Norio.

1986-01-01

Stress waves which propagate in the body are reflected at the boundary, and due to the interaction of the reflected stress waves, the focussing of stress waves will take place and a high stress level can be caused. The focussing of stress waves due to the reflection from the boundary may bring about fracture of the body, so that this is an important problem from a viewpoint of dynamic strength of structures. In this paper the process of stress wave focussing and the strain-rate dependence of constitutive equation in elastic and plastic regions are investigated. In the case where an in-plane step load uniformly acts on the straight edge of the plate with a semi-circular boundary, the propagation of stress waves in the plate was numerically analyzed by the finite element method, applying viscoelastic, elasto-plastic and elasto-viscoplastic constitutive equations. As the result, the process of focussing of stress waves due to reflection from the semi-circular boundary was observed and the difference in propagation behaviour of stress waves was discussed in materials represented by some kinds of constitutive equations. (author)

Propagation of nonlinear shock waves for the generalised Oskolkov equation and its dynamic motions in the presence of an external periodic perturbation

Science.gov (United States)

Ak, Turgut; Aydemir, Tugba; Saha, Asit; Kara, Abdul Hamid

2018-06-01

Propagation of nonlinear shock waves for the generalised Oskolkov equation and dynamic motions of the perturbed Oskolkov equation are investigated. Employing the unified method, a collection of exact shock wave solutions for the generalised Oskolkov equations is presented. Collocation finite element method is applied to the generalised Oskolkov equation for checking the accuracy of the proposed method by two test problems including the motion of shock wave and evolution of waves with Gaussian and undular bore initial conditions. Considering an external periodic perturbation, the dynamic motions of the perturbed generalised Oskolkov equation are studied depending on the system parameters with the help of phase portrait and time series plot. The perturbed generalised Oskolkov equation exhibits period-3, quasiperiodic and chaotic motions for some special values of the system parameters, whereas the generalised Oskolkov equation presents shock waves in the absence of external periodic perturbation.

E3D, 3-D Elastic Seismic Wave Propagation Code

International Nuclear Information System (INIS)

Larsen, S.; Harris, D.; Schultz, C.; Maddix, D.; Bakowsky, T.; Bent, L.

2004-01-01

1 - Description of program or function: E3D is capable of simulating seismic wave propagation in a 3D heterogeneous earth. Seismic waves are initiated by earthquake, explosive, and/or other sources. These waves propagate through a 3D geologic model, and are simulated as synthetic seismograms or other graphical output. 2 - Methods: The software simulates wave propagation by solving the elasto-dynamic formulation of the full wave equation on a staggered grid. The solution scheme is 4-order accurate in space, 2-order accurate in time

Nonlinear wave propagation through a ferromagnet with damping in ...

Indian Academy of Sciences (India)

magnetic waves in a ferromagnet can be reduced to an integro-differential equation. Keywords. Solitons; integro-differential equations; reductive perturbation method. PACS Nos 41.20 Jb; 05.45 Yv; 03.50 De; 78.20 Ls. 1. Introduction. The phenomenon of propagation of electromagnetic waves in ferromagnets are not only.

Efficient techniques for wave-based sound propagation in interactive applications

Science.gov (United States)

Mehra, Ravish

Sound propagation techniques model the effect of the environment on sound waves and predict their behavior from point of emission at the source to the final point of arrival at the listener. Sound is a pressure wave produced by mechanical vibration of a surface that propagates through a medium such as air or water, and the problem of sound propagation can be formulated mathematically as a second-order partial differential equation called the wave equation. Accurate techniques based on solving the wave equation, also called the wave-based techniques, are too expensive computationally and memory-wise. Therefore, these techniques face many challenges in terms of their applicability in interactive applications including sound propagation in large environments, time-varying source and listener directivity, and high simulation cost for mid-frequencies. In this dissertation, we propose a set of efficient wave-based sound propagation techniques that solve these three challenges and enable the use of wave-based sound propagation in interactive applications. Firstly, we propose a novel equivalent source technique for interactive wave-based sound propagation in large scenes spanning hundreds of meters. It is based on the equivalent source theory used for solving radiation and scattering problems in acoustics and electromagnetics. Instead of using a volumetric or surface-based approach, this technique takes an object-centric approach to sound propagation. The proposed equivalent source technique generates realistic acoustic effects and takes orders of magnitude less runtime memory compared to prior wave-based techniques. Secondly, we present an efficient framework for handling time-varying source and listener directivity for interactive wave-based sound propagation. The source directivity is represented as a linear combination of elementary spherical harmonic sources. This spherical harmonic-based representation of source directivity can support analytical, data

The nonlinear Schrödinger equation and the propagation of weakly nonlinear waves in optical fibers and on the water surface

Energy Technology Data Exchange (ETDEWEB)

Chabchoub, A., E-mail: achabchoub@swin.edu.au [Centre for Ocean Engineering Science and Technology, Swinburne University of Technology, Hawthorn, Victoria 3122 (Australia); Kibler, B.; Finot, C.; Millot, G. [Laboratoire Interdisciplinaire Carnot de Bourgogne (ICB), UMR 6303 CNRS, Université de Bourgogne, 21078 Dijon (France); Onorato, M. [Dipartimento di Fisica, Università degli Studi di Torino, Torino 10125 (Italy); Istituto Nazionale di Fisica Nucleare, INFN, Sezione di Torino, Torino 10125 (Italy); Dudley, J.M. [Institut FEMTO-ST, UMR 6174 CNRS- Université de Franche-Comté, 25030 Besançon (France); Babanin, A.V. [Centre for Ocean Engineering Science and Technology, Swinburne University of Technology, Hawthorn, Victoria 3122 (Australia)

2015-10-15

The dynamics of waves in weakly nonlinear dispersive media can be described by the nonlinear Schrödinger equation (NLSE). An important feature of the equation is that it can be derived in a number of different physical contexts; therefore, analogies between different fields, such as for example fiber optics, water waves, plasma waves and Bose–Einstein condensates, can be established. Here, we investigate the similarities between wave propagation in optical Kerr media and water waves. In particular, we discuss the modulation instability (MI) in both media. In analogy to the water wave problem, we derive for Kerr-media the Benjamin–Feir index, i.e. a nondimensional parameter related to the probability of formation of rogue waves in incoherent wave trains.

Effect of surface conditions on blast wave propagation

International Nuclear Information System (INIS)

Song, Seung Ho; Li, Yi Bao; Lee, Chang Hoon; Choi, Jung Il

2016-01-01

We performed numerical simulations of blast wave propagations on surfaces by solving axisymmetric two-dimensional Euler equations. Assuming the initial stage of fireball at the breakaway point after an explosion, we investigated the effect of surface conditions considering surface convex or concave elements and thermal conditions on blast wave propagations near the ground surface. Parametric studies were performed by varying the geometrical factors of the surface element as well as thermal layer characteristics. We found that the peak overpressure near the ground zero was increased due to the surface elements, while modulations of the blast wave propagations were limited within a region for the surface elements. Because of the thermal layer, the precursor was formed in the propagations, which led to the attenuation of the peak overpressure on the ground surface

Van Allen Probe observations of EMIC wave propagation in the inner magnetosphere

Science.gov (United States)

Saikin, A.; Zhang, J.; Smith, C. W.; Spence, H. E.; Torbert, R. B.; Kletzing, C.; Wygant, J. R.

2017-12-01

This study examines the propagation of inner magnetosphere (L vector, , analysis on all observed EMIC wave events to determine the direction of propagation, with bi-directionally propagating EMIC waves indicating the presence of the EMIC wave source region. EMIC waves were considered bi-directional (i.e., in the source region) if at least two wave packets exhibited opposing flux components, and (W/km2), consistently for 60 seconds. Events not observed to have opposing flux components are considered unidirectional. EMIC wave events observed at relatively high magnetic latitudes, generally, are found to propagate away from the magnetic equator (i.e., unidirectional). Bi-directionally propagating EMIC waves are preferably observed at lower magnetic latitudes. The occurrence rate, spatial distribution, and the energy propagation angle of both unidirectionally and bi-directionally propagating EMIC waves are examined with respect to L, MLT, and MLAT.

Simulation and Prediction of Weather Radar Clutter Using a Wave Propagator on High Resolution NWP Data

DEFF Research Database (Denmark)

Benzon, Hans-Henrik; Bovith, Thomas

2008-01-01

for prediction of this type of weather radar clutter is presented. The method uses a wave propagator to identify areas of potential non-standard propagation. The wave propagator uses a three dimensional refractivity field derived from the geophysical parameters: temperature, humidity, and pressure obtained from......Weather radars are essential sensors for observation of precipitation in the troposphere and play a major part in weather forecasting and hydrological modelling. Clutter caused by non-standard wave propagation is a common problem in weather radar applications, and in this paper a method...... a high-resolution Numerical Weather Prediction (NWP) model. The wave propagator is based on the parabolic equation approximation to the electromagnetic wave equation. The parabolic equation is solved using the well-known Fourier split-step method. Finally, the radar clutter prediction technique is used...

Shear wave propagation in piezoelectric-piezoelectric composite layered structure

Directory of Open Access Journals (Sweden)

Anshu Mli Gaur

Full Text Available The propagation behavior of shear wave in piezoelectric composite structure is investigated by two layer model presented in this approach. The composite structure comprises of piezoelectric layers of two different materials bonded alternatively. Dispersion equations are derived for propagation along the direction normal to the layering and in direction of layering. It has been revealed that thickness and elastic constants have significant influence on propagation behavior of shear wave. The phase velocity and wave number is numerically calculated for alternative layer of Polyvinylidene Difluoride (PVDF and Lead Zirconate Titanate (PZT-5H in composite layered structure. The analysis carried out in this paper evaluates the effect of volume fraction on the phase velocity of shear wave.

A theory of coherent propagation of light wave in semiconductors

International Nuclear Information System (INIS)

Zi-zhao, G.; Guo-zhen, Y.

1980-05-01

In this paper, we suggest a theory to describe the pheonmena of coherent propagation of light wave in semiconductors. Basing on two band system and considering the interband and intraband transitions induced by light wave and the interaction between electrons, we obtain the nonlinear equations for the description of interaction between carriers and coherent light wave. We have made use of the equations to analyse the phenomena which arise from the interaction between semiconductors and coherent light, for example, the multiphoton transitions, the saturation of light absorption of exciton, the shift of exciton line in intense light field, and the coherent propagation phenomena such as self-induced transparency, etc. (author)

Longitudinal propagation of nonlinear surface Alfven waves at a magnetic interface in a compressible atmosphere

Energy Technology Data Exchange (ETDEWEB)

Ruderman, M S

1988-08-01

Nonlinear Alfven surface wave propagation at a magnetic interface in a compressible fluid is considered. It is supposed that the magnetic field directions at both sides of the interface and the direction of wave propagation coincide. The equation governing time-evolution of nonlinear small-amplitude waves is derived by the method of multiscale expansions. This equation is similar to the equation for nonlinear Alfven surface waves in an incompressible fluid derived previously. The numerical solution of the equation shows that a sinusoidal disturbance overturns, i.e. infinite gradients arise.

Analytical and Numerical Modeling of Tsunami Wave Propagation for double layer state in Bore

Science.gov (United States)

Yuvaraj, V.; Rajasekaran, S.; Nagarajan, D.

2018-04-01

Tsunami wave enters into the river bore in the landslide. Tsunami wave propagation are described in two-layer states. The velocity and amplitude of the tsunami wave propagation are calculated using the double layer. The numerical and analytical solutions are given for the nonlinear equation of motion of the wave propagation in a bore.

True amplitude wave equation migration arising from true amplitude one-way wave equations

Science.gov (United States)

Zhang, Yu; Zhang, Guanquan; Bleistein, Norman

2003-10-01

One-way wave operators are powerful tools for use in forward modelling and inversion. Their implementation, however, involves introduction of the square root of an operator as a pseudo-differential operator. Furthermore, a simple factoring of the wave operator produces one-way wave equations that yield the same travel times as the full wave equation, but do not yield accurate amplitudes except for homogeneous media and for almost all points in heterogeneous media. Here, we present augmented one-way wave equations. We show that these equations yield solutions for which the leading order asymptotic amplitude as well as the travel time satisfy the same differential equations as the corresponding functions for the full wave equation. Exact representations of the square-root operator appearing in these differential equations are elusive, except in cases in which the heterogeneity of the medium is independent of the transverse spatial variables. Here, we address the fully heterogeneous case. Singling out depth as the preferred direction of propagation, we introduce a representation of the square-root operator as an integral in which a rational function of the transverse Laplacian appears in the integrand. This allows us to carry out explicit asymptotic analysis of the resulting one-way wave equations. To do this, we introduce an auxiliary function that satisfies a lower dimensional wave equation in transverse spatial variables only. We prove that ray theory for these one-way wave equations leads to one-way eikonal equations and the correct leading order transport equation for the full wave equation. We then introduce appropriate boundary conditions at z = 0 to generate waves at depth whose quotient leads to a reflector map and an estimate of the ray theoretical reflection coefficient on the reflector. Thus, these true amplitude one-way wave equations lead to a 'true amplitude wave equation migration' (WEM) method. In fact, we prove that applying the WEM imaging condition

Thermoelastic wave propagation in laminated composites plates

Directory of Open Access Journals (Sweden)

Verma K. L.

2012-12-01

Full Text Available The dispersion of thermoelastic waves propagation in an arbitrary direction in laminated composites plates is studied in the framework of generalized thermoelasticity in this article. Three dimensional field equations of thermoelasticity with relaxation times are considered. Characteristic equation is obtained on employing the continuity of displacements, temperature, stresses and thermal gradient at the layers’ interfaces. Some important particular cases such as of free waves on reducing plates to single layer and the surface waves when thickness tends to infinity are also discussed. Uncoupled and coupled thermoelasticity are the particular cases of the obtained results. Numerical results are also obtained and represented graphically.

Propagation of three-dimensional electron-acoustic solitary waves

International Nuclear Information System (INIS)

Shalaby, M.; El-Sherif, L. S.; El-Labany, S. K.; Sabry, R.

2011-01-01

Theoretical investigation is carried out for understanding the properties of three-dimensional electron-acoustic waves propagating in magnetized plasma whose constituents are cold magnetized electron fluid, hot electrons obeying nonthermal distribution, and stationary ions. For this purpose, the hydrodynamic equations for the cold magnetized electron fluid, nonthermal electron density distribution, and the Poisson equation are used to derive the corresponding nonlinear evolution equation, Zkharov-Kuznetsov (ZK) equation, in the small- but finite- amplitude regime. The ZK equation is solved analytically and it is found that it supports both solitary and blow-up solutions. It is found that rarefactive electron-acoustic solitary waves strongly depend on the density and temperature ratios of the hot-to-cold electron species as well as the nonthermal electron parameter. Furthermore, there is a critical value for the nonthermal electron parameter, which decides whether the electron-acoustic solitary wave's amplitude is decreased or increased by changing various plasma parameters. Importantly, the change of the propagation angles leads to miss the balance between the nonlinearity and dispersion; hence, the localized pulses convert to explosive/blow-up pulses. The relevance of this study to the nonlinear electron-acoustic structures in the dayside auroral zone in the light of Viking satellite observations is discussed.

Invertible propagator for plane wave illumination of forward-scattering structures.

Science.gov (United States)

Samelsohn, Gregory

2017-05-10

Propagation of directed waves in forward-scattering media is considered. It is assumed that the evolution of the wave field is governed by the standard parabolic wave equation. An efficient one-step momentum-space propagator, suitable for a tilted plane wave illumination of extended objects, is derived. It is expressed in terms of a propagation operator that transforms (the complex exponential of) a linogram of the illuminated object into a set of its diffraction patterns. The invertibility of the propagator is demonstrated, which permits a multiple-shot scatter correction to be performed, and makes the solution especially attractive for either projective or tomographic imaging. As an example, high-resolution tomograms are obtained in numerical simulations implemented for a synthetic phantom, with both refractive and absorptive inclusions.

On propagation of electromagnetic and gravitational waves in the expanding Universe

International Nuclear Information System (INIS)

Gladyshev, V O

2016-01-01

The purpose of this study was to obtain an equation for the propagation time of electromagnetic and gravitational waves in the expanding Universe. The velocity of electromagnetic waves propagation depends on the velocity of the interstellar medium in the observer's frame of reference. Gravitational radiation interacts weakly with the substance, so electromagnetic and gravitational waves propagate from a remote astrophysical object to the terrestrial observer at different time. Gravitational waves registration enables the inverse problem solution - by the difference in arrival time of electromagnetic and gravitational-wave signal, we can determine the characteristics of the emitting area of the astrophysical object. (paper)

Propagation of sech2-type solitary waves in higher-order KdV-type systems

International Nuclear Information System (INIS)

Ilison, O.; Salupere, A.

2005-01-01

Wave propagation in microstructured media is essentially influenced by nonlinear and dispersive effects. The simplest model governing these effects results in the Korteweg-de Vries (KdV) equation. In the present paper a KdV-type evolution equation, including the third- and fifth-order dispersive and the fourth-order nonlinear terms, is used for modelling the wave propagation in microstructured solids like martensitic-austenitic alloys. The model equation is solved numerically under localised initial conditions. Possible solution types are defined and discussed. The existence of a threshold is established. Below the threshold, the relatively small solitary waves decay in time. However, if the amplitude exceeds a certain threshold, i.e., the critical value, then such a solitary wave can propagate with nearly a constant speed and amplitude and consequently conserve the energy

Deterministic and stochastic evolution equations for fully dispersive and weakly nonlinear waves

DEFF Research Database (Denmark)

Eldeberky, Y.; Madsen, Per A.

1999-01-01

and stochastic formulations are solved numerically for the case of cross shore motion of unidirectional waves and the results are verified against laboratory data for wave propagation over submerged bars and over a plane slope. Outside the surf zone the two model predictions are generally in good agreement......This paper presents a new and more accurate set of deterministic evolution equations for the propagation of fully dispersive, weakly nonlinear, irregular, multidirectional waves. The equations are derived directly from the Laplace equation with leading order nonlinearity in the surface boundary...... is significantly underestimated for larger wave numbers. In the present work we correct this inconsistency. In addition to the improved deterministic formulation, we present improved stochastic evolution equations in terms of the energy spectrum and the bispectrum for multidirectional waves. The deterministic...

The effect of convection and shear on the damping and propagation of pressure waves

Science.gov (United States)

Kiel, Barry Vincent

Combustion instability is the positive feedback between heat release and pressure in a combustion system. Combustion instability occurs in the both air breathing and rocket propulsion devices, frequently resulting in high amplitude spinning waves. If unchecked, the resultant pressure fluctuations can cause significant damage. Models for the prediction of combustion instability typically include models for the heat release, the wave propagation and damping. Many wave propagation models for propulsion systems assume negligible flow, resulting in the wave equation. In this research the effect of flow on wave propagation was studied both numerically and experimentally. Two experiential rigs were constructed, one with axial flow to study the longitudinal waves, the other with swirling flow to study circumferential waves. The rigs were excited with speakers and the resultant pressure was measured simultaneously at many locations. Models of the rig were also developed. Equations for wave propagation were derived from the Euler Equations. The resultant resembled the wave equation with three additional terms, two for the effect of the convection and a one for the effect of shear of the mean flow on wave propagation. From the experimental and numerical data several conclusions were made. First, convection and shear both act as damping on the wave propagation, reducing the magnitude of the Frequency Response Function and the resonant frequency of the modes. Second, the energy extracted from the mean flow as a result of turbulent shear for a given condition is frequency dependent, decreasing with increasing frequency. The damping of the modes, measured for the same shear flow, also decreased with frequency. Finally, the two convective terms cause the anti-nodes of the modes to no longer be stationary. For both the longitudinal and circumferential waves, the anti-nodes move through the domain even for mean flow Mach numbers less than 0.10. It was concluded that convection

A Potential Method for Body and Surface Wave Propagation in Transversely Isotropic Half- and Full-Spaces

Directory of Open Access Journals (Sweden)

Mehdi Raoofian Naeeni

2016-12-01

Full Text Available The problem of propagation of plane wave including body and surface waves propagating in a transversely isotropic half-space with a depth-wise axis of material symmetry is investigated in details. Using the advantage of representation of displacement fields in terms of two complete scalar potential functions, the coupled equations of motion are uncoupled and reduced to two independent equations for potential functions. In this paper, the secular equations for determination of body and surface wave velocities are derived in terms of both elasticity coefficients and the direction of propagation. In particular, the longitudinal, transverse and Rayleigh wave velocities are determined in explicit forms. It is also shown that in transversely isotropic materials, a Rayleigh wave may propagate in different manner from that of isotropic materials. Some numerical results for synthetic transversely isotropic materials are also illustrated to show the behavior of wave motion due to anisotropic nature of the problem.

Breatherlike electromagnetic wave propagation in an antiferromagnetic medium with Dzyaloshinsky-Moriya interaction

International Nuclear Information System (INIS)

Kavitha, L.; Saravanan, M.; Srividya, B.; Gopi, D.

2011-01-01

We investigate the nature of propagation of electromagnetic waves (EMWs) in an antiferromagnetic medium with Dzyaloshinsky-Moriya (DM) interaction environment. The interplay of bilinear and DM exchange spin coupling with the magnetic field component of the EMW has been studied by solving Maxwell's equations coupled with a nonlinear spin equation for the magnetization of the medium. We made a nonuniform expansion of the magnetization and magnetic field along the direction of propagation of EMW, in the framework of reductive perturbation method, and the dynamics of the system is found to be governed by a generalized derivative nonlinear Schroedinger (DNLS) equation. We employ the Jacobi-elliptic function method to solve the DNLS equation, and the electromagnetic wave propagation in an antiferromagnetic medium is governed by the breatherlike spatially and temporally coherent localized modes under the influence of DM interaction parameter.

Wave propagation in a quasi-chemical equilibrium plasma

Science.gov (United States)

Fang, T.-M.; Baum, H. R.

1975-01-01

Wave propagation in a quasi-chemical equilibrium plasma is studied. The plasma is infinite and without external fields. The chemical reactions are assumed to result from the ionization and recombination processes. When the gas is near equilibrium, the dominant role describing the evolution of a reacting plasma is played by the global conservation equations. These equations are first derived and then used to study the small amplitude wave motion for a near-equilibrium situation. Nontrivial damping effects have been obtained by including the conduction current terms.

Nonlinear evolution equations for waves in random media

International Nuclear Information System (INIS)

Pelinovsky, E.; Talipova, T.

1994-01-01

The scope of this paper is to highlight the main ideas of asymptotical methods applying in modern approaches of description of nonlinear wave propagation in random media. We start with the discussion of the classical conception of ''mean field''. Then an exactly solvable model describing nonlinear wave propagation in the medium with fluctuating parameters is considered in order to demonstrate that the ''mean field'' method is not correct. We develop new asymptotic procedures of obtaining the nonlinear evolution equations for the wave fields in random media. (author). 16 refs

Nonlinear sausage-wave propagation in a magnetic slab in an incompressible fluid

International Nuclear Information System (INIS)

Ruderman, M.S.

1993-01-01

Long nonlinear sausage-wave propagation in a magnetic slab in an incompressible plasma is considered. The governing equation is derived with the aid of the reductive perturbation method. The solutions of this equation in the form of periodic waves of permanent shape are found numerically. (Author)

Earthquake wave propagation in immiscibly compressible porous soil

International Nuclear Information System (INIS)

Xue, S.; Kurita, S.; Izumi, M.

1993-01-01

This paper utilizes the formalism of the theory of immiscible compressible mixtures to formulate the wave propagation equation for the soil where the soil has been assumed as a binary mixture consisting of one solid phase and one fluid phase. The method is developed to solve the one dimensional wave equation by the above theory. The relations between the wave attenuating characteristic value Q and the volume fraction, the relative motion of two phases have been shown. It is concluded that based on such theory we can solve more precisely the soil behaviors while considering the interaction of structure and soil of immiscible mixture. (author)

The Appell transformation for the paraxial wave equation

International Nuclear Information System (INIS)

Torre, A

2011-01-01

Some issues related to the 1D heat equation are revisited and framed within the context of the free-space paraxial propagation, formally accounted for by the 2D paraxial wave equation. In particular, the Appell transformation, which is well known in the theory of the heat equation, is reformulated in optical terms, and accordingly interpreted in the light of the propagation of given source functions, which are in a definite relation with the source functions of the original wavefunctions. Basic to the discussion is the Lie-algebra-based approach, as developed in a series of seminal papers by Kalnins, Miller and Boyer, to evolutionary-type equations, ruled by Hamiltonian operators underlying a harmonic oscillator-like symmetry algebra. Indeed, both the heat equation and the paraxial wave equation are particular cases of this kind of equation. When interpreting such an approach in terms of the propagation of assigned 'source' functions, the transformations between wavefunctions may be traced back to definite relations between the respective source functions. Thus, the optical Appell transformation is seen to be a manifestation of the correspondence between wavefunctions generated by eigenstates of operators, which are linked through a Fourier-similarity transformation. As a mere consequence, one can introduce the fractional Appell transformation, thus displaying a family of symmetry transformations parameterized by a continuous parameter

Stability of negative solitary waves for an integrable modified Camassa-Holm equation

International Nuclear Information System (INIS)

Yin Jiuli; Tian Lixin; Fan Xinghua

2010-01-01

In this paper, we prove that the modified Camassa-Holm equation is Painleve integrable. We also study the orbital stability problem of negative solitary waves for this integrable equation. It is shown that the negative solitary waves are stable for arbitrary wave speed of propagation.

Wave-equation Q tomography

KAUST Repository

Dutta, Gaurav

2016-10-12

Strong subsurface attenuation leads to distortion of amplitudes and phases of seismic waves propagating inside the earth. The amplitude and the dispersion losses from attenuation are often compensated for during prestack depth migration. However, most attenuation compensation or Qcompensation migration algorithms require an estimate of the background Q model. We have developed a wave-equation gradient optimization method that inverts for the subsurface Q distribution by minimizing a skeletonized misfit function ∈, where ∈ is the sum of the squared differences between the observed and the predicted peak/centroid-frequency shifts of the early arrivals. The gradient is computed by migrating the observed traces weighted by the frequency shift residuals. The background Q model is perturbed until the predicted and the observed traces have the same peak frequencies or the same centroid frequencies. Numerical tests determined that an improved accuracy of the Q model by wave-equation Q tomography leads to a noticeable improvement in migration image quality. © 2016 Society of Exploration Geophysicists.

Wave-equation Q tomography

KAUST Repository

Dutta, Gaurav; Schuster, Gerard T.

2016-01-01

Strong subsurface attenuation leads to distortion of amplitudes and phases of seismic waves propagating inside the earth. The amplitude and the dispersion losses from attenuation are often compensated for during prestack depth migration. However, most attenuation compensation or Qcompensation migration algorithms require an estimate of the background Q model. We have developed a wave-equation gradient optimization method that inverts for the subsurface Q distribution by minimizing a skeletonized misfit function ∈, where ∈ is the sum of the squared differences between the observed and the predicted peak/centroid-frequency shifts of the early arrivals. The gradient is computed by migrating the observed traces weighted by the frequency shift residuals. The background Q model is perturbed until the predicted and the observed traces have the same peak frequencies or the same centroid frequencies. Numerical tests determined that an improved accuracy of the Q model by wave-equation Q tomography leads to a noticeable improvement in migration image quality. © 2016 Society of Exploration Geophysicists.

Pressure wave propagation in sodium loop

International Nuclear Information System (INIS)

Botelho, D.A.

1989-01-01

A study was done on the pressure wave propagation within the pipes and mixture vessel of a termohydraulic loop for thermal shock with sodium. It was used the characteristic method to solve the one-dimensional continuity and momentum equations. The numerical model includes the pipes and the effects of valves and other accidents on pressure losses. The study was based on designer informations and engineering tables. It was evaluated the pressure wave sizes, parametrically as a function of the draining valve closure times. (author) [pt

Nonlocal symmetries, solitary waves and cnoidal periodic waves of the (2+1)-dimensional breaking soliton equation

Science.gov (United States)

Zou, Li; Tian, Shou-Fu; Feng, Lian-Li

2017-12-01

In this paper, we consider the (2+1)-dimensional breaking soliton equation, which describes the interaction of a Riemann wave propagating along the y-axis with a long wave along the x-axis. By virtue of the truncated Painlevé expansion method, we obtain the nonlocal symmetry, Bäcklund transformation and Schwarzian form of the equation. Furthermore, by using the consistent Riccati expansion (CRE), we prove that the breaking soliton equation is solvable. Based on the consistent tan-function expansion, we explicitly derive the interaction solutions between solitary waves and cnoidal periodic waves.

Numerical simulation methods for wave propagation through optical waveguides

International Nuclear Information System (INIS)

Sharma, A.

1993-01-01

The simulation of the field propagation through waveguides requires numerical solutions of the Helmholtz equation. For this purpose a method based on the principle of orthogonal collocation was recently developed. The method is also applicable to nonlinear pulse propagation through optical fibers. Some of the salient features of this method and its application to both linear and nonlinear wave propagation through optical waveguides are discussed in this report. 51 refs, 8 figs, 2 tabs

WAVE: Interactive Wave-based Sound Propagation for Virtual Environments.

Science.gov (United States)

Mehra, Ravish; Rungta, Atul; Golas, Abhinav; Ming Lin; Manocha, Dinesh

2015-04-01

We present an interactive wave-based sound propagation system that generates accurate, realistic sound in virtual environments for dynamic (moving) sources and listeners. We propose a novel algorithm to accurately solve the wave equation for dynamic sources and listeners using a combination of precomputation techniques and GPU-based runtime evaluation. Our system can handle large environments typically used in VR applications, compute spatial sound corresponding to listener's motion (including head tracking) and handle both omnidirectional and directional sources, all at interactive rates. As compared to prior wave-based techniques applied to large scenes with moving sources, we observe significant improvement in runtime memory. The overall sound-propagation and rendering system has been integrated with the Half-Life 2 game engine, Oculus-Rift head-mounted display, and the Xbox game controller to enable users to experience high-quality acoustic effects (e.g., amplification, diffraction low-passing, high-order scattering) and spatial audio, based on their interactions in the VR application. We provide the results of preliminary user evaluations, conducted to study the impact of wave-based acoustic effects and spatial audio on users' navigation performance in virtual environments.

Counterstreaming magnetized plasmas. II. Perpendicular wave propagation

International Nuclear Information System (INIS)

Tautz, R.C.; Schlickeiser, R.

2006-01-01

The properties of longitudinal and transverse oscillations in magnetized symmetric counterstreaming Maxwellian plasmas with equal thermal velocities for waves propagating perpendicular to the stream direction are investigated on the basis of Maxwell equations and the nonrelativistic Vlasov equation. With the constraint of vanishing particle flux in the stream direction, three distinct dispersion relations are known, which are the ordinary-wave mode, the Bernstein wave mode, and the extraordinary electromagnetic wave mode, where the latter two are only approximations. In this article, all three dispersion relations are evaluated for a counterstreaming Maxwellian distribution function in terms of the hypergeometric function 2 F 2 . The growth rates for the ordinary-wave mode are compared to earlier results by Bornatici and Lee [Phys. Fluids 13, 3007 (1970)], who derived approximate results, whereas in this article the exact dispersion relation is solved numerically. The original results are therefore improved and show differences of up to 21% to the results obtained in this article

A phase space approach to wave propagation with dispersion.

Science.gov (United States)

Ben-Benjamin, Jonathan S; Cohen, Leon; Loughlin, Patrick J

2015-08-01

A phase space approximation method for linear dispersive wave propagation with arbitrary initial conditions is developed. The results expand on a previous approximation in terms of the Wigner distribution of a single mode. In contrast to this previously considered single-mode case, the approximation presented here is for the full wave and is obtained by a different approach. This solution requires one to obtain (i) the initial modal functions from the given initial wave, and (ii) the initial cross-Wigner distribution between different modal functions. The full wave is the sum of modal functions. The approximation is obtained for general linear wave equations by transforming the equations to phase space, and then solving in the new domain. It is shown that each modal function of the wave satisfies a Schrödinger-type equation where the equivalent "Hamiltonian" operator is the dispersion relation corresponding to the mode and where the wavenumber is replaced by the wavenumber operator. Application to the beam equation is considered to illustrate the approach.

Stress Wave Propagation in Viscoelastic-Plastic Rock-Like Materials

Directory of Open Access Journals (Sweden)

Liu Lang

2016-05-01

Full Text Available Rock-like materials are composites that can be regarded as a mixture composed of elastic, plastic, and viscous components. They exhibit viscoelastic-plastic behavior under a high-strain-rate loading according to element model theory. This paper presents an analytical solution for stress wave propagation in viscoelastic-plastic rock-like materials under a high-strain-rate loading and verifies the solution through an experimental test. A constitutive equation of viscoelastic-plastic rock-like materials was first established, and then kinematic and kinetic equations were then solved to derive the analytic solution for stress wave propagation in viscoelastic-plastic rock-like materials. An experimental test using the SHPB (Split Hopkinson Pressure Bar for a concrete specimen was conducted to obtain a stress-strain curve under a high-strain-rate loading. Inverse analysis based on differential evolution was conducted to estimate undetermined variables for constitutive equations. Finally, the relationship between the attenuation factor and the strain rate in viscoelastic-plastic rock-like materials was investigated. According to the results, the frequency of the stress wave, viscosity coefficient, modulus of elasticity, and density play dominant roles in the attenuation of the stress wave. The attenuation decreases with increasing strain rate, demonstrating strongly strain-dependent attenuation in viscoelastic-plastic rock-like materials.

On an Acoustic Wave Equation Arising in Non-Equilibrium Gasdynamics. Classroom Notes

Science.gov (United States)

Chandran, Pallath

2004-01-01

The sixth-order wave equation governing the propagation of one-dimensional acoustic waves in a viscous, heat conducting gaseous medium subject to relaxation effects has been considered. It has been reduced to a system of lower order equations corresponding to the finite speeds occurring in the equation, following a method due to Whitham. The lower…

Wave propagation in embedded inhomogeneous nanoscale plates incorporating thermal effects

Science.gov (United States)

Ebrahimi, Farzad; Barati, Mohammad Reza; Dabbagh, Ali

2018-04-01

In this article, an analytical approach is developed to study the effects of thermal loading on the wave propagation characteristics of an embedded functionally graded (FG) nanoplate based on refined four-variable plate theory. The heat conduction equation is solved to derive the nonlinear temperature distribution across the thickness. Temperature-dependent material properties of nanoplate are graded using Mori-Tanaka model. The nonlocal elasticity theory of Eringen is introduced to consider small-scale effects. The governing equations are derived by the means of Hamilton's principle. Obtained frequencies are validated with those of previously published works. Effects of different parameters such as temperature distribution, foundation parameters, nonlocal parameter, and gradient index on the wave propagation response of size-dependent FG nanoplates have been investigated.

Propagation of waves

CERN Document Server

David, P

2013-01-01

Propagation of Waves focuses on the wave propagation around the earth, which is influenced by its curvature, surface irregularities, and by passage through atmospheric layers that may be refracting, absorbing, or ionized. This book begins by outlining the behavior of waves in the various media and at their interfaces, which simplifies the basic phenomena, such as absorption, refraction, reflection, and interference. Applications to the case of the terrestrial sphere are also discussed as a natural generalization. Following the deliberation on the diffraction of the "ground? wave around the ear

Statistical Characterization of Electromagnetic Wave Propagation in Mine Environments

KAUST Repository

Yucel, Abdulkadir C.

2013-01-01

A computational framework for statistically characterizing electromagnetic (EM) wave propagation through mine tunnels and galleries is presented. The framework combines a multi-element probabilistic collocation method with a full-wave fast Fourier transform and fast multipole method accelerated surface integral equation-based EM simulator to statistically characterize fields from wireless transmitters in complex mine environments. 1536-1225 © 2013 IEEE.

Wave propagation model of heat conduction and group speed

Science.gov (United States)

Zhang, Long; Zhang, Xiaomin; Peng, Song

2018-03-01

In view of the finite relaxation model of non-Fourier's law, the Cattaneo and Vernotte (CV) model and Fourier's law are presented in this work for comparing wave propagation modes. Independent variable translation is applied to solve the partial differential equation. Results show that the general form of the time spatial distribution of temperature for the three media comprises two solutions: those corresponding to the positive and negative logarithmic heating rates. The former shows that a group of heat waves whose spatial distribution follows the exponential function law propagates at a group speed; the speed of propagation is related to the logarithmic heating rate. The total speed of all the possible heat waves can be combined to form the group speed of the wave propagation. The latter indicates that the spatial distribution of temperature, which follows the exponential function law, decays with time. These features show that propagation accelerates when heated and decelerates when cooled. For the model media that follow Fourier's law and correspond to the positive heat rate of heat conduction, the propagation mode is also considered the propagation of a group of heat waves because the group speed has no upper bound. For the finite relaxation model with non-Fourier media, the interval of group speed is bounded and the maximum speed can be obtained when the logarithmic heating rate is exactly the reciprocal of relaxation time. And for the CV model with a non-Fourier medium, the interval of group speed is also bounded and the maximum value can be obtained when the logarithmic heating rate is infinite.

Bulk elastic wave propagation in partially saturated porous solids

International Nuclear Information System (INIS)

Berryman, J.G.; Thigpen, L.; Chin, R.C.Y.

1988-01-01

The linear equations of motion that describe the behavior of small disturbances in a porous solid containing both liquid and gas are solved for bulk wave propagation. The equations have been simplified by neglecting effects due to changes in capillary pressure. With this simplifying assumption, the equations reduce to two coupled (vector) equations of the form found in Biot's equations (for full saturation) but with more complicated coefficients. As in fully saturated solids, two shear waves with the same speed but different polarizations exist as do two compressional waves with distinct speeds. Attenuation effects can be enhanced in the partially saturated solid, depending on the distribution of gas in the pore space. Two models of the liquid/gas spatial distribution are considered: a segregated-fluids model and a mixed-fluids model. The two models predict comparable attentuation when the gas saturation is low, but the segregated-fluids model predicts a more rapid roll-off of attenuation as the gas saturation increases

Nonlinear Wave Propagation and Solitary Wave Formation in Two-Dimensional Heterogeneous Media

KAUST Repository

Luna, Manuel

2011-05-01

Solitary wave formation is a well studied nonlinear phenomenon arising in propagation of dispersive nonlinear waves under suitable conditions. In non-homogeneous materials, dispersion may happen due to effective reflections between the material interfaces. This dispersion has been used along with nonlinearities to find solitary wave formation using the one-dimensional p-system. These solitary waves are called stegotons. The main goal in this work is to find two-dimensional stegoton formation. To do so we consider the nonlinear two-dimensional p-system with variable coefficients and solve it using finite volume methods. The second goal is to obtain effective equations that describe the macroscopic behavior of the variable coefficient system by a constant coefficient one. This is done through a homogenization process based on multiple-scale asymptotic expansions. We compare the solution of the effective equations with the finite volume results and find a good agreement. Finally, we study some stability properties of the homogenized equations and find they and one-dimensional versions of them are unstable in general.

Low-frequency dilatational wave propagation through unsaturated porous media containing two immiscible fluids

Energy Technology Data Exchange (ETDEWEB)

Lo, W.-C.; Sposito, G.; Majer, E.

2007-02-01

An analytical theory is presented for the low-frequency behavior of dilatational waves propagating through a homogeneous elastic porous medium containing two immiscible fluids. The theory is based on the Berryman-Thigpen-Chin (BTC) model, in which capillary pressure effects are neglected. We show that the BTC model equations in the frequency domain can be transformed, at sufficiently low frequencies, into a dissipative wave equation (telegraph equation) and a propagating wave equation in the time domain. These partial differential equations describe two independent modes of dilatational wave motion that are analogous to the Biot fast and slow compressional waves in a single-fluid system. The equations can be solved analytically under a variety of initial and boundary conditions. The stipulation of 'low frequency' underlying the derivation of our equations in the time domain is shown to require that the excitation frequency of wave motions be much smaller than a critical frequency. This frequency is shown to be the inverse of an intrinsic time scale that depends on an effective kinematic shear viscosity of the interstitial fluids and the intrinsic permeability of the porous medium. Numerical calculations indicate that the critical frequency in both unconsolidated and consolidated materials containing water and a nonaqueous phase liquid ranges typically from kHz to MHz. Thus engineering problems involving the dynamic response of an unsaturated porous medium to low excitation frequencies (e.g. seismic wave stimulation) should be accurately modeled by our equations after suitable initial and boundary conditions are imposed.

Nonlocal wave propagation in an embedded DWBNNT conveying fluid via strain gradient theory

International Nuclear Information System (INIS)

Ghorbanpour Arani, A.; Kolahchi, R.; Vossough, H.

2012-01-01

Based on the strain gradient and Eringen’s piezoelasticity theories, wave propagation of an embedded double-walled boron nitride nanotube (DWBNNT) conveying fluid is investigated using Euler-Bernoulli beam model. The elastic medium is simulated by the Pasternak foundation. The van der Waals (vdW) forces between the inner and outer nanotubes are taken into account. Since, considering electro-mechanical coupling made the nonlinear motion equations, a numerical procedure is proposed to evaluate the upstream and downstream phase velocities. The results indicate that the effect of nonlinear terms in motion equations on the phase velocity cannot be neglected at lower wave numbers. Furthermore, the effect of fluid-conveying on wave propagation of the DWBNNT is significant at lower wave numbers.

Transverse wave propagation in [ab0] direction of silicon single crystal

Energy Technology Data Exchange (ETDEWEB)

Yun, Sang Jin; Kim, Hye Jeong; Kwon, Se Ho; Kim, Young H. [Applied Acoustics Lab, Korea Science Academy of KAIST, Busan(Korea, Republic of)

2015-12-15

The speed and oscillation directions of elastic waves propagating in the [ab0] direction of a silicon single crystal were obtained by solving Christoffel's equation. It was found that the quasi waves propagate in the off-principal axis, and hence, the directions of the phase and group velocities are not the same. The maximum deviation of the two directions was 7.2 degree angle. Two modes of the pure transverse waves propagate in the [110] direction with different speeds, and hence, two peaks were observed in the pulse echo signal. The amplitude ratio of the two peaks was dependent on the initial oscillating direction of the incident wave. The pure and quasi-transverse waves propagate in the [210] direction, and the oscillation directions of these waves are perpendicular to each other. The skewing angle of the quasi wave was calculated as 7.14 degree angle, and it was measured as 9.76 degree angle. The amplitude decomposition in the [210] direction was similar to that in the [110] direction, since the oscillation directions of these waves are perpendicular to each other. These results offer useful information in measuring the crystal orientation of the silicon single crystal.

Reduction of the equation for lower hybrid waves in a plasma to a nonlinear Schroedinger equation

Science.gov (United States)

Karney, C. F. F.

1977-01-01

Equations describing the nonlinear propagation of waves in an anisotropic plasma are rarely exactly soluble. However it is often possible to make approximations that reduce the exact equations into a simpler equation. The use of MACSYMA to make such approximations, and so reduce the equation describing lower hybrid waves into the nonlinear Schrodinger equation which is soluble by the inverse scattering method is demonstrated. MACSYMA is used at several stages in the calculation only because there is a natural division between calculations that are easiest done by hand, and those that are easiest done by machine.

A frequency domain linearized Navier-Stokes equations approach to acoustic propagation in flow ducts with sharp edges.

Science.gov (United States)

Kierkegaard, Axel; Boij, Susann; Efraimsson, Gunilla

2010-02-01

Acoustic wave propagation in flow ducts is commonly modeled with time-domain non-linear Navier-Stokes equation methodologies. To reduce computational effort, investigations of a linearized approach in frequency domain are carried out. Calculations of sound wave propagation in a straight duct are presented with an orifice plate and a mean flow present. Results of transmission and reflections at the orifice are presented on a two-port scattering matrix form and are compared to measurements with good agreement. The wave propagation is modeled with a frequency domain linearized Navier-Stokes equation methodology. This methodology is found to be efficient for cases where the acoustic field does not alter the mean flow field, i.e., when whistling does not occur.

Wave Propagation in Bimodular Geomaterials

Science.gov (United States)

Kuznetsova, Maria; Pasternak, Elena; Dyskin, Arcady; Pelinovsky, Efim

2016-04-01

Observations and laboratory experiments show that fragmented or layered geomaterials have the mechanical response dependent on the sign of the load. The most adequate model accounting for this effect is the theory of bimodular (bilinear) elasticity - a hyperelastic model with different elastic moduli for tension and compression. For most of geo- and structural materials (cohesionless soils, rocks, concrete, etc.) the difference between elastic moduli is such that their modulus in compression is considerably higher than that in tension. This feature has a profound effect on oscillations [1]; however, its effect on wave propagation has not been comprehensively investigated. It is believed that incorporation of bilinear elastic constitutive equations within theory of wave dynamics will bring a deeper insight to the study of mechanical behaviour of many geomaterials. The aim of this paper is to construct a mathematical model and develop analytical methods and numerical algorithms for analysing wave propagation in bimodular materials. Geophysical and exploration applications and applications in structural engineering are envisaged. The FEM modelling of wave propagation in a 1D semi-infinite bimodular material has been performed with the use of Marlow potential [2]. In the case of the initial load expressed by a harmonic pulse loading strong dependence on the pulse sign is observed: when tension is applied before compression, the phenomenon of disappearance of negative (compressive) strains takes place. References 1. Dyskin, A., Pasternak, E., & Pelinovsky, E. (2012). Periodic motions and resonances of impact oscillators. Journal of Sound and Vibration, 331(12), 2856-2873. 2. Marlow, R. S. (2008). A Second-Invariant Extension of the Marlow Model: Representing Tension and Compression Data Exactly. In ABAQUS Users' Conference.

Nonlinear magnetoacoustic wave propagation with chemical reactions

Science.gov (United States)

Margulies, Timothy Scott

2002-11-01

The magnetoacoustic problem with an application to sound wave propagation through electrically conducting fluids such as the ocean in the Earth's magnetic field, liquid metals, or plasmas has been addressed taking into account several simultaneous chemical reactions. Using continuum balance equations for the total mass, linear momentum, energy; as well as Maxwell's electrodynamic equations, a nonlinear beam equation has been developed to generalize the Khokhlov-Zabolotskaya-Kuznetsov (KZK) equation for a fluid with linear viscosity but nonlinear and diffraction effects. Thermodynamic parameters are used and not tailored to only an adiabatic fluid case. The chemical kinetic equations build on a relaxing media approach presented, for example, by K. Naugolnukh and L. Ostrovsky [Nonlinear Wave Processes in Acoustics (Cambridge Univ. Press, Cambridge, 1998)] for a linearized single reaction and thermodynamic pressure equation of state. Approximations for large and small relaxation times and for magnetohydrodynamic parameters [Korsunskii, Sov. Phys. Acoust. 36 (1990)] are examined. Additionally, Cattaneo's equation for heat conduction and its generalization for a memory process rather than a Fourier's law are taken into account. It was introduced for the heat flux depends on the temperature gradient at an earlier time to generate heat pulses of finite speed.

Wave propagation in elastic solids

CERN Document Server

Achenbach, Jan

1984-01-01

The propagation of mechanical disturbances in solids is of interest in many branches of the physical scienses and engineering. This book aims to present an account of the theory of wave propagation in elastic solids. The material is arranged to present an exposition of the basic concepts of mechanical wave propagation within a one-dimensional setting and a discussion of formal aspects of elastodynamic theory in three dimensions, followed by chapters expounding on typical wave propagation phenomena, such as radiation, reflection, refraction, propagation in waveguides, and diffraction. The treat

Wave propagation in a piezoelectric solid bar of circular cross-section immersed in fluid

International Nuclear Information System (INIS)

Ponnusamy, P.

2013-01-01

Wave propagation in a piezoelectric solid bar of circular cross-section immersed in fluid is discussed using three-dimensional theory of piezoelectricity. The equations of motion of the cylinder are formulated using the constitutive equations of a piezoelectric material. The equations of motion of the fluid are formulated using the constitutive equations of an inviscid fluid. Three displacement potential functions are introduced to uncouple the equations of motion, electric conduction. The frequency equation of the coupled system consisting of cylinder and fluid is developed under the assumption of perfect-slip boundary conditions at the fluid–solid interfaces. The frequency equations are obtained for longitudinal and flexural modes of vibration and are studied numerically for PZT-4 material bar immersed in fluid. The computed non-dimensional wave numbers are presented in the form of dispersion curves. The secant method is used to obtain the roots of the frequency equation. -- Highlights: ► Wave propagation in a piezoelectric solid bar of circular cross-section immersed in fluid is analyzed using secant method. ► Solid–fluid interaction for piezoelectric material of PZT-4 is analyzed using the boundary conditions. ► The computed non-dimensional wave numbers are plotted in the form of dispersion curves and studied its characters. ► A comparison is made between the non-dimensional wave numbers obtained by the author with the literature results

Nonlinear elastic longitudinal strain-wave propagation in a plate with nonequilibrium laser-generated point defects

International Nuclear Information System (INIS)

Mirzade, Fikret Kh.

2005-01-01

The propagation of longitudinal strain wave in a plate with quadratic nonlinearity of elastic continuum was studied in the context of a model that takes into account the joint dynamics of elastic displacements in the medium and the concentration of the nonequilibrium laser-induced point defects. The input equations of the problem are reformulated in terms of only the total displacements of the medium points. In this case, the presence of structural defects manifests itself in the emergence of a delayed response of the system to the propagation of the strain-related perturbations, which is characteristic of media with relaxation or memory. The model equations describing the nonlinear displacement wave were derived with allowance made for the values of the relaxation parameter. The influence of the generation and relaxation of lattice defects on the propagation of this wave was analyzed. It is shown that, for short relaxation times of defects, the strain can propagate in the form of shock fronts. In the case of longer relaxation times, shock waves do not form and the strain wave propagates only in the form of solitary waves or a train of solitons. The contributions of the finiteness of the defect-recombination rate to linear and nonlinear elastic modulus, and spatial dispersion are determined

Preserving the Helmholtz dispersion relation: One-way acoustic wave propagation using matrix square roots

Science.gov (United States)

Keefe, Laurence

2016-11-01

Parabolized acoustic propagation in transversely inhomogeneous media is described by the operator update equation U (x , y , z + Δz) =eik0 (- 1 +√{ 1 + Z }) U (x , y , z) for evolution of the envelope of a wavetrain solution to the original Helmholtz equation. Here the operator, Z =∇T2 + (n2 - 1) , involves the transverse Laplacian and the refractive index distribution. Standard expansion techniques (on the assumption Z << 1)) produce pdes that approximate, to greater or lesser extent, the full dispersion relation of the original Helmholtz equation, except that none of them describe evanescent/damped waves without special modifications to the expansion coefficients. Alternatively, a discretization of both the envelope and the operator converts the operator update equation into a matrix multiply, and existing theorems on matrix functions demonstrate that the complete (discrete) Helmholtz dispersion relation, including evanescent/damped waves, is preserved by this discretization. Propagation-constant/damping-rates contour comparisons for the operator equation and various approximations demonstrate this point, and how poorly the lowest-order, textbook, parabolized equation describes propagation in lined ducts.

Propagation of nonlinear waves over submerged step: wave separation and subharmonic generation

Science.gov (United States)

Monsalve, Eduardo; Maurel, Agnes; Pagneux, Vincent; Petitjeans, Philippe

2015-11-01

Water waves can be described in simplified cases by the Helmholtz equation. However, even in these cases, they present a high complexity, among which their dispersive character and their nonlinearities are the subject of the present study. Using Fourier Transform Profilometry, we study experimentally the propagation of waves passing over a submerged step. Because of the small water depth after the step, the wave enters in a nonlinear regime. In the shallow water region, the second harmonic leads to two types of waves: bound waves which are slaves of the fundamental frequency with wavenumber 2 k (ω) , and free waves which propagate according to the usual dispersion relation with wavenumber k (2 ω) . Because of the presence of these two waves, beats are produced at the second harmonic with characteristic beat length. In this work, for the first time we extended this analysis to the third and higher harmonics. Next, the region after the step is limited to a finite size L with a reflecting wall. For certain frequencies and L- values, the spectral component becomes involved, with the appearance of sub harmonics. This regime is analyzed in more details, suggesting a transition to a chaotic and quasi-periodic wave behavior.

Mathematical problems in wave propagation theory

CERN Document Server

1970-01-01

The papers comprising this collection are directly or indirectly related to an important branch of mathematical physics - the mathematical theory of wave propagation and diffraction. The paper by V. M. Babich is concerned with the application of the parabolic-equation method (of Academician V. A. Fok and M. A, Leontovich) to the problem of the asymptotic behavior of eigenfunc­ tions concentrated in a neighborhood of a closed geodesie in a Riemannian space. The techniques used in this paper have been föund useful in solving certain problems in the theory of open resonators. The topic of G. P. Astrakhantsev's paper is similar to that of the paper by V. M. Babich. Here also the parabolic-equation method is used to find the asymptotic solution of the elasticity equations which describes Love waves concentrated in a neighborhood of some surface ray. The paper of T. F. Pankratova is concerned with finding the asymptotic behavior of th~ eigenfunc­ tions of the Laplace operator from the exact solution for the surf...

Computational study of nonlinear plasma waves. I. Simulation model and monochromatic wave propagation

International Nuclear Information System (INIS)

Matsuda, Y.; Crawford, F.W.

1975-01-01

An economical low-noise plasma simulation model originated by Denavit is applied to a series of problems associated with electrostatic wave propagation in a one-dimensional, collisionless, Maxwellian plasma, in the absence of magnetic field. The model is described and tested, first in the absence of an applied signal, and then with a small amplitude perturbation. These tests serve to establish the low-noise features of the model, and to verify the theoretical linear dispersion relation at wave energy levels as low as 10 -6 of the plasma thermal energy: Better quantitative results are obtained, for comparable computing time, than can be obtained by conventional particle simulation models, or direct solution of the Vlasov equation. The method is then used to study propagation of an essentially monochromatic plane wave. Results on amplitude oscillation and nonlinear frequency shift are compared with available theories

Analysis of Electromagnetic Wave Propagation in a Magnetized Re-Entry Plasma Sheath Via the Kinetic Equation

Science.gov (United States)

Manning, Robert M.

2009-01-01

Based on a theoretical model of the propagation of electromagnetic waves through a hypersonically induced plasma, it has been demonstrated that the classical radiofrequency communications blackout that is experienced during atmospheric reentry can be mitigated through the appropriate control of an external magnetic field of nominal magnitude. The model is based on the kinetic equation treatment of Vlasov and involves an analytical solution for the electric and magnetic fields within the plasma allowing for a description of the attendant transmission, reflection and absorption coefficients. The ability to transmit through the magnetized plasma is due to the magnetic windows that are created within the plasma via the well-known whistler modes of propagation. The case of 2 GHz transmission through a re-entry plasma is considered. The coefficients are found to be highly sensitive to the prevailing electron density and will thus require a dynamic control mechanism to vary the magnetic field as the plasma evolves through the re-entry phase.

Propagation of an ionizing surface electromagnetic wave

Energy Technology Data Exchange (ETDEWEB)

Boev, A.G.; Prokopov, A.V.

1976-11-01

The propagation of an rf surface wave in a plasma which is ionized by the wave itself is analyzed. The exact solution of the nonlinear Maxwell equations is discussed for the case in which the density of plasma electrons is an exponential function of the square of the electric field. The range over which the surface wave exists and the frequency dependence of the phase velocity are found. A detailed analysis is given for the case of a plasma whose initial density exceeds the critical density at the wave frequency. An increase in the wave amplitude is shown to expand the frequency range over which the plasma is transparent; The energy flux in the plasma tends toward a certain finite value which is governed by the effective ionization field.

Effect of small floating disks on the propagation of gravity waves

Energy Technology Data Exchange (ETDEWEB)

Santi, F De; Olla, P, E-mail: olla@dsf.unica.it [ISAC-CNR, Sez. Cagliari, I-09042 Monserrato (Italy)

2017-04-15

A dispersion relation for gravity waves in water covered by disk-like impurities embedded in a viscous matrix is derived. The macroscopic equations are obtained by ensemble-averaging the fluid equations at the disk scale in the asymptotic limit of long waves and low disk surface fraction. Various regimes are identified depending on the disk radii and the thickness and viscosity of the top layer. Semi-quantitative analysis in the close-packing regime suggests dramatic modification of the dynamics, with orders of magnitude increase in wave damping and wave dispersion. A simplified model working in this regime is proposed. Possible applications to wave propagation in an ice-covered ocean are discussed and comparison with field data is provided. (paper)

Radio wave propagation in the inhomogeneous magnetic field of the solar corona

International Nuclear Information System (INIS)

Zheleznyakov, V.V.; Zlotnik, E.Ya.

1977-01-01

Various types of linear coupling between ordinary and extra-ordinary waves in the coronal plasma with the inhomogeneous magnetic field and the effect of this phenomenon upon the polarization characteristics of solar radio emission are considered. A qualitative analysis of the wave equation indicates that in a rarefied plasma the coupling effects can be displayed in a sufficiently weak magnetic field or at the angles between the magnetic field and the direction of wave propagation close enough to zero or π/2. The wave coupling parameter are found for these three cases. The radio wave propagation through the region with a quasi-transverse magnetic field and through the neutral current sheet is discussed more in detail. A qualitative picture of coupling in such a layer is supported by a numerical solution of the ''quasi-isotropic approximation'' equations. The role of the coupling effects in formation of polarization characteristics of different components of solar radio emission has been investigated. For cm wave range, the polarization is essentially dependent on the conditions in the region of the transverse magnetic field

Scalar evolution equations for shear waves in incompressible solids: a simple derivation of the Z, ZK, KZK and KP equations

KAUST Repository

Destrade, M.

2010-12-08

We study the propagation of two-dimensional finite-amplitude shear waves in a nonlinear pre-strained incompressible solid, and derive several asymptotic amplitude equations in a simple, consistent and rigorous manner. The scalar Zabolotskaya (Z) equation is shown to be the asymptotic limit of the equations of motion for all elastic generalized neo-Hookean solids (with strain energy depending only on the first principal invariant of Cauchy-Green strain). However, we show that the Z equation cannot be a scalar equation for the propagation of two-dimensional shear waves in general elastic materials (with strain energy depending on the first and second principal invariants of strain). Then, we introduce dispersive and dissipative terms to deduce the scalar Kadomtsev-Petviashvili (KP), Zabolotskaya-Khokhlov (ZK) and Khokhlov- Zabolotskaya-Kuznetsov (KZK) equations of incompressible solid mechanics. © 2010 The Royal Society.

Scalar evolution equations for shear waves in incompressible solids: a simple derivation of the Z, ZK, KZK and KP equations

KAUST Repository

Destrade, M.; Goriely, A.; Saccomandi, G.

2010-01-01

We study the propagation of two-dimensional finite-amplitude shear waves in a nonlinear pre-strained incompressible solid, and derive several asymptotic amplitude equations in a simple, consistent and rigorous manner. The scalar Zabolotskaya (Z) equation is shown to be the asymptotic limit of the equations of motion for all elastic generalized neo-Hookean solids (with strain energy depending only on the first principal invariant of Cauchy-Green strain). However, we show that the Z equation cannot be a scalar equation for the propagation of two-dimensional shear waves in general elastic materials (with strain energy depending on the first and second principal invariants of strain). Then, we introduce dispersive and dissipative terms to deduce the scalar Kadomtsev-Petviashvili (KP), Zabolotskaya-Khokhlov (ZK) and Khokhlov- Zabolotskaya-Kuznetsov (KZK) equations of incompressible solid mechanics. © 2010 The Royal Society.

Influence of Sea Surface Roughness on the Electromagnetic Wave Propagation in the Duct Environment

OpenAIRE

Zhao, X.; Huang, S.

2010-01-01

This paper deals with a study of the influence of sea surface roughness on the electromagnetic wave propagation in the duct environment. The problem of electromagnetic wave propagation is modeled by using the parabolic equation method. The roughness of the sea surface is computed by modifying the smooth surface Fresnel reflection coefficient to account for the reduction in the specular reflection due to the roughness resulting from sea wind speed. The propagation model is solved by the mixed ...

Nonlinear propagation of intense electromagnetic waves in weakly-ionized plasmas

International Nuclear Information System (INIS)

Shukla, P.K.

1993-01-01

The nonlinear propagation of intense electromagnetic waves in weakly-ionized plasmas is considered. Stimulated scattering mechanisms involving electromagnetic and acoustic waves in an unmagnetized plasma are investigated. The growth rate and threshold for three-wave decay interactions as well as modulational and filamentation instabilities are presented. Furthermore, the electromagnetic wave modulation theory is generalized for weakly ionized collisional magnetoplasmas. Here, the radiation envelope is generally governed by a nonlinear Schroedinger equation. Accounting for the dependence of the attachment frequency on the radiation intensity, ponderomotive force, as well as the differential Joule heating nonlinearity, the authors derive the equations for the nonthermal electron density and temperature perturbations. The various nonlinear terms in the electron motion are compared. The problems of self-focusing and wave localization are discussed. The relevance of the investigation to ionospheric modification by powerful electromagnetic waves is pointed out

The surface effect on axisymmetric wave propagation in piezoelectric cylindrical shells

Directory of Open Access Journals (Sweden)

Yunying Zhou

2015-02-01

Full Text Available Based on the surface piezoelectricity theory and first-order shear deformation theory, the surface effect on the axisymmetric wave propagating in piezoelectric cylindrical shells is analyzed. The Gurtin–Murdoch theory is utilized to get the nontraditional boundary conditions and constitutive equations of the surface, in company with classical governing equations of the bulk, from which the basic formulations are obtained. Numerical results show that the surface layer has a profound effect on wave characteristics in nanostructure at a higher mode.

Statistical characterization of wave propagation in mine environments

KAUST Repository

Bakir, Onur

2012-07-01

A computational framework for statistically characterizing electromagnetic (EM) wave propagation through mine tunnels and galleries is presented. The framework combines a multi-element probabilistic collocation (ME-PC) method with a novel domain-decomposition (DD) integral equation-based EM simulator to obtain statistics of electric fields due to wireless transmitters in realistic mine environments. © 2012 IEEE.

Differential equation for Alfven ion cyclotron waves in finite-length plasma

International Nuclear Information System (INIS)

Watson, D.C.; Fateman, R.J.; Baldwin, D.E.

1977-01-01

One finds the fourth-order differential equation describing an Alfven-ion-cyclotron wave propagating along a magnetic field of varying intensity. The equation is self-adjoint and possesses non-trivial turning points. The final form of the equation is checked using MACSYMA, a system for performing algebra on a computer

On the Generalized Maxwell Equations and Their Prediction of Electroscalar Wave

Directory of Open Access Journals (Sweden)

Arbab A. I.

2009-04-01

Full Text Available We have formulated the basic laws of electromagnetic theory in quaternion form. The formalism shows that Maxwell equations and Lorentz force are derivable from just one quaternion equation that only requires the Lorentz gauge. We proposed a quaternion form of the continuity equation from which we have derived the ordinary continuity equation. We introduce new transformations that produces a scalar wave and generalize the continuity equation to a set of three equations. These equations imply that both current and density are waves. Moreover, we have shown that the current can not cir- culate around a point emanating from it. Maxwell equations are invariant under these transformations. An electroscalar wave propagating with speed of light is derived upon requiring the invariance of the energy conservation equation under the new transforma- tions. The electroscalar wave function is found to be proportional to the electric field component along the charged particle motion. This scalar wave exists with or without considering the Lorentz gauge. We have shown that the electromagnetic fields travel with speed of light in the presence or absence of free charges.

Electromagnetic wave propagating along a space curve

Science.gov (United States)

Lai, Meng-Yun; Wang, Yong-Long; Liang, Guo-Hua; Wang, Fan; Zong, Hong-Shi

2018-03-01

By using the thin-layer approach, we derive the effective equation for the electromagnetic wave propagating along a space curve. We find intrinsic spin-orbit, extrinsic spin-orbit, and extrinsic orbital angular-momentum and intrinsic orbital angular-momentum couplings induced by torsion, which can lead to geometric phase, spin, and orbital Hall effects. And we show the helicity inversion induced by curvature that can convert a right-handed circularly polarized electromagnetic wave into a left-handed polarized one, vice versa. Finally, we demonstrate that the gauge invariance of the effective dynamics is protected by the geometrically induced gauge potential.

Perfectly matched layers for radio wave propagation in inhomogeneous magnetized plasmas

International Nuclear Information System (INIS)

Gondarenko, Natalia A.; Guzdar, Parvez N.; Ossakow, Sidney L.; Bernhardt, Paul A.

2004-01-01

We present 1D and 2D numerical models of the propagation of high-frequency (HF) radio waves in inhomogeneous magnetized plasmas. The simulations allow one to describe the process of linear conversion of HF electromagnetic waves into electrostatic waves. The waves, launched from the lower boundary normally or at a specified angle on a layer of a magnetoactive plasma, can undergo linear conversion of the incident O-mode into a Z-mode at appropriate locations in an inhomogeneous prescribed plasma density. The numerical scheme for solving 2D HF wave propagation equations is described. The model employed the Maxwellian perfectly matched layers (PML) technique for approximating nonreflecting boundary conditions. Our numerical studies demonstrate the effectiveness of the PML technique for transparent boundary conditions for an open-domain problem

A theoretical analysis of the weak shock waves propagating through a bubbly flow

International Nuclear Information System (INIS)

Jun, Gu Sik; Kim, Heuy Dong; Baek, Seung Cheol

2004-01-01

Two-phase flow of liquid and gas through pipe lines are frequently encountered in nuclear power plant or industrial facility. Pressure waves which can be generated by a valve operation or any other cause in pipe lines propagate through the two-phase flow, often leading to severe noise and vibration problems or fatigue failure of pipe line system. It is of practical importance to predict the propagation characteristics of the pressure waves for the safety design for the pipe line. In the present study, a theoretical analysis is performed to understand the propagation characteristics of a weak shock wave in a bubbly flow. A wave equation is developed using a small perturbation method to analyze the weak shock wave through a bubbly flow with comparably low void fractions. It is known that the elasticity of pipe and void fraction significantly affect the propagation speed of shock wave, but the frequency of relaxation oscillation which is generated behind the shock wave is not strongly influenced by the elasticity of pipe. The present analytical results are in close agreement with existing experimental data

Wave Propagation in Jointed Geologic Media

Energy Technology Data Exchange (ETDEWEB)

Antoun, T

2009-12-17

Predictive modeling capabilities for wave propagation in a jointed geologic media remain a modern day scientific frontier. In part this is due to a lack of comprehensive understanding of the complex physical processes associated with the transient response of geologic material, and in part it is due to numerical challenges that prohibit accurate representation of the heterogeneities that influence the material response. Constitutive models whose properties are determined from laboratory experiments on intact samples have been shown to over-predict the free field environment in large scale field experiments. Current methodologies for deriving in situ properties from laboratory measured properties are based on empirical equations derived for static geomechanical applications involving loads of lower intensity and much longer durations than those encountered in applications of interest involving wave propagation. These methodologies are not validated for dynamic applications, and they do not account for anisotropic behavior stemming from direcitonal effects associated with the orientation of joint sets in realistic geologies. Recent advances in modeling capabilities coupled with modern high performance computing platforms enable physics-based simulations of jointed geologic media with unprecedented details, offering a prospect for significant advances in the state of the art. This report provides a brief overview of these modern computational approaches, discusses their advantages and limitations, and attempts to formulate an integrated framework leading to the development of predictive modeling capabilities for wave propagation in jointed and fractured geologic materials.

Three-Dimensional Coupled NLS Equations for Envelope Gravity Solitary Waves in Baroclinic Atmosphere and Modulational Instability

Directory of Open Access Journals (Sweden)

Baojun Zhao

2018-01-01

Full Text Available Envelope gravity solitary waves are an important research hot spot in the field of solitary wave. And the weakly nonlinear model equations system is a part of the research of envelope gravity solitary waves. Because of the lack of technology and theory, previous studies tried hard to reduce the variable numbers and constructed the two-dimensional model in barotropic atmosphere and could only describe the propagation feature in a direction. But for the propagation of envelope gravity solitary waves in real ocean ridges and atmospheric mountains, the three-dimensional model is more appropriate. Meanwhile, the baroclinic problem of atmosphere is also an inevitable topic. In the paper, the three-dimensional coupled nonlinear Schrödinger (CNLS equations are presented to describe the evolution of envelope gravity solitary waves in baroclinic atmosphere, which are derived from the basic dynamic equations by employing perturbation and multiscale methods. The model overcomes two disadvantages: (1 baroclinic problem and (2 propagation path problem. Then, based on trial function method, we deduce the solution of the CNLS equations. Finally, modulational instability of wave trains is also discussed.

Nonlinear propagation of electromagnetic waves in negative-refraction-index composite materials.

Science.gov (United States)

Kourakis, I; Shukla, P K

2005-07-01

We investigate the nonlinear propagation of electromagnetic waves in left-handed materials. For this purpose, we consider a set of coupled nonlinear Schrödinger (CNLS) equations, which govern the dynamics of coupled electric and magnetic field envelopes. The CNLS equations are used to obtain a nonlinear dispersion, which depicts the modulational stability profile of the coupled plane-wave solutions in left-handed materials. An exact (in)stability criterion for modulational interactions is derived, and analytical expressions for the instability growth rate are obtained.

Wave propagation in magneto-electro-elastic nanobeams via two nonlocal beam models

Science.gov (United States)

Ma, Li-Hong; Ke, Liao-Liang; Wang, Yi-Ze; Wang, Yue-Sheng

2017-02-01

This paper makes the first attempt to investigate the dispersion behavior of waves in magneto-electro-elastic (MEE) nanobeams. The Euler nanobeam model and Timoshenko nanobeam model are developed in the formulation based on the nonlocal theory. By using the Hamilton's principle, we derive the governing equations which are then solved analytically to obtain the dispersion relations of MEE nanobeams. Results are presented to highlight the influences of the thermo-electro-magnetic loadings and nonlocal parameter on the wave propagation characteristics of MEE nanobeams. It is found that the thermo-electro-magnetic loadings can lead to the occurrence of the cut-off wave number below which the wave can't propagate in MEE nanobeams.

Acoustic Wave Propagation Modeling by a Two-dimensional Finite-difference Summation-by-parts Algorithm

Energy Technology Data Exchange (ETDEWEB)

Kim, K. [Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States); Petersson, N. A. [Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States); Rodgers, A. [Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)

2016-10-25

Acoustic waveform modeling is a computationally intensive task and full three-dimensional simulations are often impractical for some geophysical applications such as long-range wave propagation and high-frequency sound simulation. In this study, we develop a two-dimensional high-order accurate finite-difference code for acoustic wave modeling. We solve the linearized Euler equations by discretizing them with the sixth order accurate finite difference stencils away from the boundary and the third order summation-by-parts (SBP) closure near the boundary. Non-planar topographic boundary is resolved by formulating the governing equation in curvilinear coordinates following the interface. We verify the implementation of the algorithm by numerical examples and demonstrate the capability of the proposed method for practical acoustic wave propagation problems in the atmosphere.

An Improved Split-Step Wavelet Transform Method for Anomalous Radio Wave Propagation Modelling

Directory of Open Access Journals (Sweden)

A. Iqbal

2014-12-01

Full Text Available Anomalous tropospheric propagation caused by ducting phenomenon is a major problem in wireless communication. Thus, it is important to study the behavior of radio wave propagation in tropospheric ducts. The Parabolic Wave Equation (PWE method is considered most reliable to model anomalous radio wave propagation. In this work, an improved Split Step Wavelet transform Method (SSWM is presented to solve PWE for the modeling of tropospheric propagation over finite and infinite conductive surfaces. A large number of numerical experiments are carried out to validate the performance of the proposed algorithm. Developed algorithm is compared with previously published techniques; Wavelet Galerkin Method (WGM and Split-Step Fourier transform Method (SSFM. A very good agreement is found between SSWM and published techniques. It is also observed that the proposed algorithm is about 18 times faster than WGM and provide more details of propagation effects as compared to SSFM.

Scattering for wave equations with dissipative terms in layered media

Directory of Open Access Journals (Sweden)

Mitsuteru Kadowaki

2011-05-01

Full Text Available In this article, we show the existence of scattering solutions to wave equations with dissipative terms in layered media. To analyze the wave propagation in layered media, it is necessary to handle singular points called thresholds in the spectrum. Our main tools are Kato's smooth perturbation theory and some approximate operators.

Propagation and scattering of electromagnetic waves by the ionospheric irregularities

International Nuclear Information System (INIS)

Ho, A.Y.; Kuo, S.P.; Lee, M.C.

1993-01-01

The problem of wave propagation and scattering in the ionosphere is particularly important in the areas of communications, remote-sensing and detection. The ionosphere is often perturbed with coherently structured (quasiperiodic) density irregularities. Experimental observations suggest that these irregularities could give rise to significant ionospheric effect on wave propagation such as causing spread-F of the probing HF sounding signals and scintillation of beacon satellite signals. It was show by the latter that scintillation index S 4 ∼ 0.5 and may be as high as 0.8. In this work a quasi-particle theory is developed to study the scintillation phenomenon. A Wigner distribution function for the wave intensity in the (k,r) space is introduced and its governing equation is derived with an effective collision term giving rise to the attenuation and scattering of the wave. This kinetic equation leads to a hierarchy of moment equations in r space. This systems of equations is then truncated to the second moment which is equivalent to assuming a cold quasi-particle distribution In this analysis, the irregularities are modeled as a two dimensional density modulation on an uniform background plasma. The analysis shows that this two dimensional density grating, effectively modulates the intensity of the beacon satellite signals. This spatial modulation of the wave intensity is converted into time modulation due to the drift of the ionospheric irregularities, which then contributes to the scintillation of the beacon satellite signals. Using the proper plasma parameters and equatorial measured data of irregularities, it is shown that the scintillation index defined by S4=( 2 >- 2 )/ 2 where stands for spatial average over an irregularity wavelength is in the range of the experimentally detected values

Propagation of internal gravity waves in the inhomogeneous atmosphere

International Nuclear Information System (INIS)

Deminov, M.G.; Ponomareva, L.I.

1988-01-01

Equations for disturbances of the density, temperature and speed of large-scale horizontally propagating internal gravity wave (IGM) wind are presented with regard to non-linearity, dispersion, molecular viscosity, thermal conductivity and background horizontal density and wind speed gradients. It is shown that values of wind speed and background atmosphere density decrease, typical of night conditions, provide for IGV amplitude increase near 250 km above the equator about 1.5 times, which with regard to the both hemispheres, fully compensates the effect of viscosity and thermal conductivity under increased solar activity. Speed and density decrease along IGW propagation can be provided both by background distribution of thermosphere parameters and by the front of a large-scale IGW on the background of which isolated IGW amplitude can grow

Modal analysis of wave propagation in dispersive media

Science.gov (United States)

Abdelrahman, M. Ismail; Gralak, B.

2018-01-01

Surveys on wave propagation in dispersive media have been limited since the pioneering work of Sommerfeld [Ann. Phys. 349, 177 (1914), 10.1002/andp.19143491002] by the presence of branches in the integral expression of the wave function. In this article a method is proposed to eliminate these critical branches and hence to establish a modal expansion of the time-dependent wave function. The different components of the transient waves are physically interpreted as the contributions of distinct sets of modes and characterized accordingly. Then, the modal expansion is used to derive a modified analytical expression of the Sommerfeld precursor improving significantly the description of the amplitude and the oscillating period up to the arrival of the Brillouin precursor. The proposed method and results apply to all waves governed by the Helmholtz equations.

Modeling of shock wave propagation in large amplitude ultrasound.

Science.gov (United States)

Pinton, Gianmarco F; Trahey, Gregg E

2008-01-01

The Rankine-Hugoniot relation for shock wave propagation describes the shock speed of a nonlinear wave. This paper investigates time-domain numerical methods that solve the nonlinear parabolic wave equation, or the Khokhlov-Zabolotskaya-Kuznetsov (KZK) equation, and the conditions they require to satisfy the Rankine-Hugoniot relation. Two numerical methods commonly used in hyperbolic conservation laws are adapted to solve the KZK equation: Godunov's method and the monotonic upwind scheme for conservation laws (MUSCL). It is shown that they satisfy the Rankine-Hugoniot relation regardless of attenuation. These two methods are compared with the current implicit solution based method. When the attenuation is small, such as in water, the current method requires a degree of grid refinement that is computationally impractical. All three numerical methods are compared in simulations for lithotripters and high intensity focused ultrasound (HIFU) where the attenuation is small compared to the nonlinearity because much of the propagation occurs in water. The simulations are performed on grid sizes that are consistent with present-day computational resources but are not sufficiently refined for the current method to satisfy the Rankine-Hugoniot condition. It is shown that satisfying the Rankine-Hugoniot conditions has a significant impact on metrics relevant to lithotripsy (such as peak pressures) and HIFU (intensity). Because the Godunov and MUSCL schemes satisfy the Rankine-Hugoniot conditions on coarse grids, they are particularly advantageous for three-dimensional simulations.

Book Review: Wave propagation in materials and structures

Science.gov (United States)

Ferguson, Neil

2018-02-01

This book's remit is to provide a very extensive and detailed coverage of many one and two dimensional wave propagating behaviours primarily in structures such as rods, beams and plates of complexity covering laminated, sandwich plates, smart configurations and complex material compositions. This is potentially where the detailed presentation, including the derivation of the governing equations of motion from first principles, i.e. Hamilton's method, for example, distracts slightly from the subsequent wave solutions, the numerical simulations showing time responses, the wave speeds and importantly the dispersion characteristics. The author introduces a number of known analytical methodologies and means to obtain wave solutions, including the spectral finite element approach and also provides numerical examples showing the approach being applied to joints and framed structures.

THE BASIS OF MATHEMATICAL DESCRIPTION FOR WAVE MODEL OF STRESSES PROPAGATION IN RAILWAY TRACK

Directory of Open Access Journals (Sweden)

D. M. Kurhan

2016-10-01

Full Text Available Purpose. Modern scientific research has repeatedly cited practical examples of the dynamic effects of railway track operation that go beyond the static calculation schemes. For the track sections where the train speed is approaching to the velocity of wave propagation in the slab track layers such issues are of particular relevance. An adequate tool for the study of such issues can be the use of the wave theory of stress propagation. The purpose of the article is the creation of a mathematical description of the basic principles of the stress propagation wave model in the railway track, which can be used as a basis for the practical development of the relevant calculation system. Methodology. The model of stress-strain states of the railway track on the basis of the stress wave propagation theory is to bring together the equations of the geometry of the outline of the space systems that is involved in the interaction at a given time, and the dynamic equilibrium equations of deformation. The solution is based on the use of the laws of the theory of elasticity. The wave front is described by an ellipsoid equation. When determining the variation in time of the surface position of the ellipsoid a vector approach is used. Findings. The geometry equations of the wave motion determine the volumes of material layers of the slab track involved in the interaction at a given time. The dynamic equilibrium determination of the deformed condition of the space bounded by the wave front makes it possible to calculate both the stresses and strains, and their changes during the time of the load perception. Thus, mathematical descriptions of the processes that occur in the perception of the load by the elements of railway track at high speeds were obtained. Originality. The simulation tasks of the track and rolling stock interaction, in particular taking into account the dynamic deflection of slab track were further developed. For the first time the article

A two dimension model of the uterine electrical wave propagation.

Science.gov (United States)

Rihana, S; Lefrançois, E; Marque, C

2007-01-01

The uterus, usually quiescent during pregnancy, exhibits forceful contractions at term leading to delivery. These contractions are caused by the synchronized propagation of electrical waves from the pacemaker cells to its neighbors inducing the whole coordinated contraction of the uterus wall leading to labor. In a previous work, we simulate the electrical activity of a single uterine cell by a set of ordinary differential equations. Then, this model has been used to simulate the electrical activity propagation. In the present work, the uterine cell tissue is assumed to have uniform and isotropic propagation, and constant electrical membrane properties. The stability of the numerical solution imposes the choice of a critical temporal step. A wave starts at a pacemaker cell; this electrical activity is initiated by the injection of an external stimulation current to the cell membrane. We observe synchronous wave propagation for axial resistance values around 0.5 GOmega or less and propoagation blocking for values greater than 0.7 GOmega. We compute the conduction velocity of the excitation, for different axial resistance values, and obtain a velocity about 10 cm/sec, approaching the one described by the literature for the rat at end of term.

Simulating Seismic Wave Propagation in Viscoelastic Media with an Irregular Free Surface

Science.gov (United States)

Liu, Xiaobo; Chen, Jingyi; Zhao, Zhencong; Lan, Haiqiang; Liu, Fuping

2018-05-01

In seismic numerical simulations of wave propagation, it is very important for us to consider surface topography and attenuation, which both have large effects (e.g., wave diffractions, conversion, amplitude/phase change) on seismic imaging and inversion. An irregular free surface provides significant information for interpreting the characteristics of seismic wave propagation in areas with rugged or rapidly varying topography, and viscoelastic media are a better representation of the earth's properties than acoustic/elastic media. In this study, we develop an approach for seismic wavefield simulation in 2D viscoelastic isotropic media with an irregular free surface. Based on the boundary-conforming grid method, the 2D time-domain second-order viscoelastic isotropic equations and irregular free surface boundary conditions are transferred from a Cartesian coordinate system to a curvilinear coordinate system. Finite difference operators with second-order accuracy are applied to discretize the viscoelastic wave equations and the irregular free surface in the curvilinear coordinate system. In addition, we select the convolutional perfectly matched layer boundary condition in order to effectively suppress artificial reflections from the edges of the model. The snapshot and seismogram results from numerical tests show that our algorithm successfully simulates seismic wavefields (e.g., P-wave, Rayleigh wave and converted waves) in viscoelastic isotropic media with an irregular free surface.

Simulation of the acoustic wave propagation using a meshless method

Directory of Open Access Journals (Sweden)

Bajko J.

2017-01-01

Full Text Available This paper presents numerical simulations of the acoustic wave propagation phenomenon modelled via Linearized Euler equations. A meshless method based on collocation of the strong form of the equation system is adopted. Moreover, the Weighted least squares method is used for local approximation of derivatives as well as stabilization technique in a form of spatial ltering. The accuracy and robustness of the method is examined on several benchmark problems.

Propagation of sound waves in ducts

DEFF Research Database (Denmark)

Jacobsen, Finn

2000-01-01

Plane wave propagation in ducts with rigid walls, radiation from ducts, classical four-pole theory for composite duct systems, and three-dimentional waves in wave guides of various cross-sectional shape are described.......Plane wave propagation in ducts with rigid walls, radiation from ducts, classical four-pole theory for composite duct systems, and three-dimentional waves in wave guides of various cross-sectional shape are described....

Effective constants for wave propagation through partially saturated porous media

International Nuclear Information System (INIS)

Berryman, J.G.; Thigpen, L.

1985-01-01

The multipole scattering coefficients for elastic wave scattering from a spherical inhomogeneity in a fluid-saturated porous medium have been calculated. These coefficients may be used to obtain estimates of the effective macroscopic constants for long-wavelength propagation of elastic waves through partially saturated media. If the volume average of the single scattering from spherical bubbles of gas and liquid is required to vanish, the resulting equations determine the effective bulk modulus, density, and viscosity of the multiphase fluid filling the pores. The formula for the effective viscosity during compressional wave excitation is apparently new

Propagation of extensional waves in a piezoelectric semiconductor rod

Directory of Open Access Journals (Sweden)

C.L. Zhang

2016-04-01

Full Text Available We studied the propagation of extensional waves in a thin piezoelectric semiconductor rod of ZnO whose c-axis is along the axis of the rod. The macroscopic theory of piezoelectric semiconductors was used which consists of the coupled equations of piezoelectricity and the conservation of charge. The problem is nonlinear because the drift current is the product of the unknown electric field and the unknown carrier density. A perturbation procedure was used which resulted in two one-way coupled linear problems of piezoelectricity and the conservation of charge, respectively. The acoustic wave and the accompanying electric field were obtained from the equations of piezoelectricity. The motion of carriers was then determined from the conservation of charge using a trigonometric series. It was found that while the acoustic wave was approximated by a sinusoidal wave, the motion of carriers deviates from a sinusoidal wave qualitatively because of the contributions of higher harmonics arising from the originally nonlinear terms. The wave crests become higher and sharper while the troughs are shallower and wider. This deviation is more pronounced for acoustic waves with larger amplitudes.

The numerical simulation of Lamb wave propagation in laser welding of stainless steel

Science.gov (United States)

Zhang, Bo; Liu, Fang; Liu, Chang; Li, Jingming; Zhang, Baojun; Zhou, Qingxiang; Han, Xiaohui; Zhao, Yang

2017-12-01

In order to explore the Lamb wave propagation in laser welding of stainless steel, the numerical simulation is used to show the feature of Lamb wave. In this paper, according to Lamb dispersion equation, excites the Lamb wave on the edge of thin stainless steel plate, and presents the reflection coefficient for quantizing the Lamb wave energy, the results show that the reflection coefficient is increased with the welding width increasing,

Terrestrial propagation of long electromagnetic waves

CERN Document Server

Galejs, Janis; Fock, V A

2013-01-01

Terrestrial Propagation of Long Electromagnetic Waves deals with the propagation of long electromagnetic waves confined principally to the shell between the earth and the ionosphere, known as the terrestrial waveguide. The discussion is limited to steady-state solutions in a waveguide that is uniform in the direction of propagation. Wave propagation is characterized almost exclusively by mode theory. The mathematics are developed only for sources at the ground surface or within the waveguide, including artificial sources as well as lightning discharges. This volume is comprised of nine chapte

Propagation of edge waves in a thinly layered laminated medium with stress couples under initial stresses

Directory of Open Access Journals (Sweden)

Pijush Pal Roy

1987-01-01

Full Text Available The propagation of edge waves in a thinly layered laminated medium with stress couples under initial stresses is examined. Based upon an approximate representation of a laminated medium by an equivalent anisotropic continuum with average initial and couple stresses, an explicit form of frequency equation is obtained to derive the phase velocity of edge waves. Edge waves exist under certain conditions. The inclusion of couple stresses increases the velocity of wave propagation. For a specific compression, the presence of couple stresses increases the velocity of wave propagation with the increase of wave number, whereas the reverse is the case when there is no couple stress. Numerical computation is performed with graphical representations. Several special cases are also examined.

A consistent formulation of wave propagation and conversion in low aspect ratio tokamaks with non-circular cross section

International Nuclear Information System (INIS)

Cuperman, S.; Bruma, C.; Komoshvili, K.

1999-01-01

The authors developed a consistent formalism for the full wave equation, appropriate for the study of propagation, absorption and wave conversion of externally launched waves in strongly toroidal, spherical tokamaks with non-circular cross-section. This includes also the formulation of rigorous regularity, boundary, gauge and periodicity conditions suitable for the exact solution of the wave equation for such devices

Sound Propagation Around Off-Shore Wind Turbines. Long-Range Parabolic Equation Calculations for Baltic Sea Conditions

Energy Technology Data Exchange (ETDEWEB)

Johansson, Lisa

2003-07-01

Low-frequency, long-range sound propagation over a sea surface has been calculated using a wide-angel Cranck-Nicholson Parabolic Equation method. The model is developed to investigate noise from off-shore wind turbines. The calculations are made using normal meteorological conditions of the Baltic Sea. Special consideration has been made to a wind phenomenon called low level jet with strong winds on rather low altitude. The effects of water waves on sound propagation have been incorporated in the ground boundary condition using a boss model. This way of including roughness in sound propagation models is valid for water wave heights that are small compared to the wave length of the sound. Nevertheless, since only low frequency sound is considered, waves up to the mean wave height of the Baltic Sea can be included in this manner. The calculation model has been tested against benchmark cases and agrees well with measurements. The calculations show that channelling of sound occurs at downwind conditions and that the sound propagation tends towards cylindrical spreading. The effects of the water waves are found to be fairly small.

Asymptotic solutions and spectral theory of linear wave equations

International Nuclear Information System (INIS)

Adam, J.A.

1982-01-01

This review contains two closely related strands. Firstly the asymptotic solution of systems of linear partial differential equations is discussed, with particular reference to Lighthill's method for obtaining the asymptotic functional form of the solution of a scalar wave equation with constant coefficients. Many of the applications of this technique are highlighted. Secondly, the methods and applications of the theory of the reduced (one-dimensional) wave equation - particularly spectral theory - are discussed. While the breadth of application and power of the techniques is emphasised throughout, the opportunity is taken to present to a wider readership, developments of the methods which have occured in some aspects of astrophysical (particularly solar) and geophysical fluid dynamics. It is believed that the topics contained herein may be of relevance to the applied mathematician or theoretical physicist interest in problems of linear wave propagation in these areas. (orig./HSI)

Investigation into stress wave propagation in metal foams

Directory of Open Access Journals (Sweden)

Li Lang

2015-01-01

Full Text Available The aim of this study is to investigate stress wave propagation in metal foams under high-speed impact loading. Three-dimensional Voronoi model is established to represent real closed-cell foam. Based on the one-dimensional stress wave theory and Voronoi model, a numerical model is developed to calculate the velocity of elastic wave and shock wave in metal foam. The effects of impact velocity and relative density of metal foam on the stress wave propagation in metal foams are explored respectively. The results show that both elastic wave and shock wave propagate faster in metal foams with larger relative density; with increasing the impact velocity, the shock wave propagation velocity increase, but the elastic wave propagation is not sensitive to the impact velocity.

Model for small arms fire muzzle blast wave propagation in air

Science.gov (United States)

Aguilar, Juan R.; Desai, Sachi V.

2011-11-01

Accurate modeling of small firearms muzzle blast wave propagation in the far field is critical to predict sound pressure levels, impulse durations and rise times, as functions of propagation distance. Such a task being relevant to a number of military applications including the determination of human response to blast noise, gunfire detection and localization, and gun suppressor design. Herein, a time domain model to predict small arms fire muzzle blast wave propagation is introduced. The model implements a Friedlander wave with finite rise time which diverges spherically from the gun muzzle. Additionally, the effects in blast wave form of thermoviscous and molecular relaxational processes, which are associated with atmospheric absorption of sound were also incorporated in the model. Atmospheric absorption of blast waves is implemented using a time domain recursive formula obtained from numerical integration of corresponding differential equations using a Crank-Nicholson finite difference scheme. Theoretical predictions from our model were compared to previously recorded real world data of muzzle blast wave signatures obtained by shooting a set different sniper weapons of varying calibers. Recordings containing gunfire acoustical signatures were taken at distances between 100 and 600 meters from the gun muzzle. Results shows that predicted blast wave slope and exponential decay agrees well with measured data. Analysis also reveals the persistency of an oscillatory phenomenon after blast overpressure in the recorded wave forms.

Modelling Acoustic Wave Propagation in Axisymmetric Varying-Radius Waveguides

DEFF Research Database (Denmark)

Bæk, David; Willatzen, Morten

2008-01-01

A computationally fast and accurate model (a set of coupled ordinary differential equations) for fluid sound-wave propagation in infinite axisymmetric waveguides of varying radius is proposed. The model accounts for fluid heat conduction and fluid irrotational viscosity. The model problem is solved...... by expanding solutions in terms of cross-sectional eigenfunctions following Stevenson’s method. A transfer matrix can be easily constructed from simple model responses of a given waveguide and later used in computing the response to any complex wave input. Energy losses due to heat conduction and viscous...

The wave equation on a curved space-time

International Nuclear Information System (INIS)

Friedlander, F.G.

1975-01-01

It is stated that chapters on differential geometry, distribution theory, and characteristics and the propagation of discontinuities are preparatory. The main matter is in three chapters, entitled: fundamental solutions, representation theorems, and wave equations on n-dimensional space-times. These deal with general construction of fundamental solutions and their application to the Cauchy problem. (U.K.)

Relativistic wave equations without the Velo-Zwanziger pathology

International Nuclear Information System (INIS)

Khalil, M.A.K.

1976-06-01

For particles described by relativistic wave equations of the form: (-iGAMMA x delta + m) psi(x) = 0 interacting with an external field B(x) it is known that the ''noncausal'' propagation characteristics are not present when (1) GAMMA 0 is diagonalizable and (2) B(x) = -eGAMMA/sub mu/A/sup mu/(x) (Amar--Dozzio). The ''noncausality''difficulties arise for the Rarita--Schwinger spin 3 / 2 equation, with nondiagonalizable GAMMA 0 , in minimal coupling (i.e., B(x) = -eGAMMA x A(x)) and the PDK spin 1 equation, with diagonalizable GAMMA 0 , in a quadrupole coupling (Velo--Zwanziger) where either (1) or (2) of the Amar--Dozzio (sufficient) conditions are violated. Some sufficient conditions are derived and explored where the Velo--Zwanziger ''noncausality'' pathology can be avoided, even though one, or the other, or both of the conditions (1) and (2) are violated. Examples with both reducible and irreducible wave equations are included

Influence of Sea Surface Roughness on the Electromagnetic Wave Propagation in the Duct Environment

Directory of Open Access Journals (Sweden)

X. Zhao

2010-12-01

Full Text Available This paper deals with a study of the influence of sea surface roughness on the electromagnetic wave propagation in the duct environment. The problem of electromagnetic wave propagation is modeled by using the parabolic equation method. The roughness of the sea surface is computed by modifying the smooth surface Fresnel reflection coefficient to account for the reduction in the specular reflection due to the roughness resulting from sea wind speed. The propagation model is solved by the mixed Fourier split-step algorithm. Numerical experiments indicate that wind-driven roughened sea surface has an impact on the electromagnetic wave propagation in the duct environment, and the strength is intensified along with the increment of sea wind speeds and/or the operating frequencies. In a fixed duct environment, however, proper disposition of the transmitter could reduce these impacts.

Effect of surface wave propagation in a four-layered oceanic crust model

Science.gov (United States)

Paul, Pasupati; Kundu, Santimoy; Mandal, Dinbandhu

2017-12-01

Dispersion of Rayleigh type surface wave propagation has been discussed in four-layered oceanic crust. It includes a sandy layer over a crystalline elastic half-space and over it there are two more layers—on the top inhomogeneous liquid layer and under it a liquid-saturated porous layer. Frequency equation is obtained in the form of determinant. The effects of the width of different layers as well as the inhomogeneity of liquid layer, sandiness of sandy layer on surface waves are depicted and shown graphically by considering all possible case of the particular model. Some special cases have been deduced, few special cases give the dispersion equation of Scholte wave and Stoneley wave, some of which have already been discussed elsewhere.

Thermo-magnetic field effects on the wave propagation behavior of smart magnetostrictive sandwich nanoplates

Science.gov (United States)

Ebrahimi, Farzad; Dabbagh, Ali

2018-03-01

In this paper, a three-variable plate model is utilized to explore the wave propagation problem of smart sandwich nanoplates made of a magnetostrictive core and ceramic face sheets while subjected to thermo-magnetic loading. Herein, the magnetostriction effect is considered and controlled via a feedback control system. The nanoplate is supposed to be embedded on a visco-Pasternak elastic substrate. The kinematic relations are derived based on the Kirchhoff plate theory; also, combining these obtained equations with Hamilton's principle, the local equations of motion are achieved. According to a nonlocal strain gradient theory (NSGT), the small-scale influences are covered precisely by introducing two scale coefficients. Afterwards, the nonlocal governing equations are derived coupling the local equations with those of the NSGT. Applying an analytical solution, the wave frequency and phase velocity of the propagated waves can be gathered solving an eigenvalue problem. On the other hand, accuracy and efficiency of the presented model are verified by setting a comparison between the obtained results with those of previous published researches. Effects of different variants are plotted in some figures and the highlights are discussed in detail.

Numerical simulation of stress wave propagation from underground nuclear explosions

Energy Technology Data Exchange (ETDEWEB)

Cherry, J T; Petersen, F L [Lawrence Radiation Laboratory, University of California, Livermore, CA (United States)

1970-05-01

This paper presents a numerical model of stress wave propagation (SOC) which uses material properties data from a preshot testing program to predict the stress-induced effects on the rock mass involved in a Plowshare application. SOC calculates stress and particle velocity history, cavity radius, extent of brittle failure, and the rock's efficiency for transmitting stress. The calculations are based on an equation of state for the rock, which is developed from preshot field and laboratory measurements of the rock properties. The field measurements, made by hole logging, determine in situ values of the rock's density, water content, and propagation velocity for elastic waves. These logs also are useful in judging the layering of the rock and in choosing which core samples to test in the laboratory. The laboratory analysis of rock cores includes determination of hydrostatic compressibility to 40 kb, triaxial strength data, tensile strength, Hugoniot elastic limit, and, for the rock near the point of detonation, high-pressure Hugoniot data. Equation-of-state data are presented for rock from three sites subjected to high explosive or underground nuclear shots, including the Hardhat and Gasbuggy sites. SOC calculations of the effects of these two shots on the surrounding rock are compared with the observed effects. In both cases SOC predicts the size of the cavity quite closely. Results of the Gasbuggy calculations indicate that useful predictions of cavity size and chimney height can be made when an adequate preshot testing program is run to determine the rock's equation of state. Seismic coupling is very sensitive to the low-pressure part of the equation of state, and its successful prediction depends on agreement between the logging data and the static compressibility data. In general, it appears that enough progress has been made in calculating stress wave propagation to begin looking at derived numbers, such as number of cracks per zone, for some insight into the

Numerical simulation of stress wave propagation from underground nuclear explosions

International Nuclear Information System (INIS)

Cherry, J.T.; Petersen, F.L.

1970-01-01

This paper presents a numerical model of stress wave propagation (SOC) which uses material properties data from a preshot testing program to predict the stress-induced effects on the rock mass involved in a Plowshare application. SOC calculates stress and particle velocity history, cavity radius, extent of brittle failure, and the rock's efficiency for transmitting stress. The calculations are based on an equation of state for the rock, which is developed from preshot field and laboratory measurements of the rock properties. The field measurements, made by hole logging, determine in situ values of the rock's density, water content, and propagation velocity for elastic waves. These logs also are useful in judging the layering of the rock and in choosing which core samples to test in the laboratory. The laboratory analysis of rock cores includes determination of hydrostatic compressibility to 40 kb, triaxial strength data, tensile strength, Hugoniot elastic limit, and, for the rock near the point of detonation, high-pressure Hugoniot data. Equation-of-state data are presented for rock from three sites subjected to high explosive or underground nuclear shots, including the Hardhat and Gasbuggy sites. SOC calculations of the effects of these two shots on the surrounding rock are compared with the observed effects. In both cases SOC predicts the size of the cavity quite closely. Results of the Gasbuggy calculations indicate that useful predictions of cavity size and chimney height can be made when an adequate preshot testing program is run to determine the rock's equation of state. Seismic coupling is very sensitive to the low-pressure part of the equation of state, and its successful prediction depends on agreement between the logging data and the static compressibility data. In general, it appears that enough progress has been made in calculating stress wave propagation to begin looking at derived numbers, such as number of cracks per zone, for some insight into the

Optimal implicit 2-D finite differences to model wave propagation in poroelastic media

Science.gov (United States)

Itzá, Reymundo; Iturrarán-Viveros, Ursula; Parra, Jorge O.

2016-08-01

Numerical modeling of seismic waves in heterogeneous porous reservoir rocks is an important tool for the interpretation of seismic surveys in reservoir engineering. We apply globally optimal implicit staggered-grid finite differences (FD) to model 2-D wave propagation in heterogeneous poroelastic media at a low-frequency range (differentiation involves solving tridiagonal linear systems of equations through Thomas' algorithm.

Propagation of sound and thermal waves in an ionizing-recombining hydrogen plasma: Revision of results

International Nuclear Information System (INIS)

Di Sigalotti, Leonardo G.; Sira, Eloy; Tremola, Ciro

2002-01-01

The propagation of acoustic and thermal waves in a heat conducting, hydrogen plasma, in which photoionization and photorecombination [H + +e - H+hν(χ)] processes are progressing, is re-examined here using linear analysis. The resulting dispersion equation is solved analytically and the results are compared with previous solutions for the same plasma model. In particular, it is found that wave propagation in a slightly and highly ionized hydrogen plasma is affected by crossing between acoustic and thermal modes. At temperatures where the plasma is partially ionized, waves of all frequencies propagate without the occurrence of mode crossing. These results disagree with those reported in previous work, thereby leading to a different physical interpretation of the propagation of small linear disturbances in a conducting, ionizing-recombining, hydrogen plasma

Modes in light wave propagating in semiconductor laser

Science.gov (United States)

Manko, Margarita A.

1994-01-01

The study of semiconductor laser based on an analogy of the Schrodinger equation and an equation describing light wave propagation in nonhomogeneous medium is developed. The active region of semiconductor laser is considered as optical waveguide confining the electromagnetic field in the cross-section (x,y) and allowing waveguide propagation along the laser resonator (z). The mode structure is investigated taking into account the transversal and what is the important part of the suggested consideration longitudinal nonhomogeneity of the optical waveguide. It is shown that the Gaussian modes in the case correspond to spatial squeezing and correlation. Spatially squeezed two-mode structure of nonhomogeneous optical waveguide is given explicitly. Distribution of light among the laser discrete modes is presented. Properties of the spatially squeezed two-mode field are described. The analog of Franck-Condon principle for finding the maxima of the distribution function and the analog of Ramsauer effect for control of spatial distribution of laser emission are discussed.

Symmetries of the triple degenerate DNLS equations for weakly nonlinear dispersive MHD waves

International Nuclear Information System (INIS)

Webb, G. M.; Brio, M.; Zank, G. P.

1996-01-01

A formulation of Hamiltonian and Lagrangian variational principles, Lie point symmetries and conservation laws for the triple degenerate DNLS equations describing the propagation of weakly nonlinear dispersive MHD waves along the ambient magnetic field, in β∼1 plasmas is given. The equations describe the interaction of the Alfven and magnetoacoustic modes near the triple umbilic point, where the fast magnetosonic, slow magnetosonic and Alfven speeds coincide and a g 2 =V A 2 where a g is the gas sound speed and V A is the Alfven speed. A discussion is given of the travelling wave similarity solutions of the equations, which include solitary wave and periodic traveling waves. Strongly compressible solutions indicate the necessity for the insertion of shocks in the flow, whereas weakly compressible, near Alfvenic solutions resemble similar, shock free travelling wave solutions of the DNLS equation

Pulse-wave propagation in straight-geometry vessels for stiffness estimation: theory, simulations, phantoms and in vitro findings.

Science.gov (United States)

Shahmirzadi, Danial; Li, Ronny X; Konofagou, Elisa E

2012-11-01

Pulse wave imaging (PWI) is an ultrasound-based method for noninvasive characterization of arterial stiffness based on pulse wave propagation. Reliable numerical models of pulse wave propagation in normal and pathological aortas could serve as powerful tools for local pulse wave analysis and a guideline for PWI measurements in vivo. The objectives of this paper are to (1) apply a fluid-structure interaction (FSI) simulation of a straight-geometry aorta to confirm the Moens-Korteweg relationship between the pulse wave velocity (PWV) and the wall modulus, and (2) validate the simulation findings against phantom and in vitro results. PWI depicted and tracked the pulse wave propagation along the abdominal wall of canine aorta in vitro in sequential Radio-Frequency (RF) ultrasound frames and estimates the PWV in the imaged wall. The same system was also used to image multiple polyacrylamide phantoms, mimicking the canine measurements as well as modeling softer and stiffer walls. Finally, the model parameters from the canine and phantom studies were used to perform 3D two-way coupled FSI simulations of pulse wave propagation and estimate the PWV. The simulation results were found to correlate well with the corresponding Moens-Korteweg equation. A high linear correlation was also established between PWV² and E measurements using the combined simulation and experimental findings (R² =  0.98) confirming the relationship established by the aforementioned equation.

Effect of parallel electric fields on the whistler mode wave propagation in the magnetosphere

International Nuclear Information System (INIS)

Gupta, G.P.; Singh, R.N.

1975-01-01

The effect of parallel electric fields on whistler mode wave propagation has been studied. To account for the parallel electric fields, the dispersion equation has been analyzed, and refractive index surfaces for magnetospheric plasma have been constructed. The presence of parallel electric fields deforms the refractive index surfaces which diffuse the energy flow and produce defocusing of the whistler mode waves. The parallel electric field induces an instability in the whistler mode waves propagating through the magnetosphere. The growth or decay of whistler mode instability depends on the direction of parallel electric fields. It is concluded that the analyses of whistler wave records received on the ground should account for the role of parallel electric fields

Propagation of SLF/ELF electromagnetic waves

CERN Document Server

Pan, Weiyan

2014-01-01

This book deals with the SLF/ELF wave propagation, an important branch of electromagnetic theory. The SLF/ELF wave propagation theory is well applied in earthquake electromagnetic radiation, submarine communication, thunderstorm detection, and geophysical prospecting and diagnostics. The propagation of SLF/ELF electromagnetic waves is introduced in various media like the earth-ionospheric waveguide, ionospheric plasma, sea water, earth, and the boundary between two different media or the stratified media. Applications in the earthquake electromagnetic radiation and the submarine communications are also addressed. This book is intended for scientists and engineers in the fields of radio propagation and EM theory and applications. Prof. Pan is a professor at China Research Institute of Radiowave Propagation in Qingdao (China). Dr. Li is a professor at Zhejiang University in Hangzhou (China).

An alternative view on the role of the β-effect in the Rossby wave propagation mechanism

Directory of Open Access Journals (Sweden)

Eyal Heifetz

2014-11-01

Full Text Available The role of the β-effect in the Rossby wave propagation mechanism is examined in the linearised shallow water equations directly in momentum–height variables, without recourse to potential vorticity (PV. Rigorous asymptotic expansion of the equations, with respect to the small non-dimensionalised β parameter, reveals in detail how the Coriolis force acting on the small ageostrophic terms translates the geostrophic leading-order solution to propagate westward in concert. This information cannot be obtained directly from the conventional PV perspective on the propagation mechanism. Furthermore, a comparison between the β-effect in planetary Rossby waves and the sloping-bottom effect in promoting topographic Rossby waves shows that the ageostrophic terms play different roles in the two cases. This is despite the fact that from the PV viewpoint whether the advection of mean PV gradient is set up by changes in planetary vorticity or by mean depth is inconsequential.

Ion stochastic heating by obliquely propagating magnetosonic waves

International Nuclear Information System (INIS)

Gao Xinliang; Lu Quanming; Wu Mingyu; Wang Shui

2012-01-01

The ion motions in obliquely propagating Alfven waves with sufficiently large amplitudes have already been studied by Chen et al.[Phys. Plasmas 8, 4713 (2001)], and it was found that the ion motions are stochastic when the wave frequency is at a fraction of the ion gyro-frequency. In this paper, with test particle simulations, we investigate the ion motions in obliquely propagating magnetosonic waves and find that the ion motions also become stochastic when the amplitude of the magnetosonic waves is sufficiently large due to the resonance at sub-cyclotron frequencies. Similar to the Alfven wave, the increase of the propagating angle, wave frequency, and the number of the wave modes can lower the stochastic threshold of the ion motions. However, because the magnetosonic waves become more and more compressive with the increase of the propagating angle, the decrease of the stochastic threshold with the increase of the propagating angle is more obvious in the magnetosonic waves than that in the Alfven waves.

Wave-equation dispersion inversion of surface waves recorded on irregular topography

KAUST Repository

Li, Jing; Schuster, Gerard T.; Lin, Fan-Chi; Alam, Amir

2017-01-01

Significant topographic variations will strongly influence the amplitudes and phases of propagating surface waves. Such effects should be taken into account, otherwise the S-velocity model inverted from the Rayleigh dispersion curves will contain significant inaccuracies. We now show that the recently developed wave-equation dispersion inversion (WD) method naturally takes into account the effects of topography to give accurate S-velocity tomograms. Application of topographic WD to demonstrates that WD can accurately invert dispersion curves from seismic data recorded over variable topography. We also apply this method to field data recorded on the crest of mountainous terrain and find with higher resolution than the standard WD tomogram.

Wave-equation dispersion inversion of surface waves recorded on irregular topography

KAUST Repository

Li, Jing

2017-08-17

Significant topographic variations will strongly influence the amplitudes and phases of propagating surface waves. Such effects should be taken into account, otherwise the S-velocity model inverted from the Rayleigh dispersion curves will contain significant inaccuracies. We now show that the recently developed wave-equation dispersion inversion (WD) method naturally takes into account the effects of topography to give accurate S-velocity tomograms. Application of topographic WD to demonstrates that WD can accurately invert dispersion curves from seismic data recorded over variable topography. We also apply this method to field data recorded on the crest of mountainous terrain and find with higher resolution than the standard WD tomogram.

Engineering equations for characterizing non-linear laser intensity propagation in air with loss.

Science.gov (United States)

Karr, Thomas; Stotts, Larry B; Tellez, Jason A; Schmidt, Jason D; Mansell, Justin D

2018-02-19

The propagation of high peak-power laser beams in real atmospheres will be affected at long range by both linear and nonlinear effects contained therein. Arguably, J. H. Marburger is associated with the mathematical characterization of this phenomenon. This paper provides a validated set of engineering equations for characterizing the self-focusing distance from a laser beam propagating through non-turbulent air with, and without, loss as well as three source configurations: (1) no lens, (2) converging lens and (3) diverging lens. The validation was done against wave-optics simulation results. Some validated equations follow Marburger completely, but others do not, requiring modification of the original theory. Our results can provide a guide for numerical simulations and field experiments.

The Green-function transform and wave propagation

Directory of Open Access Journals (Sweden)

Colin eSheppard

2014-11-01

Full Text Available Fourier methods well known in signal processing are applied to three-dimensional wave propagation problems. The Fourier transform of the Green function, when written explicitly in terms of a real-valued spatial frequency, consists of homogeneous and inhomogeneous components. Both parts are necessary to result in a pure out-going wave that satisfies causality. The homogeneous component consists only of propagating waves, but the inhomogeneous component contains both evanescent and propagating terms. Thus we make a distinction between inhomogeneous waves and evanescent waves. The evanescent component is completely contained in the region of the inhomogeneous component outside the k-space sphere. Further, propagating waves in the Weyl expansion contain both homogeneous and inhomogeneous components. The connection between the Whittaker and Weyl expansions is discussed. A list of relevant spherically symmetric Fourier transforms is given.

High Frequency Waves Propagating in Octagonal Bars: a Low Cost Computation Algorithm

Directory of Open Access Journals (Sweden)

Alessandro Marzani

2009-02-01

Full Text Available In this paper a hybrid semi-analytical Finite Element formulation is proposed to efficiently calculate the time dependent response due to stress waves propagating in a slender solid with uniform cross-section when excited by impulsive forces. The formulation takes advantage of the direct and inverse Fourier transform to formulate and solve the governing wave equation. The framework is applied to an octagonal viscoelastic isotropic steel bar.

Propagation-invariant waves in acoustic, optical, and radio-wave fields

OpenAIRE

Salo, Janne

2003-01-01

The physical phenomena considered in this thesis are associated with electromagnetic and acoustic waves that propagate in free space or in homogeneous media without diffraction. The concept of rotationally periodic wave propagation is introduced in the first journal article included in the thesis and it is subsequently used to analyse waves that avoid diffractive deterioration by repeatedly returning to their initial shape, possibly rotated around the optical axis. Such waves constitute an es...

Wave propagation in fluid-conveying viscoelastic single-walled carbon nanotubes with surface and nonlocal effects

Science.gov (United States)

Zhen, Ya-Xin

2017-02-01

In this paper, the transverse wave propagation in fluid-conveying viscoelastic single-walled carbon nanotubes is investigated based on nonlocal elasticity theory with consideration of surface effect. The governing equation is formulated utilizing nonlocal Euler-Bernoulli beam theory and Kelvin-Voigt model. Explicit wave dispersion relation is developed and wave phase velocities and frequencies are obtained. The effect of the fluid flow velocity, structural damping, surface effect, small scale effects and tube diameter on the wave propagation properties are discussed with different wave numbers. The wave frequency increases with the increase of fluid flow velocity, but decreases with the increases of tube diameter and wave number. The effect of surface elasticity and residual surface tension is more significant for small wave number and tube diameter. For larger values of wave number and nonlocal parameters, the real part of frequency ratio raises.

Wave propagation through a flexoelectric piezoelectric slab sandwiched by two piezoelectric half-spaces.

Science.gov (United States)

Jiao, Fengyu; Wei, Peijun; Li, Yueqiu

2018-01-01

Reflection and transmission of plane waves through a flexoelectric piezoelectric slab sandwiched by two piezoelectric half-spaces are studied in this paper. The secular equations in the flexoelectric piezoelectric material are first derived from the general governing equation. Different from the classical piezoelectric medium, there are five kinds of coupled elastic waves in the piezoelectric material with the microstructure effects taken into consideration. The state vectors are obtained by the summation of contributions from all possible partial waves. The state transfer equation of flexoelectric piezoelectric slab is derived from the motion equation by the reduction of order, and the transfer matrix of flexoelectric piezoelectric slab is obtained by solving the state transfer equation. By using the continuous conditions at the interface and the approach of partition matrix, we get the resultant algebraic equations in term of the transfer matrix from which the reflection and transmission coefficients can be calculated. The amplitude ratios and further the energy flux ratios of various waves are evaluated numerically. The numerical results are shown graphically and are validated by the energy conservation law. Based on these numerical results, the influences of two characteristic lengths of microstructure and the flexoelectric coefficients on the wave propagation are discussed. Copyright © 2017 Elsevier B.V. All rights reserved.

''Free-space'' boundary conditions for the time-dependent wave equation

International Nuclear Information System (INIS)

Lindman, E.L.

1975-01-01

Boundary conditions for the discrete wave equation which act like an infinite region of free space in contact with the computational region can be constructed using projection operators. Propagating and evanescent waves coming from within the computational region generate no reflected waves as they cross the boundary. At the same time arbitrary waves may be launched into the computational region. Well known projection operators for one-dimensional waves may be used for this purpose in one dimension. Extensions of these operators to higher dimensions along with numerically efficient approximations to them are described for higher-dimensional problems. The separation of waves into ingoing and outgoing waves inherent in these boundary conditions greatly facilitates diagnostics

Solitary wave for a nonintegrable discrete nonlinear Schrödinger equation in nonlinear optical waveguide arrays

Science.gov (United States)

Ma, Li-Yuan; Ji, Jia-Liang; Xu, Zong-Wei; Zhu, Zuo-Nong

2018-03-01

We study a nonintegrable discrete nonlinear Schrödinger (dNLS) equation with the term of nonlinear nearest-neighbor interaction occurred in nonlinear optical waveguide arrays. By using discrete Fourier transformation, we obtain numerical approximations of stationary and travelling solitary wave solutions of the nonintegrable dNLS equation. The analysis of stability of stationary solitary waves is performed. It is shown that the nonlinear nearest-neighbor interaction term has great influence on the form of solitary wave. The shape of solitary wave is important in the electric field propagating. If we neglect the nonlinear nearest-neighbor interaction term, much important information in the electric field propagating may be missed. Our numerical simulation also demonstrates the difference of chaos phenomenon between the nonintegrable dNLS equation with nonlinear nearest-neighbor interaction and another nonintegrable dNLS equation without the term. Project supported by the National Natural Science Foundation of China (Grant Nos. 11671255 and 11701510), the Ministry of Economy and Competitiveness of Spain (Grant No. MTM2016-80276-P (AEI/FEDER, EU)), and the China Postdoctoral Science Foundation (Grant No. 2017M621964).

Linear wave propagation in a hot axisymmetric toroidal plasma

International Nuclear Information System (INIS)

Jaun, A.

1995-03-01

Kinetic effects on the propagation of the Alfven wave are studied for the first time in a toroidal plasma relevant for experiments. This requires the resolution of a set of coupled partial differential equations whose coefficients depend locally on the plasma parameters. For this purpose, a numerical wave propagation code called PENN has been developed using either a bilinear or a bicubic Hermite finite element discretization. It solves Maxwell's equations in toroidal geometry, with a dielectric tensor operator that takes into account the linear response of the plasma. Two different models have been implemented and can be used comparatively to describe the same physical case: the first treats the plasma as resistive fluids and gives results which are in good agreement with toroidal fluid codes. The second is a kinetic model and takes into account the finite size of the Larmor radii; it has successfully been tested against a kinetic plasma model in cylindrical geometry. New results have been obtained when studying kinetic effects in toroidal geometry. Two different conversion mechanisms to the kinetic Alfven wave have been described: one occurs at toroidally coupled resonant surfaces and is the kinetic counterpart of the fluid models' resonance absorption. The other has no such correspondence and results directly from the toroidal coupling between the kinetic Alfven wave and the global wavefield. An analysis of a heating scenario suggests that it might be difficult to heat a plasma with Alfven waves up to temperatures that are relevant for a tokamak reactor. Kinetic effects are studied for three types of global Alfven modes (GAE, TAE, BAE) and a new class of kinetic eigenmodes is described which appear inside the fluid gap: it could be related to recent observations in the JET (Joint European Torus) tokamak. (author) 56 figs., 6 tabs., 58 refs

Linear wave propagation in a hot axisymmetric toroidal plasma

Energy Technology Data Exchange (ETDEWEB)

Jaun, A [Ecole Polytechnique Federale, Lausanne (Switzerland). Centre de Recherche en Physique des Plasma (CRPP)

1995-03-01

Kinetic effects on the propagation of the Alfven wave are studied for the first time in a toroidal plasma relevant for experiments. This requires the resolution of a set of coupled partial differential equations whose coefficients depend locally on the plasma parameters. For this purpose, a numerical wave propagation code called PENN has been developed using either a bilinear or a bicubic Hermite finite element discretization. It solves Maxwell`s equations in toroidal geometry, with a dielectric tensor operator that takes into account the linear response of the plasma. Two different models have been implemented and can be used comparatively to describe the same physical case: the first treats the plasma as resistive fluids and gives results which are in good agreement with toroidal fluid codes. The second is a kinetic model and takes into account the finite size of the Larmor radii; it has successfully been tested against a kinetic plasma model in cylindrical geometry. New results have been obtained when studying kinetic effects in toroidal geometry. Two different conversion mechanisms to the kinetic Alfven wave have been described: one occurs at toroidally coupled resonant surfaces and is the kinetic counterpart of the fluid models` resonance absorption. The other has no such correspondence and results directly from the toroidal coupling between the kinetic Alfven wave and the global wavefield. An analysis of a heating scenario suggests that it might be difficult to heat a plasma with Alfven waves up to temperatures that are relevant for a tokamak reactor. Kinetic effects are studied for three types of global Alfven modes (GAE, TAE, BAE) and a new class of kinetic eigenmodes is described which appear inside the fluid gap: it could be related to recent observations in the JET (Joint European Torus) tokamak. (author) 56 figs., 6 tabs., 58 refs.

Experimental and modeling analysis of fast ionization wave discharge propagation in a rectangular geometry

International Nuclear Information System (INIS)

Takashima, Keisuke; Adamovich, Igor V.; Xiong Zhongmin; Kushner, Mark J.; Starikovskaia, Svetlana; Czarnetzki, Uwe; Luggenhoelscher, Dirk

2011-01-01

Fast ionization wave (FIW), nanosecond pulse discharge propagation in nitrogen and helium in a rectangular geometry channel/waveguide is studied experimentally using calibrated capacitive probe measurements. The repetitive nanosecond pulse discharge in the channel was generated using a custom designed pulsed plasma generator (peak voltage 10-40 kV, pulse duration 30-100 ns, and voltage rise time ∼1 kV/ns), generating a sequence of alternating polarity high-voltage pulses at a pulse repetition rate of 20 Hz. Both negative polarity and positive polarity ionization waves have been studied. Ionization wave speed, as well as time-resolved potential distributions and axial electric field distributions in the propagating discharge are inferred from the capacitive probe data. ICCD images show that at the present conditions the FIW discharge in helium is diffuse and volume-filling, while in nitrogen the discharge propagates along the walls of the channel. FIW discharge propagation has been analyzed numerically using quasi-one-dimensional and two-dimensional kinetic models in a hydrodynamic (drift-diffusion), local ionization approximation. The wave speed and the electric field distribution in the wave front predicted by the model are in good agreement with the experimental results. A self-similar analytic solution of the fast ionization wave propagation equations has also been obtained. The analytic model of the FIW discharge predicts key ionization wave parameters, such as wave speed, peak electric field in the front, potential difference across the wave, and electron density as functions of the waveform on the high voltage electrode, in good agreement with the numerical calculations and the experimental results.

Nonlinear wave equation in frequency domain: accurate modeling of ultrafast interaction in anisotropic nonlinear media

DEFF Research Database (Denmark)

Guo, Hairun; Zeng, Xianglong; Zhou, Binbin

2013-01-01

We interpret the purely spectral forward Maxwell equation with up to third-order induced polarizations for pulse propagation and interactions in quadratic nonlinear crystals. The interpreted equation, also named the nonlinear wave equation in the frequency domain, includes quadratic and cubic...... nonlinearities, delayed Raman effects, and anisotropic nonlinearities. The full potential of this wave equation is demonstrated by investigating simulations of solitons generated in the process of ultrafast cascaded second-harmonic generation. We show that a balance in the soliton delay can be achieved due...

Numerical modeling of the pulse wave propagation in large blood vessels based on liquid and wall interaction

International Nuclear Information System (INIS)

Rup, K; Dróżdż, A

2014-01-01

The purpose of this article is to develop a non-linear, one-dimensional model of pulse wave propagation in the arterial cardiovascular system. The model includes partial differential equations resulting from the balance of mass and momentum for the fluid-filled area and the balance equation for the area of the wall and vessels. The considered mathematical model of pulse wave propagation in the thoracic aorta section takes into account the viscous dissipation of fluid energy, realistic values of parameters describing the physicochemical properties of blood and vessel wall. Boundary and initial conditions contain the appropriate information obtained from in vivo measurements. As a result of the numerical solution of the mass and momentum balance equations for the blood and the equilibrium equation for the arterial wall area, time- dependent deformation, respective velocity profiles and blood pressure were determined.

Heat wave propagation in a thin film irradiated by ultra-short laser pulses

International Nuclear Information System (INIS)

Yoo, Jae Gwon; Kim, Cheol Jung; Lim, C. H.

2004-01-01

A thermal wave solution of a hyperbolic heat conduction equation in a thin film is developed on the basis of the Green's function formalism. Numerical computations are carried out to investigate the temperature response and the propagation of the thermal wave inside a thin film due to a heat pulse generated by ultra-short laser pulses with various laser pulse durations and thickness of the film

Rigorous vector wave propagation for arbitrary flat media

Science.gov (United States)

Bos, Steven P.; Haffert, Sebastiaan Y.; Keller, Christoph U.

2017-08-01

Precise modelling of the (off-axis) point spread function (PSF) to identify geometrical and polarization aberrations is important for many optical systems. In order to characterise the PSF of the system in all Stokes parameters, an end-to-end simulation of the system has to be performed in which Maxwell's equations are rigorously solved. We present the first results of a python code that we are developing to perform multiscale end-to-end wave propagation simulations that include all relevant physics. Currently we can handle plane-parallel near- and far-field vector diffraction effects of propagating waves in homogeneous isotropic and anisotropic materials, refraction and reflection of flat parallel surfaces, interference effects in thin films and unpolarized light. We show that the code has a numerical precision on the order of 10-16 for non-absorbing isotropic and anisotropic materials. For absorbing materials the precision is on the order of 10-8. The capabilities of the code are demonstrated by simulating a converging beam reflecting from a flat aluminium mirror at normal incidence.

Surface Waves Propagating on Grounded Anisotropic Dielectric Slab

Directory of Open Access Journals (Sweden)

Zhuozhu Chen

2018-01-01

Full Text Available This paper investigates the characteristics of surface waves propagating on a grounded anisotropic dielectric slab. Distinct from the existing analyses that generally assume that the fields of surface wave uniformly distribute along the transverse direction of the infinitely large grounded slab, our method takes into account the field variations along the transverse direction of a finite-width slab. By solving Maxwell’s equations in closed-form, it is revealed that no pure transverse magnetic (TM or transverse electric (TE mode exists if the fields are non-uniformly distributed along the transverse direction of the grounded slab. Instead, two hybrid modes, namely quasi-TM and quasi-TE modes, are supported. In addition, the propagation characteristics of two hybrid modes supported by the grounded anisotropic slab are analyzed in terms of the slab thickness, slab width, as well as the relative permittivity tensor of the anisotropic slab. Furthermore, different methods are employed to compare the analyses, as well as to validate our derivations. The proposed method is very suitable for practical engineering applications.

Wave Propagation of Coupled Modes in the DNA Double Helix

International Nuclear Information System (INIS)

Tabi, Conrad B.; Mohamadou, Alidou; Kofane, Timoleon C.

2010-06-01

The dynamics of waves propagating along the DNA molecule is described by the coupled nonlinear Schroedinger equations. We consider both the single and the coupled nonlinear excitation modes, and we discuss their biological implications. Furthermore, the characteristics of the coupled mode solution are discussed and we show that such a solution can describe the local opening observed within the transcription and the replication phenomena. (author)

Computational study on full-wave inversion based on the acoustic wave-equation; Onkyoha hado hoteishiki full wave inversion no model keisan ni yoru kento

Energy Technology Data Exchange (ETDEWEB)

Watanabe, T; Sassa, K [Kyoto University, Kyoto (Japan); Uesaka, S [Kyoto University, Kyoto (Japan). Faculty of Engineering

1996-10-01

The effect of initial models on full-wave inversion (FWI) analysis based on acoustic wave-equation was studied for elastic wave tomography of underground structures. At present, travel time inversion using initial motion travel time is generally used, and inverse analysis is conducted using the concept `ray,` assuming very high wave frequency. Although this method can derive stable solutions relatively unaffected by initial model, it uses only the data of initial motion travel time. FWI calculates theoretical waveform at each receiver using all of observed waveforms as data by wave equation modeling where 2-D underground structure is calculated by difference calculus under the assumption that wave propagation is described by wave equation of P wave. Although it is a weak point that FWI is easily affected by noises in an initial model and data, it is featured by high resolution of solutions. This method offers very excellent convergence as a proper initial model is used, resulting in sufficient performance, however, it is strongly affected by initial model. 2 refs., 7 figs., 1 tab.

Lamb wave propagation in monocrystalline silicon wafers

OpenAIRE

Fromme, P.; Pizzolato, M.; Robyr, J-L; Masserey, B.

2018-01-01

Monocrystalline silicon wafers are widely used in the photovoltaic industry for solar panels with high conversion efficiency. Guided ultrasonic waves offer the potential to efficiently detect micro-cracks in the thin wafers. Previous studies of ultrasonic wave propagation in silicon focused on effects of material anisotropy on bulk ultrasonic waves, but the dependence of the wave propagation characteristics on the material anisotropy is not well understood for Lamb waves. The phase slowness a...

Acoustic Wave Propagation in Snow Based on a Biot-Type Porous Model

Science.gov (United States)

Sidler, R.

2014-12-01

Despite the fact that acoustic methods are inexpensive, robust and simple, the application of seismic waves to snow has been sparse. This might be due to the strong attenuation inherent to snow that prevents large scale seismic applications or due to the somewhat counterintuitive acoustic behavior of snow as a porous material. Such materials support a second kind of compressional wave that can be measured in fresh snow and which has a decreasing wave velocity with increasing density of snow. To investigate wave propagation in snow we construct a Biot-type porous model of snow as a function of porosity based on the assumptions that the solid frame is build of ice, the pore space is filled with a mix of air, or air and water, and empirical relationships for the tortuosity, the permeability, the bulk, and the shear modulus.We use this reduced model to investigate compressional and shear wave velocities of snow as a function of porosity and to asses the consequences of liquid water in the snowpack on acoustic wave propagation by solving Biot's differential equations with plain wave solutions. We find that the fast compressional wave velocity increases significantly with increasing density, but also that the fast compressional wave velocity might be even lower than the slow compressional wave velocity for very light snow. By using compressional and shear strength criteria and solving Biot's differential equations with a pseudo-spectral approach we evaluate snow failure due to acoustic waves in a heterogeneous snowpack, which we think is an important mechanism in triggering avalanches by explosives as well as by skiers. Finally, we developed a low cost seismic acquisition device to assess the theoretically obtained wave velocities in the field and to explore the possibility of an inexpensive tool to remotely gather snow water equivalent.

Spinor-electron wave guided modes in coupled quantum wells structures by solving the Dirac equation

International Nuclear Information System (INIS)

Linares, Jesus; Nistal, Maria C.

2009-01-01

A quantum analysis based on the Dirac equation of the propagation of spinor-electron waves in coupled quantum wells, or equivalently coupled electron waveguides, is presented. The complete optical wave equations for Spin-Up (SU) and Spin-Down (SD) spinor-electron waves in these electron guides couplers are derived from the Dirac equation. The relativistic amplitudes and dispersion equations of the spinor-electron wave-guided modes in a planar quantum coupler formed by two coupled quantum wells, or equivalently by two coupled slab electron waveguides, are exactly derived. The main outcomes related to the spinor modal structure, such as the breaking of the non-relativistic degenerate spin states, the appearance of phase shifts associated with the spin polarization and so on, are shown.

Scalar evolution equations for shear waves in incompressible solids: a simple derivation of the Z, ZK, KZK and KP equations

OpenAIRE

Destrade, Michel; Goriely, Alain; Saccomandi, Giuseppe

2011-01-01

We study the propagation of two-dimensional finite-amplitude shear waves in a nonlinear pre-strained incompressible solid, and derive several asymptotic amplitude equations in a simple, consistent, and rigorous manner. The scalar Zabolotskaya (Z) equation is shown to be the asymptotic limit of the equations of motion for all elastic generalized neo-Hookean solids (with strain energy depending only on the first principal invariant of Cauchy-Green strain). However, we show that the Z equation c...

Harmonic surface wave propagation in plasma

International Nuclear Information System (INIS)

Shivarova, A.; Stoychev, T.

1980-01-01

Second order harmonic surface waves generated by one fundamental high-frequency surface wave are investigated experimentally in gas discharge plasma. Two types of harmonic waves of equal frequency, associated with the linear dispersion relation and the synchronism conditions relatively propagate. The experimental conditions and the different space damping rates of the waves ensure the existence of different spatial regions (consecutively arranged along the plasma column) of a dominant propagation of each one of these two waves. Experimental data are obtained both for the wavenumbers and the space damping rates by relatively precise methods for wave investigations such as the methods of time-space diagrams and of phase shift measurements. The results are explained by the theoretical model for nonlinear mixing of dispersive waves. (author)

The propagation of nonlinear rayleigh waves in layered elastic half-space

International Nuclear Information System (INIS)

Ahmetolan, S.

2004-01-01

In this work, the propagation of small but finite amplitude generalized Rayleigh waves in an elastic half-space covered by a different elastic layer of uniform and finite thickness is considered. The constituent materials are assumed to be homogeneous, isotropic, compressible hyperelastic. Excluding the harmonic resonance phenomena, it is shown that the nonlinear self modulation of generalized Rayleigh waves is governed asymptotically by a nonlinear Schrodinger (NLS) equation. The stability of the solutions and the existence of solitary wave-type solutions a NLS are strongly depend on the sign of the product of the coefficients of the nonlinear and dipersion terms of the equation.Therefore the analysis continues with the examination of dependence of these coefficients on the nonlinear material parameters. Three different models have been considered which are nonlinear layer-nonlinear half space, linear layer-nonlinear half space and nonlinear layer-linear half space. The behavior of the coefficients of the NLS equation was also analyzed the limit as h(thickness of the layer) goes to zero and k(the wave number) is constant. Then conclusions are drawn about the effect of nonlinear material parameters on the wave modulation. In the numerical investigations both hypothetical and real material models are used

Wave propagation on a plasma media

International Nuclear Information System (INIS)

Torres-Silva, H.; Villarroel-Gonzalez, C.; Reggiani, N.; Sakanaka, P.H.

1995-01-01

Chiral-media and ferrite media have been studied over the last decade for many applications. Chiral-media have been examined as coating for reducing radar cross section, for antennas and arrays, for antenna radomes in waveguides and for microstrip substrate. Here, we examine a chiral-plasma medium, where the plasma part of the composite medium is non-reciprocal due to the external magnetic field, to find the general dispersion relation giving the ω against K behavior, vector phasor Helmholtz based equations are derived. We determine the modal eigenvalue properties in the chiral-plasma medium, which is doubly anisotropic. For the case of waves which propagate parallel to the magnetic field is a cold magnetized chiro-plasma. We compare our results with the typical results obtained for a cold plasma. Also we obtain the chiral-Faraday rotation which can be compared with the typical Faraday rotation for a pair of right-and left-handed circularly polarized waves. (author). 5 refs., 2 figs

Periodic solutions for one dimensional wave equation with bounded nonlinearity

Science.gov (United States)

Ji, Shuguan

2018-05-01

This paper is concerned with the periodic solutions for the one dimensional nonlinear wave equation with either constant or variable coefficients. The constant coefficient model corresponds to the classical wave equation, while the variable coefficient model arises from the forced vibrations of a nonhomogeneous string and the propagation of seismic waves in nonisotropic media. For finding the periodic solutions of variable coefficient wave equation, it is usually required that the coefficient u (x) satisfies ess infηu (x) > 0 with ηu (x) = 1/2 u″/u - 1/4 (u‧/u)2, which actually excludes the classical constant coefficient model. For the case ηu (x) = 0, it is indicated to remain an open problem by Barbu and Pavel (1997) [6]. In this work, for the periods having the form T = 2p-1/q (p , q are positive integers) and some types of boundary value conditions, we find some fundamental properties for the wave operator with either constant or variable coefficients. Based on these properties, we obtain the existence of periodic solutions when the nonlinearity is monotone and bounded. Such nonlinearity may cross multiple eigenvalues of the corresponding wave operator. In particular, we do not require the condition ess infηu (x) > 0.

A Full-wave Model for Wave Propagation and Dissipation in the Inner Magnetosphere Using the Finite Element Method

International Nuclear Information System (INIS)

Valeo, Ernest; Johnson, Jay R.; Kim, Eun-Hwa; Phillips, Cynthia

2012-01-01

A wide variety of plasma waves play an important role in the energization and loss of particles in the inner magnetosphere. Our ability to understand and model wave-particle interactions in this region requires improved knowledge of the spatial distribution and properties of these waves as well as improved understanding of how the waves depend on changes in solar wind forcing and/or geomagnetic activity. To this end, we have developed a two-dimensional, finite element code that solves the full wave equations in global magnetospheric geometry. The code describes three-dimensional wave structure including mode conversion when ULF, EMIC, and whistler waves are launched in a two-dimensional axisymmetric background plasma with general magnetic field topology. We illustrate the capabilities of the code by examining the role of plasmaspheric plumes on magnetosonic wave propagation; mode conversion at the ion-ion and Alfven resonances resulting from external, solar wind compressions; and wave structure and mode conversion of electromagnetic ion cyclotron waves launched in the equatorial magnetosphere, which propagate along the magnetic field lines toward the ionosphere. We also discuss advantages of the finite element method for resolving resonant structures, and how the model may be adapted to include nonlocal kinetic effects.

Wave propagation analysis of a size-dependent magneto-electro-elastic heterogeneous nanoplate

Science.gov (United States)

Ebrahimi, Farzad; Dabbagh, Ali; Reza Barati, Mohammad

2016-12-01

The analysis of the wave propagation behavior of a magneto-electro-elastic functionally graded (MEE-FG) nanoplate is carried out in the framework of a refined higher-order plate theory. In order to take into account the small-scale influence, the nonlocal elasticity theory of Eringen is employed. Furthermore, the material properties of the nanoplate are considered to be variable through the thickness based on the power-law form. Nonlocal governing equations of the MEE-FG nanoplate have been derived using Hamilton's principle. The results of the present study have been validated by comparing them with previous researches. An analytical solution of governing equations is performed to obtain wave frequencies, phase velocities and escape frequencies. The effect of different parameters, such as wave number, nonlocal parameter, gradient index, magnetic potential and electric voltage on the wave dispersion characteristics of MEE-FG nanoscale plates is studied in detail.

Solitary Wave Solutions of the Boussinesq Equation and Its Improved Form

Directory of Open Access Journals (Sweden)

Reza Abazari

2013-01-01

Full Text Available This paper presents the general case study of previous works on generalized Boussinesq equations, (Abazari, 2011 and (Kılıcman and Abazari, 2012, that focuses on the application of G′/G-expansion method with the aid of Maple to construct more general exact solutions for the coupled Boussinesq equations. In this work, the mentioned method is applied to construct more general exact solutions of Boussinesq equation and improved Boussinesq equation, which the French scientist Joseph Valentin Boussinesq (1842–1929 described in the 1870s model equations for the propagation of long waves on the surface of water with small amplitude. Our work is motivated by the fact that the G′/G-expansion method provides not only more general forms of solutions but also periodic, solitary waves and rational solutions. The method appears to be easier and faster by means of a symbolic computation.

Wave propagation in thermoelastic saturated porous medium

Indian Academy of Sciences (India)

the existence and propagation of four waves in the medium. Three of the waves are ... predicted infinite speed for propagation of ther- mal signals. Lord and ..... saturated reservoir rock (North-sea Sandstone) is chosen for the numerical model ...

Propagation of thermal and hydromagnetic waves in an ionizing-recombining hydrogen plasma

International Nuclear Information System (INIS)

Di Sigalotti, Leonardo G.; Sira, Eloy; Rendon, Otto; Tremola, Ciro; Mendoza-Briceno, Cesar A.

2004-01-01

The propagation of thermal and magnetohydrodynamic (MHD) waves in a heat-conducting, hydrogen plasma, threaded by an external uniform magnetic field (B) and in which photoionization and photorecombination [H + +e - H+hν(χ)] processes are progressing, is investigated here using linear analysis. The resulting dispersion equation is solved analytically for varied strength (β<<1 and ∼1) and orientation of the magnetic field, where β denotes the ratio of plasma to magnetic pressures. Application of this model to the interstellar medium shows that heat conduction governs the propagation of thermal waves only at relatively high frequencies regardless of the plasma temperature, strength, and orientation of the magnetic field. When the direction of wave propagation is held perpendicular to B (i.e., k perpendicular B), the magnetosonic phase velocity is closely Alfvenic for β<<1, while for β∼1 both the hydrostatic and magnetic pressures determine the wave velocity. As long as k parallel B, the fast (transverse) magnetosonic wave becomes an Alfven wave for all frequencies independent of the plasma temperature and field strength, while the slow (longitudinal) magnetosonic wave becomes a pure sound wave. Amplification of thermal and MHD waves always occur at low frequencies and preferentially at temperatures for which the plasma is either weakly or partially ionized. Compared to previous analysis for the same hydrogen plasma model with B=0, the presence of the magnetic field makes the functional dependence of the physical quantities span a longer range of frequencies, which becomes progressively longer as the field strength is increased

Controlling wave propagation through nonlinear engineered granular systems

Science.gov (United States)

Leonard, Andrea

We study the fundamental dynamic behavior of a special class of ordered granular systems in order to design new, structured materials with unique physical properties. The dynamic properties of granular systems are dictated by the nonlinear, Hertzian, potential in compression and zero tensile strength resulting from the discrete material structure. Engineering the underlying particle arrangement of granular systems allows for unique dynamic properties, not observed in natural, disordered granular media. While extensive studies on 1D granular crystals have suggested their usefulness for a variety of engineering applications, considerably less attention has been given to higher-dimensional systems. The extension of these studies in higher dimensions could enable the discovery of richer physical phenomena not possible in 1D, such as spatial redirection and anisotropic energy trapping. We present experiments, numerical simulation (based on a discrete particle model), and in some cases theoretical predictions for several engineered granular systems, studying the effects of particle arrangement on the highly nonlinear transient wave propagation to develop means for controlling the wave propagation pathways. The first component of this thesis studies the stress wave propagation resulting from a localized impulsive loading for three different 2D particle lattice structures: square, centered square, and hexagonal granular crystals. By varying the lattice structure, we observe a wide range of properties for the propagating stress waves: quasi-1D solitary wave propagation, fully 2D wave propagation with tunable wave front shapes, and 2D pulsed wave propagation. Additionally the effects of weak disorder, inevitably present in real granular systems, are investigated. The second half of this thesis studies the solitary wave propagation through 2D and 3D ordered networks of granular chains, reducing the effective density compared to granular crystals by selectively placing wave

Theory for stationary nonlinear wave propagation in complex magnetic geometry

International Nuclear Information System (INIS)

Watanabe, T.; Hojo, H.; Nishikawa, Kyoji.

1977-08-01

We present our recent efforts to derive a systematic calculation scheme for nonlinear wave propagation in the self-consistent plasma profile in complex magnetic-field geometry. Basic assumptions and/or approximations are i) use of the collisionless two-fluid model with an equation of state; ii) restriction to a steady state propagation and iii) existence of modified magnetic surface, modification due to Coriolis' force. We discuss four situations: i) weak-field propagation without static flow, ii) arbitrary field strength with flow in axisymmetric system, iii) weak field limit of case ii) and iv) arbitrary field strength in nonaxisymmetric torus. Except for case iii), we derive a simple variation principle, similar to that of Seligar and Whitham, by introducing appropriate coordinates. In cases i) and iii), we derive explicit results for quasilinear profile modification. (auth.)

Topology optimization of wave-propagation problems

DEFF Research Database (Denmark)

Jensen, Jakob Søndergaard; Sigmund, Ole

2006-01-01

Topology optimization is demonstrated as a useful tool for systematic design of wave-propagation problems. We illustrate the applicability of the method for optical, acoustic and elastic devices and structures.......Topology optimization is demonstrated as a useful tool for systematic design of wave-propagation problems. We illustrate the applicability of the method for optical, acoustic and elastic devices and structures....

Oblique Propagation of Electrostatic Waves in a Magnetized Electron-Positron-Ion Plasma in the Presence of Heavy Particles

Science.gov (United States)

Sarker, M.; Hossen, M. R.; Shah, M. G.; Hosen, B.; Mamun, A. A.

2018-06-01

A theoretical investigation is carried out to understand the basic features of nonlinear propagation of heavy ion-acoustic (HIA) waves subjected to an external magnetic field in an electron-positron-ion plasma that consists of cold magnetized positively charged heavy ion fluids and superthermal distributed electrons and positrons. In the nonlinear regime, the Korteweg-de Vries (K-dV) and modified K-dV (mK-dV) equations describing the propagation of HIA waves are derived. The latter admits a solitary wave solution with both positive and negative potentials (for K-dV equation) and only positive potential (for mK-dV equation) in the weak amplitude limit. It is observed that the effects of external magnetic field (obliqueness), superthermal electrons and positrons, different plasma species concentration, heavy ion dynamics, and temperature ratio significantly modify the basic features of HIA solitary waves. The application of the results in a magnetized EPI plasma, which occurs in many astrophysical objects (e.g. pulsars, cluster explosions, and active galactic nuclei) is briefly discussed.

Wave Equation Inversion of Skeletonized SurfaceWaves

KAUST Repository

Zhang, Zhendong

2015-08-19

We present a surface-wave inversion method that inverts for the S-wave velocity from the Rayleigh dispersion curve for the fundamental-mode. We call this wave equation inversion of skeletonized surface waves because the dispersion curve for the fundamental-mode Rayleigh wave is inverted using finite-difference solutions to the wave equation. The best match between the predicted and observed dispersion curves provides the optimal S-wave velocity model. Results with synthetic and field data illustrate the benefits and limitations of this method.

SSS: A code for computing one dimensional shock and detonation wave propagation

International Nuclear Information System (INIS)

Sun Chengwei

1986-01-01

The one-dimensional hydrodynamic code SSS for shock and detonation wave propagation in inert and reactive media is described. The elastic-plastic-hydrodynamic model and four burn techniques (the Arrhenius law, C-J volume, sharp shock and Forest Fire) are used. There are HOM and JWL options for the state equation of detonation products. Comparing with the SIN code published by LANL, the SSS code has several new options: laser effects, blast waves, diverging and instantaneous detonation waves with arbitrary initiation positions. Two examples are given to compare the SSS and SIN calculations with the experimental data

Wave propagation in nanostructures nonlocal continuum mechanics formulations

CERN Document Server

Gopalakrishnan, Srinivasan

2013-01-01

Wave Propagation in Nanostructures describes the fundamental and advanced concepts of waves propagating in structures that have dimensions of the order of nanometers. The book is fundamentally based on non-local elasticity theory, which includes scale effects in the continuum model. The book predominantly addresses wave behavior in carbon nanotubes and graphene structures, although the methods of analysis provided in this text are equally applicable to other nanostructures. The book takes the reader from the fundamentals of wave propagation in nanotubes to more advanced topics such as rotating nanotubes, coupled nanotubes, and nanotubes with magnetic field and surface effects. The first few chapters cover the basics of wave propagation, different modeling schemes for nanostructures and introduce non-local elasticity theories, which form the building blocks for understanding the material provided in later chapters. A number of interesting examples are provided to illustrate the important features of wave behav...

Wave-equation dispersion inversion

KAUST Repository

Li, Jing

2016-12-08

We present the theory for wave-equation inversion of dispersion curves, where the misfit function is the sum of the squared differences between the wavenumbers along the predicted and observed dispersion curves. The dispersion curves are obtained from Rayleigh waves recorded by vertical-component geophones. Similar to wave-equation traveltime tomography, the complicated surface wave arrivals in traces are skeletonized as simpler data, namely the picked dispersion curves in the phase-velocity and frequency domains. Solutions to the elastic wave equation and an iterative optimization method are then used to invert these curves for 2-D or 3-D S-wave velocity models. This procedure, denoted as wave-equation dispersion inversion (WD), does not require the assumption of a layered model and is significantly less prone to the cycle-skipping problems of full waveform inversion. The synthetic and field data examples demonstrate that WD can approximately reconstruct the S-wave velocity distributions in laterally heterogeneous media if the dispersion curves can be identified and picked. The WD method is easily extended to anisotropic data and the inversion of dispersion curves associated with Love waves.

Transient difference solutions of the inhomogeneous wave equation - Simulation of the Green's function

Science.gov (United States)

Baumeister, K. J.

1983-01-01

A time-dependent finite difference formulation to the inhomogeneous wave equation is derived for plane wave propagation with harmonic noise sources. The difference equation and boundary conditions are developed along with the techniques to simulate the Dirac delta function associated with a concentrated noise source. Example calculations are presented for the Green's function and distributed noise sources. For the example considered, the desired Fourier transformed acoustic pressures are determined from the transient pressures by use of a ramping function and an integration technique, both of which eliminates the nonharmonic pressure associated with the initial transient.

Transient difference solutions of the inhomogeneous wave equation: Simulation of the Green's function

Science.gov (United States)

Baumeiste, K. J.

1983-01-01

A time-dependent finite difference formulation to the inhomogeneous wave equation is derived for plane wave propagation with harmonic noise sources. The difference equation and boundary conditions are developed along with the techniques to simulate the Dirac delta function associated with a concentrated noise source. Example calculations are presented for the Green's function and distributed noise sources. For the example considered, the desired Fourier transformed acoustic pressures are determined from the transient pressures by use of a ramping function and an integration technique, both of which eliminates the nonharmonic pressure associated with the initial transient.

Attenuation compensation in least-squares reverse time migration using the visco-acoustic wave equation

KAUST Repository

Dutta, Gaurav

2013-08-20

Attenuation leads to distortion of amplitude and phase of seismic waves propagating inside the earth. Conventional acoustic and least-squares reverse time migration do not account for this distortion which leads to defocusing of migration images in highly attenuative geological environments. To account for this distortion, we propose to use the visco-acoustic wave equation for least-squares reverse time migration. Numerical tests on synthetic data show that least-squares reverse time migration with the visco-acoustic wave equation corrects for this distortion and produces images with better balanced amplitudes compared to the conventional approach. © 2013 SEG.

Theory and experiment on electromagnetic-wave-propagation velocities in stacked superconducting tunnel structures

DEFF Research Database (Denmark)

Sakai, S.; Ustinov, A. V.; Kohlstedt, H.

1994-01-01

Characteristic velocities of the electromagnetic waves propagating in vertically stacked Josephson transmission are theoretically discussed. An equation for solving n velocities of the waves in an n Josephson-junction stack is derived. The solutions of two- and threefold stacks are especially...... focused on. Furthermore, under the assumption that all parameters of the layers are equal, analytic solutions for a generic N-fold stack are presented. The velocities of the waves in two- and three-junction stacks by Nb-Al-AlOx-Nb systems are experimentally obtained by measuring the cavity resonance...

Magnetoelastic shear wave propagation in pre-stressed anisotropic media under gravity

Science.gov (United States)

Kumari, Nirmala; Chattopadhyay, Amares; Singh, Abhishek K.; Sahu, Sanjeev A.

2017-03-01

The present study investigates the propagation of shear wave (horizontally polarized) in two initially stressed heterogeneous anisotropic (magnetoelastic transversely isotropic) layers in the crust overlying a transversely isotropic gravitating semi-infinite medium. Heterogeneities in both the anisotropic layers are caused due to exponential variation (case-I) and linear variation (case-II) in the elastic constants with respect to the space variable pointing positively downwards. The dispersion relations have been established in closed form using Whittaker's asymptotic expansion and were found to be in the well-agreement to the classical Love wave equations. The substantial effects of magnetoelastic coupling parameters, heterogeneity parameters, horizontal compressive initial stresses, Biot's gravity parameter, and wave number on the phase velocity of shear waves have been computed and depicted by means of a graph. As a special case, dispersion equations have been deduced when the two layers and half-space are isotropic and homogeneous. The comparative study for both cases of heterogeneity of the layers has been performed and also depicted by means of graphical illustrations.

Heat-flow equation motivated by the ideal-gas shock wave.

Science.gov (United States)

Holian, Brad Lee; Mareschal, Michel

2010-08-01

We present an equation for the heat-flux vector that goes beyond Fourier's Law of heat conduction, in order to model shockwave propagation in gases. Our approach is motivated by the observation of a disequilibrium among the three components of temperature, namely, the difference between the temperature component in the direction of a planar shock wave, versus those in the transverse directions. This difference is most prominent near the shock front. We test our heat-flow equation for the case of strong shock waves in the ideal gas, which has been studied in the past and compared to Navier-Stokes solutions. The new heat-flow treatment improves the agreement with nonequilibrium molecular-dynamics simulations of hard spheres under strong shockwave conditions.

FDTD Simulation on Terahertz Waves Propagation Through a Dusty Plasma

Science.gov (United States)

Wang, Maoyan; Zhang, Meng; Li, Guiping; Jiang, Baojun; Zhang, Xiaochuan; Xu, Jun

2016-08-01

The frequency dependent permittivity for dusty plasmas is provided by introducing the charging response factor and charge relaxation rate of airborne particles. The field equations that describe the characteristics of Terahertz (THz) waves propagation in a dusty plasma sheath are derived and discretized on the basis of the auxiliary differential equation (ADE) in the finite difference time domain (FDTD) method. Compared with numerical solutions in reference, the accuracy for the ADE FDTD method is validated. The reflection property of the metal Aluminum interlayer of the sheath at THz frequencies is discussed. The effects of the thickness, effective collision frequency, airborne particle density, and charge relaxation rate of airborne particles on the electromagnetic properties of Terahertz waves through a dusty plasma slab are investigated. Finally, some potential applications for Terahertz waves in information and communication are analyzed. supported by National Natural Science Foundation of China (Nos. 41104097, 11504252, 61201007, 41304119), the Fundamental Research Funds for the Central Universities (Nos. ZYGX2015J039, ZYGX2015J041), and the Specialized Research Fund for the Doctoral Program of Higher Education of China (No. 20120185120012)

Two-dimensional wave propagation in layered periodic media

KAUST Repository

Quezada de Luna, Manuel

2014-09-16

We study two-dimensional wave propagation in materials whose properties vary periodically in one direction only. High order homogenization is carried out to derive a dispersive effective medium approximation. One-dimensional materials with constant impedance exhibit no effective dispersion. We show that a new kind of effective dispersion may arise in two dimensions, even in materials with constant impedance. This dispersion is a macroscopic effect of microscopic diffraction caused by spatial variation in the sound speed. We analyze this dispersive effect by using highorder homogenization to derive an anisotropic, dispersive effective medium. We generalize to two dimensions a homogenization approach that has been used previously for one-dimensional problems. Pseudospectral solutions of the effective medium equations agree to high accuracy with finite volume direct numerical simulations of the variable-coeffi cient equations.

The parabolic equation method for outdoor sound propagation

DEFF Research Database (Denmark)

Arranz, Marta Galindo

The parabolic equation method is a versatile tool for outdoor sound propagation. The present study has focused on the Cranck-Nicolson type Parabolic Equation method (CNPE). Three different applications of the CNPE method have been investigated. The first two applications study variations of the g......The parabolic equation method is a versatile tool for outdoor sound propagation. The present study has focused on the Cranck-Nicolson type Parabolic Equation method (CNPE). Three different applications of the CNPE method have been investigated. The first two applications study variations...

Paracousti-UQ: A Stochastic 3-D Acoustic Wave Propagation Algorithm.

Energy Technology Data Exchange (ETDEWEB)

Preston, Leiph [Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)

2017-09-01

Acoustic full waveform algorithms, such as Paracousti, provide deterministic solutions in complex, 3-D variable environments. In reality, environmental and source characteristics are often only known in a statistical sense. Thus, to fully characterize the expected sound levels within an environment, this uncertainty in environmental and source factors should be incorporated into the acoustic simulations. Performing Monte Carlo (MC) simulations is one method of assessing this uncertainty, but it can quickly become computationally intractable for realistic problems. An alternative method, using the technique of stochastic partial differential equations (SPDE), allows computation of the statistical properties of output signals at a fraction of the computational cost of MC. Paracousti-UQ solves the SPDE system of 3-D acoustic wave propagation equations and provides estimates of the uncertainty of the output simulated wave field (e.g., amplitudes, waveforms) based on estimated probability distributions of the input medium and source parameters. This report describes the derivation of the stochastic partial differential equations, their implementation, and comparison of Paracousti-UQ results with MC simulations using simple models.

Obliquely propagating cnoidal waves in a magnetized dusty plasma with variable dust charge

International Nuclear Information System (INIS)

Yadav, L. L.; Sayal, V. K.

2009-01-01

We have studied obliquely propagating dust-acoustic nonlinear periodic waves, namely, dust-acoustic cnoidal waves, in a magnetized dusty plasma consisting of electrons, ions, and dust grains with variable dust charge. Using reductive perturbation method and appropriate boundary conditions for nonlinear periodic waves, we have derived Korteweg-de Vries (KdV) equation for the plasma. It is found that the contribution to the dispersion due to the deviation from plasma approximation is dominant for small angles of obliqueness, while for large angles of obliqueness, the dispersion due to magnetic force becomes important. The cnoidal wave solution of the KdV equation is obtained. It is found that the frequency of the cnoidal wave depends on its amplitude. The effects of the magnetic field, the angle of obliqueness, the density of electrons, the dust-charge variation and the ion-temperature on the characteristics of the dust-acoustic cnoidal wave are also discussed. It is found that in the limiting case the cnoidal wave solution reduces to dust-acoustic soliton solution.

An FMM-FFT accelerated integral equation solver for characterizing electromagnetic wave propagation in mine tunnels and galleries loaded with conductors

KAUST Repository

Yücel, Abdulkadir C.

2014-07-01

Reliable wireless communication and tracking systems in underground mines are of paramount importance to increase miners\\' productivity while monitoring the environmental conditions and increasing the effectiveness of rescue operations. Key to the design and optimization of such systems are electromagnetic (EM) simulation tools capable of analyzing wave propagation in electromagnetically large mine tunnels and galleries loaded with conducting cables (power, telephone) and mining equipment (trolleys, rails, carts), and potentially partially obstructed by debris from a cave-in. Current tools for simulating EM propagation in mine environments leverage (multi-) modal decompositions (Emslie et. al., IEEE Trans. Antennas Propag., 23, 192-205, 1975; Sun and Akyildiz, IEEE Trans. Commun., 58, 1758-1768, 2010), ray-tracing techniques (Zhang, IEEE Tran. Vehic. Tech., 5, 1308-1314, 2003), or full wave methods. Modal approaches and ray-tracing techniques cannot accurately account for the presence of conductors, intricate details of transmitters/receivers, wall roughness, or unstructured debris from a cave-in. Classical full-wave methods do not suffer from such restrictions. However, they require prohibitively large computational resources when applied to the analysis of electromagnetically large tunnels loaded with conductors. Recently, an efficient hybrid method of moment and transmission line solver has been developed to analyze the EM wave propagation inside tunnels loaded with conductors (Brocker et. al., in Proc IEEE AP-S Symp, pp.1,2, 2012). However, the applicability of the solver is limited to the characterization of EM wave propagation at medium frequency band.

Numerical simulation of ultrasonic wave propagation in elastically anisotropic media

International Nuclear Information System (INIS)

Jacob, Victoria Cristina Cheade; Jospin, Reinaldo Jacques; Bittencourt, Marcelo de Siqueira Queiroz

2013-01-01

The ultrasonic non-destructive testing of components may encounter considerable difficulties to interpret some inspections results mainly in anisotropic crystalline structures. A numerical method for the simulation of elastic wave propagation in homogeneous elastically anisotropic media, based on the general finite element approach, is used to help this interpretation. The successful modeling of elastic field associated with NDE is based on the generation of a realistic pulsed ultrasonic wave, which is launched from a piezoelectric transducer into the material under inspection. The values of elastic constants are great interest information that provide the application of equations analytical models, until small and medium complexity problems through programs of numerical analysis as finite elements and/or boundary elements. The aim of this work is the comparison between the results of numerical solution of an ultrasonic wave, which is obtained from transient excitation pulse that can be specified by either force or displacement variation across the aperture of the transducer, and the results obtained from a experiment that was realized in an aluminum block in the IEN Ultrasonic Laboratory. The wave propagation can be simulated using all the characteristics of the material used in the experiment valuation associated to boundary conditions and from these results, the comparison can be made. (author)

Lamb wave propagation in monocrystalline silicon wafers.

Science.gov (United States)

Fromme, Paul; Pizzolato, Marco; Robyr, Jean-Luc; Masserey, Bernard

2018-01-01

Monocrystalline silicon wafers are widely used in the photovoltaic industry for solar panels with high conversion efficiency. Guided ultrasonic waves offer the potential to efficiently detect micro-cracks in the thin wafers. Previous studies of ultrasonic wave propagation in silicon focused on effects of material anisotropy on bulk ultrasonic waves, but the dependence of the wave propagation characteristics on the material anisotropy is not well understood for Lamb waves. The phase slowness and beam skewing of the two fundamental Lamb wave modes A 0 and S 0 were investigated. Experimental measurements using contact wedge transducer excitation and laser measurement were conducted. Good agreement was found between the theoretically calculated angular dependency of the phase slowness and measurements for different propagation directions relative to the crystal orientation. Significant wave skew and beam widening was observed experimentally due to the anisotropy, especially for the S 0 mode. Explicit finite element simulations were conducted to visualize and quantify the guided wave beam skew. Good agreement was found for the A 0 mode, but a systematic discrepancy was observed for the S 0 mode. These effects need to be considered for the non-destructive testing of wafers using guided waves.

Enhancing propagation characteristics of truncated localized waves in silica

KAUST Repository

Salem, Mohamed

2011-07-01

The spectral characteristics of truncated Localized Waves propagating in dispersive silica are analyzed. Numerical experiments show that the immunity of the truncated Localized Waves propagating in dispersive silica to decay and distortion is enhanced as the non-linearity of the relation between the transverse spatial spectral components and the wave vector gets stronger, in contrast to free-space propagating waves, which suffer from early decay and distortion. © 2011 IEEE.

Modeling elastic wave propagation in kidney stones with application to shock wave lithotripsy.

Science.gov (United States)

Cleveland, Robin O; Sapozhnikov, Oleg A

2005-10-01

A time-domain finite-difference solution to the equations of linear elasticity was used to model the propagation of lithotripsy waves in kidney stones. The model was used to determine the loading on the stone (principal stresses and strains and maximum shear stresses and strains) due to the impact of lithotripsy shock waves. The simulations show that the peak loading induced in kidney stones is generated by constructive interference from shear waves launched from the outer edge of the stone with other waves in the stone. Notably the shear wave induced loads were significantly larger than the loads generated by the classic Hopkinson or spall effect. For simulations where the diameter of the focal spot of the lithotripter was smaller than that of the stone the loading decreased by more than 50%. The constructive interference was also sensitive to shock rise time and it was found that the peak tensile stress reduced by 30% as rise time increased from 25 to 150 ns. These results demonstrate that shear waves likely play a critical role in stone comminution and that lithotripters with large focal widths and short rise times should be effective at generating high stresses inside kidney stones.

New travelling wave solutions of the (1 + 1-dimensional cubic nonlinear Schrodinger equation using novel (G′/G-expansion method

Directory of Open Access Journals (Sweden)

M.G. Hafez

2016-06-01

Full Text Available In this paper, the novel (G′/G-expansion method is applied to construct exact travelling wave solutions of the cubic nonlinear Schrodinger equation. This technique is straightforward and simple to use, and gives more new general solutions than the other existing methods. Various types of solitary and periodic wave solutions of this equation are derived. The obtained results may be helpful to describe the wave propagation in soliton physics, such as soliton propagation in optical fibers, modulus instability in plasma physics, etc. and provided us the firm mathematical foundation in soliton physics or any varied instances. Furthermore, three-dimensional modules plot of the solutions are also given to visualize the dynamics of the equation.

Propagation of nonlinear ion acoustic wave with generation of long-wavelength waves

International Nuclear Information System (INIS)

Ohsawa, Yukiharu; Kamimura, Tetsuo

1978-01-01

The nonlinear propagation of the wave packet of an ion acoustic wave with wavenumber k 0 asymptotically equals k sub(De) (the electron Debye wavenumber) is investigated by computer simulations. From the wave packet of the ion acoustic wave, waves with long wavelengths are observed to be produced within a few periods for the amplitude oscillation of the original wave packet. These waves are generated in the region where the original wave packet exists. Their characteristic wavelength is of the order of the length of the wave packet, and their propagation velocity is almost equal to the ion acoustic speed. The long-wavelength waves thus produced strongly affect the nonlinear evolution of the original wave packet. (auth.)

Propagation behavior of two transverse surface waves in a three-layer piezoelectric/piezomagnetic structure

Science.gov (United States)

Nie, Guoquan; Liu, Jinxi; Liu, Xianglin

2017-10-01

Propagation of transverse surface waves in a three-layer system consisting of a piezoelectric/piezomagnetic (PE/PM) bi-layer bonded on an elastic half-space is theoretically investigated in this paper. Dispersion relations and mode shapes for transverse surface waves are obtained in closed form under electrically open and shorted boundary conditions at the upper surface. Two transverse surface waves related both to Love-type wave and Bleustein-Gulyaev (B-G) type wave propagating in corresponding three-layer structure are discussed through numerically solving the derived dispersion equation. The results show that Love-type wave possesses the property of multiple modes, it can exist all of the values of wavenumber for every selected thickness ratios regardless of the electrical boundary conditions. The presence of PM interlayer makes the phase velocity of Love-type wave decrease. There exist two modes allowing the propagation of B-G type wave under electrically shorted circuit, while only one mode appears in the case of electrically open circuit. The modes of B-G type wave are combinations of partly normal dispersion and partly anomalous dispersion whether the electrically open or shorted. The existence range of mode for electrically open case is greatly related to the thickness ratios, with the thickness of PM interlayer increasing the wavenumber range for existence of B-G type wave quickly shortened. When the thickness ratio is large enough, the wavenumber range of the second mode for electrically shorted circuit is extremely narrow which can be used to remove as an undesired mode. The propagation behaviors and mode shapes of transverse surface waves can be regulated by the modification of the thickness of PM interlayer. The obtained results provide a theoretical prediction and basis for applications of PE-PM composites and acoustic wave devices.

Propagation of acoustic-gravity waves in arctic zones with elastic ice-sheets

Science.gov (United States)

Kadri, Usama; Abdolali, Ali; Kirby, James T.

2017-04-01

We present an analytical solution of the boundary value problem of propagating acoustic-gravity waves generated in the ocean by earthquakes or ice-quakes in arctic zones. At the surface, we assume elastic ice-sheets of a variable thickness, and show that the propagating acoustic-gravity modes have different mode shape than originally derived by Ref. [1] for a rigid ice-sheet settings. Computationally, we couple the ice-sheet problem with the free surface model by Ref. [2] representing shrinking ice blocks in realistic sea state, where the randomly oriented ice-sheets cause inter modal transition at the edges and multidirectional reflections. We then derive a depth-integrated equation valid for spatially slowly varying thickness of ice-sheet and water depth. Surprisingly, and unlike the free-surface setting, here it is found that the higher acoustic-gravity modes exhibit a larger contribution. These modes travel at the speed of sound in water carrying information on their source, e.g. ice-sheet motion or submarine earthquake, providing various implications for ocean monitoring and detection of quakes. In addition, we found that the propagating acoustic-gravity modes can result in orbital displacements of fluid parcels sufficiently high that may contribute to deep ocean currents and circulation, as postulated by Refs. [1, 3]. References [1] U. Kadri, 2016. Generation of Hydroacoustic Waves by an Oscillating Ice Block in Arctic Zones. Advances in Acoustics and Vibration, 2016, Article ID 8076108, 7 pages http://dx.doi.org/10.1155/2016/8076108 [2] A. Abdolali, J. T. Kirby and G. Bellotti, 2015, Depth-integrated equation for hydro-acoustic waves with bottom damping, J. Fluid Mech., 766, R1 doi:10.1017/jfm.2015.37 [3] U. Kadri, 2014. Deep ocean water transportation by acoustic?gravity waves. J. Geophys. Res. Oceans, 119, doi:10.1002/ 2014JC010234

Development of numerical methods to calculate the propagation and the absorption of the hybrid wave in tokamaks

International Nuclear Information System (INIS)

Sebelin, E.

1997-01-01

Full-wave calculations based on trial functions are carried out for solving the lower hybrid current drive problem in tokamaks. A variational method is developed and provides an efficient system to describe in a global manner both the propagation and the absorption of the electromagnetic waves in plasmas. The calculation is fully carried out in the case of circular and concentric flux surfaces. The existence and uniqueness of the solution of the wave propagation equation is mathematically proved. The first realistic simulations are performed for the high aspect ratio tokamak TRIAM-1M. It is checked that the main features of the lower-hybrid wave dynamics are well described numerically. (A.C.)

Generalization of Bateman-Hillion progressive wave and Bessel-Gauss pulse solutions of the wave equation via a separation of variables

CERN Document Server

Kiselev, A

2003-01-01

Two new families of exact solutions of the wave equation u sub x sub x + u sub y sub y + u sub z sub z - c sup - sup 2 u sub t sub t = 0 generalizing Bessel-Gauss pulses and Bateman-Hillion relatively undistorted progressive waves, respectively are presented. In each of these families new simple solutions describing localized wave propagation are found. The approach is based on a kind of separation of variables. (letter to the editor)

Influences of interfacial properties on second-harmonic generation of Lamb waves propagating in layered planar structures

International Nuclear Information System (INIS)

Deng Mingxi; Wang Ping; Lv Xiafu

2006-01-01

This paper describes influences of interfacial properties on second-harmonic generation of Lamb waves propagating in layered planar structures. The nonlinearity in the elastic wave propagation is treated as a second-order perturbation of the linear elastic response. Due to the kinematic nonlinearity and the elastic nonlinearity of materials, there are second-order bulk and surface/interface driving sources in layered planar structures through which Lamb waves propagate. These driving sources can be thought of as forcing functions of a series of double frequency lamb waves (DFLWs) in terms of the approach of modal expansion analysis for waveguide excitation. The total second-harmonic fields consist of a summation of DFLWs in the corresponding stress-free layered planar structures. The interfacial properties of layered planar structures can be described by the well-known finite interfacial stiffness technique. The normal and tangential interfacial stiffness constants can be coupled with the equation governing the expansion coefficient of each DFLW component. On the other hand, the normal and tangential interfacial stiffness constants are associated with the degree of dispersion between Lamb waves and DFLWs. Theoretical analyses and numerical simulations indicate that the efficiency of second-harmonic generation by Lamb wave propagation is closely dependent on the interfacial properties of layered structures. The potential of using the effect of second-harmonic generation by Lamb wave propagation to characterize the interfacial properties of layered structures are considered. Some experimental results are presented

Propagation and dispersion of shock waves in magnetoelastic materials

Science.gov (United States)

Crum, R. S.; Domann, J. P.; Carman, G. P.; Gupta, V.

2017-12-01

Previous studies examining the response of magnetoelastic materials to shock waves have predominantly focused on applications involving pulsed power generation, with limited attention given to the actual wave propagation characteristics. This study provides detailed magnetic and mechanical measurements of magnetoelastic shock wave propagation and dispersion. Laser generated rarefacted shock waves exceeding 3 GPa with rise times of 10 ns were introduced to samples of the magnetoelastic material Galfenol. The resulting mechanical measurements reveal the evolution of the shock into a compressive acoustic front with lateral release waves. Importantly, the wave continues to disperse even after it has decayed into an acoustic wave, due in large part to magnetoelastic coupling. The magnetic data reveal predominantly shear wave mediated magnetoelastic coupling, and were also used to noninvasively measure the wave speed. The external magnetic field controlled a 30% increase in wave propagation speed, attributed to a 70% increase in average stiffness. Finally, magnetic signals propagating along the sample over 20× faster than the mechanical wave were measured, indicating these materials can act as passive antennas that transmit information in response to mechanical stimuli.

Study on monostable and bistable reaction-diffusion equations by iteration of travelling wave maps

Science.gov (United States)

Yi, Taishan; Chen, Yuming

2017-12-01

In this paper, based on the iterative properties of travelling wave maps, we develop a new method to obtain spreading speeds and asymptotic propagation for monostable and bistable reaction-diffusion equations. Precisely, for Dirichlet problems of monostable reaction-diffusion equations on the half line, by making links between travelling wave maps and integral operators associated with the Dirichlet diffusion kernel (the latter is NOT invariant under translation), we obtain some iteration properties of the Dirichlet diffusion and some a priori estimates on nontrivial solutions of Dirichlet problems under travelling wave transformation. We then provide the asymptotic behavior of nontrivial solutions in the space-time region for Dirichlet problems. These enable us to develop a unified method to obtain results on heterogeneous steady states, travelling waves, spreading speeds, and asymptotic spreading behavior for Dirichlet problem of monostable reaction-diffusion equations on R+ as well as of monostable/bistable reaction-diffusion equations on R.

Propagation of ionization waves during ignition of fluorescent lamps

International Nuclear Information System (INIS)

Langer, R; Tidecks, R; Horn, S; Garner, R; Hilscher, A

2008-01-01

The propagation of the first ionization wave in a compact fluorescent lamp (T4 tube with standard electrodes) during ignition was investigated for various initial dc-voltages (both polarities measured against ground) and gas compositions (with and without mercury). In addition the effect of the presence of a fluorescent powder coating was studied. The propagation velocity of the initial wave was measured by an assembly of photomultipliers installed along the tube, which detected the light emitted by the wave head. The propagation was found to be faster for positive than for negative polarity. This effect is explained involving processes in the electrode region as well as in the wave head. Waves propagate faster in the presence of a fluorescent powder coating than without it and gases of lighter mass show a faster propagation than gases with higher mass

Skeletonized wave equation of surface wave dispersion inversion

KAUST Repository

Li, Jing

2016-09-06

We present the theory for wave equation inversion of dispersion curves, where the misfit function is the sum of the squared differences between the wavenumbers along the predicted and observed dispersion curves. Similar to wave-equation travel-time inversion, the complicated surface-wave arrivals in traces are skeletonized as simpler data, namely the picked dispersion curves in the (kx,ω) domain. Solutions to the elastic wave equation and an iterative optimization method are then used to invert these curves for 2D or 3D velocity models. This procedure, denoted as wave equation dispersion inversion (WD), does not require the assumption of a layered model and is less prone to the cycle skipping problems of full waveform inversion (FWI). The synthetic and field data examples demonstrate that WD can accurately reconstruct the S-wave velocity distribution in laterally heterogeneous media.

Wave propagation and scattering in random media

CERN Document Server

Ishimaru, Akira

1978-01-01

Wave Propagation and Scattering in Random Media, Volume 2, presents the fundamental formulations of wave propagation and scattering in random media in a unified and systematic manner. The topics covered in this book may be grouped into three categories: waves in random scatterers, waves in random continua, and rough surface scattering. Random scatterers are random distributions of many particles. Examples are rain, fog, smog, hail, ocean particles, red blood cells, polymers, and other particles in a state of Brownian motion. Random continua are the media whose characteristics vary randomly an

3D numerical simulation of the long range propagation of acoustical shock waves through a heterogeneous and moving medium

Energy Technology Data Exchange (ETDEWEB)

Luquet, David; Marchiano, Régis; Coulouvrat, François, E-mail: francois.coulouvrat@upmc.fr [Sorbonne Universités, UPMC Univ Paris 06, CNRS, UMR 7190, Institut Jean Le Rond d’Alembert, F-75005, Paris (France)

2015-10-28

Many situations involve the propagation of acoustical shock waves through flows. Natural sources such as lightning, volcano explosions, or meteoroid atmospheric entries, emit loud, low frequency, and impulsive sound that is influenced by atmospheric wind and turbulence. The sonic boom produced by a supersonic aircraft and explosion noises are examples of intense anthropogenic sources in the atmosphere. The Buzz-Saw-Noise produced by turbo-engine fan blades rotating at supersonic speed also propagates in a fast flow within the engine nacelle. Simulating these situations is challenging, given the 3D nature of the problem, the long range propagation distances relative to the central wavelength, the strongly nonlinear behavior of shocks associated to a wide-band spectrum, and finally the key role of the flow motion. With this in view, the so-called FLHOWARD (acronym for FLow and Heterogeneous One-Way Approximation for Resolution of Diffraction) method is presented with three-dimensional applications. A scalar nonlinear wave equation is established in the framework of atmospheric applications, assuming weak heterogeneities and a slow wind. It takes into account diffraction, absorption and relaxation properties of the atmosphere, quadratic nonlinearities including weak shock waves, heterogeneities of the medium in sound speed and density, and presence of a flow (assuming a mean stratified wind and 3D turbulent ? flow fluctuations of smaller amplitude). This equation is solved in the framework of the one-way method. A split-step technique allows the splitting of the non-linear wave equation into simpler equations, each corresponding to a physical effect. Each sub-equation is solved using an analytical method if possible, and finite-differences otherwise. Nonlinear effects are solved in the time domain, and others in the frequency domain. Homogeneous diffraction is handled by means of the angular spectrum method. Ground is assumed perfectly flat and rigid. Due to the 3D

3D numerical simulation of the long range propagation of acoustical shock waves through a heterogeneous and moving medium

International Nuclear Information System (INIS)

Luquet, David; Marchiano, Régis; Coulouvrat, François

2015-01-01

Many situations involve the propagation of acoustical shock waves through flows. Natural sources such as lightning, volcano explosions, or meteoroid atmospheric entries, emit loud, low frequency, and impulsive sound that is influenced by atmospheric wind and turbulence. The sonic boom produced by a supersonic aircraft and explosion noises are examples of intense anthropogenic sources in the atmosphere. The Buzz-Saw-Noise produced by turbo-engine fan blades rotating at supersonic speed also propagates in a fast flow within the engine nacelle. Simulating these situations is challenging, given the 3D nature of the problem, the long range propagation distances relative to the central wavelength, the strongly nonlinear behavior of shocks associated to a wide-band spectrum, and finally the key role of the flow motion. With this in view, the so-called FLHOWARD (acronym for FLow and Heterogeneous One-Way Approximation for Resolution of Diffraction) method is presented with three-dimensional applications. A scalar nonlinear wave equation is established in the framework of atmospheric applications, assuming weak heterogeneities and a slow wind. It takes into account diffraction, absorption and relaxation properties of the atmosphere, quadratic nonlinearities including weak shock waves, heterogeneities of the medium in sound speed and density, and presence of a flow (assuming a mean stratified wind and 3D turbulent ? flow fluctuations of smaller amplitude). This equation is solved in the framework of the one-way method. A split-step technique allows the splitting of the non-linear wave equation into simpler equations, each corresponding to a physical effect. Each sub-equation is solved using an analytical method if possible, and finite-differences otherwise. Nonlinear effects are solved in the time domain, and others in the frequency domain. Homogeneous diffraction is handled by means of the angular spectrum method. Ground is assumed perfectly flat and rigid. Due to the 3D

A high-order discontinuous Galerkin method for wave propagation through coupled elastic-acoustic media

International Nuclear Information System (INIS)

Wilcox, Lucas C.; Stadler, Georg; Burstedde, Carsten; Ghattas, Omar

2010-01-01

We introduce a high-order discontinuous Galerkin (dG) scheme for the numerical solution of three-dimensional (3D) wave propagation problems in coupled elastic-acoustic media. A velocity-strain formulation is used, which allows for the solution of the acoustic and elastic wave equations within the same unified framework. Careful attention is directed at the derivation of a numerical flux that preserves high-order accuracy in the presence of material discontinuities, including elastic-acoustic interfaces. Explicit expressions for the 3D upwind numerical flux, derived as an exact solution for the relevant Riemann problem, are provided. The method supports h-non-conforming meshes, which are particularly effective at allowing local adaptation of the mesh size to resolve strong contrasts in the local wavelength, as well as dynamic adaptivity to track solution features. The use of high-order elements controls numerical dispersion, enabling propagation over many wave periods. We prove consistency and stability of the proposed dG scheme. To study the numerical accuracy and convergence of the proposed method, we compare against analytical solutions for wave propagation problems with interfaces, including Rayleigh, Lamb, Scholte, and Stoneley waves as well as plane waves impinging on an elastic-acoustic interface. Spectral rates of convergence are demonstrated for these problems, which include a non-conforming mesh case. Finally, we present scalability results for a parallel implementation of the proposed high-order dG scheme for large-scale seismic wave propagation in a simplified earth model, demonstrating high parallel efficiency for strong scaling to the full size of the Jaguar Cray XT5 supercomputer.

Wave propagation of spectral energy content in a granular chain

NARCIS (Netherlands)

Shrivastava, Rohit Kumar; Luding, Stefan

2017-01-01

A mechanical wave is propagation of vibration with transfer of energy and momentum. Understanding the spectral energy characteristics of a propagating wave through disordered granular media can assist in understanding the overall properties of wave propagation through inhomogeneous materials like

Three-dimensional simulation of beam propagation and heat transfer in static gas Cs DPALs using wave optics and fluid dynamics models

Science.gov (United States)

Waichman, Karol; Barmashenko, Boris D.; Rosenwaks, Salman

2017-10-01

Analysis of beam propagation, kinetic and fluid dynamic processes in Cs diode pumped alkali lasers (DPALs), using wave optics model and gasdynamic code, is reported. The analysis is based on a three-dimensional, time-dependent computational fluid dynamics (3D CFD) model. The Navier-Stokes equations for momentum, heat and mass transfer are solved by a commercial Ansys FLUENT solver based on the finite volume discretization technique. The CFD code which solves the gas conservation equations includes effects of natural convection and temperature diffusion of the species in the DPAL mixture. The DPAL kinetic processes in the Cs/He/C2H6 gas mixture dealt with in this paper involve the three lowest energy levels of Cs, (1) 62S1/2, (2) 62P1/2 and (3) 62P3/2. The kinetic processes include absorption due to the 1->3 D2 transition followed by relaxation the 3 to 2 fine structure levels and stimulated emission due to the 2->1 D1 transition. Collisional quenching of levels 2 and 3 and spontaneous emission from these levels are also considered. The gas flow conservation equations are coupled to fast-Fourier-transform algorithm for transverse mode propagation to obtain a solution of the scalar paraxial propagation equation for the laser beam. The wave propagation equation is solved by the split-step beam propagation method where the gain and refractive index in the DPAL medium affect the wave amplitude and phase. Using the CFD and beam propagation models, the gas flow pattern and spatial distributions of the pump and laser intensities in the resonator were calculated for end-pumped Cs DPAL. The laser power, DPAL medium temperature and the laser beam quality were calculated as a function of pump power. The results of the theoretical model for laser power were compared to experimental results of Cs DPAL.

Inward propagating chemical waves in Taylor vortices.

Science.gov (United States)

Thompson, Barnaby W; Novak, Jan; Wilson, Mark C T; Britton, Melanie M; Taylor, Annette F

2010-04-01

Advection-reaction-diffusion (ARD) waves in the Belousov-Zhabotinsky reaction in steady Taylor-Couette vortices have been visualized using magnetic-resonance imaging and simulated using an adapted Oregonator model. We show how propagating wave behavior depends on the ratio of advective, chemical and diffusive time scales. In simulations, inward propagating spiral flamelets are observed at high Damköhler number (Da). At low Da, the reaction distributes itself over several vortices and then propagates inwards as contracting ring pulses--also observed experimentally.

Submillimeter wave propagation in tokamak plasmas

International Nuclear Information System (INIS)

Ma, C.H.; Hutchinson, D.P.; Staats, P.A.; Vander Sluis, K.L.; Mansfield, D.K.; Park, H.; Johnson, L.C.

1985-01-01

The propagation of submillimeter-waves (smm) in tokamak plasmas has been investigated both theoretically and experimentally to ensure successful measurements of electron density and plasma current distributions in tokamak devices. Theoretical analyses have been carried out to study the polarization of the smm waves in TFTR and ISX-B tokamaks. A multichord smm wave interferometer/polarimeter system has been employed to simultaneously measure the line electron density and poloidal field-induced Faraday rotation in the ISX-B tokamak. The experimental study on TFTR is under way. Computer codes have been developed and have been used to study the wave propagation and to reconstruct the distributions of plasma current and density from the measured data. The results are compared with other measurements

Submillimeter wave propagation in tokamak plasmas

International Nuclear Information System (INIS)

Ma, C.H.; Hutchinson, D.P.; Staats, P.A.; Vander Sluis, K.L.; Mansfield, D.K.; Park, H.; Johnson, L.C.

1986-01-01

Propagation of submillimeter waves (smm) in tokamak plasma was investigated both theoretically and experimentally to ensure successful measurements of electron density and plasma current distributions in tokamak devices. Theoretical analyses were carried out to study the polarization of the smm waves in TFTR and ISX-B tokamaks. A multichord smm wave interferometer/polarimeter system was employed to simultaneously measure the line electron density and poloidal field-induced Faraday rotation in the ISX-B tokamak. The experimental study on TFTR is under way. Computer codes were developed and have been used to study the wave propagation and to reconstruct the distributions of plasma current and density from the measured data. The results are compared with other measurements. 5 references, 2 figures

FDTD Simulation of Nonlinear Ultrasonic Pulse Propagation in ESWL Using Equations Including Lagrangian

Science.gov (United States)

Fukuhara, Keisuke; Morita, Nagayoshi

New FDTD algorithm is proposed for analyzing ultrasonic pulse propagation in the human body, the problem being connected with ESWL (Extracorporeal Shock Wave Lithotripsy). In this method, we do not use plane wave approximation but employ directly the original equations taking account of Lagrangian to derive new FDTD algorithms. This method is applied to an experimental setup and its numerical model that resemble actual treatment situation to compare sound pressure distributions obtained numerically with those obtained experimentally. It is shown that the present method gives clearly better results than the earlier method, in the viewpoint of numerical reappearance of strongly nonlinear waveform.

Wave propagation in a magnetically structured atmosphere. Pt. 2

International Nuclear Information System (INIS)

Roberts, B.

1981-01-01

Magnetic fields may introduce structure (inhomogeneity) into an otherwise uniform medium and thus change the nature of wave propagation in that medium. As an example of such structuring, wave propagation in an isolated magnetic slab is considered. It is supposed that disturbances outside the slab are laterally non-propagating. The effect of gravity is ignored. The field can support the propagation of both body and surface waves. The existence and nature of these waves depends upon the relative magnitudes of the sound speed c 0 and Alfven speed upsilonsub(A) inside the slab, and the sound speed csub(e) in the field-free environment. (orig./WL)

Propagation of mechanical waves through a stochastic medium with spherical symmetry

Science.gov (United States)

Avendaño, Carlos G.; Reyes, J. Adrián

2018-01-01

We theoretically analyze the propagation of outgoing mechanical waves through an infinite isotropic elastic medium possessing spherical symmetry whose Lamé coefficients and density are spatial random functions characterized by well-defined statistical parameters. We derive the differential equation that governs the average displacement for a system whose properties depend on the radial coordinate. We show that such an equation is an extended version of the well-known Bessel differential equation whose perturbative additional terms contain coefficients that depend directly on the squared noise intensities and the autocorrelation lengths in an exponential decay fashion. We numerically solve the second order differential equation for several values of noise intensities and autocorrelation lengths and compare the corresponding displacement profiles with that of the exact analytic solution for the case of absent inhomogeneities.

Wave propagation through an electron cyclotron resonance layer

International Nuclear Information System (INIS)

Westerhof, E.

1997-01-01

The propagation of a wave beam through an electron cyclotron resonance layer is analysed in two-dimensional slab geometry in order to assess the deviation from cold plasma propagation due to resonant, warm plasma changes in wave dispersion. For quasi-perpendicular propagation, N ' 'parallel to'' ≅ v t /c, an O-mode beam is shown to exhibit a strong wiggle in the trajectory of the centre of the beam when passing through the fundamental electron cyclotron resonance. The effects are largest for low temperatures and close to perpendicular propagation. Predictions from standard dielectric wave energy fluxes are inconsistent with the trajectory of the beam. Qualitatively identical results are obtained for the X-mode second harmonic. In contrast, the X-mode at the fundamental resonance shows significant deviations form cold plasma propagation only for strongly oblique propagation and/or high temperatures. On the basis of the obtained results a practical suggestion is made for ray tracing near electron cyclotron resonance. (Author)

Simulation of Sound Waves Using the Lattice Boltzmann Method for Fluid Flow: Benchmark Cases for Outdoor Sound Propagation.

Science.gov (United States)

Salomons, Erik M; Lohman, Walter J A; Zhou, Han

2016-01-01

Propagation of sound waves in air can be considered as a special case of fluid dynamics. Consequently, the lattice Boltzmann method (LBM) for fluid flow can be used for simulating sound propagation. In this article application of the LBM to sound propagation is illustrated for various cases: free-field propagation, propagation over porous and non-porous ground, propagation over a noise barrier, and propagation in an atmosphere with wind. LBM results are compared with solutions of the equations of acoustics. It is found that the LBM works well for sound waves, but dissipation of sound waves with the LBM is generally much larger than real dissipation of sound waves in air. To circumvent this problem it is proposed here to use the LBM for assessing the excess sound level, i.e. the difference between the sound level and the free-field sound level. The effect of dissipation on the excess sound level is much smaller than the effect on the sound level, so the LBM can be used to estimate the excess sound level for a non-dissipative atmosphere, which is a useful quantity in atmospheric acoustics. To reduce dissipation in an LBM simulation two approaches are considered: i) reduction of the kinematic viscosity and ii) reduction of the lattice spacing.

Propagation of acoustic shock waves between parallel rigid boundaries and into shadow zones

International Nuclear Information System (INIS)

Desjouy, C.; Ollivier, S.; Dragna, D.; Blanc-Benon, P.; Marsden, O.

2015-01-01

The study of acoustic shock propagation in complex environments is of great interest for urban acoustics, but also for source localization, an underlying problematic in military applications. To give a better understanding of the phenomenon taking place during the propagation of acoustic shocks, laboratory-scale experiments and numerical simulations were performed to study the propagation of weak shock waves between parallel rigid boundaries, and into shadow zones created by corners. In particular, this work focuses on the study of the local interactions taking place between incident, reflected, and diffracted waves according to the geometry in both regular or irregular – also called Von Neumann – regimes of reflection. In this latter case, an irregular reflection can lead to the formation of a Mach stem that can modify the spatial distribution of the acoustic pressure. Short duration acoustic shock waves were produced by a 20 kilovolts electric spark source and a schlieren optical method was used to visualize the incident shockfront and the reflection/diffraction patterns. Experimental results are compared to numerical simulations based on the high-order finite difference solution of the two dimensional Navier-Stokes equations

Higher-order rogue wave solutions of the three-wave resonant interaction equation via the generalized Darboux transformation

International Nuclear Information System (INIS)

Wang, Xin; Chen, Yong; Cao, Jianli

2015-01-01

In this paper, we utilize generalized Darboux transformation to study higher-order rogue wave solutions of the three-wave resonant interaction equation, which describes the propagation and mixing of waves with different frequencies in weakly nonlinear dispersive media. A general Nth-order rogue wave solution with two characteristic velocities structural parameters and 3N independent parameters under a determined plane-wave background and a specific parameter condition is derived. As an application, we show that four fundamental rogue waves with fundamental, two kinds of line and quadrilateral patterns, or six fundamental rogue waves with fundamental, triangular, two kinds of quadrilateral and circular patterns can emerge in the second-order rogue waves. Moreover, several important wave characteristics including the maximum values, the corresponding coordinate positions of the humps, and the stability problem for some special higher-order rogue wave solutions such as the fundamental and quadrilateral cases are discussed. (paper)

Simulation of laser propagation in a plasma with a frequency wave equation

International Nuclear Information System (INIS)

Desroziers, S.; Nataf, F.; Sentis, R.

2008-01-01

The aim of this work is to perform numerical simulations of the propagation of a laser in a plasma. At each time step, one has to solve a Helmholtz equation in a domain which consists in some hundreds of millions of cells. To solve this huge linear system, we use an iterative Krylov method preconditioned by a separable matrix. The corresponding linear system is solved with a block cyclic reduction method. Some enlightenments on the parallel implementation are also given. Lastly, numerical results are presented including some features concerning the scalability of the numerical method on a parallel architecture. (authors)

Waveguide propagation of electromagnetic waves in high-density ducts aligned along the geomagnetic field in the near-equatorial magnetospheric region

International Nuclear Information System (INIS)

Kaufman, R.N.

1988-01-01

Waveguide propagation of electromagnetic waves in axial symmetric ducts with increased plasma density aligned along the constant external magnetic field is considered for frequencies, being higher than low-hybrid, in the WKB approximation. In this case tunnel effects leading to captured wave damping are taken into account. Conditions for waveguide propagation and the logarithmic decrement of damping are found. Field construction is performed using the systems of axially symmetric WKB solutions of the Maxwell equations

Shock wave propagation in neutral and ionized gases

International Nuclear Information System (INIS)

Podder, N. K.; Wilson IV, R. B.; Bletzinger, P.

2008-01-01

Preliminary measurements on a recently built shock tube are presented. Planar shock waves are excited by the spark discharge of a capacitor, and launched into the neutral argon or nitrogen gas as well as its ionized glow discharge in the pressure region 1-17 Torr. For the shock wave propagation in the neutral argon at fixed capacitor charging voltage, the shock wave velocity is found to increase nonlinearly at the lower pressures, reach a maximum at an intermediate pressure, and then decrease almost linearly at the higher pressures, whereas the shock wave strength continues to increase at a nonlinear rate over the entire range of pressure. However, at fixed gas pressure the shock wave velocity increases almost monotonically as the capacitor charging voltage is increased. For the shock wave propagation in the ionized argon glow, the shock wave is found to be most influenced by the glow discharge plasma current. As the plasma current is increased, both the shock wave propagation velocity and the dispersion width are observed to increase nonlinearly

One-step leapfrog ADI-FDTD method for simulating electromagnetic wave propagation in general dispersive media.

Science.gov (United States)

Wang, Xiang-Hua; Yin, Wen-Yan; Chen, Zhi Zhang David

2013-09-09

The one-step leapfrog alternating-direction-implicit finite-difference time-domain (ADI-FDTD) method is reformulated for simulating general electrically dispersive media. It models material dispersive properties with equivalent polarization currents. These currents are then solved with the auxiliary differential equation (ADE) and then incorporated into the one-step leapfrog ADI-FDTD method. The final equations are presented in the form similar to that of the conventional FDTD method but with second-order perturbation. The adapted method is then applied to characterize (a) electromagnetic wave propagation in a rectangular waveguide loaded with a magnetized plasma slab, (b) transmission coefficient of a plane wave normally incident on a monolayer graphene sheet biased by a magnetostatic field, and (c) surface plasmon polaritons (SPPs) propagation along a monolayer graphene sheet biased by an electrostatic field. The numerical results verify the stability, accuracy and computational efficiency of the proposed one-step leapfrog ADI-FDTD algorithm in comparison with analytical results and the results obtained with the other methods.

Influence of Plasma Pressure Fluctuation on RF Wave Propagation

International Nuclear Information System (INIS)

Liu Zhiwei; Bao Weimin; Li Xiaoping; Liu Donglin; Zhou Hui

2016-01-01

Pressure fluctuations in the plasma sheath from spacecraft reentry affect radio-frequency (RF) wave propagation. The influence of these fluctuations on wave propagation and wave properties is studied using methods derived by synthesizing the compressible turbulent flow theory, plasma theory, and electromagnetic wave theory. We study these influences on wave propagation at GPS and Ka frequencies during typical reentry by adopting stratified modeling. We analyzed the variations in reflection and transmission properties induced by pressure fluctuations. Our results show that, at the GPS frequency, if the waves are not totally reflected then the pressure fluctuations can remarkably affect reflection, transmission, and absorption properties. In extreme situations, the fluctuations can even cause blackout. At the Ka frequency, the influences are obvious when the waves are not totally transmitted. The influences are more pronounced at the GPS frequency than at the Ka frequency. This suggests that the latter can mitigate blackout by reducing both the reflection and the absorption of waves, as well as the influences of plasma fluctuations on wave propagation. Given that communication links with the reentry vehicles are susceptible to plasma pressure fluctuations, the influences on link budgets should be taken into consideration. (paper)

Determination of particle size distributions from acoustic wave propagation measurements

International Nuclear Information System (INIS)

Spelt, P.D.; Norato, M.A.; Sangani, A.S.; Tavlarides, L.L.

1999-01-01

The wave equations for the interior and exterior of the particles are ensemble averaged and combined with an analysis by Allegra and Hawley [J. Acoust. Soc. Am. 51, 1545 (1972)] for the interaction of a single particle with the incident wave to determine the phase speed and attenuation of sound waves propagating through dilute slurries. The theory is shown to compare very well with the measured attenuation. The inverse problem, i.e., the problem of determining the particle size distribution given the attenuation as a function of frequency, is examined using regularization techniques that have been successful for bubbly liquids. It is shown that, unlike the bubbly liquids, the success of solving the inverse problem is limited since it depends strongly on the nature of particles and the frequency range used in inverse calculations. copyright 1999 American Institute of Physics

Smooth, cusped, and discontinuous traveling waves in the periodic fluid resonance equation

Science.gov (United States)

Kruse, Matthew Thomas

The principal motivation for this dissertation is to extend the study of small amplitude high frequency wave propagation in solutions for hyperbolic conservation laws begun by A. Majda and R. Rosales in 1984. It was then that Majda and Rosales obtained equations governing the leading order wave amplitudes of resonantly interacting weakly nonlinear high frequency wave trains in the compressible Euler equations. The equations were obtained through systematic application of multiple scales and result in a pair of nonlinear acoustic wave equations coupled through a convolution operator. The extended solutions satisfy a pair of inviscid Burgers' equations coupled via a spatial convolution operator. Since then, many mathematicians have used this technique to extend the time validity of solutions to systems of equations other than the Euler equations and have arrived at similar nonlinear non-local systems. This work attempts to look at some of the basic features of the linear and nonlinear coupled and decoupled non- local equations, offering some analytic solutions and numerical insight into the phenomena associated with these equations. We do so by examining a single non-local linear equation, and then a single equation coupling a Burgers' nonlinearity with a linear convolution operator. The linear case is completely solvable. Analytic solutions are provided along with numerical results showing the fundamental properties of the linear non- local equations. In the nonlinear case some analytic solutions, including steady state profiles and traveling wave solutions, are provided along with a battery of numerical simulations. Evidence indicates the existence of attractors for solutions of the single equation with a single mode kernel. Provided resonant interaction takes place, the profile of the attractor is uniquely dependent on the kernel alone. Hamiltonian equations are obtained for both the linear and nonlinear equations with the condition that the resonant kernel must

Wave propagation simulation of radio occultations based on ECMWF refractivity profiles

DEFF Research Database (Denmark)

von Benzon, Hans-Henrik; Høeg, Per

2015-01-01

This paper describes a complete radio occultation simulation environment, including realistic refractivity profiles, wave propagation modeling, instrument modeling, and bending angle retrieval. The wave propagator is used to simulate radio occultation measurements. The radio waves are propagated...... of radio occultations. The output from the wave propagator simulator is used as input to a Full Spectrum Inversion retrieval module which calculates geophysical parameters. These parameters can be compared to the ECMWF atmospheric profiles. The comparison can be used to reveal system errors and get...... a better understanding of the physics. The wave propagation simulations will in this paper also be compared to real measurements. These radio occultations have been exposed to the same atmospheric conditions as the radio occultations simulated by the wave propagator. This comparison reveals that precise...

Neutron wave reflexions in interface media with transport equation P1 approximation

International Nuclear Information System (INIS)

Oliveira Vellozo, S. de.

1977-01-01

The propagation of neutron waves in non multiplying media is investigated employing the Telegrapher's equation obtained from the P 1 approximation of the time, space and energy dependent Boltzmann equation. Solution of the problem of propagation of sinusoidally modulated source incident on one face of the medium is obtained by analysing the Fourier component of a pulsed source introduced, for the corresponding frequency. The amplitude and the phase of the flux are computed as a function of frequency in media consisting of one, two and three regions in order to study the effects of reflection at the interfaces. The results are compared with those from the Diffusion approximation obtained by neglecting the term involving the second order time derivative. (author)

Wave Equation Inversion of Skeletonized SurfaceWaves

KAUST Repository

Zhang, Zhendong; Liu, Yike; Schuster, Gerard T.

2015-01-01

We present a surface-wave inversion method that inverts for the S-wave velocity from the Rayleigh dispersion curve for the fundamental-mode. We call this wave equation inversion of skeletonized surface waves because the dispersion curve

Experimental study of the fast wave propagation in TFR

International Nuclear Information System (INIS)

1981-02-01

Several experiments (PLT, DIVA, ERASMUS, TFR) have shown that the heating mechanism of ICRF is dominated in Tokamaks by the presence of the ion-ion hybrid layer. The first experimental evidence of this effect came from propagation studies: a very strong damping was observed on magnetic probes since the hybrid layer was inside the plasma. Comparison with simple models which do not take into account boundary conditions have been undertaken. Recently a new theoretical model has been developped. Based on a plane, inhomogeneous, bounded plasma, it shows that the radial structure of the fast wave and hence the loading impedance of the launching coil depends on the position of the hybrid layer with respect to the plasma boundaries. This result is obtained by solving the wave equation, in the cold plasma approximation. We present here, a serie of experiments, performed in TFR. It confirms the validity of that model underlining thus the importance of radial eigenmodes, when the wave conversion layer is inside the plasma

Stability of Planar Rarefaction Wave to 3D Full Compressible Navier-Stokes Equations

Science.gov (United States)

Li, Lin-an; Wang, Teng; Wang, Yi

2018-05-01

We prove time-asymptotic stability toward the planar rarefaction wave for the three-dimensional full, compressible Navier-Stokes equations with the heat-conductivities in an infinite long flat nozzle domain {R × T^2} . Compared with one-dimensional case, the proof here is based on our new observations on the cancellations on the flux terms and viscous terms due to the underlying wave structures, which are crucial for overcoming the difficulties due to the wave propagation in the transverse directions x 2 and x 3 and its interactions with the planar rarefaction wave in x 1 direction.

Wave propagation near cyclotron resonance in the presence of large Larmor radius particles

International Nuclear Information System (INIS)

Cairns, R.A.; Lashmore-Davies, C.N.; Holt, H.; McDonald, D.C.

1995-02-01

Absorption of waves propagating across an inhomogeneous magnetic field is of crucial importance for cyclotron resonance heating. When the Larmor radius of the resonant particles is small compared to the wavelength, then the propagation can be described by differential equations. These have been derived by a considerable number of authors, but a comparatively simple method of obtaining them has recently been given by Cairns et al [Phys. Fluids B3, 2953 (1991)] and, for the relativistic case which is relevant to electron cyclotron heating, by McDonald et al [Phys. Plasmas 1, 842 (1994)]. In a fusion plasma there may be a significant number of hot ions for which the Larmor radius is comparable to or larger than the perpendicular wavelength. It is important to be able to calculate the effect of these ions on ion cyclotron phenomena. In this case the system is described by integro-differential equations, the structure of which is essentially determined by the fact that the response at a given position is determined by the wave amplitude over a region whose width is of the order of a Larmor radius. The equations describing this situation have been obtained by Sauter and Vaclavik [Theory of Fusion Plasmas, Editrice Compositori, Bologna (1990) p. 403] and by Brambilla [Plasma Physics and Controlled Fusion 33, 1029 (1991)]. Here we show how the simplified method referred to above can be adapted to this case and used to find various alternative forms for the equations. (author)

Wave-equation Qs Inversion of Skeletonized Surface Waves

KAUST Repository

Li, Jing

2017-02-08

We present a skeletonized inversion method that inverts surface-wave data for the Qs quality factor. Similar to the inversion of dispersion curves for the S-wave velocity model, the complicated surface-wave arrivals are skeletonized as simpler data, namely the amplitude spectra of the windowed Rayleigh-wave arrivals. The optimal Qs model is the one that minimizes the difference in the peak frequencies of the predicted and observed Rayleigh wave arrivals using a gradient-based wave-equation optimization method. Solutions to the viscoelastic wave-equation are used to compute the predicted Rayleigh-wave arrivals and the misfit gradient at every iteration. This procedure, denoted as wave-equation Qs inversion (WQs), does not require the assumption of a layered model and tends to have fast and robust convergence compared to full waveform inversion (FWI). Numerical examples with synthetic and field data demonstrate that the WQs method can accurately invert for a smoothed approximation to the subsurface Qs distribution as long as the Vs model is known with sufficient accuracy.

Skeletonized wave-equation Qs tomography using surface waves

KAUST Repository

Li, Jing

2017-08-17

We present a skeletonized inversion method that inverts surface-wave data for the Qs quality factor. Similar to the inversion of dispersion curves for the S-wave velocity model, the complicated surface-wave arrivals are skeletonized as simpler data, namely the amplitude spectra of the windowed Rayleigh-wave arrivals. The optimal Qs model is then found that minimizes the difference in the peak frequencies of the predicted and observed Rayleigh wave arrivals using a gradient-based wave-equation optimization method. Solutions to the viscoelastic wave-equation are used to compute the predicted Rayleigh-wave arrivals and the misfit gradient at every iteration. This procedure, denoted as wave-equation Qs tomography (WQs), does not require the assumption of a layered model and tends to have fast and robust convergence compared to Q full waveform inversion (Q-FWI). Numerical examples with synthetic and field data demonstrate that the WQs method can accurately invert for a smoothed approximation to the subsur-face Qs distribution as long as the Vs model is known with sufficient accuracy.

Wave-equation Qs Inversion of Skeletonized Surface Waves

KAUST Repository

Li, Jing; Dutta, Gaurav; Schuster, Gerard T.

2017-01-01

We present a skeletonized inversion method that inverts surface-wave data for the Qs quality factor. Similar to the inversion of dispersion curves for the S-wave velocity model, the complicated surface-wave arrivals are skeletonized as simpler data, namely the amplitude spectra of the windowed Rayleigh-wave arrivals. The optimal Qs model is the one that minimizes the difference in the peak frequencies of the predicted and observed Rayleigh wave arrivals using a gradient-based wave-equation optimization method. Solutions to the viscoelastic wave-equation are used to compute the predicted Rayleigh-wave arrivals and the misfit gradient at every iteration. This procedure, denoted as wave-equation Qs inversion (WQs), does not require the assumption of a layered model and tends to have fast and robust convergence compared to full waveform inversion (FWI). Numerical examples with synthetic and field data demonstrate that the WQs method can accurately invert for a smoothed approximation to the subsurface Qs distribution as long as the Vs model is known with sufficient accuracy.

On advanced variational formulation of the method of lines and its application to the wave propagation problems

CSIR Research Space (South Africa)

Shatalov, M

2012-09-01

Full Text Available are transformed into systems of ordinary differential equations with initial conditions. This reduction is obtained by means of application of particular finite difference schemes to the spatial derivatives. Many of the wave propagation problems describing...

Ion temperature effects on magnetotail Alfvén wave propagation and electron energization: ION TEMPERATURE EFFECTS ON ALFVÉN WAVES

Energy Technology Data Exchange (ETDEWEB)

Damiano, P. A. [Princeton Center for Heliophysics, Princeton Plasma Physics Laboratory, Princeton University, Princeton New Jersey USA; Johnson, J. R. [Princeton Center for Heliophysics, Princeton Plasma Physics Laboratory, Princeton University, Princeton New Jersey USA; Chaston, C. C. [Space Sciences Laboratory, University of California, Berkeley California USA; School of Physics, University of Sydney, Sydney New South Wales Australia

2015-07-01

A new 2-D self-consistent hybrid gyrofluid-kinetic electron model in dipolar coordinates is presented and used to simulate dispersive-scale Alfvén wave pulse propagation from the equator to the ionosphere along an L = 10 magnetic field line. The model is an extension of the hybrid MHD-kinetic electron model that incorporates ion Larmor radius corrections via the kinetic fluid model of Cheng and Johnson (1999). It is found that consideration of a realistic ion to electron temperature ratio decreases the propagation time of the wave from the plasma sheet to the ionosphere by several seconds relative to a ρi=0 case (which also implies shorter timing for a substorm onset signal) and leads to significant dispersion of wave energy perpendicular to the ambient magnetic field. Additionally, ion temperature effects reduce the parallel current and electron energization all along the field line for the same magnitude perpendicular electric field perturbation.

Skeletonized wave equation of surface wave dispersion inversion

KAUST Repository

Li, Jing; Schuster, Gerard T.

2016-01-01

We present the theory for wave equation inversion of dispersion curves, where the misfit function is the sum of the squared differences between the wavenumbers along the predicted and observed dispersion curves. Similar to wave-equation travel

Linear superposition solutions to nonlinear wave equations

International Nuclear Information System (INIS)

Liu Yu

2012-01-01

The solutions to a linear wave equation can satisfy the principle of superposition, i.e., the linear superposition of two or more known solutions is still a solution of the linear wave equation. We show in this article that many nonlinear wave equations possess exact traveling wave solutions involving hyperbolic, triangle, and exponential functions, and the suitable linear combinations of these known solutions can also constitute linear superposition solutions to some nonlinear wave equations with special structural characteristics. The linear superposition solutions to the generalized KdV equation K(2,2,1), the Oliver water wave equation, and the k(n, n) equation are given. The structure characteristic of the nonlinear wave equations having linear superposition solutions is analyzed, and the reason why the solutions with the forms of hyperbolic, triangle, and exponential functions can form the linear superposition solutions is also discussed

Variation principle for nonlinear wave propagation

International Nuclear Information System (INIS)

Watanabe, T.; Lee, Y.C.; Nishikawa, Kyoji; Hojo, H.; Yoshida, Y.

1976-01-01

Variation principle is derived which determines stationary nonlinear propagation of electrostatic waves in the self-consistent density profile. Example is given for lower-hybrid waves and the relation to the variation principle for the Lagrangian density of electromagnetic fluids is discussed

Generalized Multiscale Finite-Element Method (GMsFEM) for elastic wave propagation in heterogeneous, anisotropic media

Energy Technology Data Exchange (ETDEWEB)

Gao, Kai, E-mail: kaigao87@gmail.com [Department of Geology and Geophysics, Texas A& M University, College Station, TX 77843 (United States); Fu, Shubin, E-mail: shubinfu89@gmail.com [Department of Mathematics, Texas A& M University, College Station, TX 77843 (United States); Gibson, Richard L., E-mail: gibson@tamu.edu [Department of Geology and Geophysics, Texas A& M University, College Station, TX 77843 (United States); Chung, Eric T., E-mail: tschung@math.cuhk.edu.hk [Department of Mathematics, The Chinese University of Hong Kong, Shatin, NT (Hong Kong); Efendiev, Yalchin, E-mail: efendiev@math.tamu.edu [Department of Mathematics, Texas A& M University, College Station, TX 77843 (United States); Numerical Porous Media SRI Center (NumPor), King Abdullah University of Science and Technology, Thuwal (Saudi Arabia)

2015-08-15

It is important to develop fast yet accurate numerical methods for seismic wave propagation to characterize complex geological structures and oil and gas reservoirs. However, the computational cost of conventional numerical modeling methods, such as finite-difference method and finite-element method, becomes prohibitively expensive when applied to very large models. We propose a Generalized Multiscale Finite-Element Method (GMsFEM) for elastic wave propagation in heterogeneous, anisotropic media, where we construct basis functions from multiple local problems for both the boundaries and interior of a coarse node support or coarse element. The application of multiscale basis functions can capture the fine scale medium property variations, and allows us to greatly reduce the degrees of freedom that are required to implement the modeling compared with conventional finite-element method for wave equation, while restricting the error to low values. We formulate the continuous Galerkin and discontinuous Galerkin formulation of the multiscale method, both of which have pros and cons. Applications of the multiscale method to three heterogeneous models show that our multiscale method can effectively model the elastic wave propagation in anisotropic media with a significant reduction in the degrees of freedom in the modeling system.

Generalized Multiscale Finite-Element Method (GMsFEM) for elastic wave propagation in heterogeneous, anisotropic media

International Nuclear Information System (INIS)

Gao, Kai; Fu, Shubin; Gibson, Richard L.; Chung, Eric T.; Efendiev, Yalchin

2015-01-01

It is important to develop fast yet accurate numerical methods for seismic wave propagation to characterize complex geological structures and oil and gas reservoirs. However, the computational cost of conventional numerical modeling methods, such as finite-difference method and finite-element method, becomes prohibitively expensive when applied to very large models. We propose a Generalized Multiscale Finite-Element Method (GMsFEM) for elastic wave propagation in heterogeneous, anisotropic media, where we construct basis functions from multiple local problems for both the boundaries and interior of a coarse node support or coarse element. The application of multiscale basis functions can capture the fine scale medium property variations, and allows us to greatly reduce the degrees of freedom that are required to implement the modeling compared with conventional finite-element method for wave equation, while restricting the error to low values. We formulate the continuous Galerkin and discontinuous Galerkin formulation of the multiscale method, both of which have pros and cons. Applications of the multiscale method to three heterogeneous models show that our multiscale method can effectively model the elastic wave propagation in anisotropic media with a significant reduction in the degrees of freedom in the modeling system

Generalized Multiscale Finite-Element Method (GMsFEM) for elastic wave propagation in heterogeneous, anisotropic media

KAUST Repository

Gao, Kai

2015-04-14

It is important to develop fast yet accurate numerical methods for seismic wave propagation to characterize complex geological structures and oil and gas reservoirs. However, the computational cost of conventional numerical modeling methods, such as finite-difference method and finite-element method, becomes prohibitively expensive when applied to very large models. We propose a Generalized Multiscale Generalized Multiscale Finite-Element Method (GMsFEM) for elastic wave propagation in heterogeneous, anisotropic media, where we construct basis functions from multiple local problems for both boundaries and the interior of a coarse node support or coarse element. The application of multiscale basis functions can capture the fine scale medium property variations, and allows us to greatly reduce the degrees of freedom that are required to implement the modeling compared with conventional finite-element method for wave equation, while restricting the error to low values. We formulate the continuous Galerkin and discontinuous Galerkin formulation of the multiscale method, both of which have pros and cons. Applications of the multiscale method to three heterogeneous models show that our multiscale method can effectively model the elastic wave propagation in anisotropic media with a significant reduction in the degrees of freedom in the modeling system.

Propagation law of impact elastic wave based on specific materials

Directory of Open Access Journals (Sweden)

Chunmin CHEN

2017-02-01

Full Text Available In order to explore the propagation law of the impact elastic wave on the platform, the experimental platform is built by using the specific isotropic materials and anisotropic materials. The glass cloth epoxy laminated plate is used for anisotropic material, and an organic glass plate is used for isotropic material. The PVDF sensors adhered on the specific materials are utilized to collect data, and the elastic wave propagation law of different thick plates and laminated plates under impact conditions is analyzed. The Experimental results show that in anisotropic material, transverse wave propagation speed along the fiber arrangement direction is the fastest, while longitudinal wave propagation speed is the slowest. The longitudinal wave propagation speed in anisotropic laminates is much slower than that in the laminated thick plates. In the test channel arranged along a particular angle away from the central region of the material, transverse wave propagation speed is larger. Based on the experimental results, this paper proposes a material combination mode which is advantageous to elastic wave propagation and diffusion in shock-isolating materials. It is proposed to design a composite material with high acoustic velocity by adding regularly arranged fibrous materials. The overall design of the barrier material is a layered structure and a certain number of 90°zigzag structure.

Integral propagator solvers for Vlasov-Fokker-Planck equations

International Nuclear Information System (INIS)

Donoso, J M; Rio, E del

2007-01-01

We briefly discuss the use of short-time integral propagators on solving the so-called Vlasov-Fokker-Planck equation for the dynamics of a distribution function. For this equation, the diffusion tensor is singular and the usual Gaussian representation of the short-time propagator is no longer valid. However, we prove that the path-integral approach on solving the equation is, in fact, reliable by means of our generalized propagator, which is obtained through the construction of an auxiliary solvable Fokker-Planck equation. The new representation of the grid-free advancing scheme describes the inherent cross- and self-diffusion processes, in both velocity and configuration spaces, in a natural manner, although these processes are not explicitly depicted in the differential equation. We also show that some splitting methods, as well as some finite-difference schemes, could fail in describing the aforementioned diffusion processes, governed in the whole phase space only by the velocity diffusion tensor. The short-time transition probability offers a stable and robust numerical algorithm that preserves the distribution positiveness and its norm, ensuring the smoothness of the evolving solution at any time step. (fast track communication)

Study on the electromagnetic waves propagation characteristics in partially ionized plasma slabs

Directory of Open Access Journals (Sweden)

Zhi-Bin Wang

2016-05-01

Full Text Available Propagation characteristics of electromagnetic (EM waves in partially ionized plasma slabs are studied in this paper. Such features are significant to applications in plasma antennas, blackout of re-entry flying vehicles, wave energy injection to plasmas, and etc. We in this paper developed a theoretical model of EM wave propagation perpendicular to a plasma slab with a one-dimensional density inhomogeneity along propagation direction to investigate essential characteristics of EM wave propagation in nonuniform plasmas. Particularly, the EM wave propagation in sub-wavelength plasma slabs, where the geometric optics approximation fails, is studied and in comparison with thicker slabs where the geometric optics approximation applies. The influences of both plasma and collisional frequencies, as well as the width of the plasma slab, on the EM wave propagation characteristics are discussed. The results can help the further understanding of propagation behaviours of EM waves in nonuniform plasma, and applications of the interactions between EM waves and plasmas.

The relativistic electron wave equation

International Nuclear Information System (INIS)

Dirac, P.A.M.

1977-08-01

The paper was presented at the European Conference on Particle Physics held in Budapest between the 4th and 9th July of 1977. A short review is given on the birth of the relativistic electron wave equation. After Schroedinger has shown the equivalence of his wave mechanics and the matrix mechanics of Heisenberg, a general transformation theory was developed by the author. This theory required a relativistic wave equation linear in delta/delta t. As the Klein--Gordon equation available at this time did not satisfy this condition the development of a new equation became necessary. The equation which was found gave the value of the electron spin and magnetic moment automatically. (D.P.)

Enabling real-time ultrasound imaging of soft tissue mechanical properties by simplification of the shear wave motion equation.

Science.gov (United States)

Engel, Aaron J; Bashford, Gregory R

2015-08-01

Ultrasound based shear wave elastography (SWE) is a technique used for non-invasive characterization and imaging of soft tissue mechanical properties. Robust estimation of shear wave propagation speed is essential for imaging of soft tissue mechanical properties. In this study we propose to estimate shear wave speed by inversion of the first-order wave equation following directional filtering. This approach relies on estimation of first-order derivatives which allows for accurate estimations using smaller smoothing filters than when estimating second-order derivatives. The performance was compared to three current methods used to estimate shear wave propagation speed: direct inversion of the wave equation (DIWE), time-to-peak (TTP) and cross-correlation (CC). The shear wave speed of three homogeneous phantoms of different elastic moduli (gelatin by weight of 5%, 7%, and 9%) were measured with each method. The proposed method was shown to produce shear speed estimates comparable to the conventional methods (standard deviation of measurements being 0.13 m/s, 0.05 m/s, and 0.12 m/s), but with simpler processing and usually less time (by a factor of 1, 13, and 20 for DIWE, CC, and TTP respectively). The proposed method was able to produce a 2-D speed estimate from a single direction of wave propagation in about four seconds using an off-the-shelf PC, showing the feasibility of performing real-time or near real-time elasticity imaging with dedicated hardware.

State equations and stability of shock wave fronts in homogeneous and heterogeneous metallic medium

International Nuclear Information System (INIS)

Romain, Jean-Pierre

1977-01-01

This research thesis in physical sciences reports a theoretical and experimental study of some mechanical and thermodynamic aspects related to a shock wave propagation in homogeneous and heterogeneous metallic media: state equations, stability and instability of shock wave fronts. In the first part, the author reports the study of the Grueneisen coefficient for some metallic elements with known static and dynamic compression properties. The second part reports the experimental investigation of dynamic compressibility of some materials (lamellar Al-Cu compounds). The front shock wave propagation has been visualised, and experimental Hugoniot curves are compared with those deduced from a developed numeric model and other models. The bismuth Hugoniot curve is also determined, and the author compares the existence and nature of phase transitions obtained by static and dynamic compression

Obliquely propagating dust-density waves

International Nuclear Information System (INIS)

Piel, A.; Arp, O.; Klindworth, M.; Melzer, A.

2008-01-01

Self-excited dust-density waves are experimentally studied in a dusty plasma under microgravity. Two types of waves are observed: a mode inside the dust volume propagating in the direction of the ion flow and another mode propagating obliquely at the boundary between the dusty plasma and the space charge sheath. The dominance of oblique modes can be described in the frame of a fluid model. It is shown that the results fom the fluid model agree remarkably well with a kinetic electrostatic model of Rosenberg [J. Vac. Sci. Technol. A 14, 631 (1996)]. In the experiment, the instability is quenched by increasing the gas pressure or decreasing the dust density. The critical pressure and dust density are well described by the models

On the quantum inverse problem for a new type of nonlinear Schroedinger equation for Alfven waves in plasma

International Nuclear Information System (INIS)

Sen, S.; Roy Chowdhury, A.

1989-06-01

The nonlinear Alfven waves are governed by the Vector Derivative nonlinear Schroedinger (VDNLS) equation, which for parallel or quasi parallel propagation reduces to the Derivative Nonlinear Schroedinger (DNLS) equation for the circularly polarized waves. We have formulated the Quantum Inverse problem for a new type of Nonlinear Schroedinger Equation which has many properties similar to the usual NLS problem but the structure of classical and quantum R matrix are distinctly different. The commutation rules of the scattering data are obtained and the Algebraic Bethe Ansatz is formulated to derive the eigenvalue equation for the energy of the excited states. 10 refs

Test of a new heat-flow equation for dense-fluid shock waves.

Science.gov (United States)

Holian, Brad Lee; Mareschal, Michel; Ravelo, Ramon

2010-09-21

Using a recently proposed equation for the heat-flux vector that goes beyond Fourier's Law of heat conduction, we model shockwave propagation in the dense Lennard-Jones fluid. Disequilibrium among the three components of temperature, namely, the difference between the kinetic temperature in the direction of a planar shock wave and those in the transverse directions, particularly in the region near the shock front, gives rise to a new transport (equilibration) mechanism not seen in usual one-dimensional heat-flow situations. The modification of the heat-flow equation was tested earlier for the case of strong shock waves in the ideal gas, which had been studied in the past and compared to Navier-Stokes-Fourier solutions. Now, the Lennard-Jones fluid, whose equation of state and transport properties have been determined from independent calculations, allows us to study the case where potential, as well as kinetic contributions are important. The new heat-flow treatment improves the agreement with nonequilibrium molecular-dynamics simulations under strong shock wave conditions, compared to Navier-Stokes.

An FMM-FFT Accelerated SIE Simulator for Analyzing EM Wave Propagation in Mine Environments Loaded with Conductors

KAUST Repository

Yucel, Abdulkadir C.

2018-02-05

A fast and memory efficient 3D full wave simulator for analyzing electromagnetic (EM) wave propagation in electrically large and realistic mine tunnels/galleries loaded with conductors is proposed. The simulator relies on Muller and combined field surface integral equations (SIEs) to account for scattering from mine walls and conductors, respectively. During the iterative solution of the system of SIEs, the simulator uses a fast multipole method - fast Fourier transform (FMM-FFT) scheme to reduce CPU and memory requirements. The memory requirement is further reduced by compressing large data structures via singular value and Tucker decompositions. The efficiency, accuracy, and real-world applicability of the simulator are demonstrated through characterization of EM wave propagation in electrically large mine tunnels/galleries loaded with conducting cables and mine carts.

Slow Wave Propagation and Sheath Interaction for ICRF Waves in the Tokamak SOL

International Nuclear Information System (INIS)

Myra, J. R.; D'Ippolito, D. A.

2009-01-01

In previous work we studied the propagation of slow-wave resonance cones launched parasitically by a fast-wave antenna into a tenuous magnetized plasma. Here we extend the previous calculation to ''dense'' scrape-off-layer (SOL) plasmas where the usual slow wave is evanescent. Using the sheath boundary condition, it is shown that for sufficiently close limiters, the slow wave couples to a sheath plasma wave and is no longer evanescent, but radially propagating. A self-consistent calculation of the rf-sheath width yields the resulting sheath voltage in terms of the amplitude of the launched SW, plasma parameters and connection length.

A Simple FDTD Algorithm for Simulating EM-Wave Propagation in General Dispersive Anisotropic Material

KAUST Repository

Al-Jabr, Ahmad Ali; Alsunaidi, Mohammad A.; Ng, Tien Khee; Ooi, Boon S.

2013-01-01

In this paper, an finite-difference time-domain (FDTD) algorithm for simulating propagation of EM waves in anisotropic material is presented. The algorithm is based on the auxiliary differential equation and the general polarization formulation. In anisotropic materials, electric fields are coupled and elements in the permittivity tensor are, in general, multiterm dispersive. The presented algorithm resolves the field coupling using a formulation based on electric polarizations. It also offers a simple procedure for the treatment of multiterm dispersion in the FDTD scheme. The algorithm is tested by simulating wave propagation in 1-D magnetized plasma showing excellent agreement with analytical solutions. Extension of the algorithm to multidimensional structures is straightforward. The presented algorithm is efficient and simple compared to other algorithms found in the literature. © 2012 IEEE.

A Simple FDTD Algorithm for Simulating EM-Wave Propagation in General Dispersive Anisotropic Material

KAUST Repository

Al-Jabr, Ahmad Ali

2013-03-01

In this paper, an finite-difference time-domain (FDTD) algorithm for simulating propagation of EM waves in anisotropic material is presented. The algorithm is based on the auxiliary differential equation and the general polarization formulation. In anisotropic materials, electric fields are coupled and elements in the permittivity tensor are, in general, multiterm dispersive. The presented algorithm resolves the field coupling using a formulation based on electric polarizations. It also offers a simple procedure for the treatment of multiterm dispersion in the FDTD scheme. The algorithm is tested by simulating wave propagation in 1-D magnetized plasma showing excellent agreement with analytical solutions. Extension of the algorithm to multidimensional structures is straightforward. The presented algorithm is efficient and simple compared to other algorithms found in the literature. © 2012 IEEE.

Analysis of wave equation in electromagnetic field by Proca equation

International Nuclear Information System (INIS)

Pamungkas, Oky Rio; Soeparmi; Cari

2017-01-01

This research is aimed to analyze wave equation for the electric and magnetic field, vector and scalar potential, and continuity equation using Proca equation. Then, also analyze comparison of the solution on Maxwell and Proca equation for scalar potential and electric field, both as a function of distance and constant wave number. (paper)

2D full wave simulation on electromagnetic wave propagation in toroidal plasma

International Nuclear Information System (INIS)

Hojo, Hitoshi; Uruta, Go; Nakayama, Kazunori; Mase, Atsushi

2002-01-01

Global full-wave simulation on electromagnetic wave propagation in toroidal plasma with an external magnetic field imaging a tokamak configuration is performed in two dimensions. The temporal behavior of an electromagnetic wave launched into plasma from a wave-guiding region is obtained. (author)

Modeling and simulation of ocean wave propagation using lattice Boltzmann method

Science.gov (United States)

Nuraiman, Dian

2017-10-01

In this paper, we present on modeling and simulation of ocean wave propagation from the deep sea to the shoreline. This requires high computational cost for simulation with large domain. We propose to couple a 1D shallow water equations (SWE) model with a 2D incompressible Navier-Stokes equations (NSE) model in order to reduce the computational cost. The coupled model is solved using the lattice Boltzmann method (LBM) with the lattice Bhatnagar-Gross-Krook (BGK) scheme. Additionally, a special method is implemented to treat the complex behavior of free surface close to the shoreline. The result shows the coupled model can reduce computational cost significantly compared to the full NSE model.

Bending wave propagation of carbon nanotubes in a bi-parameter elastic matrix

International Nuclear Information System (INIS)

Wu, J.-X.; Li, X.-F.; Tang, G.-J.

2012-01-01

This article studies transverse waves propagating in carbon nanotubes (CNTs) embedded in a surrounding medium. The CNTs are modeled as a nonlocal elastic beam, whereas the surrounding medium is modeled as a bi-parameter elastic medium. When taking into account the effect of rotary inertia of cross-section, a governing equation is acquired. A comparison of wave speeds using the Rayleigh and Euler-Bernoulli theories of beams with the results of molecular dynamics simulation indicates that the nonlocal Rayleigh beam model is more adequate to describe flexural waves in CNTs than the nonlocal Euler-Bernoulli model. The influences of the surrounding medium and rotary inertia on the phase speed for single-walled and double-walled CNTs are analyzed. Obtained results turn out that the surrounding medium plays a dominant role for lower wave numbers, while rotary inertia strongly affects the phase speed for higher wave numbers.

Bending wave propagation of carbon nanotubes in a bi-parameter elastic matrix

Energy Technology Data Exchange (ETDEWEB)

Wu, J.-X. [School of Civil Engineering, Central South University, Changsha, Hunan 410075 (China); Li, X.-F., E-mail: xfli25@yahoo.com.cn [School of Civil Engineering, Central South University, Changsha, Hunan 410075 (China); Tang, G.-J. [College of Aerospace and Materials Engineering, National University of Defense Technology, Changsha 410073 (China)

2012-02-15

This article studies transverse waves propagating in carbon nanotubes (CNTs) embedded in a surrounding medium. The CNTs are modeled as a nonlocal elastic beam, whereas the surrounding medium is modeled as a bi-parameter elastic medium. When taking into account the effect of rotary inertia of cross-section, a governing equation is acquired. A comparison of wave speeds using the Rayleigh and Euler-Bernoulli theories of beams with the results of molecular dynamics simulation indicates that the nonlocal Rayleigh beam model is more adequate to describe flexural waves in CNTs than the nonlocal Euler-Bernoulli model. The influences of the surrounding medium and rotary inertia on the phase speed for single-walled and double-walled CNTs are analyzed. Obtained results turn out that the surrounding medium plays a dominant role for lower wave numbers, while rotary inertia strongly affects the phase speed for higher wave numbers.

Using the gauge condition to simplify the elastodynamic analysis of guided wave propagation

Directory of Open Access Journals (Sweden)

Md Yeasin BHUIYAN

2016-09-01

Full Text Available In this article, gauge condition in elastodynamics is explored more to revive its potential capability of simplifying wave propagation problems in elastic medium. The inception of gauge condition in elastodynamics happens from the Navier-Lame equations upon application of Helmholtz theorem. In order to solve the elastic wave problems by potential function approach, the gauge condition provides the necessary conditions for the potential functions. The gauge condition may be considered as the superposition of the separate gauge conditions of Lamb waves and shear horizontal (SH guided waves respectively, and thus, it may be resolved into corresponding gauges of Lamb waves and SH waves. The manipulation and proper choice of the gauge condition does not violate the classical solutions of elastic waves in plates; rather, it simplifies the problems. The gauge condition allows to obtain the analytical solution of complicated problems in a simplified manner.

Propagation of acoustic waves in a one-dimensional macroscopically inhomogeneous poroelastic material.

Science.gov (United States)

Gautier, G; Kelders, L; Groby, J P; Dazel, O; De Ryck, L; Leclaire, P

2011-09-01

Wave propagation in macroscopically inhomogeneous porous materials has received much attention in recent years. The wave equation, derived from the alternative formulation of Biot's theory of 1962, was reduced and solved recently in the case of rigid frame inhomogeneous porous materials. This paper focuses on the solution of the full wave equation in which the acoustic and the elastic properties of the poroelastic material vary in one-dimension. The reflection coefficient of a one-dimensional macroscopically inhomogeneous porous material on a rigid backing is obtained numerically using the state vector (or the so-called Stroh) formalism and Peano series. This coefficient can then be used to straightforwardly calculate the scattered field. To validate the method of resolution, results obtained by the present method are compared to those calculated by the classical transfer matrix method at both normal and oblique incidence and to experimental measurements at normal incidence for a known two-layers porous material, considered as a single inhomogeneous layer. Finally, discussion about the absorption coefficient for various inhomogeneity profiles gives further perspectives. © 2011 Acoustical Society of America

TWO-DIMENSIONAL MODELLING OF ACCIDENTAL FLOOD WAVES PROPAGATION

OpenAIRE

Lorand Catalin STOENESCU

2011-01-01

The study presented in this article describes a modern modeling methodology of the propagation of accidental flood waves in case a dam break; this methodology is applied in Romania for the first time for the pilot project „Breaking scenarios of Poiana Uzului dam”. The calculation programs used help us obtain a bidimensional calculation (2D) of the propagation of flood waves, taking into consideration the diminishing of the flood wave on a normal direction to the main direction; this diminishi...

A wave propagation matrix method in semiclassical theory

International Nuclear Information System (INIS)

Lee, S.Y.; Takigawa, N.

1977-05-01

A wave propagation matrix method is used to derive the semiclassical formulae of the multiturning point problem. A phase shift matrix and a barrier transformation matrix are introduced to describe the processes of a particle travelling through a potential well and crossing a potential barrier respectively. The wave propagation matrix is given by the products of phase shift matrices and barrier transformation matrices. The method to study scattering by surface transparent potentials and the Bloch wave in solids is then applied

Wave propagation and absorption in the electron cyclotron frequency range for TCA and TCV machines

International Nuclear Information System (INIS)

Cardinali, A.

1990-01-01

The main theoretical aspects of the propagation and absorption of electron cyclotron frequency waves are reviewed and applied to TCA and TCV tokamak plasmas. In particular the electromagnetic cold dispersion relation is solved analytically and numerically in order to recall the basic properties of mode propagation and to calculate the ray-trajectories by means of geometric optics. A numerical code which integrates the coupled first order differential ray-equations, has been developed and applied to the cases of interest. (author) 4 figs., 23 refs

The effects of nonlinear wave propagation on the stability of inertial cavitation

International Nuclear Information System (INIS)

Sinden, D; Stride, E; Saffari, N

2009-01-01

In the context of forecasting temperature and pressure fields generated by high-intensity focussed ultrasound, the accuracy of predictive models is critical for the safety and efficacy of treatment. In such fields 'inertial' cavitation is often observed. Classically, estimations of cavitation thresholds have been based on the assumption that the incident wave at the surface of a bubble is the same as in the far-field, neglecting the effect of nonlinear wave propagation. By modelling the incident wave as a solution to Burgers' equation using weak shock theory, the effects of nonlinear wave propagation on inertial cavitation are investigated using both numerical and analytical techniques. From radius-time curves for a single bubble, it is observed that there is a reduction in the maximum size of a bubble undergoing inertial cavitation and that the inertial collapse occurs earlier in contrast with the classical case. Corresponding stability thresholds for a bubble whose initial radius is slightly below the critical Blake radius are calculated, providing a lower bound for the onset of instability. Bifurcation diagrams and frequency-response curves are presented associated with the loss of stability. The consequences and physical implications of the results are discussed with respect to the classical results.

Nonlinear wave equation with intrinsic wave particle dualism

International Nuclear Information System (INIS)

Klein, J.J.

1976-01-01

A nonlinear wave equation derived from the sine-Gordon equation is shown to possess a variety of solutions, the most interesting of which is a solution that describes a wave packet travelling with velocity usub(e) modulating a carrier wave travelling with velocity usub(c). The envelop and carrier wave speeds agree precisely with the group and phase velocities found by de Broglie for matter waves. No spreading is exhibited by the soliton, so that it behaves exactly like a particle in classical mechanics. Moreover, the classically computed energy E of the disturbance turns out to be exactly equal to the frequency ω of the carrier wave, so that the Planck relation is automatically satisfied without postulating a particle-wave dualism. (author)

Wave energy converter effects on wave propagation: A sensitivity study in Monterey Bay, CA

Science.gov (United States)

Chang, G.; Jones, C. A.; Roberts, J.; Magalen, J.; Ruehl, K.; Chartrand, C.

2014-12-01

The development of renewable offshore energy in the United States is growing rapidly and wave energy is one of the largest resources currently being evaluated. The deployment of wave energy converter (WEC) arrays required to harness this resource could feasibly number in the hundreds of individual devices. The WEC arrays have the potential to alter nearshore wave propagation and circulation patterns and ecosystem processes. As the industry progresses from pilot- to commercial-scale it is important to understand and quantify the effects of WECs on the natural nearshore processes that support a local, healthy ecosystem. To help accelerate the realization of commercial-scale wave power, predictive modeling tools have been developed and utilized to evaluate the likelihood of environmental impact. At present, direct measurements of the effects of different types of WEC arrays on nearshore wave propagation are not available; therefore wave model simulations provide the groundwork for investigations of the sensitivity of model results to prescribed WEC characteristics over a range of anticipated wave conditions. The present study incorporates a modified version of an industry standard wave modeling tool, SWAN (Simulating WAves Nearshore), to simulate wave propagation through a hypothetical WEC array deployment site on the California coast. The modified SWAN, referred to as SNL-SWAN, incorporates device-specific WEC power take-off characteristics to more accurately evaluate a WEC device's effects on wave propagation. The primary objectives were to investigate the effects of a range of WEC devices and device and array characteristics (e.g., device spacing, number of WECs in an array) on nearshore wave propagation using SNL-SWAN model simulations. Results showed that significant wave height was most sensitive to variations in WEC device type and size and the number of WEC devices in an array. Locations in the lee centerline of the arrays in each modeled scenario showed the

Electromagnetic wave propagation in relativistic magnetized plasmas

International Nuclear Information System (INIS)

Weiss, I.

1985-01-01

An improved mathematical technique and a new code for deriving the conductivity tensor for collisionless plasmas have been developed. The method is applicable to a very general case, including both hot (relativistic) and cold magnetized plasmas, with only isotropic equilibrium distributions being considered here. The usual derivation starts from the relativistic Vlasov equation and leads to an integration over an infinite sum of Bessel functions which has to be done numerically. In the new solution the integration is carried out over a product of two Bessel functions only. This reduces the computing time very significantly. An added advantage over existing codes is our capability to perform the computations for waves propagating obliquely to the magnetic field. Both improvements greatly facilitate investigations of properties of the plasma under conditions hitherto unexplored

High Energy Laser Beam Propagation in the Atmosphere: The Integral Invariants of the Nonlinear Parabolic Equation and the Method of Moments

Science.gov (United States)

Manning, Robert M.

2012-01-01

The method of moments is used to define and derive expressions for laser beam deflection and beam radius broadening for high-energy propagation through the Earth s atmosphere. These expressions are augmented with the integral invariants of the corresponding nonlinear parabolic equation that describes the electric field of high-energy laser beam to propagation to yield universal equations for the aforementioned quantities; the beam deflection is a linear function of the propagation distance whereas the beam broadening is a quadratic function of distance. The coefficients of these expressions are then derived from a thin screen approximation solution of the nonlinear parabolic equation to give corresponding analytical expressions for a target located outside the Earth s atmospheric layer. These equations, which are graphically presented for a host of propagation scenarios, as well as the thin screen model, are easily amenable to the phase expansions of the wave front for the specification and design of adaptive optics algorithms to correct for the inherent phase aberrations. This work finds application in, for example, the analysis of beamed energy propulsion for space-based vehicles.

Wave propagation of spectral energy content in a granular chain

Directory of Open Access Journals (Sweden)

Shrivastava Rohit Kumar

2017-01-01

Full Text Available A mechanical wave is propagation of vibration with transfer of energy and momentum. Understanding the spectral energy characteristics of a propagating wave through disordered granular media can assist in understanding the overall properties of wave propagation through inhomogeneous materials like soil. The study of these properties is aimed at modeling wave propagation for oil, mineral or gas exploration (seismic prospecting or non-destructive testing of the internal structure of solids. The focus is on the total energy content of a pulse propagating through an idealized one-dimensional discrete particle system like a mass disordered granular chain, which allows understanding the energy attenuation due to disorder since it isolates the longitudinal P-wave from shear or rotational modes. It is observed from the signal that stronger disorder leads to faster attenuation of the signal. An ordered granular chain exhibits ballistic propagation of energy whereas, a disordered granular chain exhibits more diffusive like propagation, which eventually becomes localized at long time periods. For obtaining mean-field macroscopic/continuum properties, ensemble averaging has been used, however, such an ensemble averaged spectral energy response does not resolve multiple scattering, leading to loss of information, indicating the need for a different framework for micro-macro averaging.

Effects of positron density and temperature on ion-acoustic solitary waves in a magnetized electron-positron-ion plasma: Oblique propagation

International Nuclear Information System (INIS)

Esfandyari-Kalejahi, A.; Akbari-Moghanjoughi, M.; Mehdipoor, M.

2009-01-01

Ion-acoustic (IA) solitary waves are investigated in a magnetized three-component plasma consisting of cold ions, isothermal hot electrons, and positrons. The basic set of fluid equations is reduced to the Korteweg de Vries equation using the standard reductive perturbation (multiple-scale) technique. Theoretical and numerical analyses confirm significant effects of the presence of positrons and the dependence of the electron to positron temperature ratio on the amplitude and the width of IA solitary waves. It is shown that the rarefactive and compressive IA solitary excitations can propagate when the propagation angle θ satisfies 0≤θ 0 , whereas their width depends strictly on B 0 . The numerical analysis has been done based on the typical numerical data from a pulsar magnetosphere.

Electron trapping in the electrosound solitary wave for propagation of high intensity laser in a relativistic plasma

International Nuclear Information System (INIS)

Heidari, E; Aslaninejad, M; Eshraghi, H

2010-01-01

Using a set of relativistic equations for plasmas with warm electrons and cold ions, we have investigated the effects of trapped electrons in the propagation of an electrosound wave and discussed the possibility of the formation of electromagnetic solitons in a plasma. The effective potential energy and deviations of the electron and ion number densities in this relativistic model have been found. We have obtained the governing equations for the amplitude of the HF field with relativistic corrections. In order to show the destructive impact of the trapped electrons on the solitary wave, a relativistic effective potential and the governing equation have been found. It is shown that for certain values of the parameters the condition of localization of the HF amplitude is violated. In addition, it is shown that as the flow velocity of the plasma changes, the shape of the solitary wave shows two opposing behaviours, depending on whether the solitary wave velocity is larger than the flow velocity or smaller. Also, the existence of stationary solitary waves which are prohibited for nonrelativistic plasma has been predicted. Finally, we have obtained the Korteweg-de Vries equation showing the relativistic, trapping and nonlinearity effects.

A single-sided representation for the homogeneous Green's function of a unified scalar wave equation.

Science.gov (United States)

Wapenaar, Kees

2017-06-01

A unified scalar wave equation is formulated, which covers three-dimensional (3D) acoustic waves, 2D horizontally-polarised shear waves, 2D transverse-electric EM waves, 2D transverse-magnetic EM waves, 3D quantum-mechanical waves and 2D flexural waves. The homogeneous Green's function of this wave equation is a combination of the causal Green's function and its time-reversal, such that their singularities at the source position cancel each other. A classical representation expresses this homogeneous Green's function as a closed boundary integral. This representation finds applications in holographic imaging, time-reversed wave propagation and Green's function retrieval by cross correlation. The main drawback of the classical representation in those applications is that it requires access to a closed boundary around the medium of interest, whereas in many practical situations the medium can be accessed from one side only. Therefore, a single-sided representation is derived for the homogeneous Green's function of the unified scalar wave equation. Like the classical representation, this single-sided representation fully accounts for multiple scattering. The single-sided representation has the same applications as the classical representation, but unlike the classical representation it is applicable in situations where the medium of interest is accessible from one side only.

Magnetic Field Effects and Electromagnetic Wave Propagation in Highly Collisional Plasmas.

Science.gov (United States)

Bozeman, Steven Paul

The homogeneity and size of radio frequency (RF) and microwave driven plasmas are often limited by insufficient penetration of the electromagnetic radiation. To investigate increasing the skin depth of the radiation, we consider the propagation of electromagnetic waves in a weakly ionized plasma immersed in a steady magnetic field where the dominant collision processes are electron-neutral and ion-neutral collisions. Retaining both the electron and ion dynamics, we have adapted the theory for cold collisionless plasmas to include the effects of these collisions and obtained the dispersion relation at arbitrary frequency omega for plane waves propagating at arbitrary angles with respect to the magnetic field. We discuss in particular the cases of magnetic field enhanced wave penetration for parallel and perpendicular propagation, examining the experimental parameters which lead to electromagnetic wave propagation beyond the collisional skin depth. Our theory predicts that the most favorable scaling of skin depth with magnetic field occurs for waves propagating nearly parallel to B and for omega << Omega_{rm e} where Omega_{rm e} is the electron cyclotron frequency. The scaling is less favorable for propagation perpendicular to B, but the skin depth does increase for this case as well. Still, to achieve optimal wave penetration, we find that one must design the plasma configuration and antenna geometry so that one generates primarily the appropriate angles of propagation. We have measured plasma wave amplitudes and phases using an RF magnetic probe and densities using Stark line broadening. These measurements were performed in inductively coupled plasmas (ICP's) driven with a standard helical coil, a reverse turn (Stix) coil, and a flat spiral coil. Density measurements were also made in a microwave generated plasma. The RF magnetic probe measurements of wave propagation in a conventional ICP with wave propagation approximately perpendicular to B show an increase in

Effect of material parameters on stress wave propagation during fast upsetting

Institute of Scientific and Technical Information of China (English)

WANG Zhong-jin; CHENG Li-dong

2008-01-01

Based'on a dynamic analysis method and an explicit algorithm, a dynamic explicit finite element code was developed for modeling the fast upsetting process of block under drop hammer impact, in which the hammer velocity during the deformation was calculated by energy conservation law according to the operating principle of hammer equipment. The stress wave propagation and its effect on the deformation were analyzed by the stress and strain distributions. Industrial pure lead, oxygen-free high-conductivity (OFHC) copper and 7039 aluminum alloy were chosen to investigate the effect of material parameters on the stress wave propagation. The results show that the stress wave propagates from top to bottom of block, and then reflects back when it reaches the bottom surface. After that, stress wave propagates and reflects repeatedly between the upper surface and bottom surface. The stress wave propagation has a significant effect on the deformation at the initial stage, and then becomes weak at the middle-final stage. When the ratio of elastic modulus or the slope of stress-strain curve to mass density becomes larger, the velocity of stress wave propagation increases, and the influence of stress wave on the deformation becomes small.

Statistics of peak overpressure and shock steepness for linear and nonlinear N-wave propagation in a kinematic turbulence.

Science.gov (United States)

Yuldashev, Petr V; Ollivier, Sébastien; Karzova, Maria M; Khokhlova, Vera A; Blanc-Benon, Philippe

2017-12-01

Linear and nonlinear propagation of high amplitude acoustic pulses through a turbulent layer in air is investigated using a two-dimensional KZK-type (Khokhlov-Zabolotskaya-Kuznetsov) equation. Initial waves are symmetrical N-waves with shock fronts of finite width. A modified von Kármán spectrum model is used to generate random wind velocity fluctuations associated with the turbulence. Physical parameters in simulations correspond to previous laboratory scale experiments where N-waves with 1.4 cm wavelength propagated through a turbulence layer with the outer scale of about 16 cm. Mean value and standard deviation of peak overpressure and shock steepness, as well as cumulative probabilities to observe amplified peak overpressure and shock steepness, are analyzed. Nonlinear propagation effects are shown to enhance pressure level in random foci for moderate initial amplitudes of N-waves thus increasing the probability to observe highly peaked waveforms. Saturation of the pressure level is observed for stronger nonlinear effects. It is shown that in the linear propagation regime, the turbulence mainly leads to the smearing of shock fronts, thus decreasing the probability to observe high values of steepness, whereas nonlinear effects dramatically increase the probability to observe steep shocks.

Impact induced solitary wave propagation through a woodpile structure

International Nuclear Information System (INIS)

Kore, R; Waychal, A; Yadav, P; Shelke, A; Agarwal, S; Sahoo, N; Uddin, Ahsan

2016-01-01

In this paper, we investigate solitary wave propagation through a one-dimensional woodpile structure excited by low and high velocity impact. Woodpile structures are a sub-class of granular metamaterial, which supports propagation of nonlinear waves. Hertz contact law governs the behavior of the solitary wave propagation through the granular media. Towards an experimental study, a woodpile structure was fabricated by orthogonally stacking cylindrical rods. A shock tube facility has been developed to launch an impactor on the woodpile structure at a velocity of 30 m s −1 . Embedded granular chain sensors were fabricated to study the behavior of the solitary wave. The impact induced stress wave is studied to investigate solitary wave parameters, i.e. contact force, contact time, and solitary wave velocity. With the aid of the experimental setup, numerical simulations, and a theoretical solution based on the long wavelength approximation, formation of the solitary wave in the woodpile structure is validated to a reasonable degree of accuracy. The nondispersive and compact supported solitary waves traveling at sonic wave velocity offer unique properties that could be leveraged for application in nondestructive testing and structural health monitoring. (paper)

Propagation of Finite Amplitude Sound in Multiple Waveguide Modes.

Science.gov (United States)

van Doren, Thomas Walter

1993-01-01

This dissertation describes a theoretical and experimental investigation of the propagation of finite amplitude sound in multiple waveguide modes. Quasilinear analytical solutions of the full second order nonlinear wave equation, the Westervelt equation, and the KZK parabolic wave equation are obtained for the fundamental and second harmonic sound fields in a rectangular rigid-wall waveguide. It is shown that the Westervelt equation is an acceptable approximation of the full nonlinear wave equation for describing guided sound waves of finite amplitude. A system of first order equations based on both a modal and harmonic expansion of the Westervelt equation is developed for waveguides with locally reactive wall impedances. Fully nonlinear numerical solutions of the system of coupled equations are presented for waveguides formed by two parallel planes which are either both rigid, or one rigid and one pressure release. These numerical solutions are compared to finite -difference solutions of the KZK equation, and it is shown that solutions of the KZK equation are valid only at frequencies which are high compared to the cutoff frequencies of the most important modes of propagation (i.e., for which sound propagates at small grazing angles). Numerical solutions of both the Westervelt and KZK equations are compared to experiments performed in an air-filled, rigid-wall, rectangular waveguide. Solutions of the Westervelt equation are in good agreement with experiment for low source frequencies, at which sound propagates at large grazing angles, whereas solutions of the KZK equation are not valid for these cases. At higher frequencies, at which sound propagates at small grazing angles, agreement between numerical solutions of the Westervelt and KZK equations and experiment is only fair, because of problems in specifying the experimental source condition with sufficient accuracy.

Some considerations of wave propagation

Science.gov (United States)

Verdonk, P. L. F. M.

The meaning of group velocity and its relation to conserved quantities are demonstrated. The origin of wave dispersion in terms of nonlocal and relaxation phenomena are clarified. The character of a wave described by an equation with a general type of nonlinearity and general dispersion terms is explained. The steepening of a wave flank and the occurrence of stationary waves are discussed.

Waves propagating over a two-layer porous barrier on a seabed

Science.gov (United States)

Lin, Qiang; Meng, Qing-rui; Lu, Dong-qiang

2018-05-01

A research of wave propagation over a two-layer porous barrier, each layer of which is with different values of porosity and friction, is conducted with a theoretical model in the frame of linear potential flow theory. The model is more appropriate when the seabed consists of two different properties, such as rocks and breakwaters. It is assumed that the fluid is inviscid and incompressible and the motion is irrotational. The wave numbers in the porous region are complex ones, which are related to the decaying and propagating behaviors of wave modes. With the aid of the eigenfunction expansions, a new inner product of the eigenfunctions in the two-layer porous region is proposed to simplify the calculation. The eigenfunctions, under this new definition, possess the orthogonality from which the expansion coefficients can be easily deduced. Selecting the optimum truncation of the series, we derive a closed system of simultaneous linear equations for the same number of the unknown reflection and transmission coefficients. The effects of several physical parameters, including the porosity, friction, width, and depth of the porous barrier, on the dispersion relation, reflection and transmission coefficients are discussed in detail through the graphical representations of the solutions. It is concluded that these parameters have certain impacts on the reflection and transmission energy.

Parsimonious wave-equation travel-time inversion for refraction waves

KAUST Repository

Fu, Lei

2017-02-14

We present a parsimonious wave-equation travel-time inversion technique for refraction waves. A dense virtual refraction dataset can be generated from just two reciprocal shot gathers for the sources at the endpoints of the survey line, with N geophones evenly deployed along the line. These two reciprocal shots contain approximately 2N refraction travel times, which can be spawned into O(N2) refraction travel times by an interferometric transformation. Then, these virtual refraction travel times are used with a source wavelet to create N virtual refraction shot gathers, which are the input data for wave-equation travel-time inversion. Numerical results show that the parsimonious wave-equation travel-time tomogram has about the same accuracy as the tomogram computed by standard wave-equation travel-time inversion. The most significant benefit is that a reciprocal survey is far less time consuming than the standard refraction survey where a source is excited at each geophone location.

Uniqueness and existence result around lax-milgram lemma: application to electromagnetic waves propagation in tokamak plasmas

International Nuclear Information System (INIS)

Sebelin, E.; Peysson, Y.; Litaudon, X.; Moreau, D.

1997-09-01

In the context of complex Hilbert spaces is proved, around Lax-Milgram lemma, the existence and uniqueness of solutions associated to a class of stationary variational problems. This result is applied to the study of variational problems from the propagation equation of time-harmonic electromagnetic waves in confined tokamak plasmas. (author)

Integral representations of solutions of the wave equation based on relativistic wavelets

International Nuclear Information System (INIS)

Perel, Maria; Gorodnitskiy, Evgeny

2012-01-01

A representation of solutions of the wave equation with two spatial coordinates in terms of localized elementary ones is presented. Elementary solutions are constructed from four solutions with the help of transformations of the affine Poincaré group, i.e. with the help of translations, dilations in space and time and Lorentz transformations. The representation can be interpreted in terms of the initial-boundary value problem for the wave equation in a half-plane. It gives the solution as an integral representation of two types of solutions: propagating localized solutions running away from the boundary under different angles and packet-like surface waves running along the boundary and exponentially decreasing away from the boundary. Properties of elementary solutions are discussed. A numerical investigation of coefficients of the decomposition is carried out. An example of the decomposition of the field created by sources moving along a line with different speeds is considered, and the dependence of coefficients on speeds of sources is discussed. (paper)

Arbitrary amplitude electrostatic wave propagation in a magnetized dense plasma containing helium ions and degenerate electrons

Science.gov (United States)

Mahmood, S.; Sadiq, Safeer; Haque, Q.; Ali, Munazza Z.

2016-06-01

The obliquely propagating arbitrary amplitude electrostatic wave is studied in a dense magnetized plasma having singly and doubly charged helium ions with nonrelativistic and ultrarelativistic degenerate electrons pressures. The Fermi temperature for ultrarelativistic degenerate electrons described by N. M. Vernet [(Cambridge University Press, Cambridge, 2007), p. 57] is used to define ion acoustic speed in ultra-dense plasmas. The pseudo-potential approach is used to solve the fully nonlinear set of dynamic equations for obliquely propagating electrostatic waves in a dense magnetized plasma containing helium ions. The upper and lower Mach number ranges for the existence of electrostatic solitons are found which depends on the obliqueness of the wave propagation with respect to applied magnetic field and charge number of the helium ions. It is found that only compressive (hump) soliton structures are formed in all the cases and only subsonic solitons are formed for a singly charged helium ions plasma case with nonrelativistic degenerate electrons. Both subsonic and supersonic soliton hump structures are formed for doubly charged helium ions with nonrelativistic degenerate electrons and ultrarelativistic degenerate electrons plasma case containing singly as well as doubly charged helium ions. The effect of propagation direction on the soliton amplitude and width of the electrostatic waves is also presented. The numerical plots are also shown for illustration using dense plasma parameters of a compact star (white dwarf) from literature.

Closure of multi-fluid and kinetic equations for cyclotron-resonant interactions of solar wind ions with Alfvén waves

Directory of Open Access Journals (Sweden)

E. Marsch

1998-01-01

Full Text Available Based on quasilinear theory, a closure scheme for anisotropic multi-component fluid equations is developed for the wave-particle interactions of ions with electromagnetic Alfvén and ion-cyclotron waves propagating along the mean magnetic field. Acceleration and heating rates are calculated. They may be used in the multi-fluid momentum and energy equations as anomalous transport terms. The corresponding evolution equation for the average wave spectrum is established, and the effective growth/damping rate for the spectrum is calculated. Given a simple power-law spectrum, an anomalous collision frequency can be derived which depends on the slope and average intensity of the spectrum, and on the gyrofrequency and the differential motion (with respect to the wave frame of the actual ion species considered. The wave-particle interaction terms attain simple forms resembling the ones for collisional friction and temperature anisotropy relaxation (due to pitch angle scattering with collision rates that are proportional to the gyrofrequency but diminished substantially by the relative wave energy or the fluctuation level with respect the background field. In addition, a set of quasilinear diffusion equations is derived for the reduced (with respect to the perpendicular velocity component velocity distribution functions (VDFs, as they occur in the wave dispersion equation and the related dielectric function for parallel propagation. These reduced VDFs allow one to describe adequately the most prominent observed features, such as an ion beam and temperature anisotropy, in association with the resonant interactions of the particles with the waves on a kinetic level, yet have the advantage of being only dependent upon the parallel velocity component.

Uniqueness and existence result around lax-milgram lemma: application to electromagnetic waves propagation in tokamak plasmas

Energy Technology Data Exchange (ETDEWEB)

Sebelin, E.; Peysson, Y.; Litaudon, X.; Moreau, D. [Association Euratom-CEA, CEA Cadarache, 13 - Saint-Paul-lez-Durance (France). Dept. de Recherches sur la Fusion Controlee; Miellou, J.C. [Besancon Univ., 25 (France). Laboratoire d`Analyse Numerique; Lafitte, O. [CEA Limeil, 94 - Villeneuve-Saint-Georges (France)

1997-09-01

In the context of complex Hilbert spaces is proved, around Lax-Milgram lemma, the existence and uniqueness of solutions associated to a class of stationary variational problems. This result is applied to the study of variational problems from the propagation equation of time-harmonic electromagnetic waves in confined tokamak plasmas. (author) 21 refs.

Computer modeling of inelastic wave propagation in porous rock

International Nuclear Information System (INIS)

Cheney, J.A.; Schatz, J.F.; Snell, C.

1979-01-01

Computer modeling of wave propagation in porous rock has several important applications. Among them are prediction of fragmentation and permeability changes to be caused by chemical explosions used for in situ resource recovery, and the understanding of nuclear explosion effects such as seismic wave generation, containment, and site hardness. Of interest in all these applications are the distance from the source to which inelastic effects persist and the amount of porosity change within the inelastic region. In order to study phenomena related to these applications, the Cam Clay family of models developed at Cambridge University was used to develop a similar model that is applicable to wave propagation in porous rock. That model was incorporated into a finite-difference wave propagation computer code SOC. 10 figures, 1 table

Characteristics of coupled acoustic wave propagation in metal pipe

International Nuclear Information System (INIS)

Kim, Ho Wuk; Kim, Min Soo; Lee, Sang Kwon

2008-01-01

The circular cylinder pipes are used in the many industrial areas. In this paper, the acoustic wave propagation in the pipe containing gas is researched. First of all, the theory for the coupled acoustic wave propagation in a pipe is investigated. Acoustic wave propagation in pipe can not be occurred independently between the wave of the fluid and the shell. It requires complicated analysis. However, as a special case, the coupled wave in a high density pipe containing a light density medium is corresponded closely to the uncoupled in-vacuo shell waves and to the rigid-walled duct fluid waves. The coincidence frequencies of acoustic and shell modes contribute to the predominant energy transmission. The coincidence frequency means the frequency corresponding to the coincidence of the wavenumber in both acoustic and shell. In this paper, it is assumed that the internal medium is much lighter than the pipe shell. After the uncoupled acoustic wave in the internal medium and uncoupled shell wave are considered, the coincidence frequencies are found. The analysis is successfully confirmed by the verification of the experiment using the real long steel pipe. This work verifies that the coupled wave characteristic of the shell and the fluid is occurred as predominant energy transmission at the coincidence frequencies

Acoustic propagation operators for pressure waves on an arbitrarily curved surface in a homogeneous medium

Science.gov (United States)

Sun, Yimin; Verschuur, Eric; van Borselen, Roald

2018-03-01

The Rayleigh integral solution of the acoustic Helmholtz equation in a homogeneous medium can only be applied when the integral surface is a planar surface, while in reality almost all surfaces where pressure waves are measured exhibit some curvature. In this paper we derive a theoretically rigorous way of building propagation operators for pressure waves on an arbitrarily curved surface. Our theory is still based upon the Rayleigh integral, but it resorts to matrix inversion to overcome the limitations faced by the Rayleigh integral. Three examples are used to demonstrate the correctness of our theory - propagation of pressure waves acquired on an arbitrarily curved surface to a planar surface, on an arbitrarily curved surface to another arbitrarily curved surface, and on a spherical cap to a planar surface, and results agree well with the analytical solutions. The generalization of our method for particle velocities and the calculation cost of our method are also discussed.

High frequency guided wave propagation in monocrystalline silicon wafers

Science.gov (United States)

Pizzolato, Marco; Masserey, Bernard; Robyr, Jean-Luc; Fromme, Paul

2017-04-01

Monocrystalline silicon wafers are widely used in the photovoltaic industry for solar panels with high conversion efficiency. The cutting process can introduce micro-cracks in the thin wafers and lead to varying thickness. High frequency guided ultrasonic waves are considered for the structural monitoring of the wafers. The anisotropy of the monocrystalline silicon leads to variations of the wave characteristics, depending on the propagation direction relative to the crystal orientation. Full three-dimensional Finite Element simulations of the guided wave propagation were conducted to visualize and quantify these effects for a line source. The phase velocity (slowness) and skew angle of the two fundamental Lamb wave modes (first anti-symmetric mode A0 and first symmetric mode S0) for varying propagation directions relative to the crystal orientation were measured experimentally. Selective mode excitation was achieved using a contact piezoelectric transducer with a custom-made wedge and holder to achieve a controlled contact pressure. The out-of-plane component of the guided wave propagation was measured using a noncontact laser interferometer. Good agreement was found with the simulation results and theoretical predictions based on nominal material properties of the silicon wafer.

Propagation of electromagnetic waves in a weakly ionized dusty plasma

International Nuclear Information System (INIS)

Jia, Jieshu; Yuan, Chengxun; Gao, Ruilin; Wang, Ying; Liu, Yaoze; Gao, Junying; Zhou, Zhongxiang; Sun, Xiudong; Li, Hui; Wu, Jian; Pu, Shaozhi

2015-01-01

Propagation properties of electromagnetic (EM) waves in weakly ionized dusty plasmas are the subject of this study. Dielectric relation for EM waves propagating at a weakly ionized dusty plasma is derived based on the Boltzmann distribution law while considering the collision and charging effects of dust grains. The propagation properties of EM energy in dusty plasma of rocket exhaust are numerically calculated and studied, utilizing the parameters of rocket exhaust plasma. Results indicate that increase of dust radius and density enhance the reflection and absorption coefficient. High dust radius and density make the wave hardly transmit through the dusty plasmas. Interaction enhancements between wave and dusty plasmas are developed through effective collision frequency improvements. Numerical results coincide with observed results by indicating that GHz band wave communication is effected by dusty plasma as the presence of dust grains significantly affect propagation of EM waves in the dusty plasmas. The results are helpful to analyze the effect of dust in plasmas and also provide a theoretical basis for the experiments. (paper)

Modification of 2-D Time-Domain Shallow Water Wave Equation using Asymptotic Expansion Method

Science.gov (United States)

Khairuman, Teuku; Nasruddin, MN; Tulus; Ramli, Marwan

2018-01-01

Generally, research on the tsunami wave propagation model can be conducted by using a linear model of shallow water theory, where a non-linear side on high order is ignored. In line with research on the investigation of the tsunami waves, the Boussinesq equation model underwent a change aimed to obtain an improved quality of the dispersion relation and non-linearity by increasing the order to be higher. To solve non-linear sides at high order is used a asymptotic expansion method. This method can be used to solve non linear partial differential equations. In the present work, we found that this method needs much computational time and memory with the increase of the number of elements.

Investigating Alfvénic wave propagation in coronal open-field regions

Science.gov (United States)

Morton, R. J.; Tomczyk, S.; Pinto, R.

2015-01-01

The physical mechanisms behind accelerating solar and stellar winds are a long-standing astrophysical mystery, although recent breakthroughs have come from models invoking the turbulent dissipation of Alfvén waves. The existence of Alfvén waves far from the Sun has been known since the 1970s, and recently the presence of ubiquitous Alfvénic waves throughout the solar atmosphere has been confirmed. However, the presence of atmospheric Alfvénic waves does not, alone, provide sufficient support for wave-based models; the existence of counter-propagating Alfvénic waves is crucial for the development of turbulence. Here, we demonstrate that counter-propagating Alfvénic waves exist in open coronal magnetic fields and reveal key observational insights into the details of their generation, reflection in the upper atmosphere and outward propagation into the solar wind. The results enhance our knowledge of Alfvénic wave propagation in the solar atmosphere, providing support and constraints for some of the recent Alfvén wave turbulence models. PMID:26213234

Nonlinear wave equations

CERN Document Server

Li, Tatsien

2017-01-01

This book focuses on nonlinear wave equations, which are of considerable significance from both physical and theoretical perspectives. It also presents complete results on the lower bound estimates of lifespan (including the global existence), which are established for classical solutions to the Cauchy problem of nonlinear wave equations with small initial data in all possible space dimensions and with all possible integer powers of nonlinear terms. Further, the book proposes the global iteration method, which offers a unified and straightforward approach for treating these kinds of problems. Purely based on the properties of solut ions to the corresponding linear problems, the method simply applies the contraction mapping principle.

Directional bending wave propagation in periodically perforated plates

DEFF Research Database (Denmark)

Andreassen, Erik; Manktelow, Kevin; Ruzzene, Massimo

2015-01-01

We report on the investigation of wave propagation in a periodically perforated plate. A unit cell with double-C perforations is selected as a test article suitable to investigate two-dimensional dispersion characteristics, group velocities, and internal resonances. A numerical model, formulated...... using Mindlin plate elements, is developed to predict relevant wave characteristics such as dispersion, and group velocity variation as a function of frequency and direction of propagation. Experimental tests are conducted through a scanning laser vibrometer, which provides full wave field information...... for the design of phononic waveguides with directional and internal resonant characteristics....

Finite element modeling of acoustic wave propagation and energy deposition in bone during extracorporeal shock wave treatment

Science.gov (United States)

Wang, Xiaofeng; Matula, Thomas J.; Ma, Yong; Liu, Zheng; Tu, Juan; Guo, Xiasheng; Zhang, Dong

2013-06-01

It is well known that extracorporeal shock wave treatment is capable of providing a non-surgical and relatively pain free alternative treatment modality for patients suffering from musculoskeletal disorders but do not respond well to conservative treatments. The major objective of current work is to investigate how the shock wave (SW) field would change if a bony structure exists in the path of the acoustic wave. Here, a model of finite element method (FEM) was developed based on linear elasticity and acoustic propagation equations to examine SW propagation and deflection near a mimic musculoskeletal bone. High-speed photography experiments were performed to record cavitation bubbles generated in SW field with the presence of mimic bone. By comparing experimental and simulated results, the effectiveness of FEM model could be verified and strain energy distributions in the bone were also predicted according to numerical simulations. The results show that (1) the SW field will be deflected with the presence of bony structure and varying deflection angles can be observed as the bone shifted up in the z-direction relative to SW geometric focus (F2 focus); (2) SW deflection angels predicted by the FEM model agree well with experimental results obtained from high-speed photographs; and (3) temporal evolutions of strain energy distribution in the bone can also be evaluated based on FEM model, with varied vertical distance between F2 focus and intended target point on the bone surface. The present studies indicate that, by combining MRI/CT scans and FEM modeling work, it is possible to better understand SW propagation characteristics and energy deposition in musculoskeletal structure during extracorporeal shock wave treatment, which is important for standardizing the treatment dosage, optimizing treatment protocols, and even providing patient-specific treatment guidance in clinic.

Propagation of plane waves in a rotating magneto-thermoelastic fiber-reinforced medium under G-N theory

Directory of Open Access Journals (Sweden)

Maity N.

2017-06-01

Full Text Available The article is concernedwith the possibility of plane wave propagation in a rotating elastic medium under the action of magnetic and thermal fields. The material is assumed to be fibre-reinforced with increased stiffness, strength and load bearing capacity. Green and Nagdhi’s concepts of generalized thermoelastic models II and III have been followed in the governing equations expressed in tensor notation. The effects of various parameters of the applied fields on the plane wave velocity have been shown graphically.

Spherical shock-wave propagation in three-dimensional granular packings.

Science.gov (United States)

Xue, Kun; Bai, Chun-Hua

2011-02-01

We investigate numerically the spherical shock-wave propagation in an open dense granular packing perturbed by the sudden expansion of a spherical intruder in the interior of the pack, focusing on the correlation between geometrical fabrics and propagating properties. The measurements of the temporal and spatial variations in a variety of propagating properties define a consistent serrated wave substructure with characteristic length on the orders of particle diameters. Further inspection of particle packing reveals a well-defined particle layering that persists several particle diameters away from the intruder, although its dominant effects are only within one to two diameters. This interface-induced layering not only exactly coincides with the serrated wave profile, but also highlights the competition between two energy transmission mechanisms involving distinct transport speeds. The alternating dominances between these two mechanisms contribute to the nonlinear wave propagation on the particle scale. Moreover, the proliferation of intricate three-dimensional contact force networks suggests the anisotropic stress transmission, which is found to also arise from the localized packing structure in the vicinity of the intruder.

Deep currents in the Gulf of Guinea: along slope propagation of intraseasonal waves

Directory of Open Access Journals (Sweden)

C. Guiavarc'h

2009-05-01

Full Text Available In the Gulf of Guinea, intraseasonal variability is large at the equator and along the coast. Current data on the continental slope near 7.5° S show very energetic biweekly oscillations at 1300 m depth. A high resolution primitive equation numerical model demonstrates that this deep variability is forced by equatorial winds, through the generation of equatorial Yanai waves that propagate eastward and at depth, and then poleward as coastally-trapped waves upon reaching the coast of Africa. Intraseasonal variability is intensified along the coast of the Gulf of Guinea, especially in the 10–20 day period range and at depths between 500 and 1500 m. The kinetic energy distribution is well explained at first order by linear theory. Along the equator, eastward intensification of energy and bottom intensification are in qualitative agreement with vertically propagating Yanai waves, although the signal is influenced by the details of the bathymetry. Along the coast, baroclinic modes 3 to 5 are important close to the equator, and the signal is dominated by lower vertical modes farther south. Additional current meter data on the continental slope near 3° N display an energy profile in the 10–20 day period band that is strikingly different from the one at 7.5° S, with surface intensification rather than bottom intensification and a secondary maximum near 800 m. The model reproduces these features and explains them: the surface intensification in the north is due to the regional wind forcing, and the north-south asymmetry of the deep signal is due to the presence of the zonal African coast near 5° N. A 4 years time series of current measurements at 7.5° S shows that the biweekly oscillations are intermittent and vary from year to year. This intermittency is not well correlated with fluctuations of the equatorial winds and does not seem to be a simple linear response to the wind forcing.

Wave propagation in the magnetosphere of Jupiter

Science.gov (United States)

Liemohn, H. B.

1972-01-01

A systematic procedure is developed for identifying the spatial regimes of various modes of wave propagation in the Jupiter magnetosphere that may be encountered by flyby missions. The Clemmow-Mullaly-Allis (CMA) diagram of plasma physics is utilized to identify the frequency regimes in which different modes of propagation occur in the magnetoplasma. The Gledhill model and the Ioannidis and Brice model of the magnetoplasma are summarized, and configuration-space CMA diagrams are constructed for each model for frequencies from 10 Hz to 1 MHz. The distinctive propagation features, the radio noise regimes, and the wave-particle interactions are discussed. It is concluded that the concentration of plasma in the equatorial plane makes this region of vital importance for radio observations with flyby missions. Local radio noise around the electron cyclotron frequency will probably differ appreciably from its terrestrial counterpart due to the lack of field-line guidance. Hydromagnetic wave properties at frequencies near the ion cyclotron frequency and below will probably be similar to the terrestrial case.

Solitary wave solutions as a signature of the instability in the discrete nonlinear Schroedinger equation

Energy Technology Data Exchange (ETDEWEB)

Arevalo, Edward, E-mail: arevalo@temf.tu-darmstadt.d [Technische Universitaet Darmstadt, Institut fuer Theorie elektromagnetischer Felder, TEMF, Schlossgartenstr. 8, D-64289 Darmstadt (Germany)

2009-09-21

The effect of instability on the propagation of solitary waves along one-dimensional discrete nonlinear Schroedinger equation with cubic nonlinearity is revisited. A self-contained quasicontinuum approximation is developed to derive closed-form expressions for small-amplitude solitary waves. The notion that the existence of nonlinear solitary waves in discrete systems is a signature of the modulation instability is used. With the help of this notion we conjecture that instability effects on moving solitons can be qualitative estimated from the analytical solutions. Results from numerical simulations are presented to support this conjecture.

A comparative study of diffraction of shallow-water waves by high-level IGN and GN equations

Energy Technology Data Exchange (ETDEWEB)

Zhao, B.B. [College of Shipbuilding Engineering, Harbin Engineering University, 150001 Harbin (China); Ertekin, R.C. [Department of Ocean and Resources Engineering, University of Hawai' i, Honolulu, HI 96822 (United States); College of Shipbuilding Engineering, Harbin Engineering University, 150001 Harbin (China); Duan, W.Y., E-mail: duanwenyangheu@hotmail.com [College of Shipbuilding Engineering, Harbin Engineering University, 150001 Harbin (China)

2015-02-15

This work is on the nonlinear diffraction analysis of shallow-water waves, impinging on submerged obstacles, by two related theories, namely the classical Green–Naghdi (GN) equations and the Irrotational Green–Naghdi (IGN) equations, both sets of equations being at high levels and derived for incompressible and inviscid flows. Recently, the high-level Green–Naghdi equations have been applied to some wave transformation problems. The high-level IGN equations have also been used in the last decade to study certain wave propagation problems. However, past works on these theories used different numerical methods to solve these nonlinear and unsteady sets of differential equations and at different levels. Moreover, different physical problems have been solved in the past. Therefore, it has not been possible to understand the differences produced by these two sets of theories and their range of applicability so far. We are thus motivated to make a direct comparison of the results produced by these theories by use of the same numerical method to solve physically the same wave diffraction problems. We focus on comparing these two theories by using similar codes; only the equations used are different but other parts of the codes, such as the wave-maker, damping zone, discretion method, matrix solver, etc., are exactly the same. This way, we eliminate many potential sources of differences that could be produced by the solution of different equations. The physical problems include the presence of various submerged obstacles that can be used for example as breakwaters or to represent the continental shelf. A numerical wave tank is created by placing a wavemaker on one end and a wave absorbing beach on the other. The nonlinear and unsteady sets of differential equations are solved by the finite-difference method. The results are compared with different equations as well as with the available experimental data.

Oblique propagation of nonlinear hydromagnetic waves: One- and two-dimensional behavior

International Nuclear Information System (INIS)

Malara, F.; Elaoufir, J.

1991-01-01

The one- and two-dimensional behavior of obliquely propagating hydromagnetic waves is analyzed by means of analytical theory and numerical simulations. It is shown that the nonlinear evolution of a one-dimensional MHD wave leads to the formation of a rotational discontinuity and a compressive steepened quasi-linearly polarized pulse whose structure is similar to that of a finite amplitude magnetosonic simple wave. For small propagation angles, the pulse mode (fast or slow) depends on the value of β with respect to unity while for large propagation angles the wave mode is fixed by the sign of the initial density-field correlation. The two-dimensional evolution shows that an MHD wave is unstable against a small-amplitude long-wavelength modulation in the direction transverse to the wave propagation direction. A two-dimensional magnetosonic wave solution is found, in which the density fluctuation is driven by the corresponding total pressure fluctuation, exactly as in the one-dimensional simple wave. Along with the steepening effect, the wave experiences both wave front deformation and a self-focusing effect which may eventually lead to the collapse of the wave. The results compare well with observations of MHD waves in the Earth's foreshock and at comets

Radio Wave Propagation Scene Partitioning for High-Speed Rails

Directory of Open Access Journals (Sweden)

Bo Ai

2012-01-01

Full Text Available Radio wave propagation scene partitioning is necessary for wireless channel modeling. As far as we know, there are no standards of scene partitioning for high-speed rail (HSR scenarios, and therefore we propose the radio wave propagation scene partitioning scheme for HSR scenarios in this paper. Based on our measurements along the Wuhan-Guangzhou HSR, Zhengzhou-Xian passenger-dedicated line, Shijiazhuang-Taiyuan passenger-dedicated line, and Beijing-Tianjin intercity line in China, whose operation speeds are above 300 km/h, and based on the investigations on Beijing South Railway Station, Zhengzhou Railway Station, Wuhan Railway Station, Changsha Railway Station, Xian North Railway Station, Shijiazhuang North Railway Station, Taiyuan Railway Station, and Tianjin Railway Station, we obtain an overview of HSR propagation channels and record many valuable measurement data for HSR scenarios. On the basis of these measurements and investigations, we partitioned the HSR scene into twelve scenarios. Further work on theoretical analysis based on radio wave propagation mechanisms, such as reflection and diffraction, may lead us to develop the standard of radio wave propagation scene partitioning for HSR. Our work can also be used as a basis for the wireless channel modeling and the selection of some key techniques for HSR systems.

Mathematical Modeling of Torsional Surface Wave Propagation in a Non-Homogeneous Transverse Isotropic Elastic Solid Semi-Infinite Medium Under a Layer

Science.gov (United States)

Sethi, M.; Sharma, A.; Vasishth, A.

2017-05-01

The present paper deals with the mathematical modeling of the propagation of torsional surface waves in a non-homogeneous transverse isotropic elastic half-space under a rigid layer. Both rigidities and density of the half-space are assumed to vary inversely linearly with depth. Separation of variable method has been used to get the analytical solutions for the dispersion equation of the torsional surface waves. Also, the effects of nonhomogeneities on the phase velocity of torsional surface waves have been shown graphically. Also, dispersion equations have been derived for some particular cases, which are in complete agreement with some classical results.

Travelling wave solutions of generalized coupled Zakharov–Kuznetsov and dispersive long wave equations

Directory of Open Access Journals (Sweden)

M. Arshad

Full Text Available In this manuscript, we constructed different form of new exact solutions of generalized coupled Zakharov–Kuznetsov and dispersive long wave equations by utilizing the modified extended direct algebraic method. New exact traveling wave solutions for both equations are obtained in the form of soliton, periodic, bright, and dark solitary wave solutions. There are many applications of the present traveling wave solutions in physics and furthermore, a wide class of coupled nonlinear evolution equations can be solved by this method. Keywords: Traveling wave solutions, Elliptic solutions, Generalized coupled Zakharov–Kuznetsov equation, Dispersive long wave equation, Modified extended direct algebraic method

Numerical simulation of seismic wave propagation from land-excited large volume air-gun source

Science.gov (United States)

Cao, W.; Zhang, W.

2017-12-01

The land-excited large volume air-gun source can be used to study regional underground structures and to detect temporal velocity changes. The air-gun source is characterized by rich low frequency energy (from bubble oscillation, 2-8Hz) and high repeatability. It can be excited in rivers, reservoirs or man-made pool. Numerical simulation of the seismic wave propagation from the air-gun source helps to understand the energy partitioning and characteristics of the waveform records at stations. However, the effective energy recorded at a distance station is from the process of bubble oscillation, which can not be approximated by a single point source. We propose a method to simulate the seismic wave propagation from the land-excited large volume air-gun source by finite difference method. The process can be divided into three parts: bubble oscillation and source coupling, solid-fluid coupling and the propagation in the solid medium. For the first part, the wavelet of the bubble oscillation can be simulated by bubble model. We use wave injection method combining the bubble wavelet with elastic wave equation to achieve the source coupling. Then, the solid-fluid boundary condition is implemented along the water bottom. And the last part is the seismic wave propagation in the solid medium, which can be readily implemented by the finite difference method. Our method can get accuracy waveform of land-excited large volume air-gun source. Based on the above forward modeling technology, we analysis the effect of the excited P wave and the energy of converted S wave due to different water shapes. We study two land-excited large volume air-gun fields, one is Binchuan in Yunnan, and the other is Hutubi in Xinjiang. The station in Binchuan, Yunnan is located in a large irregular reservoir, the waveform records have a clear S wave. Nevertheless, the station in Hutubi, Xinjiang is located in a small man-made pool, the waveform records have very weak S wave. Better understanding of

Radiation and propagation of electromagnetic waves

CERN Document Server

Tyras, George; Declaris, Nicholas

1969-01-01

Radiation and Propagation of Electromagnetic Waves serves as a text in electrical engineering or electrophysics. The book discusses the electromagnetic theory; plane electromagnetic waves in homogenous isotropic and anisotropic media; and plane electromagnetic waves in inhomogenous stratified media. The text also describes the spectral representation of elementary electromagnetic sources; the field of a dipole in a stratified medium; and radiation in anisotropic plasma. The properties and the procedures of Green's function method of solution, axial currents, as well as cylindrical boundaries a

Scattering of lower-hybrid waves by drift-wave density fluctuations: solutions of the radiative transfer equation

International Nuclear Information System (INIS)

Andrews, P.L.; Perkins, F.W.

1983-01-01

The investigation of the scattering of lower-hybrid waves by density fluctuations arising from drift waves in tokamaks is distinguished by the presence in the wave equation of a large, random, derivative-coupling term. The propagation of the lower-hybrid waves is well represented by a radiative transfer equation when the scale size of the density fluctuations is small compared to the overall plasma size. The radiative transfer equation is solved in two limits: first, the forward scattering limit, where the scale size of density fluctuations is large compared to the lower-hybrid perpendicular wavelength, and second, the large-angle scattering limit, where this inequality is reversed. The most important features of these solutions are well represented by analytical formulas derived by simple arguments. Based on conventional estimates for density fluctuations arising from drift waves and a parabolic density profile, the optical depth tau for scattering through a significant angle, is given by tauroughly-equal(2/N 2 /sub parallel/) (#betta#/sub p/i0/#betta#) 2 (m/sub e/c 2 /2T/sub i/)/sup 1/2/ [c/α(Ω/sub i/Ω/sub e/)/sup 1/2/ ], where #betta#/sub p/i0 is the central ion plasma frequency and T/sub i/ denotes the ion temperature near the edge of the plasma. Most of the scattering occurs near the surface. The transmission through the scattering region scales as tau - 1 and the emerging intensity has an angular spectrum proportional to cos theta, where sin theta = k/sub perpendicular/xB/sub p//(k/sub perpendicular/B/sub p/), and B/sub p/ is the poloidal field

Models for seismic wave propagation in periodically layered porous media

NARCIS (Netherlands)

Kudarova, A.; Van Dalen, K.N.; Drijkoningen, G.G.

2014-01-01

Several models are discussed for seismic wave propagation in periodically layered poroelastic media where layers represent mesoscopic-scale heterogeneities that are larger than the pore and grain sizes but smaller than the wavelength. The layers behave according to Biot’s theory. Wave propagation

On the propagation of truncated localized waves in dispersive silica

KAUST Repository

Salem, Mohamed

2010-01-01

Propagation characteristics of truncated Localized Waves propagating in dispersive silica and free space are numerically analyzed. It is shown that those characteristics are affected by the changes in the relation between the transverse spatial spectral components and the wave vector. Numerical experiments demonstrate that as the non-linearity of this relation gets stronger, the pulses propagating in silica become more immune to decay and distortion whereas the pulses propagating in free-space suffer from early decay and distortion. © 2010 Optical Society of America.

Topology optimization of vibration and wave propagation problems

DEFF Research Database (Denmark)

Jensen, Jakob Søndergaard

2007-01-01

The method of topology optimization is a versatile method to determine optimal material layouts in mechanical structures. The method relies on, in principle, unlimited design freedom that can be used to design materials, structures and devices with significantly improved performance and sometimes...... novel functionality. This paper addresses basic issues in simulation and topology design of vibration and wave propagation problems. Steady-state and transient wave propagation problems are addressed and application examples for both cases are presented....

The Liouville equation for flavour evolution of neutrinos and neutrino wave packets

Energy Technology Data Exchange (ETDEWEB)

Hansen, Rasmus Sloth Lundkvist; Smirnov, Alexei Yu., E-mail: rasmus@mpi-hd.mpg.de, E-mail: smirnov@mpi-hd.mpg.de [Max-Planck-Institut für Kernphysik, Saupfercheckweg 1, 69117 Heidelberg (Germany)

2016-12-01

We consider several aspects related to the form, derivation and applications of the Liouville equation (LE) for flavour evolution of neutrinos. To take into account the quantum nature of neutrinos we derive the evolution equation for the matrix of densities using wave packets instead of Wigner functions. The obtained equation differs from the standard LE by an additional term which is proportional to the difference of group velocities. We show that this term describes loss of the propagation coherence in the system. In absence of momentum changing collisions, the LE can be reduced to a single derivative equation over a trajectory coordinate. Additional time and spatial dependence may stem from initial (production) conditions. The transition from single neutrino evolution to the evolution of a neutrino gas is considered.

TSOS and TSOS-FK hybrid methods for modelling the propagation of seismic waves

Science.gov (United States)

Ma, Jian; Yang, Dinghui; Tong, Ping; Ma, Xiao

2018-05-01

We develop a new time-space optimized symplectic (TSOS) method for numerically solving elastic wave equations in heterogeneous isotropic media. We use the phase-preserving symplectic partitioned Runge-Kutta method to evaluate the time derivatives and optimized explicit finite-difference (FD) schemes to discretize the space derivatives. We introduce the averaged medium scheme into the TSOS method to further increase its capability of dealing with heterogeneous media and match the boundary-modified scheme for implementing free-surface boundary conditions and the auxiliary differential equation complex frequency-shifted perfectly matched layer (ADE CFS-PML) non-reflecting boundaries with the TSOS method. A comparison of the TSOS method with analytical solutions and standard FD schemes indicates that the waveform generated by the TSOS method is more similar to the analytic solution and has a smaller error than other FD methods, which illustrates the efficiency and accuracy of the TSOS method. Subsequently, we focus on the calculation of synthetic seismograms for teleseismic P- or S-waves entering and propagating in the local heterogeneous region of interest. To improve the computational efficiency, we successfully combine the TSOS method with the frequency-wavenumber (FK) method and apply the ADE CFS-PML to absorb the scattered waves caused by the regional heterogeneity. The TSOS-FK hybrid method is benchmarked against semi-analytical solutions provided by the FK method for a 1-D layered model. Several numerical experiments, including a vertical cross-section of the Chinese capital area crustal model, illustrate that the TSOS-FK hybrid method works well for modelling waves propagating in complex heterogeneous media and remains stable for long-time computation. These numerical examples also show that the TSOS-FK method can tackle the converted and scattered waves of the teleseismic plane waves caused by local heterogeneity. Thus, the TSOS and TSOS-FK methods proposed in

Benney's long wave equations

International Nuclear Information System (INIS)

Lebedev, D.R.

1979-01-01

Benney's equations of motion of incompressible nonviscous fluid with free surface in the approximation of long waves are analyzed. The connection between the Lie algebra of Hamilton plane vector fields and the Benney's momentum equations is shown

Study of the influence of semiconductor material parameters on acoustic wave propagation modes in GaSb/AlSb bi-layered structures by Legendre polynomial method

Energy Technology Data Exchange (ETDEWEB)

Othmani, Cherif, E-mail: othmanicheriffss@gmail.com; Takali, Farid; Njeh, Anouar; Ben Ghozlen, Mohamed Hédi

2016-09-01

The propagation of Rayleigh–Lamb waves in bi-layered structures is studied. For this purpose, an extension of the Legendre polynomial (LP) method is proposed to formulate the acoustic wave equation in the bi-layered structures induced by thin film Gallium Antimonide (GaSb) and with Aluminum Antimonide (AlSb) substrate in moderate thickness. Acoustic modes propagating along a bi-layer plate are shown to be quite different than classical Lamb modes, contrary to most of the multilayered structures. The validation of the LP method is illustrated by a comparison between the associated numerical results and those obtained using the ordinary differential equation (ODE) method. The convergency of the LP method is discussed through a numerical example. Moreover, the influences of thin film GaSb parameters on the characteristics Rayleigh–Lamb waves propagation has been studied in detail. Finally, the advantages of the Legendre polynomial (LP) method to analyze the multilayered structures are described. All the developments performed in this work were implemented in Matlab software.

Study of the influence of semiconductor material parameters on acoustic wave propagation modes in GaSb/AlSb bi-layered structures by Legendre polynomial method

International Nuclear Information System (INIS)

Othmani, Cherif; Takali, Farid; Njeh, Anouar; Ben Ghozlen, Mohamed Hédi

2016-01-01

The propagation of Rayleigh–Lamb waves in bi-layered structures is studied. For this purpose, an extension of the Legendre polynomial (LP) method is proposed to formulate the acoustic wave equation in the bi-layered structures induced by thin film Gallium Antimonide (GaSb) and with Aluminum Antimonide (AlSb) substrate in moderate thickness. Acoustic modes propagating along a bi-layer plate are shown to be quite different than classical Lamb modes, contrary to most of the multilayered structures. The validation of the LP method is illustrated by a comparison between the associated numerical results and those obtained using the ordinary differential equation (ODE) method. The convergency of the LP method is discussed through a numerical example. Moreover, the influences of thin film GaSb parameters on the characteristics Rayleigh–Lamb waves propagation has been studied in detail. Finally, the advantages of the Legendre polynomial (LP) method to analyze the multilayered structures are described. All the developments performed in this work were implemented in Matlab software.

APPARENT CROSS-FIELD SUPERSLOW PROPAGATION OF MAGNETOHYDRODYNAMIC WAVES IN SOLAR PLASMAS

Energy Technology Data Exchange (ETDEWEB)

Kaneko, T.; Yokoyama, T. [Department of Earth and Planetary Science, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo, 113-0033 (Japan); Goossens, M.; Doorsselaere, T. Van [Centre for Mathematical Plasma Astrophysics, Katholieke Universiteit Leuven, Celestijnenlaan 200B, Bus 2400, B-3001 Herverlee (Belgium); Soler, R.; Terradas, J. [Departament de Física, Universitat de les Illes Balears, E-07122 Palma de Mallorca (Spain); Wright, A. N., E-mail: kaneko@eps.s.u-tokyo.ac.jp [School of Mathematics and Statistics, University of St Andrews, St Andrews, KY16 9SS (United Kingdom)

2015-10-20

In this paper we show that the phase-mixing of continuum Alfvén waves and/or continuum slow waves in the magnetic structures of the solar atmosphere as, e.g., coronal arcades, can create the illusion of wave propagation across the magnetic field. This phenomenon could be erroneously interpreted as fast magnetosonic waves. The cross-field propagation due to the phase-mixing of continuum waves is apparent because there is no real propagation of energy across the magnetic surfaces. We investigate the continuous Alfvén and slow spectra in two-dimensional (2D) Cartesian equilibrium models with a purely poloidal magnetic field. We show that apparent superslow propagation across the magnetic surfaces in solar coronal structures is a consequence of the existence of continuum Alfvén waves and continuum slow waves that naturally live on those structures and phase-mix as time evolves. The apparent cross-field phase velocity is related to the spatial variation of the local Alfvén/slow frequency across the magnetic surfaces and is slower than the Alfvén/sound velocities for typical coronal conditions. Understanding the nature of the apparent cross-field propagation is important for the correct analysis of numerical simulations and the correct interpretation of observations.

Explicit solutions of the Camassa-Holm equation

International Nuclear Information System (INIS)

Parkes, E.J.; Vakhnenko, V.O.

2005-01-01

Explicit travelling-wave solutions of the Camassa-Holm equation are sought. The solutions are characterized by two parameters. For propagation in the positive x-direction, both periodic and solitary smooth-hump, peakon, cuspon and inverted-cuspon waves are found. For propagation in the negative x-direction, there are solutions which are just the mirror image in the x-axis of the aforementioned solutions. Some composite wave solutions of the Degasperis-Procesi equation are given in an appendix

Propagation of Electron Acoustic Soliton, Periodic and Shock Waves in Dissipative Plasma with a q-Nonextensive Electron Velocity Distribution

Science.gov (United States)

El-Hanbaly, A. M.; El-Shewy, E. K.; Elgarayhi, A.; Kassem, A. I.

2015-11-01

The nonlinear properties of small amplitude electron-acoustic (EA) solitary and shock waves in a homogeneous system of unmagnetized collisionless plasma with nonextensive distribution for hot electrons have been investigated. A reductive perturbation method used to obtain the Kadomstev-Petviashvili-Burgers equation. Bifurcation analysis has been discussed for non-dissipative system in the absence of Burgers term and reveals different classes of the traveling wave solutions. The obtained solutions are related to periodic and soliton waves and their behavior are shown graphically. In the presence of the Burgers term, the EXP-function method is used to solve the Kadomstev-Petviashvili-Burgers equation and the obtained solution is related to shock wave. The obtained results may be helpful in better conception of waves propagation in various space plasma environments as well as in inertial confinement fusion laboratory plasmas.

Acoustic wave propagation in fluids with coupled chemical reactions

International Nuclear Information System (INIS)

Margulies, T.S.; Schwarz, W.H.

1984-08-01

This investigation presents a hydroacoustic theory which accounts for sound absorption and dispersion in a multicomponent mixture of reacting fluids (assuming a set of first-order acoustic equations without diffusion) such that several coupled reactions can occur simultaneously. General results are obtained in the form of a biquadratic characteristic equation (called the Kirchhoff-Langevin equation) for the complex propagation variable chi = - (α + iω/c) in which α is the attenuation coefficient, c is the phase speed of the progressive wave and ω is the angular frequency. Computer simulations of sound absorption spectra have been made for three different chemical systems, each comprised of two-step chemical reactions using physico-chemical data available in the literature. The chemical systems studied include: (1) water-dioxane, (2) aqueous solutions of glycine and (3) cobalt polyphosphate mixtures. Explicit comparisons are made between the exact biquadratic characteristic solution and the approximate equation (sometimes referred to as a Debye equation) previously applied to interpret the experimental data for the chemical reaction contribution to the absorption versus frequency. The relative chemical reaction and classical viscothermal contributions to the sound absorption are also presented. Several discrepancies that can arise when estimating thermodynamic data (chemical reaction heats or volume changes) for multistep chemical reaction systems when making dilute solution or constant density assumptions are discussed

Nonlinear and diffraction effects in propagation of N-waves in randomly inhomogeneous moving media.

Science.gov (United States)

Averiyanov, Mikhail; Blanc-Benon, Philippe; Cleveland, Robin O; Khokhlova, Vera

2011-04-01

Finite amplitude acoustic wave propagation through atmospheric turbulence is modeled using a Khokhlov-Zabolotskaya-Kuznetsov (KZK)-type equation. The equation accounts for the combined effects of nonlinearity, diffraction, absorption, and vectorial inhomogeneities of the medium. A numerical algorithm is developed which uses a shock capturing scheme to reduce the number of temporal grid points. The inhomogeneous medium is modeled using random Fourier modes technique. Propagation of N-waves through the medium produces regions of focusing and defocusing that is consistent with geometrical ray theory. However, differences up to ten wavelengths are observed in the locations of fist foci. Nonlinear effects are shown to enhance local focusing, increase the maximum peak pressure (up to 60%), and decrease the shock rise time (about 30 times). Although the peak pressure increases and the rise time decreases in focal regions, statistical analysis across the entire wavefront at a distance 120 wavelengths from the source indicates that turbulence: decreases the mean time-of-flight by 15% of a pulse duration, decreases the mean peak pressure by 6%, and increases the mean rise time by almost 100%. The peak pressure and the arrival time are primarily governed by large scale inhomogeneities, while the rise time is also sensitive to small scales.

The effects of variable dust size and charge on dust acoustic waves propagating in a hybrid Cairns–Tsallis complex plasma

Science.gov (United States)

El-Taibany, W. F.; El-Siragy, N. M.; Behery, E. E.; Elbendary, A. A.; Taha, R. M.

2018-05-01

The propagation characteristics of dust acoustic waves (DAWs) in a dusty plasma consisting of variable size dust grains, hybrid Cairns-Tsallis-distributed electrons, and nonthermal ions are studied. The charging of the dust grains is described by the orbital-motion-limited theory and the size of the dust grains obeys the power law dust size distribution. To describe the nonlinear propagation of the DAWs, a Zakharov-Kuznetsov equation is derived using a reductive perturbation method. It is found that the nonthermal and nonextensive parameters influence the main properties of DAWs. Moreover, our results reveal that the rarefactive waves can propagate mainly in the proposed plasma model while compressive waves can be detected for a very small range of the distribution parameters of plasma species, and the DAWs are faster and wider for smaller size dust grains. Applications of the present results to dusty plasma observations are briefly discussed.

Electron-cyclotron wave propagation, absorption and current drive in the presence of neoclassical tearing modes

Science.gov (United States)

Isliker, Heinz; Chatziantonaki, Ioanna; Tsironis, Christos; Vlahos, Loukas

2012-09-01

We analyze the propagation of electron-cyclotron waves, their absorption and current drive when neoclassical tearing modes (NTMs), in the form of magnetic islands, are present in a tokamak plasma. So far, the analysis of the wave propagation and power deposition in the presence of NTMs has been performed mainly in the frame of an axisymmetric magnetic field, ignoring any effects from the island topology. Our analysis starts from an axisymmetric magnetic equilibrium, which is perturbed such as to exhibit magnetic islands. In this geometry, we compute the wave evolution with a ray-tracing code, focusing on the effect of the island topology on the efficiency of the absorption and current drive. To increase the precision in the calculation of the power deposition, the standard analytical flux-surface labeling for the island region has been adjusted from the usual cylindrical to toroidal geometry. The propagation up to the O-point is found to be little affected by the island topology, whereas the power absorbed and the driven current are significantly enhanced, because the resonant particles are bound to the small volumes in between the flux surfaces of the island. The consequences of these effects on the NTM evolution are investigated in terms of the modified Rutherford equation.

Electron-cyclotron wave propagation, absorption and current drive in the presence of neoclassical tearing modes

International Nuclear Information System (INIS)

Isliker, Heinz; Chatziantonaki, Ioanna; Tsironis, Christos; Vlahos, Loukas

2012-01-01

We analyze the propagation of electron-cyclotron waves, their absorption and current drive when neoclassical tearing modes (NTMs), in the form of magnetic islands, are present in a tokamak plasma. So far, the analysis of the wave propagation and power deposition in the presence of NTMs has been performed mainly in the frame of an axisymmetric magnetic field, ignoring any effects from the island topology. Our analysis starts from an axisymmetric magnetic equilibrium, which is perturbed such as to exhibit magnetic islands. In this geometry, we compute the wave evolution with a ray-tracing code, focusing on the effect of the island topology on the efficiency of the absorption and current drive. To increase the precision in the calculation of the power deposition, the standard analytical flux-surface labeling for the island region has been adjusted from the usual cylindrical to toroidal geometry. The propagation up to the O-point is found to be little affected by the island topology, whereas the power absorbed and the driven current are significantly enhanced, because the resonant particles are bound to the small volumes in between the flux surfaces of the island. The consequences of these effects on the NTM evolution are investigated in terms of the modified Rutherford equation. (paper)

Wave propagation in elastic medium with heterogeneous quadratic nonlinearity

International Nuclear Information System (INIS)

Tang Guangxin; Jacobs, Laurence J.; Qu Jianmin

2011-01-01

This paper studies the one-dimensional wave propagation in an elastic medium with spatially non-uniform quadratic nonlinearity. Two problems are solved analytically. One is for a time-harmonic wave propagating in a half-space where the displacement is prescribed on the surface of the half-space. It is found that spatial non-uniformity of the material nonlinearity causes backscattering of the second order harmonic, which when combined with the forward propagating waves generates a standing wave in steady-state wave motion. The second problem solved is the reflection from and transmission through a layer of finite thickness embedded in an otherwise linearly elastic medium of infinite extent, where it is assumed that the layer has a spatially non-uniform quadratic nonlinearity. The results show that the transmission coefficient for the second order harmonic is proportional to the spatial average of the nonlinearity across the thickness of the layer, independent of the spatial distribution of the nonlinearity. On the other hand, the coefficient of reflection is proportional to a weighted average of the nonlinearity across the layer thickness. The weight function in this weighted average is related to the propagating phase, thus making the coefficient of reflection dependent on the spatial distribution of the nonlinearity. Finally, the paper concludes with some discussions on how to use the reflected and transmitted second harmonic waves to evaluate the variance and autocorrelation length of nonlinear parameter β when the nonlinearity distribution in the layer is a stochastic process.

Spin-wave propagation spectrum in magnetization-modulated cylindrical nanowires

Energy Technology Data Exchange (ETDEWEB)

Li, Zhi-xiong; Wang, Meng-ning; Nie, Yao-zhuang; Wang, Dao-wei; Xia, Qing-lin [School of Physics and Electronics, Central South University, Changsha 410083 (China); Tang, Wei [School of Physics and Electronics, Central South University, Changsha 410083 (China); Suzhou Institute of Nano-tech and Nano-bionics, Chinese Academy of Sciences, Suzhou 215123 (China); Zeng, Zhong-ming [Suzhou Institute of Nano-tech and Nano-bionics, Chinese Academy of Sciences, Suzhou 215123 (China); Guo, Guang-hua, E-mail: guogh@mail.csu.edu.cn [School of Physics and Electronics, Central South University, Changsha 410083 (China)

2016-09-15

Spin-wave propagation in periodic magnetization-modulated cylindrical nanowires is studied by micromagnetic simulation. Spin wave scattering at the interface of two magnetization segments causes a spin-wave band structure, which can be effectively tuned by changing either the magnetization modulation level or the period of the cylindrical nanowire magnonic crystal. The bandgap width is oscillating with either the period or magnetization modulation due to the oscillating variation of the spin wave transmission coefficient through the interface of the two magnetization segments. Analytical calculation based on band theory is used to account for the micromagnetic simulation results. - Highlights: • A magnetization-modulated cylindrical nanowire magnonic crystal is proposed. • Propagating characteristics of spin waves in such magnonic crystal are studied. • Spin-wave spectra can be manipulated by changing modulation level and period.

Traveling waves and conservation laws for highly nonlinear wave equations modeling Hertz chains

Science.gov (United States)

Przedborski, Michelle; Anco, Stephen C.

2017-09-01

A highly nonlinear, fourth-order wave equation that models the continuum theory of long wavelength pulses in weakly compressed, homogeneous, discrete chains with a general power-law contact interaction is studied. For this wave equation, all solitary wave solutions and all nonlinear periodic wave solutions, along with all conservation laws, are derived. The solutions are explicitly parameterized in terms of the asymptotic value of the wave amplitude in the case of solitary waves and the peak of the wave amplitude in the case of nonlinear periodic waves. All cases in which the solution expressions can be stated in an explicit analytic form using elementary functions are worked out. In these cases, explicit expressions for the total energy and total momentum for all solutions are obtained as well. The derivation of the solutions uses the conservation laws combined with an energy analysis argument to reduce the wave equation directly to a separable first-order differential equation that determines the wave amplitude in terms of the traveling wave variable. This method can be applied more generally to other highly nonlinear wave equations.

Supersonic propagation of ionization waves in an underdense, laser-produced plasma

International Nuclear Information System (INIS)

Constantin, C.; Back, C.A.; Fournier, K.B.; Gregori, G.; Landen, O.L.; Glenzer, S.H.; Dewald, E.L.; Miller, M.C.

2005-01-01

A laser-driven supersonic ionization wave propagating through a millimeter-scale plasma of subcritical density up to 2-3 keV electron temperatures was observed. Propagation velocities initially ten times the sound speed were measured by means of time-resolved x-ray imaging diagnostics. The measured ionization wave trajectory is modeled analytically and by a two-dimensional radiation-hydrodynamics code. The comparison to the modeling suggests that nonlocal heat transport effects may contribute to the attenuation of the heat-wave propagation

Modeling stress wave propagation in rocks by distinct lattice spring model

Directory of Open Access Journals (Sweden)

Gaofeng Zhao

2014-08-01

Full Text Available In this paper, the ability of the distinct lattice spring model (DLSM for modeling stress wave propagation in rocks was fully investigated. The influence of particle size on simulation of different types of stress waves (e.g. one-dimensional (1D P-wave, 1D S-wave and two-dimensional (2D cylindrical wave was studied through comparing results predicted by the DLSM with different mesh ratios (lr and those obtained from the corresponding analytical solutions. Suggested values of lr were obtained for modeling these stress waves accurately. Moreover, the weak material layer method and virtual joint plane method were used to model P-wave and S-wave propagating through a single discontinuity. The results were compared with the classical analytical solutions, indicating that the virtual joint plane method can give better results and is recommended. Finally, some remarks of the DLSM on modeling of stress wave propagation in rocks were provided.

A problem-based approach to elastic wave propagation: the role of constraints

International Nuclear Information System (INIS)

Fazio, Claudio; Guastella, Ivan; Tarantino, Giovanni

2009-01-01

A problem-based approach to the teaching of mechanical wave propagation, focused on observation and measurement of wave properties in solids and on modelling of these properties, is presented. In particular, some experimental results, originally aimed at measuring the propagation speed of sound waves in metallic rods, are used in order to deepen the role of constraints in mechanical wave propagation. Interpretative models of the results obtained in the laboratory are built and implemented by using a well-known simulation environment. The simulation results are, then, compared with experimental data. The approach has been developed and experimented in the context of a workshop on mechanical wave propagation of the two-year Graduate Program for Physics Teacher Education at University of Palermo.

Three-dimensional stability of solitary kinetic Alfven waves and ion-acoustic waves

International Nuclear Information System (INIS)

Ghosh, G.; Das, K.P.

1994-01-01

Starting from a set of equations that lead to a linear dispersion relation coupling kinetic Alfven waves and ion-acoustic waves, three-dimensional KdV equations are derived for these waves. These equations are then used to investigate the three-dimensional stability of solitary kinetic Alfven waves and ion-acoustic waves by the small-k perturbation expansion method of Rowlands and Infeld. For kinetic Alfven waves it is found that there is instability if the direction of the plane-wave perturbation lies inside a cone, and the growth rate of the instability attains a maximum when the direction of the perturbation lies in the plane containing the external magnetic field and the direction of propagation of the solitary wave. For ion-acoustic waves the growth rate of instability attains a maximum when the direction of the perturbation lies in a plane perpendicular to the direction of propagation of the solitary wave. (Author)

Simulation of non-hydrostatic gravity wave propagation in the upper atmosphere

Directory of Open Access Journals (Sweden)

Y. Deng

2014-04-01

Full Text Available The high-frequency and small horizontal scale gravity waves may be reflected and ducted in non-hydrostatic simulations, but usually propagate vertically in hydrostatic models. To examine gravity wave propagation, a preliminary study has been conducted with a global ionosphere–thermosphere model (GITM, which is a non-hydrostatic general circulation model for the upper atmosphere. GITM has been run regionally with a horizontal resolution of 0.2° long × 0.2° lat to resolve the gravity wave with wavelength of 250 km. A cosine wave oscillation with amplitude of 30 m s−1 has been applied to the zonal wind at the low boundary, and both high-frequency and low-frequency waves have been tested. In the high-frequency case, the gravity wave stays below 200 km, which indicates that the wave is reflected or ducted in propagation. The results are consistent with the theoretical analysis from the dispersion relationship when the wavelength is larger than the cutoff wavelength for the non-hydrostatic situation. However, the low-frequency wave propagates to the high altitudes during the whole simulation period, and the amplitude increases with height. This study shows that the non-hydrostatic model successfully reproduces the high-frequency gravity wave dissipation.

Propagation of inertial-gravity waves on an island shelf

Science.gov (United States)

Bondur, V. G.; Sabinin, K. D.; Grebenyuk, Yu. V.

2015-09-01

The propagation of inertial-gravity waves (IGV) at the boundary of the Pacific shelf near the island of Oahu (Hawaii), whose generation was studied in the first part of this work [1], is analyzed. It is shown that a significant role there is played by the plane oblique waves; whose characteristics were identified by the method of estimating 3D wave parameters for the cases when the measurements are available only for two verticals. It is established that along with the descending propagation of energy that is typical of IGVs, wave packets ascend from the bottom to the upper layers, which is caused by the emission of waves from intense jets of discharged waters flowing out of a diffusor located at the bottom.

Generic short-time propagation of sharp-boundaries wave packets

Science.gov (United States)

Granot, E.; Marchewka, A.

2005-11-01

A general solution to the "shutter" problem is presented. The propagation of an arbitrary initially bounded wave function is investigated, and the general solution for any such function is formulated. It is shown that the exact solution can be written as an expression that depends only on the values of the function (and its derivatives) at the boundaries. In particular, it is shown that at short times (t << 2mx2/hbar, where x is the distance to the boundaries) the wave function propagation depends only on the wave function's values (or its derivatives) at the boundaries of the region. Finally, we generalize these findings to a non-singular wave function (i.e., for wave packets with finite-width boundaries) and suggest an experimental verification.

Spatial evolution equation of wind wave growth

Institute of Scientific and Technical Information of China (English)

WANG; Wei; (王; 伟); SUN; Fu; (孙; 孚); DAI; Dejun; (戴德君)

2003-01-01

Based on the dynamic essence of air-sea interactions, a feedback type of spatial evolution equation is suggested to match reasonably the growing process of wind waves. This simple equation involving the dominant factors of wind wave growth is able to explain the transfer of energy from high to low frequencies without introducing the concept of nonlinear wave-wave interactions, and the results agree well with observations. The rate of wave height growth derived in this dissertation is applicable to both laboratory and open sea, which solidifies the physical basis of using laboratory experiments to investigate the generation of wind waves. Thus the proposed spatial evolution equation provides a new approach for the research on dynamic mechanism of air-sea interactions and wind wave prediction.

Parametric Excitations of Fast Plasma Waves by Counter-propagating Laser Beams

International Nuclear Information System (INIS)

Shvets, G.; Fisch, N.J.

2001-01-01

Short- and long-wavelength plasma waves can become strongly coupled in the presence of two counter-propagating laser pump pulses detuned by twice the cold plasma frequency. What makes this four-wave interaction important is that the growth rate of the plasma waves occurs much faster than in the more obvious co-propagating geometry

Bifurcations of nonlinear ion acoustic travelling waves in the frame of a Zakharov-Kuznetsov equation in magnetized plasma with a kappa distributed electron

International Nuclear Information System (INIS)

Kumar Samanta, Utpal; Saha, Asit; Chatterjee, Prasanta

2013-01-01

Bifurcations of nonlinear propagation of ion acoustic waves (IAWs) in a magnetized plasma whose constituents are cold ions and kappa distributed electron are investigated using a two component plasma model. The standard reductive perturbation technique is used to derive the Zakharov-Kuznetsov (ZK) equation for IAWs. By using the bifurcation theory of planar dynamical systems to this ZK equation, the existence of solitary wave solutions and periodic travelling wave solutions is established. All exact explicit solutions of these travelling waves are determined. The results may have relevance in dense space plasmas

Counter-propagating wave interaction for contrast-enhanced ultrasound imaging

Science.gov (United States)

Renaud, G.; Bosch, J. G.; ten Kate, G. L.; Shamdasani, V.; Entrekin, R.; de Jong, N.; van der Steen, A. F. W.

2012-11-01

Most techniques for contrast-enhanced ultrasound imaging require linear propagation to detect nonlinear scattering of contrast agent microbubbles. Waveform distortion due to nonlinear propagation impairs their ability to distinguish microbubbles from tissue. As a result, tissue can be misclassified as microbubbles, and contrast agent concentration can be overestimated; therefore, these artifacts can significantly impair the quality of medical diagnoses. Contrary to biological tissue, lipid-coated gas microbubbles used as a contrast agent allow the interaction of two acoustic waves propagating in opposite directions (counter-propagation). Based on that principle, we describe a strategy to detect microbubbles that is free from nonlinear propagation artifacts. In vitro images were acquired with an ultrasound scanner in a phantom of tissue-mimicking material with a cavity containing a contrast agent. Unlike the default mode of the scanner using amplitude modulation to detect microbubbles, the pulse sequence exploiting counter-propagating wave interaction creates no pseudoenhancement behind the cavity in the contrast image.

Counter-propagating wave interaction for contrast-enhanced ultrasound imaging

International Nuclear Information System (INIS)

Renaud, G; Bosch, J G; Ten Kate, G L; De Jong, N; Van der Steen, A F W; Shamdasani, V; Entrekin, R

2012-01-01

Most techniques for contrast-enhanced ultrasound imaging require linear propagation to detect nonlinear scattering of contrast agent microbubbles. Waveform distortion due to nonlinear propagation impairs their ability to distinguish microbubbles from tissue. As a result, tissue can be misclassified as microbubbles, and contrast agent concentration can be overestimated; therefore, these artifacts can significantly impair the quality of medical diagnoses. Contrary to biological tissue, lipid-coated gas microbubbles used as a contrast agent allow the interaction of two acoustic waves propagating in opposite directions (counter-propagation). Based on that principle, we describe a strategy to detect microbubbles that is free from nonlinear propagation artifacts. In vitro images were acquired with an ultrasound scanner in a phantom of tissue-mimicking material with a cavity containing a contrast agent. Unlike the default mode of the scanner using amplitude modulation to detect microbubbles, the pulse sequence exploiting counter-propagating wave interaction creates no pseudoenhancement behind the cavity in the contrast image. (fast track communication)

Solitary-wave families of the Ostrovsky equation: An approach via reversible systems theory and normal forms

International Nuclear Information System (INIS)

Roy Choudhury, S.

2007-01-01

The Ostrovsky equation is an important canonical model for the unidirectional propagation of weakly nonlinear long surface and internal waves in a rotating, inviscid and incompressible fluid. Limited functional analytic results exist for the occurrence of one family of solitary-wave solutions of this equation, as well as their approach to the well-known solitons of the famous Korteweg-de Vries equation in the limit as the rotation becomes vanishingly small. Since solitary-wave solutions often play a central role in the long-time evolution of an initial disturbance, we consider such solutions here (via the normal form approach) within the framework of reversible systems theory. Besides confirming the existence of the known family of solitary waves and its reduction to the KdV limit, we find a second family of multihumped (or N-pulse) solutions, as well as a continuum of delocalized solitary waves (or homoclinics to small-amplitude periodic orbits). On isolated curves in the relevant parameter region, the delocalized waves reduce to genuine embedded solitons. The second and third families of solutions occur in regions of parameter space distinct from the known solitary-wave solutions and are thus entirely new. Directions for future work are also mentioned

Theory of wave propagation in partially saturated double-porosity rocks: a triple-layer patchy model

Science.gov (United States)

Sun, Weitao; Ba, Jing; Carcione, José M.

2016-04-01

Wave-induced local fluid flow is known as a key mechanism to explain the intrinsic wave dissipation in fluid-saturated rocks. Understanding the relationship between the acoustic properties of rocks and fluid patch distributions is important to interpret the observed seismic wave phenomena. A triple-layer patchy (TLP) model is proposed to describe the P-wave dissipation process in a double-porosity media saturated with two immiscible fluids. The double-porosity rock consists of a solid matrix with unique host porosity and inclusions which contain the second type of pores. Two immiscible fluids are considered in concentric spherical patches, where the inner pocket and the outer sphere are saturated with different fluids. The kinetic and dissipation energy functions of local fluid flow (LFF) in the inner pocket are formulated through oscillations in spherical coordinates. The wave propagation equations of the TLP model are based on Biot's theory and the corresponding Lagrangian equations. The P-wave dispersion and attenuation caused by the Biot friction mechanism and the local fluid flow (related to the pore structure and the fluid distribution) are obtained by a plane-wave analysis from the Christoffel equations. Numerical examples and laboratory measurements indicate that P-wave dispersion and attenuation are significantly influenced by the spatial distributions of both, the solid heterogeneity and the fluid saturation distribution. The TLP model is in reasonably good agreement with White's and Johnson's models. However, differences in phase velocity suggest that the heterogeneities associated with double-porosity and dual-fluid distribution should be taken into account when describing the P-wave dispersion and attenuation in partially saturated rocks.

A new auxiliary equation and exact travelling wave solutions of nonlinear equations

International Nuclear Information System (INIS)

Sirendaoreji

2006-01-01

A new auxiliary ordinary differential equation and its solutions are used for constructing exact travelling wave solutions of nonlinear partial differential equations in a unified way. The main idea of this method is to take full advantage of the auxiliary equation which has more new exact solutions. More new exact travelling wave solutions are obtained for the quadratic nonlinear Klein-Gordon equation, the combined KdV and mKdV equation, the sine-Gordon equation and the Whitham-Broer-Kaup equations

Full-vectorial propagation model and modified effective mode area of four-wave mixing in straight waveguides

DEFF Research Database (Denmark)

Guo, Kai; Friis, Søren Michael Mørk; Christensen, Jesper Bjerge

2017-01-01

We derive from Maxwell's equations full-vectorial nonlinear propagation equations of four-wave mixing valid in straight semiconductor-on-insulator waveguides. Special attention is given to the resulting effective mode area, which takes a convenient form known from studies in photonic crystal fibers......, but has not been introduced in the context of integrated waveguides. We show that the difference between our full-vectorial effective mode area and the scalar equivalent often referred to in the literature may lead to mistakes when evaluating the nonlinear refractive index and optimizing designs of new...

Full-wave Simulations of LH Wave Propagation in Toroidal Plasma with non-Maxwellian Electron Distributions

International Nuclear Information System (INIS)

Valeo, E.J.; Phillips, C.K.; Bonoli, P.T.; Wright, J.C.; Brambilla, M.

2007-01-01

The generation of energetic tails in the electron distribution function is intrinsic to lower-hybrid (LH) heating and current drive in weakly collisional magnetically confined plasma. The effects of these deformations on the RF deposition profile have previously been examined within the ray approximation. Recently, the calculation of full-wave propagation of LH waves in a thermal plasma has been accomplished using an adaptation of the TORIC code. Here, initial results are presented from TORIC simulations of LH propagation in a toroidal plasma with non-thermal electrons. The required efficient computation of the hot plasma dielectric tensor is accomplished using a technique previously demonstrated in full-wave simulations of ICRF propagation in plasma with non-thermal ions

On a functional equation related to the intermediate long wave equation

International Nuclear Information System (INIS)

Hone, A N W; Novikov, V S

2004-01-01

We resolve an open problem stated by Ablowitz et al (1982 J. Phys. A: Math. Gen. 15 781) concerning the integral operator appearing in the intermediate long wave equation. We explain how this is resolved using the perturbative symmetry approach introduced by one of us with Mikhailov. By solving a certain functional equation, we prove that the intermediate long wave equation and the Benjamin-Ono equation are the unique integrable cases within a particular class of integro-differential equations. Furthermore, we explain how the perturbative symmetry approach is naturally extended to treat equations on a periodic domain. (letter to the editor)

Guided propagation of Alfven waves in a toroidal plasma

International Nuclear Information System (INIS)

Borg, G.G.; Brennan, M.H.; Cross, R.C.; Giannone, L.; Donnelly, I.J.

1985-01-01

Experimental results are presented which show that the Alfven wave is strongly guided by magnetic fields. The experiment was conducted in a Tokamak plasma using a small dipole loop antenna to generate a localised Alfven ray. The ray was observed, with magnetic probes, to propagate as a localised disturbance along the curved lines of the steady magnetic field without significant refraction due to the effects of finite frequency, resistivity or magnetic field gradients. These results agree with theoretical predictions and demonstrate that a localised Alfven wave may be excited, and may propagate, independently of the fast wave, as expected. The implication of these results for the Alfven wave heating scheme is discussed. (author)

Guided propagation of Alfven waves in a toroidal plasma

Energy Technology Data Exchange (ETDEWEB)

Borg, G G; Brennan, M H; Cross, R C; Giannone, L.; Donnelly, I J

1985-10-01

Experimental results are presented which show that the Alfven wave is strongly guided by magnetic fields. The experiment was conducted in a Tokamak plasma using a small dipole loop antenna to generate a localised Alfven ray. The ray was observed, with magnetic probes, to propagate as a localised disturbance along the curved lines of the steady magnetic field without significant refraction due to the effects of finite frequency, resistivity or magnetic field gradients. These results agree with theoretical predictions and demonstrate that a localised Alfven wave may be excited, and may propagate, independently of the fast wave, as expected. The implication of these results for the Alfven wave heating scheme is discussed.

24 GHz cmWave Radio Propagation Through Vegetation

DEFF Research Database (Denmark)

Rodriguez, Ignacio; Abreu, Renato Barbosa; Portela Lopes de Almeida, Erika

2016-01-01

This paper presents a measurement-based analysis of cm-wave radio propagation through vegetation at 24 GHz. A set of dedicated directional measurements were performed with horn antennas located close to street level inside a densely-vegetated area illuminated from above. The full azimuth was exam......This paper presents a measurement-based analysis of cm-wave radio propagation through vegetation at 24 GHz. A set of dedicated directional measurements were performed with horn antennas located close to street level inside a densely-vegetated area illuminated from above. The full azimuth...

Discrete Element Simulation of Elastoplastic Shock Wave Propagation in Spherical Particles

Directory of Open Access Journals (Sweden)

M. Shoaib

2011-01-01

Full Text Available Elastoplastic shock wave propagation in a one-dimensional assembly of spherical metal particles is presented by extending well-established quasistatic compaction models. The compaction process is modeled by a discrete element method while using elastic and plastic loading, elastic unloading, and adhesion at contacts with typical dynamic loading parameters. Of particular interest is to study the development of the elastoplastic shock wave, its propagation, and reflection during entire loading process. Simulation results yield information on contact behavior, velocity, and deformation of particles during dynamic loading. Effects of shock wave propagation on loading parameters are also discussed. The elastoplastic shock propagation in granular material has many practical applications including the high-velocity compaction of particulate material.

Solitary wave and periodic wave solutions for the thermally forced gravity waves in atmosphere

International Nuclear Information System (INIS)

Li Ziliang

2008-01-01

By introducing a new transformation, a new direct and unified algebraic method for constructing multiple travelling wave solutions of general nonlinear evolution equations is presented and implemented in a computer algebraic system, which extends Fan's direct algebraic method to the case when r > 4. The solutions of a first-order nonlinear ordinary differential equation with a higher degree nonlinear term and Fan's direct algebraic method of obtaining exact solutions to nonlinear partial differential equations are applied to the combined KdV-mKdV-GKdV equation, which is derived from a simple incompressible non-hydrostatic Boussinesq equation with the influence of thermal forcing and is applied to investigate internal gravity waves in the atmosphere. As a result, by taking advantage of the new first-order nonlinear ordinary differential equation with a fifth-degree nonlinear term and an eighth-degree nonlinear term, periodic wave solutions associated with the Jacobin elliptic function and the bell and kink profile solitary wave solutions are obtained under the effect of thermal forcing. Most importantly, the mechanism of propagation and generation of the periodic waves and the solitary waves is analysed in detail according to the values of the heating parameter, which show that the effect of heating in atmosphere helps to excite westerly or easterly propagating periodic internal gravity waves and internal solitary waves in atmosphere, which are affected by the local excitation structures in atmosphere. In addition, as an illustrative sample, the properties of the solitary wave solution and Jacobin periodic solution are shown by some figures under the consideration of heating interaction

Analytical solution for the transient wave propagation of a buried cylindrical P-wave line source in a semi-infinite elastic medium with a fluid surface layer

Science.gov (United States)

Shan, Zhendong; Ling, Daosheng

2018-02-01

This article develops an analytical solution for the transient wave propagation of a cylindrical P-wave line source in a semi-infinite elastic solid with a fluid layer. The analytical solution is presented in a simple closed form in which each term represents a transient physical wave. The Scholte equation is derived, through which the Scholte wave velocity can be determined. The Scholte wave is the wave that propagates along the interface between the fluid and solid. To develop the analytical solution, the wave fields in the fluid and solid are defined, their analytical solutions in the Laplace domain are derived using the boundary and interface conditions, and the solutions are then decomposed into series form according to the power series expansion method. Each item of the series solution has a clear physical meaning and represents a transient wave path. Finally, by applying Cagniard's method and the convolution theorem, the analytical solutions are transformed into the time domain. Numerical examples are provided to illustrate some interesting features in the fluid layer, the interface and the semi-infinite solid. When the P-wave velocity in the fluid is higher than that in the solid, two head waves in the solid, one head wave in the fluid and a Scholte wave at the interface are observed for the cylindrical P-wave line source.

Electron thermal conductivity from heat wave propagation in Wendelstein 7-AS