{"title":"Adaptive mesh based efficient approximations for Darcy scale precipitation–dissolution models in porous media","authors":"Shridhar Kumar, Pratibhamoy Das, Kundan Kumar","doi":"10.1002/fld.5294","DOIUrl":"10.1002/fld.5294","url":null,"abstract":"<p>In this work, we consider the Darcy scale precipitation–dissolution reactive transport 1D and 2D models in a porous medium and provide the adaptive mesh based numerical approximations for solving them efficiently. These models consist of a convection-diffusion-reaction PDE with reactions being described by an ODE having a nonlinear, discontinuous, possibly multi-valued right hand side describing precipitate concentration. The bulk concentration in the aqueous phase develops fronts and the precipitate concentration is described by a free and time-dependent moving boundary. The time adaptive moving mesh strategy, based on equidistribution principle in space and governed by a moving mesh PDE, is utilized and modified in the context of present problem for finite difference set up in 1D and finite element set up in 2D. Moreover, we use a predictor corrector based algorithm to solve the nonlinear precipitation–dissolution models. For equidistribution approach, we choose an adaptive monitor function and smooth it based on a diffusive mechanism. Numerical tests are performed to demonstrate the accuracy and efficiency of the proposed method by examples through finite difference approach for 1D and finite element approach in 2D. The moving mesh refinement accurately resolves the front location of Darcy scale precipitation–dissolution reactive transport model and reduces the computational cost in comparison to numerical simulations using a fixed mesh.</p>","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":"96 8","pages":"1415-1444"},"PeriodicalIF":1.7,"publicationDate":"2024-04-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140836393","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 extended discontinuous Galerkin shock tracking method","authors":"Jakob Vandergrift, Florian Kummer","doi":"10.1002/fld.5293","DOIUrl":"10.1002/fld.5293","url":null,"abstract":"<p>In this paper, we introduce a novel high-order shock tracking method and provide a proof of concept. Our method leverages concepts from implicit shock tracking and extended discontinuous Galerkin methods, primarily designed for solving partial differential equations featuring discontinuities. To address this challenge, we solve a constrained optimization problem aiming at accurately fitting the zero iso-contour of a level set function to the discontinuities. Additionally, we discuss various robustness measures inspired by both numerical experiments and existing literature. Finally, we showcase the capabilities of our method through a series of two-dimensional problems, progressively increasing in complexity.</p>","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":"96 8","pages":"1384-1414"},"PeriodicalIF":1.7,"publicationDate":"2024-04-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/fld.5293","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140614726","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 modified forcing approach in the Rothman–Keller method for simulations of flow phenomena at low capillary numbers","authors":"Anand Sudha, Martin Rohde","doi":"10.1002/fld.5292","DOIUrl":"10.1002/fld.5292","url":null,"abstract":"<p>The lattice-Boltzmann method (LBM) is becoming increasingly popular for simulating multi-phase flows on the microscale because of its advantages in terms of computational efficiency. Many applications of the method are restricted to relatively simple geometries. When a more complex geometry is considered—circular and inclined microchannels—some important physical phenomena may not be accurately captured, especially at low capillary numbers. A Y-Y micro-fluidic channel, widely used for a range of applications, is an example of a more complex geometry. This work aims to capture the various flow phenomena, with an emphasis on parallel flow and leakage, using the Rothman–Keller (RK) model of the LBM. To this purpose, we modify the forcing term to implement the surface tension for use at low capillary numbers. We compare the simulation results of the RK model with and without the force modification with experiments, Volume of Fluid and the phase field method and observe that the modified forcing term is an improvement over the current RK model at low capillary numbers, and it also captures parallel flow and leakage more accurately than the other simulation techniques.</p>","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":"96 8","pages":"1365-1383"},"PeriodicalIF":1.7,"publicationDate":"2024-04-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/fld.5292","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140589010","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 reduced smoothed integration scheme of the cell-based smoothed finite element method for solving fluid–structure interaction on severely distorted meshes","authors":"Tao He, Fang-Xing Lu, Xi Ma","doi":"10.1002/fld.5289","DOIUrl":"10.1002/fld.5289","url":null,"abstract":"<p>This article describes an inexpensive partitioned coupling strategy for computational fluid–structure interaction (FSI) admitting negative-Jacobian elements. The emphasis is very much on a reduced smoothed integration (RSI) scheme of the cell-based smoothed finite element method (CSFEM) using four-node quadrilateral (Q4) elements for a cost-effective solution to the Navier–Stokes (NS) equations. In the discrete fluid field, each Q4 element is considered as one single smoothing cell so as to diminish the smoothed integration loops substantially. However, the RSI scheme does not respect the stability condition of smoothed Galerkin weak-form integral in the CSFEM. To tackle this issue, a simple hourglass control is introduced to the under-integrated formulation of the NS solver. Importantly, the stabilized RSI scheme has an inbuilt advantage of its enormous tolerance towards negative-Jacobian elements. The developed technique is easy-to-implement and has been tested in various FSI examples adopting both fine and distorted meshes.</p>","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":"96 8","pages":"1337-1364"},"PeriodicalIF":1.7,"publicationDate":"2024-04-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140589005","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}
Harvesh Kumar Rana, Amit Kumar Singh, Abhay Kumar Pandey
{"title":"Therapeutic Potential of Morin Hydrate Against Rifampicin Induced Hepato and Renotoxicity in Albino Wistar Rats: Modulation of Organ Function, Oxidative Stress and Inflammatory Response.","authors":"Harvesh Kumar Rana, Amit Kumar Singh, Abhay Kumar Pandey","doi":"10.1007/s12291-023-01145-0","DOIUrl":"10.1007/s12291-023-01145-0","url":null,"abstract":"<p><p>Tuberculosis (TB) is a challenging public health issue, particularly in poor and developing countries. Rifampicin (RIF) is one of the most common first-line anti-TB drugs but it is known for its adverse effects on the hepato-renal system. The present study investigated the efficacy of morin hydrate (MH) in protecting hepato-renal damage inflicted by RIF in rats. RIF (50 mg/kg), and a combination of RIF (50 mg/kg) and MH (50 mg/kg) were administered orally for 4 weeks in rats. Silymarin (50 mg/kg) was used as a positive control. Increased levels of serological parameters such as AST, ALT, ALP, LDH, GGT, bilirubin, triglyceride, total cholesterol, urea, uric acid, creatinine, TNF-α, IFN-γ, IL-6 along with the decreased level of IL-10, total protein and albumin were used as markers of hepatic and renal injury. Oxidative damage in the tissues was measured by the increase in lipid peroxidation and decline in GSH, SOD and catalase activities. Histopathology of liver slices was used to study hepatic architecture. Four-week RIF treatment produced altered serological parameters with an increase in pro-inflammatory cytokines in serum suggesting hepatotoxicity and nephrotoxicity. The antioxidant status of the liver and kidney (increased lipid peroxidation and decline in GSH, SOD and catalase) was compromised. Cellular damage and necrosis were observed in liver slices. MH supplementation with RIF improved hepato-renal functions by restoring the serum and tissue markers towards normal values. Histological observations authenticated the results. MH supplementation also reduced the production of pro-inflammatory cytokines. Thus, the results revealed that MH provides protection against RIF-induced hepato-renal injury.</p>","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":"95 1","pages":"197-206"},"PeriodicalIF":2.1,"publicationDate":"2024-04-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC10987459/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"76652043","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}
F. I. Molina-Herrera, L. I. Quemada-Villagómez, J. L. Navarrete-Bolaños, H. Jiménez-Islas
{"title":"Comparative analysis of nondimensionalization approaches for solving the 2-D differentially heated cavity problem","authors":"F. I. Molina-Herrera, L. I. Quemada-Villagómez, J. L. Navarrete-Bolaños, H. Jiménez-Islas","doi":"10.1002/fld.5285","DOIUrl":"10.1002/fld.5285","url":null,"abstract":"<p>This work reports a numerical study on the effect of three nondimensionalization approaches that are commonly used to solve the classic problem of the 2-D differentially heated cavity. The governing equations were discretized using orthogonal collocation with Legendre polynomials, and the resulting algebraic system was solved via Newton–Raphson method with LU factorization. The simulations were performed for Rayleigh numbers between 10<sup>3</sup> and 10<sup>8</sup>, considering the Prandtl number equal to 0.71 and a geometric aspect ratio equal to 1, analyzing the convergence and the computation time on the flow lines, isotherms and the Nusselt number. The mesh size that provides independent results was 51 × 51. Approach II was the most suitable for the nondimensionalization of the differentially heated cavity problem.</p>","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":"96 7","pages":"1276-1303"},"PeriodicalIF":1.8,"publicationDate":"2024-03-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140362246","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}
Qi Jia, Jin Zhang, Wen-zhi Liang, Pei-qing Liu, Qiu-lin Qu
{"title":"A mathematical-boundary-recognition domain-decomposition Lattice Boltzmann method combined with large eddy simulation applied to airfoil aeroacoustics simulation","authors":"Qi Jia, Jin Zhang, Wen-zhi Liang, Pei-qing Liu, Qiu-lin Qu","doi":"10.1002/fld.5287","DOIUrl":"10.1002/fld.5287","url":null,"abstract":"<p>Being a direct computational aeroacoustics method, Lattice Boltzmann method (LBM) has great potential and broad application perspective in the field of numerical simulation of aerodynamic noise due to its low dispersion and low dissipation. A series of numerical algorithms and the related improvements based on the standard LBM method are proposed and developed in this paper to adapt to the airfoil noise calculation with complex grid at middle-high Reynolds number. First, a new mathematical-boundary-recognition algorithm based on Green's formula is proposed to deal with complex curved geometric models, which is validated by three-element airfoil 30P30N benchmark. Then, in order to reduce grid redundancy and improve computing efficiency, the grid refinement technique of domain decomposition model (DDM) is adopted and also improved, which is verified by calculating the flow and sound fields around 2D and 3D cylinders at Reynolds number equal to 90,000. Finally, three different LES turbulence models are combined with the standard MRT-LBM method, where different finite difference schemes are used to solve Reynolds stress tensor which is different from the traditional one. Through the direct acoustic numerical simulation of NACA0012 airfoil at Reynolds number equal to 200,000, the effects of Smagorinsky models and Wall-adapting local eddy-viscosity (WALE) model on aerodynamic noise prediction are compared and analyzed. Overall, the proposed methodology is shown to be appropriate for predicting the aerodynamic noise at low Mach number and can successfully simulate the generation and propagation of far field acoustics.</p>","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":"96 7","pages":"1250-1275"},"PeriodicalIF":1.8,"publicationDate":"2024-03-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140370669","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}
Matheus S. S. Macedo, Matheus A. Cruz, Bernardo P. Brener, Roney L. Thompson
{"title":"A data-driven turbulence modeling for the Reynolds stress tensor transport equation","authors":"Matheus S. S. Macedo, Matheus A. Cruz, Bernardo P. Brener, Roney L. Thompson","doi":"10.1002/fld.5284","DOIUrl":"10.1002/fld.5284","url":null,"abstract":"<p>The long lasting demand for better turbulence models and the still prohibitively computational cost of high-fidelity fluid dynamics simulations, like direct numerical simulations and large eddy simulations, have led to a rising interest in coupling available high-fidelity datasets and popular, yet limited, Reynolds averaged Navier–Stokes simulations through machine learning (ML) techniques. Many of the recent advances used the Reynolds stress tensor or, less frequently, the Reynolds force vector as the target for these corrections. In the present work, we considered an unexplored strategy, namely to use the modeled terms of the Reynolds stress transport equation as the target for the ML predictions, employing a neural network approach. After that, we solve the coupled set of governing equations to obtain the mean velocity field. We apply this strategy to solve the flow through a square duct. The obtained results consistently recover the secondary flow, which is not present in the baseline simulations that used the <span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mi>κ</mi>\u0000 <mo>−</mo>\u0000 <mi>ϵ</mi>\u0000 </mrow>\u0000 <annotation>$$ kappa -epsilon $$</annotation>\u0000 </semantics></math> model. The results were compared with other approaches of the literature, showing a path that can be useful in the seek of more universal models in turbulence.</p>","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":"96 7","pages":"1194-1214"},"PeriodicalIF":1.8,"publicationDate":"2024-03-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140303186","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}
Silvia Di Francesco, Sara Venturi, Jessica Padrone, Antonio Agresta
{"title":"Development of a cascaded lattice Boltzmann model for two-layer shallow water flows","authors":"Silvia Di Francesco, Sara Venturi, Jessica Padrone, Antonio Agresta","doi":"10.1002/fld.5288","DOIUrl":"10.1002/fld.5288","url":null,"abstract":"<div>\u0000 \u0000 <p>Many environmental phenomena, such as flows in rivers or in coastal region can be characterised by means of the ‘shallow approach’. A multi-layer scheme allows to extend it to density layered shallow water flows (e.g., gravity currents). Although a variety of models allowing numerical investigation of single and multi-layer shallow water flows, based on continuum and particle approaches, have been widely discussed, there are still some computational aspects that need further investigation. Focusing on the Lattice Boltzmann models (LBM), available multi-layer models generally use the standard linear collision operator (CO). In this work we adopt a multi relaxation time (MRT) cascaded collision operator to develop a two-layered liquid Lattice-Boltzmann model (CaLB-2). Specifically, the model solves the shallow water equations, taking into account two separate sets of particle distribution function (PDF), one for each layer, solved separately. Layers are connected through coupling terms, defined as external forces that model the mutual actions between the two layers. The model is validated through comparisons with experimental and numerical results from test cases available in literature. First results are very promising, highlighting a good correspondence between simulation results and literature benchmarks.</p>\u0000 </div>","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":"96 7","pages":"1230-1249"},"PeriodicalIF":1.8,"publicationDate":"2024-03-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140301794","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":"Enhanced conservative phase field method for moving contact line problems","authors":"Mingguang Shen, Ben Q. Li","doi":"10.1002/fld.5286","DOIUrl":"10.1002/fld.5286","url":null,"abstract":"<p>The traditional Allen–Cahn phase field model doesn't conserve mass and is mostly used in solidification microstructure formation. However, a recently modified Allen–Cahn phase field model has riveted the attention of the academic community. It was obtained by subtracting the curvature-driven flow term from the advective Allen–Cahn phase field model, and thus improves the boundedness of the phase field. More recently, the model has been further refined with the recovered signed distance function to compute interface normal vectors. This paper develops a three dimensional phase field model, based on the abovementioned Allen–Cahn phase field model. The model was discretized using a finite difference method on a half-staggered grid. More important, interfacial tension was expressed in a potential form. The model was tested against a number of cases and was applied to impacts in various conditions. Besides, the model was parallelized using the shared memory parallelism, OpenMP, to facilitate computation.</p>","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":"96 7","pages":"1215-1229"},"PeriodicalIF":1.8,"publicationDate":"2024-03-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140301460","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}