Search NASA⌕ Search

SEARCH · Search NASA

Results for “Liquid Breakup”

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 181 records · Page 10

Recent Developments in X-Ray Diagnostics for Cryogenic and Optically Dense Coaxial Rocket Sprays

The mixing and atomization of propellants is often characterized by optically dense flow fields and complex breakup dynamics. In the development of propulsion systems, the complexity of relevant physics and the range of spatio-temporal scales often makes computational simulation impractical for full scale injector elements; consequently, continued research into improved systems for experimental flow diagnostics is ongoing. One area of non-invasive flow diagnostics which has seen widespread growth is using synchrotron based x-ray diagostics. Over the past 3 years, a series of water and cryogenic based experiments were performed at the Advanced Photon Source, Argonne National Lab, on a NASA in-house designed swirl co-axial rocket injector, designed for operation using liquid oxygen and liquid methane in support of Project Morpheus. A range of techniques, such as x-ray fluorescence and time-averaged radiography were performed providing qualitative and quantitative mass and phase distributions, and were complemented by investigations using time-resolved radiography and white beam imaging, which provided information on breakup and mixing dynamics. Results of these investigations are presented, and conclusions regarding the viability of x-ray based diagnostics are discussed.

Radke, Christopher D.↗

Nonlinear Dynamics of a Diffusing Interface

Excitation of two miscible-viscous liquids inside a bounded enclosure in a microgravity environment has shown the evolution of quasi-stationary waves of various modes for a range of parameters. We examine computationally the nonlinear dynamics of the system as the interface breakup and bifurcates to resonance structures typified by the Rayleigh-Taylor instability mechanism. Results show that when the mean steady field is much smaller than the amplitude of the sinusoidal excitation, the system behaves linearly, and growth of quasi-stationary waves occurs through the Kelvin-Helmholtz instability mechanism. However, as the amplitude of excitation increases, nonlinearity occurs through subharmonic bifurcation prior to broadband chaos.

Duval, Walter M. B.↗

Burning of particles undergoing stripping breakup

A simple model for condensed-phase particle burning has been presented. Particles are assumed to start burning in quiescent air. Particles are then accelerated, adding convective effects. Convective effect due to the relative motion between the particle and air on heat and mass transfer is included by introducing the Nusselt number. Stripping of liquid is included at higher particle velocity. Evaporation and burning from the stripped tail are also included.

36 MATERIALS SCIENCE↗

Validation of Ullage Collapse Due to Violent Lateral Slosh

Understanding the coupling between thermodynamics and propellant sloshing is important for cryogenic upper stages, especially during the lift-off phase and in missions including multiple engine restarts, as the fluid dynamic condition of the propellant is mainly influenced by various flight maneuvers. This paper describes the development and validation of a Volume of Fluid (VOF)-based CFD solver, Loci/STREAM-VOF, to the ullage collapse under violent later slosh conditions. Loci/STREAM-VOF is shown to be able to capture detailed slosh dynamics, including slosh frequency, non-linear slosh damping, and violent surface breakup due to forced acceleration. The comparison of surface shape and small drops formed during large amplitude slosh with experimental imagery shows excellent agreement throughout the sloshing period. Comparison of temperature histories in the ullage, at the gas-liquid interface, and inside liquid with the experiment conducted at JAXA also show very good agreement. The simulation points to a stratified temperature field inside the liquid near the interface, critical to the ullage pressure after the collapse. The predicted pressure collapse matches well with the experimental value due to sloshing. Parametric sensitivity studies suggest that high-fidelity Loci/STREM-VOF simulations can play an essential role in predicting the actual thermodynamics of ullage collapse when following the recommended simulation methodology.

H Q Yang↗

Influence of Surface Aerosol Injection on Stratocumulus‐to‐Cumulus Transition: Cloud‐Surface Coupling and Background Aerosol Concentrations

The influence of surface aerosol injection on the stratocumulus‐to‐cumulus transition (SCT) is explored using large‐eddy simulations. We examine how cloud‐surface coupling (or the strength of the marine boundary layer (MBL) stratification that limits vertical turbulent mixing and convection) impacts the vertical transport of aerosols, and how injected aerosols influence cloud properties and associated cloud radiative effects during the SCT. By injecting aerosols at different stages of the SCT, noting that cloud is more decoupled from the surface over time due to entrainment warming, we find that cloud‐surface coupling significantly affects aerosol vertical transport. However, injection timing (before drizzle if any) does not notably affect the SCT and the efficiency of marine cloud brightening, because aerosol number concentrations due to injections at different times rapidly converge before the transition onset. By varying the background aerosol concentration, we find that injected aerosols can significantly extend the persistence of stratocumulus decks by suppressing precipitation in clean environments but have little impact on stratocumulus breakup with higher background aerosol concentrations due to saturated aerosol effects. In clean MBLs, the SCT‐delay‐induced increase in cloud fraction dominates the overall cooling effects in response to aerosols, followed by Twomey effects. These cooling effects are slightly offset by decreased liquid water path (LWP) due to entrainment drying. In polluted MBLs, the Twomey effect is more dominant, followed by cloud fraction adjustments, and these coolings are also partly offset by LWP adjustments. All the simulations are made in relatively small domains in which injected aerosols are homogenized over a short time scale.

Zhang, Haipeng [Univ. of Maryland, College Park, M↗

Bulk Spray and Individual Plume Characterization of LPG and Iso-Octane Sprays at Engine-Like Conditions

This study presents experimental and numerical examination of directly injected (DI) propane and iso-octane, surrogates for liquified petroleum gas (LPG) and gasoline, respectively, at various engine like conditions with the overall objective to establish the baseline with regards to fuel delivery required for future high efficiency DI-LPG fueled heavy-duty engines. Sprays for both iso-octane and propane were characterized and the results from the optical diagnostic techniques including high-speed Schlieren and planar Mie scattering imaging were applied to differentiate the liquid-phase regions and the bulk spray phenomenon from single plume behaviors. The experimental results, coupled with high-fidelity internal nozzle-flow simulations were then used to define best practices in CFD Lagrangian spray models. Optical imaging revealed that unlike iso-octane, propane’s spray propagation was fed by its flash boiling, spray collapse, and high degree of vaporization, resulting in a direct proportionality of propane’s penetration length to temperature. These unique features of propane and its variation from iso-octane’s spray pattern, contributed to its classification as an unconventional spray. Appropriate corrections to the injection and breakup models were developed to reproduce the under-expanded jet dynamics and to mimic the flash boiling-driven spray development observed with propane sprays. The simulation results were found to be sensitive to cone and inclusion angles of the blob injector. The current work represents a first assessment of the capability of the commonly available models for engine-spray simulations and highlights the fact that, despite the reasonable agreement obtained in the fuel vapor morphology, the representation of the liquid phase lacks accuracy and requires further model development.

02 PETROLEUM↗

Roche limit of a solid body

Tidal fission of both impacting and orbiting linear elastic solid bodies based on Kelvin's theory of earth tides is considered. It is shown that there can be more than one mutually exclusive modes of fracture - the particular mode in which a body fractures depending on its size and strength. The analysis gives a vivid picture of the propagation of the fracture with a decreasing distance from the planet. Expressions for the initiation and completion of fracture are obtained which are displayed graphically for a rigid body. The effect of elasticity on the breakup altitude is discussed. For orbiting solid bodies, the study gives the upper limit of the breakup altitude as 0.38R (where R is the radius of planet), which is much less than the value 1.44R used for such bodies in the past. The results presented include a previously given theory by Sekiguchi as a part. For the special case of a liquid body, comparison is made with Roche's calculation and the difference explained.

Aggarwal, H. R.↗

Viscous flow drag reduction; Symposium, Dallas, Tex., November 7, 8, 1979, Technical Papers

The symposium focused on laminar boundary layers, boundary layer stability analysis of a natural laminar flow glove on the F-111 TACT aircraft, drag reduction of an oscillating flat plate with an interface film, electromagnetic precipitation and ducting of particles in turbulent boundary layers, large eddy breakup scheme for turbulent viscous drag reduction, blowing and suction, polymer additives, and compliant surfaces. Topics included influence of environment in laminar boundary layer control, generation rate of turbulent patches in the laminar boundary layer of a submersible, drag reduction of small amplitude rigid surface waves, and hydrodynamic drag and surface deformations generated by liquid flows over flexible surfaces.

Hough, G. R.↗

Coaxial airblast atomizers

An experimental investigation was performed to quantify the characteristics of the sprays of coaxial injectors with particular emphasis on those aspects relevant to the performance of rocket engines. Measurements for coaxial air blast atomizers were obtained using air to represent the gaseous stream and water to represent the liquid stream. A wide range of flow conditions were examined for sprays with and without swirl for gaseous streams. The parameters varied include Weber number, gas flow rate, liquid flow rate, swirl, and nozzle geometry. Measurements were made with a phase Doppler velocimeter. Major conclusions of the study focused upon droplet size as a function of Weber number, effect of gas flow rate on atomization and spray spread, effect of nozzle geometry on atomization and spread, effect of swirl on atomization, spread, jet recirculation and breakup, and secondary atomization.

Hardalupas, Y.↗

Analytical and experimental studies of impinging liquid jets

Impinging injectors are a common type of injector used in liquid propellant rocket engines and are typically used in engines where both propellants are injected as a liquid, e.g., engines using LOX/hydrocarbon and storable propellant combinations. The present research program is focused on providing the requisite fundamental understanding associated with impinging jet injectors for the development of an advanced a priori combustion stability design analysis capability. To date, a systematic study of the atomization characteristics of impinging liquid jets under cold-flow conditions have been completed. Effects of orifice diameter, impingement angle, pre-impingement length, orifice length-to-diameter ratio, fabrication procedure, jet flow condition and jet velocity under steady and oscillating, and atmospheric- and high-pressure environments have been investigated. Results of these experimental studies have been compared to current models of sheet breakup and drop formation. In addition, the research findings have been scrutinized to provide a fundamental explanation for a proven empirical correlation used in the design of stable impinging injector-based rocket engines.

Ryan, H. M.↗

Multiple Environmental Influences on the Lightning of Cold-Based Continental Cumulonimbus Clouds. Part I: Description and Validation of Model

In this two-part paper, influences from environmental factors on lightning in a convective storm are assessed with a model. In Part I, an electrical component is described and applied in the Aerosol–Cloud model (AC). AC treats many types of secondary (e.g., breakup in ice–ice collisions, raindrop-freezing fragmentation, rime splintering) and primary (heterogeneous, homogeneous freezing) ice initiation. AC represents lightning flashes with a statistical treatment of branching from a fractal law constrained by video imagery. The storm simulated is from the Severe Thunderstorm Electrification and Precipitation Study (STEPS; 19/20 June 2000). The simulation was validated microphysically [e.g., ice/droplet concentrations and mean sizes, liquid water content (LWC), reflectivity, surface precipitation] and dynamically (e.g., ascent) in our 2017 paper. Predicted ice concentrations (~10 L -1 ) agreed—to within a factor of about 2—with aircraft data at flight levels (-10° to -15°C). Here, electrical statistics of the same simulation are compared with observations. Flash rates (to within a factor of 2), triggering altitudes and polarity of flashes, and electric fields, all agree with the coincident STEPS observations. The “normal” tripole of charge structure observed during an electrical balloon sounding is reproduced by AC. It is related to reversal of polarity of noninductive charging in ice–ice collisions seen in laboratory experiments when temperature or LWC are varied. Positively charged graupel and negatively charged snow at most midlevels, charged away from the fastest updrafts, is predicted to cause the normal tripole. Total charge separated in the simulated storm is dominated by collisions involving secondary ice from fragmentation in graupel–snow collisions.

54 ENVIRONMENTAL SCIENCES↗

A Nonlinear Model for Fuel Atomization in Spray Combustion

Most gas turbine combustion codes rely on ad-hoc statistical assumptions regarding the outcome of fuel atomization processes. The modeling effort proposed in this project is aimed at developing a realistic model to produce accurate predictions of fuel atomization parameters. The model involves application of the nonlinear stability theory to analyze the instability and subsequent disintegration of the liquid fuel sheet that is produced by fuel injection nozzles in gas turbine combustors. The fuel sheet is atomized into a multiplicity of small drops of large surface area to volume ratio to enhance the evaporation rate and combustion performance. The proposed model will effect predictions of fuel sheet atomization parameters such as drop size, velocity, and orientation as well as sheet penetration depth, breakup time and thickness. These parameters are essential for combustion simulation codes to perform a controlled and optimized design of gas turbine fuel injectors. Optimizing fuel injection processes is crucial to improving combustion efficiency and hence reducing fuel consumption and pollutants emissions.

Liu, Nan-Suey↗

Absolute and Convective Instability of a Liquid Jet

The existence of absolute instability in a liquid jet has been predicted for some time. The disturbance grows in time and propagates both upstream and downstream in an absolutely unstable liquid jet. The image of absolute instability is captured in the NASA 2.2 sec drop tower and reported here. The transition from convective to absolute instability is observed experimentally. The experimental results are compared with the theoretical predictions on the transition Weber number as functions of the Reynolds number. The role of interfacial shear relative to all other relevant forces which cause the onset of jet breakup is explained.

Lin, S. P.↗

Viscoelastic necking dynamics between attractive microgels

Hypothesis: Microgels can deform and interpenetrate and display colloid/polymer duality. The effective interaction of microgels in the collapsed state is governed by the interplay of polymer–solvent interfacial tension and bulk elasticity. A connecting neck is shown to mediate microgel interaction, but its temporal evolution has not been addressed. Here we hypothesize that the necking dynamics of attractive microgels exhibits liquid-like or solid-like behavior over different time and length scales. Experiments: We simulate the merging and pinching of attractive microgels with different crosslinking densities in explicit solvent using dissipative particle dynamics. The temporal coalescence dynamics of microgels is investigated and compared with simple liquid and polymeric droplets. We model the neck growth on long time scales using Maxwell model of polymer relaxation and compare the theoretical prediction with simulation data. The mechanical strength of the neck is characterized systematically via simulated pinch-off of microgels by steered molecular dynamics. Findings: We evidence a crossover in the coalescence dynamics reflecting the viscoelastic signature of microgels. In contrast to the common knowledge that viscoelastic materials respond elastically on short time scales, the early expansion of the microgel neck exhibits a linear behavior, similar to the viscous coalescence of liquid droplets. However, the late regime with arrested dynamics resembles sintering of solid particles. Through an analytical model relating microgel dynamics to neck growth, we show that the long-term behavior is governed by stress relaxation of the polymers in the neck region and predict an exponential decay in the rate of growth, which agrees favorably with the simulation. Different from coalescence, the thread thinning in microgel breakup primarily highlights its polymeric characteristics.

42 ENGINEERING↗

Growth of Hexagonal Boron Nitride from Molten Nickel Solutions: A Reactive Molecular Dynamics Study

Metal flux methods are excellent for synthesizing high-quality hexagonal boron nitride (hBN) crystals, but the atomic mechanisms of hBN nucleation and growth in these systems are poorly understood and difficult to probe experimentally. Here, we harness classical reactive molecular dynamics (ReaxFF) to unravel the mechanisms of hBN synthesis from liquid nickel solvent over time scales up to 30 ns. These simulations mimic experimental conditions by including relatively large liquid nickel slabs containing dissolved boron and a molecular nitrogen gas phase. Overall, the reaction takes place almost exclusively on the surface of the liquid nickel, owing to the low solubility of nitrogen in bulk nickel and the intermediate species’ preference for the metal–gas interface. The formation of hBN invariably begins by reaction of dinitrogen with nickel-solvated boron atoms at the surface, forming intermediate N–N–B species, which typically evolve into B–N–B units through a short-lived intermediate where a single nitrogen atom is coordinated by one nitrogen and two boron atoms. The resulting B–N–B units, in turn, coalesce with growing hBN nuclei and carry nitrogen between hBN nanocrystals in an Ostwald ripening process. The amount of hBN produced on the tens of nanosecond time scale depends critically on the boron concentration, while having a much weaker dependence on the N 2 pressure for the regime considered (N 2 pressures of 2.5–10 MPa, Ni–B solutions with 6–12% boron by atom fraction). The highest rate of hBN formation occurs at the lowest temperature considered (1750 K, just above the melting point of nickel), while no hBN sheets are formed at 2000 K or above. An analysis of the transition pathways for nitrogen atoms shows that the final step, incorporation of small B–N motifs into larger hBN sheets, is the rate-limiting step in the regimes considered. While raising the temperature from 1750 to 2000 K has little effect on the formation of intermediates (N–N–B, B–N–B, etc.), the lack of large hBN sheets at temperatures >1900 K is explained by decreased probability of the final step and increased probability of breakup of hBN into B–N motifs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Cratering on Titan: A Pre-Cassini Perspective

The NASA-ESA Cassini mission, comprising a formidably instrumented orbiter and parachute-borne probe to be launched this October, promises to reveal a crater population on Titan that has been heretofore hidden by atmospheric haze. This population on the largest remaining unexplored surface in the solar system will be invaluable in comparative planetological studies, since it introduces evidence of the atmospheric effects of cratering on an icy satellite. Here, I highlight some impact features we may hope to find and could devote some modeling effort toward. Titan in a Nutshell: Radius= 2575 km. Density= 1880 kg/cubic m consistent with rock-ice composition. Surface pressure = 1.5 bar. Surface gravity = 1.35 m/square s Atmosphere -94% N2 6% CH, Surface temperature = 94K Tropopause temperature = 70K at 40 km alt. Probable liquid hydrocarbon deposits exist on or near the surface.Titan in a Nutshell: Radius= 2575 km. Density= 1880 kg/cubic m consistent with rock-ice composition. Surface pressure = 1.5 bar. Surface gravity = 1.35 m/square s; Atmosphere about 94% N2 6% CH, Surface temperature = 94K Tropopause temperature = 70K at 40 km alt. Probable liquid hydrocarbon deposits exist on or near the surface. Titan is comparable to Callisto and Ganymede for strength/gravity, Mars/Earth/Venus for atmospheric interaction, and Hyperion, Rhea, and Iapetus for impactor distribution. The leading/trailing asymmetry of crater density from heliocentric impactors is expected to be about 5-6, in the absence of resurfacing. Any Saturnocentric impactor population is likely to alter this. In particular the impact disruption of Hyperion is noted; because of the 3:4 orbital resonance with Titan, fragments from the proto-Hyperion breakup would have rapidly accreted onto Titan. Titan's resurfacing history is of course unknown. The disruption of impactors into fragments that individually create small craters is expected to occur. A crude estimate suggests a maximum separation of about 2 km (compared with 4 km on Venus, or 0.5 km on Earth). Crater chains are unlikely on Titan, since impactors must pass close enough to Saturn to be tidally disrupted; as a result, they would suffer aerodynamic disruption. Crater counting on adjacent satellites gives densities of about 200 per 10 (exp 6) square km for 20-km-diameter craters. However, the presence of a thick atmosphere leads to atmospheric shielding, depleting the relative abundance of small craters. This has been evaluated by models, and the relative abundance of small craters may be due to a diagnostic atmospheric collapse. A number of radar-dark "splotches" have been detected on Venus; these have been attributed to the interaction of the surface with the atmospheric shockwave produced by the Tunguska-like explosion of a bolide in the atmosphere. Simple analogy suggests that similar features might occur on Titan, but the shocked mass density (which controls the momentum coupling between the surface and the shockwave) of Titan's cold N2 atmosphere is about 20x smaller than that of Venus's hot CO2 atmosphere. Unless ice is much more easily turned to rubble than is rock, such features seem less probable on Titan. When the energy deposited by an impact forms a fireball with an equilibrate greater than one scale height, the fireball expands upward and can distribute ejecta. on ballistic exoatmospheric trajectories. On Venus this process is believed to be responsible for the parabolic features; the interaction of various-sized particles falling through the atmosphere with the zonal wind field winnows the particles to form a parabolic deposit. Although such a process is possible on Titan, the large scale height at higher altitudes would make it more difficult. Comparison with craters on other icy satellites suggests that craters on Titan will be fairly shallow (depth/diameter about 0.1) and craters greater than 10 km in diameter will have central peaks or domed bases, perhaps with central pits. The formation of ejecta. blankets may involve the atmosphere in a significant way, both by restraining the expansion of the ejecta cloud and by influencing the thermal history of the ejecta. Compared with Venus, Titan's atmosphere will chill an impact melt somewhat quickly, so the long ejecta flows seen on Venus seem less likely, detailed modeling needs to be performed to determine the impact melt production. Crater topography on Titan may be highlighted by the influence of liquids forming crater lakes. Craters with central peaks will typically form ring-shaped lakes, although horseshoe-shaped takes may be common; domed craters with central pits may even form bullseye lakes with islands with central ponds. If liquids have covered a substantial part of Titan's surface for a substantial period, hydroblemes and tsunami deposits may be common.

Lorenz, R. D.↗

Cloud-Resolving ICON Simulations of Secondary Ice Production in Arctic Mixed-Phase Stratocumuli Observed during M-PACE

Field measurements and modeling studies suggest that secondary ice production (SIP) may close the gap between observed Arctic ice nucleating particle (INP) concentrations and ice crystal number concentrations ni. Here, we explore sensitivities with respect to the complexity of different INP parameterizations under the premise that ni is governed by SIP. Idealized, cloud-resolving simulations are performed for the marine cold air outbreak cloud deck sampled during the Mixed-Phase Arctic Cloud Experiment (M-PACE) with the Icosahedral Nonhydrostatic (ICON) model. The impact of the droplet shattering (DS) of raindrops and collisional breakup (BR) in addition to the existing Hallet–Mossop rime splintering mechanism were investigated. Overall, 12 different model experiments (12-h runs) were performed and analyzed. Despite the considerable amount of uncertainty remaining with regard to SIP mechanisms and their process representation in numerical models, we conclude from these experiments that (i) only simulations where DS dominates the SIP signal (potentially amplified by BR) capture observed ice-phase and liquid-phase cloud properties, and (ii) SIP events cluster around the convective outflow region and are structurally linked to mesoscale cloud organization. In addition, interactions with primary nucleation parameterizations of varied complexity were investigated. Here, our simulations show that (i) a stable long-lived mixed-phase cloud (MPC) can be maintained in the absence of primary nucleation once SIP is established, (ii) experiments using a computationally more efficient relaxation-based parameterization of primary nucleation are statistically invariant from simulations considering prognostic INP, and (iii) primary nucleation at cloud-top controls the areal extent of the mixed-phase cloud region, and reduces SIP efficacy via DS due to increased depletion of cloud liquid throughout the entire cloud column.

54 ENVIRONMENTAL SCIENCES↗

Simulation for Prediction of Entry Article Demise (SPEAD): an Analysis Tool for Spacecraft Safety Analysis and Ascent/Reentry Risk Assessment

For the purpose of performing safety analysis and risk assessment for a probable offnominal suborbital/orbital atmospheric reentry resulting in vehicle breakup, a synthesis of trajectory propagation coupled with thermal analysis and the evaluation of node failure is required to predict the sequence of events, the timeline, and the progressive demise of spacecraft components. To provide this capability, the Simulation for Prediction of Entry Article Demise (SPEAD) analysis tool was developed. This report discusses the capabilities, modeling, and validation of the SPEAD analysis tool. SPEAD is applicable for Earth or Mars, with the option for 3 or 6 degrees-of-freedom (DOF) trajectory propagation. The atmosphere and aerodynamics data are supplied in tables, for linear interpolation of up to 4 independent variables. The gravitation model can include up to 20 zonal harmonic coefficients. The modeling of a single motor is available and can be adapted to multiple motors. For thermal analysis, the aerodynamic radiative and free-molecular/continuum convective heating, black-body radiative cooling, conductive heat transfer between adjacent nodes, and node ablation are modeled. In a 6- DOF simulation, the local convective heating on a node is a function of Mach, angle-ofattack, and sideslip angle, and is dependent on 1) the location of the node in the spacecraft and its orientation to the flow modeled by an exposure factor, and 2) the geometries of the spacecraft and the node modeled by a heating factor and convective area. Node failure is evaluated using criteria based on melting temperature, reference heat load, g-load, or a combination of the above. The failure of a liquid propellant tank is evaluated based on burnout flux from nucleate boiling or excess internal pressure. Following a component failure, updates are made as needed to the spacecraft mass and aerodynamic properties, nodal exposure and heating factors, and nodal convective and conductive areas. This allows the trajectory to be propagated seamlessly in a single run, inclusive of the trajectories of components that have separated from the spacecraft. The node ablation simulates the decreasing mass and convective/reference areas, and variable heating factor. A built-in database provides the thermo-mechanical properties of For the purpose of performing safety analysis and risk assessment for a probable offnominal suborbital/orbital atmospheric reentry resulting in vehicle breakup, a synthesis of trajectory propagation coupled with thermal analysis and the evaluation of node failure is required to predict the sequence of events, the timeline, and the progressive demise of spacecraft components. To provide this capability, the Simulation for Prediction of Entry Article Demise (SPEAD) analysis tool was developed. This report discusses the capabilities, modeling, and validation of the SPEAD analysis tool. SPEAD is applicable for Earth or Mars, with the option for 3 or 6 degrees-of-freedom (DOF) trajectory propagation. The atmosphere and aerodynamics data are supplied in tables, for linear interpolation of up to 4 independent variables. The gravitation model can include up to 20 zonal harmonic coefficients. The modeling of a single motor is available and can be adapted to multiple motors. For thermal analysis, the aerodynamic radiative and free-molecular/continuum convective heating, black-body radiative cooling, conductive heat transfer between adjacent nodes, and node ablation are modeled. In a 6- DOF simulation, the local convective heating on a node is a function of Mach, angle-ofattack, and sideslip angle, and is dependent on 1) the location of the node in the spacecraft and its orientation to the flow modeled by an exposure factor, and 2) the geometries of the spacecraft and the node modeled by a heating factor and convective area. Node failure is evaluated using criteria based on melting temperature, reference heat load, g-load, or a combination of the above. The failure of a liquid propellant tank is evaluated based on burnout flux from nucleate boiling or excess internal pressure. Following a component failure, updates are made as needed to the spacecraft mass and aerodynamic properties, nodal exposure and heating factors, and nodal convective and conductive areas. This allows the trajectory to be propagated seamlessly in a single run, inclusive of the trajectories of components that have separated from the spacecraft. The node ablation simulates the decreasing mass and convective/reference areas, and variable heating factor. A built-in database provides the thermo-mechanical properties of

Ling, Lisa↗