Newtonian cosmology with a varying gravitational constant

Newtonian cosmology is developed with the assumption that the gravitational constantG diminishes with time. The functional form adopted forG(t), a modification of a suggestion of Dirac, isG=A(k+t) ^–1, wheret is the age of the Universe and a small constantk is inserted to avoid a singularity in the two-body problem. IfR is the scale factor, normalized to unity at an epoch time , the differential equation is then . Here ₀ is the mean density at the epoch time. With the above form forG(t), the solution is reducible to quadratures.The scale factorR either increases indefinitely or has one and only one maximum. LetH ₀ be the present value of Hubble's constant /R and _0c the minimum density for a maximum ofR, i.e., for closure of the Universe. The conditions for a maximum lead to a boundary curve of _0c versusH ₀ and the numbers indicate strongly that thisG-variable Newtonian model corresponds to an open universe. An upward estimate of the age of the Universe from 10¹⁰ yr to five times such a value would still lead to the same conclusion.The present Newtonian cosmology appears to refute the statement, sometimes made, that the Dirac model forG necessarily leads to the conclusion that the age of the Universe is one-third the Hubble time. Appendix B treats this point, explaining that this incorrect conclusion arises from using all the assumptions in Dirac (1938). The present paper uses only Dirac's final result, viz,G(k+t)^–1, superposing it on the differential equation .     

Shear-induced instability and arch filament eruption: A magnetohydrodynamic (MHD) numerical simulation

We investigate, via a two-dimensional (nonplanar) MHD simulation, a situation wherein a bipolar magnetic field embedded in a stratified solar atmosphere (i.e., arch-filament-like structure) undergoes symmetrical shear motion at the footpoints. It was found that the vertical plasma flow velocities grow exponentially leading to a new type of global MHD-instability that could be characterized as a Dynamic Shearing Instability, with a growth rate of about 8{ovV} _A a, where {ovV} _A is the average Alfvén speed and a ^–1 is the characteristic length scale. The growth rate grows almost linearly until it reaches the same order of magnitude as the Alfvén speed. Then a nonlinear MHD instability occurs beyond this point. This simulation indicates the following physical consequences: the central loops are pinched by opposing Lorentz forces, and the outer closed loops stretch upward with the vertically-rising mass flow. This instability may apply to arch filament eruptions (AFE) and coronal mass ejections (CMEs).To illustrate the nonlinear dynamical shearing instability, a numerical example is given for three different values of the plasma beta that span several orders of magnitude. The numerical results were analyzed using a linearized asymptotic approach in which an analytical approximate solution for velocity growth is presented. Finally, this theoretical model is applied to describe the arch filament eruption as well as CMEs.     

Deep drilling at the Siljan Ring impact structure: oxygen-isotope geochemistry of granite

The Siljan Ring is a 362-Ma-old impact structure formed in 1700-Ma-old I-type granites. A 6.8-km-deep borehole provides a vertical profile through granites and isolated horizontal diabase sills. Fluid-inclusion thermometry, and oxygen-isotope compositions of vein quartz, granite, diabase, impact melt, and pseudotachylite, reveal a complex history of fluid activity in the Siljan Ring, much of which can be related to the meteorite impact. In granites from the deep borehole, ¹⁸O values of matrix quartz increase with depth from near 8.0 at the surface to 9.5 at 5760 m depth. In contrast, feldspar ¹⁸O values decrease with depth from near 10 at the surface to 7.1 at 5760 m, forming a pattern opposite to the one defined by quartz isotopic compositions. Values of ¹⁸O for surface granites outside the impact structure are distinct from those in near-surface samples from the deep borehole. In the deep borehole, feldspar coloration varies from brick-red at the surface to white at 5760 m, and the abundances of crack-healing calcite and other secondary minerals decrease over the same interval. Superimposed on the overall decrease in alteration intensity with depth are localized fracture zones at 4662, 5415, and 6044 m depth that contain altered granites, and which provided pathways for deep penetration of surface water. The antithetic variation of quartz and feldspar ¹⁸O values, which can be correlated with mineralogical evidence of alteration, provides evidence for interaction between rocks and impact-heated fluids (100–300° C) in the upper 2 km of the pluton. Penetration of water to depths below 2 km was restricted by a general decrease in impact-fracturing with depth, and by a 60-m-thick diabase sill at 1500 m depth that may have been an aquitard. At depths below 4 km in the pluton, where water/rock ratios were low, oxygen isotopic compositions preserve evidence for limited high-temperature (>500° C) exchange between alkali feldspar and fluids. The high-temperature exchange may have been a post-impact event involving impact-heated fluids, or a post-magmatic event.     

Diffusion-deposition patterns in Martian streaks

The spreading angle of a number of light and dark Martian streaks is determined from selected Mariner 9 images. The resulting frequency distributions of spreading half-angles have maxima at ～5° for light, and ～7° for dark streaks; however the dark streaks have a secondary maximum spreading angle at ～14°. The smaller values, which include most streaks, are interpreted as crater-wake spreading phenomena. The larger value, found in only a few dark streaks or “tails,” may result from atmospheric diffusion and subsequent deposition of material from isolated sources such as vents or blowouts. An atmospheric diffusion-deposition analysis is presented, assuming this streak origin, from which it is possible to deduce the eddy diffusivity, K, in Mars' boudary layer. Calculated K values are found to agree with various theoretical estimates. They lie in the range 10⁷ and 10⁹ cm² sec^?1 and exhibit the proper scale dependence. Thus it appears that, in addition to streak-derived wind direction patterns and speed information, it is possible in a few cases to derive information on Mars' boundary-layer turbulence from streak-spreading measurements.     

Velocity-space maps and transforms of tracking observations for orbital trajectory state analysis

The orbital state of a satellite in a central force field can be uniquely described by its velocity hodograph, a circle, rather than the Keplerian conic. Also, its coordinate-frame rotation about the attracting center is definable, without singularity, by the four-parameter set of Euler parameters. A unified state model of orbital trajectory and attitude dynamics has previously been developed by use of state variables of the orbital velocity hodograph and Euler parameters. The dynamical constraint equations of this orbital state model are especially effective in advanced techniques of state estimation, used for orbit determination and prediction. External observations of orbital vehicles, such as provided by optical and radar sensors of tracking systems, are transformable into corresponding velocity state maps, as presented in this paper. These transformations and the consequent state maps are essential for development of the orbit observation matrix used with the unified state matrix, in recursive estimators such as the Kalman filters. Line-of-sight rays and range spheres (or hemispheres) of observations map conformally into orthogonal spherical surfaces in velocity space, as the result of the point-contact transformations. In bispherical coordinates, the field of observation maps for a ground-based tracking system site is shown to be a reduced (or degenerate) form of the general field of observation maps for a satellite-based tracking site. These orbital state maps and transforms are directly useful in development of observation matrices for candidate observation sets, such as range only, angle only, or range plus range-rate tracking schemes. Also, surface coverage patterns can be generated for proposed new tracking systems, in mission analysis and system synthesis studies.