An inverse analysis of fluid flow through granular media using differentiable lattice Boltzmann method

This study presents a method for the inverse analysis of fluid flow problems.The focus is put on accurately determining boundary conditions and characterizing the physical properties of granular media,such as permeability,and fluid components,like viscosity.The primary aim is to deduce either constant pressure head or pressure profiles,given the known velocity field at a steady-state flow through a conduit containing obstacles,including walls,spheres,and grains.The lattice Boltzmann method(LBM)combined with automatic differentiation(AD)(AD-LBM)is employed,with the help of the GPU-capable Taichi programming language.A lightweight tape is used to generate gradients for the entire LBM simulation,enabling end-to-end backpropagation.Our AD-LBM approach accurately estimates the boundary conditions for complex flow paths in porous media,leading to observed steady-state velocity fields and deriving macro-scale permeability and fluid viscosity.The method demonstrates significant advantages in terms of prediction accuracy and computational efficiency,making it a powerful tool for solving inverse fluid flow problems in various applications.

Three-dimensional stretched flow of Jeffrey fluid with variable thermal conductivity and thermal radiation

This article addresses the three-dimensional stretched flow of the Jeffrey fluid with thermal radiation. The thermal conductivity of the fluid varies linearly with respect to temperature. Computations are performed for the velocity and temperature fields. Graphs for the velocity and temperature are plotted to examine the behaviors with different parameters. Numerical values of the local Nusselt number are presented and discussed. The present results are compared with the existing limiting solutions, showing good agreement with each other.

Three-dimensional boundary layer flow of Maxwell nanofluid:mathematical model

The present research explores the three-dimensional boundary layer flow of the Maxwell nanofluid. The flow is generated by a bidirectional stretching surface. The mathematical formulation is carried out through a boundary layer approach with the heat source/sink, the Brownian motion, and the thermophoresis effects. The newly developed boundary conditions requiring zero nanoparticle mass flux at the boundary are employed in the flow analysis for the Maxwell fluid. The governing nonlinear boundary layer equations through appropriate transformations are reduced to the coupled nonlin- ear ordinary differential system. The resulting nonlinear system is solved. Graphs are plotted to examine the effects of various interesting parameters on the non-dimensional velocities, temperature, and concentration fields. The values of the local Nusselt number are computed and examined numerically.

Unsteady three-dimensional MHD flow of the micropolar fluid over an oscillatory disk with Cattaneo-Christov double diffusion

Cattaneo-Christov heat and mass flux models are considered rather than Fourier and Fick laws due to the presence of thermal and concentration transport hyperbolic phenomena. The generalized form of the Navier-Stokes model is considered in hydromagnetic flow. Three-dimensional (3D) unsteady fluid motion is generated by the periodic oscillations of a rotating disk. Similarity transformations are used to obtain the normalized fluid flow model. The successive over relaxation (SOR) method with finite difference schemes are accomplished for the numerical solution of the obtained partial differential non-linear system. The flow features of the velocity, microrotation, temperature, and concentration fields are discussed in pictorial forms for various physical flow parameters. The couple stresses and heat and mass transfer rates for different physical quantities are explained via tabular forms. For better insight of the physical fluid model, 3D fluid phenomena and two-dimensional (2D) contours are also plotted. The results show that the micropolar fluids contain microstructure having non-symmetric stress tensor and are useful in lubrication theory. Moreover, the thermal and concentration waves in Cattaneo-Christov models have a significance role in the laser heating and enhancement in thermal conductivity.

Three-dimensional channel flow of second grade fluid in rotating frame

An analysis is performed for the hydromagnetic second grade fluid flow between two horizontal plates in a rotating system in the presence of a magnetic field. The lower sheet is considered to be a stretching sheet, and the upper sheet is a porous solid plate. By suitable transformations, the equations of conservation of mass and momentum are reduced to a system of coupled non-linear ordinary differential equations. A series of solutions to this coupled non-linear system are obtained by a powerful analytic technique, i.e., the homotopy analysis method （HAM）. The results are presented with graphs. The effects of non-dimensional parameters R, A, M2, a, and K2 on the velocity field are discussed in detail.

Newtonian heating effects in three-dimensional flow of viscoelastic fluid

A mathematical model is constructed to investigate the three-dimensional flow of a non-Newtonian fluid. An in-compressible viscoelastic fluid is used in mathematical formulation. The conjugate convective process （in which heat the transfer rate from the bounding surface with a finite capacity is proportional to the local surface temperature） in three-dimensional flow of a differential type of non-Newtonian fluid is analyzed for the first time. Series solutions for the nonlinear differential system are computed. Plots are presented for the description of emerging parameters entering into the problem. It is observed that the conjugate heating phenomenon causes an appreciable increase in the temperature at the stretching wall.

Three-dimensional flow of Oldroyd-B fluid over surface with convective boundary conditions

The present study addresses the three-dimensional flow of an Oldroyd-B fluid over a stretching surface with convective boundary conditions. The problem formulation is presented using the conservation laws of mass, momentum, and energy. The solutions to the dimensionless problems are computed. The convergence of series solutions by the homotopy analysis method （HAM） is discussed graphically and numerically. The graphs are plotted for various parameters of the temperature profile. The series solutions are verified by providing a comparison in a limiting case. The numerical values of the local Nusselt number are analyzed.

Three-dimensional flow of Powell–Eyring nanofluid with heat and mass flux boundary conditions

This article investigates the three-dimensional flow of Powell–Eyring nanofluid with thermophoresis and Brownian motion effects. The energy equation is considered in the presence of thermal radiation. The heat and mass flux conditions are taken into account. Mathematical formulation is carried out through the boundary layer approach. The governing partial differential equations are transformed into the nonlinear ordinary differential equations through suitable variables. The resulting nonlinear ordinary differential equations have been solved for the series solutions. Effects of emerging physical parameters on the temperature and nanoparticles concentration are plotted and discussed. Numerical values of local Nusselt and Sherwood numbers are computed and examined.

Analytic homotopy solution of generalized three-dimensional channel flow due to uniform stretching of the plate

In this communication a generalized three- dimensional steady flow of a viscous fluid between two infinite parallel plates is considered. The flow is generated due to uniform stretching of the lower plate in x- and y-directions. It is assumed that the upper plate is uniformly porous and is subjected to constant injection. The governing system is fully coupled and nonlinear in nature. A complete analytic solution which is uniformly valid for all values of the dimensionless parameters β Re and λ is obtained by using a purely analytic technique, namely the homotopy analysis method. Also the effects of the parameters β Re and λ on the velocity field are discussed through graphs.

STABILIZATION EFFECT OF FRICTIONS FOR TRANSONIC SHOCKS IN STEADY COMPRESSIBLE EULER FLOWS PASSING THREE-DIMENSIONAL DUCTS

Transonic shocks play a pivotal role in designation of supersonic inlets and ramjets.For the three-dimensional steady non-isentropic compressible Euler system with frictions,we constructe a family of transonic shock solutions in rectilinear ducts with square cross-sections.In this article,we are devoted to proving rigorously that a large class of these transonic shock solutions are stable,under multidimensional small perturbations of the upcoming supersonic flows and back pressures at the exits of ducts in suitable function spaces.This manifests that frictions have a stabilization effect on transonic shocks in ducts,in consideration of previous works which shown that transonic shocks in purely steady Euler flows are not stable in such ducts.Except its implications to applications,because frictions lead to a stronger coupling between the elliptic and hyperbolic parts of the three-dimensional steady subsonic Euler system,we develop the framework established in previous works to study more complex and interesting Venttsel problems of nonlocal elliptic equations.

Numerical method of the Riemann problem for two-dimensional multi-fluid flows with general equation of state

Based on the classical Roe method, we develop an interface capture method according to the general equation of state, and extend the single-fluid Roe method to the two-dimensional （2D） multi-fluid flows, as well as construct the continuous Roe matrix for the whole flow field. The interface capture equations and fluid dynamic conservative equations are coupled together and solved by using any high-resolution schemes that usually suit for the single-fluid flows. Some numerical examples are given to illustrate the solution of 1D and 2D multi-fluid Riemann problems.

THREE-DIMENSIONAL ANALYSIS OF THE FLOW IN THE CYCLONES

This paper gives the detailed mathematical expression of the flow in the spherical coordinates system. Applying the law of conservation of mass, movement theorem of steady flow, and applying the mathematical method of stream function with consideration of the axis symmetry, the three components of velocity quantum of the flow are deduced in detail. Here the overall analysis of the flow is presented in the view of the concept of whole, and the paper gives the necessary corrections of some results of Ref.[1]

Application of He’s Variational Iterative Method for Solving Thin Film Flow Problem Arising in Non-Newtonian Fluid Mechanics

In this paper, He’s variational iteration method is successfully employed to solve a nonlinear boundary value problem arising in the study of thin film flow of a third grade fluid down an inclined plane. For comparison, the same problem is solved by the Adomian decomposition method. The results show that the difference between the two solutions is negligible. The conclusion is that this technique may be considered an alternative and efficient method for finding approximate solutions of both linear and nonlinear boundary value problems. Furthermore, the variational iteration method has an advantage over the decomposition method in that it solves the nonlinear problems without using the Adomian polynomials.

A WEIGHTED PENALTY FINITE ELEMENT METHOD FOR THE ANALYSIS OF POWER-LAW FLUID FLOW PROBLEMS

In this paper, a new finite element method for the flow analysis of the viscous incompressible power-law fluid is proposed by the use of penalty-hybrid/mixed finite element formulation and by the introduction of an alternative perturbation, which is weighted by viscosity, of the continuity equation. A numerical example is presented to exhibit the efficiency of the method.

Analytical solution of a double moving boundary problem for nonlinear flows in one-dimensional semi-infinite long porous media with low permeability

Based on Huang＇s accurate tri-sectional nonlin- ear kinematic equation （1997）, a dimensionless simplified mathematical model for nonlinear flow in one-dimensional semi-infinite long porous media with low permeability is presented for the case of a constant flow rate on the inner boundary. This model contains double moving boundaries, including an internal moving boundary and an external mov- ing boundary, which are different from the classical Stefan problem in heat conduction： The velocity of the external moving boundary is proportional to the second derivative of the unknown pressure function with respect to the distance parameter on this boundary. Through a similarity transfor- mation, the nonlinear partial differential equation （PDE） sys- tem is transformed into a linear PDE system. Then an ana- lytical solution is obtained for the dimensionless simplified mathematical model. This solution can be used for strictly checking the validity of numerical methods in solving such nonlinear mathematical models for flows in low-permeable porous media for petroleum engineering applications. Finally, through plotted comparison curves from the exact an- alytical solution, the sensitive effects of three characteristic parameters are discussed. It is concluded that with a decrease in the dimensionless critical pressure gradient, the sensi- tive effects of the dimensionless variable on the dimension- less pressure distribution and dimensionless pressure gradi- ent distribution become more serious; with an increase in the dimensionless pseudo threshold pressure gradient, the sensi- tive effects of the dimensionless variable become more serious; the dimensionless threshold pressure gradient （TPG） has a great effect on the external moving boundary but has little effect on the internal moving boundary.

The Construction Method for Solving Radial Flow Problem through the Homogeneous Reservoir

On the basis of similar structure of solutions of ordinary differential equation (ODE) boundary value problem, the similar construction method was put forward by solving problems of fluid flow in porous media through the homogeneous reservoir. It is indicate that the pressure distribution of dimensionless reservoir and bottom hole in Laplace space, which take on the radial flow, also shows similar structure, and the internal relationship between the above solutions were illustrated in detail.

A GHOST FLUID BASED FRONT TRACKING METHOD FOR MULTIMEDIUM COMPRESSIBLE FLOWS

Recent years the modify ghost fluid method （MGFM） and the real ghost fluid method （RGFM） based on Riemann problem have been developed for multimedium compressible flows. According to authors, these methods have only been used with the level set technique to track the interface. In this paper, we combine the MCFM and the RGFM respectively with front tracking method, for which the fluid interfaces are explicitly tracked by connected points. The method is tested with some one-dimensional problems, and its applicability is also studied. Furthermore, in order to capture the interface more accurately, especially for strong shock impacting on interface, a shock monitor is proposed to determine the initial states of the Riemann problem. The present method is applied to various one- dimensional problems involving strong shock-interface interaction. An extension of the present method to two dimension is also introduced and preliminary results are given.

A HYPERBOLIC SYSTEM OF CONSERVATION LAWS FOR FLUID FLOWS THROUGH COMPLIANT AXISYMMETRIC VESSELS

We are concerned with the derivation and analysis of one-dimensional hyperbolic systems of conservation laws modelling fluid flows such as the blood flow through compliant axisyminetric vessels. Early models derived are nonconservative and/or nonho- mogeneous with measure source terms, which are endowed with infinitely many Riemann solutions for some Riemann data. In this paper, we derive a one-dimensional hyperbolic system that is conservative and homogeneous. Moreover, there exists a unique global Riemann solution for the Riemann problem for two vessels with arbitrarily large Riemann data, under a natural stability entropy criterion. The Riemann solutions may consist of four waves for some cases. The system can also be written as a 3 × 3 system for which strict hyperbolicity fails and the standing waves can be regarded as the contact discontinuities corresponding to the second family with zero eigenvalue.

Similarity solutions for non-Newtonian power-law fluid flow

The problem of the boundary layer flow of power law non-Newtonian fluids with a novel boundary condition is studied. The existence and uniqueness of the solutions are examined, which are found to depend on the curvature of the solutions for different values of the power law index n. It is established with the aid of the Picard-Lindelof theorem that the nonlinear boundary value problem has a unique solution in the global domain for all values of the power law index n but with certain conditions on the curva- ture of the solutions. This is done after a suitable transformation of the dependent and independent variables. For 0 〈 n ≤ 1, the solution has a positive curvature, while for n 〉 1, the solution has a negative or zero curvature on some part of the global domain. Some solutions are presented graphically to illustrate the results and the behaviors of the solutions.

可压缩两气体流动的简化神经网络模型

实用的虚拟流体方法(practical ghost fluid method, PGFM)利用Riemann问题速度解对可压缩多介质流场界面条件进行建模。基于构造的嵌入物理约束的神经网络模型预测Riemann问题速度解的方式,给出一种两气体流动的神经网络模型简化方法。首先提出完全气体状态方程下神经网络模型输入特征采样范围从无界域到有界域的转换方法,改善模型预测不同初始条件下Riemann解的泛化性能。根据该转化方法,进一步提出一种结构更加简单的神经网络优化方法,将输入维度从5个减少到3个,有效提高神经网络的训练效果。将该神经网络代理模型应用于PGFM程序框架,通过典型的一维与二维两气体流动问题进行数值验证与对比分析。结果表明,简化的网络模型与已有研究的神经网络模型相比,能取得精度相近的计算结果。而在神经网络训练效率上,简化神经网络具有明显优势。同时因为简化神经网络采样维度少,方便尝试加密采样提高拟合精度,更具备发展潜力。