Ivonne Leonor Medina Lino, Mariana Carrasco-Teja, Ian Frigaard
{"title":"GPU computing of yield stress fluid flows in narrow gaps","authors":"Ivonne Leonor Medina Lino, Mariana Carrasco-Teja, Ian Frigaard","doi":"10.1007/s00162-023-00674-x","DOIUrl":null,"url":null,"abstract":"<p>We present a Graphic Processing Units (GPU) implementation of non-Newtonian Hele-Shaw flow that models the displacement of Herschel-Bulkley fluids along narrow eccentric annuli. This flow is characteristic of many long-thin flows that require extensive calculation due to an inherent nonlinearity in the constitutive law. A common method of dealing with such flows is via an augmented Lagrangian algorithm, which is often painfully slow. Here we show that such algorithms, although involving slow iterations, can often be accelerated via parallel implementation on GPUs. Indeed, such algorithms explicitly solve the nonlinear aspects only locally on each mesh cell (or node), which makes them ideal candidates for GPUs. Combined with other advances, the optimized GPU implementation takes <span>\\(\\approx 2.5\\%\\)</span> of the time of the original algorithm.</p>","PeriodicalId":795,"journal":{"name":"Theoretical and Computational Fluid Dynamics","volume":"37 5","pages":"661 - 680"},"PeriodicalIF":2.2000,"publicationDate":"2023-10-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Theoretical and Computational Fluid Dynamics","FirstCategoryId":"5","ListUrlMain":"https://link.springer.com/article/10.1007/s00162-023-00674-x","RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q2","JCRName":"MECHANICS","Score":null,"Total":0}
引用次数: 0
Abstract
We present a Graphic Processing Units (GPU) implementation of non-Newtonian Hele-Shaw flow that models the displacement of Herschel-Bulkley fluids along narrow eccentric annuli. This flow is characteristic of many long-thin flows that require extensive calculation due to an inherent nonlinearity in the constitutive law. A common method of dealing with such flows is via an augmented Lagrangian algorithm, which is often painfully slow. Here we show that such algorithms, although involving slow iterations, can often be accelerated via parallel implementation on GPUs. Indeed, such algorithms explicitly solve the nonlinear aspects only locally on each mesh cell (or node), which makes them ideal candidates for GPUs. Combined with other advances, the optimized GPU implementation takes \(\approx 2.5\%\) of the time of the original algorithm.
期刊介绍:
Theoretical and Computational Fluid Dynamics provides a forum for the cross fertilization of ideas, tools and techniques across all disciplines in which fluid flow plays a role. The focus is on aspects of fluid dynamics where theory and computation are used to provide insights and data upon which solid physical understanding is revealed. We seek research papers, invited review articles, brief communications, letters and comments addressing flow phenomena of relevance to aeronautical, geophysical, environmental, material, mechanical and life sciences. Papers of a purely algorithmic, experimental or engineering application nature, and papers without significant new physical insights, are outside the scope of this journal. For computational work, authors are responsible for ensuring that any artifacts of discretization and/or implementation are sufficiently controlled such that the numerical results unambiguously support the conclusions drawn. Where appropriate, and to the extent possible, such papers should either include or reference supporting documentation in the form of verification and validation studies.