Complex roots of the Stoneley-wave equation

1972 ◽  
Vol 62 (1) ◽  
pp. 285-299
Author(s):  
Walter L. Pilant

Abstract The equation governing elastic waves propagating along a solid-solid interface is found to have sixteen (16) independent roots on its eight (8) associated Riemann sheets. The range of existence (in terms of material parameters) for the real root corresponding to the propagation of Stoneley waves has long been known. It is found that outside this range there are two types of behavior. If the material of greater density has a velocity slightly greater than that of the material of lesser density, the unattenuated Stoneley waves make a transition to attenuated Interface waves, i.e., they leak energy away from the interface as they propagate along it. If the more dense material has a velocity more than about three times that of the less dense, then the Interface-wave root disappears and energy is propagated along the interface as Rayleigh waves. This Rayleigh-wave propagation is associated with a different root of the fundamental equation. On the other hand, if the material of greater density has a velocity much lower than that of the material of lower density (a case that is difficult to find physically), then no energy will be propagated along the interface at all. This result was unexpected. Some rather interesting behavior of the 16 roots was noted as the physical parameters were varied over a wide range. In addition to the normal collisions between pairs of roots, and between individual roots and branch points (with attendant Riemann sheet jumping), it was found that some roots go through the point at infinity and return with a change in sign. At least one unexpected case of a multiple root was found. Another case was noted in which a pair of complex roots change quadrants in the complex phase-velocity plane, leading to a discontinuity in root type. Finally, it was noted that, in a cyclic variation of the material parameters, it is possible to choose a path such that the roots, when followed individually, will not return to their original values. In fact, as many as five cycles in parameter space can be accomplished before the roots return. All this strange mathematical behavior seems to have no physical significance, but has been presented to increase understanding of the general behavior of the dispersion relations associated with elastic-wave propagation.

Geophysics ◽  
2011 ◽  
Vol 76 (4) ◽  
pp. E127-E139 ◽  
Author(s):  
Robert K. Mallan ◽  
Carlos Torres-Verdín ◽  
Jun Ma

A numerical simulation study has been made of borehole sonic measurements that examined shoulder-bed, anisotropy, and mud-filtrate invasion effects on frequency-dispersion curves of flexural and Stoneley waves. Numerical simulations were considered for a range of models for fast and slow formations. Computations are performed with a Cartesian 3D finite-difference time-domain code. Simulations show that presence of transverse isotropy (TI) alters the dispersion of flexural and Stoneley waves. In slow formations, the flexural wave becomes less dispersive when the shear modulus (c44) governing wave propagation parallel to the TI symmetry axis is lower than the shear modulus (c66) governing wave propagation normal to the TI symmetry axis; conversely, the flexural wave becomes more dispersive when c44 > c66. Dispersion decreases by as much as 30% at higher frequencies for the considered case where c44 < c66. Dispersion of Stoneley waves, on the other hand, increases in TI formations when c44 > c66 and decreases when c44 < c66. Dispersion increases by more than a factor of 2.5 at higher frequencies for the considered case where c44 < c66. Simulations also indicate that the impact of invasion on flexural and Stoneley dispersions can be altered by the presence of TI. For the case of a slow formation and TI, where c44 decreases from the isotropic value, separation between dispersion curves for cases with and without the presence of a fast invasion zone increases by as much as 33% for the flexural wave and by as much as a factor of 1.4 for the Stoneley wave. Lastly, presence of a shoulder bed intersecting the sonic tool at high dip angles can alter flexural dispersion significantly at low frequencies. For the considered case of a shoulder bed dipping at 80°, ambiguity in the flexural cutoff frequency might lead to shear-wave velocity errors of 8%–10%.


Geophysics ◽  
1991 ◽  
Vol 56 (4) ◽  
pp. 447-460 ◽  
Author(s):  
X. M. Tang ◽  
C. H. Cheng ◽  
M. N. Toksöz

The propagation of Stoneley waves in a fluid‐filled borehole with a vertical fracture is investigated both theoretically and experimentally. The borehole propagation excites fluid motion in the fracture and the resulting fluid flow at the fracture opening perturbs the fluid‐solid interface boundary condition at the borehole wall. By developing a boundary condition perturbation technique for the borehole situation, we studied the effect of this change in the boundary condition on the Stoneley propagation. Cases of both hard and soft formations have been investigated. The fracture has minimal effects on the Stoneley velocity, except in the very low frequency range in which the Stoneley velocity drastically decreases with decreasing frequency. Significant Stoneley‐wave attenuation is produced because of the energy dissipation into the fracture. The quantitative behavior of these effects depends not only on fracture aperture and borehole radius, but also on the acoustic properties of the formation and fluid. Ultrasonic experiments were performed to measure Stoneley propagation in laboratory fracture borehole models. Aluminum and lucite were used to simulate a hard and a soft formation, respectively. Array data for wave propagation were obtained and were processed using Prony’s method to give velocity and attenuation of Stoneley waves as a function of frequency. In both hard and soft formation cases, the experimental results agreed with the theoretical predictions. The important result of this study is that it presents a quantitative relationship between the Stoneley propagation and the fracture character in conjunction with formation and fluid properties. This relationship provides a method for estimating the characteristics of a vertical fracture by means of Stoneley wave measurements.


1990 ◽  
Vol 43 (5S) ◽  
pp. S359-S365 ◽  
Author(s):  
Arthur M. B. Braga ◽  
Paul E. Barbone ◽  
George Herrmann

Harmonic waves propagating in a fluid-loaded, composite circular cylinder are investigated. The solid in this case is composed of perfectly bonded, isotropic, linearly elastic, cylindrical layers stacked radially to form a laminated cylinder. The surface impedance of the laminated shell is calculated using a recursive invariant-imbedding algorithm. This algorithm is numerically stable for a wide range of frequencies and circumferential mode numbers. Particular attention is given to the study of subsonic interface waves which propagate unattenuated along a fluid-solid interface. Results are presented which show the effects of layering and curvature on the dispersion spectrum of the subsonic interface waves.


Sensors ◽  
2021 ◽  
Vol 21 (13) ◽  
pp. 4459
Author(s):  
José R. González ◽  
Charbel Damião ◽  
Maira Moran ◽  
Cristina A. Pantaleão ◽  
Rubens A. Cruz ◽  
...  

According to experts and medical literature, healthy thyroids and thyroids containing benign nodules tend to be less inflamed and less active than those with malignant nodules. It seems to be a consensus that malignant nodules have more blood veins and more blood circulation. This may be related to the maintenance of the nodule’s heat at a higher level compared with neighboring tissues. If the internal heat modifies the skin radiation, then it could be detected by infrared sensors. The goal of this work is the investigation of the factors that allow this detection, and the possible relation with any pattern referent to nodule malignancy. We aim to consider a wide range of factors, so a great number of numerical simulations of the heat transfer in the region under analysis, based on the Finite Element method, are performed to study the influence of each nodule and patient characteristics on the infrared sensor acquisition. To do so, the protocol for infrared thyroid examination used in our university’s hospital is simulated in the numerical study. This protocol presents two phases. In the first one, the body under observation is in steady state. In the second one, it is submitted to thermal stress (transient state). Both are simulated in order to verify if it is possible (by infrared sensors) to identify different behavior referent to malignant nodules. Moreover, when the simulation indicates possible important aspects, patients with and without similar characteristics are examined to confirm such influences. The results show that the tissues between skin and thyroid, as well as the nodule size, have an influence on superficial temperatures. Other thermal parameters of thyroid nodules show little influence on surface infrared emissions, for instance, those related to the vascularization of the nodule. All details of the physical parameters used in the simulations, characteristics of the real nodules and thermal examinations are publicly available, allowing these simulations to be compared with other types of heat transfer solutions and infrared examination protocols. Among the main contributions of this work, we highlight the simulation of the possible range of parameters, and definition of the simulation approach for mapping the used infrared protocol, promoting the investigation of a possible relation between the heat transfer process and the data obtained by infrared acquisitions.


2021 ◽  
Vol 7 (2) ◽  
pp. 48
Author(s):  
Vittorio Scardaci ◽  
Giuseppe Compagnini

Laser scribing has been proposed as a fast and easy tool to reduce graphene oxide (GO) for a wide range of applications. Here, we investigate laser reduction of GO under a range of processing and material parameters, such as laser scan speed, number of laser passes, and material coverage. We use Raman spectroscopy for the characterization of the obtained materials. We demonstrate that laser scan speed is the most influential parameter, as a slower scan speed yields poor GO reduction. The number of laser passes is influential where the material coverage is higher, producing a significant improvement of GO reduction on a second pass. Material coverage is the least influential parameter, as it affects GO reduction only under restricted conditions.


2021 ◽  
Vol 504 (2) ◽  
pp. 2325-2345
Author(s):  
Emanuel Sillero ◽  
Patricia B Tissera ◽  
Diego G Lambas ◽  
Stefano Bovino ◽  
Dominik R Schleicher ◽  
...  

ABSTRACT We present p-gadget3-k, an updated version of gadget-3, that incorporates the chemistry package krome. p-gadget3-k follows the hydrodynamical and chemical evolution of cosmic structures, incorporating the chemistry and cooling of H2 and metal cooling in non-equilibrium. We performed different runs of the same ICs to assess the impact of various physical parameters and prescriptions, namely gas metallicity, molecular hydrogen formation on dust, star formation recipes including or not H2 dependence, and the effects of numerical resolution. We find that the characteristics of the simulated systems, both globally and at kpc-scales, are in good agreement with several observable properties of molecular gas in star-forming galaxies. The surface density profiles of star formation rate (SFR) and H2 are found to vary with the clumping factor and resolution. In agreement with previous results, the chemical enrichment of the gas component is found to be a key ingredient to model the formation and distribution of H2 as a function of gas density and temperature. A star formation algorithm that takes into account the H2 fraction together with a treatment for the local stellar radiation field improves the agreement with observed H2 abundances over a wide range of gas densities and with the molecular Kennicutt–Schmidt law, implying a more realistic modelling of the star formation process.


2012 ◽  
Vol 717-720 ◽  
pp. 1101-1104 ◽  
Author(s):  
M.G. Jaikumar ◽  
Shreepad Karmalkar

4H-Silicon Carbide VDMOSFET is simulated using the Sentaurus TCAD package of Synopsys. The simulator is calibrated against measured data for a wide range of bias conditions and temperature. Material parameters of 4H-SiC are taken from literature and used in the available silicon models of the simulator. The empirical parameters are adjusted to get a good fit between the simulated curves and measured data. The simulation incorporates the bias and temperature dependence of important physical mechanisms like interface trap density, coulombic interface trap scattering, surface roughness scattering and velocity saturation.


2021 ◽  
Author(s):  
Tim Brandes ◽  
Stefano Scarso ◽  
Christian Koch ◽  
Stephan Staudacher

Abstract A numerical experiment of intentionally reduced complexity is used to demonstrate a method to classify flight missions in terms of the operational severity experienced by the engines. In this proof of concept, the general term of severity is limited to the erosion of the core flow compressor blade and vane leading edges. A Monte Carlo simulation of varying operational conditions generates a required database of 10000 flight missions. Each flight is sampled at a rate of 1 Hz. Eleven measurable or synthesizable physical parameters are deemed to be relevant for the problem. They are reduced to seven universal non-dimensional groups which are averaged for each flight. The application of principal component analysis allows a further reduction to three principal components. They are used to run a support-vector machine model in order to classify the flights. A linear kernel function is chosen for the support-vector machine due to its low computation time compared to other functions. The robustness of the classification approach against measurement precision error is evaluated. In addition, a minimum number of flights required for training and a sensible number of severity classes are documented. Furthermore, the importance to train the algorithms on a sufficiently wide range of operations is presented.


Ocean Science ◽  
2009 ◽  
Vol 5 (4) ◽  
pp. 661-684 ◽  
Author(s):  
C. Moore ◽  
A. Barnard ◽  
P. Fietzek ◽  
M. R. Lewis ◽  
H. M. Sosik ◽  
...  

Abstract. Requirements for understanding the relationships between ocean color and suspended and dissolved materials within the water column, and a rapidly emerging photonics and materials technology base for performing optical based analytical techniques have generated a diverse offering of commercial sensors and research prototypes that perform optical measurements in water. Through inversion, these tools are now being used to determine a diverse set of related biogeochemical and physical parameters. Techniques engaged include measurement of the solar radiance distribution, absorption, scattering, stimulated fluorescence, flow cytometry, and various spectroscopy methods. Selective membranes and other techniques for material isolation further enhance specificity, leading to sensors for measurement of dissolved oxygen, methane, carbon dioxide, common nutrients and a variety of other parameters. Scientists are using these measurements to infer information related to an increasing set of parameters and wide range of applications over relevant scales in space and time.


2020 ◽  
Vol 197 ◽  
pp. 01004
Author(s):  
Martina Capone ◽  
Elisa Guelpa ◽  
Vittorio Verda

As District Heating (DH) networks are experiencing an evolution towards the so-called 4th generation, there is a need to update the currently used models to take into account the ever-increasing complexity of this technology. Indeed, to further improve the reduction in energy consumption and carbon-dioxide emissions, a wide range of technologies and management strategies are being introduced within district heating, such as a large exploitation of Renewable Energy Sources (RES). As a consequence, thermal transients assume a major importance, posing the need to redefine the relevant physical parameters and to develop a model which accurately describes their behaviour. In this framework, this paper proposes a quantitative analysis of the influence of the pipe heat-capacity on the model. Moreover, an equivalent-model, which is able to take into account the two heat capacities of steel and water in just one equation, is proposed and compared with two commonly used approaches. One of the features of the proposed model is the suitability for application to large networks. To prove its capabilities, an application to the Turin district heating network, which is among the largest systems in Europe, is proposed. Results show significant improvements in terms of accuracy over computational time ratio.


Sign in / Sign up

Export Citation Format

Share Document