Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical modelling”

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 613 records · Page 34

Mesoscopic finite-size effects of unconventional electron transport in PdCoO 2

A wide range of unconventional transport phenomena has recently been observed in single-crystal delafossite metals. Here, we present a theoretical framework to elucidate electron transport using a combination of first-principles calculations and numerical modeling of the anisotropic Boltzmann transport equation. Using PdCoO 2 as a model system, we study different microscopic electron and phonon scattering mechanisms and establish the mean free path hierarchy of quasiparticles at different temperatures. We treat the anisotropic Fermi surface explicitly to numerically obtain experimentally-accessible transport observables, which bridge between the “diffusive,” “ballistic,” and “hydrodynamic” transport regime limits. Further, we illustrate that the distinction between the “quasiballistic” and “quasihydrodynamic” regimes is challenging and often needs to be quantitative in nature. From first-principles calculations, we populate the resulting transport regime plots and demonstrate how the Fermi surface orientation adds complexity to the observed transport signatures in micrometer-scale devices. Our work provides key insights into microscopic interaction mechanisms on open hexagonal Fermi surfaces and establishes their connection to the macroscopic electron transport in finite-size channels.

36 MATERIALS SCIENCE↗

Definition of the Stiesdal Offshore TetraSpar Floating Wind System for OC6 Phase IV

The Offshore Code Comparison Collaboration, Continued, with Correlation, and unCertainty (OC6) is an international project focused on verifying and validating modeling tools used in the analysis of floating offshore wind systems. The fourth phase of the project aims to assess the ability of the tools to model a novel floating wind design. The TetraSpar was designed by Stiesdal Offshore Technologies and was selected as the subject platform for the project. The platform is a unique design that features spar-like stability with a very low center of gravity, but does not require rare deep-water ports. This combination is achieved with a system made of two separate steel structures, a hull at the water's surface that supports the wind turbine tower, and a keel that is suspended below the main body. This lower keel has a very high density and drives the total platform's low center of gravity. The TetraSpar, with hull and keel, spar-type stability, and relatively slender components, is unique compared to platforms studied in earlier phases of OC3-OC6 projects, and it offers an opportunity to evaluate the applicability of modeling tools to a novel design, especially unique hydrodynamics and compliant cables. Scaled model tests of the TetraSpar design were performed in 2018-2019 by the University of Maine. The data collected in the model tests provides a detailed description of the response of the platform in various controlled and relatively well understood conditions. Comparison with this data provides a good opportunity for validation of numerical models. The OC6 Phase IV work focuses on the specific parameters used in these scaled model tests. Descriptions of the platform and the test campaign in this document are based on the report written by Allen and Fowler of the University of Maine in 2019.

17 WIND ENERGY↗

A Review of Modeling Approaches for Understanding and Monitoring the Environmental Effects of Marine Renewable Energy

Understanding the environmental effects of marine energy (ME) devices is fundamental for their sustainable development and efficient regulation. However, measuring effects is difficult given the limited number of operational devices currently deployed. Numerical modeling is a powerful tool for estimating environmental effects and quantifying risks. It is most effective when informed by empirical data and coordinated with the development and implementation of monitoring protocols. We reviewed modeling techniques and information needs for six environmental stressor–receptor interactions related to ME: changes in oceanographic systems, underwater noise, electromagnetic fields (EMFs), changes in habitat, collision risk, and displacement of marine animals. This review considers the effects of tidal, wave, and ocean current energy converters. We summarized the availability and maturity of models for each stressor–receptor interaction and provide examples involving ME devices when available and analogous examples otherwise. Models for oceanographic systems and underwater noise were widely available and sometimes applied to ME, but need validation in real-world settings. Many methods are available for modeling habitat change and displacement of marine animals, but few examples related to ME exist. Models of collision risk and species response to EMFs are still in stages of theory development and need more observational data, particularly about species behavior near devices, to be effective. We conclude by synthesizing model status, commonalities between models, and overlapping monitoring needs that can be exploited to develop a coordinated and efficient set of protocols for predicting and monitoring the environmental effects of ME.

16 TIDAL AND WAVE POWER↗

Quasi‐dynamic breathing model of the lung incorporating viscoelasticity of the lung tissue

Abstract We advanced a novel model to calculate viscoelastic lung compliance and airflow resistance in presence of mucus, accounting for the quasi‐linear viscoelastic stress–strain response of the parenchyma (alveoli) tissue. We adapted a continuum‐based numerical modeling approach for the lung, integrating the fluid mechanics of the airflow within individual generations of the bronchi and alveoli. The model accounts for elasticity of the deformable bronchioles, resistance to airflow due to the presence of mucus within the bronchioles, and subsequent mucus flow. Simulated quasi‐dynamic inhalation and expiration cycles were used to characterize the net compliance and resistance of the lung, considering the rheology of the mucus and viscoelastic properties of the parenchyma tissue. The structure and material properties of the lung were identified to have an important contribution to the lung compliance and airflow resistance. The secondary objective of this work was to assess whether a higher frequency and smaller volume of harmonic air flow rate compared to a normal ventilator breathing cycle enhanced mucus outflow. Results predict, lower mucus viscosity and higher excitation frequency of breathing are favorable for the flow of mucus up the bronchi tree, towards the trachea.

60 APPLIED LIFE SCIENCES↗

Technical Justifications for Liquid Hydrogen Exposure Distances

The previous separation distances in the National Fire Protection Association (NFPA) Hydrogen Technologies Code (NFPA 2, 2020 Edition) for bulk liquid hydrogen systems lack a well-documented basis and can be onerous. This report describes the technical justifications for revisions of the bulk liquid hydrogen storage setback distances in NFPA 2, 2023 Edition. Distances are calculated based on a leak area that is 5% of the nominal pipe flow area. Models from the open source HyRAM+ toolkit are used to justify the leak size as well as calculate consequence-based separation distances from that leak size. Validation and verification of the numerical models is provided, as well as justification for the harm criteria used for the determination of the setback distances for each exposure type. This report also reviews mitigations that could result in setback distance reduction. The resulting updates to the liquid hydrogen separation distances are well-documented, retrievable, repeatable, revisable, independently verified, and use experimental results to verify the models.

08 HYDROGEN↗

Using Parameter Sweep in WaterTAP to Analyze New Water Treatment Technologies

We describe a powerful and generalized parameter sweep tool in this report that was originally developed to analyze the performance of existing and novel water treatment models being developed in WaterTAP. Since WaterTAP is built upon IDAES and Pyomo, the parameter sweep tool can be used to systematically explore and debug the behavior of most Pyomo and IDAES numerical models. In order to enable meaningful analyses, the parameter sweep tool has been designed with the following features: 1) Model flexibility: The parameter sweep tool does not enforce any restrictions on the types of models that can be used with it. As long as a Pyomo model can be solved and the parameter is active and mutable, the tool only needs functions that describe how to run the model, the sweep parameters, and the output quantities of interest. 2) Flexible sampling: The parameter sweep tool has inbuilt functions to generate samples from a random distribution or a multidimensional Euclidean space. Furthermore, the users have to ability to supply samples generated from a tool of their choice. 3) Multiple sweep types: A user can choose from one of 3 types of parameter sweeps depending on their needs. 4) Detailed outputs: Outputs generated by the parameter sweep tool can be stored in detailed H5 file or user-friendly CSV files for post processing. 5) Parallel computing: The parameter sweep supports shared and distributed memory parallel computing to enable the use of high performance computers (HPC) for large-scale analyses. 6) Modular: The parameter sweep tool is self-contained and can easily be integrated within an outer-loop analysis or as desired by the user. 7) Ease of use: The tool is well documented and a simple sweep can be easily executed by following the online documentation in a few lines of code. We demonstrate the use of the parameter sweep tool on a simple water treatment system from the WaterTAP repository and show its parallel scaling performance on an Apple laptop and NREL's Eagle HPC. The parameter sweep tool is actively being used with models currently being developed within WaterTAP and we expect its use to grow beyond it to other IDAES and Pyomo models.

97 MATHEMATICS AND COMPUTING↗

Verification and Validation of Model-Scale Turbine Performance and Control Strategies for the IEA Wind 15 MW Reference Wind Turbine

To enable the fast growth of the floating offshore wind industry, simulation models must be validated with experimental data. Floating wind model-scale experiments in wind–wave facilities have been performed over the last two decades with varying levels of fidelity and limitations. However, the turbine controls in these experiments have considered only limited control strategies and implementations. To allow for control co-design, this research focuses on implementing and experimentally validating more advanced turbine control actions and strategies in a wind–wave basin for a 1:70-scale model of the International Energy Agency’s wind 15 MW reference wind turbine. The control strategies analyzed include torque control, collective pitch control, and transition region control (setpoint smoothing). Our experimental and numerical results include the effects of varying rotor speeds, blade pitches, and wind environments on the turbine thrust and torque. Numerical models from three different software tools are presented and compared to the experimental results. Their ability to effectively represent the aero-dynamic response of the wind turbine to the control actions is successfully validated. Finally, turbine controller tuning parameters based on the derivatives of thrust and torque are derived to allow for improved offshore wind turbine dynamics and to validate the ability of modeling tools to model the dynamics of floating offshore wind turbines with control co-design.

17 WIND ENERGY↗

2.1.3.404 - WEC Array Power Management and Output Simulation Tool

An array of wave energy converter (WEC) devices has variations in power output due to the chaotic nature of the waves. Eliminating or mitigating the power fluctuations is important for reducing the integration impacts of WEC plants in both distribution and transmission grids, and in standalone isolated power systems. Reduced variability of WEC-generated power in combination with energy storage or power management control at each WEC and at the array level will help increasing hosting capacity of distribution feeders for this type of variable renewable generation, and minimization of electric losses. To reduce risks and risk perception, gain key stakeholder acceptance, and enable developers to design effective and compatible energy plants with arrays of WEC devices, it requires the use of modeling tools to simulate the resource environment, device dynamics, utility power system response, as well as the development of an interface for an array controller. The project will create a publicly accessible numerical modeling framework to empower the wave energy sector to design projects of various scales (kW-100s MW), which are optimized on a plant performance basis and are compatible with different power systems and wave conditions. The framework will integrate with WEC-Sim, a wave environment model (SWAN-FUNWAVE), as well as established relevant electrical analysis tools, to model the grid system and interconnection, and optimize power output and power management for the WEC array.

POWER TRANSMISSION AND DISTRIBUTION,TIDAL AND WAVE↗

Development of a reduced model for energetic particle transport by sawteeth in tokamaks

We report that the sawtooth instability is known for inducing transport and loss of energetic particles (EPs), and for generating seed magnetic islands that can trigger tearing modes. Both effects degrade the overall plasma performance. Several theories and numerical models have been previously developed to quantify the expected EP transport caused by sawteeth, with various degrees of sophistication to differentiate the response of EPs at different energies and on different orbits (e.g. passing vs. trapped), although the analysis is frequently limited to a single time slice during a tokamak discharge. This work describes the development and initial benchmark of a framework that enables a reduced model for EP transport by sawteeth retaining the full EP phase-space information. The model, implemented in the ORBIT hamiltonian particle-following code, can be used either as a standalone post-processor taking input data from codes such as TRANSP, or as a pre-processor to compute transport coefficients that can be fed back to TRANSP for time-dependent simulations including the effects of sawteeth on EPs. The advantage of the latter approach is that the evolution of the EP distribution can be simulated quantitatively for sawtoothing discharges, thus enabling a more accurate modeling of sources, sinks and overall transport properties of EP and thermal plasma species for comprehensive physics studies that require detailed information of the fast-ion distribution function and its evolution over time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Development of a triple-moment ice-phase cloud microphysics scheme and its application to the Single Column Atmosphere Model

Parameterization of cloud microphysics is critical for accurate simulation of weather and climate, in which the characteristics of cloud particle spectrum inevitably further affect climate simulation by changing the cloud evolution and cloud radiation effects. The popular currently used double-moment cloud microphysics schemes in numerical models can predict the intercept (N 0 ) and slope (λ) parameters of cloud particle spectrum but cannot predict the spectral shape parameter (μ), which hinders accurate description of cloud physical processes in climate models. Therefore, in the present study, we built upon the ideas of previously developed triple-moment cloud microphysics scheme, considered radar reflectivity factor as the third predictor in addition to number and mass concentration, and introduced its related prediction equation to the Single Column Atmosphere Model Version 5.3 (SCAM5.3). Moreover, the relevant microphysical process formulas in the model were revised, and a triple-moment ice-cloud microphysics scheme was constructed to predict the μ of ice particle spectrum. Based on this model, a 29-day case of the Atmospheric Radiation Measurement Program in the summer of 1997 (ARM97) was simulated, and the differences in cloud fraction and radiation simulations between the double- and triple-moment schemes were analysed. The μ of ice particle spectrum predicted by the triple-moment scheme mainly ranged from 0 to 4; the peak value was around 2, and at least 85% of the μ values were greater than the default 0. Therefore, the μ=0 setting in the double-moment scheme is unreasonable. Moreover, the developed triple-moment ice-cloud microphysics scheme yielded a narrower ice particle spectrum, which was closer to the observation results of previous studies, than the double-moment cloud microphysics scheme (e.g., height=233 hPa). Furthermore, compared with the double-moment scheme, the triple-moment scheme achieved closer simulations of the cloud fraction observations, particularly for ice clouds in the upper levels. Considering the close association between cloud fraction and radiation, the error with the triple-moment scheme was smaller than that with the double-moment scheme regardless of the shortwave (surface downward shortwave radiation, net downward shortwave radiation at the top of the atmosphere, and net shortwave radiation at the surface) or longwave (surface downward longwave radiation and net longwave radiation at the surface) flux density. Finally, the improvement mechanism of the triple-moment scheme on cloud fraction and radiation simulations was explored. The major reason is that the triple-moment scheme weakens the autoconversion of ice crystal to form snow and enhances the growth process of ice crystal deposition, thereby increasing the mass concentration of ice crystals and ice cloud fraction. Overall, the developed triple-moment ice-cloud microphysics scheme can help improve the simulation ability of models for cloud-climate feedback and other processes.

54 ENVIRONMENTAL SCIENCES↗

Coupled Experimental and High-Temperature Discrete-Element Method Modeling Studies of Aluminosilicate Particle Handling in Concentrated Solar Power Environments

Chemically inert, aluminosilicate based particles have been investigated as both a thermal transport and sensible energy storage medium for concentrated solar power facilities. These particles will experience a wide range of operating temperatures (300-1000 K) and handling conditions (dense to dilute falling particle curtains, dense granular flows, or dense structures), requiring specially-designed and optimized infrastructures. The relative influence of collisional and frictional interactions between particles varies based on temperature-dependent particulate properties and greatly impacts the bulk, granular flow behavior. These underlying physics are captured using discrete element method modeling tools. However, this modeling method is computationally expensive as each particle position and interaction is tracked during the simulation. These modeling methods are further complicated by introducing temperature-dependent particle properties, high-temperature radiative exchange, and directional irradiation sources experienced by granular flows in concentrated solar power environments. Coupled experimental and numerical studies of aluminosilicate particles in rotary kilns and dense particle curtains were performed for bulk temperatures up to 1073 K. The three particle types investigated included Carbobead HSP 30 /60, Carbobead CP 30/60, and Granusil 4030. Temperature, spatial, and velocity profile data were extracted from experimental runs using embedded K-type thermocouple probes and particle image velocimetry techniques. Experimental and numerical studies were compared using spatial temperature profiles, velocity fields, and shape profiles of the bulk, granular flows. Numerical models were developed using commercially available discrete element method modeling software, Aspherix®. Existing Aspherix® functionality was expanded by introducing coupled radiative exchange modeling tools. The laboratory-scale rotary kiln was developed to investigate the steady-state heat and mass transfer performance of aluminosilicate particles based on particle type, bulk handling temperature, and wall roughness. The rotational speed of the rotary kiln was varied to control the relative impact of collisional and frictional effects upon the granular flow behavior. Heat and mass transfer performance was categorized based on the Froude number and the observed flow regimes of slipping, rolling, cascading, and centrifuging. Coupled discrete element method modeling studies were used to evaluate the effects of temperature-dependent, particulate mechanical properties upon bulk flow behavior and upon the relative effects of radiative, advective, and/or conductive heat transfer. A high-temperature (< 1073 K) falling particle curtain was similarly fabricated to investigate the heat and mass transfer performance of aluminosilicate particles in particle handling situations dominated by inter-particle collisions. The impact of particle type, flow preheat temperatures (< 1073K), and bulk mass flow rates were investigated upon the particle curtain shape, temperature, and velocity profiles. Coupled discrete element method modeling studies were performed to evaluate the varying impact of temperature-dependent, particulate mechanical properties on the bulk flow behavior and the temperature profile of the particle curtain.

14 SOLAR ENERGY↗

A Bottom-up Approach to Improve the Representation of Deep Convective Clouds in Weather and Climate Models (Final Report)

A fundamental yet outstanding question about deep cumulus convection regards how convective updrafts, downdrafts, and surface-based cold pools are inter-related within organized storms, and how these three convective components are modulated by external and internal factors. We considered this question through numerical model simulations as well as data collected during the MC3E field campaign. One of our basic hypotheses was that deep cold pools are forced by wide downdrafts which arise from wide updrafts. This “bottom-up” perspective required consideration of the external controls on updraft width, and then the internal controls on cold-pool characteristics including depth. Model experiments showed that environmental vertical wind shear (an external factor) exerts a particularly strong control on updraft width; analyses of linear and nonlinear dynamic pressure forcing readily explain this control. The experiments also showed that increases in updraft width are accompanied by corresponding increases in downdraft width and cold-pool depth. The updraft–downdraft connection is through precipitation: the wide area of precipitation that will be initiated in a wide updraft will in turn lead to a wide downdraft. The downdraft–cold pool connection is through mass continuity: at the downdraft terminus, the large downdraft mass flux through a wide downdraft must lead to a deep cold pool. Model experiments also showed that the cold-pool depth and other characteristics can be significantly modulated by microphysical processes (internal factors). Specifically, the cold-pool expansion rate, depth, and intensity was found to depend most strongly on the amount of rain evaporation near the ground. However, the sustenance of the cold pool itself depended upon the downdrafts strengthened by sublimation of graupel higher in the storm. Finally, field campaign data revealed compelling inter-relationships between the cold-pool characteristics, and downdraft transience in terms of the microphysics. For example, pulses of hail and graupel (also internal factors) within downdrafts resulted in short-lived increases in the cold-pool depth. However, the cold-pool speed appeared to be relatively independent of the cold-pool depth and intensity, suggesting limitations of density current theory to convective cold pools in nature.

54 ENVIRONMENTAL SCIENCES↗

Dynamic systems modeling of the spallation neutron source cryogenic moderator system to optimize transient control and prepare for power upgrades

Through support of the US Department of Energy's Office of Basic Energy Sciences, Oak Ridge National Laboratory has begun applying machine learning methods to improve accelerator and target performance of the Spallation Neutron Source (SNS). These methods are being applied to the control optimization and power upgrade of the SNS Cryogenic Moderator System (CMS). A numerical model of the CMS has been developed to study these optimizations and system modifications using EcosimPro. This paper compares steady-state and transient numerical results with experimental data. Control optimization studies focused on dampening mass flow, temperature, and pressure fluctuations during sudden losses of accelerator beam power. This analysis was conducted by adjusting five proportional-integral-derivative controllers connected to four flow control valves and one heater. Future efforts include power uprate studies focused on increasing the CMS cooling capacity. The current accelerator power is 1.4 MW; the first target station is being upgraded to 2.0 MW as part of the Proton Power Upgrade effort. The CMS cooling capacity is sufficient for 2.0 MW operation.

47 OTHER INSTRUMENTATION↗

Model development for flash boiling spray and validations with isooctane, hexane, ethanol and their binary mixtures

The objective of this study is to develop a comprehensive model for spray atomization under the consideration of flash boiling, which includes nucleation, bubble growth, and droplet breakup. Although flash boiling has been previously studied, there is still no well-established analytical and numerical model of their mechanisms, or general guidelines on the implementation of this atomization strategy. In this study, detailed formulations are elaborated for the unified spray atomization model under both non-flash and flash boiling conditions. Next, the model was implemented into CFD software and 3-D simulations were performed for model validations. Two fuel blends are tested: ethanol-isooctane blends and hexane-isooctane blends. Simulation results were compared with the experimental measurements conducted using two methodologies, DBI (Diffuse Backlight Imaging) system for spray images captured by a high-speed camera, and a PDA (Phase Doppler Anemometry) system for droplet size measurement during the spray injection process. It was found that the current flash boiling model can accurately capture the spray characteristics under different environment conditions, such as the spray plume interactions and liquid penetration. The averaged Sauter mean diameter (SMD) calculations match the measurement sufficiently well. The spray structures for both single- and multi-component fuels during the injection process are well represented by current model.

42 ENGINEERING↗

Controlled 3D nanoparticle deposition by drying of colloidal suspension in designed thin micro-porous architectures

Nanoparticle deposition by drying of colloidal suspension in thin micro-porous architectures has attracted a lot of attention in scientific research as well as industrial applications. However, the underlying mechanisms of such three-dimensional (3D) deposition are not yet fully revealed due to the complexity of the co-occurring processes of two-phase fluid flows, phase change and mass transport. Consequently, the control of 3D nanoparticle deposition remains a challenge. We use a combined experimental and numerical approach to achieve controlled 3D nanoparticle deposition by drying of colloidal suspension in two pillar-based thin micro-porous architectures. By the design of pillar layout, rectangular-spiral and circular-spiral deposition configurations are obtained globally. By varying the surface wettability, vertically symmetric and sloped nanoparticle depositions can be achieved locally. While the numerical modeling reveals the mechanisms of liquid internal flow, as well as the impact of local drying rate on nanoparticle transport, accumulation and final deposition, the experimental results of deposition configurations validate the controlling strategies. This combined experimental and numerical work provides a framework to achieve desired 3D nanoparticle deposition in thin micro-porous architectures, with a thorough understanding of the underlying mechanisms of two-phase fluid flows, phase change and mass transport.

36 MATERIALS SCIENCE↗

Discrete modeling of a transformer with ALEGRA

We report progress on a task to model transformers in ALEGRA using the “Transient Magnetics” option. We specifically evaluate limits of the approach resolving individual coil wires. There are practical limits to the number of turns in a coil that can be numerically modeled, but calculated inductance can be scaled to the correct number of turns in a simple way. Our testing essentially confirmed this “turns scaling” hypothesis. We developed a conceptual transformer design, representative of practical designs of interest, and that focused our analysis. That design includes three coils wrapped around a rectangular ferromagnetic core. The secondary and tertiary coils have multiple layers. The tertiary has three layers of 13 turns each; the secondary has five layers of 44 turns; the primary has one layer of 20 turns. We validated the turns scaling of inductance for simple (one-layer) coils in air (no core) by comparison to available independent calculations for simple rectangular coils. These comparisons quantified the errors versus reduced number of turns modeled. For more than 3 turns, the errors are <5%. The magnetic field solver failed to converge (within 5000 iterations) for >10 turns. Including the core introduced some complications. It was necessary to capture the core surfaces in thin grid sheaths to minimize errors in computed magnetic energy. We do not yet have quantitative benchmarks with which to compare, but calculated results are qualitatively reasonable.

42 ENGINEERING↗

Loop Thermosyphon Enhanced Solar Collector

A two-phase loop thermosyphon solar collector was developed to transfer thermal energy efficiently and passively from a concentrated solar collector to a thermal desalination process. The two-phase flow in the loop thermosyphon is driven by a difference in gravitational pressure head between the liquid return line and the two-phase riser located between the evaporator and condenser. A numerical model incorporating mass, energy, and pressure drop balances was developed to predict system performance and use as a design tool. Three iterations of loop thermosyphon solar collectors were designed, fabricated, and tested to validate the numerical analysis and demonstrate system integration. Additionally, a novel solar receiver design was evaluated based on volumetric absorption of the concentrated sunlight by an optically absorptive two-phase working fluid. Life tests were performed to evaluate the working fluid thermal stability and resistance to ultraviolet degradation. A glass receiver tube and associated glass-to-metal junctions were developed. Finally, a full-scale loop thermosyphon was integrated with a parabolic trough solar collector and tested on-sun.

14 SOLAR ENERGY↗

Understanding Climate and Extreme Weather Events in the Greater New York Area

The goal of this project is to create a framework for modeling extreme events using climate information for risk assessment in the Greater New York City (NYC) region. The central research question we addressed was how best can we use information from historical climate data and numerical models of the climate system to estimate changing risks of rainfall extremes at specific locations – starting with the NYC region as a testbed.

54 ENVIRONMENTAL SCIENCES↗