SEARCH · Engineering Papers

Results for “EVAPORATION MODEL”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 631 records · Page 35

All sources NASA NTRS DOE OSTI

DOE OSTI2011

Calculation of the evaporation residue cross sections for the synthesis of the superheavy element Z=119 via the {sup 50}Ti+{sup 249}Bk hot fusion reaction

The evaporation residue (ER) cross sections for 3n and 4n evaporation channels in the {sup 50}Ti + {sup 249}Bk reaction leading the formation of {sup 296}119 and {sup 295}119 isotopes are evaluated by means of a modified fusion-by-diffusion model. In the model, the dynamic evolution from dinucleus to mononucleus is taken into account with the two-dimensional coupled Langenvin equations. The calculated maximum ER cross sections in 3n and 4n evaporation channels of the {sup 50}Ti + {sup 249}Bk reaction are 0.17 and 0.57 pb, respectively. The cross section of 0.57 pb is close to the present experimental limit for the registration of the evaporation residual nuclei. Therefore, superheavy element 119 may be the most hopeful new element with Z>118 to be synthesized under somehow improved experimental conditions in the near future.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2016

Heat Pump Clothes Dryer Model Development

A heat pump clothes dryer (HPCD) is an innovative appliance that uses a vapor compression system to dry clothes. Air circulates in a closed loop through the drum, so no vent is required. The condenser heats air to evaporate moisture out of the clothes, and the evaporator condenses water out of the air stream. As a result, the HPCD can achieve 50% energy savings compared to a conventional electric resistance dryer. We developed a physics-based, quasi-steady-state HPCD system model with detailed heat exchanger and compressor models. In a novel approach, we applied a heat and mass transfer effectiveness model to simulate the drying process of the clothes load in the drum. The system model is able to simulate the inherently transient HPCD drying process, to size components, and to reveal trends in key variables (e.g. compressor discharge temperature, power consumption, required drying time, etc.) The system model was calibrated using experimental data on a prototype HPCD. In the paper, the modeling method is introduced, and the model predictions are compared with experimental data measured on a prototype HPCD.

Heat Pump Clothes Dryer↗

DOE OSTI2023

Investigation of Surface and Marine-Cloud Coupling and its Impact on Cloud Droplet Number Concentrations and Cloud Cover Over the Southern Ocean

The proposed work involves characterizing and quantifying both the thermodynamic and dynamical coupling of marine low clouds (MLC) with the sea surface using Atmospheric Radiation Measurement (ARM) observations during MARCUS field campaign and an LES model. ARM data from multiple sensors will be used (e.g., Doppler cloud radar, ceilometer, and microwave radiometer) to characterize the MLC-surface coupling by virtue of the vertical structure and integrated quantities of boundary-layer clouds, aerosols, as well as atmospheric profiles and surface meteorology. We have used a combination of case studies and statistics-based composite analyses to find any linkages between large-scale dynamics, MLC-surface coupling, and cloud and boundary layer properties, the sequence of which reflects the chain of causality. An LES model with explicit aerosol physics, which includes the cycle of aerosols by being consumed as CCN and be regenerated following cloud droplet evaporation, has been run to determine the sources and sinks of Nd and their dependence on the degree of MLC-surface coupling. We have examined the systematic differences in both Nd and cloud occurrence between the two clusters under different meteorological conditions and further examine their respective roles, as well as causal relationships by means of LES modeling. This study helped improve our understanding of the ACI by differentiating the dynamic role of the coupling and cloud physics denoted by Nd, bridging the linkage in the chain toward understanding mechanisms governing the persistence of MLC over the Southern Ocean, solving the long-lasting problem of the cloud cover underestimation over the SO by GCMs. Ample ARM data and LES model have been employed to achieve the objectives of the study.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2014

Radiative and precipitation controls on root zone soil moisture spectra

Here, we present that temporal variability in root zone soil moisture content (w) exhibits a Lorentzian spectrum with memory dictated by a damping term when forced with white-noise precipitation. In the context of regional dimming, radiation and precipitation variability are needed to reproduce w trends prompting interest in how the w memory is altered by radiative forcing. A hierarchy of models that sequentially introduce the spectrum of precipitation, net radiation, and the effect of w on evaporative and drainage losses was used to analyze the spectrum of w at subtropical and temperate forested sites. Reproducing the w spectra at long time scales necessitated simultaneous precipitation and net radiation measurements depending on site conditions. The w memory inferred from observed w spectra was 25–38 days, larger than that determined from maximum wet evapotranspiration and field capacity. Finally, the w memory can be reasonably inferred from the Lorentzian spectrum when precipitation and evapotranspiration are in phase.

58 GEOSCIENCES↗

DOE OSTI2014

Region templates: Data representation and management for high-throughput image analysis

Compartment fires exhibit unique features associated with smoke accumulation and restricted air ventilation. The objective of the present study is to evaluate the ability of CFD-based fire models to simulate these features and in particular the effects of flame extinction that results from under-ventilated fire conditions. The study is performed using FireFOAM; FireFOAM is an advanced fire modeling software developed by FM Global and is based on a general-purpose open-source software called OpenFOAM. A new flame extinction model based on the concept of a critical value of the flame Damköhler number is incorporated into FireFOAM. The performance of the extinction model is evaluated via comparisons with a previously developed experimental database corresponding to a reduced-scale, heptane-fueled, compartment fire configuration. The numerical simulations also include a description of fuel evaporation driven by the computed gas-to-liquid heat feedback. Comparisons between experimental data and numerical results provide a suitable test bed to evaluate the ability of CFD-based fire models to describe the transition from over- to under-ventilated fire conditions, as well as the transition from extinction-free conditions to conditions in which the flame experiences partial or total quenching.

Energy & Fuels↗

DOE OSTI2019

Physics design for a lithium vapor box divertor experiment on Magnum-PSI

The lithium vapor box divertor is a potential solution for power exhaust in toroidal confinement devices. The divertor plasma interacts with a localized, dense cloud of lithium vapor, leading to volumetric radiation, cooling, recombination, and detachment. To minimize contamination of the core plasma, lithium vapor is condensed on cool (300°C to 400°C) baffles upstream of the detachment point. Before implementing this in a toroidal plasma device with a slot divertor geometry, we consider an experiment with a scaled baffled-pipe geometry in the high-power linear plasma device Magnum-PSI. Three 15 cm-scale open cylinders joined by 6 cm diameter ‘nozzles’ are positioned on the plasma beam axis upstream of a target. The central box may be loaded with several tens of grams of lithium, which can be evaporated at 650°C to produce a vapor predicted, using a simple plasma-neutral interaction model, to be dense enough to cause volumetric detachment in the plasma. The power delivered to the target and box walls as measured by increases in their temperatures after a 10 s plasma pulse can be compared to determine the effectiveness of the vapor in detaching the plasma. Direct Simulation Monte Carlo simulations are performed to estimate the flow rates of lithium vapor between the boxes and to estimate the trapping of H2 delivered by the plasma in the boxes, which could inadvertently lead to detachment. Details of the geometry, simulations, and possible diagnostic techniques are presented.

divertor↗

DOE OSTI2019

Thermodynamic versus Dynamic Controls on Extreme Precipitation in a Warming Climate from the Community Earth System Model Large Ensemble

The moisture budget is evaluated as a function of the probability distribution of precipitation for the end of the twentieth century and projected end of the twenty-first century in the Community Earth System Model Large Ensemble. For a given precipitation percentile, a conditional moisture budget equation relates precipitation minus evaporation ( P − E) to vertical moisture transport, horizontal moisture advection, and moisture storage. At high percentiles, moisture advection and moisture storage cancel and evaporation is negligible, so that precipitation is approximately equal to vertical moisture transport, and likewise for projected changes. Therefore, projected changes to extreme precipitation are approximately equal to the sum of thermodynamic and dynamic tendencies, representing changes to the vertical profiles of moisture content and mass convergence, respectively. The thermodynamic tendency is uniform across percentiles and regions as an intensification of the hydrological cycle, but the dynamic tendency is more complex. For extreme events, per degree of warming, in the mid-to-high latitudes the dynamic tendency is small, so that precipitation approximately scales by the Clausius–Clapeyron 7% K −1 increase. In the subtropics, a drying tendency originating from dynamics offsets the thermodynamic wetting tendency, with the net effect on precipitation varying among regions. The effect of this dynamic drying decreases with increasing percentile. In the deep tropics, a positive dynamic tendency occurs with magnitude similar to or greater than the positive thermodynamic tendency, resulting in generally a 10%–15% K −1 precipitation increase, and with a >25% K −1 increase over the tropical east Pacific. This reinforcing dynamical tendency increases rapidly for high percentiles.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Effect of particle shape on stratification in drying films of binary colloidal mixtures

The role of particle shape in evaporation-induced auto-stratification in polydisperse colloidal suspensions is explored with molecular dynamics simulations of mixtures of spheres and aspherical particles. A unified framework based on the competition between diffusion and diffusiophoresis is proposed to understand the effects of shape and size dispersity. In general, particles diffusing more slowly (e.g., larger particles) tend to accumulate more strongly at the evaporation front. However, larger particles have larger surface areas and therefore greater diffusiophoretic mobility. Hence, they are more likely to be driven away from the evaporation front via diffusiophoresis. For a rapidly dried bidisperse suspension containing small and large spheres, the competition leads to “small-on-top” stratification. Here, we employ a computational model in which the diffusion coefficient is inversely proportional to particle mass. For a mixture of spheres and aspherical particles with similar mass, the diffusion contrast is reduced, and the spheres are always enriched at the evaporation front as they have the smallest surface area for a given mass and, therefore, the lowest diffusiophoretic mobility. Furthermore, for a mixture of solid and hollow spheres that have the same outer radius and thus the same surface area, the diffusiophoretic contrast is suppressed, and the system is dominated by diffusion. Consequently, the solid spheres, which have a larger mass and diffuse more slowly, accumulate on top of the hollow spheres. Finally, for a mixture of thin disks and long rods that differ significantly in shape but have similar mass and surface area, both diffusion and diffusiophoresis contrasts are suppressed, and the mixture does not stratify.

Classical molecular dynamic simulations↗

DOE OSTI2012

Dynamic Modeling Strategy for Flow Regime Transition in Gas-Liquid Two-Phase Flows

In modeling gas-liquid two-phase flows, the concept of flow regimes has been widely used to characterize the global interfacial structure of the flows. Nearly all constitutive relations that provide closures to the interfacial transfers in two-phase flow models, such as the two-fluid model, are flow regime dependent. Current nuclear reactor safety analysis codes, such as RELAP5, classify flow regimes using flow regime maps or transition criteria that were developed for steady-state, fully-developed flows. As twophase flows are dynamic in nature, it is important to model the flow regime transitions dynamically to more accurately predict the two-phase flows. The present work aims to develop a dynamic modeling strategy to determine flow regimes in gas-liquid two-phase flows through introduction of interfacial area transport equations (IATEs) within the framework of a two-fluid model. The IATE is a transport equation that models the interfacial area concentration by considering the creation of the interfacial area, fluid particle (bubble or liquid droplet) disintegration, boiling and evaporation, and the destruction of the interfacial area, fluid particle coalescence and condensation. For flow regimes beyond bubbly flows, a two-group IATE has been proposed, in which bubbles are divided into two groups based on their size and shapes, namely group-1 and group-2 bubbles. A preliminary approach to dynamically identify the flow regimes is discussed, in which discriminator s are based on the predicted information, such as the void fraction and interfacial area concentration. The flow regime predicted with this method shows good agreement with the experimental observations.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2016

Dissociation of doubly charged clusters of lithium acetate: Asymmetric fission and breakdown of the liquid drop model: Dissociation of doubly charged clusters of lithium acetate

Unimolecular and collision-induced dissociation of doubly charged lithium acetate clusters, (CH3COOLi)nLi22+, demonstrated that Coulomb fission via charge separation is the dominant dissociation process with no contribution from the neutral evaporation processes for all such ions from the critical limit to larger cluster ions, although latter process have normally been observed in all earlier studies. These results are clearly in disagreement with the Rayleigh’s liquid drop model that has been used successfully to predict the critical size and explain the fragmentation behavior of multiply charged clusters.

Cluster Ions↗

DOE OSTI2017

The influence of land cover on surface energy partitioning and evaporative fraction regimes in the U.S. Southern Great Plains

Land-atmosphere interactions are important to climate prediction, but the underlying effects of surface forcing of the atmosphere are not well understood. In the U.S. Southern Great Plains, grassland/pasture and winter wheat are the dominant land covers but have distinct growing periods that may differently influence land-atmosphere coupling during spring and summer. Variables that influence surface flux partitioning can change seasonally, depending on the state of local vegetation. We use surface observations from multiple sites in the U.S. Department of Energy Atmospheric Radiation Measurement Southern Great Plains Climate Research Facility and statistical modeling at a paired grassland/agricultural site within this facility to quantify land cover influence on surface energy balance and variables controlling evaporative fraction (latent heat flux normalized by the sum of sensible and latent heat fluxes). We demonstrate that the radiative balance and evaporative fraction are closely related to green leaf area at both winter wheat and grassland/pasture sites and that the early summer harvest of winter wheat abruptly shifts the relationship between evaporative fraction and surface state variables. Prior to harvest, evaporative fraction of winter wheat is strongly influenced by leaf area and soil-atmosphere temperature differences. After harvest, variations in soil moisture have a stronger effect on evaporative fraction. Finally, this is in contrast with grassland/pasture sites, where variation in green leaf area has a large influence on evaporative fraction throughout spring and summer, and changes in soil-atmosphere temperature difference and soil moisture are of relatively minor importance.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2011

Thickness dependent self limiting 1-D tin oxide nanowire arrays by nanosecond pulsed laser irradiation

Fast, sensitive and discriminating detection of hydrogen at room temperature is crucial for storage, transportation, and distribution of hydrogen as an energy source. One dimensional nanowires of SnO2 are potential candidates for improved H2 sensor performance. The single directional conducting continuous nanowires can decrease electrical noise, and their large active surface area could improve the response and recovery time of the sensor. In this work we discuss synthesis and characterization of nanowire arrays made using nanosecond ultraviolet wavelength (266 nm) laser interference processing of ultrathin SnO2 films on SiO2 substrates. The laser energy was chosen to be above the melting point of the films. The results show that the final nanowire formation is dominated by preferential evaporation as compared to thermocapillary flow. The nanowire height (and hence wire aspect ratio) increased with increasing initial film thickness ho and with increasing laser energy density Eo. Furthermore, a self-limiting effect was observed where-in the wire formation ceased at a specific final remaining thickness of SnO2 that was almost independent of ho for a given Eo. To understand these effects, finite element modeling of the nanoscale laser heating was performed. This showed that the temperature rise under laser heating was a strong non-monotonic function of film thickness. As a result, the preferential evaporation rate varies as wire formation occurs, eventually leading to a shut-off of evaporation at a characteristic thickness. This results in the stoppage of wire formation. This combination of nanosecond pulsed laser experiments and thermal modeling shows that several unique synthesis approaches can be utilized to control the nanowire characteristics.

08 HYDROGEN↗

DOE OSTI2016

Vapor shielding models and the energy absorbed by divertor targets during transient events

The erosion of divertor targets caused by high heat fluxes during transients is a serious threat to ITER operation, as it is going to be the main factor determining the divertor lifetime. Under the influence of extreme heat fluxes, the surface temperature of plasma facing components can reach some certain threshold, leading to an onset of intense material evaporation. The latter results in formation of cold dense vapor and secondary plasma cloud. This layer effectively absorbs the energy of the incident plasma flow, turning it into its own kinetic and internal energy and radiating it. This so called vapor shielding is a phenomenon that may help mitigating the erosion during transient events. In particular, the vapor shielding results in saturation of energy (per unit surface area) accumulated by the target during single pulse of heat load at some level E{sub max}. Matching this value is one of the possible tests to verify complicated numerical codes, developed to calculate the erosion rate during abnormal events in tokamaks. The paper presents three very different models of vapor shielding, demonstrating that E{sub max} depends strongly on the heat pulse duration, thermodynamic properties, and evaporation energy of the irradiated target material. While its dependence on the other shielding details such as radiation capabilities of material and dynamics of the vapor cloud is logarithmically weak. The reason for this is a strong (exponential) dependence of the target material evaporation rate, and therefore the “strength” of vapor shield on the target surface temperature. As a result, the influence of the vapor shielding phenomena details, such as radiation transport in the vapor cloud and evaporated material dynamics, on the E{sub max} is virtually completely masked by the strong dependence of the evaporation rate on the target surface temperature. However, the very same details define the amount of evaporated particles, needed to provide an effective shielding to the target, and, therefore, strongly influence resulting erosion rate. Thus, E{sub max} cannot be used for validation of shielding models and codes, aimed at the target material erosion calculations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Experiments on a vapor compression air conditioner with liquid desiccants for efficient dehumidification

Buildings require air conditioning systems that not only cool and dehumidify supply air but also provide sufficient ventilation to ensure indoor air quality and occupant comfort. However, standard recirculation systems-which introduce about a 10 % to 20 % fraction of outdoor air-often fail to deliver air that is precisely cooled and dry, particularly because 80-90 % of the ventilation cooling load is latent. Mixing humid ventilation air with recirculated indoor air increases the energy and costs required to condition the air to comfortable levels. Dedicated outdoor air systems (DOASs) are designed to handle this latent dominated ventilation load and thus need to have efficient humidity removal. Many cooling cycles can perform this task. Here we describe a liquid desiccant DOAS, which combines a vapor compression cycle and a liquid desiccant absorber and desorber pair. We present its performance at 26 operating conditions and a thermodynamic model which can accurately predict the moisture removal efficiency. The model's performance predictions have a mean percentage error of 2.5 % and a coefficient of variation of the root mean square error of 7.5 %. We also compare the performance of this vapor-compression-coupled liquid desiccant system with a standard vapor compression system with the same components but no liquid desiccant. For the 26 conditions tested in this study, this comparison shows that adding liquid desiccants lowers the required evaporator cooling load by 21 %, allows for 25 % lower compressor volumetric capacity, and 25 % lower electricity use. Future work will leverage this model to quantify the reduction in annual electricity use across different climates, including the need for a standard vapor compression system to reheat the air during some of the year.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2017

Charging and Release Mechanisms of Flexible Macromolecules in Droplets

We study systematically the charging and release mechanisms of a flexible macromolecule, modeled by poly(ethylene glycol) (PEG), in a droplet by using molecular dynamics simulations. We compare how PEG is solvated and charged by sodium Na{sup +} ions in a droplet of water (H{sub 2}O), acetonitrile (MeCN), and their mixtures. Initially, we examine the location and the conformation of the macromolecule in a droplet bearing no net charge. It is revealed that the presence of charge carriers do not affect the location of PEG in aqueous and MeCN droplets compared with that in the neutral droplets, but the location of the macromolecule and the droplet size do affect the PEG conformation. PEG is charged on the surface of a sodiated aqueous droplet that is found close to the Rayleigh limit. Its charging is coupled to the extrusion mechanism, where PEG segments leave the droplet once they coordinate a Na{sup +} ion or in a correlated motion with Na{sup +} ions. In contrast, as PEG resides in the interior of a MeCN droplet, it is sodiated inside the droplet. The compact macro-ion transitions through partially unwound states to an extended conformation, a process occurring during the final stage of desolvation and in the presence of only a handful of MeCN molecules. For charged H{sub 2}O/MeCN droplets, the sodiation of PEG is determined by the H{sub 2}O component, reflecting its slower evaporation and preference over MeCN for solvating Na{sup +} ions. We use the simulation data to construct an analytical model that suggests that the droplet surface electric field may play a role in the macro-ion–droplet interactions that lead to the extrusion of the macro-ion. This study provides the first evidence of the effect of the surface electric field by using atomistic simulations. .

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2017

Fast and slow shifts of the zonal-mean intertropical convergence zone in response to an idealized anthropogenic aerosol

Previous modeling work showed that aerosol can affect the position of the tropical rain belt, i.e., the intertropical convergence zone (ITCZ). Yet it remains unclear which aspects of the aerosol impact are robust across models, and which are not. Here we present simulations with seven comprehensive atmosphere models that study the fast and slow impacts of an idealized anthropogenic aerosol on the zonal-mean ITCZ position. The fast impact, which results from aerosol atmospheric heating and land cooling before sea-surface temperature (SST) has time to respond, causes a northward ITCZ shift. Yet the fast impact is compensated locally by decreased evaporation over the ocean, and a clear northward shift is only found for an unrealistically large aerosol forcing. The local compensation implies that while models differ in atmospheric aerosol heating, this does not contribute to model differences in the ITCZ shift. The slow impact includes the aerosol impact on the ocean surface energy balance and is mediated by SST changes. The slow impact is an order of magnitude more effective than the fast impact and causes a clear southward ITCZ shift for realistic aerosol forcing. Models agree well on the slow ITCZ shift when perturbed with the same SST pattern. However, an energetic analysis suggests that the slow ITCZ shifts would be substantially more model-dependent in interactive-SST setups due to model differences in clear-sky radiative transfer and clouds. In conclusion, we also discuss implications for the representation of aerosol in climate models and attributions of recent observed ITCZ shifts to aerosol.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Quantitative Analysis of the Effects of Natural and Human Factors on a Hydrological System in Zhangweinan Canal Basin

Drought is a water deficit state caused by large-scale climate change that cannot be avoided by a water resource management system. Water scarcity refers to the unsustainable utilization of water resources over a long time, which is the result of water management policies. However, water deficits caused by drought and water scarcity often occur simultaneously and are indistinguishable. This study proposes a model-based simulation framework that can quantitatively distinguish natural factors (drought) from human factors (water scarcity) in a hydrological system. The simulation was applied to the Zhangweinan Canal Basin, based on the runoff sequences from 1950 to 2004. The results show that the runoff curve number, soil depth, soil available water, soil evaporation compensation coefficient, base runoff α coefficient and the maximum canopy interception have the highest sensitivity to runoff, and that the calibrated and validated SWAT model can effectively simulate the runoff process in the Zhangweinan Canal Basin and similar areas. Abrupt changes in human activities in 1975 and water scarcity led to the disappearance of the summer peak runoff period in both wet and dry years. Human factors are the main reason for the change in the hydrological system in the study area; the runoff loss caused by human factors is four times that caused by natural factors according to the proposed variable threshold. This study proposes a model-based simulation framework that can help water resource managers to distinguish the effects of drought and water scarcity in water-stressed areas and adjust management accordingly.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Synthesis of observed and simulated rain microphysics to inform a new Bayesian statistical framework for microphysical parameterization in climate models (Final Scientific Report)

This project aimed to investigate a new approach to bulk microphysics schemes. We identified the need to move beyond the fixed structural assumptions and approximation of existing schemes. For example, most schemes assume some functional form for the rain drop size distribution (e.g. a gamma or exponential distribution). Most schemes then evolve some number of statistical moments of that distribution via various processes, such as evaporation, sedimentation, collision-coalescence, and collisional breakup. These microphysical processes, in turn, are typically some fixed functional form derived from either some other (more detailed) model, or using some single-particle rates that are then integrated over the size distribution. While these process rate functions may have some free parameters to adjust, in other cases doing so is impossible (e.g. when the process rates are analytical or piecewise solutions to some target function). Our goal was to build and test a scheme that assumed no size distribution form, and used a series of power laws as the basis for the process rates. Moments of the size distribution would be predicted, but no underlying size distribution would be specified. Any number or choice of size distribution moments could be used, and any number of power laws could be employed to model the process rates. Thus, our approach is flexible, and can seamlessly scale across levels of complexity, as demanded by the data. Bayesian inference then would provide the formalism to estimate the model parameters and structure, allowing for robust uncertainty quantification. We proposed testing this framework in idealized simulations using bin microphysical schemes as a data source. We also proposed using real data to inform our microphysics scheme, and also that we would integrate our scheme into WRF.

58 GEOSCIENCES↗