Search NASA⌕ Search

SEARCH · Search NASA

Results for “Code Robustness”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

HFIRCON Version 1.0.5 User Guide

The High Flux Isotope Reactor (HFIR) Controller (HFIRCON) code is a collection of python routines and C plugins that automate the workflow for fuel and single- or multicycle target depletion analyses for HFIR at Oak Ridge National Laboratory (ORNL). This code calls the LAVAMINT (LAVA Model Interrogator) parallel (MCNP) Monte Carlo N-Particle model interrogator to stochastically calculate cell volumes and bounding boxes, the ADVANTG (Automatic Variance Reduction Generation) code package for all variance reduction and source biasing calculations, the ORNL-Transformative Neutronics/MCNP5 transport solver for all transport solutions, and the MSX_DEPLETE module to perform all depletion calculations via the ORIGEN (Oak Ridge Isotope Generation) application programming interface. It also performs a robust set of postprocessing functions to automatically provide summaries of several key metrics that are common to a wide variety of typical HFIR design and safety-basis analyses.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

A robust spectral element implementation of the $k - τ$ RANS model in Nek5000/NekRS

The $k - ω$ Reynolds Averaged Navier Stokes (RANS) model is one of the industry standard approaches for modeling of turbulent flows. It performs better than the $k - ϵ$ model for low Reynolds number flows and is also more suitable for boundary layers with adverse pressure gradients. Major drawback of the model, however, is that the asymptotic value of $ω$ at the walls is singular, necessitating the use of a contrived “sufficiently” large value for $ω$ as the boundary condition for its transport equation. Here, this invariably leads to the solution being sensitive to near wall grid spacing. While an acceptable solution for low order (finite volume) methods, the excessive near wall gradients lead to persistent numerical stability issues in high order codes. To alleviate the problem, specifically in the context of the high order spectral element code Nek5000, a regularized $k - ω$ approach was formulated in our prior work (Tomboulides et al., 2018). The formulation, however, relies on the use of wall distance and its gradients for modeling the closure terms and can pose problems for simulations in complex geometries. This work presents a novel implementation of the $k - τ$ RANS model in Nek5000, where $τ = 1/ω$, eliminating the need for regularization, owing to the asymptotically bounded behavior of the source terms in the $τ$ transport equation, and also eliminating dependence on wall distance. Robustness and stability of the $k - τ$ model is ensured through implicit treatment of the source terms and their careful numerical implementation and demonstrated through several cases aimed at verification and validation. Studies include both canonical and engineering relevant problems, viz., turbulent channel flow, pipe flow, backward facing step, flow over NACA 0012 airfoil and flow in a T-junction. Results from the $k - τ$ model are shown to be consistent with regularized $k - ω$ model and also with the $k - ω$ SST model in OpenFOAM (for select studies). Comparison with experimental data is also shown, where available, to bolster validation efforts for the $k - τ$ model implementation through prediction of key turbulent quantities of interest.

Nek5000↗
DOE OSTI2020

Graph interpolating activation improves both natural and robust accuracies in data-efficient deep learning

Improving the accuracy and robustness of deep neural nets (DNNs) and adapting them to small training data are primary tasks in deep learning (DL) research. In this paper, we replace the output activation function of DNNs, typically the data-agnostic softmax function, with a graph Laplacian-based high-dimensional interpolating function which, in the continuum limit, converges to the solution of a Laplace–Beltrami equation on a high-dimensional manifold. Furthermore, we propose end-to-end training and testing algorithms for this new architecture. The proposed DNN with graph interpolating activation integrates the advantages of both deep learning and manifold learning. Compared to the conventional DNNs with the softmax function as output activation, the new framework demonstrates the following major advantages: First, it is better applicable to data-efficient learning in which we train high capacity DNNs without using a large number of training data. Second, it remarkably improves both natural accuracy on the clean images and robust accuracy on the adversarial images crafted by both white-box and black-box adversarial attacks. Third, it is a natural choice for semi-supervised learning. This paper is a significant extension of our earlier work published in NeurIPS, 2018. For reproducibility, the code is available at https://github.com/BaoWangMath/DNN-DataDependentActivation .

Mathematics↗
DOE OSTI2015

An ALMA survey of sub-millimeter galaxies in the extended Chandra deep field south: Physical properties derived from ultraviolet-to-radio modeling

The ALESS survey has followed up on a sample of 122 sub-millimeter sources in the Extended Chandra Deep Field South at 870 μm with the Atacama Large Millimeter Array (ALMA), allowing us to pinpoint the positions of sub-millimeter galaxies (SMGs) to ∼0.3 arcsec and to find their precise counterparts at different wavelengths. This enabled the first compilation of the multi-wavelength spectral energy distributions (SEDs) of a statistically reliable survey of SMGs. In this paper, we present a new calibration of the magphys SED modeling code that is optimized to fit these ultraviolet-to-radio SEDs of z>1 star-forming galaxies using an energy balance technique to connect the emission from stellar populations, dust attenuation, and dust emission in a physically consistent way. We derive statistically and physically robust estimates of the photometric redshifts and physical parameters (such as stellar masses, dust attenuation, star formation rates (SFRs), and dust masses) for the ALESS SMGs. We find that the ALESS SMGs have median stellar mass M{sub ∗}=8.9±0.1×10{sup 10} M{sub ⊙}, median SFR = 280±70 M{sub ⊙} yr{sup −1}, median overall V-band dust attenuation A{sub V}=1.9±0.2 mag, median dust mass M{sub dust}=(5.6±1.0)×10{sup 8} M{sub ⊙}, and median average dust temperature T{sub dust}≃40 K. We find that the average intrinsic SED of the ALESS SMGs resembles that of local ultra-luminous infrared galaxies in the infrared range, but the stellar emission of our average SMG is brighter and bluer, indicating lower dust attenuation, possibly because they are more extended. We explore how the average SEDs vary with different parameters (redshift, sub-millimeter flux, dust attenuation, and total infrared luminosity), and we provide a new set of SMG templates that can be used to interpret other SMG observations. To put the ALESS SMGs into context, we compare their stellar masses and SFRs with those of less actively star-forming galaxies at the same redshifts. We find that at z≃2, about half of the SMGs lie above the star-forming main sequence (with SFRs three times larger than normal galaxies of the same stellar mass), while half are consistent with being at the high-mass end of the main sequence. At higher redshifts (z≃3.5), the SMGs tend to have higher SFRs and stellar masses, but the fraction of SMGs that lie significantly above the main sequence decreases to less than a third.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

A new moment-based general-relativistic neutrino-radiation transport code: Methods and first applications to neutron star mergers

ABSTRACT We present a new moment-based energy-integrated neutrino transport code for neutron star merger simulations in general relativity. In the merger context, ours is the first code to include Doppler effects at all orders in υ/c, retaining all non-linear neutrino–matter coupling terms. The code is validated with a stringent series of tests. We show that the inclusion of full neutrino–matter coupling terms is necessary to correctly capture the trapping of neutrinos in relativistically moving media, such as in differentially rotating merger remnants. We perform preliminary simulations proving the robustness of the scheme in simulating ab-initio mergers to black hole collapse and long-term neutron star remnants up to ${\sim }70\,$ ms. The latter is the longest dynamical space-time, 3D, general relativistic simulations with full neutrino transport to date. We compare results obtained at different resolutions and using two different closures for the moment scheme. We do not find evidences of significant out-of-thermodynamic equilibrium effects, such as bulk viscosity, on the post-merger dynamics or gravitational wave emission. Neutrino luminosities and average energies are in good agreement with theory expectations and previous simulations by other groups using similar schemes. We compare dynamical and early wind ejecta properties obtained with M1 and with our older neutrino treatment. We find that the M1 results have systematically larger proton fractions. However, the differences in the nucleosynthesis yields are modest. This work sets the basis for future detailed studies spanning a wider set of neutrino reactions, binaries, and equations of state.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Was There a 3.5 keV Line?

Abstract The 3.5 keV line is a purported emission line observed in galaxies, galaxy clusters, and the Milky Way whose origin is inconsistent with known atomic transitions and has previously been suggested to arise from dark matter decay. We systematically reexamine the bulk of the evidence for the 3.5 keV line, attempting to reproduce six previous analyses that found evidence for the line. We only reproduce one of the analyses; in the other five, we find no significant evidence for a 3.5 keV line when following the described analysis procedures on the original data sets. For example, previous results claimed 4 σ evidence for a 3.5 keV line from the Perseus cluster; we dispute this claim, finding no evidence for a 3.5 keV line. We find evidence for background mismodeling in multiple analyses. We show that analyzing these data in narrower energy windows diminishes the effects of mismodeling but returns no evidence for a 3.5 keV line. We conclude that there is little robust evidence for the existence of the 3.5 keV line. Some of the discrepancy of our results from those of the original works may be due to the earlier reliance on local optimizers, which we demonstrate can lead to incorrect results. For ease of reproducibility, all code and data are publicly available.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Robust estimation of charge carrier diffusivity using transient photoluminescence microscopy

Transient microscopy has emerged as a powerful tool for imaging the diffusion of excitons and free charge carriers in optoelectronic materials. In many excitonic materials, extraction of diffusion coefficients can be simplified because of the linear relationship between signal intensity and local excited state population. However, in materials where transport is dominated by free charge carriers, extracting diffusivities accurately from multidimensional data is complicated by the nonlinear dependence of the measured signal on the local charge carrier density. To obtain accurate estimates of charge carrier diffusivity from transient microscopy data, statistically robust fitting algorithms coupled to efficient 3D numerical solvers that faithfully relate local carrier dynamics to raw experimental measurables are sometimes needed. Here, we provide a detailed numerical framework for modeling the spatiotemporal dynamics of free charge carriers in bulk semiconductors with significant solving speed reduction and for simulating the corresponding transient photoluminescence microscopy data. To demonstrate the utility of this approach, we apply a fitting algorithm using a Markov chain Monte Carlo sampler to experimental data on bulk CdS and methylammonium lead bromide (MAPbBr 3 ) crystals. Parameter analyses reveal that transient photoluminescence microscopy can be used to obtain robust estimates of charge carrier diffusivities in optoelectronic materials of interest, but that other experimental approaches should be used for obtaining carrier recombination constants. Additionally, simplifications can be made to the fitting model depending on the experimental conditions and material systems studied. Our open-source simulation code and fitting algorithm are made freely available to the scientific community.

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

Sensitivity and Uncertainty Workflow of Full System SIERRA Models Supporting High Consequence Applications (Executive Summary)

The overall goal of this work was to advance the integrated workflow capabilities to provide the weapon system analysts the tools to construct, communicate and robustly execute end-to-end computational simulation models and analysis. Specifically, this includes developing and running automated, CompSim workflows that map model inputs to responses from large-scale studies executed on a distributed, heterogenous computing environment (e.g., CAD applications on Windows laptop to Sierra FEM codes on Trinity within a single workflow). The workflows developed enable ensemble calculations in support of ND program decisions, including sensitivity analysis, optimization and uncertainty quantification.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

A Riccati Solution for the Ideal MHD Plasma Response with Applications to Real-time Stability Control

To effectuate near real-time feedback control of ideal MHD instabilities in a tokamak geometry, a rapid solution for stability analysis is a prerequisite. Toward this end, we reformulate the δW stability method with a Hamilton-Jacobi theory, elucidating analytical and numerical features of the generic tokamak ideal MHD stability problem. The plasma response matrix is demonstrated to be the solution of an ideal MHD matrix Riccati differential equation (MRDE). Since Riccati equations are prevalent in the control theory literature, such a shift in perspective brings to bear a range of numerical methods that are well-suited to the robust, fast solution of control problems. We discuss the usefulness of Riccati techniques in solving the stiff ODEs often encountered in ideal MHD stability analyses-—for example, in tokamak edge and stellarator physics. We then demonstrate the applicability of such methods to an existing 2D ideal MHD stability code—DCON—enabling its parallel operation in near real-time. Output is shown to match with high accuracy, and with wall-clock time ≪ 1s. Such speed may help enable active feedback ideal MHD stability control, especially in tokamak plasmas whose ideal MHD equilibria evolve with inductive timescale τ > 1s-—as in ITER.

Parallel↗
DOE OSTI2016

Quantitative comparison of electron temperature fluctuations to nonlinear gyrokinetic simulations in C-Mod Ohmic L-mode discharges

Long wavelength turbulent electron temperature fluctuations (k{sub y}ρ{sub s} < 0.3) are measured in the outer core region (r/a > 0.8) of Ohmic L-mode plasmas at Alcator C-Mod [E. S. Marmar et al., Nucl. Fusion 49, 104014 (2009)] with a correlation electron cyclotron emission diagnostic. The relative amplitude and frequency spectrum of the fluctuations are compared quantitatively with nonlinear gyrokinetic simulations using the GYRO code [J. Candy and R. E. Waltz, J. Comput. Phys. 186, 545 (2003)] in two different confinement regimes: linear Ohmic confinement (LOC) regime and saturated Ohmic confinement (SOC) regime. When comparing experiment with nonlinear simulations, it is found that local, electrostatic ion-scale simulations (k{sub y}ρ{sub s} ≲ 1.7) performed at r/a ∼ 0.85 reproduce the experimental ion heat flux levels, electron temperature fluctuation levels, and frequency spectra within experimental error bars. In contrast, the electron heat flux is robustly under-predicted and cannot be recovered by using scans of the simulation inputs within error bars or by using global simulations. If both the ion heat flux and the measured temperature fluctuations are attributed predominantly to long-wavelength turbulence, then under-prediction of electron heat flux strongly suggests that electron scale turbulence is important for transport in C-Mod Ohmic L-mode discharges. In addition, no evidence is found from linear or nonlinear simulations for a clear transition from trapped electron mode to ion temperature gradient turbulence across the LOC/SOC transition, and also there is no evidence in these Ohmic L-mode plasmas of the “Transport Shortfall” [C. Holland et al., Phys. Plasmas 16, 052301 (2009)].

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Chasing ghosts: characterization of artifact generation in coded aperture decoding due to experimental implementation

Coded aperture imaging is a form of lensless aperture imaging that projects multiple overlapping images of the source onto the detector, enhancing signal strength, which is advantageous for low-flux sources or high-resolution imaging. This technique requires decoding of the detector signal to reconstruct the original source, which involves convolving the detector data with the aperture pattern. When the signal is from a centered point source, the reconstructed source image is known as the point spread function (PSF). A clean PSF without artifacts is a Dirac delta function [Appl. Opt. 20, 1858 (1981)]. This paper examines the robustness of the decoding process against variations in experimental tolerances by analyzing artifact growth in the reconstructed PSF. We illustrate the effects of incorrect magnification, rotation, and detector size and find that aperture–detector rotational misalignment about the imaging axis is the most sensitive parameter, with significant artifact generation occurring with angular offsets of less than one degree. We discuss compensation methods for imperfect aperture placement, finding that small detector sizes produce uncompensatable artifact generation, and compare theoretical predictions with experimental PSF measurements of a rank , 6.8 mm thick (less than one mean free path) coded aperture with a 3.5 mm cell size, conducted at the MegaJOuLe Neutron Imaging Radiography dense plasma focus [IEEE Trans. Plasma Sci. 49, 3299 (2021)] using a 2.45 MeV neutron source. Based on our findings, we recommend using magnified coded apertures in the under-sampled regime, which allows for the inclusion of fiducial markers to characterize aperture–detector rotational offsets and the addition of mechanical coupling, where possible, to constrain rotational and magnification offsets.

Selwood, M. P. [Lawrence Livermore National Labora↗
DOE OSTI2025

A direct detection method of galaxy intrinsic ellipticity-gravitational shear correlation in non-linear regimes using self-calibration

Intrinsic alignment (IA) of galaxies is a challenging source of contamination in the Cosmic shear (GG) signals. The galaxy intrinsic ellipticity-gravitational shear (IG) correlation is generally the most dominant component of such contamination for cross-correlating redshift bins. One of the most effective techniques to mitigate such contamination is the self-calibration (SC) method which extracts the IG correlation and allows for its removal from the GG signal. In a photometric survey, the SC method first extracts the galaxy number density-galaxy intrinsic ellipticity (gI) correlation from the observed galaxy-galaxy lensing correlation using the redshift dependence of lens-source pairs. The IG correlation is computed through a scaling relation using the gI correlation and other lensing observables. The applicability of the SC method has so far been focused on the linear IA scales and the linear galaxy bias. We extend the SC method beyond the linear regime by modifying its scaling relation which can account for the non-linear galaxy bias model and various IA models. In this study, we provide a framework to detect the IG correlation for the redshift bins for source galaxies for the proposed year 1 survey of the Rubin Legacy Survey of Space and Time (LSST Y1). We tested the method for the tidal alignment and tidal torquing (TATT) model of IA and we found that the scaling relation is accurate within 10% and 20% for cross-correlating and auto-correlating redshift bins, respectively. Hence the suppression of IG contamination in observed GG correlation can be accomplished with a factor of 10 and 5, for cross-correlating and auto-correlating redshift bins, respectively. We tested the method's robustness and found that the suppression of IG contamination by a factor of 5 is still achievable for all combinations of cross-correlating bins even with the inclusion of a moderate amount of uncertainties on IA and bias parameters, respectively. We also make available, a branch of the code FAST-PT to provide gI correlations up to 1-loop order term used by the new SC method.

gravitational lensing↗
DOE OSTI2024

Polarization–frequency hyperentangled photons: generation, characterization, and manipulation

Frequency-bin encoding is massively parallelizable and robust for optical fiber transmission. When coupled with an additional degree of freedom (DoF), the expansion of the Hilbert space allows for deterministic controlled operations between two DoFs within a single photon. Such capabilities, when combined with photonic hyperentanglement, are of great value for quantum communication protocols, including dense coding and single-copy entanglement distillation. In this talk, we present an all-fiber-coupled, ultrabroadband polarization–frequency hyperentangled source and conduct comprehensive quantum state tomography across multiple dense wavelength division multiplexing channels spanning the optical C+L-band (1530–1625 nm). In addition, we design and implement a high-fidelity controlled-NOT (cnot) operation between polarization and frequency DoFs by exploiting electro-optic phase modulation within a fiber Sagnac loop. Collectively, our hyperentangled source and two-qubit gate should unlock new opportunities for harnessing polarization–frequency resources in established telecommunication fiber networks for future quantum applications.

Lu, Hsuan-Hao↗
DOE OSTI2021

Upgrade of hardware controls for the STAR experiment at RHIC

The STAR experiment has been delivering significant physics results for more than 20 years. Stable operation of the experiment was achieved by using a robust controls system based on the Experimental Physics and Industrial Control System (EPICS). Now an object-oriented approach with Python libraries, adapted for EPICS software, is going to replace the procedural-based EPICS C libraries previously used at STAR. Advantages of the new approach include stability of operation, code reduction and straightforward project documentation. The first two sections of this paper introduce the STAR experiment, give an overview of the EPICS architecture, and present the use of Python for controls software. Therefore, specific examples, as well as upgrades of user interfaces, are outlined in the following sections.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Smoothing and flattening the universe through slow contraction versus inflation

In a systematic study, we use an equivalent pair of improved numerical relativity codes based on a tetrad-formulation of the classical Einstein-scalar field equations to examine whether slow contraction or inflation (or both) can resolve the homogeneity, isotropy and flatness problems. Our finding, based on a set of gauge/frame invariant diagnostics and the models considered, is that slow contraction robustly and rapidly smooths and flattens spacetime beginning from initial conditions that are outside the perturbative regime of the flat Friedmann-Robertson-Walker metric, whereas inflation fails these tests. We present new numerical evidence supporting the conjecture that the combination of ultralocal evolution and an effective equation-of-state with pressure much greater than energy density is the key to having robust and rapid smoothing. The opposite of ultralocality occurs in expanding spacetimes, which is the leading obstruction to smoothing following a big bang.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

SciDAC – PSI: Plasma Surface Interactions involving Tungsten - A Summary of Project Accomplishments

Gaining physics understanding and predictive capabilities to describe the evolution of plasma facing components (PFC) requires simultaneously addressing complex and diverse physics occurring over a wide range of length and time scales, as well as integrating extensive physical processes across the plasma - surface - bulk materials boundaries. This requires development not only of detailed physics models and computational strategies at each of these scales, but computer science algorithms and methods to strongly couple them in a way that can be robustly validated through comparisons against available data and new experiments. Therefore, the objective of this project is to develop robust, high-fidelity simulation tools capable of predicting the PFC operating lifetime and the PFC impact on plasma contamination, recycling of hydrogenic species, and tritium retention in future magnetic fusion devices, with a focus on tungsten based material systems. Deploying these tools requires the development of a leadershipscale computational code, as well as a host of simulations that span the multiple scales needed to address complex physical and computational issues at the plasma - surface interface and the transition below the surface where neutron damage processes in the bulk material dominate behavior in multiple-component materials systems. Successful development will enable improved prediction of PFC performance needed to ensure magnetic fusion energy development beyond ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

EXTINCTION DROPLET SIZING MEASUREMENTS IN DIESEL RELEVANT SPRAYS

A new diagnostic for the quantification of Sauter Mean Diameter (SMD) in high-pressure fuel sprays has been recently developed using combined optical and x-ray measurements at the Georgia Institute of Technology and Argonne National Laboratory, respectively. This diagnostic utilizes liquid scattering extinction measurements from diffuse back-illumination (DBI) imaging, conducted at Georgia Tech, and liquid absorption measurements from x-ray radiography, conducted at Argonne’s Advanced Photon Source. The new diagnostic, entitled the Scattering Absorption Measurement Ratio (SAMR), quantifies two-dimensional distributions of path-integrated SMD, enabling construction of the spatial history of drop size development within practical fuel sprays. This technique offers unique benefits over conventional drop-sizing methods in that it can be more robust in optically dense regions of the spray, while also providing high spatial resolution of the corresponding droplet field. The spatially-resolved SMD measurements that result from the SAMR diagnostic will be especially valuable to the engine modeling community for the quantitative validation of spray submodels in engine CFD codes. The methodology for quantification of SMD distributions using the SAMR technique has been previously introduced. This thesis aims to extend the initial development of the SAMR technique by presenting in detail: the experimental methodologies used in the SAMR technique, including the development of an ideal DBI setup for this technique; the data processing methodologies developed by the author; two-dimensional Sauter Mean Diameter measurements within diesel sprays for various experimental conditions; a summary of the various sources of measurement uncertainty; an assessment for how the sources of uncertainty affect the quantified SMD; and paths to improve the technique moving forward. The SMD results show that for low ambient density conditions droplets decrease in size as radial and axial position increases. For high ambient density conditions, however, the droplets show a stable size near the spray centerline and steadily increase in size as distance from the centerline increases. Additionally, this work provides 2-D droplet sizing maps which are particularly useful for model validation.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2018

EXTINCTION DROPLET SIZING MEASUREMENTS IN DIESEL RELEVANT SPRAYS

A new diagnostic for the quantification of Sauter Mean Diameter (SMD) in high-pressure fuel sprays has been recently developed using combined optical and x-ray measurements at the Georgia Institute of Technology and Argonne National Laboratory, respectively. This diagnostic utilizes liquid scattering extinction measurements from diffuse back-illumination (DBI) imaging, conducted at Georgia Tech, and liquid absorption measurements from x-ray radiography, conducted at Argonne’s Advanced Photon Source. The new diagnostic, entitled the Scattering Absorption Measurement Ratio (SAMR), quantifies two-dimensional distributions of path-integrated SMD, enabling construction of the spatial history of drop size development within practical fuel sprays. This technique offers unique benefits over conventional drop-sizing methods in that it can be more robust in optically dense regions of the spray, while also providing high spatial resolution of the corresponding droplet field. The spatially-resolved SMD measurements that result from the SAMR diagnostic will be especially valuable to the engine modeling community for the quantitative validation of spray submodels in engine CFD codes. The methodology for quantification of SMD distributions using the SAMR technique has been previously introduced. This thesis aims to extend the initial developmen of SAMR technique by presenting in detail: the experimental methodologies used in the SAMR technique, including the development of an ideal DBI setup for this technique; the data processing methodologies developed by the author; two-dimensional Sauter Mean Diameter measurements within diesel sprays for various experimental conditions; a summary of the various sources of measurement uncertainty; an assessment for how the sources of uncertainty affect the quantified SMD; and paths to improve the technique moving forward. The SMD results show that for low ambient density conditions droplets decrease in size as radial and axial position increases. For high ambient density conditions, however, the droplets show a stable size near the spray centerline and steadily increase in size as distance from the centerline increases. Additionally, this work provides 2-D droplet sizing maps which are particularly useful for model validation.

33 ADVANCED PROPULSION SYSTEMS↗