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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Towards data-driven constitutive modelling for granular materials via micromechanics-informed deep learning

The analytical description of path-dependent elastic-plastic responses of a granular system is highly complicated because of continuously evolving microstructures and strain localisation within the system undergoing deformation. This study offers an alternative to the current analytical paradigm by developing micromechanics-informed machine-learning based constitutive modelling approaches for granular materials. A set of critical variables associated with the constitutive behaviour of granular materials are identified through an incremental stress-strain relationship analysis. Depending on the strategy to exploit the priori micromechanical knowledge, three different training strategies are explored. The first model uses only the measurable external variables to make stress predictions; the second model utilises a directed graph to link all the external strain sequences and internal microstructural evolution variables into a single prediction model comprised of a series of sub-mappings, and the third model explicitly integrates the physically important non-temporal properties with external strain paths into training through an enhanced Gated Recurrent Unit (GRU). These three models show satisfactory agreement with unseen test specimens based on multi-directional loading cases. The basic features and potential applications of each model are explained. Lastly, the key factors for constitutive training and limitations of the current work are also discussed in detail.

36 MATERIALS SCIENCE↗
NASA NTRS2005

Status of Duct Liner Technology for Application to Aircraft Engine Nacelles

Grazing flows and high acoustic intensities impose unusual design requirements on acoustic liner treatments used in aircraft engine nacelles. Increased sound absorption efficiency (requiring increased accuracy of liner impedance specification) is particularly critical in the face of ever decreasing nacelle wall area available for liner treatments in modern, high-bypass ratio engines. This paper reviews the strategy developed at Langley Research Center for achieving a robust measurement technology that is crucial for validating impedance models for aircraft liners. Specifically, the paper describes the current status of computational and data acquisition technologies for reducing impedance in a flow duct. Comparisons of reduced impedances for a "validation liner" using 1980's and 2000's measurement technology are consistent, but show significant deviations (up to 0.5 c exclusive of liner anti-resonance region) from a first principles impedance prediction model as grazing flow centerline Mach numbers increase up to 0.5. The deviations, in part, are believed related to uncertainty in the choice of grazing flow parameters (e.g. cross-section averaged, core-flow averaged, or centerline Mach number?). Also, there may be an issue with incorporating the impedance discontinuities corresponding to the hard wall to liner interface (i.e. leading and trailing edge of test liner) within the discretized finite element model.

Parrott, Tony L.↗
NASA NTRS2021

Analysis of Nonlinear Shrinkage for the Bound Metal Deposition Manufacturing using Multi-scale Approach

We consider problem of nonlinear shrinkage of the metal part during bound metal deposition manufacturing on the ground and in zero-G. To analyze this problem we developed multi-scale physics-based approach that spans atomistic dynamics at the scale of nanoseconds and the full part shrinkage at the time scale of hours. Using this approach we estimated the key parameters of the problem including grain boundary width, coefficient of surface diffusion, initial redistribution of particles during debinding stage, micro-structure evolution from round particles to densely packed grains and corresponding change of the total and chemical free energy, and sintering stress. The introduced method was used to predict shrinkage at the level of two particles, filament cross-section, sub-model, and the whole green, brown, and metal parts. To further improve accuracy and reliability of the shrinkage predictions we propose concept of intelligent additive manufacturing of metal powders in space that combines the strengths of both physics-based and data-driven methods of analysis of AM.

bound metal deposition↗
DOE OSTI2023

Implementation of extrinsic cohesive zone model (ECZM) in 2D finite-discrete element method (FDEM) using node binding scheme

The combined finite-discrete element method (FDEM) has been widely used for rock fracturing simulations. Conventionally, FDEM is realized using the intrinsic cohesive zone model (ICZM); however, it has the drawback of artificial compliance and high computational expense. As a complement, the extrinsic cohesive zone model (ECZM) is seen to be realized in FDEM recently, whereas the node splitting scheme utilized is cumbersome. Here, within the framework of ICZM-based FDEM, we propose a node binding scheme to efficiently bind the pre-discretized finite elements and thus guarantee the continuum behavior of materials in the elastic stage. The yield surfaces, controlled by ECZM, are dynamically embedded by invoking the pre-inserted cohesive elements. The effectiveness and efficiency of the proposed approach are validated and tested by performing a suite of numerical experiments. Compared with ICZM-based FDEM, the proposed approach can correctly capture material deformation and reduce the computation cost. In contrast to the existing ECZM-based FDEM, the proposed approach can overcome the frequent and complex element topology updating. Finally, this work provides a novel perspective that fully inherits the advantages of both ICZM and ECZM, but circumvents their shortcomings, which guarantees a more efficient and effective simulation of brittle material evolution from continuum to discontinuum.

58 GEOSCIENCES↗
DOE OSTI2024

DEM simulation of the compression of crushable sand: does the initial particle shape matter?

Advances in DEM modeling, combined with high-resolution X-ray tomography, opened the way for computer models based on virtual replicas of the particles which preserve nearly all facets of their geometry. This leads to simulation advantages, but also high computational costs. Here we tackle a question stemming from this trend: how accurate should particle models be to ensure accuracy? We address this question for the case of the compression of crushable sand. LS-DEM was used to generate three models of Ottawa sand (exact replicas, ellipsoids, and spheres) from digital images of its grains. Compression-induced crushing was simulated for all sets by tracking evolving size and shape distribution. The results confirm that exact replicas provide the closest match of the measurements. However, intermediate degrees of rendering (e.g. ellipsoids preserving volume and aspect ratio of the real grains) led to satisfactory results only marginally different from those of exact replicas. In conclusion, these findings provide an example of the protocols that may be followed to identify the optimal degree of particle approximation which should be regarded as mandatory to achieve a conscious, sustainable use of computational resources.

58 GEOSCIENCES↗
DOE OSTI2018

Cohesion-Induced Stabilization in Stick-Slip Dynamics of Weakly Wet, Sheared Granular Fault Gouge

We use three-dimensional discrete element calculations to study stick-slip dynamics in a weakly wet granular layer designed to simulate fault gouge. The granular gouge is constituted by 8000 spherical particles with a poly-disperse size distribution. At very low liquid content, liquids impose cohesive and viscous forces on particles. Our simulations show that by increasing the liquid content, friction increases and granular layer shows higher recurrence time between slip events. We also observe that slip events exhibit larger friction drop and layer compaction in wet system compared to dry. We demonstrate that a small volume of liquid induces cohesive forces between wet particles that are responsible for an increase in coordination number leading to a more stable arrangement of particles. This stabilization is evidenced with two orders of magnitude lower particle kinetic energy in wet system during stick phase. Similar to previous experimental studies, we observe enhanced frictional strength for wet granular layers. In experiments, the physicochemical processes are believed to be the main reason for such behavior, we show however, that at low confining stresses the hydromechanical effects of induced cohesion are sufficient for observed behavior. Our simulations illuminate the role of particle interactions and demonstrate the conditions under which induced cohesion plays a significant role in fault zone processes, including slip initiation, weakening, and failure.

58 GEOSCIENCES↗
DOE OSTI2018

Discrete numerical simulations of torpedo anchor installation in granular soils

Torpedo anchors are a viable approach for mooring marine hydrokinetic (MHK) energy devices to the seafloor. These anchors can serve to maintain station and to provide the reaction force for an MHK device. The ability of the anchor to perform these duties is a strong function of its penetration depth during installation. This is a large-strain problem not amenable to typical continuum numerical approaches. In the current work, we propose that the discrete element method (DEM) is a more appropriate tool to investigate the shallow penetration of torpedo anchors in sands. The effects of anchor mass, impact velocity, and soil interparticle friction are considered in the DEM simulations. The relative maximum penetration depths for different penetration conditions are quantified and presented. Granular material response at the microscale during penetration are used to provide insight into system response. Energy dissipation in the assembly by both friction and collision at the particle scale are considered. Results show that anchor penetration increases approximately linearly with an increase in impact velocity or anchor weight. Penetration decreases with an increase in interparticle friction (i.e., soil strength). Observations of microscale behaviors and energy calculations are used to provide insight into overall system response.

58 GEOSCIENCES↗
DOE OSTI2021

The influence of preexisting host rock damage on fault network localization

The transition from stable to unstable fracture propagation occurs when fractures begin to interact and link. Thus, fracture network coalescence controls how rocks and engineered structures fail. Here, to constrain the factors that influence localization in shear zones under brittle conditions, we build discrete element method models with a rough fault embedded in a shear zone. We add varying numbers of diffuse, randomly-placed weaknesses to examine the influence of diffuse damage on fracture network localization. The number of weaknesses controls the localization behavior of the fault network and the final fault geometry. We quantify localization using the Gini coefficient of the fracture volume, which measures the nonuniformity in a population. Each model generally increases in localization toward failure. However, models with more diffuse damage experience delocalization phases that are superimposed on the overall trend of increasing localization. The observed link between delocalization and host rock damage may help explain the varying localization of low magnitude seismicity in southern California. Models with more diffuse damage produce more complex fault geometries comprised of several parallel strands of wing cracks. The propagation of these wing cracks reduces the shear stress acting on the model boundaries, indicating that this fracture development increases the mechanical efficiency of the system.

58 GEOSCIENCES↗
NASA NTRS2022

Trade Study of DEM Software

Review of Open-source Discrete Element Method Softwares: Three open-source solutions have been investigated for integration into OceanWATERS: Yade, EsysParticle, and ProjectChrono.

DEM↗
NASA NTRSDate not supplied

Monodisperse Single-Material Granular Tribocharging Modeling and Experimental Validation

Charge transfer between insulating grains has been a topic of interest for many years as this phenomenon is extremely important to many areas in industry. Of particular interest to NASA is the behavior of electrically insulative dust grains such as those found on the lunar surface. Whether poured from a scoop during sample collection activities, agitated inside a drum via mining robots, or fluidized by a gas plume, particle-particle interactions between similar granular materials will be widespread on the Moon. Designers hope to have a predictive model for how dust grains are charged, transported, and deposited on flight systems to better understand issues electrostatics may cause in future missions. While difference in work function is often seen as a driver for charge transfer between two materials, this is an incomplete representation for insulative particles. Other asymmetries in the system can also lead to charge being transferred: one such being the difference between static vs dynamic particles. Experiments performed under vacuum in the Electrostatics and Surface Physics Laboratory (ESPL) clearly show that particles having more contacts with other particles tend to charge positively and particles having fewer contacts tend to charge negatively, despite being made from the same material. During these experiments, monodisperse (population of single diameter) spherical particles were held in a reservoir which was then elevated at one end to slope into a Faraday cup used to measure bulk charge of the grains. Care was taken to limit the number of spheres bouncing out of the cup, to isolate the beads from the container itself so that only grain-to-grain interactions occurred, and to assure neutrality of the beads before the experiment began. A discrete element method (DEM) modeling package already incorporating many granular mechanics interactions was augmented to include the electrostatic behaviors seen in experiments. The model can reproduce the experimental results above by tracking the number of contacts between pairs of individual grains and applying a charge transfer condition related to this contact parameter. The experimental results as well as their incorporation into the improved modeling suite will be discussed. Future work includes adding additional asymmetries into the experiments/model such as polydisperse populations of spherical particles and irregularly shaped grains.

granular↗
NASA NTRS1980

On sound transmission into a stiffened cylindrical shell with rings and stringers treated as discrete elements

In the context of the transmission of airborne noise into an aircraft fuselage, a mathematical model is presented for the transmission of an oblique plane sound wave into a finite cylindrical shell stiffened by stringers and ring frames. The rings and stringers are modeled as discrete structural elements. The numerical case studied was typical of a narrow-bodied jet transport fuselage. The numerical results show that the ring-frequency dip in the transmission loss curve that is present for a monocoque shell is still present in the case of a stiffened shell. The ring frequency effect is a result of the cylindrical geometry of the shell. Below the ring frequency, stiffening does not appear to have any significant effect on transmission loss, but above the ring frequency, stiffeners can enhance the transmission loss of a cylindrical shell.

Koval, L. R.↗
DOE OSTI2018

A discrete element method representation of an anisotropic elastic continuum

A method for modeling cubically anisotropic elasticity within the discrete element method is presented. The discrete element method (DEM) is an approach originally intended for modeling granular materials (sand, soil, and powders); however, recent developments have usefully extended it to model stochastic mechanical processes in monolithic solids which, to date, have been assumed to be elastically isotropic. The method presented here for efficiently capturing cubic elasticity in DEM is an important prerequisite for further extending DEM to capture the influence of elastic anisotropy on the mechanical response of polycrystals, composites, etc. The system demonstrated here uses a directionally assigned stiffness in the bonds between adjacent elements and includes separate schemes for achieving anisotropy with Zener ratios greater and smaller than one. Finally, the model framework is presented along with an analysis of the accessible space of elastic properties that can be modeled and an artificial neural network interpolation scheme for mapping input parameters to model elastic behavior.

42 ENGINEERING↗
NASA NTRS2001

Analyses for Debonding of Stitched Composite Sandwich Structures Using Improved Constitutive Models

A fracture mechanics analysis based on strain energy release rates is used to study the effect of stitching in bonded sandwich beam configurations. Finite elements are used to model the configurations. The stitches were modeled as discrete nonlinear spring elements with a compliance determined by experiment. The constitutive models were developed using the results of flatwise tension tests from sandwich material rather than monolithic material. The analyses show that increasing stitch stiffness, stitch density and debond length decrease strain energy release rates for a fixed applied load.

Glaessgen, E. H.↗
NASA NTRS1989

On equivalence of discrete-discrete and continuum-discrete design sensitivity analysis

Developments in design sensitivity analysis (DSA) method have been made using two fundamentally different approaches as shown. In the first approach, a discretized structural finite element model is used to carry out DSA. There are three different methods in the discrete DSA approach: finite difference, semi-analytical, and analytical methods. The finite difference method is a popular one due to its simplicity, but a serious shortcoming of the method is the uncertainty in the choice of a perturbation step size of design variables. In the semi-analytical method, the derivatives of stiffness matrix is computed by finite differences, whereas in the analytical method, the derivatives are obtained analytically. For the shape design variable, computation of analytical derivative of stiffness matrix is quite costly. Because of this, the semi-analytical method is a popular choice in discrete shape DSA approach. However, recently, Barthelemy and Haftka presented that the semi-analytical method can have serious accuracy problems for shape design variables in structures modeled by beam, plate, truss, frame, and solid elements. They found that accuracy problems occur even for a simple cantilever beam. In the second approach, a continuum model of the structure is used to carry out DSA.

Choi, Kyung K.↗
NASA NTRSDate not supplied

Particle-Based Fiber Models of Woven Materials for Earth Entry Thermal Protection

Applications requiring materials with layer-to-layer strength, from basketry to thermal protection systems, use interlaced, three-dimensional woven materials. NASA is developing and deploying woven material heat shields for missions, including Artemis I (3D-MAT for compression pads) and potentially Mars Sample Return - Earth Entry System. These materials are complex, hierarchical and must protect from extreme environments and phenomena, such as deformation, impact and high-enthalpy heating. Woven material performance depends on microstructure, damage and weave geometry. Therefore, fiber-specific models are needed to simulate fiber contacts within the weave hierarchical geometry (fiber to tow to yarn to weave) and the inherent directionality of fibers. Explicit fiber models can simulate how weave microstructure evolution affects thermal and mechanical properties. We parameterize a discrete element bonded particle models (DEM-BPM) of fibers to capture thermal and mechanical behavior within and between fibers, with bonded and contact forces, respectively. We study the proportion of heat transfer and stress via the contact network, fiber bonds and the weave geometry, for example, with respect to yarn warp-weft identity (whether it interlaces weave layers). Our results demonstrate the importance of explicit fiber modeling for connecting microstructure with thermal and mechanical properties.

woven material↗
NASA NTRSDate not supplied

Particle-Based Fiber Models of Woven Materials for Earth Entry Thermal Protection

Applications requiring materials with layer-to-layer strength, from basketry to thermal protection systems, use interlaced, three-dimensional woven materials. NASA is developing and deploying woven material heat shields for missions, including Artemis I (3D-MAT for compression pads) and potentially Mars Sample Return - Earth Entry System. These materials are complex, hierarchical and must protect from extreme environments and phenomena, such as deformation, impact and high-enthalpy heating. Woven material performance depends on microstructure, damage and weave geometry. Therefore, fiber-specific models are needed to simulate fiber contacts within the weave hierarchical geometry (fiber to tow to yarn to weave) and the inherent directionality of fibers. Explicit fiber models can simulate how weave microstructure evolution affects thermal and mechanical properties. We parameterize a discrete element bonded particle models (DEM-BPM) of fibers to capture thermal and mechanical behavior within and between fibers, with bonded and contact forces, respectively. We study the proportion of heat transfer and stress via the contact network, fiber bonds and the weave geometry, for example, with respect to yarn warp-weft identity (whether it interlaces weave layers). Our results demonstrate the importance of explicit fiber modeling for connecting microstructure with thermal and mechanical properties.

woven material↗
NASA NTRS1986

Calculations of axisymmetric vortex sheet roll-up using a panel and a filament model

A method for calculating the self-induced motion of a vortex sheet using discrete vortex elements is presented. Vortex panels and vortex filaments are used to simulate two-dimensional and axisymmetric vortex sheet roll-up. A straight forward application using vortex elements to simulate the motion of a disk of vorticity with an elliptic circulation distribution yields unsatisfactroy results where the vortex elements move in a chaotic manner. The difficulty is assumed to be due to the inability of a finite number of discrete vortex elements to model the singularity at the sheet edge and due to large velocity calculation errors which result from uneven sheet stretching. A model of the inner portion of the spiral is introduced to eliminate the difficulty with the sheet edge singularity. The model replaces the outermost portion of the sheet with a single vortex of equivalent circulation and a number of higher order terms which account for the asymmetry of the spiral. The resulting discrete vortex model is applied to both two-dimensional and axisymmetric sheets. The two-dimensional roll-up is compared to the solution for a semi-infinite sheet with good results.

Kantelis, J. P.↗