A Performance Portable Sea Ice Model Using the Discrete Element Method
Explore the source record for details and available documents.
Peterson, Kara J.↗
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.
A Performance Portable Discrete Element Sea Ice Model for Earth System Modeling.
Abstract not provided.
Peterson, Kara J.↗
A discrete-continuum coupled finite element modelling approach for fibre reinforced concrete
Highlights: • An easy-to-implement but effective meso-scale modelling approach is developed for FRC. • It is capable of simulating and interpreting all main failure mechanisms in FRC. • It is critically validated by experiments with single or multiple random fibres under various loading conditions. • It can be used to optimise key design factors of FRC at meso/material scale by parametric studies. - Abstract: Fibre reinforced concrete (FRC) exhibits complicated failure modes such as fibre breakage, mortar cracking, crushing and spalling, fibre-mortar interfacial debonding, depending on many material properties, geometric dimensions, boundary and loading conditions. Most existing numerical models are unable to reproduce these failure modes that may occur simultaneously or sequentially in a specimen, mainly due to difficulties in generating finite element meshes with a large number of randomly-oriented fibres. Herein we develop a discrete-continuum coupled finite element modelling approach for FRC materials capable of effectively simulating all the major failure modes. The continuum damaged plasticity model is used to simulate damage and fracture behaviour of the mortar, while debonding of fibre-mortar interfaces is modelled by nonlinear cohesive interfacial elements. Unique techniques are devised to generate conforming meshes between fibres and the surrounding mortar so that the randomly-oriented fibres are easily modelled. The modelling approach is validated by simulating single fibre pullout tests with different inclination angles, notched and non-notched direct tensile tests and three-point bending beam tests with randomly-distributed multiple fibres.
36 MATERIALS SCIENCE↗
Computational Investigation of Melt Pool Process Dynamics and Pore Formation in Laser Powder Bed Fusion
In the laser powder bed fusion additive manufacturing process, the presence of porosity may result in cracks and significantly affects the part performance. A comprehensive understanding of the melt pool process dynamics and porosity evolution can help to improve build quality. In this study, a novel multi-physics computational fluid dynamics (CFD) model has been applied to investigate the fluid dynamics in melt pools and resultant pore defects. To accurately capture the melting and solidification process, major process physics, such as the surface tension, evaporation as well as laser multi-reflection, have been considered in the model. A discrete element method is utilized to model the generation of powder spreading upon build plate by additional numerical simulations. Multiple single track experiments have been performed to obtain the melt pool shape and cross-sectional dimension information. The predicted melt pool dimensions were found to have a reasonable agreement with experimental measurements, e.g., the errors are in the range of 1.3 to 10.6% for melt pool width, while they are between 1.4 and 15.9% for melt depth. Pores are captured by both CFD simulation and x-ray computed tomography measurement for the case with a laser power of 350 W and laser speed of 100 mm/s. The formation of keyholes maybe related to the melt pool front wall angle, and it is found that the front wall angle increases with the increase in laser line energy density. In addition, a larger laser power or smaller scanning speed can help to generate keyhole-induced pores; they also contribute to produce larger sized pores.
36 MATERIALS SCIENCE↗
Integrated Process Optimization for Biochemical Conversion
This research is motivated by the challenges faced during biomass processing in bioenergy plants. It has been observed that variations in biomass characteristics, such as moisture, ash, and carbohydrate contents cause variations in feeding of the system which led to underutilization of equipment and the reactor. The objective of this research is to ensure a continuous flow of biomass to the reactor in plants that use the biochemical conversion process to generate liquid fuels. The overall goal is to lower the cost of producing biofuels, which could lead to improving US’s energy independency and growing US’s rural economy. The research team developed analytical models, such as discrete element method (DEM) models and mathematical models. The DEM models are unit-level models that explicitly capture biomass characteristics and quantify the impacts of biomass characteristics on bulk material properties and the performance of specific equipment. The mathematical models are system-level models that capture the impacts of system infeed rate, equipment processing rate, storage location and capacity, and biomass characteristics on system throughput. The functional relations predicting the bulk material properties from DEM models are incorporated to the mathematical models. The models developed were validated and evaluated using data collected at Idaho National Laboratory’s biomass processing facility. Via these models, we identified process control strategies that ensure a continuous flow of biomass to the reactor, while meeting the requirements of biochemical conversion process. Our analysis indicates that sequencing of biomass bales based on moisture level, and carbohydrate contents could have a positive impact on reducing processing time and inventory level and increasing throughput rate. Short bale sequences that repeat frequently, seem to have the greatest impact on improving system’s performance. Based on our experiments, the total annual system operating costs reduced by 20-30%, and the maximum inventory level reduced by 3 to 4 times. The operating costs include the annual equipment amortization cost and processing cost. The implementation of the models developed requires the use of standardized bale format, Radio Frequency Identification technology, sensing and real time monitoring of material attributes, automated material handling equipment, and automated process control. The scope of the model proposed can be extended to include the whole supply chain. The supply chain models help identify how many bales of different biomass feedstock to purchase given biomass availability in the region, biomass price and quality, and the biomass processing capabilities of the biorefinery. Thus, the outcomes of supply chain models can be used to inform the design of long-term contracts among farmers and the biorefinery.
09 BIOMASS FUELS↗
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.↗
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.↗
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↗
Conservative Remap and Interpolation for the Discrete Element for Sea Ice Model.
Abstract not provided.
Peterson, Kara J.↗
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.↗
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↗
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.↗
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↗
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.↗
Modeling of rock scour using coupled 3-D discrete element and Lattice Boltzmann Methods
Copyright © 2018 ARMA, American Rock Mechanics Association. Scour of rock in unlined rock spillway channels is a critical issue facing many of the world’s dams. From a modeling point of view this poses a challenging and interesting problem that combines rock mechanics and hydraulics of turbulent flow. We analyze this interaction between the blocky rock mass and water by directly modeling the solid and fluid phases—the individual polyhedral blocks are modeled using the Discrete Element Method (DEM) while the water is modeled using the Lattice Boltzmann Method (LBM). The LBM mesh is entirely independent of the DEM discretization, making it possible to refine the LBM mesh such that transient and varied fluid pressures acting of the rock surface are directly modeled. This provides the capability to investigate the effect of water pressure inside the fractured rock mass, along potential sliding planes, and can be extended to rock falls and slides into standing bodies of water such as lakes and reservoirs. Herein we present preliminary results to demonstrate the capabilities of the methodology.
Gardner, MH↗
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↗
Quantification of Powder Spreading Process for Metal Additive Manufacturing
In powder bed additive manufacturing (AM), the powder spreading process can lead to inhomogeneity in layer height, porosity level, and particle size distribution. Hence, creating homogeneous layers is the first task for optimal printing. The conventional assessment methods for powder quality are limited to a measurement of powder morphology, flowability, particle size distribution, surface roughness, and packing density. However, those techniques may not be the most appropriate to quantify the powder bed characteristics since it is not able to capture local particle configuration and does not provide any vectoral information (e.g., arrangement of contacts between particles during spreading). In this project, we: 1) developed a powder spreading model using discrete element method (DEM) for multiple-layered process, 2) experimentally measure powder spreading dynamics including static angle of repose (AOR) and dynamic AOR in order to calibrate and subsequently validate the model, and 3) investigated the effect of particle spreading process parameters on spatial powder bed structure and quality using scalar (e.g., porosity, surface roughness, coordination number, and particle size distribution) and tensorial (e.g., directionality of inter-particle contacts) methods for a comprehensive quantification of powder bed.
36 MATERIALS SCIENCE↗
Numerical modelling of impact seismic signals on regolith [Slides in French]
Abstract not provided.
58 GEOSCIENCES↗