Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “integrated modeling approach”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Modeling and Calibration of Supplier Selection Problem in Freight Agent-Based Simulations

Freight transportation modeling often struggles with data limitations, especially in accurately representing complex supplier selection processes and their impact on network flows. This research addresses this critical gap by developing a large-scale, calibrated agent-based model for supplier selection, complemented by a probabilistic heuristic for international shipments. Our approach integrates trade relationships between industry sectors, transportation costs, and a supplier-rating model adapted from existing literature. The model’s core objective is to minimize the discrepancy between modeled and observed commodity flows while ensuring a close match to regional shipping distance distributions. Implemented and tested across four major U.S. metropolitan areas—Atlanta, Chicago, Dallas–Fort Worth, and Los Angeles—the model demonstrates high fidelity in replicating observed freight patterns. Key findings reveal consistent alignment with national shipping distance trends and highlight significant spatial variations in commodity trade assignments and demand across the study regions. This behaviorally informed and transport-sensitive framework is designed to approximate real-world decision making, providing a robust tool for policymakers and planners to evaluate targeted interventions, assess infrastructure investments, and enhance supply chain resilience in the face of disruptions.

Ismael, Abdelrahman (ORCID:0000000303712110)↗
DOE OSTI2023

Autonomous Inversion of In Situ Deformation Measurement Data for Injection-Induced Stress Change

Geologic carbon storage (GCS) is likely to play a key part of the global effort to dramatically reduce CO2 emissions and perhaps even reduce atmospheric CO2 concentrations through carbon negative operations. A critical part of effort to commercialize and widely deploy this technology is developing the capability to rapidly assimilate real-time monitoring data into a form that will enable site operators to make decisions to manage the safe and efficient operations. Two of the risks associate with GCS are the risk of inducing fractures in the sealing formations that can create leakage pathways and the risk of inducing earthquakes of sufficient magnitude to cause public concern, property damage, or safety risks. To properly manage these risks the site operator needs to know the initial state of stress, the change in stress induced by injection, and the relationship between operational parameters such as injection rate and pressure and the change in stress. Current methods of estimating the change in stress require choosing the type of constitutive model and the model parameters based on core, log, and geophysical data during the characterization phase, with little feedback from operational observations to validate or refine these choices. These characterization methods interrogate the geologic formations using length scales, loading rates or magnitudes that are quite different from those encountered by the actual storage system. It is shown that errors in the assumed constitutive response, even when informed by laboratory tests on core samples, are likely to be common, large, and underestimate the magnitude of stress change caused by injection. Recent advances in borehole-based strain instruments and borehole and surface-based tilt and displacement instruments have now enabled monitoring of the deformation of the storage system throughout its operational lifespan. This data can enable validation and refinement of the knowledge of the geomechanical properties and state of the system, but brings with it a challenge to transform the raw data into actionable knowledge. We demonstrate a method that uses automatic differentiation and a finite-element based geomechanical model perform a gradient-based deterministic inversion of geomechanical monitoring data. This approach allows autonomous integration of the instrument data without the need for time consuming manual interpretation and selection of updated model parameters. Furthermore, only isotropic linear elasticity is considered in this paper, the approach presented is very flexible as to what type of geomechanical constitutive response can be used. The approach is easily adaptable to nonlinear physics-based constitutive models to account for common rock behaviors such as creep and plasticity. The approach also enables training of machine learning-based constitutive models by allowing back propagation of errors through the finite element calculations. This enables strongly enforcing known physics, such as conservation of momentum and continuity, while allowing data-driven models to learn the truly unknown physics such as the constitutive or petrophysical responses.

Burghardt, Jeffrey A.↗
DOE OSTI2021

NRAP-Open-IAM: Open Wellbore Component (V.2.0)

The Open Wellbore component of NRAP-Open-IAM is applicable to the calculation of Area of Review at a carbon storage site and fast estimation of dense gas dispersion from multiple continuous CO 2 surface leakage sources for risk assessment. The Open Wellbore Component of NRAP-Open-IAM is a lookup table model based on the drift-flux approach. The National Research Assessment Partnership’s second-generation integrated assessment model, NRAP-Open-IAM, is open-source software written in Python for use in performance and quantitative risk assessment of geologic carbon storage (GCS) systems. The look-up table for open well leakage is an updated version the previous model which provides finer resolution of reservoir depth and more accurate fluid properties than the previous open well look-up table, especially in calculating the gas (CO 2 -rich) phase properties. The component model input parameters include the reservoir transmissivity, aquifer transmissivity, brine salinity, well radius, and well top and bottom depths. The reservoir transmissivity and brine salinity were set to values consistent with the reservoir component, and the well top and bottom depths varied with location. The possible outputs from the Open Wellbore component are leakage rates of CO 2 and brine to either an aquifer or atmosphere, depending on the depth of the well top. Several python test scripts are presented that demonstrate that the Open Wellbore component works correctly. Use of the open wellbore model in NRAP-Open-IAM may result in a large CO 2 leakage rates, comparable to the leakage rates of CO 2 blow out.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Hydrological connectivity: a review and emerging strategies for integrating measurement, modeling, and management

This review synthesizes methods for measuring, modeling, and managing hydrologic connectivity, offering pathways to improve practices and address environmental challenges (e.g., climate change) and sustainability. As a key driver of water movement and nutrient cycling, hydrologic connectivity influences flood mitigation, water quality regulation, and biodiversity conservation. However, traditional field-based methods (e.g., dye tracing), indirect measurements (e.g., runoff analysis), and remote sensing techniques (e.g., InSAR) often struggle to capture the complexity of catchment-scale interactions. Similarly, modeling approaches—including process-based and percolation theory-based models, graph theory, and entropy-based metrics—face limitations in fully representing these interconnected processes. Both modeling and measurement techniques are constrained by inadequate spatial and temporal coverage, high data demands, computational complexity, and difficulties in representing subsurface connectivity. Subsequently, we critique current management practices that prioritize isolated variables (e.g., streamflow, sediment transport) over system-wide strategies and emphasize the need for adaptive, connectivity-based approaches in water resource planning and restoration. Moving forward, we highlight the importance of interdisciplinary collaboration, technological innovations (e.g., AI-driven modeling, real-time monitoring), and integrated frameworks to improve connectivity measurement, modeling, and adaptive management to restore fragmented hydrologic networks. This integrated approach sets the stage for transformative water resource management, fostering proactive policy development and stakeholder engagement.

Dwivedi, Dipankar↗
DOE OSTI2022

Physics Informed Neural Networks as Computational Physics Emulators

This report is a brief overview and evaluation of Physics Informed Neural Networks (PINNs). Karniadakis and co-workers, e.g., Karniadakis et al. (2021) assert that the PINNs approach integrates seamlessly both data and mathematical physics models, even in partially understood, uncertain and high-dimensional contexts. They further claim that PINNs are effective and efficient for ill-posed and inverse problems, and when combined with domain decomposition, are scalable to large problems and a tool to discover hidden physics. While they demonstrate the capabilities in specific academic instances, their overarching claims about PINNs seem to be an overstatement, at least at the current time. We have briefly considered a few of the limitations of PINNs in this investigation. It is not clear to us if the PINNs approach can ever be competitive with approaches that use specialized algorithms to achieve high-accuracy solutions of governing equations and other techniques that can combine observational data with such solutions. As an example of the latter, consistent with the principles of Bayesian inference, data assimilation, or more generally data-model fusion, is a process that fuses observational data typically with a computational model that respects certain constraints such as conservation laws. For example, improvements in observational network combined with data assimilation have been key in improving weather predictions over the past four decades Kalnay (2003).

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2024

Applications of Remote Sensing for Land Use Planning Scenarios with Suitability Analysis

In regions undergoing rapid urbanization, such as West Africa, land use planning (LUP) is vital to accommodate growing population and manage natural resources. Suitability analysis modeling is a widely used tool in LUP to determine the extent to which a land area is suitable for a designated purpose, but there is a gap in the integration of remote sensing time series data into land use decisions. The goal of this study was to incorporate remote sensing time series information with suitability analyses to inform LUP decisions in urban areas. In the study area of Kumasi, Ghana, land cover trends and land surface temperature (LST) from 2000 to 2019 were used to understand climate change trends. Suitability analyses determined the fitness of land areas for predetermined uses. These background processes informed a genetic algorithm to project plausible futures for three land use scenarios. One scenario represented current land use planning practices for addressing population growth, another scenario prioritized minimizing climate change impacts while also accommodating population growth, and the final scenario focused on both of these climate and population goals in addition to high density urban development. Each of these scenarios was successful in achieving population accommodation and respective climate change mitigation goals. The results for these scenarios provide insight into plausible land use distributions in 2050 based on different planning approaches. The genetic algorithm was able to effectively develop results for each scenario through the integration of remotely sensed trends and suitability models, providing a novel approach to land use decision-making.

remote sensing time series↗
NASA NTRS2013

Goal-Function Tree Modeling for Systems Engineering and Fault Management

The draft NASA Fault Management (FM) Handbook (2012) states that Fault Management (FM) is a "part of systems engineering", and that it "demands a system-level perspective" (NASAHDBK- 1002, 7). What, exactly, is the relationship between systems engineering and FM? To NASA, systems engineering (SE) is "the art and science of developing an operable system capable of meeting requirements within often opposed constraints" (NASA/SP-2007-6105, 3). Systems engineering starts with the elucidation and development of requirements, which set the goals that the system is to achieve. To achieve these goals, the systems engineer typically defines functions, and the functions in turn are the basis for design trades to determine the best means to perform the functions. System Health Management (SHM), by contrast, defines "the capabilities of a system that preserve the system's ability to function as intended" (Johnson et al., 2011, 3). Fault Management, in turn, is the operational subset of SHM, which detects current or future failures, and takes operational measures to prevent or respond to these failures. Failure, in turn, is the "unacceptable performance of intended function." (Johnson 2011, 605) Thus the relationship of SE to FM is that SE defines the functions and the design to perform those functions to meet system goals and requirements, while FM detects the inability to perform those functions and takes action. SHM and FM are in essence "the dark side" of SE. For every function to be performed (SE), there is the possibility that it is not successfully performed (SHM); FM defines the means to operationally detect and respond to this lack of success. We can also describe this in terms of goals: for every goal to be achieved, there is the possibility that it is not achieved; FM defines the means to operationally detect and respond to this inability to achieve the goal. This brief description of relationships between SE, SHM, and FM provide hints to a modeling approach to provide formal connectivity between the nominal (SE), and off-nominal (SHM and FM) aspects of functions and designs. This paper describes a formal modeling approach to the initial phases of the development process that integrates the nominal and off-nominal perspectives in a model that unites SE goals and functions of with the failure to achieve goals and functions (SHM/FM).

Johnson, Stephen B.↗
DOE OSTI2024

New ASME Section III, Division 5 Creep-Fatigue Design Rules Based on EPP and SMT Approaches

The integrated elastic-perfectly plastic (EPP) and simplified model test (SMT) creep-fatigue (CF) design methodology, referred to as the EPP-SMT method, is being developed as an alternative for CF evaluation in the design of pressure boundary components for high-temperature reactors. This report reviews the conceptual basis of the EPP-SMT methodology, summarizes the SMT experimental development efforts and results, and provides the technical basis for finalizing the EPP-SMT CF design curves for Alloy 617, based on a combined experimental and numerical approach conducted in FY 2024. This report presents the effect of hold time on the CF design curves for Alloy 617 at elevated temperatures. It includes proposed EPP-SMT CF design curves and tabulated values for continuous cycling, along with the effects of maximum hold time, for the use of this EPP-SMT CF evaluation method.

36 MATERIALS SCIENCE↗
NASA NTRS1988

A constitutive material model for nonlinear finite element structural analysis using an iterative matrix approach

A unified numerical method for the integration of stiff time dependent constitutive equations is presented. The solution process is directly applied to a constitutive model proposed by Bodner. The theory confronts time dependent inelastic behavior coupled with both isotropic hardening and directional hardening behaviors. Predicted stress-strain responses from this model are compared to experimental data from cyclic tests on uniaxial specimens. An algorithm is developed for the efficient integration of the Bodner flow equation. A comparison is made with the Euler integration method. An analysis of computational time is presented for the three algorithms.

Koenig, Herbert A.↗
DOE OSTI2026

Enabling Dynamic Probabilistic Risk Assessment of Physical Security Using EMRALD and MAAP

The optimization of physical security in nuclear power plants requires sophisticated methodologies that integrate operator actions and plant behavior through advanced simulation tools. To address this, Idaho National Laboratory [JL2.1]has developed the Modeling and Analysis for Safety and Security using the Dynamic EMRALD Framework (MASS-DEF) methodology, an approach that integrates force-on-force simulations, dynamic probabilistic risk assessment, and thermal-hydraulics modeling [JL3.1]to enhance security planning while reducing costs. We developed a tool that produces reduced order models using thermal hydraulic simulations from the Modular Accident Analysis Program (MAAP) [1]. These models can quickly evaluate reactor core behavior during attack simulations, and in so doing, address two barriers of traditional methods: (1) MAAP simulations are computationally intensive, and (2) attack scenarios must be run in a secure environment, which complicates analysis and validation. By precomputing scenario outcomes for a small number of modified parameters, the reduced order model significantly decreases the computational cost and enables offsite review of the results.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Advanced Fuels Campaign Execution Plan

The Advanced Fuels Campaign (AFC) Execution Plan details the strategy, mission, scope, and goals—both near-term and long-term—along with the structure and organization of nuclear fuels and materials research, development, and demonstration (RD&D) activities within the Fuel Cycle Technologies (FCT) program. The FCT program, tasked by the U.S. Department of Energy (DOE), employs a science-based approach to advance fuel technologies. This approach integrates theory, experiments, and multi-scale modeling and simulation (M&S) to develop a predictive understanding of fuel fabrication processes and fuel/cladding performance under irradiation, moving beyond traditional empirical methods. The long-term goals of the AFC are guided by the AFC Strategic Plan and align with the DOE Office of Nuclear Energy (NE) Roadmap [1], which outlines a multi-decade vision for demonstrating and qualifying advanced fuel forms to support diverse fuel cycle options. Near-term goals focus on enhancing accident tolerant fuels (ATF) for Light Water Reactors (LWR), a significant challenge that demands balancing immediate objectives with ongoing progress toward advanced reactor missions. Accelerating the traditional fuel qualification process to meet ATF objectives is another critical challenge. A detailed set of 5-year goals, summarized below, has been developed in line with the overarching science-based fuel development approach: • Advanced LWR Fuel Technologies: By 2027, support the development of advanced LWR fuel technologies with improved performance and enhanced accident tolerance. This includes high burnup (HBu), low enriched uranium (LEU)+, coated cladding, and doped fuel, aimed at complementing industry-led significant LWR uprates and plant refurbishments. • Tristructural Isotropic (TRISO) Fuel: Achieve qualification by 2028 and develop improved designs for emerging markets. • Metal Fuel: Achieve qualification by 2028 and develop improved designs for emerging markets. • Molten Salt Fuel: By 2027, deploy a robust program that enables fuel salt qualification technologies needed to support fuel salt research and development (R&D), focusing on emergent needs to derisk fuel salt production and utilization in advanced reactors. • Long-Term ATF: Develop fuel technologies that enable significant power uprates (~50%) in refurbished or new LWRs while optimizing fissile material utilization and waste disposal. The 5-year milestones in the AFC Execution Plan are contingent on an assumed budget. This Execution Plan will be updated annually to reflect actual funding profiles as budget guidance becomes available, ensuring milestones are adjusted accordingly. In summary, the AFC Execution Plan presents a comprehensive strategy to advance nuclear fuel technologies through a science-based approach, addressing both near-term and long-term goals while adapting to funding realities.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS2003

Automated Environment Generation for Software Model Checking

A key problem in model checking open systems is environment modeling (i.e., representing the behavior of the execution context of the system under analysis). Software systems are fundamentally open since their behavior is dependent on patterns of invocation of system components and values defined outside the system but referenced within the system. Whether reasoning about the behavior of whole programs or about program components, an abstract model of the environment can be essential in enabling sufficiently precise yet tractable verification. In this paper, we describe an approach to generating environments of Java program fragments. This approach integrates formally specified assumptions about environment behavior with sound abstractions of environment implementations to form a model of the environment. The approach is implemented in the Bandera Environment Generator (BEG) which we describe along with our experience using BEG to reason about properties of several non-trivial concurrent Java programs.

Tkachuk, Oksana↗
DOE OSTI2026

SAVY-4000 Finite-Element Drop Test Analysis

PFE Auxiliary Systems conducted drop testing on SAVY-4000 containers to evaluate structural response under 12-foot drop conditions. In support of that effort, a finite-element modeling capability was developed to simulate drop response across multiple container sizes and impact orientations. The purpose of this work was to provide a consistent analysis framework that could support interpretation of testing, compare response trends across multiple configurations, and generate quantities of interest for later comparison with experimental data. More broadly, the analysis and testing were intended to assess whether the containers continued to perform their primary function after a 12-foot drop, namely maintaining structural integrity and containment of the contents. The modeling approach combined an implicit preload analysis with an explicit drop simulation so that each drop event began from a mechanically realistic assembled condition, including compression of the silicone O-ring. Separate models were developed for 2-quart, 5-quart, 12-quart, and 10-gallon containers. The results were evaluated in terms of strain-gauge response, collar-lid gap behavior, and accumulated plastic strain. In addition, parametric studies were performed on the 2-quart container to assess sensitivity to O-ring stiffness, friction, canister thickness, geometry tolerance, and mesh density. The simulations showed that predicted drop responses depended strongly on both container size and drop orientation. Gap metrics identified cases in which the predicted collar-lid opening exceeded the nominal O-ring cross-section threshold, while plastic strain metrics identified localized regions of elevated permanent deformation. Parametric studies showed that the predicted response was especially sensitive to the assumed O-ring stiffness and contact friction, while the geometry tolerance study produced smaller changes in the cases examined. The main value of this work was that it established a repeatable modeling and simulation workflow to support drop-test implementation, evaluate effects of future configuration changes, and understand modeling assumptions that most influenced predicted response. At the current stage, the results were viewed as preliminary model predictions rather than validated predictions. The next step would be to compare drop-test data to the model so that predictive values of the workflow could be refined and used with greater confidence to assess whether the containers maintained structural integrity and containment of the contents after a 12-foot drop.

42 ENGINEERING↗
NASA NTRS2021

Mission Concept Design for Autonomous Space Missions using Mission-Level Modeling and Simulation

NASA’s Europa Lander mission is to search for biosignatures on Europa based on in-situ science using a lander architecture. This mission presents a set of challenges that requires a high level of autonomy on the lander system, leading to the need for a new operational paradigm that supports better collaboration and coordination between the lander and the ground operations team. M\&S is used for both designing the onboard system-level autonomy and the ground operations paradigm that allows effective and efficient collaboration and coordination between the lander and the ground operations team. In this paper, M\&S as applied to the design of new mission and operational concepts will be discussed. At its current early stage of the mission development for Europa Lander, the M\&S is used to explore different mission concepts and gain insights for design (formative) rather than to verify and validate fully designed mission concepts quantitatively (summative). Organically, we established a new approach to mission and operational concept exploration using high-fidelity modeling and simulation. M\&S has been an integral part of the approach of defining constraints and assertions, designing mission concepts, assessing (i.e., simulating them), and discovering insights, which feeds back to the definition an design steps. This organically-established approach provided important benefits to the project at its early phase of the development by enabling the project team to be able to build shared understanding of impacts from design characteristics, constraints, and their interactions on the mission performance.

Ye, Sean↗
DOE OSTI2026

Disruption modelling for engineering and physics design of ST-E1 fusion power plant

Plasma disruptions represent a critical challenge for high-performance tokamak operations, as they can compromise machine integrity and reduce operational availability. Although future fusion devices essentially need to incorporate strategies to minimise disruption occurrence, complete avoidance remains unattainable. Consequently, assessing and characterising unmitigated disruption consequences is fundamental for the design and qualification of next-generation fusion power plants. This work supports the pre-conceptual design of ST-E1, a low aspect-ratio Tokamak Fusion Power Plant developed by Tokamak Energy Ltd., by presenting a comprehensive disruption modelling approach applied across different design stages. The methodology integrates both physics and engineering considerations to evaluate the impact of disruptions on machine performance and structural integrity. From an engineering perspective, several ST-E1 layout options were analysed to investigate the electromagnetic response of key components under disruption-induced loads, enabling comparison between alternative design solutions. On the physics side, a broad set of disruption scenarios was explored, scanning operational space parameters, plasma-material interactions, and associated thermal loads. Furthermore, the study examined variations in disruption behaviour arising from different reference equilibria, focusing on a range starting from Double Null to Single Null configurations, reflecting the increasing up-down asymmetry consequences. The results reveal significant contrasts in plasma dynamics and structures electromagnetic behaviour between configurations, highlighting the importance of disruption modelling in guiding design choices. These analyses have proven instrumental in shaping ST-E1 development, offering critical insights for mitigating risks and optimising future fusion power plant designs.

Borowiec, Katarzyna [ORNL] (ORCID:0000000335911739↗
NASA NTRS2007

Space Shuttle Propulsion Systems Plume Modeling and Simulation for the Lift-Off Computational Fluid Dynamics Model

This paper details advances being made in the development of Reynolds-Averaged Navier-Stokes numerical simulation tools, models, and methods for the integrated Space Shuttle Vehicle at launch. The conceptual model and modeling approach described includes the development of multiple computational models to appropriately analyze the potential debris transport for critical debris sources at Lift-Off. The conceptual model described herein involves the integration of propulsion analysis for the nozzle/plume flow with the overall 3D vehicle flowfield at Lift-Off. Debris Transport Analyses are being performed using the Shuttle Lift-Off models to assess the risk to the vehicle from Lift-Off debris and appropriately prioritized mitigation of potential debris sources to continue to reduce vehicle risk. These integrated simulations are being used to evaluate plume-induced debris environments where the multi-plume interactions with the launch facility can potentially accelerate debris particles toward the vehicle.

Strutzenberg, L. L.↗
DOE OSTI2024

Manufactured solutions for an electromagnetic slot model

The accurate modeling of electromagnetic penetration is an important topic in computational electromagnetics. Electromagnetic penetration occurs through intentional or inadvertent openings in an otherwise closed electromagnetic scatterer, which prevent the contents from being fully shielded from external fields. To efficiently model electromagnetic penetration, aperture or slot models can be used with surface integral equations to solve Maxwell's equations. A necessary step towards establishing the credibility of these models is to assess the correctness of the implementation of the underlying numerical methods through code verification. Surface integral equations and slot models yield multiple interacting sources of numerical error and other challenges, which render traditional code-verification approaches ineffective. In this paper, we provide approaches to separately measure the numerical errors arising from these different error sources for the method-of-moments implementation of the electric-field integral equation with a slot model. Finally, we demonstrate the effectiveness of these approaches for a variety of cases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRSDate not supplied

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

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

Peter A Liever↗