Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “surface mass variations”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
NASA NTRS1980

Lunar cold traps and their influence on argon-40

In polar areas of the moon the maximum temperatures reached in some permanently shaded areas are well below the temperature required to retain water ice for billions of years, and cold enough to hold other volatiles for shorter periods. Aside from water, the most significant lunar volatiles are the radiogenic gases, of which argon-40 is the most easily detected, both in situ and as retrapped ions in rocks returned from the surface on the moon. Argon-40 escapes from the moon at a surprisingly high rate that is between 3% and 6% of its total production. Its brief lifetime in the lunar exosphere is marked by numerous adsorption/desorption events. Collisions with the lunar surface in cold, permanently shaded areas lead to long term storage, forming reservoirs of trapped gas that may be disturbed occasionally to produce sudden increases in atmospheric argon. It is postulated that this may explain at least part of the time variations in Apollo 17 mass spectrometer measurements of argon that were previously attributed to internal processes associated with the release of radiogenic gases from the moon.

Hodges, R. R., Jr.↗
NASA NTRS1986

Magnetic connection for Saturn's rings and atmosphere

Latitudinal variations in images of Saturn's disk, upper atmospheric temperatures, and ionospheric electron densities are found in magnetic conjugacy with features in Saturn's ring plane. It is proposed that these latitudinal variations are the result of a variable influx of water transported along magnetic field lines from sources in Saturn's ring plane. These features are thus the surface expression of an electromagnetic erosion mechanism which transports water (in the form of high charge-to-mass ratio particles) from the rings to the atmosphere.

Connerney, J. E. P.↗
NASA NTRS2019

Assessing Uncertainties of a Geophysical Approach to Estimate Surface Fine Particulate Matter Distributions from Satellite-Observed Aerosol Optical Depth

Health impact analyses are increasingly tapping the broad spatial coverage of satellite aerosol optical depth (AOD) products to estimate human exposure to fine particulate matter (PM2.5). We use a forward geophysical approach to derive ground-level PM2.5 distributions from satellite AOD at 1 km2(exp) resolution for 2011 over the northeastern US by applying relationships between surface PM2.5 and column AOD (calculated offline from speciated mass distributions) from a regional air quality model (CMAQ; 12×12 km2(exp) horizontal resolution). Seasonal average satellite-derived PM2.5 reveals more spatial detail and best captures observed surface PM2.5 levels during summer. At the daily scale, however, satellite-derived PM2.5 is not only subject to measurement uncertainties from satellite instruments, but more importantly to uncertainties in the relationship between surface PM2.5 and column AOD. Using 11 ground-based AOD measurements within 10 km of surface PM2.5 monitors, we show that uncertainties in modeled PM2.5∕AOD can explain more than 70 % of the spatial and temporal variance in the total uncertainty in daily satellite-derived PM2.5 evaluated at PM2.5 monitors. This finding implies that a successful geophysical approach to deriving daily PM2.5 from satellite AOD requires model skill at capturing day-to-day variations in PM2.5∕AOD relationships. Overall, we estimate that uncertainties in the modeled PM2.5∕AOD lead to an error of 11 µg m−3(exp) in daily satellite-derived PM2.5, and uncertainties in satellite AOD lead to an error of 8 µg m−3(exp). Using multi-platform ground, airborne, and radiosonde measurements, we show that uncertainties of modeled PM2.5∕AOD are mainly driven by model uncertainties in aerosol column mass and speciation, while model representation of relative humidity and aerosol vertical profile shape contributes some systematic biases. The parameterization of aerosol optical properties, which determines the mass extinction efficiency, also contributes to random uncertainty, with the size distribution being the largest source of uncertainty and hygroscopicity of inorganic salt the second largest. Future efforts to reduce uncertainty in geophysical approaches to derive surface PM2.5 from satellite AOD would thus benefit from improving model representation of aerosol vertical distribution and aerosol optical properties, to narrow uncertainty in satellite-derived PM2.5.

aerosol optical depth↗
DOE OSTI2022

WHONDRS Surface Water Chemistry and Organic Matter Characterization along the St. Lawrence River's Inland to Coastal Gradient, Eastern North America (v2)

This dataset supports a broader study examining the inland (Lake Ontario) to coastal (North Atlantic Ocean) geochemistry gradient along the St. Lawrence River in Canada and the United States. The St. Lawrence River is unique in that it contains the convergence of multiple water masses with distinct water signatures that mix only slightly as the river flows downstream. The dataset provides dissolved organic carbon (DOC) and organic matter characterization data generated from surface water. Samples were collected by researchers on board the Lampsilis research vessel (l’Université du Québec à Trois-Rivières) and small boats (St. Lawrence River Institute of Environmental Sciences, Cornwall) at 94 locations across and along the St. Lawrence River to capture longitudinal and transverse variation. Related data were collected and will be published separately in collaboration with the MicrEAU Laboratory (François Guillemette; l’Université du Québec à Trois-Rivières) and the Exploration of Coastal Hydrobiogeochemistry Across a Network of Gradients and Experiments (EXCHANGE) program. This dataset is comprised of one main data folder containing (1) file-level metadata; (2) data dictionary; (3) field metadata; (4) dissolved organic carbon (DOC, measured as non-purgeable organic carbon, NPOC) data (5) surface water sampling protocol; (6) readme; (7) methods codes; (8) international geo-sample number (IGSN) mapping file; and (9) folder of high resolution characterization of organic matter via 12 Tesla Fourier transform ion cyclotron resonance mass spectrometry (FTICR-MS) through the Environmental Molecular Sciences Laboratory (EMSL; https://www.pnnl.gov/environmental-molecular-sciences-laboratory). The FTICR folder contains two subfolders, one containing the .xml data files and the other containing instructions for using Formularity (https://omics.pnl.gov/software/formularity) and an R script to process the data based on the user's specific needs. All files are .csv, .pdf, .R, .ref, or .xml. The data package was originally published in November 202. It was updated in April 2025 (v2; modified files). See the change history section in the readme for details.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1973

A first look at atmospheric dynamics and temperature variations on Titan.

Pollack (1973) has used a radiative equilibrium model to match radiometric data for Titan and infers the atmospheric mass, composition, opacity, and gross vertical thermal structure. These results are used to estimate the atmospheric temperature variations by means of scaling analysis, taking into account dynamics both for a baroclinic wave regime and for an axially symmetric circulation regime. Horizontal temperature variations of the atmosphere and surface are found to be very small, and the circulation is found to be weak and probably axially symmetric. The small temperature variations appear to preclude the storage of volatiles in polar caps, so that the present atmospheric methane content may be due to a balance between outgassing and photodissociation.

Leovy, C. B.↗
NASA NTRS1981

Surface photometry of spiral galaxies. I - Theoretical color variation and surface brightness across spiral arms

A method for the construction of color variations and surface brightness across a galactic spiral arm, using the concepts and results of the density-wave theory, is developed. The color variation measured by the index Q varies with the ratio of the light of young stars to that of the old stars, the initial mass function, and the pattern speed. The fact that the variation of brightness and color across a spiral arm is relatively smooth and symmetrical and not sharp as in the galactic shock picture is due to the smearing-out effects of star formation inside a cloud in the regions behind the galactic shock, and the tendency of new stars to fall back to small galactocentric distances after their formation because of their lower angular velocity at birth. Both effects arise naturally in the context of the density-wave theory.

Yuan, C.↗
NASA NTRS1986

Satellite microwave and in situ observations of the Weddell Sea ice cover and its marginal ice zone

The radiative and physical characteristics of the Weddell Sea ice cover and its marginal ice zone are analyzed using multichannel satellite passive microwave data and ship and helicopter observations obtained during the 1983 Antarctic Marine Ecosystem Research. Winter and spring brightness temperatures are examined; spatial variability in the brightness temperatures of consolidated ice in winter and spring cyclic increases and decrease in brightness temperatures of consolidated ice with an amplitude of 50 K at 37 GHz and 20 K at 18 GHz are observed. The roles of variations in air temperature and surface characteristics in the variability of spring brightness temperatures are investigated. Ice concentrations are derived using the frequency and polarization techniques, and the data are compared with the helicopter and ship observations. Temporal changes in the ice margin structure and the mass balance of fresh water and of biological features of the marginal ice zone are studied.

Comiso, J. C.↗
NASA NTRS1993

The influence of thermal inertia on Mars' seasonal pressure variation and the effect of the weather component

Using a Leighton-Murray type diurnal and seasonal Mars thermal model, we found that it is possible to reproduce the seasonal variation in daily-averaged pressures (approximately 680-890 Pa) measured by Viking Lander 1 (VL1), during years without global dust storms, with a standard deviation of less than 5 Pa. In this simple model, surface CO2, frost condensation, and sublimation rates at each latitude are determined by the net effects of radiation, latent heat, and heat conduction in subsurface soil layers. An inherent assumption of our model is that the seasonal pressure variation is due entirely to the exchange of mass between the atmosphere and polar caps. However, the results of recent Mars GCM modeling have made it clear that there is a significant dynamical contribution to the seasonal pressure variation. This 'weather' component is primarily due to large-scale changes in atmospheric circulation, and its magnitude depends somewhat on the dust content of the atmosphere. The overall form of the theoretical weather component at the location of VL1, as calculated by the AMES GCM, remains the same over the typical range of Mars dust opacities.

Wood, S. E.↗
NASA NTRS2016

Lunar Phase Function at 1064 Nm from Lunar Orbiter Laser Altimeter Passive and Active Radiometry

We present initial calibration and results of passive radiometry collected by the Lunar Orbiter Laser Altimeter onboard the Lunar Reconnaissance Orbiter over the course of 12 months. After correcting for time- and temperature-dependent dark noise and detector responsivity variations, the LOLA passive radiometry measurements are brought onto the absolute radiance scale of the SELENE Spectral Profiler. The resulting photometric precision is estimated to be 5%. We leverage the unique ability of LOLA to measure normal albedo to explore the 1064 nm phase function's dependence on various geologic parameters. On a global scale, we find that iron abundance and optical maturity (quantified by FeO and OMAT) are the dominant controlling parameters. Titanium abundance (TiO2), surface roughness on decimeter to decameter scales, and soil thermo- physical properties have a smaller effect, but the latter two are correlated with OMAT, indicating that exposure age is the driving force behind their effects in a globally-averaged sense. The phase function also exhibits a dependence on surface slope at approximately 300 m baselines, possibly the result of mass wasting exposing immature material and/or less space weathering due to reduced sky visibility. Modeling the photometric function in the Hapke framework, we find that, relative to the highlands, the maria exhibit decreased backscattering, a smaller opposition effect (OE) width, and a smaller OE amplitude. Immature highlands regolith has a higher backscattering fraction and a larger OE width compared to mature highlands regolith. Within the maria, the backscattering fraction and OE width show little dependence on TiO2 and OMAT. Variations in the phase function shape at large phase angles are observed in and around the Copernican-aged Jackson crater, including its dark halo, a putative impact melt deposit. Finally, the phase function of the Reiner Gamma Formation behaves more optically immature than is typical for its composition and OMAT, suggesting the visible-to-near-infrared spectrum and phase function respond differently to the unusual regolith evolution and properties at this location. The phase difference map revealed additional geologically-influenced variations in the phase function's shape. In particular, variations were observed associated with the dark halo around Jackson crater, the impact ejecta of Copernicus and Giordano Bruno, and the Reiner Gamma Formation. For the latter, we 915 found that the phase function behaves more optically immature than the global phase function for its composition and OMAT, suggesting a difference in how the visible-to-near-IR spectrum and phase function respond to the unusual regolith evolution and properties at this location.

moon↗
NASA NTRS1989

Preliminary design of a large tetrahedral truss/hexagonal heatshield panel aerobrake

An aerobrake structural concept is introduced which consists of two primary components: (1) a lightweight erectable tetrahedral support truss; and (2) sandwich hexagonal heatshield panels which, when attached to the truss, form a continuous impermeable aerobraking surface. Generic finite element models and a general analysis procedure to design tetrahedral truss/hexagonal heatshield panel aerobrakes is developed, and values of the aerobrake design parameters which minimize mass and packaging volume for a 120-foot-diameter aerobrake are determined. Sensitivity of the aerobrake design to variations in design parameters is also assessed. The results show that a 120-foot-diameter aerobrake is viable using the concept presented (i.e., the aerobrake mass is less than or equal to 15 percent of the payload spacecraft mass). Minimizing the aerobrake mass (by increasing the number of rings in the support truss) however, leads to aerobrakes with the highest part count.

Dorsey, John T.↗
NASA NTRS2014

Development of a Linear Ion Trap Mass Spectrometer (LITMS) Investigation for Future Planetary Surface Missions

Future surface missions to Mars and other planetary bodies will benefit from continued advances in miniature sensor and sample handling technologies that enable high-performance chemical analyses of natural samples. Fine-scale (approx.1 mm and below) analyses of rock surfaces and interiors, such as exposed on a drill core, will permit (1) the detection of habitability markers including complex organics in association with their original depositional environment, and (2) the characterization of successive layers and gradients that can reveal the time-evolution of those environments. In particular, if broad-based and highly-sensitive mass spectrometry techniques could be brought to such scales, the resulting planetary science capability would be truly powerful. The Linear Ion Trap Mass Spectrometer (LITMS) investigation is designed to conduct fine-scale organic and inorganic analyses of short (approx.5-10 cm) rock cores such as could be acquired by a planetary lander or rover arm-based drill. LITMS combines both pyrolysis/gas chromatograph mass spectrometry (GCMS) of sub-sampled core fines, and laser desorption mass spectrometry (LDMS) of the intact core surface, using a common mass analyzer, enhanced from the design used in the Mars Organic Molecule Analyzer (MOMA) instrument on the 2018 ExoMars rover. LITMS additionally features developments based on the Sample Analysis at Mars (SAM) investigation on MSL and recent NASA-funded prototype efforts in laser mass spectrometry, pyrolysis, and precision subsampling. LITMS brings these combined capabilities to achieve its four measurement objectives: (1) Organics: Broad Survey Detect organic molecules over a wide range of molecular weight, volatility, electronegativity, concentration, and host mineralogy. (2) Organic: Molecular Structure Characterize internal molecular structure to identify individual compounds, and reveal functionalization and processing. (3) Inorganic Host Environment Assess the local chemical/mineralogical makeup of organic host phases to help determine deposition and preservation factors. (4) Chemical Stratigraphy Analyze the fine spatial distribution and variation of key species with depth.

Mars↗
NASA NTRS2000

Feeling Gravity's Pull: Gravity Modeling. The Gravity Field of Mars

Most people take the constant presence of gravitys pull for granted. However, the Earth's gravitational strength actually varies from location to location. This variation occurs because mass, which influences an object's gravitational pull, is not evenly distributed within the planet. Changes in topography, such as glacial movement, an earthquake, or a rise in the ocean level, can subtly affect the gravity field. An accurate measurement of the Earth's gravity field helps us understand the distribution of mass beneath the surface. This insight can assist us in locating petroleum, mineral deposits, ground water, and other valuable substances. Gravity mapping can also help notice or verify changes in sea surface height and other ocean characteristics. Such changes may indicate climate change from polar ice melting and other phenomena. In addition, gravity mapping can indicate how land moves under the surface after earthquakes and other plate tectonic processes. Finally, changes in the Earth's gravity field might indicate a shift in water distribution that could affect agriculture, water supplies for population centers, and long-term weather prediction. Scientists can map out the Earth's gravity field by watching satellite orbits. When a satellite shifts in vertical position, it might be passing over an area where gravity changes in strength. Gravity is only one factor that may shape a satellite's orbital path. To derive a gravity measurement from satellite movement, scientists must remove other factors that might affect a satellite's position: 1. Drag from atmospheric friction. 2. Pressure from solar radiation as it heads toward Earth and. as it is reflected off the surface of the Earth 3. Gravitational pull from the Sun, the Moon, and other planets in the Solar System. 4. The effect of tides. 5. Relativistic effects. Scientists must also correct for the satellite tracking process. For example, the tracking signal must be corrected for refraction through the atmosphere of the Earth. Supercomputers can calculate the effect of gravity for specific locations in space following a mathematical process known as spherical harmonics, which quantifies the gravity field of a planetary body. The process is based on Laplace's fundamental differential equation of gravity. The accuracy of a spherical harmonic solution is rated by its degree and order. Minute variations in gravity are measured against the geoid, a surface of constant gravity acceleration at mean sea level. The geoid reference gravity model strength includes the central body gravitational attraction (9.8 m/sq s) and a geopotential variation in latitude partially caused by the rotation of the Earth. The rotational effect modifies the shape of the geoid to be more like an ellipsoid, rather than a perfect, circle. Variations of gravity strength from the ellipsoidal reference model are measured in units called milli-Galileos (mGals). One mGal equals 10(exp -5) m/sq s. Research projects have also measured the gravity fields of other planetary bodies, as noted in the user profile that follows. From this information, we may make inferences about our own planet's internal structure and evolution. Moreover, mapping the gravity fields of other planets can help scientists plot the most fuel-efficient course for spacecraft expeditions to those planets.

Lemoine, Frank↗
NASA NTRS1982

Crater densities and geological histories of Rhea, Dione, Mimas and Tethys

The Saturn icy satellites Rhea, Dione, Mimas and Tethys were imaged by Voyager 1 with sufficiently high resolution to permit crater counts that will determine the character and magnitude of crater density variations and the relative ages of the different terrain types. The data show that the surfaces of these satellites are not uniformly cratered, and therefore not of uniform age. Density variations are found to occur at all diameters, especially in the distribution of large-diameter craters (greater than 50 km). It is concluded that this satellite group, despite their small size and mass and low density, have evolved beyond their original, accretionary state. The alteration processes appear in light of crater density determinations to have been active beyond the period of heavy bombardment. The energy source for the evolution of these bodies has yet to be determined.

Plescia, J. B.↗
NASA NTRS1992

Two-dimensional viscous accretion disk models. I - On meridional circulations in radiative regions

It has been suggested that thin Keplerian disks, like rotating stars, cannot attain hydrostatic and thermal equilibrium simultaneously if the angular velocity is a function of radius only. To address the problem of vertical structure and stability of thin accretion disks, we performed 2D axisymmetric calculations by solving the set of fully nonlinear hydrodynamic equations, including radiation transport. The results indicate that the variation of the rotational velocity with height is too small to sustain hydrostatic equilibrium. Instead, the variation of the radial (inflow) velocity with height suffices to establish a stationary state. In particular, we find that for low values of the viscosity parameter alpha, there is mass outflow in the central parts of the disk close to the equatorial plane, and inflow only near the surface. Only for modestly high values of alpha, the flow throughout the disk is directed inward. Thus, the flow within the disk is the result of viscous radial inflow plus a circulatory flow directed outward in the midplane and inward near the surface.

Kley, W.↗
DOE OSTI2024

DRIFT: Diffusivity Regulation of Isotopes by nanoconFinementT

Spontaneous isotope fractionation has been reported under nanoconfinement conditions in naturally occurring systems, but the origin of this phenomena is currently unknown. Two existing hypotheses have been proposed, one based on changes in the solvation environment of the isotopes that reduces the non-mass dependent hydrodynamics contribution to diffusion. The other is that isotopes have mass-dependent surface adsorption, varying their total diffusion through nanoconfined channels. To investigate these hypotheses, benchtop experiments, nuclear magnetic resonance (NMR) spectroscopy, and molecule scale modeling were applied. Classical molecular dynamics simulations identified that the Na + and Cl - hydration shells across the three different salt solutions ( 22 Na 35 Cl, 23 Na 35 Cl, 24 Na 35 Cl) did not vary as a function of the Na + isotope, but that there was a significant pore size effect, with larger hydration shells at larger pore sizes. Additionally, while total adsorption times did not vary as a function of the Na + isotope or pore size, the free ion concentration, or those adsorbed on the surface for <5% of the simulation time did exhibit isotope dependence. Experimentally, challenges occurred developing a repeatable experiment, but NMR characterization of water diffusion rates through ordered alumina membranes was able to identify the existence of two distinct water environments associated with water inside and outside the pore. Further NMR studies could be used to confirm variation in hydration shells and diffusion rates of dissolved ions in water. Ultimately, mass-dependence adsorption is a primary driver of variations in isotope diffusion rates, rather than variation in hydration shells that occur under nanoconfinement.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1998

An Integrated Tool for the Coupled Thermal and Mechanical Analysis of Pyrolyzing Heatshield Materials

Materials that pyrolyze at elevated temperature have been commonly used as thermal protection materials in hypersonic flight, and advanced pyrolyzing materials for this purpose continue to be developed. Because of the large temperature gradients that can arise in thermal protection materials, significant thermal stresses can develop. Advanced applications of pyrolytic materials are calling for more complex heatshield configurations, making accurate thermal stress analysis more important, and more challenging. For non-pyrolyzing materials, many finite element codes are available and capable of performing coupled thermal-mechanical analyses. These codes do not, however, have a built-in capability to perform analyses that include pyrolysis effects. When a pyrolyzing material is heated, one or more components of the original virgin material pyrolyze and create a gas. This gas flows away from the pyrolysis zone to the surface, resulting in a reduction in surface heating. A porous residue, referred to as char, remains in place of the virgin material. While the processes involved can be complex, it has been found that a simple physical model in which virgin material reacts to form char and pyrolysis gas, will yield satisfactory analytical results. Specifically, the effects that must be modeled include: (1) Variation of thermal properties (density, specific heat, thermal conductivity) as the material composition changes; (2) Energy released or absorbed by the pyrolysis reactions; (3) Energy convected by the flow of pyrolysis gas from the interior to the surface; (4) The reduction in surface heating due to surface blowing; and (5) Chemical and mass diffusion effects at the surface between the pyrolysis gas and edge gas Computational tools for the one-dimensional thermal analysis these materials exist and have proven to be reliable design tools. The objective of the present work is to extend the analysis capabilities of pyrolyzing materials to axisymmetric configurations, and to couple thermal and mechanical analyses so that thermal stresses may be efficiently and accurately calculated.

Pronchick, Stephen W.↗
NASA NTRS2012

Modeling of Melt Growth During Carbothermal Processing of Lunar Regolith

The carbothermal processing of lunar regolith has been proposed as a means to produce carbon monoxide and ultimately oxygen to support human exploration of the moon. In this process, gaseous methane is pyrolyzed as it flows over the hot surface of a molten zone of lunar regolith and is converted to carbon and hydrogen. Carbon gets deposited on the surface of the melt, and mixes and reacts with the metal oxides in it to produce carbon monoxide that bubbles out of the melt. Carbon monoxide is further processed in other reactors downstream to ultimately produce oxygen. The amount of oxygen produced crucially depends on the amount of regolith that is molten. In this paper we develop a model of the heat transfer in carbothermal processing. Regolith in a suitable container is heated by a heat flux at its surface such as by continuously shining a beam of solar energy or a laser on it. The regolith on the surface absorbs the energy and its temperature rises until it attains the melting point. The energy from the heat flux is then used for the latent heat necessary to change phase from solid to liquid, after which the temperature continues to rise. Thus a small melt pool appears under the heated zone shortly after the heat flux is turned on. As time progresses, the pool absorbs more heat and supplies the energy required to melt more of the regolith, and the size of the molten zone increases. Ultimately, a steady-state is achieved when the heat flux absorbed by the melt is balanced by radiative losses from the surface. In this paper, we model the melting and the growth of the melt zone with time in a bed of regolith when a portion of its surface is subjected to a constant heat flux. The heat flux is assumed to impinge on a circular area. Our model is based on an axisymmetric three-dimensional variation of the temperature field in the domain. Heat transfer occurs only by conduction, and effects of convective heat transport are assumed negligible. Radiative heat loss from the surface of the melt and the regolith to the surroundings is permitted. We perform numerical computations to determine the shape and the mass of the melt at steady state and its time evolution. We first neglect the volume change upon melting, and subsequently perform calculations including it. Predictions from our model are compared to test data to determine the effective thermal conductivities of the regolith and the melt that are compatible with the data

Balasubramaniam, R.↗
DOE OSTI2026

A stress-based fracture model for reacting metal ejecta

The evolution of reacting metal ejecta continues to be a topic of interest at the forefront of metals in reactive and extreme environments. Ejecta are small particles formed when the surface of a metal undergoes Richtmyer–Meshkov instability from a strong shock. Experiments have shown that in the case where ejecta are in ambient conditions that induce a reaction, the ejecta behave irregularly. The ejecta temperature rises and then plateaus, and the acceleration profile shows unexpected jumps. These variations are assumed to be related to the exothermic heat release and particle mass loss caused by the reaction. To explain this phenomenon, efforts to model this in simulations have increased. While current models can capture many of these physical processes, they currently assign a constant reaction shell thickness with little physical reasoning. This work remedies this problem by assigning a dynamic physically informed shell thickness to the reacting particles, using solid analysis. The shell thickness of the particles impacts the rate of change of reacted mass in the system, as well as the rate at which the particles react. The model is based on a simple stress–strain relationship and gives a dynamic assignment for when the reacting particle should begin to fracture. We compare our model to the previous computational and simulation data to analyze the effects of different model parameters.

42 ENGINEERING↗