Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userAnalytic and numerical solutions to the seismic wave equation in continuous media
106
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
This paper presents two approaches to mathematical modelling of a synthetic seismic pulse, and a comparison between them. First, a new analytical model is developed in two-dimensional Cartesian coordinates. Combined with an initial condition of sufficient symmetry, this provides a valuable check for the validity of the numerical method that follows. A particular initial condition is found which allows for a new closed-form solution. A numerical scheme is then presented which combines a spectral (Fourier) representation for displacement components and wave-speed parameters, a fourth order Runge-Kutta integration method, and an absorbing boundary layer. The resulting large system of differential equations is solved in parallel on suitable enhanced performance desktop hardware in a new software implementation. This provides an alternative approach to forward modelling of waves within isotropic media which is efficient, and tailored to rapid and flexible developments in modelling seismic structure, for example, shallow depth environmental applications. Visual comparisons of the analytic solution and the numerical scheme are presented.
rate research
Read More
Seismology was developed on Earth and shaped our model of the Earths interior over the 20th century. With the exception of the Philae lander, all in situ extraterrestrial seismological effort to date was limited to other terrestrial planets. All have in common a rigid crust above a solid mantle. The coming years may see the installation of seismometers on Europa, Titan and Enceladus, so it is necessary to adapt seismological concepts to the setting of worlds with global oceans covered in ice. Here we use waveform analyses to identify and classify wave types, developing a lexicon for icy ocean world seismology intended to be useful to both seismologists and planetary scientists. We use results from spectral-element simulations of broadband seismic wavefields to adapt seismological concepts to icy ocean worlds. We present a concise naming scheme for seismic waves and an overview of the features of the seismic wavefield on Europa, Titan, Ganymede and Enceladus. In close connection with geophysical interior models, we analyze simulated seismic measurements of Europa and Titan that might be used to constrain geochemical parameters governing the habitability of a sub-ice ocean.
Seismic wave propagation forms the basis for most aspects of seismological research, yet solving the wave equation is a major computational burden that inhibits the progress of research. This is exaspirated by the fact that new simulations must be performed when the velocity structure or source location is perturbed. Here, we explore a prototype framework for learning general solutions using a recently developed machine learning paradigm called Neural Operator. A trained Neural Operator can compute a solution in negligible time for any velocity structure or source location. We develop a scheme to train Neural Operators on an ensemble of simulations performed with random velocity models and source locations. As Neural Operators are grid-free, it is possible to evaluate solutions on higher resolution velocity models than trained on, providing additional computational efficiency. We illustrate the method with the 2D acoustic wave equation and demonstrate the methods applicability to seismic tomography, using reverse mode automatic differentiation to compute gradients of the wavefield with respect to the velocity structure. The developed procedure is nearly an order of magnitude faster than using conventional numerical methods for full waveform inversion.
We obtain time dependent $q$-Gaussian wave-packet solutions to a non linear Schrodinger equation recently advanced by Nobre, Rego-Montero and Tsallis (NRT) [Phys. Rev. Lett. 106 (2011) 10601]. The NRT non-linear equation admits plane wave-like solutions ($q$-plane waves) compatible with the celebrated de Broglie relations connecting wave number and frequency, respectively, with energy and momentum. The NRT equation, inspired in the $q$-generalized thermostatistical formalism, is characterized by a parameter $q$, and in the limit $q to 1$ reduces to the standard, linear Schrodinger equation. The $q$-Gaussian solutions to the NRT equation investigated here admit as a particular instance the previously known $q$-plane wave solutions. The present work thus extends the range of possible processes yielded by the NRT dynamics that admit an analytical, exact treatment. In the $q to 1$ limit the $q$-Gaussian solutions correspond to the Gaussian wave packet solutions to the free particle linear Schrodinger equation. In the present work we also show that there are other families of nonlinear Schrodinger-like equations, besides the NRT one, exhibiting a dynamics compatible with the de Broglie relations. Remarkably, however, the existence of time dependent Gaussian-like wave packet solutions is a unique feature of the NRT equation not shared by the aforementioned, more general, families of nonlinear evolution equations.
Crack microgeometries pose a paramount influence on effective elastic characteristics and sonic responses. Geophysical exploration based on seismic methods are widely used to assess and understand the presence of fractures. Numerical simulation as a promising way for this issue, still faces some challenges. With the rapid development of computers and computational techniques, discrete-based numerical approaches with desirable properties have been increasingly developed, but have not yet extensively applied to seismic response simulation for complex fractured media. For this purpose, we apply the coupled LSM-DFN model (Liu and Fu, 2020b) to examining the validity in emulating elastic wave propagation and scattering in naturally-fractured media. By comparing to the theoretical values, the implement of the schema is validated with input parameters optimization. Moreover, dynamic elastic moduli from seismic responses are calculated and compared with static ones from quasi-static loading of uniaxial compression tests. Numerical results are consistent with the tendency of theoretical predictions and available experimental data. It shows the potential for reproducing the seismic responses in complex fractured media and quantitatively investigating the correlations and differences between static and dynamic elastic moduli.
We consider the radial wave equation in similarity coordinates within the semigroup formalism. It is known that the generator of the semigroup exhibits a continuum of eigenvalues and embedded in this continuum there exists a discrete set of eigenvalues with analytic eigenfunctions. Our results show that, for sufficiently regular data, the long time behaviour of the solution is governed by the analytic eigenfunctions. The same techniques are applied to the linear stability problem for the fundamental self--similar solution $chi_T$ of the wave equation with a focusing power nonlinearity. Analogous to the free wave equation, we show that the long time behaviour (in similarity coordinates) of linear perturbations around $chi_T$ is governed by analytic mode solutions. In particular, this yields a rigorous proof for the linear stability of $chi_T$ with the sharp decay rate for the perturbations.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
