Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Slow dynamics models”

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.

1,197 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS2026

The Significance of Water Vapor Isotopes in Improving Weather Prediction

Water vapor isotopes carry the integrated history of evaporation, condensation, mixing, and transport. Although previous studies have shown potential to improve forecasts under controlled conditions, real-world applications have been limited by systematic biases in both models and satellite retrievals arising from sparse measurements in the free troposphere. Here we assimilate mid‑tropospheric δD retrievals (peak sensitivity ~4.2 km) from the Infrared Atmospheric Sounding Interferometer into the Isotope‑incorporated Global Spectral Model and evaluate the added value beyond co‑assimilated temperature and specific humidity with identical spatial and temporal coverage. Assimilating δD improves 0–120 h forecasts of wind, temperature, specific humidity, and geopotential height, with the largest gains in the midlatitudes; heavy‑precipitation skill also increases for thresholds >3 mm per 6 h. Demonstrated in a coarse‑resolution configuration with limited observations, the results indicate that isotopic information strengthens transport tracking and hydrological constraints, motivating evaluation in operational high‑resolution forecasting systems.

Hydrology↗
NASA NTRS2026

Organic Inventory - Planetary Protection on the Moon

All spacecraft generate and carry contaminants, i.e., unwanted and potentially harmful material. When a spacecraft lands and operates in vacuum, as onto Earth’s Moon, it introduces contaminants into its environment that may compromise mission science objectives and engineering performance. Contamination may degrade sites of unique value to planetary science or in situ resource utilization. This presentation will identify and compare source terms and transport vectors for contaminants – in particular, organic material – generated by landed spacecraft. An integrated modeling framework for the organic contamination footprint of spacecraft missions will be described and presented.

Planetary Protection↗
NASA NTRS2026

Beyond Melting: Amorphous Bonding for Joining and Consolidation

Crystallization may be the hidden constraint in thermoplastic composite manufacturing. It requires tightly controlled cooling, induces residual stresses through shrinkage, and introduces path-dependent behavior that complicates predictive modeling yet remains essential for structural performance. This work asks: can bonding be achieved without relying on melt-driven crystallization? To address this, thin (5–20 μm) polyetherimide (PEI) interlayers are pre-healed to slow-cooled polyaryletherketone (PAEK) in two contexts. The first, Thermabond®, is sub-melt joining of low melt-PAEK laminates. Results show that bond quality is governed primarily by processing (i.e., adequate healing and film handling) rather than modest changes in interlayer thickness. This concept is then extended to laminate-scale manufacturing through an architecture known as OATMEAL (Out-of-autoclave Amorphous/semicrystalline Thermoplastic Material for Energy-efficient Aerospace-grade Laminates). PEI is healed to carbon fiber reinforced polyetheretherketone (PEEK) at the prepreg and excess PEI is then ablated from the surface. Crystallinity is developed off-line during prepreg fabrication, while subsequent consolidation occurs below the melt temperature to preserve it. Cross-ply warpage experiments show that, contrary to intuition, repeated amorphous interfaces reduce global curvature by lowering the effective stress lock-in temperature and eliminating crystallization shrinkage from the lamina response. Correspondingly, laminate behavior is accurately predicted using classical laminate theory (CLT) with a single effective stress-free temperature, whereas conventional CF/PEEK requires accounting for crystallization-driven effects. By decoupling interfacial healing from crystallization, OATMEAL enables sub-melt consolidation, reduces energy consumption by up to 75%, and increases manufacturing throughput by fivefold. These results demonstrate that amorphous bonding is not only a joining strategy, but a pathway to more predictable and scalable thermoplastic composite manufacturing.

solidification↗
DOE OSTI2021

Tradeoffs and Synergies in Tropical Forest Root Traits and Dynamics for Nutrient and Water Acquisition: Field and Modeling Advances

Vegetation processes are fundamentally limited by nutrient and water availability, the uptake of which is mediated by plant roots in terrestrial ecosystems. While tropical forests play a central role in global water, carbon, and nutrient cycling, we know very little about tradeoffs and synergies in root traits that respond to resource scarcity. Tropical trees face a unique set of resource limitations, with rock-derived nutrients and moisture seasonality governing many ecosystem functions, and nutrient versus water availability often separated spatially and temporally. Root traits that characterize biomass, depth distributions, production and phenology, morphology, physiology, chemistry, and symbiotic relationships can be predictive of plants’ capacities to access and acquire nutrients and water, with links to aboveground processes like transpiration, wood productivity, and leaf phenology. In this review, we identify an emerging trend in the literature that tropical fine root biomass and production in surface soils are greatest in infertile or sufficiently moist soils. We also identify interesting paradoxes in tropical forest root responses to changing resources that merit further exploration. For example, specific root length, which typically increases under resource scarcity to expand the volume of soil explored, instead can increase with greater base cation availability, both across natural tropical forest gradients and in fertilization experiments. Also, nutrient additions, rather than reducing mycorrhizal colonization of fine roots as might be expected, increased colonization rates under scenarios of water scarcity in some forests. Efforts to include fine root traits and functions in vegetation models have grown more sophisticated over time, yet there is a disconnect between the emphasis in models characterizing nutrient and water uptake rates and carbon costs versus the emphasis in field experiments on measuring root biomass, production, and morphology in response to changes in resource availability. Closer integration of field and modeling efforts could connect mechanistic investigation of fine-root dynamics to ecosystem-scale understanding of nutrient and water cycling, allowing us to better predict tropical forest-climate feedbacks.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Formation and Shape Changing of Conductive Helical Ribbons via Deposition of Highly Stressed Films on Mechanically Responsive Substrates

Abstract This work demonstrates that the electrodeposition of highly stressed films on compliant ribbons is a robust process to obtain helical structures with excellent mechanical stability and potentially high thermal and electrical conductance. Electrodeposition on end‐tethered ribbons alters their axial and bending stiffness while imparting mechanical stress to drive the formation of a helix with a microscale diameter and pitch in a controlled and scalable manner. The process generates helices with diameters and pitches between 80 and 200 µm and lengths as large as several millimeters. The approach is amenable to parallel processing a large number of 3D structures on any substrate, including large‐area semiconductor wafers. This phenomenon is explained in terms of the change of stress gradients as material is added. Applications of the fabricated helices include antennas, metamaterials, and slow‐wave structures in frequency ranges not previously attainable.

Chemistry↗
DOE OSTI2025

Rethinking 𝛼−RuCl 3 : Parameters, models, and phase diagram

RuCl 3 was likely the first ever deliberately synthesized ruthenium compound, following the discovery of the 44 Ru element in 1844. For a long time it was known as an oxidation catalyst, with its physical properties being discrepant and confusing, until a decade ago when its allotropic form 𝛼−RuCl 3 rose to exceptional prominence. This “rediscovery” of 𝛼−RuCl 3 has not only reshaped the hunt for a material manifestation of the Kitaev spin liquid, but it has opened the floodgates of theoretical and experimental research in the many unusual phases and excitations that the anisotropic-exchange magnets as a class of compounds have to offer. Given its importance for the field of Kitaev materials, it is astonishing that the low-energy spin model that describes this compound and its possible proximity to the much-desired spin-liquid state is still a subject of significant debate ten years later. In the present study, we argue that the existing key phenomenological observations put strong natural constraints on the effective microscopic spin model of 𝛼−RuCl 3 , and specifically on its spin-orbit-induced anisotropic-exchange parameters that are responsible for the nontrivial physical properties of this material. These constraints allow one to focus on the relevant region of the multidimensional phase diagram of the 𝛼−RuCl 3 model, suggest an intuitive description of it via a different parametrization of the exchange matrix, offer a unifying view on the earlier assessments of its parameters, and bring closer together several approaches to the derivation of anisotropic-exchange models. We explore extended phase diagrams relevant to the 𝛼−RuCl 3 parameter space using quasiclassical, Luttinger-Tisza, exact diagonalization, and density-matrix renormalization-group methods, demonstrating a remarkably close quantitative accord between them on the general structure and hierarchy of the phases, with the zigzag, ferromagnetic, and incommensurate phases that are proximate to each other. As a result, one of the highlights is the detailed agreement on the nature of the incommensurate phases that realize two distinct counterrotating helical states.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2026

CFD-Assisted Nodal Modeling of Sloshing in a Cryogenic Propellant Tank

During autogenous pressurization, tank sloshing causes a significant increase in pressurant consumption to maintain constant ullage pressure during draining of the tank. This increase in pressurant consumption is caused by a significant increase in condensation at the liquid-vapor interface. Sloshing strongly affects the liquid side heat transfer coefficient and thereby the condensation rate. Traditionally, sloshing is modeled by CFD code using the VOF (Volume of Fluid) method to track the liquid vapor interface during sloshing. CFD calculations require a very fine grid to accurately compute the heat and mass transfer at the interface. Therefore, computations are time consuming and prohibit performing many parametric studies often needed during the design of a new system. This paper describes an alternative approach by developing a CFD-assisted nodal model to predict system parameters more economically with reasonable accuracy. In this approach, a multi-node model of tank pressurization was developed using GFSSP. A multi-node model was needed to account for stratification. The model computes heat and mass transfer at the interface to calculate the condensation rate. The liquid side heat transfer is computed using the parameters of sloshing dynamics such as frequency, wave amplitude, and interface area. The parameters of sloshing dynamics are computed by the CFD code LOCI-Stream. The model predictions were compared with test data for several cases.

Cryogenic Tank Sloshing↗
DOE OSTI2026

Divergent Rickettsia species exhibit distinct mechanisms of actin-based motility

Many Rickettsia species undergo actin-based motility to promote cell-cell spread during infection. Rickettsial genomes often encode two motility effectors, RickA and Sca2, which in the spotted fever group I species, Rickettsia parkeri act by activating the host Arp2/3 complex and by mimicking eukaryotic formins, respectively. The function of RickA and Sca2 orthologs in the distantly related species Rickettsia bellii was unclear. We report that R. bellii RickA activates the host Arp2/3 complex but has no discernible role in bacterial motility. The R. bellii Sca2 ortholog, Sca2/6, nucleates and elongates actin with a flexible structure and an unusual actin monomer-binding motif in a mechanism distinct from formins or other microbial actin nucleators. R. bellii motility is solely correlated with Sca2/6 localization and, compared with R. parkeri motility, is slow and meandering, generating distinctly organized actin tails. The evolutionary flexibility in the mechanism and regulation of rickettsial actin-based motility suggests similar adaptability for other microbes.

Bacher, Meghan C↗
DOE OSTI2026

Vanishing dynamic strength measured in a transient high-pressure phase of tin

This work reports two independent experimental estimations of the dynamic strength of the body-centered tetragonal (bct) 𝛾 phase of tin, which occurs at pressures above about 9 GPa and transforms promptly back to the ambient 𝛽 phase upon pressure release. Measuring strength in such a transient high-pressure phase is challenging. One measurement used free-surface Richtmyer-Meshkov instabilities generated with gas gun impact. Another set of measurements used ramp-release loading in a pulsed-power facility. Strength estimations came from comparison to simulations using a comprehensive multiphase modeling framework that includes equation of state, shear moduli, and strength separately for each phase, and treats mixed-phase regions. Both experimental methods found the 𝛾-phase strength to be very low, within experimental uncertainty of zero. The existing literature on other materials, by contrast, almost universally reports higher strength in transient high-pressure phases compared with ambient phases. Recent Molecular-Dynamics simulations on tin in the literature showed almost zero deviatoric strength during a deformation-induced bct → bct transformation in which one of the 𝛾-phase 𝑎 axes “flips” to the 𝑐 axis. This reorientation currently provides the most plausible explanation for the low observed strength in 𝛾-phase tin.

36 MATERIALS SCIENCE↗
DOE OSTI2026

SODAs: sparse optimization for the discovery of differential and algebraic equations

Differential-algebraic equations (DAEs) integrate ordinary differential equations (ODEs) with algebraic constraints, providing a fundamental framework for developing models of dynamical systems characterized by time-scale separation, conservation laws and physical constraints. While sparse optimization has revolutionized model development by allowing data-driven discovery of parsimonious models from a library of possible equations, existing approaches for dynamical systems assume DAEs can be reduced to ODEs by eliminating variables before model discovery. This assumption limits the applicability of such methods for DAE systems with unknown constraints and time scales. We introduce sparse optimization for differential-algebraic systems (SODAs), a data-driven method for the identification of DAEs in their explicit form. By discovering the algebraic and dynamic components sequentially without prior identification of the algebraic variables, this approach leads to a sequence of convex optimization problems. It has the advantage of discovering interpretable models that preserve the structure of the underlying physical system. To this end, SODAs improves since SODAs is singular numerical stability when handling high correlations between library terms, caused by near-perfect algebraic relationships, by iteratively refining the conditioning of the candidate library. We demonstrate the performance of our method on biological, mechanical and electrical systems, showcasing its robustness to noise in both simulated time series and real-time experimental data.

DAE↗
DOE OSTI2026

NanoPSD: A software for automatic detection of Nano-Particle Shape Distribution in electron microscopy images

Accurate quantification of the size and morphology of nanoparticles from electron microscopy (EM) images is essential to understand growth mechanisms, surface reactivity, and functional behavior in nanoscale materials. Manual analysis remains slow, subjective, and difficult to reproduce in large datasets. We introduce NanoPSD (Nano-Particle Shape Distribution), an open-source and fully automated framework for quantitative particle detection and morphology analysis from EM images. NanoPSD integrates adaptive contrast enhancement, polarity-agnostic scale-bar detection, Optical Character Recognition (OCR)-based calibration, and classical segmentation via Otsu thresholding with morphological refinement. Particle contours are used to extract geometric descriptors, including equivalent circular diameter, aspect ratio, circularity, and solidity, enabling automated classification into spherical, rod-like, and aggregate morphologies. The framework supports both single-image and batch processing, generating publication-quality visualizations, LaTeX-ready tables, and structured comma-separated values (CSV) datasets. As a demonstration, we applied NanoPSD to plasma-synthesized nanoparticle samples diagnosed via transmission electron microscopy (TEM). The code produced statistically robust size and morphology distributions spanning a few to tens of nanometers with minimal user supervision. The pipeline demonstrates high reproducibility and scalability, processing large image collections with consistent calibration and output formatting. Its modular design enables seamless integration of future deep-learning-based segmentation models, providing a pathway toward intelligent, data-driven electron microscopy analysis.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Mitigating electrochemical degradation in CsPbBr{sub 3} gamma detectors by organic and inorganic encapsulation.

CsPbBr3 perovskite semiconductors have emerged as a leading candidate for nextgeneration radiation detectors because of their exceptional charge transport properties, defect tolerance, and record-breaking sensitivity and energy resolution. Their long-term stability, however, is hindered by electrode-driven electrochemical decomposition, which is accelerated by moisture- and oxygen-assisted ion migration during operation. Here, we investigated organic and inorganic encapsulation strategies as both environmental barriers and means to suppress interfacial degradation pathways. Atomic layer deposition (ALD) of Al2O3 provided a conformal passivation layer that blocked environmental ingress, suppressed ionic diffusion, reduced leakage current, enhanced energy resolution and expanded the operational electric-field window beyond 5 kV∙cm1 . By contrast, organic encapsulants such as paraffin wax and polystyrene slowed moisture diffusion but did not suppress interfacial reactions, with wax extending stability to over 90 days. These results show that ALD-Al2O3 suppresses dominant interfacial degradation pathways, enabling stable, high-field operation and advancing the practical deployment of CsPbBr3 γ-ray detectors.

Unal, Mustafa↗
NASA NTRSDate not supplied

Advances in Design Capabilities for Planetary Missions from the NASA Entry Systems Modeling and Instrumentation Portfolio

The Entry Systems Modeling project (ESM) is supported by both the NASA Space Technology and the Science Mission Directorates and focuses on developing simulation tools and validated models for characterizing the performance of entry systems tailored to planetary destinations across the Solar System. ESM is organized into six technical capability areas that together address all relevant factors related to spacecraft entry, as well as some aspects of descent: Thermal Protection System (TPS) Materials; Aerothermodynamics; Entry & Descent Vehicle Dynamics; Guidance, Navigation, and Control; Vehicle Systems Analysis; and Advanced Tools and Numerical Methods. Development within the capability areas is undertaken explicitly with a focus on transition and infusion to science missions, human exploration missions, and commercial space activities. The present talk details developments that specifically impact science missions, including simulation tool capabilities that aid in mission design and model development to understand entry system performance at a given destination. Examples of the successful infusion and transition of such project outcomes to science missions also are provided. Several simulation tool development efforts within ESM have resulted in new design capabilities for missions. One such outcome is improved toolsets for mission trajectory and concept of operations design. Specifically, an initiative to couple a leading tool for entry, ascent/descent, and orbital trajectory optimization (Program to Optimize Simulated Trajectories II or POST2) to those used within the Agency for interplanetary trajectory optimization (Copernicus and Monte) has made substantial progress, with the outcomes to date promising to allow efficient trajectory optimization across mission phases. Additionally, toolchains for the evaluation of vehicle performance during entry and descent have been developed that allow assessment of multi-dimensional aeroheating on detailed vehicle geometries, characterization of deployment and inflation of parachutes, and assessment of vehicle dynamic stability during descent. These capabilities are achieved by coupling diverse sets of physics together – material response, computational fluid dynamics, radiation, and vehicle dynamics – to suitably describe complex entry and descent phenomena. Several model development and validation efforts for specific destinations and entry regimes also are underway within the ESM project. For instance, new experimental capabilities to validate radiation models at low densities/high altitudes recently have been established with project support, specifically the Low-Density Shock Tube (LDST) at the NASA Ames Research Center Electric Arc Shock Tube (EAST) facility. The LDST is being leveraged to develop improved models of shock layer kinetics and radiation in Titan atmospheres, while future studies will be conducted in the LDST and the existing high velocity shock tube to provide validation data for radiation models of Venus, Ice Giants, and Mars atmospheres. Models describing the aerothermal and thermo-structural performance of Thermal Protection System (TPS) materials has been another focus, with multiscale modeling activities on-going for the two leading TPS materials applicable to a range of entry conditions and science missions: the Phenolic-Impregnated Carbon Ablator (PICA) and woven materials like 3D Mid-Density Carbon Phenolic (3MDCP). A continual effort is made to infuse and transition outcomes from ESM simulation tool and model development activities into relevant science missions. Significant progress has been made on this front, with missions such as Dragonfly, DAVINCI, and Mars Missions benefitting from project outcomes. The groundwork also is being laid to provide insights into forward looking missions to Gas/Ice Giants as well as for potential sample returns.

Justin Haskins↗
NASA NTRS2014

Analytical Model for Steady Flow through a Finite Channel with One Porous Wall with Arbitrary Variable Suction or Injection

This paper presents an exact solution of two-dimensional laminar flow through a finite length channel with one porous wall. It improves upon previous solutions by (1) satisfying the no-slip boundary condition at the channel dead end, (2) adding a turbulent term to the porous wall boundary condition, (3) allowing for arbitrary variable suction or injection across the porous wall, and (4) model validation against new cryogenic liquid hydrogen and oxygen experimental data. Of particular interest in the current work is the modeling of cryogenic propellant flow through a porous liquid acquisition device (LAD) screen and channel inside a propellant tank. First, a detailed review of the literature is presented for previously attempted solutions to channel flow with one porous wall. Next, the governing equations, boundary conditions, and model assumptions are used to derive the analytical flow solution and present general model results for pressure and velocity fields within the channel. Then, the model solution is compared with horizontal LAD channel flow data in liquid oxygen as well as vertical LAD channel flow data in an inverted outflow configuration in liquid hydrogen. Model results are used to update the static cryogenic bubble point pressure model with a dynamic bubble point term which factors in enhanced convection and cooling at the screen during propellant outflow. Convective heat transfer at the LAD screen during outflow is also quantified by comparing model and data. The new analytical flow solution with the dynamic bubble point model is shown to compare well with available cryogenic experimental data

Navier Stokes Equations↗
DOE OSTI2026

Adaptation of virtual synchronous generators to dynamic conditions in power grids

Virtual synchronous generators (VSGs) are widely adopted as grid-forming controls for inverter-based resources. However, when grid conditions vary significantly as characterized by changes in short-circuit ratio (SCR) and the reactance-to-resistance (X/R) ratio, fixed-gain designs and the commonly used P–Q decoupling assumption can become inaccurate. Such conditions can degrade transient power performance, leading to oscillations, prolonged settling, and overshoot, particularly in stiff-grid operating points. This paper quantifies how grid strength and impedance-dependent coupling affect the active–reactive power dynamics of a conventional VSG over a broad range of SCR and X/R values. An adaptive VSG tuning framework is then developed by combining (i) a coupling-explicit, impedance-parameterized state-space model to enable systematic controller synthesis, (ii) a full-state-feedback law designed via pole placement to meet prescribed damping and settling-time specifications, and (iii) a physics-informed neural network (PINN)–based online grid-impedance estimator that updates controller gains in real time as grid conditions vary. Offline simulations in MATLAB/Simulink and real-time validation on an OPAL-RT platform show that the proposed method preserves consistent damping and settling behavior with reduced overshoot across wide SCR and X/R ranges, compared with fixed-gain VSG baselines.

Adaptive control↗
DOE OSTI2025

Generative learning of densities on manifolds

A generative modeling framework is proposed that combines diffusion models and manifold learning to efficiently sample data densities on manifolds. The approach utilizes Diffusion Maps to uncover possible low-dimensional underlying (latent) spaces in the high-dimensional data (ambient) space. Two approaches for sampling from the latent data density are described. The first is a score-based diffusion model, which is trained to map a standard normal distribution to the latent data distribution using a neural network. The second one involves solving an Itô stochastic differential equation in the latent space. Additional realizations of the data are generated by lifting the samples back to the ambient space using Double Diffusion Maps , a recently introduced technique typically employed in studying dynamical system reduction; here the focus lies in sampling densities rather than system dynamics. The proposed approaches enable sampling high dimensional data densities restricted to low-dimensional, a priori unknown manifolds. The efficacy of the proposed framework is demonstrated through a benchmark problem and a material with multiscale structure.

Double diffusion maps↗
DOE OSTI2026

Future Projection of Tropical Upper-Tropospheric Troughs and Implications for Tropical Cyclone Activity

Summertime tropical upper-tropospheric troughs (TUTTs) provide a unified framework to better understand how extratropical and tropical forcings jointly modulate basin-scale tropical cyclone (TC) activity. In this study, we examine future changes in TUTTs and their implications for TC activity. Multimodel ensemble-mean projections from 45 Coupled Model Intercomparison Project phase 6 (CMIP6) models suggest a contraction of the Pacific TUTT and an expansion of the Atlantic TUTT as the climate warms. Consistently, future changes in environment-based TC indices indicate that the large-scale conditions will become more favorable for TC genesis and intensification over the central North Pacific but less favorable over the tropical North Atlantic and Gulf of Mexico. Utilizing a TC-permitting large-ensemble dataset [i.e., the Database for Policy Decision-Making for Future Climate Change (d4PDF)] that adequately captures the observed interannual TUTT–TC relationships, we further confirm the impacts of projected TUTT changes on the TC activity in a warmer climate. In contrast, the TUTT–TC relationship is poorly represented in most CMIP6 High-Resolution Model Intercomparison Project (HighResMIP) models; such deficiencies call for caution when assessing future TC risk based on explicitly tracked TCs in these models. Additionally, CMIP6 projections show large intermodel spread in TUTT changes, implying uncertainty in projected TC activity, especially over the central-to-eastern Pacific and the North Atlantic. This intermodel spread is associated with interhemispheric sea surface temperature warming asymmetry, which leads to a meridional shift of the intertropical convergence zone (ITCZ) and the simultaneous weakening or strengthening of TUTTs in the North Pacific and North Atlantic. The potential contributions of anthropogenic aerosol forcing and oceanic circulation to this interhemispheric warming asymmetry are briefly examined.

Climate Change↗
NASA NTRS2004

Ejecta Generation and Redistribution on 433 Eros: Modeling Ejecta Launch Conditions

The NEAR-Shoemaker mission to asteroid 433 Eros presents an unprecedented opportunity to gain fundamental new knowledge about the processes governing regolith formation and redistribution on small bodies. NEAR-Shoemaker’s high-resolution imaging of the surface of Eros makes the asteroid a valuable and heretofore unparalleled laboratory for the detailed study of impact ejecta reaccretion and regolith redistribution on low-gravity (of order 10 -3 g) objects. Regolith is produced on asteroids by impact cratering, and the existence of regolith on the smallest solar system bodies supports the view that some of the ejecta from impact events on such objects may be retained. Impact craters and retained ejecta on low-gravity objects like Eros represent valuable natural laboratories for evaluating various models of impact cratering processes, since they may present crater structures or ejecta features that either do not form or are hidden on higher-gravity bodies like the Moon. Further, quantifying the extent to which impact processes generate and redistribute regoliths on small body surfaces (excavation depths, retained fraction, turnover timescales, etc.) is pivotal to the issue of how to relate meteoritical samples to their asteroidal parent bodies when surficial processes ( i.e., “space weathering”) may disguise or cover up underlying material and confound the ability of remote sensing techniques to provide reliable mineralogical assays of the parent objects. The rich variety of data on Eros’ regolith properties and distribution returned by NEAR-Shoemaker now require detailed analysis in order to take full advantage of the clues these observations offer for elucidating details of the impact cratering process on small bodies. Complicating simple interpretations of crater and ejecta morphology are dynamical effects on ejecta emplacement resulting from Eros’ irregular shape, rapid (5.27 hr) rotation, and low gravity. Figure 1 shows the very different ejecta deposit morphology that can result if the effects of rotation alone are neglected. Considering the additional complicating factors of Eros’ irregular shape and complex gravitational field, simple calculations of the extent and thickness of ejecta blankets and the spatial distribution of ejecta blocks from basic crater scaling laws or numerical hydrocodes alone do not suffice. In order to fully interpret the suite of NEAR-Shoemaker observations of regolith features across the surface of Eros and to evaluate various impact models for specific craters on the asteroid, detailed dynamical modeling of the deposition of crater ejecta from those craters is required . Here, I describe some modifications and improvements to the dynamical model being used for these studies.

D D Durda↗