Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Statistical techniques”

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 451 records · Page 25

Inverse regression for ridge recovery: a data-driven approach for parameter reduction in computer experiments

Parameter reduction can enable otherwise infeasible design and uncertainty studies with modern computational science models that contain several input parameters. In statistical regression, techniques for sufficient dimension reduction (SDR) use data to reduce the predictor dimension of a regression problem. A computational scientist hoping to use SDR for parameter reduction encounters a problem: a computer prediction is best represented by a deterministic function of the inputs, so data comprised of computer simulation queries fail to satisfy the SDR assumptions. To address this problem, here we interpret SDR methods sliced inverse regression (SIR) and sliced average variance estimation (SAVE) as estimating the directions of a ridge function, which is a composition of a low-dimensional linear transformation with a nonlinear function. Within this interpretation, SIR and SAVE estimate matrices of integrals whose column spaces are contained in the ridge directions’ span; we analyze and numerically verify convergence of these column spaces as the number of computer model queries increases. Moreover, we show example functions that are not ridge functions but whose inverse conditional moment matrices are low-rank. Consequently, the computational scientist should beware when using SIR and SAVE for parameter reduction, since SIR and SAVE may mistakenly suggest that truly important directions are unimportant.

97 MATHEMATICS AND COMPUTING↗

A systematic feature extraction and selection framework for data-driven whole-building automated fault detection and diagnostics in commercial buildings

In data-driven automated fault detection and diagnostics (AFDD) modeling for building energy systems, feature engineering is a critical process of extracting information from high-dimensional and noisy sensor measurement and turning it into informative and representative inputs or features for data-driven modeling. However, few studies specifically discuss the feature engineering, especially the interactions between feature extraction and feature selection in whole-building AFDD. We developed a systematic feature extraction and selection framework for whole-building AFDD. In this framework, features are aggressively extracted from raw sensor data using statistical feature extraction techniques with various window sizes and statistics. With many features extracted, a hybrid feature selection algorithm that combines the filter and wrapper method then selects the best feature set. The framework considers diversity in the duration of fault behavior among fault types in whole-building AFDD, thus achieving high model generalization. We implemented our developed framework in a virtual testbed calibrated with measured data from Oak Ridge National Laboratory's Flexible Research Platform designed to mimic the operation of a typical small commercial building. The AFDD model is trained by the simulation data generated from the virtual testbed. The results show that (1) the developed framework improves the generalization of the AFDD model by 10.7% compared with literature-reported feature extraction and selection methods and (2) features with diverse window sizes and statistics are selected, providing insight into physical systems beyond the current understanding of buildings and faults and improving the detection and diagnostics of multiple fault types.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Micromechanical characterization of shales through nanoindentation and energy dispersive x-ray spectrometry

Shales are heterogeneous sedimentary rocks which typically comprise a variable mineralogy (including compacted clay particles sub-micrometer in size), silt grains, and nanometer sized pores collectively arranged with transversely isotropic symmetry. Moreover, a detailed understanding of the micro- and sub-microscale geomechanics of these minerals is required to improve models of shale strength and stiffness properties. In this paper, we propose a linked experimental–computational approach and validate a combination of grid nanoindentation and Scanning Electron Microscopy (SEM) with Energy and Wavelength Dispersive X-ray Spectrometry (EDS/WDS) at the same spatial locations to identify both the nano-mechanical morphology and local mineralogy of these nanocomposites. The experimental parameters of each method are chosen to assess a similar volume of material. By considering three different shales of varying mineralogy and mechanical diversity, we show through the EMMIX statistical iterative technique that the constituent phases, including highly compacted plate- or sheet-like clay particles, carbonates, silicates, and sulfides, have distinct nano-mechanical morphologies and associated indentation moduli and hardness. Nanoindentation-based strength homogenization analysis determines an average clay packing density, friction coefficient, and solid cohesion for each tested shale sample. Comparison of bulk to microscale geomechanical properties, through bulk porosimetry measurements, reveals a close correspondence between bulk and microscale clay packing densities. Determining the mechanical microstructure and material properties is useful for predictive microporomechanical models of the stiffness and strength properties of shale. Furthermore, the experimental and computational approaches presented here also apply to other chemically and mechanically complex materials exhibiting nanogranular, composite behavior.

36 MATERIALS SCIENCE↗

Quantifying CO 2 storage efficiency factors in hydrocarbon reservoirs. A detailed look at CO 2 Enhanced Oil Recovery

Carbon dioxide (CO 2 ) enhanced oil recovery (EOR) will likely be the primary means of geologic CO 2 storage during the early stages of commercial-scale carbon capture and storage (CCS) because of the inherent economic incentives as well as the abundant experience and demonstrated success in the United States, where CO 2 EOR has been employed since 1974. The work presented here estimates CO 2 storage efficiency factors in CO 2 EOR operations using a unique industry database of CO 2 EOR sites and 12 different reservoir simulation models. The simulation models encompass fluvial clastic and shallow shelf carbonate depositional environments for reservoir depths of 1219 and 2438 m (4000 and 8000 feet) and 7.6-, 20-, and 64-m (25-, 66-, and 209-foot)-thick pay zones. A novel statistical modeling technique incorporating the Michaelis–Menten function is used to generate empirical percentile estimates of CO 2 storage efficiency factors. West Texas San Andres dolomite water alternating gas (WAG) CO 2 flood performance data were used to derive P10, P50, and P90 CO 2 storage efficiency factors of 0.76, 1.28, and 1.74 Mscf/STB (stock tank barrel) of original oil in place. Median CO 2 storage efficiency factors from continuous CO 2 injection following conventional waterflood varied from 15% to 61% and 8% to 40% for fluvial clastic and shallow shelf carbonate simulation models, respectively, while those from WAG injection varied from 14% to 42% and 8% to 31%, respectively. Variation in the CO 2 storage efficiency factors was largely attributable to reservoir depth (a surrogate for reservoir pressure and temperature) and lithology (clastic versus carbonate). The results of this work provide practical information that can be used to quantify CO 2 storage resource estimates in oil reservoirs during CO 2 EOR operations (as opposed to storage following depletion) and the uncertainty associated with those estimates.

54 ENVIRONMENTAL SCIENCES↗

Dependence of Achievable Plan Quality on Treatment Technique and Planning Goal Refinement: A Head-and-Neck Intensity Modulated Radiation Therapy Application

Purpose: To develop a practical workflow for retrospectively analyzing target and normal tissue dose–volume endpoints for various intensity modulated radiation therapy (IMRT) delivery techniques; to develop technique-specific planning goals to improve plan consistency and quality when feasible. Methods and Materials: A total of 165 consecutive head-and-neck patients from our patient registry were selected and retrospectively analyzed. All IMRT plans were generated using the same dose–volume guidelines for TomoTherapy (Tomo, Accuray), TrueBeam (TB, Varian) using fixed-field IMRT (TB-IMRT) or RAPIDARC (TB-RAPIDARC), or Siemens Oncor (Siemens-IMRT, Siemens). A MATLAB-based dose–volume extraction and analysis tool was developed to export dosimetric endpoints for each patient. With a fair stratification of patient cohort, the variation of achieved dosimetric endpoints was analyzed among different treatment techniques. Upon identification of statistically significant variations, technique-specific planning goals were derived from dynamically accumulated institutional data. Results: Retrospective analysis showed that although all techniques yielded comparable target coverage, the doses to the critical structures differed. The maximum cord doses were 34.1 ± 2.6, 42.7 ± 2.1, 43.3 ± 2.0, and 45.1 ± 1.6 Gy for Tomo, TB-IMRT, TB-RAPIDARC, and Siemens-IMRT plans, respectively. Analyses of variance showed significant differences for the maximum cord doses but no significant differences for other selected structures among the investigated IMRT delivery techniques. Subsequently, a refined technique-specific dose–volume guideline for maximum cord dose was derived at a confidence level of 95%. The dosimetric plans that failed the refined technique-specific planning goals were reoptimized according to the refined constraints. We observed better cord sparing with minimal variations for the target coverage and other organ at risk sparing for the Tomo cases, and higher parotid doses for C-arm linear accelerator–based IMRT and RAPIDARC plans. Conclusion: Patient registry–based processes allowed easy and systematic dosimetric assessment of treatment plan quality and consistency. Our analysis revealed the dependence of certain dosimetric endpoints on the treatment techniques. Technique-specific refinement of planning goals may lead to improvement in plan consistency and plan quality.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

A general solution strategy of modified power method for higher mode solutions

A general solution strategy of the modified power iteration method for calculating higher eigenmodes has been developed and applied in continuous energy Monte Carlo simulation. The new approach adopts four features: 1) the eigen decomposition of transfer matrix, 2) weight cancellation for higher modes, 3) population control with higher mode weights, and 4) stabilization technique of statistical fluctuations using multi-cycle accumulations. The numerical tests of neutron transport eigenvalue problems successfully demonstrate that the new strategy can significantly accelerate the fission source convergence with stable convergence behavior while obtaining multiple higher eigenmodes at the same time. The advantages of the new strategy can be summarized as 1) the replacement of the cumbersome solution step of high order polynomial equations required by Booth's original method with the simple matrix eigen decomposition, 2) faster fission source convergence in inactive cycles, 3) more stable behaviors in both inactive and active cycles, and 4) smaller variances in active cycles. Advantages 3 and 4 can be attributed to the lower sensitivity of the new strategy to statistical fluctuations due to the multi-cycle accumulations. The application of the modified power method to continuous energy Monte Carlo simulation and the higher eigenmodes up to 4th order are reported for the first time in this paper. -- Graphical abstract: -- Highlights: •Modified power method is applied to continuous energy Monte Carlo simulation. •Transfer matrix is introduced to generalize the modified power method. •All mode based population control is applied to get the higher eigenmodes. •Statistic fluctuation can be greatly reduced using accumulated tally results. •Fission source convergence is accelerated with higher mode solutions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A general CFD framework for fault-resilient simulations based on multi-resolution information fusion

We develop a general CFD framework for multi-resolution simulations to target multiscale problems but also resilience in exascale simulations, where faulty processors may lead to gappy, in space-time, simulated fields. We combine approximation theory and domain decomposition together with statistical learning techniques, e.g. coKriging, to estimate boundary conditions and minimize communications by performing independent parallel runs. To demonstrate this new simulation approach, we consider two benchmark problems. First, we solve the heat equation (a) on a small number of spatial “patches” distributed across the domain, simulated by finite differences at fine resolution and (b) on the entire domain simulated at very low resolution, thus fusing multi-resolution models to obtain the final answer. Second, we simulate the flow in a lid-driven cavity in an analogous fashion, by fusing finite difference solutions obtained with fine and low resolution assuming gappy data sets. We investigate the influence of various parameters for this framework, including the correlation kernel, the size of a buffer employed in estimating boundary conditions, the coarseness of the resolution of auxiliary data, and the communication frequency across different patches in fusing the information at different resolution levels. In addition to its robustness and resilience, the new framework can be employed to generalize previous multiscale approaches involving heterogeneous discretizations or even fundamentally different flow descriptions, e.g. in continuum-atomistic simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Evolution of microstructure and grain boundary character distribution of a tin bronze annealed at different temperatures

Specimens cut from a rolled tin bronze sheet were annealed at 400–800 °C for 1 h and evolution of their microstructures was then characterized in details by electron channeling contrast imaging and electron backscatter diffraction techniques. Particularly, statistics on special boundaries (SBs) with Σ ≤ 29 and network connectivity of random high angle boundaries (HABs) in the annealed specimens were examined to probe optimization potentials of grain boundary character distribution (GBCD) for this material. Results show that the deformed microstructure in the as-received material begins to be recrystallized when the annealing temperature increase to 500 °C and average grain sizes surge with further increasing temperatures. As a result of the recrystallization, a large number of annealing twins (with Σ3 misorientation) are produced, leading to remarkably increased fractions of SBs (f{sub SBs}). Thanks to preexisting dense low angle boundaries, the majority of SBs in the 500 °C specimen with only partial recrystallization are Σ3{sub ic} (incoherent) boundaries, which effectively disrupt connectivity of random HABs network. Although the f{sub SBs} can be further increased (up to 72.5%) in specimens with full recrystallization (at higher temperatures), the Σ3{sub ic} boundaries would be replaced to some extent by Σ3{sub c} (coherent) boundaries which do not contribute directly to optimizing the GBCD. This work should be able to provide clear suggestions on applying the concept of grain boundary engineering to tin bronze alloys. - Highlights: • The rolled tin bronze begins to be recrystallized as temperature increases to 500 °C. • A lot of SBs are produced after recrystallization and the highest f{sub SBs} is 72.5%. • Partially recrystallized specimen has the optimum GBCD due to more Σ3{sub ic} boundaries. • The Σ3{sub ic} boundaries are replaced by Σ3{sub c} boundaries after full recrystallization.

36 MATERIALS SCIENCE↗

Entropy-driven melting point depression in fcc HEAs

High Entropy Alloys (HEAs) are an increasingly dominant alloy design paradigm. The premise of entropic stabilization of single-phase alloys has motivated much of the research on HEAs. Chemical complexity may indeed help stabilize single alloy phases relative to other lower-entropy competing solid phases. Paradoxically, this complexity may de-stabilize these alloys against the liquid phase, potentially limiting the application space of HEAs at elevated temperatures. In this work, we carry out a comprehensive investigation of the phase stability in the fcc CoCrFeMnNiV-Al HEA space using a state of the art CALPHAD database. By using modern visualization techniques and statistical analysis we examine the trade-off between chemical complexity and stability against the liquid state and identify a potentially difficult to overcome barrier for development of high temperature alloys, at least within the conventional fcc HEA space. Here, limited experimental data seem to be consistent with this analysis.

36 MATERIALS SCIENCE↗

Comparing smoothing techniques for extracting MOSFET threshold voltage

Measurement noise acts as a barrier to the accurate calculation of threshold voltage by derivative-based extraction methods. We examined several smoothing techniques and their effects on measurement noise in threshold voltage extracted by the linear extrapolation method and the g m /I d method. We applied these techniques on a set of SiC power MOSFET devices before and after they had undergone accelerated thermal testing. The smoothing methods examined are least-squares polynomial fitting, low-pass filtering using Butterworth and Chebyshev fillers, and differential smoothing. These methods are compared by their shift in threshold voltage compared to unsmoothed data and their standard deviation among the measurements of each device. All methods performed poorly at smoothing when extracting via the g m /I d due to the large amounts of high-magnitude subthreshold noise, mostly tending to extract threshold voltages far below the expected value. Differential smoothing provided results closest to the expected value with an absolute shift of 1.2 V on average. The smoothing methods performed better with the LE method. On a failed device, all smoothing methods except polynomial smoothing preserved the large threshold voltage anomaly for four out of five measurements. These results show that choice of smoothing method, if used must be evaluated based on the device and extraction method.

42 ENGINEERING↗

Data Science Meets Physical Organic Chemistry

At the heart of synthetic chemistry is the holy grail of predictable catalyst design. In particular, researchers involved in reaction development in asymmetric catalysis have pursued a variety of strategies toward this goal. This is driven by both the pragmatic need to achieve high selectivities and the inability to readily identify why a certain catalyst is effective for a given reaction. While empiricism and intuition have dominated the field of asymmetric catalysis since its inception, enantioselectivity offers a mechanistically rich platform to interrogate catalyst-structure response patterns that explain the performance of a particular catalyst or substrate. In the early stages of an asymmetric reaction development campaign, the overarching mechanism of the reaction, catalyst speciation, the turnover limiting step, and many other details are unknown or posited based on related reactions. Considering the unclear details leading to a successful reaction, initial enantioselectivity data are often used to intuitively guide the ultimate direction of optimization. However, if the conditions of the Curtin-Hammett principle are satisfied, then measured enantioselectivity can be directly connected to the ensemble of diastereomeric transition states (TSs) that lead to the enantiomeric products, and the associated free energy difference between competing TSs (ΔΔ G ‡ = - RT ln[( S )/( R )], where ( S ) and ( R ) represent the concentrations of the enantiomeric products). We, and others, speculated that this important piece of information can be leveraged to guide reaction optimization in a quantitative way. Although traditional linear free energy relationships (LFERs), such as Hammett plots, have been used to illuminate important mechanistic features, we sought to develop data science derived tools to expand the power of LFERs in order to describe complex reactions frequently encountered in modern asymmetric catalysis. Specifically, we investigated whether enantioselectivity data from a reaction can be quantitatively connected to the attributes of reaction components, such as catalyst and substrate structural features, to harness data for asymmetric catalyst design. In this context, we developed a workflow to relate computationally derived features of reaction components to enantioselectivity using data science tools. The mathematical representation of molecules can incorporate many aspects of a transformation, such as molecular features from substrate, product, catalyst, and proposed transition states. Statistical models relating these features to reaction outputs can be used for various tasks, such as performance prediction of untested molecules. Perhaps most importantly, statistical models can guide the generation of mechanistic hypotheses that are embedded within complex patterns of reaction responses. Overall, merging traditional physical organic experiments with statistical modeling techniques creates a feedback loop that enables both evaluation of multiple mechanistic hypotheses and future catalyst design. In this Account, we highlight the evolution and application of this approach in the context of a collaborative program based on chiral phosphoric acid catalysts (CPAs) in asymmetric catalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Multi-objective Bayesian active learning for MeV-ultrafast electron diffraction

Ultrafast electron diffraction using MeV energy beams(MeV-UED) has enabled unprecedented scientific opportunities in the study of ultrafast structural dynamics in a variety of gas, liquid and solid state systems. Broad scientific applications usually pose different requirements for electron probe properties. Due to the complex, nonlinear and correlated nature of accelerator systems, electron beam property optimization is a time-taking process and often relies on extensive hand-tuning by experienced human operators. Algorithm based efficient online tuning strategies are highly desired. Here, we demonstrate multi-objective Bayesian active learning for speeding up online beam tuning at the SLAC MeV-UED facility. The multi-objective Bayesian optimization algorithm was used for efficiently searching the parameter space and mapping out the Pareto Fronts which give the trade-offs between key beam properties. Such scheme enables an unprecedented overview of the global behavior of the experimental system and takes a significantly smaller number of measurements compared with traditional methods such as a grid scan. This methodology can be applied in other experimental scenarios that require simultaneously optimizing multiple objectives by explorations in high dimensional, nonlinear and correlated systems.

43 PARTICLE ACCELERATORS↗

J-PLUS: Systematic impact of metallicity on photometric calibration with the stellar locus

Here, we present the photometric calibration of the twelve optical passbands for the Javalambre Photometric Local Universe Survey (J-PLUS) second data release (DR2), comprising 1088 pointings of two square degrees, and study the systematic impact of metallicity on the stellar locus technique. The [Fe/H] metallicity from the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) for 146 184 high-quality calibration stars, defined with signal-to-noise ratio larger than ten in J-PLUS passbands and larger than three in Gaia parallax, was used to compute the metallicity-dependent stellar locus (ZSL). The initial homogenization of J-PLUS photometry, performed with a unique stellar locus, was refined by including the metallicity effect in colors via the ZSL. The variation of the average metallicity along the Milky Way produces a systematic offset in J-PLUS calibration. This effect is well above 1% for the bluer passbands and amounts 0.07, 0.07, 0.05, 0.03, and 0.02 mag in u, J0378, J0395, J0410, and J0430, respectively. We modeled this effect with the Milky Way location of the J-PLUS pointing, also providing an updated calibration for those observations without LAMOST information. The estimated accuracy in the calibration after including the metallicity effect is at 1% for the bluer J-PLUS passbands and below for the rest. Photometric calibration with the stellar locus technique is prone to significant systematic bias in the Milky Way for passbands bluer than λ = 4500 Å. The calibration method for J-PLUS DR2 reaches 1–2% precision and 1% accuracy for 12 optical filters within an area of 2176 square degrees.

79 ASTRONOMY AND ASTROPHYSICS↗

Implementation of Feldman-Cousins Corrections and Oscillation Calculations in the HPC Environment for the NOvA Experiment

Analysis of neutrino oscillation data involves a combination of complex fitting procedures and statistical correction techniques that are used to determine the full three-flavor PMNS parameters and constraint contours. These techniques rely on computationally intensive “multi-universe” stochastic modeling. The process of calculating these contours and corrections can dominate final stages of the data analysis and become a bottleneck for examining the ef-fect of systematic variations on the final results. As part of the DOE SciDAC-4 sponsored research program, we present a new implementation of a neutrino oscillation fitting and framework to carry out calculations of Feldman-Cousins corrections. This implementation is specifically designed and optimized to operate on modern High-Performance Computing (HPC) facilities. We present the performance of the system in calculating allowed regions for the NOvA experiment based on 8.9 x 1020 and 6.9 x 1020 protons-on-target (POT) neutrino and antineutrino datasets, and compare the performance of this new implementation run at the NERSC supercomputing facility with that from methods used previously by the NOvA collaboration running on grid computing facilities.

Sousa, Alexandre↗

Notes on power of normality tests of error terms in regression models

Normality is one of the basic assumptions in applying statistical procedures. For example in linear regression most of the inferential procedures are based on the assumption of normality, i.e. the disturbance vector is assumed to be normally distributed. Failure to assess non-normality of the error terms may lead to incorrect results of usual statistical inference techniques such as t-test or F-test. Thus, error terms should be normally distributed in order to allow us to make exact inferences. As a consequence, normally distributed stochastic errors are necessary in order to make a not misleading inferences which explains a necessity and importance of robust tests of normality. Therefore, the aim of this contribution is to discuss normality testing of error terms in regression models. In this contribution, we introduce the general RT class of robust tests for normality, and present and discuss the trade-off between power and robustness of selected classical and robust normality tests of error terms in regression models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Enhanced configurational sampling with hybrid non-equilibrium molecular dynamics–Monte Carlo propagator

Molecular dynamics (MD) trajectories based on classical equations of motion can be used to sample the configurational space of complex molecular systems. However, brute-force MD often converges slowly due to the ruggedness of the underlying potential energy surface. Several schemes have been proposed to address this problem by effectively smoothing the potential energy surface. However, in order to recover the proper Boltzmann equilibrium probability distribution, these approaches must then rely on statistical reweighting techniques or generate the simulations within a Hamiltonian tempering replica-exchange scheme. The present work puts forth a novel hybrid sampling propagator combining Metropolis-Hastings Monte Carlo (MC) with proposed moves generated by non-equilibrium MD (neMD). This hybrid neMD-MC propagator comprises three elementary elements: (i) an atomic system is dynamically propagated for some period of time using standard equilibrium MD on the correct potential energy surface; (ii) the system is then propagated for a brief period of time during what is referred to as a “boosting phase,” via a time-dependent Hamiltonian that is evolved toward the perturbed potential energy surface and then back to the correct potential energy surface; (iii) the resulting configuration at the end of the neMD trajectory is then accepted or rejected according to a Metropolis criterion before returning to step 1. A symmetric two-end momentum reversal prescription is used at the end of the neMD trajectories to guarantee that the hybrid neMD-MC sampling propagator obeys microscopic detailed balance and rigorously yields the equilibrium Boltzmann distribution. The hybrid neMD-MC sampling propagator is designed and implemented to enhance the sampling by relying on the accelerated MD and solute tempering schemes. It is also combined with the adaptive biased force sampling algorithm to examine. Illustrative tests with specific biomolecular systems indicate that the method can yield a significant speedup.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Accelerating quantum optics experiments with statistical learning

Quantum optics experiments, involving the measurement of low-probability photon events, are known to be extremely time-consuming. We present a methodology for accelerating such experiments using physically motivated ansatzes together with simple statistical learning techniques such as Bayesian maximum a posteriori estimation based on few-shot data. We show that it is possible to reconstruct time-dependent data using a small number of detected photons, allowing for fast estimates in under a minute and providing a one-to-two order of magnitude speed-up in data acquisition time. We test our approach using real experimental data to retrieve the second order intensity correlation function, G (2) ($τ$), as a function of time delay τ between detector counts, for thermal light as well as anti-bunched light emitted by a quantum dot driven by periodic laser pulses. The proposed methodology has a wide range of applicability and has the potential to impact the scientific discovery process across a multitude of domains.

36 MATERIALS SCIENCE↗

Space-filling designs for computer experiments: A review

Improving the quality of a product/process using a computer simulator is a much less expensive option than the real physical testing. However, simulation using computationally intensive computer models can be time consuming and therefore, directly doing the optimization on the computer simulator can be infeasible. Experimental design and statistical modeling techniques can be used for overcoming this problem. This article reviews experimental designs known as space-filling designs that are suitable for computer simulations. In the review, a special emphasis is given for a recently developed space-filling design called maximum projection design. Furthermore, its advantages are illustrated using a simulation conducted for optimizing a milling process.

97 MATHEMATICS AND COMPUTING↗