Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Heat conduction modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Assessing the Impact of Effective Thermal Conductivity on Gas-cooled Reactor Transients in RELAP5-3D

High-temperature gas-cooled reactors (HTGRs) are a relatively mature advanced reactor concept that are of interest to industry and government for near-term deployment. These systems feature passive safety through low power density, coated particle fuels, and large graphite volumes that lead to slow heatup. In loss of flow transients, cooldown is achieved through conduction and radiation heat transfer. In prismatic block-type reactors, such as the mHTGR-350, the presence of coolant holes and fuel compacts alters the flow of heat through the blocks. To capture this effect in systems codes like RELAP5-3D, relationships for effective thermal conductivity (ETC) must be used to appropriately degrade the thermal conductivity. This work presents an assessment of the impact of ETC on performance of block-type gas-cooled reactors in a pressurized conduction cooldown (PCC) and depressurized conduction cooldown (DCC) transients. Models for PCC and DCC with bulk material thermal conductivity and ETC were run for both the mHTGR-350 and the High-Temperature Test Facility (HTTF) to determine the impact of ETC on transient performance.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Lattice Anisotropy-Driven Reduction of Phonon Velocities in Black Phosphorus

Phonon dynamics and transport determine how heat is utilized and dissipated in materials. In 2D systems for optoelectronics and thermoelectrics, the impact of nanoscale material structure on phonon propagation is central to controlling thermal conduction. Here, we directly observe in-plane coherent acoustic phonon propagation in black phosphorus (BP) using ultrafast electron microscopy. We identify a significant reduction of the group velocities in directions intermediate to the armchair and zigzag lattice directions. Using a machine learning-based model with an >8000 atom supercell, we find that this slowing results from the mixing of in-plane transverse and longitudinal acoustic phonons and is independent of broken symmetries of edge reconstructions. In conclusion, this work demonstrates how coherent phonon transport is sensitive to propagation direction in the lattice plane.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Impact of Urbanization on Convection, Lightning, and Precipitation Over the Houston Metropolitan Area: Case Study Simulation From the TRACER Campaign

This study investigates the effects of urbanization, specifically land use change and anthropogenic emissions (AE), on convection, lightning, and surface precipitation for a case of summertime sea‐breeze convection observed over the Houston metropolitan area. The unique capabilities of the NASA‐Unified Weather Research and Forecasting model allows us to conduct a series of sensitivity experiments with complex configurations, in particular including multi‐year land model spin‐up simulations, treatment of aerosols and their precursors, and explicit cloud charging and lightning. The simulation results show that urban land use primarily alters the temporal evolution of convection, lightning, and surface precipitation, leading to late afternoon thunderstorm development. The decrease in latent heat flux from the land surface caused by urbanization weakens convection in the early afternoon, while a condition suitable for convection development is maintained in the late afternoon due to less stabilization of the lower troposphere by the weaker convection development and high sensible heat flux from the surface. On the other hand, anthropogenic aerosols directly enhance convection, lightning, and surface precipitation by increasing convective updrafts due to the aerosol‐induced convective invigoration. The combined effects of urban land use and AE lead to even stronger thunderstorms in the late afternoon, mostly consistent with observations. These results indicate that urbanization increases the probability of late afternoon thunderstorms over the Houston area during the summer season. Advanced weather forecasting models that incorporate these urbanization effects might support sustainable urban planning to better mitigate the impacts of urbanization on local weather and public safety.

Iguchi, T. [Univ. of Maryland, College Park, MD (U↗
DOE OSTI2024

Reversed-direction 2-point modelling applied to divertor conditions in DIII-D *

Abstract A predictive form of the extended 2-point model known as the ‘reverse 2-point model’, Rev2PM, is applied to a range of detachment levels in the open lower divertor of DIII-D, showing that the experimentally measured electron temperature ( T e ) and pressure ( p e ) at the divertor entrance can be calculated within 50% from target measurements, if and only if a posteriori corrections for convective heat flux are included in the model. Unlike the standard 2-point model, the Rev2PM calculates upstream scrape-off layer (SOL) quantities (such as separatrix T e and p e ) from target conditions (such as T e and parallel heat flux), with volumetric power and momentum losses depending solely on target T e . The Rev2PM is tested against a database of DIII-D inter-ELM divertor Thomson scattering measurements, built from a series of 6 MW, 1.3 MA, LSN H-mode discharges with varied main ion density, drift direction, and nitrogen puffing rate. Measured target T e ranged from 0.4–25 eV over this database, and upstream T e ranged from 5–60 eV. Poor agreement is found between upstream measurements and Rev2PM calculations that assume purely conductive parallel heat transport. However, introducing a posteriori corrections to account for convective heat transport brings the Rev2PM calculations within 50% of the measured upstream values across the dataset. These corrections imply that up to 99% of the parallel heat flux is carried by convection in detached conditions in the DIII-D open lower divertor, though further work is required to assess any potential dependencies on device size or divertor closure.

2-point model↗
DOE OSTI2021

A Flexible and Generic Functional Mock-up Unit Based Threat Injection Framework for Grid-interactive Efficient Buildings: A Case Study in Modelica

Grid-interactive efficient buildings (GEBs) have been considered as an important asset to support the power grid reliability by utilizing the demand flexibility offered by GEBs. GEBs are enabled by advances in sensors and controls, and the communication between building equipment, whole buildings, and the grid. The integration of different building technologies and network-based communication system makes GEBs vulnerable to passive threats such as equipment failure and active threats such as cyber-attacks. Modeling and simulation is an effective way to evaluate the impact of threats on the system performance. This paper proposes a generic and flexible threat injection framework for commonly-used building energy simulators such as EnergyPlus and Modelica to support threat modeling and evaluation. This framework leverages functional mock-up unit (FMU) to develop a general modeling interface for threat injection and simulation. A numerical case study using Modelica as a building energy simulator is conducted to demonstrate the capability of the framework for supporting single/multiple-order threat modeling and simulation of a GEB. Four threats and their combinations are injected on a Modelica-based threat-free building energy and control system, including operating supply fan at its full speed, remotely cycling the chiller on and off, blocking the chiller from receiving the chilled water supply temperature setpoints, and hijacking the global zone air temperature setpoint. Simulation results show that the cyber-attack that leads to short-term signal blocking has small effects on the system operation due to the "self-healing" feature of the heating, ventilation, and air-conditioning (HVAC) interactive control system. The threat that takes control of resetting the global zone air temperature setpoints has the most adverse impact on the system energy use, peak power demand, thermal comfort and the provision of demand flexibility. The combination of four threats have aggregative effects on the system but the effects are less than the additive effects of the individual threat.

Fu, Yanyang↗
DOE OSTI2023

Thermohydraulic experiments on a supercritical carbon dioxide–air microtube heat exchanger

Heat exchangers are critical components in supercritical CO 2 Brayton cycles and provide necessary heat input, recovery, and dissipation. Tubular heat exchangers with unconventionally small tube sizes (tube diameters less than 5 mm) are promising components for supercritical CO 2 cycles and potentially provide excellent structural stability with wide scope of application. This paper provides essential design and fabrication guidelines for a compact shell-and-tube heat exchanger with microtubes (with an inner diameter equal to 1.75 mm). A heat exchanger test rig is used to evaluate the thermohydraulic performance of this heat exchanger with supercritical CO 2 and air as working fluids. Thermohydraulic data are reported for more than forty sets of experiments with varying Reynolds numbers for shell and tube flows. Critical performance metrics are calculated from the data and compared with predictions from a previously described numerical model. The average deviations between the experimental and model results fall within 10% for all critical metrics. This excellent agreement validates the numerical model for supercritical CO 2 heat exchanger optimization and scale-up.

42 ENGINEERING↗
DOE OSTI2020

Validation of Hygrothermal Simulations with Wall Performance Experiments in an Environmental Chamber

Oak Ridge National Laboratory is developing a web tool, built on a rule-based expert system, that aids stakeholders in designing energy-efficient moisture-durable walls. Currently, the tool’s expert system database is mostly populated with expert opinion, but work is being done to incorporate durability assessment based on stochastic hygrothermal modeling. Chamber experiments have been conducted to validate the hygrothermal models. The measured temperature and humidity have been compared with those predicted by the hygrothermal model. Experiments and comparison to one-dimensional hygrothermal modeling results were completed with two walls, a Structural Insulated Panel (SIP) based wall and a Concrete Masonry Unit (CMU) based wall. The two walls were succumbed to typical Chicago weather in ORNL’s Heat, Air and Moisture chamber. The walls were sequentially exposed to different scenarios, including diffusion, wetting, air leakage and solar radiation. For the most part, the hygrothermal simulations can be successfully used to predict the performance of these wall systems for the moisture transport phenomenon. Errors between measured and simulated values decreased as measurements got closer to the interior side of the wall. The root mean square error was larger for relative humidity (up to 17.5%-RH for CMU) than for temperature (up to 4.5°C for CMU wall). The errors were larger for the CMU wall than the wood frame wall. The phenomenon, including liquid water, caused large discrepancies between measurement and simulation results, and simulated results showed slower drying of materials than measured results. The one-dimensional nature of the simulation model made simulating air leaks difficult but not impossible.

Salonvaara, Mikael↗
DOE OSTI2022

Thermal conductivity of multilayer polymer-nanocomposite thin films

The development of electrical insulators that are thermally conducting is critical for thermal management applications in many advanced electronics and electrical devices. Here, we synthesized polymer nanocomposite (PNC) films composed of polymers [polyethylenimine, poly(vinylamine), poly(acrylic acid), and poly(ethylene oxide)] and dielectric fillers (montmorillonite clay and hexagonal boron nitride) by layer-by-layer technique. The cross-plane thermal conductivity (k⊥) of the film was measured by the 3ω method. The effect of various factors such as film growth, filler type, filler volume fraction, polymer chemical structures, and temperature on the thermal conductivity is reported. The k⊥ of PNCs with thickness from 37 nm to 1.34 μm was found to be in the range of 0.11 to 0.21 ± 0.02 W m−1 K−1. The k⊥ values were found to be lower than the constituent polymer matrix. The experimental result is compared with existing theoretical models of nanocomposite systems to get insight into heat transfer behavior in such layered films composed of dielectrics and polymers.

Physics↗
DOE OSTI2023

Quality Ranking of Unary Fluoride Salt Property Data in MSTDB-TP

Molten salt reactor developers rely on thermal property data to design, license and operate the reactors they are developing. The Molten Salt Thermal Database-Thermophysical Properties (MSTDB-TP) is being developed under the DOE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program and managed by Oak Ridge National Laboratory. The MSTDB-TP is intended to serve as a single source of measured thermophysical property values for a wide variety of molten salt systems available to researchers, molten salt reactor developers, and regulators. These properties include melting and boiling points, heat capacity, density, viscosity and thermal diffusivity and conductivity. Published measurements of molten salt properties are lacking for many salts and the data that are available are often inconsistent. This creates a challenge for MSR developers in determining which property values to select when designing their reactors. The MSDTDB-TP collects all available property data and down-selects to preferred data sets or correlations. It is the purpose of this work to apply a recently developed ranking system that indicates the quality of property values listed in the database to inform decisions about how to use the available data. Quality assessments and rankings are being applied to data in MSTDB-TP to provide an indication of the quality of the data available. A previous report detailed the ranking system that was followed.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Thermal performance and energy consumption validation of an occupied local government office building outfitted with ceiling tile phase change materials

Buildings present an opportunity for energy conservation and the modulation of peak energy demand through controlled Heating, Ventilation, and Air Conditioning (HVAC) energy use. The administrative and office building stock in the United States holds potential to achieve energy and demand savings through retrofits such as insulation, weatherization, and thermal energy storage. Specifically, there is a need to validate passive phase change material (PCM) applications in full scale in aging administrative buildings in the US to evaluate the energy benefits. Aim of this study was to conduct a whole building level thermal and energy validation of an operational building and explore an alternative method for evaluating energy efficiency. To accomplish this, the study employed PCMs in the drop ceiling and carry out an energy audit and on-site measurement of HVAC systems' energy demand and consumption. A full-scale EnergyPlus energy model, modeled by the authors, served as a baseline for evaluation. The results show that calibrated model's envelope temperature measures fall within the accepted errors. HVAC energy simulation results also fall within the accepted errors for monthly and hourly pre- and post- PCM retrofit electricity and natural gas data. The novelty of this study is that it employees energy scales per Heating Degree Hour and Cooling Degree Hour, in contrast to the commonly used Heating Degree Days and Cooling Degree Days as reported in the literature to analyze energy savings. These findings underscore the pivotal role of a calibrated model in assessing the efficacy of a singular energy measure, like a PCM-retrofitted ceiling, in an occupied office building.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2026

COGIHTO: Comprehensive Geothermal Integrated HVAC Retrofits in Disadvantaged Communities

Geothermal heating and cooling (GHC) systems offer one of the most efficient approaches to provide space conditioning and water heating in buildings. Due to the relatively stable temperature of the ground year-round, these systems can provide an efficient solution to decarbonize heating systems. Yet this technology has faced limited market adoption due to its high first cost. This project seeks to address this barrier by designing a community GHC system at a public housing development in New York, taking advantage of load diversity (in terms of the timing and amount of heating / cooling used per person) that a larger community offers to reduce the size and thus first cost of a GHC system. Project success would result in societal benefit from reduced energy and fossil fuel use, reduced local emissions, and improved comfort, which would positively impact the health and well-being this community. This report presents the findings from Phase 1 project supported by the Department of Energy’s Geothermal Technologies Office and includes the technical design of a GHC at a public housing development in New York. To support the success of a potential deployment of a GHC system, the team conducted community engagement and prepared an assessment of the local workforce needs. The project team, led by EPRI, considered a subset of NYCHA housing developments that were previously identified by NYCHA as suitable for GHC retrofits. This report presents the site selection criteria the team developed and the overall score assigned to each of the eight developments considered. One site was selected for focus of the technical design: Stapleton Houses in Staten Island, a mid-size development with roughly 700 apartments across 6 buildings. The site has an existing hot water hydronic loop served by natural gas boilers for space heating and water heating. The technical design partner assessed baseline energy consumption at the site and established thermal models of the apartment buildings to assist in the design. The team considered several GHC system configurations and selected a hybrid geothermal design fed from water-source heat pumps with auxiliary gas boilers in a central plant. This design best leveraged the existing site configuration and minimized costs due to additional electrical infrastructure upgrades. The report also details the feedback received from community engagement efforts and the findings from the workforce needs assessment and workforce transition plan. Members of the site maintenance staff and the resident advisory council were interviewed for their input and feedback about both the site’s existing heating/cooling infrastructure and how a construction project would impact resident life. This feedback was taken into consideration when completing the geothermal system design. Research was conducted into existing workforce development opportunities available to NYCHA residents. A workforce development plan was established using information from NYCHA’s resident economic empowerment and sustainability group, along with training material and curriculum from the International Ground Source Heat Pump Association (IGSHPA).

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Data-Enabled Predictive Control for Building HVAC Systems

Model predictive control is widely used as a control technology for the computation of optimal control inputs of building heating, ventilating, and air conditioning (HVAC) systems. However, both the benefits and widespread adoption of model predictive control (MPC) are hindered by the effort of model creation, calibration, and accuracy of the predictions. In this paper, we apply the data-enabled predictive control (DeePC) algorithm for designing controls for building HVAC systems. The algorithm solely depends on input/output data from the system to predict future state trajectories without the need for system identification. The algorithm relies on the idea that a vector space of all input–output trajectories of a discrete-time linear time-invariant (LTI) system is spanned by time-shifts of a single measured trajectory, given the input signal is persistently exciting. Closed-loop simulations using EnergyPlus are performed to demonstrate the approach. The simulated building modeled in EnergyPlus is a modified commercial large office prototype building served by an air handling unit-variable air volume HVAC system. Temperature setpoints of zones are used as control variables to minimize the HVAC energy cost of the building considering a time-of-use electricity rate structure. Furthermore, sensitivity analysis is conducted to gain insights into the effect of parameter tuning on DeePC performance. Simulation results are used to illustrate the performance of the algorithm and compare the algorithm with model-based MPC and occupancy-based setpoint controller. Overall, DeePC achieves similar performance compared to MPC for lower engineering effort.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Morphological Ordering of the Organic Layer for High-Performance Hybrid Thermoelectrics

Inorganic-organic hybrids, such as Te-PEDOT:PSS core/shell nanowires, have emerged as a class of promising thermoelectric materials with combined attributes of mechanical flexibility and low cost. However, the poorly understood structure-property relationship calls for further investigation for performance enhancement. Here, through precise treatments of focused electron beam irradiation and thermal annealing on individual Te-PEDOT:PSS nanowires, new, nonchemical mechanisms are introduced to specifically engineer the organic phase, and the measured results provide an unprecedented piece of evidence, confirming the dominant role of organic shell in charge transport. Paired with the Kang-Snyder model and molecular dynamics simulations, this work provides mechanistic insights in terms of heating-enabled morphological ordering of the polymer chains. The measured results show that thermal annealing on the 42 nm nanowire results in a ZT value of 0.78 at 450 K. In conclusion, through leveraging the interfacial self-assembly of the organic phase to construct a high electrical conductivity domain, this work lays out a clear framework for the development of next-generation soft thermoelectrics.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Update of ASHRAE Standard 140 Section 5.2 and Related Sections (BESTEST Building Thermal Fabric Test Cases)

This report documents new work intended to update the current building thermal fabric software-to-software comparative test specifications and example results of ANSI/ASHRAE Standard 140-2017, Sections 5.2.1, 5.2.2, 5.2.3, and related supporting sections. The Building Energy Simulation Test and Diagnostic Method (BESTEST) test specification and example results that provide the basis for these Standard 140 sections were originally published in 1995. The state of the art in building energy simulation modeling has advanced since then. Additionally, comments on the test specification over the years highlighted ambiguities, further necessitating the update to the original work. This new work was conducted for the U.S. Department of Energy (DOE) by Argonne National Laboratory (Argonne) and National Renewable Energy Laboratory (NREL), in collaboration with a working group of international experts associated with the American Society of Heating, Refrigerating, and Air-Conditioning Engineers (ASHRAE) Standing Standard Project Committee (SSPC) 140. SSPC 140 is responsible for ANSI/ASHRAE Standard 140. The new work, presented in this report, defined the test cases and specifications such that building energy simulation programs with time steps of one hour or less can perform the tests. Revisions to the test specification and update of example results follow the vetting process of iterative field trials established for previous BESTEST work that has been adapted for ASHRAE Standard 140.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Evaluating simplifications of subsurface process representations for field-scale permafrost hydrology models

Abstract. Permafrost degradation within a warming climate poses a significant environmental threat through both the permafrost carbon feedback and damage to human communities and infrastructure. Understanding this threat relies on better understanding and numerical representation of thermo-hydrological permafrost processes and the subsequent accurate prediction of permafrost dynamics. All models include simplified assumptions, implying a tradeoff between model complexity and prediction accuracy. The main purpose of this work is to investigate this tradeoff when applying the following commonly made assumptions: (1) assuming equal density of ice and liquid water in frozen soil, (2) neglecting the effect of cryosuction in unsaturated freezing soil, and (3) neglecting advective heat transport during soil freezing and thaw. This study designed a set of 62 numerical experiments using the Advanced Terrestrial Simulator (ATS v1.2) to evaluate the effects of these choices on permafrost hydrological outputs, including both integrated and pointwise quantities. Simulations were conducted under different climate conditions and soil properties from three different sites in both column- and hillslope-scale configurations. Results showed that amongst the three physical assumptions, soil cryosuction is the most crucial yet commonly ignored process. Neglecting cryosuction, on average, can cause 10 %–20 % error in predicting evaporation, 50 %–60 % error in discharge, 10 %–30 % error in thaw depth, and 10 %–30 % error in soil temperature at 1 m beneath the surface. The prediction error for subsurface temperature and water saturation is more obvious at hillslope scales due to the presence of lateral flux. By comparison, using equal ice–liquid density has a minor impact on most hydrological metrics of interest but significantly affects soil water saturation with an averaged 5 %–15 % error. Neglecting advective heat transport presents the least error, 5 % or even much lower, in most metrics of interest for a large-scale Arctic tundra system without apparent influence caused by localized groundwater flow, and it can decrease the simulation time at hillslope scales by 40 %–80 %. By challenging these commonly made assumptions, this work provides permafrost hydrology scientists an important context for understanding the underlying physical processes, including allowing modelers to better choose the appropriate process representation for a given modeling experiment.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Simulation of electron Bernstein waves using FullWave with a 2D non-local hot plasma model

Hot plasma wave simulation capability is expanded in the FullWave code by updating the hybrid iterative solver in the code with a semi-implicit time stepping method. The new approach is used to simulate Electron Bernstein Wave (EBW) heating in over-dense spherical tokamak plasmas. The code’s hybrid iterative solver circumvents the prohibitive memory cost of direct methods by combining a time evolution of Maxwell’s equations with frequency-domain relaxation, while the conductivity kernel, calculated via 3D particle tracking, captures the essential non-local wave–particle interactions. One-dimensional EBW simulations verify the algorithm’s accuracy by demonstrating mode conversion from X-mode wave to EBW at the upper hybrid resonance and a strong cyclotron damping near the plasma core. Two-dimensional simulation reproduces the predicted short EBW wavelength and quantitatively matches the hot-plasma dispersion relation. This study demonstrates the fidelity of the hybrid solver for the electron cyclotron frequency range.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Turbulence regimes in the nocturnal roughness sublayer: Interaction with deep convection and tree mortality in the Amazon

We investigated the influence of seasonality and proximity to the forest canopy on nocturnal turbulence regimes in the roughness sublayer of a Central Amazon forest. Since convective systems of different scales are common in this region, we also analyzed the effect of extreme wind gusts (propagated from convective downdrafts) on the organization of the turbulence regimes, and their potential to cause the mortality of canopy trees. Our data include high-frequency winds, temperature and ozone concentration at different heights during the dry and wet seasons of 2014. In addition, we used critical wind-speed data derived from a tree-winching experiment and a modeling study conducted in the same study site. Two different turbulence regimes were identified at three heights above the canopy: a weakly stable (WS) and a very stable regime (VS). The threshold wind speeds that mark the transition between turbulence regimes were larger during the dry season and increased as a function of the height above the canopy. The turbulent fluxes of sensible heat and momentum during the WS accounted for 88% of the entire nighttime flux. Downdrafts occurred only in the WS and favored a fully coupled state of wind flow along the canopy profile. The destructive potential of winds was four times higher than on nights without downdrafts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

What can the GRACE seasonal cycle tell us about lake-aquifer interactions?

Lake-aquifer interactions have been the subject of investigation and debate for decades. Traditional investigation methods include measurement of water flux across the groundwater-surface water interface, application of heat and environmental tracer methods, conducting numerical simulations of the water flow, and mass balance-based approaches. We first review, evaluate and describe the different approaches that have been applied for examining lake-aquifer interactions and propose an additional complementary approach. While the leakage from lakes and its slow movement through porous media (laminar flow) is well characterized by existing methods, modeling the rapid and turbulent flow through preferred pathways (faults/karst) remains a challenging task. To better understand the nature, and the full scale, of the lake aquifer interactions, and using Lake Nasser and the underlying Dakhla subbasin of the Nubian Sandstone Aquifer System (NSAS) in the Western Desert of Egypt (area: 0.66 x 10 6 km 2 ) as a test site, we examine, model, and correlate temporally and spatially, the variations in Gravity Recovery and Climate Experiment terrestrial water storage (GRACE TWS ) solutions, with precipitation, lake levels, area, and water volume. We review current GRACE applications in hydrology and present our novel approach. Findings include: firstly, large seasonal fluctuations (peak: Nov./Dec.; trough: July/Aug.) in Lake Nasser's surface water levels (average annual fluctuations: from 3.7 m to 7.5 m), area (average area: 3622 km 2 but up to 4530 km 2 ), and volume (annual average: 13.4 km 3 , and up to 34 km 3 ) are observed throughout years 2006 to 2015. These fluctuations are accompanied by an increase in GRACE TWS (average: 50 ± 13 mm/yr, up to 77 ± 18 mm/yr) over Lake Nasser and by a progression of a front of increasing GRACE TWS values (> 50 ± 13 mm) away from the lake reaching distances of up to 700 km some 3 to 5 months following peak lake level periods. The areas witnessing the seasonal increase in GRACE TWS display a progressive increase in phase and decrease in amplitude with distance from the lake. Secondly, the negligible precipitation over the Dakhla subbasin cannot account for the observed seasonal GRACE TWS patterns and neither can the leakage signal from Lake Nasser. Thirdly, overall similarities in the distribution of GRACE TWS seasonal spatial patterns are observed. During high lake level periods (e.g., Dec. 2007; Feb. 2015) additional, and more distant, areas from Lake Nasser saw increases in GRACE TWS and vice versa during low lake level periods (e.g., July 2006; Aug. 2010). These observations are consistent with Lake Nasser being the main source of modern recharge for the Dakhla subbasin and suggest a new conceptual model for the subbasin: a slow groundwater flow through a porous matrix and a rapid groundwater flow along a network of faults, fractures, and karst topography across the matrix. We suggest that the proposed conceptual model for the interaction between Lake Nasser and the Dakhla subbasin could be applicable to aquifers of similar geologic, climatic, and hydrologic settings worldwide and that approaches similar to those advocated here could be used to investigate the validity of this suggestion.

54 ENVIRONMENTAL SCIENCES↗