Chamber Simulation Research Articles

Optimizing calibration settings for accurate water equivalent path length assessment using flat panel proton radiography

Objective: Proton range uncertainties can compromise the effectiveness of proton therapy treatments. Water equivalent path length (WEPL) assessment by flat panel detector proton radiography (FP-PR) can provide means of range uncertainty detection. Since WEPL accuracy intrinsically relies on the FP-PR calibration parameters, the purpose of this study is to establish an optimal calibration procedure that ensures high accuracy of WEPL measurements. To that end, several calibration settings were investigated. Approach: FP-PR calibration datasets were obtained simulating PR fields with different proton energies, directed towards water-equivalent material slabs of increasing thickness. The parameters investigated were the spacing between energy layers (ΔE) and the increment in thickness of the water-equivalent material slabs (ΔX) used for calibration. 30 calibrations were simulated, as a result of combining ΔE = 9, 7, 5, 3, 1 MeV and ΔX = 10, 8, 5, 3, 2, 1 mm. FP-PRs through a CIRS electron density phantom were simulated, and WEPL images corresponding to each calibration were obtained. Ground truth WEPL values were provided by range probing multi-layer ionization chamber simulations on each insert of the phantom. Relative WEPL errors between FP-PR simulations and ground truth were calculated for each insert. Mean relative WEPL errors and standard deviations across all inserts were computed for WEPL images obtained with each calibration. Main results: Large mean and standard deviations were found in WEPL images obtained with large ΔE values (ΔE = 9 or 7 MeV), for any ΔX. WEPL images obtained with ΔE ≤ 5 MeV and ΔX ≤ 5 mm resulted in a WEPL accuracy with mean values within ±0.5% and standard deviations around 1%. Significance: An optimal FP calibration in the framework of this study was established, characterized by 3 MeV ≤ ΔE ≤ 5 MeV and 2 mm ≤ ΔX ≤ 5 mm. Within these boundaries, highly accurate WEPL acquisitions using FP-PR are feasible and practical, holding the potential to assist future online range verification quality control procedures.

Development and implementation of a time-based signal generation scheme for the muon chamber simulation of the CBM experiment at FAIR

The Compressed Baryonic Matter (CBM) experiment in the upcoming Facility for Antiproton and Ion Research (FAIR), designed to take data in nuclear collisions at very high interaction rates of up to 10 MHz, will employ a free-streaming data acquisition with self-triggered readout electronics, without any hardware trigger. A simulation framework with a realistic digitization of the detectors in the muon chamber (MuCh) subsystem in CBM has been developed to provide a realistic simulation of the time-stamped data stream. In this article, we describe the implementation of the free-streaming detector simulation and the basic data related effects on the detector with respect to the interaction rate.

Response of the peatland carbon dioxide sink function to future climate change scenarios and water level management.

Stress factors such as climate change and drought may switch the role of temperate peatlands from carbon dioxide (CO2 ) sinks to sources, leading to positive feedback to global climate change. Water level management has been regarded as an important climate change mitigation strategy as it can sustain the natural net CO2 sink function of a peatland. Little is known about how resilient peatlands are in the face of future climate change scenarios, as well as how effectively water level management can sustain the CO2 sink function to mitigate global warming. The authors assess the effect of climate change on CO2 exchange of south Swedish temperate peatlands, which were either unmanaged or subject to water level regulation. Climate chamber simulations were conducted using experimental peatland mesocosms exposed to current and future representative concentration pathway (RCP) climate scenarios (RCP 2.6, 4.5 and 8.5). The results showed that all managed and unmanaged systems under future climate scenarios could serve as CO2 sinks throughout the experimental period. However, the 2018 extreme drought caused the unmanaged mesocosms under the RCP 4.5 and RCP 8.5 switch from a net CO2 sink to a source during summer. Surprisingly, the unmanaged mesocosms under RCP 2.6 benefited from the warmer climate, and served as the best sink among the other unmanaged systems. Water level management had the greatest effect on the CO2 sink function under RCP 8.5 and RCP 4.5, which improved their CO2 sink capability up to six and two times, respectively. Under the current climate scenario, water level management had a negative effect on the CO2 sink function, and it had almost no effect under RCP 2.6. Therefore, the researchers conclude that water level management is necessary for RCP 8.5, beneficial for RCP 4.5 and unimportant for RCP 2.6 and the current climate.

Highly oxygenated organic molecule (HOM) formation in the isoprene oxidation by NO<sub>3</sub> radical

Abstract. Highly oxygenated organic molecules (HOM) are found to play an important role in the formation and growth of secondary organic aerosol (SOA). SOA is an important type of aerosol with significant impact on air quality and climate. Compared with the oxidation of volatile organic compounds by ozone (O3) and hydroxyl radical (OH), HOM formation in the oxidation by nitrate radical (NO3), an important oxidant at nighttime and dawn, has received less attention. In this study, HOM formation in the reaction of isoprene with NO3 was investigated in the SAPHIR chamber (Simulation of Atmospheric PHotochemistry In a large Reaction chamber). A large number of HOM, including monomers (C5), dimers (C10), and trimers (C15), both closed-shell compounds and open-shell peroxy radicals (RO2), were identified and were classified into various series according to their formula. Their formation pathways were proposed based on the peroxy radicals observed and known mechanisms in the literature, which were further constrained by the time profiles of HOM after sequential isoprene addition to differentiate first- and second-generation products. HOM monomers containing one to three N atoms (1–3N-monomers) were formed, starting with NO3 addition to carbon double bond, forming peroxy radicals, followed by autoxidation. 1N-monomers were formed by both the direct reaction of NO3 with isoprene and of NO3 with first-generation products. 2N-monomers (e.g., C5H8N2On(n=7–13), C5H10N2On(n=8–14)) were likely the termination products of C5H9N2On⚫, which was formed by the addition of NO3 to C5-hydroxynitrate (C5H9NO4), a first-generation product containing one carbon double bond. 2N-monomers, which were second-generation products, dominated in monomers and accounted for ∼34 % of all HOM, indicating the important role of second-generation oxidation in HOM formation in the isoprene + NO3 reaction under our experimental conditions. H shift of alkoxy radicals to form peroxy radicals and subsequent autoxidation (“alkoxy–peroxy” pathway) was found to be an important pathway of HOM formation. HOM dimers were mostly formed by the accretion reaction of various HOM monomer RO2 and via the termination reactions of dimer RO2 formed by further reaction of closed-shell dimers with NO3 and possibly by the reaction of C5–RO2 with isoprene. HOM trimers were likely formed by the accretion reaction of dimer RO2 with monomer RO2. The concentrations of different HOM showed distinct time profiles during the reaction, which was linked to their formation pathway. HOM concentrations either showed a typical time profile of first-generation products, second-generation products, or a combination of both, indicating multiple formation pathways and/or multiple isomers. Total HOM molar yield was estimated to be 1.2 %-0.7%+1.3%, which corresponded to a SOA yield of ∼3.6 % assuming the molecular weight of C5H9NO6 as the lower limit. This yield suggests that HOM may contribute a significant fraction to SOA yield in the reaction of isoprene with NO3.

Numerical simulation of hydraulic transients in hydropower station with complicated water conveyance system

With the construction of Jinping II hydropower station with extra-large diversion tunnel, Baihetan hydropower station with complicated tailwater system, Yangfanggou hydropower station with complicated surge chamber, High requirements are put forward for the calculation and analysis of hydraulic transient process. In the water transmission and power generation system design of these projects, the following technical problems was solved by Power China Huadong Engineering Corporation Limited: simulation of variable impedance coefficient of differential surge chamber, interface tracking analysis of free-surface-pressure flow system, narrow upper chamber open channel effect and simulation of trigeminy surge chamber. The solutions to those technical problems and were summarized, and the hydraulic transients simulation and calculation software of complex waterway system - Hysim was developed. Based on the field test of Jinping II Hydropower Station, Baihetan hydropower station and Yangfanggou hydropower station, the simulation calculations were carried out. The accuracy of the hydraulic transients simulation software is verified by comparing the inversion results with the test values.

Rapid sulfate formation from synergetic oxidation of SO2 by O3 and NO2 under ammonia-rich conditions: Implications for the explosive growth of atmospheric PM2.5 during haze events in China

Extremely high levels of atmospheric sulfate aerosols have still frequently occurred in China especially in winter haze periods and often been underestimated by models due to some missing formation mechanisms. Here we investigated the heterogeneous reaction dynamics of SO2 oxidation by the abundantly co-existing O3 and NO2 in the urban atmosphere of China by using a laboratory smog chamber simulation technique. Our results showed that with an increase of NH3 concentrations from 0.05 ppm to 1.5 ppm, SO2 oxidation by O3 can be greatly promoted and lead to an exponential increase of diameter growth factor (GF) of particles in the chamber from 1.29 to 1.98 for NaCl seeds and from 1.20 to 1.60 for (NH4)2SO4 seeds, along with an increasing uptake coefficient (γ) of SO2 from 4.47 × 10−5 to 1.52 × 10−4 on NaCl seeds and from 2.32 × 10−5 to 5.74 × 10−5 on (NH4)2SO4 seeds, respectively. The heterogeneous production of sulfate from oxidation of SO2 under NH3-rich conditions by O3 and NO2 mixture in the chamber was 2.0–3.5 times the sum of sulfate from SO2 oxidations by O3 and NO2, suggesting a strongly synergetic effect of the mixed oxidants on the heterogeneous oxidation of SO2, which can cause rapid formation of (NH4)2SO4 and NH4NO3 and is responsible for the explosive growth of PM2.5 in the winter haze period of China. Our chamber results further showed that such synergetic process is only efficient under NH3-rich conditions, clearly indicating that the combined controls on O3, NOx and NH3 are necessary for further mitigating the PM2.5 pollution in China.

Efficient Heterogeneous Formation of Ammonium Nitrate on the Saline Mineral Particle Surface in the Atmosphere of East Asia during Dust Storm Periods.

To understand the chemical evolution of dust in the current East Asian atmosphere, the chemistry of PM2.5 and size-resolved aerosols in Shanghai, China, during the 2019 dust storm event was investigated. Our results showed that concentrations of SO42- in the city during the event highly correlated with Ca2+ and Na+ due to the direct emissions of CaSO4 and Na2SO4 from the upwind deserts. In contrast, during the event, NO3- linearly correlated with NH4+ at a molar ratio close to 1:1, and both almost entirely stayed in coarse particles, suggesting they accumulated on the dust surface as NH4NO3. Based on the field observations and laboratory smog chamber simulations, we found that NO2 and O3 in Shanghai during the dust period reacted to form N2O5, which subsequently hydrolyzed into HNO3 on the surface of saline mineral dusts (e.g., CaSO4 and Na2SO4) and was further neutralized by NH3 as NH4NO3. The relative abundances of NO3- and NH4+ in Shanghai during the dust event were notably higher than those a decade ago, indicating that this heterogeneous formation of NH4NO3 on dust was enhanced by the abundantly coexisting NOx, O3, and NH3 in the current East Asian atmosphere, which should be considered in future modeling studies.

On the mechanisms of flame kernel extinction or survival during aeronautical ignition sequences: Experimental and numerical analysis

A detailed study of ignition mechanisms in aeronautical burners is presented. Even in global lean conditions, liquid fuel injection leads to a strong stratification of the mixture and the whole range of equivalence ratio may be encountered in the burner. The observation of ignition sequences in a representative lab-scale combustion chamber, from both experiment and numerical simulation, reveals a variety of scenarios leading to success or failure. In particular the occurrence of quenching events after the successful creation of a kernel flame is a key mechanism for the outcome of the sequence. A first analysis leads to a classification of ignition/extinction scenarios, based on similar time evolutions and trajectories of the flame kernel. It is found that this classification is much dependent on the sparking location. Then a deeper analysis allows to decompose all scenarios in a succession of more simple, generic mechanisms which are independent of the geometry. This decomposition is a useful tool to describe, understand and predict ignition sequences, being successful or not, in any combustion chamber geometry.

A fully decoupled iterative method with three-dimensional anisotropic immersed finite elements for Kaufman-type discharge problems

In order to simulate the Kaufman-type discharge problems, a fully decoupled iterative method with anisotropic immersed finite elements on Cartesian meshes is proposed, especially for a three-dimensional (3D) non-axisymmetric anisotropic hybrid model which is more difficult than the axisymmetric or isotropic models. The classical hybrid model, which describes the important plasma distribution of the Kaufman-type discharge problems, couples several difficult equations together to form a large scale system. The 3D non-axisymmetric and anisotropic properties will further increase the complexity of this system. Hence it generally needs to be solved in the decoupled way for significantly reducing the computational cost. Based on the Particle-in-Cell Monte Carlo collision (PIC–MCC) method and the immersed finite element (IFE) method, we propose a fully decoupled iterative method for solving this complex system. The IFE method allows Cartesian meshes for general interface problems, while the traditional finite element methods require body-fitting meshes which are often unstructured. Compared with traditional finite element methods, this feature significantly improves the efficiency of the proposed 3D fully decoupled iterative method, while maintaining the optimal accuracy of the chosen finite elements. Numerical simulations of traditional Kaufman ion thruster and annular ion thruster discharge chambers are provided and compared with the corresponding lab experiment results to illustrate the features of the proposed method.

New particle formation (NPF) events in China urban clusters given by sever composite pollution background

New Particle Formation (NPF) refers to transformation of gaseous precursors in the atmosphere due to nucleation and subsequent growth process through physicochemical interaction. It has generated a lot of interest due to its profound impact on global and regional environment, climate and human health. We reviewed the studies on NPF in three city clusters of China: the North China Plain, the Yangtze River Delta and the Pearl River Delta obtained through experiment simulations (e.g., chamber simulation, flow-tube simulation, etc.), field observations, and numerical simulations. Due to its atmospheric background pollution and strong oxidation capacities resulting in high source rate of precursors, China’s atmosphere possesses challenges different from those evaluated in previous studies on cleaning sites and other developing countries. Hence, NPF events can simultaneously exhibit high condensable sink, formation rate and growth rate. In addition, the high intensity of anthropogenic emissions in urban China has led to greater diversity of pollutant species involved in NPF nucleation and subsequent growth, compared to the dominant role of biogenic precursors at cleaning sites. Differences in geographical location and industrial structure also lead to significant distinctions in NPF characteristics of the three city clusters. Consequently, the lack of understanding of nucleation mechanism of complexly polluted background sites makes the global and regional climate models with submodels based on clean background have enormous uncertainty when applied to urban China. The establishment of a mature research ecosystem including field observations, laboratory simulations and numerical simulations is the key to the breakthrough of NPF research in China.

English

The high shear rate thrombus formation was only recently recognized as another way of thrombosis. Models proposed in Weller (2008), (2010) take into account this type of thrombosis. This work uses the phase-field method to model these evolving interface problems. A loosely coupled iterative procedure is introduced to solve the coupled system of equations. Convergence behavior on two levels of refinement of perfusion chamber geometry and cylinder geometry is then studied. The perfusion chamber simulations show good agreement with the original results of Weller. The code is implemented in FEM-library deal.ii Alzeta et al. (2018), which enables distribution of computations to large number of processing units. A scalability and numerical performance study of the loosely coupled iterative procedure is performed, combined with several preconditioners for the linear subproblems.

Experimental measurement of ionization chamber angular response and associated magnetic field correction factors in MR-linac.

To measure ionization chamber dose response as a function of the angle between magnetic field direction and ionization chamber orientation in magnetic resonance-guided radiation therapy (MRgRT) system, and to evaluate angular dependence of magnetic field correction factor for reference dosimetry. Measurements were performed on an Elekta MR-linac that integrates a 1.5-T Philips MRI and a 7-MV FFF photon beam accelerator. The response of four reference class chambers (Exradin-A19, A1SL, IBA FC65-G, and CC13, paired with a PTW UE electrometer) was studied. An in-house built MR-compatible water tank and an accompanying cylindrical insert that allowed chamber rotation around the cylinder's axis was used. The EPID onboard imaging was used to center chamber at the MR-linac isocenter (143.5cm, SAD), as well as to verify position at each datapoint. A clear angular dependence of dose response for all chambers has been measured. The most significant effect of magnetic field on relative chamber response in the presence of magnetic field was observed in the orientation when chamber axis is perpendicular to the direction of magnetic field with the tip pointing in the same direction as Lorentz force. This effect is more pronounced for larger volume chambers; the maximum relative variation in the chamber response (between the setup described above and the one where chamber and magnetic field are parallel) is a 5.3% and 4.6% increase for A19 and FC65-G, respectively, and only 2.0% and 1.9% for smaller volume A1SL and CC13 chamber, respectively. We measured the absolute magnitude of the magnetic field correction factor for the Exradin-A19, A1SL, IBA FC65-G, and CC13 to be 0.938±1.13%, 0.968±0.99%, 0.950±1.13%, and 0.975±1.13%, respectively. The values are for perpendicular orientation of the chamber relative to magnetic field and parallel to the Lorentz force. Experimental measurements carried out in this study have verified the optimal orientation of ionization chamber in terms of minimizing effect of magnetic field on the chamber dose response. This study provides a detailed high-resolution measurement of absolute values for four reference class chambers as a function of the angle between ionization chamber's central axis and the direction of strong magnetic field over a full 360° rotation. The experimental results of this study can further be used for optimization of the actual sensitive volume of the chamber (and analysis of dead volume) in future Monte Carlo chamber simulations in the presence of strong magnetic fields. In addition, it will provide some necessary data for future reference dosimetry protocols for MR-linac.

COSICAF, a fission chamber simulation tool for academic purposes

Since a few years, simulation codes were built at CEA Cadarache to predict the signal of ionisation chambers and taylor detectors for specific applications. It is proposed here to present COSICAF, a tool developed for mainly academic purpose and rapid fission chamber prototyping. This numerical simulation, mostly based on semi-empirical models and Monte-Carlo method will help students to understand how ionisation chambers work. Through the paper, models and their numerical implementation will be discussed. A focus is made on recently implemented features like charge multiplication and correlated source which make the simulation of proportional counter possible. To demonstrate the interest of the code, simulations of a planar fission chamber is proposed.

COSICAF, a fission chamber simulation tool for academic purposes

Nuclear instrumentation is a complex topic since it involves a wide range of physics phenomena like nuclear reactions, heavy ion interactions with matter, electrostatic, charge creation etc. Understanding and modelling fission chambers is a difficult task usually performed with Monte-Carlo and finite element simulations. Since a few years, analytical and simplified Monte Carlo models were introduced at the French Atomic Energy Commission to easily design detectors. It is proposed here to present the derivation of such model, called COSICAF, for academic purposes; this numerical model provided with this article, will help students and researchers to understand and design fission chambers. To demonstrate the interest and the limitation of proposed work in research field, the model is applied to simulate two real miniature fission chamber designs.

Simulations of a parallel plate Fission Chamber using the MCNPX simulation code

MCNPX in its latest version is able to simulate the transportation of Fission Fragments. It opens the door to Fission Chamber simulations. Such simulations are not straightforward and comparisons with experimental spectra often failed. A procedure is described in the present paper to perform such simulation and to process the result to obtain realistic spectra. Simulated spectra are compared with experimental ones in various conditions to validate the method and to present its limitations.

Aerosol mass yields of selected biogenic volatile organic compounds – a theoretical study with nearly explicit gas-phase chemistry

Abstract. In this study we modeled secondary organic aerosol (SOA) mass loadings from the oxidation (by O3, OH and NO3) of five representative biogenic volatile organic compounds (BVOCs): isoprene, endocyclic bond-containing monoterpenes (α-pinene and limonene), exocyclic double-bond compound (β-pinene) and a sesquiterpene (β-caryophyllene). The simulations were designed to replicate an idealized smog chamber and oxidative flow reactors (OFRs). The Master Chemical Mechanism (MCM) together with the peroxy radical autoxidation mechanism (PRAM) were used to simulate the gas-phase chemistry. The aim of this study was to compare the potency of MCM and MCM + PRAM in predicting SOA formation. SOA yields were in good agreement with experimental values for chamber simulations when MCM + PRAM was applied, while a stand-alone MCM underpredicted the SOA yields. Compared to experimental yields, the OFR simulations using MCM + PRAM yields were in good agreement for BVOCs oxidized by both O3 and OH. On the other hand, a stand-alone MCM underpredicted the SOA mass yields. SOA yields increased with decreasing temperatures and NO concentrations and vice versa. This highlights the limitations posed when using fixed SOA yields in a majority of global and regional models. Few compounds that play a crucial role (>95 % of mass load) in contributing to SOA mass increase (using MCM + PRAM) are identified. The results further emphasized that incorporating PRAM in conjunction with MCM does improve SOA mass yield estimation.

An Efficiency Studying of an Ion Chamber Simulation Using Vriance Reduction Techniques with EGSnrc.

Background:Radiotherapy is an important technique of cancer treatment using ionizing radiation. The determination of total dose in reference conditions is an important contribution to uncertainty that could achieve 2%. The source of this uncertainty comes from cavity theory that relates the in-air cavity dose and the dose to water. These correction factors are determined from Monte Carlo calculations of ionization chambers. The main problem of this type of calculation is the extremely long computation time to achieve reasonable statistics. Objective:The main purpose of this work is to present a combination with variance reduction techniques for the case of an ionization chamber in water. Methods:The egs_chamber code allows for very efficient computation of ionization chamber doses and dose ratios by using various variance reduction techniques, and also permits realistic simulations of the experimental setup due to the use of EGSnrc C++ library. Russian roulette and Photon Cross Section Enhancement were used with egs_chamber code. Tests were performed to obtain the parameters of variance reduction techniques resulting in a maximum efficiency. Results:It can be seen that the parameters which result in improved Monte Carlo calculation of the efficiency values are XCSE 64 and Russian Roulette (RR) 128. Conclusion:This study determines the parameters of variance reduction techniques that result in an optimal computational efficiency.

An Efficiency Studying of an Ion Chamber Simulation Using Vriance Reduction Techniques with EGSnrc

Background: Radiotherapy is an important technique of cancer treatment using ionizing radiation. The determination of total dose in reference conditions is an important contribution to uncertainty that could achieve 2%. The source of this uncertainty comes from cavity theory that relates the in-air cavity dose and the dose to water. These correction factors are determined from Monte Carlo calculations of ionization chambers. The main problem of this type of calculation is the extremely long computation time to achieve reasonable statistics.Objective: The main purpose of this work is to present a combination with variance reduction techniques for the case of an ionization chamber in water.Methods: The egs_chamber code allows for very efficient computation of ionization chamber doses and dose ratios by using various variance reduction techniques, and also permits realistic simulations of the experimental setup due to the use of EGSnrc C++ library. Russian roulette and Photon Cross Section Enhancement were used with egs_chamber code. Tests were performed to obtain the parameters of variance reduction techniques resulting in a maximum efficiency.Results: It can be seen that the parameters which result in improved Monte Carlo calculation of the efficiency values are XCSE 64 and Russian Roulette (RR) 128.Conclusion: This study determines the parameters of variance reduction techniques that result in an optimal computational efficiency.

Penumbra Measurements and Comparison of In-House and Standard Circular Cones by the Gafchoromic Film, Pinpoint Ion Chamber, and MCNPX Monte Carlo Simulation

Introduction: Penumbra is an important property of the radiation beam to obtain a suitable margin surrounding the target volume. Therefore, the precise penumbra width determination in stereotactic radiotherapy is necessary for treatment planning. This study aimed to compare the obtained results of penumbra width by in-house and standard circular cones by different dosimeters, as well as evaluating the function of EBT3 for dosimetric properties of the small field radiation. Material and Methods: Different circular cones were mounted on the head of the accelerator to produce 12, 20, and 40 mm field sizes at isocenter. Dosimetric measurements were performed with the EBT3 film, PinPoint ion chamber. Afterwards, MCNPX Monte Carlo simulation was used to evaluate the dosimetric parameters. Results: According to the obtained results, the penumbra width was increased by larger diameters of circular cones. The obtained measured data by PinPoint ion chamber showed a larger penumbra width compared to those calculated by Monte Carlo at all field sizes. The gamma index analysis revealed distance-to-agreement and dose-difference of 2 mm /2%/ at all points. The results of this study showed that source to diaphragm distance had a major role in penumbra size determination of small field dosimetry with PinPoint ion chamber, EBT3 film, and Monte Carlo simulation. Conclusion: As findings of this study reported, EBT3 films are reliable detectors for relative dosimetry due to high spatial resolution for small field sizes. Furthermore, they can be used for measuring beam profile and percentage depth dose curves.

Field Test and Analysis of Energy-Saving Effects of Energy-Recovery Ventilators on Heat-Pump Electricity Consumption in a Classroom

Energy-recovery ventilators (ERVs) are regarded as important energy-saving systems in buildings. It has been reported that they have high energy-saving rates compared with conventional ventilators that operate without energy recovery, but the saving rates have been obtained typically by employing chamber tests and simulations. In this work, a field-test method is proposed that uses a single test room but alternates the tested ventilation modes hourly. This proposed method is useful because an additional comparison room is not always available and can be a source of uncertainty for field tests. The test is performed in a classroom during a heating period, and the results are calibrated to account for different experimental conditions during the test period. The calibrated energy-saving rates indicate the effectiveness of the ERV; however, they are lower in the early hours of the system operation, for two reasons: (1) the maximum power control schemes of the heat pumps are applied for cases where the indoor temperatures are far lower than the set-point temperature; (2) the ventilation load seemingly represents a decreasing proportion of the total heating load in early hours owing to the thermal-capacity effects for the building, which was cooled for many hours. The findings are verified via a chamber test and simulations. As a consequence, it is important to account for actual system characteristics affected by the thermal behaviors of classrooms when the overall performance of a system is evaluated.