Search NASA⌕ Search

SEARCH · Search NASA

Results for “Energy Code Field Study”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Reassessing energy deposition for the ITER 5 MA vertical displacement event with an improved DINA model

The beryllium (Be) main chamber wall interaction during a 5 MA/1.8 T upward, unmitigated VDE scenario, first analysed in [J. Coburn et al., Phys. Scr. T171 (2020) 014076] for ITER, has been re-evaluated using the latest energy deposition analysis software. Updates to the DINA disruption model are summarized, including an improved numerical convergence for the 0D power balance, limitations on the safety factor within the plasma core, and the choice to maintain a constant plasma +halo poloidal cross-section. Such updates result in a broad halo region and higher radiated power fractions compared to previous models. The new scenario lasts for ~75 ms and deposits ~29 MJ of energy, with the radial distribution of parallel heat flux q‖(r)resembling an exponential falloff with an effective λ_q=75-198 mm. A maximum halo width w_h of 0.52 m at the outboard midplane is observed. SMITER field line tracing and energy deposition simulations calculate a q_{⊥,max} of ~83 MW/m^2 on the upper first wall panels (FWP). Heat transfer calculations with the MEMOS-U code show that the FWP surface temperature reaches ~1000 K, well below the Be melt threshold. Variations of this 5 MA scenario with Be im-purity densities from 0 to 3∙10^19 m^-3 also remain below the melt threshold despite differences in energy deposition and duration. These results are in contrast to the early study which predicted melt damage to the first wall [J. Coburn et al., Phys. Scr. T171 (2020) 014076], and emphasize the importance of accurate models for the halo width w_h and the heat flux distribution q‖(r)within that halo width. The 2020 halo model in DINA has been compared with halo current experiments on COMPASS, JET, and Alcator C-Mod, and the preliminary results build confidence in the broad halo width predictions. Results for the 5 MA VDE are compared with those for a 15 MA equivalent, generated using the new DINA model. At the higher current, significant melting of the upper FWP is to be expected.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2004

Swift Burst Alert Telescope (BAT) Instrument Response

The Burst Alert Telescope (BAT), a large coded aperture instrument with a wide field-of-view (FOV), provides the gamma-ray burst triggers and locations for the Swift Gamma-Ray Burst Explorer. In addition to providing this imaging information, BAT will perform a 15 keV - 150 keV all-sky hard x-ray survey based on the serendipitous pointings resulting from the study of gamma-ray bursts and will also monitor the sky for transient hard x-ray sources. For BAT to provide spectral and photometric information for the gamma-ray bursts, the transient sources and the all-sky survey, the BAT instrument response must be determined to an increasingly greater accuracy. In this talk, we describe the BAT instrument response as determined to an accuracy suitable for gamma-ray burst studies. We will also discuss the public data analysis tools developed to calculate the BAT response to sources at different energies and locations in the FOV. The level of accuracy required for the BAT instrument response used for the hard x-ray survey is significantly higher because this response must be used in the iterative clean algorithm for finding fainter sources. Because the bright sources add a lot of coding noise to the BAT sky image, fainter sources can be seen only after the counts due to the bright sources are removed. The better we know the BAT response, the lower the noise in the cleaned spectrum and thus the more sensitive the survey. Since the BAT detector plane consists of 32768 individual, 4 mm square CZT gamma-ray detectors, the most accurate BAT response would include 32768 individual detector response functions to separate mask modulation effects from differences in detector efficiencies! We describe OUT continuing work to improve the accuracy of the BAT instrument response and will present the current results of Monte Carlo simulations as well as BAT ground calibration data.

Parsons, A.↗
NASA NTRS2015

Solid Rocket Motor Combustion Instability Modeling in COMSOL Multiphysics

Combustion instability modeling of Solid Rocket Motors (SRM) remains a topic of active research. Many rockets display violent fluctuations in pressure, velocity, and temperature originating from the complex interactions between the combustion process, acoustics, and steady-state gas dynamics. Recent advances in defining the energy transport of disturbances within steady flow-fields have been applied by combustion stability modelers to improve the analysis framework. Employing this more accurate global energy balance requires a higher fidelity model of the SRM flow-field and acoustic mode shapes. The current industry standard analysis tool utilizes a one dimensional analysis of the time dependent fluid dynamics along with a quasi-three dimensional propellant grain regression model to determine the SRM ballistics. The code then couples with another application that calculates the eigenvalues of the one dimensional homogenous wave equation. The mean flow parameters and acoustic normal modes are coupled to evaluate the stability theory developed and popularized by Culick. The assumption of a linear, non-dissipative wave in a quiescent fluid remains valid while acoustic amplitudes are small and local gas velocities stay below Mach 0.2. The current study employs the COMSOL Multiphysics finite element framework to model the steady flow-field parameters and acoustic normal modes of a generic SRM. This work builds upon previous efforts to verify the use of the acoustic velocity potential equation (AVPE) laid out by Campos. The acoustic velocity potential (psi) describing the acoustic wave motion in the presence of an inhomogeneous steady high-speed flow is defined by, del squared psi - (lambda/c) squared psi - M x [M x del((del)(psi))] - 2((lambda)(M)/c + M x del(M) x (del)(psi) - 2(lambda)(psi)[M x del(1/c)] = 0. with M as the Mach vector, c as the speed of sound, and λ as the complex eigenvalue. The study requires one way coupling of the CFD High Mach Number Flow (HMNF) and mathematics module. The HMNF module evaluates the gas flow inside of a SRM using St. Robert's law to model the solid propellant burn rate, slip boundary conditions, and the supersonic outflow condition. Results from the HMNF model are verified by comparing the pertinent ballistics parameters with the industry standard code outputs (i.e. pressure drop, axial velocity, exit velocity). These results are then used by the coefficient form of the mathematics module to determine the complex eigenvalues of the AVPE. The mathematics model is truncated at the nozzle sonic line, where a zero flux boundary condition is self-satisfying. The remaining boundaries are modeled with a zero flux boundary condition, assuming zero acoustic absorption on all surfaces. The one way coupled analysis is perform four times utilizing geometries determined through traditional SRM modeling procedures. The results of the steady-state CFD and AVPE analyses are used to calculate the linear acoustic growth rate as is defined by Flandro and Jacob. In order to verify the process implemented within COMSOL we first employ the Culick theory and compare the results with the industry standard. After the process is verified, the Flandro/Jacob energy balance theory is employed and results displayed.

Fischbach, S. R.↗
NASA NTRS1996

Studies of Westward Electrojets and Field-Aligned Currents in the Magnetotail During Substorms: Implications for Magnetic Field Models

This section outlines those tasks undertaken in the final year that contribute integrally to the overarching project goals. Fast, during the final year, it is important to note that the project benefited greatly with the addition of a Boston University graduate student, Ms. Karen Hirsch. Jointly, we made substantial progress on the development of and improvements to magnetotail magnetic field and plasma models. The ultimate aim of this specific task was to assess critically the utility of such models for mapping low-altitude phenomena into the magnetotail (and vice-versa). The bulk of this effort centered around the finite-width- magnetotail convection model developed by and described by Spence and Kivelson (J. Geophys. Res., 98, 15,487, 1993). This analytic, theoretical model specifies the bulk plasma characteristics of the magnetotail plasma sheet (number density, temperature, pressure) across the full width of the tail from the inner edge of the plasma sheet to lunar distances. Model outputs are specified by boundary conditions of the source particle populations as well as the magnetic and electric field configuration. During the reporting period, we modified this code such that it can be interfaced with the auroral particle precipitation model developed by Dr. Terry Onsager. Together, our models provide a simple analytic specification of the equatorial distribution of fields and plasma along with their low-altitude consequences. Specifically, we have built a simple, yet powerful tool which allows us to indirectly 'map' auroral precipitation signatures (VDIS, inverted-V's, etc.) measured by polar orbiting spacecraft in the ionosphere, to the magnetospheric equatorial plane. The combined models allow us to associate latitudinal gradients measured in the ion energy fluxes at low-altitudes with the large-scale pressure gradients in the equatorial plane. Given this global, quasi-static association, we can then make fairly strong statements regarding the location of discrete features in the context of the global picture. We reported on our initial study at national and international meetings and published the results of our predictions of the low-altitude signatures of the plasma sheet. In addition, the PI was invited to contribute a publication to the so-called 'Great Debate in Space Physics' series that is a feature of EOS. The topic was on the nature of magnetospheric substorms. Specific questions of the when and where a substorm occurs and the connection between the auroral and magnetospheric components were discussed in that paper. This paper therefore was derived exclusively from the research supported by this grant. Attachment: Empirical modeling of the quite time nightside magnetosphere.' 'CRRES observations of particle flux dropout event.' The what, where, when, and why of magnetospheric substorm triggers'. and 'Low altitude signature of the plasma sheet: model prediction of local time dependence'.

Spence, Harlan E.↗
DOE OSTI2022

Efficient reconstruction and validation of heterogeneous microstructures for energy applications

The digital reconstruction of microstructures is necessary for simulations in fields ranging from geology to electrochemistry, but the state-of-the-art digital reconstruction techniques often compromise between resolution and field of view. It is challenging to retain detailed microstructure information in large-scale reconstructions. Here, this study investigates different aspects of the Yeong-Torquato algorithm based on correlation functions to make it more efficient. We achieve this goal by reducing the computational complexity of the chord-length distribution function and the two-point correlation function, applying the random sphere-packing method as the initial condition, and restricting potential voxel swaps to interfaces. In addition, a novel superposition parallel scheme is introduced to aid in searching for potential voxel swaps. The algorithm proposed is validated by comparing the pore-size distributions of reconstructed 3D custom battery electrodes from a sample dataset obtained from transmission X-ray microscopy. From a sample image with 200 x 200 pixels, the code can reconstruct a 300 x 300 x 300 structure in under 22 h and reconstruct a 400 x 400 x 400 structure in 43 h with eight cores.

42 ENGINEERING↗
DOE OSTI2025

Solenoidal Fields from Permanent Magnets for a Polarized Electron Preinjector

This technical note describes the design of a compact focusing solenoid constructed from permanent bar magnets, intended for the low-energy section of a spin-polarized electron injector. Typically, quadrupole triplets provide appropriate focusing for the low-energy spin rotator section; however, their large footprint precludes the insertion of sufficient beam diagnostics. The primary motivation of this study is to provide beam focusing while preserving the vertical spin orientation of the electron beam. We propose a solenoidal field composed of multiple permanent magnets arranged to naturally yield a zero longitudinal magnetic field integral ($\int B_zdz = 0$). This design serves as a compact, passive, and non-coupling alternative to traditional quadrupole triplets. We present the magnet arrangement, field calculation methodology, spin tracking simulations using space charge codes, and calculations demonstrating that the required field strength is achievable.

43 PARTICLE ACCELERATORS↗
NASA NTRS1985

Numerical techniques for large cosmological N-body simulations

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

Efstathiou, G.↗
DOE OSTI2020

Continued Validation Studies using the MOOSE Framework for Plasma Simulation with Electromagnetics

Resources and tools for the modeling and simulation of low-temperature plasma (LTP) discharges are increasingly vital to progress in the field in order to properly characterize and study complex source designs and plasma chemistries beyond the scope of traditional diagnostics. Open-source software provides powerful platforms for this work and can enable community-driven LTP R&D. Within the Multiphysics Object-Oriented Simulation Environment (MOOSE) open-source framework [1], capabilities have been demonstrated in the areas of plasma fluids (Zapdos [2]), plasma chemistry (CRANE [3]), and general electromagnetic wave theory (Electromagnetic Library for Kinetics & fluids [ELK] [4]). ELK has since been coupled to Zapdos/CRANE to enable fully coupled electromagnetic plasma simulations. This talk will detail the continued validation efforts and discuss Zapdos-ELK-CRANE code coupling with various-low temperature plasma sources. The impact of fully coupled electromagnetics on process parameters (e.g., temperature and electron/ion energy) versus an electrostatic description will also be discussed. [1] Permann et al., SoftwareX, 11 (2020) 100430. [2] Lindsay et al., J. Phys. D: Appl. Phys. 49 (2016) 235204. [3] Keniley et al., Talk, GEC Session DT3.00002 (2019). [4] Icenhour et al., Talk, GEC Session RR1.00006 (2019). Acknowledgements: INL Graduate Fellowship Program, NSF SI2 Grant 1740300, U.S. Dept. of Energy Office of Science Graduate Student Research Program

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Sensitivities of time-dependent temperature profile predictions for NSTX with the multi-mode model

The Multi-Mode Model (MMM) for turbulent transport was applied to a large set of well-analyzed discharges from the National Spherical Torus Experiment (NSTX) in order to evaluate its sensitivities to a wide range of plasma conditions. MMM calculations were performed for hundreds of milliseconds in each discharge by performing time-dependent predictive simulations with the 1.5D tokamak integrated modeling code TRANSP. A closely related study (Lestz et al 2025 Plasma Phys. Control. Fusion 67 105029) concluded that MMM predicted electron and ion temperature profiles that were in reasonable agreement with NSTX observations, generally outperforming a different reduced transport model, TGLF. This finding motivates the more thorough investigation of the characteristics of the MMM predictions conducted in this work. The simulations with MMM have electron energy transport dominated by electron temperature gradient modes in the examined discharges with relatively low plasma β (ratio of kinetic plasma pressure to magnetic field pressure) and high collisionality, transitioning to a mixture of different modes for higher β and lower collisionality. The thermal ion diffusivity predicted by MMM is much smaller than the neoclassical contribution, in line with previous experimental analysis of NSTX. Nonetheless, the electron and ion temperature profiles are coupled via collisional energy exchange and thus sensitive to which transport channels are predicted. The time-dependent simulations with MMM are robust to the simulation start time, converging to remarkably similar temperature profiles later during the discharge. MMM typically overpredicts confinement relative to NSTX observations, leading to the prediction of overly steep temperature profiles. Plasmas with spatially broader temperature profiles, higher plasma β, and longer energy confinement times tend to be predicted by MMM with better agreement with the experiment. As a result, these findings provide useful context for understanding the regime-dependent tendencies of MMM in anticipation of self-consistent, time-dependent predictive simulations of NSTX-U discharges with these same modeling tools.

MMM↗
DOE OSTI2023

Investigation of Nanoparticle Degradation in Hydrogen Fuel Cell Systems through Automated Electron Microscopy

Proton exchange membrane fuel cells (PEMFC) are promising devices for the deployment of hydrogen-powered heavy-duty vehicles, providing a higher efficiency for similar driving range and fueling time than the existing ones. However, PEMFCs still encounter durability challenges mainly due to catalyst degradation in the cathode. Mitigating these performance losses requires a better understanding of the degradation mechanisms under heavy-duty accelerated stress tests (ASTs) [1]. Scanning transmission electron microscopy (STEM) combined with energy dispersive X-ray spectroscopy (EDS) are key tools for the analysis of Pt and PtCo nanoparticle size, spatial distribution and composition [2]. Electron tomography is also used to determine the rate and type of degradation of catalyst nanoparticles as a function of their position on the carbon support. In this work, automated data acquisition software, paired with a custom Python code, have been used to study the effect of different accelerated stress tests (ASTs) on nanoparticle coarsening [2]. Figure 1 shows high-angle annular dark-field (HAADF)-STEM images and EDS maps comparing the cathodes of membrane electrode assemblies (MEAs) following an electrocatalyst AST performed under H2/N2 with that of the heavy-duty AST performed under H2/air. We will discuss how AST conditions affect considerably the spatial distribution of the nanoparticles across the electrode between the membrane and microporous layer. Although the median particle size increased more in the MEA aged under the heavy-duty AST, as determined using a high-throughput image analysis, the quantitative EDS measurements demonstrate that the electrocatalyst AST resulted in more Pt and Co dissolution from the cathode, which is another important indicator of electrocatalyst degradation. We will further present the impact of the relative humidity (% RH) on the degradation mechanisms demonstrated using the same approach. Electron tomography has been used to distinguish the Pt nanoparticles residing on the carbon support surface (exterior) from those within the pore structure (interior) in order to determine the relative stability of interior and exterior nanoparticles. As shown in Figure 2, we will compare the Pt catalyst particle size at the beginning of test (BOT) and end of test (EOT), and discuss the importance of automating the electron tomography workflow, i.e. acquisition, reconstruction, and visualization, to increase sampling and determine the standard deviation of these measurement. The outlook for utilizing low-dose cryo-tomography for limiting damage to the catalyst, support, and especially proton-conducting ionomer will also be discussed [3].

Amichi, Lynda↗
NASA NTRS2018

Computational Model Prediction and Biological Validation Using Simplified Mixed Field Exposures for the Development of a GCR Reference Field

The yield of chromosomal aberrations has been shown to increase in the lymphocytes of astronauts after long-duration missions of several months in space. Chromosome exchanges, especially translocations, are positively correlated with many cancers and are therefore a potential biomarker of cancer risk associated with radiation exposure. Although extensive studies have been carried out on the induction of chromosomal aberrations by low- and high-LET radiation in human lymphocytes, fibroblasts, and epithelial cells exposed in vitro, there is a lack of data on chromosome aberrations induced by low dose-rate chronic exposure and mixed field beams such as those expected in space. Chromosome aberration studies at NSRL will provide the biological validation needed to extend the computational models over a broader range of experimental conditions (more complicated mixed fields leading up to the galactic cosmic rays (GCR) simulator), helping to reduce uncertainties in radiation quality effects and dose-rate dependence in cancer risk models. These models can then be used to answer some of the open questions regarding requirements for a full GCR reference field, including particle type and number, energy, dose rate, and delivery order. In this study, we designed a simplified mixed field beam with a combination of proton, helium, oxygen, and iron ions with shielding or proton, helium, oxygen, and titanium without shielding. Human fibroblasts cells were irradiated with these mixed field beam as well as each single beam with acute and chronic dose rate, and chromosome aberrations (CA) were measured with 3-color fluorescent in situ hybridization (FISH) chromosome painting methods. Frequency and type of CA induced with acute dose rate and chronic dose rates with single and mixed field beam will be discussed. A computational chromosome and radiation-induced DNA damage model, BDSTRACKS (Biological Damage by Stochastic Tracks), was updated to simulate various types of CA induced by acute exposures of the mixed field beams used for the experiments. The chromosomes were simulated by a polymer random walk algorithm with restrictions to their respective domains in the nucleus [1]. The stochastic dose to the nucleus was calculated with the code RITRACKS [2]. Irradiation of a target volume by a mixed field of ions was implemented within RITRACKs, and the fields of ions can be delivered over specific periods of time, allowing the simulation of dose-rate effects. Similarly, particles of various types and energies extracted from a pre-calculated spectra of galactic cosmic rays (GCR) can be used in RITRACKS. The number and spatial location of DSBs (DNA double-strand breaks) were calculated in BDSTRACKS using the simulated chromosomes and local (voxel) dose. Assuming that DSBs led to chromosome breaks, and simulating the rejoining of damaged chromosomes occurring during repair, BDSTRACKS produces the yield of various types of chromosome aberrations as a function of time (only final yields are presented). A comparison between experimental and simulation results will be shown.

Hada, M.↗
DOE OSTI2022

Evaluation of Gamma Ray Production - Quarterly Report FY22Q2

The LLNL team has performed the first systematic study of the effect of two-body weak currents in β decay rates across the entire mass table. In light and medium-mass nuclei, two-body weak currents from chiral effective field theory account for a significant portion of the phenomenological quenching of Gamow-Teller transition matrix elements. We examined the systematic effects of two-body axial currents on Gamow-Teller strength and β-decay rates in heavy nuclei within energy-density functional theory. We extended the charge-changing finiteamplitude pnFAM code to incorporate the contributions of two-body currents to the usual onebody linear response in the Gamow-Teller channel. We found that the two-body currents, as expected, usually quench both summed Gamow-Teller strength and decay rates, but by an amount that decreases as the neutron excess grows; see Fig. 1. In addition, they can enhance individual low-lying transitions, leading to decay rates that are quite different from those that an energy-independent quenching would produce, particularly in neutron-rich nuclei. This new physics capability is unique and provides the highest-fidelity description of β decay on the market. The underlying computational implementation can be reused to calculate γ-strength functions for finite-range nuclear potentials such as the Gogny force, which will allow us to benchmark some of our future developments against published results. The paper is in print in Phys. Rev. C.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1988

Numerical investigations in three-dimensional internal flows

Renewed interest in simulation of high enthalpy flows was brought about by interest in the design and testing of hypersonic propulsion components. This renewed interest has led to refurbishment of high speed ground based wind tunnel facilities. Simulation of the high enthalpy flows occurring in the NASP propulsion system is of particular interest. Studying the combustor presentes some of the more challenging problems in the ground based simulation of the fluid flow. A difficulty in simulating the required conditions at the entrance to the combustor is related to achieveing a sufficiently high static pressure such that combustion would be possible at a high Mach number. An analysis of the flow within the original nozzle was not performed previously using numerical simulations, although it was operated experimentally. A Navier-Stokes code was chosen in the present study to examine analytically the flow within representative selected nozzles for use in combustion experiments. The numerical simulation of these fluid flows is popular and may be able to shed light on the direction that redesign of the shock tube nozzles should take. The numerical codes were modified to account for the axially symmetric flow field, for the inclusion of equilibrium real gas effects and for ideal gases whose thermodynamic properties are frozen at the stagnation conditions. Numerical simulation of flows within the original and two proposed nozzles, including the energy lost to the relatively cold walls, are investigated.

Rose, William C.↗
NASA NTRS1992

Direct numerical simulation of a combusting droplet with convection

The evaporation and combustion of a single droplet under forced and natural convection was studied numerically from first principles using a numerical scheme that solves the time-dependent multiphase and multispecies Navier-Stokes equations and tracks the sharp gas-liquid interface cutting across an arbitrary Eulerian grid. The flow fields both inside and outside of the droplet are resolved in a unified fashion. Additional governing equations model the interphase mass, energy, and momentum exchange. Test cases involving iso-octane, n-hexane, and n-propanol droplets show reasonable comparison rate, and flame stand-off distance. The partially validated code is, thus, readied to be applied to more demanding droplet combustion situations where substantial drop deformation render classical models inadequate.

Liang, Pak-Yan↗
NASA NTRS2004

Research on the Natural Variability of Climate

The scheme for ocean vertical mixing presented in Canuto, Howard, Cheng & Dubovikov and Canuto, Howard, Hogan, Cheng, Dubovikov & Montenegro [GISS mixing scheme] was extended to include the latitudinal dependence reported by Gregg, Sanford & Winkel of the input to the interior ocean turbulence field of energy from internal waves. The resulting latitude dependence was implemented in our module of code for the GISS vertical mixing scheme and tested in the global NCAR-CSM ocean model with 3 degree X 3 degree resolution and 25 levels in the vertical. With the latitude dependence, the background diffusivity decreases from approx. 0.1 sq cm/sec at midlatitudes to approx. less than 0.01 sq cm/sec at the equator. A significant improvement was seen in the tropical ocean, with tropical thermoclines being sharpened in agreement with the requirements of observational studies and the needs of ENS0 modeling. At the same time global measures nf performance such as the meridional overturning and northward heat transport were not adversely affected. Results were presented at the 2004 Layered Ocean Modeling meeting and a paper has been prepared and submitted to Geophysical Research Letters.

Source record↗
NASA NTRS1986

Computer simulations of electromagnetic cool ion beam instabilities

Electromagnetic ion beam instabilities driven by cool ion beams at propagation parallel or antiparallel to a uniform magnetic field are studied using computer simulations. The elements of linear theory applicable to electromagnetic ion beam instabilities and the simulations derived from a one-dimensional hybrid computer code are described. The quasi-linear regime of the right-hand resonant ion beam instability, and the gyrophase bunching of the nonlinear regime of the right-hand resonant and nonresonant instabilities are examined. It is detected that in the quasi-linear regime the instability saturation is due to a reduction in the beam core relative drift speed and an increase in the perpendicular-to-parallel beam temperature; in the nonlinear regime the instabilities saturate when half the initial beam drift kinetic energy density is converted to fluctuating magnetic field energy density.

Gary, S. P.↗
DOE OSTI2024

Single-Particle Soot Photometer (SP2) Black Carbon Number and Mass Concentrations

The single-particle soot photometer (SP2) records particle-by-particle measurements of the intensity of both the scattering signature and incandescence signature of particles that enter its laser beam. These intensities are then used to calculate refractory black carbon (rBC) masses and particle diameters. In previous U.S. Department of Energy Atmospheric Radiation Measurement (ARM) user facility field campaigns, the SP2 data was difficult to process because the manufacturer-supplied code was not scalable to distributed machines, making it unusable for the large amounts of data output by the SP2. This prohibited the SP2 from becoming an operational instrument for ARM. Therefore, PySP2 was developed to solve this issue and enable SP2 to be an operational instrument for ARM. This technical document summarizes the test data sets from the ARM North Slope of Alaska (NSA) site and the Multidisciplinary Drifting Observatory for the Study of Arctic Climate (MOSAiC) field campaign that were used to develop and test PySP2.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2016

Enhancement of the Open National Combustion Code (OpenNCC) and Initial Simulation of Energy Efficient Engine Combustor

In this paper, we present the recent enhancement of the Open National Combustion Code (OpenNCC) and apply the OpenNCC to model a realistic combustor configuration (Energy Efficient Engine (E3)). First, we perform a series of validation tests for the newly-implemented advection upstream splitting method (AUSM) and the extended version of the AUSM-family schemes (AUSM+-up). Compared with the analytical/experimental data of the validation tests, we achieved good agreement. In the steady-state E3 cold flow results using the Reynolds-averaged Navier-Stokes(RANS), we find a noticeable difference in the flow fields calculated by the two different numerical schemes, the standard Jameson- Schmidt-Turkel (JST) scheme and the AUSM scheme. The main differences are that the AUSM scheme is less numerical dissipative and it predicts much stronger reverse flow in the recirculation zone. This study indicates that two schemes could show different flame-holding predictions and overall flame structures.

Code validation↗