Visualization of elastic wavefields computed with a finite difference code
Larsen, S.; Harris, D.
1994-11-15
The authors have developed a finite difference elastic propagation model to simulate seismic wave propagation through geophysically complex regions. To facilitate debugging and to assist seismologists in interpreting the seismograms generated by the code, they have developed an X Windows interface that permits viewing of successive temporal snapshots of the (2D) wavefield as they are calculated. The authors present a brief video displaying the generation of seismic waves by an explosive source on a continent, which propagate to the edge of the continent then convert to two types of acoustic waves. This sample calculation was part of an effort to study the potential of offshore hydroacoustic systems to monitor seismic events occurring onshore.
3D Elastic Wavefield Tomography
NASA Astrophysics Data System (ADS)
Guasch, L.; Warner, M.; Stekl, I.; Umpleby, A.; Shah, N.
2010-12-01
Wavefield tomography, or waveform inversion, aims to extract the maximum information from seismic data by matching trace by trace the response of the solid earth to seismic waves using numerical modelling tools. Its first formulation dates from the early 80's, when Albert Tarantola developed a solid theoretical basis that is still used today with little change. Due to computational limitations, the application of the method to 3D problems has been unaffordable until a few years ago, and then only under the acoustic approximation. Although acoustic wavefield tomography is widely used, a complete solution of the seismic inversion problem requires that we account properly for the physics of wave propagation, and so must include elastic effects. We have developed a 3D tomographic wavefield inversion code that incorporates the full elastic wave equation. The bottle neck of the different implementations is the forward modelling algorithm that generates the synthetic data to be compared with the field seismograms as well as the backpropagation of the residuals needed to form the direction update of the model parameters. Furthermore, one or two extra modelling runs are needed in order to calculate the step-length. Our approach uses a FD scheme explicit time-stepping by finite differences that are 4th order in space and 2nd order in time, which is a 3D version of the one developed by Jean Virieux in 1986. We chose the time domain because an explicit time scheme is much less demanding in terms of memory than its frequency domain analogue, although the discussion of wich domain is more efficient still remains open. We calculate the parameter gradients for Vp and Vs by correlating the normal and shear stress wavefields respectively. A straightforward application would lead to the storage of the wavefield at all grid points at each time-step. We tackled this problem using two different approaches. The first one makes better use of resources for small models of dimension equal
Modeling and Processing of Continuous 3D Elastic Wavefield Data
NASA Astrophysics Data System (ADS)
Milkereit, B.; Bohlen, T.
2001-12-01
Continuous seismic wavefields are excited by earthquake clustering, induced seismicity in reservoirs, and mining. In hydrocarbon reservoirs, for example, pore pressure changes and fluid flow (mass transfer) will cause incremental deviatoric stresses sufficient to trigger and sustain seismic activity. Here we address three aspects of seismic wavefields in three-dimensional heterogeneous media triggered by distributed sources in space and time: forward modeling, multichannel data processing, and source location imaging. A power law distribution of seismic sources (such as the Gutenberg-Richter law) is used for the modeling of viscoelastic/elastic wave propagation through a realistic earth model. 3D modeling provides new insight in the interaction of multi-source wavefields and the role of scale-dependend elastic model parameters on transmitted and reflected/back-scattered wavefields. There exists a strong correlation between the spatial properties of the compressional, shear wave and density perturbations and the lateral correlation length of the resulting reflected or transmitted seismic wavefields. Modeling is based on the implementation of 3D elastic/viscoelastic FD codes on massive parallel and/or distributed computing resources using MPI (message passing interface). For parallelization, large grid 3D earth models are decomposed into subvolume processing elements whereby each processing element is updating the wavefield within its portion of the grid. Processing of continuous seismic wavefields excited by multiple distributed sources is based on a combination of crosscorrelated or slowness-transformed array data and Kirchhoff or reverse time migration for source location or source volume imaging. The appearance of slowness in both migration and array data processing suggests the possibility of combining them into a single process. In order to place further constraints on the migration, the directivity properties of 3-component receiver arrays can be included in
ELASTIC-WAVEFIELD SEISMIC STRATIGRAPHY: A NEW SEISMIC IMAGING TECHNOLOGY
Bob A. Hardage
2004-05-06
The focus of elastic-wavefield seismic stratigraphy research shifted from onshore prospects to marine environments during this report period. Four-component ocean-bottom-cable (4-C OBC) seismic data acquired in water depths of 2400 to 2500 feet across Green Canyon Block 237 in the Gulf of Mexico were processed and analyzed. The P-P and P-SV images of strata immediately below the seafloor exhibit amazing differences in P-P and P-SV seismic facies. These data may be one of the classic examples of the basic concepts of elastic-wavefield seismic stratigraphy.
Elastic-Wavefield Seismic Stratigraphy: A New Seismic Imaging Technology
Bob A. Hardage; Milo M. Backus; Michael V. DeAngelo; Sergey Fomel; Khaled Fouad; Robert J. Graebner; Paul E. Murray; Randy Remington; Diana Sava
2006-07-31
The purpose of our research has been to develop and demonstrate a seismic technology that will provide the oil and gas industry a better methodology for understanding reservoir and seal architectures and for improving interpretations of hydrocarbon systems. Our research goal was to expand the valuable science of seismic stratigraphy beyond the constraints of compressional (P-P) seismic data by using all modes (P-P, P-SV, SH-SH, SV-SV, SV-P) of a seismic elastic wavefield to define depositional sequences and facies. Our objective was to demonstrate that one or more modes of an elastic wavefield may image stratal surfaces across some stratigraphic intervals that are not seen by companion wave modes and thus provide different, but equally valid, information regarding depositional sequences and sedimentary facies within that interval. We use the term elastic wavefield stratigraphy to describe the methodology we use to integrate seismic sequences and seismic facies from all modes of an elastic wavefield into a seismic interpretation. We interpreted both onshore and marine multicomponent seismic surveys to select the data examples that we use to document the principles of elastic wavefield stratigraphy. We have also used examples from published papers that illustrate some concepts better than did the multicomponent seismic data that were available for our analysis. In each interpretation study, we used rock physics modeling to explain how and why certain geological conditions caused differences in P and S reflectivities that resulted in P-wave seismic sequences and facies being different from depth-equivalent S-wave sequences and facies across the targets we studied.
NASA Astrophysics Data System (ADS)
Ha, Jiho; Shin, Sungryul; Shin, Changsoo; Chung, Wookeen
2015-05-01
Because complex mixed waves are typically generated in elastic media, wavefield decomposition is required for such media to obtain migration images accurately. In isotropic media, this is achieved according to the Helmholtz decomposition theorem; in particular, the divergence operator is commonly applied to P-wavefield decomposition. In this study, two types of elastic reverse-time migration algorithms are proposed for decomposition of the P-wavefield without requiring the divergence operator. The first algorithm involves formulation of the stress tensor by spatially differentiated displacement according to the stress-strain relationship and is utilized to construct an imaging condition for the decomposed P-wavefield. We demonstrate this approach through numerical testing. The second algorithm allows us to obtain emphasized interfaces through the application of the absolute value function to decomposed wavefield in imaging condition. Because reverse-time migration can be defined by a zero-lag cross-correlation relationship between the partial-derivative wavefield and the observed wavefield data, we derive the virtual source to construct the partial-derivative wavefield based on a 2D staggered-grid finite-difference modeling method in the time domain. The explicitly computed partial-derivative wavefield from virtual sources with the stress tensor is in agreement with the partial-derivative wavefield directly computed from residual by between with and without a perturbation point in the subsurface. Moreover, the back-propagation technique is used to enhance the computational efficiency. To validate our two types of imaging conditions, numerical tests are conducted. The migration images created according to our imaging conditions can represent the subsurface structure accurately. Thus, we can confirm the feasibility of obtaining migration images of the decomposed P-wavefield without requiring the application of the divergence operator.
Wavefield separation and polarity reversal correction in elastic reverse time migration
NASA Astrophysics Data System (ADS)
Li, Zhiyuan; Ma, Xiaona; Fu, Chao; Liang, Guanghe
2016-04-01
In elastic reverse time migration (RTM), one of the problems that are often encountered is the cross-talk between P- and S-wavefields. A useful processing technique to reduce the cross-talk is separating the P- and S-wavefields by using divergence and curl operators before applying an elastic imaging condition. However, the separated wavefields lose their physical meaning because their phase and amplitude are changed. In this paper, we modify the divergence and curl operators to give the separated wavefields a clear physical meaning: the separated wavefield is the first derivative of the input wavefield with respect to time. Another problem often encountered is polarity reversals in PS and SP images, which can cause destructive interference in the final stacked image and thus destroy the migrated events. In this paper we also develop a procedure for polarity reversal correction based on the polarization vectors of the P- and S-wavefields in the common-shot domain. The correction factor is first calculated at every imaging point during the wavefield reconstruction and is then multiplied by the PS and SP images at each time step when an elastic imaging condition is applied. Numerical examples with synthetic data have shown that the modified wavefield separation method is correct, and the procedure of polarity reversal correction is effective for a complex model.
Pre-stack full wavefield inversion for elastic parameters of TI media
NASA Astrophysics Data System (ADS)
Zhang, Meigen; Huang, Zhongyu; Li, Xiaofan; Wang, Miaoyue; Xu, Guangyin
2006-03-01
Pre-stack full wavefield inversion for the elastic parameters of transversely isotropical media is implemented. The Jacobian matrix is derived directly with the finite element method, just like the full wavefield forward modelling. An absorbing boundary scheme combining Liao's transparent boundary condition with Sarma's attenuation boundary condition is applied to the forward modelling and Jacobian calculation. The input data are the complete ground-recorded wavefields containing full kinematic and dynamic information for the seismic waves. Inversion with such data is desirable as it should improve the accuracy of the estimated parameters and also reduce data pre-processing, such as wavefield identification and separation. A scheme called energy grading inversion is presented to deal with the instability caused by the large energy difference between different arrivals in the input data. With this method, parameters in the shallow areas, which mainly affect wave patterns with strong energy, converge before those of deeper media. Thus, the number of unknowns in each inversion step is reduced, and the stability and reliability of the inversion process is greatly improved. As a result, the scheme is helpful to reduce the non-uniqueness in the inversion. Two synthetic examples show that the inversion system is reliable and accurate even though initial models deviate significantly from the actual models. Also, the system can accurately invert for transversely isotropic model parameters even with the introduction of strong random noise.
Implementation of elastic reverse-time migration using wavefield separation in the frequency domain
NASA Astrophysics Data System (ADS)
Chung, Wookeen; Pyun, Sukjoon; Bae, Ho Seuk; Shin, Changsoo; Marfurt, Kurt J.
2012-06-01
Considerable effort has been devoted to the migration of multicomponent data in elastic media with wavefield separation techniques being the most successful. Most of this work has been carried out in the time domain. In this paper, we formulate a multicomponent migration technique in the frequency domain. Reverse-time migration can be viewed as the zero-lag cross-correlation between virtual source and back-propagated wavefields. Cross-correlating the Helmholtz decomposed wavefields rather than directly correlating the vector displacement fields results in sharper, more interpretable images, contaminated by fewer crosstalk artefacts. The end products are separate P and S wave (and if desired, PS and SP) migration images. We test our migration algorithm on synthetic seismic data generated using the SEG/EAGE salt-dome, Overthrust and Marmousi-2 models. We correctly image the location and shape of the target zone for oil exploration using these data sets. Furthermore, we demonstrate that our new migration technique provides good images even when the initial velocity model is only approximate.
Wavefield imaging using the energy norm
NASA Astrophysics Data System (ADS)
Rocha, Daniel Carvalho, Junior
For various wavefield types, one can formulate a measurement of the mechanical energy that satisfies conservation laws. Based on this formulation, one can derive an energy norm that is applicable to wavefield imaging. Extending the concept of the norm to an inner product, one can compare two related wavefields. Therefore, an imaging condition can be defined as the inner product between the source and receiver wavefields at every spatial location. In this regard, the imaging condition based on the energy inner product accounts for wavefield directionality in space and time, overcoming some problems present in the conventional imaging condition. I exploit the wavefield directionality information from the energy imaging condition to attenuate unwanted events in reverse time migrated (RTM) images. For acoustic wavefields, these unwanted events are characterized by the collinearity of the source and receiver raypaths, and they are described as RTM backscattering artifacts. For elastic wavefields, these events are characterized by the fact that source and receiver displacement fields have the same polarization and wave propagation directions. In both acoustic and elastic cases, one can to attenuate these artifacts and produce high quality images. Another application that uses the wavefield directionality is to enhance the full waveform inversion (FWI) gradient for acoustic wavefields. By enhancing wave events that are collinear and suppressing all other wave events, I am able to compute gradients that are more suitable for the inversion process. Numerical experiments show the efficacy of these applications for synthetic models that emulate the complexity of subsurface structures found in exploration seismology, such as salt bodies, diffractors, dipping layers and faults.
NASA Astrophysics Data System (ADS)
MacLean, L. S.; Romanowicz, B. A.; French, S.
2015-12-01
Seismic wavefield computations using the Spectral Element Method are now regularly used to recover tomographic images of the upper mantle and crust at the local, regional, and global scales (e.g. Fichtner et al., GJI, 2009; Tape et al., Science 2010; Lekic and Romanowicz, GJI, 2011; French and Romanowicz, GJI, 2014). However, the heaviness of the computations remains a challenge, and contributes to limiting the resolution of the produced images. Using source stacking, as suggested by Capdeville et al. (GJI,2005), can considerably speed up the process by reducing the wavefield computations to only one per each set of N sources. This method was demonstrated through synthetic tests on low frequency datasets, and therefore should work for global mantle tomography. However, the large amplitudes of surface waves dominates the stacked seismograms and these cases can no longer be separated by windowing in the time domain. We have developed a processing approach that helps address this issue and demonstrate its usefulness through a series of synthetic tests performed at long periods (T >60 s) on toy upper mantle models. The summed synthetics are computed using the CSEM code (Capdeville et al., 2002). As for the inverse part of the procedure, we use a quasi-Newton method, computing Frechet derivatives and Hessian using normal mode perturbation theory.
NASA Astrophysics Data System (ADS)
Oh, Ju-Won; Alkhalifah, Tariq
2016-07-01
Multi-parameter full waveform inversion (FWI) applied to an elastic orthorhombic model description of the subsurface requires in theory a nine-parameter representation of each pixel of the model. Even with optimal acquisition on the Earth surface that includes large offsets, full azimuth, and multi component sensors, the potential for tradeoff between the elastic orthorhombic parameters are large. The first step to understanding such trade-off is analysing the scattering potential of each parameter, and specifically, its scattering radiation patterns. We investigate such radiation patterns for diffraction and for scattering from a horizontal reflector considering a background isotropic model. The radiation patterns show considerable potential for trade-off between the parameters and the potentially limited resolution in their recovery. The radiation patterns of C11, C22 and C33 are well separated so that we expect to recover these parameters with limited trade-offs. However, the resolution of their recovery represented by recovered range of model wavenumbers varies between these parameters. We can only invert for the short wavelength components (reflection) of C33 while we can mainly invert for the long wavelength components (transmission) of the elastic coefficients C11 and C22 if we have large enough offsets. The elastic coefficients C13, C23 and C12 suffer from strong trade-offs with C55, C44 and C66, respectively. The trade-offs between C13 and C55, as well as C23 and C44, can be partially mitigated if we acquire P-SV and SV-SV waves. However, to reduce the trade-offs between C12 and C66, we require credible SH-SH waves. The analytical radiation patterns of the elastic constants are supported by numerical gradients of these parameters.
NASA Astrophysics Data System (ADS)
Oh, Ju-Won; Alkhalifah, Tariq
2016-09-01
Multiparameter full waveform inversion (FWI) applied to an elastic orthorhombic model description of the subsurface requires in theory a nine-parameter representation of each pixel of the model. Even with optimal acquisition on the Earth surface that includes large offsets, full azimuth, and multicomponent sensors, the potential for trade-off between the elastic orthorhombic parameters are large. The first step to understanding such trade-off is analysing the scattering potential of each parameter, and specifically, its scattering radiation patterns. We investigate such radiation patterns for diffraction and for scattering from a horizontal reflector considering a background isotropic model. The radiation patterns show considerable potential for trade-off between the parameters and the potentially limited resolution in their recovery. The radiation patterns of C11, C22, and C33 are well separated so that we expect to recover these parameters with limited trade-offs. However, the resolution of their recovery represented by recovered range of model wavenumbers varies between these parameters. We can only invert for the short wavelength components (reflection) of C33 while we can mainly invert for the long wavelength components (transmission) of the elastic coefficients C11 and C22 if we have large enough offsets. The elastic coefficients C13, C23, and C12 suffer from strong trade-offs with C55, C44, and C66, respectively. The trade-offs between C13 and C55, as well as C23 and C44, can be partially mitigated if we acquire P-SV and SV-SV waves. However, to reduce the trade-offs between C12 and C66, we require credible SH-SH waves. The analytical radiation patterns of the elastic constants are supported by numerical gradients of these parameters.
Source-independent full wavefield converted-phase elastic migration velocity analysis
NASA Astrophysics Data System (ADS)
Shabelansky, A. H.; Malcolm, A. E.; Fehler, M. C.; Shang, X.; Rodi, W. L.
2015-02-01
Converted phase (CP) elastic seismic signals are comparable in amplitude to the primary signals recorded at large offsets and have the potential to be used in seismic imaging and velocity analysis. We present an approach for CP elastic wave equation velocity analysis that does not use source information and is applicable to surface-seismic, microseismic, teleseismic and vertical seismic profile (VSP) studies. Our approach is based on the cross-correlation between reflected or transmitted PP and CP PS (and/or SS and CP PS) waves propagated backward in time, and is formulated as an optimization problem with a differential semblance criterion objective function for the simultaneous update of both P- and S-wave velocity models. The merit of this approach is that it is fully data-driven, uses full waveform information, and requires only one elastic backward propagation to form an image rather than the two (one forward and one backward) propagations needed for standard reverse-time migration. Moreover, as the method does not require forward propagation, it does not suffer from migration operator source aliasing when a small number of shots are used. We present a derivation of the method and test it with a synthetic model and field micro-seismic data.
EXPLORING FOR SUBTLE MISSION CANYON STRATIGRAPHIC TRAPS WITH ELASTIC WAVEFIELD SEISMIC TECHNOLOGY
John Beecherl
2004-02-01
The 9C3D seismic data that will form the principal data base needed for this research program have been successfully acquired. The seismic field data exhibit a good signal-to-noise (S/N) ratio for all elastic-wave modes. Thus the major hurdle of acquiring optimal-quality 9-C seismic data has been cleared. The stratigraphic oil-reservoir target that will be the imaging objective of the seismic data-processing effort is described in this report to indicate the challenge that now confronts the data-processing phase of the project.
Pollitz, F.F.
2002-01-01
I present a new algorithm for calculating seismic wave propagation through a three-dimensional heterogeneous medium using the framework of mode coupling theory originally developed to perform very low frequency (f < ???0.01-0.05 Hz) seismic wavefield computation. It is a Greens function approach for multiple scattering within a defined volume and employs a truncated traveling wave basis set using the locked mode approximation. Interactions between incident and scattered wavefields are prescribed by mode coupling theory and account for the coupling among surface waves, body waves, and evanescent waves. The described algorithm is, in principle, applicable to global and regional wave propagation problems, but I focus on higher frequency (typically f ??????0.25 Hz) applications at regional and local distances where the locked mode approximation is best utilized and which involve wavefields strongly shaped by propagation through a highly heterogeneous crust. Synthetic examples are shown for P-SV-wave propagation through a semi-ellipsoidal basin and SH-wave propagation through a fault zone.
NASA Astrophysics Data System (ADS)
Pierre, C.
2015-12-01
The Earthscope TA deployment across the continental United-State (US) has reached its eastern part, providing the opportunity for high-resolution 3D seismic velocity imaging of both lithosphere and asthenosphere across the entire north-American continent (NA). Previously (Yuan et al., 2014), we presented a 3D radially anisotropic shear wave (Vs) model of North America (NA) lithospheric mantle based on full waveform tomography, combining teleseismic and regional distance data sampling the NA. Regional wavefield computations were performed numerically, using a regional Spectral Element code (RegSEM, Cupillard et al., 2012), while teleseismic computations were performed approximately, using non-linear asymptotic coupling theory (NACT, Li and Romanowicz, 1995). For both datasets, the inversion was performed iteratively, using a Gauss-Newton scheme, with kernels computed using either NACT or the surface wave, path average approximation (PAVA), depending on the source-station distance. We here present a new radially anisotropic lithospheric/asthenospheric model of Vs for NA based entirely on SEM-based numerical waveforms from an augmented dataset of 155 regional events and 70 teleseismic events. The forward wavefield computations are performed using RegSEM down to 40s, starting from our most recent whole mantle 3D radially anisotropic Vs model (SEMUCB-wm1, French and Romanowicz, 2014). To model teleseismic wavefields within our regional computational domain, we developed a new modeling technique which allows us to replace a distant source by virtual sources at the boundary of the computational domain (Masson et al., 2014). Computing virtual sources requires one global simulation per teleseismic events.We then compare two models obtained: one using NACT/PAVA kernels as in our previous work, and another using hybrid kernels, where the Hessian is computed using NACT/PAVA, but the gradient is computed numerically from the adjoint wavefield, providing more accurate kernels
Finite-difference modelling of wavefield constituents
NASA Astrophysics Data System (ADS)
Robertsson, Johan O. A.; van Manen, Dirk-Jan; Schmelzbach, Cedric; Van Renterghem, Cederic; Amundsen, Lasse
2015-11-01
The finite-difference method is among the most popular methods for modelling seismic wave propagation. Although the method has enjoyed huge success for its ability to produce full wavefield seismograms in complex models, it has one major limitation which is of critical importance for many modelling applications; to naturally output up- and downgoing and P- and S-wave constituents of synthesized seismograms. In this paper, we show how such wavefield constituents can be isolated in finite-difference-computed synthetics in complex models with high numerical precision by means of a simple algorithm. The description focuses on up- and downgoing and P- and S-wave separation of data generated using an isotropic elastic finite-difference modelling method. However, the same principles can also be applied to acoustic, electromagnetic and other wave equations.
NASA Astrophysics Data System (ADS)
Pierre, C.; Masson, Y.; Romanowicz, B. A.; French, S. W.; Yuan, H.
2014-12-01
The Earthscope TA deployment across the continental US now has reached the eastern part of the United States, providing the opportunity for high-resolution 3D seismic velocity imaging of both lithosphere and asthenosphere across the entire north-American continent (NA). Previously (Yuan et al., 2014), we presented a 3D radially anisotropic shear wave model of North America (NA) lithospheric mantle based on full waveform tomography, combining teleseismic and regional distance data sampling the NA. Regional wavefield computations were performed numerically, using a regional Spectral Element code (RegSEM, Cupillard et al., 2012), while teleseismic computations were performed approximately, using non-linear asymptotic coupling theory (NACT, Li and Romanowicz, 1995). For both datasets, the inversion was performed iteratively, using a Gauss-Newton scheme, with kernels computed using either NACT or the surface wave, path average approximation (PAVA), depending on the source-station distance. Building upon our previous work, we here present a new radially anisotropic lithospheric/asthenospheric model of shear velocity for North America based entirely on regional waveforms from an augmented dataset of ~150 events contained and observed inside the study region, with forward wavefield computations performed using RegSEM down to 40s, starting from our most recent whole mantle 3D radially anisotropic shear velocity model (SEMUCB-wm1, French and Romanowicz, 2014). Several iterations of inversion are performed using a Gauss-Newton scheme. We present and compare two models obtained, on the one hand, using NACT/PAVA kernels as in our previous work, and on the other, using hybrid kernels, where the Hessian is computed using NACT/PAVA, but the gradient is computed numerically from the adjoint wavefield, providing more accurate kernels while preserving the fast convergence properties of the Gauss-Newton inversion scheme. We also present an update to our azimuthally anisotropic shear
NASA Technical Reports Server (NTRS)
Schmidt, H.; Tango, G. J.; Werby, M. F.
1985-01-01
A new matrix method for rapid wave propagation modeling in generalized stratified media, which has recently been applied to numerical simulations in diverse areas of underwater acoustics, solid earth seismology, and nondestructive ultrasonic scattering is explained and illustrated. A portion of recent efforts jointly undertaken at NATOSACLANT and NORDA Numerical Modeling groups in developing, implementing, and testing a new fast general-applications wave propagation algorithm, SAFARI, formulated at SACLANT is summarized. The present general-applications SAFARI program uses a Direct Global Matrix Approach to multilayer Green's function calculation. A rapid and unconditionally stable solution is readily obtained via simple Gaussian ellimination on the resulting sparsely banded block system, precisely analogous to that arising in the Finite Element Method. The resulting gains in accuracy and computational speed allow consideration of much larger multilayered air/ocean/Earth/engineering material media models, for many more source-receiver configurations than previously possible. The validity and versatility of the SAFARI-DGM method is demonstrated by reviewing three practical examples of engineering interest, drawn from ocean acoustics, engineering seismology and ultrasonic scattering.
NASA Astrophysics Data System (ADS)
Schuberth, Bernhard; Zaroli, Christophe; Nolet, Guust
2015-04-01
Of particular interest for the tectonic evolution of the Atlantic region is the influence of lower mantle structure under Africa on flow in the upper mantle beneath the ocean basin. Along with its Pacific counterpart, the large African anomaly in the lowermost mantle with strongly reduced seismic velocities has received considerable attention in seismological and geodynamic studies. Several seismological observations are typically taken as an indication that these two anomalies are being caused by large-scale compositional variations and that they are piles of material with higher density than normal mantle rock. This would imply negative buoyancy in the lowermost mantle under Africa, which has important implications for the flow at shallower depth and inferences on the processes that led to the formation of the Atlantic Ocean basin. However, a large number of recent studies argue for a strong thermal gradient across the core-mantle boundary that might provide an alternative explanation for the lower mantle anomaly through the resulting large lateral temperature variations. Recently, we developed a new joint forward modeling approach to test such geodynamic hypotheses directly against the seismic observations: Seismic heterogeneity is predicted by converting the temperature field of a high-resolution 3-D mantle circulation model into seismic velocities using thermodynamic models of mantle mineralogy. 3-D global wave propagation in the synthetic elastic structures is then simulated using a spectral element method. Being based on forward modelling only, this approach allows us to generate synthetic wavefields and seismograms independently of seismic observations. The statistics of observed long-period body wave traveltime variations show a markedly different behaviour for P- and S-waves: the standard deviation of P-wave delay times stays almost constant with ray turning depth, while that of the S-wave delay times increases strongly throughout the mantle. In an
Wavefield Compression for Full-Waveform Inversion
NASA Astrophysics Data System (ADS)
Boehm, Christian; Fichtner, Andreas; de la Puente, Josep; Hanzich, Mauricio
2015-04-01
We present compression techniques tailored to iterative nonlinear minimization methods that significantly reduce the memory requirements to store the forward wavefield for the computation of sensitivity kernels. Full-waveform inversion on 3d data sets requires massive computing and memory capabilities. Adjoint techniques offer a powerful tool to compute the first and second derivatives. However, due to the asynchronous nature of forward and adjoint simulations, a severe bottleneck is introduced by the necessity to access both wavefields simultaneously when computing sensitivity kernels. There exist two opposing strategies to deal with this challenge. On the one hand, conventional approaches save the whole forward wavefield to the disk, which yields a significant I/O overhead and might require several terabytes of storage capacity per seismic event. On the other hand, checkpointing techniques allow to trade an almost arbitrary amount of memory requirements for a - potentially large - number of additional forward simulations. We propose an alternative approach that strikes a balance between memory requirements and the need for additional computations. Here, we aim at compressing the forward wavefield in such a way that (1) the I/O overhead is reduced substantially without the need for additional simulations, (2) the costs for compressing/decompressing the wavefield are negligible, and (3) the approximate derivatives resulting from the compressed forward wavefield do not affect the rate of convergence of a Newton-type minimization method. To this end, we apply an adaptive re-quantization of the displacement field that uses dynamically adjusted floating-point accuracies - i.e., a locally varying number of bits - to store the data. Furthermore, the spectral element functions are adaptively downsampled to a lower polynomial degree. In addition, a sliding-window cubic spline re-interpolates the temporal snapshots to recover a smooth signal. Moreover, a preprocessing step
NASA Astrophysics Data System (ADS)
Komatitsch, Dimitri; Monteiller, Vadim; Chevrot, Sébastien; Wang, Yi; Durochat, Clément
2015-04-01
We present a method for high-resolution imaging of lithospheric structures based on full waveform inversion of teleseismic wavefields. We model the propagation of seismic waves using our recently developed Direct Solution Method (DSM) / Spectral-Element Method (SEM) hybrid technique, which allows us to simulate the propagation of short period teleseismic waves through a regional 3-D model. We implement an iterative quasi-Newton method based upon the L-BFGS algorithm, with a gradient of the misfit function computed with the adjoint-state method. Compared to gradient or conjugate-gradient methods, the L-BFGS algorithm finds solutions that better explain the observed waveforms, and at a much faster convergence rate. We illustrate the potential of this method on a synthetic test case that consists in a crustal model with a crustal discontinuity at 25 km depth and a sharp Moho jump. This simple model contains short and long wavelength heterogeneities along both the lateral and vertical dimensions. In order to do that successfully we resort to high-performance computing on supercomputing clusters using an improved version of our SPECFEM3D open-source software package, which exhibits excellent scalability on parallel machines.
True-amplitude wavefield separation using staggered-grid interpolation in the wavenumber domain
NASA Astrophysics Data System (ADS)
Du, Qi-Zhen; Zhang, Ming-Qiang; Chen, Xiao-Ran; Gong, Xu-Fei; Guo, Cheng-Feng
2014-12-01
Wavefield separation of multicomponent seismic data to image subsurface structures can be realized in either the space domain or the wavenumber domain. However, as the particle velocity components used in the wavenumber-domain wavefield separation are not defined at the same grid point with the staggered-grid finite-difference method for elastic wavefield simulation, we propose the wavenumber-domain interpolation method to estimate the required values at the common grid points prior to the wavenumber-domain true-amplitude wavefield separation. Moreover, numerical experiments show that the wavenumber-domain interpolation method has high interpolation accuracy and the trueamplitude wavefield separation method shows good amplitude preservation. The application of the proposed methodology to elastic reverse-time migration can obtain good amplitudepreserved images even in the case of some velocity error.
Guided wavefield reconstruction from sparse measurements
NASA Astrophysics Data System (ADS)
Mesnil, Olivier; Ruzzene, Massimo
2016-02-01
Guided wave measurements are at the basis of several Non-Destructive Evaluation (NDE) techniques. Although sparse measurements of guided wave obtained using piezoelectric sensors can efficiently detect and locate defects, extensive informa-tion on the shape and subsurface location of defects can be extracted from full-field measurements acquired by Laser Doppler Vibrometers (LDV). Wavefield acquisition from LDVs is generally a slow operation due to the fact that the wave propagation to record must be repeated for each point measurement and the initial conditions must be reached between each measurement. In this research, a Sparse Wavefield Reconstruction (SWR) process using Compressed Sensing is developed. The goal of this technique is to reduce the number of point measurements needed to apply NDE techniques by at least one order of magnitude by extrapolating the knowledge of a few randomly chosen measured pixels over an over-sampled grid. To achieve this, the Lamb wave propagation equation is used to formulate a basis of shape functions in which the wavefield has a sparse representation, in order to comply with the Compressed Sensing requirements and use l1-minimization solvers. The main assumption of this reconstruction process is that every material point of the studied area is a potential source. The Compressed Sensing matrix is defined as being the contribution that would have been received at a measurement location from each possible source, using the dispersion relations of the specimen computed using a Semi-Analytical Finite Element technique. The measurements are then processed through an l1-minimizer to find a minimum corresponding to the set of active sources and their corresponding excitation functions. This minimum represents the best combination of the parameters of the model matching the sparse measurements. Wavefields are then reconstructed using the propagation equation. The set of active sources found by minimization contains all the wave
Prestack elastic generalized-screen migration for multicomponent data
NASA Astrophysics Data System (ADS)
Kim, Byoung Yeop; Seol, Soon Jee; Lee, Ho-Young; Byun, Joongmoo
2016-03-01
An efficient prestack depth migration method based on the elastic one-way wave equation was developed using an improved elastic generalized-screen propagator, which effectively describes the behavior of elastic waves with mode conversion at the interfaces and efficiently computes wave propagation in media with lateral velocity variations. The elastic propagator presented in this study was improved from the elastic generalized-screen propagator. Several terms of the vertical slowness right symbol with orders are corrected from the original formulation, and the vertical slowness operator in the propagator was expanded up to the 2nd order which yields a more accurate approximation. In each screen propagation step, the multicomponent wavefields are automatically separated into the P and S wavefields by the P-S decomposition operator included in the elastic generalized-screen propagator. This process facilitates the imaging of the P and S waves separately without any additional P and S separation process after the wavefield extrapolation. Impulse response tests of the improved elastic generalized-screen propagator in a uniform-property medium proved that propagation accuracy increases with order, even when large medium perturbations are assigned. Migration tests using the developed algorithm on a simple layered model and two complex models (the SEG/EAGE salt and elastic Marmousi-2 model) demonstrated the functional advantages and capabilities of the algorithm compared with elastic migration using the scalar wave equation.
Explicit Fourier wavefield operators
NASA Astrophysics Data System (ADS)
Ferguson, R. J.; Margrave, G. F.
2006-04-01
Explicit wavefield extrapolators are based on direct analytic mathematical formulae that express the output as an extrapolation operator acting on the input, while implicit methods usually require the calculation of the numerical inverse of a matrix to obtain the output. Typically, explicit methods are faster than implicit methods, and they often give more insight into the physics of the wave propagation, but they often suffer from instability. Four different explicit extrapolators based on Fourier theory are presented and analysed. They are: PS (ordinary phase shift), GPSPI (generalized phase shift plus interpolation), NSPS (non-stationary phase shift) and SNPS (symmetric non-stationary phase shift). A formal proof is given that NSPS in a direction orthogonal to the velocity gradient is the mathematical adjoint process to GPSPI in the opposite direction. This motivates the construction of SNPS that combines NSPS and GPSPI in a symmetric fashion. This symmetry (under interchange of input and output lateral coordinates) is required by reciprocity arguments. PS and SNPS are symmetric while NSPS and GPSPI are not. A numerical stability study using SVD (singular value decomposition) shows that all of these extrapolators can become unstable for strong lateral velocity gradients. Unstable operators allow amplitudes to grow non-physically in a recursion. Stability is enhanced by introducing a small (~3 per cent) imaginary component to the velocities. This causes a numerical attenuation that tends to stabilize the operators but does not address the cause of the instability. For the velocity model studied (a very challenging case) GPSPI and NSPS have exactly the same instability while SNPS is always more stable. Instability manifests in a complicated way as a function of extrapolation step size, frequency, velocity gradient, and strength of numerical attenuation. The SNPS operator can be stabilized over a wide range of conditions with considerably less attenuation than is
An excitation potential imaging condition for elastic reverse time migration
NASA Astrophysics Data System (ADS)
Gu, Bingluo; Liu, Youshan; Li, Zhiyuan; Ma, Xiaona; Liang, Guanghe
2014-09-01
Elastic reverse time migration (ERTM) has been demonstrated to be more accurate than scalar RTM. However, low efficiency (large storage and heavy calculated amount) and strong artifacts caused by the crosstalk between different wave modes are the two primary barriers to the application of the ERTM during the processing of real data. The scalar (P) and vector (S) potentials of the elastic wavefield and the arrival times corresponding to the first energy extremum of the wavefield are saved at each grid point during the forward modeling of the source wavefield. The angle-dependent reflection coefficient images are subsequently obtained by dividing the scalar and vector potentials of the backward extrapolated receiver wavefield by the saved scalar and vector potentials at the grid points that satisfy the image time at each time step, respectively. The proposed imaging condition does not need to store the snapshots of the source wavefield, while it can considerably improve the computational efficiency and decrease the amount of storage and Input/Output manipulation (compared with the cross-correlation imaging condition) in addition to suppressing the crosstalk between compressive and shear wave modes. Compared with the excitation time imaging condition, the proposed imaging condition reduces the energy loss caused by the opposite polarity of the horizontal component at opposite sides of the source in stacked images. Numerical tests with synthetic data of the Sigsbee2a model have demonstrated that this imaging condition is a cost-effective and practical imaging condition for use in prestack ERTM.
An Expert Fitness Diagnosis System Based on Elastic Cloud Computing
Tseng, Kevin C.; Wu, Chia-Chuan
2014-01-01
This paper presents an expert diagnosis system based on cloud computing. It classifies a user's fitness level based on supervised machine learning techniques. This system is able to learn and make customized diagnoses according to the user's physiological data, such as age, gender, and body mass index (BMI). In addition, an elastic algorithm based on Poisson distribution is presented to allocate computation resources dynamically. It predicts the required resources in the future according to the exponential moving average of past observations. The experimental results show that Naïve Bayes is the best classifier with the highest accuracy (90.8%) and that the elastic algorithm is able to capture tightly the trend of requests generated from the Internet and thus assign corresponding computation resources to ensure the quality of service. PMID:24723842
Angle-domain imaging condition for elastic reverse time migration
NASA Astrophysics Data System (ADS)
Yan, R.; Xie, X.; Wu, R.
2010-12-01
In exploration seismology, elastic reverse time migration (RTM) has spurred much interest in recent years because of the increased imaging challenges posed by complex subsurface targets and affordable computer resources such as Linux clusters. Elastic RTM reconstructs the source wavefields forward in time and receiver wavefields backward in time by finite difference method. It then applies an imaging condition that evaluates the match between source and receiver wavefields. To construct the image which describes the physical property of the medium, we prefer to separate the wavefields into P and S modes, and implement the imaging condition as cross-correlation of pure wave mode rather than of Cartesian component of the displacement wavefields. However, simple cross-correlation can cause some image problems which impede the further seismic interpretation. For example, PP image is contaminated by strong artifacts resulted from unwanted cross-correlation between diving wave, turning wave or back-scattered wave and their time-reversed counterpart; PS image suffers from polarization problem due to the polarity reversal of converted S-wave. We found it necessary to apply certain intrinsic properties (e.g., the wave propagation directions and particle motion information) in reconstructed wavefields in order to solve the problems exhibited in elastic RTM image. This paper presents a procedure to decompose the source and receiver wavefields into local plane waves in pure P and S modes. We generate the partial PP and PS images in angle domain by cross-correlating any combination of decomposed plane wave component, and then formulate an imaging condition as a product of an angle-domain operator and the partial images. The new angle-domain imaging condition substantially reduces the artifacts in the PP image and produces the PS image with correct polarizations. Synthetic examples demonstrate that the imaging condition works very well on producing clean and consistent image
Computational aspects of dispersive computational continua for elastic heterogeneous media
NASA Astrophysics Data System (ADS)
Fafalis, Dimitrios; Fish, Jacob
2015-12-01
The present manuscript focusses on computational aspects of dispersive computational continua (C^2) formulation previously introduced by the authors. The dispersive C^2 formulation is a multiscale approach that showed strikingly accurate dispersion curves. However, the seemingly theoretical advantage may be inconsequential due to tremendous computational cost involved. Unlike classical dispersive methods pioneered more than a half a century ago where the unit cell is quasi-static and provides effective mechanical and dispersive properties to the coarse-scale problem, the dispersive C^2 gives rise to transient problems at all scales and for all microphases involved. An efficient block time-integration scheme is proposed that takes advantage of the fact that the transient unit cell problems are not coupled to each other, but rather to a single coarse-scale finite element they are positioned in. We show that the computational cost of the method is comparable to the classical dispersive methods for short load durations.
Wavefield tomography using extended images
NASA Astrophysics Data System (ADS)
Yang, Tongning
Estimating an accurate velocity model is crucial for seismic imaging to obtain a good understanding of the subsurface structure. The objective of this thesis is to investigate methods of velocity analysis by optimizing seismic images. A conventional seismic image is obtained by zero-lag crosscorrelation of wavefields extrapolated from a source wavelet and recorded data on the surface using a velocity model. The velocity model provides the kinematic information needed by the imaging algorithm to position the reflectors at correct locations and to focus the image. In complex geology, wave-equation migration is a powerful tool for accurately imaging the earth's interior; the quality of the output image, however, depends on the accuracy of the velocity model. Given such a dependency between the image and model, analyzing the velocity information from the image is still not intuitive and often ambiguous. If the nonzero space- and time-lags information are preserved in the crosscorrelation, the output are image hypercube defined as extended images. Compared to the conventional image, the extended images provide a straightforward way to analyze the image quality and to characterize the velocity model accuracy. Understanding the reflection moveout is the key to developing velocity model building methods using extended images. In the extended image space, reflections form coherent objects which depend on space (lags) and time (lags). These objects resemble cones which ideally have their apex at zero space and time lags. The symmetry axis of the cone lies along the time-lag axis. The apex of the cone is located at zero lags only if the velocity model is accurate. This corresponds to the situation when reflection energy focuses at origin in both the space- and time-lag common-image gathers (the slices at zero time and space lags, respectively). When the velocity model is inaccurate, the cone shifts along the time-lag axis. This results in residual moveout in space-lag gathers
Simplified computational methods for elastic and elastic-plastic fracture problems
NASA Technical Reports Server (NTRS)
Atluri, Satya N.
1992-01-01
An overview is given of some of the recent (1984-1991) developments in computational/analytical methods in the mechanics of fractures. Topics covered include analytical solutions for elliptical or circular cracks embedded in isotropic or transversely isotropic solids, with crack faces being subjected to arbitrary tractions; finite element or boundary element alternating methods for two or three dimensional crack problems; a 'direct stiffness' method for stiffened panels with flexible fasteners and with multiple cracks; multiple site damage near a row of fastener holes; an analysis of cracks with bonded repair patches; methods for the generation of weight functions for two and three dimensional crack problems; and domain-integral methods for elastic-plastic or inelastic crack mechanics.
Simultaneous source wavefield separation by signal apparition
NASA Astrophysics Data System (ADS)
Robertsson, Johan O. A.; Amundsen, Lasse; Pedersen, Åsmund Sjøen
2016-06-01
A new method for discrete sampling of signals is presented with specific applications to the reconstruction of recorded interfering wavefields from two or more sources excited simultaneously at discrete positions along lines. By utilizing a periodic sequence of source signatures along one of the source lines, the corresponding wavefield becomes separately visible in a part of the spectral domain where it can be isolated, processed, and subtracted from the interfering wavefields. As a result, interfering wavefields from multiple sources recorded at a single location can be fully separated from each other. The concept is referred to as signal apparition which we suggest refers to "the act of becoming visible". It may find applications in a wide range of disciplines relying on wave experimentation such as acoustic, seismic, and electromagnetic imaging of the Earth's interior for instance to significantly enhance resolution of sub-surface images.
Passive imaging in nondiffuse acoustic wavefields.
Mulargia, Francesco; Castellaro, Silvia
2008-05-30
A main property of diffuse acoustic wavefields is that, taken any two points, each of them can be seen as the source of waves and the other as the recording station. This property is shown to follow simply from array azimuthal selectivity and Huygens principle in a locally isotropic wavefield. Without time reversal, this property holds approximately also in anisotropic azimuthally uniform wavefields, implying much looser constraints for undistorted passive imaging than those required by a diffuse field. A notable example is the seismic noise field, which is generally nondiffuse, but is found to be compatible with a finite aperture anisotropic uniform wavefield. The theoretical predictions were confirmed by an experiment on seismic noise in the mainland of Venice, Italy. PMID:18518643
Passive Imaging in Nondiffuse Acoustic Wavefields
Mulargia, Francesco; Castellaro, Silvia
2008-05-30
A main property of diffuse acoustic wavefields is that, taken any two points, each of them can be seen as the source of waves and the other as the recording station. This property is shown to follow simply from array azimuthal selectivity and Huygens principle in a locally isotropic wavefield. Without time reversal, this property holds approximately also in anisotropic azimuthally uniform wavefields, implying much looser constraints for undistorted passive imaging than those required by a diffuse field. A notable example is the seismic noise field, which is generally nondiffuse, but is found to be compatible with a finite aperture anisotropic uniform wavefield. The theoretical predictions were confirmed by an experiment on seismic noise in the mainland of Venice, Italy.
Signal apparition for simultaneous source wavefield separation
NASA Astrophysics Data System (ADS)
Robertsson, Johan O. A.; Amundsen, Lasse; Pedersen, Åsmund Sjøen
2016-08-01
A new method for discrete sampling of signals is presented with specific applications to the reconstruction of recorded interfering wavefields from two or more sources excited simultaneously at discrete positions along lines. By utilizing a periodic sequence of source signatures along one of the source lines, the corresponding wavefield becomes separately visible in a part of the spectral domain where it can be isolated, processed and subtracted from the interfering wavefields. As a result, interfering wavefields from multiple sources recorded at a single location can be fully separated from each other. The concept is referred to as signal apparition which we suggest refers to `the act of becoming visible'. It may find applications in a wide range of disciplines relying on wave experimentation, such as acoustic, seismic and electromagnetic imaging of the Earth's interior for instance to significantly enhance resolution of subsurface images.
Computational Model of Three Dimensional Elastic Wing Driven by Muscles
NASA Astrophysics Data System (ADS)
Wang, Z. Jane; Cowen, Nathaniel; Peskin, Charles S.; Childress, Stephen W.
2003-11-01
The flapping wing motion observed in nature results from couplings of muscles, flexible wing structures, and unsteady flows. Previously we have studied the unsteady flows and forces of a rigid two dimensional wing undergoing prescribed motion similar to kinematics observed in insects, as a means of understanding basic unsteady aerodynamic mechanisms. In this talk, we describe our recent progress in constructing a more realistic model insect, which consists of a pair of elastic wings immersed in fluids, and is driven by periodically contracting 'muscles'. A natural computational framework for such a system is the immersed boundary method, which is used here. We present simulations of flapping flight at Reynolds number 10^2, in the same range as that of fruitflies and butterflies.
Seismic wavefield propagation in 2D anisotropic media: Ray theory versus wave-equation simulation
NASA Astrophysics Data System (ADS)
Bai, Chao-ying; Hu, Guang-yi; Zhang, Yan-teng; Li, Zhong-sheng
2014-05-01
Despite the ray theory that is based on the high frequency assumption of the elastic wave-equation, the ray theory and the wave-equation simulation methods should be mutually proof of each other and hence jointly developed, but in fact parallel independent progressively. For this reason, in this paper we try an alternative way to mutually verify and test the computational accuracy and the solution correctness of both the ray theory (the multistage irregular shortest-path method) and the wave-equation simulation method (both the staggered finite difference method and the pseudo-spectral method) in anisotropic VTI and TTI media. Through the analysis and comparison of wavefield snapshot, common source gather profile and synthetic seismogram, it is able not only to verify the accuracy and correctness of each of the methods at least for kinematic features, but also to thoroughly understand the kinematic and dynamic features of the wave propagation in anisotropic media. The results show that both the staggered finite difference method and the pseudo-spectral method are able to yield the same results even for complex anisotropic media (such as a fault model); the multistage irregular shortest-path method is capable of predicting similar kinematic features as the wave-equation simulation method does, which can be used to mutually test each other for methodology accuracy and solution correctness. In addition, with the aid of the ray tracing results, it is easy to identify the multi-phases (or multiples) in the wavefield snapshot, common source point gather seismic section and synthetic seismogram predicted by the wave-equation simulation method, which is a key issue for later seismic application.
NASA Astrophysics Data System (ADS)
Warner, Michael; Shah, Nikhil; Guasch, Lluís.; Umpleby, Adrian; Stekl, Ivan; Morgan, Joanna
2010-05-01
Wavefield tomography (AKA full waveform inversion) is a method of inverting geophysical field data that seeks to find a quantitative model of physical properties in the subsurface that can be used to generate synthetic data that match field data "wiggle-for-wiggle". It is most often applied to active-source seismic data, but it can also be applied to passive seismic data and to controlled-source electro-magnetic data. The method has significantly higher spatial resolution and fidelity than can be achieved using conventional imaging methods. The first practical methods were developed in 2D, mostly in the frequency domain. Here we report recent algorithmic developments which, coupled with hardware advances, make these methods tractable in three dimensions. We have implemented finite-difference computer codes for 3D acoustic and elastic wavefield tomography by explicit time-stepping in the time-domain, and for visco-acoustic tomography by iterative solution of the implicit matrix equations in the frequency-domain. Unlike the situation in two dimensions, where frequency-domain methods have proven to be far more efficient, in three dimensions both methods require approximately similar computational resource, and have largely complementary properties in terms of their effectiveness. We have applied these methods to a variety of synthetic and real-world problems taken from petroleum, mining and academic field datasets with a variety of 3D acquisition geometries and target depths. These include conventional marine multi-streamer acquisition, multi-azimuth marine OBC, high-resolution land surveys, and deep-ocean single-streamer acquisition. In each case, wavefield tomography was able to obtain a high-resolution high-fidelity velocity model of the heterogeneous overburden, and consequently to improve subsequent depth imaging of an underlying target. One of the serious practical limitations on the wider applicability of wavefield tomography is the necessity to have low
The seismic noise wavefield is not diffuse.
Mulargia, Francesco
2012-04-01
Passive seismology is burgeoning under the apparent theoretical support of diffuse acoustics. However, basic physical arguments suggest that this theory may not be applicable to seismic noise. A procedure is developed to establish the applicability of the diffuse field paradigm to a wavefield, based on testing the latter for azimuthal isotropy and spatial homogeneity. This procedure is then applied to the seismic noise recorded at 65 sites covering a wide variety of environmental and subsoil conditions. Considering the instantaneous oscillation vector measured at single triaxial stations, the hypothesis of azimuthal isotropy is rejected in all cases with high confidence, which makes the spatial homogeneity test unnecessary and leads directly to conclude that the seismic noise wavefield is not diffuse. However, such a conclusion has no practical effect on passive imaging, which is also possible in non-diffuse wavefields. PMID:22501063
Using quasiphotons to compute wave fields in an elastic medium
Kachalov, A.P.
1987-07-10
Quasiphoton solutions are constructed for longitudinal and transversal waves in an elastic medium. The quasiphotons are then applied to determine the fields of nonstationary high-frequency point sources in a medium with parameters dependent on two Euclidean coordinates.
Computation of graphene elastic moduli at low temperature
Zubko, I. Yu. Kochurov, V. I.
2015-10-27
Finding the values of parameters for the simplest Mie’s family potentials is performed in order to estimate elastic moduli of graphene monolayers using lattice statics approach. The coincidence criterion of the experimentally determined Poisson’s ratio with the estimated value is taken in order to select dimensionless power parameters of the Mie-type potential. It allowed obtaining more precise estimation of elastic properties in comparison with variety of other potentials for carbon atoms in graphene monolayer.
NASA Astrophysics Data System (ADS)
Minato, Shohei; Ghose, Ranajit
2014-05-01
The inverse scattering of seismic waves can reveal the spatial distribution of the elastic compliances along a non-welded interface, such as a fracture surface. The spatial heterogeneity along the surface of a fracture is a key determinant for fracture-associated hydraulic properties. In this paper, we demonstrate that the inverse scattering solution can be successfully applied to the point source response of a subhorizontal fracture. In the scale of seismic exploration, it is more appropriate to consider spherical waves from a point source than plane waves. Further, from only the P-wave point source response it is possible to estimate both normal and tangential fracture compliances. The synthetic seismic wavefield due to a P-wave point source in a 2-D elastic medium was computed using a time-domain finite difference approach. On this spherical wave data set, the correct estimation of the position and dip of the non-welded interface was possible through reverse-time migration followed by least-square fitting of the maximum amplitude of the P-P reflection. In order to estimate the heterogeneity along the non-welded interface, we first extract the elastic wavefield at the interface position. The extrapolated wavefield is then rotated such that the horizontal axis aligns along the fracture plane. Next, using this extrapolated and rotated wavefield, we solve the linear-slip boundary condition to obtain the distribution of normal and tangential compliances. Our result shows that the estimates of normal compliance are very accurate around the dominant frequency of the incident seismic wavefield. At lower frequencies, the estimated compliance distribution is less accurate and rather smooth due to the presence of evanescent waves. Extracting the distribution of the tangential compliance requires a larger stabilization factor. For a correct estimation of the tangential compliance, one needs S-wave sources or multiple sources providing more grazing angles to avoid the shadow
Using Amazon's Elastic Compute Cloud to dynamically scale CMS computational resources
NASA Astrophysics Data System (ADS)
Evans, D.; Fisk, I.; Holzman, B.; Melo, A.; Metson, S.; Pordes, R.; Sheldon, P.; Tiradani, A.
2011-12-01
Large international scientific collaborations such as the Compact Muon Solenoid (CMS) experiment at the Large Hadron Collider have traditionally addressed their data reduction and analysis needs by building and maintaining dedicated computational infrastructure. Emerging cloud computing services such as Amazon's Elastic Compute Cloud (EC2) offer short-term CPU and storage resources with costs based on usage. These services allow experiments to purchase computing resources as needed, without significant prior planning and without long term investments in facilities and their management. We have demonstrated that services such as EC2 can successfully be integrated into the production-computing model of CMS, and find that they work very well as worker nodes. The cost-structure and transient nature of EC2 services makes them inappropriate for some CMS production services and functions. We also found that the resources are not truely "on-demand" as limits and caps on usage are imposed. Our trial workflows allow us to make a cost comparison between EC2 resources and dedicated CMS resources at a University, and conclude that it is most cost effective to purchase dedicated resources for the "base-line" needs of experiments such as CMS. However, if the ability to use cloud computing resources is built into an experiment's software framework before demand requires their use, cloud computing resources make sense for bursting during times when spikes in usage are required.
Computing elastic moduli on 3-D X-ray computed tomography image stacks
NASA Astrophysics Data System (ADS)
Garboczi, E. J.; Kushch, V. I.
2015-03-01
A numerical task of current interest is to compute the effective elastic properties of a random composite material by operating on a 3D digital image of its microstructure obtained via X-ray computed tomography (CT). The 3-D image is usually sub-sampled since an X-ray CT image is typically of order 10003 voxels or larger, which is considered to be a very large finite element problem. Two main questions for the validity of any such study are then: can the sub-sample size be made sufficiently large to capture enough of the important details of the random microstructure so that the computed moduli can be thought of as accurate, and what boundary conditions should be chosen for these sub-samples? This paper contributes to the answer of both questions by studying a simulated X-ray CT cylindrical microstructure with three phases, cut from a random model system with known elastic properties. A new hybrid numerical method is introduced, which makes use of finite element solutions coupled with exact solutions for elastic moduli of square arrays of parallel cylindrical fibers. The new method allows, in principle, all of the microstructural data to be used when the X-ray CT image is in the form of a cylinder, which is often the case. Appendix A describes a similar algorithm for spherical sub-samples, which may be of use when examining the mechanical properties of particles. Cubic sub-samples are also taken from this simulated X-ray CT structure to investigate the effect of two different kinds of boundary conditions: forced periodic and fixed displacements. It is found that using forced periodic displacements on the non-geometrically periodic cubic sub-samples always gave more accurate results than using fixed displacements, although with about the same precision. The larger the cubic sub-sample, the more accurate and precise was the elastic computation, and using the complete cylindrical sample with the new method gave still more accurate and precise results. Fortran 90
Inter-station coda wavefield studies using a novel icequake database on Erebus volcano
NASA Astrophysics Data System (ADS)
Chaput, J. A.; Campillo, M.; Roux, P.; Aster, R. C.
2013-12-01
Recent theoretical advances pertaining to the properties of multiply scattered wavefields have yielded a plethora of numerical and controlled source studies aiming to better understand what information may be derived from these otherwise chaotic signals. Practically, multiply scattered wavefields are difficult to compare to numerically derived models due to a combination of source paucity/directionality and array density limitations, particularly in passive seismology scenarios. Furthermore, in situations where data quantities are abundant, such as for ambient noise correlations, it remains very difficult to recover pseudo-Green's function symmetry in the ballistic components of the wavefield, let alone in the coda of the correlations. In this study, we use a large network of short period and broadband instruments on Erebus volcano to show that actual Green's function recovery is indeed possible in some cases. We make use of a large database of small impulsive icequakes distributed randomly on the summit plateau and, using fundamental theoretical properties of equipartitioned wavefields and interstation icequake coda correlations, are able to directly derive notoriously difficult quantities such as the bulk elastic mean free path for the volcano, demonstrations of correlation coda symmetry and its dependence on the number of icequakes used, and a theoretically predicted coherent backscattering amplification factor associated with weak localization. We furthermore show that stable equipartition and H^2/V^2 ratios may be consistently observed for icequake coda, and we perform simple depth inversions of these frequency dependent quantities to compare with known structures.
NASA Technical Reports Server (NTRS)
Levy, R.
1991-01-01
Post-processing algorithms are given to compute the vibratory elastic-rigid coupling matrices and the modal contributions to the rigid-body mass matrices and to the effective modal inertias and masses. Recomputation of the elastic-rigid coupling matrices for a change in origin is also described. A computational example is included. The algorithms can all be executed by using standard finite-element program eigenvalue analysis output with no changes to existing code or source programs.
NASA Astrophysics Data System (ADS)
Nguyen, Bao D.
Imaging with prestack reverse-time migration (RTM) is typically approached via a zero-lag crosscorrelation between source and receiver wavefields, which imposes unnecessarily stringent requirements for computational resources and disk storage. The imaging principle for reflectivity is analyzed and we demonstrate that a single maximal energy arrival event is often sufficient for migration imaging. Methods to alleviate the cost of crosscorrelation imaging are proposed and categorized into reconstructive and non-reconstructive schemes. Source wavefield reconstruction treats the source extrapolation as a method of providing the auxiliary conditions for an initial-boundary value problem. A first-pass (forward-time) extrapolation for the source wavefield identifies the boundary and/or initial values necessary to uniquely reconstruct it using a second (reverse-time) backward propagation. Mixed value, or hybrid, reconstruction is proposed as the most accurate alternative to storing the source wavefield time history. Reconstructing the source wavefield reduces storage costs by up to two orders of magnitude without an appreciable loss of image quality. Boundary value and initial value reconstruction methods are extended from acoustic to elastic RTM. Non-reconstructive approaches deviate from the conventional imaging paradigm, as only the most salient information required for imaging is kept. A maximal energy arrival event (termed the `excitation amplitude') imaging condition is explored as the direct analog for the theoretical reflection coefficient for acoustic isotropic media, and extended for elastic RTM. Sparse crosscorrelation is proposed as an equivalent method to standard crosscorrelation where the migrated image is now represented with a minimized data set. Time-binning is dynamic sorting algorithm with linear time complexity proposed for use with both excitation amplitude and sparse crosscorrelation approches to further expedite imaging. These parsimonious imaging
NASA Astrophysics Data System (ADS)
Maeda, Takuto; Nishida, Kiwamu; Takagi, Ryota; Obara, Kazushige
2016-04-01
The high-sensitive seismograph network (Hi-net) operated by National Research Institute for Earth Science and Disaster Prevention (NIED) has about 800 stations with average separation of 20 km all over the Japanese archipelago. Although it is equipped with short-period seismometers, we also can observe long-period seismic wave up to 100 s in periods for significantly large earthquakes. In this case, we may treat long-period seismic waves as a 2D wavefield with station separations shorter than wavelength rather than individual traces at stations. In this study, we attempt to reconstruct 2D wavefield and obtain its propagation properties from seismic gradiometry (SG) method. The SG estimates the wave amplitude and its spatial derivative coefficients from discrete station record by the Taylor series approximation with an inverse problem. By using spatial derivatives in horizontal directions, we can obtain properties of propagating wave packet such as the arrival direction, slowness, geometrical spreading and radiation pattern. In addition, by using spatial derivatives together with free-surface boundary condition, we may decompose the vector elastic 2D wavefield estimated by the SG into divergence and rotation components. First, we applied the seismic gradiometry to a synthetic long-period (20-50 s) seismogram dataset computed by numerical simulation in realistic 3D medium at the Hi-net station layout as a feasibility test. We confirmed that the wave amplitude and its spatial derivatives are very well reproduced with average correlation coefficients higher than 0.99 in this period range. Applications to a real large earthquakes show that the amplitude and phase of the wavefield are well reconstructed with additional information of arrival direction and its slowness. The reconstructed wavefield contained a clear contrast in slowness between body and surface waves, regional non-great-circle-path wave propagation which may be attributed to scattering. Slowness
Modeling scattering from azimuthally symmetric bathymetric features using wavefield superposition.
Fawcett, John A
2007-12-01
In this paper, an approach for modeling the scattering from azimuthally symmetric bathymetric features is described. These features are useful models for small mounds and indentations on the seafloor at high frequencies and seamounts, shoals, and basins at low frequencies. A bathymetric feature can be considered as a compact closed region, with the same sound speed and density as one of the surrounding media. Using this approach, a number of numerical methods appropriate for a partially buried target or facet problem can be applied. This paper considers the use of wavefield superposition and because of the azimuthal symmetry, the three-dimensional solution to the scattering problem can be expressed as a Fourier sum of solutions to a set of two-dimensional scattering problems. In the case where the surrounding two half spaces have only a density contrast, a semianalytic coupled mode solution is derived. This provides a benchmark solution to scattering from a class of penetrable hemispherical bosses or indentations. The details and problems of the numerical implementation of the wavefield superposition method are described. Example computations using the method for a simple scattering feature on a seabed are presented for a wide band of frequencies. PMID:18247740
Multiscale simulation of 2D elastic wave propagation
NASA Astrophysics Data System (ADS)
Zhang, Wensheng; Zheng, Hui
2016-06-01
In this paper, we develop the multiscale method for simulation of elastic wave propagation. Based on the first-order velocity-stress hyperbolic form of 2D elastic wave equation, the particle velocities are solved first ona coarse grid by the finite volume method. Then the stress tensor is solved by using the multiscale basis functions which can represent the fine-scale variation of the wavefield on the coarse grid. The basis functions are computed by solving a local problem with the finite element method. The theoretical formulae and description of the multiscale method for elastic wave equation are given in more detail. The numerical computations for an inhomogeneous model with random scatter are completed. The results show the effectiveness of the multiscale method.
NASA Astrophysics Data System (ADS)
Arnoux, A.; Batou, A.; Soize, C.; Gagliardini, L.
2013-08-01
This paper is devoted to the construction of a stochastic reduced order computational model of structures having numerous local elastic modes in low frequency dynamics. We are particularly interested in automotive vehicles which are made up of stiff parts and flexible components. This type of structure is characterized by the fact that it exhibits, in the low frequency range, not only the classical global elastic modes but also numerous local elastic modes which cannot easily be separated from the global elastic modes. To solve this difficult problem, an innovative method has recently been proposed for constructing a reduced order computational dynamical model adapted to this particular situation for the low frequency range. Then a new adapted generalized eigenvalue problem is introduced and allows a global vector basis to be constructed for the global displacements space. This method requires to decompose the domain of the structure into sub-domains. Such a decomposition is carried out using the Fast Marching Method. This global vector basis is then used to construct the reduced order computational model. Since there are model uncertainties induced by modeling errors in the computational model, the nonparametric probabilistic approach of uncertainties is used and implemented in the reduced order computational model. The methodology is applied to a complex computational model of an automotive vehicle.
Adaptive finite difference for seismic wavefield modelling in acoustic media.
Yao, Gang; Wu, Di; Debens, Henry Alexander
2016-01-01
Efficient numerical seismic wavefield modelling is a key component of modern seismic imaging techniques, such as reverse-time migration and full-waveform inversion. Finite difference methods are perhaps the most widely used numerical approach for forward modelling, and here we introduce a novel scheme for implementing finite difference by introducing a time-to-space wavelet mapping. Finite difference coefficients are then computed by minimising the difference between the spatial derivatives of the mapped wavelet and the finite difference operator over all propagation angles. Since the coefficients vary adaptively with different velocities and source wavelet bandwidths, the method is capable to maximise the accuracy of the finite difference operator. Numerical examples demonstrate that this method is superior to standard finite difference methods, while comparable to Zhang's optimised finite difference scheme. PMID:27491333
Adaptive finite difference for seismic wavefield modelling in acoustic media
NASA Astrophysics Data System (ADS)
Yao, Gang; Wu, Di; Debens, Henry Alexander
2016-08-01
Efficient numerical seismic wavefield modelling is a key component of modern seismic imaging techniques, such as reverse-time migration and full-waveform inversion. Finite difference methods are perhaps the most widely used numerical approach for forward modelling, and here we introduce a novel scheme for implementing finite difference by introducing a time-to-space wavelet mapping. Finite difference coefficients are then computed by minimising the difference between the spatial derivatives of the mapped wavelet and the finite difference operator over all propagation angles. Since the coefficients vary adaptively with different velocities and source wavelet bandwidths, the method is capable to maximise the accuracy of the finite difference operator. Numerical examples demonstrate that this method is superior to standard finite difference methods, while comparable to Zhang’s optimised finite difference scheme.
Adaptive finite difference for seismic wavefield modelling in acoustic media
Yao, Gang; Wu, Di; Debens, Henry Alexander
2016-01-01
Efficient numerical seismic wavefield modelling is a key component of modern seismic imaging techniques, such as reverse-time migration and full-waveform inversion. Finite difference methods are perhaps the most widely used numerical approach for forward modelling, and here we introduce a novel scheme for implementing finite difference by introducing a time-to-space wavelet mapping. Finite difference coefficients are then computed by minimising the difference between the spatial derivatives of the mapped wavelet and the finite difference operator over all propagation angles. Since the coefficients vary adaptively with different velocities and source wavelet bandwidths, the method is capable to maximise the accuracy of the finite difference operator. Numerical examples demonstrate that this method is superior to standard finite difference methods, while comparable to Zhang’s optimised finite difference scheme. PMID:27491333
Managing competing elastic Grid and Cloud scientific computing applications using OpenNebula
NASA Astrophysics Data System (ADS)
Bagnasco, S.; Berzano, D.; Lusso, S.; Masera, M.; Vallero, S.
2015-12-01
Elastic cloud computing applications, i.e. applications that automatically scale according to computing needs, work on the ideal assumption of infinite resources. While large public cloud infrastructures may be a reasonable approximation of this condition, scientific computing centres like WLCG Grid sites usually work in a saturated regime, in which applications compete for scarce resources through queues, priorities and scheduling policies, and keeping a fraction of the computing cores idle to allow for headroom is usually not an option. In our particular environment one of the applications (a WLCG Tier-2 Grid site) is much larger than all the others and cannot autoscale easily. Nevertheless, other smaller applications can benefit of automatic elasticity; the implementation of this property in our infrastructure, based on the OpenNebula cloud stack, will be described and the very first operational experiences with a small number of strategies for timely allocation and release of resources will be discussed.
NASA Astrophysics Data System (ADS)
Liu, B.; Arsenlis, A.; Aubry, S.
2016-06-01
Driven by the growing interest in numerical simulations of dislocation–interface interactions in general crystalline materials with elastic anisotropy, we develop algorithms for the integration of interface tractions needed to couple dislocation dynamics with a finite element or boundary element solver. The dislocation stress fields in elastically anisotropic media are made analytically accessible through the spherical harmonics expansion of the derivative of Green’s function, and analytical expressions for the forces on interface elements are derived by analytically integrating the spherical harmonics series recursively. Compared with numerical integration by Gaussian quadrature, the newly developed analytical algorithm for interface traction integration is highly beneficial in terms of both computation precision and speed.
Numerical Modeling of Ocean Acoustic Wavefields
NASA Astrophysics Data System (ADS)
Tappert, Frederick
1997-08-01
The U.S. Navy requires real-time ``acoustic performance prediction'' models in order to optimize sonar tactics in naval combat situations. The need for numerical models that solve the acoustic wave equation in realistic ocean environments is being met by a collaborative effort between university researchers, industrial contractors, and navy laboratory workers. This paper discusses one particularly successful numerical model, called the PE/SSF model, that was originally developed by the author. Here PE stands for Parabolic Equation, a good approximation to the elliptic Helmholtz equation; and SSF stands for the Split-Step Fourier algorithm, a highly efficient marching algorithm for solving parabolic type equations. These techniques are analyzed, and examples are displayed of ocean acoustic wavefields generated by the PE/SSF model.
Isolation of ultrasonic scattering by wavefield baseline subtraction
NASA Astrophysics Data System (ADS)
Dawson, Alexander J.; Michaels, Jennifer E.; Michaels, Thomas E.
2016-03-01
Wavefield imaging generally refers to the measurement of signals over a two-dimensional rectilinear grid that originate from a spatially fixed source. Subtraction of such wavefields is investigated as a means of separating scattered signals from the total wavefield; that is, baseline wavefield data acquired from a defect-free specimen are subtracted from analogous data acquired after introduction of a defect. The wavefields considered here are generated by a 5 MHz angle-beam probe and measured over an area of the accessible specimen surface using a laser vibrometer. The primary challenge in isolating the scattered waves is imperfect temporal and spatial alignment of the two wavefields. Two methods for aligning the wavefields in space and time prior to performing baseline subtraction are presented and their efficacy is evaluated using data acquired before and after introducing notches that originate from a through-hole. Although perfect baseline subtraction is not achieved, the improvement in performance after alignment using either method allows for scattered waves from small defects to be separated and visualized, even when their amplitudes are much smaller than those of the incident waves.
Elastic Cloud Computing Infrastructures in the Open Cirrus Testbed Implemented via Eucalyptus
NASA Astrophysics Data System (ADS)
Baun, Christian; Kunze, Marcel
Cloud computing realizes the advantages and overcomes some restrictionsof the grid computing paradigm. Elastic infrastructures can easily be createdand managed by cloud users. In order to accelerate the research ondata center management and cloud services the OpenCirrusTM researchtestbed has been started by HP, Intel and Yahoo!. Although commercialcloud offerings are proprietary, Open Source solutions exist in the field ofIaaS with Eucalyptus, PaaS with AppScale and at the applications layerwith Hadoop MapReduce. This paper examines the I/O performance ofcloud computing infrastructures implemented with Eucalyptus in contrastto Amazon S3.
NASA Technical Reports Server (NTRS)
Poole, L. R.
1972-01-01
A computer program is presented by which the effects of nonlinear suspension-system elastic characteristics on parachute inflation loads and motions can be investigated. A mathematical elastic model of suspension-system geometry is coupled to the planar equations of motion of a general vehicle and canopy. Canopy geometry and aerodynamic drag characteristics and suspension-system elastic properties are tabular inputs. The equations of motion are numerically integrated by use of an equivalent fifth-order Runge-Kutta technique.
Numerical Resolution of Seismic Wavefield Simulations in Southern California
NASA Astrophysics Data System (ADS)
Tape, C.; Casarotti, E.
2014-12-01
Seismic wavefield simulations can provide accurate solutions to the wave equation, even for three-dimensional seismic velocity models with topography, basin structures, anisotropy, attenuation, and other complexities. But how accurate are these numerical solutions? In many cases the effect of numerical dispersion on the synthetic seismograms looks quite similar to the effects of structural complexities. Therefore it is important to know the numerical resolution of the synthetic seismograms, that is, the minimum period that provides a quantifiably numerically accurate solution to the wave equation. Numerical resolution can be discussed in terms of a combined mesh and velocity model, or in terms of an individual source-station path within the same mesh and velocity model. Here we discuss two approaches for quantifying the numerical resolution. In the first approach we estimate the minimum resolvable period of each element within the finite element mesh. This calculation, performed within SPECFEM3D, requires no wavefield simulations and is based only on the size of each element and the minimum velocity within each element. The calculation produces a volumetric field that shows the estimated minimum resolvable period within each element of the (generally unstructured) mesh. In the second approach we choose a set of earthquakes to be used within a tomographic inversion. For each earthquake we compute one simulation using a fine discretization of gridpoints on the finite-element mesh and another simulation using a coarse discretization. We filter both sets of seismograms over a range of periods and then quantify the waveform differences. The minimum resolvable period (or numerical resolution) is identified by the minimum period for which the synthetic seismogram from the coarse-mesh and fine-mesh simulations is quantifiably the same. This calculation provides a path-specific minimum resolvable period that can be used to guide the choice of measurement filtering for a
Computing dislocation stress fields in anisotropic elastic media using fast multipole expansions
NASA Astrophysics Data System (ADS)
Yin, Jie; Barnett, D. M.; Fitzgerald, S. P.; Cai, Wei
2012-06-01
The calculation of stress fields due to dislocations and hence the forces they exert on each other is the most time consuming step in dislocation dynamics (DD) simulations. The fast multipole method (FMM) can reduce the computational cost at each simulation step from { O}(N^2) to { O}(N) for an ensemble of N dislocation segments. However, FMM has not yet been applied to three-dimensional DD simulations which take into account anisotropic elasticity. We demonstrate a systematic procedure to establish this capability by first obtaining the derivatives of the elastic Green's function to arbitrary order for a medium of general anisotropy. We then compute the stress field of a dislocation loop using multipole expansions based on these derivatives, and analyze the dependence of numerical errors on the expansion order. This method can be implemented in large scale DD simulations when the consideration of elastic anisotropy is necessary, for example the technologically important cases of iron and ferritic steels at high temperatures.
Computation of Temperature-Dependent Legendre Moments of a Double-Differential Elastic Cross Section
Arbanas, Goran; Dunn, Michael E; Larson, Nancy M; Leal, Luiz C; Williams, Mark L; Becker, B.; Dagan, R
2011-01-01
A general expression for temperature-dependent Legendre moments of a double-differential elastic scattering cross section was derived by Ouisloumen and Sanchez [Nucl. Sci. Eng. 107, 189-200 (1991)]. Attempts to compute this expression are hindered by the three-fold nested integral, limiting their practical application to just the zeroth Legendre moment of an isotropic scattering. It is shown that the two innermost integrals could be evaluated analytically to all orders of Legendre moments, and for anisotropic scattering, by a recursive application of the integration by parts method. For this method to work, the anisotropic angular distribution in the center of mass is expressed as an expansion in Legendre polynomials. The first several Legendre moments of elastic scattering of neutrons on U-238 are computed at T=1000 K at incoming energy 6.5 eV for isotropic scattering in the center of mass frame. Legendre moments of the anisotropic angular distribution given via Blatt-Biedenharn coefficients are computed at ~1 keV. The results are in agreement with those computed by the Monte Carlo method.
A computational framework for polyconvex large strain elasticity for geometrically exact beam theory
NASA Astrophysics Data System (ADS)
Ortigosa, Rogelio; Gil, Antonio J.; Bonet, Javier; Hesch, Christian
2016-02-01
In this paper, a new computational framework is presented for the analysis of nonlinear beam finite elements subjected to large strains. Specifically, the methodology recently introduced in Bonet et al. (Comput Methods Appl Mech Eng 283:1061-1094, 2015) in the context of three dimensional polyconvex elasticity is extended to the geometrically exact beam model of Simo (Comput Methods Appl Mech Eng 49:55-70, 1985), the starting point of so many other finite element beam type formulations. This new variational framework can be viewed as a continuum degenerate formulation which, moreover, is enhanced by three key novelties. First, in order to facilitate the implementation of the sophisticated polyconvex constitutive laws particularly associated with beams undergoing large strains, a novel tensor cross product algebra by Bonet et al. (Comput Methods Appl Mech Eng 283:1061-1094, 2015) is adopted, leading to an elegant and physically meaningful representation of an otherwise complex computational framework. Second, the paper shows how the novel algebra facilitates the re-expression of any invariant of the deformation gradient, its cofactor and its determinant in terms of the classical beam strain measures. The latter being very useful whenever a classical beam implementation is preferred. This is particularised for the case of a Mooney-Rivlin model although the technique can be straightforwardly generalised to other more complex isotropic and anisotropic polyconvex models. Third, the connection between the two most accepted restrictions for the definition of constitutive models in three dimensional elasticity and beams is shown, bridging the gap between the continuum and its degenerate beam description. This is carried out via a novel insightful representation of the tangent operator.
Local guided wavefield analysis for characterization of delaminations in composites
NASA Astrophysics Data System (ADS)
Rogge, M. D.; Leckey, C. A. C.
2013-01-01
Delaminations in composite laminates resulting from impact events may be accompanied by minimal indication of damage at the surface. As such, inspection techniques are required to ensure defects are within allowable limits. Conventional ultrasonic scanning techniques have been shown to effectively characterize the size and depth of delaminations but require physical contact with the structure. Alternatively, a noncontact scanning laser vibrometer may be used to measure guided wave propagation in the laminate structure. A local Fourier domain analysis method is presented for processing guided wavefield data to estimate spatially-dependent wavenumber values, which can be used to determine delamination depth. The technique is applied to simulated wavefields and results are analyzed to determine limitations of the technique with regards to determining defect size and depth. Finally, experimental wavefield data obtained in quasi-isotropic carbon fiber reinforced polymer (CFRP) laminates with impact damage is analyzed and wavenumber is measured to an accuracy of 8.5% in the region of shallow delaminations.
Elasticity and phase stability of pyrope garnet from ab initio computation
NASA Astrophysics Data System (ADS)
Kawai, Kenji; Tsuchiya, Taku
2015-03-01
We study the high-pressure stability and elastic properties of Mg3Al2Si3O12 pyrope garnet using the density functional first principles computation method. Pyrope garnet is found to dissociate into an assemblage of MgSiO3 Mg-perovskite (Pv) and Al2O3 corundum (Cor) solid solutions at ∼19.7 GPa at static conditions. Then, this assemblage undergoes a phase transition to pyropic (Al-bearing) Pv at ∼65 GPa. The enthalpy of an assemblage of MgAl2O4 calcium ferrite (CF), MgPv, and stishovite (St) is less stable than that of MgPv plus Cor. A continuous reaction in the MgSiO3-Al2O3 system suggested by this study is consistent with previous experimental and computational studies but not with a recently modeled phase diagram. This implies that the formation of pyropic Pv could not cause any seismic scatterers in the mid-lower mantle. The investigated anisotropy of elastic velocities further indicates that pyrope garnet is a very isotropic mineral. Our results suggest that seismological anisotropy inferred in the upper mantle could not be due to garnet.
Three-dimensional full-wavefield seismic tomography on field data (Invited)
NASA Astrophysics Data System (ADS)
Warner, M.; Umpleby, A.; Stekl, I.; Guasch, L.
2010-12-01
In contrast to conventional seismic tomography, where we minimise the mismatch between observed and calculated seismic travel times, in full-wavefield tomography we seek a model that is able to match the entire observed wavefield, wiggle-for-wiggle. Wavefield tomography has a long history, but it is only recently that advances in algorithms and in hardware have made the technique feasible on realistic-sized datasets in three dimensions. With sponsorship from the petroleum industry, we have developed full 3D codes for anisotropic acoustic, and isotropic elastic, wavefield tomography in the time-domain, and for visco-acoustic tomography in the frequency domain. In both domains, we solve the wave equation using finite differences on a regular mesh; we use explicit time-stepping in the time domain, and use an implicit iterative solver in the frequency domain. The codes are parallelised to run on a cluster of multi-core nodes, and they are able to deal with large irregular 3D datasets efficiently. We report here the results of applying these codes to a 3D ocean-bottom seismic dataset acquired over the Tommeliten oil field in the North Sea. The field data are composed of 1920 four-component ocean-bottom receivers, recording about 30,000 air-gun sources over an area of 12 x 9 km. A low-velocity, high-attenuation gas cloud is located at a depth of 1 to 2 km; this gas cloud partially obscures the geology of the underlying oil field. There is significant anisotropy within the section; vertical and horizontal p-wave velocities can differ by more than 15%. Wavefield tomography is successful in imaging the complex velocity structure in 3D within this gas cloud with a lateral resolution of about 25 m. This resolution is much better than that obtained using reflection travel-time tomography or migration velocity analysis. Subsequent pre-stack reverse-time depth migration of the underlying reflection data demonstrates that the recovered velocity structure is real. It is necessary
GPU computing with OpenCL to model 2D elastic wave propagation: exploring memory usage
NASA Astrophysics Data System (ADS)
Iturrarán-Viveros, Ursula; Molero-Armenta, Miguel
2015-01-01
Graphics processing units (GPUs) have become increasingly powerful in recent years. Programs exploring the advantages of this architecture could achieve large performance gains and this is the aim of new initiatives in high performance computing. The objective of this work is to develop an efficient tool to model 2D elastic wave propagation on parallel computing devices. To this end, we implement the elastodynamic finite integration technique, using the industry open standard open computing language (OpenCL) for cross-platform, parallel programming of modern processors, and an open-source toolkit called [Py]OpenCL. The code written with [Py]OpenCL can run on a wide variety of platforms; it can be used on AMD or NVIDIA GPUs as well as classical multicore CPUs, adapting to the underlying architecture. Our main contribution is its implementation with local and global memory and the performance analysis using five different computing devices (including Kepler, one of the fastest and most efficient high performance computing technologies) with various operating systems.
Signal and image processing algorithm performance in a virtual and elastic computing environment
NASA Astrophysics Data System (ADS)
Bennett, Kelly W.; Robertson, James
2013-05-01
The U.S. Army Research Laboratory (ARL) supports the development of classification, detection, tracking, and localization algorithms using multiple sensing modalities including acoustic, seismic, E-field, magnetic field, PIR, and visual and IR imaging. Multimodal sensors collect large amounts of data in support of algorithm development. The resulting large amount of data, and their associated high-performance computing needs, increases and challenges existing computing infrastructures. Purchasing computer power as a commodity using a Cloud service offers low-cost, pay-as-you-go pricing models, scalability, and elasticity that may provide solutions to develop and optimize algorithms without having to procure additional hardware and resources. This paper provides a detailed look at using a commercial cloud service provider, such as Amazon Web Services (AWS), to develop and deploy simple signal and image processing algorithms in a cloud and run the algorithms on a large set of data archived in the ARL Multimodal Signatures Database (MMSDB). Analytical results will provide performance comparisons with existing infrastructure. A discussion on using cloud computing with government data will discuss best security practices that exist within cloud services, such as AWS.
Impact Induced Delamination Detection and Quantification With Guided Wavefield Analysis
NASA Technical Reports Server (NTRS)
Tian, Zhenhua; Leckey, Cara A. C.; Yu, Lingyu; Seebo, Jeffrey P.
2015-01-01
This paper studies impact induced delamination detection and quantification by using guided wavefield data and spatial wavenumber imaging. The complex geometry impact-like delamination is created through a quasi-static indentation on a CFRP plate. To detect and quantify the impact delamination in the CFRP plate, PZT-SLDV sensing and spatial wavenumber imaging are performed. In the PZT-SLDV sensing, the guided waves are generated from the PZT, and the high spatial resolution guided wavefields are measured by the SLDV. The guided wavefield data acquired from the PZT-SLDV sensing represent guided wave propagation in the composite laminate and include guided wave interaction with the delamination damage. The measured guided wavefields are analyzed through the spatial wavenumber imaging method, which generates an image containing the dominant local wavenumber at each spatial location. The spatial wavenumber imaging result for the simple single layer Teflon insert delamination provided quantitative information on delamination damage size and location. The location of delamination damage is indicated by the area with larger wavenumbers in the spatial wavenumber image. The impact-like delamination results only partially agreed with the damage size and shape. The results also demonstrated the dependence on excitation frequency. Future work will further investigate the accuracy of the wavenumber imaging method for real composite damage and the dependence on frequency of excitation.
Local Guided Wavefield Analysis for Characterization of Delaminations in Composites
NASA Technical Reports Server (NTRS)
Rogge, Matthew D.; Campbell Leckey, Cara A.
2012-01-01
Delaminations in composite laminates resulting from impact events may be accompanied by minimal indication of damage at the surface. As such, inspection techniques are required to ensure defects are within allowable limits. Conventional ultrasonic scanning techniques have been shown to effectively characterize the size and depth of delaminations but require physical contact with the structure. Alternatively, a noncontact scanning laser vibrometer may be used to measure guided wave propagation in the laminate structure. A local Fourier domain analysis method is presented for processing guided wavefield data to estimate spatially-dependent wavenumber values, which can be used to determine delamination depth. The technique is applied to simulated wavefields and results are analyzed to determine limitations of the technique with regards to determining defect size and depth. Finally, experimental wavefield data obtained in quasi-isotropic carbon fiber reinforced polymer (CFRP) laminates with impact damage is analyzed and wavenumber is measured to an accuracy of 8.5% in the region of shallow delaminations. Keywords: Ultrasonic wavefield imaging, Windowed Fourier transforms, Guided waves, Structural health monitoring, Nondestructive evaluation
Seismic interferometry by multidimensional deconvolution without wavefield separation
NASA Astrophysics Data System (ADS)
Ravasi, Matteo; Meles, Giovanni; Curtis, Andrew; Rawlinson, Zara; Yikuo, Liu
2015-07-01
Seismic interferometry comprises a suite of methods to redatum recorded wavefields to those that would have been recorded if different sources (so-called virtual sources) had been activated. Seismic interferometry by cross-correlation has been formulated using either two-way (for full wavefields) or one-way (for directionally decomposed wavefields) representation theorems. To obtain improved Green's function estimates, the cross-correlation result can be deconvolved by a quantity that identifies the smearing of the virtual source in space and time, the so-called point-spread function. This type of interferometry, known as interferometry by multidimensional deconvolution (MDD), has so far been applied only to one-way directionally decomposed fields, requiring accurate wavefield decomposition from dual (e.g. pressure and velocity) recordings. Here we propose a form of interferometry by multidimensional deconvolution that uses full wavefields with two-way representations, and simultaneously invert for pressure and (normal) velocity Green's functions, rather than only velocity responses as for its one-way counterpart. Tests on synthetic data show that two-way MDD improves on results of interferometry by cross-correlation, and generally produces estimates of similar quality to those obtained by one-way MDD, suggesting that the preliminary decomposition into up- and downgoing components of the pressure field is not required if pressure and velocity data are jointly used in the deconvolution. We also show that constraints on the directionality of the Green's functions sought can be added directly into the MDD inversion process to further improve two-way multidimensional deconvolution. Finally, as a by-product of having pressure and particle velocity measurements, we adapt one- and two-way representation theorems to convert any particle velocity receiver into its corresponding virtual dipole/gradient source by means of MDD. Thus data recorded from standard monopolar (e
Barkaoui, Abdelwahed; Chamekh, Abdessalem; Merzouki, Tarek; Hambli, Ridha; Mkaddem, Ali
2014-03-01
The complexity and heterogeneity of bone tissue require a multiscale modeling to understand its mechanical behavior and its remodeling mechanisms. In this paper, a novel multiscale hierarchical approach including microfibril scale based on hybrid neural network (NN) computation and homogenization equations was developed to link nanoscopic and macroscopic scales to estimate the elastic properties of human cortical bone. The multiscale model is divided into three main phases: (i) in step 0, the elastic constants of collagen-water and mineral-water composites are calculated by averaging the upper and lower Hill bounds; (ii) in step 1, the elastic properties of the collagen microfibril are computed using a trained NN simulation. Finite element calculation is performed at nanoscopic levels to provide a database to train an in-house NN program; and (iii) in steps 2-10 from fibril to continuum cortical bone tissue, homogenization equations are used to perform the computation at the higher scales. The NN outputs (elastic properties of the microfibril) are used as inputs for the homogenization computation to determine the properties of mineralized collagen fibril. The mechanical and geometrical properties of bone constituents (mineral, collagen, and cross-links) as well as the porosity were taken in consideration. This paper aims to predict analytically the effective elastic constants of cortical bone by modeling its elastic response at these different scales, ranging from the nanostructural to mesostructural levels. Our findings of the lowest scale's output were well integrated with the other higher levels and serve as inputs for the next higher scale modeling. Good agreement was obtained between our predicted results and literature data. PMID:24123969
AxiSEM and instaseis: Fast simulation of global wavefields across the frequency band
NASA Astrophysics Data System (ADS)
Nissen-Meyer, T.; van Driel, M.; Krischer, L.; Stähler, S. C.; Hosseini, K.; Leng, K.
2015-12-01
We present our seismic modeling methods AxiSEM and instaseis. These methods exploit recent developments in high-performance computing and suitable numerical methods for seismic wave propagation, while operating efficiently across the vast observable frequency spectrum of global waves in sparse yet realistic structures. AxiSEM (www.axisem.info and geodynamics.org) relies upon axisymmetric (including spherically symmetric) models, thereby satisfying a large fraction of observable data. The benefit of this method lies in the resultant dimensional collapse to two numerical dimensions, whereby the third azimuthal dimension is tackled analytically. For high-frequency wave propagation, this leads to 3-4 orders of magnitude speedup in computational cost compared to 3D domain discretizations. AxiSEM is highly scalable and accommodates efficient implementations of viscoelasticity and anisotropy. We will present benchmarks, data comparisons, a diverse range of applications from inner-core anisotropy to noise modeling and lowermost mantle structures, and wavefields for sensitivity kernels. We also touch upon ongoing efforts for linking computational cost to structural complexity in the vein of Occam's razor, eventually allowing for an adaptive rendition of 1D, 2D and 3D structures at optimally low computational cost, as well as 1D/3D hybrid approaches. Instaseis (www.instaseis.net) is a methodology to extract full, broadband and accurate waveforms instantaneously from wavefield databases computed with AxiSEM. This "once-and-for-all solution" relies on reciprocity and requires only two AxiSEM simulations to construct the databases, while allowing for arbitrary parameter changes (e.g. source, processing, structure) instantaneously with modest computational cost and storage requirements. The instaseis python package is integrated with ObsPy, contains a graphical user interface, and can be used for source inversion, noise simulations, finite-fault modeling, waveform tomography
NASA Astrophysics Data System (ADS)
Sidler, Rolf; Rubino, J. Germán; Holliger, Klaus
2013-04-01
There is increasing evidence to suggest that the presence of mesoscopic heterogeneities constitutes an important seismic attenuation mechanism in porous rocks. As a consequence, centimetre-scale perturbations of the rock physical properties should be taken into account for seismic modelling whenever detailed and accurate responses of specific target structures are desired, which is, however, computationally prohibitive. A convenient way to circumvent this problem is to use an upscaling procedure to replace each of the heterogeneous porous media composing the geological model by corresponding equivalent visco-elastic solids and to solve the visco-elastic equations of motion for the inferred equivalent model. While the overall qualitative validity of this procedure is well established, there are as of yet no quantitative analyses regarding the equivalence of the seismograms resulting from the original poro-elastic and the corresponding upscaled visco-elastic models. To address this issue, we compare poro-elastic and visco-elastic solutions for a range of marine-type models of increasing complexity. We found that despite the identical dispersion and attenuation behaviour of the heterogeneous poro-elastic and the equivalent visco-elastic media, the seismograms may differ substantially due to diverging boundary conditions, where there exist additional options for the poro-elastic case. In particular, we observe that at the fluid/porous-solid interface, the poro- and visco-elastic seismograms agree for closed-pore boundary conditions, but differ significantly for open-pore boundary conditions. This is an important result which has potentially far-reaching implications for wave-equation-based algorithms in exploration geophysics involving fluid/porous-solid interfaces, such as, for example, wavefield decomposition.
Guided Wave Delamination Detection and Quantification With Wavefield Data Analysis
NASA Technical Reports Server (NTRS)
Tian, Zhenhua; Campbell Leckey, Cara A.; Seebo, Jeffrey P.; Yu, Lingyu
2014-01-01
Unexpected damage can occur in aerospace composites due to impact events or material stress during off-nominal loading events. In particular, laminated composites are susceptible to delamination damage due to weak transverse tensile and inter-laminar shear strengths. Developments of reliable and quantitative techniques to detect delamination damage in laminated composites are imperative for safe and functional optimally-designed next-generation composite structures. In this paper, we investigate guided wave interactions with delamination damage and develop quantification algorithms by using wavefield data analysis. The trapped guided waves in the delamination region are observed from the wavefield data and further quantitatively interpreted by using different wavenumber analysis methods. The frequency-wavenumber representation of the wavefield shows that new wavenumbers are present and correlate to trapped waves in the damage region. These new wavenumbers are used to detect and quantify the delamination damage through the wavenumber analysis, which can show how the wavenumber changes as a function of wave propagation distance. The location and spatial duration of the new wavenumbers can be identified, providing a useful means not only for detecting the presence of delamination damage but also allowing for estimation of the delamination size. Our method has been applied to detect and quantify real delamination damage with complex geometry (grown using a quasi-static indentation technique). The detection and quantification results show the location, size, and shape of the delamination damage.
Polarity-consistent excitation amplitude imaging condition for elastic reverse time migration
NASA Astrophysics Data System (ADS)
Du, Qizhen; Zhang, Mingqiang; Gong, Xufei; Chen, Xiaoran
2015-02-01
Imaging conditions can not only affect the computational efficiency and storage cost of reverse time migration (RTM) but determine the quality of the final migrated images. This paper extends the idea of the well amplitude-preserved and highly-efficient excitation amplitude imaging condition from acoustic RTM to elastic RTM. For elastic RTM, the maximum amplitude of the separated P-wave and the corresponding image time of each grid point are saved during the forward modeling of the source wavefield and then PP and PS images are obtained by dividing the separated P- and S-waves of the backward-propagated receiver wavefield by the precomputed P-waves at each grid point that satisfies the image time. However, polarity reversals of the PS image will cause destructive interference when the stacked image is needed. In order to solve this problem, we propose the polarity-consistent excitation amplitude imaging condition by combining the excitation amplitude imaging condition with a shot-domain polarity reversal correction method. Then we provide the detailed realization process of this imaging condition in elastic RTM. By utilizing the relatively stable and well amplitude-preserved source-normalized cross-correlation imaging condition as a comparison, we testify to the feasibility and validity of the proposed imaging condition in the aspects of amplitude preservation property, imaging capability of complex structures, storage cost and computational efficiency. Considering the balance between the efficiency and image quality, the polarity-consistent excitation amplitude imaging condition can be a good choice for elastic RTM.
Computation of elastic properties of 3D digital cores from the Longmaxi shale
NASA Astrophysics Data System (ADS)
Zhang, Wen-Hui; Fu, Li-Yun; Zhang, Yan; Jin, Wei-Jun
2016-06-01
The dependence of elastic moduli of shales on the mineralogy and microstructure of shales is important for the prediction of sweet spots and shale gas production. Based on 3D digital images of the microstructure of Longmaxi black shale samples using X-ray CT, we built detailed 3D digital images of cores with porosity properties and mineral contents. Next, we used finite-element (FE) methods to derive the elastic properties of the samples. The FE method can accurately model the shale mineralogy. Particular attention is paid to the derived elastic properties and their dependence on porosity and kerogen. The elastic moduli generally decrease with increasing porosity and kerogen, and there is a critical porosity (0.75) and kerogen content (ca. ≤3%) over which the elastic moduli decrease rapidly and slowly, respectively. The derived elastic moduli of gas- and oil-saturated digital cores differ little probably because of the low porosity (4.5%) of the Longmaxi black shale. Clearly, the numerical experiments demonstrated the feasibility of combining microstructure images of shale samples with elastic moduli calculations to predict shale properties.
Challenges in the separation and analysis of scattered waves in angle-beam wavefield data
Dawson, Alexander J.; Michaels, Jennifer E.; Michaels, Thomas E.
2015-03-31
The measurement of ultrasonic signals on a 2-D rectilinear grid resulting from a fixed source, referred to as wavefield imaging, is a powerful tool for visualizing wave propagation and scattering. Wavefield imaging provides a more complete picture of wave propagation than conventional single-point measurements, but creates more challenges for analysis. This work considers the development of wavefield-based methods for analyzing angle-beam wave propagation and scattering in plates. Methods of analysis focus on the separation of scattered waves from the total wavefield with the eventual goal of quantitative scatterer characterization in a laboratory environment. Two methods for wave separation are considered: frequency-wavenumber filtering and wavefield baseline subtraction. Frequency-wavenumber filtering is applied to wavefield data that are finely sampled in both space and time, whereas baseline subtraction is a technique that has typically been applied to individual signals recorded from fixed transducers rather than to full wavefield data. Baseline subtraction of wavefields, particularly for the frequency range considered here, is sensitive to both specimen alignment and temperature variations, whereas frequency-wavenumber methods are limited in their ability to separate waves traveling in the same direction. Results are shown for both methods with a focus on investigating and overcoming the challenges to full wavefield baseline subtraction.
Direct and Inverse Problems in Statistical Wavefields
Wolf, Emil
2002-09-01
In this report account is presented of research carried out during the period September 1, 1999-August 31, 2002 under the sponsorship of the Department of Energy, grant DE-FG02-90ER14119. The research covered several areas of modern optical physics, particularly propagation of partially coherent light and its interaction with deterministic and with random media, spectroscopy with partially coherent light, polarization properties of statistical wave fields, effects of moving diffusers on coherence and on the spectra of light transmitted and scattered by them, reciprocity inequalities involving spatial and angular correlations of partially coherent beams, spreading of partially coherent beams in-random media, inverse source problems, computed and diffraction tomography and partially coherent solitons. We have discovered a new phenomenon in an emerging field of physical optics, known as singular optics; specifically we found that the spectrum of light changes drastically in the neighborhood of points where the intensity has zero value and where, consequently, the phase becomes singular, We noted some potential applications of this phenomenon. The results of our investigations were reported in 39 publications. They are listed on pages 3 to 5. Summaries of these publications are given on pages 6-13. Scientists who have participated in this research are listed on page 14.
Ilegbusi, Olusegun; Li, Ziang; Min, Yugang; Meeks, Sanford; Kupelian, Patrick; Santhanam, Anand P
2012-01-01
The aim of this paper is to model the airflow inside lungs during breathing and its fluid-structure interaction with the lung tissues and the lung tumor using subject-specific elastic properties. The fluid-structure interaction technique simultaneously simulates flow within the airway and anisotropic deformation of the lung lobes. The three-dimensional (3D) lung geometry is reconstructed from the end-expiration 3D CT scan datasets of humans with lung cancer. The lung is modeled as a poro-elastic medium with anisotropic elastic property (non-linear Young's modulus) obtained from inverse lung elastography of 4D CT scans for the same patients. The predicted results include the 3D anisotropic lung deformation along with the airflow pattern inside the lungs. The effect is also presented of anisotropic elasticity on both the spatio-temporal volumetric lung displacement and the regional lung hysteresis. PMID:22356987
Imaging of offset VSP data with an elastic iterative migration scheme
Mittet, R.; Hokstad, K.; Helgesen, J.; Canadas, G.
1997-03-01
VSP data are usually acquired in order to obtain high-resolution images of complex structures in reservoirs and near boreholes. The authors present an elastic iterative migration scheme which has few limitations regarding the complexity of the geology, and where the macromodel for both P- and S-wave velocities is automatically improved and updated at each iteration. They avoid wavefield separation (up/down and P/S) and the simplifying assumptions of small dips underlying most such methods. The migration scheme is based on elastic inversion theory. The wavefield extrapolation is based on a high-order, coarse-grid, finite-difference solution to the elastic two-way wave equation. At each iteration, the macromodel is updated using a gradient method, in which the gradient is computed by correlation of forward-modelled fields with back-propagated residual fields. The first iteration of the migration scheme is equivalent to elastic reverse-time migration with an imaging condition similar to Claerbout`s principle. Both P- and S-wave reflections contribute to the images. In numerical examples with both synthetic and real offset VSP data, they find that increasing the number of iterations improves the image quality. Images based on both P- and S-wave energy give more near-well information and higher spatial resolution than images on only acoustic energy. In the real data example the authors show that the iterative migration scheme can image a relatively complex geological structure. A fault and a small graben intersecting the well can be identified.
Frequency-domain elastic full-waveform multiscale inversion method based on dual-level parallelism
NASA Astrophysics Data System (ADS)
Li, Yuan-Yuan; Li, Zhen-Chun; Zhang, Kai; Zhang, Xuan
2015-12-01
The complexity of an elastic wavefield increases the nonlinearity of inversion. To some extent, multiscale inversion decreases the nonlinearity of inversion and prevents it from falling into local extremes. A multiscale strategy based on the simultaneous use of frequency groups and layer stripping method based on damped wave field improves the stability of inversion. A dual-level parallel algorithm is then used to decrease the computational cost and improve practicability. The seismic wave modeling of a single frequency and inversion in a frequency group are computed in parallel by multiple nodes based on multifrontal massively parallel sparse direct solver and MPI. Numerical tests using an overthrust model show that the proposed inversion algorithm can effectively improve the stability and accuracy of inversion by selecting the appropriate inversion frequency and damping factor in lowfrequency seismic data.
Frequency-wavenumber domain analysis of guided wavefields.
Michaels, Thomas E; Michaels, Jennifer E; Ruzzene, Massimo
2011-05-01
Full wavefield measurements obtained with either an air-coupled transducer mounted on a scanning stage or a scanning laser vibrometer can be combined with effective signal and imaging processing algorithms to support characterization of guided waves as well as detection, localization and quantification of structural damage. These wavefield images contain a wealth of information that clearly shows details of guided waves as they propagate outward from the source, reflect from specimen boundaries, and scatter from discontinuities within the structure. The analysis of weaker scattered waves is facilitated by the removal of source waves and the separation of wave modes, which is effectively achieved via frequency-wavenumber domain filtering in conjunction with the subsequent analysis of the resulting residual signals. Incident wave removal highlights the presence and the location of weak scatterers, while the separation of individual guided wave modes allows the characterization of their separate contribution to the scattered field and the evaluation of mode conversion phenomena. The effectiveness of these methods is demonstrated through their application to detection of a delamination in a composite plate and detection of a crack emanating from a hole. PMID:21190706
NASA Astrophysics Data System (ADS)
Petrov, P.; Newman, G. A.
2013-12-01
Full Waveform Inversion (FWI) is a promising seismic imaging method which aims to compute quantitative estimates of the subsurface parameters (bulk wave velocity, shear wave velocity, rock density) from local measurements of the seismic wavefield. It based on a modeling of wave propagation from controlled seismic sources and consists in minimizing iteratively the difference between the predicted wavefields at the receivers and the recorded data. This amounts to solving a strongly nonlinear large-scale inverse problem. We have formulated a 3D inverse solution for the elastic wave problem in Laplace-Fourier domain using the non-linear conjugate gradient methods. FWI in the Laplace-Fourier domain has seen considerable interest over the last several years as an effective approach to imaging seismic data, where conventional FWI schemes have difficulties in converging to acceptable solutions. Finite difference methods, employing staggered grids are used to compute predicted data and objective functional gradients. Typically three forward modeling applications per frequency are required to produce the model update at each iteration. Because the solution's realism and complexity are still limited by the speed and memory of serial processors the code has been implemented on a massive parallel computing platform where hundreds to thousands of processors operate on the problem simultaneously. The inversion scheme is tested by inverting data produced with a forward modelling code for a few models.
Rapid acquisition of high resolution full wave-field borehole seismic data
Sleefe, G.E.; Harding, R.S. Jr.; Fairborn, J.W.; Paulsson, B.N.P.
1993-04-01
An essential requirement for both Vertical Seismic Profiling (VSP) and Cross-Hole Seismic Profiling (CHSP) is the rapid acquisition of high resolution borehole seismic data. Additionally, full wave-field recording using three-component receivers enables the use of both transmitted and reflected elastic wave events in the resulting seismic images of the subsurface. To this end, an advanced three- component multi-station borehole seismic receiver system has been designed and developed by Sandia National Laboratory (SNL) and OYO Geospace. The system requires data from multiple three-component wall-locking accelerometer packages and telemeters digital data to the surface in real-time. Due to the multiplicity of measurement stations and the real-time data link, acquisition time for the borehole seismic survey is significantly reduced. The system was tested at the Chevron La Habra Test Site using Chevron`s clamped axial borehole vibrator as the seismic source. Several source and receiver fans were acquired using a four-station version of the advanced system. For comparison purposes, an equivalent data set was acquired using a standard analog wall-locking geophone receiver. The test data indicate several enhancements provided by the multi-station receiver relative to the standard, drastically improved signal-to-noise ratio, increased signal bandwidth, the detection of multiple reflectors, and a true 4:1 reduction in survey time.
Traction force microscopy on soft elastic substrates: A guide to recent computational advances.
Schwarz, Ulrich S; Soiné, Jérôme R D
2015-11-01
The measurement of cellular traction forces on soft elastic substrates has become a standard tool for many labs working on mechanobiology. Here we review the basic principles and different variants of this approach. In general, the extraction of the substrate displacement field from image data and the reconstruction procedure for the forces are closely linked to each other and limited by the presence of experimental noise. We discuss different strategies to reconstruct cellular forces as they follow from the foundations of elasticity theory, including two- versus three-dimensional, inverse versus direct and linear versus non-linear approaches. We also discuss how biophysical models can improve force reconstruction and comment on practical issues like substrate preparation, image processing and the availability of software for traction force microscopy. This article is part of a Special Issue entitled: Mechanobiology. PMID:26026889
A remark on the computation of shear-horizontal and torsional modes in elastic waveguides.
Gravenkamp, Hauke
2016-07-01
When modeling the propagation of elastic guided waves in plates or cylinders, Finite Element based numerical methods such as the Scaled Boundary Finite Element Method (SBFEM) or the Semi-Analytical Finite Element (SAFE) Method lead to an eigenvalue problem to be solved at each frequency. For the particular case of shear horizontal modes in a homogeneous plate or torsional modes in a homogeneous cylinder, the problem can be drastically simplified. The eigenvalues become simple functions of the frequency, while the eigenvectors are constant. The current contribution discusses how this behavior is represented in the numerical formulation and derives the expressions for the eigenvalues and eigenvectors as well as the dynamic stiffness matrix of infinite elastic waveguides. PMID:27014855
A fourth order accurate finite difference scheme for the computation of elastic waves
NASA Technical Reports Server (NTRS)
Bayliss, A.; Jordan, K. E.; Lemesurier, B. J.; Turkel, E.
1986-01-01
A finite difference for elastic waves is introduced. The model is based on the first order system of equations for the velocities and stresses. The differencing is fourth order accurate on the spatial derivatives and second order accurate in time. The model is tested on a series of examples including the Lamb problem, scattering from plane interf aces and scattering from a fluid-elastic interface. The scheme is shown to be effective for these problems. The accuracy and stability is insensitive to the Poisson ratio. For the class of problems considered here it is found that the fourth order scheme requires for two-thirds to one-half the resolution of a typical second order scheme to give comparable accuracy.
Analytical and Computational Study of One-Dimensional Impact of Graded Elastic Solids
NASA Astrophysics Data System (ADS)
Scheidler, Mike; Gazonas, George
2002-07-01
Some of the recent efforts to improve the ballistic performance of lightweight armors utilize functionally graded materials to provide a continuous transition in properties between dissimilar materials. It has been conjectured that the elimination of abrupt acoustic impedance changes may result in beneficial stress wave attenuation. We examine this issue for some idealized one-dimensional problems in which all materials are linear elastic. Exact solutions are compared with DYNA3D simulations of the same problems.
A study of self-propelled elastic cylindrical micro-swimmers using modeling and computation
NASA Astrophysics Data System (ADS)
Shi, Lingling; Čanić, Sunčica; Quaini, Annalisa; Pan, Tsorng-Whay
2016-06-01
We study propulsion of micro-swimmers in 3D creeping flow. The swimmers are assumed to be made of elastic cylindrical hollow tubes. The swimming is generated by the contractions of the tube's elastic membrane walls producing a traveling wave in the form of a "step-function" traversing the swimmer from right to left, propelling the swimmer from left to right. The problem is motivated by medical applications such as drug delivery. The influence of several non-dimensional design parameters on the velocity of the swimmer is investigated, including the swimmer aspect ratio, and the amplitude of the traveling wave relative to the swimmer radius. An immersed boundary method based on a finite element method approach is successfully combined with an elastic spring network model to simulate the two-way fluid-structure interaction coupling between the elastic cylindrical tube and the flow of a 3D viscous, incompressible fluid. To gain a deeper insight into the influence of various parameters on the swimmer speed, a reduced 1D fluid-structure interaction model was derived and validated. It was found that fast swimmers are those with large tube aspect ratios, and with the amplitude of the traveling wave which is roughly 50% of the reference swimmer radius. It was shown that the speed of our "optimal swimmer" is around 1.5 swimmer lengths per second, which is at the top of the class of all currently manufactured micro-swimmers swimming in low Reynolds number flows (Re =10-6), reported in [11].
Uniform asymptotic formula for the Aharonov Bohm wavefield
NASA Astrophysics Data System (ADS)
Hannay, J. H.
2016-06-01
A uniform asymptotic formula for the Aharonov–Bohm wavefield (that of a plane quantum wave scattered by a thin straight solenoid) far away from the solenoid is obtained in a direct way. Actually quite good accuracy is achieved even down to one wavelength away. The error is numerically of order radius^(‑3/2) for all values of polar angle, including directly forwards. Several previous formulas, uniform and otherwise, for the far field limit exist in the literature. All contain the essential ingredient: the Fresnel integral (complex error function), but ordinarily the error in these formulas is of order radius^(‑1/2) in the forwards direction where the Fresnel contribution is most important.
Multicomponent wavefield characterization with a novel scanning laser interferometer.
Blum, Thomas E; van Wijk, Kasper; Pouet, Bruno; Wartelle, Alexis
2010-07-01
The in-plane component of the wavefield provides valuable information about media properties from seismology to nondestructive testing. A new compact scanning laser ultrasonic interferometer collects light scattered away from the angle of incidence to provide the absolute ultrasonic displacement for both the out-of-plane and an in-plane components. This new system is tested by measuring the radial and vertical polarization of a Rayleigh wave in an aluminum half-space. The estimated amplitude ratio of the horizontal and vertical displacement agrees well with the theoretical value. The phase difference exhibits a small bias between the two components due to a slightly different frequency response between the two processing channels of the prototype electronic circuitry. PMID:20687699
Multicomponent wavefield characterization with a novel scanning laser interferometer
Blum, Thomas E.; Wijk, Kasper van; Pouet, Bruno; Wartelle, Alexis
2010-07-15
The in-plane component of the wavefield provides valuable information about media properties from seismology to nondestructive testing. A new compact scanning laser ultrasonic interferometer collects light scattered away from the angle of incidence to provide the absolute ultrasonic displacement for both the out-of-plane and an in-plane components. This new system is tested by measuring the radial and vertical polarization of a Rayleigh wave in an aluminum half-space. The estimated amplitude ratio of the horizontal and vertical displacement agrees well with the theoretical value. The phase difference exhibits a small bias between the two components due to a slightly different frequency response between the two processing channels of the prototype electronic circuitry.
Elasticity and geometry: a computational model of the Heineke-Mikulicz strictureplasty.
Tsamis, Alkiviadis; Pocivavsek, Luka; Vorp, David A
2014-11-01
Crohn's disease is a challenging inflammatory process with a propensity for focal gastro-intestinal tract inflammation and stricture. Surgically, Crohn's is often treated with resection. However, a subtype of diffuse disease with multiple strictures is treated by strictureplasty procedures in hope of avoiding short-gut syndrome. Prior work by Pocivavsek et al. defined the geometry of a Heineke-Mikulicz strictureplasty. Here, we bring this analysis one step closer to clinical and biological relevance by calculating the mechanical stresses and strains that the strictureplasty deformation generates on a model intestinal wall. The small bowel is simulated as a linearly elastic isotropic deformable cylindrical shell using finite element modeling. Data show a divergence in elastic response between the anti-mesenteric and mesenteric halves. The anti-mesenteric surface shows a bending dominated elastic response that correlates with the prior purely geometric analysis. However, the mesenteric half is not a neutral bystander during strictureplasty formation, as geometric arguments predict. Strong in-plane stretching strains develop in a rim around the image of the transverse closure, which may impact local perfusion and serve as sites of disease recurrence. Lastly, nearly all the deformation energy is stored in the central vertex stitch, placing this part at highest risk of dehiscence. This study enhances our understanding of mechanical response in complex nonlinear cylindrical geometries like the surgically manipulated intestinal tract. The developed framework serves as a platform for future addition of more complex clinically relevant parameters to our model, including real tissue properties, anisotropy, blood supply modeling, and patient deriver anatomic factors. PMID:24671519
Analytical and Computational Study of One-Dimensional Impact of Graded Elastic Materials
NASA Astrophysics Data System (ADS)
Scheidler, Mike; Gazonas, George
2001-06-01
Traditional armor designs consist of discrete layers of homogeneous materials with different properties. Some recent efforts to improve the ballistic performance of lightweight armors involve replacing several or all of these layers by a monolithic material with continuous property variations. The stress waves generated by ballistic impact in these functionally graded materials tend to be much more complex than in the corresponding layered materials. This fact has led us to examine whether existing hydrocodes can capture the complex wave structures in graded materials and, in particular, whether these codes can accurately predict the peak stresses that may lead to material failure. In this paper we report on the initial phase of this program. We consider the one-dimensional impact of a homogeneous projectile on a graded target, with or without a homogeneous backing plate. All materials are assumed to be linear elastic. Smooth gradings in the density and the elastic modulus are chosen for which exact analytical solutions (including multiple reflections) can be derived using Laplace transforms. These exact solutions are compared with DYNA3D simulations of the same problems. The effects artificial viscosity, nodal spacing (i.e., mesh grading), and number of elements on the DYNA3D solutions are examined.
3D elastic full waveform inversion: case study from a land seismic survey
NASA Astrophysics Data System (ADS)
Kormann, Jean; Marti, David; Rodriguez, Juan-Esteban; Marzan, Ignacio; Ferrer, Miguel; Gutierrez, Natalia; Farres, Albert; Hanzich, Mauricio; de la Puente, Josep; Carbonell, Ramon
2016-04-01
Full Waveform Inversion (FWI) is one of the most advanced processing methods that is recently reaching a mature state after years of solving theoretical and technical issues such as the non-uniqueness of the solution and harnessing the huge computational power required by realistic scenarios. BSIT (Barcelona Subsurface Imaging Tools, www.bsc.es/bsit) includes a FWI algorithm that can tackle with very complex problems involving large datasets. We present here the application of this system to a 3D dataset acquired to constrain the shallow subsurface. This is where the wavefield is the most complicated, because most of the wavefield conversions takes place in the shallow region and also because the media is much more laterally heterogeneous. With this in mind, at least isotropic elastic approximation would be suitable as kernel engine for FWI. The current study explores the possibilities to apply elastic isotropic FWI using only the vertical component of the recorded seismograms. The survey covers an area of 500×500 m2, and consists in a receivers grid of 10 m×20 m combined with a 250 kg accelerated weight-drop as source on a displaced grid of 20 m×20 m. One of the main challenges in this case study is the costly 3D modeling that includes topography and substantial free surface effects. FWI is applied to a data subset (shooting lines 4 to 12), and is performed for 3 frequencies ranging from 15 to 25 Hz. The starting models are obtained from travel-time tomography and the all computation is run on 75 nodes of Mare Nostrum supercomputer during 3 days. The resulting models provide a higher resolution of the subsurface structures, and show a good correlation with the available borehole measurements. FWI allows to extend in a reliable way this 1D knowledge (borehole) to 3D.
Quaini, A.; Canic, S.; Glowinski, R.; Igo, S.; Hartley, C.J.; Zoghbi, W.; Little, S.
2011-01-01
This work presents a validation of a fluid-structure interaction computational model simulating the flow conditions in an in vitro mock heart chamber modeling mitral valve regurgitation during the ejection phase during which the trans-valvular pressure drop and valve displacement are not as large. The mock heart chamber was developed to study the use of 2D and 3D color Doppler techniques in imaging the clinically relevant complex intra-cardiac flow events associated with mitral regurgitation. Computational models are expected to play an important role in supporting, refining, and reinforcing the emerging 3D echocardiographic applications. We have developed a 3D computational fluid-structure interaction algorithm based on a semi-implicit, monolithic method, combined with an arbitrary Lagrangian-Eulerian approach to capture the fluid domain motion. The mock regurgitant mitral valve corresponding to an elastic plate with a geometric orifice, was modeled using 3D elasticity, while the blood flow was modeled using the 3D Navier-Stokes equations for an incompressible, viscous fluid. The two are coupled via the kinematic and dynamic conditions describing the two-way coupling. The pressure, the flow rate, and orifice plate displacement were measured and compared with numerical simulation results. In-line flow meter was used to measure the flow, pressure transducers were used to measure the pressure, and a Doppler method developed by one of the authors was used to measure the axial displacement of the orifice plate. The maximum recorded difference between experiment and numerical simulation for the flow rate was 4%, the pressure 3.6%, and for the orifice displacement 15%, showing excellent agreement between the two. PMID:22138194
Quaini, A; Canic, S; Glowinski, R; Igo, S; Hartley, C J; Zoghbi, W; Little, S
2012-01-10
This work presents a validation of a fluid-structure interaction computational model simulating the flow conditions in an in vitro mock heart chamber modeling mitral valve regurgitation during the ejection phase during which the trans-valvular pressure drop and valve displacement are not as large. The mock heart chamber was developed to study the use of 2D and 3D color Doppler techniques in imaging the clinically relevant complex intra-cardiac flow events associated with mitral regurgitation. Computational models are expected to play an important role in supporting, refining, and reinforcing the emerging 3D echocardiographic applications. We have developed a 3D computational fluid-structure interaction algorithm based on a semi-implicit, monolithic method, combined with an arbitrary Lagrangian-Eulerian approach to capture the fluid domain motion. The mock regurgitant mitral valve corresponding to an elastic plate with a geometric orifice, was modeled using 3D elasticity, while the blood flow was modeled using the 3D Navier-Stokes equations for an incompressible, viscous fluid. The two are coupled via the kinematic and dynamic conditions describing the two-way coupling. The pressure, the flow rate, and orifice plate displacement were measured and compared with numerical simulation results. In-line flow meter was used to measure the flow, pressure transducers were used to measure the pressure, and a Doppler method developed by one of the authors was used to measure the axial displacement of the orifice plate. The maximum recorded difference between experiment and numerical simulation for the flow rate was 4%, the pressure 3.6%, and for the orifice displacement 15%, showing excellent agreement between the two. PMID:22138194
Scattered Wavefield Within the San Andreas Fault System, California
NASA Astrophysics Data System (ADS)
Taira, T.; Silver, P. G.; Niu, F.; Nadeau, R. M.
2004-12-01
Since transient aseismic deformation at seismogenic depth is one of the key phenomena related to earthquake occurrence, it is important to estimate the physical characteristics of such stress/strain transients within the deeper structure of the fault zone. Analysis of coda (scattered) waves has the potential for identifying such transients because the scattered wavefield is attributed mainly to small-scale heterogeneity that is likely formed by these transient events. In addition, the sampling area of scattered waves is concentrated within the fault zone, compared to the area sampled by direct waves. We have begun a program to map the spatial distribution of fault-zone scatterers and their time dependence within two regions of San Andreas Fault system: the Hayward Fault and the Parkfield segment of the San Andreas Fault. In order to most reliably evaluate the scattered wavefields, we limited our analysis to records from borehole seismographs recorded by the Hayward Fault Network (HFN) and the High-Resolution Seismic Network (HRSN) in each area. For the Hayward fault, we mapped the spatial distribution of scatterers by analysis of the S-wave coda amplification factor (CAF) in a manner similar to Taira and Yomogida (2003). CAF is defined as the amplitude ratio of coda waves among different stations after corrections for source, station, and overall propagation effects (e.g., coda Q). This parameter allows for a statistical characterization of the distribution of scatterers. The station effect for each station and the coda Q averaged over all the seismograms in this area were estimated by the coda-normalization and maximum likelihood methods, respectively, using five regional earthquakes (epicentral distance > 50 km). We evaluated the CAF value of each source-station pair for the transverse-component, using 294 seismograms for 39 local earthquakes (epicentral distance < 50 km) recorded by 14 stations of the HFN. A map of CAF values for stations near the Hayward Fault
IRIS Arrays: Observing Wavefields at Multiple Scales and Frequencies
NASA Astrophysics Data System (ADS)
Sumy, D. F.; Woodward, R.; Frassetto, A.
2014-12-01
The Incorporated Research Institutions for Seismology (IRIS) provides instruments for creating and operating seismic arrays at a wide range of scales. As an example, for over thirty years the IRIS PASSCAL program has provided instruments to individual Principal Investigators to deploy arrays of all shapes and sizes on every continent. These arrays have ranged from just a few sensors to hundreds or even thousands of sensors, covering areas with dimensions of meters to thousands of kilometers. IRIS also operates arrays directly, such as the USArray Transportable Array (TA) as part of the EarthScope program. Since 2004, the TA has rolled across North America, at any given time spanning a swath of approximately 800 km by 2,500 km, and thus far sampling 2% of the Earth's surface. This achievement includes all of the lower-48 U.S., southernmost Canada, and now parts of Alaska. IRIS has also facilitated specialized arrays in polar environments and on the seafloor. In all cases, the data from these arrays are freely available to the scientific community. As the community of scientists who use IRIS facilities and data look to the future they have identified a clear need for new array capabilities. In particular, as part of its Wavefields Initiative, IRIS is exploring new technologies that can enable large, dense array deployments to record unaliased wavefields at a wide range of frequencies. Large-scale arrays might utilize multiple sensor technologies to best achieve observing objectives and optimize equipment and logistical costs. Improvements in packaging and power systems can provide equipment with reduced size, weight, and power that will reduce logistical constraints for large experiments, and can make a critical difference for deployments in harsh environments or other situations where rapid deployment is required. We will review the range of existing IRIS array capabilities with an overview of previous and current deployments and examples of data and results. We
EXPLORING FOR SUBTLE MISSION CANYON STRATIGRAPHIC TRAPS WITH ELASTIC WAVEFIELD SEISMIC TECHNOLOGY
John Beecherl
2003-10-01
A source-receiver geometry was designed for a 9C3D seismic survey in Montrail County, North Dakota, that will involve the largest number of active 3-component stations (1,800 to 2,100) ever attempted in an onshore U.S. multicomponent seismic survey. To achieve the data-acquisition objectives, 3-component geophone strings will be provided by the Bureau of Economic Geology, Dawson Geophysical, and Vecta Technology. Data acquisition will commence in late October 2003. The general objective of this study is to demonstrate the value of multicomponent seismic technology for exploring for subtle oolitic-bank reservoirs in the Mission Canyon Formation of the Williston Basin. The work tasks done during this report period concentrated on developing an optimal design for the seismic survey. This first semiannual report defines the geographical location and geometrical shape of the survey and documents the key acquisition parameters that will be implemented to yield high-fold, high-resolution 9-component seismic data.
NASA Astrophysics Data System (ADS)
Gilabert, F. A.; Roux, J.-N.; Castellanos, A.
2008-09-01
The quasistatic behavior of a simple two-dimensional model of a cohesive powder under isotropic loads is investigated by discrete element simulations. We ignore contact plasticity and focus on the effect of geometry and collective rearrangements on the material behavior. The loose packing states, as assembled and characterized in a previous numerical study [Gilabert, Roux, and Castellanos, Phys. Rev. E 75, 011303 (2007)], are observed, under growing confining pressure P , to undergo important structural changes, while solid fraction Φ irreversibly increases (typically, from 0.4-0.5 to 0.75-0.8). The system state goes through three stages, with different forms of the plastic consolidation curve, i.e., Φ as a function of the growing reduced pressure P*=Pa/F0 , defined with adhesion force F0 and grain diameter a . In the low-confinement regime (I), the system undergoes negligible plastic compaction, and its structure is influenced by the assembling process. In regime II the material state is independent of initial conditions, and the void ratio varies linearly with lnP [i.e., Δ(1/Φ)=λΔ(lnP*) ], as described in the engineering literature. Plasticity index λ is reduced in the presence of a small rolling resistance (RR). In the last stage of compaction (III), Φ approaches an asymptotic, maximum solid fraction Φmax , as a power law Φmax-Φ∝(P*)-α , with α≃1 , and properties of cohesionless granular packs are gradually retrieved. Under consolidation, while the range ξ of fractal density correlations decreases, force patterns reorganize from self-balanced clusters to force chains, with correlative evolutions of force distributions, and elastic moduli increase by a large amount. Plastic deformation events correspond to very small changes in the network topology, while the denser regions tend to move like rigid bodies. Elastic properties are dominated by the bending of thin junctions in loose systems. For growing RR those tend to form particle chains, the
NASA Astrophysics Data System (ADS)
French, S. W.; Lekic, V.; Romanowicz, B. A.
2010-12-01
As global waveform-modeling schemes rooted in perturbation theory are supplanted by fully numerical alternatives, such as the Spectral Element Method (e.g. SEM: Komatitsch and Tromp, 2002), the improved wavefield accuracy for complex 3D structures also carries increased computational cost. Lekic and Romanowicz (2010) inverted waveforms of fundamental and higher mode surface waves for a radially anisotropic upper-mantle Vs model using SEM (SEMum). The SEM computations were made feasible by an appropriate choice of cutoff period (T≥ 60 s.), as well as the implementation of a homogenized anisotropic crustal layer based on fitting of short period group velocity dispersion curves. These choices allowed for an efficient SEM mesh undeformed by true Moho topography. Further, instead of homogenization of a possibly biased a priori crustal model, Lekic and Romanowicz jointly inverted for the crustal layer, constrained by surface wave group velocity dispersion maps for T≥ 25 s. We are currently developing a radially anisotropic Vs model of the whole mantle using SEM, following an approach broadly similar to that employed in SEMum. Extension of this methodology to imaging of lower-mantle structure requires the inclusion of a body wave dataset, and thus shorter-period modeling of the global wavefield (T≥ 32 s.). While this period range dictates finer sampling of our SEM mesh, reduced computational cost is still possible through the crustal homogenization scheme. Here, we first discuss the development of an analogous homogenized crustal model and its validity for both the fundamental and higher mode surface wave and the body wave datasets. We focus on maintaining a simplified Moho topography, thus obviating expensive deformation of the SEM mesh, while accurately treating valuable surface-reflected body wave phases (ex: multiple ScS). Second, we discuss implications of treating the crust in this manner for the overall inversion methodology. In particular, we intend to
Plietsch, R; Bourauel, C; Drescher, D; Nellen, B
1994-04-01
Metals are the most commonly used materials in the construction of orthodontic appliances designed for the correction of malocclusions. Knowledge of the force systems at work is a prerequisite for judging the functionality of these appliances. The elasticity parameters (Young's E-moduli, strain limits) of the alloys employed can be drawn upon to calculate numerically forces and torsional moments. Both tensile tests and bending experiments are used to determine the E-moduli and strain limits of standard steel and highly flexible NiTi wires frequently used in orthodontics. However, parameters obtained by tensile tests are less suited for studying the mechanical properties of orthodontic appliances. Since bending deformation prevails, bending experiments should be preferred method for ascertaining the relevant parameters. This study, therefore, presents a new experimental method for testing the bend ability of highly flexible materials and the determination of the underlying material parameters. A comparison of calculated force systems with direct measurements revealed that bending parameters lead to an appropriate description of forces and moments generated during clinical treatment, whereas calculations based on tensile test parameters differ substantially. The bending test proposed here is, thus, a suitable means for dependably predicting the force systems produced by an orthodontic appliance and the test therefore can contribute to an accurate design of new types of therapeutic devices. PMID:8194813
NASA Astrophysics Data System (ADS)
Bagnasco, S.; Berzano, D.; Guarise, A.; Lusso, S.; Masera, M.; Vallero, S.
2015-12-01
The INFN computing centre in Torino hosts a private Cloud, which is managed with the OpenNebula cloud controller. The infrastructure offers Infrastructure-as-a-Service (IaaS) and Platform-as-a-Service (PaaS) services to different scientific computing applications. The main stakeholders of the facility are a grid Tier-2 site for the ALICE collaboration at LHC, an interactive analysis facility for the same experiment and a grid Tier-2 site for the BESIII collaboration, plus an increasing number of other small tenants. The dynamic allocation of resources to tenants is partially automated. This feature requires detailed monitoring and accounting of the resource usage. We set up a monitoring framework to inspect the site activities both in terms of IaaS and applications running on the hosted virtual instances. For this purpose we used the ElasticSearch, Logstash and Kibana (ELK) stack. The infrastructure relies on a MySQL database back-end for data preservation and to ensure flexibility to choose a different monitoring solution if needed. The heterogeneous accounting information is transferred from the database to the ElasticSearch engine via a custom Logstash plugin. Each use-case is indexed separately in ElasticSearch and we setup a set of Kibana dashboards with pre-defined queries in order to monitor the relevant information in each case. For the IaaS metering, we developed sensors for the OpenNebula API. The IaaS level information gathered through the API is sent to the MySQL database through an ad-hoc developed RESTful web service. Moreover, we have developed a billing system for our private Cloud, which relies on the RabbitMQ message queue for asynchronous communication to the database and on the ELK stack for its graphical interface. The Italian Grid accounting framework is also migrating to a similar set-up. Concerning the application level, we used the Root plugin TProofMonSenderSQL to collect accounting data from the interactive analysis facility. The BESIII
Finite-Difference Algorithm for Simulating 3D Electromagnetic Wavefields in Conductive Media
NASA Astrophysics Data System (ADS)
Aldridge, D. F.; Bartel, L. C.; Knox, H. A.
2013-12-01
Electromagnetic (EM) wavefields are routinely used in geophysical exploration for detection and characterization of subsurface geological formations of economic interest. Recorded EM signals depend strongly on the current conductivity of geologic media. Hence, they are particularly useful for inferring fluid content of saturated porous bodies. In order to enhance understanding of field-recorded data, we are developing a numerical algorithm for simulating three-dimensional (3D) EM wave propagation and diffusion in heterogeneous conductive materials. Maxwell's equations are combined with isotropic constitutive relations to obtain a set of six, coupled, first-order partial differential equations governing the electric and magnetic vectors. An advantage of this system is that it does not contain spatial derivatives of the three medium parameters electric permittivity, magnetic permeability, and current conductivity. Numerical solution methodology consists of explicit, time-domain finite-differencing on a 3D staggered rectangular grid. Temporal and spatial FD operators have order 2 and N, where N is user-selectable. We use an artificially-large electric permittivity to maximize the FD timestep, and thus reduce execution time. For the low frequencies typically used in geophysical exploration, accuracy is not unduly compromised. Grid boundary reflections are mitigated via convolutional perfectly matched layers (C-PMLs) imposed at the six grid flanks. A shared-memory-parallel code implementation via OpenMP directives enables rapid algorithm execution on a multi-thread computational platform. Good agreement is obtained in comparisons of numerically-generated data with reference solutions. EM wavefields are sourced via point current density and magnetic dipole vectors. Spatially-extended inductive sources (current carrying wire loops) are under development. We are particularly interested in accurate representation of high-conductivity sub-grid-scale features that are common
NASA Astrophysics Data System (ADS)
Galiana-Merino, J. J.; Rosa-Cintas, S.; Rosa-Herranz, J. L.; Molina-Palacios, S.; Martinez-Espla, J. J.
2011-12-01
Microzonation studies using ambient noise measurements constitute an extended and useful procedure for determine the local soil characteristics and its response due to an earthquake. Several methods exist for analyzing the noise measurements, being the most popular the horizontal-to-vertical spectral ratio (H/V) and the array techniques, i.e. the frequency-wavenumber (F-K) transform. Many works exist about this topic and it stills being an ongoing debate about ambient noise composition, whether body or surface waves constitute most of it, showing the importance of identifying the different kinds of waves presents in a seismic record. In this work we utilize a new method of time-frequency polarization analysis, based on the stationary wavelet packet transform, to investigate how the polarization characteristics of the wavefield influence in the application of ambient noise techniques. The signals are divided in different bands, according to their reciprocal ellipticity values and then the H/V method and the F-K array analysis are computed for each band. The qualitative and quantitative comparison between the original curve and the obtained for the analyzed intervals provide information about the signals composition, showing that the major components of the seismic noise present reciprocal ellipticity values lower than 0.5. The efficient application of the studied techniques by using just the main a part of the entire signal, [0 - 0.5], is also evaluated, showing favorable results.
NASA Astrophysics Data System (ADS)
Sinescu, Cosmin; Negrutiu, Meda; Topala, Florin; Ionita, Ciprian; Negru, Radu; Fabriky, Mihai; Marcauteanu, Corina; Bradu, Adrian; Dobre, George; Marsavina, Liviu; Rominu, Mihai; Podoleanu, Adrian
2011-10-01
Dental onlays are restorations used to repair rear teeth that have a mild to moderate amount of decay. They can also be used to restore teeth that are cracked or fractured if the damage is not severe enough to require a dental crown. The use of onlays requires less tooth reduction than does the use of metal fillings. This allows dentists to conserve more of a patient's natural tooth structure in the treatment process. The aims of this study are to evaluate the biomechanical comportment of the dental onlays, by using the 3D photo elasticity method and to investigate the integrity of the structures and their fitting to the dental support. For this optical coherence tomography and micro-computed tomography were employed. Both methods were used to investigate 37 dental onlays, 17 integral polymeric and 20 integral ceramic. The results permit to observe materials defects inside the ceramic or polymeric onlays situate in the biomechanically tensioned areas that could lead to fracture of the prosthetic structure. Marginal fitting problems of the onlays related to the teeth preparations were presented in order to observe the possibility of secondary cavities. The resulted images from the optical coherence tomography were verified by the micro-computed tomography. In conclusion, the optical coherence tomography can be used as a clinical method in order to evaluate the integrity of the dental ceramic and polymeric onlays and to investigate the quality of the marginal fitting to the teeth preparations.
VFLOW2D - A Vorte-Based Code for Computing Flow Over Elastically Supported Tubes and Tube Arrays
WOLFE,WALTER P.; STRICKLAND,JAMES H.; HOMICZ,GREGORY F.; GOSSLER,ALBERT A.
2000-10-11
A numerical flow model is developed to simulate two-dimensional fluid flow past immersed, elastically supported tube arrays. This work is motivated by the objective of predicting forces and motion associated with both deep-water drilling and production risers in the oil industry. This work has other engineering applications including simulation of flow past tubular heat exchangers or submarine-towed sensor arrays and the flow about parachute ribbons. In the present work, a vortex method is used for solving the unsteady flow field. This method demonstrates inherent advantages over more conventional grid-based computational fluid dynamics. The vortex method is non-iterative, does not require artificial viscosity for stability, displays minimal numerical diffusion, can easily treat moving boundaries, and allows a greatly reduced computational domain since vorticity occupies only a small fraction of the fluid volume. A gridless approach is used in the flow sufficiently distant from surfaces. A Lagrangian remap scheme is used near surfaces to calculate diffusion and convection of vorticity. A fast multipole technique is utilized for efficient calculation of velocity from the vorticity field. The ability of the method to correctly predict lift and drag forces on simple stationary geometries over a broad range of Reynolds numbers is presented.
Acquisition and analysis of angle-beam wavefield data
Dawson, Alexander J.; Michaels, Jennifer E.; Levine, Ross M.; Chen, Xin; Michaels, Thomas E.
2014-02-18
Angle-beam ultrasonic testing is a common practical technique used for nondestructive evaluation to detect, locate, and characterize a variety of material defects and damage. Greater understanding of the both the incident wavefield produced by an angle-beam transducer and the subsequent scattering from a variety of defects and geometrical features is anticipated to increase the reliability of data interpretation. The focus of this paper is on acquiring and analyzing propagating waves from angle-beam transducers in simple, defect-free plates as a first step in the development of methods for flaw characterization. Unlike guided waves, which excite the plate throughout its thickness, angle-beam bulk waves bounce back and forth between the plate surfaces, resulting in the well-known multiple “skips” or “V-paths.” The experimental setup consists of a laser vibrometer mounted on an XYZ scanning stage, which is programmed to move point-to-point on a rectilinear grid to acquire waveform data. Although laser vibrometry is now routinely used to record guided waves for which the frequency content is below 1 MHz, it is more challenging to acquire higher frequency bulk waves in the 1–10 MHz range. Signals are recorded on the surface of an aluminum plate that were generated from a 5 MHz, 65° refracted angle, shear wave transducer-wedge combination. Data are analyzed directly in the x-t domain, via a slant stack Radon transform in the τ-p (offset time-slowness) domain, and via a 2-D Fourier transform in the ω-k domain, thereby enabling identification of specific arrivals and modes. Results compare well to those expected from a simple ray tracing analysis except for the unexpected presence of a strong Rayleigh wave.
Exploring the Seismic Wavefield with EarthScope's Transportable Array
NASA Astrophysics Data System (ADS)
Ammon, C. J.; Lay, T.
2007-12-01
The unprecedented aperture and dense station spacing of EarthScope's Transportable Array (TA) provide opportunities to explore subtle features of the seismic wave field using array methods, including stacking to enhance signals from small, distant events, and using spatial animations and slowness spectra to investigate the origin and nature of long-period (tens to hundreds of seconds) seismic coda. We illustrate the TA's capabilities using teleseismic short-period P waves and complete broadband ground motion recordings. The high quality and uniform installation procedures of the TA produce many stations with superb short-period signal- to-noise ratios enabling direct detection of tiny amplitude (~2nm) P-waves from the mb 4.2, 09 October, 2006 North Korean nuclear test. Stacking of many TA signals resolves subtle complexity in the source signals. Guided by seismic wavefield animations from large earthquakes, which display many wavefronts sweeping across the array, visual inspection and seismic slowness spectra of the extended coda resolve off-great circle Rayleigh-wave arrivals. The scattered waves approach the array from preferred directions depending on the source location relative to the TA. For example, scattered signals from the 01 April, 2007 Solomon Islands earthquakes are preferentially located to the west and northwest, with timing that suggests scattering from subduction zones along the margin of the western and north Pacific ocean. For the 16 August, 2007 Central Peru earthquake, Rayleigh wave coda contains significant energy arriving from the north and northeast, scattered near the North American margin. Analysis of scattered arrivals may lead to improved constraints on lateral gradients in lithospheric structures.
On the reliability and limitations of the SPAC method with a directional wavefield
NASA Astrophysics Data System (ADS)
Luo, Song; Luo, Yinhe; Zhu, Lupei; Xu, Yixian
2016-03-01
The spatial autocorrelation (SPAC) method is one of the most efficient ways to extract phase velocities of surface waves from ambient seismic noise. Most studies apply the method based on the assumption that the wavefield of ambient noise is diffuse. However, the actual distribution of sources is neither diffuse nor stationary. In this study, we examined the reliability and limitations of the SPAC method with a directional wavefield. We calculated the SPAC coefficients and phase velocities from a directional wavefield for a four-layer model and characterized the limitations of the SPAC. We then applied the SPAC method to real data in Karamay, China. Our results show that, 1) the SPAC method can accurately measure surface wave phase velocities from a square array with a directional wavefield down to a wavelength of twice the shortest interstation distance; and 2) phase velocities obtained from real data by the SPAC method are stable and reliable, which demonstrates that this method can be applied to measure phase velocities in a square array with a directional wavefield.
Anderson, C.A.
1981-01-01
Six years ago the Reactor Safety Research Division of the Nuclear Regulatory Commission (NRC) approached the Los Alamos National Laboratory to develop a comprehensive concrete structural analysis code to predict the static and dynamic behavior of Prestressed Concrete Reactor Vessels (PCRVs) that serve as the containment structure of a High-Temperature Gas-Cooled Reactor. The PCRV is a complex concrete structure that must be modeled in three dimensions and posseses other complicating features such as a steel liner for the reactor cavity and woven cables embedded vertically in the PCRV and wound circumferentially on the outside of the PCRV. The cables, or tendons, are used for prestressing the reactor vessel. In addition to developing the computational capability to predict inelastic three dimensional concrete structural behavior, the code response was verified against documented experiments on concrete structural behavior. This code development/verification effort is described.
NASA Astrophysics Data System (ADS)
Borisov, Dmitry; Singh, Satish C.; Fuji, Nobuaki
2015-09-01
Seismic full waveform inversion is an objective method to estimate elastic properties of the subsurface and is an important area of research, particularly in seismic exploration community. It is a data-fitting approach, where the difference between observed and synthetic data is minimized iteratively. Due to a very high computational cost, the practical implementation of waveform inversion has so far been restricted to a 2-D geometry with different levels of physics incorporated in it (e.g. elasticity/viscoelasticity) or to a 3-D geometry but using an acoustic approximation. However, the earth is three-dimensional, elastic and heterogeneous and therefore a full 3-D elastic inversion is required in order to obtain more accurate and valuable models of the subsurface. Despite the recent increase in computing power, the application of 3-D elastic full waveform inversion to real-scale problems remains quite challenging on the current computer architecture. Here, we present an efficient method to perform 3-D elastic full waveform inversion for time-lapse seismic data using a finite-difference injection method. In this method, the wavefield is computed in the whole model and is stored on a surface above a finite volume where the model is perturbed and localized inversion is performed. Comparison of the final results using the 3-D finite-difference injection method and conventional 3-D inversion performed within the whole volume shows that our new method provides significant reductions in computational time and memory requirements without any notable loss in accuracy. Our approach shows a big potential for efficient reservoir monitoring in real time-lapse experiments.
Korsa, Radim; Lukes, Jaroslav; Sepitka, Josef; Mares, Tomas
2015-08-01
Knowledge of the anisotropic elastic properties of osteon and osteonal lamellae provides a better understanding of various pathophysiological conditions, such as aging, osteoporosis, osteoarthritis, and other degenerative diseases. For this reason, it is important to investigate and understand the elasticity of cortical bone. We created a bidirectional micromechanical model based on inverse homogenization for predicting the elastic properties of osteon and osteonal lamellae of cortical bone. The shape, the dimensions, and the curvature of osteon and osteonal lamellae are described by appropriately chosen curvilinear coordinate systems, so that the model operates close to the real morphology of these bone components. The model was used to calculate nine orthotropic elastic constants of osteonal lamellae. The input values have the elastic properties of a single osteon. We also expressed the dependence of the elastic properties of the lamellae on the angle of orientation. To validate the model, we performed nanoindentation tests on several osteonal lamellae. We compared the experimental results with the calculated results, and there was good agreement between them. The inverted model was used to calculate the elastic properties of a single osteon, where the input values are the elastic constants of osteonal lamellae. These calculations reveal that the model can be used in both directions of homogenization, i.e., direct homogenization and also inverse homogenization. The model described here can provide either the unknown elastic properties of a single lamella from the known elastic properties at the level of a single osteon, or the unknown elastic properties of a single osteon from the known elastic properties at the level of a single lamella. PMID:25901781
NASA Astrophysics Data System (ADS)
Zhou, X.; Sha, D.; Tamma, K. K.
2006-10-01
Of interest here are the class of static/dynamic finite deformation problems that arise in computational mechanics, and the question of the suitability in employing the total strain measure for this class of problems is raised. An attempt to resolve the problem by proposing a new arbitrary reference configuration (ARC) framework is described in this exposition. The ARC framework consists of the ARC elasticity, which bridges the Truesdell stress rate hypo-elasticity and the St. Venant-Kirchhoff hyperelasticity, and the ARC Lagrangian formulation, which bridges the updated Lagrangian formulation and the total Lagrangian formulation. The ARC framework serves as a generalized computational framework to handle both the computational infinitesimal and the finite deformation/strain deformation applications in a consistent and unified manner. In part II of the paper [1], we further extend the ARC framework to elasto-plasticity.
Imaging the North Anatolian Fault using the scattered teleseismic wavefield
NASA Astrophysics Data System (ADS)
Thompson, D. A.; Rost, S.; Houseman, G. A.; Cornwell, D. G.; Turkelli, N.; Teoman, U.; Kahraman, M.; Altuncu Poyraz, S.; Gülen, L.; Utkucu, M.; Frederiksen, A. W.; Rondenay, S.
2013-12-01
The North Anatolian Fault Zone (NAFZ) is a major continental strike-slip fault system, similar in size and scale to the San Andreas system, that extends ˜1200 km across Turkey. In 2012, a new multidisciplinary project (FaultLab) was instigated to better understand deformation throughout the entire crust in the NAFZ, in particular the expected transition from narrow zones of brittle deformation in the upper crust to possibly broader shear zones in the lower crust/upper mantle and how these features contribute to the earthquake loading cycle. This contribution will discuss the first results from the seismic component of the project, a 73 station network encompassing the northern and southern branches of the NAFZ in the Sakarya region. The Dense Array for North Anatolia (DANA) is arranged as a 6×11 grid with a nominal station spacing of 7 km, with a further 7 stations located outside of the main grid. With the excellent resolution afforded by the DANA network, we will present images of crustal structure using the technique of teleseismic scattering tomography. The method uses a full waveform inversion of the teleseismic scattered wavefield coupled with array processing techniques to infer the properties and location of small-scale heterogeneities (with scales on the order of the seismic wavelength) within the crust. We will also present preliminary results of teleseismic scattering migration, another powerful method that benefits from the dense data coverage of the deployed seismic network. Images obtained using these methods together with other conventional imaging techniques will provide evidence for how the deformation is distributed within the fault zone at depth, providing constraints that can be used in conjunction with structural analyses of exhumed fault segments and models of geodetic strain-rate across the fault system. By linking together results from the complementary techniques being employed in the FaultLab project, we aim to produce a comprehensive
NASA Technical Reports Server (NTRS)
Schmerr, Nicholas C.; Weber, Renee C.; Lin, Pei-Ying Patty; Thorne, Michael Scott; Garnero, Ed J.
2011-01-01
Lunar seismograms are distinctly different from their terrestrial counterparts. The Apollo lunar seismometers recorded moonquakes without distinct P- or S-wave arrivals; instead waves arrive as a diffuse coda that decays over several hours making the identification of body waves difficult. The unusual character of the lunar seismic wavefield is generally tied to properties of the megaregolith: it consists of highly fractured and broken crustal rock, the result of extensive bombardment of the Moon. The megaregolith extends several kilometers into the lunar crust, possibly into the mantle in some regions, and is covered by a thin coating of fine-scale dust. These materials possess very low seismic velocities that strongly scatter the seismic wavefield at high frequencies. Directly modeling the effects of the megaregolith to simulate an accurate lunar seismic wavefield is a challenging computational problem, owing to the inherent 3-D nature of the problem and the high frequencies (greater than 1 Hz) required. Here we focus on modeling the long duration code, studying the effects of the low velocities found in the megaregolith. We produce synthetic seismograms using 1-D slowness integration methodologies, GEMINI and reflectivity, and a 3-D Cartesian finite difference code, Wave Propagation Program, to study the effect of thin layers of low velocity on the surface of a planet. These codes allow us generate seismograms with dominant frequencies of approximately 1 Hz. For background lunar seismic structure we explore several models, including the recent model of Weber et al., Science, 2011. We also investigate variations in megaregolithic thickness, velocity, attenuation, and seismogram frequency content. Our results are compared to the Apollo seismic dataset, using both a cross correlation technique and integrated envelope approach to investigate coda decay. We find our new high frequency results strongly support the hypothesis that the long duration of the lunar seismic
Unsupervised pattern recognition in continuous seismic wavefield records using Self-Organizing Maps
NASA Astrophysics Data System (ADS)
Köhler, Andreas; Ohrnberger, Matthias; Scherbaum, Frank
2010-09-01
Modern acquisition of seismic data on receiver networks worldwide produces an increasing amount of continuous wavefield recordings. In addition to manual data inspection, seismogram interpretation requires therefore new processing utilities for event detection, signal classification and data visualization. The use of machine learning techniques automatises decision processes and reveals the statistical properties of data. This approach is becoming more and more important and valuable for large and complex seismic records. Unsupervised learning allows the recognition of wavefield patterns, such as short-term transients and long-term variations, with a minimum of domain knowledge. This study applies an unsupervised pattern recognition approach for the discovery, imaging and interpretation of temporal patterns in seismic array recordings. For this purpose, the data is parameterized by feature vectors, which combine different real-valued wavefield attributes for short time windows. Standard seismic analysis tools are used as feature generation methods, such as frequency-wavenumber, polarization and spectral analysis. We use Self-Organizing Maps (SOMs) for a data-driven feature selection, visualization and clustering procedure. The application to continuous recordings of seismic signals from an active volcano (Mount Merapi, Java, Indonesia) shows that volcano-tectonic and rockfall events can be detected and distinguished by clustering the feature vectors. Similar results are obtained in terms of correctly classifying events compared to a previously implemented supervised classification system. Furthermore, patterns in the background wavefield, that is the 24-hr cycle due to human activity, are intuitively visualized by means of the SOM representation. Finally, we apply our technique to an ambient seismic vibration record, which has been acquired for local site characterization. Disturbing wavefield patterns are identified which affect the quality of Love wave dispersion
Wavefield properties of a shallow long-period event and tremor at Kilauea Volcano, Hawaii
Saccorotti, G.; Chouet, B.; Dawson, P.
2001-01-01
The wavefields of tremor and a long-period (LP) event associated with the ongoing eruptive activity at Kilauea Volcano, Hawaii, are investigated using a combination of dense small-aperture (300 m) and sparse large-aperture (5 km) arrays deployed in the vicinity of the summit caldera. Measurements of azimuth and slowness for tremor recorded on the small-aperture array indicate a bimodal nature of the observed wavefield. At frequencies below 2 Hz, the wavefield is dominated by body waves impinging the array with steep incidence. These arrivals are attributed to the oceanic microseismic noise. In the 2-6 Hz band, the wavefield is dominated by waves propagating from sources located at shallow depths (<1 km) beneath the eastern edge of the Halemaumau pit crater. The hypocenter of the LP event, determined from frequency-slowness analyses combined with phase picks, appears to be located close to the source of tremor but at a shallower depth (<0.1 km). The wavefields of tremor and LP event are characterized by a complex composition of body and surface waves, whose propagation and polarization properties are strongly affected by topographic and structural features in the summit caldera region. Analyses of the directional properties of the wavefield in the 2-6 Hz band point to the directions of main scattering sources, which are consistent with pronounced velocity contrasts imaged in a high-resolution three-dimensional velocity model of the caldera region. The frequency and Q of the dominant peak observed in the spectra of the LP event may be explained as the dominant oscillation mode of a crack with scale length 20-100 m and aperture of a few centimeters filled with bubbly water. The mechanism driving the shallow tremor appears to be consistent with a sustained excitation originating in the oscillations of a bubbly cloud resulting from vesiculation and degassing in the magma. ?? 2001 Elsevier Science B.V. All rights reserved.
Delamination detection in foam core composite structures using transient flexural wavefields
NASA Astrophysics Data System (ADS)
Lamboul, B.; Osmont, D.
2016-03-01
This paper investigates a health monitoring technique for foam sandwich structures based on flexural wavefields imaged by Laser Doppler vibrometry. A study of calibrated artificial defect responses in harmonic regime demonstrates that the use of a low frequency regime (below 30 kHz) makes it possible to excite defects in their first flexural resonance modes. The analysis performed in harmonic regime is used to interpret signature patterns obtained with accumulated energy maps from transient wavefield recordings. The potential of the technique is demonstrated on a real impact-induced defect. The robustness of the method relatively to the excitation center frequency selection and to the presence of wave reverberation is demonstrated.
NASA Astrophysics Data System (ADS)
Hu, H.; Guan, W.; Wang, Z.
2014-12-01
Borehole electrokinetic wavefields have been theoretically simulated and experimentally recorded. However, it is still challenging to explain some of the signals in the full seismoeletric waveform. Similarly, while earthquake coseismic electric and magnetic signals were recorded and theoretically modeled, there are some basic questions to be answered regarding the formulation of the earthquake electrokinetic field. First, an electromagnetic signal appears at the same time in all recorded full waveforms when an acoustic wave is incident on the borehole wall or an interface between two porous media. Is it a traveling electromagnetic wave or a field? This is explained by a comparison between the waveforms obtained by solving the full Pride equations and those by a quasi-static approximation to the seismic-to-electric conversion. Second, a magnetic signal accompanies the borehole P-wave. Does that contradicts to Pride's prediction that no magnetic signal travels with a P-wave? We will show that the borehole P-wave consists of plane fast-P, slow-P and shear waves. It is the plane S-wave that brings about the magnetic field. Thirdly, it was proposed that there were no seismoelectric signal accompanying the collar wave during seismoelectric logging while drilling, because the electrokinetic conversion occurs only in the porous formation. Why there is an electric signal accompanies the acoustic collar wave? A detailed study of the acoustic field in the formation reveals that there is a wave propagates with the collar wave speed in the formation. This wave is present in the calculated full waveforms, either by the discrete wavenumber method or by the finite-difference-time-domain algorithm. That explains the existence of a noise signal with collar-wave speed in the full waveform of the electric field recorded during seismoelectric logging while drilling. Finally, an earthquake is usually modeled by a double couple in an elastic medium, and the displacement field is
NASA Astrophysics Data System (ADS)
Minato, S.; Ghose, R.
2013-12-01
The elastdynamic response of a fracture is often modeled using the linear-slip model (LSM) for the fracture compliance. In earlier theoretical and laboratory studies, the distribution of compliance along the plane of a fracture has generally been assumed to be homogeneous. However, naturally occurring fractures are spatially heterogeneous, with the microscale properties varying along the fracture plane. The spatial heterogeneity of the microscale parameters along the fracture plane, e.g., roughness, contact area and distribution of fluid filled aperture, controls significantly the mechanical and hydraulic response of a fracture. When the fracture compliance is spatially heterogeneous, an incident elastic wavefield will be scattered at the fracture plane. This scattered wavefield contains information of the spatial heterogeneity of fracture compliance. In this study, we show through numerical modeling that the scattered elastic wavefield is sensitive to the spatial heterogeneity in compliance distribution. We find that the back-scattered elastic wavefield from a spatially heterogeneous fracture appears as the coda of the specular reflection, with amplitude differing from that for a homogenous fracture compliance. An analysis of the scattered wavefield does reveal the spatial heterogeneity along the fracture plane. In order to estimate the spatially heterogeneous compliance distribution, we have developed an inversion scheme. The scheme has the following two steps: (1) extrapolating the recorded back-scattered elastic wavefield and estimating the stress field at the fracture plane, and (2) solving the boundary condition of LSM using the estimated stress field. We illustrate this new method through numerical examples mimicking laboratory-scale measurements (Figure). In the low frequency, the estimated compliance distribution is smooth and inaccurate because of the presence of the evanescent waves. However, at the peak frequency, the compliance distribution can be
Effect of strong elastic contrasts on the propagation of seismic wave in hard-rock environments
NASA Astrophysics Data System (ADS)
Saleh, R.; Zheng, L.; Liu, Q.; Milkereit, B.
2013-12-01
Understanding the propagation of seismic waves in a presence of strong elastic contrasts, such as topography, tunnels and ore-bodies is still a challenge. Safety in mining is a major concern and seismic monitoring is the main tool here. For engineering purposes, amplitudes (peak particle velocity/acceleration) and travel times of seismic events (mostly blasts or microseismic events) are critical parameters that have to be determined at various locations in a mine. These parameters are useful in preparing risk maps or to better understand the process of spatial and temporal stress distributions in a mine. Simple constant velocity models used for monitoring studies in mining, cannot explain the observed complexities in scattered seismic waves. In hard-rock environments modeling of elastic seismic wavefield require detailed 3D petrophysical, infrastructure and topographical data to simulate the propagation of seismic wave with a frequencies up to few kilohertz. With the development of efficient numerical techniques, and parallel computation facilities, a solution for such a problem is achievable. In this study, the effects of strong elastic contrasts such as ore-bodies, rough topography and tunnels will be illustrated using 3D modeling method. The main tools here are finite difference code (SOFI3D)[1] that has been benchmarked for engineering studies, and spectral element code (SPECFEM) [2], which was, developed for global seismology problems. The modeling results show locally enhanced peak particle velocity due to presence of strong elastic contrast and topography in models. [1] Bohlen, T. Parallel 3-D viscoelastic finite difference seismic modeling. Computers & Geosciences 28 (2002) 887-899 [2] Komatitsch, D., and J. Tromp, Introduction to the spectral-element method for 3-D seismic wave propagation, Geophys. J. Int., 139, 806-822, 1999.
Astillero, Antonio; Santos, Andrés
2005-09-01
In the preceding paper, we have conjectured that the main transport properties of a dilute gas of inelastic hard spheres (IHSs) can be satisfactorily captured by an equivalent gas of elastic hard spheres (EHSs), provided that the latter are under the action of an effective drag force and their collision rate is reduced by a factor (1+alpha)/2 (where alpha is the constant coefficient of normal restitution). In this paper we test the above expectation in a paradigmatic nonequilibrium state, namely, the simple or uniform shear flow, by performing Monte Carlo computer simulations of the Boltzmann equation for both classes of dissipative gases with a dissipation range 0.5 < or = alpha < or = 0.95 and two values of the imposed shear rate a. It is observed that the evolution toward the steady state proceeds in two stages: a short kinetic stage (strongly dependent on the initial preparation of the system) followed by a slower hydrodynamic regime that becomes increasingly less dependent on the initial state. Once conveniently scaled, the intrinsic quantities in the hydrodynamic regime depend on time, at a given value of alpha, only through the reduced shear rate a*(t) is proportional to a/square root(T(t)), until a steady state, independent of the imposed shear rate and of the initial preparation, is reached. The distortion of the steady-state velocity distribution from the local equilibrium state is measured by the shear stress, the normal stress differences, the cooling rate, the fourth and sixth cumulants, and the shape of the distribution itself. In particular, the simulation results seem to be consistent with an exponential overpopulation of the high-velocity tail. These properties are common to both the IHS and EHS systems. In addition, the EHS results are in general hardly distinguishable from the IHS ones if alpha approximately > 0.7, so that the distinct signature of the IHS gas (higher anisotropy and overpopulation) only manifests itself at relatively high
INVERSION OF FULL ACOUSTIC WAVEFIELD IN LOCAL HELIOSEISMOLOGY: A STUDY WITH SYNTHETIC DATA
Cobden, L. J.; Warner, M. R.; Tong, C. H.
2011-02-01
We present the first results from the inversion of full acoustic wavefield in the helioseismic context. In contrast to time-distance helioseismology, which involves analyzing the travel times of seismic waves propagating into the solar interior, wavefield tomography models both the travel times and amplitude variations present in the entire seismic record. Unlike the use of ray-based, Fresnel-zone, Born, or Rytov approximations in previous time-distance studies, this method does not require any simplifications to be made to the sensitivity kernel in the inversion. In this study, the acoustic wavefield is simulated for all iterations in the inversion. The sensitivity kernel is therefore updated while lateral variations in sound-speed structure in the model emerge during the course of the inversion. Our results demonstrate that the amplitude-based inversion approach is capable of resolving sound-speed structures defined by relatively sharp vertical and horizontal boundaries. This study therefore provides the foundation for a new type of subsurface imaging in local helioseismology that is based on the inversion of the entire seismic wavefield.
NASA Astrophysics Data System (ADS)
Schissele, E.; Cansi, Y.; Gaffet, S.
Many observations and studies as well as numerical simulations have been done in order to completely understand the whole seismogram recorded during an earthquake. At regional distances, the seismic wavefield is strongly influenced by crustal hetero- geneities. The primary wavefield constituted by Pn, Pg, Sn, Sg, Rg, Lg.... phases is diffracted and refracted by these heterogeneities and hence forms the coda of the seis- mogram. But the different mechanisms of propagation in a heterogeneous medium are not fully understood. The identification of the different phases contributing to the coda seems to be essential to progress in the comprehension of the seismic wavefield propagation. Seismic arrays are then well-adapted tools since they provide the spatio-temporal evo- lution of the wavefield. In 1998, 4 small-scales arrays were deployed for 2 months around the Annot region, located in the southern French Alps. Each array was constituted by 9 short-period seismometers, recording frequencies greater than 0.2 Hz. Its aperture was 250 meters, with a minimal distance between 2 adjacent sensors of 20 meters. That allows us to study the seismic wavefield for very low wavelength without any problem of spatial aliasing. It will be interesting to characterize in terms of wavefield deformation the signature of the different kinds of heterogeneities (fault system, topographic relief, impedance contrast...) surrounding this area. We expect the primary wavefield to be diffracted or refracted by all these heterogeneities. A time-frequency-wavenumber technique which allows us to characterize the whole coherent part of the energy which prop- agates through the seismic array has been derived. Such a characterization involves, for each coherent wavelet, an estimate of: (i) an arrival time and a frequency content and (ii) an azimuth and an apparent velocity. This way, the principal phases will be described. What will be more interesting, is the extraction of the deterministic part of the
Doyle, Heather; Lohfeld, Stefan; McHugh, Peter
2014-03-01
This study assesses the ability of finite element (FE) models to capture the mechanical behaviour of sintered orthopaedic scaffold materials. Individual scaffold struts were fabricated from a 50:50 wt% poly-ε-caprolactone (PCL)/β-tricalcium phosphate (β-TCP) blend, using selective laser sintering. The tensile elastic modulus of single struts was determined experimentally. High resolution FE models of single struts were generated from micro-CT scans (28.8 μm resolution) and an effective strut elastic modulus was calculated from tensile loading simulations. Three material assignment methods were employed: (1) homogeneous PCL elastic constants, (2) composite PCL/β-TCP elastic constants based on rule of mixtures, and (3) heterogeneous distribution of micromechanically-determined elastic constants. In comparison with experimental results, the use of homogeneous PCL properties gave a good estimate of strut modulus; however it is not sufficiently representative of the real material as it neglects the β-TCP phase. The rule of mixtures method significantly overestimated strut modulus, while there was no significant difference between strut modulus evaluated using the micromechanically-determined elastic constants and experimentally evaluated strut modulus. These results indicate that the multi-scale approach of linking micromechanical modelling of the sintered scaffold material with macroscale modelling gives an accurate prediction of the mechanical behaviour of the sintered structure. PMID:24057867
NASA Astrophysics Data System (ADS)
Gaudot, Ianis; Beucler, Éric; Mocquet, Antoine; Schimmel, Martin; Le Feuvre, Mathieu
2016-04-01
In order to detect possible signal redundancies in the seismic ambient wavefield, we develop a new method based on pairwise comparisons among a set of synchronous time-series. This approach is based on instantaneous phase coherence statistics. The first and second moments of the pairwise phase coherence distribution are used to characterize the phase randomness. Both theory and synthetic experiments show that, for perfect phase randomness, the theoretical values of the mean and variance are equal to 0 and 1 ‑ 2/π, respectively. As a consequence, any deviation from these values indicates the presence of a redundant phase in the raw continuous signal. Using the ergodicity property of a random signal, we split an initial time-series into a set of synchronous signals. This allows us to detect and to quantify the repetitiveness of any possible temporally persistent and spatially localized source, during a given period of observation. In the case of the detection of a redundant phase, individual coherences (one trace against all others) quantify the contribution of each time-series independently. A previously detected 26 s period microseismic source located near the Gulf of Guinea is used to illustrate one of the possible ways of handling phase coherence statistics. We use the continuous vertical component data recorded during the month of 2004 August by four broad-band stations of the Federation of Digital Seismography Network. To compute coherence statistics among a set composed of a sufficient number of synchronous traces, the raw seismic signal is split into 372 2-hr sliding time windows. Only the basic signal processing steps (including removing the mean, trend and the instrumental response) are applied. After bandpass filtering the data between 23 and 32 s periods, the 2-hr time-series are cross-correlated, leading to a set of 372 synchronous cross-correlations for each station pair. We observe that, for all station pairs, the mean overall coherence value is
Biskri, Zine Elabidine; Rached, Habib; Bouchear, Merzoug; Rached, Djamel
2014-04-01
The objective of this study is to investigate theoretically the structural, elastic and electronic properties of Lithium Disilicate (LD) crystal (Li2Si2O5), using the pseudo potential method based on Density Functional Theory (DFT) with the Local Density Approximation (LDA) and the Generalized Gradient Approximation (GGA). The calculated structural properties namely the equilibrium lattice parameters and cell volume are in good agreement with the available experimental results. However, for the LD crystal elastic moduli: Shear modulus G, Young's modulus E and Poisson's ratio ν we have found a discrepancy between our theoretical values and experimental ones reported in polycrystalline sample containing LD crystals. The calculated elastic properties show that LD is more rigid compared with other components. We also investigated the mechanical stability of Li2Si2O5 compound and we have noticed that this compound is stable against elastic deformations. On the basis of shear to bulk modulus ratio analysis, we inferred that Li2Si2O5 compound is brittle in nature. In order to complete the fundamental characteristics of this compound we have measured the elastic anisotropy. Our results for the energy band structure and Density of States (DOS) show that Li2Si2O5 compound has an insulator characteristic. PMID:24411692
Boundary element modeling of earthquake site effects including the complete incident wavefield
NASA Astrophysics Data System (ADS)
Kim, Kyoung-Tae
Numerical modeling of earthquake site effects in realistic, three-dimensional structures, including high frequencies, low surface velocities and surface topography, has not been possible simply because the amount of computer memory constrains the number of grid points available. In principle, this problem is reduced in the Boundary Element Method (BEM) since only the surface of the velocity discontinuity is discretized; wave propagation both inside and outside this boundary is computed analytically. Equivalent body forces are determined on the boundary by solving a matrix equation containing frequency-domain displacement and stress Green's functions from every point on the boundary to every other point. This matrix problem has imposed a practical limit on the size or maximum frequency of previous BEM models. Although the matrix can be quite large, it also seems to be fairly sparse. We have used iterative matrix algorithms of the PETSc package and direct solution algorithms of the ScaLAPACK on the massively parallel supercomputers at Cornell, San Diego and Michigan. Preconditioning has been applied using blockwise ILU decomposition for the iterative approach or LU decomposition for the direct approach. The matrix equation is solved using the GMRES method for the iterative approach and a tri-diagonal solver for the direct approach. Previous BEM applications typically have assumed a single, incident plane wave. However, it is clear that for more realistic ground motion simulations, we need to consider the complete incident wavefield. If we assume that the basin or three-dimensional structure of interest is embedded in a surrounding plane-layered medium, we may use the propagator matrix method to solve for the displacements and stresses at depth on the boundary. This is done in the frequency domain with integration over wavenumber so that all P, S, mode conversions, reverberations and surface waves are included. The Boundary Element Method succeeds in modeling
The Propagation of Seismic Waves in the Presence of Strong Elastic Property Contrasts
NASA Astrophysics Data System (ADS)
Saleh, R.; Jeyaraj, R.; Milkereit, B.; Liu, Q.; Valley, B.
2012-12-01
In an active underground mine there are many seismic activities taking place, such as seismic noises, blasts, tremors and microseismic events. In between the activities, the microseismic events are mainly used for monitoring purposes. The frequency content of microseismic events can be up to few KHz, which can result in wavelengths on the order of a few meters in hard rock environment. In an underground mine, considering the presence of both small wavelength and strong elastic contrasts, the simulation of seismic wave propagation is a challenge. With the recent availability of detailed 3D rock property models of mines, in addition to the development of efficient numerical techniques (such as Spectral Element Method (SEM)), and parallel computation facilities, a solution for such a problem is achievable. Most seismic wave scattering studies focus on large scales (>1 km) and weak elastic contrasts (velocity perturbations less than 10%). However, scattering in the presence of small-scale heterogeneities and large elastic contrasts is an area of ongoing research. In a mine environment, the presence of strong contrast discontinuities such as massive ore bodies, tunnels and infrastructure lead to discontinuities of displacement and/or stress tensor components, and have significant impact on the propagation of seismic waves. In order to obtain an accurate image of wave propagation in such a complex media, it is necessary to consider the presence of these discontinuities in numerical models. In this study, the effects of such a contrast are illustrated with 2D/3D modeling and compared with real broadband 3-component seismic data. The real broadband 3-component seismic data will be obtained in one of the Canadian underground mines in Ontario. One of the possible scenarios investigated in this study that may explain the observed complexity in seismic wavefield pattern in hard rock environments is the effect of near field displacements rather than far field. Considering the
Mitri, F. G.
2015-09-15
The standard Resonance Scattering Theory (RST) of plane waves is extended for the case of any two-dimensional (2D) arbitrarily-shaped monochromatic beam incident upon an elastic cylinder with arbitrary location using an exact methodology based on Graf’s translational addition theorem for the cylindrical wave functions. The analysis is exact as it does not require numerical integration procedures. The formulation is valid for any cylinder of finite size and material that is immersed in a nonviscous fluid. Partial-wave series expansions (PWSEs) for the incident, internal and scattered linear pressure fields are derived, and the analysis is further extended to obtain generalized expressions for the on-axis and off-axis acoustic radiation force components. The wave-fields are expressed using generalized PWSEs involving the beam-shape coefficients (BSCs) and the scattering coefficients of the cylinder. The off-axial BSCs are expressed analytically in terms of an infinite PWSE with emphasis on the translational offset distance d. Numerical computations are considered for a zeroth-order quasi-Gaussian beam chosen as an example to illustrate the analysis. Acoustic resonance scattering directivity diagrams are calculated by subtracting an appropriate background from the expression of the scattered pressure field. In addition, computations for the radiation force exerted on an elastic cylinder centered on the axis of wave propagation of the beam, and shifted off-axially are analyzed and discussed.
Elastic modeling and migration in Earth models
NASA Astrophysics Data System (ADS)
Filhocunha, Carlos Alves
Migration and inversion of marine seismic data using the elastic wave equation requires the transformation of the recorded pressure data into a vector particle-displacement field. This can be done easily when the recording geometry samples the wavefield both horizontally and vertically. However, only experimental surveys have cables located at different depths. Using a few assumptions, I derive a method for performing this transformation, which is applicable to standard surveys. The assumptions are: smooth water surface, cable nearly parallel to water surface, and perfect seismic-reflection at the water surface. Results in a realistic example, where these assumptions are only partially fulfilled, demonstrate that the method is robust. Elastic, reverse-time migration/inversion schemes in the space-time domain are usually implemented by finite-difference or finite-element methods. When imaging beyond structures, a dynamically accurate scheme must be used. For models characterized by layers with sharp boundaries, traditional finite-difference methods fail to correctly describe the dynamics of the propagation process. Failure comes from the lack of distinction between model and field variables; the same difference operator is applied to discontinuous (model) and continuous (wavefield) components. The problem is solved with a modified finite-difference scheme (dual-operator), which uses long operators for wave-fields, short operators for elastic parameters, Shoenberg-Muir (1989) equivalence relations, and a modified Virieux (1984) staggered grid scheme. Tests show that the dual-operator is dynamically more accurate than traditional finite-difference schemes and comparable to Haskell-Thomson schemes. In structurally complex media, accurate recovery of angle-dependent reflectivities requires elastic prestack migration. Mode separation can be done before or after depth extrapolation. Though more complex, the latter is more complete because it images mode-converted waves
Elastic Modeling and Migration in Earth Models
NASA Astrophysics Data System (ADS)
Cunha, Carlos Alves, Filho
Migration and inversion of marine seismic data using the elastic wave equation requires the transformation of the recorded pressure data into a vector particle-displacement field. This can be done easily when the recording geometry samples the wavefield both horizontally and vertically. However, only experimental surveys have cables located at different depths. Using a few assumptions, I derive a method for performing this transformation, which is applicable to standard surveys. The assumptions are: smooth water surface, cable nearly parallel to water surface, and perfect seismic -reflection at the water surface. Results in a realistic example, where these assumptions are only partially fulfilled, demonstrate that the method is robust. Elastic, reverse-time migration/inversion schemes in the space-time domain are usually implemented by finite -difference or finite-element methods. When imaging beyond structures, a dynamically accurate scheme must be used. For models characterized by layers with sharp boundaries traditional finite-difference methods fail to correctly describe the dynamics of the propagation process. Failure comes from the lack of distinction between model and field variables; the same difference operator is applied to discontinuous (model) and continuous (wavefield) components. The problem is solved with a modified finite-difference scheme (dual -operator), which uses long operators for wavefields, short operators for elastic parameters, Shoenberg-Muir (1989) equivalence relations and a modified Virieux (1984) staggered grid scheme. Tests show that the the dual-operator is dynamically more accurate than traditional finite-difference schemes and comparable to Haskell-Thomson schemes. In structurally complex media, accurate recovery of angle-dependent reflectivities requires elastic prestack migration. Mode separation can be done before or after depth extrapolation. Though more complex, the latter is more complete because it images mode
On the seismic noise wavefield in a weakly dissipative layered Earth
NASA Astrophysics Data System (ADS)
Lunedei, Enrico; Albarello, Dario
2009-06-01
Proceeding on the assumption that the seismic noise wavefield mainly consists of surface waves generated by a random distribution of independent, point-like and harmonic sources, we modelled a set of relationships that allows us to interpret empirical H/V and apparent dispersion curves in terms of the dynamic properties of a viscoelastic layered Earth. This formalization can be used to determine the effect of material damping on the expected noise wavefield in cases where the sources are relatively distant from the receivers and where the frequencies taken into account are higher than the fundamental resonance frequency of the local subsoil. Numerical experiments indicate that such effects are significant on the H/V ratio curve, whereas the apparent dispersion curve appears to be less strongly influenced by damping. This finding opens new perspectives for the use of noise measurements to infer the in situ anelastic properties of subsoil from passive seismic measurements made at the surface.
Michaels, T.E.; Michaels, J.E.; Mi, B.; Ruzzene, M.
2005-04-09
A methodology is presented for health monitoring and subsequent inspection of critical structures. Algorithms have been developed to detect and approximately locate damaged regions by analyzing signals recorded from a permanently mounted, sparse array of transducers. Followup inspections of suspected flaw locations are performed using a dual transducer ultrasonic approach where a permanently mounted transducer is the source and an externally scanned transducer is the receiver. Scan results are presented as snapshots of the propagating ultrasonic wavefield radiating out from the attached transducers. This method, referred to here as Acoustic Wavefield Imaging (AWI), provides an excellent visual representation of the interaction of propagating ultrasonic waves with the structure. Pre-flaw and post-flaw ultrasonic waveforms are analyzed from an aluminum plate specimen with artificially induced damage, and the AWI results show the location and spatial extent of all of the defects.
NASA Astrophysics Data System (ADS)
Sollberger, David; Schmelzbach, Cedric; Robertsson, Johan O. A.; Greenhalgh, Stewart A.; Nakamura, Yosio; Khan, Amir
2016-04-01
We present a new seismic velocity model of the shallow lunar crust, including, for the first time, shear wave velocity information. So far, the shear wave velocity structure of the lunar near-surface was effectively unconstrained due to the complexity of lunar seismograms. Intense scattering and low attenuation in the lunar crust lead to characteristic long-duration reverberations on the seismograms. The reverberations obscure later arriving shear waves and mode conversions, rendering them impossible to identify and analyze. Additionally, only vertical component data were recorded during the Apollo active seismic experiments, which further compromises the identification of shear waves. We applied a novel processing and analysis technique to the data of the Apollo 17 lunar seismic profiling experiment (LSPE), which involved recording seismic energy generated by several explosive packages on a small areal array of four vertical component geophones. Our approach is based on the analysis of the spatial gradients of the seismic wavefield and yields key parameters such as apparent phase velocity and rotational ground motion as a function of time (depth), which cannot be obtained through conventional seismic data analysis. These new observables significantly enhance the data for interpretation of the recorded seismic wavefield and allow, for example, for the identification of S wave arrivals based on their lower apparent phase velocities and distinct higher amount of generated rotational motion relative to compressional (P-) waves. Using our methodology, we successfully identified pure-mode and mode-converted refracted shear wave arrivals in the complex LSPE data and derived a P- and S-wave velocity model of the shallow lunar crust at the Apollo 17 landing site. The extracted elastic-parameter model supports the current understanding of the lunar near-surface structure, suggesting a thin layer of low-velocity lunar regolith overlying a heavily fractured crust of basaltic
NASA Astrophysics Data System (ADS)
Liu, Xin; Pavlis, Gary L.
2013-03-01
We develop a generic method to appraise the reliability of wavefield imaging methods and use it to validate some novel observations on the 410-km discontinuity. The core concept of the error appraisal method is to produce a simulated data set that replicates the geometry of the real data. Here we implemented two simulation methods: (1) flat layer primary P to S conversions, and (2) a point source scattering model for P to S conversion data based on the Born approximation and ray theory propagators. We show how the approach can be extended for any simulation algorithm. We apply this new approach to appraise recent results using a 3-D, three-component P to S conversion imaging method applied to data collected by the USArray. Multiple metrics show that the amplitude of P to S converted energy scattered from the 410-km discontinuity varies by 18 dB with a systematically lower amplitude in an irregular band running from Idaho through northern Arizona. In addition, we observe strong lateral changes in the ratio of amplitudes recovered on the radial versus the transverse component. We compute point resolution functions and a checkerboard test to demonstrate we can reliably recover relative amplitudes with a lateral scale of the order of 200 km and a vertical scale of approximately 10 km. Irregular coverage locally distorts the amplitudes recovered in the checkerboard, but a 156 km scale checkerboard pattern is recovered. Flat layer simulations show we can recover relative amplitudes to within a range of 1 dB and the reconstructed transverse to radial amplitude is everywhere less than 0.1. A model with north-south oriented ridges with a 3° wavelength and 12.5 km amplitude shows of the order of ±6 dB amplitude variations and small, but clear correlation of the transverse/radial amplitude ratio topography in the model. Finally, we model the 410-km discontinuity as a rough surface characterized by variations in amplitude and depth derived from the USArray data. The rough
Fracture-related wavefield polarization and seismic anisotropy across the Greendale Fault
NASA Astrophysics Data System (ADS)
Pischiutta, M.; Savage, M. K.; Holt, R. A.; Salvini, F.
2015-10-01
We investigate seismic signatures of fracturing in a newly ruptured strike-slip fault by determining the wavefield polarization in the New Zealand Canterbury Plains area and across the Greendale Fault, which was responsible for the 3 September 2010 Darfield Mw 7.1 earthquake. Previous studies suggested that fractured rocks in fault damage zones cause directional amplification and ground motion polarization in the fracture-perpendicular direction as an effect of stiffness anisotropy, and cause velocity anisotropy with shear wave velocity larger in the fracture-parallel component. An array of 14 stations was installed following the Darfield earthquake. We assess polarization both in the frequency and time domains through the individual-station horizontal-to-vertical spectral ratio and covariance matrix analysis, respectively, and compare the results to previously reported anisotropy measurements from shear wave splitting. Stations installed in the Canterbury Plains have an amplification peak between 0.1 and 0.3 Hz for both earthquakes and ambient noise. We relate the amplification to the resonance of a considerable thickness (c. 1 km) of soft sediments lying over the metamorphic bedrock. Analysis of seismic events revealed the existence of another peak in amplification between 2 and 5 Hz at two on-fault stations, which was not visible in the noise analysis. In contrast to the lower frequency peak, the ones between 2 and 5 Hz are more strongly anisotropic, attaining amplitudes up to a factor of 4 in the N52° direction. To interpret this effect we model the fracture pattern in the fault damage zone produced by the fault kinematics. We conclude that the horizontal polarization is orthogonal to extensional fractures, which predominate in the shallow layers (<2 km) with an expected strike of N139°. Fracture orientation is consistent with coseismic surface rupture observations, confirming the reliability of the model. S wave splitting is produced by velocity anisotropy
NASA Technical Reports Server (NTRS)
Wu, R. W.; Witmer, E. A.
1972-01-01
A user-oriented FORTRAN 4 computer program, called JET 3, is presented. The JET 3 program, which employs the spatial finite-element and timewise finite-difference method, can be used to predict the large two-dimensional elastic-plastic transient Kirchhoff-type deformations of a complete or partial structural ring, with various support conditions and restraints, subjected to a variety of initial velocity distributions and externally-applied transient forcing functions. The geometric shapes of the structural ring can be circular or arbitrarily curved and with variable thickness. Strain-hardening and strain-rate effects of the material are taken into account.
Quantitatively understanding the imprint of fractures in the seismic wave-field
NASA Astrophysics Data System (ADS)
Vsemirnova, Ekaterina; Roberts, Alan; Long, Jon; Jones, Richard; McCaffrey, Ken; Hobbs, Richard
2015-04-01
Understanding fracture connectivity in the shallow crust is of major importance for the development and production of hydrocarbon fields. Fracture datasets collected from wells have limited spatial coverage compared to remote sensing methods such as seismic imaging, Ground Penetrating Radar (GPR), electromagnetic recording, Terrestrial Laser Scanning (TLS), and Unmanned Aerial Vehicles ("drones"). In this study we focus on quantitatively understanding the imprint of several classes of realistic fracture network on the seismic wave-field. The thin, often rough sheet-like form of fractures poses challenges for reliable imaging of fracture networks using seismic methods, and the seismic response can be significantly altered by the highly variable dip of the fractures. A number of studies have been published showing the effect of the presence of simple fracture configurations on the synthetic seismic wave-field. At present, however, due to the inherent complexity of real fracture networks, there is limited understanding regarding the extraction of network characteristics from seismic data. Our work involves forward seismic wave-field simulation of a range of complex fracture networks derived from detailed quantitative characterisation of fractures in outcrop. We aim to build a library of calibrated examples from which to both develop understanding of the information contained in a seismic dataset related to the fracture network, and further research into the quantitative inversion and imaging of such information.
An analytical solution to separate P-waves and S-waves in the VSP wavefield
Amano, Hiroshi
1994-12-31
An analytical solution to separate P-waves and S-waves in the VSP wavefield is derived with combinations of the formal solution of a forward VSP modeling. Some practical applications of this method to synthetic seismograms and field data are investigated and evaluated. Little wave distortion is recognized and the weak wavefield masked by dominant wave trains can be extracted with this method. The decomposed wavefield is expressed in frequency-depth (f-z) domain as a linear combination of up to the third order differential of traces, which is approximated by trace difference sin the practical separation process. In general, five traces with single-component data are required in this process, but the same process is implemented with only three traces in the acoustic case. Two-trace extrapolation is applied to each edge of data gather in order to enhance the accuracy of trace difference. Since the formulas are developed in f-z domain, the influence of anelasticity is taken into account with simplicity and the calculation is carried out fast enough with the benefit of fast Fourier transform (FFT).
2013-01-01
Background The resistance of the bone against damage by repairing itself and adapting to environmental conditions is its most important property. These adaptive changes are regulated by physiological process commonly called the bone remodeling. Better understanding this process requires that we apply the theory of elastic-damage under the hypothesis of small displacements to a bone structure and see its mechanical behavior. Results The purpose of the present study is to simulate a two dimensional model of a proximal femur by taking into consideration elastic-damage and mechanical stimulus. Here, we present a mathematical model based on a system of nonlinear ordinary differential equations and we develop the variational formulation for the mechanical problem. Then, we implement our mathematical model into the finite element method algorithm to investigate the effect of the damage. Conclusion The results are consistent with the existing literature which shows that the bone stiffness drops in damaged bone structure under mechanical loading. PMID:23663260
NASA Astrophysics Data System (ADS)
Shin, Hosop; Park, Jonghyun; Han, Sangwoo; Sastry, Ann Marie; Lu, Wei
2015-03-01
The mechanical instability of the Solid Electrolyte Interphase (SEI) layer in lithium ion (Li-ion) batteries causes significant side reactions resulting in Li-ion consumption and cell impedance rise by forming further SEI layers, which eventually leads to battery capacity fade and power fade. In this paper, the composition-/structure-dependent elasticity of the SEI layer is investigated via Atomic Force Microscopy (AFM) measurements coupled with X-ray Photoelectron Spectroscopy (XPS) analysis, and atomistic calculations. It is observed that the inner layer is stiffer than the outer layer. The measured Young's moduli are mostly in the range of 0.2-4.5 GPa, while some values above 80 GPa are also observed. This wide variation of the observed elastic modulus is elucidated by atomistic calculations with a focus on chemical and structural analysis. The numerical analysis shows the Young's moduli range from 2.4 GPa to 58.1 GPa in the order of the polymeric, organic, and amorphous inorganic components. The crystalline inorganic component (LiF) shows the highest value (135.3 GPa) among the SEI species. This quantitative observation on the elasticity of individual components of the SEI layer must be essential to analyzing the mechanical behavior of the SEI layer and to optimizing and controlling it.
Homogenization of Electromagnetic and Seismic Wavefields for Joint Inverse Modeling
NASA Astrophysics Data System (ADS)
Newman, G. A.; Commer, M.; Petrov, P.; Um, E. S.
2011-12-01
) possibilities to image across multiple scale lengths, incorporating different types of geophysical data and attributes in the process. Important numerical details of 3D seismic wave field simulation in the Laplace-Fourier domain for both acoustic and elastic cases will also be discussed.
NASA Astrophysics Data System (ADS)
Husebye, E. S.; Kebeasy, T. R.
2004-12-01
The Luroy earthquake of August 31, 1819 with MS around 6.0 is, by many colleagues, rated as the largest in NW Europe in historical times (pre-1900) and even up to present. Local shaking manifestations were most spectacular with rock and mud avalanches, mast-high waves in nearby Rana fjord and even liquefaction was reported. Most surprisingly, at epicentral distances exceeding 100 km except for Stockholm 800 km away, very few macroseismic observations are available. Another peculiarity was the lack of any significant housing damage even in the Luroy parish itself. In a recent paper, we postulated that the earthquake was of moderate size, reestimated at MS = 5.1, but of shallow depth between 5 - 10 km causing the intense local shaking. In this article, we add a new dimension to the many of Luroy earthquake studies namely simulating the seismic wavefield response of Luroy itself and adjacent areas characterized by steep topographic relief. We use a 3D finite difference scheme and compute ground motion for a point source. We used a shear waves source with a focal depth of 5 km. Water covered areas are replaced by crystalline crust due to the dearth of dense bathymetric data. Main results are that the topography of the Luroy, close to the mountain peak at 685 meter, cause wavefield amplification by a factor of 20 and even stronger. Further away in the Rana fjord and surrounding areas, we also got strong amplification in particular where the relief is sharpest thus explaining triggering of avalanches in a quantitative manner. In other words, macroseismic observations would be biased upward due to the topographic focusing effects and unless properly corrected for also will increase the final earthquake magnitude estimate. We take these results to strongly support our claim that the historic Luroy earthquake was of moderate size of MS = 5.1 and not at MS = 6.0 class as claimed by many colleagues. Finally, downscaling of maximum earthquake magnitude would also lower the
Elastic Velocity Updating through Image-Domain Tomographic Inversion of Passive Seismic Data
NASA Astrophysics Data System (ADS)
Witten, B.; Shragge, J. C.
2014-12-01
Seismic monitoring at injection sites (e.g., CO2sequestration, waste water disposal, hydraulic fracturing) has become an increasingly important tool for hazard identification and avoidance. The information obtained from this data is often limited to seismic event properties (e.g., location, approximate time, moment tensor), the accuracy of which greatly depends on the estimated elastic velocity models. However, creating accurate velocity models from passive array data remains a challenging problem. Common techniques rely on picking arrivals or matching waveforms requiring high signal-to-noise data that is often not available for the magnitude earthquakes observed over injection sites. We present a new method for obtaining elastic velocity information from earthquakes though full-wavefield wave-equation imaging and adjoint-state tomography. The technique exploits images of the earthquake source using various imaging conditions based upon the P- and S-wavefield data. We generate image volumes by back propagating data through initial models and then applying a correlation-based imaging condition. We use the P-wavefield autocorrelation, S-wavefield autocorrelation, and P-S wavefield cross-correlation images. Inconsistencies in the images form the residuals, which are used to update the P- and S-wave velocity models through adjoint-state tomography. Because the image volumes are constructed from all trace data, the signal-to-noise in this space is increased when compared to the individual traces. Moreover, it eliminates the need for picking and does not require any estimation of the source location and timing. Initial tests show that with reasonable source distribution and acquisition array, velocity anomalies can be recovered. Future tests will apply this methodology to other scales from laboratory to global.
Visco-elastic controlled-source full waveform inversion without surface waves
NASA Astrophysics Data System (ADS)
Paschke, Marco; Krause, Martin; Bleibinhaus, Florian
2016-04-01
We developed a frequency-domain visco-elastic full waveform inversion for onshore seismic experiments with topography. The forward modeling is based on a finite-difference time-domain algorithm by Robertsson that uses the image-method to ensure a stress-free condition at the surface. The time-domain data is Fourier-transformed at every point in the model space during the forward modeling for a given set of frequencies. The motivation for this approach is the reduced amount of memory when computing kernels, and the straightforward implementation of the multiscale approach. For the inversion, we calculate the Frechet derivative matrix explicitly, and we implement a Levenberg-Marquardt scheme that allows for computing the resolution matrix. To reduce the size of the Frechet derivative matrix, and to stabilize the inversion, an adapted inverse mesh is used. The node spacing is controlled by the velocity distribution and the chosen frequencies. To focus the inversion on body waves (P, P-coda, and S) we mute the surface waves from the data. Consistent spatiotemporal weighting factors are applied to the wavefields during the Fourier transform to obtain the corresponding kernels. We test our code with a synthetic study using the Marmousi model with arbitrary topography. This study also demonstrates the importance of topography and muting surface waves in controlled-source full waveform inversion.
The Broad-Band Seismic Noise Wavefield at the Larderello-Travale Geothermal Field (Italy)
NASA Astrophysics Data System (ADS)
Zupo, M.; Saccorotti, G.; Piccinini, D.
2013-12-01
Cross-correlation of ambient noise wave-field between a pair of receivers (NCF), provides an estimate of the Green's Function between the two sites, thus allowing extraction of the associated group velocity dispersion curve. This is valid under the assumption that noise sources and/or scatterers are isotropically distributed and uncorrelated each other. These conditions are usually met once the cross-correlations are averaged over long time intervals. At frequencies lower than 1 Hz, ambient noise wavefield is essentially composed by surface waves that are mostly associated with oceanic sources; as a consequence, the noise wavefield may exhibit marked directional properties over short (day) to intermediate (weeks) time scales. A detailed assessment of the nature and duration of these sources is therefore required in order to define the optimal conditions for retrieving the Green's functions from NCF analysis. This study presents ambient noise analysis for the Larderello-Travale Geothermal Field (Italy). We use data collected by a temporary seismic array consisting of 20 broad-band instruments, with station spacing ranging from 6 to 50 Km. Below 1 Hz, the most energetic sources are those associated with both primary and secondary microseisms, with dominant spectral peaks spanning the 0.05-0.5 Hz frequency range. Focusing on the secondary microseism sources (f > 0.1 Hz), we test the validity of the isotropic-wavefield assumption by determining the kinematic properties of the wavefield using frequency-domain beamforming. For the May-November 2012 time span, our results show that the most energetic and persistent wavefield components propagate from WNW (Gulf of Marseille and Genova) and SW (Sardinia channel). In the late part of the observation period, additional wavefield components are observed to propagate from the NE-SE azimuthal range, corresponding to sources located throughout the Adriatic sea. This suggests that the conditions for the application of the NCF
NASA Astrophysics Data System (ADS)
Zubov, V. I.; Sanchez, J. F.; Tretiakov, N. P.; Yusef, A. E.
Based on the correlative method of an unsymmetrized self-consistent field,16-23 we have derived expressions for elastic constant tensors of strongly anharmonic crystals of cubic symmetry. Each isothermal elastic constant consists of four terms. The first one is the zeroth approximation containing the main anharmonicity (up to the fourth order). The second term is the quantum correction. It is important at temperatures below the De-bye characteristic temperature. Finally, the third and fourth terms are the perturbation theory corrections which take into account the influence of the correlations in atomic displacements from the lattice points and that of the high-order anharmonicity respectively. These corrections appear to be small up to the melting temperatures. It is sufficient for a personal computer to perform all our calculations with just a little computer time. A comparison with certain Monte Carlo simulations and with experimental data for Ar and Kr is made. For the most part, our results are between. The quasi-harmonic approximation fails at high temperatures, confirming once again the crucial role of strong anharmonicity.
NASA Astrophysics Data System (ADS)
de Hoop, Adrianus T.
2004-06-01
A uniqueness theorem for the (analytic or computational) time-domain modeling of the elastic wave motion in a scattering configuration that consists of inhomogeneous, anisotropic solids with arbitrary relaxation properties, occupying a bounded subdomain in an unbounded homogeneous, isotropic, perfectly elastic embedding, is presented. No direct time-domain uniqueness proof seems to exist for this kind of configuration. As an intermediate step, the one-to-one correspondence between the causal time-domain wavefield components and the constitutive material response functions on the one hand, and their time Laplace-transform counterparts for (a sequence of) real, positive values of the transform parameter on the other hand, seems a necessary tool. It is shown that such an approach leads to simple, explicit, sufficiency conditions on the inertial loss and compliance relaxation tensors describing the solid's constitutive behavior for uniqueness to hold. In it, the property of causality plays an essential role. In Christensen [Theory of Viscoelasticity-An Introduction (Academic, New York, 1971)] a similar approach is applied to the problem of uniqueness of the elastodynamic initial-/boundary-value problem associated with a viscoelastic object of bounded extent, the surface of which is subject to an admissible set of explicit boundary values. In the scattering configuration of unbounded extent, no explicit boundary values occur and the far-field compressional and shear wave radiation characteristics at ``infinity'' in the embedding play a key role in the proof.
NASA Astrophysics Data System (ADS)
Pusey, Jason L.; Yoo, Jin-Hyeong
2014-06-01
We document the design and preliminary numerical simulation study of a high fidelity model of Canid, a recently introduced bounding robot. Canid is a free-standing, power-autonomous quadrupedal machine constructed from standard commercially available electromechanical and structural elements, incorporating compliant C-shaped legs like those of the decade old RHex design, but departing from that standard (and, to the best of our knowledge, from any prior) robot platform in its parallel actuated elastic spine. We have used a commercial modeling package to develop a finite-element model of the actuated, cable-driven, rigid-plate-reinforced harness for the carbon-fiber spring that joins the robot's fore- and hind-quarters. We compare a numerical model of this parallel actuated elastic spine with empirical data from preliminary physical experiments with the most important component of the spine assembly: the composite leaf spring. Specifically, we report our progress in tuning the mechanical properties of a standard modal approximation to a conventional compliant beam model whose boundary conditions represent constraints imposed by the actuated cable driven vertebral plates that comprise the active control affordance over the spine. We conclude with a brief look ahead at near-term future experiments that will compare predictions of this fitted composite spring model with data taken from the physical spine flexed in isolation from the actuated harness.
Sahariah, Munima B; Ghosh, Subhradip; Singh, Chabungbam S; Gowtham, S; Pandey, Ravindra
2013-01-16
The structural stabilities, elastic, electronic and magnetic properties of the Heusler-type shape memory alloy Ni(2)FeGa are calculated using density functional theory. The volume conserving tetragonal distortion of the austenite Ni(2)FeGa find an energy minimum at c/a = 1.33. Metastable behaviour of the high temperature cubic austenite phase is predicted due to elastic softening in the [110] direction. Calculations of the total and partial magnetic moments show a dominant contribution from Fe atoms of the alloy. The calculated density of states shows a depression in the minority spin channel of the cubic Ni(2)FeGa just above the Fermi level which gets partially filled up in the tetragonal phase. In contrast to Ni(2)MnGa, the transition metal spin-down states show partial hybridization in Ni(2)FeGa and there is a relatively high electron density of states near the Fermi level in both phases. PMID:23186622
Wang, Yuxiang; Gerling, Gregory J.
2015-01-01
Our capability to discriminate object compliance is based on cues both tactile and proprioceptive, in addition to visual. To understand how the mechanics of the fingertip skin and bone might encode such information, we used finite element models to simulate the task of differentiating spherical indenters of radii (4, 6 and 8 mm) and elasticity (initial shear modulus of 10, 50 and 90 kPa). In particular, we considered two response variables, the strain energy density (SED) at the epidermal-dermal interface where Merkel cell end-organs of slowly adapting type I afferents reside, and the displacement of the fingertip bone necessary to achieve certain surface contact force. The former variable ties to tactile cues while the latter ties to proprioceptive cues. The results indicate that distributions of SED are clearly distinct for most combinations of object radii and elasticity. However, for certain combinations – e.g., between 4 mm spheres of 10 kPa and 8 mm of 90 kPa – spatial distributions of SED are nearly identical. In such cases where tactile-only cues are non-differentiable, we may rely on proprioceptive cues to discriminate compliance. PMID:25699293
NASA Astrophysics Data System (ADS)
Johnson, S. E.; Vel, S. S.; Cook, A. C.; Song, W. J.; Gerbi, C. C.; Okaya, D. A.
2014-12-01
Owing to the abundance of highly anisotropic minerals in the crust, the Voigt and Reuss bounds on the seismic velocities can be separated by more than 1 km/s. These bounds are determined by modal mineralogy and crystallographic preferred orientations (CPO) of the constituent minerals, but where the true velocities lie between these bounds is determined by other fabric parameters such as the shapes, shape-preferred orientations, and spatial arrangements of grains. Thus, the calculation of accurate bulk stiffness relies on explicitly treating the grain-scale heterogeneity, and the same principle applies at larger scales, for example calculating accurate bulk stiffness for a crustal volume with varying proportions and distributions of folds or shear zones. We have developed stand-alone GUI software - ESP Toolbox - for the calculation of 3D bulk elastic and seismic properties of heterogeneous and polycrystalline materials using image or EBSD data. The GUI includes a number of different homogenization techniques, including Voigt, Reuss, Hill, geometric mean, self-consistent and asymptotic expansion homogenization (AEH) methods. The AEH method, which uses a finite element mesh, is most accurate since it explicitly accounts for elastic interactions of constituent minerals/phases. The user need only specify the microstructure and material properties of the minerals/phases. We use the Toolbox to explore changes in bulk elasticity and related seismic anisotropy caused by specific variables, including: (a) the quartz alpha-beta phase change in rocks with varying proportions of quartz, (b) changes in modal mineralogy and CPO fabric that occur during progressive deformation and metamorphism, and (c) shear zones of varying thickness, abundance and geometry in continental crust. The Toolbox allows rapid sensitivity analysis around these and other variables, and the resulting bulk stiffness matrices can be used to populate volumes for synthetic wave propagation experiments that
Metzler, Adam M; Collis, Jon M
2013-04-01
Shallow-water environments typically include sediments containing thin or low-shear layers. Numerical treatments of these types of layers require finer depth grid spacing than is needed elsewhere in the domain. Thin layers require finer grids to fully sample effects due to elasticity within the layer. As shear wave speeds approach zero, the governing system becomes singular and fine-grid spacing becomes necessary to obtain converged solutions. In this paper, a seismo-acoustic parabolic equation solution is derived utilizing modified difference formulas using Galerkin's method to allow for variable-grid spacing in depth. Propagation results are shown for environments containing thin layers and low-shear layers. PMID:23556690
NASA Astrophysics Data System (ADS)
Kumar, Hemant; Er, Dequan; Dong, Liang; Li, Junwen; Shenoy, Vivek B.
2015-06-01
Recent technological advances in the isolation and transfer of different 2-dimensional (2D) materials have led to renewed interest in stacked Van der Waals (vdW) heterostructures. Interlayer interactions and lattice mismatch between two different monolayers cause elastic strains, which significantly affects their electronic properties. Using a multiscale computational method, we demonstrate that significant in-plane strains and the out-of-plane displacements are introduced in three different bilayer structures, namely graphene-hBN, MoS2-WS2 and MoSe2-WSe2, due to interlayer interactions which can cause bandgap change of up to ~300 meV. Furthermore, the magnitude of the elastic deformations can be controlled by changing the relative rotation angle between two layers. Magnitude of the out-of-plane displacements in graphene agrees well with those observed in experiments and can explain the experimentally observed bandgap opening in graphene. Upon increasing the relative rotation angle between the two lattices from 0° to 10°, the magnitude of the out-of-plane displacements decrease while in-plane strains peaks when the angle is ~6°. For large misorientation angles (>10°), the out-of-plane displacements become negligible. We further predict the deformation fields for MoS2-WS2 and MoSe2-WSe2 heterostructures that have been recently synthesized experimentally and estimate the effect of these deformation fields on near-gap states.
1-D seismic inversion of dual wavefield data. Part 2: A Gulf of Mexico example
Hildebrand, S.T. ); McMechan, G.A. . Center for Lithospheric Studies)
1994-05-01
Stable estimates for angle-dependent acoustic impedance are obtained from a dual wavefield data set that consists of pressure and vertical particle velocity measurements. Using these acoustic impedance estimates, a nonlocal inversion procedure predicts sediment material impedance, interval velocity, density, and their associated relative errors as a function of depth. The method is successfully applied to data from a marine seismic experiment conducted by Halliburton Geophysical Services in the Gulf of Mexico in 1989. The inverted values for the seafloor sediment density and interval velocity lie close to those obtained independently by other methods, implying that the inversion illustrated is valid for use with field data.
NASA Astrophysics Data System (ADS)
Van Der Kruk, J.; Yang, X.; Klotzsche, A.; von Hebel, C.; Busch, S.; Mester, A.; Huisman, J. A.; Vereecken, H.
2014-12-01
Ray-based or approximate forward modeling techniques have been often used to reduce the computational demands for inversion purposes. Due to increasing computational power and possible parallelization of inversion algorithms, accurate forward modeling can be included in advanced inversion approaches such that the full-wavefield content can be exploited. Here, recent developments of large-scale quantitative electromagnetic induction (EMI) inversion and full-waveform ground penetrating radar (GPR) inversions are discussed that yield higher resolution of quantitative medium properties compared to conventional approaches due to the use of accurate modeling tools that are based on Maxwell's equations. For a limited number of parameters, a combined global and local search using the simplex search algorithm or the shuffled complex evolution (SCE) can be used for inversion. Examples will be shown where calibrated large-scale multi-configuration EMI data measured with new generation multi-offset EMI systems are inverted for a layered electrical conductivity earth, and quantitative permittivity and conductivity values of a layered subsurface can be obtained using on-ground GPR full-waveform inversion that includes the estimation of the unknown source wavelet. For a large number of unknowns, gradient-based optimization methods are commonly used that need a good start model to prevent it from being trapped in a local minimum. Examples will be shown where the non-linearity invoked by the presence of high contrast media can be tamed by using a novel combined frequency-time-domain full-waveform inversion, and a low-velocity waveguide layer can be imaged by using crosshole GPR full-waveform inversion, after adapting the starting model using waveguide identification in the measured data. Synthetic data calculated using the inverted permittivity and conductivity models show similar amplitudes and phases as observed in the measured data, which indicates the reliability of the
NASA Technical Reports Server (NTRS)
Bartels, Robert E.
1999-01-01
This paper presents a modification of the spring analogy scheme which uses axial linear spring stiffness with selective spring stiffening/relaxation. An alternate approach to solving the geometric conservation law is taken which eliminates the need for storage of metric Jacobians at previous time steps. Efficiency and verification are illustrated with several unsteady 2-D airfoil Euler computations. The method is next applied to the computation of the turbulent flow about a 2-D airfoil and wing with two and three- dimensional moving spoiler surfaces, and the results compared with Benchmark Active Controls Technology (BACT) experimental data. The aeroelastic response at low dynamic pressure of an airfoil to a single large scale oscillation of a spoiler surface is computed. This study confirms that it is possible to achieve accurate solutions with a very large time step for aeroelastic problems using the fluid solver and aeroelastic integrator as discussed in this paper.
NASA Astrophysics Data System (ADS)
Cassereau, Didier; Mézière, Fabien; Muller, Marie; Bossy, Emmanuel; Derode, Arnaud
In this paper, we are interested in the 2D computation of the pressure scattered by an elliptic scatterer using a semi-analytical method based on a decomposition of the solutions on a basis of cylindrical waves. This approach is perfectly adapted to circular scatterers, and has been extended to scatterers of arbitrary shape [F. Chati et al. (2004)]. We will see that this extended formulation yields some very difficult numerical issues, particularly in our context of a flat and small elliptic scatterer. The use of arbitrary precision mathematics appears as a possible workaround, even if the cost in terms of the computation time may be prohibitive.
NASA Astrophysics Data System (ADS)
Gaudot, I.; Beucler, É.; Mocquet, A.; Schimmel, M.; Le Feuvre, M.
2016-02-01
In order to detect possible signal redundancies in the ambient seismic wavefield, we develop a new method based on pairwise comparisons among a set of synchronous time-series. This approach is based on instantaneous phase coherence statistics. The first and second moments of the pairwise phase coherence distribution are used to characterize the phase randomness. For perfect phase randomness, the theoretical values of the mean and variance are equal to 0 and √{1-2/π }, respectively. As a consequence, any deviation from these values indicates the presence of a redundant phase in the raw continuous signal. A previously detected microseismic source in the Gulf of Guinea is used to illustrate one of the possible ways of handling phase coherence statistics. The proposed approach allows us to properly localize this persistent source, and to quantify its contribution to the overall seismic ambient wavefield. The strength of the phase coherence statistics relies in its ability to quantify the redundancy of a given phase among a set of time-series with various useful applications in seismic noise-based studies (tomography and/or source characterization).
NASA Astrophysics Data System (ADS)
An, Yun-Kyu; Song, Homin; Sohn, Hoon
2014-09-01
This paper presents a wireless ultrasonic wavefield imaging (WUWI) technique for detecting hidden damage inside a steel box girder bridge. The proposed technique allows (1) complete wireless excitation of piezoelectric transducers and noncontact sensing of the corresponding responses using laser beams, (2) autonomous damage visualization without comparing against baseline data previously accumulated from the pristine condition of a target structure and (3) robust damage diagnosis even for real structures with complex structural geometries. First, a new WUWI hardware system was developed by integrating optoelectronic-based signal transmitting and receiving devices and a scanning laser Doppler vibrometer. Next, a damage visualization algorithm, self-referencing f-k filter (SRF), was introduced to isolate and visualize only crack-induced ultrasonic modes from measured ultrasonic wavefield images. Finally, the performance of the proposed technique was validated through hidden crack visualization at a decommissioned Ramp-G Bridge in South Korea. The experimental results reveal that the proposed technique instantaneously detects and successfully visualizes hidden cracks even in the complex structure of a real bridge.
Kim, Min Hyeok; Kim, Young Jin; Kim, Hee Ryung; Jeon, Tae-Joon; Choi, Jae Boong; Chung, Ka Young; Kim, Moon Ki
2016-01-01
Agonist-activated G protein-coupled receptors (GPCRs) interact with GDP-bound G protein heterotrimers (Gαβγ) promoting GDP/GTP exchange, which results in dissociation of Gα from the receptor and Gβγ. The GTPase activity of Gα hydrolyzes GTP to GDP, and the GDP-bound Gα interacts with Gβγ, forming a GDP-bound G protein heterotrimer. The G protein cycle is allosterically modulated by conformational changes of the Gα subunit. Although biochemical and biophysical methods have elucidated the structure and dynamics of Gα, the precise conformational mechanisms underlying the G protein cycle are not fully understood yet. Simulation methods could help to provide additional details to gain further insight into G protein signal transduction mechanisms. In this study, using the available X-ray crystal structures of Gα, we simulated the entire G protein cycle and described not only the steric features of the Gα structure, but also conformational changes at each step. Each reference structure in the G protein cycle was modeled as an elastic network model and subjected to normal mode analysis. Our simulation data suggests that activated receptors trigger conformational changes of the Gα subunit that are thermodynamically favorable for opening of the nucleotide-binding pocket and GDP release. Furthermore, the effects of GTP binding and hydrolysis on mobility changes of the C and N termini and switch regions are elucidated. In summary, our simulation results enabled us to provide detailed descriptions of the structural and dynamic features of the G protein cycle. PMID:27483005
Sewell, T. D.; Bedrov, D.; Menikoff, Ralph; Smith, G. D.
2001-01-01
Atomistic molecular dynamics simulations have been used to calculate isothermal elastic properties for {beta}-, {alpha}-, and {delta}-HMX. The complete elastic tensor for each polymorph was determined at room temperature and pressure via analysis of microscopic strain fluctuations using formalism due to Rahman and Parrinello [J. Chem. Phys. 76,2662 (1982)]. Additionally, the isothermal compression curve was computed for {beta}-HMX for 0 {le} p {le} 10.6 GPa; the bulk modulus K and its pressure derivative K{prime} were obtained from two fitting forms employed previously in experimental studies of the {beta}-HMX equation of state. Overall, the results indicate good agreement between the bulk modulus predicted from the measured and calculated compression curves. The bulk modulus determined directly from the elastic tensor of {beta}-HMX is in significant disagreement with the compression curve-based results. The explanation for this discrepancy is an area of current research.
Elastic membranes in confinement.
Bostwick, J B; Miksis, M J; Davis, S H
2016-07-01
An elastic membrane stretched between two walls takes a shape defined by its length and the volume of fluid it encloses. Many biological structures, such as cells, mitochondria and coiled DNA, have fine internal structure in which a membrane (or elastic member) is geometrically 'confined' by another object. Here, the two-dimensional shape of an elastic membrane in a 'confining' box is studied by introducing a repulsive confinement pressure that prevents the membrane from intersecting the wall. The stage is set by contrasting confined and unconfined solutions. Continuation methods are then used to compute response diagrams, from which we identify the particular membrane mechanics that generate mitochondria-like shapes. Large confinement pressures yield complex response diagrams with secondary bifurcations and multiple turning points where modal identities may change. Regions in parameter space where such behaviour occurs are then mapped. PMID:27440257
NASA Astrophysics Data System (ADS)
Moraru, Gheorghe; Mursa, Condrat
2006-12-01
In this book we present the basic concepts of the theory of elasticity: stress and deformation states (plane and three-dimensional) and generalized Hooke's law. We present a number of problems which have applications in strength analysis. The book includes a synthesis of the theory of elasticity and modern methods of applied mathematics. This book is designed for students, post graduate students and specialists in strength analysis. the book contains a number of appendixes which includes: elements of matrix-calculation, concepts of tensorial calculation, the Fourier transform, the notion of improper integrals,singular and hypersingular integrals, generalized functions, the Dirac Delta function
On the anisotropic elastic properties of hydroxyapatite.
NASA Technical Reports Server (NTRS)
Katz, J. L.; Ukraincik, K.
1971-01-01
Experimental measurements of the isotropic elastic moduli on polycrystalline specimens of hydroxyapatite and fluorapatite are compared with elastic constants measured directly from single crystals of fluorapatite in order to derive a set of pseudo single crystal elastic constants for hydroxyapatite. The stiffness coefficients thus derived are given. The anisotropic and isotropic elastic properties are then computed and compared with similar properties derived from experimental observations of the anisotropic behavior of bone.
Development of Active Seismic Vector-Wavefield Imaging Technology for Geothermal Applications
B. A. Hardage; J. L. Simmons, Jr.; M. DeAngelo
1999-10-01
This report describes the development and testing of vector-wavefield seismic sources that can generate shear (S) waves that may be valuable in geothermal exploration and reservoir characterization. Also described is a 3-D seismic data-processing effort to create images of Rye Patch geothermal reservoir from 3-D sign-bit data recorded over the geothermal prospect. Two seismic sources were developed and tested in this study that can be used to illuminate geothermal reservoirs with S-waves. The first was an explosive package that generates a strong, azimuth-oriented, horizontal force vector when deployed in a conventional shot hole. This vector-explosive source has never been available to industry before. The second source was a dipole formed by operating two vertical vibrators in either a force or phase imbalance. Field data are shown that document the strong S-wave modes generated by these sources.
NASA Astrophysics Data System (ADS)
Yoshimitsu, Nana; Furumura, Takashi; Maeda, Takuto
2016-09-01
The coda part of a waveform transmitted through a laboratory sample should be examined for the high-resolution monitoring of the sample characteristics in detail. However, the origin and propagation process of the later phases in a finite-sized small sample are very complicated with the overlap of multiple unknown reflections and conversions. In this study, we investigated the three-dimensional (3D) geometric effect of a finite-sized cylindrical sample to understand the development of these later phases. This study used 3D finite difference method simulation employing a free-surface boundary condition over a curved model surface and a realistic circular shape of the source model. The simulated waveforms and the visualized 3D wavefield in a stainless steel sample clearly demonstrated the process of multiple reflections and the conversions of the P and S waves at the side surface as well as at the top and bottom of the sample. Rayleigh wave propagation along the curved side boundary was also confirmed, and these waves dominate in the later portion of the simulated waveform with much larger amplitudes than the P and S wave reflections. The feature of the simulated waveforms showed good agreement with laboratory observed waveforms. For the simulation, an introduction of an absorbing boundary condition at the top and bottom of the sample made it possible to efficiently separate the contribution of the vertical and horizontal boundary effects in the simulated wavefield. This procedure helped to confirm the additional finding of vertically propagating multiple surface waves and their conversion at the corner of the sample. This new laboratory-scale 3D simulation enabled the appearance of a variety of geometric effects that constitute the later phases of the transmitted waves.
Saccorotti, G.; Chouet, B.; Dawson, P.
2003-01-01
The properties of the surface wavefield at Kilauea Volcano are analysed using data from small-aperture arrays of short-period seismometers deployed in and around the Kilauea caldera. Tremor recordings were obtained during two Japan-US cooperative experiments conducted in 1996 and 1997. The seismometers were deployed in three semi-circular arrays with apertures of 300, 300 and 400 m, and a linear array with length of 1680 m. Data are analysed using a spatio-temporal correlation technique well suited for the study of the stationary stochastic wavefields of Rayleigh and Love waves associated with volcanic activity and scattering sources distributed in and around the summit caldera. Spatial autocorrelation coefficients are obtained as a function of frequency and are inverted for the dispersion characteristics of Rayleigh and Love waves using a grid search that seeks phase velocities for which the L-2 norm between data and forward modelling operators is minimized. Within the caldera, the phase velocities of Rayleigh waves range from 1400 to 1800 m s-1 at 1 Hz down to 300-400 m s-1 at 10 Hz, and the phase velocities of Love waves range from 2600 to 400 m s-1 within the same frequency band. Outside the caldera, Rayleigh wave velocities range from 1800 to 1600 m s-1 at 1 Hz down to 260-360 m s-1 at 10 Hz, and Love wave velocities range from 600 to 150 m s-1 within the same frequency band. The dispersion curves are inverted for velocity structure beneath each array, assuming these dispersions represent the fundamental modes of Rayleigh and Love waves. The velocity structures observed at different array sites are consistent with results from a recent 3-D traveltime tomography of the caldera region, and point to a marked velocity discontinuity associated with the southern caldera boundary.
NASA Astrophysics Data System (ADS)
Lunedei, Enrico; Albarello, Dario
2016-03-01
Synthetic dispersion curves are here computed in the frame of an ambient-vibration full-wavefield model, which relies on the description of both ambient-vibration ground displacement and its sources as stochastic fields defined on the Earth's surface, stationary in time and homogeneous in space. In this model, previously developed for computing synthetic Horizontal-to-Vertical Spectral Ratio curves, the power spectral density function and the spatial autocorrelation of the displacement are naturally described as functions of the power spectral density function of the generating forces and of the subsoil properties (via the relevant Green's function), by also accounting for spatial correlation of these forces. Dispersion curves are computed from the displacement power spectral density function and from the spatial autocorrelation according with the well-known f-k and SPAC techniques, respectively. Two examples illustrate the way this new ambient-vibration model works, showing its possible use in better understanding the role of the surface waves in forming the dispersion curves, as well as its capability to capture some remarkable experimental findings.
NASA Technical Reports Server (NTRS)
Oline, L.; Medaglia, J.
1972-01-01
The dynamic finite element method was used to investigate elastic stress waves in a plate. Strain displacement and stress strain relations are discussed along with the stiffness and mass matrix. The results of studying point load, and distributed load over small, intermediate, and large radii are reported. The derivation of finite element matrices, and the derivation of lumped and consistent matrices for one dimensional problems with Laplace transfer solutions are included. The computer program JMMSPALL is also included.
Elasticity of crystalline molecular explosives
Hooks, Daniel E.; Ramos, Kyle J.; Bolme, C. A.; Cawkwell, Marc J.
2015-04-14
Crystalline molecular explosives are key components of engineered explosive formulations. In precision applications a high degree of consistency and predictability is desired under a range of conditions to a variety of stimuli. Prediction of behaviors from mechanical response and failure to detonation initiation and detonation performance of the material is linked to accurate knowledge of the material structure and first stage of deformation: elasticity. The elastic response of pentaerythritol tetranitrate (PETN), cyclotrimethylene trinitramine (RDX), and cyclotetramethylene tetranitramine (HMX), including aspects of material and measurement variability, and computational methods are described in detail. Experimental determinations of elastic tensors are compared, andmore » an evaluation of sources of error is presented. Furthermore, computed elastic constants are also compared for these materials and for triaminotrinitrobenzene (TATB), for which there are no measurements.« less
Elasticity of crystalline molecular explosives
Hooks, Daniel E.; Ramos, Kyle J.; Bolme, C. A.; Cawkwell, Marc J.
2015-04-14
Crystalline molecular explosives are key components of engineered explosive formulations. In precision applications a high degree of consistency and predictability is desired under a range of conditions to a variety of stimuli. Prediction of behaviors from mechanical response and failure to detonation initiation and detonation performance of the material is linked to accurate knowledge of the material structure and first stage of deformation: elasticity. The elastic response of pentaerythritol tetranitrate (PETN), cyclotrimethylene trinitramine (RDX), and cyclotetramethylene tetranitramine (HMX), including aspects of material and measurement variability, and computational methods are described in detail. Experimental determinations of elastic tensors are compared, and an evaluation of sources of error is presented. Furthermore, computed elastic constants are also compared for these materials and for triaminotrinitrobenzene (TATB), for which there are no measurements.
1-D seismic inversion of dual wavefield data. Part 1: Nonuniqueness and stability
Hildebrand, S.T. ); McMechan, G.A. . Center for Lithospheric Studies)
1994-05-01
The inversion problem for determining seismic impedance is nonunique and nonstable because of limited recording aperture, data bandwidth, and data noise. For large reflection angles, small errors in the reflection coefficients give rise to arbitrarily large errors in the seismic impedance estimates. Spatial resolution of the seismic impedance response is controlled by the dominant wavelength corresponding to the source time wavelet; aperture limitations control the resolution of material impedance and interval velocity. Analysis of a linearization-approximation approach shows that this method degenerates into a single-parameter estimator for material impedance when using only small-offset data and for velocity when using only far-offset data. A nonlocal inversion method is introduced to estimate the material impedance and interval velocity by exploring interval velocity space and computing an associated variance estimate surface. Using this method, the resolution of the material impedance and compressional and shear interval velocities is shown to be poor in the elastic case because of a ''valley'' feature in the variance estimate surface; in the acoustic problem, resolution of the material impedance and interval velocity is excellent.
NASA Astrophysics Data System (ADS)
Khaniani, Hassan
boundary condition of the wave equation is set up along reflection surfaces. Hence, the surface integral Kirchhoff approximation is used as a mathematical framework instead of the volume integral of the Born approximation. In addition, I study the feasibility of iterative coupling of ray theory with the Kirchhoff approximation for inversion. For the amplitude considerations, the direct relationship between the scattering potential of the Born approximation with the reflectivity function of the asymptotic Kirchhoff approximation for elastic waves is used. Therefore, I use the linearized Zoeppritz approximation of Aki and Richards (1980) for computation of the forward modeling and migration operators as well as gradient function from Amplitude vs Offset (AVO) inversion. The multiparameter elastic inversion approach is applicable to all types of reflected wavefields such as P-to-P, P-to-S, S-to-S and S-to-P. Traveltime estimation of forward modeling and migration/inversion operators are based on the DSR equation. All operators involved in inversion, including the background model for DSR and AVO are updated at each iteration. The migration/inversion procedure maps the mode converted waves to the traveltime of incident waves which fixes the registration problem of events that travel from source to scatter point. The inversion of the reflected P-to-P and P-to-S synthetic and field data are provided for the numerical examples. This approach is applicable for complex structures however, to estimate the traveltime of scatterpoints, ray tracing can be added to the algorithm. For such a medium, the scatterpoint traveltime approximations from the PSTM, is compared to the PSDM approach using numerical analysis of ray- and FDTD-based modeling. In part of this thesis, I further improve the conventional velocity analysis of Common Scatter Point (CSP) gathers by including the tilt effects. I show that travel time response of scatter points beneath a dipping interface experiences an
NASA Technical Reports Server (NTRS)
Roskam, J.; Lan, C.; Mehrotra, S.
1972-01-01
The computer program used to determine the rigid and elastic stability derivatives presented in the summary report is listed in this appendix along with instructions for its use, sample input data and answers. This program represents the airplane at subsonic and supersonic speeds as (a) thin surface(s) (without dihedral) composed of discrete panels of constant pressure according to the method of Woodward for the aerodynamic effects and slender beam(s) for the structural effects. Given a set of input data, the computer program calculates an aerodynamic influence coefficient matrix and a structural influence coefficient matrix.
Mitri, F G
2015-09-01
The classical Resonance Scattering Theory (RST) for plane waves in acoustics is generalized for the case of a 2D arbitrarily-shaped beam incident upon an elastic cylinder with arbitrary location that is immersed in a nonviscous fluid. The formulation is valid for an elastic (or viscoelastic) cylinder (or a cylindrical shell, a layered cylinder/shell, or a multilayered cylindrical shell, etc.) of any size and material. Partial-wave series expansions (PWSEs) for the incident, internal and scattered fields are derived, and numerical examples illustrate the theory. The wave-fields are expressed using a generalized PWSE involving the beam-shape coefficients (BSCs) and the scattering coefficients of the cylinder. When the beam is shifted off the center of the cylinder, the off-axial BSCs are evaluated by performing standard numerical integration. Acoustic resonance scattering directivity diagrams are calculated by subtracting an appropriate background from the expression of the scattered pressure field. The properties related to the arbitrary scattering of a zeroth-order quasi-Gaussian cylindrical beam (chosen as an example) by an elastic brass cylinder centered on the axis of wave propagation of the beam, and shifted off-axially are analyzed and discussed. Moreover, the total and resonance backscattering form function moduli are numerically computed, and the results discussed with emphasis on the contribution of the surface waves circumnavigating the cylinder circular surface to the resonance backscattering. Furthermore, the analysis is extended to derive general expressions for the axial and transverse acoustic radiation force functions for the cylinder in any 2D beam of arbitrary shape. Examples are provided for a zeroth-order quasi Gaussian cylindrical beam with different waist. Potential applications are in underwater and physical acoustics, however, ongoing research in biomedical ultrasound, non-destructive evaluation, imaging, manufacturing, instrumentation, and
Subcritical scattering from buried elastic shells.
Lucifredi, Irena; Schmidt, Henrik
2006-12-01
Buried objects have been largely undetectable by traditional high-frequency sonars due to their insignificant bottom penetration. Further, even a high grazing angle sonar approach is vastly limited by the coverage rate dictated by the finite water depth, making the detection and classification of buried objects using low frequency, subcritical sonar an interesting alternative. On the other hand, such a concept would require classification clues different from the traditional high-resolution imaging and shadows to maintain low false alarm rates. A potential alternative, even for buried targets, is classification based on the acoustic signatures of man-made elastic targets. However, the elastic responses of buried and proud targets are significantly different. The objective of this work is to identify, analyze, and explain some of the effects of the sediment and the proximity of the seabed interface on the scattering of sound from completely and partially buried elastic shells. The analysis was performed using focused array processing of data from the GOATS98 experiment carried out jointly by MIT and SACLANTCEN, and a new hybrid modeling capability combining a virtual source-or wave-field superposition-approach with an exact spectral integral representation of the Green's functions for a stratified ocean waveguide, incorporating all multiple scattering between the object and the seabed. Among the principal results is the demonstration of the significant role of structural circumferential waves in converting incident, evanescent waves into backscattered body waves, emanating to the receivers at supercritical grazing angles, in effect making the target appear closer to the sonar than predicted by traditional ray theory. PMID:17225387
Full-Wavefield Modeling and Pre-Stack Depth Migration of Common-Source Seismic Data
NASA Astrophysics Data System (ADS)
Chen, How-Wei
1992-01-01
The goal of solving geophysical problems can be thought of as a data mapping or transform procedure. Through various techniques, the observed seismic data can be transformed into the solution domain to estimate the Earth properties. Seismic wave field simulation is a forward process used to synthesize the seismic responses of an Earth model. Seismic wave field imaging is an inverse process used to estimate the Earth parameters from observed seismic data. In this dissertation, finite-difference and pseudo-spectral computations, in two- and three-dimensional space, are used for full wave field simulations and imaging of common-source data. Numerical simulation is developed for seismic sources and multi-attribute wave fields in two-dimensional acoustic and elastic media. P- and S-waves can be primarily separated in the resulting seismograms by vector operators in simulated surface survey, Vertical Seismic Profile (VSP) and cross-hole recording geometries. Three-component displacement seismograms can be approximately simulated by treating the acoustic field as a scalar potential field. The algorithm is applied to a complex multi-component, multi-offset walkaway circular VSP data from offshore California. Numerical modeling of large-scale, wide-aperture 3-D seismic data volumes is performed using a 3-D pseudo-spectral approach. Asymmetrical source and 3-D wave propagation effects in physical model data are identified and interpreted through iterative numerical modeling. Prestack reverse-time migration algorithms based on finite-difference and pseudo-spectral wave field extrapolators are developed for acoustic media in two- and three-dimensions. The excitation time imaging condition is computed by ray tracing and by finite-difference solution of the Eikonal equation. I generalize the concept of reverse-time migration and apply it for the correction of near-surface static effects. The feasibility of using very large scale, very wide-aperture 3-D seismic data recorded on a
Array analysis of regional Pn and Pg wavefields from the Nevada Test Site
Leonard, M.A. . Dept. of Geology and Geophysics Lawrence Berkeley Lab., CA )
1991-06-01
Small-aperture high-frequency seismic arrays with dimensions of a few kilometers or less, can improve our ability to seismically monitor compliance with a low-yield Threshold Test Ban Treaty. This work studies the characteristics and effectiveness of array processing of the regional Pn and Pg wavefields generated by underground nuclear explosions at the Nevada Test Site. Waveform data from the explosion HARDIN (m{sub b} = 5.5) is recorded at a temporary 12-element, 3-component, 1.5 km-aperture array sited in an area of northern Nevada. The explosions VILLE (m{sub b} = 4.4) and SALUT (m{sub b} = 5.5) are recorded at two arrays sited in the Mojave desert, one a 96-element vertical-component 7 km-aperture array and the other a 155-element vertical-component 4 km-aperture array. Among the mean spectra for the m{sub b} = 5.5 events there are significant differences in low-frequency spectral amplitudes between array sites. The spectra become nearly identical beyond about 6 Hz. Spectral ratios are used to examine seismic source properties and the partitioning of energy between Pn and Pg. Frequency-wavenumber analysis at the 12-element array is used to obtain estimates of signal gain, phase velocity, and source azimuth. This analysis reveals frequency-dependent biases in velocity and azimuth of the coherent Pn and Pg arrivals. Signal correlation, the principal factor governing array performance, is examined in terms of spatial coherence estimates. The coherence is found to vary between the three sites. In all cases the coherence of Pn is greater than that for Pg. 81 refs., 92 figs., 5 tabs.
The Utility of the Extended Images in Ambient Seismic Wavefield Migration
NASA Astrophysics Data System (ADS)
Girard, A. J.; Shragge, J. C.
2015-12-01
Active-source 3D seismic migration and migration velocity analysis (MVA) are robust and highly used methods for imaging Earth structure. One class of migration methods uses extended images constructed by incorporating spatial and/or temporal wavefield correlation lags to the imaging conditions. These extended images allow users to directly assess whether images focus better with different parameters, which leads to MVA techniques that are based on the tenets of adjoint-state theory. Under certain conditions (e.g., geographical, cultural or financial), however, active-source methods can prove impractical. Utilizing ambient seismic energy that naturally propagates through the Earth is an alternate method currently used in the scientific community. Thus, an open question is whether extended images are similarly useful for ambient seismic migration processing and verifying subsurface velocity models, and whether one can similarly apply adjoint-state methods to perform ambient migration velocity analysis (AMVA). Herein, we conduct a number of numerical experiments that construct extended images from ambient seismic recordings. We demonstrate that, similar to active-source methods, there is a sensitivity to velocity in ambient seismic recordings in the migrated extended image domain. In synthetic ambient imaging tests with varying degrees of error introduced to the velocity model, the extended images are sensitive to velocity model errors. To determine the extent of this sensitivity, we utilize acoustic wave-equation propagation and cross-correlation-based migration methods to image weak body-wave signals present in the recordings. Importantly, we have also observed scenarios where non-zero correlation lags show signal while zero-lags show none. This may be a valuable missing piece for ambient migration techniques that have yielded largely inconclusive results, and might be an important piece of information for performing AMVA from ambient seismic recordings.
A hybrid method for the computation of quasi-3D seismograms.
NASA Astrophysics Data System (ADS)
Masson, Yder; Romanowicz, Barbara
2013-04-01
The development of powerful computer clusters and efficient numerical computation methods, such as the Spectral Element Method (SEM) made possible the computation of seismic wave propagation in a heterogeneous 3D earth. However, the cost of theses computations is still problematic for global scale tomography that requires hundreds of such simulations. Part of the ongoing research effort is dedicated to the development of faster modeling methods based on the spectral element method. Capdeville et al. (2002) proposed to couple SEM simulations with normal modes calculation (C-SEM). Nissen-Meyer et al. (2007) used 2D SEM simulations to compute 3D seismograms in a 1D earth model. Thanks to these developments, and for the first time, Lekic et al. (2011) developed a 3D global model of the upper mantle using SEM simulations. At the local and continental scale, adjoint tomography that is using a lot of SEM simulation can be implemented on current computers (Tape, Liu et al. 2009). Due to their smaller size, these models offer higher resolution. They provide us with images of the crust and the upper part of the mantle. In an attempt to teleport such local adjoint tomographic inversions into the deep earth, we are developing a hybrid method where SEM computation are limited to a region of interest within the earth. That region can have an arbitrary shape and size. Outside this region, the seismic wavefield is extrapolated to obtain synthetic data at the Earth's surface. A key feature of the method is the use of a time reversal mirror to inject the wavefield induced by distant seismic source into the region of interest (Robertsson and Chapman 2000). We compute synthetic seismograms as follow: Inside the region of interest, we are using regional spectral element software RegSEM to compute wave propagation in 3D. Outside this region, the wavefield is extrapolated to the surface by convolution with the Green's functions from the mirror to the seismic stations. For now, these
Full Elastic Waveform Search Engine for Near Surface Imaging
NASA Astrophysics Data System (ADS)
Zhang, J.; Zhang, X.
2014-12-01
For processing land seismic data, the near-surface problem is often very complex and may severely affect our capability to image the subsurface. The current state-of-the-art technology for near surface imaging is the early arrival waveform inversion that solves an acoustic wave-equation problem. However, fitting land seismic data with acoustic wavefield is sometimes invalid. On the other hand, performing elastic waveform inversion is very time-consuming. Similar to a web search engine, we develop a full elastic waveform search engine that includes a large database with synthetic elastic waveforms accounting for a wide range of interval velocity models in the CMP domain. With each CMP gather of real data as an entry, the search engine applies Multiple-Randomized K-Dimensional (MRKD) tree method to find approximate best matches to the entry in about a second. Interpolation of the velocity models at CMP positions creates 2D or 3D Vp, Vs, and density models for the near surface area. The method does not just return one solution; it gives a series of best matches in a solution space. Therefore, the results can help us to examine the resolution and nonuniqueness of the final solution. Further, this full waveform search method can avoid the issues of initial model and cycle skipping that the method of full waveform inversion is difficult to deal with.
NASA Astrophysics Data System (ADS)
Williams, Westin B.; Michaels, Jennifer E.; Michaels, Thomas E.
2016-04-01
Detection, localization, and characterization of impact damage in composites using in situ transducers are important objectives for the aerospace industry to both reduce maintenance costs and prevent failures. A network of piezoelectric transducers spatially distributed over an area of interest is one practical configuration for utilizing guided waves to accomplish these objectives. Detecting and localizing barely visible impact damage with such a sparse array has been demonstrated in prior work, and improvements in localization were demonstrated by incorporating fairly crude estimates of scattering patterns in the imaging algorithms. Here we obtain more estimates of scattering patterns from a simulated defect by employing baseline subtraction of wavefield data recorded in a circle centered at the scatterer. Scattering patterns are estimated from the wavefield residual signals before and after simulated damage is introduced and the estimated scattering patterns are then incorporated into sparse array imaging via the minimum variance imaging method. Images created with different scattering patterns are compared and the efficacy of the methodology is assessed.
Loewenthal, M.; Loseke, K.; Dow, T.A.; Scattergood, R.O.
1988-12-01
Elastic emission polishing, also called elastic emission machining (EEM), is a process where a stream of abrasive slurry is used to remove material from a substrate and produce damage free surfaces with controlled surface form. It is a noncontacting method utilizing a thick elasto-hydrodynamic film formed between a soft rotating ball and the workpiece to control the flow of the abrasive. An apparatus was built in the Center, which consists of a stationary spindle, a two-axis table for the workpiece, and a pump to circulate the working fluid. The process is controlled by a programmable computer numerical controller (CNC), which presently can operate the spindle speed and movement of the workpiece in one axis only. This apparatus has been used to determine material removal rates on different material samples as a function of time, utilizing zirconium oxide (ZrO{sub 2}) particles suspended in distilled water as the working fluid. By continuing a study of removal rates the process should become predictable, and thus create a new, effective, yet simple tool for ultra-precision mechanical machining of surfaces.
NASA Astrophysics Data System (ADS)
Hoshiba, M.
2013-05-01
In this presentation, I propose a new approach for real-time prediction of seismic ground motion which is applicable to Earthquake Early Waning (EEW), in which hypocentral location and magnitude are not required. . Many methods of EEW are based on a network method in which hypocenter and magnitude (source parameters) are quickly determined, and then the ground motions are predicted, and warnings are issued depending on the strength of the predicted ground motion. In this method, it is necessary to determine the hypocenter and magnitude at first, and error of the source parameters leads directly to the error of the prediction. It is not easy to take the effects of rupture directivity and source extent into account, and it is impossible to fully reproduce the current wavefield from the interpreted source parameters. Time evolutional prediction is a method in which future wavefield is iteratively predicted from the wavefield at the certain time, that is u(x, t+Δt)=P(u(x, t)), where u is the wave motion at location x at lapse time t, and P is the prediction operator. Future wave motion, u(x, t+Δt), is predicted from the distribution of the current wave motion u(x, t) using P. For P, finite difference technique or boundary integral equation method, such as Kirchhoff integral, is used. Kirchhoff integral is qualitatively approximated by Huygens principle. The real time monitoring of wavefield are important for this method, but it is possible to predict ground motion without a hypocentral location and magnitude. In the time evolutional prediction, determination of detailed distribution of current wavefield is an important key, so that dense seismic observation network is required. Data assimilation is a technique to produce artificially denser network, which is widely used for numerical weather forecast and oceanography. Distribution of current wave motion is estimated from not only the current real observation of u(xi, t) where xi is the location of the i-th site, but
NASA Astrophysics Data System (ADS)
Reading, A. M.; Heckscher, N.; Graham, L.; Arroucau, P.; Rawlinson, N.
2009-12-01
Techniques using ambient seismic energy have become powerful tools in the imaging of the Earth's lithosphere. In this work we investigate the potential for using ambient energy to image the upper crust at the highest possible resolution. Such images are sought-after in deducing regional-scale 3D geological structure. They are critical to activities which use the Earth's crust (e.g. hot dry rock geothermal energy production, carbon sequestration) as well as placing fundamental constraints on the crustal architecture which defines tectonic structure and evolution. We have deployed a line of variously spaced seismic stations in a region of relatively well-constrained regional geology in Eastern Tasmania. Located south of mainland Australia, in the 'roaring 40's' southerly latitudes, the ambient seismic energy is of considerable amplitude. We are able to characterise the ambient seismic wavefield: investigating the influence of station deployment geometry, deployment and processing protocols, and ocean and atmospheric conditions on the amplitude and frequency content of the signals derived from correlated waveforms. We also investigate the potential of using diffuse seismic sources from a highway and railway using ambient energy techniques. We find that, in this ocean-dominated island setting, a correlated signal of sufficient strength to model for structure is obtained in a few days. The relationship between station separation and the dominant wavelengths in the correlated signals places a minimum limit on station separations which are usable with standard modelling techniques. Hence, in this environment, crustal imaging may be best carried out using deployments with overlapping, frequently moved (or numerous and short-deployed), sets of station pairs. The potential for using cultural noise sources is limited by the dominance of the natural noise sources in the ambient wavefield.
Elastic Wave Radiation from a Line Source of Finite Length
Aldridge, D.F.
1998-11-04
Straightforward algebraic expressions describing the elastic wavefield produced by a line source of finite length are derived in circular cylindrical coordinates. The surrounding elastic medium is assumed to be both homogeneous and isotropic, anc[ the source stress distribution is considered axisymmetic. The time- and space-domain formulae are accurate at all distances and directions from the source; no fa-field or long-wavelength assumptions are adopted for the derivation. The mathematics yield a unified treatment of three different types of sources: an axial torque, an axial force, and a radial pressure. The torque source radiates only azirnuthally polarized shear waves, whereas force and pressure sources generate simultaneous compressional and shear radiation polarized in planes containing the line source. The formulae reduce to more familiar expressions in the two limiting cases where the length of the line source approaches zero and infinity. Far-field approximations to the exact equations indicate that waves radiated parallel to the line source axI.s are attenuated relative to those radiated normal to the axis. The attenuation is more severe for higher I?equencies and for lower wavespeeds. Hence, shear waves are affected more than compressional waves. This fi-equency- and directiondependent attenuation is characterized by an extremely simple mathematical formula, and is readily apparent in example synthetic seismograms.
Polycrystalline gamma plutonium's elastic moduli versus temperature
Migliori, Albert; Betts, J; Trugman, A; Mielke, C H; Mitchell, J N; Ramos, M; Stroe, I
2009-01-01
Resonant ultrasound spectroscopy was used to measure the elastic properties of pure polycrystalline {sup 239}Pu in the {gamma} phase. Shear and longitudinal elastic moduli were measured simultaneously and the bulk modulus was computed from them. A smooth, linear, and large decrease of all elastic moduli with increasing temperature was observed. They calculated the Poisson ratio and found that it increases from 0.242 at 519 K to 0.252 at 571 K. These measurements on extremely well characterized pure Pu are in agreement with other reported results where overlap occurs.
Elastically Decoupling Dark Matter
NASA Astrophysics Data System (ADS)
Kuflik, Eric; Perelstein, Maxim; Lorier, Nicolas Rey-Le; Tsai, Yu-Dai
2016-06-01
We present a novel dark matter candidate, an elastically decoupling relic, which is a cold thermal relic whose present abundance is determined by the cross section of its elastic scattering on standard model particles. The dark matter candidate is predicted to have a mass ranging from a few to a few hundred MeV, and an elastic scattering cross section with electrons, photons and/or neutrinos in the 10-3- 1 fb range.
Elastically Decoupling Dark Matter.
Kuflik, Eric; Perelstein, Maxim; Lorier, Nicolas Rey-Le; Tsai, Yu-Dai
2016-06-01
We present a novel dark matter candidate, an elastically decoupling relic, which is a cold thermal relic whose present abundance is determined by the cross section of its elastic scattering on standard model particles. The dark matter candidate is predicted to have a mass ranging from a few to a few hundred MeV, and an elastic scattering cross section with electrons, photons and/or neutrinos in the 10^{-3}-1 fb range. PMID:27314712
Boyd, O.S.
2006-01-01
We have created a second-order finite-difference solution to the anisotropic elastic wave equation in three dimensions and implemented the solution as an efficient Matlab script. This program allows the user to generate synthetic seismograms for three-dimensional anisotropic earth structure. The code was written for teleseismic wave propagation in the 1-0.1 Hz frequency range but is of general utility and can be used at all scales of space and time. This program was created to help distinguish among various types of lithospheric structure given the uneven distribution of sources and receivers commonly utilized in passive source seismology. Several successful implementations have resulted in a better appreciation for subduction zone structure, the fate of a transform fault with depth, lithospheric delamination, and the effects of wavefield focusing and defocusing on attenuation. Companion scripts are provided which help the user prepare input to the finite-difference solution. Boundary conditions including specification of the initial wavefield, absorption and two types of reflection are available. ?? 2005 Elsevier Ltd. All rights reserved.
Soukup, Darin J; Odom, Robert I; Park, Jeffrey
2013-07-01
Theoretical and numerical results are presented for modal characteristics of the seismo-acoustic wavefield in anisotropic range-independent media. General anisotropy affects the form of the elastic-stiffness tensor, particle-motion polarization, the frequency and angular dispersion curves, and introduces near-degenerate modes. Horizontally polarized particle motion (SH) cannot be ignored when anisotropy is present for low-frequency modes having significant bottom interaction. The seismo-acoustic wavefield has polarizations in all three coordinate directions even in the absence of any scattering or heterogeneity. Even weak anisotropy may have a significant impact on seismo-acoustic wave propagation. Unlike isotropic and transversely isotropic media with a vertical symmetry axis where acoustic signals comprise P-SV modes alone (in the absence of any scattering), tilted TI media allow both quasi-P-SV and quasi-SH modes to carry seismo-acoustic energy. Discrete modes for an anisotropic medium are best described as generalized P-SV-SH modes with polarizations in all three Cartesian directions. Conversion to SH is a loss that will mimic acoustic attenuation. An in-water explosion will excite quasi-SH. PMID:23862797
NASA Astrophysics Data System (ADS)
Boyd, Oliver S.
2006-03-01
We have created a second-order finite-difference solution to the anisotropic elastic wave equation in three dimensions and implemented the solution as an efficient Matlab script. This program allows the user to generate synthetic seismograms for three-dimensional anisotropic earth structure. The code was written for teleseismic wave propagation in the 1-0.1 Hz frequency range but is of general utility and can be used at all scales of space and time. This program was created to help distinguish among various types of lithospheric structure given the uneven distribution of sources and receivers commonly utilized in passive source seismology. Several successful implementations have resulted in a better appreciation for subduction zone structure, the fate of a transform fault with depth, lithospheric delamination, and the effects of wavefield focusing and defocusing on attenuation. Companion scripts are provided which help the user prepare input to the finite-difference solution. Boundary conditions including specification of the initial wavefield, absorption and two types of reflection are available.
NASA Astrophysics Data System (ADS)
Cherry, Matthew R.; Aldrin, John C.; Boehnlein, Thomas; Blackshire, James L.
2013-01-01
In this work, a combined grid/FEM method that is capable of using parallelepiped or tetragonal mesh elements, as well as a combination of the two, is investigated. A formulation was developed that leverages the architecture of GPUs with irregular grids to efficiently address complex structures and heterogeneous materials. Benchmark studies are presented comparing the computational time, memory requirements, and simulation accuracy for GPU and CPU solvers with several challenge NDE problems.
Linear elastic fracture mechanics primer
NASA Astrophysics Data System (ADS)
Wilson, Christopher D.
1992-07-01
This primer is intended to remove the blackbox perception of fracture mechanics computer software by structural engineers. The fundamental concepts of linear elastic fracture mechanics are presented with emphasis on the practical application of fracture mechanics to real problems. Numerous rules of thumb are provided. Recommended texts for additional reading, and a discussion of the significance of fracture mechanics in structural design are given. Griffith's criterion for crack extension, Irwin's elastic stress field near the crack tip, and the influence of small-scale plasticity are discussed. Common stress intensities factor solutions and methods for determining them are included. Fracture toughness and subcritical crack growth are discussed. The application of fracture mechanics to damage tolerance and fracture control is discussed. Several example problems and a practice set of problems are given.
Elastic Flows Of Ellipsoidal Particles
NASA Astrophysics Data System (ADS)
Campbell, Charles S.
2009-06-01
Granular flow rheology can be divided into two global regimes, the Elastic, which is dominated by force chains and the inertial which are nearly free of force chains. The propensity of a material to form force chains is strongly influenced by particle shape. This paper is an attempt to assess the effect of particle shape on flow regime transitions, through computer simulations of shear flow of ellipsoidal particles. On one hand, the results show that at a given concentration, ellipsoidal particles generate smaller quasistatic stress than spheres, likely a result of their ability to form denser static packings. But at the same time, large aspect ratio ellipsoids more readily form force chains and demonstrate Elastic behavior at smaller concentrations than spheres.
Linear elastic fracture mechanics primer
NASA Technical Reports Server (NTRS)
Wilson, Christopher D.
1992-01-01
This primer is intended to remove the blackbox perception of fracture mechanics computer software by structural engineers. The fundamental concepts of linear elastic fracture mechanics are presented with emphasis on the practical application of fracture mechanics to real problems. Numerous rules of thumb are provided. Recommended texts for additional reading, and a discussion of the significance of fracture mechanics in structural design are given. Griffith's criterion for crack extension, Irwin's elastic stress field near the crack tip, and the influence of small-scale plasticity are discussed. Common stress intensities factor solutions and methods for determining them are included. Fracture toughness and subcritical crack growth are discussed. The application of fracture mechanics to damage tolerance and fracture control is discussed. Several example problems and a practice set of problems are given.
NASA Astrophysics Data System (ADS)
Perez, G. L.; Larour, E. Y.; Halkides, D. J.; Cheng, D. L. C.
2015-12-01
The Virtual Ice Sheet Laboratory(VISL) is a Cryosphere outreach effort byscientists at the Jet Propulsion Laboratory(JPL) in Pasadena, CA, Earth and SpaceResearch(ESR) in Seattle, WA, and the University of California at Irvine (UCI), with the goal of providing interactive lessons for K-12 and college level students,while conforming to STEM guidelines. At the core of VISL is the Ice Sheet System Model(ISSM), an open-source project developed jointlyat JPL and UCI whose main purpose is to model the evolution of the polar ice caps in Greenland and Antarctica. By using ISSM, VISL students have access tostate-of-the-art modeling software that is being used to conduct scientificresearch by users all over the world. However, providing this functionality isby no means simple. The modeling of ice sheets in response to sea and atmospheric temperatures, among many other possible parameters, requiressignificant computational resources. Furthermore, this service needs to beresponsive and capable of handling burst requests produced by classrooms ofstudents. Cloud computing providers represent a burgeoning industry. With majorinvestments by tech giants like Amazon, Google and Microsoft, it has never beeneasier or more affordable to deploy computational elements on-demand. This isexactly what VISL needs and ISSM is capable of. Moreover, this is a promisingalternative to investing in expensive and rapidly devaluing hardware.
Elastic internal flywheel gimbal
Rabenhorst, D.W.
1981-01-13
An elastic joint mounting and rotatably coupling a rotary inertial energy storage device or flywheel, to a shaft, the present gimbal structure reduces vibration and shock while allowing precession of the flywheel without the need for external gimbal mounts. The present elastic joint usually takes the form of an annular elastic member either integrally formed into the flywheel as a centermost segment thereof or attached to the flywheel or flywheel hub member at the center thereof, the rotary shaft then being mounted centrally to the elastic member.
NASA Technical Reports Server (NTRS)
Wu, R. W. H.; Stagliano, T. R.; Witmer, E. A.; Spilker, R. L.
1978-01-01
These structural ring deflections lie essentially in one plane and, hence, are called two-dimensional (2-d). The structural rings may be complete or partial; the former may be regarded as representing a fragment containment ring while the latter may be viewed as a 2-d fragment-deflector structure. These two types of rings may be either free or supported in various ways (pinned-fixed, locally clamped, elastic-foundation supported, mounting-bracket supported, etc.). The initial geometry of each ring may be circular or arbitrarily curved; uniform-thickness or variable-thickness rings may be analyzed. Strain-hardening and strain-rate effects of initially-isotropic material are taken into account. An approximate analysis utilizing kinetic energy and momentum conservation relations is used to predict the after-impact velocities of each fragment and of the impact-affected region of the ring; this procedure is termed the collision-imparted velocity method (CIVM) and is used in the CIVM-JET 5 B program. This imparted-velocity information is used in conjunction with a finite-element structural response computation code to predict the transient, large-deflection, elastic-plastic responses of the ring. Similarly, the equations of motion of each fragment are solved in small steps in time. Provisions are made in the CIVM-JET 5B code to analyze structural ring response to impact attack by from 1 to 3 fragments, each with its own size, mass, translational velocity components, and rotational velocity. The effects of friction between each fragment and the impacted ring are included.
Estimation of In vivo Cancellous Bone Elasticity
NASA Astrophysics Data System (ADS)
Otani, Takahiko; Mano, Isao; Tsujimoto, Toshiyuki; Yamamoto, Tadahito; Teshima, Ryota; Naka, Hiroshi
2009-07-01
The effect of decreasing bone density (a symptom of osteoporosis) is greater for cancellous bone than for dense cortical bone, because cancellous bone is metabolically more active. Therefore, the bone density or bone mineral density of cancellous bone is generally used to estimate the onset of osteoporosis. Elasticity or elastic constant is a fundamental mechanical parameter and is directly related to the mechanical strength of bone. Accordingly, elasticity is a preferable parameter for assessing fracture risk. A novel ultrasonic bone densitometer LD-100 has been developed to determine the mass density and elasticity of cancellous bone with a spatial resolution comparable to that of peripheral quantitative computed tomography. Bone density and bone elasticity are evaluated using ultrasonic parameters based on fast and slow waves in cancellous bone by modeling the ultrasonic wave propagation path. Elasticity is deduced from the measured bone density and the propagation speed of the fast wave. Thus, the elasticity of cancellous bone is approximately expressed by a cubic equation of bone density.
Elastic properties of minerals
Aleksandrov, K.S.; Prodaivoda, G.T.
1993-09-01
Investigations of the elastic properties of the main rock-forming minerals were begun by T.V. Ryzhova and K.S. Aleksandrov over 30 years ago on the initiative of B.P. Belikov. At the time, information on the elasticity of single crystals in general, and especially of minerals, was very scanty. In the surveys of that time there was information on the elasticity of 20 or 30 minerals. These, as a rule, did not include the main rock-forming minerals; silicates were represented only by garnets, quartz, topaz, tourmaline, zircon, beryl, and staurolite, which are often found in nature in the form of large and fairly high-quality crystals. Then and even much later it was still necessary to prove a supposition which now seems obvious: The elastic properties of rocks, and hence the velocities of elastic (seismic) waves in the earth`s crust, are primarily determined by the elastic characteristics of the minerals composing these rocks. Proof of this assertion, with rare exceptions of mono-mineralic rocks (marble, quartzite, etc.) cannot be obtained without information on the elasticities of a sufficiently large number of minerals, primarily framework, layer, and chain silicates which constitute the basis of most rocks. This also served as the starting point and main problem of the undertakings of Aleksandrov, Ryzhova, and Belikov - systematic investigations of the elastic properties of minerals and then of various rocks. 108 refs., 7 tabs.
Postinstability models in elasticity
NASA Technical Reports Server (NTRS)
Zak, M.
1984-01-01
It is demonstrated that the instability caused by the failure of hyperbolicity in elasticity and associated with the problem of unpredictability in classical mechanics expresses the incompleteness of the original model of an elastic medium. The instability as well as the ill-posedness of the Cauchy problem are eliminated by reformulating the original model.
Peng, Qing; De, Suvranu
2014-10-21
Silicane is a fully hydrogenated silicene-a counterpart of graphene-having promising applications in hydrogen storage with capacities larger than 6 wt%. Knowledge of its elastic limit is critical in its applications as well as tailoring its electronic properties by strain. Here we investigate the mechanical response of silicane to various strains using first-principles calculations based on density functional theory. We illustrate that non-linear elastic behavior is prominent in two-dimensional nanomaterials as opposed to bulk materials. The elastic limits defined by ultimate tensile strains are 0.22, 0.28, and 0.25 along armchair, zigzag, and biaxial directions, respectively, an increase of 29%, 33%, and 24% respectively in reference to silicene. The in-plane stiffness and Poisson ratio are reduced by a factor of 16% and 26%, respectively. However, hydrogenation/dehydrogenation has little effect on its ultimate tensile strengths. We obtained high order elastic constants for a rigorous continuum description of the nonlinear elastic response. The limitation of second, third, fourth, and fifth order elastic constants are in the strain range of 0.02, 0.08, and 0.13, and 0.21, respectively. The pressure effect on the second order elastic constants and Poisson's ratio were predicted from the third order elastic constants. Our results could provide a safe guide for promising applications and strain-engineering the functions and properties of silicane monolayers. PMID:25190587
NASA Astrophysics Data System (ADS)
Stachnik, Joshua C.
Surface wave analysis of both earthquake and ambient noise seismic data from arrays of broadband seismic stations provides new high resolution images of shear wave velocity of the crust and upper mantle in western North America. In the Yellowstone Hotspot region, new constraints are shown on the high velocity midcrustal layer of the eastern Snake River Plain that represents approximately 10 km of magmatic thickening and subsequent forcing of lower crustal outflow. In the Coast Mountain Batholith area of western British Columbia, the lack of a significant region of high velocities in the lower crust indicates that the foundering of negatively buoyant eclogitic lower crust has been efficient. A high resolution shear velocity model of the Sierra Nevada batholith region finds crustal thickening beneath the batholith, sinking material beneath the central Sierras with adjacent upwelling asthenospheric mantle, and the new image suggests that the Isabella (San Joaquin Valley) anomaly has a quasi-planar NW-SE striking geometry perhaps more consistent with being a Monterey plate slab remnant than an eclogite dominated feature. In addition to the surface wave results, new constraints are found on the sharpness of the 410-km velocity discontinuity via the wave-field continuation approach applied to five regional earthquakes in western North America. The 410-km discontinuity gradient ranges from 7-25 km, indicating the presence of water atop the transition zone in the two regions with large 410 km discontinuity gradient widths.
Proton Nucleus Elastic Scattering Data.
Energy Science and Technology Software Center (ESTSC)
1993-08-18
Version 00 The Proton Nucleus Elastic Scattering Data file PNESD contains the numerical data and the related bibliography for the differential elastic cross sections, polarization and integral nonelastic cross sections for elastic proton-nucleus scattering.
Wavefield Analysis of Rayleigh Waves for Near-Surface Shear-Wave Velocity
NASA Astrophysics Data System (ADS)
Zeng, Chong
2011-12-01
Shear (S)-wave velocity is a key property of near-surface materials and is the fundamental parameter for many environmental and engineering geophysical studies. Directly acquiring accurate S-wave velocities from a seismic shot gather is usually difficult due to the poor signal-to-noise ratio. The relationship between Rayleigh-wave phase velocity and frequency has been widely utilized to estimate the S-wave velocities in shallow layers using the multichannel analysis of surface waves (MASW) technique. Hence, Rayleigh wave is a main focus of most near-surface seismic studies. Conventional dispersion analysis of Rayleigh waves assumes that the earth is laterally homogeneous and the free surface is horizontally flat, which limits the application of surface-wave methods to only 1D earth models or very smooth 2D models. In this study I extend the analysis of Rayleigh waves to a 2D domain by employing the 2D full elastic wave equation so as to address the lateral heterogeneity problem. I first discuss the accurate simulation of Rayleigh waves through finite-difference method and the boundary absorbing problems in the numerical modeling with a high Poisson's ratio (> 0.4), which is a unique near-surface problem. Then I develop an improved vacuum formulation to generate accurate synthetic seismograms focusing on Rayleigh waves in presence of surface topography and internal discontinuities. With these solutions to forward modeling of Rayleigh waves, I evaluate the influence of surface topography to conventional dispersion analysis in 2D and 3D domains by numerical investigations. At last I examine the feasibility of inverting waveforms of Rayleigh waves for shallow S-wave velocities using a genetic algorithm. Results of the study show that Rayleigh waves can be accurately simulated in near surface using the improved vacuum formulation. Spurious reflections during the numerical modeling can be efficiently suppressed by the simplified multiaxial perfectly matched layers. The
Locating the Origin of Scattered Waves By Simulating Time Reversal of the Seismic Wavefield
NASA Astrophysics Data System (ADS)
Myers, S. C.; Pitarka, A.; Sjogreen, B.; Petersson, A.; Simmons, N. A.; Johannesson, G.
2014-12-01
The Source Physics Experiment (SPE) is a series of underground chemical explosions at the Nevada National Security Site (NNSS) that are improving our physical understanding how explosion sources generate seismic waves. Better understanding the origin of S-waves from explosions is a primary goal of the SPE. Even at distances of a few kilometers from the SPE sources, seismic recordings include arrivals of unknown origin that could originate as S-waves at the explosive source or from topographic and subsurface scatterers. Back propagation of time reversed seismograms has been used to determine the location of seismic events (e.g. Tromp et al., 2005; Larmat et al., 2006), and Myers et al. (2007) demonstrated that the time-reversal method can be used to determine the origin of direct and scattered waves in seismic simulations. In this study we identify the origin of distinct features in synthetic seismograms that are generated by elastic, finite-difference simulation of seismic propagation from SPE explosions through a model that has been developed specifically for the SPE. The SPE model includes 3-dimensional velocity discontinuities at geologic boundaries, as well as free-surface topography. Although the largest arrivals in the synthetic seismograms are expected to originate at the explosion source, other prominent features are likely to originate as scattered energy from model discontinuities. Scattering sources in the SPE model that are needed in order to match synthetic seismograms to field recordings of SPE shots will be identified. Conversely, model structures may be removed if they result in disagreement between synthetic seismograms and field recordings. Ultimately, we plan to constrain the origin of prominent features in field recordings of SPE shots by directly using the field recordings as inputs to time reversal simulations. Direct use of field recordings will require development of methods that account for the uncertainty of the seismic model through which
Sub-wavelength energy trapping of elastic waves in a metamaterial.
Colombi, Andrea; Roux, Philippe; Rupin, Matthieu
2014-08-01
Deep sub-wavelength focusing has been demonstrated for locally resonant metamaterials using electromagnetic and acoustic waves. The elastic equivalents of such objects are made of sub-wavelength resonating beams fixed to a two-dimensional plate, as presented here. Independent of a random or regular arrangement of the resonators, the metamaterial shows large bandgaps that are independent of the incident wave direction. Numerical simulations demonstrate that the insertion of a defect in the layout, as a shorter resonator, creates strong amplification of the wave-field on the defect. This energy trapping, which is localized on a spatial scale that is much smaller than the wavelength in the two-dimensional plate, leads to a >1 factor in terms of the local density of energy. PMID:25096146
Sun, Qicheng; Jin, Feng; Wang, Guangqian; Song, Shixiong; Zhang, Guohua
2015-01-01
Mesoscopic structures form in dense granular materials due to the self-organisation of the constituent particles. These structures have internal structural degrees of freedom in addition to the translational degree of freedom. The resultant granular elasticity, which exhibits intrinsic variations and inevitable relaxation, is a key quantity that accounts for macroscopic solid- or fluid-like properties and the transitions between them. In this work, we propose a potential energy landscape (PEL) with local stable basins and low elastic energy barriers to analyse the nature of granular elasticity. A function for the elastic energy density is proposed for stable states and is further calibrated with ultrasonic measurements. Fluctuations in the elastic energy due to the evolution of internal structures are proposed to describe a so-called configuration temperature Tc as a counterpart of the classical kinetic granular temperature Tk that is attributed to the translational degrees of freedom. The two granular temperatures are chosen as the state variables, and a fundamental equation is established to develop non-equilibrium thermodynamics for granular materials. Due to the relatively low elastic energy barrier in the PEL, granular elasticity relaxes more under common mechanical loadings, and a simple model based on mean-field theory is developed to account for this behaviour. PMID:25951049
Elastic membranes in confinement
NASA Astrophysics Data System (ADS)
Bostwick, Joshua; Miksis, Michael; Davis, Stephen
2014-11-01
An elastic membrane stretched between two walls takes a shape defined by its length and the volume of fluid it encloses. Many biological structures, such as cells, mitochondria and DNA, have finer internal structure in which a membrane (or elastic member) is geometrically ``confined'' by another object. We study the shape stability of elastic membranes in a ``confining'' box and introduce repulsive van der Waals forces to prevent the membrane from intersecting the wall. We aim to define the parameter space associated with mitochondria-like deformations. We compare the confined to `unconfined' solutions and show how the structure and stability of the membrane shapes changes with the system parameters.
Mechanism of Resilin Elasticity
Qin, Guokui; Hu, Xiao; Cebe, Peggy; Kaplan, David L.
2012-01-01
Resilin is critical in the flight and jumping systems of insects as a polymeric rubber-like protein with outstanding elasticity. However, insight into the underlying molecular mechanisms responsible for resilin elasticity remains undefined. Here we report the structure and function of resilin from Drosophila CG15920. A reversible beta-turn transition was identified in the peptide encoded by exon III and for full length resilin during energy input and release, features that correlate to the rapid deformation of resilin during functions in vivo. Micellar structures and nano-porous patterns formed after beta-turn structures were present via changes in either the thermal or mechanical inputs. A model is proposed to explain the super elasticity and energy conversion mechanisms of resilin, providing important insight into structure-function relationships for this protein. Further, this model offers a view of elastomeric proteins in general where beta-turn related structures serve as fundamental units of the structure and elasticity. PMID:22893127
NASA Astrophysics Data System (ADS)
Thomson, C. J.
2005-10-01
Several observations are made concerning the numerical implementation of wide-angle one-way wave equations, using for illustration scalar waves obeying the Helmholtz equation in two space dimensions. This simple case permits clear identification of a sequence of physically motivated approximations of use when the mathematically exact pseudo-differential operator (PSDO) one-way method is applied. As intuition suggests, these approximations largely depend on the medium gradients in the direction transverse to the main propagation direction. A key point is that narrow-angle approximations are to be avoided in the interests of accuracy. Another key consideration stems from the fact that the so-called `standard-ordering' PSDO indicates how lateral interpolation of the velocity structure can significantly reduce computational costs associated with the Fourier or plane-wave synthesis lying at the heart of the calculations. A third important point is that the PSDO theory shows what approximations are necessary in order to generate an exponential one-way propagator for the laterally varying case, representing the intuitive extension of classical integral-transform solutions for a laterally homogeneous medium. This exponential propagator permits larger forward stepsizes. Numerical comparisons with Helmholtz (i.e. full) wave-equation finite-difference solutions are presented for various canonical problems. These include propagation along an interfacial gradient, the effects of a compact inclusion and the formation of extended transmitted and backscattered wave trains by model roughness. The ideas extend to the 3-D, generally anisotropic case and to multiple scattering by invariant embedding. It is concluded that the method is very competitive, striking a new balance between simplifying approximations and computational labour. Complicated wave-scattering effects are retained without the need for expensive global solutions, providing a robust and flexible modelling tool.
NASA Astrophysics Data System (ADS)
Quilliet, Catherine; Quemeneur, François; Marmottant, Philippe; Imhof, Arnout; Pépin-Donat, Brigitte; van Blaaderen, Alfons
2010-03-01
The deflation of elastic spherical surfaces has been numerically investigated, and show very different types of deformations according the range of elastic parameters, some of them being quantitatively explained through simple calculations. This allows to retrieve various shapes observed on hollow shells (from colloidal to centimeter scale), on lipid vesicles, or on some biological objects. The extension of this process to other geometries allows to modelize vegetal objects such as the ultrafast trap of carnivorous plants.
NASA Astrophysics Data System (ADS)
Khaniani, Hassan
boundary condition of the wave equation is set up along reflection surfaces. Hence, the surface integral Kirchhoff approximation is used as a mathematical framework instead of the volume integral of the Born approximation. In addition, I study the feasibility of iterative coupling of ray theory with the Kirchhoff approximation for inversion. For the amplitude considerations, the direct relationship between the scattering potential of the Born approximation with the reflectivity function of the asymptotic Kirchhoff approximation for elastic waves is used. Therefore, I use the linearized Zoeppritz approximation of Aki and Richards (1980) for computation of the forward modeling and migration operators as well as gradient function from Amplitude vs Offset (AVO) inversion. The multiparameter elastic inversion approach is applicable to all types of reflected wavefields such as P-to-P, P-to-S, S-to-S and S-to-P. Traveltime estimation of forward modeling and migration/inversion operators are based on the DSR equation. All operators involved in inversion, including the background model for DSR and AVO are updated at each iteration. The migration/inversion procedure maps the mode converted waves to the traveltime of incident waves which fixes the registration problem of events that travel from source to scatter point. The inversion of the reflected P-to-P and P-to-S synthetic and field data are provided for the numerical examples. This approach is applicable for complex structures however, to estimate the traveltime of scatterpoints, ray tracing can be added to the algorithm. For such a medium, the scatterpoint traveltime approximations from the PSTM, is compared to the PSDM approach using numerical analysis of ray- and FDTD-based modeling. In part of this thesis, I further improve the conventional velocity analysis of Common Scatter Point (CSP) gathers by including the tilt effects. I show that travel time response of scatter points beneath a dipping interface experiences an
Elastic Collisions and Gravity
NASA Astrophysics Data System (ADS)
Ball, Steven
2009-04-01
Elastic collisions are fascinating demonstrations of conservation principles. The mediating force must be conservative in an elastic collision. Truly elastic collisions take place only when the objects in collision do not touch, e.g. magnetic bumpers on low friction carts. This requires that we define a collision as a momentum transfer. Elastic collisions in 1-D can be solved in general and the implications are quite remarkable. For example, a heavy object moving initially towards a light object followed by an elastic collision results in a final velocity of the light object greater than either initial velocity. This is easily demonstrated with low friction carts. Gravitational elastic collisions involving a light spacecraft and an extremely massive body like a moon or planet can be approximated as 1-D collisions, such as the ``free return'' trajectory of Apollo 13 around the moon. The most fascinating gravitational collisions involve the gravitational slingshot effect used to boost spacecraft velocities. The maximum gravitational slingshot effect occurs when approaching a nearly 1-D collision, revealing that the spacecraft can be boosted to greater than twice the planet velocity, enabling the spacecraft to travel much further away from the Sun.
NASA Astrophysics Data System (ADS)
Chen, H.; Li, K.
2012-12-01
We applied a wave-equation based adjoint wavefield method for seismic illumination/resolution analyses and full waveform inversion. A two-way wave-equation is used to calculate directional and diffracted energy fluxes for waves propagating between sources and receivers to the subsurface target. The first-order staggered-grid pressure-velocity formulation, which lacks the characteristic of being self-adjoint is further validated and corrected to render the modeling operator before its practical application. Despite most published papers on synthetic kernel research, realistic applications to two field experiments are demonstrated and emphasize its practical needs. The Fréchet sensitivity kernels are used to quantify the target illumination conditions. For realistic illumination measurements and resolution analyses, two completely different survey geometries and nontrivial pre-conditioning strategies based on seismic data type are demonstrated and compared. From illumination studies, particle velocity responses are more sensitive to lateral velocity variations than pressure records. For waveform inversion, the more accurately estimated velocity model obtained the deeper the depth of investigation would be reached. To achieve better resolution and illumination, closely spaced OBS receiver interval is preferred. Based on the results, waveform inversion is applied for a gas hydrate site in Taiwan for shallow structure and BSR detection. Full waveform approach potentially provides better depth resolution than ray approach. The quantitative analyses, a by-product of full waveform inversion, are useful for quantifying seismic processing and depth migration strategies.llumination/resolution analysis for a 3D MCS/OBS survey in 2008. Analysis of OBS data shows that pressure (top), horizontal (middle) and vertical (bottom) velocity records produce different resolving power for gas hydrate exploration. ull waveform inversion of 8 OBS data along Yuan-An Ridge in SW Taiwan
Elasticity of plagioclase feldspars
NASA Astrophysics Data System (ADS)
Brown, J. Michael; Angel, Ross J.; Ross, Nancy L.
2016-02-01
Elastic properties are reported for eight plagioclase feldspars that span compositions from albite (NaSi3AlO8) to anorthite (CaSi2Al2O8). Surface acoustic wave velocities measured using Impulsive Stimulated Light Scattering and compliance sums from high-pressure X-ray compression studies accurately determine all 21 components of the elasticity tensor for these triclinic minerals. The overall pattern of elasticity and the changes in individual elastic components with composition can be rationalized on the basis of the evolution of crystal structures and chemistry across this solid-solution join. All plagioclase feldspars have high elastic anisotropy; a* (the direction perpendicular to the b and c axes) is the softest direction by a factor of 3 in albite. From albite to anorthite the stiffness of this direction undergoes the greatest change, increasing twofold. Small discontinuities in the elastic components, inferred to occur between the three plagioclase phases with distinct symmetry (C1>¯, I1>¯, and P1>¯), appear consistent with the nature of the underlying conformation of the framework-linked tetrahedra and the associated structural changes. Measured body wave velocities of plagioclase-rich rocks, reported over the last five decades, are consistent with calculated Hill-averaged velocities using the current moduli. This confirms long-standing speculation that previously reported elastic moduli for plagioclase feldspars are systematically in error. The current results provide greater assurance that the seismic structure of the middle and lower crusts can be accurately estimated on the basis of specified mineral modes, chemistry, and fabric.
NASA Astrophysics Data System (ADS)
Gu, Bingluo; Li, Zhiyuan; Ma, Xiaona; Liang, Guanghe
2015-01-01
The vector P- and S-seismograms in an elastic common-shot gather generated with a P-wave source in a two-dimensional model using a finite difference solution of the P- and S-wave separated velocity-stress equations can be imaged by two independent ERTMs based on the same equations. The inputs as boundary conditions for reverse-time extrapolation are the pure vector P- and S-waves, respectively. The vector P-wave image components can be obtained by the normalized correlation operation between the vector P-source wavefields and receiver wavefields, which are obtained by extrapolating the vector P-seismograms in reverse-time using the same equations as the forward modeling. The vector S-wave image components can be obtained by a similar method. Compared with the conventional ERTM, this method can minimize the artifacts caused by the crosstalk between different wave modes and can preserve the phase and amplitude attributes of migration images very well. Furthermore, the polarity-reversal of the vector S-wave data can be corrected automatically during the imaging process, so destructive interferences between data from adjacent sources do not exist. Numerical examples with synthetic data have demonstrated the feasibility and effectiveness of this method for complex structures.
NASA Astrophysics Data System (ADS)
Zacharias, Mario; Paul, Indranil; Garst, Markus
2015-07-01
We discuss elastic instabilities of the atomic crystal lattice at zero temperature. Because of long-range shear forces of the solid, at such transitions the phonon velocities vanish, if at all, only along certain crystallographic directions, and, consequently, the critical phonon fluctuations are suppressed to a lower dimensional manifold and governed by a Gaussian fixed point. In the case of symmetry-breaking elastic transitions, a characteristic critical phonon thermodynamics arises that is found, e.g., to violate Debye's T3 law for the specific heat. We point out that quantum critical elasticity is triggered whenever a critical soft mode couples linearly to the strain tensor. In particular, this is relevant for the electronic Ising-nematic quantum phase transition in a tetragonal crystal as discussed in the context of certain cuprates, ruthenates, and iron-based superconductors.
Norris, Andrew N.
2014-01-01
We consider a periodic lattice structure in d=2 or 3 dimensions with unit cell comprising Z thin elastic members emanating from a similarly situated central node. A general theoretical approach provides an algebraic formula for the effective elasticity of such frameworks. The method yields the effective cubic elastic constants for three-dimensional space-filling lattices with Z=4, 6, 8, 12 and 14, the last being the ‘stiffest’ lattice proposed by Gurtner & Durand (Gurtner & Durand 2014 Proc. R. Soc. A 470, 20130611. (doi:10.1098/rspa.2013.0611)). The analytical expressions provide explicit formulae for the effective properties of pentamode materials, both isotropic and anisotropic, obtained from the general formulation in the stretch-dominated limit for Z=d+1. PMID:25484608
Yu, Betty; Kang, Soo-Young; Akthakul, Ariya; Ramadurai, Nithin; Pilkenton, Morgan; Patel, Alpesh; Nashat, Amir; Anderson, Daniel G; Sakamoto, Fernanda H; Gilchrest, Barbara A; Anderson, R Rox; Langer, Robert
2016-08-01
We report the synthesis and application of an elastic, wearable crosslinked polymer layer (XPL) that mimics the properties of normal, youthful skin. XPL is made of a tunable polysiloxane-based material that can be engineered with specific elasticity, contractility, adhesion, tensile strength and occlusivity. XPL can be topically applied, rapidly curing at the skin interface without the need for heat- or light-mediated activation. In a pilot human study, we examined the performance of a prototype XPL that has a tensile modulus matching normal skin responses at low strain (<40%), and that withstands elongations exceeding 250%, elastically recoiling with minimal strain-energy loss on repeated deformation. The application of XPL to the herniated lower eyelid fat pads of 12 subjects resulted in an average 2-grade decrease in herniation appearance in a 5-point severity scale. The XPL platform may offer advanced solutions to compromised skin barrier function, pharmaceutical delivery and wound dressings. PMID:27159017
Peselnick, L.; Robie, R.A.
1962-01-01
The recent measurements of the elastic constants of calcite by Reddy and Subrahmanyam (1960) disagree with the values obtained independently by Voigt (1910) and Bhimasenachar (1945). The present authors, using an ultrasonic pulse technique at 3 Mc and 25??C, determined the elastic constants of calcite using the exact equations governing the wave velocities in the single crystal. The results are C11=13.7, C33=8.11, C44=3.50, C12=4.82, C13=5.68, and C14=-2.00, in units of 1011 dyncm2. Independent checks of several of the elastic constants were made employing other directions and polarizations of the wave velocities. With the exception of C13, these values substantially agree with the data of Voigt and Bhimasenachar. ?? 1962 The American Institute of Physics.
Zacharias, Mario; Paul, Indranil; Garst, Markus
2015-07-10
We discuss elastic instabilities of the atomic crystal lattice at zero temperature. Because of long-range shear forces of the solid, at such transitions the phonon velocities vanish, if at all, only along certain crystallographic directions, and, consequently, the critical phonon fluctuations are suppressed to a lower dimensional manifold and governed by a Gaussian fixed point. In the case of symmetry-breaking elastic transitions, a characteristic critical phonon thermodynamics arises that is found, e.g., to violate Debye's T(3) law for the specific heat. We point out that quantum critical elasticity is triggered whenever a critical soft mode couples linearly to the strain tensor. In particular, this is relevant for the electronic Ising-nematic quantum phase transition in a tetragonal crystal as discussed in the context of certain cuprates, ruthenates, and iron-based superconductors. PMID:26207483
NASA Astrophysics Data System (ADS)
Yu, Betty; Kang, Soo-Young; Akthakul, Ariya; Ramadurai, Nithin; Pilkenton, Morgan; Patel, Alpesh; Nashat, Amir; Anderson, Daniel G.; Sakamoto, Fernanda H.; Gilchrest, Barbara A.; Anderson, R. Rox; Langer, Robert
2016-08-01
We report the synthesis and application of an elastic, wearable crosslinked polymer layer (XPL) that mimics the properties of normal, youthful skin. XPL is made of a tunable polysiloxane-based material that can be engineered with specific elasticity, contractility, adhesion, tensile strength and occlusivity. XPL can be topically applied, rapidly curing at the skin interface without the need for heat- or light-mediated activation. In a pilot human study, we examined the performance of a prototype XPL that has a tensile modulus matching normal skin responses at low strain (<40%), and that withstands elongations exceeding 250%, elastically recoiling with minimal strain-energy loss on repeated deformation. The application of XPL to the herniated lower eyelid fat pads of 12 subjects resulted in an average 2-grade decrease in herniation appearance in a 5-point severity scale. The XPL platform may offer advanced solutions to compromised skin barrier function, pharmaceutical delivery and wound dressings.
NASA Astrophysics Data System (ADS)
Habibi, M.; Ribe, N. M.; Bonn, Daniel
2007-10-01
A rope falling onto a solid surface typically forms a series of regular coils. Here, we study this phenomenon using laboratory experiments (with cotton threads and softened spaghetti) and an asymptotic “slender-rope” numerical model. The excellent agreement between the two with no adjustable parameters allows us to determine a complete phase diagram for elastic coiling comprising three basic regimes involving different force balances (elastic, gravitational, and inertial) together with resonant “whirling string” and “whirling shaft” eigenmodes in the inertial regime.
Dynamic Elasticity Model of Resilin Biopolymers
NASA Astrophysics Data System (ADS)
Hu, Xiao; Duki, Solomon
2013-03-01
Resilin proteins are `super elastic rubbers' in the flight and jumping systems of most insects, and can extend and retract millions of times. Natural resilin exhibits high resilience (> 95%) under high-frequency conditions, and could be stretched to over 300% of its original length with a low elastic modulus of 0.1-3 MPa. However, insight into the underlying molecular mechanisms responsible for resilin elasticity remains undefined. We report on the dynamic structure transitions and functions of full length resilin from fruit fly (D. melanogaster CG15920) and its different functional domains. A dynamic computational model is proposed to explain the super elasticity and energy conversion mechanisms of resilin, providing important insight into structure-function relationships for resilins, as well as other elastomeric proteins. A strong beta-turn transition was experimentally identified in the full length resilin and its non-elastic domains (Exon III). Changes in periodic long-range order were demonstrated during this transition, induced either by thermal or mechanical inputs, to confirm the universality of proposed mechanism. Further, this model offers new options for designing protein-based biopolymers with tunable material applications.
Vascular elastic photoacoustic tomography in humans
NASA Astrophysics Data System (ADS)
Hai, Pengfei; Zhou, Yong; Liang, Jinyang; Li, Chiye; Wang, Lihong V.
2016-03-01
Quantification of vascular elasticity can help detect thrombosis and prevent life-threatening conditions such as acute myocardial infarction or stroke. Here, we propose vascular elastic photoacoustic tomography (VE-PAT) to measure vascular elasticity in humans. VE-PAT was developed by incorporating a linear-array-based photoacoustic computed tomography system with a customized compression stage. By measuring the deformation of blood vessels under uniaxial loading, VE-PAT was able to quantify the vascular compliance. We first demonstrated the feasibility of VE-PAT in blood vessel phantoms. In large vessel phantoms, VE-PAT detected a decrease in vascular compliance due to simulated thrombosis, which was validated by a standard compression test. In small blood vessel phantoms embedded 3 mm deep in gelatin, VE-PAT detected elasticity changes at depths that are difficult to image using other elasticity imaging techniques. We then applied VE-PAT to assess vascular compliance in a human subject and detected a decrease in vascular compliance when an occlusion occurred downstream from the measurement point, demonstrating the potential of VE-PAT in clinical applications such as detection of deep venous thrombosis.
NASA Astrophysics Data System (ADS)
Pavlis, G. L.; Wang, Y.
2015-12-01
A significant number of P and S wave tomography models have been produced in the past decade using various subsets of data from the Earthscope USArray and different inversion algorithms. We focus here on published tomography results that span large portions of the final footprint of the USArray. We use 3D visualization techniques to search for common features in different tomography models. We also compare tomography results to features seen in our current generation wavefield images. Recent innovations of our plane wave migration method have yielded what is arguably the highest resolution image ever produced of the mantle in the vicinity of the transition zone. The new results reveal a rich collection of coherent, dipping structures seen throughout the upper mantle and transition zone. These dipping interfaces are judged significant according to a coherence metric. We treat these surfaces as strain markers to assess proposed models for geometry of the 3D geometry of the Farallon Slab under North America. We find the following geologic interpretations are well supported by independent results: 1. The old Farallon under eastern North America and below the base of transition zone is universally seen as a high velocity anomaly. 2. All results support a simple, 3D kinematic model of the updip limit of the Farallon slab window that follows a track from Cape Mendocino, across Nevada, and northern Arizona and New Mexico. 3. All models show a strong low-velocity mantle under the southwestern U.S. 4. A low-velocity features is universally seen related to the Yellowstone-Snake River system. Shorter wavelength features observed in different tomography models are inconsistent showing that the theme of this session is very important to understand what features are in current results are real. Isopach maps of the thickness of the transition show a systematic difference in transition zone thickness in the western and eastern US. The transition zone thickens in the eastern US in
McHugh, Stuart
1976-01-01
The material in this report can be grouped into two categories: 1) programs that compute tilts produced by a vertically oriented expanding rectangular dislocation loop in an elastic or viscoelastic material and 2) programs that compute the shear stresses, strains, and shear displacements in a three-phase half-space (i.e. a half-space containing a vertical slab). Each section describes the relevant theory, and provides a detailed guide to the operation of the programs. A series of examples is provided at the end of each section.
Elastic and Inelastic Collisions
ERIC Educational Resources Information Center
Gluck, Paul
2010-01-01
There have been two articles in this journal that described a pair of collision carts used to demonstrate vividly the difference between elastic and inelastic collisions. One cart had a series of washers that were mounted rigidly on a rigid wooden framework, the other had washers mounted on rubber bands stretched across a framework. The rigidly…
ERIC Educational Resources Information Center
Gordon, Warren B.
2006-01-01
This paper examines the elasticity of demand, and shows that geometrically, it may be interpreted as the ratio of two simple distances along the tangent line: the distance from the point on the curve to the x-intercept to the distance from the point on the curve to the y-intercept. It also shows that total revenue is maximized at the transition…
ERIC Educational Resources Information Center
Cocco, Alberto; Masin, Sergio Cesare
2010-01-01
Participants estimated the imagined elongation of a spring while they were imagining that a load was stretching the spring. This elongation turned out to be a multiplicative function of spring length and load weight--a cognitive law analogous to Hooke's law of elasticity. Participants also estimated the total imagined elongation of springs joined…
Elastically tailored composite structures
NASA Technical Reports Server (NTRS)
2000-01-01
Elastically tailored composite structures using out-of-autoclave processes. Several unsymetric autoclave-cured and electron-beam-cured composite laminates are compared. Cantilevered beam (unbalanced/asymetric laminate) used to demonstrate bend-twist coupling effects. Photographed in building 1145, photographic studio.
Hydrodynamic Elastic Magneto Plastic
Energy Science and Technology Software Center (ESTSC)
1985-02-01
The HEMP code solves the conservation equations of two-dimensional elastic-plastic flow, in plane x-y coordinates or in cylindrical symmetry around the x-axis. Provisions for calculation of fixed boundaries, free surfaces, pistons, and boundary slide planes have been included, along with other special conditions.
Renormalization of curvature elastic constants for elastic and fluid membranes
NASA Astrophysics Data System (ADS)
Ami, S.; Kleinert, H.
1987-02-01
We study the fluctuations of membranes with area and curvature elasticity and calculate the renormalization of the curvature elastic constants due to thermal fluctuations. For the mean curvature elastic constant the result is the same as obtained previously for “ideal membranes” which resist only to curvature deformations. The renormalization of the gaussian curvature, on the other hand, depends on the elastic contants. In an incompressible membrane, it is five times weaker than in an ideal membrane.
Elastic wave turbulence and intermittency
NASA Astrophysics Data System (ADS)
Chibbaro, Sergio; Josserand, Christophe
2016-07-01
We investigate the onset of intermittency for vibrating elastic plate turbulence in the framework of the weak wave turbulence theory using a numerical approach. The spectrum of the displacement field and the structure functions of the fluctuations are computed for different forcing amplitudes. At low forcing, the spectrum predicted by the theory is observed, while the fluctuations are consistent with Gaussian statistics. When the forcing is increased, the spectrum varies at large scales, corresponding to the oscillations of nonlinear structures made of ridges delimited by d cones. In this regime, the fluctuations exhibit small-scale intermittency that can be fitted via a multifractal model. The analysis of the nonlinear frequency shows that the intermittency is linked to the breakdown of the weak turbulence at large scales only.
Elastic wave turbulence and intermittency.
Chibbaro, Sergio; Josserand, Christophe
2016-07-01
We investigate the onset of intermittency for vibrating elastic plate turbulence in the framework of the weak wave turbulence theory using a numerical approach. The spectrum of the displacement field and the structure functions of the fluctuations are computed for different forcing amplitudes. At low forcing, the spectrum predicted by the theory is observed, while the fluctuations are consistent with Gaussian statistics. When the forcing is increased, the spectrum varies at large scales, corresponding to the oscillations of nonlinear structures made of ridges delimited by d cones. In this regime, the fluctuations exhibit small-scale intermittency that can be fitted via a multifractal model. The analysis of the nonlinear frequency shows that the intermittency is linked to the breakdown of the weak turbulence at large scales only. PMID:27575068
NASA Astrophysics Data System (ADS)
Qi, Yusheng; Yue, Dick
2015-11-01
We use direct nonlinear phase-resolved simulations based on a High-Order Spectral (HOS) method (Dommermuth & Yue 1987) to understand and quantify wave-breaking dissipation in the evolution of general irregular short-crested wave-fields. We achieve this by incorporating a robust phenomenological-based wave breaking model in HOS simulations to account for energy dissipation. This model can automatically simulate the onset of wave breaking, and the simulated wave-breaking dissipation strength differentiates corresponding to different wave breaking type (such as spilling or plunging breaking waves). The efficacy of this model is confirmed by direct comparisons against measurements for the energy loss in 2D and 3D breaking events. By comparing simulated wave-fields with and without the dissipation model in HOS, we obtain the dissipation field, which provides the times, locations and intensity of wave breaking events. From the dissipation field we further calculate the distribution of total length of breaking wave front per unit surface area per unit increment of breaking velocity (Phillips 1985), and obtain qualitative agreement with Phillips theoretical power-law.
Frequency dependent elastic impedance inversion for interstratified dispersive elastic parameters
NASA Astrophysics Data System (ADS)
Zong, Zhaoyun; Yin, Xingyao; Wu, Guochen
2016-08-01
The elastic impedance equation is extended to frequency dependent elastic impedance equation by taking partial derivative to frequency. With this equation as the forward solver, a practical frequency dependent elastic impedance inversion approach is presented to implement the estimation of the interstratified dispersive elastic parameters which makes full use of the frequency information of elastic impedances. Three main steps are included in this approach. Firstly, the elastic Bayesian inversion is implemented for the estimation of elastic impedances from different incident angle. Secondly, with those estimated elastic impedances, their variations are used to estimate P-wave velocity and S-wave velocity. Finally, with the prior elastic impedance and P-wave and S-wave velocity information, the frequency dependent elastic variation with incident angle inversion is presented for the estimation of the interstratified elastic parameters. With this approach, the interstratified elastic parameters rather than the interface information can be estimated, making easier the interpretation of frequency dependent seismic attributes. The model examples illustrate the feasibility and stability of the proposed method in P-wave velocity dispersion and S-wave velocity dispersion estimation. The field data example validates the possibility and efficiency in hydrocarbon indication of the estimated P-wave velocity dispersion and S-wave velocity dispersion.
Elastic properties of hedenbergite
NASA Astrophysics Data System (ADS)
Kandelin, John; Weidner, Donald J.
1988-02-01
The single-crystal elastic moduli of hedenbergite (CaFeSi2O6) hare been measured at 20°C and 1 bar using Brillouin spectroscopy. The moduli are (in gigapascals): C11 = 222, C22 = 176, C23 = 249, C44 = 55, C55 = 63, C66 = 60, C12 = 69, C13 = 79, C33, = 86, C15 = 12, C25 = 13, C35 = 26, C46 = -10. The comparison of elastic properties among Mg-Fe-Ca bearing pyroxenes, known as quadrilateral pyroxenes, reveals only weak variations with changes in composition. Of the four quadrilateral pyroxenes, orthoferrosilite has elastic properties distinctive from the others. The principal differences among these pyroxenes are due to subtle structural differences. In particular, the mechanical linkage between the M2 polyhedral chains in clinopyroxenes enhances the importance of the cation in this site. In contrast to the orthopyroxenes, the aggregate shear modulus μ of the calcium-bearing clinopyroxenes (diopside and hedenbergite) exhibits no dependence on the amount of iron (Fe2+) present in the structure, while the ratio K/μ does. As a result, the compressional and shear acoustic velocities of the calcium-bearing clinopyroxenes show a smaller dependency on iron content than do the orthopyroxenes.
Elastic properties of gamma-Pu by resonant ultrasound spectroscopy
Migliori, Albert; Betts, J; Trugman, A; Mielke, C H; Mitchell, J N; Ramos, M; Stroe, I
2009-01-01
Despite intense experimental and theoretical work on Pu, there is still little understanding of the strange properties of this metal. We used resonant ultrasound spectroscopy method to investigate the elastic properties of pure polycrystalline Pu at high temperatures. Shear and longitudinal elastic moduli of the {gamma}-phase of Pu were determined simultaneously and the bulk modulus was computed from them. A smooth linear and large decrease of all elastic moduli with increasing temperature was observed. We calculated the Poisson ratio and found that it increases from 0.242 at 519K to 0.252 at 571K.
Elastic regimes of subisostatic athermal fiber networks
NASA Astrophysics Data System (ADS)
Licup, A. J.; Sharma, A.; MacKintosh, F. C.
2016-01-01
Athermal models of disordered fibrous networks are highly useful for studying the mechanics of elastic networks composed of stiff biopolymers. The underlying network architecture is a key aspect that can affect the elastic properties of these systems, which include rich linear and nonlinear elasticity. Existing computational approaches have focused on both lattice-based and off-lattice networks obtained from the random placement of rods. It is not obvious, a priori, whether the two architectures have fundamentally similar or different mechanics. If they are different, it is not clear which of these represents a better model for biological networks. Here, we show that both approaches are essentially equivalent for the same network connectivity, provided the networks are subisostatic with respect to central force interactions. Moreover, for a given subisostatic connectivity, we even find that lattice-based networks in both two and three dimensions exhibit nearly identical nonlinear elastic response. We provide a description of the linear mechanics for both architectures in terms of a scaling function. We also show that the nonlinear regime is dominated by fiber bending and that stiffening originates from the stabilization of subisostatic networks by stress. We propose a generalized relation for this regime in terms of the self-generated normal stresses that develop under deformation. Different network architectures have different susceptibilities to the normal stress but essentially exhibit the same nonlinear mechanics. Such a stiffening mechanism has been shown to successfully capture the nonlinear mechanics of collagen networks.