Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dissipative systems”

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 1,009 records · Page 56

On the prediction of equilibrium states in homogeneous turbulence

A comparison of several commonly used turbulence models (including the Kappa-epsilon and two second-order closures) is made for the test problem of homogeneous turbulent shear flow in a rotating frame. The time evolution of the turbulent kinetic energy and dissipation rate is calculated for a variety of models and comparisons are made with previously published experiments and numerical simulations. Particular emphasis is placed on examining the ability of each model to accurately predict equilibruim states for a range of the parameter Omega/S (the ratio of the rotation rate to the shear rate). It is found that none of the commonly used second-order closure models yield substantially improved predictions for the time evolution of the turbulent kinetic energy and dissipation rate over the somewhat defective results obtained from the simpler Kappa-epsilon model for the turbulent flow regime. There is also a problem with the equilibrium states predicted by the various models. For example, the Kappa-epsilon model erroneously yields equilibrium states that are independent of Omega/S while the Launder, Reece, and Rodi model predicts a flow relaminarization when Omega/S is greater than 0.39 - a result which is contrary to numerical simulations and linear spectral analysis which indicate flow instability for at least the range 0 less than or = Omega/S less than or = 0.5. The physical implications of the results obtained from the various turbulence models considered here are discussed in detail along with proposals to remedy the deficiencies based on a dynamical systems approach.

Speziale, Charles G.↗

On the prediction of equilibrium states in homogeneous turbulence

A comparison of several commonly used turbulence models (including the Kappa-epsilon and two second-order closures) is made for the test problem of homogeneous turbulent shear flow in a rotating frame. The time evolution of the turbulent kinetic energy and dissipation rate is calculated for a variety of models and comparisons are made with previously published experiments and numerical simulations. Particular emphasis is placed on examining the ability of each model to accurately predict equilibrium states for a range of the parameter Omega/S (the ratio of the rotation rate to the shear rate). It is found that none of the commonly used second-order closure models yield substantially improved predictions for the time evolution of the turbulent kinetic energy and dissipation rate over the somewhat defective results obtained from the simpler Kappa-epsilon model for the turbulent flow regime. There is also a problem with the equilibrium states predicted by the various models. For example, the Kappa-epsilon model erroneously yields equilibrium states that are independent of Omega/S while the Launder, Reece, and Rodi model predicts a flow relaminarization when Omega/S is greater than 0.39 - a result which is contrary to numerical simulations and linear spectral analysis which indicate flow instability for at least the range 0 less than or = Omega/S less than or = 0.5. The physical implications of the results obtained from the various turbulence models considered here are discussed in detail along with proposals to remedy the deficiencies based on a dynamical systems approach.

Speziale, C. G.↗

System for the installation and replacement of components in hostile environments

Service equipment for use in hostile environments is presented. The equipment includes a detachable service unit secured to a stationary service unit. The detachable service unit includes a housing with an exterior plate, a power control interface for connection to an exterior power source, locating pins located in said exterior plate, an electrical connector in the exterior plate electrically coupled to said power control interface, and a pair of clamping receptacles formed in the exterior plate and located on adjacent opposite edges of the exterior plate. The stationary unit includes an electrical connector for connection to the electrical connector of the detachable service unit, a clamping apparatus for clamping and unclamping the detachable service unit from the stationary unit, a base clamp assembly for mounting the clamping apparatus onto the stationary unit, and locating pin holes for receiving the locating pins and aligning the detachable service unit onto the stationary unit. The detachable service unit have mating scalloped faces which aid in alignment and provide a mechanism for heat dissipation.

Gordan, Andrew L.↗

Investigation of advanced counterrotation blade configuration concepts for high speed turboprop systems. Task 5: Unsteady counterrotation ducted propfan analysis. Computer program user's manual

The primary objective of this study was the development of a time-marching three-dimensional Euler/Navier-Stokes aerodynamic analysis to predict steady and unsteady compressible transonic flows about ducted and unducted propfan propulsion systems employing multiple blade rows. The computer codes resulting from this study are referred to as ADPAC-AOACR (Advanced Ducted Propfan Analysis Codes-Angle of Attack Coupled Row). This report is intended to serve as a computer program user's manual for the ADPAC-AOACR codes developed under Task 5 of NASA Contract NAS3-25270, Unsteady Counterrotating Ducted Propfan Analysis. The ADPAC-AOACR program is based on a flexible multiple blocked grid discretization scheme permitting coupled 2-D/3-D mesh block solutions with application to a wide variety of geometries. For convenience, several standard mesh block structures are described for turbomachinery applications. Aerodynamic calculations are based on a four-stage Runge-Kutta time-marching finite volume solution technique with added numerical dissipation. Steady flow predictions are accelerated by a multigrid procedure. Numerical calculations are compared with experimental data for several test cases to demonstrate the utility of this approach for predicting the aerodynamics of modern turbomachinery configurations employing multiple blade rows.

Hall, Edward J.↗

Investigation of advanced counterrotation blade configuration concepts for high speed turboprop systems. Task 5: Unsteady counterrotation ducted propfan analysis

The primary objective of this study was the development of a time-marching three-dimensional Euler/Navier-Stokes aerodynamic analysis to predict steady and unsteady compressible transonic flows about ducted and unducted propfan propulsion systems employing multiple blade rows. The computer codes resulting from this study are referred to as ADPAC-AOAR\CR (Advanced Ducted Propfan Analysis Codes-Angle of Attack Coupled Row). This document is the final report describing the theoretical basis and analytical results from the ADPAC-AOACR codes developed under task 5 of NASA Contract NAS3-25270, Unsteady Counterrotating Ducted Propfan Analysis. The ADPAC-AOACR Program is based on a flexible multiple blocked grid discretization scheme permitting coupled 2-D/3-D mesh block solutions with application to a wide variety of geometries. For convenience, several standard mesh block structures are described for turbomachinery applications. Aerodynamic calculations are based on a four-stage Runge-Kutta time-marching finite volume solution technique with added numerical dissipation. Steady flow predictions are accelerated by a multigrid procedure. Numerical calculations are compared with experimental data for several test cases to demonstrate the utility of this approach for predicting the aerodynamics of modern turbomachinery configurations employing multiple blade rows.

Hall, Edward J.↗

Demonstration of Capability to Simulate Particle Irregular Shape and Poly-Disperse Mixtures Within Lunar Lander Plume-Surface Interaction

Plume-Surface Interaction (PSI) between lander engine plumes and regolith soil creates hazards in obscuration and contamination by particle clouds, high-energy ejecta streams, and landing area cratering damage. The MSFC Fluid Dynamics Branch is developing simulation tools to offer a predictive PSI capability to NASA customers such as the Human Lander System (HLS) and Commercial Lunar Payload Services (CLPS). The Gas-Granular Flow Solver (GGFS) is the main application tool for coupled gas-particle two-phase flow simulations to predict the range of PSI effects from onset of surface erosion to deep crater formation. GGFS features an Eulerian-Eulerian modeling approach, treating both gas and granular material as interacting continuum phases. Modeling the lunar regolith granular material fluidic characteristics poses special challenges due to complex particle shapes and mixture composition. The lunar regolith is poorly sorted with broad particle size distributions and large fines content. It has significant cohesion, due to interlocking jagged particle shapes. Eulerian granular material flow modeling requires closure formulations for the granular material constitutive models (stress, friction, collisional and kinetic energy dissipation, drag, etc.). While closure models for spherical particles are available from particle kinetic theory, closure models for realistic non-spherical particles must be extracted from unit physics Discrete Element Model (DEM) particle interaction simulations and provided in the form of tabular datasets. The effects of particle irregular shape (non-spherical shape factors, angular particle surface roughness, and interlocking features) are simulated by approximating the particle features in the form of grouped elemental spheres to form composite particles in the DEM simulations. The effects of the wide range of regolith mixture particle sizes and the strong effects of the presence of the small particle sizes results in high cohesion and low porosity of the regolith mixture. The range of particle sizes is simulated by binning the particle sizes into an appropriate finite number of particle-size species and solving the problem as a species mixture. Combining these two modeling approaches enables simulations to capture both, the contributions of the irregular particle shape and the particle size distribution. The integration and maturation of the DEM-based constitutive model database generation process and poly-disperse mixture binning approach into the GGFS simulation framework are proceeding under funding by the NASA Game Changing Development program. The status of current capabilities will be presented in comparisons of crater characteristics resulting for spherical and irregular shape particles, and for mono-, bi-, and tri-disperse mixture simulations of Apollo LM plume-surface interaction. The computational results confirm the significance of including the particle shape and mixture effects. Going forward plans for the full implementation of the general poly-disperse regolith modeling capability and maturation towards NASA project application readiness under the GCD program will be presented.

Peter A Liever↗

Demonstration of Capability to Simulate Particle Irregular Shape and Poly-Disperse Mixtures Within Lunar Lander Plume-Surface Interaction

Plume-Surface Interaction (PSI) between lander engine plumes and regolith soil creates hazards in obscuration and contamination by particle clouds, high-energy ejecta streams, and landing area cratering damage. The MSFC Fluid Dynamics Branch is developing simulation tools to offer a predictive PSI capability to NASA customers such as the Human Lander System (HLS) and Commercial Lunar Payload Services (CLPS). The Gas-Granular Flow Solver (GGFS) is the main application tool for coupled gas-particle two-phase flow simulations to predict the range of PSI effects from onset of surface erosion to deep crater formation. GGFS features an Eulerian-Eulerian modeling approach, treating both gas and granular material as interacting continuum phases. Modeling the lunar regolith granular material fluidic characteristics poses special challenges due to complex particle shapes and mixture composition. The lunar regolith is poorly sorted with broad particle size distributions and large fines content. It has significant cohesion, due to interlocking jagged particle shapes. Eulerian granular material flow modeling requires closure formulations for the granular material constitutive models (stress, friction, collisional and kinetic energy dissipation, drag, etc.). While closure models for spherical particles are available from particle kinetic theory, closure models for realistic non-spherical particles must be extracted from unit physics Discrete Element Model (DEM) particle interaction simulations and provided in the form of tabular datasets. The effects of particle irregular shape (non-spherical shape factors, angular particle surface roughness, and interlocking features) are simulated by approximating the particle features in the form of grouped elemental spheres to form composite particles in the DEM simulations. The effects of the wide range of regolith mixture particle sizes and the strong effects of the presence of the small particle sizes results in high cohesion and low porosity of the regolith mixture. The range of particle sizes is simulated by binning the particle sizes into an appropriate finite number of particle-size species and solving the problem as a species mixture. Combining these two modeling approaches enables simulations to capture both, the contributions of the irregular particle shape and the particle size distribution. The integration and maturation of the DEM-based constitutive model database generation process and poly-disperse mixture binning approach into the GGFS simulation framework are proceeding under funding by the NASA Game Changing Development program. The status of current capabilities will be presented in comparisons of crater characteristics resulting for spherical and irregular shape particles, and for mono-, bi-, and tri-disperse mixture simulations of Apollo LM plume-surface interaction. The computational results confirm the significance of including the particle shape and mixture effects. Going forward plans for the full implementation of the general poly-disperse regolith modeling capability and maturation towards NASA project application readiness under the GCD program will be presented.

Peter A. Liever↗

Direct Numerical Simulation of Wetting and Spreading Behavior on Heterogeneous and Roughened Substrates

A method of calculation is presented that allows the simulation of the time-dependent three-dimensional motion of thin liquid layers on solid substrates for systems with finite equilibrium contact angles. The contact angle is a prescribed function of position on the substrate. Similar mathematical models are constructed for substrates with a pattern of roughness. Evolution equations are given, using the lubrication approximation, that include viscous, capillary and disjoining forces. Motion to and from dry substrate regions is made possible by use of a thin energetically-stable wetting layer. We simulate motion on heterogeneous substrates with periodic arrays of high contact-angle patches. Two different problems are treated for heterogenous substrates. The first is spontaneous motion driven only by wetting forces. If the contact-angle difference is sufficiently high, the droplet can find several different stable positions, depending on the previous history of the motion. A second simulation treats a forced cyclical motion. Energy dissipation per cycle for a heterogeneous substrate is found to be larger than for a uniform substrate with the same total energy. The Landau-Levich solution for plate removal from a liquid bath is extended to account for a pattern of roughness on the plate.

Schwartz, Leonard W.↗

Strategies for Determining the Cascade Rate in MHD Turbulence: Isotropy, Anisotropy, and Spacecraft Sampling

Exact laws for evaluating cascade rates, tracing back to the Kolmogorov "4/5" law, have been extended to many systems of interest including magnetohydrodynamics (MHD), and compressible flows of the magnetofluid and ordinary fluid types. It is understood that implementations may be limited by the quantity of available data and by the lack of turbulence symmetry. Assessment of the accuracy and feasibility of such third-order (or Yaglom) relations is most effectively accomplished by examining the von Kármán–Howarth equation in increment form, a framework from which the third-order laws are derived as asymptotic approximations. Using this approach, we examine the context of third-order laws for incompressible MHD in some detail. The simplest versions rely on the assumption of isotropy and the presence of a well-defined inertial range, while related procedures generalize the same idea to arbitrary rotational symmetries. Conditions for obtaining correct and accurate values of the dissipation rate from these laws based on several sampling and fitting strategies are investigated using results from simulations. The questions we address are of particular relevance to sampling of solar wind turbulence by one or more spacecraft.

Yanwen Wang↗

Attenuation of Gas Turbulence by a Nearly Stationary Dispersion of Fine Particles

Turbulence attenuation by greater than a factor of two has been observed in many practical gas flows carrying volume fractions as small as 0.01% of dispersed particles. Particles which cause such attenuation usually are smaller than the smallest scales of the turbulence and have time constants 5 to 10 times greater than the time scale of a typical turbulent eddy. That is, strongly attenuating particles usually have Stokes numbers in the range of 5 to 10, indicating that they do not respond to the turbulent fluctuations, but instead just fall through the flow responding only to the mean flow. There are two mechanisms by which free falling particles may attenuate turbulence. First, the unresponsive particles act as a drag on the turbulent eddies, passing energy from the turbulent eddies to the small scale wakes of the particles where it is quickly dissipated by viscosity. The second mechanism is more complicated. Particles falling under gravity convert gravitational potential energy to turbulent velocity fluctuations. If the particles are large, this mechanism increases the overall turbulence level. However, with moderate size particles, the small scale turbulence generated apparently distorts the turbulent eddies leading to more rapid dissipation. Unfortunately, this conclusion is supported only by circumstantial evidence to date. The objectives of the experiment are to use microgravity to separate the two mechanisms. A region of nearly-isotropic decaying turbulence with zero mean flow will be formed in a box in the microgravity environment. Different sets of particles with Stokes numbers in the range of 2 to 20 will be dispersed in the flow. With zero gravity and no mean fluid velocity the particles will have zero mean velocity. With the large Stokes numbers, the fluctuating velocities will also be small. Therefore, the only attenuation mechanism will be the direct action of the particles on the turbulence. Control experiments will also be done in which the particles fall through the measurement volume. Measurements will be acquired using a high resolution image velocimetry (PIV) system being developed specifically for work in particle-laden flows. The measurements will include the decay of the turbulence kinetic energy under various particle loadings. The spatial spectra of the turbulence will also be measured. In a second set of experiments, the interaction of a single eddy with a collection of nearly stationary particles will be examined. The eddy will be a vortex ring emitted by a jet pulse through an orifice. The distortion of the vortex under the influence of the particles will be examined to gain a better understanding of how fine particles can cause such large reductions in turbulence levels. This experiment could not be conducted in terrestrial gravity because the high particle velocities would overwhelm the relatively low speed motion of the vortex ring. This experimental program is just getting underway. The initial challenge is to build a closed facility containing reasonably homogeneous and isotropic turbulence with zero mean velocity. Our approach is to use a set of synthetic jets mounted on the periphery of a transparent plexiglass box to create the turbulence. A synthetic jet is a plenum chamber with an orifice open to the volume of interest. The volume of the chamber fluctuates periodically so alternately a jet is ejected from the volume or flow is drawn back in as a sink. The asymmetry of this situation results in a net transport of momentum and kinetic energy into the volume of interest. The present apparatus includes eight synthetic jets each powered independently by a six inch loudspeaker. The synthetic jets discharge through ejector tubes to increase the scale of the turbulence. Construction of the apparatus is now complete and preliminary flow visualization studies have been conducted. The PIV system is also under development. A compact dual-pulse YAG laser has been acquired as the light source and special software is under development to allow simultaneous measurements of both the particle phase and the fluid phase (marked by fine tracers).

Paris, A. D.↗

Assimilation of SMAP Observations Over Land Improves the Simulation and Prediction of Tropical Cyclone Idai

This work is focused on the role of soil moisture in the prediction of tropical cyclones (TCs) approaching land and after landfall. Soil moisture conditions can impact the circulation and structure of an existing tropical cyclone (TC) when part or all of the circulation is over land. For example, dry land surface conditions may lead to faster dissipation of a TC over land (often associated with changes in precipitation structure), whereas very wet conditions may help sustain or in rare cases re-intensify a TC. Moreover, the presence of strong soil moisture gradients may affect the symmetry and development of the TC circulation leading to changes in its over-land track. While the link between soil moisture conditions and TC evolution in proximity to land is relatively well understood in theory, applications of these findings in the context of numerical weather prediction (NWP) have been limited. Here we present a case study that explores the potential of improving TC predictions through an improved soil moisture initialization in an NWP framework. Specifically, we examine the impact of assimilating observations from the NASA Soil Moisture Active Passive (SMAP) mission into the NASA Goddard Earth Observing System (GEOS) global weather model on the prediction of South-West Indian Ocean TC Idai (2019). SMAP provides accurate L-band (1.4 GHz) brightness temperatures (Tb) observations that are sensitive to soil moisture globally and at high revisit times of 2-3 days. It has previously been shown that the assimilation of SMAP Tb observations significantly improves modeled land surface states. Thus, it is expected that SMAP can be used to constrain land surface initial conditions and potentially benefit TC forecasts. Here we present two sets of retrospective forecasts of TC Idai that are compared in an Observing System Experiment framework at ¼ degree resolution: (i) forecasts initialized from an analysis that is comparable to the GEOS operational analysis (without SMAP Tb assimilation) and (ii) forecasts initialized from an analysis that additionally assimilates SMAP brightness temperature observations over land using a weakly-coupled land analysis. We find that the assimilation of SMAP meaningfully improves the representation of TC Idai’s structure as well as the prediction of its intensity and track. The analyzed TC size, as measured by the wind speed radius, is improved by up to 18% in the analysis with SMAP assimilation relative to the control run. The forecast intensity error, measured against the observed intensity, is reduced by up to 23%. At the 1/4-degree resolution used here, GEOS unavoidably under-estimates TC intensity and over-estimates TC size. The SMAP assimilation therefore corrects the model in the right direction, leading to a storm that is more energetic and more compact. Furthermore, we find that the along-track forecast error is reduced by up to 34%, indicating a more accurate propagation speed, which is consistent with the fact that TC speed over land is strongly affected by surface processes. The impact of SMAP assimilation on the forecast cross-track error is neutral. Across the TC forecast skill metrics used here, the improvements from SMAP DA are largest at lead times of 36 to 72 hours, suggesting that the predictability of forecasts at shorter lead times may be dominated by short-term convective processes, while the land and its longer memory gains in importance as a source of predictability on a 2-3 day timescale. We further investigated the underlying mechanisms leading to the skill improvements from SMAP data assimilation by isolating the land areas that directly influence TC Idai using a back trajectory analysis. We find that the assimilation of SMAP leads to wetter soil moisture conditions that cause an increased latent heat flux, which ultimately results in TC analyzed representation that has higher column-integrated total moisture content and total energy compared to the analysis in the control run without SMAP assimilation. Overall, the results highlight that the assimilation of SMAP observations into a global numerical weather prediction model can lead to pronounced improvements of TC predictions. This is a crucial step towards a better mitigation of the socio-economic impact of landfalling TCs and thus safeguarding human lives. Finally, our study presents an event-based approach that assesses the impact of land data assimilation for a particular weather event rather than by globally averaging differences in skill. We argue that global skill assessments – while necessary – can mute the impact of land data assimilation, because the land’s influence on the atmosphere is constrained to certain locations and certain times. Instead, the event-based approach better highlights the true potential of land data assimilation in the context of NWP, especially for extreme events when accurate predictions are critical.

Jana Kolassa↗

Assimilation of SMAP Observations Over Land Improves the Simulation and Prediction of Tropical Cyclone Idai

This work is focused on the role of soil moisture in the prediction of tropical cyclones (TCs) approaching land and after landfall. Soil moisture conditions can impact the circulation and structure of an existing tropical cyclone (TC) when part or all of the circulation is over land. For example, dry land surface conditions may lead to faster dissipation of a TC over land (often associated with changes in precipitation structure), whereas very wet conditions may help sustain or in rare cases re-intensify a TC. Moreover, the presence of strong soil moisture gradients may affect the symmetry and development of the TC circulation leading to changes in its over-land track. While the link between soil moisture conditions and TC evolution in proximity to land is relatively well understood in theory, applications of these findings in the context of numerical weather prediction (NWP) have been limited. Here we present a case study that explores the potential of improving TC predictions through an improved soil moisture initialization in an NWP framework. Specifically, we examine the impact of assimilating observations from the NASA Soil Moisture Active Passive (SMAP) mission into the NASA Goddard Earth Observing System (GEOS) global weather model on the prediction of South-West Indian Ocean TC Idai (2019). SMAP provides accurate L-band (1.4 GHz) brightness temperatures (Tb) observations that are sensitive to soil moisture globally and at high revisit times of 2-3 days. It has previously been shown that the assimilation of SMAP Tb observations significantly improves modeled land surface states. Thus, it is expected that SMAP can be used to constrain land surface initial conditions and potentially benefit TC forecasts. Here we present two sets of retrospective forecasts of TC Idai that are compared in an Observing System Experiment framework at ¼ degree resolution: (i) forecasts initialized from an analysis that is comparable to the GEOS operational analysis (without SMAP Tb assimilation) and (ii) forecasts initialized from an analysis that additionally assimilates SMAP brightness temperature observations over land using a weakly-coupled land analysis. We find that the assimilation of SMAP meaningfully improves the representation of TC Idai’s structure as well as the prediction of its intensity and track. The analyzed TC size, as measured by the wind speed radius, is improved by up to 18% in the analysis with SMAP assimilation relative to the control run. The forecast intensity error, measured against the observed intensity, is reduced by up to 23%. At the 1/4-degree resolution used here, GEOS unavoidably under-estimates TC intensity and over-estimates TC size. The SMAP assimilation therefore corrects the model in the right direction, leading to a storm that is more energetic and more compact. Furthermore, we find that the along-track forecast error is reduced by up to 34%, indicating a more accurate propagation speed, which is consistent with the fact that TC speed over land is strongly affected by surface processes. The impact of SMAP assimilation on the forecast cross-track error is neutral. Across the TC forecast skill metrics used here, the improvements from SMAP DA are largest at lead times of 36 to 72 hours, suggesting that the predictability of forecasts at shorter lead times may be dominated by short-term convective processes, while the land and its longer memory gains in importance as a source of predictability on a 2-3 day timescale. We further investigated the underlying mechanisms leading to the skill improvements from SMAP data assimilation by isolating the land areas that directly influence TC Idai using a back trajectory analysis. We find that the assimilation of SMAP leads to wetter soil moisture conditions that cause an increased latent heat flux, which ultimately results in TC analyzed representation that has higher column-integrated total moisture content and total energy compared to the analysis in the control run without SMAP assimilation. Overall, the results highlight that the assimilation of SMAP observations into a global numerical weather prediction model can lead to pronounced improvements of TC predictions. This is a crucial step towards a better mitigation of the socio-economic impact of landfalling TCs and thus safeguarding human lives. Finally, our study presents an event-based approach that assesses the impact of land data assimilation for a particular weather event rather than by globally averaging differences in skill. We argue that global skill assessments – while necessary – can mute the impact of land data assimilation, because the land’s influence on the atmosphere is constrained to certain locations and certain times. Instead, the event-based approach better highlights the true potential of land data assimilation in the context of NWP, especially for extreme events when accurate predictions are critical.

Jana Kolassa↗

Iapetus' Geophysics: Rotation Rate, Shape, and Equatorial Ridge

Iapetus has preserved evidence that constrains the modeling of its geophysical history from the time of its accretion until now. The evidence is (a) its present 79.33-day rotation or spin rate, (b) its shape that corresponds to the equilibrium figure for a hydrostatic body rotating with a period of approximately 16 h, and (c) its high, equatorial ridge, which is unique in the Solar System. This paper reports the results of an investigation into the coupling between Iapetus' thermal and orbital evolution for a wide range of conditions including the spatial distributions with time of composition, porosity, short-lived radioactive isotopes (SLRI), and temperature. The thermal model uses conductive heat transfer with temperature-dependent conductivity. Only models with a thick lithosphere and an interior viscosity in the range of about the water ice melting point can explain the observed shape. Short-lived radioactive isotopes provide the heat needed to decrease porosity in Iapetus? early history. This increases thermal conductivity and allows the development of the strong lithosphere that is required to preserve the 16-h rotational shape and the high vertical relief of the topography. Long-lived radioactive isotopes and SLRI raise internal temperatures high enough that significant tidal dissipation can start, and despin Iapetus to synchronous rotation. This occurred several hundred million years after Iapetus formed. The models also constrain the time when Iapetus formed because the successful models are critically dependent upon having just the right amount of heat added by SLRI decay in this early period. The amount of heat available from short-lived radioactivity is not a free parameter but is fixed by the time when Iapetus accreted, by the canonical concentration of Al-26, and, to a lesser extent, by the concentration of Fe-60. The needed amount of heat is available only if Iapetus accreted between 2.5 and 5.0Myr after the formation of the calcium aluminum inclusions as found in meteorites. Models with these features allow us to explain Iapetus? present synchronous rotation, its fossil 16-h shape, and the context within which the equatorial ridge arose.

geophysics↗

A speculation about comets and the earth

Assuming that some 100 earth masses of excess icy conglomerate material accumulated in the Uranus-Neptune region, it is suggested that dispersal by perturbations could have transferred a large fraction of that material into Jupiter-crossing orbits, with a time constant for dissipation of the order of 10 to 100 million years. The possible formation of a 'cometary nebula' within Jupiter's orbit is considered along with the accumulation of a significant amount of cometary material on earth. Calculations are presented which show that the cometary flux within Jupiter may have attained a million times its present value for a period of 10 to 100 million years, that the current flux on earth adds only about one ten-millionth to one hundred-millionth of an earth mass in 4.6 billion years, and that the accretion rate scales only with the local density of the cometary material being collisionally destroyed. The possibility is discussed that the solar wind was unable to cope with the large cometary flux within Jupiter, resulting in a 'choke up' of cometary debris in the inner solar system and providing the cometary nebula with a density many orders of magnitude higher than the present value.

Whipple, F. L.↗

Energy coupling in catastrophic collisions

The prediction of events leading to the catastrophic collisions and disruption of solar system bodies is fraught with the same difficulties as are other theories of impact events; since one simply cannot perform experiments in the regime of interest. In the catastrophic collisions of asteroids that regime involves bodies of a few tons to hundred of kilometers in diameter, and velocities of several kilometers pre second. For hundred kilometer bodies, gravitational stresses dominate material fracture strengths, but those gravitational stresses are essentially absent for laboratory experiments. Only numerical simulations using hydrocodes can in principle analyze the true problems, but they have their own major uncertainties about the correctness of the physical models and properties. The question of the measure of the impactor and its energy coupling is investigated using numerical code calculations. The material model was that of a generic silicate rock, including high pressure melt and vapor phases, and includes material nonlinearity and dissipation via a Mie-Gruniesen model. A series of calculations with various size ratios and impact velocities are reported.

Holsapple, K. A.↗

High Altitude Platform System (HAPS) Communication Support for Wildland Firefighting

High Altitude Platform Systems (HAPS) are emerging aircraft and balloon-type technology that can host payloads and provide services from the stratosphere. One potential HAPS use case is to provide wireless communication services for mobile devices, such as LTE, to wildland firefighters who often operate in locations without terrestrial wireless communications coverage. In this research we analyze historical wildland fire data to provide estimates of the annual number of HAPS required to support a fire season. We apply agglomerative clustering to group historical daily satellite-based fire observations where each cluster is analogous to a required HAPS vehicle. Our lower and upper bound estimates span a range of years, communication payload footprints, the minimum days of clusters prior to launch, and categories of fires. Additionally, we consider a case where HAPS vehicles can be transferred between fires after the initial fire has dissipated. In our specific case study from 2022 “Significant” fires (greater than 40,000 acres), we approximate that either 8 balloon HAPS vehicles without considering overprovisioning for station-keeping limitations or 23 fixed-wing aircraft would be required. Overprovisioning can scale the estimate for balloon vehicles based on reader preference, and for reference, Google Loon overprovisioned by 5-10x. Furthermore, in the case where budgets are constrained and not all of the estimated HAPS vehicles can be acquired, we provide operational insight on where to deploy HAPS vehicles. Generally, in the Spring months we see that HAPS vehicles are needed in the south and southeast of the US which transitions to the north and west as the fire season progresses.

Aaron J. Burns↗

Formation Of the Giant Planets By Concurrent Accretion Of Solids And Gas

New numerical simulations of the formation of the giant planets are presented, in which for the first time both the gas and planetesimal accretion rates are calculated in a self-consistent, interactive fashion. The simulations combine three elements: 1) three-body accretion cross-sections of solids onto an isolated planetary embryo, 2) a stellar evolution code for the planet's gaseous envelope, and 3) a planetesimal dissolution code within the envelope, used to evaluate the planet's effective capture radius and the energy deposition profile of accreted material. Major assumptions include: The planet is embedded in a disk of gas and small planetesimals with locally uniform initial surface mass density, and planetesimals are not allowed to migrate into or out of the planet's feeding zone. All simulations are characterized by three major phases. During the first phase, the planet's mass consists primarily of solid material. The planetesimal accretion rate, which dominates that of gas, rapidly increases owing to runaway accretion, then decreases as the planet's feeding zone is depleted. During the second phase, both solid and gas accretion rates are small and nearly independent of time. The third phase, marked by runaway gas accretion, starts when the solid and gas masses are about equal. It is engendered by a strong positive feedback on the gas accretion rates, driven by the rapid contraction of the gaseous envelope and the rapid expansion of the outer boundary, which depends on the planet's total mass. The overall evolutionary time scale is generally determined by the length of the second phase. The actual rates at which the giant planets accreted small planetesimals is probably intermediate between the constant rates assumed in most previous studies and the highly variable rates that we have used. Within the context, of the adopted model of planetesimal accretion, the joint constraints of the time scale for dissipation of the solar nebula and the current high-Z masses of the giant planets lead to estimates of the initial surface density (sigma(sub init)) of planetesimals in the outer region of the solar nebula. The results show sigma(sub init) approx. = 10 g/sq cm near Jupiter's orbit and that sigma(sub init) proportional to alpha(sup -2), where alpha is the distance from the Sun. These values are a factor of 3 - 4 times as high as that of the "minimum mass" solar nebula at Jupiter's distance and a factor of 2 - 3 times as high it Saturn's distance. Our estimates for the formation time of Jupiter and Saturn are 1 - 10 million years while those for Uranus fall in the range of 2 - 16 million years. These estimates follow from the properties of our Solar System and do not necessarily apply to giant planets in other planetary systems.

Pollack, James B.↗

Exploring Lightning Jump Characteristics

This study is concerned with the characteristics of storms exhibiting an abrupt temporal increase in the total lightning flash rate (i.e., lightning jump, LJ). An automated storm tracking method is used to identify storm "clusters" and total lightning activity from three different lightning detection systems over Oklahoma, northern Alabama and Washington, D.C. On average and for different employed thresholds, the clusters that encompass at least one LJ (LJ1) last longer, relate to higher Maximum Expected Size of Hail, Vertical Integrated Liquid and lightning flash rates (area-normalized) than the clusters that did not exhibit any LJ (LJ0). The respective mean values for LJ1 (LJ0) clusters are 80 min (35 min), 14 mm (8 mm), 25 kg per square meter (18 kg per square meter) and 0.05 flash per min per square kilometer (0.01 flash per min per square kilometer). Furthermore, the LJ1 clusters are also characterized by slower decaying autocorrelation functions, a result that implies a less "random" behavior in the temporal flash rate evolution. In addition, the temporal occurrence of the last LJ provides an estimate of the time remaining to the storm's dissipation. Depending of the LJ strength (i.e., varying thresholds), these values typically range between 20-60 min, with stronger jumps indicating more time until storm decay. This study's results support the hypothesis that the LJ is a proxy for the storm's kinematic and microphysical state rather than a coincidental value.

lightning jump↗