共查询到20条相似文献，搜索用时 46 毫秒
1.
Steven D. London 《地球物理与天体物理流体动力学》2013,107(4):317-333
We consider an electrically conducting fluid in rotating cylindrical coordinates in which the Elsasser and magnetic Reynolds numbers are assumed to be large while the Rossby number is assumed to vanish in an appropriate limit. This may be taken as a simple model for the Earth's outer core. Fully nonlinear waves dominated by the nonlinear Lorentz forces are studied using the method of geometric optics (essentially WKB). These waves are assumed to be of the form of an asymptotic series expanded about ambient magnetic and velocity fields which vanish on the equatorial plane. They take the form of short wave, slowly varying wave trains. The first-order approximation is sinusoidal and basically the same as in the linear problem, with a dispersion relation modified by the appearance of mean terms. These mean terms, as well the undetermined amplitude functions, are found by suppressing secular terms in a “fast” variable in the second-order approximation. The interaction of the mean terms with the dispersion relation is the primary cause of behaviors which differ from the linear case. In particular, new singularities appear in the wave amplitude functions and an initial value problem results in a singularity in one of the mean terms which propagates through the fluid. The singularities corresponding to the linear ones are shown to develop when the corresponding waves propagate toward the equatorial plane. 相似文献
2.
S. D. London 《地球物理与天体物理流体动力学》2013,107(3-4):251-268
Abstract Numerical work indicates that resistive instability may be the dominant mode of instability in the Earth's outer core for realistic core parameter regimes. In this paper, we assume that the Elsasser number is large in order to obtain an asymptotic analysis of resistive instability in an electrically conducting fluid confined to a rotating cylindrical shell of infinite extent in the axial direction. The dimensionless equations of motion are linearized about an ambient magnetic field which is purely azimuthal and depends only on the cylindrical radial variable. Applying the theory of ordinary differential equations with a large parameter, we obtain an asymptotic approximation to the solution. Relatively simple analytic expressions for the complex frequencies are obtained by applying the boundary conditions for insulating boundaries at the cylindrical sidewalls and then assuming that the ambient magnetic field vanishes at one or both of those sidewalls. The results appear to be consistent with previous numerical work. 相似文献
3.
David P Thambiratnam 《地震工程与结构动力学》1986,14(3):475-485
Herein the propagation of transient waves in a rod due to impulsive end loading is investigated using the method of wavefront expansion. The Mindlin–Herrmann equations which include the Poisson effect are used to model the rod. Since this rod theory is restricted to lower values of Poisson's ratio (< 0.283), the present work also has the same restriction. The impulsive end loading on the rod can be prescribed in the form of stress, strain, velocity or acceleration boundary conditions. The analysis is based on the concept of a wave as a carrier of discontinuities in the field variables and their derivatives. These discontinuities are determined from a set of recurrence relations which are in turn generated by the use of time harmonic asymptotic series solutions to the equations of motion. Transients due to both velocity and strain (or stress) boundary conditions are treated. Numerical examples are presented to illustrate the method of solution. The results indicate the distortion of the transient waves due to the dispersive nature of the governing equations. Transients due to other boundary conditions at the rod end can be treated by using the results presented in the paper. 相似文献
4.
Thermal instabilities in the form of oscillatory magnetoconvection representing diffusively modified Alfvén waves in an electrically-conducting Bénard fluid layer of rigid walls in the presence of a vertical magnetic field are investigated. Emphasis of the article is on the transition from a nearly undamped Alfvén wave to diffusively modified Alfvén waves, and on the effect of physically realisable magnetic field boundary conditions on magnetoconvection. It is found that the extra magnetic dissipation in the magnetic Hartmann boundary layers can enhance oscillatory magnetoconvection in the form of strongly modified Alfvén waves. Oscillatory magnetoconvection produced solely by the Alfvén wave mechanism can be the most unstable mode even in the presence of a strong viscous effect. This article also represents the first study on the effect of an electrically conducting wall on magnetoconvection which is associated with a nonlinear eigenvalue problem. We find that the electrically perfectly conducting condition does not yield a good approximation for magnetoconvection with an electrically highly conducting wall. The size of oscillation frequency with an electrically highly conducting wall can be more than a factor of 2 larger than that obtained using the perfectly conducting condition. 相似文献
5.
An absorbing boundary for saturated porous media is developed that can be used for transient analyses in the time domain. The elastic constitutive equations for the saturated porous media follow Bowen's formulation. The method consists of applying viscous tractions along the artificial boundary. The absorbing boundary behaviour is assumed linear and isotropic. Hadamard's conditions provide the speeds of the dilatational and shear waves that propagate in saturated porous media. Since these expressions are frequency independent, the intensities of the viscous tractions are evaluated in the time domain, and the two dilatational waves are accounted for. The viscous tractions are defined from the drained characteristics, assuming an infinite permeability, at variance with the traditional ‘undrained’ method based on undrained characteristics and a null permeability. Solid media and materials with low permeability are also retrieved as subcases. The results show that, at no additional cost, this ‘drained’ method is more accurate for all permeabilities than the ‘undrained’ method, which disregards the existence of the second dilatational wave. Copyright © 2003 John Wiley & Sons, Ltd. 相似文献
6.
7.
Scattering of elastic waves by three-dimensional canyons embedded within an elastic half-space is investigated by using a wave function expansion technique. The geometry of the canyon is assumed to be non-axisymmetric. The canyon is subjected to incident plane Rayleigh waves and oblique incident SH, SV and P waves. The unknown scattered wavefield is expressed in terms of spherical wave functions which satisfy the equations of motion and radiation conditions at infinity, but they do not satisfy stress-free boundary conditions at the half-space surface. The boundary conditions are imposed locally in the least-squares sense at several points on the surface of the canyon and the half-space. Through a comparative study the validity and limitations of two-dimensional approximations (antiplane strain and plane strain models) have been examined. It is shown that scattering of waves by three-dimensional canyons may cause substantial change in the surface displacement patterns in comparison to the two-dimensional models. These results emphasize the need for three-dimensional modelling of realistic problems of interest in strong ground motion seismology and earthquake engineering. 相似文献
8.
In this work, we develop the indirect boundary element method (IBEM) to simulate the seismic site response in a realistic, large-scale 3-D sedimentary basin. Most previous applications of boundary element method have used full-space Green’s functions for wave propagation between element points. We use half-space Green’s functions, which include the seismic wavefield interactions at the free surface and require only the boundary elements of the basin interface. In this way, the size of the matrix equation for solution in the IBEM can be reduced to approximately a quarter of that using full-space Green’s functions. The site response modeling of the Granada basin in southern Spain using the IBEM shows that the basin-induced scattering waves were identified as propagating back and forth inside the basin. The scattered waves also generate surface waves that are weakly propagated outside of the basin. The wave propagation inside and outside of the basin shows different patterns. We observe that the scattered wave is locally amplified, and its propagation direction deviates from that of the incident waves propagation direction. Therefore, the computed seismic response in the basin could provide us with good estimates of the seismic motion. 相似文献
9.
Two-layer equatorial primitive equations for the free troposphere in the presence of a thin atmospheric boundary layer and thermal dissipation are developed here. An asymptotic theory for the resonant nonlinear interaction of long equatorial baroclinic and barotropic Rossby waves is derived in the presence of such dissipation. In this model, a self-consistent asymptotic derivation establishes that boundary layer flows are generated by meridional pressure gradients in the lower troposphere and give rise to degenerate equatorial Ekman friction. That is to say, the asymptotic model has the property that the dissipation matrix has one eigenvalue which is nearly zero: therefore the dynamics rapidly dissipates flows with pressure at the base of the troposphere and creates barotropic/baroclinic spin up/spin down. The simplified asymptotic equations for the amplitudes of the dissipative equatorial barotropic and baroclinic waves are studied by linear theory and integrated numerically. The results indicate that although the dissipation slightly weakens the tropics to midlatitude connection, strong localized wave packets are nonetheless able to exchange energy between barotropic and baroclinic waves on intraseasonal timescales in the presence of baroclinic mean shear. Interesting dissipation balanced wave-mean flow states are discovered through numerical simulations. In general, the boundary layer dissipation is very efficient for flows in which the barotropic and baroclinic components are of the same sign at the base of the free troposphere whereas the boundary layer dissipation is less efficient for flows whose barotropic and baroclinic components are of opposite sign at the base of the free troposphere. 相似文献
10.
M. R. Sarangi 《地球物理与天体物理流体动力学》2013,107(3):303-325
ABSTRACTA hydro-elastic frame has been considered to investigate the proliferation of waves over small base deformation on an infinitely extended flexible seabed. The flexible base surface is assumed as a thin elastic plate of very small thickness and it depends on the Euler–Bernoulli beam equation. For any particular frequency, there are two different modes of time-harmonic propagating wave exists rather than one mode of propagating wave along the positive horizontal direction. The waves with smaller wavenumber spread along the free-surface of the sea (say, free-surface mode) and the waves with higher wavenumber spread along the flexible base surface (say, flexural mode). A simplified perturbation approach is utilised to bring down the entire equations which govern the original boundary value problem (bvp) to a less complex bvp for the first-order velocity potential function. The first-order potential function along with the first-order reflection and transmission coefficients for both modes are calculated by a procedure based upon Fourier transform approach. A shape of sinusoidal swells flexible base surface is taken as an example to approve the scientific results. It is observed that when the train of normal incident propagating wave spreads over base distortion because of either the free-surface unsettling influence or the flexural wave movement in the sea, the reflected and transmitted energy are always feasible to be exchanged from one particular wave mode to another wave mode. Furthermore, we notice that the realistic changes in the flexural rigidity behaviour on the flexible base surface of the sea have a significant effect on the problem of water wave proliferation over small base deformation. Moreover, the energy conservation equation is derived with the help of the Green's integral theorem. The results for the values of reflection and transmission coefficients obtained for both the free-surface unsettling influence as well as flexural wave movement in the fluid are found to satisfy the energy conservation equation almost accurately. 相似文献
11.
A velocity formula is proposed for flow over a mobile sediment bed induced by velocity-skewed waves and current. The formula is obtained by a separation of waves and current velocities and requires seven free variables related to free stream velocity and sediment characteristics. The formula includes two parts:(1) a wave part consisting of the free stream velocity and defect function, which considers phase lead, wave boundary layer thickness, and mobile bed level, and(2) a current part, which ch... 相似文献
12.
Steven D. London 《地球物理与天体物理流体动力学》2013,107(5):397-411
The equations for an electrically conducting fluid in cylindrical coordinates are linearized assuming that the inertial terms in the momentum equation can be ignored (small Rossby number), and that the ratio of the Elsasser number and magnetic Reynolds number is one. After these assumptions, the governing equations are linearized about an ambient solution which vanishes at the the equator. Upon assuming large Elsasser and magnetic Reynolds number, the solutions to the linearized equations are approximated by wave trains having very short wave length (relative to the core radius) but which vary slowly (on a scale of the core radius). The period of the waves is much longer than a day but much shorter than the period of the slow hydromagnetic oscillations. These waves are found to be trapped in a region about the equator and away from the axis of rotation. The waves break at a latitudinal wave region boundary, in the sense that the waves become exponentially large in a boundary layer, having as an exponent some positive power of the large azimuthal wave number. This behavior is amplified as the Elsasser number becomes smaller while still remaining relatively large. Waves in more Earth-like parameter regimes are discussed briefly. 相似文献
13.
Diffraction of plane P waves by a canyon of arbitrary shape in poroelastic half-space (Ⅰ): Formulation 
下载免费PDF全文
下载免费PDF全文 This paper presents an indirect boundary integration equation method for diffraction of plane P waves by a two-dimensional canyon of arbitrary shape in poroelastic half-space. The Green's functions of compressional and shear wave sources in poroelastic half-space are derived based on Biot's theory. The scattered waves are constructed using the fictitious wave sources close to the boundary of the canyon, and magnitude of the fictitious wave sources are determined by the boundary conditions. The precision of the method is verified by the satisfaction extent of boundary conditions, the comparison between the degenerated solutions of single-phased half-space and the well-known solutions, and the numerical stability of the method. 相似文献
14.
Effect of underground cavities on surface earthquake ground motion under SH wave propagation 总被引：2，自引：0，他引：2
A theoretical approach is presented to study the antiplane seismic response of underground structures, subjected to the incidence of both plane and cylindrical waves. The structure is assumed to be a circular inclusion embedded in a homogenous, isotropic and linear visco‐elastic halfspace. The inclusion may consist either of a cavity, with or without a ring‐shaped boundary, or it may be filled in with a linear‐elastic material, without loss of generality. The analytical solution is obtained using expansions of wave functions in terms of Bessel and Hankel functions, relying on the technique of images and the use of Graf's addition theorem to enforce the boundary conditions. The effects of underground cavities on surface earthquake ground motion are studied as a function of the size of the cavity, its embedment depth, the frequency content of the excitation, the incidence angle and the distance from the axis of symmetry of the cavity itself. A simple application of Rayleigh's method allows us to verify that the ground surface response is dominated by the fundamental vibration mode of the portion of soil between the cavity and ground surface itself, in the frequency range of interest for engineering purposes. A simple relationship to estimate the fundamental natural frequency as a function of the embedment depth of the cavity is given. Finally, amplification factors on response spectra are obtained, to provide a practical insight into the effect of an underground cavity on surface ground motion during real earthquakes. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献
15.
16.
S. Mohapatra 《地球物理与天体物理流体动力学》2013,107(5):347-374
Scattering of oblique incident waves by small bottom undulation in a two-layer fluid, where the upper layer has a thin ice-cover while the lower one has the undulation, is investigated within the framework of linearized water wave theory. The ice-cover is being modeled as an elastic plate of very small thickness. There exist two modes of time-harmonic waves–one with lower wave number propagating along the ice-cover (ice-cover mode) and the other with higher wave number along the interface (interfacial mode). A perturbation analysis is employed to solve the corresponding boundary value problem governed by modified Helmholtz equation and thereby evaluating the reflection and transmission coefficients approximately up to first order for both modes. A patch of sinusoidal ripples, having two different wave numbers over two consecutive stretches, is considered as an example and the related coefficients are determined. It is observed that when the wave is incident on the ice-cover surface we always find energy transfer to the interface, but for interfacial incident waves there are parameter ranges for which no energy transfer to the ice-cover surface is possible. Also it is observed that for small angles of incidence, the reflected energy is more as compared to the other angles of incidence. These results are demonstrated in graphical form. From the derived results, the solutions for problems with free surface can be obtained as particular cases. 相似文献
17.
In this work, a hybrid boundary integral equation method (BIEM) is developed, based on both displacement and hypersingular traction formulations, for the analysis of time-harmonic seismic waves propagating through cracked, multi-layered geological regions with surface topography and under plane strain conditions. Specifically, the displacement-based BIEM is used for a multi-layered deposit with interface cracks, while the regularized, traction-based BIEM is used when internal cracks are present within the layers. The standard uni-dimensional boundary element with parabolic shape functions is employed for discretizing the free surface and the layer interfaces, while special discontinuous boundary elements are placed near the crack tips to model the asymptotic behaviour of both displacements and tractions. This formulation yields displacement amplitudes and phase angles on the free surface of a geological deposit, as well as stress intensity factors near the tips of the cracks. Finally, in the companion paper, numerical results are presented which show that both scattered wave and stress concentration fields are sensitive to the incidence seismic wave parameters and to specific site conditions such as surface topography, layering, the presence of cracks and crack interaction. 相似文献
18.
Scattering of incident plane harmonic pseudo P‐, SH‐, and SV‐waves by a two‐dimensional basin of arbitrary shape is investigated by using an indirect boundary integral equation approach. The basin and surrounding half‐space are assumed to be generally anisotropic, homogeneous, linearly elastic solids. No material symmetries are assumed. The unknown scattered waves are expressed as linear combinations of full‐space time‐harmonic two‐dimensional Green functions. Using the Radon transform, the Green functions are obtained in the form of finite integrals over a unit circle. An algorithm for the accurate and efficient numerical evaluation of the Green functions is discussed. A detailed convergence and parametric analysis of the problem is presented. Excellent agreement is obtained with isotropic results available in the literature. Steady‐state surface ground motion is presented for semi‐circular basins with generally anisotropic material properties. The results show that surface motion strongly depends upon the material properties of the basin as well as the angle of incidence and frequency of the incident wave. Significant mode conversion can be observed for general triclinic materials which are not present in isotropic models. Comparison with an isotropic basin response demonstrates that anisotropy is very important for assessing the nature of surface motion atop basins. Copyright © 2001 John Wiley & Sons, Ltd. 相似文献
19.
In the present paper, effects of non-uniform excitation due to spatially variation of seismic waves under the reservoir bottom on linear and nonlinear responses of arch dams are studied. Foundation is assumed to be massed and infinite elements are utilized to model semi-infinite medium via the far-end boundary of the foundation FE model. A continuum crack propagation model based on the smeared crack approach in 3D space is introduced. Reservoir's water is assumed compressible and the coupled system is solved using the staggered method. As a case study, Amir-Kabir double curvature arch dam in Iran is selected to investigate seismic behavior of the system. Two cases are analyzed in which wave travelling velocities are taken as 650 m/s and infinity. According to the results, non-uniform excitation leads to reduction in seismic response of the system and, in addition, frequency content of the response is different completely in comparison with the system under the uniform excitation. 相似文献
20.
This article commences by surveying the basic dynamics of Earth's core and their impact on various mechanisms of core-mantle coupling. The physics governing core convection and magnetic field production in the Earth is briefly reviewed. Convection is taken to be a small perturbation from a hydrostatic, “adiabatic reference state” of uniform composition and specific entropy, in which thermodynamic variables depend only on the gravitational potential. The four principal processes coupling the rotation of the mantle to the rotations of the inner and outer cores are analyzed: viscosity, topography, gravity and magnetic field. The gravitational potential of density anomalies in the mantle and inner core creates density differences in the fluid core that greatly exceed those associated with convection. The implications of the resulting “adiabatic torques” on topographic and gravitational coupling are considered. A new approach to the gravitational interaction between the inner core and the mantle, and the associated gravitational oscillations, is presented. Magnetic coupling through torsional waves is studied. A fresh analysis of torsional waves identifies new terms previously overlooked. The magnetic boundary layer on the core-mantle boundary is studied and shown to attenuate the waves significantly. It also hosts relatively high speed flows that influence the angular momentum budget. The magnetic coupling of the solid core to fluid in the tangent cylinder is investigated. Four technical appendices derive, and present solutions of, the torsional wave equation, analyze the associated magnetic boundary layers at the top and bottom of the fluid core, and consider gravitational and magnetic coupling from a more general standpoint. A fifth presents a simple model of the adiabatic reference state. 相似文献