Galactic civilizations: Population dynamics and interstellar diffusion

The interstellar diffusion of galactic civilizations is reexamined by potential theory; both numerical and analytical solutions are derived for the nonlinear partial differential equations which specify a range of relevant models, drawn from blast wave physics, soil science, and, especially, population biology. An essential feature of these models is that, for all civilizations, population growth must be limited by the carrying capacity of the environment. Dispersal is fundamentally a diffusion process; a density-dependent diffusivity describes interstellar emigration. Two models are considered: the first describing zero population growth (ZPG), and the second which also includes local growth and saturation of a planetary population, and for which an asymptotic traveling wave solution is found

Limb-darkening and the structure of the Jovian atmosphere

By observing the transit of various cloud features across the Jovian disk, limb-darkening curves were constructed for three regions in the 4.6 to 5.1 mu cm band. Several models currently employed in describing the radiative or dynamical properties of planetary atmospheres are here examined to understand their implications for limb-darkening. The statistical problem of fitting these models to the observed data is reviewed and methods for applying multiple regression analysis are discussed. Analysis of variance techniques are introduced to test the viability of a given physical process as a cause of the observed limb-darkening

Stress-intensity factor calculations using the boundary force method

The Boundary Force Method (BFM) was formulated for the three fundamental problems of elasticity: the stress boundary value problem, the displacement boundary value problem, and the mixed boundary value problem. Because the BFM is a form of an indirect boundary element method, only the boundaries of the region of interest are modeled. The elasticity solution for the stress distribution due to concentrated forces and a moment applied at an arbitrary point in a cracked infinite plate is used as the fundamental solution. Thus, unlike other boundary element methods, here the crack face need not be modeled as part of the boundary. The formulation of the BFM is described and the accuracy of the method is established by analyzing a center-cracked specimen subjected to mixed boundary conditions and a three-hole cracked configuration subjected to traction boundary conditions. The results obtained are in good agreement with accepted numerical solutions. The method is then used to generate stress-intensity solutions for two common cracked configurations: an edge crack emanating from a semi-elliptical notch, and an edge crack emanating from a V-notch. The BFM is a versatile technique that can be used to obtain very accurate stress intensity factors for complex crack configurations subjected to stress, displacement, or mixed boundary conditions. The method requires a minimal amount of modeling effort

Boundary force method for analyzing two-dimensional cracked bodies

The Boundary Force Method (BFM) was formulated for the two-dimensional stress analysis of complex crack configurations. In this method, only the boundaries of the region of interest are modeled. The boundaries are divided into a finite number of straight-line segments, and at the center of each segment, concentrated forces and a moment are applied. This set of unknown forces and moments is calculated to satisfy the prescribed boundary conditions of the problem. The elasticity solution for the stress distribution due to concentrated forces and a moment applied at an arbitrary point in a cracked infinite plate are used as the fundamental solution. Thus, the crack need not be modeled as part of the boundary. The formulation of the BFM is described and the accuracy of the method is established by analyzing several crack configurations for which accepted stress-intensity factor solutions are known. The crack configurations investigated include mode I and mixed mode (mode I and II) problems. The results obtained are, in general, within + or - 0.5 percent of accurate numerical solutions. The versatility of the method is demonstrated through the analysis of complex crack configurations for which limited or no solutions are known

Bounds for the time to failure of hierarchical systems of fracture

For years limited Monte Carlo simulations have led to the suspicion that the time to failure of hierarchically organized load-transfer models of fracture is non-zero for sets of infinite size. This fact could have a profound significance in engineering practice and also in geophysics. Here, we develop an exact algebraic iterative method to compute the successive time intervals for individual breaking in systems of height $n$ in terms of the information calculated in the previous height $n-1$. As a byproduct of this method, rigorous lower and higher bounds for the time to failure of very large systems are easily obtained. The asymptotic behavior of the resulting lower bound leads to the evidence that the above mentioned suspicion is actually true.Comment: Final version. To appear in Phys. Rev. E, Feb 199

Probabilistic Approach to Time-Dependent Load-Transfer Models of Fracture

A probabilistic method for solving time-dependent load-transfer models of fracture is developed. It is applicable to any rule of load redistribution, i.e, local, hierarchical, etc. In the new method, the fluctuations are generated during the breaking process (annealed randomness) while in the usual method, the random lifetimes are fixed at the beginning (quenched disorder). Both approaches are equivalent.Comment: 13 pages, 4 figures. To appear in Phys.Rev.

Growth and form of the mound in Gale Crater, Mars: Slope wind enhanced erosion and transport

Ancient sediments provide archives of climate and habitability on Mars. Gale Crater, the landing site for the Mars Science Laboratory (MSL), hosts a 5-km-high sedimentary mound (Mount Sharp/Aeolis Mons). Hypotheses for mound formation include evaporitic, lacustrine, fluviodeltaic, and aeolian processes, but the origin and original extent of Gale’s mound is unknown. Here we show new measurements of sedimentary strata within the mound that indicate ∼3° outward dips oriented radially away from the mound center, inconsistent with the first three hypotheses. Moreover, although mounds are widely considered to be erosional remnants of a once crater-filling unit, we find that the Gale mound’s current form is close to its maximal extent. Instead we propose that the mound’s structure, stratigraphy, and current shape can be explained by growth in place near the center of the crater mediated by wind-topography feedbacks. Our model shows how sediment can initially accrete near the crater center far from crater-wall katabatic winds, until the increasing relief of the resulting mound generates mound-flank slope winds strong enough to erode the mound. The slope wind enhanced erosion and transport (SWEET) hypothesis indicates mound formation dominantly by aeolian deposition with limited organic carbon preservation potential, and a relatively limited role for lacustrine and fluvial activity. Morphodynamic feedbacks between wind and topography are widely applicable to a range of sedimentary and ice mounds across the Martian surface, and possibly other planets

A re-evaluation of finite-element models and stress-intensity factors for surface cracks emanating from stress concentrations

A re-evaluation of the 3-D finite-element models and methods used to analyze surface crack at stress concentrations is presented. Previous finite-element models used by Raju and Newman for surface and corner cracks at holes were shown to have ill-shaped elements at the intersection of the hole and crack boundaries. These ill-shaped elements tended to make the model too stiff and, hence, gave lower stress-intensity factors near the hole-crack intersection than models without these elements. Improved models, without these ill-shaped elements, were developed for a surface crack at a circular hole and at a semi-circular edge notch. Stress-intensity factors were calculated by both the nodal-force and virtual-crack-closure methods. Both methods and different models gave essentially the same results. Comparisons made between the previously developed stress-intensity factor equations and the results from the improved models agreed well except for configurations with large notch-radii-to-plate-thickness ratios. Stress-intensity factors for a semi-elliptical surface crack located at the center of a semi-circular edge notch in a plate subjected to remote tensile loadings were calculated using the improved models. The ratio of crack depth to crack length ranged form 0.4 to 2; the ratio of crack depth to plate thickness ranged from 0.2 to 0.8; and the ratio of notch radius to the plate thickness ranged from 1 to 3. The models had about 15,000 degrees-of-freedom. Stress-intensity factors were calculated by using the nodal-force method

A Variational Principle for the Asymptotic Speed of Fronts of the Density Dependent Diffusion--Reaction Equation

We show that the minimal speed for the existence of monotonic fronts of the equation $u_t = (u^m)_{xx} + f(u)$ with $f(0) = f(1) = 0$, $m >1$ and $f>0$ in $(0,1)$ derives from a variational principle. The variational principle allows to calculate, in principle, the exact speed for arbitrary $f$. The case $m=1$ when $f'(0)=0$ is included as an extension of the results.Comment: Latex, postcript figure availabl