Assessing the Structure-Based Turbulence Model Performance for Thermal Striping Applications Using Symmetric Jet Experiments
Not Available
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.
Not Available
Not Available
Not provided.
A high-fidelity computational fluid dynamics (CFD) analysis was performed using the Large Eddy Simulation (LES) model for the lower plenum of the High–Temperature Test Facility (HTTF), a ¼ scale test facility of the modular high temperature gas-cooled reactor (MHTGR) managed by Oregon State University. In most next–generation nuclear reactors, thermal stress due to thermal striping is one of the risks to be curiously considered. This is also true for HTGRs, especially since the exhaust helium gas temperature is high. In order to evaluate these risks and performance, organizations in the United States led by the OECD NEA are conducting a thermal hydraulic code benchmark for HTGR, and the test facility used for this benchmark is HTTF. HTTF can perform experiments in both normal and accident situations and provide high-quality experimental data. However, it is difficult to provide sufficient data for benchmarking through experiments, and there is a problem with the reliability of CFD analysis results based on Reynolds–averaged Navier–Stokes to analyze thermal hydraulic behavior without verification. To solve this problem, high-fidelity 3-D CFD analysis was performed using the LES model for HTTF. It was also verified that the LES model can properly simulate this jet mixing phenomenon via a unit cell test that provides experimental information. As a result of CFD analysis, the lower the dependency of the sub-grid scale model, the closer to the actual analysis result. In the case of unit cell test CFD analysis and HTTF CFD analysis, the volume-averaged sub-grid scale model dependency was calculated to be 13.0% and 9.16%, respectively. As a result of HTTF analysis, quantitative data of the fluid inside the HTTF lower plenum was provided in this paper. As a result of qualitative analysis, the temperature was highest at the center of the lower plenum, while the temperature fluctuation was highest near the edge of the lower plenum wall. The power spectral density of temperature was analyzed via fast Fourier transform (FFT) for specific points on the center and side of the lower plenum. FFT results did not reveal specific frequency-dominant temperature fluctuations in the center part. It was confirmed that the temperature power spectral density (PSD) at the top increased from the center to the wake. The vortex was visualized using the well-known scalar Q-criterion, and as a result, the closer to the outlet duct, the greater the influence of the mainstream, so that the inflow jet vortex was dissipated and mixed at the top of the lower plenum. Additionally, FFT analysis was performed on the support structure near the corner of the lower plenum with large temperature fluctuations, and as a result, it was confirmed that the temperature fluctuation of the flow did not have a significant effect near the corner wall. In addition, the vortices generated from the lower plenum to the outlet duct were identified in this paper. It is considered that the quantitative and qualitative results presented in this paper will serve as reference data for the benchmark.
This report reviews the modeling and simulation capabilities of Argonne National Laboratory’s DASSH code that is used in present reactor analysis activities. These capabilities will be used to establish the set of verification tasks necessary to verify DASSH for use on commercial projects. A similar approach was taken for the PERSENT, REBUS and DIF3D software packages. The DASSH program is a thermal analysis code designed to rapidly allow a reactor design engineer to obtain flow rates requirements that satisfy peak temperature constraints in the domain. DASSH is a follow-on development to the SE2-ANL software and SUPERENERGY-2 software that it is based upon. DASSH was designed to account for both neutron and gamma heating and is inherently connected to the GAMSOR part of the ARC suite of fast reactor analysis software. SE2-ANL is a developed piece of software from the 1980s while DASSH is a modern implementation with notable improvements in geometry handling. The most important upgrade of DASSH relative to SE2-ANL is that it can analyze multiple time points in a single run where SE2-ANL can only treat a single time point. This allows the user to understand the impact of and search the flow distribution for the entire operational period of a reactor design considering pressure drop, peak coolant and fuel temperatures, and thermal striping. DASSH has three input paths that have to be verified. The first input path builds the geometry and power distribution based upon the DIF3D model but ignores the gamma heating aspects of the problem. The second input path also builds the geometry from the DIF3D model but it takes the neutron and gamma heating distributions from GAMSOR. The third input path is to take the geometry and power distribution directly from user input (i.e. not coupled to DIF3D or GAMSOR). DASSH also has many built in correlations for material properties along with a user defined specification of the fuel, structure, and coolant properties. There are correlations for flow split, mixing, pressure drop, and heat transfer coefficients (subchannel rather than a direct methodology). In total, verification of DASSH will require an extensive testing to cover all possible user features of the software.
We report the results of synchrotron Bragg Coherent x-ray Diffraction Imaging (BCDI) experiments to investigate domain formation in a micron-sized magnetite crystal undergoing the Verwey transition at low temperature. A strong splitting of the measured 311 Bragg reflection was observed in the low-temperature phase, indicating the formation of domains. BCDI revealed pronounced strain distributions, characterized by a clear layered stripe domain structure in real-space. Stripes were seen only along the [001] crystallographic direction, normal to the substrate surface direction, breaking the symmetry of the cubic high-temperature phase. It is argued that other domain directions were suppressed by the sample mounting orientation. More surprisingly, only a single domain orientation was observed, suggesting an additional symmetry-breaking influence originating from the shape of the crystal. Gaining insight into how thermal effects induce the formation of layered or striped phases offers a valuable framework for understanding the development of mesoscopic domains and strain patterns in functional materials.
Context.Pulsars are detected over the whole electromagnetic spectrum, from radio wavelengths up to very high energies, in the GeV-TeV range. While the radio emission site for young pulsars is well constrained to occur at altitudes about several percent of the light-cylinder radius andγ-ray emission is believed to be produced in the striped wind, outside the light cylinder, the non-thermal X-ray production site remains unknown. Aims.The aim of this Letter is to localize the non-thermal X-ray emission region based on a multi-wavelength pulse profile fitting for PSR J2229+6114, which stands as a particularly good candidate given its high X-ray brightness. Methods.Based on the geometry deduced from the joint radio andγ-ray pulse profiles, we fixed the magnetic axis inclination angle and the line-of-sight (LoS) inclination angle. However, we left the region of X-ray emission unlocalized, setting it somewhere between the surface and the light cylinder. We localized this region and its extension by fitting the X-ray pulse profile as observed by the NICER,NuSTAR, and RXTE telescopes in the ranges of 2–7 keV, 3–10 keV, and 9.4–22.4 keV, respectively. Results.We constrained the non-thermal X-ray emission to arise from altitudes between 0.2 r L and 0.55 r L wherer L is the light-cylinder radius. The magnetic obliquity is approximatelyα ≈ 45° −50° and the LoS inclination angle isζ ≈ 32° −48°. Conclusions.This Letter is among the first works to offer a tight constraint on the location of non-thermal X-ray emission from pulsars. We plan to apply this procedure to several other promising candidates to confirm this new result.
The University of Missouri Research Reactor (MURR) is expected to be converted from highly enriched uranium (HEU, ≥ 20 wt% U-235) U-Al x dispersion fuel to low-enriched uranium (LEU, < 20 wt% U-235) with U-10Mo monolithic fuel. This work introduces high-fidelity irradiation thermo–mechanical (T-M) analysis of the MURR LEU focusing on changes in coolant channel gap thickness. Three-dimensional (3D) finite element (FE) models were developed to simulate the irradiation T-M behavior of the MURR LEU element with all 23 curved fuel plates, the two side plates, and the combs. It was shown that channel gap thickness changes were influenced not only by plate thickness variations due to fuel swelling and creep but also by the radial displacement of consecutive MURR LEU plates. Modeling the fuel element assembly captured side plate displacements, which were shown to reduce radial fuel plate displacements towards the convex side. The maximum local radial displacement in the element was predicted at the end of life (EOL) as 23.7 mil (602.0 µm) on the lateral centerline of plate 23 towards the convex side. The maximum stripe-averaged reduction in channel gap thickness, particularly relevant for thermal hydraulics (TH) safety analysis, was calculated as 15.9 mil (403.9 µm) in single-side heated channel 24 (the outermost channel). These results account for the thermal resistance from the oxide build-up on cladding surfaces which was shown to be up to 0.82 mil (20.8 µm) thick. It was demonstrated that accounting for oxide layer thermal resistance led to a 10 °C higher peak fuel temperature and a 4.4 mil (111.8 µm) greater maximum local radial displacement. The impact of the calculated channel gap thickness changes on the MURR LEU TH safety analysis is evaluated in Part II.
Intercalation of metal atoms at the SiC(0001)-graphene (Gr) interface can provide confined 2D metal layers with interesting electronic properties. The intercalated Pb monolayer (ML) has shown the coexistence of the Gr(10 x 10)-moiré and a stripe phase, which still lacks understanding. Using density functional theory calculation and thermal annealing with ab initio molecular dynamics as motivated by experiment, we have studied the formation energy of Gr/Pb/SiC(0001) for different Pb coverages. Near the coverage of a Pb(111)-like ML mimicking the (10 x 10)-moiré, we find a slightly more stable stripe structure, where one half of the structure has compressive strain with Pb occupying the Si-top sites and the other half has tensile strain with Pb off the Si-top sites. This stripe structure along the Gr zigzag direction has a periodicity of 2.3 nm across the [1$\overline{2}$10] direction agreeing with the previous observations using scanning tunneling microscopy. Analysis with electron density difference and density of states show the tensile region has a more metallic character than the compressive region, while both are dominated by charge transfer from Pb ML to SiC(0001). As a result, the small energy difference between the stripe and Pb(111)-like structures means the two phases are almost degenerate and can coexist, which explains the experimental observations.
The impetus for this work is the Proliferation Resistance Optimization (PRO-X) program, which was created by the National Nuclear Security Administration (NNSA) of the US Department of Energy (DOE) to provide a framework for developing reactor designs that minimize the production of special nuclear materials while maximizing performance for peaceful uses. The Argentine company INVAP and Argonne National Laboratory collaborated to share knowledge and expertise regarding the thermal analysis of research reactors. Working together, they designed a series of test problems based on the Open Pool Australian Light Water (OPAL) reactor, whose fuel is in the form of flat plates cooled by water flowing through narrow rectangular channels. ANL used their PLTEMP/ANL code. INVAP used their TERMIC code (and its new multi-plate version, TERMIC-MP). Then, INVAP and ANL compared the two sets of results. 1-D models were used during the collaboration since for research reactor analysis, largely 1-D thermal-hydraulics models are commonly used because they are simple to apply and technically defendable. The collaboration between the two organizations aims for improvements in thermal-hydraulic models. Improvements to methods and models can reduce modeling uncertainties and excessive conservatism, allowing greater reactor performance without reducing predicted safety margins. Therefore, these improvements potentially expand the design space, making additional designs feasible, which, in turn, may achieve better fuel utilization and proliferation resistance in the designs. This report provides justification for the lateral node size (stripe width) chosen for the 1-D computer models during the collaboration. The lateral heat conduction in OPAL fuel plates is evaluated using the two-dimensional conduction capability of PLTEMP/ANL to determine a representative hot-stripe width that can be conservatively used for one-dimensional steady-state thermal-hydraulic analysis. Basically, a representative hot-stripe width needs to be selected so that the actual heat flux averaged over this stripe width is greater than or equal to the maximum heat flux obtained by the two-dimensional analysis. Detailed results in this report show that a hot-stripe width of 8.125 mm can be conservatively used for the one-dimensional thermal-hydraulics analysis of all plates. Therefore, eight equal stripes (8.125 mm each) are recommended to be modeled in the MCNP and PLTEMP/ANL full core models of the OPAL reactor.
A two-dimensional van der Waals ferromagnet Fe 5–x GeTe 2 (F5GT) exhibits room-temperature magnetic transition and magnetic anisotropy. Studies have identified diverse magnetic states, including stripe domains, skyrmionic type-I and topologically trivial type-II bubbles, depending on the magnetic and thermal history. Yet, the underlying micromagnetic energetics driving these states remains unclear. Here, we establish the magnetic phase diagram of bulk F5GT using cryogenic Lorentz four-dimensional scanning transmission electron microscopy. We reveal that stripe domains spontaneously form upon zero-field cooling, while metastable magnetic bubbles are stabilized by external fields. Transitions from type-I to type-II bubbles are driven by an oblique external field. Micromagnetic simulations confirm that these transitions arise from the interplay of uniaxial anisotropy, dipolar interactions, and external fields, without requiring significant Dzyaloshinskii–Moriya Interaction (DMI). In conclusion, these findings clarify the micromagnetic origin of spin textures of F5GT and establish it as a tunable platform for nanoscale topological magnetism.
We report the first numerical prediction of a “spin microemulsion”—a phase with undulating spin domains resembling classical bicontinuous oil-water-surfactant emulsions—in two-dimensional systems of spinor Bose-Einstein condensates with isotropic Rashba spin-orbit coupling. Using field-theoretic numerical simulations, we investigated the melting of a low-temperature stripe phase with supersolid character and find that the stripes lose their superfluidity at elevated temperature and undergo a Kosterlitz-Thouless-like transition into a spin microemulsion. Momentum distribution calculations highlight a thermally broadened occupation of the Rashba circle of low-energy states with macroscopic and isotropic occupation around the ring. Here, we provide a finite-temperature phase diagram that positions the emulsion as an intermediate, structured isotropic phase with residual quantum character before transitioning at higher temperature into a structureless normal fluid.
The need for reliable bonded interface materials is critical to realize the performance benefits of wide-bandgap devices in power electronic modules, especially in operating temperatures greater than 150 degrees Celsius. In this paper, we investigate the thermomechanical performance of sintered copper (Cu) as a large-area attachment, bonded between Cu baseplates and active-metal-bonded substrates, under accelerated thermal shock (-40 degrees Celsius to 200 degrees Celsius) conditions. In the fabrication phase of the samples, we used different stencil patterns and found out that the grid and stripe patterns resulted in a better outgassing of the residual organics during the sintering process, thereby ensuring a substantially improved bond quality than a full-area print. The paste consisted of Cu microflakes and we performed sintering using a Budatec SP300 sintering press at 275 degrees Celsius with 15 MPa of bonding pressure for 5 minutes in a nitrogen atmosphere. Under accelerated experiments, we monitored the degradation of the sintered Cu bond in the samples through C-mode scanning acoustic microscope (C-SAM) images. To quantify the defect percentage in C-SAM images, we investigated image denoising techniques to exclude the pattern prints. Finally, we cross-sectioned a sample and obtained digital microscope images, which revealed adhesive fracture as the dominant failure mechanism.
The need for reliable bonded interface materials is critical to realize the performance benefits of wide-bandgap devices in power electronics modules, especially in operating temperatures greater than 150 degrees C. In this paper, we investigate the thermomechanical performance of sintered copper (Cu) as a large-area attachment, bonded between Cu baseplates and active-metal-bonded substrates, under accelerated thermal shock (-40 degrees C to 200 degrees C) conditions. In the fabrication phase of the samples, we used different stencil patterns and found that the grid and stripe patterns resulted in a better outgassing of the residual organics during the sintering process, thereby ensuring a substantially improved bond quality compared to a full-area print. The paste consisted of Cu microflakes, and we performed sintering using a Budatec SP300 sintering press at 275 degrees C with 15 MPa of bonding pressure for 5 minutes in a nitrogen atmosphere. Under accelerated experiments, we monitored the degradation of the sintered Cu bond in the samples through C-mode scanning acoustic microscope (C-SAM) images. To quantify the defect percentage in C-SAM images, we investigated image denoising techniques to exclude the pattern prints. Finally, we cross-sectioned a sample and obtained scanning electron microscope images, which revealed adhesive fracture as the dominant failure mechanism.
Real-space Bragg coherent x-ray-diffractive imaging unveils the formation of bubble and stripe antiferromagnetic (AFM) phase domains on the surface of Ni 2 CoTeO 6 single crystal. The stripe domains exhibit dislocation-type topological defects. The defects typically form as defect-antidefect pairs and can be created or annihilated by temperature changes and local heating. Thermal fluctuations of the AFM domain walls are observed near the Néel temperature. Topological defect pairs dynamically form and disappear in the fluctuating state. These observations provide a real-space perspective on the dynamics of the AFM phase transition in a helical antiferromagnet. Finally, the remarkable tunability of AFM domain walls in a helical antiferromagnet underscores their potential for AFM spintronics applications.
2D van der Waals (vdW) magnets open landmark horizons in the development of innovative spintronic device architectures. However, their fabrication with large scale poses challenges due to high synthesis temperatures (>500 °C) and difficulties in integrating them with standard complementary metal-oxide semiconductor (CMOS) technology on amorphous substrates such as silicon oxide (SiO 2 ) and silicon nitride (SiN x ). Here, a seeded growth technique for crystallizing CrTe 2 films on amorphous SiN x /Si and SiO 2 /Si substrates with a low thermal budget is presented. This fabrication process optimizes large-scale, granular atomic layers on amorphous substrates, yielding a substantial coercivity of 11.5 kilo-oersted, attributed to weak intergranular exchange coupling. Field-driven Néel-type stripe domain dynamics explain the amplified coercivity. Moreover, the granular CrTe 2 devices on Si wafers display significantly enhanced magnetoresistance, more than doubling that of single-crystalline counterparts. Current-assisted magnetization switching, enabled by a substantial spin–orbit torque with a large spin Hall angle (85) and spin Hall conductivity (1.02 × 10 7 ℏ/2e Ω -1 m -1 ), is also demonstrated. These observations underscore the proficiency in manipulating crystallinity within integrated 2D magnetic films on Si wafers, paving the way for large-scale batch manufacturing of practical magnetoelectronic and spintronic devices, heralding a new era of technological innovation.
Heterogeneous landscapes can influence the development of convection through the generation of thermally driven mesoscale circulations. To assess the impacts of these circulations and their interaction with sea breezes, we simulated convection in an idealized coastal environment using the Regional Atmospheric Modeling System (RAMS). We compared simulations with striped patterns of surface vegetation to those of uniform vegetation to identify the importance of vegetation heterogeneity in impacting convective development. Under dry soil conditions representative of those during the Tracking Aerosol Convection Interactions Experiment (TRACER) and Experiment of Sea Breeze Convection, Aerosols, Precipitation, and Environment (ESCAPE) campaigns in June 2022, we found that these vegetation-induced circulations, referred to in the literature as “forest breezes,” are more important than the sea breeze in determining the location of convection initiation. Convection and precipitation are also found to be favored over forests and suppressed over pasture and suburban landscapes as a result of greater surface sensible heat flux over the forest. Our findings also indicate that forest breezes are important for initiating convection along the boundaries of the forest, but that cold pools may play a key role in propagating the forest breezes toward the center of the forest stripe. In our simulations, the collisions of these breezes in the center of the forest stripe lead to uplift and strong convection there; however, a different width of the forest stripe would alter when the forest breezes collide or whether they collide at all. The presence of these cold pools may therefore impact the “ideal stripe width,” the width of each vegetation stripe which maximizes domain-wide precipitation.