Search NASA⌕ Search

SEARCH · Search NASA

Results for “wall modeling”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Toward Verification of RANS Simulations of the T-Tube Modular Divertor Using Large Eddy Simulations of Impinging Turbulent Plane Jets

Turbulent impinging jets have been proposed to cool high heat flux plasma-facing components such as the solid tungsten target plates of the divertor in long-pulse magnetic fusion energy reactors. In particular, the T-tube modular divertor, originally developed by the ARIES Team, consists of two concentric cylindrical tubes where helium flows through a slot in the inner tube, forming an approximately planar jet that impinges upon and cools the inner surface of the pressure boundary (namely, the outer tube) and the ~15-cm 2 plasma-facing W target. The objective of this work is to demonstrate that large eddy simulations (LESs) accurately simulate the thermal transport in canonical flows that comprise the cooling flow in the T-tube, as well as validate temperatures from LES with experimental measurements in a simplified T-tube geometry. Wall‑resolved LESs, validated by experimental data and verified by direct numerical simulations (DNSs), provide benchmark data for two canonical flows in the T‑tube, namely, planar impinging and wall jets, for Reynolds numbers Re B = 4 × 10 3 to 2 × 10 4 . Our LES results are within 4% to 12% root-mean-square error (RMSE) of surface Nusselt number distributions (Nu) from experiments and DNSs. The validated LES results are then used as the ground truth to evaluate four Reynolds‑averaged Navier-Stokes (RANS) turbulence closures, namely, the k‑ω SST, realizable k‑ε, GEKO, and γ‑SST models. The k‑ω SST model has the best overall performance in terms of heat transfer, giving surface Nu within 12% RMSE of the LES results for high‑ReB impinging jets and reduced overprediction in the wall‑jet region. The GEKO model with default constants has the next best performance, providing slightly better Nu predictions for low ReB impinging jets (versus k-ω SST) but worse overall performance over the full range of ReB studied here. The realizable k‑ε turbulence model significantly overestimates turbulence near the stagnation point, while the γ‑SST model suppresses near‑wall production, biasing the simulations toward simulating laminar surface heat transfer. Simulations of the simplified T‑tube show that LES and RANS simulations with the k‑ω SST model give nearly identical average heat transfer coefficients (HTCs) over the impingement surface. The realizable k‑ε model predicts significantly lower wall temperatures due to overestimation of HTC in the outlet flow.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Domain-wall-mediated polarization switching in ferroelectric AlScN: Strain relief and field-dependent dynamics

While scandium-doped aluminum nitride (AlScN) exhibits robust ferroelectricity and excellent thermal stability, its utility is limited by an exceptionally high coercive field (𝐸 𝑐 ) for polarization switching. Unraveling the atomistic switching dynamics is therefore critical for tailoring 𝐸 𝑐 . Here, in this study, we combine density functional theory and machine-learning molecular dynamics to elucidate the polarization switching mechanisms in AlScN over various Sc concentrations and applied electric fields. We find that excessive lattice strain strictly prohibits collective polarization switching, but the preexisting domain walls relieve strain and lead to a distinct switching dynamics—dictating a field-dependent switching mechanism. At low electric fields, switching occurs via gradual domain-wall propagation consistent with the Kolmogorov-Avrami-Ishibashi model. In contrast, high fields stimulate additional nucleation, driving a rapid, homogeneous reversal process described by the simultaneous nonlinear nucleation and growth model. These findings highlight the critical role of domain-wall dynamics and suggest domain engineering as a viable strategy to tailor coercive fields in AlScN and related ferroelectrics.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Wall heating by subcritical energetic electrons generated by the runaway electron avalanche source *

Abstract Subcritical energetic electrons (SEEs) produced by the runaway electron (RE) avalanche source at energies below the runaway threshold are found to be the primary contributor to surface heating of plasma-facing components (PFCs) during final loss events. This finding is supported by theoretical analysis, computational modeling with the Kinetic Orbit Runaway electrons Code (KORC), and qualitative agreement with DIII-D experimental observations. The avalanche source generates significantly more secondary electrons below the runaway threshold, which thermalize rapidly when well-confined. However, during a final loss event, the RE beam impacts the first wall, and SEEs are deconfined before they can thermalize. Additionally, because the energy deposition length decreases faster than energy, the deposited energy density, and thus the maximum PFC surface temperature change, is larger for SEEs than REs. KORC simulations employ an analytic first wall to model particle deconfinement onto a non-axisymmetric wall composed of individual tiles. PFC surface heating is calculated using a 1D model extended to include an energy-dependent deposition length scale. Simulations of DIII-D qualitatively agree with infrared (IR) imaging only when SEEs from the avalanche source are included. These results demonstrate that SEEs are the dominant contributor to PFC surface heating and indicate that the avalanche source plays a critical role in the PFC damage caused during final loss events. The prominence of SEEs also has important implications for interpreting IR imaging, one of the primary diagnostics for RE-wall interaction diagnosis, despite REs dominating the energy and current density. This result improves predictions of wall damage due to post-disruption REs to estimate material lifetime and design RE mitigation systems for ITER and future reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Axion domain walls, small instantons, and non-invertible symmetry breaking

Non-invertible global symmetry often predicts degeneracy in axion potentials and carries important information about the global form of the gauge group. When these symmetries are spontaneously broken they can lead to the formation of stable axion domain wall networks which support topological degrees of freedom on their worldvolume. Such non-invertible symmetries can be broken by embedding into appropriate larger UV gauge groups where small instanton contributions lift the vacuum degeneracy, and provide a possible solution to the domain wall problem. We explain these ideas in simple illustrative examples and then apply them to the Standard Model, whose gauge algebra and matter content are consistent with several possible global structures. Each possible global structure leads to different selection rules on the axion couplings, and various UV completions of the Standard Model lead to more specific relations. As a proof of principle, we also present an example of a UV embedding of the Standard Model which can solve the axion domain wall problem. The formation and annihilation of the long-lived axion domain walls can lead to observables, such as gravitational wave signals. Observing such signals, in combination with the axion coupling measurements, can provide valuable insight into the global structure of the Standard Model, as well as its UV completion.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

A boundary value “reservoir problem” and boundary conditions for multi-moment multifluid simulations of sheaths

Multifluid simulations of plasma sheaths are increasingly used to model a wide variety of problems in plasma physics ranging from global magnetospheric flows around celestial bodies to plasma–wall interactions in thrusters and fusion devices. For multifluid problems, accurate boundary conditions to model an absorbing wall that resolves a classical sheath remain an open research area. This work justifies the use of vacuum boundary conditions for absorbing walls to show comparable accuracy between a multifluid sheath and lower moments of a continuum-kinetic sheath.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Recalculation of Soil Bulk Density Used in the Scorpius MCNP6 Model

The soil composition used in the Scorpius MCNP6.2 (Ref. 1) model was recently recalculated. Originally, it was “determined using the United States Geological Survey (USGS) Mercury Core Library and the Nevada National Security Site U.S. Geological Survey Databases (NNSS USGS) and Nevada National Security Site Petrographic, Geochemical, and Geophysical Database (NNSS PGG),” but the details had been lost. Reference 2 documents the updated soil composition determination. Reference 2 estimated the grain density of the soil instead of the bulk density. This report corrects that error. I am indebted to Garrett Euler (EES-17) for reading Ref. 2, pointing out its errors, and guiding me through the correct calculations. This report is organized as follows. For completeness, the composition calculations from Ref. 2 are repeated: Sec. II discusses the data that are reported in the PGG Access database, and Sec. III uses the data in the PGG Access database to compute the composition of the soil. Section IV estimates the bulk density of the soil. Section V evaluates the soil wall thickness in the Scorpius model. Section VI estimates the effect of the soil density change on previously calculated results. Section VII presents a summary and conclusions. Input files and output files are listed in appendices.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2023

Mesh and model requirements for capturing deep-stall aerodynamics in low-Mach-number flows

Here, this work presents a comprehensive computational fluid dynamics investigation of the effects of grid resolution and turbulence-model choice for capturing the unsteady three-dimensional aerodynamic performance of NACA 0012 and 0021 airfoils, with specific focus on the deep-stall regime. At high angles of attack (α), wind turbine blades routinely experience vortex-induced vibrations, which can cause significant structural damages. Accurate predictions of post-stall aerodynamics can identify the frequencies at which such vibrations maybe triggered. In this context, the NACA 0012 airfoil simulations are conducted at a chord-based Reynolds number, Re c =2×10 6 , with the k-ω Shear-Stress Transport Reynolds-Averaged Navier-Stokes (RANS) and Improved Delayed Detached Eddy Simulation (IDDES) hybrid RANS-Large Eddy Simulation turbulence models. The effect of mesh resolution both in the wall-normal and spanwise directions is investigated. Only the IDDES model with a minimum spanwise resolution of 24 cells per chord length correctly predicts the aerodynamic forces. Spectral analysis shows the peak primary shedding frequency at α=30°, which signifies the end of the stall region. In the post-stall regime, both lift and drag frequencies drop asymptotically with increasing α. The Strouhal number, based on normalised chord length, remains nearly constant in this region. Based on this study, NACA 0021 airfoil runs are performed with IDDES for Re c =2.7×10 5 and 2.0×10 6 on the finest wall-normal mesh and three spanwise grids. Simulations conducted on the finer spanwise grids demonstrate grid independence and show good agreement with experiments. The effect of varying Rec on the airfoil frequency statistics is investigated. Additionally, comparison studies are presented to investigate the impact of airfoil thickness on the frequency content at Re c =2.0×10 6 . The results from the study provide guidance on the choice of mesh resolution with the IDDES model to accurately capture aerodynamic quantities for complex industrial applications.

17 WIND ENERGY↗
DOE OSTI2020

Integrated modelling: Coupling of surface evolution and plasma-impurity transport

During the interaction of the scrape off layer (SOL) plasma with the first wall the evolution of both wall and plasma are tightly coupled: The erosion of the first wall leads to an impurity concentration in the plasma which affects the particle and power balance in the plasma. In turn the impurities, when leaving the plasma via transport, can form deposits and mixed materials, far away from their initial source location. These deposits can be eroded, allowing the impurities to stepwise migrate through the fusion device until they end up at a location where the plasma at the wall is cold enough and no further erosion occurs. To describe these processes an integrated model of surface evolution and plasma transport of impurities is needed. The WallDYN code achieves this required coupling of processes by parameterising the output of surface evolution- and plasma impurity-migration-codes by analytical models. For a given fixed background plasma it evolves the surface composition, derives impurity flux into and from the plasma and can from this derive the impurity densities in the plasma. This paper will show the importance of including this recycling of impurities at the wall in impurity migration modelling: The 13 CH 4 seeding experiment performed in the Wendelstein 7-X Stellerator is modelled using the recent extension of the WallDYN code to 3D plasma and wall geometries. A comparison with post mortem analysis of the 13 C deposition shows both qualitative and quantitative agreement with the WallDYN calculations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Analysis of Compressibility Effects and Nonlinear Property Variations in a Supercritical CO2 Mixing Layer

Calculations are performed using the large eddy simulation technique to qualitatively assess the effects of fluid compressibility and thermodynamic nonlinearities on the dynamics of a supercritical flow field. Resulting data are also used to interrogate instantaneous subfilter-scale fields, as modeled by the mixed-dynamic Smagorinsky closure and gradient diffusion model. A three-dimensional wall-resolved calculation of a spatially evolving mixing layer composed solely of carbon dioxide is performed. This is followed by a coarse grid calculation to examine model performance as a function of resolution. Theoretical analysis indicates that the partial derivatives of density with respect to pressure and temperature are important modulators of the temperature and pressure field evolution, respectively. Results pertaining to subfilter velocity-velocity and velocity-temperature fields indicate that the models are most stressed in the regions where compressibility and thermo-fluid nonlinearities dominate the physics.

42 ENGINEERING↗
DOE OSTI2025

A Model Intercomparison Study of Aerosol‐Cloud‐Turbulence Interactions in a Cloud Chamber: 1. Model Results

This study presents the first model intercomparison of aerosol‐cloud‐turbulence interactions in a controlled cloudy Rayleigh‐Bénard Convection chamber environment, utilizing the Pi Chamber at Michigan Technological University. We analyzed simulated cloud chamber‐averaged statistics of microphysics and thermodynamics in a warm‐phase, cloudy environment under steady‐state conditions at varying aerosol injection rates. Simulation results from seven distinct models (DNS, LES, and a 1D turbulence model) were compared. Our findings demonstrate that while all models qualitatively capture observed trends in droplet number concentration, mean radius, and droplet size distributions at both high and low aerosol injection rates, significant quantitative differences were observed. Notably, droplet number concentrations varied by over two orders of magnitude between models for the same injection rates, indicating sensitivities to the model treatments in droplet activation and removal and wall fluxes. Furthermore, inconsistencies in vertical relative humidity profiles and in achieving steady‐state liquid water content suggest the need for further investigation into the mechanisms driving these variations. Despite these discrepancies, the models generally reproduced consistent power‐law relationships between the microphysical variables. This model intercomparison underscores the importance of controlled cloud chamber experiments for validating and improving cloud microphysical parameterizations. Recommendations for future modeling studies are also highlighted, including constraining wall conditions and processes, investigating droplet/aerosol removal (including sidewall losses), and conducting simplified experiments to isolate specific processes contributing to model divergence and reduce model uncertainties.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Drift kinetic effects on plasma response to resonant magnetic perturbation for EU DEMO design

Abstract A systematic investigation of the plasma response to an externally applied resonant magnetic perturbation (RMP) field, for the purpose of controlling edge localized modes, is carried out for an EU DEMO reference plasma. Particular emphasis is placed on the role of kinetic effects associated with both thermal particles and fusion-born alphas. The single fluid, resistive model predicts a large peak amplification of the n = 1 ( n is the toroidal mode number) plasma response to the target equilibrium, which is found to be close to the Troyon no-wall limit. A more advanced response model, including kinetic resonances between the RMP perturbation and drift motions of thermal and energetic particles, on the other hand, finds a strong suppression of the n = 1 field amplification. A major role is played by the precessional drift resonance of fusion-born alphas. A strong parallel sound wave damping model is found to well reproduce the full kinetic response results for the DEMO plasma, in terms of both the resonant field response amplitude and the plasma displacement. Finally, both fluid and kinetic models produce similar responses for the n = 2 and 3 RMP fields for the considered DEMO plasma, whilst kinetic effects again become important for the n = 4 RMP due to proximity of the reference plasma to the no-wall limit for the n = 4 ideal kink instability.

Physics↗
DOE OSTI2020

A model for radiative heating of a high-Z pusher

Several approaches to Inertial Confinement Fusion (ICF), including double-shell, pushered-single-shell, and the Revolver designs, have fuel surrounded by pushers made from high-Z materials. An advantage of these designs is that radiation emitted by the hot fuel will be absorbed and re-radiated into the fuel to reduce cooling. This process is referred to as radiation-trapping, and it lowers the fuel temperature required for ignition. To elucidate the physics of radiation trapping, a simplified model comprising a set of coupled differential equations has been developed to model the arbitrary time-dependent wall temperature from the Hammer and Rosen solution [J. H. Hammer and M. D. Rosen, Phys. Plasmas 10, 1829 (2003)] to a Marshak wave. The derivation of the model and a set of analytical power-law solutions will be presented. A discussion on numerical implementation of the differential equations into a model for burn in ICF capsules is also included.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

A Simplified Numerical Approach to Characterize the Thermal Response of a Moving Bed Solar Reactor

Abstract Concentrated solar thermochemical storage in the form of a zero-emission fuel is a promising option to produce long-duration energy storage. Solar fuel is produced using a cavity reactor that captures concentrated solar radiation from a solar field of heliostats. In this paper, heat transfer model of a tubular plug-flow reactor designed and manufactured for a solar fuel production is presented. Experimental data collected from a fixed bed tubular reactor testing are used for model comparison. The system consists of an externally heated tube with counter-current flowing gas and moving solid particles as the heated media. The proposed model simulates the dynamic behavior of temperature profiles of the tube wall, gas, and particles under various gas flowrates and residence times. The heat transfer between gas–wall, solid particle–wall, and gas–solid particle is numerically studied. The model results are compared with the results of experiments done using a 4 kW furnace with a 150 mm heating zone surrounding a horizontal alumina tube (reactor) with 50.8 mm outer diameter and thickness of 3.175 mm. Solid fixed particles of magnesium manganese oxide (MgMn2O4) with the size of 1 mm are packed within the length of 250 mm at the center of the tube length. Simulation results are assessed with respect to fixed bed experimental data for four different gas flowrates, namely, 5, 10, 15, and 20 standard liters per minute of air, and furnace temperatures in the range of 200–1200 °C. The simulation results showed good agreement with maximum steady state error that is less than 6% of those obtained from the experiments for all runs. The proposed model can be implemented as a low-order physical model for the control of temperature inside plug-flow reactors for thermochemical energy storage applications.

Engineering↗
DOE OSTI2025

Simplified Model and Approach to Transform Infrared Surface Temperature to Film Effectiveness in a Conjugate Heat Transfer Experiment

This paper describes a simplified engineering model based on a one-dimensional thermal resistance network. The model is used to develop a new method to relate film cooling effectiveness and heat transfer augmentation to local overall cooling effectiveness in a conjugate flat plate experiment. This paper presents experimental proof-of-concept data to demonstrate the potential for this model. In contrast to previous approaches, neither the wall heat flux nor the adiabatic wall temperature is required to estimate the local film cooling performance parameters. The model predicts surface temperatures that are within the experimental uncertainties over the range for which the model is trained, and to within five percent when the model is extrapolated to higher coolant channel Reynolds numbers. This paper is relevant to conjugate test rigs that can measure the hot surface temperature distribution with and without film cooling. This information may also be relevant to designers as a method to approximate surface temperatures or used as an approximate heat transfer model for optimization studies.

advanced gas turbines↗
DOE OSTI2024

Flow reversal benchmark of a one-sided heated narrow rectangular channel with CATHARE and RELAP5

Flow reversal in narrow coolant channels can be a crucial phenomenon for the safety of research reactors with a downward nominal flow direction. During a loss of forced flow accident, the downward flow stagnates briefly before transitioning into an upward natural circulation flow. The fuel may be damaged if dryout occurs and threshold fuel and/or cladding temperatures are exceeded. A comprehensive study is provided for flow reversal in narrow rectangular channels by examining experimental data and conducting software model analyses. The literature on flow reversal was reviewed, and selected experimental datasets were used to benchmark against CATHARE and RELAP5 models and also compare the code calculations with each other. The experimental data comes from flow reversal tests conducted with a narrow rectangular channel with one-sided heating. The results were compared with experimental data for successful flow reversal tests and predicted dryout power for dryout conditions. Also, the study examined the effects of the pump coastdown period, inlet liquid temperature, system pressure, and localized pressure drops. The experimental results showed that shorter coastdown periods, reduced pressure drops, and lower coolant inlet temperatures increased the dryout power. However, the system pressure did not noticeably affect the results. The simulation results showed that both CATHARE and RELAP5 agreed with experimental data, capturing the trends of the experimental results. Slight differences between each code calculation, as well as the predicted and measured dryout powers, were attributed to experimental uncertainties and the modeling of physical phenomena such as wall nucleation, interfacial heat transfer, drag coefficients, and critical heat flux. Overall, this study provides an understanding of flow reversal and the prediction capabilities of thermal-hydraulics software models. In conclusion, a future study of the flow reversal benchmark of a narrow rectangular channel with two-sided heating may provide additional valuable insights.

CATHARE↗
DOE OSTI2021

Interlayer Interactions in 1D Van der Waals Moiré Superlattices.

Studying two-dimensional (2D) van der Waals (vdW) moiré superlattices and their interlayer interactions have received surging attention after recent discoveries of many new phases of matter that are highly tunable. Different atomistic registry between layers forming the inner and outer nanotubes can also form one-dimensional (1D) vdW moiré superlattices. In this review, experimental observations and theoretical perspectives related to interlayer interactions in 1D vdW moiré superlattices are summarized. The discussion focuses on double-walled carbon nanotubes (DWNTs), a model 1D vdW moiré system, and the authors highlight the new optical features emerging from the non-trivial strong interlayer coupling effect and the unique physics in 1D DWNTs. Future directions and questions in probing the intriguing physical phenomena in 1D vdW moiré superlattices such as, correlated physics in different 1D moiré systems beyond DWNTs are proposed and discussed.

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

Elucidating key reducing species beyond ions in hydrogen plasma smelting reduction of iron ore

Hydrogen plasma smelting reduction (HPSR) of iron ore has attracted significant attention over the past decade due to its high-temperature operation, rapid plasma mediated reduction kinetics, and simpler density-based separation of molten iron product, compared to H2-based solid-state reduction. All of these attributes enable processing of low-grade ores for downstream use in electric-arc furnaces, as virgin iron with low gangue content is required for high quality steel and improved furnace operation. While positive ions exist within the plasma arc, this work demonstrates that near the anodic ore surface, hydrogen radicals and vibrationally excited hydrogen species dominate and their densities correlate well with observed reduction rates. Species concentrations in the transferred plasma arc and at the plasma-ore interface are evaluated using coupled thermal plasma and near-wall non-equilibrium plasma models. The thermal plasma model is validated against experimental voltage data and spectroscopic measurements of plasma temperature and density for varying current inputs. Modeling of the near surface thermochemical non-equilibrium and micrometer scale anode sheath layer reveals, in addition to the expected H + , significant concentrations of ArH + and H$^+_3$ ions, typically not observed in thermal plasmas under thermodynamic equilibrium. Our results show that the inverted sheath structure at the anodic ore surface strongly suppresses reactive positive ion fluxes, while non-equilibrium electron-impact processes generate abundant hydrogen radicals and vibrationally excited species. These findings highlight the critical role of non-equilibrium effects in hydrogen arc-driven iron ore reduction and advance understanding beyond prevailing hypotheses centered on hydrogen ion-driven mechanisms.

08 HYDROGEN↗
DOE OSTI2020

Prediction of the thermal-hydraulic-mechanical response of a geological repository at large scale and sensitivity analyses

To better understand the Thermal-Hydraulic-Mechanical (THM) response of the rock caused by the high-level radioactive waste (HLW) released heat in the Callovo-Oxfordian formation (COx) in the near-field and far-field areas, a series of coupled THM modelling is performed. This study presents a case study based on the French geological repository for HLW (Cigéo project) to better assess modelling of deep geological repositories (DGR) within the DECOVALEX-2019 framework. In this study, the proposed coupled THM model, implemented in the platform of COMSOL, is validated by comparing the modelled THM results for a point heat source in an infinite rock mass with an analytical solution. It is also validated by comparing its thermal component for the Base Case with accurate results calculated using another numerical method. To shorten the calculation time, the DGR is modelled using six cuboid blocks representing 162 HLW cells and only six central HLW cells of interest modelled in detail. This simplification is validated. This model can easily be used to model not only the near-field THM response of the rock but also the far-field THM response. The influence of the boundary conditions applied on the gallery wall, HLW cell walls, and external surfaces of the model and the influence of different vertical dimensions of the geometry are studied. The sensitivity analyses of the THM parameters of COx on the THM response at different locations are performed. Finally, the influence of using a 2-dimensional (2D) model to represent a 3-dimensional (3D) repository is also investigated.

58 GEOSCIENCES↗