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 703 records · Page 39

Granular Flow in Novel Octet Shape–Based Lattice Frame Material

Moving packed-bed heat exchangers in concentrated solar power (CSP) plants involves heat transfer between heated falling particles and supercritical carbon dioxide. The overall effective thermal conductivity of the moving packed bed and particle-side channel contact resistances are still the bottlenecks in achieving the desirable thermal transport levels. To this end, a novel moving packed bed heat exchanger consisting of an Octet lattice packed between the walls of the particle-side channel is proposed in this study. Granular flow analysis in Octet lattice moving packed bed heat exchanger (OLHX) was conducted through experiments and discrete element method (DEM)-based numerical simulations. The experimental images clearly demonstrated stagnation regions upstream of lattice fibers, void regions downstream of the fiber junctions, and wavy-type unobstructed flow on the lateral sides of the fibers. DEM simulations were successful in capturing all these critical flow phenomena. Larger flow velocities were observed on the lateral sides of the fibers in the simulations. Also, when the particles in the silo were emptied, the final images showed an accumulation of particles on the inter-fiber as well as fiber–channel wall junctions. Moreover, the fiber connections resulted in some regions devoid of particle contact on the channel endwall, which means that these regions would suffer from poor thermal exchange. Finally, the overall mass flowrate increased with increasing porosity for a fixed particle diameter.

14 SOLAR ENERGY↗

Two-Fluid and Discrete Element Modeling of a Parallel Plate Fluidized Bed Heat Exchanger for Concentrating Solar Power

A novel high-temperature particle solar receiver is developed using a light trapping planar cavity configuration. As particles fall through the cavity, the concentrated solar radiation warms the boundaries of the receiver and in turn heats the particles. Particles flow through the system, forming a fluidized bed at the lower section, leaving the system from the bottom at a constant flowrate. Air is introduced to the system as the fluidizing medium to improve particle heat transfer and mixing. A laboratory scale cavity receiver is built by collaborators at the Colorado School of Mines and their data are used for model validation. In this experimental setup, near IR quartz lamp is used to provide flux to the vertical wall of the heat exchanger. The system is modeled using the discrete element method and a continuum two-fluid method. The computational model matches the experimental system size and the particle size distribution is assumed monodisperse. A new continuum conduction model that accounts for the effects of solid concentration is implemented, and the heat flux boundary condition matches the experimental setup. Radiative heat transfer is estimated using a widely used correlation during the post-processing step to determine an overall heat transfer coefficient. The model is validated against testing data and achieves less than 30% discrepancy and a heat transfer coefficient greater than 1000 W/m2 K.

concentrating solar power↗

Comparison of Radiation Models for a Turbulent Piloted Methane/Air Jet Flame: A Frozen-Field Study

Numerical modeling of radiative transfer in nongray reacting media is a challenging problem in computational science and engineering. The choice of radiation models is important for accurate and efficient high-fidelity combustion simulations. Different applications usually involve different degrees of complexity, so there is yet no consensus in the community. In this paper, the performance of different radiative transfer equation (RTE) solvers and spectral models for a turbulent piloted methane/air jet flame are studied. The flame is scaled from the Sandia Flame D with a Reynolds number of 22,400. Three classes of RTE solvers, namely the discrete ordinates method, spherical harmonics method, and Monte Carlo method, are examined. The spectral models include the Planck-mean model, the full-spectrum k-distribution (FSK) method, and the line-by-line (LBL) calculation. The performances of different radiation models in terms of accuracy and computational cost are benchmarked. The results have shown that both RTE solvers and spectral models are critical in the prediction of radiative heat source terms for this jet flame. The trade-offs between the accuracy, the computational cost, and the implementation difficulty are discussed in detail. The results can be used as a reference for radiation model selection in combustor simulations.

David, Chloe↗

Comparison of Spectroscopic Data with Cluster Calculations of Plutonium, Plutonium Dioxide and Uranium Dioxide

Using spectroscopic data produced in the experimental investigations of bulk systems, including X-Ray Absorption Spectroscopy (XAS), Photoelectron Spectroscopy (PES) and Bremstrahlung Isochromat Spectroscopy (BIS), the theoretical results within for UO{sub 2}{sup 6}, PuO{sub 2}{sup 6} and Pu{sup 7} clusters have been evaluated. The calculations of the electronic structure of the clusters have been performed within the framework of the Relativistic Discrete-Variational Method (RDV). The comparisons between the LLNL experimental data and the Russian calculations are quite favorable. The cluster calculations may represent a new and useful avenue to address unresolved questions within the field of actinide electron structure, particularly that of Pu. Observation of the changes in the Pu electronic structure as a function of size suggests interesting implications for bulk Pu electronic structure.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Error Estimates for the Optimal Control of a Parabolic Fractional PDE

In this work, we consider the integral definition of the fractional Laplacian and analyze a linear-quadratic optimal control problem for the so-called fractional heat equation; control constraints are also considered. We derive existence and uniqueness results, first order optimality conditions, and regularity estimates for the optimal variables. To discretize the state equation we propose a fully discrete scheme that relies on an implicit finite difference discretization in time combined with a piecewise linear finite element discretization in space. We derive stability results and a novel $L^2(0,T;L^2(\Omega))$ a priori error estimate. On the basis of the aforementioned solution technique, we propose a fully discrete scheme for our optimal control problem that discretizes the control variable with piecewise constant functions, and we derive a priori error estimates for it. We illustrate the theory with one- and two-dimensional numerical experiments.

97 MATHEMATICS AND COMPUTING↗

Granular Flow

LIGGGHTS-INL is a capability-extended version of LIGGGHTS. LIGGGHTS is an Open Source Discrete Element Method Particle Simulation Software (https://www.cfdem.com/media/DEM/docu/Manual.html). LIGGGHTS stands for LAMMPS improved for general granular and granular heat transfer simulations. LAMMPS is a classical molecular dynamics simulator. It is widely used in the field of Molecular Dynamics (https://lammps.sandia.gov/doc/Manual.html). Thanks to physical and algorithmic analogies, LAMMPS is a very good platform for DEM simulations. LAMMPS offers a GRANULAR package to perform these kinds of simulations. LIGGGHTS aims to improve those capability with the goal to apply it to industrial applications.

Xia, Yidong↗

Quandary

Quandary numerically simulates and optimizes the time-evolution of open quantum systems. The underlying dynamics are modelled by Lindblad's master equation, a linear ordinary differential equation (ODE) describing quantum systems interacting with the environment. Quandary solves this ODE numerically by applying a time-stepping integration scheme, and utilizes a gradient-based optimization approach to determine optimal control pulses that drive the quantum system to a desired target state. Two optimization objectives are considered: (a) Unitary gate optimization that finds controls to realize a unitary gate transformation, and (b) optimal reset that aims to drive the quantum system to the ground states. Gradient-based optimization schemes utilizing Petsc's Tao optimization package are applied to generate control pulses that minimize the respective measure. To evaluate the gradient of the objective function, the discrete adjoint method is used while leveraging techniques from Algorithmic Differentiation to produce exact and consistent gradients. To mitigate excessive execution run times, the software can be build together with the XBraid software library which provides a parallelization strategy to distribute the time-evolution of the underlying dynamics onto multiple processor.

Petersson, NilsA.↗

BDEM (Discrete-element-simulator for high-solids granular flows) [SWR-22-72]

BDEM is a discrete element method based simulation tool developed specifically for modeling high-solids granular flows that include polydispersity, heat-transfer, moving boundaries and chemistry. Our solver provides facilities for simulating spherical/non-spherical particles with modified contact and friction models in complex dynamic geometries defined using level-sets or triangulated files. The solver is developed on top of NREL's open-source performance portable library, AMReX, providing parallel execution capabilities on current and upcoming high-performance-computing (HPC) architectures. Simulations at the scale of several millions to billion particles have been performed using this software on large scale computing resources. This software can be applied to non-reacting solids dominant flows in silos, hoppers and screw conveyors as well as in high temperature reacting systems such as screw kilns and auger reactors.

Sitaraman, Hariswaran↗

BARC Element Classifier (barc)

SAND2024-02120O This program enumerates and classifies all possible functional elements in several different categories within the Ballistic Asynchronous Reversible Computing (barc) model of computation. barc software is a research tool that is being used to help document the possible digital behaviors of primitive functional elements in the ABRC a.k.a. BARC model of computation. The method of operation of this software leverages elementary concepts and methods from discrete mathematics: combinatorics (permutations), partitions and equivalence classes, symmetry transformations, and the theory of finite groups. Implementation is in Python programming language. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Frank, Michael↗

Towards performance-portability of the Albany Finite Element analysis code using the Kokkos library of Trilinos

Performance portability on heterogeneous high-performance computing (HPC) systems is a major challenge faced today by code developers: parallel code needs to execute correctly as well as with high performance on machines with different architectures, operating systems, and software libraries. The Finite Element Method (FEM) is a popular and flexible method for discretizing partial differential equations arising in a wide variety of scientific, engineering, and industry applications that require HPC. This paper presents some preliminary results pertaining to our development of a performance portable implementation of the FEM-based Albany code. Performance portability is achieved using the Kokkos library of Trilinos. We present performance results for two different physics simulations modules in Albany: the Aeras global atmosphere dynamical code and the FELIX land-ice solver. As a result, numerical experiments show that our single code implementation gives reasonable performance across two multi-core/many-core architectures: NVIDIA GPUs and multi-core CPUs.

97 MATHEMATICS AND COMPUTING↗

Interaction Between Proppant Compaction and Single-/Multiphase Flows in a Hydraulic Fracture

Understanding proppant-transport and -deposition patterns in a hydraulic fracture is vital for effective and economical production of petroleum hydrocarbons. In this research, a numerical-modeling approach, combining the discrete-element method (DEM) with single-/multiphase lattice Boltzmann (LB) simulation, was adopted to advance understanding of the interaction between reservoir depletion, proppant-particle compaction, and single-/multiphase flows in a hydraulic fracture. DEM was used to simulate effective-stress increase and the resultant proppant-particle movement and rearrangement during the process of reservoir depletion caused by hydrocarbon production. Simulated pore structure of the proppant pack was extracted and used as boundary conditions in the LB simulation to calculate the time-dependent permeability of the proppant pack. We first validated the DEM-LB numerical work flow, and the simulated proppant-pack permeabilities as functions of effective stress were in good agreement with laboratory measurements. Furthermore, three proppant packs with the same average particle diameter but different diameter distributions were generated to study the role of proppant-size heterogeneity (variation in particle diameter). Specifically, we used the coefficient of variation (COV) of diameter, defined as the ratio of standard deviation of diameter to mean diameter, to characterize the heterogeneity of particle size. We obtained proppant-pack porosity, permeability, and fracture-width reduction (closure distance) as functions of effective stress. Under the same effective stress, a proppant pack with a higher diameter COV had lower porosity and permeability and larger fracture-width reduction. This was because the high diameter COV gave rise to a wider diameter distribution of proppant particles; smaller particles were compressed into the pore space between larger particles with increasing stress, leading to larger closure distance and lower porosity and permeability. With multiphase LB simulation, relative permeability curves were obtained, which are critical for larger-scale reservoir simulations under various effective stresses. The relative permeability of the oil phase was more sensitive to changes in geometry and stress, compared with the wetting phase. This was because the oil phase occupied larger pores; compaction of the proppant pack affected the structure of the pores, because the pores were farther from the grain contacts and thus more susceptible to collapse. It is also interesting to note that when effective stress increased continuously, oil relative permeability increased first and then decreased. This nonmonotonic behavior was the result of the nonmonotonic development of pore structure and oil connectivity under increasing stress.

Fan, Ming↗

Studies of elasticity, sound propagation and attenuation of acoustic modes in granular media: final report

This is the final report describing the results of DOE Grant # DE-FG02-03ER15458 with original termination date of April 31, 2013, which has been extended to April 31, 2014. The goal of this project is to develop a theoretical and experimental understanding of sound propagation, elasticity and dissipation in granular materials. The topic is relevant for the efficient production of hydrocarbon and for identifying and characterizing the underground formation for storage of either CO 2 or nuclear waste material. Furthermore, understanding the basic properties of acoustic propagation in granular media is of importance not only to the energy industry, but also to the pharmaceutical, chemical and agricultural industries. We employ a set of experimental, theoretical and computational tools to develop a study of acoustics and dissipation in granular media. These include the concept effective mass of granular media, normal modes analysis, statistical mechanics frameworks and numerical simulations based on Discrete Element Methods. Effective mass measurements allow us to study the mechanisms of the elastic response and attenuation of acoustic modes in granular media. We perform experiments and simulations under varying conditions, including humidity and vacuum, and different interparticle force-laws to develop a fundamental understanding of the mechanisms of damping and acoustic propagation in granular media. A theoretical statistical approach studies the necessary phase space of configurations in pressure, volume fraction to classify granular materials.

58 GEOSCIENCES↗

Implementation and Re nement of a Comprehensive Model for Dense Granular Flows

Dense granular ows are ubiquitous in both natural and industrial processes. They manifest three di erent ow regimes, each exhibiting its own dependence on solids volume fraction, shear rate, and particle-level properties. This research project sought to develop continuum rheological models for dense granular ows that bridges multiple regimes of ow, implement them in open-source platforms for gas-particle ows and perform test simulations. The rst phase of the research covered in this project involved implementation of a steady- shear rheological model that bridges quasi-static, intermediate and inertial regimes of ow into MFIX (Multiphase Flow with Interphase eXchanges - a general purpose computer code developed at the National Energy Technology Laboratory). MFIX simulations of dense granular ows in hourglass-shaped hopper were then performed as test examples. The second phase focused on formulation of a modi ed kinetic theory for frictional particles that can be used over a wider range of particle volume fractions and also apply for dynamic, multi- dimensional ow conditions. To guide this work, simulations of simple shear ows of identical mono-disperse spheres were also performed using the discrete element method. The third phase of this project sought to develop and implement a more rigorous treatment of boundary e ects. Towards this end, simulations of simple shear ows of identical mono-disperse spheres con ned between parallel plates were performed and analyzed to formulate compact wall boundary conditions that can be used for dense frictional ows at at frictional boundaries. The fourth phase explored the role of modest levels of cohesive interactions between particles on the dense phase rheology. The nal phase of this project focused on implementation and testing of the modi ed kinetic theory in MFIX and running bin-discharge simulations as test examples.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Numerical simulation of a full-loop circulating fluidized bed under different operating conditions

Both experimental and computational studies of the fluidization of high-density polyethylene (HDPE) particles in a small-scale full-loop circulating fluidized bed are conducted. Experimental measurements of pressure drop are taken at different locations along the bed. The solids circulation rate is measured with an advanced Particle Image Velocimetry (PIV) technique. The bed height of the quasi-static region in the standpipe is also measured. Comparative numerical simulations are performed with a Computational Fluid Dynamics solver utilizing a Discrete Element Method (CFD-DEM). This paper reports a detailed and direct comparison between CFD-DEM results and experimental data for realistic gas-solid fluidization in a full-loop circulating fluidized bed system. The comparison reveals good agreement with respect to system component pressure drop and inventory height in the standpipe. In addition, the effect of different drag laws applied within the CFD simulation is examined and compared with experimental results.

01 COAL, LIGNITE, AND PEAT↗

FY18 Progress on Improvements and V&V of MC 2 -3 and PROTEUS

The PROTEUS code is a high-fidelity capable deterministic neutron transport code based on unstructured finite element meshes, which solves the steady-state and transient neutron transport problem using the 2nd-order discrete ordinate method (SN2ND), the method of characteristics (MOC), and the advanced nodal transport method (NODAL). The PROTEUS system in a broader sense consists of cross section generation tools, mesh generation tools, and output post processing tools as well as neutron transport solvers.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Thermal Test Suite for Code Maintenance and Validation of Diablo

This document describes a series of test problems that ensure that the Thermal capabilities of Diablo are working correctly. Diablo utilizes the Finite Element Method to discretize the problem in space and a trapezoidal rule finite difference scheme to discretize the problem in time. This suite will test the codes capabilities to handle various boundary conditions, source terms, and initial conditions. It will also test the coupling between thermal and solid effects, and the use of various linear and nonlinear solvers.

42 ENGINEERING↗

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↗

IC reports (2021)

We requested HPC support to continue research on the seismic waves generated by impacts. We used the following codes: (1) the Hybrid Optimization Software Suite (HOSS), developed at LANL. HOSS is based on a combined Finite and Discrete Element Method (FDEM). New material models are developed for the sedimentary rocks. HOSS has been recently benchmarked to iSale and FLAG codes. SPECFEM3D is an open-source code developed since the last 90s. It won the Gordon Bell award for best performance in 2003, was finalist again in 2008 for a run at 0.16 petaflops on 149,784 cores on the ‘Jaguar’ Cray system at Oak Ridge National Laboratory. It also won the BULL Joseph Fourier supercomputing award in 2010.; SW4 is a 4th-order finite difference code developed at LLNL which is currently actively developed to handle complex 3D models and to be ported on future exascale platforms. We assessed our need to a total of 3.9M CPU-hrs for year 1 and and 3.1 M for year 2.

79 ASTRONOMY AND ASTROPHYSICS↗