Search NASA⌕ Search

SEARCH · Search NASA

Results for “Discrete element models”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

A Review of Computational Models for the Flow of Milled Biomass Part I: Discrete-Particle Models

Biomass is a renewable and sustainable energy resource. Current design of biomass handling and feeding equipment leverage both experiments and numerical modeling. This paper reviews the state-of-the-art discrete element methods (DEM) for the flow of milled biomass (Part I), accompanied by a comprehensive review on continuum-based computational models (Part II). The present review on DEM is primarily focused on the features and suitability of various particle shape models for different types of milled biomass because particle shape is the predominant attribute controlling the flow behavior of complex-shaped granular material. The general strengths and weaknesses in the applicability of those models for the milled biomass modeling are summarized. In particular, comments are provided to balance the numerical model capabilities and the computational cost for the development of DEM models. To our best knowledge, this is the first-of-its-kind review on DEM specifically for biomass. Our study indicates that the current DEM models require further development, calibration, and validation based on a deep understanding of biomass particle contact mechanics and experimental data support before they can be reliably used for predictive simulations in handling and feeding systems.

42 ENGINEERING↗
DOE OSTI2020

An Integrated Multiscale Experimental-Numerical Analysis on Reconsolidation of Salt-Clay Mixture for Disposal of Heat-Generating Waste (Final NEUP Technical Report)

The overall purpose of this research is to improve understanding of THMC coupling effect on the reconsolidation of granular (or crushed) salt-clay mixture used for seal systems of shafts and drifts in salt repositories. This proposed work is partially motivated by the recent work on the Waste Isolation Pilot Plant (WIPP) that shows the promising sealing capability of clay-salt mixture compared to crushed salt. In particular, primary emphasis is to develop a fully integrated multiscale experiment-numerical study to determine and explain what leads to the superior sealing ability of the clay-salt mixture. These research activities are designed to seek further understanding of (1) why clay additives may enhance the fluid trapping and (2) whether this flow barrier effect may prevail under different combinations of temperature, confining pressure, deviatoric stress and other foreseeable environmental factors. If successful, this enhanced flow trapping ability of the seal provides significant improvement to the seal and repository performance and therefore make the repository safer in the long-term. The experiment component includes microstructural investigation and macroscopic tests on a reconsolidated salt-clay mixture. In the microstructural study, the goal is to (1) characterize microscopic distributions of distinct phases (e.g., clay, salt crystal boundaries, trapped brine, and pore) to examine the connectivity of the pore network inside the salt-clay mixture with different amounts of clay additive and moisture content and (2) analyze multiscale imaging data to reconstruct the polycrystalline microstructures for numerical simulations. Meanwhile, macroscopic tests are performed to analyze how clay alters the failure/creep mechanisms in the salt-clay mixture. Microscopic and macroscopic experimental observations will both be used to calibrate and validate a multiscale model that explicitly simulates the capillary and multiphase flow in the connected pores and the deformation due to the presence of intra-crystalline brine at the pore scale via a new polyhedral discrete element–lattice Boltzmann method (DEM-LBM) coupling model. The pore-scale simulations are homogenized via an upscaling procedure that converts pore-scale information (e.g. force exerted on grain boundary, sliding, pressure-solution) to continuum measures (e.g. Cauchy stress, Darcy’s flow) at each integration point in the macroscopic multiphase TMHC model. This multiscale scheme will allow coupling be- tween high-fidelity simulations of brine-salt-clay interaction and the macroscopic TMHC model. The multiscale model helps the understanding of how the trapped brine inclusion affects the pressure-solution mechanism with the presence of clay and moisture. This work brings new insight into the sealing capacity of salt-clay mixture under elevated temperature over a long period of time - a key to evaluating the potential of salt-clay mixture usage for salt repositories.

42 ENGINEERING↗
DOE OSTI2022

A coupled DEM-IMB-LBM model for simulating methane hydrate exploitation involving particle dissolution

The coupled discrete element and lattice Boltzmann method using an immersed moving boundary scheme was extended to simulate methane hydrate exploitation involving mass transport and particle dissolution. In this coupled DEM-IMB-LBM model, a new Dirichlet-type thermal boundary condition is extended to simulate moving curved boundaries with constant concentration. A novel periodic boundary including an efficient searching algorithm for particle contact is proposed to reduce the computational cost and boundary effect. So, this model is validated by two numerical examples: a circular particle with concentration convection-diffusion moving in a horizontal channel and mass transport from a cylinder particle in a simple shear flow. The numerical results obtained from the proposed model agree well with previous studies. To further demonstrate the capacity of the proposed model, simulations of methane hydrate exploitation including two formations in marine sediments are carried out. The numerical results indicate that the coupled DEM-IMB-LBM is not only capable of simulating the dissolution of hydrate particles at the grain level, but also recover the sand erosion and migration process in a fundamental perspective during the methane hydrate exploitation process.

42 ENGINEERING↗
DOE OSTI2021

The mechanics of brittle granular materials with coevolving grain size and shape

The influence of particle shape on the mechanics of sand is widely recognized, especially in mineral processing and geomechanics. However, most existing continuum theories for engineering applications do not encompass the morphology of the grains and its evolution during comminution. Similarly, the relatively few engineering models accounting for grain-scale processes tend to idealize particles as spheres, with their diameters considered as the primary and sole geometric descriptor. This paper inspires a new generation of constitutive laws for crushable granular continua with arbitrary, yet evolving, particle morphology. We explore the idea of introducing multiple grain shape descriptors into Continuum Breakage Mechanics (CBM), a theory originally designed to track changes in particle size distributions during confined comminution. We incorporate the influence of these descriptors on the elastic strain energy potential and treat them as dissipative state variables. In analogy with the original CBM, and in light of evidence from extreme fragmentation in nature, the evolution of the additional shape descriptors is postulated to converge towards an attractor. Comparisons with laboratory experiments, discrete element analyses and particle-scale fracture models illustrate the encouraging performance of the theory. The theory provides insights into the feedback among particle shape, compressive yielding and inelastic deformation in crushable granular continua. These results inspire new questions that should guide future research into crushable granular systems using particle-scale imaging and computations.

36 MATERIALS SCIENCE↗
NASA NTRSDate not supplied

Understanding Variability in Mechanical Behavior of Adhesively Bonded Sandwich Joints due to Stochastic Material Damage

A mesoscale progressive damage analysis was applied to the analysis of recently tested adhesively bonded sandwich panels, in which ten nominally identical tension loaded panels failed in three different failure modes, with peakloads ranging from 13.0 –16.1 kips. A parametric finite element model for progressive damage analysis of an adhesively bonded joint (ABJ) considering matrix and fiber damage, and delamination cracks in unidirectional and fabric composite sections of the joint, as well as core crushing, and adhesive fracture was developed. The model was developed using Abaqus/Explicit and utilized continuum damage mechanics material models for intralaminar damage in the facesheet and doubler plies, cohesive elements for delamination and adhesive debonding, and discretized cells considering plasticity to model crushing of the aluminum honeycomb core. Analysis predictions obtained with a pristine model agreed well with the average experimental peak load, strain in the joint, and post-test damage states. Then, various methods to incorporate stochastic material variability into ABJ model were evaluated. These methods included simulation of pre-existing matrix damage, pre-existing adhesive porosity, and uniform modification of adhesive properties. These material and manufacturing defects changed the damage mode predictions in some cases. The analysis of stochastic strength due to the considered defects revealed unintuitive damage mode interaction in the ABJ, explaining some structure-property relationships observed experimentally.

Richard Larson↗
NASA NTRS1980

Prediction of the interior noise levels of high-speed propeller-driven aircraft

The theoretical basis for an analytical model developed to predict the interior noise levels of high-speed propeller-driven airplanes is presented. Particular emphasis is given to modeling the transmission of discrete tones through a fuselage element into a cavity, estimates for the mean and standard deviation of the acoustic power flow, the coupling between a non-homogeneous excitation and the fuselage vibration response, and the prediction of maximum interior noise levels. The model allows for convenient examination of the various roles of the excitation and fuselage structural characteristics on the fuselage vibration response and the interior noise levels, as is required for the design of model or prototype noise control validation tests.

Rennison, D. C.↗
DOE OSTI2025

Multiscale Modeling Framework Using Element‐Based Galerkin Methods for Moist Atmospheric Limited‐Area Simulations

This paper presents a multiscale modeling framework (MMF) to model moist atmospheric limited-area weather. The MMF resolves large-scale convection using a coarse grid while simultaneously resolving local features through numerous fine local grids and coupling them seamlessly. Both large- and small-scale processes are modeled using the compressible Navier-Stokes equations within the Nonhydrostatic Unified Model of the Atmosphere (NUMA), and are discretized using a continuous element-based Galerkin method (spectral elements) with high-order basis functions. Consequently, the large-scale and small-scale models share the same dynamical core but have the flexibility to be adjusted individually. The proposed MMF method is tested in 2D and 3D idealized limited-area weather problems involving storm clouds produced by squall line and supercell simulations. Numerical results from the MMF showed enhanced representation of cloud processes compared to the coarse model.

Kang, Soonpil [Naval Postgraduate School, Monterey↗
NASA NTRS1995

Micromechanical models for textile composites

Numerical and analytical micromechanical models are presented to predict the thermoelastic behavior of a textile composite. In the numerical model, the unit-cell is discretized with finite elements, and periodic boundary conditions are imposed between opposite faces of the unit-cell. For a thin textile composite, stress gradients effects through the thickness are demonstrated. The consequent difference in the stiffness and strength behavior of thick and thin composites are discussed. The numerical model is implemented to predict 3-D thermo-elastic constants for a thick textile composite, and the plate thermo-mechanical properties for a thin textile composite. The numerical model is extended to compute the thermal residual microstresses due to processing to predict the composite failure envelopes. An analytical model - Selective Averaging Method (SAM) - is proposed, which is based on a judicious combination of stiffness and compliance averaging to estimate the 3-D elastic constants. Both the models are tested and verified for several examples by comparing the stiffness properties with elasticity solutions and available results.

Sankar, Bhavani V.↗
NASA NTRS1982

Finite element thermal-structural modeling of orbiting truss structures

A description of an integrated finite element (FE) thermal-structural approach for accurate and efficient modeling of large space structures is presented. A geometric model with a common discretization for all analyses is employed. It uses improved thermal elements and the results from the thermal analysis directly in the structural analysis without any intervening data processing. The differences between the conventional FE approach as implemented in large programs and an integrated FE approach currently under development are described. Considerations for thermal modeling of truss members is discussed and three thermal truss finite elements are presented. The performance of these elements was evaluated for typical truss members neglecting joint effects. A simple truss with metallic joints and composite members was studied to evaluate the effectiveness of the approach for realistic truss designs. A study of the effects of aluminum joints on the thermal deformations of a simple, plane truss with composite members showed that joint effects may be significant. Further study is needed to assess the role of joint effects on the deformation of large trusses.

Thornton, E. A.↗
DOE OSTI2017

Simulations of the effects of proppant placement on the conductivity and mechanical stability of hydraulic fractures

Here, we generate a wide range of models of proppant-packed fractures using discrete element simulations, and measure fracture conductivity using finite element flow simulations. This allows for a controlled computational study of proppant structure and its relationship to fracture conductivity and stress in the proppant pack. For homogeneous multi-layered packings, we observe the expected increase in fracture conductivity with increasing fracture aperture, while the stress on the proppant pack remains nearly constant. This is consistent with the expected behavior in conventional proppant-packed fractures, but the present work offers a novel quantitative analysis with an explicit geometric representation of the proppant particles. In single-layered packings (i.e. proppant monolayers), there is a drastic increase in fracture conductivity as the proppant volume fraction decreases and open flow channels form. However, this also corresponds to a sharp increase in the mechanical stress on the proppant pack, as measured by the maximum normal stress relative to the side crushing strength of typical proppant particles. We also generate a variety of computational geometries that resemble highly heterogeneous proppant packings hypothesized to form during channel fracturing. In some cases, these heterogeneous packings show drastic improvements in conductivity with only moderate increase in the stress on the proppant particles, suggesting that in certain applications these structures are indeed optimal. We also compare our computer-generated structures to micro computed tomography imaging of a manually fractured laboratory-scale shale specimen, and find reasonable agreement in the geometric characteristics.

58 GEOSCIENCES↗
DOE OSTI2019

From Stress Chains to Acoustic Emission

A numerical scheme using the combined finite-discrete element method is employed to study a model of an earthquake system comprising a granular layer embedded in a formation. When the formation is driven so as to shear the granular layer, a system of stress chains emerges. The stress chains endow the layer with resistance to shear and on failure launch broadcasts into the formation. Furthermore, these broadcasts, received as acoustic emission, provide a remote monitor of the state of the granular layer of the earthquake system.

58 GEOSCIENCES↗
DOE OSTI2024

High-Order Hybrid RANS-LES Study of NACA0012 Wing Sections

We develop hybrid RANS-LES strategies within Nek5000 for application to airfoil sections at small flight configurations. We present a validation and verification study of $k \space – \space \tau$ SST applied to a NACA 0012 wing section in a pure RANS and in a hybrid RANS-LES setup. The study shows good corroboration with existing experimental and numerical datasets. We also analyze some of the observed discrepancies with the experiments by evaluating the side wall “blocking” effect. We demonstrate that for the hybrid turbulence modeling approach a high-order spectral- element discretization converges faster (i.e., with less resolution) than a representative low-order finite-volume-based approach.

42 ENGINEERING↗
NASA NTRS1982

Combined discrete network. Continuous control modelling of operator behavior

The situation in which an operator is faced with a continuous control task plus one or more discrete information processing tasks is discussed. A modelling approach which has the capability of realistically representing both types of tasks and the resulting interactions is sought. The modelling approach utilizes discrete network models for the cognitive tasks and elements of an open-loop - closed-loop control representation for the continuous task. the approach is demonstrated through its application to a simulated Digital Avionics Information System in which subjects were required to perform retrieval and processing tasks as well as flight control. This model differs from conventional model in that system status sampling is not necessarily continuous or periodic. Rather, the pilot is assumed to read system status displays only as time permits and operate in a so-called open-loop preprogrammed fashion between sampling.

Miller, R. A.↗
DOE OSTI2023

Modelling and simulation of brinicle formation

Below the Arctic sea ice, under the right conditions, a flux of icy brine flows down into the sea. The icy brine has a much lower fusion point and is denser than normal seawater. As a result, it sinks while freezing everything around it, forming an ice channel called a brinicle (also known as ice stalactite). In this paper, we develop a mathematical model for this phenomenon, assuming cylindrical symmetry. The fluid is considered to be viscous and quasi-stationary. The heat and salt transport are weakly coupled to the fluid motion and are modelled with the corresponding conservation equations, accounting for diffusive and convective effects. Finite-element discretization is employed to solve the coupled system of partial differential equations. We find that the model can capture the general behaviour of the physical system and generate brinicle-like structures while also recovering dendrite composition, which is a physically expected feature aligned with previous experimental results. This represents, to our knowledge, the first complete model proposed that captures the global structure of the physical phenomenon even though it has some discrepancies, such as brine accumulation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

GPU accelerated Discrete Element Method (DEM) molecular dynamics for conservative, faceted particle simulations

Faceted shapes, such as polyhedra, are commonly found in systems of nanoscale, colloidal, and granular particles. Many interesting physical phenomena, like crystal nucleation and growth, vacancy motion, and glassy dynamics are challenging to model in these systems because they require detailed dynamical information at the individual particle level. Within the granular materials community the Discrete Element Method has been used extensively to model systems of anisotropic particles under gravity, with friction. We provide an implementation of this method intended for simulation of hard, faceted nanoparticles, with a conservative Weeks–Chandler–Andersen (WCA) interparticle potential, coupled to a thermodynamic ensemble. This method is a natural extension of classical molecular dynamics and enables rigorous thermodynamic calculations for faceted particles.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Investigating Powder Spreading Dynamics in Additive Manufacturing Processes by In-situ High-speed X-ray Imaging

The quality of the powder bed is known to be one of the main factors that influence the quality of the part being manufactured by powderbed- based additive manufacturing processes. For example, in the laser powder bed fusion process, it has been demonstrated that improving the powder bed density and homogeneity improves the final quality of the part. Understanding powder spreading dynamics in the powder-bed-based additive manufacturing processes is critical for designing spreading structures and spreading parameters, identifying defect formation mechanisms in the powder bed, and guiding the design of optimized feedstock powders. Extensive efforts have been devoted to understand the fundamental mechanism of the powder flow and spreading process. Various powder flow testers were used to study the flow behavior of powders used in additive manufacturing. However, the powder flow conditions in flow meters are different from those in real additive manufacturing conditions. Models have also been developed to simulate and study the spreading process. The discrete element method (DEM) is currently the main tool for computationally simulating the spreading process. However, these models are based on many physical assumptions. Visible light cameras have been used recently to conduct in-situ study of powder spreading in additive manufacturing. Due to the limited resolution and penetration depth of visible light, it is extremely difficulty to study the particle-scale powder spreading dynamics. Here, to overcome the challenges of studying powder spreading dynamics under real additive manufacturing conditions, we developed an experimental approach to study the particle-scale powder spreading process with high spatial and temporal resolution using the high-speed high-energy X-ray imaging facility at beamline 32-ID-B of the APS. Figure 1 shows an overview of the research setup. It is composed of a vacuum chamber (building environment), an automatic powder spreading system, a laser system, an X-ray source (X-ray outlet), and an X-ray detection/imaging system.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1984

The numerical solution of the Helmholtz equation for wave propagation problems in underwater acoustics

The Helmholtz Equation (-delta-K(2)n(2))u=0 with a variable index of refraction, n, and a suitable radiation condition at infinity serves as a model for a wide variety of wave propagation problems. A numerical algorithm was developed and a computer code implemented that can effectively solve this equation in the intermediate frequency range. The equation is discretized using the finite element method, thus allowing for the modeling of complicated geometrices (including interfaces) and complicated boundary conditions. A global radiation boundary condition is imposed at the far field boundary that is exact for an arbitrary number of propagating modes. The resulting large, non-selfadjoint system of linear equations with indefinite symmetric part is solved using the preconditioned conjugate gradient method applied to the normal equations. A new preconditioner is developed based on the multigrid method. This preconditioner is vectorizable and is extremely effective over a wide range of frequencies provided the number of grid levels is reduced for large frequencies. A heuristic argument is given that indicates the superior convergence properties of this preconditioner.

Bayliss, A.↗
NASA NTRS1986

A finite element overlay technique for modeling pinned composite joints

A finite element technique using overlaid plane stress elements is presented for modeling a pinned composite joint. The technique allows for discrete modeling of the pin and web regions of the joint in a two-dimensional (2-D) finite element analysis. Thus, predictive capability is substantially increased without the added complexity of a three-dimensional (3-D) analysis. The overlay technique requires a thorough evaluation of the model constraints between the joint components. Bearing panels are used to account for a reduced compressive modulus under bearing loads. Application of the plane stress overlay technique to the Filament Wound Case (FWC) Program is discussed. The experimental approach to determine bearing panel moduli using a double lap shear test is detailed, and finite element results are correlated with full-scale FWC joint test results.

Colvin, G. E.↗