Search NASA⌕ Search

SEARCH · Search NASA

Results for “generative diffusion model”

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 577 records · Page 32

Application of the pressure sensitive paint technique to steady and unsteady flow

Pressure sensitive paint is a newly-developed optical measurement technique with which one can get a continuous pressure distribution in much shorter time and lower cost than a conventional pressure tap measurement. However, most of the current pressure sensitive paint applications are restricted to steady pressure measurement at high speeds because of the small signal-to-noise ratio at low speed and a slow response to pressure changes. In the present study, three phases of work have been completed to extend the application of the pressure sensitive paint technique to low-speed testing and to investigate the applicability of the paint technique to unsteady flow. First the measurement system using a commercially available PtOEP/GP-197 pressure sensitive paint was established and applied to impinging jet measurements. An in-situ calibration using only five pressure tap data points was applied and the results showed good repeatability and good agreement with conventional pressure tap measurements on the whole painted area. The overall measurement accuracy in these experiments was found to be within 0.1 psi. The pressure sensitive paint technique was then applied to low-speed wind tunnel tests using a 60 deg delta wing model with leading edge blowing slots. The technical problems encountered in low-speed testing were resolved by using a high grade CCD camera and applying corrections to improve the measurement accuracy. Even at 35 m/s, the paint data not only agreed well with conventional pressure tap measurements but also clearly showed the suction region generated by the leading edge vortices. The vortex breakdown was also detected at alpha=30 deg. It was found that a pressure difference of 0.2 psi was required for a quantitative pressure measurement in this experiment and that temperature control or a parallel temperature measurement is necessary if thermal uniformity does not hold on the model. Finally, the pressure sensitive paint was applied to a periodically changing pressure field with a 12.8s time period. A simple first-order pole model was applied to deal with the phase lag of the paint. The unsteady pressure estimated from the time-changing pressure sensitive paint data agreed well with the pressure transducer data in regions of higher pressure and showed the possibility of extending the technique to unsteady pressure measurements. However, the model still needs further refinement based on the physics of the oxygen diffusion into the paint layer and the oxygen quenching on the paint luminescence.

Shimbo, Y.↗

Li1−xNiO2 Many-body DMC Benchmark Dataset

The dataset contains all numerical data generated in support of the manuscript “Many‑body Benchmark of Electronic Charge and Spin Densities for Li1–xNiO2​” (Journal of Chemical Theory and Computation, DOI: 10.1021/acs.jctc.5c02097, URL: https://pubs.acs.org/doi/10.1021/acs.jctc.5c02097). The materials included in this repository are: 1. Data files used to produce all figures and tables in the main manuscript and supporting information. 2. Benchmark density‑functional theory (DFT) datasets used for the charge‑ and spin‑density analyses. 3. Reference many‑body diffusion Monte Carlo (DMC) calculations and associated input/output files.

36 MATERIALS SCIENCE↗

Nonadiabatic Force Matching for Alchemical Free-Energy Estimation

We propose a method to compute free-energy differences from nonadiabatic alchemical transformations by using flow-based generative models. The method, nonadiabatic force matching, hinges on estimating the dissipation along an alchemical switching process in terms of a nonadiabatic force field that can be learned through stochastic flow matching. The learned field can be used in conjunction with short-time trajectory data to evaluate upper and lower bounds on the alchemical free energy that variationally converge to the exact value if the field is optimal. Applying the method to evaluate the alchemical free energy of atomistic models shows that it can substantially reduce the simulation cost of a free-energy estimate at a negligible loss of accuracy when compared with thermodynamic integration.

Computational chemistry↗

Tailoring composition and deformation modes at the microstructural level for next generation low-cost high-strength austenitic stainless steels

The objective of this project is to enable deliberate development of cost-effective, hydrogen resistant alloys by establishing detailed relationships specific to the effects of alloy composition, short-range order (SRO), and microsegregation in the presence of hydrogen on the transition between homogeneous deformation and localized plasticity in shear bands. In collaboration with the International Institute for Carbon-Neutral Energy Research, I2CNER, at Kyushu University in Japan, we conceptualized, designed, and manufactured four austenitic alloys that maintain corrosion resistance and ensure lower cost relative to baseline commercial alloys. The mechanical properties and deformation modes of the novel alloys (KU alloys) were assessed in the presence of hydrogen (H). Correlations between composition and performance revealed that two of the KU alloys are suitable replacements for 316 steel, while another is a viable replacement for 304 steel at room temperature. We found that, in the presence of other austenite stabilizing elements namely Mn and N, replacing Ni with Cu does not lead to martensite formation as has been previously reported.1–3 Furthermore, we found that the addition of Cu leads to an earlier onset of multiple slip resulting in an relative earlier onset of a higher work hardening rate (WHR). Greater understanding of the relationships between alloy composition and SRO required the development of a novel advanced electron diffraction methodology to characterize SRO in complex FCC alloys. This innovative approach, which combines fluctuation and correlation analyses of diffuse-scattering signals, successfully differentiated between SRO and long-range ordering (LRO). Further investigations into annealed austenitic stainless steels could provide insights into manipulating SRO and its effects on material properties. Atomistic simulations provided understanding of SRO behavior that was difficult to capture experimentally. This project created the first spin cluster expansion model that is able to capture and describe SRO effects in Fe-Ni-Cr FCC alloys, accounting for the non-negligible effects of magnetism. An automated computational workflow was established to provide reliable predictions of SRO in Fe-Ni-Cr austenitic alloys, both with and without the presence of H atoms. Analysis of the propensity for SRO in Fe-Ni-Cr alloys revealed that H tends to cluster with specific, well-defined SRO domains. The computational framework is general purpose and can be extended to realistic stainless steels across diverse composition ranges. With confidence that SRO is possible in austenitic stainless steels, we developed a discrete dislocation finite element code to understand the interaction of dislocations with SRO in the presence of H. By incorporating H effects on the dislocation emission and SRO stress field we show that the critical stress for the dislocation pileup to breakthrough the SRO domain decreases in the presence of H, which directly contributes localized deformation at the macroscale. Through the simulation of a uniaxial tension test, we demonstrated that H-induced weakening of SRO stress field and H-enhanced dislocation emission can lead to the onset of shear localization at lower macroscopic strains. As a whole, this project identified three novel alloys that show improvements in performance and cost efficiency for H-facing applications by studying correlations between alloy chemistry and deformation behavior. We also made significant advancements to experimental and computational methodologies necessary to study the chemistry and distribution of SRO across a range of alloys, which in turn allowed us to demonstrate how deformation mechanisms change due to the contributions of SRO in austenitic alloys in the presence of H. The combined advancements in fundamental understanding with novel alloy development in this project has increased the viability of next generation H-technologies for the broader public through accessible low-cost alloys and accelerated development towards future H-infrastructure.

08 HYDROGEN↗

Measuring the neutrino-oxygen neutral current quasielastic cross section using the accelerator neutrino neutron interaction experiment

The Accelerator Neutrino Neutron Interaction Experiment (ANNIE) is a 26-ton gadolinium-doped water Cherenkov detector located on-axis to Fermilab’s Booster Neutrino Beam (BNB). ANNIE is uniquely positioned to perform high-statistics measurements of neutrino-nucleus interactions in water, benefiting from a large neutrino flux due to a short (100-meter) baseline. A central focus of ANNIE’s physics program is the measurement of both charged current (CC) and neutral current (NC) cross sections on water, including neutral current quasielastic (NCQE) and CC-inclusive channels. The NCQE measurement is particularly critical for constraining uncertainties in rare-event searches such as the Diffuse Supernova Neutrino Background (DSNB), where atmospheric $\nu$NCQE interactions constitute a significant and poorly constrained background. This dissertation presents a measurement of the flux-averaged neutrino-oxygen neutral current quasielastic ($\nu$NCQE) cross section using $2.573 \times 10^{20}$~POT of BNB exposure from the 2022 and 2023 beam years. The $\nu$NCQE interaction is identified through the primary $\gamma$-rays produced by nuclear de-excitation of the residual $^{15}$N$^*$ or $^{15}$O$^*$ nucleus following nucleon knockout from $^{16}$O. A dedicated Monte Carlo (MC) re-tuning campaign was conducted using an americium-beryllium (AmBe) calibration source, Michel electrons from stopped muons, and throughgoing dirt muons originating upstream of the detector. This multi-sample approach provided a wide-ranging $\mathcal{O}(\text{MeV})$--$\mathcal{O}(\text{GeV})$ dataset for tuning the simulated detector response, which was subsequently validated against AmBe neutron and Michel electron data for use in the $\nu$NCQE analysis. A dedicated laser calibration campaign was carried out to reduce timing uncertainties across the PMT system, enabling reconstruction of the BNB bunch substructure with sufficient resolution to serve as a background rejection tool. By selecting events in-time with individual neutrino bunches, beam-correlated $\nu$NCQE events are separated from diffuse and accelerator-induced backgrounds, notably skyshine neutrons and externally-originating events, that would otherwise dominate traditional charge-based selections within a small-scale, surface-level, short-baseline detector. A data-driven estimation of the skyshine neutron and external background rates was performed and incorporated into the systematic uncertainty budget. The flux-averaged $\nu$NCQE cross section on oxygen is measured to be $1.57 \pm 0.06\,(\text{stat.})$ $^{+0.91}_{-0.67}\,(\text{syst.})$ $\times 10^{-38}\ \text{cm}^{2}$. A full systematic budget is constructed by propagating uncertainties in the secondary hadronic interaction modeling, background cross section normalizations, detector response, neutrino flux, and the primary $\gamma$-ray emission probabilities from oxygen nuclear de-excitation. An idealized de-excitation model, constructed from existing measurements in the literature is developed to benchmark the predictions of the \textsc{GENIE} event generator. A comparison reveals that \textsc{GENIE} systematically overpredicts the primary $\gamma$-ray emission probability from oxygen de-excitation by a factor of $1.49\times$ for $E_\gamma > 6$~MeV and $3.07\times$ in the $3$--$6$~MeV band. This comparison motivates the dominant systematic uncertainty in this analysis, where a conservative uncertainty of $^{+39.9\%}_{-0\%}$ on the primary $\gamma$-ray signal prediction is assigned. The ANNIE result is consistent with and complementary to existing flux-averaged $\nu$NCQE cross section measurements from T2K and Super-Kamiokande, providing an independent measurement with a different detector, neutrino beam, and analysis methodology. Looking ahead, an upgrade to the ANNIE DAQ infrastructure enabling continuous extended readout will allow a complementary $\nu$NCQE neutron multiplicity measurement, directly relevant to constraining the NCQE background in DSNB searches, competitive with the recent T2K measurement at SK-Gd. The planned Super-SANDI upgrade, deploying a large Water-based Liquid Scintillator (WbLS) volume, will further extend ANNIE's reach to hadronic final states and exclusive NC channels, and enable joint measurements with liquid argon detectors sharing the BNB beamline ahead of DUNE and Hyper-Kamiokande.

Doran, Steven [Iowa State U.]↗

Satellite Beach Energy - Restructuring the Energy Balance in Satellite Beach, Florida, by Quantifying Solar Energy Production Potential using NASA POWER Data Products and LiDAR

The City of Satellite Beach, Florida, has committed to supplying 100% of its energy use from renewable energy, primarily solar, by the year 2050. The team created a methodology for estimating rooftop solar power potential using a high-resolution Light Detection and Ranging (LiDAR) dataset and the NASA Prediction of Worldwide Energy Resources (POWER) dataset to assist Satellite Beach in reaching their solar renewable energy goals. The POWER dataset provides information on direct and diffuse solar irradiation on horizontal surfaces, surface albedo, and effects of local meteorology, such as clouds. The team integrated the solar irradiance data with the LiDAR data to model slope, aspect, and shadowing in the 7 km2 study area to find suitable roof segments for solar panel installation and estimate the solar potential of each segment. This process was supplemented by an analysis of land surface temperature and urban greenness measured through the Normalized Difference Vegetation Index (NDVI) from Landsat 8 Operational Land Imager and Thermal Infrared Sensor (OLI/TIRS) observations. These metrics serve to target areas for cooling initiatives aimed at reducing Satellite Beach’s overall energy consumption. The team found the total rooftop solar potential throughout the city to be 221,919,330KWh per year with an average annual rooftop photovoltaic, or PV, potential of 55,647KWh per building. As such, the average building could generate over five times the annual energy needs for an average household if PV panels were installed on all viable areas of its roof.

Spencer Nelson↗

TPSAS-NF1676L-16988-DND

There has been renewed interest to uniformly recalibrate historical geostationary (GEO) data records to aid in climate monitoring. GEO sensors have annual repeatable angular sampling over a given location. The view angle is fixed and the imaging schedule is usually constant through its lifetime. Given the fact that colocated GEOs are always share the same sub-satellite point and maintain their imaging schedules provides repeatable angular sampling over decades. One of the biggest challenges in transferring a reference sensor calibration using invariant desert targets to another sensor is the accuracy of the bidirectional reflectance distribution function (BDRF). A well-calibrated GEO can be used to predict the daily exoatmospheric radiance model (DERM) over a desert target for a given GMT that is valid for any GEO sensor at the same location. The advantage of this method is that a BRDF is not needed. Another challenge of invariant target calibration is the unique spectra signature of the desert. However, since most GEOs are built in batches, the spectral response functions (SRF) are very similar for most historical GEOs, the spectral band adjustment factor (SBAF) between GEO sensors is much smaller than for MODIS and GEO sensors. Since the water vapor burden over the desert is seasonal, both the TOA and desert surface can be considered invariant for a given day of the year. The reference GEO can be inter-calibrated with MODIS or VIIRS, which have onboard visible calibration using solar diffusers, using other methods, such as ray-matching or deep convective clouds. Also the next generation GEOs will have onboard visible calibration, which will increase the accuracy of this method. Three Meteosats over the Libyan desert will be used to illustrate the DERM method. The reference Meteosat will be inter-calibrated against Aqua-MODIS. The reference GEO DERM will be constructed and used to calibrate the remaining Meteosats. The DERM calibration will be validated by comparing the calibration using Aqua-MODIS ray-matching. Similarly, two GOES sensors using the Sonoran desert will also be highlighted. An uncertainty analysis will also be performed with emphasis on the SBAF, derived over the desert targets using both SCIAMACHY and Hyperion hyper-spectral radiances.

David Doelling↗

Operating envelope charts for the Langley 0.3-meter transonic cryogenic wind tunnel

To take full advantage of the unique Reynolds number capabilities of the 0.3-meter Transonic Cryogenic Tunnel (0.3-m TCT) at the NASA Langley Research Center, it was designed to accommodate test sections other than the original, octagonal, three-dimensional test section. A 20- by 60-cm two-dimensional test section was installed in 1976 and was extensively used, primarily for airfoil testing, through the fall of 1984. The tunnel was inactive during 1985 so that a new test section and improved high speed diffuser could be installed in the tunnel circuit. The new test section has solid adaptive top and bottom walls to reduce or eliminate wall interference for two-dimensional testing. The test section is 33- by 33-cm in cross section at the entrance and is 142 cm long. In the planning and running of past airfoil tests in the 0.3-m TCT, the use of operating envelope charts have proven very useful. These charts give the variation of total temperature and pressure with Mach number and Reynolds number. The operating total temperature range of the 0.3-m TCT is from about 78 K to 327 K with total pressures ranging from about 17.5 psia to 88 psia. This report presents the operating envelope charts for the 0.3-m TCT with the adaptive wall tes t section installed. They were all generated based on a 1-foot chord model. The Mach numbers vary from 0.1 to 0.95.

Rallo, R. A.↗

The Turbulent/Non-Turbulent Interface Bounding a Far-Wake

The velocity fields of a turbulent wake behind a flat plate obtained from the direct numerical simulations of Moser et al. are used to study the structure of the flow in the intermittent zone where there are, alternately, regions of fully turbulent flow and non-turbulent velocity fluctuations either side of a thin randomly moving interface. Comparisons are made with a wake that is 'forced' by amplifying initial velocity fluctuations. There is also a random temperature field T in the flow; T varies between constant values of 0.0 and 1.0 on the sides of the wake. The value of the Reynolds number based on the centreplane mean velocity defect and halfwidth b of the wake is Re approx. = 2000. It is found that the thickness of the continuous interface is about equal to 0.07b, whereas the amplitude of fluctuations of the instantaneous interface displacement y(sub I)(t) is an order of magnitude larger, being about 0.5b. This explains why the mean statistics of vorticity in the intermittent zone can be calculated in terms of the probability distribution of y(sub I) and the instantaneous discontinuity in vorticity across the interface. When plotted as functions of y - y(sub I), the conditional mean velocity (U) and temperature (T) profiles show sharp jumps Delta(U) and Delta(T) at the interface adjacent to a thick zone where (U) and (T) vary much more slowly. Statistics for the vorticity and velocity variances, available in such detail only from DNS data, show how streamwise and spanwise components of vorticity are generated by vortex stretching in the bulges of the interface. Flow fields around the interface, analyzed in terms of the local streamline pattern, confirm previous results that the advancement of the vortical interface into the irrotational flow is driven by large-scale eddy motion. It is argued that because this is an inviscid mechanism the entrainment process is not sensitive to the value of Re, and that small-scale nibbling only plays a subsidiary role. While mean Reynolds stresses decrease gradually in the intermittent zone, conditional stresses are found to decrease sharply towards zero at the interface. Using one-point turbulence models applied to either unconditional or conditional statistics for the turbulent region and then averaged, the entrainment rate E(sub b) would, if calculated exactly, be zero. But if computed with standard computational methods, E(sub b) would be non-zero because of numerical diffusion. It is concluded that the current practice in statistical models of approximating entrainment by a diffusion process is computationally arbitrary and physically incorrect. An analysis shows how E(sub b) is related to Delta(U) and the jump in shear stress at the interface, and correspondingly to Delta(T) and the heat flux.

Bisset, David K.↗

Coupled effects of conduction in the crystal and thermo-solutal convection in a rectangular inclined enclosure

To date modeling of crystal growth of optoelectronic materials using Physical Vapor Transport has been limited to the study of the fluid phase. To achieve it, the equations of coupled heat, mass and momentum transfer in the gas have to be solved. The first objective of this study is to examine the effect of heat conduction in the crystal on the fluid flow in the neighborhood of the interface. Heat transfer boundary conditions on both interfaces were modified to take into account the additional heat flux between gas and solid. It is proved that heat conduction does not affect the fluid flow. In the presence of gravity, density gradients in the fluid phase generate convection responsible for the problem of a nonplanar growth of the interface. The second objective is to study systematically under one-g the different possible flows in order to solve this problem. Depending on the parameters, a diffusive mode and three convective modes (thermal, solutal and thermo-solutal) are observed. The competition between thermal and solutal convections leads to a mathematical condition which can be used to achieve a planar growth. It is proven that, under the physical conditions chosen, this mathematical condition cannot be thermodynamically satisfied.

Mennetrier, Christophe↗

Spatiotemporal 4D Whole-cell Modeling of a Minimal Autotroph Reveals Central Carbon Metabolism Regulated Locally by Protein Megacomplexes via Post-translational Modifications under Light Disturbance

Photosynthetic microorganisms rely on multiple pathways in central carbon metabolism to adapt to fluctuating light and energy availability across diel cycles. Mechanistic insight into the regulatory dynamics of this adaptation requires integrating processes spanning disparate timescales, from rapid redox-dependent post-translational modifications (PTMs) to slower changes in protein expression and metabolic pathway usage. To address this complexity beyond genome-based inference and traditional modeling, we develop a whole-cell four-dimensional (3D + time) model of the marine cyanobacterium Prochlorococcus marinus MED4 that explicitly represents the spatial organization of enzymatic and molecular processes in central carbon metabolism under light perturbation. We employ a perturbation-based research design to experimentally generate time-series, multi-omics measurements that provide molecular descriptors and cryo-ET derived 3D segmented volumes as constraints for this dynamic 4D framework. The integration of experiments and modeling across defined light regimes enables quantitative validation of system-level responses and forecasting under distinct light disturbances. We test the hypothesis that light-dependent redox PTMs regulating the structural assembly of a protein megacomplex, the “dark complex,” modulate metabolic flux at a conserved regulatory node of the Calvin–Benson cycle (CBC) in cyanobacteria. Our model shows that subcellular spatial organization buffers rapid light-induced changes in thylakoid reaction rates, which are followed by redox-PTM-mediated sequestration or release of CBC enzymes in the dark complex, ultimately impacting carbon fixation dynamics within carboxysomes. Comparison with an equivalently parameterized well-mixed stochastic model demonstrates that post-translational regulation not only buffers transcriptional noise and diffusion-driven fluctuations but also stabilizes phenotypic outcomes, underscoring the importance of spatial heterogeneity in phenotypic robustness. This ability to probe adaptive, spatiotemporally resolved mechanisms in photosynthetic machinery and central carbon metabolism addresses a critical gap in genotype-to-phenotype inference and expands modeling and design capabilities for understudied or genetically intractable autotrophs such as P. marinus MED4.

Johnson, Connah G.↗

Challenge Problem 1: Preliminary Results of the Direct Numerical Simulation of Transient Flows

This report presents the first direct numerical simulations (DNS) of transient mixed convection in an idealized downcomer-like channel (Challenge Problem 1, Phase II). Using the GPU-accelerated NekRS solver, we modeled a sudden decay in driving pressure, mimicking loss-of-flow events, and tracked the resulting evolution of Reynolds number, boundary-layer structure, turbulence statistics, and heat-transfer metrics. Key findings include the systematic thickening and eventual asymmetry of velocity and thermal boundary layers under buoyant deceleration; minimal “memory” lag in Reynolds shear stress and TKE profiles when sampled at matching Re, yet clear shifts of peak locations toward the cooled wall; overshoots in transient eddy-viscosity and eddy-diffusivity (and corresponding sub-unity turbulent Prandtl numbers) on the cooled side; and a pronounced transient Nusselt-number enhancement driven by wall-temperature inertia and residual eddy mixing. These effects combined to offer a temporary cooling margin above steady-state predictions during reactor LOF transients. Future work will extend this work to a more complex “Case II” geometry (90° turn + lower plenum) and generate multi-Re/Pr datasets for data-driven turbulence closures.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Annual Research Briefs

This report contains the 1997 annual progress reports of the research fellows and students supported by the Center for Turbulence Research (CTR). Titles include: Invariant modeling in large-eddy simulation of turbulence; Validation of large-eddy simulation in a plain asymmetric diffuser; Progress in large-eddy simulation of trailing-edge turbulence and aeronautics; Resolution requirements in large-eddy simulations of shear flows; A general theory of discrete filtering for LES in complex geometry; On the use of discrete filters for large eddy simulation; Wall models in large eddy simulation of separated flow; Perspectives for ensemble average LES; Anisotropic grid-based formulas for subgrid-scale models; Some modeling requirements for wall models in large eddy simulation; Numerical simulation of 3D turbulent boundary layers using the V2F model; Accurate modeling of impinging jet heat transfer; Application of turbulence models to high-lift airfoils; Advances in structure-based turbulence modeling; Incorporating realistic chemistry into direct numerical simulations of turbulent non-premixed combustion; Effects of small-scale structure on turbulent mixing; Turbulent premixed combustion in the laminar flamelet and the thin reaction zone regime; Large eddy simulation of combustion instabilities in turbulent premixed burners; On the generation of vorticity at a free-surface; Active control of turbulent channel flow; A generalized framework for robust control in fluid mechanics; Combined immersed-boundary/B-spline methods for simulations of flow in complex geometries; and DNS of shock boundary-layer interaction - preliminary results for compression ramp flow.

Spinks, Debra↗

The First Fermi-LAT SNR Catalog and Cosmic Ray Implications

While supernova remnants (SNRs) are widely thought to be powerful cosmic-ray accelerators, indirect evidence comes from a small number of well-studied cases. Here we systematically determine the gamma-ray emission detected by the Fermi Large Area Telescope (LAT) from all known Galactic SNRs, disentangling them from the sea of cosmic-ray generated photons in the Galactic plane. Using LAT data we have characterized the 1-100 GeV emission in 279 regions containing SNRs, accounting for systematic uncertainties caused by source misattribution and instrumental response. We classified 30 sources as SNRs, using spatial overlap with the radio emission position. For all the remaining regions we evaluated upper limits on SNRs' emission. In the First Fermi-LAT SNR Catalog there is a study of the common characteristics of these SNRs, such as comparisons between GeV, radio and TeV quantities. We show that previously satisfactory models of SNRs' GeV emission no longer adequately describe the data. To address the question of cosmic ray (CR) origins, we also examine the SNRs' maximal CR contribution assuming the GeV emission arises solely from proton interactions. Improved breadth and quality of multiwavelength (MW) data, including distances and local densities, and more, higher resolution gamma-ray data with correspondingly improved Galactic diffuse models will strengthen this constraint.

de Palma, F.↗

The First Fermi-LAT SNR Catalog and Cosmic Ray Implications

While supernova remnants (SNRs) are widely thought to be powerful cosmicray accelerators, indirect evidence comes from a small number of well-studied cases. Here we systematically determine the gamma-ray emission detected by the ermi Large Area Telescope (LAT) from all known Galactic SNRs, disentangling them from the sea of cosmic-ray generated photons in the Galactic plane. Using LAT data we have characterized the 1-100 GeV emission in 279 regions containing SNRs, accounting for systematic uncertainties caused by source misattribution and instrumental response. We classified 30 sources as SNRs, using spatial overlap with the radio emission position. For all the remaining regions we evaluated upper limits on SNRs' emission. In the First Fermi-LAT SNR Catalog there is a study of the common characteristics of these SNRs, such as comparisons between GeV, radio and TeV quantities. We show that previously satisfactory models of SNRs' GeV emission no longer adequately describe the data. To address the question of cosmic ray (CR) origins, we also examine the SNRs' maximal CR contribution assuming the GeV emission arises solely from proton interactions. Improved breadth and quality of multiwavelength (MW) data, including distances and local densities, and more, higher resolution gamma-ray data with correspondingly improved Galactic diffuse models will strengthen this constraint.

de Palma, F.↗

An Overview of Inter-comparison Methodologies for Terra and Aqua MODIS Calibration

With increasing efforts on data fusion and long-term climate data records (CDR) using observations made by multiple sensors, on the same or different platforms, their cross-calibration and validation work has become more and more important. The uncertainties of the climate models and data records depend not only on the calibration quality of individual sensors but also on their calibration consistency. This paper provides an overview of methodologies used by the MODIS Characterization Support Team (MCST) at NASA GSFC for the inter-comparison studies of Terra and Aqua MODIS on-orbit calibration. Each MODIS was built with a set of onboard calibrators (OBC) that include a blackbody (BB), a solar diffuser (SD), and a solar diffuser stability monitor (SDSM). The BB is primarily used for the thermal emissive bands (TEB) calibration and the SD/SDSM system for the reflective solar bands (RSB) calibration. Although the-instrument design and calibration approach are nearly identical for both Terra and Aqua MODIS and they all went through an extensive and similar pre-launch calibration and characterization process, still their on-orbit calibration consistency needs to be carefully examined and validated as many science products have been generated from observations made by both instruments. Methodologies discussed in this paper include inter-comparison studies using the Moon, a third sensor, and ground targets. Our results show that Terra and Aqua reflective solar bands and thermal emissive bands have been calibrated consistently with excellent long-term stabilities. For the 11 and 12 micrometers sea surface temperature (SST) bands, the calibration difference of Terra and Aqua MODIS is less than 0.2K.

Xiong, X.↗

Numerical investigations in three-dimensional internal flows

In the present reporting period, the 3D version of the OVERFLOW code was used to solve the flow within the internal portion of the supersonic inlet. The internal portion of this inlet is bounded by an inflow plane containing the leading edge of the sidewalls, the sidewalls, the ramp and cowl surfaces and an outflow plane just downstream of the minimum geometric area of the inlet. Boundary layer bleed was used in the two-dimensional calculations discussed in the previous progress report and that same bleed was applied in the present study. For reference, this bleed corresponds to locations designated as R2 and R3 in the Mach 5 inlet model test. Using the GRIDGEN code, a three dimensional grid was generated that accounted for the viscous effects expected to occur on the sidewall, as well as those known to occur on the ramp and cowl surfaces. The internal flow grid size was 141 streamwise by 101 cross stream by 71 in the lateral direction between sidewalls. Since the flow entering the inlet was not symmetrical, the inlet was solved from sidewall to sidewall (without using a symmetry plane). In addition to the short sidewalls proposed in the Langley geometry database, a set of shorter sidewalls was also investigated in the present study and was shown to have beneficial effects with respect to the flow distortion exiting the supersonic inlet. In addition to these calculations, additional 3D solutions using the OVERFLOW code were obtained for the flow downstream of the throat of the supersonic inlet, including a terminal shock wave system produced by a backpressured subsonic diffuser.

Rose, William C.↗

Griffin: A MOOSE-based reactor physics application for multiphysics simulation of advanced nuclear reactors

Griffin is a Multiphysics Object-Oriented Simulation Environment (MOOSE) based reactor physics application for multiphysics simulations of advanced reactor designs jointly developed by Idaho National Laboratory and Argonne National Laboratory. This paper summarizes the motivation, significance, architecture, design, and features of Griffin. Griffin offers flexible and extensible features to address the challenges associated with advanced reactor designs. These features range from fundamental particle transport to specific reactor physics tasks. The features cover a wide range including on-the-fly and traditional two-step cross-section generation methods, steady-state and transient transport solvers suitable for both heterogeneous and homogeneous models, high-fidelity depletion where thousands of isotopes can be tracked and low-fidelity depletion characterized by burnup, etc. The most fundamental aspect that sets Griffin apart from other reactor analysis codes is that it is developed based on the MOOSE framework. A modular development approach is strongly enforced, with multiphysics being an essential element considered since the beginning of Griffin’s development. Griffin links various MOOSE physics modules and couples to other MOOSE-based applications and non-MOOSE-based applications for multiphyiscs simulations. Griffin includes three modules: ISOXML for preparing and managing multigroup cross sections, radiation transport for solving the neutron transport equation, and reactor analysis for user-oriented reactor physics analysis functionalities. Griffin uses various finite element methods for spatial discretization, multigroup approximation for energy discretization and discrete ordinates method, spherical harmonics expansion method, and diffusion approximation for streaming direction discretization to solve the neutron transport equation. Griffin’s flexibility is evidenced through Griffin’s various applications to fast reactor, high-temperature reactor, pebble bed reactor, molten salt reactor, and microreactor designs. Griffin development follows the software quality assurance procedure for MOOSE-based applications and with software requirements consistent with the ASME NQA-1 standard. Griffin has been adopted into the reactor analysis system for the U.S. NRC and is in use at U.S. companies, universities and national laboratories.

97 MATHEMATICS AND COMPUTING↗