Skip to main content Accessibility help
×
Home

Reynolds-averaged Navier–Stokes equations with explicit data-driven Reynolds stress closure can be ill-conditioned

Published online by Cambridge University Press:  29 April 2019

Jinlong Wu
Affiliation:
Kevin T. Crofton Department of Aerospace and Ocean Engineering, Virginia Tech, Blacksburg, VA 24060, USA
Heng Xiao
Affiliation:
Kevin T. Crofton Department of Aerospace and Ocean Engineering, Virginia Tech, Blacksburg, VA 24060, USA
Rui Sun
Affiliation:
Kevin T. Crofton Department of Aerospace and Ocean Engineering, Virginia Tech, Blacksburg, VA 24060, USA
Qiqi Wang
Affiliation:
Department of Aeronautics and Astronautics, MIT, Cambridge, MA 02139, USA
Corresponding
E-mail address:

Abstract

Reynolds-averaged Navier–Stokes (RANS) simulations with turbulence closure models continue to play important roles in industrial flow simulations. However, the commonly used linear eddy-viscosity models are intrinsically unable to handle flows with non-equilibrium turbulence (e.g. flows with massive separation). Reynolds stress models, on the other hand, are plagued by their lack of robustness. Recent studies in plane channel flows found that even substituting Reynolds stresses with errors below 0.5 % from direct numerical simulation databases into RANS equations leads to velocities with large errors (up to 35 %). While such an observation may have only marginal relevance to traditional Reynolds stress models, it is disturbing for the recently emerging data-driven models that treat the Reynolds stress as an explicit source term in the RANS equations, as it suggests that the RANS equations with such models can be ill-conditioned. So far, a rigorous analysis of the condition of such models is still lacking. As such, in this work we propose a metric based on local condition number function for a priori evaluation of the conditioning of the RANS equations. We further show that the ill-conditioning cannot be explained by the global matrix condition number of the discretized RANS equations. Comprehensive numerical tests are performed on turbulent channel flows at various Reynolds numbers and additionally on two complex flows, i.e. flow over periodic hills, and flow in a square duct. Results suggest that the proposed metric can adequately explain observations in previous studies, i.e. deteriorated model conditioning with increasing Reynolds number and better conditioning of the implicit treatment of the Reynolds stress compared to the explicit treatment. This metric can play critical roles in the future development of data-driven turbulence models by enforcing the conditioning as a requirement on these models.

Type
JFM Papers
Copyright
© 2019 Cambridge University Press 

Access options

Get access to the full version of this content by using one of the access options below.

References

Basara, B. & Jakirlic, S. 2003 A new hybrid turbulence modelling strategy for industrial CFD. Intl J. Numer. Meth. Fluids 42 (1), 89116.CrossRefGoogle Scholar
Bernardini, M., Pirozzoli, S. & Orlandi, P. 2014 Velocity statistics in turbulent channel flow up to Re = 4000. J. Fluid Mech. 742, 171191.CrossRefGoogle Scholar
Breuer, M., Peller, N., Rapp, C. & Manhart, M. 2009 Flow over periodic hills: numerical and experimental study in a wide range of Reynolds numbers. Comput. Fluids 38 (2), 433457.CrossRefGoogle Scholar
Chandrasekaran, S. & Ipsen, I. C. F. 1995 On the sensitivity of solution components in linear systems of equations. SIAM J. Matrix Anal. Applics. 16 (1), 93112.CrossRefGoogle Scholar
Debnath, L. & Mikusiński, P. 2005 Hilbert Spaces with Applications. Academic Press.Google Scholar
Del Alamo, J. C. & Jiménez, J. 2003 Spectra of the very large anisotropic scales in turbulent channels. Phys. Fluids 15 (6), L41L44.CrossRefGoogle Scholar
Gamahara, M. & Hattori, Y. 2017 Searching for turbulence models by artificial neural network. Phys. Rev. Fluids 2 (5), 054604.CrossRefGoogle Scholar
Geneva, N. & Zabaras, N. 2019 Quantifying model form uncertainty in Reynolds-averaged turbulence models with Bayesian deep neural networks. J. Comput. Phys. 383, 125147.CrossRefGoogle Scholar
Hamlington, P. E. & Dahm, W. J. A. 2008 Reynolds stress closure for nonequilibrium effects in turbulent flows. Phys. Fluids 20 (11), 115101.CrossRefGoogle Scholar
Hamlington, P. E. & Ihme, M. 2014 Modeling of non-equilibrium homogeneous turbulence in rapidly compressed flows. Flow Turbul. Combust. 93 (1), 93124.CrossRefGoogle Scholar
King, R. N., Hamlington, P. E. & Dahm, W. J. A. 2016 Autonomic closure for turbulence simulations. Phys. Rev. E 93 (3), 031301.Google ScholarPubMed
Laizet, S. & Lamballais, E. 2009 High-order compact schemes for incompressible flows: a simple and efficient method with quasi-spectral accuracy. J. Comput. Phys. 228 (16), 59896015.CrossRefGoogle Scholar
Laizet, S. & Li, N. 2011 Incompact3d: a powerful tool to tackle turbulence problems with up to O (105) computational cores. Intl J. Numer. Meth. Fluids 67 (11), 17351757.CrossRefGoogle Scholar
Lanczos, C. 1996 Linear Differential Operators. SIAM.CrossRefGoogle Scholar
Launder, B. E. & Sharma, B. I. 1974 Application of the energy-dissipation model of turbulence to the calculation of flow near a spinning disc. Lett. Heat Mass Transfer 1 (2), 131137.CrossRefGoogle Scholar
Lee, M. & Moser, R. D. 2015 Direct numerical simulation of turbulent channel flow up to Re = 5200. J. Fluid Mech. 774, 395415.CrossRefGoogle Scholar
Ling, J., Kurzawski, A. & Templeton, J. 2016 Reynolds averaged turbulence modelling using deep neural networks with embedded invariance. J. Fluid Mech. 807, 155166.CrossRefGoogle Scholar
Maduta, R. & Jakirlic, S. 2017 Improved RANS computations of flow over the 25° -slant-angle Ahmed body. SAE Intl J. Passenger Cars – Mech. Syst. 10 (2), 649661.CrossRefGoogle Scholar
Maulik, R. & San, O. 2017 A neural network approach for the blind deconvolution of turbulent flows. J. Fluid Mech. 831, 151181.CrossRefGoogle Scholar
Maulik, R., San, O., Rasheed, A. & Vedula, P. 2019 Sub-grid modelling for two-dimensional turbulence using neural networks. J. Fluid Mech. 858, 122144.CrossRefGoogle Scholar
Menter, F. R. 1994 Two-equation eddy-viscosity turbulence models for engineering applications. AIAA J. 32 (8), 15981605.CrossRefGoogle Scholar
Parish, E. J. & Duraisamy, K. 2016 A paradigm for data-driven predictive modeling using field inversion and machine learning. J. Comput. Phys. 305, 758774.CrossRefGoogle Scholar
Pinelli, A., Uhlmann, M., Sekimoto, A. & Kawahara, G. 2010 Reynolds number dependence of mean flow structure in square duct turbulence. J. Fluid Mech. 644, 107122.CrossRefGoogle Scholar
Pope, S. B. 1975 A more general effective-viscosity hypothesis. J. Fluid Mech. 72 (2), 331340.CrossRefGoogle Scholar
Pope, S. B. 2000 Turbulent Flows. Cambridge University Press.CrossRefGoogle Scholar
Poroseva, S. V., Colmenares, F., Juan, D. & Murman, S. M. 2016 On the accuracy of RANS simulations with DNS data. Phys. Fluids 28 (11), 115102.CrossRefGoogle ScholarPubMed
Singh, A. P. & Duraisamy, K. 2016 Using field inversion to quantify functional errors in turbulence closures. Phys. Fluids 28, 045110.CrossRefGoogle Scholar
Singh, A. P., Medida, S. & Duraisamy, K. 2017 Machine-learning-augmented predictive modeling of turbulent separated flows over airfoils. AIAA J. 55 (7), 22152227.CrossRefGoogle Scholar
Slotnick, J., Khodadoust, A., Alonso, J., Darmofal, D., Gropp, W., Lurie, E. & Mavriplis, D.2014 CFD Vision 2030 Study: a path to revolutionary computational aerosciences. Tech. Rep. National Aeronautics and Space Administration, Langley Research Center, Hampton, VA.Google Scholar
Spalart, P. R. & Allmaras, S. R. 1994 A one-equation turbulence model for aerodynamic flows. Rech. Aerosp, 1, 521.Google Scholar
Speziale, C. G. & Xu, X.-H. 1996 Towards the development of second-order closure models for nonequilibrium turbulent flows. Intl J. Heat Fluid Flow 17 (3), 238244.CrossRefGoogle Scholar
Steele, J. M. 2004 The Cauchy–Schwarz Master Class: An Introduction to the Art of Mathematical Inequalities. Cambridge University Press.CrossRefGoogle Scholar
Strang, G. 2016 Introduction to Linear Algebra, 5th edn. Wellesley–Cambridge Press.Google Scholar
Thompson, R. L., Sampaio, L. E. B., de Bragança Alves, F. A. V., Thais, L. & Mompean, G. 2016 A methodology to evaluate statistical errors in DNS data of plane channel flows. Comput. Fluids 130, 17.CrossRefGoogle Scholar
Vollant, A., Balarac, G. & Corre, C. 2017 Subgrid-scale scalar flux modelling based on optimal estimation theory and machine-learning procedures. J. Turbul. 18 (9), 854878.CrossRefGoogle Scholar
Wang, J.-X., Wu, J.-L. & Xiao, H. 2017 Physics-informed machine learning approach for reconstructing Reynolds stress modeling discrepancies based on DNS data. Phys. Rev. Fluids 2 (3), 034603.CrossRefGoogle Scholar
Weatheritt, J. & Sandberg, R. 2016 A novel evolutionary algorithm applied to algebraic modifications of the RANS stress–strain relationship. J. Comput. Phys. 325, 2237.CrossRefGoogle Scholar
Weller, H. G., Tabor, G., Jasak, H. & Fureby, C. 1998 A tensorial approach to computational continuum mechanics using object-oriented techniques. Comput. Phys. 12 (6), 620631.CrossRefGoogle Scholar
Wilcox, D. C. 1988 Reassessment of the scale-determining equation for advanced turbulence models. AIAA J. 26 (11), 12991310.CrossRefGoogle Scholar
Wu, J.-L., Xiao, H. & Paterson, E. 2018 Physics-informed machine learning approach for augmenting turbulence models: a comprehensive framework. Phys. Rev. Fluids 3 (7), 074602.CrossRefGoogle Scholar
Zhu, L., Zhang, W., Kou, J. & Liu, Y. 2019 Machine learning methods for turbulence modeling in subsonic flows around airfoils. Phys. Fluids 31 (1), 015105.CrossRefGoogle Scholar

Full text views

Full text views reflects PDF downloads, PDFs sent to Google Drive, Dropbox and Kindle and HTML full text views.

Total number of HTML views: 55
Total number of PDF views: 657 *
View data table for this chart

* Views captured on Cambridge Core between 29th April 2019 - 16th January 2021. This data will be updated every 24 hours.

Hostname: page-component-77fc7d77f9-vchrx Total loading time: 0.334 Render date: 2021-01-16T06:00:34.675Z Query parameters: { "hasAccess": "0", "openAccess": "0", "isLogged": "0", "lang": "en" } Feature Flags last update: Sat Jan 16 2021 05:52:35 GMT+0000 (Coordinated Universal Time) Feature Flags: { "metrics": true, "metricsAbstractViews": false, "peerReview": true, "crossMark": true, "comments": true, "relatedCommentaries": true, "subject": true, "clr": true, "languageSwitch": true, "figures": false, "newCiteModal": false, "shouldUseShareProductTool": true, "shouldUseHypothesis": true, "isUnsiloEnabled": true }

Send article to Kindle

To send this article to your Kindle, first ensure no-reply@cambridge.org is added to your Approved Personal Document E-mail List under your Personal Document Settings on the Manage Your Content and Devices page of your Amazon account. Then enter the ‘name’ part of your Kindle email address below. Find out more about sending to your Kindle. Find out more about sending to your Kindle.

Note you can select to send to either the @free.kindle.com or @kindle.com variations. ‘@free.kindle.com’ emails are free but can only be sent to your device when it is connected to wi-fi. ‘@kindle.com’ emails can be delivered even when you are not connected to wi-fi, but note that service fees apply.

Find out more about the Kindle Personal Document Service.

Reynolds-averaged Navier–Stokes equations with explicit data-driven Reynolds stress closure can be ill-conditioned
Available formats
×

Send article to Dropbox

To send this article to your Dropbox account, please select one or more formats and confirm that you agree to abide by our usage policies. If this is the first time you use this feature, you will be asked to authorise Cambridge Core to connect with your <service> account. Find out more about sending content to Dropbox.

Reynolds-averaged Navier–Stokes equations with explicit data-driven Reynolds stress closure can be ill-conditioned
Available formats
×

Send article to Google Drive

To send this article to your Google Drive account, please select one or more formats and confirm that you agree to abide by our usage policies. If this is the first time you use this feature, you will be asked to authorise Cambridge Core to connect with your <service> account. Find out more about sending content to Google Drive.

Reynolds-averaged Navier–Stokes equations with explicit data-driven Reynolds stress closure can be ill-conditioned
Available formats
×
×

Reply to: Submit a response


Your details


Conflicting interests

Do you have any conflicting interests? *