Dynamic anti-plane behavior of rare earth giant magnetostrictive medium with a circular cavity defect in semi-space

An analytical solution to the anti-plane dynamics problem of semi-space rare earth giant magnetostrictive media with circular cavity defects near the horizontal boundary under the action of SH wave is studied. Based on the Helmholtz theorem and the theory of complex function, the elastic-magnetic dynamic equation of magnetostrictive medium is established, and the semi-space incident wave field is written. In addition, based on the theory of complex function and the method of wave function expansion, the expression of the wave function of the scattered displacement field and the corresponding magnetic potential of the scattered wave under the condition of no stress and magnetic insulation of the horizontal boundary are obtained. Then, based on the conditions of free boundary stress, continuous magnetic induction intensity and continuous magnetic potential around the circular cavity, the infinite linear algebraic equations are established. Finally, the analytical expressions of dynamic stress concentration factor and magnetic field intensity concentration factor around circular cavity in semi-space rare earth giant magnetostrictive medium are obtained. Numerical examples show that the analysis results depend on the following parameters: permeability, dimensional-piezomagnetic coefficient, frequency of the incident wave, incident angle, distance between the circular cavity and horizontal boundary. These results have certain reference value for the study of non-destructive testing and failure analysis of rare earth giant magnetostrictive materials.

Dynamic Anti-plane Behavior of Rare Earth Giant Magnetostrictive Medium With a Circular Cavity Defect in Semi-space

An analytical solution to the anti-plane dynamics problem of semi-space rare earth giant magnetostrictive media with circular cavity defects near the horizontal boundary under the action of SH wave is studied. Based on the Helmholtz theorem and the theory of complex function, the elastic-magnetic dynamic equation of magnetostrictive medium is established, and the semi-space incident wave field is written. In addition, the scattered displacement field and the corresponding magnetic potential of the scattered wave under the condition of no stress and magnetic insulation of the horizontal boundary are obtained by image method. Then, based on the conditions of free boundary stress, continuous magnetic induction intensity and continuous magnetic potential around the circular cavity, the infinite linear algebraic equations are established. Finally, the analytical expressions of dynamic stress concentration factor and magnetic field intensity concentration factor around circular cavity in semi-space rare earth giant magnetostrictive medium are obtained. Numerical examples show that the analysis results depend on the following parameters: permeability, dimensional-piezomagnetic coefficient, frequency of the incident wave, incident angle, distance between the circular cavity and horizontal boundary. These results have certain reference value for the study of non-destructive testing and failure analysis of rare earth giant magnetostrictive materials.

Stress and deformation mechanisms at a subduction zone: insights from 2-D thermomechanical numerical modelling

SUMMARY Numerous processes such as metamorphic reactions, fluid and melt transfer and earthquakes occur at a subducting zone, but are still incompletely understood. These processes are affected, or even controlled, by the magnitude and distribution of stress and deformation mechanism. To eventually understand subduction zone processes, we quantify here stresses and deformation mechanisms in and around a subducting lithosphere, surrounded by asthenosphere and overlain by an overriding plate. We use 2-D thermomechanical numerical simulations based on the finite difference and marker-in-cell method and consider a 3200 km wide and 660 km deep numerical domain with a resolution of 1 km by 1 km. We apply a combined visco-elasto-plastic deformation behaviour using a linear combination of diffusion creep, dislocation creep and Peierls creep for the viscous deformation. We consider two end-member subduction scenarios: forced and free subduction. In the forced scenario, horizontal velocities are applied to the lateral boundaries of the plates during the entire simulation. In the free scenario, we set the horizontal boundary velocities to zero once the subducted slab is long enough to generate a slab pull force large enough to maintain subduction without horizontal boundary velocities. A slab pull of at least 1.8 TN m–1 is required to continue subduction in the free scenario. We also quantify along-profile variations of gravitational potential energy (GPE). We evaluate the contributions of topography and density variations to GPE variations across a subduction system. The GPE variations indicate large-scale horizontal compressive forces around the trench region and extension forces on both sides of the trench region. Corresponding vertically averaged differential stresses are between 120 and 170 MPa. Furthermore, we calculate the distribution of the dominant deformation mechanisms. Elastoplastic deformation is the dominant mechanism in the upper region of the lithosphere and subducting slab (from ca. 5 to 60 km depth from the top of the slab). Viscous deformation dominates in the lower region of the lithosphere and in the asthenosphere. Considering elasticity in the calculations has an important impact on the magnitude and distribution of deviatoric stress; hence, simulations with increased shear modulus, in order to reduce elasticity, exhibit considerably different stress fields. Limiting absolute stress magnitudes by decreasing the internal friction angle causes slab detachment so that slab pull cannot be transmitted anymore to the horizontal lithosphere. Applying different boundary conditions shows that forced subduction simulations are stronger affected by the applied boundary conditions than free subduction simulations. We also compare our modelled topography and gravity anomaly with natural data of seafloor bathymetry and free-air gravity anomalies across the Mariana trench. Elasticity and deviatoric stress magnitudes of several hundreds of MPa are required to best fit the natural data. This agreement suggests that the modelled flexural behaviour and density field are compatible with natural data. Moreover, we discuss potential applications of our results to the depth of faulting in a subducting plate and to the generation of petit-spot volcanoes.

Rock Stress Around Noncircular Tunnel: a New Simple Mathematical Method

A new simple mathematical method has been proposed to predict rock stress around a noncircular tunnel and the method is calibrated and validated with a numerical model. It can be found that the tunnel shapes and polar angles affect the applicable zone of the theoretical model significantly and the applicable zone of a rectangular tunnel was obtained using this method. The method can be used to predict the values of the concentrated stress, and to analyze the change rate of rock stress and back to calculate the mechanical boundary condition in the applicable zone. The results of the stress change rate indicate that the horizontal stress is negatively related to the vertical boundary load and positively related to the horizontal boundary load. The vertical stress is negatively related to the horizontal boundary load and positively related to the vertical boundary load. These findings can be used to explain the evolution of the vertical increment in stress obtained with field-based borehole stress monitoring.

Turbulent horizontal convection under spatially periodic forcing: a regime governed by interior inertia

Differential heating applied at a single horizontal boundary forces ‘horizontal convection’, even when there is no net heat flux through the boundary. However, almost all studies of horizontal convection have been limited to a special class of problem in which temperature or heat flux differences were applied in only one direction and over the horizontal length of a box (the Rossby problem; Rossby, Deep-Sea Res., vol. 12, 1965, pp. 9–16). These conditions strongly constrain the flow. Here we report laboratory experiments and direct numerical simulations (DNS) extending the results of Griffiths & Gayen (Phys. Rev. Lett., vol. 115, 2015, 204301) for horizontal convection forced by boundary conditions imposed in a two-dimensional periodic array at a horizontal boundary. The experiments use saline and freshwater fluxes at a permeable base with the imposed boundary salinity having a horizontal length scale one quarter of the width of the box. The flow reaches a state in which the net boundary buoyancy flux vanishes and the bulk of the fluid shows an inertial range of turbulence length scales. A regime transition is seen for increasing water depth, from an array of individual coherent plumes on the forcing scale to convection dominated by emergent larger scales of overturning. The DNS explore the analogous thermally forced case with sinusoidal boundary temperature of wavenumber $n=4$, and are used to examine the Rayleigh number ($Ra$) dependence for shallow- and deep-water cases. For shallow water the flow transitions with increasing $Ra$ from laminar to turbulent boundary layer regimes that are familiar from the Rossby problem and which have normalised heat transport scaling as $Nu\sim Ra^{1/5}$ and $Nu\sim (Ra\,Pr)^{1/5}$, with $Nu$ the Nusselt number and $Pr$ the Prandtl number, in this case maintaining a stable array of coherent turbulent plumes. For deep-water and large $Ra$ the laminar scaling transitions to $Nu\sim (Ra\,Pr)^{1/4}$, with the scales of turbulence extending to the dimensions of the box. The $1/4$ power law regime is explained in terms of the momentum of symmetric, inviscid large scales of motion in the interior coupled to diffusive loss of heat through stabilised parts of the boundary layer. The turbulence production is predominantly by shear instability rather than convection, with viscous dissipation distributed throughout the bulk of the fluid. These conditions are not seen in the highly asymmetric flow in the Rossby problem even at Rayleigh numbers up to six orders of magnitude greater than the transition found here. The new inertial interior regime has the rate of supply of available potential energy, and its removal by mixing of density, increasing as $Ra^{5/4}$, which is faster than $Ra^{6/5}$ in the Rossby problem. Irreversible mixing is confined close to the forcing boundary and is very much larger than the viscous dissipation, which is proportional to $Ra$.