3D Elastic Seismic Wave Propagation Code
Energy Science and Technology Software Center (ESTSC)
1998-09-23
E3D is capable of simulating seismic wave propagation in a 3D heterogeneous earth. Seismic waves are initiated by earthquake, explosive, and/or other sources. These waves propagate through a 3D geologic model, and are simulated as synthetic seismograms or other graphical output.
Spatial parallelism of a 3D finite difference, velocity-stress elastic wave propagation code
Minkoff, S.E.
1999-12-01
Finite difference methods for solving the wave equation more accurately capture the physics of waves propagating through the earth than asymptotic solution methods. Unfortunately, finite difference simulations for 3D elastic wave propagation are expensive. The authors model waves in a 3D isotropic elastic earth. The wave equation solution consists of three velocity components and six stresses. The partial derivatives are discretized using 2nd-order in time and 4th-order in space staggered finite difference operators. Staggered schemes allow one to obtain additional accuracy (via centered finite differences) without requiring additional storage. The serial code is most unique in its ability to model a number of different types of seismic sources. The parallel implementation uses the MPI library, thus allowing for portability between platforms. Spatial parallelism provides a highly efficient strategy for parallelizing finite difference simulations. In this implementation, one can decompose the global problem domain into one-, two-, and three-dimensional processor decompositions with 3D decompositions generally producing the best parallel speedup. Because I/O is handled largely outside of the time-step loop (the most expensive part of the simulation) the authors have opted for straight-forward broadcast and reduce operations to handle I/O. The majority of the communication in the code consists of passing subdomain face information to neighboring processors for use as ghost cells. When this communication is balanced against computation by allocating subdomains of reasonable size, they observe excellent scaled speedup. Allocating subdomains of size 25 x 25 x 25 on each node, they achieve efficiencies of 94% on 128 processors. Numerical examples for both a layered earth model and a homogeneous medium with a high-velocity blocky inclusion illustrate the accuracy of the parallel code.
Spatial Parallelism of a 3D Finite Difference, Velocity-Stress Elastic Wave Propagation Code
MINKOFF,SUSAN E.
1999-12-09
Finite difference methods for solving the wave equation more accurately capture the physics of waves propagating through the earth than asymptotic solution methods. Unfortunately. finite difference simulations for 3D elastic wave propagation are expensive. We model waves in a 3D isotropic elastic earth. The wave equation solution consists of three velocity components and six stresses. The partial derivatives are discretized using 2nd-order in time and 4th-order in space staggered finite difference operators. Staggered schemes allow one to obtain additional accuracy (via centered finite differences) without requiring additional storage. The serial code is most unique in its ability to model a number of different types of seismic sources. The parallel implementation uses the MP1 library, thus allowing for portability between platforms. Spatial parallelism provides a highly efficient strategy for parallelizing finite difference simulations. In this implementation, one can decompose the global problem domain into one-, two-, and three-dimensional processor decompositions with 3D decompositions generally producing the best parallel speed up. Because i/o is handled largely outside of the time-step loop (the most expensive part of the simulation) we have opted for straight-forward broadcast and reduce operations to handle i/o. The majority of the communication in the code consists of passing subdomain face information to neighboring processors for use as ''ghost cells''. When this communication is balanced against computation by allocating subdomains of reasonable size, we observe excellent scaled speed up. Allocating subdomains of size 25 x 25 x 25 on each node, we achieve efficiencies of 94% on 128 processors. Numerical examples for both a layered earth model and a homogeneous medium with a high-velocity blocky inclusion illustrate the accuracy of the parallel code.
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
NASA Astrophysics Data System (ADS)
Wang, S.; De Hoop, M. V.; Xia, J.; Li, X.
2011-12-01
We consider the modeling of elastic seismic wave propagation on a rectangular domain via the discretization and solution of the inhomogeneous coupled Helmholtz equation in 3D, by exploiting a parallel multifrontal sparse direct solver equipped with Hierarchically Semi-Separable (HSS) structure to reduce the computational complexity and storage. In particular, we are concerned with solving this equation on a large domain, for a large number of different forcing terms in the context of seismic problems in general, and modeling in particular. We resort to a parsimonious mixed grid finite differences scheme for discretizing the Helmholtz operator and Perfect Matched Layer boundaries, resulting in a non-Hermitian matrix. We make use of a nested dissection based domain decomposition, and introduce an approximate direct solver by developing a parallel HSS matrix compression, factorization, and solution approach. We cast our massive parallelization in the framework of the multifrontal method. The assembly tree is partitioned into local trees and a global tree. The local trees are eliminated independently in each processor, while the global tree is eliminated through massive communication. The solver for the inhomogeneous equation is a parallel hybrid between multifrontal and HSS structure. The computational complexity associated with the factorization is almost linear with the size of the Helmholtz matrix. Our numerical approach can be compared with the spectral element method in 3D seismic applications.
NASA Astrophysics Data System (ADS)
Petrov, P.; Newman, G. A.
2010-12-01
-Fourier domain we had developed 3D code for full-wave field simulation in the elastic media which take into account nonlinearity introduced by free-surface effects. Our approach is based on the velocity-stress formulation. In the contrast to conventional formulation we defined the material properties such as density and Lame constants not at nodal points but within cells. This second order finite differences method formulated in the cell-based grid, generate numerical solutions compatible with analytical ones within the range errors determinate by dispersion analysis. Our simulator will be embedded in an inversion scheme for joint seismic- electromagnetic imaging. It also offers possibilities for preconditioning the seismic wave propagation problems in the frequency domain. References. Shin, C. & Cha, Y. (2009), Waveform inversion in the Laplace-Fourier domain, Geophys. J. Int. 177(3), 1067- 1079. Shin, C. & Cha, Y. H. (2008), Waveform inversion in the Laplace domain, Geophys. J. Int. 173(3), 922-931. Commer, M. & Newman, G. (2008), New advances in three-dimensional controlled-source electromagnetic inversion, Geophys. J. Int. 172(2), 513-535. Newman, G. A., Commer, M. & Carazzone, J. J. (2010), Imaging CSEM data in the presence of electrical anisotropy, Geophysics, in press.
NASA Astrophysics Data System (ADS)
Obermann, Anne; Planès, Thomas; Hadziioannou, Céline; Campillo, Michel
2016-07-01
In the context of seismic monitoring, recent studies made successful use of seismic coda waves to locate medium changes on the horizontal plane. Locating the depth of the changes, however, remains a challenge. In this paper, we use 3-D wavefield simulations to address two problems: firstly, we evaluate the contribution of surface and body wave sensitivity to a change at depth. We introduce a thin layer with a perturbed velocity at different depths and measure the apparent relative velocity changes due to this layer at different times in the coda and for different degrees of heterogeneity of the model. We show that the depth sensitivity can be modelled as a linear combination of body- and surface-wave sensitivity. The lapse-time dependent sensitivity ratio of body waves and surface waves can be used to build 3-D sensitivity kernels for imaging purposes. Secondly, we compare the lapse-time behavior in the presence of a perturbation in horizontal and vertical slabs to address, for instance, the origin of the velocity changes detected after large earthquakes.
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
3D Ultrasonic Wave Simulations for Structural Health Monitoring
NASA Technical Reports Server (NTRS)
Campbell, Leckey Cara A/; Miler, Corey A.; Hinders, Mark K.
2011-01-01
Structural health monitoring (SHM) for the detection of damage in aerospace materials is an important area of research at NASA. Ultrasonic guided Lamb waves are a promising SHM damage detection technique since the waves can propagate long distances. For complicated flaw geometries experimental signals can be difficult to interpret. High performance computing can now handle full 3-dimensional (3D) simulations of elastic wave propagation in materials. We have developed and implemented parallel 3D elastodynamic finite integration technique (3D EFIT) code to investigate ultrasound scattering from flaws in materials. EFIT results have been compared to experimental data and the simulations provide unique insight into details of the wave behavior. This type of insight is useful for developing optimized experimental SHM techniques. 3D EFIT can also be expanded to model wave propagation and scattering in anisotropic composite materials.
Numerical simulation of 3D breaking waves
NASA Astrophysics Data System (ADS)
Fraunie, Philippe; Golay, Frederic
2015-04-01
Numerical methods dealing with two phase flows basically can be classified in two ways : the "interface tracking" methods when the two phases are resolved separately including boundary conditions fixed at the interface and the "interface capturing" methods when a single flow is considered with variable density. Physical and numerical properties of the two approaches are discussed, based on some numerical experiments performed concerning 3D breaking waves. Acknowledgements : This research was supported by the Modtercom program of Region PACA.
Elastically deformable 3D organs for haptic surgical simulation.
Webster, Roger; Haluck, Randy; Ravenscroft, Rob; Mohler, Betty; Crouthamel, Eric; Frack, Tyson; Terlecki, Steve; Sheaffer, Jeremy
2002-01-01
This paper describes a technique for incorporating real-time elastically deformable 3D organs in haptic surgical simulators. Our system is a physically based particle model utilizing a mass-springs-damper connectivity with an implicit predictor to speed up calculations during each time step. The solution involves repeated application of Newton's 2ndd Law of motion: F = ma using an implicit solver for numerically solving the differential equations. PMID:15458154
Laplace-domain waveform modeling and inversion for the 3D acoustic-elastic coupled media
NASA Astrophysics Data System (ADS)
Shin, Jungkyun; Shin, Changsoo; Calandra, Henri
2016-06-01
Laplace-domain waveform inversion reconstructs long-wavelength subsurface models by using the zero-frequency component of damped seismic signals. Despite the computational advantages of Laplace-domain waveform inversion over conventional frequency-domain waveform inversion, an acoustic assumption and an iterative matrix solver have been used to invert 3D marine datasets to mitigate the intensive computing cost. In this study, we develop a Laplace-domain waveform modeling and inversion algorithm for 3D acoustic-elastic coupled media by using a parallel sparse direct solver library (MUltifrontal Massively Parallel Solver, MUMPS). We precisely simulate a real marine environment by coupling the 3D acoustic and elastic wave equations with the proper boundary condition at the fluid-solid interface. In addition, we can extract the elastic properties of the Earth below the sea bottom from the recorded acoustic pressure datasets. As a matrix solver, the parallel sparse direct solver is used to factorize the non-symmetric impedance matrix in a distributed memory architecture and rapidly solve the wave field for a number of shots by using the lower and upper matrix factors. Using both synthetic datasets and real datasets obtained by a 3D wide azimuth survey, the long-wavelength component of the P-wave and S-wave velocity models is reconstructed and the proposed modeling and inversion algorithm are verified. A cluster of 80 CPU cores is used for this study.
3-D study of texture and elastic anisotropy on rocks from NW Italy Ivrea zone
NASA Astrophysics Data System (ADS)
Pros, Z.; Lokajicek, T.; Prikryl, R.; Klima, K.; Nikitin, A. N.; Ivankina, T. I.; Martinkova, M.
2003-04-01
The direct measurement of physical properties of lower crustal and upper mantle rocks, which can be found on the Earth's surface, could be used for the improving of our knowledge of deep rocks. These results could be used mainly for the correction of geological and geophysical models based on the indirect data. Elastic properties of rocks are one of the most important parameters studied and could be applied in many fields of Earth sciences. In this study several quite different methods were applied to determine elastic properties. P-wave ultrasonic sounding of mafic and ultrabasic rock samples in 132 independent directions at several levels of confining pressure enable to determine elastic anisotropy of P-wave velocity. The samples were collected in nearby of Balmuccia ultra basic massif (Ivrea zone, southern Alps, NW Italy). This method revealed large directional variance of maximum P-wave velocity and different symmetric (orthorhombic vs. transversal isotropic) of elastic waves 3-D distribution, that has not been found on these rocks before. Identical samples were studied by means of neutron diffraction. Neutron diffraction provide data on CPO orientation in identical spherical samples, on which was measured P-wave velocity. Laboratory 3-D measurement of P-wave velocity thus present powerful method for detection of magmatic fabric features not visible by naked eye. One dunite sample exhibits P-wave velocity approaching to that of olivine crystal 9.8 km/s due to the strong CPO of olivine in this sample. Such observation was not done before on the natural olivine-rich rocks. It follows from the comparison of measured and calculated P-wave velocities, that these values are more reliable than data obtained from measurement in few directions only. This project was supported by Grant Agency of the Czech Republic No.: 205/01/1430.
Simulation of 3D Global Wave Propagation Through Geodynamic Models
NASA Astrophysics Data System (ADS)
Schuberth, B.; Piazzoni, A.; Bunge, H.; Igel, H.; Steinle-Neumann, G.
2005-12-01
This project aims at a better understanding of the forward problem of global 3D wave propagation. We use the spectral element program "SPECFEM3D" (Komatitsch and Tromp, 2002a,b) with varying input models of seismic velocities derived from mantle convection simulations (Bunge et al., 2002). The purpose of this approach is to obtain seismic velocity models independently from seismological studies. In this way one can test the effects of varying parameters of the mantle convection models on the seismic wave field. In order to obtain the seismic velocities from the temperature field of the geodynamical simulations we follow a mineral physics approach. Assuming a certain mantle composition (e.g. pyrolite with CMASF composition) we compute the stable phases for each depth (i.e. pressure) and temperature by system Gibbs free energy minimization. Elastic moduli and density are calculated from the equations of state of the stable mineral phases. For this we use a mineral physics database derived from calorimetric experiments (enthalphy and entropy of formation, heat capacity) and EOS parameters.
Subduction zone guided waves: 3D modelling and attenuation effects
NASA Astrophysics Data System (ADS)
Garth, T.; Rietbrock, A.
2013-12-01
Waveform modelling is an important tool for understanding complex seismic structures such as subduction zone waveguides. These structures are often simplified to 2D structures for modelling purposes to reduce computational costs. In the case of subduction zone waveguide affects, 2D models have shown that dispersed arrivals are caused by a low velocity waveguide, inferred to be subducted oceanic crust and/or hydrated outer rise normal faults. However, due to the 2D modelling limitations the inferred seismic properties such as velocity contrast and waveguide thickness are still debated. Here we test these limitations with full 3D waveform modelling. For waveguide effects to be observable the waveform must be accurately modelled to relatively high frequencies (> 2 Hz). This requires a small grid spacing due to the high seismic velocities present in subduction zones. A large area must be modelled as well due to the long propagation distances (400 - 600 km) of waves interacting with subduction zone waveguides. The combination of the large model area and small grid spacing required means that these simulations require a large amount of computational resources, only available at high performance computational centres like the UK National super computer HECTOR (used in this study). To minimize the cost of modelling for such a large area, the width of the model area perpendicular to the subduction trench (the y-direction) is made as small as possible. This reduces the overall volume of the 3D model domain. Therefore the wave field is simulated in a model ';corridor' of the subduction zone velocity structure. This introduces new potential sources of error particularly from grazing wave side reflections in the y-direction. Various dampening methods are explored to reduce these grazing side reflections, including perfectly matched layers (PML) and more traditional exponential dampening layers. Defining a corridor model allows waveguide affects to be modelled up to at least 2
Towards an Anisotropic Whole Mantle 3D Elastic Velocity Model
NASA Astrophysics Data System (ADS)
Panning, M. P.; Romanowicz, B.; Gung, Y.
2001-12-01
Many studies have documented the existence of anisotropy in the earth's upper mantle, concentrated in the top 200 km. This evidence comes from the study of surface waves as well as shear wave splitting. There is also evidence for shear wave splitting in D", at least in well sampled regions. There are some hints of anisotropy at the base of the transition zone. Tomographic models of the upper mantle have been developed with simplifying assumptions about the nature of the anisotropy, in order to minimize the number of free parameters in the inversions. Some assume transverse isotropy (e.g Ekström and Dziewonski, 1997), others include additional degrees of freedom with some realistic constraints on mineralogy (e.g. Montagner and Tanimoto, 1991). Our goal is to investigate anisotropy in the whole mantle, using the framework of waveform inversion, and the nonlinear asymptotic mode coupling theory (NACT), previously developed and applied to the construction of whole-mantle SH velocity models (Li and Romanowicz, 1996; Mégnin and Romanowicz, 2000). For this we require a 3 component dataset, and we have extended our automatic transverse (T) component wavepicking procedures to the vertical (Z) and longitudinal (L) component - a non-trivial task given the large number of phases present in the coupled P-SV system. A useful initial assumption, for which the theory has been readily adapted, is that of transverse isotropy. As a first step towards this, we have been investigating inversions using T component and Z,L component data separately. In particular, this allows us to explore the sampling that can be achieved with Z,L component data alone in the deepest part of the mantle. Indeed, D" is in general much better sampled in SH than in SV, owing to the availability of SHdiff at large distances, while SVdiff decays more rapidly due to mantle-core coupling. We present the results of our resolution experiments and discuss the differences between the 3D SV model obtained in well
3-D FDTD simulation of shear waves for evaluation of complex modulus imaging.
Orescanin, Marko; Wang, Yue; Insana, Michael
2011-02-01
The Navier equation describing shear wave propagation in 3-D viscoelastic media is solved numerically with a finite differences time domain (FDTD) method. Solutions are formed in terms of transverse scatterer velocity waves and then verified via comparison to measured wave fields in heterogeneous hydrogel phantoms. The numerical algorithm is used as a tool to study the effects on complex shear modulus estimation from wave propagation in heterogeneous viscoelastic media. We used an algebraic Helmholtz inversion (AHI) technique to solve for the complex shear modulus from simulated and experimental velocity data acquired in 2-D and 3-D. Although 3-D velocity estimates are required in general, there are object geometries for which 2-D inversions provide accurate estimations of the material properties. Through simulations and experiments, we explored artifacts generated in elastic and dynamic-viscous shear modulus images related to the shear wavelength and average viscosity. PMID:21342824
3D Guided Wave Motion Analysis on Laminated Composites
NASA Technical Reports Server (NTRS)
Tian, Zhenhua; Leckey, Cara; Yu, Lingyu
2013-01-01
Ultrasonic guided waves have proved useful for structural health monitoring (SHM) and nondestructive evaluation (NDE) due to their ability to propagate long distances with less energy loss compared to bulk waves and due to their sensitivity to small defects in the structure. Analysis of actively transmitted ultrasonic signals has long been used to detect and assess damage. However, there remain many challenging tasks for guided wave based SHM due to the complexity involved with propagating guided waves, especially in the case of composite materials. The multimodal nature of the ultrasonic guided waves complicates the related damage analysis. This paper presents results from parallel 3D elastodynamic finite integration technique (EFIT) simulations used to acquire 3D wave motion in the subject laminated carbon fiber reinforced polymer composites. The acquired 3D wave motion is then analyzed by frequency-wavenumber analysis to study the wave propagation and interaction in the composite laminate. The frequency-wavenumber analysis enables the study of individual modes and visualization of mode conversion. Delamination damage has been incorporated into the EFIT model to generate "damaged" data. The potential for damage detection in laminated composites is discussed in the end.
Wave-CAIPI for Highly Accelerated 3D Imaging
Bilgic, Berkin; Gagoski, Borjan A.; Cauley, Stephen F.; Fan, Audrey P.; Polimeni, Jonathan R.; Grant, P. Ellen; Wald, Lawrence L.; Setsompop, Kawin
2014-01-01
Purpose To introduce the Wave-CAIPI (Controlled Aliasing in Parallel Imaging) acquisition and reconstruction technique for highly accelerated 3D imaging with negligible g-factor and artifact penalties. Methods The Wave-CAIPI 3D acquisition involves playing sinusoidal gy and gz gradients during the readout of each kx encoding line, while modifying the 3D phase encoding strategy to incur inter-slice shifts as in 2D-CAIPI acquisitions. The resulting acquisition spreads the aliasing evenly in all spatial directions, thereby taking full advantage of 3D coil sensitivity distribution. By expressing the voxel spreading effect as a convolution in image space, an efficient reconstruction scheme that does not require data gridding is proposed. Rapid acquisition and high quality image reconstruction with Wave-CAIPI is demonstrated for high-resolution magnitude and phase imaging and Quantitative Susceptibility Mapping (QSM). Results Wave-CAIPI enables full-brain gradient echo (GRE) acquisition at 1 mm isotropic voxel size and R=3×3 acceleration with maximum g-factors of 1.08 at 3T, and 1.05 at 7T. Relative to the other advanced Cartesian encoding strategies 2D-CAIPI and Bunched Phase Encoding, Wave-CAIPI yields up to 2-fold reduction in maximum g-factor for 9-fold acceleration at both field strengths. Conclusion Wave-CAIPI allows highly accelerated 3D acquisitions with low artifact and negligible g-factor penalties, and may facilitate clinical application of high-resolution volumetric imaging. PMID:24986223
Direct measurement of 3D elastic anisotropy on rocks from the Ivrea zone (Southern Alps, NW Italy)
NASA Astrophysics Data System (ADS)
Pros, Z.; Lokajíček, T.; Přikryl, R.; Klíma, K.
2003-07-01
Lower crustal and upper mantle rocks exposed at the earth's surface present direct possibility to measure their physical properties that must be, in other cases, interpreted using indirect methods. The results of these direct measurements can be then used for the corrections of models based on the indirect data. Elastic properties are among the most important parameters studied in geophysics and employed in many fields of earth sciences. In laboratory, dynamic elastic properties are commonly tested in three mutually perpendicular directions. The spatial distribution of P- and S-wave velocities are then computed using textural data, modal composition, density and elastic constants. During such computation, it is virtually impossible to involve all microfabric parameters like different types of microcracking, micropores, mineral alteration or quality of grain boundaries. In this study, complete 3D ultrasonic transmission of spherical samples in 132 independent directions at several levels of confining pressure up to 400 MPa has been employed for study of selected mafic and ultrabasic rocks sampled in and nearby Balmuccia ultrabasic massif (Ivrea zone, Southern Alps, NW Italy). This method revealed large directional variance of maximum P-wave velocity and different symmetries (orthorhombic vs. transversal isotropic) of elastic waves 3D distribution that has not been recorded on these rocks before. Moreover, one dunite sample exhibits P-wave velocity approaching to that of olivine single crystal being interpreted as influence of CPO.
NASA Astrophysics Data System (ADS)
Saxena, Nishank; Mavko, Gary
2016-03-01
Estimation of elastic rock moduli using 2D plane strain computations from thin sections has several numerical and analytical advantages over using 3D rock images, including faster computation, smaller memory requirements, and the availability of cheap thin sections. These advantages, however, must be weighed against the estimation accuracy of 3D rock properties from thin sections. We present a new method for predicting elastic properties of natural rocks using thin sections. Our method is based on a simple power-law transform that correlates computed 2D thin section moduli and the corresponding 3D rock moduli. The validity of this transform is established using a dataset comprised of FEM-computed elastic moduli of rock samples from various geologic formations, including Fontainebleau sandstone, Berea sandstone, Bituminous sand, and Grossmont carbonate. We note that using the power-law transform with a power-law coefficient between 0.4-0.6 contains 2D moduli to 3D moduli transformations for all rocks that are considered in this study. We also find that reliable estimates of P-wave (Vp) and S-wave velocity (Vs) trends can be obtained using 2D thin sections.
Robust elastic 2D/3D geometric graph matching
NASA Astrophysics Data System (ADS)
Serradell, Eduard; Kybic, Jan; Moreno-Noguer, Francesc; Fua, Pascal
2012-02-01
We present an algorithm for geometric matching of graphs embedded in 2D or 3D space. It is applicable for registering any graph-like structures appearing in biomedical images, such as blood vessels, pulmonary bronchi, nerve fibers, or dendritic arbors. Our approach does not rely on the similarity of local appearance features, so it is suitable for multimodal registration with a large difference in appearance. Unlike earlier methods, the algorithm uses edge shape, does not require an initial pose estimate, can handle partial matches, and can cope with nonlinear deformations and topological differences. The matching consists of two steps. First, we find an affine transform that roughly aligns the graphs by exploring the set of all consistent correspondences between the nodes. This can be done at an acceptably low computational expense by using parameter uncertainties for pruning, backtracking as needed. Parameter uncertainties are updated in a Kalman-like scheme with each match. In the second step we allow for a nonlinear part of the deformation, modeled as a Gaussian Process. Short sequences of edges are grouped into superedges, which are then matched between graphs. This allows for topological differences. A maximum consistent set of superedge matches is found using a dedicated branch-and-bound solver, which is over 100 times faster than a standard linear programming approach. Geometrical and topological consistency of candidate matches is determined in a fast hierarchical manner. We demonstrate the effectiveness of our technique at registering angiography and retinal fundus images, as well as neural image stacks.
Importance of a 3D forward modeling tool for surface wave analysis methods
NASA Astrophysics Data System (ADS)
Pageot, Damien; Le Feuvre, Mathieu; Donatienne, Leparoux; Philippe, Côte; Yann, Capdeville
2016-04-01
Since a few years, seismic surface waves analysis methods (SWM) have been widely developed and tested in the context of subsurface characterization and have demonstrated their effectiveness for sounding and monitoring purposes, e.g., high-resolution tomography of the principal geological units of California or real time monitoring of the Piton de la Fournaise volcano. Historically, these methods are mostly developed under the assumption of semi-infinite 1D layered medium without topography. The forward modeling is generally based on Thomson-Haskell matrix based modeling algorithm and the inversion is driven by Monte-Carlo sampling. Given their efficiency, SWM have been transfered to several scale of which civil engineering structures in order to, e.g., determine the so-called V s30 parameter or assess other critical constructional parameters in pavement engineering. However, at this scale, many structures may often exhibit 3D surface variations which drastically limit the efficiency of SWM application. Indeed, even in the case of an homogeneous structure, 3D geometry can bias the dispersion diagram of Rayleigh waves up to obtain discontinuous phase velocity curves which drastically impact the 1D mean velocity model obtained from dispersion inversion. Taking advantages of high-performance computing center accessibility and wave propagation modeling algorithm development, it is now possible to consider the use of a 3D elastic forward modeling algorithm instead of Thomson-Haskell method in the SWM inversion process. We use a parallelized 3D elastic modeling code based on the spectral element method which allows to obtain accurate synthetic data with very low numerical dispersion and a reasonable numerical cost. In this study, we choose dike embankments as an illustrative example. We first show that their longitudinal geometry may have a significant effect on dispersion diagrams of Rayleigh waves. Then, we demonstrate the necessity of 3D elastic modeling as a forward
3D modeling of ultrasonic wave interaction with disbonds and weak bonds
NASA Astrophysics Data System (ADS)
Leckey, C.; Hinders, M.
2012-05-01
Ultrasonic techniques, such as the use of guided waves, can be ideal for finding damage in the plate and pipe-like structures used in aerospace applications. However, the interaction of waves with real flaw types and geometries can lead to experimental signals that are difficult to interpret. 3-dimensional (3D) elastic wave simulations can be a powerful tool in understanding the complicated wave scattering involved in flaw detection and for optimizing experimental techniques. We have developed and implemented parallel 3D elastodynamic finite integration technique (3D EFIT) code to investigate Lamb wave scattering from realistic flaws. This paper discusses simulation results for an aluminum-aluminum diffusion disbond and an aluminum-epoxy disbond and compares results from the disbond case to the common artificial flaw type of a flat-bottom hole. The paper also discusses the potential for extending the 3D EFIT equations to incorporate physics-based weak bond models for simulating wave scattering from weak adhesive bonds.
3D Modeling of Ultrasonic Wave Interaction with Disbonds and Weak Bonds
NASA Technical Reports Server (NTRS)
Leckey, C.; Hinders, M.
2011-01-01
Ultrasonic techniques, such as the use of guided waves, can be ideal for finding damage in the plate and pipe-like structures used in aerospace applications. However, the interaction of waves with real flaw types and geometries can lead to experimental signals that are difficult to interpret. 3-dimensional (3D) elastic wave simulations can be a powerful tool in understanding the complicated wave scattering involved in flaw detection and for optimizing experimental techniques. We have developed and implemented parallel 3D elastodynamic finite integration technique (3D EFIT) code to investigate Lamb wave scattering from realistic flaws. This paper discusses simulation results for an aluminum-aluminum diffusion disbond and an aluminum-epoxy disbond and compares results from the disbond case to the common artificial flaw type of a flat-bottom hole. The paper also discusses the potential for extending the 3D EFIT equations to incorporate physics-based weak bond models for simulating wave scattering from weak adhesive bonds.
NASA Astrophysics Data System (ADS)
Duru, Kenneth; Dunham, Eric M.
2016-01-01
Dynamic propagation of shear ruptures on a frictional interface in an elastic solid is a useful idealization of natural earthquakes. The conditions relating discontinuities in particle velocities across fault zones and tractions acting on the fault are often expressed as nonlinear friction laws. The corresponding initial boundary value problems are both numerically and computationally challenging. In addition, seismic waves generated by earthquake ruptures must be propagated for many wavelengths away from the fault. Therefore, reliable and efficient numerical simulations require both provably stable and high order accurate numerical methods. We present a high order accurate finite difference method for: a) enforcing nonlinear friction laws, in a consistent and provably stable manner, suitable for efficient explicit time integration; b) dynamic propagation of earthquake ruptures along nonplanar faults; and c) accurate propagation of seismic waves in heterogeneous media with free surface topography. We solve the first order form of the 3D elastic wave equation on a boundary-conforming curvilinear mesh, in terms of particle velocities and stresses that are collocated in space and time, using summation-by-parts (SBP) finite difference operators in space. Boundary and interface conditions are imposed weakly using penalties. By deriving semi-discrete energy estimates analogous to the continuous energy estimates we prove numerical stability. The finite difference stencils used in this paper are sixth order accurate in the interior and third order accurate close to the boundaries. However, the method is applicable to any spatial operator with a diagonal norm satisfying the SBP property. Time stepping is performed with a 4th order accurate explicit low storage Runge-Kutta scheme, thus yielding a globally fourth order accurate method in both space and time. We show numerical simulations on band limited self-similar fractal faults revealing the complexity of rupture dynamics
Instability and Wave Propagation in Structured 3D Composites
NASA Astrophysics Data System (ADS)
Kaynia, Narges; Fang, Nicholas X.; Boyce, Mary C.
2014-03-01
Many structured composites found in nature possess undulating and wrinkled interfacial layers that regulate mechanical, chemical, acoustic, adhesive, thermal, electrical and optical functions of the material. This research focused on the complex instability and wrinkling pattern arising in 3D structured composites and the effect of the buckling pattern on the overall structural response. The 3D structured composites consisted of stiffer plates supported by soft matrix on both sides. Compression beyond the critical strain led to complex buckling patterns in the initially straight plates. The motivation of our work is to elaborate the formation of a system of prescribed periodic scatterers (metamaterials) due to buckling, and their effect to interfere wave propagation through the metamaterial structures. Such metamaterials made from elastomers enable large reversible deformation and, as a result, significant changes of the wave propagation properties. We developed analytical and finite element models to capture various aspects of the instability mechanism. Mechanical experiments were designed to further explore the modeling results. The ability to actively alter the 3D composite structure can enable on-demand tunability of many different functions, such as active control of wave propagation to create band-gaps and waveguides.
Upscaling small heterogeneities for seismic wave propagation in 3D complex media
NASA Astrophysics Data System (ADS)
Cupillard, P.; Capdeville, Y.
2012-04-01
Seismic waves propagating in the Earth are affected by different sizes of heterogeneities. When modelling these waves (using numerical methods such as the SEM), taking into account heterogeneities that are much smaller than the minimum wavelength is a challenge because meshing small heterogeneities often requires important efforts and leads to high numerical costs. In this work, we present a technique which allows to upscale the small heterogeneities that can lie in an elastic medium. This technique yields a smooth effective medium and effective equations. We describe its implementation in the 3D case and we show relevant examples.
Protrusive waves guide 3D cell migration along nanofibers
Guetta-Terrier, Charlotte; Monzo, Pascale; Zhu, Jie; Long, Hongyan; Venkatraman, Lakshmi; Zhou, Yue; Wang, PeiPei; Chew, Sing Yian; Mogilner, Alexander
2015-01-01
In vivo, cells migrate on complex three-dimensional (3D) fibrous matrices, which has made investigation of the key molecular and physical mechanisms that drive cell migration difficult. Using reductionist approaches based on 3D electrospun fibers, we report for various cell types that single-cell migration along fibronectin-coated nanofibers is associated with lateral actin-based waves. These cyclical waves have a fin-like shape and propagate up to several hundred micrometers from the cell body, extending the leading edge and promoting highly persistent directional movement. Cells generate these waves through balanced activation of the Rac1/N-WASP/Arp2/3 and Rho/formins pathways. The waves originate from one major adhesion site at leading end of the cell body, which is linked through actomyosin contractility to another site at the back of the cell, allowing force generation, matrix deformation and cell translocation. By combining experimental and modeling data, we demonstrate that cell migration in a fibrous environment requires the formation and propagation of dynamic, actin based fin-like protrusions. PMID:26553933
3D simulation of seismic wave propagation around a tunnel using the spectral element method
NASA Astrophysics Data System (ADS)
Lambrecht, L.; Friederich, W.
2010-05-01
We model seismic wave propagation in the environment of a tunnel for later application to reconnaissance. Elastic wave propagation can be simulated by different numerical techniques such as finite differences and pseudospectral methods. Their disadvantage is the lack of accuracy on free surfaces, numerical dispersion and inflexibility of the mesh. Here we use the software package SPECFEM3D_SESAME in an svn development version, which is based on the spectral element method (SEM) and can handle complex mesh geometries. A weak form of the elastic wave equation leads to a linear system of equations with a diagonal mass matrix, where the free surface boundary of the tunnel can be treated under realistic conditions and can be effectively implemented in parallel. We have designed a 3D external mesh including a tunnel and realistic features such as layers and holes to simulate elastic wave propagation in the zone around the tunnel. The source is acting at the tunnel surface so that we excite Rayleigh waves which propagate to the front face of the tunnel. A conversion takes place and a high amplitude S-wave is radiated in the direction of the tunnel axis. Reflections from perturbations in front of the tunnel can be measured by receivers implemented on the tunnel face. For a shallow tunnel the land surface has high influence on the wave propagation. By implementing additional receivers at this surface we intent to improve the prediction. It shows that the SEM is very capable to handle the complex geometry of the model and especially incorporates the free surfaces of the model.
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.
NASA Astrophysics Data System (ADS)
Bognet, B.; Leygue, A.; Chinesta, F.; Poitou, A.
2011-01-01
In this paper, we focus on the simulation of linear elastic behaviour of plates using a 3D approach which numerical cost only scales like a 2D one. In the case of plates, the kinematic hypothesis introduced in plate theories to go from 3D to 2D is usually unsatisfactory where one cannot rely on St Venant's principle (usually close to the plate edges). We propose to apply the PGD (Proper Generalized Decomposition) method [1] to the simulation of the linear elastic behavior of plates. This method allows us to separately search for the in-plane and the out-of plane contributions to the 3D solution, yielding significant savings in computational cost. The method is validated on a simple case and its full potential is then presented for the simulation of the behavior of laminated composite plates.
An Application of the Method of Arbitrary Lines to 3D Elastic Stress Analysis
NASA Astrophysics Data System (ADS)
Kaminishi, Ken; Ando, Ryuma
The MAL (Method of Arbitrary Lines) is a technique of reducing a partial differential equation to a system of ordinary differential equations. It is known that relevant use of this procedure yields high accuracy in some problems of two-dimensional elasticity and elastoplasticity. Since the basic concept of MAL is simple and based on generality, it is expected that many problems in other fields will be effectively solvable by this method. In this study, we consider the application of MAL to 3D (three-dimensional) elasticity analysis. We first give a MAL formulation of 3D elasticity problems, and demonstrate its effectiveness and accuracy for a typical problem. The reported numerical results are compared with the exact solution or that of the finite element method (FEM).
ZIP3D: An elastic and elastic-plastic finite-element analysis program for cracked bodies
NASA Technical Reports Server (NTRS)
Shivakumar, K. N.; Newman, J. C., Jr.
1990-01-01
ZIP3D is an elastic and an elastic-plastic finite element program to analyze cracks in three dimensional solids. The program may also be used to analyze uncracked bodies or multi-body problems involving contacting surfaces. For crack problems, the program has several unique features including the calculation of mixed-mode strain energy release rates using the three dimensional virtual crack closure technique, the calculation of the J integral using the equivalent domain integral method, the capability to extend the crack front under monotonic or cyclic loading, and the capability to close or open the crack surfaces during cyclic loading. The theories behind the various aspects of the program are explained briefly. Line-by-line data preparation is presented. Input data and results for an elastic analysis of a surface crack in a plate and for an elastic-plastic analysis of a single-edge-crack-tension specimen are also presented.
Evaluation of Fracture Azimuth by EM Wave and Elastic Wave
NASA Astrophysics Data System (ADS)
Feng, X.; Wang, Q.; Liu, C.; Lu, Q.; Zeng, Z.; Liang, W.; Yu, Y.; Ren, Q.
2013-12-01
Fracture system plays an important role in the development of underground energy, for example enhanced geothermal system (EGS), oil shale and shale gas, etc. Therefore, it becomes more and more important to detect and evaluate the fracture system. Geophysical prospecting is an useful method to evaluate the characteristics of the subsurface fractures. Currently, micro-seismology, multi-wave seismic exploration, and electromagnetic (EM) survey are reported to be used for the purpose. We are studying a method using both elastic wave and EM wave to detect and evaluate the fracture azimuth in laboratory. First, we build a 3D horizontal transverse isotropy (HTI) model, shown in the figure 1, by dry parallel fractures system, which was constructed by plexiglass plates and papers. Then, we used the ultrasonic system to obtain reflected S-wave data. Depending on the shear wave splitting, we evaluated the fracture azimuth by the algorithm of Pearson correlation coefficient. In addition, we used the full Polarimetric ultra wide band electromagnetic (FP-UWB-EM) wave System, shown in the figure 2, to obtain full polarimetric reflected EM-wave data. Depending on the rotation of the EM wave polarimetry, we evaluated the fracture azimuth by the the ration between maximum amplitude of co-polarimetric EM wave and maximum amplitude of cross-polarimetric EM wave. Finally, we used both EM-wave data and S-wave data to evaluate the fracture azimuth by the method of cross plot and statistical mathematics. To sum up, we found that FP-UWB-EM wave can be used to evaluated the fracture azimuth and is more accurate than ultrasound wave. Also joint evaluation using both data could improve the precision.
Simulation of 3D Seismic Wave Propagation with Volcano Topography
NASA Astrophysics Data System (ADS)
Ripperger, J.; Igel, H.; Wassermann, J.
2001-12-01
We investigate the possibilities of using three-dimensional finite difference (FD) methods for numerical simulation of the seismic wave field at active volcanoes. We put special emphasis on the implementation of the boundary conditions for free surface topography. We compare two different approaches to solve the free surface boundary conditions. The algorithms are implemented on parallel hardware and have been tested for correctness and stability. We apply them to smooth artificial topographies and to the real topography of Mount Merapi, Indonesia. We conclude, that grid stretching type methods (e.g. Hestholm & Ruud, 1994) are not well suited for realistic volcano topography as they tend to become unstable for large topographic gradients. The representation of topography through staircase shaped grids (Ohminato & Chouet, 1997) results in stable calculations, while demanding very fine gridding. The simulations show the effects of a three-dimensional surface topography on elastic wave propagation. Ground motion at the surface is severely affected by topography. If neglected, this may jeopardize attempts to determine source location by analyzing particle motion. Numerical studies like this can help to understand wave propagation phenomena observed on field recordings in volcano seismology. Future studies will aim at separating the wave effects of internal scattering, topography and sources (tremors, tectonic events, pyroclastic flows).
ATHENA 3D: A finite element code for ultrasonic wave propagation
NASA Astrophysics Data System (ADS)
Rose, C.; Rupin, F.; Fouquet, T.; Chassignole, B.
2014-04-01
The understanding of wave propagation phenomena requires use of robust numerical models. 3D finite element (FE) models are generally prohibitively time consuming. However, advances in computing processor speed and memory allow them to be more and more competitive. In this context, EDF R&D developed the 3D version of the well-validated FE code ATHENA2D. The code is dedicated to the simulation of wave propagation in all kinds of elastic media and in particular, heterogeneous and anisotropic materials like welds. It is based on solving elastodynamic equations in the calculation zone expressed in terms of stress and particle velocities. The particularity of the code relies on the fact that the discretization of the calculation domain uses a Cartesian regular 3D mesh while the defect of complex geometry can be described using a separate (2D) mesh using the fictitious domains method. This allows combining the rapidity of regular meshes computation with the capability of modelling arbitrary shaped defects. Furthermore, the calculation domain is discretized with a quasi-explicit time evolution scheme. Thereby only local linear systems of small size have to be solved. The final step to reduce the computation time relies on the fact that ATHENA3D has been parallelized and adapted to the use of HPC resources. In this paper, the validation of the 3D FE model is discussed. A cross-validation of ATHENA 3D and CIVA is proposed for several inspection configurations. The performances in terms of calculation time are also presented in the cases of both local computer and computation cluster use.
Nonlinear dynamics of Airy-vortex 3D wave packets: emission of vortex light waves.
Driben, Rodislav; Meier, Torsten
2014-10-01
The dynamics of 3D Airy-vortex wave packets is studied under the action of strong self-focusing Kerr nonlinearity. Emissions of nonlinear 3D waves out of the main wave packets with the topological charges were demonstrated. Because of the conservation of the total angular momentum, charges of the emitted waves are equal to those carried by the parental light structure. The rapid collapse imposes a severe limitation on the propagation of multidimensional waves in Kerr media. However, the structure of the Airy beam carrier allows the coupling of light from the leading, most intense peak into neighboring peaks and consequently strongly postpones the collapse. The dependence of the critical input amplitude for the appearance of a fast collapse on the beam width is studied for wave packets with zero and nonzero topological charges. Wave packets carrying angular momentum are found to be much more resistant to the rapid collapse. PMID:25360922
Numerical study of elastic turbulence in a 3D curvilinear micro-channel
NASA Astrophysics Data System (ADS)
Zhang, Hongna; Kunugi, Tomoaki; Li, Fengchen
2012-11-01
Elastic turbulence is an intriguing phenomenon of viscoelastic fluid flow, and dominated by the strong nonlinear elasticity due to the existence of flexible microstructures. It implies the possibility to generate a turbulent state (so-called an elastic turbulence) in the micro-scale devices by introducing the viscoelastic fluids, which could significantly enhance the mixing efficiency therein. Several experiments have been carried out to study its characteristics and underlying physics. However, the difficulty in measuring the flow information and behaviors of the microstructures, especially in the cross section normal to the mean flow direction, limits our current understanding and controlling. In the present study, the nondimensionalization method in which the characteristic velocity is defined as the ratio of the solution viscosity to the width of the channel was adopted to simulate the elastic turbulence in the micro-scale devices. And the elastic turbulent flow was obtained numerically in the 3D curvilinear micro-channel. Therein, the characteristics of the velocity field and polymer's behavior are discussed. Moreover, the energy transfer between the kinetic energy and the polymer's elastic energy is also investigated to understand its physical mechanism. Supported by the Japan Society for the Promotion of Science research fellowship and the Ministry of Education, Culture, Sports, Science and Technology via `Energy Science in the Age of Global Warming' of Global Center of Excellence (G-COE) program (J-051).
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
Faraday wave lattice as an elastic metamaterial
NASA Astrophysics Data System (ADS)
Domino, L.; Tarpin, M.; Patinet, S.; Eddi, A.
2016-05-01
Metamaterials enable the emergence of novel physical properties due to the existence of an underlying subwavelength structure. Here, we use the Faraday instability to shape the fluid-air interface with a regular pattern. This pattern undergoes an oscillating secondary instability and exhibits spontaneous vibrations that are analogous to transverse elastic waves. By locally forcing these waves, we fully characterize their dispersion relation and show that a Faraday pattern presents an effective shear elasticity. We propose a physical mechanism combining surface tension with the Faraday structured interface that quantitatively predicts the elastic wave phase speed, revealing that the liquid interface behaves as an elastic metamaterial.
Faraday wave lattice as an elastic metamaterial.
Domino, L; Tarpin, M; Patinet, S; Eddi, A
2016-05-01
Metamaterials enable the emergence of novel physical properties due to the existence of an underlying subwavelength structure. Here, we use the Faraday instability to shape the fluid-air interface with a regular pattern. This pattern undergoes an oscillating secondary instability and exhibits spontaneous vibrations that are analogous to transverse elastic waves. By locally forcing these waves, we fully characterize their dispersion relation and show that a Faraday pattern presents an effective shear elasticity. We propose a physical mechanism combining surface tension with the Faraday structured interface that quantitatively predicts the elastic wave phase speed, revealing that the liquid interface behaves as an elastic metamaterial. PMID:27300815
Elastic registration using 3D ChainMail: application to virtual colonoscopy
NASA Astrophysics Data System (ADS)
Castro-Pareja, Carlos R.; Daly, Barry; Shekhar, Raj
2006-03-01
We present an elastic registration algorithm based on local deformations modeled using cubic B-splines and controlled using 3D ChainMail. Our algorithm eliminates the appearance of folding artifacts and allows local rigidity and compressibility control independent of the image similarity metric being used. 3D ChainMail propagates large internal deformations between neighboring B-Spline control points, thereby preserving the topology of the transformed image without requiring the addition of penalty terms based on rigidity of the transformation field to the equation used to maximize image similarity. A novel application to virtual colonoscopy is presented where the algorithm is used to significantly improve cross-localization between colon locations in prone and supine CT images.
Noninvasive 3D elasticity mapping using phase-stabilized optical coherence elastography
NASA Astrophysics Data System (ADS)
Singh, Manmohan; Li, Jiasong; Wang, Shang; Twa, Michael; Larin, Kirill V.
2015-03-01
We demonstrate a novel method for noninvasive elasticity mapping in three dimensions using phase stabilized swept source optical coherence elastography (PhS-SSOCE). By calculating the velocity in all radial directions from the origin of the induced shear wave, a volumetric elasticity map of the sample was generated. Due to the submicrometer spatial sensitivity of PhS-SSOCE, the loading force and the induced deformation amplitude can be minimal, thus preserving the structure and function of delicate tissues such as the cornea and sclera of the eye. Tissue mimicking agar phantoms were utilized for proof of concept testing and the results show that this method can noninvasively provide a three dimensional estimation of sample elasticity.
NASA Astrophysics Data System (ADS)
Okaya, D. A.; Van Avendonk, H. J.
2013-12-01
Recent anisotropy studies at scales ranging from crust to full mantle have recognized the importance of 3D anisotropy geometry and heterogeneity as well as variability in anisotropic symmetry and orientation (tilt) of the Earth. The strong relationship between seismic anisotropy and geodynamic processes highlights the need to construct realistic Earth models that can explain observations of anisotropy in modern seismic data sets. For example, ray paths through a mantle slab window or a mountain belt may show that the crust or mantle exhibits low-order anisotropy due to a history of deformation and the development of tectonic fabrics. Observed traveltimes might not be fit with simple Transverse Isotropy (TI), so realistic calculations require an Earth model that accurately describes the wave speeds of compressional and shear waves. We have developed an anisotropic traveltime solver that allows for full 3D heterogeneity of anisotropy tensors, degrees of symmetry, and arbitrary orientation. This traveltime solver is based on the robust shortest path method (SPM) and a ray-bending algorithm that were previously applied to isotropic media (e.g., Van Avendonk et al., 2001). Instead of using an isotropic description of the seismic wave velocity, we define the full elastic tensor at each location in the model. The directional seismic velocity can subsequently be extracted using solutions of the Christoffel equations. For computational efficiency, we calculate all directional seismic velocities at each model node before the start of ray tracing. As we calculate a new ray segment, this information is quickly retrieved. We use these directional velocity maps to separately describe the propagation of compressional (P) and shear (S) body waves in anisotropic media and to subsequently calculate their traveltimes. Patterns within the velocity maps represent tensor symmetries and tilts, allowing for the construction of discretized large-scale 3D LPO flow fields or fabric
Numerical homogenization for seismic wave propagation in 3D geological media
NASA Astrophysics Data System (ADS)
Cupillard, P.; Capdeville, Y.; Botella, A.
2014-12-01
Despite the important increase of the computational power in the last decades, simulating the seismic wave propagation through realistic geological models is still a challenge. By realistic models we here mean 3D media in which a broad variety (in terms of amplitude and extent) of heterogeneities lies, including discontinuities with complex geometry such as faulted and folded horizons, intrusive geological contacts and fault systems. To perform accurate numerical simulations, these discontinuities require complicated meshes which usually contain extremely small elements, yielding large, sometimes prohibitive, computation costs. Fortunately, the recent development of the non-periodic homogenization technique now enables to overcome this problem by computing smooth equivalent models for which a coarse mesh is sufficient to get an accurate wavefield. In this work, we present an efficient implementation of the technique which now allows for the homogenization of large 3D geological models. This implementation relies on a tetrahedral finite-element solution of the elasto-static equation behind the homogenization problem. Because this equation is time-independent, solving it is numerically cheaper than solving the wave equation, but it nevertheless requires some care because of the large size of the stiffness matrix arising from the fine mesh of realistic geological structures. A domain decomposition is therefore adopted. In our strategy, the obtained sub-domains overlap but they are independent so the solution within each of them can be computed either in series or in parallel. In addition, well-balanced loads, efficient search algorithms and multithreading are implemented to speed up the computation. The resulting code enables the homogenization of 3D elastic media in a time that is neglectable with respect to the simulation time of the wave propagation within. This is illustrated through a sub-surface model of the Furfooz karstic region, Belgium.
Jammed elastic shells - a 3D experimental soft frictionless granular system
NASA Astrophysics Data System (ADS)
Jose, Jissy; Blab, Gerhard A.; van Blaaderen, Alfons; Imhof, Arnout
2015-03-01
We present a new experimental system of monodisperse, soft, frictionless, fluorescent labelled elastic shells for the characterization of structure, universal scaling laws and force networks in 3D jammed matter. The interesting fact about these elastic shells is that they can reversibly deform and therefore serve as sensors of local stress in jammed matter. Similar to other soft particles, like emulsion droplets and bubbles in foam, the shells can be packed to volume fractions close to unity, which allows us to characterize the contact force distribution and universal scaling laws as a function of volume fraction, and to compare them with theoretical predictions and numerical simulations. However, our shells, unlike other soft particles, deform rather differently at large stresses. They deform without conserving their inner volume, by forming dimples at contact regions. At each contact one of the shells buckled with a dimple and the other remained spherical, closely resembling overlapping spheres. We conducted 3D quantitative analysis using confocal microscopy and image analysis routines specially developed for these particles. In addition, we analysed the randomness of the process of dimpling, which was found to be volume fraction dependent.
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.
Simulation domain size requirements for elastic response of 3D polycrystalline materials
NASA Astrophysics Data System (ADS)
Ozturk, Tugce; Stein, Clayton; Pokharel, Reeju; Hefferan, Christopher; Tucker, Harris; Jha, Sushant; John, Reji; Lebensohn, Ricardo A.; Kenesei, Peter; Suter, Robert M.; Rollett, Anthony D.
2016-01-01
A fast Fourier transform (FFT) based spectral algorithm is used to compute the full field mechanical response of polycrystalline microstructures. The field distributions in a specific region are used to determine the sensitivity of the method to the number of surrounding grains through quantification of the divergence of the field values from the largest simulation domain, as successively smaller surrounding volumes are included in the simulation. The analysis considers a mapped 3D structure where the location of interest is taken to be a particular pair of surface grains that enclose a small fatigue crack, and synthetically created statistically representative microstructures to further investigate the effect of anisotropy, loading condition, loading direction, and texture. The synthetic structures are generated via DREAM3D and the measured material is a cyclically loaded, Ni-based, low solvus high refractory (LSHR) superalloy that was characterized via 3D high energy x-ray diffraction microscopy (HEDM). Point-wise comparison of distributions in the grain pairs shows that, in order to obtain a Pearson correlation coefficient larger than 99%, the domain must extend to at least the third nearest neighbor. For an elastic FFT calculation, the stress-strain distributions are not sensitive to the shape of the domain. The main result is that convergence can be specified in terms of the number of grains surrounding a region of interest.
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.
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
Fracture imaging with converted elastic waves
Nihei, K.T.; Nakagawa, S.; Myer, L.R.
2001-05-29
This paper examines the seismic signatures of discrete, finite-length fractures, and outlines an approach for elastic, prestack reverse-time imaging of discrete fractures. The results of this study highlight the importance of incorporating fracture-generated P-S converted waves into the imaging method, and presents an alternate imaging condition that can be used in elastic reverse-time imaging when a direct wave is recorded (e.g., for crosswell and VSP acquisition geometries).
A Full-Wave Approach to Elastic and Q Tomography
NASA Astrophysics Data System (ADS)
Zhao, L.; Chen, P.
2006-12-01
Phase delays and traveltimes of seismic waves have been used extensively in seismic tomography to image the laterally heterogeneous elastic structures of the Earth. However, the amplitudes of seismic waves have not been as fully exploited. The difficulties in utilizing amplitudes in structural studies are two folds. The amplitudes of seismic waves are often affected by structural variations in a very nonlinear fashion and as a result the amplitudes are not robust data for tomography inversions. Moreover, the amplitudes of seismic waves are affected by not only the elastic structures through focusing/defocusing and scattering, but also the anelastic structures through attenuation. We propose a consistent and comprehensive approach to phase- delay and amplitude tomography inversion for the Earth's elastic and anelastic structures. We adopt a consistent definition for the phase-delay and amplitude anomalies and measure both from the same cross- correlation between synthetic and recorded seismograms. Frequency-dependent anomalies can be obtained from narrow-band filtered cross-correlagrams. We also assure consistency in interpreting the measurements in terms of structural variations by linearly relating the frequency-dependent phase-delay anomalies to both the elastic parameters to account for scattering and the Q values to account for physical dispersion; and at the same time linearly relating the frequency-dependent amplitude anomalies to the same elastic parameters and Q values to account for scattering and attenuation. We present examples of full-wave 3D sensitivity kernels for these linear relationships computed by coupled normal-mode summations, as well as results of an experimental Q tomography using regional Rayleigh waves in East Asia.
Embedding SAS approach into conjugate gradient algorithms for asymmetric 3D elasticity problems
Chen, Hsin-Chu; Warsi, N.A.; Sameh, A.
1996-12-31
In this paper, we present two strategies to embed the SAS (symmetric-and-antisymmetric) scheme into conjugate gradient (CG) algorithms to make solving 3D elasticity problems, with or without global reflexive symmetry, more efficient. The SAS approach is physically a domain decomposition scheme that takes advantage of reflexive symmetry of discretized physical problems, and algebraically a matrix transformation method that exploits special reflexivity properties of the matrix resulting from discretization. In addition to offering large-grain parallelism, which is valuable in a multiprocessing environment, the SAS scheme also has the potential for reducing arithmetic operations in the numerical solution of a reasonably wide class of scientific and engineering problems. This approach can be applied directly to problems that have global reflexive symmetry, yielding smaller and independent subproblems to solve, or indirectly to problems with partial symmetry, resulting in loosely coupled subproblems. The decomposition is achieved by separating the reflexive subspace from the antireflexive one, possessed by a special class of matrices A, A {element_of} C{sup n x n} that satisfy the relation A = PAP where P is a reflection matrix (symmetric signed permutation matrix).
3D dynamic simulation of crack propagation in extracorporeal shock wave lithotripsy
NASA Astrophysics Data System (ADS)
Wijerathne, M. L. L.; Hori, Muneo; Sakaguchi, Hide; Oguni, Kenji
2010-06-01
Some experimental observations of Shock Wave Lithotripsy(SWL), which include 3D dynamic crack propagation, are simulated with the aim of reproducing fragmentation of kidney stones with SWL. Extracorporeal shock wave lithotripsy (ESWL) is the fragmentation of kidney stones by focusing an ultrasonic pressure pulse onto the stones. 3D models with fine discretization are used to accurately capture the high amplitude shear shock waves. For solving the resulting large scale dynamic crack propagation problem, PDS-FEM is used; it provides numerically efficient failure treatments. With a distributed memory parallel code of PDS-FEM, experimentally observed 3D photoelastic images of transient stress waves and crack patterns in cylindrical samples are successfully reproduced. The numerical crack patterns are in good agreement with the experimental ones, quantitatively. The results shows that the high amplitude shear waves induced in solid, by the lithotriptor generated shock wave, play a dominant role in stone fragmentation.
Simulation of the elastic wave propagation in anisotropic microstructures
NASA Astrophysics Data System (ADS)
Bryner, Juerg; Vollmann, Jacqueline; Profunser, Dieter M.; Dual, Jurg
2007-06-01
For the interpretation of optical Pump-Probe Measurements on microstructures the wave propagation in anisotropic 3-D structures with arbitrary geometries is numerically calculated. The laser acoustic Pump-Probe technique generates bulk waves in structures in a thermo-elastic way. This method is well established for non-destructive measurements of thin films with an indepth resolution in the order of 10 nm. The Pump-Probe technique can also be used for measurements, e.g. for quality inspection of three-dimensional structures with arbitrary geometries, like MEMS components. For the interpretation of the measurements it is necessary that the wave propagation in the specimen to be inspected can be calculated. Here, the wave propagation for various geometries and materials is investigated. In the first part, the wave propagation in isotropic axisymmetric structures is simulated with a 2-D finite difference formulation. The numerical results are verified with measurements of macroscopic specimens. In a second step, the simulations are extended to 3-D structures with orthotopic material properties. The implemented code allows the calculation of the wave propagation for different orientations of the material axes (orientation of the orthotropic axes relative to the geometry of the structure). Limits of the presented algorithm are discussed and future directions of the on-going research project are presented.
A dispersive wave equation using nonlocal elasticity
NASA Astrophysics Data System (ADS)
Challamel, Noël; Rakotomanana, Lalaonirina; Le Marrec, Loïc
2009-08-01
Nonlocal continuum mechanics allows one to account for the small length scale effect that becomes significant when dealing with micro- or nano-structures. This Note investigates a model of wave propagation in a nonlocal elastic material. We show that a dispersive wave equation is obtained from a nonlocal elastic constitutive law, based on a mixture of a local and a nonlocal strain. This model comprises both the classical gradient model and the Eringen's integral model. The dynamic properties of the model are discussed, and corroborate well some recent theoretical studies published to unify both static and dynamics gradient elasticity theories. Moreover, an excellent matching of the dispersive curve of the Born-Kármán model of lattice dynamics is obtained with such nonlocal model. To cite this article: N. Challamel et al., C. R. Mecanique 337 (2009).
Wave propagation in polar elastic superlattices
NASA Astrophysics Data System (ADS)
Green, W. A.; Green, E. Rhian
1994-08-01
This paper examines the passband and stop band regions for time-periodic waves travelling normal to the layering through an infinite medium composed of alternating layers of two different elastic materials. The materials are such that the elastic energy density is a function of the strains and the strain gradients and, in consequence, a deformation gives rise to both the usual Cauchy stress and to a hyperstress or couple-stress. Such materials can exhibit a non-uniform wrinkling deformation at a free surface and similar non-uniform deformations can arise at interfaces between two different media. The presence of the strain derivatives in the elastic energy function introduces a natural length scale l into the material and the depth of the non-uniform deformation is of the order of this length scale. This model can give rise to enhanced elastic response when the layer depths are comparable with l and it is of interest as a possible mathematical model of nanolayered structures. The model also includes a non-standard set of continuity conditions at material interfaces. These arise from the elastic interaction energy of the two materials at the boundary and their effect is localized in a boundary layer whose depth is of order l. The periodic layering gives rise to displacements which are periodic with a frequency-dependent wave number, the Floquet wave number. Dispersion curves, relating circular frequency to the Floquet wave number, are obtained for different ratios of the layer depth to the natural length l and for different values of the elastic interface coupling parameters.
Simultaneous elastic parameter inversion in 2-D/3-D TTI medium combined later arrival times
NASA Astrophysics Data System (ADS)
Bai, Chao-ying; Wang, Tao; Yang, Shang-bei; Li, Xing-wang; Huang, Guo-jiao
2016-04-01
Traditional traveltime inversion for anisotropic medium is, in general, based on a "weak" assumption in the anisotropic property, which simplifies both the forward part (ray tracing is performed once only) and the inversion part (a linear inversion solver is possible). But for some real applications, a general (both "weak" and "strong") anisotropic medium should be considered. In such cases, one has to develop a ray tracing algorithm to handle with the general (including "strong") anisotropic medium and also to design a non-linear inversion solver for later tomography. Meanwhile, it is constructive to investigate how much the tomographic resolution can be improved by introducing the later arrivals. For this motivation, we incorporated our newly developed ray tracing algorithm (multistage irregular shortest-path method) for general anisotropic media with a non-linear inversion solver (a damped minimum norm, constrained least squares problem with a conjugate gradient approach) to formulate a non-linear inversion solver for anisotropic medium. This anisotropic traveltime inversion procedure is able to combine the later (reflected) arrival times. Both 2-D/3-D synthetic inversion experiments and comparison tests show that (1) the proposed anisotropic traveltime inversion scheme is able to recover the high contrast anomalies and (2) it is possible to improve the tomographic resolution by introducing the later (reflected) arrivals, but not as expected in the isotropic medium, because the different velocity (qP, qSV and qSH) sensitivities (or derivatives) respective to the different elastic parameters are not the same but are also dependent on the inclination angle.
Geodynamic background of the 2008 Wenchuan earthquake based on 3D visco-elastic numerical modelling
NASA Astrophysics Data System (ADS)
Liu, Chang; Zhu, Bojing; Yang, Xiaolin; Shi, Yaolin
2016-03-01
The 2008 Wenchuan earthquake (Mw7.9) occurred in the Longmen Shan fault zone. The stress change and crustal deformation during the accumulation period is computed using 3D finite element modelling assuming visco-elastic rheology. Our results support that the eastward movement of the Tibetan Plateau resulting from the India-Eurasia collision is obstructed at the Longmen Shan fault zone by the strong Yangtze craton. In response, the Tibetan ductile crust thickens and accumulates at the contact between the Tibetan Plateau and the Sichuan Basin. This process implies a strong uplift with the rate of about 1.8 mm/a of the upper crust and induces a stress concentration nearly at the bottom of the Longmen Shan fault zone. We believe that the stress concentration in the Longmen Shan fault zone provides a very important geodynamic background of the 2008 Wenchuan earthquake. Using numerical experiments we find that the key factor controlling this stress concentration process is the large viscosity contrast in the middle and lower crusts between the Tibetan Plateau and the Sichuan Basin. The results show that large viscosity contrast in the middle and lower crusts accelerates the stress concentration in the Longmen Shan fault zone. Fast moving lower crustal flow accelerates this stress accumulation process. During the inter-seismic period, spatially the maximum stress accumulation rate of the eastern margin of the Tibetan Plateau is located nearly at the bottom of the brittle upper crust of the Longmen Shan fault zone. The spatial distribution of the stress accumulation along the strike of the Longmen Shan fault zone is as follows: the normal stress decreases while the shear stress increases from southwest to northeast along the Longmen Shan fault zone. This stress distribution explains the thrust motion in the SW and strike-slip motion in the NE during the 2008 Wenchuan earthquake.
Characterizing the propagation of gravity waves in 3D nonlinear simulations of solar-like stars
NASA Astrophysics Data System (ADS)
Alvan, L.; Strugarek, A.; Brun, A. S.; Mathis, S.; Garcia, R. A.
2015-09-01
Context. The revolution of helio- and asteroseismology provides access to the detailed properties of stellar interiors by studying the star's oscillation modes. Among them, gravity (g) modes are formed by constructive interferences between progressive internal gravity waves (IGWs), propagating in stellar radiative zones. Our new 3D nonlinear simulations of the interior of a solar-like star allows us to study the excitation, propagation, and dissipation of these waves. Aims: The aim of this article is to clarify our understanding of the behavior of IGWs in a 3D radiative zone and to provide a clear overview of their properties. Methods: We use a method of frequency filtering that reveals the path of individual gravity waves of different frequencies in the radiative zone. Results: We are able to identify the region of propagation of different waves in 2D and 3D, to compare them to the linear raytracing theory and to distinguish between propagative and standing waves (g-modes). We also show that the energy carried by waves is distributed in different planes in the sphere, depending on their azimuthal wave number. Conclusions: We are able to isolate individual IGWs from a complex spectrum and to study their propagation in space and time. In particular, we highlight in this paper the necessity of studying the propagation of waves in 3D spherical geometry, since the distribution of their energy is not equipartitioned in the sphere.
Wave anisotropy of shear viscosity and elasticity
NASA Astrophysics Data System (ADS)
Rudenko, O. V.; Sarvazyan, A. P.
2014-11-01
The paper presents the theory of shear wave propagation in a "soft solid" material possessing anisotropy of elastic and dissipative properties. The theory is developed mainly for understanding the nature of the low-frequency acoustic characteristics of skeletal muscles, which carry important diagnostic information on the functional state of muscles and their pathologies. It is shown that the shear elasticity of muscles is determined by two independent moduli. The dissipative properties are determined by the fourth-rank viscosity tensor, which also has two independent components. The propagation velocity and attenuation of shear waves in muscle depend on the relative orientation of three vectors: the wave vector, the polarization vector, and the direction of muscle fiber. For one of the many experiments where attention was distinctly focused on the vector character of the wave process, it was possible to make a comparison with the theory, estimate the elasticity moduli, and obtain agreement with the angular dependence of the wave propagation velocity predicted by the theory.
Piezoresistive Sensor with High Elasticity Based on 3D Hybrid Network of Sponge@CNTs@Ag NPs.
Zhang, Hui; Liu, Nishuang; Shi, Yuling; Liu, Weijie; Yue, Yang; Wang, Siliang; Ma, Yanan; Wen, Li; Li, Luying; Long, Fei; Zou, Zhengguang; Gao, Yihua
2016-08-31
Pressure sensors with high elasticity are in great demand for the realization of intelligent sensing, but there is a need to develope a simple, inexpensive, and scalable method for the manufacture of the sensors. Here, we reported an efficient, simple, facile, and repeatable "dipping and coating" process to manufacture a piezoresistive sensor with high elasticity, based on homogeneous 3D hybrid network of carbon nanotubes@silver nanoparticles (CNTs@Ag NPs) anchored on a skeleton sponge. Highly elastic, sensitive, and wearable sensors are obtained using the porous structure of sponge and the synergy effect of CNTs/Ag NPs. Our sensor was also tested for over 2000 compression-release cycles, exhibiting excellent elasticity and cycling stability. Sensors with high performance and a simple fabrication process are promising devices for commercial production in various electronic devices, for example, sport performance monitoring and man-machine interfaces. PMID:27482721
Bulk solitary waves in elastic solids
NASA Astrophysics Data System (ADS)
Samsonov, A. M.; Dreiden, G. V.; Semenova, I. V.; Shvartz, A. G.
2015-10-01
A short and object oriented conspectus of bulk solitary wave theory, numerical simulations and real experiments in condensed matter is given. Upon a brief description of the soliton history and development we focus on bulk solitary waves of strain, also known as waves of density and, sometimes, as elastic and/or acoustic solitons. We consider the problem of nonlinear bulk wave generation and detection in basic structural elements, rods, plates and shells, that are exhaustively studied and widely used in physics and engineering. However, it is mostly valid for linear elasticity, whereas dynamic nonlinear theory of these elements is still far from being completed. In order to show how the nonlinear waves can be used in various applications, we studied the solitary elastic wave propagation along lengthy wave guides, and remarkably small attenuation of elastic solitons was proven in physical experiments. Both theory and generation for strain soliton in a shell, however, remained unsolved problems until recently, and we consider in more details the nonlinear bulk wave propagation in a shell. We studied an axially symmetric deformation of an infinite nonlinearly elastic cylindrical shell without torsion. The problem for bulk longitudinal waves is shown to be reducible to the one equation, if a relation between transversal displacement and the longitudinal strain is found. It is found that both the 1+1D and even the 1+2D problems for long travelling waves in nonlinear solids can be reduced to the Weierstrass equation for elliptic functions, which provide the solitary wave solutions as appropriate limits. We show that the accuracy in the boundary conditions on free lateral surfaces is of crucial importance for solution, derive the only equation for longitudinal nonlinear strain wave and show, that the equation has, amongst others, a bidirectional solitary wave solution, which lead us to successful physical experiments. We observed first the compression solitary wave in the
NASA Astrophysics Data System (ADS)
Lee, H.; Min, D.; Lim, S.; Yang, J.; Kwon, B.; Yoo, H.
2009-12-01
In a conventional marine seismic data analysis, pressure data have been usually interpreted on the basis of acoustic wave equation. The acoustic wave equation, however, only deals with P-wave propagation, and it cannot correctly describe the wave propagation in acoustic-elastic (fluid-solid) coupled media. Recently, in 4C OBC survey (4-component ocean bottom cable), it is possible to acquire both pressure and 3-component displacements (measured at the sea-bottom). Combining pressure and displacement data allows us to interpret subsurface structures more accurately. In order to accurately simulate wave propagation in fluid-solid coupled media, we need an acoustic-elastic coupled modeling algorithm, which deals with displacements in elastic region and pressure in acoustic region. For waveform inversion and reverse-time migration that require a great number of forward modeling, it is essential to develop an efficient scheme that reduces computing time and computer core memory. In this study, we present a 3D time-domain acoustic-elastic coupled modeling algorithm on the basis of the cell-based finite difference method. The cell-based method has proven to properly describe the free-surface boundary, which indicates that it will also properly describe the fluid-solid interface boundaries. In the acoustic-elastic coupled modeling, we first compose cell-based finite differences individually for the 3D acoustic and elastic media, and then combine the differences using the fluid-solid interface boundary conditions. Considering that the 2D acoustic-elastic coupled modeling algorithm gives numerical solutions comparable to analytic solutions, we expect that the 3D acoustic-elastic coupled modeling will correctly describe wave propagation in the fluid-solid coupled media. We apply our algorithm to 3D horizontal two- and three-layer models. Numerical experiments show that the cell-based coupled modeling algorithm properly describes S- and converted waves as well as P-waves. The
Wave-induced fluid flow in random porous media: attenuation and dispersion of elastic waves.
Müller, Tobias M; Gurevich, Boris
2005-05-01
A detailed analysis of the relationship between elastic waves in inhomogeneous, porous media and the effect of wave-induced fluid flow is presented. Based on the results of the poroelastic first-order statistical smoothing approximation applied to Biot's equations of poroelasticity, a model for elastic wave attenuation and dispersion due to wave-induced fluid flow in 3-D randomly inhomogeneous poroelastic media is developed. Attenuation and dispersion depend on linear combinations of the spatial correlations of the fluctuating poroelastic parameters. The observed frequency dependence is typical for a relaxation phenomenon. Further, the analytic properties of attenuation and dispersion are analyzed. It is shown that the low-frequency asymptote of the attenuation coefficient of a plane compressional wave is proportional to the square of frequency. At high frequencies the attenuation coefficient becomes proportional to the square root of frequency. A comparison with the 1-D theory shows that attenuation is of the same order but slightly larger in 3-D random media. Several modeling choices of the approach including the effect of cross correlations between fluid and solid phase properties are demonstrated. The potential application of the results to real porous materials is discussed. PMID:15957744
NASA Astrophysics Data System (ADS)
Chen, Wei-Qiu
2015-10-01
Significant progress has been made in mixed boundary-value problems associated with three-dimensional (3D) crack and contact analyses of advanced materials featuring more complexities compared to the conventional isotropic elastic materials. These include material anisotropy and multifield coupling, two typical characteristics of most current multifunctional materials. In this paper we try to present a state-of-the-art description of 3D exact/analytical solutions derived for crack and contact problems of elastic solids with both transverse isotropy and multifield coupling in the latest decade by the potential theory method in the spirit of V. I. Fabrikant, whose ingenious breakthrough brings new vigor and vitality to the old research subject of classical potential theory. We are particularly interested in crack and contact problems with certain nonlinear features. Emphasis is also placed on the coupling between the temperature field (or the like) and other physical fields (e.g., elastic, electric, and magnetic fields). We further highlight the practical significance of 3D contact solutions, in particular in applications related to modern scanning probe microscopes.
NASA Astrophysics Data System (ADS)
Petrov, P.; Newman, G. A.
2014-12-01
An application of the 3D elastic full-waveform inversion (FWI) to wide-aperture seismic data obtained for a complex geological setting is presented. Imaging is implemented in the Fourier domain, exploiting damped wave fields. The modeling part is solved with a finite-difference method. The non-linear conjugate gradient method is used for the inverse problem solution. The nonlinearity of FWI leads to the presence of local and multiple minima in the least-squares error functional especially for large offset problems. That leads to the shutdown of the inverse problem convergence and uncertainty in the solution. An accurate starting velocity model can avoid this problem, but in many cases may not be available. Hence other strategies are necessary to address the problem. We propose a robust inversion process for an arbitrary starting velocity model, which allows avoiding local minima and obtaining acceptable images of the deep seated structures defined by large offset data. We proceed from the assumption that decreasing data offset reduces local minima problems but decreases the depth of the recovered image. So, the inversion process is realized sequentially from small to large offsets, allowing recovery of geological structures over the entire depth range of interest from the near surface to deeper depths sensed only by large aperture offsets. Increasing of data offset is first performed at the lowest frequency and then proceeding with treatment of all data offsets from low to high frequencies. A reverse loop is also implemented in the laddering of frequencies, where after the inversion at high frequencies and all offsets we return to the lower frequencies data to continue the IP. Returning to lower frequency data provides helping to ameliorate multiple minima encountered in the inversion. The inversion then concludes by sweeping over higher frequency data, at all offsets. We demonstrate our strategies for treating wide aperture offset data on the Marmousi model, using
Field structure of collapsing wave packets in 3D strong Langmuir turbulence
NASA Technical Reports Server (NTRS)
Newman, D. L.; Robinson, P. A.; Goldman, M. V.
1989-01-01
A simple model is constructed for the electric fields in the collapsing wave packets found in 3D simulations of driven and damped isotropic strong Langmuir turbulence. This model, based on a spherical-harmonic decomposition of the electrostatic potential, accounts for the distribution of wave-packet shapes observed in the simulations, particularly the predominance of oblate wave packets. In contrast with predictions for undamped and undriven subsonic collapse of scalar fields, oblate vector-field wave packets do not flatten during collapse but, instead, remain approximately self-similar and rigid.
Bubbles attenuate elastic waves at seismic frequencies
NASA Astrophysics Data System (ADS)
Tisato, Nicola; Quintal, Beatriz; Chapman, Samuel; Podladchikov, Yury; Burg, Jean-Pierre
2016-04-01
The vertical migration of multiphase fluids in the crust can cause hazardous events such as eruptions, explosions, pollution and earthquakes. Although seismic tomography could potentially provide a detailed image of such fluid-saturated regions, the interpretation of the tomographic signals is often controversial and fails in providing a conclusive map of the subsurface saturation. Seismic tomography should be improved considering seismic wave attenuation (1/Q) and the dispersive elastic moduli which allow accounting for the energy lost by the propagating elastic wave. In particular, in saturated media a significant portion of the energy carried by the propagating wave is dissipated by the wave-induced-fluid-flow and the wave-induced-gas-exsolution-dissolution (WIGED) mechanisms. The WIGED mechanism describes how a propagating wave modifies the thermodynamic equillibrium between different fluid phases causing the exsolution and the dissolution of the gas in the liquid, which in turn causes a significant frequency dependent 1/Q and moduli dispersion. The WIGED theory was initially postulated for bubbly magmas but only recently was extended to bubbly water and experimentally demonstrated. Here we report these theory and laboratory experiments. Specifically, we present i) attenuation measurements performed by means of the Broad Band Attenuation Vessel on porous media saturated with water and different gases, and ii) numerical experiments validating the laboratory observations. Finally, we will extend the theory to fluids and to pressure-temperature conditions which are typical of phreatomagmatic and hydrocarbon domains and we will compare the propagation of seismic waves in bubble-free and bubble-bearing subsurface domains. With the present contribution we extend the knowledge about attenuation in rocks which are saturated with multiphase fluid demonstrating that the WIGED mechanism could be extremely important to image subsurface gas plumes.
3D extension of Tensorial Polar Decomposition. Application to (photo-)elasticity tensors
NASA Astrophysics Data System (ADS)
Desmorat, Rodrigue; Desmorat, Boris
2016-06-01
The orthogonalized harmonic decomposition of symmetric fourth-order tensors (i.e. having major and minor indicial symmetries, such as elasticity tensors) is completed by a representation of harmonic fourth-order tensors H by means of two second-order harmonic (symmetric deviatoric) tensors only. A similar decomposition is obtained for non-symmetric tensors (i.e. having minor indicial symmetry only, such as photo-elasticity tensors or elasto-plasticity tangent operators) introducing a fourth-order major antisymmetric traceless tensor Z. The tensor Z is represented by means of one harmonic second-order tensor and one antisymmetric second-order tensor only. Representations of totally symmetric (rari-constant), symmetric and major antisymmetric fourth-order tensors are simple particular cases of the proposed general representation. Closed-form expressions for tensor decomposition are given in the monoclinic case. Practical applications to elasticity and photo-elasticity monoclinic tensors are finally presented. xml:lang="fr"
3D WKB solution for fast magnetoacoustic wave behaviour around an X-line
NASA Astrophysics Data System (ADS)
McLaughlin, J. A.; Botha, G. J. J.; Régnier, S.; Spoors, D. L.
2016-06-01
Context. We study the propagation of a fast magnetoacoustic wave in a 3D magnetic field created from two magnetic dipoles. The magnetic topology contains an X-line. Aims: We aim to contribute to the overall understanding of MHD wave propagation within inhomogeneous media, specifically around X-lines. Methods: We investigate the linearised, 3D MHD equations under the assumptions of ideal and cold plasma. We utilise the WKB approximation and Charpit's method during our investigation. Results: It is found that the behaviour of the fast magnetoacoustic wave is entirely dictated by the local, inhomogeneous, equilibrium Alfvén speed profile. All parts of the wave experience refraction during propagation, where the magnitude of the refraction effect depends on the location of an individual wave element within the inhomogeneous magnetic field. The X-line, along which the Alfvén speed is identically zero, acts as a focus for the refraction effect. There are two main types of wave behaviour: part of the wave is either trapped by the X-line or escapes the system, and there exists a critical starting region around the X-line that divides these two types of behaviour. For the set-up investigated, it is found that 15.5% of the fast wave energy is trapped by the X-line. Conclusions: We conclude that linear, β = 0 fast magnetoacoustic waves can accumulate along X-lines and thus these will be specific locations of fast wave energy deposition and thus preferential heating. The work here highlights the importance of understanding the magnetic topology of a system. We also demonstrate how the 3D WKB technique described in this paper can be applied to other magnetic configurations.
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
Elastic waves along a fracture intersection
NASA Astrophysics Data System (ADS)
Abell, Bradley Charles
Fractures and fracture networks play a significant role in the subsurface hydraulic connectivity within the Earth. While a significant amount of research has been performed on the seismic response of single fractures and sets of fractures, few studies have examined the effect of fracture intersections on elastic wave propagation. Intersections play a key role in the connectivity of a fracture network that ultimately affects the hydraulic integrity of a rock mass. In this dissertation two new types of coupled waves are examined that propagate along intersections. 1) A coupled wedge wave that propagates along a surface fracture with particle motion highly localized to the intersection of a fracture with a free surface, and 2) fracture intersection waves that propagate along the intersection between two orthogonal fractures. Theoretical formulations were derived to determine the particle motion and velocity of intersection waves. Vibrational modes calculated from the theoretical formulation match those predicted by group theory based on the symmetry of the problem. For the coupled wedge wave, two vibrational modes exist that range in velocity between the wedge wave and Rayleigh wave velocity and exhibit either wagging or breathing motion depending on the Poisson's ratio. For the intersection waves, the observed modes depend on the properties of the fractures forming the intersection. If both fractures have equal stiffness four modes exist, two with wagging and two with breathing motion. If the fractures have unequal stiffness, four modes also exist, but the motion depends on the Poisson's ratio. The velocity of intersection waves depends on the coupling or stiffness of the intersection and frequency of the signal. In general, the different modes travel with speeds between the wedge wave and bulk shear wave velocity. Laboratory experiments were performed on isotropic and anisotropic samples to verify the existence of these waves. For both waves, the observed signals
Fast Wave Trains Associated with Solar Eruptions: Insights from 3D Thermodynamic MHD Simulations
NASA Astrophysics Data System (ADS)
Downs, C.; Liu, W.; Torok, T.; Linker, J.; Mikic, Z.; Ofman, L.
2015-12-01
EUV imaging observations during the SDO/AIA era have provided new insights into a variety of wave phenomena occurring in the low solar corona. One example is the observation of quasi-periodic, fast-propagating wave trains that are associated with solar eruptions, including flares and CMEs. While there has been considerable progress in understanding such waves from both an observational and theoretical perspective, it remains a challenge to pin down their physical origin. In this work, we detail our results from a case-study 3D thermodynamic MHD simulation of a coronal mass ejection where quasi-periodic wave trains are generated during the simulated eruption. We find a direct correlation between the onset of non-steady reconnection in the flare current sheet and the generation of quasi-periodic wave train signatures when patchy, collimated downflows interact with the flare arcade. Via forward modeling of SDO/AIA observables, we explore how the appearance of the wave trains is affected by line-of-sight integration and the multi-thermal nature of the coronal medium. We also examine how the wave trains themselves are channeled by natural waveguides formed in 3D by the non-uniform background magnetic field. While the physical association of the reconnection dynamics to the generation of quasi-periodic wave trains appears to be a compelling result, unanswered questions posed from recent observations as well as future prospects will be discussed.
NASA Astrophysics Data System (ADS)
Montiel, F.; Squire, V. A.
2013-12-01
A new ocean wave/sea-ice interaction model is proposed that simulates how a directional wave spectrum evolves as it travels through a realistic marginal ice zone (MIZ), where wave/ice dynamics are entirely governed by coherent conservative wave scattering effects. Field experiments conducted by Wadhams et al. (1986) in the Greenland Sea generated important data on wave attenuation in the MIZ and, particularly, on whether the wave spectrum spreads directionally or collimates with distance from the ice edge. The data suggest that angular isotropy, arising from multiple scattering by ice floes, occurs close to the edge and thenceforth dominates wave propagation throughout the MIZ. Although several attempts have been made to replicate this finding theoretically, including by the use of numerical models, none have confronted this problem in a 3D MIZ with fully randomised floe distribution properties. We construct such a model by subdividing the discontinuous ice cover into adjacent infinite slabs of finite width parallel to the ice edge. Each slab contains an arbitrary (but finite) number of circular ice floes with randomly distributed properties. Ice floes are modeled as thin elastic plates with uniform thickness and finite draught. We consider a directional wave spectrum with harmonic time dependence incident on the MIZ from the open ocean, defined as a continuous superposition of plane waves traveling at different angles. The scattering problem within each slab is then solved using Graf's interaction theory for an arbitrary incident directional plane wave spectrum. Using an appropriate integral representation of the Hankel function of the first kind (see Cincotti et al., 1993), we map the outgoing circular wave field from each floe on the slab boundaries into a directional spectrum of plane waves, which characterizes the slab reflected and transmitted fields. Discretizing the angular spectrum, we can obtain a scattering matrix for each slab. Standard recursive
NASA Astrophysics Data System (ADS)
Parisi, Laura; Ferreira, Ana M. G.
2016-04-01
The surface wave full ray theory (FRT) is an efficient tool to calculate synthetic waveforms of surface waves. It combines the concept of local modes with exact ray tracing as a function of frequency, providing a more complete description of surface wave propagation than the widely used great circle approximation (GCA). The purpose of this study is to evaluate the ability of the FRT approach to model teleseismic long-period surface waveforms (T ˜ 45-150 s) in the context of current 3-D Earth models to empirically assess its validity domain and its scope for future studies in seismic tomography. To achieve this goal, we compute vertical and horizontal component fundamental mode synthetic Rayleigh waveforms using the FRT, which are compared with calculations using the highly accurate spectral element method. We use 13 global earth models including 3-D crustal and mantle structure, which are derived by successively varying the strength and lengthscale of heterogeneity in current tomographic models. For completeness, GCA waveforms are also compared with the spectral element method. We find that the FRT accurately predicts the phase and amplitude of long-period Rayleigh waves (T ˜ 45-150 s) for almost all the models considered, with errors in the modelling of the phase (amplitude) of Rayleigh waves being smaller than 5 per cent (10 per cent) in most cases. The largest errors in phase and amplitude are observed for T ˜ 45 s and for the three roughest earth models considered that exhibit shear wave anomalies of up to ˜20 per cent, which is much larger than in current global tomographic models. In addition, we find that overall the GCA does not predict Rayleigh wave amplitudes well, except for the longest wave periods (T ˜ 150 s) and the smoothest models considered. Although the GCA accurately predicts Rayleigh wave phase for current earth models such as S20RTS and S40RTS, FRT's phase errors are smaller, notably for the shortest wave periods considered (T ˜ 45 s and
Elastic waves in structurally chiral composites
Yang, Shiuhkuang.
1990-01-01
Elastic wave propagation through structurally chiral (handed) media was studied. The primary objectives are to construct structurally chiral composites and to characterize their properties. Structurally chiral composites are constructed by stacking identical uniaxial plates, whose consecutive symmetric axes describe either a right- or a left-handed spiral. A matrix representation method is used to solve the elastic wave propagation in such layered composites. Numerical computation of the plane wave reflection and transmission characteristics for chiral arrangements are compared with those for the non-chiral one. It is concluded that the co-polarized characteristics are unaffected by the structural chirality, while the cross-polarized reflected and transmitted fields are greatly influenced by it. Numerical modeling is also applied for the real samples. The polarization ellipse of the transmitted field of each sample is calculated. To verify the form chirality, four glass-reinforced chiral and non-chiral composite samples are made from helix tape, molded, debulked, and cured individually under identical temperature and pressure histories. The spiral composites are characterized using shear and longitudinal wave transducers in ultrasonic experiments. Both the material properties and the polarization ellipse of the transmitted field of each sample are measured. It is proved conclusively that left and right handedness in the microstructures of a material rotates the plane of polarization of a propagating shear wave in the opposite directions. Thus it is now possible to say that by reducing the length scale of the handed microstructures tone more appropriate to its propagating wavelength, a medium is obtained that gives rise to effects similar to optical radar and optical dichroism.
Elastic Wave Propagation and Generation in Seismology
NASA Astrophysics Data System (ADS)
Lees, Jonathan M.
The majority of mature seismologists of my generation were introduced to theoretical seismology via classic textbooks written in the early 1980s. Since this generation has matured and taken the mantle of teaching seismology to a new generation, several new books have been put forward as replacements, or alternatives, to the original classical texts. The target readers of the new texts range from beginner through intermediate to more advanced, although all have been attempts to improve upon what is now considered standard convention in quantitative seismology. To this plethora of choices we now have a new addition by Jose Pujol, titledElastic Wave Propagation and Generation in Seismology.
Topologically protected elastic waves in phononic metamaterials
Mousavi, S. Hossein; Khanikaev, Alexander B.; Wang, Zheng
2015-01-01
Surface waves in topological states of quantum matter exhibit unique protection from backscattering induced by disorders, making them ideal carriers for both classical and quantum information. Topological matters for electrons and photons are largely limited by the range of bulk properties, and the associated performance trade-offs. In contrast, phononic metamaterials provide access to a much wider range of material properties. Here we demonstrate numerically a phononic topological metamaterial in an elastic-wave analogue of the quantum spin Hall effect. A dual-scale phononic crystal slab is used to support two effective spins for phonons over a broad bandwidth, and strong spin–orbit coupling is realized by breaking spatial mirror symmetry. By preserving the spin polarization with an external load or spatial symmetry, phononic edge states are shown to be robust against scattering from discrete defects as well as disorders in the continuum, demonstrating topological protection for phonons in both static and time-dependent regimes. PMID:26530426
Topologically protected elastic waves in phononic metamaterials.
Mousavi, S Hossein; Khanikaev, Alexander B; Wang, Zheng
2015-01-01
Surface waves in topological states of quantum matter exhibit unique protection from backscattering induced by disorders, making them ideal carriers for both classical and quantum information. Topological matters for electrons and photons are largely limited by the range of bulk properties, and the associated performance trade-offs. In contrast, phononic metamaterials provide access to a much wider range of material properties. Here we demonstrate numerically a phononic topological metamaterial in an elastic-wave analogue of the quantum spin Hall effect. A dual-scale phononic crystal slab is used to support two effective spins for phonons over a broad bandwidth, and strong spin-orbit coupling is realized by breaking spatial mirror symmetry. By preserving the spin polarization with an external load or spatial symmetry, phononic edge states are shown to be robust against scattering from discrete defects as well as disorders in the continuum, demonstrating topological protection for phonons in both static and time-dependent regimes. PMID:26530426
NASA Astrophysics Data System (ADS)
Paroni, Roberto; Sili, Ali
2016-02-01
We first consider an elastic thin heterogeneous cylinder of radius of order ε: the interior of the cylinder is occupied by a stiff material (fiber) that is surrounded by a soft material (matrix). By assuming that the elasticity tensor of the fiber does not scale with ε and that of the matrix scales with ε2, we prove that the one dimensional model is a nonlocal system. We then consider a reference configuration domain filled out by periodically distributed rods similar to those described above. We prove that the homogenized model is a second order nonlocal problem. In particular, we show that the homogenization problem is directly connected to the 3D-1D dimensional reduction problem.
Force sensing using 3D displacement measurements in linear elastic bodies
NASA Astrophysics Data System (ADS)
Feng, Xinzeng; Hui, Chung-Yuen
2016-07-01
In cell traction microscopy, the mechanical forces exerted by a cell on its environment is usually determined from experimentally measured displacement by solving an inverse problem in elasticity. In this paper, an innovative numerical method is proposed which finds the "optimal" traction to the inverse problem. When sufficient regularization is applied, we demonstrate that the proposed method significantly improves the widely used approach using Green's functions. Motivated by real cell experiments, the equilibrium condition of a slowly migrating cell is imposed as a set of equality constraints on the unknown traction. Our validation benchmarks demonstrate that the numeric solution to the constrained inverse problem well recovers the actual traction when the optimal regularization parameter is used. The proposed method can thus be applied to study general force sensing problems, which utilize displacement measurements to sense inaccessible forces in linear elastic bodies with a priori constraints.
A 3-D elasticity theory based model for acoustic radiation from multilayered anisotropic plates.
Shen, C; Xin, F X; Lu, T J
2014-05-01
A theoretical model built upon three-dimensional elasticity theory is developed to investigate the acoustic radiation from multilayered anisotropic plates subjected to a harmonic point force excitation. Fourier transform technique and stationary phase method are combined to predict the far-field radiated sound pressure of one-side water immersed plate. Compared to equivalent single-layer plate models, the present model based on elasticity theory can differentiate radiated sound pressure between dry-side and wet-side excited cases, as well as discrepancies induced by different layer sequences for multilayered anisotropic plates. These results highlight the superiority of the present theoretical model especially for handling multilayered anisotropic structures. PMID:24815294
Force sensing using 3D displacement measurements in linear elastic bodies
NASA Astrophysics Data System (ADS)
Feng, Xinzeng; Hui, Chung-Yuen
2016-04-01
In cell traction microscopy, the mechanical forces exerted by a cell on its environment is usually determined from experimentally measured displacement by solving an inverse problem in elasticity. In this paper, an innovative numerical method is proposed which finds the "optimal" traction to the inverse problem. When sufficient regularization is applied, we demonstrate that the proposed method significantly improves the widely used approach using Green's functions. Motivated by real cell experiments, the equilibrium condition of a slowly migrating cell is imposed as a set of equality constraints on the unknown traction. Our validation benchmarks demonstrate that the numeric solution to the constrained inverse problem well recovers the actual traction when the optimal regularization parameter is used. The proposed method can thus be applied to study general force sensing problems, which utilize displacement measurements to sense inaccessible forces in linear elastic bodies with a priori constraints.
Approximate method for controlling solid elastic waves by transformation media
NASA Astrophysics Data System (ADS)
Hu, Jin; Chang, Zheng; Hu, Gengkai
2011-11-01
By idealizing a general mapping as a series of local affine ones, we derive approximately transformed material parameters necessary to control solid elastic waves within classical elasticity theory. The transformed elastic moduli are symmetric, and can be used with Navier's equation to manipulate elastic waves. It is shown numerically that the method can provide a powerful tool to control elastic waves in solids in case of high frequency or small material gradient. Potential applications can be anticipated in nondestructive testing, structure impact protection, petroleum exploration, and seismology.
A 3D Orthotropic Strain-Rate Dependent Elastic Damage Material Model.
English, Shawn Allen
2014-09-01
A three dimensional orthotropic elastic constitutive model with continuum damage and cohesive based fracture is implemented for a general polymer matrix composite lamina. The formulation assumes the possibility of distributed (continuum) damage followed b y localized damage. The current damage activation functions are simply partially interactive quadratic strain criteria . However, the code structure allows for changes in the functions without extraordinary effort. The material model formulation, implementation, characterization and use cases are presented.
Interface waves in almost incompressible elastic materials
NASA Astrophysics Data System (ADS)
Virta, Kristoffer; Kreiss, Gunilla
2015-12-01
We study the problem of two elastic half-planes in contact and the Stoneley interface wave that may exist at the interface between two different elastic materials, emphasis being put on the case when the half-planes are almost incompressible. We show that numerical simulations involving interface waves require an unexpectedly high number of grid points per wavelength as the materials become more incompressible. Let λ, μ, ρ and λ‧, μ‧, ρ‧ be the Lamé parameters and densities of the first and second half-plane, respectively. A theoretical study shows that if K is a real constant, λ‧ = Kλ, μ‧ = Kμ, ρ‧ = Kρ and μ → 0, then for an accurate solution the required number of grid points per wavelength scales as (μ / λ) - 1 / p, where p is the order of accuracy of the numerical method. This requirement becomes very restrictive close to the incompressible limit μ ≪ λ, especially for lower order methods i.e., a small p. The theoretical findings are supported by numerical experiments that illustrate the demanding resolution requirement as well as the superiority of higher order methods. The scaling is also seen to hold for a more general choice of Lamé parameters. Numerical experiments when one of the half-planes is a vacuum are also presented, where the higher resolution requirement is illustrated in a numerical solution of Lamb's problem.
Gravitational Wave Signals from 2D and 3D Core Collapse Supernova Explosions
NASA Astrophysics Data System (ADS)
Yakunin, Konstantin; Mezzacappa, Anthony; Marronetti, Pedro; Bruenn, Stephen; Hix, W. Raphael; Lentz, Eric J.; Messer, O. E. Bronson; Harris, J. Austin; Endeve, Eirik; Blondin, John
2016-03-01
We study two- and three-dimensional (2D and 3D) core-collapse supernovae (CCSN) using our first-principles CCSN simulations performed with the neutrino hydrodynamics code CHIMERA. The following physics is included: Newtonian hydrodynamics with a nuclear equation of state capable of describing matter in both NSE and non-NSE, MGFLD neutrino transport with realistic neutrino interactions, an effective GR gravitational potential, and a nuclear reaction network. Both our 2D and 3D models achieve explosion, which in turn enables us to determine their complete gravitational wave signals. In this talk, we present them, and we analyze the similarities and differences between the 2D and 3D signals.
NASA Astrophysics Data System (ADS)
Chakraborty, Bidisha; Heyde, Brecht; Alessandrini, Martino; D'hooge, Jan
2016-04-01
Image registration techniques using free-form deformation models have shown promising results for 3D myocardial strain estimation from ultrasound. However, the use of this technique has mostly been limited to research institutes due to the high computational demand, which is primarily due to the computational load of the regularization term ensuring spatially smooth cardiac strain estimates. Indeed, this term typically requires evaluating derivatives of the transformation field numerically in each voxel of the image during every iteration of the optimization process. In this paper, we replace this time-consuming step with a closed-form solution directly associated with the transformation field resulting in a speed up factor of ~10-60,000, for a typical 3D B-mode image of 2503 and 5003 voxels, depending upon the size and the parametrization of the transformation field. The performance of the numeric and the analytic solutions was contrasted by computing tracking and strain accuracy on two realistic synthetic 3D cardiac ultrasound sequences, mimicking two ischemic motion patterns. Mean and standard deviation of the displacement errors over the cardiac cycle for the numeric and analytic solutions were 0.68+/-0.40 mm and 0.75+/-0.43 mm respectively. Correlations for the radial, longitudinal and circumferential strain components at end-systole were 0.89, 0.83 and 0.95 versus 0.90, 0.88 and 0.92 for the numeric and analytic regularization respectively. The analytic solution matched the performance of the numeric solution as no statistically significant differences (p>0.05) were found when expressed in terms of bias or limits-of-agreement.
3D printed elastic honeycombs with graded density for tailorable energy absorption
NASA Astrophysics Data System (ADS)
Bates, Simon R. G.; Farrow, Ian R.; Trask, Richard S.
2016-04-01
This work describes the development and experimental analysis of hyperelastic honeycombs with graded densities, for the purpose of energy absorption. Hexagonal arrays are manufactured from thermoplastic polyurethane (TPU) via fused filament fabrication (FFF) 3D printing and the density graded by varying cell wall thickness though the structures. Manufactured samples are subject to static compression tests and their energy absorbing potential analysed via the formation of energy absorption diagrams. It is shown that by grading the density through the structure, the energy absorption profile of these structures can be manipulated such that a wide range of compression energies can be efficiently absorbed.
Filtering of elastic waves by opal-based hypersonic crystal.
Salasyuk, Alexey S; Scherbakov, Alexey V; Yakovlev, Dmitri R; Akimov, Andrey V; Kaplyanskii, Alexander A; Kaplan, Saveliy F; Grudinkin, Sergey A; Nashchekin, Alexey V; Pevtsov, Alexander B; Golubev, Valery G; Berstermann, Thorsten; Brüggemann, Christian; Bombeck, Michael; Bayer, Manfred
2010-04-14
We report experiments in which high quality silica opal films are used as three-dimensional hypersonic crystals in the 10 GHz range. Controlled sintering of these structures leads to well-defined elastic bonding between the submicrometer-sized silica spheres, due to which a band structure for elastic waves is formed. The sonic crystal properties are studied by injection of a broadband elastic wave packet with a femtosecond laser. Depending on the elastic bonding strength, the band structure separates long-living surface acoustic waves with frequencies in the complete band gap from bulk waves with band frequencies that propagate into the crystal leading to a fast decay. PMID:20232893
The small data solutions of general 3-D quasilinear wave equations. II
NASA Astrophysics Data System (ADS)
Ding, Bingbing; Witt, Ingo; Yin, Huicheng
2016-07-01
This paper is a continuation of the work in [8], where the authors established the global existence of smooth small data solutions to the general 3-D quasilinear wave equation ∑ i , j = 0 3 gij (u , ∂ u) ∂ij2 u = 0 when the weak null condition holds. In the present paper, we show that the smooth small data solutions of equation ∑ i , j = 0 3 gij (u , ∂ u) ∂ij2 u = 0 will blow up in finite time when the weak null condition does not hold and a generic nondegenerate condition on the initial data is satisfied, moreover, a precise blowup time is completely determined. Therefore, collecting the main results in this paper and [8], we have given a basically complete study on the blowup or global existence of small data solutions to the 3-D quasilinear wave equation ∑ i , j = 0 3 gij (u , ∂ u) ∂ij2 u = 0.
The Vajont disaster: a 3D numerical simulation for the slide and the waves
NASA Astrophysics Data System (ADS)
Rubino, Angelo; Androsov, Alexey; Vacondio, Renato; Zanchettin, Davide; Voltzinger, Naum
2016-04-01
A very high resolution O(5 m), 3D hydrostatic nonlinear numerical model was used to simulate the dynamics of both the slide and the surface waves produced during the Vajont disaster (north Italy, 1963), one of the major landslide-induced tsunamis ever documented. Different simulated wave phenomena like, e.g., maximum run-up on the opposite shore, maximum height, and water velocity were analyzed and compared with data available in literature, including the results of a fully 3D simulation obtained with a Smoothed Particle Hydrodynamic code. The difference between measured and simulated after-slide bathymetries was calculated and used in an attempt to quantify the relative magnitude and extension of rigid and fluid motion components during the event.
A Self-Consistent Beam Loaded Travelling Wave Accelerator Model for use in TRACE-3D
NASA Astrophysics Data System (ADS)
Lampel, M. C.
1997-05-01
An optics model of a constant gradient traveling wave (CGTW) accelerator structure has been implemented for TRACE-3D. TRACE-3D is an envelope code including space charge that is used to model bunched beams in magnetic transport systems and radio frequency (rf) accelerators when the effects of beam current might be significant. The new matrix model has been developed to allow incorporation of particle beam loading (current) effects on the accelerator gradient and the accelerator structure's beam focusing properties in a self-consistent manner. The beam loaded electric field for a CGTW accelerator structure is constant for only a particular design current (e.g., 0 current), otherwise it can be written as a function of accelerator attenuation and axial position along the structure. The variation of the electric field through the structure has been taken into account in the new model. CGTW structures differ substantially in focusing properties and beam loading properties from standing wave structures. Examples will be presented using the new TW model, propagating electron beams with different currents through the Stanford Linear Accelerator Center's 3 m structure. The results will be compared to the zero current TW structure model in TRANSPORT and the Tank model (a standing wave structure model) in TRACE-3D. A computer demonstration of the code with the new element will also be presented.
NASA Astrophysics Data System (ADS)
Lieske, Mike; Schlurmann, Torsten
2016-04-01
INTRODUCTION & MOTIVATION The design of structures in coastal and offshore areas and their maintenance are key components of coastal protection. Usually, assessments of processes and loads on coastal structures are derived from experiments with flow and wave parameters in separate physical models. However, Peregrin (1976) already points out that processes in natural shallow coastal waters flow and sea state processes do not occur separately, but influence each other nonlinearly. Kemp & Simons (1982) perform 2D laboratory tests and study the interactions between a turbulent flow and following waves. They highlight the significance of wave-induced changes in the current properties, especially in the mean flow profiles, and draw attention to turbulent fluctuations and bottom shear stresses. Kemp & Simons (1983) also study these processes and features with opposing waves. Studies on the wave-current interaction in three-dimensional space for a certain wave height, wave period and water depth were conducted by MacIver et al. (2006). The research focus is set on the investigation of long-crested waves on obliquely opposing and following currents in the new 3D wave-current basin. METHODOLOGY In a first step the flow analysis without waves is carried out and includes measurements of flow profiles in the sweet spot of the basin at predefined measurement positions. Five measuring points in the water column have been delineated in different water depths in order to obtain vertical flow profiles. For the characterization of the undisturbed flow properties in the basin, an uniformly distributed flow was generated in the wave basin. In the second step wave analysis without current, the unidirectional wave propagation and wave height were investigated for long-crested waves in intermediate wave conditions. In the sweet spot of the wave basin waves with three different wave directions, three wave periods and uniform wave steepness were examined. For evaluation, we applied a common
Sagar, Nitin; Khanna, Kunal; Sardesai, Varda S; Singh, Atul K; Temgire, Mayur; Kalita, Mridula Phukan; Kadam, Sachin S; Soni, Vivek P; Bhartiya, Deepa; Bellare, Jayesh R
2016-12-01
Bioactive 3D composites play an important role in advanced biomaterial design to provide molecular coupling and improve integrity with the cellular environment of the native bone. In the present study, a hybrid lyophilized polymer composite blend of anionic charged sodium salt of carboxymethyl chitin and gelatin (CMChNa-GEL) reinforced with nano-rod agglomerated hydroxyapatite (nHA) has been developed with enhanced biocompatibility and tunable elasticity. The scaffolds have an open, uniform and interconnected porous structure with an average pore diameter of 157±30μm and 89.47+0.03% with four dimensional X-ray. The aspect ratio of ellipsoidal pores decrease from 4.4 to 1.2 with increase in gelatin concentration; and from 2.14 to 1.93 with decrease in gelling temperature. The samples were resilient with elastic stain at 1.2MPa of stress also decreased from 0.33 to 0.23 with increase in gelatin concentration. The crosslinker HMDI (hexamethylene diisocyanate) yielded more resilient samples at 1.2MPa in comparison to glutaraldehyde. Increased crosslinking time from 2 to 4h in continuous compression cycle show no improvement in maximum elastic stain of 1.2MPa stress. This surface elasticity of the scaffold enables the capacity of these materials for adherent self renewal and cultivation of the NTERA-2 cL.D1 (NT2/D1), pluripotent embryonal carcinoma cell with biomechanical surface, as is shown here. Proliferation with MG-63, ALP activity and Alizarin red mineralization assay on optimized scaffold demonstrated ***p<0.001 between different time points thus showing its potential for bone healing. In pre-clinical study histological bone response of the scaffold construct displayed improved activity of bone regeneration in comparison to self healing of control groups (sham) up to week 07 after implantation in rabbit tibia critical-size defect. Therefore, this nHA-CMChNa-GEL scaffold composite exhibits inherent and efficient physicochemical, mechanical and biological
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.
NASA Astrophysics Data System (ADS)
Pathak, Himanshu; Singh, Akhilendra; Singh, Indra Vir
2016-06-01
In this work, finite element method (FEM) and element free Galerkin method (EFGM) are coupled for solving 3D crack domains subjected to cyclic thermal load of constant amplitude. Crack growth contours and fatigue life have been obtained for each of the considered numerical examples. Thermo-elastic problems are decoupled into thermal and elastic problems . Firstly, the unknown temperature field is obtained by solving heat conduction equation, then, it is used as the input load in the elastic problem to calculate the displacement and stress fields. The geometrical discontinuity across crack surface is modelled by extrinsically enriched EFGM and the remaining part of the domain is approximated by standard finite element method. At the crack interface, a ramp function based interpolation scheme has been implemented. This coupled approach combines the advantages of both EFGM and FEM. A linear successive crack increment approach is used to model crack growth. The growing crack surface is traced by level set function. Standard Paris law is used for life estimation of the three-dimensional crack models. Different cases of planar and non-planar crack problems have been solved and their results are compared with the results obtained using extended finite element method to check accuracy, efficiency and robustness of the coupled FE-EFG approach implemented in this study.
NASA Astrophysics Data System (ADS)
Belashov, Vasily
We study the formation, structure, stability and dynamics of the multidimensional soliton-like beam structures forming on the low-frequency branch of oscillation in the ionospheric and magnetospheric plasma for cases when beta=4pinT/B(2) <<1 and beta>1. In first case with the conditions omega
3D resolution tests of two-plane wave approach using synthetic seismograms
NASA Astrophysics Data System (ADS)
Ceylan, S.; Larmat, C. S.; Sandvol, E. A.
2012-12-01
Two-plane wave tomography (TPWT) is becoming a standard approach to obtain fundamental mode Rayleigh wave phase velocities for a variety of tectonic settings. A recent study by Ceylan et al. (2012) has applied this method to eastern Tibet, using data from INDEPTH-IV and Namche-Barwa seismic experiments. The TPWT assumes that distortion of wavefronts at each station can be expressed as the sum of two plane waves. However, there is currently no robust or complete resolution test for TPWT, to address its limitations such as wavefront healing. In this study, we test the capabilities of TPWT and resolution of INDEPTH-IV seismic experiment, by performing 3D resolution tests using synthetic seismograms. Utilizing SPECFEM3D software, we compute synthetic data sets resolving periods down to ~30 s. We implement a checkerboard upper mantle (for depths between 50 and 650 km) with variable cell sizes, superimposed to PREM as the background model. We then calculate fundamental mode surface wave phase velocities using TPWT for periods between 33-143 seconds, using synthetic seismograms computed from our three dimensional hypothetical model. Assuming a constant Poisson's ratio, we use partial derivatives from Saito (1988) to invert for shear wave velocities. We show that the combination of TPWT and Saito (1988) methods is capable of retrieving anomalies down to depths of ~200 km for Rayleigh waves. Below these depths, we observe evidence of both lateral and vertical smearing. We also find that the traditional method for estimating the resolution of TPWT consistently overestimates phase velocity resolutions. Love waves exhibit adequate resolution down to depths of ~100 km. At depths greater than 100 km, smearing is more evident in SH wave results than those of SV waves. Increased smearing of SH waves is most probably due to propagation characteristics and shallower sensitivity of Love waves. Our results imply that TPWT can be applied to Love waves, making future investigations of
Multiple-mode Lamb wave scattering simulations using 3D elastodynamic finite integration technique.
Leckey, Cara A C; Rogge, Matthew D; Miller, Corey A; Hinders, Mark K
2012-02-01
We have implemented three-dimensional (3D) elastodynamic finite integration technique (EFIT) simulations to model Lamb wave scattering for two flaw-types in an aircraft-grade aluminum plate, a rounded rectangle flat-bottom hole and a disbond of the same shape. The plate thickness and flaws explored in this work include frequency-thickness regions where several Lamb wave modes exist and sometimes overlap in phase and/or group velocity. For the case of the flat-bottom hole the depth was incrementally increased to explore progressive changes in multiple-mode Lamb wave scattering due to the damage. The flat-bottom hole simulation results have been compared to experimental data and are shown to provide key insight for this well-defined experimental case by explaining unexpected results in experimental waveforms. For the rounded rectangle disbond flaw, which would be difficult to implement experimentally, we found that Lamb wave behavior differed significantly from the flat-bottom hole flaw. Most of the literature in this field is restricted to low frequency-thickness regions due to difficulties in interpreting data when multiple modes exist. We found that benchmarked 3D EFIT simulations can yield an understanding of scattering behavior for these higher frequency-thickness regions and in cases that would be difficult to set up experimentally. Additionally, our results show that 2D simulations would not have been sufficient for modeling the complicated scattering that occurred. PMID:21908011
3-D P Wave Velocity Structure of Marmara Region Using Local Earthquake Tomography
NASA Astrophysics Data System (ADS)
Işık, S. E.; Gurbuz, C.
2014-12-01
The 3D P wave velocity model of upper and lower crust of the Marmara Region between 40.200- 41.200N and 26.500- 30.500E is obtained by tomographic inversion (Simulps) of 47034 P wave arrivals of local earthquakes recorded at 90 land stations between October 2009 and December 2012 and 30 OBO stations and 14162 shot arrivals recorded at 35 OBO stations (Seismarmara Survey, 2001). We first obtained a 1D minimum model with Velest code in order to obtain an initial model for 3D inversion with 648 well located earthquakes located within the study area. After several 3D inversion trials we decided to create a more adequate initial model for 3D inversion. Choosing the initial model we estimated the 3D P wave velocity model representing the whole region both for land and sea. The results are tested by making Checkerboard , Restoring Resolution and Characteristic Tests, and the reliable areas of the resulting model is defined in terms of RDE, DWS, SF and Hit count distributions. By taking cross sections from the resulting model we observed the vertical velocity change along profiles crossing both land and sea. All the profiles crossing the basins showed that the high velocities of lower crust make extensions towards the basin area which looks like the force that gives a shape to the basins. These extensions of lower crust towards the basins appeared with an average velocity of 6.3 km/s which might be the result of the deformation due the shearing in the region. It is also interpreted that the development of these high velocities coincide with the development of the basins. Thus, both the basins and the high velocity zones around them might be resulted from the entrance of the NAF into the Marmara Sea and at the same time a shear regime was dominated due to the resistance of the northern Marmara Region (Yılmaz, 2010). The seismicity is observed between 5 km and 15 km after the 3D location of the earthquakes. The locations of the earthquakes improved and the seismogenic zone
Study of nonlinear 3-D evolution of kinetic Alfvén wave and fluctuation spectra
NASA Astrophysics Data System (ADS)
Sharma, Prachi; Yadav, Nitin; Sharma, R. P.
2015-11-01
Waves and instabilities play a very crucial role in astrophysical plasmas e.g. solar wind, Geospace etc. The main objective of current study is to investigate the importance of nonlinear processes associated with kinetic Alfvén waves (KAWs) in order to understand the physical mechanism behind the magnetopause turbulence. Numerical simulation of the coupled equations guiding the dynamics of three dimensionally propagating kinetic Alfvén wave (KAW) and slow magnetosonic wave has been performed for intermediate beta plasma (i.e. me/mi ≪ β < 1, where β is thermal to magnetic pressure ratio) applicable to the magnetopause. A simplified semi-analytical model based on paraxial approach has also been developed. We have examined the field localization and associated power spectrum of 3-D kinetic Alfvén wave for this nonlinear interaction. Governing dynamical equations of KAW and slow magnetosonic wave get coupled when the ponderomotive force arising due to pump KAW is taken into account while studying the slow magnetosonic wave dynamics. The numerical prediction of power law scaling is just consistent with the observation of THEMIS spacecraft in the magnetopause.
NASA Astrophysics Data System (ADS)
Wilder, F. D.; Ergun, R.; Goodrich, K.; Malaspina, D.; Eriksson, S.; Stawarz, J. E.; Sturner, A. P.; Holmes, J.; Burch, J. L.; Torbert, R. B.; Phan, T.; Le Contel, O.; Goldman, M. V.; Newman, D. L.; Lindqvist, P. A.; Khotyaintsev, Y. V.; Strangeway, R. J.; Russell, C. T.; Giles, B. L.; Pollock, C. J.
2015-12-01
The phenomenon of magnetic reconnection, especially at electron scales, is still poorly understood. One process that warrants further investigation is the role of wave phenomenon in mediating magnetic reconnection. Previous observations have shown the presence of electrostatic solitary waves (ESWs) as well as whistler mode waves near the dayside reconnection site. Additionally, recent simulations have suggested that whistler waves might be generated by electron phase space holes associated with ESWs as they propagate along the magnetic separatrix towards the diffusion region. Other observations have shown ESWs with distinct speeds and time scales, suggesting that different instabilities generate the ESWs. NASA's recently launched Magnetospheric Multiscale (MMS) mission presents a unique opportunity to investigate the roles of wave phenomena, such as ESWs and whistlers, in asymmetric reconnection at the dayside magnetopause. We will present 3-D electric and magnetic field data from magnetopause crossings by MMS during its first dayside science phase. Burst mode wave data and electron distributions from all four spacecraft will be analyzed to investigate the origin of these wave phenomena, as well as their impact on the reconnection electric field.
Numerical investigation of wave attenuation by vegetation using a 3D RANS model
NASA Astrophysics Data System (ADS)
Marsooli, Reza; Wu, Weiming
2014-12-01
Vegetation has been recognized as an important natural shoreline protection against storm surges and waves. Understanding of wave-vegetation interaction is essential for assessing the ability of vegetation patches, such as wetlands, to mitigate storm damages. In this study the wave attenuation by vegetation is investigated numerically using a 3-D model which solves the Reynolds-Averaged Navier-Stokes equations (RANS) by means of a finite-volume method based on collocated hexahedron mesh. A mixing length model is used for turbulence closure of the RANS equations. The water surface boundary is tracked using the Volume-of-Fluid (VOF) method with the Compressive Interface Capturing Scheme for Arbitrary Meshes (CICSAM) to solve the VOF advection equation. The presence of vegetation is taken into account by adding the vegetation drag and inertia forces to the momentum equations. The model is validated by several laboratory experiments of short wave propagation through vegetation over flat and sloping beds. The comparisons show good agreement between the measured data and calculated results, but the swaying motion of flexible vegetation which is neglected in this study can influence the accuracy of the wave height predictions. The model is then applied to one of the validation tests with different vegetation properties, revealing that the wave height attenuation by vegetation depends not only on the wave conditions, but also the vegetation characteristics such as vegetation height and density.
Support minimized inversion of acoustic and elastic wave scattering
Safaeinili, A.
1994-04-24
This report discusses the following topics on support minimized inversion of acoustic and elastic wave scattering: Minimum support inversion; forward modelling of elastodynamic wave scattering; minimum support linearized acoustic inversion; support minimized nonlinear acoustic inversion without absolute phase; and support minimized nonlinear elastic inversion.
Simulating Seismic Wave Propagation in 3-D Structure: A Case Study For Istanbul City
NASA Astrophysics Data System (ADS)
Yelkenci, Seda; Aktar, Mustafa
2013-04-01
Investigation of the wave propagation around the Marmara Sea, in particular for the city of Istanbul is critical because this target area is identified as one of the megacities with the highest seismic risk in the world. This study makes an attempt for creating an integrated 3D seismic/geologic model and precise understanding of 3-D wave propagation in the city of Istanbul. The approach is based on generating synthetic seismograms using realistic velocity structures as well as accurate location, focal mechanism and source parameters of reference earthquakes. The modarate size reference earthquakes occured in the Marmara Sea and were recorded by the National Seismic Network of Turkey as well as the network of Istanbul Early Warning and Rapid Response System. The seismograms are simulated by means of a 3-D finite difference method operated on parallel processing environment. In the content of creating a robust velocity model; 1D velocity models which are derived fom previous crustal studies of Marmara region such as refraction seismic and receiver functions have been conducted firstly for depths greater than 1km. Velocity structure in shallower part of the study region is then derived from recent geophysical and geotechnical surveys. To construct 3-D model from the obtained 1-D model data, a variety of interpolation methods are considered. According to the observations on amplitude and arrival time based on comparison of simulated seismograms, the considered velocity model is refined the way that S delay times are compensated. Another important task of this work is an application of the finite difference method to estimate three-dimensional seismic responses for a specified basin structure including soft sediments with low shear velocities in respect of the surrounded area in the Asian part of Istanbul. The analysis performed both in the time and frequency domain, helps in understanding of the comprehensive wave propagation characteristics and the distribution of
Characteristics of Oceanic Waves Caused by Landslides: Insights from 3D-Hydrocode Modeling
NASA Astrophysics Data System (ADS)
Elbeshausen, D.; Wünnemann, K.; Weiss, R.
2008-12-01
In an experimental framework, the generation of tsunami waves can be considered as a two-dimensional or three-dimensional problem by finding the respective geometry of the experimental set up. In nature, of course, it is a fully three-dimensional problem. The generation of tsunami waves caused by landslides (submarine and subaerial) must be approached as a dynamical problem. The understanding of the slide body's dynamics plays a key role in understanding the generated waves. Numerical calculations are a standard tool in tsunami science as the propagation of long waves can be tackled with depth-averaged equations. These classical models have often been used for modeling the propagation and run up of those tsunami waves caused by earthquakes. Tsunamis generated during slide motion are different. They are shorter and have larger amplitudes. It could be demonstrated in respective laboratory experiments and two dimensional numerical studies that in the initial phase the behavior of the waves is very complex, resulting in wave breaking and plunging. We conducted hydrocode simulations to model the generation of tsunami by slide in two dimensions. Results could show the complexity of the initial wave evolution as well as the development of the slide body itself. As an extension to these two-dimensional simulations, we now consider the three-dimensional problem and reveal some differences to the two-dimensional results. For this purpose we are using iSALE-3D, a multi- material, multi-rheology hydrocode capable of studying landslide processes in both two and three dimensions. iSALE-3D has been originally developed to study shock waves and high pressure scenarios like meteorite impacts or explosions and has been successfully validated against theoretical and experimental results as well as other numerical codes. Previous studies revealed differences in the formation and propagation of oceanic waves caused by meteorite impacts at different angles of incidence. Hence, for
3D P-Wave Velocity Structure of the Deep Galicia Rifted Margin
NASA Astrophysics Data System (ADS)
Bayrakci, Gaye; Minshull, Timothy; Davy, Richard; Sawyer, Dale; Klaeschen, Dirk; Papenberg, Cord; Reston, Timothy; Shillington, Donna; Ranero, Cesar
2015-04-01
The combined wide-angle reflection-refraction and multi-channel seismic (MCS) experiment, Galicia 3D, was carried out in 2013 at the Galicia rifted margin in the northeast Atlantic Ocean, west of Spain. The main geological features within the 64 by 20 km (1280 km²) 3D box investigated by the survey are the peridotite ridge (PR), the fault bounded, rotated basement blocks and the S reflector, which has been interpreted to be a low angle detachment fault. 44 short period four-component ocean bottom seismometers and 28 ocean bottom hydrophones were deployed in the 3D box. 3D MCS profiles sampling the whole box were acquired with two airgun arrays of 3300 cu.in. fired alternately every 37.5 m. We present the results from 3D first-arrival time tomography that constrains the P-wave velocity in the 3D box, for the entire depth sampled by reflection data. Results are validated by synthetic tests and by the comparison with Galicia 3D MCS lines. The main outcomes are as follows: 1- The 3.5 km/s iso-velocity contour mimics the top of the acoustic basement observed on MCS profiles. Block bounding faults are imaged as velocity contrasts and basement blocks exhibit 3D topographic variations. 2- On the southern profiles, the top of the PR rises up to 5.5 km depth whereas, 20 km northward, its basement expression (at 6.5 km depth) nearly disappears. 3- The 6.5 km/s iso-velocity contour matches the topography of the S reflector where the latter is visible on MCS profiles. Within a depth interval of 0.6 km (in average), velocities beneath the S reflector increase from 6.5 km/s to 7 km/s, which would correspond to a decrease in the degree of serpentinization from ~45 % to ~30 % if these velocity variations are caused solely by variations in hydration. At the intersections between the block bounding normal faults and the S reflector, this decrease happens over a larger depth interval (> 1 km), suggesting that faults act as conduit for the water flow in the upper mantle.
Stratification effects on nonlinear elastic surface waves
NASA Astrophysics Data System (ADS)
Parker, D. F.
1988-01-01
On a homogeneous elastic half-space, linear surface waves are nondispersive. In each direction, waves having any profile travel without distortion. Nonlinearity causes intermodulation between the various wavelengths so that the signal distorts. Even when nonlinearity is small, sinusoidal profiles do not remain approximately sinusoidal. The absence of dispersion means that profiles suffer cumulative distortion, until the surface slope and strain become locally unbounded. Although this behaviour is typical of many signals, there are some signals for which intermodulation is constructive. These signals can travel coherently over large distances. For seismological applications, it is important to study the effects due to stratification. Dependence of the material constants on depth modifies the nonlinear evolution equations previously derived for homogeneous media. It has a smaller effect on higher frequencies than on lower frequencies. An approximate theory for short wavelength (high frequency) signals is introduced. Calculations show that when nonlinearity is no more important than dispersion, initially sinusoidal profiles propagate with surface slope remaining finite. When dispersion is small compared to nonlinearity, certain sharp peaked profiles can travel large distances while suffering little distortion.
Concealed threat detection with the IRAD sub-millimeter wave 3D imaging radar
NASA Astrophysics Data System (ADS)
Robertson, Duncan A.; Cassidy, Scott L.; Jones, Ben; Clark, Anthony
2014-06-01
Sub-millimeter wave 3D imaging radar is a promising technology for the stand-off detection of threats concealed on people. The IRAD 340 GHz 3D imaging radar uses polarization intensity information to identify signatures associated with concealed threats. We report on an extensive trials program which has been carried out involving dozens of individual subjects wearing a variety of different clothing to evaluate the detection of a wide range of threat and benign items. We have developed an automatic algorithm to run on the radar which yields a level of anomaly indication in real time. Statistical analysis of the large volume of recorded data has enabled performance metrics for the radar system to be evaluated.
A support-operator method for viscoelastic wave modelling in 3-D heterogeneous media
NASA Astrophysics Data System (ADS)
Ely, Geoffrey P.; Day, Steven M.; Minster, Jean-Bernard
2008-01-01
We apply the method of support operators (SOM) to solve the 3-D, viscoelastic equations of motion for use in earthquake simulations. SOM is a generalized finite-difference method that can utilize meshes of arbitrary structure and incorporate irregular geometry. Our implementation uses a 3-D, logically rectangular, hexahedral mesh. Calculations are second-order in space and time. A correction term is employed for suppression of spurious zero-energy modes (hourglass oscillations). We develop a free surface boundary condition, and an absorbing boundary condition using the method of perfectly matched layers (PML). Numerical tests using a layered material model in a highly deformed mesh show good agreement with the frequency-wavenumber method, for resolutions greater than 10 nodes per wavelength. We also test a vertically incident P wave on a semi-circular canyon, for which results match boundary integral solutions at resolutions greater that 20 nodes per wavelength. We also demonstrate excellent parallel scalability of our code.
An efficient flexible-order model for 3D nonlinear water waves
Engsig-Karup, A.P. Bingham, H.B.; Lindberg, O.
2009-04-01
The flexible-order, finite difference based fully nonlinear potential flow model described in [H.B. Bingham, H. Zhang, On the accuracy of finite difference solutions for nonlinear water waves, J. Eng. Math. 58 (2007) 211-228] is extended to three dimensions (3D). In order to obtain an optimal scaling of the solution effort multigrid is employed to precondition a GMRES iterative solution of the discretized Laplace problem. A robust multigrid method based on Gauss-Seidel smoothing is found to require special treatment of the boundary conditions along solid boundaries, and in particular on the sea bottom. A new discretization scheme using one layer of grid points outside the fluid domain is presented and shown to provide convergent solutions over the full physical and discrete parameter space of interest. Linear analysis of the fundamental properties of the scheme with respect to accuracy, robustness and energy conservation are presented together with demonstrations of grid independent iteration count and optimal scaling of the solution effort. Calculations are made for 3D nonlinear wave problems for steep nonlinear waves and a shoaling problem which show good agreement with experimental measurements and other calculations from the literature.
3D Numerical Simulation on the Sloshing Waves Excited by the Seismic Shacking
NASA Astrophysics Data System (ADS)
Zhang, Lin; Wu, Tso-Ren
2016-04-01
In the event of 2015 Nepal earthquake, a video clip broadcasted worldwide showed a violent water spilling in a hotel swimming pool. This sloshing phenomenon indicates a potential water loss in the sensitive facilities, e.g. the spent fuel pools in nuclear power plant, has to be taken into account carefully under the consideration of seismic-induced ground acceleration. In the previous studies, the simulation of sloshing mainly focused on the pressure force on the structure by using a simplified Spring-Mass Method developed in the field of solid mechanics. However, restricted by the assumptions of plane water surface and limited wave height, significant error will be made in evaluating the amount of water loss in the tank. In this paper, the computational fluid dynamical model, Splash3D, was adopted for studying the sloshing problem accurately. Splash3D solved 3D Navier-Stokes Equation directly with Large-Eddy Simulation (LES) turbulent closure. The Volume-of-fluid (VOF) method with piecewise linear interface calculation (PLIC) was used to track the complex breaking water surface. The time series acceleration of a design seismic was loaded to excite the water. With few restrictions from the assumptions, the accuracy of the simulation results were improved dramatically. A series model validations were conducted by compared to a 2D theoretical solution, and a 3D experimental data. Good comparisons can be seen. After the validation, we performed the simulation for considering a sloshing case in a rectangular water tank with a dimension of 12 m long, 8 m wide, 8 m deep, which contained water with 7 m in depth. The seismic movement was imported by considering time-series acceleration in three dimensions, which were about 0.5 g to 1.2 g in the horizontal directions, and 0.3 g to 1 g in the vertical direction. We focused the discussions on the kinematics of the water surface, wave breaking, velocity field, pressure field, water force on the side walls, and, most
3D frequency-domain finite-difference modeling of acoustic wave propagation
NASA Astrophysics Data System (ADS)
Operto, S.; Virieux, J.
2006-12-01
We present a 3D frequency-domain finite-difference method for acoustic wave propagation modeling. This method is developed as a tool to perform 3D frequency-domain full-waveform inversion of wide-angle seismic data. For wide-angle data, frequency-domain full-waveform inversion can be applied only to few discrete frequencies to develop reliable velocity model. Frequency-domain finite-difference (FD) modeling of wave propagation requires resolution of a huge sparse system of linear equations. If this system can be solved with a direct method, solutions for multiple sources can be computed efficiently once the underlying matrix has been factorized. The drawback of the direct method is the memory requirement resulting from the fill-in of the matrix during factorization. We assess in this study whether representative problems can be addressed in 3D geometry with such approach. We start from the velocity-stress formulation of the 3D acoustic wave equation. The spatial derivatives are discretized with second-order accurate staggered-grid stencil on different coordinate systems such that the axis span over as many directions as possible. Once the discrete equations were developed on each coordinate system, the particle velocity fields are eliminated from the first-order hyperbolic system (following the so-called parsimonious staggered-grid method) leading to second-order elliptic wave equations in pressure. The second-order wave equations discretized on each coordinate system are combined linearly to mitigate the numerical anisotropy. Secondly, grid dispersion is minimized by replacing the mass term at the collocation point by its weighted averaging over all the grid points of the stencil. Use of second-order accurate staggered- grid stencil allows to reduce the bandwidth of the matrix to be factorized. The final stencil incorporates 27 points. Absorbing conditions are PML. The system is solved using the parallel direct solver MUMPS developed for distributed
Mach-wave coherence in 3D media with random heterogeneities
NASA Astrophysics Data System (ADS)
Vyas, Jagdish C.; Mai, P. Martin; Galis, Martin; Dunham, Eric M.; Imperatori, Walter
2016-04-01
We investigate Mach-waves coherence for complex super-shear ruptures embedded in 3D random media that lead to seismic scattering. We simulate Mach-wave using kinematic earthquake sources that include fault-regions over which the rupture propagates at super-shear speed. The local slip rate is modeled with the regularized Yoffe function. The medium heterogeneities are characterized by Von Karman correlation function. We consider various realizations of 3D random media from combinations of different values of correlation length (0.5 km, 2 km, 5 km), standard deviation (5%, 10%, 15%) and Hurst exponent (0.2). Simulations in a homogeneous medium serve as a reference case. The ground-motion simulations (maximum resolved frequency of 5 Hz) are conducted by solving the elasto-dynamic equations of motions using a generalized finite-difference method, assuming a vertical strike-slip fault. The seismic wavefield is sampled at numerous locations within the Mach-cone region to study the properties and evolution of the Mach-waves in scattering media. We find that the medium scattering from random heterogeneities significantly diminishes the coherence of Mach-wave in terms of both amplitude and frequencies. We observe that Mach-waves are considerably scattered at distances RJB > 20 km (and beyond) for random media with standard deviation 10%. The scattering efficiency of the medium for small Hurst exponents (H <= 0.2) is mainly controlled by the standard deviation of the velocity heterogeneities, rather than their correlation length, as both theoretical considerations and numerical experiments indicate. Based on our simulations, we propose that local super-shear ruptures may be more common in nature then reported, but are very difficult to detect due to the strong seismic scattering. We suggest that if an earthquake is recorded within 10-15 km fault perpendicular distance and has high PGA, then inversion should be carried out by allowing rupture speed variations from sub
NASA Astrophysics Data System (ADS)
Woodbury, D.; Kubota, S.; Johnson, I.
2014-10-01
Computer simulations of electromagnetic wave propagation in magnetized plasmas are an important tool for both plasma heating and diagnostics. For active millimeter-wave and microwave diagnostics, accurately modeling the evolution of the beam parameters for launched, reflected or scattered waves in a toroidal plasma requires that calculations be done using the full 3-D geometry. Previously, we reported on the application of GPGPU (General-Purpose computing on Graphics Processing Units) to a 3-D vacuum Maxwell code using the FDTD (Finite-Difference Time-Domain) method. Tests were done for Gaussian beam propagation with a hard source antenna, utilizing the parallel processing capabilities of the NVIDIA K20M. In the current study, we have modified the 3-D code to include a soft source antenna and an induced current density based on the cold plasma approximation. Results from Gaussian beam propagation in an inhomogeneous anisotropic plasma, along with comparisons to ray- and beam-tracing calculations will be presented. Additional enhancements, such as advanced coding techniques for improved speedup, will also be investigated. Supported by U.S. DoE Grant DE-FG02-99-ER54527 and in part by the U.S. DoE, Office of Science, WDTS under the Science Undergraduate Laboratory Internship program.
3D Simulation of an Audible Ultrasonic Electrolarynx Using Difference Waves
Mills, Patrick; Zara, Jason
2014-01-01
A total laryngectomy removes the vocal folds which are fundamental in forming voiced sounds that make speech possible. Although implanted prosthetics are commonly used in developed countries, simple handheld vibrating electrolarynxes are still common worldwide. These devices are easy to use but suffer from many drawbacks including dedication of a hand, mechanical sounding voice, and sound leakage. To address some of these drawbacks, we introduce a novel electrolarynx that uses vibro-acoustic interference of dual ultrasonic waves to generate an audible fundamental frequency. A 3D simulation of the principles of the device is presented in this paper. PMID:25401965
Joint inversion of 3D crustal structure with ambient noise and earthquake body wave travel time
NASA Astrophysics Data System (ADS)
Li, Z.; Ni, S.; Chong, J.; Wang, X.
2012-12-01
Surface wave tomography based on the noise correlation function of seismic ambient noise has been widely used in studies of crustal and mantle structure . However, the periods of surface wave dispersions in the ambient noise tomography are typically less than 40 s, which limits its resolution on the lower crust. Travel times of earthquake body waves, such as Sg and SmS, could provide additional constraints to the crustal structure, especially to the lower crust due to the ray paths of SmS traveling through the lower crust twice. Here, we proposed a joint inversion method for 3D crustal structure with ambient noise and earthquake body wave travel time data, with the goal of providing better constraints and resolutions on the whole crust. We constructed the linear equations for joint inversion of crustal S velocity structure with the surface wave dispersion and body wave travel time data, and solved the equations with LSQR algorithm. Different weighting and damping factors, together with smoothing constraints, are adopted for surface wave dispersion and body wave travel time data to fit both dataset simultaneously. Synthetics experiments showed that the joint inversion could resolve the crust structure better than sole tomography of ambient noise or body wave travel time. We conducted the joint inversion around the Yangtze block in the eastern China. Rayleigh wave dispersions are extracted from the seismic ambient noise tomography by Zheng et al (2011) in this area. The body waves (e.g., Sg, SmS, Sn) are coherent to be identified and their travel times are measured with accuracy from high quality waveforms of some recent local earthquakes in this area. In order to minimize the travel time uncertainties, the focal depth and epicenter of these local earthquakes were resolved by depth phases and temporary aftershock observations. The result from joint inversion suggests that the crustal velocity structure, especially the lower crust, was well improved, which not only
NASA Astrophysics Data System (ADS)
Tago, J.; Cruz-Atienza, V. M.; Etienne, V.; Virieux, J.; Benjemaa, M.; Sanchez-Sesma, F. J.
2010-12-01
Simulating any realistic seismic scenario requires incorporating physical basis into the model. Considering both the dynamics of the rupture process and the anelastic attenuation of seismic waves is essential to this purpose and, therefore, we choose to extend the hp-adaptive Discontinuous Galerkin finite-element method to integrate these physical aspects. The 3D elastodynamic equations in an unstructured tetrahedral mesh are solved with a second-order time marching approach in a high-performance computing environment. The first extension incorporates the viscoelastic rheology so that the intrinsic attenuation of the medium is considered in terms of frequency dependent quality factors (Q). On the other hand, the extension related to dynamic rupture is integrated through explicit boundary conditions over the crack surface. For this visco-elastodynamic formulation, we introduce an original discrete scheme that preserves the optimal code performance of the elastodynamic equations. A set of relaxation mechanisms describes the behavior of a generalized Maxwell body. We approximate almost constant Q in a wide frequency range by selecting both suitable relaxation frequencies and anelastic coefficients characterizing these mechanisms. In order to do so, we solve an optimization problem which is critical to minimize the amount of relaxation mechanisms. Two strategies are explored: 1) a least squares method and 2) a genetic algorithm (GA). We found that the improvement provided by the heuristic GA method is negligible. Both optimization strategies yield Q values within the 5% of the target constant Q mechanism. Anelastic functions (i.e. memory variables) are introduced to efficiently evaluate the time convolution terms involved in the constitutive equations and thus to minimize the computational cost. The incorporation of anelastic functions implies new terms with ordinary differential equations in the mathematical formulation. We solve these equations using the same order
3D rendering of passive millimeter-wave scenes using modified open source software
NASA Astrophysics Data System (ADS)
Murakowski, Maciej; Wilson, John; Murakowski, Janusz; Schneider, Garrett; Schuetz, Christopher; Prather, Dennis
2011-05-01
As millimeter-wave imaging technology becomes more mature, several applications are emerging for which this technology may be useful. However, effectively predicting the nuances of millimeter-wave phenomenology on the usefulness for a given application remains a challenge. To this end, an accurate millimeter-wave scene simulator would have tremendous value in predicting imager requirements for a given application. Herein, we present a passive millimeter-wave scene simulator built on the open-source 3d modeling software Blender. We describe the changes made to the Blender rendering engine to make it suitable for this purpose, including physically accurate reflections at each material interface, volumetric absorption and scattering, and tracking of both s and p polarizations. In addition, we have incorporated a mmW material database and world model that emulates the effects of cold sky profiles for varying weather conditions and frequencies of operation. The images produced by this model have been validated against calibrated experimental imagery captured by a passive scanning millimeter-wave imager for maritime, desert, and standoff detection applications.
Wave optics theory and 3-D deconvolution for the light field microscope
Broxton, Michael; Grosenick, Logan; Yang, Samuel; Cohen, Noy; Andalman, Aaron; Deisseroth, Karl; Levoy, Marc
2013-01-01
Light field microscopy is a new technique for high-speed volumetric imaging of weakly scattering or fluorescent specimens. It employs an array of microlenses to trade off spatial resolution against angular resolution, thereby allowing a 4-D light field to be captured using a single photographic exposure without the need for scanning. The recorded light field can then be used to computationally reconstruct a full volume. In this paper, we present an optical model for light field microscopy based on wave optics, instead of previously reported ray optics models. We also present a 3-D deconvolution method for light field microscopy that is able to reconstruct volumes at higher spatial resolution, and with better optical sectioning, than previously reported. To accomplish this, we take advantage of the dense spatio-angular sampling provided by a microlens array at axial positions away from the native object plane. This dense sampling permits us to decode aliasing present in the light field to reconstruct high-frequency information. We formulate our method as an inverse problem for reconstructing the 3-D volume, which we solve using a GPU-accelerated iterative algorithm. Theoretical limits on the depth-dependent lateral resolution of the reconstructed volumes are derived. We show that these limits are in good agreement with experimental results on a standard USAF 1951 resolution target. Finally, we present 3-D reconstructions of pollen grains that demonstrate the improvements in fidelity made possible by our method. PMID:24150383
Solving tolerancing and 3D beam shaping problems by multifunctional wave optical design
NASA Astrophysics Data System (ADS)
Buehling, Sven; Wyrowski, Frank
2000-10-01
A strategy for designing optical systems that are optimized for multiple optical functions on the basis of wave optics is presented. Each optical function is composed of an input field, a set of fixed system parameters, and a merit function. A design algorithm is proposed which is applicable for optical systems consisting of an transmission operator followed by an arbitrary linear operator. The goal is to find the transmission operator which is optimal for all optical functions simultaneously. In later design steps, the found transmission operator can be transformed to real optical elements, for instance by using the thin element approximation. It is shown that the algorithm is efficiently applicable by investigating two sample applications for multifunctional wave optical design: the design of tolerant systems and 3D beam shaping.
NASA Astrophysics Data System (ADS)
Crosta, G.; Imposimato, S.; Roddeman, D.; Frattini, P.
2012-04-01
Fast moving landslides can be originated along slopes in mountainous terrains with natural and artificial lakes, or fjords at the slope foot. This landslides can reach extremely high speed and the impact with the immobile reservoir water can be influenced by the local topography and the landslide mass profile. The impact can generate large impulse waves and landslide tsunami. Initiation, propagation and runup are the three phases that need to be considered. The landslide evolution and the consequent wave can be controlled by the initial mass position (subaerial, partially or completely submerged), the landslide speed, the type of material, the subaerial and subaqueous slope geometry, the landslide depth and length at the impact, and the water depth. Extreme events have been caused by subaerial landslides: the 1963 Vajont rockslide (Italy), the 1958 Lituya Bay event (Alaska), the Tafjord and the Loen multiple events event (Norway), also from volcanic collapses (Hawaii and Canary islands). Various researchers completed a systematic experimental work on 2D and 3D wave generation and propagation (Kamphuis and Bowering, 1970; Huber, 1980; Müller, 1995; Huber and Hager, 1997; Fritz, 2002; Zweifel, 2004; Panizzo et al., 2005; Heller, 2007; Heller and Kinnear, 2010; Sælevik et al., 2009), using both rigid blocks and deformable granular" masses. Model data and results have been used to calibrate and validate numerical modelling tools (Harbitz, 1992; Jiang and LeBlond, 1993; Grilli et al., 2002; Grilli and Watts, 2005; Lynett and Liu, 2005; Tinti et al., 2006; Abadie et al., 2010) generally considering simplified rheologies (e.g. viscous rheologies) for subaerial subaqueous spreading. We use a FEM code (Roddeman, 2011; Crosta et al., 2006, 2009, 2010, 2011) adopting an Eulerian-Lagrangian approach to give accurate results for large deformations. We model both 2D and fully 3D events considering different settings. The material is considered as a fully deformable elasto
Development of a GPU-Accelerated 3-D Full-Wave Code for Reflectometry Simulations
NASA Astrophysics Data System (ADS)
Reuther, K. S.; Kubota, S.; Feibush, E.; Johnson, I.
2013-10-01
1-D and 2-D full-wave codes used as synthetic diagnostics in microwave reflectometry are standard tools for understanding electron density fluctuations in fusion plasmas. The accuracy of the code depends on how well the wave properties along the ignored dimensions can be pre-specified or neglected. In a toroidal magnetic geometry, such assumptions are never strictly correct and ray tracing has shown that beam propagation is inherently a 3-D problem. Previously, we reported on the application of GPGPU's (General-Purpose computing on Graphics Processing Units) to a 2-D FDTD (Finite-Difference Time-Domain) code ported to utilize the parallel processing capabilities of the NVIDIA C870 and C1060. Here, we report on the development of a FDTD code for 3-D problems. Initial tests will use NVIDIA's M2070 GPU and concentrate on the launching and propagation of Gaussian beams in free space. If available, results using a plasma target will also be presented. Performance will be compared with previous generations of GPGPU cards as well as with NVIDIA's newest K20C GPU. Finally, the possibility of utilizing multiple GPGPU cards in a cluster environment or in a single node will also be discussed. Supported by U.S. DoE Grants DE-FG02-99-ER54527 and DE-AC02-09CH11466 and the DoE National Undergraduate Fusion Fellowship.
Rayleigh wave effects in an elastic half-space.
NASA Technical Reports Server (NTRS)
Aggarwal, H. R.
1972-01-01
Consideration of Rayleigh wave effects in a homogeneous isotropic linearly elastic half-space subject to an impulsive uniform disk pressure loading. An approximate formula is obtained for the Rayleigh wave effects. It is shown that the Rayleigh waves near the center of loading arise from the portion of the dilatational and shear waves moving toward the axis, after they originate at the edge of the load disk. A study is made of the vertical displacement due to Rayleigh waves at points on the axis near the surface of the elastic half-space.
NASA Astrophysics Data System (ADS)
Plotnikov, Illya; Vourlidas, Angelos; Tylka, Allan J.; Pinto, Rui; Rouillard, Alexis; Tirole, Margot
2016-07-01
Identifying the physical mechanisms that produce the most energetic particles is a long-standing observational and theoretical challenge in astrophysics. Strong pressure waves have been proposed as efficient accelerators both in the solar and astrophysical contexts via various mechanisms such as diffusive-shock/shock-drift acceleration and betatron effects. In diffusive-shock acceleration, the efficacy of the process relies on shock waves being super-critical or moving several times faster than the characteristic speed of the medium they propagate through (a high Alfven Mach number) and on the orientation of the magnetic field upstream of the shock front. High-cadence, multipoint imaging using the NASA STEREO, SOHO and SDO spacecrafts now permits the 3-D reconstruction of pressure waves formed during the eruption of coronal mass ejections. Using these unprecedented capabilities, some recent studies have provided new insights on the timing and longitudinal extent of solar energetic particles, including the first derivations of the time-dependent 3-dimensional distribution of the expansion speed and Mach numbers of coronal shock waves. We will review these recent developments by focusing on particle events that occurred between 2011 and 2015. These new techniques also provide the opportunity to investigate the enigmatic long-duration gamma ray events.
3D laboratory experiments on a system of low-crested breakwaters under oblique wave attack
NASA Astrophysics Data System (ADS)
Papacharalampous, Georgia; Karantinos, Michalis; Giantsi, Theodora; Moutzouris, Constantinos
2016-04-01
Low-crested breakwaters are being increasingly used for shore protection. Hydrodynamics around coastal structures are complicated and have not been fully understood. A series of large scale (1:40) 3D laboratory experiments were carried out in the Laboratory of Harbour Works, National Technical University of Athens to investigate the wave disturbance around a system of two non-parallel to the shoreline breakwaters. The structures were of the type of low-crested, permeable and attacked by obliquely incident waves. Three different water depths were tested in the basin with a range of various different spectra. The transmission and reflection coefficients were measured in the middle of each breakwater. For this purpose, 1 gauge and 4 gauges (in line) were placed on the landward and seaward side of each breakwater respectively. The effect of diffraction is incorporate at the measured wave heights. The measured coefficients are being compared to their corresponding estimated using existing empirical formulas. Most of those formulas neglect wave obliquity.
Active elastic metamaterials for subwavelength wave propagation control
NASA Astrophysics Data System (ADS)
Chen, Y. Y.; Huang, G. L.
2015-06-01
Recent research activities in elastic metamaterials demonstrate a significant potential for subwavelength wave propagation control owing to their interior locally resonant mechanism. The growing technological developments in electro/magnetomechanical couplings of smart materials have introduced a controlling degree of freedom for passive elastic metamaterials. Active elastic metamaterials could allow for a fine control of material physical behavior and thereby induce new functional properties that cannot be produced by passive approaches. In this paper, two types of active elastic metamaterials with shunted piezoelectric materials and electrorheological elastomers are proposed. Theoretical analyses and numerical validations of the active elastic metamaterials with detailed microstructures are presented for designing adaptive applications in band gap structures and extraordinary waveguides. The active elastic metamaterial could provide a new design methodology for adaptive wave filters, high signal-to-noise sensors, and structural health monitoring applications.
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.
Wilcox, Lucas C.; Stadler, Georg; Burstedde, Carsten; Ghattas, Omar
2010-12-10
We introduce a high-order discontinuous Galerkin (dG) scheme for the numerical solution of three-dimensional (3D) wave propagation problems in coupled elastic-acoustic media. A velocity-strain formulation is used, which allows for the solution of the acoustic and elastic wave equations within the same unified framework. Careful attention is directed at the derivation of a numerical flux that preserves high-order accuracy in the presence of material discontinuities, including elastic-acoustic interfaces. Explicit expressions for the 3D upwind numerical flux, derived as an exact solution for the relevant Riemann problem, are provided. The method supports h-non-conforming meshes, which are particularly effective at allowing local adaptation of the mesh size to resolve strong contrasts in the local wavelength, as well as dynamic adaptivity to track solution features. The use of high-order elements controls numerical dispersion, enabling propagation over many wave periods. We prove consistency and stability of the proposed dG scheme. To study the numerical accuracy and convergence of the proposed method, we compare against analytical solutions for wave propagation problems with interfaces, including Rayleigh, Lamb, Scholte, and Stoneley waves as well as plane waves impinging on an elastic-acoustic interface. Spectral rates of convergence are demonstrated for these problems, which include a non-conforming mesh case. Finally, we present scalability results for a parallel implementation of the proposed high-order dG scheme for large-scale seismic wave propagation in a simplified earth model, demonstrating high parallel efficiency for strong scaling to the full size of the Jaguar Cray XT5 supercomputer.
NASA Astrophysics Data System (ADS)
Wilcox, Lucas C.; Stadler, Georg; Burstedde, Carsten; Ghattas, Omar
2010-12-01
We introduce a high-order discontinuous Galerkin (dG) scheme for the numerical solution of three-dimensional (3D) wave propagation problems in coupled elastic-acoustic media. A velocity-strain formulation is used, which allows for the solution of the acoustic and elastic wave equations within the same unified framework. Careful attention is directed at the derivation of a numerical flux that preserves high-order accuracy in the presence of material discontinuities, including elastic-acoustic interfaces. Explicit expressions for the 3D upwind numerical flux, derived as an exact solution for the relevant Riemann problem, are provided. The method supports h-non-conforming meshes, which are particularly effective at allowing local adaptation of the mesh size to resolve strong contrasts in the local wavelength, as well as dynamic adaptivity to track solution features. The use of high-order elements controls numerical dispersion, enabling propagation over many wave periods. We prove consistency and stability of the proposed dG scheme. To study the numerical accuracy and convergence of the proposed method, we compare against analytical solutions for wave propagation problems with interfaces, including Rayleigh, Lamb, Scholte, and Stoneley waves as well as plane waves impinging on an elastic-acoustic interface. Spectral rates of convergence are demonstrated for these problems, which include a non-conforming mesh case. Finally, we present scalability results for a parallel implementation of the proposed high-order dG scheme for large-scale seismic wave propagation in a simplified earth model, demonstrating high parallel efficiency for strong scaling to the full size of the Jaguar Cray XT5 supercomputer.
Rayleigh–Bloch waves along elastic diffraction gratings
Colquitt, D. J.; Craster, R. V.; Antonakakis, T.; Guenneau, S.
2015-01-01
Rayleigh–Bloch (RB) waves in elasticity, in contrast to those in scalar wave systems, appear to have had little attention. Despite the importance of RB waves in applications, their connections to trapped modes and the ubiquitous nature of diffraction gratings, there has been no investigation of whether such waves occur within elastic diffraction gratings for the in-plane vector elastic system. We identify boundary conditions that support such waves and numerical simulations confirm their presence. An asymptotic technique is also developed to generate effective medium homogenized equations for the grating that allows us to replace the detailed microstructure by a continuum representation. Further numerical simulations confirm that the asymptotic scheme captures the essential features of these waves. PMID:25568616
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.
3D surface-wave tomography in the central Baltic Shield
NASA Astrophysics Data System (ADS)
Bruneton, M.; Pedersen, H. A.; Farra, V.; Svekalapko Seismic Tomography Working Group
2003-04-01
The main objective of the SVEKALAPKO deep seismic experiment was to image in details the lithosphere-asthenosphere system of the central Baltic Shield, therefore enhancing our knowledge of the structure and evolution of cratonic lithosphere. During the experiment a regular 2D grid of 46 broad-band stations covered the southern part of Finland. This exceptional stations distribution made it possible to undertake a high precision surface-wave tomography. We developed a technique based on paraxial ray tracing to obtain 2D phase-velocity maps as a function of frequency which can subsequently be inverted for the 3D structure. The major improvement is that we jointly inverted for the velocity model under the array and the shape of incoming wave fronts, therefore reducing artifacts due to structure outside the study region. The data set included carefully selected fundamental mode Rayleigh wave arrival times of 69 teleseismic events, computed using Wiener filtering. An average dispersion curve was obtained imposing the phase-velocity to be quasi constant. It leads to shear-wave velocities for the lithospheric mantle 4% faster than standard Earth model ak135. The inversion of the same data set was also conducted using weaker constraints to obtain the lateral variations of the phase-velocity at each frequency and subsequently of the shear-wave velocity as a function of depth. Three Vs profiles were computed respectively in the Karelian Archean province, in the Proterozoic Svekofennia, and at the suture between the two domains. They showed significant variations, the higher lithospheric velocities were seen in the proterozoic domain, a low velocity zone was necessary only in the suture zone. Our results showed that chemical changes are maintained within the lithosphere over extended periods of time.
System-in-package LTCC platform for 3D RF to millimeter wave
NASA Astrophysics Data System (ADS)
Vähä-Heikkilä, T.; Lahti, M.
2011-04-01
This presentation shows recent trends and results in 3D Low Temperature Co-Fired Ceramics (LTCC) modules in applications from RF to millimeter waves. The system-in-package LTCC platform is a true three dimensional module technology. LTCC is a lightweight multi-layer technology having typically 6-20 ceramic layers and metallizations between. The metallization levels i.e different metal layers can be patterned and connected together with metal vias. Passive devices can also be fabricated on LTCC while active devices and other chips are connected with flip-chip, wire bonding or soldering. In addition to passives directly fabricated to LTCC, several different technologies/ chips can be hybrid integrated to the same module. LTCC platform is also well suited for the realization of antenna arrays for microwave and millimeter wave applications. Potential applications are ranging from short range communications to space and radars. VTT has designed, fabricated and characterized microwave and millimeter wave packages for Radio Frequency (RF) Micro Electro Mechanical Systems (MEMS) as well as active devices. Also, several types of system-in-package modules have been realized containing hybrid integrated CMOS and GaAs MMICs and antenna arrays.
Energy in elastic fiber embedded in elastic matrix containing incident SH wave
NASA Technical Reports Server (NTRS)
Williams, James H., Jr.; Nagem, Raymond J.
1989-01-01
A single elastic fiber embedded in an infinite elastic matrix is considered. An incident plane SH wave is assumed in the infinite matrix, and an expression is derived for the total energy in the fiber due to the incident SH wave. A nondimensional form of the fiber energy is plotted as a function of the nondimensional wavenumber of the SH wave. It is shown that the fiber energy attains maximum values at specific values of the wavenumber of the incident wave. The results obtained here are interpreted in the context of phenomena observed in acousto-ultrasonic experiments on fiber reinforced composite materials.
Full 3D dispersion curve solutions for guided waves in generally anisotropic media
NASA Astrophysics Data System (ADS)
Hernando Quintanilla, F.; Lowe, M. J. S.; Craster, R. V.
2016-02-01
Dispersion curves of guided waves provide valuable information about the physical and elastic properties of waves propagating within a given waveguide structure. Algorithms to accurately compute these curves are an essential tool for engineers working in non-destructive evaluation and for scientists studying wave phenomena. Dispersion curves are typically computed for low or zero attenuation and presented in two or three dimensional plots. The former do not always provide a clear and complete picture of the dispersion loci and the latter are very difficult to obtain when high values of attenuation are involved and arbitrary anisotropy is considered in single or multi-layered systems. As a consequence, drawing correct and reliable conclusions is a challenging task in the modern applications that often utilize multi-layered anisotropic viscoelastic materials. These challenges are overcome here by using a spectral collocation method (SCM) to robustly find dispersion curves in the most complicated cases of high attenuation and arbitrary anisotropy. Solutions are then plotted in three-dimensional frequency-complex wavenumber space, thus gaining much deeper insight into the nature of these problems. The cases studied range from classical examples, which validate this approach, to new ones involving materials up to the most general triclinic class for both flat and cylindrical geometry in multi-layered systems. The apparent crossing of modes within the same symmetry family in viscoelastic media is also explained and clarified by the results. Finally, the consequences of the centre of symmetry, present in every crystal class, on the solutions are discussed.
Utilization of 3-D elastic transformation in the registration of chest x-ray CT and whole body PET
Tai, Yuan-Chuan; Hoh, C.K.; Hoffman, E.J.
1996-12-31
X-ray CT is widely used for detection and localization of lesions in the thorax. Whole Body PET with 18-FDG is becoming accepted for staging of cancer because of its ability to detect malignancy. Combining information from these two modalities has a significant value to improve lung cancer staging and treatment planning. Due to the non-rigid nature of the thorax and the differences in the acquisition conventions, the subject is stretched non-uniformly and the images of these two modalities requires non-rigid transformation for proper registration. Techniques to register chest x-ray CT and Whole Body PET images were developed and evaluated. Accuracy of 3-D elastic transformation was tested by phantom study. Studies on patients with lung carcinoma were used to validate the technique in localizing the 18-FDG uptake and in correlating PET to x-ray CT images. The fused images showed an accurate alignment and provided confident identification of the detailed anatomy of the CT with the functional information of the PET images.
3D P and S Wave Velocity Structure and Tremor Locations in the Parkfield Region
NASA Astrophysics Data System (ADS)
Zeng, X.; Thurber, C. H.; Shelly, D. R.; Bennington, N. L.; Cochran, E. S.; Harrington, R. M.
2014-12-01
We have assembled a new dataset to refine the 3D seismic velocity model in the Parkfield region. The S arrivals from 184 earthquakes recorded by the Parkfield Experiment to Record MIcroseismicity and Tremor array (PERMIT) during 2010-2011 were picked by a new S wave picker, which is based on machine learning. 74 blasts have been assigned to four quarries, whose locations were identified with Google Earth. About 1000 P and S wave arrivals from these blasts at permanent seismic network were also incorporated. Low frequency earthquakes (LFEs) occurring within non-volcanic tremor (NVT) are valuable for improving the precision of NVT location and the seismic velocity model at greater depths. Based on previous work (Shelley and Hardebeck, 2010), waveforms of hundreds of LFEs in same family were stacked to improve signal qualify. In a previous study (McClement et al., 2013), stacked traces of more than 30 LFE families at the Parkfileld Array Seismic Observatory (PASO) have been picked. We expanded our work to include LFEs recorded by the PERMIT array. The time-frequency Phase Weight Stacking (tf-PWS) method was introduced to improve the stack quality, as direct stacking does not produce clear S-wave arrivals on the PERMIT stations. This technique uses the coherence of the instantaneous phase among the stacked signals to enhance the signal-to-noise ratio (SNR) of the stack. We found that it is extremely effective for picking LFE arrivals (Thurber et al., 2014). More than 500 P and about 1000 S arrivals from 58 LFE families were picked at the PERMIT and PASO arrays. Since the depths of LFEs are much deeper than earthquakes, we are able to extend model resolution to lower crustal depths. Both P and S wave velocity structure have been obtained with the tomoDD method. The result suggests that there is a low velocity zone (LVZ) in the lower crust and the location of the LVZ is consistent with the high conductivity zone beneath the southern segment of the Rinconada fault that
NASA Astrophysics Data System (ADS)
Bayrakci, Gaye; Minshull, Timothy A.; Davy, Richard G.; Karplus, Marianne S.; Kaeschen, Dirk; Papenberg, Cord; Krabbenhoeft, Anne; Sawyer, Dale; Reston, Timothy J.; Shillington, Donna J.; Ranero, César R.
2014-05-01
Galicia 3D, a reflection-refraction and long offset seismic experiment was carried out from May through September 2013, at the Galicia rifted margin (in the northeast Atlantic Ocean, west of Spain) as a collaboration between US, UK, German and Spanish groups. The 3D multichannel seismic acquisition conducted by R/V Marcus Langseth covered a 64 km by 20 km (1280 km2) zone where the main geological features are the Peridotite Ridge (PR), composed of serpentinized peridotite and thought be upper mantle exhumed to the seafloor during rifting, and the S reflector which has been interpreted to be a low angle detachment fault overlain by fault bounded, rotated, continental crustal blocks. In the 3D box, two airgun arrays of 3300 cu.in. were fired alternately (in flip-flop configuration) every 37.5 m. All shots are recorded by 44 short period four component ocean bottom seismometers (OBS) and 26 ocean bottom hydrophones (OBH) deployed and recovered by R/V Poseidon, as well as four 6 km hydrophone streamers with 12.5 m channel spacing towed by R/V Marcus Langseth. We present the preliminary results of the first arrival time tomography study which is carried out with a subset of the wide-angle dataset, in order to generate a 3D P-wave velocity volume for the entire depth sampled by the reflection data. After the relocation of OBSs and OBHs, an automatic first-arrival time picking approach is applied to a subset of the dataset, which comprises more than 5.5 million source-receiver pairs. Then, the first-arrival times are checked visually, in 3-dimensions. The a priori model used for the first-arrival time tomography is built up using information from previous seismic surveys carried out at the Galicia margin (e.g. ISE, 1997). The FAST algorithm of Zelt and Barton (1998) is used for the first-arrival time inversion. The 3D P-wave velocity volume can be used in interpreting the reflection dataset, as a starting point for migration, to quantify the thinning of the crustal layers
Hydrodynamic analysis of elastic floating collars in random waves
NASA Astrophysics Data System (ADS)
Bai, Xiao-dong; Zhao, Yun-peng; Dong, Guo-hai; Li, Yu-cheng
2015-06-01
As the main load-bearing component of fish cages, the floating collar supports the whole cage and undergoes large deformations. In this paper, a mathematical method is developed to study the motions and elastic deformations of elastic floating collars in random waves. The irregular wave is simulated by the random phase method and the statistical approach and Fourier transfer are applied to analyze the elastic response in both time and frequency domains. The governing equations of motions are established by Newton's second law, and the governing equations of deformations are obtained based on curved beam theory and modal superposition method. In order to validate the numerical model of the floating collar attacked by random waves, a series of physical model tests are conducted. Good relationship between numerical simulation and experimental observations is obtained. The numerical results indicate that the transfer function of out-of-plane and in-plane deformations increase with the increasing of wave frequency. In the frequency range between 0.6 Hz and 1.1 Hz, a linear relationship exists between the wave elevations and the deformations. The average phase difference between the wave elevation and out-of-plane deformation is 60° with waves leading and the phase between the wave elevation and in-plane deformation is 10° with waves lagging. In addition, the effect of fish net on the elastic response is analyzed. The results suggest that the deformation of the floating collar with fish net is a little larger than that without net.
Intersymbol Interference Investigations Using a 3D Time-Dependent Traveling Wave Tube Model
NASA Technical Reports Server (NTRS)
Kory, Carol L.; Andro, Monty; Downey, Alan (Technical Monitor)
2001-01-01
For the first time, a physics based computational model has been used to provide a direct description of the effects of the TWT (Traveling Wave Tube) on modulated digital signals. The TWT model comprehensively takes into account the effects of frequency dependent AM/AM and AM/PM conversion; gain and phase ripple; drive-induced oscillations; harmonic generation; intermodulation products; and backward waves. Thus, signal integrity can be investigated in the presence of these sources of potential distortion as a function of the physical geometry of the high power amplifier and the operational digital signal. This method promises superior predictive fidelity compared to methods using TWT models based on swept amplitude and/or swept frequency data. The fully three-dimensional (3D), time-dependent, TWT interaction model using the electromagnetic code MAFIA is presented. This model is used to investigate assumptions made in TWT black box models used in communication system level simulations. In addition, digital signal performance, including intersymbol interference (ISI), is compared using direct data input into the MAFIA model and using the system level analysis tool, SPW (Signal Processing Worksystem).
Intersymbol Interference Investigations Using a 3D Time-Dependent Traveling Wave Tube Model
NASA Technical Reports Server (NTRS)
Kory, Carol L.; Andro, Monty
2002-01-01
For the first time, a time-dependent, physics-based computational model has been used to provide a direct description of the effects of the traveling wave tube amplifier (TWTA) on modulated digital signals. The TWT model comprehensively takes into account the effects of frequency dependent AM/AM and AM/PM conversion; gain and phase ripple; drive-induced oscillations; harmonic generation; intermodulation products; and backward waves. Thus, signal integrity can be investigated in the presence of these sources of potential distortion as a function of the physical geometry and operating characteristics of the high power amplifier and the operational digital signal. This method promises superior predictive fidelity compared to methods using TWT models based on swept- amplitude and/or swept-frequency data. First, the TWT model using the three dimensional (3D) electromagnetic code MAFIA is presented. Then, this comprehensive model is used to investigate approximations made in conventional TWT black-box models used in communication system level simulations. To quantitatively demonstrate the effects these approximations have on digital signal performance predictions, including intersymbol interference (ISI), the MAFIA results are compared to the system level analysis tool, Signal Processing Workstation (SPW), using high order modulation schemes including 16 and 64-QAM.
Modeling Recent Large Earthquakes Using the 3-D Global Wave Field
NASA Astrophysics Data System (ADS)
Hjörleifsdóttir, V.; Kanamori, H.; Tromp, J.
2003-04-01
We use the spectral-element method (SEM) to accurately compute waveforms at periods of 40 s and longer for three recent large earthquakes using 3D Earth models and finite source models. The M_w~7.6, Jan~26, 2001, Bhuj, India event had a small rupture area and is well modeled at long periods with a point source. We use this event as a calibration event to investigate the effects of 3-D Earth models on the waveforms. The M_w~7.9, Nov~11, 2001, Kunlun, China, event exhibits a large directivity (an asymmetry in the radiation pattern) even at periods longer than 200~s. We used the source time function determined by Kikuchi and Yamanaka (2001) and the overall pattern of slip distribution determined by Lin et al. to guide the wave-form modeling. The large directivity is consistent with a long fault, at least 300 km, and an average rupture speed of 3±0.3~km/s. The directivity at long periods is not sensitive to variations in the rupture speed along strike as long as the average rupture speed is constant. Thus, local variations in rupture speed cannot be ruled out. The rupture speed is a key parameter for estimating the fracture energy of earthquakes. The M_w~8.1, March~25, 1998, event near the Balleny Islands on the Antarctic Plate exhibits large directivity in long period surface waves, similar to the Kunlun event. Many slip models have been obtained from body waves for this earthquake (Kuge et al. (1999), Nettles et al. (1999), Antolik et al. (2000), Henry et al. (2000) and Tsuboi et al. (2000)). We used the slip model from Henry et al. to compute SEM waveforms for this event. The synthetic waveforms show a good fit to the data at periods from 40-200~s, but the amplitude and directivity at longer periods are significantly smaller than observed. Henry et al. suggest that this event comprised two subevents with one triggering the other at a distance of 100 km. To explain the observed directivity however, a significant amount of slip is required between the two subevents
Liu, Fengming; Liu, Zhengyou
2015-10-23
We theoretically investigate elastic waves propagating in metamaterials with simultaneous zero indices for both the longitudinal and transverse waves. With scattering objects (here cylinders) present in the metamaterial slabs, while the elastic waves can mostly transmit through the metamaterial slabs perfectly, exhibiting the well-known cloaking effect of zero-index metamaterials, they nevertheless become totally blocked at resonances, indicating strong elastic wave scattering by the objects in the cases. However, despite the occurrence of the elastic wave scattering, there is, counterintuitively, no mode conversion between the longitudinal and transverse waves in the process, completely in contrast to that in conventional elastic media. A design of a two-dimensional phononic crystal with these peculiar properties is presented. PMID:26551124
NASA Astrophysics Data System (ADS)
Miah, Khalid; Bellefleur, Gilles
2014-05-01
The global demand for base metals, uranium and precious metals has been pushing mineral explorations at greater depth. Seismic techniques and surveys have become essential in finding and extracting mineral rich ore bodies, especially for deep VMS mining camps. Geophysical parameters collected from borehole logs and laboratory measurements of core samples provide preliminary information about the nature and type of subsurface lithologic units. Alteration halos formed during the hydrothermal alteration process contain ore bodies, which are of primary interests among geologists and mining industries. It is known that the alteration halos are easier to detect than the ore bodies itself. Many 3D geological models are merely projection of 2D surface geology based on outcrop inspections and geochemical analysis of a small number of core samples collected from the area. Since a large scale 3D multicomponent seismic survey can be prohibitively expensive, performance analysis of such geological models can be helpful in reducing exploration costs. In this abstract, we discussed challenges and constraints encountered in geophysical modelling of ore bodies and surrounding geologic structures from the available coarse 3D geological models of the Lalor Lake mining camp, located in northern Manitoba, Canada. Ore bodies in the Lalor lake VMS camp are rich in gold, zinc, lead and copper, and have an approximate weight of 27 Mt. For better understanding of physical parameters of these known ore bodies and potentially unknown ones at greater depth, we constructed a fine resolution 3D seismic model with dimensions: 2000 m (width), 2000 m (height), and 1500 m (vertical depth). Seismic properties (P-wave, S-wave velocities, and density) were assigned based on a previous rock properties study of the same mining camp. 3D finite-difference elastic wave propagation simulation was performed in the model using appropriate parameters. The generated synthetic 3D seismic data was then compared to
Attenuation of Elastic Waves due to Scattering from Spherical Cavities and Elastic Inclusions.
NASA Astrophysics Data System (ADS)
Hinders, Mark Karl
1990-01-01
The attenuation of elastic waves due to scattering from a spherical inclusion of arbitrary size in an infinitely extended medium is investigated. The spherical scatterer and the exterior medium are isotropic, homogeneous, and linearly elastic, but of arbitrarily differing material parameters, with compressional and shear waves supported in both media. Exact expressions for scattered and transmitted fields caused by an incident plane compressional or shear wave of unit amplitude are calculated analytically and general expressions for extinction and scattering cross -sections are derived for both lossy and lossless scattering. Application to ultrasonic determination of porosity in cast aluminum is investigated.
Laser Excitation of a Fracture Source for Elastic Waves
Blum, Thomas E.; Wijk, Kasper van; Snieder, Roel; Willis, Mark E.
2011-12-30
We show that elastic waves can be excited at a fracture inside a transparent sample by focusing laser light directly onto this fracture. The associated displacement field, measured by a laser interferometer, has pronounced waves that are diffracted at the fracture tips. We confirm that these are tip diffractions from direct excitation of the fracture by comparing them with tip diffractions from scattered elastic waves excited on the exterior of the sample. Being able to investigate fractures - in this case in an optically transparent material - via direct excitation opens the door to more detailed studies of fracture properties in general.
NASA Astrophysics Data System (ADS)
Liu, Xiaolong; Eckert, Andreas; Connolly, Peter
2016-06-01
Buckle folds of sedimentary strata commonly feature a variety of different fracture sets. Some fracture sets including outer arc tensile fractures and inner arc shear fractures at the fold hinge zones are well understood by the extensional and compressional strain/stress pattern. However, other commonly observed fracture sets, including tensile fractures parallel to the fold axis, tensile fractures cutting through the limb, extensional faults at the fold hinge, and other shear fractures of various orientations in the fold limb, fail to be intuitively explained by the strain/stress regimes during the buckling process. To obtain a better understanding of the conditions for the initiation of the various fractures sets associated with single-layer cylindrical buckle folds, a 3D finite element modeling approach using a Maxwell visco-elastic rheology is utilized. The influences of three model parameters with significant influence on fracture initiation are considered: burial depth, viscosity, and permeability. It is concluded that these parameters are critical for the initiation of major fracture sets at the hinge zone with varying degrees. The numerical simulation results further show that the buckling process fails to explain most of the fracture sets occurring in the limb unless the process of erosional unloading as a post-fold phenomenon is considered. For fracture sets that only develop under unrealistic boundary conditions, the results demonstrate that their development is realistic for a perclinal fold geometry. In summary, a more thorough understanding of fractures sets associated with buckle folds is obtained based on the simulation of in-situ stress conditions during the structural development of buckle folds.
Quantitative analysis of accuracy of seismic wave-propagation codes in 3D random scattering media
NASA Astrophysics Data System (ADS)
Galis, Martin; Imperatori, Walter; Mai, P. Martin
2013-04-01
Several recent verification studies (e.g. Day et al., 2001; Bielak et al., 2010, Chaljub et al., 2010) have demonstrated the importance of assessing the accuracy of available numerical tools at low frequency in presence of large-scale features (basins, topography, etc.). The fast progress in high-performance computing, including efficient optimization of numerical codes on petascale supercomputers, has permitted the simulation of 3D seismic wave propagation at frequencies of engineering interest (up to 10Hz) in highly heterogeneous media (e.g. Hartzell et al, 2010; Imperatori and Mai, 2013). However, high frequency numerical simulations involving random scattering media, characterized by small-scale heterogeneities, are much more challenging for most numerical methods, and their verification may therefore be even more crucial than in the low-frequency case. Our goal is to quantitatively compare the accuracy and the behavior of three different numerical codes for seismic wave propagation in 3D random scattering media at high frequency. We deploy a point source with omega-squared spectrum, and focus on the near-source region, being of great interest in strong motion seismology. We use two codes based on finite-difference method (FD1 and FD2) and one code based on support-operator method (SO). Both FD1 and FD2 are 4-th order staggered-grid finite-difference codes (for FD1 see Olsen et al., 2009; for FD2 see Moczo et al., 2007). The FD1 and FD2 codes are characterized by slightly different medium representations, since FD1 uses point values of material parameters in each FD-cell, while FD2 uses the effective material parameters at each grid-point (Moczo et al., 2002). SO is 2-nd order support-operator method (Ely et al., 2008). We considered models with random velocity perturbations described by van Karman correlation function with different correlation lengths and different standard deviations. Our results show significant variability in both phase and amplitude as
Anomalously low amplitude of S waves produced by the 3D structures in the lower mantle
NASA Astrophysics Data System (ADS)
To, Akiko; Capdeville, Yann; Romanowicz, Barbara
2016-07-01
Direct S and Sdiff phases with anomalously low amplitudes are recorded for the earthquakes in Papua New Guinea by seismographs in northern America. According to the prediction by a standard 1D model, the amplitudes are the lowest at stations in southern California, at a distance and azimuth of around 95° and 55°, respectively, from the earthquake. The amplitude anomaly is more prominent at frequencies higher than 0.03 Hz. We checked and ruled out the possibility of the anomalies appearing because of the errors in the focal mechanism used in the reference synthetic waveform calculations. The observed anomaly distribution changes drastically with a relatively small shift in the location of the earthquake. The observations indicate that the amplitude reduction is likely due to the 3D shear velocity (Vs) structure, which deflects the wave energy away from the original ray paths. Moreover, some previous studies suggested that some of the S and Sdiff phases in our dataset are followed by a prominent postcursor and show a large travel time delay, which was explained by placing a large ultra-low velocity zone (ULVZ) located on the core-mantle boundary southwest of Hawaii. In this study, we evaluated the extent of amplitude anomalies that can be explained by the lower mantle structures in the existing models, including the previously proposed ULVZ. In addition, we modified and tested some models and searched for the possible causes of low amplitudes. Full 3D synthetic waveforms were calculated and compared with the observations. Our results show that while the existing models explain the trends of the observed amplitude anomalies, the size of such anomalies remain under-predicted especially at large distances. Adding a low velocity zone, which is spatially larger and has less Vs reduction than ULVZ, on the southwest side of ULVZ, contributes to explain the low amplitudes observed at distances larger than 100° from the earthquake. The newly proposed low velocity zone
Wave Phase-Sensitive Transformation of 3d-Straining of Mechanical Fields
NASA Astrophysics Data System (ADS)
Smirnov, I. N.; Speranskiy, A. A.
2015-11-01
It is the area of research of oscillatory processes in elastic mechanical systems. Technical result of innovation is creation of spectral set of multidimensional images which reflect time-correlated three-dimensional vector parameters of metrological, and\\or estimated, and\\or design parameters of oscillations in mechanical systems. Reconstructed images of different dimensionality integrated in various combinations depending on their objective function can be used as homeostatic profile or cybernetic image of oscillatory processes in mechanical systems for an objective estimation of current operational conditions in real time. The innovation can be widely used to enhance the efficiency of monitoring and research of oscillation processes in mechanical systems (objects) in construction, mechanical engineering, acoustics, etc. Concept method of vector vibrometry based on application of vector 3D phase- sensitive vibro-transducers permits unique evaluation of real stressed-strained states of power aggregates and loaded constructions and opens fundamental innovation opportunities: conduct of continuous (on-line regime) reliable monitoring of turboagregates of electrical machines, compressor installations, bases, supports, pipe-lines and other objects subjected to damaging effect of vibrations; control of operational safety of technical systems at all the stages of life cycle including design, test production, tuning, testing, operational use, repairs and resource enlargement; creation of vibro-diagnostic systems of authentic non-destructive control of anisotropic characteristics of materials resistance of power aggregates and loaded constructions under outer effects and operational flaws. The described technology is revolutionary, universal and common for all branches of engineering industry and construction building objects.
NASA Astrophysics Data System (ADS)
Polkowski, Marcin; Grad, Marek
2015-04-01
The 3D P-wave seismic velocity model was obtained by combining data from multiple studies during past 50 years. Data sources included refraction seismology, reflection seismology, geological boreholes, vertical seismic profiling, magnetotellurics and gravimetry. Use of many data sources allowed creation of detailed 3D P-wave velocity model that reaches to depth of 60 km and includes 6-layers of sediments and 3-layers of the crust. Purpose of this study is to analyze how 3D model influences local (accuracy of location and source time estimation for local events), regional (identification of wide-angle seismic phases) and global (teleseismic tomography) seismic travel times. Additionally we compare results of forward seismic wave propagation with signals observed on short period and broadband stations. National Science Centre Poland provided financial support for this work by NCN grant DEC-2011/02/A/ST10/00284.
3D Anisotropic structure of the south-central Mongolia from Rayleigh and Love wave tomography
NASA Astrophysics Data System (ADS)
Yu, D.; Wu, Q.; Montagner, J. P.
2014-12-01
A better understanding of the geodynamics of the crust and mantle below Baikal-Mongolia is required to identify the role of mantle processes versus that of far-field tectonic effects from India-Asia collision. Anisotropy tomography can provide new perspective to the continental growth mechanism. In order to study the 3D anisotropic structure of the upper mantle in the south-central Mongolia, we collected the vertical and transverse components of seismograms recorded at 69 broadband seismic stations. We have measured inter-station phase velocities of 7181 Rayleigh waves and 901 Love waves using the frequency-time analysis of wavelet transformation method for the fundamental mode at period range 10~80s. The lateral phase velocity variations are computed by using a regionalization method. These phase velocities have been inverted to obtain the first anisotropic model including Sv velocities, azimuthal and radial anisotropy. The Middle Gobi is associated with low velocity. Based on the distribution of the Cenozoic basalts in the Middle Gobi, it refers that the low velocity anomaly is related to the Cenozoic volcanism. In the northern domain, near to Baikal zone, the azimuthal anisotropy is normal to the Baikal rift and consistent with the fast direction of previous SKS splitting measurements. In the South Gobi, north to Main Mongolian Lineament, the azimuthal anisotropy is NEE-SWW in the crust and NW-SE in the mantle. It indicates that the crust and mantle are decoupled. We propose that the crustal deformation is related to the far-field effects of India-Asia collision and that the mantle flow is correlated with the Baikal rift activity. Further study in process will provide more evidence and insight to better understand the geodynamics in this region.
Seismic waves in 3-D: from mantle asymmetries to reliable seismic hazard assessment
NASA Astrophysics Data System (ADS)
Panza, Giuliano F.; Romanelli, Fabio
2014-10-01
A global cross-section of the Earth parallel to the tectonic equator (TE) path, the great circle representing the equator of net lithosphere rotation, shows a difference in shear wave velocities between the western and eastern flanks of the three major oceanic rift basins. The low-velocity layer in the upper asthenosphere, at a depth range of 120 to 200 km, is assumed to represent the decoupling between the lithosphere and the underlying mantle. Along the TE-perturbed (TE-pert) path, a ubiquitous LVZ, about 1,000-km-wide and 100-km-thick, occurs in the asthenosphere. The existence of the TE-pert is a necessary prerequisite for the existence of a continuous global flow within the Earth. Ground-shaking scenarios were constructed using a scenario-based method for seismic hazard analysis (NDSHA), using realistic and duly validated synthetic time series, and generating a data bank of several thousands of seismograms that account for source, propagation, and site effects. Accordingly, with basic self-organized criticality concepts, NDSHA permits the integration of available information provided by the most updated seismological, geological, geophysical, and geotechnical databases for the site of interest, as well as advanced physical modeling techniques, to provide a reliable and robust background for the development of a design basis for cultural heritage and civil infrastructures. Estimates of seismic hazard obtained using the NDSHA and standard probabilistic approaches are compared for the Italian territory, and a case-study is discussed. In order to enable a reliable estimation of the ground motion response to an earthquake, three-dimensional velocity models have to be considered, resulting in a new, very efficient, analytical procedure for computing the broadband seismic wave-field in a 3-D anelastic Earth model.
Propagation of Electromagnetic Waves in 3D Opal-based Magnetophotonic Crystals
NASA Astrophysics Data System (ADS)
Pardavi-Horvath, Martha; Makeeva, Galina S.; Golovanov, Oleg A.; Rinkevich, Anatolii B.
2013-03-01
Opals, a class of self-organized 3D nanostructures, are typical representatives of photonic bandgap structures. The voids inside of the opal structure of close packed SiO2 spheres can be infiltrated by a magnetic material, creating magnetically tunable magnetophotonic crystals with interesting and potentially useful properties at GHz and THz frequencies. The propagation of electromagnetic waves at microwave frequencies was investigated numerically in SiO2 opal based magnetic nanostructures, using rigorous mathematical models to solve Maxwell's equations complemented by the Landau-Lifshitz equation with electrodynamic boundary conditions. The numerical approach is based on Galerkin's projection method using the decomposition algorithm on autonomous blocks with Floquet channels. The opal structure consists of SiO2 nanospheres, with inter-sphere voids infiltrated with nanoparticles of Ni-Zn ferrites. Both the opal matrix and the ferrite are assumed to be lossy. A model, taking into account the real structure of the ferrite particles in the opal's voids was developed to simulate the measured FMR lineshape of the ferrite infiltrated opal. The numerical technique shows an excellent agreement when applied to model recent experimental data on similar ferrite opals.
Nonhydrostatic granular flow over 3-D terrain: New Boussinesq-type gravity waves?
NASA Astrophysics Data System (ADS)
Castro-Orgaz, Oscar; Hutter, Kolumban; Giraldez, Juan V.; Hager, Willi H.
2015-01-01
granular mass flow is a basic step in the prediction and control of natural or man-made disasters related to avalanches on the Earth. Savage and Hutter (1989) pioneered the mathematical modeling of these geophysical flows introducing Saint-Venant-type mass and momentum depth-averaged hydrostatic equations using the continuum mechanics approach. However, Denlinger and Iverson (2004) found that vertical accelerations in granular mass flows are of the same order as the gravity acceleration, requiring the consideration of nonhydrostatic modeling of granular mass flows. Although free surface water flow simulations based on nonhydrostatic depth-averaged models are commonly used since the works of Boussinesq (1872, 1877), they have not yet been applied to the modeling of debris flow. Can granular mass flow be described by Boussinesq-type gravity waves? This is a fundamental question to which an answer is required, given the potential to expand the successful Boussinesq-type water theory to granular flow over 3-D terrain. This issue is explored in this work by generalizing the basic Boussinesq-type theory used in civil and coastal engineering for more than a century to an arbitrary granular mass flow using the continuum mechanics approach. Using simple test cases, it is demonstrated that the above question can be answered in the affirmative way, thereby opening a new framework for the physical and mathematical modeling of granular mass flow in geophysics, whereby the effect of vertical motion is mathematically included without the need of ad hoc assumptions.
Characterization of an SRF gun: a 3D full wave simulation
Wang, E.; Ben-Zvi, I.; Wang, J.
2011-03-28
We characterized a BNL 1.3GHz half-cell SRF gun is tested for GaAs photocathode. The gun already was simulated several years ago via two-dimensional (2D) numerical codes (i.e., Superfish and Parmela) with and without the beam. In this paper, we discuss our investigation of its characteristics using a three dimensional (3D) full-wave code (CST STUDIO SUITE{trademark}).The input/pickup couplers are sited symmetrically on the same side of the gun at an angle of 180{sup o}. In particular, the inner conductor of the pickup coupler is considerably shorter than that of the input coupler. We evaluated the cross-talk between the beam (trajectory) and the signal on the input coupler compared our findings with published results based on analytical models. The CST STUDIO SUITE{trademark} also was used to predict the field within the cavity; particularly, a combination of transient/eigenmode solvers was employed to accurately construct the RF field for the particles, which also includes the effects of the couplers. Finally, we explored the beam's dynamics with a particle in cell (PIC) simulation, validated the results and compare them with 2D code result.
Interaction of highly nonlinear solitary waves with linear elastic media
NASA Astrophysics Data System (ADS)
Yang, Jinkyu; Silvestro, Claudio; Khatri, Devvrath; de Nardo, Luigi; Daraio, Chiara
2011-04-01
We study the interaction of highly nonlinear solitary waves propagating in granular crystals with an adjacent linear elastic medium. We investigate the effects of interface dynamics on the reflection of incident waves and on the formation of primary and secondary reflected waves. Experimental tests are performed to correlate the linear medium geometry, materials, and mass with the formation and propagation of reflected waves. We compare the experimental results with theoretical analysis based on the long-wavelength approximation and with numerical predictions obtained from discrete particle models. Experimental results are found to be in agreement with theoretical analysis and numerical simulations. This preliminary study establishes the foundation for utilizing reflected solitary waves as novel information carriers in nondestructive evaluation of elastic material systems.
Interaction of highly nonlinear solitary waves with linear elastic media.
Yang, Jinkyu; Silvestro, Claudio; Khatri, Devvrath; De Nardo, Luigi; Daraio, Chiara
2011-04-01
We study the interaction of highly nonlinear solitary waves propagating in granular crystals with an adjacent linear elastic medium. We investigate the effects of interface dynamics on the reflection of incident waves and on the formation of primary and secondary reflected waves. Experimental tests are performed to correlate the linear medium geometry, materials, and mass with the formation and propagation of reflected waves. We compare the experimental results with theoretical analysis based on the long-wavelength approximation and with numerical predictions obtained from discrete particle models. Experimental results are found to be in agreement with theoretical analysis and numerical simulations. This preliminary study establishes the foundation for utilizing reflected solitary waves as novel information carriers in nondestructive evaluation of elastic material systems. PMID:21599325
Theoretical relationship between elastic wave velocity and electrical resistivity
NASA Astrophysics Data System (ADS)
Lee, Jong-Sub; Yoon, Hyung-Koo
2015-05-01
Elastic wave velocity and electrical resistivity have been commonly applied to estimate stratum structures and obtain subsurface soil design parameters. Both elastic wave velocity and electrical resistivity are related to the void ratio; the objective of this study is therefore to suggest a theoretical relationship between the two physical parameters. Gassmann theory and Archie's equation are applied to propose a new theoretical equation, which relates the compressional wave velocity to shear wave velocity and electrical resistivity. The piezo disk element (PDE) and bender element (BE) are used to measure the compressional and shear wave velocities, respectively. In addition, the electrical resistivity is obtained by using the electrical resistivity probe (ERP). The elastic wave velocity and electrical resistivity are recorded in several types of soils including sand, silty sand, silty clay, silt, and clay-sand mixture. The appropriate input parameters are determined based on the error norm in order to increase the reliability of the proposed relationship. The predicted compressional wave velocities from the shear wave velocity and electrical resistivity are similar to the measured compressional velocities. This study demonstrates that the new theoretical relationship may be effectively used to predict the unknown geophysical property from the measured values.
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.
Coal Thickness Gauging Using Elastic Waves
NASA Technical Reports Server (NTRS)
Nazarian, Soheil; Bar-Cohen, Yoseph
1999-01-01
The efforts of a mining crew can be optimized, if the thickness of the coal layers to be excavated is known before excavation. Wave propagation techniques can be used to estimate the thickness of the layer based on the contrast in the wave velocity between coal and rock beyond it. Another advantage of repeated wave measurement is that the state of the stress within the mine can be estimated. The state of the stress can be used in many safety-related decisions made during the operation of the mine. Given these two advantages, a study was carried out to determine the feasibility of the methodology. The results are presented herein.
Scattering and coupling effects of electromagnetic waves in 3D networks of spheres
NASA Astrophysics Data System (ADS)
Defos Du Rau, M.; Pessan, F.; Ruffie, G.; Vignéras-Lefebvre, V.; Parneix, J. P.
1998-01-01
In this paper, the problem of electromagnetic scattering from a 3D system of spheres is considered and an iterative solution that accounts for multiple scattering is proposed. The Mie formalism used for a single sphere is extended to account for multiple scattered fields between several particles. The translational addition theorems for spherical wave functions are used to express the electromagnetic field scattered by a sphere S_i in terms of an incident field for a sphere S_k in a spherical coordinates system attached to the sphere S_k. In this work, the numerical convergence of the method is discussed and associated computational times are given. Numerical computations including Radar Cross Section (RCS) and radiation patterns for various 3D configurations are presented. Some of them are compared with free-space measurements made in the 8 to 100 GHz frequency band using vectorial network analyzers. 11.55.-m S-matrix theory; analytic structure of amplitudes Cet article étudie la diffusion des ondes électromagnétiques par des réseaux tridimensionnels de sphères et propose une méthode itérative pour prendre en compte les effets de multidiffusion. Le formalisme de Mie utilisé dans le cas d'une sphère est étendu pour calculer les champs "multidiffusés" entre plusieurs particules. Les théorèmes d'addition et de translation des fonctions d'onde sphériques sont utilisés pour exprimer le champ diffusé par une sphère S_i comme étant incident sur une sphère S_k, dans un système de coordonnées sphériques lié au centre de S_k. La convergence numérique de la méthode est discutée et des temps de calcul sont donnés. Des résultats numériques tels que des Surfaces Équivalentes Radar (SER) et des diagrammes de rayonnement pour différentes configurations tridimensionnelles sont montrés. Certains d'entre eux sont comparés à des mesures en espace libre faites à l'aide d'analyseurs de réseaux vectoriels dans la bande de fréquence 8{-}100 GHz.
Modulated pressure waves in large elastic tubes.
Mefire Yone, G R; Tabi, C B; Mohamadou, A; Ekobena Fouda, H P; Kofané, T C
2013-09-01
Modulational instability is the direct way for the emergence of wave patterns and localized structures in nonlinear systems. We show in this work that it can be explored in the framework of blood flow models. The whole modified Navier-Stokes equations are reduced to a difference-differential amplitude equation. The modulational instability criterion is therefore derived from the latter, and unstable patterns occurrence is discussed on the basis of the nonlinear parameter model of the vessel. It is found that the critical amplitude is an increasing function of α, whereas the region of instability expands. The subsequent modulated pressure waves are obtained through numerical simulations, in agreement with our analytical expectations. Different classes of modulated pressure waves are obtained, and their close relationship with Mayer waves is discussed. PMID:24089964
Modulated pressure waves in large elastic tubes
NASA Astrophysics Data System (ADS)
Mefire Yone, G. R.; Tabi, C. B.; Mohamadou, A.; Ekobena Fouda, H. P.; Kofané, T. C.
2013-09-01
Modulational instability is the direct way for the emergence of wave patterns and localized structures in nonlinear systems. We show in this work that it can be explored in the framework of blood flow models. The whole modified Navier-Stokes equations are reduced to a difference-differential amplitude equation. The modulational instability criterion is therefore derived from the latter, and unstable patterns occurrence is discussed on the basis of the nonlinear parameter model of the vessel. It is found that the critical amplitude is an increasing function of α, whereas the region of instability expands. The subsequent modulated pressure waves are obtained through numerical simulations, in agreement with our analytical expectations. Different classes of modulated pressure waves are obtained, and their close relationship with Mayer waves is discussed.
Waves guided by a thin viscoelastic layer between elastic solids
NASA Technical Reports Server (NTRS)
Claus, R. O.; Roges, R. T.
1981-01-01
The propagation of ultrasonic waves guided along a viscoelastic layer which separates two dissimilar elastic solid half spaces is described. For the limiting cases of rigid and soft bonding by a layer thin compared to acoustic wavelength, Stoneley and independent Rayleigh wave characteristic equations, respectively, result. For combinations of layer rigidity and wave length that correspond to guided wave propagation, external cyclic loading of the layer produces mechanical hysteresis and a resulting dynamic change in bond properties. Measurements of velocity hysteresis in an aluminum polymer adhesive aluminum system are described.
Elastic waves in discontinuous media: Three-dimensional scattering
NASA Astrophysics Data System (ADS)
Molino, F. R.; Sabatier, P. C.
1994-09-01
This report contains an exact study of elastic wave propagation and its scattering in discontinuous media where hard reflectors are onionlike sets of surfaces. In order to reformulate the problem as a finite set of boundary integral equations, the wave motion between reflectors is represented by means of elastic potentials which involve vectorial densities on the surfaces. In the external medium, an outgoing asymptotic condition generalizes the Silver-Müller (and the Sommerfeld) condition to the case of coupled waves (S and P waves) moving with different velocities. The uniqueness of the Green's function, which guarantees the uniqueness of the direct problem solution, is proven. For any incident wave and arbitrary number of surfaces, the transmission and scattering problems are studied, with and without the simplification obtained by assuming constant Poisson ratios. According to the parameter ranges, the equations which are obtained are well posed, either as second kind Fredholm equations, or because they reduce to the inverse of the sum of the identity operator and a ``small norm'' bounded operator. The results can be used to describe rigorously the three-dimensional scattering of elastic waves in the frequency domain for any kind of incident wave function (P,S,...) as well as the response to a localized source.
High Frequency Elastic Wave Propagation in Media with a Microstructure
NASA Astrophysics Data System (ADS)
Tie, B.; Aubry, D.; Mouronval, A.-S.; Solas, D.; Thébault, J.; Tian, B.-Y.
2010-05-01
This contribution deals with the theoretical analysis and numerical modeling of elastic wave propagation in media with a microstructure. Two kinds of media are considered: polycrystalline material and honeycomb core sandwich shells, in which elastic waves are triggered by transient signals that result in large frequency ranges including high frequencies. Our theoretical and numerical investigations aim at understanding and simulating the interactions between the microstructure of those media and the wave propagation phenomena, when the characteristic lengths of the microstructure and the involved shortest wavelengths have roughly the same scale. In this paper, some key mechanisms of interaction between the considered microstructures and the elastic waves are highlighted. In polycrystalline superalloys, the misorientation distribution and the average grain size are considered, as they can alter pressure/shear wave propagation and also the permeability to ultrasonic waves monitored to perform non-destructive testing. For the flexure behavior of honeycomb core sandwich shells, the fundamental role played by the honeycomb cells, especially in high frequency domain, is analyzed. Relevant numerical modeling that provides a promising way to quantify micro-structure/wave interactions is presented. The important issue of how to take into account these micro-scale interactions in a homogenized macro-scale modeling is also discussed.
Measurements of radiated elastic wave energy from dynamic tensile cracks
NASA Technical Reports Server (NTRS)
Boler, Frances M.
1990-01-01
The role of fracture-velocity, microstructure, and fracture-energy barriers in elastic wave radiation during a dynamic fracture was investigated in experiments in which dynamic tensile cracks of two fracture cofigurations of double cantilever beam geometry were propagating in glass samples. The first, referred to as primary fracture, consisted of fractures of intact glass specimens; the second configuration, referred to as secondary fracture, consisted of a refracture of primary fracture specimens which were rebonded with an intermittent pattern of adhesive to produce variations in fracture surface energy along the crack path. For primary fracture cases, measurable elastic waves were generated in 31 percent of the 16 fracture events observed; the condition for radiation of measurable waves appears to be a local abrupt change in the fracture path direction, such as occurs when the fracture intersects a surface flaw. For secondary fractures, 100 percent of events showed measurable elastic waves; in these fractures, the ratio of radiated elastic wave energy in the measured component to fracture surface energy was 10 times greater than for primary fracture.
SHEAR WAVE SEISMIC STUDY COMPARING 9C3D SV AND SH IMAGES WITH 3C3D C-WAVE IMAGES
John Beecherl; Bob A. Hardage
2004-07-01
The objective of this study was to compare the relative merits of shear-wave (S-wave) seismic data acquired with nine-component (9-C) technology and with three-component (3-C) technology. The original proposal was written as if the investigation would be restricted to a single 9-C seismic survey in southwest Kansas (the Ashland survey), on the basis of the assumption that both 9-C and 3-C S-wave images could be created from that one data set. The Ashland survey was designed as a 9-C seismic program. We found that although the acquisition geometry was adequate for 9-C data analysis, the source-receiver geometry did not allow 3-C data to be extracted on an equitable and competitive basis with 9-C data. To do a fair assessment of the relative value of 9-C and 3-C seismic S-wave data, we expanded the study beyond the Ashland survey and included multicomponent seismic data from surveys done in a variety of basins. These additional data were made available through the Bureau of Economic Geology, our research subcontractor. Bureau scientists have added theoretical analyses to this report that provide valuable insights into several key distinctions between 9-C and 3-C seismic data. These theoretical considerations about distinctions between 3-C and 9-C S-wave data are presented first, followed by a discussion of differences between processing 9-C common-midpoint data and 3-C common-conversion-point data. Examples of 9-C and 3-C data are illustrated and discussed in the last part of the report. The key findings of this study are that each S-wave mode (SH-SH, SV-SV, or PSV) involves a different subsurface illumination pattern and a different reflectivity behavior and that each mode senses a different Earth fabric along its propagation path because of the unique orientation of its particle-displacement vector. As a result of the distinct orientation of each mode's particle-displacement vector, one mode may react to a critical geologic condition in a more optimal way than do
NASA Astrophysics Data System (ADS)
Tong, Ping; Komatitsch, Dimitri; Tseng, Tai-Lin; Hung, Shu-Huei; Chen, Chin-Wu; Basini, Piero; Liu, Qinya
2014-10-01
We present a three-dimensional (3-D) hybrid method that interfaces the spectral-element method (SEM) with the frequency-wave number (FK) technique to model the propagation of teleseismic plane waves beneath seismic arrays. The accuracy of the resulting 3-D SEM-FK hybrid method is benchmarked against semianalytical FK solutions for 1-D models. The accuracy of 2.5-D modeling based on 2-D SEM-FK hybrid method is also investigated through comparisons to this 3-D hybrid method. Synthetic examples for structural models of the Alaska subduction zone and the central Tibet crust show that this method is capable of accurately capturing interactions between incident plane waves and local heterogeneities. This hybrid method presents an essential tool for the receiver function and scattering imaging community to verify and further improve their techniques. These numerical examples also show the promising future of the 3-D SEM-FK hybrid method in high-resolution regional seismic imaging based on waveform inversions of converted/scattered waves recorded by seismic array.
Finite-difference staggered grids in GPUs for anisotropic elastic wave propagation simulation
NASA Astrophysics Data System (ADS)
Rubio, Felix; Hanzich, Mauricio; Farrés, Albert; de la Puente, Josep; María Cela, José
2014-09-01
The 3D elastic wave equations can be used to simulate the physics of waves traveling through the Earth more precisely than acoustic approximations. However, this improvement in quality has a counterpart in the cost of the numerical scheme. A possible strategy to mitigate that expense is using specialized, high-performing architectures such as GPUs. Nevertheless, porting and optimizing a code for such a platform require a deep understanding of both the underlying hardware architecture and the algorithm at hand. Furthermore, for very large problems, multiple GPUs must work concurrently, which adds yet another layer of complexity to the codes. In this work, we have tackled the problem of porting and optimizing a 3D elastic wave propagation engine which supports both standard- and fully-staggered grids to multi-GPU clusters. At the single GPU level, we have proposed and evaluated many optimization strategies and adopted the best performing ones for our final code. At the distributed memory level, a domain decomposition approach has been used which allows for good scalability thanks to using asynchronous communications and I/O.
Transformation ray method: controlling high frequency elastic waves (L).
Chang, Zheng; Liu, Xiaoning; Hu, Gengkai; Hu, Jin
2012-10-01
Elastic ray theory is a high frequency asymptotic approximation of solution of elastodynamic equation, and is widely used in seismology. In this paper, the form invariance under a general spatial mapping and high frequency wave control have been examined by transformation method. It is showed that with the constraint of major and minor symmetry of the transformed elastic tensor, the eikonal equation keeps its form under a general mapping, however, the transport equation loses its form except for conformal mapping. Therefore, the elastic ray path can be controlled in an exact manner by a transformation method, whereas energy distribution along the ray is only approximately controlled. An elastic rotator based on ray tracing method is also provided to illustrate the method and to access the approximation. PMID:23039561
Elastic responses of a flotation ring in water waves
NASA Astrophysics Data System (ADS)
Dong, Guo-Hai; Hao, Shuang-Hu; Zhao, Yun-Peng; Zong, Zhi; Gui, Fu-Kun
2010-01-01
The gravity-type fish cage is extensively applied in open-sea fishery aquaculture. Its practicality is closely related to the reliability of the flotation ring which is its main load-bearing component. Therefore, it is necessary to study the elastic responses of the flotation ring in ocean waves. Here, an analytical method is proposed to analyze the elastic deformations of a circular ring subjected to water waves. The governing equations of six degree-of-freedom motions and elastic deformations are obtained according to Euler's laws and curved beam theory. In order to examine the method, a series of physical model tests were carried out. The surge and heave displacements of the ring between the predicted results and experimental measurements are compared, and good correlation is represented. The effects of the propagation directions of the incident wave on elastic responses of the ring are then analyzed. It is concluded that small deformations of the ring occur when the configuration of the mooring cables is symmetrically arranged along the propagation direction of the incident waves. Additionally, the out-of-plane stiffness is suggested to be strengthened in order to diminish the corresponding deformations.
3D Hot Test Simulations of a 220 GHz Folded Waveguide Traveling Wave Tube Using a CFDTD PIC Method
NASA Astrophysics Data System (ADS)
Lin, Ming-Chieh; Song, Heather
2015-11-01
Millimeter or sub-THz wave sources centered at 220 GHz is of interest due to the potential for its commercial and military applications including high resolution radar, remote sensing, and high-data-rate communications. It has been demonstrated via 3D cold test finite element method (FEM) simulations that a folded waveguide traveling wave tube (FWTWT) can be designed and optimized at this frequency range with a small signal gain of 18 dB over a comparatively broad (-3 dB) bandwidth of ~ 10%. On the other hand, 3D hot test simulations of a V-band ladder TWT have been successfully demonstrated using a conformal finite-difference time-domain (CFDTD) particle-in-cell (PIC) method for center frequency of 50 GHz. In the present work, the 220 GHz FWTWT designs have been reviewed and studied. 3D Cold test simulations using both the CFDTD and FEM methods have been carried out and compared with each other as basis for 3D hot test CFDTD PIC simulations. The preliminary simulation result shows that the gain-bandwidth features at 220 GHz are achievable while carefully avoiding beam interceptions. Our study shows that the interaction characteristics are very sensitive to the operating beam parameters. Detail simulation results and discussions will be presented.
3D millimeter wave imaging of vertical cracks and its application for the inspection of HDPE pipes
NASA Astrophysics Data System (ADS)
Ghasr, Mohammad Tayeb; Ying, Kuang; Zoughi, Reza
2014-02-01
Robust detection of vertical cracks in high-density polyethylene (HDPE) pipes is a challenging task for the majority of nondestructive testing (NDT) techniques. Vertical cracks are specifically referred to those whose largest planar view is parallel to the signal direction of propagation, leaving very little signal to be scattered for detection. In such pipes this commonly occurs between two pipes sections when thermally or adhesively joined. This work presents the utility and efficacy of three-dimensional (3D) millimeter wave holographical imaging based on synthetic aperture radar (SAR) algorithm for imaging such cracks. Such a 3D millimeter wave image can readily represent the type, size, and location of various flaws within a structure. Two-dimensional (2D) slices of the 3D image, at different orientations, can also be readily produced showing the cross-sectional views of the structure and flaws, further aiding in identifying, and sizing a flaw or vertical crack. Imaging results for planner and curved (pipe section) specimen with machined flaws are presented. These images are produced using a novel field-portable, small, and low-cost wideband millimeter-wave reflectometer capable of rapid 3D image production.
CIP-MOC Modeling of Seismic Wave Propagation in Elastic Media
NASA Astrophysics Data System (ADS)
Yoshimi, M.
2004-12-01
In many fields such as hydrodynamics and MHD, the CIP method, an upwind difference hyperbolic equation solver, has widely been employed for advection calculation. The CIP scheme was constructed considering that an advected property and its spatial derivative follow same advection equation. This effects low numerical dispersion and relaxed CFL condition in the advection calculation. In the present work, we developed a CIP-MOC (CIP with method of characteristics) scheme for seismic wave propagation in 3D elastic heterogeneous media with flat free surface. 3D elastic wave equations in velocity-stress formulation and their spatial derivatives, as well, are converted into sets of 1D advection equations and non-advection equations for each direction (x,y,z in Cartesian coodinate system) with the method of characteristics. Since the Riemann invariant of each advection equation consists of stress and velocity, updatings of velocity and stress are simultaneous and a collocated grid system is employed. A free surface is modeled as a zero-stress surface. A reflection free boundary is installed by considering no incident wave comes from outside of the boundary. A double coupled seismic point source is introduced as external point stresses. Overall scheme is made up of multiphases employing time-splitting and directional-splitting techniques. Each time step is composed of three directional updating phases each for wave propagation in x, y and z direction. Each directional updating phase is made up of advection phase and non-advection phase. In the advection phase, advection equations are solved with the CIP method. In the non-advection phases, non-advetion equations and boundary conditions are evaluated with central finite differences. We conducted CIP-MOC seismic wave propagation simulations in a half-space, layered and fully heterogeneous media for embedded point source. By comparing our products with those produced with discrete wavenumber method and finite difference method
NASA Astrophysics Data System (ADS)
Velioǧlu, Deniz; Cevdet Yalçıner, Ahmet; Zaytsev, Andrey
2016-04-01
Tsunamis are huge waves with long wave periods and wave lengths that can cause great devastation and loss of life when they strike a coast. The interest in experimental and numerical modeling of tsunami propagation and inundation increased considerably after the 2011 Great East Japan earthquake. In this study, two numerical codes, FLOW 3D and NAMI DANCE, that analyze tsunami propagation and inundation patterns are considered. Flow 3D simulates linear and nonlinear propagating surface waves as well as long waves by solving three-dimensional Navier-Stokes (3D-NS) equations. NAMI DANCE uses finite difference computational method to solve 2D depth-averaged linear and nonlinear forms of shallow water equations (NSWE) in long wave problems, specifically tsunamis. In order to validate these two codes and analyze the differences between 3D-NS and 2D depth-averaged NSWE equations, two benchmark problems are applied. One benchmark problem investigates the runup of long waves over a complex 3D beach. The experimental setup is a 1:400 scale model of Monai Valley located on the west coast of Okushiri Island, Japan. Other benchmark problem is discussed in 2015 National Tsunami Hazard Mitigation Program (NTHMP) Annual meeting in Portland, USA. It is a field dataset, recording the Japan 2011 tsunami in Hilo Harbor, Hawaii. The computed water surface elevation and velocity data are compared with the measured data. The comparisons showed that both codes are in fairly good agreement with each other and benchmark data. The differences between 3D-NS and 2D depth-averaged NSWE equations are highlighted. All results are presented with discussions and comparisons. Acknowledgements: Partial support by Japan-Turkey Joint Research Project by JICA on earthquakes and tsunamis in Marmara Region (JICA SATREPS - MarDiM Project), 603839 ASTARTE Project of EU, UDAP-C-12-14 project of AFAD Turkey, 108Y227, 113M556 and 213M534 projects of TUBITAK Turkey, RAPSODI (CONCERT_Dis-021) of CONCERT
NASA Technical Reports Server (NTRS)
Temmer, M.; Veronig, A. M.; Gopalswamy, N.; Yashiro, S.
2011-01-01
We study the kinematical characteristics and 3D geometry of a large-scale coronal wave that occurred in association with the 26 April 2008 flare-CME event. The wave was observed with the EUVI instruments aboard both STEREO spacecraft (STEREO-A and STEREO-B) with a mean speed of approx 240 km/s. The wave is more pronounced in the eastern propagation direction, and is thus, better observable in STEREO-B images. From STEREO-B observations we derive two separate initiation centers for the wave, and their locations fit with the coronal dimming regions. Assuming a simple geometry of the wave we reconstruct its 3D nature from combined STEREO-A and STEREO-B observations. We find that the wave structure is asymmetric with an inclination toward East. The associated CME has a deprojected speed of approx 750 +/- 50 km/s, and it shows a non-radial outward motion toward the East with respect to the underlying source region location. Applying the forward fitting model developed by Thernisien, Howard, and Vourlidas we derive the CME flux rope position on the solar surface to be close to the dimming regions. We conclude that the expanding flanks of the CME most likely drive and shape the coronal wave.
Propagation of elastic waves through textured polycrystals: application to ice
Maurel, Agnès; Lund, Fernando; Montagnat, Maurine
2015-01-01
The propagation of elastic waves in polycrystals is revisited, with an emphasis on configurations relevant to the study of ice. Randomly oriented hexagonal single crystals are considered with specific, non-uniform, probability distributions for their major axis. Three typical textures or fabrics (i.e. preferred grain orientations) are studied in detail: one cluster fabric and two girdle fabrics, as found in ice recovered from deep ice cores. After computing the averaged elasticity tensor for the considered textures, wave propagation is studied using a wave equation with elastic constants c=〈c〉+δc that are equal to an average plus deviations, presumed small, from that average. This allows for the use of the Voigt average in the wave equation, and velocities are obtained solving the appropriate Christoffel equation. The velocity for vertical propagation, as appropriate to interpret sonic logging measurements, is analysed in more details. Our formulae are shown to be accurate at the 0.5% level and they provide a rationale for previous empirical fits to wave propagation velocities with a quantitative agreement at the 0.07–0.7% level. We conclude that, within the formalism presented here, it is appropriate to use, with confidence, velocity measurements to characterize ice fabrics. PMID:27547099
Topological Phononic Crystals with One-Way Elastic Edge Waves
NASA Astrophysics Data System (ADS)
Wang, Pai; Lu, Ling; Bertoldi, Katia
2015-09-01
We report a new type of phononic crystals with topologically nontrivial band gaps for both longitudinal and transverse polarizations, resulting in protected one-way elastic edge waves. In our design, gyroscopic inertial effects are used to break the time-reversal symmetry and realize the phononic analogue of the electronic quantum (anomalous) Hall effect. We investigate the response of both hexagonal and square gyroscopic lattices and observe bulk Chern numbers of 1 and 2, indicating that these structures support single and multimode edge elastic waves immune to backscattering. These robust one-way phononic waveguides could potentially lead to the design of a novel class of surface wave devices that are widely used in electronics, telecommunication, and acoustic imaging.
Elastic wave velocities of Apollo 14, 15, and 16 rocks
NASA Technical Reports Server (NTRS)
Mizutani, H.; Newbigging, D. F.
1973-01-01
Elastic wave velocities of two Apollo 14 rocks, 14053 and 14321, three Apollo 15 rocks, 15058, 15415, and 15545, and one Apollo 16 rock 60315 have been determined at pressures up to 10 kb. For sample 14321, the variation of the compressional wave velocities with temperature has been measured over the temperature range from 27 to 200 C. Overall elastic properties of these samples except sample 15415 are very similar to those of Apollo 11, 12, and 14 rocks and are concordant with Toksoz et al.'s (1972) interpretation that lunar upper crust is of basaltic composition. Temperature derivative of the P wave velocity for sample 14321 is a half to one order of magnitude larger than that for single crystalline minerals. This suggests that the seismic velocity in the lunar crust may be affected significantly by the temperature distribution.
Joint evaluation of fracture azimuth by electromagnetic wave and elastic wave
NASA Astrophysics Data System (ADS)
Feng, Xuan; Liu, Cai; Wang, Qiao; Wang, Kai; Lu, Qi; Xue, Jian; Liang, Wenjing; Yu, Yue; Ren, Qianci
2013-12-01
With the multi-wave, multi-component seismic wave exploration, one can apply the anisotropy of fracture media to analyze the attributes of the fracture media, including the fracture azimuth. In the meantime, the techniques of full-polarimetric electromagnetic wave, including full-polarimetric borehole radar, can also be used to analyze the attributes of the fracture. However, the analysis precision of both the multi-component elastic wave exploration and full-polarimetric electromagnetic wave exploration is prone to the influence of noise and other factors. So far, some researchers have conducted studies on the joint inversion of electromagnetic waves and seismic waves. This paper develops evaluation techniques of fracture azimuth by electromagnetic wave, elastic wave, and joint analysis of coincident elastic reflection and electromagnetic data. Firstly, based on the shear wave splitting of elastic waves, this paper develops a statistical analysis technique which applies Pearson correlation coefficient to count and analyze the azimuth angle of fracture. Secondly, based on the information of electromagnetic polarization rotated by fracture, this paper develops a statistical analysis method of full-polarimetric electromagnetic waves which applies the maximum amplitude ratio between the co-polarization and cross-polarization to analyze the azimuth angle of fracture. Furthermore, based on the analysis result of the elastic wave and full-polarimetric electromagnetic wave, this paper develops a joint analysis technique which adopts the standard deviation. At last, authors in this study conduct joint detection experiments on the coincident fracture medium by using the ultrasonic and full-polarimetric ground penetrating radar. The experimental result indicates that both single geophysical methods are capable of analyzing the fracture azimuth angle, but the joint analysis is more accurate.
NASA Astrophysics Data System (ADS)
Gumerov, Nail A.; Karavaev, Alexey V.; Surjalal Sharma, A.; Shao, Xi; Papadopoulos, Konstantinos D.
2011-04-01
Efficient spectral and pseudospectral algorithms for simulation of linear and nonlinear 3D whistler waves in a cold electron plasma are developed. These algorithms are applied to the simulation of whistler waves generated by loop antennas and spheromak-like stationary waves of considerable amplitude. The algorithms are linearly stable and show good stability properties for computations of nonlinear waves over tens of thousands of time steps. Additional speedups by factors of 10-20 (comparing single core CPU and one GPU) are achieved by using graphics processors (GPUs), which enable efficient numerical simulation of the wave propagation on relatively high resolution meshes (tens of millions nodes) in personal computing environment. Comparisons of the numerical results with analytical solutions and experiments show good agreement. The limitations of the codes and the performance of the GPU computing are discussed.
Internal waves patterns in the wake of a 3D body towed in a two-layer fluid
NASA Astrophysics Data System (ADS)
Lacaze, Laurent; Mercier, Matthieu; Thual, Olivier; Paci, Alexandre
2014-11-01
Stratified flows over obstacles are important features in meteorology and oceanography. The characterization of these flows is crucial in order to propose models of geophysical processes such as mixing and ocean circulation or orographic drag in the atmosphere. For some specific stratification profiles, the energy of internal waves generated by the obstacle can be trapped at a given depth, at the base of the oceanic mixing layer or at the top of the atmospheric boundary layer for instance. This scenario can be modelled by a two-layer stratified fluid for which gravity waves spread at the interface between the two layers. The work presented here focuses on a two-layer flow over a 3D obstacle, or equivalently, an obstacle towed in a fluid at rest. Experiments performed both in the large-scale flume of CNRM-GAME Toulouse (METEO-FRANCE & CNRS) and in a smaller tank apparatus, are presented with a specific attention on the measurement of the 3D wave patterns. A non-hydrostatic linear analysis is used to describe the observed wave patterns. The experiments highlight the strong influence of the Froude number on the generated waves. More specifically, we investigate the nature of the wake angle obtained from the wave pattern, and discuss a transition from Kelvin to Mach angle.
Acoustic and elastic waves in metamaterials for underwater applications
NASA Astrophysics Data System (ADS)
Titovich, Alexey S.
Elastic effects in acoustic metamaterials are investigated. Water-based periodic arrays of elastic scatterers, sonic crystals, suffer from low transmission due to the impedance and index mismatch of typical engineering materials with water. A new type of acoustic metamaterial element is proposed that can be tuned to match the acoustic properties of water in the quasi-static regime. The element comprises a hollow elastic cylindrical shell fitted with an optimized internal substructure consisting of a central mass supported by an axisymmetric distribution of elastic stiffeners, which dictate the shell's effective bulk modulus and density. The derived closed form scattering solution for this system shows that the subsonic flexural waves excited in the shell by the attachment of stiffeners are suppressed by including a sufficiently large number of such stiffeners. As an example of refraction-based wave steering, a cylindrical-to-plane wave lens is designed by varying the bulk modulus in the array according to the conformal mapping of a unit circle to a square. Elastic shells provide rich scattering properties, mainly due to their ability to support highly dispersive flexural waves. Analysis of flexural-borne waves on a pair of shells yields an analytical expression for the width of a flexural resonance, which is then used with the theory of multiple scattering to accurately predict the splitting of the resonance frequency. This analysis leads to the discovery of the acoustic Poisson-like effect in a periodic wave medium. This effect redirects an incident acoustic wave by 90° in an otherwise acoustically transparent sonic crystal. An unresponsive "deaf" antisymmetric mode locked to band gap boundaries is unlocked by matching Bragg scattering with a quadrupole flexural resonance of the shell. The dynamic effect causes normal unidirectional wave motion to strongly couple to perpendicular motion, analogous to the quasi-static Poisson effect in solids. The Poisson
NASA Astrophysics Data System (ADS)
Perton, Mathieu; Contreras-Zazueta, Marcial A.; Sánchez-Sesma, Francisco J.
2016-04-01
A new implementation of IBEM allows simulating the elastic wave propagation in complex configurations made of embedded regions that are or homogeneous with irregular boundaries or flat layers. In an older implementation, each layer of a flat layered region would have been treated as a separated homogeneous region without taking into account the flat boundary information. For both types of regions, the scattered field results from fictitious sources positioned along their boundaries. For the homogeneous regions, the fictitious sources emit as in a full-space and the wave field is given by analytical Green's functions. For flat layered regions, fictitious sources emit as in an unbounded flat layered region and the wave field is given by Green's functions obtained from the Discrete Wave Number (DWN) method. The new implementation allows then reducing the length of the discretized boundaries but DWN Green's functions require much more computation time than the full space Green's functions. Several optimization steps are then implemented and commented. Validations are presented for 2D and 3D problems. Higher efficiency is achieved in 3D.
NASA Astrophysics Data System (ADS)
dot Zak, A.; Ostachowicz, W.; Krawczuk, M.
2011-07-01
Damage of aircraft structural elements in any form always present high risks. Failures of these elements can be caused by various reasons including material fatigue or impact leading to damage initiation and growth. Detection of these failures at their earliest stage of development, estimation of their size and location, are one of the most crucial factors for each damage detection method. Structural health monitoring strategies based on propagation of guided elastic waves in structures and wave interaction with damage related discontinuities are very promising tools that offer not only damage detection capabilities, but are also meant to provide precise information about the state of the structures and their remaining lifetime. Because of that various techniques are employed to simulate and mimic the wave-discontinuity interactions. The use of various types of sensors, their networks together with sophisticated contactless measuring techniques are investigated both numerically and experimentally. Certain results of numerical simulations obtained by the use of the spectral finite element method are presented by the authors and related with propagation of guided elastic waves in shell-type aircraft structures. Two types of structures are considered: flat 2D panels with or without stiffeners and 3D shell structures. The applicability of two different damage detection approaches is evaluated in order to detect and localise damage in these structures. Selected results related with the use of laser scanning vibrometry are also presented and discussed by the authors.
Scattering resonance of elastic wave and low-frequency equivalent slow wave
NASA Astrophysics Data System (ADS)
Meng, X.; Liu, H.; Hu, T.; Yang, L.
2015-12-01
Transmitted wave occurs as fast p-wave and slow p-wave in certain conditions when seismic waves travel through inhomogeneous layers. Energy of slow p-waves is strongest at some frequency band, but rather weak at both high frequency band and low frequency band, called scattering resonance. For practical seismic exploration, the frequency of slow p-wave occurs is below 10Hz, which cannot be explained by Biot's theory which predicts existence of the slow p-wave at ultrasonic band in the porous media. The slow p-wave equation have been derived, but which only adapted to explaining slow p-wave in the ultrasonic band. Experimental observations exhibit that slow p-wave also exists in nonporous media but with enormous low-velocity interbeds. When vertical incidence, elastic wave is simplified as compressing wave, the generation of slow waves is independent on shear wave. In the case of flat interbed and gas bubble, Liu (2006) has studied the transmission of acoustic waves, and found that the slow waves below the 10Hz frequency band can be explained. In the case of general elastic anisotropy medium, the tiheoretical research on the generation of slow waves is insufficient. Aiming at this problem, this paper presents an exponential mapping method based on transmitted wave (Magnus 1954), which can successfully explain the generation of the slow wave transmission in that case. Using the prediction operator (Claerbout 1985) to represent the transmission wave, this can be derived as first order partial differential equation. Using expansions in the frequency domain and the wave number domain, we find that the solutions have different expressions in the case of weak scattering and strong scattering. Besides, the method of combining the prediction operator and the exponential map is needed to extend to the elastic wave equation. Using the equation (Frazer and Fryer 1984, 1987), we derive the exponential mapping solution for the prediction operator of the general elastic medium
NASA Astrophysics Data System (ADS)
Uhrig, Matthias P.; Kim, Jin-Yeon; Jacobs, Laurence J.
2016-02-01
This research presents a 3D numerical finite element (FE) model which, previously developed, precisely simulates non-contact, air-coupled measurements of nonlinear Rayleigh wave propagation. The commercial FE-solver ABAQUS is used to perform the simulations. First, frequency dependent pressure wave attenuation is investigated numerically to reconstruct the sound pressure distribution along the active surface of the non-contact receiver. Second, constitutive law and excitation source properties are optimized to match nonlinear ultrasonic experimental data. Finally, the FE-model data are fit with analytical solutions showing a good agreement and thus, indicating the significance of the study performed.
NASA Astrophysics Data System (ADS)
Greenhalgh, Stewart; Zhou, Bing; Maurer, Hansruedi
2010-05-01
We have developed a modified version of the spectral element method (SEM), called the Gaussian Quadrature Grid (GQG) approach, for frequency domain 3D seismic modelling in arbitrary heterogeneous, anisotropic media. The model may incorporate an arbitrary free-surface topography and irregular subsurface interfaces. Unlike the SEM ,it does not require a powerful mesh generator such as the Delauney Triangular or TetGen. Rather, the GQG approach replaces the element mesh with Gaussian quadrature abscissae to directly sample the physical properties of the model parameters and compute the weighted residual or variational integral. This renders the model discretisation simple and easily matched to the model topography, as well as direct control of the model paramterisation for subsequent inversion. In addition, it offers high accuracy in numerical modelling provided that an appropriate density of the Gaussian quadrature abscissae is employed. The second innovation of the GQG is the incorporation of a new implementation of perfectly matched layers to suppress artificial reflections from the domain margins. We employ PML model parameters (specified complex valued density and elastic moduli) rather than explicitly solving the governing wave equation with a complex co-ordinate system as in conventional approaches. Such an implementation is simple, general, effective and easily extendable to any class of anisotropy and other numerical modelling methods. The accuracy of the GQG approach is controlled by the number of Gaussian quadrature points per minimum wavelength, the so-called sampling density. The optimal sampling density should be the one which enables high definition of geological characteristics and high precision of the variational integral evaluation and spatial differentiation. Our experiments show that satisfactory results can be obtained using sampling densities of 5 points per minimum wavelength. Efficiency of the GQG approach mainly depends on the linear
Wave-induced response of poro-elastic offshore foundations
Toha, F.X.
1983-01-01
A plane strain analysis based on Biot's theory of consolidation is utilized to investigate the pore-water pressures and displacements induced by steady-state linear planar waves beneath offshore gravity structures placed on poro-elastic seabed of finite thickness. The response is characterized by three controlling parameters, i.e., poroelasticity factor, fluid compressibility factor, and Poisson's Ratio. A parametric study, using a range of the controlling parameters encountered in practice, indicates that the response is governed primarily by the poro-elasticity factor and to a lesser extent by Poisson's Ratio (fluid compressibility factor has negligible influence for saturated soils). For offshore gravity structures, the wave-induced moment and horizontal force exerted by the structure on the foundation contribute most to the overall response while the influence of sea water pressure penetration through the sea bottom is, in general, negligible. For smaller structures, however, the latter influence dominates the response for soils as fine as silty sand. The maximum shear strain amplitudes evaluated imply that the soil response remains primarily in the linear elastic range, even for severe storm loadings. An experimental procedure was developed for the measurement of a composite poro-elastic property which combines all the pertinent soil properties. Two test configurations were found to be suitable depending on the hydraulic conductivity of the soil considered. The results of test on nine specimens reconstituted from three soil samples conform with the theoretical estimates and indicate that the developed procedure is viable.
Elastic Wave Propagation in Concrete and Continuous Wavelet Transform
Chiang, C.-H.; Gi, Y.-F.; Pan, C.-L.; Cheng, C.-C.
2005-04-09
Elastic wave methods, such as the ultrasonic pulse velocity and the impact echo, are often subject to multiple reflections at the boundaries of various constituents of concrete. Current study aims to improve the feature identification of elastic wave propagation due to buried objects in concrete slabs and cylinders. Embedded steel reinforcement, steel and PVC tubes, wooden disks, and rubber spheres are tested. The received signals are analyzed using continuous wavelet transform. As a result, signals are decomposed into distinctive frequency bands with transient information preserved. The interpretation of multiple reflections at different boundary conditions thus becomes more straightforward. Features related to reflections from steel bar, PVC tube, and steel tube can be readily identified in the magnitude plot of wavelet coefficients. Vibration modes of the concrete slab corresponding to different buried objects can also be separated based on corresponding time duration.
Elastic Wave Propagation in Concrete and Continuous Wavelet Transform
NASA Astrophysics Data System (ADS)
Chiang, Chih-Hung; Gi, Yu-Fung; Pan, Chi-Ling; Cheng, Chia-Chi
2005-04-01
Elastic wave methods, such as the ultrasonic pulse velocity and the impact echo, are often subject to multiple reflections at the boundaries of various constituents of concrete. Current study aims to improve the feature identification of elastic wave propagation due to buried objects in concrete slabs and cylinders. Embedded steel reinforcement, steel and PVC tubes, wooden disks, and rubber spheres are tested. The received signals are analyzed using continuous wavelet transform. As a result, signals are decomposed into distinctive frequency bands with transient information preserved. The interpretation of multiple reflections at different boundary conditions thus becomes more straightforward. Features related to reflections from steel bar, PVC tube, and steel tube can be readily identified in the magnitude plot of wavelet coefficients. Vibration modes of the concrete slab corresponding to different buried objects can also be separated based on corresponding time duration.
Decoupling Nonclassical Nonlinear Behavior of Elastic Wave Types
NASA Astrophysics Data System (ADS)
Remillieux, Marcel C.; Guyer, Robert A.; Payan, Cédric; Ulrich, T. J.
2016-03-01
In this Letter, the tensorial nature of the nonequilibrium dynamics in nonlinear mesoscopic elastic materials is evidenced via multimode resonance experiments. In these experiments the dynamic response, including the spatial variations of velocities and strains, is carefully monitored while the sample is vibrated in a purely longitudinal or a purely torsional mode. By analogy with the fact that such experiments can decouple the elements of the linear elastic tensor, we demonstrate that the parameters quantifying the nonequilibrium dynamics of the material differ substantially for a compressional wave and for a shear wave. This result could lead to further understanding of the nonlinear mechanical phenomena that arise in natural systems as well as to the design and engineering of nonlinear acoustic metamaterials.
Decoupling Nonclassical Nonlinear Behavior of Elastic Wave Types.
Remillieux, Marcel C; Guyer, Robert A; Payan, Cédric; Ulrich, T J
2016-03-18
In this Letter, the tensorial nature of the nonequilibrium dynamics in nonlinear mesoscopic elastic materials is evidenced via multimode resonance experiments. In these experiments the dynamic response, including the spatial variations of velocities and strains, is carefully monitored while the sample is vibrated in a purely longitudinal or a purely torsional mode. By analogy with the fact that such experiments can decouple the elements of the linear elastic tensor, we demonstrate that the parameters quantifying the nonequilibrium dynamics of the material differ substantially for a compressional wave and for a shear wave. This result could lead to further understanding of the nonlinear mechanical phenomena that arise in natural systems as well as to the design and engineering of nonlinear acoustic metamaterials. PMID:27035309
Prediction of Tsunami Waves and Runup Generated by 3d Granular Landslides
NASA Astrophysics Data System (ADS)
Mohammed, F.; Fritz, H. M.
2008-12-01
Subaerial and submarine landslides can trigger tsunamis with locally high amplitudes and runup, which can cause devastating effects in the near field region. The 50th anniversary of the Lituya Bay 1958 landslide impact generated mega tsunami recalls the largest tsunami runup of 524m in recorded history. In contrast to earthquake generated tsunamis, landslide generated tsunami sources are not confined to active tectonic regions and therefore are of particular importance for the Atlantic Ocean. Landslide generated tsunamis were studied in the three dimensional NEES tsunami wave basin TWB at OSU based on the generalized Froude similarity. A novel pneumatic landslide generator was deployed to control the landslide geometry and kinematics. Granular materials were used to model deformable landslides. Measurement techniques such as particle image velocimetry (PIV), multiple above and underwater video cameras, multiple acoustic transducer arrays (MTA), as well as resistance wave and runup gauges were applied. The wave generation was characterized by an extremely unsteady three phase flow consisting of the slide granulate, water and air entrained into the flow. The underwater cameras and the MTA provide data on the landslide deformation as it impacts the water surface, penetrates the water and finally deposits on the bottom of the basin. The influence of the landslide volume, shape and the impact speed on the generated tsunami wave characteristics were extensively studied. The experimental data provides prediction models for the generated tsunami wave characteristics based on the initial landslide characteristics and the final slide deposits. PIV provided instantaneous surface velocity vector fields, which gave insight into the kinematics of the landslide and wave generation process. At high impact velocities flow separation occurred on the slide shoulder resulting in a hydrodynamic impact crater. The recorded wave profiles yielded information on the wave propagation and
3D finite element modelling of guided wave scattering at delaminations in composites
NASA Astrophysics Data System (ADS)
Murat, Bibi Intan Suraya; Fromme, Paul
2016-02-01
Carbon fiber laminate composites are increasingly used for aerospace structures as they offer a number of advantages including a good strength to weight ratio. However, impact during the operation and servicing of the aircraft can lead to barely visible and difficult to detect damage. Depending on the severity of the impact, delaminations can occur, reducing the load carrying capacity of the structure. Efficient nondestructive testing of composite panels can be achieved using guided ultrasonic waves propagating along the structure. The guided wave (A0 Lamb wave mode) scattering at delaminations was modeled using full three-dimensional Finite Element (FE) simulations. The influence of the delamination size was systematically investigated from a parameter study. A significant influence of the delamination width on the guided wave scattering was found, especially on the angular dependency of the scattered guided wave amplitude. The sensitivity of guided ultrasonic waves for the detection of delamination damage in composite panels is discussed.
NASA Technical Reports Server (NTRS)
Silva, Walter A.; Sanetrik, Mark D.; Chwalowski, Pawel; Connolly, Joseph; Kopasakis, George
2016-01-01
An overview of recent applications of the FUN3D CFD code to computational aeroelastic, sonic boom, and aeropropulsoservoelasticity (APSE) analyses of a low-boom supersonic configuration is presented. The overview includes details of the computational models developed including multiple unstructured CFD grids suitable for aeroelastic and sonic boom analyses. In addition, aeroelastic Reduced-Order Models (ROMs) are generated and used to rapidly compute the aeroelastic response and utter boundaries at multiple flight conditions.
Moulding and shielding flexural waves in elastic plates
NASA Astrophysics Data System (ADS)
Antonakakis, T.; Craster, R. V.; Guenneau, S.
2014-03-01
Platonic crystals (PlCs) are the elastic plate analogue of the photonic crystals widely used in optics, and are thin structured elastic plates along which flexural waves cannot propagate within certain stop band frequency intervals. The practical importance of PlCs is twofold: These can be used either in the design of microstructured acoustic metamaterials or as an approximate model for surface elastic waves propagating in meter scale seismic metamaterials. Here, we make use of the band spectrum of PlCs created by an array of either very small or densely packed clamped circles to achieve surface wave reflectors at very large wavelengths, a flat lens, a waveguide effect, a directive antenna near the stop band frequencies. The limit in which the circles reduce to points is particularly appealing as there is an exact dispersion relation available so the origin of these phenomena can be explained and interpreted using Fourier series and high-frequency homogenization (HFH). We then enlarge the radius of clamped circles, which both makes the zero-frequency stop band up to five times wider and flattens the dispersion curves. Here, HFH notably captures the essence of localized modes, one of which appears in the zero-frequency stop band and is used in the design of a highly directive waveguide.
Making and Propagating Elastic Waves: Overview of the new wave propagation code WPP
McCandless, K P; Petersson, N A; Nilsson, S; Rodgers, A; Sjogreen, B; Blair, S C
2006-05-09
We are developing a new parallel 3D wave propagation code at LLNL called WPP (Wave Propagation Program). WPP is being designed to incorporate the latest developments in embedded boundary and mesh refinement technology for finite difference methods, as well as having an efficient portable implementation to run on the latest supercomputers at LLNL. We are currently exploring seismic wave applications, including a recent effort to compute ground motions for the 1906 Great San Francisco Earthquake. This paper will briefly describe the wave propagation problem, features of our numerical method to model it, implementation of the wave propagation code, and results from the 1906 Great San Francisco Earthquake simulation.
Morita, Yasuyuki; Kawase, Naoki; Ju, Yang; Yamauchi, Takashi
2016-07-01
Cells maintain homeostasis and perform various functions by interacting mechanically with a cell-adhesive matrix. Regarding cellular differentiation, it has been found that matrix elasticity can determine the differentiation lineage of mesenchymal stem cells (MSCs). Direct quantitative measurements of the mechanical interaction between MSCs and matrix for differentiation, however, have yet to be reported. Herein, the displacement field of the cell-adhesive matrix was observed quantitatively using a digital volume correlation (DVC) method. Maximum displacement and cellular traction stress were analyzed when the MSC differentiated into a neuron-like cell or an osteoblast-like cell on a soft or hard elastic matrix, respectively. The function of non-muscle myosin II (NMM II), which plays an important role in intracellular cytoskeletal dynamics, was investigated during cellular differentiation. The mechanical interaction (maximum displacement and subjected area of the matrix) between the cell and matrix was dependent on matrix elasticity. It has also been shown that the mechanical interaction between the intracellular cytoskeleton and cell-adhesion matrix is indispensable for cellular differentiation. This work provides the first quantitative visualization of the mechanical interaction between MSCs and the cell-adhesion matrix for differentiation. PMID:26945874
WaveQ3D: Fast and accurate acoustic transmission loss (TL) eigenrays, in littoral environments
NASA Astrophysics Data System (ADS)
Reilly, Sean M.
This study defines a new 3D Gaussian ray bundling acoustic transmission loss model in geodetic coordinates: latitude, longitude, and altitude. This approach is designed to lower the computation burden of computing accurate environmental effects in sonar training application by eliminating the need to transform the ocean environment into a collection of Nx2D Cartesian radials. This approach also improves model accuracy by incorporating real world 3D effects, like horizontal refraction, into the model. This study starts with derivations for a 3D variant of Gaussian ray bundles in this coordinate system. To verify the accuracy of this approach, acoustic propagation predictions of transmission loss, time of arrival, and propagation direction are compared to analytic solutions and other models. To validate the model's ability to predict real world phenomena, predictions of transmission loss and propagation direction are compared to at-sea measurements, in an environment where strong horizontal refraction effect have been observed. This model has been integrated into U.S. Navy active sonar training system applications, where testing has demonstrated its ability to improve transmission loss calculation speed without sacrificing accuracy.
Analysis of non linear partially standing waves from 3D velocity measurements
NASA Astrophysics Data System (ADS)
Drevard, D.; Rey, V.; Svendsen, Ib; Fraunie, P.
2003-04-01
Surface gravity waves in the ocean exhibit an energy spectrum distributed in both frequency and direction of propagation. Wave data collection is of great importance in coastal zones for engineering and scientific studies. In particular, partially standing waves measurements near coastal structures and steep or barred beaches may be a requirement, for instance for morphodynamic studies. The aim of the present study is the analysis of partially standing surface waves icluding non-linear effects. According to 1st order Stokes theory, synchronous measurements of horizontal and vertical velocity components allow calculation of rate of standing waves (Drevard et al, 2003). In the present study, it is demonstrated that for deep water conditions, partially standing 2nd order Stokes waves induced velocity field is still represented by the 1st order solution for the velocity potential contrary to the surface elevation which exhibits harmonic components. For intermediate water depth, harmonic components appear not only in the surface elevation but also in the velocity fields, but their weight remains much smaller, because of the vertical decreasing wave induced motion. For irregular waves, the influence of the spectrum width on the non-linear effects in the analysis is discussed. Keywords: Wave measurements ; reflection ; non-linear effects Acknowledgements: This work was initiated during the stay of Prof. Ib Svendsen, as invited Professor, at LSEET in autumn 2002. This study is carried out in the framework of the Scientific French National Programmes PNEC ART7 and PATOM. Their financial supports are acknowledged References: Drevard, D., Meuret, A., Rey, V. Piazzola, J. And Dolle, A.. (2002). "Partially reflected waves measurements using Acoustic Doppler Velocimeter (ADV)", Submitted to ISOPE 03, Honolulu, Hawaii, May 2003.
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
NASA Astrophysics Data System (ADS)
Deng, Y.; Ebert-Uphoff, I.; Chen, J.
2015-12-01
Causal discovery seeks to discover potential cause-effect relationships from observational data. Here we adopt the idea of interpreting large-scale atmospheric dynamical processes, particularly those tied to propagation of large-scale waves, as information flow around the globe, which can then be calculated using causal discovery methods. We apply a well-established causal discovery algorithm - based on constraint-based structure learning of probabilistic graphical models - toward 51 years of 6-hourly, atmospheric isobaric-level geopotential height data to construct the first-ever graphs of 3D information flow in the atmosphere. These graphs are created globally for different seasons and their connection to phase/energy propagation of atmospheric waves are investigated. Specifically, we examine the information flows 1) in the topical region that represent horizontal and vertical propagations of Kelvin and Rossby-gravity waves whose associated momentum transfer are known to play a key role in the Quasi-Biennial Oscillation (QBO), and 2) in the northern extratropics that represent propagations of planetary-scale waves whose heat/momentum fluxes are responsible for vacillations in the polar stratospheric vortex and occurrences of extreme events such as the stratospheric sudden warming. The sensitivity of the constructed graphs of 3D information flow to data resolution and pre-processing methods (e.g., spatial and temporal filtering) will be discussed.
Wang, Tao; Green, Ryan; Nair, Rajesh Ramakrishnan; Howell, Mark; Mohapatra, Subhra; Guldiken, Rasim; Mohapatra, Shyam Sundar
2015-01-01
Detection and quantification of cell viability and growth in two-dimensional (2D) and three-dimensional (3D) cell cultures commonly involve harvesting of cells and therefore requires a parallel set-up of several replicates for time-lapse or dose–response studies. Thus, developing a non-invasive and touch-free detection of cell growth in longitudinal studies of 3D tumor spheroid cultures or of stem cell regeneration remains a major unmet need. Since surface acoustic waves (SAWs) permit mass loading-based biosensing and have been touted due to their many advantages including low cost, small size and ease of assembly, we examined the potential of SAW-biosensing to detect and quantify cell growth. Herein, we demonstrate that a shear horizontal-surface acoustic waves (SH-SAW) device comprising two pairs of resonators consisting of interdigital transducers and reflecting fingers can be used to quantify mass loading by the cells in suspension as well as within a 3D cell culture platform. A 3D COMSOL model was built to simulate the mass loading response of increasing concentrations of cells in suspension in the polydimethylsiloxane (PDMS) well in order to predict the characteristics and optimize the design of the SH-SAW biosensor. The simulated relative frequency shift from the two oscillatory circuit systems (one of which functions as control) were found to be concordant to experimental data generated with RAW264.7 macrophage and A549 cancer cells. In addition, results showed that SAW measurements per se did not affect viability of cells. Further, SH-SAW biosensing was applied to A549 cells cultured on a 3D electrospun nanofiber scaffold that generate tumor spheroids (tumoroids) and the results showed the device's ability to detect changes in tumor spheroid growth over the course of eight days. Taken together, these results demonstrate the use of SH-SAW device for detection and quantification of cell growth changes over time in 2D suspension cultures and in 3D cell
Wang, Tao; Green, Ryan; Nair, Rajesh Ramakrishnan; Howell, Mark; Mohapatra, Subhra; Guldiken, Rasim; Mohapatra, Shyam Sundar
2015-01-01
Detection and quantification of cell viability and growth in two-dimensional (2D) and three-dimensional (3D) cell cultures commonly involve harvesting of cells and therefore requires a parallel set-up of several replicates for time-lapse or dose-response studies. Thus, developing a non-invasive and touch-free detection of cell growth in longitudinal studies of 3D tumor spheroid cultures or of stem cell regeneration remains a major unmet need. Since surface acoustic waves (SAWs) permit mass loading-based biosensing and have been touted due to their many advantages including low cost, small size and ease of assembly, we examined the potential of SAW-biosensing to detect and quantify cell growth. Herein, we demonstrate that a shear horizontal-surface acoustic waves (SH-SAW) device comprising two pairs of resonators consisting of interdigital transducers and reflecting fingers can be used to quantify mass loading by the cells in suspension as well as within a 3D cell culture platform. A 3D COMSOL model was built to simulate the mass loading response of increasing concentrations of cells in suspension in the polydimethylsiloxane (PDMS) well in order to predict the characteristics and optimize the design of the SH-SAW biosensor. The simulated relative frequency shift from the two oscillatory circuit systems (one of which functions as control) were found to be concordant to experimental data generated with RAW264.7 macrophage and A549 cancer cells. In addition, results showed that SAW measurements per se did not affect viability of cells. Further, SH-SAW biosensing was applied to A549 cells cultured on a 3D electrospun nanofiber scaffold that generate tumor spheroids (tumoroids) and the results showed the device's ability to detect changes in tumor spheroid growth over the course of eight days. Taken together, these results demonstrate the use of SH-SAW device for detection and quantification of cell growth changes over time in 2D suspension cultures and in 3D cell
NASA Astrophysics Data System (ADS)
Soboleva, O. N.; Kurochkina, E. P.
2016-01-01
The effective coefficients in the problem of the acoustic wave propagation have been calculated for a multiscale 3D isotropic medium using a subgrid modeling approach. The density and the elastic stiffness have been represented mathematically by the Kolmogorov multiplicative cascades, which, to date, appear to be the only mechanisms for generating a stationary multifractal fields with a log-stable probability distribution. The fields with the stable distribution are described with the help of linear combination random values ?, ? and weight coefficients ?, ?, which satisfy certain conditions in the nodes of spatial grid ?. The parameters of the stable distribution of the random values ?, ? are equal: ?, ?, ?, ?. The wavelength is assumed to be large as compared with the scale of heterogeneities of the medium. We consider the regime in which the waves propagate over a distance of the typical wave length in source. The theoretical results obtained in this paper are compared with the results of a direct 3D numerical simulation.
NASA Astrophysics Data System (ADS)
Moritake, Hiroshi; Seike, Tadaaki; Toda, Kohji
1999-05-01
An elastic wave delay line with a glass plate as a propagation medium is investigated in relations to acoustooptic effects of nematic liquid crystals. A nematic liquid crystal cell is mounted on the central region of a glass plate for elastic wave propagation. The elastic wave propagating in the glass plate interacts with the liquid crystal in a wide frequency range. Two types of periodical domain structures are observed in the nematic liquid crystal cell under the existence of the elastic wave. One exists in both homeotropically and homogeneously aligned cells, and depends not on the kind of liquid crystal but on the carrier frequency of the elastic wave. The other is recognized only in homogeneously aligned cells and depends on the layer thickness and the kind of liquid crystal, but not on the carrier frequency of the elastic wave. Both periodical domain structures are induced by the elastic wave propagating in the glass plate.
Vdovin V.L.
2005-08-15
In this report we describe theory and 3D full wave code description for the wave excitation, propagation and absorption in 3-dimensional (3D) stellarator equilibrium high beta plasma in ion cyclotron frequency range (ICRF). This theory forms a basis for a 3D code creation, urgently needed for the ICRF heating scenarios development for the operated LHD, constructed W7-X, NCSX and projected CSX3 stellarators, as well for re evaluation of ICRF scenarios in operated tokamaks and in the ITER . The theory solves the 3D Maxwell-Vlasov antenna-plasma-conducting shell boundary value problem in the non-orthogonal flux coordinates ({Psi}, {theta}, {var_phi}), {Psi} being magnetic flux function, {theta} and {var_phi} being the poloidal and toroidal angles, respectively. All basic physics, like wave refraction, reflection and diffraction are self consistently included, along with the fundamental ion and ion minority cyclotron resonances, two ion hybrid resonance, electron Landau and TTMP absorption. Antenna reactive impedance and loading resistance are also calculated and urgently needed for an antenna -generator matching. This is accomplished in a real confining magnetic field being varying in a plasma major radius direction, in toroidal and poloidal directions, through making use of the hot dense plasma wave induced currents with account to the finite Larmor radius effects. We expand the solution in Fourier series over the toroidal ({var_phi}) and poloidal ({theta}) angles and solve resulting ordinary differential equations in a radial like {Psi}-coordinate by finite difference method. The constructed discretization scheme is divergent-free one, thus retaining the basic properties of original equations. The Fourier expansion over the angle coordinates has given to us the possibility to correctly construct the ''parallel'' wave number k{sub //}, and thereby to correctly describe the ICRF waves absorption by a hot plasma. The toroidal harmonics are tightly coupled with each
NASA Astrophysics Data System (ADS)
Perton, Mathieu; Contreras-Zazueta, Marcial A.; Sánchez-Sesma, Francisco J.
2016-06-01
A new implementation of indirect boundary element method allows simulating the elastic wave propagation in complex configurations made of embedded regions that are homogeneous with irregular boundaries or flat layers. In an older implementation, each layer of a flat layered region would have been treated as a separated homogeneous region without taking into account the flat boundary information. For both types of regions, the scattered field results from fictitious sources positioned along their boundaries. For the homogeneous regions, the fictitious sources emit as in a full-space and the wave field is given by analytical Green's functions. For flat layered regions, fictitious sources emit as in an unbounded flat layered region and the wave field is given by Green's functions obtained from the discrete wavenumber (DWN) method. The new implementation allows then reducing the length of the discretized boundaries but DWN Green's functions require much more computation time than the full-space Green's functions. Several optimization steps are then implemented and commented. Validations are presented for 2-D and 3-D problems. Higher efficiency is achieved in 3-D.
NASA Astrophysics Data System (ADS)
Sidler, Rolf; Carcione, José M.; Holliger, Klaus
2014-02-01
We present a novel approach for the comprehensive, flexible and accurate simulation of poroelastic wave propagation in 3-D cylindrical coordinates. An important application of this method is the realistic modelling of complex seismic wave phenomena in fluid-filled boreholes, which represents a major, as of yet largely unresolved, problem in exploration geophysics. To this end, we consider a numerical mesh consisting of three concentric domains representing the borehole fluid in the centre followed by the mudcake and/or casing, and the surrounding porous formation. The spatial discretization is based on a Chebyshev expansion in the radial direction and Fourier expansions in the vertical and azimuthal directions as well as a Runge-Kutta integration scheme for the time evolution. Trigonometric interpolation and a domain decomposition method based on the method of characteristics are used to match the boundary conditions at the fluid/porous-solid and porous-solid/porous-solid interfaces as well as to reduce the number of gridpoints in the innermost domain for computational efficiency. We apply this novel modelling approach to the particularly challenging scenario of near-surface borehole environments. To this end, we compare 3-D heterogeneous and corresponding rotationally invariant simulations, assess the sensitivity of Stoneley waves to formation permeability in the presence of a casing and evaluate the effects of an excavation damage zone behind a casing on sonic log recordings. Our results indicate that only first arrival times of fast modes are reasonably well described by rotationally invariant approximations of 3-D heterogenous media. We also find that Stoneley waves are indeed remarkably sensitive to the average permeability behind a perforated PVC casing, and that the presence of an excavation damage zone behind a casing tends to dominate the overall signature of recorded seismograms.
Rayleigh scattering and nonlinear inversion of elastic waves
Gritto, R.
1995-12-01
Rayleigh scattering of elastic waves by an inclusion is investigated and the limitations determined. In the near field of the inhomogeneity, the scattered waves are up to a factor of 300 stronger than in the far field, excluding the application of the far field Rayleigh approximation for this range. The investigation of the relative error as a function of parameter perturbation shows a range of applicability broader than previously assumed, with errors of 37% and 17% for perturbations of {minus}100% and +100%, respectively. The validity range for the Rayleigh limit is controlled by large inequalities, and therefore, the exact limit is determined as a function of various parameter configurations, resulting in surprisingly high values of up to k{sub p}R = 0.9. The nonlinear scattering problem can be solved by inverting for equivalent source terms (moments) of the scatterer, before the elastic parameters are determined. The nonlinear dependence between the moments and the elastic parameters reveals a strong asymmetry around the origin, which will produce different results for weak scattering approximations depending on the sign of the anomaly. Numerical modeling of cross hole situations shows that near field terms are important to yield correct estimates of the inhomogeneities in the vicinity of the receivers, while a few well positioned sources and receivers considerably increase the angular coverage, and thus the model resolution of the inversion parameters. The pattern of scattered energy by an inhomogeneity is complicated and varies depending on the object, the wavelength of the incident wave, and the elastic parameters involved. Therefore, it is necessary to investigate the direction of scattered amplitudes to determine the best survey geometry.
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.
Non-linear interaction of elastic waves in rocks
NASA Astrophysics Data System (ADS)
Kuvshinov, B. N.; Smit, T. J. H.; Campman, X. H.
2013-09-01
We study theoretically the interaction of elastic waves caused by non-linearities of rock elastic moduli, and assess the possibility to use this phenomenon in hydrocarbon exploration and in the analysis of rock samples. In our calculations we use the five-constant model by Gol'dberg. It is shown that the interaction of plane waves in isotropic solids is completely described by five coupling coefficients, which have the same order of magnitude. By considering scattering of compressional waves generated by controlled sources at the Earth surface from a non-linear layer at the subsurface, we conclude that non-linear signals from deep formations are unlikely to be measured with the current level of technology. Our analysis of field tests where non-linear signals were measured, suggests that these signals are generated either in the shallow subsurface or in the vicinity of sources. Non-linear wave interaction might be observable in lab tests with focused ultrasonic beams. In this case, the non-linear response is generated in the secondary parametric array formed by linear beams scattered from inclusions. Although the strength of this response is controlled by non-linearity of the surrounding medium rather than by non-linearity of inclusions, its measurement can help to obtain better images of rock samples.
Skin-friction measurements in a 3-D, supersonic shock-wave/boundary-layer interaction
NASA Astrophysics Data System (ADS)
Wideman, Jeffrey Kenneth
An experimental study has been conducted in a three-dimensional, supersonic shockwave/boundary-layer interaction (3-D SW/BLI) with the intent of providing accurate experimental data for turbulence modeling and computational fluid dynamics (CFD) code validation. The experiment was performed in the High Reynolds Channel 1 (HRCI) wind tunnel at NASA Ames Research Center. The test was conducted at a Mach number of M(sub infinity) = 2.89 and at a Reynolds number of Re = 15 x 106/m. The model consisted of a sting-supported cylinder aligned with the tunnel axis and a 20 deg half-angle conical flare offset 1.27 cm from the cylinder centerline. The generated shock system was verified to be steady by schlieren visualization. The highlight of the study was the acquisition of 3-D skin-friction data by a laser interferometric skin friction (LISF) meter. Surface pressure measurements were obtained in 15 deg intervals around the cylinder and flare. Additional measurements included surface oil flow and laser light sheet illumination which were used to document the flow topology. Skin-friction measurements are proving to be a very challenging test of a CFD code predictive capability. However, at the present time there is a very limited amount of accurate skin-friction data in complex flows such as in 3-D SW/BLI. The LISF technique is advantageous as compared to other skin-friction measurement techniques for application in complex flows like the present since it is non-intrusive and is capable of performing measurements in flows with large shear and pressure gradients where the reliability of other techniques is questionable. Thus, the prevent skin-friction data will prove valuable to turbulence modeling and CFD code validation efforts.
Propagation of 3D nonlinear waves over complex bathymetry using a High-Order Spectral method
NASA Astrophysics Data System (ADS)
Gouin, Maïté; Ducrozet, Guillaume; Ferrant, Pierre
2016-04-01
Scattering of regular and irregular surface gravity waves propagating over a region of arbitrary three-dimensional varying bathymetry is considered here. The three-dimensional High-Order Spectral method (HOS) with an extension to account for a variable bathymetry is used. The efficiency of the model has been proved to be conserved even with this extension. The method is first applied to a bathymetry consisting of an elliptical lens, as used in the Vincent and Briggs (1989) experiment. Incident waves passing across the lens are transformed and a strong convergence region is observed after the elliptical mound. The wave amplification depends on the incident wave. Numerical results for regular and irregular waves are analysed and compared with other methods and experimental data demonstrating the efficiency and practical applicability of the present approach. Then the method is used to model waves propagating over a real bathymetry: the canyons of Scripps/La Jolla in California. The implementation of this complex bathymetry in the model is presented, as well as the first results achieved. They will be compared to the ones obtained with another numerical model.
NASA Astrophysics Data System (ADS)
Bennis, A.; Ardhuin, F.; Dumas, F.; Bonneton, P.
2010-12-01
The interaction of waves with three-dimensional current structure is investigated using a two-way coupled modelling system combining MARS3D (Lazure and Dumas 2008) with WAVEWATCH III (Tolman 2008, Ardhuin et al. 2009) , a wave model (NOAA/NCEP, Tolman 2008). After a basic validation in two dimensions, the flow model MARS3D was adapted with three options that solve for the total momentum (Mellor 2003, 2008) or the quasi-Eulerian momentum (Ardhuin et al. 2008b). Adiabatic model results show that, as expected from theory (Ardhuin et al. 2008a), the total momentum fluxes parameterized by Mellor are not self-consistent and can lead to very large errors (Bennis and Ardhuin 2010). We thus use the model option to solve for the quasi-Eulerian momentum, including sources of momentum and turbulent kinetic energy (TKE). The influence of these TKE sources is investigated in the case of the NSTS experiment (Thornton and Guza, 1986). The feedback of the currents on the waves is negligible in this case. The sources of TKE from wave breaking and wave bottom friction are found to have strong influence on the bottom friction, in a way consistent with the parameterizations by Longuet-Higgins (1970) and Mellor (2002). The complete model is then applied to a real case of a large rip current on the South-West coast of France (Bruneau et al., Cont. Shelf Res. 2009). The breaking of waves on the opposed current generates a strong coupling on the rip current that partially controls the strength of the current and it three-dimensional shape.
Conducting a 3D Converted Shear Wave Project to Reduce Exploration Risk at Wister, CA
Matlick, Skip; Walsh, Patrick; Rhodes, Greg; Fercho, Steven
2015-06-30
Ormat sited 2 full-size exploration wells based on 3D seismic interpretation of fractures, prior drilling results, and temperature anomaly. The wells indicated commercial temperatures (>300 F), but almost no permeability, despite one of the wells being drilled within 820 ft of an older exploration well with reported indications of permeability. Following completion of the second well in 2012, Ormat undertook a lengthy program to 1) evaluate the lack of observed permeability, 2) estimate the likelihood of finding permeability with additional drilling, and 3) estimate resource size based on an anticipated extent of permeability.
Display depth analyses with the wave aberration for the auto-stereoscopic 3D display
NASA Astrophysics Data System (ADS)
Gao, Xin; Sang, Xinzhu; Yu, Xunbo; Chen, Duo; Chen, Zhidong; Zhang, Wanlu; Yan, Binbin; Yuan, Jinhui; Wang, Kuiru; Yu, Chongxiu; Dou, Wenhua; Xiao, Liquan
2016-07-01
Because the aberration severely affects the display performances of the auto-stereoscopic 3D display, the diffraction theory is used to analyze the diffraction field distribution and the display depth through aberration analysis. Based on the proposed method, the display depth of central and marginal reconstructed images is discussed. The experimental results agree with the theoretical analyses. Increasing the viewing distance or decreasing the lens aperture can improve the display depth. Different viewing distances and the LCD with two lens-arrays are used to verify the conclusion.
Bending analysis of a general cross-ply laminate using 3D elasticity solution and layerwise theory
NASA Astrophysics Data System (ADS)
Yazdani Sarvestani, H.; Naghashpour, A.; Heidari-Rarani, M.
2015-12-01
In this study, the analytical solution of interlaminar stresses near the free edges of a general (symmetric and unsymmetric layups) cross-ply composite laminate subjected to pure bending loading is presented based on Reddy's layerwise theory (LWT) for the first time. First, the reduced form of displacement field is obtained for a general cross-ply composite laminate subjected to a bending moment by elasticity theory. Then, first-order shear deformation theory of plates and LWT is utilized to determine the global and local deformation parameters appearing in the displacement fields, respectively. One of the main advantages of the developed solution based on the LWT is exact prediction of interlaminar stresses at the boundary layer regions. To show the accuracy of this solution, three-dimensional elasticity bending problem of a laminated composite is solved for special set of boundary conditions as well. Finally, LWT results are presented for edge-effect problems of several symmetric and unsymmetric cross-ply laminates under the bending moment. The obtained results indicate high stress gradients of interlaminar stresses near the edges of laminates.
Elastic wave scattering to characterize heterogeneities in the borehole environment
NASA Astrophysics Data System (ADS)
Tang, Xiao-Ming; Li, Zhen; Hei, Chuang; Su, Yuan-Da
2016-04-01
Scattering due to small-scale heterogeneities in the rock formation surrounding a wellbore can significantly change the acoustic waveform from a logging measurement which in turn can be used to characterize the formation heterogeneities. This study simulates the elastic heterogeneity scattering in monopole and dipole acoustic logging and analyse the resulting effects on the waveforms. The results show that significant coda waves are generated in both monopole and dipole waveforms and the dipole coda is dominated by S-to-S scattering, which can be effectively utilized to diagnose the heterogeneity in the rock formation. The coda wave modelling and analysis were used to characterize dipole acoustic data logged before and after fracturing a reservoir interval, with significant coda wave in the after-fracturing data indicating fracturing-induced heterogeneous property change in the rock volume surrounding the borehole.
STUDY ON DISPERSION OF LONGITUDINAL ELASTIC WAVES IN ROCK SPECIMENS
NASA Astrophysics Data System (ADS)
Nishiyama, Satoshi; Ohnishi, Yuzo; Yano, Takao; Takahashi, Manabu; Yoshimura, Kimitaka; Ando, Ken-Ichi
In order to estimate the permeable properties and pore fluid properties in the ground, survey methods using dispersion and attenuation of elastic waves have been developed. The Biot theory have been applied to the frequency-dependent dispersion data, but the influence of geological properties on velocity dispersion and the relation between the permeable properties and velocity dispersion in rock specimens have not become clear. Sedimentary rock and granite specimens were tested using longitudinal waves for the difference of velocity dispersion phenomena observed in each specimen, and we examine whether the Biot theory can be applied to observed experimental data. Moreover, we tried to estimate the permeability of the rock specimen based on the theory and show that the Biot-squirt theory can be applied to estimating rock permeability using the seismic wave dispersion characteristics.
SAFE-3D analysis of a piezoelectric transducer to excite guided waves in a rail web
NASA Astrophysics Data System (ADS)
Ramatlo, Dineo A.; Long, Craig S.; Loveday, Philip W.; Wilke, Daniel N.
2016-02-01
Our existing Ultrasonic Broken Rail Detection system detects complete breaks and primarily uses a propagating mode with energy concentrated in the head of the rail. Previous experimental studies have demonstrated that a mode with energy concentrated in the head of the rail, is capable of detecting weld reflections at long distances. Exploiting a mode with energy concentrated in the web of the rail would allow us to effectively detect defects in the web of the rail and could also help to distinguish between reflections from welds and cracks. In this paper, we will demonstrate the analysis of a piezoelectric transducer attached to the rail web. The forced response at different frequencies is computed by the Semi-Analytical Finite Element (SAFE) method and compared to a full three-dimensional finite element method using ABAQUS. The SAFE method only requires the rail track cross-section to be meshed using two-dimensional elements. The ABAQUS model in turn requires a full three-dimensional discretisation of the rail track. The SAFE approach can yield poor predictions at cut-on frequencies associated with other modes in the rail. Problematic frequencies are identified and a suitable frequency range identified for transducer design. The forced response results of the two methods were found to be in good agreement with each other. We then use a previously developed SAFE-3D method to analyse a practical transducer over the selected frequency range. The results obtained from the SAFE-3D method are in good agreement with experimental measurements.
Nonlinear Evolution of a 3D Inertial Alfvén Wave and Its Implication in Particle Acceleration
NASA Astrophysics Data System (ADS)
Sharma, Prachi; Yadav, Nitin; Sharma, R. P.
2016-03-01
A simulation based on a pseudo-spectral method has been performed in order to study particle acceleration. A model for the acceleration of charged particles by field localization is developed for the low-β plasma. For this purpose, a fractional diffusion approach has been employed. The nonlinear interaction between a 3D inertial Alfvén wave and a slow magnetosonic wave has been examined, and the dynamical equations of these two waves in the presence of ponderomotive nonlinearity have been solved numerically. The nonlinear evolution of the inertial Alfvén wave in the presence of slow magnetosonic wave undergoes a filamentation instability and results in field intensity localization. The results obtained show the localization and power spectrum of inertial Alfvén wave due to nonlinear coupling. The scaling obtained after the first break point of the magnetic power spectrum has been used to calculate the formation of the thermal tail of energetic particles in the solar corona.
Rigas, Fotis; Sklavounos, Spyros
2005-05-20
Accidental blast wave generation and propagation in the surroundings poses severe threats for people and property. The prediction of overpressure maxima and its change with time at specified distances can lead to useful conclusions in quantitative risk analysis applications. In this paper, the use of a computational fluid dynamics (CFD) code CFX-5.6 on dense explosive detonation events is described. The work deals with the three-dimensional simulation of overpressure wave propagation generated by the detonation of a dense explosive within a small-scale branched tunnel. It also aids at validating the code against published experimental data as well as to study the way that the resulting shock wave propagates in a confined space configuration. Predicted overpressure histories were plotted and compared versus experimental measurements showing a reasonably good agreement. Overpressure maxima and corresponding times were found close to the measured ones confirming that CFDs may constitute a useful tool in explosion hazard assessment procedures. Moreover, it was found that blast wave propagates preserving supersonic speed along the tunnel accompanied by high overpressure levels, and indicating that space confinement favors the formation and maintenance of a shock rather than a weak pressure wave. PMID:15885402
Effect of spatial dispersion on transient acoustic wave propagation in 3D.
Every, A G
2006-12-22
Spatial dispersion is the variation of wave speed with wavelength. It sets in when the acoustic wavelength approaches the natural scale of length of the medium, which could, for example, be the lattice constant of a crystal, the repeat distance in a superlattice, or the grain size in a granular material. In centrosymmetric media, the first onset of dispersion is accommodated by the introduction of fourth order spatial derivatives into the wave equation. These lead to a correction to the phase velocity which is quadratic in the spatial frequency. This paper treats the effect of spatial dispersion on the point force elastodynamic Green's functions of solids. The effects of dispersion are shown to be most pronounced in the vicinity of wave arrivals. These lose their singular form, and are transformed into wave trains known as quasi-arrivals. The step and ramp function wave arrivals are treated, and it is shown that their unfolded quasi-arrival forms can be expressed in terms of integrals involving the Airy function. PMID:16828830
NASA Astrophysics Data System (ADS)
Schuberth, B.; Piazzoni, A.; Bunge, H.; Igel, H.; Steinle-Neumann, G.; Moder, C.; Oeser, J.
2007-12-01
Our current understanding of mantle structure and dynamics is to a large part based on inversion of seismic data resulting in tomographic images and on direct analysis of a wide range of seismic phases such as Pdiff, PcP, ScS SdS etc. For solving inverse problems, forward modeling is needed to obtain a synthetic dataset for a given set of model parameters. In this respect, great progress has been made over the last years in the developement of sophisticated numerical full waveform modeling tools. However, the main limitation in the application of this new class of techniques for the forward problem of seismology is the lack of accurate predictions of mantle heterogeneity that allow us to test hypotheses about Earth's mantle. Such predictive models should be based on geodynamic and mineralogical considerations and derived independently of seismological observations. Here, we demonstrate the feasibility of joining forward simulations from geodynamics, mineral physics and seismology to obtain earth-like seismograms. 3D global wave propagation is simulated for dynamically consistent thermal structures derived from 3D mantle circulation modeling (e.g. Bunge et al. 2002), for which the temperatures are converted to seismic velocities using a recently published, thermodynamically self-consistent mineral physics approach (Piazzoni et al. 2007). Assuming a certain, fixed mantle composition (e.g. pyrolite) our mineralogic modeling algorithm computes the stable phases at mantle pressures for a wide range of temperatures by system Gibbs free energy minimization. Through the same equations of state that model the Gibbs free energy, we compute elastic moduli and density for each stable phase assemblage at the same P-T conditions. One straightforward application of this approach is the study of the seismic signature of synthetic mantle discontinuities arising in such models, as the temperature dependent phase transformations occuring at around 410 Km and 660 Km depth are
NASA Astrophysics Data System (ADS)
Vidal, A.; San-Blas, A. A.; Quesada-Pereira, F. D.; Pérez-Soler, J.; Gil, J.; Vicente, C.; Gimeno, B.; Boria, V. E.
2015-07-01
A novel technique for the full-wave analysis of 3-D complex waveguide devices is presented. This new formulation, based on the Boundary Integral-Resonant Mode Expansion (BI-RME) method, allows the rigorous full-wave electromagnetic characterization of 3-D arbitrarily shaped metallic structures making use of extremely low CPU resources (both time and memory). The unknown electric current density on the surface of the metallic elements is represented by means of Rao-Wilton-Glisson basis functions, and an algebraic procedure based on a singular value decomposition is applied to transform such functions into the classical solenoidal and nonsolenoidal basis functions needed by the original BI-RME technique. The developed tool also provides an accurate computation of the electromagnetic fields at an arbitrary observation point of the considered device, so it can be used for predicting high-power breakdown phenomena. In order to validate the accuracy and efficiency of this novel approach, several new designs of band-pass waveguides filters are presented. The obtained results (S-parameters and electromagnetic fields) are successfully compared both to experimental data and to numerical simulations provided by a commercial software based on the finite element technique. The results obtained show that the new technique is specially suitable for the efficient full-wave analysis of complex waveguide devices considering an integrated coaxial excitation, where the coaxial probes may be in contact with the metallic insets of the component.
Real-time 3D millimeter wave imaging based FMCW using GGD focal plane array as detectors
NASA Astrophysics Data System (ADS)
Levanon, Assaf; Rozban, Daniel; Kopeika, Natan S.; Yitzhaky, Yitzhak; Abramovich, Amir
2014-03-01
Millimeter wave (MMW) imaging systems are required for applications in medicine, communications, homeland security, and space technology. This is because there is no known ionization hazard for biological tissue, and atmospheric attenuation in this range of the spectrum is relatively low. The lack of inexpensive room temperature imaging systems makes it difficult to give a suitable MMW system for many of the above applications. 3D MMW imaging system based on chirp radar was studied previously using a scanning imaging system of a single detector. The system presented here proposes to employ a chirp radar method with a Glow Discharge Detector (GDD) Focal Plane Array (FPA) of plasma based detectors. Each point on the object corresponds to a point in the image and includes the distance information. This will enable 3D MMW imaging. The radar system requires that the millimeter wave detector (GDD) will be able to operate as a heterodyne detector. Since the source of radiation is a frequency modulated continuous wave (FMCW), the detected signal as a result of heterodyne detection gives the object's depth information according to value of difference frequency, in addition to the reflectance of the image. In this work we experimentally demonstrate the feasibility of implementing an imaging system based on radar principles and FPA of GDD devices. This imaging system is shown to be capable of imaging objects from distances of at least 10 meters.
NASA Astrophysics Data System (ADS)
Xie, J.; Yang, Y.; Ni, S.; Zhao, K.
2015-12-01
In the past decade, ambient noise tomography (ANT) has become an estimated method to construct the earth's interior structures thanks to its advantage in extracting surface waves from cross-correlations of ambient noise without using earthquake data. However, most of previous ambient noise tomography studies concentrate on short and intermediate periods (<50sec) due to the dominant energy of the microseism at these periods. Studies of long period surface waves from cross-correlation of ambient noise are limited. In this study, we verify the accuracy of the long period (50-250sec) surface wave (Rayleigh wave) from ambient noise by comparing both dispersion curves and seismic waveforms from ambient noise with those from earthquake records quantitatively. After that, we calculate vertical-vertical cross-correlation functions among more than1800 USArray Transportable Array stations and extract high quality interstation phase velocity dispersion curves from them at 10-200 sec periods. Then, we adopt a finite frequency ambient noise tomography method based on Born approximation to obtain high resolution phase velocity maps using the obtained dispersion measurements at 10-150 sec periods. Afterward, we extract local dispersion curves from these dispersion maps and invert them for 1D shear wave velocity profiles at individual grids using a Bayesian Monte Carlo method. Finally, a 3D shear velocity model is constructed by assembling all the 1D Vs profiles. Our 3D model is overall similar to other models constructed using earthquake surface waves and body waves. In summary, we demonstrate that the long period surface waves can be extracted from ambient noise, and the long period dispersion measurements from ambient noise are as accurate as those from earthquake data and can be used to construct 3D lithospheric structure from surface down to lithosphere/asthenosphere depths.
Chen, Jiankang; Wang, Wencai; Wang, Ji; Yang, Zengtao; Yang, Jiashi
2008-08-01
We studied thickness vibration of 2 elastic layers with an elastic interface mounted on a plate piezoelectric resonator. The effect of the interface elasticity on resonant frequencies was examined. The result obtained suggests an acoustic wave sensor for measuring the elastic property of an interface between 2 materials. PMID:18986911
Finite Difference Elastic Wave Field Simulation On GPU
NASA Astrophysics Data System (ADS)
Hu, Y.; Zhang, W.
2011-12-01
Numerical modeling of seismic wave propagation is considered as a basic and important aspect in investigation of the Earth's structure, and earthquake phenomenon. Among various numerical methods, the finite-difference method is considered one of the most efficient tools for the wave field simulation. However, with the increment of computing scale, the power of computing has becoming a bottleneck. With the development of hardware, in recent years, GPU shows powerful computational ability and bright application prospects in scientific computing. Many works using GPU demonstrate that GPU is powerful . Recently, GPU has not be used widely in the simulation of wave field. In this work, we present forward finite difference simulation of acoustic and elastic seismic wave propagation in heterogeneous media on NVIDIA graphics cards with the CUDA programming language. We also implement perfectly matched layers on the graphics cards to efficiently absorb outgoing waves on the fictitious edges of the grid Simulations compared with the results on CPU platform shows reliable accuracy and remarkable efficiency. This work proves that GPU can be an effective platform for wave field simulation, and it can also be used as a practical tool for real-time strong ground motion simulation.
Fatigue crack detection in metallic structures with Lamb waves and 3D laser vibrometry
NASA Astrophysics Data System (ADS)
Staszewski, W. J.; Lee, B. C.; Traynor, R.
2007-03-01
The paper presents the application of ultrasonic guided waves for fatigue crack detection in metallic structures. The study involves a simple fatigue test performed to introduce a crack into an aluminium plate. Lamb waves generated by a low-profile, surface-bonded piezoceramic transducer are sensed using a tri-axis, multi-position scanning laser vibrometer. The results demonstrate the potential of laser vibrometry for simple, rapid and robust detection of fatigue cracks in metallic structures. The method could be used in quality inspection and in-service maintenance of metallic structures in aerospace, civil and mechanical engineering industries.
Non-linear 3D Born Shear Wave Tomography in Southeastern Asia
NASA Astrophysics Data System (ADS)
Cao, A.; Panning, M.; Kim, A.; Romanowicz, B.
2007-12-01
We have developed a 3D radially anisotropic shear velocity model of the upper mantle in southeastern Asia from the inversion of long period seismic multimode waveforms. Our approach is based on normal mode perturbation theory, specifically, on a recent modification of the Born approximation, which we call "N-Born", and which includes a non-linear term that allows the accurate inclusion of accumulated phase shifts which arise when the wavepath traverses a spatially extended region with a smooth velocity anomaly of constant sign. We apply the N-Born approximation in the forward modeling part and calculate linear 3D Born kernels in the inverse part. Our starting model is a 3D radially anisotropic model which we derived from a large dataset of teleseismic multimode long period waveforms in the period range 60 to 400 s, using a finite-frequency 2D approximation (NACT, Li and Romanowicz, 1995). This model covered a larger region of East Asia (longitude 30 to 150 degrees and latitude -10 to 60 degrees), while our N-Born model is restricted to a smaller subregion (longitude 75 to 150 degrees and latitude 0 to 45 degrees) for computational efficiency. In this subregion, our N-Born isotropic and anisotropic models are both parameterized at relatively short wavelengths corresponding to a spherical spline level 6 (~200km). Our N-Born model can fit waveforms as well as the NACT model, with up to ~ 83% variance reduction. While the models agree in general, the N-Born isotropic model shows a stronger fast velocity anomaly beneath the Tibetan plateau in the depth range of 150 km to 250 km, which disappears at greater depth, consistent with other studies. More importantly, the N-Born anisotropic model can recover well the downwelling structure associated with subducted slabs. Beneath the Tibet plateau, radial anisotropy shows VSH>VSV, which is indicative of horizontal rather than vertical flow and may help distinguish between end member models of the tectonics of Tibet.
Luquet, David; Marchiano, Régis; Coulouvrat, François
2015-10-28
Many situations involve the propagation of acoustical shock waves through flows. Natural sources such as lightning, volcano explosions, or meteoroid atmospheric entries, emit loud, low frequency, and impulsive sound that is influenced by atmospheric wind and turbulence. The sonic boom produced by a supersonic aircraft and explosion noises are examples of intense anthropogenic sources in the atmosphere. The Buzz-Saw-Noise produced by turbo-engine fan blades rotating at supersonic speed also propagates in a fast flow within the engine nacelle. Simulating these situations is challenging, given the 3D nature of the problem, the long range propagation distances relative to the central wavelength, the strongly nonlinear behavior of shocks associated to a wide-band spectrum, and finally the key role of the flow motion. With this in view, the so-called FLHOWARD (acronym for FLow and Heterogeneous One-Way Approximation for Resolution of Diffraction) method is presented with three-dimensional applications. A scalar nonlinear wave equation is established in the framework of atmospheric applications, assuming weak heterogeneities and a slow wind. It takes into account diffraction, absorption and relaxation properties of the atmosphere, quadratic nonlinearities including weak shock waves, heterogeneities of the medium in sound speed and density, and presence of a flow (assuming a mean stratified wind and 3D turbulent ? flow fluctuations of smaller amplitude). This equation is solved in the framework of the one-way method. A split-step technique allows the splitting of the non-linear wave equation into simpler equations, each corresponding to a physical effect. Each sub-equation is solved using an analytical method if possible, and finite-differences otherwise. Nonlinear effects are solved in the time domain, and others in the frequency domain. Homogeneous diffraction is handled by means of the angular spectrum method. Ground is assumed perfectly flat and rigid. Due to the 3D
NASA Astrophysics Data System (ADS)
Luquet, David; Marchiano, Régis; Coulouvrat, François
2015-10-01
Many situations involve the propagation of acoustical shock waves through flows. Natural sources such as lightning, volcano explosions, or meteoroid atmospheric entries, emit loud, low frequency, and impulsive sound that is influenced by atmospheric wind and turbulence. The sonic boom produced by a supersonic aircraft and explosion noises are examples of intense anthropogenic sources in the atmosphere. The Buzz-Saw-Noise produced by turbo-engine fan blades rotating at supersonic speed also propagates in a fast flow within the engine nacelle. Simulating these situations is challenging, given the 3D nature of the problem, the long range propagation distances relative to the central wavelength, the strongly nonlinear behavior of shocks associated to a wide-band spectrum, and finally the key role of the flow motion. With this in view, the so-called FLHOWARD (acronym for FLow and Heterogeneous One-Way Approximation for Resolution of Diffraction) method is presented with three-dimensional applications. A scalar nonlinear wave equation is established in the framework of atmospheric applications, assuming weak heterogeneities and a slow wind. It takes into account diffraction, absorption and relaxation properties of the atmosphere, quadratic nonlinearities including weak shock waves, heterogeneities of the medium in sound speed and density, and presence of a flow (assuming a mean stratified wind and 3D turbulent ? flow fluctuations of smaller amplitude). This equation is solved in the framework of the one-way method. A split-step technique allows the splitting of the non-linear wave equation into simpler equations, each corresponding to a physical effect. Each sub-equation is solved using an analytical method if possible, and finite-differences otherwise. Nonlinear effects are solved in the time domain, and others in the frequency domain. Homogeneous diffraction is handled by means of the angular spectrum method. Ground is assumed perfectly flat and rigid. Due to the 3D
Internal wave attractors examined using laboratory experiments and 3D numerical simulations
NASA Astrophysics Data System (ADS)
Brouzet, C.; Sibgatullin, I. N.; Scolan, H.; Ermanyuk, E. V.; Dauxois, T.
2016-04-01
In the present paper, we combine numerical and experimental approaches to study the dynamics of stable and unstable internal wave attractors. The problem is considered in a classic trapezoidal setup filled with a uniformly stratified fluid. Energy is injected into the system at global scale by the small-amplitude motion of a vertical wall. Wave motion in the test tank is measured with the help of conventional synthetic schlieren and PIV techniques. The numerical setup closely reproduces the experimental one in terms of geometry and the operational range of the Reynolds and Schmidt numbers. The spectral element method is used as a numerical tool to simulate the nonlinear dynamics of a viscous salt-stratified fluid. We show that the results of three-dimensional calculations are in excellent qualitative and quantitative agreement with the experimental data, including the spatial and temporal parameters of the secondary waves produced by triadic resonance instability. Further, we explore experimentally and numerically the effect of lateral walls on secondary currents and spanwise distribution of velocity amplitudes in the wave beams. Finally, we test the assumption of a bidimensional flow and estimate the error made in synthetic schlieren measurements due to this assumption.
A 3D MPI-Parallel GPU-accelerated framework for simulating ocean wave energy converters
NASA Astrophysics Data System (ADS)
Pathak, Ashish; Raessi, Mehdi
2015-11-01
We present an MPI-parallel GPU-accelerated computational framework for studying the interaction between ocean waves and wave energy converters (WECs). The computational framework captures the viscous effects, nonlinear fluid-structure interaction (FSI), and breaking of waves around the structure, which cannot be captured in many potential flow solvers commonly used for WEC simulations. The full Navier-Stokes equations are solved using the two-step projection method, which is accelerated by porting the pressure Poisson equation to GPUs. The FSI is captured using the numerically stable fictitious domain method. A novel three-phase interface reconstruction algorithm is used to resolve three phases in a VOF-PLIC context. A consistent mass and momentum transport approach enables simulations at high density ratios. The accuracy of the overall framework is demonstrated via an array of test cases. Numerical simulations of the interaction between ocean waves and WECs are presented. Funding from the National Science Foundation CBET-1236462 grant is gratefully acknowledged.
Summary of work on shock wave feature extraction in 3-D datasets
NASA Technical Reports Server (NTRS)
Hesselink, Lambertus (Principal Investigator)
1996-01-01
A method for extracting and visualizing shock waves from three dimensional data-sets is discussed. Issues concerning computation time, robustness to numerical perturbations, and noise introduction are considered and compared with other methods. Finally, results using this method are discussed.
NASA Astrophysics Data System (ADS)
Moulik, P.; Ekstrom, G.
2012-12-01
We have developed a framework that can be used to investigate anisotropic velocity, density and anelastic heterogeneity in the Earth's mantle using a wide spectrum (0.3-50 mHz) of seismological observables. We start with the extensive dataset of surface-wave phase anomalies, long-period waveforms, and body-wave travel times collected by Kustowski et al. (2008) for the development of the global model S362ANI. The additional data included in our analysis are splitting functions of spheroidal and toroidal modes, which are analogous to phase velocity maps at low frequencies. We include in this set of observations a new dataset containing the splitting functions of 56 spheroidal fundamental modes and overtones, measured by Deuss et al. (2011, 2012) using data from large recent earthquakes. Apart from providing unique constraints on the long-wavelength elastic and density structure in the mantle, the overtone splitting data are especially sensitive to the velocity (and anisotropic) structure in the transition zone and in the deeper mantle. The detection of anisotropy, a marker of flow, in the transition zone has implications for our understanding of mantle convection. Our forward modeling of the splitting functions, like the other types of data, includes the effects of radial anisotropy (Mochizuki, 1986). We show that the upper-mantle shear-wave anisotropy of S362ANI generates a clear contribution to the splitting functions of the modes that are sensitive to the upper-mantle structure. We explore the tradeoffs between fitting the mode splitting functions and the travel-times of body waves that turn in the transition zone or in the lower mantle (e.g. SS), while observing that the waveforms and the surface wave phase-anomalies provide complementary information about the mantle. Our experiments suggest that the splitting data are sufficiently sensitive to the anisotropy in the mantle such that their inclusion may provide a better depth resolution of the anisotropic shear
Elastic guided waves in a coated spherical shell
NASA Astrophysics Data System (ADS)
Qiao, Song; Shang, Xinchun; Pan, Ernian
2016-04-01
Elastic-guided wave inspection technique is important in non-destructive detection of coated shell structures. It is based on the wave propagation characteristics and various factors which influences it. In this paper, the dispersion equations of the spherical shell are derived by the decomposition approach in order to investigate the influences of the coating thickness and viscoelastic damping on the dispersion characteristics. The viscoelastic properties of the coating layer are modelled by the standard linear solid with two damping factors in the Láme constants. The dispersion equation of the coated shell is deduced by the transfer matrix method, and the dispersion and attenuation curves for different thicknesses and damping factors are calculated. The frequency range which is less affected by coating is identified by comparing the dispersion curves of the bare shell to those of the coated shell with different coating thicknesses. The effect of damping factors on the mode shapes is also examined. The present numerical results on the elastic guided wave in coated spherical shell would provide a theoretical basis for non-destructive inspections in layered spherical shell structures.
Miles, A R; Edwards, M J; Greenough, J A
2004-11-08
Perturbations on an interface driven by a strong blast wave grow in time due to a combination of Rayleigh-Taylor, Richtmyer-Meshkov, and decompression effects. In this paper, results from three-dimensional numerical simulations of such a system under drive conditions to be attainable on the National Ignition Facility [E. M. Campbell, Laser Part. Beams, 9(2), 209 (1991)] are presented. Using the multi-physics, adaptive mesh refinement, higher order Godunov Eulerian hydrocode, Raptor [L. H. Howell and J.A. Greenough, J. Comp. Phys. 184, 53 (2003)], the late nonlinear instability evolution, including transition to turbulence, is considered for various multimode perturbation spectra. The 3D post-transition state differs from the 2D result, but the process of transition proceeds similarly in both 2D and 3D. The turbulent mixing transition results in a reduction in the growth rate of the mixing layer relative to its pre-transition value and, in the case of the bubble front, relative to the 2D result. The post-transition spike front velocity is approximately the same in 2D and 3D. Implications for hydrodynamic mixing in core-collapse supernova are discussed.