Sampling low-fidelity outputs for estimation of high-fidelity density and its tails

Minji Kim,Vladas Pipiras,Kevin O'Connor,Themistoklis Sapsis

2024-02-28

Abstract:In a multifidelity setting, data are available under the same conditions from two (or more) sources, e.g. computer codes, one being lower-fidelity but computationally cheaper, and the other higher-fidelity and more expensive. This work studies for which low-fidelity outputs, one should obtain high-fidelity outputs, if the goal is to estimate the probability density function of the latter, especially when it comes to the distribution tails and extremes. It is suggested to approach this problem from the perspective of the importance sampling of low-fidelity outputs according to some proposal distribution, combined with special considerations for the distribution tails based on extreme value theory. The notion of an optimal proposal distribution is introduced and investigated, in both theory and simulations. The approach is motivated and illustrated with an application to estimate the probability density function of record extremes of ship motions, obtained through two computer codes of different fidelities.

Methodology

What problem does this paper attempt to address?

The paper attempts to address the problem of how to utilize low-cost, low-fidelity outputs to estimate the probability density function (PDF) of high-fidelity outputs in a multi-fidelity setting, particularly focusing on the tail distribution. Specifically:

1. **Research Background**: In ship motion simulation, there are two types of computer codes: one is LAMP, which has high fidelity but high computational cost; the other is SC, which has low fidelity but high computational efficiency. Data generated by these two codes under the same conditions can be used to study the extreme values of ship motion.

2. **Core Issue**: How to select low-fidelity outputs (such as SC) to generate corresponding high-fidelity outputs (such as LAMP) in order to better estimate the probability density function of high-fidelity outputs, especially in the tail distribution.

3. **Methodology**: The paper proposes an importance sampling method that selects specific samples from low-fidelity outputs and combines them with extreme value theory to optimize the estimation of the probability density function of high-fidelity outputs. This method aims to find the optimal importance sampling distribution to improve the accuracy of the estimation.

4. **Objective**: The paper aims to develop an effective method to estimate the probability density function of high-fidelity outputs, particularly to improve the estimation accuracy when dealing with extreme values. This helps to better understand the extreme conditions of ship motion in practical applications.