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 127 records · Page 7

Atmospheric Dissolved Iron Depostiion to the Global Oceans: Effects of Oxalate-Promoted Fe Dissolution, Photochemical Redox Cycling, and Dust Mineralogy

Mineral dust deposition is suggested to be a significant atmospheric supply pathway of bioavailable iron (Fe) to Fe-depleted surface oceans. In this study, mineral dust and dissolved Fe (Fed) deposition rates are predicted for March 2009 to February 2010 using the 3-D chemical transport model GEOS-Chem implemented with a comprehensive dust-Fe dissolution scheme. The model simulates Fed production during the atmospheric transport of mineral dust taking into account inorganic and organic (oxalate)-promoted Fe dissolution processes, photochemical redox cycling between ferric (Fe(III)) and ferrous (Fe(II)) forms of Fe, dissolution of three different Fe-containing minerals (hematite, goethite, and aluminosilicates), and detailed mineralogy of windblown dust from the major desert regions. Our calculations suggest that during the yearlong simulation is approximately 0.26 Tg (1 Tg = 1012 g) of Fed was deposited to global oceanic regions. Compared to simulations only taking into account proton-promoted Fe dissolution, the addition of oxalate to the dust-Fe mobilization scheme increased total annual model-predicted Fed deposition to global oceanic regions by approximately 75%. The implementation of Fe(II)/Fe(III) photochemical redox cycling in the model allows for the distinction between different oxidation states of deposited Fed. Our calculations suggest that during the daytime, large fractions of Fed deposited to the global oceans is likely to be in Fe(II) form, while nocturnal fluxes of Fed are largely in Fe(III) form. Model simulations also show that atmospheric fluxes of Fed can be strongly influenced by the mineralogy of Fe-containing compounds. This study shows that Fed deposition to the oceans is controlled by total dust-Fe mass concentrations, mineralogy, the surface area of dust particles, atmospheric chemical composition, cloud processing, and meteorological parameters and exhibits complex and spatiotemporally variable patterns. Our study suggests that the explicit model representation of individual processes leading to Fed production within mineral dust are needed to improve the understanding of the atmospheric Fe cycle, and quantify the effect of dust-Fe on ocean biological productivity, carbon cycle, and climate.

Johnson, M. S.↗

Jovian plasma modeling for mission design

The purpose of this report is to address uncertainties in the plasma models at Jupiter responsible for surface charging and to update the jovian plasma models using the most recent data available. The updated plasma environment models were then used to evaluate two proposed Europa mission designs for spacecraft charging effects using the Nascap-2k code. The original Divine/Garrett jovian plasma model (or "DG1", T. N. Divine and H. B. Garrett, "Charged particle distributions in Jupiter's magnetosphere," J. Geophys. Res., vol. 88, pp. 6889-6903,1983) has not been updated in 30 years, and there are known errors in the model. As an example, the cold ion plasma temperatures between approx.5 and 10 Jupiter radii (Rj) were found by the experimenters who originally published the data to have been underestimated by approx.2 shortly after publication of the original DG1 model. As knowledge of the plasma environment is critical to any evaluation of the surface charging at Jupiter, the original DG1 model needed to be updated to correct for this and other changes in our interpretation of the data so that charging levels could beproperly estimated using the Nascap-2k charging code. As an additional task, the Nascap-2k spacecraft charging tool has been adapted to incorporate the so-called Kappa plasma distribution function--an important component of the plasma model necessary to compute the particle fluxes between approx.5 keV and 100 keV (at the outset of this study,Nascap-2k did not directly incorporate this common representation of the plasma thus limiting the accuracy of our charging estimates). The updating of the DG1 model and its integration into the Nascap-2k design tool means that charging concerns can now be more efficiently evaluated and mitigated. (We note that, given the subsequent decision by the Europa project to utilize solar arrays for its baseline design, surface charging effects have becomeeven more of an issue for its mission design). The modifications and results of those modifications to the DG1 model to produce the new DG2 model presented here and the steps taken to integrate the DG2 predictions into Nascap-2k are described in this report.

Voyager 2↗

Jovian Plasma Modeling for Mission Design

The purpose of this report is to address uncertainties in the plasma models at Jupiter responsible for surface charging and to update the jovian plasma models using the most recent data available. The updated plasma environment models were then used to evaluate two proposed Europa mission designs for spacecraft charging effects using the Nascap-2k code. The original Divine/Garrett jovian plasma model (or "DG1", T. N. Divine and H. B. Garrett, "Charged particle distributions in Jupiter's magnetosphere," J. Geophys. Res., vol. 88, pp. 6889-6903,1983) has not been updated in 30 years, and there are known errors in the model. As an example, the cold ion plasma temperatures between approx.5 and 10 Jupiter radii (Rj) were found by the experimenters who originally published the data to have been underestimated by approx.2 shortly after publication of the original DG1 model. As knowledge of the plasma environment is critical to any evaluation of the surface charging at Jupiter, the original DG1 model needed to be updated to correct for this and other changes in our interpretation of the data so that charging levels could beproperly estimated using the Nascap-2k charging code. As an additional task, the Nascap-2k spacecraft charging tool has been adapted to incorporate the so-called Kappa plasma distribution function--an important component of the plasma model necessary to compute the particle fluxes between approx.5 keV and 100 keV (at the outset of this study,Nascap-2k did not directly incorporate this common representation of the plasma thus limiting the accuracy of our charging estimates). The updating of the DG1 model and its integration into the Nascap-2k design tool means that charging concerns can now be more efficiently evaluated and mitigated. (We note that, given the subsequent decision by the Europa project to utilize solar arrays for its baseline design, surface charging effects have becomeeven more of an issue for its mission design). The modifications and results of those modifications to the DG1 model to produce the new DG2 model presented here and the steps taken to integrate the DG2 predictions into Nascap-2k are described in this report

Garrett, Henry B.↗

Multi-Scale Modeling of Low-Density Carbon-Phenolic Ablators

Protecting a spacecraft during atmospheric entry is one of highest risk factors that needs to be mitigated during design of a space exploration mission. At entry speeds from space, air turns into high-temperature plasma, and spacecraft Thermal Protection Systems (TPS) are needed to protect the vehicle payload. Modern successful material architectures of spacecraft shields use a porous carbon fiber substrate impregnated with phenolic as an ablator material. In the lecture, efforts to build a Predictive Material Modeling framework for porous ablators from micro-scale to macro-scale will be presented. Several numerical methods and techniques will be summarized that use voxelized images to compute geometrical properties of the porous substrate. These computed properties include porosity, specific surface area and tortuosity that are otherwise indirectly measured through experimental techniques. Direct simulation Monte Carlo (DSMC), a particle-based method for approximating the Boltzmann equation, is used to compute the permeability coefficient of the porous substrate based on its digitized representation. The method computes the flow within the microstructure, where the size of the pores may approach the mean-free-path of the flow. Finally, a high-fidelity model implemented in PATO (Porous-material Analysis Toolbox) is discussed, and some examples of ablative material response are presented, including for the first time 3D simulations of the full tiled heat shield for the Mars Science Laboratory (MSL) capsule.

Ablators↗

Theory of multicolor lattice gas - A cellular automaton Poisson solver

The present class of models for cellular automata involving a quiescent hydrodynamic lattice gas with multiple-valued passive labels termed 'colors', the lattice collisions change individual particle colors while preserving net color. The rigorous proofs of the multicolor lattice gases' essential features are rendered more tractable by an equivalent subparticle representation in which the color is represented by underlying two-state 'spins'. Schemes for the introduction of Dirichlet and Neumann boundary conditions are described, and two illustrative numerical test cases are used to verify the theory. The lattice gas model is equivalent to a Poisson equation solution.

Chen, H.↗

Remarks about massive and massless particles in supersymmetry

The internal space-time symmetry and simple supersymmetry of relativistic particles are briefly discussed in terms of the little group of the Poincare group. The little group generators in a finite-dimensional matrix representation of the N = 1 super-Poincare algebra are explicitly constructed. The supergeometry of a massive case continuously becomes that of a massless case in the infinite-momentum limit. The origin of the gage transformations associated with the massless supermultiplets becomes transparent in that limit.

Ketov, S. V.↗

Binary collision rates of relativistic thermal plasmas. II - Spectra

Spectra of importance for the analysis of relativistic thermal plasmas are numerically calculated assuming a thermal form for the particle distribution functions. Complete sets of optically thin thermal electron-proton, electron-electron, and electron-positron bremsstrahlung spectra are calculated throughout the transrelativistic regime of electron temperatures and compared with approximate expressions for the spectra in the nonrelativistic and extreme relativistic regimes of temperature. A method for calculating accurate secondary particle production spectra in proton-proton collisions from threshold to the highest energies is presented based on an isobaric model near threshold and scaling representations at high energies. The production spectra of charged and neutral pions resulting from proton-proton collisions in relativistic proton plasmas are calculated, and the resultant electron, positron, and gamma-ray spectra from the decay of secondary pions are presented.

Dermer, C. D.↗

Using XR for Improving Scientific Discovery With Numerical Weather Models

Earth science (ES) digital twins will help us understand the complex interactions and interrelationships that make up our Earth system and the impacts of earth science phenomena on it. Our work addresses two underdeveloped areas in current ES digital twin work: improving the understanding and interaction with ES model outputs by using Virtual and Mixed Reality (XR) tools and improving the non-intuitive mapping of continuous ES natural phenomena to gridded reference frames in current numerical models. Traditionally, scientists working on ES view and analyze the results of calculated or measured observables with static 1-dimensional (1D), 2D or 3D plots displayed on flat computer screens or paper. Using such limited mediums, it can be very difficult to identify, track and understand the evolution of key features due to poor viewing angles and the nature of flat computer screens. In addition, numerical models, such as the NASA Goddard Earth Observing System (GEOS) ES model, are almost exclusively formulated, visualized and analyzed in an Eulerian reference frame with fixed grid points in space and time. However, ES phenomena such as convective clouds, hurricanes and wildfire smoke plumes are visualized and analyzed in a Lagrangian reference frame: therefore it is often difficult and unnatural to understand these phenomena in relation to each other, visualized either in an Eulerian or Lagrangian context. In 3D visualizations, data generally takes one of three forms: gridded (e.g., voxelized) data, where space is divided into regions; point clouds, where data is represented as a set of points; and meshes, where objects are rendered as surfaces composed of small polygons (usually triangles). A gridded, Eulerian reference frame has been the default representation for the 2D visual analysis of atmospheric data in part because the numerical methods used to generate atmospheric model data in the first place use a gridded approach, with equations defining the relationships between the physical variables in each of a grid's cells across successive timesteps. In our work, we are particularly interested in data from GEOS. Another reason why gridded representations tend to be used for visualizing data from such models is because trajectories are difficult to interpret from representations on 2D surfaces, due to line-of-sight ambiguity. Instead of a fixed grid from GEOS, we embed a trajectory model to simulate particles' movement throughout a GEOS run. We then ingest these particle trajectories as animated point clouds with a NASA open source XR toolkit, the Mixed Reality Exploration Toolkit (MRET), and merge GEOS data with ES phenomena data onto one combined visualization that the user can intuitively interact with. Efficient rendering of arbitrarily large point clouds is an ongoing challenge being addressed by the computer science community, with the GPU-based optimizations and efficient GPU memory utilization a common theme of recent advances, especially for XR, where sustained high frame rate is mandatory to save the user from suffering due to simulation sickness. In this work, we describe and evaluate our progress in choosing and implementing appropriate methods for rendering arbitrarily large point clouds within MRET for XR. While tracking the XR headset enables the immersion of a user within a 3D scene of a data visualization, tracking of XR handheld controllers or user’s hands enables us to implement intuitive user interactions with the visualized datasets. Conventional tools require a user working with an ES visualization to conduct many interactions to commit their intended selections or manipulations with a visualized dataset; for example to specify a set of points in 3D space. Doing so in a 2D flat screen interface has traditionally required specifying a set of points in three distinct 2D coordinate systems (XY, XZ, and YZ), which is cumbersome. In other scientific domains, it has been shown that specifying or selecting a location or volume in XR using handheld controllers or tracked hands allows for greater speed and accuracy. We anticipate the same will hold true for atmospheric data, and we will share initial results of measuring the utility of such an interface. Notably, as the data being visualized is generated by GEOS as a prediction based on initial conditions, an intended application of our tool is to serve as part of an iterative feedback loop. Through XR, a scientist will review and manipulate a GEOS model run, modifying the conditions as needed to do subsequent runs of GEOS. Thereby, XR-based improvements to speed and accuracy of 3D tagging of points minimizes the effort required by both the scientist and the computer cluster conducting the necessary calculations.

Thomas Grubb↗

MATRIX-VBS (v1.0): Implementing an Evolving Organic Aerosol Volatility in an Aerosol Microphysics Model

The gas-particle partitioning and chemical aging of semi-volatile organic aerosol are presented in a newly developed box model scheme, where its effect on the growth, composition, and mixing state of particles is examined. The volatility-basis set (VBS) framework is implemented into the aerosol microphysical scheme MATRIX (Multiconfiguration Aerosol TRacker of mIXing state), which resolves mass and number aerosol concentrations and in multiple mixing-state classes. The new scheme, MATRIX-VBS, has the potential to significantly advance the representation of organic aerosols in Earth system models by improving upon the conventional representation as non-volatile particulate organic matter, often also with an assumed fixed size distribution. We present results from idealized cases representing Beijing, Mexico City, a Finnish forest, and a southeastern US forest, and investigate the evolution of mass concentrations and volatility distributions for organic species across the gas and particle phases, as well as assessing their mixing state among aerosol populations. Emitted semi-volatile primary organic aerosols evaporate almost completely in the intermediate-volatility range, while they remain in the particle phase in the low-volatility range. Their volatility distribution at any point in time depends on the applied emission factors, oxidation by OH radicals, and temperature. We also compare against parallel simulations with the original scheme, which represented only the particulate and non-volatile component of the organic aerosol, examining how differently the condensed-phase organic matter is distributed across the mixing states in the model. The results demonstrate the importance of representing organic aerosol as a semi-volatile aerosol, and explicitly calculating the partitioning of organic species between the gas and particulate phases.

volatility-basis set↗

Spacecraft Attitude Determination Methods

This document is presentation in viewgraph form, which outlines the methods of determining spacecraft attitude. The presentation reviews several parameterizations relating to spacecraft attitude, such as Euler's Theorem, Rodriques parameters, and Euler-Rodriques parameters or Quaternion. Onboard attitude determination is the norm, using either single frame or filtering methods. The presentation reviews several mathematical representations of attitude. The mechanisms for determining attitude on board the Hubble Space Telescope, the Tropical Rainfall and Measuring Mission and the Solar Anomalous and Magnetospheric Particle Explorer are reviewed. Wahba's problem, Procrustes Problem, and some solutions are also summarized.

Markley, F. Landis↗

Analysis of the charging of the SCATHA (P78-2) satellite

The charging of a large object in polar Earth orbit was investigated in order to obtain a preliminary indication of the response of the shuttle orbiter to such an environment. Two NASCAP (NASA Charging Analyzer Program) models of SCATHA (Satellite Charging at High Altitudes) were used in simulations of charging events. The properties of the satellite's constituent materials were compiled and representations of the experimentally observed plasma spectra were constructed. Actual charging events, as well as those using test environments, were simulated. Numerical models for the simulation of particle emitters and detectors were used to analyze the operation of these devices onboard SCATHA. The effect of highly charged surface regions on the charging conductivity within a photosheath was used to interpret results from the onboard electric field experiment. Shadowing calculations were carried out for the satellite and a table of effective illuminated areas was compiled.

Stannard, P. R.↗

Swelling of Transported Smoke from Savanna Fires over the Southeast Atlantic Ocean

We use the recently released Cloud Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) Version 4.1 (V4) lidar data to study the smoke plumes transported from Southern African biomass burning areas. Significant improvements in the CALIPSO V4 Level 1 calibration and V4 Level 2 algorithms lead to a better representation of their optical properties, with the aerosol subtype improvements being particularly relevant to smoke over this area. For the first time, we show evidence of smoke particles increasing in size, evidenced in their particulate color ratios, as they are transported over the South Atlantic Ocean from the source regions over Southern Africa. We hypothesize that this is due to hygroscopic swelling of the smoke particles and is reflected in the higher relative humidity in the middle troposphere for profiles with smoke. This finding may have implications for radiative forcing estimates over this area and is also relevant to the ORACLES field mission.

Kar, J.↗

Swelling of Transported Smoke from Savanna Fires over the Southeast Atlantic Ocean

We use the recently released Cloud Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) Version 4.1 (V4) lidar data to study the smoke plumes transported from Southern African biomass burning areas. Significant improvements in the CALIPSO V4 Level 1 calibration and V4 Level 2 algorithms lead to a better representation of their optical properties, with the aerosol subtype improvements being particularly relevant to smoke over this area. For the first time, we show evidence of smoke particles increasing in size, as demonstrated by their particulate color ratios, as they are transported over the South Atlantic Ocean from the source regions over Southern Africa. We hypothesize that this is due to hygroscopic swelling of the smoke particles and is reflected in the higher relative humidity in the middle troposphere for profiles with smoke. This finding may have implications for radiative forcing estimates over this area and is also relevant to the ORACLES field mission.

J. Kar↗

Apparent dust size discrepancy in aerosol reanalysis in north African dust after long-range transport

North African dust reaches the southeastern United States every summer. Size-resolved dust mass measurements taken in Miami, Florida, indicate that more than one-half of the surface dust mass concentrations reside in particles with geometric diameters less than 2.1 µm, while vertical profiles of micropulse lidar depolarization ratios show dust reaching above 4 km during pronounced events. These observations are compared to the representation of dust in the Modern-Era Retrospective analysis for Research and Applications, version 2 (MERRA-2) aerosol reanalysis and closely related Goddard Earth Observing System model version 5 (GEOS-5) Forward Processing (FP) aerosol product, both of which assimilate satellite-derived aerosol optical depths using a similar protocol and inputs. These capture the day-to-day variability in aerosol optical depth well, in a comparison to an independent sun-photometer-derived aerosol optical depth dataset. Most of the modeled dust mass resides in diameters between 2 and 6 µm, in contrast to the measurements. Model-specified mass extinction efficiencies equate light extinction with approximately 3 times as much aerosol mass, in this size range, compared to the measured dust sizes. GEOS-5 FP surface-layer sea salt mass concentrations greatly exceed observed values, despite realistic winds and relative humidities. In combination, these observations help explain why, despite realistic total aerosol optical depths, (1) free-tropospheric model volume extinction coefficients are lower than those retrieved from the micro-pulse lidar, suggesting too-low model dust loadings in the free troposphere, and (2) model dust mass concentrations near the surface can be higher than those measured. The modeled vertical distribution of dust, when captured, is reasonable. Large, aspherical particles exceeding the modeled dust sizes are also occasionally present, but dust particles with diameters exceeding 10 µm contribute little to the measured total dust mass concentrations after such long-range transport. Remaining uncertainties warrant a further integrated assessment to confirm this study's interpretations.

MEERA-2 aerosols reanalysis↗

Polarized Radiative Transfer of a Cirrus Cloud Consisting of Randomly Oriented Hexagonal Ice Crystals: The 3 x 3 Approximation for Non-Spherical Particles

The reflection and transmission of polarized light for a cirrus cloud consisting of randomly oriented hexagonal columns were calculated by two very different vector radiative transfer models. The forward peak of the phase function for the ensemble-averaged ice crystals has a value of order 6 x 10(exp 3) so a truncation procedure was used to help produce numerically efficient yet accurate results. One of these models, the Vectorized Line-by-Line Equivalent model (VLBLE), is based on the doubling- adding principle, while the other is based on a vector discrete ordinates method (VDISORT). A comparison shows that the two models provide very close although not entirely identical results, which can be explained by differences in treatment of single scattering and the representation of the scattering phase matrix. The relative differences in the reflected I and Q Stokes parameters are within 0.5 for I and within 1.5 for Q for all viewing angles. In 1971 Hansen showed that for scattering by spherical particles the 3 x 3 approximation is sufficient to produce accurate results for the reflected radiance I and the degree of polarization (DOP), and he conjectured that these results would hold also for non-spherical particles. Simulations were conducted to test Hansen's conjecture for the cirrus cloud particles considered in this study. It was found that the 3 x 3 approximation also gives accurate results for the transmitted light, and for Q and U in addition to I and DOP. For these non-spherical ice particles the 3 x 3 approximation leads to an absolute error 2 x 10(exp -6) for the reflected and transmitted I, Q and U Stokes parameters. Hence, it appears to be an excellent approximation, which significantly reduces the computational complexity and burden required for multiple scattering calculations.

Polarization↗

Three-dimensional separated flows past ground vehicles

The separated flow past a thick body is calculated using a theoretical model based on a discrete wake-vortex representation. Rehbach's unsteady inviscid scheme (1977) was used in the three-dimensional analysis. The obstacle is represented by singularities, and the wake is modeled by vortex particles emitted at the separation line. Calculated results for an axisymmetric body are compared with flow visualization results.

Chometon, F.↗

Grid-free simulation of diffusion using random walk methods

The simulation of the diffusion of a continuum field by the random walk (RW) displacement of a set of particles is considered. Elements of the gradients of the diffusive concentration are transported by computational particles. It is demonstrated that, by the use of concentration gradients in the RW process, statistical errors are reduced and each realization of the numerical solution is a representation of the exact solution. The algorithm is grid-free, and the computational elements move to follow the gradients; hence, the algorithm is self-adaptive, and uniform resolution is achieved for all times.

Ghoniem, A. F.↗