SEARCH · Engineering Papers

Results for “variable geometry”

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 379 records · Page 21

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Refractory-based thermal energy storage for industrial process heat: one-dimensional modeling, control, and optimization

The variable and weather-dependent output of wind and solar power plants present a substantial challenge for planning and operating electricity-systems, particularly in the absence of cost-effective and dispatchable energy storage technologies. This study investigates a high-temperature, electrically heated, refractory-based thermal energy storage (RTES) system that stores electrical energy as sensible heat in dense ceramic bricks over the 950–1800 °C range. The stored heat can be discharged as a controlled hot-gas stream for industrial heating, fuel substitution in high-temperature processes, or electricity generation. The main novelty is a comprehensive modelling, control, mapping, and optimization framework that integrates one-dimensional transient gas–solid heat transfer, fan-assisted discharge, bypass-flow regulation, reheating logic, fan-power evaluation, insulation-loss assessment, and genetic-algorithm-based design optimization. The model uses feedback from outlet temperature and delivered power to regulate discharge, while a two-stage genetic algorithm optimizes brick-channel geometry, gas-flow operation, and multilayer insulation thicknesses. Storage capacities below 50 MWh and discharge powers of 5–30 MW are analyzed to evaluate hold time, thermal delivery, fan-power penalty, heat loss, state-of-charge evolution, and indicative capital cost. Results demonstrate that optimized and well-insulated refractory-based thermal energy storage units can provide stable, efficient, and repeatable heat delivery over multiple discharge cycles. The generated performance and cost maps support modular refractory thermal energy storage as a practical option for large-scale integration of wind and solar generation and for high-temperature industrial process heat.

25 ENERGY STORAGE↗

DOE OSTI2020

Single-ion magnetism in the extended solid-state: insights from X-ray absorption and emission spectroscopy

Large single-ion magnetic anisotropy is observed in lithium nitride doped with iron. The iron sites are two-coordinate, putting iron doped lithium nitride amongst a growing number of two coordinate transition metal single-ion magnets (SIMs). Uniquely, the relaxation times to magnetisation reversal are over two orders of magnitude longer in iron doped lithium nitride than other 3d-metal SIMs, and comparable with high-performance lanthanide-based SIMs. To understand the origin of these enhanced magnetic properties a detailed characterisation of electronic structure is presented. Access to dopant electronic structure calls for atomic specific techniques, hence a combination of detailed single-crystal X-ray absorption and emission spectroscopies are applied. Together K-edge, L 2,3 -edge and Kβ X-ray spectroscopies probe local geometry and electronic structure, identifying iron doped lithium nitride to be a prototype, solid-state SIM, clean of stoichiometric vacancies where Fe lattice sites are geometrically equivalent. Extended X-ray absorption fine structure and angular dependent single-crystal X-ray absorption near edge spectroscopy measurements determine Fe I dopant ions to be linearly coordinated, occupying a D 6h symmetry pocket. The dopant engages in strong 3dπ-bonding, resulting in an exceptionally short Fe–N bond length (1.873(7) Å) and rigorous linearity. It is proposed that this structure protects dopant sites from Renner–Teller vibronic coupling and pseudo Jahn–Teller distortions, enhancing magnetic properties with respect to molecular-based linear complexes. The Fe ligand field is quantified by L 2,3 -edge XAS from which the energy reduction of 3d z 2 due to strong 4s mixing is deduced. Quantification of magnetic anisotropy barriers in low concentration dopant sites is inhibited by many established methods, including far-infrared and neutron scattering. We deduce variable temperature L 3 -edge XAS can be applied to quantify the J = 7/2 magnetic anisotropy barrier, 34.80 meV (~280 cm –1 ), that corresponds with Orbach relaxation via the first excited, M J = ±5/2 doublet. The results demonstrate that dopant sites within solid-state host lattices could offer a viable alternative to rare-earth bulk magnets and high-performance SIMs, where the host matrix can be tailored to impose high symmetry and control lattice induced relaxation effects.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Holographic Space-Time and Quantum Information

The formalism of Holographic Space-time (HST) is a translation of the principles of Lorentzian geometry into the language of quantum information. Intervals along time-like trajectories, and their associated causal diamonds, completely characterize a Lorentzian geometry. The Bekenstein-Hawking-Gibbons-’t Hooft-Jacobson-Fischler-Susskind-Bousso Covariant Entropy Principle, equates the logarithm of the dimension of the Hilbert space associated with a diamond to one quarter of the area of the diamond’s holographic screen, measured in Planck units. The most convincing argument for this principle is Jacobson’s derivation of Einstein’s equations as the hydrodynamic expression of this entropy law. In that context, the null energy condition (NEC) is seen to be the analog of the local law of entropy increase. The quantum version of Einstein’s relativity principle is a set of constraints on the mutual quantum information shared by causal diamonds along different time-like trajectories. The implementation of this constraint for trajectories in relative motion is the greatest unsolved problem in HST. The other key feature of HST is its claim that, for non-negative cosmological constant or causal diamonds much smaller than the asymptotic radius of curvature for negative c.c., the degrees of freedom localized in the bulk of a diamond are constrained states of variables defined on the holographic screen. This principle gives a simple explanation of otherwise puzzling features of BH entropy formulae, and resolves the firewall problem for black holes in Minkowski space. It motivates a covariant version of the CKN [1] bound on the regime of validity of quantum field theory (QFT) and a detailed picture of the way in which QFT emerges as an approximation to the exact theory

tensor network↗

DOE OSTI2020

Koopman-Based Approach to Nonintrusive Reduced Order Modeling: Application to Aerodynamic Shape Optimization and Uncertainty Propagation

A methodology for nonintrusive projection-based nonlinear model reduction originally presented by Renganathan et al. is further extended here toward parametric systems with a focus on application to aerospace design. Specifically, the method is extended to address static systems with parametric geometry (that deforms the mesh) in addition to parametric freestream boundary conditions. The main idea is to first perform a transformation on the governing equations such that it is lifted to a higher-dimensional but linear underdetermined system. This enables one to extract the system matrices easily as compared to that of the original nonlinear system. The underdetermined system is closed with a set of model-dependent nonlinear constraints upon which the model reduction is finally performed. The methodology is validated on the subsonic and transonic inviscid flows past the NACA0012 and the RAE2822 airfoils with parametrized shapes. The utility of the approach is further demonstrated by applying it to two common problems in aerospace design, namely, derivative-free global optimization and parametric uncertainty quantification with Monte Carlo sampling. Overall, the methodology is shown to achieve accuracy up to 5% and computational speedup of two to three orders of magnitude relative to the full order model. Comparison against another nonintrusive model reduction method revealed that the proposed approach is more robust, accurate, and retains the consistency between the state variables.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Versailles Project on Advanced Materials and Standards interlaboratory study on intensity calibration for x-ray photoelectron spectroscopy instruments using low-density polyethylene

We report the results of a Versailles Project on Advanced Materials and Standards interlaboratory study on the intensity scale calibration of x-ray photoelectron spectrometers using low-density polyethylene (LDPE) as an alternative material to gold, silver, and copper. An improved set of LDPE reference spectra, corrected for different instrument geometries using a quartz-monochromated Al Kα x-ray source, was developed using data provided by participants in this study. Using these new reference spectra, a transmission function was calculated for each dataset that participants provided. When compared to a similar calibration procedure using the NPL reference spectra for gold, the LDPE intensity calibration method achieves an absolute offset of ~3.0% and a systematic deviation of ±6.5% on average across all participants. For spectra recorded at high pass energies (≥90 eV), values of absolute offset and systematic deviation are ~5.8% and ±5.7%, respectively, whereas for spectra collected at lower pass energies (<90 eV), values of absolute offset and systematic deviation are ~4.9% and ±8.8%, respectively; low pass energy spectra perform worse than the global average, in terms of systematic deviations, due to diminished count rates and signal-to-noise ratio. Differences in absolute offset are attributed to the surface roughness of the LDPE induced by sample preparation. We further assess the usability of LDPE as a secondary reference material and comment on its performance in the presence of issues such as variable dark noise, x-ray warm up times, inaccuracy at low count rates, and underlying spectrometer problems. In response to participant feedback and the results of the study, we provide an updated LDPE intensity calibration protocol to address the issues highlighted in the interlaboratory study. Finally, we also comment on the lack of implementation of a consistent and traceable intensity calibration method across the community of x-ray photoelectron spectroscopy (XPS) users and, therefore, propose a route to achieving this with the assistance of instrument manufacturers, metrology laboratories, and experts leading to an international standard for XPS intensity scale calibration.

36 MATERIALS SCIENCE↗

DOE OSTI2024

An explainable variational autoencoder model for three-dimensional acoustic emission source localization in hollow cylindrical structures

We introduce an explainable variational autoencoder for three-dimensional (3D) localization of acoustic emission sources in hollow cylindrical structures, with an unsupervised approach. This research capitalizes on multi-arrival waveforms generated by helical path propagation in cylindrical geometries to enable efficient two-receiver localization. By integrating the modal characteristics of Lamb modes under multi-path conditions, we demonstrate that two sets of time-of-arrival differences and peak amplitudes extracted from one receiver can serve as effective localization features. This initial approach identifies four potential source locations, highlighting the feasibility of two-receiver source localization using traditional feature extraction methods. However, direct extraction can be challenging when mode overlaps occur, complicating the localization process. To address this, our work proposes a novel waveform-based method. This method leverages the consistent dispersion characteristics within isotropic materials, where each unique combination of mode arrival times and peak amplitudes constructs a distinct waveform. This distinctiveness overcomes the ambiguities associated with mode overlaps, significantly enhancing the method’s precision and robustness. Our approach adopts a data-driven strategy for waveform-based localization using variational autoencoder (VAE). VAE discerns waveform patterns for localization, while also addressing data uncertainties. The VAE’s encoder and decoder networks capture the localization process and the source’s influence on waveform generation, respectively, guiding latent variables to segregate waveforms by source in the latent space. The design of the learning process focuses on specific localization characteristics to enhance result explainability. Localization predictions are generated by projecting test waveforms, not included in the training set, onto a trained latent space. The prediction is determined using a nearest-neighbor approach based on the closest latent representation of a source. Validation with pencil-lead-break tests on a metallic pipe confirmed our method’s effectiveness, achieving an averaged 3D localization accuracy of 0.84.

Lee, Guan-Wei↗

DOE OSTI2026

Measurements of Quasar Proximity Zones with the Lyα Forest of DESI Y1 Quasars

Abstract The intergalactic medium (IGM) around a quasar is shaped by its dense environment and by its excess ionizing radiation, which form a “quasar proximity zone” whose size and anisotropy depend on the quasar’s halo mass, luminosity, age, and radiation geometry. Using over 10,000 quasar pairs from the Dark Energy Spectroscopic Instrument (DESI) Year 1 data, with projected comoving separations r ⊥ < 2 h −1 Mpc, we investigate how the proximity zone of foreground quasars at z ∼ 2–3.5 affects Lyα absorption in their background quasars. The large DESI sample enables unprecedented precision in measuring this “transverse proximity” effect, allowing a detailed investigation of the signal’s dependence on the projected separation of quasar pairs and the luminosity of the foreground quasar. We find that enhanced gas clustering near quasars dominates over their ionizing effect, leading to stronger absorption on neighboring sightlines. Under the assumption that quasar ionizing luminosity is isotropic and steady, we infer the IGM overdensity profile in the vicinity of quasars, finding overdensities as high as Δ ∼ 10 at comoving distance ∼1 h −1 Mpc from the most luminous systems. Surprisingly, however, we find no significant dependence of the proximity profile on the luminosity of the foreground quasar. This lack of luminosity dependence could reflect a cancellation between higher ionizing flux and higher gas overdensity, or it could indicate that quasar emission is highly time-variable or anisotropic, so that the observed luminosity does not trace the ionizing flux on nearby sightlines.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Data-Efficient Dimensionality Reduction and Surrogate Modeling of High-Dimensional Stress Fields

Tensor datatypes representing field variables like stress, displacement, velocity, etc., have increasingly become a common occurrence in data-driven modeling and analysis of simulations. Numerous methods [such as convolutional neural networks (CNNs)] exist to address the meta-modeling of field data from simulations. As the complexity of the simulation increases, so does the cost of acquisition, leading to limited data scenarios. Modeling of tensor datatypes under limited data scenarios remains a hindrance for engineering applications. Here, in this article, we introduce a direct image-to-image modeling framework of convolutional autoencoders enhanced by information bottleneck loss function to tackle the tensor data types with limited data. The information bottleneck method penalizes the nuisance information in the latent space while maximizing relevant information making it robust for limited data scenarios. The entire neural network framework is further combined with robust hyperparameter optimization. We perform numerical studies to compare the predictive performance of the proposed method with a dimensionality reduction-based surrogate modeling framework on a representative linear elastic ellipsoidal void problem with uniaxial loading. The data structure focuses on the low-data regime (fewer than 100 data points) and includes the parameterized geometry of the ellipsoidal void as the input and the predicted stress field as the output. The results of the numerical studies show that the information bottleneck approach yields improved overall accuracy and more precise prediction of the extremes of the stress field. Additionally, an in-depth analysis is carried out to elucidate the information compression behavior of the proposed framework.

artificial intelligence↗

DOE OSTI2021

Floodplain Wetland Channel Planform, Cross-Sectional Morphology, and Sediment Characteristics Along an Estuarine to Tidal River Gradient

The formation of wetland channel networks on coastal river floodplains is affected by spatial transitions in tidal-fluvial processes. This study evaluated sedimentation patterns in tidal marshes on the 147-km tidal Columbia River and its 87-km estuary and characterized the cross-sectional geometry and planform morphometry of tidal wetlands along the estuarine to riverine gradient. Tidal marshes were predominantly depositional (median 0.7 cm yr-1, primarily silt), consistent with late-Holocene rates. Marsh sediment accretion rate increased with annual maximum 7-day average daily maximum water depth. Elevation was negatively correlated with accretion rate in the tidal river and positively correlated with total organic carbon (TOC%) in the tidal river and estuary. TOC% was greater in the estuary and, like percent fines, greater farther from the main-channel Columbia River. For four single-channel and six channel-network planforms of tidal marsh, shrub-dominated, and forested wetlands: (1) Spatial patterns of cross-sectional morphology conformed to typical coastal wetland morphology in estuarine reaches up to the tidal river. (2) Reach-based and common slope linear models for channel perimeter on wetland area were predictive in the estuary (R2 > 0.81) and differed between island (n = 144) and mainland (n = 164) wetlands. No planform linear models were predictive for the entire study area or the tidal river. No models for channel outlets were acceptable. (3) Channel perimeter, surface area, and outlets were all highly variable within and among estuarine–riverine reaches. Analyses of sediment, cross-section, and planform support the use of a “natural network paradigm” for wetland-restoration design on tidal river floodplains.

channel morphology↗

DOE OSTI2026

Modulation induced by very-large-scale motions on the inclination angle of wall-attached eddies: an atmospheric surface layer experiment

The forward leaning inclination angle, 𝛾, of coherent turbulent structures is a well-known feature of wall-bounded turbulent flows. Although invariant across friction Reynolds numbers within the range 𝑅𝑒 𝜏 =10 3 −10 6 , 𝛾 can vary significantly across turbulent scales within a high-Reynolds-number flow. Very-large-scale motions (VLSMs) are known to induce significant changes in the instantaneous shear profile, which is a conditioning event that could trigger variability in the inclination angle of smaller coherent turbulent structures. Although this aspect has been extensively studied via numerical and laboratory experiments, few studies have explored this feature for a very-high-Reynolds-number atmospheric flow. In this work, the inclination angle of turbulent structures within the atmospheric surface layer at a very high Reynolds number (𝑅𝑒 𝜏 =7.9 ×10 5 ) is investigated by deploying a scanning Doppler light detection and ranging and a super large particle image velocimetry (SLPIV) apparatus. The inclination angle of wall-attached eddies is inferred either from the two-point correlation of streamwise velocity (𝛾 =41.1∘) or with a scale-dependent approach through the spectral linear stochastic estimator (SLSE). The SLSE (and, thus, the scale-dependent inclination angle) is conditionally evaluated based on the high- and low-momentum events induced by VLSMs, both in the streamwise (𝑢′𝑉𝐿𝑆𝑀) and in the vertical (𝑤′ 𝑉𝐿𝑆𝑀) velocity components. As a result, lower inclination angles (𝛾 =30° −50°) are found for 𝑢$^{'}_{VLSM}$ >0 (𝑤$^{'}_{VLSM}$ <0), while higher values (50° −85°) are ascribed to 𝑢$^{'}_{VLSM}$ <0 (𝑤$^{'}_{VLSM}$ >0). This result emphasises the primary role that VLSMs play in shaping the wall-attached eddy geometry, which, in turn, is crucial to determine the Reynolds stress balance within the wall-attached eddy range.

42 ENGINEERING↗

DOE OSTI2021

HD 63021: Chromospheric Activity and Mass Transfer in a Close Binary

Prompted by X-ray detections from multiple surveys, we investigated the A-type star HD 63021 and found that it is a double-lined spectroscopic binary with highly variable emission associated with the primary star. Analysis of our multiepoch spectroscopic observations, the majority of which were carried out on small-aperture telescopes, indicates a very short orbital period of just 2.9 days and a mass ratio M {sub 2}/M {sub 1} of 0.23. The A1 V star is a slow rotator, with a rotational speed of ∼34 km s{sup −1}. Assuming that its mass is 2.3 M {sub ⊙}, the present-day secondary is an evolved star of ∼0.5 M {sub ⊙} that nearly fills its Roche lobe. This secondary star rotates comparatively rapidly at ∼44 km s{sup −1}, and we see evidence that it is chromospherically active. Analysis of a photometric light curve from TESS reveals two strong periods, one at the orbital period for the system and another at half the orbital period. These findings suggest that HD 63021 is a close binary system undergoing mass transfer from the secondary star onto the primary star—in all ways like an Algol eclipsing binary system, except without the eclipse. We discuss the system’s mass transfer, which is not steady but seems to run in fits and bursts, and infer the system’s basic physical properties from an orbital parameter study, the Roche lobe geometry, and its extant X-ray emission.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2020

Annihilation-to-nothing: a quantum gravitational boundary condition for the Schwarzschild black hole

The interior of a static Schwarzschild metric can be written in terms of two functions, similar to some models of anisotropic cosmology. With a suitable choice of canonical variables, we solve the Wheeler-DeWitt equation (WDW) inside the horizon of a Schwarzschild black hole. By imposing classicality near the horizon, and requiring boundedness of the wave function, we get a rather generic solution of the WDW equation, whose steepest-descent solution, i.e., the ridge of the wave function, coincides nicely with the classical trajectory. However, there is an ambiguity in defining the arrow of time which leads to two possible interpretations—(i) if there is only one arrow of time, one can infer that the steepest-descent of the wave function follows the classical trajectory throughout: coming from the event horizon and going all the way down to the singularity, while (ii) if there are two different arrows of time in two separate regimes, it can be inferred that the steepest-descent of the wave function comes inwards from the event horizon in one region while it moves outwards from the singularity in the other region, and there exists an annihilation process of these two parts of the wave function inside the horizon. Additionally, adopting the second interpretation could shed light on the information loss paradox: as time goes on, probabilities for histories that include black holes and singularities decay to zero and eventually only trivial geometries dominate.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Microstructure-sensitive mechanical behavior of an additively manufactured psuedoelastic shape memory alloy

The additive manufacturing of shape memory alloys into complex geometries enables fabrication of advanced functional systems across a variety of fields and domains. This work presents results focused on the mechanical behavior of additively manufactured shape memory pseudoelastic NiTi. The deformation induced solid state phase transformation from austenite to martensite allows this system to accommodate large recoverable strains. This deformation behavior is fundamentally driven by crystal-scale transformation physics. Laser powder bed fusion processing reveals that the resulting microstructure, both grain morphology and crystallographic texture, is strongly dependent on the manufacturing processing history. Exhaustive mechanical testing demonstrates that these microstructural factors strongly impact both tensile and cyclic stress–strain behavior. Cyclic dissipative behavior, however, is similar across all tested microstructures following an initial transient period. Remarkably, analysis of spatial strain fields during tensile loading reveals two distinctly different localization “modes”. The first is initiation of localized deformation bands which continuously propagate through the tensile bar during loading. In the second mode localization is observed but lacks propagation; instead additional localization cites nucleate during subsequent loading. The latter phenomena is suspected to be driven by grain-scale deformation physics as the localized band morphologies coincide with grain morphologies. These phenomena strongly impact the resulting aggregate stress–strain behavior. Hence, manufacturers and designers of psuedoelastic functional components must at the very least consider the potential variability in properties when considering additive manufacturing processing. More ideally the process–structure–property relations can be used to further tailor and optimize final functional performance.

Additive manufacturing↗

DOE OSTI2025

Microreactor Optimization Using Simulation And Economics (mouse)

Microreactor Optimization Using Simulation and Economics (MOUSE) is a tool that integrates both nuclear microreactor design and reactor economics to provide comprehensive evaluations and optimizations. This tool enables stakeholders to explore the interplay between technical and economic variables, guiding them towards effective and competitive microreactor solutions. For the reactor core simulations, MOUSE leverages the OpenMC Monte Carlo Particle Transport Code to perform detailed core simulations for various microreactor designs. The included OpenMC models are 2D core designs of a Liquid Metal Thermal Microreactor (LMTR), a Gas-Cooled TRISO-Fueled Microreactor (GCMR), and a Heat Pipe Microreactor. Beyond core design, MOUSE includes simplified calculations for: - Calculating the masses of heat exchangers within the system. - Mechanical power of pumps. - Estimating the area occupied by various buildings within the nuclear plant. For the economic analysis, MOUSE provides detailed bottom-up cost estimates, encompassing a wide range of costs including preconstruction costs, direct costs, indirect costs, training costs, financial costs, operation & maintenance (O&M) costs, and fuel costs. These cost estimations are developed using data from the MARVEL project and additional literature sources, enabling the calculation of total capital costs and levelized cost of energy for both first-of-a-kind and nth-of-a-kind microreactors. MOUSE also enables analysis of the cost drivers and competitiveness in the electricity market. MOUSE allows users to modify a wide array of technical and economic parameters to evaluate different scenarios and their impacts. Examples of these parameters include: Fuels, coolants, or reflector materials Enrichment levels Control drum materials and geometry Fuel pin geometry and materials Moderator pin geometry and materials Reactor core and reflector dimensions Packing factor for the TRISO particles Nuclear reactor power and reactor burnup Number of sensors Shielding thickness Reactor vessel and guard vessel dimensions Operational staff requirements Number of emergency shutdowns Levelization period Interest rate Construction duration Since MOUSE is powered by the WATTS toolkit, it supports optimization studies, parametric analyses, and uncertainty calculations/propagation. The optimization techniques enable users to identify optimal design and economic configurations. The parametric analysis tools allow users to explore the sensitivity of various parameters, while uncertainty propagation helps quantify the impact of uncertainties on overall performance and cost. User Interface and Workflow: Currently, MOUSE is a command-line-based tool. Users can input various reactor design or economic parameters, modify the designs, run simulations, and visualize results through comprehensive data visualization and reporting capabilities. The typical workflow involves setting up the reactor model, defining economic parameters, running simulations, and analyzing the results to make informed decisions. By combining advanced design calculations with detailed economic modeling, MOUSE provides a robust framework for optimizing nuclear microreactor technologies, enhancing their competitiveness, and guiding stakeholders towards innovative and cost-effective solutions.

Hanna, Botros [Idaho National Laboratory (INL), Id↗

DOE OSTI2026

Cluster bootstrap for cosmological correlators

We show that cosmological wavefunction coefficients associated with n-site chain and loop graphs for a cubic scalar theory in de Sitter spacetime have symbol alphabets given by subsets of A 2n−2 and B 2n−1 cluster variables, respectively, and satisfy the associated cluster adjacency properties. The key step in proving this is identifying a precise connection between graph “tubings” that appear in the kinematic flow equation and polygon “triangulations” that encode the combinatorics of cluster compatibility. Our results imply that cosmological wavefunction coefficients in a general power-law FRW cosmology satisfy cluster adjacency to all orders in the ϵ expansion around the de Sitter limit. We use this information as bootstrap input to show that de Sitter symbols for n ≤ 4 are uniquely determined by simple physical constraints.

differential and algebraic geometry↗

DOE OSTI2023

Fuel Injection Dynamics and Composition Effects on RDE Performance

Rotating detonation engines (RDEs) provide a promising route to substantially increasing cycle efficiency in stationary gas turbines. Much of this increase relies on the ability to achieve consistent pressure gain within the combustor. In particular, the design of injectors that feed fuel and air into the detonation channel plays a crucial role. Such injectors have to ensure proper mixing of fuel and oxidizer, while minimizing backflow of detonation products into the feed plenums, and reduce susceptibility to the complex wave structures that exist within the combustor. From a practical perspective, such RDEs also need to operate with variable fuel composition. When fuel mixtures with components that possess vastly different oxidation pathways and time-scales are used, there could be additional losses through deflagrative burning instead of detonation-driven heat release. Such sensitivity to the complete flow path is akin to the physics of thermoacoustic instabilities in conventional gas turbines. In this sense, RDEs pose a unique research challenge: the performance of the device relies on the small-scale heat release process, which is highly dependent on the flow interactions within the full-scale system. As a result, canonical flow configurations, instrumented with detailed diagnostics or modeled using high-fidelity tools, but only focus on the small-scale processes will not contain the key system-level interactions. At the same time, macroscopic measurements and models that only capture system-level performance will not provide insight into the key sources of pressure losses. These couplings and sensitivities provide a formidable challenge to both experimental and simulation studies of the effects Thus, a joint experimental/computational program designed specifically to address these challenges was undertaken in this program. The focus of this program was on two key topics: a) the interaction between injector flow and the overall wave dynamics within the combustor, and b) the deflagration/detonation structure in multi-component fuels that are of practical interest. Both topics involve interaction of small-scale heat release processes with the geometry-dependent wave structure. Studies focused on the study of full-scale RDE systems, based on a 6-inch conventional annular geometry. Experimentally RDEs were studied using a combination of diagnostics. A combination of optical diagnostics and aero-thermo-acoustic analysis based on a combination of spectral and mode decomposition analysis was used to identify the dynamics of the detonation wave and other secondary waves that exist in the system. These studies have helped the identification and investigation of injector and detonation dynamics arising from coupling, and how they affect RDE mixing, detonation structure, operability and performance. Performance of RDEs was investigated through thrust stand measurements, which was used to evaluate the effective pressure gain generated by the system through the concept of equivalent available pressure. Optical diagnostics were developed and implemented to investigate the distribution of heat release, across the detonation wave. Novel optical diagnostics of NIR imaging was also developed and applied to investigate the high temperature / high pressure distribution across the detonation wave. In order to complement the experiments, the computational tools were geared to simulate the full experimental setup. GPU-based acceleration of the models and computations were developed to enable rapid simulation of the full system. In addition, the use of adaptive mesh refinement, and unstructured grid formulation, enabled the investigation of realistic geometries studied in the laboratory. The simulations produced a wealth of detail on the structure of the detonation wave under different operating conditions. Emphasis was placed on quantifying mixture pre-burning and the impact on wave propagation and structure.

03 NATURAL GAS↗

DOE OSTI2022

Optimization of the moderators in the STS preliminary design

This report details the results for an optimization of the dimensions of the moderators in the preliminary design of the Spallation Neutron Source Second Target Station (STS). This study uses the optimization algorithms of Dakota and an unstructured mesh model for the moderators in MCNP. More details on the unstructured mesh model and the automated mesh generation can be found in [3]. Parallel to this effort, the same moderator geometries have been optimized using a constructive solid geometry (CSG) MCNP model. More details on this model and its results can be found in [4]. Three optimal designs are selected for each moderator: one that is optimized for maximum peak brightness, one for maximum time-integrated brightness, and one for a combination of peak and time-integrated brightness. The backbone of the optimization work flow is provided by Dakota. For each set of design parameters requested by Dakota, a new solid geometry is automatically built in Creo and SpaceClaim, and subsequently exported to Attila4MC to generate an unstructured mesh geometry for MCNP. After the MCNP calculation is finished, the objective function (e.g., brightness metric) is returned to Dakota. After the new design has been evaluated, a result-file is written, and Dakota proposes the next set of design parameters to be evaluated. The loop continues until a specified convergence criterion has been met. The design parameters of the cylindrical (upper) moderator include the hydrogen radius, the premoderator thickness (top, bottom, radial), the beryllium radius and the horizontal position of the moderator. The crucial design choice is the hydrogen radius. A radius of 62 mm is shown to provide the maximum time-integrated brightness. The maximum peak brightness occurs with a radius of 40 mm. A combined (middle) design, which balances peak and time-integrated brightnesses, is obtained with a hydrogen radius of 50 mm. The premoderator thicknesses and the beryllium radius are slightly larger in the design optimized for time-integrated brightness than in the design optimized for peak brightness. The sensitivity to these two parameters is relatively small close to the optimal configurations. The hydrogen vessel and vacuum vessel wall thicknesses are dependent on the radius of the liquid hydrogen due to structural integrity requirements. The increased wall thicknesses for larger vessels significantly penalize the time-integrated brightness, with the maximum obtainable value reduced by more than 10% relative to earlier studies which used fixed vessel wall thicknesses. The impact of the variable wall thicknesses is much less for the peak brightness and combined brightness designs. The design parameters of the tube (lower) moderator selected for the optimization are the tube length, the annular premoderator thickness, the beryllium radius and the horizontal position of the moderator. The tube length is the crucial parameter and is chosen large (210 mm) and small (125 mm) in the designs optimized for time-integrated and peak brightness respectively. A combined optimal design has a tube length of 170 mm. The premoderator thickness and the beryllium radius are chosen larger in the design optimized for time-integrated brightness.

42 ENGINEERING↗

DOE OSTI2024

The effect of plasma shaping on high density H-mode SOL profiles and fluctuations in TCV

Abstract The impact of plasma shaping on the properties of high density H-mode scrape-off layer (SOL) profiles and transport at the outer midplane has been investigated on Tokamakà configuration variable. The experimental dataset has been acquired by evolving the upper triangularity while keeping the other parameters constant. The scan comprises δ u p values between 0.0 and 0.6, excluding negative triangularity scenarios. Within this study, a transition from type-I edge localised modes to the quasi-continuous exhaust regime takes place from low to high δ u p . The modification of the upstream SOL profiles has been assessed, in terms of separatrix quantities, within the α t turbulence control parameter theoretical framework (Eich et al 2020 Nucl. Fusion 60 056016). The target parallel heat load and the upstream near-SOL density profiles have been shown to broaden significantly for increasing α t . Correspondingly, in the far SOL a density shoulder formation is observed when moving from low to high δ u p . These behaviours have been correlated with an enhancement of the SOL fluctuation level, as registered by wall-mounted Langmuir probes as well as the thermal helium beam diagnostic. Specifically, both the background and the filamentary-induced fluctuating parts of the first wall ion saturation current signal are larger at higher δ u p , with filaments being ejected more frequently into the SOL. Comparison of two pulses at the extremes of the δ u p scan range, but with otherwise same input parameters, shows that the midplane neutral pressure does not change much during the H-mode phase of the discharge. This indicates that indirect effects of the change in geometry, linked to first wall recycling sources, should not play a significant role. The total core radiation increases at high δ u p , on account of a stronger plasma–wall interaction and resulting larger carbon impurity intake from the first wall. This is likely associated to the enhanced first wall fluctuations, as well as a smaller outer gap and the close-to-double-null magnetic topology at high shaping.

Stagni, A. (ORCID:0000000180841544)↗