NASA Technical Reports Server (NTRS)
Dahl, Milo D.
2000-01-01
An acoustic source inside of a 2-D jet excites an instability wave in the shear layer resulting in sound radiating away from the shear layer. Solve the linearized Euler equations to predict the sound radiation outside of the jet. The jet static pressure is assumed to be constant. The jet flow is parallel and symmetric about the x-axis. Use a symmetry boundary condition along the x-axis.
A scanning-mode 2D shear wave imaging (s2D-SWI) system for ultrasound elastography.
Qiu, Weibao; Wang, Congzhi; Li, Yongchuan; Zhou, Juan; Yang, Ge; Xiao, Yang; Feng, Ge; Jin, Qiaofeng; Mu, Peitian; Qian, Ming; Zheng, Hairong
2015-09-01
Ultrasound elastography is widely used for the non-invasive measurement of tissue elasticity properties. Shear wave imaging (SWI) is a quantitative method for assessing tissue stiffness. SWI has been demonstrated to be less operator dependent than quasi-static elastography, and has the ability to acquire quantitative elasticity information in contrast with acoustic radiation force impulse (ARFI) imaging. However, traditional SWI implementations cannot acquire two dimensional (2D) quantitative images of the tissue elasticity distribution. This study proposes and evaluates a scanning-mode 2D SWI (s2D-SWI) system. The hardware and image processing algorithms are presented in detail. Programmable devices are used to support flexible control of the system and the image processing algorithms. An analytic signal based cross-correlation method and a Radon transformation based shear wave speed determination method are proposed, which can be implemented using parallel computation. Imaging of tissue mimicking phantoms, and in vitro, and in vivo imaging test are conducted to demonstrate the performance of the proposed system. The s2D-SWI system represents a new choice for the quantitative mapping of tissue elasticity, and has great potential for implementation in commercial ultrasound scanners. PMID:26025508
2D instabilities of surface gravity waves on a linear shear current
NASA Astrophysics Data System (ADS)
Francius, Marc; Kharif, Christian
2016-04-01
Periodic 2D surface water waves propagating steadily on a rotational current have been studied by many authors (see [1] and references therein). Although the recent important theoretical developments have confirmed that periodic waves can exist over flows with arbitrary vorticity, their stability and their nonlinear evolution have not been much studied extensively so far. In fact, even in the rather simple case of uniform vorticity (linear shear), few papers have been published on the effect of a vertical shear current on the side-band instability of a uniform wave train over finite depth. In most of these studies [2-5], asymptotic expansions and multiple scales method have been used to obtain envelope evolution equations, which allow eventually to formulate a condition of (linear) instability to long modulational perturbations. It is noted here that this instability is often referred in the literature as the Benjamin-Feir or modulational instability. In the present study, we consider the linear stability of finite amplitude two-dimensional, periodic water waves propagating steadily on the free surface of a fluid with constant vorticity and finite depth. First, the steadily propagating surface waves are computed with steepness up to very close to the highest, using a Fourier series expansions and a collocation method, which constitutes a simple extension of Fenton's method [6] to the cases with a linear shear current. Then, the linear stability of these permanent waves to infinitesimal 2D perturbations is developed from the fully nonlinear equations in the framework of normal modes analysis. This linear stability analysis is an extension of [7] to the case of waves in the presence of a linear shear current and permits the determination of the dominant instability as a function of depth and vorticity for a given steepness. The numerical results are used to assess the accuracy of the vor-NLS equation derived in [5] for the characteristics of modulational
Generation and Radiation of Acoustic Waves from a 2D Shear Layer
NASA Technical Reports Server (NTRS)
Dahl, Milo D.
2000-01-01
A thin free shear layer containing an inflection point in the mean velocity profile is inherently unstable. Disturbances in the flow field can excite the unstable behavior of a shear layer, if the appropriate combination of frequencies and shear layer thicknesses exists, causing instability waves to grow. For other combinations of frequencies and thicknesses, these instability waves remain neutral in amplitude or decay in the downstream direction. A growing instability wave radiates noise when its phase velocity becomes supersonic relative to the ambient speed of sound. This occurs primarily when the mean jet flow velocity is supersonic. Thus, the small disturbances in the flow, which themselves may generate noise, have generated an additional noise source. It is the purpose of this problem to test the ability of CAA to compute this additional source of noise. The problem is idealized such that the exciting disturbance is a fixed known acoustic source pulsating at a single frequency. The source is placed inside of a 2D jet with parallel flow; hence, the shear layer thickness is constant. With the source amplitude small enough, the problem is governed by the following set of linear equations given in dimensional form.
Generation and Radiation of Acoustic Waves from a 2-D Shear Layer using the CE/SE Method
NASA Technical Reports Server (NTRS)
Loh, Ching Y.; Wang, Xiao Y.; Chang, Sin-Chung; Jorgenson, Philip C. E.
2000-01-01
In the present work, the generation and radiation of acoustic waves from a 2-D shear layer problem is considered. An acoustic source inside of a 2-D jet excites an instability wave in the shear layer, resulting in sound Mach radiation. The numerical solution is obtained by solving the Euler equations using the space time conservation element and solution element (CE/SE) method. Linearization is achieved through choosing a small acoustic source amplitude. The Euler equations are nondimensionalized as instructed in the problem statement. All other conditions are the same except that the Crocco's relation has a slightly different form. In the following, after a brief sketch of the CE/SE method, the numerical results for this problem are presented.
Periodically sheared 2D Yukawa systems
Kovács, Anikó Zsuzsa; Hartmann, Peter; Donkó, Zoltán
2015-10-15
We present non-equilibrium molecular dynamics simulation studies on the dynamic (complex) shear viscosity of a 2D Yukawa system. We have identified a non-monotonic frequency dependence of the viscosity at high frequencies and shear rates, an energy absorption maximum (local resonance) at the Einstein frequency of the system at medium shear rates, an enhanced collective wave activity, when the excitation is near the plateau frequency of the longitudinal wave dispersion, and the emergence of significant configurational anisotropy at small frequencies and high shear rates.
Generation and Radiation of Acoustic Waves from a 2-D Shear Layer
NASA Technical Reports Server (NTRS)
Agarwal, Anurag; Morris, Philip J.
2000-01-01
A parallel numerical simulation of the radiation of sound from an acoustic source inside a 2-D jet is presented in this paper. This basic benchmark problem is used as a test case for scattering problems that are presently being solved by using the Impedance Mismatch Method (IMM). In this technique, a solid body in the domain is represented by setting the acoustic impedance of each medium, encountered by a wave, to a different value. This impedance discrepancy results in reflected and scattered waves with appropriate amplitudes. The great advantage of the use of this method is that no modifications to a simple Cartesian grid need to be made for complicated geometry bodies. Thus, high order finite difference schemes may be applied simply to all parts of the domain. In the IMM, the total perturbation field is split into incident and scattered fields. The incident pressure is assumed to be known and the equivalent sources for the scattered field are associated with the presence of the scattering body (through the impedance mismatch) and the propagation of the incident field through a non-uniform flow. An earlier version of the technique could only handle uniform flow in the vicinity of the source and at the outflow boundary. Scattering problems in non-uniform mean flow are of great practical importance (for example, scattering from a high lift device in a non-uniform mean flow or the effects of a fuselage boundary layer). The solution to this benchmark problem, which has an acoustic wave propagating through a non-uniform mean flow, serves as a test case for the extensions of the IMM technique.
Gersak, Mariana M; Badea, Radu; Lenghel, Lavinia M; Vasilescu, Dan; Botar-Jid, Carolina; Dudea, Sorin M
2016-06-01
Transient elastography and Acoustic Radiation Force Impulse imaging are useful non-invasive methods for liver stiffness estimation, although both are influenced by food intake. The aim of the work described here was to identify liver stiffness variation after a standardized meal using 2-D shear wave elastography. Liver stiffness was estimated in 31 apparently healthy subjects, under fasting conditions and after a standardized meal (20, 40, 60, 80, 100 and 120 min after food intake). In most of the cases, liver stiffness values increased between 20 and 40 min after the meal (p < 0.05) and then significantly decreased between 60 and 80 min (p < 0.05). At 120 min after food intake, liver stiffness values were significantly lower compared with liver stiffness values under fasting conditions (p < 0.05). Gender, but not body mass index, had an important role in liver stiffness variation after food intake (p < 0.01). In conclusion, to avoid the influence of food intake on liver stiffness estimation, 2-D shear wave elastography should be performed only under fasting conditions. PMID:26947447
Li, Changtian; Zhang, Changsheng; Li, Junlai; Cao, Xiaolin; Song, Danfei
2016-07-01
2-D Shear wave elastography (SWE) imaging is widely used in clinical practice, and some researchers have applied this technique in the evaluation of neonatal brains. However, the immediate and long-term impacts of dynamic radiation force exposure on the neonatal central nervous system remain unknown. In this study, we exposed neonatal mice to 2-D SWE scanning for 10 min, 20 min and 30 min under diagnostic mode (mechanical index [MI]: 1.3; thermal index [TI]: 0.5), respectively. For the control group, the neonatal mice were sham irradiated for 30 min with the machine powered off. Their brains were collected and analyzed using histologic staining and western blot analysis at 24 h and 3 mo after the 2-D SWE scanning. The Morris water maze (MWM) test was used to assess learning and memory function of the mice at 3 mo of age. The results indicated that using 2-D SWE in evaluating brains of neonatal mice does not cause detectable histologic changes, nor does it have long-term effects on their learning and memory abilities. However, the PI3 K/AKT/mTOR pathway was disturbed when the 2-D SWE scanning lasted for more than 30 min, and the expression of p-PKCa was suppressed by 10 min or more in 2-D SWE scanning. Although these injuries may be self-repaired as the mice grow, more attention should be paid to the scanning duration when applying 2-D-SWE elastography in the assessment of neonatal brains. PMID:27112914
Tian, Wen-Shuo; Lin, Man-Xia; Zhou, Lu-Yao; Pan, Fu-Shun; Huang, Guang-Liang; Wang, Wei; Lu, Ming-De; Xie, Xiao-Yan
2016-09-01
The goal of the work described here was to evaluate the diagnostic efficacy of 2-D shear wave elastography (2-D SWE) in differentiating malignancy from benign focal liver lesions (FLLs). The maxima, minima, means and the standard deviations of 2-D SWE measurements, expressed in kilopascals (Emax, Emin, Emean, ESD), were obtained for 221 patients with 229 FLLs. Receiver operating characteristic curve analysis was performed to evaluate the diagnostic performance of 2-D SWE. The Mann-Whitney U-test was used to assess inter-group differences. Emax, Emin, Emean and ESD were significantly higher in the 164 malignant lesions than in the 65 benign lesions (p < 0.001). For identification of malignant FLLs, the areas under receiver operating characteristic curves for Emax, Emin, Emean and ESD were 0.920, 0.710, 0.879 and 0.915, respectively. Emax was 96.21 ± 35.40 for 19 intrahepatic cholangiocarcinomas and 90.32 ± 54.71 for 35 liver metastatic lesions, which were significantly higher than 61.83 ± 28.87 for 103 hepatocellular carcinomas (p < 0.0001 and p = 0.0237). Emax was 38.72 ± 18.65 for 15 focal nodular hyperplasias, which was significantly higher than 20.56 ± 10.74 for 37 hemangiomas (p = 0.0009). The Emax values for adjacent liver parenchyma of hepatocellular carcinomas and intrahepatic cholangiocarcinomas were significantly higher than those for the other three lesion types (p < 0.005). In conclusion, Emax values of FLLs and adjacent liver parenchyma could help in differentiating malignant from benign FLLs. PMID:27283039
Arnal, B; Pinton, G; Garapon, P; Pernot, M; Fink, M; Tanter, M
2013-10-01
Shear wave imaging (SWI) maps soft tissue elasticity by measuring shear wave propagation with ultrafast ultrasound acquisitions (10 000 frames s(-1)). This spatiotemporal data can be used as an input for an inverse problem that determines a shear modulus map. Common inversion methods are local: the shear modulus at each point is calculated based on the values of its neighbour (e.g. time-of-flight, wave equation inversion). However, these approaches are sensitive to the information loss such as noise or the lack of the backscattered signal. In this paper, we evaluate the benefits of a global approach for elasticity inversion using a least-squares formulation, which is derived from full waveform inversion in geophysics known as the adjoint method. We simulate an acoustic waveform in a medium with a soft and a hard lesion. For this initial application, full elastic propagation and viscosity are ignored. We demonstrate that the reconstruction of the shear modulus map is robust with a non-uniform background or in the presence of noise with regularization. Compared to regular local inversions, the global approach leads to an increase of contrast (∼+3 dB) and a decrease of the quantification error (∼+2%). We demonstrate that the inversion is reliable in the case when there is no signal measured within the inclusions like hypoechoic lesions which could have an impact on medical diagnosis. PMID:24018867
Wang, Michael; Byram, Brett; Palmeri, Mark; Rouze, Ned; Nightingale, Kathryn
2013-01-01
A 2D matrix ultrasound array is used to monitor acoustic radiation force impulse (ARFI) induced shear wave propagation in 3D in excised canine muscle. From a single acquisition, both the shear wave phase and group velocity can be calculated to estimate the shear wave speed (SWS) along and across the fibers, as well as the fiber orientation in 3D. The true fiber orientation found using the 3D Radon Transform on B-mode volumes of the muscle was used to verify the fiber direction estimated from shear wave data. For the simplified imaging case when the ARFI push can be oriented perpendicular to the fibers, the error in estimating the fiber orientation using phase and group velocity measurements was 3.5 ±2.6° and 3.4 ±1.4° (mean ± standard deviation), respectively, over six acquisitions in different muscle samples. For the more general case when the push is oblique to the fibers, the angle between the push and the fibers is found using the dominant orientation of the shear wave displacement magnitude. In 30 acquisitions on six different muscle samples with oblique push angles up to 40°, the error in the estimated fiber orientation using phase and group velocity measurements was 5.4±2.9° and 5.3±3.2°, respectively, after estimating and accounting for the additional unknown push angle. Either the phase or group velocity measurements can be used to estimate fiber orientation and SWS along and across the fibers. Although it is possible to perform these measurements when the push is not perpendicular to the fibers, highly oblique push angles induce lower shear wave amplitudes which can cause inaccurate SWS measurements. PMID:23686942
NASA Astrophysics Data System (ADS)
Louie, J. N.; Pancha, A.; Pullammanappallil, S. K.
2014-12-01
Refraction microtermor routinely assesses 1D and 2D velocity-depth profiles to shallow depths of approximately 100 m primarily for engineering applications. Estimation of both shallow and deep (>100 m) shear-velocity structure are key elements in the assessment of urban areas for potential earthquake ground shaking, damage, and the calibration of recorded ground motions. Three independent studies investigated the ability of the refraction microtremor technology to image deep velocity structure, to depths exceeding 1 km (Deep ReMi). In the first study, we were able to delineate basin thicknesses of up to 900 m and the deep-basin velocity structure beneath the Reno-area basin. Constraints on lateral velocity changes in 3D as well as on velocity profiles extended down to 1500 m, and show a possible fault offset. This deployment used 30 stand-alone wireless instruments mated to 4.5 Hz geophones, along each of five arrays 2.9 to 5.8 km long. Rayleigh-wave dispersion was clear at frequencies as low as 0.5 Hz using 120 sec ambient urban noise records. The results allowed construction of a 3D velocity model, vetted by agreement with gravity studies. In a second test, a 5.8 km array delimited the southern edge of the Tahoe Basin, with constraints from gravity. There, bedrock depth increased by 250 m in thickness over a distance of 1600 m, with definition of the velocity of the deeper basin sediments. The third study delineated the collapse region of an underground nuclear explosion within a thick sequence of volcanic extrusives, using a shear-wave minivibe in a radial direction, and horizontal geophones. Analysis showed the cavity extends to 620 m depth, with a width of 180 m and a height of 220 m. Our results demonstrate that deep velocity structure can be recovered using ambient noise. This technique offers the ability to define 2D and 3D structural representations essential for seismic hazard evaluation.
Song, Pengfei; Mellema, Daniel C.; Sheedy, Shannon P.; Meixner, Duane D.; Karshen, Ryan M.; Urban, Matthew W.; Manduca, Armando; Sanchez, William; Callstrom, Matthew R.; Greenleaf, James F.; Chen, Shigao
2015-01-01
Objective To investigate the correlation between 2-D ultrasound shear wave elastography (SWE) and magnetic resonance elastography (MRE) in liver stiffness measurement and the diagnostic performance of 2-D SWE for liver fibrosis when imaging from different intercostal spaces and using MRE as the reference standard. Methods 2-D SWE was performed on 47 patients (22 females and 25 males, age 19–77) using the GE LOGIQ E9 scanner. Each of the 47 patients had same day MRE obtained for clinical purposes. The study was HIPAA-compliant and approved by the institutional review board. Informed consent was obtained from each subject. 2-D SWE measurements were acquired from the 9th, 8th, and 7th intercostal spaces. Correlation with MRE was calculated at each intercostal space and multiple intercostal spaces combined. The performance of 2-D SWE in diagnosing liver fibrosis was evaluated with receiver operating characteristic (ROC) curve analysis using MRE as the standard. Results The highest correlation between 2-D SWE and MRE was from the 8th and 7th intercostal spaces (r = 0.68 – 0.76). The range of the areas under the ROC curve for separating normal or inflamed livers from fibrotic livers using MRE as the clinical reference were 0.84 – 0.92 when using 8th and 7th intercostal spaces individually, and 0.89 –0.9 when combined. Conclusion The results suggest that 2-D SWE and MRE are well correlated when SWE is performed at the 8th and 7th intercostal spaces. The 9th intercostal space is less reliable for diagnosing fibrosis using 2-D SWE. Combining measurements from multiple intercostal spaces does not significantly improve 2-D SWE performance for the detection of fibrosis. PMID:26782164
Luo, Y.; Xia, J.; Liu, J.; Xu, Y.; Liu, Q.
2008-01-01
Multichannel Analysis of Surface Waves utilizes a multichannel recording system to estimate near-surface shear (S)-wave velocities from high-frequency Rayleigh waves. A pseudo-2D S-wave velocity (vS) section is constructed by aligning 1D models at the midpoint of each receiver spread and using a spatial interpolation scheme. The horizontal resolution of the section is therefore most influenced by the receiver spread length and the source interval. The receiver spread length sets the theoretical lower limit and any vS structure with its lateral dimension smaller than this length will not be properly resolved in the final vS section. A source interval smaller than the spread length will not improve the horizontal resolution because spatial smearing has already been introduced by the receiver spread. In this paper, we first analyze the horizontal resolution of a pair of synthetic traces. Resolution analysis shows that (1) a pair of traces with a smaller receiver spacing achieves higher horizontal resolution of inverted S-wave velocities but results in a larger relative error; (2) the relative error of the phase velocity at a high frequency is smaller than at a low frequency; and (3) a relative error of the inverted S-wave velocity is affected by the signal-to-noise ratio of data. These results provide us with a guideline to balance the trade-off between receiver spacing (horizontal resolution) and accuracy of the inverted S-wave velocity. We then present a scheme to generate a pseudo-2D S-wave velocity section with high horizontal resolution using multichannel records by inverting high-frequency surface-wave dispersion curves calculated through cross-correlation combined with a phase-shift scanning method. This method chooses only a pair of consecutive traces within a shot gather to calculate a dispersion curve. We finally invert surface-wave dispersion curves of synthetic and real-world data. Inversion results of both synthetic and real-world data demonstrate that
Energy Science and Technology Software Center (ESTSC)
2004-08-01
AnisWave2D is a 2D finite-difference code for a simulating seismic wave propagation in fully anisotropic materials. The code is implemented to run in parallel over multiple processors and is fully portable. A mesh refinement algorithm has been utilized to allow the grid-spacing to be tailored to the velocity model, avoiding the over-sampling of high-velocity materials that usually occurs in fixed-grid schemes.
Diagnostic Accuracy of 2D-Shear Wave Elastography for Liver Fibrosis Severity: A Meta-Analysis
Jiang, Tian’an; Tian, Guo; Zhao, Qiyu; Kong, Dexing; Cheng, Chao; Zhong, Liyun; Li, Lanjuan
2016-01-01
Purpose To evaluate the accuracy of shear wave elastography (SWE) in the quantitative diagnosis of liver fibrosis severity. Methods The published literatures were systematically retrieved from PubMed, Embase, Web of science and Scopus up to May 13th, 2016. Included studies reported the pooled sensitivity, specificity, positive and negative predictive values, as well as the diagnostic odds ratio of SWE in populations with liver fibrosis. A bivariate mixed-effects regression model was used, which was estimated by the I2 statistics. The quality of articles was evaluated by quality assessment of diagnostic accuracy studies (QUADAS). Results Thirteen articles including 2303 patients were qualified for the study. The pooled sensitivity and specificity of SWE for the diagnosis of liver fibrosis are as follows: ≥F1 0.76 (p<0.001, 95% CI, 0.71–0.81, I2 = 75.33%), 0.92 (p<0.001, 95% CI, 0.80–0.97, I2 = 79.36%); ≥F2 0.84 (p = 0.35, 95% CI, 0.81–0.86, I2 = 9.55%), 0.83 (p<0.001, 95% CI, 0.77–0.88, I2 = 86.56%); ≥F3 0.89 (p = 0.56, 95% CI, 0.86–0.92, I2 = 0%), 0.86 (p<0.001, 95% CI, 0.82–0.90, I2 = 75.73%); F4 0.89 (p = 0.24, 95% CI, 0.84–0.92, I2 = 20.56%), 0.88 (p<0.001, 95% CI, 0.84–0.92, I2 = 82.75%), respectively. Sensitivity analysis showed no significant changes if any one of the studies was excluded. Publication bias was not detected in this meta-analysis. Conclusions Our study suggests that SWE is a helpful method to appraise liver fibrosis severity. Future studies that validate these findings would be appropriate. PMID:27300569
Gerber, Ludmila; Kasper, Daniela; Fitting, Daniel; Knop, Viola; Vermehren, Annika; Sprinzl, Kathrin; Hansmann, Martin L; Herrmann, Eva; Bojunga, Joerg; Albert, Joerg; Sarrazin, Christoph; Zeuzem, Stefan; Friedrich-Rust, Mireen
2015-09-01
Two-dimensional shear wave elastography (2-D SWE) is an ultrasound-based elastography method integrated into a conventional ultrasound machine. It can evaluate larger regions of interest and, therefore, might be better at determining the overall fibrosis distribution. The aim of this prospective study was to compare 2-D SWE with the two best evaluated liver elastography methods, transient elastography and acoustic radiation force impulse (point SWE using acoustic radiation force impulse) imaging, in the same population group. The study included 132 patients with chronic hepatopathies, in which liver stiffness was evaluated using transient elastography, acoustic radiation force impulse imaging and 2-D SWE. The reference methods were liver biopsy for the assessment of liver fibrosis (n = 101) and magnetic resonance imaging/computed tomography for the diagnosis of liver cirrhosis (n = 31). No significant difference in diagnostic accuracy, assessed as the area under the receiver operating characteristic curve (AUROC), was found between the three elastography methods (2-D SWE, transient elastography, acoustic radiation force impulse imaging) for the diagnosis of significant and advanced fibrosis and liver cirrhosis in the "per protocol" (AUROCs for fibrosis stages ≥2: 0.90, 0.95 and 0.91; for fibrosis stage [F] ≥3: 0.93, 0.95 and 0.94; for F = 4: 0.92, 0.96 and 0.92) and "intention to diagnose" cohort (AUROCs for F ≥2: 0.87, 0.92 and 0.91; for F ≥3: 0.91, 0.93 and 0.94; for F = 4: 0.88, 0.90 and 0.89). Therefore, 2-D SWE, ARFI imaging and transient elastography seem to be comparably good methods for non-invasive assessment of liver fibrosis. PMID:26116161
Ultrasonic shear wave couplant
Kupperman, D.S.; Lanham, R.N.
1984-04-11
Ultrasonically testing of an article at high temperatures is accomplished by the use of a compact layer of a dry ceramic powder as a couplant in a method which involves providing an ultrasonic transducer as a probe capable of transmitting shear waves, coupling the probe to the article through a thin compact layer of a dry ceramic powder, propagating a shear wave from the probe through the ceramic powder and into the article to develop echo signals, and analyzing the echo signals to determine at least one physical characteristic of the article.
Ultrasonic shear wave couplant
Kupperman, David S.; Lanham, Ronald N.
1985-01-01
Ultrasonically testing of an article at high temperatures is accomplished by the use of a compact layer of a dry ceramic powder as a couplant in a method which involves providing an ultrasonic transducer as a probe capable of transmitting shear waves, coupling the probe to the article through a thin compact layer of a dry ceramic powder, propagating a shear wave from the probe through the ceramic powder and into the article to develop echo signals, and analyzing the echo signals to determine at least one physical characteristic of the article.
Shear wave transmissivity measurement by color Doppler shear wave imaging
NASA Astrophysics Data System (ADS)
Yamakoshi, Yoshiki; Yamazaki, Mayuko; Kasahara, Toshihiro; Sunaguchi, Naoki; Yuminaka, Yasushi
2016-07-01
Shear wave elastography is a useful method for evaluating tissue stiffness. We have proposed a novel shear wave imaging method (color Doppler shear wave imaging: CD SWI), which utilizes a signal processing unit in ultrasound color flow imaging in order to detect the shear wave wavefront in real time. Shear wave velocity is adopted to characterize tissue stiffness; however, it is difficult to measure tissue stiffness with high spatial resolution because of the artifact produced by shear wave diffraction. Spatial average processing in the image reconstruction method also degrades the spatial resolution. In this paper, we propose a novel measurement method for the shear wave transmissivity of a tissue boundary. Shear wave wavefront maps are acquired by changing the displacement amplitude of the shear wave and the transmissivity of the shear wave, which gives the difference in shear wave velocity between two mediums separated by the boundary, is measured from the ratio of two threshold voltages required to form the shear wave wavefronts in the two mediums. From this method, a high-resolution shear wave amplitude imaging method that reconstructs a tissue boundary is proposed.
A new method for shear wave speed estimation in shear wave elastography.
Engel, Aaron J; Bashford, Gregory R
2015-12-01
Visualization of mechanical properties of tissue can aid in noninvasive pathology diagnosis. Shear wave elastography (SWE) measures the elastic properties of soft tissues by estimation of local shear wave propagation speed. In this paper, a new robust method for estimation of shear wave speed is introduced which has the potential for simplifying continuous filtering and real-time elasticity processing. Shear waves were generated by external mechanical excitation and imaged at a high frame rate. Three homogeneous phantoms of varying elastic moduli and one inclusion phantom were imaged. Waves propagating in separate directions were filtered and shear wave speed was estimated by inversion of the 1-D first-order wave equation. Final 2-D shear wave speed maps were constructed by weighted averaging of estimates from opposite traveling directions. Shear wave speed results for phantoms with gelatin concentrations of 5%, 7%, and 9% were 1.52 ± 0.10 m/s, 1.86 ± 0.10 m/s, and 2.37 ± 0.15 m/s, respectively, which were consistent with estimates computed from three other conventional methods, as well as compression tests done with a commercial texture analyzer. The method was shown to be able to reconstruct a 2-D speed map of an inclusion phantom with good image quality and variance comparable to conventional methods. Suggestions for further work are given. PMID:26670851
NASA Astrophysics Data System (ADS)
Di Fiore, V.; Cavuoto, G.; Tarallo, D.; Punzo, M.; Evangelista, L.
2016-05-01
A joint analysis of down-hole (DH) and multichannel analysis of surface waves (MASW) measurements offers a complete evaluation of shear wave velocity profiles, especially for sites where a strong lateral variability is expected, such as archeological sites. In this complex stratigraphic setting, the high "subsoil anisotropy" (i.e., sharp lithological changes due to the presence of anthropogenic backfill deposits and/or buried man-made structures) implies a different role for DH and MASW tests. This paper discusses some results of a broad experimental program conducted on the Palatine Hill, one of the most ancient areas of the city of Rome (Italy). The experiments were part of a project on seismic microzoning and consisted of 20 MASW and 11 DH tests. The main objective of this study was to examine the difficulties related to the interpretation of the DH and MASW tests and the reliability limits inherent in the application of the noninvasive method in complex stratigraphic settings. As is well known, DH tests provide good determinations of shear wave velocities (Vs) for different lithologies and man-made materials, whereas MASW tests provide average values for the subsoil volume investigated. The data obtained from each method with blind tests were compared and were correlated to site-specific subsurface conditions, including lateral variability. Differences between punctual (DH) and global (MASW) Vs measurements are discussed, quantifying the errors by synthetic comparison and by site response analyses. This study demonstrates that, for archeological sites, VS profiles obtained from the DH and MASW methods differ by more than 15 %. However, the local site effect showed comparable results in terms of natural frequencies, whereas the resolution of the inverted shear wave velocity was influenced by the fundamental mode of propagation.
Calculating tissue shear modulus and pressure by 2D Log-Elastographic methods
McLaughlin, Joyce R; Zhang, Ning; Manduca, Armando
2010-01-01
Shear modulus imaging, often called elastography, enables detection and characterization of tissue abnormalities. In this paper the data is two displacement components obtained from successive MR or ultrasound data sets acquired while the tissue is excited mechanically. A 2D plane strain elastic model is assumed to govern the 2D displacement, u. The shear modulus, μ, is unknown and whether or not the first Lamé parameter, λ, is known the pressure p = λ∇ · u which is present in the plane strain model cannot be measured and is unreliably computed from measured data and can be shown to be an order one quantity in the units kPa. So here we present a 2D Log-Elastographic inverse algorithm that: (1) simultaneously reconstructs the shear modulus, μ, and p, which together satisfy a first order partial differential equation system, with the goal of imaging μ; (2) controls potential exponential growth in the numerical error; and (3) reliably reconstructs the quantity p in the inverse algorithm as compared to the same quantity computed with a forward algorithm. This work generalizes the Log-Elastographic algorithm in [20] which uses one displacement component, is derived assuming the component satisfies the wave equation, and is tested on synthetic data computed with the wave equation model. The 2D Log-Elastographic algorithm is tested on 2D synthetic data and 2D in-vivo data from Mayo Clinic. We also exhibit examples to show that the 2D Log-Elastographic algorithm improves the quality of the recovered images as compared to the Log-Elastographic and Direct Inversion algorithms. PMID:21822349
Calculating tissue shear modulus and pressure by 2D Log-Elastographic methods.
McLaughlin, Joyce R; Zhang, Ning; Manduca, Armando
2010-01-01
Shear modulus imaging, often called elastography, enables detection and characterization of tissue abnormalities. In this paper the data is two displacement components obtained from successive MR or ultrasound data sets acquired while the tissue is excited mechanically. A 2D plane strain elastic model is assumed to govern the 2D displacement, u. The shear modulus, μ, is unknown and whether or not the first Lamé parameter, λ, is known the pressure p = λ∇ · u which is present in the plane strain model cannot be measured and is unreliably computed from measured data and can be shown to be an order one quantity in the units kPa. So here we present a 2D Log-Elastographic inverse algorithm that: (1) simultaneously reconstructs the shear modulus, μ, and p, which together satisfy a first order partial differential equation system, with the goal of imaging μ; (2) controls potential exponential growth in the numerical error; and (3) reliably reconstructs the quantity p in the inverse algorithm as compared to the same quantity computed with a forward algorithm. This work generalizes the Log-Elastographic algorithm in [20] which uses one displacement component, is derived assuming the component satisfies the wave equation, and is tested on synthetic data computed with the wave equation model. The 2D Log-Elastographic algorithm is tested on 2D synthetic data and 2Din-vivo data from Mayo Clinic. We also exhibit examples to show that the 2D Log-Elastographic algorithm improves the quality of the recovered images as compared to the Log-Elastographic and Direct Inversion algorithms. PMID:21822349
NASA Technical Reports Server (NTRS)
Bechert, D. W.
1982-01-01
The generation of instability waves in free shear layers is investigated. The model assumes an infinitesimally thin shear layer shed from a semi-infinite plate which is exposed to sound excitation. The acoustical shear layer excitation by a source further away from the plate edge in the downstream direction is very weak while upstream from the plate edge the excitation is relatively efficient. A special solution is given for the source at the plate edge. The theory is then extended to two streams on both sides of the shear layer having different velocities and densities. Furthermore, the excitation of a shear layer in a channel is calculated. A reference quantity is found for the magnitude of the excited instability waves. For a comparison with measurements, numerical computations of the velocity field outside the shear layer were carried out.
Shear wave logging using guided waves
Winbow, G.A.; Chen, S.T.; Rice, J.A.
1988-09-27
This patent describes a method for acoustically logging an earth formation surrounding a borehole which contains a liquid where the approximate shear wave velocity v of the formation is known. The method consists of: vibrating a dipole source in the liquid to generate in the liquid a guided wave the frequencies of which include a critical frequency f given by zeta = ..nu..12a where a is the borehole radius, so that the fastest component of the guided wave has velocity substantially equal to ..nu..; and detecting the arrival of the fastest component of the guided wave at least one location in the liquid spaced longitudinally along the borehole from the dipole source.
WFR-2D: an analytical model for PWAS-generated 2D ultrasonic guided wave propagation
NASA Astrophysics Data System (ADS)
Shen, Yanfeng; Giurgiutiu, Victor
2014-03-01
This paper presents WaveFormRevealer 2-D (WFR-2D), an analytical predictive tool for the simulation of 2-D ultrasonic guided wave propagation and interaction with damage. The design of structural health monitoring (SHM) systems and self-aware smart structures requires the exploration of a wide range of parameters to achieve best detection and quantification of certain types of damage. Such need for parameter exploration on sensor dimension, location, guided wave characteristics (mode type, frequency, wavelength, etc.) can be best satisfied with analytical models which are fast and efficient. The analytical model was constructed based on the exact 2-D Lamb wave solution using Bessel and Hankel functions. Damage effects were inserted in the model by considering the damage as a secondary wave source with complex-valued directivity scattering coefficients containing both amplitude and phase information from wave-damage interaction. The analytical procedure was coded with MATLAB, and a predictive simulation tool called WaveFormRevealer 2-D was developed. The wave-damage interaction coefficients (WDICs) were extracted from harmonic analysis of local finite element model (FEM) with artificial non-reflective boundaries (NRB). The WFR-2D analytical simulation results were compared and verified with full scale multiphysics finite element models and experiments with scanning laser vibrometer. First, Lamb wave propagation in a pristine aluminum plate was simulated with WFR-2D, compared with finite element results, and verified by experiments. Then, an inhomogeneity was machined into the plate to represent damage. Analytical modeling was carried out, and verified by finite element simulation and experiments. This paper finishes with conclusions and suggestions for future work.
Magnetized stratified rotating shear waves
NASA Astrophysics Data System (ADS)
Salhi, A.; Lehner, T.; Godeferd, F.; Cambon, C.
2012-02-01
solution at infinite vertical wavelength (k3=0): There is an oscillatory behavior for τ>1+|K2/k1|, where τ=St is a dimensionless time and K2 is the radial component of the wave vector at τ=0. The model is suitable to describe instabilities leading to turbulence by the bypass mechanism that can be relevant for the analysis of magnetized stratified Keplerian disks with a purely azimuthal field. For initial isotropic conditions, the time evolution of the spectral density of total energy (kinetic + magnetic + potential) is considered. At k3=0, the vertical motion is purely oscillatory, and the sum of the vertical (kinetic + magnetic) energy plus the potential energy does not evolve with time and remains equal to its initial value. The horizontal motion can induce a rapid transient growth provided K2/k1≫1. This rapid growth is due to the aperiodic velocity vortex mode that behaves like Kh/kh where kh(τ)=[k12+(K2-k1τ)2]1/2 and Kh=kh(0). After the leading phase (τ>K2/k1≫1), the horizontal magnetic energy and the horizontal kinetic energy exhibit a similar (oscillatory) behavior yielding a high level of total energy. The contribution to energies coming from the modes k1=0 and k3=0 is addressed by investigating the one-dimensional spectra for an initial Gaussian dense spectrum. For a magnetized Keplerian disk with a purely vertical field, it is found that an important contribution to magnetic and kinetic energies comes from the region near k1=0. The limit at k1=0 of the streamwise one-dimensional spectra of energies, or equivalently, the streamwise two-dimensional (2D) energy, is then computed. The comparison of the ratios of these 2D quantities with their three-dimensional counterparts provided by previous direct numerical simulations shows a quantitative agreement.
Magnetized stratified rotating shear waves.
Salhi, A; Lehner, T; Godeferd, F; Cambon, C
2012-02-01
stability of the solution at infinite vertical wavelength (k(3) = 0): There is an oscillatory behavior for τ > 1+|K(2)/k(1)|, where τ = St is a dimensionless time and K(2) is the radial component of the wave vector at τ = 0. The model is suitable to describe instabilities leading to turbulence by the bypass mechanism that can be relevant for the analysis of magnetized stratified Keplerian disks with a purely azimuthal field. For initial isotropic conditions, the time evolution of the spectral density of total energy (kinetic + magnetic + potential) is considered. At k(3) = 0, the vertical motion is purely oscillatory, and the sum of the vertical (kinetic + magnetic) energy plus the potential energy does not evolve with time and remains equal to its initial value. The horizontal motion can induce a rapid transient growth provided K(2)/k(1)>1. This rapid growth is due to the aperiodic velocity vortex mode that behaves like K(h)/k(h) where k(h)(τ)=[k(1)(2) + (K(2) - k(1)τ)(2)](1/2) and K(h) =k(h)(0). After the leading phase (τ > K(2)/k(1)>1), the horizontal magnetic energy and the horizontal kinetic energy exhibit a similar (oscillatory) behavior yielding a high level of total energy. The contribution to energies coming from the modes k(1) = 0 and k(3) = 0 is addressed by investigating the one-dimensional spectra for an initial Gaussian dense spectrum. For a magnetized Keplerian disk with a purely vertical field, it is found that an important contribution to magnetic and kinetic energies comes from the region near k(1) = 0. The limit at k(1) = 0 of the streamwise one-dimensional spectra of energies, or equivalently, the streamwise two-dimensional (2D) energy, is then computed. The comparison of the ratios of these 2D quantities with their three-dimensional counterparts provided by previous direct numerical simulations shows a quantitative agreement. PMID:22463311
Explosion Shear Wave Generation and Scattering
NASA Astrophysics Data System (ADS)
Baker, G. E.; Stevens, J. L.; Xu, H.
2004-12-01
We use observations of explosion-generated Lg together with three separate types of numerical models to determine how underground nuclear explosions generate shear wave phases. This question is fundamental to how Lg phases are interpreted for use in explosion yield estimation and earthquake/explosion discrimination. A simple point explosion in a uniform medium generates no shear waves, so the Lg phase is generated entirely by non-spherical components of the source and conversions through reflections and scattering. Our results indicate that the most important sources of high frequency explosion shear waves are P to S conversions at the free surface and S waves generated directly by a realistic distributed explosion source including nonlinear effects due to the free surface and gravity. In addition, Rg scattering may contribute to lower frequency Lg. Near source S is observed on both radial and tangential component records from a diverse set of explosion data. The data sets include 1) Degelen Mountain explosions recorded at distances less than 100 km and corresponding recordings at Borovoye (BOR) at 650 km; 2) recordings from Russian deep seismic sounding experiments; 3) Nevada Test Site (NTS) explosion sources including the Nonproliferation Experiment (NPE) and nuclear tests covering a range of source depths and media properties. We model the overburied NPE, and underburied and overburied Degelen explosions, using point sources and two-dimensional nonlinear finite difference calculations to quantify the source effects. We use energy conservation to determine an upper bound on Rg to Lg scattering. Results indicate that Rg to Lg scattering may be important at frequencies less than 1 Hz, and in Lg coda, but is less than Lg generated directly by the explosion at higher frequencies. We use 2D and 3D finite difference calculations, using the known topography and velocity structure at Degelen Mt. and lateral heterogeneities within the crust, to estimate the effect of
Maximum likelihood estimation of shear wave speed in transient elastography.
Audière, Stéphane; Angelini, Elsa D; Sandrin, Laurent; Charbit, Maurice
2014-06-01
Ultrasonic transient elastography (TE), enables to assess, under active mechanical constraints, the elasticity of the liver, which correlates with hepatic fibrosis stages. This technique is routinely used in clinical practice to assess noninvasively liver stiffness. The Fibroscan system used in this work generates a shear wave via an impulse stress applied on the surface of the skin and records a temporal series of radio-frequency (RF) lines using a single-element ultrasound probe. A shear wave propagation map (SWPM) is generated as a 2-D map of the displacements along depth and time, derived from the correlations of the sequential 1-D RF lines, assuming that the direction of propagation (DOP) of the shear wave coincides with the ultrasound beam axis (UBA). Under the assumption of pure elastic tissue, elasticity is proportional to the shear wave speed. This paper introduces a novel approach to the processing of the SWPM, deriving the maximum likelihood estimate of the shear wave speed when comparing the observed displacements and the estimates provided by the Green's functions. A simple parametric model is used to interface Green's theoretical values of noisy measures provided by the SWPM, taking into account depth-varying attenuation and time-delay. The proposed method was evaluated on numerical simulations using a finite element method simulator and on physical phantoms. Evaluation on this test database reported very high agreements of shear wave speed measures when DOP and UBA coincide. PMID:24835213
Li, Sinan; Cheng, Yi; Eckersley, Robert J; Elson, Daniel S; Tang, Meng-Xing
2015-01-01
Shear wave speed is quantitatively related to tissue viscoelasticity. Previously we reported shear wave tracking at centimetre depths in a turbid optical medium using laser speckle contrast detection. Shear wave progression modulates displacement of optical scatterers and therefore modulates photon phase and changes the laser speckle patterns. Time-resolved charge-coupled device (CCD)-based speckle contrast analysis was used to track shear waves and measure the time-of-flight of shear waves for speed measurement. In this manuscript, we report a new observation of the laser speckle contrast difference signal for dual shear waves. A modulation of CCD speckle contrast difference was observed and simulation reproduces the modulation pattern, suggesting its origin. Both experimental and simulation results show that the dual shear wave approach generates an improved definition of temporal features in the time-of-flight optical signal and an improved signal to noise ratio with a standard deviation less than 50% that of individual shear waves. Results also show that dual shear waves can correct the bias of shear wave speed measurement caused by shear wave reflections from elastic boundaries. PMID:26114021
Stress dynamics of a 2D dense granular system near shear jamming
NASA Astrophysics Data System (ADS)
Ren, Jie; Dijksman, Joshua; Behringer, Robert
2013-03-01
We study the dynamics of pressure and shear stress in a frictional 2D dense granular system using a novel apparatus that can provide fixed-volume shear without generating inhomogeneities. Under increasing shear strain, the system's pressure shows a strong increase with strain, characterized by a ``Reynolds coefficient,'' R =d2 P / dγ2 . R depends only on packing fraction ϕ, and shows a strong increase as ϕ approaches ϕJ from below. In the meantime, the system's shear stress shows a non-monotonic behavior with increasing strain. It first increases with strain as the system is in ``fragile'' states and builds up long force chains along the compression direction. After a certain amount of strain, force chains along the dilation direction starts to build up, and the system transfers into a ``shear-jammed'' state and the shear stress starts to decrease with strain. Under oscillatory shear, both pressure and shear stress show limit-cycle behavior and reach steady states after many cycles. However, the limit cycles of pressure and shear stress are very different: the pressure exhibits a hysteresis-free parabolic curve, while the shear stress exhibits a strongly hysteretic loop. This work is funded by NSF grants: DMR0906908, DMS0835571, NASA grant NNX10AU01G and ARO grant W911NF-11-1-0110.
A new shear wave imaging system for ultrasound elastography.
Qiu, Weibao; Wang, Congzhi; Xiao, Yang; Qian, Ming; Zheng, Hairong
2015-08-01
Ultrasound elastography is able to provide a non-invasive measurement of tissue elasticity properties. Shear wave imaging (SWI) technique is a quantitative method for tissue stiffness assessment. However, traditional SWI implementations cannot acquire 2D quantitative images of tissue elasticity distribution. In this study, a new shear wave imaging system is proposed and evaluated. Detailed delineation of hardware and image processing algorithms are presented. Programmable devices are selected to support flexible control of the system and the image processing algorithms. Analytic signal based cross-correlation method and a Radon transform based shear wave speed determination method are proposed with parallel computation ability. Tissue mimicking phantom imaging, and in vitro imaging measurements are conducted to demonstrate the performance of the proposed system. The system has the ability to provide a new choice for quantitative mapping of the tissue elasticity, and has good potential to be implemented into commercial ultrasound scanner. PMID:26737133
Continuous wave laser for wind shear detection
NASA Technical Reports Server (NTRS)
Nelson, Loren
1991-01-01
Details of the design and development of a continuous-wave heterodyne carbon dioxide laser which has wind shear detection capabilities are given in viewgraph form. The goal of the development was to investigate the lower cost CW (rather than pulsed) lidar option for look-ahead wind shear detection from aircraft. The device has potential utility for ground based wind shear detection at secondary airports where the high cost of a Terminal Doppler Weather Radar system is not justifiable.
AnisWave2D: User's Guide to the 2d Anisotropic Finite-DifferenceCode
Toomey, Aoife
2005-01-06
This document describes a parallel finite-difference code for modeling wave propagation in 2D, fully anisotropic materials. The code utilizes a mesh refinement scheme to improve computational efficiency. Mesh refinement allows the grid spacing to be tailored to the velocity model, so that fine grid spacing can be used in low velocity zones where the seismic wavelength is short, and coarse grid spacing can be used in zones with higher material velocities. Over-sampling of the seismic wavefield in high velocity zones is therefore avoided. The code has been implemented to run in parallel over multiple processors and allows large-scale models and models with large velocity contrasts to be simulated with ease.
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.
Song, Pengfei; Macdonald, Michael; Behler, Russell; Lanning, Justin; Wang, Michael; Urban, Matthew; Manduca, Armando; Zhao, Heng; Callstrom, Matthew; Alizad, Azra; Greenleaf, James; Chen, Shigao
2015-02-01
Two-dimensional shear-wave elastography presents 2-D quantitative shear elasticity maps of tissue, which are clinically useful for both focal lesion detection and diffuse disease diagnosis. Realization of 2-D shear-wave elastography on conventional ultrasound scanners, however, is challenging because of the low tracking pulse-repetition-frequency (PRF) of these systems. Although some clinical and research platforms support software beamforming and plane-wave imaging with high PRF, the majority of current clinical ultrasound systems do not have the software beamforming capability, which presents a critical challenge for translating the 2-D shear-wave elastography technique from laboratory to clinical scanners. To address this challenge, this paper presents a time-aligned sequential tracking (TAST) method for shear-wave tracking on conventional ultrasound scanners. TAST takes advantage of the parallel beamforming capability of conventional systems and realizes high-PRF shear-wave tracking by sequentially firing tracking vectors and aligning shear wave data in the temporal direction. The comb-push ultrasound shear elastography (CUSE) technique was used to simultaneously produce multiple shear wave sources within the field-of-view (FOV) to enhance shear wave SNR and facilitate robust reconstructions of 2-D elasticity maps. TAST and CUSE were realized on a conventional ultrasound scanner. A phantom study showed that the shear-wave speed measurements from the conventional ultrasound scanner were in good agreement with the values measured from other 2-D shear wave imaging technologies. An inclusion phantom study showed that the conventional ultrasound scanner had comparable performance to a state-of-the-art shear-wave imaging system in terms of bias and precision in measuring different sized inclusions. Finally, in vivo case analysis of a breast with a malignant mass, and a liver from a healthy subject demonstrated the feasibility of using the conventional ultrasound
Monitoring of thermal therapy based on shear modulus changes: I. shear wave thermometry.
Arnal, Bastien; Pernot, Mathieu; Tanter, Mickael
2011-02-01
The clinical applicability of high-intensity focused ultrasound (HIFU) for noninvasive therapy is today hampered by the lack of robust and real-time monitoring of tissue damage during treatment. The goal of this study is to show that the estimation of local tissue elasticity from shear wave imaging (SWI) can lead to the 2-D mapping of temperature changes during HIFU treatments. This new concept of shear wave thermometry is experimentally implemented here using conventional ultrasonic imaging probes. HIFU treatment and monitoring were, respectively, performed using a confocal setup consisting of a 2.5-MHz single-element transducer focused at 30 mm on ex vivo samples and an 8-MHz ultrasound diagnostic probe. Thermocouple measurements and ultrasound-based thermometry were used as a gold standard technique and were combined with SWI on the same device. The SWI sequences consisted of 2 successive shear waves induced at different lateral positions. Each wave was created using 100-μs pushing beams at 3 depths. The shear wave propagation was acquired at 17,000 frames/s, from which the elasticity map was recovered. HIFU sonications were interleaved with fast imaging acquisitions, allowing a duty cycle of more than 90%. Elasticity and temperature mapping was achieved every 3 s, leading to realtime monitoring of the treatment. Tissue stiffness was found to decrease in the focal zone for temperatures up to 43°C. Ultrasound-based temperature estimation was highly correlated to stiffness variation maps (r² = 0.91 to 0.97). A reversible calibration phase of the changes of elasticity with temperature can be made locally using sighting shots. This calibration process allows for the derivation of temperature maps from shear wave imaging. Compared with conventional ultrasound-based approaches, shear wave thermometry is found to be much more robust to motion artifacts. PMID:21342822
NASA Astrophysics Data System (ADS)
Pasternack, Gregory B.; Gilbert, Andrew T.; Wheaton, Joseph M.; Buckland, Evan M.
2006-08-01
SummaryResource managers, scientists, government regulators, and stakeholders are considering sophisticated numerical models for managing complex environmental problems. In this study, observations from a river-rehabilitation experiment involving gravel augmentation and spawning habitat enhancement were used to assess sources and magnitudes of error in depth, velocity, and shear velocity predictions made at the 1-m scale with a commercial two-dimensional (depth-averaged) model. Error in 2D model depth prediction averaged 21%. This error was attributable to topographic survey resolution, which at 1 pt per 1.14 m 2, was inadequate to resolve small humps and depressions influencing point measurements. Error in 2D model velocity prediction averaged 29%. More than half of this error was attributable to depth prediction error. Despite depth and velocity error, 56% of tested 2D model predictions of shear velocity were within the 95% confidence limit of the best field-based estimation method. Ninety percent of the error in shear velocity prediction was explained by velocity prediction error. Multiple field-based estimates of shear velocity differed by up to 160%, so the lower error for the 2D model's predictions suggests such models are at least as accurate as field measurement. 2D models enable detailed, spatially distributed estimates compared to the small number measurable in a field campaign of comparable cost. They also can be used for design evaluation. Although such numerical models are limited to channel types adhering to model assumptions and yield predictions only accurate to ˜20-30%, they can provide a useful tool for river-rehabilitation design and assessment, including spatially diverse habitat heterogeneity as well as for pre- and post-project appraisal.
Breakdown of wave diffusion in 2D due to loops.
Haney, Matthew; Snieder, Roel
2003-08-29
The validity of the diffusion approximation for the intensity of multiply scattered waves is tested with numerical simulations in a strongly scattering 2D medium of finite extent. We show that the diffusion equation underestimates the intensity and attribute this to both the neglect of recurrent scattering paths and interference within diffusion theory. We present a theory to quantify this discrepancy based on counting all possible scattering paths between point scatterers. Interference phenomena, due to loop paths, are incorporated in a way similar to coherent backscattering. PMID:14525183
Coded Excitation Plane Wave Imaging for Shear Wave Motion Detection
Song, Pengfei; Urban, Matthew W.; Manduca, Armando; Greenleaf, James F.; Chen, Shigao
2015-01-01
Plane wave imaging has greatly advanced the field of shear wave elastography thanks to its ultrafast imaging frame rate and the large field-of-view (FOV). However, plane wave imaging also has decreased penetration due to lack of transmit focusing, which makes it challenging to use plane waves for shear wave detection in deep tissues and in obese patients. This study investigated the feasibility of implementing coded excitation in plane wave imaging for shear wave detection, with the hypothesis that coded ultrasound signals can provide superior detection penetration and shear wave signal-to-noise-ratio (SNR) compared to conventional ultrasound signals. Both phase encoding (Barker code) and frequency encoding (chirp code) methods were studied. A first phantom experiment showed an approximate penetration gain of 2-4 cm for the coded pulses. Two subsequent phantom studies showed that all coded pulses outperformed the conventional short imaging pulse by providing superior sensitivity to small motion and robustness to weak ultrasound signals. Finally, an in vivo liver case study on an obese subject (Body Mass Index = 40) demonstrated the feasibility of using the proposed method for in vivo applications, and showed that all coded pulses could provide higher SNR shear wave signals than the conventional short pulse. These findings indicate that by using coded excitation shear wave detection, one can benefit from the ultrafast imaging frame rate and large FOV provided by plane wave imaging while preserving good penetration and shear wave signal quality, which is essential for obtaining robust shear elasticity measurements of tissue. PMID:26168181
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
Seismic shear waves as Foucault pendulum
NASA Astrophysics Data System (ADS)
Snieder, Roel; Sens-Schönfelder, Christoph; Ruigrok, Elmer; Shiomi, Katsuhiko
2016-03-01
Earth's rotation causes splitting of normal modes. Wave fronts and rays are, however, not affected by Earth's rotation, as we show theoretically and with observations made with USArray. We derive that the Coriolis force causes a small transverse component for P waves and a small longitudinal component for S waves. More importantly, Earth's rotation leads to a slow rotation of the transverse polarization of S waves; during the propagation of S waves the particle motion behaves just like a Foucault pendulum. The polarization plane of shear waves counteracts Earth's rotation and rotates clockwise in the Northern Hemisphere. The rotation rate is independent of the wave frequency and is purely geometric, like the Berry phase. Using the polarization of ScS and ScS2 waves, we show that the Foucault-like rotation of the S wave polarization can be observed. This can affect the determination of source mechanisms and the interpretation of observed SKS splitting.
NASA Astrophysics Data System (ADS)
Yamakoshi, Yoshiki; Yamamoto, Atsushi; Kasahara, Toshihiro; Iijima, Tomohiro; Yuminaka, Yasushi
2015-07-01
We have proposed a quantitative shear wave imaging technique for continuous shear wave excitation. Shear wave wavefront is observed directly by color flow imaging using a general-purpose ultrasonic imaging system. In this study, the proposed method is applied to experiments in vivo, and shear wave maps, namely, the shear wave phase map, which shows the shear wave propagation inside the medium, and the shear wave velocity map, are observed for the skeletal muscle in the shoulder. To excite the shear wave inside the skeletal muscle of the shoulder, a hybrid ultrasonic wave transducer, which combines a small vibrator with an ultrasonic wave probe, is adopted. The shear wave velocity of supraspinatus muscle, which is measured by the proposed method, is 4.11 ± 0.06 m/s (N = 4). This value is consistent with those obtained by the acoustic radiation force impulse method.
Estimation of seabed shear-wave velocity profiles using shear-wave source data.
Dong, Hefeng; Nguyen, Thanh-Duong; Duffaut, Kenneth
2013-07-01
This paper estimates seabed shear-wave velocity profiles and their uncertainties using interface-wave dispersion curves extracted from data generated by a shear-wave source. The shear-wave source generated a seismic signature over a frequency range between 2 and 60 Hz and was polarized in both in-line and cross-line orientations. Low-frequency Scholte- and Love-waves were recorded. Dispersion curves of the Scholte- and Love-waves for the fundamental mode and higher-order modes are extracted by three time-frequency analysis methods. Both the vertically and horizontally polarized shear-wave velocity profiles in the sediment are estimated by the Scholte- and Love-wave dispersion curves, respectively. A Bayesian approach is utilized for the inversion. Differential evolution, a global search algorithm is applied to estimate the most-probable shear-velocity models. Marginal posterior probability profiles are computed by Metropolis-Hastings sampling. The estimated vertically and horizontally polarized shear-wave velocity profiles fit well with the core and in situ measurements. PMID:23862796
2D modeling of electromagnetic waves in cold plasmas
Crombé, K.; Van Eester, D.; Koch, R.; Kyrytsya, V.
2014-02-12
The consequences of sheath (rectified) electric fields, resulting from the different mobility of electrons and ions as a response to radio frequency (RF) fields, are a concern for RF antenna design as it can cause damage to antenna parts, limiters and other in-vessel components. As a first step to a more complete description, the usual cold plasma dielectric description has been adopted, and the density profile was assumed to be known as input. Ultimately, the relevant equations describing the wave-particle interaction both on the fast and slow timescale will need to be tackled but prior to doing so was felt as a necessity to get a feeling of the wave dynamics involved. Maxwell's equations are solved for a cold plasma in a 2D antenna box with strongly varying density profiles crossing also lower hybrid and ion-ion hybrid resonance layers. Numerical modelling quickly becomes demanding on computer power, since a fine grid spacing is required to capture the small wavelengths effects of strongly evanescent modes.
Shear wave speed and dispersion measurements using crawling wave chirps.
Hah, Zaegyoo; Partin, Alexander; Parker, Kevin J
2014-10-01
This article demonstrates the measurement of shear wave speed and shear speed dispersion of biomaterials using a chirp signal that launches waves over a range of frequencies. A biomaterial is vibrated by two vibration sources that generate shear waves inside the medium, which is scanned by an ultrasound imaging system. Doppler processing of the acquired signal produces an image of the square of vibration amplitude that shows repetitive constructive and destructive interference patterns called "crawling waves." With a chirp vibration signal, successive Doppler frames are generated from different source frequencies. Collected frames generate a distinctive pattern which is used to calculate the shear speed and shear speed dispersion. A special reciprocal chirp is designed such that the equi-phase lines of a motion slice image are straight lines. Detailed analysis is provided to generate a closed-form solution for calculating the shear wave speed and the dispersion. Also several phantoms and an ex vivo human liver sample are scanned and the estimation results are presented. PMID:24658144
Fan-structure waves in shear ruptures
NASA Astrophysics Data System (ADS)
Tarasov, Boris
2016-04-01
This presentation introduces a recently identified shear rupture mechanism providing a paradoxical feature of hard rocks - the possibility of shear rupture propagation through the highly confined intact rock mass at shear stress levels significantly less than frictional strength. According to the fan-mechanism the shear rupture propagation is associated with consecutive creation of small slabs in the fracture tip which, due to rotation caused by shear displacement of the fracture interfaces, form a fan-structure representing the fracture head. The fan-head combines such unique features as: extremely low shear resistance (below the frictional strength), self-sustaining stress intensification in the rupture tip (providing easy formation of new slabs), and self-unbalancing conditions in the fan-head (making the failure process inevitably spontaneous and violent). An important feature of the fan-mechanism is the fact that for the initial formation of the fan-structure an enhanced local shear stress is required, however, after completion of the fan-structure it can propagate as a dynamic wave through intact rock mass at shear stresses below the frictional strength. Paradoxically low shear strength of pristine rocks provided by the fan-mechanism determines the correspondingly low transient strength of the lithosphere, which favours generation of new earthquake faults in the intact rock mass adjoining pre-existing faults in preference to frictional stick-slip instability along these faults. The new approach reveals an alternative role of pre-existing faults in earthquake activity: they represent local stress concentrates in pristine rock adjoining the fault where special conditions for the fan-mechanism nucleation are created, while further dynamic propagation of the new fault (earthquake) occurs at low field stresses even below the frictional strength.
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.
Spatial correlation of shear-wave velocity within San Francisco Bay Sediments
Thompson, E.M.; Baise, L.G.; Kayen, R.E.
2006-01-01
Sediment properties are spatially variable at all scales, and this variability at smaller scales influences high frequency ground motions. We show that surface shear-wave velocity is highly correlated within San Francisco Bay Area sediments using shear-wave velocity measurements from 210 seismic cone penetration tests. We use this correlation to estimate the surface sediment velocity structure using geostatistics. We find that the variance of the estimated shear-wave velocity is reduced using ordinary kriging, and that including this velocity structure in 2D ground motion simulations of a moderate sized earthquake improves the accuracy of the synthetics. Copyright ASCE 2006.
Shear wave velocities in the earth's mantle.
NASA Technical Reports Server (NTRS)
Robinson, R.; Kovach, R. L.
1972-01-01
Direct measurement of the travel time gradient for S waves together with travel time data are used to derive a shear velocity model for the earth's mantle. In order to satisfy the data it is necessary to discard the usual assumption of lateral homogeneity below shallow depths. A shear velocity differential is proposed for a region between western North America and areas of the Pacific Ocean. Distinctive features of the velocity model for the upper mantle beneath western North America are a low-velocity zone centered at 100 km depth and zones of high velocity gradient beginning at 400, 650, and 900 km.
Waves in Turbulent Stably Stratified Shear Flow
NASA Technical Reports Server (NTRS)
Jacobitz, F. G.; Rogers, M. M.; Ferziger, J. H.; Parks, John W. (Technical Monitor)
2002-01-01
Two approaches for the identification of internal gravity waves in sheared and unsheared homogeneous stratified turbulence are investigated. First, the phase angle between the vertical velocity and density fluctuations is considered. It was found, however, that a continuous distribution of the phase angle is present in weakly and strongly stratified flow. Second, a projection onto the solution of the linearized inviscid equations of motion of unsheared stratified flow is investigated. It was found that a solution of the fully nonlinear viscous Navier-Stokes equations can be represented by the linearized inviscid solution. The projection yields a decomposition into vertical wave modes and horizontal vortical modes.
ML shear wave velocity tomography for the Iranian Plateau
NASA Astrophysics Data System (ADS)
Maheri-Peyrov, Mehdi; Ghods, Abdolreza; Abbasi, Madjid; Bergman, Eric; Sobouti, Farhad
2016-04-01
Iranian Plateau reflects several different tectonic styles of collision, and large-scale strike-slip faults. We calculate a high-resolution 2-D ML shear velocity map for the Iranian Plateau to detect lateral crustal thickness changes associated with different tectonic boundaries. The ML velocity is very sensitive to strong lateral variations of crustal thickness and varies between the velocity of Lg and Sn phases. Our data set consists of 65 795 ML amplitude velocity measurements from 2531 precisely relocated events recorded by Iranian networks in the period 1996-2014. Using a constrained least-squares inversion scheme, we inverted the ML velocities for a 2-D shear velocity map of Iran. Our results show that the Zagros and South Caspian Basin (SCB) have shear wave velocities close to the Sn phase, and are thus Lg-blocking regions. High velocities in the High Zagros and the Simply Folded Belt imply significant crustal undulations within these zones. We note that in the central and south Zagros, the velocity border between the Zagros and central Iran is not coincident with the Zagros suture line that marks underthrusting of the Arabian plate beneath central Iran. The low plains of Gilan and Gorgan to the south of the Caspian Sea show high shear velocities similar to the SCB, implying that they are either underlain by an oceanic type crust or a transitional crust with a strong lateral crustal thickness gradient. The Lut block is an Lg-passing block implying that it is not surrounded by any sudden crustal thickness changes along its borders with central Iran. In the Alborz, NW Iran, Kopeh-Dagh, Binalud and most of the central Iran, low shear velocity near the Lg velocity is attributed to smooth or minor Moho undulations within these regions.
Horizontal Shear Wave Imaging of Large Optics
Quarry, M J
2007-09-05
When complete the National Ignition Facility (NIF) will be the world's largest and most energetic laser and will be capable of achieving for the first time fusion ignition in the laboratory. Detecting optics features within the laser beamlines and sizing them at diameters of 0.1 mm to 10 mm allows timely decisions concerning refurbishment and will help with the routine operation of the system. Horizontally polarized shear waves at 10 MHz were shown to accurately detect, locate, and size features created by laser operations from 0.5 mm to 8 mm by placing sensors at the edge of the optic. The shear wave technique utilizes highly directed beams. The outer edge of an optic can be covered with shear wave transducers on four sides. Each transducer sends a pulse into the optic and any damage reflects the pulse back to the transmitter. The transducers are multiplexed, and the collected time waveforms are enveloped and replicated across the width of the element. Multiplying the data sets from four directions produces a map of reflected amplitude to the fourth power, which images the surface of the optic. Surface area can be measured directly from the image, and maximum depth was shown to be correlated to maximum amplitude of the reflected waveform.
The stability of Rossby waves in a stratified shear fluid
NASA Astrophysics Data System (ADS)
Tan, Benkui
1990-11-01
An investigation is undertaken of the stability of linear Rossby waves in a stratified shear fluid by means of a qualitative theory employing ordinary differential equations. It is noted that, while the basic current has no detectable shear, the Rossby waves are always stable. If the basic current possesses only horizontal shear, the unstable criterion for waves takes one form, but it takes entirely another in the case where the basic current possesses only vertical shear.
Wavefield Analysis of Rayleigh Waves for Near-Surface Shear-Wave Velocity
NASA Astrophysics Data System (ADS)
Zeng, Chong
2011-12-01
Shear (S)-wave velocity is a key property of near-surface materials and is the fundamental parameter for many environmental and engineering geophysical studies. Directly acquiring accurate S-wave velocities from a seismic shot gather is usually difficult due to the poor signal-to-noise ratio. The relationship between Rayleigh-wave phase velocity and frequency has been widely utilized to estimate the S-wave velocities in shallow layers using the multichannel analysis of surface waves (MASW) technique. Hence, Rayleigh wave is a main focus of most near-surface seismic studies. Conventional dispersion analysis of Rayleigh waves assumes that the earth is laterally homogeneous and the free surface is horizontally flat, which limits the application of surface-wave methods to only 1D earth models or very smooth 2D models. In this study I extend the analysis of Rayleigh waves to a 2D domain by employing the 2D full elastic wave equation so as to address the lateral heterogeneity problem. I first discuss the accurate simulation of Rayleigh waves through finite-difference method and the boundary absorbing problems in the numerical modeling with a high Poisson's ratio (> 0.4), which is a unique near-surface problem. Then I develop an improved vacuum formulation to generate accurate synthetic seismograms focusing on Rayleigh waves in presence of surface topography and internal discontinuities. With these solutions to forward modeling of Rayleigh waves, I evaluate the influence of surface topography to conventional dispersion analysis in 2D and 3D domains by numerical investigations. At last I examine the feasibility of inverting waveforms of Rayleigh waves for shallow S-wave velocities using a genetic algorithm. Results of the study show that Rayleigh waves can be accurately simulated in near surface using the improved vacuum formulation. Spurious reflections during the numerical modeling can be efficiently suppressed by the simplified multiaxial perfectly matched layers. The
Cascade properties of shear Alfven wave turbulence
NASA Technical Reports Server (NTRS)
Bondeson, A.
1985-01-01
Nonlinear three-wave interactions of linear normal modes are investigated for two-dimensional incompressible magnetohydrodynamics and the weakly three-dimensional Strauss equations in the case where a strong uniform background field B0 is present. In both systems the only resonant interaction affecting Alfven waves is caused by the shear of the background field plus the zero frequency components of the perturbation. It is shown that the Alfven waves are cascaded in wavenumber space by a mechanism equivalent to the resonant absorption at the Alfven resonance. For large wavenumbers perpendicular to B0, the cascade is described by Hamilton's ray equations, dk/dt = -(first-order) partial derivative of omega with respect to vector r, where omega includes the effects of the zero frequency perturbations.
Ion acoustic wave collapse via two-ion wave decay: 2D Vlasov simulation and theory
NASA Astrophysics Data System (ADS)
Chapman, Thomas; Berger, Richard; Banks, Jeffrey; Brunner, Stephan
2015-11-01
The decay of ion acoustic waves (IAWs) via two-ion wave decay may transfer energy from the electric field of the IAWs to the particles, resulting in a significant heating of resonant particles. This process has previously been shown in numerical simulations to decrease the plasma reflectivity due to stimulated Brillouin scattering. Two-ion wave decay is a fundamental property of ion acoustic waves that occurs over most if not all of the parameter space of relevance to inertial confinement fusion experiments, and can lead to a sudden collapse of IAWs. The treatment of all species kinetically, and in particular the electrons, is required to describe the decay process correctly. We present fully kinetic 2D+2V Vlasov simulations of IAWs undergoing decay to a highly nonlinear turbulent state using the code LOKI. The scaling of the decay rate with characteristic plasma parameters and wave amplitude is shown. A new theory describing two-ion wave decay in 2D, that incorporates key kinetic properties of the electrons, is presented and used to explain quantitatively for the first time the observed decay of IAWs. Work performed under auspices of U.S. DoE by LLNL, Contract DE-AC52-07NA2734. Funded by LDRD 15-ERD-038 and supported by LLNL Grand Challenge allocation.
Focusing surface wave imaging with flexible 2D array
NASA Astrophysics Data System (ADS)
Zhou, Shiyuan; Fu, Junqiang; Li, Zhe; Xu, Chunguang; Xiao, Dingguo; Wang, Shaohan
2016-04-01
Curved surface is widely exist in key parts of energy and power equipment, such as, turbine blade cylinder block and so on. Cycling loading and harsh working condition of enable fatigue cracks appear on the surface. The crack should be found in time to avoid catastrophic damage to the equipment. A flexible 2D array transducer was developed. 2D Phased Array focusing method (2DPA), Mode-Spatial Double Phased focusing method (MSDPF) and the imaging method using the flexible 2D array probe are studied. Experiments using these focusing and imaging method are carried out. Surface crack image is obtained with both 2DPA and MSDPF focusing method. It have been proved that MSDPF can be more adaptable for curved surface and more calculate efficient than 2DPA.
Hammering Yucca Flat, Part Two: Shear-Wave Velocity
NASA Astrophysics Data System (ADS)
Finlay, T. S.; Abbott, R. E.; Knox, H. A.; Tang, D. G.; James, S. R.; Haney, M. M.; Hampshire, J. B., II
2015-12-01
In preparation for the next phase of the Source Physics Experiment (SPE), we conducted an active-source seismic survey of Yucca Flat, Nevada, on the Nevada National Security Site. Results from this survey will be used to inform the geologic models associated with the SPE project. For this study, we used a novel 13,000 kilogram weight-drop seismic source to interrogate an 18-km North-South transect of Yucca Flat. Source points were spaced every 200 meters and were recorded by 350 to 380 3-component 2-Hz geophones with variable spacings of 10, 20, and 100 meters. We utilized the Refraction-Microtremor (ReMi) technique to create multiple 1D dispersion curves, which were then inverted for shear-wave velocity profiles using the Dix inversion method (Tsai and Haney, 2015). Each of these 1D velocity models was subsequently stitched together to create a 2D profile over the survey area. The dispersion results indicate a general decrease in surface-wave phase velocity to the south. This result is supported by slower shear-wave velocity sediments and increasing basin depth towards the survey's southern extent. Sandia National Laboratories is a multi-program laboratory managed and operated by Sandia Corporation, a wholly owned subsidiary of Lockheed Martin Corporation, for the U.S. Department of Energy's National Nuclear Security Administration under contract DE-AC04-94AL85000.
Shear waves in inhomogeneous, compressible fluids in a gravity field.
Godin, Oleg A
2014-03-01
While elastic solids support compressional and shear waves, waves in ideal compressible fluids are usually thought of as compressional waves. Here, a class of acoustic-gravity waves is studied in which the dilatation is identically zero, and the pressure and density remain constant in each fluid particle. These shear waves are described by an exact analytic solution of linearized hydrodynamics equations in inhomogeneous, quiescent, inviscid, compressible fluids with piecewise continuous parameters in a uniform gravity field. It is demonstrated that the shear acoustic-gravity waves also can be supported by moving fluids as well as quiescent, viscous fluids with and without thermal conductivity. Excitation of a shear-wave normal mode by a point source and the normal mode distortion in realistic environmental models are considered. The shear acoustic-gravity waves are likely to play a significant role in coupling wave processes in the ocean and atmosphere. PMID:24606251
From supersonic shear wave imaging to full-field optical coherence shear wave elastography
NASA Astrophysics Data System (ADS)
Nahas, Amir; Tanter, Mickaël; Nguyen, Thu-Mai; Chassot, Jean-Marie; Fink, Mathias; Claude Boccara, A.
2013-12-01
Elasticity maps of tissue have proved to be particularly useful in providing complementary contrast to ultrasonic imaging, e.g., for cancer diagnosis at the millimeter scale. Optical coherence tomography (OCT) offers an endogenous contrast based on singly backscattered optical waves. Adding complementary contrast to OCT images by recording elasticity maps could also be valuable in improving OCT-based diagnosis at the microscopic scale. Static elastography has been successfully coupled with full-field OCT (FF-OCT) in order to realize both micrometer-scale sectioning and elasticity maps. Nevertheless, static elastography presents a number of drawbacks, mainly when stiffness quantification is required. Here, we describe the combination of two methods: transient elastography, based on speed measurements of shear waves induced by ultrasonic radiation forces, and FF-OCT, an en face OCT approach using an incoherent light source. The use of an ultrafast ultrasonic scanner and an ultrafast camera working at 10,000 to 30,000 images/s made it possible to follow shear wave propagation with both modalities. As expected, FF-OCT is found to be much more sensitive than ultrafast ultrasound to tiny shear vibrations (a few nanometers and micrometers, respectively). Stiffness assessed in gel phantoms and an ex vivo rat brain by FF-OCT is found to be in good agreement with ultrasound shear wave elastography.
2-D traveling-wave patterns in binary fluid convection
Surko, C.M.; Porta, A.L.
1996-12-31
An overview is presented of recent experiments designed to study two-dimensional traveling-wave convection in binary fluid convection in a large aspect ratio container. Disordered patterns are observed when convection is initiated. As time proceeds, they evolve to more ordered patterns, consisting of several domains of traveling-waves separated by well-defined domain boundaries. The detailed character of the patterns depends sensitively on the Rayleigh number. Numerical techniques are described which were developed to provide a quantitative characterization of the traveling-wave patterns. Applications of complex demodulation techniques are also described, which make a detailed study of the structure and dynamics of the domain boundaries possible.
Wave-forced reconfiguration of a 2D artificial canopy
NASA Astrophysics Data System (ADS)
Barsu, Sylvie; Doppler, Delphine; Rivière, Nicolas; Lance, Michel
2015-11-01
Blades inside aquatic vegetation canopies show collective motion when submitted to a water flow. Coherent deformation waves might be observed under given flow conditions, which might enhance mass and sediment transfers between the canopy and surrounding flow, thus impacting the plants development. However, most studies have been focused on the flow velocity while the cover motion has been far less studied. Here we present experimental results about the dynamic reconfiguration of a single array of PVC blades in a wave flume. The oscillations of the blades are imaged while the water level is separately measured using resistive probes. A delayed coherent wave motion is observed within the canopy, as a response to the oscillatory flow. The associated transfer function (amplitude, phase, wave speed) is built by correlating blade displacements and water local velocity time series. The canopy-flow interaction is then modelled by a simple linear damped oscillator chain whose parameters are deduced from experiments.
Standing shear waves in anisotropic viscoelastic media
NASA Astrophysics Data System (ADS)
Krit, T.; Golubkova, I.; Andreev, V.
2015-10-01
We studied standing shear waves in anisotropic resonator represented by a rectangular parallelepiped (layer) fixed without slipping between two wooden plates of finite mass. The viscoelastic layer with edges of 70 mm × 40 mm × 15 mm was made of a rubber-like polymer plastisol with rubber bands inside. The bands were placed vertical between the top and the bottom plate. Mechanical properties of the plastisol itself were carefully measured previously. It was found that plastisol shows a cubic nonlinear behavior, i.e. the stress-strain curve could be represented as: σ = μɛ + βμɛ3, where ɛ stands for shear strain and σ is an applied shear stress. The value of shear modulus μ depends on frequency and was found to be several kilopascals which is common for such soft solids. Nonlinear parameter β is frequency dependent too and varies in range from tenths to unity at 1-100 Hz frequency range, decreasing with frequency growth. Stretching the rubber bands inside the layer leads to change of elastic properties in resonator. Such effect could be noticed due to frequency response of the resonator. The numerical model of the resonator was based on finite elements method (FEM) and performed in MatLab. The resonator was cut in hundreds of right triangular prisms. Each prism was provided with viscoelastic properties of the layer except for the top prisms provided with the wooden plate properties and the prisms at the site of the rubber bands provided with the rubber properties. The boundary conditions on each prism satisfied the requirements that resonator is inseparable and all its boundaries but bottom are free. The bottom boundary was set to move horizontally with constant acceleration amplitude. It was shown numerically that the resonator shows anisotropic behavior expressed in different frequency response to oscillations applied to a bottom boundary in different directions.
Shear Wave Splitting Beneath the Galapagos Archipelago
NASA Astrophysics Data System (ADS)
Fontaine, F. R.; Burkett, P. G.; Hooft, E. E.; Toomey, D. R.; Solomon, S. C.; Silver, P. G.
2004-12-01
We report measurements of teleseismic shear wave splitting in the Galápagos Archipelago. The inferred lateral variations in azimuthal anisotropy allow us to examine the dynamics of an evolving hotspot-ridge system. The data are from SKS and SKKS phases, as well as S waves from deep sources, recorded by a relatively dense network of 10 portable broadband seismometers deployed from 1999 to 2003 for the IGUANA (Imaging Galápagos Upwelling and Neotectonics of the Archipelago) experiment and from the GSN broadband station in Santa Cruz (PAYG). We find a delay time between fast and slow shear waves of 0.4 to 0.9 s and fast polarization directions of N85-90° E beneath five stations at the leading and southern edge of the archipelago. Despite clear seismic signals, we did not find any anisotropy at the six stations located in the interior of the archipelago. For those stations that show shear wave splitting, there is an increase in the delay time toward the expected location of the Galápagos hotspot at the western edge of the archipelago. With the exception of Española, fast polarization directions (N85-90° E) are close to the current direction of absolute plate motion of the overlying Nazca plate (N91° E). The lack of azimuthal anisotropy in the interior of the archipelago is interpreted as an absence of strongly oriented mantle fabric beneath these stations. The apparent isotropy in this dynamic region, where we expect considerable mantle strain, is surprising. It is not likely that the olivine a-axis is oriented vertically beneath the interior of the archipelago as the Galápagos plume is thought to lie at the western edge. It is also unlikely that there are two layers of perpendicularly-oriented anisotropy which are solely confined to the center of the archipelago. However, there appears to be some correlation between the region of apparent isotropy and a zone of anomalously low upper mantle velocities imaged beneath Santiago and Marchena from surface waves by
Origin of crustal anisotropy: Shear wave splitting studies in Japan
Kaneshima, Satoshi )
1990-07-10
Shear wave splitting manifested as leading shear wave polarization, that is, parallel alignment of leading shear wave particle motions from a variety of sources, has been observed at a number of seismograph stations in Japan. Detected on shear wave seismograms from crustal earthquakes over a wide range of source zones and source-receiver azimuths, the shear wave splitting can be attributed to crustal anisotropy. This paper discusses the relation between leading shear wave polarization directions and tectonic features of Japan. To explain the observed shear wave splitting, the author proposes that at least three phenomena should be taken into account: stress-induced microcracks primarily aligned in vertical or subvertical planes; cracks or fractures in the vicinity of active faults having their orientation parallel to the fault planes; and intrinsic rock anisotropy resulting from preferred orientation of minerals. Travel time differences between leading and slower split shear waves from crustal and upper mantle earthquakes analyzed for about one third of the stations suggest that the crustal anisotropy which causes the observed shear wave splitting may be limited to the upper 15-25 km. This implies that the density of nonhorizontally aligned cracks or fractures below 15-25 km and into the upper mantle is much smaller than that in the crust above 15-25 km.
A proposed experimental test to distinguish waves from 2-D turbulence
NASA Technical Reports Server (NTRS)
Dewan, E. M.
1986-01-01
A theory of buoyancy range turbulence that leads to a unique scale, K sub B, that allows one to differentiate between waves and turbulence for the special case of theta = 0 (i.e., horizontally propagating waves) is discussed. The theory does not seem to lead to a practical empirical distinction for the general situation. This is due to the fact that, as theta is increased, one has the ever-increasing presence of BRT for longer wavelengths. The fact that the numerical values of epsilon prime are not yet available compounds the difficulty. In addition, it does not appear possible to encompass true 2-D turbulence in the theory. We are thus driven to a test which circumvents all these difficulties. A proposed test is based on the idea that waves are coherent and propagate, while in turbulence we have the opposite situation. In particular, the test is suggested by the following quotation from MULLER (1984), on the nature of such turbulence: The turbulence in each horizontal plane is independent from the turbulence in the other planes. If this statement were to be taken literally, it would imply that the temporal coherence between horizontal speeds, separated only in altitude, would be zero. Any vertical separation would be forced to take into account the effects of viscosity: that is to say, a specific finite vertical separation would be needed to destroy coherence. In order to estimate this distance, L, one can use L = C(v/S) (1/2) were v is the kinematic viscosity, S is the shear scale, and C is a constant of order unity.
Shear wavelength estimation based on inverse filtering and multiple-point shear wave generation
NASA Astrophysics Data System (ADS)
Kitazaki, Tomoaki; Kondo, Kengo; Yamakawa, Makoto; Shiina, Tsuyoshi
2016-07-01
Elastography provides important diagnostic information because tissue elasticity is related to pathological conditions. For example, in a mammary gland, higher grade malignancies yield harder tumors. Estimating shear wave speed enables the quantification of tissue elasticity imaging using time-of-flight. However, time-of-flight measurement is based on an assumption about the propagation direction of a shear wave which is highly affected by reflection and refraction, and thus might cause an artifact. An alternative elasticity estimation approach based on shear wavelength was proposed and applied to passive configurations. To determine the elasticity of tissue more quickly and more accurately, we proposed a new method for shear wave elasticity imaging that combines the shear wavelength approach and inverse filtering with multiple shear wave sources induced by acoustic radiation force (ARF). The feasibility of the proposed method was verified using an elasticity phantom with a hard inclusion.
2D wave-front shaping in optical superlattices using nonlinear volume holography.
Yang, Bo; Hong, Xu-Hao; Lu, Rong-Er; Yue, Yang-Yang; Zhang, Chao; Qin, Yi-Qiang; Zhu, Yong-Yuan
2016-07-01
Nonlinear volume holography is employed to realize arbitrary wave-front shaping during nonlinear processes with properly designed 2D optical superlattices. The concept of a nonlinear polarization wave in nonlinear volume holography is investigated. The holographic imaging of irregular patterns was performed using 2D LiTaO_{3} crystals with fundamental wave propagating along the spontaneous polarization direction, and the results agree well with the theoretical predictions. This Letter not only extends the application area of optical superlattices, but also offers an efficient method for wave-front shaping technology. PMID:27367067
Application of a laser/EMAT system for using shear and LS mode converted waves.
Murray, P R; Dewhurst, R J
2002-05-01
Quantitative time-of-flight analysis of laser-generated shear waves and longitudinal-shear mode-converted waves has demonstrated an effective method for non-contact monitoring of the thickness of metal plates. Q-switched Nd:YAG laser pulses with energies of approximately 18 mJ, delivered to the material surface via an optical fibre and focused to a line source by a cylindrical lens, excited surface waves, longitudinal and shear waves. Bulk waves propagated through the plate to be reflected from the far surface. Returning waves were detected using an electro-magnetic acoustic transducer (EMAT) sensitive to in-plane motion. The compilation of B-scans generated as the sensor head was moved along the material's surface to produce a 2-D intensity profile made any changes in the plate thickness easy to visualise. The longitudinal-shear (L-S) and shear-longitudinal (S-L) mode-converted waves provided a method of simultaneously monitoring two different points on the far surface enabling any changes in the material thickness to be clearly identified. This method was used to determine the thickness of aluminium samples ranging in from 5 to 70 mm. PMID:12160043
Piezoelectric shear wave resonator and method of making same
Wang, J.S.; Lakin, K.M.; Landin, A.R.
1985-05-20
An acoustic shear wave resonator comprising a piezoelectric film having its C-axis substantially inclined from the film normal such that the shear wave coupling coefficient significantly exceeds the longitudinal wave coupling coefficient, whereby the film is capable of shear wave resonance, and means for exciting said film to resonate. The film is prepared by deposition in a dc planar magnetron sputtering system to which a supplemental electric field is applied. The resonator structure may also include a semiconductor material having a positive temperature coefficient of resonance such that the resonator has a temperature coefficient of resonance approaching 0 ppM//sup 0/C.
Piezoelectric shear wave resonator and method of making same
Wang, Jin S.; Lakin, Kenneth M.; Landin, Allen R.
1988-01-01
An acoustic shear wave resonator comprising a piezoelectric film having its C-axis substantially inclined from the film normal such that the shear wave coupling coefficient significantly exceeds the longitudinal wave coupling coefficient, whereby the film is capable of shear wave resonance, and means for exciting said film to resonate. The film is prepared by deposition in a dc planar magnetron sputtering system to which a supplemental electric field is applied. The resonator structure may also include a semiconductor material having a positive temperature coefficient of resonance such that the resonator has a temperature coefficient of resonance approaching 0 ppm/.degree.C.
Piezoelectric shear wave resonator and method of making same
Wang, J.S.; Lakin, K.M.; Landin, A.R.
1983-10-25
An acoustic shear wave resonator comprising a piezoelectric film having its C-axis substantially inclined from the film normal such that the shear wave coupling coefficient significantly exceeds the longitudinal wave coupling coefficient, whereby the film is capable of shear wave resonance, and means for exciting said film to resonate. The film is prepared by deposition in a dc planar magnetron sputtering system to which a supplemental electric field is applied. The resonator structure may also include a semiconductor material having a positive temperature coefficient of resonance such that the resonator has a temperature coefficient of resonance approaching 0 ppM//sup 0/C.
Method of making a piezoelectric shear wave resonator
Wang, Jin S.; Lakin, Kenneth M.; Landin, Allen R.
1987-02-03
An acoustic shear wave resonator comprising a piezoelectric film having its C-axis substantially inclined from the film normal such that the shear wave coupling coefficient significantly exceeds the longitudinal wave coupling coefficient, whereby the film is capable of shear wave resonance, and means for exciting said film to resonate. The film is prepared by deposition in a dc planar magnetron sputtering system to which a supplemental electric field is applied. The resonator structure may also include a semiconductor material having a positive temperature coefficient of resonance such that the resonator has a temperature coefficient of resonance approaching 0 ppm/.degree.C.
Rouze, Ned C.; Palmeri, Mark L.; Nightingale, Kathryn R.
2015-01-01
Recent measurements of shear wave propagation in viscoelastic materials have been analyzed by constructing the two-dimensional Fourier transform (2D-FT) of the spatial-temporal shear wave signal and using an analysis procedure derived under the assumption the wave is described as a plane wave, or as the asymptotic form of a wave expanding radially from a cylindrically symmetric source. This study presents an exact, analytic expression for the 2D-FT description of shear wave propagation in viscoelastic materials following asymmetric Gaussian excitations and uses this expression to evaluate the bias in 2D-FT measurements obtained using the plane or cylindrical wave assumptions. A wide range of biases are observed depending on specific values of frequency, aspect ratio R of the source asymmetry, and material properties. These biases can be reduced significantly by weighting the shear wave signal in the spatial domain to correct for the geometric spreading of the shear wavefront using a factor of xp. The optimal weighting power p is found to be near the theoretical value of 0.5 for the case of a cylindrical source with R = 1, and decreases for asymmetric sources with R > 1. PMID:26328717
Nonlinear Alfvén wave dynamics at a 2D magnetic null point: ponderomotive force
NASA Astrophysics Data System (ADS)
Thurgood, J. O.; McLaughlin, J. A.
2013-07-01
Context. In the linear, β = 0 MHD regime, the transient properties of magnetohydrodynamic (MHD) waves in the vicinity of 2D null points are well known. The waves are decoupled and accumulate at predictable parts of the magnetic topology: fast waves accumulate at the null point; whereas Alfvén waves cannot cross the separatricies. However, in nonlinear MHD mode conversion can occur at regions of inhomogeneous Alfvén speed, suggesting that the decoupled nature of waves may not extend to the nonlinear regime. Aims: We investigate the behaviour of low-amplitude Alfvén waves about a 2D magnetic null point in nonlinear, β = 0 MHD. Methods: We numerically simulate the introduction of low-amplitude Alfvén waves into the vicinity of a magnetic null point using the nonlinear LARE2D code. Results: Unlike in the linear regime, we find that the Alfvén wave sustains cospatial daughter disturbances, manifest in the transverse and longitudinal fluid velocity, owing to the action of nonlinear magnetic pressure gradients (viz. the ponderomotive force). These disturbances are dependent on the Alfvén wave and do not interact with the medium to excite magnetoacoustic waves, although the transverse daughter becomes focused at the null point. Additionally, an independently propagating fast magnetoacoustic wave is generated during the early stages, which transports some of the initial Alfvén wave energy towards the null point. Subsequently, despite undergoing dispersion and phase-mixing due to gradients in the Alfvén-speed profile (∇cA ≠ 0) there is no further nonlinear generation of fast waves. Conclusions: We find that Alfvén waves at 2D cold null points behave largely as in the linear regime, however they sustain transverse and longitudinal disturbances - effects absent in the linear regime - due to nonlinear magnetic pressure gradients.
Spatial variations in Achilles tendon shear wave speed
DeWall, Ryan J.; Slane, Laura C.; Lee, Kenneth S.; Thelen, Darryl G.
2014-01-01
Supersonic shear imaging (SSI) is an ultrasound imaging modality that can provide insight into tissue mechanics by measuring shear wave propagation speed, a property that depends on tissue elasticity. SSI has previously been used to characterize the increase in Achilles tendon shear wave speed that occurs with loading, an effect attributable to the strain-stiffening behavior of the tissue. However, little is known about how shear wave speed varies spatially, which is important, given the anatomical variation that occurs between the calcaneus insertion and the gastrocnemius musculotendon junction. The purpose of this study was to investigate spatial variations in shear wave speed along medial and lateral paths of the Achilles tendon for three different ankle postures: resting ankle angle (R, i.e. neutral), plantarflexed (P; R − 15 deg), and dorsiflexed (D; R + 15 deg). We observed significant spatial and posture variations in tendon shear wave speed in ten healthy young adults. Shear wave speeds in the Achilles free tendon averaged 12 ± 1.2 m/s in a resting position, but decreased to 7.2 ± 1.8 m/s with passive plantarflexion. Distal tendon shear wave speeds often reached the maximum tracking limit (16.3 m/s) of the system when the ankle was in the passively dorsiflexed posture (+15 deg from R). At a fixed posture, shear wave speeds decreased significantly from the free tendon to the gastrocnemius musculotendon junction, with slightly higher speeds measured on the medial side than on the lateral side. Shear wave speeds were only weakly correlated with the thickness and depth of the tendon, suggesting that the distal-to-proximal variations may reflect greater compliance in the aponeurosis relative to the free tendon. The results highlight the importance of considering both limb posture and transducer positioning when using SSI for biomechanical and clinical assessments of the Achilles tendon. PMID:24933528
Large amplitude compression and shear wave propagation in an elastomer
NASA Astrophysics Data System (ADS)
Gupta, Y. M.; Murri, W. J.; Henley, D.
1982-04-01
Experimental techniques have been developed to measure the high strain-rate compression and shear response of Solithane 113. Compression and shear wave profiles have been measured in specimens compressed to 20% (compressive stresses ˜1.2 GPa). The compressive profiles are nearly steady and the compressive stress-strain response is typical of a compliant material. The shear wave profiles are dispersive and show attenuation with propagation. Analyses of these wave profiles will be presented. Shear moduli vary from 0.35 GPa to 0.8 GPa for the compression range examined to date. These values are within a factor of two of the static shear moduli in the glassy state. The data described here have been used to calculate the high strain rate compressive and shear stress-strain curves for Solithane 113.
Triad resonance between gravity and vorticity waves in vertical shear
NASA Astrophysics Data System (ADS)
Drivas, Theodore D.; Wunsch, Scott
2016-07-01
Weakly nonlinear theory is used to explore the effect of vertical shear on surface gravity waves in three dimensions. An idealized piecewise-linear shear profile motivated by wind-driven profiles and ambient currents in the ocean is used. It is shown that shear may mediate weakly nonlinear resonant triad interactions between gravity and vorticity waves. The triad results in energy exchange between gravity waves of comparable wavelengths propagating in different directions. For realistic ocean shears, shear-mediated energy exchange may occur on timescales of minutes for shorter wavelengths, but slows as the wavelength increases. Hence this triad mechanism may contribute to the larger angular spreading (relative to wind direction) for shorter wind-waves observed in the oceans.
Shear waves in vegetal tissues at ultrasonic frequencies
NASA Astrophysics Data System (ADS)
Fariñas, M. D.; Sancho-Knapik, D.; Peguero-Pina, J. J.; Gil-Pelegrín, E.; Gómez Álvarez-Arenas, T. E.
2013-03-01
Shear waves are investigated in leaves of two plant species using air-coupled ultrasound. Magnitude and phase spectra of the transmission coefficient around the first two orders of the thickness resonances (normal and oblique incidence) have been measured. A bilayer acoustic model for plant leaves (comprising the palisade parenchyma and the spongy mesophyll) is proposed to extract, from measured spectra, properties of these tissues like: velocity and attenuation of longitudinal and shear waves and hence Young modulus, rigidity modulus, and Poisson's ratio. Elastic moduli values are typical of cellular solids and both, shear and longitudinal waves exhibit classical viscoelastic losses. Influence of leaf water content is also analyzed.
Excitation of fundamental shear horizontal wave by using face-shear (d36) piezoelectric ceramics
NASA Astrophysics Data System (ADS)
Miao, Hongchen; Dong, Shuxiang; Li, Faxin
2016-05-01
The fundamental shear horizontal (SH0) wave in plate-like structures is extremely useful for non-destructive testing (NDT) and structural health monitoring (SHM) as it is non-dispersive. However, currently, the SH0 wave is usually excited by electromagnetic acoustic transducers (EMAT) whose energy conversion efficiency is fairly low. The face-shear ( d 36 ) mode piezoelectrics is more promising for SH0 wave excitation, but this mode cannot appear in conventional piezoelectric ceramics. Recently, by modifying the symmetry of poled PbZr1-xTixO3 (PZT) ceramics via ferroelastic domain engineering, we realized the face-shear d 36 mode in both soft and hard PZT ceramics. In this work, we further improved the face-shear properties of PZT-4 and PZT-5H ceramics via lateral compression under elevated temperature. It was found that when bonded on a 1 mm-thick aluminum plate, the d 36 type PZT-4 exhibited better face-shear performance than PZT-5H. We then successfully excite SH0 wave in the aluminum plate using a face-shear PZT-4 square patch and receive the wave using a face-shear 0.72[Pb(Mg1/3Nb2/3)O3]-0.28[PbTiO3] (PMN-PT) patch. The frequency response and directionality of the excited SH0 wave were also investigated. The SH0 wave can be dominated over the Lamb waves (S0 and A0 waves) from 160 kHz to 280 kHz. The wave amplitude reaches its maxima along the two main directions (0° and 90°). The amplitude can keep over 80% of the maxima when the deviate angle is less than 30°, while it vanishes quickly at the 45° direction. The excited SH0 wave using piezoelectric ceramics could be very promising in the fields of NDT and SHM.
Shear viscosity of quasi-2D dipolar Bose-Fermi mixtures with long-range 1/r interactions
NASA Astrophysics Data System (ADS)
Darsheshdar, E.; Yavari, H.; Moniri, S. M.
2016-05-01
Low-temperature shear viscosity of a spin polarized two-component quasi-2D dipolar Fermi gas with long-range 1/ r interaction in the Bose-Einstein condensation (BEC) limit, where the system can be considered as dimers and the unpaired fermions, is calculated by means of the Kubo formalism. By taking into account the dimer-atom, dimer-dimer, and atom-atom interactions in the self-energies the viscous relaxation time (τ_{η}= (τ_{DA}^{-1}+τ_{DD}^{-1}+ τ_{AA}^{-1})^{-1}) is determined. Since the relaxation rates due to these interactions τ_{DA}^{-1} , τ_{DD}^{ -1} and τ_{AA}^{-1} varies, respectively, as T , T2 , and T in the low-temperature limit T→0 , the dimer-atom and atom-atom interactions play the dominant role to the shear viscosity and the shear viscosity varies as T^{-1} . For small polarization the effect of dimer-dimer interaction is important (τ_{DA},τ_{AA}≫τ_{DD}) , and the shear viscosity changes as the standard T^{-2} behviour. In this case, the temperature behavior of the dimer relaxation rate unaffected by 1/ r interaction and the contact, dipole-dipole, and 1/ r interactions play the same role in the temperature dependence of the shear viscosity. Our results have important consequences for developing experiments and theoretical researches on the transport properties of ultracold gases with repulsive or attractive long range 1/ r interaction.
Flipping and scooping of curved 2D rigid fibers in simple shear: The Jeffery equations
NASA Astrophysics Data System (ADS)
Crowdy, Darren
2016-05-01
The dynamical system governing the motion of a curved rigid two-dimensional circular-arc fiber in simple shear is derived in analytical form. This is achieved by finding the solution for the associated low-Reynolds-number flow around such a fiber using the methods of complex analysis. Solutions of the dynamical system display the "flipping" and "scooping" recently observed in computational studies of three-dimensional fibers using linked rigid rod and bead-shell models [J. Wang et al., "Flipping, scooping, and spinning: Drift of rigid curved nonchiral fibers in simple shear flows," Phys. Fluids 24, 123304 (2012)]. To complete the Jeffery-type equations for a curved fiber in a linear flow field we also derive its evolution equations in an extensional flow. It is expected that the equations derived here also govern the motion of slender, curved, three-dimensional rigid fibers when they evolve purely in the plane of shear or strain.
NASA Astrophysics Data System (ADS)
Köhn, D.; Wilken, D.; Rabbel, W.
2012-04-01
The FINO3 project is aimed at the construction of an offshore research platform in the north-sea, hosting research projects dealing with offshore wind energy topics. As part of FINO3 our sub-project deals with the development of new seismic acquisition and inversion concepts for offshore-building foundation soil analysis. We are focussed on the determination of seismic parameters and structural information of the building plot of the platform. Possible changes of the shear modulus of the sediments in the vicinity of the FINO3 monopile due to mechanic loads on the platform are estimated by a tube-waveform tomography. The tube-waves are excited by a hammer blow at the internal wall of the FINO3 monopile above the water line. The tube-waves are propagating through the water column and the sediments and are measured in situ by hydrophones at the external wall of the monopile. Homogenous long wavelength starting models for the waveform tomography are estimated using simple 2D finite difference models. Possible shear-wave velocity starting models range from 150-300 m/s. The resolution of the tube-waveform tomography is estimated by simple chequerboard and random media models. Additionally first results of the data application in the vicinity of the monopile are presented.
Effect of shear on failure waves in soda lime glass
Clifton, R. J.; Mello, M.; Brar, N. S.
1998-07-10
By means of in-material stress gauges, failure waves in shock-compressed soda lime glass have been shown to be distinguished by a marked reduction in shear stress. To explore further the relation between failure waves and shearing resistance, a series of pressure-shear impact experiments have been performed involving the impact of a glass plate by a steel flyer plate and vice versa. The latter configuration is designed to allow direct measurements of the shearing resistance of the failed material. In both configurations, the normal and transverse motion of the free surface of the target is monitored using laser interferometry. The transverse velocity-time profiles show a pronounced loss in shearing resistance of the glass at impact velocities above the threshold for failure waves to occur.
On acoustic wave generation in uniform shear flow
NASA Astrophysics Data System (ADS)
Gogoberidze, G.
2016-07-01
The linear dynamics of acoustic waves and vortices in uniform shear flow is studied. For flows with very low shear rates, the dynamics of perturbations is adiabatic and can be described by the WKB approximation. However, for flows with moderate and high shear rates the WKB approximation is not appropriate, and alternative analysis shows that two important phenomena occur: acoustic wave over-reflection and wave generation by vortices. The later phenomenon is a known linear mechanisms for sound generation in shear flows, a mechanism that is related to the continuous spectrum that arises in linear shear flow dynamics. A detailed analytical study of these phenomena is performed and the main quantitative and qualitative characteristics of the radiated acoustic field are obtained and analyzed.
Shear wave speed dispersion and attenuation in granular marine sediments.
Kimura, Masao
2013-07-01
The reported compressional wave speed dispersion and attenuation could be explained by a modified gap stiffness model incorporated into the Biot model (the BIMGS model). In contrast, shear wave speed dispersion and attenuation have not been investigated in detail. No measurements of shear wave speed dispersion have been reported, and only Brunson's data provide the frequency characteristics of shear wave attenuation. In this study, Brunson's attenuation measurements are compared to predictions using the Biot-Stoll model and the BIMGS model. It is shown that the BIMGS model accurately predicts the frequency dependence of shear wave attenuation. Then, the shear wave speed dispersion and attenuation in water-saturated silica sand are measured in the frequency range of 4-20 kHz. The vertical stress applied to the sample is 17.6 kPa. The temperature of the sample is set to be 5 °C, 20 °C, and 35 °C in order to change the relaxation frequency in the BIMGS model. The measured results are compared with those calculated using the Biot-Stoll model and the BIMGS model. It is shown that the shear wave speed dispersion and attenuation are predicted accurately by using the BIMGS model. PMID:23862793
Convertion Shear Wave Velocity to Standard Penetration Resistance
NASA Astrophysics Data System (ADS)
Madun, A.; Tajuddin, S. A. A.; Abdullah, M. E.; Abidin, M. H. Z.; Sani, S.; Siang, A. J. L. M.; Yusof, M. F.
2016-07-01
Multichannel Analysis Surface Wave (MASW) measurement is one of the geophysics exploration techniques to determine the soil profile based on shear wave velocity. Meanwhile, borehole intrusive technique identifies the changes of soil layer based on soil penetration resistance, i.e. standard penetration test-number of blows (SPT-N). Researchers across the world introduced many empirical conversions of standard penetration test blow number of borehole data to shear wave velocity or vice versa. This is because geophysics test is a non-destructive and relatively fast assessment, and thus should be promoted to compliment the site investigation work. These empirical conversions of shear wave velocity to SPT-N blow can be utilised, and thus suitable geotechnical parameters for design purposes can be achieved. This study has demonstrated the conversion between MASW and SPT-N value. The study was conducted at the university campus and Sejagung Sri Medan. The MASW seismic profiles at the University campus test site and Sejagung were at a depth of 21 m and 13 m, respectively. The shear wave velocities were also calculated empirically using SPT-N value, and thus both calculated and measured shear wave velocities were compared. It is essential to note that the MASW test and empirical conversion always underestimate the actual shear wave velocity of hard layer or rock due to the effect of soil properties on the upper layer.
2-D Path Corrections for Local and Regional Coda Waves: A Test of Transportability
Mayeda, K M; Malagnini, L; Phillips, W S; Walter, W R; Dreger, D S; Morasca, P
2005-07-13
Reliable estimates of the seismic source spectrum are necessary for accurate magnitude, yield, and energy estimation. In particular, how seismic radiated energy scales with increasing earthquake size has been the focus of recent debate within the community and has direct implications on earthquake source physics studies as well as hazard mitigation. The 1-D coda methodology of Mayeda et al. [2003] has provided the lowest variance estimate of the source spectrum when compared against traditional approaches that use direct S-waves, thus making it ideal for networks that have sparse station distribution. The 1-D coda methodology has been mostly confined to regions of approximately uniform complexity. For larger, more geophysically complicated regions, 2-D path corrections may be required. We will compare performance of 1-D versus 2-D path corrections in a variety of regions. First, the complicated tectonics of the northern California region coupled with high quality broadband seismic data provides for an ideal ''apples-to-apples'' test of 1-D and 2-D path assumptions on direct waves and their coda. Next, we will compare results for the Italian Alps using high frequency data from the University of Genoa. For Northern California, we used the same station and event distribution and compared 1-D and 2-D path corrections and observed the following results: (1) 1-D coda results reduced the amplitude variance relative to direct S-waves by roughly a factor of 8 (800%); (2) Applying a 2-D correction to the coda resulted in up to 40% variance reduction from the 1-D coda results; (3) 2-D direct S-wave results, though better than 1-D direct waves, were significantly worse than the 1-D coda. We found that coda-based moment-rate source spectra derived from the 2-D approach were essentially identical to those from the 1-D approach for frequencies less than {approx}0.7-Hz, however for the high frequencies (0.7 {le} f {le} 8.0-Hz), the 2-D approach resulted in inter-station scatter
Shear Wave Propagation Across Filled Joints with the Effect of Interfacial Shear Strength
NASA Astrophysics Data System (ADS)
Li, J. C.; Liu, T. T.; Li, H. B.; Liu, Y. Q.; Liu, B.; Xia, X.
2015-07-01
The thin-layer interface model for filled joints is extended to analyze shear wave propagation across filled rock joints when the interfacial shear strength between the filling material and the rocks is taken into account. During the wave propagation process, the two sides of the filled joint are welded with the adjacent rocks first and slide on each other when the shear stress on the joint is greater than the interfacial shear strength. By back analysis, the relation between the shear stress and the relative tangential deformation of the filled joints is obtained from the present approach, which is shown as a cycle parallelogram. Comparison between the present approach and the existing method based on the zero-thickness interface model indicates that the present approach is efficient to analyze shear wave propagation across rock joints with slippery behavior. The calculation results show that the slippery behavior of joints is related to the interfacial failure. In addition, the interaction between the shear stress wave and the two sides of the filling joint influences not only the wave propagation process but also the dynamic response of the filled joint.
NASA Astrophysics Data System (ADS)
Luo, Xin; Lu, Xin; Cong, Chunxiao; Yu, Ting; Xiong, Qihua; Ying Quek, Su
2015-10-01
2D layered materials have recently attracted tremendous interest due to their fascinating properties and potential applications. The interlayer interactions are much weaker than the intralayer bonds, allowing the as-synthesized materials to exhibit different stacking sequences, leading to different physical properties. Here, we show that regardless of the space group of the 2D materials, the Raman frequencies of the interlayer shear modes observed under the typical configuration blue shift for AB stacked materials, and red shift for ABC stacked materials, as the number of layers increases. Our predictions are made using an intuitive bond polarizability model which shows that stacking sequence plays a key role in determining which interlayer shear modes lead to the largest change in polarizability (Raman intensity); the modes with the largest Raman intensity determining the frequency trends. We present direct evidence for these conclusions by studying the Raman modes in few layer graphene, MoS2, MoSe2, WSe2 and Bi2Se3, using both first principles calculations and Raman spectroscopy. This study sheds light on the influence of stacking sequence on the Raman intensities of intrinsic interlayer modes in 2D layered materials in general, and leads to a practical way of identifying the stacking sequence in these materials.
Three Dimensional Shear Wave Elastographic Reconstruction of Ablations*
Ingle, Atul; Varghese, Tomy
2014-01-01
This paper presents an algorithm for three dimensional (3D) reconstruction of tumor ablations using ultrasound electrode vibration elastography. Shear wave velocity, which is used as a surrogate for tissue stiffness, is estimated by perturbing the ablation needle and tracking frame-to-frame displacements using radiofrequency ultrasound echo data. This process is repeated over many imaging planes that share a common axis of intersection collinear with needle. A 3D volume is reconstructed by solving an optimization problem which smoothly approximates shear wave velocities on a stack of transverse planes. The mean shear wave velocity estimates obtained in the phantom experiments are within 20% of those measured using a commercial shear wave imaging system. PMID:25570587
NASA Astrophysics Data System (ADS)
Zvietcovich, Fernando; Yao, Jianing; Chu, Ying-Ju; Meemon, Panomsak; Rolland, Jannick P.; Parker, Kevin J.
2016-03-01
Optical Coherence Elastography (OCE) is a widely investigated noninvasive technique for estimating the mechanical properties of tissue. In particular, vibrational OCE methods aim to estimate the shear wave velocity generated by an external stimulus in order to calculate the elastic modulus of tissue. In this study, we compare the performance of five acquisition and processing techniques for estimating the shear wave speed in simulations and experiments using tissue-mimicking phantoms. Accuracy, contrast-to-noise ratio, and resolution are measured for all cases. The first two techniques make the use of one piezoelectric actuator for generating a continuous shear wave propagation (SWP) and a tone-burst propagation (TBP) of 400 Hz over the gelatin phantom. The other techniques make use of one additional actuator located on the opposite side of the region of interest in order to create an interference pattern. When both actuators have the same frequency, a standing wave (SW) pattern is generated. Otherwise, when there is a frequency difference df between both actuators, a crawling wave (CrW) pattern is generated and propagates with less speed than a shear wave, which makes it suitable for being detected by the 2D cross-sectional OCE imaging. If df is not small compared to the operational frequency, the CrW travels faster and a sampled version of it (SCrW) is acquired by the system. Preliminary results suggest that TBP (error < 4.1%) and SWP (error < 6%) techniques are more accurate when compared to mechanical measurement test results.
Compressive and Shear Wave Velocity Profiles using Seismic Refraction Technique
NASA Astrophysics Data System (ADS)
Aziman, M.; Hazreek, Z. A. M.; Azhar, A. T. S.; Haimi, D. S.
2016-04-01
Seismic refraction measurement is one of the geophysics exploration techniques to determine soil profile. Meanwhile, the borehole technique is an established way to identify the changes of soil layer based on number of blows penetrating the soil. Both techniques are commonly adopted for subsurface investigation. The seismic refraction test is a non-destructive and relatively fast assessment compared to borehole technique. The soil velocities of compressive wave and shear wave derived from the seismic refraction measurements can be directly utilised to calculate soil parameters such as soil modulus and Poisson’s ratio. This study investigates the seismic refraction techniques to obtain compressive and shear wave velocity profile. Using the vertical and horizontal geophones as well as vertical and horizontal strike directions of the transient seismic source, the propagation of compressive wave and shear wave can be examined, respectively. The study was conducted at Sejagung Sri Medan. The seismic velocity profile was obtained at a depth of 20 m. The velocity of the shear wave is about half of the velocity of the compression wave. The soil profiles of compressive and shear wave velocities were verified using the borehole data and showed good agreement with the borehole data.
Surface wave phase velocities from 2-D surface wave tomography studies in the Anatolian plate
NASA Astrophysics Data System (ADS)
Arif Kutlu, Yusuf; Erduran, Murat; Çakır, Özcan; Vinnik, Lev; Kosarev, Grigoriy; Oreshin, Sergey
2014-05-01
We study the Rayleigh and Love surface wave fundamental mode propagation beneath the Anatolian plate. To examine the inter-station phase velocities a two-station method is used along with the Multiple Filter Technique (MFT) in the Computer Programs in Seismology (Herrmann and Ammon, 2004). The near-station waveform is deconvolved from the far-station waveform removing the propagation effects between the source and the station. This method requires that the near and far stations are aligned with the epicentre on a great circle path. The azimuthal difference of the earthquake to the two-stations and the azimuthal difference between the earthquake and the station are restricted to be smaller than 5o. We selected 3378 teleseismic events (Mw >= 5.7) recorded by 394 broadband local stations with high signal-to-noise ratio within the years 1999-2013. Corrected for the instrument response suitable seismogram pairs are analyzed with the two-station method yielding a collection of phase velocity curves in various period ranges (mainly in the range 25-185 sec). Diffraction from lateral heterogeneities, multipathing, interference of Rayleigh and Love waves can alter the dispersion measurements. In order to obtain quality measurements, we select only smooth portions of the phase velocity curves, remove outliers and average over many measurements. We discard these average phase velocity curves suspected of suffering from phase wrapping errors by comparing them with a reference Earth model (IASP91 by Kennett and Engdahl, 1991). The outlined analysis procedure yields 3035 Rayleigh and 1637 Love individual phase velocity curves. To obtain Rayleigh and Love wave travel times for a given region we performed 2-D tomographic inversion for which the Fast Marching Surface Tomography (FMST) code developed by N. Rawlinson at the Australian National University was utilized. This software package is based on the multistage fast marching method by Rawlinson and Sambridge (2004a, 2004b). The
A pitfall in shallow shear-wave refraction surveying
Xia, J.; Miller, R.D.; Park, C.B.; Wightman, E.; Nigbor, R.
2002-01-01
The shallow shear-wave refraction method works successfully in an area with a series of horizontal layers. However, complex near-surface geology may not fit into the assumption of a series of horizontal layers. That a plane SH-wave undergoes wave-type conversion along an interface in an area of nonhorizontal layers is theoretically inevitable. One real example shows that the shallow shear-wave refraction method provides velocities of a converted wave rather than an SH- wave. Moreover, it is impossible to identify the converted wave by refraction data itself. As most geophysical engineering firms have limited resources, an additional P-wave refraction survey is necessary to verify if velocities calculated from a shear-wave refraction survey are velocities of converted waves. The alternative at this time may be the surface wave method, which can provide reliable S-wave velocities, even in an area of velocity inversion (a higher velocity layer underlain by a lower velocity layer). ?? 2002 Elsevier Science B.V. All rights reserved.
Galerkin Spectral Method for the 2D Solitary Waves of Boussinesq Paradigm Equation
Christou, M. A.; Christov, C. I.
2009-10-29
We consider the 2D stationary propagating solitary waves of the so-called Boussinesq Paradigm equation. The fourth- order elliptic boundary value problem on infinite interval is solved by a Galerkin spectral method. An iterative procedure based on artificial time ('false transients') and operator splitting is used. Results are obtained for the shapes of the solitary waves for different values of the dispersion parameters for both subcritical and supercritical phase speeds.
On System-Dependent Sources of Uncertainty and Bias in Ultrasonic Quantitative Shear-Wave Imaging.
Deng, Yufeng; Rouze, Ned C; Palmeri, Mark L; Nightingale, Kathryn R
2016-03-01
Ultrasonic quantitative shear-wave imaging methods have been developed over the last decade to estimate tissue elasticity by measuring the speed of propagating shear waves following acoustic radiation force excitation. This work discusses eight sources of uncertainty and bias arising from ultrasound system-dependent parameters in ultrasound shear-wave speed (SWS) measurements. Each of the eight sources of error is discussed in the context of a linear, isotropic, elastic, homogeneous medium, combining previously reported analyses with Field II simulations, full-wave 2-D acoustic propagation simulations, and experimental studies. Errors arising from both spatial and temporal sources lead to errors in SWS measurements. Arrival time estimation noise, speckle bias, hardware fluctuations, and phase aberration cause uncertainties (variance) in SWS measurements, while pulse repetition frequency (PRF) and beamforming errors, as well as coupling medium sound speed mismatch, cause biases in SWS measurements (accuracy errors). Calibration of the sources of bias is an important step in the development of shear-wave imaging systems. In a well-calibrated system, where the sources of bias are minimized, and averaging over a region of interest (ROI) is employed to reduce the sources of uncertainty, an SWS error can be expected. PMID:26886980
Analysis of vegetation effect on waves using a vertical 2-D RANS model
Technology Transfer Automated Retrieval System (TEKTRAN)
A vertical two-dimensional (2-D) model has been applied in the simulation of wave propagation through vegetated water bodies. The model is based on an existing model SOLA-VOF which solves the Reynolds-Averaged Navier-Stokes (RANS) equations with the finite difference method on a staggered rectangula...
NASA Astrophysics Data System (ADS)
Hong, T.-K.
2009-04-01
Understanding the shear-wave excitation mechanism is a key issue for effective seismic monitoring of underground nuclear explosions (UNEs). We often observe strong shear waves from UNEs, which causes difficulty in prompt discrimination of nuclear explosions from natural earthquakes. Various mechanisms have been proposed to explain the shear-wave excitation from the UNEs. Consensus on dominant mechanism of shear-wave excitation has not been made. To constrain the shear-wave excitation mechanism, we examine the consistency in shear-wave radiation pattern using a source-array slowness-wavenumber (F-K) analysis, which allows us to check the time-invariant feature in the shear waves. We examine regional and teleseismic waveforms for the UNEs of the Balapan test site and Nevada test site along with the Indian and North Korean UNEs. We observe consistent radiation pattern in both regional and teleseismic shear waves. The observed radiation pattern suggests that the shear waves were not excited azimuthally-isotropic. Shear waves observed in teleseismic distances are far weak compared to those in regional distances, which implies that shear waves are excited stronger at high takeoff angles. Also, spectra of shear waves display significantly low overshoot feature that is different from those of P phases. The time-invariant anisotropic radiation pattern, strong excitation in high takeoff angle and low overshoot feature allow us to constrain the shear-wave excitation mechanism.
WAVE ACTION AND BOTTOM SHEAR STRESSES IN LAKE ERIE
For Lake Erie, the amplitudes and periods of wind-driven, surface gravity waves were calculated by means of the SMB hindcasting method. Bottom orbital velocities and bottom shear stresses were then calculated using linear wave theory and Kajiura's (1968) turbulent oscillating bou...
Measurement of Oblique Impact-generated Shear Waves
NASA Technical Reports Server (NTRS)
Dahl, J. M.; Schultz, P. H.
2001-01-01
Experimental strain measurements reveal that oblique impacts can generate shear waves with displacements as large as those in the P-wave. Large oblique impacts may thus be more efficient sources of surface disruption than vertical impacts. Additional information is contained in the original extended abstract.
Ismail, A.; Anderson, N.
2007-01-01
Shear-wave velocity (Vs) as a function of soil stiffness is an essential parameter in geotechnical characterization of the subsurface. In this study, multichannel analysis of surface wave (MASW) and downhole methods were used to map the shear-wave velocity-structure and depth to the bed-rock surface at a 125m ?? 125m geotechnical site in Missouri. The main objective was to assess the suitability of the site for constructing a large, heavy building. The acquired multichannel surface wave data were inverted to provide 1D shear-wave velocity profile corresponding to each shot gather. These 1D velocity profiles were interpolated and contoured to generate a suite of 2D shear-wave velocity sections. Integrating the shear-wave velocity data from the MASW method with the downhole velocity data and the available borehole lithologic information enabled us to map shear-wave velocity-structure to a depth on the order of 20m. The bedrock surface, which is dissected by a significant cut-and-fill valley, was imaged. The results suggest that the study site will require special consideration prior to construction. The results also demonstrate the successful use of MASW methods, when integrated with downhole velocity measurements and borehole lithologic information, in the characterization of the near surface at the geotechnical sites. ?? 2007 European Association of Geoscientists & Engineers.
Analytical solution of boundary integral equations for 2-D steady linear wave problems
NASA Astrophysics Data System (ADS)
Chuang, J. M.
2005-10-01
Based on the Fourier transform, the analytical solution of boundary integral equations formulated for the complex velocity of a 2-D steady linear surface flow is derived. It has been found that before the radiation condition is imposed, free waves appear both far upstream and downstream. In order to cancel the free waves in far upstream regions, the eigensolution of a specific eigenvalue, which satisfies the homogeneous boundary integral equation, is found and superposed to the analytical solution. An example, a submerged vortex, is used to demonstrate the derived analytical solution. Furthermore, an analytical approach to imposing the radiation condition in the numerical solution of boundary integral equations for 2-D steady linear wave problems is proposed.
Well-posedness and generalized plane waves simulations of a 2D mode conversion model
NASA Astrophysics Data System (ADS)
Imbert-Gérard, Lise-Marie
2015-12-01
Certain types of electro-magnetic waves propagating in a plasma can undergo a mode conversion process. In magnetic confinement fusion, this phenomenon is very useful to heat the plasma, since it permits to transfer the heat at or near the plasma center. This work focuses on a mathematical model of wave propagation around the mode conversion region, from both theoretical and numerical points of view. It aims at developing, for a well-posed equation, specific basis functions to study a wave mode conversion process. These basis functions, called generalized plane waves, are intrinsically based on variable coefficients. As such, they are particularly adapted to the mode conversion problem. The design of generalized plane waves for the proposed model is described in detail. Their implementation within a discontinuous Galerkin method then provides numerical simulations of the process. These first 2D simulations for this model agree with qualitative aspects studied in previous works.
Arnal, Bastien; Pernot, Mathieu; Tanter, Mickael
2011-08-01
The clinical applicability of high-intensity focused ultrasound (HIFU) for noninvasive therapy is currently hampered by the lack of robust and real-time monitoring of tissue damage during treatment. The goal of this study is to show that the estimation of local tissue elasticity from shear wave imaging (SWI) can lead to a precise mapping of the lesion. HIFU treatment and monitoring were respectively performed using a confocal setup consisting of a 2.5-MHz single element transducer focused at 34 mm on ex vivo samples and an 8-MHz ultrasound diagnostic probe. Ultrasound-based strain imaging was combined with shear wave imaging on the same device. The SWI sequences consisted of 2 successive shear waves induced at different lateral positions. Each wave was created with pushing beams of 100 μs at 3 depths. The shear wave propagation was acquired at 17,000 frames/s, from which the elasticity map was recovered. HIFU sonications were interleaved with fast imaging acquisitions, allowing a duty cycle of more than 90%. Thus, elasticity and strain mapping was achieved every 3 s, leading to real-time monitoring of the treatment. When thermal damage occurs, tissue stiffness was found to increase up to 4-fold and strain imaging showed strong shrinkages that blur the temperature information. We show that strain imaging elastograms are not easy to interpret for accurate lesion characterization, but SWI provides a quantitative mapping of the thermal lesion. Moreover, the concept of shear wave thermometry (SWT) developed in the companion paper allows mapping temperature with the same method. Combined SWT and shear wave imaging can map the lesion stiffening and temperature outside the lesion, which could be used to predict the eventual lesion growth by thermal dose calculation. Finally, SWI is shown to be robust to motion and reliable in vivo on sheep muscle. PMID:21859579
Numerical modeling of electromagnetic waves scattering from 2D coastal breaking sea waves
NASA Astrophysics Data System (ADS)
Khairi, Refzul; Coatanhay, Arnaud; Khenchaf, Ali; Scolan, Yves Marie
2013-11-01
The aim of this work is to model the interaction of L-band electromagnetic waves with coastal breaking sea waves. The breaking sea waves' profiles are generated using the desingularized technique and the electromagnetic waves scattering is computed using the high-order method of moments (HO-MoM) combined with non uniform rational basis spline (NURBS) geometry. Our study mainly focuses upon the electromagnetic waves behavior in the crest and the cavity of breaking sea waves. Contribution to the Topical Issue "Numelec 2012", Edited by Adel Razek.
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.
Lithology and shear-wave velocity in Memphis, Tennessee
Gomberg, J.; Waldron, B.; Schweig, E.; Hwang, H.; Webbers, A.; Van Arsdale, R.; Tucker, K.; Williams, R.; Street, R.; Mayne, P.; Stephenson, W.; Odum, J.; Cramer, C.; Updike, R.; Hutson, S.; Bradley, M.
2003-01-01
We have derived a new three-dimensional model of the lithologic structure beneath the city of Memphis, Tennessee, and examined its correlation with measured shear-wave velocity profiles. The correlation is sufficiently high that the better-constrained lithologic model may be used as a proxy for shear-wave velocities, which are required to calculate site-amplification for new seismic hazard maps for Memphis. The lithologic model and its uncertainties are derived from over 1200 newly compiled well and boring logs, some sampling to 500 m depth, and a moving-least-squares algorithm. Seventy-six new shear-wave velocity profiles have been measured and used for this study, most sampling to 30 m depth or less. All log and velocity observations are publicly available via new web sites.
Shear wave velocity structures of the Arabian Peninsula
NASA Astrophysics Data System (ADS)
Mokhtar, Talal A.; Al-Saeed, Mohammed M.
1994-02-01
The shear velocity structures of the different tectonic provinces of the Arabian Peninsula has been studied using surface wave data recorded by the RYD (Riyadh) station. The inversion of Rayleigh wave group velocities indicates that the Arabian shield can be modeled by two layers, each of which is 20 km thick with a shear velocity of 3.61 km/s in the upper crust and 3.88 km/s in the lower crust. The underlying upper mantle velocity is 4.61 km/s. Inversion of both Love and Rayleigh waves group velocities shows that the Arabian platform upper and lower crusts are comparable in their thicknesses to those of the shield, but with shear velocities of 3.4 and 4 km/s, respectively. The upper mantle velocity beneath the platform is 4.4 km/s and the average total thickness of the crust is 45 km.
Terrane-controlled crustal shear wave splitting in Taiwan
NASA Astrophysics Data System (ADS)
Okaya, David; Christensen, Nikolas I.; Ross, Zachary E.; Wu, Francis T.
2016-01-01
Taiwan is the result of arc-continent collision associated with the convergence of the Philippine Sea plate with the eastern Eurasian plate continental margin. The locus of deformation is found in eastern Taiwan in the form of mountain building (Central Range) with underlying thickened lithosphere. Rapid tectonic exhumation in the Central Range has uncovered low-to-high-grade metamorphic rocks marked by steep cleavage. We carried out a crustal seismic anisotropy study across Taiwan, producing a database of over 27,000 local earthquake shear wave splitting measurements. Additionally, we carried out rock physics measurements of metamorphic outcrop samples to quantify shear wave rock anisotropy. We produced a map of station-averaged splitting measurements across Taiwan. Patterns of fast shear wave directions correlate with tectonic terranes produced by plate convergence. Deformation-related mineral-preferred orientation in the metamorphic rocks produces a significant amount of the crustal anisotropy in the Taiwan collision zone.
Could linear hysteresis contribute to shear wave losses in tissues?
Parker, Kevin J
2015-04-01
For nearly 100 y in the study of cyclical motion in materials, a particular phenomenon called "linear hysteresis" or "ideal hysteretic damping" has been widely observed. More recently in the field of shear wave elastography, the basic mechanisms underlying shear wave losses in soft tissues are in question. Could linear hysteresis play a role? An underlying theoretical question must be answered: Is there a real and causal physical model that is capable of producing linear hysteresis over a band of shear wave frequencies used in diagnostic imaging schemes? One model that can approximately produce classic linear hysteresis behavior, by examining a generalized Maxwell model with a specific power law relaxation spectrum, is described here. This provides a theoretical plausibility for the phenomenon as a candidate for models of tissue behavior. PMID:25701527
Wave Transformation and Breaking on a Sheared Current
NASA Astrophysics Data System (ADS)
Zippel, S.; Thomson, J. M.; Rusch, C.
2014-12-01
Waves shoaling against tidal currents at river inlets have long been a hazard to navigation. We present measurements of waves, currents, and turbulence from SWIFT drifters at the Columbia River Mouth to diagnose wave transformation, breaking, and the resulting turbulence. In particular, down-looking velocity profiles, measured onboard the drifters, allow for evaluation of wave transformation on a vertically sheared current, for which theory exists but few in situ measurements are available. One consequence of wave transformation is steepening and breaking, which is identified using visual images, increased near surface turbulence, and gradients in wave energy flux. Vertical turbulent dissipation profiles measured during breaking are compared to existing scalings developed for deep and shallow water and expanded to the intermediate depth conditions common at the Columbia River Mouth. The analysis is intended to improve hydrodynamic models, especially two-way coupled wave-current models, and to aid navigation by better predicting dangerous wave conditions.
Crustal Shear Wave Anisotropy in the Taiwan Orogen
NASA Astrophysics Data System (ADS)
Rau, R.; Yang, C.
2002-12-01
Crustal shear wave anisotropy is analyzed in seismograms from local earthquakes (1.5 < M < 4.5) recorded at 75 permanent network stations in Taiwan during the period between 1991 and 2000. We investigate the origin of Taiwan crustal anisotropy by analyzing splitting in recorded shear waves, which can be characterized by a fast polarization direction and a time delay between fast and slow shear waves. Particle motion analysis and cross-correlation method are applied to estimate the splitting parameters of local shear waves. Clear evidence of shear wave splitting with split times of 0.02 to 0.2 s is found in about 75% of the stations studied. Stations in Coastal Plain, the foreland basin, show consistent alignment of the fast polarization directions that are parallel to the directions of local maximum horizontal compressive stress. Around the Chukou fault in the foothills region, fast shear-wave polarization directions of over 20 earthquakes recorded at station TWL are normal to the direction of local maximum horizontal compressive stress and parallel to the NE-SW trend of this east-dipping thrust fault. Except TWL, all the stations in foothills show large scatter in measured shear-wave polarizations. In the southern Central Range, two stations with a stable polarization direction of NNE-SSW, which is consistent with the local preferential mineral orientation, are found near the ChaoChou fault system. Two persistent polarization directions of fast shear wave are observed in a station (STY) near the boundary between the foothills and the southern Central Range: NW-SE direction for earthquakes located beneath the foothills and NE-SW direction for earthquakes coming under the Central Range. Our study indicates that the crustal anisotropy beneath the Taiwan orogen cannot be simply explained by the hypothesis of extensive dilatancy anisotropy (EDA), where parallel alignment of fluid-filled fractures produces the anisotropy. Other factors, such as intrinsic rock anisotropy
2D Traveling Wave Array Employing a Trapezoidal Dielectric Wedge for Beam Steering
NASA Technical Reports Server (NTRS)
Host, Nicholas K.; Chen, Chi-Chih; Volakis, John L.; Miranada, Felix A.
2014-01-01
This presentation addresses the progress made so far in the development of an antenna array with reconfigurable transmission line feeds connecting each element in series. In particular, 2D traveling wave array employing trapezoidal Dielectric Wedge for Beam Steering will be discussed. The presentation includes current status of the effort and suggested future work. The work is being done as part of the NASA Office of the Chief Technologist's Space Technology Research Fellowship (NSTRF).
The parametric decay of Alfven waves into shear Alfven waves and dust lower hybrid waves
Jamil, M.; Shah, H. A.; Zubia, K.; Zeba, I.; Uzma, Ch.; Salimullah, M.
2010-07-15
The parametric decay instability of Alfven wave into low-frequency electrostatic dust-lower-hybrid and electromagnetic shear Alfven waves has been investigated in detail in a dusty plasma in the presence of external/ambient uniform magnetic field. Magnetohydrodynamic fluid equations of plasmas have been employed to find the linear and nonlinear response of the plasma particles for this three-wave nonlinear coupling in a dusty magnetoplasma. Here, relatively high frequency electromagnetic Alfven wave has been taken as the pump wave. It couples with other two low-frequency internal possible modes of the dusty magnetoplasma, viz., the dust-lower-hybrid and shear Alfven waves. The nonlinear dispersion relation of the dust-lower-hybrid wave has been solved to obtain the growth rate of the parametric decay instability. The growth rate is maximum for small value of external magnetic field B{sub s}. It is noticed that the growth rate is proportional to the unperturbed electron number density n{sub oe}.
A nearly analytic exponential time difference method for solving 2D seismic wave equations
NASA Astrophysics Data System (ADS)
Zhang, Xiao; Yang, Dinghui; Song, Guojie
2014-02-01
In this paper, we propose a nearly analytic exponential time difference (NETD) method for solving the 2D acoustic and elastic wave equations. In this method, we use the nearly analytic discrete operator to approximate the high-order spatial differential operators and transform the seismic wave equations into semi-discrete ordinary differential equations (ODEs). Then, the converted ODE system is solved by the exponential time difference (ETD) method. We investigate the properties of NETD in detail, including the stability condition for 1-D and 2-D cases, the theoretical and relative errors, the numerical dispersion relation for the 2-D acoustic case, and the computational efficiency. In order to further validate the method, we apply it to simulating acoustic/elastic wave propagation in multilayer models which have strong contrasts and complex heterogeneous media, e.g., the SEG model and the Marmousi model. From our theoretical analyses and numerical results, the NETD can suppress numerical dispersion effectively by using the displacement and gradient to approximate the high-order spatial derivatives. In addition, because NETD is based on the structure of the Lie group method which preserves the quantitative properties of differential equations, it can achieve more accurate results than the classical methods.
Observations of wave shear stress on a steep beach
NASA Astrophysics Data System (ADS)
Wilson, G. W.; Hay, A. E.; Bowen, A. J.
2014-11-01
Observations are presented of the wave shear stress on a steeply sloping beach. Above the wave boundary layer (WBL), positive values of were observed and are attributed to a combination of both wave shoaling due to the large-scale bed slope, and dissipation due to wave breaking, in agreement with the wave theory of Zou et al. (2003). Within the WBL, observed vertical profiles of were also in good agreement with theory, in cases where the wave height was small. As wave heights increased, however, the WBL profile of generally did not agree with theory. Near-simultaneous rotary sonar observations of the bed suggest the disagreement with theory was due to the presence of orbital-scale ripples, which the present theory does not accommodate.
Wave Propagation in 2-D Granular Matrix and Dust Mitigation of Fabrics for Space Exploration Mission
NASA Technical Reports Server (NTRS)
Thanh, Phi Hung X.
2004-01-01
Wave Propagation study is essential to exploring the soil on Mars or Moon and Dust Mitigation is a necessity in terms of crew's health in exploration missions. The study of Dust Mitigation has a significant impact on the crew s health when astronauts track dust back into their living space after exploration trips. We are trying to use piezoelectric fiber to create waves and vibrations at certain critical frequencies and amplitudes so that we can shake the particles off from the astronaut s fabrics. By shaking off the dust and removing it, the astronauts no longer have to worry about breathing in small and possibly hazardous materials, when they are back in their living quarters. The Wave Propagation in 2-D Granular Matrix studies how the individual particles interact with each other when a pressure wave travels through the matrix. This experiment allows us to understand how wave propagates through soils and other materials. By knowing the details about the interactions of particles when they act as a medium for waves, we can better understand how wave propagates through soils and other materials. With this experiment, we can study how less gravity effects the wave propagation and hence device a way to study soils in space and on Moon or Mars. Some scientists treat the medium that waves travel through as a "black box", they did not pay much attention to how individual particles act as wave travels through them. With this data, I believe that we can use it to model ways to measure the properties of different materials such as density and composition. In order to study how the particles interact with each other, I have continued Juan Agui's experiment of the effects of impacts on a 2-D matrix. By controlling the inputs and measuring the outputs of the system, I will be able to study now the particles in that system interact with each other. I will also try to model this with the software called PFC2D in order to obtain theoretical data to compare with the experiment
Instability of plasma waves during relaxation of 2D turbulent flows
NASA Astrophysics Data System (ADS)
Kabantsev, A. A.; Drsicoll, C. F.
2015-11-01
We observe strong excitation of novel low-frequency z-dependent plasma waves (mθ = 0 ,kz = 1) , occurring during the nominally 2D relaxation of turbulent initial conditions (10 -100 interacting vortices) in strongly magnetized electron columns. This initial relaxation often results in ``2D vortex crystal'' states. Here we describe experiments showing the concomitant growth of ill-understood low-frequency plasma waves, probably due to ``leakage'' of 2D turbulent potential energy into z-dependent fluctuations. With plasma injection, the lowest regular Trivelpiece- Gould mode (mθ = 0 ,kz = 1) is observed at fTG (t) ~ 2 . 8 MHz and exponential decay time τTG ~ 1 msec. Also, we observe rapid exponential growth of a novel low-frequency mode with fLF (t) ~ 0 . 3 MHz, nominally also with mθ = 0 ,kz = 1 . In a few milliseconds (several tens of rotation times at B = 10kG), the LF-mode becomes highly nonlinear, developing up to a dozen temporal harmonics. When a LF-harmonic resonates with the decaying TG-mode, LF-mode energy is transferred into the TG-mode, and both modes remain at moderate amplitudes until the 2D turbulent relaxation abates (hundreds of rotation times). The ill-understood fLF is independent of B, even though the growth and duration times follow scale as B1 from the 2D flows. Supported by National Science Foundation Grant PHY-1414570, Department of Energy Grants DE-SC0008693.
Shear wave elastography using phase sensitive optical coherence tomography
NASA Astrophysics Data System (ADS)
Song, Shaozhen; Huang, Zhihong; Nguyen, Thu-Mai; Wong, Emily Y.; Arnal, Bastien; O'Donnell, Matthew; Wang, Ruikang K.
2014-03-01
Optical coherence tomography (OCT) provides high spatial resolution and sensitivity that are ideal for imaging the cornea and lens. Quantifying the biomechanical properties of these tissues could add clinically valuable information. Thus, we propose a dynamic elastography method combining OCT detection and a mechanical actuator to map the shear modulus of soft tissues. We used a piezoelectric actuator driven in the kHz range and we used phase-sensitive OCT (PhS-OCT) to track the resulting shear waves at an equivalent frame rate of 47 kHz. We mapped the shear wave speed of anesthetized mice cornea using monochromatic excitations. We found a significant difference between a group of knock-out (3.92 +/- 0.35 m/s, N=4) and wild-type mice (5.04 +/- 0.51 m/s, N=3). These preliminary results demonstrate the feasibility of using PhS-OCT to perform in vivo shear wave elastography of the cornea. We then implemented a shear pulse compression approach on ex vivo human cornea. For that purpose, frequency- modulated excitations were used and the resulting displacement field was digitally compressed in a short broadband pulse with a 7 dB gain in signal-to-noise ratio (SNR).
A Hammer-Impact, Aluminum, Shear-Wave Seismic Source
Haines, Seth S.
2007-01-01
Near-surface seismic surveys often employ hammer impacts to create seismic energy. Shear-wave surveys using horizontally polarized waves require horizontal hammer impacts against a rigid object (the source) that is coupled to the ground surface. I have designed, built, and tested a source made out of aluminum and equipped with spikes to improve coupling. The source is effective in a variety of settings, and it is relatively simple and inexpensive to build.
Scattering of elastic waves by a 2-D crack using the Indirect Boundary Element Method (IBEM)
NASA Astrophysics Data System (ADS)
Iturrarán-Viveros, Ursula; Vai, Rossana; Sánchez-Sesma, Francisco J.
2005-09-01
The scattering of elastic waves by cracks is an old problem and various ways to solve it have been proposed in the last decades. One approach is using dual integral equations, another useful and common formulation is the Boundary Element Method (BEM). With the last one, the boundary conditions of the crack lead to hyper-singularities and particular care should be taken to regularize and solve the resulting integral equations. In this work, instead, the Indirect Boundary Element Method (IBEM) is applied to study problems of zero-thickness 2-D cracks. The IBEM yields the Crack Opening Displacement (COD) which is used to evaluate the solution away from the crack. We use a multiregional approach which consists of splitting a boundary S into two identical boundaries S+ and S- chosen such that the cracks lie in the interface. The resulting integral equations are not hyper-singular and wave propagation within media that contain zero-thickness cracks can be rigorously solved. In order to validate the method, we deal with the scalar case, namely the scattering of antiplane SH waves by a 2-D crack. We compare results against a recently published analytic solution, obtaining an excellent agreement. This comparison gives us confidence to study cases where no analytic solutions exist. Some examples of incidence of P- or SV waves are depicted and the salient aspects of the method are also discussed.
Laser probe for measuring 2-D wave slope spectra of ocean capillary waves
NASA Technical Reports Server (NTRS)
Palm, C. S.; Anderson, R. C.; Reece, A. M.
1977-01-01
A laser-optical instrument for use in determining the two-dimensional wave-slope spectrum of ocean capillary waves is described. The instrument measures up to a 35-deg tip angle of the surface normal by measuring the position of a refracted laser beam directed vertically upward through a water surface. A telescope, a continuous two-dimensional Schottky barrier photodiode, and a pair of analog dividers render the signals independent of water height and insensitive to laser-beam intensity fluctuations. Calibration is performed entirely in the laboratory before field use. Sample records and wave-slope spectra are shown for one-dimensional wave-tank tests and for two-dimensional ocean tests. These are presented along with comparison spectra for calm and choppy water conditions. A mechanical wave follower was used to adjust the instrument position in the presence of large ocean swell and tides.
Density gradient effects on transverse shear driven lower hybrid waves
DuBois, Ami M.; Thomas, Edward; Amatucci, William E.; Ganguli, Gurudas
2014-06-15
Shear driven instabilities are commonly observed in the near-Earth space, particularly in boundary layer plasmas. When the shear scale length (L{sub E}) is much less than the ion gyro-radius (ρ{sub i}) but greater than the electron gyro-radius (ρ{sub e}), the electrons are magnetized in the shear layer, but the ions are effectively un-magnetized. The resulting shear driven instability, the electron-ion hybrid (EIH) instability, is investigated in a new interpenetrating plasma configuration in the Auburn Linear EXperiment for Instability Studies. In order to understand the dynamics of magnetospheric boundary layers, the EIH instability is studied in the presence of a density gradient located at the boundary layer between two plasmas. This paper reports on a recent experiment in which electrostatic lower hybrid waves are identified as the EIH instability, and the effect of a density gradient on the instability properties are investigated.
Shear Wave Elastography and Cervical Lymph Nodes: Predicting Malignancy.
Azizi, Ghobad; Keller, James M; Mayo, Michelle L; Piper, Kelé; Puett, David; Earp, Karly M; Malchoff, Carl D
2016-06-01
This prospective study evaluates the accuracy of virtual touch imaging quantification (VTIQ), a non-invasive shear wave elastography method for measuring cervical lymph nodes (LN) stiffness in differentiating benign from malignant LN. The study evaluated 270 LN in 236 patients with both conventional B-mode ultrasound and VTIQ shear wave elastography before fine-needle aspiration biopsy (FNAB). LN stiffness was measured as shear wave velocity (SWV) in m/s. Surgical resection was advised for FNAB results that were not clearly benign. Surgical pathology confirmed 54 malignant LN. The receiver operating curve (ROC) identified a single cut-off value of 2.93 m/s as the maximum SWV for predicting a malignant cervical LN. The sensitivity and specificity were 92.59% and 75.46%, respectively. Positive predictive value (PPV) was 48.54% and negative predictive value (NPV) was 97.60%. LN stiffness measured by VTIQ-generated shear wave elastography is an independent predictor of malignancy. PMID:26976785
Fast Acceleration of 2D Wave Propagation Simulations Using Modern Computational Accelerators
Wang, Wei; Xu, Lifan; Cavazos, John; Huang, Howie H.; Kay, Matthew
2014-01-01
Recent developments in modern computational accelerators like Graphics Processing Units (GPUs) and coprocessors provide great opportunities for making scientific applications run faster than ever before. However, efficient parallelization of scientific code using new programming tools like CUDA requires a high level of expertise that is not available to many scientists. This, plus the fact that parallelized code is usually not portable to different architectures, creates major challenges for exploiting the full capabilities of modern computational accelerators. In this work, we sought to overcome these challenges by studying how to achieve both automated parallelization using OpenACC and enhanced portability using OpenCL. We applied our parallelization schemes using GPUs as well as Intel Many Integrated Core (MIC) coprocessor to reduce the run time of wave propagation simulations. We used a well-established 2D cardiac action potential model as a specific case-study. To the best of our knowledge, we are the first to study auto-parallelization of 2D cardiac wave propagation simulations using OpenACC. Our results identify several approaches that provide substantial speedups. The OpenACC-generated GPU code achieved more than speedup above the sequential implementation and required the addition of only a few OpenACC pragmas to the code. An OpenCL implementation provided speedups on GPUs of at least faster than the sequential implementation and faster than a parallelized OpenMP implementation. An implementation of OpenMP on Intel MIC coprocessor provided speedups of with only a few code changes to the sequential implementation. We highlight that OpenACC provides an automatic, efficient, and portable approach to achieve parallelization of 2D cardiac wave simulations on GPUs. Our approach of using OpenACC, OpenCL, and OpenMP to parallelize this particular model on modern computational accelerators should be applicable to other computational models of wave propagation in
Cosmic shear from scalar-induced gravitational waves
Sarkar, Devdeep; Serra, Paolo; Cooray, Asantha; Ichiki, Kiyotomo; Baumann, Daniel
2008-05-15
Weak gravitational lensing by foreground density perturbations generates a gradient mode in the shear of background images. In contrast, cosmological tensor perturbations induce a nonzero curl mode associated with image rotations. In this note, we study the lensing signatures of both primordial gravitational waves from inflation and second-order gravitational waves generated from the observed spectrum of primordial density fluctuations. We derive the curl mode for galaxy lensing surveys at redshifts of 1-3 and for lensing of the cosmic microwave background at a redshift of 1100. We find that the curl mode angular power spectrum associated with secondary tensor modes for galaxy lensing surveys dominates over the corresponding signal generated by primary gravitational waves from inflation. However, both tensor contributions to the shear curl mode spectrum are below the projected noise levels of upcoming galaxy and cosmic microwave background lensing surveys and therefore are unlikely to be detectable.
NASA Astrophysics Data System (ADS)
Martowicz, A.; Ruzzene, M.; Staszewski, W. J.; Rimoli, J. J.; Uhl, T.
2014-03-01
The work deals with the reduction of numerical dispersion in simulations of wave propagation in solids. The phenomenon of numerical dispersion naturally results from time and spatial discretization present in a numerical model of mechanical continuum. Although discretization itself makes possible to model wave propagation in structures with complicated geometries and made of different materials, it inevitably causes simulation errors when improper time and length scales are chosen for the simulations domains. Therefore, by definition, any characteristic parameter for spatial and time resolution must create limitations on maximal wavenumber and frequency for a numerical model. It should be however noted that expected increase of the model quality and its functionality in terms of affordable wavenumbers, frequencies and speeds should not be achieved merely by denser mesh and reduced time integration step. The computational cost would be simply unacceptable. The authors present a nonlocal finite difference scheme with the coefficients calculated applying a Fourier series, which allows for considerable reduction of numerical dispersion. There are presented the results of analyses for 2D models, with isotropic and anisotropic materials, fulfilling the planar stress state. Reduced numerical dispersion is shown in the dispersion surfaces for longitudinal and shear waves propagating for different directions with respect to the mesh orientation and without dramatic increase of required number of nonlocal interactions. A case with the propagation of longitudinal wave in composite material is studied with given referential solution of the initial value problem for verification of the time-domain outcomes. The work gives a perspective of modeling of any type of real material dispersion according to measurements and with assumed accuracy.
NASA Astrophysics Data System (ADS)
Vanneste, Maarten; Madshus, Christian; Socco, Valentina L.; Maraschini, Margherita; Sparrevik, Per M.; Westerdahl, Harald; Duffaut, Kenneth; Skomedal, Eiliv; Bjørnarâ, Tore I.
2011-04-01
Pure shear wave data are only very rarely acquired for offshore site investigations and exploration. Here, we present details of a novel, seabed-coupled, shear wave vibrator and field data recorded by a densely populated, multicomponent ocean-bottom cable, to improve shallow soil characterization. The prototype shear wave vibrator uses vibroseis technology adopted for marine environments through its instalment on top of a suction anchor, assuring seabed coupling in combination with self-weight penetration. The prototype is depth rated to 1500 m water depth, and can be rotated while installed in the seabed. The philosophy is to acquire fully complementary seismic data to conventional P- and P-to-S-converted waves, in particular for 2-D profiling, VSP (vertical seismic profiling) or monitoring purposes, thereby exploiting advantages of shear waves over compressional waves for determining, for example, anisotropy, small-strain shear modulus and excess pore pressures/effective stress. The source was primarily designed for reservoir depths. However, significant energy is emitted as surface waves, which provide detailed geotechnical information through mapping of shear wave velocities in potentially high resolution of the upper soil units. To fully utilize pure shear wave content, a proper analysis of surface waves is paramount, due to the proximity of surface wave propagation speed with shear wave velocities. The experiment was carried out in the northern North Sea in 364 m water depth. Cable dragging was necessary to obtain close receiver spacing (2.5 m effective spacing), with total line length of 600 m. Frequency-waveform transforms reveal both Scholte and Love waves. Up to six surface wave modes are identified, that is, fundamental mode and several higher surface wave modes. The occurrence of these two different dispersive surface wave types with well-resolved higher modes allows for a unique analysis and inversion scheme for high-resolution mapping of physical
2-D Coda and Direct Wave Attenuation Tomography in Northern Italy
Morasca, P; Mayeda, K; Gok, R; Phillips, W S; Malagnini, L
2007-10-17
A 1-D coda method was proposed by Mayeda et al. (2003) in order to obtain stable seismic source moment-rate spectra using narrowband coda envelope measurements. That study took advantage of the averaging nature of coda waves to derive stable amplitude measurements taking into account all propagation, site, and Sto-coda transfer function effects. Recently this methodology was applied to micro earthquake data sets from three sub-regions of northern Italy (i.e., western Alps, northern Apennines and eastern Alps). Since the study regions were small, ranging between local-to-near-regional distances, the simple 1-D path assumptions used in the coda method worked very well. The lateral complexity of this region would suggest, however, that a 2-D path correction might provide even better results if the datasets were combined, especially when paths traverse larger distances and complicated regions. The structural heterogeneity of northern Italy makes the region ideal to test the extent to which coda variance can be reduced further by using a 2-D Q tomography technique. The approach we use has been developed by Phillips et al. (2005) and is an extension of previous amplitude ratio techniques to remove source effects from the inversion. The method requires some assumptions such as isotropic source radiation which is generally true for coda waves. Our results are compared against direct Swave inversions for 1/Q and results from both share very similar attenuation features that coincide with known geologic structures. We compare our results with those derived from direct waves as well as some recent results from northern California obtained by Mayeda et al. (2005) which tested the same tomographic methodology applied in this study to invert for 1/Q. We find that 2-D coda path corrections for this region significantly improve upon the 1-D corrections, in contrast to California where only a marginal improvement was observed. We attribute this difference to stronger lateral
Designing of sparse 2D arrays for Lamb wave imaging using coarray concept
NASA Astrophysics Data System (ADS)
Ambroziński, Łukasz; Stepinski, Tadeusz; Uhl, Tadeusz
2015-03-01
2D ultrasonic arrays have considerable application potential in Lamb wave based SHM systems, since they enable equivocal damage imaging and even in some cases wave-mode selection. Recently, it has been shown that the 2D arrays can be used in SHM applications in a synthetic focusing (SF) mode, which is much more effective than the classical phase array mode commonly used in NDT. The SF mode assumes a single element excitation of subsequent transmitters and off-line processing the acquired data. In the simplest implementation of the technique, only single multiplexed input and output channels are required, which results in significant hardware simplification. Application of the SF mode for 2D arrays creates additional degrees of freedom during the design of the array topology, which complicates the array design process, however, it enables sparse array designs with performance similar to that of the fully populated dense arrays. In this paper we present the coarray concept to facilitate synthesis process of an array's aperture used in the multistatic synthetic focusing approach in Lamb waves-based imaging systems. In the coherent imaging, performed in the transmit/receive mode, the sum coarray is a morphological convolution of the transmit/receive sub-arrays. It can be calculated as the set of sums of the individual sub-arrays' elements locations. The coarray framework will be presented here using a an example of a star-shaped array. The approach will be discussed in terms of beampatterns of the resulting imaging systems. Both simulated and experimental results will be included.
Designing of sparse 2D arrays for Lamb wave imaging using coarray concept
Ambroziński, Łukasz Stepinski, Tadeusz Uhl, Tadeusz
2015-03-31
2D ultrasonic arrays have considerable application potential in Lamb wave based SHM systems, since they enable equivocal damage imaging and even in some cases wave-mode selection. Recently, it has been shown that the 2D arrays can be used in SHM applications in a synthetic focusing (SF) mode, which is much more effective than the classical phase array mode commonly used in NDT. The SF mode assumes a single element excitation of subsequent transmitters and off-line processing the acquired data. In the simplest implementation of the technique, only single multiplexed input and output channels are required, which results in significant hardware simplification. Application of the SF mode for 2D arrays creates additional degrees of freedom during the design of the array topology, which complicates the array design process, however, it enables sparse array designs with performance similar to that of the fully populated dense arrays. In this paper we present the coarray concept to facilitate synthesis process of an array’s aperture used in the multistatic synthetic focusing approach in Lamb waves-based imaging systems. In the coherent imaging, performed in the transmit/receive mode, the sum coarray is a morphological convolution of the transmit/receive sub-arrays. It can be calculated as the set of sums of the individual sub-arrays’ elements locations. The coarray framework will be presented here using a an example of a star-shaped array. The approach will be discussed in terms of beampatterns of the resulting imaging systems. Both simulated and experimental results will be included.
NASA Astrophysics Data System (ADS)
Nguyen, Thu-Mai; Song, Shaozhen; Arnal, Bastien; Wong, Emily Y.; Shen, Tueng T.; Wang, Ruikang K.; O'Donnell, Matthew
2015-03-01
Tissue stiffness can be measured from the propagation speed of shear waves. Acoustic radiation force (ARF) can generate shear waves by focusing ultrasound in tissue for ~100 μs. Safety considerations and electronics abilities limit ultrasound pressures. We previously presented shear wave elastography combining ARF and phase-sensitive optical coherence tomography (PhS-OCT) [1]. Here, we use amplitude-modulated ARF to enhance shear wave signal-to-noise ratio (SNR) at low pressures. Experiments were performed on tissue-mimicking phantoms. ARF was applied using a single-element transducer, driven by a 7.5 MHz, 3-ms, sine wave modulated in amplitude by a linear-swept frequency (1 to 7 kHz). Pressures between 1 to 3 MPa were tested. Displacements were tracked using PhS-OCT and numerically compressed using pulse compression methods detailed in previous work [2]. SNR was compared to that of 200-μs bursts. Stiffness maps were reconstructed using time-of-flight computations. 200-μs bursts give barely detectable displacements at 1 MPa (3.7 dB SNR). Pulse compression gives 36.2 dB at 1.5 MPa. In all cases with detectable displacements, shear wave speeds were determined in 5%-gelatin and 10%-gelatin phantoms and compared to literature values. Applicability to ocular tissues (cornea, intraocular lens) is under investigation.
Ship waves on uniform shear current at finite depth: wave resistance and critical velocity
NASA Astrophysics Data System (ADS)
Li, Yan; Ellingsen, Simen Å.
2016-03-01
We present a comprehensive theory for linear gravity-driven ship waves in the presence of a shear current with uniform vorticity, including the effects of finite water depth. The wave resistance in the presence of shear current is calculated for the first time, containing in general a non-zero lateral component. While formally apparently a straightforward extension of existing deep water theory, the introduction of finite water depth is physically non-trivial, since the surface waves are now affected by a subtle interplay of the effects of the current and the sea bed. This becomes particularly pronounced when considering the phenomenon of critical velocity, the velocity at which transversely propagating waves become unable to keep up with the moving source. The phenomenon is well known for shallow water, and was recently shown to exist also in deep water in the presence of a shear current [Ellingsen, J.~Fluid Mech.\\ {\\bf 742} R2 (2014)]. We derive the exact criterion for criticality as a function of an intrinsic shear Froude number $S\\sqrt{b/g}$ ($S$ is uniform vorticity, $b$ size of source), the water depth, and the angle between the shear current and the ship's motion. Formulae for both the normal and lateral wave resistance force are derived, and we analyse its dependence on the source velocity (or Froude number $Fr$) for different amounts of shear and different directions of motion. The effect of the shear current is to increase wave resistance for upstream ship motion and decrease it for downstream motion. Also the value of $Fr$ at which $R$ is maximal is lowered for upstream and increased for downstream directions of ship motion. For oblique angles between ship motion and current there is a lateral wave resistance component which can amount to $10$-$20\\%$ of the normal wave resistance for side-on shear and $S\\sqrt{b/g}$ of order unity. (Continues...)
Parametric analysis of 2D guided-wave photonic band gap structures
NASA Astrophysics Data System (ADS)
Ciminelli, C.; Peluso, F.; Armenise, M. N.
2005-11-01
The parametric analysis of the electromagnetic properties of 2D guided wave photonic band gap structures is reported with the aim of providing a valid tool for the optimal design. The modelling approach is based on the Bloch-Floquet method. Different lattice configurations and geometrical parameters are considered. An optimum value for the ratio between the hole (or rod) radius and the lattice constant does exist and the calculation demonstrated that it is almost independent from the etching depth, only depending on the lattice type. The results are suitable for the design optimisation of photonic crystal reflectors to be used in integrated optical devices.
Parametric analysis of 2D guided-wave photonic band gap structures.
Ciminelli, C; Peluso, F; Armenise, M
2005-11-28
The parametric analysis of the electromagnetic properties of 2D guided wave photonic band gap structures is reported with the aim of providing a valid tool for the optimal design. The modelling approach is based on the Bloch-Floquet method. Different lattice configurations and geometrical parameters are considered. An optimum value for the ratio between the hole (or rod) radius and the lattice constant does exist and the calculation demonstrated that it is almost independent from the etching depth, only depending on the lattice type. The results are suitable for the design optimisation of photonic crystal reflectors to be used in integrated optical devices. PMID:19503180
Multipacting Simulation Study for 56 MHz Quarter Wave Resonator using 2D Code
Naik,D.; Ben-Zvi, I.
2009-01-02
A beam excited 56 MHz Radio Frequency (RF) Niobium Quarter Wave Resonator (QWR) has been proposed to enhance RHIC beam luminosity and bunching. Being a RF cavity, multipacting is expected; therefore an extensive study was carried out with the Multipac 2.1 2D simulation code. The study revealed that multipacting occurs in various bands up to peak surface electric field 50 kV/m and is concentrated mostly above the beam gap and on the outer conductor. To suppress multipacting, a ripple structure was introduced to the outer conductor and the phenomenon was successfully eliminated from the cavity.
Hα Moreton waves observed on December 06, 2006. A 2D case study
NASA Astrophysics Data System (ADS)
Francile, C.; Costa, A.; Luoni, M. L.; Elaskar, S.
2013-04-01
Context. We present high temporal resolution observations of a Moreton wave event detected with the Hα Solar Telescope for Argentina (HASTA) in the Hα line 656.3 nm, on December 6, 2006. Aims: The aim is to contribute to the discussion about the nature and triggering mechanisms of Moreton wave events. Methods: We describe the HASTA telescope capabilities and the observational techniques. We carried out a detailed analysis to determine the flare onset, the radiant point location, the kinematics of the disturbance and the activation time of two distant filaments. We used a 2D reconstruction of the HASTA and corresponding TRACE observations, together with conventional techniques, to analyze the probable origin of the phenomenon. Results: The kinematic parameters and the probable onset time of the Moreton wave event are determined. A small-scale ejectum and the winking of two remote filaments are analyzed to discuss their relation with the Moreton disturbance. Conclusions: The analysis of the Moreton wave event favors the hypothesis that the phenomenon can be described as the chromospheric imprint of a single fast coronal shock triggered from a single source in association with a coronal mass ejection. Its onset time is concurrent with a Lorentz force peak measured in the photosphere, as stated by other authors. However, the existence of multiple shock waves that were generated almost simultaneously cannot be discarded.
Bottom shear stress and pressure perturbations under an internal solitary wave
NASA Astrophysics Data System (ADS)
Rivera, Gustavo; Diamessis, Peter
2014-11-01
The bottom boundary layer (BBL) under a mode-1 internal solitary wave (ISW) of depression propagating against an oncoming model barotropic current is examined using 2-D direct numerical simulation based on a spectral multidomain penalty method model. Use of a postprocessing projection onto a modified set of divergence-free basis functions enables investigation of wave-based Reynolds numbers within the range [105 ,106 ] . At sufficiently high ISW amplitude, the BBL undergoes a global instability which produces intermittent vortex shedding from within the separation bubble in the lee of the wave. The interplay between the bottom shear stress field and pressure perturbations during vortex ejection events and the subsequent evolution of the vortices is the focus of this presentation. Implications for resuspension of bottom particulate matter are discussed in the context of specific sediment transport models. Support from the Cornell Sloan Diversity Fellowship program is gratefully acknowledged.
High-frequency shear-horizontal surface acoustic wave sensor
Branch, Darren W
2013-05-07
A Love wave sensor uses a single-phase unidirectional interdigital transducer (IDT) on a piezoelectric substrate for leaky surface acoustic wave generation. The IDT design minimizes propagation losses, bulk wave interferences, provides a highly linear phase response, and eliminates the need for impedance matching. As an example, a high frequency (.about.300-400 MHz) surface acoustic wave (SAW) transducer enables efficient excitation of shear-horizontal waves on 36.degree. Y-cut lithium tantalate (LTO) giving a highly linear phase response (2.8.degree. P-P). The sensor has the ability to detect at the pg/mm.sup.2 level and can perform multi-analyte detection in real-time. The sensor can be used for rapid autonomous detection of pathogenic microorganisms and bioagents by field deployable platforms.
High-frequency shear-horizontal surface acoustic wave sensor
Branch, Darren W
2014-03-11
A Love wave sensor uses a single-phase unidirectional interdigital transducer (IDT) on a piezoelectric substrate for leaky surface acoustic wave generation. The IDT design minimizes propagation losses, bulk wave interferences, provides a highly linear phase response, and eliminates the need for impedance matching. As an example, a high frequency (.about.300-400 MHz) surface acoustic wave (SAW) transducer enables efficient excitation of shear-horizontal waves on 36.degree. Y-cut lithium tantalate (LTO) giving a highly linear phase response (2.8.degree. P-P). The sensor has the ability to detect at the pg/mm.sup.2 level and can perform multi-analyte detection in real-time. The sensor can be used for rapid autonomous detection of pathogenic microorganisms and bioagents by field deployable platforms.
Shear Wave Structure in the Crust and Upper Mantle Beneath the Wyoming Craton
NASA Astrophysics Data System (ADS)
Dave, R.; Li, A.
2013-12-01
The Wyoming craton was pervasively deformed during the Laramide and has been encroached by the Yellowstone hotspot. It is therefore an ideal place to study the survival and evolution of cratonic lithosphere. We have analyzed Rayleigh wave data recorded by the USArray and applied the two-plane-wave tomography method to determine average and 2-D phase velocities in the Wyoming craton area. These phase velocities were inverted to construct 3-D shear wave structure of the crust and upper mantle. Vertical component seismograms from 82 events at 103 stations were filtered at 18 center frequencies with a narrow bandwidth of 10 mHz. The average phase velocity varies from 3.48 km/s at 20 s to 4.36 km/s at 166 s. Compared to the average global phase velocities, these values are lower at longer periods, indicating a generally slower upper mantle. 2-D variation of phase velocity shows significantly low velocities at the hotspot and relatively high velocities in the north of the stable craton. The same pattern of velocity variation is also imaged in the 3-D shear wave model. The lowest velocity tilted to north and east to the current hotspot location with depth, reflecting shearing of the hotspot material due to the relative motion of the North American plate. In addition, a low velocity anomaly is imaged at central-eastern portion of the craton at depths greater than 150 km, correlated with a thin crust. Relative high velocity anomaly is found in between this slow anomaly and the hotspot. This slow-fast-slow velocity variation in the upper mantle could be associated with small-scale mantle convection of upwelling-downwelling-upwelling trigged by the hotspot, which could thermally erode the cratonic lithosphere.
NASA Astrophysics Data System (ADS)
Chanishvili, R.; Chagelishvili, G.; Uchava, E.; Kharshiladze, O.
2016-04-01
Our goal is to gain new insight into the physics of wave dynamics in ionospheric zonal shear flows. We study the shear flow non-normality induced linear coupling of planetary scale (slow) modified Rossby waves and westward propagating fast magnetized (Khantadze) waves using an approach different from the existing one to the linear wave dynamics. The performed analysis allows us to separate from each other different physical processes, grasp their interplay, and, by this way, construct the basic physics of the linear coupling of the slow and fast waves in an ionospheric zonal flow with linear shear of mean velocity, U0=(S y ,0 ) . It should be noted from the beginning that we consider incompressible flow and the classified "slow" and "fast" waves are not connected with the similarly labeled magnetosonic waves in compressible heliosphere. We show that: the modified Rossby waves generate fast magnetized waves due to the coupling for a quite wide range of ionospheric and shear flow parameters; the linear transient processes are highly anisotropic in wavenumber plane; the generation of the magnetized waves/oscillations is most efficient/optimal for S ≃0.1 (S is the shear rate normalized to the combination of the angular velocity and latitude, Ω0 cos θ0 ); the streamwise wave number of the optimally generated magnetized wave harmonics decreases (the length scale increases) with increasing the Hall parameter, α. At the end, we discuss nonlinear consequences of the described anisotropic linear dynamics—they should lead to an anisotropy of nonlinear cascade processes (in wavenumber plane). In turn, an interplay of the analyzed quite strong transient growth of the fast magnetic waves with anisotropic nonlinear processes should ensure self-sustenance of (stochastic or regular) magnetic perturbations.
The Hartle-Hawking wave function in 2D causal set quantum gravity
NASA Astrophysics Data System (ADS)
Glaser, Lisa; Surya, Sumati
2016-03-01
We define the Hartle-Hawking no-boundary wave function for causal set theory (CST) over the discrete analogs of spacelike hypersurfaces. Using Markov Chain Monte Carlo and numerical integration methods we analyze the wave function in non-perturbative 2D CST. We find that in the low-temperature regime it is dominated by causal sets which have no continuum counterparts but possess physically interesting geometric properties. Not only do they exhibit a rapid spatial expansion with respect to the discrete proper time, but a high degree of spatial homogeneity. The latter is due to the extensive overlap of the causal pasts of the elements in the final discrete hypersurface and corresponds to high graph connectivity. Our results thus suggest new possibilities for the role of quantum gravity in the observable Universe.
Instability of subharmonic resonances in magnetogravity shear waves
NASA Astrophysics Data System (ADS)
Salhi, A.; Nasraoui, S.
2013-12-01
We study analytically the instability of the subharmonic resonances in magnetogravity waves excited by a (vertical) time-periodic shear for an inviscid and nondiffusive unbounded conducting fluid. Due to the fact that the magnetic potential induction is a Lagrangian invariant for magnetohydrodynamic Euler-Boussinesq equations, we show that plane-wave disturbances are governed by a four-dimensional Floquet system in which appears, among others, the parameter ɛ representing the ratio of the periodic shear amplitude to the vertical Brunt-Väisälä frequency N3. For sufficiently small ɛ and when the magnetic field is horizontal, we perform an asymptotic analysis of the Floquet system following the method of Lebovitz and Zweibel [Astrophys. J. 609, 301 (2004), 10.1086/420972]. We determine the width and the maximal growth rate of the instability bands associated with subharmonic resonances. We show that the instability of subharmonic resonance occurring in gravity shear waves has a maximal growth rate of the form Δm=(3√3 /16)ɛ. This instability persists in the presence of magnetic fields, but its growth rate decreases as the magnetic strength increases. We also find a second instability involving a mixing of hydrodynamic and magnetic modes that occurs for all magnetic field strengths. We also elucidate the similarity between the effect of a vertical magnetic field and the effect of a vertical Coriolis force on the gravity shear waves considering axisymmetric disturbances. For both cases, plane waves are governed by a Hill equation, and, when ɛ is sufficiently small, the subharmonic instability band is determined by a Mathieu equation. We find that, when the Coriolis parameter (or the magnetic strength) exceeds N3/2, the instability of the subharmonic resonance vanishes.
Instability of subharmonic resonances in magnetogravity shear waves.
Salhi, A; Nasraoui, S
2013-12-01
We study analytically the instability of the subharmonic resonances in magnetogravity waves excited by a (vertical) time-periodic shear for an inviscid and nondiffusive unbounded conducting fluid. Due to the fact that the magnetic potential induction is a Lagrangian invariant for magnetohydrodynamic Euler-Boussinesq equations, we show that plane-wave disturbances are governed by a four-dimensional Floquet system in which appears, among others, the parameter ɛ representing the ratio of the periodic shear amplitude to the vertical Brunt-Väisälä frequency N(3). For sufficiently small ɛ and when the magnetic field is horizontal, we perform an asymptotic analysis of the Floquet system following the method of Lebovitz and Zweibel [Astrophys. J. 609, 301 (2004)]. We determine the width and the maximal growth rate of the instability bands associated with subharmonic resonances. We show that the instability of subharmonic resonance occurring in gravity shear waves has a maximal growth rate of the form Δ(m)=(3√[3]/16)ɛ. This instability persists in the presence of magnetic fields, but its growth rate decreases as the magnetic strength increases. We also find a second instability involving a mixing of hydrodynamic and magnetic modes that occurs for all magnetic field strengths. We also elucidate the similarity between the effect of a vertical magnetic field and the effect of a vertical Coriolis force on the gravity shear waves considering axisymmetric disturbances. For both cases, plane waves are governed by a Hill equation, and, when ɛ is sufficiently small, the subharmonic instability band is determined by a Mathieu equation. We find that, when the Coriolis parameter (or the magnetic strength) exceeds N(3)/2, the instability of the subharmonic resonance vanishes. PMID:24483566
2D aperture synthesis for Lamb wave imaging using co-arrays
NASA Astrophysics Data System (ADS)
Ambrozinski, Lukasz; Stepinski, Tadeusz; Uhl, Tadeusz
2014-03-01
2D ultrasonic arrays in Lamb wave based SHM systems can operate in the phased array (PA) or synthetic focusing (SF) mode. In the real-time PA approach, multiple electronically delayed signals excite transmitting elements to form the desired wave-front, whereas receiving elements are used to sense scattered waves. Due to that, the PA mode requires multi channeled hardware and multiple excitations at numerous azimuths to scan the inspected region of interest. To the contrary, the SF mode, assumes a single element excitation of subsequent transmitters and off-line processing of the acquired data. In the simplest implementation of the SF technique, a single multiplexed input and output channels are required, which results in significant hardware simplification. Performance of a 2D imaging array depends on many parameters, such as, its topology, number of its transducers and their spacing in terms of wavelength as well as the type of weighting function (apodization). Moreover, it is possible to use sparse arrays, which means that not all array elements are used for transmitting and/ or receiving. In this paper the co-array concept is applied to facilitate the synthesis process of an array's aperture used in the multistatic synthetic focusing approach in Lamb waves-based imaging systems. In the coherent imaging, performed in the transmit/receive mode, the sum co-array is a morphological convolution of the transmit/receive sub-arrays. It can be calculated as the set of sums of the individual elements' locations in the sub-arrays used for imaging. The coarray framework will be presented here using two different array topologies, aID uniform linear array and a cross-shaped array that will result in a square coarray. The approach will be discussed in terms of array patterns and beam patterns of the resulting imaging systems. Both, theoretical and experimental results will be given.
Accelerating numerical modeling of wave propagation through 2-D anisotropic materials using OpenCL.
Molero, Miguel; Iturrarán-Viveros, Ursula
2013-03-01
We present an implementation of the numerical modeling of elastic waves propagation, in 2D anisotropic materials, using the new parallel computing devices (PCDs). Our study is aimed both to model laboratory experiments and explore the capabilities of the emerging PCDs by discussing performance issues. In the experiments a sample plate of an anisotropic material placed inside a water tank is rotated and, for every angle of rotation it is subjected to an ultrasonic wave (produced by a large source transducer) that propagates in the water and through the material producing some reflection and transmission signals that are recording by a "point-like" receiver. This experiment is numerically modeled by running a finite difference code covering a set of angles θ∈[-50°, 50°], and recorded the signals for the transmission and reflection results. Transversely anisotropic and weakly orthorhombic materials are considered. We accelerated the computation using an open-source toolkit called PyOpenCL, which lets one to easily access the OpenCL parallel computation API's from the high-level programming environment of Python. A speedup factor over 19 using the GPU is obtained when compared with the execution of the same program in parallel using a CPU multi-core (in this case we use the 4-cores that has the CPU). The performance for different graphic cards and operating systems is included together with the full 2-D finite difference code with PyOpenCL. PMID:23290584
Probing the shear-band formation in granular media with sound waves.
Khidas, Y; Jia, X
2012-05-01
We investigate the mechanical responses of dense granular materials, using a direct shear box combined with simultaneous acoustic measurements. Measured shear wave speeds evidence the structural change of the material under shear, from the jammed state to the flowing state. There is a clear acoustic signature when the shear band is formed. Subjected to cyclic shear, both shear stress and wave speed show the strong hysteretic dependence on the shear strain, likely associated with the geometry change in the packing structure. Moreover, the correlation function of configuration-specific multiply scattered waves reveals an intermittent behavior before the failure of material. PMID:23004745
Prediction of the Shear Wave Velocity from Compressional Wave Velocity for Gachsaran Formation
NASA Astrophysics Data System (ADS)
Parvizi, Saeed; Kharrat, Riyaz; Asef, Mohammad R.; Jahangiry, Bijan; Hashemi, Abdolnabi
2015-10-01
Shear and compressional wave velocities, coupled with other petrophysical data, are very important for hydrocarbon reservoir characterization. In situ shear wave velocity (Vs) is measured by some sonic logging tools. Shear velocity coupled with compressional velocity is vitally important in determining geomechanical parameters, identifying the lithology, mud weight design, hydraulic fracturing, geophysical studies such as VSP, etc. In this paper, a correlation between compressional and shear wave velocity is obtained for Gachsaran formation in Maroon oil field. Real data were used to examine the accuracy of the prediction equation. Moreover, the genetic algorithm was used to obtain the optimal value for constants of the suggested equation. Furthermore, artificial neural network was used to inspect the reliability of this method. These investigations verify the notion that the suggested equation could be considered as an efficient, fast, and cost-effective method for predicting Vs from Vp.
Shear waves in an inhomogeneous strongly coupled dusty plasma
Janaki, M. S.; Banerjee, D.; Chakrabarti, N.
2011-09-15
The properties of electrostatic transverse shear waves propagating in a strongly coupled dusty plasma with an equilibrium density gradient are examined using the generalized hydrodynamic (GH) equation. In the usual kinetic limit, the resulting equation has similarity to zero energy Schrodinger's equation. This has helped in obtaining some exact eigenmode solutions in both Cartesian and cylindrical geometries for certain nontrivial density profiles. The corresponding velocity profiles and the discrete eigenfrequencies are obtained for several interesting situations and their physics discussed.
Monitoring temporal changes with shear wave splitting: testing the methodology
NASA Astrophysics Data System (ADS)
Walsh, E.; Savage, M. K.; Arnold, R.; Brenguier, F.; Rivemale, E.
2012-12-01
Changes in stress or fluid content of cracks should cause changes in seismic anisotropy, which can be monitored via shear wave splitting. However, shear wave splitting measurements often yield larger variation between the results than is expected from the formal error bars generated by computer programs. Furthermore, recent measurements on multiplets of nearly identical events have yielded large variations and "jumps" in apparent splitting parameters as a function of time rather than smooth variations as expected for slow changes in path properties. We carry out a systematic analysis of the Silver and Chan method. We recently reported on variations based on synthetic seismograms with various types of noise added. Here we extend the analysis to examine the variations of shear wave splitting on a family of 513 earthquakes that was recorded on station BOR on Piton de la Fournaise Volcano on La Reunion Island. Earthquakes were assigned to families based on cross correlation of the P waveform of higher than 90%. Of these 513 events, only 268 high quality events were used. A fixed 2 pole Butterworth bandpass filter with corner frequencies of 4 and 10Hz was applied. Several clusters of fast directions appeared so we chose to analyse the cluster with the most (146) events, which had fast directions between 0 and 45 degrees. The analysis window was then fixed at 3.88-4.39 seconds for all events. In the process of examination, we discovered that the cause of the jumps is most likely related to cycle skipping. Cycle skipping in shear wave splitting is a phenomenon that is known to lead to delay time shifts that are half integer multiples of the dominant period and can create 90 degree errors in the fast direction. Using plots of particle motion, the polarisation components and delay times against incoming polarisation, our preliminary findings show that 90 degree errors in the incoming polarisation can lead to correctly identifying the fast direction, but produces jumps in
Shear-driven Dynamo Waves in the Fully Nonlinear Regime
NASA Astrophysics Data System (ADS)
Pongkitiwanichakul, P.; Nigro, G.; Cattaneo, F.; Tobias, S. M.
2016-07-01
Large-scale dynamo action is well understood when the magnetic Reynolds number (Rm) is small, but becomes problematic in the astrophysically relevant large Rm limit since the fluctuations may control the operation of the dynamo, obscuring the large-scale behavior. Recent works by Tobias & Cattaneo demonstrated numerically the existence of large-scale dynamo action in the form of dynamo waves driven by strongly helical turbulence and shear. Their calculations were carried out in the kinematic regime in which the back-reaction of the Lorentz force on the flow is neglected. Here, we have undertaken a systematic extension of their work to the fully nonlinear regime. Helical turbulence and large-scale shear are produced self-consistently by prescribing body forces that, in the kinematic regime, drive flows that resemble the original velocity used by Tobias & Cattaneo. We have found four different solution types in the nonlinear regime for various ratios of the fluctuating velocity to the shear and Reynolds numbers. Some of the solutions are in the form of propagating waves. Some solutions show large-scale helical magnetic structure. Both waves and structures are permanent only when the kinetic helicity is non-zero on average.
Laboratory measurements of compressional and shear wave speeds through methane hydrate
Waite, W.F.; Helgerud, M.B.; Nur, A.; Pinkston, J.C.; Stern, L.A.; Kirby, S.H.; Durham, W.B.
2000-01-01
Simultaneous measurements of compressional and shear wave speeds through polycrystalline methane hydrate have been made. Methane hydrate, grown directly in a wave speed measurement chamber, was uniaxially compacted to a final porosity below 2%. At 277 K, the compacted material had a compressional wave speed of 3650 ?? 50 m/s. The shear wave speed, measured simultaneously, was 1890 ?? 30 m/s. From these wave speed measurements, we derive V(p)/V(s), Poisson's ratio, bulk, shear, and Young's moduli.
Shear wave velocity structure in West Java, Indonesia as inferred from surface wave dispersion
NASA Astrophysics Data System (ADS)
Anggono, Titi; Syuhada
2016-02-01
We investigated the crust and upper mantle of West Java, Indonesia by measuring the group velocity dispersion of surface waves. We analyzed waveform from four teleseismic earthquake recorded at three 3-component broadband seismometers. We analyzed fundamental mode of Rayleigh and Love waves from vertical, radial, and transverse components using multiple filter technique. We inverted the measured group velocity to obtain shear wave velocity profile down to 200 km depth. We observed low shear wave velocity zone at depth of about 20 km. Shear velocity reduction is estimated to be 18% compared to the upper and lower velocity layer. The low velocity zone might be associated with the subducting slab of Indo-Australian Plate as similar characteristics of low velocity zones also observed at other subducting regions.
NASA Technical Reports Server (NTRS)
Marvin, Joseph G.; Brown, James L.; Gnoffo, Peter A.
2013-01-01
A database compilation of hypersonic shock-wave/turbulent boundary layer experiments is provided. The experiments selected for the database are either 2D or axisymmetric, and include both compression corner and impinging type SWTBL interactions. The strength of the interactions range from attached to incipient separation to fully separated flows. The experiments were chosen based on criterion to ensure quality of the datasets, to be relevant to NASA's missions and to be useful for validation and uncertainty assessment of CFD Navier-Stokes predictive methods, both now and in the future. An emphasis on datasets selected was on surface pressures and surface heating throughout the interaction, but include some wall shear stress distributions and flowfield profiles. Included, for selected cases, are example CFD grids and setup information, along with surface pressure and wall heating results from simulations using current NASA real-gas Navier-Stokes codes by which future CFD investigators can compare and evaluate physics modeling improvements and validation and uncertainty assessments of future CFD code developments. The experimental database is presented tabulated in the Appendices describing each experiment. The database is also provided in computer-readable ASCII files located on a companion DVD.
Zhu, Jiang; Qu, Yueqiao; Ma, Teng; Li, Rui; Du, Yongzhao; Huang, Shenghai; Shung, K Kirk; Zhou, Qifa; Chen, Zhongping
2015-05-01
We report on a novel acoustic radiation force orthogonal excitation optical coherence elastography (ARFOE-OCE) technique for imaging shear wave and quantifying shear modulus under orthogonal acoustic radiation force (ARF) excitation using the optical coherence tomography (OCT) Doppler variance method. The ARF perpendicular to the OCT beam is produced by a remote ultrasonic transducer. A shear wave induced by ARF excitation propagates parallel to the OCT beam. The OCT Doppler variance method, which is sensitive to the transverse vibration, is used to measure the ARF-induced vibration. For analysis of the shear modulus, the Doppler variance method is utilized to visualize shear wave propagation instead of Doppler OCT method, and the propagation velocity of the shear wave is measured at different depths of one location with the M scan. In order to quantify shear modulus beyond the OCT imaging depth, we move ARF to a deeper layer at a known step and measure the time delay of the shear wave propagating to the same OCT imaging depth. We also quantitatively map the shear modulus of a cross-section in a tissue-equivalent phantom after employing the B scan. PMID:25927794
Zhu, Jiang; Qu, Yueqiao; Ma, Teng; Li, Rui; Du, Yongzhao; Huang, Shenghai; Shung, K. Kirk; Zhou, Qifa; Chen, Zhongping
2015-01-01
We report on a novel acoustic radiation force orthogonal excitation optical coherence elastography (ARFOE-OCE) technique for imaging shear wave and quantifying shear modulus under orthogonal acoustic radiation force (ARF) excitation using the optical coherence tomography (OCT) Doppler variance method. The ARF perpendicular to the OCT beam is produced by a remote ultrasonic transducer. A shear wave induced by ARF excitation propagates parallel to the OCT beam. The OCT Doppler variance method, which is sensitive to the transverse vibration, is used to measure the ARF-induced vibration. For analysis of the shear modulus, the Doppler variance method is utilized to visualize shear wave propagation instead of Doppler OCT method, and the propagation velocity of the shear wave is measured at different depths of one location with the M scan. In order to quantify shear modulus beyond the OCT imaging depth, we move ARF to a deeper layer at a known step and measure the time delay of the shear wave propagating to the same OCT imaging depth. We also quantitatively map the shear modulus of a cross-section in a tissue-equivalent phantom after employing the B scan. PMID:25927794
2D time-domain finite-difference modeling for viscoelastic seismic wave propagation
NASA Astrophysics Data System (ADS)
Fan, Na; Zhao, Lian-Feng; Xie, Xiao-Bi; Ge, Zengxi; Yao, Zhen-Xing
2016-07-01
Real Earth media are not perfectly elastic. Instead, they attenuate propagating mechanical waves. This anelastic phenomenon in wave propagation can be modeled by a viscoelastic mechanical model consisting of several standard linear solids. Using this viscoelastic model, we approximate a constant Q over a frequency band of interest. We use a four-element viscoelastic model with a tradeoff between accuracy and computational costs to incorporate Q into 2D time-domain first-order velocity-stress wave equations. To improve the computational efficiency, we limit the Q in the model to a list of discrete values between 2 and 1000. The related stress and strain relaxation times that characterize the viscoelastic model are pre-calculated and stored in a database for use by the finite-difference calculation. A viscoelastic finite-difference scheme that is second-order in time and fourth-order in space is developed based on the MacCormack algorithm. The new method is validated by comparing the numerical result with analytical solutions that are calculated using the generalized reflection/transmission coefficient method. The synthetic seismograms exhibit greater than 95 per cent consistency in a two-layer viscoelastic model. The dispersion generated from the simulation is consistent with the Kolsky-Futterman dispersion relationship.
2D fluid simulations of acoustic waves in pulsed ICP discharges: Comparison with experiments
NASA Astrophysics Data System (ADS)
Despiau-Pujo, Emilie; Cunge, Gilles; Sadeghi, Nader; Braithwaite, N. St. J.
2012-10-01
Neutral depletion, which is mostly caused by gas heating under typical material processing conditions, is an important phenomenon in high-density plasmas. In low pressure pulsed discharges, experiments show that additional depletion due to electron pressure (Pe) may have a non-negligible influence on radical transport [1]. To evaluate this effect, comparisons between 2D fluid simulations and measurements of gas convection in Ar/Cl2 pulsed ICP plasmas are reported. In the afterglow, Pe drops rapidly by electron cooling which generates a neutral pressure gradient between the plasma bulk and the reactor walls. This in turn forces the cold surrounding gas to move rapidly towards the center, thus launching an acoustic wave in the reactor. Time-resolved measurements of atoms drift velocity and gas temperature by LIF and LAS in the early afterglow are consistent with gas drifting at acoustic wave velocity followed by rapid gas cooling. Similar results are predicted by the model. The ion flux at the reactor walls is also shown to oscillate in phase with the acoustic wave due to ion-neutral friction forces. Finally, during plasma ignition, experiments show opposite phenomena when Pe rises.[4pt] [1] Cunge et al, APL 96, 131501 (2010)
Optimal implicit 2-D finite differences to model wave propagation in poroelastic media
NASA Astrophysics Data System (ADS)
Itzá, Reymundo; Iturrarán-Viveros, Ursula; Parra, Jorge O.
2016-08-01
Numerical modeling of seismic waves in heterogeneous porous reservoir rocks is an important tool for the interpretation of seismic surveys in reservoir engineering. We apply globally optimal implicit staggered-grid finite differences (FD) to model 2-D wave propagation in heterogeneous poroelastic media at a low-frequency range (<10 kHz). We validate the numerical solution by comparing it to an analytical-transient solution obtaining clear seismic wavefields including fast P and slow P and S waves (for a porous media saturated with fluid). The numerical dispersion and stability conditions are derived using von Neumann analysis, showing that over a wide range of porous materials the Courant condition governs the stability and this optimal implicit scheme improves the stability of explicit schemes. High-order explicit FD can be replaced by some lower order optimal implicit FD so computational cost will not be as expensive while maintaining the accuracy. Here, we compute weights for the optimal implicit FD scheme to attain an accuracy of γ = 10-8. The implicit spatial differentiation involves solving tridiagonal linear systems of equations through Thomas' algorithm.
Optimal implicit 2-D finite differences to model wave propagation in poroelastic media
NASA Astrophysics Data System (ADS)
Itzá, Reymundo; Iturrarán-Viveros, Ursula; Parra, Jorge O.
2016-05-01
Numerical modeling of seismic waves in heterogeneous porous reservoir rocks is an important tool for the interpretation of seismic surveys in reservoir engineering. We apply globally optimal implicit staggered-grid finite-differences to model 2-D wave propagation in heterogeneous poroelastic media at a low-frequency range (<10kHz). We validate the numerical solution by comparing it to an analytical-transient solution obtaining clear seismic wavefields including fast P, slow P and S waves (for a porous media saturated with fluid). The numerical dispersion and stability conditions are derived using von Neumann analysis, showing that over a wide range of porous materials the Courant condition governs the stability and this optimal implicit scheme improves the stability of explicit schemes. High order explicit finite-differences (FD) can be replaced by some lower order optimal implicit FD so computational cost will not be as expensive while maintaining the accuracy. Here we compute weights for the optimal implicit FD scheme to attain an accuracy of γ = 10-8. The implicit spatial differentiation involves solving tridiagonal linear systems of equations through Thomas' algorithm.
Schmidt, J L; Tweten, D J; Benegal, A N; Walker, C H; Portnoi, T E; Okamoto, R J; Garbow, J R; Bayly, P V
2016-05-01
Mechanical anisotropy is an important property of fibrous tissues; for example, the anisotropic mechanical properties of brain white matter may play a key role in the mechanics of traumatic brain injury (TBI). The simplest anisotropic material model for small deformations of soft tissue is a nearly incompressible, transversely isotropic (ITI) material characterized by three parameters: minimum shear modulus (µ), shear anisotropy (ϕ=µ1µ-1) and tensile anisotropy (ζ=E1E2-1). These parameters can be determined using magnetic resonance elastography (MRE) to visualize shear waves, if the angle between the shear-wave propagation direction and fiber direction is known. Most MRE studies assume isotropic material models with a single shear (µ) or tensile (E) modulus. In this study, two types of shear waves, "fast" and "slow", were analyzed for a given propagation direction to estimate anisotropic parameters µ, ϕ, and ζ in two fibrous soft materials: turkey breast ex vivo and aligned fibrin gels. As expected, the speed of slow shear waves depended on the angle between fiber direction and propagation direction. Fast shear waves were observed when the deformations due to wave motion induced stretch in the fiber direction. Finally, MRE estimates of anisotropic mechanical properties in turkey breast were compared to estimates from direct mechanical tests. PMID:26920505
Validation of Shear Wave Elastography in Skeletal Muscle
Eby, Sarah F.; Song, Pengfei; Chen, Shigao; Chen, Qingshan; Greenleaf, James F.; An, Kai-Nan
2013-01-01
Skeletal muscle is a very dynamic tissue, thus accurate quantification of skeletal muscle stiffness throughout its functional range is crucial to improve the physical functioning and independence following pathology. Shear wave elastography (SWE) is an ultrasound-based technique that characterizes tissue mechanical properties based on the propagation of remotely induced shear waves. The objective of this study is to validate SWE throughout the functional range of motion of skeletal muscle for three ultrasound transducer orientations. We hypothesized that combining traditional materials testing (MTS) techniques with SWE measurements will show increased stiffness measures with increasing tensile load, and will correlate well with each other for trials in which the transducer is parallel to underlying muscle fibers. To evaluate this hypothesis, we monitored the deformation throughout tensile loading of four porcine brachialis whole-muscle tissue specimens, while simultaneously making SWE measurements of the same specimen. We used regression to examine the correlation between Young's modulus from MTS and shear modulus from SWE for each of the transducer orientations. We applied a generalized linear model to account for repeated testing. Model parameters were estimated via generalized estimating equations. The regression coefficient was 0.1944, with a 95% confidence interval of (0.1463 – 0.2425) for parallel transducer trials. Shear waves did not propagate well for both the 45° and perpendicular transducer orientations. Both parallel SWE and MTS showed increased stiffness with increasing tensile load. This study provides the necessary first step for additional studies that can evaluate the distribution of stiffness throughout muscle. PMID:23953670
Shape waves in 2D Josephson junctions: exact solutions and time dilation.
Gulevich, D R; Kusmartsev, F V; Savel'ev, Sergey; Yampol'skii, V A; Nori, Franco
2008-09-19
We predict a new class of excitations propagating along a Josephson vortex in two-dimensional Josephson junctions. These excitations are associated with the distortion of a Josephson vortex line and have an analogy with shear waves in solid mechanics. Their shapes can have an arbitrary profile, which is retained when propagating. We derive a universal analytical expression for the energy of arbitrary shape excitations, investigate their influence on the dynamics of a vortex line, and discuss conditions where such excitations can be created. Finally, we show that such excitations play the role of a clock for a relativistically moving Josephson vortex and suggest an experiment to measure a time dilation effect analogous to that in special relativity. PMID:18851404
Shape Waves in 2D Josephson Junctions: Exact Solutions and Time Dilation
Gulevich, D. R.; Savel'ev, Sergey; Kusmartsev, F. V.; Yampol'skii, V. A.; Nori, Franco
2008-09-19
We predict a new class of excitations propagating along a Josephson vortex in two-dimensional Josephson junctions. These excitations are associated with the distortion of a Josephson vortex line and have an analogy with shear waves in solid mechanics. Their shapes can have an arbitrary profile, which is retained when propagating. We derive a universal analytical expression for the energy of arbitrary shape excitations, investigate their influence on the dynamics of a vortex line, and discuss conditions where such excitations can be created. Finally, we show that such excitations play the role of a clock for a relativistically moving Josephson vortex and suggest an experiment to measure a time dilation effect analogous to that in special relativity.
Unbounded wall flow with free surface waves and horizontal shear
NASA Astrophysics Data System (ADS)
Lapham, Gary; McHugh, John
2015-11-01
Free surface waves in the presence of a non-uniform shear flow are treated. The shear flow of interest varies with both the transverse and vertical coordinates, U (y , z) . Initial results treat a mean flow varying only with the transverse, U (y) . The domain is bounded on one side by a flat rigid vertical wall and is unbounded on the other side. The mean flows considered here are nonzero near the vertical wall and approach zero far from the wall, e.g. U =e-γy . The flowfield is treated as inviscid but rotational. Linear solutions are obtained using a nonuniform coordinate transformation that converts the free surface boundary condition into a modified Bessel equation. Velocity components are expanded in modified Bessel functions of the first kind of purely imaginary order. The dispersion relation for steady waves are found with wavespeeds outside the range of U, matching previous results for a flow bounded on both sides. Corresponding eigenvectors show a sequence of wave profiles of increasing complexity near the wall. The wave amplitude approaches zero far from the wall.
The stability of freely-propagating ion acoustic waves in 2D systems
NASA Astrophysics Data System (ADS)
Chapman, Thomas; Berger, Richard; Banks, Jeffrey; Brunner, Stephan
2014-10-01
The stability of a freely-propagating ion acoustic wave (IAW) is a basic science problem that is made difficult by the need to resolve electron kinetic effects over a timescale that greatly exceeds the IAW period during numerical simulation. Recent results examining IAW stability using a 1D+1V Vlasov-Poisson solver indicate that instability is a fundamental property of IAWs that occurs over most if not all of the parameter space of relevance to ICF experiments. We present here new results addressing the fundamental question of IAW stability across a broad range of plasma conditions in a 2D+2V system using LOKI, ranging from a regime of relatively weak to a regime of relatively strong ion kinetic effects. Work performed under the auspices of the U.S. DOE by LLNL (DE-AC52-07NA27344) and funded by the LDRD Program at LLNL (12-ERD-061).
2D full wave modeling for a synthetic Doppler backscattering diagnostic
Hillesheim, J. C.; Schmitz, L.; Kubota, S.; Rhodes, T. L.; Carter, T. A.; Holland, C.
2012-10-15
Doppler backscattering (DBS) is a plasma diagnostic used in tokamaks and other magnetic confinement devices to measure the fluctuation level of intermediate wavenumber (k{sub {theta}}{rho}{sub s}{approx} 1) density fluctuations and the lab frame propagation velocity of turbulence. Here, a synthetic DBS diagnostic is described, which has been used for comparisons between measurements in the DIII-D tokamak and predictions from nonlinear gyrokinetic simulations. To estimate the wavenumber range to which a Gaussian beam would be sensitive, a ray tracing code and a 2D finite difference, time domain full wave code are used. Experimental density profiles and magnetic geometry are used along with the experimental antenna and beam characteristics. An example of the effect of the synthetic diagnostic on the output of a nonlinear gyrokinetic simulation is presented.
Assessment of local pulse wave velocity in arteries using 2D distension waveforms.
Meinders, J M; Kornet, L; Brands, P J; Hoeks, A P
2001-10-01
The reciprocal of the arterial pulse wave velocity contains crucial information about the mechanical characteristics of the arterial wall but is difficult to assess noninvasively in vivo. In this paper, a new method to assess local pulse wave velocity (PWV) is presented. To this end, multiple adjacent distension waveforms are determined simultaneously along a short arterial segment, using a single 2D-vessel wall tracking system with a high frame rate (651 Hz). Each B-mode image consists of 16 echo lines spanning a total width of 15.86 mm. Dedicated software has been developed to extract the end-diastolic diameter from the B-mode image and the distension waveforms from the underlying radiofrequency (rf) information for each echo-line. The PWV is obtained by determining the ratio of the temporal and spatial gradient of adjacent distension velocity waveforms. The proposed method is verified in a phantom and in the common carotid artery (CCA) of humans. Phantom experiments show a high concordance between the PWV obtained from 2D distension velocity waveforms (4.21 +/- 0.02 m/s) and the PWV determined using two pressure catheters (4.26 +/- 0.02 m/s). Assuming linear spatial gradients, the PWV can also be obtained in vivo for CCA and averages to 5.5 +/- 1.5 m/s (intersubject variation, n = 23), which compares well to values found in literature. Furthermore, intrasubject PWV compares well with those calculated using the Bramwell-Hill equation. It can be concluded that the PWV can be obtained from the spatial and temporal gradient if the spatial gradient is linear over the observed length of the artery, i.e. the artery should be homogenous in diameter and distension and the influence of reflections must be small. PMID:12051275
Helfenstein-Didier, C; Andrade, R J; Brum, J; Hug, F; Tanter, M; Nordez, A; Gennisson, J-L
2016-03-21
The shear wave velocity dispersion was analyzed in the Achilles tendon (AT) during passive dorsiflexion using a phase velocity method in order to obtain the tendon shear modulus (C 55). Based on this analysis, the aims of the present study were (i) to assess the reproducibility of the shear modulus for different ankle angles, (ii) to assess the effect of the probe locations, and (iii) to compare results with elasticity values obtained with the supersonic shear imaging (SSI) technique. The AT shear modulus (C 55) consistently increased with the ankle dorsiflexion (N = 10, p < 0.05). Furthermore, the technique showed a very good reproducibility (all standard error of the mean values <10.7 kPa and all coefficient of variation (CV) values ⩽0.05%). In addition, independently from the ankle dorsiflexion, the shear modulus was significantly higher in the proximal location compared to the more distal one. The shear modulus provided by SSI was always lower than C55 and the difference increased with the ankle dorsiflexion. However, shear modulus values provided by both methods were highly correlated (R = 0.84), indicating that the conventional shear wave elastography technique (SSI technique) can be used to compare tendon mechanical properties across populations. Future studies should determine the clinical relevance of the shear wave dispersion analysis, for instance in the case of tendinopathy or tendon tear. PMID:26948399
NASA Astrophysics Data System (ADS)
Helfenstein-Didier, C.; Andrade, R. J.; Brum, J.; Hug, F.; Tanter, M.; Nordez, A.; Gennisson, J.-L.
2016-03-01
The shear wave velocity dispersion was analyzed in the Achilles tendon (AT) during passive dorsiflexion using a phase velocity method in order to obtain the tendon shear modulus (C 55). Based on this analysis, the aims of the present study were (i) to assess the reproducibility of the shear modulus for different ankle angles, (ii) to assess the effect of the probe locations, and (iii) to compare results with elasticity values obtained with the supersonic shear imaging (SSI) technique. The AT shear modulus (C 55) consistently increased with the ankle dorsiflexion (N = 10, p < 0.05). Furthermore, the technique showed a very good reproducibility (all standard error of the mean values <10.7 kPa and all coefficient of variation (CV) values ⩽0.05%). In addition, independently from the ankle dorsiflexion, the shear modulus was significantly higher in the proximal location compared to the more distal one. The shear modulus provided by SSI was always lower than C55 and the difference increased with the ankle dorsiflexion. However, shear modulus values provided by both methods were highly correlated (R = 0.84), indicating that the conventional shear wave elastography technique (SSI technique) can be used to compare tendon mechanical properties across populations. Future studies should determine the clinical relevance of the shear wave dispersion analysis, for instance in the case of tendinopathy or tendon tear.
The Role of Shear and Interface Waves in the Excitation of T-waves
NASA Astrophysics Data System (ADS)
Odom, R. I.
2007-12-01
T-waves are late arriving phases on a seismogram which travel a significant portion of their journey from the source to the receiver along a slow oceanic path. Park, Soukup and Odom [2001] proposed a modal scattering mechanism which permits energy from steeply arriving rays impinging on the ocean bottom to be converted to shallow grazing angle rays corresponding to low order modes known to comprise the T-wave signal. The Scholte interface waves are crucial to this energy transfer as they have anti-nodes nearly coincident with the ocean bottom. Any seafloor roughness acts as a secondary source located right on the Scholte wave anti-node. This allows energy to be efficiently transferred from high energy modes to lower energy, and slower, T-wave modes. In fact the presence of the Scholte waves is crucial to the existence of T-waves. Some finite shear strength in the bottom sediments and/or upper ocean crust is essential for the existence of the Scholte waves. Elastic interface waves do not exist the boundary between two fluids. The effect of the the shear modulus of the ocean bottom sediment and ocean crust on the excitation of the interface waves and T-waves is discussed.
Shear Wave Splitting Inversion in a Complex Crust
NASA Astrophysics Data System (ADS)
Lucas, A.
2015-12-01
Shear wave splitting (SWS) inversion presents a method whereby the upper crust can be interrogated for fracture density. It is caused when a shear wave traverses an area of anisotropy, splits in two, with each wave experiencing a different velocity resulting in an observable separation in arrival times. A SWS observation consists of the first arrival polarization direction and the time delay. Given the large amount of data common in SWS studies, manual inspection for polarization and time delay is considered prohibitively time intensive. All automated techniques used can produce high amounts of observations falsely interpreted as SWS. Thus introducing error into the interpretation. The technique often used for removing these false observations is to manually inspect all SWS observations defined as high quality by the automated routine, and remove false identifications. We investigate the nature of events falsely identified compared to those correctly identified. Once this identification is complete we conduct a inversion for crack density from SWS time delay. The current body of work on linear SWS inversion utilizes an equation that defines the time delay between arriving shear waves with respect to fracture density. This equation makes the assumption that no fluid flow occurs as a result of the passing shear wave, a situation called squirt flow. We show that the assumption is not applicable in all geological situations. When it is not true, its use in an inversion produces a result which is negatively affected by the assumptions. This is shown to be the case at the test case of 6894 SWS observations gathered in a small area at Puna geothermal field, Hawaii. To rectify this situation, a series of new time delay formulae, applicable to linear inversion, are derived from velocity equations presented in literature. The new formula use a 'fluid influence parameter' which indicates the degree to which squirt flow is influencing the SWS. It is found that accounting for
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)
Song, Shaozhen; Joy, Joyce; Wang, Ruikang K.; Huang, Zhihong
2015-03-01
A quantitative measurement of the mechanical properties of biological tissue is a useful assessment of its physiologic conditions, which may aid medical diagnosis and treatment of, e.g., scleroderma and skin cancer. Traditional elastography techniques such as magnetic resonance elastography and ultrasound elastography have limited scope of application on skin due to insufficient spatial resolution. Recently, dynamic / transient elastography are attracting more applications with the advantage of non-destructive measurements, and revealing the absolute moduli values of tissue mechanical properties. Shear wave optical coherence elastography (SW-OCE) is a novel transient elastography method, which lays emphasis on the propagation of dynamic mechanical waves. In this study, high speed shear wave imaging technique was applied to a range of soft-embalmed mouse skin, where 3 kHz shear waves were launched with a piezoelectric actuator as an external excitation. The shear wave velocity was estimated from the shear wave images, and used to recover a shear modulus map in the same OCT imaging range. Results revealed significant difference in shear modulus and structure in compliance with gender, and images on fresh mouse skin are also compared. Thiel embalming technique is also proven to present the ability to furthest preserve the mechanical property of biological tissue. The experiment results suggest that SW-OCE is an effective technique for quantitative estimation of skin tissue biomechanical status.
Simulations of SH wave scattering due to cracks by the 2-D finite difference method
NASA Astrophysics Data System (ADS)
Suzuki, Y.; Kawahara, J.; Okamoto, T.; Miyashita, K.
2006-05-01
We simulate SH wave scattering by 2-D parallel cracks using the finite difference method (FDM), instead of the popularly used boundary integral equation method (BIEM). Here special emphasis is put on simplicity; we apply a standard FDM (fourth-order velocity-stress scheme with a staggered grid) to media in cluding traction-freecracks, which are expressed by arrays of grid points with zero traction. Two types of accuracy tests based oncomparison with a reliable BIEM, suggest that the present method gives practically sufficient accuracy, except for the wavefields in the vicinity of cracks, which can be well handled if the second-order FDM is used instead. As an application of this method, we also simulate wave propagation in media with randomly distributed cracks of the same length. We experimentally determine the attenuation and velocity dispersion induced by scattering from the synthetic seismograms, using a waveform averaging technique. It is shown that the results are well explained by a theory based on the Foldy approximation for crack densities of up to about 01. The presence of a free surface does not affect the validity of the theory. A preliminary experiment also suggests that the validity will not change even for multi-scale cracks.
Electrostatic drift waves in a 2D magnetic current sheet - a new kinetic theory
NASA Astrophysics Data System (ADS)
Fruit, G.; Louarn, P.; Tur, A.
2015-12-01
In the general context of understanding the possible destabilization of the magnetotail before a substorm, a kinetic model for electromagnetic instabilities in resonant interaction with trapped bouncing electrons has been proposed for several years. Fruit et al. 2013 already used it to investigate the possibilities for electrostatic instabilities. Tur et al. 2014 generalizes the model for full electromagnetic perturbations.It turns out that some corrections should be added to the electrostatic version of Fruit et al. 2013. We propose to revist the theory in this present paper.Starting with a modified 2D Harris sheet as equilibrium state, the linearized gyrokinetic Vlasov equation is solved for electrostatic fluctuations with period of the order of the electron bounce period (a few seconds). The particle motion is restricted to its first Fourier component along the magnetic field and this allows the complete time integration of the non local perturbed distribution functions. The dispersion relation for electrostatic modes is finally obtained through the quasineutrality condition.The new feature of the present model is the inclusion of diamagnetic drift effects due to the density gradient in the tail. It is well known in MHD theory that drift waves are driven unstable through collisions or other dissipative effects. Here electrostatic drift waves are revisited in this more complete kinetic model including bouncing electrons and finite Larmor radius effects. A new mode has been found with original propagation proprieties. It is moreover mildly unstable due to electron or ion damping (dissipative instability).
1D and 2D simulations of seismic wave propagation in fractured media
NASA Astrophysics Data System (ADS)
Möller, Thomas; Friederich, Wolfgang
2016-04-01
Fractures and cracks have a significant influence on the propagation of seismic waves. Their presence causes reflections and scattering and makes the medium effectively anisotropic. We present a numerical approach to simulation of seismic waves in fractured media that does not require direct modelling of the fracture itself, but uses the concept of linear slip interfaces developed by Schoenberg (1980). This condition states that at an interface between two imperfectly bonded elastic media, stress is continuous across the interface while displacement is discontinuous. It is assumed that the jump of displacement is proportional to stress which implies a jump in particle velocity at the interface. We use this condition as a boundary condition to the elastic wave equation and solve this equation in the framework of a Nodal Discontinuous Galerkin scheme using a velocity-stress formulation. We use meshes with tetrahedral elements to discretise the medium. Each individual element face may be declared as a slip interface. Numerical fluxes have been derived by solving the 1D Riemann problem for slip interfaces with elastic and viscoelastic rheology. Viscoelasticity is realised either by a Kelvin-Voigt body or a Standard Linear Solid. These fluxes are not limited to 1D and can - with little modification - be used for simulations in higher dimensions as well. The Nodal Discontinuous Galerkin code "neXd" developed by Lambrecht (2013) is used as a basis for the numerical implementation of this concept. We present examples of simulations in 1D and 2D that illustrate the influence of fractures on the seismic wavefield. We demonstrate the accuracy of the simulation through comparison to an analytical solution in 1D.
Error in Estimates of Tissue Material Properties from Shear Wave Dispersion Ultrasound Vibrometry
Urban, Matthew W.; Chen, Shigao; Greenleaf, James F.
2009-01-01
Shear wave velocity measurements are used in elasticity imaging to find the shear elasticity and viscosity of tissue. A technique called shear wave dispersion ultrasound vibrometry (SDUV) has been introduced to use the dispersive nature of shear wave velocity to locally estimate the material properties of tissue. Shear waves are created using a multifrequency ultrasound radiation force, and the propagating shear waves are measured a few millimeters away from the excitation point. The shear wave velocity is measured using a repetitive pulse-echo method and Kalman filtering to find the phase of the harmonic shear wave at 2 different locations. A viscoelastic Voigt model and the shear wave velocity measurements at different frequencies are used to find the shear elasticity (μ1) and viscosity (μ2) of the tissue. The purpose of this paper is to report the accuracy of the SDUV method over a range of different values of μ1 and μ2. A motion detection model of a vibrating scattering medium was used to analyze measurement errors of vibration phase in a scattering medium. To assess the accuracy of the SDUV method, we modeled the effects of phase errors on estimates of shear wave velocity and material properties while varying parameters such as shear stiffness and viscosity, shear wave amplitude, the distance between shear wave measurements (Δr), signal-to-noise ratio (SNR) of the ultrasound pulse-echo method, and the frequency range of the measurements. We performed an experiment in a section of porcine muscle to evaluate variation of the aforementioned parameters on the estimated shear wave velocity and material property measurements and to validate the error prediction model. The model showed that errors in the shear wave velocity and material property estimates were minimized by maximizing shear wave amplitude, pulse-echo SNR, Δr, and the bandwidth used for shear wave measurements. The experimental model showed optimum performance could be obtained for Δr = 3-6 mm
NASA Astrophysics Data System (ADS)
Markus Schmalholz, Stefan; Jaquet, Yoann
2016-04-01
We study the formation of an orogenic wedge during lithospheric shortening with 2D numerical simulations. We consider a viscoelastoplastic rheology, thermo-mechanical coupling by shear heating and temperature-dependent viscosities, gravity and erosion. In the initial model configuration there is either a lateral temperature variation at the model base or a lateral variation in crustal thickness to generate slight stress variations during lithospheric shortening. These stress variations can trigger the formation of shear zones which are caused by thermal softening associated with shear heating. We do not apply any kind of strain softening, such as reduction of friction angle with progressive plastic strain. The first major shear zone that appears during shortening crosscuts the entire crust and initiates the asymmetric subduction/underthrusting of mainly the mechanically strong lower crust. After some deformation, the first shear zone in the upper crust is abandoned, the deformation propagates towards the foreland and a new shear zone forms only in the upper crust. The shear zone propagation occurs several times where new shear zones form in the upper crust and the mechanically strong top of the lower crust acts as detachment horizon. We calculate the magnitudes of the maximal and minimal principal stresses and of the mean stress (or dynamic pressure), and we record also the temperature for several marker points in the upper and lower crust. We analyse the evolution of stresses and temperature with burial depth and time. Deviatoric stresses (half the differential stress) in the upper crust are up to 200 MPa and associated shear heating in shear zones ranges between 40 - 80 °C. Lower crustal rocks remain either at the base of the orogenic wedge at depths of around 50 km or are subducted to depths of up to 120 km, depending on their position when the first shear zone formed. Largest deviatotric stresses in the strong part of the lower crust are about 1000 MPa and
Long-wave shear instability of fluid interfaces
NASA Astrophysics Data System (ADS)
Cherniavski, Vladimir
2010-05-01
The earlier oceanographic works largely focused on the case of very large density differences between the two fluids sepa¬rated by an interface. The aim of this investigation is to extend the shear-flow analysis to more wide range of density ratios for the log profile of the "wind". Long wave asymptotic leads to the analytic determination of the stability characteristics of the flow. The present work is originally motivated by a laboratory experiment and can be useful for an astrophysical problem.
Explicit wave action conservation for water waves on vertically sheared flows
NASA Astrophysics Data System (ADS)
Quinn, Brenda; Toledo, Yaron; Shrira, Victor
2016-04-01
Water waves almost always propagate on currents with a vertical structure such as currents directed towards the beach accompanied by an under-current directed back toward the deep sea or wind-induced currents which change magnitude with depth due to viscosity effects. On larger scales they also change their direction due to the Coriolis force as described by the Ekman spiral. This implies that the existing wave models, which assume vertically-averaged currents, is an approximation which is far from realistic. In recent years, ocean circulation models have significantly improved with the capability to model vertically-sheared current profiles in contrast with the earlier vertically-averaged current profiles. Further advancements have coupled wave action models to circulation models to relate the mutual effects between the two types of motion. Restricting wave models to vertically-averaged non-turbulent current profiles is obviously problematic in these cases and the primary goal of this work is to derive and examine a general wave action equation which accounts for these shortcoming. The formulation of the wave action conservation equation is made explicit by following the work of Voronovich (1976) and using known asymptotic solutions of the boundary value problem which exploit the smallness of the current magnitude compared to the wave phase velocity and/or its vertical shear and curvature. The adopted approximations are shown to be sufficient for most of the conceivable applications. This provides correction terms to the group velocity and wave action definition accounting for the shear effects, which are fitting for application to operational wave models. In the limit of vanishing current shear, the new formulation reduces to the commonly used Bretherton & Garrett (1968) no-shear wave action equation where the invariant is calculated with the current magnitude taken at the free surface. It is shown that in realistic oceanic conditions, the neglect of the vertical
Berryman, J G
2004-02-24
Layered earth models are well justified by experience, and provide a simple means of studying fairly general behavior of the elastic and poroelastic characteristics of seismic waves in the earth. Thomsen's anisotropy parameters for weak elastic and poroelastic anisotropy are now commonly used in exploration, and can be conveniently expressed in terms of the layer averages of Backus. Since our main interest is usually in the fluids underground, it would be helpful to have a set of general equations relating the Thomsen parameters as directly as possible to the fluid properties. This end can be achieved in a rather straightforward fashion for these layered earth models, and the present paper develops and then discusses these relations. Furthermore, it is found that, although there are five effective shear moduli for any layered VTI medium, one and only one effective shear modulus for the layered system contains all the dependence of pore fluids on the elastic or poroelastic constants that can be observed in vertically polarized shear waves in VTI media. The effects of the pore fluids on this effective shear modulus can be substantial - an increase of shear wave speed on the order of 10% is shown to be possible when circumstances are favorable -when the medium behaves in an undrained fashion, and the shear modulus fluctuations are large (resulting in strong anisotropy). These effects are expected to be seen at higher frequencies such as sonic and ultrasonic waves for well-logging or laboratory experiments, or at seismic wave frequencies for low permeability regions of reservoirs, prior to hydrofracing. Results presented are strictly for velocity analysis.
The Gaussian Shear Wave in a Dispersive Medium
Parker, Kevin J.; Baddour, Natalie
2014-01-01
Within the field of “imaging the biomechanical properties of tissues,” a number of approaches analyze shear wave propagation initiated by a short radiation force push. Unfortunately, it is experimentally observed that the displacement vs. time curves in lossy tissues are rapidly damped and distorted in ways that confound any simple tracking approach. This paper addresses the propagation, decay, and distortion of pulses in lossy and dispersive media, in order to derive closed form analytic expressions for the propagating pulses. The theory identifies key terms that drive the distortion and broadening of the pulse. Furthermore, the approach taken is not dependent on any particular viscoelastic model of tissue, but instead takes a general first order approach to dispersion. Examples with a Gaussian beam pattern and realistic dispersion parameters are given along with general guidelines for identifying the features of the distorting wave that are the most compact. PMID:24412170
Monitoring the lesion formation during histotripsy treatment using shear wave imaging
NASA Astrophysics Data System (ADS)
Arnal, Bastien; Lee, Wei-Ning; Pernot, Mathieu; Fink, Mathias; Tanter, Mickael
2012-11-01
Monitoring the lesion formation induced by histotripsy has mainly relied on the quantitative change in backscatter intensity using ultrasound B-mode imaging. However, how the mechanical properties of the histotripsy-treated tissue region alter during the procedure is yet to be fully investigated. We thus proposed here to monitor such a therapeutic process based on shear modulus estimated by shear wave imaging (SWI). In the therapeutic procedure, a single-element piezo-composite focused transducer (Imasonic, Besançon, France) with a center frequency of 660 kHz, a focal length of 45 mm, and an fnumber of 1 was driven by a function generator (AFG 3101, Tektronix, Beaverton, OR) and a gated RF power amplifier (GA-2500A, RITEC Inc., USA) to generate ultrasound histotripsy pulses. Histotripsy pulses were delivered for 20 seconds and then followed by a 30-second pause and a rapid monitoring step. Such a treatment and monitoring scheme was repeated for 10 mins. Both the reference measurement and monitoring were realized by SWI, where plane shear waves were generated by an 8 MHz linear array probe connected to a prototype ultrasound scanner, and acquired at a frame rate of 10000 Hz. Shear modulus was estimated and mapped in 2D through a time-of-flight algorithm. Gelatin (8%)-agar (2%) phantoms and ex-vivo porcine liver samples were tested. Regions of interests (ROI's) of 2 mm-by-2 mm in both untreated and treated regions were selected to compute the contrast-to-noise ratio (CNR). In all three scenarios where different PD's and PRF's were implemented, during the first 100 seconds of the treatment, 50% decrease in the shear modulus within the histotripsy-targeted zone was already observed, and the CNR of the shear modulus increased by 18 dB. In contrast, the backscatter intensity began to reduce and the corresponding CNR was found to increase by 6 dB only after 120 seconds of treatment. The results demonstrated that SWI can map quantitatively the change of mechanical
Solitary waves of permanent form in a deep fluid with weak shear
NASA Astrophysics Data System (ADS)
Derzho, Oleg G.; Velarde, Manuel G.
1995-06-01
The Benjamin-Davis-Acrivos-Ono equation is generalized to account for finite, large amplitude solitary waves in a sheared deep fluid. It is shown how fine structure of stratification and weak noncritical shear in such geophysical flows do affect length (shape), wave (phase) velocity, and even stability of finite amplitude solitary waves.
Shear Wave Splitting Observations Beneath Uturuncu Volcano, Bolivia
NASA Astrophysics Data System (ADS)
Sims, N. E.; Christensen, D. H.; Moore-Driskell, M. M.
2015-12-01
Anisotropy in the upper mantle is often associated with mantle flow direction through the lattice preferred orientation of anisotropic minerals such as olivine in the upper mantle material. The flow of the mantle around subduction zones can be particularly complex, and thus difficult to explain. Because of its relationship to anisotropy, analysis of shear wave splitting measurements can help to answer questions regarding the upper mantle flow that surrounds subducting slabs. Here we present SK(K)S shear wave splitting measurements from a temporary broadband network (PLUTONS) of 33 stations deployed from April 2009 to October 2012 on the Altiplano plateau around Uturuncu volcano in Bolivia. The stations are spaced 10-20 km apart, providing a high spatial resolution of the region of the mantle directly below Uturuncu volcano. Despite the lack of numerous splitting results to analyze, preliminary measurements indicate a relatively consistent pattern of fast-polarization directions in a NW-SE orientation of about N80ºW. We think that it is likely that these observations come from anisotropy in the mantle wedge above the subducting Nazca plate indicating a direction of flow in the mantle wedge that is sub-parallel to the subduction direction of the Nazca plate. Although W-E flow beneath the subducting Nazca plate cannot be completely ruled out, these results appear to be consistent with the simple model of two-dimensional corner flow in the mantle wedge and slab-entrained mantle flow beneath the slab.
Predicting prognostic factors of breast cancer using shear wave elastography.
Choi, Woo Jung; Kim, Hak Hee; Cha, Joo Hee; Shin, Hee Jung; Kim, Hyunji; Chae, Eun Young; Hong, Min Ji
2014-02-01
The purpose of the study described here was to investigate the correlation between histologic factors, including immunohistochemical factors, related to the prognosis of breast cancer and shear wave elastography (SWE) measurements. One hundred twenty-two breast cancers from 116 women were subjected to sonoelastography. Of the SWE features, mean and maximum elasticity and SWE ratio were extracted. The SWE ratio was calculated as the ratio of the stiffness of a portion of the lesion to that of a similar region of interest in fatty tissue. High ratios indicate stiffer lesions. The Mann-Whitney U-test, Kruskal-Wallis test and receiver operating characteristic (ROC) curve were used for statistical analysis. Estrogen receptor negativity, progesterone receptor negativity, p53 positivity, Ki-67 positivity, high nuclear grade, high histologic grade and large tumor (invasive) size were associated with a significantly high SWE ratio (p < 0.05). ROC curve analysis yielded SWE ratio cutoff values of 2.74-3.69 for significant immunohistochemical factors and 4.21 for the basal-like subtype by maximizing specificity while ensuring more than 80% sensitivity. Breast cancers with aggressive histologic features had high SWE ratios. Shear wave elastography may provide useful information for determining prognosis. PMID:24268451
2-D modeling of laterally acoustically coupled thin film bulk acoustic wave resonator filters.
Pensala, Tuomas; Meltaus, Johanna; Kokkonen, Kimmo; Ylilammi, Markku
2010-11-01
A 2-D model is developed for calculating lateral acoustical coupling between adjacent thin film BAW resonators forming an electrical N-port. The model is based on solution and superposition of lateral eigenmodes and eigenfrequencies in a structure consisting of adjacent regions with known plate wave dispersion properties. Mechanical and electrical response of the device are calculated as a superposition of eigenmodes according to voltage drive at one electrical port at a time while extracting current induced in the other ports, leading to a full Y-parameter description of the device. Exemplary cases are simulated to show the usefulness of the model in the study of the basic design rules of laterally coupled thin film BAW resonator filters. Model predictions are compared to an experimental 1.9-GHz band-pass filter based on aluminum nitride thin film technology and lateral acoustical coupling. Good agreement is obtained in prediction of passband behavior. The eigenmode-based model forms a useful tool for fast simulation of laterally coupled acoustic devices. It allows one to gain insight into basic device physics in a very intuitive fashion compared with more detailed but heavier finite element method. Shortcomings of this model and possible improvements are discussed. PMID:21041141
GPU computing with OpenCL to model 2D elastic wave propagation: exploring memory usage
NASA Astrophysics Data System (ADS)
Iturrarán-Viveros, Ursula; Molero-Armenta, Miguel
2015-01-01
Graphics processing units (GPUs) have become increasingly powerful in recent years. Programs exploring the advantages of this architecture could achieve large performance gains and this is the aim of new initiatives in high performance computing. The objective of this work is to develop an efficient tool to model 2D elastic wave propagation on parallel computing devices. To this end, we implement the elastodynamic finite integration technique, using the industry open standard open computing language (OpenCL) for cross-platform, parallel programming of modern processors, and an open-source toolkit called [Py]OpenCL. The code written with [Py]OpenCL can run on a wide variety of platforms; it can be used on AMD or NVIDIA GPUs as well as classical multicore CPUs, adapting to the underlying architecture. Our main contribution is its implementation with local and global memory and the performance analysis using five different computing devices (including Kepler, one of the fastest and most efficient high performance computing technologies) with various operating systems.
Wentland, Andrew L.; Wieben, Oliver; Shanmuganayagam, Dhanansayan; Krueger, Christian G.; Meudt, Jennifer J.; Consigny, Daniel; Rivera, Leonardo; McBride, Patrick E.; Reed, Jess D.; Grist, Thomas M.
2014-01-01
PURPOSE To assess measurements of pulse wave velocity (PWV) and wall shear stress (WSS) in a swine model of atherosclerosis. MATERIALS AND METHODS Nine familial hypercholesterolemic (FH) swine with angioplasty balloon catheter-induced atherosclerotic lesions to the abdominal aorta (injured group) and ten uninjured FH swine were evaluated with a 4D phase contrast (PC) MRI acquisition, as well as with radial and Cartesian 2D PC acquisitions, on a 3T MR scanner. PWV values were computed from the 2D and 4D PC techniques, compared between the injured and uninjured swine, and were validated against reference standard pressure probe-based PWV measurements. WSS values were also computed from the 4D PC MRI technique and compared between injured and uninjured groups. RESULTS PWV values were significantly greater in the injured than in the uninjured groups with the 4D PC MRI technique (p=0.03) and pressure probes (p=0.02). No significant differences were found in PWV between groups using the 2D PC techniques (p=0.75–0.83). No significant differences were found for WSS values between the injured and uninjured groups. CONCLUSION The 4D PC MRI technique provides a promising means of evaluating PWV and WSS in a swine model of atherosclerosis, providing a potential platform for developing the technique for the early detection of atherosclerosis. PMID:24964097
Near surface shear wave velocity in Bucharest, Romania
NASA Astrophysics Data System (ADS)
von Steht, M.; Jaskolla, B.; Ritter, J. R. R.
2008-12-01
Bucharest, the capital of Romania with nearly 2 1/2 million inhabitants, is endangered by the strong earthquakes in the Vrancea seismic zone. To obtain information on the near surface shear-wave velocity Vs structure and to improve the available microzonations we conducted seismic refraction measurements in two parks of the city. There the shallow Vs structure is determined along five profiles, and the compressional-wave velocity (Vp) structure is obtained along one profile. Although the amount of data collected is limited, they offer a reasonable idea about the seismic velocity distribution in these two locations. This knowledge is useful for a city like Bucharest where seismic velocity information so far is sparse and poorly documented. Using sledge-hammer blows on a steel plate and a 24-channel recording unit, we observe clear shear-wave arrivals in a very noisy environment up to a distance of 300 m from the source. The Vp model along profile 1 can be correlated with the known near surface sedimentary layers. Vp increases from 320 m/s near the surface to 1280 m/s above 55 65 m depth. The Vs models along all five profiles are characterized by low Vs (<350 m/s) in the upper 60 m depth and a maximum Vs of about 1000 m/s below this depth. In the upper 30 m the average Vs30 varies from 210 m/s to 290 m/s. The Vp-Vs relations lead to a high Poisson's ratio of 0.45 0.49 in the upper ~60 m depth, which is an indication for water-saturated clayey sediments. Such ground conditions may severely influence the ground motion during strong Vrancea earthquakes.
Wave blocking phenomenon of surface waves on a shear flow with a constant vorticity
NASA Astrophysics Data System (ADS)
Maïssa, Philippe; Rousseaux, Germain; Stepanyants, Yury
2016-03-01
Propagation of gravity-capillary surface waves on a background shear flow with a constant vorticity is studied and compared with the case when the background flow is uniform in depth. Under the assumption that the background flow gradually varies in the horizontal direction, the primary attention is paid to the wave blocking phenomenon; the effect of vorticity on this phenomenon is studied in detail. The conditions for wave blocking are obtained and categorized for different values of the governing dimensionless parameters: Froude number, dimensionless vorticity, and surface tension.
Ouared, Abderrahmane; Montagnon, Emmanuel; Cloutier, Guy
2015-10-21
A method based on adaptive torsional shear waves (ATSW) is proposed to overcome the strong attenuation of shear waves generated by a radiation force in dynamic elastography. During the inward propagation of ATSW, the magnitude of displacements is enhanced due to the convergence of shear waves and constructive interferences. The proposed method consists in generating ATSW fields from the combination of quasi-plane shear wavefronts by considering a linear superposition of displacement maps. Adaptive torsional shear waves were experimentally generated in homogeneous and heterogeneous tissue mimicking phantoms, and compared to quasi-plane shear wave propagations. Results demonstrated that displacement magnitudes by ATSW could be up to 3 times higher than those obtained with quasi-plane shear waves, that the variability of shear wave speeds was reduced, and that the signal-to-noise ratio of displacements was improved. It was also observed that ATSW could cause mechanical inclusions to resonate in heterogeneous phantoms, which further increased the displacement contrast between the inclusion and the surrounding medium. This method opens a way for the development of new noninvasive tissue characterization strategies based on ATSW in the framework of our previously reported shear wave induced resonance elastography (SWIRE) method proposed for breast cancer diagnosis. PMID:26439616
New Measurements of Shear-wave Splitting in Saudi Arabia
NASA Astrophysics Data System (ADS)
Chen, S.; Mooney, W. D.; Suzuki, J.; Zahran, H. M.; El-Hadidy, S. Y.
2015-12-01
The Saudi Geological Survey (SGS) operates a nationwide digital seismic network with more than 160 broadband seismometers that transmit to a central location at the SGS. These seismic data have been used to measure shear-wave splitting in infer anisotropy within and beneath the Arabian plate. We selected for analysis more than 300 teleseismic recordings between January, 2008 and February, 2015. Individual seismometers located on the crystalline rock of the Arabian shield provide 20 to 30 shear-wave splitting results, whereas seismometers located on volcanic rocks provide 2 to 14 reliable measurements. Here we summarize results obtained from the Tertiary volcanic fields ("harrats") of western Saudi Arabia, in particular Harrat Lunayyir and Harrat Rahat. Both of these volcanic fields have been active in historic times. Eighteen seismic stations with an average inter-station spacing of 10 km are located within Harrat Lunayyir. Seismic stations there have consistent shear-wave splitting directions ranging from N2°E to N20°W and delay times from 0.7 s to 1.6 s. This volcanic field is of particular interest because in 2009 it experienced abundant seismic activity and measureable crustal deformation that was associated with a dike intrusion into the upper crust (Pallister et al., 2010, Nature Geoscience). However, our analysis does not reveal any anomalous splitting results beneath this harrat. Fifteen seismic stations with an average inter-station spacing of 30 km are located in or adjacent to Harrat Rahat. These show very similar splitting directions to Harrat Lunayyir, ranging from N1°W to N16°W, with delay times of 1.0 s to 1.4 s. Following previous studies, we assume that these delay times are dominantly due to mantle anisotropy, with crustal anisotropy being secondary. Our results indicate a highly uniform fast-direction of anisotropy oriented approximately N10°W beneath these two volcanic fields. The measured orientation is inconsistent with the N40
Teleseismic shear wave tomography of the Japan subduction zone
NASA Astrophysics Data System (ADS)
Asamori, Koichi; Zhao, Dapeng
2015-12-01
We present a high-resolution shear wave tomography of the Japan subduction zone down to a depth of 700 km, which is determined by inverting a large number of high-quality S-wave arrival-time data from local earthquakes and teleseismic events. The subducting Pacific and Philippine Sea (PHS) slabs are revealed clearly as high-velocity (high-V) zones, whereas low-velocity (low-V) anomalies are revealed in the mantle wedge above the two slabs. The PHS slab has subducted aseismically down to a depth of 480 km under the Japan Sea and to a depth of 540 km under the Tsushima Strait. A window is revealed within the aseismic PHS slab, being consistent with P-wave tomography. Prominent low-V and high-Poisson's ratio (σ) anomalies exist below the PHS slab and above the Pacific slab, which reflect hot and wet mantle upwelling caused by the joint effect of deep dehydration of the Pacific slab and convective circulation process in the mantle wedge above the Pacific slab. The hot and wet mantle upwelling has caused the complex geometry and structure of the PHS slab in SW Japan, and contributed to the Quaternary volcanism along the Japan Sea coast. In eastern Japan, low-V zones are revealed at depths of 200-700 km below the Pacific slab, which may reflect hot upwelling from the lower mantle or even the core-mantle boundary.
Three-Dimensional Shear Wave Velocity Structure of the Peru Flat Slab Subduction Segment
NASA Astrophysics Data System (ADS)
Knezevic Antonijevic, S.; Wagner, L. S.; Beck, S. L.; Zandt, G.; Long, M. D.
2012-12-01
Recent studies focused on flat slab subduction segments in central Chile (L. S. Wagner, 2006) and Alaska (B. R. Hacker and G. A. Aber, 2012) suggest significant differences in seismic velocity structures, and hence, composition in the mantle wedge between flat and normal "steep" subducting slabs. Instead of finding the low velocities and high Vp/Vs ratios common in normal subduction zones, these studies find low Vp, high Vs, and very low Vp/Vs above flat slabs. This may indicate the presence of dry, cold material in the mantle wedge. In order to investigate the seismic velocities of the upper mantle above the Peruvian flat segment, we have inverted for 2D Rayleigh wave phase velocity maps using data from the currently deployed 40 station PULSE seismic network and some adjacent stations from the CAUGHT seismic network. We then used the sensitivity of surface waves to shear wave velocity structure with depth to develop a 3D shear wave velocity model. This model will allow us to determine the nature of the mantle lithosphere above the flat slab, and how this may have influenced the development of local topography. For example, dry conditions (high Vs velocities) above the flat slab would imply greater strength of this material, possibly making it capable of causing further inland overthrusting, while wet conditions (low Vs) would imply weaker material. This could provide some insight into the ongoing debate over whether the Fitzcarrald arch (along the northern most flank of the Altiplano) could be a topographical response to the subducted Nazca ridge hundred kilometers away from the trench (N. Espurt, 2012, P. Baby, 2005, V. A. Ramos, 2012) or not (J. Martinod, 2005, M. Wipf, 2008, T. Gerya, 2008).
Transition from 1D to 2D Laser-Induced Ultrasonic Wave Propagation in an Extended Plate
NASA Astrophysics Data System (ADS)
Laloš, Jernej; Požar, Tomaž; Možina, Janez
2016-05-01
Optodynamic interaction between a laser pulse and the surface of an opaque, solid elastic object produces transient waves that propagate and reverberate within the object. They can be, in general, categorized into three distinctive types which are all formed through different mechanisms: ablation-induced waves, light-pressure-induced waves, and thermoelastic waves. In this paper, out-of-plane displacements of such waves are simulated at the epicentral position on the opposite side of an extended plane-parallel elastic plate. Wave propagation is mathematically described by Green's transfer functions convolved with suitable time profiles of the incoming laser pulses. The simulated size of the circularly symmetric laser-illuminated area on the plate surface is varied to show the limit-to-limit transition of the displacement waveforms: from a 2D point source to an infinite 1D source.
Shear Wave Splitting Analysis of Aftershocks of the 2013 Mw6.6 Lushan Earthquake, China
NASA Astrophysics Data System (ADS)
Liu, Y.; Zhang, H.
2013-12-01
Shear wave splits into faster and slower shear waves that are nearly perpendicular when it travels through an anisotropic medium. There are two important parameters of shear wave splitting, one is the fast polarization direction of the fast shear wave and the other one is the time delay of the slow shear wave. The mechanisms for anisotropy in the upper crust can be divided into two categories. The first category is stress-induced anisotropy related to alignment of cracks in response to the in situ stress field. The second category is structural anisotropy associated with aligned planar features such as fault zone fabrics, sedimentary bedding planes and aligned minerals. We can characterize anisotropy around fault zone by shear wave splitting analysis. We used cross-correlation method for the shear wave splitting analysis. Since the faster shear wave and the slower shear wave are from the same source, they will correlate well after the time delay correction. We rotated two horizontal seismograms at a 10 increment of azimuth α from 00 to 1800. For each α, the cross-correlation coefficients between the two orthogonal seismograms are calculated for a range of time delays τ. When the absolute value of cross-correlation coefficient reaches a maximum, the corresponding values of α and τ are chosen as the fast polarization direction of the faster shear wave and the time delay of the slower shear wave, respectively. We chose 200 aftershocks observed at a temporary array consisting of 29 stations in the Lushan region. Shear wave arrivals were first picked for setting up the time window for the shear wave splitting analysis using the cross-correlation method. Because these 200 events are shallower than 20km, we can infer that the shear wave splitting is caused by crustal anisotropy. The rose diagram of the fast polarization directions of the fast shear waves showed two major directions. One is nearly parallel to the south-north trending fault system in this region, and
Razani, Marjan; Mariampillai, Adrian; Sun, Cuiru; Luk, Timothy W. H.; Yang, Victor X. D.; Kolios, Michael C.
2012-01-01
In this work, we explored the potential of measuring shear wave propagation using optical coherence elastography (OCE) based on a swept-source optical coherence tomography (OCT) system. Shear waves were generated using a 20 MHz piezoelectric transducer (circular element 8.5 mm diameter) transmitting sine-wave bursts of 400 μs, synchronized with the OCT swept source wavelength sweep. The acoustic radiation force (ARF) was applied to two gelatin phantoms (differing in gelatin concentration by weight, 8% vs. 14%). Differential OCT phase maps, measured with and without the ARF, demonstrate microscopic displacement generated by shear wave propagation in these phantoms of different stiffness. We present preliminary results of OCT derived shear wave propagation velocity and modulus, and compare these results to rheometer measurements. The results demonstrate the feasibility of shear wave OCE (SW-OCE) for high-resolution microscopic homogeneous tissue mechanical property characterization. PMID:22567590
Razani, Marjan; Mariampillai, Adrian; Sun, Cuiru; Luk, Timothy W H; Yang, Victor X D; Kolios, Michael C
2012-05-01
In this work, we explored the potential of measuring shear wave propagation using optical coherence elastography (OCE) based on a swept-source optical coherence tomography (OCT) system. Shear waves were generated using a 20 MHz piezoelectric transducer (circular element 8.5 mm diameter) transmitting sine-wave bursts of 400 μs, synchronized with the OCT swept source wavelength sweep. The acoustic radiation force (ARF) was applied to two gelatin phantoms (differing in gelatin concentration by weight, 8% vs. 14%). Differential OCT phase maps, measured with and without the ARF, demonstrate microscopic displacement generated by shear wave propagation in these phantoms of different stiffness. We present preliminary results of OCT derived shear wave propagation velocity and modulus, and compare these results to rheometer measurements. The results demonstrate the feasibility of shear wave OCE (SW-OCE) for high-resolution microscopic homogeneous tissue mechanical property characterization. PMID:22567590
Drift-wave transport in the velocity shear layer
NASA Astrophysics Data System (ADS)
Rosalem, K. C.; Roberto, M.; Caldas, I. L.
2016-07-01
Particle drift driven by electrostatic wave fluctuations is numerically computed to describe the transport in a gradient velocity layer at the tokamak plasma edge. We consider an equilibrium plasma in large aspect ratio approximation with E × B flow and specified toroidal plasma velocity, electric field, and magnetic field profiles. A symplectic map, previously derived for infinite coherent time modes, is used to describe the transport dependence on the electric, magnetic, and plasma velocity shears. We also show that resonant perturbations and their correspondent islands in the Poincaré maps are much affected by the toroidal velocity profiles. Moreover, shearless transport barriers, identified by extremum values of the perturbed rotation number profiles of the invariant curves, allow chaotic trajectories trapped into the plasma. We investigate the influence of the toroidal plasma velocity profile on these shearless transport barriers.
Shear wave velocities of unconsolidated shallow sediments in the Gulf of Mexico
Lee, Myung W.
2013-01-01
Accurate shear-wave velocities for shallow sediments are important for a variety of seismic applications such as inver-sion and amplitude versus offset analysis. During the U.S. Department of Energy-sponsored Gas Hydrate Joint Industry Project Leg II, shear-wave velocities were measured at six wells in the Gulf of Mexico using the logging-while-drilling SonicScope acoustic tool. Because the tool measurement point was only 35 feet from the drill bit, the adverse effect of the borehole condition, which is severe for the shallow unconsolidated sediments in the Gulf of Mexico, was mini-mized and accurate shear-wave velocities of unconsolidated sediments were measured. Measured shear-wave velocities were compared with the shear-wave velocities predicted from the compressional-wave velocities using empirical formulas and the rock physics models based on the Biot-Gassmann theory, and the effectiveness of the two prediction methods was evaluated. Although the empirical equation derived from measured shear-wave data is accurate for predicting shear-wave velocities for depths greater than 500 feet in these wells, the three-phase Biot-Gassmann-theory -based theory appears to be optimum for predicting shear-wave velocities for shallow unconsolidated sediments in the Gulf of Mexico.
Luo, Xin; Lu, Xin; Cong, Chunxiao; Yu, Ting; Xiong, Qihua; Ying Quek, Su
2015-01-01
2D layered materials have recently attracted tremendous interest due to their fascinating properties and potential applications. The interlayer interactions are much weaker than the intralayer bonds, allowing the as-synthesized materials to exhibit different stacking sequences, leading to different physical properties. Here, we show that regardless of the space group of the 2D materials, the Raman frequencies of the interlayer shear modes observed under the typical configuration blue shift for AB stacked materials, and red shift for ABC stacked materials, as the number of layers increases. Our predictions are made using an intuitive bond polarizability model which shows that stacking sequence plays a key role in determining which interlayer shear modes lead to the largest change in polarizability (Raman intensity); the modes with the largest Raman intensity determining the frequency trends. We present direct evidence for these conclusions by studying the Raman modes in few layer graphene, MoS2, MoSe2, WSe2 and Bi2Se3, using both first principles calculations and Raman spectroscopy. This study sheds light on the influence of stacking sequence on the Raman intensities of intrinsic interlayer modes in 2D layered materials in general, and leads to a practical way of identifying the stacking sequence in these materials. PMID:26469313
Luo, Xin; Lu, Xin; Cong, Chunxiao; Yu, Ting; Xiong, Qihua; Quek, Su Ying
2015-01-01
2D layered materials have recently attracted tremendous interest due to their fascinating properties and potential applications. The interlayer interactions are much weaker than the intralayer bonds, allowing the as-synthesized materials to exhibit different stacking sequences, leading to different physical properties. Here, we show that regardless of the space group of the 2D materials, the Raman frequencies of the interlayer shear modes observed under the typical z(xx)z configuration blue shift for AB stacked materials, and red shift for ABC stacked materials, as the number of layers increases. Our predictions are made using an intuitive bond polarizability model which shows that stacking sequence plays a key role in determining which interlayer shear modes lead to the largest change in polarizability (Raman intensity); the modes with the largest Raman intensity determining the frequency trends. We present direct evidence for these conclusions by studying the Raman modes in few layer graphene, MoS2, MoSe2, WSe2 and Bi2Se3, using both first principles calculations and Raman spectroscopy. This study sheds light on the influence of stacking sequence on the Raman intensities of intrinsic interlayer modes in 2D layered materials in general, and leads to a practical way of identifying the stacking sequence in these materials. PMID:26469313
[INVITED] Laser generation and detection of ultrafast shear acoustic waves in solids and liquids
NASA Astrophysics Data System (ADS)
Pezeril, Thomas
2016-09-01
The aim of this article is to provide an overview of the up-to-date findings related to ultrafast shear acoustic waves. Recent progress obtained for the laser generation and detection of picosecond shear acoustic waves in solids and liquids is reviewed. Examples in which the transverse isotropic symmetry of the sample structure is broken in order to permit shear acoustic wave generation through sudden laser heating are described in detail. Alternative photo-induced mechanisms for ultrafast shear acoustic generation in metals, semiconductors, insulators, magnetostrictive, piezoelectric and electrostrictive materials are reviewed as well. With reference to key experiments, an all-optical technique employed to probe longitudinal and shear structural dynamics in the GHz frequency range in ultra-thin liquid films is described. This technique, based on specific ultrafast shear acoustic transducers, has opened new perspectives that will be discussed for ultrafast shear acoustic probing of viscoelastic liquids at the nanometer scale.
Analysis shear wave velocity structure obtained from surface wave methods in Bornova, Izmir
NASA Astrophysics Data System (ADS)
Pamuk, Eren; Özdaǧ, Özkan Cevdet; Akgün, Mustafa
2016-04-01
Properties of the soil from the bedrock is necessary to describe accurately and reliably for the reduction of earthquake damage. Because seismic waves change their amplitude and frequency content owing to acoustic impedance difference between soil and bedrock. Firstly, shear wave velocity and depth information of layers on bedrock is needed to detect this changing. Shear wave velocity can be obtained using inversion of Rayleigh wave dispersion curves obtained from surface wave methods (MASW- the Multichannel Analysis of Surface Waves, ReMi-Refraction Microtremor, SPAC-Spatial Autocorrelation). While research depth is limeted in active source study, a passive source methods are utilized for deep depth which is not reached using active source methods. ReMi method is used to determine layer thickness and velocity up to 100 m using seismic refraction measurement systems.The research carried out up to desired depth depending on radius using SPAC which is utilized easily in conditions that district using of seismic studies in the city. Vs profiles which are required to calculate deformations in under static and dynamic loads can be obtained with high resolution using combining rayleigh wave dispersion curve obtained from active and passive source methods. In the this study, Surface waves data were collected using the measurements of MASW, ReMi and SPAC at the İzmir Bornova region. Dispersion curves obtained from surface wave methods were combined in wide frequency band and Vs-depth profiles were obtained using inversion. Reliability of the resulting soil profiles were provided by comparison with theoretical transfer function obtained from soil paremeters and observed soil transfer function from Nakamura technique and by examination of fitting between these functions. Vs values are changed between 200-830 m/s and engineering bedrock (Vs>760 m/s) depth is approximately 150 m.
Shear wave velocities from noise correlation at local scale
De Nisco, G.; Nunziata, C.; Vaccari, F.; Panza, G. F.
2008-07-08
Cross correlations of ambient seismic noise recordings have been studied to infer shear seismic velocities with depth. Experiments have been done in the crowded and noisy historical centre of Napoli over inter-station distances from 50 m to about 400 m, whereas active seismic spreadings are prohibitive, even for just one receiver. Group velocity dispersion curves have been extracted with FTAN method from the noise cross correlations and then the non linear inversion of them has resulted in Vs profiles with depth. The information of near by stratigraphies and the range of Vs variability for samples of Neapolitan soils and rocks confirms the validity of results obtained with our expeditious procedure. Moreover, the good comparison of noise H/V frequency of the first main peak with 1D and 2D spectral amplifications encourages to continue experiments of noise cross-correlation. If confirmed in other geological settings, the proposed approach could reveal a low cost methodology to obtain reliable and detailed Vs velocity profiles.
Passive elastography: shear-wave tomography from physiological-noise correlation in soft tissues.
Gallot, Thomas; Catheline, Stefan; Roux, Philippe; Brum, Javier; Benech, Nicolas; Negreira, Carlos
2011-06-01
Inspired by seismic-noise correlation and time reversal, a shear-wave tomography of soft tissues using an ultrafast ultrasonic scanner is presented here. Free from the need for controlled shear-wave sources, this passive elastography is based on Green's function retrieval and takes advantage of the permanent physiological noise of the human body. PMID:21693392
Modelling the impulse diffraction field of shear waves in transverse isotropic viscoelastic medium.
Chatelin, Simon; Gennisson, Jean-Luc; Bernal, Miguel; Tanter, Mickael; Pernot, Mathieu
2015-05-01
The generation of shear waves from an ultrasound focused beam has been developed as a major concept for remote palpation using shear wave elastography (SWE). For muscular diagnostic applications, characteristics of the shear wave profile will strongly depend on characteristics of the transducer as well as the orientation of muscular fibers and the tissue viscoelastic properties. The numerical simulation of shear waves generated from a specific probe in an anisotropic viscoelastic medium is a key issue for further developments of SWE in fibrous soft tissues. In this study we propose a complete numerical tool allowing 3D simulation of a shear wave front in anisotropic viscoelastic media. From the description of an ultrasonic transducer, the shear wave source is simulated by using Field's II software and shear wave propagation described by using the Green's formalism. Finally, the comparison between simulations and experiments are successively performed for both shear wave velocity and dispersion profile in a transverse isotropic hydrogel phantom, in vivo forearm muscle and in vivo biceps brachii. PMID:25880794
Modelling the impulse diffraction field of shear waves in transverse isotropic viscoelastic medium
NASA Astrophysics Data System (ADS)
Chatelin, Simon; Gennisson, Jean-Luc; Bernal, Miguel; Tanter, Mickael; Pernot, Mathieu
2015-05-01
The generation of shear waves from an ultrasound focused beam has been developed as a major concept for remote palpation using shear wave elastography (SWE). For muscular diagnostic applications, characteristics of the shear wave profile will strongly depend on characteristics of the transducer as well as the orientation of muscular fibers and the tissue viscoelastic properties. The numerical simulation of shear waves generated from a specific probe in an anisotropic viscoelastic medium is a key issue for further developments of SWE in fibrous soft tissues. In this study we propose a complete numerical tool allowing 3D simulation of a shear wave front in anisotropic viscoelastic media. From the description of an ultrasonic transducer, the shear wave source is simulated by using Field’s II software and shear wave propagation described by using the Green’s formalism. Finally, the comparison between simulations and experiments are successively performed for both shear wave velocity and dispersion profile in a transverse isotropic hydrogel phantom, in vivo forearm muscle and in vivo biceps brachii.
NASA Astrophysics Data System (ADS)
Kimura, Masao
2005-11-01
Shear-wave velocity is one of the important parameters that characterize the physical properties of marine sediments. In this study, a new method is proposed for measuring shear-wave velocity in marine sediments by using radiation impedance. Shear-wave velocities for three kinds of urethane rubber with different Japanese Industrial Standards hardness values were obtained by radiation impedance and time-of-flight measurement techniques. It was shown that the values of the shear-wave velocity measured by the radiation impedance method were consistent with those of time-of-flight measurements. It was then shown that the shear-wave velocities for air- and water-saturated beach sands are different. It was also found that the indicated shear-wave velocity is dependent on the vibrating plate radius because the instrument measures an average shear-wave velocity within a depth window beneath the plate; the larger the plate radius, the deeper the averaging window. Finally, measurements were made on two-layered media in which air-saturated beach sand or urethane rubber was covered with air-saturated clay, and the relationship between the thickness of the clay layer and the indicated shear-wave velocity was investigated.
Anomalous shear wave delays and surface wave velocities at Yellowstone Caldera, Wyoming
Daniel, R.G.; Boore, D.M.
1982-04-10
To investigate the effects of a geothermal area on the propagation of intermediate-period (1--30 s) teleseismic body waves and surface waves, a specially designed portable seismograph system was operated in Yellowstone Caldera, Wyoming. Travel time residuals, relative to a station outside the caldera, of up to 2 s for compressional phases are in agreement with short-period residuals for P phases measured by other investigators. Travel time delays for shear arrivals in the intermediate-period band range from 2 to 9 s and decrease with increasing dT/d..delta... Measured Rayleigh wave phase velocities are extremely low, ranging from 3.2 km/s at 27-s period to 2.0 km/s at 7-s period; the estimated uncertainty associated with these values is 15%. We propose a model for compressional and shear velocities and Poisson's ratio beneath the Yellowstone caldera which fits the teleseismic body and surface wave data: it consists of a highly anomalous crust with an average shear velocity of 3.0 km/s overlying an upper mantle with average velocity of 4.1 km/s. The high average value of Poisson's ratio in the crust (0.34) suggests the presence of fluids there; Poisson's ratio in the mantle between 40 and approximately 200 km is more nearly normal (0.29) than in the crust. A discrepancy between normal values of Poisson's ratio in the crust calculated from short-period data and high values calculated from teleseismic data can be resolved by postulating a viscoelastic crustal model with frequency-dependent shear velocity and attenuation.
Finite-difference modelling to evaluate seismic P-wave and shear-wave field data
NASA Astrophysics Data System (ADS)
Burschil, T.; Beilecke, T.; Krawczyk, C. M.
2015-01-01
High-resolution reflection seismic methods are an established non-destructive tool for engineering tasks. In the near surface, shear-wave reflection seismic measurements usually offer a higher spatial resolution in the same effective signal frequency spectrum than P-wave data, but data quality varies more strongly. To discuss the causes of these differences, we investigated a P-wave and a SH-wave seismic reflection profile measured at the same location on the island of Föhr, Germany and applied seismic reflection processing to the field data as well as finite-difference modelling of the seismic wave field. The simulations calculated were adapted to the acquisition field geometry, comprising 2 m receiver distance (1 m for SH wave) and 4 m shot distance along the 1.5 km long P-wave and 800 m long SH-wave profiles. A Ricker wavelet and the use of absorbing frames were first-order model parameters. The petrophysical parameters to populate the structural models down to 400 m depth were taken from borehole data, VSP (vertical seismic profile) measurements and cross-plot relations. The simulation of the P-wave wave-field was based on interpretation of the P-wave depth section that included a priori information from boreholes and airborne electromagnetics. Velocities for 14 layers in the model were derived from the analysis of five nearby VSPs (vP =1600-2300 m s-1). Synthetic shot data were compared with the field data and seismic sections were created. Major features like direct wave and reflections are imaged. We reproduce the mayor reflectors in the depth section of the field data, e.g. a prominent till layer and several deep reflectors. The SH-wave model was adapted accordingly but only led to minor correlation with the field data and produced a higher signal-to-noise ratio. Therefore, we suggest to consider for future simulations additional features like intrinsic damping, thin layering, or a near-surface weathering layer. These may lead to a better understanding of
Acoustic-radiation-force-induced shear wave propagation in cardiac tissue
NASA Astrophysics Data System (ADS)
Bouchard, Richard R.; Wolf, Patrick D.; Hsu, Stephen J.; Dumont, Douglas M.; Trahey, Gregg E.
2009-02-01
Shear wave elasticity imaging (SWEI) was employed to track acoustic radiation force (ARF)-induced shear waves in the myocardium of a beating heart. Shear waves were generated in and tracked through the myocardium of the left ventricular free wall (LVFW) in an in vivo heart that was exposed through a thoracotomy; matched studies were also preformed on an ex vivo myocardial specimen. Average shear wave velocities ranged from 2.22 to 2.53 m/s for the ex vivo specimen and 1.5 to 2.9 m/s (1.5-2.09 m/s during diastole; 2.9 m/s during systole) for in vivo specimens. Despite the known rotation of myocardial fiber orientation with tissue depth, there was no statistically significant shear wave velocity depth dependence observed in any of the experimental trials.
Shear Wave Splitting Analysis to Estimate Fracture Orientation and Frequency Dependent Anisotropy
NASA Astrophysics Data System (ADS)
Gholami, Raoof; Moradzadeh, Ali; Rasouli, Vamegh; Hanachi, Javid
2016-02-01
Shear wave splitting is a well-known method for indication of orientation, radius, and length of fractures in subsurface layers. In this paper, a three component near offset VSP data acquired from a fractured sandstone reservoir in southern part of Iran was used to analyse shear wave splitting and frequency-dependent anisotropy assessment. Polarization angle obtained by performing rotation on radial and transverse components of VSP data was used to determine the direction of polarization of fast shear wave which corresponds to direction of fractures. It was shown that correct implementation of shear wave splitting analysis can be used for determination of fracture direction. During frequency-dependent anisotropy analysis, it was found that the time delays in shear-waves decrease as the frequency increases. It was clearly demonstrated throughout this study that anisotropy may have an inverse relationship with frequency. The analysis presented in this paper complements the studied conducted by other researchers in this field of research.
A shear-wave polarization study in the Wellington region New Zealand
Gledhill, K.R. )
1990-08-01
A month of digital data from two three component seismograph stations near Wellington, New Zealand, was analyzed as part of a feasibility study for a major project to investigate shear-wave splitting. Although the total number of earthquakes studies was small (14), some suggestive results were obtained. Almost all events recorded within the shear wave window showed a phase reversal of the horizontal components after one or two shear wave cycles, suggesting that there are actually two shear-wave arrivals. The measured polarization of the first shear wave arrivals was N (31 {plus minus} 11) E. This polarization alignment cannot be explained by focal mechanisms, and it is unlikely to be due to topography because of the station-to-station correlation. The present evidence suggests the most likely cause is crustal anistropy due to the geological structure at shallow depth, rather than stress aligned micro-cracks.
NASA Astrophysics Data System (ADS)
Nguyen, Thu-Mai; Song, Shaozhen; Arnal, Bastien; Wong, Emily Y.; Wang, Ruikang K.; O'Donnell, Matthew
2014-02-01
Shear wave elastography measures the stiffness of soft tissues from the speed of propagating shear waves induced in tissue. Optical coherence tomography (OCT) is a promising detection modality given its high sensitivity and spatial resolution, making it suitable for elastic characterization of skin, peripheral vasculature or ocular tissues. For clinical applications, it would be valuable to use a non-contact shear source. Thus, we propose acoustic radiation force as a remote shear source combined with OCT for visualization. A single-element focused transducer (central frequency 7.5 MHz) was used to apply a maximal pressure of ~3 MPa for 100 μs in agar phantoms. It induced shear waves with an amplitude of several hundreds of nanometers and a broadband spectrum in the kilohertz range. Phasesensitive OCT was used to track shear waves at an equivalent frame rate of 47 kHz. We reconstructed shear modulus maps in a heterogeneous phantom. In addition, we use 3-ms long coded excitation to increase the displacement signal-to-noise ratio. We applied digital pulse compression to the resulting displacement field to obtain a gain of ~15 dB compared to standard pulse excitation while maintaining the US pressure level and the shear wave spatial and temporal resolution. This is a promising result for shear wave generation at low US pressures (~ 1 MPa).
Anisotropy across the Sorgenfrei Tornquist Zone from shear wave splitting
NASA Astrophysics Data System (ADS)
Wylegalla, K.; Bock, G.; Gossler, J.; Hanka, W.; TOR Working Group
1999-12-01
During the TOR-1 passive seismic experiment in 1996/97, a maximum of 139 temporary seismograph stations were operating over the Sorgenfrei-Tornquist Zone (STZ) in an area extending from northern Germany through Denmark to central Sweden. One of the objectives was to study horizontal anisotropy directions in the subcrustal lithosphere and asthenosphere across the Trans-European Suture Zone. To achieve this goal, broad-band and intermediate-period (5 s) data of the TOR-1 stations and additional stations of permanent networks (GRSN, GEOFON) were analysed for splitting of SKS and SKKS phases. As a result of the relatively dense station spacing, the method offers good lateral resolution of anisotropy. Preliminary results suggest that the directions of the fast horizontal S wave velocity are affected by the STZ. In central Europe and southern Sweden, far away from the STZ, fast S wave directions are approximately E-W while they turn more northerly closer to the STZ where they are approximately parallel to the trend of the STZ. No significant shear wave splitting was observed north of 57°N and east of 14°E. Small delay times between 0.2 and 0.5 s observed at the northernmost TOR-1 station T40S and T60S may be controlled by anisotropy in a thickened crust. The mantle contribution of horizontal anisotropy within the STZ is probably constrained to an approximately 60-km-thick zone in the depth range between 70 and 300 km. The observations are consistent with a model where azimuthal anisotropy is not governed by present-day mantle flow in the asthenosphere, but rather is frozen into the subcrustal lithosphere during the last episode of tectonic activity.
Modelling study of challenges in sinkhole detection with shear wave reflection seismics
NASA Astrophysics Data System (ADS)
Burschil, Thomas; Krawczyk, CharLotte M.
2016-04-01
The detection of cavities with reflection seismics is a difficult task even if high impedance contrasts are assumed. Especially the shear wave reflection method with a higher resolution potential trough lower velocities and short wavelength has come into focus of investigation. But shear wave propagation fails if material exists that partially has no shear strength. The shear wave does not propagate into or through those voids. Here, we evaluate the influence of a possible fracture zone above a cavity. We simulate shear wave propagation with finite difference modelling for two reference models, with and without cavity, and various sets of input models with a fracture zone above the cavity. Reflections and multiples of the reference models image the subsidence structure and the cavity. For the fracture input models, we implemented a fracture network, derived from numerical crack propagation modelling (Schneider-Löbens et al., 2015). The cracks possess the minimum possible aperture of one grid point (i.e. 0.1 m) and no shear stiffness. The seismic modelling exhibits that the shear wave does not pass through the fracture zone and shadows the subjacent cavity. Sequences of randomly discontinuous cracks, cf. displacement discontinuity model with zero crack stiffness, approximate partially seismic connected rock on both sides of the crack. The amount of these seismic pathways determines whether a reflection of the cavity can be detected at the surface or not. Cracks with higher aperture, e.g. two or three grid points, need a higher amount of intact rock/defective cracks, since more connected grid points are necessary to create seismic pathways. Furthermore, it turns out that the crack filling is important for shear wave transmission. While a mineralized fracture zone, implemented with high velocity, facilitate shear wave propagation, water or air-filled cracks avoid shear wave transmission. Crack orientation affects the shear wave propagation through the geometry. A
Flow shear suppression of turbulence using externally driven ion Bernstein and Alfven waves
Biglari, H.; Ono, M. . Plasma Physics Lab.); Diamond, P.H. . Dept. of Physics); Craddock, G.G. )
1991-01-01
The utilization of externally-launched radio-frequency waves as a means of active confinement control through the generation of sheared poloidal flows is explored. For low-frequency waves, kinetic Alfven waves are proposed, and are shown to drive sheared E {times} B flows as a result of the radial variation in the electromagnetic Reynolds stress. In the high frequency regime, ion Bernstein waves are considered, and shown to generate sheared poloidal rotation through the ponderomotive force. In either case, it is shown that modest amounts of absorbed power ({approximately} few 100 kW) are required to suppress turbulence in a region of several cm radial width. 9 refs.
Estimation of viscoelastic parameters in Prony series from shear wave propagation
NASA Astrophysics Data System (ADS)
Jung, Jae-Wook; Hong, Jung-Wuk; Lee, Hyoung-Ki; Choi, Kiwan
2016-06-01
When acquiring accurate ultrasonic images, we must precisely estimate the mechanical properties of the soft tissue. This study investigates and estimates the viscoelastic properties of the tissue by analyzing shear waves generated through an acoustic radiation force. The shear waves are sourced from a localized pushing force acting for a certain duration, and the generated waves travel horizontally. The wave velocities depend on the mechanical properties of the tissue such as the shear modulus and viscoelastic properties; therefore, we can inversely calculate the properties of the tissue through parametric studies.
Fluid Effects on Shear for Seismic Waves in Finely Layered Porous Media
Berryman, J G
2004-07-22
Although there are five effective shear moduli for any layered VTI medium, one and only one effective shear modulus of the layered system (namely the uniaxial shear) contains all the dependence of pore fluids on the elastic or poroelastic constants that can be observed in vertically polarized shear waves. Pore fluids can increase the magnitude the shear energy stored in this modulus by an amount that ranges from the smallest to the largest effective shear moduli of the VTI system. But, since there are five shear moduli in play, the overall increase in shear energy due to fluids is reduced by a factor of about 5 in general. We can therefore give definite bounds on the maximum increase of overall shear modulus, being about 20% of the allowed range as liquid is fully substituted for gas. An attendant increase of density (depending on porosity and fluid density) by approximately 5 to 10% decreases the shear wave speed and, thereby, partially offsets the effect of this shear modulus increase. The final result is an increase of shear wave speed on the order of 5 to 10%. This increase is shown to be possible under most favorable circumstances - i.e. when the shear modulus fluctuations are large (resulting in strong anisotropy) and the medium behaves in an undrained fashion due to fluid trapping. At frequencies higher than seismic (such as sonic and ultrasonic waves for well-logging or laboratory experiments), resulting short response times also produce the requisite undrained behavior and, therefore, fluids also affect shear waves at high frequencies by increasing rigidity.
The shear wave velocity underneath Bucharest city, Romania, from the analysis of Love waves
NASA Astrophysics Data System (ADS)
Sèbe, Olivier; Forbriger, Thomas; Ritter, Joachim R. R.
2009-03-01
From the dispersion of Love waves, we infer models of shear wave velocity structure underneath Bucharest (Romania) at depths down to 2km that can contribute to seismic hazard estimation. Waves from eight regional events recorded during 10months with a network of 34 seismic broad-band stations of the URban Seismology (URS) experiment are used. Although these events provide poor azimuthal coverage the data reliably constrain a shear wave velocity model with an interface between the Neogene and the Cretaceous sediments that is dipping northwards towards the Carpathian mountains. Array processing techniques that account for non-uniform wave propagation are used to estimate the dispersion of structural phase velocity. From this, we infer subsurface structure at three different latitudes. The Neogene sediments are represented by a gradient layer with no significant lateral variation. Shear wave velocity increases from approximately 400ms-1 near the surface to 1kms-1 at 1km depth and 5km in the south, and to 1.35kms-1 at 1.5km depth and 5km in the north from the centre of Bucharest, respectively. For the half-space representing the Cretaceous sediments, we obtain shear wave velocities of 2.7-2.9kms-1. The results are consistent with results from boreholes and shallow seismics for the near-surface structure and results from receiver function studies and crustal refraction seismic studies for the deeper structure. The details of the Neogene layer comprising a vertical gradient fill a gap in existing models of the subsurface structure of Bucharest and can contribute to modelling of seismic hazard for the city. Since the signal-to-noise ratio restricted useful data to the frequency range from 90 to 290mHz, the inversion could not constrain the near-surface velocity independently. Due to strong trade-off between near-surface velocity and depth of half-space, the latter had to be introduced as a priori data from previous studies.
Ren, X; Domier, C W; Kramer, G; Luhmann, N C; Muscatello, C M; Shi, L; Tobias, B J; Valeo, E
2014-11-01
A synthetic microwave imaging reflectometer (MIR) diagnostic employing the full-wave reflectometer code (FWR2D) has been developed and is currently being used to guide the design of real systems, such as the one recently installed on DIII-D. The FWR2D code utilizes real plasma profiles as input, and it is combined with optical simulation tools for synthetic diagnostic signal generation. A detailed discussion of FWR2D and the process to generate the synthetic signal are presented in this paper. The synthetic signal is also compared to a prescribed density fluctuation spectrum to quantify the imaging quality. An example is presented with H-mode-like plasma profiles derived from a DIII-D discharge, where the MIR focal is located in the pedestal region. It is shown that MIR is suitable for diagnosing fluctuations with poloidal wavenumber up to 2.0 cm(-1) and fluctuation amplitudes less than 5%. PMID:25430276
Xia, J.; Miller, R.D.; Park, C.B.; Hunter, J.A.; Harris, J.B.; Ivanov, J.
2002-01-01
Recent field tests illustrate the accuracy and consistency of calculating near-surface shear (S)-wave velocities using multichannel analysis of surface waves (MASW). S-wave velocity profiles (S-wave velocity vs. depth) derived from MASW compared favorably to direct borehole measurements at sites in Kansas, British Columbia, and Wyoming. Effects of changing the total number of recording channels, sampling interval, source offset, and receiver spacing on the inverted S-wave velocity were studied at a test site in Lawrence, Kansas. On the average, the difference between MASW calculated Vs and borehole measured Vs in eight wells along the Fraser River in Vancouver, Canada was less than 15%. One of the eight wells was a blind test well with the calculated overall difference between MASW and borehole measurements less than 9%. No systematic differences were observed in derived Vs values from any of the eight test sites. Surface wave analysis performed on surface data from Wyoming provided S-wave velocities in near-surface materials. Velocity profiles from MASW were confirmed by measurements based on suspension log analysis. ?? 2002 Elsevier Science Ltd. All rights reserved.
Measurement of sound speed vs. depth in South Pole ice: pressure waves and shear waves
IceCube Collaboration; Klein, Spencer
2009-06-04
We have measured the speed of both pressure waves and shear waves as a function of depth between 80 and 500 m depth in South Pole ice with better than 1% precision. The measurements were made using the South Pole Acoustic Test Setup (SPATS), an array of transmitters and sensors deployed in the ice at the South Pole in order to measure the acoustic properties relevant to acoustic detection of astrophysical neutrinos. The transmitters and sensors use piezoceramics operating at {approx}5-25 kHz. Between 200 m and 500 m depth, the measured profile is consistent with zero variation of the sound speed with depth, resulting in zero refraction, for both pressure and shear waves. We also performed a complementary study featuring an explosive signal propagating vertically from 50 to 2250 m depth, from which we determined a value for the pressure wave speed consistent with that determined for shallower depths, higher frequencies, and horizontal propagation with the SPATS sensors. The sound speed profile presented here can be used to achieve good acoustic source position and emission time reconstruction in general, and neutrino direction and energy reconstruction in particular. The reconstructed quantities could also help separate neutrino signals from background.
High speed all optical shear wave imaging optical coherence elastography (Conference Presentation)
NASA Astrophysics Data System (ADS)
Song, Shaozhen; Hsieh, Bao-Yu; Wei, Wei; Shen, Tueng; O'Donnell, Matthew; Wang, Ruikang K.
2016-03-01
Optical Coherence Elastography (OCE) is a non-invasive testing modality that maps the mechanical property of soft tissues with high sensitivity and spatial resolution using phase-sensitive optical coherence tomography (PhS-OCT). Shear wave OCE (SW-OCE) is a leading technique that relies on the speed of propagating shear waves to provide a quantitative elastography. Previous shear wave imaging OCT techniques are based on repeated M-B scans, which have several drawbacks such as long acquisition time and repeated wave stimulations. Recent developments of Fourier domain mode-locked high-speed swept-source OCT system has enabled enough speed to perform KHz B-scan rate OCT imaging. Here we propose ultra-high speed, single shot shear wave imaging to capture single-shot transient shear wave propagation to perform SW-OCE. The frame rate of shear wave imaging is 16 kHz, at A-line rate of ~1.62 MHz, which allows the detection of high-frequency shear wave of up to 8 kHz. The shear wave is generated photothermal-acoustically, by ultra-violet pulsed laser, which requires no contact to OCE subjects, while launching high frequency shear waves that carries rich localized elasticity information. The image acquisition and processing can be performed at video-rate, which enables real-time 3D elastography. SW-OCE measurements are demonstrated on tissue-mimicking phantoms and porcine ocular tissue. This approach opens up the feasibility to perform real-time 3D SW-OCE in clinical applications, to obtain high-resolution localized quantitative measurement of tissue biomechanical property.
McAleavey, Stephen A
2014-05-01
Shear wave induced phase encoding (SWIPE) imaging generates ultrasound backscatter images of tissue-like elastic materials by using traveling shear waves to encode the lateral position of the scatters in the phase of the received echo. In contrast to conventional ultrasound B-scan imaging, SWIPE offers the potential advantages of image formation without beam focusing or steering from a single transducer element, lateral resolution independent of aperture size, and the potential to achieve relatively high lateral resolution with low frequency ultrasound. Here a Fourier series description of the phase modulated echo signal is developed, demonstrating that echo harmonics at multiples of the shear wave frequency reveal target k-space data at identical multiples of the shear wavenumber. Modulation transfer functions of SWIPE imaging systems are calculated for maximum shear wave acceleration and maximum shear constraints, and compared with a conventionally focused aperture. The relative signal-to-noise ratio of the SWIPE method versus a conventionally focused aperture is found through these calculations. Reconstructions of wire targets in a gelatin phantom using 1 and 3.5 MHz ultrasound and a cylindrical shear wave source are presented, generated from the fundamental and second harmonic of the shear wave modulation frequency, demonstrating weak dependence of lateral resolution with ultrasound frequency. PMID:24815265
McAleavey, Stephen A.
2014-01-01
Shear wave induced phase encoding (SWIPE) imaging generates ultrasound backscatter images of tissue-like elastic materials by using traveling shear waves to encode the lateral position of the scatters in the phase of the received echo. In contrast to conventional ultrasound B-scan imaging, SWIPE offers the potential advantages of image formation without beam focusing or steering from a single transducer element, lateral resolution independent of aperture size, and the potential to achieve relatively high lateral resolution with low frequency ultrasound. Here a Fourier series description of the phase modulated echo signal is developed, demonstrating that echo harmonics at multiples of the shear wave frequency reveal target k-space data at identical multiples of the shear wavenumber. Modulation transfer functions of SWIPE imaging systems are calculated for maximum shear wave acceleration and maximum shear constraints, and compared with a conventionally focused aperture. The relative signal-to-noise ratio of the SWIPE method versus a conventionally focused aperture is found through these calculations. Reconstructions of wire targets in a gelatin phantom using 1 and 3.5 MHz ultrasound and a cylindrical shear wave source are presented, generated from the fundamental and second harmonic of the shear wave modulation frequency, demonstrating weak dependence of lateral resolution with ultrasound frequency. PMID:24815265
[Effective control of excitable waves in 2D cardiac excitable media].
Li, Li; Liu, Li; Zhang, Guangcai; Wang, Guangrui; Qu, Zhi
2005-12-01
We propose a method for effective control of patter on dissipative system by use of little perturbation analysis, and apply this nonuniform feed back method to control the polarization wave of heart represented by FHN equation. In the numerical experiment, we successfully alter the propagating direction of planar waves and drift the spiral waves to boundary without resonant repulsion. The effective control of the excited system will be used to study the mechanism of defibrillation and that is our interesting work. PMID:16422076
Acoustic wave flow sensor using quartz thickness shear mode resonator.
Qin, Lifeng; Zeng, Zijing; Cheng, Hongbin; Wang, Qing-Ming
2009-09-01
A quartz thickness shear mode (TSM) bulk acoustic wave resonator was used for in situ and real-time detection of liquid flow rate in this study. A special flow chamber made of 2 parallel acrylic plates was designed for flow measurement. The flow chamber has a rectangular flow channel, 2 flow reservoirs for stabilizing the fluid flow, a sensor mounting port for resonator holding, one inlet port, and one outlet port for pipe connection. A 5-MHz TSM quartz resonator was edge-bonded to the sensor mounting port with one side exposed to the flowing liquid and other side exposed to air. The electrical impedance spectra of the quartz resonator at different volumetric flow rate conditions were measured by an impedance analyzer for the extraction of the resonant frequency through a data-fitting method. The fundamental, 3rd, 5th, 7th, and 9th resonant frequency shifts were found to be around 920, 3572, 5947, 8228, and 10,300 Hz for flow rate variation from 0 to 3000 mL/min, which had a corresponding Reynolds number change from 0 to 822. The resonant frequency shifts of different modes are found to be quadratic with flow rate, which is attributed to the nonlinear effect of quartz resonator due to the effective normal pressure imposing on the resonator sensor by the flowing fluid. The results indicate that quartz TSM resonators can be used for flow sensors with characteristics of simplicity, fast response, and good repeatability. PMID:19811997
Mapping an aquitard breach using shear-wave seismic reflection
NASA Astrophysics Data System (ADS)
Waldron, B. A.; Harris, J. B.; Larsen, D.; Pell, A.
2009-05-01
In multi-layered hydrostratigraphic systems, aquitard breaches caused by faulting or paleo-erosion can allow substantial quantities of water of differing quality to be exchanged between aquifers. Seismic reflection technology was used to map the extent and orientation of an aquitard breach connecting a shallow alluvial aquifer to the deeper semi-confined Memphis aquifer in southwestern Tennessee, USA. Geophysical well logs indicate the presence of the aquitard at borehole locations that define the beginning and end points on two seismic survey lines, which intersect at a borehole where the aquitard is absent. A SE-NW-oriented paleochannel, 350 m wide and approximately 35-40 m deep, is interpreted from the seismic reflection surveys. The paleochannel cuts through the aquitard and into the upper part of the Memphis aquifer, thus creating a hydraulic connection between the shallow unconfined and deeper, semi-confined aquifers. The results indicate the potential of the shear-wave seismic reflection methods to resolve shallow breaches through fine-grained aquitards given availability of sufficient well control.
Statistical Analysis of Shear Wave Speed in the Uterine Cervix
Carlson, Lindsey C.; Feltovich, Helen; Palmeri, Mark L.; del Rio, Alejandro Muñoz; Hall, Timothy J.
2014-01-01
Although cervical softening is critical in pregnancy, there currently is no objective method for assessing the softness of the cervix. Shear wave speed (SWS) estimation is a noninvasive tool used to measure tissue mechanical properties such as stiffness. The goal of this study was to determine the spatial variability and assess the ability of SWS to classify ripened vs. unripened tissue samples. Ex vivo human hysterectomy samples (n = 22) were collected, a subset (n = 13) were ripened. SWS estimates were made at 4–5 locations along the length of the canal on both anterior and posterior halves. A linear mixed model was used for a robust multivariate analysis. Receiver operating characteristic (ROC) analysis and the area under the ROC curve (AUC) were calculated to describe the utility of SWS to classify ripened vs. unripened tissue samples. Results showed that all variables used in the linear mixed model were significant (p<0.05). Estimates at the mid location for the unripened group were 3.45 ± 0.95 m/s (anterior) and 3.56 ± 0.92 m/s (posterior), and 2.11 ± 0.45 m/s (anterior) and 2.68 ± 0.57 m/s (posterior) for the ripened (p < 0.001). The AUC’s were 0.91 and 0.84 for anterior and posterior respectively suggesting SWS estimates may be useful for quantifying cervical softening. PMID:25392863
Statistical analysis of shear wave speed in the uterine cervix.
Carlson, Lindsey C; Feltovich, Helen; Palmeri, Mark L; del Rio, Alejandro Muñoz; Hall, Timothy J
2014-10-01
Although cervical softening is critical in pregnancy, there currently is no objective method for assessing the softness of the cervix. Shear wave speed (SWS) estimation is a noninvasive tool used to measure tissue mechanical properties such as stiffness. The goal of this study was to determine the spatial variability and assess the ability of SWS to classify ripened versus unripened tissue samples. Ex vivo human hysterectomy samples (n = 22) were collected; a subset (n = 13) were ripened. SWS estimates were made at 4 to 5 locations along the length of the canal on both anterior and posterior halves. A linear mixed model was used for a robust multivariate analysis. Receiver operating characteristic (ROC) analysis and the area under the ROC curve (AUC) were calculated to describe the utility of SWS to classify ripened versus unripened tissue samples. Results showed that all variables used in the linear mixed model were significant ( p < 0.05). Estimates at the mid location for the unripened group were 3.45 ± 0.95 m/s (anterior) and 3.56 ± 0.92 m/s (posterior), and 2.11 ± 0.45 m/s (anterior) and 2.68 ± 0.57 m/s (posterior) for the ripened ( p < 0.001). The AUCs were 0.91 and 0.84 for anterior and posterior, respectively, suggesting that SWS estimates may be useful for quantifying cervical softening. PMID:25392863
Shear Wave Splitting Across Eastern, Western and Southern Africa
NASA Astrophysics Data System (ADS)
Nyblade, A.; Ramirez, C.; Bagley, B. C.; Mulibo, G. D.; Tugume, F.; Wysession, M. E.; Wiens, D. A.
2014-12-01
The expansion of the AfricaArray network across eastern, western and southern Africa, in conjunction with seismic data from many PASSCAL deployments over the past 20 years, is helping to fill in major gaps in the global coverage of shear wave splitting measurements. New results from stations in Ghana, Nigeria, Mozambique, Botswana, Angola, Namibia and South Africa are presented in this study that when combined with previously published measurements help to map the pattern of seismic anisotropy over much of the African continent. A general pattern of fast polarization directions, characterized by NE orientations, is found, and superimposed on this subcontinental-scale pattern is local and regional variability, most notably around the Archean Tanzania craton in eastern Africa. The subcontinental-scale pattern, as well as local and regional variations in this pattern, are interpreted in terms of large-scale mantle flow from the African Superplume, fossil anisotropy in the lithosphere, and shape anisotropy in magmatic regions of the East African rift system.
NASA Astrophysics Data System (ADS)
Weese, Juergen; Goecke, Roland; Penney, Graeme P.; Desmedt, Paul; Buzug, Thorsten M.; Schumann, Heidrun
1999-05-01
2D/3D registration makes it possible to use pre-operative CT scans for navigation purposes during X-ray fluoroscopy guided interventions. We present a fast voxel-based method for this registration task, which uses a recently introduced similarity measure (pattern intensity). This measure is especially suitable for 2D/3D registration, because it is robust with respect to structures such as a stent visible in the X-ray fluoroscopy image but not in the CT scan. The method uses only a part of the CT scan for the generation of digitally reconstructed radiographs (DRRs) to accelerate their computation. Nevertheless, computation time is crucial for intra-operative application and a further speed-up is required, because numerous DRRs must be computed. For that reason, the suitability of different volume rendering methods for 2D/3D registration has been investigated. A method based on the shear-warp factorization of the viewing transformation turned out to be especially suitable and builds the basis of the registration algorithm. The algorithm has been applied to images of a spine phantom and to clinical images. For comparison, registration results have been calculated using ray-casting. The shear-warp factorization based rendering method accelerates registration by a factor of up to seven compared to ray-casting without degrading registration accuracy. Using a vertebra as feature for registration, computation time is in the range of 3-4s (Sun UltraSparc, 300 MHz) which is acceptable for intra-operative application.
In Vivo Cardiac, Acoustic-Radiation-Force-Driven, Shear Wave Velocimetry
Hsu, Stephen J.; Wolf, Patrick D.; Trahey, Gregg E.
2009-01-01
Shear wave elasticity imaging (SWEI) was employed to track acoustic radiation force impulse (ARFI) -induced shear waves in the mid-myocardium of the left ventricular free wall (LVFW) of a beating canine heart. Shear waves were generated and tracked with a linear ultrasound transducer that was placed directly on the exposed epicardium. Acquinsition was ECG-gated arid coincided with the mid-diastolic portion of the cardiac cycle. Axial displacement profiles consistent with shear wave propagation were clearly evident in all SWEI acquisitions (i.e., those including an ARFI excitation); displacement data from control cases (i.e., sequences lacking an ARFI excitation) offered no evidence of shear wave propagation and yielded a peak absolute mean displacement below 0.31 μm after motion filtering. Shear wave velocity estimates ranged from 0.82 to 2.65 m/s and were stable across multiple heartbeats for the same interrogation region, with coefficients of variation less than 19% for all matched acquisitions. Variations in velocity estimates suggest a spatial dependence of shear wave velocity through the mid-myocardium of the LVFW, with velocity estimates changing, in limited cases, through depth and lateral position. PMID:19771962
Surface and downhole shear wave seismic methods for thick soil site investigations
Hunter, J.A.; Benjumea, B.; Harris, J.B.; Miller, R.D.; Pullan, S.E.; Burns, R.A.; Good, R.L.
2002-01-01
Shear wave velocity-depth information is required for predicting the ground motion response to earthquakes in areas where significant soil cover exists over firm bedrock. Rather than estimating this critical parameter, it can be reliably measured using a suite of surface (non-invasive) and downhole (invasive) seismic methods. Shear wave velocities from surface measurements can be obtained using SH refraction techniques. Array lengths as large as 1000 m and depth of penetration to 250 m have been achieved in some areas. High resolution shear wave reflection techniques utilizing the common midpoint method can delineate the overburden-bedrock surface as well as reflecting boundaries within the overburden. Reflection data can also be used to obtain direct estimates of fundamental site periods from shear wave reflections without the requirement of measuring average shear wave velocity and total thickness of unconsolidated overburden above the bedrock surface. Accurate measurements of vertical shear wave velocities can be obtained using a seismic cone penetrometer in soft sediments, or with a well-locked geophone array in a borehole. Examples from thick soil sites in Canada demonstrate the type of shear wave velocity information that can be obtained with these geophysical techniques, and show how these data can be used to provide a first look at predicted ground motion response for thick soil sites. ?? 2002 Published by Elsevier Science Ltd.
Johnston, J.E.; Christensen, N.I. . Dept. of Earth and Atmospheric Sciences)
1992-01-01
The physical properties of a sequence of Paleozoic sedimentary rocks have been examined in detail, with an emphasis on laboratory measurements of density, shear wave velocity, shear wave splitting, and Vp/Vs ratios. Seismic properties of 147 cores from 49 rock samples collected from the thorn hill sedimentary sequence of eastern Tennessee are examined in terms of implications for future seismic studies in the southern Appalachians. The shear wave velocities of these rocks are strongly influenced by the relatively high shear wave velocity of quartz. Shear wave velocity anisotropy is present in most of the lithologic groups: it is highest in the shales while being almost insignificant in the dolostones. The related phenomenon of shear wave splitting occurs to some degree in all of the lithologies studied and at high pressures originates from mineral orientation. Compressional to shear velocity (Vp/Vs) ratios of approximately 1.82 (dolostones) and 1.95 (limestones) effectively characterize the carbonates while other lithologies display wider ranges of Vp/Vs, primarily due to the influence of accessory minerals such as quartz. Densities of the sample suite range from 2.34 g/cm[sup 3] (shale) to 2.86 g/cm[sup 3] (dolostone). Normal incidence shear and compressional wave synthetic seismograms of the entire Thorn Hill section indicate that three zones of high amplitude reflections would be seen on reflection records obtained over this 3,327 meter thick sequence. differences are seen at some interfaces in the Mississippian-Devonian interval, which are more reflective to shear waves, and in the Ordovician Martinsburg Formation, which appears more reflective to compressional waves.
Coupling of an acoustic wave to shear motion due to viscous heating
NASA Astrophysics Data System (ADS)
Liu, Bin; Goree, J.
2016-07-01
Viscous heating due to shear motion in a plasma can result in the excitation of a longitudinal acoustic wave, if the shear motion is modulated in time. The coupling mechanism is a thermal effect: time-dependent shear motion causes viscous heating, which leads to a rarefaction that can couple into a longitudinal wave, such as an acoustic wave. This coupling mechanism is demonstrated in an electrostatic three-dimensional (3D) simulation of a dusty plasma, in which a localized shear flow is initiated as a pulse, resulting in a delayed outward propagation of a longitudinal acoustic wave. This coupling effect can be profound in plasmas that exhibit localized viscous heating, such as the dusty plasma we simulated using parameters typical of the PK-4 experiment. We expect that a similar phenomenon can occur with other kinds of plasma waves.
The propagation of horizontally polarized shear waves in plates bordered with viscous liquid.
Gitis, Alexander; Sauer, Dirk Uwe
2016-09-01
Requirements for ultrasonic horizontally polarized shear waves based viscosity sensors and their applicability for continuous in-line measurement are presented and discussed. The results reveal, that sensors using non-piezoelectric plates as well as wave guides and sensing surface have application-oriented advantages in corrosive and hot liquids. For such non-piezoelectric plate sensors, the dispersion relations are found and the linking equation among propagation velocity as well as attenuation coefficient and Newtonian liquid parameters are obtained. The findings show that in presence of viscous liquids the propagation parameters of horizontally polarized shear waves (HPSW) in non-piezoelectric plate change and a viscosity depending attenuation occurs. It is shown that the measurement sensitivity, in physical terms, of the investigated device highly depends on plate thickness, shear wave impedance of the plate material, and the shear wave impedance of the ambient liquid. Further, reasonable geometrical optimizations and suited plate materials are discussed. PMID:27423968
A hybrid wave-mode formulation for the vibro-acoustic analysis of 2D periodic structures
NASA Astrophysics Data System (ADS)
Droz, C.; Zhou, C.; Ichchou, M. N.; Lainé, J.-P.
2016-02-01
In the framework of vibrational analysis of 2D periodic waveguides, Floquet-Bloch theorem is widely applied for the determination of wave dispersion characteristics. In this context, the Wave Finite Element Method (WFEM) combines Periodic Structure Theory (PST) with standard FE packages, enabling wave dispersion analysis of waveguides involving structurally realistic unit-cells. For such applications, the computational efficiency of the WFEM depends on the choice of the formulation and can lead to numerical issues, worsen by extensive computational cost. This paper presents a coupled wave-mode approach for the determination of wave dispersion characteristics in structurally advanced periodic structures. It combines two scales of model order reduction. At the unit-cell's scale, Component Mode Synthesis (CMS) provides the displacement field associated with local resonances of the periodic structure, while the free wave propagation is considered using a spectral problem projection on a reduced set of shape functions associated with propagating waves, thus providing considerable reduction of the computational cost. An application is provided for a bi-directionally stiffened panel and the influence of reduction parameters is discussed, as well as the robustness of the numerical results.
Nonlinear shear wave in a non Newtonian visco-elastic medium
Banerjee, D.; Janaki, M. S.; Chakrabarti, N.
2012-06-15
An analysis of nonlinear transverse shear wave has been carried out on non-Newtonian viscoelastic liquid using generalized hydrodynamic model. The nonlinear viscoelastic behavior is introduced through velocity shear dependence of viscosity coefficient by well known Carreau-Bird model. The dynamical feature of this shear wave leads to the celebrated Fermi-Pasta-Ulam problem. Numerical solution has been obtained which shows that initial periodic solutions reoccur after passing through several patterns of periodic waves. A possible explanation for this periodic solution is given by constructing modified Korteweg de Vries equation. This model has application from laboratory to astrophysical plasmas as well as in biological systems.