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 .     

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.     

Formation Mechanism of Mud Bank Along the Southwest Coast of India

Mud bank formation during the southwest monsoon along the southwest coast of India remains an enigma to the researchers and coastal community in spite of several earlier studies. The present study attempts to unravel the mystery through a high-frequency, season-long time-series observation at Alappuzha, located at the southern part of the west coast of India, a region of frequent occurrence of mud bank. Using 7-month-long weekly time-series observation, we identified strong winds and high waves associated with onset of the southwest monsoon and subsequent three episodic atmospheric low-pressure events (LPEs).With the help of in situ time-series data, we show that the strong winds and high waves associated with southwest monsoon pre-conditions the near shore bottom sediment to bring it into suspension. The high amplitude waves associated with the southwest monsoon, while propagating from the deep water to shallow water region, interact with the bottom initiating bottom-sediment movement and its suspension due to wave refraction and shoaling. The sporadic occurrence of the atmospheric LPEs enhances the process of suspension of bottom sediment in the near shore region leading to the formation of fluid mud. Simulations with a cohesive sediment transport model yielded realistic estimates of sediment transport, in the presence of an onshore current, a pre-requisite for transporting the fluid mud toward the coast. The prevailing onshore upwelling current during the southwest monsoon provides the favorable pre-requisite conditions for transporting the fluid mud through depression channel network towards the coast. Once sufficient quantity and thickness of fluid mud is accumulated in the near shore region, it acts as a wave damper for subsequent high monsoon waves, as indicated by the time-series wave data, leading to the formation of tranquil mud bank region. Depression channel networks extending from the shelf to the coast off Alappuzha, Kochi, Ponnani, Beypore, and Ullal were found in the bathymetric charts, thus explaining why mud banks occur only at few locations in spite of the prevalence of similar monsoon conditions.     

An inter-comparison of tidal solutions computed with a range of unstructured grid models of the Irish and Celtic Sea Regions

Three finite element codes, namely TELEMAC, ADCIRC and QUODDY, are used to compute the spatial distributions of the M₂, M₄ and M₆ components of the tide in the sea region off the west coast of Britain. This region is chosen because there is an accurate topographic dataset in the area and detailed open boundary M₂ tidal forcing for driving the model. In addition, accurate solutions (based upon comparisons with extensive observations) using uniform grid finite difference models forced with these open boundary data exist for comparison purposes. By using boundary forcing, bottom topography and bottom drag coefficients identical to those used in an earlier finite difference model, there is no danger of comparing finite element solutions for “untuned unoptimised solutions” with those from a “tuned optimised solution”. In addition, by placing the open boundary in all finite element calculations at the same location as that used in a previous finite difference model and using the same M₂ tidal boundary forcing and water depths, a like with like comparison of solutions derived with the various finite element models was possible. In addition, this open boundary was well removed from the shallow water region, namely the eastern Irish Sea where the higher harmonics were generated. Since these are not included in the open boundary, forcing their generation was determined by physical processes within the models. Consequently, an inter-comparison of these higher harmonics generated by the various finite element codes gives some indication of the degree of variability in the solution particularly in coastal regions from one finite element model to another. Initial calculations using high-resolution near-shore topography in the eastern Irish Sea and including “wetting and drying” showed that M₂ tidal amplitudes and phases in the region computed with TELEMAC were in good agreement with observations. The ADCIRC code gave amplitudes about 30 cm lower and phases about 8° higher. For the M₄ tide, in the eastern Irish Sea amplitudes computed with TELEMAC were about 4 cm higher than ADCIRC on average, with phase differences of order 5°. For the M₆ component, amplitudes and phases showed significant small-scale variability in the eastern Irish Sea, and no clear bias between the models could be found. Although setting a minimum water depth of 5 m in the near-shore region, hence removing wetting and drying, reduced the small-scale variability in the models, the differences in M₂ and M₄ tide between models remained. For M₆, a significant reduction in variability occurred in the eastern Irish Sea when a minimum 5-m water depth was specified. In this case, TELEMAC gave amplitudes that were 1 cm higher and phases 30° lower than ADCIRC on average. For QUODDY in the eastern Irish Sea, average M₂ tidal amplitudes were about 10 cm higher and phase 8° higher than those computed with TELEMAC. For M₄, amplitudes were approximately 2 cm higher with phases of order 15° higher in the northern part of the region and 15° lower in the southern part. For M₆ in the north of the region, amplitudes were 2 cm higher and about 2 cm lower in the south. Very rapid M₆ tidal-phase changes occurred in the near-shore regions. The lessons learned from this model inter-comparison study are summarised in the final section of the paper. In addition, the problems of performing a detailed model–model inter-comparison are discussed, as are the enormous difficulties of conducting a true model skill assessment that would require detailed measurements of tidal boundary forcing, near-shore topography and precise knowledge of bed types and bed forms. Such data are at present not available.     

A robust algorithm for ellipse-based discrete element modelling of granular materials

With recent research indicating the importance of the rolling mechanism of deformation in granular systems consisting of perfectly round particles, it has become popular to use ellipse-shaped particles in the Discrete Element Method (DEM) numerical model. Inherent in this technique is the need for accurately computing ellipse to ellipse intersection, in order to properly detect contact formation and compute relative contact velocities. However, the commonly used algorithms for computing ellipse-ellipse intersection are generally poorly conditioned and can be inaccurate. An alternate method for computing ellipse-ellipse intersection is developed and presented which results in a well-conditioned, stable and accurate contact detection method. These modification are incorporated into the general DEM algorithm.