Source author record

Yaomin Zhao

Yaomin Zhao appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

physics.flu-dyn physics.comp-ph

Catalog footprint

What is connected

6works

2topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2026arXiv

Low-wavenumber wall pressure fluctuations in turbulent flows within concentric annular ducts

Compressible direct numerical simulations of turbulent channel flows in concentric annular ducts of height $2δ$ are performed to study the low-wavenumber wall pressure fluctuations (WPF) over cylindrical walls at a bulk Mach number $M_b = 0.4$ and bulk Reynolds number $Re_b=3000$. The radius of the inner cylinder $R$ is varied between $0.2δ$, $δ$, $2δ$ and $\infty$. As $R$ decreases, the one-point power spectral density of the WPF decreases at intermediate but increases at high frequencies. When $R$ decreases, the 1D (streamwise) wavenumber-frequency spectrum of the WPF decreases at high wavenumbers. At low wavenumbers, however, as $R$ reduces to $0.2δ$ the 1D wavenumber-frequency spectrum exhibits multiple spectral peaks whose strengths increase with frequency. Examination of the 2D wavenumber-frequency spectra shows that these represent acoustic duct modes that closely match theoretical predictions. The acoustic modes of higher radial orders exhibit increasingly high amplitude on the inner than on the outer walls. The low-wavenumber components of the $0$th-order (azimuthal) 2D wavenumber-frequency spectrum are of great importance in practice, and their magnitude increases as $R$ reduces; this increase is increasingly pronounced at higher frequencies. Analytical modelling and numerical validation show that this increase appears to arise from the ``geometric'' effects connected with the Green's function, and they are generated mainly by radial and azimuthal disturbances. Disturbances closer to the wall are shown to be increasingly important in WPF generation as $R$ reduces, which highlights a potential in WPF control using wall treatment on thin cylinders.

preprint2026arXiv

Progressive Mixture-of-Experts with autoencoder routing for continual RANS turbulence modelling

Developing Reynolds-averaged Navier-Stokes (RANS) turbulence models that remain accurate across diverse flow regimes remains a long-standing challenge. In this work, we propose a novel framework, termed the progressive mixture-of-experts (PMoE), designed to enable continual learning for RANS turbulence modelling. The framework employs a modular autoencoder-based router to associate each flow scenario with a specialised turbulence model, referred to as an expert. When an unseen flow regime cannot be adequately represented by the existing router and expert set, a new expert together with its routing component can be introduced at low cost, without modifying or degrading previously trained ones, thereby naturally avoiding catastrophic forgetting. The framework is applied to a range of flows with distinct physical characteristics, including baseline airfoil wakes, wall-attached flows, separated flows and corner-induced secondary flows. The resulting PMoE model effectively integrates multiple experts and achieves improved predictive accuracy across both seen and unseen test cases. Owing to sparse activation, model expansion does not incur additional computational cost during inference. The proposed framework therefore provides a scalable pathway towards lifelong-learning turbulence models for industrial computational fluid dynamics.

preprint2023arXiv

Turbulence Model Development based on a Novel Method Combining Gene Expression Programming with an Artificial Neural Network

Data-driven methods are widely used to develop physical models, but there still exist limitations that affect their performance, generalizability and robustness. By combining gene expression programming (GEP) with artificial neural network (ANN), we propose a novel method for symbolic regression called the gene expression programming neural network (GEPNN). In this method, candidate expressions generated by evolutionary algorithms are transformed between the GEP and ANN structures during training iterations, and efficient and robust convergence to accurate models is achieved by combining the GEP's global searching and the ANN's gradient optimization capabilities. In addition, sparsity-enhancing strategies have been introduced to GEPNN to improve the interpretability of the trained models. The GEPNN method has been tested for finding different physical laws, showing improved convergence to models with precise coefficients. Furthermore, for large-eddy simulation of turbulence, the subgrid-scale stress model trained by GEPNN significantly improves the prediction of turbulence statistics and flow structures over traditional models, showing advantages compared to the existing GEP and ANN methods in both a priori and a posteriori tests.

preprint2021arXiv

Multi-Objective CFD-Driven Development of Coupled Turbulence Closure Models

This paper introduces two novel concepts in data-driven turbulence modeling that enable the simultaneous development of multiple closure models and the training towards multiple objectives. The concepts extend the evolutionary framework by Weatheritt and Sandberg (2016), which derives interpretable and implementation-ready expressions from high-fidelity simulation data. By assigning a shared fitness value to the evolved closure models and utilizing the CFD-driven training approach by Zhao et al. (2020), the multi-expression training concept introduced here is able to account for the coupling between the trained models, i.e. Reynolds stress anisotropy, turbulent heat flux and turbulence production correction models. As a second concept, a multi-objective optimization algorithm is applied to the framework. The extension yields a diverse set of candidate models and allows a trade-off between the training objectives after analyzing the training results. In this study, the novel concepts are applied to a benchmark periodic hills case and a vertical natural convection flow. The predictions of mean flow quantities are improved compared to decoupled training strategies with distinct and robust improvements for strongly coupled momentum and thermal fields. The coupled training of closure models and the balancing of multiple training objectives are considered important capabilities on the path towards generalized data-driven turbulence models.

preprint2020arXiv

RANS Turbulence Model Development using CFD-Driven Machine Learning

This paper presents a novel CFD-driven machine learning framework to develop Reynolds-averaged Navier-Stokes (RANS) models. The CFD-driven training is an extension of the gene expression programming method (Weatheritt and Sandberg, 2016), but crucially the fitness of candidate models is now evaluated by running RANS calculations in an integrated way, rather than using an algebraic function. Unlike other data-driven methods that fit the Reynolds stresses of trained models to high-fidelity data, the cost function for the CFD-driven training can be defined based on any flow feature from the CFD results. This extends the applicability of the method especially when the training data is limited. Furthermore, the resulting model, which is the one providing the most accurate CFD results at the end of the training, inherently shows good performance in RANS calculations. To demonstrate the potential of this new method, the CFD-driven machine learning approach is applied to model development for wake mixing in turbomachines. A new model is trained based on a high-pressure turbine case and then tested for three additional cases, all representative of modern turbine nozzles. Despite the geometric configurations and operating conditions being different among the cases, the predicted wake mixing profiles are significantly improved in all of these a posteriori tests. Moreover, the model equation is explicitly given and available for analysis, thus it could be deduced that the enhanced wake prediction is predominantly due to the extra diffusion introduced by the CFD-driven model.

preprint2016arXiv

Evolution of material surfaces in the temporal transition in channel flow

We report a Lagrangian study on the evolution of material surfaces in the K-type temporal transitional channel flow. Based on the Eulerian velocity field from the DNS, a backward-particle-tracking method is applied to solve the transport equation of the Lagrangian scalar field, and then the iso-surfaces of the Lagrangian field can be extracted as material surfaces in the evolution. Three critical issues for Lagrangian investigations on the evolution of coherent structures using material surfaces are addressed. First, the initial scalar field is uniquely determined based on proposed criteria, so that the initial material surfaces can be approximated as vortex surfaces, and keep invariant in the initial laminar state. Second, the evolution of typical material surfaces initially from different wall distances is presented, and then the influential material surface with the maximum deformation is identified. Large vorticity variations with the maximum curvature growth of vortex lines are also observed on this surface. Moreover, crucial events in the transition can be characterized in a Lagrangian approach by conditional statistics on the material surfaces. Finally, the influential material surface, which is initially a vortex surface, is demonstrated as a surrogate of the vortex surface before significant topological changes of vortical structures. Therefore, this material surface can be used to elucidate the continuous evolution of vortical structures in transitional wall-bounded flows in a Lagrangian perspective. The evolution of the influential material surface is divided into three stages: the formation of a triangular bulge from an initially disturbed streamwise-spanwise sheet, rolling up of the vortex sheet near the bulge ridges with the vorticity intensification, and the generation and evolution of signature hairpin-like structures with self-induced dynamics of vortex filaments.