Search NASA⌕ Search

SEARCH · Search NASA

Results for “Complex Framework Materials”

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 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS2016

Tangle-Free Mesh Motion for Ablation Simulations

Problems involving mesh motion-which should not be mistakenly associated with moving mesh methods, a class of adaptive mesh redistribution techniques-are of critical importance in numerical simulations of the thermal response of melting and ablative materials. Ablation is the process by which material vaporizes or otherwise erodes due to strong heating. Accurate modeling of such materials is of the utmost importance in design of passive thermal protection systems ("heatshields") for spacecraft, the layer of the vehicle that ensures survival of crew and craft during re-entry. In an explicit mesh motion approach, a complete thermal solve is first performed. Afterwards, the thermal response is used to determine surface recession rates. These values are then used to generate boundary conditions for an a posteriori correction designed to update the location of the mesh nodes. Most often, linear elastic or biharmonic equations are used to model this material response, traditionally in a finite element framework so that complex geometries can be simulated. A simple scheme for moving the boundary nodes involves receding along the surface normals. However, for all but the simplest problem geometries, evolution in time following such a scheme will eventually bring the mesh to intersect and "tangle" with itself, inducing failure. This presentation demonstrates a comprehensive and sophisticated scheme that analyzes the local geometry of each node with help from user-provided clues to eliminate the tangle and enable simulations on a wide-class of difficult problem geometries. The method developed is demonstrated for linear elastic equations but is general enough that it may be adapted to other modeling equations. The presentation will explicate the inner workings of the tangle-free mesh motion algorithm for both two and three-dimensional meshes. It will show abstract examples of the method's success, including a verification problem that demonstrates its accuracy and correctness. The focus of the presentation will be on the algorithm; specifics on how the techniques may be used in spacecraft design will be not discussed.

Droba, Justin↗
NASA NTRS2018

Gas-Granular Simulation Framework for Spacecraft Landing Plume-Surface Interaction and Debris Transport Analysis

The Gas-Granular Flow Solver (GGFS) multi-phase flow computational framework has been developed to enable simulations of particle flows complex extra-terrestrial regolith materials. Particle flows of interest include the damage of unprepared landing sites from rocket plume impingement on Moon, Mars, and asteroids. The flow solver implements an Eulerian-Eulerian two-fluid model with fluid representation of the gas phase and granular phase to avoid the need to model billions of particle interactions. The granular phase is modeled as an Eulerian fluid with constituent physics closure models derived from first-principle Discrete Element Model (DEM) particle interaction simulations that capture the complex, non-linear granular particle interaction effects. Granular phase constituent models have been developed and integrated that address the complex, non-linear granular material mechanics complexities resulting from both: the irregular, jagged particle shapes and poly-disperse mixture effects encountered in extra-terrestrial regolith, with lunar regolith as the extreme. The GGFS capabilities are being integrated into a proven NASA plume-surface interaction and debris transport simulation framework featuring the Loci/CHEM CFD program and Debris Transport Analysis (DTA) post-processing tools for applications in robotic and human Moon and Mars lander development. Integration of the three simulation tool components. Loci/CHEM, GGFS, and DTA, into a coordinated simulation framework will enable time-accurate spacecraft landing simulations that account for the alteration of the landing surface through plume-induced cratering and the resulting redirection of plume impingement flow and debris transport. Initial implementation of this simulation framework and application examples will be presented.

Liever, Peter A.↗
NASA NTRSDate not supplied

Affordable Development Strategy for NEP Nuclear Systems

One nuclear electric propulsion (NEP) reactor systems under consideration is a hydride moderated thermal spectrum reactor fueled by high assay low enriched uranium (HALEU). While such a reactor is expected to yield the lightest HALEU reactor design, its development challenges grow exponentially with increasing mission demands, most notably power output, specific weight of the overall system (which may require operation at temperatures exceeding 1200 K), service lifetime, and human-rated reliability. Two of the greatest cost drivers are full-powered nuclear demonstrations and extensive material development campaigns, so it is important to consider options that can minimize the need for or complexity of such tasks. This paper discusses a structured framework being developed for assessing how NEP design choices, such as materials selection, neutronic features, and heat-removal technologies, can translate into project risk and how project performance goals can be traded with development cost.Reactors operating at high temperatures often require cutting-edge heat transfer technologies and creep-resistant materials. Use of new materials in high temperature reactors brings additional complication beyond those common to any new space materials development campaign. For example, such materials may not possess necessary neutronic cross-sectional or neutronic irradiation data. Similarly, use of new materials may significantly influence core neutronics; in some extreme cases, neutronic reactivity feed-back of certain new materials can vary during their service life as radiation damage impacts the scattering cross-section. In an affordable development approach, high fidelity modeling and simulation tools are used to identify and characterize potential ‘knees-in-the-curves’ in the relationship that exists between the mission characteristics and the project risk. Of particular significance is use of modern uncertainty management and variance reduction methods to perform gap analyses that feed into phenomena identification and ranking tables (PIRT) commonly used to communicate nuclear readiness levels. Model-based measurements techniques are used to design sub-scale experiments as a substitute to minimize orcompletely eliminate the need for nuclear demonstrations.This paper will describe the approach and present preliminary results. It will lay the groundwork for developing a set of metrics that can be broadly characterized as system nuclear readiness levels and advancement degree of difficulty for nuclear systems. Equally importantly, a goal of this paper is to initiate a dialogue among stakeholders on what is the sufficient level of maturity that is required for launching a demonstration unit.

Dasari V Rao↗
NASA NTRS2016

Revisiting of Multiscale Static Analysis of Notched Laminates Using the Generalized Method of Cells

Composite material systems generally exhibit a range of behavior on different length scales (from constituent level to macro); therefore, a multiscale framework is beneficial for the design and engineering of these material systems. The complex nature of the observed composite failure during experiments suggests the need for a three-dimensional (3D) multiscale model to attain a reliable prediction. However, the size of a multiscale three-dimensional finite element model can become prohibitively large and computationally costly. Two-dimensional (2D) models are preferred due to computational efficiency, especially if many different configurations have to be analyzed for an in-depth damage tolerance and durability design study. In this study, various 2D and 3D multiscale analyses will be employed to conduct a detailed investigation into the tensile failure of a given multidirectional, notched carbon fiber reinforced polymer laminate. Threedimensional finite element analysis is typically considered more accurate than a 2D finite element model, as compared with experiments. Nevertheless, in the absence of adequate mesh refinement, large differences may be observed between a 2D and 3D analysis, especially for a shear-dominated layup. This observed difference has not been widely addressed in previous literature and is the main focus of this paper.

Materials Engineering↗
NASA NTRS1990

CFD spinoff - Computational electromagnetics for radar cross section (RCS) studies

A finite-volume discretization procedure derived from proven CFD methods is used to solve the conservation form of the time-domain Maxwell's equations, in order to compute EM scattering from layered objects. This time-domain approach handles both single-frequency/continuous wave and broadband-frequency/pulse incident excitation. Arbitrarily shaped objects are modeled by means of a body-fitted coordinate transformation; complex internal/external structures with many material layers are treated through the implementation of a multizone framework capable of handling any type of zonal boundary condition. Results are presented for various two- and three-dimensional problems.

Shankar, Vijaya↗
NASA NTRS2007

Quantum Approximate Methods for the Atomistic Modeling of Multicomponent Alloys

This chapter describes the role of quantum approximate methods in the understanding of complex multicomponent alloys at the atomic level. The need to accelerate materials design programs based on economical and efficient modeling techniques provides the framework for the introduction of approximations and simplifications in otherwise rigorous theoretical schemes. As a promising example of the role that such approximate methods might have in the development of complex systems, the BFS method for alloys is presented and applied to Ru-rich Ni-base superalloys and also to the NiAI(Ti,Cu) system, highlighting the benefits that can be obtained from introducing simple modeling techniques to the investigation of such complex systems.

Bozzolo, Guillermo↗
NASA NTRS2009

The Role of Guidance, Navigation, and Control in Hypersonic Vehicle Multidisciplinary Design and Optimization

Airbreathing hypersonic systems offer distinct performance advantages over rocket-based systems for space access vehicles. However, these performance advantages are dependent upon advances in current state-of-the-art technologies in many areas such as ram/scramjet propulsion integration, high temperature materials, aero-elastic structures, thermal protection systems, transition to hypersonics and hypersonic control elements within the framework of complex physics and new design methods. The complex interactions between elements of an airbreathing hypersonic vehicle represent a new paradigm in vehicle design to achieve the optimal performance necessary to meet space access mission objectives. In the past, guidance, navigation, and control (GNC) analysis often follows completion of the vehicle conceptual design process. Individual component groups design subsystems which are then integrated into a vehicle configuration. GNC is presented the task of developing control approaches to meet vehicle performance objectives given that configuration. This approach may be sufficient for vehicles where significant performance margins exist. However, for higher performance vehicles engaging the GNC discipline too late in the design cycle has been costly. For example, the X-29 experimental flight vehicle was built as a technology demonstrator. One of the many technologies to be demonstrated was the use of light-weight material composites for structural components. The use of light-weight materials increased the flexibility of the X- 29 beyond that of conventional metal alloy constructed aircraft. This effect was not considered when the vehicle control system was designed and built. The impact of this is that the control system did not have enough control authority to compensate for the effects of the first fundamental structural mode of the vehicle. As a result, the resulting pitch rate response of the vehicle was below specification and no post-design changes could recover the desired capability.

Ouzts, Peter J.↗
NASA NTRS2001

Reverse Micelle Based Synthesis of Microporous Materials in Microgravity

Microporous materials include a large group of solids of varying chemical composition as well as porosity. These materials are characterized by channels and cavities of molecular dimensions. The framework structure is made up of interconnecting T-O-T' bonds, where T and T' can be Si, Al, P, Ga, Fe, Co, Zn, B and a host of other elements. Materials with Si-O-Al bonding in the framework are called zeolites and are extensively used in many applications. Ion-exchange properties of these materials are exploited in the consumer and environmental industries. Chemical and petroleum industries use zeolites as catalysts in hydrocarbon transform ations. Synthesis of new microporous frameworks has led to the development of new technologies, and thus considerable effort worldwide is expended in their discovery. Microporous materials are typically made under hydrothermal conditions. Influence of nature of starting reactants, structure directing agents, pH, temperature, and aging all have profound influence on the synthesis process. This is primarily because the most interesting open frameworks are not necessarily the stable structures in the reaction medium. Thus, the discovery of new frameworks is often tied to finding the right composition and synthesis conditions that allow for kinetic stabilization of the structure. This complexity of the synthesis process and limited understanding of it has made it difficult to develop directed is of microporous materials and most advances in this field have been made by trial and error. The basic issues in crystal growth of these materials include: (1) Nature of the nucleation process; (2) Molecular structure and assembly of nuclei; (3) Growth of nuclei into crystals; (4) Morphology control; and (5) Transformation of frameworks into other structures. The NASA-funded research described in this paper focuses on all the above issues and has been described in several publications. We present the highlights of our program, especially with the focus on possible experiments in microgravity.

Dutta, Prabir K.↗
NASA NTRS2002

Fundamental Studies of Crystal Growth of Microporous Materials

Microporous materials are framework structures with well-defined porosity, often of molecular dimensions. Zeolites contain aluminum and silicon atoms in their framework and are the most extensively studied amongst all microporous materials. Framework structures with P, Ga, Fe, Co, Zn, B, Ti and a host of other elements have also been made. Typical synthesis of microporous materials involve mixing the framework elements (or compounds, thereof) in a basic solution, followed by aging in some cases and then heating at elevated temperatures. This process is termed hydrothermal synthesis, and involves complex chemical and physical changes. Because of a limited understanding of this process, most synthesis advancements happen by a trial and error approach. There is considerable interest in understanding the synthesis process at a molecular level with the expectation that eventually new framework structures will be built by design. The basic issues in the microporous materials crystallization process include: (1) Nature of the molecular units responsible for the crystal nuclei formation; (2) Nature of the nuclei and nucleation process; (3) Growth process of the nuclei into crystal; (4) Morphological control and size of the resulting crystal; (5) Surface structure of the resulting crystals; (6) Transformation of frameworks into other frameworks or condensed structures. The NASA-funded research described in this report focuses to varying degrees on all of the above issues and has been described in several publications. Following is the presentation of the highlights of our current research program. The report is divided into five sections: (1) Fundamental aspects of the crystal growth process; (2) Morphological and Surface properties of crystals; (3) Crystal dissolution and transformations; (4) Modeling of Crystal Growth; (5) Relevant Microgravity Experiments.

Dutta, P.↗
NASA NTRS2003

Fundamental Studies of Crystal Growth of Microporous Materials

Microporous materials are framework structures with well-defined porosity, often of molecular dimensions. Zeolites contain aluminum and silicon atoms in their framework and are the most extensively studied amongst all microporous materials. Framework structures with P, Ga, Fe, Co, Zn, B, Ti and a host of other elements have also been made. Typical synthesis of microporous materials involve mixing the framework elements (or compounds, thereof) in a basic solution, followed by aging in some cases and then heating at elevated temperatures. This process is termed hydrothermal synthesis, and involves complex chemical and physical changes. Because of a limited understanding of this process, most synthesis advancements happen by a trial and error approach. There is considerable interest in understanding the synthesis process at a molecular level with the expectation that eventually new framework structures will be built by design. The basic issues in the microporous materials crystallization process include: (a) Nature of the molecular units responsible for the crystal nuclei formation; (b) Nature of the nuclei and nucleation process; (c) Growth process of the nuclei into crystal; (d) Morphological control and size of the resulting crystal; (e) Surface structure of the resulting crystals; and (f) Transformation of frameworks into other frameworks or condensed structures.

Singh, Ramsharan↗
NASA NTRS2019

Probabilistic Risk Assessment of Aging Layered Pressure Vessels

The National Aeronautics and Space Administration (NASA) operates approximately 300 aging layered pressure vessels that were designed and manufactured prior to ASME Boiler and Pressure Vessel (B&PV) code requirements. In order to make decisions regarding the continued fitness-for-service of these non-code carbon steel vessels, it is necessary to perform a relative risk of failure assessment for each vessel. However, risk assessment of these vessels is confounded by uncertainties and variabilities related to the use of proprietary materials in fabrication, missing construction records, geometric discontinuities, weld residual stresses, and complex service stress gradients in and around the welds. Therefore, a probabilistic framework that can capture these uncertainties and variabilities has been developed to assess the fracture risk of flaws in regions of interest, such as longitudinal and circumferential welds, using the NESSUS® probabilistic modeling software and NASGRO® fracture mechanics software. In this study, the probabilistic framework was used to predict variability in the stress intensity factor associated with different reference flaws located in the head-to-shell circumferential welds of a 4-layer and 14-layer pressure vessel. The probabilistic studies predict variability in flaw behavior and the important uncertain parameters for each reference flaw location.

Riha, D. S.↗
NASA NTRS2010

Spinoff 2010

Topics covered include: Burnishing Techniques Strengthen Hip Implants; Signal Processing Methods Monitor Cranial Pressure; Ultraviolet-Blocking Lenses Protect, Enhance Vision; Hyperspectral Systems Increase Imaging Capabilities; Programs Model the Future of Air Traffic Management; Tail Rotor Airfoils Stabilize Helicopters, Reduce Noise; Personal Aircraft Point to the Future of Transportation; Ducted Fan Designs Lead to Potential New Vehicles; Winglets Save Billions of Dollars in Fuel Costs; Sensor Systems Collect Critical Aerodynamics Data; Coatings Extend Life of Engines and Infrastructure; Radiometers Optimize Local Weather Prediction; Energy-Efficient Systems Eliminate Icing Danger for UAVs; Rocket-Powered Parachutes Rescue Entire Planes; Technologies Advance UAVs for Science, Military; Inflatable Antennas Support Emergency Communication; Smart Sensors Assess Structural Health; Hand-Held Devices Detect Explosives and Chemical Agents; Terahertz Tools Advance Imaging for Security, Industry; LED Systems Target Plant Growth; Aerogels Insulate Against Extreme Temperatures; Image Sensors Enhance Camera Technologies; Lightweight Material Patches Allow for Quick Repairs; Nanomaterials Transform Hairstyling Tools; Do-It-Yourself Additives Recharge Auto Air Conditioning; Systems Analyze Water Quality in Real Time; Compact Radiometers Expand Climate Knowledge; Energy Servers Deliver Clean, Affordable Power; Solutions Remediate Contaminated Groundwater; Bacteria Provide Cleanup of Oil Spills, Wastewater; Reflective Coatings Protect People and Animals; Innovative Techniques Simplify Vibration Analysis; Modeling Tools Predict Flow in Fluid Dynamics; Verification Tools Secure Online Shopping, Banking; Toolsets Maintain Health of Complex Systems; Framework Resources Multiply Computing Power; Tools Automate Spacecraft Testing, Operation; GPS Software Packages Deliver Positioning Solutions; Solid-State Recorders Enhance Scientific Data Collection; Computer Models Simulate Fine Particle Dispersion; Composite Sandwich Technologies Lighten Components; Cameras Reveal Elements in the Short Wave Infrared; Deformable Mirrors Correct Optical Distortions; Stitching Techniques Advance Optics Manufacturing; Compact, Robust Chips Integrate Optical Functions; Fuel Cell Stations Automate Processes, Catalyst Testing; Onboard Systems Record Unique Videos of Space Missions; Space Research Results Purify Semiconductor Materials; and Toolkits Control Motion of Complex Robotics.

Source record↗
NASA NTRS1992

Continuation of tailored composite structures of ordered staple thermoplastic material

The search for the cost effective composite structure has motivated the investigation of several approaches to develop composite structure from innovative material forms. Among the promising approaches is the conversion of a planar sheet to components of complex curvature through sheet forming or stretch forming. In both cases, the potential for material stretch in the fiber direction appears to offer a clear advantage in formability over continuous fiber systems. A framework was established which allows the simulation of the anisotropic mechanisms of deformation of long discontinuous fiber laminates wherein the matrix phase is a viscous fluid. Predictions for the effective viscosities of a hyper-anisotropic medium consisting of collimated, discontinuous fibers suspended in viscous matrix were extended to capture the characteristics of typical polymers including non-Newtonian behavior and temperature dependence. In addition, the influence of fiber misorientation was also modeled by compliance averaging to determine ensemble properties for a given orientation distribution. A design tool is presented for predicting the effect of material heterogeneity on the performance of curved composite beams such as those used in aircraft fuselage structures. Material heterogeneity can be induced during manufacturing processes such as sheet forming and stretch forming of thermoplastic composites. This heterogeneity can be introduced in the form of fiber realignment and spreading during the manufacturing process causing radial and tangential gradients in material properties. Two analysis procedures are used to solve the beam problems. The first method uses separate two-dimensional elasticity solutions for the stresses in the flange and web sections of the beam. The separate solutions are coupled by requiring that forces and displacements match section boundaries. The second method uses an approximate Rayleigh-Ritz technique to find the solutions for more complex beams. Analyses are performed for curved beams of various cross-sections loaded in pure bending and with a uniform distributed load. Preliminary results show that the geometry of the beam dictates the effect of heterogeneity on performance. The role of heterogeneity is larger in beams with a small average radius-to-depth ration, R/t, where R is the average radius of the beam and t is the difference between the inside and outside radii. Results of the anlysis are in the form of stresses and displacements and are compared to both mechanics of materials and numerical solutions obtained using finite element analysis.

Santare, Michael H.↗
NASA NTRSDate not supplied

CFD Predictions of Boiling Regime Transitions during Line Chilldown validated against a 1G LN2 Experiment

Introduction Before filling a propellant tank on the ground or in Space, the transfer line between the donor and receiver tanks must be cooled down preferably by sacrificing a minimum amount of the cryogenic fluid. The cryogenic line chill-down process involves a transition between different flow boiling regimes, namely, film boiling, transition film boiling, and nucleate boiling which are complex and may be quite gravity-dependent. Capturing these boiling phenomena and predicting the transition between them in a CFD framework is new and challenging both for 1g and microgravity applications. Materials & Methods The present work addresses this challenge by employing a two-phase Eulerian approach in the context of a homogeneous fluid mixture together with the Lee phase change model to capture the film boiling regime of the chill-down process using ANSYS Fluent®. The nucleate boiling regime is predicted by incorporating an in-house developed sub-grid model that accounts for bubble nucleation, bubble growth, bubble departure diameter, and their shedding frequency. The sub-grid model is encoded and implemented into Fluent via a user-defined function for the wall-fluid heat flux calculations. The mathematical formulation and numerical implementation of the CFD model are described in detail. The coupled CFD-Subgrid model is validated against published experimental data for liquid nitrogen chill-down of a heated stainless-steel pipe in 1g. Results Numerical simulation results show good agreements between the CFD predictions of the wall temperature evolution, rewetting temperature, and transition between film and nucleate boiling, with the experimental measurements published by Darr et al [2] for several different LN2 flowrates in the vertical pipe orientation. The CFD predictions for the wall temperature distribution indicate a rapid quenching of the wall at two upstream and downstream temperature sensing locations as compared to the experimental measurement. The only tuning parameter in the CFD model is the Lee mass transfer coefficient. The CFD Model predicts the Liedenfrost rewetting temperature in close agreement with the experiment. This marks a transition between stable and transitionary flow boiling regimes. The CFD-predicted boiling curve for the downstream sensor location is also compared against its experimental counterpart and indicates that the model is able to predict all the key temperature and heat flux parameters during the transitions from stable to transitionary film boiling to nucleate boiling in close agreement with the experiment. A sequence of predicted volume fraction, and temperature contours depicting these transitions will be presented.

Evaporation Condensation↗
NASA NTRS1973

Gravitation

This textbook on gravitation physics (Einstein's general relativity or geometrodynamics) is designed for a rigorous full-year course at the graduate level. The material is presented in two parallel tracks in an attempt to divide key physical ideas from more complex enrichment material to be selected at the discretion of the reader or teacher. The full book is intended to provide competence relative to the laws of physics in flat space-time, Einstein's geometric framework for physics, applications with pulsars and neutron stars, cosmology, the Schwarzschild geometry and gravitational collapse, gravitational waves, experimental tests of Einstein's theory, and mathematical concepts of differential geometry.

Misner, C. W.↗
NASA NTRS1983

Complexities of high temperature metal fatigue: Some steps toward understanding

After pointing out many of the complexities that attend high temperature metal fatigue beyond those already studied in the sub-creep range, a description of the micromechanisms of deformation and fracture is presented for several classes of materials that were studied over the past dozen years. Strainrange Partitioning (SRP) is used as a framework for interpreting the results. Several generic types of behavior were observed with regard both to deformation and fracture and each is discussed in the context of the micromechanisms involved. Treatment of cumulative fatigue damage and the possibility of ""healing'' of damage in successive loading loops, has led to a new interpretation of the Interaction Damage Rule of SRP. Using the concept of ""equivalent micromechanistic damage'' -- that the same damage on a microscopic scale is induced if the same hysteresis loops are generated, element for element -- it turns out the Interaction Damage Rule essentially compounds a number of variants of hysteresis loops, all of which have the same damage according to SRP concepts, into a set of loops each containing only one of the generic SRP strainranges. Thus the damage associcated with complex loops comprising several types of strainrange is analyzed by considering a combination of loops each containing only one type of strainrange. This concept is expanded to show how several independent loops can combine to ""heal'' creep damage in a complex loading history.

Manson, S. S.↗
NASA NTRS1997

Parallel adaptive mesh refinement techniques for plasticity problems

The accurate modeling of the nonlinear properties of materials can be computationally expensive. Parallel computing offers an attractive way for solving such problems; however, the efficient use of these systems requires the vertical integration of a number of very different software components, we explore the solution of two- and three-dimensional, small-strain plasticity problems. We consider a finite-element formulation of the problem with adaptive refinement of an unstructured mesh to accurately model plastic transition zones. We present a framework for the parallel implementation of such complex algorithms. This framework, using libraries from the SUMAA3d project, allows a user to build a parallel finite-element application without writing any parallel code. To demonstrate the effectiveness of this approach on widely varying parallel architectures, we present experimental results from an IBM SP parallel computer and an ATM-connected network of Sun UltraSparc workstations. The results detail the parallel performance of the computational phases of the application during the process while the material is incrementally loaded.

Barry, W. J.↗
NASA NTRSDate not supplied

Implementation of Charged Particle Behavior in Discrete Element Method (DEM) Simulations

Lunar landers will agitate the surface of the Moon with an exhaust plume during descent which will, due to the particulate nature of the lunar regolith, loosen and eject grains from the surface. This ejection is not only coupled with the charged plume gas, but also results in significant particle-particle interactions. Settling of these grains after plume effects have subsided takes much longer than expected in a ballistic trajectory. The prevailing hypothesis attributes this behavior to the accumulated charge on the particles. We are thus developing a discrete element method (DEM) approach to explore these charged particle interactions on the lunar surface. The Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) Improved for General Granular and Granular Heat Transfer Simulations (LIGGGHTS) software package provides a DEM modeling framework for granular interactions. It includes many complexities such as non-spherical particle shapes, cohesion and frictional forces, and heat transfer, but has no provision for inter-particle electrostatic forces and charge transfer that are important to examine in the lunar environment. In this work, a standard Coulomb potential and a Yukawa potential are integrated into the LIGGGHTS framework to provide a basis for particle-particle electrostatic interactions, as well as a gravitational potential to enable inter-grain gravitational interactions. A preliminary approach to charge transfer between grains incorporating properties such as work function and electrical conductivity to the library of available material characteristics will be presented. Several scenarios have been simulated that include charged particle interactions within a diffuse granular gas, settling of charged grains into a regolith bed, sliding of granular material along an incline, and vibration of settled grains to produce a behavior similar to fluidization. There are numerous challenges to incorporate realistic interactions between complex lunar particles. Currently, grains are modeled to behave as if the entirety of the charge acts at the center of mass, such as conductors with spherical symmetry and insulators with homogeneously distributed charge. We are developing improvements that will include the use of non-spherical particle geometries, as well as reasonable approximations of insulating/dielectric materials that have non-uniform charge distributions. The cases simulated thus far will be examined in a relevant environment within a vacuum chamber to validate the simulations. These simulations will be bounded by experiments utilizing high-speed camera observations of the motion for validation. The grains in the experiment will exchange charge during their motion and this can be quantified by collection within a charge measurement device such as a Faraday cup. Such a device may be modeled within the software by defining an integration region and computing the contained charge as a function of simulation time, allowing for side-by-side comparison of simulated and measured bulk charging results. Any differences will be reconciled by updating the mathematical mechanisms described within the simulation suite. Successfully combining results from experiments within a relevant environment into the LIGGGHTS framework will improve modeling of the charged grain dynamics experienced on the Moon to provide insights into dust behavior for future lunar exploration missions.

Electrostatics↗