Search NASA⌕ Search

SEARCH · Search NASA

Results for “SCALING LAW”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Molecular concentration field design using closed-form steady-state solutions

Control over spatial concentration fields represents a fundamental challenge in designing synthetic biological systems and programmable soft materials. While nature creates morphogen gradients that orchestrate complex developmental processes, synthetic approaches have largely relied on empirical optimization and computationally intensive simulations. Here, we present an analytical framework for steady-state concentration fields generated by finite-sized localized sources in diffusion–degradation systems and derive closed-form solutions for one-, two-, and three-dimensional geometries. By expressing these solutions in dimensionless form, we show that gradient steepness and spatial structure are organized by the Thiele modulus, which captures the competition between diffusion and degradation length scales. The analysis reveals distinct design regimes: in degradation-dominated systems, gradient shape is governed by exponential decay and becomes dimension-independent, whereas in diffusion-dominated systems, gradient magnitude and extent follow dimension-dependent power-law scaling. Building on these results, we introduce a quantitative design strategy that uses threshold-based criteria to program concentration ranges by tuning physically accessible parameters, most directly the production rate, while holding transport and degradation properties fixed. Comparisons with numerical solutions and reported experimental systems demonstrate consistency with the predicted scaling behavior. Together, this work provides a generalizable and physically transparent framework for designing steady-state concentration fields in synthetic biological and soft matter systems, enabling predictive control of gradient-mediated organization without reliance on extensive numerical optimization.

Kim, Dong Woo [Johns Hopkins University, Baltimore↗
DOE OSTI2024

Self-similar Reynolds-averaged mechanical–scalar turbulence models for Rayleigh–Taylor mixing induced by power-law accelerations in the small Atwood number limit

Analytical self-similar solutions to two-, three-, and four-equation Reynolds-averaged mechanical–scalar turbulence models describing turbulent Rayleigh–Taylor mixing driven by a temporal power-law acceleration are derived in the small Atwood number (Boussinesq) limit. The solutions generalize those previously derived for constant acceleration Rayleigh–Taylor mixing for models based on the turbulent kinetic energy K and its dissipation rate ε, together with the scalar variance S and its dissipation rate χ [O. Schilling, “Self-similar Reynolds-averaged mechanical–scalar turbulence models for Rayleigh–Taylor, Richtmyer–Meshkov, and Kelvin–Helmholtz instability-induced mixing in the small Atwood number limit,” Phys. Fluids 33, 085129 (2021)]. The turbulent fields are expressed in terms of the model coefficients and power-law exponent, with their temporal power-law scalings obtained by requiring that the self-similar equations are explicitly time-independent. Mixing layer growth parameters and other physical observables are obtained explicitly as functions of the model coefficients and parameterized by the exponent of the power-law acceleration. Values for physical observables in the constant acceleration case are used to calibrate the two-, three-, and four-equation models, such that the self-similar solutions are consistent with experimental and numerical simulation data corresponding to a canonical (i.e., constant acceleration) Rayleigh–Taylor turbulent flow. The calibrated four-equation model is then used to numerically reconstruct the mean and turbulent fields, and turbulent equation budgets across the mixing layer for several values of the power-law exponent. Finally, the reference solutions derived here can be used to understand the model predictions for strongly accelerated or decelerated Rayleigh–Taylor mixing in the large Reynolds number limit.

42 ENGINEERING↗
DOE OSTI2025

A study of resistive peeling–ballooning modes across spherical tokamaks

We investigate how non-ideal-magnetohydrodynamics (MHD) effects, in particular plasma resistivity, impact the peeling–ballooning stability thresholds in spherical tokamaks. This analysis follows the discovery of resistive kink-peeling modes in ELMing National Spherical Torus Experiment (NSTX) discharges. In the present study we extend this modeling to ELMing pulses in the Mega Ampere Spherical Tokamak (MAST) and MAST—Upgrade (MAST-U), where we find a clear resistive scaling for peeling–ballooning modes. While in NSTX ideal-MHD predicts stability for ELMing discharges, in MAST-U we find that the plasma is slightly unstable to peeling–ballooning modes, but is fully stabilized once diamagnetic effects are considered in terms of a growth rate normalization. A resistive power law scaling is calculated for these modes on MAST-U, which lies in between that of tearing modes and resistive interchange modes. A comparison between M3D-C1 and NIMROD shows reasonable agreement for this scaling. Resistivity destabilizes the modes and the peeling–ballooning unstable domain is considerably expanded in both, MAST and MAST-U. In addition to the MAST/-U pulses we also analyze resistive PB stability in a NSTX-similarity discharge on DIII-D. While having a different aspect ratio from NSTX, this discharge uses NSTX-like shaping parameters, toroidal field and plasma current. By considering these discharges alongside NSTX cases, we identify conditions influencing the onset of resistive peeling–ballooning modes. Furthermore, our findings indicate that magnetic shear in the pedestal region is closely linked to the emergence of resistive edge modes.

MAST-U↗
DOE OSTI2024

Interplay of intercalation dynamics and lithium plating in monolithic and architectured graphite anodes during fast charging

Fast charging of high-capacity anodes is challenging due to lithium plating reactions, which lead to poor cycling performance and safety concerns. Thus, accurate predictions of plating onset and an understanding of this electrochemical process are crucial for robust battery design. However, the most commonly used models, based on porous electrode theory (e.g., the pseudo-2D model), are notoriously difficult to calibrate due to their complexity, limiting their predictive power. This work studies the process of lithium plating during fast charging of (small-particle) graphite half-cells by measuring local reaction progression and plating behavior using optical operando techniques. These experiments employ a realistic 1D graphite electrode geometry with commercially-relevant mass loading charged at fast charge rates. It is demonstrated that the local reaction progression and plating onset can not only be predicted accurately with a p2D numerical model, but that these processes follow a simple scaling law. Remarkably, the entire reaction histories of different electrodes charged at different rates (e.g., 160 μm thickness at 0.5C, 111 μm at 1C or 66 μm at 4C) were observed to have self-similar intercalation profiles. It is demonstrated that plating onset is in turn governed by the reaction profile which explains why both processes exhibit the same scaling behavior. Finally, operando measurements of local reaction dynamics are conducted for the first time in electrodes with channeled architectures, quantitatively determining how channels affect reaction uniformity and plating onset. Together, these results reveal underlying simplicity in the complex electrochemical environment of fast charging and lithium plating, improving understanding of this process. These fundamental insights are broadly applicable for design processes, modeling and experimental evaluation of lithium ion batteries.

25 ENERGY STORAGE↗
DOE OSTI2024

Generalized vs. fractional: a comparative analysis of Maxwell models applied to entangled polymer solutions

Fractional viscoelastic models provide an excellent description of rheological data for polymer systems with power-law behaviour. However, the physical interpretation of their model parameters, which carry fractional units of time, often remains elusive. We show that for poly(ethylene oxide) (PEO) solutions, the fractional Maxwell model (FMM) requires much fewer model parameters than the classical generalized Maxwell model for a good description of the data and that it can be applied consistently to solutions with varying degrees of viscoelasticity. Despite their fractional units, the parameters exhibit scaling laws similar to classical parameters as a function of polymer concentration.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

ν-point energy correletors with F AST EEC: Small-x physics from LHC jets

In recent years, energy correlators have emerged as a powerful tool for studying jet substructure, with promising applications such as probing the hadronization transition, analyzing the quark-gluon plasma, and improving the precision of top quark mass measurements. The projected N-point correlator measures correlations between N final-state particles by tracking the largest separation between them, showing a scaling behavior related to DGLAP splitting functions. These correlators can be analytically continued in N, commonly referred to as ν-correlators, allowing access to non-integer moments of the splitting functions. Of particular interest is the ν → 0 limit, where the small momentum fraction behavior of the splitting functions requires resummation. Originally, the computational complexity of evaluating ν-correlators for M particles scaled as 2 2M , making it impractical for real-world analyses. However, by using recursion, we reduce this to M 2M , and through the FastEEC method of dynamically resolving subjets, M is replaced by the number of subjets. This breakthrough enables, for the first time, the computation of ν-correlators for LHC data. In practice, limiting the number of subjets to 16 is sufficient to achieve percent-level precision, which we validate using known integer-ν results and convergence tests for non-integer ν. We have implemented this in an update to FastEEC and conducted an initial study of power-law scaling in the perturbative regime as a function of ν, using CMS Open Data on jets. The results agree with DGLAP evolution, except at small ν, where the anomalous dimension saturates to a value that matches the BFKL anomalous dimension.

Energy correlators↗
DOE OSTI2025

Cosmological quasiparticles and the cosmological collider

The interplay between cosmology and strongly coupled dynamics can yield transient spectral features that vanish at late times, but which may leave behind phenomenological signatures in the spectrum of primordial fluctuations. Of particular interest are strongly coupled extensions of the standard model featuring approximate conformal invariance. In flat space, the spectral density for a scalar operator in a conformal field theory is characterized by a continuum with scaling law governed by the dimension of the operator, and is otherwise featureless. Anti–de Sitter/conformal field theory arguments suggest that for large 𝑁, in an inflationary background with Hubble rate 𝐻, this continuum is gapped. We demonstrate that there can be additional peak structures that become sharp and particlelike at phenomenologically interesting regions in parameter space, and we estimate their contribution to cosmological observables. We find phenomena that are potentially observable in future experiments that are unique to these models, including displaced oscillatory features in the squeezed limit of the bispectrum. These particles can be either fundamental, and localized to a UV brane, or composite at the Hubble scale, 𝐻, and bound to a horizon in the bulk of the 5D geometry. We comment on how stabilization of conformal symmetry breaking vacua can be correlated with these spectral features and their phenomenology.

conformal field theory↗
DOE OSTI2024

Shock wave formation in radiative plasmas

The temporal evolution of weak shocks in radiative media is theoretically investigated in this work. The structure of radiative shocks has traditionally been studied in a stationary framework. Their systematic classification is complex because layers of optically thick and thin regions alternate to form a radiatively-driven precursor and a temperature-relaxation layer, between which the hydrodynamic shock is embedded. In this work, we analyze the formation of weak shocks when two radiative plasmas with different pressures are put in contact. Applying a reductive perturbative method yields a Burgers-type equation that governs the temporal evolution of the perturbed variables including the radiation field. The conditions upon which optically thick and thin solutions exist have been derived and expressed as a function of the shock strength and Boltzmann number. Below a certain Boltzmann number threshold, weak shocks always become optically thick asymptotically in time, while thin solutions appear as transitory structures. The existence of an optically thin regime is related to the presence of an overdense layer in the compressed material. Scaling laws for the characteristic formation time and shock width are provided for each regime. The theoretical analysis is supported by FLASH simulations, and a comprehensive testcase has been designed to benchmark radiative hydrodynamic codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Are turbulence effects on droplet collision–coalescence a key to understanding observed rain formation in clouds?

Rain formation is a critical factor governing the lifecycle and radiative forcing of clouds and therefore it is a key element of weather and climate. Cloud microphysics–turbulence interactions occur across a wide range of scales and are challenging to represent in atmospheric models with limited resolution. Based on past experiments and idealized numerical simulations, it has been postulated that cloud turbulence accelerates rain formation by enhancing drop collision–coalescence. We provide substantial evidence for significant impacts of turbulence on the evolution of cloud droplet size distributions and rain formation by comparing high-resolution observations of cumulus congestus clouds with state-of-the-art large-eddy simulations coupled with a Lagrangian particle-based microphysics scheme. Turbulent coalescence must be included in the model to accurately represent the observed drop size distributions, especially for drizzle drop sizes at lower heights in the cloud. Turbulence causes earlier rain formation and greater rain accumulation compared to simulations with gravitational coalescence only. The observed rain size distribution tail just above cloud base follows a power law scaling that deviates from theoretical scalings considering either a purely gravitation collision kernel or a turbulent kernel neglecting droplet inertial effects, providing additional evidence for turbulent coalescence in clouds. In contrast, large aerosols acting as cloud condensation nuclei (“giant CCN”) do not significantly impact rain formation owing to their long timescale to reach equilibrium wet size relative to the lifetime of rising cumulus thermals. Overall, turbulent drop coalescence exerts a dominant influence on rain initiation in warm cumulus clouds, with limited impacts of giant CCN.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Fundamental bandwidth limits and shaping of frequency-modulated combs

Frequency-modulated (FM) combs based on active cavities like quantum cascade lasers have recently emerged as promising light sources in many spectral regions. Unlike passive modelocking, which generates amplitude modulation using the field’s amplitude, FM comb formation relies on the generation of phase modulation from the field’s phase. They can therefore be regarded as a phase-domain version of passive modelocking. However, while the ultimate scaling laws of passive modelocking have long been known—Haus showed in 1975 that pulses modelocked by a fast saturable absorber have a bandwidth proportional to effective gain bandwidth—the limits of FM combs have been much less clear. Here, we show that FM combs based on fast gain media are governed by the same fundamental limits, producing combs whose bandwidths are linear in the effective gain bandwidth. Not only do we show theoretically that the diffusive effect of gain curvature limits comb bandwidth, but we also show experimentally how this limit can be increased. By adding carefully designed resonant-loss structures that are evanescently coupled to the cavity of a terahertz laser, we reduce the curvature and increase the effective gain bandwidth of the laser, demonstrating bandwidth enhancement. Our results can better enable the creation of active chip-scale combs and be applied to a wide array of cavity geometries.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Overview of ASDEX upgrade results in view of ITER and DEMO

Experiments on ASDEX Upgrade (AUG) in 2021 and 2022 have addressed a number of critical issues for ITER and EU DEMO. A major objective of the AUG programme is to shed light on the underlying physics of confinement, stability, and plasma exhaust in order to allow reliable extrapolation of results obtained on present day machines to these reactor-grade devices. Concerning pedestal physics, the mitigation of edge localised modes (ELMs) using resonant magnetic perturbations (RMPs) was found to be consistent with a reduction of the linear peeling-ballooning stability threshold due to the helical deformation of the plasma. Conversely, ELM suppression by RMPs is ascribed to an increased pedestal transport that keeps the plasma away from this boundary. Candidates for this increased transport are locally enhanced turbulence and a locked magnetic island in the pedestal. The enhanced D-alpha (EDA) and quasi-continuous exhaust (QCE) regimes have been established as promising ELM-free scenarios. Here, the pressure gradient at the foot of the H-mode pedestal is reduced by a quasi-coherent mode, consistent with violation of the high-n ballooning mode stability limit there. This is suggestive that the EDA and QCE regimes have a common underlying physics origin. In the area of transport physics, full radius models for both L- and H-modes have been developed. These models predict energy confinement in AUG better than the commonly used global scaling laws, representing a large step towards the goal of predictive capability. A new momentum transport analysis framework has been developed that provides access to the intrinsic torque in the plasma core. In the field of exhaust, the X-Point Radiator (XPR), a cold and dense plasma region on closed flux surfaces close to the X-point, was described by an analytical model that provides an understanding of its formation as well as its stability, i.e., the conditions under which it transitions into a deleterious MARFE with the potential to result in a disruptive termination. With the XPR close to the divertor target, a new detached divertor concept, the compact radiative divertor, was developed. Here, the exhaust power is radiated before reaching the target, allowing close proximity of the X-point to the target. No limitations by the shallow field line angle due to the large flux expansion were observed, and sufficient compression of neutral density was demonstrated. With respect to the pumping of non-recycling impurities, the divertor enrichment was found to mainly depend on the ionisation energy of the impurity under consideration. In the area of MHD physics, analysis of the hot plasma core motion in sawtooth crashes showed good agreement with nonlinear 2-fluid simulations. This indicates that the fast reconnection observed in these events is adequately described including the pressure gradient and the electron inertia in the parallel Ohm’s law. Concerning disruption physics, a shattered pellet injection system was installed in collaboration with the ITER International Organisation. Thanks to the ability to vary the shard size distribution independently of the injection velocity, as well as its impurity admixture, it was possible to tailor the current quench rate, which is an important requirement for future large devices such as ITER. Progress was also made modelling the force reduction of VDEs induced by massive gas injection on AUG. The H-mode density limit was characterised in terms of safe operational space with a newly developed active feedback control method that allowed the stability boundary to be probed several times within a single discharge without inducing a disruptive termination. Regarding integrated operation scenarios, the role of density peaking in the confinement of the ITER baseline scenario (high plasma current) was clarified. The usual energy confinement scaling ITER98(p,y) does not capture this effect, but the more recent H20 scaling does, highlighting again the importance of developing adequate physics based models. Advanced tokamak scenarios, aiming at large non-inductive current fraction due to non-standard profiles of the safety factor in combination with high normalised plasma pressure were studied with a focus on their access conditions. A method to guide the approach of the targeted safety factor profiles was developed, and the conditions for achieving good confinement were clarified. Based on this, two types of advanced scenarios (‘hybrid’ and ‘elevated’

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Long‐Range Non‐Coulombic Coupling at the LaAlO 3 /SrTiO 3 Interface

The LaAlO 3 /SrTiO 3 interface hosts a plethora of gate‐tunable electronic phases. Gating of LaAlO 3 /SrTiO 3 interfaces is usually assumed to occur electrostatically. However, increasing evidence suggests that non‐local interactions can influence and, in some cases, dominate the coupling between applied gate voltages and electronic properties. Here, quasi‐1D ballistic electron waveguides are sketched at the LaAlO 3 /SrTiO 3 interface as a probe to understand how gate tunability varies as a function of spatial separation. Gate tunability measurements reveal the scaling law to be at odds with the pure electrostatic coupling observed in traditional semiconductor systems. The non‐Coulombic gating at the interface is attributed to a long‐range nanoelectromechanical coupling between the gate and electron waveguide, possibly mediated by the ferroelastic domains in SrTiO 3 . The long‐range interactions at the LaAlO 3 /SrTiO 3 interface add unexpected richness and complexity to this correlated electron system.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Shear rate dependency on flowing granular biomass material

The commercialization of bioenergy has been significantly limited by various material handling issues due to the poor flowability of granular biomass materials. A good understanding of flow physics and robust constitutive models to predict flow behavior across multiple regimes are essential to address these issues. In this study, we investigated the multi-regime flow behavior of loblolly pine chips, a widely used bioenergy feedstock, through comprehensive inclined plane flow experiments and simulations. A quasi-static hypoplastic model and a cross-regime Drucker-Prager-µ(I) model were calibrated and validated against the physical experiments to investigate the quasi-static and dense flow behavior. The results show that for granular biomass, 1) plane flow (iso-thickness along the plane) exists within a smaller range of inclination while heap flow (varying thickness along the plane) exists within a broader range of inclination, as compared with conventional granular materials (e.g., glass beads); 2) the scaling law of granular biomass flowing on an inclined plane (Froude Number versus dimensionless thickness) forms a bi-linear trend with the turning point governed by the quasi-static and dense flow regimes; 3) the multi-regime DP-µ(I) model can capture the flow behavior in both regimes well at the cost of extra calibration. In conclusion, these findings advance the scientific understanding of the multi-regime flow behavior of granular biomass materials and shed light on formulating novel constitutive models to assist granular biomass handling in the bioenergy industry.

09 BIOMASS FUELS↗
DOE OSTI2024

Impulse-driven release of gas-encapsulated drops

Gas-encapsulated drops, much like antibubbles, are drops enclosed in a bubble within a liquid. They show potential as payload carriers in fluid transport and mixing techniques where sound waves can be leveraged to induce the collapse of the gas core and the subsequent release of the drop. Here, the interaction of millimetre-sized gas-encapsulated drops with impulsive laser-induced shock waves is investigated to gain fundamental insights into the release process. Experimental synchrotron X-ray phase contrast imaging, which allows the drop dynamics to be visualised inside the encapsulating bubble, is complemented by numerical simulations to study the intricate physics at play. Three drop dynamical release regimes are discovered, namely the drop impact , partial deposition and jet impact regimes. The regime type is mainly dependent on the shape of the bubble interface impacting the drop and the associated Weber and Reynolds numbers. The drop dynamics of the drop impact and partial deposition regimes show similarities with the canonical configuration of drops impacting flat liquid surfaces, whereas the jet impact regime resembles binary drop collisions, which allows existing scaling laws to be applied to describe the underlying processes. The release of the drop is investigated numerically. The time evolution of the drop dissemination within the surrounding liquid discloses enhanced mixing for dynamics involving high Weber and Reynolds numbers such as the drop impact and jet impact regimes.

42 ENGINEERING↗
DOE OSTI2024

Magneto-Stokes flow in a shallow free-surface annulus

In this study, we analyse ‘magneto-Stokes’ flow, a fundamental magnetohydrodynamic (MHD) flow that shares the cylindrical-annular geometry of the Taylor–Couette cell but uses applied electromagnetic forces to circulate a free-surface layer of electrolyte at low Reynolds numbers. The first complete, analytical solution for time-dependent magneto-Stokes flow is presented and validated with coupled laboratory and numerical experiments. Three regimes are distinguished (shallow-layer, transitional and deep-layer flow regimes), and their influence on the efficiency of microscale mixing is clarified. The solution in the shallow-layer limit belongs to a newly identified class of MHD potential flows, and thus induces mixing without the aid of axial vorticity. We show that these shallow-layer magneto-Stokes flows can still augment mixing in distinct Taylor dispersion and advection-dominated mixing regimes. The existence of enhanced mixing across all three distinguished flow regimes is predicted by asymptotic scaling laws and supported by three-dimensional numerical simulations. Mixing enhancement is initiated with the least electromagnetic forcing in channels with order-unity depth-to-gap-width ratios. If the strength of the electromagnetic forcing is not a constraint, then shallow-layer flows can still yield the shortest mixing times in the advection-dominated limit. Our robust description of momentum evolution and mixing of passive tracers makes the annular magneto-Stokes system fit for use as an MHD reference flow.

58 GEOSCIENCES↗
DOE OSTI2025

Shock-induced bubble jets: a dual perspective of bubble collapse and interfacial instability theory

Interactions between shock waves and gas bubbles in a liquid can lead to bubble collapse and high-speed liquid jet formation, relevant to biomedical applications such as shock wave lithotripsy and targeted drug delivery. This study reveals a complex interplay between acceleration-induced instabilities that drive jet formation and radial accelerations causing overall bubble collapse under shock wave pressure. Using high-speed synchrotron X-ray phase contrast imaging, the dynamics of micrometre-sized air bubbles interacting with laser-induced underwater shock waves are visualised. These images offer full optical access to phase discontinuities along the X-ray path, including jet formation, its propagation inside the bubble, and penetration through the distal side. Jet formation from laser-induced shock waves is suggested to be an acceleration-driven process. A model predicting jet speed based on the perturbation growth rate of a single-mode Richtmyer–Meshkov instability shows good agreement with experimental data, despite uncertainties in the jet-driving mechanisms. The jet initially follows a linear growth phase, transitioning into a nonlinear regime as it evolves. To capture this transition, a heuristic model bridging the linear and nonlinear growth phases is introduced, also approximating jet shape as a single-mode instability, again matching experimental observations. Upon piercing the distal bubble surface, jets can entrain gas and form a toroidal secondary bubble. Linear scaling laws are identified for the pinch-off time and volume of the ejected bubble relative to the jet’s Weber number, characterising the balance of inertia and surface tension. At low speeds, jets destabilise due to capillary effects, resulting in ligament pinch-off.

Drops and Bubbles: Bubble dynamics↗
DOE OSTI2025

Radiative Heat Transfer and 2D Transition Metal Dichalcogenide Materials

Here we study the radiative heat transfer power in the family of transition metal dichalcogenide monolayers in their H- and T-symmetries. For this purpose, the electronic and optical properties computed from first-principles are used in effective models to understand the emerging scaling laws for metals and semiconductors as well as specific material signatures as control knobs for radiative heat transfer. Our combined approach of analytical modeling with properties from ab initio simulations can be used for other materials families to build a materials database for radiative heat transfer.

36 MATERIALS SCIENCE↗
DOE OSTI2026

A Panspermia Origin for Venus Cloud Life

Decades of study have hinted at the astrobiological potential of Venus's cloud layers. This potential is often cast as stemming from the idea that the Venusian surface was clement in the past. As the climate changed, life then remained in, or perhaps evolved and migrated to, the last habitable niche: the altitudes above ∼50 km with Earth-like temperatures and pressures today. Here we explore an alternative scenario where life was delivered to Venus' clouds from Earth or Mars (“panspermia”). This process requires a life-containing bolide to enter the atmosphere, without experiencing complete sterilization, and then be dispersed at high altitude in fragments small enough to dwell in the clouds. We adapt a widely used model of bolide-atmosphere interaction to investigate the fate of bolides delivered to Venus from Earth and Mars. Starting at the top of the atmosphere, bolides ablate and fragment. Aerodynamic drag spreads these fragments horizontally, forming a “pancake” with an increased effective cross-section, causing rapid deceleration. An airburst occurs when the bolide deposits its highest amount of kinetic energy in the atmosphere. Observations of terrestrial meteorites provide a scaling law for the distribution of post-airburst fragment sizes. Inspired by the “Venus Life Equation,” we present a framework for calculating the rate at which panspermia delivers microbial life to the clouds of Venus. Our best estimate is an average of ∼100 cells dispersed in the clouds per Earth-year. Whether this life can survive and thrive in its new home remains an open question.

Guinan, Emma [Arizona State Univ., Tempe, AZ (Unit↗