Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userAnalysis of dynamic ruptures generating seismic waves in a self-gravitating planet: an iterative coupling scheme and well-posedness
97
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We study the solution of the system of equations describing the dynamical evolution of spontaneous ruptures generated in a prestressed elastic-gravitational deforming body and governed by rate and state friction laws. We propose an iterative coupling scheme based on a weak formulation with nonlinear interior boundary conditions, both for continuous time and with implicit discretization (backward Euler) in time. We regularize the problem by introducing viscosity. This guarantees the convergence of the scheme for solutions of the regularized problems in both cases. We also make precise the conditions on the relevant coefficients for convergence to hold.
rate research
Read More
In this work, an efficient numerical scheme is presented for seismic blind deconvolution in a multichannel scenario. The proposed method iterate with wo steps: first, wavelet estimation across all channels and second, refinement of the reflectivity estimate simultaneously in all channels using sparse deconvolution. The reflectivity update step is formulated as a basis pursuit denoising problem and a sparse solution is obtained with the spectral projected-gradient algorithm - faithfulness to the recorded traces is constrained by the measured noise level. Wavelet re-estimation has a closed form solution when performed in the frequency domain by finding the minimum energy wavelet common to all channels. Nothing is assumed known about the wavelet apart from its time duration. In tests with both synthetic and real data, the method yields sparse reflectivity series and stable wavelet estimates results compared to existing methods with significantly less computational effort.
The Reynolds-Averaged Navier-Stokes equations and the Large-Eddy Simulation equations can be coupled using a transition function to switch from a set of equations applied in some areas of a domain to the other set in the other part of the domain. Following this idea, different time integration schemes can be coupled. In this context, we developed a hybrid time integration scheme that spatially couples the explicit scheme of Heun and the implicit scheme of Crank and Nicolson using a dedicated transition function. This scheme is linearly stable and second-order accurate. In this paper, an extension of this hybrid scheme is introduced to deal with a temporal adaptive procedure. The idea is to treat the time integration procedure with unstructured grids as it is performed with Cartesian grids with local mesh refinement. Depending on its characteristic size, each mesh cell is assigned a rank. And for two cells from two consecutive ranks, the ratio of the associated time steps for time marching the solutions is $2$. As a consequence, the cells with the lowest rank iterate more than the other ones to reach the same physical time. In a finite-volume context, a key ingredient is to keep the conservation property for the interfaces that separate two cells of different ranks. After introducing the different schemes, the paper recalls briefly the coupling procedure, and details the extension to the temporal adaptive procedure. The new time integration scheme is validated with the propagation of 1D wave packet, the Sods tube, and the transport of a bi-dimensional vortex in an uniform flow.
Seismology is the main tool for inferring the deep interior structures of Earth and potentially also of other planetary bodies in the solar system. Terrestrial seismology is influenced by the presence of the ocean-generated microseismic signal, which sets a lower limit on the earthquake detection capabilities but also provides a strong energy source to infer the interior structure on scales from local to continental. Titan is the only other place in the solar system with permanent surface liquids and future lander missions there might carry a seismic package. Therefore, the presence of microseisms would be of great benefit for interior studies, but also for detecting storm-generated waves on the lakes remotely. We estimated the strength of microseismic signals on Titan, based on wind speeds predicted from modeled global circulation models interior structure. We find that storms of more than 2 m/s wind speed, would create a signal that is globally observable with a high-quality broadband sensor and observable to a thousand kilometer distance with a space-ready seismometer, such as the InSight instruments currently operating on the surface of Mars.
We develop an approach for simulating acousto-elastic wave phenomena, including scattering from fluid-solid boundaries, where the solid is allowed to be anisotropic, with the Discontinuous Galerkin method. We use a coupled first-order elastic strain-velocity, acoustic velocity-pressure formulation, and append penalty terms based on interior boundary continuity conditions to the numerical (central) flux so that the consistency condition holds for the discretized Discontinuous Galerkin weak formulation. We incorporate the fluid-solid boundaries through these penalty terms and obtain a stable algorithm. Our approach avoids the diagonalization into polarized wave constituents such as in the approach based on solving elementwise Riemann problems.
Fireballs are infrequently recorded by seismic sensors on the ground. If recorded, they are usually reported as one-off events. This study is the first seismic bulk analysis of the largest single fireball data set, observed by the Desert Fireball Network (DFN) in Australia in the period 2014-2019. The DFN typically observes fireballs from cm-m scale impactors. We identified 25 fireballs in seismic time series data recorded by the Australian National Seismograph Network (ANSN). This corresponds to 1.8% of surveyed fireballs, at the kinetic energy range of 10$^6$ to 10$^{10}$ J. The peaks observed in the seismic time series data were consistent with calculated arrival times of the direct airwave or ground-coupled Rayleigh wave caused by shock waves by the fireball in the atmosphere (either due to fragmentation or the passage of the Mach cone). Our work suggests that identification of fireball events in the seismic time series data depends both on physical properties of a fireball (such as fireball energy and entry angle in the atmosphere) and the sensitivity of a seismic instrument. This work suggests that fireballs are likely detectable within 200 km direct air distance between a fireball and seismic station, for sensors used in the ANSN. If each DFN observatory had been accompanied by a seismic sensor of similar sensitivity, 50% of surveyed fireballs could have been detected. These statistics justify the future consideration of expanding the DFN camera network into the seismic domain.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
