Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “building load”

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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Harmonic Signals Dataset

The Harmonic Signals Dataset (HSD) consists of one-dimensional time series data collected from current sensors with a goal of providing data for research and development of Non-Intrusive Load Monitoring (NILM) with high sample rate (800kHz) sensors. NILM seeks to detect and characterize electrical equipment operating within a facility by sensing changes in electrical current on the building power system associated with the equipment. To provide data with a known ground truth, but with the realism of a signal introduced within a building facility, a series of known ground truth signals were injected as voltage into the power system of a building and measurements of the signals at multiple locations across the building power system were made with electrical current sensors.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Peak cooling load shift capability of a thermal energy storage system integrated with an active insulation system in US climate zones

This study aims to evaluate a thermal energy storage (TES) system integrated with an active insulation system (AIS) to form a TES + AIS integrated wall system as a partition and as a secondary cooling system to shift the peak load and reduce cooling energy consumption. To understand and demonstrate its cooling performance, the TES + AIS integrated wall system was installed in an office building in Oak Ridge, Tennessee. To investigate the effect of the TES + AIS integrated wall system in a typical office building and various climate zones, the US Department of Energy’s prototype office building model was modified to accommodate the proposed system. In this work, results showed that the minimum size of the proposed system to achieve energy savings varies depending on climate conditions. The minimum size of the proposed system for cooling energy saving and shifting peak cooling demand in climate zones 2, 3, and 5 is 29.7 m 2 , whereas it is 44.6 m 2 in zones 4 and 6. By installing the minimum size of the proposed system, 11.3 % to 16.4 % of cooling energy can be shifted during discharge hours in a representative summer day.

25 ENERGY STORAGE↗
DOE OSTI2023

Cost analysis of distributed storage in AC and DC microgrids

Building and microgrid designs with highly-distributed electrical storage have potential advantages over today's conventional topologies with centralized storage. Herein this paper studies the capital cost benefits of several residential behind-the-meter distributed-storage topologies, including AC and DC versions of systems with load-packaged batteries and resilient sub-networks. The study begins by defining the block configuration of each topology. This work then develops a model for the cost of the power electronics necessary to interface with the storage elements. Finally, the analysis develops a model for the total cost of each storage topology, incorporating the installation and soft costs. The results suggest that while the cost of power electronics is lower in centralized topologies, the total cost is lower for distributed storage due to the avoided costs of installation and permitting. This paper also explores the benefits of load-packaged batteries for savings in electrical infrastructure.

25 ENERGY STORAGE↗
NASA NTRS2005

Simulating Vibrations in a Complex Loaded Structure

The Dynamic Response Computation (DIRECT) computer program simulates vibrations induced in a complex structure by applied dynamic loads. Developed to enable rapid analysis of launch- and landing- induced vibrations and stresses in a space shuttle, DIRECT also can be used to analyze dynamic responses of other structures - for example, the response of a building to an earthquake, or the response of an oil-drilling platform and attached tanks to large ocean waves. For a space-shuttle simulation, the required input to DIRECT includes mathematical models of the space shuttle and its payloads, and a set of forcing functions that simulates launch and landing loads. DIRECT can accommodate multiple levels of payload attachment and substructure as well as nonlinear dynamic responses of structural interfaces. DIRECT combines the shuttle and payload models into a single structural model, to which the forcing functions are then applied. The resulting equations of motion are reduced to an optimum set and decoupled into a unique format for simulating dynamics. During the simulation, maximum vibrations, loads, and stresses are monitored and recorded for subsequent analysis to identify structural deficiencies in the shuttle and/or payloads.

Cao, Tim T.↗
DOE OSTI2021

Hardware-Based Emulator with Deep Learning Model for Building Energy Control and Prediction Based on Occupancy Sensors’ Data

Heating, ventilation, and air conditioning (HVAC) is the largest source of residential energy consumption. Occupancy sensors’ data can be used for HVAC control since it indicates the number of people in the building. HVAC and sensors form a typical cyber-physical system (CPS). In this paper, we aim to build a hardware-based emulation platform to study the occupancy data’s features, which can be further extracted by using machine learning models. In particular, we propose two hardware-based emulators to investigate the use of wired/wireless communication interfaces for occupancy sensor-based building CPS control, and the use of deep learning to predict the building energy consumption with the sensor data. We hypothesize is that the building energy consumption may be predicted by using the occupancy data collected by the sensors, and question what type of prediction model should be used to accurately predict the energy load. Another hypothesis is that an in-lab hardware/software platform could be built to emulate the occupancy sensing process. The machine learning algorithms can then be used to analyze the energy load based on the sensing data. To test the emulator, the occupancy data from the sensors is used to predict energy consumption. The synchronization scheme between sensors and the HVAC server will be discussed. We have built two hardware/software emulation platforms to investigate the sensor/HVAC integration strategies, and used an enhanced deep learning model—which has sequence-to-sequence long short-term memory (Seq2Seq LSTM)—with an attention model to predict the building energy consumption with the preservation of the intrinsic patterns. Because the long-range temporal dependencies are captured, the Seq2Seq models may provide a higher accuracy by using LSTM architectures with encoder and decoder. Meanwhile, LSTMs can capture the temporal and spatial patterns of time series data. The attention model can highlight the most relevant input information in the energy prediction by allocating the attention weights. The communication overhead between the sensors and the HVAC control server can also be alleviated via the attention mechanism, which can automatically ignore the irrelevant information and amplify the relevant information during CNN training. Our experiments and performance analysis show that, compared with the traditional LSTM neural network, the performance of the proposed method has a 30% higher prediction accuracy.

Ye, Zhijing↗
DOE OSTI2026

Consumer safety-oriented scheduling of rotating power outages during heat waves

Extreme heat events have widespread effects on power systems, reducing available generation capacity, limiting transmission capabilities, and causing unusual demand patterns on the consumer side. As these combined effects expose bulk transmission systems to potential large-scale blackouts, utilities may be required to schedule and apply rotating outages, by temporarily and alternately disconnecting distribution substations to reduce overload. However, utilities lack mechanisms to inform these events, exacerbating the negative effects of heat waves on affected communities. This paper introduces a novel framework for scheduling rotating outages during heat waves while considering impacts on consumers’ safety. Instead of random sequential load shedding, we propose a methodology to rotate power outages considering a metric that quantifies the indoor overheating risk of groups of consumers during a power outage. The overheating risk is derived from a detailed building simulation using CityBES, where the buildings are modeled based on available data—use type, year built, floor area, number of stories, location—while presence of air conditioning and occupancy are calibrated from smart meter data. Based on the metric, an algorithm to schedule the rotating outages is applied to prioritize feeders for disconnection at each hour according to their overheating risk to meet a utility load reduction target. Applied to two substations and seven feeders in the Portland General Electric territory, the results show that this approach effectively leads to the lowest overheating risk during the resulting outage schedules, with an average 10.1% lower overheating compared to uninformed schedules.

Building thermal simulation↗
DOE OSTI2023

KG-Hub—building and exchanging biological knowledge graphs

Knowledge graphs (KGs) are a powerful approach for integrating heterogeneous data and making inferences in biology and many other domains, but a coherent solution for constructing, exchanging, and facilitating the downstream use of KGs is lacking. Here we present KG-Hub, a platform that enables standardized construction, exchange, and reuse of KGs. Features include a simple, modular extract–transform–load pattern for producing graphs compliant with Biolink Model (a high-level data model for standardizing biological data), easy integration of any OBO (Open Biological and Biomedical Ontologies) ontology, cached downloads of upstream data sources, versioned and automatically updated builds with stable URLs, web-browsable storage of KG artifacts on cloud infrastructure, and easy reuse of transformed subgraphs across projects. Current KG-Hub projects span use cases including COVID-19 research, drug repurposing, microbial–environmental interactions, and rare disease research. KG-Hub is equipped with tooling to easily analyze and manipulate KGs. KG-Hub is also tightly integrated with graph machine learning (ML) tools which allow automated graph ML, including node embeddings and training of models for link prediction and node classification.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Case Study: Resilience Benefits of Distributed Wind Against Fuel and Weather Hazards in Alaska

In this case study of St. Mary’s Village, Alaska, we present a resilience evaluation exercise. A resilience framework is employed to identify system characteristics, relevant metrics, and resilience hazards and to assess the performance against the hazards with and without a distributed wind system installed. The results show the resilience benefits provided by the distributed wind installation against fuel shortage hazards and cold weather hazards. The resilience benefits can be assigned monetary values, which will be highly dependent on actual circumstances of the hazard, but provide insight into value streams of distributed wind that are not usually considered. For example, the single 900 kW turbine was found to prevent an average of 14,643 kWh of load from being dropped during a two-day diesel fuel shortage event, which saved the community $447,592 from the prevented outages. This case study serves as an example for novel power system resilience analysis and builds understanding of resilience hazards that are common across many power systems.

Culler, Megan J.↗
DOE OSTI2021

Tri-Level Scheduling Model Considering Residential Demand Flexibility of Aggregated HVACs and EVs Under Distribution LMP

Residential loads, especially heating, ventilation and air conditioners (HVACs) and electric vehicles (EVs), have great potentials to provide demand flexibility which is an attribute of grid-interactive efficient buildings (GEB). Under this new paradigm, first, EV and HVAC aggregator models are developed in this paper to represent the fleet of GEBs, in which the aggregated parameters are obtained based on a new approach of data generation and least squares parameter estimation (DG-LSPE), which can deal with heterogeneous HVACs. Then, a tri-level bidding and dispatching framework is established based on competitive distribution operation with distribution locational marginal price (DLMP). Furthermore, the first two levels form a bilevel model to optimize the aggregators’ payment and to represent the interdependency between load aggregators and the distribution system operator (DSO) using DLMP, and the third level is to dispatch the optimal load aggregation to all residents by the proposed priority list-based demand dispatching algorithm. Finally, case studies on a modified IEEE 33-Bus system illustrate three main technical reasons of payment reduction due to demand flexibility: load shift, DLMP step changes, and power losses. They can be used as general guidelines for better decision-making for future planning and operation of demand response programs.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRSDate not supplied

Wind Tunnel to Full Scale Mapping of Winds and Loads for Launch-Vehicle Ground Wind Loads

A launch vehicle ground-wind-loads program was conducted at the NASA Langley Transonic Dynamics Tunnel. The objectives were to quantify key aerodynamic and structural characteristics that impact the occurrence of large wind-induced oscillations of a launch vehicle when exposed to ground winds prior to launch. Of particular interest is the dynamic response of a launch vehicle when a von Kármán vortex street forms in the wake of the vehicle resulting in quasiperiodic lift and drag forces.Vehicle response to these quasiperiodic forces can become quite large when the frequency of vortex shedding nears that of a lowly-damped structural mode thereby exciting a resonant response. The study of ground-wind-loads presents unique challenges to quantify significant characteristics of the approaching wind and to relate the wind-tunnel acquired static, dynamic, and gravitational loads to full-scale vehicles. This paper will explain the correlation process between model-scale and full-scale wind characteristics and resulting structural loads. Characterization of the wind at a launch site requires knowledge of the vehicle dynamics and the anemometer performance. Vehicle dynamics dictate the frequency range of interest in characterization of the atmospheric turbulence, and in defining the required data window-length to model dynamic oscillation build-up during wind-gust analysis. Knowledge of the anemometer performance is required to infer peak wind magnitudes since most anemometers cannot directly measure instantaneous peak values. Following proper wind-speed correlation, static and dynamic wind-tunnel loads must be converted to full-scale equivalent values to derive design and operational guidance. In addition to the direct computation of full-scale equivalent loads using relevant scaling parameters, the corrections applied to dynamic loads accounting for differences in structural damping will also be discussed. Furthermore, a phenomenon that manifests in many slender and flexible launch vehicles is a load magnification due to gravity resulting from mass offset under deflection. A method for scaling this contribution from model-scale to full-scale is presented. Finally, a comparison of wind-tunnel derived loads using the described methodologies to launch vehicle measured loads is presented.

ground wind loads↗
NASA NTRSDate not supplied

Wind-Tunnel to Full-Scale Mapping of Winds and Loads for Launch-Vehicle Ground Wind Loads

A launch vehicle ground-wind-loads program was conducted at the NASA Langley Transonic Dynamics Tunnel. The objectives were to quantify key aerodynamic and structural characteristics that impact the occurrence of large wind-induced oscillations of a launch vehicle when exposed to ground winds prior to launch. Of particular interest is the dynamic response of a launch vehicle when a von Kármán vortex street forms in the wake of the vehicle resulting in quasiperiodic lift and drag forces. Vehicle response to these quasiperiodic forces can become quite large when the frequency of vortex shedding nears that of a lowly-damped structural mode thereby exciting a resonant response. The study of ground-wind-loads presents unique challenges to quantify significant characteristics of the approaching wind and to relate the wind-tunnel acquired static, dynamic, and gravitational loads to full-scale vehicles. This paper will explain the correlation process between model-scale and full-scale wind characteristics and resulting structural loads. Characterization of the wind at a launch site requires knowledge of the vehicle dynamics and the anemometer performance. Vehicle dynamics dictate the frequency range of interest in characterization of the atmospheric turbulence, and in defining the required data window-length to model dynamic oscillation build-up during wind-gust analysis. Knowledge of the anemometer performance is required to infer peak wind magnitudes since most anemometers cannot directly measure instantaneous peak values. Following proper wind-speed correlation, static and dynamic wind-tunnel loads must be converted to full-scale equivalent values to derive design and operational guidance. In addition to the direct computation of full-scale equivalent loads using relevant scaling parameters, the corrections applied to dynamic loads accounting for differences in structural damping will also be discussed. Furthermore, a phenomenon that manifests in many slender and flexible launch vehicles is a load magnification due to gravity resulting from mass offset under deflection. A method for scaling this contribution from model-scale to full-scale is presented. Finally, a comparison of wind-tunnel derived loads using the described methodologies to launch vehicle measured loads is presented.

Ground wind loads↗
DOE OSTI2020

Prototype energy models for data centers

Data centers in the United States consume about two percent of the nation’s electricity. Because heat gains from IT equipment drive cooling demand, data centers offer unique opportunities for energy savings. However, no prototype energy model for data centers is available in the suite of existing U.S. Department of Energy’s Commercial Prototype Building Models. Here we present the development of two new data center prototype models and their implementation in OpenStudio and EnergyPlus. The small-size data center model represents a computer room in a building served by computer room air conditioners (CRACs); while the large-sized model represents stand-alone data centers served by computer room air handlers (CRAHs) with a central chiller plant. For each data center model, two levels of IT equipment (ITE) load density were considered, to cover the wide range of IT power density of data centers: 40 and 100 W/ft 2 (430 and 1076 W/m 2 ) for the computer room, and 100 and 500 W/ft 2 (1076 and 5382 W/m 2 ) for the stand-alone data center. All other assumptions, such as building envelope, lighting, HVAC efficiencies and schedules, were based on the minimal requirements of ASHRAE Standard 90.1 at various vintages. We introduced a novel concept of supply and return air approach temperatures to capture the essential effects of non-uniform airflow and temperature distribution in data centers. The approach temperatures were pre-computed by computational fluid dynamics (CFD) simulations for various configurations of ITE loads and airflow containment management in data centers. A new feature was developed in EnergyPlus to implement the approach temperature method. A case study was conducted to demonstrate the use of the data center models. The two data center models cover all U.S. climate zones and can be used to evaluate energy saving measures for data centers, as well as to support development of data center energy efficiency codes and standards.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRS2000

Fracture Tests of Etched Components Using a Focused Ion Beam Machine

Many optical MEMS device designs involve large arrays of thin (0.5 to 1 micron components subjected to high stresses due to cyclic loading. These devices are fabricated from a variety of materials, and the properties strongly depend on size and processing. Our objective is to develop standard and convenient test methods that can be used to measure the properties of large numbers of witness samples, for every device we build. In this work we explore a variety of fracture test configurations for 0.5 micron thick silicon nitride membranes machined using the Reactive Ion Etching (RIE) process. Testing was completed using an FEI 620 dual focused ion beam milling machine. Static loads were applied using a probe. and dynamic loads were applied through a piezo-electric stack mounted at the base of the probe. Results from the tests are presented and compared, and application for predicting fracture probability of large arrays of devices are considered.

Kuhn, Jonathan, L.↗
NASA NTRS2011

On Calculation Methods and Results for Straight Cylindrical Roller Bearing Deflection, Stiffness, and Stress

The purpose of this study was to assess some calculation methods for quantifying the relationships of bearing geometry, material properties, load, deflection, stiffness, and stress. The scope of the work was limited to two-dimensional modeling of straight cylindrical roller bearings. Preparations for studies of dynamic response of bearings with damaged surfaces motivated this work. Studies were selected to exercise and build confidence in the numerical tools. Three calculation methods were used in this work. Two of the methods were numerical solutions of the Hertz contact approach. The third method used was a combined finite element surface integral method. Example calculations were done for a single roller loaded between an inner and outer raceway for code verification. Next, a bearing with 13 rollers and all-steel construction was used as an example to do additional code verification, including an assessment of the leading order of accuracy of the finite element and surface integral method. Results from that study show that the method is at least first-order accurate. Those results also show that the contact grid refinement has a more significant influence on precision as compared to the finite element grid refinement. To explore the influence of material properties, the 13-roller bearing was modeled as made from Nitinol 60, a material with very different properties from steel and showing some potential for bearing applications. The codes were exercised to compare contact areas and stress levels for steel and Nitinol 60 bearings operating at equivalent power density. As a step toward modeling the dynamic response of bearings having surface damage, static analyses were completed to simulate a bearing with a spall or similar damage.

Krantz, Timothy L.↗
DOE OSTI2021

Transactive HVAC Agent - Design and Performance Evaluation

Transactive energy systems are playing an increasingly important role in the efficient and reliable marketbased operation of the power grid. Since a significant portion of the residential building energy consumption is from heating ventilation and air conditioning (HVAC) systems, HVAC is one of the most promising resources to provide load flexibility. However, utilizing HVAC flexibility to provide various grid services while simultaneously maintaining consumer comfort and cost-reductions is challenging. This paper presents a design of a transactive HVAC agent (T-HVAC) to be used as a supervisory control for the HVAC system that can simultaneously ensure comfort and cost-reduction. In particular, the T-HVAC a) estimates HVAC thermal dynamics, b) ensures optimal operations of the HVAC system, and c) participates into markets, and d) implements a market-based control via controlling the thermostat temperature set-point. The T-HVAC performance is demonstrated through multiple scenarios and illustrations.

Demand flexibility, distribution system, HVAC, Tra↗
DOE OSTI2021

Maximum Cs-137 Curie Loading onto Crystalline Silicotitanate for the Documented Safety Analysis of the Tank Side Cesium Removal Platform

The Tank Side Cesium Removal (TSCR) system is currently being constructed to process Hanford tank waste supernates for vitrification. TSCR incorporates a filtration system and cesium (Cs) removal system using columns filled with crystalline silicotitanate (CST) ion exchanger, produced by Honeywell UOP, LLC. The documented safety analysis (DSA) developed for TSCR limits a single column curie loading to 141,600 Ci; given a 137 Cs isotopic mass fraction of 20% and the planned CST bed size of a TSCR column, this equates to 0.10 mmole Cs per g CST. Factors that influence 137Cs loading onto the CST include, but are not limited to, CST production lot (different production lots behave differently), contact temperature, contact duration, 137Cs mass fraction, and competitors in the tank waste feed. Seventeen tank waste feeds (compositions) were identified by Washington River Protection Solutions to be processed through TSCR. These feed compositions were used to develop a simulant (referred to herein as Stage 1) that would provide an upper bound to the Cs loading onto CST based on maximizing the Cs/Na activity coefficient ratios in solution while maintaining Na at no less than 5.0 M. Building upon this Stage 1 simulant, a series of four additional simulants were developed based on the cationic/anionic species that impact Cs exchange, with each successive formulation relaxing one or more matrix component concentration constraints as show in Table S.1

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

DECOVALEX-2023: Task C Final Report

The Full-scale Emplacement (FE) heater experiment at the Mont Terri Underground Rock Laboratory (URL) was designed and conducted by Nagra to replicate an emplacement tunnel of Nagra’s reference repository design at 1:1 scale. Alongside testing the technical feasibility of constructing disposal tunnels, emplacing waste containers in the tunnels and then backfilling them, the main goals of the FE experiment are (1) to obtain a better understanding of the coupled effects of induced thermo-hydro-mechanical (THM) processes that may occur and (2) to validate existing coupled THM models (Müller et al., 2017). A key aspect of ensuring safety for repositories located in low-permeability rock involves minimizing any damage to the rock itself, thereby preserving its integrity and promoting a stable environment Amongst a number of processes that could damage the rock is the increase in pore pressure due to thermal loading caused by heat emitted from the waste. To reduce the potential damage of the rock, it is important to analyse the evolution of heat over time due to the heat load of the containers and assess possible consequences by coupled THM models. The aim of Task C of DECOVALEX-2023 was to build 3D numerical models of the FE experiment, focussing in particular on the heating induced pore pressure change in the Opalinus Clay. Data from a large number of sensors were available from the FE experiment for model comparison. These sensors measured temperature and relative humidity in the bentonite around the heaters, and temperature, pressure and displacement/strain in the surrounding Opalinus clay. Data were available from the start of excavation (April 2012) up to August 2020 for most sensors (more than 5 years from the start of heating in December 2014). To fulfil the overall aim of the task, the work was broken down into a number of steps, starting with simpler models to build confidence in each team’s approach and then moving to more complex models that better represent the FE experiment. Step 0 consisted of 2D benchmark models, gradually increasing the number of processes that are represented from thermal (T) only models in Step 0a, to coupled thermal hydraulic (TH) models in Step 0b with a representation of changing porosity, to coupled thermo-hydro-mechanical (THM) models in Step 0c, where porosity changes are calculated by the mechanical model. A detailed specification of processes, parameters, initial and boundary conditions was provided for this step, with the ambition that all teams would work towards close agreement in their model results, thus building confidence in the model implementations. vi It was not straightforward to achieve agreement between the teams, so additional steps (Step 0b2, 0b3, 0c2, 0c3) were added along with derivation of some analytical solutions against which the models could be compared. The reasons for the differences between teams were investigated and found to be caused primarily by different conceptual model assumptions (including temperature dependence of the thermal expansion of water), different model formulations (including porosity evolution) and differences in modelled domain sizes, boundary conditions and grid discretisation. This demonstrates that comparisons between multiple modelling teams and/or comparison with analytical results and experimental data are highly beneficial in providing an indication of uncertainty in model predictions. At the conclusion of Step 0, almost all teams had achieved a close agreement in model results and those that had not achieved an agreement knew the reason for this. Step 1 moved from 2D models to 3D models of the FE experiment without adding technical features like shotcrete or EDZ, and only considering the heating phase. Initially the 3D model was tightly specified to continue to build confidence in the model implementations (Step 1a). The results of Step 1a were compared to the data from the FE-experiment without the teams seeing the data. The teams were then provided with a sub-set of the data from the FE-experiment and invited to consider how best to use the large dataset for model comparison (Step 1b). Teams were then asked to use the data provided to calibrate their models, only changing material property values rather than adding features or processes to their models (Step 1c). In Step 1, teams were asked to only model the heating phase of the experiment, so pressure in the Opalinus Clay was reported as change in pressure since the initial conditions were specified rather than modelled. The change from 2D to 3D models was accompanied by an increase in the dispersion of results between the teams. Some of this was resolved during the task, but some remained and is potentially due to model discretisation. Calibration of parameters was useful in improving the fit of the models to the data but the remaining differences indicated that the models were missing features or processes. In Step 2, the teams were asked to update their models with additional features and processes as well as calibrating parameters to try and improve the fit of the models to the data. Teams were encouraged to represent ventilation of the open FE tunnel prior to backfilling with heaters and bentonite and in Step 2, the absolute pressure in the Opalinus Clay was compared between the teams. Teams took different approaches, but there was consideration of adding shotcrete and an EDZ into the model, representing stress change during excavation and different approaches to modelling ventilation of the FE tunnel. Overall, the documented results showed a very good agreement for temperature. The results for porewater pressure evolution showed a significant improvement for most teams compared to Step 1c with a good agreement to the measurements for several teams whereas some teams overpredicted the pressure increase and others overpredicted the drainage effect especially for the sensors close to the heater. Step 3 was an opportunity for teams to use the models developed in Step 1 and Step 2 to make predictions about the temperature and pressure changes that will be expected at the FE experiment over the next few years in light of the planned changes in thermal output of the heaters.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗