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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

An Adaptive Timestepping Algorithm for Particle Time Integration in Coupled CFD-DEM Simulations

Coupled Computational Fluid Dynamics(CFD) - Discrete Element Method (DEM) models provide an accurate description of multiphase physical systems where a granular phase exists in an underlying continuous medium. The time integration of the granular phase in these simulations are typically handled using an explicit time integration scheme with a global time step that is much smaller than fluid time scales, in order to resolve inter-particle collisions. An adapative time integration technique that avoids this global constraint for the granular phase is presented in this work. The algorithm involves the use of a local time stepping approach to resolve collisional time scales for only a subset of particles that are in close proximity to potential collision partners, thereby resulting in substantial reduction of computational cost. This approach is observed to be 2-3X faster than traditional explicit methods for problems that involve both dense and dilute regions, while maintaining the same level of accuracy.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Islet: interpolation semi-Lagrangian element-based transport

Abstract. Advection of trace species, or tracers, also called tracer transport, in models of the atmosphere and other physical domains is an important and potentially computationally expensive part of a model's dynamical core. Semi-Lagrangian (SL) advection methods are efficient because they permit a time step much larger than the advective stability limit for explicit Eulerian methods without requiring the solution of a globally coupled system of equations as implicit Eulerian methods do. Thus, to reduce the computational expense of tracer transport, dynamical cores often use SL methods to advect tracers. The class of interpolation semi-Lagrangian (ISL) methods contains potentially extremely efficient SL methods. We describe a finite-element ISL transport method that we call the interpolation semi-Lagrangian element-based transport (Islet) method, such as for use with atmosphere models discretized using the spectral element method. The Islet method uses three grids that share an element grid: a dynamics grid supporting, for example, the Gauss–Legendre–Lobatto basis of degree three; a physics parameterizations grid with a configurable number of finite-volume subcells per element; and a tracer grid supporting use of Islet bases with particular basis again configurable. This method provides extremely accurate tracer transport and excellent diagnostic values in a number of verification problems.

58 GEOSCIENCES↗
NASA NTRS1998

Debond Analyses for Stitched Composite Structures

The effect of stitching on mode I and mode II strain energy release rates for debond configurations is studied using an analysis based on plate finite elements and the virtual crack closure technique. The stitches were modeled as discrete nonlinear fastener elements with a compliance determined by experiment. The axial and shear behavior of the stitches was considered with both the compliances and failure loads assumed to be independent. The mode I strain energy release rate, G(sub I), was shown to decrease once the debond had grown beyond the first row of stitches and was reduced to zero for long debonds, however, the mode II strain energy release rate, G(sub II), continued to be of significant magnitude over the range of debond lengths considered.

Glaessgen, E. H.↗
DOE OSTI2020

A numerical study of particle jetting in a dense particle bed driven by an air-blast

In this work, the particle jetting behavior in a blast-driven dense particle bed is studied at early times. Four-way coupled Euler–Lagrange simulations are performed using a high-order discontinuous Galerkin spectral element solver coupled with a high-order Lagrangian particle solver, wherein the inter-particle collisions are resolved using a discrete element method collision model. Following the experiments of Rodriguez et al. [“Formation of particle jetting in a cylindrical shock tube,” Shock Waves 23(6), 619–634 (2013)] and the simulations of Osnes et al. [“Numerical simulation of particle jet formation induced by shock wave acceleration in a Hele-Shaw cell,” Shock Waves 28(3), 451–461 (2018)], the simulations are performed in a quasi-two-dimensional cylindrical geometry (Hele-Shaw cell). Parametric studies are carried out to assess the impact of the coefficient of restitution and the strength of the incident shock on the particle jetting behavior. The deposition of vorticity through a multiphase (gas–particle) analog of Richtmyer–Meshkov instability is observed to play a crucial role in channeling the particles into well-defined jets at the outer edge of the particle bed. This is confirmed by the presence of vortex pairs around the outer jets. Furthermore, the effect of the relaxation of the relative velocity between the two phases on the vorticity generation is explored by analyzing the correlation between the radial velocity of particles and the radial velocity of the gas at the particle location.

42 ENGINEERING↗
DOE OSTI2020

An error-controlled adaptive time-stepping method for particle advancement in coupled CFD-DEM simulations

Coupled Computational-Fluid-Dynamics (CFD) and Discrete-Element-Method (DEM) models provide an accurate description of multiphase physical systems where a solid granular particle phase exists in an underlying gaseous continuous medium. The time integration of the granular phase in these simulations is typically handled using an explicit scheme with a constant time-step among all particles that is invariant in time to resolve inter-particle collisions. A locally third-order accurate adaptive time integration technique for particles that employs an embedded locally second-order scheme for error determination is presented in this work. The particle time-step size is dynamically adapted based on solution error, thus leading to significant savings in computational time. The efficacy of our scheme is quantified using four test cases of varying complexity (binary collision, homogeneous cooling system, fluidized bed and hopper discharge). The adaptive time-stepping method exhibits improved performance (~ 2–3 times in most of the cases studied) compared to three commonly used non-adaptive time-step methods (first-order Euler-explicit, second-order Adams-Bashforth and third-order Runge-Kutta schemes), while maintaining the same level of accuracy and parallel scalability.

42 ENGINEERING↗
DOE OSTI2020

Numerical analyses of high temperature dense, granular flows coupled to high temperature flow property measurements for solar thermal energy storage

High temperature particle flow properties necessary to predict granular flow behavior for solar thermal energy storage applications were measured and calculated for Carbobead CP 30/60 up to 800 °C. Here, the measured properties included elastic and shear moduli, particle-particle coefficients of static sliding and rolling friction, and particle-particle coefficients of restitution. Poisson’s ratio was calculated with elastic and shear moduli. The flow properties were used as inputs for a numerical model using the discrete element method to examine granular flows along an inclined plane at high temperature. The flow behavior was strongly influenced by the coefficients of static friction, which impacted the particle residence time, shear effects from the side walls, and particle flow mass flux. An 8.7%, 15.6%, and 8.5% increase and 37.9% decrease in steady state mass flow rate was observed for 200 °C, 400 °C, 600 °C, and 800 °C, respectively, when compared to room temperature simulations. A 52%, 59%, and 33% decrease in the time to reach steady state was observed for 200 °C, 400 °C, and 600 °C, respectively, while a 53% increase in time was observed for 800 °C. A significant delay in the flow development at 800 °C was observed due to significantly higher frictional forces.

14 SOLAR ENERGY↗
DOE OSTI2017

Effects of defects in composite wind turbine blades – Part 2: Progressive damage modeling of fiberglass-reinforced epoxy composites with manufacturing-induced waves

Abstract. Composite wind turbine blades are typically reliable; however, premature failures are often in regions of manufacturing defects. While the use of damage modeling has increased with improved computational capabilities, they are often performed for worst-case scenarios in which damage or defects are replaced with notches or holes. To better understand and predict these effects, an effects-of-defects study has been undertaken. As a portion of this study, various progressive damage modeling approaches were investigated to determine if proven modeling capabilities could be adapted to predict damage progression of composite laminates with typical manufacturing flaws commonly found in wind turbine blades. Models were constructed to match the coupons from, and compare the results to, the characterization and material testing study presented as a companion. Modeling methods were chosen from established methodologies and included continuum damage models (linear elastic with Hashin failure criteria, user-defined failure criteria, nonlinear shear criteria), a discrete damage model (cohesive elements), and a combined damage model (nonlinear shear with cohesive elements). A systematic, combined qualitative–quantitative approach was used to compare consistency, accuracy, and predictive capability for each model to responses found experimentally. Results indicated that the Hashin and combined models were best able to predict material response to be within 10 % of the strain at peak stress and within 10 % of the peak stress. In both cases, the correlation was not as accurate as the wave shapes were changed in the model; correlation was still within 20 % in many cases. The other modeling approaches did not correlate well within the comparative framework. Overall, the results indicate that this combined approach may provide insight into blade performance with known defects when used in conjunction with a probabilistic flaw framework.

Nelson, Jared W. (ORCID:0000000207736091)↗
NASA NTRS1978

Effect of longitudinal stringers on sound transmission into a thin cylindrical shell

In the context of the transmission of airborne noise into an aircraft fuselage, a mathematical model is presented for the transmission of airborne noise into a stiffened cylindrical shell. The stiffeners are longitudinal stringers and are modeled as discrete structural elements. The numerical cases examined were typical of a narrow-bodied jet transport fuselage. The stringers appeared to raise the cylinder transmission loss in the mass-controlled region, although they produced dips at the stringer resonances. The ring-frequency dip in transmission loss, which is characteristic of monocoque shells, was found to still be present. There appeared to be a small increase in transmission loss as the number of stringers was increased. The effect of high damping was also investigated.

Koval, L. R.↗
NASA NTRS1986

Parameter sensitivity in the dynamics of rotor-bearing systems

When designing a rotor system it is frequently desirable to have at hand a set of design sensitivity coefficients which quantitatively predict a change in specific system characteristics to changes in design parameters. This paper presents eigenvalue sensitivity coefficients for the damped natural frequencies of whirl of general linear rotor system modelled by finite element discretization. In addition, a simple and direct method for calculation of the damped critical speeds is presented, which utilizes the eigenvalue sensitivity with respect to the spin speed. It is shown that the combination of design parameter and spin speed whirl frequency sensitivity coefficients may be used to also evaluate the damped critical speed sensitivity coefficients.

Rajan, M.↗
NASA NTRS1993

Model reduction for the dynamics and control of large structural systems via neutral network processing direct numerical optimization

Three neural network processing approaches in a direct numerical optimization model reduction scheme are proposed and investigated. Large structural systems, such as large space structures, offer new challenges to both structural dynamicists and control engineers. One such challenge is that of dimensionality. Indeed these distributed parameter systems can be modeled either by infinite dimensional mathematical models (typically partial differential equations) or by high dimensional discrete models (typically finite element models) often exhibiting thousands of vibrational modes usually closely spaced and with little, if any, damping. Clearly, some form of model reduction is in order, especially for the control engineer who can actively control but a few of the modes using system identification based on a limited number of sensors. Inasmuch as the amount of 'control spillover' (in which the control inputs excite the neglected dynamics) and/or 'observation spillover' (where neglected dynamics affect system identification) is to a large extent determined by the choice of particular reduced model (RM), the way in which this model reduction is carried out is often critical.

Becus, Georges A.↗
NASA NTRSDate not supplied

OceanWATERS Lander Robotic Arm Operation

Ocean Worlds Autonomy Testbed for Exploration Research and Simulation (OceanWATERS) is an open-source simulator for developing onboard autonomy software for robotic exploration of ocean worlds, such as Europa, Enceladus, and Titan, built on the Robot Operating System (ROS) and Gazebo simulation environment. Inevitable ground communication delays increase demand for a high degree of autonomy during excavation, collection and transfer of samples to scientific instruments for in-situ analysis. This paper offers a detailed discussion of the robotic arm design and operation for such autonomous surface exploration, taking as reference the Europa Lander mission. The lander arm, which is designed primarily to acquire icy surface and subsurface samples within the arm’s workspace, is a 6-degree-of-freedom manipulator with two end effectors: a sample excavation tool and a trenching end-effector. The robotic arm’s modes and operations can be summarized as follows: stowed arm, intended as the lander arm default configuration characterized by zero-power consumption; un-stowed arm, target arm configuration after its first deployment; selection and deployment of the end-effector to use next; guarded move, to detect ground level at the desired trenching location; drill ice using the grinder; dig trench at a particular location using the scoop; deliver sample to the sample transfer dock; discard redundant samples. The motion planning tool used for the lander arm is MoveIt, a ROS package. MoveIt uses sampling-based planning and collision checking libraries to determine safe paths. The Rapidly Exploring Random Trees* (RRT*) has been chosen as default planning algorithm as it provides optimal plans with an exponential speed and is guaranteed to find a solution, if feasible solutions exist. Furthermore, this work quantifies and discusses the energy requirements for excavating and collecting samples. In OceanWATERS, force feedback from the terrain, which influences the arm dynamics, is modelled using a discrete element method (DEM) simulation. The DEM and Gazebo software run in parallel and communicate through a co-simulation plugin. This paper presents an analysis and comparison of three DEM open source software (YADE, ESyS-Particle, Project Chrono) for implementation in OceanWATERS and motivates the choice of YADE as most suitable candidate.

Damiana Catanoso↗
NASA NTRS2024

Experimental Validation of Charged Lunar Dust Dynamics Simulants

The goal of this work is to experimentally verify the electrostatic interaction physics that have been incorporated into the discrete element method (DEM) modeling software package during the FY20 CIF project titled “Charged Particle Dynamics in the Lunar Environment.” Current state-of-the-art (SOA) granular gas dynamics models used to explain rocket plume impingement physics have not taken into account the natural or induced electrostatic environment of the lunar surface, nor the effect of charged regolith grains being present in the plume plasma. This work improved upon the SOA by adding long-range and contact inter-particle electrostatic interactions to a granular mechanics DEM modeling software package. Tribocharging mechanics between spherical particles has been experimentally examined along with the natural electrostatic plasma environments produced via emission from an electron gun and ultraviolet lamp.

Jay Phillips↗
DOE OSTI2018

Flow Behavior of Particulate Pine Forest Residues and Corn Stover: A Comparison of Experiments and Simulations

Feeding and handling problems have been recognized as a major contributor to low production of renewable cellulosic ethanol in the US, which achieved only 7% of nameplate capacity in 2016. In this work we compare the flow behavior of pine forest residues and corn stover at various particle sizes and moisture contents in instrumented pilot-scale hopper tests and instrumented shear testers to continuum finite element method (FEM) and discrete element method (DEM) models. The shear testers include a Schulze ring shear tester, a custom large Peschl-style tester, and a modified Johanson direct shear tester. In addition to measurements of normal force and torque, the shear testers are also instrumented with arrays of pressure transducers. Modeling approaches include a Drucker-Prager model, a micropolar hypoelasticity model, and DEM simulations employing non-spherical, deformable particles. Results from physical experiments indicate that continuum elastoplastic models, such as Drucker-Prager models, do not account for the full complexity of biomass flow, and that more advanced rheological models are better able to predict nuances of the flow behavior, especially the compressive stress field.

09 BIOMASS FUELS↗
DOE OSTI2018

Effects of curvature-related DEM contact model on the macro- and micro-mechanical behaviours of granular soils

The presented paper conducted a comprehensive comparison between the Hertz-Mindlin model and the linear spring model in true triaxial shear simulations of granular soils using the discrete element method (DEM). The no-slip Hertz-Mindlin model for general elastic non-spherical particles with smooth surfaces was revisited, and implemented for superellipsoidal particles in an in-house DEM code. Three groups of specimens with a grain size distribution of Ottawa 20-30 sands, consisting of spheres, ellipsoids, and superellipsoids, respectively, were subjected to triaxial shear DEM simulations with the Hertz-Mindlin model and the linear spring model. The corresponding mechanical behaviours were examined in terms of a series of macro-and micro-parameters. It is found that the linear spring model is able to resemble the Hertz-Mindlin model on aspects of both microscopic and macroscopic mechanical behaviours of granular media with spherical and/or non-spherical particles. This finding suggests that the linear spring model can be used to investigate micro-mechanical behaviours of granular soils even with complex particle shapes.

42 ENGINEERING↗
DOE OSTI2018

New Mechanistic Models of Long Term Evolution of Microstructure and Mechanical Properties of Nickel Based Alloys

The report describes the development of a discrete element method (DEM) based modeling approach to quantitatively predict deformation and failure of typical nickel based superalloys. A series of experimental data, including microstructure and mechanical property characterization at 600°C, was collected for a relatively simple, model solid solution Ni-20Cr alloy (Nimonic 75) to determine inputs for the model and provide data for model validation. Nimonic 75 was considered ideal for this study because it is a certified tensile and creep reference material. A series of new DEM modeling approaches were developed to capture the complexity of metal deformation, including cubic elastic anisotropy and plastic deformation both with and without strain hardening. Our model approaches were implemented into a commercially available DEM code, PFC3D, that is commonly used by engineers. It is envisioned that once further developed, this new DEM modeling approach can be adapted to a wide range of engineering applications.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Study the thermal radiation effects in gas-solid flows with gray and non-gray P1 models implemented in MFiX

Thermal radiation is a dominant mode of heat transfer in combustion/gasification, packed/circulating bed reactors, and energy storage/conversion devices. In hydrocarbon fuel combustion, the high absorption and emission of product gases like CO 2 and H 2 O significantly affect the heat transfer characteristics. It is essential to understand the radiative energy propagation in such applications. In the present study, we developed a framework to implement the P-1 radiation model in MFiX for Eulerian-Eulerian (Two fluid-TFM) and Eulerian-Lagrangian (Discrete element-DEM and Particle in cell-PIC) models. The P-1 is a simplification of the spherical harmonics method. The framework has been verified on various flue gas conditions. The verified solver results are validated against the benchmark results available in the literature. The verified and validated gray and non-gray weighted sum of gray gas models is evaluated for dry and wet flue gas conditions. Further, they are also used to study gas-solid flows in a fluidized bed/riser.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2020

Absorber Column CFD model validation against PNNL’s device-scale absorber column on the LCFS unit

Absorber column has been widely used for CO 2 capture in the coal-fired power plants. High-fidelity CFD models play an important role in absorber column design and solvents optimization, which help enhance the CO 2 capture efficiency and reduce the operation cost. This report provides a comprehensive description of the development of CFD absorber models at two different levels, the design and implementation of PNNL’s device-scale absorber column experiment, and the methodology to combine the CFD results, experiment data, and Aspen model for a better understanding of the interface area in packed column. A composite model is firstly proposed in the Discrete Element Method (DEM) packing process, which can model the complex geometry of various packing elements. This generates a realistic packing pattern and accurate packing porosity ε and specific area a p compared to the actual values for absorber column used in LCFS. Two level of CFD absorber models were developed, namely the full-size column model (FCM) to simulate the entire packed column with a focus on the wall/entrance effects, and the representative column model (RCM) to simulate a section of column with a focus on the sensitivity study of interface area. A Design of Experiment (DoE) plan was developed to guide the collection of 100 run CFD data and 12 experiment runs. The 100 CFD runs were carried out in the RCM with Pro-Pak packing and cover a wide operation range and solvent properties. Impact of influential factors on the interface area in packed column were investigated in details. The CFD interface area was then combined with the experimental data and Aspen model to infer some information of the effective contact angle in the column. The accuracy and uncertainties in the interface area and contact angle can be quantified.

01 COAL, LIGNITE, AND PEAT↗
NASA NTRS2012

Implementation of a Smeared Crack Band Model in a Micromechanics Framework

The smeared crack band theory is implemented within the generalized method of cells and high-fidelity generalized method of cells micromechanics models to capture progressive failure within the constituents of a composite material while retaining objectivity with respect to the size of the discretization elements used in the model. An repeating unit cell containing 13 randomly arranged fibers is modeled and subjected to a combination of transverse tension/compression and transverse shear loading. The implementation is verified against experimental data (where available), and an equivalent finite element model utilizing the same implementation of the crack band theory. To evaluate the performance of the crack band theory within a repeating unit cell that is more amenable to a multiscale implementation, a single fiber is modeled with generalized method of cells and high-fidelity generalized method of cells using a relatively coarse subcell mesh which is subjected to the same loading scenarios as the multiple fiber repeating unit cell. The generalized method of cells and high-fidelity generalized method of cells models are validated against a very refined finite element model.

Pineda, Evan J.↗