Electromagnetic scattering by a perfectly conducting half-plane in a conductive half-space

FollowingDmitriev (1960) a rigorous theoretical solution for the problem of scattering by a perfectly conducting inclined half-plane buried in a uniform conductive half-space has been obtained for plane wave excitation. The resultant integral equation for the Laplace transform of scattering current in the half-plane is solved numerically by the method of successive approximation. The scattered fields at the surface of the half-space are found by integrating the half-space Green's function over the transform of the scattering current.The effects of depth of burial and inclination, of the half-plane on the scattered fields are studied in detail. An increase in the depth of burial leads to attenuation of the fields. Inclination introduces asymmetry in the field profiles beside affecting its magnitude. Depth of exploration is greater for quadrature component. An interpretation scheme based on a phasor diagram is presented for the VLF-EM method of exploration for rich vein deposits in a conductive terrain.List of symbols x, y, z Space co-ordinates - Half-space conductivity - ₀ Free-space permeability - Excitation frequency (angular) - T Time - h Depth of the half-plane - a Inclination of the half-plane - E _x x-Directed total electric field - E _x ^p x-Directed primary electric field - E _xo ^p x-Directed primary electric field atz=0 directly over the half-plane - H _y y-Component of total magnetic field - H _y ^p y-Component of primary magnetic field - H _y0 ^p y-Component of primary magnetic field atz=0 directly over the half-plane - H _z z-Component of total magnetic field - H _z ^p z-Component of primary magnetic field - J _x Surface density ofx-directed scattering current - G Green's function - k ₀,K Wave numbers - u,u ₀,u ₁,u ₂ Functions - Space co-ordinate - s Variable in transform domain - Variable of integration - Normalized scattering current - Laplace transform of - N Normalized - , ₀, ₁, ₂ Functions - t Variable of integration - Skin depth - H Total magnetic field - H ^p Primary magnetic field - H ₀ ^p Primary magnetic field atz=0 directly over the half-plane - M,Q,R,S,U,V Functions - N ₁,N ₂ Functions     

Source parameters of some large earthquakes in the eastern Mediterranean region based on an itertive maximum entropy technique

TeleseismicP-waves of some large earthquakes that occurred in the eastern Mediterranean region have been analysed by using an iterative maximum entropy technique in order to obtain the independent spectral parameters, the long-period spectral level ₀ and the corner frequencyf ₀ of the far-field displacement spectra.Based on these parameters, the seismic source parameters seismic momentM ₀, source dimensionr, fault lengthl, average displacement u, shear stress drop , radiated energyE _s and apparent stressn are calculated for the considered earthquakes by using Brune's and Madariaga's models.The striking feature of the source parameters obtained in this study is the low stress drop value which varies between 5 and 15 bars. If Madariaga's model had been used, higher stress drop values would have been obtained.The low stress drop earthquakes in the eastern Mediterranean region might be interpreted either by the possible presence of low strength material near the source or by the partial stress drop model.     

Propagation of surface waves on a nonhomogenous spherically aeolotropic shell

Summary This paper studies the propagation of Surface Waves on a spherically aeolotropic shell surrounded by vacuum. The elastic constantsc _ij and density of the material of the shell are assumed to be of the form _ij r ^l and _o r ^m respectively, where _{ij o} are constants andl, m are any integers.     

A numerical study of the heat transfer through a fluid layer with recirculating flow between concentric and eccentric spheres

A numerical study has been made of the heat transfer through a fluid layer with recirculating flow. The outer fluid surface was assumed to be spherical, while the inner surface consisted of a sphere concentrically or eccentrically located with respect to the outer spherical surface. The recirculating flow was assumed to be driven by a gas flow creating stress on the fluid's outer surface so that creeping (low Reynolds number) flow developed in its interior. The present study solves the Stokes equation of motion and the convective diffusion equation in bispherical coordinates and presents the streamline and isotherm patterns.Nomenclature a _i inner sphere radius - a _d outer sphere radius - A ₁ defined by equation (5) - A ₂ defined by equation (6) - B ₁ defined by equation (7) - B ₂ defined by equation (8) - c dimensional factor for bispherical coordinates - C constant in equation (4) - d narrowest distance between the two eccentric spheres - E ² operator defined by equation (1) in spherical coordinates and by equation (21) in bispherical coordinates - G modified vorticity, defined in equation (22) - G ^* non-dimensional modified vorticity, defined in equation (28) - h metric coefficient of bispherical coordinate system, defined in equation (18) - k _w thermal conductivity of water - K ₁ defined by equation (9) - K ₂ defined by equation (10) - N _Re Reynolds number=2a _dU/gn - N _Pe,h Peclet number=2a _dU/ - n integer counter - q heat flux - r radius - r ^* non-dimensional radius=r/a _d - S surface area - t time - t ^* non-dimensional time=t/a _d ² - T temperature - T _o temperature at inner sphere surface - T _a temperature at outer sphere surface - T ^* non-dimensional temperature;=(T–T _o)/(T_a–T_o) - u velocity - u _r radial velocity in spherical coordinates - u angular velocity in spherical coordinates - u radial velocity in bispherical coordinates - u angular velocity in bispherical coordinates - U free stream velocity - u _r ^* =u _r/U - u ^* =u /U - u ^* =u /U - u ^* =u /U Greek symbols a ₁ small displacement - vorticity, defined in equation (17) - ^* non-dimensional vorticity, defined in equation (27) - radial bispherical coordinates - _o bispherical coordinate of inner sphere - _a bispherical coordinate of outer sphere - angular coordinate in spherical coordinates - thermal diffusivity - _w thermal diffusivity of water - kinematic viscosity - angular bispherical coordinate - spherical coordinate - streamfunction - non-dimensional streamfunction for spherical coordinates, = /(U a _d ² ) - ^* non-dimensional streamfunction for bispherical coordinates, defined in equation (26)     

Aspherical correction for free vibration of elastically anisotropic solid by means of theXYZ algorithm

By applying the perturbation theory to theXYZ algorithm (a kind of variational method), the difference f in free vibration frequencies between sphere and ellipsoid was approximated as , where i and _i (i = x,y andz) (i=x, y andz) are aspherical coefficients and asphericities of the ellipsoid, respectively. We developed an analytic method to compute the aspherical coefficients719-4 by using theXYZ algorithm. A numerical example was given for an ellipsoidal olivine, and an attempt was made to estimate the asphericities of the specimen by a least-squares method, based on the relationship between frequency shift and asphericity.     

On the propagation of finite Love Waves

Summary The problem of the propagation of finite Love Waves in a heterogeneous elastic half space lying over a homogeneous elastic half space, using the quasilinear stress-strain relation due toS. Ferhst [4] is considered in detail. The variations of the parameter in the layer assumed to be of the form 1= 0e ^z, 0e ^z where is a constant andz is distance measured from the surface into the layer.     

Rayleigh waves in an inhomogeneous medium

Summary Dispersion in Rayleigh waves is discussed for semi-infinite media with = ₁(1 ± cos s z) and = ₁(1 ± cosh s z), being the rigidity of the medium. A few workers tried with the above Fourier type of model but failed to find the dispersive nature. Because they neglected s due to the complexity of the calculation they arrived at a non dispersive frequency equation. This difficulty is removed in this paper and a dispersive frequency equation is obtained which shows both direct and inverse dispersion. The second model leads to non-convergent solution forz but shows many interesting results which are also discussed.     

On the estimation of the heat flow at the earth's crust lower boundary

Summary The paper deals with a comparison of two models used in determining the heat flow q_M at the lower boundary of two-dimensional geothermal models of the Earth's crust. Method I is based on a linear relationship between the component of the surface heat flow, which is not generated by heat sources inside the model, and the heat flow q_M. This method uses a regularization process, in which the variation of the sought function q_M is limited from above. Method II is the most frequently used iterative method, in which the (i + 1)-th approximation of q_M is determined from the surface heat flow, corresponding to the i-th approximation of q_M used as a boundary condition in solving the direct problem. The comparison of both methods has revealed that the solution obtained by method I satisfies the supposed reality better than the solution obtained by method II. Method II is attractive especially for its simplicity. To eliminate the local variations of the estimated heat flow q_M, which are due to automatic transmission from the surface heat flow, a combination of method II with some smoothing procedure could be applied.

---

¶rt;am ¶rt;a m¶rt;a uu mn nma q_M a uau ¶rt; muu ¶rt; . m¶rt; unm u ¶rt; ma nm nma, ma uaa umuau mna mu ¶rt;u, u nm q_M. m m¶rt; nu n uauu, auuau auau u uu q_M . m m¶rt;, mu nm ¶rt; u n a , umau: (i + 1)-a annuau nma q_M um n auma nm nma, ma i- annuauu q_M. au ¶rt;m, m nm q_M, a¶rt; n m¶rt;, mam n¶rt;naa ¶rt;mummu , u, n m m¶rt;. m m¶rt; nuam nmm. ¶rt; mau a auau a¶rt; nma q_M, uu ¶rt;mu amamu na u auau nm nma, ¶rt;a uuam au-u¶rt; m¶rt; auau.

     

Some exact solutions for the flow of fluid through tubes with uniformly porous walls

Summary This paper considers an incompressible fluid flowing through a straight, circular tube whose walls are uniformly porous. The flow is steady and one dimensional. The loss of fluid through the wall is proportional to the mean static pressure in the tube. Several formulations of the wall shear stress are considered; these formulations were motivated by the results from Hamel's radial flow problem, boundary layer flows/and boundary layer suction profiles. For each of these formulations exact solutions for the mean axial velocity and the mean static pressure of the fluid are obtained. Sample results are plotted on graphs. For the constant wall shear stress problem, the theoretical solutions compare favorably with some experimental results.Notations A, B, D, E constant parameters - a, b constant parameters - Ai(z), Bi(z) Airy functions - Ai, Bi derivatives of Airy functions - k constant of proportionality betweenV andp - L length of pores - p,p mean static pressure - p ₀ static pressure outside the tube - p ₀ value ofp atx=0 - Q constant exponent - R inside radius of the tube - T wall shear stress - T ₀ shear parameter - t wall thickness - U free stream velocity - ,u mean axial velocity - u ₀ value ofu atx=0 - V,V mean seepage velocity through the wall - v ₀ mean seepage velocity - x,x axial distance along the tube - z transformed axial distance - z ₀ value ofz atx=0 - mean outflow angle through the wall - cos - density of the fluid - wall shear stress - dynamic viscosity of the fluid - over-bar dimensional terms - no bar nondimensional terms The National Center for Atmospheric Research is sponsored by the National Science Foundation     

Determination of the geopotential scale factor from satellite altimetry

Summary The geopotential scale factor R ₀ =GM/W ₀ has been determined on the basis of satellite altimetry as R _{0=(6 363 672·5±0·3) m} and/or the geopotential value on the geoid W ₀ =(62 636 256·5±3) m ² s ^–2 . It has been stated that R ₀ and/or W ₀ is independent of the tidal distortion of surface W=W ₀ due to the zero frequency tide.

---

¶rt;a nmu amumuu u ama amnmuaa R ₀ =GM/W ₀ =(6 363 672,5±0,3) m u/uu aunmuaa a nmuu¶rt;a W ₀ =(62 636 256,5±3) m² s^–2. m, m R ₀ u/uu W ₀ auum m nm amu a a nuu ¶rt;au nmu W=W ₀ .

     

Transmission of Rayleigh waves over the surface of a heterogeneous medium

Summary The frequency equation of Rayleigh waves propagating over the free surface of an isotropic, perfectly elastic, heterogeneous semi-infinite medium with material properties varying as = ₀ e ^az , = ₀ e ^az , = ₀ e ^az (a>0) has been obtained. Solution of the frequency equation in closed form is obtained in two cases (i) =0, (ii) =, and the Rayleigh wave dispersion curves for phase and group velocities drawn. In both the cases the medium yields single Rayleigh modes which cannot propagate below certain cut-off frequencies. It is found that in case (i), <c<c ₀ and 0.87500 <c _g <c ₀, and in case (ii), 1.03082 <c<c ₁ and 0.90850 <c _g <c ₁, wherec andc _g denote phase nad group velocities respectively, is the constant shear wave velocity of the mediumc ₀ andc ₁ are the corresponding Rayleigh wave velocities of the homogeneous medium of the same Poisson's ratio. The motion of the surface particles is found to be retrograde elliptical as in the homogeneous case, but the ratic of the major and minor axes now becomes frequency dependent and is plotted against frequency. In both the cases (i) and (ii), the ratio starts at a lower value at the cut-off frequency and approaches the corresponding value of the homogeneous medium at high frequencies.     

Die Seismizität der Türkei

Summary Using the fromulae given byGutenberg andRichter, the writer has computed the magnitude and energy of 1804 earthquakes which occurred in Turkey during the period 1850–1960. For drawing the Isenerget, the formula =log₁₀ S has been used in accordance with the definitions given byToperczer andTrapp, whereS=e _i/F·p represents the energy in erg/m² h corresponding to the surface element of 0.5° Lat. x 0.5° Long. Also the relationship between the seismicity and the tectonics of Turkey has been studied by drawing the maps of the epicenters, the focus-depths and the frequences of the earthquakes with various intensities.     

The mean energetic level. theory

Summary One of the important atmospheric levels, the mean energetic level (MEL), which in a sense reflects the energetics of the whole atmosphere, is defined. Its fundamental properties are shown. In order to describe the MEL correctly a new vertical coordinate is introduced and discussed. The new coordinate, , is defined as the ratio of height and temperature. The MEL is shown to be a level with constant value of . Some incorrect conclusions concerning the MEL, derived in the past, have been corrected.List of symbols used c _p specific heat of air at constant pressure - c _v specific heat of air at constant volume - e base of natural logarithms - E total potential energy - f Coriolis parameter - g acceleration of gravity - i specific internal energy - I internal energy - J enthalpy - k unit vector pointing upwards - p pressure - Q diabatic heating rate - R gas constant of the air - t time - T temperature - v horizontal velocity - v ⁽³⁾ three-dimensional velocity - w vertical velocity in thez-system - z height - temperature growth rate (T/z) - Pechala's vertical coordinate (z/T) - generalized vertical velocity in the -system (d/dt) - specific potential energy - potential energy - density of the air - Ruppert function - T(1–)^–1 - ( ) _S quantity at the sea level - ( )* quantity at the MEL     

Experimental investigation of the elastic modulus of a fractal system—A model of fractured rocks

The elastic properties of a physical model representing a damaged rock matrix were studied using a square lattice deformed under tensile stress. The elastic modulusM of such a system varies in agreement with percolation theory as|x–x _c | ^f , wherex is the damage parameter andx _c the threshold value of the damage parameter,f3.6. Atxx _c the scale dependence ofM can be expressed asML ^–f/v , whereL is the size of the sample andv the correlation exponent in percolation theory.The experimental results are of interest in assessing elastic properties in earthquake focal zones and fault zones in general.     

Towards realistic flow modelling. Creation and evaluation of two-dimensional simulated porous media: An image analysis approach

This work is part of an attempt to quantify the relationship between the permeability tensor (K) and the micro-structure of natural porous media. A brief account is first provided of popular theories used to relate the micro-structure toK. Reasons for the lack of predictive power and restricted generality of current models are discussed. An alternative is an empirically based implicit model whereinK is expressed as a consequence of a few pore-types arising from the dynamics of depositional processes. The analytical form of that implicit model arises from evidence of universal association between pore-type and throat size in sandstones and carbonates. An explicit model, relying on the local change of scale technique is then addressed. That explicit model allows, from knowledge of the three-dimensional micro-geometry to calculateK explicitly without having recourse to any constitutive assumptions. The predictive and general character of the explicit model is underlined. The relevance of the change of scale technique is recalled to be contingent on the availability of rock-like three-dimensional synthetic media. A random stationary ergodic process is developed, that allows us to generate three-dimensional synthetic media from a two-dimensional autocorrelation functionr( _x , _y ) and associated probability density function measured on a single binary image. The focus of this work is to ensure the rock-like character of those synthetic media. This is done first through a direct approach:n two-dimensional synthetic media, derived from single set ( ,r( _x , _y )) yieldn permeability tensorsK _i-1,n ⁱ (calculated by the local change of scale) of the same order. This is a necessary condition to ensure thatr( _x , _y ) and carry all structural information relevant toK. The limits of this direct approach, in terms of required Central Process Unit time and Memory is underlined, raising the need for an alternative. This is done by comparing the pore-type content of a sandstone sample andn synthetic media derived fromr( _x , _y ) and measured on that sandstone-sample. Achievement of a good match ensures that the synthetic media comprise the fundamental structural level of all natural sandstones, that is a domainal structure of well-packed clusters of grains bounded by loose-packed pores.     

The effect of the number of charges carried by the argon ion on the darkening of ilford Q2 emulsion as ion detectors in mass-spectrometer

Summary The darkening (S) of Illford Q₂ photographic plates as ion detectors in mass spectrometer has been investigated. The dependence of the darkening (S) on the ion density (n=ions/mm²) i.e.S=S(n)E for constant energy (E)=z U ranging from 4U20 Kv of the impinging⁴⁰A⁺¹-,⁴⁰A⁺²- and⁴⁰A⁺³-ions whenS does not exceed the value 0.15 and the second relationn=n(z U) _S for darkening 0.05S0.15 constructed from the above relationS=S(n) _E has been determined. The darkening was found to increase with increasing ion-density which inturn decreases with the ionenergy. For⁴⁰A⁺¹-,⁴⁰A⁺²-, and⁴⁰A⁺³-ion of equal energy and ion-density the darkening effect was independent of the number of the charges carried by the argon ion.     

(ϖμ/ϖT) P of the lower mantle

We estimate (/T) _P of the lower mantle at seismic frequencies using two distinct approaches by combining ambient laboratory measurements on lower mantle minerals with seismic data. In the first approach, an upper bound is estimated for |(/T) _P | by comparing the shear modulus () profile of PREM with laboratory room-temperature data of extrapolated to high pressures. The second approach employs a seismic tomography constraint ( lnV _S / lnV _P ) _P =1.8–2, which directly relates (/T) _P with (K _S /T) _P . An average (K _S /T) _P can be obtained by comparing the well-established room-temperature compression data for lower mantle minerals with theK _S profile of PREM along several possible adiabats. Both (K _S /T) and (/T) depend on silicon content [or (Mg+Fe)/Sil of the model. For various compositions, the two approaches predict rather distinct (/T) _P vs. (K _S /T) _P curves, which intersect at a composition similar to pyrolite with (/T) _P =–0.02 to –0.035 and (K _S /T) _P =–0.015 to –0.020 GPa/K. The pure perovskite model, on the other hand, yields grossly inconsistent results using the two approaches. We conclude that both vertical and lateral variations in seismic velocities are consistent with variation due to pressure, temperature, and phase transformations of a uniform composition. Additional physical properties of a pyrolite lower mantle are further predicted. Lateral temperature variations are predicted to be about 100–250 K, and the ratio of ( lnp/ lnV _S ) _P around 0.13 and 0.26. All of these parameters increase slightly with depth if the ratio of ( lnV _S / lnV _P ) _P remains constant throughout the lower mantle. These predicted values are in excellent agreement with geodynamic analyses, in which the ratios ( ln / lnV _S ) _P and ( / lnV _S ) _P are free parameters arbitrarily adjusted to fit the tomography and geoid data.     

A note on the stress-dilatancy relation for simulated fault gouge

Theoretical constraints on the stress-dilation relation for a deforming Coulomb material requirev ifC=0 andv sin^-1( _m / _m ) always, wherev is the dilation angle, is the friction angle,C is cohesion, _m is the maximum shear stress, and _m is the mean effective stress. Recent laboratory measurements of friction and dilatancy of simulated fault gouge show that small amplitude shear-load cycling causes compaction and consolidation. Comparison of the data with theory indicates that such load cycling produces: (1) increased coefficient of friction (or friction angle), (2) increased cohesion, and (3) increased dilatancy rate (or dilation angle). Under certain conditions of load cycling without significant plastic shear strain accumulation ( ^p <0.005) we find thatv exceeds both and, in contrast to theory, sin^-1( _m / _m ). This result is interpreted in terms of enhanced cohesion and overconsolidation, which lead to residual stresses within the gouge. An analogy is drawn between these special loading conditions and those extant on natural faults. In particular, our results imply that jostling and minor stress variations associated with microearthquakes may produce strengthening of fault gouge and changes in the fault zone's stress-dilatancy relation. Hence, compaction associated with microseismicity may lead to subsequent dilation of fault gouge, even for faults with large displacement rates and large net offsets (e.g., San Andreas). In regions where such dilation persists over sufficient displacements (on the order of the critical slip distance for seismic faulting) it may tend to inhibit unstable slip.     

Non-adiabatic motion of charged particles traversing a weak magnetic field: Pitch angle scattering

Summary A charged particle moves with velocityv in a constant non-uniform magnetic fieldH, spiralling with Larmor radiusR. IfR is small compared with the scale lengthL of the field, the magnetic moment associated with the Larmor motion of the particle is nearly constant. Consequently , the (pitch) angle betweenv andH, varies as arcsinH ^1/2. Hence in such adiabatic motion is approximately the same at points on the path whereH has the same value. But the magnetic moment and the pitch angle may differ materially at two such points, each in the region whereR/L is small, if between them the particle traverses a region whereR/L is not small. This region of non-adiabatic motion scatters the pitch angles.Such scattering is investigated for regions of weak field (R large), with and without the presence of a neutral line along whichH=0. Either type of region, it is found, can scatter the pitch angles. This gives support to the theory proposed byAkasofu andChapman to explain why auroral arcs and bands are very thin.The scattering here examined is of interest also in connection with magnetic mirror devices for nuclear energy transformation. It may also have applications to phenomena of solar and stellar atmospheres.     

On the accuracy of location of local shocks on the territory of Czechoslovakia and adjacent areas

Summary The calculation procedures for determining epicentre parameters of weak near shocks with foci in Poland are discussed and tested for explosions with known epicentres.

---

m m¶rt; ¶rt; n¶rt;u num a uu m num nmua n auauu, u mu a mumuu u, n muu ¶rt;au uu mau. au mam nam (a. 4) nu nuuu na 71 u m ¶rt;u n¶rt; ¶rt; a auu ¶rt;a [11].