SEARCH · Search NASA
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.
We have developed a mathematical model that predicts that hydrodynamic slip can eliminate the diffusioosmotic flow reversal due to electrostatic correlations.
Thermoacoustic instability (TAI) is a critical issue in modern lean-burn gas-turbine combustors, which is induced by a strong coupling between the resonant combustor acoustics and fluctuations in the heat release rate. This instability may lead to high-amplitude pressure waves that generate undesirable noise levels as well as fatigue stresses in mechanical structures of the combustor. The intense pressure fluctuations due to TAI may also cause large flow perturbations and possibly flow reversal that may lead to flame oscillations, flame liftoff, and even flame blow-out. Hence, there is a strong need for exercising control actions in a timely fashion to mitigate the TAI phenomena. Anomaly detection is an essential prerequisite to the design of a good controller and such a detector must be able to reliably predict a forthcoming TAI. To detect and predict the onset of a TAI from an ensemble of pressure time series, this paper investigates three data-driven methods: Fast Fourier transform (FFT), symbolic time series analysis (STSA), and hidden Markov modeling (HMM). The main focus of the paper is to make a comparative evaluation of these three anomaly detection methods for classification of the current regime of operation into stable and unstable categories as well as for real-time identification of precursors to impending instabilities with short-length time series of measured variables (e.g., pressure oscillations). Finally, the results, generated on experimental data from a multi-nozzle combustor apparatus, have been compared to evaluate the performance of FFT, STSA, and HMM methods for TAI analysis.
In fusion neutral gas analysis, such as with the Diagnostic Residual Gas Analyzer (DRGA) for ITER, the primary measurement range of interest comprises the low-amu species (1 to 6), especially deuterium and helium. The challenge in successfully obtaining accurate measurements is two-fold. First, the sensitivity of the method must be sufficient to resolve trace amounts accurately; typically, one percent or less. Second, the gas signal from the fusion processes must be free of bias caused by the latent presence (from system outgassing and/or vacuum backstreaming) of these gases to enable accurate interpretation of the measured signal. This latter criterion can be problematic for the lightest gases since there is a propensity for some fraction of the pumped gas load to undergo a phenomenon known as backstreaming. This behavior is manifested in pumping systems for gas properties related to relative atomic weight (lightest) and size (smallest). Backstreaming results in a significant amount of the pumped gas undertaking a reverse flow and re-entering the measurement region; thus, contaminating the forward, real-time measurement. To fully eliminate this adverse effect, a conductance-limiting device – or orifice – has been installed in the high-vacuum pumping system of the present ITER DRGA prototype. The system was already equipped with a secondary turbomolecular pump (TMP), but with limited effectiveness against backstreaming in the inter-pump volume (IPV). This orifice is placed within the suction inlet coupling of the secondary TMP, which is downstream of the IPV. Its objective is to eliminate the backstreaming phenomenon by increasing the back pressure in the IPV. However, the orifice sizing must take into consideration other factors, such as the diagnostic measurement objectives. For example, in the ITER DRGA, one of the measurement requirements is a dynamic response time of ~1s. Fortunately, an added benefit of the pumping restriction created by the orifice is that the upstream pressure increase is beneficial for the DRGA’s optical gas analysis (OGA) sensors. These sensors are attached to the IPV in the present design. The glow discharges, when used as an OGA light source, will typically have a brighter light emission with increasing plasma cell pressure. In addition to the fusion machine research sector, there are other potential applications of this pumping technique where the monitoring of lighter gas concentrations is essential, such as the photolithography process for the semiconductor fabrication of integrated circuits. This presentation will describe the vacuum system used to demonstrate a process to eliminate backstreaming as well as show test results to verify the accomplishment of this critical objective.
We evaluate two leading-edge-based dynamic stall-onset criteria (namely, the maximum magnitudes of the leading-edge suction parameter and the boundary enstrophy flux) for mixed and trailing-edge stall. These criteria have been shown to successfully predict the onset of leading-edge stall at Reynolds numbers of [Formula: see text], where the leading-edge suction drops abruptly. However, for mixed/trailing-edge stall, leading-edge suction tends to persist even when there is significant trailing-edge reversed flow and stall is underway, necessitating further investigation into the effectiveness of these criteria. Using wall-resolved large-eddy simulations and the unsteady Reynolds-averaged Navier–Stokes method, we simulate one leading-edge stall and three mixed/trailing-edge stall cases at Reynolds numbers of [Formula: see text] and [Formula: see text]. We contrast the progression of flow features such as trailing-edge separation and vortex formation across different stall types and evaluate the stall-onset criteria relative to critical points in the flow. We find that the criteria nearly coincide with the instance of leading-edge suction collapse and are reached in advance of dynamic stall vortex formation and lift stall for all four cases. We conclude that the two criteria effectively signal dynamic stall onset in cases where the dynamic stall vortex plays a prominent role.
This is an M2 report for the NEAMS program. This work focuses on advancing the modeling and simulation capabilities for Molten Salt Reactors (MSRs) within Pronghorn, addressing key thermal-hydraulic challenges that must be addressed for the safe and efficient operation of these advanced reactors. This report summarizes the development and validation of the following new models in Pronghorn: (i) two-equation RANS turbulence models, (ii) near-wall turbulence corrections due to curvature and roughness, and (iii) multidimensional two-phase drift-flux models. For each of these cases, along with the experimental validation, an application study of interest to MSRs is presented. In addition, this report also presents recent improvements in the domain-overlapping coupling between Pronghorn and SAM to address the cases of weakly compressible flow and flow reversal. Finally, an MSR action that greatly simplifies the setup of MSR models using NEAMS tools is described.
We study the structure of reverse shocks formed by the collision of supersonic, magnetized plasma flows driven by an inverse (or exploding) wire array with a planar conducting obstacle. We observe that the structure of these reverse shocks varies dramatically with wire material, despite the similar upstream flow velocities and mass densities. For aluminum wire arrays, the shock is sharp and well-defined, consistent with magneto-hydrodynamic theory. In contrast, we do not observe a well-defined shock using tungsten wires, and instead we see a broad region dominated by density fluctuations on a wide range of spatial scales. We diagnose these two very different interactions using interferometry, Thomson scattering, shadowgraphy, and a newly developed imaging refractometer that is sensitive to small deflections of the probing laser corresponding to small-scale density perturbations. We conclude that the differences in shock structure are most likely due to radiative cooling instabilities, which create small-scale density perturbations elongated along magnetic field lines in the tungsten plasma. These instabilities grow more slowly and are smoothed by thermal conduction in the aluminum plasma.
The rapid adoption of photovoltaic (PV) systems combined with falling prices of energy storage is paving the way for a future in which customers could locally supply their energy needs and export surplus power to the grid. However, reverse power flow from multiple sites in a network can cause overvoltage that decreases system reliability and the utilization level of PV system. In this paper, we present a two-stage approach to solve the energy coordination problem in distribution networks. In the first stage, a household with PV and energy storage is optimized to allow customers to adopt different control strategies such as load following and tariff arbitrage. We show that such customer-led control strategies can cause reverse power flows during peak PV generation periods, thus resulting in overvoltage. In such cases, the second stage of the proposed approach is called upon to solve a fairness-based energy coordination problem that can maintain voltages within limits while maximizing the PV generation while fairly distributing the actions among all the PV systems. As a result, each of the resources in the network has equal opportunity to export power and shares the responsibility of voltage regulation. Simulation studies have been carried out on the IEEE 13-bus test feeder to verify the effectiveness of the proposed approach.
During photoheterotrophic growth on organic substrates, purple nonsulfur photosynthetic bacteria like Rhodospirillum rubrum can acquire electrons by multiple means, including oxidation of organic substrates, oxidation of inorganic electron donors (e.g., H 2 ), and by reverse electron flow from the photosynthetic electron transport chain. These electrons are stored as reduced electron-carrying cofactors (e.g., NAD(P)H and ferredoxin). The overall ratio of oxidized to reduced cofactors (e.g., NAD(P)+:NAD(P)H), or ’redox poise’, is difficult to understand or predict, as are the cellular processes for dissipating these reducing equivalents. Using physics-based models that capture mass action kinetics consistent with the thermodynamics of reactions and pathways, a range of redox conditions for heterophototrophic growth are evaluated, from conditions in which the NADP+/NADPH levels approach thermodynamic equilibrium to conditions in which the NADP+/NADPH ratio is far above the typical physiological values. Modeling predictions together with experimental measurements indicate that the redox poise of the cell results in large-scale changes in the activity of biosynthetic pathways and, thus, changes in cell macromolecule levels (DNA, RNA, proteins, and fatty acids). Furthermore, modeling predictions indicate that during phototrophic growth, reverse electron flow from the quinone pool is a minor contributor to the production of reduced cofactors (e.g., NAD(P)H) compared to other oxidative processes (H 2 and carbon substrate oxidation). Instead, the quinone pool primarily operates to aid ATP production. The high level of ATP, in turn, drives reduction processes even when NADPH levels are relatively low compared to NADP+ by coupling ATP hydrolysis to the reductive processes. The model, in agreement with experimental measurements of macromolecule ratios of cells growing on different carbon substrates, indicates that the dynamics of nucleotide versus lipid and protein production is likely a significant mechanism of balancing oxidation and reduction in the cell.
During photoheterotrophic growth on organic substrates, purple nonsulfur photosynthetic bacteria like Rhodospirillum rubrum can acquire electrons by multiple means, including oxidation of organic substrates, oxidation of inorganic electron donors (e.g., H2), and by reverse electron flow from the photosynthetic electron transport chain. These electrons are stored as reduced electron-carrying cofactors (e.g., NAD(P)H and ferredoxin). The overall ratio of oxidized to reduced cofactors (e.g., NAD(P)+:NAD(P)H), or ’redox poise’, is difficult to understand or predict, as are the cellular processes for dissipating these reducing equivalents. Using physics-based models that capture mass action kinetics consistent with the thermodynamics of reactions and pathways, a range of redox conditions for heterophototrophic growth are evaluated, from conditions in which the NADP+/NADPH levels approach thermodynamic equilibrium to conditions in which the NADP+/NADPH ratio is far above the typical physiological values. Modeling predictions together with experimental measurements indicate that the redox poise of the cell results in large-scale changes in the activity of biosynthetic pathways and, thus, changes in cell macromolecule levels (DNA, RNA, proteins, and fatty acids). Furthermore, modeling predictions indicate that during phototrophic growth, reverse electron flow from the quinone pool is a minor contributor to the production of reduced cofactors (e.g., NAD(P)H) compared to other oxidative processes (H2 and carbon substrate oxidation). Instead, the quinone pool primarily operates to aid ATP production. The high level of ATP, in turn, drives reduction processes even when NADPH levels are relatively low compared to NADP+ by coupling ATP hydrolysis to the reductive processes. The model, in agreement with experimental measurements of macromolecule ratios of cells growing on different carbon substrates, indicates that the dynamics of nucleotide versus lipid and protein production is likely a significant mechanism of balancing oxidation and reduction in the cell.
Solid-state transformer (SST) and hybrid transformer (HT) are promising alternatives to the line-frequency transformer (LFT) in smart grids. The SST features medium-frequency isolation, full controllability for voltage regulation, reactive power compensation, and the capability of battery energy storage system (BESS) integration with multiport configuration. The HT has a partially-rated converter for fractional controllability and can integrate a small BESS. Fast grid-edge voltage fluctuations from increased solar photovoltaic (PV) and electric vehicle (EV) penetration are difficult to manage for mechanical load tap changers. Hence, along with the trend towards more BESS in the grid, the controllability and the storage integration capability of the SST and HT are of strong interest. However, a review of literature shows existing SST and HT research is mostly at converter level, while system-level assessments are scarce. Assessing technical and economic impacts is critical to understanding the benefits and role of the SST and HT to guide future research, which is presented for the first time in this article. Experimental results from medium-voltage (MV) SST and MV HT prototypes are shown to confirm equipment-level feasibility, where the voltage controllability waveforms of a MV HT prototype are reported for the first time. Comparative simulations are performed on a modified IEEE 34-bus system. Here, a grid-model-less decentralized grid-edge voltage control method and a day-ahead BESS dispatch method are proposed for the SST and HT. The simulations show that the SST and HT with integrated storage can host more PV, achieve peak shaving, mitigate voltage fluctuation and reverse power flow, and support energy arbitrage for operational cost reduction, as compared to the LFT. Moreover, comprehensive analyses of net present value (NPV) and internal rate of return (IRR) are performed under different installed PV capacities, HT’s partial converter ratings, and BESS capacities. Sensitivities to future cost reductions of the PV and BESS are studied. Although the NPV and IRR are currently negative, 60% capital cost reduction or 150% revenue increase will make the SST and HT economically viable in the use case studied.
An electron temperature bifurcation is observed in the small angle slot divertor, which has been developed to enhance neutral cooling across the divertor target by coupling a closed slot structure with appropriate target shaping. Experiments in the DIII-D tokamak and associated SOLPS-ITER modeling with full drifts find a strong interplay between drifts and divertor geometry on divertor dissipation. The coupling of divertor geometry and drift flows can strongly affect the path towards divertor detachment onset as the plasma density is raised. With the strike point on the inner slanted surface and ion B × ∇B away from the magnetic X-point, bifurcative transitions were observed with sharp decrease of T e towards detachment onset both experimentally and computationally. This differs from the situation for the open divertor where the T e cliff was only observed for ion B × ∇B towards the X-point. SOLPS-ITER modeling with full drifts demonstrates that the magnitude of the E × B drift flow is comparable with the main plasma flow. The reversal of both the poloidal and radial E × B flows near the strike point leads to rapid density accumulation right near the separatrix, which results in bifurcative step transition of divertor conditions with cold plasma across the entire divertor target plate. Furthermore, these results indicate that the interplay between geometry and drifts should be fully taken into account in future fusion reactor divertor designs.
We studied a variety of absorber thicknesses in Cu(In,Ga)Se 2 devices to examine the ability of the cell to withstand reverse current flow without damage. Absorbers were varied from 2.5 um to 0.4 um, with the expectation that larger reverse current would be allowed to flow at lower voltages as absorbers were made thinner. Our initial experiments showed promise. However, as the efficiency of thin devices was improved, the reverse current allowed at a given voltage decreased. Here we present our negative results on the ability of thin devices to allow higher reverse currents at lower voltages. We also give details on our procedures for reverse-bias testing small-area solar cells. We hope this information will be useful for those who want to study reverse-bias at the cell level.
The Prototype Material Plasma Exposure eXperiment (Proto-MPEX) at Oak Ridge National Laboratory is developing a helicon radio frequency (RF) source for plasma production along with microwave electron heating and RF ion heating. Proto-MPEX has several magnetic mirrors that affect plasma transport. The magnetic field near the helicon source is ~ 0.05 T, and the peak magnetic field in Proto-MPEX is ~ 1.8 T, which results in a large magnetic well for particles sourced at the helicon transporting toward the target region. Keeping the helicon source field and the target field constant, experiments have been conducted by methodically changing the magnetic field, at a local axial location, within the source and the target. Plasma flow measurements using Mach probes have shown a decrease in the flow as a function of the magnetic mirror ratio. Plasma flows were further investigated with ion cyclotron resonant heating (ICH). Results with the addition of ICH measured upstream of the antenna have shown a decrease in Mach number as a function of mirror ratio when compared to helicon discharges. Such flow behavior demonstrated slowing down of the plasma due to the increased reverse plasma flow toward the helicon plasma source. However, in high magnetic fields, when the resonant magnetic field moved away from the ICH antenna, no effects on the plasma flow due to ICH were observed.
The multi-region relaxed magnetohydrodynamics (MRxMHD) has been successful in the construction of equilibria in three-dimensional (3D) configurations. In MRxMHD, the plasma is sliced into sub-volumes separated by ideal interfaces, each undergoing relaxation, allowing the formation of islands and chaos. The resulting equilibrium has a stepped pressure profile across sub-volumes. The stepped pressure equilibrium code (SPEC) was developed to calculate MRxMHD equilibria numerically. In this work, we have extended the SPEC code to compute MRxMHD equilibria with field-aligned flow and rotation, following the theoretical development to incorporate cross-helicity and angular momentum constraints. The code has been verified for convergence and compared to a Grad–Shafranov solver in 2D. We apply our new tool to study the flow profile change before and after the sawtooth crash of a reversed-field pinch discharge, in which data of the parallel flow is available. We find the promising result that under the constraints of cross-helicity and angular momentum, the parallel flow profile in post-crash SPEC equilibrium is flat in the plasma core and the amplitude of the flow matches experimental observations. Lastly, we provide an example equilibrium with a 3D helical field structure as the favored lower energy state. This will be the first 3D numerical equilibrium in which the flow effects are self-consistently calculated.
Network protector units (NPUs) are crucial parts of the protection of secondary networks to effectively isolate faults occurring on the primary feeders. When a fault occurs on the primary feeder, there is a path of the fault current going through the service transformers that causes a negative flow of current on the NPU connected to the faulted feeder. Conventionally, NPUs rely on the direction of current with respect to the voltage to detect faults and make a correct trip decision. However, the conventional NPU logic does not allow the reverse power flow caused by distributed energy resources installed on secondary networks. The communication-assisted direct transfer trip logic for NPUs can be used to address this challenge. However, the communication-assisted scheme is exposed to some vulnerabilities arising from the disruption or corruption of the communicated data that can endanger the reliable operation of NPUs. This paper evaluates the impact of the malfunction of the communication system on the operation of communication-assisted NPU logic. To this end, the impact of packet modification and denial-of-service cyberattacks on the communication-assisted scheme are evaluated. The evaluation was performed using a hardware-in-the-loop (HIL) co-simulation testbed that includes both real-time power system and communication network digital simulators. This paper evaluates the impact of the cyberattacks for different fault scenarios and provides a list of recommendations to improve the reliability of communication-assisted NPU protection.
The accompanying poster to this short paper presents a combination of reverse mode AD and formal methods to enable efficient differentiation of (or backpropagation through) shared-memory parallel code. Compared to the state of the art, our approach can more often avoid the need for atomic updates or private data copies during the parallel derivative computation, even in the presence of unstructured or data-dependent data access patterns. This is achieved by gathering information about the memory access patterns from the input program, which is assumed to be correctly parallelized. This information is then used to build a model of assertions in a theorem prover, which can be used to check the safety of shared memory accesses during the parallel derivative computation
Here, we present an efficient level set method to track an arbitrary number of materials. The algorithm is optimal in the sense that it only needs to store a single unsigned distance-like function and a single integer indicator function, independent of the number of materials or distinct regions being tracked. Furthermore, for smooth velocity fields and smooth interface shape, arbitrarily high order solutions can be demonstrated. For interfaces that are or become kinked, the solution is limited to second-order convergence rates in the L 1 norm and first-order in the L ∞ norm.
Extrusion-based additive manufacturing has emerged as a powerful platform for designing shape-morphing materials through controlled orientation. However, existing approaches primarily rely on a single mode of flow-induced alignment, limiting orientation programmability. Herein, we present a direct-ink-writing approach for smectic liquid crystal elastics that exploits two distinct alignment modes within a single ink. The smectic ink exhibits shear- and temperature-dependent orientation switching, enabling molecular alignment either perpendicular or parallel to the print direction. Combined rheological, X-ray, and molecular dynamics analyses reveal that this alignment inversion arises from the preservation or collapse of smectic layers under flow. This reversible switching encodes both contractile and elongational actuation within individual filaments, greatly expanding the design freedom of printed liquid crystal elastomers. We demonstrate 2D and 3D structures with diverse programmed shape transformations, highlighting the potential of this platform for adaptive soft actuators and architected functional materials.