Efficient low rank model order reduction of vibroacoustic problems under stochastic loads

Yannik Hüpel,Ulrich Römer,Matthias Bollhöfer,Sabine Langer
2024-08-16
Abstract:This contribution combines a low-rank matrix approximation through Singular Value Decomposition (SVD) with second-order Krylov subspace-based Model Order Reduction (MOR), in order to efficiently propagate input uncertainties through a given vibroacoustic model. The vibroacoustic model consists of a plate coupled to a fluid into which the plate radiates sound due to a turbulent boundary layer excitation. This excitation is subject to uncertainties due to the stochastic nature of the turbulence and the computational cost of simulating the coupled problem with stochastic forcing is very high. The proposed method approximates the output uncertainties in an efficient way, by reducing the evaluation cost of the model in terms of DOFs and samples by using the factors of the SVD low-rank approximation directly as input for the MOR algorithm. Here, the covariance matrix of the vector of unknowns can efficiently be approximated with only a fraction of the original number of evaluations. Therefore, the approach is a promising step to further reducing the computational effort of large-scale vibroacoustic evaluations.
Computational Engineering, Finance, and Science
What problem does this paper attempt to address?
### What problem does this paper attempt to solve? This paper aims to solve the computational efficiency problem of input uncertainty propagation in vibro - acoustic models under random loads. Specifically, the research background involves the evaluation of aircraft cabin noise. Due to the randomness and complexity of excitation by the turbulent boundary layer (TBL), traditional high - fidelity wave - analytical vibro - acoustic models are very computationally expensive when considering uncertainties and material property variations. #### Main problems: 1. **High computational cost**: When uncertainties in excitation or material properties need to be considered, the computational amount of traditional methods is very large, making vibro - acoustic evaluation in the early design stage infeasible. 2. **Efficient uncertainty quantification**: In order to effectively evaluate the impact of these uncertainties on the system response, a method that can significantly reduce the computational cost needs to be developed. 3. **Multi - query tasks**: For tasks involving multiple solutions (such as different frequency and parameter combinations), efficient numerical methods are required to accelerate the computational process. #### Solutions: The paper proposes a method that combines low - rank matrix approximation (through singular value decomposition SVD) and a model order reduction method (MOR) based on the second - order Krylov subspace to efficiently propagate input uncertainties. The specific steps are as follows: 1. **Low - rank approximation**: Perform low - rank approximation on input uncertainties through SVD to reduce the degrees of freedom and the number of samples required for model evaluation. 2. **Model order reduction**: Use the factors of the low - rank approximation as the input of the MOR algorithm to further reduce the computational cost. 3. **Efficient uncertainty propagation**: Through the above methods, the output uncertainty can be efficiently approximated and the computational burden of the original model can be significantly reduced. This method not only reduces the number of system solutions at each frequency point but also can significantly reduce the computational complexity while maintaining sufficient accuracy. The ultimate goal is to provide an efficient and reliable solution for uncertainty quantification in large - scale vibro - acoustic models. ### Formula summary: - **SVD decomposition**: \[ \Sigma_f = U S V^H \approx U_l S_l V_l^H \] where \(U\) and \(V\) represent the left and right singular vectors respectively, \(S\) is the singular value matrix, and the subscript \(l\) represents low - rank approximation. - **Covariance matrix calculation**: \[ \Sigma_x \approx U_x S_x V_x^T \] where \(S_x\) comes from the SVD decomposition of \(\Sigma_f\), and satisfies: \[ A U_x = U, \quad A V_x = V \] - **Error measurement**: \[ \text{err} = \max(||p_{\text{FOM}}| - |p_{\text{ROM}}||) \] This is an absolute error indicator used to quantitatively evaluate the accuracy of the ROM. Through these methods, the paper shows how to significantly improve the computational efficiency of uncertainty propagation in vibro - acoustic models while ensuring accuracy.