Simulating Photoelectric Radiation Transport in Cylindrical Gas Filled Photoemission Driven Cavities
Explore the source record for details and available documents.
SEARCH · Search NASA
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.
Explore the source record for details and available documents.
Fast and accurate approximations of gas cooling and heating functions are needed for hydrodynamic galaxy simulations. We use machine learning to analyze atomic gas cooling and heating functions in the presence of a generalized incident local radiation field computed by Cloudy. We characterize the radiation field through binned radiation field intensities instead of the photoionization rates used in our previous work. We find a set of 6 energy bins whose intensities exhibit relatively low correlation. We use these bins as features to train machine learning models to predict Cloudy cooling and heating functions at fixed metallicity. We compare the relative SHapley Additive exPlanation (SHAP) value importance of the features. From the SHAP analysis, we identify a feature subset of 3 energy bins (0.5-1, 1-4, and 13-16Ry) with the largest importance and train additional models on this subset. We compare the mean squared errors and distribution of errors on both the entire training data table and a randomly selected 20% test set withheld from model training. The machine learning models trained with 3 and 6 bins, as well as 3 and 4 photoionization rates, have comparable accuracy everywhere, with errors ≳10 times smaller than for the interpolation table of Gnedin and Hollon (2012). We conclude that 3 energy bins (or 3 analogous photoionization rates: molecular hydrogen photodissociation, neutral hydrogen HI, and fully ionized carbon CVI) are sufficient to characterize the dependence of the gas cooling and heating functions on our assumed incident radiation field model.
Fast and accurate approximations of gas cooling and heating functions are needed for hydrodynamic galaxy simulations. We use machine learning to analyze atomic gas cooling and heating functions computed by Cloudy in the presence of a generalized incident local radiation field. We characterize the radiation field through binned radiation field intensities instead of the photoionization rates used in our previous work. We find a set of 6 energy bins whose intensities exhibit relatively low correlation. We use these bins as features to train machine learning models to predict Cloudy cooling and heating functions at fixed metallicity. We compare the relative SHapley Additive exPlanation (SHAP) value importance of the features. From the SHAP analysis, we identify a feature subset of 3 energy bins ($0.5-1, 1-4$, and $13-16 \, \mathrm{Ry}$) with the largest importance and train additional models on this subset. We compare the mean squared errors and distribution of errors on both the entire training data table and a randomly selected 20% test set withheld from model training. The machine learning models trained with 3 and 6 bins, as well as 3 and 4 photoionization rates, have comparable accuracy everywhere, with errors $\gtrsim 10$ times smaller than for the interpolation table of Gnedin and Hollon (2012). We conclude that 3 energy bins (or 3 analogous photoionization rates: molecular hydrogen photodissociation, neutral hydrogen HI, and fully ionized carbon CVI) are sufficient to characterize the dependence of the gas cooling and heating functions on our assumed incident radiation field model.
Accurate and efficient interatomic potentials are essential for molecular dynamics (MD) simulations of radiation damage, gas diffusion, and phase stability in complex ceramics such as LiAlO 2 , especially under extreme conditions relevant to tritium production. Here, we evaluate the performance of six machine-learned interatomic potentials (MLIPs), moment tensor potential (MTP), Gaussian approximation potential, deep potential (DeePMD), atomic cluster expansion (ACE), message-passing ACE (multilayer atomic cluster expansion (MACE) pretrained) and MACE (trained from-scratch), all trained on the same density functional theory dataset with inclusion of tritium. The MLIPs are benchmarked against traditional Buckingham and ReaxFF potentials in terms of energy accuracy, density predictions, thermal equilibration behavior, threshold displacement energy (E d ), tritium diffusivity, and computational cost. Among the models, MTP shows the best overall balance between efficiency and accuracy, with low force and energy errors and realistic E d values for Li and Al. The ACE and MACE (pretrained and trained from scratch) models exhibit high E d (>200 eV) and unphysical pair interactions. DeePMD underestimates Ed due to overly repulsive behavior even at equilibrium distances. All models over-estimate tritium diffusion but the pretrained MACE model behaves well during tritium-diffusion simulations up to 500 K, maintaining diffusivities in the physically consistent 10 −11 m 2 /s range. Finally, we quantify the computational cost of each potential in large-scale atomic/molecular massively parallel simulator, finding that only MTP is more efficient than traditional empirical potentials, while others are significantly more expensive. These findings explain the trade-offs between accuracy and computational cost in MLIP development and provide essential guidance for use in high-throughput radiation damage and gas diffusion simulations in nuclear ceramics.
Synchrotron radiation from hot gas near a black hole results in a polarized image. The image polarization is determined by effects including the orientation of the magnetic field in the emitting region, relativistic motion of the gas, strong gravitational lensing by the black hole, and parallel transport in the curved spacetime. We explore these effects using a simple model of an axisymmetric, equatorial accretion disk around a Schwarzschild black hole. By using an approximate expression for the null geodesics derived by Beloborodov and conservation of the Walker–Penrose constant, we provide analytic estimates for the image polarization. We test this model using currently favored general relativistic magnetohydrodynamic simulations of M87*, using ring parameters given by the simulations. For a subset of these with modest Faraday effects, we show that the ring model broadly reproduces the polarimetric image morphology. Our model also predicts the polarization evolution for compact flaring regions, such as those observed from Sgr A* with GRAVITY. With suitably chosen parameters, our simple model can reproduce the EVPA pattern and relative polarized intensity in Event Horizon Telescope images of M87*. Under the physically motivated assumption that the magnetic field trails the fluid velocity, this comparison is consistent with the clockwise rotation inferred from total intensity images.
ABSTRACT Gas cooling and heating functions play a crucial role in galaxy formation. But, it is computationally expensive to exactly compute these functions in the presence of an incident radiation field. These computations can be greatly sped up by using interpolation tables of pre-computed values, at the expense of making significant and sometimes even unjustified approximations. Here, we explore the capacity of machine learning to approximate cooling and heating functions with a generalized radiation field. Specifically, we use the machine learning algorithm XGBoost to predict cooling and heating functions calculated with the photoionization code cloudy at fixed metallicity, using different combinations of photoionization rates as features. We perform a constrained quadratic fit in metallicity to enable a fair comparison with traditional interpolation methods at arbitrary metallicity. We consider the relative importance of various photoionization rates through both a principal component analysis (PCA) and calculation of SHapley Additive exPlanation (shap) values for our XGBoost models. We use feature importance information to select different subsets of rates to use in model training. Our XGBoost models outperform a traditional interpolation approach at each fixed metallicity, regardless of feature selection. At arbitrary metallicity, we are able to reduce the frequency of the largest cooling and heating function errors compared to an interpolation table. We find that the primary bottleneck to increasing accuracy lies in accurately capturing the metallicity dependence. This study demonstrates the potential of machine learning methods such as XGBoost to capture the non-linear behaviour of cooling and heating functions.
Miniaturized photoionization detectors (PIDs) are used in conjunction with gas chromatography systems to detect volatile compounds in gases by collecting the current from the photoionized gas analytes. PIDs should be inexpensive and compatible with a wide range of analyte species. One such PID is based on the formation of a He plasma in a dielectric barrier discharge (DBD), which generates vacuum UV (VUV) photons from excited states of He to photoionize gas analytes. There are several design parameters that can be leveraged to increase the ionizing photon flux to gas analytes to increase the sensitivity of the PID. To that end, the methods to maximize the photon flux from a pulsed He plasma in a DBD-PID were investigated using a two-dimensional plasma hydrodynamics model. The ionizing photon flux originated from the resonance states of helium, He(3P) and He(2 1 P), and from the dimer excimer He 2 *. While the photon flux from the resonant states was modulated over the voltage pulse, the photon flux from He 2 * persisted long after the voltage pulse passed. Several geometrical optimizations were investigated, such as using an array of pointed electrodes. However, increasing the capacitance of the dielectric enclosing the plasma chamber had the largest effect on increasing the VUV photon fluence to gas analytes.
Uranium mononitride (UN) is one of the ceramic nuclear fuel alternatives to oxide fuels considered for light water reactor and advanced reactor designs, as it presents significant advantages such as high uranium density (better economics) and high thermal conductivity and melting point (increased safety). Self- and fission gas diffusivities need to be better understood, given that they influence key fuel performance phenomena like swelling and fission gas release. Recently, radiation enhanced diffusivity was investigated in UN by means of cluster dynamics simulations relying on empirical potential-based parameterizations, the reliability of which highly depends on the interatomic potential accuracy. Here, in this work, we refine this approach by determining, using ab-initio calculations, the properties of defect clusters containing vacancies, self-interstitials and Xe impurities. We also consider larger clusters than previous studies. The obtained dataset (formation enthalpies, entropies, and migration barriers) is used to parameterize a cluster dynamics model of mobile clusters, and to calculate the defect cluster concentrations under irradiation. This gives us access to the radiation enhanced self- and fission gas diffusivities. Although the resulting diffusivities are close to the values reported in the literature, we find important qualitative differences in the diffusion mechanisms. Capturing the correct mechanisms is crucial to properly describe the chemistry and fission rate dependence of the model.
Published data by Task 2 (Oxyhydroxide Layer Radiolytic Gas Generation Resolution) for Technical Considerations and Challenges for Extended (>50 yrs) Dry Storage of Aluminum Cladded Spent Nuclear Fuel (ASNF), demonstrated that radiolytic molecular hydrogen (H2) production from gamma irradiated aluminum alloy 1100 (Al-1100) coupons exhibited significant dependence on backfill gaseous environment conditions: air, due to the presence of oxygen, completely inhibited H2 production; nitrogen promoted H2 production; and argon yielded more H2 than nitrogen environments. The concern here is that helium has been proposed as the backfill gas for extended storage and is more inert than argon, which may translate into significantly more H2 production than current Task 3 argon-based models predict. However, the measurement of H2 in helium media was not possible using the previously establish gas chromatography (GC) flame ionization detector, due to similarities in thermal conductivity. This milestone was initiated to evaluate an alternative GC approach that employed a mercuric oxide (HgO) reduction gas detector (RGD). Using the HgO RGD setup, H2 was successfully calibrated in the presence of a helium carrier gas, and then subsequently measured in control and irradiated crush-tube sample vials, consistent with previous sample measurements. Overall, the HgO RGD approach was found to be sufficient for future helium environment irradiations.
Intergovernmental Panel on Climate Change (IPCC) assessments are the trusted source of scientific evidence for climate negotiations taking place under the United Nations Framework Convention on Climate Change (UNFCCC). Evidence-based decision-making needs to be informed by up-to-date and timely information on key indicators of the state of the climate system and of the human influence on the global climate system. However, successive IPCC reports are published at intervals of 5–10 years, creating potential for an information gap between report cycles. We follow methods as close as possible to those used in the IPCC Sixth Assessment Report (AR6) Working Group One (WGI) report. We compile monitoring datasets to produce estimates for key climate indicators related to forcing of the climate system: emissions of greenhouse gases and short-lived climate forcers, greenhouse gas concentrations, radiative forcing, the Earth's energy imbalance, surface temperature changes, warming attributed to human activities, the remaining carbon budget, and estimates of global temperature extremes. The purpose of this effort, grounded in an open-data, open-science approach, is to make annually updated reliable global climate indicators available in the public domain. As they are traceable to IPCC report methods, they can be trusted by all parties involved in UNFCCC negotiations and help convey wider understanding of the latest knowledge of the climate system and its direction of travel. The indicators show that, for the 2014–2023 decade average, observed warming was 1.19 [1.06 to 1.30] °C, of which 1.19 [1.0 to 1.4] °C was human-induced. For the single-year average, human-induced warming reached 1.31 [1.1 to 1.7] °C in 2023 relative to 1850–1900. The best estimate is below the 2023-observed warming record of 1.43 [1.32 to 1.53] °C, indicating a substantial contribution of internal variability in the 2023 record. Human-induced warming has been increasing at a rate that is unprecedented in the instrumental record, reaching 0.26 [0.2–0.4] °C per decade over 2014–2023. This high rate of warming is caused by a combination of net greenhouse gas emissions being at a persistent high of 53±5.4 Gt CO 2 e yr -1 over the last decade, as well as reductions in the strength of aerosol cooling. Despite this, there is evidence that the rate of increase in CO 2 emissions over the last decade has slowed compared to the 2000s, and depending on societal choices, a continued series of these annual updates over the critical 2020s decade could track a change of direction for some of the indicators presented here.
Classical molecular dynamics simulations show that production, accumulation, and evolution of defects in monolayer phosphorene can be precisely controlled by varying fluence of noble gas ion radiation.
We present the first comprehensive search for high-energy neutrino emission from high- and low-mass X-ray binaries conducted by IceCube. Galactic X-ray binaries are long-standing candidates for the source of Galactic hadronic cosmic rays and neutrinos. The compact object in these systems can be the site of cosmic-ray acceleration, and neutrinos can be produced by interactions of cosmic rays with radiation or gas, in the jet of a microquasar, in the stellar wind, or in the atmosphere of the companion star. We study X-ray binaries using 7.5 yr of IceCube data with three separate analyses. In the first, we search for periodic neutrino emission from 55 binaries in the Northern Sky with known orbital periods. In the second, the X-ray light curves of 102 binaries across the entire sky are used as templates to search for time-dependent neutrino emission. Finally, we search for time-integrated emission of neutrinos for a list of 4 notable binaries identified as microquasars. In the absence of a significant excess, we place upper limits on the neutrino flux for each hypothesis and compare our results with theoretical predictions for several binaries. In addition, we evaluate the sensitivity of the next generation neutrino telescope at the South Pole, IceCube-Gen2, and demonstrate its power to identify potential neutrino emission from these binary sources in the Galaxy.
In the past few decades, significant efforts have been made to improve the theoretical understanding of thermal transport mechanisms in thermal insulation materials and push the thermal conductivity's lower limits. However, most works focused singularly on specific types of materials, and the models used for thermal conductivity predictions are diverse - a model that fits one material might not fit others. Here, we improve and unify the gas and solid thermal conductivity models for porous materials. Through experimental characterization of several different materials as well as literature data for other materials, these models are validated. Further, we have also found that the pressure-dependent gas thermal conductivity of most materials can be well fitted by using one or two pore sizes without using a complex pore size distribution. With the refined models, we decompose the effective thermal conductivity of several thermal insulation materials into gas, solid, and radiation contributions. For cellular (polystyrene and polyurethane) foams, the relative contributions from air, solid, and radiation are 58–75%, 3–11%, 16–38%, respectively. For granular porous materials (polyurethane and silica in this work), the contributions from air, solid, and radiation are 45–66%, 34–46%, and 0–8%, respectively. This work is expected to provide guidance on the design and optimization of the next generation of thermal insulation materials, for example, through the effort of reducing gas conduction and radiation in foams and suppressing gas and solid conduction in aerogels.
Thermal radiation of the gaseous and particle phases in a pilot-scale pressurized oxy-combustor is computationally studied. In particular, the radiation characteristics of gases and particles are estimated by employing the statistical narrow-band model and the large-particle model. It is found that thermal radiation of the particle cloud dominates in the combustor under a furnace temperature of 1500 K and when there is no substantial loss of particles to the walls. Another important observation is that radiation from the gas and particles can be approximately treated as a graybody under these conditions. More specifically, the results on the spectral radiation intensity of a gas comprising 40% (vol) H 2 O and 60% CO 2 show that when the pressure is increased to 15 bar, and the radiation pathlength is 100 cm, the spectral radiation profile of the gas phase approaches that of a blackbody at the respective temperature. In addition, the emissivity of the particulate cloud has been evaluated as a function of the particle concentration and diameter by employing the large-particle model. It is shown that the emissivity grows with the particle concentration but decreases with the particle size for the same mass of the particles. Finally, this outcome of the present study is expected to be used to validate the assumption of the gray-gas model adopted in the numerical simulations of pressurized oxy-combustion.
Infrared thermography of the first wall in DIII-D is used to show the toroidal peaking of plasma radiation during mitigated disruptions with shattered pellet injection (SPI) and with massive gas injection (MGI). During MGI, the radiation peak location is shown to be due to the n = 1 magnetohydrodynamic (MHD) activity and continuously variable in toroidal phase based on externally applied error fields, consistent with previous experimental results. Furthermore, the measured toroidal peaking factor (TPF) is in agreement with the previous estimates based on radiometry, with a total TPF of 1.07 ± 0.05 when integrated over the entire duration of the disruption, and a value of 1.3 ± 0.1 during the thermal quench (TQ). For SPI, the location of the radiation peak is found to be determined predominantly by the particle source location, and thus can be varied on DIII-D by injecting from either of two toroidally separated injectors. 3D non-linear MHD simulations with the NIMROD code support this picture. The localized SPI particle source leads to higher peaking than for MGI, with a total TPF of 1.3 ± 0.1 over the full disruption and a TQ TPF of 1.9 +0.5/−0.3, a value also consistent with the NIMROD simulations. This TQ value is close to or potentially exceeding the allowable limit for the previously planned Be first wall in ITER, pointing to the importance of radiation asymmetries depending on the choice of first wall material.
The spark discharge of an aircraft plasma jet igniter is studied using high-fidelity numerical simulations and X-ray radiography measurements. The target problem here features the thermal expansion of hot gas introduced by the electric spark within a confined igniter cavity, which eventually evolves into a pulsed jet of a high-temperature kernel. A comprehensive set of models adapted from existing strategies for internal combustion engine spark plug discharge is extended to the target problem, including the modeling of energy deposition, plasma reactions, thermodynamic properties, and heat losses. A series of validation and parameter studies are performed and presented. The kernel size is found to be sensitive to heat losses arising from radiation and hot gas remained within the discharge cavity, rather than heat conduction to the wall in the discharge cavity. Depending on the enforced shape of the post-breakdown electric arc, the spark kernel can be off-centered, tilted, and considerably asymmetric. These features have been previously not considered when studying such igniter configurations and may have a first-order impact on the ignition process. Provided a proper setup of the heat loss models and electric arc shape, the numerical results are quantitatively comparable to the experimental results in terms of the kernel size, shape, and velocity throughout different stages after the spark discharge.
One of the most pressing questions in cosmology is how the black holes (BHs) powering quasars at high redshift grow to supermassive scales within a billion years of the Big Bang. Here we show that sustained super-Eddington accretion can be achieved for BHs with Eddington fractions f Edd ≳ 2/ϵ, where ϵ is the efficiency with which radiation is generated in the accretion process. In this regime, the radiation carries too little momentum to halt the accretion flow and the infalling gas traps the radiation. The BH growth then proceeds unimpeded until the gas supply is exhausted, in contrast to accretion at lower rates, which is limited by the radiation generated in the accretion process. The large gas supply available in massive high-redshift quasar host galaxies may be readily accreted onto seed BHs via this supply-limited mode of accretion, providing an explanation for how such supermassive BHs are assembled in the early universe. This sustained super-Eddington growth may also explain the short lifetimes inferred for the H ii regions surrounding high-redshift quasars, if the bulk of the BH growth occurs without the associated radiation escaping to ionize the intergalactic medium. It furthermore implies that a population of obscured rapidly growing BHs may be difficult to detect, perhaps explaining why so few quasars with Eddington fractions higher than a few have been observed. Finally, this simple condition for sustained super-Eddington growth can easily be implemented in cosmological simulations that can be used to assess in which environments it occurs.
To present a methodology to use pulmonary gas exchange maps to guide functional avoidance treatment planning in radiation therapy (RT) and evaluate its efficacy compared with ventilation-guided treatment planning.