Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “normalizing flows”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

A block-spectral adaptive H-/$p$-refinement strategy for shock-dominated problems

An adaptive H-/p-refinement strategy using a novel sensor is devised and tested in a block-spectral compressible Euler code equipped with adaptive-mesh refinement (AMR) and high-order flux-reconstruction numerics. At each Gauss quadrature point (or solution point) within each spectral block (or mesh element) the discrete velocity jump ΔU = ∂U/∂y 1 Δy 1 + ∂V/∂y 2 Δy 2 + ∂W/∂y 3 Δy 3 is calculated and normalized by the local speed of sound, a. Here, the grid spacing, Δx i , is calculated in each direction as the distance between auxiliary Gauss-Lobatto points, staggered relative to the solution points. The polynomial order is increased from p = 0 to p = p max in regions of weak compression, (ΔU/a) crit < ΔU/a < 0 and kept at p = p max in regions of flow expansion ΔU/a ≥ 0, while staying at the H = 0 base mesh level. Regions experiencing strong compressions, i.e. ΔU/a < (ΔU/a) crit , are H-refined up to H = H max where H max is applied at the location of maximum compression, ΔU/a = min(ΔU/a) in the domain, while keeping p = 0 to guarantee robustness and monotonicity of the solution in the H refined region. The critical value of (ΔU/a) crit = -0.06 is found to effectively separate smooth and non-smooth solution regions, supported by a 1D detonation initiation test case in ideal gas and a shock-to-detonation transition in high explosives. Using this value, the Sod shock tube, Shu-Osher problem, double Mach reflection and a 2D detonation in a high-explosive are simulated with the proposed adaptive H-/p-refinement. In the Sod shock tube case, p-refinement resolves the (weak) contact discontinuity while H-refinement enhances the grid resolution in the shock exploiting the monotonicity of the p = 0 reconstruction. For the Shu-Osher problem, p-refinement captures the small-scale oscillations trailing the shock that would be otherwise attenuated, while H-refinement triggered by the ΔU-sensor appropriately tracks the shock. In the double Mach reflection problem, H-refinement confines the numerical diffusion around the reflected shock while p-refinement recaptures many physical features trailing the shock. Finally, in the 2D high-explosive detonation case, H-refinement follows the leading shock and resolves the curvature of the detonation wave, while p-refinement adds resolution to the trailing reaction zone. Finally, the proposed methodology is tested in a detonation-wave propagation test case in high-explosives with numerical predictions comparing favorably against experiments.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

The corrosion behavior of borosilicate glass in the presence of cementitious waste forms

Borosilicate glasses are widely used for radioactive waste disposal due to their ability to incorporate a variety of contaminants and radionuclides while exhibiting high durability in various disposal scenarios. This research evaluated the dissolution of borosilicate glass using both single-pass-flow-through (ASTM C1662-18) and product consistency test (ASTM C1285-21) methods with different solutions, including a cementitious-contacted water (called grout-contacted, GC, from this point) and solutions with varying levels of dissolved cementitious species such as Si, Ca, Al. Here, the results indicated that the presence of Ca plays a crucial role in suppressing glass corrosion, as evidenced by the slower normalized dissolution rates, which were one order of magnitude lower for boron and two orders of magnitude lower for rhenium, observed in both Ca-amended and GC solutions compared to the pH 12 buffer solution. This effect is attributed to the formation of a dense, low-porosity, and strongly bonded calcium silicate hydrate (CSH) layer on the glass surface, which implies that a glass corrosion process is influenced by ion exchange involving alkali ions Na + , K + , Ca 2+ , and hydrogen-containing species. A small number of glass particles treated in the GC solution showed minor corrosion pits in the form of shallow craters with an average diameter of approximately 500 μm. This observation is correlated with a significant reduction, 2000 to 3000 times lower, in the cumulative volume of glass pores, indicating that smaller pore voids were “sealed” in the presence of Ca 2+ ions, likely attributed to the formation of CSH precipitation or other corrosion products such as calcium carbonate saturated from the grout solution. These findings suggest that the presence of dissolved Ca in the GC solution can slow down the dissolution of borosilicate glass, contrary to the expected trend of higher dissolution rates resulting from exposure to high alkaline and thus higher pH solutions.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Mechanisms and models of the turbulent boundary layers at transcritical conditions

Computational models for high-pressure transcritical turbulence in wall-modeled large-eddy simulation typically rely on wall-function models to overcome the need for resolving the boundary layer structures. However, the mechanisms and models of turbulent boundary layer at transcritical conditions remain poorly understood since the near-wall flow and heat transfer are significantly affected by the enhanced fluctuations and steep gradients in thermodynamic properties. Here, to address this issue, we study the mechanisms of transcritical turbulent boundary layers and wall-attached models at transcritical conditions. It is shown that the real-fluid variable-property effects are associated with the coupling between the near-wall cycle and the wall-normal coherent motions in the outer layer, resulting in the amplification of turbulent energy in the outer layer and noticeable energy transfer between the log-layer and the outer layer; hence, turbulence in the log-layer is modulated by the outer layer. Based on this underlying physical principle, we propose the characteristic velocity and length scales for the attached eddy at transcritical conditions, and extend the attached eddy model to transcritical turbulent boundary layers by introducing a mixed scaling that incorporates both inner and outer scalings. We show that the new characteristic scales and extended attached-eddy model perform well in characterizing the structures in transcritical turbulent boundary layers.

Li, Fangbo↗
DOE OSTI2024

Measurement of flow coefficients in high-multiplicity 𝑝 + Au, 𝑑 + Au, and 3 He + Au collisions at $\sqrt{s_{NN}}$ = 200 GeV

Flow coefficients (𝑣 2 and 𝑣 3 ) are measured in high-multiplicity 𝑝 + Au, 𝑑 + Au, and 3 He + Au collisions at a center-of-mass energy of $\sqrt{s_{NN}}$ = 200 GeV using the STAR detector. The measurements utilize two-particle correlations with a pseudorapidity requirement of |𝜂| < 0.9 and a pair gap of |Δ⁢𝜂| > 1.0. The primary focus is on analysis methods, particularly the subtraction of nonflow contributions. Four established nonflow subtraction methods are applied to determine 𝑣 𝑛 , validated using the HIJING event generator. 𝑣 𝑛 values are compared across the three collision systems at similar multiplicities; this comparison cancels the final-state effects and isolates the impact of initial geometry. While 𝑣 2 values show differences among these collision systems, 𝑣 3 values are largely similar, consistent with expectations of subnucleon fluctuations in the initial geometry. The ordering of 𝑣 𝑛 differs quantitatively from previous measurements using two-particle correlations with a larger rapidity gap, which, according to model calculations, can be partially attributed to the effects of longitudinal flow decorrelations. As a result, the prospects for future measurements to improve our understanding of flow decorrelation and subnucleonic fluctuations are also discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Distributed Coordination of Networked Microgrids for Voltage Support in Bulk Power Grids

The increasing deployment of distributed energy resources (DERs) and microgrids (MGs) in power distribution systems has enabled the adjustment of reactive power consumption as seen at the substation, which can be used to provide voltage support for the bulk power system (BPS). Leveraging this new capability will provide greater resiliency to the power system as a whole. Here, the goal of this paper is to develop and compare three different algorithms, namely distributed optimal power flow, distributed consensus algorithm, and fully decentralized collaborative autonomy for unbalanced distribution systems for microgrid coordination. These algorithms use networked MGs to support the BPS voltage when a contingency at the bulk grid results in abnormally low voltages, which may be a precursor to voltage collapse. Our comparative analysis includes both qualitative and quantitative assessments of the three algorithms and a discussion of the trade-offs between the decentralized and distributed methods in normal and disrupted conditions. Each algorithm was evaluated on the modified IEEE 13-bus system and a real power distribution system at Chattanooga, Tennessee, that encompasses more than 4500 buses. Each algorithms excels differently and may be suited for different scenarios depending on the condition, operations, and priorities of the power and communication systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Data-driven closure modeling for hypersonic turbulent flows

The Reynolds-averaged Navier–Stokes (RANS) equations remain a workhorse technology for simulating compressible fluid flows of practical interest. Due to model-form errors, however, RANS models can yield erroneous predictions that preclude their use on mission-critical problems. This report summarizes work performed from FY22-FY24 focused on improving RANS models for hypersonic flows using data-driven modeling and scientific machine learning. In this work we: 1. Investigate the current capabilities of RANS models in Sandia’s parallel aerodynamics and re-entry code (SPARC) for hypersonic flows with a focus on shock boundary layer interactions (SBLIs), 2. Assess several established corrections that exist in the literature aimed at improving predictions for SBLIs, 3. Develop improved models for the Reynolds stress tensor using tensor-basis neural networks, 4. Develop a neural-network-based variable turbulent Prandtl number model to reduce errors in wall heating in SBLIs. 5. Begin future investigations including employing the LIFE framework to improve wall heating predictions in SBLIs as well as the ensemble Kalman filter. We find that current RANS models in SPARC are deficient for complex SBLI flows. In particular, no current model jointly predicts wall heat flux, wall shear stress, and wall pressure with reasonable accuracy. Existing corrections help, but do not alleviate this issue altogether. The development of improved models for the Reynolds stress tensor via tensor-basis neural networks results in more predictive RANS models across a suite of low-speed and high-speed cases. For hypersonic boundary layers, the inclusion of the wall-normal Reynolds stress via TBNNs has an appreciable impact on the wall-normal momentum balance and wall quantities. However, we find that improvements to the Reynolds stress tensor do not address the over-prediction in wall heat flux in SBLIs. We find that a neural-network-based variable turbulent Prandtl number model systematically and substantially improves wall heating predictions for a range of SBLI cases.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Modeling Radionuclide Inventories in MSR Off-Gas Systems with Radiochemical Transport Analysis

Molten salt reactors (MSR) contain unique characteristics that may require enhancements to modeling tools to accurately predict phenomena. One characteristic that may be advantageous to leverage during normal operation is on-line processing of the circulating fuel salt, such as an off-gas system (OGS) to remove volatile fission products. Therefore, new modeling tools must be developed to integrate spatial resolution and chemistry effects into fuel depletion tools to be able to account for these non-core sources of radioactivity. Such types of radiochemical transport analysis tools were used to estimate the removal rates for 12 elements within a flow model of the Molten Salt Reactor Experiment (MSRE) by optimizing against legacy experimental data of the gas-borne (GB) percentages of 12 nuclides. The removal rates were used in a depletion model to calculate the FP inventory that enters the MSRE OGS. Calculations are in good agreement with the empirical GB percentages reported for the 12 nuclides, which validates the approach and verifies each tool’s treatment of the radiochemical flow effects. The OGS inventory is discussed in terms of the largest nuclide contributors to activity, dose consequence, decay heat, and elemental composition. Finally, insights from the study allow recommendations to be made for future code development activities.

Shahbazi, Shayan↗
DOE OSTI2025

A comparative TEA of a two-step process using chemical solvents for producing an ultra-sweet natural gas

Here, a comprehensive Techno-Economic Analysis (TEA) was performed to evaluate the economic feasibility of a novel two-step process (TSP) developed in Aspen Plus V12.1 to desulfurize and decarbonize a raw natural gas containing (2 mol% H 2 S and 5 mol% CO 2 ) into an ultra-sweet natural gas containing (1.72 ppmv H 2 S and 4.19 ppmv CO 2 ). The raw natural gas flow rate used in the TSP was 117.74 kg/s at 60 °C and 50 bar. The TSP combines an H 2 S desulfurization step using potassium carbonate (K 2 CO 3 ) and a CO 2 capture step using 3 different chemical solvents, monoethanolamine (MEA), sodium glycinate (SGS), and potassium glycinate (PGS). Both steps employ fixed-bed absorbers packed with Mellapak 250Y structured packing. The hydraulics and mass transfer characteristics for the TSP were calculated, indicating normal operation with higher gas-side (k G ) than liquid-side (k L ) mass transfer coefficients. The TEA of TSP indicated that PGS had the most promising economic feasibility among the 3 solvents as it exhibited the lowest Levelized Cost of CO 2 capture (LCOC) of $\$$47.54/ton.CO 2 at a Capital Expenditure (CAPEX) of $\$$24.98 million, and an Operating Expenditure (OPEX) of $\$$12.20 million/year. Also, the TSP could produce one MMSCF of ultra-sweet natural gas at a total cost of $\$$339.55.

03 NATURAL GAS↗
DOE OSTI2026

Mechanistic Modeling of TEG Dehydrator Emissions in Oil and Gas Industry

This work presents a mechanistic modeling approach for simulating methane emissions from triethylene glycol (TEG) dehydrators used in oil & gas (O&G) operations. The model was developed as a modular component of the Mechanistic Air Emissions Simulator (MAES) tool, incorporating species-specific absorption and emission dynamics through two-level, second-order polynomial regression (PR) models trained on ProMax simulation data: (1) species-level regression models that track the transfer rates of individual gas species within the dehydrator unit streams, and (2) outlet flow stream regression models that predict the fraction of inlet gas distributed among the outlet streams of the dehydrator unit. These behaviors were characterized over a range of glycol circulation ratios, wet gas pressures, and temperatures. The model was validated using root mean square error (RMSE) analysis. The species-level PR achieved low root mean square error (RMSE) values (<0.03) for light hydrocarbon species across all dehydrator components, ranging from 0.0009 for methane to 0.029 for normal pentane. Similarly, the outlet-level PR yielded RMSE values below 0.002 for the dry gas fraction, 0.001 for the flash tank fraction, and 0.002 for the still vent fraction, demonstrating strong agreement between predicted and reference ProMax values. When deployed at field facilities, the model significantly improved MAES-simulated dehydrator emissions, revealing that gas-assisted glycol pump emissions are the dominant contributors to both dehydrator-level and site-level methane emissions under uncontrolled conditions. Further analysis of the 154 dehydrator units reported by operators under the AMI 2024 project showed that 54 units (31%) used gas-driven glycol pumps, of which 6 units (11%) operated with uncontrolled flash tanks, and 22 units (40.7%) were identified as potentially oversized. Of the six dehydrator units with uncontrolled gas-assisted pumps, pump emissions accounted for 90.25% of total dehydrator emissions and 63.10% of total site-level emissions. These findings highlight substantial opportunities for emissions mitigation through equipment upgrades.

MAES↗
DOE OSTI2024

Optimization of a Welding Procedure for Making Critical Aluminum Welds on the LBNF Absorber Core Block

The LBNF Absorber consists of thirteen 6061-T6 aluminum core blocks. The core blocks are water cooled with de-ionized (DI) water which becomes radioactive during beam operations. The cooling water flows through gun-drilled channels in the core blocks. A weld quality optimization was performed to produce National Aeronautical Standard 1514 Class I quality welds on the aluminum core blocks. This was not successful in all cases. An existing Gas Tungsten Arc Welding Procedure Specification was fine tuned to minimize, in most cases, and eliminate detect-able tungsten inclusions in the welds. All the weld coupons, however passed welding inspection as per the piping code: ASME B31.3 Normal Fluid Service. Tungsten electrode diameter, type, and manufacturer were varied. Some of the samples were pre-heated and others were not. It was observed that larger diameter electrodes, 5/32 in., with pre-heated joints resulted in welds with the least number of tungsten inclusions. It is hypothesized that thinner electrodes breakdown easily and get lodged into the weld pool during the welding process. This breakdown is further enhanced by the large temperature differential between the un-preheated sample and the hot electrode.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Multi-X-point radiation and its dynamics in up/down asymmetry in plasma detachment regimes in Wendelstein 7-X

In Wendelstein 7-X (W7-X), stable plasma detachment in ECRH hydrogen plasmas is routinely achieved with the standard magnetic configuration, containing five magnetic islands in the scrape-off layer (SOL). Plasma detachment induced by intrinsic carbon impurities is characterized by a high radiation fraction (ƒ rad ) with significantly reduced divertor heat load and particle flux. In this work, we present the radiation dynamics during the detached plasma (DP) phase, such as the shift of radiation zones in radial and poloidal directions towards the X-points as ƒ rad increases (∼0.6–0.9), as well as the penetration of the radiation layer into the confinement region in the deep DP phase with ƒ rad approaching unity. In particular, a structure of multi-X-point radiation (multi-XPR) with an up/down asymmetry in the DP phase is highlighted, which is revealed by bolometer tomography and is further confirmed by video diagnostics. The multi-XPR structure forms helical 3D bands aligned with W7-X field periodicity. Field reversal experiments show that the brightest XPR displaces between the upper and lower SOL regions, suggesting the E× B drift effects. This paper presents the multi-XPR structure observed in the W7-X plasma in detail for the first time. A simplified model considering the influence of the poloidal E× B drift ( V d ) on the impurity flow in the SOL shows that the downstream drift toward the target or target shadow region ( V d < 0) decreases the impurity content, while the upstream drift toward the LCFS ( V d > 0) increases the impurity content. On this basis, the poloidal drift potentially leads to an up/down asymmetry of impurity density in the SOL despite the symmetry magnetic topology. The dynamics of the up/down asymmetry in the multi-XPR structure is also related to the magnitude V d /D (normalized to the impurity diffusivity), with an additional effect owing to the radial inward shift of the emission zone. These results provide new insights into impurity-induced detachment dynamics and provide a basis for improving 3D modeling of impurity transport, which typically does not consider drifts.

E × B drift↗
DOE OSTI2025

Expansion of the Fast Neutron Hodoscope at TREAT to Support Fuel Safety Experiments

In 2024 the Fuel Motion Monitoring System (FMMS), or Hodoscope, at the Transient Reactor Test Facility, is being expanded from 96 viewing channels to 192 viewing channels, effectively doubling the FMMS field of view. Here, this increase in capability will allow the FMMS to support larger scale fuel tests, encompassing height-of-core test devices, test devices with multi-pin fuel assemblies, and test devices with recirculating coolant flow. Work supporting the expansion included refurbishing 96 additional proton recoil scintillator (PRS) detectors, doubling the data acquisition system (DAS) installed architecture, improving time synchronization in the DAS, and new research to measure the PRS detector energy-dependent, fast-neutron detection efficiency. In addition, laboratory activities have produced an improved benchtop testing capability for assessing the DAS, time synchronization, and external start triggering, along with an updated capability to scan PRS detectors to develop a preliminary flat-field normalization prior to deployment to TREAT.

46 - INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AN↗
DOE OSTI2024

Optimization of a Welding Procedure for Making Critical Aluminum Welds on the LBNF Absorber Core Block

he LBNF Absorber consists of thirteen 6061-T6 aluminum core blocks. The core blocks are water cooled with de-ionized (DI) water which becomes radioactive during beam operations. The cooling water flows through gun-drilled channels in the core blocks. A weld quality optimization was performed to produce National Aeronautical Standard (NAS) 1514 Class I quality welds on the aluminum core blocks. This was not successful in all cases. An existing Gas Tungsten Arc Welding (GTAW) Welding Procedure Specification (WPS) was fine tuned to minimize, in most cases, and eliminate detectable tungsten inclusions in the welds. All the weld coupons however, passed welding inspection as per the piping code: ASME B31.3 Normal Fluid Service. Tungsten electrode diameter, type, and manufacturer were varied. Some of the samples were pre-heated and others were not. It was observed that larger diameter electrodes, 5/32 in., with pre-heated joints resulted in welds with the least number of tungsten inclusions. It is hypothesized that thinner electrodes breakdown easily and get lodged into the weld pool during the welding process. This breakdown is further enhanced by the large temperature differential between the un-preheated sample and the hot electrode.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Assessment of multi-stage telescopic control rods in a pebble bed reactor. Part II: corrosion and wear performance

To support High-Temperature Gas-Cooled Reactor designs, the effects of corrosion and wear on a SS316L telescopic control rod prototype were investigated in helium at 775 °C with up to 5 ppm oxygen for durations up to 500 h. Three deployment conditions were tested to assess their impact on wear and oxide formation: severe cycles, daily cycles, and gravity-driven deployments. Here, a chromium-rich oxide layer developed on regions exposed to the helium flow, reaching up to 2 µm in the most affected areas. No carburization or decarburization was observed. Oxide breakthrough occurred mainly at grooves and contact points during severe cycling of the telescopic rod, while in the other tests it was limited to contact points. Friction forces increased due to oxide growth and thermal expansion, the latter being the dominant factor, but never enough to hinder rod deployment significantly. Coaxial misalignment of stages was identified as the main potential risk to reliable operation under normal conditions.

Control rod↗
DOE OSTI2024

Investigation of two-dimensional radio-frequency sheath properties using a microscale fluid model

Abstract In previous work (Kohno and Myra 2023 Comput. Phys. Commun. 291 108841), we developed a numerical scheme based on a two-dimensional microscale radio-frequency (RF) sheath model with periodically curved wall boundaries. Here, we expand the capability of this scheme through modification of the boundary conditions (BCs) on the conducting walls, which allows the ion flow to turn back to the plasma at locations on the walls where the electromagnetic force on the ions is reversed from its usual direction. Numerical simulations are carried out to investigate the dependences of the surface-integrated admittances on the wall bump height, ion magnetization, ion mobility, and the magnetic field angle, and to visualize the sheath structures in several cases. One of the main results is the ion cyclotron admittance resonance observed under the condition of low ion mobility (high normalized frequency). It is shown that the amplitude of the resonance peak depends on the wall bump height and the ion velocity is reversed on the sides of the bump in an RF cycle for the resonance cases. Furthermore, the differences in the admittances between the one- and two-dimensional microscale models are assessed for the purpose of understanding non-locality of the sheath near the wall surface for the parameters considered in this study. This information will be essential for improving the sheath BC for macroscale calculations in the future.

Kohno, H. (ORCID:0000000332690010)↗
DOE OSTI2025

Development and performance of high-temperature superconducting CORC ® cables with CFD REBCO tapes

Abstract Increasing the normal zone propagation velocity (NZPV) in superconducting cables based on rare earth barium copper oxide (REBCO) tapes is expected to cause a paradigm shift by enabling the development of faster and more reliable quench detection systems to mitigate the development of destructive hot spots. Furthermore, a higher NZPV is beneficial in terms of fault current limitation capabilities by accelerating the homogenization of the quench in superconducting power devices. One way envisioned to increase the NZPV of cables based on REBCO tapes is the current flow diverter (CFD) concept. A Conductor on Round Core (CORC ® ) cable made with CFD REBCO tapes, called CFD CORC ® cable, and a CORC ® cable made with regular REBCO tapes, were fabricated and tested. The critical current of the CFD and regular CORC ® cables were obtained at temperatures ranging from 67 to 77 K in self-field. Measurements showed that the NZPV was increased by a factor of 4–4.5 in the case of the CFD CORC ® cable. Furthermore, the results suggest that the NZPV depends only on the applied current, similar to what has been observed previously on single REBCO tapes. These results demonstrate the successful integration of CFD REBCO tapes in CORC ® cables, without compromising their superconducting properties or their enhanced NZPV. The possibility of enhancing the NZPV of REBCO-based cables could facilitate quench detection, which remains an important issue in HTS magnet applications. This advancement also holds promise for the fault current limitation functionality of CORC ® cables.

Ben Saâd, Haïfa (ORCID:0000000216205364)↗
DOE OSTI2025

Evolution of in-medium baryon-baryon scattering cross sections and stiffness of dense nuclear matter from Bayesian analyses of FOPI proton-flow excitation functions

Within a Bayesian statistical framework using a Gaussian process (GP) emulator for an isospin-dependent Boltzmann-Uehling-Uhlenbeck (IBUU) transport model simulator of heavy-ion reactions with momentum-independent Skyrme interactions, we infer from the proton directed and elliptical flow in mid-central Au + Au reactions, measured at beam energies from 150 to 1200 MeV/nucleon by the FOPI Collaboration, the posterior probability distribution functions (PDFs) of the in-medium baryon-baryon scattering cross section (BBSCS) modification factor 𝑋 (with respect to their free-space values) and the stiffness parameter 𝐾 of dense nuclear matter. Here, we find that the most probable value of 𝑋 evolves from around 0.7 to 1.0 as the beam energy 𝐸 beam /𝐴 increases. On the other hand, the posterior PDF⁡(𝐾) may have dual peaks having roughly the same height or extended shoulders at high 𝐾 values. More quantitatively, the posterior PDF⁡(𝐾) changes from having a major peak around 220 MeV, characterizing a soft EOS in the reaction at 𝐸 beam /𝐴 = 150 MeV, to one that peaks around 320 MeV, indicating a stiff EOS in the reactions at 𝐸 beam /𝐴 higher than about 600 MeV. The transition from soft to stiff happens in mid-central Au + Au reactions at beam energies around 250 MeV/nucleon in which 𝐾 = 220 MeV and 𝐾 = 320 MeV are approximately equally probable. Altogether, the FOPI proton-flow excitation function data indicate a gradual hardening of hot and dense nuclear matter as its density and temperature increase in reactions with higher beam energies.

nuclear matter↗
DOE OSTI2025

Advances in understanding vacuum break dynamics in liquid helium-cooled tubes for accelerator beamline applications

Understanding air propagation and condensation following a catastrophic vacuum break in particle accelerator beamlines cooled by liquid helium is essential for ensuring operational safety. This review summarizes experimental and theoretical work conducted in our cryogenics lab to address this issue. Systematic measurements were performed to study nitrogen gas propagation in uniform copper tubes cooled by both normal liquid helium (He I) and superfluid helium (He II). These experiments revealed a nearly exponential deceleration of the gas front, with stronger deceleration observed in He II-cooled tubes. To interpret these results, a one-dimensional (1D) theoretical model was developed, incorporating gas dynamics, heat transfer, and condensation mechanisms. The model successfully reproduced key experimental observations in the uniform tube system. However, recent experiments involving a bulky copper cavity designed to mimic the geometry of a superconducting radiofrequency (SRF) cavity revealed strong anisotropic flow patterns of nitrogen gas within the cavity, highlighting limitations in extrapolating results from simplified tube geometries to real accelerator beamlines. To address these complexities, we outline plans for systematic studies using tubes with multiple bulky cavities and the development of a two-dimensional (2D) model to simulate gas dynamics in these more intricate configurations. As a result, these efforts aim to provide a comprehensive understanding of vacuum breaks in particle accelerators and improve predictive capabilities for their operational safety.

Beamline tube↗