Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “particle 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,405 records · Page 78
All sourcesNASA NTRSDOE OSTI
NASA NTRS2023

Oxidation of Silicon Particles Suspended in Mullite Matrix With Dual-Mode Vacancy Diffusion of Oxygen

Mullite bond coats filled with silicon particles have recently been proposed for the next generation of environmental barrier coatings (EBCs). The sacrificial oxidation of the silicon particles eliminates some of the environmental durability issues associated with traditional EBC designs. To investigate this novel design approach, a mathematical formulation and resulting numerical method have been developed for simulating the oxidation of silicon particles embedded in a mullite matrix. The focus is on the oxidation of the silicon in a dry oxygen environment where the oxidant is presumably oxygen. Oxygen transport from the exterior environment through the mullite is assumed to occur via two lattice defect diffusion mechanisms. As a result, the rate of oxygen transport is dictated by the self-diffusivity of the lattice defects. The numerical method is implemented using the COMSOL Multiphysics® program (COMSOL, Inc.). The method is applied to simulate the oxidation of silicon particles in a mullite pellet in recent oxidation experiments. Numerical results compare favorably with the oxidation behavior observed in the oxidation experiments. The numerical method is also used to investigate the effects of various features of the mullite-silicon particle system on the oxidation behavior.

silicon oxidation↗
DOE OSTI2022

Hybrid-gyrokinetic simulations of low- n toroidal Alfvén eigenmodes using gKPSP

Here, we report a benchmark study of toroidal Alfvén eigenmode (TAE) simulation using the hybrid-gyrokinetic code GyroKinetic Plasma Simulation Program (gKPSP). A simulation capability for energetic particles based on the gyrokinetic δf method has been newly implemented in the gKPSP code. Benchmark simulations have been performed in both circular and realistic tokamak geometries. Good agreement has been found with previously reported results, demonstrating the new capability of the gKPSP code. We have investigated the effects of the distribution function on TAE stability by examining both isotropic and anisotropic slowing-down distributions of energetic particles. The slowing-down distribution produces a higher linear growth rate than a Maxwellian distribution, while a growth rate scan with the anisotropy shows an opposite trend. This can be attributed to competition between Landau damping and the linear drive, which are correlated with the fraction of resonant passing particles and their distribution in phase space.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Spatially Aware Linear Transformer (SAL-T) for Particle Jet Tagging

Transformers are very effective in capturing both global and local correlations within high-energy particle collisions, but they present deployment challenges in high-data-throughput environments, such as the CERN LHC. The quadratic complexity of transformer models demands substantial resources and increases latency during inference. In order to address these issues, we introduce the Spatially Aware Linear Transformer (SAL-T), a physics-inspired enhancement of the linformer architecture that maintains linear attention. Our method incorporates spatially aware partitioning of particles based on kinematic features, thereby computing attention between regions of physical significance. Additionally, we employ convolutional layers to capture local correlations, informed by insights from jet physics. In addition to outperforming the standard linformer in jet classification tasks, SAL-T also achieves classification results comparable to full-attention transformers, while using considerably fewer resources with lower latency during inference. Experiments on a generic point cloud classification dataset (ModelNet10) further confirm this trend. Our code is available at https://github.com/aaronw5/SAL-T4HEP.

Wang, Aaron [Illinois U., Chicago] (ORCID:00000003↗
NASA NTRS1986

Kinetic theory based new upwind methods for inviscid compressible flows

Two new upwind methods called the Kinetic Numerical Method (KNM) and the Kinetic Flux Vector Splitting (KFVS) method for the solution of the Euler equations have been presented. Both of these methods can be regarded as some suitable moments of an upwind scheme for the solution of the Boltzmann equation provided the distribution function is Maxwellian. This moment-method strategy leads to a unification of the Riemann approach and the pseudo-particle approach used earlier in the development of upwind methods for the Euler equations. A very important aspect of the moment-method strategy is that the new upwind methods satisfy the entropy condition because of the Boltzmann H-Theorem and suggest a possible way of extending the Total Variation Diminishing (TVD) principle within the framework of the H-Theorem. The ability of these methods in obtaining accurate wiggle-free solution is demonstrated by applying them to two test problems.

Deshpande, S. M.↗
DOE OSTI2022

Solar receivers and methods for capturing solar energy

Thermal receivers, systems, and methods are disclosed that efficiently capture concentrated solar energy into a plurality of heat absorption bodies for conversion into thermal energy. In an embodiment, the thermal receivers, systems, and methods enable simultaneous electricity conversion and thermal energy capture. The receiver design enables a high penetration of concentrated sunlight deep into the thermal receiver to increase light trapping and reduce thermal losses. The thermal receiver is integrated with a photovoltaic (PV) receiver platform that converts some of the incident light to electricity while passing the remaining light to the thermal receiver. In another embodiment, other thermal receivers, systems, and methods are disclosed that efficiently capture concentrated solar energy into a sheet of falling particles. In an embodiment, the thermal receivers, systems, and methods enable simultaneous electricity conversion and thermal energy capture.

Ho, Clifford K.↗
DOE OSTI2025

Impacts of Biomass Feedstock Pre-Processing on Heat and Mass Transfer During Pyrolysis Using X-Ray Computed Tomography and Multiscale Modeling

Knowledge of the transport properties of biomass particles such as porosity, tortuosity, and permeability is paramount for high-fidelity modeling of biomass pyrolysis due to the heat and mass transfer limitations imposed by particle microstructure. X-ray computed tomography (XCT) is a non-destructive imaging method that enables full 3D reconstructions of the biomass particle microstructure with high resolution, permitting direct calculation of porosity, tortuosity, and permeability from real particle geometries. In this study, XCT imaging revealed the 3D microstructures of particles and chars from pyrolytic conversion of cylindrically cut or milled/pelletized loblolly pine samples. The porosity, tortuosity, and permeability were calculated directly from the XCT geometries via open-source microstructural analysis tool MATBOX+TauFactor (https://github.com/NREL/MATBOX_Microstructure_analysis_toolbox) and computational fluid dynamics (CFD) simulations using our solver, Mesoflow (https://github.com/NREL/mesoflow). These properties were used in a reactor scale model developed in COMSOL of the single particle reactor at NREL to investigate the impact of feedstock pre-processing on biomass conversion during pyrolysis with rigorous experimental validation.

biomass↗
NASA NTRS1976

The partially averaged field approach to cosmic ray diffusion

The kinetic equation for particles interacting with turbulent fluctuations is derived by a new nonlinear technique which successfully corrects the difficulties associated with quasilinear theory. In this new method the effects of the fluctuations are evaluated along particle orbits which themselves include the effects of a statistically averaged subset of the possible configurations of the turbulence. The new method is illustrated by calculating the pitch angle diffusion coefficient D sub Mu Mu for particles interacting with slab model magnetic turbulence, i.e., magnetic fluctuations linearly polarized transverse to a mean magnetic field. Results are compared with those of quasilinear theory and also with those of Monte Carlo calculations. The major effect of the nonlinear treatment in this illustration is the determination of D sub Mu Mu in the vicinity of 90 deg pitch angles where quasilinear theory breaks down. The spatial diffusion coefficient parallel to a mean magnetic field is evaluated using D sub Mu Mu as calculated by this technique. It is argued that the partially averaged field method is not limited to small amplitude fluctuating fields and is hence not a perturbation theory.

Jones, F. C.↗
NASA NTRS1990

Modeling the dust size distribution in comets with dust fragmentation

A hydrodynamic model was developed of a spherically symmetric dusty gas flow in a cometary atmosphere assuming a single fluid, inviscid, perfect gas. The hydrodynamics for gas and dust, which involves the gas drag force (momentum transfer), heat exchange between gas and dust, photodissociation energy for H2O gas, and radiative heating and cooling terms for dust particles, are solved using the Gear method for stiff, coupled differential equations. Calculations were done with a dust size distribution for radii alpha = 0.01 micron to 10 cm with densities variable with the size. A nucleus size of 4.0 km radius with a density of 0.5 g/cu cm and a total dust-to-gas mass ratio chi = 1 was adopted. There are indications from in situ observations that dust particle fragment into smaller ones. Fragmentation of dust particles was incorporated into the model. This is done by adding source and sink terms in the continuity equations for the dust. Lifetimes for the decay of dust particles were assumed as a function of particle size. It is also assumed that dust particles always fragment only into the next smaller size.

Konno, Ichishiro↗
NASA NTRSDate not supplied

AM Powder Flowability Capabilities at NASA Marshall Space Flight Center (MSFC)

Additive Manufacturing (AM) powder flowability is critical to metal 3D printing, because the more fluid the powder is, the better powder spreads. The NASA MSFC’s Contamination Control Team (CCT) studies density, flow, particle size distribution, and morphology for programs such as Moon-to-Mars Planetary Autonomous Construction Technology (MMPACT), ASTM Proficiency Testing for AM and Powder Metallurgy, and MSFC’s AM team. Examining and measuring powder characteristics is essential for improving flowability of AM powder and ensuring lot-to-lot consistency, which will help to prevent defects in manufactured parts. The CCT uses several methods to characterize AM powder. Optical particle size distribution and morphology analyses determine particle parameters including size, circularity, convexity, and dimensions. The Carney and Hall Flowmeter Funnel measures time it takes powder to flow through a funnel to compare relative flowability (free-flowing and non-free-flowing). In addition to this, the device utilizes a density cup which determines apparent density. Lastly, the Revolution Powder Analyzer measures dynamic powder flowability and behavior over time via digital imaging. Using the CCT’s lab, the team has produced ASTM Powder Proficiency Testing results with flow rates, apparent density, particle size distribution, and tap density on titanium-based and nickel-based powder and compared data with other companies/labs. Morphology data was also collected for various regolith simulant powder for the MMPACT program. Lastly, the CCT used morphology data to examine the flowability of Inconel 718 powder for Selective Laser Melting (SLM) and Directed Energy Deposition (DED) printers at MSFC. The CCT’s essential work in studying and enabling AM powder characterization has helped and will continue to help study new and refined powders in the AM industry.

Additive Manufacturing↗
DOE OSTI2025

Extending TOUGH + HYDRATE with a parallel particle transport simulator: numerical investigation of sand production during gas production from hydrate deposits

A new parallel code for simulating particle transport in porous media is integrated with the TOUGH + HYDRATE simulator to investigate sand production associated with gas production from unconsolidated gas hydrate-bearing sediments (HBS). Here, the parallel coupled simulator is named THMPT and uses the integral finite difference method to describe the Darcian and non-Darcian flow of fluids and heat transport, the finite element method to describe the associated geomechanical changes, and the discrete element method to track the trajectory of individual sand particles within the HBS. The THMPT simulator is written in Fortran, incorporates multiple optimized algorithms, and can comprehensively address the coupled flow, thermal, chemical, geomechanical, and particle transport processes that characterize the system behaviors during gas production from HBS. The simulator can capture all processes involved in sand particle transport in porous media, including sand detachment, collision, clogging (i.e., bridging), and migration. A benchmark case study of sand production in the course of depressurization-induced gas production from a representative HBS reveals various distinct microscopic particle migration mechanisms and the adverse impact of sand particle detachment, transport, and clogging. The numerical investigation also examines the effect of bottomhole pressure on mitigating sand production. The simulation results indicate that sand clogging near the wellbore significantly reduces permeability, decreasing gas production by at least 50%. Lastly, the efficiency of gravel packing in mitigating sand production is numerically evaluated, revealing that the structure of the porous media appears to profoundly influence the macroscopic motion behavior of sand particles and sand clogging characteristics.

discrete element method↗
DOE OSTI2024

Low Energy Neutron-induced Charged-particle (Z) (LENZ) instrument development with a focus on pulse shape discrimination for low-energy charged particles

To study neutron-induced charged-particle reactions with high precision, the Low Energy Neutron-induced Charged-particle (Z) instrument (LENZ) was developed at the Los Alamos Neutron Science Center. For the interest of measuring (n,p) and (n, α) reactions simultaneously, Pulse Shape Discrimination methods were investigated to identify different charged particles in the energy range of 3 - 20 MeV and improve signal-to-background ratios using Double-sided Silicon Strip Detectors and waveform digitizers. The risetime and pulse shape properties of detected charged particles were characterized for various silicon detector’s thickness with different orientations. During the post-processing of waveforms, we implemented different digital filters for effective particle identifications and improved energy- and timing- resolutions. We validated the optimized digital filters and pulse shape analyses, via measurements with 228 Th and 229 Th calibration sources, proton-induced reactions on a 7 LiF target, and neutron-induced reactions on CH 2 , Ta 2 O 5 , 58 Ni, and 6 LiF/ 59 Ni targets at the time-of-flight facility, LANSCE. In conclusion, the summary of effective thresholds and Figure of Merits on separating different charged particles is reported.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS2003

Numerical Simulation of One- and Two-Phase Flows in Propulsion Systems

Four subprojects were conducted using analytical methods, numerical simulation and experimental tests: (A) Shock wave mitigation by spike-shaped blunt bodies with application for the purpose of drag, lift and longitudinal momentum optimization. The main result in this subproject is: application of a single needle against a supersonic flow provides higher benefits for blunt body drag reduction and heat transfer to the body than the application of multiple needles. (B) Solid particles, liquid and air jet injection through the front of a blunt body against a supersonic flow. In this case, the research conducted and analysis of multiple previous investigations in this area have shown essential benefits and preferable application of solid particle injection. (C) Comparison of different methods of fuel injection into supersonic duct flows. Preliminary numerical simulations and theoretical analysis show promising results for Telescope-shaped inlet applications in SCRAMJET; and (D) Development of an acoustic source location method for different applications including propulsion systems.

Gilinsky, Mikhail↗
NASA NTRS2004

MESTRN: A Deterministic Meson-Muon Transport Code for Space Radiation

A safe and efficient exploration of space requires an understanding of space radiations, so that human life and sensitive equipment can be protected. On the way to these sensitive sites, the radiation fields are modified in both quality and quantity. Many of these modifications are thought to be due to the production of pions and muons in the interactions between the radiation and intervening matter. A method used to predict the effects of the presence of these particles on the transport of radiation through materials is developed. This method was then used to develop software, which was used to calculate the fluxes of pions and muons after the transport of a cosmic ray spectrum through aluminum and water. Software descriptions are given in the appendices.

Blattnig, Steve R.↗
DOE OSTI2024

Studying fast-ion populations using oscillations in solid-state neutral-particle analyzer signal and neutral-beam injection power

A method for determining the fast-ion population density in magnetically confined plasmas as a function of pitch-radius, (λ, R), using a solid-state neutral-particle analyzer (ssNPA) signal and neutral-beam injection (NBI) power-output data has been developed. Oscillations in the NBI power output are replicated only in the active part of the ssNPA signal, allowing this to be separated from the passive and background signals, which usually complicate data from this diagnostic. Results obtained using this method are compared with those from standard techniques using data from the Mega-Amp Spherical Tokamak Upgrade spherical tokamak.

Instruments & Instrumentation↗
NASA NTRSDate not supplied

AM Powder Flowability Capabilities at NASA Marshall Space Flight Center (MSFC)

Additive Manufacturing (AM) powder flowability is critical to metal 3D printing, because the more fluid the powder is, the better powder spreads. AM Powder Flowability Capabilities at NASA MSFC studies density, flow, particle size distribution, and morphology for programs such as Moon-to-Mars Planetary Autonomous Construction Technology (MMPACT), ASTM Proficiency Testing for AM and Powder Metallurgy, and MSFC’s AM team. Examining and measuring powder characteristics is essential for improving flowability of AM powder and ensuring lot-to-lot consistency, which will help to prevent defects in manufactured parts. The Contamination Control Team (CCT) uses several methods to characterize AM powder. Optical particle size distribution and morphology analyses determine particle parameters including size, circularity, convexity, and dimensions. The Carney and Hall Flowmeter Funnel measures time it takes powder to flow through a funnel to compare relative flowability (free-flowing and non-free-flowing). In addition to this, the device utilizes a density cup which determines apparent density. Lastly, the Revolution Powder Analyzer measures dynamic powder flowability and behavior over time via digital imaging. Using the CCT’s lab, the team has produced ASTM Powder Proficiency Testing results with flow rates, apparent density, particle size distribution, and tap density on titanium-based and nickel-based powder and compared data with other companies/labs. Morphology data was also collected for various regolith simulant powder for the MMPACT program. Lastly, the CCT used morphology data to examine the flowability of Inconel 718 powder for Selective Laser Melting (SLM) and Directed Energy Deposition (DED) printers at MSFC. The CCT’s essential work in studying and enabling AM powder characterization has helped and will continue to help study new and refined powders in the AM industry.

Additive Manufacturing↗
DOE OSTI2025

Identification of uranium oxidation states using oxygen K-edge scanning transmission X-ray microscopy

The field of nuclear forensics is growing in importance, and the increasing capabilities at synchrotron radiation light sources enable non-destructive characterization of oxide particles with better spatial, compositional, and oxidation state speciation resolution than ever before. Here, uranium oxide particles derived from multiple wet chemical processing methods were examined using a scanning transmission X-ray microscope (STXM), and a weakly-supervised method was developed to automatically analyze the collected data. Multiple uranium oxidation states were observed and quantified within and between samples, yielding information about differences between particles produced via the various processing routes.

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

A large expert-curated cryo-EM image dataset for machine learning protein particle picking

Cryo-electron microscopy (cryo-EM) is a powerful technique for determining the structures of biological macromolecular complexes. Picking single-protein particles from cryo-EM micrographs is a crucial step in reconstructing protein structures. However, the widely used template-based particle picking process is labor-intensive and time-consuming. Though machine learning and artificial intelligence (AI) based particle picking can potentially automate the process, its development is hindered by lack of large, high-quality labelled training data. To address this bottleneck, we present CryoPPP, a large, diverse, expert-curated cryo-EM image dataset for protein particle picking and analysis. It consists of labelled cryo-EM micrographs (images) of 34 representative protein datasets selected from the Electron Microscopy Public Image Archive (EMPIAR). The dataset is 2.6 terabytes and includes 9,893 high-resolution micrographs with labelled protein particle coordinates. The labelling process was rigorously validated through 2D particle class validation and 3D density map validation with the gold standard. The dataset is expected to greatly facilitate the development of both AI and classical methods for automated cryo-EM protein particle picking.

42 ENGINEERING↗
DOE OSTI2026

Direct cross section measurement of 102 Pd ⁢(𝛾,𝑝) and 102 Pd ⁢(𝛾,𝛼) for the astrophysical 𝑝 process

Background: A handful of neutron-deficient stable nuclei, known as the “p nuclei,” cannot be produced through astrophysical neutron capture processes. Instead, some of these nuclei are proposed to be produced by 𝛾-induced reactions on existing r- and s-process seeds. The specific astrophysical site or sites are not yet identified, however, with uncertainties in the cross sections of these 𝛾-induced reactions playing a role. Databases of reaction rates for astrophysical simulations often rely on theoretical statistical model calculations, such as Hauser-Feshbach, for rates where no experimental information is known. However, reasonable variations in the choice of parametrizations of various nuclear properties can create order-of-magnitude variations in the final predicted cross sections and reaction rates, which are then propagated through the models to the predicted final abundances. Purpose: To better constrain these statistical model calculations and ultimately reduce the uncertainties from the nuclear physics on our understanding of the p nuclei, a measurement of the cross sections of 𝛾-induced reactions on the p-nucleus 102 Pd was undertaken. This work represents the first measurement of its kind, using segmented silicon detectors to measure prompt charged particle emission from 𝛾-induced reactions. Methods: Quasimonoenergetic gamma beams from the High Intensity 𝛾 Source facility bombarded an enriched 102 Pd target. A segmented silicon array was arranged to detect the particles emitted from (𝛾,𝑝) and (𝛾,𝛼) reactions. Results: Reaction cross sections were deduced at multiple 𝛾-beam energies between 10 and 19 MeV, and compared to statistical model calculations using talys-1.96. The 102 Pd ⁢(𝛾,𝑝)⁢ 101 Rh reaction cross section was reasonably well reproduced by a subset of photon strength functions and level densities, though the strength to the ground state of 101 Rh was underestimated at higher incident gamma energies. The 102 Pd ⁢(𝛾,𝛼)⁢ 98 Ru was in general overpredicted by the various alpha-nucleus optical model potentials. Conclusions: While the theoretical cross sections used to model the (𝛾,𝑝) reactions for the p process may be reasonable, a more careful approach is needed in the case of (𝛾,𝛼). Further work to probe gamma-induced reaction cross sections at and near the p nuclei is warranted.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗