SEARCH · Engineering Papers
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.
Why Austenitic Stainless Steel for Tritium Service? • High strength and ductility for structural components and pressure boundaries • Austenitic SS is highly resistant to hydrogen isotope embrittlement • Not immune though
Abstract Boron (B) powder injection is a potential alternative to glow discharge boronization as a wall conditioning method for tokamaks. This technique is currently being studied in WEST experiments, during which B powder is injected by an Impurity Powder Dropper developed by PPPL. In order to interpret and analyse experimental trends, and to help develop future experiments, a modelling workflow using a boundary plasma simulation (SOLEDGE-EIRENE) and powder ablation simulation (Dust Injection Simulator) was developed and tested. The effect of adding a B neutral source to simulated deuterium + oxygen (D + O) plasmas was compared to experimental data from the WEST C5 campaign, where B powder was injected in a dedicated experiment. While the impact of B injection on radiated power P rad measurements at the upper divertor was similar, there were significant differences in measurements of P rad , outer strike point electron temperature T e OSP and O-II line intensity at the lower divertor between experiment and simulation. This discrepancy suggests that those parameters were affected by phenomena not present in the simulations, with the most likely candidates being reduced D recycling and a reduced O sourcing from the divertor.
The termolecular, association reaction between OH and NO is a source of nitrous acid (HONO), an important atmospheric trace gas. Rate coefficients for the title reaction as recommended by evaluation panels differ substantially at the temperatures and pressures that prevail in the Earth’s boundary layer where the reaction is in the fall-off regime between low- and high-pressure limiting rate coefficients. Using pulsed laser methods for generation and detection of OH, we have reinvestigated the kinetics of the title reaction at pressures of 22–743 Torr (1 Torr = 1.333 hPa) and temperatures (273, 298, and 333 K) in pure N 2 and in N 2 –H 2 O bath gases. In situ optical absorption measurements were used to rule out any bias due to NO 2 or HONO impurities. Our rate coefficients (k 1 ) in N 2 bath gas are parametrized in terms of low-pressure (k o ) and high-pressure (k ∞ ) rate coefficients and a fall-off parameter (F C ) with k 1,0 N 2 = 7.24 × 10 –31 (T/300 K) –2.17 cm 6 molecule –2 s –1 , k 1,∞ = 3.3 × 10 –12 (T/300 K) –3 cm 3 molecule –1 s –1 , and F C = 0.53. Used with the “Troe” expression for termolecular reactions, these parameters accurately reproduce the current data in the fall-off regime and also capture literature rate coefficients at extrapolated temperatures. The presence of water vapor was found to enhance the rate coefficients of the title reaction significantly. The low-pressure limiting rate coefficient in H 2 O bath gas is a factor 5–6 larger than in N 2 , at room temperature (k 1,o H 2 O = 4.55 × 10 –30 (T/300 K) –4.85 cm 6 molecule –2 s –1 ) indicating that H 2 O is much more efficient in quenching the association complex HONO* through collisional energy transfer. Based on measurements in N 2 –H 2 O mixtures, a parametrization of k 1 including both N 2 and H 2 O as third-body quenchers was derived. Neglecting the effect of H 2 O results, e.g., in an underestimation of k 1 by >10% in the tropical boundary layer.
Thin ferroelectric materials hold great promise for compact nonvolatile memory and nonlinear optical and optoelectronic devices. Herein, an ultrathin in-plane ferroelectric material that exhibits a giant nonlinear optical effect, group-IV monochalcogenide SnSe, is reported. Nanometer-scale ferroelectric domains with ≈90°/270° twin boundaries or ≈180° domain walls are revealed in physical-vapor-deposited SnSe by lateral piezoresponse force microscopy. Atomic structure characterization reveals both parallel and antiparallel stacking of neighboring van der Waals ferroelectric layers, leading to ferroelectric or antiferroelectric ordering. Ferroelectric domains exhibit giant nonlinear optical activity due to coherent enhancement of second-harmonic fields and the as-resulted second-harmonic generation was observed to be 100 times more intense than monolayer WS 2 . This work demonstrates in-plane ferroelectric ordering and giant nonlinear optical activity in SnSe, which paves the way for applications in on-chip nonlinear optical components and nonvolatile memory devices.
Purpose: A preclinical MV-energy photon FLASH radiotherapy system is being designed at Stanford and SLAC National Accelerator Laboratory. Because of the higher energy and dose rate compared to conventional kV-energy photon laboratory-scale irradiators, adequate shielding in a stand-alone cabinet form factor is more challenging to achieve. We present a Monte Carlo simulation of multilayered shielding for a compact self-shielding system without the need for a radiation therapy vault. Methods: A multilayered shielding approach using multiple alternating layers of high-Z and low-Z materials is applied to the self-shielded cabinet to effectively mitigate the secondary radiation produced and to allow the device to be housed in a Controlled Radiation Area outside of a radiation vault. The multilayered shielding approach takes advantage of the properties of high-Z and low-Z radiation shielding materials such as density, cross-section, atomic number of the shielding elements, and products of radiation interactions within each layer. The Monte Carlo radiation transport code, FLUKA, is used to simulate the total effective dose produced by the operation. Results: The multilayered shielding designs proposed and simulated produced effective dose rates significantly lower than monolayer designs with the same total material thickness at the regulatory boundary; this is accomplished through the manipulation of the locations where secondary radiation is produced and reactions due to material properties such as neutron back reflection in hydrogen. Borated polyethylene at 5 wt% significantly increased the shielding performance as compared to regular polyethylene, with the magnitude of the reduction depending upon the order of the shielding material. Conclusions: The multilayered shielding provides a path for shielding preclinical FLASH systems that deliver MV-energy bremsstrahlung photons. Furthermore, this approach promises to be more efficient with respect to the shielding material mass and space claim as compared to shielded vaults typically required for clinical radiation therapy with MV photons.
We conduct pore-scale simulations of two-phase flow using the 2D Rothman–Keller colour gradient lattice Boltzmann method to study the effect of wettability on saturation at breakthrough (sweep) when the injected fluid first passes through the right boundary of the model. We performed a suite of 189 simulations in which a “red” fluid is injected at the left side of a 2D porous model that is initially saturated with a “blue” fluid spanning viscosity ratios M=ν r /ν b ∈[0.001,100] and wetting angles θ w ∈[0°,180°]. As expected, at low-viscosity ratios M=ν r /ν b $\ll$1 we observe viscous fingering in which narrow tendrils of the red fluid span the model, and for high-viscosity ratios M$\gg$1, we observe stable displacement. The viscous finger morphology is affected by the wetting angle with a tendency for more rounded fingers when the injected fluid is wetting. However, rather than the expected result of increased saturation with increasing wettability, we observe a complex saturation landscape at breakthrough as a function of viscosity ratio and wetting angle that contains hills and valleys with specific wetting angles at given viscosity ratios that maximize sweep. This unexpected result that sweep does not necessarily increase with wettability has major implications to enhanced oil recovery and suggests that the dynamics of multiphase flow in porous media has a complex relationship with the geometry of the medium and the hydrodynamical parameters.
Hydrogen is a promising component of a future energy-secure and efficient economy, but its competitiveness depends on reducing production costs. One strategy is to operate alkaline water electrolyzers at higher current densities to increase output. However, this intensifies performance losses due to gas bubble accumulation, which blocks transport pathways and deactivates electrochemically active surfaces. Enhancing bubble evacuation through electrode design is therefore essential. Previous studies have explored various approaches — such as modifying surface morphology, applying sonication or pressure modulation, and introducing surfactants — but these efforts have addressed a limited range of conditions due to the complexity of two-phase flow and electrode geometries. Experiments have also largely been focused on either cell level improvements, which lack the information necessary to isolate each contributing factor, or on modified geometries that are not relevant to practical cell operation. From a modeling perspective, conventional Eulerian multiphase models do not track the complex gas–liquid interfacial dynamics and often neglect surface tension and contact angle effects, reducing their predictive accuracy. To provide insights on the effects of different electrode geometries on the performance of alklaine water electrolyzers this work employs an immersed boundary volume-of-fluid method to simulate bubble behavior in 3D porous electrodes. Multiple base electrode geometries, typically used in practice, with varying porosity are evaluated under a constant surface gas generation rate. Simulation data is analyzed to quantify electrode gas coverage, bubble size dynamics and other relevant metrics. Results show that porosity strongly influences bubble accumulation on electrode surfaces, with higher porosity reducing gas coverage, and its not strictly dependent on the electrode geometry. However, the electrode’s base geometry significantly affects gas accumulation at the separator gap, independent of porosity. A foam electrode geometry resulted in the lowest gas coverage of all electrodes with a median volumetric gas coverage of 11%, but at the cost of a 70% reduction in active area compared with the largest surface area electrode, while gyroid electrodes showed the best trade-off between gas coverage, particularly at the separator surface, and electrochemically active area. In conclusion, the results highlight the need for holistic electrode design strategies.
This paper describes a numerical algorithm for solving the classic problem of a plane strain (KGD) fracture propagating in an impermeable elastic medium with zero toughness. The method, which takes advantage of the self-similar nature of the solution, combines a domain-based scheme to solve the elasticity equations and a finite volume method to solve the nonlinear lubrication equation. This work represents a first step towards developing a model able to account for pore pressure diffusion in the medium and corresponding poroelastic effects, noting that these processes are more efficiently solved using a domain-based rather than a boundary integral method. To enhance the efficiency and accuracy of the numerical scheme, the far-field crack asymptotics is embedded in the discretized elastic relationship between the fluid pressure and the crack opening, while the coupled fluid-solid tip asymptote is enforced in a weak form when solving the nonlinear lubrication equation. The proposed technique yields results that closely match the analytical solution, even with a coarse mesh. This approach offers potential for addressing more complex hydraulic fracturing problems in the future.
Mechanical wave transmission through a material is influenced by the mechanical discontinuity in the material. The propagation of embedded discontinuities can be monitored by analyzing the wave-transmission measurements recorded by a multipoint sensor system placed on the surface of the material. The proposed workflow monitors the propagation of mechanical discontinuity through three stages, namely initial, intermediate, and final stages, by using supervised learning followed by data-driven causal discovery. To the end, the workflow processes the multipoint waveform measurements resulting from a single impulse source, while considering the effects of wave attenuation, dispersion and multiple wave-propagation modes due to the discontinuity and material boundaries. Among various feature reduction techniques ranging from decomposition methods to manifold approximation methods, the features derived based on statistical parameterizations of the measured waveforms lead to reliable monitoring that is robust to changes in precision, resolution, and signal-to-noise ratio of the multipoint sensor measurements. The numbers of zero-crossing, negative-turning, and positive turning in the waveforms are the strongest causal signatures of the propagation of mechanical discontinuity. Higher order moments of the waveforms, such as variance, skewness and kurtosis, are also strong causal signatures of the propagation. Finally, the newly discovered causal signatures confirm that the statistical correlations and conventional feature rankings are not always statistically significant indicators of causality.
In many ferroelectrics, large electromechanical strains are observed near regions of composition- or temperature- driven phase coexistence. Phenomenologically, this is attributed to easy re-orientation of the polarization vector and/or phase transition, although their effects are highly convoluted and difficult to distinguish experimentally. Here, we used synchrotron X-ray scattering and digital image correlation to differentiate between the microscopic mechanisms leading to large electrostrains in an exemplary Pb-free piezoceramic Sn-doped barium calcium zirconate titanate. Large electrostrains of ~0.2% measured at room-temperature are attributed to an unconventional effect, wherein polarization switching is aided by a reversible phase transition near the tetragonal-orthorhombic phase boundary. Additionally, electrostrains of ~0.1% or more could be maintained from room temperature to 140°C due to a succession of different microscopic mechanisms. In situ X-ray diffraction elucidates that while 90° domain reorientation is pertinent below the Curie temperature ( T C ), isotropic distortion of polar clusters is the dominant mechanism above T C .
In this work, we examine the role of image charge effects on the electrostatic potential fluctuations experienced by ionic species in the vicinity of an electrode surface. We combine simulation and theory to quantify these fluctuations and how they vary with distance from the electrode surface. We observe that the potential distribution narrows significantly for species within a few electrolyte screening lengths of the electrode. We attribute this narrowing to the effects of image charge fluctuations originating from the polarization response of the electrode. We show that the physical consequences of these image charge effects can be captured in the context of a simple analytical field theory with anti-symmetric boundary conditions. We contextualize these results by discussing their implications for rates of Marcus-like outer-sphere interfacial electron transfer.
Chemical ordering is a common phenomenon and highly correlated with the properties of solid materials. By means of the redistribution of atoms and chemical bonds, it invokes an effective lattice adjustment and tailors corresponding physical properties. To date, however, directly probing the 3D interfacial interactions of chemical ordering remains a big challenge. In this work, we deciphered the interlaced distribution of nanosized domains with chemical order/disorder in Fe 3 Pt bulk alloy. HAADF-STEM images evidence the existence of such nanodomains. The reverse Monte Carlo method with the X-ray pair distribution function data reveal the 3D distribution of local structures and the tensile effect in the disordered domains at the single-atomic level. The chemical bonding around the domain boundary changes the bonding feature in the disordered side and reduces the local magnetic moment of Fe atoms. This results in a suppressed negative thermal expansion and extended temperature range in Fe 3 Pt bulk alloy with nanodomains. Our study demonstrates a local revelation for the chemical order/disorder nanodomains in bulk alloy. The understanding gained from atomic short-range interactions within the domain boundaries provides useful insights with regard to designing new functional compounds.
In this work, we study ribbons of the dice two-dimensional lattice (that we call “dice ladders”) known to have nontrivial topological properties, such as Chern numbers 2. Our main results are two fold. (1) Analyzing the tight-binding model in the presence of Rashba spin-orbit coupling and an external magnetic field, we observed that dice ladders qualitatively display properties similar to their two-dimensional counterpart all the way to the limit of only two legs in the short direction. This includes flat bands near the Fermi level, edge currents and edge charge localization near zero energy when open boundary conditions are used, two chiral edge modes, and a nonzero Hall conductance. (2) We studied the effect of Hubbard correlation U in the two-leg dice ladder using Lanczos and density matrix renormalization group techniques. We show that by increasing U the flat bands split without the need of introducing external fields. Moreover, robust ferrimagnetic order develops. Overall, our work establishes dice ladders as a promising playground to study the combined effect of topology and correlation effects, one of the frontiers in quantum materials.
We derive the soft effective action in (d+2)-dimensional Abelian gauge theories from the on-shell action obeying Neumann boundary conditions at timelike and null infinity and Dirichlet boundary conditions at spatial infinity. This allows us to identify the on-shell degrees of freedom on the boundary with the soft modes living on the celestial sphere. Following the work of Donnelly and Wall, this suggests that we can interpret soft modes as entanglement edge modes on the celestial sphere and study entanglement properties of soft modes in Abelian gauge theories.
This study explores the added value (AV) of a regional Earth system model (ESM) compared to an atmosphere-only regional climate model (RCM) in simulating West African monsoon (WAM) rainfall. The primary goals are to foster discussions on the suitability of coupled RCMs for WAM projections and deepen our understanding of ocean–atmosphere coupling’s influence on the WAM system. The study employs results from dynamical downscaling of the ERA-Interim reanalysis and Max Plank Institute ESM, low resolution (MPI-ESM-LR), by two RCMs, atmosphere only (REMO) and REMO coupled with Max Planck Institute Ocean Model (MPIOM) (ROM), at ~25-km horizontal resolution. Results show that in regions distant from coupling domain boundaries such as West Africa (WA), constraint conditions from ERA-Interim are more beneficial than coupling effects. REMO, reliant on oceanic sea surface temperatures (SSTs) from observations and influenced by ERA-Interim, is biased under coupling conditions, although coupling offers potential advantages in representing heat and mass fluxes. Contrastingly, as intended, coupling improves SSTs and monsoon fluxes’ relationships under ESM-forced conditions. In this latter case, the coupling features a dipole-like spatial structure of AV, improving precipitation over the Guinea Coast but degrading precipitation over half of the Sahel. Our extensive examination of physical processes and mechanisms underpinning the WAM system supports the plausibility of AV. Additionally, we found that the monsoonal dynamics over the ocean respond to convective activity, with the Sahara–Sahel surface temperature gradient serving as the maintenance mechanism. While further efforts are needed to enhance the coupled RCM, we advocate for its use in the context of WAM rainfall forecasts and projections.
This report documents the findings of a project aiming to improve the initial conditions of km-scale simulations of deep convective storms by assimilating cloud-scale weather radar observations. More accurate numerical analyses will increase the effectiveness of research efforts using LES cloud models as tool to better understand land-atmosphere coupling, boundary layer turbulence, and cloud processes, each used for model parameterization development. Work reported on herein includes: i) assessment of radar data sets for use in data assimilation (‘DA’) experiments, ii) quality control of the radar data set and format conversion to one acceptable to DA schemes utilized by the Weather Research and Forecasting (WRF) model, and iii) examination of the sensitivity of the numerical representation of cloud-scale wind and microphysical features to a variety of tunable DA parameters.
This report presents the microstructural characterization using scanning transmission electron microscopy (STEM) of two creep-tested ESR1150-PT specimens, one short-term and one long-term. The results show that during creep testing, the microstructure of the A709 ESR1150-PT material continues to evolve as expected due to the synergistic effect of strain, temperature and time. Key features include dislocation forest and cell structures, precipitates on grain boundaries, precipitates randomly within grains, and precipitates associated with dislocations. In addition, the STEM study of one ESR1100-SA specimen that has been tested to rupture at a relatively low temperature (550°C) for a relatively long term (16226 h) is also presented. The results show that there is a high density of very fine nano precipitates form on dislocations.
Alumina-forming austenitic (AFA) stainless steels have emerged as a candidate alloy because of their high-temperature strength, formability, cost, and compatibility with primary coolants for lead-cooled fast reactors (LFRs). This class of steels has exceptional high-temperature oxidation performance; however, a high concentration of Ni is required to stabilize the austenite phase and to provide sufficient creep strength. AFA stainless steels are also susceptible to liquid metal embrittlement (LME). Additionally, under neutron irradiation, Ni will enrich at grain boundaries due to radiation-induced segregation (RIS). Nickel RIS can increase the LME under these coupled effects. Oak Ridge National Laboratory (ORNL) and the NSUF program have leveraged its High Flux Isotope Reactor (HFIR) and experience with complex irradiation experiments to design experiment capsules that test the aforementioned coupled effects. These capsules are designed for insertion in the central flux trap, the highest flux region, of HFIR. The experiment capsules will be filled with Pb, designed to passively melt from the gamma heating in HFIR. The specimens were fabricated into miniature tensile specimens from two different alloys, GA05-25Ni and GA05-20Ni, varying Ni concentrations. The experiment capsules are designed to achieve target temperatures of 400 °C and 650 °C with accumulated dosage of 3 dpa. This report documents the specimen alloy characterization, experimental design, and expected performance of the capsules.