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 1,369 records · Page 76
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Coarse-grained discrete particle simulations of particle segregation in rotating fluidized beds in vortex chambers [Discrete particle simulations of particle segregation in rotating fluidized beds in vortex chambers]

Vortex chambers allow the generation of rotating fluidized beds, offering high-G intensified gas-solid contact, gas-solids separation and solids-solids segregation. Focusing on binary particle mixtures and fixing the density and diameter of the heavy/large particles, transient batch CFD-coarse-grained DPM simulations were carried out with varying densities or sizes of the light/small particles to evaluate to what extent combining these three functionalities is possible within a vortex chamber of given design. Both the rate and quality of segregation were analyzed. Within a relatively wide density and size range, fast and efficient segregation takes place, with an inner and slower rotating bed of the lighter/small particles forming within the outer and faster rotating bed of the heavier/large particles. Simulations show that the contamination of the outer bed with lighter particles occurs more easily than contamination of the inner bed with heavier particles and increases with decreasing difference in size or density of the particles. Bubbling in the inner bed is observed with an inner bed of very low density or small particles. Porosity plots show that vortex chambers with a sufficient number of gas inlet slots have to be used to guarantee a uniform gas distribution and particle bed. Lastly, the flexibility of particle segregation in vortex chambers with respect to the gas flow rate is demonstrated.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2018

Classification of particle height in a hopper bin from limited discharge data using convolutional neural network models

Hopper bin discharge processes are ubiquitous in many industrial applications. The preponderance of previous studies of this system has been focused on the forward discharge process given an initial particle packing configuration. Motivated by the needs in reaction design and optimization, we develop here an inverse reconstruction procedure that enables one to obtain the particle height information in the initial packing configuration in the hopper bin from limited discharge dynamics data and particle characteristics. Our procedure is based on convolutional neural network (CNN) models, which take experimentally measurable particle residence time, diameter and density as input, and provide a classification of particle height in the initial packing as output. Using MFIX-DEM simulations with enhanced physical modeling capabilities, we generate extensive discharge data for hoppers containing two distinct solid phase particles with four distinct classes of initial packing geometries. The CNN reconstruction models are subsequently trained and tested using the discharge data. Here, we find that the reconstruction (i.e., particle height classification) accuracy strongly depends on the initial packing geometry. For hopper-specific CNN models exclusively trained using discharge data from a single hopper, the classification can attain an accuracy up to almost 90%. The accuracy decreases for generic reconstruction models trained using a collection of discharge data for multiple hoppers. Finally, we apply the CNN model to accurately reconstruct a hopper containing four distinct solid phases with a layered configuration to further demonstrate its utility.

42 ENGINEERING↗
DOE OSTI2021

A nonlinear elasto-plastic bond model for the discrete element modeling of woody biomass particles

This work presents an experiment-informed, semi-empirical, elasto-plastic bond model for the discrete element modeling of woody biomass particles. The model renders nonlinear plastic deformation of materials when subjected to compression/tension, bending, and twisting, essential for accurately simulating the behavior of biomass in comminution. The model is implemented in an open-source DEM package LIGGGHTS and is assessed in a number of verification tests. The model is applied to simulate the fracturing test of notched loblolly pine blocks. Reduced bond strength is prescribed to the weak planes along the growth rings. It is found that the square-root formulation for describing the elasto-plastic normal bond contact delivers the best agreement with the experimental data over its linear and quadratic counterparts. Furthermore, the model predicts a strong anisotropic behavior of wood blocks with regard to the grain orientations, indicating that the cutting angle should be an important attribute for the optimization of biomass comminution.

42 ENGINEERING↗
DOE OSTI2023

Particle size and shape effect of Crumbler® rotary shear-milled granular woody biomass on the performance of Acrison® screw feeder: A computational and experimental investigation

Physical experiments and discrete element model (DEM) simulations are conducted to evaluate particle characteristics and operation parameter effects on screw feeding performance for rotary shear-milled Douglas fir. Three performance metrics are used: mass flow rate, shaft driving torque, and specific energy consumption. The impact of particle size, particle size distribution (PSD), shaft rotational speed (rpm), and hopper dimensions on the performance are investigated. All employed performance metrics reveal the superior flowability of the 2-mm particles in contrast to the larger 6-mm counterpart. Remarkably, wider PSD results in poorer flowability than the two mono-sized particles, proving the flowability enhancement achieved by narrower PSD of the woody feedstock. More importantly, DEM simulations unveil PSD-induced degradation in flowability is attributed to mechanical interlocking and particle segregation effects. Furthermore, higher shaft rpm causes higher mass flow rate at the cost of higher specific energy consumption due to viscous dissipation and changes in flow pattern.

09 BIOMASS FUELS↗
DOE OSTI2011

Selective adsorption of atomic hydrogen on a h-BN thin film

The adsorption of atomic hydrogen on hexagonal boron nitride (h-BN) is studied using two element-specific spectroscopies, i.e., near-edge x-ray absorption fine structure (NEXAFS) spectroscopy and x-ray photoelectron spectroscopy (XPS). B K-edge NEXAFS spectra show a clear change in the energy region of the {pi}{sup *} band before and after reaction with atomic deuterium. On the other hand, N K-edge NEXAFS spectra show only a little change. B 1s XPS spectra show a distinct component at the low binding energy side of a main component, while N 1s XPS spectra show peak broadening at the high binding energy side. These experimental results are analyzed by the discrete variational X{alpha} method with a core-hole effect and are explained by a model in which hydrogen atoms are preferentially adsorbed on the B sites of h-BN. Based on the experimental and theoretical results, we propose a site-selective property of BN material on adsorption of atomic hydrogen.

36 MATERIALS SCIENCE↗
DOE OSTI2013

Dynamic and impact contact mechanics of geologic materials: Grain-scale experiments and modeling

High fidelity treatments of the generation and propagation of seismic waves in naturally occurring granular materials is becoming more practical given recent advancements in our ability to model complex particle shapes and their mechanical interaction. Of particular interest are the grain-scale processes that are activated by impact events and the characteristics of force transmission through grain contacts. To address this issue, we have developed a physics based approach that involves laboratory experiments to quantify the dynamic contact and impact behavior of granular materials and incorporation of the observed behavior indiscrete element models. The dynamic experiments do not involve particle damage and emphasis is placed on measured values of contact stiffness and frictional loss. The normal stiffness observed in dynamic contact experiments at low frequencies (e.g., 10 Hz) are shown to be in good agreement with quasistatic experiments on quartz sand. The results of impact experiments - which involve moderate to extensive levels of particle damage - are presented for several types of naturally occurring granular materials (several quartz sands, magnesite and calcium carbonate ooids). Implementation of the experimental findings in discrete element models is discussed and the results of impact simulations involving up to 5 Multiplication-Sign 105 grains are presented.

58 GEOSCIENCES↗
DOE OSTI2015

Free surface modeling in OWC chamber with parabolic side walls using 3D BEM

In this paper, BEM was used for free surface modeling in OWC chamber and out of it. Linear kinematic and dynamic boundary conditions were used for free surface out of OWC chamber and nonlinear forms were used for free surface in the chamber. These boundary conditions were discretized by finite differences method. Also, some thermodynamics relations were applied for trapped air behavior modeling in OWC chamber. Wave specifications in Chabahar region were used in modeling because these waves have an acceptable power for electricity generation. The results show a good agreement with results of other researches.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Electron-impact vibrational excitation of cyclopropane

We report a very detailed test of the ab initio discrete momentum representation (DMR) method of calculating vibrational excitation of polyatomic molecules by electron impact, by comparison of its results with an extensive set of experimental data, covering the entire range of scattering angles from 10{sup ∘} to 180{sup ∘} and electron energies from 0.4 to 20 eV. The DMR calculations were carried out by solving the two-channel Lippmann-Schwinger equation in the momentum space, and the interaction between the scattered electron and the target molecule was described by exact static-exchange potential corrected by a density functional theory (DFT) correlation-polarization interaction that models target’s response to the field of incoming electron. The theory is found to quantitatively reproduce the measured spectra for all normal modes, even at the difficult conditions of extreme angles and at low energies, and thus provides full understanding of the excitation mechanism. It is shown that the overlap of individual vibrational bands caused by limited experimental resolution and rotational excitation must be properly taken into account for correct comparison of experiment and theory. By doing so, an apparent discrepancy between published experimental data could be reconciled. A substantial cross section is found for excitation of the non-symmetric HCH twisting mode ν{sub 4} of A{sub 1}{sup ″} symmetry by the 5.5 eV A{sub 2}{sup ′} resonance, surprisingly because the currently accepted selection rules predict this process to be forbidden. The DMR theory shows that the excitation is caused by an incoming electron in an f-wave of A{sub 2}{sup ′} symmetry which causes excitation of the non-symmetric HCH twisting mode ν{sub 4} of the A{sub 1}{sup ″} symmetry and departs in p- and f-waves of A{sub 2}{sup ″} symmetry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

Microscopic study of {sup 6}He elastic scattering around the Coulomb barrier

We investigate {sup 6}He scattering on {sup 27}Al, {sup 58}Ni, {sup 120}Sn, and {sup 208}Pb in a microscopic version of the Continuum Discretized Coupled Channel (CDCC) method. We essentially focus on energies around the Coulomb barrier. The {sup 6}He nucleus is described by an antisymmetric 6-nucleon wave function, defined in the Resonating Group Method. The {sup 6}He continuum is simulated by square-integrable positive-energy states. The model does not depend on any adjustable parameter as it is based only on well known nucleon-target potentials. We show that experimental elastic cross sections are fairly well reproduced. The calculation suggests that breakup effects increase for high target masses. For a light system such as {sup 6}He+{sup 27}Al, breakup effects are small, and a single-channel approximation provides fair results.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Energy-conserving contact dynamics of nonspherical rigid-body particles

Understanding the contact dynamics of nonspherical particles is crucial for accurately modeling colloidal and granular systems where shape anisotropy dictates structural organization and transport properties. We here introduce an energy-conserving contact dynamics framework for arbitrary convex rigid-body particles by implementing vertex–boundary interactions in 2D and vertex–surface and edge–edge detection in 3D. The established formulation enables continuous force evaluation and prevents particle overlap while conserving total energy during translational and rotational motion. We demonstrate the framework’s stability and its utility to capture packing behavior, anisotropic diffusion, and equations of state of polygonal and polyhedral particles as examples. The framework establishes a robust and extensible foundation for investigating nonequilibrium dynamics of complex nonspherical particulate systems, enabling enhanced understanding of phenomena across spatiotemporal scales in self- and directed-assembly, granular flows, and hydrodynamics, potentially coupled with interactions that represent underlying physical mechanisms.

Discrete element method↗
DOE OSTI2021

Skyshine Calculations for a Large Spent Nuclear Fuel Storage Facility with SCALE 6.2.3

The SCALE code system developed at Oak Ridge National Laboratory includes state-of-the-art capabilities for radiation source term and radiation transport simulations that can be used in numerous applications, including dose rate analyses of complex consolidated interim storage facilities (CISFs). A licensed CISF could be used to store tens of thousands of tonnes of spent nuclear fuel discharged from commercial power reactors using various cask and storage pad designs. A CISF design must comply with the regulatory requirements provided in 10 CFR Part 72, including requirements related to annual dose limits applicable to real individuals located beyond the area controlled by the licensee. Therefore, calculating a dose to the public is a necessary part of the licensing process for the construction of a CISF. These calculations are very challenging because of the complexity of the CISF design and the low magnitude of dose rate at large distances from the facility. This paper describes detailed far-field dose rate calculations performed for a proposed CISF using MAVRIC, the Monte Carlo radiation shielding sequence in SCALE 6.2.3, with automated variance reduction based on discrete ordinates calculations. The method presented in this paper uses a detailed Monte Carlo radiation transport simulation in one step from source to dose rate. A series of independent simulations was made using the complete site geometry (all casks present), but with only one cask containing radiation sources to obtain the dose rate maps produced by each storage cask. The CISF dose rate map was obtained by adding the dose rate maps produced by the independent individual cask simulations. Ample volumes of air and soil extending beyond the location of interest for dose rate calculation were included in the calculation model to properly simulate important radiation attenuation and scattering events that affect far-field dose rates. Furthermore, a comprehensive sensitivity study is included in this paper to illustrate the importance of selecting appropriate air volume, mass density, and composition for CISF skyshine dose rate calculations. Dry soil and soil containing water were analyzed to determine their effects on groundshine radiation.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Optimizing training trajectories in variational autoencoders via latent Bayesian optimization approach *

Unsupervised and semi-supervised ML methods such as variational autoencoders (VAE) have become widely adopted across multiple areas of physics, chemistry, and materials sciences due to their capability in disentangling representations and ability to find latent manifolds for classification and/or regression of complex experimental data. Like other ML problems, VAEs require hyperparameter tuning, e.g. balancing the Kullback–Leibler and reconstruction terms. However, the training process and resulting manifold topology and connectivity depend not only on hyperparameters, but also their evolution during training. Because of the inefficiency of exhaustive search in a high-dimensional hyperparameter space for the expensive-to-train models, here we have explored a latent Bayesian optimization (zBO) approach for the hyperparameter trajectory optimization for the unsupervised and semi-supervised ML and demonstrated for joint-VAE with rotational invariances. We have demonstrated an application of this method for finding joint discrete and continuous rotationally invariant representations for modified national institute of standards and technology database (MNIST) and experimental data of a plasmonic nanoparticles material system. The performance of the proposed approach has been discussed extensively, where it allows for any high dimensional hyperparameter trajectory optimization of other ML models.

42 ENGINEERING↗
DOE OSTI2011

Effective radii of deuteron-induced reactions

The continuum-discretized coupled-channels method (CDCC) for exclusive reactions and the eikonal reaction theory (ERT) as an extension of CDCC to inclusive reactions are applied to deuteron-induced reactions. The CDCC result reproduces experimental data on the reaction cross section for d+{sup 58}Ni scattering at 200 MeV/nucleon, and ERT provides data on the neutron-stripping cross section for inclusive {sup 7}Li(d,n) reaction at 40 MeV. For deuteron-induced reactions at 200 MeV/nucleon, target-dependence of the reaction, elastic-breakup, nucleon-stripping, nucleon-removal, and complete- and incomplete-fusion cross sections is clearly explained by simple formulas. Accuracy of the Glauber model is also investigated.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2011

Elastic scattering and total reaction cross sections for the {sup 8}B, {sup 7}Be, and {sup 6}Li +{sup 12}C systems

Angular distributions for the elastic scattering of {sup 8}B, {sup 7}Be, and {sup 6}Li on a {sup 12}C target have been measured at E{sub lab}=25.8, 18.8, and 12.3 MeV, respectively. The analyses of these angular distributions have been performed in terms of the optical model using Woods-Saxon and double-folding type potentials. The effect of breakup in the elastic scattering of {sup 8}B+{sup 12}C is investigated by performing coupled-channels calculations with the continuum discretized coupled-channel method and cluster-model folding potentials. Total reaction cross sections were deduced from the elastic-scattering analysis and compared with published data on elastic scattering of other weakly and tightly bound projectiles on {sup 12}C, as a function of energy. With the exception of {sup 4}He and {sup 16}O, the data can be described using a universal function for the reduced cross sections.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Physical Layer Security at a Point-to-Point MIMO System With 1-Bit DACs and ADCs

Low-resolution converters reduce the power consumption of wireless devices and help enable large energy-efficient arrays. In this letter, we evaluate the secrecy rate for the extreme case of 1-bit single-user multiple-input multiple-output systems following two signaling approaches. In the first method, the system uses discrete signaling. In the second, the system combines Gaussian signaling with artificial noise alongside power allocation. In the first case, we give a closed-form expression of the secrecy rate and propose an algorithm that provides near-optimal results for low signal-to-noise ratios (SNRs). Further, in the second scenario, depending on the channel state information, we suggest different precoding plans and derive a secrecy rate lower bound. At low SNR, power is allocated to the main signal. Only at high SNR is optimal power allocation effective. Moreover, knowing the eavesdropper channel greatly improves the system’s secrecy.

1-bit digital-to-analog (DAC) and analog-to-digita↗
DOE OSTI2016

SU-C-207A-01: A Novel Maximum Likelihood Method for High-Resolution Proton Radiography/proton CT

Purpose: Multiple Coulomb scattering is the largest contributor to blurring in proton imaging. Here we tested a maximum likelihood least squares estimator (MLLSE) to improve the spatial resolution of proton radiography (pRad) and proton computed tomography (pCT). Methods: The object is discretized into voxels and the average relative stopping power through voxel columns defined from the source to the detector pixels is optimized such that it maximizes the likelihood of the proton energy loss. The length spent by individual protons in each column is calculated through an optimized cubic spline estimate. pRad images were first produced using Geant4 simulations. An anthropomorphic head phantom and the Catphan line-pair module for 3-D spatial resolution were studied and resulting images were analyzed. Both parallel and conical beam have been investigated for simulated pRad acquisition. Then, experimental data of a pediatric head phantom (CIRS) were acquired using a recently completed experimental pCT scanner. Specific filters were applied on proton angle and energy loss data to remove proton histories that underwent nuclear interactions. The MTF10% (lp/mm) was used to evaluate and compare spatial resolution. Results: Numerical simulations showed improvement in the pRad spatial resolution for the parallel (2.75 to 6.71 lp/cm) and conical beam (3.08 to 5.83 lp/cm) reconstructed with the MLLSE compared to averaging detector pixel signals. For full tomographic reconstruction, the improved pRad were used as input into a simultaneous algebraic reconstruction algorithm. The Catphan pCT reconstruction based on the MLLSE-enhanced projection showed spatial resolution improvement for the parallel (2.83 to 5.86 lp/cm) and conical beam (3.03 to 5.15 lp/cm). The anthropomorphic head pCT displayed important contrast gains in high-gradient regions. Experimental results also demonstrated significant improvement in spatial resolution of the pediatric head radiography. Conclusion: The proposed MLLSE shows promising potential to increase the spatial resolution (up to 244%) in proton imaging.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2026

Spatially Accelerated Winding Numbers for Curved Geometry

The generalized winding number (GWN) is a scalar field that supports robust containment queries on curved geometry, including non-watertight, overlapping, and nested boundary representations. While queries can be easily parallelized over samples, direct evaluation on parametric curves and surfaces remains costly for large and complex models. Fast, state-of-the-art GWN approaches leverage a spatial index to approximate the GWN, typically coupled with a Taylor expansion which approximates the GWN contribution for far clusters of geometric primitives. However, such methods operate only on discrete inputs such as triangle meshes and point clouds, and would introduce containment errors near boundaries if applied to curved input. We extend support for fast GWN evaluation over arbitrary collections of NURBS curves in 2D and trimmed NURBS patches in 3D via a Bounding Volume Hierarchy that stores efficiently precomputed moment data in the hierarchy nodes. When querying the hierarchy, approximations for far clusters are used alongside direct evaluation for nearby NURBS primitives, achieving sub-linear complexity while preserving the geometric features in the vicinity of the query point. Central to our performance improvements is an adaptive subdivision strategy for NURBS primitives during a preprocessing phase, creating better spatial partitions while retaining the same accuracy for containment decisions as a direct evaluation. We demonstrate the performance and accuracy of our approach across a large collection of 2D and 3D datasets.

Computer science↗