Investigation of Three Highly Loaded Subsonic-inlet-stage Axial-flow Compressors Employing Varying Radial Gradients of Energy Addition
Explore the source record for details and available documents.
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.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
The influence of current density and pressure gradient profiles in the pedestal on the access to the regimes free from edge localized modes (ELMs) like quiescent H-mode in ITER is investigated. Using the simulator of MHD modes localized near plasma boundary based on the KINX code, calculations of the ELM stability were performed for the ITER plasma in scenarios 2 and 4 under variations of density and temperature profiles with the self-consistent bootstrap current in the pedestal. Low pressure gradient values at the separatrix, the same position of the density and temperature pedestals and high poloidal beta values facilitate reaching high current density in the pedestal and a potential transition into the regime with saturated large scale kink modes. New version of the localized MHD mode simulator allows one to compute the growth rates of ideal peeling-ballooning modes with different toroidal mode numbers and to determine the stability region taking into account diamagnetic stabilization. The edge stability diagrams computations and sensitivity studies of the stability limits to the value of diamagnetic frequency show that diamagnetic stabilization of the modes with high toroidal mode numbers can help to access the quiescent H-mode even with high plasma density but only with low pressure gradient values at the separatrix. The limiting pressure at the top of the pedestal increases for higher plasma density. With flat density profile the access to the quiescent H-mode is closed even with diamagnetic stabilization taken into account, while toroidal mode numbers of the most unstable peeling-ballooning mode decrease from n = 10−40 to n = 3−20.
The space-time conservation element solution element (CESE) method is modified to address the robustness issues of high-aspect-ratio, viscous, near-wall meshes. In this new approach, the dependent variable gradients are evaluated using element edges and the corresponding neighboring solution elements while keeping the original flux integration procedure intact. As such, the excellent flux conservation property is retained and the new edge-based gradients evaluation significantly improves the robustness for high-aspect ratio meshes frequently encountered in three-dimensional, Navier-Stokes calculations. The order of accuracy of the proposed method is demonstrated for oblique acoustic wave propagation, shock-wave interaction, and hypersonic flows over a blunt body. The confirmed second-order convergence along with the enhanced robustness in handling hypersonic blunt body flow calculations makes the proposed approach a very competitive CFD framework for 3D Navier-Stokes simulations.
Additive manufacturing (AM) plays a significant role in the fabrication of regeneratively cooled rocket engine combustion chambers, where internal channels are used to transport high pressure propellants to keep the engine from overheating. Cooling channels are designed with thin walls between the cooling and combustion sides, experiencing high thermal and mechanical stresses from large temperature gradients and internal pressures. Hence, selecting materials suited for this environment and maintaining a high-quality processing for this critical application is of utmost importance. GRCop-alloys (Cu-Cr-Nb) were developed as an ideal material for combustion chambers due to its high conductivity and high strength at temperature, being employed as a structural member for the hotwall with excellent heat transfer capabilities. GRCop-42 (Cu, 4 at% Cr, 2 at% Nb) has significantly matured using the Laser Powder Bed Fusion (L-PBF) process through microstructure, mechanical and thermophysical properties and evaluations using hot-fire testing. Additionally, Laser Powder Directed Energy Deposition (LP-DED) has demonstrated successful deposition and geometry samples using GRCop-42. While the maturation in properties has focused on bulk material properties from standard test specimens, these cannot be extrapolated to be used in thin walls, as different features (i.e., microstructure, porosity, and surface texture) have a greater effect on thin walls, especially on AM structures built using L-PBF and LP-DED. This presentation provides an overview of the effects introduced by reducing the thickness on the fabrication, microstructure, surface metrology and subsequent mechanical properties of GRCop-42 using L-PBF and LP-DED as well as future work planned for this research.
In this paper, we propose a frozen Gaussian approximation (FGA)-based multi-level particle swarm optimization (MLPSO) method for seismic inversion of high-frequency wave data. The method addresses two challenges in it: First, the optimization problem is highly non-convex, which makes hard for gradient-based methods to reach global minima. This is tackled by MLPSO which can escape from undesired local minima. Second, the character of high-frequency of seismic waves requires a large number of grid points in direct computational methods, and thus renders an extremely high computational demand on the simulation of each sample in MLPSO. We overcome this difficulty by three steps: First, we use FGA to compute high-frequency wave propagation based on asymptotic analysis on phase plane; Then we design a constrained full waveform inversion problem to prevent the optimization search getting into regions of velocity where FGA is not accurate; Last, we solve the constrained optimization problem by MLPSO that employs FGA solvers with different fidelity. The performance of the proposed method is demonstrated by a two-dimensional full-waveform inversion example of the smoothed Marmousi model.
The quest for high-efficiency heat-to-electricity conversion has been one of the major driving forces toward renewable energy production for the future. Efficient thermoelectric devices require high voltage generation from a temperature gradient and a large electrical conductivity while maintaining a low thermal conductivity. For a given thermal conductivity and temperature, the thermoelectric power factor is determined by the electronic structure of the material. Low dimensionality (1D and 2D) opens new routes to a high power factor due to the unique density of states (DOS) of confined electrons and holes. The 2D transition metal dichalcogenide (TMDC) semiconductors represent a new class of thermoelectric materials not only due to such confinement effects but especially due to their large effective masses and valley degeneracies. Here, we report a power factor of MoS2 as large as 8.5 mW m –1 K –2 at room temperature, which is among the highest measured in traditional, gapped thermoelectric materials. To obtain these high power factors, we perform thermoelectric measurements on few-layer MoS 2 in the metallic regime, which allows us to access the 2D DOS near the conduction band edge and exploit the effect of 2D confinement on electron scattering rates, resulting in a large Seebeck coefficient. We conclude the demonstrated high, electronically modulated power factor in 2D TMDCs holds promise for efficient thermoelectric energy conversion.
One of the risks of prolonged manned space flight is the exposure of astronauts to radiation from galactic cosmic rays, which contain heavy ions such as (56)Fe. To study the effects of such exposures, experiments were conducted at the Brookhaven National Laboratory by exposing Wistar rats to high-mass, high-Z, high-energy (HZE) particles using the Alternating Gradient Synchrotron (AGS). The biological effectiveness of (56)Fe ions (1000 MeV/nucleon) relative to low-LET gamma rays and high-LET alpha particles for the induction of chromosome damage and micronuclei was determined. The mitotic index and the frequency of chromosome aberrations were evaluated in bone marrow cells, and the frequency of micronuclei was measured in cells isolated from the trachea and the deep lung. A marked delay in the entry of cells into mitosis was induced in the bone marrow cells that decreased as a function of time after the exposure. The frequencies of chromatid aberrations and micronuclei increased as linear functions of dose. The frequency of chromosome aberrations induced by HZE particles was about 3.2 times higher than that observed after exposure to (60)Co gamma rays. The frequency of micronuclei in rat lung fibroblasts, lung epithelial cells, and tracheal epithelial cells increased linearly, with slopes of 7 x 10(-4), 12 x 10(-4), and 11 x 10(-4) micronuclei/binucleated cell cGy(-1), respectively. When genetic damage induced by radiation from (56)Fe ions was compared to that from exposure to (60)Co gamma rays, (56)Fe-ion radiation was between 0.9 and 3.3 times more effective than (60)Co gamma rays. However, the HZE-particle exposures were only 10-20% as effective as radon in producing micronuclei in either deep lung or tracheal epithelial cells. Using microdosimetric techniques, we estimated that 32 cells were hit by delta rays for each cell that was traversed by the primary HZE (56)Fe particle. These calculations and the observed low relative effectiveness of the exposure to HZE particles suggest that at least part of the cytogenetic damage measured was caused by the delta rays. Much of the energy deposited by the primary HZE particles may result in cell killing and may therefore be "wasted" as far as production of detectable micronuclei is concerned. The role of wasted energy in studies of cancer induction may be important in risk estimates for exposure to HZE particles.
In the interest of improving the predictability of high-lift systems at maximum lift conditions, a series of fundamental experiments were conducted to study the effects of adverse pressure gradient on a wake flow. Mean and fluctuating velocities were measured with a two-component laser-Doppler velocimeter. Data were obtained for several cases of adverse pressure gradient, producing flows ranging from no reversed flow to massively reversed flow. While the turbulent Reynolds stresses increase with increasing size of the reversed flow region, the gradient of Reynolds stress does not. Computations using various turbulence models were unable to reproduce the reversed flow.
Accurate reconstruction of the neutron flux distribution within a reactor core is essential for safe and efficient reactor operation. Traditional power shape synthesis in Light Water Reactors relies on hundreds of in-core detectors. However, this approach becomes impractical for Advanced Reactors and Microreactors due to limited space and harsh environments. To address this challenge, we propose a data-driven methodology that combines high-fidelity modeling with real-time ex-core sensor measurements, enabling the reconstruction of core power distribution while minimizing the reliance on intrusive in-core instrumentation. This project began in FY24 and achieved two initial milestones: (1) the definition of a three-year development plan for a Digital Twin framework and (2) the development of high-fidelity neutronics models of the Purdue University Reactor One (PUR-1) using both MCNP6 and OpenMC. The PUR-1 reactor, a zero-power facility, was selected due to its suitability for neutronics-focused modeling and the availability of experimental data for validation. Both models were benchmarked using neutron flux measurements obtained from irradiated gold foils, which were strategically placed within the core during a dedicated campaign in July 2024. This report marks the continuation and completion of those foundational tasks. The OpenMC model has been refined (improved geometric accuracy, expanded cross-section libraries, and refined sampling) and validated using additional experimental data. An updated sensor design—based on quadrupole configuration—was designed to measure both ex-core flux and its spatial gradient. These measurements will serve as inputs to a neural network-based reconstruction algorithm. Finally, the methodology was demonstrated on a two-dimensional test case representative of the heterogeneous material composition of the PUR-1 reactor core. A neural network implementation of the Kirchhoff-Helmholtz integral equation was employed to solve the boundary value problem using peripheral sensor measurements. The preliminary results confirm the strong potential of the proposed approach for accurate and minimally invasive neutron flux reconstruction.
MOSCATO (Molten Salt Chemistry and Transport) is a multiphysics code that provides high-fidelity, coupled simulations of fluid flow, heat transfer, mass transfer, chemistry, electrochemical phenomena, and alloy corrosion for molten salt systems. In FY25, significant developments were made to the code package, enhancing its capabilities for modeling all relevant phenomena within operating moltens salt reactors (MSRs). The developments and activities in FY25 included: 1. Implementation of Level-Set methods to enable modeling of single-bubble behavior in molten salts. In FY25, the Level-Set two-phase flow modeling implementation was improved to simulate single bubble behavior with molten salt media. The large density and viscosity ratios between typical gases and molten salt liquids present challenges for these types of numerical solvers. With enhancements to the pressure projection method, MOSCATO’s Level-Set solver was able to be successfully validated to experiments related to helium bubble rise in stagnant molten salt. The simulated bubble rising velocity showed reasonable good agreement with experimental measurements. The bubble shape and dynamics were also visually compared with experimental snapshots, demonstrating a good qualitative match. 2. Generation of mass transfer correlations for multiphase flow systems. To enable calculations of the tritium transport across the interface between gas bubbles and salt, we modeled high- Schmidt-number mass transfer around a sphere across a broad range of Reynolds numbers. The mesh near the sphere surface was highly refined to resolve steep concentration gradients caused by the low diffusion coefficient. Literature-based mass transfer correlations were compared with the numerical results, and modifications were proposed to improve agreement, particularly at higher Schmidt numbers. These mass transfer correlations were subsequently provided to other national laboratories to help enable high quality mass transfer simulations using lower-order solvers under development within the NEAMS program. 3. Preliminary implementation of a bubbly flow solver. To model bubbly flow in molten salt, we implemented a bubbly flow solver for void fractions less than 5%. To do so, an algebraic relative velocity model that assumes small bubbles with rapid momentum equilibration was added to MOSCATO to compute bubble velocities. Preliminary comparisons with experimental data showed reasonable agreement, and further improvements are underway. 4. Generation of mass transfer correlations for MSRE subchannel The Molten-Salt Reactor Experiment (MSRE) was a landmark historical project that demonstrated the feasibility of molten-salt reactor technology. The MSRE campaign also generated a significant body of experimental data and reports that continue to support molten-salt–related research. In this report, the MSRE core subchannel was used as the reference geometry for a mass transfer study performed with MOSCATO. The geometry and computational mesh were adapted from a previous study, providing adequate resolution for the relatively low Reynolds number in this case. Additional mesh refinement was applied to reach higher Schmidt numbers, enabling the derivation of a reliable mass-transfer correlation for the present scenario. 5. Simulations of oxygen ingressions into molten salt. In the previous fiscal year, we initiated a study to simulate oxygen ingression in stagnant salt. As oxygen enters the salt through its surface, it reacts with Ce 3+ to form solid CeO 2 and other reaction products. To more fully capture the complex diffusion-convection-reaction mechanisms, capabilities for modeling natural convection in the salt vessel were added. These were needed as the flow of the ingressed gas induced flow in the salt caused by surface shear and non-isothermal effects. With these updated physics in place, we were able to successfully reproduce the experimental results for the rate of change of CeCl 3 concentrations versus time. 6. Flow corrosion model validation. In FY24, MOSCATO’s corrosion model was validated against static corrosion experiments. In FY25, this work was extended to a flow corrosion experiment, where FLiNaK salt was driven by natural convection, with initial salt impurities to initiate corrosion. Despite uncertainties in parameters such as elemental diffusion coefficients in the alloy and unknown H + concentrations, the simulations achieved good agreement with experimental results, especially in predicting sample mass losses.
Transition metal dichalcogenides such as molybdenum disulfide (MoS2) may see service in the heart of next-generation nanoelectronic devices, where highly localized power dissipation can produce nontrivial temperature gradients over nanometer-scale distances. Here, we demonstrate that MoS2 is a promising target for plasmon energy expansion thermometry (PEET), a high-spatial resolution temperature mapping technique employed in a scanning transmission electron microscope (STEM) equipped with electron energy loss spectroscopy (EELS). We first use a calibrated, commercial MEMS-style TEM sample heater chip to measure the temperature dependence of the MoS2 bulk plasmon. We corroborate the chip's temperature calibration with Raman thermometry and determine the bulk thermal expansion coefficient (TEC) of MoS2 in the temperature range of 300–1100 K. Applying this TEC value to PEET measurements on a suspended MoS2 flake, we map 70–90 K/μm temperature gradients with a submicrometer spatial resolution.
Designing materials with high Li to Na selectivity is the key for direct extraction of Li from unconventional sources, such as seawater, without complex separations. However, material discovery is hindered by an insufficient understanding on how a material's composition and structure relate to Li selectivity during electrochemical co-intercalation. Here, we synthesized a core-shell structured (NaLi) 1-x CoO 2 with a core Li phase (Li 0.94 CoO 2 ), a shell Na phase (Na 0.51 CoO 2 ), and a transition intermediate phase via Na ion-exchange of parent Li 1-x CoO 2 . We reveal that the chemical compositions and spatial relation of these phases play critical roles in governing high Li to Na selectivity. The core Li phase restricts large layer spacing expansion and inhibits Na intercalation. The shell Na phase serves to maintain the electrode stability against further Na ion exchange. With the layer spacing gradient design, we achieved a high Li selectivity of 1.5 x 10 5 with a recovery of similar to 7.6:1 Li to Na from 1:20,000 initial ratio.
In 3D seismic interpretation, curvature is a popular attribute that depicts the geometry of seismic reflectors and has been widely used to detect faults in the subsurface; however, it provides only part of the solutions to subsurface structure analysis. This study extends the curvature algorithm to a new curvature gradient algorithm, and integrates both algorithms for fracture detection using a 3D seismic test data set over Teapot Dome (Wyoming). In fractured reservoirs at Teapot Dome known to be formed by tectonic folding and faulting, curvature helps define the crestal portion of the reservoirs that is associated with strong seismic amplitude and high oil productivity. In contrast, curvature gradient helps better define the regional northwest-trending and the cross-regional northeast-trending lineaments that are associated with weak seismic amplitude and low oil productivity. In concert with previous reports from image logs, cores, and outcrops, the current study based on an integrated seismic curvature and curvature gradient analysis suggests that curvature might help define areas of enhanced potential to form tensile fractures, whereas curvature gradient might help define zones of enhanced potential to develop shear fractures. In certain fractured reservoirs such as at Teapot Dome where faulting and fault-related folding contribute dominantly to the formation and evolution of fractures, curvature and curvature gradient attributes can be potentially applied to differentiate fracture mode, to predict fracture intensity and orientation, to detect fracture volume and connectivity, and to model fracture networks.
Transport aircraft performance is strongly influenced by the effectiveness of high-lift systems. Developing wakes generated by the airfoil elements are subjected to strong pressure gradients and can thicken very rapidly, limiting maximum lift. This paper focuses on the effects of various pressure gradients on developing symmetric wakes and on the ability of a linear eddy viscosity model and a non-linear explicit algebraic stress model to accurately predict their downstream evolution. In order to reduce the uncertainties arising from numerical issues when assessing the performance of turbulence models, three different numerical codes with the same turbulence models are used. Results are compared to available experimental data to assess the accuracy of the computational results.
Abstract. Emissions of anthropogenic aerosol and their precursors are often prescribed in global aerosol models. Most of these emissions are spatially heterogeneous at model grid scales. When remapped from low-resolution data, the spatial heterogeneity in emissions can be lost, leading to large errors in the simulation. It can also cause the conservation problem if non-conservative remapping is used. The default anthropogenic emission treatment in the Energy Exascale Earth System Model (E3SM) is subject to both problems. In this study, we introduce a revised emission treatment for the E3SM Atmosphere Model (EAM) that ensures conservation of mass fluxes and preserves the original emission heterogeneity at the model-resolved grid scale. We assess the error estimates associated with the default emission treatment and the impact of improved heterogeneity and mass conservation in both globally uniform standard-resolution (∼ 165 km) and regionally refined high-resolution (∼ 42 km) simulations. The default treatment incurs significant errors near the surface, particularly over sharp emission gradient zones. Much larger errors are observed in high-resolution simulations. It substantially underestimates the aerosol burden, surface concentration, and aerosol sources over highly polluted regions, while it overestimates these quantities over less-polluted adjacent areas. Large errors can persist at higher elevation for daily mean estimates, which can affect aerosol extinction profiles and aerosol optical depth (AOD). We find that the revised treatment significantly improves the accuracy of the aerosol emissions from surface and elevated sources near sharp spatial gradient regions, with significant improvement in the spatial heterogeneity and variability of simulated surface concentration in high-resolution simulations. In the next-generation E3SM running at convection-permitting scales where the resolved spatial heterogeneity is significantly increased, the revised emission treatment is expected to better represent the aerosol emissions as well as their lifecycle and impacts on climate.
Lithium‐ion batteries with high gravimetric capacity density and improved cycle life performance under fast‐charging conditions are crucial for widespread electric vehicle (EV) adoption. This study investigates how designing graphite anode microstructure, specifically porosity, and particle‐size gradients, improves lithium‐ion (Li + ) transport during fast‐charging conditions. Three‐layered graphite anodes with varying porosity (24%, 36%, 46%) and particle size gradients (3, 5, 10 μm) were compared to a conventional single‐layered electrode in half‐cell configurations. At room temperature and high discharge rate (2C), both gradient structures showed significantly enhanced capacity retention (80% and 67% vs. 50%) compared to the conventional electrode, highlighting the effectiveness of microstructure engineering for fast charging. The study also investigated the temperature's impact on cycle life. After 200 cycles at 2°C and 45°C, all gradient structures demonstrated superior capacity retention (≈80%) compared to the conventional electrode (35%), suggesting the gradients mitigate degradation rate at high temperatures. Electrochemical impedance spectroscopy confirmed superior Li+ diffusion and lower resistivity in gradient electrodes. Simulations explored the influence of gradient profiles on reaction kinetics across the electrode thickness. Overall, this research demonstrates that the fast‐charging capability of graphite electrodes can be greatly enhanced by engineering the electrode microstructure, thereby making EV technology more accessible and appealing.
The infra-red video bolometer (IRVB) is a diagnostic equipped with an infra-red camera that measures the total radiated power in thousands of lines of sight within a large field of view. Recently validated in MAST-U [Fderici et al., Rev. Sci. Instrum. 94, 033502 (2023)], it offers a high spatial resolution map of the radiated power in the divertor region, where large gradients are expected. The IRVB’s sensing element comprises a thin layer of high Z absorbing material, typically platinum, usually coated with carbon to reduce reflections [Peterson et al., Rev. Sci. Instrum. 79, 10E301 (2008)].Here, the possibility of using a relatively inert material such as titanium, is explored that can be produced in layers up to 1 μm compared to 2.5 μm for Pt and then coat it with Pt of the desired thickness (0.3 μm per side here) and carbon. This leads to a higher temperature signal (about 3 times) and better spatial resolution (about 4 times), resulting in higher accuracy in the measured power [Peterson et al., Rev. Sci. Instrum. 79, 10E301 (2008)]. This assembly is also expected to improve foil uniformity, as the Pt layer is obtained via deposition rather than mechanical processes [Mukai et al., Rev. Sci. Instrum. 87, 2014 (2016)].Given its multi-material composition, measuring the thermal properties of the foil assembly is vital. Various methods using a calibrated laser as a heat source have been developed, analyzing the temperature profile shape [Sano et al., Plasma and Fusion Res. 7, 2405039 (2012)] and [Mukai et al., Rev. Sci. Instrum. 89, 10E114 (2018)] or fitting the calculated laser power for different intensities and frequencies [Fderici et al., Rev. Sci. Instrum. 94, 033502 (2023)]. Here, a simpler approach is presented, which relies on analyzing the separate components of the foil heat equation for a single laser exposure in a given area. This can then be iterated over the entire foil to capture local deviations.