SEARCH · Engineering Papers

Results for “energy use behavior”

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.

SourceResearch connectionFrom yearThrough year

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

At least 1,315 records · Page 73

All sources NASA NTRS DOE OSTI

DOE OSTI2025

The Effect of Flow on CO2 Corrosion of Self-Healing Metallic Coatings

Internal corrosion is an issue that affects natural gas pipelines, a significant part of the United States’ energy infrastructure. Over time, this corrosion has worn away longstanding pipelines due to the original construction materials used for the lines and impurities in the gas and liquid streams flowing through them. The main impurities in natural gas are H2O, CO2, H2S, and O2. One of the leading causes of corrosion is the CO2 dissolved in water, giving rise to carbonic acid formation, which can further dissolve the steel pipe. The National Energy Technology Laboratory (NETL) has been studying different solutions to this problem. One potential answer is a self-healing sacrificial metallic coating applied by a novel cold spray technique. This study explored the impact of flow on corrosion with this coating and carbon steel when exposed to a saturated CO2 environment by simulating the pipeline flow profile in a small-scale lab setting. The samples were evaluated using multiple electrochemical techniques that found corrosion rates and were backed up with surface analysis to conclude the behavior of the corrosion mechanisms. The cold spray coating exhibited steel corrosion protection under flow conditions. This study has significantly advanced our understanding of CO2 corrosion, particularly in the context of natural gas pipelines and coating design.

carbon steel↗

NASA NTRS1988

An experimental and theoretical study of the ice accretion process during artificial and natural icing conditions

Real-time measurements of ice growth during artificial and natural icing conditions were conducted using an ultrasonic pulse-echo technique. This technique allows ice thickness to be measured with an accuracy of + or - 0.5 mm; in addition, the ultrasonic signal characteristics may be used to detect the presence of liquid on the ice surface and hence discern wet and dry ice growth behavior. Ice growth was measured on the stagnation line of a cylinder exposed to artificial icing conditions in the NASA Lewis Icing Research Tunnel (IRT), and similarly for a cylinder exposed in flight to natural icing conditions. Ice thickness was observed to increase approximately linearly with exposure time during the initial icing period. The ice accretion rate was found to vary with cloud temperature during wet ice growth, and liquid runback from the stagnation region was inferred. A steady-state energy balance model for the icing surface was used to compare heat transfer characteristics for IRT and natural icing conditions. Ultrasonic measurements of wet and dry ice growth observed in the IRT and in flight were compared with icing regimes predicted by a series of heat transfer coefficients. The heat transfer magnitude was generally inferred to be higher for the IRT than for the natural icing conditions encountered in flight. An apparent variation in the heat transfer magnitude was also observed for flights conducted through different natural icing-cloud formations.

Kirby, Mark S.↗

DOE OSTI2020

Understanding ICF hohlraums using NIF gated laser-entrance-hole images

The newly available ns-gated laser-entrance-hole (LEH) imager on the National Ignition Facility provides routine, non-perturbative measurements of the x-ray emission from laser-heated plasmas inside the hohlraum as viewed at 19° to the hohlraum axis through one of its LEHs. Multiple images are acquired for a series of times and filter-selected x-ray energy bands within a single shot. The images provide time dependent data on phenomena including the effective radius of the LEH, the length of the gold-plasma “bubble” evolving off the interior wall surface heated by the outer beams, the evolving radius of the x-ray heated hohlraum wall, and the radius of the ablation front of the fusion capsule. Here, these measurements are explained and illustrated with sample data. These techniques are then applied to understand hohlraum behavior as a function of gas fill. For hohlraums with helium gas fill densities of 0.15 to 0.30 mg/cm 3 , synthetic images computed from simulations agree well with experimental gated LEH images when an inhibited heat transport model [Jones et al., Phys. Plasmas 24, 056312 (2017)] is used. This model can be adjusted to reproduce the expansion rate of the laser-heated plasma bubble in such a way as to improve agreement with the images. At the higher 0.6 mg/cc gas fill, the experimental images show more pronounced 3D features, resulting in slightly less good agreement with the 2D simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Inverse reinforcement learning control for building energy management

Reinforcement learning (RL) based control is widely considered a promising approach in building automation and control as it has demonstrated the potential to deal with complex objectives in adjacent domains like robotics, autonomous vehicles, gaming applications, and advertisement recommendations. When applied to any environment, model-free RL learns to improve its control performance over time without requiring a control model, by receiving and then analyzing feedback from the building environment after each control action. Operational objectives are becoming increasingly complex through the simultaneous consideration of thermal comfort, carbon emissions, grid services, and indoor air quality. In this context, conventional rule-based control approaches are proving sub-optimal, mostly heuristic, and inadequate. The model-free and self learning nature of RL appears promising and attractive as it may address the scalability issues associated with advanced control approaches. However, it suffers from long training times and unstable control behavior during the early stages of its learning process, which makes it unsuitable to be applied directly to buildings. This paper addresses these issues using an inverse reinforcement learning approach (IRL), a technique utilized to learn the objective of a controller agent which is considered an expert in its respective domain. Here, we consider a rule-based control as the expert demonstrator. IRL is different from a direct imitation (i.e., direct mapping of states to actions) of control actions as it tries to find the underlying intent of an expert's policy, providing the controller with a better-generalized policy for unseen states or environments with slightly different dynamics. This approach propels the RL controller's policy to levels similar to or better than that of a rule-based policy before it starts learning by interacting with the building. This makes RL for building energy management applications more practical as it prevents the erratic and exploratory behavior in the initial training period, simultaneously speeding up the learning process when compared to applying an untrained RL agent directly to a building environment.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2020

Nanostructures for Electrical Energy Storage (NEES) (2020 Final Technical Report)

Nanostructures for Electrical Energy Storage (NEES, www.efrc.umd.edu) was an Energy Frontier Research Center supported by the DOE Office of Science, Basic Energy Sciences, from 8/1/2009 to 7/31/2020. Led by the University of Maryland, NEES enjoyed extensive collaborations with its funded partners, including two DOE Laboratories and six universities. The NEES vision has been to reveal a set of scientific insights and design principles that can underpin a next-generation electrical energy storage approach, building on advances in nanoscale science and technology to achieve simultaneous high power and high energy over extended charge/discharge cycling. The vision is motivated by the recognition that scaling into the nano regime opens the door to new physical phenomena and that the tools enlisted in nanoscale research provide major new opportunities for the synthesis not only of materials at molecular scale but for structures at nano scale and above. NEES has translated this vision into its research program based on two observations. First, while the behavior of ions and electrons in electrolytes and in electrode materials is crucial to electrical energy storage (or more appropriately electrochemical energy storage), it is the transport of ion and electron charge between different structural components of a storage device that ultimately determine its performance. With it well recognized that the choice of electrode materials typically constrain ion transport kinetics as well as maximum ion concentration, the search for better electrode materials has been a primary driver of battery research. At the same time the synthesis of electrodes is typically based on aggregation of particles with varying size, shape, and orientation in the electrode. Together with the presence of additional materials to impart electrical conductivity and cohesion to the composite electrode, change in electrode materials is necessarily accompanied by structural changes at the nano/micro scale that are difficult to categorize and manage. From the beginning, NEES’ vision has been to create and study simpler, highly controlled spatial arrangements of known materials as battery components (electrodes, current collectors, and electrolyte) and to understand how design and structure above the molecular scale determines the energy storage performance available from known materials. Second, advances in nanoscience dramatically expanded the portfolio of synthesis methods, structural motifs, and new phenomena available for research. Some of these gave rapid access to new building blocks at the deep nanoscale (e.g., carbon nanotubes grown by self-assembly, nanoscale arrays formed by electrochemical self-alignment, monolayer films controlled by self-limiting reaction). Such advances served as the enabler for the NEES vision to be pursued experimentally through study of 3D structures created and controlled at the nano, micro, and meso scales. Here, we use meso as in the BES MESO Report, implying not only intermediate or varying length scales, but very much the way behavior is influenced by other factors including aggregation of nanocomponents at different densities and spatial configurations, statistical variations in the aggregates, hierarchical architectures in which they can be assembled, or local 3D configurations that result from the architectures. Over its life cycle, NEES has pursued two overarching goals: (1) to understand the scientific fundamentals of electrochemical storage from the nanoscale to the mesoscale; and (2) to create and learn from innovative, controlled, heterogeneous nanostructures, where such nanostructures can enable the first goal and serve as models for future paradigms in energy storage. Specific goals have included: Synthesize heterogeneous nanostructures comprised of multiple materials arranged in controlled fashion and characterize their behavior; Demonstrate and elucidate design principles for achieving simultaneous high power and high energy; Develop materials processes which enable precision control of thin layers and 3D structures; Investigate the impact of artificial interphases on electrode stability during ion insertion/deinsertion; Create dense arrays of nanostructures to understand how the architecture of these assemblies, along with nanostructure design, influences energy storage behavior at the mesoscale; Identify and understand the consequences of nanoconfinement and local inhomogeneities in 3D mesoscale arrays; Develop and apply computational models to stimulate, guide and interpret experiments.

25 ENERGY STORAGE↗

DOE OSTI2021

Process model for multilayer slide coating of polymer electrolyte membrane fuel cells

Slide coating is a precision method suitable for depositing multiple liquid-film layers simultaneously. Originally developed in the photographic film industry, it has been deployed for manufacturing of other products that benefit from multilayer coatings. One emerging application is the manufacture of polymer electrolyte membrane fuel cells (PEMFCs), which are used to produce electricity through electrochemical reactions of hydrogen and oxygen gas. The membrane-electrode assembly (MEA), in which key electrochemical reactions occur, consists of three layers that are typically deposited separately in serial fashion and then laminated together to form the three-layer MEA, i.e., three sequential steps of coat and dry. Adapting the process to simultaneous, multilayer slide coating of all three layers will save equipment cost and space while minimizing possible exposure to contaminants during transition between the steps. We are developing a three-layer slide coating model to aid the manufacturing process design of PEMFC. The model accounts for rheology of each layer, which typically exhibit shear thinning behavior. Model predictions are used to investigate simultaneous coatability of catalyst inks and to determine the best layer-by-layer ink selection.

25 ENERGY STORAGE↗

DOE OSTI2022

Optimization and Commercialization of the Juvenile Eel/Lamprey Acoustic Transmitter and Micro-battery (Abstract)

Enhance studies to track in 3D and sub-meter accuracy the movements of sensitive species and early life stages of fish to advance understanding of migration timing and behaviors, habitat use, and survival rates, resulting in more informed management decisions regarding new and existing hydroelectric facilities and better designs of new hydropower systems that minimize or avoid environmental impacts.

13 HYDRO ENERGY↗

DOE OSTI2025

Strategies for the in-orbit gain tracking using the modulated X-ray sources for the Resolve microcalorimeter spectrometer on the X-ray Imaging and Spectroscopy Mission

Accurate and precise correction of the gain drift is the key to achieve the required energy resolution of the Resolve microcalorimeter spectrometer on the X-ray Imaging and Spectroscopy Mission (XRISM). Therefore, Resolve is equipped with highly configurable X-ray sources called the modulated X-ray source (MXS). The pulsed nature allows us to separate calibration and astrophysical X-rays by time interval selections. However, undesirable characteristics of the MXS, such as the afterglow X-rays, restrict the allowed configuration range. Moreover, the nonlinear and discontinuous behaviors of the calibration line count rate make the determination of the optimal setting highly complex. The MXS count rate model has been established using measurements in the spacecraft thermal vacuum test with the flight detector and MXS. A trade-off study using the model enables us to choose a few settings for the continuous use of the MXS optimized for different ranges of target gain tracking intervals. An alternative approach, where the MXS is used only intermittently, has also been developed and implemented. This new mode enables us to reconstruct the drift without having most of the undesirable effects in science data. This also forms the basis of the gain tracking under the current Resolve configuration with the closed gate valve.

Astronomy and AstroPhysics↗

DOE OSTI2026

Evaluating hydrogen permeation through medium-density and high-density polyethylene pipes to enable a hydrogen-compatible infrastructure

As the United States and the global community pursue energy abundance and resilience, hydrogen is becoming a critical energy carrier. The future of hydrogen transport depends on one thing: verifying that today’s pipelines are hydrogen ready. This study investigates the permeation behavior of medium-density polyethylene (MDPE) and high-density polyethylene (HDPE) pipelines under pure hydrogen gas environment at various pressures and temperatures. Using thermal desorption analysis (TDA) and direct pipe permeability measurements, hydrogen loss rates are quantified, revealing minimal permeation losses (< 10.5 g/km/day) under realistic operating conditions. Effects of temperature and pressure with Standard Dimension Ratios (SDR), were evaluated to determine their influence on permeation rates of MDPE and HDPE pipelines. In addition, ex-situ density and degree of crystallinity (DOC) of MDPE and HDPEs after hydrogen exposure were also investigated to understand pressure-dependency of polymer properties. The findings provide essential data for assessing polymer-based pipeline materials for safe and reliable hydrogen transport.

08 HYDROGEN↗

DOE OSTI2026

Effects of Concentration, Salinity and Temperature on the Conformations of Zwitterionic Poly(2-Vinylpyridine- N -Oxide) Chains in Semidilute Solutions Probed by Small-Angle X-Ray and Neutron Scattering

Polyzwitterions, composed of repeating units with equal numbers of anionic and cationic groups, have drawn considerable interest for applications ranging from antimicrobial coatings and antifouling membranes to low-friction interfaces, energy storage media, and actuators. A characteristic feature of many polyzwitterions is the anti-polyelectrolyte effect salt-induced chain expansion in aqueous solution a phenomenon often implicated in their functional behavior. In this study, we investigated the conformational behavior of poly(2-vinylpyridine-N-oxide) (P2VPNO) using smallangle X-ray and neutron scattering (SAXS and SANS). Small-angle X-ray scattering (SAXS) revealed that, in salt-free aqueous solution, P2VPNO adopts an expanded wormlike conformation stabilized by hydration. With increasing concentration, the chains contract due to screening of intrachain excluded-volume interactions, qualitatively consistent with de Gennes’ scaling predictions for neutral polymers. Small-angle neutron scattering (SANS) further demonstrated that P2VPNO exhibits the characteristic anti-polyelectrolyte response to added salt, as observed in many other polyzwitterions. At elevated temperatures, chain flexibility increases, leading to a shorter Kuhn length and reduced radius of gyration. Notably, however, salt-induced chain expansion persists, indicating that the anti-polyelectrolyte effect remains operative under elevated thermal conditions. These findings provide the first experimental evidence of scaling behavior in polyzwitterions, as well as the first observation of their altered antipolyelectrolyte response at elevated temperatures, offering new insights into the solution physics of this important class of polymers.

Conformation↗

DOE OSTI2026

Self-Assembly Driven Microlithography via Near-Infrared Light Activation

Current vat photopolymerization (VP) relies on UV or visible light to start the rapid crosslinking of liquid photocurable resins into 3D-printed structures. Here, we develop a self-assembly-driven photopatterning approach to photothermally generate polymeric solids by combining thermoplasmonic nanoparticles and thermoresponsive polymers, in which near-infrared (NIR) light activates thermoplasmonic heating of nanoparticles, triggering the irreversible self-assembly of thermoresponsive polymers into insoluble mesoglobules. A small amount of thermal initiator and crosslinker leads to irreversible self-assembly of polymer nanocomposites. NIR light offers deeper penetration and reduced scattering compared to UV, enabling more uniform curing of thicker or filled materials and expanded process control for composites or opaque systems. Thermoplasmonic heat generation is achieved using surface-modified gold nanorods (AuNRs) with a longitudinal localized surface plasmon resonance peak in the NIR region. Key variables such as polymer composition, molecular weight, physical interactions at the nanoparticle–polymer interface, which can be tuned by surface functionalization, AuNR concentration, and pH can be used to tailor the assembly behavior of these systems, including photothermal effect, flocculation, and cloud point temperature, and the mechanical properties of the final structures. Collectively, these results highlight a platform for photothermally-driven microlithography of polymer solids with diverse, tunable macroscopic properties, enabled by low-energy NIR light-activated self-assembly.

36 MATERIALS SCIENCE↗

DOE OSTI2022

From Fluid Flow to Coupled Processes in Fractured Rock: Recent Advances and New Frontiers

Abstract Quantitative predictions of natural and induced phenomena in fractured rock is one of the great challenges in the Earth and Energy Sciences with far‐reaching economic and environmental impacts. Fractures occupy a very small volume of a subsurface formation but often dominate fluid flow, solute transport and mechanical deformation behavior. They play a central role in CO 2 sequestration, nuclear waste disposal, hydrogen storage, geothermal energy production, nuclear nonproliferation, and hydrocarbon extraction. These applications require predictions of fracture‐dependent quantities of interest such as CO 2 leakage rate, hydrocarbon production, radionuclide plume migration, and seismicity; to be useful, these predictions must account for uncertainty inherent in subsurface systems. Here, we review recent advances in fractured rock research covering field‐ and laboratory‐scale experimentation, numerical simulations, and uncertainty quantification. We discuss how these have greatly improved the fundamental understanding of fractures and one's ability to predict flow and transport in fractured systems. Dedicated field sites provide quantitative measurements of fracture flow that can be used to identify dominant coupled processes and to validate models. Laboratory‐scale experiments fill critical knowledge gaps by providing direct observations and measurements of fracture geometry and flow under controlled conditions that cannot be obtained in the field. Physics‐based simulation of flow and transport provide a bridge in understanding between controlled simple laboratory experiments and the massively complex field‐scale fracture systems. Finally, we review the use of machine learning‐based emulators to rapidly investigate different fracture property scenarios and accelerate physics‐based models by orders of magnitude to enable uncertainty quantification and near real‐time analysis.

58 GEOSCIENCES↗

DOE OSTI2023

In situ X-ray area detector flat-field correction at an operating photon energy without flat illumination

Flat-field calibration of X-ray area detectors is a challenge due to the inability to generate an X-ray flat-field at the selected photon energy the beamline operates at, which has a strong influence on the measurement behavior of the detector. A method is presented in which a simulated flat-field correction is calculated without flat-field measurements. Instead, a series of quick scattering measurements from an amorphous scatterer is used to calculate a flat-field response. The ability to rapidly obtain a flat-field response allows for recalibration of an X-ray detector as needed without significant expenditure of either time or effort. Area detectors on the beamlines used, such as the Pilatus 2M CdTe, PE XRD1621 and Varex XRD 4343CT, were found to have detector responses that drift slightly over timescales of several weeks or after exposure to high photon flux, suggesting the need to more frequently recalibrate with a new flat-field correction map.

47 OTHER INSTRUMENTATION↗

DOE OSTI2025

Micromechanical Surrogate Machine Learning Model for Creep Deformation Modeling

Process variability during the manufacture of gas turbine engine hot section components can significantly affect the material’s resulting microstructure. In casting, for instance, geometric variation within a component (thin sections versus thick sections, radial location) influences cooling rates and the resulting grain size. The high temperature creep response is known to be sensitive to grain size owing to a diffusional creep mechanism which occurs more readily along grain boundaries. Microstructural variation correspondingly drives mechanical behavior which propagates into component scale performance uncertainty. These factors are essential when planning inspection, maintenance, and repair strategies within a reliability framework. These benefits provide opportunities to increase overall energy efficiency through refined margins. Critically, there is an opportunity to bolster existing data-driven reliability models using physics-driven process-structure-property relations. Here we present recent work establishing a framework for evaluating the probabilistic creep performance of high-temperature materials. A novel microstructure-sensitive crystal plasticity finite element model is established that captures both grain boundary and crystallographic deformation effects. The computationally expensive physics model is calibrated using a statistical approach and this high-fidelity model is subsequently used to train a computationally efficient machine learning surrogate model. The surrogate model is essential for sampling a large ensemble of simulated structure-property pair results. The ensemble data are then mined to extract salient trends to be incorporated into a microstructure-sensitive reliability model. The proposed approach represents a novel way to capture microstructure-sensitive trends from physics-based models within a modern reliability framework.

Fernandez-Zelaia, Patxi [ORNL]↗

DOE OSTI2024

Bioblendstocks to Optimize Mixing Controlled Compression Ignition (MCCI) Engines

In this project, a team of researchers from the University of Massachusetts Lowell, the University of Maine, and Mainstream Engineering developed an integrated process for the product of bioblendstocks to optimize mixing controlled compression ignition (MCCI) engines. The objective was to improve the energy density, sooting propensity, and cetane number of base diesel fuel while maintaining cold weather behavior. The process converts woody biomass (e.g. sawmill residues) into bio-oil through selective fast pyrolysis; the bio-oil is then selectively upgraded to form selectively oxygenated, minimally-branched hydrocarbons using non-noble metal catalysts in combination with metal-catalyzed hydrogenation. Advanced predictive models, in conjunction with existing property databases, and experimental testing are used to evaluate overall bioblendstock properties and their impact on base diesel fuel. An iterative, targeted upgrading approach was implemented to optimize the proposed bioblendstock’s properties. Assessment methodologies included techno-economic analysis, life-cycle assessment, property testing, and engine testing. Ultimately, the project team successfully produced a viable bioblendstock while identifying critical process points related to scale-up efforts. It was found that producing pyrolysis oils at 500 degrees C and with pine particle sizes of 1-2 mm led to bio-oil with a higher yield (of approximately 45 wt%) and rich amounts of aromatic alcohols. The resultant pyrolysis oil was then upgraded using a sequence of mild hydrotreating, followed by catalytic etherification and esterification, followed by another final mild hydrotreating to produce a blendstock containing saturated species with a limited, but non-zero, amount of oxygen. The aromatic alcohols produced by pyrolysis were especially helpful in this regard, as the resulting bicycloethers and derivatives exhibited high cetane numbers. While most bulk properties of the bioblendstock met or exceeded targeted thresholds, viscosity and cloud point notably fell outside the expected range; this could be addressed by blending limits and/or through the use of additives that are commonplace in current refinding practices. Identification of a bioblendstock that can be produced economically at scale while improving the performance and emissions characteristics of internal combustion engines positively affects the economy by boosting domestic fuel production and the environment by decreasing harmful emissions and increasing efficiency.

09 BIOMASS FUELS↗

DOE OSTI2025

Utilization of Existing Pipelines in Hydrogen Transport: Literature Review Report

This report critically reviews the flow behavior of hydrogen-natural gas (H 2 -NG) mixtures in pipelines and examines the critical factors of hydrogen integration into existing natural gas infrastructure. It addresses the choking behavior characterized by velocity increase and pressure drop, as well as the effects of flow restrictions and pressure losses during hydrogen transport. Computational and analytical models are used to investigate these effects, and their effects on thermodynamic properties and system performance are evaluated. The study also reviews the energy efficiency and flow dynamics of hydrogen and methane-hydrogen mixtures and optimizes the hydrogen flow rate. In addition, the effects of these mixtures on the flow characteristics are discussed in detail, with special emphasis on the compressibility factor (z factor) and fluid properties based on equations of state for hydrogen-natural gas mixtures. The study also analyzes the mixture ratios and highlights the thermophysical properties, flow dynamics, and hydrogen-blended natural gas application potential. These investigations assess flow stability, material interactions, and operational feasibility of transporting hydrogen mixtures through natural gas pipelines, which contribute to developing sustainable and efficient energy systems.

08 HYDROGEN↗

NASA NTRS1987

Transient radiative cooling of a droplet-filled layer

A proposed lightweight radiator system for waste heat dissipation in space would eject streams of coolant in the form of small, hot liquid droplets. The droplets would lose radiative energy by direct exposure to the very low-temperature environment of space, and would then be collected for reuse. The cooling behavior of a layer composed of many small droplets was studied by numerical solution of the radiative integral equations. Since there is mutual interference for radiative energy dissipation, an array droplet will cool more slowly than if each drop is exposed individually. Since liquid metal droplets may be used, the study includes results for conditions with high scattering. For optically thin regions, especially with high scattering, the temperature distribution is sufficiently uniform that the cooling can be computed using the approximation of a constant layer emittance. For optically thick layers starting at uniform temperature, the temperature distributions become nonuniform with time. It was found that the cooling process goes through a starting transient; a constant emittance condition is then achieved where the emittance is lower than that for a layer at uniform temperature.

Siegel, R.↗

NASA NTRS1995

Pressure updating methods for the steady-state fluid equations

We consider the steady state equations for a compressible fluid. Since we wish to solve for a range of speeds we must consider the equations in conservation form. For transonic speeds these equations are of mixed type. Hence, the usual approach is to add time derivatives to the steady state equations and then march these equations in time. One then adds a time derivative of the density to the continuity equation, a derivative of the momentum to the momentum equation and a derivative of the total energy to the energy equation. This choice is dictated by the time consistent equations. However, since we are only interested in the steady state this is not necessary. Thus we shall consider the possibility of adding a time derivative of the pressure to the continuity equation and similar modifications for the energy equation. This can then be generalized to adding combinations of time derivatives to each equation since these vanish in the steady state. When using acceleration techniques such as residual smoothing, multigrid, etc. these are applied to the pressure rather than the density. Hence, the code duplicates the behavior of the incompressible equations for low speeds.

Fiterman, A.↗