CFD-driven Symbolic Identification of Algebraic Reynolds-Stress Models

TL;DR

This paper introduces a CFD-driven symbolic identification algorithm for learning explicit algebraic Reynolds-stress models from high-fidelity data, enabling more accurate turbulence modeling in complex flows.

Contribution

The paper presents a novel CFD-driven optimization approach for discovering explicit algebraic Reynolds-stress models, overcoming limitations of previous data-restricted methods.

Findings

Discovered models outperform baseline linear eddy viscosity models in predicting separated flows.

The method effectively reduces the search space using sensitivity analysis and response surface approximation.

Models generalize well to unseen flow configurations with higher Reynolds numbers.

Abstract

A CFD-driven deterministic symbolic identification algorithm for learning explicit algebraic Reynolds-stress models (EARSM) from high-fidelity data is developed building on the frozen-training SpaRTA algorithm of [1]. Corrections for the Reynolds stress tensor and the production of transported turbulent quantities of a baseline linear eddy viscosity model (LEVM) are expressed as functions of tensor polynomials selected from a library of candidate functions. The CFD-driven training consists in solving a blackbox optimization problem in which the fitness of candidate EARSM models is evaluated by running RANS simulations. Unlike the frozen-training approach, the proposed methodology is not restricted to data sets for which full fields of high-fidelity data are available. However, the solution of a high-dimensional expensive blackbox function optimization problem is required. Several steps are then undertaken to reduce the associated computational burden. First, a sensitivity analysis is used to identify the most influential terms and to reduce the dimensionality of the search space. Afterwards, the Constrained Optimization using Response Surface (CORS) algorithm, which approximates the black-box cost function using a response surface constructed from a limited number of CFD solves, is used to find the optimal model parameters. Model discovery and cross-validation is performed for three configurations of 2D turbulent separated flows in channels of variable section using different sets of training data to show the flexibility of the method. The discovered models are then applied to the prediction of an unseen 2D separated flow with higher Reynolds number and different geometry. The predictions for the new case are shown to be not only more accurate than the baseline LEVM, but also of a multi-purpose EARSM model derived from purely physical arguments.