Search

Within the context of machine learning-based closure mappings for RANS turbulence modelling, physical realizability is often enforced using ad-hoc postprocessing of the predicted anisotropy tensor. In this study, we address the realizability issue via a new physics-based loss function that penalizes non-realizable results during training, thereby embedding a preference for realizable predictions into the model. Additionally, we propose a new framework for data-driven turbulence modelling which retains the stability and conditioning of optimal eddy viscosity-based approaches while embedding equivariance. Several modifications to the tensor basis neural network to enhance training and testing stability are proposed. We demonstrate the conditioning, stability, and generalization of the new framework and model architecture on three flows: flow over a flat plate, flow over periodic hills, and flow through a square duct. The realizability-informed loss function is demonstrated to significantly increase the number of realizable predictions made by the model when generalizing to a new flow configuration. Altogether, the proposed framework enables the training of stable and equivariant anisotropy mappings, with more physically realizable predictions on new data. We make our code available for use and modification by others., Comment: 44 pages, 21 figures

In this paper, we present novel identification strategies to develop a unified framework for vortex-induced vibration (VIV) prediction based on the general semi-empirical wake oscillator. Greybox nonlinear system identification method accompanying high-fidelity computational fluid dynamics (CFD) and/or experimental data could be applied for the identification process. The proposed template of general wake oscillators contains low- to high-order damping terms to be identified for characterizing the possible flow dynamics. Two different strategies, including individual identification of single wake oscillator and overall identification of coupled VIV control equations, are proposed. VIV system consisting of an elastically-mounted circular cylinder submerged in laminar flow at Reynold number of 100 is considered. Both strategies have been tested and have exhibited high accuracy. The second strategy, i.e., overall identification of coupled VIV control equations, would be more suitable for the future framework owing that its training process considers the effect of fluid damping. A detailed mathematical introduction to future works on framework development covering the wide Reynold number range is addressed. The proposed unified framework is a landmark update of past wake oscillators both in terms of prediction accuracy and physical principles and has considerable research significance and practical engineering value., Comment: This version is not intended to be sent for peer review and is merely a summary of the preliminary work

Industrial simulations of turbulent flows often rely on Reynolds-averaged Navier-Stokes (RANS) turbulence models, which contain numerous closure coefficients that need to be calibrated. In this work, we address this issue by proposing a semi-automated calibration of these coefficients using a new framework (referred to as turbo-RANS) based on Bayesian optimization. We introduce the generalized error and default coefficient preference (GEDCP) objective function, which can be used with integral, sparse, or dense reference data for the purpose of calibrating RANS turbulence closure model coefficients. Then, we describe a Bayesian optimization-based algorithm for conducting the calibration of these model coefficients. An in-depth hyperparameter tuning study is conducted to recommend efficient settings for the turbo-RANS optimization procedure. We demonstrate that the performance of the $k$-$\omega$ shear stress transport (SST) and Generalized $k$-$\omega$ (GEKO) turbulence models can be efficiently improved via turbo-RANS, for three example cases: predicting the lift coefficient of an airfoil; predicting the velocity and turbulent kinetic energy fields for a separated flow; and, predicting the wall pressure coefficient distribution for flow through a converging-diverging channel. This work is the first to propose and provide an open-source black-box calibration procedure for turbulence model coefficients based on Bayesian optimization. We propose a data-flexible objective function for the calibration target. Our open-source implementation of the turbo-RANS framework includes OpenFOAM, Ansys Fluent, STAR-CCM+, and solver-agnostic templates for user application., Comment: 24 pages, 14 figures. Published in International Journal of Numerical Methods for Heat & Fluid Flow (https://doi.org/10.1108/HFF-12-2023-0726)

Several machine learning frameworks for augmenting turbulence closure models have been recently proposed. However, the generalizability of an augmented turbulence model remains an open question. We investigate this question by systematically varying the training and test sets of several models. An optimal three-term tensor basis expansion is used to develop a model-agnostic data-driven turbulence closure approximation. Then, hyperparameter optimization was performed for a random forest, a neural network, and an eXtreme Gradient Boosting (XGBoost) model. We recommend XGBoost for data-driven turbulence closure modelling owing to its low-tuning cost and good performance. We also find that machine learning models generalize well to new parametric variations of flows seen in the training dataset, but lack generalizability to new flow types. This generalizability gap suggests that machine learning methods are most suited for developing specialized models for a given flow type, a problem often encountered in industrial applications., Comment: Article: 40 pages, 15 figures, supplemental (ancillary) file: 29 pages, 60 figures. Accepted in International Journal of Computational Fluid Dynamics

Despite well-known limitations of Reynolds-averaged Navier-Stokes (RANS) simulations, this methodology remains the most widely used tool for predicting many turbulent flows, due to computational efficiency. Machine learning is a promising approach to improve the accuracy of RANS simulations. One major area of improvement is using machine learning models to represent the complex relationship between the mean flow field gradients and the Reynolds stress tensor. In the present work, modifications to improve the stability of previous optimal eddy viscosity approaches for RANS simulations are presented and evaluated. The optimal eddy viscosity is reformulated with a non-negativity constraint, which promotes numerical stability. We demonstrate that the new formulation of the optimal eddy viscosity improves the conditioning of the RANS equations for a periodic hills test case. To demonstrate the suitability of this proportional/orthogonal tensor decomposition for use in a physics-informed data-driven turbulence closure, we use two neural networks (structured on this specific tensor decomposition which is incorporated as an inductive bias into the network design) to predict the newly reformulated linear and non-linear parts of the Reynolds stress tensor. Injecting these network model predictions for the Reynolds stresses into a RANS simulation improves predictions of the velocity field, even when compared to a sophisticated (state of the art) physics-based turbulence closure model. Finally, we apply SHAP (SHapley Additive exPlanations) values to obtain insights from the learned representation for the inner workings of the neural network used to predict the optimal eddy viscosity from the input feature data., Comment: 18 pages, 23 figures. Article submitted to Physics of Fluids

The recent surge in machine learning augmented turbulence modelling is a promising approach for addressing the limitations of Reynolds-averaged Navier-Stokes (RANS) models. This work presents the development of the first open-source dataset, curated and structured for immediate use in machine learning augmented turbulence closure modelling. The dataset features a variety of RANS simulations with matching direct numerical simulation (DNS) and large-eddy simulation (LES) data. Four turbulence models are selected to form the initial dataset: $k$-$\varepsilon$, $k$-$\varepsilon$-$\phi_t$-$f$, $k$-$\omega$, and $k$-$\omega$ SST. The dataset consists of 29 cases per turbulence model, for several parametrically sweeping reference DNS/LES cases: periodic hills, square duct, parametric bumps, converging-diverging channel, and a curved backward-facing step. At each of the 895,640 points, various RANS features with DNS/LES labels are available. The feature set includes quantities used in current state-of-the-art models, and additional fields which enable the generation of new feature sets. The dataset reduces effort required to train, test, and benchmark new models. The dataset is available at https://doi.org/10.34740/kaggle/dsv/2044393 ., Comment: 21 pages, 20 figures, 15 tables. Article yet to be submitted to a journal

In this paper, a simple field function is presented for facilitating the solution of complex and dynamic fluid-solid systems on Cartesian grids with interface-resolved fluid-fluid, fluid-solid, and solid-solid interactions. For a Cartesian-grid-discretized computational domain segmented by a set of solid bodies, this field function explicitly tracks each subdomain with multiple resolved interfacial node layers. As a result, the presented field function enables low-memory-cost multidomain node mapping, efficient node remapping, fast collision detection, and expedient surface force integration. Implementation algorithms for the field function and its described functionalities are also presented. Equipped with a deterministic multibody collision model, numerical experiments involving fluid-solid systems with flow conditions ranging from subsonic to supersonic states are conducted to validate and illustrate the applicability of the proposed field function.

In this paper, a novel immersed boundary method is developed, validated, and applied. Through devising a second-order three-step flow reconstruction scheme, the proposed method is able to enforce the Dirichlet, Neumann, Robin, and Cauchy boundary conditions in a straightforward and consistent manner. Equipped with a fluid-solid coupling framework that integrates high-order temporal and spatial discretization schemes, numerical experiments concerning flow involving stationary and moving objects, convex and concave geometries, no-slip and slip wall boundary conditions, as well as subsonic and supersonic motions are conducted to validate the method. It is demonstrated that the proposed method can provide efficient, accurate, and robust boundary treatment for solving flow with arbitrarily irregular and moving geometries on Cartesian grids.