SEARCH · Engineering Papers

Results for “water 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 OSTI2025

FAST Irradiations, Postirradiation Examinations, and Modeling of U-Mo for Light Water Reactor Applications

Many next generation light water reactor (LWR) concepts, such as mobile small modular reactors, are seeking to use smaller core dimensions than conventional reactor types. Smaller reactor cores require an increase in fissile material to maintain reactivity. For non-proliferation purposes, enrichment increases are limited to less than 20% (high assay low enriched uranium, [HALEU]) and so higher uranium density fuels than UO 2 must be considered. To this end, uranium-molybdenum alloys were tested using the Fission Accelerated Steady-state Test (FAST) approach. The experiment test matrix is focused on identifying the temperature transition between low swelling and high fission gas retention to break away swelling and low fission gas retention. This paper documents the results of irradiation tests and post-irradiation examinations (PIE) including neutron radiography, rodlet profilometry, fission gas collection analysis, and optical metallography. The results of these tests showed that unconstrained U-Mo fuels (solid, Na-bonded rodlets) have a swelling threshold between 400-450°C with minimal fission gas release below this point. Higher temperature solid fuel showed microstructural zoning with small pore networks while lower temperature solid fuels have a uniform microstructure with large pore networks. U-Annular Mo fuels where swelling had some self-constraint imposed upon it, were shown to have much reduced swelling compared to their solid counterparts as well as very low fission gas release for irradiation temperatures up to 500°C. These initial results show that the use of U-Mo in constrained fuel geometries could be used as a high uranium density HALEU fuel for LWRs.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2023

Modeling Assessment of Residential Air-to-Water Heat Pumps Coupled with Cooling Thermal Storage

This study explored the performance and operating cost viability of air-to-water heat pumps (AWHPs) coupled with thermal energy storage (TES) in efficient new residential construction.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2023

Advanced Modeling and Simulation to Characterize Advanced Boiling Water Reactor Source Terms to Support a Regulatory Approval Pathway for Right-Sized Emergency Planning Zone. Advanced Boiling Water Reactor (A-BWR) Project Phase 2

The United States Department of Energy (US DOE) is currently supporting the development of various small modular reactor (SMR), microreactor and advanced reactor designs. These reactors have improved safety features as compared to conventional Large Light Water Reactors (LLWRs), which include features that improve potential reduction of radiological source terms in the event of design and beyond-design basis accidents. Specifically, some reactors feature a relatively smaller containment volume with respect to the available surface area for fission product deposition; these include Integrated Pressurized Water Reactors (iPWRs), such as the NuScale SMR and SMR-type Boiling Water Reactors (BWRs) such as General Electric-Hitachi’s BWRX-300 design.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2022

Modeling and Analysing the Impact of Heat Pump Water Heaters on Distribution Systems Using GridLAB-D

With the constant increase in energy demand, finding ways to reduce peak load and the energy-costs factors has become more imperative. Domestic water heating showcases a significant opportunity for such applications. Water heating is the second-highest energy consumer in the residential sector across the United States. Electric Water Heaters (EWHs), in particular, constitute nearly 43% of American household water heating energy consumption. Heat Pump Water Heater (HPWH), on the other hand, are an advanced water heating technology that has recently emerged in the United States residential market. The objectives of this work are to develop a HPWH model and build a case study that evaluates various penetration levels of HPWH in providing reduced peak load and cost-effective energy savings for both utilities and customers. The HPWH model was developed and integrated within the GridLAB-D simulation environment. The model behavior was then validated against a real HPWH unit at Portland State University (PSU). The case studies incorporated five HPWH penetration levels, ranging from 20% to 100%. In each case, EWHs were replaced with HPWHs. The results showed that a high penetration level of HPWHs can reduce the energy consumption on a distribution system to 38%.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2021

Condensed Phase Water Molecular Multipole Moments from Deep Neural Network Models Trained on Ab Initio Simulation Data

Ionic solvation phenomena in liquids involve intense interactions in the inner solvation shell. For interactions beyond the first shell, the ion–solvent interaction energies result from the sum of many smaller-magnitude contributions that can still include polarization effects. Deep neural network (DNN) methods have recently found wide application in developing efficient molecular models that maintain near-quantum accuracy. In this work, we extend the DeePMD-kit code to produce accurate molecular multipole moments in the bulk and near interfaces. The new method is validated by comparing the DNN moments with those generated by ab initio simulations. The moments are used to compute the electrostatic potential at the center of a molecular-sized hydrophobic cavity in water. The results show that the fields produced by the DNN models are in quantitative agreement with the AIMD-derived values. These efficient methods will open the door to more accurate solvation models for large solutes such as proteins.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Surface Salinity Under Transitioning Ice Cover in the Canada Basin: Climate Model Biases Linked to Vertical Distribution of Fresh Water

Abstract The Canada Basin has exhibited a significant trend toward a fresher surface layer and thus a more stratified upper‐ocean over the past three decades. State‐of‐the‐art ice‐ocean models, by contrast, tend to simulate a surface layer that is saltier and less stratified than observed. Here, we examine decadal changes to seasonal processes that may contribute to this wide‐reaching model bias using climate model simulations from the Community Earth System Model and below‐ice observations from the Arctic Ice Dynamics Joint Experiment in 1975 and Ice Tethered Profilers in 2006–2012. In contrast to the observations, the models simulate salinity profiles that show relatively little variation between 1975 and 2012. We demonstrate that this bias can be mainly attributed to unrealistically deep vertical mixing in the model, creating a surface layer that is saltier than observed. The results provide insight for climate model improvement with broad implications for Arctic sea ice and ecosystem dynamics.

58 GEOSCIENCES↗

DOE OSTI2024

A More Transparent Infrared Window

The infrared window region (780–1,250 cm –1 , 12.8 to 8.0 μm) is of great importance to Earth's climate due to its high transparency and thermal energy. We present here a new investigation of the transparency of this spectral region based on observations by interferometers of downwelling surface radiance at two DOE Atmospheric Radiation Measurement program sites. We focus on the dominant source of absorption in this region, the water vapor continuum, and derive updated values of spectral absorption coefficients for both the self and foreign continua. Our results show that the self continuum is too strong in the previous version of Mlawer-Tobin_Clough-Kneizys-Davies (MT_CKD) water vapor continuum model, a result that is consistent with other recent analyses, while the foreign continuum is too weak in MT_CKD. In general, the weaker self continuum derived in this study results in an overall increase in atmospheric transparency in the window, although in atmospheres with low amounts of water vapor the transparency may slightly decrease due to the increase in foreign continuum absorption. These continuum changes lead to a significant decrease in downwelling longwave flux at the surface for moist atmospheres and a modest increase in outgoing longwave radiation. The increased fraction of surface-leaving radiation that escapes to space leads to a notable increase (~5–10%) in climate feedback, implying that climate simulations that use the new infrared window continuum will show somewhat less warming than before. This study also points out the possibly important role that aerosol absorption may play in the longwave radiative budget.

58 GEOSCIENCES↗

DOE OSTI2025

Approaches for Water Removal in Direct-Fired sCO 2 Power Cycles

There is interest in investigation of water removal processes in direct fired sCO 2 flows, as this may potentially lead to greater system efficiency as the removal of this contaminant will result in sCO 2 behaving close to idealized behaviors. Water removal should be split into a two-step process, condensation of the water, followed by separation of the liquid phase water from the sCO 2 . The two main avenues of condensation are manipulation of pressure and temperature for phase change. For this paper, temperature-based phase change is the primary focus through the implementation of heat exchangers. Of the heat exchangers investigated it was found that printed circuit heat exchangers (PCHEs) could be an alternative for this use case, though the specific design of flow channel geometry and flow direction depends on the specific system case and cannot be determined at this point. For water separation there were four processes identified, all of which already assume water is in liquid phase at that point in the system. Of these separation avenues the best candidate is the hydrocyclone as it has a proven history of separating liquid-liquid phase mixtures with small density differences in oilfield use, in addition they have been investigated and modeled specifically for water separation for sCO 2 flows and the footprint is relatively small.

20 FOSSIL-FUELED POWER PLANTS↗

DOE OSTI2024

Data/code for Mlawer et al. 2024 (A more transparent infrared window)

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Progress Towards the Validation of a new RELAP5-3D model of the High Temperature Test Facility

Validation is a key step in the development of any type of systems model. As the next generation of reactors approaches, the need for codes that have been validated for these new types of systems continues to grow. An example of a prominent option is the Reactor Excursion Leak Analysis Program (RELAP5-3D), developed by Idaho National Laboratory. This code was developed for the purpose of systems level thermal-hydraulic modeling of light water reactors (LWRS) and postulated transients that can occur in LWRS.RELAP5-3D has been substantially validated against LWR data. Due to its long history as a reactor safety analysis tool, there has been an effort to adapt RELAP5-3D for the purposes of advanced reactor concepts such as prismatic high-temperature gas-cooled reactors (HTGRs). However, RELAP5-3D has not nearly been validated and verified for HTGRs to the degree of LWRs, warranting verification and validation opportunities with computational benchmarks and existing experimental facilities. Examples of such facilities include the modular high-temperature gas-cooled reactor (MHTGR) 350 and the high temperature engineering test reactor (HTTR) from Japan. The MHTGR 350 is a benchmark design concept for code-to-code verification purposes; therefore, it does not provide any experimental data for validation opportunities The HTTR provides useful multiphysics validation data but does not have the in-core instruments to generate thermal-hydraulic experimental data to help with RELAP5-3D validation. Consequently, a facility that could provide key in-core temperatures for thermal-hydraulic validation was still needed. The High Temperature Test Facility (HTTF) is an integral effects facility for HTGR thermal hydraulics developed and operated by Oregon State University. HTTF represents ¼ length scale of the General Atomics MHTGR and is rated for a total power of 2.2 MW. Axially, the core consists of an upper and lower reflector and 10 blocks, numbered from bottom to top (Block 1 is right above lower reflector). The core is heated via graphite resistive heater rods, with respective channels distributed throughout the core. The primary coolant is helium and heat can radiate out of the core to the reactor cavity cooling system (RCCS), which is cooled by water. The primary purpose of the facility is to investigate pressurized conduction cooldown (PCC) and depressurized conduction cooldown (DCC) transients, which are also referred to as the pressurized and depressurized loss of forced cooling respectively. Two experiments were chosen to perform the validation study with a RELAP5-3D model of HTTF. These experiments are PG-27 (PCC) and PG-29 (DCC). These were chosen based off of the quality of available experimental data before and during the experiment which led to their inclusion in the HTGR Thermal Hydraulics Benchmark.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2026

Knowledge-guided graph machine learning for spatially distributed prediction of daily discharge and nitrogen export dynamics

Spatially distributed prediction of streamflow and nitrogen export dynamics is essential for precision management of agricultural watersheds. While temporal deep learning models such as Long Short-Term Memory (LSTM) have shown strong performance at basin scales, their ability to generalize spatially is limited by insufficient representation of spatial dependencies and flow paths, particularly under data-scarce conditions. To address this gap, we propose HydroGraphNet, a knowledge-guided graph machine learning framework that integrates process-based knowledge and explicit spatial learning into temporal modeling. This framework incorporates directed graph topology to encode watershed connectivity and upstream inflows, with mass balance constraints to improve physical consistency. To enhance generalization in sparsely monitored regions, HydroGraphNet is pretrained on synthetic data generated by the SWAT+ (Soil and Water Assessment Tool Plus) model. We evaluated HydroGraphNet in the Upper Sangamon River Basin (44 HUC-12 subwatersheds, 2001–2020) against two LSTM baselines: a lumped basin-level model and a distributed variant. When benchmarked on SWAT+ simulations in pretraining, HydroGraphNet improved test NSEs by 8.9% (discharge) and 13.7% (NO₃–N load) in temporal extrapolation, and by 27.1% and 34.7% in spatial extrapolation, relative to the Lumped LSTM baseline. After fine-tuning with USGS monitoring data, the model achieved mean test NSE (KGE) scores of 0.768 (0.861) for discharge and 0.626 (0.664) for NO₃–N load, substantially outperforming baselines. Attribution analysis further highlighted the importance of upstream inflow representation and graph-based spatial learning in capturing cross-subwatershed dependencies. The model also reproduced seasonal hydrological and biogeochemical patterns consistent with known processes, demonstrating its robustness and process fidelity for spatially distributed prediction. Altogether, HydroGraphNet advances the integration of physical knowledge and spatially explicit learning in hydrological modeling, offering a generalizable framework for distributed modeling to support spatially targeted water quality management in data-scarce watersheds.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

How Well Do U.S. Western Water Markets Convey Economic Information?

An efficient market implies that potential gains from trade are fully captured. Achieving this requires a well-functioning market where prices reflect all available information. In the case of water rights markets, this implies that the permanent water rights transfer price reflects the sum of discounted returns to this asset (i.e., the lease price), the market interest rate, and a risk premium that reflects potential future water scarcity. Here, the purpose of this study is to assess the efficiency of western U.S. water markets by using the asset pricing model to measure how well prices reflect long-run returns to permanent water rights.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

DOE OSTI2026

Improving coastal water level estimation by merging nadir-only satellite altimetry data into a hydrodynamic model

Providing robust real time flood warnings is of paramount importance to coastal communities. Although state-of-the-art hydrodynamic models are capable of robustly predicting Coastal Water Levels (CWL), unresolved drivers affecting level fluctuations are often not represented by the model governing equations. This work evaluates a novel method to improve the performance of the ADvanced CIRCulation (ADCIRC) hydrodynamic model by assimilating observations from four nadir-only satellite altimetry missions against a set of National Oceanic and Atmospheric Administration (NOAA) gauge stations located across the entire U.S. East Coast. Two different types of simulations were performed – Open Loop (OL) and Data Assimilation (DA). Five different simulations were performed where four different satellite altimetry observations were assimilated individually and combined with two different scenarios – with and without considering the data quality flags. Results indicate that, despite their limited spatial coverage, merging nadir-only observations into ADCIRC from the newly launched Surface Water and Ocean Topography (SWOT)’s nadir altimeter can improve the model performance at 76% of the gauge locations, whereas Sentinel-6 improves 73% of the total locations, Jason-3 74%, and SARAL 21%. Furthermore, combining observations from SWOT-nadir, Jason-3, and Sentinel-6 can improve the ADCIRC performance at more than 80% of the gauge locations for 107-day simulation. Nadir-only satellite altimetry observations can be useful for improving the model performance even if flagged as “poor quality” near the coast. When the flagged data are disregarded, SWOT can improve ADCIRC at 78%, Sentinel-6 at 73%, Jason-3 at 53%, and SARAL at 21% of the gauge locations. The ability to improve the model simulations largely depends on the availability of a satellite overpass nearby. Therefore, model performance can be further enhanced if satellite observations are available during a storm surge event, stressing the importance of frequent satellite overpasses.

Aafnan Bhuiyan, Soelem↗

DOE OSTI2026

Accelerated Carbon and Water Cycles in the Amazon and Congo Basins Revealed From TRENDY Models and Remote Sensing Products

Tropical forests play a vital role in the global carbon cycle and land–atmosphere interactions. Estimating tropical forest carbon–water dynamics is challenging due to observational and modeling uncertainties. This study leverages the “Trends and drivers of the regional scale terrestrial sources and sinks of carbon dioxide” (TRENDY) project models and satellite observations to assess changes (2003–2021) in vegetation carbon, gross primary production (GPP), evapotranspiration (ET), and net biosphere production (NBP) in the Amazon and Congo. Atmospheric CO 2 , climate variability, and land use and land cover changes constrain these variables between 1700 and 2021 with the overall increasing trends of carbon stock and fluxes. The models overestimate vegetation carbon and GPP, while ET and NBP are consistent with observations. Fire-activated models predict lower values for vegetation carbon and GPP, ET, and NBP, aligning more closely with observations. The higher ET from fire-activated models may result from enhanced soil evaporation due to increased canopy openings. Fire-inactivated models could well estimate the magnitudes of NBP. The high vegetation carbon in nitrogen-enabled models points to simulation uncertainties and imbalance in model numbers regarding the nitrogen cycle. Although the nitrogen cycle enhances water use efficiency in both the Amazon and Congo, the models show a higher sensitivity to the nitrogen cycle in the Congo. This study highlights the challenges in accurately representing tropical biogeochemical cycles and the values of satellite products in model evaluations, underscoring the need for standard modeling protocols that address biogeochemical components (e.g., nutrient cycles) to better resolve process-based representations.

Shi, Mingjie [Pacific Northwest National Laborator↗

DOE OSTI2020

FY20 Verification of BISON Using Analytic and Manufactured Solutions

In 2010, the U.S. Department of Energy created its first Energy Innovation Hub, which is focused on developing high-fidelity and high-resolution Modeling and Simulation (M&S) tools for modeling of Light Water Reactors (LWRs). This hub, Consortium for Advanced Simulation of LWRs (CASL), has developed an LWR simulation tool called Virtual Environment for Reactor Applications (VERA). The multi-physics capability of VERA is achieved through the coupling of single-physics codes, including BISON, CTF, MPACT, and MAMBA. BISON is a fuel performance code which models the thermo-mechanical behavior of nuclear fuel using high performance M&S. It is capable of modeling traditional LWR fuel rods, fuel plates, and TRi-structural ISOtropic (TRISO) fuel particles. It can employ three-dimensional Cartesian, two-dimensional axisymmetric cylindrical, or one-dimensional radial spherical geometry. It includes empirical models for a large variety of fuel physics: temperature- and burnup-dependent thermal properties, fuel swelling and densification, fission gas production, cladding creep, fracture, cladding plasticity, and gap/plenum models. This document details a series of code verification test problems that are used to test BISON. These problems add confidence that the BISON code is a faithful representation of its underlying mathematical model. The suite of verification tests are mapped to the underlying conservation equations solved by the code: heat conduction, mechanics, and species conservation. Twenty-two problems are added for the heat conduction solution, two for the mechanics solution, and none for species conservation. Method of Manufactured Solutions (MMS) capability is demonstrated with three problems, and temperature drops across the fuel gap are tested.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2022

Methods for Computing Physically Realistic Estimates of Electric Water Heater Demand Response Resource Suitable for Bulk Power System Planning Models

Demand response is commonly called on to reduce load during system peak times or to respond to contingency events. In future power systems with higher shares of wind and solar generation (which we describe together as variable generation [VG]), demand response could have more opportunities to provide energy shifting or operating reserve services. This report evaluates the ability of residential electric water heaters, both electric resistance water heaters (ERWHs) and heat pump water heaters (HPWHs), to provide such services starting from detailed whole-building energy models that realistically represent New England single family home stock. We use a parsimonious surrogate model to represent operational flexibility in a form suitable for linear and mixed integer programming. This enables relatively fast determination of aggregate contingency reserve resource, price-taking energy shifting outcomes, and in some cases the determination of aggregate models at the megawatt (MW) scale that can be directly included in large-scale grid models. After selecting modeling methods and parameters through various computational experiments, we find interquartile ranges of contingency reserve resource in ISO-NE for about 603,400 ERWHs of 45 MW - 69 MW for Claim10 (50 minute responses provided with 10 minutes of advanced notification) and 65 MW - 102 MW for Claim30 (30 minute responses provided with 30 minutes of advanced notification), and for about 619,000 HPWHs of 48 MW - 88 MW for Claim10 and 52 MW - 90 MW for Claim30. The overall reserve resource is up to 32% of total load for ERWHs providing Claim10 service, 47% for ERWHs providing Claim30 service, 93% for HPWHs providing Claim10 service, and 97% for HPWHs providing Claim30 service. More work is required to determine if HPWHs are inherently more suitable than ERWHs for providing contingency reserve or if these results reflect idiosyncrasies of the single family home stock model used in this study. The value of this contingency resource in a Near-term VG model of ISO-NE is $\$ 0.40$ to $\$1.20$ per water heater-year, and significantly larger, $\$ 3.80$ to $\$ 5.30$ per water heater-year in a Mid-term VG model of ISONE. Aggregating surrogate models to the MW-scale for energy shifting service is more challenging than for contingency service and we only present such results for ERWHs, because we were unable to determine satisfactory ways to deal with HPWHs' time-varying and path dependent operational characteristics. Individual surrogate models suitable for evaluating the energy shifting resource from both ERWHs and HPWHs are created, however, and dispatched against day-ahead prices from the Near-Term VG and Mid-Term VG models of ISO-NE. The individual surrogate models are able to access and potentially shift all 640 GWh of HPWH load and 1,547 GWh of ERWH load we modeled in two different single family home stock models. In contrast, the most effective model of aggregate ERWH shifting resource we created only captured 34.7% of the total ERWH load. Energy shifting affected by price-taking dispatch against modeled day-ahead energy prices produces per water heater year profits of $\$19.44$ - $\$22.93$ for individual HPWHs, $\$39.11$ - $\$40.54$ for individual ERWHs, and up to $\$4.00$ - $\$4.24$ for aggregated ERWHs, with the variations mainly due to grid conditions (more or less VG). When the supply-side response to these changes is accounted for, the per water heater year production cost savings for ISO-NE are $\$7.50$ to $\$17.70$ for the most effective set of endogenously dispatched aggregate ERWHs, $\$15.60$ to $\$15.70$ for individual ERWHs dispatched against the DA prices, and $\$10.70$ to $\$11.20$ for individual HPWHs dispatched against DA prices. Those ranges primarily represent the difference between Near-Term VG and Mid-Term VG grid conditions.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2024

Asymmetric Bipolar Membrane for High Current Density Electrodialysis Operation with Exceptional Stability

Bipolar membranes (BPMs) enable isolated acidic/alkaline regions in electrochemical devices, facilitating optimized environments for electrochemical separations and catalysis. For economic viability, BPMs must attain stable, high current density operation with low overpotentials in a freestanding configuration. We report an asymmetric, graphene oxide (GrOx)-catalyzed BPM capable of freestanding electrodialysis operation at 1 A cm–2 with overpotentials <250 mV. Use of a thin anion-exchange layer improves water transport while maintaining near unity Faradaic efficiency for acid and base generation. Voltage stability exceeding 1100 h with an average drift of 70 μV/h at 80 mA cm–2 and 100 h with an average drift of −300 μV/h at 500 mA cm–2 and implementation in an electrodialysis stack demonstrate real-world applicability. Continuum modeling reveals that water dissociation in GrOx BPMs is both catalyzed and electric-field enhanced, where low pK a moieties on GrOx enhance local electric fields and high pK a moieties serve as active sites for surface-catalyzed water dissociation. These results establish commercially viable BPM electrodialysis and provide fundamental insight to advance design of next-generation devices.

Lucas, Éowyn↗

DOE OSTI2024

Induced Polarization of Clayey Rocks and Soils: Non‐Linear Complex Conductivity Models

Abstract The past decades have witnessed the increased applications of induced polarization (IP) method in the critical zone studies with ubiquitous clay minerals. Although IP outperforms traditional electrical and electromagnetic methods through its unique ability to measure quadrature conductivity, the nonlinearity that quadrature conductivity behaves with salinities and frequencies greatly tortures IP practitioners, as (a) salinity‐dependency makes the quadrature conductivity a varyingly unstable parameter to quantitatively estimate hydraulic properties and clay content; (b) frequency‐dependent Cole‐Cole and Debye/Warburg decomposition models, although mathematically sound, physically mingle the properties of pore water and clay minerals and are empirical in nature. From basic principles, we demonstrate that quadrature conductivity remains a hybrid property involving both clay and water, and develop relevant models to distinguish them. Our models are validated by theories, experiments, simulations, and comparisons, all of which proclaim considerable advantages over previous models and offer the prospect of quantitative applications.

58 GEOSCIENCES↗