Investigation of the Effectiveness of Multi-Fidelity Surrogates on Extrapolation

Author(s):  
Chanyoung Park ◽  
Raphael T. Haftka ◽  
Nam H. Kim

Surrogates have been used as an approximate tool to emulate simulation responses based on a handful of response samples. However, for high fidelity simulations, often only a small number samples are affordable, and this increases the risk of extrapolation using surrogates. Frequently, most of the sampling domain is not in the interpolation domain (called coverage), usually defined as the convex hull of these samples. For example, when we build a surrogate with 20 samples in six-dimensional space, the coverage is merely 2% of the sampling domain. Multi-fidelity surrogates (MFS) may mitigate this problem, because they use large number of low fidelity simulations, so that most of the domain is covered with at least some simulations. This paper explores the extrapolation capability of MFS frameworks through examples including algebraic functions. To examine the effects of different MFS frameworks, we consider six MFS frameworks in terms of their functional forms and frameworks for fitting the forms to data. We consider three different functional forms based on different approaches: 1) a model discrepancy function, 2) model calibration, and 3) both. Bayesian MFS frameworks based on the functional forms are considered. We include also their counterparts in simple frameworks, which have the same functional form but can be built with ready-made surrogates. We examined the effect of the high fidelity sample coverage on extrapolation while the number of high fidelity samples remains the same. The root mean square errors (RMSE) of the interpolation and extrapolation domains are calculated to see their effectiveness on the overall RMSE of whole MFS. For the examples considered, we found that the presence of a regression scalar could be important to extrapolation. Bayesian framework is useful to find a good regression scalar, which simplifies the discrepancy function.

2021 ◽  
Vol 4 ◽  
pp. 28-37
Author(s):  
Alexander Nakonechnyi ◽  
◽  
Grigoriy Kudin ◽  
Taras Zinko ◽  
Petr Zinko ◽  
...  

The issues of parameter estimation in linear regression problems with random matrix coefficients were researched. Given that random linear functions are observed from unknown matrices with random errors that have unknown correlation matrices, the problems of guaranteed mean square estimation of linear functions of matrices were investigated. The estimates of the upper and lower guaranteed standard errors of linear estimates of observations of linear functions of matrices were obtained in the case when the sets are found, for which the unknown matrices and correlation matrices of observation errors are known. It was proved that for some partial cases such estimates are accurate. Assuming that the sets are bounded, convex and closed, more accurate two-sided estimates have been gained for guaranteed errors. The conditions when the guaranteed mean squared errors approach zero as the number of observations increases were found. The necessary and sufficient conditions for the unbiasedness of linear estimates of linear functions of matrices were provided. The notion of quasi-optimal estimates for linear functions of matrices was introduced, and it was proved that in the class of unbiased estimates, quasi-optimal estimates exist and are unique. For such estimates, the conditions of convergence to zero of the guaranteed mean-square errors were obtained. Also, for linear estimates of unknown matrices, the concept of quasi-minimax estimates was introduced and it was confirmed that they are unbiased. For special sets, which include an unknown matrix and correlation matrices of observation errors, such estimates were expressed through the solution of linear operator equations in a finite-dimensional space. For quasi-minimax estimates under certain assumptions, the form of the guaranteed mean squared error of the unknown matrix was found. It was shown that such errors are limited by the sum of traces of the known matrices. An example of finding a minimax unbiased linear estimation was given for a special type of random matrices that are included in the observation equation.


Foods ◽  
2021 ◽  
Vol 10 (4) ◽  
pp. 885
Author(s):  
Sergio Ghidini ◽  
Luca Maria Chiesa ◽  
Sara Panseri ◽  
Maria Olga Varrà ◽  
Adriana Ianieri ◽  
...  

The present study was designed to investigate whether near infrared (NIR) spectroscopy with minimal sample processing could be a suitable technique to rapidly measure histamine levels in raw and processed tuna fish. Calibration models based on orthogonal partial least square regression (OPLSR) were built to predict histamine in the range 10–1000 mg kg−1 using the 1000–2500 nm NIR spectra of artificially-contaminated fish. The two models were then validated using a new set of naturally contaminated samples in which histamine content was determined by conventional high-performance liquid chromatography (HPLC) analysis. As for calibration results, coefficient of determination (r2) > 0.98, root mean square of estimation (RMSEE) ≤ 5 mg kg−1 and root mean square of cross-validation (RMSECV) ≤ 6 mg kg−1 were achieved. Both models were optimal also in the validation stage, showing r2 values > 0.97, root mean square errors of prediction (RMSEP) ≤ 10 mg kg−1 and relative range error (RER) ≥ 25, with better results showed by the model for processed fish. The promising results achieved suggest NIR spectroscopy as an implemental analytical solution in fish industries and markets to effectively determine histamine amounts.


2021 ◽  
pp. 1-14
Author(s):  
Noura Hamze ◽  
Lukas Nocker ◽  
Nikolaus Rauch ◽  
Markus Walzthöni ◽  
Matthias Harders ◽  
...  

BACKGROUND: Accurate segmentation of connective soft tissues in medical images is very challenging, hampering the generation of geometric models for bio-mechanical computations. Alternatively, one could predict ligament insertion sites and then approximate the shapes, based on anatomical knowledge and morphological studies. OBJECTIVE: In this work, we describe an integrated framework for automatic modelling of human musculoskeletal ligaments. METHOD: We combine statistical shape modelling with geometric algorithms to automatically identify insertion sites, based on which geometric surface/volume meshes are created. As clinical use case, the framework has been applied to generate models of the forearm interosseous membrane. Ligament insertion sites in the statistical model were defined according to anatomical predictions following a published approach. RESULTS: For evaluation we compared the generated sites, as well as the ligament shapes, to data obtained from a cadaveric study, involving five forearms with 15 ligaments. Our framework permitted the creation of models approximating ligaments’ shapes with good fidelity. However, we found that the statistical model trained with the state-of-the-art prediction of the insertion sites was not always reliable. Average mean square errors as well as Hausdorff distances of the meshes could increase by an order of magnitude, as compared to employing known insertion locations of the cadaveric study. Using those, an average mean square error of 0.59 mm and an average Hausdorff distance of less than 7 mm resulted, for all ligaments. CONCLUSIONS: The presented approach for automatic generation of ligament shapes from insertion points appears to be feasible but the detection of the insertion sites with a SSM is too inaccurate, thus making a patient-specific approach necessary.


Author(s):  
Pavle Šćepanović ◽  
Frederik A. Döring

AbstractFor a broad range of applications, flight mechanics simulator models have to accurately predict the aircraft dynamics. However, the development and improvement of such models is a difficult and time consuming process. This is especially true for helicopters. In this paper, two rapidly applicable and implementable methods to derive linear input filters that improve the simulator model are presented. The first method is based on model inversion, the second on feedback control. Both methods are evaluated in the time domain, compared to recorded helicopter flight test data, and assessed based on root mean square errors and the Qualification Test Guide bounds. The best results were achieved when using the first method.


1944 ◽  
Vol 7 (53) ◽  
pp. 279-294
Author(s):  
G. H. Menzies

2002 ◽  
Vol 12 (9) ◽  
pp. 315-315
Author(s):  
K. Cicak ◽  
K. O'Neill ◽  
R. E. Thorne

Below T=40 K, charge-density wave (CDW) transport in NbSe3 is characterized by two well-defined driving force thresholds ET and ET*. Between these thresholds the CDW moves extremely slowly with creep-like temperature and driving force dependencies. At the same time, the CDW exhibits coherent oscillations with a frequency proportional to the CDW current and having very narrow spectral widths, suggesting that the collective motion is temporally ordered. We have extended our initial work to doped crystals containing isoelectronic (Ta) and nonisoelectronic (Ti) impurities, and to crystals of different thicknesses. These experiments show that the qualitative features are extremely robust, and that the functional form of the creep velocity versus driving force and temperature is consistent across all samples for currents ranging over five orders of magnitude. The temperature dependence is consistent with processes having an energy comparable to the CDW gap, but the field and impurity dependencies are inconsistent with all predicted functional forms for creep in CDWs and related systems, and with our earlier picture of amplitude collapse at each impurity. We compare our results to measurements of creep-like behavior in other CDW and SDW systems, and discuss possible mechanisms.


2017 ◽  
Vol 2017 ◽  
pp. 1-12
Author(s):  
Lin Lin ◽  
Fang Wang ◽  
Shisheng Zhong

Prediction technology for aeroengine performance is significantly important in operational maintenance and safety engineering. In the prediction of engine performance, to address overfitting and underfitting problems with the approximation modeling technique, we derived a generalized approximation model that could be used to adjust fitting precision. Approximation precision was combined with fitting sensitivity to allow the model to obtain excellent fitting accuracy and generalization performance. Taking the Grey model (GM) as an example, we discussed the modeling approach of the novel GM based on fitting sensitivity, analyzed the setting methods and optimization range of model parameters, and solved the model by using a genetic algorithm. By investigating the effect of every model parameter on the prediction precision in experiments, we summarized the change regularities of the root-mean-square errors (RMSEs) varying with the model parameters in novel GM. Also, by analyzing the novel ANN and ANN with Bayesian regularization, it is concluded that the generalized approximation model based on fitting sensitivity can achieve a reasonable fitting degree and generalization ability.


2017 ◽  
Vol 10 (1) ◽  
pp. 155-165 ◽  
Author(s):  
Wengang Zhang ◽  
Guirong Xu ◽  
Yuanyuan Liu ◽  
Guopao Yan ◽  
Dejun Li ◽  
...  

Abstract. This paper is to investigate the uncertainties of microwave radiometer (MWR) retrievals in snow conditions and also explore the discrepancies of MWR retrievals in zenith and off-zenith observations. The MWR retrievals were averaged in a ±15 min period centered at sounding times of 00:00 and 12:00 UTC and compared with radiosonde observations (RAOBs). In general, the MWR retrievals have a better correlation with RAOB profiles in off-zenith observations than in zenith observations, and the biases (MWR observations minus RAOBs) and root mean square errors (RMSEs) between MWR and RAOB are also clearly reduced in off-zenith observations. The biases of temperature, relative humidity, and vapor density decrease from 4.6 K, 9 %, and 1.43 g m−3 in zenith observations to −0.6 K, −2 %, and 0.10 g m−3 in off-zenith observations, respectively. The discrepancies between MWR retrievals and RAOB profiles by altitude present the same situation. Cases studies show that the impact of snow on accuracies of MWR retrievals is more serious in heavy snowfall than in light snowfall, but off-zenith observation can mitigate the impact of snowfall. The MWR measurements become less accurate in snowfall mainly due to the retrieval algorithm, which does not consider the effect of snow, and the accumulated snow on the top of the radome increases the signal noise of MWR measurements. As the snowfall drops away by gravity on the sides of the radome, the off-zenith observations are more representative of the atmospheric conditions for RAOBs.


Symmetry ◽  
2021 ◽  
Vol 14 (1) ◽  
pp. 23
Author(s):  
Yuping Li ◽  
Brady K. Quinn ◽  
Johan Gielis ◽  
Yirong Li ◽  
Peijian Shi

Many natural radial symmetrical shapes (e.g., sea stars) follow the Gielis equation (GE) or its twin equation (TGE). A supertriangle (three triangles arranged around a central polygon) represents such a shape, but no study has tested whether natural shapes can be represented as/are supertriangles or whether the GE or TGE can describe their shape. We collected 100 pieces of Koelreuteria paniculata fruit, which have a supertriangular shape, extracted the boundary coordinates for their vertical projections, and then fitted them with the GE and TGE. The adjusted root mean square errors (RMSEadj) of the two equations were always less than 0.08, and >70% were less than 0.05. For 57/100 fruit projections, the GE had a lower RMSEadj than the TGE, although overall differences in the goodness of fit were non-significant. However, the TGE produces more symmetrical shapes than the GE as the two parameters controlling the extent of symmetry in it are approximately equal. This work demonstrates that natural supertriangles exist, validates the use of the GE and TGE to model their shapes, and suggests that different complex radially symmetrical shapes can be generated by the same equation, implying that different types of biological symmetry may result from the same biophysical mechanisms.


Sign in / Sign up

Export Citation Format

Share Document