Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “channel estimation”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Cross-Layered Cyber-Physical Power System State Estimation towards a Secure Grid Operation

In the Smart Grid paradigm, this critical infrastructure operation is increasingly exposed to cyber-threats due to the increased dependency on communication networks. An adversary can launch an attack on a power grid operation through False Data Injection into system measurements and/or through attacks on the communication network, such as flooding the communication channels with unnecessary data or intercepting messages. A cross-layered strategy that combines power grid data, communication grid monitoring and Machine Learning based processing is a promising solution for detecting cyberthreats. In this paper, an implementation of an integrated solution of a cross-layer framework is presented. The advantage of such a framework is the augmentation of valuable data that enhances the detection of anomalies in the operation of power grid. IEEE 118-bus system is built in Simulink to provide a power grid testing environment and communication network data is emulated using SimComponents. The performance of the framework is investigated under various FDI and communication attacks.

cyber security, network security, cyber-physical s↗
DOE OSTI2023

Urban Flood Modeling: Uncertainty Quantification and Physics‐Informed Gaussian Processes Regression Forecasting

Abstract Estimating uncertainty in flood model predictions is important for many applications, including risk assessment and flood forecasting. We focus on uncertainty in physics‐based urban flooding models. We consider the effects of the model's complexity and uncertainty in key input parameters. The effect of rainfall intensity on the uncertainty in water depth predictions is also studied. As a test study, we choose the Interconnected Channel and Pond Routing (ICPR) model of a part of the city of Minneapolis. The uncertainty in the ICPR model's predictions of the floodwater depth is quantified in terms of the ensemble variance using the multilevel Monte Carlo (MC) simulation method. Our results show that uncertainties in the studied domain are highly localized. Model simplifications, such as disregarding the groundwater flow, lead to overly confident predictions, that is, predictions that are both less accurate and uncertain than those of the more complex model. We find that for the same number of uncertain parameters, increasing the model resolution reduces uncertainty in the model predictions (and increases the MC method's computational cost). We employ the multilevel MC method to reduce the cost of estimating uncertainty in a high‐resolution ICPR model. Finally, we use the ensemble estimates of the mean and covariance of the flood depth for real‐time flood depth forecasting using the physics‐informed Gaussian process regression method. We show that even with few measurements, the proposed framework results in a more accurate forecast than that provided by the mean prediction of the ICPR model.

Kohanpur, Amir H.↗
DOE OSTI2026

Quantifying the Impact of Metal Population Distribution in MFI-Supported Mo Catalysts for Methane Dehydroaromatization

Precise evaluation of intrinsic kinetic behavior in Mo/MFI catalysts for methane dehydroaromatization (MDA) is confounded by variations in Mo dispersion and speciation, which are influenced by metal loading and zeolite acidity. This work advances the utility of H 2 -temperature programmed reduction (H 2 -TPR) for characterizing Mo/MFI catalysts by enabling quantitative comparison of MoO x populations distinguished by reduction behavior and linked to initial catalytic performance. A comprehensive H 2 reduction pathway is established through systematic H 2 - TPR studies varying catalyst composition (1–10 wt % Mo, Si/Al = 15, 40, ∞), supplemented by UV-Raman spectroscopy, X-ray powder diffraction, N 2 physisorption, NH 3 -TPD, and advanced spectroscopic analysis (in situ XAS with principal component analysis/multivariate curve resolution—alternating least squares). Two low-temperature H 2 -TPR regions capture distinct Mo populations undergoing initial Mo(VI)→Mo(IV) reduction: a lower-temperature population (Mo-RI) associated primarily with highly dispersed, anchored MoO x species expected to predominantly reside within MFI channels, and a higher-temperature population (Mo-RII) corresponding to a broader set of MoO x species that becomes increasingly bulk-like/extrazeolitic at higher Mo loading. Quantification of these populations provides practical, kinetically relevant descriptors for comparing initial MDA rates across catalysts with varying Mo loading, Si/Al ratio, and MoO x heterogeneity. Application of lower-temperature Mo-RI estimates to kinetic measurements reveals a minimum threshold of ∼0.12 × 10 –3 mol Mo-RI/g cat , above which initial forward benzene rates normalized to this population converge despite differences in metal loading and zeolite Brønsted acidity. This threshold coincides with a transition toward a common C 2 -mediated benzene-forming regime, as indicated by approach-to-equilibrium analysis of methane-to-ethane, ethane-to-ethylene, and ethylene-to-benzene reaction steps. Above this threshold, normalized initial benzene rates are nearly invariant with increasing Mo-RII/Mo-RI population ratio, indicating that excess Mo-RII populations, including bulk-like/extrazeolitic MoO x domains present at higher loading, do not measurably suppress benzene formation associated with anchored and mostly channel-confined Mo population under the initial-rate conditions examined.

Mo/MFI↗
DOE OSTI2024

Weak Gravitational Lensing around Low Surface Brightness Galaxies in the DES Year 3 Data

We present galaxy-galaxy lensing measurements using a sample of low surface brightness galaxies (LSBGs) drawn from the Dark Energy Survey Year 3 (Y3) data as lenses. LSBGs are diffuse galaxies with a surface brightness dimmer than the ambient night sky. These dark-matter-dominated objects are intriguing due to potentially unusual formation channels that lead to their diffuse stellar component. Given the faintness of LSBGs, using standard observational techniques to characterize their total masses proves challenging. Weak gravitational lensing, which is less sensitive to the stellar component of galaxies, could be a promising avenue to estimate the masses of LSBGs. Our LSBG sample consists of 23,790 galaxies separated into red and blue color types at g - i ≥ 0.60 and g - i < 0.60 , respectively. Combined with the DES Y3 shear catalog, we measure the tangential shear around these LSBGs and find signal-to-noise ratios of 6.67 for the red sample, 2.17 for the blue sample, and 5.30 for the full sample. We use the clustering redshifts method to obtain redshift distributions for the red and blue LSBG samples. Assuming all red LSBGs are satellites, we fit a simple model to the measurements and estimate the host halo mass of these LSBGs to be log(M host /M ⊙ ) = $12.98^{+0.10}_{-0.11}$. We place a 95% upper bound on the subhalo mass at log(M sub /M ⊙ ) < 11.51. By contrast, we assume the blue LSBGs are centrals, and place a 95% upper bound on the halo mass at log(M host /M ⊙ ) < 11.84. We find that the stellar-to-halo mass ratio of the LSBG samples is consistent with that of the general galaxy population. This work illustrates the viability of using weak gravitational lensing to constrain the halo masses of LSBGs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Properties and Astrophysical Implications of the 150 M ⊙ Binary Black Hole Merger GW190521

The gravitational-wave signal GW190521 is consistent with a binary black hole (BBH) merger source at redshift 0.8 with unusually high component masses, 85-14+21 M o˙ and 66-18+17 M ⊙ , compared to previously reported events, and shows mild evidence for spin-induced orbital precession. The primary falls in the mass gap predicted by (pulsational) pair-instability supernova theory, in the approximate range 65-120 M ⊙ . The probability that at least one of the black holes in GW190521 is in that range is 99.0%. The final mass of the merger (142-16+28 M ⊙ ) classifies it as an intermediate-mass black hole. Under the assumption of a quasi-circular BBH coalescence, we detail the physical properties of GW190521's source binary and its post-merger remnant, including component masses and spin vectors. Three different waveform models, as well as direct comparison to numerical solutions of general relativity, yield consistent estimates of these properties. Tests of strong-field general relativity targeting the merger-ringdown stages of the coalescence indicate consistency of the observed signal with theoretical predictions. We estimate the merger rate of similar systems to be 0.13-0.11+0.30 Gpc -3 yr -1 . We discuss the astrophysical implications of GW190521 for stellar collapse and for the possible formation of black holes in the pair-instability mass gap through various channels: via (multiple) stellar coalescences, or via hierarchical mergers of lower-mass black holes in star clusters or in active galactic nuclei. We find it to be unlikely that GW190521 is a strongly lensed signal of a lower-mass black hole binary merger. We also discuss more exotic possible sources for GW190521, including a highly eccentric black hole binary, or a primordial black hole binary.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Toward AlGaN channel HEMTs on AlN: Polarization-induced 2DEGs in AlN/AlGaN/AlN heterostructures

Due to its high breakdown electric field, the ultra-wide bandgap semiconductor AlGaN has garnered much attention recently as a promising channel material for next-generation high electron mobility transistors (HEMTs). A comprehensive experimental study of the effects of Al composition x on the transport and structural properties is lacking. We report the charge control and transport properties of polarization-induced 2D electron gases (2DEGs) in strained AlGaN quantum well channels in molecular-beam-epitaxy-grown AlN/Al x Ga 1-x N/AlN double heterostructures by systematically varying the Al content from x = 0 (GaN) to x = 0.74, spanning energy bandgaps of the conducting HEMT channels from 3.49 to 4.9 eV measured by photoluminescence. This results in a tunable 2DEG density from 0 to 3.7 × 10 13 cm 2 . The room temperature mobilities of x ≥ 0.25 AlGaN channel HEMTs were limited by alloy disorder scattering to below 50 cm 2 /(V.s) for these 2DEG densities, leaving ample room for further heterostructure design improvements to boost mobilities. A characteristic alloy fluctuation energy of ≥1.8 eV for electron scattering in AlGaN alloy is estimated based on the temperature dependent electron transport experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2026

BSEC Air Quality: Stationary Black and Brown Carbon Properties and Mass Concentration

Data is 1-min absorption/attenuation of aerosol collected on filter and sampled at seven wavelengths (950, 880, 660, 590, 520, 470, and 370 nm) with a Magee AE-33 dual spot aethalometer. The wavelength value for black carbon, as recommended by the manufacturer, is the 880 nm channel. The derived mass concentrations, using the factory default mass absorption cross-sections, are reported for all seven wavelengths. The absorption at all seven wavelengths are reported by dividing the derived mass concentrations by the factory default mass absorption cross-sections. The estimated mass concentration contribution from biomass burning and fossil fuel black carbon, using factory default values, are reported. The estimated brown carbon mass concentration is derived by finding the difference between the recommended black carbon wavelength (880 nm) and the ultra violet wavelength (370 nm), assuming that brown carbon absorption is strongest at 370 nm and minimal at 880 nm. Finally, the estimated absorption Angstrom exponent for black and brown carbon is derived. For black carbon, the slope between each minute log(AbsXXX), where XXX is the wavelength used versus log(XXX) is used, and the wavelengths are 590, 660, 880, and 950 nm. For brown carbon, the slope between each minute log(AbsXXX), where XXX is the wavelength used versus log(XXX) is used, and the wavelengths are 370, 470, and 520 nm.

aerosol↗
DOE OSTI2023

On the convergence of statistics in simulations of stationary incompressible turbulent flows

When reporting statistics from simulations of statistically stationary chaotic phenomenon, it is important to verify that the simulations are time-converged. This condition is connected with the statistical error or number of digits with which statistics can be reliably reported. In this work we consider numerical experiments of low Reynolds number incompressible homogeneous and isotropic turbulence as a model problem to investigate statistical convergence over finite simulation times. Specifically, we investigate the time integration requirements that allow meaningful reporting of the statistical error associated with finiteness of the temporal domain. We address two key questions: (1) How long should a simulation be performed in terms of large eddy time, and (2) How should the simulation time be divided among temporal windows over which a quantity of interest is estimated so that its statistical error could be reliably reported? We find that reliable reporting of statistical errors requires simulations on the order of 10 4 large eddy times, which is orders of magnitude longer than typically performed. Additionally, data post-processing should employ windows of at least ten times the large eddy time scale, with the most robust computation of statistical error of the mean requiring window sizes of an additional factor of ten. For practical simulations, we demonstrate that it is possible to estimate the statistical error within a factor of two under a less stringent condition in which a minimum of four windows with size at least ten large eddy times are used. In conclusion, our observations for homogeneous isotropic turbulence are also shown to hold in turbulent channel flow.

42 ENGINEERING↗
DOE OSTI2023

Muon capture on the deuteron in chiral effective field theory

Here, we consider the capture of a muon on a deuteron. An uncertainty analysis of the dominant channels is important for a careful analysis of forthcoming experimental data. We quantify the theoretical uncertainties of chiral effective-field-theory predictions of the muon-deuteron capture rate from the relevant neutron-neutron partial wave channels in the final state. We study the dependence on the cutoff used to regularize the interactions, low-energy constants calibrated using different fitting data and strategies, and truncation of the effective-field-theory expansion of the currents. Combining these approaches gives as an estimate of $Γ^{1/2}_{μd}$ = 399.1 ± 7.6 ± 4.4 s –1 for capture from the atomic doublet state, and $Γ^{3/2}_{μd}$= 12.31 ± 0.47 ± 0.04 s –1 for capture from the quartet state and the first and second uncertainties given here are due to the effective field theory truncation error and the uncertainty in the axial radius, respectively.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

A workflow for automatic generation and efficient refinement of individual pressure-dependent networks

Manual chemical kinetic model construction often requires modelists to at least implicitly guess all possible decay pathways and their relative fluxes for all important species. We present a workflow and associated tool for automatic generation and re finement of pressure-dependent networks that should enable modelists to efficiently and comprehensively identify decay pathways and estimate respective parameters for chemical species. This tool combines the capabilities of the Reaction Mechanism Generator (RMG) software to generate possible reaction paths, determine thermochemistry, approximate rate coefficients and estimate frequencies with the capabilities of the Arkane software to use quantum chemical parameters to compute thermochemistry and pressure dependent rates. A flux-based algorithm is used to decide which isomers to add to the network. Isomers added to the network are reacted to form other channels. When enough of the flux is accounted for by the isomers and bimolecular product channels, the generation process terminates to yield a comprehensive network. Network sensitivity analysis is applied to the network to identify the most important wells and barriers that can be refined using quantum chemistry calculations. Iterative refinement can then be used to achieve accuracy. The comprehensive network can be easily reduced using energy and flux-based algorithms to the most important channels and isomers in the network.

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

Lithium indium diselenide — An advanced material for neutron detection

Here, this paper describes the synthesis, crystal growth, detector fabrication, radiation hardening studies, MCNP modeling, and characterization of lithium indium diselenide or LiInSe 2 . This newly-developed room temperature thermal neutron detector has semiconducting and scintillating properties and it is suitable for neutron detection application. LiInSe 2 was synthesized starting from elemental Li, In, Se in two steps due to high reactivity of Li. A single crystal of LiInSe 2 was grown using the Vertical Bridgman method. The room temperature band gap was found to be 2.8 eV using optical absorption measurements. Bulk resistivity was measured at ∼5 × 10 11 Ωcm. Photoconductivity measurements of LiInSe 2 wafers identified a peak in the photocurrent around 445 nm. Nuclear radiation detectors were fabricated from single crystal wafer and the responses to alpha particles at various biases were measured. The mobility-lifetime product was estimated. Gamma irradiation studies were performed with calculated absorbed doses ranging from 0.2126 to 21,262 Gy. The characterization of the two wafers for their scintillator performance was conducted after each irradiation. The gamma irradiation produced a reduction of the light yield that translated to a lower channel number for the centroid of alpha detection spectra. It also showed a considerable reduction of the decay time after the first irradiation. These are the first studies on gamma radiation hardening with this material.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Evidence for pressure induced unconventional quantum criticality in the coupled spin ladder antiferromagnet C 9 H 18 N 2 CuBr 4

Quantum phase transitions in quantum matter occur at zero temperature between distinct ground states by tuning a nonthermal control parameter. Often, they can be accurately described within the Landau theory of phase transitions, similarly to conventional thermal phase transitions. However, this picture can break down under certain circumstances. Here, we present a comprehensive study of the effect of hydrostatic pressure on the magnetic structure and spin dynamics of the spin-1/2 ladder compound C 9 H 18 N 2 CuBr 4 . Single-crystal heat capacity and neutron diffraction measurements reveal that the Néel-ordered phase breaks down beyond a critical pressure of Pc~1.0 GPa through a continuous quantum phase transition. Estimates of the critical exponents suggest that this transition may fall outside the traditional Landau paradigm. The inelastic neutron scattering spectra at 1.3 GPa are characterized by two well-separated gapped modes, including one continuum-like and another resolution-limited excitation in distinct scattering channels, which further indicates an exotic quantum-disordered phase above Pc.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Role of low-lying resonances for the 10 Be(p,α) 7 Li reaction rate and implications for the formation of the Solar System

Evidence for the presence of short-lived radioactive isotopes when the Solar System formed is preserved in meteorites, providing insights into the conditions at the birth of our Sun. A low-mass core-collapse supernova had been postulated as a candidate for the origin of 10Be, reinforcing the idea that a supernova triggered the formation of the Solar System. We present a detailed study of the production of 10Be by the ν process in supernovae, which is very sensitive to the reaction rate of the major destruction channel, 10 Be(p,α) 7 Li. With data from recent nuclear experiments that show the presence of a resonant state in 11B at ≈193 keV, we derive new values for the 1 0 Be(p,α) 7 Li reaction rate, which are significantly higher than previous estimates. We show that, with the new 10 Be(p,α) 7 Li reaction rate, a low-mass CCSN is unlikely to produce enough 10 Be to explain the observed 10 Be/ 9 Be ratio in meteorites, even for a wide range of neutrino spectra considered in our models.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Poles and poltergeists in e + e − → D D ¯ data

A recent report of e + e − → D D ¯ events by the BESIII collaboration suggests the presence of a structure R at 3900 MeV. We argue that this structure, called G ( 3900 ) in the past, is not in fact due to a new c c ¯ resonance but rather naturally emerges due to a combination of interference between nearby resonances and the opening of the D * D ¯ channel. We further find that the appearance of this structure does not require suppression because of a radial node in the ψ ( 4040 ) wave function, although a node improves fit quality. The measured e + e − coupling of ψ ( 4040 ) is found to be substantially smaller than previously estimated. In addition, we report new corrections to the measured cross section σ ( e + e − → D D ¯ ) at energies near ψ ( 3770 ) . Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2021

A Novel Machine Learning Algorithm for Planetary Boundary Layer Height Estimation Using AERI Measurement Data

Accurately determining the height of the planetary boundary layer (PBL) is important since it can affect the climate, weather, and air quality. Ground-based infrared hyperspectral remote sensing is an effective way to obtain this parameter. Compared with radiosonde measurements, its temporal resolution is much higher. In this study, a method to retrieve the PBL height (PBLH) from the ground-based infrared hyperspectral radiance data is proposed based on machine learning. In this method, the channels that are sensitive to temperature and humidity profiles are selected as the feature vectors, and the PBLHs derived from radiosonde are taken as the true values. The support vector machine (SVM) is applied to train and test the data set, and the parameters are optimized in the process. The data set collected at the Atmospheric Radiation Measurement (ARM) program Southern Great Plains (SGP) from 2012 to 2015 is analyzed. The instruments used in this letter include Atmospheric Emitted Radiance Interferometer (AERI), Vaisala CL31 ceilometer, and radiosonde. It shows that the root mean square error (RMSE) between the PBLHs calculated by the proposed method using AERI data and those from radiosonde data can be within 370 m, and the square correlation coefficient (SCC) is greater than 0.7. Compared with the PBLHs derived from the ceilometer, it can be found that the new method is more stable and less affected by clouds.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Quantifying Thermal Output of Energetic Materials (LDRD Final Report)

We present the results of an LDRD project, funded by the Nuclear Deterrence IA, to develop capabilities for quantitative assessment of pyrotechnic thermal output. The thermal battery igniter is used as an exemplar system. Experimental methodologies for thermal output evaluation are demonstrated here, which can help designers and engineers better specify pyrotechnic components , provide thermal output guidelines for new formulations, and generate new metrics for assessing component performance and margin given a known failure condition. A heat-transfer analysis confirms that the dominant mode of energy transfer from the pyrotechnic output plume to the heat pellet is conduction via deposition of hot titanium particles. A simple lumped-parameter model of titanium particle heat transfer and a detailed multi-phase model of deposition heat transfer are discussed. Pyrotechnic function, as defined by "go/no-go" standoff testing of a heat pellet, is correlated with experimentally measured igniter plume temperature, titanium metal particle temperature, and energy deposition. Three high-speed thermal diagnostics were developed for this task. A three-color imaging pyrometer, acquiring 100k images per second on three color channels, is deployed for measurement of titanium particle temperatures. Complimentary measurements of the overall igniter plume emission ("color") temperature were conducted using a transmission-grating spectrograph in line-imaging mode. Heat flux and energy deposition to a cold wall at the heat-pellet location were estimated using an eroding thermocouple probe, with a frequency response of ~5 kHz. Ultimate "go/no-go" function in the igniter/heat-pellet system was correlated with quantitative thermal metrics, in particular surface energy deposition and plume color temperature. Titanium metal-particle and plume color temperatures both experience an upper bound approximated by the 3245-K boiling point of TiO 2 . Average metal-particle temperatures remained nearly constant for all standoff distances at T = 2850 K, ± 300 K, while plume color temperature and heat flux decay with standoff—suggesting that heat-pellet failure results from a drop in metal-particle flux and not particle temperature. At 50% likelihood of heat-pellet failure, peak time-resolved plume color temperatures drop well below TiO 2 boiling to ~2000 - 2200 K, near the TiO 2 melting point. Estimates of peak heat flux decline from up to 1 GW/m 2 for near-field standoffs to below 320 MW/m 2 at 50% failure likelihood.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Investigations of Relative Stability and Distribution of Mo Species in Pore Space of Mo/ZSM-5 Systems

Reduced forms of molybdenum (Mo) carbides anchoring in pore space of zeolites are considered as activated catalysts for the methane dehydroaromatization reaction. However, the atomic structure of the active sites and their stability under reactive conditions are still not clearly understood. Herein, we employ a synergistic theoretical and experimental approach to investigate this research gap. The Raman data shows the existence of the distributed binuclear and mononuclear Mo oxides in the pores that are corroborated by the hydrogen temperature-programmed reduction experiments. The binding free energies of mono- and bi-nuclear Mo carbides in straight and sinusoidal channels and their intersections as a function of temperature and anchoring sites have been determined. The catalysts with Mo coordinated at two anchoring sites are more stable compared to those located at the rings with one Al substitution. The activation energies of the catalysts migration to adjacent rings were estimated as a function of temperature and kinetics parameters of such migrations were determined. The migration is facilitated if there is one anchoring site shared between two rings. This process can contribute to an agglomeration mechanism leading to loss of catalytic activity.

Myshakin, Evgeniy↗
DOE OSTI2023

Development of experimental and computational frameworks to predict subcooled flow boiling in the LANL Isotope Production Facility

Cooling is crucial to maintain the integrity of target systems in isotope production facilities. At Los Alamos National Laboratory (LANL)’s Isotope Production Facility (IPF), multiple encapsulated targets are stacked and irradiated in tandem with a 100 MeV, ~250μA proton beam. To facilitate effective heat removal, these stacked targets are separated and cooled via a series of water channels. At these beam currents, this high-energy proton beam heats the target system, likely initiating subcooled flow boiling in the cooling channels. However, in-beam monitoring of the IPF target system is not possible due to the extreme radiation environment, and the necessarily significant shielding. To better understand high-power target performance, we developed ex-situ experimental and computational frameworks to predict the behavior of subcooled flow boiling at IPF. Subcooled flow boiling experiments on Inconel 625 samples under IPF conditions (2 bar pressure, 10 GPM flow rate (i.e., 2249 kg/m 2 /s), 85 K subcooling) revealed that IPF's average operating power is at the early stage of boiling with a heat transfer coefficient of 48,000 W/m 2 /s. The proposed modeling framework enables us to predict a complete boiling curve, i.e., single-phase heat transfer, onset of nucleate boiling, two-phase heat transfer, and critical heat flux (CHF), with specification of input boiling parameters up to intermediate heat flux levels. The estimated CHF under IPF conditions is 5.2 MW/m 2 . Experimental data under reduced conditions (2 bar pressure, 1.5 GPM flow rate (i.e., 337 kg/m 2 /s), 45 K subcooling) served as validation cases for the computational modeling. This computational model can be further extended to more complicated systems replicating the real IPF configuration, for instance, to study void distribution as a function of the incident proton beam profile and coolant velocity profile of multiple cooling channels. Finally, the proposed experimental and computational frameworks provide a means to better understand cooling systems in the isotope production facilities at different accelerators, where in-beam monitoring of the cooling process is not available.

07 ISOTOPE AND RADIATION SOURCES↗