Next Article in Journal
Sustainable Integration of a Solar Heating System into a Single-Family House in the Climate of Central Europe—A Case Study
Next Article in Special Issue
Analyzing Atrium Volume Designs for Hot and Humid Climates
Previous Article in Journal
Strength Development Characteristics of SBR-Modified Cementitious Mixtures for 3-Demensional Concrete Printing
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

A Review of CFD Analysis Methods for Personalized Ventilation (PV) in Indoor Built Environments

1
School of Thermal Engineering, Shandong Jianzhu University, Jinan 250101, China
2
Department of Mechanical Engineering, University of Maryland, College Park, MD 20742, USA
3
Department of Architecture, Xi’an Jiaotong-Liverpool University, Suzhou 215123, China
*
Authors to whom correspondence should be addressed.
Sustainability 2019, 11(15), 4166; https://doi.org/10.3390/su11154166
Submission received: 7 July 2019 / Revised: 26 July 2019 / Accepted: 28 July 2019 / Published: 1 August 2019

Abstract

:
Computational fluid dynamics (CFD) is an effective analysis method of personalized ventilation (PV) in indoor built environments. As an increasingly important supplement to experimental and theoretical methods, the quality of CFD simulations must be maintained through an adequately controlled numerical modeling process. CFD numerical data can explain PV performance in terms of inhaled air quality, occupants’ thermal comfort, and building energy savings. Therefore, this paper presents state-of-the-art CFD analyses of PV systems in indoor built environments. The results emphasize the importance of accurate thermal boundary conditions for computational thermal manikins (CTMs) to properly analyze the heat exchange between human body and the microenvironment, including both convective and radiative heat exchange. CFD modeling performance is examined in terms of effectiveness of computational grids, convergence criteria, and validation methods. Additionally, indices of PV performance are suggested as system-performance evaluation criteria. A specific utilization of realistic PV air supply diffuser configurations remains a challenging task for further study. Overall, the adaptable airflow characteristics of a PV air supply provide an opportunity to achieve better thermal comfort with lower energy use based on CFD numerical analyses.

1. Introduction

Conventional heating, ventilation, and air conditioning (HVAC) systems are used to maintain a uniform indoor environment in compliance with design standards. However, a certain population of occupants is more sensitive to local air movement, temperature, and relative humidity due to individual preferences as well as differences in clothing and activity levels. Therefore, thermal discomfort exists for a percentage of occupants even when the thermal environment complies with recommendations in the standards [1,2,3]. Moreover, symptoms of sick-building syndrome have resulted in increased attention to inhaled air quality based on individual needs and preferences [4]. Personalized ventilation (PV) aims to supply clean air directly to the breathing zone of each occupant to improve air quality. Several recent studies have introduced PV as a viable ventilation concept for office environments, and not just vehicles [3,5,6,7]. A paradigm shift from the design of indoor uniform environments to the design of individually controlled environments is advantageous even with respect to energy consumption in buildings. Several recent studies have found that individually controlled environments combined with a wider thermostat setpoint range for the background indoor environment result in greater thermal comfort and enable building energy savings [8,9,10,11].
The original PV studies for indoor built environments, such as offices and hospitals, (i) focused on delivering a small amount of personalized clean air directly to the breathing zone of each occupant [3,12]; (ii) improved thermal comfort [13,14]; and (iii) improved energy efficiency [11,15,16]. Numerous laboratory studies have demonstrated PV performance in different buildings and climate zones, including hot and humid climates [15,17,18] and cold climates [19]. PV is typically used in conjunction with total volume ventilation [10,20], such as mixing ventilation (MV) [21,22], displacement ventilation (DV) [16,23], and underfloor air distribution (UFAD) systems [24,25,26], or combined with liquid-based cooling systems, such as chilled ceilings [27,28]. PV is also considered as a standalone system with no assistive ventilation or cooling systems. For instance, a chair-based PV has been proposed for theatres, cinemas, and lecture halls [29]. The performance of PV on improving inhaled air quality, thermal comfort, and energy use has been widely studied.
In addition to experimental methods, computational fluid dynamics (CFD), as a simulation technique that can efficiently predict both spatial and temporal fields in the research of ventilation and indoor air science, has been adopted and studied [30,31,32]. Despite the dramatic increase in its application, CFD has not replaced experimental and theoretical analyses, but it has become increasingly important as a supplement to them [33,34,35]. In the PV field, a considerable number of studies have used CFD to analyze indoor air quality, thermal comfort, ventilation effectiveness, and energy saving in association with turbulence modeling, numerical approximations, and boundary conditions [36,37,38,39]. Compared to experimental methods, CFD can provide precise information regarding the distribution of flow and concentration fields in the whole simulation domain, rather than just targeted areas for data collection [40,41,42]. However, the quality of CFD simulation requires adequate grid generation, selection of numerical models, assessment of numerical errors, and control over other numerical parameters for proper verification and validation [43,44,45,46]. Therefore, as a system suitable for CFD analysis, PV requires a careful consideration of numerical parameters.
In recent years, several literature-review papers have been published in the PV research field. These reviews focused on different ventilation types supplemented by PV [47], personalized conditioning influencing thermal comfort and energy consumption [48], quantification of the ability of personal comfort systems to produce comfort [49], and advanced personal comfort systems in the workplace [50]. A number of review papers regarding CFD techniques in indoor built environments have been also published, including a CFD analysis procedure [51], quality control [52], application aspects and trends [53], turbulence models [54,55], verification and validation [56], pressure–velocity decoupling algorithms [57], CFD coupled with network modeling [58], large eddy simulation (LES), and Reynolds-averaged Navier–Stokes (RANS) approaches in building simulations [59], as well as reviews of CFD and indoor ventilation [30,33,60,61]. However, previous review studies have not conducted a PV review focusing on the CFD modeling approach, hence this paper targets this important knowledge base. The most relevant papers include a review of CFD studies to analyze and design the microenvironment around the human body [62], and a review of CFD studies of the thermal environment around a human body accounting for the PV system [63]. Nevertheless, with the increasingly fast development of the CFD technique and PV systems, a timely, critical literature review of CFD applications to PV systems is needed. Most importantly, there has been no comprehensive overview with respect to PV and its effect on the microenvironment around the human body.
This study reviews the existing knowledge base of CFD analyses focused on PV systems in indoor built environments. This review first provides an overview of CFD applications for PV in indoor built environments, starting with a brief introduction of PV research topics, CFD codes, PV types, and their development benefits from the perspective of CFD technical advancements. CFD setups of PV simulations are then reviewed, including turbulence models and model performance. This study specifically emphasizes the heat transfer around the human body. PV systems’ parameters and their performance are analyzed, with a focus on PV air supply diffusers, PV air supply parameters, and evaluation indices. Through state-of-the-art analyses of current CFD applications in studying PV systems, future CFD procedures and research interests are suggested to provide a preliminary guideline of PV analyses for HVAC engineers and researchers.

2. Research Methodology

A series of scientific databases—including ScienceDirect, SpringerLink, Taylor and Francis Online, Wiley Online Library, SAGE journals, and Web of Science, as well as the Google Scholar search engine—were employed for this scientific literature review. In combination with “CFD” or “numerical simulation”, the following search keywords were used: “personalized ventilation”, “task ambient conditioning”, “air terminal diffuser”, “task ambient air conditioning”, “personalized conditioning”, “personal thermal comfort”, “personal cooling”, and “personal heating”. The searched articles were examined as they related to the aim of this review paper. Moreover, relevant literature known to the authors or reviewed through the above references was included. Additional studies that examine the capacity and accuracy of the CFD technique in simulating indoor built environments were also included.
A total 60 journal articles that studied PV’s application in indoor built environments by CFD methods were reviewed in this study. Moreover, this review covered the reviews for CFD studies on indoor ventilation, which referred to many aspects of CFD simulation for indoor built environments, therefore were considered beneficial to this study [30,33,60,61]. Some important studies on CTM development were also included in the review, even they did not involve PV systems. However, those CFD studies on PV’s application in vehicles, such as aircraft [64,65] and vehicle compartments [66,67], are excluded in this review as they are out of the scope.

3. Overview of Previous CFD Applications in PV in Indoor Built Environments

Compared to CFD studies of indoor ventilation, a comprehensive PV study using CFD simulations must account not only just for the geometric model, grid generation, turbulence model, boundary condition, numerical scheme, and post-processing, but also for the computational thermal manikin (CTM). Airflow interactions with the CTM surface significantly affect inhaled indoor air quality and thermal comfort induced by convective human body plume or forced jet flow [68]. Therefore, CTM in PV studies is reviewed. Sixty journal articles focus on PV and CFD. Table 1 shows an overview of CFD analysis of PV, with the following categories: (1) author(s) and publication year, (2) PV type, (3) CFD code, (4) turbulence model, (5) radiation model, (6) CTM surface boundary condition, (7) near wall treatment, and (8) target parameters. The primary overview covers the yearly publication distribution, CFD codes, research topics, and PV types with background HVAC systems, to provide a development status of PV studies using CFD.

3.1. Yearly Publication Distribution

Published in 2004 [36], the first article in the scope of PV and CFD studied the microenvironment around the human body resulting from a PV system operation. Figure 1a shows the yearly distribution of CFD studies of PV systems, which are increasingly popular. It is found that the number of CFD studies of PV in the last three years already exceeds 45% of all of the studies.

3.2. CFD Codes

Figure 1b shows the CFD code distribution of CFD studies of PV systems. The figure shows that of all of the publications reviewed in this study, 38 of 60 articles (63%) employed ANSYS Fluent software (including the early version, called Fluent). Other commonly used programs are STAR-CCM+, Airpak, PHOENICS, STACH-3, and ANSYS. A study compared ANSYS Fluent and STAR-CCM+, which are currently the most popular, in predicting indoor air parameter distributions [125]. The study pointed out that the two programs produced almost the same results with similar computing effort, and emphasized that ANSYS Fluent was more user-friendly, owing to its user-defined functions. Authors do not promote any commercial CFD package, and encourage diverse application of CFD in PV studies, including the open-source OpenFOAM, which has powerful grid capabilities and a wide range of applications [126]. Note that Airpak is no longer updated, but it has been used in recent studies [123].

3.3. Research Topics of PV Studies

Research topics of PV using CFD are categorized into three major types: (i) local ventilation and inhaled air quality, (ii) thermal comfort, and (iii) energy savings, based on the target parameters summarized in Table 1. Local ventilation and inhaled air quality are evaluated based on the air change rate (h−1) (ACH), carbon dioxide concentration (ppm) (CO2), concentration asymmetry (-) (ΔC), deposited fraction (DFr), inhaled contaminant concentration (ppm) (Cinh), intake fraction (-) (IF), particle concentration (kg/m3) (PC), pollutant exposure reduction (%) (PER), ventilation effectiveness (-) (ES), personal exposure effectiveness (-) (PEE), scale for ventilation efficiency 3 (-) (SVE3), new scale for ventilation efficiency 3 (-) (SVE3*), scale for ventilation efficiency 4 (-) (SVE4), re-inhaled exposure index (-) (RIE), and tracer gas concentration (ppm) (TGC). Thermal comfort studies use local thermal comfort (-) (LTC), local thermal sensation (-) (LTS), overall thermal comfort (-) (OTC), overall thermal sensation (-) (OTS), predicted mean vote (-) (PMV), predicted percentage dissatisfied (%) (PPD), draft risk (%) (DR), CTM-based equivalent temperature (°C) (Teq), operative temperature (°C) (Top), vertical air temperature difference at the head and foot level (°C) (VATD), and mean radiant temperature (°C) (MRT). Energy saving studies evaluate the cooling efficiency (-) (CE), energy utilization coefficient (-) (EUC), and calculated energy consumption (W) (Qc).
It is important to note that the target parameters used in Table 1 are not the recommended evaluation indices, but instead are the main parameters used in literature focused on validation, comparison, or performance evaluation PV systems with CFD. Being consistent with the original aim of studying PV systems, the research topic category of local ventilation and air quality constitutes almost 80% (46 of 60) of the reviewed studies. Thermal comfort and energy saving were investigated by 28 studies and 18 studies, respectively. Additionally, in the last three years, the three major types of research topic—local ventilation and air quality, thermal comfort, and energy saving—have the proportions of 20:18:16 in a total of 27 studies, which shows almost equal importance. This distribution coincides with two of the current challenges in building science: (i) the growing need for comfort improvements and (ii) increased concern for energy consumption [48]. The commonly used evaluation indices will be described in the next section.

3.4. PV Types with Background HVAC Systems

A comprehensive study presented five air terminal devices (ATDs), namely, a movable panel (MP), computer monitor panel (CMP), vertical desk grill (VDG), horizontal desk grill (HDG), and personal environment module (PEM) [6]. Another study redesigned a round movable panel (RMP) [127], used in several review studies [38,80,92]. The commonly used PV type also includes ceiling-mounted personalized ventilation (CMPV) [88,89], chair-based personalized ventilation (CBPV) [29,37], desk-mounted personalized ventilation (DMPV) [120], desk fan (DF) and/or chair fan (CF) [103,104], and bed-based task/ambient air conditioning (BTAC) [105,109,110]. Figure 2 shows the commonly used PV types. Note that several PV types included in relatively few articles are not presented in this figure and can be found in the references. These less popular PV types include ductless personalized ventilation (DPV) [118], local exhaust ventilation for offices (LEVO) [111,112], movable slot diffuser plate (MSDP) [82], natural personalized ventilation (NPV) [83], robotic personal air conditioning (RPAC) [107], targeted personalized ventilation (TPV) [122], and wide-cover personalized ventilation (WCPV) [101].
Depending on individual preferences, especially in large offices with several workstations, occupants may adjust their PV system supply flow rate or temperature, or even switch it off. Therefore, the background environment must maintain an acceptable comfort level. The background HVAC system can be used in combination with PV, resulting in a whole space with non-uniform velocity and temperature fields and unevenly distributed pollutants [3]. Background HVAC systems are classified into two types: (i) background ventilation and (ii) radiant cooling/heating. The former is a convection-dominated system, i.e., total volume system ventilation [10], including MV (14 studies), DV (14 studies), and UFAD (8 studies). The latter is radiation-dominated and uses radiant panels [114,128] or a chilled ceiling [116,119]. It should be noted that some PV systems recirculate air in addition to using outside air by utilizing a low-mixing PV nozzle integrated with a four-sided slot diffuser [89,93], as shown in Figure 2c. The supply PV nozzle creates a jet that can penetrate the thermal plume around the human body and reach an occupant’s breathing zone with the necessary amount of outside air to maintain acceptable inhaled air quality.

4. Modeling Computational Thermal Manikin

In reviewed PV studies, modeling of CTMs plays a significant role in determining the distributions of airflow temperature, velocity, and pollutant concentrations. These distributions directly affect inhaled air quality, occupant thermal comfort, and energy consumption. The personal microenvironment is highly influenced by the convective heat transfer between the body surface and ambient environment, radiant heat transfer between the body surface and wall surface, and buoyancy-driven thermal plume around the body. Therefore, it is important to identify best CFD practices for implementation of realistic CTMs. This review explores related topics, including turbulence models, heat transfer around CTMs, and modeling performance.

4.1. Turbulence Models

PV studies primarily employ RANS methods. The three most commonly used RANS turbulence models are the renormalization group k-ε turbulence model (RNG) [129] (20 of 60 studies), realizable k-ε turbulence model (RKE) [130] (16 of 60 studies), and standard k-ε turbulence model (SKE) [131] (13 of 60 studies). Shear stress transport k-ω (SST) [132] was used in four studies, zero-equation turbulence model (ZEQ) [133] in three studies, and Reynolds stress models (RSMs) [134] in one study. Table 1 presents additional details on the turbulence models used in studies of PV systems.
For RANS-based models, selecting a proper turbulence model is a prerequisite for conducting performance evaluations of PV systems. In general, a PV system represents a mixed convection heat-transfer mode, including human body buoyancy-driven thermal plumes, natural convection induced by background ventilation, and forced convection due to PV nozzle supply airflow. Several studies have evaluated the performance of turbulence models. A study simulated a low-mixing co-flow nozzle directing a PV fresh-air jet toward the breathing zone of a seated occupant using SKE, RKE, RNG, the k-omega model (SKW), and SST turbulence models [75]. The results showed that the k-ε family of turbulence models presented similar results and agreed well with experimental data, while the SKW and SST turbulence models performed poorly at predicting jet deflection/distortion and jet spreading. For instance, for average convective heat flux, the k-ε family of turbulence models overestimated the values by up to 1%, while the k-ω family of turbulence models over-predicted the values by up to 4%. Similarly, the results of SKE and RNG are in good agreement with experimental data, while performance of the k-ω group models does not agree well with experiments in forced and mixed convection cases [135]. Another study validated a conventional round nozzle delivering fresh air directly to an occupant’s face using STK, RKE, and RNG [118]. RKE achieved the best results based on the air temperature, ventilation effectiveness index, and air velocity between the PV outlet and the occupant’s face. Another study evaluated the performance of the SKE, RNG, and SST models for predicting vertical profiles of air temperature and air velocity distribution around a sleeping occupant using bed-based task/ambient air conditioning; SST achieved a better agreement than the other models with measured data [105]. A study compared SKE, RKE, SST, and RSM for predicting the air velocity around the human head using a personalized square diffuser, and pointed out that RSM performed better than the other models [122].
Besides the reviewed comparisons for different turbulence models, most studies followed recommendations outlined in CFD program manuals and/or guidelines for CFD applications in indoor/enclosure environments, to choose appropriate turbulence models among a wide range of popular ones. For instance, among the RANS models, RNG was extensively used due to its stable performance in simulations of both forced and mixed convection around the human body [107], accuracy and good predictability of the jet spreading rate and the behavior of recirculated air [113], and better overall performance in terms of accuracy, computing efficiency, and robustness [81]. RKE is also popular because of its high accuracy in predicting the spreading rate of planar and round jets, and its superior predictions of flows involving separation and recirculation in indoor environments [88,103,117]. Moreover, the RKE model is effective for both high and low Reynolds number turbulent flows, and can therefore produce good-quality results even when the flow is not fully turbulent [74]. Finally, SKE is widely used for CFD simulations of PV airflows in enclosed environments because it can simulate convective heat transfer of buoyancy-driven airflows with reasonably acceptable results [36,79], despite its many deficiencies, such as inaccurate predictions of turbulence kinetic energy and adverse pressure gradient flows [71].
High Reynolds number turbulence models are almost always applied in PV studies due to the development of enhanced wall treatment (EWT, 25 of 54 studies) in ANSYS fluent and/or two-layer y+ wall treatment (TLWT, 2 of 54 studies) in STAR CCM+. The EWT approach uses a near-wall modeling method that combines a two-layer model with enhanced wall functions. If the near-wall mesh is fine enough (typically y+≈1), then EWT will automatically resolve the viscous sublayer. Additional details on the grid distribution adjacent to human body surfaces are provided in the following section. Furthermore, low Reynolds k-ε turbulence models (LRN) were typically used to account for convective heat transfer with the locally-fine grid adjacent to the body surface before applicable wall functions were adopted [62,136,137]. LRN can predict the transport process throughout the boundary layer with a large number of computational cells, in contrast to the standard k-ε model, which is valid only for fully developed turbulent flows [52]. However, LRN models converge with difficulty due to the complicated geometry associated with extremely fine grids near the human body [63], so their popularity is waning.
It should be noted that few studies utilize the modified v2f-dav model, detached eddy simulation model (DES), or the more robust but computationally demanding LES, which have been found to perform better for force convection and mixed convection in ventilated spaces [54,55]. A study concluded that v2f-dav and RNG showed the best overall performance compared to other models in terms of accuracy, computing efficiency, and robustness [138]. The v2f-dav model is not commonly applied in indoor environments, probably because it is unavailable in commercial CFD software. Another assessment of turbulence models for transitional flows in an enclosed environment [139], demonstrated that the RNG and RSM models could achieve a good overall performance to characterizing a flow before it became too complex, especially when these two models were used to simulate mixed convection in a room with a heated box. DES with sub-grid-scale turbulence kinetic energy and LES still has the most accurate and stable performance in enclosed environments [54,55]. LES provided the most detailed flow features and accuracy, but used significant computing time [61].
Although no single turbulence model can optimally and economically handle all flow elements—e.g., jet flow, momentum-driven flow, stratified flow, and buoyancy-driven flow [30]—there is still a great necessity to utilize more accurate turbulence models or predict the human body microenvironment characterized by a complex flow characteristic. In fact, LES has been applied to predict the microenvironment around a manikin, mostly using simplified geometry [140,141,142]. With increasing computer power, LES may gain popularity for PV studies. Further study is required to investigate the use of more advanced models, to improve the balance between CFD simulation accuracy and required computational effort.

4.2. Heat Transfer around the Human Body

Accurate predictions of heat transfer around the human body have direct implications on accuracy of predicted thermal comfort and inhaled air quality with PV systems. When introducing CTMs in CFD modeling of PV systems, it is important to consider CTM’s geometry, heat exchange models, boundary conditions, and thermoregulation models.

4.2.1. Geometry of CTMs

There are no relevant standards for the sizes, shapes, and postures of CTMs, and their usage mostly depends on the purposes of research studies. CTMs used a surface area of 1.57 m2 [36,92], 1.59 m2 [69], or 1.8 m2 [115], and a simple shape, e.g., rectangular, or exact-shape, as well as sleeping, seated, or standing postures. Most studies used complex CTM geometry, the exception being a rectangular person used to assess the performance of a CMPV system in regard to inhaled air quality influenced by a moving person [94]. This is mainly due to a complex flow field and the extensive computation required. Note that in Airpak software, the CTM geometry is predefined as a rectangular shape of some limited body parts and can only be changed according to the posture or height/length, because its library function stores or retrieves groups of objects in a parts library [143].
Computational requirements are mostly driven by the geometric complexity of the human body, resulting in different levels of findings adjacent to the body surface. Theoretically, the more realistic the CTM geometry the more accurate the results. Although existing literature has not examined the effect of CTM geometry simplification on PV performance, a few studies have qualified the influence of geometric complexity on an indoor environment, e.g., the shape and size of CTM [144]. One study pointed out that a detailed manikin shape can provide more accurate predictions (4–10%) at some locations, especially those close to the body [141]. Another study revealed that although human body simplification only affected airflow field prediction of the thermal plume regions, it could increase the predictive error of contaminant transport in the whole computational domain [145]. Furthermore, a study proposed a surface smoothing approach that can decrease the number of grids on the CTM surface to reduce the mesh independence effect, improving computational efficiency while obtaining a reliable evaluation. This illustrates the importance of keeping the key body features when simplifying human body [146]. Therefore, it might be possible to utilize more accurate representation of critical body surfaces, such as the head and face, while simplifying the other human body surfaces.

4.2.2. Heat Exchange between Human Body and Microenvironment

The current CFD modeling technique only accounts for heat exchange between the human body and the microenvironment, instead of within the human body, i.e., CTMs do not have the thermoregulation functions that can respond to the microenvironment as human bodies do [147]. The CTMs are mainly used for heat dissipation from the human body, including heat loss by radiation, sensible heat loss by convection, and latent heat loss by convection [62]. For instance, the human body releases its metabolic heat production to the surrounding environment, 49.3 W by convection, 65.1 W by radiation, 41.3 W by evaporation, and 14.8 W by respiration [136]. Moreover, a total of 115 W for a seated person with very light work was generated, including 70 W of sensible heat loss due to convection and radiation, and 45 W of latent heat loss induced by evaporation and respiration, representing the amount of heat released to the surrounding environment to achieve thermal balance [115].
As the PV flow involves mixed convection in most articles, the buoyancy effect should be accounted for. The commonly used method to include thermal buoyancy utilizes the Boussinesq approximation for incompressible flow, and an incompressible ideal gas model. The Boussinesq approximation is applied in indoor airflows as long as the density variations are small or β ( T T O ) 1 [135], which cannot be combined with the species transport equation in the RANS model [148]. The study found that buoyancy tends to suppress turbulence at a stable stratification but promotes it at an unstable stratification [149]. Turbulence production or dissipation by buoyancy can be incorporated when using SKE, RNG, SKE, SST (using a user-defined function), and RSM models. Therefore, the full buoyancy-effects option may be selected to consider the effect of buoyancy on the turbulent dissipation rate (ε), since it is not well understood [148]. The predictive accuracy of the two models above is evaluated in mixed ventilation except in the case of a high-temperature heat source [150].
For radiation models in 20 of 60 studies, the surface-to-surface (S2S) and discrete ordinate radiation (DO) models were mostly applied in the reviewed papers in ANSYS Fluent [148], while S2S and the participating radiation model were selected in the Star CCM+ program [151]. We do not intend to explain the radiation models, whose details can be found in program tutorials. Instead, we provide the following analysis for setting boundary conditions once the radiation is considered. A review evaluated three thermal radiation approaches when simulating normal indoor thermal flows in terms of predictive accuracy and computational time in the ANSYS CFX program [152], demonstrating the importance of including thermal radiation in the CFD model [153].

4.2.3. Boundary Conditions

The thermal boundary condition for CTMs representing the heat source can be modeled in three ways: constant surface temperature (35 of 60 studies), constant surface heat flux (12 of 60 studies), and total heat power output through the volume called a constant volumetric heat generation for the human body (5 of 60 studies). The last is popularly used in the Airpak and PHOENICS programs, and is employed in the ANSYS Fluent program, although the study did not consider the PV case [154].
In general, the CTM surface temperature/heat flux can be selected from the measured data [74,79,92,110] or estimated from the literature and/or experience for simplification [38,80]. For simulations disregarding radiation, the convective heat flux per unit area was determined for the CTM surface boundary condition, e.g., as 39 W/m2 [117,120]. When accounting for radiation effects, some articles still utilized surface heat flux [38,80,111,112]. As the studies simplified the boundary condition setup with a uniform heat flux value, it was necessary to determine the convection-to-radiation ratio (C:R) to properly calculate the convective heat flux. Many values of C:R were found in the reviewed literature, such as 43:57 for a seated person with a surface area of 1.69 m2 [136], 40:60 for a seated person with a surface area of 1.59 m2 and sensible heat loss of 89 W [52], 42:58 for a seated person with 1.45 m2 surface area in a uniform radiation environment [155,156], 40:60 for a standing person with 1.48 m2 surface area and 76 W [56], 30:70 for a standing person with 1.46 m2 surface area [144], 38:62 for a seated person with 1.48 m2 surface area [157], and 60:40 for a seated person using CFD radiation [158]. Note that those variable values emphasize the importance of modeling radiative heat flux between human body and the microenvironment together with the convective heat fluxes.
Note that it is expensive to calculate radiative heat transfer from a human body to room walls by using a surface-to-surface radiation model, because of the large number of segment surfaces accounting for the irregular human body shape [63,159]. Therefore, radiation heat transfer is commonly avoided by assuming that a certain percentage (usually half) of CTMs’ heat loss was caused by convection. The heat flux boundary condition is then used by calculating the ratio of convective heat loss at CTM surface area. If the research is only aimed to globally capture the effect of the thermal plume, the convection-to-radiation ratio might be sufficient [160]. Review found that 11 studies utilized one uniform heat flux value. However, research has shown that it was erroneous to neglect radiation modeling when heat flux boundary conditions were used [157,161]. Using heat flux cannot represent the exact heat-loss portions of convection, radiation, and evaporation, because of their nonlinear relationships, and their local distributions in the segments. Therefore, it is not realistic to only consider convective heat transfer around the human body in an indoor environment, especially when radiation plays a significant role.
A total of 12 studies adopted one uniform temperature value for whole body surfaces because they were focused on PV system properties that were not significantly affected by non-uniform heat transfer around the human body. A study pointed out that when a predefined temperature was used for the body surface, the heat transfer was independent of thermal radiation, which might be disregarded in the simulations [52]. To avoid neglecting about half of the heat loss, using a temperature boundary condition is crucial for proper specification of CTM heat exchange between surfaces and the ambient environment [115]. In addition, even when the PV is used to provide clean air to inhalation, instead of providing local cooling or heating, it is also important to take into consideration the non-uniform temperature distribution over a human body, because of its impact on thermal plume around the human body. Before being inhaled, personalized air needs to penetrate thermal plume first. Therefore, it is relatively correct to define the actual surface temperatures in a realistic environment, and then the effects of radiation can be included without a radiation model [157]. However, obtaining a detailed description of the CTM surface temperature distribution is difficult in an experiment. The commonly used method is to obtain the average temperature for the major segment surface and set it as the constant temperature boundary condition [160]. In addition, coupling CFD with a thermoregulation model was used to define the thermal boundary condition, which will be presented in next section. The representative surfaces, including the indoor floor and wall temperature, are best measured in advance. This can be a drawback because the surface temperature should first be determined. Overall, it is more accurate to set non-uniform thermal boundary conditions when seeking local evaluation indices, e.g., local thermal comfort and specified segment temperatures. There is no standard for the number of segments required for indoor environment studies. The number of segments may be consistent with the manikin number used in the experiments. Studies of segment numbers are listed in Table 1.
Although it was demonstrated that breathing of a sedentary manikin had a very small impact on the occupant thermal plume [162], discrepancies in the healthy manikin’s exposure to exhaled contaminated air still existed between simulations and experiments, which could possibly be attributed to disuse of the transient inhalation and exhalation [91]. Therefore, some studies (10 of 60) considered the breathing process of CTMs, although most only took into account the steady-state inhalation [36,37,69,71] or exhalation process [38,86,87], or both [101], instead of a completely transient breathing process [82,85]. The constant inhaled flow rate was about 0.14 L/s [37,71], and the exhaled flow rate using the nose or mouth ranged from 0.02 L/s [87] to 0.14 L/s [37,91,103], with an exhaled air temperature of 34 °C [38,87,101] or 35 °C [86], depending on the activity level. Note that owing to the research topics that are mainly predicting particle transport and assessing ventilation effectiveness, most studies presented that particle-generation characteristics were considered as horizontal generation steady rate with mass flow rate of 5 × 10−5 kg/s from normal breathing of the infected person [103,104].
For the transient breathing process, a study described the inhalation and exhalation process as a sinusoidal curve function consisting of inhalation (2.5 s), exhalation (2.5 s), and a pause (0.5 s), with an amplitude flow rate of 0.3 L/s [82]. This study demonstrated the usefulness of CFD in investigating the influence of the breathing process on inhaled air quality, and the necessity of proper breathing simulation to understand the complex airflow interactions in the breathing zone. Another study compared steady and transient breathing methods, and recommended the use of steady-state inhalation as the breathing method instead of exhalation, as it does not increase the complexity of the simulation compared to unsteady methods [85]. In addition, it is necessary to model the transient breathing process when studying re-inhaled air and the transport of exhaled contaminants, although exhalation concentration was found to have little effect on the inhaled concentration because of the strong momentum of the PV air supply to carry the pollutant away from the breathing zone. Given the limited number of studies addressing this issue, more studies are suggested to evaluate the effect of a steady or transient breathing process. Note that the detailed respiration cycle and its flow rate were reviewed by Gao and Niu [63].

4.2.4. Coupling of CFD and Thermoregulation Models

CTMs do not include thermoregulation models to define the heat output of the body by default. However, non-uniform microenvironments and transient environments are typically encountered in buildings and in situations where a person’s job involves exposure to high- or low-temperature surfaces or sources [159]. Therefore, knowledge of human thermal response to non-uniform and transient environments is important to accurately provide thermal boundary conditions to a CFD solver using thermoregulatory models. The coupled process includes the data exchange with respect to body segment heat fluxes and temperatures, and detailed local microenvironment parameters [163]. An early PV study coupled CFD and the thermoregulation model to investigate PV in regard to inhaled air quality and local thermal comfort [37]. This model consisted of 16 body segments and considered heat and moisture transfer as well as convection, conduction, and radiation [164,165]. Similarly, another study examined the thermal adaptive effect by changing the human body’s orientation or position relative to spot airflow by means of CFD coupled with Fanger’s neutral model [166]. A CFD study for PV integrated with a transient bio-heat model included the computation of segmental skin temperatures, overall thermal sensation, and overall thermal comfort [113]. In addition to PV studies, coupling CFD and a human-body thermoregulation model has been extensively employed in common studies in indoor environments [147,155,156,161,167,168].
There are already several comprehensive reviews of thermoregulation models [163,169,170,171]. The accuracy of transferred parameters of the environment around the human body in PV studies is vital to the accuracy of CFD simulation results. The quality of coupling is highly dependent upon CTM geometry, heat transfer around CTMs, clothing, and computational resources [161]. In addition, the coupling process might hinder the analyses of PV system performance if convergence is difficult to achieve [172]. Thus, the primary task for the application of coupled CFDs with thermoregulation models is still to ensure the accuracy and efficiency of CTM simulations.

4.3. Modeling Performance

Any CFD study, including studies of PV systems, needs attention to modeling performance focused on computational grids, convergence criteria, and validation.

4.3.1. Computational Grids

The computational grid has attracted more attention with respect to quality control [52], verification [56], evaluation errors [158], and sensitivity [173] because the computational results depend greatly on the grid quality, particularly due to complex CTM geometry [63]. This section outlines recent progress in the use of computational grids for PV studies, and suggests grid-independence analyses using evaluation indices.
Table 2 summarizes the grids in the computational domain and around CTMs for a few studies. Due to the complex human body geometry, hybrid grids—including tetrahedral, hexahedral, and prismatic grids—were commonly adopted in the simulation domain. Figure 3 shows different types of hybrid grids created in a room and around a human body. The prismatic grids were used in the human body boundary layer. In early PV studies, in addition to the boundary layer (prismatic grid), hexahedral combined with tetrahedral grids were typically used for the ambient domain, owing to the utilization of Gambit [174]. Tetrahedral and prismatic grids were mainly used to study complex CTM shapes and indoor configurations [175]. Honeycomb and hexahedral grids were used for the majority of PV studies in Star CCM+ [39,115].
The total grid number reached 5.77 million and 6.5 million, representing 0.27 million grids and 0.20 million grids per cubic meter, respectively, in two recent studies [118,120]. However, two early studies used only 1.81 million and 1.75 million grids, resulting in 0.12 million/m3 and 0.11 million/m3, respectively [36,69]. Future studies will employ larger grid numbers in line with increasing computational capacity. However, CFD studies may become more inefficient than experimental studies, even with the development of computer technology, especially when considering dynamic simulation scenarios. Therefore, a balance between efficiency and accuracy of simulations must be considered.
The range for the initial inflated prismatic height adjacent to the human body varied from 0.4 mm [114] to 2.2 mm [91], and y+ values on the majority of the body’s surface were lower than 4.5 in all cases, which is consistent with non-PV studies [56,158]. Most studies utilized a maximum CTM surface grid size of up to 20 mm. Based on this literature review, a prismatic layer number of at least 3 is suggested, with the expansion ratio under 1.2, which results in sufficiently good grid quality represented by equi-volume skewness within the recommended value of CFD solvers [56,160]. Moreover, y+ varies with the location, and the maximum value occurs mostly in the breathing zone close to the body, so the local grid refinement is required to place the first grid point within the viscous sublayer [113].
Verification of CFD includes the grid convergence index (GCI), which is used to compare discrete solutions at different grid spacing [44,176]. Note that GCI for validation of CFD was suggested for the grid refinement analysis [52]. However, this approach received little attention in PV studies. This is probably due to the difficulty in selecting variables and locations for comparisons, so a new normalized root mean square error index (RMSE*, less than 10%) developed from the GCI index, and a numerical viscosity analysis method may be applied for a grid independence analysis in indoor environments [173]. Grid independence analysis must be sufficiently considered, instead of simply comparing absolute errors or relative errors for limited variables in an ambient environment or on the body surface.

4.3.2. Convergence Criteria

The default convergence criteria in CFD programs are not always convincible when the residual root mean square error value is below the criterion (typically 10−4 or 10−5) [52]. Therefore, besides the condition when the residual values reach an acceptable level for a steady-state condition, two additional assessment indices are required to ensure convergence. One is that the overall imbalance for the specific variables in the domain should be less than 1%, and the other is that the monitored point for the values of interest should reach a steady solution.
In PV studies, the convergence criterion was reached not only when the scaled residuals met the requirement, but also when velocity magnitude at specific points in the thermal plume and breathing zone had stabilized [74,75]. A net heat flux of less than 1% of the total heat gained in the domain was used as the convergence criterion [113,120]. Another study examined the sensitivity of iterative convergence in an indoor environment, and indicated that iterative convergence criteria have a large impact on the results of the age of air [154]. It was suggested that convergence can be reached when the net heat flux imbalance of the domain boundary is less than 1%. In addition, the criteria were employed when the convective and sensible heat losses from the human body were stabilized with the iteration [107]. Another study judged numerical convergence using scaled residuals, net heat flux of total heat gain, and stabilization of CO2 concentration in the breathing zone [93]. Overall, PV studies have clearly taken into account the overall convergence criteria.

4.3.3. Validation

The fundamental strategy of validation is to assess how accurately the computational results compare with the experimental data, with quantified error and uncertainty estimates for both [45]. Unfortunately, no sufficient estimates of uncertainty or error were provided for CTMs until a detailed validation study was conducted in an indoor environment [56,160]. Through the validation, the aspects of grid convergence, turbulence model, and radiation model were identified and qualified, accounting for errors associated with both the computational model and experiment. Practically, it is not possible for every study to pursue a comprehensive validation. Note that the majority of studies selected a computational model based on research experience or references. In 2003, a study established a benchmark test case including a CTM in an MV or DV system for validation of CFD models [177]. A subsequent study provided another benchmark test case involving PV combined with a general ventilation system in a room from the experiment [178,179]. The latter test case has been applied in several publications [39,74,76,77,85]. Therefore, it is a useful database to compare simulated results with a PV benchmark test to validate new programs, turbulence models, grid dependency, and numerical schemes [30].

5. PV Air Supply, Thermal Comfort, and Energy Savings

The most important design considerations for a PV system include air supply characterized by air supply diffusers and air supply parameters. The performance of PV air supply defines how successful such a system is in delivering desired local thermal comfort and inhaled air quality with potential energy savings for the central HVAC unit.

5.1. PV Air Supply Diffusers

An experimental study validated that the air terminal device within three different layouts of blinds with the same diffuser size and effective area had a significant influence on the occupant’s thermal sensations [180]. Therefore, correct modeling of a PV air supply diffuser is an important prerequisite to accurately predict the microenvironment around a human body. The commonly used models involve the simplified boundary condition method, box-method, prescribed velocity method, and momentum method [30]. Note that the air supply diffusers for the background ventilation system are not discussed in detail, as extensive studies have investigated diffuser boundary conditions without PV diffusers [181,182,183]. This literature review reveals that the majority of PV studies used the momentum method, which is specified as a momentum source term in the conservation equation and is calculated from the mass flow rate and effective diffuser area. This is probably because the relatively small air supply opening area is incapable of modeling the detailed diffuser geometry because it requires a large number of grids. Another explanation is that the momentum method is more effective and saves time in setting diffuser boundary conditions. It is noticeable that few studies have paid specific attention to the effect of real diffuser geometry, which has mostly occurred in general CFD cases of indoor microenvironments in recent years, e.g., inserted lobes [184], blade of louver [185], and swirl blades [186].
An exception for PV studies utilized the measured data in a PV and UFAD air supply diffuser as a boundary condition, as shown in Figure 4, and reported that the inlet velocities have non-uniform distributions [39]. For instance, the largest velocity occurred in the central two slots of the PV diffuser, especially for large flow rates. Thus, given that PV flow directly affects the occupant’s body within a short distance, the shape, configuration, or inner components in air supply diffusers—e.g., guide vanes, perforated plates, sound dampers, and curved surfaces—can have a positive influence on thermal comfort [180,187] and ventilation effectiveness [188]. Future studies using PV air supply diffusers with their actual geometries are necessary to evaluate their effect on the microenvironment around the human body [30].

5.2. PV Air Supply Parameters

To evaluate the performance of PV systems, CFD studies require PV air supply parameters, such as air velocity or flow rate, temperature, air humidity, and trace gas concentration. This review aims to provide an overview of CFD applications in PV studies, suggestions for CFD simulation setups, and further research interests. Thus, air supply parameters for typical PV with a background ventilation system are not reviewed in detail. Instead, we compare commonly used parameters, such as flow rate, air temperature, and turbulence intensity.
Table 3 summarizes air supply parameters in a typical PV and background ventilation system. It shows that the majority of studies selected an airflow rate from 2.4 to 20 L/s, and an air temperature from 16 to 24 °C to meet occupant needs. For desk fans or chair fans, the supply boundary conditions can be set as fan boundary conditions with predefined pressure changes according to the flow rates [99,100], or inlet boundary conditions with an assumed temperature of ambient air [79]. Background ventilation systems typically have flow rates of 14 to 60 L/s per person and air temperatures from 16 to 24 °C. The air turbulence intensity mostly varied from 2.5% to 10% for PV and 2.5% to 15% for background ventilation. The distances between PV supply openings to an occupant’s face/head ranged from 0.2 to 0.6 m for an air terminal device around a desk, or about 1.4 m for PV mounted in the ceiling.
Figure 5 shows two special PV nozzles, including a reduced-mixing personal ventilation jet [74] and low-mixing coaxial nozzle mounted in the ceiling (CMPV) [89,93]. For example, CMPV has two nozzles, the primary nozzle delivering outside air at a flow rate ranging from 5 to 11 L/s, and the secondary nozzle supplying recirculated air at a flow rate ranging from 10 to 20 L/s. The review for PV and background ventilation air supply parameters are given as further references. Not all of the studies provide the minimum information that should be included in CFD-related papers, as required [43].
As a powerful tool, CFD has a significant advantage in capturing flow characteristics and particle trajectories. However, the majority of studies, especially in the early stage of development, only accounted for the steady-state condition. For the indoor environment, a study pointed out that airflow cannot be expected to be in steady-state condition, and this can only be predicted for a given flow with time-dependent equations [30]. Another study found that air-movement parameters, such as turbulence intensity and fluctuation frequency, can significantly affect human thermal perception, and suggested that dynamic airflows can achieve stronger cooling effects and better thermal comfort evaluation without negatively impacting occupants’ work performance [189]. A few well-known unsteady studies utilized CFD to assess the performance of PV under the influence of a moving person [73,94], the breathing process [82,85], and sinusoidal airflow with different fluctuation frequencies [113,119,120]. One aspect of dynamic airflow is characterized by perceived air movement with a range of fluctuation frequency, e.g., 0.5 and 1.0 Hz [190], where the range of associated amplitudes and mean velocities is determined by minimum and maximum jet flows [113]. Suitable air velocities enable the penetration of the thermal plume around the human body to reach the breathing zone. Furthermore, an acceptable temperature range to deliver cold/warm air to the human body with ambient air is guaranteed to maintain better thermal comfort.
The fluctuation frequency in PV studies using CFD has only gained attention in the last two years, while its effect on human thermal response was extensively examined in experiments starting in the early 21st century [191,192,193]. Compared to constant flow, a sinusoidal airflow with the right frequency can have a greater cooling effect on the same average velocity [194]. An average flow rate combined with a specific frequency can guarantee good breathable air quality and comfort, providing a good compromise between thermal comfort and air quality, and better protection from direct contamination [120]. Therefore, to some extent, to use average flow rate solely in steady state is not an adequate simulation method. Moreover, other aspects of dynamic airflow involve the frequent change of an occupant’s posture and position, transport of simulated natural wind, and coughing and sneezing behavior. All of these aspects have an important impact on PV performance, which can be fulfilled by CFD technique. Therefore, further CFD work is required to account for the use of PV in more realistic environments.

5.3. Performance Evaluation Indices

Several evaluation criteria emerged to assess performance of PV systems, including inhaled air quality, thermal comfort, and energy savings.

5.3.1. Inhaled Air Quality

PV considers the occupants’ health, comfort, and productivity around a human microenvironment. Its performance can be assessed based on occupant-related performance criteria, e.g., improvement of inhaled air quality, improvement of thermal comfort, and protection from and minimization of contaminant transmission [3]. Evaluation indices for ventilation effectiveness and inhaled air quality have been introduced in the literature [6]. However, these criteria were developed for experimental studies, and recent studies employ new indices specifically developed for CFD studies [88,103,113,119]. Table 4 summarizes commonly used evaluation indices for ventilation effectiveness and inhaled air quality. It is found that, due to the powerful post-processing functions in CFD programs, relevant parameters in specific locations/regions can be monitored. For instance, a study calculated the ventilation effectiveness based on tracer gas concentrations at the breathing zone that was defined as a small sphere of 1 cm radius at 2.5 cm from the occupant’s nose [88,103]. Moreover, the cross-infection level can be assessed by monitoring the number of particles deposited in the vicinity of an exposed person and the number of particles generated by an infected person using deposited fraction (DFr) [103,120]. A few other evaluation indices have been applied in PV studies using CFD, e.g., new scale for ventilation efficiency 3 (SVE3*) and scale for ventilation efficiency 4 (SVE4) [101], age of air (AGE) [81], and concentration asymmetry (ΔC) [117]. Note that SVE3 and SVE4 respectively represent the mean age and residual life of air, without considering where air is from and where it is exhausted, while SVE3* and SVE4* respectively represent the mean age and residual life of air from a specific inlet or outlet [195]. They have been used in few PV studies; details can be found in Table 1.

5.3.2. Thermal Comfort

Regarding thermal comfort involving air temperature, radiant temperature, air velocity, humidity, clothes insulation, and activity, it is recommended to consider all six factors for human response, e.g., by using an integrated thermal comfort model. A recent advanced thermal comfort (ATC) model [196,197,198] is widely accepted for assessing segmental/overall body thermal comfort and sensation using CFD simulation for PV studies [98,118]. A thermal comfort scale from −4 (very uncomfortable) to +4 (very comfortable), and thermal sensation scale from −4 (very cold) to +4 (very hot) are utilized. This comfort model is more comprehensive and more robust than the previous model [199] for PV studies [37], as it accounts for transient variation of the skin and core temperatures. CFD is more likely coupled with a physiological regulation model through exchanging ambient environment data and manikin segmental skin temperature and heat flux in order to obtain input data to predict thermal sensation and comfort with an ATC model [89,93,99,100,113,119].
The CTM-based equivalent temperature (teq) [200,201] has been widely used to assess the thermal comfort of individual segments for PV studies [39,79,82], and is included in ISO 14505 [202]. It is defined by thermal resistance and sensible heat loss of body segments located in a thermal equilibrium environment, accounting for radiative heat release in some cases [39]. Note that this method does not include the thermoregulation model, and its reliability should be questioned when evaporation from skin is involved, but it has been shown to be in agreement with experimental studies and the ATC model [172]. In addition, the draft risk indicates the local draft discomfort at face exposure to the airflow from a PV nozzle [102,107,108], which has already been adopted in the ASHRAE and ISO thermal comfort standards [203,204]. The airflow interaction in the vicinity of the human body is complex and important, especially when PV nozzles supply fresh air to the breathing zone and penetrate the convective boundary layer around the body, so as to provide more acceptable inhaled air quality without draft complaints [205,206].
To assess the influence of a PV system on thermal comfort, mean skin surface temperature, as an important factor influencing a human’s thermal sensation, was employed by calculating Fanger’s thermal comfort model parameters and by considering the heat generation and balance for a sleeping person [109,110]. Finally, the commonly used comfort indices—e.g., predicted mean vote (PMV) and predicated percentage of dissatisfaction (PPD), as well as the improvement of overall thermal sensation (DPMV) [107]—were used to examine an occupant’s thermal comfort [123]. The impact of PV on CTM thermal comfort is mostly evaluated using the ATC model.

5.3.3. Energy Savings

Significant energy saving can result from maintaining comfort in a wider range by increasing the cooling temperature setpoint or decreasing the heating temperature setpoint [11,48]. CFD can generally not be directly applied to evaluate energy saving, but can provide sufficient information to calculate energy-saving potential compared to traditional mechanical ventilation. For instance, the evaluation criteria are established in terms of the improved temperature setpoint and reduced outdoor airflow rate. The energy-saving potential of PV has been widely investigated by means of energy simulation programs [18,19,50], which are not within the scope of this study. Therefore, only CFD-based valuation indices for energy saving from PV are reviewed in this study.
The direct energy-saving evaluation index mainly accounts for the reduction of cooling loads in the occupied zone using PV devices [96,109,207], or a difference in power consumption of a PV or MV air supply fan due to different flow rates [113]. A cooling load calculation method originally designed for stratified air distribution systems was used [208] by ventilating only the occupied zone and then calculating the cooling energy saving in the coil capacity of DV compared to MV [123]. The coil capacity can be assessed by multiplying the exhaust mass flow rate and the difference between the exhaust temperature and indoor setpoint temperature [111,112]. Similarly, the cooling capacity can be calculated by taking into account the supply flow rate and temperature difference between outside air and supplied air [96,207]. As for PV assisting with a chilled ceiling, the radiant heat load removal and convective heat load removal should both be considered [119]. The energy consumption for dehumidifying thermal manikins should be considered in addition to the convective cooling load [110].
The reduction in the power consumption of an integrated air supply fan has been assessed in several studies [88,98,117]. Review found that CMPV performs better than conventional MV at maintaining occupants’ thermal sensation and inhaled air quality [88,89,90]. Owing to the advanced design of coaxial nozzles in CMPV, a fan with about 15 W of nominal power supplied fresh air through the ceiling-mounted primary nozzle, and an additional fan with the same power supplied recirculated air through a secondary nozzle, with a varied fresh airflow rate and total supply flow rate. The resulting reduction in power consumption of the supply fan demonstrated the potential energy saving. An evaluation of energy savings with PV used an indirect valuation index called the cooling efficiency (CE) [107]. It is defined as the ratio of additional sensible heat loss removed from the human body to the PRAC’s cooling capacity. Note that a portion of a device’s cooling capacity is fulfilled by dissipating heat with convective airflow around the human body, resulting in a compromise to achieve high cooling performance and acceptable thermal comfort. In addition, the energy utilization coefficient (EUC), reflecting the difference in average operative temperatures between the occupied and unoccupied zone was used in BTAC with a radiant panel [114]. When EUC is greater than one, it can save energy because less cooling energy is used to remove the heat in the unoccupied zone.
A major finding from reviewing energy saving evaluations for PV systems is that most evaluation indices focus on the terminal device-based energy consumption without considering the effect of the building envelope. Although two articles studied the energy use due to the variation of the envelope heat gain in a bedroom at night [109,110], variations of the envelope heat gain were specified in advance instead of being utilized as an evaluation index in a wider analysis framework. Moreover, the weather data file, equipment schedules, and more detailed control strategies were not considered, resulting in rough evaluations of potential energy savings. Therefore, regardless of whether the direct or indirect method was used, the indices provided a preliminary evaluation of the energy-saving potential of PV systems, and future studies should account for more comprehensive system changes due to the PV presence.

6. Conclusions

This paper has reviewed studies of CFD analysis methods for PV systems in indoor built environments. First, this study found that a significant increase in the number of PV system studies was presented using CFD with diverse PV system types associated with different background HVAC systems. Then, besides the basic elements of CFD simulation for PV study, such as turbulence models, boundary condition, computational grids, convergence criteria, and validation of CFD model, this study specifically investigated the modeling of human body and personalized air supply. For human body, this study focused on CTM, deployed in CFD studies on PV to account for its impacts on local thermal state. This study reviewed the influences of CTM body shape, heat exchange over skin surface, and inclusion/exclusion of heat radiation. Overall, the inclusion of radiative heat transfer is fairly important, especially in spaces with radiant floor/panel or external building walls and windows. For personalized air supply, this study reviewed the set up for PV system parameters and evaluation indices for PV performance on inhaled air quality, thermal comfort, and energy-saving. Overall, reviewed studies have demonstrated that CFD is a powerful tool for predicting inhaled air quality, occupants’ thermal comfort, and energy-saving potential in buildings. Nevertheless, current studies do not examine interactions of an integrated system that would include thermal comfort of individual occupants, distributed PVs, and other building systems. This integrated approach could uncover major PV system advantages and disadvantages prior to costly implementations in actual buildings. Furthermore, accurate CFD results require careful consideration of the numerical setup, and understanding this is a prerequisite for optimizing PV performance, integration with other building systems, and promoting its application. Through state-of-the-art analyses of current CFD applications in studying PV systems, future CFD procedures and research interests are suggested to provide a preliminary guideline of PV analyses for building system engineers and researchers.

Author Contributions

Conceptualization, J.L. and J.S.; methodology, J.L., S.Z. and J.S.; investigation, J.L. and M.K.K.; supervision, J.S.; writing – original draft, J.L.; writing – review & editing, S.Z., M.K.K. and J.S.; funding acquisition, J.L. and M.K.K.

Funding

This research was funded by National Natural Science Foundation of China (51608310) and Science and Technology Plan Project of University in Shandong Province (J16LG07).

Acknowledgments

This study is supported by Cluster of Sustainability at the University of Maryland (CITY@UMD). Moreover, Shandong Province Green Building Collaborative Innovation Center Innovation Team Support Program is specially acknowledged.

Conflicts of Interest

The authors declared that they have no conflicts of interest to this work.

References

  1. Fanger, P.O. Indoor Air Quality in the 21st Century: Search for Excellence. Indoor Air 2000, 10, 68–73. [Google Scholar] [CrossRef] [PubMed]
  2. Fanger, P.O. Human requirements in future air-conditioned environments. Int. J. Refrig. 2001, 24, 148–153. [Google Scholar] [CrossRef]
  3. Melikov, A.K. Personalized ventilation. Indoor Air 2004, 14, 157–167. [Google Scholar] [CrossRef] [PubMed]
  4. Wargocki, P.; Wyon, D.P.; Baik, Y.K.; Clausen, G.; Fanger, P.O. Perceived Air Quality, Sick Building Syndrome (SBS) Symptoms and Productivity in an Office with Two Different Pollution Loads. Indoor Air 1999, 9, 165–179. [Google Scholar] [CrossRef] [PubMed]
  5. Melikov, A.K.; Cermak, R.; Majer, M. Personalized ventilation: Performance of different air terminal devices. In Proceedings of the CLIMA 2000, Napoli, Italy, 15–18 September 2001. [Google Scholar]
  6. Melikov, A.K.; Cermak, R.; Majer, M. Personalized ventilation: Evaluation of different air terminal devices. Energy Build. 2002, 34, 829–836. [Google Scholar] [CrossRef]
  7. Bolashikov, Z.D.; Nikolaev, L.; Melikov, A.K.; Kaczmarczyk, J.; Fanger, P.O. Personalized ventilation: Air terminal devices with high efficiency. In Proceedings of the 7th Healthy Buildings 2003, Singapore, 7–11 September 2003. [Google Scholar]
  8. Hoyt, T.; Arens, E.; Zhang, H. Extending air temperature setpoints: Simulated energy savings and design considerations for new and retrofit buildings. Build. Environ. 2015, 88, 89–96. [Google Scholar] [CrossRef] [Green Version]
  9. Hoyt, T.; Lee, K.H.; Zhang, H.; Arens, E.; Webster, T. In Energy savings from extended air temperature setpoints and reductions in room air mixing. In Proceedings of the International Conference on Environmental Ergonomics, Boston, MA, USA, 2–7 August 2009. [Google Scholar]
  10. Melikov, A.K. Advanced air distribution: Improving health and comfort while reducing energy use. Indoor Air 2016, 26, 112–124. [Google Scholar] [CrossRef]
  11. Heidarinejad, M.; Dalgo, D.A.; Mattise, N.W.; Srebric, J. Personalized cooling as an energy efficiency technology for city energy footprint reduction. J. Clean Prod. 2018, 171, 491–505. [Google Scholar] [CrossRef]
  12. Melikov, A.K.; Skwarczynski, M.A.; Kaczmarczyk, J.; Zabecky, J. Use of personalized ventilation for improving health, comfort, and performance at high room temperature and humidity. Indoor Air 2012, 23, 250–263. [Google Scholar] [CrossRef]
  13. Amai, H.; Tanabe, S.-I.; Akimoto, T.; Genma, T. Thermal sensation and comfort with different task conditioning systems. Build. Environ. 2007, 42, 3955–3964. [Google Scholar] [CrossRef]
  14. Bogdan, A.; Chludzinska, M. Assessment of Thermal Comfort Using Personalized Ventilation. HVAC&R Res. 2010, 16, 529–542. [Google Scholar]
  15. Chen, Y.; Raphael, B.; Sekhar, S.C. Experimental and simulated energy performance of a personalized ventilation system with individual airflow control in a hot and humid climate. Build. Environ. 2016, 96, 283–292. [Google Scholar] [CrossRef]
  16. Makhoul, A.; Ghali, K.; Ghaddar, N. The Energy Saving Potential and the Associated Thermal Comfort of Displacement Ventilation Systems Assisted by Personalised Ventilation. Indoor Built Environ. 2012, 22, 508–519. [Google Scholar] [CrossRef]
  17. Sekhar, S.C.; Gong, N.; Tham, K.W.; Cheong, K.W.; Melikov, A.K.; Wyon, D.P.; Fanger, P.O. Findings of Personalized Ventilation Studies in a Hot and Humid Climate. HVAC&R Res. 2005, 11, 603–620. [Google Scholar]
  18. Schiavon, S.; Melikov, A.K.; Sekhar, C. Energy analysis of the personalized ventilation system in hot and humid climates. Energy Build. 2010, 42, 699–707. [Google Scholar] [CrossRef]
  19. Schiavon, S.; Melikov, A.K. Energy-saving strategies with personalized ventilation in cold climates. Energy Build. 2009, 41, 543–550. [Google Scholar] [CrossRef] [Green Version]
  20. Melikov, A.K. Advanced air distribution. ASHRAE J. 2011, 53, 73–77. [Google Scholar]
  21. Chen, Y.; Raphael, B.; Sekhar, S.C. Individual control of a personalized ventilation system integrated with an ambient mixing ventilation system. HVAC&R Res. 2012, 18, 1136–1152. [Google Scholar]
  22. Cermak, R.; Melikov, A.K.; Forejt, L.; Kovar, O. Performance of Personalized Ventilation in Conjunction with Mixing and Displacement Ventilation. HVAC&R Res. 2006, 12, 295–311. [Google Scholar]
  23. Alain, M.; Kamel, G.; Nesreen, G. A simplified combined displacement and personalized ventilation model. HVAC&R Res. 2012, 18, 737–749. [Google Scholar]
  24. Cermak, R.; Melikov, A.K. Air Quality and Thermal Comfort in an Office with Underfloor, Mixing and Displacement Ventilation. Int. J. Vent. 2006, 5, 323–352. [Google Scholar] [CrossRef]
  25. Cermak, R.; Melikov, A.K. Protection of Occupants from Exhaled Infectious Agents and Floor Material Emissions in Rooms with Personalized and Underfloor Ventilation. HVAC&R Res. 2007, 13, 23–38. [Google Scholar]
  26. Li, R.; Sekhar, S.C.; Melikov, A.K. Thermal comfort and indoor air quality in rooms with integrated personalized ventilation and under-floor air distribution systems. HVAC&R Res. 2011, 17, 829–846. [Google Scholar]
  27. Chakroun, W.; Ghaddar, N.; Ghali, K. Chilled ceiling and displacement ventilation aided with personalized evaporative cooler. Energy Build. 2011, 43, 3250–3257. [Google Scholar] [CrossRef]
  28. Lipczynska, A.; Kaczmarczyk, J.; Melikov, A.K. Thermal environment and air quality in office with personalized ventilation combined with chilled ceiling. Build. Environ. 2015, 92, 603–614. [Google Scholar] [CrossRef]
  29. Niu, J.; Gao, N.; Phoebe, M.; Huigang, Z. Experimental study on a chair-based personalized ventilation system. Build. Environ. 2007, 42, 913–925. [Google Scholar] [CrossRef]
  30. Nielsen, P.V. Fifty years of CFD for room air distribution. Build. Environ. 2015, 91, 78–90. [Google Scholar] [CrossRef]
  31. Liu, J.; Heidarinejad, M.; Gracik, S.; Srebric, J. The impact of exterior surface convective heat transfer coefficients on the building energy consumption in urban neighborhoods with different plan area densities. Energy Build. 2015, 86, 449–463. [Google Scholar] [CrossRef]
  32. Liu, J.; Heidarinejad, M.; Pitchurov, G.; Zhang, L.; Srebric, J. An extensive comparison of modified zero-equation, standard k-ε, and LES models in predicting urban airflow. Sust. Cities Soc. 2018, 40, 28–43. [Google Scholar] [CrossRef]
  33. Li, Y.; Nielsen, P.V. CFD and ventilation research. Indoor Air 2011, 21, 442–453. [Google Scholar] [CrossRef]
  34. Liu, J.; Heidarinejad, M.; Gracik, S.; Srebric, J.; Yu, N. An indirect validation of convective heat transfer coefficients (CHTCs) for external building surfaces in an actual urban environment. Build. Simul. 2015, 8, 337–352. [Google Scholar] [CrossRef]
  35. Srebric, J.; Heidarinejad, M.; Liu, J. Building neighborhood emerging properties and their impacts on multi-scale modeling of building energy and airflows. Build. Environ. 2015, 91, 246–262. [Google Scholar] [CrossRef]
  36. Gao, N.; Niu, J. CFD study on micro-environment around human body and personalized ventilation. Build. Environ. 2004, 39, 795–805. [Google Scholar] [CrossRef]
  37. Gao, N.; Niu, J.; Zhang, H. Coupling CFD and Human Body Thermoregulation Model for the Assessment of Personalized Ventilation. HVAC&R Res. 2006, 12, 497–518. [Google Scholar]
  38. Shen, C.; Gao, N.; Wang, T. CFD study on the transmission of indoor pollutants under personalized ventilation. Build. Environ. 2013, 63, 69–78. [Google Scholar] [CrossRef]
  39. Kong, M.; Zhang, J.; Wang, J. Air and air contaminant flows in office cubicles with and without personal ventilation: A CFD modeling and simulation study. Build. Simul. 2015, 8, 381–392. [Google Scholar] [CrossRef]
  40. Chen, Q.; Lee, K.; Mazumdar, S.; Poussou, S.; Wang, L.; Wang, M.; Zhang, Z. Ventilation performance prediction for buildings: Model assessment. Build. Environ. 2010, 45, 295–303. [Google Scholar] [CrossRef] [Green Version]
  41. Hayashi, T.; Ishizu, Y.; Kato, S.; Murakami, S. CFD analysis on characteristics of contaminated indoor air ventilation and its application in the evaluation of the effects of contaminant inhalation by a human occupant. Build. Environ. 2002, 37, 219–230. [Google Scholar] [CrossRef]
  42. Liu, J.; Srebric, J.; Yu, N. Numerical simulation of convective heat transfer coefficients at the external surfaces of building arrays immersed in a turbulent boundary layer. Int. J. Heat Mass Transf. 2013, 61, 209–225. [Google Scholar] [CrossRef]
  43. Sørensen, D.N.; Nielsen, P.V. Guest editorial: CFD In Indoor Air. Indoor Air 2004, 13, 1. [Google Scholar] [CrossRef]
  44. Roache, P.J. Quantification of uncertainty in computational fluid dynamics. Annual Review of Fluid Mechanics 1997, 29, 123–160. [Google Scholar] [CrossRef]
  45. Oberkampf, W.L.; Trucano, T.G. Verification and validation in computational fluid dynamics. Progress Aerosp. Sci. 2002, 38, 209–272. [Google Scholar] [CrossRef] [Green Version]
  46. Stamou, A.; Katsiris, I. Verification of a CFD model for indoor airflow and heat transfer. Build. Environ. 2006, 41, 1171–1181. [Google Scholar] [CrossRef]
  47. Cao, G.; Awbi, H.; Yao, R.; Fan, Y.; Sirén, K.; Kosonen, R.; Zhang, J. A review of the performance of different ventilation and airflow distribution systems in buildings. Build. Environ. 2014, 73, 171–186. [Google Scholar] [CrossRef]
  48. Veselý, M.; Zeiler, W. Personalized conditioning and its impact on thermal comfort and energy performance—A review. Renew. Sust. Energ. Rev. 2014, 34, 401–408. [Google Scholar] [CrossRef]
  49. Zhang, H.; Arens, E.; Zhai, Y. A review of the corrective power of personal comfort systems in non-neutral ambient environments. Build. Environ. 2015, 91, 15–41. [Google Scholar] [CrossRef] [Green Version]
  50. Shahzad, S.; Calautit, J.K.; Calautit, K.; Hughes, B.; Aquino, A.I. Advanced Personal Comfort System (APCS) for the workplace: A review and case study. Energy Build. 2018, 173, 689–709. [Google Scholar] [CrossRef]
  51. Chen, Q.; Srebric, J. A Procedure for Verification, Validation, and Reporting of Indoor Environment CFD Analyses. HVAC&R Res. 2002, 8, 201–216. [Google Scholar]
  52. Sørensen, D.N.; Nielsen, P.V. Quality control of computational fluid dynamics in indoor environments. Indoor Air 2003, 13, 2–17. [Google Scholar] [CrossRef]
  53. Zhai, Z. Application of Computational Fluid Dynamics in Building Design: Aspects and Trends. Indoor Built Environ. 2006, 15, 305–313. [Google Scholar] [CrossRef]
  54. Zhai, Z.J.; Zhang, Z.; Zhang, W.; Chen, Q.Y. Evaluation of Various Turbulence Models in Predicting Airflow and Turbulence in Enclosed Environments by CFD: Part 1—Summary of Prevalent Turbulence Models. HVAC&R Res. 2007, 13, 853–870. [Google Scholar]
  55. Zhang, Z.; Zhang, W.; Zhai, Z.J.; Chen, Q.Y. Evaluation of Various Turbulence Models in Predicting Airflow and Turbulence in Enclosed Environments by CFD: Part 2—Comparison with Experimental Data from Literature. HVAC&R Res. 2007, 13, 871–886. [Google Scholar]
  56. Sideroff, C.N.; Dang, T.Q. Verification and Validation of CFD for the Personal Micro-Environment. ASHRAE Trans. 2008, 114, 45–56. [Google Scholar]
  57. Wang, H.; Wang, H.; Gao, F.; Zhou, P.; Zhai, Z. Literature review on pressure–velocity decoupling algorithms applied to built-environment CFD simulation. Build. Environ. 2018, 143, 671–678. [Google Scholar] [CrossRef]
  58. Kato, S. Review of airflow and transport analysis in building using CFD and network model. Japan Archit. Rev. 2018, 1, 299–309. [Google Scholar] [CrossRef]
  59. Blocken, B. LES over RANS in building simulation for outdoor and indoor applications: A foregone conclusion? Build. Simul. 2018, 11, 821–870. [Google Scholar] [CrossRef] [Green Version]
  60. Nielsen, P.V. Computational fluid dynamics and room air movement. Indoor Air 2004, 14, 134–143. [Google Scholar] [CrossRef] [Green Version]
  61. Chen, Q. Ventilation performance prediction for buildings: A method overview and recent applications. Build. Environ. 2009, 44, 848–858. [Google Scholar] [CrossRef]
  62. Murakami, S. Analysis and design of micro-climate around the human body with respiration by CFD. Indoor Air 2004, 14, 144–156. [Google Scholar] [CrossRef]
  63. Gao, N.; Niu, J. CFD Study of the Thermal Environment around a Human Body: A Review. Indoor Built Environ. 2005, 14, 5–16. [Google Scholar] [CrossRef]
  64. You, R.; Zhang, Y.; Zhao, X.; Lin, C.-H.; Wei, D.; Liu, J.; Chen, Q. An innovative personalized displacement ventilation system for airliner cabins. Build. Environ. 2018, 137, 41–50. [Google Scholar] [CrossRef]
  65. Zítek, P.; Vyhlídal, T.; Simeunović, G.; Nováková, L.; Čížek, J. Novel personalized and humidified air supply for airliner passengers. Build. Environ. 2010, 45, 2345–2353. [Google Scholar] [CrossRef]
  66. Oh, M.S.; Ahn, J.H.; Kim, D.W.; Jang, D.S.; Kim, Y. Thermal comfort and energy saving in a vehicle compartment using a localized air-conditioning system. Appl. Energy 2014, 133, 14–21. [Google Scholar] [CrossRef]
  67. Zhang, H.; Dai, L.; Xu, G.; Li, Y.; Chen, W.; Tao, W.-Q. Studies of air-flow and temperature fields inside a passenger compartment for improving thermal comfort and saving energy. Part I: Test/numerical model and validation. Appl. Therm. Eng. 2009, 29, 2022–2027. [Google Scholar] [CrossRef]
  68. Croitoru, C.; Nastase, I.; Bode, F.; Meslem, A.; Dogeanu, A. Thermal comfort models for indoor spaces and vehicles—Current capabilities and future perspectives. Renew. Sust. Energ. Rev. 2015, 44, 304–318. [Google Scholar] [CrossRef]
  69. Gao, N.; Niu, J. Modeling the Performance of Personalized Ventilation under Different Conditions of Room Air and Personalized Air. HVAC&R Res. 2005, 11, 587–602. [Google Scholar]
  70. Nielsen, P.V.; Hyldgaard, C.E.; Melikov, A.; Andersen, H.; Soennichsen, M. Personal Exposure Between People in a Room Ventilated by Textile Terminals—with and without Personalized Ventilation. HVAC&R Res. 2007, 13, 635–643. [Google Scholar]
  71. Gao, N.; Zhang, H.; Niu, J. Investigating Indoor Air Quality and Thermal Comfort Using a Numerical Thermal Manikin. Indoor Built Environ. 2007, 16, 7–17. [Google Scholar] [CrossRef]
  72. Zhao, B.; Guan, P. Modeling particle dispersion in personalized ventilated room. Build. Environ. 2007, 42, 1099–1109. [Google Scholar] [CrossRef]
  73. Yang, B.; Sekhar, S.C. The Influence of Evenly Distributed Ceiling Mounted Personalized Ventilation Devices on the Indoor Environment. Int. J. Vent. 2008, 7, 99–112. [Google Scholar]
  74. Russo, J.S.; Dang, T.Q.; Khalifa, H.E. Computational analysis of reduced-mixing personal ventilation jets. Build. Environ. 2009, 44, 1559–1567. [Google Scholar] [CrossRef]
  75. Dygert, R.K.; Russo, J.S.; Dang, T.Q.; Khalifa, H.E. Modeling of the Human Body to Study the Personal Micro Environment. ASHRAE Trans. 2009, 115, 407–420. [Google Scholar]
  76. Russo, J.S.; Khalifa, H.E. CFD assessment of intake fraction in the indoor environment. Build. Environ. 2010, 45, 1968–1975. [Google Scholar] [CrossRef]
  77. Russo, J.; Khalifa, H.E. CFD Analysis of Personal Ventilation with Volumetric Chemical Reactions. HVAC&R Res. 2010, 16, 799–812. [Google Scholar]
  78. Conceição, E.Z.E.; Lúcio, M.M.J.R.; Rosa, S.P.; Custódio, A.L.V.; Andrade, R.L.; Meira, M.J.P.A. Evaluation of comfort level in desks equipped with two personalized ventilation systems in slightly warm environments. Build. Environ. 2010, 45, 601–609. [Google Scholar] [CrossRef]
  79. Tham, K.W.; Pantelic, J. Performance evaluation of the coupling of a desktop personalized ventilation air terminal device and desk mounted fans. Build. Environ. 2010, 45, 1941–1950. [Google Scholar] [CrossRef]
  80. He, Q.; Niu, J.; Gao, N.; Zhu, T.; Wu, J. CFD study of exhaled droplet transmission between occupants under different ventilation strategies in a typical office room. Build. Environ. 2011, 46, 397–408. [Google Scholar] [CrossRef]
  81. Zhai, Z.; Metzger, I.D. Taguchi-Method-Based CFD Study and Optimisation of Personalised Ventilation Systems. Indoor Built Environ. 2011, 21, 690–702. [Google Scholar]
  82. Mazej, M.; Butala, V. Investigation in the Characteristics of the Personal Ventilation Using Computational Fluid Dynamics. Indoor Built Environ. 2011, 21, 749–771. [Google Scholar] [CrossRef]
  83. Adamu, Z.A.; Cook, M.J.; Price, A.D.F. Natural Personalised Ventilation—A Novel Approach. Int. J. Vent. 2011, 10, 263–275. [Google Scholar] [CrossRef]
  84. Ishiguro, R.; Chikamoto, T.; Hashimoto, S.; Inada, R.; Nishino, A.; Akimoto, T. Airflow Control for Personal Air-Conditioning in aPartly Unoccupied Zone using a Multi-Flow Ceiling Cassette Type Packaged Air-Conditioner. Int. J. Vent. 2011, 9, 393–403. [Google Scholar] [CrossRef]
  85. Russo, J.S.; Khalifa, E. Computational study of breathing methods for inhalation exposure. HVAC&R Res. 2011, 17, 419–431. [Google Scholar]
  86. Li, X.; Niu, J.; Gao, N. Co-occupant’s exposure to exhaled pollutants with two types of personalized ventilation strategies under mixing and displacement ventilation systems. Indoor Air 2012, 23, 162–171. [Google Scholar] [CrossRef]
  87. Kanaan, M.; Ghaddar, N.; Ghali, K. Quality of inhaled air in displacement ventilation systems assisted by personalized ventilation. HVAC&R Res. 2012, 18, 500–514. [Google Scholar]
  88. Makhoul, A.; Ghali, K.; Ghaddar, N. Desk fans for the control of the convection flow around occupants using ceiling mounted personalized ventilation. Build. Environ. 2013, 59, 336–348. [Google Scholar] [CrossRef]
  89. Makhoul, A.; Ghali, K.; Ghaddar, N. Thermal comfort and energy performance of a low-mixing ceiling-mounted personalized ventilator system. Build. Environ. 2013, 60, 126–136. [Google Scholar] [CrossRef]
  90. Makhoul, A.; Ghali, K.; Ghaddar, N.; Chakroun, W. Investigation of particle transport in offices equipped with ceiling-mounted personalized ventilators. Build. Environ. 2013, 63, 97–107. [Google Scholar] [CrossRef]
  91. Yang, J.; Sekhar, C.; Wai, D.C.K.; Raphael, B. Computational fluid dynamics study and evaluation of different personalized exhaust devices. HVAC&R Res. 2013, 19, 934–946. [Google Scholar]
  92. Cheong, D.K.W.; Huang, S. Performance evaluation of personalized ventilation system with two types of air terminal devices coupled with displacement ventilation in a mock-up office. HVAC&R Res. 2013, 19, 974–985. [Google Scholar]
  93. Makhoul, A.; Ghali, K.; Ghaddar, N. Low-mixing coaxial nozzle for effective personalized ventilation. Indoor Built Environ. 2013, 24, 225–243. [Google Scholar] [CrossRef]
  94. Yang, B.; Sekhar, C. Interaction of dynamic indoor environment with moving person and performance of ceiling mounted personalized ventilation system. Indoor Built Environ. 2014, 23, 920–932. [Google Scholar] [CrossRef]
  95. Russo, J.; Ezzat Khalifa, H. Effect of Chemical Reactions in the Personal Micro-environment on Inhaled Air Quality. ASHRAE Trans. 2014, 120, 1–8. [Google Scholar]
  96. Naumov, A.L.; Tabunshchikov, I.A.; Kapko, D.V.; Brodach, M.M. Research of the microclimate formed by the local DCV. Energy Build. 2015, 90, 1–5. [Google Scholar] [CrossRef]
  97. Shao, X.; Li, X. Evaluating the potential of airflow patterns to maintain a non-uniform indoor environment. Renew. Energy 2015, 73, 99–108. [Google Scholar] [CrossRef]
  98. Antoun, S.; Ghaddar, N.; Ghali, K. Coaxial personalized ventilation system and window performance for human thermal comfort in asymmetrical environment. Energy Build. 2016, 111, 253–266. [Google Scholar] [CrossRef]
  99. Abou Hweij, W.; Ghaddar, N.; Ghali, K.; Habchi, C. Optimized performance of displacement ventilation aided with chair fans for comfort and indoor air quality. Energy Build. 2016, 127, 907–919. [Google Scholar] [CrossRef]
  100. El-Fil, B.; Ghaddar, N.; Ghali, K. Optimizing performance of ceiling-mounted personalized ventilation system assisted by chair fans: Assessment of thermal comfort and indoor air quality. Sci. Technol. Built Environ. 2016, 22, 412–430. [Google Scholar] [CrossRef]
  101. Zhu, S.; Cai, W.; Spengler, J.D. Control of sleep environment of an infant by wide-cover type personalized ventilation. Energy Build. 2016, 129, 69–80. [Google Scholar] [CrossRef]
  102. Conceição, E.Z.E.; Santiago, C.I.M.; Awbi, H.B. Numerical study of different ceiling-mounted air distribution systems for a virtual classroom environment. Indoor Built Environ. 2016, 26, 1382–1396. [Google Scholar] [CrossRef]
  103. Habchi, C.; Chakroun, W.; Alotaibi, S.; Ghali, K.; Ghaddar, N. Effect of shifts from occupant design position on performance of ceiling personalized ventilation assisted with desk fan or chair fans. Energy Build. 2016, 117, 20–32. [Google Scholar] [CrossRef]
  104. Habchi, C.; Ghali, K.; Ghaddar, N.; Chakroun, W.; Alotaibi, S. Ceiling personalized ventilation combined with desk fans for reduced direct and indirect cross-contamination and efficient use of office space. Energy Conv. Manag. 2016, 111, 158–173. [Google Scholar] [CrossRef]
  105. Mao, N.; Song, M.; Deng, S.; Pan, D.; Chen, S. Experimental and numerical study on air flow and moisture transport in sleeping environments with a task/ambient air conditioning (TAC) system. Energy Build. 2016, 133, 596–604. [Google Scholar] [CrossRef]
  106. Taheri, M.; Schuss, M.; Fail, A.; Mahdavi, A. A performance assessment of an office space with displacement, personal, and natural ventilation systems. Build. Simul. 2016, 9, 89–100. [Google Scholar] [CrossRef]
  107. Zhu, S.; Dalgo, D.; Srebric, J.; Kato, S. Cooling efficiency of a spot-type personalized air-conditioner. Build. Environ. 2017, 121, 35–48. [Google Scholar] [CrossRef]
  108. Conceição, E.Z.E.; Farinho, J.P.; Lúcio, M.M.J.R. Numerical and experimental study of personalized ventilation installed in a double occupation desk placed nearby a window subjected to solar radiation. Int. J. Vent. 2017, 17, 1–21. [Google Scholar] [CrossRef]
  109. Mao, N.; Pan, D.; Li, Z.; Xu, Y.; Song, M.; Deng, S. A numerical study on influences of building envelope heat gain on operating performances of a bed-based task/ambient air conditioning (TAC) system in energy saving and thermal comfort. Appl. Energy 2017, 192, 213–221. [Google Scholar] [CrossRef]
  110. Mao, N.; Pan, D.; Song, M.; Li, Z.; Xu, Y.; Deng, S. Operating optimization for improved energy consumption of a TAC system affected by nighttime thermal loads of building envelopes. Energy 2017, 133, 491–501. [Google Scholar] [CrossRef]
  111. Ahmed, A.Q.; Gao, S.; Kareem, A.K. Energy saving and indoor thermal comfort evaluation using a novel local exhaust ventilation system for office rooms. Appl. Therm. Eng. 2017, 110, 821–834. [Google Scholar] [CrossRef] [Green Version]
  112. Ahmed, A.Q.; Gao, S. Numerical investigation of height impact of local exhaust combined with an office work station on energy saving and indoor environment. Build. Environ. 2017, 122, 194–205. [Google Scholar] [CrossRef]
  113. Al Assaad, D.; Ghali, K.; Ghaddar, N.; Habchi, C. Mixing ventilation coupled with personalized sinusoidal ventilation: Optimal frequency and flow rate for acceptable air quality. Energy Build. 2017, 154, 569–580. [Google Scholar] [CrossRef]
  114. Du, J.; Chan, M.; Pan, D.; Deng, S. A numerical study on the effects of design/operating parameters of the radiant panel in a radiation-based task air conditioning system on indoor thermal comfort and energy saving for a sleeping environment. Energy Build. 2017, 151, 250–262. [Google Scholar] [CrossRef]
  115. Kong, M.; Dang, T.Q.; Zhang, J.; Khalifa, H.E. Micro-environmental control for efficient local cooling. Build. Environ. 2017, 118, 300–312. [Google Scholar] [CrossRef] [Green Version]
  116. Sun, L.; Tan, S.; Zhen, L. Working parameter selection and analysis of radiant floor cooling combined with desk-based task/ambient conditioning (TAC) system. Build. Serv. Eng. Res. Technol. 2017, 38, 104–120. [Google Scholar]
  117. Alotaibi, S.; Chakroun, W.; Habchi, C.; Ghali, K.; Ghaddar, N. Effectiveness of contaminant confinement in office spaces equipped with ceiling personalized ventilation system. Build. Simul. 2018, 11, 773–786. [Google Scholar] [CrossRef]
  118. Alsaad, H.; Voelker, C. Performance assessment of a ductless personalized ventilation system using a validated CFD model. J. Build. Perf. Simul. 2018, 11, 689–704. [Google Scholar] [CrossRef]
  119. Al Assaad, D.; Ghali, K.; Ghaddar, N. Effectiveness of intermittent personalized ventilation assisting a chilled ceiling for enhanced thermal comfort and acceptable indoor air quality. Build. Environ. 2018, 144, 9–22. [Google Scholar] [CrossRef]
  120. Al Assaad, D.; Habchi, C.; Ghali, K.; Ghaddar, N. Effectiveness of intermittent personalized ventilation in protecting occupant from indoor particles. Build. Environ. 2018, 128, 22–32. [Google Scholar] [CrossRef]
  121. Conceição, Z.E.; Santiago, I.C.; Lúcio, M.M.; Awbi, B.H. Predicting the Air Quality, Thermal Comfort and Draught Risk for a Virtual Classroom with Desk-Type Personalized Ventilation Systems. Buildings 2018, 8, 35. [Google Scholar] [CrossRef]
  122. Gao, R.; Wang, C.; Li, A.; Yu, S.; Deng, B. A novel targeted personalized ventilation system based on the shooting concept. Build. Environ. 2018, 135, 269–279. [Google Scholar] [CrossRef]
  123. Rahmati, B.; Mohammad Jadidi, A.; Heidarian, A. Investigation in performance of a hybrid under-floor air distribution with improved desk displacement ventilation system in a small office. Appl. Therm. Eng. 2018, 138, 861–872. [Google Scholar] [CrossRef]
  124. Sekhar, C.; Zheng, L. Study of an integrated personalized ventilation and local fan-induced active chilled beam air conditioning system in hot and humid climate. Build. Simul. 2018, 11, 787–801. [Google Scholar] [CrossRef]
  125. Zou, Y.; Zhao, X.; Chen, Q. Comparison of STAR-CCM+ and ANSYS Fluent for simulating indoor airflows. Build. Simul. 2018, 11, 165–174. [Google Scholar] [CrossRef]
  126. Greenshields, C.J. OpenFOAM User Guide. Version 3; Openfoam Foundation Ltd.: London, UK, 2015. [Google Scholar]
  127. Kaczmarczyk, J.; Melikov, A.; Bolashikov, Z.; Nikolaev, L.; Fanger, P.O. Human Response to Five Designs of Personalized Ventilation. HVAC&R Res. 2006, 12, 367–384. [Google Scholar]
  128. Liu, J.; Xie, X.; Qin, F.; Song, S.; Lv, D. A case study of ground source direct cooling system integrated with water storage tank system. Build. Simul. 2016, 9, 659–668. [Google Scholar] [CrossRef]
  129. Yakhot, V.; Orszag, S.A. Renormalization group analysis of turbulence. I. Basic theory. J. Sci. Comput. 1986, 1, 3–51. [Google Scholar] [CrossRef]
  130. Shih, T.-H.; Liou, W.W.; Shabbir, A.; Yang, Z.; Zhu, J. A new k-ϵ eddy viscosity model for high reynolds number turbulent flows. Comput. Fluids 1995, 24, 227–238. [Google Scholar] [CrossRef]
  131. Jones, W.P.; Launder, B.E. The prediction of laminarization with a two-equation model of turbulence. Int. J. Heat Mass Transf. 1972, 15, 301–314. [Google Scholar] [CrossRef]
  132. Menter, F.R. Two-equation eddy-viscosity turbulence models for engineering applications. AIAA J. 1994, 32, 1598–1605. [Google Scholar] [CrossRef] [Green Version]
  133. Chen, Q.; Xu, W. A zero-equation turbulence model for indoor airflow simulation. Energy Build. 1998, 28, 137–144. [Google Scholar] [CrossRef]
  134. Leschziner, M.A. Modelling engineering flows with Reynolds stress turbulence closure. J. Wind Eng. Ind. Aerodyn. 1990, 35, 21–47. [Google Scholar] [CrossRef]
  135. Yuce, B.E.; Pulat, E. Forced, natural and mixed convection benchmark studies for indoor thermal environments. Int. Commun. Heat Mass Transf. 2018, 92, 1–14. [Google Scholar] [CrossRef]
  136. Murakami, S.; Kato, S.; Zeng, J. Combined simulation of airflow, radiation and moisture transport for heat release from a human body. Build. Environ. 2000, 35, 489–500. [Google Scholar] [CrossRef]
  137. Zhu, S.; Kato, S.; Murakami, S.; Hayashi, T. Study on inhalation region by means of CFD analysis and experiment. Build. Environ. 2005, 40, 1329–1336. [Google Scholar] [CrossRef]
  138. Davidson, L.; Nielsen, P.; Sveningsson, A. Modification of the v2f model for computing the flow in a 3D wall jet. Turbulence Heat and Mass Transfer 2003, 4, 577–584. [Google Scholar]
  139. Wang, M.; Chen, Q. Assessment of Various Turbulence Models for Transitional Flows in an Enclosed Environment (RP-1271). HVAC&R Res. 2009, 15, 1099–1119. [Google Scholar]
  140. Berrouk, A.S.; Lai, A.C.K.; Cheung, A.C.T.; Wong, S.L. Experimental measurements and large eddy simulation of expiratory droplet dispersion in a mechanically ventilated enclosure with thermal effects. Build. Environ. 2010, 45, 371–379. [Google Scholar] [CrossRef]
  141. Taghinia, J.H.; Rahman, M.M.; Lu, X. Effects of different CFD modeling approaches and simplification of shape on prediction of flow field around manikin. Energy Build. 2018, 170, 47–60. [Google Scholar] [CrossRef]
  142. Taghinia, J.; Rahman, M.M.; Siikonen, T. Numerical simulation of airflow and temperature fields around an occupant in indoor environment. Energy Build. 2015, 104, 199–207. [Google Scholar] [CrossRef]
  143. ANSYS Inc. ANSYS Airpak User Guide; ANSYS Inc.: Canonsburg, PA, USA, 2007. [Google Scholar]
  144. Srebric, J.; Vukovic, V.; He, G.; Yang, X. CFD boundary conditions for contaminant dispersion, heat transfer and airflow simulations around human occupants in indoor environments. Build. Environ. 2008, 43, 294–303. [Google Scholar] [CrossRef]
  145. Li, X.; Yan, Y.; Tu, J. The simplification of computer simulated persons (CSPs) in CFD models of occupied indoor spaces. Build. Environ. 2015, 93, 155–164. [Google Scholar] [CrossRef]
  146. Yan, Y.; Li, X.; Tu, J. Effects of manikin model simplification on CFD predictions of thermal flow field around human bodies. Indoor Built Environ. 2017, 26, 1185–1197. [Google Scholar] [CrossRef]
  147. Yang, J.; Ni, S.; Weng, W. Modelling heat transfer and physiological responses of unclothed human body in hot environment by coupling CFD simulation with thermal model. Int. J. Therm. Sci. 2017, 120, 437–445. [Google Scholar] [CrossRef]
  148. ANSYS Inc. ANSYS FLUENT Theory Guide, Release 16.1; ANSYS Inc.: Canonsburg, PA, USA, 2014. [Google Scholar]
  149. Visser, D.C.; Houkema, M.; Siccama, N.B.; Komen, E.M.J. Validation of a FLUENT CFD model for hydrogen distribution in a containment. Nucl. Eng. Des. 2012, 245, 161–171. [Google Scholar] [CrossRef] [Green Version]
  150. Zhou, Y.; Wang, M.; Wang, M.; Wang, Y. Predictive accuracy of Boussinesq approximation in opposed mixed convection with a high-temperature heat source inside a building. Build. Environ. 2018, 144, 349–356. [Google Scholar] [CrossRef]
  151. CD-Adapco Inc. STAR-CCM+ 11.0 User Guide; CD-Adapco Inc.: Melville, NY, USA, 2016. [Google Scholar]
  152. ANSYS Inc. ANSYS CFX-Solver Theory Guide (Release 15.0); ANSYS Inc.: Canonsburg, PA, USA, 2013. [Google Scholar]
  153. Li, X.; Yan, Y.; Tu, J. Evaluation of models and methods to simulate thermal radiation in indoor spaces. Build. Environ. 2018, 144, 259–267. [Google Scholar] [CrossRef]
  154. Gilani, S.; Montazeri, H.; Blocken, B. CFD simulation of stratified indoor environment in displacement ventilation: Validation and sensitivity analysis. Build. Environ. 2016, 95, 299–313. [Google Scholar] [CrossRef]
  155. Zhu, S.; Kato, S.; Ooka, R.; Sakoi, T. Development of a Computational Thermal Manikin Applicable in a Nonuniform Thermal Environment—Part 1: Coupled Simulation of Convection, Radiation, and Smith’s Human Thermal Physiological Model for Sensible Heat Transfer from a Seated Human Body in Radiant Environment. HVAC&R Res. 2007, 13, 661–679. [Google Scholar]
  156. Zhu, S.; Kato, S.; Ooka, R.; Sakoi, T.; Tsuzuki, K. Development of a Computational Thermal Manikin Applicable in a Non-Uniform Thermal Environment—Part 2: Coupled Simulation Using Sakoi’s Human Thermal Physiological Model. HVAC&R Res. 2008, 14, 545–564. [Google Scholar]
  157. Russo, J. A Detailed and Systematic Investigation of Personal Ventilation Systems. Ph.D. Thesis, Syracuse University, Syracuse, NY, USA, 2011. [Google Scholar]
  158. Martinho, N.; Lopes, A.; Gameiro da Silva, M. Evaluation of errors on the CFD computation of air flow and heat transfer around the human body. Build. Environ. 2012, 58, 58–69. [Google Scholar] [CrossRef]
  159. Al-Othmani, M.; Ghaddar, N.; Ghali, K. A multi-segmented human bioheat model for transient and asymmetric radiative environments. Int. J. Heat Mass Transf. 2008, 51, 5522–5533. [Google Scholar] [CrossRef]
  160. Sideroff, C.N. Detailed Examinations of the Human Micro-Environment by CFD. Ph.D. Thesis, Syracuse University, Syracuse, NY, USA, 2009. [Google Scholar]
  161. Voelker, C.; Alsaad, H. Simulating the human body’s microclimate using automatic coupling of CFD and an advanced thermoregulation model. Indoor Air 2018, 28, 415–425. [Google Scholar] [CrossRef]
  162. Rim, D.; Novoselac, A. Transport of particulate and gaseous pollutants in the vicinity of a human body. Build. Environ. 2009, 44, 1840–1849. [Google Scholar] [CrossRef]
  163. Fu, M.; Weng, W.; Chen, W.; Luo, N. Review on modeling heat transfer and thermoregulatory responses in human body. J. Therm. Biol. 2016, 62, 189–200. [Google Scholar] [CrossRef]
  164. Zhang, H.; Huizenga, C.; Arens, E.; Yu, T. Considering individual physiological differences in a human thermal model. J. Therm. Biol. 2001, 26, 401–408. [Google Scholar] [CrossRef] [Green Version]
  165. Huizenga, C.; Hui, Z.; Arens, E. A model of human physiology and comfort for assessing complex thermal environments. Build. Environ. 2001, 36, 691–699. [Google Scholar] [CrossRef] [Green Version]
  166. Zhu, S.; Kato, S.; Yang, L. Examination of Thermal Adaptive Effect of Postural and Positional Adjustment of a Seated Human Body Exposed to Spot Airflow. ASHRAE Trans. 2007, 113, 300–312. [Google Scholar]
  167. Dixit, A.; Gade, U. A case study on human bio-heat transfer and thermal comfort within CFD. Build. Environ. 2015, 94, 122–130. [Google Scholar] [CrossRef]
  168. Cropper, P.C.; Yang, T.; Cook, M.; Fiala, D.; Yousaf, R. Coupling a model of human thermoregulation with computational fluid dynamics for predicting human–environment interaction. J. Build. Perf. Simul. 2010, 3, 233–243. [Google Scholar] [CrossRef]
  169. Katić, K.; Li, R.; Zeiler, W. Thermophysiological models and their applications: A review. Build. Environ. 2016, 106, 286–300. [Google Scholar] [CrossRef] [Green Version]
  170. Koelblen, B.; Psikuta, A.; Bogdan, A.; Annaheim, S.; Rossi, R.M. Human simulator—A tool for predicting thermal sensation in the built environment. Build. Environ. 2018, 143, 632–644. [Google Scholar] [CrossRef]
  171. Psikuta, A.; Allegrini, J.; Koelblen, B.; Bogdan, A.; Annaheim, S.; Martínez, N.; Derome, D.; Carmeliet, J.; Rossi, R.M. Thermal manikins controlled by human thermoregulation models for energy efficiency and thermal comfort research—A review. Renew. Sust. Energ. Rev. 2017, 78, 1315–1330. [Google Scholar] [CrossRef]
  172. Cheng, Y.; Niu, J.; Gao, N. Thermal comfort models: A review and numerical investigation. Build. Environ. 2012, 47, 13–22. [Google Scholar] [CrossRef]
  173. Wang, H.; Zhai, Z. Analyzing grid independency and numerical viscosity of computational fluid dynamics for indoor environment applications. Build. Environ. 2012, 52, 107–118. [Google Scholar] [CrossRef]
  174. ANSYS Inc. Gambit Version 2.4; ANSYS Inc.: Canonsburg, PA, USA, 2007. [Google Scholar]
  175. ANSYS Inc. ICEM CFD Guide, Release 16.1; ANSYS Inc.: Canonsburg, PA, USA, 2014. [Google Scholar]
  176. Roache, P.J. Perspective: A Method for Uniform Reporting of Grid Refinement Studies. J. Fluids Eng. 1994, 116, 405–413. [Google Scholar] [CrossRef]
  177. Nielsen, P.V.; Murakami, S.; Kato, S.; Topp, C.; Yang, J.-H. Benchmark Tests for a Computer Simulated Person; Aalborg University: Aalborg, Denmark; Indoor Environmental Engineering: San Francisco, CA, USA, 2003; Available online: https://Www.Cfd-Benchmarks.Com/ (accessed on 7 July 2019).
  178. Khalifa, H.E.; Janos, M.; Russo, J. Benchmark Tests for a Computer Simulated Person with Personal Ventilation; Aalborg University: Aalborg, Denmark; Indoor Environmental Engineering: San Francisco, CA, USA, 2008; Available online: https://Www.Cfd-Benchmarks.Com/ (accessed on 7 July 2019).
  179. Khalifa, H.E.; Janos, M.I.; Dannenhoffer, J.F. Experimental investigation of reduced-mixing personal ventilation jets. Build. Environ. 2009, 44, 1551–1558. [Google Scholar] [CrossRef]
  180. Bogdan, A.; Chludzińska, M.; Koelblen, B.; Ogłodzinski, K.; Młynarczyk, M. Sensitivity of personalized ventilation air terminal device geometry with regard to a human thermal sensation. Sci. Technol. Built Environ. 2017, 23, 1105–1112. [Google Scholar] [CrossRef]
  181. Zhang, T.; Lee, K.; Chen, Q. A simplified approach to describe complex diffusers in displacement ventilation for CFD simulations. Indoor Air 2009, 19, 255–267. [Google Scholar] [CrossRef]
  182. Srebric, J.; Chen, Q. Simplified numerical models for complex air supply diffusers. HVAC&R Res. 2002, 8, 277–294. [Google Scholar]
  183. Nielsen, P.V.; Allard, F.; Awbi, H.B.; Davidson, L.; Schälin, A. Computational Fluid Dynamics in Ventilation Design REHVA Guidebook No 10; Taylor & Francis: Didcot, UK, 2007. [Google Scholar]
  184. Bragança, P.; Sodjavi, K.; Meslem, A.; Serres, L. Airflow characteristics and thermal comfort generated by a multi-cone ceiling diffuser with and without inserted lobes. Build. Environ. 2016, 108, 143–158. [Google Scholar] [CrossRef]
  185. Awwad, A.; Mohamed, M.H.; Fatouh, M. Optimal design of a louver face ceiling diffuser using CFD to improve occupant’s thermal comfort. J. Build. Eng. 2017, 11, 134–157. [Google Scholar] [CrossRef]
  186. Yau, Y.H.; Poh, K.S.; Badarudin, A. A numerical airflow pattern study of a floor swirl diffuser for UFAD system. Energy Build. 2018, 158, 525–535. [Google Scholar] [CrossRef]
  187. Chludzińska, M.; Bogdan, A. The role of the front pattern shape in modelling personalized airflow and its capacity to affect human thermal comfort. Build. Environ. 2017, 126, 373–381. [Google Scholar] [CrossRef]
  188. Zhang, T.; Wang, S.; Sun, G.; Xu, L.; Takaoka, D. Flow impact of an air conditioner to portable air cleaning. Build. Environ. 2010, 45, 2047–2056. [Google Scholar] [CrossRef]
  189. Zhu, Y.; Luo, M.; Ouyang, Q.; Huang, L.; Cao, B. Dynamic characteristics and comfort assessment of airflows in indoor environments: A review. Build. Environ. 2015, 91, 5–14. [Google Scholar] [CrossRef]
  190. Huang, L.; Ouyang, Q.; Zhu, Y. Perceptible airflow fluctuation frequency and human thermal response. Build. Environ. 2012, 54, 14–19. [Google Scholar] [CrossRef]
  191. Yang, J.; Melikov, A.K.; Fanger, P.O.; Li, X.; Yan, Q. Impact of Personalized Ventilation on Human Response: Comparison between Constant and Fluctuating Airflows under Warm Condition. In Proceedings of the Roomvent-8th International Conference on Air Distribution in Rooms, Copenhagen, Denmark, 8–11 September 2002; pp. 305–308. [Google Scholar]
  192. Zhou, G.; Melikov, A.K. Equivalent frequency-a new parameter for description of frequency characteristics of airflow fluctuations. In Proceedings of the Roomvent-8th International Conference on Air Distribution in Rooms, Copenhagen, Denmark, 8–11 September 2002; pp. 357–363. [Google Scholar]
  193. Melikov, A.; Popiolek, Z. Comparison of different methods for the determination of dynamic characteristics of low velocity thermal anemometers. Meas. Sci. Technol. 2004, 15, 1709. [Google Scholar] [CrossRef]
  194. Uğursal, A.; Culp, C.H. The effect of temperature, metabolic rate and dynamic localized airflow on thermal comfort. Appl. Energy 2013, 111, 64–73. [Google Scholar] [CrossRef]
  195. Kato, S.; Yang, J.-H. Study on inhaled air quality in a personal air-conditioning environment using new scales of ventilation efficiency. Build. Environ. 2008, 43, 494–507. [Google Scholar] [CrossRef]
  196. Zhang, H.; Arens, E.; Huizenga, C.; Han, T. Thermal sensation and comfort models for non-uniform and transient environments: Part I: Local sensation of individual body parts. Build. Environ. 2010, 45, 380–388. [Google Scholar] [CrossRef] [Green Version]
  197. Zhang, H.; Arens, E.; Huizenga, C.; Han, T. Thermal sensation and comfort models for non-uniform and transient environments, part II: Local comfort of individual body parts. Build. Environ. 2010, 45, 389–398. [Google Scholar] [CrossRef] [Green Version]
  198. Zhang, H.; Arens, E.; Huizenga, C.; Han, T. Thermal sensation and comfort models for non-uniform and transient environments, part III: Whole-body sensation and comfort. Build. Environ. 2010, 45, 399–410. [Google Scholar] [CrossRef] [Green Version]
  199. Zhang, H. Human Thermal Sensation and Comfort in Transient and Non-Uniform Thermal Environments; University of California: Berkeley, CA, USA, 2003. [Google Scholar]
  200. Tanabe, S.; Arens, E.A.; Bauman, F.; Zhang, H.; Madsen, T. Evaluating thermal environments by using a thermal manikin with controlled skin surface temperature. ASHRAE Trans. 1994, 100, 39–48. [Google Scholar]
  201. Nilsson, H.; Holmér, I. Comfort climate evaluation with thermal manikin methods and computer simulation models. Indoor Air 2003, 13, 28. [Google Scholar] [CrossRef]
  202. ISO/TS. Ergonomics of the Thermal Environment-Evaluation of Thermal Environments in Vehicles; ISO/TS 14505; ISO: Geneva, Switzerland, 2007. [Google Scholar]
  203. ASHRAE Standard 55. Thermal Environmental Conditions for Human Occupancy; American Society of Heating, Refrigerating and Air-Conditioning Engineers: Atlanta, GA, USA, 2004. [Google Scholar]
  204. ISO Standard 7730. Ergonomics of the Thermal Environment e Analytical Determination and Interpretation of Thermal Comfort Using Calculation of the PMV and PPD Indices and Local Thermal Comfort Criteria; International Organization for Standardization: Geneva, Switzerland, 2005. [Google Scholar]
  205. Melikov, A.; Ivanova, T.; Stefanova, G. Seat headrest-incorporated personalized ventilation: Thermal comfort and inhaled air quality. Build. Environ. 2012, 47, 100–108. [Google Scholar] [CrossRef]
  206. Melikov, A.K. Human body micro-environment: The benefits of controlling airflow interaction. Build. Environ. 2015, 91, 70–77. [Google Scholar] [CrossRef]
  207. Al Assaad, D.; Habchi, C.; Ghali, K.; Ghaddar, N. Simplified model for thermal comfort, IAQ and energy savings in rooms conditioned by displacement ventilation aided with transient personalized ventilation. Energy Conv. Manag. 2018, 162, 203–217. [Google Scholar] [CrossRef]
  208. Cheng, Y.; Niu, J.; Gao, N. Stratified air distribution systems in a large lecture theatre: A numerical method to optimize thermal comfort and maximize energy saving. Energy Build. 2012, 55, 515–525. [Google Scholar] [CrossRef]
Figure 1. Distribution of reviewed journal papers on PV studies using CFD. (a) Year of publication; (b) CFD code.
Figure 1. Distribution of reviewed journal papers on PV studies using CFD. (a) Year of publication; (b) CFD code.
Sustainability 11 04166 g001
Figure 2. Commonly used PV types. (a) ATDs (modified from [6]), (b) CBPV [29,37], (c) CMPV [88,89] and DMPV [120], (d) DF and CF [103,104].
Figure 2. Commonly used PV types. (a) ATDs (modified from [6]), (b) CBPV [29,37], (c) CMPV [88,89] and DMPV [120], (d) DF and CF [103,104].
Sustainability 11 04166 g002
Figure 3. Grid distribution in room and around the human body. (a) Hybrid grids (prismatic, tetrahedral, and hexahedral) [69], (b) Hybrid grids (tetrahedral and prismatic) [117], (c) Hybrid grids (honeycomb and hexahedral) [115].
Figure 3. Grid distribution in room and around the human body. (a) Hybrid grids (prismatic, tetrahedral, and hexahedral) [69], (b) Hybrid grids (tetrahedral and prismatic) [117], (c) Hybrid grids (honeycomb and hexahedral) [115].
Sustainability 11 04166 g003
Figure 4. (a) UFAD and PV configuration view, and (b) velocity distribution [39].
Figure 4. (a) UFAD and PV configuration view, and (b) velocity distribution [39].
Sustainability 11 04166 g004
Figure 5. (a) Low mixing co-flow nozzle [179] and (b) low-mixing ceiling-mounted nozzle [88].
Figure 5. (a) Low mixing co-flow nozzle [179] and (b) low-mixing ceiling-mounted nozzle [88].
Sustainability 11 04166 g005
Table 1. CFD analyses of PV in indoor built environments.
Table 1. CFD analyses of PV in indoor built environments.
Author & YearPV with BVCFD CodeTurb.Rad.BC (no.)NWTTarget Parameters
Gao and Niu 2004 [36]PV with DVNSSKENOTsk(1)SWFAT, AV, CHTC, CO2, PER, PPD, PUE
Gao and Niu 2005 [69]PV with DVFluentRNGNSTsk(1)EWTAT, AV, CHTC, CO2, PER
Gao et al. 2006 [37]DMPV or CBPV with DV, RMPFluentRNGYESTsk(16)EWTACE, AT, AV, AVV, CHTC, LTC, LTS, OTC, OTS, PER
Nielsen et al. 2007 [70]PV with MVNSk-εNSq(1)NSACH, AV, AVV, PEI
Gao et al. 2007 [71]PV with DV or MVFluentSKEYESTsk(16)NSAT, AV, AVV, Cinh, LTC, LTS, OTC, OTS, PER
Zhao and Guan 2007 [72]PVSTACH-3ZEQNSq(1)NoAT, AV, AVV, PC
Yang and Sekhar 2008 [73]PV or CMPV with MVFluentSKENONONSAV, AVV, FAP, PEE
Russo et al. 2009 [74]PV with UFADANSYS FluentRKENSTsk(20)EWTAT, AV, PEE, TI
Dygert et al. 2009 [75]PV with UFADFluentSKE, RKE, RNG, SKW, SSTNOTsk(1)EWTAQI, AV, AT
Russo and Ezzat Khalifa 2010 [76]PVNSNSNSTsk(1)EWTIF, PEE, Sct
Russo and Khalifa 2010 [77]PVFluentRKENSTsk(1)EWTIF, Ozone
Conceição et al. 2010 [78]DMPVNSk-εNSTsk(15)NSAT, AV, DR, PMV
Tham and Pantelic 2010 [79]DMPV and DF with MVFluentSKENOTsk(26)EWTAV, CHF, PEE, Teq
He et al. 2011 [80]RMP with DV, MV, UFADFluentRNGS2Sq(1)SWFAV, AVV, IF, PC, Rc, TGC
Zhai and Metzger 2011 [81]PV with MVPHOENICSRNGNONSNSAGE, ES, PMV
Mazej and Butala 2011 [82]MSDP with UFADANSYS FluentRNGS2STsk(16)EWTAV, AT, AVV, CHTC, CO2, PEE, RHTC, RIE, Teq, TGC
Adamu et al. 2011 [83]NPVPHOENICSRNGNSQ(1)NSAGE, AT, ATT, AV, CRE, LACH, MACE
Ishiguro et al. 2011 [84]CMPVNSNSNOQ(1)NSAV, AT, AVV, PMV, Qc
Russo and Khalifa 2011 [85]PVNSNSNSTsk(1)NSAV, AT, AQI, IF, TGC
Li et al. 2012 [86]CBPV or DMPV with MV or DVFluentRNGNSq(1)SWFACE, AT, AV, AVV, IF, PC
Kanaan et al. 2012 [87]PV with DVAirpakSKEYESQ(1)NSAT, AV, AVV, CO2, PEE
Shen et al. 2013 [38]RMP or VDG with DV, MV or UFADANSYS FluentSKES2Sq(1)SWFAT, AV, AVV, MFP, TGC, VE
Makhoul et al. 2013 [88]CMPV and DFANSYS FluentRKENSTsk(11)EWTAV, AT, AVV, CO2, ES, LTC, OTC, PEE
Makhoul et al. 2013 [89]CMPVANSYS FluentRKENSTsk(11)EWTLTC, LTS, OTC, OTS, PEE, Qc
Makhoul et al. 2013 [90]CMPVANSYS FluentRKENSTsk(-)NSAT, IF, PC, PVV
Yang et al. 2013 [91]DMPV or VDG with MVANSYS FluentSKENSTsk(1)EWTCO2, IF, PEE
Cheong and Huang 2013 [92]RMP or DMPV with DVFluentRNGNSTsk(26)NSAV, AT, AVV, PEIc
Makhoul et al. 2013 [93]CMPVANSYS FluentRKENSTsk(11)EWTAV, AT, CO2, PEE
Yang and Sekhar 2014 [94]CMPV with MVFluentSKENONONSAV, AVV, PEE, PAP, TGC
Russo and Ezzat Khalifa 2014 [95]PV with UFADFluentRKENOTsk(1)EWTOzone, IF
Kong et al. 2015 [39]PV with UFADSTAR-CCM+RKENSTsk(20)TLWTAT, AV, AVV, PEE, Teq, TGC
Naumov et al. 2015 [96]PVANSYSNSNSTsk(-)NSAV, AT, EXT, ES
Shao and Li 2015[97]PVAirpakZEQNSQ(1)NODPSA, TACS, TAIC, TASA
Antoun et al. 2016 [98]CMPVANSYS FluentRKES2STsk(11)EWTAV, LTC, LTS, OTS, OTC, Qc, RHF, TI, Tsk
Abou Hweij et al. 2016 [99]CF with DVANSYS FluentSKENSTsk(11)EWTAT, AV, LTC, LTS, OTC, OTS, Tsk, PEE
El-Fil et al. 2016 [100]CMPV and CFANSYF FluentRKENOTsk(11)EWTAT, AV, CO2, LTC, LTS, OTC, OTS, PEE, Qc
Zhu et al. 2016 [101]WCPV with MVSTAR-CDRNGNOq(1)&Tsk(1)SWFAT, AV, DR, PEE, SVE3, SVE3*, SVE4
Conceição et al. 2016 [102]CMPVNSRNGYESTsk(25)NSAQN, AT, AV, CO2, ADI, DR, PPD, Tsk
Habchi et al. 2016 [103]CMPV and CF or DFANSYS FluentRKENOTsk(11)EWTAT, AV, CO2, DFr, IF, LTC, OTC, PC, PEE, Qc
Habchi et al. 2016 [104]CMPV and DFANSYS FluentRKENOq(20)EWTAT, AV, DFr, IF, PC, Qc
Mao et al. 2016 [105]BTACANSYS FluentSSTS2STsk(16)NSAV, AT, AVV, RH, Tsk
Taheri et al. 2016 [106]PV with DVNSNSNSNSNSAT, AV, CO2, PMV, TSV
Zhu et al. 2017 [107]RPAC with MVFluentRNGYESTsk(17)EWTAT, AV, AVV, CE, DR, OTS, Tsk, ΔPMV
Conceição et al. 2017 [108]DMPVNSRNGNSTsk(15)NSAV, AT, CO2, DR, MRT, PMV, PPD, Tclo, Tsk, TI
Mao et al. 2017 [109]BTACANSYS FluentSSTS2STsk(1)NSAT, AV, PMV, Qc
Mao et al. 2017 [110]BTACANSYS FluentSSTS2STsk(1)NSAT, AV, PMV, Qc
Ahmed et al. 2017 [111]LEVO with DVANSYS FluentRNGDOq(1)EWTAV, AT, PC, PMV, PPD, Qc, VATD
Ahmed and Gao 2017 [112]LEVO with DVANSYS FluentRNGDOq(1)EWTAV, AT, DR, PC, Qc, VTAD
Al Assaad et al. 2017 [113]PV with MVANSYS FluentRNGNSTsk(11)EWTAT, AV, CO2, FRF, OTC, PEE, Qc, Tsk
Du et al. 2017 [114]BTAC with radiant panelANSYS FluentSSTS2STsk(1)NSAT, AV, EUC, DR, PMV, To
Kong et al. 2017 [115]RMP or DMPV with MVSTAR-CCM+SKEYESTsk(20)TLWTAT, AV, CHF, CHTC, To, Tsk
Sun et al. 2017 [116]DMPV with cooling ceilingAirpakSKEDOQ(1)NSAT, DR, PMV, PPD, To, ES
Alotaibi et al. 2018 [117]CMPV and DF or CFANSYS FluentRKENSq(1)NSAV, AT, CI, PC, Qc, ΔC
Alsaad and Voelker 2018 [118]DPV with DVANSYS FluentSKE, RKE, RNGS2STsk(16)EWTAT, AV, CO2, LTC, LTS, OTC, OTS, PEIc,
Al Assaad et al. 2018 [119]PV with cooling ceilingANSYS FluentRNGS2STsk(10)EWTAT, CO2, OTC, OTS, PEE, Qc, TI
Al Assaad et al. 2018 [120]DMPV with MVANSYS FluentRNGNSq(1)EWTAT, AV, CO2, DFr, FRF, IF, PC
Conceição et al. 2018 [121]DMPVNSRNGNSTsk(25)NSADI, AQN, AT, AV, CO2, DR, PPD
Gao et al. 2018 [122]TPVNSRSMNSNSNSAV, AT, Qc, Ts, Tr
Rahmat et al. 2018 [123]DF and DMPV with UFADAirpakZEQS2SQ(1)NOAT, AV, PMV, PPD, Qc, VATD
Sekhar and Zheng 2018 [124]PV with ACBANSYS FluentRKENSq(1)SWFAT, AV, AVV, PPD, PMV, Qc,
Abbreviations: PV type with background ventilation type (PV with BV), new wall treatment (NWT), radiation model (Rad.), turbulence model (Turb.), CTM boundary condition type with segment number (BC (no.)); PV with BV: bed-based task/ambient air conditioning (BTAC), chair fan (CF), chair based personalized ventilation (CBPV), ceiling mounted personalized ventilation (CMPV), desk mounted personalized ventilation (DMPV), desk fan (DF), ductless personalized ventilation (DPV), local exhaust ventilation for offices (LEVO), movable slot diffuser plate (MSDP), natural personalized ventilation (NPV), robotic personal air conditioning (RPAC), round movable panel (RMP), targeted personalized ventilation (TPV), vertical desk grille (VDG), wide-cover personalized ventilation (WCPV); Background HVAC system type: displacement ventilation (DV), mixing ventilation (MV), underfloor air distribution (UFAD), local fan-induced active chilled beam air conditioning system (ACB); Turbulence modes: low Reynolds number k-ε (LRNKE), realizable k-ε (RKE), Reynolds stress models (RSM), renormalization group k-ε (RNG), SST k-ω (SST), standard k-ε (SKE), zero-equation turbulence model (ZEQ); CFD codes: not specified (NS); Near wall treatment: enhanced wall treatment (EWT); standard wall function (SWF), two-Layer all y+ wall treatment (TLWT); Radiation model: discrete ordinates radiation model (DO), surface-to-surface radiation model (S2S); Target parameters: (i) local ventilation and air quality: air change rate (h−1) (ACH), carbon dioxide concentration (ppm) (CO2), concentration asymmetry (-) (ΔC), deposited fraction (DFr), inhaled contaminant concentration (ppm) (Cinh), intake fraction (IF), particle concentration (kg/m3) (PC), re-inhaled exposure index (-) (RIE), scale for ventilation efficiency 3 (-) (SVE3), new scale for ventilation efficiency 3 (-) (SVE3*), scale for ventilation efficiency 4 (-) (SVE4), tracer gas concentration (ppm) (TGC), personal exposure index (-) (PEI), personal exposure effectiveness (-) (PEE); (ii) Thermal comfort: change of predicted mean vote (-) (ΔPMV), local thermal comfort (-) (LTC), local thermal sensation (-) (LTS), overall thermal comfort (-) (OTC); overall thermal sensation (-) (OTS), predicted percentage dissatisfied (%) (PPD), draught risk (%) (DR), predicted mean vote (-) (PMV), CTM-based equivalent temperature (°C) (Teq), vertical air temperature difference at the head and foot level (°C) (VATD), mean radiant temperature (°C) (MRT), operative temperature (°C) (To); (iii) Energy saving: cooling efficiency (-) (CE), energy utilization coefficient (-) (EUC), calculated energy consumption (W) (Qc); Others: age of air (s) (AGE), air distribution index (-) (ADI), air quality number (-) (AQN), air temperature (°C) (AT), air velocity (m/s) (AV), air velocity vector (-) (AVV), air turnover time (s) (ATT), concentration uniformity index (-) (Rc), convective heat flux (W/m2) (CHF), confinement index (%) (CI), convective heat transfer coefficient (W/m2K) (CHTC), contaminant removal efficiency (CRE), cloth surface temperature (°C) (Tclo), derivation of the index of difference potential by supply air (-) (DPSA), excess temperature (°C) (EXT), flow rate frequency (Hz) (FRF), fresh air percentage (%) (FAP), local air change index (%) (LACH), mean air exchange efficiency (%) (MACE), mass fraction of pollutants (%) (MFP), ozone concentration (ppm) (Ozone), particle velocity vector (-) (PVV), radiative heat flux (W/m2) (RHF), radiative heat transfer coefficient (W/m2K) (RHTC), relative humidity (%) (RH), skin surface temperature (°C) (Tsk), transient accessibility of contaminant source (-) (TACS), transient accessibility of initial condition (-) (TAIC), transient accessibility of supply air (-) (TASA), turbulent Schmidt number (-) (Sct).
Table 2. Grid distribution in the computational domain and around CTMs
Table 2. Grid distribution in the computational domain and around CTMs
Authors and YearCDDTCN(LCN)TCNPFLH(NBL)MSCSy+MSAGrid Type
Gao and Niu 2004 [36]2.6 × 2.2 × 2.71.81(1.59)0.12 1.57Hex/Tet
Gao and Niu 2005 [69]2.6 × 2.2 × 2.71.75(1.06)0.11 most < 11.59Hex/Tet
Russo, Dang et al. 2009 [74]2.0 × 2.6 × 2.54.20.321.5(4)<12<3 Hex/Tet
Tham and Pantelic 2010 [79]5.55 × 3.7 × 2.62.64(2.08)0.052(15) <1 Hex/Tet
Shen, Gao et al. 2013 [38]5.4 × 4.8 × 2.61.080.02
Makhoul, Ghali et al. 2013 [88]1.7 × 3.4 × 2.8(half)1.20.071.5(4)100.8–41.78
Makhoul, Ghali et al. 2013 [89]1.7 × 3.4 × 2.8(half)1.210.071.5(4)200.8–41.78
Makhoul, Ghali et al. 2013 [90]1.7 × 3.4 × 2.6(half)1.210.08 --
Makhoul, Ghali et al. 2013 [93]1.7 × 3.4 × 2.6(half)1.210.081.5 0.8–41.78
Cheong and Huang 2013 [92]6.6 × 3.7 × 2.72.71(2.16)0.04 most < 11.57Hex/Tet
Yang, Sekhar et al. 2013 [91]6.6 × 3.7 × 2.64.610.072.2(2) 0.5–4
Kong, Zhang et al. 2015 [39]1.9 × 1.8 × 1.70.340.06(5)501.58 Hon
Antoun, Ghaddar et al. 2016 [98]1.7 × 2.7 × 2.6 (half)1.420.12(3) 0.8–4
Abou Hweij, Ghaddar et al. 2016 [99]2.5 × 2.75 × 2.83.520.18 0.8–4
El-Fil, Ghaddar et al. 2016 [100]3.4 × 3.4 × 2.82.830.091.8(3)200.9–4.5
Al Assaad, Ghali et al. 2017 [113]2.8 × 2.75 × 2.51.060.061.5
Ahmed, Gao et al. 2017 [111]4.0 × 2.7 × 3.02.750.081.5 0.7–4.5
Ahmed and Gao 2017 [112]4.0 × 2.7 × 3.02.750.081.5 0.7–4.5
Du, Chan et al. 2017 [114]3.7 × 2.6 × 3.292.44(1.45)0.080.4 Hex/Tet
Kong, Dang et al. 2017 [115]4.8 × 3.66 × 3.055.00.091(10)2.5-10<0.81.8Hon
Alsaad and Voelker 2018 [118]3.0 × 3.0 × 2.45.770.27 most < 1
Al Assaad, Habchi et al. 2018 [120]3.4 × 3.4 × 2.86.50.20 150.8–4 Tet
Al Assaad, Ghali et al. 2018 [119]2.5 × 2.75 × 2.82.110.11 150.8–4 Tet
Alotaibi, Chakroun et al. 2018 [117]3.4 × 3.4 × 2.62.840.09
Gao, Wang et al. 2018 [122]3 × 1.5 × 3 (half)1.60.12
Abbreviations: computational domain dimension (m × m × m) (CDD), total cells number (-) (TCN), local zone cell number around the human body (-) (LCN), first layer height (mm) (FLH), number of inflated prismatic layers (-) (NBL), manikin surface cell size (mm) (MSCS), non-dimensional wall distance in manikin surface (-) (y+), manikin surface area (m2) (MSA), grid type (GY), hexahedral grid (Hex), tetrahedral grid (Tet), honeycomb grid (Hon), total cell number per volume unit (TCNP).
Table 3. Summary of air supply parameters in a typical PV and background ventilation system.
Table 3. Summary of air supply parameters in a typical PV and background ventilation system.
Authors and YearPV SystemBackground Ventilation SystemDistance (m)
Flow Rate (L/s)Tin (°C)Ti (%)Flow Rate (L/s)Tin (°C)Ti (%)
Gao and Niu 2004 [36]PV0~3200.5DV14225<0.1
Kong, Zhang et al. 2015 [39]PV7.1, 9.421.7 UFAD48.824.3
Gao and Niu 2005 [69]PV0.8~1.618~245,10,20DV142210<0.1
Gao, Zhang et al. 2007 [71]PV10~2017~2120MV/DV/31~5117~21MV30, DV15
Russo, Dang et al. 2009 [74]Pri0.6~4.8, Se1.7~13.4Pri23.5Pri/Se1.7~10UFAD0.7~16.621 0.4
Russo and Khalifa 2010 [76,77]PV2.4, Sec6.7211.7UFAD10.1211.70.4
He, Niu et al. 2011 [80]RMP7~1520 DV/MV/UFAD25~4020 0.586
Mazej and Butala 2011 [82]PV5,1020,235DV12020,23,2610
Russo and Khalifa 2011 [85]PV2.4, Sec6.7211.7UFAD10.1211.70.4
Li, Niu et al. 2012 [86]CBPV/DMPV0.8~6.518 DV/MV22~2820
Kanaan, Ghaddar et al. 2012 [87]PV4~1018~22 DV11618 0.3~0.5
Cheong and Huang 2013 [92]PV5, 10, 15, 2019~235, 10DV602310
Makhoul, Ghali et al. 2013 [89,93]Pri5~10, Sec10~20Pri16~24, Sec26~28 Dif10~20, 40~57162.71.49
Shen, Gao et al. 2013 [38]RMP/VDG7,151710DV/MV/UFAD10~2317DV5, MV/UFAD15
Yang, Sekhar et al. 2013 [91]PV4, 82310MV(-)23100.15~0.55
Antoun, Ghaddar et al. 2016 [98]Pri8.5, Sec(-)16, Sec(-)3Dif50, 55, 6016, 18, 1921.6
El-Fil, Ghaddar et al. 2016 [100]CMPV8.5, CF2.516, CF(-)2.5Dif50162.51.4
Habchi, Chakroun et al. 2016 [103]Pri8.5,11, CF10, DF1016 Dif3516 1.4
Habchi, Ghali et al. 2016 [104]Pri10, DF5,10,1516 Dif3516 1.4 m
Zhu, Cai et al. 2016 [101]WCPV2.5,518,1910AC242410
Al Assaad, Ghali et al. 2017 [113]PV3.5,5,7.522 MV6320 0.4
Kong, Dang et al. 2017 [115]PV8~3821~2510MV7527.8, 28.3100.2~0.61
Zhu, Dalgo et al. 2017 [107]RPAC11.8~5919~24.610MV52.52610
Alsaad and Voelker 2018 [118]DPV5,6.5 DV16,24,4319,22 0.4
Alotaibi, Chakroun et al. 2018 [117]CMPV8.5, DF10, CF10162.5Dif35162.51.4
Al Assaad, Ghali et al. 2018 [119]PV3.5,5,7.522 MV6315 0.4
Abbreviations: personalized ventilation (PV), ceiling mounted personalized ventilation (CMPV), desk mounted personalized ventilation (DMPV), chair based personalized ventilation (CBPV), desk fan (DF), chair fan (CF), round movable panel (RMP), vertical desk grille (VDG), ductless personalized ventilation (DPV), robotic personal air conditioning (RPAC), wide-cover personalized ventilation (WCPV), displacement ventilation (DV), mixing ventilation (MV), underfloor air distribution (UFAD), primary nozzle (Pri), secondary nozzle (Sec), diffuser in the CMPV (Dif).
Table 4. Evaluation indices for ventilation effectiveness and inhaled air quality
Table 4. Evaluation indices for ventilation effectiveness and inhaled air quality
IndexDefinitionExplanation for ParametersPurpose of the Evaluation IndexTypical References
EV1 ε E V = C e C S C j C S Ce tracer gas concentration at the exhaust duct
CS tracer gas concentration at the supply duct
Cj tracer gas concentration at the simulated point (e.g., mouth of the exposure manikin)
Assess the air distribution efficiency in rooms and around a human body[80,118]
EV2 ε E V 2 = C p C S C e C S Cp pollutant concentration at simulated point (inhalation zone)
Ce pollutant concentration at the exhaust opening
CS pollutant concentration at the supply air
Evaluate the effect of PV on inhaled air quality and trace gas pollutant transport characteristics around a polluting occupant[38]
EV3 ε E V 3 = C r C b C r C f Cr tracer gas concentration of the recirculated air or at outlet
Cf tracer gas concentration of the fresh air
Cb tracer gas concentration at the breathing zone defined as a small sphere of 1 cm radius, at 2.5 cm from the occupant nose
Assess the mixing level of recirculated air and delivered fresh air supplied by second nozzle and diffuser and primary nozzle, respectively, using the tracer gas[88,103,113,119]
PEI ε P E I = C j C S C s f C S Cj pollutant concentration at the simulated point
CS pollutant concentration at the supply air
Csf pollutant concentration on the surface of pollutant source
Assess the extent of the measuring location affected by the polluting source[92]
PER ε P E R = C L C a C f C a CL tracer gas concentration of inhaled air
Cf tracer gas concentration of personalized air
Ca tracer gas concentration of ambient air.
Expressed as the fraction of personalized air in inhaled air[36,37,69,71]
PEE ε P E E = C I , 0 C I C I , 0 C P V CI,0 concentration of pollution in the inhaled air without PV
CI concentration of pollution in the inhaled air
CPV concentration of pollution in personalized air
Estimate personal exposure effectiveness for PV in inhaled air, C P V = 0 means that PV provides clean air with no pollutants[82,87,91,101]
PEEc ε P E E c = C R , S F 6 C I , S F 6 C R , S F 6 C P V , S F 6 CRS6 SF6 concentration in the re-circulated air duct
CPV,SF6 SF6 concentration of the supplying fresh air
CI,SF6 SF6 concentration at the breathing zone
Estimate personal exposure effectiveness for SF6 in inhaled air [73,94]
AQI A Q I = C b C e C P r i C e Cb tracer gas concentration at a point in the breathing zone
Cpri tracer gas concentration at the primary nozzle exit
Ce tracer gas concentration in the exhaust
Evaluate inhaled air quality in the floor diffuser and the secondary nozzle that supplied tracer gas, while primary air was kept free of tracer gas[39,74,76,85]
IF I F = M b r M g e Mbr Particle concentration at breathing level of healthy person, or pollutant mass exposed to person
Mge Particle generation concentration, or pollutant mass emitted from a source/polluting person
Represent the proportion of contaminants generated by the infected occupant that is inhaled by the exposed occupant[77,85,95,103,104,120]
DFr D F r = N b r N g e Nbr Number of particles deposited at the vicinity of the exposed person
Nge Number of particles generated by the infected person
Assess the cross-infection by the monitoring the particles deposited at the human body of the exposed person and other facilities in the room[103,104,120]
Abbreviations: personal exposure effectiveness (PEE), pollutant exposure index (PEI), ventilation effectiveness (VE1, VE2, or VE3), intake fraction (IF), deposited fraction (DFr), air quality index (AQI), pollutant exposure reduction (PER).

Share and Cite

MDPI and ACS Style

Liu, J.; Zhu, S.; Kim, M.K.; Srebric, J. A Review of CFD Analysis Methods for Personalized Ventilation (PV) in Indoor Built Environments. Sustainability 2019, 11, 4166. https://doi.org/10.3390/su11154166

AMA Style

Liu J, Zhu S, Kim MK, Srebric J. A Review of CFD Analysis Methods for Personalized Ventilation (PV) in Indoor Built Environments. Sustainability. 2019; 11(15):4166. https://doi.org/10.3390/su11154166

Chicago/Turabian Style

Liu, Jiying, Shengwei Zhu, Moon Keun Kim, and Jelena Srebric. 2019. "A Review of CFD Analysis Methods for Personalized Ventilation (PV) in Indoor Built Environments" Sustainability 11, no. 15: 4166. https://doi.org/10.3390/su11154166

APA Style

Liu, J., Zhu, S., Kim, M. K., & Srebric, J. (2019). A Review of CFD Analysis Methods for Personalized Ventilation (PV) in Indoor Built Environments. Sustainability, 11(15), 4166. https://doi.org/10.3390/su11154166

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop