Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Stability boundary”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

NSTX-U theory, modeling and analysis results

Here, the mission of the low aspect ratio spherical tokamak NSTX-U is to advance the physics basis and technical solutions required for optimizing the configuration of next-step steady-state tokamak fusion devices. NSTX-U will ultimately operate at up to 2 MA of plasma current and 1 T toroidal field on axis for 5 s, and has available up to 15 MW of neutral beam injection power at different tangency radii and 6 MW of high harmonic fast wave heating. With these capabilities NSTX-U will develop the physics understanding and control tools to ramp-up and sustain high performance fully non-inductive plasmas with large bootstrap fraction and enhanced confinement enabled via the low aspect ratio, high beta configuration. With its unique capabilities, NSTX-U research also supports ITER and other critical fusion development needs. Super-Alfvénic ions in beam-heated NSTX-U plasmas access energetic particle (EP) parameter space that is relevant for both α-heated conventional and low aspect ratio burning plasmas. NSTX-U can also generate very large target heat fluxes to test conventional and innovative plasma exhaust and plasma facing component solutions. This paper summarizes recent analysis, theory and modelling progress to advance the tokamak physics basis in the areas of macrostability and 3D fields, EP stability and fast ion transport, thermal transport and pedestal structure, boundary and plasma material interaction, RF heating, scenario optimization and real-time control.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Energy transport analysis of NSTX plasmas with the TGLF turbulent and NEO neoclassical transport models

This work presents a study of plasma transport at low aspect ratio on the National Spherical Torus Experiment tokamak, where the turbulent and neoclassical energy fluxes calculated by the quasilinear Trapped Gyro Landau Fluid (TGLF) model and the multi species drift-kinetic Neoclassical solver (NEO) are validated against experimental data. The turbulent energy transport of two plasma discharges, one in the L-mode confinement regime and another in the H-mode regime, is dominated by electrostatic drift-wave instabilities, while the ion heat transport has a significant neoclassical contribution. The data analysis workflow is described in detail to understand how the variations of mapping and fitting of experimental data affect the power balance solution and subsequent flux-matching plasma profile predictions with the TGYRO solver. On average, the predicted plasma profiles are consistent with experimental data. However, the solutions are sensitive to various input parameters, including boundary conditions, and the electron-ion coupling. Linear gyrokinetic stability analysis demonstrates close agreement of the real frequencies of unstable modes between TGLF and CGYRO gyrokinetic simulations, but higher growth rates are predicted by TGLF, especially for the H-mode case. Estimates of the low-k modes' contributions to the total flux are consistent with linear stability analysis and the E × B suppression of turbulence in TGLF simulations with the SAT1 saturation model, while the SAT2 saturation model over-predicts the low-k modes' contribution in the H-mode case. Moreover, the results with SAT1 model are consistent with power balance analysis, which indicates only neoclassical ion energy fluxes inside ρ < 0.4 in the L-mode case and $\rho \unicode{x2A7D} 0.7$ in the H-mode case. The presence of multi-scale turbulence and ion-scale driven zonal flow mixing effects are also observed in TGLF scans of the electron turbulent heat flux over a range of temperature gradients and the electron-ion temperature ratio, which could explain the strong model sensitivity to variations of input parameters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Correlation of the L-mode density limit with edge collisionality

The "density limit'' is one of the fundamental bounds on tokamak operating space, and is commonly estimated via the empirical Greenwald scaling. This limit has garnered renewed interest in recent years as it has become clear that ITER and many tokamak pilot plant concepts must operate near or above the widely-used Greenwald limit to achieve their objectives. Evidence has also grown that the Greenwald scaling - in its remarkable simplicity - may not capture the full complexity of the disruptive density limit. In this study, we assemble a multi-machine database to quantify the effectiveness of the Greenwald limit as a predictor of the L-mode density limit and identify alternative stability metrics. We find that a two-parameter dimensionless boundary in the plasma edge, $\nu_{*\rm, edge}^{\rm limit} = 3.0 \beta_{T,{\rm edge}}^{-0.4}$, achieves significantly higher accuracy (true negative rate of 97.7\% at a true positive rate of 95\%) than the Greenwald limit (true negative rate 86.1\% at a true positive rate of 95\%) across a multi-machine dataset including metal- and carbon-wall tokamaks (AUG, C-Mod, DIII-D, and TCV). The collisionality boundary presented here can be applied for density limit avoidance in current devices and in ITER, where it can be measured and responded to in real time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Zonal magnetic fields regulate nonlinear edge-localized-mode dynamics via self-consistent force balance

Edge-localized modes (ELMs) eject intense bursts of heat and particles that threaten plasma-facing components in fusion reactors. Nonlinear full-torus BOUT++ simulations show that turbulence-driven zonal magnetic fields (ZMFs) play an essential role in nonlinear ELM evolution by maintaining self-consistent force balance. Zonal flows mitigate the initial crash through shear but do not prevent continued radial transport. When ZMFs are self-consistently included, turbulence-driven zonal currents modify the parallel current distribution and magnetic tension and are associated with a reduction of the axisymmetric (𝑛 = 0) perturbed radial force imbalance. This coincides with a transition from convective, bursty propagation to more localized, diffusive transport. Similar behavior is observed across the regimes considered, including both resistive-ballooning and peeling-ballooning cases. Finally, associated signatures, including radial electric field shear and parallel current redistribution, provide experimentally accessible diagnostics for present devices and ITER-relevant conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

In situ ion irradiation of amorphous TiO 2 nanotubes

Understanding of structural and morphological evolution in nanomaterials is critical in tailoring their functionality for applications such as energy conversion and storage. Here, we examine irradiation effects on the morphology and structure of amorphous TiO 2 nanotubes in comparison with their crystalline counterpart, anatase TiO 2 nanotubes, using high-resolution transmission electron microscopy (TEM), in situ ion irradiation TEM, and molecular dynamics (MD) simulations. Anatase TiO 2 nanotubes exhibit morphological and structural stability under irradiation due to their high concentration of grain boundaries and surfaces as defect sinks. On the other hand, amorphous TiO 2 nanotubes undergo irradiation-induced crystallization, with some tubes remaining only partially crystallized. The partially crystalline tubes bend due to internal stresses associated with densification during crystallization as suggested by MD calculations. Furthermore, these results present a novel irradiation-based pathway for potentially tuning structure and morphology of energy storage materials.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Assessing Strain Rate Sensitivity of Nanotwinned Al–Zr Alloys through Nanoindentation

Nanotwinned metals have exhibited many enhanced physical and mechanical properties. Twin boundaries have recently been introduced into sputtered Al alloys in spite of their high stacking fault energy. These twinned Al alloys possess unique microstructures composed of vertically aligned Σ3(112) incoherent twin boundaries (ITBs) and have demonstrated remarkable mechanical strengths and thermal stability. However, their strain rate sensitivity has not been fully assessed. A modified nanoindentation method has been employed here to accurately determine the strain rate sensitivity of nanotwinned Al–Zr alloys. The hardness of these alloys reaches 4.2 GPa while simultaneously exhibiting an improved strain rate sensitivity. The nanotwinned Al–Zr alloys have shown grain size-dependent strain rate sensitivity, consistent with previous findings in the literature. This work provides insight into a previously unstudied aspect of nanotwinned Al alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Distinct regional meteorological influences on low-cloud albedo susceptibility over global marine stratocumulus regions

Abstract. Marine stratocumuli cool the Earth effectively due to their high reflectance of incoming solar radiation and persistent occurrence. The susceptibility of cloud albedo to droplet number concentration perturbations depends strongly on large-scale meteorological conditions. Studies focused on the meteorological dependence of cloud adjustments often overlook the covariability among meteorological factors and their geographical and temporal variability. We use 8 years of satellite observations sorted by day and geographical location to show the global distribution of marine low-cloud albedo susceptibility. We find an overall cloud brightening potential for most marine warm clouds, which is more pronounced over subtropical coastal regions. A weak cloud darkening potential in the annual mean is evident over the remote SE Pacific and SE Atlantic. We show that large-scale meteorological fields from the ERA5 reanalysis data, including lower-tropospheric stability, free-tropospheric relative humidity, sea surface temperature, and boundary layer depth, have distinct covariabilities over each of the eastern subtropical ocean basins where marine stratocumuli prevail. This leads to a markedly different annual cycle in albedo susceptibility over each basin. Moreover, we find that basin-specific regional relationships between key meteorological factors and albedo susceptibilities are absent in a global analysis. Our results stress the importance of considering the geographical distinctiveness of temporal meteorological covariability when scaling up the local-to-global response of cloud albedo to aerosol perturbations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Investigation of onshore wind farm wake recovery with in situ aircraft measurements during AWAKEN

The share of wind power for electricity supply is increasing worldwide. This highly variable resource requires the improved prediction of power output for network stability. The interaction between wind farm wakes and the atmospheric boundary layer (ABL) introduces uncertainties in power production that warrant detailed investigation. The flow downwind of wind farms is characterized by a reduction in wind speed and an increase in turbulence, which both vary with atmospheric conditions. During the American WAKE experimeNt (AWAKEN), the Technische Universität Braunschweig conducted measurement flights with a research aircraft upwind and downwind of onshore wind farms in the southern Great Plains in Oklahoma in the USA. This study utilizes data from 20 flights conducted at approximately hub height in September 2023 to investigate the wind field variability downwind of the wind farms and vertical profiles to observe atmospheric stratification. The flights were aligned perpendicular to the main wind direction downwind of the King Plains and Armadillo Flats wind farms. Additionally, lidar data from both upwind and downwind ground-based measurement sites and sonic anemometer data were used for comprehensive analysis. Results indicate that under stable ABL conditions, the wake persists at greater downwind distances with a higher velocity deficit in the wake relative to the undisturbed flow compared to unstable stratification. In homogeneous terrain under stable conditions, wake recovery to 95 % occurs between a distance of 4.5 and 9 km downwind of the wind farm. In the semi-complex terrain characterized by shallow hills, slopes, and valleys, the wake exhibits a higher velocity deficit compared to homogeneous terrain, while in some cases the wake was amplified by the terrain resulting in higher velocity deficit 10 km downwind of the wind farm compared to the measurements closer to the wind farm. The turbulent kinetic energy (TKE) and “TKE difference” was found to be a valuable measure in understanding wakes in a semi-complex terrain, showing a clear wake recovery and formation depending on the stratification of the ABL.

17 WIND ENERGY↗
DOE OSTI2025

Bulk nanocrystalline Al–Mg–Y alloys with amorphous grain boundary complexions display high strength and compressive plasticity

Although nanocrystalline alloys regularly exhibit high strengths, their use in structural applications often face challenges due to sample size limitations, unstable microstructures, and the limited ability to plastically deform. The incorporation of amorphous grain boundary complexions has been proposed to address these issues, by simultaneously stabilizing nanocrystalline grain structures for scale-up processing and improving alloy toughness. In the present study, the mechanical behavior of bulk nanocrystalline Al–Mg–Y is examined with macroscale compression testing, probing a length scale that is relevant to real-world structural applications. Bulk samples were fabricated via a simple powder metallurgy approach, with different hot-pressing temperatures and durations employed for consolidation in order to investigate microstructural and property evolution. All of the specimens contained primary face-centered cubic Al and secondary Al 4 C 3 and Al 3 Y phases, with the Al 3 Y particles exhibiting two populations of small equiaxed and larger elongated particles. Appreciable plasticity was measured along with high ultimate stresses over 800 MPa due to the presence of amorphous grain boundary complexions. Microstructural characterization of fracture surfaces revealed that the area fraction of dimpled regions increased with longer hot-pressing time. Most importantly, the elongated Al 3 Y particles formed regular cellular patterns with increasing hot-pressing time, delaying shear localization and significantly enhancing plasticity. The hierarchy present in the microstructure of the Al–Mg–Y alloy, from amorphous grain boundary complexions to secondary phases, gives rise to excellent bulk mechanical properties, which are attractive for structural applications.

Bulk nanocrystalline alloy↗
DOE OSTI2026

An implicit-explicit time splitting strategy for the far SOL plasma fluid model with DG-FEM discretization

We consider a far scrape-off layer (SOL) plasma fluid model of ions that is governed by a Braginskiitype model: a one-dimensional, nonlinear system of advection-diffusion equations coupled with a diffusion equation for neutral particles. Our motivation for studying this system arises from the coupling between the edge plasma and radio-frequency (RF) heating, where solving a far SOL plasma fluid model provides critical insights into edge plasma dynamics. Numerical simulations of plasma fluid models require advanced computational techniques to achieve both efficiency and accuracy, especially when resolving the boundary layer in magnetically confined plasmas. In this work, we propose an implicit-explicit time operator splitting strategy that allows for an efficient solution algorithm, where the diffusive terms are treated semi-implicitly requiring only a linear solve, while the advection part is handled explicitly using a strong-stability-preserving Runge-Kutta (SSP-RK3) scheme. This leads to a fully decoupled system in which the diffusion and advection sub-problems can be solved separately, simplifying the overall solution procedure and allowing for efficient parallelization, which is particularly relevant for exploring the impact of RF heating on the SOL plasma. The main challenge of the discretization is due to the strong coupling between diffusion and advection, particularly through the boundary conditions. This makes implementation of such a scheme in an accurate and stable manner nontrivial. We discuss in detail how to split the equations and manage boundary conditions to maintain stability and well-posedness for each subsystem. We also describe a spatial discretization approach, based on the discontinuous Galerkin finite element method (DG-FEM) and present numerical results for a one-dimensional system.

Burkovska, Olena [ORNL] (ORCID:0000000163101130)↗
DOE OSTI2020

Precipitation Trade Wind Cumuli within EUREC 4 A Field Campaign Report

The campaign Trades Cumuli Precipitation Initiation and Characterization (TCPIC) within the overarching Elucidating the Role of Clouds-Circulation Coupling in Climate EUREC 4 A campaign took place from 17 January 2020 to 20 February 2020, a total of 35 days. The campaign complemented the more extensive EUREC 4 A campaign focusing on cumulus clouds' role in the Earth's climate. TCPIC concentrated on understanding what triggers the onset and development of precipitation in trade wind cumulus clouds. The campaign requested the deployment of the U.S. Department of Energy Atmospheric Radiation Measurement (ARM) Mobile Facility (AMF) roll, pitch, and heave (RPH) stabilization platform to stabilize two radars on the Maria S. Merian research vessel (RV). The advanced remote-sensing techniques from stabilized radars characterized the entire precipitation life cycle, from the formation of larger cloud drops to precipitation evaporating on its way to the surface. Retrievals of evaporation rates, to track how rain affects the trade wind boundary layer, will be based on the zenith-pointing Doppler observations collected from the stabilized micro rain radar (MRR). Figure 1 shows the instrumentation's deployment on the RV: we installed the AMF RPH stabilization platform and the radars on the uppermost deck (peil deck) of the ship. For installing the devices, we coordinated a test phase of the equipment on the Maria S. Merian RV from 25 to 29 October 2019 in Emden harbor, Germany. In this testing phase, we set up all the connections, tested how to mount the W-band and the micro rain radars firmly on the table, and experimented with guided movements of the table to test the resistance of the whole system. Figure 1 shows the final tests conducted in Emden with the stable platform and the radars mounted.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

An Eulerian Vlasov-Fokker–Planck algorithm for spherical implosion simulations of inertial confinement fusion capsules

Here, we present a numerical algorithm that enables a phase-space adaptive Eulerian Vlasov–Fokker–Planck (VFP) simulation of inertial confinement fusion (ICF) capsule implosions. The approach relies on extending a recent mass, momentum, and energy conserving phase-space moving-mesh adaptivity strategy to spherical geometry. In configuration space, we employ a mesh motion partial differential equation (MMPDE) strategy while, in velocity space, the mesh is expanded/contracted and shifted with the plasma’s evolving temperature and drift velocity. The mesh motion is dealt with by transforming the underlying VFP equations into a computational (logical) coordinate, with the resulting inertial terms carefully discretized to ensure conservation. To deal with the spatial and temporally varying dynamics in a spherically imploding system, we have developed a novel nonlinear stabilization strategy for MMPDE in the configuration space. The strategy relies on a nonlinear optimization procedure that optimizes between mesh quality and the volumetric rate change of the mesh to ensure both accuracy and stability of the solution. Implosions of ICF capsules are driven by several boundary conditions: (1) an elastic moving wall boundary; (2) a time-dependent Maxwellian Dirichlet boundary; and (3) a pressure-driven Lagrangian boundary. Of these, the pressure-driven Lagrangian boundary driver is new to our knowledge. The implementation of our strategy is verified through a set of test problems, including the Guderley and Van-Dyke implosion problems — the first-ever reported using a Vlasov–Fokker–Planck model.

1D2V↗
DOE OSTI2025

Speciation and diffusive dynamics in hydrated grain boundaries of complex oxide Gd2Ti2O7

Abstract Grain boundaries in polycrystalline materials significantly affect their properties, such as ionic transport, corrosion, and chemical durability. The pyrochlore compound (Gd 2 Ti 2 O 7 ) is employed as a model for complex oxides and is known for its diverse applications, including nuclear waste immobilization. Density functional theory-based first-principles molecular dynamics simulations were performed at different temperatures on the hydrated grain boundary system. The results show extensive transformations within the grain boundaries among hydrous water species (OH − , H 2 O, and H 3 O + ). The temperature dependence of self-diffusion coefficients follows Arrhenius behavior, with an activation energy of 35.9 kJ/mol for hydrogen and 46.3 kJ/mol for oxygen. The lifetime of OH − is about three to four times longer than that of H 2 O at temperatures from 800 to 2100 K, suggesting the greater stability of OH − over H 2 O, a unique characteristic of the grain boundaries. The estimated lifetime of the hydrous species decreases as the temperature increases, with an activation energy of 9.9 kJ/mol for OH − and 13.4 kJ/mol for H 2 O. While Gd 3 + is more mobile than Ti 4+ , both the Gd 3 + and Ti 4+ cations are orders of magnitude less mobile than the water species. The results suggest that water species are much more mobile within grain boundaries than in the bulk crystal and have the potential to penetrate deep into polycrystalline materials through grain boundaries, leading to grain boundary degradation and dissolution. The different mobilities of cations in complex oxides can lead to leaching of certain cations and incongruent dissolution during the chemical weathering of Earth and industrial materials.

B. Ghosh, Dipta↗
DOE OSTI2022

Magnetic configuration scans during divertor operation of Wendelstein 7-X

Wendelstein 7-X (W7-X) (Greifswald, Germany) is an advanced stellarator, which uses the modular coil concept to realize a magnetic configuration optimized for fusion-relevant plasma properties. The magnet system of the machine allows a variation of the rotational transform (iota) at the boundary. In the latest W7-X operational phase a dedicated configuration scan has been performed varying the rotational transform between magnetic configurations with iota = 5/4 and iota = 5/5 at the boundary. This paper presents an overview of the experiments and of the main results with respect to confinement and stability. The main observation is an increase of the plasma energy in several intermediate configurations of the scan when the 5/5-islands are close to the plasma boundary but still inside the last-closed-flux-surface. In addition, these configurations showed marked MHD-activity with a crashing behavior related to the 5/5-islands. The corresponding mode amplitude was correlated with the size of the internal 5/5 islands.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Effects of interstitial oxygen on ω transformations and twin formation in bcc NbTaTiHf multi-principal element alloy from first-principles

Transformation- and twinning-induced plasticity (TRIP and TWIP) have been reported to contribute to the low-temperature deformation of some body-centered cubic (bcc) multi-principal element alloys (MPEAs) containing large fractions of group IV transition metals. The influence of interstitial solutes on the mechanisms underlying these forms of plasticity, however, remains unclear. Using first-principles calculations, we study the effects of interstitial O atoms on the relative stability of bcc and ω phases and on unstable and twin boundary stacking fault energy profiles in a representative bcc MPEA with high group-IV elemental fraction: NbTaTiHf. We find that O additions generally promote the relaxation of ω configurations back to their parent bcc structure, therefore inhibiting ω transformation. Calculations of the Rice parameter for bulk bcc and phases, as well as bcc-ω interfaces, further show that ω formation is a potent embrittlement factor, an effect that is enhanced by O additions, suggesting that the formation of bcc-ω interfaces is energetically preferred over the formation of the bulk ω phase. By contrast, the Rice parameter for twin boundaries indicates that these interfaces do not embrittle the material, even with O atoms at twin boundaries, providing a more favorable pathway for plastic deformation compared to ω transformation.

Density functional theory↗
DOE OSTI2025

Flow instabilities in helical-coil steam generators for small modular reactors: A review

Here, this study covers the research and discoveries in two-phase flow-boiling instabilities available in the literature—specifically for a helical-coil steam generator (HCSG), including experimental findings, theoretical research, computational models, and system code analyses—supporting research and development of representative small modular reactors (SMRs). Like other new and advanced reactor systems, water-cooled SMRs require experimental data from both integral and separate thermal-hydraulics test facilities for the verification and validation (V&V) of the computational models and computer codes in order to design and obtain regulatory approval. The complex dynamics of two-phase flow-boiling instabilities includes flow regimes physics phenomena, flow-channel geometries, heat-transfer behavior, and interactions among the solid–liquid-gas within the system boundary, all of which are pivotal for understanding the design and operational challenges of SMRs. This study focuses on identifying the relevant knowledge gaps on boiling instabilities—specifically for a HCSG—and provides insights about future research direction optimizing the transport of thermal energy, mass-flow rates, and boundary conditions that ensure the adequate heat-transfer performance, operational stability, and safety associated with SMR systems.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2021

Implications of the iron oxide phase transition on the interiors of rocky exoplanets

The discovery of an extraordinary number of extrasolar planets, characterized by an unexpected variety of sizes, masses and orbits, challenges our understanding of the formation and evolution of the planets in the Solar System and the perception of the Earth as the prototypical habitable world. Many exoplanets appear to be rocky and yet more massive than the Earth, with expected pressures and temperatures of hundreds of gigapascal and thousands of Kelvin in their deep interiors. At these conditions, the properties of bridgmanite and ferropericlase, expected to dominate their mantles, are largely unconstrained, limiting our knowledge of their interior structure. Here we used nano-second X-ray diffraction and dynamic compression to experimentally investigate the atomic structure and density of iron oxide (FeO), one of the end-members of the (Mg,Fe)O ferropericlase solid solution, up to 700 GPa, a pressure exceeding the core–mantle boundary of a 5 Earth masses planet. Our data document the stability of the high-pressure cesium-chloride B2 structure above 300 GPa, well below the pressure required for magnesium oxide (MgO) to adopt the same phase. These observations, complemented by the calculation of the binary MgO–FeO phase diagram, reveal complex stratification and rheology inside large terrestrial exoplanets.

58 GEOSCIENCES↗
DOE OSTI2021

Analysis of the SBP-SAT Stabilization for Finite Element Methods Part II: Entropy Stability

In the hyperbolic research community, there exists the strong belief that a continuous Galerkin scheme is notoriously unstable and additional stabilization terms have to be added to guarantee stability. In the first part of the series, the application of simultaneous approximation terms for linear problems is investigated where the boundary conditions are imposed weakly. By applying this technique, the authors demonstrate that a pure continuous Galerkin scheme is indeed linearly stable if the boundary conditions are imposed in the correct way. In this work, we extend this investigation to the nonlinear case and focus on entropy conservation. Here, by switching to entropy variables, we provide an estimation of the boundary operators also for nonlinear problems, that guarantee conservation. In numerical simulations, we verify our theoretical analysis.

97 MATHEMATICS AND COMPUTING↗