Ahmed Jan, Muhammad Mushtaq, Muhammad Imran Khan, Umer Farooq
{"title":"Integrated artificial intelligence and non-similar analysis for forced convection of radially magnetized ternary hybrid nanofluid of Carreau-Yasuda fluid model over a curved stretching surface","authors":"Ahmed Jan, Muhammad Mushtaq, Muhammad Imran Khan, Umer Farooq","doi":"10.1002/fld.5325","DOIUrl":"10.1002/fld.5325","url":null,"abstract":"<p>The current study investigates the boundary layer flow of Carreau-Yasuda (C-Y) ternary hybrid nanofluid model in a porous medium across curved surface stretching at linear rate under the influence of applied radial magnetic field. <span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mi>A</mi>\u0000 <msub>\u0000 <mi>l</mi>\u0000 <mn>2</mn>\u0000 </msub>\u0000 <msub>\u0000 <mi>O</mi>\u0000 <mn>3</mn>\u0000 </msub>\u0000 </mrow>\u0000 <annotation>$$ A{l}_2{O}_3 $$</annotation>\u0000 </semantics></math>, <span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mi>F</mi>\u0000 <msub>\u0000 <mi>e</mi>\u0000 <mn>3</mn>\u0000 </msub>\u0000 <msub>\u0000 <mi>O</mi>\u0000 <mn>4</mn>\u0000 </msub>\u0000 </mrow>\u0000 <annotation>$$ F{e}_3{O}_4 $$</annotation>\u0000 </semantics></math> and <span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mi>Si</mi>\u0000 <msub>\u0000 <mi>O</mi>\u0000 <mn>2</mn>\u0000 </msub>\u0000 </mrow>\u0000 <annotation>$$ Si{O}_2 $$</annotation>\u0000 </semantics></math> are nanoparticles and ethylene glycol is considered as base fluid. The effects of viscous dissipation and ohmic heating are present in the energy equation. The governing partial differential equation (PDEs) is nondimensionalized using non-similarity transformations. They can be treated as ordinary differential equations (ODEs) using local non-similarity method and solutions are obtained via bvp4c MATLAB tools. The results are evaluated by introducing computational intelligence approach utilizing the AI-based Levenberg–Marquardt scheme with a backpropagation neural network (LMS-BPNN) to investigate flow stability. The authors intend to use AI-based LMS-BPNN is to optimize the behavior of the hybrid nanofluid (HNF) flow of Carreau-Yasuda fluid across a stretching curved sheet. Initial/reference solutions are obtained through bvp4c function (an embedded MATLAB function designed to solve systems of ODEs) by systematically adjusting input parameters as demonstrated in Scenarios 1–5. There are three options to divide the numerical data: 80% for training, 10% for testing, and an additional 10% for validation. The LMS-BPNN is used for approximate solutions of Scenario 1–5. The efficiency and reliability of LMS-BPNN are validated through fitness curves based on correlation index (R), error, and regression analysis. The velocity and temperature profiles asymptotically satisfy boundary conditions of Scenario 1–5 with LMS-BPNN.</p>","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":"96 12","pages":"1864-1882"},"PeriodicalIF":1.7,"publicationDate":"2024-08-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141941392","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"An efficient direct-forcing immersed boundary method for flow around a pair of spheres","authors":"Der Chang Lo, Katherine Lee, Pao-Lan Shen","doi":"10.1002/fld.5326","DOIUrl":"10.1002/fld.5326","url":null,"abstract":"<p>The numerical study of flow around a pair of spheres and a square array of spheres is investigated by using a direct-forcing immersed boundary method. Using high resolution three-dimensional computations, we analyzed the flow around several configurations: a sphere, a pair of spheres in a tandem arrangement with center-to-center streamwise ratio L/D ranging from 1 to 6, and a square array with 9 spheres in a uniform arrangement. In the latter case, we explore the ratio of array diameter (<i>D</i><sub>G</sub>) to sphere diameter (D) at 4, 5, 6 and 7. The center-to-center streamwise and transverse pitch is the same, varied from L/D = 1.5, 2, 2.5 to 3, and they were arranged in a square periodic array to allow uniform distribution within the array. Based on the effective direct-forcing immersed boundary projection method, the fractional time marching methodology is applied for solving four field variables involving three velocities and one pressure component. The pressure Poisson equation is advanced in space by using the fast Fourier transform (FFT) and a tridiagonal matrix algorithm (TDMA), effectively solving for the diagonally dominant tridiagonal matrix equations. A direct-forcing immersed boundary method is involved to treat the interfacial terms by adding the appropriate sources as force function at the boundary, separating the phases. Geometries featuring the stationary solid obstacles in the flow are embedded in the Cartesian grid with special discretizations near the embedded boundary using a discrete Dirac delta function to ensure the accuracy of the solution in the cut cells. An important characteristic of flow over the multiple spheres is devised by comparing with the drag and lift coefficients, as well as vortex shedding.</p>","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":"96 12","pages":"1830-1863"},"PeriodicalIF":1.7,"publicationDate":"2024-07-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141866391","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"A two-phase volume of fluid approach to model rigid-perfectly plastic granular materials","authors":"W. Düsterhöft-Wriggers, S. Schubert, T. Rung","doi":"10.1002/fld.5323","DOIUrl":"10.1002/fld.5323","url":null,"abstract":"<p>Granular flow problems characterized by large deformations are widespread in various applications, including coastal and geotechnical engineering. The paper deals with the application of a rigid-perfectly plastic two-phase model extended by the Drucker–Prager yield criterion to simulate granular media with a finite volume flow solver (FV). The model refers to the combination of a Bingham fluid and an Eulerian strain measure to assess the failure region of granular dam slides. A monolithic volume-of-fluid (VoF) method is used to distinguish between the air and granular phases, both governed by the incompressible Navier–Stokes equations. The numerical framework enables modeling of large displacements and arbitrary shapes for large-scale applications. The displayed validation and verification focuses on the rigid-perfectly plastic material model for noncohesive and cohesive materials with varying angles of repose. Results indicate a good agreement of the predicted soil surface and strain results with experimental and numerical data.</p>","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":"96 12","pages":"1813-1829"},"PeriodicalIF":1.7,"publicationDate":"2024-07-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/fld.5323","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141866463","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"A spectral element discretization for quasi-static magnetohydrodynamic flows","authors":"Mattias Brynjell-Rahkola","doi":"10.1002/fld.5321","DOIUrl":"10.1002/fld.5321","url":null,"abstract":"<p>The classical staggered <span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <msub>\u0000 <mrow>\u0000 <mi>ℙ</mi>\u0000 </mrow>\u0000 <mrow>\u0000 <mi>N</mi>\u0000 </mrow>\u0000 </msub>\u0000 <mi>-</mi>\u0000 <msub>\u0000 <mrow>\u0000 <mi>ℙ</mi>\u0000 </mrow>\u0000 <mrow>\u0000 <mi>N</mi>\u0000 <mo>−</mo>\u0000 <mn>2</mn>\u0000 </mrow>\u0000 </msub>\u0000 </mrow>\u0000 <annotation>$$ {mathbb{P}}_Nhbox{-} {mathbb{P}}_{N-2} $$</annotation>\u0000 </semantics></math> spectral element method (SEM) is revisited and extended to quasi-static magnetohydrodynamic (MHD) flows. In this realm, which is valid in the limit of vanishing magnetic Reynolds number, the evaluation of the Lorentz force in the momentum equation requires the electric current density, governed by Ohm's law and a charge conservation condition derived from Ampère's law, to be determined. Once discretized with the SEM, this translates into solving one additional problem for the electric potential involving the so-called consistent Poisson operator. The method is well suited for fully three-dimensional flows in complex geometries. Changes in resolution requirements aside, consideration of the electromagnetic quantities is estimated to increase the computational cost associated with MHD by about 40% relative to hydrodynamics. The accuracy and the capabilities of the scheme is demonstrated on a set of common flows from the MHD literature. Exponential convergence with polynomial order is confirmed for the electric current density.</p>","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":"96 11","pages":"1795-1812"},"PeriodicalIF":1.7,"publicationDate":"2024-07-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/fld.5321","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141744841","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Large Courant–Friedrichs–Lewy explicit scheme for one-dimensional hyperbolic conservation laws","authors":"Vincent Guinot, Antoine Rousseau","doi":"10.1002/fld.5322","DOIUrl":"10.1002/fld.5322","url":null,"abstract":"<p>A large Courant–Friedrichs–Lewy (CFL) algorithm is presented for the explicit, finite volume solution of hyperbolic systems of conservation laws, with a focus on the shallow water equations. The Riemann problems used in the flux computation are determined using averaging kernels that extend over several computational cells. The usual CFL stability constraint is replaced with a constraint involving the kernel support size. This makes the method unconditionally stable with respect to the size of the computational cells, allowing the computational mesh to be refined locally to an arbitrary degree without altering solution stability. The practical implementation of the method is detailed for the shallow water equations with topographical source term. Computational examples report applications of the method to the linear advection, Burgers and shallow water equations. In the case of sharp bottom discontinuities, the need for improved, well-balanced discretisations of the geometric source term is acknowledged.</p>","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":"96 11","pages":"1760-1794"},"PeriodicalIF":1.7,"publicationDate":"2024-07-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141578024","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"A hybrid marching cubes based IsoAlpha method for interface reconstruction","authors":"G.S. Abhishek, Shyamprasad Karagadde","doi":"10.1002/fld.5320","DOIUrl":"10.1002/fld.5320","url":null,"abstract":"<div>\u0000 \u0000 <p>In modelling two-phase flows, accurate representation of interfaces is crucial. A class of methods for interface reconstruction are based on isosurface extraction, which involves a non-iterative, interpolation based approach. These approaches have been shown to be faster by an order of magnitude than the conventional PLIC schemes. In this work, we present a new isosurface extraction based interface reconstruction scheme based on the Marching Cubes algorithm (MC), which is commonly used in computer graphics for visualizing isosurfaces. The MC algorithm apriori lists and categorizes all possible interface configurations in a single grid cell into a Look Up Table (LUT), which makes this approach fast and robust. We also show that for certain interface configurations, the inverse problem of obtaining the isovalue from the cell volume fraction is not surjective, and a special treatment is required while handling these cases. We then demonstrate the capabilities of the method through benchmark cases for 2D and 3D static/dynamic interface reconstruction.</p>\u0000 </div>","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":"96 11","pages":"1741-1759"},"PeriodicalIF":1.7,"publicationDate":"2024-07-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141550771","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Marina G. Kontou, Xenofon S. Trompoukis, Varvara G. Asouti, Kyriakos C. Giannakoglou
{"title":"The continuous adjoint method to the \u0000 \u0000 \u0000 γ\u0000 \u0000 −\u0000 \u0000 \u0000 \u0000 R\u0000 \u0000 ˜\u0000 \u0000 \u0000 \u0000 e\u0000 \u0000 \u0000 θ\u0000 t\u0000 \u0000 \u0000 \u0000 $$ gamma -tilde{R}{e}_{theta t} $$\u0000 transition model coupled with the Spalart–Allmaras model for compressible flows","authors":"Marina G. Kontou, Xenofon S. Trompoukis, Varvara G. Asouti, Kyriakos C. Giannakoglou","doi":"10.1002/fld.5319","DOIUrl":"10.1002/fld.5319","url":null,"abstract":"<p>The continuous adjoint method for transitional flows of compressible fluids is developed and assessed, for the first time in the literature. The gradient of aerodynamic objective functions (aerodynamic forces) with respect to design variables, in problems governed by the compressible Navier–Stokes equations coupled with the Spalart–Allmaras turbulence model and the <span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mi>γ</mi>\u0000 <mspace></mspace>\u0000 <mo>−</mo>\u0000 <mspace></mspace>\u0000 <mover>\u0000 <mrow>\u0000 <mi>R</mi>\u0000 </mrow>\u0000 <mo>˜</mo>\u0000 </mover>\u0000 <msub>\u0000 <mrow>\u0000 <mi>e</mi>\u0000 </mrow>\u0000 <mrow>\u0000 <mi>θ</mi>\u0000 <mi>t</mi>\u0000 </mrow>\u0000 </msub>\u0000 </mrow>\u0000 <annotation>$$ gamma -tilde{R}{e}_{theta t} $$</annotation>\u0000 </semantics></math> transition model (in three, non-smooth and smooth, variants of it), is computed based on the continuous adjoint method. The development of the adjoint to the smooth transition model variant proved to be beneficial. The accuracy of the computed sensitivity derivatives is verified against finite differences. Programming is performed in an in-house, vertex-centered finite-volume code, efficiently running on GPUs. The proposed continuous adjoint method is used in 2D and 3D aerodynamic shape optimization problems, namely the constrained optimization of the NLF(1)–0416 isolated airfoil and that of the ONERA M6 wing. The impact of “frozen transition” (assumption according to which the adjoint to the transition model equations are not solved) or “frozen turbulence” (by additionally ignoring the adjoint to the turbulence model) are evaluated; it is shown that both lead to inaccurate sensitivities.</p>","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":"96 11","pages":"1715-1740"},"PeriodicalIF":1.7,"publicationDate":"2024-07-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/fld.5319","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141503919","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}