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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
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.↗
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↗
NASA NTRS2013

Representation of Nucleation Mode Microphysics in a Global Aerosol Model with Sectional Microphysics

In models, nucleation mode (1 nm<Dp <10 nm) particle microphysics can be represented explicitly with aerosol microphysical processes or can be parameterized to obtain the growth and survival of nuclei to the model's lower size boundary. This study investigates how the representation of nucleation mode microphysics impacts aerosol number predictions in the TwO-Moment Aerosol Sectional (TOMAS) aerosol microphysics model running with the GISS GCM II-prime by varying its lowest diameter boundary: 1 nm, 3 nm, and 10 nm. The model with the 1 nm boundary simulates the nucleation mode particles with fully resolved microphysical processes, while the model with the 10 nm and 3 nm boundaries uses a nucleation mode dynamics parameterization to account for the growth of nucleated particles to 10 nm and 3 nm, respectively.We also investigate the impact of the time step for aerosol microphysical processes (a 10 min versus a 1 h time step) to aerosol number predictions in the TOMAS models with explicit dynamics for the nucleation mode particles (i.e., 3 nm and 1 nm boundary). The model with the explicit microphysics (i.e., 1 nm boundary) with the 10 min time step is used as a numerical benchmark simulation to estimate biases caused by varying the lower size cutoff and the time step. Different representations of the nucleation mode have a significant effect on the formation rate of particles larger than 10 nm from nucleated particles (J10) and the burdens and lifetimes of ultrafinemode (10 nm≤Dp ≤70 nm) particles but have less impact on the burdens and lifetimes of CCN-sized particles. The models using parameterized microphysics (i.e., 10 nm and 3 nm boundaries) result in higher J10 and shorter coagulation lifetimes of ultrafine-mode particles than the model with explicit dynamics (i.e., 1 nm boundary). The spatial distributions of CN10 (Dp ≥10 nm) and CCN(0.2 %) (i.e., CCN concentrations at 0.2%supersaturation) are moderately affected, especially CN10 predictions above ∼700 hPa where nucleation contributes most strongly to CN10 concentrations. The lowermost-layer CN10 is substantially improved with the 3 nm boundary (compared to 10 nm) in most areas. The overprediction in CN10 with the 3 nm and 10 nm boundaries can be explained by the overprediction of J10 or J3 with the parameterized microphysics, possibly due to the instantaneous growth rate assumption in the survival and growth parameterization. The errors in CN10 predictions are sensitive to the choice of the lower size boundary but not to the choice of the time step applied to the microphysical processes. The spatial distribution of CCN(0.2 %) with the 3 nm boundary is almost identical to that with the 1 nm boundary, but that with the 10 nm boundary can differ more than 10-40% in some areas. We found that the deviation in the 10 nm simulations is partly due to the longer time step (i.e., 1 h time step used in the 10 nm simulations compared to 10 min time step used in the benchmark simulations), but, even with the same time step, the 10 nm cutoff showed noticeably higher errors than the 3 nm cutoff. In conclusion, we generally recommend using a lower diameter boundary of 3 nm for studies focused on aerosol indirect effects but down to 1 nm boundary for studies focused on CN10 predictions or nucleation.

microphysics↗
NASA NTRS1971

An explicit form of the Mie phase matrix for multiple scattering calculations in the I, Q, U, and V representation

An explicit expression is obtained for the phase matrix in the I, Q, U, and V Stokes vector representation for a system containing a polydispersion of spherical particles. All of the symmetry relations derived by Hovenier using general arguments are established explicitly. Convenient algorithms are given for the computation of the phase matrix for a spherical polydispersion. Since this theory is so vitally important in radiative transfer, many researchers will need to compute these functions for realistic aerosols distributions. Therefore, results are presented for a haze L distribution so that other researchers will have a way of checking their programs which compute these quantities.

Kattawar, G. W.↗
NASA NTRS1991

A diagnostic stratospheric aerosol size distribution inferred from SAGE II measurements

An aerosol size distribution model for the stratosphere is inferred based on 5 years of Stratospheric Aerosol and Gas Experiment (SAGE) II measurements of multispectral aerosol and water vapor extinction. The SAGE II aerosol and water vapor extinction data strongly suggest that there is a critical particle radius below which there is a relatively weak dependence of particle number density with size and above which there are few, if any, particles. A segmented power law model, as a simple representation of this dependence, is used in theoretical calculations and intercomparisons with a variety of aerosol measurements including dustsondes, longwave lidar, and wire impactors and shows a consistently good agreement.

Thomason, Larry W.↗
NASA NTRS2018

Overview of Methods for Deriving the Radiative Transfer Theory from the Maxwell Equations. I: Approach Based on the Far-Field Foldy Equations

In this paper, we revisit, with further enhancements and clarifications, the self-consistent first-principles approach developed previously for deriving the vector radiative transfer theory for a discrete random medium with a sparse concentration of particles. We specifically consider the case of a plane-parallel particulate layer embedded in an otherwise homogeneous unbounded medium. The solution method is based on the far-field Foldy equations, an order-of-scattering expansion for the total field derived under the Twersky approximation, the computation of the coherent field by assuming that the positions of the particles are uncorrelated, and the ladder approximation for the coherency dyadic. The latter yields an integral equation for the diffuse specific coherency dyadic, defined through an angular spectrum repre- sentation for the coherency dyadic, which in turn, gives the vector radiative transfer equation for the diffuse specific intensity column vector. We analyze specifically the computation of the coherent field for inhomogeneous particulate media and multiple species of particles, the continuous extension of the far-field representation to the near field, the Foldy approximation, and the Foldy integral equation for the coherent field. Finally, we discuss the transition from the vector to the scalar radiative transfer equation.

Doicu, Adrian↗
NASA NTRS1983

Computer simulation of surface and film processes

Adequate computer methods, based on interactions between discrete particles, provide information leading to an atomic level understanding of various physical processes. The success of these simulation methods, however, is related to the accuracy of the potential energy function representing the interactions among the particles. The development of a potential energy function for crystalline SiO2 forms that can be employed in lengthy computer modelling procedures was investigated. In many of the simulation methods which deal with discrete particles, semiempirical two body potentials were employed to analyze energy and structure related properties of the system. Many body interactions are required for a proper representation of the total energy for many systems. Many body interactions for simulations based on discrete particles are discussed.

Tiller, W. A.↗
NASA NTRS1976

Quantum principles and free particles

The quantum principles that establish the energy levels and degeneracies needed to evaluate the partition functions are explored. The uncertainty principle is associated with the dual wave-particle nature of the model used to describe quantized gas particles. The Schroedinger wave equation is presented as a generalization of Maxwell's wave equation; the former applies to all particles while the Maxwell equation applies to the special case of photon particles. The size of the quantum cell in phase space and the representation of momentum as a space derivative operator follow from the uncertainty principle. A consequence of this is that steady-state problems that are space-time dependent for the classical model become only space dependent for the quantum model and are often easier to solve. The partition function is derived for quantized free particles and, at normal conditions, the result is the same as that given by the classical phase integral. The quantum corrections that occur at very low temperatures or high densities are derived. These corrections for the Einstein-Bose gas qualitatively describe the condensation effects that occur in liquid helium, but are unimportant for most practical purposes otherwise. However, the corrections for the Fermi-Dirac gas are important because they quantitatively describe the behavior of high-density conduction electron gases in metals and explain the zero point energy and low specific heat exhibited in this case.

Source record↗
NASA NTRS1973

Atmospheric secondary gamma radiation

A theoretical model in the form of a three-dimensional Monte Carlo calculation is developed that starts with known cosmic ray flux at the top of the atmosphere and follows particles as they interact with the atmosphere for a schematic representation of how atmospheric secondary gamma rays are produced. The model predicts atmospheric secondary gamma radiation as a function of altitude, energy, and zenith angle.

Thompson, D. J.↗
NASA NTRS1974

Lagrangian methods in plasma dynamics. II - Construction of Lagrangians for plasmas

Consideration of the construction of suitable Lagrangian functions for the dynamics of a cold plasma in such a way as to retain the relativistically covariant formalism. In one method, this is achieved by the introduction of a set of three variables which label the world lines of the particles. A second method results in a Clebsch-type representation. Sturrock's relativistic Lagrangian and Low's hot plasma Lagrangian are also briefly discussed in the context of the present work. The behavior of the canonical stress tensor is considered. The applicability of many of the general results in part I (Dougherty, 1970) is ensured by establishing the existence of the Lagrangian function.

Dougherty, J. P.↗
NASA NTRS2013

Evaluation of Aerosol Mixing State Classes in the GISS Modele-matrix Climate Model Using Single-particle Mass Spectrometry Measurements

Aerosol particles in the atmosphere are composed of multiple chemical species. The aerosol mixing state, which describes how chemical species are mixed at the single-particle level, provides critical information on microphysical characteristics that determine the interaction of aerosols with the climate system. The evaluation of mixing state has become the next challenge. This study uses aerosol time-of-flight mass spectrometry (ATOFMS) data and compares the results to those of the Goddard Institute for Space Studies modelE-MATRIX (Multiconfiguration Aerosol TRacker of mIXing state) model, a global climate model that includes a detailed aerosol microphysical scheme. We use data from field campaigns that examine a variety of air mass regimens (urban, rural, and maritime). At all locations, polluted areas in California (Riverside, La Jolla, and Long Beach), a remote location in the Sierra Nevada Mountains (Sugar Pine) and observations from Jeju (South Korea), the majority of aerosol species are internally mixed. Coarse aerosol particles, those above 1 micron, are typically aged, such as coated dust or reacted sea-salt particles. Particles below 1 micron contain large fractions of organic material, internally-mixed with sulfate and black carbon, and few external mixtures. We conclude that observations taken over multiple weeks characterize typical air mass types at a given location well; however, due to the instrumentation, we could not evaluate mass budgets. These results represent the first detailed comparison of single-particle mixing states in a global climate model with real-time single-particle mass spectrometry data, an important step in improving the representation of mixing state in global climate models.

climate models↗
NASA NTRS2019

Analysis of the Global Microwave Polarization Data of Clouds

Two NASA microwave radiometers, the satellite-borne GPM (Global Precipitation Measurement) Microwave Imager (GMI) and the aircraft-borne CoSMIR (the Conical Scanning Millimeter-wave Imaging Radiometer), measure vertically- and horizontally-polarized microwaves emitted by cloud particles and the Earth below, providing unique information on ice crystal properties in clouds. Their data reveal that non-spherical ice crystals are common and they fall in a preferred horizontally aligned orientation in convective and optically thick clouds especially near cloud top. A bin (particle-size-resolving) microphysical model with an ice crystal shape representation is developed to simulate the evolution of ice crystal properties (i.e., size, shape and orientation), where the radiation effect on microphysics (REM) is taken into account. Since REM represents the effect of all (e.g., both infrared and solar) radiation on ice crystal temperature, it relies upon the ice crystal proprieties that determine how an ice crystal receives radiation. Definitely, REM is different from the radiative effects that cause sensitivity at the microwave frequencies in the GPM and CoSMIR observations. Model results show that horizontally-oriented ice crystals grow faster than vertically-oriented (or spherical) ones due to REM. When both horizontally- and vertically-oriented ice crystals coexist in an air parcel, the model results show that the former grow by vapor deposition whereas the latter shrink by sublimation and disappear eventually. These modeling results are supported by the GMI data and the CoSMIR observations from MC3E (Midlatitude Continental Convective Clouds Experiment) and OLYMPEX (Olympic Mountains Experiment) on the prevalence of horizontally-oriented ice crystals. Moreover, the REM-induced precipitation explains the CloudSat observations of rare thin clouds in the tropical mid-troposphere as well as the common diamond dust in the high latitudes.

Xiping Zeng↗
NASA NTRS2022

An Evaluation of Biomass Burning Aerosol Mass, Extinction, and Size Distribution in Geos Using Observations From Camp2EX

Biomass burning aerosol impacts aspects of the atmosphere and Earth system through radiative forcing, serving as cloud condensation nuclei, and air quality. Despite its importance, the representation of biomass burning aerosol is not always accurate in numerical weather prediction and climate models or reanalysis products. Using observations collected as part of the Cloud, Aerosol and Monsoon Processes Philippines Experiment (CAMP2Ex) in August through October of 2019, aerosol concentration and optical properties are evaluated within the Goddard Earth Observing System (GEOS) and its underlying aerosol module, GOCART. In the operational configuration, GEOS assimilates aerosol optical depth observations at 550 nm to constrain aerosol fields. Particularly for biomass burning aerosol, without the assimilation of aerosol optical depth, aerosol extinction is underestimated compared to observations collected in the Philippines region during the CAMP2Ex campaign. The assimilation process adds excessive amounts of carbon to account for the underestimated extinction, resulting in positive biases in the mass of black and organic carbon, especially within the boundary layer, relative to in situ observations from the Langley Aerosol Research Group Experiment. Counteracting this, GEOS is deficient in sulphate and nitrate aerosol just above the boundary layer. Aside from aerosol mass, extinction within GEOS is a function of ambient relative humidity and an assumed particle size distribution. The relationship between dry and ambient extinction in GEOS reveals that hygroscopic growth is too aggressive within the model for biomass burning aerosol. An additional concern lies in the assumed particle size distribution for GEOS, which has a mode radius that is too small for organic carbon. Variability in the observed particle size distribution for biomass burning aerosol within a single flight also illuminates the fact that a single assumed particle size distribution is not sufficient and that for a proper representation, a more advanced aerosol module with GEOS may be necessary.

Aerosol Mass↗
NASA NTRSDate not supplied

Evaluation of GEOS Aerosols Using CAMP2Ex Observations

Biomass burning aerosol impacts aspects of the atmosphere and Earth system through direct and semi-direct effects, as well as influencing air quality. Despite its importance, the representation of biomass burning aerosol is not always accurate in numerical weather prediction and climate models or reanalysis products. Using observations collected as part of the Cloud, Aerosol and Monsoon Processes Philippines Experiment (CAMP2Ex) in August through October of 2019, aerosol concentration and optical properties are evaluated within the Goddard Earth Observing System (GEOS) and its underlying aerosol module, GOCART. In the operational configuration, GEOS assimilates aerosol optical depth observations at 550 nm from AERONET and MODIS to constrain aerosol fields. Particularly for biomass burning aerosol, without the assimilation of aerosol optical depth, aerosol extinction is underestimated compared to observations collected in the Philippines region during the CAMP2Ex campaign. The assimilation process adds excessive amounts of carbon to account for the underestimated extinction, resulting in positive biases in the mass of black and organic carbon, especially within the boundary layer, relative to in situ observations from the Langley Aerosol Research Group Experiment. Counteracting this, GEOS is deficient in sulfate and nitrate aerosol just above the boundary layer. Aerosol extinction within GEOS is a function of the mass of different aerosol species, the ambient relative humidity, the assumed spectral optical properties, and particle size distribution per species. The relationship between dry and ambient extinction in GEOS reveals that hygroscopic growth is too high within the model for biomass burning aerosol. An additional concern lies in the assumed particle size distribution for GEOS, which has a single mode radius that is too small for organic carbon. Variability in the observed particle size distribution for biomass burning aerosol within a single flight also illuminates the fact that a single assumed particle size distribution is not sufficient and that for a proper representation, a more advanced aerosol module within GEOS may be necessary.

Allison Collow↗
NASA NTRS2020

Models Transport Saharan Dust too Low in the Atmosphere: a Comparison of the MetUM and CAMS Forecasts with Observations

We investigate the dust forecasts from two operational global atmospheric models in comparison with in situ and remote sensing measurements obtained during the AERosol properties – Dust (AER-D) field campaign. Airborne elastic backscatter lidar measurements were performed on board the Facility for Airborne Atmospheric Measurements during August 2015 over the eastern Atlantic, and they permitted us to characterise the dust vertical distribution in detail, offering insights on transport from the Sahara. They were complemented with airborne in situ measurements of dust size distribution and optical properties, as well as datasets from the Cloud–Aerosol Transport System (CATS) spaceborne lidar and the Moderate Resolution Imaging Spectroradiometer (MODIS). We compare the airborne and spaceborne datasets to operational predictions obtained from the Met Office Unified Model (MetUM) and the Copernicus Atmosphere Monitoring Service (CAMS). The dust aerosol optical depth predictions from the models are generally in agreement with the observations but display a low bias. However, the predicted vertical distribution places the dust lower in the atmosphere than highlighted in our observations. This is particularly noticeable for the MetUM, which does not transport coarse dust high enough in the atmosphere or far enough away from the source.We also found that both model forecasts underpredict coarse-mode dust and at times overpredict fine-mode dust, but as they are fine-tuned to represent the observed optical depth, the fine mode is set to compensate for the underestimation of the coarse mode. As aerosol–cloud interactions are dependent on particle numbers rather than on the optical properties, this behaviour is likely to affect their correct representation. This leads us to propose an augmentation of the set of aerosol observations available on a global scale for constraining models, with a better focus on the vertical distribution and on the particle size distribution. Mineral dust is a major component of the climate system; therefore, it is important to work towards improving how models reproduce its properties and transport mechanisms.

Debbie OSullivan↗
NASA NTRS2018

Gas-Granular Simulation Framework for Spacecraft Landing Plume-Surface Interaction and Debris Transport Analysis

The Gas-Granular Flow Solver (GGFS) multi-phase flow computational framework has been developed to enable simulations of particle flows complex extra-terrestrial regolith materials. Particle flows of interest include the damage of unprepared landing sites from rocket plume impingement on Moon, Mars, and asteroids. The flow solver implements an Eulerian-Eulerian two-fluid model with fluid representation of the gas phase and granular phase to avoid the need to model billions of particle interactions. The granular phase is modeled as an Eulerian fluid with constituent physics closure models derived from first-principle Discrete Element Model (DEM) particle interaction simulations that capture the complex, non-linear granular particle interaction effects. Granular phase constituent models have been developed and integrated that address the complex, non-linear granular material mechanics complexities resulting from both: the irregular, jagged particle shapes and poly-disperse mixture effects encountered in extra-terrestrial regolith, with lunar regolith as the extreme. The GGFS capabilities are being integrated into a proven NASA plume-surface interaction and debris transport simulation framework featuring the Loci/CHEM CFD program and Debris Transport Analysis (DTA) post-processing tools for applications in robotic and human Moon and Mars lander development. Integration of the three simulation tool components. Loci/CHEM, GGFS, and DTA, into a coordinated simulation framework will enable time-accurate spacecraft landing simulations that account for the alteration of the landing surface through plume-induced cratering and the resulting redirection of plume impingement flow and debris transport. Initial implementation of this simulation framework and application examples will be presented.

Liever, Peter A.↗
NASA NTRS2012

Difference in Fractional Occurrences of Precipitation Categories in Terms of Cloud Properties

The relationship between the precipitation stage and water cloud properties was studied in the mid-latitudes using space-borne radar and radiometer data, with a focus on a comparison between continental (over China) and oceanic (over the northwest Pacific) clouds. More specifically, fractional occurrences of the precipitation categories were investigated in terms of two-dimensional representations of cloud parameters, such as combinations of optical depth - effective particle radius, and liquid water path - droplet number density. As a result, systematic changes in the transition regarding pairs of the cloud parameters were well captured for all precipitation categories. The transition pattern was similar for both land and oceanic clouds. Notably, the ranges of all values of cloud parameters (except cloud droplet number density) were considerably larger for oceanic clouds for the precipitating category.

drizzle↗
NASA NTRS2012

Statistical Orbit Determination using the Particle Filter for Incorporating Non-Gaussian Uncertainties

The tracking of space objects requires frequent and accurate monitoring for collision avoidance. As even collision events with very low probability are important, accurate prediction of collisions require the representation of the full probability density function (PDF) of the random orbit state. Through representing the full PDF of the orbit state for orbit maintenance and collision avoidance, we can take advantage of the statistical information present in the heavy tailed distributions, more accurately representing the orbit states with low probability. The classical methods of orbit determination (i.e. Kalman Filter and its derivatives) provide state estimates based on only the second moments of the state and measurement errors that are captured by assuming a Gaussian distribution. Although the measurement errors can be accurately assumed to have a Gaussian distribution, errors with a non-Gaussian distribution could arise during propagation between observations. Moreover, unmodeled dynamics in the orbit model could introduce non-Gaussian errors into the process noise. A Particle Filter (PF) is proposed as a nonlinear filtering technique that is capable of propagating and estimating a more complete representation of the state distribution as an accurate approximation of a full PDF. The PF uses Monte Carlo runs to generate particles that approximate the full PDF representation. The PF is applied in the estimation and propagation of a highly eccentric orbit and the results are compared to the Extended Kalman Filter and Splitting Gaussian Mixture algorithms to demonstrate its proficiency.

Mashiku, Alinda↗