ترغب بنشر مسار تعليمي؟ اضغط هنا

Numerical simulation of pulmonary airway reopening by the EOS-based multiphase lattice Boltzmann method

55   0   0.0 ( 0 )
 نشر من قبل Binghai Wen
 تاريخ النشر 2021
  مجال البحث فيزياء
والبحث باللغة English




اسأل ChatGPT حول البحث

The aerosol formation is associated with the rupture of the liquid plug during the pulmonary airway reopening. The fluid dynamics of this process is difficult to predict because the rupture involved complex liquid-gas transition. Equation of state (EOS) plays a key role in the thermodynamic process of liquid-gas transition. Here, we propose an EOS-based multiphase lattice Boltzmann model, in which the nonideal force is directly evaluated by EOSs. This multiphase model is used to model the pulmonary airway reopening and study aerosol formation during exhalation. The numerical model is first validated with the simulations of Fujioka et al.(2008). and the result is in reasonable agreement with their study. Furthermore, two rupture cases with and without aerosol formation are contrasted and analyzed. It is found that the injury on the epithelium in the case with aerosol formation is essentially the same that of without aerosol formation even while the pressure drop in airway increases by about 67%. Then extensive simulations are performed to investigate the effects of pressure drop, thickness of liquid plug and film on aerosol size and the mechanical stresses. The results show that aerosol size and the mechanical stresses increase as the pressure drop enlarges and thickness of liquid plug become thicken, while aerosol size and the mechanical stresses decrease as thickness of liquid film is thicken. The present multiphase model can be extended to study the generation and transmission of bioaerosols which can carry the bioparticles of influenza or coronavirus.

قيم البحث

اقرأ أيضاً

This work presents a new multiphase SPH model that includes the shifting algorithm and a variable smoothing length formalism to simulate multi-phase flows with accuracy and proper interphase management. The implementation was performed in the DualSPH ysics code and validated for different canonical experiments, such as the single-phase and multiphase Poiseuille and Couette test cases. The method is accurate even for the multiphase case for which two phases are simulated. The shifting algorithm and the variable smoothing length formalism has been applied in the multiphase SPH model to improve the numerical results at the interphase even when it is highly deformed and non-linear effects become important. The obtained accuracy in the validation tests and the good interphase definition in the instability cases indicate an important improvement in the numerical results compared with single-phase and multiphase models where the shifting algorithm and the variable smoothing length formalism are not applied.
Many numerical studies of blood flow impose a rigid wall assumption due to the simplicity of its implementation compared to a full coupling to a solid mechanics model. In this paper, we present a localised method for incorporating the effects of elas tic walls into blood flow simulations using the lattice Boltzmann method. We demonstrate that our approach is able to more accurately capture the flow behaviour expected in an elastic walled vessel than a rigid wall model and achieves this without a loss of computational performance. We also demonstrate that our approach can capture trends in wall shear stress distribution captured by fully coupled models in personalised vascular geometries.
Simulating inhomogeneous flows with different characteristic scales in different coordinate directions using the collide-and-stream based lattice Boltzmann methods (LBM) can be accomplished efficiently using rectangular lattice grids. We develop and investigate a new rectangular central moment LBM based on non-orthogonal moment basis (referred to as RC-LBM). The equilibria to which the central moments relax under collision in this approach are obtained from matching with those corresponding to the continuous Maxwell distribution. A Chapman-Enskog analysis is performed to derive the correction terms to the second order moment equilibria involving the grid aspect ratio and velocity gradients that restores the isotropy of the viscous stress tensor and eliminates the non-Galilean invariant cubic velocity terms of the resulting hydrodynamical equations. A special case of this rectangular formulation involving the raw moments (referred to as the RNR-LBM) is also constructed. The resulting schemes represent a considerable simplification, especially for the transformation matrices and isotropy corrections, and improvement over the existing MRT-LB schemes on rectangular lattice grids that use orthogonal moment basis. Numerical validation study of both the RC-LBM and RNR-LBM for a variety of benchmark flow problems are performed that show good accuracy at various grid aspect ratios. The ability of our proposed schemes to simulate flows using relatively lower grid aspect ratios than considered in prior rectangular LB approaches is demonstrated. Furthermore, simulations reveal the superior stability characteristics of the RC-LBM over RNR-LBM in handling shear flows at lower viscosities and/or higher characteristic velocities. In addition, computational advantages of using our rectangular LB formulation in lieu of that based on the square lattice is shown.
We present a multi-scale lattice Boltzmann scheme, which adaptively refines particles velocity space. Different velocity sets, i.e., higher- and lower-order lattices, are consistently and efficiently coupled, allowing us to use the higher-order latti ce only when and where needed. This includes regions of either high Mach number or high Knudsen number. The coupling procedure of different lattices consists of either projection of the moments of the higher-order lattice onto the lower-order lattice or lifting of the lower-order lattice to the higher-order velocity space. Both lifting and projection are local operations, which enable a flexible adaptive velocity set. The proposed scheme can be formulated both in a static and an optimal, co-moving reference frame, in the spirit of the recently introduced Particles on Demand method. The multi-scale scheme is first validated through a convected athermal vortex and also studied in a jet flow setup. The performance of the proposed scheme is further investigated through the shock structure problem and a high Knudsen Couette flow, typical examples of highly non-equilibrium flows in which the order of the velocity set plays a decisive role. The results demonstrate that the proposed multi-scale scheme can operate accurately, with flexibility in terms of the underlying models and with reduced computational requirements.
A hybrid computational method coupling the lattice-Boltzmann (LB) method and a Langevin-dynamics (LD) method is developed to simulate nanoscale particle and polymer (NPP) suspensions in the presence of both thermal fluctuation and long-range many-bod y hydrodynamic interactions (HI). Brownian motion of the NPP is explicitly captured by a stochastic forcing term in the LD method. The LD method is two-way coupled to the non-fluctuating LB fluid through a discrete LB forcing source distribution to capture the long-range HI. To ensure intrinsically linear scalability with respect to the number of particles, an Eulerian-host algorithm for short-distance particle neighbor search and interaction is developed and embedded to LB-LD framework. The validity and accuracy of the LB-LD approach are demonstrated through several sample problems. The simulation results show good agreements with theory and experiment. The LB-LD approach can be favorably incorporated into complex multiscale computational frameworks for efficiently simulating multiscale, multicomponent particulate suspension systems such as complex blood suspensions.
التعليقات
جاري جلب التعليقات جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
mircosoft-partner

هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا