Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “GRID MODELING”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

An adaptive grid, unsteady model for two-dimensional magnetohydrodynamic (MHD) flow

An adaptive grid finite difference method for solving multi-dimensional, time-dependent, magnetohydrodynamic (MHD) equations is developed. The method is capable of solving problems that include high gradients due to geometry, propagation of shock waves and unsteady boundary conditions. The grid generation technique is based on variational principles with direct control over grid concentration, smoothness and skewness. An example for a two-dimensional MHD simulation of the propagation of a solar-flare-generated shock wave in solar wind flow in the heliographic equatorial plane is selected for illustration of this method.

Panitchob, Supat↗
DOE OSTI2026

Flexible AI Models for Grid Resilience

The rapid growth in size and complexity of artificial intelligence (AI) and machine learning (ML) models has led to increased energy demands, posing a threat to the reliability of the existing power grid. This project addresses the challenge of highly intermittent and energy-intensive inference workloads by (1) developing fidelity-adaptive neural networks capable of dynamic response to grid conditions and (2) integrating these networks with power flow simulations to assess their impact on power grid reliability. We will explore both top-down and bottom-up approaches to create hierarchies of submodels that provide a controlled trade-off between power draw and prediction accuracy. The top-down method utilizes NN pruning to reduce a flagship model into progressively smaller, energy-efficient variants. The bottom-up approach employs geometrically principled weight setting strategies to construct depth-efficient models from the ground up. A real-time hardware-in-the-loop (HIL) platform will be developed to simulate a scaled AC power grid, integrating live AI workload power draw and enabling dynamic model switching in response to grid feedback. This work will provide a novel framework for evaluating the impact of flexible AI/ML workloads on grid performance and establish new methodologies for energy-aware computing in data centers. The outcomes will demonstrate that adaptive AI/ML can play a critical role in improving grid stability while advancing NREL's leadership in energy-efficient computing research.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Grid Enhancing Technologies Modeling Methodology: Telos Approach

Grid Enhancing Technologies Modeling Methodology Introduction -Grid Enhancing Technology -PLEXOS and TARA Modeling Methodology Overview -PLEXOS/TARA agreement -Flowgate identification -DLR and PFC Location Identification Key Lessons Learned and Possible Improvements

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Quantification of the importance of the modeling fidelity in the Best Estimate Plus Uncertainty analysis

In this work, the multi-physics coupling of CTF-PARCS for both steady-state and transient simulations is evaluated focusing on the quantification of the importance of the modeling fidelity when applied to the Best Estimate Plus Uncertainty safety analysis. For the steady-state, a Picard iteration approach is adopted in the coupled code while for the transient an explicit approach has been developed as a first step. The CTF-PARCS coupling is applied to the REA at HFP BOC scenario proposed in the exercise I-2c of LWR-UAM benchmark Phase III. Two different modeling choices are investigated for the nominal calculations. The first consists of including the 3D burnup information in CTF, and the second of modeling the spacers grids in both PARCS and CTF. It is observed that due to the degradation of the fuel thermal conductivity, the 3D burnup inclusion in CTF leads to significant higher centerline temperatures than the base case, with a difference up to 300 K. The spacer grids modeling shows the expected power and temperature decrease in the axial positions of the spacer grids. For the uncertainty quantification study, the multi-physics modeling including the burnup information in CTF is used. The results show that the inserted reactivity has a relative standard deviation of 13.6% mainly attributed to the nuclear data. The maximum linear power has a relative standard deviation of 16.2% due to the nuclear data as well, while the maximum centerline temperature has a relative standard deviation of 3.5% primarily caused by the fuel thermal conductivity. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2013

Sub-grid drag models for horizontal cylinder arrays immersed in gas-particle multiphase flows

Immersed cylindrical tube arrays often are used as heat exchangers in gas-particle fluidized beds. In multiphase computational fluid dynamics (CFD) simulations of large fluidized beds, explicit resolution of small cylinders is computationally infeasible. Instead, the cylinder array may be viewed as an effective porous medium in coarse-grid simulations. The cylinders' influence on the suspension as a whole, manifested as an effective drag force, and on the relative motion between gas and particles, manifested as a correction to the gas-particle drag, must be modeled via suitable sub-grid constitutive relationships. In this work, highly resolved unit-cell simulations of flow around an array of horizontal cylinders, arranged in a staggered configuration, are filtered to construct sub-grid, or `filtered', drag models, which can be implemented in coarse-grid simulations. The force on the suspension exerted by the cylinders is comprised of, as expected, a buoyancy contribution, and a kinetic component analogous to fluid drag on a single cylinder. Furthermore, the introduction of tubes also is found to enhance segregation at the scale of the cylinder size, which, in turn, leads to a reduction in the filtered gas-particle drag.

Computational fluid dynamics (CFD), Cylinders, Fil↗
DOE OSTI2021

Sensor enabled data-driven predictive analytics for modeling and control with high penetration of DERs in distribution systems

The electric power grid is undergoing a tremendous transformation due to the increasing penetration of renewable energy resources beginning with wind and more recently with the distributed energy resources (DERs) such as solar and battery storage. DERs have dramatically changed the role of the distribution systems in the overall power grid, and they are expected to contribute a significant portion of power generation in the future. If current trends for DERs continue, system operation and control will need to change dramatically for improved grid reliability and resiliency. As renewable resources increase in penetration, new and challenging operational, planning, and design problems are expected to emerge. Some of the key challenges that arise in the planning and operation of the future grid are: 1) Quantifying the impact of high DER penetration in distribution systems on bulk grid behavior over multiple time scales. 2) Identifying whether a particular DER configuration/settings have a large impact on the overall grid behavior. These challenges can be addressed in an offline manner using detailed T&D grid models and they can also be addressed in an online manner using sensor measurements. In particular, the advancement and planned growth in sensor technology in power grid over various voltage levels provide us with a unique opportunity to tackle these challenges from a data analytic perspective without needing detailed T&D grid models. A few questions that naturally arise when addressing the challenges from DERs using sensor data are: 1) How can we use limited sensor measurements to monitor & control voltage stability and small signal stability of the bulk system? 2) How can we ensure that the developed data analytic methods are robust to data availability and quality issues? 3) How can we compute the developed analytics in a scalable manner using streaming measurements? In this project, we addressed the aforementioned challenges arising from DERs and answered the questions raised above on how to effectively use the sensor measurements to enhance the reliability and performance of the electric grid. Thus, the overarching goal of this project is to develop effective reduced/representative system models from data that make the computational complexity sufficiently manageable so as to be useful to simulate, analyze, and even control complex non-linear power systems dynamics with large penetrations of DERs. In order to achieve the objective, the project team established a four-fold technical approach 1) Formulated a combined transmission-distribution co-simulation framework for data generation and validation, 2) Derived reduced/representative models of power systems based on data-driven methods for efficient computation and appropriate representation of system behavior, 3) Developed data driven characterization of power system behavior based on transfer operator theory, machine learning and optimization for model estimation, 4) Incorporated a scalable data management and processing architecture using distributed Kafka streaming applications that coordinate input data streams to the developed data analytics. The key accomplishments of the project are: 1) Development of a scalable multi-timescale T&D co-simulation framework (both for steady state and for dynamic co-simulation) using commercial solvers (PSSE and GridLAB-D). The steady-state T&D co-simulation interface is shared with our industry partner (PJM). 2) A structured reduced order dynamic model of distribution systems that can represent partial motor stalling along with a systematic procedure to derive the model parameters. 3) A PMU based online method to monitor, localize and mitigate fault-induced delayed voltage recovery using DER reactive support and load control in distribution systems. 4) Development of linear operator based robust methodologies for dynamic state estimation, uncertainty quantification, system identification and trajectory prediction for power system dynamics. 5) An adaptive damping control for utilizing wind energy resources to provide oscillation damping and system stability. 6) Implementation of Kafka-based framework for efficient processing of streaming data using Linux-based local virtual environment.

DER integration↗
DOE OSTI2017

Reducing Cascading Failure Risk by Increasing Infrastructure Network Interdependence

Increased coupling between critical infrastructure networks, such as power and communication systems, has important implications for the reliability and security of these systems. To understand the effects of power-communication coupling, several researchers have studied models of interdependent networks and reported that increased coupling can increase vulnerability. However, these conclusions come largely from models that have substantially different mechanisms of cascading failure, relative to those found in actual power and communication networks, and that do not capture the benefits of connecting systems with complementary capabilities. In order to understand the importance of these details, this paper compares network vulnerability in simple topological models and in models that more accurately capture the dynamics of cascading in power systems. First, we compare a simple model of topological contagion to a model of cascading in power systems and find that the power grid model shows a higher level of vulnerability, relative to the contagion model. Second, we compare a percolation model of topological cascading in coupled networks to three different models of power networks coupled to communication systems. Again, the more accurate models suggest very different conclusions than the percolation model. In all but the most extreme case, the physics-based power grid models indicate that increased power-communication coupling decreases vulnerability. This is opposite from what one would conclude from the percolation model, in which zero coupling is optimal. Only in an extreme case, in which communication failures immediately cause grid failures, did we find that increased coupling can be harmful. Together, these results suggest design strategies for reducing the risk of cascades in interdependent infrastructure systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2017

1–2.4 μm Near-IR Spectrum of the Giant Planet β Pictoris b Obtained with the Gemini Planet Imager

Using the Gemini Planet Imager located at Gemini South, we measured the near-infrared (1.0–2.4 μm) spectrum of the planetary companion to the nearby, young star β Pictoris. We compare the spectrum obtained with currently published model grids and with known substellar objects and present the best matching models as well as the best matching observed objects. Comparing the empirical measurement of the bolometric luminosity to evolutionary models, we find a mass of 12.9 ± 0.2 M{sub Jup}, an effective temperature of 1724 ± 15 K, a radius of 1.46 ± 0.01 R{sub Jup}, and a surface gravity of logg=4.18±0.01 [dex] (cgs). The stated uncertainties are statistical errors only, and do not incorporate any uncertainty on the evolutionary models. Using atmospheric models, we find an effective temperature of 1700–1800 K and a surface gravity of logg=3.5–4.0 [dex] depending upon the model. These values agree well with other publications and with “hot-start” predictions from planetary evolution models. Further, we find that the spectrum of β Pic b best matches a low surface gravity L2 ± 1 brown dwarf. Finally, comparing the spectrum to field brown dwarfs, we find the the spectrum best matches 2MASS J04062677–381210 and 2MASS J03552337+1133437.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Future Grid State Modeling for Transactive Systems

The simulation design and evaluation of a transactive energy system is contingent on the setting in which it would exist, which could vary significantly with the many drivers that impact the potential future states of the power system. This report summarizes the development of the Transactive Future Grid State (TFGS) model, which is a lightweight capacity expansion model that will support transactive energy research by allowing users to explore alternative future states. The TFGS model determines the long-term generation capacity needs of a given power system that satisfy specified constraints and support a least cost hourly dispatch. This is accomplished by solving a single optimization problem that simultaneously finds the generation fleet with capacity sufficient to meet peak load and the flexibility to meet hourly operational dispatch, all at lowest cost. The TFGS model was applied to a larger ongoing study taking place that will be simulating how distribution system operators (DSOs) could use a transactive energy system to manage distributed energy resources. The larger study, commonly referred to as the DSO+T study, involves two parallel analyses; one that represents the system today and one that represents a future system with a high degree of renewable energy penetration. The TFGS model provided the generation mix of the high renewable base case. As the TFGS model was being used to provide inputs to the DSO+T study, two issues became apparent. First, additional functionality would have been useful, and second, a larger-than-anticipated amount of user inputs and data were required to successfully use the TFGS model. These discoveries have lead to the recommendation of specific functionality that can be added in future development of the tool.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS1999

Navier-Stokes Simulation of Several High-Lift Reference H Configurations

The subsonic flow field was numerically simulated around several High Speed Research Reference H configurations at various pitch and yaw angles. A sequence of structured-viscous grids were generated; the first grid modeled the wing-body high-lift geometry, and the second grid incorporated the nacelles and the horizontal tail. The third grid modeled the full-span geometry for sideslip calculations, and was obtained by mirroring a coarser version of the second grid. The CFL3D code, a Reynolds averaged, thin-layer Navier-Stokes flow solver for structural grids, was used for the flow solver and modeled the free-air Reference H high-lift configuration at wind tunnel conditions of Mach number 0.24 and Reynolds number of 1.4 x 10(exp 5) per in. Pitch sweeps were performed at angles of attack from 6 deg to 15 deg. Sideslip angle sweeps at 0 deg <= Beta <= +18 deg were performed at an angle of attack of 8 deg. The lateral and longitudinal performance characteristics were well predicted and very good force and moment comparisons were obtained. A very complex multiple vortical system develops at the higher angles of attack, and detailed postprocessing of the solutions provided a comprehensive three-dimensional understanding of the flow which helps to correlate and interpret the wind tunnel data.

Lessard, Wendy B.↗
NASA NTRS2022

The Effects of Surface Fossil Magnetic Fields on Massive Star Evolution: IV. Grids of Models at Solar, LMC, and SMC Metallicities

Magnetic fields can drastically change predictions of evolutionary models of massive stars via mass-loss quenching, magnetic braking, and efficient angular momentum transport, which we aim to quantify in this work. We use the MESA software instrument to compute an e xtensiv e main-sequence grid of stellar structure and evolution models, as well as isochrones, accounting for the effects attributed to a surface fossil magnetic field. The grid is densely populated in initial mass (3–60 M ), surface equatorial magnetic field strength (0–50 kG), and metallicity (representative of the Solar neighbourhood and the Magellanic Clouds). We use two magnetic braking and two chemical mixing schemes and compare the model predictions for slowly rotating, nitrogen- enriched (‘Group 2’) stars with observations in the Large Magellanic Cloud. We quantify a range of initial field strengths that allow for producing Group 2 stars and find that typical values (up to a few kG) lead to solutions. Between the subgrids, we find notable departures in surface abundances and evolutionary paths. In our magnetic models, chemical mixing is al w ays less efficient compared to non-magnetic models due to the rapid spin-down. We identify that quasi-chemically homogeneous main sequence evolution by efficient mixing could be prevented by fossil magnetic fields. We recommend comparing this grid of evolutionary models with spectropolarimetric and spectroscopic observations with the goals of (i) revisiting the derived stellar parameters of known magnetic stars, and (ii) observationally constraining the uncertain magnetic braking and chemical mixing schemes.

stars: abundances↗
NASA NTRS2013

A Storm Surge and Inundation Model of the Back River Watershed at NASA Langley Research Center

This report on a Virginia Institute for Marine Science project demonstrates that the sub-grid modeling technology (now as part of Chesapeake Bay Inundation Prediction System, CIPS) can incorporate high-resolution Lidar measurements provided by NASA Langley Research Center into the sub-grid model framework to resolve detailed topographic features for use as a hydrological transport model for run-off simulations within NASA Langley and Langley Air Force Base. The rainfall over land accumulates in the ditches/channels resolved via the model sub-grid was tested to simulate the run-off induced by heavy precipitation. Possessing both the capabilities for storm surge and run-off simulations, the CIPS model was then applied to simulate real storm events starting with Hurricane Isabel in 2003. It will be shown that the model can generate highly accurate on-land inundation maps as demonstrated by excellent comparison of the Langley tidal gauge time series data (CAPABLE.larc.nasa.gov) and spatial patterns of real storm wrack line measurements with the model results simulated during Hurricanes Isabel (2003), Irene (2011), and a 2009 Nor'easter. With confidence built upon the model's performance, sea level rise scenarios from the ICCP (International Climate Change Partnership) were also included in the model scenario runs to simulate future inundation cases.

Loftis, Jon Derek↗
DOE OSTI2016

Detailed high-resolution three-dimensional simulations of OMEGA separated reactants inertial confinement fusion experiments

Here, we present results from the comparison of high-resolution three-dimensional (3D) simulations with data from the implosions of inertial confinement fusion capsules with separated reactants performed on the OMEGA laser facility. Each capsule, referred to as a “CD Mixcap,” is filled with tritium and has a polystyrene (CH) shell with a deuterated polystyrene (CD) layer whose burial depth is varied. In these implosions, fusion reactions between deuterium and tritium ions can occur only in the presence of atomic mix between the gas fill and shell material. The simulations feature accurate models for all known experimental asymmetries and do not employ any adjustable parameters to improve agreement with experimental data. Simulations are performed with the RAGE radiation-hydrodynamics code using an Implicit Large Eddy Simulation (ILES) strategy for the hydrodynamics. We obtain good agreement with the experimental data, including the DT/TT neutron yield ratios used to diagnose mix, for all burial depths of the deuterated shell layer. Additionally, simulations demonstrate good agreement with converged simulations employing explicit models for plasma diffusion and viscosity, suggesting that the implicit sub-grid model used in ILES is sufficient to model these processes in these experiments. In our simulations, mixing is driven by short-wavelength asymmetries and longer-wavelength features are responsible for developing flows that transport mixed material towards the center of the hot spot. Mix material transported by this process is responsible for most of the mix (DT) yield even for the capsule with a CD layer adjacent to the tritium fuel. Consistent with our previous results, mix does not play a significant role in TT neutron yield degradation; instead, this is dominated by the displacement of fuel from the center of the implosion due to the development of turbulent instabilities seeded by long-wavelength asymmetries. Through these processes, the long-wavelength asymmetries degrade TT yield more than the DT yield and thus bring DT/TT neutron yield ratios into agreement with experiment. Finally, we present a detailed comparison of the flows in 2D and 3D simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1991

Comparison of 2-D and 3-D models of grid erosion in an ion thruster

Numerical results of particle-in-cell/Monte Carlo calculations of accelerator grid erosion in an ion thruster are presented. Specifically, it is shown that a three-dimensional model is required to account for the experimentally observed pitting of the accelerator grid between grid apertures. Some comparisons with earlier two-dimensional, axisymmetric model are made, and it is shown that, for identical operating conditions of the thruster, the wear-through time in the three-dimensional model is about two to three times higher than that obtained previously with the two-dimensional model, namely on the order of 10,000 hours for sample calculation.

Peng, Xiaohang↗
NASA NTRS2004

Surface Modeling and Grid Generation for Iced Airfoils (SmaggIce)

Many of the troubles associated with problem solving are alleviated when there is a model that can be used to represent the problem. Through the Advanced Graphics and Visualization (G-VIS) Laboratory and other facilities located within the Research Analysis Center, the Computer Services Division (CSD) is able to develop and maintain programs and software that allow for the modeling of various situations. For example, the Icing Research Branch is devoted to investigating the effect of ice that forms on the wings and other airfoils of airplanes while in flight. While running tests that physically generate ice and wind on airfoils within the laboratories and wind tunnels on site are done, it would be beneficial if most of the preliminary work could be done outside of the lab. Therefore, individuals from within CSD have collaborated with Icing Research in order to create SmaggIce. This software allows users to create ice patterns on clean airfoils or open files containing a variety of icing situations, manipulate and measure these forms, generate, divide, and merge grids around these elements for more explicit analysis, and specify and rediscretize subcurves. With the projected completion date of Summer 2005, the majority of the focus of the Smagglce team is user-functionality and error handling. My primary responsibility is to test the Graphical User Interface (GUI) in SmaggIce in order to ensure the usability and verify the expected results of the events (buttons, menus, etc.) within the program. However, there is no standardized, systematic way in which to test all the possible combinations or permutations of events, not to mention unsolicited events such as errors. Moreover, scripting tests, if not done properly and with a view towards inevitable revision, can result in more apparent errors within the software and in effect become useless whenever the developers of the program make a slight change in the way a specific process is executed. My task therefore requires a brief yet intense study into GUI coverage criteria and creating algorithms for GUI implementation. Nevertheless, there are still heavily graphical features of SmaggIceSmaggIce that must be either corrected or redesigned before its release. A particular feature of SmaggIce is the ability to smooth out curves created by control points that form an arbitrary shape into something more acquiescent to gridding (while maintaining the integrity of the data). This is done by a mathematical model known as Non-Uniform Rational B-Spline (NURBS) curves. Existing NURBS code is written in FORTRAN-77 with static arrays for holding information. My new assignment is to allow for dynamic memory allocation within the code and to make it possible for the developers to call out functions from the NURBS code using C.

Hammond, Brandy M.↗
DOE OSTI2024

WTK-LED: The WIND Toolkit Long-Term Ensemble Dataset

To satisfy a wide group of stakeholders across various wind energy disciplines, including but not limited to stakeholders in the distributed and utility scale wind industry, the new emerging airborne wind energy field, grid integration, power systems modeling, environmental modeling, and researchers in academia, and to close some of the gaps that current public datasets have, we aimed at developing an updated version of the meteorological WIND Toolkit, named WIND Toolkit Long-term Ensemble Dataset (WTK-LED), which is a meteorological dataset providing time series every 5 min and 2 km, including model uncertainty of wind speed at every modeling grid point so that users are provided with a range of possible wind speeds every 2 km. The data were produced using the Weather Research and Forecasting Model (WRF). The vertical grid used in WTK-LED includes many vertical layers in the atmospheric boundary layer to provide information of atmospheric quantities across the rotor layer of utility scale and distributed wind turbines. The WTK-LED includes: 1) Numerical simulations covering the continental United States, Alaska, and Hawaii, with high-resolution data being available for 3 years (2018-2020). 2) Climate simulations from Argonne National Laboratories covering the North American continent, including Alaska, Canada, and most of Mexico and the Caribbean Islands. These simulations complement the new WTK-LED to offer a 4-km dataset covering 20 years, from 2001-2020. 3) Specific long-term,high-resolution offshore simulations have been conducted separately for the US coasts, Hawaii, and the Great Lakes, leading to the 2023 National Offshore Wind data set. This report focuses on a description of the land-based WTK-LED for CONUS, Hawaii, and Alaska, for the 3-year 2-km/5-min dataset and the 20-year 4-km/hourly dataset, as well as the uncertainty quantification method. We also provide limited validation results. Based on our results to date, we suggest use cases and applications for each dataset of the WTK-LED.

17 WIND ENERGY↗
DOE OSTI2026

HydraGNN_OPF_GFM_2026 - Ensemble of predictive graph foundation models for power grid applications

This dataset supports research on graph foundation models for optimal power flow (OPF) on electric grids using HydraGNN. It contains heterogeneous graph representations of PGLib-OPF cases spanning systems from 14 to 13,659 buses, together with packed HDF5 datasets for pretraining, feasibility classification, and N-1 contingency analysis. The release includes OPF solution data, downstream fine-tuning datasets, pretrained HeteroSAGE and HeteroHEAT model checkpoints, hyperparameter-optimization summaries across multiple heterogeneous GNN architectures, and aggregated fine-tuning results for sample-efficiency studies. The dataset is designed to enable scalable training, evaluation, and transfer-learning studies for OPF surrogate modeling, including node-level AC-OPF solution prediction, graph-level prediction, feasibility classification, operating-condition generalization, and contingency-response tasks.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Explicit quantification of residence and flushing times in the Salish Sea using a sub-basin scale shoreline resolving model

The Salish Sea, located in the Pacific Northwest region of North America has complex currents and circulation features distributed over numerous interconnected deep basins with islands. Increased risk of exposure to oil spills and untreated wastewater from maritime emergencies and treatment plant failures have led to a need for quantifying residence and flushing characteristics at a sub-basin scale using a shoreline resolving hydrodynamic model. An unstructured grid model of the Salish Sea was developed using FVCOM with a ≈75-100m shoreline resolution. In addition to 7 tides and 23 salinity and temperature monitoring stations, an extensive currents data set from 135 stations collected over a span of three years was used for skill assessment and validation. Explicit forward computations were then conducted to define and quantify residence and flushing times in various sub-basins of interest using (a) Lagrangian particles and (b) Numerical/virtual dye experiments. The results in most basins show expected seasonal variability with longer flushing time associated with summer lower tides and lower freshwater inflows. However, contrary to expectation, flushing time is significantly longer in wintertime in large fjord-like basins such as Hood Canal (≈138 days), likely due to increased stratification and reduced mixing. The flushing time for the Puget Sound region of the Salish Sea is ≈115 days, while Georgia Basin is 240 days when analyzed as stand-alone basins with zero background concentrations. When examined as part of the flushing of the entire Puget Sound filled with virtual dye, the compounded flushing times for embedded sub-basins can be significantly longer in order of magnitudes and largely dictated by the flushing time of Puget Sound. The computed residence and flushing time scales tabulated over 36 sub-basins provide an improved understanding of water renewal in the system, informing pollution management actions.

54 ENVIRONMENTAL SCIENCES↗