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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

A High-Performance Discrete-Element Framework for Simulating Flow and Jamming of Moisture Bearing Biomass Feedstocks

We developed and verified a high-performance open-source discrete element method (DEM) solver with simultaneously-supported feedstock-specific interaction models, including bonded-sphere, liquid bridge, cohesion, and non-linear contact models. Our solver uses parallel data structures on hybrid central and graphics processing unit (CPU/GPU) architectures, with favorable strong scaling performance observed for large problem sizes comprised of (100 M particles), and 4X single-node GPU speedup. The particles for corn stover feedstock were conceptualized and calibrated based on experimental measurements and results. Sensitivity analyses demonstrate that the mass flow rate from a wedge hopper is governed primarily by moisture content, friction coefficient, and cohesion energy density. The model is used to reproduce experimentally observed hopper jamming results, highlighting that the experimental no-flow trends can only be achieved by using non-spherical particles, liquid bridge and cohesion models, highlighting the importance of using concurrent feedstock specialized models for the effective representation of biomass material handling problems.

bioenergy↗
DOE OSTI2025

Effect of Ice Number Concentration on the Evolution of Boundary Layer Clouds During Arctic Marine Cold‐Air Outbreaks

Abstract Marine cold‐air outbreaks (MCAOs) are crucial for Arctic Ocean heat loss, featuring convective cloud rolls that transition into convection cells downstream. Understanding factors controlling this transformation is the key for improving MCAO cloud representation in climate models. This study employs large‐eddy simulations to investigate how cloud ice number concentrations () affect cloud evolution using a case from the Cold‐Air Outbreaks in the Marine Boundary Layer Experiment (COMBLE) campaign. The simulations, performed in a Lagrangian framework following an air mass trajectory, are driven by ERA5 reanalysis data. Initially, all simulations produce similar cloud patterns, but higher leads to earlier breakup of cloud rolls. Between 4 and 10 hr, surface precipitation rates are similar across simulations, but precipitation initiates earlier, and the cloud‐base precipitation rates are higher when is higher. The stronger precipitation evaporation leads to increased stability of the boundary layer and reduced intensity of vertical mixing between the surface and cloud layer. An increased sink of cloud layer moisture via precipitation and decreased source through diminished vertical transport result in earlier cloud breakup in higher conditions. Simulations with different sea surface temperatures (SST) indicate that this cloud breakup mechanism remains valid for MCAOs of different strengths, although the cloud organization is more sensitive to SST changes in low environments. This work highlights the importance of accurate representations of ice processes in simulating MCAO clouds and suggests the need for observational constraints of ice nucleating particles and over the mixed‐phase cloud regimes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

mesoflow [SWR-22-56]

Mesoflow is a continuum scale simulation tool developed specifically for modeling transport and chemistry at the mesoscale. Our solver utilizes Cartesian block-structured adaptive mesh refinement to resolve complex surface morphologies (of catalysts/biomass particles among others) directly obtained from X-ray tomography data. An immersed boundary based formulation enables rapid representation of complex geometries prevalent in most mesoporous interfaces. The solver is developed on top of open-source performance portable library, AMReX, providing parallel execution capabilities on current and upcoming high-performance-computing (HPC) architectures. Our flexible software framework enables integration of complex chemical mechanisms at heterogenous interfaces and time-split algorithms for circumventing highly disparate reaction and flow time-scales. Our current studies indicate a ten-fold performance gain by using graphics-processing-units (GPU) compared to a single processor for representative problem sizes (2 million cell mesh).

Sitaraman, Hariswaran↗
NASA NTRS1988

Correlation studies on surface particle detection methods

The accurate determination of dust levels on optical surfaces is necessary to assess sensor system performance. A comparison study was made on several particle measurement methods including those based on direct imaging and light scattering. The effectiveness of removing the particles from the surface prior to determining particle size distributions was also assessed. These studies revealed that some methods, especially those requiring particle removal before analysis, are subject to large systematic errors affecting particle size distributions. Thus, an understanding of the particle measurement methods employed is necessary before any surface cleanliness or obstruction value assignments are accepted as true representations of an optical surface contamination condition.

Peterson, Ronald V.↗
NASA NTRS2013

Influence of Computational Drop Representation in LES of a Droplet-Laden Mixing Layer

Multiphase turbulent flows are encountered in many practical applications including turbine engines or natural phenomena involving particle dispersion. Numerical computations of multiphase turbulent flows are important because they provide a cheaper alternative to performing experiments during an engine design process or because they can provide predictions of pollutant dispersion, etc. Two-phase flows contain millions and sometimes billions of particles. For flows with volumetrically dilute particle loading, the most accurate method of numerically simulating the flow is based on direct numerical simulation (DNS) of the governing equations in which all scales of the flow including the small scales that are responsible for the overwhelming amount of dissipation are resolved. DNS, however, requires high computational cost and cannot be used in engineering design applications where iterations among several design conditions are necessary. Because of high computational cost, numerical simulations of such flows cannot track all these drops. The objective of this work is to quantify the influence of the number of computational drops and grid spacing on the accuracy of predicted flow statistics, and to possibly identify the minimum number, or, if not possible, the optimal number of computational drops that provide minimal error in flow prediction. For this purpose, several Large Eddy Simulation (LES) of a mixing layer with evaporating drops have been performed by using coarse, medium, and fine grid spacings and computational drops, rather than physical drops. To define computational drops, an integer NR is introduced that represents the ratio of the number of existing physical drops to the desired number of computational drops; for example, if NR=8, this means that a computational drop represents 8 physical drops in the flow field. The desired number of computational drops is determined by the available computational resources; the larger NR is, the less computationally intensive is the simulation. A set of first order and second order flow statistics, and of drop statistics are extracted from LES predictions and are compared to results obtained by filtering a DNS database. First order statistics such as Favre averaged stream-wise velocity, Favre averaged vapor mass fraction, and the drop stream-wise velocity, are predicted accurately independent of the number of computational drops and grid spacing. Second order flow statistics depend both on the number of computational drops and on grid spacing. The scalar variance and turbulent vapor flux are predicted accurately by the fine mesh LES only when NR is less than 32, and by the coarse mesh LES reasonably accurately for all NR values. This is attributed to the fact that when the grid spacing is coarsened, the number of drops in a computational cell must not be significantly lower than that in the DNS.

Bellan, Josette↗
DOE OSTI2020

Wave function methods for canonical ensemble thermal averages in correlated many-fermion systems

We present a wave function representation for the canonical ensemble thermal density matrix by projecting the thermofield double state against the desired number of particles. Furthermore, the resulting canonical thermal state obeys an imaginary-time evolution equation. Starting with the mean-field approximation, where the canonical thermal state becomes an antisymmetrized geminal power (AGP) wave function, we explore two different schemes to add correlation: by number-projecting a correlated grand-canonical thermal state and by adding correlation to the number-projected mean-field state. As benchmark examples, we use number-projected configuration interaction and an AGP-based perturbation theory to study the hydrogen molecule in a minimal basis and the six-site Hubbard model.

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

Properties and Controlling Processes of Aerosol and Cloud Condensation Nuclei in Marine Boundary Layer Over Eastern North Atlantic (Final Report)

Remote marine low cloud systems have large spatial and temporal coverages. Because of their relatively low optical thickness and background aerosol concentrations, marine low clouds are particularly susceptible to perturbations in aerosols associated with anthropogenic emissions. Recent studies show that a large fraction of the global aerosol indirect forcing can be attributed to changes in marine low clouds, despite their relatively long distance from most anthropogenic sources. Currently, there remain large uncertainties in the magnitude of the global aerosol indirect radiative forcing. One major contribution to this large uncertainty derives from poor understanding of the aerosols under natural conditions, the perturbation by anthropogenic emissions, and the processes that drive them. In this project, we took advantage of the comprehensive, multi-dimensional characterization of meteorological parameters, trace gases, aerosol, cloud, and drizzle fields during the Aerosol and Cloud Experiments in Eastern North Atlantic (ACE-ENA) and the long-term measurements at the Atmospheric Radiation Measurement (ARM) Eastern North Atlantic (ENA) site on Graciosa Island. By synergistically combining in situ measurements, remote sensing retrievals, and model simulations, we quantified aerosol properties in the ENA, their vertical profiles, and variations with season, synoptic state, and airmass. We investigated key aerosol processes, including entrainment of free troposphere aerosol into marine boundary layer, new particle formation, growth of Aitken mode particles, and the influences of these processes on the cloud condensation nuclei (CCN) population in marine boundary layer. The representation of key aerosol processes in the Energy Exascale Earth System Model (E3SM) was evaluated using the observations. This project has led to 26 publications, which are listed at the end of this report.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1972

Interplanetary gas dynamics.

Observations of the interplanetary medium are considered together with the structure of the solar atmosphere and the kinetic properties of the interplanetary plasma. Methods of gas kinetics and continuum flows are examined, taking into account microscopic and macroscopic representations, the kinetic flow equations, the formal solution of the Vlasov equation, and the continuum flow equations. The collective particle behavior of that interplanetary gas is discussed along with the hydrodynamic coronal expansion, a wave-pump problem, the free expansion phenomenon, and a generalized free expansion problem.

Liu, V. C.↗
NASA NTRS1984

A Lagrangian method for a two layer shallow water model

One of the biggest drawbacks of Eulerian methods is the presence of the nonlinear convective terms which leads to inaccurate representations of advection and fronts. An alternative and conceptually simpler approach is the Lagrangian approach in which the fluid particles themselves are tracked and equations are derived based on local spatial interaction. In this formulation the nonlinear convective terms do not appear. The main drawback with the Lagrangian formulation, however, is that the local spatial interactions are time dependent. That is, at each time level, one must know the current neighbors of a given particle to accurately compute the forces acting on the particle. The Voroni Mesh which is continuously deforming and always linking nearest neighbors is needed. The effectiveness of a simple Lagrangian method for simulating the onset of a Kelvin-Helmholtz instability in a 2 layer fluid is demonstrated.

Augenbaum, J. M.↗
NASA NTRS2018

Impact of the Ozone Monitoring Instrument Row Anomaly on the Long-Term Record of Aerosol Products

Since about three years after the launch the Ozone Monitoring Instrument (OMI) on the EOS-Aura satellite, the sensor's viewing capability has been affected by what is believed to be an internal obstruction that has reduced OMI's spatial coverage. It currently affects about half of the instrument's 60 viewing positions. In this work we carry out an analysis to assess the effect of the reduced spatial coverage on the monthly average values of retrieved aerosol optical depth (AOD), single scattering albedo (SSA) and the UV Aerosol Index (UVAI) using the 2005-2007 three-year period prior to the onset of the row anomaly. Regional monthly average values calculated using viewing positions 1 through 30 were compared to similarly obtained values using positions 31 through 60, with the expectation of finding close agreement between the two calculations. As expected, mean monthly values of AOD and SSA obtained with these two scattering-angle dependent subsets of OMI observations agreed over regions where carbonaceous or sulphate aerosol particles are the predominant aerosol type. However, over arid regions, where desert dust is the main aerosol type, significant differences between the two sets of calculated regional mean values of AOD were observed. As it turned out, the difference in retrieved desert dust AOD between the scattering-angle dependent observation subsets was due to the incorrect representation of desert dust scattering phase function. A sensitivity analysis using radiative transfer calculations demonstrated that the source of the observed AOD bias was the spherical shape assumption of desert dust particles. A similar analysis in terms of UVAI yielded large differences in the monthly mean values for the two sets of calculations over cloudy regions. On the contrary, in arid regions with minimum cloud presence, the resulting UVAI monthly average values for the two sets of observations were in very close agreement. The discrepancy under cloudy conditions was found to be caused by the parameterization of clouds as opaque Lambertian reflectors. When properly accounting for cloud scattering effects using Mie theory, the observed UVAI angular bias was significantly reduced. The analysis discussed here has uncovered important algorithmic deficiencies associated with the model representation of the angular dependence of scattering effects of desert dust aerosols and cloud droplets. The resulting improvements in the handling of desert dust and cloud scattering have been incorporated in an improved version of the OMAERUV algorithm.

Atmospheric Chemistry↗
DOE OSTI2024

Robust Tensor Hypercontraction of the Particle–Particle Ladder Term in Equation-of-Motion Coupled Cluster Theory

One method of representing a high-rank tensor as a (hyper-)product of lower-rank tensors is the tensor hypercontraction (THC) method of Hohenstein et al. This strategy has been found to be useful for reducing the polynomial scaling of coupled-cluster methods by representation of a four-dimensional tensor of electron-repulsion integrals in terms of five two-dimensional matrices. Pierce et al. have already shown that the application of a robust form of THC to the particle–particle ladder (PPL) term reduces the cost of this term in couple-cluster singles and doubles (CCSD) from O(N 6 ) to O(N 5 ) with negligible errors in energy with respect to the density-fitted variant. In this work, we have implemented the least-squares variant of THC (LS-THC) which does not require a nonlinear tensor factorization, including the robust form (R-LS-THC), for the calculation of the excitation and electron attachment energies using equation-of-motion coupled cluster methods EOMEE-CCSD and EOMEA-CCSD, respectively. We have benchmarked the effect of the R-LS-THC-PPL approximation on excitation energies using the comprehensive QUEST database and the accuracy of electron attachment energies using the NAB22 database. Here, we find that errors on the order of 1 meV are achievable with a reduction in total calculation time of approximately 5x.

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

The viscous magnetopause

Magnetopause representation by mixing region of plasma streams with different velocities and magnetic fields, assuming hydromagnetic viscosity caused by wave-particle interactions

Cassen, P.↗
NASA NTRS2020

Global Airborne Sampling Reveals a Previously Unobserved Dimethyl Sulfide Oxidation Mechanism in the Marine Atmosphere

Dimethyl sulfide (DMS), emitted from the oceans, is the most abundant biological source of sulfur to the marine atmosphere. Atmospheric DMS is oxidized to condensable products that form secondary aerosols that affect Earth’s radiative balance by scattering solar radiation and serving as cloud condensation nuclei. We report the atmospheric discovery of a previously unquantified DMS oxidation product, hydroperoxymethyl thioformate (HPMTF, HOOCH2SCHO), identified through global-scale airborne observations that demonstrate it to be a major reservoir of marine sulfur. Observationally constrained model results show that more than 30% of oceanic DMS emitted to the atmosphere forms HPMTF. Coincident particle measurements suggest a strong link between HPMTF concentration and new particle formation and growth. Analyses of these observations show that HPMTF chemistry must be included in atmospheric models to improve representation of key linkages between the biogeochemistry of the ocean, marine aerosol formation and growth, and their combined effects on climate.

DMS gas-phase oxidation product↗
NASA NTRS2024

A Multi‐Probe Automated Classification of Ice Crystal Habits During the IMPACTS Campaign

Although all ice crystals are unique, many can be grouped together by shape or habit, with members of a habit class sharing similar representations of properties such as fall velocity and growth rate. A decision tree algorithm designed to be adaptable to any particle imaging probe, thus enabling the creation of habit size distributions over a size range larger than that of any probe on its own, is used to classify ice crystals imaged by three airborne cloud probes in mid-latitude winter cyclones during the Investigation of Microphysics and Precipitation for Atlantic Coast-Threatening Snowstorms (IMPACTS) field campaign. Crystals are sorted into seven habit classes based on their morphological properties: sphere, column/needle, plate, graupel, dendrite, aggregate, and irregular. Although adaptability was its primary goal, the algorithm was found to be moderately skillful for identifying idealized habit images. Quantitative tests of the algorithm’s adaptability displayed mixed results, as Two-Dimensional Stereo Probe (2DS) classifications showed moderate correlation with Particle Habit Imaging and Polar Scattering Probe (PHIPS) classifications, but only weak correlation with High Volume Precipitation Spectrometer (HVPS) classifications. The algorithm was applied to random sets of images from each probe in a case study of a mesoscale snow band sampled on 7 February 2020. In the case study, qualitative analysis of particle images revealed general agreement on classifications among the probes, supporting the algorithm’s applicability to multiple cloud probes. Most classifications appeared correct upon manual inspection, suggesting that in practical use, the algorithm is reasonably able to classify non-idealized images.

Julian Schima↗
DOE OSTI2021

Role of quadrupole deformation and continuum effects in the “island of inversion” nuclei 28,29,31 F

The peculiar properties of nuclei in the so-called “island of inversion” around Z = 10 and N = 20 are the focus of current nuclear physics research. Recent studies showed that 28 F has a negative-parity ground state and thus lies within the southern shore of the island of inversion, and 29 F presents a halo structure in its ground state, but it is unclear which effects, such as deformation, shell evolution due to tensor forces, or couplings to the continuum, lead to this situation. Here, we investigate the role of quadrupole deformation and continuum effects on the single-particle structure of 28,29,31 F from a relativistic mean-field approach and show how both phenomena can lead to a negative-parity ground state in 28 F and halo structures in 29,31 F. We solve the Dirac equation in the complex-momentum (Berggren) representation for a potential with quadrupole deformation at the first order obtained from relativistic mean-field calculations using the NL3 interaction and calculate the continuum level densities using the Green’s function method. We extract single-particle energies and widths from the continuum level densities to construct the Nilsson diagrams of 28,29,31 F in the continuum and analyze the evolution of both the widths and occupation probabilities of relevant Nilsson orbitals in 28 F and find that some amount of prolate deformation must be present. In addition, we calculate the density distributions for bound Nilsson orbitals near the Fermi surface in 29,31 F and reveal that, for a quadrupole deformation 0.3 ≤ β 2 ≤ 0.45 (prolate), characteristic halo tails appear at large distances. Using the relativistic mean-field approach in the complex-momentum representation with the Green’s function method, we demonstrate that in neutron-rich fluorine isotopes, while in the spherical case the pf shells are already inverted and close to the neutron emission threshold, a small amount of quadrupole deformation can dramatically reduce the gap between positive- and negative-parity states and increase the role of continuum states, ultimately leading to the negative parity in the ground state of 28 F and the halo structures in 29,31 F.

20 ≤ A ≤ 38↗
DOE OSTI2024

Cluster-dynamics-based parameterization for sulfuric acid–dimethylamine nucleation: comparison and selection through box and three-dimensional modeling

Clustering of gaseous sulfuric acid (SA) enhanced by dimethylamine (DMA) is a major mechanism for new particle formation (NPF) in polluted atmospheres. However, uncertainty remains regarding the SA–DMA nucleation parameterization that reasonably represents cluster dynamics and is applicable across various atmospheric conditions. This uncertainty hinders accurate three-dimensional (3-D) modeling of NPF and the subsequent assessment of its environmental and climatic impacts. Here we extensively compare different cluster-dynamics-based parameterizations for SA–DMA nucleation and identify the most reliable one through a combination of box model simulations, 3-D modeling, and in situ observations. Results show that the parameterization derived from Atmospheric Cluster Dynamic Code (ACDC) simulations, incorporating the latest theoretical insights (DLPNO-CCSD(T)/aug-cc-pVTZ//ωB97X-D/6-311++G(3df,3pd) level of theory) and adequate representation of cluster dynamics, exhibits dependable performance in 3-D NPF simulation for both winter and summer conditions in Beijing and shows promise for application in diverse atmospheric conditions. Another ACDC-derived parameterization, replacing the level of theory with RI-CC2/aug-cc-pV(T+d)Z//M06-2X/6–311++G(3df,3pd), also performs well in NPF modeling at relatively low temperatures around 280 K but exhibits limitations at higher temperatures due to inappropriate representation of SA–DMA cluster thermodynamics. Additionally, a previously reported parameterization incorporating simplifications is applicable for simulating NPF in polluted atmospheres but tends to overestimate particle formation rates under conditions of elevated temperature (>∼300 K) and low-condensation sink (<∼3×10 -3 s −1 ). Our findings highlight the applicability of the new ACDC-derived parameterization, which couples the latest SA–DMA nucleation theory and holistic cluster dynamics, in 3-D NPF modeling. The ACDC-derived parameterization framework provides a valuable reference for developing parameterizations for other nucleation systems.

Shen, Jiewen↗
DOE OSTI2025

Measurements and Analyses to Enable Science for the Eastern Pacific Cloud Aerosol Precipitation Experiment (EPCAPE)

Coastal cities provide the opportunity to characterize the substantial effects of manmade particles on marine cloud properties and processes. La Jolla lies to the north of San Diego, California, but it is often about a day directly downwind of the major pollution sources located in the ports of Los Angeles and Long Beach. The large dynamic range of aerosol particle concentrations combined with the multi-hour to multi-day persistence of stratocumulus cloud layers makes the site ideal for investigating the seasonal changes in cloud and aerosol properties as well as the quantitative relationships between cloud and aerosol properties. The Eastern Pacific Cloud Aerosol Precipitation Experiment (EPCAPE) characterized the extent, radiative properties, aerosol interactions, and precipitation characteristics of stratocumulus clouds in the Eastern Pacific across all four seasons at two coastal sites in La Jolla. This project was designed to enhance and expand the scientific uses of the ARM AMF1 measurements from EPCAPE. The goal was to ensure meaningful observations were collected that would enable science. The project included the following objectives: (1) Reviewing the ARM AMF1 measurements [ARM, 2021a; b] and circulating a summary of the measurements each week, (2) Comparing AMF1 measurements to those provided by collaborators (including filter measurements at the pier), (3) Collecting and analyzing filter samples from Scripps Pier by Fourier Transform Infrared spectroscopy (FTIR) and X-ray Fluorescence (XRF), (4) Assisting in operations of instruments provided by Guest PIs when possible, and (5) Providing an initial compilation of EPCAPE aerosol and cloud seasonal differences. The expected outcomes of these objectives were enhanced proposals and publications using EPCAPE measurements by helping to identify instrument issues, expanded data access and awareness by distributing weekly plots and related summaries, improved source-related attribution of aerosols with elemental tracers, additional observations provided by Guest PIs, and accelerated ACI studies enabled by the compiled seasonal summaries of aerosol and cloud properties. Three examples of the science enabled by this project are findings that (i) aerosol and cloud aqueous production contributes more than half of sulfate particle mass concentration, (ii) upwind sources make chemical composition very similar at nearby sites despite local differences in meteorology, and (iii) most of the large mass concentration of semi-volatile organic components is co-emitted and co-evaporated with nitrate. Together these findings illustrate how ARM extended field campaigns in coastal regions can be used to constrain ACI processes with direct observations. By using the unique ARM suite of cloud radiative products in addition to the measured aerosol properties at a coastal location, we were able to address the more specific question of which aerosol particles cause how much of the effects on clouds. Identifying this signature in coastal areas provides an opportunity to test the representation of aerosol sources by global models in a range of clean and urban-influenced conditions.

Russell, Lynn [Univ. of California, San Diego, CA ↗