In this work, we present a statistical treatment of stress-life (S-N) data drawn from a collection of records of fatigue experiments that were performed on 75S-T6 aluminum alloys. Our main objective is to predict the fatigue life of materials by providing a systematic approach to model calibration, model selection and model ranking with reference to S-N data. To this purpose, we consider fatigue-limit models and random fatigue-limit models that are specially designed to allow the treatment of the run-outs (right-censored data). We first fit the models to the data by maximum likelihood methods and estimate the quantiles of the life distribution of the alloy specimen. To assess the robustness of the estimation of the quantile functions, we obtain bootstrap confidence bands by stratified resampling with respect to the cycle ratio. We then compare and rank the models by classical measures of fit based on information criteria. We also consider a Bayesian approach that provides, under the prior distribution of the model parameters selected by the user, their simulation-based posterior distributions. We implement and apply Bayesian model comparison methods, such as Bayes factor ranking and predictive information criteria based on cross-validation techniques under various a priori scenarios

The assessment of the thermal properties of walls is essential for accurate building energy simulations that are needed to make effective energy-saving policies. These properties are usually investigated through in-situ measurements of temperature and heat flux over extended time periods. The one-dimensional heat equation with unknown Dirichlet boundary conditions is used to model the heat transfer process through the wall. In [F. Ruggeri, Z. Sawlan, M. Scavino, R. Tempone, A hierarchical Bayesian setting for an inverse problem in linear parabolic PDEs with noisy boundary conditions, Bayesian Analysis 12 (2) (2017) 407-433], it was assessed the uncertainty about the thermal diffusivity parameter using different synthetic data sets. In this work, we adapt this methodology to an experimental study conducted in an environmental chamber, with measurements recorded every minute from temperature probes and heat flux sensors placed on both sides of a solid brick wall over a five-day period. The observed time series are locally averaged, according to a smoothing procedure determined by the solution of a criterion function optimization problem, to fit the required set of noise model assumptions. Therefore, after preprocessing, we can reasonably assume that the temperature and the heat flux measurements have stationary Gaussian noise and we can avoid working with full covariance matrices. The results show that our technique reduces the bias error of the estimated parameters when compared to other approaches. Finally, we compute the information gain under two experimental setups to recommend how the user can efficiently determine the duration of the measurement campaign and the range of the external temperature oscillation.

In this work we develop a Bayesian setting to infer unknown parameters in initial-boundary value problems related to linear parabolic partial differential equations. We realistically assume that the boundary data are noisy, for a given prescribed initial condition. We show how to derive the joint likelihood function for the forward problem, given some measurements of the solution field subject to Gaussian noise. Given Gaussian priors for the time-dependent Dirichlet boundary values, we analytically marginalize the joint likelihood using the linearity of the equation. Our hierarchical Bayesian approach is fully implemented in an example that involves the heat equation. In this example, the thermal diffusivity is the unknown parameter. We assume that the thermal diffusivity parameter can be modeled a priori through a lognormal random variable or by means of a spacedependent stationary lognormal random field. Synthetic data are used to test the inference. We exploit the behavior of the non-normalized log posterior distribution of the thermal diffusivity. Then, we use the Laplace method to obtain an approximated Gaussian posterior and therefore avoid costly Markov Chain Monte Carlo computations. Expected information gains and predictive posterior densities for observable quantities are numerically estimated using Laplace approximation for different experimental setups.

In this work we propose a stochastic model for estimating the occurrence of crack initiations on the surface of metallic specimens in fatigue problems that can be applied to a general class of geometries. The stochastic model is based on spatial Poisson processes with intensity function that combines stresslife (S-N) curves with averaged effective stress, σ ∆ eff (x), which is computed after solving numerically the linear elasticity equations on the specimen domains using finite element methods. Here, ∆ is a parameter that characterizes the size of the neighbors covering the domain boundary. The averaged effective stress, parameterized by ∆, maps the stress tensor to a scalar field upon the specimen domain. Data from fatigue experiments on notched and unnotched sheet specimens of 75S-T6 aluminum alloys are used to calibrate the model parameters for the individual data sets and their combination. Bayesian and classical approaches are applied to estimate the survival-probability function for any specimen tested under a prescribed fatigue experimental setup. Our proposed model can predict the initiation of cracks in specimens made from the same material with new geometries.

In this work, we present the ensemble-marginalized Kalman filter (EnMKF), a sequential algorithm analogous to our previously proposed approach [1,2], for estimating the state and parameters of linear parabolic partial differential equations in initial-boundary value problems when the boundary data are noisy. We apply EnMKF to infer the thermal properties of building walls and to estimate the corresponding heat flux from real and synthetic data. Compared with a modified Ensemble Kalman Filter (EnKF) that is not marginalized, EnMKF reduces the bias error, avoids the collapse of the ensemble without needing to add inflation, and converges to the mean field posterior using 50% or less of the ensemble size required by EnKF. According to our results, the marginalization technique in EnMKF is key to performance improvement with smaller ensembles at any fixed time.