Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Averages”

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 1,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

The mechanism of twin thickening and the elastic strain state of TWIP steel nanotwins

A Twinning Induced Plasticity (TWIP) steel with a nominal composition of Fe-16.4Mn-0.9C-0.5Si-0.05Nb-0.05V was deformed to an engineering strain of 6%. The strain around the deformation twins were mapped using the 4D-STEM technique. Strain mapping showed a large average elastic strain of approximately 6% in the directions parallel and perpendicular to the twinning direction. However, the large average strain comprised of several hot spots of even larger strains of up to 12%. These hot spots could be attributed to a high density of sessile Frank dislocations on the twin boundary and correspond to shear stresses of 1–1.5 GPa. The strain and therefore stress fields are significantly larger than other materials known to twin and are speculated to be responsible for the early thickness saturation of TWIP steel nanotwins. The ability to keep twins extremely thin helps improve grain fragmentation, i.e. the dynamic Hall–Petch effect, and underpins the large elongations and strain hardening rates in TWIP steels.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A lightweight Fe–Mn–Al–C austenitic steel with ultra-high strength and ductility fabricated via laser powder bed fusion

Lightweight Fe–Mn–Al–C steels have become a topic of significant interest for the defense and automotive industries. These alloys can maintain high strength and ductility while also reducing weight in structural applications. Conventionally processed Fe–Mn–Al–C austenitic steels with high Al content (~9 wt%) demonstrate greater than 1.5 GPa strength with 35% elongation. Several recent studies have demonstrated success in fabricating steel parts using laser powder bed fusion (L-PBF) additive manufacturing (AM), which can generate near-net-shape components with complex geometries and is capable of local microstructural control. However, studies on L-PBF processing of Fe–Mn–Al–C alloys have focused on low Al content (<5 wt%) compositional regimes representing alloys that undergo transformation-induced plasticity (TRIP) and twinning-induced plasticity (TWIP). Here, in this study, we present the effects of L-PBF processing on the microstructure and mechanical properties of an Fe–30Mn–9Al–1Si-0.5Mo-0.9C austenitic steel. A process optimization framework is employed to determine an ideal L-PBF processing space that will result in >99% density parts. Implementing this framework resulted in near-fully dense specimens fabricated over a broad range of process parameters. Additionally, two bi-directional scan rotation strategies (90° and 67°) were applied to understand their effects on texture and anisotropy in this material. As-printed specimens displayed considerable work-hardening characteristics with average strengths of up to 1.3 GPa and 36% elongation in the build direction. However, solidification microcracks oriented in the build direction resulted in anisotropy in tensile strength and ductility resulting in average strengths of 1.1 GPa and 20% elongation perpendicular to the build direction. The successful L-PBF fabrication of Fe–30Mn–9Al–1Si-0.5Mo-0.9C presented here is expected to open new avenues for weight reduction in structural applications with a high degree of control over part topology.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Deformation behavior of metallic glass composites and plasticity accommodation at microstructural length-scales

Bulk metallic glass matrix composites represent a unique microstructural design strategy for overcoming the strength/ductility trade-off in structural alloys. Site-specific mechanical behavior of a Ti-based bulk metallic glass composite was evaluated at microstructural length-scale. The micro-pillars on the amorphous matrix showed on average a yield point of 1.9 GPa followed by serrated plastic deformation characteristic of shear banding. In contrast, micro-pillars on the crystalline dendritic phase showed a much lower yield strength of 0.72 GPa on average for the four different crystallographic orientations chosen, smooth plastic flow with no recognizable load burst, and stable strain-hardening after yield point. A number of micro-pillars were made at the interface between the glassy matrix and the crystalline dendrite. The mixed micro-pillars showed homogeneous deformation and greater plasticity which was attributed to their smaller stored elastic energy. Shear bands initiated in the amorphous matrix were arrested by the crystalline dendrite, which accommodated plasticity through slip bands and dislocations pile-ups. The interface remained intact after plastic deformation, with no observable signs of devitrification in the amorphous phase.

Bulk metallic glass matrix composite, Micro-pillar↗
DOE OSTI2026

A comprehensive first-principles study of the effects of the exchange-correlation functional and magnetism on defect and diffusion properties of the CoCrNi medium-entropy alloy

The present work is a novel, systematic study of the effect of density functional theory input parameters on the vacancy formation energy (VFE), migration barrier for diffusion, and electronic structure for each element in the CoCrNi medium-entropy alloy (MEA). In particular, the novelties include: (1) calculating the aforementioned properties of Co, Cr, or Ni, in the CoCrNi MEA using magnetic and non-magnetic states, and two versions of the generalized gradient approximation: Perdew, Burke, and Ernzerhof (PBE) and the PBE version for solids (PBEsol), and (2) a detailed comparison of 0 K activation energy to experimental creep activation energies. First-principles calculations at 0 K are performed using the Vienna ab-initio simulation package. Special quasirandom structures (SQS) and Widom-type substitution are employed. For each element, Co, Cr, or Ni, non-magnetic calculations result in a higher VFE and larger range of calculated values for the configurations studied. The averaged migration barrier is the highest for Co in the CoCrNi for three of four sets of calculation parameters in the configurations studied. Finally, the results indicate that the average 0 K activation energy for diffusion makes up 70–80% of the experimental creep activation energy, depending on the exchange-correlation functional employed.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Amorphous complexion-aided sintering enables scalable processing of bulk nanocrystalline Cu-Zr with high strength and compressive plasticity

Nanocrystalline alloys can have exceptional strengths, yet due to limited microstructural stability it is difficult to fabricate bulk pieces through traditional processing routes that retain nanosized grains. In this study, centimeter-sized Cu-Zr alloy pellets were fabricated via a simple and improved powder metallurgy processing route. Different consolidation temperatures and times were employed to investigate the effect of amorphous grain boundary complexions on densification and the resulting mechanical properties. Bulk compression tests were carried out, with the samples that were hot pressed at 900 °C for 10 h exhibiting an excellent combination of average yield strength of 722 ± 45 MPa and average failure strain of 25.3 ± 2.4 %. Furthermore, we find that a powder processing route which enables amorphous complexion-assisted sintering leads to samples that (1) reach full density without requiring quenching treatments or other complex processing, (2) demonstrate appreciable plasticity, and (3) have strength that competes with commercially available high-strength Cu alloys.

Bulk nanocrystalline alloy↗
DOE OSTI2023

Investigation of thermal hydraulic behavior of the High Temperature Test Facility's lower plenum via large eddy simulation

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.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

A suite of neutron time-of-flight detectors to measure hot-spot motion in direct-drive inertial confinement fusion experiments on OMEGA

A suite of six neutron time of flight detectors capable of measuring the primary DT neutron energy spectrum has recently been completed on the OMEGA Laser System. The detectors are positioned along multiple quasi-orthogonal lines of sight (LOS’s) in either a stand-alone, collinear, or antipodal configuration. The collinear detector configuration enables one to measure the neutron velocity along the detector LOS independent of the neutron-production history, while the antipodal detector configuration enables one to directly measure of the hot-spot velocity along the detector axis and Gamow velocity shift due to the plasma ion temperature. Using these six measurements of the primary DT neutron energy spectrum, the neutron-averaged hot-spot velocity and Gamow velocity shift of DT fusion neutrons have been measured in cryogenic experiments on OMEGA for the first time. Finally, the detector suite is operable between DT neutron yields of 1 × 10 13 to 2 × 10 14 with an uncertainty of ~14 km/s in the neutron-averaged hot-spot velocity measurement.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Anisotropic neutron response of trans-stilbene and impact on a handheld dual particle imager

Here this manuscript details measurements of the anisotropic light-output response of trans-stilbene, henceforth referred to as stilbene, along the three primary crystal planes and analyzes how this anisotropic response impacts neutron image reconstruction in a handheld dual particle imager (H2DPI). A prototype H2DPI composed of stilbene pillars (6 x 6 x 50 mm 3 ) coupled to silicon photomultipliers was built and is capable of imaging both fast neutrons and gamma rays from kilogram quantities of special nuclear material. Stilbene was chosen as the scintillating medium for the imager because of its pulse shape discrimination capability, time resolution and relatively high light output. A drawback to using stilbene, however, is the anisotropic response. A recoiling proton from a neutron elastic-scattering event in stilbene will yield different amounts of scintillation light depending on the direction of the recoiling proton with respect to the crystal lattice. This manuscript analyzes how this anisotropic response impacts neutron image reconstruction in 4π. The light output in the three primary crystal planes of stilbene were independently measured 6–7 times for proton recoil energies ranging from 0.5–5.0 MeV using a quasi-monoenergetic time-of-flight neutron source. The measured light-output data in the three crystal planes were fit with a semi-empirical function based on Birks’ formula. These fits were used to reconstruct the location of a 252 Cf source in front of the imager (0°, -0.8°), to the left of the imager (-90°, -0.8°) and directly above the imager (0°, 90°) to determine if it is necessary to apply the directionally-dependent response of stilbene to accurately reconstruct source locations. A bootstrapping technique was applied to the measured data sets to produce 1000 images composed of 2000 cone projections for each source location and each measured light-output curve. List-mode maximum likelihood expectation maximization was applied to each image; the highest pixel location and the full width at half maximum (FWHM) were then recorded. The average and standard deviations of these parameters were taken for each set of images. The maximum difference in the average azimuthal pixel location when varying between the appropriate and incorrect light-output curves for the source locations at (0°, -0.8°) and (-90°, -0.8°) were respectively found to be 1.85 ± 1.39° and 2.03 ± 1.78°. All other pixel locations and FWHMs were within a single standard deviation of uncertainty. Taken together, the anisotropic response of stilbene has negligible impact on the neutron image reconstruction capability of the H2DPI.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Machine learning framework for predicting uranium enrichments from M400 CZT gamma spectra

A machine learning framework was developed for predicting uranium enrichments from M400 CZT gamma spectra. This framework leverages the availability of a large amount of measured M400 gamma spectra and uses a recently updated version of Gamma Detector Response and Analysis Software (GADRAS) for gamma spectrum analysis and generation. It also leverages the existing machine learning modules in Python for gamma spectrum data processing, curation, model training, benchmarking, and optimization of the deep machine learning models. The framework is used to develop a deep learning model to analyze gamma spectra from a set of U 3 O 8 samples with enrichments ranging from 0.31 to 93.17% and UF 6 cylinders with enrichments ranging from 0.2 to 4.95%, and the model performance is tested using a set of measured spectra and the respective declared enrichment values. Results show that the model can correctly classify 99.35% of the U 3 O 8 sample enrichments, and can predict the samples’ enrichments within an average absolute error of 0.099% (in percentage points of enrichment). For the UF 6 cylinders, the average absolute error was approximately 0.03%, with an accuracy of 98% in classifying discrete enrichment values of UF 6 samples. Finally, the results also show that the model has performed significantly better in terms of predicting enrichments in UF 6 cylinders based on measured gamma spectra than the GEM code, with a standard deviation (of the relative errors) of 2.23% (compared with the 11.51% value for the GEM code) based on results from a set of test data.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Integration of high confinement, high poloidal beta plasma with dual radiated power and detachment controls for divertor protection and ELM suppression

Divertor detachment without serious core confinement quality loss in DIII-D’s high poloidal β scenario has been combined with impurity-induced ELM mitigation without disruption. Use of Ne previously granted access to a detached, non-ELMing regime that retained high confinement quality due to stimulation of Internal Transport Barrier (ITB) growth, but suffered from on-going core fuel dilution and high disruptivity. Excess Ne accumulation in the core plasma has now been avoided by feeding back core radiated power ($P_{rad}$) measurements to control Ne seeding, rather than using attachment fraction ($A_{frac}$) control with Ne; this also reduces disruptivity. At the same time, N 2 seeding is used in a feedback loop with $A_{frac}$ measurements, which previously posed low disruption risk. In this way, the effect of Ne in the core is managed while avoiding excess seeding, and N 2 acts to correct for any excess heat exhaust that might interfere with detachment. The average Ne flow rate was 38% of what was used in pure Ne $A_{frac}$ control, plus average N 2 flow that was 43% of pure N 2 $A_{frac}$ control, even while meeting an even deeper detachment target. Meanwhile, a steep pressure gradient in the core plasma reduces sensitivity to impurity-induced degradation of the pedestal and reasonable confinement quality was maintained despite operational challenges that blocked formation of an ITB in these experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

The effectiveness of D 2 pellet and gas injection in reducing intra-ELM tungsten erosion and heat flux in the DIII-D small angle slot divertor

Edge localized modes (ELMs) in H-mode plasmas can melt and erode plasma-facing components (PFCs) and lead to impurities in the core, reducing confinement. This study analyzes the use of D 2 pellet and gas injection for ELM mitigation on the DIII-D tokamak during the 2022 Small Angle Slot V-shaped Tungsten (W) (SAS-VW) divertor campaign, reducing W erosion and heat flux during ELMs. D α (656 nm) and WI (400.9 nm) filterscopes and Langmuir probes provide photon emission and electron density/temperature to estimate W atom erosion using the S/XB method, while surface eroding thermocouples measured ELM peak heat flux at the outer strike point (OSP). Thomson Scattering measurements of pedestal T e and n e provided input to predict W divertor erosion and heat flux during ELMs via the Free-Streaming plus Recycling Model (FSRM). While a greater D 2 mass injection rate decreased the ELM peak heat flux, the impact on W erosion was not monotonic. The average ‘large’ intra-ELM W erosion was lower for plasma shots with D 2 mass injection compared to the plasma shots without any D 2 mass injection in the SAS-VW divertor. However, plasma shots did not see significant changes in ‘large’ intra-ELM W divertor erosion as the D 2 mass injection rate increased. On average, the FSRM overestimated SAS-VW experimental intra-ELM W erosion by a factor of 6.9, but the overestimation is reduced to a factor of 3.6 with the implication of C deposition effects. Generally, experimentally measured and predicted quantities were worse for plasma shots with a lower D 2 mass injection rate and corresponding higher plasma carbon (C) impurity percentage (f C ). The discrepancies are postulated to be due to C/W material mixing, for which a simple analytic mixed-material model is presented. These results highlight the importance of incorporating and improving the robustness of a mixed-material layer model in the analysis of PFC erosion on present and future tokamak devices.

Carbon (C)↗
DOE OSTI2020

Development and demonstration of a methodology to evaluate high burnup fuel susceptibility to pulverization under a loss of coolant transient

For economic reasons, the US nuclear industry is renewing efforts to build a technical basis to extend rod average burnup limits above the current regulatory burnup limit of 62 GWd/MTU. The primary driver is to increase pressurized water reactor cycle lengths to 24 months, reducing the number of fresh fuel assemblies and core design constraints, thereby making core energy utilization more efficient. However, fuel pellet fragmentation and pulverization, termed high burnup fuel fragmentation (HBFF), has been observed in the high burnup (>90 GWd/MTU) Halden loss-of-coolant-accident (LOCA) integral test series. The issue gained attention when fuel fragmentation and pulverization were also observed closer to the current US regulatory limit during the US Nuclear Regulatory Commission (NRC) sponsored out-of-core integral test at Studsvik Nuclear in early 2011. This led to NRC concerns with potential changes to fuel and core designs relative to fuel pellet pulverization. In a letter to the NRC Commissioners, the staff specifically identified a need to “…define the boundary of safe operation for key fuel design and operating parameters,” stating that “the staff is challenged to evaluate the acceptability of future fuel design advancements and fuel utilization changes.” As such, it can be concluded that HBFF and potential dispersal into the reactor coolant system introduces additional complications in light-water reactor (LWR) fuel safety evaluations. However, it is not clear how much fuel will be susceptible to HBFF; nor has there been a methodology developed to evaluate fuel susceptibility to HBFF. To that end, this paper proposes an analysis methodology to assess fuel susceptibility to HBFF during LOCA scenarios. The work presented here uses the BISON fuel performance code to evaluate a representative pressurized water reactor fuel rod exposed to a rod average burnup of 75 GWd/MTU. Sensitivity studies investigated the impact of the peak cladding temperature, transient fission gas released, and pre-transient fission gas release on cladding ballooning and burst timing. Subsequently, a methodology to assess fuel susceptibility to HBFF will be developed based on experimental data published in the open literature. The methodology will then be demonstrated by calculating the mass of fuel susceptibility to HBFF. The BISON results conclude that increasing peak cladding temperature drastically decreased time to failure, and decreased balloon size both of which have been confirmed experimentally. Additionally, the effect of pre-transient and transient fission gas release affected cladding balloon size and burst timing. Finally, fuel susceptibility to HBFF significantly decreased as a function of peak cladding temperature.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

An improved pressure drop correlation for modeling localized effects in a pebble bed reactor

Advances in the development of pebble bed reactors (PBRs) has created a desire for accurate and cost-effective simulation tools for design scoping studies and safety analysis. The current state-of-the-art for these simulations is the use of porous media models, although these models rely on correlations to capture the effects of flow features that are not explicitly modeled. One of the areas where correlation accuracy is currently lacking is in the near-wall region of the bed. In this region, the presence of the wall causes the pebbles to pack more orderly, drastically changing the geometry and flow behavior in this region. This work presents a new generalized pressure drop correlation for PBRs based on the KTA equation. A high-to-low methodology is applied, where large eddy simulation (LES) is performed on two beds of 1568 and 1700 pebbles to generate a high-fidelity dataset. The flow fields are then averaged in time and separated into concentric rings of 0.05 D peb width. Average porosity, velocity, and pressure drop are extracted for each ring and the friction and form losses are calculated. The Reynolds number range for this study is 625–10,000, and thus the form losses are dominant over the friction losses. The form losses across the rings are investigated, and a correction term for the form loss calculation is determined and applied to the KTA equation to drastically improve the capability of modeling localized porosity effects in a porous media code. Finally, the improved correlation reduces near-wall velocity prediction error from over 50% with the KTA correlation to around 5%. Agreement in pressure drop prediction between LES and porous media simulations is also improved.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Reactivity-initiated accidents in two pressurized water reactor high burnup core designs

Here, this paper presents a safety analysis of two proposed core loadings for 24-month Pressurized Water Reactor (PWR) fuel cycles. This analysis focuses on reactivity-initiated accidents (RIAs) and evaluates core safety performance impacts of rod-averaged burnup limits up to 75 GWD/MTU and less than 7 % enriched UO 2 . The capabilities of Polaris, PARCS, and RELAP5-3D are leveraged to evaluate the core neutronic and thermal–hydraulic behavior for normal-operation, uncontrolled control rod withdrawal (CRW) transients, and control rod ejection (CRE) accidents. The two core designs are compared to identify features of realistic high burnup/extended enrichment core design approaches which have significant safety impact, identify experimental data needs for high-fidelity predictive modeling, and provide recommendations for future high burnup core designs. The first core design evaluated in this study was developed by Southern Nuclear Company and used an ZrB 2 Integral Fuel Burnable Absorber (IFBA) and B 4 C Wet-Annular Burnable Absorber (WABA)-based burnable poison strategy. The second core design assessed in this work used a Gd 2 O 3 -doped UO 2 burnable poison, similar to that used in boiling water reactors or French PWRs. Results indicate that fuel thermal limits are maintained for limiting CRW and hot full power (HFP) CRE transients. Cladding failure is predicted for the highest energy deposition rods in each core during limiting hot zero power (HZP) CRE accidents (where maximum radially averaged enthalpy exceeds 120 cal/g), though licensing may be permissible with a limited number of failed rods. While concerns exist regarding high critical boron concentration during steady state for the IFBA core and large plenum pressures for the gadolinia core design, the analysis demonstrates adequate safety performance during limiting RIA accident scenarios for two representative high burnup core designs. Design changes limiting plenum pressures and implementation of accident tolerant fuel (ATF) cladding features which minimize hydriding and susceptibility to pellet-cladding mechanical interaction (PCMI) are recommended for future high burnup fuel concepts. To support the technical basis for burnup limit increases, high-fidelity fuel performance models are needed to address physical effects not considered in this analysis, and high burnup irradiated fuel tests are required to extend applicability of the fuel failure limits and validate existing and future models.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Evaluation of the spatial self-shielding impact for TRISO-based nuclear fuel depletion

Reactor physics analyses of nuclear cores with nuclear fuel concepts containing tristructural isotropic (TRISO) particles, such as pebbles or compact fuel elements, rely on various degrees of simplification to keep these highly heterogeneous problems computationally tractable. One such limitation regards the level of spatial discretization employed during burnup calculations, where traditionally only a limited number of spatial zones are modeled at the full core level and assume that the spectrum is constant within the fuel elements and TRISO particles in this depletion zone. This type of assumption neglects the impact of spatial self-shielding effect within the kernels (microscale level) as well as within the compact or pebbles (mesoscale level). Furthermore, the Monte Carlo code Serpent 2 contains many relevant features for efficiently modeling this type of geometry, including a collision-based domain decomposition intended for very large burnup calculations, which we leveraged for this work to quantify the impact of capturing neutron flux variations occurring at the micro- and mesoscale level on a series of high-temperature gas-cooled reactor fuel element depletion problems. While spatial self-shielding is observed at both scales, with differences from a volume-averaged burnup of ±7% within the kernels and ±2% between TRISO particles within the fuel element, the conjugated effect on nuclide inventories and multiplication factor are negligible, hence confirming that assuming a single average spectrum value may be sufficient for most applications.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Evaluation of RANS vs. LES simulation of fluid flow through 3 × 3 rod bundle with a simple spacer grid as a precursor to coupled fluid–structure interaction simulations

The research literature on Computational Fluid Dynamics (CFD) of coolant flow through rod bundles with spacer-grids and mixing vanes is replete, ranging from high fidelity Large Eddy Simulation (LES)/Direct Numerical Simulation (DNS) simulations to Reynolds-Averaged Navier–Stokes (RANS) modeled studies. The mixing of flow between subchannels and the pressure drop through the bundle are fundamental quantities useful for comparing and evaluating CFD methods. Less commonly observed and compared are the forces exerted onto the structure by the fluid. The present study seeks to evaluate the use of RANS simulations for predicting the structural response to fluid flow. Wall resolved RANS simulations are benchmarked against LES simulations of fluid flow at a Reynolds number of 15,000 through a 3 × 3 fuel rod bundle with a simple spacer grid. Velocity line-plots are compared showing good agreement between RANS and LES results, ascertaining that the former is capable of capturing the essential time-averaged velocity profile. Additionally, the distribution of forces on the spacer grid and fuel rods are collected as a function of time and space. The RANS methods are evaluated using the frequency and magnitude of the fluctuating forces on various portions of the structure as compared to LES. In conclusion, the power spectral density evaluation of the models reveal underprediction of force amplitude on the rod walls by RANS and also discrepancy in the prediction of high frequency spectra, especially in the immediate vicinity of spacer-grid structure, which may be attributed to the lack of random turbulence fluctuation or insufficient modeling of small-scale eddies in RANS simulation.

FIV↗
DOE OSTI2026

Comparison of URANS and LES predictions for the open phase of the OECD NEA CSNI fluid structure interaction CFD benchmark

The OECD NEA CSNI WGAMA CFD Task Group ran a benchmark in 2020 and 2021 to assess the predictive capabilities of coupled fluid structure interaction (FSI) CFD analysis methods. This paper presents the predictions made for the open phase of the benchmark using URANS and LES turbulence modelling approaches, and a comparison of the results to the experimental data. The benchmark comprised a channel containing two inline cylinders in cross-flow. The cylinders were fixed at one end, free at the other, and had measured resonant frequencies and damping properties. The URANS modelling used ANSYS Fluent 2-way coupled to ANSYS Mechanical. The LES modelling used Nek5000, 1-way coupled to Diablo. Comparisons with cross-channel velocity profiles are presented, both for the mean flow and its RMS. Comparisons are also made to the frequency spectra for point measurements of fluid velocity and pressure, and for the accelerations of the free end of each cylinder. URANS predicts the average velocity profiles relatively well, and is able to predict the velocity and acceleration spectra at the shedding frequency. However, the frequency content at the 4th harmonic of the shedding frequency is low in the URANS flow fields, and so does not excite accelerations at the resonant frequency of the cylinders. LES makes better predictions of the average profiles, and the velocity spectra agree well at both the shedding frequency and at higher frequencies. In conclusion, the 1-way coupled LES results show good agreement for acceleration spectra.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Polysiloxane coatings on microspheres based on Multiphase Flow with Interface Exchange-Discrete Element Modelling

Polysiloxane coatings on yttria stabilized zirconia (YSZ) microspheres of 500 μm were simulated in a spout fluidized bed coater using Multiphase Flow with Interface Exchange-Discrete Element Modelling (MFiX-DEM). In this work, two different coater configurations were developed to study the influence of gas velocity and its distribution on particle dynamics. The presence of the Wurster tube not only enhances the distribution but also increases the overall residence time of the particles. Investigations were also carried out with different Wurster tube positions (normal, 10% and 20% lowered from its initial position). Among these, 20% lowered Wurster tube position demonstrated the most effective coating process. The effects of gas inlet pressure on the average gas velocity and the distribution of particles were analyzed. More than 97% of the particles can be retained. The derived results, including average gas velocity, particle retention percentage, and distribution of particles with gas velocity, are being used to guide the experimental work in obtaining defect-free coatings for YSZ microspheres.

36 MATERIALS SCIENCE↗