SEARCH · Engineering Papers

Results for “Wall Model”

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 73 records · Page 4

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Model Development for the ACME Lab Modular Wall and A-Frame Assemblies

The structural frame used for conducting large structural tests was replaced by the ACME (Applied Combined Mechanical Environments) Lab in Building 860 of Sandia National Laboratories in Albuquerque, NM. This new design required simulation and analysis to verify maximum loads and corresponding factors of safety. Two major assemblies which are part of the ACME Lab, the Modular Wall and A-Frame, were designed to operate with a multitude of components in various configurations and with varying loads. Simulations of the Modular Wall and A-Frame structures were performed through explicit analysis to determine safety factors for numerous loading scenarios. The Modular Wall had safety factors of 1.1 for the lateral, diagonal, and gusset loading conditions. The A-Frame had safety factors of 2 for all loading conditions simulated.

42 ENGINEERING↗

DOE OSTI2026

First Wall Design of a Tokamak Pilot Plant Using a Monte Carlo Model for 3-D Heat Flux Deposition

We present a method for calculating the heat fluxes deposited on nonaxisymmetric tokamak first wall components, allowing for a first-of-its-kind model for power handling in the tokamak far scrape-off layer (SOL). The DIV3D Monte Carlo model features strict global power conservation and can calculate the finite cross-field plasma transport into magnetically-shadowed regions, which is significant when dealing with meter-scale shadows introduced by components such as poloidal limiters or antennas. As a case study, we apply the DIV3D model to inform the distribution of first wall poloidal limiters in an ARC-class reactor device. We demonstrate that discrete protection limiters can efficiently reduce peak heat fluxes on recessed breeder wall components in the presence of significant far-SOL plasma fluxes. By varying the toroidal periodicity and radial standoff depth of the limiters, we demonstrate one of the tradeoffs that must be considered in first wall design: more limiters provide greater protection, but at the cost of reduced breeding performance. We also present the impact that radial misalignments between limiters would have on first wall power loading.

Monte Carlo methods↗

DOE OSTI2024

Validation of the ERO2.0 code using W7-X and JET experiments and predictions for ITER operation

Abstract The paper provides an overview of recent modelling of global material erosion and deposition in the fusion devices Wendelstein 7-X (W7-X), JET and ITER using the Monte-Carlo code ERO2.0. For validating the modelling tool in a three-dimensional environment, W7-X simulations are performed to describe carbon erosion from the graphite test divertor units, which were equipped in operational phase OP 1.2 and analysed post-mortem. Synthetic spectroscopy of carbon line emission is compared with experimental results from the divertor spectrometer measurement system, showing a good agreement in the e-folding lengths in the radial intensity profiles of carbon. In the case of metallic wall materials, earlier modelling of the Be/W environment in JET and ITER is revisited and extended with an updated set of sputtering and reflection data, as well as including the mixing model for describing the Be/W dynamics in the divertor. Motivated by recent H/D/T isotope experiments in JET, limited and diverted configuration pulses are modelled, showing the expected trend of both Be and W erosion increasing with isotope mass. For the JET diverted configuration pulses, it is shown that Be migrates predominantly to the upper part of the inner divertor where it initially leads to strong W erosion. With longer exposure time, the growth of a Be deposited layer leads to a reduction of W erosion in that region. A similar trend is observed in simulations of the ITER baseline Q = 10 scenario, however with a more symmetric Be migration pattern leading to deposition also on the outer divertor.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Computational Fluid Dynamics Analysis of the Molten Salt Tritium Transport Experiment Test Section

Tritium, a radionuclide produced through neutron capture by lithium and other elements (beryllium and fluoride) in molten salts, presents unique challenges to radionuclide release. This is true for both fusion energy breeder blankets and molten salt fission reactors. The fundamental understanding of tritium transport is crucial to the safe design and operation of these reactors. The Molten Salt Tritium Transport Experiment (MSTTE), currently under construction at Idaho National Laboratory, aims to investigate tritium transport phenomena using a forced-convection fluoride salt loop. This loop is designed to study various transport mechanisms, such as permeation through metals and gas-liquid interactions, and is intended to support future research on tritium extraction units. A critical aspect of the MSTTE loop design is ensuring a fully developed velocity profile before the fluid reaches the permeation test section where measurements are made. This study employs computational fluid dynamics to model the salt flow behavior within the MSTTE permeation test section. A realizable k-ε turbulent model with enhanced wall treatment is used to simulate the single-phase, vertical upward flow of molten salt FLiNaK under isothermal conditions. The simulation results indicated flow distortion and underdeveloped profiles at all planned flow rates within the test section due to the 85-deg sharp bend. To address this issue, a reduced diameter with a reducer and expander and a flow conditioner are investigated to achieve fully developed flow. The analysis showed that the flow conditioner successfully corrected the flow profile, achieving fully developed behavior at a flow rate of 50 liters per minute (LPM). In conclusion, this research enhances our understanding of flow dynamics in molten salt systems and contributes to optimizing tritium transport control technologies.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2026

Composite-dimensional topological codes with boundaries and defects

We introduce new algorithms and provide example constructions of stabilizer models for the gapped boundaries, domain walls, and 0D defects of Abelian composite-dimensional twisted quantum doubles. Using the physically intuitive concept of condensation, our algorithm explicitly describes how to construct the boundary and domain-wall stabilizers starting from the bulk model. This extends the utility of Pauli stabilizer models in describing nontranslationally invariant topological orders with gapped boundaries. To highlight this utility, we provide a series of examples, including a new family of quantum error-correcting codes where the double of ℤ4 is coupled to instances of the double semion (DS) phase. We discuss the codes' utility in the burgeoning area of quantum error correction with an emphasis on the interplay between deconfined anyons, logical operators, error rates, and decoding. We also augment our construction, built using algorithmic tools to describe the properties of explicit stabilizer layouts at the microscopic lattice level, with dimensional counting arguments and macroscopic-level constructions building on pants decompositions. The latter outlines how such codes' representation and design can be automated. Our results are validated by a series of error-correcting threshold calculations comparing our codes' performance with that of standard surface codes. To do so, we introduce a composite-dimensional belief-propagation decoder with ordered statistics that utilizes combination sweeps. Going beyond our worked-out examples, we expect our explicit step-by-step algorithms to pave the path for higher-dimensional codes to be discovered and implemented in near-future architectures that take advantage of various hardware platforms.

Mousa, Mohamad [Purdue University]↗

DOE OSTI2025

Outflow Boundary Conditions for Turbine-Integrated Rotating Detonation Combustors

This study examines outflow boundary conditions (BCs) in computational fluid dynamics (CFD) simulations of a transition duct with and without guide vanes that converts supersonic flow exiting a rotating detonation combustor (RDC) to subsonic flow to drive a turbine. Since the flow exiting the transition duct has swirling shock waves with significant spatial and temporal variations in pressure, temperature, and Mach number, imposing proper BCs poses a challenge. To ensure all swirling shock waves exit the transition duct without creating non-physical reflected waves at its outlet, this study examined three outflow BCs: (1) the average pressure imposed at the duct’s outlet, (2) a nonreflecting BC (NRBC) with a specified average pressure imposed at the duct’s outlet, (3) the average pressure imposed at the outlet of an extension duct made up of a buffer layer and a sponge layer. This study is based on the three-dimensional, unsteady density-weighted-ensemble-averaged continuity, Navier–Stokes, and energy equations for a thermally perfect gas closed by the realizable k–ε model and “enhanced” wall functions. The results obtained show that imposing an average pressure at the transition duct’s outlet produces spurious waves that degrade the physical meaningfulness of the solution. When the NRBC was applied, swirling shock waves exited the duct’s outlet without creating spurious waves. However, its usage requires the gas to be thermally, as well as calorically, perfect, which this study shows could be a concern. By imposing the average pressure at the outlet of an extension duct, the gas does not need to be calorically perfect. The results obtained show the effects of the sponge layer’s length and coarsening ratio on damping nonuniformities in non-physical reflected waves to ensure the flow exiting the transition duct’s outlet can do so as if there are no boundaries present and has the desired average pressure—even though the BC is applied at the extension duct’s outlet.

gas turbines↗

DOE OSTI2024

Advancing Thermal-Hydraulic Modeling Capabilities for Molten Salt Reactors in Pronghorn

This is an M2 report for the NEAMS program. This work focuses on advancing the modeling and simulation capabilities for Molten Salt Reactors (MSRs) within Pronghorn, addressing key thermal-hydraulic challenges that must be addressed for the safe and efficient operation of these advanced reactors. This report summarizes the development and validation of the following new models in Pronghorn: (i) two-equation RANS turbulence models, (ii) near-wall turbulence corrections due to curvature and roughness, and (iii) multidimensional two-phase drift-flux models. For each of these cases, along with the experimental validation, an application study of interest to MSRs is presented. In addition, this report also presents recent improvements in the domain-overlapping coupling between Pronghorn and SAM to address the cases of weakly compressible flow and flow reversal. Finally, an MSR action that greatly simplifies the setup of MSR models using NEAMS tools is described.

42 ENGINEERING↗

DOE OSTI2025

Proximity ferroelectricity in wurtzite heterostructures

Proximity ferroelectricity is an interface-associated phenomenon in electric-field-driven polarization reversal in a non-ferroelectric polar material induced by one or more adjacent ferroelectric materials. Here, in this paper, we report proximity ferroelectricity in wurtzite ferroelectric heterostructures. In the present case, the non-ferroelectric layers are AlN and ZnO, whereas the ferroelectric layers are Al 1−x B x N, Al 1−x Sc x N and Zn 1−x Mg x O. The layered structures include nitride–nitride, oxide–oxide and nitride–oxide stacks that feature two-layer (asymmetric) and three-layer (symmetric) configurations. Ferroelectric switching in both layers is validated by multimodal characterization methods, including polarization hysteresis, anisotropic chemical etching, second harmonic generation, piezo response force microscopy, electromechanical testing and atomic resolution polarization orientation imaging in real space by scanning transmission electron microscopy. We present a physical switching model in which antipolar nuclei originate in the ferroelectric layer and propagate towards the internal non-ferroelectric interface. The domain wall leading edge produces elastic and electric fields that extend beyond the interface at close proximity, reducing the switching barrier in the non-ferroelectric layer, and allowing complete domain propagation without breakdown. Density functional theory calculations of polymorph energies, reversal barriers and domain wall energies support this model. Proximity ferroelectricity enables polarization reversal in wurtzites without the chemical or structural disorder that accompanies elemental substitution, opening new questions and opportunities regarding interface-based ferroelectricity.

36 MATERIALS SCIENCE↗

DOE OSTI2024

A robust spectral element implementation of the $k - τ$ RANS model in Nek5000/NekRS

The $k - ω$ Reynolds Averaged Navier Stokes (RANS) model is one of the industry standard approaches for modeling of turbulent flows. It performs better than the $k - ϵ$ model for low Reynolds number flows and is also more suitable for boundary layers with adverse pressure gradients. Major drawback of the model, however, is that the asymptotic value of $ω$ at the walls is singular, necessitating the use of a contrived “sufficiently” large value for $ω$ as the boundary condition for its transport equation. Here, this invariably leads to the solution being sensitive to near wall grid spacing. While an acceptable solution for low order (finite volume) methods, the excessive near wall gradients lead to persistent numerical stability issues in high order codes. To alleviate the problem, specifically in the context of the high order spectral element code Nek5000, a regularized $k - ω$ approach was formulated in our prior work (Tomboulides et al., 2018). The formulation, however, relies on the use of wall distance and its gradients for modeling the closure terms and can pose problems for simulations in complex geometries. This work presents a novel implementation of the $k - τ$ RANS model in Nek5000, where $τ = 1/ω$, eliminating the need for regularization, owing to the asymptotically bounded behavior of the source terms in the $τ$ transport equation, and also eliminating dependence on wall distance. Robustness and stability of the $k - τ$ model is ensured through implicit treatment of the source terms and their careful numerical implementation and demonstrated through several cases aimed at verification and validation. Studies include both canonical and engineering relevant problems, viz., turbulent channel flow, pipe flow, backward facing step, flow over NACA 0012 airfoil and flow in a T-junction. Results from the $k - τ$ model are shown to be consistent with regularized $k - ω$ model and also with the $k - ω$ SST model in OpenFOAM (for select studies). Comparison with experimental data is also shown, where available, to bolster validation efforts for the $k - τ$ model implementation through prediction of key turbulent quantities of interest.

Nek5000↗

DOE OSTI2024

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

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

Kohno, H. (ORCID:0000000332690010)↗

DOE OSTI2026

Microstructure Validation of Graph Theory Model-Derived Cooling Rates in the Wire Arc Additive Manufacturing of ER70S-6 Steel

Wire arc additive manufacturing (WAAM) enables high-rate fabrication of large metallic components, but spatial variations in thermal history can lead to microstructural heterogeneity that requires efficient process models to evaluate. This study evaluates whether cooling rates extracted from a graph theory model (GTM)-based thermal simulation are consistent with the microstructural evolution observed in an ER70S-6 WAAM wall. Thermal histories from the model were analyzed at selected build heights, and cooling rates were extracted from the final thermal excursion through the austenite phase field. Microstructures at corresponding locations were characterized using electron backscatter diffraction (EBSD) to quantify grain size distributions, and pearlite interlamellar spacing was used as an additional indicator of cooling behavior. The modeled cooling rates were highest near the substrate and generally decreased with build height, consistent with the observed reduction in the fine grain fraction and the progressive shift in the grain size distribution as build height increased. Pearlite spacing trends also supported the modeled cooling rate variation. These results indicate that GTM-derived thermal histories can be post-processed into metallurgically meaningful cooling rate estimates for WAAM steel builds and linked to dataset specific empirical grain size distribution relationships for process–thermal history–microstructure assessment.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Thermally anisotropic building envelope for thermal management: finite element model calibration using field evaluation data

The thermally anisotropic building envelope (TABE) is an active building envelope that redistributes thermal loads in response to weather conditions and building energy demand. Conductive layers throughout the TABE distribute low-grade heat among hydronic loops, altering heat flow direction and intensity. Finite element models of TABE roof and wall panels were developed and calibrated using field evaluation data. The calibration results showed that heat flux differences between the experimental data and finite element models averaged –0.42% and 3.57%, with a maximum mean square error of 1.78 and 3.96 for roof and wall panels, respectively. A reduction in heat flux from the environment to the building living space over the entire testing period (weeks in July/August) was found to be 85% for roof panels and 335% (load reversed) for wall panels. Finally, these results indicate TABE can effectively harness low-grade thermal energy sources to achieve high energy efficiency and promote demand-side management.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2026

Towards fully predictive gyrokinetic full- f simulations: validation and triangularity studies in TCV

Designing economical magnetic confinement fusion power plants motivates computational tools that can estimate plasma behavior from engineering parameters without direct reliance on experimental measurement of the plasma profiles. In this work, we present full-f global long-wavelength gyrokinetic simulations of edge and scrape-off layer turbulence in tokamaks that use only magnetic geometry, heating power, and particle inventory as inputs. Unlike many modeling approaches that employ free parameters fitted to experimental data, raising uncertainties when extrapolating to reactor scales. This approach directly simulates turbulence and resulting profiles through gyrokinetics without such empirical adjustments. This is achieved via an adaptive sourcing algorithm in Gkeyll that strictly controls energy injection and emulates particle sourcing due to neutral recycling. We show that the simulated kinetic profiles compare reasonably well with Thomson scattering and Langmuir probe data for Tokamak á Configuration Variable (TCV) discharge #65125, and that the simulations reproduce characteristic features such as blob transport and self-organized electric fields. Applying the same framework to study triangularity effects suggests mechanisms contributing to the improved confinement reported for negative triangularity (NT). Simulations of TCV discharges #65125 and #65130 indicate that NT increases the E x B flow shear (by about 20% in these cases), which correlates with reduced turbulent losses and a modest change in the distribution of power exhaust to the vessel wall. While the physical models contain approximations that can be refined in future work, the predictive capability demonstrated here, evolving multiple profile relaxation times with kinetic electron and ion models in hundreds of GPU hours, indicates the feasibility of using Gkeyll to support design studies of fusion devices.

Hoffmann, Antoine Cyril David [Princeton Plasma Ph↗

DOE OSTI2025

Mechanical Roles of Polysaccharide Assembly and Interactions in Plant Cell Walls

Plants synthesize polysaccharide-based primary cell walls that possess unique microstructures and mechanical properties to accommodate plant growth and provide protection. Here, it remains challenging to assess the role of polysaccharide organization and interactions in the mechanical behavior of primary cell walls owing to their complex microstructure and highly nonlinear mechanical responses. Employing a coarse-grained molecular dynamics model developed for onion epidermal walls, this work explores the conditions under which polysaccharide assembly and interactions might play a significant role in primary cell wall mechanics. Cellulose–cellulose adhesion plays a dominant role in the wall load-bearing capacity, but when cellulose–cellulose adhesion was disrupted computationally, cellulose–xyloglucan adhesion could influence the wall load-bearing capacity. Contrary to the common concept that xyloglucans mechanically tether well-separated cellulose microfibrils, xyloglucans functioned in this case as interfibrillar adhesives capable of transmitting tensile forces between cellulose microfibrils. Our findings may inform design criteria of new materials inspired by plant cell walls.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2025

Integrative path modeling and QTL mapping identify maturity, stem strength, and cell wall composition driving lettuce resistance to Sclerotinia minor

Lettuce ( Lactuca sativa ) is highly vulnerable to Sclerotinia minor , the pathogen causing lettuce drop. Breeding for resistance is the most effective control strategy; however, full resistance has not been achieved, and current partial resistance sources are often linked with undesirable traits, such as early bolting. This study aimed to unravel the genetic basis of partial resistance to S. minor and its relationship with plant maturity (bolting), stem mechanical strength (SMS), and cell wall composition (CWC) using a recombinant inbred line (RIL) population derived from a cross between the susceptible iceberg cv. ‘Salinas’ and the resistant oil-seed accession PI 251246. Field evaluations indicated that resistance was linked to earlier bolting, stronger stems, and higher pentose content. Path analysis demonstrated that earlier-maturing plants exhibited increased resistance through enhanced SMS and modified CWC, particularly with higher xylose and lower arabinose levels. Further analysis indicated a significant relationship between syringyl lignin content and resistance, especially in plants with varying bolting responses. Three key quantitative trait loci (QTLs) on linkage groups (LG) 2, 6, and 7 were consistently associated with resistance, bolting, and SMS. Importantly, residual QTL analysis revealed that the resistance locus on LG7 acted independently of maturity, suggesting a distinct resistance mechanism. Callose synthase emerged as a key candidate gene within the LG7 resistance QTL, located near - but distinct from - genes associated with plant maturity and flowering. These findings provide valuable insights into decoupling resistance from early bolting, suggesting a pathway for breeding lettuce cultivars with improved disease resistance and delayed bolting.

Lactuca↗

DOE OSTI2025

Hybrid SO(10) Axion Model without Quality Problem

Invisible axion models that solve the strong C P problem via the Peccei-Quinn (PQ) mechanism typically have a quality problem that arises from quantum gravity effects which violate all global symmetries. These models therefore require extreme fine-tuning of parameters for consistency. We present a new solution to the quality problem in a unified SO ( 10 ) × U ( 1 ) a gauge model, where U ( 1 ) a is an anomaly free axial gauge symmetry. PQ symmetry emerges as an accidental symmetry in this setup, which admits a PQ breaking scale as large as 7 × 10 11 GeV , allowing for the axion to be the cosmological dark matter. We call this a hybrid axion model due to its unique feature that it interpolates between the popular KSVZ and DFSZ axion models. Its predictions for the experimentally measurable axion couplings to the nucleon and electron are distinct from those of the usual models, a feature that can be used to test it. Furthermore, the model has no domain wall problem and it provides a realistic and predictive framework for fermion masses and mixings. Published by the American Physical Society 2025

Babu, K. S. (ORCID:0000000161475155)↗

DOE OSTI2026

Analytic model for neutral penetration and plasma fueling

Neutral atoms recycled from wall interaction interact with confined plasma, thereby refuelling it, most strongly in the region closest to the wall. This occurs near the X-point in diverted configurations, or else near the wall itself in limited configurations. A progression of analytic models is developed for neutral density in the vicinity of a planar or linear source in an ionising domain. First-principles neutral transport simulations with DEGAS2 are used throughout to test the validity and limits of the model when using equivalent sources. The model is further generalised for strong plasma gradients or the inclusion of charge exchange. An important part of the problem of neutral fuelling from recycling is thereby isolated and solved with a closed-form analytic model. A key finding is that charge exchange with the confined plasma can be significantly simplified with a reasonable sacrifice of accuracy by treating it as a loss. The several assumptions inherent to the model (and the simulations with which it is compared) can be adapted according to the particular behaviour of neutrals in the divertor and the manner in which they cross the separatrix.

fusion plasma↗

DOE OSTI2024

Role of Cell Wall Polysaccharides in Water Distribution During Seed Imbibition of Hymenaea courbaril L.

Seed water imbibition is critical to seedling establishment in tropical forests. The seeds of the neotropical tree Hymenaea courbaril have no oil reserves and have been used as a model to study storage cell wall polysaccharide (xyloglucan - XyG) mobilization.We studied pathways of water imbibition in Hymenaea seeds. To understand seed features, we performed carbohydrate analysis and scanning electron microscopy. We found that the seed coat comprises a palisade of lignified cells, below which are several cell layers with cell walls rich in pectin. The cotyledons are composed mainly of storage XyG. From a single point of scarification on the seed surface, we followed water imbibition pathways in the entire seed using fluorescent dye and NMRi spectroscopy. We constructed composites of cellulose with Hymenaea pectin or XyG. In vitro experiments demonstrated cell wall polymer capacity to imbibe water, with XyG imbibition much slower than the pectin-rich layer of the seed coat.We found that water rapidly crosses the lignified layer and reaches the pectin-rich palisade layer so that water rapidly surrounds the whole seed. Water travels very slowly in cotyledons (most of the seed mass) because it is imbibed in the XyG-rich storage walls. However, there are channels among the cotyledon cells through which water travels rapidly, so the primary cell walls containing pectins will retain water around each storage cell.The different seed tissue dynamic interactions between water and wall polysaccharides (pectins and XyG) are essential to determining water distribution and preparing the seed for germination.

arabinoxylan↗