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)     

L'ordre du coefficient de récombinaison dans la coucheE

Résumé On propose une méthode de détermination de l'ordre du coefficient équivalent de récombinaison en n'utilisant que les données ionosphériques. A cette fin on évalue l'apport de la vitesse des changements passagersdN _m E/dt dans la densité électroniqueN _m E. L'évaluation se fait en présence d'une série de valeurs choisies de , on observe laquelle des courbes définis _m (t) dans cette étude coïncidera de plus près avec _m (t) mesurée. Les changements de montrés sur la figure 4 sont obtenus de cette manière. La valeur diurne de varie de (0.5÷1) 10^–7 cm³ sec^–1; peu après le lever et peu avant le coucher du soleil devient >10^–7 cm³ sec^–1. Lorsque cos>0, d'après la règle décrcit et devient 10^–8 cm³ sec^–1. Les explications des changements de obtenus de cette manière sont données en rendant compte des changements par rapport aux ions atomiques et moléculaires, à la disparition rapide des ions positifs pendant le coucher du soleil en présence d'un coefficient de récombinaison plus grand et des changements de température eventuels. La symétrie ou l'asymétrie des fréquences critiquesf ₀ E quand cos sont égaux permettent dans la marche diurne de juger de l'ordre de . Toutes ces évaluations indiquent également des valeurs de environ 10^–7 cm³ sec^–1. En précisant des mesures def ₀ E il est possible de définir non seulement l'ordre, mais aussi la valeur de elle-même.

---

Summary A method for the determination of the order of the equivalent recombination coefficient is suggested, by using ionospheric data only. The increase in the speed of temporary changes indN _m /dt in the electronic densityN _m E is estimated for this purpose. The estimation is done with series of selected for the sake of expediency values of , following at the same time which curves _m (t), determined in the course of work, will coincide most closely with the _m (t) measured. The changes in , shown in figure 4, have been obtained in this way. The diurnal value of is in the range of (0.5 to 1) 10^–7 cm³ sec, being >10^–7 cm³ sec a little after sunrise and a little prior to sunset. At cos<0, by rule decreases and becomes 10^–8 cm³ sec. Explanations of the changes thus obtained in are indicated, taking into consideration the relation of atomic and molecular ions, the rapid disappearance of positive ions at sunset with a higher recombination coefficient and eventual temperature changes. From the symmetry or asymmetry of the critical frequencies off ₀ E at equal cos in the course of the day it can also be judged for the order of . All those estimates show values of in the range of 10^–7 cm³ sec. In the case of precise measurements off ₀ E, not only the determination of the order but also the real value of is possible.

     

Stratospheric ozone variability near the IMF sector boundary

Summary We search for the effects of the interplanetary magnetic field (IMF) sector boundary crossing (SBC) in upper stratospheric ozone. The SBUV data (Nimbus-7) at the 10, 3 and 1 hPa levels are analysed for latitudes 45° N and 55° N for winters of the period December 1979 to December 1982. An effect of the IMF SBC wos only found at the 10 hPa level. These first results concerning the IMF SBC effect in upper stratospheric ozone are rather preliminary.

---

¶rt;m uu nu mau () nam aum n () a mam. SBUV ¶rt;a (u-7) a nm ¶rt;au 10, 3 u 1a aauum ¶rt; um 45° u 55° . . ua nu¶rt;a ¶rt;a 1979 – ¶rt;a 1982. m uu a¶rt; m a 10a. mu n mam n uuu a mam m n¶rt;aumu.

     

Rotatory vibration of a thick spherical shell of isotropic non-homogeneous elastic material

Summary Rotatory vibrations of a thick spherical shell of isotropic non-homogeneous material with rigidity and density given by (i) = ₀ r ^-2 withQ =Q ₀ r ^-2 e ^2mr and (ii) = ₀ r ^m with =Q ₀ r ⁿ have been discussed and the frequency equation is derived with numerical enumeration of frequency in each case.     

Planar charged-particle trajectories in multipole magnetic fields

This paper provides a complete generalization of the classic result that the radius of curvature () of a charged-particle trajectory confined to the equatorial plane of a magnetic dipole is directly proportional to the cube of the particles equatorial distance () from the dipole (i.e. ³). Comparable results are derived for the radii of curvature of all possible planar chargedparticle trajectories in an individual static magnetic multipole of arbitrary order m and degree n. Such trajectories arise wherever there exists a plane (or planes) such that the multipole magnetic field is locally perpendicular to this plane (or planes), everywhere apart from possibly at a set of magnetic neutral lines. Therefore planar trajectories exist in the equatorial plane of an axisymmetric (m = 0), or zonal, magnetic multipole, provided n is odd: the radius of curvature varies directly as ⁿ⁼². This result reduces to the classic one in the case of a zonal magnetic dipole (n = 1). Planar trajectories exist in 2m meridional planes in the case of the general tesseral (0 < m < n) magnetic multipole. These meridional planes are defined by the 2m roots of the equation cos[m()–_n^m)] = 0, where _n^m = (1/m) arctan (h_n^m/g_n^m); g_n^m and h_n^m denote the spherical harmonic coefficients. Equatorial planar trajectories also exist if (n – m) is odd. The polar axis ( = O,) of a tesseral magnetic multipole is a magnetic neutral line if m > I. A further 2_m(n – m) neutral lines exist at the intersections of the 2_m meridional planes with the (n – m) cones defined by the (n – m) roots of the equation P_n^m(cos ) = 0 in the range 0 < 9 < , where P_n^m(cos ) denotes the associated Legendre function. If (n – m) is odd, one of these cones coincides with the equator and the magnetic field is then perpendicular to the equator everywhere apart from the 2m equatorial neutral lines. The radius of curvature of an equatorial trajectory is directly proportional to ⁿ⁼² and inversely proportional to cos[m(–)]. Since this last expression vanishes at the 2m equatorial neutral ines, the radius of curvature becomes infinitely large as the particle approaches any one of these neutral lines. The radius of curvature of a meridional trajectory is directly proportional to rⁿ⁺², where r denotes radial distance from the multiple, and inversely proportional to P_n^m(cos )/sin . Hence the radius of curvature becomes infinitely large if the particle approaches the polar magnetic neutral ine (m > 1) or any one of the 2m(n – m) neutral ines located at the intersections of the 2m meridional planes with the (n – m) cones. Illustrative particle trajectories, derived by stepwise numerical integration of the exact equations of particle motion, are pressented for low-degree (n 3) magnetic multipoles. These computed particle trajectories clearly demonstrate the non-adiabatic scattering of charged particles at magnetic neutral lines. Brief comments are made on the different regions of phase space defined by regular and irregular trajectories.Also Visiting Reader in Physics, University of Sussex, Falmer, Brighton, BN1 9QH, UK     

Some possibilities of predicting radio wave absorption in the lower ionosphere

Summary An attempt is made to show possible ways of predicting radio wave absorption in the midlatitude lower ionosphere using relations between absorption and the intensity of solar ionizing radiation and/or common solar activity indices, and between absorption and f₀F2.

---

aa mu nuau nu a¶rt;u ¶rt;um u u a mu ¶rt; nu u umum uuu uu (uu uu u¶rt;au amumu) u ¶rt; nu u f₀F2.

     

Effects of the imf sector structure in the midlatitude troposphere and stratosphere

Summary The data on geopotential heights and temperatures at 7 pressure levels between 1000-10 hPa above Berlin(52.5 °N, 13.4 °E) are analysed for the winters of 1963–1973. No demonstrable effect of the interplanetary magnetic field sector boundary crossing (IMF SBC) is found in the lower and middle stratosphere, but there is a demonstrable effect in the middle troposphere at the 500 hPa level. This effect is less important than the IMF SBC effect in the tropospheric vorticity area index and seems to be of a different type.

---

auum ¶rt;a nnmua m u mnam a 7 nm ¶rt;au ¶rt; 1000-10 a a¶rt; u(52,5 °.., 13,4 °.¶rt;.) ¶rt; u 1963–1973. ua ¶rt;aam m nu mau nam aum n( ) ¶rt;a amu u u ¶rt; mam, ma m a¶rt; ¶rt; mn a 500 a. mm m a, m u¶rt; na¶rt;u aumu am, u am m ¶rt; muna.

     

Characteristics of the mid-latitude ionosphere suitable for modelling ionospheric filtration of ULF waves

Summary An approximate method of one-dimensional modelling of the plasma of the Earth's ionosphere is demonstrated for purposes of studying the ionospheric filtration of ULF waves (micropulsations). Apart from the basic local parameters, characterizing the plasma, also derived local characteristics have been defined, i.e. the mass of the so-called effective ion and its effective collision frequency . Drawing on existing empirical models of the mid-latitude ionosphere, vertical profiles (50 km h 1000 km) were determined of the characteristics N_e N_i, v_e, and for the daytime and nighttime mid-latitude ionosphere under low and enhanced solar activity, which can be used to study the ionospheric ULF filter.

---

aa nu uum m¶rt; ¶rt;a ¶rt;uau na ¶rt; u ¶rt; nmm uu u umauu (unau). ¶rt; u mu naamau na n¶rt; mu m aamumuu — aa m. a. mu ua, , u mua amma mu, . a mu nuuu ¶rt; ¶rt;um u u n¶rt; mua nuu (50 h 1000 ) aamumu N_e N_i, v_e, u ¶rt; u u ¶rt;u um nu u u amumu. ¶rt;u nam nu uuu u uma.

     

The ionosphere of higher geomagnetic latitudes (L=3 andL=5) as a ulf wave filter

Summary Eight models of the ionosphere in the MHD approximation over the whole range of required heights(50–3000 km) are synthesized on the basis of other theoretical studies and a number of empirical methods for the purposes of modelling the ionospheric filtration of the natural micropulsation (Pc1) signal. The models represent a vertically inhomogeneous and dissipative daytime and nighttime ionosphere of higher latitudes in two regions (L=3 and L=5) of the Northern Hemisphere in summer under low (F _10.7 =70) and high (F _10.7 =200) solar activity. The higher ionosphere (h>200 km) is approximated by a quasineutral oxygen-hydrogen plasma(O ⁺,H ⁺) taking into account the wave-dissipating effects of the neutral component(O, He). The procedure of computing the physical parameters, locally characterizing the medium of the ionosphere, is demonstrated.

---

a ¶rt;a mmuu am u a nuuu m¶rt; ¶rt;a um u ¶rt; u ( nuuuu) n a m (50–3000 ), m ¶rt; ¶rt;uau u umauu mm unau (1) uaa. ¶rt;u n¶rt;mam m-¶rt; ¶rt;, na ¶rt; u m u u um ¶rt; am (L=3 u L=5) n, uu nu u, (F _10,7 =70), (F _10,7 =200) amumu. a ua (h>200 ) n¶rt;maa nuuuu auma u¶rt;-¶rt; ¶rt;(O ⁺,H ⁺) na, uau nuuam nu u-a ma maQu(O, He). aa m¶rt; uu uuu naam, m aamuu u ¶rt;.

     

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.     

The effect of the velocity of the centre of a cyclone on the generation of microseisms

Summary The influence of the velocity of the movement of the centre of the cycloneV _c.c. on the rate of amplitudes' change A/t and periods' change T/t of storm microseisms is investigated. The dependence A/t=k V _c.c. and T/t=k ¹ V _c.c. is obtained. Unmovable depression (V _c.c. =0) does not stipulate the change of A/t and T/t.

---

u V _c.c. A/t T/t . A/t=k V _c.c. T/t=¹ V _c.c. . (V _c.c. =0) A/t T/t.

---

---

Presented as a scientific communication to the IASPEI Assembly in Madrid, Sept. 1969.     

Anisotropic magnetic susceptibility of rocks under mechanical stresses

Summary The magnetic susceptibility of a rock under a uniaxial compression () decreases along the axis of compression and increases along the direction perpendicular to the axis, with an increase of . Thus, the magnetic susceptibility of a compressed rock becomes anisotropic.The decrease of longitudinal susceptibility,K (), and the increase of transverse susceptibility,K (), are theoretically derived from a model of rock which assumes the uniaxial anisotropy and the isotropic magnetostriction of magnetic minerals in rocks and a random orientation of the minerals. Results show thatK () decreases toward zero whereasK () increases and approaches a finite asymptotic value with an increase of , and –(/)K () is twice as large as /K () for small values of . These results are in good agreement with experimental data.

---

Zusammenfassung Die magnetische Suszeptibilität eines Steines unter zunehmender uniachsigen Druckspannung () nimmt ab längs der Achse der Druckspannung und nimmt zu längs der Richtung senkrecht der Achse. Somit wird die magnetische Suszeptibilität des gedrückten Steines anisotrop.Die Abnahme der longitudinalen Suszeptibilität,K (), und die Zunahme der transversalen Suszeptibilität,K (), werden theoretisch von einem Modell eines Steines hergeleitet, das die uniachsige Anisotropie, die isotrope Magnetostriktion, und eine nichtbevorzugte Orientierung der magnetischen Minerals im Stein annimmt. Die Ergebnisse zeigen, dass mit einer Zunahme des ,K () gegen Null abnimmt, währendK () zunimmt und sich einem begrenzten asymtotitschen Wert nähert und, dass für kleine Werte des , –(/)K () zweimal so gross wie /K () ist. Diese Ergebnisse stimmen gut mit den Versuchangaben überein.

     

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.     

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

     

La fonction d'appui dans les formules de la géodésie mathématique

na ma¶rt;uu n ¶rt;a amamu¶rt;uu ¶rt; una, unu¶rt;a au u m unu¶rt;a — u nua na u n u, m. . nmu, ¶rt; amamu uu.     

Simple model for long-range forecasts of AT500 hPa level fluctuations

Summary A forecast model to predict the fluctuations of level AT 500 hPa in a selected grid of points is derived. The solution of the compensation equation is sought in the form of a trigonometric polynomial in three idependent variables. It constitutes the basis of a numerical solution of the prediction problem with the use of a high-speed computer. A three-month forecast of the altitude fluctuation of level 500 hPa is evaluated by means of the daily values of the correlation coefficient. The results are satisfactory and the general evaluation shows the model to be prospective.

---

¶rt;um nmua ¶rt; ¶rt; na au m AT 500a a uum m m. u au nauu um u¶rt; mu mu nua m m auu n. m u u nmu a¶rt;au a . u¶rt;um a m na au m 500a nu nu m au uuma uu. mam m anumu, u nu ¶rt; n¶rt;mam nnmu.

---

List of symbols used * coupling coefficient between two conjugate atmospheric elements - * ageostrophicity coefficient of the atmospheric system - Coriolis' parameter (=2 sin) - , * geographic latitude, geographic longitude - *, geopotential reference and pressure levels - , * compensation and coupling frequencies - integration field over the entire atmospheric system - A** constant (A*=2 ²(*+*)) - A _r ,B _r ,C _r ,D _r constants related to subsystem -r- - A _s ,B _s ,C _s ,D _s constants related to subsystem -s- - B** constant (B**= ⁴(*²+2**)) - C _r constant (i=1, 2, 3, 4) - E _k ,E _v ,E _p energies of the atmospheric system: kinetic, internal and potential - K transformation constant - m total number of generalized frequencies - R(T) frequency characteristics of the numerical band filter - r0(t),r1(t) daily values of the correlation coefficient - Q heat - x, y coordinates in the reference plane - t time - _p ² Laplace's differential operator in thep-system     

Seismic response of weakly laterally inhomogeneous structure

Summary The seismic response (the transfer properties) of a model consisting of two layers above a homogeneous halfspace, with curved interfaces and continuously laterally varying velocity within the layers is investigated assuming vertical plane SH-wave incidence from below. The frequency response and the synthetic seismograms are computed by the ray and matrix methods. The dependence of the seismic response on the prevailing frequencies of the incident wave, the role of the individual layers in forming the seismic response, as well as the effect of causal absorption are studied.

---

aamuam ua au (n¶rt;am ma) ¶rt;u, m u n¶rt;mua ¶rt;¶rt; nnmam ¶rt; , uuuauau a¶rt;a u n u anau m mu , n¶rt;nuu mua na¶rt;a n SH-. amm auu u ummuu a auma u amu m¶rt;au. ¶rt;am auum u auu m na¶rt;a amm na¶rt;a , m¶rt; nu uauu u auu u ma m nu.

     

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     

Q measurements for PhaseX overtones

Linear stacking procedures are used to retrieve the attenuation of 91 modes belonging to the 3rd, 4th and 5th Rayleigh overtones branches in the 80–160 s period range, and contributing to the so-called PhaseX wave group. Our data show in general slightly less attenuation than expected from available models. Data space inversion shows that, when combined with previously measured fundamental modeQ's, this new dataset improves resolution significantly in the 1000–2000 km depth range. Based on this remark, we carry out a number of parameter space inversions. Our results suggest a narrow (80–200 km) zone of high attenuation (Q =75–90), low attenuation in the intermediate mantle (670–1500 km); (Q 350), and lower values in the deeper mantle (Q 200).     

Palaeointensity of the geomagnetic field during upper cainozoic derived from palaeo-slags and porcellanites in north Bohemia

Summary Porcellanites and palaeo-slags from North Bohemia are natural materials which can be used to derive the palaeomagnetic directions and palaeointensity of the geomagnetic field active at the time of caustic alteration. The origin of these rocks, called erdbrands, was due to the caustic alteration of predominantly pelitic sediments as a result of underground fires conditioned by spontaneous ignition of coal seams. The caustic alteration occurred during the Upper Pliocene to the Quaternary. Three procedure based on the methods by Thellier and Nagata are presented in the paper. The newly developed apparatus MAVACS (Magnetic Vacuum Control System) was used for the thermal demagnetization of samples. A procedure based on multi-component analysis was also proposed and tested. Besides some methodic results, it was found that the geomagnetic field intensity varied during the respective period within the limits of 48%±4% to 154%±32% of the present geomagnetic field intensity.

---

aum u na au a mumuu uu n¶rt;mam nu¶rt; amua, m n¶rt;¶rt;um ¶rt; ¶rt;u naaum anau u naumumuaum n, ¶rt;m amu uu. mu n¶rt;, aa ¶rt;a¶rt;, uu n¶rt; uu amu uu num num a¶rt; n¶rt; ¶rt;mu aau . amu uu u m nu¶rt; m nua ¶rt; mmu nu¶rt;a. am n¶rt; mu umnmau nua, n¶rt;¶rt;u ¶rt; ¶rt;u naumumu, nuau a m¶rt; u aama. a aamaa annaama MAVACS (Magnetic Vacuum Control System) a unaa ¶rt; mu aauuau ¶rt; amu aa. n¶rt; u n n¶rt;¶rt;, a a munm aau amuauu. nu m¶rt;uu au, ma ma, m umumaum n u¶rt; nu¶rt; a n¶rt;a 48%±4%-154%±32% au umumu aum n.