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.

Study of hybrid nanofluid flow in a stationary cone-disk system with temperature-dependent fluid properties

Cone-disk systems find frequent use such as conical diffusers,medical devices,various rheometric,and viscosimetry applications.In this study,we investigate the three-dimensional flow of a water-based Ag-Mg O hybrid nanofluid in a static cone-disk system while considering temperature-dependent fluid properties.How the variable fluid properties affect the dynamics and heat transfer features is studied by Reynolds's linearized model for variable viscosity and Chiam's model for variable thermal conductivity.The single-phase nanofluid model is utilized to describe convective heat transfer in hybrid nanofluids,incorporating the experimental data.This model is developed as a coupled system of convective-diffusion equations,encompassing the conservation of momentum and the conservation of thermal energy,in conjunction with an incompressibility condition.A self-similar model is developed by the Lie-group scaling transformations,and the subsequent self-similar equations are then solved numerically.The influence of variable fluid parameters on both swirling and non-swirling flow cases is analyzed.Additionally,the Nusselt number for the disk surface is calculated.It is found that an increase in the temperature-dependent viscosity parameter enhances heat transfer characteristics in the static cone-disk system,while the thermal conductivity parameter has the opposite effect.

Film Flow of Nano-Micropolar Fluid with Dissipation Effect

The physical problem of the thin film flow of a micropolar fluid over a dynamic and inclined substrate under the influence of gravitational and thermal forces in the presence of nanoparticles is formulated.Five different types of nanoparticle samples are accounted for in this current study,namely gold Au,silver Ag,molybdenum disulfide MoS_(2),aluminum oxide Al_(2)O_(3),and silicon dioxide SiO_(2).Blood,a micropolar fluid,serves as the common base fluid.An exact closed-form solution for this problem is derived for the first time in the literature.The results are particularly validated against those for the Newtonian fluid and show excellent agreement.It was found that increasing values of the spin boundary condition and micropolarity lead to a reduction in both the thermal and momentum boundary layers.A quantitative decay in the Nusselt number for a micropolar fluid,as compared to a Newtonian one for all the tested nanoparticles,is anticipated.Gold and silver nanoparticles(i)intensify in the flow parameter as the concentration of nanoparticles increases(ii)yield a higher thermal transfer rate,whereas molybdenum disulfide,aluminum oxide,and silicon dioxide exhibit a converse attitude for both Newtonian and micropolar fluids.The reduction in film thickness for fluid comprising gold particles,as compared to the rest of the nanoparticles,is remarkable.

Flow Field Characteristics of Multi-Trophic Artificial Reef Based on Computation Fluid Dynamics

On the basis of computational fluid dynamics,the flow field characteristics of multi-trophic artificial reefs,including the flow field distribution features of a single reef under three different velocities and the effect of spacing between reefs on flow scale and the flow state,were analyzed.Results indicate upwelling,slow flow,and eddy around a single reef.Maximum velocity,height,and volume of upwelling in front of a single reef were positively correlated with inflow velocity.The length and volume of slow flow increased with the increase in inflow velocity.Eddies were present both inside and backward,and vorticity was positively correlated with inflow velocity.Space between reefs had a minor influence on the maximum velocity and height of upwelling.With the increase in space from 0.5 L to 1.5 L(L is the reef lehgth),the length of slow flow in the front and back of the combined reefs increased slightly.When the space was 2.0 L,the length of the slow flow decreased.In four different spaces,eddies were present inside and at the back of each reef.The maximum vorticity was negatively correlated with space from 0.5 L to 1.5 L,but under 2.0 L space,the maximum vorticity was close to the vorticity of a single reef under the same inflow velocity.

Numerical Solutions of the Classical and Modified Buckley-Leverett Equations Applied to Two-Phase Fluid Flow

Climate change is a reality. The burning of fossil fuels from oil, natural gas and coal is responsible for much of the pollution and the increase in the planet’s average temperature, which has raised discussions on the subject, given the emergencies related to climate. An energy transition to clean and renewable sources is necessary and urgent, but it will not be quick. In this sense, increasing the efficiency of oil extraction from existing sources is crucial, to avoid waste and the drilling of new wells. The purpose of this work was to add diffusive and dispersive terms to the Buckley-Leverett equation in order to incorporate extra phenomena in the temporal evolution between the water-oil and oil-water transitions in the pipeline. For this, the modified Buckley-Leverett equation was discretized via essentially weighted non-oscillatory schemes, coupled with a three-stage Runge-Kutta and a fourth-order centered finite difference methods. Then, computational simulations were performed and the results showed that new features emerge in the transitions, when compared to classical simulations. For instance, the dispersive term inhibits the diffusive term, adding oscillations, which indicates that the absorption of the fluid by the porous medium occurs in a non-homogeneous manner. Therefore, based on research such as this, decisions can be made regarding the replacement of the porous medium or the insertion of new components to delay the replacement.

Particulate flow modelling in a spiral separator by using the Eulerian multi-fluid VOF approach

The Euler-Euler model is less effective in capturing the free surface of flow film in the spiral separator,and thus a Eulerian multi-fluid volume of fluid(VOF)model was first proposed to describe the particulate flow in spiral separators.In order to improve the applicability of the model in the high solid concentration system,the Bagnold effect was incorporated into the modelling framework.The capability of the proposed model in terms of predicting the flow film shape in a LD9 spiral separator was evaluated via comparison with measured flow film thicknesses reported in literature.Results showed that sharp air–water and air-pulp interfaces can be obtained using the proposed model,and the shapes of the predicted flow films before and after particle addition were reasonably consistent with the observations reported in literature.Furthermore,the experimental and numerical simulation of the separation of quartz and hematite were performed in a laboratory-scale spiral separator.When the Bagnold lift force model was considered,predictions of the grade of iron and solid concentration by mass for different trough lengths were more consistent with experimental data.In the initial development stage,the quartz particles at the bottom of the flow layer were more possible to be lifted due to the Bagnold force.Thus,a better predicted vertical stratification between quartz and hematite particles was obtained,which provided favorable conditions for subsequent radial segregation.

Influences of double diffusion upon radiative flow of thin film Maxwell fluid through a stretching channel

This work explores the influence of double diffusion over thermally radiative flow of thin film hybrid nanofluid and irreversibility generation through a stretching channel.The nanoparticles of silver and alumina have mixed in the Maxwell fluid(base fluid).Magnetic field influence has been employed to channel in normal direction.Equations that are going to administer the fluid flow have been converted to dimension-free notations by using appropriate variables.Homotopy analysis method is used for the solution of the resultant equations.In this investigation it has pointed out that motion of fluid has declined with growth in magnetic effects,thin film thickness,and unsteadiness factor.Temperature of fluid has grown up with upsurge in Brownian motion,radiation factor,and thermophoresis effects,while it has declined with greater values of thermal Maxwell factor and thickness factor of the thin film.Concentration distribution has grown up with higher values of thermophoresis effects and has declined for augmentation in Brownian motion.

Venous Doppler flow patterns,venous congestion,heart disease and renal dysfunction:A complex liaison

The World Journal of Cardiology published an article written by Kuwahara et al that we take the pleasure to comment on.We focused our attention on venous congestion.In intensive care settings,it is now widely accepted that venous congestion is an important clinical feature worthy of investigation.Evaluating venous Doppler profile abnormalities at multiple sites could suggest adequate treatment and monitor its efficacy.Renal dysfunction could trigger or worsen fluid overload in heart disease,and cardio-renal syndrome is a well-characterized spectrum of disorders describing the complex interactions between heart and kidney diseases.Fluid overload and venous congestion,including renal venous hypertension,are major determinants of acute and chronic renal dysfunction arising in heart disease.Organ congestion from venous hypertension could be involved in the development of organ injury in several clinical situations,such as critical diseases,congestive heart failure,and chronic kidney disease.Ultrasonography and abnormal Doppler flow patterns diagnose clinically significant systemic venous congestion.Cardiologists and nephrologists might use this valuable,noninvasive,bedside diagnostic tool to establish fluid status and guide clinical choices.

Non-invasive assessment for intratumoural distribution of interstitial fluid flow

Interstitial fluid plays a vital role in drug delivery and tumour treatment.However,few non-invasive measurement methods are available for measuring low-velocity biological fluid flow.Therefore,this study aimed to develop a novel technology called interstitial flow velocity-MRI.The interstitial flow velocity-MRI sequence consists of a dual inversion recovery preparation and an improved stimulated echo sequence(ISTE)combined with phase-contrast MRI.A homemade flow phantom was used to assess the feasibility and sensitivity of interstitial flow velocity-MRI.In addition,xenografts of female BALB/c mouse models of 4T1 breast cancer administered losartan(40 mg/kg)or saline(n?6)were subjected to imaging on a 7.0 T scanner to assess the in vivo interstitial fluid flow velocity.The results showed a significant correlation(P<0.001)between the theoretical velocities and velocities measured using the flow phantom.Interstitial flow velocity-MRI could detect a velocity as low as 10.21±2.65 mm/s with a spatial resolution of 0.313 mm.The losartan group had a lower mean interstitial fluid velocity than the control group(85±16 vs 113±24 mm/s).In addition,compared to the saline treatment,losartan treatment reduced the proportion of collagen fibres by 10%and 12%in the Masson and Sirius red staining groups,respectively.Interstitial flow velocity-MRI has the potential to determine interstitial fluid flow velocity non-invasively and exhibits an intuitive velocity map.

Influence of the earth rotation on the interfacial wave solutions and wave-induced tangential stress in a two-layer fluid system

Interfacial waves and wave-induced tangential stress are studied for geostrophic small amplitude waves of two-layer fluid with a top free surface and a fiat bottom. The solutions were deduced from the general form of linear fluid dynamic equations of two-layer fluid under the f-plane approximation, and wave-induced tangential stress were estimated based on the solutions obtained. As expected, the solutions derived from the present work include as special cases those obtained by Sun et al. （2004. Science in China, Ser. D, 47（12）： 1147-1154） for geostrophic small amplitude surface wave solutions and wave-induced tangential stress if the density of the upper layer is much smaller than that of the lower layer. The results show that the interface and the surface will oscillate synchronously, and the influence of the earth＇s rotation both on the surface wave solutions and the interfacial wave solutions should be considered.

Averaged Dynamics of Fluids near the Oscillating Interface in a Hele-Shaw Cell

The steady flow in a Hele-Shaw cell filled with fluids with a high viscosity contrast in the presence of fluid oscillations is experimentally studied.The control of oscillatory dynamics of multiphase systems with interfaces is a challenging technological problem.We consider miscible(water and glycerol)and immiscible(water and high-viscosity silicone oil PMS-1000)fluids under subsonic oscillations perpendicular to the interface.Observations show that the interface shape depends on the amplitude and frequency of oscillations.The interface is undisturbed only in the absence of oscillations.Under small amplitudes,the interface between water and glycerol widens due to mixing.When the critical amplitude is reached,the interface becomes unstable to the fingering instability:Aqueous fingers penetrate the high-viscosity glycerol and induce intensive mixing of miscible fluids and associated decay of the instability.After the disappearance of the fingers,the interface takes a U-shape in the central part of the cell.A similar effect is observed for immiscible fluids:The oscillating interface tends to bend to the side of a high-viscosity fluid.Again,when the critical amplitude is reached,the fingering instability arises at the convex interface.This paper focuses on the causes of bending of the initially undisturbed interface between miscible or immiscible fluids.For this purpose,we measure the steady flow velocity near the interface and in the bulk of a high-viscosity fluid using Particle Image Velocimetry(PIV).

A seismic elastic moduli module for the measurements of low-frequency wave dispersion and attenuation of fluid-saturated rocks under different pressures

Knowledge about the seismic elastic modulus dispersion,and associated attenuation,in fluid-saturated rocks is essential for better interpretation of seismic observations taken as part of hydrocarbon identification and time-lapse seismic surveillance of both conventional and unconventional reservoir and overburden performances.A Seismic Elastic Moduli Module has been developed,based on the forced-oscillations method,to experimentally investigate the frequency dependence of Young's modulus and Poisson's ratio,as well as the inferred attenuation,of cylindrical samples under different confining pressure conditions.Calibration with three standard samples showed that the measured elastic moduli were consistent with the published data,indicating that the new apparatus can operate reliably over a wide frequency range of f∈[1-2000,10^(6)]Hz.The Young's modulus and Poisson's ratio of the shale and the tight sandstone samples were measured under axial stress oscillations to assess the frequency-and pressure-dependent effects.Under dry condition,both samples appear to be nearly frequency independent,with weak pressure dependence for the shale and significant pressure dependence for the sandstone.In particular,it was found that the tight sandstone with complex pore microstructure exhibited apparent dispersion and attenuation under brine or glycerin saturation conditions,the levels of which were strongly influenced by the increased effective pressure.In addition,the measured Young's moduli results were compared with the theoretical predictions from a scaled poroelastic model with a reasonably good agreement,revealing that the combined fluid flow mechanisms at both mesoscopic and microscopic scales possibly responsible for the measured dispersion.

Towards implementation of alloy-specific thermo-fluid modelling for laser powder-bed fusion of Mg alloys

Multi-physics thermo-fluid modeling has been extensively used as an approach to understand melt pool dynamics and defect formation as well as optimizing the process-related parameters of laser powder-bed fusion(L-PBF).However,its capabilities for being implemented as a reliable tool for material design,where minor changes in material-related parameters must be accurately captured,is still in question.In the present research,first,a thermo-fluid computational fluid dynamics(CFD)model is developed and validated against experimental data.Considering the predicted material properties of the pure Mg and commercial ZK60 and WE43 Mg alloys,parametric studies are done attempting to elucidate how the difference in some of the material properties,i.e.,saturated vapor pressure,viscosity,and solidification range,can influence the melt pool dynamics.It is found that a higher saturated vapor pressure,associated with the ZK60 alloy,leads to a deeper unstable keyhole,increasing the keyhole-induced porosity and evaporation mass loss.Higher viscosity and wider solidification range can increase the non-uniformity of temperature and velocity distribution on the keyhole walls,resulting in increased keyhole instability and formation of defects.Finally,the WE43 alloy showed the best behavior in terms of defect formation and evaporation mass loss,providing theoretical support to the extensive use of this alloy in L-PBF.In summary,this study suggests an approach to investigate the effect of materials-related parameters on L-PBF melting and solidification,which can be extremely helpful for future design of new alloys suitable for L-PBF.

Pressure transient characteristics of non-uniform conductivity fractured wells in viscoelasticity polymer flooding based on oil-water two-phase flow

Polymer flooding in fractured wells has been extensively applied in oilfields to enhance oil recovery.In contrast to water,polymer solution exhibits non-Newtonian and nonlinear behavior such as effects of shear thinning and shear thickening,polymer convection,diffusion,adsorption retention,inaccessible pore volume and reduced effective permeability.Meanwhile,the flux density and fracture conductivity along the hydraulic fracture are generally non-uniform due to the effects of pressure distribution,formation damage,and proppant breakage.In this paper,we present an oil-water two-phase flow model that captures these complex non-Newtonian and nonlinear behavior,and non-uniform fracture characteristics in fractured polymer flooding.The hydraulic fracture is firstly divided into two parts:high-conductivity fracture near the wellbore and low-conductivity fracture in the far-wellbore section.A hybrid grid system,including perpendicular bisection(PEBI)and Cartesian grid,is applied to discrete the partial differential flow equations,and the local grid refinement method is applied in the near-wellbore region to accurately calculate the pressure distribution and shear rate of polymer solution.The combination of polymer behavior characterizations and numerical flow simulations are applied,resulting in the calculation for the distribution of water saturation,polymer concentration and reservoir pressure.Compared with the polymer flooding well with uniform fracture conductivity,this non-uniform fracture conductivity model exhibits the larger pressure difference,and the shorter bilinear flow period due to the decrease of fracture flow ability in the far-wellbore section.The field case of the fall-off test demonstrates that the proposed method characterizes fracture characteristics more accurately,and yields fracture half-lengths that better match engineering reality,enabling a quantitative segmented characterization of the near-wellbore section with high fracture conductivity and the far-wellbore section with low fracture conductivity.The novelty of this paper is the analysis of pressure performances caused by the fracture dynamics and polymer rheology,as well as an analysis method that derives formation and fracture parameters based on the pressure and its derivative curves.

Numerical Study on the Aerodynamic and Fluid−Structure Interaction of An NREL-5MW Wind Turbine

A 5-MW wind turbine has been modeled and analyzed for fluid-structure interaction and aerodynamic performance.In this study, a full-scale model of a 5-MW wind turbine is first developed based on a computational fluid dynamics(CFD) approach, in which the unsteady, noncompressible Reynolds Averaged Navier-Stokes(RANS) method is used. The main focus of the study is to analyze the tower shadow effect on the aerodynamic performance of the wind turbine under different inlet flow conditions. Subsequently, the finite element model is established by considering fluid/structure interactions to study the structural stress, displacement, strain distributions and flow field information of the structure under the uniform wind speed. Finally, the fluid-structure interaction model is established by considering turbulent wind and the tower shadow effect. The variation rules of the dynamic response of the one-way and two-way fluid-structure interaction(FSI) models under different wind speeds are analyzed, and the numerical calculation results are compared with those of the centralized mass model. The results show that the tower shadow effect and structural deformation are the main factors affecting the aerodynamic load fluctuation of the wind turbine, which in turn affects the aerodynamic performance and structural stability of the blades. The structural dynamic response of the coupled model shows significant similarity, while the structural displacement response of the former exhibits less fluctuation compared with the conventional centralized mass model. The one-way fluid-structure interaction(FSI)model shows a higher frequency of stress-strain and displacement oscillations on the blade compared with the two-way FSI model.

Flow Rate Measurement of Gravity Infusion Set and Functional Evaluation of Drop Counter: A Pilot Study

Developing a novel drop counter by introducing the Internet of Things concept has been vigorously conducted in recent years. Understanding the newly introduced drop counter’s flow rate control accuracy and flow rate count feature is essential for improving safety in infusion management. This study aimed to verify if the new drop counters could secure accurate flow rate and drip count by conducting actual flow rate measurements using gravimetry and functional evaluation. A drop counter was attached to each drip chamber of the infusion set, and an IV drip was conducted at the 100 ml/h flow rate. The weight of discharged physiological saline was measured to plot trumpet curves. Next, three different types of drop counters were evaluated to determine if they maintained drip count accuracy according to the changes in their position angles. The flow rate errors in all conditions indicated trumpet-like curves, exhibiting an overall error range within ±10% in all observation windows. Although every drop counter successfully detected and measured dripping, it was challenging in some counters to detect dripping when the drip chamber was tilted. In comparing adult and pediatric IV sets, the adult IV set was found to be less likely to detect dripping in the angled position. No significant differences in results were confirmed between high and low flow rates, suggesting that the drop count function would not be affected by the flow rate in the ranges of typical infusion practices. Doppler sensors have a wide range of measurements and high sensitivity;the dripping was detected successfully even when the drip chamber was tilted, probably due to the advantages of these sensors. In contrast, miscounts occurred in those equipped with infrared sensors, which could not detect light intensity changes in tilted positions. Understanding the tendencies in flow rate errors in infusion can be valuable information for infusion management.

基于Flow-Simulation的气液分离器设计仿真分析研究

为解决燃气发电机组瓦斯气源含水量高的问题,对燃气发电机组气液分离器进行了设计与仿真分析。首先根据功能要求进行气液分离器设计参数计算,用Solidworks三维建模,然后利用Flow-Simulation进行CFD流体动力学仿真,得到气液分离器内的压力云图、速度云图和运动轨迹图;从压力云图中看出压力损失为784 Pa,能够满足压力损失≤1 kPa的要求。后经现场试验气液分离器分离效率达到93.6%,有效地提高了燃气发电机组的发电效率。该款设计优化后的旋风式气液分离器能够满足设计使用要求。

Three-dimensional simulation of pore scale fluid flow in granular ore media with realistic geometry

The images of granular ore media were captured by X-ray CT scanner. Combined with digital image processing and finite element techniques, the three-dimensional geometrical model, which represents the realistic pore structure of the media, was constructed. With this model, three dimensional pore scale fluid flow among particles was simulated. Then the distributions of fluid flow velocity and pressure were analyzed and the hydraulic conductivity was calculated. The simulation results indicate the fluid flow behaviors are mainly dominated by the volume and topological structure of pore space. There exist obvious preferential flow and leaching blind zones simultaneously in the medium. The highest velocities generally occur in those narrow pores with high pressure drops. The hydraulic conductivity obtained by simulation is the same order of magnitude as the laboratory test result, which denotes the validity of the model. The pore-scale and macro-scale are combined and the established geometrical model can be used for the simulations of other phenomena during heap leaching process.

3D topographic correction of the BSR heat flow and detection of focused fluid flow

The bottom-simulating reflector（BSR） is a seismic indicator of the bottom of a gas hydrate stability zone. Its depth can be used to calculate the seafloor surface heat flow. The calculated BSR heat flow variations include disturbances from two important factors：（1） seafloor topography, which focuses the heat flow over regions of concave topography and defocuses it over regions of convex topography, and（2） the focused warm fluid flow within the accretionary prism coming from depths deeper than BSR. The focused fluid flow can be detected if the contribution of the topography to the BSR heat flow is removed. However, the analytical equation cannot solve the topographic effect at complex seafloor regions. We prove that 3D finite element method can model the topographic effect on the regional background heat flow with high accuracy, which can then be used to correct the topographic effect and obtain the BSR heat flow under the condition of perfectly flat topography. By comparing the corrected BSR heat flow with the regional background heat flow, focused fluid flow regions can be detected that are originally too small and cannot be detected using present-day equipment. This method was successfully applied to the midslope region of northern Cascadia subducting margin. The results suggest that the Cucumber Ridge and its neighboring area are positive heat flow anomalies, about 10%–20% higher than the background heat flow after 3D topographic correction. Moreover, the seismic imaging associated the positive heat flow anomaly areas with seabed fracture–cavity systems. This suggests flow of warm gas-carrying fluids along these high-permeability pathways, which could result in higher gas hydrate concentrations.

Relationship between the transition frequency of local fluid flow and the peak frequency of attenuation

Local fluid flow（LFF） at the mesoscopic scale is the main dissipation mechanism of seismic waves in heterogeneous porous media within the seismic frequency band.LFF is easily influenced by the structure and boundary conditions of the porous media,which leads to different behaviors of the peak frequency of attenuation.The associated transition frequency can provide detailed information about the trend of LFF;therefore,research on the transition frequency of LFF and its relationship with the peak frequency of the corresponding attenuation（i.e.,inverse of quality factor） facilitates the detailed understanding of the effect of inner structures and boundary conditions in porous media.In this study,we firstly obtain the transition frequency of fluid flux based on Biot＇s theory of poroelasticity and the fast Fourier transform algorithm in a sample containing one repeating unit cell（RUC）.We then analyze changes of these two frequencies in porous media with different porous properties.Finally,we extend our analysis to the influence of the undrained boundary condition on the transition frequency and peak frequency in porous media with multiple RUCs.This setup can facilitate the understanding of the effect from the undrained boundary condition.Results demonstrate that these two frequencies have the same trend at low water saturation,but amplitude variations differ between the frequencies as the amount of saturation increases.However,for cases of high water saturation,both the trend and the amplitude variation of these two frequencies fit well with each other.