Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Particle representation”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Simulating energetic particle losses in JET plasmas with a reverse integration biasing scheme

We report that an integrated energetic particle transport model has been constructed in JET plasmas constrained by experimental fast ion loss measurements. The model incorporates a synthetic fast ion loss detector identical to JET’s thin-foil Faraday cup fast ion loss detector array. The loss model combines analyses from the TRANSP and ORBIT-kick codes with enhanced features for producing the synthetic diagnostic. Extensions to the ORBIT code framework allow a full-orbit representation within the vacuum region that can map particles directly to an installed detector geometry. Since synthetic fast ion loss detectors are plagued by weak loss statistics, a novel reverse integration biasing scheme has been implemented to boost computational efficiency. The model is validated against experimental loss measurements induced by long-lived kink modes and is found to be in good agreement. This confirms the development of a fully integrated transport/loss model which can be quantitatively verified against experiment allowing for future validation and predictive studies. The model is particularly useful for more complicated plasma scenarios that involve multiple fast ion species such as JET’s 2021 DT-campaign.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Effects of the method of active site characterization for determining structure-sensitivity in Ag-catalyzed ethylene epoxidation

We report catalysts have been prepared on a low surface area α-Al 2 O 3 support used commercially for previous generation olefin epoxidation catalysts. Prescreening of the low surface area alumina (0.73 m 2 /g) indicated the absence of acid catalyzed isomerization of EO at an evaluation temperature of 210 °C. A 0.1 wt% Ag base material synthesized by incipient wetness impregnation of AgNO 3 was used as a base material for electroless deposition (ED) of additional Ag to increase particle sizes by controlled reduction of Ag + directly onto the preexisting Ag surface to form weight loadings between 0.3 and 5.0 wt% metal. A 12 wt% Ag/α-Al 2 O 3 using Ag 2 C 2 O 4 as the Ag precursor was also prepared to compare performance of the ED samples with a catalyst more typical of industrial formulations. Characterization by SEM, STEM, and hydrogen titration of oxygen precovered Ag characterized before and after catalytic evaluation indicated that microscopy is required to accurately represent distributions of Ag particle sizes, but H 2 titration of O-precovered Ag gives the best representation of active sites since it directly counts the number of Ag surface sites. Larger particles >100 nm are resistant to both Ag sintering and carbon foulant; TOF values were relatively insensitive to particle size with only a 2.2 × difference between the best and worst performing samples. Selectivity, which is not a function of TOF, shows the most significant structure sensitivity effect where particle sizes follow the trend 67 nm ≈ 92 nm (58% EO) < 157 nm (67% EO) < 211–542 nm (73% EO). The lower EO selectivities were also correlated with increased fouling for the smaller Ag sizes, suggesting that more strongly bound EO precursor(s) leads to combustion and CO 2 /H 2 O formation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1985

Numerical techniques for large cosmological N-body simulations

Techniques for carrying out large N-body simulations of the gravitational evolution of clustering in the fundamental cube of an infinite periodic universe are described and compared. The accuracy of the forces derived from several commonly used particle mesh schemes is examined, showing how submesh resolution can be achieved by including short-range forces between particles by direct summation techniques. The time integration of the equations of motion is discussed, and the accuracy of the codes for various choices of 'time' variable and time step is tested by considering energy conservation as well as by direct analysis of particle trajectories. Methods for generating initial particle positions and velocities corresponding to a growing mode representation of a specified power spectrum of linear density fluctuations are described. The effects of force resolution are studied and different simulation schemes are compared. An algorithm is implemented for generating initial conditions by varying the number of particles, the initial amplitude of density fluctuations, and the initial peculiar velocity field.

Efstathiou, G.↗
NASA NTRS1993

The giant impact produced a precipitated Moon

The author's current simulations of Giant Impacts on the protoearth show the development of large hot rock vapor atmospheres. The Balbus-Hawley mechanism will pump mass and angular momentum outwards in the equatorial plane; upon cooling and expansion the rock vapor will condense refractory material beyond the Roche distance, where it is available for lunar formation. During the last seven years, the author together with several colleagues has carried out a series of numerical investigations of the Giant Impact theory for the origin of the Moon. These involved three-dimensional simulations of the impact and its aftermath using Smooth Particle Hydrodynamics (SPH), in which the matter in the system is divided into discrete particles whose motions and internal energies are determined as a result of the imposed initial conditions. Densities and pressures are determined from the combined overlaps of the particles, which have a bell-shaped density distribution characterized by a smoothing length. In the original series of runs all particle masses and smoothing lengths had the same values; the matter in the colliding bodies consisted of initial iron cores and rock (dunite) mantles. Each of 41 runs used 3,008 particles, took several weeks of continuous computation, and gave fairly good representations of the ultimate state of the post-collision body or bodies but at best crude and qualitative information about individual particles in orbit. During the last two years an improved SPH program was used in which the masses and smoothing lengths of the particles are variable, and the intent of the current series of computations is to investigate the behavior of the matter exterior to the main parts of the body or bodies subsequent to the collisions. These runs are taking times comparable to a year of continuous computation in each case; they use 10,000 particles with 5,000 particles in the target and 5,000 in the impactor, and the particles thus have variable masses and smoothing lengths (the latter are dynamically adjusted so that a particle typically overlaps a few tens of its neighbors). Since the matter in the impactor provides the majority of the mass left in orbit after the collision, and since the masses of the particles that originated in the impactor are smaller than those in the target, the mass resolution in the exterior parts of the problem is greatly improved and the exterior particles properly simulate atmospheres in hydrostatic equilibrium.

Cameron, A. G. W.↗
NASA NTRS2016

Aerosol Complexity and Implications for Predictability and Short-Term Forecasting

There are clear NWP and climate impacts from including aerosol radiative and cloud interactions. Changes in dynamics and cloud fields affect aerosol lifecycle, plume height, long-range transport, overall forcing of the climate system, etc. Inclusion of aerosols in NWP systems has benefit to surface field biases (e.g., T2m, U10m). Including aerosol affects has impact on analysis increments and can have statistically significant impacts on, e.g., tropical cyclogenesis. Above points are made especially with respect to aerosol radiative interactions, but aerosol-cloud interaction is a bigger signal on the global system. Many of these impacts are realized even in models with relatively simple (bulk) aerosol schemes (approx.10 -20 tracers). Simple schemes though imply simple representation of aerosol absorption and importantly for aerosol-cloud interaction particle-size distribution. Even so, more complex schemes exhibit a lot of diversity between different models, with issues such as size selection both for emitted particles and for modes. Prospects for complex sectional schemes to tune modal (and even bulk) schemes toward better selection of size representation. I think this is a ripe topic for more research -Systematic documentation of benefits of no vs. climatological vs. interactive (direct and then direct+indirect) aerosols. Document aerosol impact on analysis increments, inclusion in NWP data assimilation operator -Further refinement of baseline assumptions in model design (e.g., absorption, particle size distribution). Did not get into model resolution and interplay of other physical processes with aerosols (e.g., moist physics, obviously important), chemistry

Predictability↗
DOE OSTI2020

The O(N) S-matrix monolith

We consider the scattering matrices of massive quantum field theories with no bound states and a global O(N) symmetry in two spacetime dimensions. In particular we explore the space of two-to-two S-matrices of particles of mass m transforming in the vector representation as restricted by the general conditions of unitarity, crossing, analyticity and O(N) symmetry. We found a rich structure in that space by using convex maximization and in particular its convex dual minimization problem. At the boundary of the allowed space special geometric points such as vertices were found to correspond to integrable models. The dual convex minimization problem provides a novel and useful approach to the problem allowing, for example, to prove that generically the S-matrices so obtained saturate unitarity and, insome cases, that they are at vertices of the allowed space.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Graviton topology

Over the past three decades, it has been shown that discrete and continuous media can support topologically nontrivial waves. Recently, it was shown that the same is true of the vacuum, in particular, right (R) and left (L) circularly polarized photons are topologically nontrivial. Here, we study the topology of another class of massless particles, namely gravitons. We show that the collection of all gravitons forms a topologically trivial vector bundle over the lightcone, allowing us to construct a globally smooth basis for gravitons. The graviton bundle also has a natural geometric splitting into two topologically nontrivial subbundles, consisting of the R and L gravitons. The R and L gravitons are unitary irreducible bundle representations of the Poincaré group, and are thus elementary particles; their topology is characterized by the Chern numbers ∓4. This nontrivial topology obstructs the splitting of graviton angular momentum into spin and orbital angular momentum.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS2004

Inner Radiation Belt Representation of the Energetic Electron Environment: Model and Data Synthesis Using the Salammbo Radiation Belt Transport Code and Los Alamos Geosynchronous and GPS Energetic Particle Data

The highly energetic electron environment in the inner magnetosphere (GEO inward) has received a lot of research attention in resent years, as the dynamics of relativistic electron acceleration and transport are not yet fully understood. These electrons can cause deep dielectric charging in any space hardware in the MEO to GEO region. We use a new and novel approach to obtain a global representation of the inner magnetospheric energetic electron environment, which can reproduce the absolute environment (flux) for any spacecraft orbit in that region to within a factor of 2 for the energy range of 100 KeV to 5 MeV electrons, for any levels of magnetospheric activity. We combine the extensive set of inner magnetospheric energetic electron observations available at Los Alamos with the physics based Salammbo transport code, using the data assimilation technique of "nudging". This in effect input in-situ data into the code and allows the diffusion mechanisms in the code to interpolate the data into regions and times of no data availability. We present here details of the methods used, both in the data assimilation process and in the necessary inter-calibration of the input data used. We will present sample runs of the model/data code and compare the results to test spacecraft data not used in the data assimilation process.

Friedel, R. H. W.↗
DOE OSTI2025

Subject-specific modeling framework for particle deposition using computational fluid dynamics

Quantifying particle deposition and dose in the respiratory tract requires a physiologically realistic representation and reproducible computational workflows. However, existing modeling frameworks, such as the International Commission on Radiological Protection (ICRP) compartmental models and the Multiple Path Particle Dosimetry (MPPD) tool, lack detailed deposition profiles and subject-specific capabilities. The combination of advances in computer vision algorithms applied to the respiratory tract and Computational Fluid and Particle Dynamics (CFPD) allows high-fidelity simulations of particle behavior in anatomically accurate geometries derived from individual CT scans. The segmentation, preprocessing, and file preparation task for a CFPD simulation was often time-consuming, and no prior studies to-date have yet presented a fully automated framework. This work presents a fully automated workflow to obtain individualized particle deposition profiles in the human respiratory tract. The pipeline starts with segmenting upper and lower airway geometries using morphological and deep learning-based methods, generating three-dimensional (3D) models from CT imaging data. Next, a series of algorithms are presented to quality check and prepare the 3D geometry for a CFD or CFPD simulation. The preprocessing step includes correcting geometric artifacts, enforcing a physically consistent mesh, and automatically identifying and capping multiple outlets, which is required for CFD/CFPD simulations. These processed models are then input into open-source (OpenFOAM) or commercial (StarCCM+) CFD solvers, where flow and transient particle transport equations — including turbulence and particle–wall interactions are solved under realistic breathing conditions. Finally, the resulting particle deposition profiles can be integrated with Monte Carlo radiation transport codes and state-of-the-art computational phantoms to assess organ-specific absorbed doses in scenarios of radioactive aerosol inhalation. The presented work streamlines respiratory tract segmentation, preprocessing for CFD/CFPD simulations, and integration with dose assessment workflows, reducing manual intervention and improving access to high-fidelity, subject-specific modeling. The high precision in predicted particle deposition and dose distributions can improve personalized treatment strategies in respiratory medicine and refine dose estimates for radiation protection.

AI↗
DOE OSTI2024

Anti-symmetric and positivity preserving formulation of a spectral method for Vlasov-Poisson equations

We analyze the anti-symmetric properties of a spectral discretization for the one-dimensional Vlasov-Poisson equations. The discretization is based on a spectral expansion in velocity with the symmetrically weighted Hermite basis functions, central finite differencing in space, and an implicit Runge Kutta integrator in time. The proposed discretization preserves the anti-symmetric structure of the advection operator in the Vlasov equation, resulting in a stable numerical method. We apply such discretization to two formulations: the canonical Vlasov-Poisson equations and their continuously transformed square-root representation. The latter preserves the positivity of the particle distribution function. We derive analytically the conservation properties of both formulations, including particle number, momentum, and energy, which are verified numerically on the following benchmark problems: manufactured solution, linear and nonlinear Landau damping, two-stream instability, bump-on-tail instability, and ion-acoustic wave.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Coarse-grained methods for heterogeneous vesicles with phase-separated domains: Elastic mechanics of shape fluctuations, plate compression, and channel insertion

We develop coarse-grained particle approaches for studying the elastic mechanics of vesicles with heterogeneous membranes having phase-separated domains. We perform simulations both of passive shape fluctuations and of active systems where vesicles are subjected to compression between two plates or subjected to insertion into narrow channels. Analysis methods are developed for mapping particle configurations to continuum fields with spherical harmonics representations. Heterogeneous vesicles are found to exhibit rich behaviors where the heterogeneity can amplify surface two-point correlations, reduce resistance during compression, and augment vesicle transport times in channels. The developed methods provide general approaches for characterizing the mechanics of coarse-grained heterogeneous systems taking into account the roles of thermal fluctuations, geometry, and phase separation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Employing ternary fission of Pu 242 as a probe of very neutron-rich matter

Detailed assessments of the ability of recent theoretical approaches to modeling existing experimental data for ternary fission confirm earlier indications that the dominant mode of cluster formation in ternary fission is clusterization in very neutron rich, very low density, essentially chemically equilibrated, nucleonic matter. An extended study and comparison of these approaches applied to ternary fission yields in the thermal neutron induced reaction 241 Pu(n th ,f) has been undertaken to refine the characterization of the source matter. Here, the resonance gas approximation has been improved taking in-medium effects on the binding energies into account. A temperature of 1.29 MeV, density of 6.7 × 10 -5 nucleons/fm 3 and proton fraction Y p = 0.035 are found to provide a good representation of yields of the ternary emitted light particles and clusters. In particular, results for Z = 1 and 2 isotopes are presented. Isotopes with larger Z are discussed, and the roles of medium and continuum effects, even at very low density are illustrated.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1989

Scattering by fluffy grains

There are indications (Greenberg et al., 1988), that fluffy (i.e., porous) particles are responsible for the observed 3.4 and 10 micron emissions of comet Halley. The absorption characteristics of small particles both solid and fluffy are needed in order to explain the Halley emissions. How isolated small solid particles react to an external radiation field is well known - the Rayleigh approximation. How these same small particles emit when assembled as fluffy aggregates in another question. To what degree are the emission spectra of isolated and aggregated particles comparable. In order to quantify the assertion that fluffy particles produce the observed Halley infrared emission features, the authors are performing calculations to determine the effect of porosity on the absorption characteristics of aggregates of interstellar grain-type particles. The calculations are based on an integral representation of the scattered electromagnetic field. Results are given with application to comet Halley.

Hage, J. I.↗
NASA NTRS2015

Use of Existing CAD Models for Radiation Shielding Analysis

The utility of a radiation exposure analysis depends not only on the accuracy of the underlying particle transport code, but also on the accuracy of the geometric representations of both the vehicle used as radiation shielding mass and the phantom representation of the human form. The current NASA/Space Radiation Analysis Group (SRAG) process to determine crew radiation exposure in a vehicle design incorporates both output from an analytic High Z and Energy Particle Transport (HZETRN) code and the properties (i.e., material thicknesses) of a previously processed drawing. This geometry pre-process can be time-consuming, and the results are less accurate than those determined using a Monte Carlo-based particle transport code. The current work aims to improve this process. Although several Monte Carlo programs (FLUKA, Geant4) are readily available, most use an internal geometry engine. The lack of an interface with the standard CAD formats used by the vehicle designers limits the ability of the user to communicate complex geometries. Translation of native CAD drawings into a format readable by these transport programs is time consuming and prone to error. The Direct Accelerated Geometry -United (DAGU) project is intended to provide an interface between the native vehicle or phantom CAD geometry and multiple particle transport codes to minimize problem setup, computing time and analysis error.

Lee, K. T.↗
NASA NTRS1989

Three dimensional discrete particle simulation of an AOTV

A discrete Particle Simulation method, recently formulated by Baganoff, is discussed in the context of its application to the simulation of the flow field about the Aeroassisted Flight Experiment (AFE). As a basis for discussion the current algorithm is first described. Because of the use of a cubic Cartesian mesh, the representation of the geometry is different than that of other particle methods and an algorithm for its generation is discussed. The method is applied to test problems and then to the AFE calculation with the use of 9.52 million particles and 432,000 cells.

Feiereisen, William J.↗
DOE OSTI2023

Variational Neural-Network Ansatz for Continuum Quantum Field Theory

Physicists dating back to Feynman have lamented the difficulties of applying the variational principle to quantum field theories. In nonrelativistic quantum field theories, the challenge is to parametrize and optimize over the infinitely many n-particle wave functions comprising the state’s Fock-space representation. Here we approach this problem by introducing neural-network quantum field states, a deep learning ansatz that enables application of the variational principle to nonrelativistic quantum field theories in the continuum. Our ansatz uses the Deep Sets neural network architecture to simultaneously parametrize all of the n-particle wave functions comprising a quantum field state. We employ our ansatz to approximate ground states of various field theories, including an inhomogeneous system and a system with long-range interactions, thus demonstrating a powerful new tool for probing quantum field theories.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

A 1D Model for Nucleation of Ice From Aerosol Particles: An Application to a Mixed‐Phase Arctic Stratus Cloud Layer

Abstract Mixed‐phase clouds (MPCs) have been identified as significant contributors to uncertainties in climate projections, attributable to model representation of processes controlling the formation and loss of supercooled water droplets and ice particles from the atmosphere. Arctic MPCs are commonly widespread and long‐lived, with sustained ice crystal formation processes that challenge current understanding. This study examines the ice‐nucleating particle (INP) reservoir dynamics governing immersion‐mode heterogeneous freezing in an observed case of Arctic MPCs using a simplified 1D aerosol‐cloud model. The model setup includes prescribed dynamical forcings and thermodynamic profiles, and represents INPs as multicomponent and polydisperse particle size distributions. Diagnostic and prognostic approaches to immersion freezing parameterization are compared, including time‐independent (singular) number‐ and surface area‐based descriptions and a time‐dependent description following classical nucleation theory (CNT). The choice of freezing parameterization defines the size of the INP reservoir. The CNT‐based description yields an orders of magnitude larger INP reservoir than the singular parameterizations, which is the dominant factor for sustained ice crystal formation. The efficiency of the freezing process and cloud cooling are of secondary importance. A diagnostic treatment neglecting INP loss is only accurate when the INP reservoir size is large and INP depletion weak. Since a larger INP reservoir sustains ice crystal formation substantially longer, and ice water path scales with ice crystal concentrations for the conditions considered, resolving the source of differences in INP reservoir dynamics due to model implementation is a high priority for advancing climate model physics.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2023

A 1D Model for Nucleation of Ice From Aerosol Particles: An Application to a Mixed-Phase Arctic Stratus Cloud Layer

Mixed-phase clouds (MPCs) have been identified as significant contributors to uncertainties in climate projections, attributable to model representation of processes controlling the formation and loss of supercooled water droplets and ice particles from the atmosphere. Arctic MPCs are commonly widespread and long-lived, with sustained ice crystal formation processes that challenge current understanding. This study examines the ice-nucleating particle (INP) reservoir dynamics governing immersion-mode heterogeneous freezing in an observed case of Arctic MPCs using a simplified 1D aerosol-cloud model. The model setup includes prescribed dynamical forcings and thermodynamic profiles, and represents INPs as multicomponent and polydisperse particle size distributions. Diagnostic and prognostic approaches to immersion freezing parameterization are compared, including time-independent (singular) number- and surface area-based descriptions and a time-dependent description following classical nucleation theory (CNT). The choice of freezing parameterization defines the size of the INP reservoir. The CNT-based description yields an orders of magnitude larger INP reservoir than the singular parameterizations, which is the dominant factor for sustained ice crystal formation. The efficiency of the freezing process and cloud cooling are of secondary importance. A diagnostic treatment neglecting INP loss is only accurate when the INP reservoir size is large and INP depletion weak. Since a larger INP reservoir sustains ice crystal formation substantially longer, and ice water path scales with ice crystal concentrations for the conditions considered, resolving the source of differences in INP reservoir dynamics due to model implementation is a high priority for advancing climate model physics.

aerosol↗