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.     

Cyclic variations of the initial height of the sporadicE-layer at middle geographical latitudes

Summary The variations of the initialh E _s height are investigated in the solar cycle 1957–1968, deriving the regressive dependency:h E _s =121.4–6·10^–2 R referring to the median monthly values at a solar zenith angle =75°. The similar variations ofh E _s (R) during the day and night are interpreted as a domination of the sporadic layer formation from a redistribution of the day-time ionization and secondary participation of nightly ionizing sources. The analogous cyclich E _s andh E variations confirm this conclusion while the seasonal variations in the state of the sporadic layer show outlined dynamical effects. The comparatively not big cyclic variation in the spatial state of theE-region are considered to confirm the predominating ionizing action of the ultraviolet range (933–1038 Å) in the lower part of theE-region, while the soft X-radiation influences mainly the near maximum part of this region.     

Quasi-static thermal deformations of a sphere by radiation from a point source

Summary In this paper the quasi-static temperature and stress distributions set up in an elastic sphere by radiation from a point source at a finite distance from the centre of the sphere and out-side it, have been discussed. The temperature boundary condition has been taken in the general form involving an arbitrary function of time. The final solutions have been obtained in terms of series involving Legendre polynomials. Numerical calculations have been done on IBM 1620 Computer and a desk calculator. The results have been represented in graphs.Notation the del operator - u the displacement vector - T the excess of temperature over that at state of zero stress and strain - , Lamé's constants - /2(+) Poisson's ratio - coefficient of linear expansion - 2(1+) - a radius of the sphere - d distance of the point source from the centre of the sphere - d _o a/d - K coefficient of thermal conductivity - h heat transfer coefficient of the surface     

Relations entre le creusement et le comblement des perturbations et leur déplacement

Résumé On commence par définir le creusement et le comblement d'une fonctionp(, t) du tempst et des points (, ) d'une surface régulière fermée en se donnant, sur cette surface, un vecteur vitesse d'advection ou de transfert tangent à . Le creusement (ou le comblement) est la variation dep sur les particules fictives se déplaçant constamment et partout à la vitesse , A chaque vecteur et pour un mêmep(, ,t) correspond naturellement une fonction creusementC (, ,t) admissible a priori; mais une condition analytique très générale (l'intégrale du creusement sur toute la surface fermée du champ est nulle à chaque instant), à laquelle satisfont les fonctions de perturbation sur les surfaces géopotentielles, permet de restreindre beaucoup la généralité des vecteurs d'advection admissibles a priori et conduit à des vecteurs de la forme: , oùT est un scalaire régulier, () une fonction régulière de la latitude , le vecteur unitaire des verticales ascendantes etR/2 une constante. Ces vecteurs sont donc une généralisation naturelle des vitesses géostrophiques attachées à tout scalaire régulier. Dans le cas oùp(, ,t) est la perturbation de la pression sur la surface du géoïde, le vecteur d'advection par rapport auquel on doit définir le creusement est précisément une vitesse géostrophique: on a alors ()=sin etT un certain champ bien défini de température moyenne.On déduit ensuite une formule générale de géométrie et de cinématique différentielles reliant la vitesse de déplacement d'un centre ou d'un col d'un champp(, ,t) à son champ de creusementC (, ,t) et au vecteur d'advection correspondant. Cette formule peut être transformée et prend la forme d'une relation générale entre le creusement (ou le comblement) d'un centre ou d'un col et la vitesse de son déplacement, sans que le vecteur d'advection intervienne explicitement. On analyse alors les conséquences de ces formules dans les cas suivants: 1^o) perturbations circulaires dans le voisinage du centre; 2^o) perturbations ayant, dans le voisinage du centre, un axe de symétrie normal ou tangent à la vitesse du centre; 3^o) évolution normale des cyclones tropicaux.Finalement, on examine les relations qui existent entre le creusement ou le comblement d'un champ, le vecteur d'advection et la configuration des iso-lignes du champ dans le voisinage d'un centre.Ces considérations permettent d'expliquer plusieurs propriétés bien connues du comportement des perturbations dans différentes régions.

---

Summary The deepening and filling (development) of a functionp(, ,t) of the timet and the points (, ) of a regular closed surface is first of all defined, in respect to a given advection or transfer velocity field tangent to , as the variation ofp on any fictitious particle moving constantly and everywhere with the velocity . For a givenp(, ,t) and to any there corresponds a well defined development fieldC (, ,t). All theseC fields are a priori admissible, but a very general analytical condition of the perturbation fields in synoptic meteorology (the integral of the development fieldC (, ,t) on any geopotential surface vanishes at any moment), leads to an important restriction to advection vectors of the form: , whereT is any regular scalar, () any regular function of latitude, the unit vector of the ascending verticals andR/2 a constant. These vectors are a natural generalisation of the geostrophic velocities attached to any regular scalar. Whenp(, ,t) is the pressure perturbation at sea level, its development must be defined in respect to a geostrophic advection vector belonging to the above defined class of vectors with ()=sin andT a well defined mean temperature field.A general formula of the differential geometry and kinematics ofp(, ,t) is then derived, giving the velocity of any centre and col of ap(, ,t) as a function of the advection vector and the corresponding development fieldC (, ,t). This formula can be transformed and takes the form of a general relation between the deepening (and filling) of a centre (or a col) of ap(, ,t) and its displament velocity, the advection vector appearing no more explicitly. A detailed analysis of the consequences of these formulae is then given for the following cases: 1^o) circular perturbations in the vicinity of a centre; 2^o) perturbations having, in the vicinity of a centre, an axis of symmetry normal or tangent to the velocity of the centre; 3^o) normal evolution of the tropical cyclones.Finally, the relations between the developmentC (, ,t) of a fieldp(, ,t), the advection velocity vector and the configuration of the iso-lines in the vicinity of a centre are analysed.These theoretical results give a rational explanation of several well known properties of the behaviour of the perturbations in different geographical regions.

---

---

Communication à la 2ème Assemblée de la «Società Italiana di Geofisica e Meteorologia» (Gênes, 23–25 Avril 1954).     

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.     

Generalized thermoelastic longitudinal waves in unbounded medium

Summary In this paper the generalized thermoelastic longitudinal waves and the temperature field set up due to coupling of the displacement and the temperature fields, with heat wave travelling with certain finite velocity, in an unbounded medium are studied. The thermoelastic displacement potential and the temperature field at any point are obtained in terms of the surface integrals involving the potential, the temperature and their normal derivatives.Notation x _i the cartesian coordinate system,i=1,2,3 - n (u _i) the displacement vector - (/x _i) the del operator - _t /t the derivative with respect to time - T _o the temperature corresponding to the natural stat of zero stress and strain - T Absolute temperature - c _e the specific heat - , Lamé's constants - ₀ the density - coefficient of linear thermal expansion - K thermal conductivity coefficient - _kk u - ₀ the relaxation time     

Détermination des paramètres d'Ångström par des observations actionmétriques courantes

Summary By comparing the equalities expressing the solar flux within a wide spectral region on the assumption of an extintion within the atmospheric aerosol given by ^– and that of an extintion represented by ₁^–1, a relationship is established by means of which: 1) The error is discussed which is done when the air opacity is expressed by ₁, in the case when 1; 2) A simple and quick method is worked out for determining the parameters and from actinometric observations, carried out within the spectral regions <525 m and 525 m<<625 m.     

(ϖμ/ϖ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.     

Turbidometric estimations in Mexico city using the Volz sun photometer

Summary Measuring, with the aid of two filters, the instantaneous intensity of the solar radiation in two wave lengths ( _B = 0.44 , _R = 0.64 ) by means of a sun photometer designed byVolz, we carried out determinations of the decadic turbidity coefficientB (=0.5 ) and the wave length exponent of the haze extinction for Mexico City. Observations were made for almost two and a half years (1960 to 1962 period). A seasonal size distribution in both parameters was found. Although the data thus obtained are provenient of a contaminated atmosphere, comparison of our data is made with those found for higher latitudes ofÅngström, Schüepp andVolz. The height of the homogeneous haze layerH _D was calculated showing pronounced variations for a given wind direction. The maximum and minimum values ofB enable us to get, by the first approximation, the aerosol size distribution ofJunge for our latitudes. However, for exceptional very clear days having maximum actinometric intensity of the solar radiation the sensitivity of the microamperimeter in theVolz sun photometer fails.     

Reservoir lifetime and heat recovery factor in geothermal aquifers used for urban heating

Simple models are discussed to evaluate reservoir lifetime and heat recovery factor in geothermal aquifers used for urban heating. By comparing various single well and doublet production schemes, it is shown that reinjection of heat depleted water greatly enhances heat recovery and reservoir lifetime, and can be optimized for maximum heat production. It is concluded that geothermal aquifer production should be unitized, as is already done in oil and gas reservoirs.Nomenclature a distance between doublets in multi-doublet patterns, meters - A area of aquifer at base temperature, m² drainage area of individual doublets in multidoublet patterns, m² - D distance between doublet wells, meters - h aquifer thickness, meters - H water head, meters - Q production rate, m³/sec. - r _e aquifer radius, meters - r _w well radius, meters - R _g heat recovery factor, fraction - S water level drawdown, meters - t producing time, sec. - T aquifer transmissivity, m²/sec. - v stream-channel water velocity, m/sec. - actual temperature change, °C - theoretical temperature change, °C - water temperature, °C - heat conductivity, W/m/°C - _r rock heat conductivity, W/m/°C - _aC_a aquifer heat capacity, J/m³/°C - _aC_r rock heat capacity, J/m³/°C - _WC_W water heat capacity, J/m³/°C - aquifer porosity, fraction     

On estimating frequency-magnitude relations from heterogeneous catalogs

The paper studies the effect of magnitude errors on heterogeneous catalogs, by applying the apparent magnitude theory (seeTinti andMulargia, 1985a), which proves to be the most natural and rigorous approach to the problem. Heterogeneities in seismic catalogs are due to a number of various sources and affect both instrumental as well as noninstrumental earthquake compilations.The most frequent basis of heterogeneity is certainly that the recent instrumental records are to be combined with the historic and prehistoric event listings to secure a time coverage, considerably longer than the recurrence time of the major earthquakes. Therefore the case which attracts the greatest attention in the present analysis is that of a catalog consisting of a subset of higher quality data, generallyS ₁, spanning the interval T ₁ (the instrumental catalog), and of a second subset of more uncertain magnitude determination, generallyS ₂, covering a vastly longer interval T ₂ (the historic and/or the geologic catalog). The magnitude threshold of the subcatalogS ₁ is supposedly smaller than that ofS ₂, which, as we will see, is one of the major causes of discrepancy between the apparent magnitude and the true magnitude distributions. We will further suppose that true magnitude occurrences conform to theGutenberg-Richter (GR) law, because the assumption simplified the analysis without reducing the relevancy of our findings.The main results are: 1) the apparent occurrence rate exceeds the true occurrence rate from a certain magnitude onward, saym _GR; 2) the apparent occurrence rate shows two distinct GR regimes separated by an intermediate transition region. The offset between the two regimes is the essential outcome ofS ₁ being heterogeneous with respect toS ₂. The most important consequences of this study are that: 1) it provides a basis to infer the parameters of the true magnitude distribution, by correcting the bias deriving from heterogeneous magnitude errors; 2) it demonstrates that the double GR decay, that several authors have taken as the incontestable proof of the failure of the GR law and of the experimental evidence of the characteristic earthquake theory, is instead perfectly consistent with a GR-type seismicity.     

Strong Potential Magnetic Field Influence on Slightly Differentially Rotating Spherical Shell

An electrically conducting viscous fluid-filled spherical shell is permeated by an axisymmetric strong potential magnetic field with large Elssaser number ² 1. We describe analytically the steady flow driven by a slightly faster rotation of the conducting inner boundary of the shell. The main flow is controlled by Ekman-Hartmann boundary layers with a small thickness /, where ² is the Ekman number. Asymptotics based on small ^–1 1 reveal the nature of a free shear layer O((/)^1/2) and a super-rotation that allows a part of the fluid to rotate faster than the inner sphere. The free shear is following an imposed field line that is tangent to the inner or outer sphere. Meridional flux is concentrated in the shear and boundary layers. Fluid tends to rotate with the inner sphere and to expel azimuthal magnetic field from an -region restricted by the free shear in the spherical shell. For an imposed axial uniform magnetic field, this -region is outside the cylinder tangent to the inner sphere and rotates with the outer sphere. Weak differential rotation O(/) is inside the cylinder, while almost all difference in rotation rates between spheres is accommodated in the thin O((/)^1/2) free shear. For an imposed dipole magnet, the region has a shape of a lobe touching the outer equator. Inside a super-rotation exists; this is the common case for such when the source of the imposed field is inside.     

A shot noise model for the generation of soil moisture data

An attempt is made to estimate the expected contribution of rainfall to soil moisture during the irrigation season. Effective rainfall and evapotranspiration are the parameters considered in the water balance carried out in the root zone. Rainfall occurrence is simulated by a Poisson process whereas evapotranspiration is described by a simple deterministic function of potential evapotranspiration and soil moisture in the root zone. Using the theory of shot noise models a closed form solution is derived from the expected soil moisture in the root zone at the end of the time interval (0,t]. For illustration purposes the proposed model is applied to a series of data from Mikra meteorological station in Greece.List of symbols x change in water storage in the root zone during the time interval t - X infiltrated rainfall of thei th storm event - ET actual evapotranspiration during thej th day - Poisson rate - number of storm events in (0,t] - t _r duration of rainfall - t _b interarrival time - h _i rainfall depth of thei th storm event - i _m mean rainfall intensity - i(t) instantaneous rainfall intensity - x(0),x(t) available soil moisture in the root zone at time 0 andt, respectively - PET potential evapotranspiration rate - x _F available soil moisture in the root zone at field capacity - soil moisture depletion rate (=PET/x _F ) - w impulse shape of filtered Poisson processes - E[·] mean value - S _i time of thei th rainfall event - N(t) time of storm events in (0,t] - estimated standard deviation The following symbols were used in this paper     

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 relationship between phonon conductivity and seismic parameter Φ for silicate minerals

Summary The relationship between the phonon conductivity at room temperature (K _N ) and the seismic parameter () for silicate minerals is suggested. The considerations are based on the Debye model of thermal energy transport phenomena in solids and on the seismic equation of state for silicates and oxides given byAnderson (1967). The semiempirical relationship is the formK _N = 0.43^0.82 where is in km²/s² andK _N in mcal/cm s K, and the empirical relationship isK _N =(0.528±0.006) –(8.18±2.11). The laboratory data on thermal and elastic properties for several silicates were taken fromHorai andSimmons (1970).     

Nonuniform friction as a physical basis for earthquake mechanics

A review of simple models and observations suggests that the main first-order features of active faulting-mechanical instability, the frequency-magnitude relations, seismic and aseismie slip, seismic radiation, incoherency and rupture stoppage — may be explained by a single characteristic of crustal faults: the spatial variation of the effective frictional stress, which resists slippage on faults. Faultoffset data suggest that rupture propagation ceases in regions of high resistance which act, as barriers. In these regions slippage is associated with negative stress drop. The spacing and the amplitudeA() of the barriers, as inferred from the frequency-magnitude and moment relation for earthquakes, obeys a simple statistical relationA()^p. On the scale of particle motion, this variability of frictional stress provides a mechanical instability which may be associated with the concept of dynamic friction. Invariably, the rapid particle motion in the model is always preceded by accelerated creep. The particle acceleration is highly irregular, giving rise to an almost random acceleration record on the fault. The particle displacement is relatively smooth, giving rise to simple displacement time function in the far field. Rupture propagation time is approximately proportional to the gradient of frictional stress along the fault. Consequently sharp changes of this stress may cause multiple events and other long period irregularities in the fault motion.The power density spectrum associated with the frictional stress implies that stress may be related to a Poisson distribution of lengths. The autocorrelation of such type of distribution yields a correlation lengthk _L ^–1 , similar perhaps toHaskell's (1964) andAki's (1967) correlation lengths inferred from spectral analysis of seismic waves. The partial incoherency of faulting implies that preseismic deformation may be significantly incoherent, consequently the prediction of small moderate earthquakes may be subject to inherent uncertainties. We conclude that frictional stress heterogeneities may be necessary and sufficient to explain active faulting associated with small and moderate earthquakes.     

Efficient Methods to Simulate Planar Free Surface in the 3D 4th-Order Staggered-Grid Finite-Difference Schemes

We numerically tested accuracy of two formulations of Levander's (1988) stress-imaging technique for simulating a planar free surface in the 4^th-order staggered-grid finite-difference schemes. We have found that both formulations (one with normal stress-tensor components at the surface, the other with shear stress-tensor components at the surface) require at least 10 grid spacings per minimum wavelength ( _min÷h = 10) if Rayleigh waves are to be propagated without significant grid dispersion in the range of epicentral distances up to 15 _dom ^S.Because interior 4^th-order staggered-grid schemes usually do not require more than 6 grid spacings per minimum wavelength, in the considered range of epicentral distances, it was desirable to find alternative techniques to simulate a planar free surface, which would not require denser spatial sampling than _min÷h = 6. Therefore, we have developed and tested new techniques: 1. Combination of the stress imaging (with the shear stress-tensor components at the surface) with Rodrigues' (1993) vertically refined grid near the free surface. 2. Application of the adjusted finite-difference approximations to the z-derivatives at the grid points at and below the surface that uses no virtual values above the surface and no stress imaging. The normal stress-tensor components are at the surface in one formulation, while the shear stress-tensor components are at the surface in the other formulation.The three developed formulations give for the spatial sampling _min÷h = 6 results very close to those obtained by the discrete-wavenumber method. Because, however, the technique with the vertically refined grid near the free surface requires 3 times smaller time step (due to the refined grid), the technique with adjusted finite-difference approximations is the most accurate and efficient technique from the examined formulations in the homogeneous halfspace.     

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 ₀ .

     

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     

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.