Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variable flow”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

A regional comparison of sub-daily flow variability in regulated and unregulated rivers in the United States

Regulating rivers for hydropower or other purposes can dramatically alter river flow patterns, including creating substantial changes in flow over short, minutes-to-hours-long timespans known as sub-daily flow variability (SDFV). The impacts of flexible hydropower production on flow and aquatic organisms are increasingly documented in research. However, the degree to which flow alteration relates to different hydropower operational modes in distinct geographical regions and seasons is not well understood. This study offers a methodology for regional- and species-appropriate evaluations of potential impacts of flow on fish based on sub-daily flow characteristics of hydropower operational modes. We analyzed 15-min discharge data between 2018 and 2021 from 69 USGS stream gages to compare SDFV in hydropeaking, run-of-river, and unregulated systems in the US Southeast and Pacific Northwest. Regulated systems exhibited significant SDFV downstream from hydropower facilities relative to unregulated systems, but specific impacts differed between regions. Regulated systems in the Southeast were characterized by high flow coefficients of variation and ratios (hydropeaking only) and extended durations of daily upramping flow phases. Regulated systems in the Pacific Northwest were characterized by many short flow phases per day and large portions of the day spent upramping. Pacific Northwest unregulated systems displayed the strongest seasonal flow patterns while Southeastern hydropeaking systems displayed the greatest SDFV. Given that SDFV impacts multiple dimensions of fish ecology, region-specific sub-daily flow signatures have important implications for understanding and mitigating potential community-, species-, and age-specific effects on fish in different parts of the country.

Fish↗
DOE OSTI2024

The environmental impact of hydropower: a systematic review of the ecological effects of sub-daily flow variability on riverine fish

Hydropower can help facilitate power grid decarbonization because it can respond to short-term changes in power demand and is comparatively more reliable than intermittent wind and solar. However, flexible hydropower operations can create rapid and abnormal fluctuations in downstream flow conditions, which can negatively impact aquatic ecosystems. Accordingly, we conducted a systematic review on the ecological effects of hydropower-driven sub-daily flow variability (SDFV) on riverine fishes. We reviewed and synthesized 109 articles relevant to fish-SDFV relationships from seven sources, most of which focused on Salmonids in North America and northern and western Europe and were published in the last 15 years. We found strong agreement in the literature that SDFV increases fish stranding risk, destabilizes habitat, and decreases production and diversity. We found moderate agreement that SDFV interrupts fish reproduction, increases or has no impact on condition, and prompts or discourages movement depending on local channel conditions. We found little to no agreement for relationships between SDFV and mortality, physiology, and behavior. The effects of SDFV on riverine fish ecology are intertwined in the complex suite of biotic and abiotic characteristics that structure aquatic ecosystems and are highly site-, species-, and life stage-specific. Assessments of the impact of SDFV on fish ecology should first characterize local habitat and channel quality and fish community composition to identify specific, measurable ecological outcomes to sustain or enhance, and then design mitigation strategies tailored to those ecological objectives.

13 HYDRO ENERGY↗
DOE OSTI2025

Investigating the relationship between simulation parameters and flow variables in simulating atmospheric gravity waves for wind energy applications

Wind farms, particularly offshore clusters, are becoming larger than ever before. Besides influencing the surface wind flow and the inflow for downstream wind farms, large wind farms can trigger atmospheric gravity waves in the inversion layer and the free atmosphere aloft. Wind-farm-induced gravity waves can cause adverse pressure gradients upstream of the wind farm, which contribute to the global blockage effect, and can induce favorable pressure gradients above and downstream of the wind farm that enhance wake recovery. Numerical modeling is a powerful means of studying these wind-farm-induced atmospheric gravity waves, but it comes with the challenge of handling spurious reflections of these waves from domain boundaries. Typically, approaches which employ radiation boundary conditions and forcing zones are used to avoid these reflections. However, the simulation setup of these approaches relies heavily on ad hoc processes. For instance, the widely used Rayleigh damping method requires ad hoc tuning to produce a setup that may only produce satisfactory results for a particular case. To provide more systematic guidance on setting up realistic simulations of atmospheric gravity waves, we conduct a large-eddy simulation (LES) study of flow over a 2D hill and through a wind farm canopy that explores the optimum domain size and damping layer setup depending on the fundamental parameters which determine the flow characteristics. In this work, we only consider linearly stratified conditions (i.e., no inversion layer), thereby focusing on internal gravity waves in the free atmosphere and their reflections from the domain boundaries. This type of flow is governed by a single Froude number, which dictates most of the internal wave properties, such as wavelength, amplitude, and direction. This, in turn, will dictate the optimum domain size and Rayleigh damping layer setup. We find the effective horizontal and vertical wavelengths (the representative wavelengths of the entire wave spectrum) to be the appropriate length scales to size the domain and damping layer thickness, and the optimal Rayleigh damping coefficient scales with the Brunt–Väisälä frequency. Considering Froude numbers seen in wind farm applications, we propose recommendations to limit the reflections to less than 10 % of the total upward-propagating wave energy. Typically, damping is done at the top boundary, but given the non-periodic lateral boundary conditions of practical wind farm simulation domains, we find that damping the inflow–outflow boundaries is of equal importance to damping the top boundary. The Brunt–Väisälä frequency-normalized damping coefficient should be between 1 and 10. The damping layer thickness should be at least one effective vertical wavelength; damping layers exceeding 1.5 times the vertical wavelength are found to be unnecessary. The domain length and height should accommodate at least one effective horizontal and vertical wavelength, respectively. Moreover, Rayleigh damping does not damp the waves completely, and the non-damped energy might accumulate over the simulation time.

17 WIND ENERGY↗
DOE OSTI2023

Emergence of Unstable Focused Flow Induced by Variable‐Density Flows in Vertical Fractures

Abstract Fluids with different densities often coexist in subsurface fractures and lead to variable‐density flows that control subsurface processes such as seawater intrusion, contaminant transport, and geologic carbon sequestration. In nature, fractures have dip angles relative to gravity, and density effects are maximized in vertical fractures. However, most studies on flow and transport through fractures are often limited to horizontal fractures. Here, we study the mixing and transport of variable‐density fluids in vertical fractures by combining three‐dimensional (3D) pore‐scale numerical simulations and visual laboratory experiments. Two miscible fluids with different densities are injected through two inlets at the bottom of a fracture and exit from an outlet at the top of the fracture. Laboratory experiments show the emergence of an unstable focused flow path, which we term a “runlet.” We successfully reproduce the unstable runlet using 3D numerical simulations and elucidate the underlying mechanisms triggering the runlet. Dimensionless number analysis shows that the runlet instability arises due to the Rayleigh‐Taylor instability (RTI), and flow topology analysis is applied to identify 3D vortices that are caused by the RTI. Even under laminar flow regimes, fluid inertia is shown to control the runlet instability by affecting the size and movement of vortices. Finally, we confirm the emergence of a runlet in rough‐walled fractures. Since a runlet dramatically affects fluid distribution, residence time, and mixing, the findings in this study have direct implications for the management of groundwater resources and subsurface applications.

Environmental Sciences & Ecology↗
DOE OSTI2025

Variable Refrigerant Flow Performance: Humidity and Efficiency in Multifamily Buildings in New York, NY

Variable refrigerant flow (VRF) systems have become common in large multifamily buildings. The upfront cost is appealing to developers, as a single system can provide all heating and cooling. VRF systems also use no on-site fossil fuels which can help avoid greenhouse gas emissions. However, there have been very few data on how VRF systems actually perform in buildings. Researchers monitored three new, multifamily buildings with a focus on VRF system sizing, humidity control, and energy efficiency.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

End-Use Savings Shapes Measure Documentation: Variable Refrigerant Flow With 25% Upsizing Allowance for Heating

This documentation focuses on a single end-use savings shape upgrade - a variable refrigerant flow with heat recovery (VRF HR) heating and cooling system coupled with a dedicated outdoor air system (DOAS) for ventilation and where 25% upsizing (or 125% of the original size) is allowed for heating in colder climates (i.e., heating dominant regions). This document will primarily discuss the additional changes to the sizing algorithm and modeling approach, while a comprehensive overview of the fundamental modeling methodology and background of the VRF HR DOAS upgrade, including applicability and other key assumptions can be found in the original documentation: Variable Refrigerant Flow with Heat Recovery and Dedicated Outdoor Air System. To provide more context on the 25% upsizing algorithm, if the building is a cooling dominant (e.g., design cooling load higher than design heating load ), then the outdoor unit capacity of the VRF heat pump is sized based on the design cooling load. However, if the building is heating dominant, then this building is a candidate for the 25% upsizing allowance. Once the 25% upsizing is allowed, if the 25% upsized capacity (or 125% from the original size) represented with design condition exceeds the design heating load, then the design heating load is used for calculating the rated capacity of the outdoor unit. And if the 25% upsized capacity represented with the design condition does not exceed the design heating load, then the 25% upsized capacity represented with rated condition is used for the capacity of the outdoor unit while the remaining heating load is handled with the supplemental/backup electric resistance coil.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2026

Multi-Split Variable Refrigerant Flow (VRF) System Building Energy Simulations Using Performance Maps

Multi-split variable refrigerant flow (VRF) systems are highly energy-efficient HVAC (heating, ventilation and air conditioning) technologies that connect a single outdoor unit to multiple independent indoor terminal units using a common refrigerant circuit and a variable-speed compressor. Building energy simulations that incorporate VRF systems help model their unique operational characteristics and predict energy consumption in specific building designs. Traditionally, EnergyPlus models these systems by employing multiple sets of performance curves to characterize both individual terminal units and the outdoor unit. However, producing these curves is labor intensive and error prone, and they often do not capture all the key input and output variables. This paper introduces a novel approach that uses multi-dimensional performance maps to model VRF systems in building environments for space cooling. In this approach, performance maps are developed at the component level—separately for the outdoor unit and for each indoor terminal. The new modeling method is validated within EnergyPlus via a Python plug-in that contains a simple solver loop to coordinate the component-level, indoor, and outdoor unit maps. Furthermore, because performance maps can span more variables than traditional performance curves, they offer the opportunity to implement advanced controls, such as enhanced dehumidification and compressor modulation. A VRF air conditioner’s hardware system was modeled using the DOE/ORNL Heat Pump Design Model, which was automated to produce extensive performance maps for both the indoor and outdoor units.

Shen, Bo [ORNL] (ORCID:0000000336600393)↗
DOE OSTI2024

Modeling and Validation of a Residential Multi-Functional Variable Refrigerant Flow Heat Pump System with Heat Recovery

To bridge the existing gap in modeling the variable refrigerant flow heat pump systems with heat recovery (VRFHR), we developed a suite of dynamic VRFHR system models in Modelica. These models are specifically tailored for residential multi-functional VRFHR (MF-VRFHR) applications, including space conditioning and domestic hot water (DHW) heating, utilizing both the TIL library for HVAC equipment and the Buildings library for thermal load calculations. The development comprises essential component models, including the newly developed heat recovery unit (HRU), along with system models that integrate the heat pump system and building envelope. These system models accommodate various operational modes such as heating-only, cooling-only, and heating-recovery (including heating-dominant and cooling-dominant) modes. Furthermore, we propose an efficient optimization-based model calibration method that identifies critical model parameters while utilizing a small amount of data obtained from either real systems or manufacturer's specifications. We demonstrate the effectiveness of these models and the proposed calibration method for a MF-VRFHR system installed in Richland, WA. The developed models are calibrated and validated using data collected under different operational modes during both heating and cooling seasons. The results show that the models capture the system dynamics and achieve high accuracy, with the coefficient of the variation of the root-mean-square-error less than 15% for variables such as outdoor unit power consumption, compressor speed, space temperature and DHW temperature. The validated models serve as a reliable representation of the MF-VRFHR system, facilitating the development and validation of optimized controls needed to realize the full benefits of integrated heat pump systems. Future research will utilize these models to develop advanced controls and optimize system performance for improved energy efficiency and demand flexibility.

Modeling, Variable refrigerant flow (VRF) systems,↗
DOE OSTI2024

A numerical method for simulating variable density flows in membrane desalination systems

Here, we present a novel method for simulating unsteady, variable density, fluid flows in membrane desalination systems. By assuming the density varies only with concentration and temperature, the scheme decouples the solution of the governing equations into two sequential blocks. The first solves the governing equations for the temperature and concentration fields, which are used to compute all thermophysical properties. The second block solves the conservation of mass and momentum equations for the velocity and pressure. We show that this is computationally more efficient than schemes that iterate over the full coupled equations in one block. We verify that the method achieves second-order spatial-temporal accuracy, and we use the method to investigate buoyancy-driven convection in a desalination process called vacuum membrane distillation. Specifically, we show that with gravity properly oriented, variations in temperature and concentration can trigger a double-diffusive instability that enhances mixing and improves water recovery. We also show that the instability can be strengthened by providing external heating.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

End-Use Savings Shapes Measure Documentation: Variable Refrigerant Flow with Heat Recovery and Dedicated Outdoor Air System

Building on the successfully completed effort to calibrate and validate the U.S. Department of Energy's ResStock™ and ComStock™ models over the past 3 years, the objective of this work is to produce national data sets that empower analysts working for federal, state, utility, city, and manufacturer stakeholders to answer a broad range of analysis questions. The goal of this work is to develop energy efficiency, electrification, and demand flexibility end-use load shapes (electricity, gas, propane, or fuel oil) that cover a majority of the high-impact, market-ready (or nearly market-ready) measures. "Measures" refers to energy efficiency variables that can be applied to buildings during modeling. An end-use savings shape is the difference in energy consumption between a baseline building and a building with an energy efficiency, electrification, or demand flexibility measure applied. It results in a time-series profile that is broken down by end use and fuel (electricity or on-site gas, propane, or fuel oil use) at each time step. ComStock is a highly granular, bottom-up model that uses multiple data sources, statistical sampling methods, and advanced building energy simulations to estimate the annual subhourly energy consumption of the commercial building stock across the United States. The baseline model intends to represent the U.S. commercial building stock as it existed in 2018. The methodology and results of the baseline model are discussed in the final technical report of the End-Use Load Profiles project. This documentation focuses on a single heating, ventilation, and air-conditioning (HVAC) end-use savings shape measure - a variable refrigerant flow with heat recovery (VRF HR) heating and cooling system coupled with a dedicated outdoor air system (DOAS) for ventilation. This measure replaces existing multi-zone variable air volume (VAV) systems or single-zone rooftop units (RTU) with a VRF HR system coupled with a DOAS that includes an energy/heat recovery ventilator (E/HRV). The measure covers 53% of exisiting building stock's floor area and is not applicable to HVAC system types using district heating or cooling or buildings/spaces that include high-ventilation spaces such as kitchens where the amount of exhaust air is large. A DOAS with E/HRV is used to provide required outdoor ventilation air to spaces since ventilation air is generally not supplied by a VRF HR system. An exhaust air energy recovery ventilator (ERV ) with sensible and latent heat exchange is added to humid climate zones while a heat recovery ventilator (HRV ) with sensible only exchange is added to drier climate zones. The ERV is modeled as a fixed membrane plate counterflow heat exchanger, while the HRV is modeled as a sensible-only fixed aluminum plate counterflow heat exchanger. Both systems include a bypass (for temperature control and economizer lockout) and minimum exhaust temperature control for frost prevention.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

Influence of the Kuroshio Intrusion on Deep Flow Intraseasonal Variability in the Northern South China Sea

Abstract Interactions between the open ocean and marginal seas have been suggested to be critical to the redistribution and dissipation of global energy. Here, we propose a mechanism for the upper open ocean influencing the deep flow in marginal seas that hinges on the formation and propagation of topographic Rossby waves (TRWs). Observations and high‐resolution simulations suggest substantial intraseasonal variability with periods of 5–60 days associated with the deep flow over continental slopes of the northern South China Sea (NSCS). These fluctuations generally account for over 40% of the total deep kinetic energy variability and the number can reach 70% over the slopes to the west of the Luzon Strait, southwest of the Dongsha Islands, and northeast of the Xisha Islands. By examining the spatiotemporal features of the fluctuations, we demonstrate that the intraseasonal variability of the deep flow in the NSCS is closely associated with TRWs. Using a recently developed multiscale energetics analysis in combination with a ray tracing model, we show that the energy sources of TRWs can be traced back to the east of the Dongsha Islands, where the Kuroshio intrusion and related eddies energize the TRWs primarily through pressure work. These waves propagate westward across the NSCS and drive the intraseasonal variability of the deep flow over continental slopes. Our findings highlight an energy pathway from the open ocean western boundary current to the abyssal marginal sea that could modulate regional circulation as well as exchanges between major ocean basins.

Quan, Qi↗
DOE OSTI2020

Overland flow numerical model prediction, Lower Triangle Region in East River Watershed, Colorado, 3 days

This data package contains numerical simulation results of surface flow variables such as flow velocity and water depth in Lower Triangle Region in East River Watershed, Colorado. The surface flow is a consequence of a high intensity rainfall event with a total duration of 3 days, available at a resolution of 10 minutes. The results are computed on triangular multiresolution meshes with resolutions ranging from 10 meter to 80 meter. The data package also contains a simulation on a uniform triangular mesh with a resolution of 10 meter. The simulations consider surface flow only and neglect subsurface flow, infiltration, and evapotranspiration. The purpose of the data is to assess the quality of a mesh refinement strategy.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Kinetics of solid–liquid interface motion in molecular dynamics and phase-field models: crystallization of chromium and silicon

The results of molecular dynamics (MD) simulations of the crystallization process in one-component materials and solid solution alloys reveal a complex temperature dependence of the velocity of the crystal–liquid interface featuring an increase up to a maximum at 10–30% undercooling below the equilibrium melting temperature followed by a gradual decrease of the velocity at deeper levels of undercooling. At the qualitative level, such non-monotonous behaviour of the crystallization front velocity is consistent with the diffusion-controlled crystallization process described by the Wilson–Frenkel model, where the almost linear increase of the interface velocity in the vicinity of melting temperature is defined by the growth of the thermodynamic driving force for the phase transformation, while the decrease in atomic mobility with further increase of the undercooling drives the velocity through the maximum and into a gradual decrease at lower temperatures. At the quantitative level, however, the diffusional model fails to describe the results of MD simulations in the whole range of temperatures with a single set of parameters for some of the model materials. The limited ability of the existing theoretical models to adequately describe the MD results is illustrated in the present work for two materials, chromium and silicon. It is also demonstrated that the MD results can be well described by the solution following from the hodograph equation, previously found from the kinetic phase-field model (kinetic PFM) in the sharp interface limit. The ability of the hodograph equation to describe the predictions of MD simulation in the whole range of temperatures is related to the introduction of slow (phase field) and fast (gradient flow) variables into the original kinetic PFM from which the hodograph equation is obtained. The slow phase-field variable is responsible for the description of data at small undercoolings and the fast gradient flow variable accounts for local non-equilibrium effects at high undercoolings. The introduction of these two types of variables makes the solution of the hodograph equation sufficiently flexible for a reliable description of all nonlinearities of the kinetic curves predicted in MD simulations of Cr and Si. This article is part of the theme issue ‘Transport phenomena in complex systems (part 1)’.

Science & Technology - Other Topics↗
DOE OSTI2020

Identifying Characteristics of Wildfire Towers and Troughs

Wildfire behavior is dictated by the complex interaction of numerous physical phenomena including dynamic ambient and fire-induced winds, heat transfer, aerodynamic drag on the wind by the fuel and combustion. These phenomena create complex feedback effects between the fire and its surroundings. In this study, we aim to study the mechanisms by which buoyant flame dynamics along with vortical motions and instabilities control wildfire propagation. Specifically, this study employs a suite of simulations conducted with the physics-based coupled fire-atmosphere behavior model (FIRETEC). The simulations are initialized with a fire line and the fires are allowed to propagate on a grass bed, where the fuel heights and wind conditions are varied systematically. Flow variables are extracted to identify the characteristics of the alternating counter-rotational vortices, called towers and troughs, that drive convective heat transfer and fire spread. These vortices have previously been observed in wildfires and laboratory fires, and have also been observed to arise spontaneously in FIRETEC due to the fundamental physics incorporated in the model. However, these past observations have been qualitative in nature and no quantitative studies can be found in the literature which connected these coherent structures fundamental to fire behavior with the constitutive flow variables. To that end, a variety of state variables are examined in the context of these coherent structures under various wind profile and grass height conditions. Identification of various correlated signatures and fire-atmosphere feedbacks in simulations provides a hypothesis that can be tested in future observational or experimental efforts, potentially assisting experimental design, and can aid in the interpretation of data from in situ detectors.

54 ENVIRONMENTAL SCIENCES↗