A general transformation for compact waveguide coupler by using homogeneous media

Full-field numerical solutions for a crack which lies along the interface of an elastic-plastic medium and a rigid substrate are presented. The solutions are obtained using a small strain version of the J2-deformation theory with power-law strain hardening. In the present article, results for loading causing only small scale yielding at the crack tip are described; in subsequent articles the mathematical structure of the crack-tip fields under small scale yielding and results for contained yielding and fully plastic behavior will be presented. We find that although the near-tip fields do not appear to have a separable singular form, of the HRR-type fields as in homogeneous media, they do, however, bear interesting similarities to certain mixed-mode HRR fields. Under small scale yielding, where the remote elastic fields are specified by a complex stress-concentration vector Q = |Q|eiφ with φ being the phase angle between the two in-plane stress modes, we find that the plastic fields are members of a family parameterized by a new phase angle ξ, ≡ φ + εln(QQ/σ02L), and the fields nearly scale with the well-defined energy release rate as evaluated by the J-integral. Here σ0 is the reference yield stress and L is the total crack length (or a relevant length of the crack geometry). Numerical procedures appropriate for solving a general class of interface crack problems are also presented. A description of a numerical method for extracting the mixed mode stress intensities for cracks at interfaces and in homogeneous isotropic or anisotropic media, is included.

Deformable, rigid, and inviscid elliptical inclusions in a homogeneous incompressible anisotropic viscous fluid

This paper investigated numerically a natural convection in a porous cavity saturated by nanofluide. The left and right wall of the cavity are maintained at the hot-cold temperature respectively, the other walls are adiabatic. The two-phase Buongiorno model has been adopted to take account Brownian and thermophoretic diffusion in order to demonstrate the spatial distribution of the local nanoparticles concentration. After following the temporal evolution of the different structures (0<τ<5), Numerical simulations are performed to explore the effect of density buoyancy (104<Ra<106), the permeability of the homogeneous porous medium (10-5<Da<10-2), on the hydrodynamic, thermal and mass behavior. An original motivation was also introduced in our work to examine the effect of linearly variable permeability along the opposite direction of the cavity by varying the initial Darcy number from (10-5<Da<10-2) and fixing the final Darcy number Daf =10-5. The dimensionless partial differential equations are solved using the finite element method. The effects of the governing parameters on heat transfer are analyzed. Results indicate that the stationary regime is formed after the unsteady regime at dimensionless time τ = 0.3. The movement of nanofluide is strongly influenced by thermal buoyancy forces and depends on the Darcy number, heat transfer is accentuated for the homogenous medium compared to this one with variable permeability. It is found that the convective flow in a homogeneous porous medium is considerably affected by the variation of permeability and consequently the heat transfer is reduced.

Elastic study on singularities interacting with interfaces using alternating technique

A hyper-lens with planar object/image surfaces and uniform magnification factor is designed based on directional projecting property of optic-null medium (ONM). The proposed planar hyper-lens only requires one homogeneous anisotropic medium, i.e. ONM. We use metallic plates inside ϵ-near-zero medium to realize the proposed hyper-lens. 0.191λ 0 resolution is observed (0.238λ 0 can be clearly resolved) in experimental measurement, which consists very well with 3D numerical simulations and verifies the performance of the proposed planar hyper-lens. The planar hyper-lens with uniform pre-designed magnification factor may have some applications in super-resolution imaging technology.

The attenuation vector plays an important role in investigating waves propagating in dissipative, isotropic or anisotropic media. It is defined as the imaginary part of the complex-valued traveltime gradient. The real part of the complex-valued traveltime gradient is referred to as the propagation vector. In this paper, simple expressions for the attenuation vector are derived for seismic body waves propagating in heterogeneous, weakly dissipative, anisotropic media. The ray-theory perturbation method is used, in which the dissipative medium is considered to be a small perturbation of a perfectly elastic medium. The general perturbation procedure proposed by Klimeš (2002) is very suitable for this purpose, as it does not require ray tracing in dissipative media; the computation of rays in the reference perfectly elastic medium, followed by quadratures along them, is quite sufficient. It is shown that the waves propagating in a heterogeneous, weakly dissipative, anisotropic or isotropic medium are, in general, inhomogeneous. This means that their attenuation vector is not parallel to the propagation vector. This holds even for waves generated by a point-source in a homogeneous anisotropic weakly dissipative medium. An exception is a wave generated by a point-source in a homogeneous isotropic dissipative medium, where the generated wave is homogeneous. Thus, the commonly used concept of homogeneous waves can be applied neither to heterogeneous nor to anisotropic dissipative media. The situation is different for plane waves propagating in homogeneous dissipative isotropic or anisotropic media. In this case, the homogeneity or inhomogeneity of the plane wave may be chosen freely. Besides the attenuation vector, we also study the complex-valued ray-velocity vector. In heterogeneous media, the complex-valued ray-velocity vector is generally inhomogeneous, that is, its real and imaginary parts are not parallel. It is also shown that twice the scalar product of the attenuation vector with the ray-velocity vector in the reference medium yields 1/Q, where Q is the position- and direction-dependent quality factor. Quality factor Q does not depend on the inhomogeneity of the wave under consideration and offers a convenient measure of intrinsic material dissipation.

Nearly all dip‐moveout correction (DMO) implementations to date assume isotropic homogeneous media. Usually, this has been acceptable considering the tremendous cost savings of homogeneous isotropic DMO and considering the difficulty of obtaining the anisotropy parameters required for effective implementation. In the presence of typical anisotropy, however, ignoring the anisotropy can yield inadequate results. Since anisotropy may introduce large deviations from hyperbolic moveout, accurate transformation to zero‐offset in anisotropic media should address such nonhyperbolic moveout behavior of reflections. Artley and Hale’s v(z) ray‐tracing‐based DMO, developed for isotropic media, provides an attractive approach to treating such problems. By using a ray‐tracing procedure crafted for anisotropic media, I modify some aspects of their DMO so that it can work for v(z) anisotropic media. DMO impulse responses in typical transversely isotropic (TI) models (such as those associated with shales) deviate substantially from the familiar elliptical shape associated with responses in homogeneous isotropic media (to the extent that triplications arise even where the medium is homogeneous). Such deviations can exceed those caused by vertical inhomogeneity, thus emphasizing the importance of taking anisotropy into account in DMO processing. For isotropic or elliptically anisotropic media, the impulse response is an ellipse; but as the key anisotropy parameter η varies, the shape of the response differs substantially from elliptical. For typical η > 0, the impulse response in TI media tends to broaden compared to the response in an isotropic homogeneous medium, a behavior opposite to that encountered in typical v(z) isotropic media, where the response tends to be squeezed. Furthermore, the amplitude distribution along the DMO operator differs significantly from that for isotropic media. Application of this anisotropic DMO to data from offshore Africa resulted in a considerably better alignment of reflections from horizontal and dipping reflectors in common‐midpoint gather than that obtained using an isotropic DMO. Even the presence of vertical inhomogeneity in this medium could not eliminate the importance of considering the shale‐induced anisotropy.

Temporal changes in seismic anisotropy can be interpreted as variations in the orientation of cracks in seismogenic zones, and thus as variations in the stress field. Such temporal changes have been observed in seismogenic zones before and after earthquakes, although they are still not well understood. In this study, we investigate the azimuthal polarization of surface waves in anisotropic media with respect to the orientation of anisotropy, from a numerical point of view. This technique is based on the observation of the signature of anisotropy on the nine-component cross-correlation tensor (CCT) computed from seismic ambient noise recorded on pairs of three-component sensors. If noise sources are spatially distributed in a homogeneous medium, the CCT allows the reconstruction of the surface wave Green's tensor between the station pairs. In homogeneous, isotropic medium, four off-diagonal terms of the surface wave Green's tensor are null, but not in anisotropic medium. This technique is applied to three-component synthetic seismograms computed in a transversely isotropic medium with a horizontal symmetry axis, using a spectral element code. The CCT is computed between each pair of stations and then rotated, to approximate the surface wave Green's tensor by minimizing the off-diagonal components. This procedure allows the calculation of the azimuthal variation of quasi-Rayleigh and quasi-Love waves. In an anisotropic medium, in some cases, the azimuth of seismic anisotropy can induce a large variation in the horizontal polarization of surface waves. This variation depends on the relative angle between a pair of stations and the direction of anisotropy, the amplitude of the anisotropy, the frequency band of the signal and the depth of the anisotropic layer.

A 2D finite-difference, frequency-domain method was developed for modeling viscoacoustic seismic waves in transversely isotropic media with a tilted symmetry axis. The medium is parameterized by the P-wave velocity on the symmetry axis, the density, the attenuation factor, Thomsen's anisotropic parameters δ and ϵ, and the tilt angle. The finite-difference discretization relies on a parsimonious mixed-grid approach that designs accurate yet spatially compact stencils. The system of linear equations resulting from discretizing the time-harmonic wave equation is solved with a parallel direct solver that computes monochromatic wavefields efficiently for many sources. Dispersion analysis shows that four grid points per P-wavelength provide sufficiently accurate solutions in homogeneous media. The absorbing boundary conditions are perfectly matched layers (PMLs). The kinematic and dynamic accuracy of the method wasassessed with several synthetic examples which illustrate the propagation of S-waves excited at the source or at seismic discontinuities when ϵ<δ. In frequency-domain modeling with absorbing boundary conditions, the unstable S-wave mode is not excited when ϵ<δ, allowing stable simulations of the P-wave mode for such anisotropic media. Some S-wave instabilities are seen in the PMLs when the symmetry axis is tilted and ϵ>δ. These instabilities are consistent with previous theoretical analyses of PMLs in anisotropic media; they are removed if the grid interval is matched to the P-wavelength that leads to dispersive S-waves. Comparisons between seismograms computed with the frequency-domain acoustic TTI method and a finite-difference, time-domain method for the vertical transversely isotropic elastic equation show good agreement for weak to moderate anisotropy. This suggests the method can be used as a forward problem for viscoacoustic anisotropic full-waveform inversion.

Even though the propagation of elastic/acoustic waves in inhomogeneous and layered media has been an active research topic for many decades already, new problems and challenges continue to be posed even up to now. In fact, during the last few years, renewed interests have been witnessed by researchers in the various fields of acoustics, such as acoustic mirrors, filters, resonators, waveguides, and other kinds of acoustic devices, in relation to wave propagation in periodic elastic media. In acoustics and applied mechanics, these developments have been triggered by the need for new acoustic devices in order to obtain quality control of elastic/acoustic waves. What sort of material can allow us to have complete control over the elastic/acoustic wave’s propagation? We would like to discuss and answer this question in this chapter. It is well known that the successful applications of photonic band-gap materials have hastened the related researches on phononic band-gap materials. Analysis of Acoustic Wave in Homogeneous and Inhomogeneous Media Using Finite Element Method explores the theoretical road leading to the possible applications of phononic band gaps. It should quickly bring the elastic/acoustic professionals and engineers up to speed in this field of study where elastic/acoustic waves and solid-state physics meet. It will also provide an excellent overview to any course in elastic/acoustic media. Previous research on photonic crystals (Johnson & Joannopoulos, 2001, 2003; Joannopoulos et al., 1995; Leung & Liu, 1990) has sparked rapidly growing interest in the analogous acoustic effects of phononic crystals and periodic elastic structures. The various techniques for band structure calculations were introduced (Hussein, 2009). There are many wellknown methods of calculating the band structures of photonic and phononic crystals in addition to the reduced Bloch mode expansion method: the plane-wave expansion (PWE) method (Huang & Wu, 2005; Kushwaha et al., 1993; Laude et al., 2005; Tanaka & Tamura, 1998; Wu et al., 2004 ; Wu & Huang, 2004), the multiple-scattering theory (MST) (Leung & Qiu, 1993; Kafesaki & Economou, 1999; Psarobas & Stefanou, 2000; Wang et al., 1993), the finite-difference (FD) method (Garica-Pabloset et al., 2000; Sun & Wu, 2005; Yang, 1996), the transfer matrix method (Pendry & MacKinnon, 1992), the meshless method (Jun et al., 2003), the multiple multipole method (Moreno et al., 2002), the wavelet method (Checoury & Lourtioz, 2006; Yan & Wang, 2006), the pseudospectral method (Chiang et al., 2007), the finite element method (FEM) (Axmann & Kuchment, 1999; Dobson, 1999; Huang & Chen,

This study develops a thermal homogenization model for an aluminum honeycomb panel using the representative volume element (RVE) concept, considering the orthotropic nature of the structure. The RVE thermal homogenization method is a numerical approach for analyzing heterogeneous materials. It employs a constitutive model based on RVE performance to represent thermal behavior. Effective parameters are determined through averaging techniques, and the finite element method solves the thermal problem, accounting for structure topology and material behavior. The resulting heat conduction problem is solved using the finite element method (FEM) to evaluate the effective thermal characteristics. A 3D RVE is generated based on the honeycomb panel's geometry, evaluating thermal conductivity tensor and describing the medium's thermal performance. Numerical tests validate the model by comparing it with the real honeycomb structure under sinusoidal heat flux. Results show good correlation, with maximum temperatures of 1101.9 °C in the real structure and 1096.4 °C in the medium. The homogeneous medium is further used to investigate thermal performance under convective conditions with varying panel thicknesses, achieving over 77 °C temperature reduction with the thickest panel. Natural vibration behavior is considered, demonstrating strong correlation between modal responses and natural frequencies. This modeling approach efficiently analyzes thermal behavior in large honeycomb structures, reducing computational time significantly.

Discrete element method (DEM) is an interesting alternative to classical approaches as the finite element method (FEM) to simulate homogeneous and heterogeneous materials. Indeed, although DEM was initially developed to simulate granular systems in motion, it also enables to model a continuous medium with the help of a dense granular packing in which the cohesion is introduced between each pair of particles in contact using beam or spring elements. However, among other issues, the local stress field obtained using DEM is heterogeneous, even in the case where it is theoretically homogeneous. In the present contribution, we investigate the stress field distribution in 3D DEM mechanical simulation, and propose an approach we named Halo, to evaluate the level of stress distribution. The idea of the method is to evaluate the stress at the scale of every discrete element (DE) taking into account the contribution of its neighbors, inside the Halo of the DE. The effect of the number of DEs per Halo on the stress distribution is discussed in homogeneous and heterogeneous media. The approach results for homogeneous media are compared to FEM calculations and to the theoretical Hertz solution via a Brazilian test. For the homogeneous media, we use the unidirectional flax/bio based epoxy composite, for which macroscopic longitudinal Young’s modulus is experimentally quantified. Furthermore, comparisons with FEM are performed in this context. Results highlight a good agreement between both approaches in terms of stress fields including extremal values.

Rotational motions in homogeneous anisotropic elastic media are studied under the assumption of plane wave propagation. The main goal is to investigate the influences of anisotropy in the behavior of the rotational wavefield. The focus is on P-waves that theoretically do not generate rotational motion in isotropic media. By using the Kelvin–Christoffel equation, expressions are obtained of the rotational motions of body waves as a function of the propagation direction and the coefficients of the elastic modulus matrix. As a result, the amplitudes of the rotation rates and their radiation patterns are quantified and it is concluded that (1) for strong local earthquakes and typical reservoir situations quasi P-rotation rates induced by anisotropy are significant, recordable, and can be used for inverse problems; and (2) for teleseismic wavefields, anisotropic effects are unlikely to be responsible for the observed rotational energy in the P coda.

Thin inclusions of large or small conductivity in an anisotropic homogeneous medium

Metamaterials are compounds of periodic sub-wavelength elements, whose properties are strongly influenced by the applied external electromagnetic (EM) waves within a typically narrow frequency band. Metamaterial structures could be designed to produce very thin films with high losses manifested by perfect absorption. This work discusses a metamaterial (MTM) based dual band absorber operating at mid infrared (IR) regime. The idea of this absorber originates from combining the absorption responses of the cross shaped and cut wire (CW) resonators. The strong coupling between the two resonators is enhanced by introducing a U-shaped resonator interdigitated between the CW and the cross-shaped resonators. This setup would improve the coupling between the neighboring cells and hence, reduce the overall dimensions. The simulated results are obtained using the full electromagnetic wave simulation package CST MWS. The absorber is designed to possess two distinct absorption bands centered at 9.57μm and 4 μm. The two bands are well separated with a first to second band ratio of ∼ 2.39. The sub-wavelength period size is found to be 0.19X and 0A5X with respect to the first and second bands respectively. High absorption is also observed for wide range of incident angles. Moreover, circuit elements based on circuit approximation are calculated by analyzing the electric and magnetic fields to provide a physical meaningful circuit representation of the electromagnetic response, which can be very helpful in the design and optimization process. The extracted circuit model fully describes and simplifies the behavior of the periodically distributed physical structure and translates it into equivalent discrete entities. Therefore, it offers an approximate prediction for the resonance frequencies. Given that the size of the sub-wavelength elements are much smaller than the incident EM wave operating wavelength, the MTM compound behaves as a homogeneous medium, whose properties can be expressed by effective constitutive parameters. Effective constitutive parameters are fitted using Drude-Lorentz oscillators. This provides insight into the nature of the microscopic behavior of the absorber elements when they are excited by external fields. Finally, samples of this design were fabricated and tested. The measured response shows two absorption peaks with better than 85% absorption for both peaks. This design shows significant improvement over the existing designs, which makes it appealing for various planar applications in the mid-IR frequencies.