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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Thermal Hydraulics Analysis of a Divertor Monoblock Using SALAMANDER

Divertors are critical components in magnetic confinement fusion devices. One of the Divertor's crictial roles is absorbing the highest heat flux from the plasma. The Divertor accomplishes this through the use of divertor monoblocks and cooling channels. This work presents a thermal hydraulic analysis of a divertor monoblock using the Software for Advanced Large-scale Analysis of MAgnetic confinement for Numerical Design, Engineering and Research (SALAMANDER). SALAMANDER is an open source tool developed at Idaho National Laboratory for conducting multiphysics and multscale analysis of Plasma Facing Components (PFCs). For the monoblock, solid heat conduction is modeled in the block with a constant heat flux of 1E7 W/m^2 applied to the top of the block to represent peak heat flux from the plasma. The cooling channel (which is pressurized water at a Reynolds number of ~1,000,000) is modeled using a k-epsilon turbulence model with standard wall functions. The heat transfer between the solid and fluid domain is modeled using the Dittus-Boelter correlation for the convective heat transfer coefficient. This work aims to show the importance of a multiphysics modeling approach when performing simulations on PFCs.

70 - PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Challenging Common Assumptions of Thick-Wall Chamber Dynamics in Inertial Fusion Systems Using MOOSE-based Multiphysics Simulations

As an increasing number of companies look toward commercial Inertial Fusion Energy (IFE) designs, there is a pressing need to understand the physics of thick-wall chamber gas dynamics. The thick liquid wall approach implements a renewable wall to mitigate the fusion target emissions, thereby reducing the radiation damage rate and significantly extending the lifetime of chamber structures, leading to increased plant availability and reduced waste streams in comparison to dry wall chamber designs. It is necessary, however, to assess the critical performance and safety aspects of these systems. For example, it is crucial to predict (1) where the mass ablated from the liquid walls will vent, which determines the placement of condensing surfaces; (2) debris propagation up the beam lines, which provides essential information for design and protection requirements; (3) peak pressures and impulse on chamber walls, which affect chamber structural design; and (4) momentum transfer to the liquid jets, which constrains the shape and positioning of the jets. In turn, the chamber design and its liquid walls affect shielding requirements, material activation, and tritium fuel cycle. Currently available simulation tools, however, are unable to accurately capture key thick-wall chamber dynamics. Significant assumptions are often made to simplify the system and reduce computational cost and modeling capability needs, but the impact of these assumptions on simulation predictions has not been evaluated. For example, no three-dimensional simulations can be found in the open literature to evaluate gas venting and momentum transfer to the jets with simulations using two-dimensional domains to represent complex three-dimensional geometries. Moreover, limited studies have been dedicated to jet breakup due to both turbulence and neutron heating, and no studies have been found that evaluate how jet breakup can impact shock-jet interaction. Furthermore, effects of radiative heat transfer have rarely been included for the hydrodynamic phase of shock propagation, and integration of proper equations of state in shock dynamics codes has been mostly exploratory. In this study, we use the flexible, high-fidelity Multiphysics Object-Oriented Simulation Environment (MOOSE) to model these complex phenomena and inform design and safety studies. Capabilities to model thick-wall chamber gas dynamics are being developed, and the impact of the assumptions listed above (i.e., two-dimensional vs three-dimensional, absence of jet breakout, no radiative heat transfer, and ideal gas behavior) are being quantified.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Challenging Common Assumptions of Thick-Wall Chamber Dynamics in Inertial Fusion Systems using MOOSE

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Prevention of resistive wall tearing mode major disruptions with feedback

Resistive wall tearing modes (RWTMs) can cause major disruptions. A signature of RWTMs is that the rational surface is sufficiently close to the wall to interact with it. For (m,n) = (2,1) modes, a RWTM requires normalized minor radius of the rational surface ρ q2 ≥ 0.75, which can also be expressed as q 75 ≤ 2. Major disruptions can occur when the criterion is satisfied. This is confirmed in simulations and theory and in a DIII-D locked mode disruption database. The q 75 < 2 criterion is valid at high β as well as at low β. A very important feature of RWTMs is that they can be feedback stabilized. If the ρ q2 criterion is not satisfied, or if the wall is ideally conducting, then the mode does not produce a major disruption, although it can produce a minor disruption. Feedback, or rotation of the mode at the wall by complex feedback, can emulate an ideal wall, preventing major disruptions. The ρ q2 criterion depends weakly on the wall radius. A simple geometric model of its dependence on wall radius is given.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Simulation of Channel Flow with Square Ribs for Blanket First-Wall Cooling: Geometry-Specific Tuning of k-ω Model Using Adjoint Method

Cooling of the plasma-facing first wall is challenging in the design of blanket components because of the high heat flux (on the order of 𝑀𝑊/𝑚2) from the plasma, especially when a low thermal mass medium like helium is chosen as the coolant. Therefore, heat transfer enhancement in which the convective heat transfer rate is augmented by the addition of turbulence-promoting structures becomes a key initiative for providing sufficient cooling capability with helium. Previously, computational fluid dynamics simulations had been performed on pipe flows with different transverse and longitudinal ribbed geometries at Oak Ridge National Laboratory to compare the enhancement performance among different ribbed geometries. Rib shape morphing had been conducted to obtain an optimized rib profile. In the work presented here, the adjoint method is adopted in the ANSYS Fluent solver for turbulence model augmentation, and the Generalized k-ω (GEKO) turbulence model is employed because of its ability of tuning the turbulence model. The Nusselt number and pressure drop obtained from the channel flow with bottom ribbed wall experiments are used as the targets. Sensitivity analysis provides information as guidance to improve the turbulence model accuracy. The augmented GEKO model is tuned for the studied ribbed channel geometry and flow conditions, providing improved predictive accuracy within this context. Extension to other configurations offers potential but may require additional tuning and validation.

Xu, Tracy [ORNL] (ORCID:0009000193700887)↗

DOE OSTI2026

CFD Modeling of an 81" Single Slope Barrier Wall Drainage Window for on Grade and Sag Conditions

The main goal of this study is to design a drainage window in an 81" single slope barrier wall for on-grade and sag conditions that could replace the standard ODOT I-3 barrier catch basin. The ODOT design guidelines are as follows: a) Examine configurations of the opening to the flume behind the wall under the wall section to maximize conveyance of flow to a roadside ditch section. b) Examine window lengths and configurations to maximize hydraulic efficiency on a grade up to an absolute maximum window length of 20 ft. c) Use a minimum opening height of 4 in. and consider tapering to an increased height at the outlet to the flume to reduce the potential for debris clogging beneath the barrier section. On-grade locations are analyzed for longitudinal slope up to 5% and spread up to 10 ft. Sag locations, represented by pavement at zero longitudinal slope, are analyzed for water depths up to 18 in. The roadway cross section consists of a 10-ft shoulder at a 4.0% cross slope and pavement at a 1.6% cross slope. All surfaces are represented with Manning’s n = 0.015. Computational fluid dynamics (CFD) is used to perform the analysis. The modeling approach follows previous studies by the authors. The new design is similar to a slotted inlet: a rectangular opening cut into the bottom of a barrier (e.g., a temporary barrier or single slope barrier wall), which typically has a uniform cross section. Findings from laboratory tests of flow under Florida DOT barrier walls were used for additional validation of the modeling approach used in this study.

42 ENGINEERING↗

DOE OSTI2020

A simple model to scope out parameter space for indirect drive designs on NIF

In this work, we present a simple model to scope out parameter space for indirect-drive, inertial confinement fusion designs for the National Ignition Facility laser. Because the parameter space is large, simple models can be used to identify regions of parameter space for further study with more sophisticated models and experiments. We include a model for Hohlraum radiation drive and symmetry—both based on empirical scalings from the data. The model for radiation drive is based on assuming that the high atomic number (Z) Hohlraum wall dominates the energy balance during the high power, peak of the pulse (≳ 300 TW). We find that the time-dependent radiation drive flux can be described by the running integral of the laser energy divided by the Hohlraum area multiplied by constant slopes in two distinct time periods. The first period is when the laser power rises rapidly, so the radiation temperature increases due to changes in laser power and wall albedo. The second period is during peak power—here, the laser power is typically held constant—so, the radiation temperature increases only due to changes in the wall albedo. This model is applied to several NIF designs with different Hohlraum sizes, laser pulse length durations, and peak powers and energies. Drive and symmetry models can be combined to find regions of parameter space that have high capsule absorbed energy while maintaining a symmetric implosion. We propose a new metric for evaluating designs based on minimizing the radius at which the maximum implosion kinetic energy is achieved.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Earthquake-induced impact of base-isolated buildings: theory, numerical modeling, and design solutions

Earthquake shaking more intense than that used to size the horizontal clearance between a base–isolated building and near–rigid perimeter moat wall will result in hard impact, producing high–frequency, high–amplitude acceleration response in the structure and supported equipment. Herein this paper provides a design solution for the damaging effects of hard impact by installing a compliant engineered element in the load path between the base–isolated building and the moat wall, resulting in soft impact and a much smaller acceleration response. The engineered element assumed herein is a commercial–off–the–shelf marine fender with mechanical properties determined by physical testing. The attachment of a flexible engineered element, with well–defined stiffness and damping, to a near–rigid moat wall, simplifies the numerical modeling of the building–moat wall system and eliminates the need to bound the lateral stiffness of the wall for impact calculations. The simple model of the engineered element can be implemented in commercial finite element codes. Theory is developed for two–sided impact of a single–degree–of–freedom oscillator. Analytical solutions are derived for the shifted first–mode frequency of the impacted oscillator and for its free–vibration response. The shifted first–mode frequency is a function of the composite lateral stiffness of the isolator–engineered element assembly and its earthquake–induced displacement. Local peaks in the spectral response of the impacted oscillator form at odd integer multiples of the shifted first–mode frequency. The analytical solutions can be used to verify, in part, the numerical model used for impact analysis.

42 ENGINEERING↗

DOE OSTI2023

QCD-collapsed domain walls: QCD phase transition and gravitational wave spectroscopy

Abstract For a discrete symmetry that is anomalous under QCD, the domain walls produced in the early universe from its spontaneous breaking can naturally annihilate due to QCD instanton effects. The gravitational waves generated from wall annihilation have their amplitude and frequency determined by both the discrete symmetry breaking scale and the QCD scale. The evidence of stochastic gravitational waves at nanohertz observed by pulsar timing array experiments suggests that the discrete-symmetry-breaking scale is around 100 TeV, assuming the domain-wall explanation. The annihilation temperature is about 100 MeV, which could naturally be below the QCD phase transition temperature. We point out that the QCD phase transition within some domains with an effective large QCDθangle could be a first-order one. To derive the phase diagram inθand temperature, we adopt a phenomenological linear sigma model with three quark flavors. The domain-wall explanation for the NANOGrav, EPTA, PPTA and CPTA results hints at a first-order QCD phase transition, which predicts additional gravitational waves at higher frequencies. If the initial formation of domain walls is also a first-order process, this class of domain-wall models predicts an interesting gravitational wave spectroscopy with frequencies spanning more than ten orders of magnitude, from nanohertz to 100 Hz.

Physics↗

DOE OSTI2021

Model-based real-time surface heat flux and temperature estimation for the DIII-D tokamak

A control-oriented model for monitoring of wall power flux densities on the DIII-D tokamak has been successfully implemented and validated experimentally. Future reactors will have to withstand severe steady state high heat flux loads on plasma-facing components (PFCs). Due to the difficulty of directly-measuring local heat fluxes on these components, monitoring and protection of PFCs during the plasma discharge can benefit from simplified physics-based real-time (RT) functional models to estimate and guide heat load control. As a first step into the development, a control-oriented model for monitoring of wall power flux densities and temperatures on DIII-D tokamak has been successfully implemented. The paper discusses the experimental demonstration and comparison of the 2-D model-based wall heat flux algorithm on the DIII-Dinner wall limiter (IWL) against infra-red camera heat flux measurements for limited plasma configurations. The paper also reports on the benchmarking of the field line tracing environment, SMITER, developed at ITER organization on DIII-D tokamak against experimental IR diagnostic data and the derivation of the component shaping weighting factors for the 2-D model-based approach. Here, the extension of the model-based approach for surface temperature estimation on the DIII-D IWL is also presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Striations in moderate pressure dc driven nitrogen glow discharge

Plasma stratification has been studied for more than a century. Despite the many experimental studies reported on this topic, theoretical analyses and numerical modeling of this phenomenon have been mostly limited to rare gases. In this work, a one-dimensional fluid model with detailed kinetics of electrons and vibrationally excited molecules is employed to simulate moderate-pressure (i.e. a few Torrs) dc discharge in nitrogen in a 15.5 cm long tube of radius 0.55 cm. The model also considers ambipolar diffusion to account for the radial loss of ions and electrons to the wall. The proposed model predicts self-excited standing striations in nitrogen for a range of discharge currents. The impact of electron transport parameters and reaction rates obtained from a solution of local two-term and a multi-term Boltzmann equation on the predictions are assessed. In-depth kinetic analysis indicates that the striations result from the undulations in electron temperature caused due to the interaction between ionization and vibrational reactions. Furthermore, the vibrationally excited molecules associated with the lower energy levels are found to influence nitrogen plasma stratification and the striation pattern strongly. A balance between ionization processes and electron energy transport allows the formation of the observed standing striations. Simulations were conducted for a range of discharge current densities from ~0.018 to 0.080 mA cm –2 , for an operating pressure of 0.7 Torr. Parametric studies show that the striation length decreases with increasing discharge current. Furthermore, the predictions from the model are compared against experimental measurements and are found to agree favorably.

42 ENGINEERING↗

DOE OSTI2024

Mineral Deposition on the Rough Walls of a Fracture

Modeling carbonate growth in fractures and pores is important for understanding carbon sequestration in the environment or when supersaturated solutions are injected into rocks. Here, we study the simple but nontrivial problem of calcite growth on fractures with rough walls of the same mineral using kinetic Monte Carlo simulations of attachment and detachment of molecules and scaling approaches. First, we consider wedge-shaped fracture walls whose upper terraces are in the same low-energy planes and show that the valleys are slowly filled by the propagation of parallel monolayer steps in the wedge sides. The growth ceases when the walls reach these low-energy configurations so that a gap between the walls may not be filled. Second, we consider fracture walls with equally separated monolayer steps (vicinal surfaces with roughness below 1 nm) and show that growth by step propagation will eventually clog the fracture gap. In both cases, scaling approaches predict the times to attain the final configurations as a function of the initial geometry and the step-propagation velocity, which is set by the saturation index. The same reasoning applied to a random wall geometry shows that step propagation leads to lateral filling of surface valleys until the wall reaches the low-energy crystalline plane that has the smallest initial density of molecules. Thus, the final configurations of the fracture walls are much more sensitive to the crystallography than to the roughness or the local curvature. The framework developed here may be used to determine those configurations, the times to reach them, and the mass of deposited mineral. Effects of transport limitations are discussed when the fracture gap is significantly narrowed.

calcite↗

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 OSTI2025

Néel domain walls with bistable chirality in a perpendicularly magnetized ferrimagnetic insulator

Field-free spin-orbit torque-driven domain wall motion in magnetic thin films with perpendicular magnetic anisotropy (PMA) requires the domain walls to have Néel character. Conventionally, Néel domain walls are stabilized by the Dzyaloshinskii-Moriya interaction (DMI) in ultrathin films. Here, in a europium iron garnet thin film with PMA and an additional uniaxial in-plane anisotropy, we demonstrate two bistable Néel domain wall states in the absence of DMI, and the capability to toggle the wall states with an in-plane field pulse and consequently their directions of motion under a current pulse. We present a phase diagram for the bistable Néel domain wall states as a function of in-plane field pulse width and amplitude. By fitting the experimental data to an analytical model of Néel wall reversal through the nucleation and propagation of Bloch lines, we extract the length of the initial reversed domain wall segment and Bloch line nucleation energy barrier. Current-driven motion of in-plane anisotropy stabilized Néel walls is qualitatively different from that of DMI-stabilized ones owing to the different symmetry of the effective fields that stabilize the Néel configuration. Furthermore, we present a proof of principle demonstration for 2-bit random number generation based on the stochastic reversal of domain wall chirality. These results provide critical insight into the topological energy barrier of Bloch lines and identify paths towards domain wall-based memory and computing devices.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Discrete element simulation of Pebble Bed Reactors on graphics processing units

Prediction of pebble positions in a Pebble Bed Reactor (PBR) is necessary for both reactor physics and thermal hydraulics simulations as the arrangement of pebbles has a significant impact on the resulting core power, coolant flow, and fuel temperature. Knowledge of pebble movement as the fuel is cycled through the core is also critical for predicting the fuel residence time and subsequently, the fuel burnup. Simulation with the Discrete Element Method (DEM) can provide knowledge of both the fuel packing and the fuel movement during cycling. Previous works that have performed 3D full-core DEM simulation of PBRs have used simplified models that neglect reflector wall features. This work employs a graphics processing unit (GPU)-enabled DEM code, Project Chrono, to analyze the differences in pebble packing and pebble velocities between a simplified smooth PBR reflector and a more realistic reflector that includes circular wall features. Additionally, a sensitivity study is performed on the depth of the wall features to ensure that crystallization is prevented. Project Chrono is also validated for PBR cycling applications using experimental data. It is found that wall features with a depth of at least 0.5 pebble diameters significantly reduce crystallization in the near-wall region, leading to discrepancies in both packing fraction and pebble velocity in this region compared to the simplified reflector models. These discrepancies are found to lead to roughly a 5–10% difference in the prediction of the near-wall porosity and a 10% difference in the prediction of the velocity of pebbles near the wall. As a result of these discrepancies, it is suggested that future DEM simulations of PBRs include wall features to reduce modeling errors.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2022

Empower Wall: Active insulation system leveraging additive manufacturing and model predictive control

Buildings are one of the largest energy consumers worldwide, using large amounts of energy during their construction and for climate control during operation. Active insulation systems (AIS) have been shown to reduce the energy needed for climate control in buildings by dynamically regulating the heat transferred between a building’s interior and exterior. Infrastructure-scale additive manufacturing (AM) has the potential to reduce the resources needed for building construction. Combining these two technologies into a single building envelope would create a path towards more sustainable buildings. A test was conducted for the Federal Energy Management Program (FEMP) Energy Exchange training and trade show, in August 2021, to investigate a new building envelope design, termed the Empower Wall, that utilized an AIS and was constructed using AM. Model predictive control was implemented to manage operation of the Empower Wall in concert with the existing HVAC system. Finally, the prototype system demonstrated that the Empower Wall lowered total energy consumption and reduced the cost of energy used.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Internal energy balance and aerodynamic heating predictions for hypersonic turbulent boundary layers

The elemental equation governing heat transfer in aerodynamic flows is the internal energy equation. For a boundary layer flow, a double integration of the Reynolds-averaged form of this equation provides an expression of the wall heat flux in terms of the integrated effects, over the boundary layer, of various physical processes: turbulent dissipation, mean dissipation, turbulent heat flux, etc. Recently available direct numerical simulation data for a Mach 11 cold-wall turbulent boundary layer allows a comparison of the exact contributions of these terms in the energy equation to the wall heat flux with their counterparts modeled in the Reynolds-averaged Navier-Stokes (RANS) framework. Various approximations involved in RANS, both closure models as well as approximations involved in adapting incompressible RANS models to a compressible form, are assessed through examination of the internal energy balance. There are a number of potentially problematic assumptions and terms identified through this analysis. Here, the effect of compressibility corrections of the dilatational dissipation type is explored, as is the role of the modeled turbulent dissipation, in the context of wall heat flux predictions. The results indicate several potential avenues for RANS model improvement for hypersonic cold-wall boundary-layer flows.

42 ENGINEERING↗

DOE OSTI2025

High-Temperature Gas-Cooled Reactors Multiphysics Simulation Demonstration and Code Validation

This study presents a comprehensive benchmarking and verification effort of several thermal-hydraulic and multiphysics capabilities for high-temperature gas-cooled reactor applications. The first part of this effort focuses on the running-in verification of Griffin’s multiphysics capabilities, specifically for simulating the evolution of pebble-bed reactor cores from startup to equilibrium. Since Fiscal Year 2024, improvements and enhancements have been implemented in Griffin, including simplifying the process to specify streamlines and developing the online cross-section generation capability. In the absence of validation data, code-to-code comparisons are conducted with kugelpy, showing good agreement for integral quantities like k-eff predictions and predictions for maximum power density. However, accuracy issues are noted for more detailed quantities like the spatial distribution of fission rate densities which will require further work to address. The second part of this report presents an improved System Analysis Module (SAM) core channel model where the effects of cross flow are considered during the pressurized loss of forced cooling transient, resulting in an improved agreement of the predicted pebble temperature with respect to the predictions from the SAM 2D porous media model. Additionally, the wall channeling effect due to variable porosity at the near wall region of the core is also investigated. Furthermore, to demonstrate Griffin’s online cross-section generation capability, a Multiphysics simulation is performed by coupling Griffin to the SAM core channel model. In the third part of the report, as a part of the Organisation for Economic Co-operation and Development/Nuclear Energy Agency (OECD/NEA) thermal-hydraulic code validation benchmark activity for a high-temperature gas-cooled reactor, the High Temperature Test Facility (HTTF) is investigated first using the NekRS computational fluid dynamics (CFD) code to study the flow mixing phenomenon in the lower plenum of the facility. Then, code-to-code and code-to-data comparisons are performed for Test PG27, which is a pressurized conduction cooldown (PCC) test, using five different codes by six organizations from five countries. The different simulations show good agreements in terms of the general trend but there are differences in some results such as the peak temperatures of different regions and heat removal rate.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗