Analytical solutions for movement of cold water thermal front in a heterogeneous geothermal reservoir

In the present study an analytical model has been presented to describe the transient temperature distribution and advancement of the thermal front generated due to the reinjection of heat depleted water in a heterogeneous geothermal reservoir. One dimensional heat transport equation in porous media with advection and longitudinal heat conduction has been solved analytically using Laplace transform technique in a semi infinite medium. The heterogeneity of the porous medium is expressed by the spatial variation of the flow velocity and the longitudinal effective thermal conductivity of the medium. A simpler solution is also derived afterwards neglecting the longitudinal conduction depending on the situation where the contribution to the transient heat transport phenomenon in the porous media is negligible. Solution for a homogeneous aquifer with constant values of the rock and fluid parameters is also derived with an aim to compare the results with that of the heterogeneous one. The effect of some of the parameters involved, on the transient heat transport phenomenon is assessed by observing the variation of the results with different magnitudes of those parameters. Results prove the heterogeneity of the medium, the flow velocity and the longitudinal conductivity to have great influence and porosity to have negligible effect on the transient temperature distribution. (C) 2013 Elsevier Inc. All rights reserved.

Majorana edge modes in the Kitaev model

We study the Majorana modes, both equilibrium and Floquet, which can appear at the edges of the Kitaev model on the honeycomb lattice. We first present the analytical solutions known for the equilibrium Majorana edge modes for both zigzag and armchair edges of a semi-infinite Kitaev model and chart the parameter regimes in which they appear. We then examine how edge modes can be generated if the Kitaev coupling on the bonds perpendicular to the edge is varied periodically in time as periodic delta-function kicks. We derive a general condition for the appearance and disappearance of the Floquet edge modes as a function of the drive frequency for a generic d-dimensional integrable system. We confirm this general condition for the Kitaev model with a finite width by mapping it to a one-dimensional model. Our numerical and analytical study of this problem shows that Floquet Majorana modes can appear on some edges in the kicked system even when the corresponding equilibrium Hamiltonian has no Majorana mode solutions on those edges. We support our analytical studies by numerics for a finite sized system which show that periodic kicks can generate modes at the edges and the corners of the lattice.

Traveling waves in trimer granular lattice I: Bifurcation structure of traveling waves in the unit-cell model

Present paper is the first one in the series devoted to the dynamics of traveling waves emerging in the uncompressed, tri-atomic granular crystals. This work is primarily concerned with the dynamics of one-dimensional periodic granular trimer (tri-atomic) chains in the state of acoustic vacuum. Each unit cell consists of three spherical particles of different masses subject to periodic boundary conditions. Hertzian interaction law governs the mutual interaction of these particles. Under the assumption of zero pre-compression, this interaction is modeled as purely nonlinear, which means the absence of linear force component. The dynamics of such chains is governed by the two system parameters that scale the mass ratios between the particles of the unit cell. Such a system supports two different classes of periodic solutions namely the traveling and standing waves. The primary objective of the present study is the numerical analysis of the bifurcation structure of these solutions with emphasis on the dynamics of traveling waves. In fact, understanding of the bifurcation structure of the traveling wave solutions emerging in the unit-cell granular trimer is rather important and can shed light on the more complex nonlinear wave phenomena emerging in semi-infinite trimer chains. (c) 2016 Elsevier B.V. All rights reserved.

Dislocation theory of the fracture criterion for anisotropic solids

A dislocation theory of fracture criterion for the mixed dislocation emission and cleavage process in an anisotropic solid is developed in this paper. The complicated cases involving mixed-mode loading are considered here. The explicit formula for dislocations interaction with a semi-infinite crack is obtained. The governing equation for the critical condition of crack cleavage in an anisotropic solid after a number dislocation emissions is established. The effects of elastic anisotropy, crack geometry and load phase angle on the critical energy release rate and the total number of the emitted dislocations at the onset of cleavage are analysed in detail. The analyses revealed that the critical energy release rates can increase to one or two magnitudes larger than the surface energy because of the dislocation emission. It is also found elastic anisotropy and crystal orientation have significant effects on the critical energy release rates. The anisotropic values can be several times the isotropic value in one crack orientation. The values may be as much as 40% less than the isotropic value in another crack orientation and another anisotropy parameter. Then the theory is applied to a fee single crystal. An edge dislocation can emit from the crack tip along the most highly shear stressed slip plane. Crack cleavage can occur along the most highly stressed slip plane after a number of dislocation emissions. Calculation is carried out step by step. Each step we should judge by which slip system is the most highly shear stressed slip system and which slip system has the largest energy release rate. The calculation clearly shows that the crack orientation and the load phase angle have significant effects on the crystal brittle-ductile behaviours.

The Displacement and Strain Field of Three-Dimensional Rheologic Model of Earthquake Preparation

Based on the three-dimensional elastic inclusion model proposed by Dobrovolskii, we developed a rheological inclusion model to study earthquake preparation processes. By using the Corresponding Principle in the theory of rheologic mechanics, we derived the analytic expressions of viscoelastic displacement U(r, t) , V(r, t) and W(r, t), normal strains epsilon(xx) (r, t), epsilon(yy) (r, t) and epsilon(zz) (r, t) and the bulk strain theta (r, t) at an arbitrary point (x, y, z) in three directions of X axis, Y axis and Z axis produced by a three-dimensional inclusion in the semi-infinite rheologic medium defined by the standard linear rheologic model. Subsequent to the spatial-temporal variation of bulk strain being computed on the ground produced by such a spherical rheologic inclusion, interesting results are obtained, suggesting that the bulk strain produced by a hard inclusion change with time according to three stages (alpha, beta, gamma) with different characteristics, similar to that of geodetic deformation observations, but different with the results of a soft inclusion. These theoretical results can be used to explain the characteristics of spatial-temporal evolution, patterns, quadrant-distribution of earthquake precursors, the changeability, spontaneity and complexity of short-term and imminent-term precursors. It offers a theoretical base to build physical models for earthquake precursors and to predict the earthquakes.

Models for acoustic propagation through turbofan exhaust flows-lined ducts

Turbomachinery noise radiating into the rearward arc is an important problem. This noise is scattered by the trailing edges of the nacelle and the jet exhaust, and interacts with the shear layers between the external flow, bypass stream and jet, en route to the far field. In the past a range of relevant model problems involving semi-infinite cylinders have been solved. However, one limitation of these previous solutions is that they do not allow for the jet nozzle protruding a finite distance beyond the end of the nacelle (or in certain configurations being buried a finite distance upstream). With this in mind, we have used the matrix Wiener-Hopf technique to allow precisely this finite nacelle-jet nozzle separation to be included. We have previously reported results for the case of hard-walled ducts, which requires factorisation of a 2 × 2 matrix. In this paper we extend this work by allowing one of the duct walls, in this case the outer wall of the jet pipe, to be acoustically lined. This results in the need to factorise a 3 × 3 matrix, which is completed by use of a combination of pole-removal and Pad́e approximant techniques. Sample results are presented, investigating in particular the effects of exit plane stagger and liner impedance. Here we take the mean flow to be zero, but extension to nonzero Mach numbers in the core and bypass flow has also been completed. Copyright © 2009 by Nigel Peake & Ben Veitch.

The investigation on dynamic fracture behavour of materials under compressive-shear combined stress waves

In this paper, an improved plate impact experimental technique is presented for studying dynamic fracture mechanism of materials, under the conditions that the impacting loading is provided by a single pulse and the loading time is in the sub-microsecond range. The impacting tests are carried out on the pressure-shear gas gun. The loading rate achieved is dK/dt similar to 10(8) MPa m(1/2) s(-1). With the elimination of influence of the specimen boundary, the plane strain state of a semi-infinite crack in an infinite elastic plate is used to simulate the deformation fields of crack tip. The single pulses are obtained by using the "momentum trap" technique. Therefore, the one-time actions of the single pulse are achieved by eradicating the stress waves reflected from the specimen boundary or diffracted from the crack surfaces. In the current study, some important phenomena have been observed. The special loading of the single pulse can bring about material damage around crack tip, and affect the material behavior, such as kinking and branching of the crack propagation. Failure mode transitions from mode I to mode II crack are observed under asymmetrical impact conditions. The mechanisms of the dynamic crack propagation are consistent with the damage failure model.

Role Of Vertical Component Of Surface Tension Of The Droplet On The Elastic Deformation Of Pdms Membrane

Poly(dimethylsiloxane) (PDMS) has been widely used in lab-on-a-chip and micro- total analysis systems (mu-TAS), thus wetting and electrowetting behaviors of PDMS are of great importance in these devices. PDMS is a kind of soft polymer material, so the elastic deformation of PDMS membrane by a droplet cannot be neglected due to the vertical component of the interfacial tension between the liquid and vapor, and this vertical component of liquid-vapor surface tension is also balanced by the stress distribution within the PDMS membrane. Such elastic deformation and stress distribution not only affect the exact measurement of contact angle, but also have influence on the micro-fluidic behavior of the devices. Using ANSYS code, we simulated numerically the elastic deformation and stress distribution of PDMS membrane on a rigid substrate due to the liquid-vapor surface tension. It is found that the vertical elastic deformation of the PDMS membrane is on the order of several tens of nanometers due to the application of a droplet with a diameter of 2.31 mm, which is no longer negligible for lab-on-a-chip and mu-TAS. The vertical elastic deformation increases with the thickness of the PDMS membrane, and there exists a saturated membrane thickness, regarded as a semi-infinite membrane thickness, and the vertical elastic deformation reaches a limiting value when the membrane thickness is equal to or thicker than such saturated thickness. (C) Koninklijke Brill NV, Leiden, 2008.

On The Thermally Induced Interfacial Delamination Of A Segmented Coating

The thermally induced interfacial delamination problem of a segmented coating is investigated using finite element method (FEM). The coating-substrate system, modeled as a coated semi-infinite medium with periodic segmentation cracks within coating, is assumed to be exposed to convective cooling from surface. The failure criterion based on the interfacial fracture toughness is adopted, in which the energy release rate for an interface crack is considered to be the driving force for interfacial delamination extension. The results confirm that a segmented coating has higher delamination resistance than an intact one under the same thermal transients, as the segmentation crack spacing is smaller than a critical value. Based on dimensional analysis, sensitivity analyses of the crack driving force are also obtained as a function of various dimensionless parameters such as time, convection severity and material constants. These results may provide some helpful references for the integrity of coating-substrate systems under thermal loading. (C) 2007 Elsevier B.V. All rights reserved.

Stress intensity factor histories of a steadily propagating crack under 3-D combined mode traction

Elastodynamic stress intensity factor histories of an unbounded solid containing a semi-infinite plane crack that propagates at a constant velocity under 3-D time-independent combined mode loading are considered. The fundamental solution, which is the response of point loading, is obtained. Then, stress intensity factor histories of a general loading system are written out in terms of superposition integrals. The methods used here are the Laplace transform methods in conjunction with the Wiener-Hopf technique.

A discussion about a class of stress intensity factors and its verification

In this paper, new formulae of a class of stress intensity factors for an infinite plane with two collinear semi-infinite cracks are presented. The formulae differ from those gathered in several handbooks used all over the world. Some experiments and finite element calculations have been developed to verify the new formulae and the results have shown its reliability. Finally, the new formulae and the old are listed to show the differences between them.

The Displacement and Strain Field of Three-Dimensional Rheologic Model of Earthquake Preparation

Based on the three-dimensional elastic inclusion model proposed by Dobrovolskii, we developed a rheological inclusion model to study earthquake preparation processes. By using the Corresponding Principle in the theory of rheologic mechanics, we derived the analytic expressions of viscoelastic displacement U(r, t) , V(r, t) and W(r, t), normal strains epsilon(xx) (r, t), epsilon(yy) (r, t) and epsilon(zz) (r, t) and the bulk strain theta (r, t) at an arbitrary point (x, y, z) in three directions of X axis, Y axis and Z axis produced by a three-dimensional inclusion in the semi-infinite rheologic medium defined by the standard linear rheologic model. Subsequent to the spatial-temporal variation of bulk strain being computed on the ground produced by such a spherical rheologic inclusion, interesting results are obtained, suggesting that the bulk strain produced by a hard inclusion change with time according to three stages (alpha, beta, gamma) with different characteristics, similar to that of geodetic deformation observations, but different with the results of a soft inclusion. These theoretical results can be used to explain the characteristics of spatial-temporal evolution, patterns, quadrant-distribution of earthquake precursors, the changeability, spontaneity and complexity of short-term and imminent-term precursors. It offers a theoretical base to build physical models for earthquake precursors and to predict the earthquakes.

Subsonic interface crack with crack face contact

A steady-state subsonic interface crack propagating between an elastic solid and a rigid substrate with crack face contact is studied. Two cases with respective to the contact length are considered, i.e., semi-infinite and finite crack face contact. Different from a stationary or an open subsonic interface crack, stress singularity at the crack tip in the present paper is found to be non-oscillatory. Furthermore, in the semi-infinite contact case, the singularity of the stress field near the crack tip is less than 1/2. In the finite contact case, no singularity exists near the crack tip, but less than 1/2 singularity does at the end of the contact zone. In both cases, the singularity depends on the linear contact coefficient and the crack speed. Asymptotic solutions near the crack tip are given and analyzed. In order to satisfy the contact conditions, reasonable region of the linear contact coefficient is found. In addition, the solution predicts a non-zero-energy dissipation rate due to crack face contact.

Some approximate solutions of dynamic problems in the linear theory of thin elastic shells

Some aspects of wave propagation in thin elastic shells are considered. The governing equations are derived by a method which makes their relationship to the exact equations of linear elasticity quite clear. Finite wave propagation speeds are ensured by the inclusion of the appropriate physical effects.

The problem of a constant pressure front moving with constant velocity along a semi-infinite circular cylindrical shell is studied. The behavior of the solution immediately under the leading wave is found, as well as the short time solution behind the characteristic wavefronts. The main long time disturbance is found to travel with the velocity of very long longitudinal waves in a bar and an expression for this part of the solution is given.

When a constant moment is applied to the lip of an open spherical shell, there is an interesting effect due to the focusing of the waves. This phenomenon is studied and an expression is derived for the wavefront behavior for the first passage of the leading wave and its first reflection.

For the two problems mentioned, the method used involves reducing the governing partial differential equations to ordinary differential equations by means of a Laplace transform in time. The information sought is then extracted by doing the appropriate asymptotic expansion with the Laplace variable as parameter.

A study of nonlinear phenomena in the propagation of electromagnetic waves in a weakly ionized gas

This thesis is a study of nonlinear phenomena in the propagation of electromagnetic waves in a weakly ionized gas externally biased with a magnetostatic field. The present study is restricted to the nonlinear phenomena rising from the interaction of electromagnetic waves in the ionized gas. The important effects of nonlinearity are wave-form distortion leads to cross modulation of one wave by a second amplitude-modulated wave.

The nonlinear effects are assumed to be small so that a perturbation method can be used. Boltzmann’s kinetic equation with an appropriate expression for the collision term is solved by expanding the electron distribution function into spherical harmonics in velocity space. In turn, the electron convection current density and the conductivity tensors of the nonlinear ionized gas are found from the distribution function. Finally, the expression for the current density and Maxwell’s equations are employed to investigate the effects of nonlinearity on the propagation of electromagnetic waves in the ionized gas, and also on the reflection of waves from an ionized gas of semi-infinite extent.