skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Von Kármán spatial correlation function to describe wave propagation in polycrystalline media
Analytical functions that describe the spatial heterogeneity in polycrystalline media are highly desirable. These mathematically tractable descriptors can be readily implemented in physical models of static and dynamic material behavior, including wave propagation. This paper explores the suitability of von Kármán spatial correlation functions (SCFs) to describe polycrystalline media with a distribution of grain sizes. The empirical two-point statistics are compared to the von Kármán and other commonly reported SCFs. The von Kármán function is shown to be more accurate than the exponential function and more tractable than the sum of exponentials form. The impact of the SCF on wave propagation and scattering is studied by employing a well-defined analytical model for attenuation. The attenuation varies by over a factor of two for the aluminum case considered. These results provide preliminary insights into the suitability of a closed-form von Kármán SCF to describe polycrystalline media with increasingly complex microstructures.  more » « less
Award ID(s):
2029142
PAR ID:
10334466
Author(s) / Creator(s):
Date Published:
Journal Name:
Journal of Applied Physics
Volume:
131
Issue:
22
ISSN:
0021-8979
Page Range / eLocation ID:
225109
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. The impact of vertical wind shear on the land–sea-breeze circulation at the equator is explored using idealized 2D numerical simulations and a simple 2D linear analytical model. Both the idealized and linear analytical models indicate Doppler shifting and attenuation effects coexist under the effect of vertical wind shear for the propagation of gravity waves that characterize the land–sea-breeze circulation. Without a background wind, the idealized sea breeze has two ray paths of gravity waves that extend outward and upward from the coast. A uniform background wind causes a tilting of the two ray paths due to Doppler shifting. With vertical shear in the background wind, the downstream ray path of wave propagation can be rapidly attenuated near a certain level, whereas the upstream ray path is not attenuated and the amplitudes even increase with height. The downstream attenuation level is found to descend with increasing linear wind shear. The present analytical model establishes that the attenuation level corresponds to the critical level where the background wind is equal to the horizontal gravity wave phase speed. The upstream gravity wave ray path can propagate upward without attenuation as there is no critical level there. 
    more » « less
  2. Abstract Dislocation‐based dissipation mechanisms potentially control the viscoelastic response of Earth's upper mantle across a variety of geodynamic contexts, including glacial isostatic adjustment, postseismic creep, and seismic‐wave attenuation. However, there is no consensus on which dislocation‐based, microphysical process controls the viscoelastic behavior of the upper mantle. Although both intergranular (plastic anisotropy) and intragranular (backstress) mechanisms have been proposed, there is currently insufficient laboratory data to discriminate between those mechanisms. Here, we present the results of forced‐oscillation experiments in a deformation‐DIA apparatus at confining pressures of 3–7 GPa and temperatures of 298–1370 K. Our experiments tested the viscoelastic response of polycrystalline olivine—the main constituent of the upper mantle—at stress amplitudes from 70 to 2,800 MPa. Mechanical data are complemented by microstructural analyses of grain size, crystallographic preferred orientation, and dislocation density. We observe amplitude‐ and frequency‐dependent attenuation and modulus relaxation and find that numerical solutions of the backstress model match our results well. Therefore, we argue that interactions among dislocations, rather than intergranular processes (e.g., plastic anisotropy or grain boundary sliding), control the viscoelastic behavior of polycrystalline olivine in our experiments. In addition, we present a linearized version of the constitutive equations of the backstress model and extrapolate it to conditions typical of seismic‐wave propagation in the upper mantle. Our extrapolation demonstrates that the backstress model can explain the magnitude of seismic‐wave attenuation in the upper mantle, although some modification is required to explain the weak frequency dependence of attenuation observed in nature and in previous experimental work. 
    more » « less
  3. Abstract Early research in aerodynamics and biological propulsion was dramatically advanced by the analytical solutions of Theodorsen, von Kármán, Wu and others. While these classical solutions apply only to isolated swimmers, the flow interactions between multiple swimmers are relevant to many practical applications, including the schooling and flocking of animal collectives. In this work, we derive a class of solutions that describe the hydrodynamic interactions between an arbitrary number of swimmers in a two‐dimensional inviscid fluid. Our approach is rooted in multiply‐connected complex analysis and exploits several recent results. Specifically, the transcendental (Schottky–Klein) prime function serves as the basic building block to construct the appropriate conformal maps and leading‐edge‐suction functions, which allows us to solve the modified Schwarz problem that arises. As such, our solutions generalize classical thin aerofoil theory, specifically Wu's waving‐plate analysis, to the case of multiple swimmers. For the case of a pair of interacting swimmers, we develop an efficient numerical implementation that allows rapid computations of the forces on each swimmer. We investigate flow‐mediated equilibria and find excellent agreement between our new solutions and previously reported experimental results. Our solutions recover and unify disparate results in the literature, thereby opening the door for future studies into the interactions between multiple swimmers. 
    more » « less
  4. null (Ed.)
    Seismic waves in earth media usually undergo attenuation, causing energy losses and phase distortions. In the regime of high-frequency asymptotics, a complex-valued eikonal is an essential ingredient for describing wave propagation in attenuating media, where the real and imaginary parts of the eikonal function capture dispersion effects and amplitude attenuation of seismic waves, respectively. Conventionally, such a complex-valued eikonal is mainly computed either by tracing rays exactly in complex space or by tracing rays approximately in real space so that the resulting eikonal is distributed irregularly in real space. However, seismic data processing methods, such as prestack depth migration and tomography, usually require uniformly distributed complex-valued eikonals. Therefore, we have developed a unified framework to Eulerianize several popular approximate real-space ray-tracing methods for complex-valued eikonals so that the real and imaginary parts of the eikonal function satisfy the classic real-space eikonal equation and a novel real-space advection equation, respectively, and we dub the resulting method the Eulerian partial-differential-equation method. We further develop highly efficient high-order methods to solve these two equations by using the factorization idea and the Lax-Friedrichs weighted essentially nonoscillatory schemes. Numerical examples demonstrate that our method yields highly accurate complex-valued eikonals, analogous to those from ray-tracing methods. Our methods can be useful for migration and tomography in attenuating media. 
    more » « less
  5. This paper investigates the validity of two different analytical homogenization methods: the Mori–Tanaka mean-field theory and Milton’s correlation function-dependent bounds. We focus on biphase linearly elastic transversely isotropic composites. The composites consist of a matrix reinforced with long fibers of either circular or irregular cross section shapes formed by overlapping circles, with different degrees of radius polydispersity. The Mori–Tanaka effective stiffness depends on the phase moduli, volume fractions, and on a few geometric descriptors of the fibers that can be readily evaluated. In contrast, the computation of Milton’s bounds requires finer knowledge of the microstructure, in terms of two and three-point spatial correlation functions, which are not always analytically tractable. We thus consider very specific random microstructure geometries with known correlation functions. The effective moduli estimates of the two methods are validated against the results of numerical homogenization using the finite element method. It is shown that the Mori– Tanaka predictions of the effective transverse bulk modulus are significantly more accurate than those of the transverse and axial shear moduli. In addition, the predictions of the scheme generally deteriorate with an increasing fiber volume fraction. By contrast, the average of Milton’s upper and lower bounds provides a highly accurate estimate for all three independent effective moduli, without any limitation on the fiber concentration. This study highlights the indisputable effect of the spatial correlation functions on the effective properties of composites, and aspires to pave the way towards the development of more predictive, correlation function-dependent homogenization methods. 
    more » « less