Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variable energy generation”

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 1,153 records · Page 64
All sourcesNASA NTRSDOE OSTI
NASA NTRS2009

The Sodium Tail of the Moon

During the few days centered about new Moon, the lunar surface is optically hidden from Earth-based observers. However, the Moon still offers an observable: an extended sodium tail. The lunar sodium tail is the escaping "hot" component of a coma-like exosphere of sodium generated by photon-stimulated desorption, solar wind sputtering and meteoroid impact. Neutral sodium atoms escaping lunar gravity experience solar radiation pressure that drives them into the anti-solar direction forming a comet-like tail. During new Moon time, the geometry of the Sun, Moon and Earth is such that the anti-sunward sodium flux is perturbed by the terrestrial gravitational field resulting in its focusing into a dense core that extends beyond the Earth. An all-sky camera situated at the El Leoncito Observatory (CASLEO) in Argentina has been successfully imaging this tail through a sodium filter at each lunation since April 2006. This paper reports on the results of the brightness of the lunar sodium tail spanning 31 lunations between April 2006 and September 2008. Brightness variability trends are compared with both sporadic and shower meteor activity, solar wind proton energy flux and solar near ultra violet (NUV) patterns for possible correlations. Results suggest minimal variability in the brightness of the observed lunar sodium tail, generally uncorrelated with any single source, yet consistent with a multi-year period of minimal solar activity and non-intense meteoric fluxes.

Matta, M.↗
DOE OSTI2020

Probabilistic Forecast of All-Sky Solar Radiation Using Enhanced WRF-Solar: Preprint

This study presents enhancements of the Weather Research and Forecasting model with solar extensions (WRF-Solar) to provide probabilistic forecasts of solar radiation. Our approach builds ensemble WRF-Solar runs by introducing stochastic perturbations of variables that produce the largest uncertainties in predicting surface irradiance and clouds. The key variables are identified using tangent linear sensitivity analysis of six physics packages responsible for all-sky irradiance variability. An optimal strategy to stochastically perturb the selected variables is developed and applied to WRF-Solar to generate ensemble members for day-ahead solar prediction. The National Solar Radiation Database (NSRDB) is used to validate the ensemble forecast at arbitrary locations on the model grid. Preliminary results indicate that the proposed technique can potentially produce WRF-Solar ensembles providing reliable information of solar prediction uncertainty. This study describes the implemented methodology and initial results as well as future research to improve ensemble-based probabilistic forecasts with WRF-Solar.

all-sky↗
DOE OSTI2020

Accelerating Large-Format Metal Additive Manufacturing: How Controls R&D is Driving Speed, Scale, and Efficiency

This article highlights work at Oak Ridge National Laboratory’s Manufacturing Demonstration Facility to develop closed-loop, feedback control for laser-wire based Directed Energy Deposition, a form of metal Big Area Additive Manufacturing (m-BAAM), a process being developed in partnership with GKN Aerospace specifically for the production of Ti-6Al-4V pre-forms for aerospace components. A large-scale structural demonstrator component is presented as a case-study in which not just control, but the entire 3D printing workflow for m-BAAM is discussed in detail, including design principles for large-format metal AM, toolpath generation, parameter development, process control, and system operation, as well as post-print net-shape geometric analysis and finish machining. In terms of control, a multi-sensor approach has been utilized to measure both layer height and melt pool size, and multiple modes of closed-loop control have been developed to manipulate process parameters (laser power, print speed, deposition rate) to control these variables. Layer height control and melt pool size control have yielded excellent local (intralayer) and global (component-level) geometry control, and the impact of melt pool size control in particular on thermal gradients and material properties is the subject of continuing research. Further, these modes of control have allowed the process to advance to higher deposition rates (exceeding 7.5 lb/hr), larger parts (1-meter scale), shorter build times, and higher overall efficiency. The control modes are examined individually, highlighting their development, demonstration, and lessons learned, and it is shown how they operate concurrently to enable the printing of a large-scale, near net shape Ti-6Al-4V component.

Gibson, Brian↗
DOE OSTI2025

Environmental Conditions Affecting Global Mesoscale Convective System Occurrence

Abstract The ERA5 environments of mesoscale convective systems (MCSs), tracked from satellite observations, are assessed over a 20-yr period. The use of a large set of MCS tracks allows us to robustly test the sensitivity of the results to factors such as region, latitude, and diurnal cycle. We aim to provide novel information on environments of observed MCSs for assessments of global atmospheric models and to improve their ability to simulate MCSs. Statistical analysis of all tracked MCSs is performed in two complementary ways. First, we investigate the environments when an MCS has occurred at different spatial scales before and after MCS formation. Several environmental variables are found to show marked changes before MCS initiation, particularly over land. The vertically integrated moisture flux convergence shows a robust signal across different regions and when considering MCS initiation diurnal cycle. We also found spatial scale dependence of the environments between 200 and 500 km, providing new evidence of a natural length scale for use with MCS parameterization. In the second analysis, the likelihood of MCS occurrence for given environmental conditions is evaluated, by considering all environments and determining the probability of being in an MCS core or shield region. These are compared to analogous non-MCS environments, allowing discrimination between conditions suitable for MCS and non-MCS occurrence. Three environmental variables are found to be useful predictors of MCS occurrence: total column water vapor, midlevel relative humidity, and total column moisture flux convergence. Such relations could be used as trigger conditions for the parameterization of MCSs, thereby strengthening the dependence of the MCS scheme on the environment. Significance Statement Large storm systems called mesoscale convective systems form across Earth. These are collections of thunderstorms, with associated high-level clouds that produce substantial, lighter rainfall and modulate Earth’s energy balance. They produce hazardous weather conditions, such as floods and high winds, and are responsible for a high percentage of rainfall in many regions globally. We investigate the environmental conditions under which they form, so that we can understand the spatial extent of the environment which is important for their formation, and also where and when the effects of these storms might be felt. The novel information generated here should help improve the representation of these storms in weather and climate models, improving the prediction of rainfall, thunderclouds, and high-level clouds.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Fluorogenic aptamers resolve the flexibility of RNA junctions using orientation-dependent FRET

To further understand the transcriptome, new tools capable of measuring folding, interactions, and localization of RNA are needed. Although Förster resonance energy transfer (FRET) is an angle- and distance-dependent phenomenon, the majority of FRET measurements have been used to report distances, by assuming rotationally averaged donor–acceptor pairs. Angle-dependent FRET measurements have proven challenging for nucleic acids due to the difficulties in incorporating fluorophores rigidly into local substructures in a biocompatible manner. Fluorescence turn-on RNA aptamers are genetically encodable tags that appear to rigidly confine their cognate fluorophores, and thus have the potential to report angular-resolved FRET. In this work, we use the fluorescent aptamers Broccoli and Mango-III as donor and acceptor, respectively, to measure the angular dependence of FRET. Joining the two fluorescent aptamers by a helix of variable length allowed systematic rotation of the acceptor fluorophore relative to the donor. FRET oscillated in a sinusoidal manner as a function of helix length, consistent with simulated data generated from models of oriented fluorophores separated by an inflexible helix. Analysis of the orientation dependence of FRET allowed us to demonstrate structural rigidification of the NiCo riboswitch upon transition metal-ion binding. This application of fluorescence turn-on aptamers opens the way to improved structural interpretation of ensemble and single-molecule FRET measurements of RNA.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1996

Rayleigh Light Scattering for Concentration Measurements in Turbulent Flows

Despite intensive research over a number of years, an understanding of scalar mixing in turbulent flows remains elusive. An understanding is required because turbulent mixing has a pivotal role in a wide variety of natural and technologically important processes. As an example, the mixing and transport of pollutants in the atmosphere and in bodies of water are often dependent on turbulent mixing processes. Turbulent mixing is also central to turbulent combustion which underlies most hydrocarbon energy use in modern societies as well as in unwanted fire behavior. Development of models for combusting flows is therefore crucial, however, an understanding of scalar mixing is required before useful models of turbulent mixing and, ultimately, turbulent combustion can be developed. An important subset of turbulent flows is axisymmetric turbulent jets and plumes because they are relatively simple to generate, and because the provide an appropriate test bed for the development of general theories of turbulent mixing which can be applied to more complex geometries and flows. This paper focuses on a number of experimental techniques which have been developed at the National Institute of Standards and Development for measuring concentration in binary axisymmetric turbulent jets. In order to demonstrate the value of these diagnostics, some of the more important results from earlier and on-going investigations are summarized. Topics addressed include the similarity behavior of variable density axisymmetric jets, the behavior of absolutely unstable axisymmetric helium jets, and the role of large scale structures and scalar dissipation in these flows.

Pitts, William M.↗
DOE OSTI2021

UNH TDP - Concurrent Measurements of Inflow, Power Performance, and Loads for a Grid-Synchronized Vertical Axis Cross-Flow Turbine Operating in a Tidal Estuary

This data was collected between October 12 and December 15 of 2021 at the University of New Hampshire (UNH) and Atlantic Marine Energy Center (AMEC) turbine deployment platform (TDP). This data set includes over 29 days of grid connected turbine operation during this 65 day time frame. The priority for this measurement campaign was to collect data while the turbine was electrically connected to the grid by means of a rectifier and inverter. The Fall_2021_UNH_Measurement_Timeline.png highlights when each instrument was functioning and the Fall_2021_UNH_Test_Log.jpg indicates the four main regions for analysis available from this measurement campaign. The TDP is a floating structure moored on the Portsmouth facing side of Memorial Bridge pier #2, which spans the Piscataqua River between Portsmouth, NH and Kittery, ME. The Piscataqua River connects the Great Bay Estuary to the Gulf of Maine and the river currents are dominated by tidal forcing with water velocities exceeding 2.5 m/s during spring ebb tides at this site which were previously characterized by Kaelin Chancey (Assessment Of The Localized Flow And Tidal Energy Conversion System At An Estuarine Bridge - UNH MS Thesis 2019). The turbine under test was a modified New Energy Corporation (Calgary, CA) model EVG-025 4-blade H-Darrius type vertical axis cross flow turbine that rotates in the clockwise direction with a rotor diameter of 3.2m and blade length of 1.7m. The hydro-foil profile was a NACA 0021 with a 10 inch chord length and a blade preset pitch angle of +4deg with a positive angle corresponding with the toe in direction. The standard EVG-025 has a rotor diameter of 3.4m and its rated power output is 25kW at 3 m/s. The rotor diameter was reduced to accommodate the size of the existing TDP moon-pool. This project was pursued to quantify device performance for cross flow turbines operating in a marine environment. Accurate physical models, to characterize cross flow turbine performance, require real operational data sets due to the complexity of blade fluid interactions. This data can help support model development which will help predict turbine performance when analyzing perspective project locations in the future. Instrumentation was deployed to measure; water speed/direction, electrical power output, turbine shaft speed, turbine thrust force, and platform motion. Concurrent measurements of these parameters allow for correlations (cause and affect) to be inferred, allowing for characterization of device performance over a range of operating conditions. Water currents were measured using Acoustic Doppler Current Profilers (ADCP's) and Acoustic Doppler Velocimeters (ADV's) directly upstream and downstream of the turbine for inflow, wake and turbulence measurements. Electrical power output was measured using the Voltsys rectifier and the Shark power meter. Shaft speed was calculated based on the Voltsys measurements of the permanent magnet three phase generator AC generation frequency, coupled directly to the cross flow turbine under test (i.e., no gear box). Platform motions were captured using a Yost IMU (inertial measurement unit). Turbine thrust loading was measured using a reaction arm about the turbine deployment platform spanning beam, where two bi-directional load cells were connected to the system via a pinned connection. This submission includes zipped folders for each instrument containing quality controlled (QC'd) data in daily .csv files for the relevant duration specific to each instrument, along with separate .csv file that contains the units for each variable. Some instrument daily files are quite large and can pose a challenge for a visual spreadsheet editor to open. A processing software like MATLAB or Python is recommended. Note the degree of QC varied between each instrument due to time constraints. Particular time and attention was given to perform quality control tests on the acoustic based instruments that are particularly suscep...

16 TIDAL AND WAVE POWER↗
DOE OSTI2026

Structural Design of Bismuth Telluride Nanoplates through Process Variables

Binary pnictogen chalcogen compounds, primarily bismuth tellurides and selenides, are of great interest due to their applications in emerging quantum devices, as well as thermoelectric generators. The performance of bismuth telluride in these roles depends on its structure at the nanoscale, particularly the size, shape, and crystallinity of its nanocrystalline forms. However, current methods for controlling these features are often slow, inconsistent, or difficult to scale. Here, we demonstrate that through a solvothermal synthesis and hot injection process, precise control over the morphology of bismuth telluride nanoplates is possible with independent tuning of process variables, such as temperature and reaction time. We find that the nanoplate shape and internal porosity vary systematically with synthesis temperature and that the same morphological outcomes can be rapidly achieved at a fixed temperature by adjusting reaction duration. These results reveal that both the temperature and time can independently direct bismuth telluride morphological features, allowing for rapid, tunable synthesis strategies. Our approach offers a scalable framework, not only for bismuth telluride but also for related layered chalcogenides used in energy harvesting and quantum technologies.

Ackley, Jordan [Boise State Univ., ID (United Stat↗
DOE OSTI2020

Deep Reinforcement Learning Based Volt-VAR Optimization in Smart Distribution Systems

This paper develops a model-free volt-VAR optimization (VVO) algorithm via multi-agent deep reinforcement learning (DRL) in unbalanced distribution systems. This method is novel since we cast the VVO problem in distribution networks to an intelligent deep Q-network (DQN) framework, which avoids solving a specific optimization model directly when facing time-varying operating conditions in the systems. We consider statuses/ratios of switchable capacitors, voltage regulators, and smart inverters installed at distributed generators as the action variables of the agents. A delicately designed reward function guides these agents to interact with the distribution system, in the direction of reinforcing voltage regulation and power loss reduction simultaneously. The forward-backward sweep method for radial three-phase distribution systems provides accurate power flow results within a few iterations to the DRL environment. The proposed method realizes the dual goals for VVO. We test this algorithm on the unbalanced IEEE 13-bus and 123-bus systems. Numerical simulations validate the excellent performance of this method in voltage regulation and power loss reduction.

41 EE - Solar Energy Technologies Office (EE-4S)↗
DOE OSTI2024

MW-Class SOFC Pilot System Development

The overall objective of this project is to develop a low-cost, efficient, and reliable MWe-class Solid Oxide Fuel Cell (SOFC) power system towards commercial deployment in natural-gas fueled distributed generation applications. This report details the successful completion of all project objectives. The main goals of the project were to develop the conceptual design of a MWe-class SOFC power system, and to complete a techno-economic analysis (TEA) to determine system costs at low-volume production levels of >100MW per year. The combination of the BoP equipment cost and SOFC stack module factory cost provided the net fabrication cost of the selected 1MW SOFC system configuration. The results of the study showed that the cost target of less than $6000/kW is achievable for a FOAK 1MW SOFC system, exclusive of the first-time engineering cost. Also, the project included an exploratory analysis of the factory cost of next generation stack technology based on the Compact Solid-oxide Architecture (CSA) design at a high annual manufacturing rate of one-GW produced in a giga-factory. The BoP capital cost for one giga-watt power production was estimated using learning-curve analysis based on the historical manufacturing cost data related to similar high temperature Molten Carbonate fuel cell systems built by FuelCell Energy (FCE). The giga-factory manufactured SOFC cost estimate combined with the BoP cost at high production rates led to the estimation of entire system capital cost. Ultimately, a Cost of Electricity (COE) analysis was performed including natural gas fuel prices, installation and O&M costs, and stack replacement cost. The results of the techno-economic analysis were summarized by parametric estimation of COE for 1MW SOFC NOAK (Nth-of-a-Kind) system with natural gas prices as a variable parameter.

08 HYDROGEN↗
DOE OSTI2025

High-shade dryland agrivoltaic conditions enhanced carbon uptake and water-use efficiency in zucchini ( Cucurbita pepo )

Introduction: The increasing global demand for food and energy is intensifying land-use competition. Agrivoltaic systems are a multifunctional land-use approach that vertically integrates the production of agricultural crops and solar power on the same land area. Most food crops are adapted to full-sun conditions, and the physiological responses of these crops to the novel microclimate under solar panels remain poorly understood. We hypothesized that the microclimate beneath the high-density photovoltaic system would influence carbon uptake, water use, and yield outcomes of zucchini summer squash.Methods: We conducted a field experiment in a hot, semi-arid climate on zucchini (Cucurbita pepo). Plants were grown under an agrivoltaic system with a 75% ground cover ratio (GCR) and in a full-sun control plot, each with two irrigation regimes (100 and 50%). We measured leaf-level photosynthesis, microclimate variables, and fruit yield at plant maturity and throughout the growing season.Results: The agrivoltaic array reduced photosynthetically active radiation (PAR) by ~79%, resulting in a cooler (−1.1 °C), more humid environment with higher soil moisture. These microclimatic conditions enhanced midday photosynthesis and daily cumulative carbon uptake. However, fruit yield was consistently lower under the panels, indicating a shift in carbon allocation toward vegetative growth. Photosynthesis was primarily driven by PAR across treatments, while soil moisture significantly influenced photosynthesis only in the control plots, suggesting water limitation was alleviated under the panels.Discussion: These findings highlight a trade-off between improved physiological performance and reduced yield under high-density agrivoltaics. While the system buffered heat and drought stress and improved overall plant function, excessive shade reduced reproductive output. Optimizing panel density or selecting crops cultivated for non-fruit yields will be essential for balancing food production and energy generation in dryland agrivoltaic settings.

14 SOLAR ENERGY↗
DOE OSTI2025

Tailoring Nitinol for elastocaloric application

This study focuses on tailoring commercial Nitinol, the most commonly used elastocaloric material, for near-room-temperature cooling applications. Short heat treatments near 500 °C were used to fine-tune the material’s transition temperature, resulting in austenite finish temperatures ranging from 6.0 to 25.5°C and altered superelastic and elastocaloric properties. Plateau stresses decreased while temperature changes rose from 21.5°C up to 27.9°C at 6% strain. Significant variability in the Nitinol response when testing below its austenite finish temperature was observed. In conclusion, the effect of mechanical cycling on transition temperatures was also evaluated, demonstrating an increase for all the samples.

Efficiency↗
DOE OSTI2024

Development of an Optimal Variable-Pitch Controller for Floating Axial-Flow Marine Hydrokinetic Turbines

This article discusses the development of an optimal variable-pitch controller for floating, axial-flow marine turbines. Recently, OpenFAST, an open-source wind turbine modeling tool, has been extended to model marine turbines. A controller is necessary to simulate marine turbines for different load cases using OpenFAST, which greatly impacts the performance of the energy system. Previous studies have designed controllers using a linearized model of the marine turbine, which can be time-consuming and require the expertise of a control engineer. In this study, we use an automated approach that uses generic models of the marine turbine to identify the controller gains, which can expedite the process of designing a controller. Using an optimizer to identify the control system parameters can additionally improve the controller's performance. The optimal controller tuned using such an approach results in a 20% reduction in the tower-base damage equivalent loading and better tracking of the rated generator speed and power.

axial flow↗
NASA NTRS1991

Final design of a free-piston hydraulic advanced Stirling conversion system

Under the US Department of Energy's (DOEs) Solar Thermal Technology Program, Sandia National Laboratories is evaluating heat engines for solar distributed receiver systems. The final design is described of an engineering prototype advanced Stirling conversion system (ASCS) with a free-piston hydraulic engine output capable of delivering about 25 kW of electric power to a utility grid. The free-piston Stirling engine has the potential for a highly reliable engine with long life because it has only a few moving parts, has noncontacting bearings, and can be hermetically sealed. The ASCS is designed to deliver maximum power per year over a range of solar input with a design life of 30 years (60,000 h). The system includes a liquid Nak pool boiler heat transport system and a free-piston Stirling engine with high-pressure hydraulic output, coupled with a bent axis variable displacement hydraulic motor and a rotary induction generator.

Wallace, D. A.↗
NASA NTRS2013

Associating Long-term Gamma-ray Variability with the Superorbital Period of LS I + 61 Deg. 303

Gamma-ray binaries are stellar systems for which the spectral energy distribution (discounting the thermal stellar emission) peaks at high energies. Detected from radio to TeV gamma rays, the gamma-ray binary LS I + 61◦303 is highly variable across all frequencies. One aspect of this system's variability is the modulation of its emission with the timescale set by the approx. 26.4960 day orbital period. Here we show that, during the time of our observations, the gamma-ray emission of LS I + 61 deg. 303 also presents a sinusoidal variability consistent with the previously known superorbital period of 1667 days. This modulation is more prominently seen at orbital phases around apastron, whereas it does not introduce a visible change close to periastron. It is also found in the appearance and disappearance of variability at the orbital period in the power spectrum of the data. This behavior could be explained by a quasi-cyclical evolution of the equatorial outflow of the Be companion star, whose features influence the conditions for generating gamma rays. These findings open the possibility to use gamma-ray observations to study the outflows of massive stars in eccentric binary systems.

Be↗
NASA NTRS2015

The Two-Column Aerosol Project: Phase I - Overview and Impact of Elevated Aerosol Layers on Aerosol Optical Depth

The Two-Column Aerosol Project (TCAP), conducted from June 2012 through June 2013, was a unique study designed to provide a comprehensive data set that can be used to investigate a number of important climate science questions, including those related to aerosol mixing state and aerosol radiative forcing. The study was designed to sample the atmosphere be tween and within two atmospheric columns; one fixed near the coast of North America (over Cape Cod, MA) and a second moveable column over the Atlantic Ocean several hundred kilometers from the coast. The U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) Mobile Facility (AMF) was deployed at the base of the Cape Cod column, and the ARM Aerial Facility was utilized for the summer and winter intensive observation periods. One important finding from TCAP is that four of six nearly cloud-free flight days had aerosol layers aloft in both the Cape Cod and maritime columns that were detected using the nadir pointing second-generation NASA high-spectral resolution lidar (HSRL-2).These layer s contributed up to 60 of the total observed aerosol optical depth (AOD). Many of these layers were also intercepted by the aircraft configured for in situ sampling, and the aerosol in the layers was found to have increased amounts of biomass burning material and nitrate compared to aerosol found near the surface. In addition, while there was a great deal of spatial and day-to-day variability in the aerosol chemical composition and optical properties, no systematic differences between the two columns were observed.

aerosols↗
DOE OSTI2022

Integration of cryogenic energy storage with renewables and power plants: Optimal strategies and cost analysis

Energy storage is critical for overcoming challenges associated with the intermittency and the variable availability of renewable sources for decarbonizing the energy sector. Cryogenic energy storage (CES) is of interest due to its high technology readiness level, no geographical limitations, and moderate round-trip efficiency. The time-varying nature of demands and renewable availability needs to be considered at the design and integration stages of energy storage. We develop a mixed-integer nonlinear program (MINLP) model to obtain the energy storage costs on a daily basis for different scenarios that typically arise over an entire year. Using this optimization-based framework, we address key decision-making questions towards energy transition: What is the energy cost when CES is integrated with renewables and power plants? How does each scenario affect the overall energy cost? How much storage is needed for complete transition to renewables? What is the optimal integration towards 100% renewable energy? What are the optimal storage designs for both renewables and fossil-based power generation with current and future energy demands? Here, we discuss different scenarios and solutions to these questions.

25 ENERGY STORAGE↗
DOE OSTI2020

A simple and effective five-equation two-phase numerical model for liquid-vapor phase transition in cavitating flows

Numerical difficulties, notably the non-monotonic behavior of the Wood speed of sound and the volume fraction positivity, associated with the reduced five-equation two-phase flow model of Kapila et al. (2001) [A.K. Kapila, R. Menikoff, J.B. Bdzil, S.F. Son, D.S. Stewart, 2001. Two-phase modeling of deflagration-to-detonation transition in granular materials: reduced equations, Physics of Fluids 13(10), 3002–3024] have been resolved in the past through the introduction of a frozen speed of sound and an algebraic approach for mechanical relaxation afforded by a pressure non-equilibrium six-equation model proposed by [R. Saurel, F. Petitpas, R.A. Berry, 2009. Simple and efficient relaxation methods for interfaces separating compressible fluids, cavitating flows and shocks in multiphase mixture, J. Comput. Phys. 228, 1678–1712]. By contrast, it is explored and demonstrated in this work that these difficulties can in fact still be resolved within the numerical scheme for solving the reduced five-equation model by numerically replacing the Wood speed of sound for the estimates of wave speeds in the approximate Riemann solver HLLC with the monotonic mixture speed of sound for a transport five-equation model. For shock interface (artificial mixture separating pure or nearly pure fluids) interaction problems, with the apparent advantage of monotonic behavior of the speed of sound in the interface, the effect of the numerical replacement is also confined to the interface. Differences other than the behavior of the speed of sound within the interface in the solutions due to the replacement diminish with increasing resolution when reasonable solution can be obtained with Wood speed of sound. For cavitating/expansion problems in physical fluid mixture, it is pointed out and explained why acoustics in the numerical solutions still propagate at the Wood speed of sound (therefore consistent with the reduced five-equation model) even though in some cases a much higher speed of sound like the numerical replacement above for solving the reduced five-equation model or the frozen speed of sound for solving a six-equation model is used for the estimates of wave speeds in the HLLC scheme. A variant of the five-equation two-phase flow model by Saurel et al. (2008) [R. Saurel, F. Petitpas, R. Abgrall, 2008. Modelling phase transition in metastable liquids: application to cavitating and flashing flows, J. Fluid Mech. 607, 313–350] is then constructed for liquid-vapor phase transition in cavitating flows. The relaxation toward thermo-chemical equilibrium during phase transition is achieved by solving a simple system of algebraic equations for the equilibrium state variables for better efficiency, following Pelanti and Shyue (2014) [M. Pelanti, K.-M. Shyue, 2014. A mixture-energy-consistent six-equation two-phase numerical model for fluid with interfaces, cavitation and evaporation waves. J. Comput. Phys. 259, 331–357]. Therefore, the current model retains both the simplicity afforded by the five-equation model and the efficiency of the algebraic relaxation solver. An alternative algebraic approach for handling the non-conservative term (the so-called K∇ · u term) in the reduced five-equation model for mechanical equilibrium of a liquid-vapor mixture is also explored by enforcing the thermal equilibrium at the same time. Finally, numerical results of sample tests in both one and two dimensions in the literature as well as that in three dimensions demonstrate the effectiveness and ability of the proposed model to simulate cavitating flows. An interesting mechanism of shock generation by acoustics in water due to phase transition is then found by the numerical simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗