Search NASA⌕ Search

SEARCH · Search NASA

Results for “discretization methods”

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
DOE OSTI2017

Energy stable discontinuous Galerkin methods for Maxwell's equations in nonlinear optical media

The propagation of electromagnetic waves in general media is modeled by the time-dependent Maxwell's partial differential equations (PDEs), coupled with constitutive laws that describe the response of the media. In this work, we focus on nonlinear optical media whose response is modeled by a system of first order nonlinear ordinary differential equations (ODEs), which include a single resonance linear Lorentz dispersion, and the nonlinearity comes from the instantaneous electronic Kerr response and the residual Raman molecular vibrational response. To design efficient, accurate, and stable computational methods, we apply high order discontinuous Galerkin discretizations in space to the hybrid PDE-ODE Maxwell system with several choices of numerical fluxes, and the resulting semi-discrete methods are shown to be energy stable. Under some restrictions on the strength of the nonlinearity, error estimates are also established. When we turn to fully discrete methods, the challenge to achieve provable stability lies in the temporal discretizations of the nonlinear terms. To overcome this, novel strategies are proposed to treat the nonlinearity in our model within the framework of the second-order leap-frog and implicit trapezoidal time integrators. The performance of the overall algorithms are demonstrated through numerical simulations of kink and antikink waves, and third-harmonic generation in soliton propagation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRSDate not supplied

Implementation of Charged Particle Behavior in Discrete Element Method (DEM) Simulations

To understand the behavior of charged lunar regolith when perturbed by lunar landers, it is important to couple the grain dynamics with mechanical and electrical particle interactions. To accomplish this, improvements have been made to a discrete element method (DEM) software package to include both short- and long-range interactions between spherical particles. Short-range interactions rely on contact between the particles, such as electrical conduction and triboelectric charge transfer driven by work functions. Long-range interactions act at a distance between every pairing of particles, such as electrostatic forces and gravitational forces. Results from simulations between a few particles are compared with theory to verify these added behaviors prior to scaling up to more complex scenarios. The radii, initial charges, electrical conductivities, work functions, and separation of the particles are varied and the resultant charges as well as the time required to reach the final state are determined.

Electrostatics↗
DOE OSTI2015

Modeling the Interaction Between Hydraulic and Natural Fractures Using Dual-Lattice Discrete Element Method

The interaction between hydraulic fractures (HF) and natural fractures (NF) will lead to complex fracture networks due to the branching and merging of natural and hydraulic fractures in unconventional reservoirs. In this paper, a newly developed hydraulic fracturing simulator based on discrete element method is used to predict the generation of complex fracture network in the presence of pre-existing natural fractures. By coupling geomechanics and reservoir flow within a dual lattice system, this simulator can effectively capture the poro-elastic effects and fluid leakoff into the formation. When HFs are intercepting single or multiple NFs, complex mechanisms such as direct crossing, arresting, dilating and branching can be simulated. Based on the model, the effects of injected fluid rate and viscosity, the orientation and permeability of NFs and stress anisotropy on the HF-NF interaction process are investigated. Combined impacts from multiple parameters are also examined in the paper. The numerical results show that large values of stress anisotropy, intercepting angle, injection rate and viscosity will impede the opening of NFs.

02 PETROLEUM↗
NASA NTRS2015

Calculated Low-Speed Steady and Time-Dependent Aerodynamic Derivatives for Some Airfoils Using a Discrete Vortex Method

This paper contains a collection of some results of four individual studies presenting calculated numerical values for airfoil aerodynamic stability derivatives in unseparated inviscid incompressible flow due separately to angle-of-attack, pitch rate, flap deflection, and airfoil camber using a discrete vortex method. Both steady conditions and oscillatory motion were considered. Variables include the number of vortices representing the airfoil, the pitch axis / moment center chordwise location, flap chord to airfoil chord ratio, and circular or parabolic arc camber. Comparisons with some experimental and other theoretical information are included. The calculated aerodynamic numerical results obtained using a limited number of vortices provided in each study compared favorably with thin airfoil theory predictions. Of particular interest are those aerodynamic results calculated herein (such as induced drag) that are not readily available elsewhere.

Riley, Donald R.↗
DOE OSTI2016

High-order solution methods for grey discrete ordinates thermal radiative transfer

This work presents a solution methodology for solving the grey radiative transfer equations that is both spatially and temporally more accurate than the canonical radiative transfer solution technique of linear discontinuous finite element discretization in space with implicit Euler integration in time. We solve the grey radiative transfer equations by fully converging the nonlinear temperature dependence of the material specific heat, material opacities, and Planck function. The grey radiative transfer equations are discretized in space using arbitrary-order self-lumping discontinuous finite elements and integrated in time with arbitrary-order diagonally implicit Runge–Kutta time integration techniques. Iterative convergence of the radiation equation is accelerated using a modified interior penalty diffusion operator to precondition the full discrete ordinates transport operator.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

A discrete element method-based approach to predict the breakage of coal

Pulverization is an essential pre-combustion technique employed for solid fuels, such as coal, to reduce particle sizes. Smaller particles ensure rapid and complete combustion, leading to low carbon emissions. Traditionally, the resulting particle size distributions from pulverizers have been informed by empirical or semi-empirical approaches that rely on extensive data gathered over several decades during operations or experiments. However, the predictive capabilities for new coals and processes are limited. This work presents a Discrete Element Method based computational framework to predict particle size distribution resulting from the breakage of coal particles characterized by the coal’s physical properties. The effect of certain operating parameters on the breakage behavior of coal particles also is examined.

Brazilian Test.↗
DOE OSTI2024

Investigation of thermal damage in explosive bridgewire detonators via discrete element method simulations

Exploding bridgewire (EBW) detonators are used to rapidly and reliably initiate energetic reactions by exploding a bridgewire via Joule heating. While the mechanisms of EBW detonators have been studied extensively in nominal conditions, comparatively few studies have addressed thermally damaged detonator operability. We present a mesoscale simulation study of thermal damage in a representative EBW detonator, using discrete element method (DEM) simulations that explicitly account for individual particles in the pressed explosive powder. We use a simplified model of melting, where solid spherical particles undergo uniform shrinking, and fluid dynamics are ignored. The subsequent settling of particles results in the formation of a gap between the solid powder and the bridgewire, which we study under different conditions. In particular, particle cohesion has a significant effect on gap formation and settling behavior, where sufficiently high cohesion leads to coalescence of particles into a free-standing pellet. This behavior is qualitatively compared to experimental visualization data, and simulations are shown to capture several key changes in pellet shape. We derive a minimum and maximum limit on gap formation during melting using simple geometric arguments. In the absence of cohesion, results agree with the maximum gap size. With increasing cohesion, the gap size decreases, eventually saturating at the minimum limit. In conclusion, we present results for different combinations of interparticle cohesion and detonator orientations with respect to gravity, demonstrating the complex behavior of these systems and the potential for DEM simulations to capture a range of scenarios.

cohesive powders↗
NASA NTRS2020

Analysis of Sample Acquisition Dynamics Using Discrete Element Method

The analysis presented in this paper is conducted in the framework of the Ocean Worlds Autonomy Testbed for Exploration Research and Simulation (OceanWATERS) project, currently under development at NASA Ames Research Center. OceanWATERS aims at designing a simulation environment which allows for testing autonomy of scientific lander missions to the icy moons of our solar system. Mainly focused on reproducing the end effector interaction with the inherent terrain, this paper introduces a novel discrete element method (DEM)-based approach to determine forces and torques acting on the lander’s scoop during the sample acquisition process. An accurate force feedback from the terrain on the scoop is required by fault-detection and autonomous decision-making algorithms to identify when the requested torque on the robotic arm’s joints exceeds the maximum available torque. Knowledge of the terrain force feedback significantly helps evaluating the arm’s links structural properties and properly selecting actuators for the joints. Models available in literature constitute a partial representation of the dynamics of the interaction. As an example, Balovnev derived an analytical expression of the vertical and horizontal force acting on a bucket while collecting a sample as a function of its geometry and velocity, soil parameters and reached depth. Although the model represents an adequate approximation of the two force components, it ignores the direction orthogonal to the scoop motion and neglects the torque. This work relies on DEM analysis to compensate for analytical models’ deficiencies and inaccuracies, i. e. provide force and torque 3D vectors, defined in the moving reference (body) frame attached to the scoop, at each instant of the sample collection process. Results from the first presented analysis relate to the specific OceanWATERS sampling strategy, which consists of collecting the sample through five consecutive passes with increasing depth, each pass following the same circularlinear- circular trajectory. Data is collected given a specific scoop design interacting with two types of bulk materials, which may characterize the surface of icy planetary bodies: snow and ice. Although specifically concerned with the OceanWATERS design, this first analysis provides the expected force trends for similar sampling strategies and allows to deduce phenomenological information about the general scooping process. In order to further instruct the community on the use of DEM tools as a solution to the sampling collection problem, two more analyses have been carried out, mainly focused on reducing the DEM computation time, which increases with a decrease in particle size. After running a set of identical simulations, where the only changing parameter is the size of the spherical particle, it is observed that the resulting force trajectories, starting from a given particle size, converge to the true trend. It is deducible that a further decrease in size yields negligible improvements in the accuracy, while it sensibly increases computation time. A final analysis aims at discussing limitations of approximating bulk material particles having a complex shape, e. g. ice fragments, with spheres, by comparing force trends resulting in the two cases for the same simulation scenario.

Catanoso, Damiana↗
NASA NTRS1971

Comparison of a discrete steepest ascent method with the continuous steepest ascent method for optimal programing

A discrete steepest ascent method which allows controls which are not piecewise constant (for example, it allows all continuous piecewise linear controls) was derived for the solution of optimal programming problems. This method is based on the continuous steepest ascent method of Bryson and Denham and new concepts introduced by Kelley and Denham in their development of compatible adjoints for taking into account the effects of numerical integration. The method is a generalization of the algorithm suggested by Canon, Cullum, and Polak with the details of the gradient computation given. The discrete method was compared with the continuous method for an aerodynamics problem for which an analytic solution is given by Pontryagin's maximum principle, and numerical results are presented. The discrete method converges more rapidly than the continuous method at first, but then for some undetermined reason, loses its exponential convergence rate. A comparsion was also made for the algorithm of Canon, Cullum, and Polak using piecewise constant controls. This algorithm is very competitive with the continuous algorithm.

Childs, A. G.↗
DOE OSTI2018

Assessment of different coarse graining strategies to simulate polydisperse gas-solids flow

Continuum methods require the additional development of solids stress closures for polydisperse powders based on complex kinetic theories that are non-trivial to develop, code, and numerically converge for the wide range of fluidization regimes from very dilute to dense/frictional flow limit. On the other hand, it is straightforward to model the flow of polydisperse granular materials by treating particles as discrete rigid bodies that are tracked following simple physical laws of motion. The coarsening of these discrete methods by lumping several particles in a parcel alleviates the significant computational cost associated with these discrete methods while introducing some inaccuracies in the numerical results. In this research, we explore two different coarse graining methods that can be applied to polydisperse powders, namely the same statistic weight method (SSW) and the same size parcel method (SSP), and assess their accuracy by comparison with the finest simulation results obtained with a discrete element method (DEM). For Geldart group B powders fluidized at a relative low superficial velocity, the numerical results indicate that the SSW is more accurate than the SSP method. For type A powders fluidized at relatively high velocity, these two methods predict similar results. Interestingly, up to four times increase in the speed of simulation of the SSP method was obtained because the original polydisperse powder is scaled to a mono-disperse system in terms of particle-particle collision. Lastly, these results suggest that the SSP method is more favorable for the simulation of fluidized beds due to its accuracy and efficiency while the SSW method may be used for granular flow and dense fluidized bed systems where capturing the size segregation of particles due to collision is important.

42 ENGINEERING↗
NASA NTRSDate not supplied

Implementation of Charged Particle Behavior in Discrete Element Method (DEM) Simulations

Lunar landers will agitate the surface of the Moon with an exhaust plume during descent which will, due to the particulate nature of the lunar regolith, loosen and eject grains from the surface. This ejection is not only coupled with the charged plume gas, but also results in significant particle-particle interactions. Settling of these grains after plume effects have subsided takes much longer than expected in a ballistic trajectory. The prevailing hypothesis attributes this behavior to the accumulated charge on the particles. We are thus developing a discrete element method (DEM) approach to explore these charged particle interactions on the lunar surface. The Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) Improved for General Granular and Granular Heat Transfer Simulations (LIGGGHTS) software package provides a DEM modeling framework for granular interactions. It includes many complexities such as non-spherical particle shapes, cohesion and frictional forces, and heat transfer, but has no provision for inter-particle electrostatic forces and charge transfer that are important to examine in the lunar environment. In this work, a standard Coulomb potential and a Yukawa potential are integrated into the LIGGGHTS framework to provide a basis for particle-particle electrostatic interactions, as well as a gravitational potential to enable inter-grain gravitational interactions. A preliminary approach to charge transfer between grains incorporating properties such as work function and electrical conductivity to the library of available material characteristics will be presented. Several scenarios have been simulated that include charged particle interactions within a diffuse granular gas, settling of charged grains into a regolith bed, sliding of granular material along an incline, and vibration of settled grains to produce a behavior similar to fluidization. There are numerous challenges to incorporate realistic interactions between complex lunar particles. Currently, grains are modeled to behave as if the entirety of the charge acts at the center of mass, such as conductors with spherical symmetry and insulators with homogeneously distributed charge. We are developing improvements that will include the use of non-spherical particle geometries, as well as reasonable approximations of insulating/dielectric materials that have non-uniform charge distributions. The cases simulated thus far will be examined in a relevant environment within a vacuum chamber to validate the simulations. These simulations will be bounded by experiments utilizing high-speed camera observations of the motion for validation. The grains in the experiment will exchange charge during their motion and this can be quantified by collection within a charge measurement device such as a Faraday cup. Such a device may be modeled within the software by defining an integration region and computing the contained charge as a function of simulation time, allowing for side-by-side comparison of simulated and measured bulk charging results. Any differences will be reconciled by updating the mathematical mechanisms described within the simulation suite. Successfully combining results from experiments within a relevant environment into the LIGGGHTS framework will improve modeling of the charged grain dynamics experienced on the Moon to provide insights into dust behavior for future lunar exploration missions.

Electrostatics↗
DOE OSTI2018

Effects of Particle Diameter and Coke Layer Thickness on Solid Flow and Stress Distribution in BF by 3D Discrete Element Method

In order to reduce ironmaking-related CO{sub 2} emissions, hydrogen-enriched blast furnace (BF) operation is currently under development. In hydrogen-enriched BF operation, coke layer thickness can be decreased to reduce CO{sub 2} emissions. However, BFs operating with thin coke layers may experience instability or discontinuous phenomena such as particle slip and gas channeling problems, so it is important to optimize the particle diameter and coke layer thickness for optimal BF operation. In this study, the effects of particle diameter and coke layer thickness on the solid flow and stress distribution in a BF were analyzed using a three-dimensional discrete element method. Furthermore, the effects of particle diameter and coke layer thickness on the burden layer stabilities, particle velocities, and particle stress distributions have been investigated. The results show that decreasing the coke layer thickness caused instability owing to the mixing of the coke and ore layers in the BF-cohesive zone and slight increases in both the average particle velocities and the average normal particle stress magnitudes. In addition, the average particle velocities and average normal particle stresses were higher for the smaller particles than for the larger ones during the simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Fish Passage though Hydropower Turbines: Simulating Blade Strike using the Discrete Element Method

mong the hazardous hydraulic conditions affecting anadromous and resident fish during their passage though turbine flows, two are believed to cause considerable injury and mortality: collision on moving blades and decompression. Several methods are currently available to evaluate these stressors in installed turbines, i.e. using live fish or autonomous sensor devices, and in reduced-scale physical models, i.e. registering collisions from plastic beads. However, a priori estimates with computational modeling approaches applied early in the process of turbine design can facilitate the development of fish-friendly turbines. In the present study, we evaluated the frequency of blade strike and nadir pressure environment by modeling potential fish trajectories with the Discrete Element Method (DEM) applied to fish-like composite particles. In the DEM approach, particles are subjected to realistic hydraulic conditions simulated with computational fluid dynamics (CFD), and particle-structure interactions—representing fish collisions with turbine blades—are explicitly recorded and accounted for in the calculation of particle trajectories. We conducted transient CFD simulations by setting the runner in motion and allowing for better turbulence resolution, a modeling improvement over the conventional practice of simulating the system in steady state which was also done here. While both schemes yielded comparable bulk hydraulic performance, transient conditions exhibited a visual improvement in describing flow variability. We released streamtraces (steady flow solution) and DEM particles (transient solution) at the same location from where sensor fish (SF) have been released in field studies of the modeled turbine unit. The streamtrace-based results showed a better agreement with SF data than the DEM-based nadir pressures did because the former accounted for the turbulent dispersion at the intake but the latter did not. However, the DEM-based strike frequency is more representative of blade-strike probability than the steady solution is, mainly because DEM particles accounted for the full fish length, thus resolving (instead of modeling) the collision event.

CFD↗
DOE OSTI2016

A non-locking composite tetrahedron element for the combined finite discrete element method

Here, a number of validation cases considering different loading and boundary conditions and different materials for the proposed element are presented. A practical application of the use of the composite tetrahedral finite element is presented by quantitative comparisons of numerical model results against simple theoretical estimates and results from acrylic fracturing experiments. All of these examples clearly show the capability of the composite element in eliminating volumetric locking.

42 ENGINEERING↗
DOE OSTI2019

Yearly report for the period Jan. 2018 – Jan. 2019 IC Project: w17_unedetection - “Simulation of underground nuclear explosions using the combined finite discrete element method” and “Fracture Formation and Permeability Evolution at in situ Pressure, Temperature and Stress Conditions”

Understanding the underlying physical processes associated with an underground nuclear explosion (UNE) can help with detection and identification of a suspected event. The Underground Nuclear Explosion Signatures Experiment (UNESE) is a collaborative effort between national laboratories aimed towards improving our understanding of physical and radiochemical signatures associated with a UNE. One primary objective of this effort is to perform hydrodynamic modeling of a UNE, with the purpose of estimating the extent of damage caused by such an event. Estimating the damage associated with an event is important for: a) detection of possible physical signatures at the surface of the earth such as fractures, spallation, and cratering and (b) damage to subsurface structures such as tunnels. Additionally, the damage produced by the hydrodynamic models play can be used to inform subsurface flow and transport models. The Combined Finite-Discrete Element Method (FDEM), implemented as an in-house code known as HOSS (Hybrid Optimization Software Suite), is used to model the damage associated with a UNE in various geological settings.

42 ENGINEERING↗
DOE OSTI2026

Optimal Design and Techno-Economic Analysis of 3D-Printed, Intensified Packings for Absorbers and Strippers in Solvent-Based CO 2 Capture

A potential technology for the CO 2 absorption process is utilizing intensified structured packing with embedded cooling/heating channels for continuous heat exchange, which can overcome limitations of discrete methods, such as discrete intercooling and centralized reboilers, to aid in reducing energy consumption and decreasing costs. This work investigates the modeling of intensified packing (IP) for the stripper tower, extending on previous work for the absorber, which distributes heat internally within the column, improving the thermodynamics for the solvent regeneration process. The model includes submodels for steam turbine extraction to produce steam at various qualities as well as a surrogate model for calculating steam enthalpy. A cost model for a plant-scale absorption capture process was developed, allowing for the design of the plant to be optimized, subject to minimizing capture cost using two different power plant flue gas sources. In this optimization, the placement of IP in both towers is optimized to balance the trade-off between enhanced heat transfer and reduced mass transfer volume. For natural gas combined cycle flue gas, the standard process configuration had a minimum cost of $\$$65.40/tonne CO 2 , and considering IP, the minimum capture cost is reduced to $\$$62.73/tonne, with utilization in the stripper column, which reduces yearly costs by up to $\$$2.67 MM/yr. Cooling the absorber through IP, or intercoolers, was only found to be beneficial at higher capture rates, with IP in both towers having a cost of capture of $\$$68.08/tonne at 99.9% capture, a reduction of $\$$12.64/tonne when using only intercoolers at the same capture rate. When capturing from pulverized-coal power plants, the minimum cost of capture when using IP in both towers is $\$$44.18/tonne (at 97% capture), while the standard configuration with and without intercoolers was $\$$45.69 and $\$$47.22 per tonne, respectively. This results in a reduction in yearly costs of $\$$16.98 MM/yr from the base-case configuration. At this higher CO 2 concentration, cooling in the absorber from the IP becomes extremely beneficial, reducing energy consumption by up to 6%.

20 FOSSIL-FUELED POWER PLANTS↗