Search NASA⌕ Search

SEARCH · Search NASA

Results for “discrete particle method”

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
DOE OSTI2020

Three-Dimensional Evaluation of Sand Particle Fracture Using Discrete-Element Method and Synchrotron Microcomputed Tomography Images

Recent research showed that fracture of sand particles plays a significant role in determining the plastic bulk volumetric changes of granular materials under different loading conditions. One of the major tools used to better understand the influence of particle fracture on the behavior of granular materials is discrete-element modeling (DEM). This paper employed the bonded block model (BBM) to simulate the fracture behavior of sand. Each sand particle is modeled as an agglomerate of rigid blocks bonded at their contacts using the linear-parallel contact model, which can transmit both moment and force. DEM simulated particles closely matched the actual three-dimensional (3D) shape of sand particles acquired using high-resolution 3D synchrotron microcomputed tomography (SMT). Results from unconfined one-dimensional (1D) compression of a single synthetic silica cube were used to calibrate the model parameters. Particle fracture was investigated for specimens composed of three sand particles that were loaded under confined 1D compression. Breakage energy measured from DEM models matched well with that measured experimentally. The paper studied the effects of contact loading condition and particle interaction on the fracture mode of particles using BBM that can closely capture the 3D shape of real sand particles.

58 GEOSCIENCES↗
NASA NTRS1990

Hypergolic bipropellant spray combustion and flow modelling in rocket engines

A predictive tool for hypergolic bipropellant spray combustion and flow evolution in small rocket combustion chambers is described. It encompasses a computational technique for the gas-phase governing equations, a discrete particle method for liquid bipropellant sprays, and constitutive models for combustion chemistry, interphase exchanges, and unlike impinging hypergolic spray interactions. Emphasis is placed on the phenomenological modeling of the hypergolic liquid bipropellant gasification processes. Sample computations with the N2H4-N2O4 propellant system are given in order to show some of the capabilities and inadequacies of this tool.

Larosiliere, Louis M.↗
NASA NTRS1990

Combustion and flow modelling applied to the OMV VTE

A predictive tool for hypergolic bipropellant spray combustion and flow evolution in the OMV VTE (orbital maneuvering vehicle variable thrust engine) is described. It encompasses a computational technique for the gas phase governing equations, a discrete particle method for liquid bipropellant sprays, and constitutive models for combustion chemistry, interphase exchanges, and unlike impinging liquid hypergolic stream interactions. Emphasis is placed on the phenomenological modelling of the hypergolic liquid bipropellant gasification processes. An application to the OMV VTE combustion chamber is given in order to show some of the capabilities and inadequacies of this tool.

Larosiliere, Louis M.↗
NASA NTRS2009

Hall-Effect Thruster Simulations with 2-D Electron Transport and Hydrodynamic Ions

A computational approach that has been used extensively in the last two decades for Hall thruster simulations is to solve a diffusion equation and energy conservation law for the electrons in a direction that is perpendicular to the magnetic field, and use discrete-particle methods for the heavy species. This "hybrid" approach has allowed for the capture of bulk plasma phenomena inside these thrusters within reasonable computational times. Regions of the thruster with complex magnetic field arrangements (such as those near eroded walls and magnets) and/or reduced Hall parameter (such as those near the anode and the cathode plume) challenge the validity of the quasi-one-dimensional assumption for the electrons. This paper reports on the development of a computer code that solves numerically the 2-D axisymmetric vector form of Ohm's law, with no assumptions regarding the rate of electron transport in the parallel and perpendicular directions. The numerical challenges related to the large disparity of the transport coefficients in the two directions are met by solving the equations in a computational mesh that is aligned with the magnetic field. The fully-2D approach allows for a large physical domain that extends more than five times the thruster channel length in the axial direction, and encompasses the cathode boundary. Ions are treated as an isothermal, cold (relative to the electrons) fluid, accounting for charge-exchange and multiple-ionization collisions in the momentum equations. A first series of simulations of two Hall thrusters, namely the BPT-4000 and a 6-kW laboratory thruster, quantifies the significance of ion diffusion in the anode region and the importance of the extended physical domain on studies related to the impact of the transport coefficients on the electron flow field.

plasma↗
NASA NTRS2010

A Review of Discrete Element Method (DEM) Particle Shapes and Size Distributions for Lunar Soil

As part of ongoing efforts to develop models of lunar soil mechanics, this report reviews two topics that are important to discrete element method (DEM) modeling the behavior of soils (such as lunar soils): (1) methods of modeling particle shapes and (2) analytical representations of particle size distribution. The choice of particle shape complexity is driven primarily by opposing tradeoffs with total number of particles, computer memory, and total simulation computer processing time. The choice is also dependent on available DEM software capabilities. For example, PFC2D/PFC3D and EDEM support clustering of spheres; MIMES incorporates superquadric particle shapes; and BLOKS3D provides polyhedra shapes. Most commercial and custom DEM software supports some type of complex particle shape beyond the standard sphere. Convex polyhedra, clusters of spheres and single parametric particle shapes such as the ellipsoid, polyellipsoid, and superquadric, are all motivated by the desire to introduce asymmetry into the particle shape, as well as edges and corners, in order to better simulate actual granular particle shapes and behavior. An empirical particle size distribution (PSD) formula is shown to fit desert sand data from Bagnold. Particle size data of JSC-1a obtained from a fine particle analyzer at the NASA Kennedy Space Center is also fitted to a similar empirical PSD function.

Lane, John E.↗
DOE OSTI2021

Thermal conductivity modeling of monodispersed microspheres using discrete element method

Particle beds are widely used in various systems and processes, such as particle heat exchangers, granular flow reactors, and additive manufacturing. Accurate modeling of the thermal conductivity of particle beds and understanding of their heat transfer mechanisms are important. However, previous models were based on simple cubic packing of particles, which could not accurately represent the actual heat transfer processes under certain conditions. Here, we examine the effect of the packing structure on the thermal conductivity of particle beds. We use monodispersed silica microspheres with average particle sizes ranging from 23 to 330 μm as a model material. We employ a transient hot-wire technique to measure the thermal conductivity of the particle beds with packing density of 43%–57% within a temperature range of room temperature to 500 °C and under N2 gaseous pressures of 20–760 Torr. We then use a discrete element method (DEM) to obtain the realistic packing structure of the particles, which is then fed into a finite-element model (FEM) to calculate the thermal conductivity, with the consideration of solid conduction, gas conduction, and radiation heat transfer. Our results show that the thermal conductivity model based on the more realistic random packing structure derived from the DEM shows better agreement with the experimental data compared to that based on the simple cubic-packing structure. The combined DEM and FEM methodology can serve as a useful tool to predict the effective thermal conductivity of particle beds and to quantify different heat transfer mechanisms under various conditions.

14 SOLAR ENERGY↗
NASA NTRSDate not supplied

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

To understand the behavior of charged lunar regolith when perturbed by lunar landers, it is important to couple the grain dynamics with mechanical and electrical particle interactions. To accomplish this, improvements have been made to a discrete element method (DEM) software package to include both short- and long-range interactions between spherical particles. Short-range interactions rely on contact between the particles, such as electrical conduction and triboelectric charge transfer driven by work functions. Long-range interactions act at a distance between every pairing of particles, such as electrostatic forces and gravitational forces. Results from simulations between a few particles are compared with theory to verify these added behaviors prior to scaling up to more complex scenarios. The radii, initial charges, electrical conductivities, work functions, and separation of the particles are varied and the resultant charges as well as the time required to reach the final state are determined.

Electrostatics↗
DOE OSTI2021

Applicability of discrete element method with spherical and clumped particles for constitutive study of granular materials

Discrete element method (DEM) has been intensively used to study the constitutive behaviour of granular materials. However, to what extent a real granular material can be reproduced by virtual DEM simulations remains unclear. This study attempts to answer this question by comparing DEM simulations with typical features of experimental granular materials. Three groups of models with spherical and clumped particles are investigated from four perspectives: (i) deviatoric stress and volumetric behaviour; (ii) critical state behaviour; (iii) stress-dilatancy relationship; and (iv) the evolution of principal stress ratio against axial strain. The results demonstrate that DEM with spherical or clumped particles is capable of qualitatively describing macroscopic deviatoric stress responses, volumetric behaviour, and critical state behaviour observed in experiments for granular materials. On the other hand, some qualitative deviations between experiments and the investigated DEM simulations are also observed, in terms of the stress-dilatancy behaviour and principal stress ratio against axial strain, which are proven to be critical for constitutive modelling. The results demonstrate that DEM with spherical or clumped particles may not necessarily fully capture experimental features of granular materials even from a qualitative perspective. It is thus encouraged to thoroughly validate DEM with experiments when developing constitutive models based on DEM observations.

58 GEOSCIENCES↗
DOE OSTI2023

Investigation of Best-Practices and Computationally Inexpensive Radiative Exchange Models for Discrete Element Method Modeling of Aluminosilicate Particles in Concentrating Solar Power Environments

Chemically inert, aluminosilicate based particles have been investigated as both a thermal transport and sensible energy storage medium for concentrating solar power facilities. These particles will experience a wide range of operating temperatures (300-1000 K) and handling conditions (dense to dilute falling particle curtains, dense granular flows, or dense structures), requiring specially-designed and optimized infrastructures. The relative influence of collisional and frictional interactions between particles varies based on temperature-dependent particulate properties and greatly impacts the bulk, granular flow behavior. These underlying physics are captured using discrete element method modeling tools. However, this modeling method is computationally expensive as each particle position and interaction is tracked during the simulation. These modeling methods are further complicated by introducing temperature-dependent particle properties, high-temperature radiative exchange, and directional irradiation sources experienced by granular flows in concentrating solar power environments. In this study, coupled experimental and numerical slump testing of aluminosilicate particles was performed and computationally efficient radiative exchange models were evaluated to establish best-practices for discrete element method models for concentrating solar power environments. The three particle types investigated included Carbobead HSP 30 /60, Carbobead CP 30/60, and Granusil 4030. Existing modeling limitations and computationally-efficient multi-modal heat transfer models were evaluated using Aspherix®, a commercial discrete element method software. High-temperature (< 1073 K) slump testing of aluminosilicate particles was performed to investigate the deviation between experimentally-observed and numerically-predicted angles of repose introduced by computation-time reduction practices including the relaxation of the particle elastic modulus and coarse-graining. Coarse-graining is used to use a single modeled particle that is representative of a collection of smaller particles, decreasing the computational cost at the expense of geometric accuracy. Additionally, relaxation of the elastic modulus is used to reduce computational time at the expense of an increased, modeled particle overlap. Prior studies have determined that aluminosilicate particles retain a high elastic modulus at high temperatures (< 1073 K), requiring small simulation timesteps to ensure resolved contact forces resemble appropriate solid mechanics. A parametric study was performed to evaluate the influence of computation time improvements on the deviation between experimental and modeled angle of repose across high temperatures < 1073 K. Additionally, numerical case studies were performed on candidate particle systems at varying porosities and temperatures. These studies were performed to investigate the influence of computationally-efficient radiative-exchange modeling methods coupled to Aspherix® on modeled accuracy and computation time. The recently-developed distance-based approximation was evaluated in estimating radiative exchange between particles and participating surfaces located in close proximity. The distance based approximation was developed to use tabulated estimates of the radiative distribution factor between individual particles and surfaces in close proximity (< 40 particle radii). These methods were expanded to the aluminosilicate particles of interest, including the influence of particle size distributions. To capture radiative exchange between particles and surfaces not in close proximity (> 40 particle radii) and to capture the absorption of directional irradiation from concentrating solar resources, a volumetrically-averaged radiative distribution factor was calculated between the modeled granular flow and surfaces using Monte Carlo ray-tracing for participating media. Volume-averaged absorption and scattering coefficients were predicted using a volumetric discretization of the modeled domain with monodisperse approximations based on geometric optics and experimentally-determined scattering phase functions for aluminosilicate particles.

14 SOLAR ENERGY↗
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↗
DOE OSTI2025

Coarse Graining Discrete Element Method Information in Particle-in-Cell Length Scales Using a Machine Learning Approach

This report details the development of a machine learning (ML)-driven framework to coarse-grain inter-particle collision dynamics from high-fidelity Discrete Element Method (DEM) simulations to Particle-in-Cell (PIC) scales for gas-solid systems. Traditional PIC models, while computationally efficient, rely on empirical granular stress formulations that fail to capture the full complexity of collision physics, particularly the heterogeneity in particle dynamics. This study adopts a bottom-up approach, integrating insights from DEM simulations to improve the physical fidelity and interpretability of PIC-scale models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1991

High-performance multiprocessor architecture for a 3-D lattice gas model

The lattice gas method has recently emerged as a promising discrete particle simulation method in areas such as fluid dynamics. We present a very high-performance scalable multiprocessor architecture, called ALGE, proposed for the simulation of a realistic 3-D lattice gas model, Henon's 24-bit FCHC isometric model. Each of these VLSI processors is as powerful as a CRAY-2 for this application. ALGE is scalable in the sense that it achieves linear speedup for both fixed and increasing problem sizes with more processors. The core computation of a lattice gas model consists of many repetitions of two alternating phases: particle collision and propagation. Functional decomposition by symmetry group and virtual move are the respective keys to efficient implementation of collision and propagation.

Lee, F.↗
DOE OSTI2023

Coupled Experimental and High-Temperature Discrete-Element Method Modeling Studies of Aluminosilicate Particle Handling in Concentrated Solar Power Environments

Chemically inert, aluminosilicate based particles have been investigated as both a thermal transport and sensible energy storage medium for concentrated solar power facilities. These particles will experience a wide range of operating temperatures (300-1000 K) and handling conditions (dense to dilute falling particle curtains, dense granular flows, or dense structures), requiring specially-designed and optimized infrastructures. The relative influence of collisional and frictional interactions between particles varies based on temperature-dependent particulate properties and greatly impacts the bulk, granular flow behavior. These underlying physics are captured using discrete element method modeling tools. However, this modeling method is computationally expensive as each particle position and interaction is tracked during the simulation. These modeling methods are further complicated by introducing temperature-dependent particle properties, high-temperature radiative exchange, and directional irradiation sources experienced by granular flows in concentrated solar power environments. Coupled experimental and numerical studies of aluminosilicate particles in rotary kilns and dense particle curtains were performed for bulk temperatures up to 1073 K. The three particle types investigated included Carbobead HSP 30 /60, Carbobead CP 30/60, and Granusil 4030. Temperature, spatial, and velocity profile data were extracted from experimental runs using embedded K-type thermocouple probes and particle image velocimetry techniques. Experimental and numerical studies were compared using spatial temperature profiles, velocity fields, and shape profiles of the bulk, granular flows. Numerical models were developed using commercially available discrete element method modeling software, Aspherix®. Existing Aspherix® functionality was expanded by introducing coupled radiative exchange modeling tools. The laboratory-scale rotary kiln was developed to investigate the steady-state heat and mass transfer performance of aluminosilicate particles based on particle type, bulk handling temperature, and wall roughness. The rotational speed of the rotary kiln was varied to control the relative impact of collisional and frictional effects upon the granular flow behavior. Heat and mass transfer performance was categorized based on the Froude number and the observed flow regimes of slipping, rolling, cascading, and centrifuging. Coupled discrete element method modeling studies were used to evaluate the effects of temperature-dependent, particulate mechanical properties upon bulk flow behavior and upon the relative effects of radiative, advective, and/or conductive heat transfer. A high-temperature (< 1073 K) falling particle curtain was similarly fabricated to investigate the heat and mass transfer performance of aluminosilicate particles in particle handling situations dominated by inter-particle collisions. The impact of particle type, flow preheat temperatures (< 1073K), and bulk mass flow rates were investigated upon the particle curtain shape, temperature, and velocity profiles. Coupled discrete element method modeling studies were performed to evaluate the varying impact of temperature-dependent, particulate mechanical properties on the bulk flow behavior and the temperature profile of the particle curtain.

14 SOLAR ENERGY↗
NASA NTRS1989

Assessment of impingement effects in the isentropic core of a small satellite control thruster plume

Two computational techniques commonly employed in the calculation of rocket and thruster expansion plumes are assessed. These are the method of characteristics (MOC), which is derived from the continuum Euler equations, and the direct simulation Monte Carlo (DSMC) method, which adopts a discrete particle approach. These techniques vary both in the computational expense and in the accuracy and detail of the solutions that they provide, depending upon the regime of application. The assessment is made with reference to the plume expanding from a small monopropellant hydrazine thruster and concentrates on the isentropic core of the jet for the flow regime lying between the continuum and free molecular limits. It is found that the more numerically intensive DSMC method offers the better correspondence to the available experimental data. In addition, large differences in typical impingement effects such as drag force and heat transfer are found at the free molecular limit of the plume expansion for the two predictive techniques. It is concluded that accurate estimation of impingement potential may only be achieved through application of the discrete particle method.

Boyd, I. D.↗
DOE OSTI2022

Protein drift-diffusion dynamics and phase separation in curved cell membranes and dendritic spines: Hybrid discrete-continuum methods

We develop methods for investigating protein drift-diffusion dynamics in heterogeneous cell membranes and the roles played by geometry, diffusion, chemical kinetics, and phase separation. Our hybrid stochastic numerical methods combine discrete particle descriptions with continuum-level models for tracking the individual protein drift-diffusion dynamics when coupled to continuum fields. We show how our approaches can be used to investigate phenomena motivated by protein kinetics within dendritic spines. The spine geometry is hypothesized to play an important biological role regulating synaptic strength, protein kinetics, and self-assembly of clusters. We perform simulation studies for model spine geometries varying the neck size to investigate how phase-separation and protein organization is influenced by different shapes. We also show how our methods can be used to study the roles of geometry in reaction-diffusion systems including Turing instabilities. Furthermore, our methods provide general approaches for investigating protein kinetics and drift-diffusion dynamics within curved membrane structures.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

A projection method for particle resampling

Particle discretizations of partial differential equations are advantageous for high-dimensional kinetic models in phase-space due to their better scalability than continuum approaches with respect to dimension. Complex processes collectively referred to as particle noise hamper long time simulations with particle methods. One approach to address this problem is particle mesh adaptivity, or remapping, known as particle resampling and remeshing. Here, this work introduces a resampling method that projects particles to and from a (finite element) function space. The method is simple, using standard sparse linear algebra and finite element techniques, and it preserves all moments up to the order of a polynomial represented exactly by the continuum function space. It is distinguished from most other mesh-based methods in that new particle positions and number are decoupled from the mesh, allowing particle and continuum meshes to be adapted relatively independently. While this work is developed with structured particle and continuum phase-space grids on 1X + 1V Vlasov-Poisson models of Landau damping and two-stream instability, the method is well-suited to unstructured grids. Stable long time dynamics are demonstrated up to time T = 500. Reproducibility artifacts and data are publicly available.

Kinetic methods↗
DOE OSTI2021

Granular Flow

LIGGGHTS-INL is a capability-extended version of LIGGGHTS. LIGGGHTS is an Open Source Discrete Element Method Particle Simulation Software (https://www.cfdem.com/media/DEM/docu/Manual.html). LIGGGHTS stands for LAMMPS improved for general granular and granular heat transfer simulations. LAMMPS is a classical molecular dynamics simulator. It is widely used in the field of Molecular Dynamics (https://lammps.sandia.gov/doc/Manual.html). Thanks to physical and algorithmic analogies, LAMMPS is a very good platform for DEM simulations. LAMMPS offers a GRANULAR package to perform these kinds of simulations. LIGGGHTS aims to improve those capability with the goal to apply it to industrial applications.

Xia, Yidong↗
DOE OSTI2026

Onset of Fluidization in MP-PIC Simulations using MFIX-Exa

Fluidized bed reactors are used across a variety of industries, including for energy processes like pyrolysis that result in low-cost energy products. Design and scale-up of fluidized beds is de-risked by modeling and simulation, utilizing tools like NETL’s MFIX-Exa High-Performance Computing (HPC) code for reacting multiphase flow. This report summarizes an investigation into the breadth of problems to which MFIX-Exa may be applied, specifically with regard to low fluid velocities and the onset of fluidization. A simple fluidization study is conducted both experimentally and numerically for particles of interest, then reactor simulations are compared to cold flow experiments for uniform distributor plates. Approaches for modeling bubble caps are also presented.

discrete particle method↗