Double-heater-wire circuits and heat-and-moisture exchangers and the risk of ventilator-associated pneumonia

Objective: To compare the incidence of ventilator-associated pneumonia (VAP) in patients ventilated in intensive care by means of circuits humidified with a hygroscopic heat-and-moisture exchanger with a bacterial viral filter (HME) or hot-water humidification with a heater wire in both inspiratory and expiratory circuit limbs (DHW) or the inspiratory limb only (SHW). Design: A prospective, randomized trial. Setting: A metropolitan teaching hospital's general intensive care unit. Patients: Three hundred eighty-one patients requiring a minimum period of mechanical ventilation of 48 hrs. Interventions: Patients were randomized to humidification with use of an HME (n = 190), SHW (n = 94), or DHW (n = 97). Measurements and Main Results. Study end points were VAP diagnosed on the basis of Clinical Pulmonary Infection Score (CPIS) (1), HME resistance after 24 hrs of use, endotracheal tube resistance, and HME use per patient. VAP occurred with similar frequency in all groups (13%, HME; 14%, DHW; 10%, SHW; p = 0.61) and was predicted only by current smoking (adjusted odds ratio [AOR], 2.1; 95% confidence interval [CI], 1.1-3.9; p =.03) and ventilation days (AOR, 1.05; 95% Cl, 1.0-1.2; p =.001); VAP was less likely for patients with an admission diagnosis of pneumonia (AOR, 0.40; 95% Cl, 0.4-0.2; p =.04). HME resistance after 24 hrs of use measured at a gas flow of 50 L/min was 0.9 cm H2O (0.4-2.9). Endotracheal tube resistance was similar for all three groups (16-19 cm H2O min/L; p =.2), as were suction frequency, secretion thickness, and blood on suctioning (p =.32, p =.06, and p =.34, respectively). The HME use per patient per day was 1.13. Conclusions: Humidification technique does not influence either VAP incidence or secretion characteristics, but HMEs may have air-flow resistance higher than manufacturer specifications after 24 hrs of use.

A techno-economic assessment of the use of fast pyrolysis bio-oil from UK energy crops in the production of electricity and combined heat and power

This thesis investigates the cost of electricity generation using bio-oil produced by the fast pyrolysis of UK energy crops. The study covers cost from the farm to the generator’s terminals. The use of short rotation coppice willow and miscanthus as feedstocks was investigated. All costs and performance data have been taken from published papers, reports or web sites. Generation technologies are compared at scales where they have proved economic burning other fuels, rather than at a given size. A pyrolysis yield model was developed for a bubbling fluidised bed fast pyrolysis reactor from published data to predict bio-oil yields and pyrolysis plant energy demands. Generation using diesel engines, gas turbines in open and combined cycle (CCGT) operation and steam cycle plants was considered. The use of bio-oil storage to allow the pyrolysis and generation plants to operate independently of each other was investigated. The option of using diesel generators and open cycle gas turbines for combined heat and power was examined. The possible cost reductions that could be expected through learning if the technology is widely implemented were considered. It was found that none of the systems analysed would be viable without subsidy, but with the current Renewable Obligation Scheme CCGT plants in the 200 to 350 MWe range, super-critical coal fired boilers co-fired with bio-oil, and groups of diesel engine based CHP schemes supplied by a central pyrolysis plant would be viable. It was found that the cost would reduce with implementation and the planting of more energy crops but some subsidy would still be needed to make the plants viable.

Heat and mass transfer in membrane distillation used for desalination with slip flow

A theoretical model for the transport phenomena in an air gap membrane distillation is presented. The model is based on the conservation equations for the mass, momentum, energy and species within the feed water solution as well as on the mass and energy balances on the membrane sides. The slip flow occurs due to the hydrophobic properties of the membrane. The slip boundary condition applied on the feed saline solution-membrane interface is taken into consideration showing its effects on process parameters particularly permeate flow, heat transfer coefficient and thermal efficiency. The theoretical model was validated with available experimental data and was found to be in good agreement especially when the slip condition is introduced. Increasing slip length from zero to 200 μm was found to increase the permeate flux and the thermal efficiency by 33% and 1.7% respectively.

Cleavage initiation in the intercritically reheated coarse-grained heat affected zone:Part II. Failure criteria and statistical effects

In part 1 of this article, cleavage initiation in the intercritically reheated coarse-grained heat affected zone (IC CG HAZ) of high-strength low-alloy (HSLA) steels was determined to occur between two closely spaced blocky MA particles. Blunt notch, crack tip opening displacement (CTOD), and precracked Charpy testing were used in this investigation to determine the failure criteria required for cleavage initiation to occur by this mechanism in the IC CG HAZ. It was found that the attainment of a critical level of strain was required in addition to a critical level of stress. This does not occur in the case of high strain rate testing, for example, during precracked Charpy testing. A different cleavage initiation mechanism is then found to operate. The precise fracture criteria and microstructural requirements (described in part I of this article) result in competition between potential cleavage initiation mechanisms in the IC CG HAZ.

Controls on sensible heat and latent energy fluxes from a short-hydroperiod Florida Everglades marsh

Little is known of energy balance in low latitude wetlands where there is a year-round growing season and a climate best defined by wet and dry seasons. The Florida Everglades is a highly managed and extensive subtropical wetland that exerts a substantial influence on the hydrology and climate of the south Florida region. However, the effects of seasonality and active water management on energy balance in the Everglades ecosystem are poorly understood. An eddy covariance and micrometeorological tower was established in a short-hydroperiod Everglades marsh to examine the dominant environmental controls on sensible heat (H) and latent energy (LE) fluxes, as well as the effects of seasonality on these parameters. Seasonality differentially affected H and LE fluxes in this marsh, such that H was principally dominant in the dry season and LE was strongly dominant in the wet season. The Bowen ratio was high for much of the dry season (1.5–2.4), but relatively low (H and LE fluxes across nearly all seasons and years (). However, the 2009 dry season LE data were not consistent with this relationship () because of low seasonal variation in LE following a prolonged end to the previous wet season. In addition to net radiation, H and LE fluxes were significantly related to soil volumetric water content (VWC), water depth, air temperature, and occasionally vapor pressure deficit. Given that VWC and water depth were determined in part by water management decisions, it is clear that human actions have the ability to influence the mode of energy dissipation from this ecosystem. Impending modifications to water management under the Comprehensive Everglades Restoration Plan may shift the dominant turbulent flux from this ecosystem further toward LE, and this change will likely affect local hydrology and climate.

Developing an enhanced model for combined heat and air infiltration energy simulation

The need for efficient, sustainable, and planned utilization of resources is ever more critical. In the U.S. alone, buildings consume 34.8 Quadrillion (1015) BTU of energy annually at a cost of $1.4 Trillion. Of this energy 58% is utilized for heating and air conditioning. ^ Several building energy analysis tools have been developed to assess energy demands and lifecycle energy costs in buildings. Such analyses are also essential for an efficient HVAC design that overcomes the pitfalls of an under/over-designed system. DOE-2 is among the most widely known full building energy analysis models. It also constitutes the simulation engine of other prominent software such as eQUEST, EnergyPro, PowerDOE. Therefore, it is essential that DOE-2 energy simulations be characterized by high accuracy. ^ Infiltration is an uncontrolled process through which outside air leaks into a building. Studies have estimated infiltration to account for up to 50% of a building's energy demand. This, considered alongside the annual cost of buildings energy consumption, reveals the costs of air infiltration. It also stresses the need that prominent building energy simulation engines accurately account for its impact. ^ In this research the relative accuracy of current air infiltration calculation methods is evaluated against an intricate Multiphysics Hygrothermal CFD building envelope analysis. The full-scale CFD analysis is based on a meticulous representation of cracking in building envelopes and on real-life conditions. The research found that even the most advanced current infiltration methods, including in DOE-2, are at up to 96.13% relative error versus CFD analysis. ^ An Enhanced Model for Combined Heat and Air Infiltration Simulation was developed. The model resulted in 91.6% improvement in relative accuracy over current models. It reduces error versus CFD analysis to less than 4.5% while requiring less than 1% of the time required for such a complex hygrothermal analysis. The algorithm used in our model was demonstrated to be easy to integrate into DOE-2 and other engines as a standalone method for evaluating infiltration heat loads. This will vastly increase the accuracy of such simulation engines while maintaining their speed and ease of use characteristics that make them very widely used in building design.^

Under-ice export flux measurements by short-term drifting sediment traps at 27 stations in the Arctic Ocean during summer 1995, 1997, and 2012

A critical question regarding the organic carbon cycle in the Arctic Ocean is whether the decline in ice extent and thickness and the associated increase in solar irradiance in the upper ocean will result in increased primary production and particulate organic carbon (POC) export. To assess spatial and temporal variability in POC export, under-ice export fluxes were measured with short-term sediment traps in the northern Laptev Sea in July-August-September 1995, north of the Fram Strait in July 1997, and in the Central Arctic in August-September 2012. Sediment traps were deployed at 2-5 m and 20-25 m under ice for periods ranging from 8.5 to 71 h. In addition to POC fluxes, total particulate matter, chlorophyll a, biogenic particulate silica, phytoplankton, and zooplankton fecal pellet fluxes were measured to evaluate the amount and composition of the material exported in the upper Arctic Ocean. Whereas elevated export fluxes observed on and near the Laptev Sea shelf were likely the combined result of high primary production, resuspension, and release of particulate matter from melting ice, low export fluxes above the central basins despite increased light availability during the record minimum ice extent of 2012 suggest that POC export was limited by nutrient supply during summer. These results suggest that the ongoing decline in ice cover affects export fluxes differently on Arctic shelves and over the deep Arctic Ocean and that POC export is likely to remain low above the central basins unless additional nutrients are supplied to surface waters.

Convective Cold Pools over the Atlas Mountains and Their Influence on the Saharan Heat Low

The West African Monsoon (WAM) and its representation in numerical models are strongly influenced by the Saharan Heat Low (SHL), a low-pressure system driven by radiative heating over the central Sahara and ventilated by the cold and moist inflow from adjacent oceans. It has recently been shown that a significant part of the southerly moisture flux into the SHL originates from convective cold pools over the Sahel. These density currents driven by evaporation of rain are largely absent in models with parameterized convection. This crucial issue has been hypothesized to contribute to the inability of many climate models to reproduce the variability of the WAM. Here, the role of convective cold pools approaching the SHL from the Atlas Mountains, which are a strong orographic trigger for deep convection in Northwest Africa, is analyzed. Knowledge about the frequency of these events, as well as their impact on large-scale dynamics, is required to understand their contribution to the variability of the SHL and to known model uncertainties. The first aspect is addressed through the development of an objective and automated method for the generation of multi-year climatologies not available before. The algorithm combines freely available standard surface observations with satellite microwave data. Representativeness of stations and influence of their spatial density are addressed by comparison to a satellite-only climatology. Applying this algorithm to data from automated weather stations and manned synoptic stations in and south of the Atlas Mountains reveals the frequent occurrence. On the order of 6 events per month are detected from May to September when the SHL is in its northernmost position. The events tend to cluster into several-days long convectively active periods, often with strong events on consecutive days. This study is the first to diagnose dynamical impacts of such periods on the SHL, based on simulations of two example cases using the Weather Research and Forecast (WRF) model at convection-permitting resolution. Sensitivity experiments with artificially removed cold pools as well as different resolutions and parameterizations are conducted. Results indicate increases in surface pressure of more than 1 hPa and significant moisture transports into the desert over several days. This moisture affects radiative heating and thus the energy balance of the SHL. Even though cold pool events north of the SHL are less frequent when compared to their Sahelian counterparts, it is shown that they gain importance due to their temporal clustering on synoptic timescale. Together with studies focusing on the Sahel, this work emphasizes the need for improved parameterization schemes for deep convection in order to produce more reliable climate projections for the WAM.

DIRECT NUMERICAL SIMULATION OF INCOMPRESSIBLE MULTIPHASE FLOW WITH PHASE CHANGE

Phase change problems arise in many practical applications such as air-conditioning and refrigeration, thermal energy storage systems and thermal management of electronic devices. The physical phenomenon in such applications are complex and are often difficult to be studied in detail with the help of only experimental techniques. The efforts to improve computational techniques for analyzing two-phase flow problems with phase change are therefore gaining momentum. The development of numerical methods for multiphase flow has been motivated generally by the need to account more accurately for (a) large topological changes such as phase breakup and merging, (b) sharp representation of the interface and its discontinuous properties and (c) accurate and mass conserving motion of the interface. In addition to these considerations, numerical simulation of multiphase flow with phase change introduces additional challenges related to discontinuities in the velocity and the temperature fields. Moreover, the velocity field is no longer divergence free. For phase change problems, the focus of developmental efforts has thus been on numerically attaining a proper conservation of energy across the interface in addition to the accurate treatment of fluxes of mass and momentum conservation as well as the associated interface advection. Among the initial efforts related to the simulation of bubble growth in film boiling applications the work in \cite{Welch1995} was based on the interface tracking method using a moving unstructured mesh. That study considered moderate interfacial deformations. A similar problem was subsequently studied using moving, boundary fitted grids \cite{Son1997}, again for regimes of relatively small topological changes. A hybrid interface tracking method with a moving interface grid overlapping a static Eulerian grid was developed \cite{Juric1998} for the computation of a range of phase change problems including, three-dimensional film boiling \cite{esmaeeli2004computations}, multimode two-dimensional pool boiling \cite{Esmaeeli2004} and film boiling on horizontal cylinders \cite{Esmaeeli2004a}. The handling of interface merging and pinch off however remains a challenge with methods that explicitly track the interface. As large topological changes are crucial for phase change problems, attention has turned in recent years to front capturing methods utilizing implicit interfaces that are more effective in treating complex interface deformations. The VOF (Volume of Fluid) method was adopted in \cite{Welch2000} to simulate the one-dimensional Stefan problem and the two-dimensional film boiling problem. The approach employed a specific model for mass transfer across the interface involving a mass source term within cells containing the interface. This VOF based approach was further coupled with the level set method in \cite{Son1998}, employing a smeared-out Heaviside function to avoid the numerical instability related to the source term. The coupled level set, volume of fluid method and the diffused interface approach was used for film boiling with water and R134a at the near critical pressure condition \cite{Tomar2005}. The effect of superheat and saturation pressure on the frequency of bubble formation were analyzed with this approach. The work in \cite{Gibou2007} used the ghost fluid and the level set methods for phase change simulations. A similar approach was adopted in \cite{Son2008} to study various boiling problems including three-dimensional film boiling on a horizontal cylinder, nucleate boiling in microcavity \cite{lee2010numerical} and flow boiling in a finned microchannel \cite{lee2012direct}. The work in \cite{tanguy2007level} also used the ghost fluid method and proposed an improved algorithm based on enforcing continuity and divergence-free condition for the extended velocity field. The work in \cite{sato2013sharp} employed a multiphase model based on volume fraction with interface sharpening scheme and derived a phase change model based on local interface area and mass flux. Among the front capturing methods, sharp interface methods have been found to be particularly effective both for implementing sharp jumps and for resolving the interfacial velocity field. However, sharp velocity jumps render the solution susceptible to erroneous oscillations in pressure and also lead to spurious interface velocities. To implement phase change, the work in \cite{Hardt2008} employed point mass source terms derived from a physical basis for the evaporating mass flux. To avoid numerical instability, the authors smeared the mass source by solving a pseudo time-step diffusion equation. This measure however led to mass conservation issues due to non-symmetric integration over the distributed mass source region. The problem of spurious pressure oscillations related to point mass sources was also investigated by \cite{Schlottke2008}. Although their method is based on the VOF, the large pressure peaks associated with sharp mass source was observed to be similar to that for the interface tracking method. Such spurious fluctuation in pressure are essentially undesirable because the effect is globally transmitted in incompressible flow. Hence, the pressure field formation due to phase change need to be implemented with greater accuracy than is reported in current literature. The accuracy of interface advection in the presence of interfacial mass flux (mass flux conservation) has been discussed in \cite{tanguy2007level,tanguy2014benchmarks}. The authors found that the method of extending one phase velocity to entire domain suggested by Nguyen et al. in \cite{nguyen2001boundary} suffers from a lack of mass flux conservation when the density difference is high. To improve the solution, the authors impose a divergence-free condition for the extended velocity field by solving a constant coefficient Poisson equation. The approach has shown good results with enclosed bubble or droplet but is not general for more complex flow and requires additional solution of the linear system of equations. In current thesis, an improved approach that addresses both the numerical oscillation of pressure and the spurious interface velocity field is presented by featuring (i) continuous velocity and density fields within a thin interfacial region and (ii) temporal velocity correction steps to avoid unphysical pressure source term. Also I propose a general (iii) mass flux projection correction for improved mass flux conservation. The pressure and the temperature gradient jump condition are treated sharply. A series of one-dimensional and two-dimensional problems are solved to verify the performance of the new algorithm. Two-dimensional and cylindrical film boiling problems are also demonstrated and show good qualitative agreement with the experimental observations and heat transfer correlations. Finally, a study on Taylor bubble flow with heat transfer and phase change in a small vertical tube in axisymmetric coordinates is carried out using the new multiphase, phase change method.

COMPREHENSIVE ELECTRICAL/OPTICAL/THERMAL CHARACTERIZATIONS OF HIGH POWER LIGHT EMITTING DIODES AND LASER DIODES

Thermal characterizations of high power light emitting diodes (LEDs) and laser diodes (LDs) are one of the most critical issues to achieve optimal performance such as center wavelength, spectrum, power efficiency, and reliability. Unique electrical/optical/thermal characterizations are proposed to analyze the complex thermal issues of high power LEDs and LDs. First, an advanced inverse approach, based on the transient junction temperature behavior, is proposed and implemented to quantify the resistance of the die-attach thermal interface (DTI) in high power LEDs. A hybrid analytical/numerical model is utilized to determine an approximate transient junction temperature behavior, which is governed predominantly by the resistance of the DTI. Then, an accurate value of the resistance of the DTI is determined inversely from the experimental data over the predetermined transient time domain using numerical modeling. Secondly, the effect of junction temperature on heat dissipation of high power LEDs is investigated. The theoretical aspect of junction temperature dependency of two major parameters – the forward voltage and the radiant flux – on heat dissipation is reviewed. Actual measurements of the heat dissipation over a wide range of junction temperatures are followed to quantify the effect of the parameters using commercially available LEDs. An empirical model of heat dissipation is proposed for applications in practice. Finally, a hybrid experimental/numerical method is proposed to predict the junction temperature distribution of a high power LD bar. A commercial water-cooled LD bar is used to present the proposed method. A unique experimental setup is developed and implemented to measure the average junction temperatures of the LD bar. After measuring the heat dissipation of the LD bar, the effective heat transfer coefficient of the cooling system is determined inversely. The characterized properties are used to predict the junction temperature distribution over the LD bar under high operating currents. The results are presented in conjunction with the wall-plug efficiency and the center wavelength shift.

Critical core mass for enriched envelopes: the role of H2O condensation

Context. Within the core accretion scenario of planetary formation, most simulations performed so far always assume the accreting envelope to have a solar composition. From the study of meteorite showers on Earth and numerical simulations, we know that planetesimals must undergo thermal ablation and disruption when crossing a protoplanetary envelope. Thus, once the protoplanet has acquired an atmosphere, not all planetesimals reach the core intact, i.e. the primordial envelope (mainly H and He) gets enriched in volatiles and silicates from the planetesimals. This change of envelope composition during the formation can have a significant effect on the final atmospheric composition and on the formation timescale of giant planets. Aims. We investigate the physical implications of considering the envelope enrichment of protoplanets due to the disruption of icy planetesimals during their way to the core. Particular focus is placed on the effect on the critical core mass for envelopes where condensation of water can occur. Methods. Internal structure models are numerically solved with the implementation of updated opacities for all ranges of metallicities and the software Chemical Equilibrium with Applications to compute the equation of state. This package computes the chemical equilibrium for an arbitrary mixture of gases and allows the condensation of some species, including water. This means that the latent heat of phase transitions is consistently incorporated in the total energy budget. Results. The critical core mass is found to decrease significantly when an enriched envelope composition is considered in the internal structure equations. A particularly strong reduction of the critical core mass is obtained for planets whose envelope metallicity is larger than Z approximate to 0.45 when the outer boundary conditions are suitable for condensation of water to occur in the top layers of the atmosphere. We show that this effect is qualitatively preserved even when the atmosphere is out of chemical equilibrium. Conclusions. Our results indicate that the effect of water condensation in the envelope of protoplanets can severely affect the critical core mass, and should be considered in future studies.

Heat fluxes between the Guadalquivir river and the Gulf of Cádiz Continental Shelf

An 18-year time series of daily sea surface temperature of Gulf of Cadiz and an 18-month time series of temperature collected in the vicinity of the Guadalquivir estuary mouth have been analyzed to investigate the heat exchange between the estuary and the adjacent continental shelf. The first time identifies a continental shelf area where seasonal thermal oscillation signal (amplitudes and phase) changes abruptly. In order to explain this anomaly, the second data set allows a description of thermal fluctuations in a wide range of frequencies and an estimation of the upstream heat budget of the Guadalquivir estuary. Results show that high frequency thermal signal, diurnal and semidiurnal, and water flux signal through Guadalquivir mouth, mainly semidiurnal, apparently interact randomly to give a small exchange of thermal energy at high frequency. There is no trace, at the estuary's mouth, of daily heat exchanges with intertidal mudflats probably because it tends to cancel on daily time scales. Results also show that fluctuations of estimated air-sea fluxes force fluctuations of temperature in a quite homogeneous estuarine, with a delay of 20 days. The along-channel thermal energy gradient reaches magnitudes of 300-400 J m-4 near the mouth during the summer and winter and drives the estuary-shelf exchange of thermal energy at seasonal scale. Particularly, the thermal heat imported by the estuary from the shelf area during late fall-winter-early spring of 2008/2009 is balanced by the thermal heat that the estuary exports to the shelf area during late spring-summer of 2008. In summary, Guadalquivir river removes/imports excess of thermal energy towards/from the continental shelf seasonally, as a mechanism to accommodate excess of heat from one side respect to the other side.

Treatment of textile wastewater with membrane bioreactor: a critical review

Membrane bioreactor (MBR) technology has been used widely for various industrial wastewater treatments due to its distinct advantages over conventional bioreactors. Treatment of textile wastewater using MBR has been investigated as a simple, reliable and cost-effective process with a significant removal of contaminants. However, a major drawback in the operation of MBR is membrane fouling, which leads to the decline in permeate flux and therefore requires membrane cleaning. This eventually decreases the lifespan of the membrane. In this paper, the application of aerobic and anaerobic MBR for textile wastewater treatment as well as fouling and control of fouling in MBR processes have been reviewed. It has been found that long sludge retention time increases the degradation of pollutants by allowing slow growing microorganisms to establish but also contributes to membrane fouling. Further research aspects of MBR for textile wastewater treatment are also considered for sustainable operations of the process.

A survey on mechanisms and critical parameters on solidification of selective laser melting during fabrication of Ti-6Al-4V prosthetic acetabular cup

Additive Manufacturing (AM) includes a range of approaches that correlate with computer aided design (CAD) and manufacturing by fabrication via precise layers and is a promising method for the production of medical tools. In this study, different aspects and mechanisms of solidification for curved surfaces based on equilibrium at curved interfaces, Monge patch, interfacial and Gibbs energy will be discussed. Also, the effect of capillarity, geometry, substrate temperature, cooling rate and scanning parameters in the solidification of a prosthetic acetabular cup (PAC) using selective laser melting (SLM) is analysed. The contributions of this work are analysing solidification and effective factors in this process to produce parts with a higher quality and mechanical properties such as strength, strain, porosity, relative density and hardness. Results indicate that due to the surface to volume (S/V) ratio, and the increasing effect of the radius on Monge patch, thermal stresses and surface forces are more prevalent on outer surfaces. Moreover, solidification and mechanical properties are related to capillarity, geometry, substrate temperature, cooling rate, scanning power and speed. The results also indicate the interaction of solute diffusion and heat transfer with interatomic forces in large S/V ratio and at small scales tend to improve solidification.

Developing an Enhanced Model for Combined Heat and Air Infiltration Energy Simulation

The need for efficient, sustainable, and planned utilization of resources is ever more critical. In the U.S. alone, buildings consume 34.8 Quadrillion (1015) BTU of energy annually at a cost of $1.4 Trillion. Of this energy 58% is utilized for heating and air conditioning. Several building energy analysis tools have been developed to assess energy demands and lifecycle energy costs in buildings. Such analyses are also essential for an efficient HVAC design that overcomes the pitfalls of an under/over-designed system. DOE-2 is among the most widely known full building energy analysis models. It also constitutes the simulation engine of other prominent software such as eQUEST, EnergyPro, PowerDOE. Therefore, it is essential that DOE-2 energy simulations be characterized by high accuracy. Infiltration is an uncontrolled process through which outside air leaks into a building. Studies have estimated infiltration to account for up to 50% of a building’s energy demand. This, considered alongside the annual cost of buildings energy consumption, reveals the costs of air infiltration. It also stresses the need that prominent building energy simulation engines accurately account for its impact. In this research the relative accuracy of current air infiltration calculation methods is evaluated against an intricate Multiphysics Hygrothermal CFD building envelope analysis. The full-scale CFD analysis is based on a meticulous representation of cracking in building envelopes and on real-life conditions. The research found that even the most advanced current infiltration methods, including in DOE-2, are at up to 96.13% relative error versus CFD analysis. An Enhanced Model for Combined Heat and Air Infiltration Simulation was developed. The model resulted in 91.6% improvement in relative accuracy over current models. It reduces error versus CFD analysis to less than 4.5% while requiring less than 1% of the time required for such a complex hygrothermal analysis. The algorithm used in our model was demonstrated to be easy to integrate into DOE-2 and other engines as a standalone method for evaluating infiltration heat loads. This will vastly increase the accuracy of such simulation engines while maintaining their speed and ease of use characteristics that make them very widely used in building design.