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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Data-driven closure modeling for hypersonic turbulent flows

The Reynolds-averaged Navier–Stokes (RANS) equations remain a workhorse technology for simulating compressible fluid flows of practical interest. Due to model-form errors, however, RANS models can yield erroneous predictions that preclude their use on mission-critical problems. This report summarizes work performed from FY22-FY24 focused on improving RANS models for hypersonic flows using data-driven modeling and scientific machine learning. In this work we: 1. Investigate the current capabilities of RANS models in Sandia’s parallel aerodynamics and re-entry code (SPARC) for hypersonic flows with a focus on shock boundary layer interactions (SBLIs), 2. Assess several established corrections that exist in the literature aimed at improving predictions for SBLIs, 3. Develop improved models for the Reynolds stress tensor using tensor-basis neural networks, 4. Develop a neural-network-based variable turbulent Prandtl number model to reduce errors in wall heating in SBLIs. 5. Begin future investigations including employing the LIFE framework to improve wall heating predictions in SBLIs as well as the ensemble Kalman filter. We find that current RANS models in SPARC are deficient for complex SBLI flows. In particular, no current model jointly predicts wall heat flux, wall shear stress, and wall pressure with reasonable accuracy. Existing corrections help, but do not alleviate this issue altogether. The development of improved models for the Reynolds stress tensor via tensor-basis neural networks results in more predictive RANS models across a suite of low-speed and high-speed cases. For hypersonic boundary layers, the inclusion of the wall-normal Reynolds stress via TBNNs has an appreciable impact on the wall-normal momentum balance and wall quantities. However, we find that improvements to the Reynolds stress tensor do not address the over-prediction in wall heat flux in SBLIs. We find that a neural-network-based variable turbulent Prandtl number model systematically and substantially improves wall heating predictions for a range of SBLI cases.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Alternating Direction Decomposition with Strong Bounding and Convexification (ADDSBC) for Solving Security Constrained AC Unit Commitment Problems

This project aims to develop efficient and robust computational methods for solving the security-constrained unit commitment and alternating current optimal power flow problem (SC-UC-ACOPF). The SC-UC-ACOPF problem is at the center of the short-term operation of the U.S. Power Grid. It is solved every week, every day, and every 10 minutes to plan for the optimal action of electricity generation and consumption by minimizing the generation cost and maintaining power system reliability against potential disruptions of equipment failures. In mathematical terms, SC-UC-ACOPF is a challenging large-scale mixed-integer nonlinear optimization model. This means that the decisions involve both discrete variables, e.g. the turning on and off of generators and switching of transmission lines and transformers, and continuous decisions, e.g. the amount of energy generated by each generator and the power flows in the power grid. The physics of the power flow is described by nonlinear equations involving real and reactive power and bus voltages. Another key feature is the large number of contingencies, i.e. the system needs to stay reliable in face of failure of any one equipment, such as transmission lines and generators. The U.S. power grids are extremely complicated and large scale with more than 5,000 generators, 50,000 buses, and 100,000 high-voltage transmission lines, making the SC-UC-ACOPF a very large-scale computation challenge. The research developed in this project aims to solve the SC-UC-ACOPF problems in the three timescales, i.e. weekly, daily, and every 10-min. The proposed computational methods are built on a principled algorithmic approach of decomposition and penalization. More specifically, the algorithm develops spatial and temporal decomposition by exploiting the strong temporal coupling and weak spatial coupling of the UC problem and the complementary feature, i.e. weak temporal coupling and strong spatial coupling of the ACOPF problem. The algorithm also leverages recent progresses in strong convex relaxation of ACOPF. A unique feature of the proposed approach is that it generates a valid, global upper bound on the optimal maximum profit. In this way, a global optimality gap is available to measure the quality of the solution. To further speed up computation, the research team has developed a plethora of effective heuristics to strengthen the iterative penalty-based decomposition framework. For instance, a heuristic is developed to construct inner approximations of the time coupling constraints within the time decoupled problems. Contingencies are pre-screened and low-rank matrix computation is exploited to find the almost unique solution to each contingency. A novel heuristic for line switching is proposed and tested with positive impacts on instances where line switching is beneficial. Taking a systematic approach and carefully handling every detail of the problem pays off. The TIM-GO’s performance throughout the trials and the final event was stellar. TIM-GO garnered the second highest total prize money and is ranked in the top three positions across all categories of comparison.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Deep Learning for Subsurface Flow: A Comparative Study of U‐Net, Fourier Neural Operators, and Transformers in Underground Hydrogen Storage

Subsurface flow research is essential for the sustainable management of natural resources and the environment. Deep learning (DL) has significantly advanced this field by developing efficient and accurate surrogate models to replace computationally expensive physics‐based simulations. These surrogate models are commonly used to predict the spatiotemporal evolution of state variables, such as gas saturation and reservoir pressure, in heterogeneous geological formations. Despite the various DL models applied to this task, there is a lack of studies systematically comparing their performance. This absence of comparative analysis leads to somewhat arbitrary DL model selection in subsurface flow research, resulting in suboptimal performance and potentially inaccurate predictions. To bridge this gap, we conduct a systematic comparison study of three popular DL architectures—U‐Net, Fourier Neural Operators (FNO), and Segmentation Transformer (SETR)—in surrogate modeling of underground hydrogen storage (UHS). We focus on UHS due to its promise of enhancing clean energy resilience and its cyclic operational conditions that represent common scenarios in various subsurface applications. We evaluate the models based on accuracy, training cost, and inference speed. The comparison shows that U‐Net achieves the highest accuracy, followed by SETR and FNO. Despite its lower accuracy, FNO has the highest inference speed. SETR offers competitive accuracy with the least training memory usage, demonstrating the potential of transformers in learning subsurface flow. Our results provide guidance for selecting DL models for surrogate modeling in a wide range of subsurface flow problems.

42 ENGINEERING↗
DOE OSTI2022

Warming Pattern over the Northern Hemisphere Midlatitudes in Boreal Summer 1979–2020

Abstract Significant surface air temperature warming during summer 1979–2020 is not uniformly distributed in the northern midlatitudes over land but rather is confined to several longitudinal sectors including Europe, central Siberia and Mongolia, and both coasts of North America. These hot spots are accompanied by a chain of high pressure ridges from an anomalous, circumglobal Rossby wave train in the upper troposphere. From reanalysis data and several baseline experiments from phase 6 of the Coupled Model Intercomparison Project (CMIP6), we find that the circulation trend pattern is associated with fluctuations of the Atlantic multidecadal variability (AMV) and the interdecadal Pacific oscillation. The phase shift of AMV in the 1990s is particularly noteworthy for accelerating warming averaged over the northern midlatitude land. The amplitude of the observed trend in both surface air temperature and the upper-level geopotential height generally falls beyond the range of multidecadal trends simulated by the CMIP6 preindustrial control runs, supporting the likelihood that anthropogenic forcing played a critical role in the observed trend. On the other hand, the fidelity of the simulated low-frequency modes of variability and their teleconnections, especially on multidecadal time scales, is difficult to assess because of the relatively short observational records. Our mechanistic modeling results indicate that synoptic eddy–mean flow interaction is a key to the formation of the anomalous wave train but how the multidecadal modes can modulate the synoptic eddies through atmosphere–ocean and atmosphere–land interactions remains poorly understood. This gap in our knowledge makes it challenging to quantify the roles of the low-frequency modes and external forcings in causing the observed multidecadal trends.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Enhancing ACPF Analysis: Integrating Newton-Raphson Method with Gradient Descent and Computational Graphs

This paper presents a new method for enhancing Alternating Current Power Flow (ACPF) analysis. The method integrates the Newton-Raphson (NR) method with Enhanced-Gradient Descent (GD) and computational graphs. The integration of renewable energy sources in power systems introduces variability and unpredictability, and this method addresses these challenges. It leverages the robustness of NR for accurate approximations and the flexibility of GD for handling variable conditions, all without requiring Jacobian matrix inversion. Furthermore, computational graphs provide a structured and visual framework that simplifies and systematizes the application of these methods. The goal of this fusion is to overcome the limitations of traditional ACPF methods and improve the resilience, adaptability, and efficiency of modern power grid analyses. We validate the effectiveness of our advanced algorithm through comprehensive testing on established IEEE benchmark systems. Furthermore, our findings demonstrate that our approach not only speeds up the convergence process but also ensures consistent performance across diverse system states, representing a significant advancement in power flow computation.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Computationally efficient CFD prediction of bubbly flow using physics-guided deep learning

To realize efficient computational fluid dynamics (CFD) prediction of two-phase flow, a multi-scale framework was proposed in this paper by applying a physics-guided data-driven approach. Instrumental to this framework, Feature Similarity Measurement (FSM) technique was developed for error estimation in two-phase flow simulation using coarse-mesh CFD, to achieve a comparable accuracy as fine-mesh simulations with fast-running feature. In this work, by defining physics-guided parameters and variable gradients as physical features, FSM has the capability to capture the underlying local patterns in the coarse-mesh CFD simulation. Massive low-fidelity data and respective high-fidelity data are used to explore the underlying information relevant to the main simulation errors and the effects of phenomenological scaling. By learning from previous simulation data, a surrogate model using deep feedforward neural network (DFNN) can be developed and trained to estimate the simulation error of coarse-mesh CFD. In a demonstration case of two-phase bubbly flow, the DFNN model well captured and corrected the unphysical “peaks” in the velocity and void fraction profiles near the wall in the coarse-mesh configuration, even for extrapolative predictions. The research documented supports the feasibility of the physics-guided deep learning methods for coarse mesh CFD simulations which has a potential for the efficient industrial design.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Design and execution of a verification, validation, and uncertainty quantification plan for a numerical model of left ventricular flow after LVAD implantation

Left ventricular assist devices (LVADs) are implantable pumps that act as a life support therapy for patients with severe heart failure. Despite improving the survival rate, LVAD therapy can carry major complications. Particularly, the flow distortion introduced by the LVAD in the left ventricle (LV) may induce thrombus formation. While previous works have used numerical models to study the impact of multiple variables in the intra-LV stagnation regions, a comprehensive validation analysis has never been executed. The main goal of this work is to present a model of the LV-LVAD system and to design and follow a verification, validation and uncertainty quantification (VVUQ) plan based on the ASME V&V40 and V&V20 standards to ensure credible predictions. The experiment used to validate the simulation is the SDSU cardiac simulator, a bench mock-up of the cardiovascular system that allows mimicking multiple operation conditions for the heart-LVAD system. The numerical model is based on Alya, the BSC’s in-house platform for numerical modelling. Alya solves the Navier-Stokes equation with an Arbitrary Lagrangian-Eulerian (ALE) formulation in a deformable ventricle and includes pressure-driven valves, a 0D Windkessel model for the arterial output and a LVAD boundary condition modeled through a dynamic pressure-flow performance curve. The designed VVUQ plan involves: (a) a risk analysis and the associated credibility goals; (b) a verification stage to ensure correctness in the numerical solution procedure; (c) a sensitivity analysis to quantify the impact of the inputs on the four quantities of interest (QoIs) (average aortic root flow $Q^{avg}_{Ao}$, maximum aortic root flow $Q^{max}_{Ao}$, average LVAD flow $Q^{avg}_{VAD}$, and maximum LVAD flow $Q^{max}_{VAD}$; (d) an uncertainty quantification using six validation experiments that include extreme operating conditions. Numerical code verification tests ensured correctness of the solution procedure and numerical calculation verification showed a grid convergence index (GCI)95% <3.3%. The total Sobol indices obtained during the sensitivity analysis demonstrated that the ejection fraction, the heart rate, and the pump performance curve coefficients are the most impactful inputs for the analysed QoIs. The Minkowski norm is used as validation metric for the uncertainty quantification. It shows that the midpoint cases have more accurate results when compared to the extreme cases. The total computational cost of the simulations was above 100 [core-years] executed in around three weeks time span in Marenostrum IV supercomputer. This work details a novel numerical model for the LV-LVAD system, that is supported by the design and execution of a VVUQ plan created following recognised international standards. We present a methodology demonstrating that stringent VVUQ according to ASME standards is feasible but computationally expensive.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

A Deep Learning-Based Workflow for Fast Prediction of 3D State Variables in Geological Carbon Storage: A Dimension Reduction Approach

In this study, we used deep learning techniques, which are a form of artificial intelligence, to create fast and effective models for predicting how fluids flow in underground geological formations. This is important for managing geological carbon storage, a method used to fight climate change by storing carbon dioxide underground. The challenge lies in the complex nature of these underground spaces and the large amount of data needed to accurately simulate them. To overcome these issues, we developed a new workflow that reduces the data’s complexity before training the deep learning model and then reconstructs the predicted results in their original form. We also proposed a unique approach to handle the specific complexities found in 3D saturation fields, a crucial aspect of fluid flow prediction. We tested our method using real-world data from the Gulf of Mexico. Our results show that our approach not only accurately predicts fluid behavior but also significantly reduces computation time. This will greatly improve real-time decision-making and risk assessment in large-scale geological carbon storage operations.

Wang, Hongsheng↗
DOE OSTI2025

Spatial Replication Is Important for Developing Landscape Genetic Inferences for a Wetland Salamander

Habitat fragmentation is a pressing threat to wildlife populations, and maintenance of gene flow between populations is an essential goal of conservation. Resistance surfaces have emerged as an important tool for modelling connectivity and developing management strategies to mitigate effects of habitat fragmentation. However, recent studies have noted inconsistencies in the factors most strongly associated with connectivity across different landscapes. Thus, replication of genetic-based resistance surface optimisation across landscapes may be necessary for making robust conclusions about the influence of environmental variables. Accordingly, replication represents a substantive challenge and opportunity in the field of landscape genetics. In this study, we conducted replicated landscape genetic analyses across five landscapes in Tennessee and Kentucky for a threatened wetland amphibian, the four-toed salamander (Hemidactylium scutatum). We tested multiple hypotheses of how different landscape features that could directly affect small, desiccation-intolerant amphibians (e.g., canopy cover) influenced gene flow and assessed the appropriate scale at which to model different features. We found some concordance in the landscape features that influenced gene flow (e.g., a common importance of forest cover and topography), but also some differences—potentially owing to the difference in variability of predictors across landscapes. We also found discordance in the scale of effect of different features across landscapes. In conclusion, our work emphasises that flat areas of moist forest not bisected by roads may be important for H. scutatum conservation, and our replicated design allows us to identify relationships that would have been missed if only using one study site.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Radial Wall Jet Flow Over Sigmoidal Surfaces

Radial wall jet flows across flat smooth surfaces have been studied for decades. These studies show that the radial velocity of these jets decays inversely with distance from the nozzle with modest contribution from friction (Poreh, et al., 1967; Rajaratnam, 1976). However, the extent to which flat surface results apply to curved surfaces remains unclear. In this article we explore the influence of settled particle bed slope on these velocity profiles. We model the step change in thickness as a sigmoidal curve of variable steepness and use conservation of momentum to evaluate the velocity profile. We show that the velocity profile attenuates because of curvature. Next, we evaluate the influence of surface curvature on radial wall jet flow. Jet flows over particle beds often introduce curvature in the particle bed profile but the influence of the developed curvature on the velocity profile has not been explored. Here we develop a solution for steady fixed beds based on conservation of momentum. We find that surface curvature has a significant influence on the velocity decay coefficients and celerity of radial wall jets, provided there is velocity slip in the vicinity of the particle bed interface, which is strictly true for particle surfaces. Conservation of momentum predicts conditions where the forward momentum of the flow is directed completely upward. The solution identifies two new dimensionless groups that determine whether a curved surface is sufficient to block radial flow and force flow vertically.

Surface curvature, jet flow, multiphase flow↗
DOE OSTI2020

A physics-informed operator regression framework for extracting data-driven continuum models

The application of deep learning toward discovery of data-driven models requires careful application of inductive biases to obtain a description of physics which is both accurate and robust. We present here a framework for discovering continuum models from high fidelity molecular simulation data. Our approach applies a neural network parameterization of governing physics in modal space, allowing a characterization of differential operators while providing structure which may be used to impose biases related to symmetry, isotropy, and conservation form. Here, we demonstrate the effectiveness of our framework for a variety of physics, including local and nonlocal diffusion processes and single and multiphase flows. For the flow physics we demonstrate this approach leads to a learned operator that generalizes to system characteristics not included in the training sets, such as variable particle sizes, densities, and concentration.

42 ENGINEERING↗
DOE OSTI2025

Intrinsic Kinetics of Polyethylene Terephthalate Pyrolysis via Micropyrolysis and Multivariate Chromatographic Analysis

This study provides an in-depth investigation of the primary decomposition of polyethylene terephthalate (PET) via pyrolysis, employing an experimental-analytic workflow that integrates design of experiments (DoE), micropyrolysis coupled with comprehensive two-dimensional gas chromatography (GC×GC), and multivariate data analysis to verify intrinsic kinetic conditions and elucidate evolving product distributions for mapping key reaction pathways. Peaks that could not be identified using commercial spectral libraries were assigned using Mass Frontier simulations, enabling the identification of divinyl terephthalate, ethyl vinyl terephthalate, and 2-(benzoyloxy)ethyl vinyl terephthalate. A polar×polar (non-orthogonal) column set tailored for the detection of carboxylic acids enhanced the quantification of benzoic acid, 4-vinylbenzoic acid, 4-ethylbenzoic acid, and methylbenzoic acid by up to 6-fold relative to an orthogonal column combination (non-polar×mid-polar). Moreover, pyrolysis variables were systematically evaluated using a Box- Behnken design (BBD), encompassing pyrolysis temperature (500−600 °C), sample weight (50−150 μg), and carrier gas flow rate (100−300 mL min −1 ). Among these, pyrolysis temperature was the only statistically significant factor influencing product yields, ranging from 58.78 to 84.26 wt %. In contrast, neither the sample weight nor the carrier gas flow rate had a significant effect on product yields within the evaluated experimental space. At 600 °C, the major pyrolysis products were benzoic acid (up to 20.20 ± 1.46 wt %) and CO 2 (up to 21.28 ± 1.46 wt %), which can be produced through decarboxylation reactions. These findings underscore the critical importance of selecting appropriate analytical columns for the accurate quantification of heteroatomcontaining products such as carboxylic acids, which may otherwise be underestimated or undetected due to their reactivity with the stationary phase of non-polar and mid-polar columns, as well as other GC components. They also highlight the importance of selecting pyrolysis conditions for investigating the primary decomposition of PET under an isothermal kinetically limited regime.

aromatic compounds↗
DOE OSTI2023

A Gradient Based Subgrid-Scale Parameterization for Ocean Mesoscale Eddies

Mesoscale eddies play an important role in transport of heat and biogeochemical tracers in the global ocean circulation. Resolving these energetic eddies, however, is challenging in ocean general circulation models (OGCM) because it requires a horizontal grid spacing of ≲1/8° that is computationally expensive. As a result, we are required to parameterize mesoscale eddy effects on large-scale ocean flows. In this work, we introduce a new subgrid-scale (SGS) model that is developed based on a Taylor series expansion of resolved variables to parameterize subgrid mesoscale eddy transports and momentum fluxes in OGCM. We have performed an a priori study to evaluate the performance of our new gradient model using high-resolution ocean simulations. Our results show that the gradient model well predicts the actual SGS thickness fluxes in the zonal and meridional directions in coarse-resolution simulations with the grid spacing ≳1/4°. The unresolved kinetic energy at the ocean surface is also skillfully estimated. More importantly, unlike current mesoscale eddy parameterizations, which are mainly developed based on an assumption of flat bottom ocean, our new SGS model can capture the structure of unresolved standing meanders at the ocean surface. We have also developed a dynamic procedure for setting in non-dimensional parameters in our new parameterization through a non-ad hoc and tuning-free method. Overall, this work suggests that implementing the gradient model in OGCM can improve the model accuracy with an affordable computational cost in eddy-permitting and non-eddying simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Tracking blobs in the turbulent edge plasma of a tokamak fusion device

Abstract The analysis of turbulence in plasmas is fundamental in fusion research. Despite extensive progress in theoretical modeling in the past 15 years, we still lack a complete and consistent understanding of turbulence in magnetic confinement devices, such as tokamaks. Experimental studies are challenging due to the diverse processes that drive the high-speed dynamics of turbulent phenomena. This work presents a novel application of motion tracking to identify and track turbulent filaments in fusion plasmas, called blobs, in a high-frequency video obtained from Gas Puff Imaging diagnostics. We compare four baseline methods (RAFT, Mask R-CNN, GMA, and Flow Walk) trained on synthetic data and then test on synthetic and real-world data obtained from plasmas in the Tokamak à Configuration Variable (TCV). The blob regime identified from an analysis of blob trajectories agrees with state-of-the-art conditional averaging methods for each of the baseline methods employed, giving confidence in the accuracy of these techniques. By making a dataset and benchmark publicly available, we aim to lower the entry barrier to tokamak plasma research, thereby greatly broadening the community of scientists and engineers who might apply their talents to this endeavor.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

A High-Fidelity Modeling Tool to Support the Design of Oxy-Combustors for Direct-Fired Supercritical CO2 Cycles

The challenge in the design of oxy-combustors for direct-fired supercritical CO 2 (sCO 2 ) cycles is in addressing disparate performance metrics and objectives. Key design parameters to consider include, among others, injector design for mixing and flame stability, split of recycled CO 2 diluent between injectors and cooling films, target flame temperatures to control noncondensable products, and strategies to inject the diluent CO 2 for film cooling and thermal control. In order to support novel oxy-combustor designs, a high-fidelity yet numerically efficient modeling framework based on the CRUNCH CFD® flow solver is presented, featuring key physics-based submodels relevant in this regime. For computational efficiency in modeling large kinetic sets, a flamelet/progress variable (FPV) based tabulated-chemistry approach is utilized featuring a three-stream extension to allow for the simulation of the CO 2 film cooling stream in addition to the fuel and oxidizer streams. Finite rate chemistry effects are modeled in terms of multiple progress variables for the primary flame as well as for slower-evolving chemical species such as NO x and SO x contaminants. Real fluid effects are modeled using advanced equations of states. The predictive capabilities of this computationally tractable design support tool are demonstrated on a conceptual injector design for an oxy-combustor operating near 30 MPa. Simulations results provide quantitative feedback on the effectiveness of the film cooling as well as the level of contaminants (CO, NO, and N) in the exhaust due to impurities entering from the injectors. These results indicate that this framework would be a useful tool for refining and optimizing the oxy-combustor designs as well as risk mitigation analyses.

Engineering↗
DOE OSTI2024

The Four-Dimensional Carbon Cycle of the Southern Ocean

The Southern Ocean plays a fundamental role in the global carbon cycle, dominating the oceanic uptake of heat and carbon added by anthropogenic activities and modulating atmospheric carbon concentrations in past, present, and future climates. However, the remote and extreme conditions found there make the Southern Ocean perpetually one of the most difficult places on the planet to observe and to model, resulting in significant and persistent uncertainties in our knowledge of the oceanic carbon cycle there. The flow of carbon in the Southern Ocean is traditionally understood using a zonal mean framework, in which the meridional overturning circulation drives the latitudinal variability observed in both air–sea flux and interior ocean carbon concentration. However, recent advances, based largely on expanded observation and modeling capabilities in the region, reveal the importance of processes acting at smaller scales, including basin-scale zonal asymmetries in mixed-layer depth, mesoscale eddies, and high-frequency atmospheric variability. Assessing the current state of knowledge and remaining gaps emphasizes the need to move beyond the zonal mean picture and embrace a four-dimensional understanding of the carbon cycle in the Southern Ocean.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Microgrid Assisted Design for Remote Areas

In this work, we present a three-stage multiobjective mixed-integer linear programming (MILP) for the optimal expansion planning and operation of isolated multienergy microgrids in remote areas. By selecting the optimal distributed generators (DGs) and energy storage systems (ESSs) mix selection, siting, sizing, and scheduling in the remote microgrid, the proposed model is targeted to minimize the annualized total cost of microgrids while enhancing the performance of the system, i.e., minimizing the voltage deviations and line power loss. To represent the electricity and heat flow between generation resources and various electrical, heating, and cooling loads in the isolated microgrid, linearized power flow, and heat flow constraints are employed in the proposed optimization model. The available capacity of DGs and ESSs are modeled as discrete constants instead of continuous variables for practical purpose. Numerical simulation results on a remote microgrid consisting of DGs, ESSs, and various loads validate the proposed method.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Shallow-Water Habitat in the Lower Columbia River Estuary: A Highly Altered System

Abstract Decreases in shallow-water habitat area (SWHA) in the Lower Columbia River and Estuary (LCRE) have adversely affected salmonid populations. We investigate the causes by hindcasting SWHA from 1928 to 2004, system-wide, based on daily higher high water (HHW) and system hypsometry. Physics-based regression models are used to represent HHW along the system as a function of river inflow, tides, and coastal processes, and hypsometry is used to estimate the associated SWHA. Scenario modeling is employed to attribute SWHA losses to levees, flow regulation, diversion, navigational development, and climate-induced hydrologic change, for subsidence scenarios of up to 2 m, and for 0.5 m fill. For zero subsidence, the system-wide annual-average loss of SWHA is 55 ± 5%, or 51 × 10 5 ha/year; levees have caused the largest decrease ( $${54}_{-14}^{+5}$$ 54 - 14 + 5 %, or ~ 50 × 10 5 ha/year). The loss in SWHA due to operation of the hydropower system is small, but spatially and seasonally variable. During the spring freshet critical to juvenile salmonids, the total SWHA loss was $${63}_{-3}^{+2}$$ 63 - 3 + 2 %, with the hydropower system causing losses of 5–16% (depending on subsidence). Climate change and navigation have caused SWHA losses of $${5}_{-5}^{+16}$$ 5 - 5 + 16 % and $${4}_{-6}^{+14}$$ 4 - 6 + 14 %, respectively, but with high spatial variability; irrigation impacts have been small. Uncertain subsidence causes most of the uncertainty in estimates; the sum of the individual factors exceeds the total loss, because factors interact. Any factor that reduces mean or peak flows (reservoirs, diversion, and climate change) or alters tides and along-channel slope (navigation) becomes more impactful as assumed historical elevations are increased to account for subsidence, while levees matter less.

59 BASIC BIOLOGICAL SCIENCES↗