Numerical simulation of wave propagation in a porous elastic solid using biot theory

1989 ◽  
Author(s):  
Xianhuai Zhu ◽  
George A. McMechan
Geophysics ◽  
2021 ◽  
pp. 1-39
Author(s):  
Josué G. Gonzalez ◽  
Pratap N. Sahay ◽  
Tobias M. Mueller

The shear motion in Newtonian fluids, i.e., the fluid vorticity, represents an intrinsic loss mechanism governed by a diffusion equation. Its description involves the trace-free part of the fluid viscous stress tensor. This part is missing in the Biot theory of poroelasticity. As a result, the fluid vorticity is not captured, and only one S-wave is predicted. The missing fluid vorticity has implications for the propagation of S-waves across discontinuities. This becomes most apparent in the problem of S-wave propagation across the welded contact of an elastic solid with a porous medium. At such a contact, the no-slip condition between the elastic solid and the constituent parts of the porous medium, the solid-frame, and the pore-fluid, must hold. This requirement translates into a vanishing relative motion of the fluid with respect to the solid-frame, i.e., filtration field, at the contact. Nevertheless, our analysis shows that for the Biot theory, in the low-frequency regime, a non-zero, although insignificantly small filtration field exists at the contact. But, more importantly, the filtration field is noticeable when the transition to the high-frequency regime occurs. This constitutes a disagreement with the requirement of a no-slip boundary condition and renders the prediction unphysical. This shortcoming is circumvented by including the fluid viscous stress tensor into the poroelastic constitutive relations, as stipulated by the de la Cruz-Spanos poroelasticity theory. Then, a second S-wave is predicted which manifests as the fluid vorticity at macroscale. This process is distinct from the fast S-wave, the other predicted S-wave akin to the Biot S-wave. We find that the generation of this process at the contact induces a filtration field equal and opposite to that associated with the fast S-wave. Therefore, the no-slip condition is satisfied, and the S-wave reflection/transmission across a discontinuity becomes physically meaningful.


Materials ◽  
2021 ◽  
Vol 14 (13) ◽  
pp. 3553
Author(s):  
Dengwang Wang ◽  
Yong Gao ◽  
Sheng Wang ◽  
Jie Wang ◽  
Haipeng Li

Carbon/Phenolic (C/P), a typical anisotropic material, is an important component of aerospace and often used to protect the thermodynamic effects of strong X-ray radiation. In this paper, we establish the anisotropic elastic-plastic constitutive model, which is embedded in the in-house code “RAMA” to simulate a two-dimensional thermal shock wave induced by X-ray. Then, we compare the numerical simulation results with the thermal shock wave stress generated by the same strong current electron beam via experiment to verify the correctness of the numerical simulation. Subsequently, we discuss and analyze the rules of thermal shock wave propagation in C/P material by further numerical simulation. The results reveal that the thermal shock wave represents different shapes and mechanisms by the radiation of 1 keV and 3 keV X-rays. The vaporization recoil phenomenon appears as a compression wave under 1 keV X-ray irradiation, and X-ray penetration is caused by thermal deformation under 3 keV X-ray irradiation. The thermal shock wave propagation exhibits two-dimensional characteristics, the energy deposition of 1 keV and 3 keV both decays exponentially, the energy deposition of 1 keV-peak soft X-ray is high, and the deposition depth is shallow, while the energy deposition of 3 keV-peak hard X-ray is low, and the deposition depth is deep. RAMA can successfully realize two-dimensional orthotropic elastoplastic constitutive relation, the corresponding program was designed and checked, and the calculation results for inspection are consistent with the theory. This study has great significance in the evaluation of anisotropic material protection under the radiation of intense X-rays.


Author(s):  
Z. Y. Song ◽  
C. Cheng ◽  
F. M. Xu ◽  
J. Kong

Based on the analytical solution of one-dimensional simplified equation of damping tidal wave and Heuristic stability analysis, the precision of numerical solution, computational time and the relationship between the numerical dissipation and the friction dissipation are discussed with different numerical schemes in this paper. The results show that (1) when Courant number is less than unity, the explicit solution of tidal wave propagation has higher precision and requires less computational time than the implicit one; (2) large time step is allowed in the implicit scheme in order to reduce the computational time, but the precision of the solution also reduce and the calculation precision should be guaranteed by reducing the friction factor: (3) the friction factor in the implicit solution is related to Courant number, presented as the determined friction factor is smaller than the natural value when Courant number is larger than unity, and their relationship formula is given from the theoretical analysis and the numerical experiments. These results have important application value for the numerical simulation of the tidal wave.


1974 ◽  
Vol 41 (1) ◽  
pp. 83-88 ◽  
Author(s):  
M. Kurashige

A study is made of the radial propagation of axial shear waves in an incompressible elastic solid under finite radial deformation. Basic equations are derived on the basis of Biot’s mechanics of incremental deformations, and analysis is made by the method of characteristics. Numerical examples are given by specializing the initial deformation to two cases: (a) an infinite solid with a cylindrical bore is inflated by all-around tension, and (b) a cuboid is rounded into a ring and its ends are bonded to each other. The influence of the inhomogeneities of such deformations upon the laws of shear wave propagation is presented in the form of curves.


Sign in / Sign up

Export Citation Format

Share Document