Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Building thermal dynamics”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Investigating shock-induced chemical reactions in Ni+Al multilayers: A continuum-based mesoscale approach with Arrhenius kinetics and artificial thermal conduction

A new continuum-based mesoscale modeling approach for shock-induced chemical reactions (SICRs) in Ni+Al multilayers is demonstrated in Sandia’s shock physics hydrocode, CTH. The approach utilizes Arrhenius-type kinetics and artificial thermal conduction. Our work builds upon previous efforts to parameterize equations of state for Ni x Al y [J. Appl. Phys. 137, 075102 (2025)], as well as simulations of inert shocks in realistic 2D microstructures [J. Appl. Phys. 137, 225301 (2025)]. To calibrate the reaction kinetics, pairs of the reaction coordinate, R′, vs time are extracted from the molecular dynamics (MD) literature. Here, the MD-informed kinetics are used to simulate the dynamic evolution of pressure and temperature in 2D mesoscale simulations. Overall, the MD-informed kinetics obtained for planar interfaces are too slow, as initial reaction is not observed on a nanosecond time scale. Even with quasi-periodic shock focusing leading to the formation of so-called “hot-spots,” the hot spots are unable to grow and coalesce using the fitted Arrhenius rate constants. However, by increasing the rate constants by two orders of magnitude, SICRs are observed at a shock pressure near 30 GPa, which is supported by experiments. Consequently, these mesoscale simulations suggest that unresolved shear-based mechanical mixing might possibly account for the discrepancies in kinetic rates, with shock-generated intense perturbations, interfacial vortical flows, and elevated temperatures serving as favorable reaction conditions. Future work will calibrate a shear-dependent reaction rate from the MD simulations with realistic interfaces that are wavy, diffuse, and disordered.

Kittell, David E. [Sandia National Laboratories (S↗
DOE OSTI2025

Decision and Control of Complex Systems – A Data-Drive Framework

During the project period, we have collaborated with other team members and developed novel algorithms for novelty detection, continual learning, and graph learning algorithms for dynamic systems. The results are documented in publications and meeting notes. Moreover, we leverage virtual collaboration tools (such as Basecamp, Microsoft Teams and Zoom) for technical exchanges. Our research on novelty detection was published at AAAI 2022 and Lecture Notes in Artificial Intelligence, Springer Nature. The newly developed algorithms were successfully applied to realistic cases, including thermal data from buildings at Pacific Northwest National Lab and microelectronic data provided by GlobalFoundries. Multiple publications have been produced from this project, in collaboration with other team members. Three PhD students were supported in this project to conduct their research.

42 ENGINEERING↗
DOE OSTI2021

Thermal transport in defective and disordered materials

With significant recent advancements in thermal sciences—such as the development of new theoretical and experimental techniques, and the discovery of new transport mechanisms—it is helpful to revisit the fundamentals of vibrational heat conduction to formulate an updated and informed physical understanding. The increasing maturity of simulation and modeling methods sparks the desire to leverage these techniques to rapidly improve and develop technology through digital engineering and multi-scale, electro-thermal models. With that vision in mind, this review attempts to build a holistic understanding of thermal transport by focusing on the often unaddressed relationships between subfields, which can be critical for multi-scale modeling approaches. For example, we outline the relationship between mode-specific (computational) and spectral (analytical) models. We relate thermal boundary resistance models based on perturbation approaches and classic transmissivity based models. We discuss the relationship between lattice dynamics and molecular dynamics approaches along with two-channel transport frameworks that have emerged recently and that connect crystal-like and amorphous-like heat conduction. Throughout, we discuss best practices for modeling experimental data and outline how these models can guide material-level and system-level design.

42 ENGINEERING↗
DOE OSTI2022

Numerical modeling and parametric study of a dual purpose underground thermal battery

We report integrating thermal energy storage with building energy systems can enable flexible building electric demands at buildings to help mitigate the mismatch between electricity supply and demand. A novel building heating and cooling system that integrates a dual-source heat pump with hybrid thermal storage named dual-purpose underground thermal battery (DPUTB) has been developed for reshaping building electric demands. The proposed DPUTB integrated geothermal heat pump system is an original innovation that enables Grid-interactive Energy Efficient buildings. This paper focuses on the study of the novel DPUTB. The DPUTB works as both a thermal storage tank (an inner tank) and a ground heat exchanger (an outer tank separated from the inner tank by the insulation material). High fidelity and computationally effective models are needed to predict the performance of the novel DPUTB. This study has developed a simplified dynamic model for the DPUTB according to heat transfer and energy conservation principles and validated it by using experimental data obtained from testing a small-scale DPUTB apparatus. A parametric study was conducted to identify a design that can achieve the target thermal storage performance of load shift and energy efficiency. The parametric study results show that the inner tank shell thermal conductivity and the phase change material's melting point are the two most influencing factors on the performance of the DPUTB. One single full-size DPUTB with the identified design could provide 1-ton cooling (3.51 kW) with the supply temperature lower than 11 °C for 4 hours in summer after being fully charged in 8 hours. The inner tank filled with phase change material is for cooling thermal storage as a latent tank in the design. However, its capacity can be as high as 60 MJ as a sensible water tank for heating storage in winter. In the future study, the DPUTB model will be incorporated into the dual-source heat pump system for evaluating the overall system performance of demand side management in the long term.

25 ENERGY STORAGE↗
DOE OSTI2023

Degradation of Insulating Glass Units: Thermal Performance, Measurements and Energy Impacts

Insulating glass unit (IGU) degradation has been studied extensively. However, there is limited understanding of how present durability evaluation standards relate to product lifetime. Furthermore, there is debate on how to quantify performance of installed windows over time to better understand degradation processes. More knowledge on these topics is required to link durability evaluation to product lifetime projections based on energy performance. Energy models provide helpful estimations of total annual building energy consumption. However, most models are based on “as installed” performance of envelope components and fail to account for performance degradation. This can lead to an underestimation of building lifetime energy consumption. A better understanding of the relationship between durability and energy performance can inform integration of degradation dynamics into energy modeling software. This will improve lifetime building energy consumption estimations as well as inform appropriate retrofit strategies and timing. This paper reviews current durability literature, various standards for window performance ratings and weathering methods, existing in situ IGU energy performance measurement techniques, and whole-building energy effects. The challenges and disparities among various studies are analyzed and discussed. The authors hope that further work in this area will lead to the development of improved in situ test methods to assess IGU degradation in the field and link this knowledge to improved energy performance modeling approaches.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Building Energy Modeling Enhancements to Identify Least-Cost Pathways to Net-Zero Carbon Homes: Preprint

Residential grid-interactive efficient buildings (GEBs) can utilize high levels of energy efficiency and demand flexibility to deliver value to occupants, the grid, and society. However, without the ability to analyze, design, and optimize residential GEBs there is a risk that homes will not realize their full potential value as efficient and dynamic resources, which could ultimately lead to higher than necessary energy costs for U.S. households and lower realization of potential energy efficiency and demand flexibility benefits such as increased convenience/automation, thermal comfort, durability, indoor air quality, and resilience. In 2018, the National Renewable Energy Laboratory (NREL) Residential Modeling Team developed a vision for an open source Residential GEB Analytics Platform built within DOE's EnergyPlus and OpenStudio modeling environment that would enable the design and optimization of residential GEBs, including the identification of least-cost pathways to highly energy-efficient and energy-flexible homes (e.g., net zero carbon homes). We created a detailed workplan for the development of new and enhanced residential GEB component-level models for EnergyPlus and OpenStudio needed to progress toward the vision for the analytics platform while delivering near-term benefits to industry. This paper 1) presents the vision for the platform, summarizing the workplan for residential GEB modeling enhancements; 2) highlights significant advancements that have been achieved between 2018 and 2022, including stochastic residential occupancy modeling, flexible water heater modeling, detailed lithium-ion stationary battery modeling, and realistic residential HVAC modeling; and 3) outlines ongoing efforts and next steps toward the full vision for the platform.

building energy modeling↗
DOE OSTI2020

Elucidating the effect of preheating temperature on melt pool morphology variation in Inconel 718 laser powder bed fusion via simulation and experiment

In laser powder bed fusion (L-PBF) additive manufacturing, the mechanical performance, microstructure and defects of fabricated parts are closely associated with the melt pool morphology, e.g., its dimension and shape through the building process. Past studies have largely focused on how the process parameters such as laser power and scan speed affect melt pool characteristics. In this study, the melt pool morphology variation as a function of preheating temperature in the conduction, transition, and keyhole regimes and the underlying mechanisms in each regime are investigated through ex-situ sample characterization and computation thermal fluid dynamics (CtFD) simulation. Single tracks with different combinations of laser power and scan speed are deposited on an Inconel 718 bare plate preheated to a temperature range of 100–500 °C in the experiment. Significant changes are observed in melt pool morphology as a function of preheating temperature from optical measurements of melt track cross sections. The depth of melt pool in the three regimes increases monotonically with preheating temperature, e.g., at 500 °C, the experimental melt pool depth is increased by 49% in conduction regime, 34% in transition regime and 33% in keyhole regime, respectively, while the variation of melt pool width in each regime does not all follow an increasing trend but depends on the melt pool regimes. Melt pool width variation in the conduction and transition regimes is found to depend on the enhanced heat conduction directly related to temperature dependent thermal properties. Through validated CtFD simulations, it is found that in the keyhole regime the evaporation mass, recoil pressure, and laser drilling effect is enhanced with higher preheating temperature, which gives rise to a deeper melt pool. The simulations also reveal that preheating temperature significantly elongates the melt track length due to the increased flow rate and strong recoil pressure that accelerates the backward flow.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Implementing Ordinary Differential Equation Solvers in Rust Programming Language for Modeling Vehicle Powertrain Systems: Preprint

Efficient and accurate ordinary differential equation (ODE) solvers are necessary for powertrain and vehicle dynamics modeling. However, current commercial ODE solvers can be financially prohibitive, leading to a need for accessible, effective, open-source ODE solvers designed for powertrain modeling. Rust is a compiled programming language that has the potential to be used for fast and easy-to-use powertrain models, given its exceptional computational performance, robust package ecosystem, and short time required for modelers to become proficient. However, of the three commonly used (>3,000 downloads) packages in Rust with ODE solver capabilities, only one has more than four numerical methods implemented, and none are designed specifically for modeling physical systems. Therefore, the goal of the Differential Equation System Solver (DESS) was to implement accurate ODE solvers in Rust designed for the component-based problems often seen in powertrain modeling. DESS is a text-based software package that provides a flexible framework for building and solving systems of ODEs. This allows DESS to be included as a dependency for automotive powertrain models that require a variety of solvers and solver configurations. Seven explicit ODE solver methods have been implemented in DESS: Euler’s, Heun’s, midpoint, Ralston’s, classic Runge-Kutta, Bogacki-Shampine, and Cash-Karp. These represent five fixed-step methods and two adaptive-step methods. This paper shows that the solver implementations increase accuracy and computational efficiency compared to Euler's method when modeling a system of three thermal masses in Rust. DESS also includes features designed for modeling component-based physical systems. Users can define relationships between nodes in their system, which the package then translates into a system of equations, leading to simpler and more intuitive code. In the case of a three-thermal-mass system, the user can specify node thermal properties (e.g., thermal capacitance), how nodes are interconnected, and thermal conductance between nodes rather than providing a system of equations. The core contribution from this work is an open-source, text-based Rust package with ODE solvers for automotive powertrain modeling to support cost-free, fast, and accurate simulation.

ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2025

Stress evolution and creep deformation in solid-oxide electrolysis cell systems – Dynamic modeling and multi-objective optimization to maximize stack life and efficiency

Here, this study develops a thermal stress model of solid-oxide electrolysis cells (SOECs) including a model for creep strain and failure probability that is integrated with a dynamic plant-wide model of a hydrogen production process. Uncertainties in key material properties of the cell are quantified to assess their impact on stress profile variability. The oxygen electrode is found to have about 10 times higher failure probability compared to the fuel electrode. The study shows that if the stack operation is not optimized, cycling operation would lead to stress build-up eventually leading to catastrophic failure. A dynamic optimization problem is set up for obtaining the optimal operational profile considering a variable hydrogen production rate. Due to the tradeoff between the efficiency and stress build-up, the dynamic optimization problem is multi-objective. It is observed that the optimizer can considerably reduce the stress build-up (i.e., can increase the stack life) albeit at the cost of a lower efficiency thus exhibiting strong tradeoffs between capital and operating costs. For example, if the stack would be replaced in 0.5 yr, specific energy requirement would be 48.5 kWh/kg H 2 while for a stack replacement time of about 6 yr, the specific energy requirement rises by about 4.2 %.

SOEC↗
DOE OSTI2026

Systematic analysis of melt pool dynamics in laser processing of mixed powder feedstocks

Functionally graded materials (FGMs) fabricated via additive manufacturing of blended powders offer the potential to spatially tailor properties for new technologies, such as fusion first-wall systems, turbine blades, and spacecraft. However, processing these materials is difficult due to the multiplicity of processing parameters to optimize, all of which must be changed as substrate material, powder feedstock compositions, and melt pool dynamics evolve. Here, this work systematically evaluates the qualitative and quantitative effects of these variables on the melt pool size, shape, composition, and particle distribution in an exemplar Ti-Ta system, and connects the experimental results to Marangoni flow behavior and phenomena observed in other systems. Increasing laser power linearly increases melt pool size and layer thickness, driving engineering considerations such as part/geometrical tolerances. Decreasing laser velocity changes the melt pool shape from lenticular to convex and reduces chemical homogeneity due to extreme thermal and compositional gradients between the melt pool center and boundaries. Thermophysical property differences between the powder feedstock and substrate material, as well as the directionality of the gradient, affect dilution and melt pool dynamics, which in turn affect the melt pool boundary characteristics, shape, and uniformity. Mixed powder feedstocks of intermediate compositions do not behave according to linear interpolations between single-material endpoints, instead building taller and wider melt pools. As such, it is recommended to quantify process maps for at least one intermediate composition in the FGM or multi-material system of interest to ensure optimized processing parameters, predictable melt pool sizes and shapes, and compositional and spatial precision.

Dissimilar↗
DOE OSTI2026

Prime Time for Model-Predictive Control? Assessing the Technical and Market Readiness of Advanced Controls in Buildings

Despite three decades of extensive research and field testing that have consistently validated the benefits of Model Predictive Control (MPC) in building applications, the technology has seen limited market adoption. This paper evaluates the readiness of MPC for widespread deployment, showcases recent demonstrations and field tests across diverse building types, including residential, small commercial, large commercial, and campus settings. Our results demonstrate that MPC can optimize system operations to achieve load shifting, minimize curtailment of on-site generation, and reduce energy costs by up to 80 %, while maintaining or improving occupant comfort. We also show that MPC can effectively control large assets, such as MW-sized thermal storage systems, and respond to dynamic pricing signals. However, achieving scale remains difficult due to labor-intensive workflows, reliance on a “PhD-in-the-loop” for MPC design and maintenance, susceptibility to fragile data infrastructure, and persistent workforce education and acceptance barriers. To bridge this gap, we outline a transition from bespoke, labor intensive prototypes toward streamlined, segment-targeted deployment strategies that leverage model templates, semantic tools, and generative AI. By automating control configuration and reducing engineering effort, these recommendations provide a pathway for transforming successful research demonstrations into scalable, market ready solutions for MPC-based controls.

Pritoni, Marco↗
DOE OSTI2022

Protein Vesicles with pH-Responsive Disassembly

We report protein biomaterials offer several advantages over those made from other components because their amino acid sequence can be precisely controlled with genetic engineering to produce a diverse set of material building blocks. In this work, three different elastin-like polypeptide (ELP) sequences were designed to synthesize pH-responsive protein vesicles. ELPs undergo a thermally induced hydrophobic transition that enables self-assembly of different kinds of protein biomaterials. The transition can be tuned by the composition of the guest residue, X, within the ELP pentapeptide repeat unit, VPGXG. When the guest residue is substituted with an ionizable amino acid, such as histidine, the ELP undergoes a pH-dependent hydrophobic phase transition. We used pH-responsive ELPs with different levels of histidine substitution, in combination with leucine zippers and globular, functional proteins, to fabricate protein vesicles. We demonstrate pH-dependent self-assembly, diameter, and disassembly of the vesicles using a combination of turbidimetry, dynamic light scattering, microscopy, and small angle X-ray scattering. As the ELP transition is dependent on the sequence, the vesicle properties also depend on the histidine content in the ELP building blocks. These results demonstrate the tunability of protein vesicles endowed with pH responsiveness, which expands their potential in drug-delivery applications.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Frequency Control and Dynamics (Part 2) [Slides]

This presentation provides an introductory overview of system dynamics and frequency control in electric power systems, with a focus on concepts relevant to small and interconnected grids such as those in Malawi. It explains foundational principles of AC system frequency, the relationship between generation-demand balance and frequency deviations, and the operational limits of generators and end-use equipment. The deck discusses frequency stability within broader system stability classifications and illustrates how inertia and turbine-governor dynamics shape system response to disturbances. It then outlines the tiered approach to frequency control - primary, secondary, and tertiary - detailing the roles, characteristics, timescales, and response mechanisms of each. Special emphasis is placed on hydro and thermal unit behavior, area control error (ACE), automatic generation control (AGC), and the operational implications of interconnecting small systems with larger grids. The material was developed to support Malawi's electricity sector and the establishment of the Southern Africa Battery Energy Storage Center of Excellence (SABESS CoE).

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Failure investigation of nuclear grade POCO graphite target in high energy neutrino physics through numerical simulation

In many high energy physics neutrino beamlines, targets are made up of isotropic POCO graphite grade for production of neutrinos for high energy particle physics research and are bombarded with highly energetic pulsed proton beam. The pulsed proton beam with a small beam size creates thermal stress waves as well as radiation damage such as displacement damage, void formation, swelling, and gas formation amongst others. As a result, after a few years of operation one of the longest operated targets appears to have undergone bulk swelling and fractured along the beam centerline. A complex interaction of dynamic loading due to beam, material degradation due to irradiation are thought to have caused such a failure. Here, we present a numerical simulation to capture the combined effects of swelling and the dynamic effect of thermal stress waves due to beam heating, to explain the failure of target. An empirical formula has been developed to take into account swelling as a function of temperature and proton fluence and has been implemented in a commercial finite element code. Extensive X-ray diffraction (XRD) were performed on graphite samples from the failed target to understand the lattice parameter changes due to irradiation and build a valid empirical model. Simulation results show that swelling plays a major role in elevating the stress state in the material exceeding the failure strength while the pulsed beam heating also introduces fatigue loading that would explain brittle fracture emanating from inside of the material and propagating outward. Such empirical formulations will help in identifying the threshold values of critical parameters to extend the service life of future multi-megawatt neutrino targets and targets in other accelerator environments.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

One System, Many Models: Designing a Surrogate Model for Sulfur Thermal Energy Storage: Preprint

Industrial process heating (IPH) relies primarily on thermal energy generated by fossil fuel combustion to produce, treat, and alter manufactured goods. Thermal energy storage (TES) helps reduce the carbon footprint of IPH systems by facilitating the utilization of renewable and waste heat sources. A promising new TES technology uses elemental sulfur as the heat-storage medium. The design of sulfur TES systems can be evaluated with the aid of computational fluid dynamics (CFD). However, the computational cost of such CFD efforts is prohibitive to comprehensive optimizations over design parameters. To reduce this computational cost, machine learning (ML) models can be developed to act as surrogates for CFD. In this paper, we describe the process of building and evaluating surrogate ML models for facilitating optimization of sulfur TES systems for IPH. To enforce the thermodynamic relationship between the two modeled quantities, we develop a hybrid model for sulfur temperature using both direct predictions of temperature and calculations of temperature from predictions of the heat transfer coefficient. The hybrid model enforces this constraint at the expense of the slightly reduced accuracy compared to two disjoint models. The overall high accuracy observed in our model evaluation demonstrates the usefulness of such surrogate modeling for studying TES systems. This work contributes to the field of TES surrogate modeling by offering a novel accurate hybrid approach to predicting simultaneously the heat transfer coefficient and temperature of the heat storage medium.

computational fluid dynamics↗
DOE OSTI2022

Reconfiguring band-edge states and charge distribution of organic semiconductor–incorporated 2D perovskites via pressure gating

Two-dimensional (2D) semiconductor heterostructures are key building blocks for many electronic and optoelectronic devices. Reconfiguring the band-edge states and modulating their interplay with charge carriers at the interface in a continuous manner have long been sought yet are challenging. Here, using organic semiconductor–incorporated 2D halide perovskites as the model system, we realize the manipulation of band-edge states and charge distribution via mechanical—rather than chemical or thermal—regulation. Compression induces band-alignment switching and charge redistribution due to the different pressure responses of organic and inorganic building blocks, giving controllable emission properties of 2D perovskites. We propose and demonstrate a “pressure gating” strategy that enables the control of multiple emission states within a single material. We also reveal that band-alignment transition at the organic-inorganic interface is intrinsically not well resolved at room temperature owing to the thermally activated transfer and shuffling of band-edge carriers. This work provides important fundamental insights into the energetics and carrier dynamics of hybrid semiconductor heterostructures.

42 ENGINEERING↗
DOE OSTI2025

Beam Shaping and Oscillation for Metals Additive Manufacturing

Additive Manufacturing (AM) is pivotal for next-generation industrial manufacturing but is hampered by process control challenges leading to unwanted defects and sub-optimal material properties. Our research builds on prior work that utilized beam oscillation techniques to mitigate these issues (Wu et al., 2023). We focus on investigating the influence of beam oscillation frequency on melt pool dynamics, utilizing in-situ synchrotron material characterization from the ESRF (Chen et al., 2026). Preliminary findings indicated that beam oscillation results in unique fluid flow patterns and enhanced melting stability. These findings suggested the possibility of optimized AM process control, including minimized defects and tailored thermal histories. We proposed further research to explore the effects of beam oscillation frequencies comprehensively, using the in-situ characterization capabilities at the ESRF. The goal was to quantify this relationship, offering a pathway to more precise AM process control and optimization. Our research is expected to offer new scientific insights into the fluid mechanics and thermal behaviour in AM, contributing to its broader applicability and commercial viability.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A Survey of Multizone Airflow Simulation Software

H-Canyon simulates the building and process ventilation to ensure positive contamination control for any modifications as well as evaluating off-normal conditions. Simulations can also be used to identify air quality of personnel occupied areas. The models can be classified into two groups, single and multi-zone models. As the name states, a single zone model assumes that the building can be described by a single, well-mixed zone. Whereas large buildings have structures that must be described as multizone. Multizone airflow models deal with the difficulty of flows in a building by distinguishing the effects of internal flow restrictions. Several models have been developed that take internal partitions into account. CONTAM (currently used) is a multizone airflow program created by NIST to determine airflows including infiltration, exfiltration, wind pressures acting outside the building, buoyancy effects caused by the indoor and outdoor air temperature difference. COMIS is a multizone air infiltration, ventilation, and contaminant transport model that was developed at Lawrence Berkeley National Laboratory. COMIS deals with applied pressure distribution, airflow components, HVAC systems, schedules, and contaminant transport. EnergyPlus is designed by DOE. It supports flexible component level configuration of HVAC models and simulates sub-hourly time steps to control fast dynamic systems. The program features include providing accurate environmental performance data including energy consumption and carbon emissions, and reporting solar gains on surfaces. Heavent is used to model industrial ventilation systems to control airborne contaminants. It's designed based on the American Conference of Governmental Industrial Hygienists (ACGIH) Ventilation manual. DesignBuilder is the most recognized and advanced user interface to EnergyPlus. It gives access to all of the most frequently required simulation capabilities covering building fabric, thermal mass, glazing, shading renewables, and HVAC analysis. The objective of this study is to identify all relevant software available including commercially licensed, freely licensed, and open source, that can be used to simulate large, complex process facility ventilation systems. All discovered software tools are evaluated for salient features and compared to the software tool currently used to simulate the H-Canyon ventilation system. The design criteria of a ventilation system has two basic requirements: Supplying fresh air for the occupants and changing the air in the room sufficiently; therefore, the smells, fumes, and contaminants will be removed. In this study, first, whether the software simulates airflow (infiltration, natural ventilation, and interzonal) or considers them to be inputs is discussed. Second, if indoor contaminant levels are considered; and finally suggestions on how airflow analysis can be more appropriately accounted for when performing airflow simulation of a building. All of the airflow software reviewed for this study can account for constant infiltration and natural ventilation rates that are not affected by changes in indoor and outdoor conditions. Modelling infiltration in EnergyPlus or DesignBuilder are applicable for small buildings and not appropriate for larger buildings like the H-Canyon. The effect of wind on external pressures, infiltration, and natural and mechanical ventilation, can also be calculated using all of these software. These programs have control options to switch between natural to mechanical ventilation modes or both (hybrid). Interzonal airflow can be applied which means one zone can be specified to receive a fixed amount of air from another zone. The following ventilation modeling software is recommended for further testing: CONTAM, COMIS, Heavent. The next step in evaluating the software will consist of creating test facility simulations using evaluation copies of the software.

42 ENGINEERING↗