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.

Volumetric lattice Boltzmann method for pore-scale mass diffusionadvection process in geopolymer porous structures

Porous materials present significant advantages for absorbing radioactive isotopes in nuclear waste streams.To improve absorption efficiency in nuclear waste treatment,a thorough understanding of the diffusion-advection process within porous structures is essential for material design.In this study,we present advancements in the volumetric lattice Boltzmann method(VLBM)for modeling and simulating pore-scale diffusion-advection of radioactive isotopes within geopolymer porous structures.These structures are created using the phase field method(PFM)to precisely control pore architectures.In our VLBM approach,we introduce a concentration field of an isotope seamlessly coupled with the velocity field and solve it by the time evolution of its particle population function.To address the computational intensity inherent in the coupled lattice Boltzmann equations for velocity and concentration fields,we implement graphics processing unit(GPU)parallelization.Validation of the developed model involves examining the flow and diffusion fields in porous structures.Remarkably,good agreement is observed for both the velocity field from VLBM and multiphysics object-oriented simulation environment(MOOSE),and the concentration field from VLBM and the finite difference method(FDM).Furthermore,we investigate the effects of background flow,species diffusivity,and porosity on the diffusion-advection behavior by varying the background flow velocity,diffusion coefficient,and pore volume fraction,respectively.Notably,all three parameters exert an influence on the diffusion-advection process.Increased background flow and diffusivity markedly accelerate the process due to increased advection intensity and enhanced diffusion capability,respectively.Conversely,increasing the porosity has a less significant effect,causing a slight slowdown of the diffusion-advection process due to the expanded pore volume.This comprehensive parametric study provides valuable insights into the kinetics of isotope uptake in porous structures,facilitating the development of porous materials for nuclear waste treatment applications.

线性Poisson-Boltzmann方程的虚单元L^(2)误差估计

针对一类线性Poisson-Boltzmann方程的虚单元L^(2)误差估计进行分析.首先引入正则化的线性Poisson-Boltzmann方程,将原问题转化为非奇性Poisson-Boltzmann方程.然后给出L^(2)范数的误差估计.最后在四边形和五边形混合多边形网格上进行数值实验,数值结果验证了理论分析的正确性.

THE SMOOTHING EFFECT IN SHARP GEVREY SPACE FOR THE SPATIALLY HOMOGENEOUS NON-CUTOFF BOLTZMANN EQUATIONS WITH A HARDPOTENTIAL

In this article, we study the smoothing effect of the Cauchy problem for the spatially homogeneous non-cutoff Boltzmann equation for hard potentials. It has long been suspected that the non-cutoff Boltzmann equation enjoys similar regularity properties as to whose of the fractional heat equation. We prove that any solution with mild regularity will become smooth in Gevrey class at positive time, with a sharp Gevrey index, depending on the angular singularity. Our proof relies on the elementary L^(2) weighted estimates.

On the spreading behavior of a droplet on a circular cylinder using the lattice Boltzmann method

The study of a droplet spreading on a circular cylinder under gravity was carried out using the pseudo-potential lattice Boltzmann high-density ratios multiphase model with a non-ideal Peng–Robinson equation of state. The calculation results indicate that the motion of the droplet on the cylinder can be divided into three stages: spreading, sliding, and aggregating.The contact length and contact time of a droplet on a cylindrical surface can be affected by factors such as the wettability gradient of the cylindrical wall, the Bond number, and droplet size. Furthermore, phase diagrams showing the relationship between Bond number, cylinder wall wettability gradient, and contact time as well as maximum contact length for three different droplet sizes are given. A theoretical foundation for additional research into the heat and mass transfer process between the droplet and the cylinder can be established by comprehending the variable rules of maximum contact length and contact time.

分数阶Navier-Stokes方程的格子Boltzmann方法

对于分数阶Navier-Stokes方程的数值求解问题,首先将该方程进行离散化处理,然后构造D1Q3格子Boltzmann模型,并采用Taylor展开和Chapman-Enskog多尺度展开等技术恢复宏观方程,同时推导出该模型平衡态分布函数的表达式。最后根据一维的两个数值算例对方程进行数值模拟以及误差分析,并将得到的数值解与精确解进行比较,从而验证格子Boltzmann方法的准确性与有效性。

Investigation of Acoustomagnetoelectric Effect in Bandgap Graphene by the Boltzmann Transport Equation

We study the acoustomagnetoelectric (AME) effect in two-dimensional graphene with an energy bandgap using the semiclassical Boltzmann transport equation within the hypersound regime, (where represents the acoustic wavenumber and is the mean free path of the electron). The Boltzmann transport equation and other relevant equations were solved analytically to obtain an expression for the AME current density, consisting of longitudinal and Hall components. Our numerical results indicate that both components of the AME current densities display oscillatory behaviour. Furthermore, geometric resonances and Weiss oscillations were each defined using the relationship between the current density and Surface Acoustic Wave (SAW) frequency and the inverse of the applied magnetic field, respectively. Our results show that the AME current density of bandgap graphene, which can be controlled to suit a particular electronic device application, is smaller than that of (gapless) graphene and is therefore, more suited for nanophotonic device applications.

Exploration of the coupled lattice Boltzmann model based on a multiphase field model:A study of the solid-liquid-gas interaction mechanism in the solidification process

A multiphase field model coupled with a lattice Boltzmann(PF-LBM)model is proposed to simulate the distribution mechanism of bubbles and solutes at the solid-liquid interface,the interaction between dendrites and bubbles,and the effects of different temperatures,anisotropic strengths and tilting angles on the solidified organization of the SCN-0.24wt.%butanedinitrile alloy during the solidification process.The model adopts a multiphase field model to simulate the growth of dendrites,calculates the growth motions of dendrites based on the interfacial solute equilibrium;and adopts a lattice Boltzmann model(LBM)based on the Shan-Chen multiphase flow to simulate the growth and motions of bubbles in the liquid phase,which includes the interaction between solid-liquid-gas phases.The simulation results show that during the directional growth of columnar dendrites,bubbles first precipitate out slowly at the very bottom of the dendrites,and then rise up due to the different solid-liquid densities and pressure differences.The bubbles will interact with the dendrite in the process of flow migration,such as extrusion,overflow,fusion and disappearance.In the case of wide gaps in the dendrite channels,bubbles will fuse to form larger irregular bubbles,and in the case of dense channels,bubbles will deform due to the extrusion of dendrites.In the simulated region,as the dendrites converge and diverge,the bubbles precipitate out of the dendrites by compression and diffusion,which also causes physical phenomena such as fusion and spillage of the bubbles.These results reveal the physical mechanisms of bubble nucleation,growth and kinematic evolution during solidification and interaction with dendrite growth.

A High-Accuracy Curve Boundary Recognition Method Based on the Lattice Boltzmann Method and Immersed Moving Boundary Method

Applying numerical simulation technology to investigate fluid-solid interaction involving complex curved bound-aries is vital in aircraft design,ocean,and construction engineering.However,current methods such as Lattice Boltzmann(LBM)and the immersion boundary method based on solid ratio(IMB)have limitations in identifying custom curved boundaries.Meanwhile,IBM based on velocity correction(IBM-VC)suffers from inaccuracies and numerical instability.Therefore,this study introduces a high-accuracy curve boundary recognition method(IMB-CB),which identifies boundary nodes by moving the search box,and corrects the weighting function in LBM by calculating the solid ratio of the boundary nodes,achieving accurate recognition of custom curve boundaries.In addition,curve boundary image and dot methods are utilized to verify IMB-CB.The findings revealed that IMB-CB can accurately identify the boundary,showing an error of less than 1.8%with 500 lattices.Also,the flow in the custom curve boundary and aerodynamic characteristics of the NACA0012 airfoil are calculated and compared to IBM-VC.Results showed that IMB-CB yields lower lift and drag coefficient errors than IBM-VC,with a 1.45%drag coefficient error.In addition,the characteristic curve of IMB-CB is very stable,whereas that of IBM-VC is not.For the moving boundary problem,LBM-IMB-CB with discrete element method(DEM)is capable of accurately simulating the physical phenomena of multi-moving particle flow in complex curved pipelines.This research proposes a new curve boundary recognition method,which can significantly promote the stability and accuracy of fluid-solid interaction simulations and thus has huge applications in engineering.

Lattice Boltzmann simulation study of anode degradation in solid oxide fuel cells during the initial aging process

For present solid oxide fuel cells(SOFCs),rapid performance degradation is observed in the initial aging process,and the dis-cussion of the degradation mechanism necessitates quantitative analysis.Herein,focused ion beam-scanning electron microscopy was em-ployed to characterize and reconstruct the ceramic microstructures of SOFC anodes.The lattice Boltzmann method(LBM)simulation of multiphysical and electrochemical processes in the reconstructed models was performed.Two samples collected from industrial-size cells were characterized,including a reduced reference cell and a cell with an initial aging process.Statistical parameters of the reconstructed microstructures revealed a significant decrease in the active triple-phase boundary and Ni connectivity in the aged cell compared with the reference cell.The LBM simulation revealed that activity degradation is dominant compared with microstructural degradation during the initial aging process,and the electrochemical reactions spread to the support layer in the aged cell.The microstructural and activity de-gradations are attributed to Ni migration and coarsening.

A combined method using Lattice Boltzmann Method(LBM)and Finite Volume Method(FVM)to simulate geothermal reservoirs in Enhanced Geothermal System(EGS)

With the development of industrial activities,global warming has accelerated due to excessive emission of CO_(2).Enhanced Geothermal System(EGS)utilizes deep geothermal heat for power generation.Although porous medium theory is commonly employed to model geothermal reservoirs in EGS,Hot Dry Rock(HDR)presents a challenge as it consists of impermeable granite with zero porosity,potentially distorting the physical interpretation.To address this,the Lattice Boltzmann Method(LBM)is employed to simulate CO_(2)flow within geothermal reservoirs and the Finite Volume Method(FVM)to solve the energy conservation equation for temperature distribution.This combined method of LBM and FVM is imple-mented using MATLAB.The results showed that the Reynolds numbers(Re)of 3,000 and 8,000 lead to higher heat extraction rates from geothermal reservoirs.However,higher Re values may accelerate thermal breakthrough,posing challenges to EGS operation.Meanwhile,non-equilibrium of density in fractures becomes more pronounced during the system's life cycle,with non-Darcy's law becoming significant at Re values of 3,000 and 8,000.Density stratification due to buoyancy effects significantly impacts temperature distribution within geothermal reservoirs,with buoyancy effects at Re=100 under gravitational influence being noteworthy.Larger Re values(3,000 and 8,000)induce stronger forced convection,leading to more uniform density distribution.The addition of proppant negatively affects heat transfer performance in geothermal reservoirs,especially in single fractures.Practical engineering considerations should determine the quantity of proppant through detailed numerical simulations.

How to Predict the Temperature of the CMB Directly Using the Hubble Parameter and the Planck Scale Using the Stefan-Boltzmann Law

Based on recent progress in quantum gravity and quantum cosmology, we are also presenting a way to estimate the temperature in the cosmos, the Hubble sphere, from a relation between the Planck temperature and the Hubble scale. Our analysis predicts the Hubble sphere temperature of 2.72 K with the one standard deviation confidence interval between 2.65 K and 2.80 K, which corresponds well with the measured temperature observed from the cosmic microwave background (CMB) of about 2.72 K. This adds evidence that there is a close connection between the Planck scale, gravity, and the cosmological scales as anticipated by Eddington already in 1918.1.

Applying the Stefan-Boltzmann Law to a Cosmological Model (a Brief Note)

This brief note brings the reader up-to-date with the recent successes of the new Haug-Tatum cosmology model. In particular, the significance of recent proof that the Stefan-Boltzmann law applies to such a model is emphasized and a rationale for this is given. Remarkably, the proposed solutions of this model have incorporated all 580 supernova redshifts in the Union2 database. Therefore, one can usefully apply this thermodynamic law in the form of a continually expanding black-body universe model. To our knowledge, no other cosmological model has achieved such high-precision observational correlation.

S型异质结光催化剂ZnFe_(2)O_(4)/WO_(3)的构筑及光催化还原CO_(2)性能

通过在WO_(3)纳米片表面负载ZnFe_(2)O_(4)纳米颗粒,构建了一系列S型异质结光催化剂ZnFe_(2)O_(4)/WO_(3),并研究了其光催化CO_(2)还原性能。在没有助催化剂和牺牲剂的条件下,所制备的ZnFe_(2)O_(4)/WO_(3)复合材料可对CO_(2)与水蒸汽进行光催化反应。优化后的材料光照5 h后CO_(2)还原产物CO和CH_(4)的产量分别为7.87和4.88μmol·g^(-1)。相对于单相组分,CO和CH_(4)的产量明显提高。光催化活性的提高,归因于ZnFe_(2)O_(4)和WO_(3)异质结的形成以及光生载流子的S型电荷传输模式。

2022年四川芦山M_(S)6.1地震前应力状态研究

为研究2022年6月1日四川芦山M_(S)6.1地震的孕育和发生过程,采用CAP方法反演了2013年芦山M_(S)7.0主震及M_(S)≥5.0余震的震源机制解,并基于应力张量方差与b值时空分布特征,探讨了芦山M_(S)6.1地震的力学机制和震源区的应力状态。结果表明:2022年芦山M_(S)6.1地震震源机制表现出与2013年芦山M_(S)7.0主震和5级余震相似的逆冲型破裂特征,压应力轴方位与龙门山断裂带南段区域应力场一致。2013年芦山M_(S)7.0地震后震中及附近的应力张量方差和b值长期处于低值状态,2022年芦山M_(S)6.1地震前震中及附近出现了应力张量方差和b值的低值异常,表明芦山余震区处于较高的应力水平。分析认为:巴颜喀拉块体持续东向运动受到华南块体的阻挡,震中所在区域长期受挤压逆冲作用,从而使芦山余震区长期处于应力积累的状态,芦山M_(S)6.1地震也是在这种动力学背景下发生的。

阴道卷曲乳杆菌S层蛋白多样性研究

目的:分析比较西安市女性阴道卷曲乳杆菌S层蛋白的多样性及NCBI数据库中阴道来源卷曲乳杆菌S层蛋白的多样性。方法:选取100例妇科患者的阴道分泌物标本作为研究对象,从其中10例标本中分离得到480株卷曲乳杆菌,其中两株卷曲乳杆菌S层蛋白有明显差异,对其进行全基因组测序;并从NCBI数据库中下载人阴道源的卷曲乳杆菌全基因组数据,分析并比较西安市来源与数据库来源的卷曲乳杆菌S层蛋白的多样性。结果:S层蛋白注释基因具有菌株内和菌株间的多样性,同一套基因组中有多个S层蛋白注释基因,菌株间的相似性偏差较大,通过系统发育分析发现,中国人阴道卷曲乳杆菌的S层蛋白基因与其他人种的具有系统发育相似性。结论:阴道卷曲乳杆菌S层蛋白具有菌种间和菌株间多样性的特征,不同人种的卷曲乳杆菌具有相似性。

结合主动光源和改进YOLOv5s模型的夜间柑橘检测方法

【目的】解决夜间环境下遮挡和较小柑橘难以准确识别的问题,实现采摘机器人全天候智能化作业。【方法】提出一种结合主动光源的夜间柑橘识别方法。首先,通过分析主动光源下颜色特征不同的夜间柑橘图像,选择最佳的光源色并进行图像采集。然后,提出一种夜间柑橘检测模型BI-YOLOv5s,该模型采用双向特征金字塔网络(Bi-FPN)进行多尺度交叉连接和加权特征融合,提高对遮挡和较小果实的识别能力;引入Coordinate attention(CA)注意力机制模块,进一步加强对目标位置信息的提取;采用融入Transformer结构的C3TR模块,在减少计算量的同时更好地提取全局信息。【结果】本文提出的BI-YOLOv5s模型在测试集上的精准率、召回率、平均准确率分别为93.4%、92.2%和97.1%,相比YOLOv5s模型分别提升了3.2、1.5和2.3个百分点。在所采用的光源色环境下,模型对夜间柑橘识别的正确率为95.3%,相比白光环境下提高了10.4个百分点。【结论】本文提出的方法对夜间环境下遮挡和小目标柑橘的识别具有较高的准确性,可为夜间果蔬智能化采摘的视觉精准识别提供技术支持。

基于改进YOLOv5s的不同成熟度苹果目标检测方法

[目的]本文旨在解决在自然环境下不同成熟度苹果目标检测精度较低的问题。[方法]提出了一种改进的YOLOv5s模型SODSTR-YOLOv5s(YOLOv5s with small detection layer and omni-dimensional dynamic convolution and swin transformer block),用于不同成熟度苹果检测。首先改进YOLOv5s的多尺度目标检测层,在Prediction中构建检测160×160特征图的检测头,提高小尺寸的不同成熟度苹果的检测精度;其次在Backbone结构中融合Swin Transformer Block,加强同级成熟度的苹果纹理特征融合,弱化纹理特征分布差异带来的消极影响,提高模型泛化能力;最后将Neck结构的Conv模块替换为动态卷积模块ODConv,细化局部特征映射,实现局部苹果细粒度特征的充分提取。基于不同成熟度苹果数据集进行试验,验证改进模型的性能。[结果]改进模型SODSTR-YOLOv5s检测的精确率、召回率、平均精度均值分别为89.1%、95.5%、93.6%,高、中、低成熟度苹果平均精度均值分别为94.1%、93.1%、93.7%,平均检测时间为16 ms,参数量为7.34 M。相比于YOLOv5s模型,改进模型SODSTR-YOLOv5s精确率、召回率、平均精度均值分别提高了3.8%、5.0%、2.9%,参数量和平均检测时间分别增加了0.32 M和5 ms。[结论]改进模型SODSTR-YOLOv5s提升了在自然环境下对不同成熟度苹果的检测能力,能较好地满足实际采摘苹果的检测要求。

献血者HBVs基因变异的生物信息学分析

目的分析HBV DNA+/HBsAg-献血者的传染性指标、S区基因序列突变情况及生物信息学特征变化。方法通过PCR方法筛查出1份HBV DNA+/HBsAg-的标本,应用酶联免疫和化学发光方法检测该标本乙肝病毒5项指标,对该标本HBVs区基因片段测序并分析变异情况,应用分析软件对所测序列进行生物信息学分析。结果该标本HBV DNA+、HBsAg-、HBsAb+、HBeAg+、HBeAb-、HBcAb+;HBVs区基因序列内发生氨基酸单位点突变(P151L),空间结构发生改变。结论HBVs区基因序列空间结构改变,可能是导致检测结果出现血清学HBsAg-/HBsAb+/HBV DNA+的原因。

在片S参数计量比对结果浅析

中国电子科技集团公司第十三研究所作为主导实验室开展了在片S参数计量比对工作,对参比实验室提交的在片S参数测量结果进行了汇总分析,并用E_n值对各参比实验室测量结果进行了评价。通过在片S参数计量比对,确保了量值传递的准确、可靠,特别是对在片S参数测量不确定度的主要来源统一了认识。同时也为业内提供了在片S参数测量一致性的比较平台。