3D-solitons in a strongly magnetized plasma: Transition from flat to spherical waves of the Zakharov-Kuznetsov equation

The Zakharov-Kuznetsov equation has been used to describe ion-acoustic wave propagation in a strong magnetic plasma. An initial-value problem has been solved for this equation on the basis of a numerical method that uses the fast Fourier transform technique for calculating space derivatives and a fourth order Runge-Kutta method for the time scheme. Numerical simulations have shown that the disturbed flat solitary waves can break up into spherical ones.     

An unstaggered constrained transport method for the 3D ideal magnetohydrodynamic equations

Numerical methods for solving the ideal magnetohydrodynamic (MHD) equations in more than one space dimension must either confront the challenge of controlling errors in the discrete divergence of the magnetic field, or else be faced with nonlinear numerical instabilities. One approach for controlling the discrete divergence is through a so-called constrained transport method, which is based on first predicting a magnetic field through a standard finite volume solver, and then correcting this field through the appropriate use of a magnetic vector potential. In this work we develop a constrained transport method for the 3D ideal MHD equations that is based on a high-resolution wave propagation scheme. Our proposed scheme is the 3D extension of the 2D scheme developed by Rossmanith [J.A. Rossmanith, An unstaggered, high-resolution constrained transport method for magnetohydrodynamic flows, SIAM J. Sci. Comput. 28 (2006) 1766], and is based on the high-resolution wave propagation method of Langseth and LeVeque [J.O. Langseth, R.J. LeVeque, A wave propagation method for threedimensional hyperbolic conservation laws, J. Comput. Phys. 165 (2000) 126]. In particular, in our extension we take great care to maintain the three most important properties of the 2D scheme: (1) all quantities, including all components of the magnetic field and magnetic potential, are treated as cell-centered; (2) we develop a high-resolution wave propagation scheme for evolving the magnetic potential; and (3) we develop a wave limiting approach that is applied during the vector potential evolution, which controls unphysical oscillations in the magnetic field. One of the key numerical difficulties that is novel to 3D is that the transport equation that must be solved for the magnetic vector potential is only weakly hyperbolic. In presenting our numerical algorithm we describe how to numerically handle this problem of weak hyperbolicity, as well as how to choose an appropriate gauge condition. The resulting scheme is applied to several numerical test cases.     

A method of solving the stiffness problem in Biot＇s poroelastic equations using a staggered high-order finite-difference

In modelling elastic wave propagation in a porous medium, when the ratio between the fluid viscosity and the medium permeability is comparatively large, the stiffness problem of Blot＇s poroelastic equations will be encountered. In the paper, a partition method is developed to solve the stiffness problem with a staggered high-order finite-difference. The method splits the Biot equations into two systems. One is stiff, and solved analytically, the other is nonstiff, and solved numerically by using a high-order staggered-grid finite-difference scheme. The time step is determined by the staggered finite-difference algorithm in solving the nonstiff equations, thus a coarse time step may be employed. Therefore, the computation efficiency and computational stability are improved greatly. Also a perfect by matched layer technology is used in the split method as absorbing boundary conditions. The numerical results are compared with the analytical results and those obtained from the conventional staggered-grid finite-difference method in a homogeneous model, respectively. They are in good agreement with each other. Finally, a slightly more complex model is investigated and compared with related equivalent model to illustrate the good performance of the staggered-grid finite-difference scheme in the partition method.     

From the paddle to the beach – A Boussinesq shallow water numerical wave tank based on Madsen and Sørensen’s equations

This article describes a one-dimensional numerical model of a shallow-water flume with an in-built piston paddle moving boundary wavemaker. The model is based on a set of enhanced Boussinesq equations and the nonlinear shallow water equations. Wave breaking is described approximately, by locally switching to the nonlinear shallow water equations when a critical wave steepness is reached. The moving shoreline is calculated as part of the solution. The piston paddle wavemaker operates on a movable grid, which is Lagrangian on the paddle face and Eulerian away from the paddle. The governing equations are, however, evolved on a fixed mapped grid, and the newly calculated solution is transformed back onto the moving grid via a domain mapping technique. Validation test results are compared against analytical solutions, confirming correct discretisation of the governing equations, wave generation via the numerical paddle, and movement of the wet/dry front. Simulations are presented that reproduce laboratory experiments of wave runup on a plane beach and wave overtopping of a laboratory seawall, involving solitary waves and compact wave groups. In practice, the numerical model is suitable for simulating the propagation of weakly dispersive waves and can additionally model any associated inundation, overtopping or inland flooding within the same simulation.     

Plasmons and polaritons in a semi-infinite plasma and a plasma slab

Plasmon and polariton modes are derived for an ideal semi-infinite (half-space) plasma and an ideal plasma slab by using a general, unifying procedure, based on equations of motion, Maxwell's equations and suitable boundary conditions. Known results are re-obtained in much a more direct manner and new ones are derived. The approach consists of representing the charge disturbances by a displacement field in the positions of the moving particles (electrons). The dielectric response and the electron energy loss are computed. The surface contribution to the energy loss exhibits an oscillatory behaviour in the transient regime near the surfaces. The propagation of an electromagnetic wave in these plasmas is treated by using the retarded electromagnetic potentials. The resulting integral equations are solved and the reflected and refracted waves are computed, as well as the reflection coefficient. For the slab we compute also the transmitted wave and the transmission coefficient. Generalized Fresnel's relations are thereby obtained for any incidence angle and polarization. Bulk and surface plasmon-polariton modes are identified. As it is well known, the field inside the plasma is either damped (evanescent) or propagating (transparency regime), and the reflection coefficient for a semi-infinite plasma exhibits an abrupt enhancement on passing from the propagating regime to the damped one (total reflection). Similarly, apart from characteristic oscillations, the reflection and transmission coefficients for a plasma slab exhibit an appreciable enhancement in the damped regime.     

Exact boundary condition for time-dependent wave equation based on boundary integral

An exact non-reflecting boundary conditions based on a boundary integral equation or a modified Kirchhoff-type formula is derived for exterior three-dimensional wave equations. The Kirchhoff-type non-reflecting boundary condition is originally proposed by L. Ting and M.J. Miksis [J. Acoust. Soc. Am. 80 (1986) 1825] and numerically tested by D. Givoli and D. Cohen [J. Comput. Phys. 117 (1995) 102] for a spherically symmetric problem. The computational advantage of Ting–Miksis boundary condition is that its temporal non-locality is limited to a fixed amount of past information. However, a long-time instability is exhibited in testing numerical solutions by using a standard non-dissipative finite-difference scheme. The main purpose of this work is to present a new exact boundary condition and to eliminate the long-time instability. The proposed exact boundary condition can be considered as a limit case of Ting–Miksis boundary condition when the two artificial boundaries used in their method approach each other. Our boundary condition is actually a boundary integral equation on a single artificial boundary for wave equations, which is to be solved in conjunction with the interior wave equation. The new boundary condition needs only one artificial boundary, which can be of any shape, i.e., sphere, cubic surface, etc. It keeps all merits of the original Kirchhoff boundary condition such as restricting the temporal non-locality, free of numerical evaluation of any special functions and so on. Numerical approximation to the artificial boundary condition on cubic surface is derived and three-dimensional numerical tests are carried out on the cubic computational domain.     

Electron exchange-correlation effects on dispersion properties of electrostatic waves in degenerate quantum plasmas

Based on the quantum Magnetohydrodynamic (QMHD) model, the obliquely propagation of electrostatic waves in degenerate magnetized quantum plasmas with electron exchange-correlation effects are theoretically investigated. The modified linear dispersion relations of electrostatic waves are obtained and discussed in some specific cases. The analytical results clearly show that the dispersion properties of the high frequency electron waves (including the Langmuir wave and upper-hybrid wave) and the low frequency ion acoustic wave are modified by the quantum effects together with the electron exchange-correlation effects. The numerical results depict that the Langmuir wave and upper-hybrid wave can be unstable in the presence of the electron exchange-correlation effects, and it is also evidently indicated that the electron exchange-correlation effects can reduce the phase velocity of the waves, especially in the high wave number region. The corresponding results should be of relevance for identifying electrostatic fluctuations which transport in an inhomogeneous and magnetized quantum plasmas.     

Wide-angle one-way wave equations

A one-way wave equation, also known as a paraxial or parabolic wave equation, is a differential equation that permits wave propagation in certain directions only. Such equations are used regularly in underwater acoustics, in geophysics, and as energy-absorbing numerical boundary conditions. The design of a one-way wave equation is connected with the approximation of (1-s2)1/2 on [-1,1] by a rational function, which has usually been carried out by Padé approximation. This article presents coefficients for L2, L infinity, and other alternative classes of approximants that have better wide-angle behavior. For theoretical results establishing the well posedness of these wide-angle equations, see the work of Trefethen and Halpern ["Well-posedness of one-way wave equations and absorbing boundary conditions," Math. Comput. 47, 421-435 (1986)].     

Poroelastic modeling of seismic boundary conditions across a fracture

Permeability of a fracture can affect how the fracture interacts with seismic waves. To examine this effect, a simple mathematical model that describes the poroelastic nature of wave-fracture interaction is useful. In this paper, a set of boundary conditions is presented which relate wave-induced particle velocity (or displacement) and stress including fluid pressure across a compliant, fluid-bearing fracture. These conditions are derived by modeling a fracture as a thin porous layer with increased compliance and finite permeability. Assuming a small layer thickness, the boundary conditions can be derived by integrating the governing equations of poroelastic wave propagation. A finite jump in the stress and velocity across a fracture is expressed as a function of the stress and velocity at the boundaries. Further simplification for a thin fracture yields a set of characteristic parameters that control the seismic response of single fractures with a wide range of mechanical and hydraulic properties. These boundary conditions have potential applications in simplifying numerical models such as finite-difference and finite-element methods to compute seismic wave scattering off nonplanar (e.g., curved and intersecting) fractures.     

Well-posedness of continuum models for weakly ionized plasmas

Continuum fluid models of weakly ionized plasmas are useful in the design and control of plasma-assisted deposition and etching processes. The equations in these models are numerically stiff. Their stiffness is affected by the imposed boundary conditions. In this work, a DC discharge model is studied and the effect of the boundary conditions on the model solution is investigated. It is established, both analytically and numerically, that depending on the choice of boundary conditions the model may range from being ill-posed to being solvable with standard software. It is also established that excessive truncation error maybe present in numerical simulations which appear to qualitatively capture plasma structure. Accurate numerical simulations of the considered model, with alternate boundary conditions, are shown to capture many characteristics of a DC discharge, albeit at lower values of applied voltage than those reported in the literature. Finally, model shortcomings are discussed     

The contribution of a lateral wave in simulating low-frequency sound fields in an irregular waveguide with a liquid bottom

A further development of a previously proposed approach to calculating the sound field in an arbitrarily irregular ocean is presented. The approach is based on solving the first-order causal mode equations, which are equivalent to the boundary-value problem for acoustic wave equations in terms of the cross-section method. For the mode functions depending on the horizontal coordinate, additional terms are introduced in the cross-section equations to allow for the multilayer structure of the medium. A numerical solution to the causal equations is sought using the fundamental matrix equation. For the modes of the discrete spectrum and two fixed low frequencies, calculations are performed for an irregular two-layer waveguide model with fluid sediments, which is close to the actual conditions of low-frequency sound propagation in the coastal zone of the oceanic shelf. The calculated propagation loss curves are used as an example for comparison with results that can be obtained for the given waveguide model with the use of adiabatic and one-way propagation approximations.     

Numerical modeling of urban sound fields by a diffusion process

In this paper, we present the numerical implementation of a sound field model used in urban acoustics. The mathematical model being based on a classic diffusion equation for the sound energy, a simple finite difference scheme is applied. We give also some finite difference equations for simple boundary conditions, like absorption by a wall and at building edges. The two-dimensional numerical scheme is then compared to analytical solutions of the sound field propagation in a rectangular street with a good agreement, both in the steady state and in the time varying state. Finally it is suggested that the adjustment of usual softwares for heat transfer could be an interesting and low cost way to develop powerful acoustic softwares for the prediction of noise in urban areas.     

Sound field prediction of ultrasonic lithotripsy in water with spheroidal beam equations

With converged shock wave,extracorporeal shock wave lithotripsy(ESWL)has become a preferable way to crush human calculi because of its advantages of efficiency and non-intrusion.Nonlinear spheroidal beam equations(SBE)are employed to illustrate the acoustic wave propagation for transducers with a wide aperture angle.To predict the acoustic field distribution precisely,boundary conditions are obtained for the SBE model of the monochromatic wave when the source is located on the focus of an ESWL transducer.Numerical results of the monochromatic wave propagation in water are analyzed and the influences of half-angle,fundamental frequency,and initial pressure are investigated.According to our results,with optimization of these factors,the pressure focal gain of ESWL can be enhanced and the effectiveness of treatment can be improved.     

基于广义交替数值通量的局部间断Galerkin方法求解二维波动方程

波的传播往往在复杂的地质结构中进行,如何有效地求解非均匀介质中的波动方程一直是研究的热点.本文将局部间断Galekin(local discontinuous Galerkin, LDG)方法引入到数值求解波动方程中.首先引入辅助变量,将二阶波动方程写成一阶偏微分方程组,然后对相应的线性化波动方程和伴随方程构造间断Galerkin格式;为了保证离散格式满足能量守恒,在单元边界上选取广义交替数值通量,理论证明该方法满足能量守恒性.在时间离散上,采用指数积分因子方法,为了提高计算效率,应用Krylov子空间方法近似指数矩阵与向量的乘积.数值实验中给出了带有精确解的算例,验证了LDG方法的数值精度和能量守恒性;此外,也考虑了非均匀介质和复杂计算区域的计算,结果表明LDG方法适合模拟具有复杂结构和多尺度结构介质中的传播.     

An iterative method to solve a regularized model for strongly nonlinear long internal waves

We present a simple iterative scheme to solve numerically a regularized internal wave model describing the large amplitude motion of the interface between two layers of different densities. Compared with the original strongly nonlinear internal wave model of Miyata [10] and Choi and Camassa [2], the regularized model adopted here suppresses shear instability associated with a velocity jump across the interface, but the coupling between the upper and lower layers is more complicated so that an additional system of coupled linear equations must be solved at every time step after a set of nonlinear evolution equations are integrated in time. Therefore, an efficient numerical scheme is desirable. In our iterative scheme, the linear system is decoupled and simple linear operators with constant coefficients are required to be inverted. Through linear analysis, it is shown that the scheme converges fast with an optimum choice of iteration parameters. After demonstrating its effectiveness for a model problem, the iterative scheme is applied to solve the regularized internal wave model using a pseudo-spectral method for the propagation of a single internal solitary wave and the head-on collision between two solitary waves of different wave amplitudes.     

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

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.     

Numerical simulation of nonlinear propagation of sound waves in a finite horn

Nonlinear acoustic propagation generated by a piston in a finite horn is numerically studied.A quasi-one-dimensional nonlinear model with varying cross-section uses high-order low-dispersion numerical schemes to solve the governing equation.Because of the nonlinear wave distortion and reflected sound waves at the mouth,broadband time-domain impedance boundary conditions are employed.The impedance approximation can be optimized to identify the complex-conjugate pole-residue pairs of the impedance functions,which can be calculated by fast and efficient recursive convolution.The numerical results agree very well with experimental data in the situations of weak nonlinear wave propagation in an exponential horn,it is shown that the model can describe the broadband characteristics caused by nonlinear distortion.Moreover,finite-amplitude acoustic propagation in types of horns is simulated,including hyperbolic,conical,exponential and sinusoidal horns.It is found that sound pressure levels at the horn mouth are strongly affected by the horn profiles,the driving velocity and frequency of the piston.The paper also discusses the influence of the horn geometry on nonlinear waveform distortion.     

A study of propagation of nonlinear acoustic-gravity waves in the middle and upper atmosphere by means of numerical modeling

numerical model of the vertical propagation and decay of nonlinear acoustic-gravity waves (AGW) from the Earth surface to the upper atmosphere is described. Monochromatic vertical velocity variations at the Earth surface are used as the AGW source in the model. The numerical method for solving three-dimensional hydrodynamic equations is based on finite-difference representation of the fundamental laws of conservation, which makes it possible to calculate not only smooth, but also physically correct generalized solutions of the hydrodynamic equations. The equations are solved in a range of altitudes from the ground up to 500 km. The background temperature, density, molecular viscosity and thermal conductivity coefficient are specified according to standard atmosphere models. The dependence of the characteristics of the waves on the amplitude of the wave source at the lower boundary is examined. The amplitudes of the AGW increase with the altitude, and the waves can break down due to nonlinear effects in the middle and upper atmosphere, depending on the amplitude of the source.     

A numerical method for computation of sound radiation from an unflanged duct

Some issues involved in establishing a numerical model for sound radiation from a straight duct are addressed in this paper. The main ingredient of the numerical method is solutions of linearized Euler equations using a high order compact scheme. Farfield directivity is estimated through an integral solution of Ffowcs-Williams Hawkings equations. A generic test case of planar wave radiation from an unflanged duct is studied. The sound pressure level and wave propagation in the nearfield are analyzed, together with the farfield directivity. Comparison with analytical solutions shows good agreement. The effect of grid resolution on the sound radiation pattern and the construction of integration surface on the farfield directivity are discussed.