Search NASA⌕ Search

SEARCH · Search NASA

Results for “unsteady”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Mixing in Low Reynolds Number Reacting Impinging Jets in Crossflow

Previous efforts to model uranyl fluoride formation in an impinging jet gas reactor underpredicted spatial mixing and overpredicted chemical conversion into particulates. The previous fluid dynamics model was based on the solution of the Reynolds Averaged Navier Stokes equations. After simulating fluid dynamics, aerosol dynamics were superimposed onto CFD-simulated gas reactant species concentrations. The current work explores the influence of complex unsteady flow features on the overall flow physics and chemistry for a low Reynolds number, opposed flow, impinging jet gas reactor where there is a low Reynolds number cross flow. The objective of this study was to assess the impact of model formulation on scalar mixing and transport. Here, transient flow simulations were performed using Scale Resolving Simulations. Large-Eddy Simulations with the dynamic Smagorinsky turbulence model were performed along with simulations which directly resolved the flow. Average and root-mean-square (RMS) velocities and species concentrations were computed along with modeled and resolved turbulence kinetic energy (TKE), modeled turbulence dissipation, and modeled turbulent viscosity. Lagrangian flow tracers were also used to quantify species concentrations along path lines emanating from the jet tips. Transient simulation data were compared to results from RANS simulations using the k-ω shear stress transport (SST) model and Reynolds Stress Model (RSM). Transient simulations showed spatial mixing patterns which were more consistent with experimental data and helped elucidate the process of particle formation observed in experiments.

42 ENGINEERING↗
DOE OSTI2022

Optimum Injector Parameters for Thermoacoustic Stability in a Multi-Nozzle Can Combustion System

Abstract High-frequency transverse instabilities are an important concern in can combustor configurations. In these configurations which are typically operated with multiple injectors around a central injector, each injector is subjected to different parts of the acoustic mode shape and thus respond differently for the same instability mode. Recent work by the author has modeled the response of premixed flames to excitation by natural high-frequency transverse modes in a can combustor both in the center and outer nozzles. The stability of these acoustically non-compact flames was assessed using the Rayleigh criterion (Rayleigh Integral denoted as RI) and not the overall unsteady heat release as is the case for compact flames. Several key control parameters were studied, namely — flame angle, swirling strength, nozzle location. For non-axisymmetric modes such as the commonly occurring 1-T mode, both radial and azimuthal offsets of the nozzle location affected stability. The framework was applied to an optimization study to identify the optimal combination of parameters that minimizes RI for the different nozzles in the multi-nozzle system. In this study, a N-around-1 configuration was studied, and the results indicated that the different nozzles needed to be operated at different flame angles and swirl numbers to result in an overall minimum RI. However, the specific response of the different injectors was not considered. The helical mode distribution at each injector varies as we azimuthally go around the combustor’s injector distribution and thus the most amplified mode and the resulting flame response would be different. To minimize RI, it is important to determine the injector configurations that result in a hydrodynamic profile that minimizes the individual RI for each nozzle. The resulting relationship between the injector’s flow and local hydrodynamics can then be used in a hydrodynamics study of an individual injector so that the most optimal injector is chosen depending on its location in the combustor dump plane.

Acharya, Vishal↗
DOE OSTI2025

Corrosion of Rare-Earth Disilicate Environmental Barrier Coatings in Gas Turbine Environments—Part 2: Microstructure Model

A model for the volatization of ytterbium disilicate (YbDS) by water vapor is developed. The new model is based on the Deal & Grove theory of the thermal oxidation of silicon [1], generalized here to account for the multi-step, second-order chemical reactions for silica depletion, multi-component mixture effects, and the unsteady porous microstructural evolution of the corroded inner zone. With the new model, the fundamental kinetic rate and standard porous media parameters are determined using the atmospheric steam jet impingement experiment of Ridley & Opila [2] and compared to results from the original Deal & Grove model. Fixing the input parameters in each respective model, a validation study (blind prediction) is presented for the high Reynolds number (1 × 105) and high pressure (14.7 atm) burner rig experiment of Wan et al. [3].

Cha, Chong M [Rolls Royce]↗
DOE OSTI2023

Computational Fluid Dynamics Study of a Cross-Flow Marine Hydrokinetic Turbine and the Combined Influence of Struts and Helical Blades

A computational fluid dynamics study was performed for a cross-flow marine hydrokinetic turbine. The analysis was done in three dimensions and used the unsteady Reynolds-averaged Navier-Stokes solver in the commercial code STAR-CCM+. The base turbine configuration is the RivGen(R) Turbine, designed by the Ocean Renewable Power Company. A convergence and uncertainty analysis was performed for both the spatial and temporal discretization; this was done using the base configuration, which features support struts and helical foils. Both struts and helical blades introduce three-dimensional flow effects, influencing the complex flow phenomenon of dynamic stall. The study compares the relative impact of struts on power performance and blade loading for both helical and straight blades, and found that for this turbine the relative loss in power from struts was lower with helical blades.

CFD↗
DOE OSTI2024

High- and Mid-Fidelity Modeling Comparison for a Floating Marine Turbine System

There is a lack of suitable numerical tools, particularly open-source tools, that can be used for designing and optimizing marine turbine systems. The National Renewable Energy Laboratory has added features to their widely used mid-fidelity wind turbine modeling code, OpenFAST, to enable modeling of axial-flow marine turbines. This necessitated the addition of several physical effects relevant to marine turbines that are neglected for wind turbines. These include buoyancy, added mass and inertial loads, wave-current superposition, and changes to the coordinate systems. This updated version of OpenFAST allows for the modeling of both fixed and floating marine turbines at a speed comparable to real time. While efficient for large sets of load cases and design studies, mid-fidelity codes make simplifying assumptions that may impact their accuracy. High-fidelity computational fluid dynamics (CFD) simulations can capture more flow effects with fewer assumptions and provide detailed body pressure mapping and flow-field information. It is important to compare predictions between mid-fidelity and high-fidelity codes, both to verify the models and to understand the limitations. A floating marine turbine system was modeled both with OpenFAST and with the commercial CFD code STAR-CCM+. The CFD model used a three-dimensional unsteady Reynolds-averaged Navier-Stokes solver for a volume-of-fluid numerical wave and current tank. The blade-resolved simulations used the sliding-interface technique for the spinning rotor and an overset grid to accommodate the rigid-body motion of the floating system. The mooring system was modeled with a custom coupling of the CFD solver with the open-source code MoorDyn. This improves upon the existing quasi-static catenary solver in STAR-CCM+, which lacks seabed contact or line-to-line connections. Simulation results for a floating marine turbine are compared between OpenFAST and CFD, highlighting the capabilities of the mid-fidelity code and identifying the areas where a high-fidelity approach is needed.

CFD↗
DOE OSTI2025

Artificial to Spiking Neural Networks Conversion with Calibration in Scientific Machine Learning

Here, we introduce a method to convert physics-informed neural networks (PINNs), commonly used in scientific machine learning, to spiking neural networks (SNNs), which are expected to have higher energy efficiency compared to traditional artificial neural networks (ANNs). We first extend the calibration technique of SNNs to arbitrary activation functions beyond ReLU, making it more versatile, and we prove a theorem that ensures the effectiveness of the calibration. We successfully convert PINNs to SNNs, enabling computational efficiency for diverse regression tasks in solving multiple differential equations, including the unsteady Navier–Stokes equations. We demonstrate great gains in terms of overall efficiency, including separable PINNs (SPINNs), which accelerate the training process. Overall, this is the first work of this kind and the proposed method achieves relatively good accuracy with low spike rates.

PINN↗
DOE OSTI2025

Mapping Spatiotemporal Solvent Velocity from Measured Concentration Gradients in a Polarized Electrolyte

The electric-field induced motion of neutral species impedes the efficacy of electrochemical devices. By combining operando X-ray transmission measurements with continuum mechanics, we have developed a methodology for determining the velocity of neutral solvent molecules under an applied field. The X-ray transmission experiments were used to determine ion concentration profiles as a function of space and time in a polymer electrolyte. The unsteady state solvent mass balance equation was solved numerically with experimental concentration profiles to map spatiotemporal solvent velocities. We compare our experimentally derived results with predictions made with concentrated solution theory. We use the cation transference number as the only adjustable parameter to match experimental measurements of both concentration and solvent velocity. Our approach may be used to determine solvent velocity with any operando technique used to measure time-dependent ion concentration profiles.

Abdo, Emily E↗
DOE OSTI2020

Scalability of OpenFOAM Density-Based Solver with Runge–Kutta Temporal Discretization Scheme

Compressible density-based solvers are widely used in OpenFOAM, and the parallel scalability of these solvers is crucial for large-scale simulations. In this paper, we report our experiences with the scalability of OpenFOAM’s native rhoCentralFoam solver, and by making a small number of modifications to it, we show the degree to which the scalability of the solver can be improved. The main modification made is to replace the first-order accurate Euler scheme in rhoCentralFoam with a third-order accurate, four-stage Runge-Kutta or RK4 scheme for the time integration. The scaling test we used is the transonic flow over the ONERA M6 wing. This is a common validation test for compressible flows solvers in aerospace and other engineering applications. Numerical experiments show that our modified solver, referred to as rhoCentralRK4Foam, for the same spatial discretization, achieves as much as a 123.2% improvement in scalability over the rhoCentralFoam solver. As expected, the better time resolution of the Runge–Kutta scheme makes it more suitable for unsteady problems such as the Taylor–Green vortex decay where the new solver showed a 50% decrease in the overall time-to-solution compared to rhoCentralFoam to get to the final solution with the same numerical accuracy. Finally, the improved scalability can be traced to the improvement of the computation to communication ratio obtained by substituting the RK4 scheme in place of the Euler scheme. All numerical tests were conducted on a Cray XC40 parallel system, Theta, at Argonne National Laboratory.

Li, Sibo↗
DOE OSTI2026

UnsteadyKineticRotorDynamics.jl

SAND2026-22948O UnsteadyKineticRotorDynamics.jl is a Julia-based tool for time-domain unsteady rotor dynamics in wind and hydrokinetic rotors. The tool combines steady CCBlade aerodynamic calculations with low-pass dynamic state updates to model transient rotor response. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Michelen Strofer, Carlos [Sandia National Lab. (SN↗
DOE OSTI2022

Unified Wind-Wave Growth and Spectrum Functions for All Water Depths: Field Observations and Model Results

Abstract Wind-wave development is governed by the fetch- or duration-limited growth principle that is expressed as a pair of similarity functions relating the dimensionless elevation variance (wave energy) and spectral peak frequency to fetch or duration. Combining the pair of similarity functions, the fetch or duration variable can be removed to form a dimensionless function of elevation variance and spectral peak frequency, which is interpreted as the wave energy evolution with wave age. The relationship is initially developed for quasi-neural stability and quasi-steady wind forcing conditions. Further analyses show that the same fetch, duration, and wave-age similarity functions are applicable to unsteady wind forcing conditions, including rapidly accelerating and decelerating mountain gap wind episodes and tropical cyclone (TC) wind fields. Here it is shown that with the dimensionless frequency converted to dimensionless wavenumber using the surface wave dispersion relationship, the same similarity function is applicable in all water depths. Field data collected in shallow to deep waters and mild to TC wind conditions and synthetic data generated by spectrum model computations are assembled to illustrate the applicability. For the simulation work, the finite-depth wind-wave spectrum model and its shoaling function are formulated for variable spectral slopes. Given wind speed, wave age, and water depth, the measured and spectrum-computed significant wave heights and the associated growth parameters are in good agreement in forcing conditions from mild to TC winds and in all depths from deep ocean to shallow lake. Significance Statement This paper presents a growth function and spectrum model to describe wind-wave development in all water depths. Their applicability covers a wide range of wind forcing conditions including steady, accelerating, decelerating, and tropical cyclone events. Support for the unified spectrum model and growth function is presented with field observations and numerical computations.

Hwang, Paul A.↗
DOE OSTI2020

Multi-line FLEET by imaging periodic masks

A simple linear configuration for multi-line FLEET velocimetry is used for the first time to image an overexpanded unsteady supersonic jet. The FLEET lines are spaced 0.5-1.0 mm apart and up to six lines can simultaneously be used to visualize the flowfield. These lines are created using periodic masks, despite the mask blocking 25 — 30% of the 10 mJ incident beam. Maps of mean single-component velocity in the direction along the principal flow axis, and turbulence intensity in that same direction, are created using multi-line FLEET, and computed velocities agree well with that obtained from single-line (traditional) FLEET. Compared to traditional FLEET, multi-line FLEET offers increased simultaneous spatial coverage and the ability to produce spatial correlations in the streamwise direction. This FLEET permutation is especially well-suited for short-duration test facilities

42 ENGINEERING↗
DOE OSTI2020

Long-term persistence of infectious Zika virus: Inflammation and behavioral sequela in mice

The neurodevelopmental defects associated with ZIKV infections early in pregnancy are well documented, however the potential defects and long-term consequences associated with milder infections in late pregnancy and perinatal period are less well understood. To model these, we challenged 1 day old (P1) immunocompetent C57BL/6 mice with ZIKV. The animals developed a transient neurological syndrome including unsteady gait, kinetic tremors, severe ataxia and seizures 10–15 days post-infection (dpi) but symptoms subsided after a week, and most animals survived. Despite apparent recovery, MRI of convalescent mice show reduced cerebellar volume that correlates with altered coordination and motor function as well as hyperactivity and impulsivity. Persistent mRNA levels of pro-inflammatory genes including Cd80 , Il-1α , and Ifn-γ together with Cd3 , Cd8 and perforin ( PrfA) , suggested persistence of low-grade inflammation. Surprisingly, the brain parenchyma of convalescent mice harbor multiple small discrete foci with viral antigen, active apoptotic processes in neurons, and cellular infiltrates, surrounded by activated astrocytes and microglia as late as 1-year post-infection. Detection of negative-sense strand viral RNA and isolation of infectious virus derived from these convalescent mice by blinded passage in Vero cells confirmed long-term persistence of replicating ZIKV in CNS of convalescent mice. Although the infection appears to persist in defined reservoirs within CNS, the resulting inflammation could increase the risk of neurodegenerative disorders. This raises concern regarding possible long-term effects in asymptomatic children exposed to the virus and suggests that long-term neurological and behavioral monitoring as well as anti-viral treatment to clear virus from the CNS may be useful in patients exposed to ZIKV at an early age.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Validation of simulated wave energy converter responses to focused waves

Focused waves are increasingly used as a design method to predict extreme loads for offshore structures, such as wave energy converters (WECs). However, the ability of computational methods to accurately simulate the prescribed focused waves and resulting wave structure interactions is not well validated. This study uses two computational fluid dynamic methods, WEC-Sim and STAR-CCM+, to evaluate the responses of two taut-moored WEC-like buoys to three different focused waves. WEC-Sim is a computationally efficient, mid-fidelity model, based on linearised potential flow theory. In contrast, STAR-CCM+ is a high-fidelity, three-dimensional, unsteady, Reynolds-averaged Navier–Stokes-based model. The two buoy geometries considered are a hemispherical-bottom cylinder and a truncated cylinder with a cylindrical moon pool. Three prescribed focused waves are generated in both codes, and the simulated displacements and mooring loads are obtained for each of the hydrodynamic bodies. The resulting STAR-CCM+ generated focused waves have approximately the same accuracy, in comparison with the analytic solution, as experimentally generated focused waves. The WEC-Sim simulated responses, when using viscous coefficients tuned to the measured displacements, are, on average, within 11.5% of experimentally measured values, whereas the STAR-CCM+ simulated responses are, on average, within 14.2% of the experimental values.

16 TIDAL AND WAVE POWER↗
DOE OSTI2020

Pre-Sinkhole Seismicity at the Napoleonville Salt Dome: Implications for Local Seismic Monitoring of Underground Caverns

The formation of a sinkhole at the Napoleonville salt dome located near Bayou Corne in Assumption Parish, Louisiana, was preceded by an increase in seismicity, as observed 10 km away by a continuously recording broadband seismic sensor. The seismicity rate in the days prior to the sinkhole formation contained hundreds of long-period events, some of which could be detected using a power detector (short-term to long-term energy ratio) tuned to the duration and frequency content of the seismicity. A characteristic of some of these events is a resonance near 0.4 Hz, which could be due to either an acoustic resonance within the brine-filled cavity or possibly due to unsteady fluid flow through a crack. In this work, a more exhaustive catalog was obtained using some of the events detected by the power detector as templates in a correlation detector. Furthermore, the full pre-sinkhole sequence of detections could offer lessons in remote monitoring of underground caverns.

58 GEOSCIENCES↗
DOE OSTI2020

Analysis of the performance of different culvert boundary conditions in 2D shallow flow models

Culverts allow roads to safely traverse small streams or drainage ditches, and their proper design is critical to ensure a safe and reliable transportation network. A correct modelization of these hydraulic structures becomes crucial in the assessment of flood footprints or discharge peak estimation in a risk evaluation plan. The question of how to include culverts comes up frequently when assembling a hydraulic model that requires the presence of as many singular elements as possible. In this work, three different culvert integrations with the surface domain are studied and compared in the context of a 2D shallow water (SW) model. All of them are based on the Federal Highway Administration (FHWA) formulation for the culvert discharge estimation but differ in complexity and in the interaction with the numerical model for surface flow, some of them as internal boundary conditions. Several steady and unsteady validation test cases are presented and the numerical results are compared with the predictions from HEC-RAS 1D and HY-8 software. Furthermore, the culvert area, shape and their sensitivity to the 2D computational mesh is also analyzed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Understanding Transient Combustion Phenomena in Low-NO x Gas Turbines

The goal of this work was to understand, quantify, and eventually predict combustion instability during operational transients in low-NO x gas turbine combustors. Transient operation is defined as a mode where one or more key operating parameters, including equivalence ratio, fuel composition, or fuel splitting between nozzles in a combustor, are varying in time. The timescale over which this variation occurs will be referred to as the transient timescale, the magnitude of the change in the transient parameter will be referred to as the transient amplitude, and the direction of the variation will be referred to as the transient direction. Transients occur regularly in power generation gas turbine operation. For example, load-following gas turbines vary load to meet power demands while accounting for variation in unsteady grid contributors like wind and solar power. Operators achieve this by varying equivalence ratio and fuel splitting inside the combustors to produce more or less power. Another example of transient operation is during fuel-switching operations, a particular issue for turbines operating with syngas fuels in a co-generation facility. In this case, the composition of the fuel stream may vary in time either as a result of fluctuations in the pipeline, or in the case of a co-generation plant, variations in feedstock or gasification processes. This fluctuation in fuel composition can have significant implications for combustion processes, particularly with high-hydrogen fuels. A small number of studies have investigated transient behavior in gas turbine combustors, particularly in the context of flame blow-off; here, investigators have typically focused on issues of transient amplitude, though without the depth and physics-based experimental design that is proposed in this program. To our knowledge, there have been no studies that have quantified the impact of the transient timescale on combustion. The work helped us understand the impact of several types of transients on combustion stability. Both static and dynamic combustion stability will be considered. The transients will be quantified using their amplitude, direction, and timescale, and we will consider transients in equivalence ratio, fuel composition, and fuel splitting. There were two main objectives of this work: 1. Provide a fundamental understanding of combustion stability during transient operation using data obtained from a number of high-speed diagnostic techniques on experimental configurations that capture important gas turbine combustor features and operating conditions. This goal was met in conjunction with our industrial partner, GE Global Research, with whom we worked closely throughout the program. 2. Develop a stability prediction framework that will allow for the prediction, and possible control, of instabilities as a function of the transient amplitude, direction, and timescale. This input/output model should be physics-based and allow for screening of operating conditions and prediction of unstable behavior at the end of a transient by monitoring combustor health and precursor signals during operation.

03 NATURAL GAS↗
DOE OSTI2020

Collaborative Research: The Evolution of Magnetic Skeletons During 3D Reconnection. Final report

We investigated how the plasmoid instability affects the evolution of topological features of a magnetic field and the process of magnetic reconnection with three sub-projects. (i) We extensively studied the criterion for plasmoid-mediated current sheet disruption and the onset of fast reconnection. Scaling relations with respect to plasma parameters have been established, using both direct numerical simulations and analytic theories, for the disruption time, the current sheet width, the instability linear growth rate, and the dominant mode wavenumber. (ii) We applied the theories of plasmoid-mediated current sheet disruption to magnetohydrodynamic turbulence and predicted a steepening of the turbulence energy spectrum. This steepening has been confirmed by direct numerical simulations. (iii) We performed numerical simulations of ultraviolet burst events in the solar atmosphere with a spine-dome magnetic topology. The plasmoid instability leads to unsteady reconnection, which is consistent with solar observations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

ATDM AD projects to ECP Report (Q3 FY2020)

The SPARC team has completed much of the workflow development needed for the coupled and uncoupled analyses needed to support our flight test validation work. The SPARC team has continued to make progress on development activities to support unsteady, full reentry vehicle aero analysis, which have focused on turbulence modeling, uniform mesh refinement, in-situ visualization (with the Paraview/Catalyst team), and performance analysis. The SPARC team is making solid progress towards our Q4 goal of demonstrating an improved Wall Modeled Large Eddy Simulation (WMLES) capability for standard geometries using high-order finite difference, discontinuous Galerkin, and low-dissipation finite volume methods. The SPARC team has met the minimum completion criteria for parts of the L1, including documenting performance and scaling of SPARC on Trinity, Sierra and Astra, and performing runs for the flight test validation case.

97 MATHEMATICS AND COMPUTING↗