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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
NASA NTRS1975

The inversion of aureole measurements to derive aerosol size distributions

An iterative method to invert size distributions from simulated scattered radiance measurements at small angles from the sun has been investigated. The inferred size distributions were represented by piecewise linear and cubic spline functions. Various relevant characteristics were investigated and it was found that: (1) the inverted size distribution was insensitive to the number of knots in the piecewise linear spline; (2) within the range of sensitivity, the choice of initial guess had little effect on the inverted size distribution; (3) five per cent random noise in the simulated radiances appreciably deteriorated the result but variations are still tolerable when compared with other methods for determining size distributions; (4) the inverted distribution was insensitive to the index of refraction used in the kernel for particle radii greater than 1 micron; (5) the choice of wavelength between 0.40 and 0.70 microns has a negligible effect on the inverted distribution; (6) a range of tropospheric aerosol size distributions gives acceptable inverted results; and (7) the cubic spline representation can give reasonable inverted distributions, but may become unstable.

Twitty, J. T.↗
NASA NTRS1978

Remote sensing of environmental particulate pollutants - Optical methods for determinations of size distribution and complex refractive index

A unifying approach, based on a generalization of Pearson's differential equation of statistical theory, is proposed for both the representation of particulate size distribution and the interpretation of radiometric measurements in terms of this parameter. A single-parameter gamma-type distribution is introduced, and it is shown that inversion can only provide the dimensionless parameter, r/ab (where r = particle radius, a = effective radius, b = effective variance), at least when the distribution vanishes at both ends. The basic inversion problem in reconstructing the particle size distribution is analyzed, and the existing methods are reviewed (with emphasis on their capabilities) and classified. A two-step strategy is proposed for simultaneously determining the complex refractive index and reconstructing the size distribution of atmospheric particulates.

Fymat, A. L.↗
NASA NTRS1982

On the stability of Schwarzschild's triaxial galaxy model

The numerical model for a triaxial galaxy constructed by Schwartzschild (1979) has been tested by means of fully self-consistent three-dimensional n-body numerical experiments that use a representation of this model to start the integration. The integration, study models, initial load, and correcting potential of the experiment are detailed. Preliminary checks included density values, potential values, and integration in potential of an ideal model. The results showed some particle loss and no significant shape change. The model is rugged and robust, and no growing disturbances with growth rate in excess of 0.5 per crossing time could be detected. Runs were done at 50 percent excess halo and with tube orbits reversed.

Smith, B. F.↗
NASA NTRS1996

Ionospheric convection during different phases of magnetospheric substorms

Observations of ionospheric convection flows at a range of local times during the various phases of the substorm cycle are reported on with the aim of investigating the convection behavior during a range of times and phases. The ionospheric flow observations are from the EISCAT and DMSP satellites. The substorm phases are identified from energetic particle measurements from geosynchronous satellites. The growth phase convection indicates an initial expansion of the polar cap. There is an unexplained poleward motion of the flow reversal boundary (FRB). It is concluded that this motion does not necessarily provide a true representation of the balance between reconnection at the dayside and in the tail. The expansion phase flows do not show any evidence for tail reconnection until late in the phase. The convection during the recovery phase is indicative of tail reconnection as there is evidence that there is only a lobe cell driving convection on the dayside.

Lester, M.↗
NASA NTRS2022

Importance of Different Parameterization Changes for the Updated Dust Cycle Modeling in the Community Atmosphere Model (Version 6.1)

The Community Atmosphere Model (CAM6.1), the atmospheric component of the Community Earth System Model (CESM; version 2.1), simulates the life cycle (emission, transport, and deposition) of mineral dust and its interactions with physio-chemical components to quantify the impacts of dust on climate and the Earth system. The accuracy of such quantifications relies on how well dust-related processes are represented in the model. Here we update the parameterizations for the dust module, including those on the dust emission scheme, the aerosol dry deposition scheme, the size distribution of transported dust, and the treatment of dust particle shape. Multiple simulations were undertaken to evaluate the model performance against diverse observations, and to understand how each update alters the modeled dust cycle and the simulated dust direct radiative effect. The model–observation comparisons suggest that substantially improved model representations of the dust cycle are achieved primarily through the new more physically-based dust emission scheme. In comparison, the other modifications induced small changes to the modeled dust cycle and model–observation comparisons, except the size distribution of dust in the coarse mode, which can be even more influential than that of replacing the dust emission scheme. We highlight which changes introduced here are important for which regions, shedding light on further dust model developments required for more accurately estimating interactions between dust and climate.

mineral dust parameterizations↗
DOE OSTI2022

The MP{sub N} method: a new angular discretization method based on piecewise polynomial interfaces fluxes

In transport calculations, it is well known how S{sub N} method is extremely inefficient in problems where the particle physics is dominated by streaming. The ray-effect eventually produced by the insufficient angular discretization, appears to be extremely persistent with respect to the refinement of the angular quadrature. The MP{sub N} method, that relies on continuous angular representation, offers a robust remedy to such an issue. MP{sub N} is based on the decomposition of the unit sphere into solid angles and on a piecewise continuous definition of interface fluxes, which are expanded in polynomials in each solid angle. This allows propagating more than one angular degree of freedom simultaneously while maintaining unaltered the block-diagonal pattern of the displacement plus removal operator. The method is therefore well suited for the flux resolution by means of a conventional sweep algorithm. Furthermore, unlike the S{sub N} method, MP{sub N} does not rely on discrete directions and, thus, on an angular quadrature formula, but rather constructs a set of linear equations solving for the angular moments of the flux for all discrete solid angles within the sweep. MP{sub N} shows an error convergence rate higher than S{sub N} at the expense of an increased size of the coefficient matrices, so of the computational cost. Although MP{sub N} is not free from ray-effect, the latter is effectively mitigated and less persistent with respect to the increase of the angular refinement order. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Semi-Abelian gauge theories, non-invertible symmetries, and string tensions beyond N -ality

We study a 3d lattice gauge theory with gauge group U(1) N– 1 ⋉ S N , which is obtained by gauging the S N global symmetry of a pure U(1) N– 1 gauge theory, and we call it the semi-Abelian gauge theory. We compute mass gaps and string tensions for both theories using the monopole-gas description. We find that the effective potential receives equal contributions at leading order from monopoles associated with the entire SU( N ) root system. Even though the center symmetry of the semi-Abelian gauge theory is given by ℤ N , we observe that the string tensions do not obey the N -ality rule and carry more detailed information on the representations of the gauge group. We find that this refinement is due to the presence of non-invertible topological lines as a remnant of U(1) N– 1 one-form symmetry in the original Abelian lattice theory. Upon adding charged particles corresponding to W -bosons, such non-invertible symmetries are explicitly broken so that the N -ality rule should emerge in the deep infrared regime.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Topological heavy fermions in magnetic field

The recently introduced topological heavy fermion model (THFM) provides a means for interpreting the low-energy electronic degrees of freedom of the magic angle twisted bilayer graphene as hybridization amidst highly dispersing topological conduction and weakly dispersing localized heavy fermions. In order to understand the Landau quantization of the ensuing electronic spectrum, a generalization of THFM to include the magnetic field B is desired, but currently missing. Here we provide a systematic derivation of the THFM in B and solve the resulting model to obtain the interacting Hofstadter spectra for single particle charged excitations. While naive minimal substitution within THFM fails to correctly account for the total number of magnetic subbands within the narrow band i.e., its total Chern number, our method—based on projecting the light and heavy fermions onto the irreducible representations of the magnetic translation group— reproduces the correct total Chern number. Analytical results presented here offer an intuitive understanding of the nature of the (strongly interacting) Hofstadter bands.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

MontePy: a Python library for reading, editing, and writing MCNP input files.

The Monte Carlo N-Particle (MCNP) radiation transport code is a highly capable and accurate code with a long legacy. MCNP uses the Monte Carlo simulation process to simulate the path of particles (e.g., neutrons, photons, charged particles, etc.), and their interaction with materials. It is widely used in nuclear engineering, high-energy physics, and other fields. Its origins in the mid-twentieth century predate many modern software conventions. MCNP users provide an input file to MCNP, which it then uses to create an internal representation of the simulation problem. These input files originally had to be stored as punchcard decks, and the user manual still uses the terminology of cards and decks, despite moving beyond punchcards. MCNP predates nearly all modern human readable markup or data serialization languages, such as the extensible Markup Language (XML), the Standard Generalized Markup Language (SGML), YAML (YAML Ain’t Markup Language), and Javascript Object Notation (JSON). Due to this, MCNP uses an entirely custom defined syntax language for its input, making off-the-shelf libraries for XML, YAML, and JSON impossible to use for scripting various operations on MCNP input files (Kulesza et al., 2022).

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Scheme for quantitative description of longitudinal drifts in the Fermilab Linac and their correction

The Fermilab Linac accepts the 0.75 MeV H- ions from the front end and accelerates them to 400 MeV for injection into the Booster. Day-to-day drifts of the longitudinal trajectory in the Linac, reconstructed from phase readings of Beam Position Monitors, are at the level of several degrees. They are believed to cause additional losses both in the Linac and Booster, and are addressed by empirically adjusting the phases of Linac cavities. This work explores the option of expressing these drifts in terms of phase shifts in two cavities at the low-energy part of the Linac. Such a description allows for a simplified visual representation of the drifts, suggests a clear algorithm for their compensation, and provides a tool for estimating efficiency of such compensation.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Quantitative Description and Correction of Longitudinal Drifts in the Fermilab Linac

The Fermilab Linac accepts the 0.75 MeV H- ions from the front end and accelerates them to 400 MeV for injection into the Booster. Day-to-day drifts of the longitudinal trajectory in the Linac, reconstructed from phase readings of Beam Position Monitors, are at the level of several degrees. They are believed to cause additional losses both in the Linac and Booster, and are addressed by empirically adjusting the phases of Linac cavities. This work explores the option of expressing these drifts in terms of phase shifts in two cavities at the low-energy part of the Linac. Such description allows for a simplified visual representation of the drifts, suggest a clear algorithm for their compensation, and provides a tool for estimating efficiency of such compensation.

43 PARTICLE ACCELERATORS↗
NASA NTRS1981

Changes in cosmic ray propagation induced by corotating interaction regions

By simulating the trajectories for scatter free and diffusive propagation of relativistic cosmic rays in a model of the heliospheric magnetic fields containing a representation of corotating interaction regions (CIR's), it is found that, due to the large gradients associated with these compression regions, the motion is strongly affected and differs substantially from the predictions of current modulation theory. For positive (outward) northern hemisphere polarity, particles do not stream purely from high latitudes but can come from almost any latitude in the outer heliosphere; for negative polarity, many particles come along the current sheet (as predicted) but a second, equally important, population exists comprising particles that do not start on the current sheet but are brought to low latitudes by their interaction with CIR's. Thus, it is concluded that CIR's (and other large scale structures) cannot be ignored in analyses of cosmic ray modulation.

Gall, R.↗
NASA NTRS1994

The stability of perfect elliptic disks. 1: The maximum streaming case

Self-consistent distribution functions are constructed for two-dimensional perfect elliptic disks (for which the potential is exactly integrable) in the limit of maximum streaming; these are tested for stability by N-body integration. To obtain a discrete representation for each model, simulated annealing is used to choose a set of orbits which sample the distribution function and reproduce the required density profile while carrying the greatest possible amount of angular momentum. A quiet start technique is developed to place particles on these orbits uniformly in action-angle space, making the initial conditions as smooth as possible. The roundest models exhibit spiral instabilities similar to those of cold axisymmetric disks; the most elongated models show bending instabilities like those seen in prolate systems. Between these extremes, there is a range of axial ratios 0.25 approximately less than b/a approximately less than 0.6 within which these models appear to be stable. All the methods developed in this investigation can easily be extended to integrable potentials in three dimensions.

Levine, Stephen E.↗
NASA NTRSDate not supplied

Experimental and Computational Analysis of a Rotating Detonation Combustor

Rotating detonation combustor (RDC) research has progressed along parallel experimental and computational paths with limited opportunities for validation. Validation is a crucial step to ensure RDC simulations provide an accurate representation of the physical phenomena, and thus, can be used with confidence for design optimization in different applications. In this study, data from a quasi-two-dimensional computational fluid dynamic simulation in an annular RDC are compared with experimental measurements including particle image velocimetry at 100 kHz. Results show good agreement between measured and simulated dynamic pressures, Chapman-Jouguet speeds, and static pressure distributions. Velocity measurements at the annulus exit show good temporal and quantitative agreement with simulations including spikes from the passing shock.

Detonation↗
DOE OSTI2022

Design Optimization for EUROPA Critical Experiment

The Experiment for Unresolved Region Of Plutonium Actinides (EUROPA) is an integral critical experiment currently being designed to target the intermediate energy region of plutonium. Intermediate energies, those between 0.625 eV - 100 keV, contain the end of resolved and beginning of the unresolved resonance region, making experiments in this energy range prudent for validating the representation of cross sections in this transition region. Below is the preliminary design of the experiment to be performed at the National Criticality Experiments Research Center (NCERC). Experiment optimization required an exhaustive down selection of several moderating, absorbing, and reflecting materials using Particle Swarm Optimization (PSO) in order to achieve maximum sensitivity to the intermediate energy region.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2024

Quantitative description and correction of longitudinal drifts in the Fermilab linac

The Fermi National Accelerator Laboratory (Fermilab) Linac accepts 750 keV H- ions from the front end and accelerates them to 400 MeV for injection into the Booster rapid cycling synchrotron. Day-to-day drifts in the beam longitudinal trajectory during regular operation are of the order of several degrees. They are believed to cause additional losses in both the Linac and the Booster and are addressed by empirically adjusting cavity phases of front end and Linac RF cavities. This work explores a scheme for expressing these drifts in terms of phase shifts in the low-energy part of the Linac. Such a description allows for a simplified visual representation of the drifts, suggests a clear algorithm for their compensation, and provides a tool for estimating efficiency of such compensation.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Object detection with deep learning for rare event search in the GADGET II TPC

In the pursuit of identifying rare two-particle events within the GADGET II Time Projection Chamber (TPC), this paper presents a comprehensive approach for leveraging Convolutional Neural Networks (CNNs) and various data processing methods. To address the inherent complexities of 3D TPC track reconstructions, the data is expressed in 2D projections and 1D quantities. This approach capitalizes on the diverse data modalities of the TPC, allowing for the efficient representation of the distinct features of the 3D events, with no loss in topology uniqueness. Additionally, it leverages the computational efficiency of 2D CNNs and benefits from the extensive availability of pre-trained models. Given the scarcity of real training data for the rare events of interest, simulated events are used to train the models to detect real events. To account for potential distribution shifts when predominantly depending on simulations, significant perturbations are embedded within the simulations. This produces a broad parameter space that works to account for potential physics parameter and detector response variations and uncertainties. These parameter-varied simulations are used to train sensitive 2D CNN object detectors. When combined with 1D histogram peak detection algorithms, this multi-modal detection framework is highly adept at identifying rare, two-particle events in data taken during experiment 21072 at the Facility for Rare Isotope Beams (FRIB), demonstrating a 100% recall for events of interest. Here, we present the methods and outcomes of our investigation and discuss the potential future applications of these techniques.

Convolutional neural network↗