No Arabic abstract
We develop a fast multi-fidelity modeling method for very complex correlations between high- and low-fidelity data by working in modal space to extract the proper correlation function. We apply this method to infer the amplitude of motion of a flexible marine riser in cross-flow, subject to vortex-induced vibrations (VIV). VIV are driven by an absolute instability in the flow, which imposes a frequency (Strouhal) law that requires a matching with the impedance of the structure; this matching is easily achieved because of the rapid parametric variation of the added mass force. As a result, the wavenumber of the riser spatial response is within narrow bands of uncertainty. Hence, an error in wavenumber prediction can cause significant phase-related errors in the shape of the amplitude of response along the riser, rendering correlation between low- and high-fidelity data very complex. Working in modal space as outlined herein, dense data from low-fidelity data, provided by the semi-empirical computer code VIVA, can correlate in modal space with few high-fidelity data, obtained from experiments or fully-resolved CFD simulations, to correct both phase and amplitude and provide predictions that agree very well overall with the correct shape of the amplitude response. We also quantify the uncertainty in the prediction using Bayesian modeling and exploit this uncertainty to formulate an active learning strategy for the best possible location of the sensors providing the high fidelity measurements.
Slender marine structures such as deep-water marine risers are subjected to currents and will normally experience Vortex Induced Vibrations (VIV), which can cause fast accumulation of fatigue damage. The ocean current is often three-dimensional (3D), i.e., the direction and magnitude of the current vary throughout the water column. Today, semi-empirical tools are used by the industry to predict VIV induced fatigue on risers. The load model and hydrodynamic parameters in present VIV prediction tools are developed based on two-dimensional (2D) flow conditions, as it is challenging to consider the effect of 3D flow along the risers. Accordingly, the current profiles must be purposely made 2D during the design process, which leads to significant uncertainty in the prediction results. Further, due to the limitations in the laboratory, VIV model tests are mostly carried out under 2D flow conditions and thus little experimental data exist to document VIV response of riser subjected to varying directions of the current. However, a few experiments have been conducted with 3D current. We have used results from one of these experiments to investigate how well 1) traditional and 2) an alternative method based on a data driven prediction can describe VIV in 3D currents. Data driven modelling is particularly suited for complicated problems with many parameters and non-linear relationships. We have applied a data clustering algorithm to the experimental 3D flow data in order to identify measurable parameters that can influence responses. The riser responses are grouped based on their statistical characteristics, which relate to the direction of the flow. Furthermore we fit a random forest regression model to the measured VIV response and compare its performance with the predictions of existing VIV prediction tools (VIVANA-FD).
We present a scale-bridging approach based on a multi-fidelity (MF) machine-learning (ML) framework leveraging Gaussian processes (GP) to fuse atomistic computational model predictions across multiple levels of fidelity. Through the posterior variance of the MFGP, our framework naturally enables uncertainty quantification, providing estimates of confidence in the predictions. We used Density Functional Theory as high-fidelity prediction, while a ML interatomic potential is used as the low-fidelity prediction. Practical materials design efficiency is demonstrated by reproducing the ternary composition dependence of a quantity of interest (bulk modulus) across the full aluminum-niobium-titanium ternary random alloy composition space. The MFGP is then coupled to a Bayesian optimization procedure and the computational efficiency of this approach is demonstrated by performing an on-the-fly search for the global optimum of bulk modulus in the ternary composition space. The framework presented in this manuscript is the first application of MFGP to atomistic materials simulations fusing predictions between Density Functional Theory and classical interatomic potential calculations.
Spatio-temporally extended nonlinear systems often exhibit a remarkable complexity in space and time. In many cases, extensive datasets of such systems are difficult to obtain, yet needed for a range of applications. Here, we present a method to generate synthetic time series or fields that reproduce statistical multi-scale features of complex systems. The method is based on a hierarchical refinement employing transition probability density functions (PDFs) from one scale to another. We address the case in which such PDFs can be obtained from experimental measurements or simulations and then used to generate arbitrarily large synthetic datasets. The validity of our approach is demonstrated at the example of an experimental dataset of high Reynolds number turbulence.
Calibrating a Reynolds-averaged Navier-Stokes (RANS) model against data leads to an improvement. Determining {it a priori} if such an improvement generalizes to flows outside the calibration data is an outstanding challenge. This work attempts to address this challenge via global epistemic Uncertainty Quantification (UQ). Unlike the available epistemic UQ methods that are local and tell us a models uncertainty at one specific flow condition, the global epistemic UQ method presented in this work tells us also whether a perturbation of the original model would generalize. Specifically, the global epistemic UQ method evaluates a potential improvement in terms of its effectiveness and inconsistency. Any improvement can be put in one of the following four quadrants: first, high effectiveness, low inconsistency; second, high effectiveness, high inconsistency; third, low effectiveness, low inconsistency; and fourth, low effectiveness, high inconsistency. An improvement would generalize if and only if it is in the high effectiveness and low inconsistency quadrant. To demonstrate the concept, we apply the global epistemic UQ to full Reynolds stress modeling of a stratified shear layer. The global epistemic UQ results point to a model coefficient in the pressure-strain correlation closure (among others) as effective and consistent for predicting the quantity of interest of shear layers growth. We calibrate the model coefficient such that our RANS matches direct numerical simulation data at one flow condition. We show that the calibrated model generalizes to several other test flow conditions. On the other hand, when calibrating a high inconsistency term, we get a model that works at only the calibrated condition.
We present an accurate and robust numerical method to track quantized vortex lines in a superfluid described by the Gross-Pitaevskii equation. By utilizing the pseudo-vorticity field of the associated complex scalar order parameter of the superfluid, we are able to track the topological defects of the superfluid and reconstruct the vortex lines which correspond to zeros of the field. Throughout, we assume our field is periodic to allow us to make extensive use of the Fourier representation of the field and its derivatives in order to retain spectral accuracy. We present several case studies to test the precision of the method which include the evaluation of the curvature and torsion of a torus vortex knot, and the measurement of the Kelvin wave spectrum of a vortex line and a vortex ring. The method we present makes no a-priori assumptions on the geometry of the vortices and is therefore applicable to a wide range of systems such as a superfluid in a turbulent state that is characterised by many vortex rings coexisting with sound waves. This allows us to track the positions of the vortex filaments in a dense turbulent vortex tangle and extract statistical information about the distribution of the size of the vortex rings and the inter-vortex separations. In principle, the method can be extended to track similar topological defects arising in other physical systems.