Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Systematic error”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,531 records · Page 85

A COMPARISON OF STELLAR ELEMENTAL ABUNDANCE TECHNIQUES AND MEASUREMENTS

Stellar elemental abundances are important for understanding the fundamental properties of a star or stellar group, such as age and evolutionary history, as well as the composition of an orbiting planet. However, as abundance measurement techniques have progressed, there has been little standardization between individual methods and their comparisons. As a result, different stellar abundance procedures determine measurements that vary beyond the quoted error for the same elements within the same stars. The purpose of this paper is to better understand the systematic variations between methods and offer recommendations for producing more accurate results in the future. We invited a number of participants from around the world (Australia, Portugal, Sweden, Switzerland, and the United States) to calculate 10 element abundances (C, O, Na, Mg, Al, Si, Fe, Ni, Ba, and Eu) using the same stellar spectra for four stars (HD 361, HD 10700, HD 121504, and HD 202206). Each group produced measurements for each star using (1) their own autonomous techniques, (2) standardized stellar parameters, (3) a standardized line list, and (4) both standardized parameters and a line list. We present the resulting stellar parameters, absolute abundances, and a metric of data similarity that quantifies the homogeneity of the data. We conclude that standardization of some kind, particularly stellar parameters, improves the consistency between methods. However, because results did not converge as more free parameters were standardized, it is clear there are inherent issues within the techniques that need to be reconciled. Therefore, we encourage more conversation and transparency within the community such that stellar abundance determinations can be reproducible as well as accurate and precise.

79 ASTRONOMY AND ASTROPHYSICS↗

Optimizing Long-term Monitoring of Radiation Air Dose Rates - 19169

Radiation air dose rates near the Fukushima Daiichi Nuclear Power Plant (NPP) have been steadily decreasing over the past seven years since the accident in March 2011. After seven years, the radiation monitoring program is expected to transition from short-term emergency responses to long-term monitoring. The objectives of long-term monitoring are often different from early emergency responses, since extensive data accumulation has led to an understanding of contaminant distributions and mobility. The main long-term monitoring objectives are (1) to confirm system stability and the continuing reduction of contaminant and hazard levels, (2) to provide assurance for the public, (3) to accumulate the basic datasets for scientific knowledge and future preparation, and (4) to detect changes or anomalies in contaminant mobility (if they occur), or any unexpected processes or events. In this work, we have developed a methodology for optimizing the sampling locations of radiation air dose-rate monitoring. Three steps are taken in order to determine sampling locations in a systematic manner: (1) prioritizing the critical locations, such as schools or regulatory requirement locations, (2) diversifying locations across the key environmental controls that are known to influence contaminant mobility and distributions, and (3) capturing the heterogeneity of radiation air dose rates across the domain. In the second step, we use a Gaussian mixture model to identify the representative locations in multiple environmental variables, such as elevation and land-cover types. In the third step, we use a Gaussian process model to capture and estimate the heterogeneity of air dose rates across the domain. We use the integrated dose-rate map from the Bayesian geostatistical method as the reference map and distribute the sampling in such a way as to capture the heterogeneity of the reference map. Although we assume hypothetical sets of measurement locations to evaluate the performance of the methodology and algorithms, we have demonstrated this approach using real, spatially extensive datasets, including the integrated air-dose-rate map of 2016. Our results have shown that this approach allows us to add sampling locations in a systematic manner. Increasing the number of sampling locations can better capture the heterogeneity of dose rates, although the estimation error does not decrease further after a certain number of samples. We have also found that when there are restrictions such as pre-existing monitoring locations or the ones along roads, the spatial estimation becomes poor even with the same number of monitoring locations. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Two-step fabrication technique of gold tips for use in point-contact spectroscopy

For a successful point-contact spectroscopy (PCS) measurement, metallic tips of proper shape and smoothness are essential to ensure the ballistic nature of a point-contact junction. Until recently, the fabrication of Au tips suitable for use in point-contact spectroscopy has remained more of an art involving a trial and error method rather than an automated scientific process. To address these issues, we have developed a technique with which one can prepare high quality Au tips reproducibly and systematically. It involves an electronic control of the driving voltages used for an electrochemical etching of a gold wire in a HCl-glycerol mixture or a HCl solution. We find that a stopping current, below which the circuit is set to shut off, is a single very important parameter to produce an Au tip of desired shape. We present detailed descriptions for a two-step etching process for Au tips and also test results from PCS measurements using them.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

BeyondPlanck: I. Global Bayesian analysis of the Planck Low Frequency Instrument data

We describe the BEYONDPLANCK project in terms of our motivation, methodology, and main products, and provide a guide to a set of companion papers that describe each result in more detail. Building directly on experience from ESA’s Planck mission, we implemented a complete end-to-end Bayesian analysis framework for the Planck Low Frequency Instrument (LFI) observations. The primary product is a full joint posterior distribution P(ω | d), where ω represents the set of all free instrumental (gain, correlated noise, bandpass, etc.), astrophysical (synchrotron, free-free, thermal dust emission, etc.), and cosmological (cosmic microwave background – CMB – map, power spectrum, etc.) parameters. Some notable advantages of this approach compared to a traditional pipeline procedure are seamless end-to-end propagation of uncertainties; accurate modeling of both astrophysical and instrumental effects in the most natural basis for each uncertain quantity; optimized computational costs with little or no need for intermediate human interaction between various analysis steps; and a complete overview of the entire analysis process within one single framework. As a practical demonstration of this framework, we focus in particular on low-ℓ CMB polarization reconstruction with Planck LFI. In this process, we identify several important new effects that have not been accounted for in previous pipelines, including gain over-smoothing and time-variable and non-1/f correlated noise in the 30 and 44 GHz channels. Modeling and mitigating both previously known and newly discovered systematic effects, we find that all results are consistent with the ΛCDM model, and we constrained the reionization optical depth to τ = 0.066 ± 0.013, with a low-resolution CMB-based χ 2 probability to exceed of 32%. This uncertainty is about 30% larger than the official pipelines, arising from taking a more complete instrumental model into account. The marginal CMB solar dipole amplitude is 3362.7 ± 1.4 μK, where the error bar was derived directly from the posterior distribution without the need of any ad hoc instrumental corrections. We are currently not aware of any significant unmodeled systematic effects remaining in the Planck LFI data, and, for the first time, the 44 GHz channel is fully exploited in the current analysis. We argue that this framework can play a central role in the analysis of many current and future high-sensitivity CMB experiments, including LiteBIRD, and it will serve as the computational foundation of the emerging community-wide COSMOGLOBE effort, which aims to combine state-of-the-art radio, microwave, and submillimeter data sets into one global astrophysical model.

79 ASTRONOMY AND ASTROPHYSICS↗

Enhanced reaction rates in NDP analysis with neutron scattering

Neutron depth profiling (NDP) makes accessible quantitative information on a few isotopic concentration profiles ranging from the surface into the sample a few micrometers. Because the candidate analytes for NDP are few, there is little interference encountered. Furthermore, neutrons have no charge so mixed chemical states in the sample are of no direct concern. There are a few nuclides that exhibit large probabilities for neutron scattering. The effect of neutron scattering on NDP measurements has not previously been evaluated as a basis for either enhancing the reaction rates or as a source of measurement error. Hydrogen is a common element exhibiting large neutron scattering probability found in or around sample volumes being analyzed by NDP. A systematic study was conducted to determine the degree of signal change when neutron scattering occurs during analysis. The relative signal perturbation was evaluated for materials of varied neutron scattering probability, concentration, total mass, and geometry. Signal enhancements up to 50% are observed when the hydrogen density is high and in close proximity to the region of analysis with neutron beams of sub thermal energies. Greater signal enhancements for the same neutron number density are reported for thermal neutron beams. Even adhesive tape used to position the sample produces a measureable signal enhancement. Because of the shallow volume, negligible distortion of the NDP measured profile shape is encountered from neutron scattering.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Rapid systematic assessment of the detection and attribution of regional anthropogenic climate change

Despite being a well-established research field, the detection and attribution of observed climate change to anthropogenic forcing is not yet provided as a climate service. One reason for this is the lack of a methodology for performing tailored detection and attribution assessments on a rapid time scale. Here we develop such an approach, based on the translation of quantitative analysis into the “confidence” language employed in recent Assessment Reports of the Intergovernmental Panel on Climate Change. While its systematic nature necessarily ignores some nuances examined in detailed expert assessments, the approach nevertheless goes beyond most detection and attribution studies in considering contributors to building confidence such as errors in observational data products arising from sparse monitoring networks. When compared against recent expert assessments, the results of this approach closely match those of the existing assessments. Where there are small discrepancies, these variously reflect ambiguities in the details of what is being assessed, reveal nuances or limitations of the expert assessments, or indicate limitations of the accuracy of the sort of systematic approach employed here. Deployment of the method on 116 regional assessments of recent temperature and precipitation changes indicates that existing rules of thumb concerning the detectability of climate change ignore the full range of sources of uncertainty, most particularly the importance of adequate observational monitoring.

54 ENVIRONMENTAL SCIENCES↗

High-throughput cancer hypothesis testing with an integrated PhysiCell-EMEWS workflow

BackgroundCancer is a complex, multiscale dynamical system, with interactions between tumor cells and non-cancerous host systems. Therapies act on this combined cancer-host system, sometimes with unexpected results. Systematic investigation of mechanistic computational models can augment traditional laboratory and clinical studies, helping identify the factors driving a treatment's success or failure. However, given the uncertainties regarding the underlying biology, these multiscale computational models can take many potential forms, in addition to encompassing high-dimensional parameter spaces. Therefore, the exploration of these models is computationally challenging. We propose that integrating two existing technologiesone to aid the construction of multiscale agent-based models, the other developed to enhance model exploration and optimizationcan provide a computational means for high-throughput hypothesis testing, and eventually, optimization.ResultsIn this paper, we introduce a high throughput computing (HTC) framework that integrates a mechanistic 3-D multicellular simulator (PhysiCell) with an extreme-scale model exploration platform (EMEWS) to investigate high-dimensional parameter spaces. We show early results in applying PhysiCell-EMEWS to 3-D cancer immunotherapy and show insights on therapeutic failure. We describe a generalized PhysiCell-EMEWS workflow for high-throughput cancer hypothesis testing, where hundreds or thousands of mechanistic simulations are compared against data-driven error metrics to perform hypothesis optimization.ConclusionsWhile key notational and computational challenges remain, mechanistic agent-based models and high-throughput model exploration environments can be combined to systematically and rapidly explore key problems in cancer. These high-throughput computational experiments can improve our understanding of the underlying biology, drive future experiments, and ultimately inform clinical practice.

Agent-based model↗

Uncertainty of Rotating Shadowband Irradiometers and Si-Pyranometers Including the Spectral Irradiance Error

Concentrating solar power projects require accurate direct normal irradiance (DNI) data including uncertainty specifications for plant layout and cost calculations. Ground measured data are necessary to obtain the required level of accuracy and are often obtained with Rotating Shadowband Irradiometers (RSI) that use photodiode pyranometers and correction functions to account for systematic effects. The uncertainty of Si-pyranometers has been investigated, but so far basically empirical studies were published or decisive uncertainty influences had to be estimated based on experience in analytical studies. One of the most crucial estimated influences is the spectral irradiance error because Si-photodiode-pyranometers only detect visible and color infrared radiation and have a spectral response that varies strongly within this wavelength interval. Furthermore, analytic studies did not discuss the role of correction functions and the uncertainty introduced by imperfect shading. In order to further improve the bankability of RSI and Si-pyranometer data, a detailed uncertainty analysis following the Guide to the Expression of Uncertainty in Measurement (GUM) has been carried out. The study defines a method for the derivation of the spectral error and spectral uncertainties and presents quantitative values of the spectral and overall uncertainties. Data from the PSA station in southern Spain was selected for the analysis. Average standard uncertainties for corrected 10 min data of 2% for global horizontal irradiance (GHI), and 2.9% for DNI (for GHI and DNI over 300 W/m2) were found for the 2012 yearly dataset when separate GHI and DHI calibration constants were used. Also the uncertainty in 1 min resolution was analyzed. The effect of correction functions is significant. The uncertainties found in this study are consistent with results of previous empirical studies.

14 SOLAR ENERGY↗

Performance Analysis of OpenSHMEM Applications with TAU Commander

The TAU Performance System® (TAU) is a powerful and highly versatile profiling and tracing tool ecosystem for performance engineering of parallel programs. Developed over the last twenty years, TAU has evolved with each new generation of HPC systems and scales efficiently to hundreds of thousands of cores. TAU’s organic growth has resulted in a loosely coupled software toolbox such that novice users first encountering TAU’s complexity and vast array of features are often intimidated and easily frustrated. To lower the barrier to entry for novice TAU users, ParaTools and the US Department of Energy have developed “TAU Commander,” a performance engineering workflow manager that facilitates a systematic approach to performance engineering, guides users through common profiling and tracing workflows, and offers constructive feedback in case of error. This work compares TAU and TAU Commander workflows for common performance engineering tasks in OpenSHMEM applications and demonstrates workflows targeting two different SHMEM implementations, Intel Xeon “Haswell” and “Knights Landing” processors, direct and indirect measurement methods, callsite, profiles, and traces.

97 MATHEMATICS AND COMPUTING↗

Emergent Constraints for Cloud Feedbacks and Climate Sensitivity

Emergent constraints are physically explainable empirical relationships between characteristics of the current climate and long-term climate prediction that emerge in collections of climate model simulations. With the prospect of constraining long-term climate prediction, scientists have recently uncovered several emergent constraints related to long-term cloud feedbacks. We review these proposed emergent constraints, many of which involve the behavior of low-level clouds, and discuss criteria to assess their credibility. With further research, some of the cases we review may eventually become confirmed emergent constraints, provided they are accompanied by credible physical explanations. Because confirmed emergent constraints identify a source of model error that projects onto climate predictions, they deserve extra attention from those developing climate models and climate observations. While a systematic bias cannot be ruled out, it is noteworthy that the promising emergent constraints suggest larger cloud feedback and hence climate sensitivity.

58 GEOSCIENCES↗

Equilibrium clumped-isotope effects in doubly substituted isotopologues of ethane

We combine path-integral Monte Carlo methods with a new intramolecular potential energy surface to quantify the equilibrium enrichment of doubly substituted ethane isotopologues due to clumped-isotope effects. Ethane represents the simplest molecule to simultaneously exhibit 13C–13C, 13C–D, and D–D clumped-isotope effects, and the analysis of corresponding signatures may provide useful geochemical and biogeochemical proxies of formation temperatures or reaction pathways. Utilizing path-integral statistical mechanics, we predict equilibrium fractionation factors that fully incorporate nuclear quantum effects, such as anharmonicity and rotational-vibrational coupling which are typically neglected by the widely used Urey model. The magnitude of the calculated fractionation factors for the doubly substituted ethane isotopologues indicates that isotopic clumping can be observed if rare-isotope substitutions are separated by up to three chemical bonds, but the diminishing strength of these effects suggests that enrichment at further separations will be negligible. The Urey model systematically underestimates enrichment due to 13C–D and D–D clumped-isotope effects in ethane, leading to small relative errors in the apparent equilibrium temperature, ranging from 5 K at 273.15 K to 30 K at 873.15 K. We additionally note that the rotameric dependence of isotopologue enrichment must be carefully considered when using the Urey model, whereas the path-integral calculations automatically account for such effects due to configurational sampling. These findings are of direct relevance to future clumped-isotope studies of ethane, as well as studies of 13C–13C, 13C–D, and D–D clumped-isotope effects in other hydrocarbons.

Geochemistry & Geophysics↗

A spectroscopic study of the extreme black widow PSR J1311–3430

We report on a series of spectroscopic observations of PSR J1311–3430, an extreme black-widow gamma-ray pulsar with a helium-star companion. In a previous study we estimated the neutron star mass as M{sub NS}=2.68±0.14 M{sub ⊙} (statistical error), based on limited spectroscopy and a basic (direct heating) light-curve model; however, much larger model-dependent systematics dominate the mass uncertainty. Our new spectroscopy reveals a range of complex source behavior. The variable He i companion wind emission lines can dominate broadband photometry, especially in red filters or near minimum brightness, and the wind flux should complete companion evaporation in a spin-down time. The heated companion face also undergoes dramatic flares, reaching ∼40,000 K over ∼20% of the star; this is likely powered by a magnetic field generated in the companion. The companion center-of-light radial velocity is now well measured with K{sub CoL}=615.4±5.1 km s{sup −1}. We detect non-sinusoidal velocity components due to the heated face flux distribution. Using our spectra to excise flares and wind lines, we generate substantially improved light curves for companion continuum fitting. We show that the inferred inclination and neutron star mass, however, remain sensitive to the poorly constrained heating pattern. The neutron star’s mass, M{sub NS}, is likely less than the direct heating value and could range as low as 1.8 M{sub ⊙} for extreme equatorial heating concentration. While we cannot yet pin down M{sub NS}, our data imply that an intrabinary shock reprocesses the pulsar emission and heats the companion. Improved spectra and, especially, models that include such shock heating are needed for precise parameter measurement.

79 ASTRONOMY AND ASTROPHYSICS↗

TH-AB-207A-11: Anatomical and Functional Assessment of Coronary Artery Disease Using Low-Dose Whole-Organ Dynamic Computed Tomography

Purpose: To evaluate the accuracy of a low-dose, first-pass-analysis (FPA) dynamic computed tomography angiography and perfusion (CTAP) technique, for whole-organ anatomical and functional assessment of coronary artery disease (CAD). Methods: An angioplasty balloon was advanced into the left anterior descending (LAD) coronary artery of five swine (35–45 kg) to induce several levels of stenosis at maximal hyperemia (intracoronary adenosine, 240 µg/min). Reference fluorescence microspheres and intravenous contrast (370 mg/mL iodine, 25 mL, 7 mL/s) were injected centrally and dynamic imaging was performed using a 320-slice CT scanner at 100 kVp and 200 mA. Twenty volume scans were acquired per stenosis level to capture complete aortic and myocardial enhancement curves, but only two volume scans were used for whole-organ dynamic FPA CTAP measurement. All CTAP measurements in the LAD were compared to the reference microsphere perfusion measurements using linear regression, concordance correlation, and Bland-Altman analysis. Results: The result of dynamic FPA CTAP measurement in the LAD was in good agreement with the reference microsphere perfusion measurement (P-CTAP = 1.01 P-MICRO + 0.16, R{sup 2} = 0.95). The root mean square error (RMSE) and difference (RMSD) of measurement were 0.51 mL/min/g and 0.47 mL/min/g, respectively. Bland-Altman analysis demonstrated negligible systematic measurement bias. Additionally, the concordance correlation coefficient (CCC) was found to be ρ = 0.97, indicating excellent agreement between dynamic FPA CTAP measurement and the reference microsphere perfusion measurement. Lastly, the effective dose of the proposed technique using the “simulated” two-volume scan CTAP acquisition protocol was 2.6 mSv; much lower than the ∼10 mSv effective dose of current dynamic CTP techniques alone. Conclusion: The results indicate the potential for significant improvements in CAD assessment through low-dose, quantitative dynamic FPA CTAP. Such improvements are afforded through whole-organ CT scanning technology, and have the potential to improve patient outcomes and lead to healthier patient lives. Conflict of Interest (only if applicable): Grant funding from Toshiba America Medical Systems.

60 APPLIED LIFE SCIENCES↗

AI to Predict Glass Compositions Satisfying Property and Cooling Rate Criteria

This project aimed to develop a predictive, artificial intelligence/machine learning-based model to identify glass compositions satisfying specified property requirements. Such a model would provide a systematic approach for narrowing down the nearly infinite range of possible compositions for glasses and minimize unnecessary experimental trial and error. A large empirical data set for training and testing the algorithm was obtained from the SciGlass database. It contains glass compositions and corresponding property data from a wide range of literature sources. However, the currently available form of this data, recently released under an open database license, is not conducive to easy querying and use. The data structure was deciphered and a customized parsing code developed to make this data more usable for the current and future work. Neural network models were developed and trained on viscosity data from the database and demonstrated potential for improving prediction accuracy over a traditional regression model.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Economic optimization of operations for hybrid energy systems under variable markets

We prosed a hybrid energy systems (HES) which is an important element to enable increasing penetration of clean energy. Our paper investigates the operations flexibility of HES, and develops a methodology for operations optimization for maximizing economic value based on predicted renewable generation and market information. A multi-environment computational platform for performing such operations optimization is also developed. In order to compensate for prediction error, a control strategy is accordingly designed to operate a standby energy storage element (ESE) to avoid energy imbalance within HES. The proposed operations optimizer allows systematic control of energy conversion for maximal economic value. Simulation results of two specific HES configurations are included to illustrate the proposed methodology and computational capability. These results demonstrate the economic viability of HES under proposed operations optimizer, suggesting the diversion of energy for alternative energy output while participating in the ancillary service market. Economic advantages of such operations optimizer and associated flexible operations are illustrated by comparing the economic performance of flexible operations against that of constant operations. Sensitivity analysis with respect to market variability and prediction error, are also performed.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Lu-Hf and Sm-Nd geochronological constraints on the influence of subduction metamorphism in controlling the Hf-Nd terrestrial array: Evidence from the world’s orogenic belts

To create a better understanding of the influence of subduction metamorphism on the Hf-Nd terrestrial array, this study examines all of the available Lu-Hf and Sm-Nd garnet isochron ages from the world’s orogenic eclogites and amphi bolites. The garnet isochron y-intercept values provide a record of each sample’s epsilon Hf ( ε Hf) and Nd ( ε Nd) at the age recorded by the isochron, which in the case of most orogenic eclogites and amphibolites is the time of prograde metamorphism during subduction. It is further possible to assess the extent to which each sample’s Lu-Hf and Sm-Nd isotope systematics have remained undisturbed since the date recorded by the isochron by comparison of the bulk-rock isotope systematics with the isochron initial value. This approach provides a means of understanding the Lu-Hf and Sm-Nd systematics of post-subduction metamorphic oceanic crust. On a global scale, the coupled Lu-Hf and Sm-Nd isotope systematics of this crust closely align with the Hf-Nd terrestrial array and have average 176 Lu/ 177 Hf and 147 Sm/ 144 Nd ratios within error of the chondritic uniform reservoir values. This is an important finding because it suggests that a reservoir of deeply subducted oceanic crust, possibly residing at the core-mantle boundary, is unlikely to be significantly fractionated from the Hf-Nd terrestrial array. In contrast to the samples of subduction-metamorphosed oceanic crust, metamorphic rocks with a continental crustal protolith demonstrate garnet isochron initial values and bulk-rock Lu-Hf and Sm-Nd isotope systematics indicating significant decoupling from the Hf-Nd terrestrial array. However, because this material is not readily subducted, this strongly decoupled signature is unlikely to be transferred into deep Earth. Mixing calculations aimed at developing a clearer understanding of what factors control decoupling of a sample’s Lu-Hf and Sm-Nd systematics from the Hf-Nd terrestrial array highlight the complexity of mass-transfer regimes operating in the subduction environment.

54 ENVIRONMENTAL SCIENCES↗

Evaluating Climate Models’ Cloud Feedbacks Against Expert Judgment

Abstract The persistent and growing spread in effective climate sensitivity (ECS) across global climate models necessitates rigorous evaluation of their cloud feedbacks. Here we evaluate several cloud feedback components simulated in 19 climate models against benchmark values determined via an expert synthesis of observational, theoretical, and high‐resolution modeling studies. We find that models with smallest feedback errors relative to these benchmark values generally have moderate total cloud feedbacks (0.4–0.6 W m −2 K −1 ) and ECS (3–4 K). Those with largest errors generally have total cloud feedback and ECS values that are too large or too small. Models tend to achieve large positive total cloud feedbacks by having several cloud feedback components that are systematically biased high rather than by having a single anomalously large component, and vice versa. In general, better simulation of mean‐state cloud properties leads to stronger but not necessarily better cloud feedbacks. The Python code base provided herein could be applied to developmental versions of models to assess cloud feedbacks and cloud errors and place them in the context of other models and of expert judgment in real‐time during model development.

58 GEOSCIENCES↗

Fast and self-learning indoor airflow simulation based on in situ adaptive tabulation

Fast simulation for stratified indoor airflow distributions is desired for various applications, such as design of advanced indoor environments, emergency management, and coupled annual energy simulation for buildings with stratified air distributions. Reduced order models trained by pre-computed computational fluid dynamics results are fast, but their prediction may be inaccurate when applied for conditions outside the training domain. To overcome this limitation, we propose in this paper a fast and self-learning model based on an in situ adaptive tabulation (ISAT) algorithm, which is trained by a fast fluid dynamics (FFD) model as an example. The idea is that the ISAT will retrieve the solutions from an existing data set if the estimated prediction error is within a pre-defined tolerance. Otherwise, the ISAT will execute the FFD simulation, which is accelerated by running in parallel on a graphics processing unit, for a full-scale simulation. This paper systematically investigates the feasibility of the ISAT for indoor airflow simulations by presenting the ISAT-FFD implementation alongside results related to its overall performance. Using a stratified indoor airflow as an example, we evaluated how the training time of ISAT was impacted by four factors (training methods, error tolerances, number of inputs, and number of outputs). Then we demonstrated that a trained ISAT model can predict the key information for inputs both inside and outside the training domain. The ISAT was able to answer query points both inside and close to training domain using retrieve actions within a time less than 0.001 s for each query. In conclusion, we provided suggestions for using the ISAT for building applications.

24 POWER TRANSMISSION AND DISTRIBUTION↗