Numerical determination of families of three-dimensional double-symmetric periodic orbits in the restricted three-body problem,II

New families of three-dimensional double-symmetric periodic orbits are determined numerically in the Sun-Jupiter case of the restricted three-body problem. These families bifurcate from the vertical-critical orbits ( _v = – 1, b _v – 0) of the basic plane familiesi, g ₁, g₂, c andI. Further, the predictor-corrector procedure employed to reveal these families has been described and interesting numerical results have been pointed out. Also, computer plots of the orbits of these families have been shown in conical projections.     

Families of three-dimensional axisymmetric periodic orbits in the restricted three-body problem: Sun-Jupiter case

Families of three-dimensional axisymmetric periodic orbits are determined numerically in the Sun-Jupiter case of the restricted three-body problem. These families bifurcate from the vertical-critical orbits (_v = 1,b _v = 0) of the basic plane familiesi andI. Further the predictor-corrector procedure employed to reveal these families has been described and interesting numerical results have been pointed out. Also, computer plots of the orbits of these families have been shown in conical projections.     

Families of three-dimensional plane-symmetric periodic orbits in the restricted three-body problem: Sun-Jupiter case

Two new families of three-dimensional simple-symmetric periodic orbits are determined numerically in the Sun-Jupiter case of the restricted three-body problem. These families emanate from the vertical-critical orbits (_v = 1,c _v = 0)of the familiesi andl of plane symmetric simpleperiodic orbits direct around the Sun and the Sun-Jupiter respectively. Further, the numerical technique employed in the determination of these families has been described and interesting results have been pointed out. Also, computer plots of the orbits of these families have been shown in conical projections.     

Parameter values for stable low-inclination periodic motion in the restricted three-body problem with oblateness

The known intervals of possible stability, on the mgr-axis, of basicfamilies of 3D periodic orbits in the restricted three-body problem areextended into -A₁ regions for oblate larger primary, A ₁ beingthe oblateness coefficient. Eight regions, corresponding to the basicstable bifurcation orbits l1v, l1v, l2v, l3v, m1v, m1v,m2v, i1v are determined and related branching 3D periodic orbits arecomputed systematically and tested for stability. The regions for l1v,m1v and m2v survive the test emerging as the regions allowing thesimplest types of stable low inclination 3D motion. For l1v, l2v,l3v, m1v and m2v oblateness seems to have a stabilising effect,while stability of i1v survives only for a very small range of A ₁values.     

Multiple expansion of the tidal potential

The Earth tidal deformation causes an additional gravitational potential. Its effect on the Moon orbital motion has been studied by several authors.In this contribution, we develop this additional potential without specifying the inertial frame chosen.For this purpose, we use the properties of the representation of rotation groups in 3 dimensions space. We finally obtain the interaction potential between the distorted Earth and the Moon which is a necessary preliminary to the study of the evolution of the Earth-Moon system.Nomenclature T.R.O Tide raising object - (, , ) Spherical coordinates of the T.R.O. - (J, _E ) Earth spin axis orientation. _E is the longitude of the ascending node of Earth's equator on thexy-plane - (a ,I ,e , , ,M ) Elliptics elements of the T.R.O     

Effect of perturbations in Coriolis and centrifugal forces on the stability of libration points in the restricted problem

The location and the stability of the libration points in the restricted problem have been studied when small perturbation and are given to the Coriolis and the centrifugal forces respectively. It is seen that the pointsL ₄ andL ₅ form nearly equilateral triangles with the primaries and the pointsL ₁,L ₂,L ₃ remain collinear. It is further observed that for the pointsL ₄ andL ₅, the range of stability increases or decreases depending upon whether the point (, ) lies in one or the other of the two parts in which the (, ) plane is divided by the line 36-19=0 and the stability of the collinear points is not influenced by the perturbations and they remain unstable.     

The motion of a satellite in an axi-symmetric gravitational field

The equations for the variation of the osculating elements of a satellite moving in an axi-symmetric gravitational field are integrated to yield the complete first-order perturbations for the elements of the orbit. The expressions obtained include the effects produced by the second to eighth spherical harmonics. The orbital elements are presented in the most general form of summations by means of Hansen coefficients. Due to their general forms it is a simple matter to estimate the perturbations of any higher harmonic by simply increasing the index of summation. Finally, this paper gives the respective general expressions for the secular perturbations of the orbital elements. The formulae presented should be useful for the reductions of Earth-satellite observations and geopotential studies based on them.List of Symbols semi-major axis - C _jk ⁿ (, ) cosine functions of and - e eccentricity of the orbit - f acceleration vector of perturbing force - f sin²t - i inclination of the orbit - J _n coefficients in the potential expansion - M mean anomaly - n mean motion - p semi-latus rectum of the orbit - R, S, andW components of the perturbing acceleration - r radius-vector of satellite - r magnitude ofr - S _jk ⁿ (, ) sine functions of and - T time of perigee passage - u argument of latitude - U gravitational potential - true anomaly - V perturbing potential - G(M₊+m) (gravitational constant times the sum of the masses of Earth and satellite) - _n,k coefficients ofR component of disturbing acceleration (funtions off) - _n,k coefficients ofS andW components of disturbing acceleration (functions off) - mean anomaly at timet=0 - X ₀ ^n,m zero-order Hansen coefficients - argument of perigee - right ascension of the ascending node     

The photometric solution of cc com by use of the

The contact binary system CC Com (=12^h09^m33^s.8, =+22°4339, (1950);V _max=11.31, (B-V)_max=1.24) is a W UMa-type system with the shortest known period. The photometric solution of CC Com is presented using the Wilson and Devinney method. The results show that the CC Com belongs to the late-type eclipsing binary with the spectral type K5V and K6V, low temperatureT ₁=4300 K,T ₂=4265 K, the mass ratioq=0.5873±0.0021, and the inclinationi=87°.719±1°.44. The best regions of the gravity darkening exponents , the bolomotric albedov, and the limb-darkening coefficients are tested. It is found that 0.1250.065, 0.1v0.5, =0.5 are better regions for CC Com. The third body ofl ₃ is not found to be significant. The results are combined with the spectroscopic results of Rucinski to provide an estimate of the absolute parameters.     

Effects of Hall current on hydromagnetic free-convective flow through a porous medium

An analysis of the effects of Hall current on hydromagnetic free-convective flow through a porous medium bounded by a vertical plate is theoretically investigated when a strong magnetic field is imposed in a direction which is perpendicular to the free stream and makes an angle to the vertical direction. The influence of Hall currents on the flow is studied for various values of .Nomenclature c _p specific heat at constant pressure - e electrical charge - E Eckert number - E electrical field intensity - g acceleration due to gravity - G Grashof number - H ₀ applied magnetic field - H magnetic field intensity - (j _x , j _y , j _z ) components of current densityJ - J current density - K permeability of porous medium - M magnetic parameter - m Hall parameter - n _e electron number density - P Prandtl number - q velocity vector - (T, T _w , T ) temperature - t time - (u, v, w) components of the velocity vectorq - U ₀ uniform velocity - v ₀ suction velocity - (x, y, z) Cartesian coordinates Greek Symbols angle - coefficient of volume expansion - _e cyclotron frequency - frequency - dimensionless temperature - thermal conductivity - coefficient of viscosity - magnetic permeability - kinematic viscosity - mass density of fluid - _e charge density - electrical conductivity - _e electron collision time     

The structure of periodic solutions in the restricted problem of three bodies

We consider the basic families of plane-symmetric simply-periodic orbits in the Sun-Jupiter case of the plane restricted three-body problem and we study their horizontal and vertical stabilities. We give the critical orbits of these families, corresponding to the vertical stability parameter = 1 and in future communications we shall give the three-dimensional families which emanate from these plane bifurcations.     

Possible effects of anisotropy ofG on celestial orbits

Will (1971) has discussed a possible anisotropy in the gravitational constantG. Suppose that the attractive gravitational force between two particles of massesm ₁ andm ₂ is given by the usual expressionF=–Gm ₁ m ₂ r/r ³, wherer is the separation vector. Ifc is the velocity of light in vacuo and if 1 _r r/r, he expresses the anisotropy byG=G [1+(v·1 r/c)²], whereG is a constant,v is identified practically as the velocity of the Sun around the galaxy, and 1. Will's suggestion is to look for such an effect in the laboratory.The purpose of the present paper is to look for such an effect in the solar system, wherem ₁ andm ₂ become the masses of the Sun and a planet or of the Earth and the Moon. For simplicity I consider only those planets whose orbits are close to the ecliptic, so that the angle betweenv and the plane of the ecliptic is about 59°.With the above force, the resulting two-body problem is completely solvable. The results are these. If =1, there is an increase in mean motion of 7 parts in 10⁸, a periodic fluctuation in true longitude with period half that of the orbit and amplitude ranging possibly from 0.01 to 0.02, and periodic fluctuations in the radius vector, with period also one half that for the orbit. The amplitudes are: 2.7 km for Mercury, 5.1 km for Venus, 7.0 km for Mars, 18 m for the Moon about the Earth, and 28 cm for a close artificial satellite with inclination 23°. The more conservative estimate <0.0115 would reduce these values by the factor 70.     

Effects of metal abundances on the evolutionary period changes of classical Cepheids

Some peculiarities in the behaviour of a model self-gravitating system described by hydrodynamical equations and isothermal equation of state connected with the presence of thermodynamical fluctuations in real systems were investigated in numerical experiment. The values of density and velocity , , respectively, were computed by numerical code perturbed on each time-step and in each computational cell by random values , for modeling such fluctuations. Perturbed values _i = _i + _i ,v _i = _i + v _i were used to initiate the next step of computations. This procedure is equivalent to an introduction into original hydrodynamical equations of Langevin sources which are random functions. It is shown that these small fluctuations (= v =0,² =v ² = 10^–8) grow many times in marginally-stable state.     

Effects of mass transfer on the unsteady free convection flow past an accelerated vertical porous plate with variable suction

An exact analysis of the effects of mass transfer on the flow of a viscous incompressible fluid past an uniformly accelerated vertical porous and non-porous plate has been presented on taking into account the free convection currents. The results are discussed with the effects of the Grashof number Gr, the modified Grashof number Sc, the Schmidt number Sc, and the suction parametera for Pr (the Prandtl number)=0.71 representating air at 20°C.Nomenclature a suction parameter - C species concentration - C species concentration at the free stream - g acceleration due gravity - Gc modified Grashof number (vg^*(C –C )/U ₀ ³ ) - Pr Prandtl number (C _p/K) - T temperature of the fluid near the plate - T dimensionless temperature near the plate ((T-T )/(T -T )) - U(t) dimensionless velocity of the plate (U/U ₀) - v normal velocity component - v ₀ suction/injection velocity - x, y coordinate along and normal to the plate - v kinematic viscosity (/gr) - C _p specific heat at constant pressure - C _w species concentration at the plate - C non-dimensional species concentration ((C-C )/(C _w -C )) - Gr Grashof number (g(T _w -T )/U ₀ ³ ) - D chemical molecular diffusivity - K thermal conductivity - Sc Schmidt number (/D) - T _w temperature of the plate - T free stream temperature - t time variable - t dimensionless time (tU ₀ ² /) - U ₀ reference velocity - u velocity of the fluid near the plate - u non-dimensional velocity (u/U ₀) - v dimensionless velocity (v/U ₀) - v ₀ non-dimensionalv ₀ (v ₀ /U₀)=–at^–1/2 - y dimensionless ordinate (yU ₀/) - density of the fluid - coefficient of viscosity     

Drift in toroidal configurations

It is found that charged particles of positive energiesE, when constrained on axisymmetric isoflux surfaces , execute sinusoidal motions with typical frequencies =(2E/m)^1/2). In general, it was found that under equilibrium condition p=J ^B/cthe particles develop a non-ambipolar drift velocityv _d =(cµ/eb)[1+q ² +2(q/)²]p.     

Implications of a strongly peaked polar magnetic field

Using the flux-transport equation in the absence of sources, we study the relation between a highly peaked polar magnetic field and the poleward meridional flow that concentrates it. If the maximum flow speed _m greatly exceeds the effective diffusion speed /R, then the field has a quasi-equilibrium configuration in which the poleward convection of flux via meridional flow approximately balances the equatorward spreading via supergranular diffusion. In this case, the flow speed () and the magnetic field B() are related by the steady-state approximation () (/R)B()/B() over a wide range of colatitudes from the poles to midlatitudes. In particular, a general flow profile of the form sin ^p cos ^q which peaks near the equator (q p) will correspond to a cos ⁿ magnetic field at high latitudes only if p = 1 and _m = n /R. Recent measurements of n 8 and 600 km² s^–1 would then give _m 7 m s^–1.     

New kinds of asymmetric periodic orbits in the restricted three-body problem

The procedure of numerical ascent from families of planar to three-dimensional periodic orbits and the subsequent descent to the plane is proved efficient in determining new families of planar asymmetric periodic orbits in the restricted three-body problem. Two such families are computed and described for values of the mass parameter for which it has been found that they exist. Two new families of three-dimensional asymmetric periodic orbits are also presented in this paper.     

Newtonian and relativistic periodic orbits around two fixed black holes

We compare families of simple periodic orbits of test particles in the Newtonian and relativistic problems of two fixed centers (black holes). The Newtonian problem is integrable, while the relativistic problem is highly non-integrable.The orbits are calculated on the meridian plane through the fixed centersM ₁ (atz=+1) andM ₂ (atz=–1) for energies smaller than the escape energyE=1. We use prolate spheroidal coordinates (, , =const) and also the variables =cosh and =–cos . The orbits are inside a curve of zero velocity (CZV). The Newtonian orbits are also limited by an ellipse and a hyperbola, or by two eillipses. There are 3 main types of periodic orbits (1) elliptic type (around both centers), (2) hyperbolic-type, and (3) resonant-type.The elliptic type orbits are stable in the Newtonian case and both stable and unstable in the relativistic case. From the stable orbits bifurcate double period orbits both symmetric and asymmetric with respect to thez-axis. There are also higher order bifurcations. The hyperbolic-type orbits are unstable. The Newtonian resonant orbits are defined by the ratiot _µ/t =n/m of oscillations along and during one period, and they are all marginally unstable. The corresponding relativistic orbits are stable, or unstable. The main families are figure eight orbits aroundM ₁, or aroundM ₂ (3/1 orbits); gamma, or inverse gamma orbits (4/2); higher resonant families 5/1,7/1,...,8/2,12/2,...;, more complicated orbits, like 5/3, and bifurcations from the above orbits. Satellite orbits aroundM ₁, orM ₂, and their bifurcations (e.g. double period) exist in the relativistic case but not in the Newtonian case. The characteristics of the various families are quite different in the Newtonian and the relativistic cases. The sizes of the orbits and their stabilities are also quite different in general. In the Appendix we study the various types of straight line orbits and prove that some subcases introduced by Charlier (1902) are impossible.     

Free-convection flow of water at 4°C past an infinite vertical porous plate with constant heat flux

An analysis of the two-dimensional flow of water at 4°C past an infinite porous plate is presented, when the plate is subjected to a normal suction velocity and the heat flux at the plate is constant. Approximate solutions are derived for the velocity and temperature fields and the skin-friction. The effects ofG (Grashof number) andE (Eckert number) on the velocity and temperature fields are discussed.Nomenclature u, v velocity components of the fluid inx, y direction - g acceleration due to gravity - coefficient of thermal expansion of water at 4°C - v kinematic viscosity - density - T temperature inside thermal boundary layer - T free-stream temperature - k thermal conductivity - C _p specific heat at constant pressure