skip to main content

Title: Green's function for anisotropic dispersive poroelastic media based on the Radon transform and eigenvector diagonalization
A compact Green's function for general dispersive anisotropic poroelastic media in a full-frequency regime is presented for the first time. First, starting in a frequency domain, the anisotropic dispersion is exactly incorporated into the constitutive relationship, thus avoiding fractional derivatives in a time domain. Then, based on the Radon transform, the original three-dimensional differential equation is effectively reduced to a one-dimensional system in space. Furthermore, inspired by the strategy adopted in the characteristic analysis of hyperbolic equations, the eigenvector diagonalization method is applied to decouple the one-dimensional vector problem into several independent scalar equations. Consequently, the fundamental solutions are easily obtained. A further derivation shows that Green's function can be decomposed into circumferential and spherical integrals, corresponding to static and transient responses, respectively. The procedures shown in this study are also compatible with other pertinent multi-physics coupling problems, such as piezoelectric, magneto-electro-elastic and thermo-elastic materials. Finally, the verifications and validations with existing analytical solutions and numerical solvers corroborate the correctness of the proposed Green's function.
; ; ; ;
Award ID(s):
Publication Date:
Journal Name:
Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences
Page Range or eLocation-ID:
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract The paper extends the recent work (JAM, 88, 061002, 2021) of the Eshelby's tensors for polynomial eigenstrains from a two dimensional (2D) to three dimensional (3D) domain, which provides the solution to the elastic field with continuously distributed eigenstrain on a polyhedral inclusion approximated by the Taylor series of polynomials. Similarly, the polynomial eigenstrain is expanded at the centroid of the polyhedral inclusion with uniform, linear and quadratic order terms, which provides tailorable accuracy of the elastic solutions of polyhedral inhomogeneity by using Eshelby's equivalent inclusion method. However, for both 2D and 3D cases, the stress distribution in the inhomogeneity exhibits a certain discrepancy from the finite element results at the neighborhood of the vertices due to the singularity of Eshelby's tensors, which makes it inaccurate to use the Taylor series of polynomials at the centroid to catch the eigenstrain at the vertices. This paper formulates the domain discretization with tetrahedral elements to accurately solve for eigenstrain distribution and predict the stress field. With the eigenstrain determined at each node, the elastic field can be predicted with the closed-form domain integral of Green's function. The parametric analysis shows the performance difference between the polynomial eigenstrain by the Taylor expansionmore »at the centroid and the 𝐶0 continuous eigenstrain by particle discretization. Because the stress singularity is evaluated by the analytical form of the Eshelby's tensor, the elastic analysis is robust, stable and efficient.« less
  2. Abstract We outline and interpret a recently developed theory of impedance matching or reflectionless excitation of arbitrary finite photonic structures in any dimension. The theory includes both the case of guided wave and free-space excitation. It describes the necessary and sufficient conditions for perfectly reflectionless excitation to be possible and specifies how many physical parameters must be tuned to achieve this. In the absence of geometric symmetries, such as parity and time-reversal, the product of parity and time-reversal, or rotational symmetry, the tuning of at least one structural parameter will be necessary to achieve reflectionless excitation. The theory employs a recently identified set of complex frequency solutions of the Maxwell equations as a starting point, which are defined by having zero reflection into a chosen set of input channels, and which are referred to as R-zeros. Tuning is generically necessary in order to move an R-zero to the real frequency axis, where it becomes a physical steady-state impedance-matched solution, which we refer to as a reflectionless scattering mode (RSM). In addition, except in single-channel systems, the RSM corresponds to a particular input wavefront, and any other wavefront will generally not be reflectionless. It is useful to consider the theory asmore »representing a generalization of the concept of critical coupling of a resonator, but it holds in arbitrary dimension, for arbitrary number of channels, and even when resonances are not spectrally isolated. In a structure with parity and time-reversal symmetry (a real dielectric function) or with parity–time symmetry, generically a subset of the R-zeros has real frequencies, and reflectionless states exist at discrete frequencies without tuning. However, they do not exist within every spectral range, as they do in the special case of the Fabry–Pérot or two-mirror resonator, due to a spontaneous symmetry-breaking phenomenon when two RSMs meet. Such symmetry-breaking transitions correspond to a new kind of exceptional point, only recently identified, at which the shape of the reflection and transmission resonance lineshape is flattened. Numerical examples of RSMs are given for one-dimensional multimirror cavities, a two-dimensional multiwaveguide junction, and a multimode waveguide functioning as a perfect mode converter. Two solution methods to find R-zeros and RSMs are discussed. The first one is a straightforward generalization of the complex scaling or perfectly matched layer method and is applicable in a number of important cases; the second one involves a mode-specific boundary matching method that has only recently been demonstrated and can be applied to all geometries for which the theory is valid, including free space and multimode waveguide problems of the type solved here.« less
  3. Consider the elastic scattering of a time-harmonic wave by multiple well-separated rigid particles with smooth boundaries in two dimensions. Instead of using the complex Green's tensor of the elastic wave equation, we utilize the Helmholtz decomposition to convert the boundary value problem of the elastic wave equation into a coupled boundary value problem of the Helmholtz equation. Based on single, double, and combined layer potentials with the simpler Green's function of the Helmholtz equation, we present three different boundary integral equations for the coupled boundary value problem. The well-posedness of the new integral equations is established. Computationally, a scattering matrix based method is proposed to evaluate the elastic wave for arbitrarily shaped particles. The method uses the local expansion for the incident wave and the multipole expansion for the scattered wave. The linear system of algebraic equations is solved by GMRES with fast multipole method (FMM) acceleration. Numerical results show that the method is fast and highly accurate for solving elastic scattering problems with multiple particles.
  4. The chief objective of this paper is to explore energy transfer mechanism between the sub-systems that are coupled by a nonlinear elastic path. In the proposed model (via a minimal order, two degree of freedom system), both sub-systems are defined as damped harmonic oscillators with linear springs and dampers. The first sub-system is attached to the ground on one side but connected to the second sub-system on the other side. In addition, linear elastic and dissipative characteristics of both oscillators are assumed to be identical, and a harmonic force excitation is applied only on the mass element of second oscillator. The nonlinear spring (placed in between the two sub-systems) is assumed to exhibit cubic, hardening type nonlinearity. First, the governing equations of the two degree of freedom system with a nonlinear elastic path are obtained. Second, the nonlinear differential equations are solved with a semi-analytical (multi-term harmonic balance) method, and nonlinear frequency responses of the system are calculated for different path coupling cases. As such, the nonlinear path stiffness is gradually increased so that the stiffness ratio of nonlinear element to the linear element is 0.01, 0.05, 0.1, 0.5 and 1.0 while the absolute value of linear spring stiffness ismore »kept intact. In all solutions, it is observed that the frequency response curves at the vicinity of resonant frequencies bend towards higher frequencies as expected due to the hardening effect. However, at moderate or higher levels of path coupling (say 0.1, 0.5 and 1.0), additional branches emerge in the frequency response curves but only at the first resonant frequency. This is due to higher displacement amplitudes at the first resonant frequency as compared to the second one. Even though the oscillators move in-phase around the first natural frequency, high amplitudes increase the contribution of the stored potential energy in the nonlinear spring to the total mechanical energy. The out-of-phase motion around the second natural frequency cannot significantly contribute due to very low motion amplitudes. Finally, the governing equations are numerically solved for the same levels of nonlinearity, and the motion responses of both sub-systems are calculated. Both in-phase and out-of-phase motion responses are successfully shown in numerical solutions, and phase portraits of the system are generated in order to illustrate its nonlinear dynamics. In conclusion, a better understanding of the effect of nonlinear elastic path on two damped harmonic oscillators is gained.« less
  5. We investigate aspects of the spherical squirmer model employing both large-scale numerical simulations and asymptotic methods when the squirmer is placed in weakly elastic fluids. The fluids are modelled by differential equations, including the upper-convected Maxwell (UCM)/Oldroyd-B, finite-extensibility nonlinear elastic model with Peterlin approximation (FENE-P) and Giesekus models. The squirmer model we examine is characterized by two dimensionless parameters related to the fluid velocity at the surface of the micro-swimmer: the slip parameter $\xi $ and the swirl parameter $\zeta $ . We show that, for swimming in UCM/Oldroyd-B fluids, the elastic stress becomes singular at a critical Weissenberg number, Wi , that depends only on $\xi$ . This singularity for the UCM/Oldroyd-B models is independent of the domain exterior to the swimmer, or any other forces considered in the momentum balance for the fluid – we believe that this is the first time such a singularity has been explicitly demonstrated. Moreover, we show that the behaviour of the solution at the poles is purely extensional in character and is the primary reason for the singularity in the Oldroyd-B model. When the Giesekus or the FENE-P models are utilized, the singularity is removed. We also investigate the mechanism behind themore »speed and rotation rate enhancement associated with the addition of swirl in the swimmer's gait. We demonstrate that, for all models, the speed is enhanced by swirl, but the mechanism of enhancement depends intrinsically on the rheological model employed. Special attention is paid to the propulsive role of the pressure and elucidated upon. We also study the region of convergence of the perturbation solutions in terms of Wi . When techniques that accelerate the convergence of series are applied, transformed solutions are derived that are in very good agreement with the results obtained by simulations. Finally, both the analytical and numerical results clearly indicate that the low- Wi region is more important than one would expect and demonstrate all the major phenomena observed when swimming with swirl in a viscoelastic fluid.« less