Coupled bio-chemo-hydro-mechanical modeling of microbially induced calcite precipitation process considering biomass encapsulation using a micro-scale relationship

Geomaterials with inferior hydraulic and strength characteristics often need improvement to enhance their engineering behaviors.Traditional ground improvement techniques require enormous mechanical effort or synthetic chemicals.Sustainable stabilization technique such as microbially induced calcite precipitation(MICP)utilizes bacterial metabolic processes to precipitate cementitious calcium carbonate.The reactive transport of biochemical species in the soil mass initiates the precipitation of biocement during the MICP process.The precipitated biocement alters the hydro-mechanical performance of the soil mass.Usually,the flow,deformation,and transport phenomena regulate the biocementation technique via coupled bio-chemo-hydro-mechanical(BCHM)processes.Among all,one crucial phenomenon controlling the precipitation mechanism is the encapsulation of biomass by calcium carbonate.Biomass encapsulation can potentially reduce the biochemical reaction rate and decelerate biocementation.Laboratory examination of the encapsulation process demands a thorough analysis of associated coupled effects.Despite this,a numerical model can assist in capturing the coupled processes influencing encapsulation during the MICP treatment.However,most numerical models did not consider biochemical reaction rate kinetics accounting for the influence of bacterial encapsulation.Given this,the current study developed a coupled BCHM model to evaluate the effect of encapsulation on the precipitated calcite content using a micro-scale semiempirical relationship.Firstly,the developed BCHM model was verified and validated using numerical and experimental observations of soil column tests.Later,the encapsulation phenomenon was investigated in the soil columns of variable maximum calcite crystal sizes.The results depict altered reaction rates due to the encapsulation phenomenon and an observable change in the precipitated calcite content for each maximum crystal size.Furthermore,the permeability and deformation of the soil mass were affected by the simultaneous precipitation of calcium carbonate.Overall,the present study comprehended the influence of the encapsulation of bacteria on cement morphology-induced permeability,biocement-induced stresses and displacements.

On the analytical solution of transient friction in channel flows

The transient friction in channel mean flows is the sum of two contributions,i.e.,the underlying laminar flow(ULF)and the purely turbulent component(PTC),and the contributions are analyzed separately by theoretical experiments.It is found that,the transient friction may be higher or remarkably lower than that in equal-Reynolds number steady-state flows.The universal time constant for plane-parallel laminar flows is reported,and the role of the time constant in a turbulent mean flow is examined.It is shown that the time constant is related to the turbulence's frozen time.Finally,a study of the logarithmic layer during the transient flow is accomplished,which shows that the logarithmic layer is destroyed.

Debris Fan Produced by Failure of Canyon-Blocking Pyroclastic Flows

Ash-rich pyroclastic flows from the cataclysmic eruption of Mount Mazama (~7700 yr. B. P.), Cascade volcanic arc, Oregon, entered and blocked the narrow, bedrock-lined canyon of the Williamson River approximately 35 to 44 km from the source volcano. The blockage impounded a body of water which then released producing four stratigraphic units in the downstream debris fan. The four stratigraphic units are a boulder core comprised of locally sourced bedrock boulders and three sand-rich units including a fine-grained sand unit, a sandy pumice gravel (±basalt/hydrovolcanic tuff) unit, and a pumice pebble-bearing, crystal-rich sand unit. Hand-drilled auger holes up to ~1.6 m deep were used to obtain samples of the sand-rich units. Units were delimited using surface and down-hole observations, composition and texture, estimated density, statistical parameters of grain size, and vertical and lateral distribution of properties. Overtopping followed by rapid incision into the ash-rich pyroclastic flows progressively cleared the canyon, but a bedrock knickpoint near the head of the canyon limited the volume of debris available for transport to about 0.04 km3 to 0.08 km3. Co-deposition of bedrock boulders and lithic-rich sand was followed by rapid deposition with minimal reworking of remobilized pyroclastics. Continued draining of the impounded lake sent hyperconcentrated flows onto the debris fan depositing pumice-rich gravels that graded upward to crystal-rich sands.

Verification and Validation of High-Resolution Inviscid and Viscous Conical Nozzle Flows

Capturing elaborated flow structures and phenomena is required for well-solved numerical flows.The finite difference methods allow simple discretization of mesh and model equations.However,they need simpler meshes,e.g.,rectangular.The inverse Lax-Wendroff(ILW)procedure can handle complex geometries for rectangular meshes.High-resolution and high-order methods can capture elaborated flow structures and phenomena.They also have strong mathematical and physical backgrounds,such as positivity-preserving,jump conditions,and wave propagation concepts.We perceive an effort toward direct numerical simulation,for instance,regarding weighted essentially non-oscillatory(WENO)schemes.Thus,we propose to solve a challenging engineering application without turbulence models.We aim to verify and validate recent high-resolution and high-order methods.To check the solver accuracy,we solved vortex and Couette flows.Then,we solved inviscid and viscous nozzle flows for a conical profile.We employed the finite difference method,positivity-preserving Lax-Friedrichs splitting,high-resolution viscous terms discretization,fifth-order multi-resolution WENO,ILW,and third-order strong stability preserving Runge-Kutta.We showed the solver is high-order and captured elaborated flow structures and phenomena.One can see oblique shocks in both nozzle flows.In the viscous flow,we also captured a free-shock separation,recirculation,entrainment region,Mach disk,and the diamond-shaped pattern of nozzle flows.

Depth-Guided Vision Transformer With Normalizing Flows for Monocular 3D Object Detection

Monocular 3D object detection is challenging due to the lack of accurate depth information.Some methods estimate the pixel-wise depth maps from off-the-shelf depth estimators and then use them as an additional input to augment the RGB images.Depth-based methods attempt to convert estimated depth maps to pseudo-LiDAR and then use LiDAR-based object detectors or focus on the perspective of image and depth fusion learning.However,they demonstrate limited performance and efficiency as a result of depth inaccuracy and complex fusion mode with convolutions.Different from these approaches,our proposed depth-guided vision transformer with a normalizing flows(NF-DVT)network uses normalizing flows to build priors in depth maps to achieve more accurate depth information.Then we develop a novel Swin-Transformer-based backbone with a fusion module to process RGB image patches and depth map patches with two separate branches and fuse them using cross-attention to exchange information with each other.Furthermore,with the help of pixel-wise relative depth values in depth maps,we develop new relative position embeddings in the cross-attention mechanism to capture more accurate sequence ordering of input tokens.Our method is the first Swin-Transformer-based backbone architecture for monocular 3D object detection.The experimental results on the KITTI and the challenging Waymo Open datasets show the effectiveness of our proposed method and superior performance over previous counterparts.

A Novel Model for the Prediction of Liquid Film Thickness Distribution in Pipe Gas-Liquid Flows

A model is proposed for liquid film profile prediction in gas-liquid two-phase flow,which is able to provide the film thickness along the circumferential direction and the pressure gradient in the flow direction.A two-fluid model is used to calculate both gas and liquid phases’flow characteristics.The secondary flow occurring in the gas phase is taken into account and a sailing boat mechanism is introduced.Moreover,energy conservation is applied for obtaining the liquid film thickness distribution along the circumference.Liquid film thickness distribution is calculated accordingly for different cases;its values are compared with other models and available experimental data.As a result,the newly proposed model is tested and good performances are demonstrated.The liquid film thickness distribution in small pipes and inclined pipes is also studied,and regime transition is revealed by liquid film profile evolution.The observed inflection point demonstrates that the liquid film thickness decreases steeply along the circumference,when the circle angle ranges between 30°and 50°for gas-liquid stratified flow with small superficial velocities.

THE NONLINEAR STABILITY OF PLANE PARALLEL SHEAR FLOWS WITH RESPECT TO TILTED PERTURBATIONS

The nonlinear stability of plane parallel shear flows with respect to tilted perturbations is studied by energy methods.Tilted perturbation refers to the fact that perturbations form an angleθ∈(0,π/2)with the direction of the basic flows.By defining an energy functional,it is proven that plane parallel shear flows are unconditionally nonlinearly exponentially stable for tilted streamwise perturbation when the Reynolds number is below a certain critical value and the boundary conditions are either rigid or stress-free.In the case of stress-free boundaries,by taking advantage of the poloidal-toroidal decomposition of a solenoidal field to define energy functionals,it can be even shown that plane parallel shear flows are unconditionally nonlinearly exponentially stable for all Reynolds numbers,where the tilted perturbation can be either spanwise or streamwise.

Fixed-Time Gradient Flows for Solving Constrained Optimization: A Unified Approach

The accelerated method in solving optimization problems has always been an absorbing topic.Based on the fixedtime(FxT)stability of nonlinear dynamical systems,we provide a unified approach for designing FxT gradient flows(FxTGFs).First,a general class of nonlinear functions in designing FxTGFs is provided.A unified method for designing first-order FxTGFs is shown under Polyak-Łjasiewicz inequality assumption,a weaker condition than strong convexity.When there exist both bounded and vanishing disturbances in the gradient flow,a specific class of nonsmooth robust FxTGFs with disturbance rejection is presented.Under the strict convexity assumption,Newton-based FxTGFs is given and further extended to solve time-varying optimization.Besides,the proposed FxTGFs are further used for solving equation-constrained optimization.Moreover,an FxT proximal gradient flow with a wide range of parameters is provided for solving nonsmooth composite optimization.To show the effectiveness of various FxTGFs,the static regret analyses for several typical FxTGFs are also provided in detail.Finally,the proposed FxTGFs are applied to solve two network problems,i.e.,the network consensus problem and solving a system linear equations,respectively,from the perspective of optimization.Particularly,by choosing component-wisely sign-preserving functions,these problems can be solved in a distributed way,which extends the existing results.The accelerated convergence and robustness of the proposed FxTGFs are validated in several numerical examples stemming from practical applications.

Fiber Bundle Topology Optimization for Surface Flows

This paper presents a topology optimization approach for the surface flows on variable design domains.Via this approach,the matching between the pattern of a surface flow and the 2-manifold used to define the pattern can be optimized,where the 2-manifold is implicitly defined on another fixed 2-manifold named as the base manifold.The fiber bundle topology optimization approach is developed based on the description of the topological structure of the surface flow by using the differential geometry concept of the fiber bundle.The material distribution method is used to achieve the evolution of the pattern of the surface flow.The evolution of the implicit 2-manifold is realized via a homeomorphous map.The design variable of the pattern of the surface flow and that of the implicit 2-manifold are regularized by two sequentially implemented surface-PDE filters.The two surface-PDE filters are coupled,because they are defined on the implicit 2-manifold and base manifold,respectively.The surface Navier-Stokes equations,defined on the implicit 2-manifold,are used to describe the surface flow.The fiber bundle topology optimization problem is analyzed using the continuous adjoint method implemented on the first-order Sobolev space.Several numerical examples have been provided to demonstrate this approach,where the combination of the viscous dissipation and pressure drop is used as the design objective.

Stochastic Analysis and Modeling of Velocity Observations in Turbulent Flows

Highly turbulent water flows,often encountered near human constructions like bridge piers,spillways,and weirs,display intricate dynamics characterized by the formation of eddies and vortices.These formations,varying in sizes and lifespans,significantly influence the distribution of fluid velocities within the flow.Subsequently,the rapid velocity fluctuations in highly turbulent flows lead to elevated shear and normal stress levels.For this reason,to meticulously study these dynamics,more often than not,physical modeling is employed for studying the impact of turbulent flows on the stability and longevity of nearby structures.Despite the effectiveness of physical modeling,various monitoring challenges arise,including flow disruption,the necessity for concurrent gauging at multiple locations,and the duration of measurements.Addressing these challenges,image velocimetry emerges as an ideal method in fluid mechanics,particularly for studying turbulent flows.To account for measurement duration,a probabilistic approach utilizing a probability density function(PDF)is suggested to mitigate uncertainty in estimated average and maximum values.However,it becomes evident that deriving the PDF is not straightforward for all turbulence-induced stresses.In response,this study proposes a novel approach by combining image velocimetry with a stochastic model to provide a generic yet accurate description of flow dynamics in such applications.This integration enables an approach based on the probability of failure,facilitating a more comprehensive analysis of turbulent flows.Such an approach is essential for estimating both short-and long-term stresses on hydraulic constructions under assessment.

Exploring the Element Integration of Urban and Rural“Space of Flows”in the Context of Consumption Upgrading:A Case Study of Changsi Village in Anhui Province

The mobility and interaction between urban and rural areas are becoming more and more intensive,and their links and exchanges are increasingly closer due to constant flow of factors such as information,capital,personnel and technology.In this context,the element integration of urban and rural“space of flows”can promote the integrated development of urban and rural areas,improve the consumption environment and experience,and promote the industrial upgrading and technological progress.To realize the element integration of urban and rural“space of flows”,it is necessary to explore and innovate in infrastructure construction,information technology application,industrial cooperation and cultural exchanges.Government departments,enterprises and social organizations also need to work together to give play to their respective advantages and jointly promote the process of element integration of urban and rural“space of flows”.

Arbitrary High-Order Fully-Decoupled Numerical Schemes for Phase-Field Models of Two-Phase Incompressible Flows

Due to the coupling between the hydrodynamic equation and the phase-field equation in two-phase incompressible flows,it is desirable to develop efficient and high-order accurate numerical schemes that can decouple these two equations.One popular and efficient strategy is to add an explicit stabilizing term to the convective velocity in the phase-field equation to decouple them.The resulting schemes are only first-order accurate in time,and it seems extremely difficult to generalize the idea of stabilization to the second-order or higher version.In this paper,we employ the spectral deferred correction method to improve the temporal accuracy,based on the first-order decoupled and energy-stable scheme constructed by the stabilization idea.The novelty lies in how the decoupling and linear implicit properties are maintained to improve the efficiency.Within the framework of the spatially discretized local discontinuous Galerkin method,the resulting numerical schemes are fully decoupled,efficient,and high-order accurate in both time and space.Numerical experiments are performed to validate the high-order accuracy and efficiency of the methods for solving phase-field models of two-phase incompressible flows.

Bound-Preserving Discontinuous Galerkin Methods with Modified Patankar Time Integrations for Chemical Reacting Flows

In this paper,we develop bound-preserving discontinuous Galerkin(DG)methods for chemical reactive flows.There are several difficulties in constructing suitable numerical schemes.First of all,the density and internal energy are positive,and the mass fraction of each species is between 0 and 1.Second,due to the rapid reaction rate,the system may contain stiff sources,and the strong-stability-preserving explicit Runge-Kutta method may result in limited time-step sizes.To obtain physically relevant numerical approximations,we apply the bound-preserving technique to the DG methods.Though traditional positivity-preserving techniques can successfully yield positive density,internal energy,and mass fractions,they may not enforce the upper bound 1 of the mass fractions.To solve this problem,we need to(i)make sure the numerical fluxes in the equations of the mass fractions are consistent with that in the equation of the density;(ii)choose conservative time integrations,such that the summation of the mass fractions is preserved.With the above two conditions,the positive mass fractions have summation 1,and then,they are all between 0 and 1.For time discretization,we apply the modified Runge-Kutta/multi-step Patankar methods,which are explicit for the flux while implicit for the source.Such methods can handle stiff sources with relatively large time steps,preserve the positivity of the target variables,and keep the summation of the mass fractions to be 1.Finally,it is not straightforward to combine the bound-preserving DG methods and the Patankar time integrations.The positivity-preserving technique for DG methods requires positive numerical approximations at the cell interfaces,while Patankar methods can keep the positivity of the pre-selected point values of the target variables.To match the degree of freedom,we use polynomials on rectangular meshes for problems in two space dimensions.To evolve in time,we first read the polynomials at the Gaussian points.Then,suitable slope limiters can be applied to enforce the positivity of the solutions at those points,which can be preserved by the Patankar methods,leading to positive updated numerical cell averages.In addition,we use another slope limiter to get positive solutions used for the bound-preserving technique for the flux.Numerical examples are given to demonstrate the good performance of the proposed schemes.

Numerical Simulation of Dam-Break Flows Using Radial Basis Functions: Application to Urban Flood Inundation

Dam-break flows pose significant threats to urban areas due to their potential for causing rapid and extensive flooding. Traditional numerical methods for simulating these events struggle with complex urban landscapes. This paper presents an alternative approach using Radial Basis Functions to simulate dam-break flows and their impact on urban flood inundation. The proposed method adapts a new strategy based on Particle Swarm Optimization for variable shape parameter selection on meshfree formulation to enhance the numerical stability and convergence of the simulation. The method’s accuracy and efficiency are demonstrated through numerical experiments, including well-known partial and circular dam-break problems and an idealized city with a single building, highlighting its potential as a valuable tool for urban flood risk management.

A flexible multiscale algorithm based on an improved smoothed particle hydrodynamics method for complex viscoelastic flows

Viscoelastic flows play an important role in numerous engineering fields,and the multiscale algorithms for simulating viscoelastic flows have received significant attention in order to deepen our understanding of the nonlinear dynamic behaviors of viscoelastic fluids.However,traditional grid-based multiscale methods are confined to simple viscoelastic flows with short relaxation time,and there is a lack of uniform multiscale scheme available for coupling different solvers in the simulations of viscoelastic fluids.In this paper,a universal multiscale method coupling an improved smoothed particle hydrodynamics(SPH)and multiscale universal interface(MUI)library is presented for viscoelastic flows.The proposed multiscale method builds on an improved SPH method and leverages the MUI library to facilitate the exchange of information among different solvers in the overlapping domain.We test the capability and flexibility of the presented multiscale method to deal with complex viscoelastic flows by solving different multiscale problems of viscoelastic flows.In the first example,the simulation of a viscoelastic Poiseuille flow is carried out by two coupled improved SPH methods with different spatial resolutions.The effects of exchanging different physical quantities on the numerical results in both the upper and lower domains are also investigated as well as the absolute errors in the overlapping domain.In the second example,the complex Wannier flow with different Weissenberg numbers is further simulated by two improved SPH methods and coupling the improved SPH method and the dissipative particle dynamics(DPD)method.The numerical results show that the physical quantities for viscoelastic flows obtained by the presented multiscale method are in consistence with those obtained by a single solver in the overlapping domain.Moreover,transferring different physical quantities has an important effect on the numerical results.

An Improved Coupled Level Set and Continuous Moment-of-Fluid Method for Simulating Multiphase Flows with Phase Change

An improved algorithm for computing multiphase flows is presented in which the multimaterial Moment-of-Fluid(MOF)algorithm for multiphase flows,initially described by Li et al.(2015),is enhanced addressing existing MOF difficulties in computing solutions to problems in which surface tension forces are crucial for understanding salient flow mechanisms.The Continuous MOF(CMOF)method is motivated in this article.The CMOF reconstruction method inherently removes the"checkerboard instability"that persists when using the MOF method on surface tension driven multiphase(multimaterial)flows.The CMOF reconstruction algorithm is accelerated by coupling the CMOF method to the level set method and coupling the CMOF method to a decision tree machine learning(ML)algorithm.Multiphase flow examples are shown in the two-dimensional(2D),three-dimensional(3D)axisymmetric"RZ",and 3D coordinate systems.Examples include two material and three material multiphase flows:bubble formation,the impingement of a liquid jet on a gas bubble in a cryogenic fuel tank,freezing,and liquid lens dynamics.

Impacts of monsoon break events in the western North Pacific on the cross-equatorial flows over the Maritime Continent

本文利用ERA5逐日再分析资料,探讨了1979-2020年间西北太平洋季风中断事件对海洋性大陆越赤道气流的影响.合成结果表明,西北太平洋季风中断事件会造成高,低空越赤道气流减弱,即高层南风异常,低层北风异常,与此相关的环流异常表现为西北太平洋高层气旋,低层反气旋的斜压结构。特别的是,西北太平洋季风中断对高空越赤道气流的影响更为显著,92%的季风中断事件都导致高空越赤道气流减弱,而只有70%的事件造成低空越赤道气流减弱,这是由于低空越赤道气流同时还受到赤道中东太平洋海温异常的调控.

4D-Flow MRI在肥厚型心肌病左室流出道血流评估中的价值探索

目的 探索四维血流(four-dimensional flow,4D-Flow)磁共振成像(magnetic resonance imaging,MRI)技术在左心室腔内应用的可行性。材料与方法 本研究为前瞻性、横断面研究,纳入2022年8月至2023年1月于我院接受心脏MRI检查的21例肥厚型心肌病患者,采用3.0 T MRI扫描仪进行二维血流(tow-dimensional flow,2D-Flow)及4D-Flow成像,收集患者一周内进行的超声心动图检查结果。采用组内相关系数(inter-class correlation coefficient,ICC)、变异系数(coefficients of variation,COV)及Bland-Altman分析比较2D-Flow、4D-Flow评估左室流出道峰值流速的可重复性及一致性,并通过Pearson相关性分析探究二者与超声心动图测量结果的关系。结果 2D-Flow和4D-Flow观察者内/观察者间的ICC分别为0.999/0.999和0.995/0.992,COV分别为0.5%/0.5%和2.4%/2.6%。4D-Flow与超声心动图的测量结果呈中度相关,相关系数r值为0.574(P=0.006),但一致性较差,ICC为0.375(P=0.013)。2D-Flow与4D-Flow和超声心动图间无显著的一致性及相关性。结论 4D-Flow技术能够可视化心腔内血流模式,对左室流出道峰值流速的测量具有高度可重复性,且与超声心动图的测量结果具有显著的一致性。

针对无线控制环境的OpenFlow流表申请优化机制

在无线控制环境下,受限于控制通道的带宽和时延,OpenFlow交换机通过Packet-In消息申请流表的过程耗时较长,这个时段内交换机可能收到同一条流的大量后续报文,并分别产生与首包报文完全相同的Packet-In消息发往控制器,明显浪费网络资源。同时,针对每条流独立发送一个Packet-In报文的方式容易导致控制通道负载以小包为主、包头开销占比过大、网络资源利用率不高的问题。因此,针对无线控制条件下报文粒度Packet-In消息触发和离散化Packet-In消息发送机制存在的问题,提出了一种基于流粒度实现Packet-In消息触发和Packet-In消息聚合发送的方案,并通过仿真试验,验证了方案的效果。

基于 Moldflow 的汽车中控台框架翘曲变形分析及优化

以某汽车中控台框架为研究对象,测量试模样品发现其翘曲变形量超过了装配要求。通过Moldflow软件模拟了该塑件实际的注塑过程,翘曲变形量的模拟值与实测平均值的最大误差为5.98%,发现该塑件翘曲变形的主要因素为冷却不均和收缩不均。本文在原物料中添加质量分数为25%的玻璃纤维以及优化工艺参数后,翘曲变形量的模拟值与初始方案相比降低了86.22%。试模验证表明,优化后的翘曲变形量模拟值与实测平均值的最大误差为4.35%,证明了Moldflow软件模拟分析的准确性。试模后各检测点的最大翘曲变形量降到了1.6 mm以下,较优化之前降低了80%以上,为类似大型复杂注塑件的翘曲变形分析及优化提供了思路。