Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “discrete element modeling”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Modeling bidirectional reflectance of forests and woodlands using Boolean models and geometric optics

Geometric-optical discrete-element mathematical models for forest canopies have been developed using the Boolean logic and models of Serra. The geometric-optical approach is considered to be particularly well suited to describing the bidirectional reflectance of forest woodland canopies, where the concentration of leaf material within crowns and the resulting between-tree gaps make plane-parallel, radiative-transfer models inappropriate. The approach leads to invertible formulations, in which the spatial and directional variance provides the means for remote estimation of tree crown size, shape, and total cover from remotedly sensed imagery.

Strahler, Alan H.↗
NASA NTRS1982

Approximating a finite element model of a hybrid deployable truss by a discrete mass model

The NASTRAN finite element model of the MSFC/hybrid deployable truss is approximated by a model consisting of four bodies serially connected by a spring hinge suspension. The truss is decomposed into four serially connected modules determining the total mass and rotational inertias for a rigid body with its mass center on the axis of the truss for each module. A five-body cantilever truss model is developed along with an undamped form of this model. Spring coefficient values are generated to match eigenvalues of the latter model with those of the NASTRAN model. Damping coefficients are then determined that provide n percent of critical damping for the model.

Chichester, F. D.↗
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 NTRS2012

Elastic Model Transitions Using Quadratic Inequality Constrained Least Squares

A technique is presented for initializing multiple discrete finite element model (FEM) mode sets for certain types of flight dynamics formulations that rely on superposition of orthogonal modes for modeling the elastic response. Such approaches are commonly used for modeling launch vehicle dynamics, and challenges arise due to the rapidly time-varying nature of the rigid-body and elastic characteristics. By way of an energy argument, a quadratic inequality constrained least squares (LSQI) algorithm is employed to e ect a smooth transition from one set of FEM eigenvectors to another with no requirement that the models be of similar dimension or that the eigenvectors be correlated in any particular way. The physically unrealistic and controversial method of eigenvector interpolation is completely avoided, and the discrete solution approximates that of the continuously varying system. The real-time computational burden is shown to be negligible due to convenient features of the solution method. Simulation results are presented, and applications to staging and other discontinuous mass changes are discussed

Orr, Jeb S.↗
DOE OSTI2023

Application of a Distributed Element Roughness Model to Additively Manufactured Internal Cooling Channels

Abstract Design for cooling effectiveness in turbine blades relies on accurate models for dynamic losses and heat transfer of internal cooling passages. Metal additive manufacturing (AM) has expanded the design space for these configurations, but can give rise to large-scale roughness features. The range of roughness length scales in these systems makes morphology resolved computational fluid dynamics (CFD) impractical. However, volumetric roughness models can be leveraged, as they have computational costs orders of magnitude lower. In this work, a discrete element roughness model (DERM), based on the double-averaged Navier–Stokes equations, is presented and applied to additively manufactured rough channels, representative of gas turbine blade cooling passages. Unique to this formulation of DERM is a generalized sheltering-based treatment of drag, a two-layer model for spatially averaged Reynolds stresses, and explicit treatment of dispersion. Six different AM rough surface channel configurations are studied, with roughness trough to peak sizes ranging from 15% to 60% nominal channel passage half-width, and the roughness Reynolds number ranges from Rek = 60 to 300. DERM predictions for spatially and temporally averaged mean flow quantities are compared to previously reported direct numerical simulation results. Good agreement in the mean velocity profiles, stress balances, and drag partitions are observed. While DERM models are typically calibrated to specific deterministic roughness morphologies at comparatively small roughness Reynolds numbers, the present more generalized DERM formulation has wider applicability. Here, it is demonstrated that the model can accommodate random roughness of large scale, typical of AM.

Engineering↗
NASA NTRS1983

An analytical investigation of shape control of large space structures by applied temperatures

An analytical procedure for the static shape control of flexible space structures subjected to thermal distortions is developed which is based on prescribing temperatures in control elements having much higher coefficients of thermal expansion than the main structure. The temperatures at the control elements are defined so as to minimize the overall thermal distortion of the structure from its ideal shape, and a matrix equation is obtained which can be solved for the set of optimum control temperatures. A formulation of the procedure for continuous structures governed by differential equations and a formulation for discrete (finite element modeled) structures governed by matrix equations are presented. The equations from the continuous formulation are employed for the shape control of a simple beam distorted by nonuniform heating, and the discrete formulation is applied in a general purpose finite-element structural analysis computer program for the shape control of a 750 m radiometer antenna reflector dish subjected to orbital heating. A reduction in thermal distortion by a factor of nearly 50 was obtained with the use of only seven control elements. Results for four different sets of control locations for the antenna are presented in which reductions in distortion of up to a factor of four were obtained.

Haftka, R. T.↗
DOE OSTI2025

Dense granular flows with MFIX-Exa

This report extends the linear spring dashpot collision model the discrete element method available in MFIX-Exa to include static a static tangential friction force. Additionally, two rolling friction models frequently used in the literature are also implemented. The governing equations are provided with an emphasis on the new terms. The new model is validated by comparison to existing experimental data of single particle oblique collisions. The model is then tested on three dense granular flow problems: the formation of static piles, the discharge from a flat-bottom hopper and the self-induced granular Rayleigh-Taylor instability.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1985

Modeling of convection phenomena in Bridgman-Stockbarger crystal growth

Thermal convection phenomena in a vertically oriented Bridgman-Stockbarger apparatus were modeled by computer simulations for different gravity conditions, ranging from earth conditions to extremely low gravity, approximate space conditions. The modeling results were obtained by the application of a state-of-the art, transient, multi-dimensional, completely densimetrically coupled, discrete-element computational model which was specifically developed for the simulation of flow, temperature, and species concentration conditions in two-phase (solid-liquid) systems. The computational model was applied to the simulation of the flow and the thermal conditions associated with the convection phenomena in a modified Germanium-Silicon charge enclosed in a stationary fused-silica ampoule. The results clearly indicated that the gravitational field strength influences the characteristics of the coherent vortical flow patterns, interface shape and position, maximum melt velocity, and interfacial normal temperature gradient.

Carlson, F. M.↗
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 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 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↗