High‐order gas kinetic flux solver for viscous compressible flow simulationsJiang, Lan; Wu, Jie; Yang, Liming; Dong, Hao
doi: 10.1002/fld.5272pmid: N/A
Although the gas kinetic schemes (GKS) have emerged as one of the powerful tools for simulating compressible flows, they exhibit several shortcomings. Since the local solution of continuous Boltzmann equation with the Maxwellian distribution function is used to calculate the numerical fluxes at the cell interface, the flux expression in GKS is usually more complicated. In this paper, a high‐order simplified gas kinetic flux solver (GKFS) is presented for simulating two‐dimensional compressible flows. Circular function‐based GKFS (C‐GKFS), which simplifies the Maxwellian distribution function into the circular function, combined with an improved weighted essentially non‐oscillatory (WENO‐Z) scheme is applied to capture more details of the flow fields with fewer grids. As a result, a simple high‐order accurate C‐GKFS is obtained, which improves the computing efficiency and reduce its complexity to facilitate the practical application of engineering. A series of benchmark‐test problems are simulated and good agreement can be obtained compared with the references, which demonstrate that the high‐order C‐GKFS can achieve the desired accuracy.
A finite element model for concentration polarization and osmotic effects in a membrane channelCarro, Nicolás; Mora, David; Vellojin, Jesus
doi: 10.1002/fld.5252pmid: N/A
In this article, we study a mathematical model that represents the concentration polarization and osmosis effects in a reverse osmosis cross‐flow channel with dense membranes at some of its boundaries. The fluid is modeled using the Navier–Stokes equations and the solution‐diffusion is used to impose the momentum balance on the membrane. The scheme consist of a conforming finite element method with the velocity–pressure formulation for the Navier–Stokes equations, together with a primal scheme for the convection–diffusion equations. The Nitsche's method is used to impose the permeability condition across the membrane. Several numerical experiments are performed to show the robustness of the method. The resulting model accurately replicates the analytical models and predicts similar results to previous works. It is found that the submerged configuration has the highest permeate production, but also has the greatest pressure loss of all three configurations studied.
Multiphase lattice Boltzmann flux solver with conservative Allen‐Cahn model for modeling high‐density‐ratio flowsChen, Z.; Sun, Y. H.
doi: 10.1002/fld.5265pmid: N/A
In this paper, the Allen‐Cahn‐Multiphase lattice Boltzmann flux solver (AC‐MLBFS) is proposed as a new and effective numerical simulation method for multiphase flows with high density ratios. The MLBFS resolves the macroscopic governing equations with the finite volume method and reconstructs numerical fluxes on the cell interface from local solutions to the lattice Boltzmann equation, which combines the advantages of conventional Navier–Stokes solvers and lattice Boltzmann methods for simulating incompressible multiphase flows while alleviating their limitations. Previous MLBFS‐based multiphase solvers performed poorly in mass conservation, which might be caused by the excessive numerical diffusion in the Cahn‐Hilliard (CH) model used as the interface tracking algorithm. To resolve this problem, the present method proposes using the conservative Allen‐Cahn (AC) model as the interfacial tracking algorithm, which can ease the numerical implementation by removing high order derivative terms and alleviate mass leakage by enforcing local mass conservation in the physical model. Numerical validations will be carried out through benchmark tests at high density ratios and in extreme conditions with large Reynolds or Weber numbers. Through these examples, the accuracy and robustness as well as the mass conservation characteristics of the proposed method are demonstrated.
Unconditionally stable fully‐discrete finite element numerical scheme for active fluid modelWang, Bo; Zhang, Yuxing; Zou, Guang‐an
doi: 10.1002/fld.5260pmid: N/A
In this paper, we propose a linear, decoupled, unconditionally stable fully‐discrete finite element scheme for the active fluid model, which is derived from the gradient flow approach for an effective non‐equilibrium free energy. The developed scheme is employed by an implicit‐explicit treatment of the nonlinear terms and a second‐order Gauge–Uzawa method for the decoupling of computations for the velocity and pressure. We rigorously prove the unique solvability and unconditional stability of the proposed scheme. Several numerical tests are presented to verify the accuracy, stability, and efficiency of the proposed scheme. We also simulate the self‐organized motion under the various external body forces in 2D and 3D cases, including the motion direction of active fluid from disorder to order. Numerical results show that the scheme has a good performance in accurately capturing and handling the complex dynamics of active fluid motion.
A comparative investigation of a time‐dependent mesh method and physics‐informed neural networks to analyze the generalized Kolmogorov–Petrovsky–Piskunov equationSultan, Saad; Zhang, Zhengce
doi: 10.1002/fld.5259pmid: N/A
The Kolmogorov–Petrovsky–Piskunov (KPP) partial differential equation (PDE) is solved in this article using the moving mesh finite difference technique (MMFDM) in conjunction with physics‐informed neural networks (PINNs). We construct a time‐dependent mesh to obtain approximate solutions for the KPP problem. The temporal derivative is discretized using a backward Euler, while the spatial derivatives are discretized using a central implicit difference scheme. Depending on the error measure, several moving mesh partial differential equations (MMPDEs) are employed along the arc‐length and curvature mesh density functions (MDF). The proposed strategy has been suggested to yield remarkably precise and consistent results. To find the approximate solution, we additionally employ physics‐informed neural networks (PINNs) to compare the outcomes of the adaptive moving mesh approach. It has been observed that solutions obtained using the moving mesh method (MMM) are sufficiently accurate, and the absolute error is also much lower than the PINNs.
Assessment of implicit adaptive mesh‐free CFD modellingZhang, Tao; Barakos, George N.
doi: 10.1002/fld.5266pmid: N/A
This work presents details and assesses implicit and adaptive mesh‐free CFD modelling approaches, to alleviate laborious mesh generation in modern CFD processes. A weighted‐least‐squares‐based, mesh‐free, discretisation scheme was first derived for the compressible RANS equations, and the implicit dual‐time stepping was adopted for improved stability and convergence. A novel weight balancing concept was introduced to improve the mesh‐free modelling on highly irregular point clouds. Automatic point cloud generations based on strand and level‐set points were also discussed. A novel, polar selection approach, was also introduced to establish high‐quality point collocations. The spatial accuracy and convergence properties were validated using 2D and 3D benchmark cases. The impact of irregular point clouds and various point collocation search methods were evaluated in detail. The proposed weight balancing and the polar selection approaches were found capable of improving the mesh‐free modelling on highly irregular point clouds. The mesh‐free flexibility was then exploited for adaptive modelling. Various adaptation strategies were assessed using simulations of an isentropic vortex, combining different point refinement mechanisms and collocation search methods. The mesh‐free modelling was then successfully applied to transonic aerofoil simulations with automated point generation. A weighted pressure gradient metric prioritising high gradient regions with large point sizes was introduced to drive the adaptation. The mesh‐free adaptation was found to effectively improve the shock resolution. The results highlight the potential of mesh‐free methods in alleviating the meshing bottleneck in modern CFD.
A two‐stage reliable computational scheme for stochastic unsteady mixed convection flow of Casson nanofluidNawaz, Yasir; Arif, Muhammad Shoaib; Nazeer, Amna; Abbasi, Javeria Nawaz; Abodayeh, Kamaleldin
doi: 10.1002/fld.5264pmid: N/A
Researchers can incorporate uncertainties in computational fluid dynamics (CFD) that go beyond the inaccuracies caused by numerical discretization thanks to stochastic simulations. This study confirms the validity of current stochastic modeling tools by providing examples of stochastic simulations in conjunction with numerical solutions for incompressible flows. A numerical technique for solving deterministic and stochastic models is developed in this work. Our approach employs the Euler‐Maruyama method for stochastic modeling, representing a stochastic version of the third‐order explicit‐implicit scheme. For the deterministic model, the scheme is third‐order accurate. The consistency and stability of the constructed scheme are provided in the mean square sense. The scheme is the predictor–corrector type that is built on two time levels. Moreover, a mathematical model of the Casson nanofluid flow with variable thermal conductivity is given with the effect of the chemical reaction. The appropriate transformations are used to condense the set of partial differential equations (PDEs) down to one that is dimensionless. The scheme is applied for the deterministic and stochastic models of dimensionless flow problems. The velocity profile's deterministic and stochastic behavior are shown using contour plots. Results show that growing values of the thermal mixed convection parameter enhance the velocity profile. This article presents the progress made in stochastic computational fluid dynamics (SCFD) and highlights the energy‐related aspects of our discoveries. Our computational approach and stochastic modeling techniques provide new insights into the energy properties of Casson nanofluid flow, specifically regarding the variability of thermal conductivity and chemical processes. Our objective is to clarify the complex interaction of these factors on energy dynamics. This article presents a contemporary summary of the latest SCFD advancements. Additionally, it highlights potential directions for future research and unresolved issues that require attention from the members of the field of computational mathematics.