Analysis of a new multispecies tumor growth model coupling 3D phase-fields with a 1D vascular network

In this work, we present and analyze a mathematical model for tumor growth incorporating ECM erosion, interstitial flow, and the effect of vascular flow and nutrient transport. The model is of phase-field or diffused-interface type in which multiple phases of cell species and other constituents are separated by smooth evolving interfaces. The model involves a mesoscale version of Darcy’s law to capture the flow mechanism in the tissue matrix. Modeling flow and transport processes in the vasculature supplying the healthy and cancerous tissue, one-dimensional (1D) equations are considered. Since the models governing the transport and flow processes are defined together with cell species models on a three-dimensional (3D) domain, we obtain a 3D–1D coupled model.     

Imposing accurate wall boundary conditions in corrective‐matrix‐based moving particle semi‐implicit method for free surface flow

Corrective matrix that is derived to restore consistency of discretization schemes can significantly enhance accuracy for the inside particles in the Moving Particle Semi‐implicit method. In this situation, the error due to free surface and wall boundaries becomes dominant. Based on the recent study on Neumann boundary condition (Matsunaga et al, CMAME, 2020), the corrective matrix schemes in MPS are generalized to straightforwardly and accurately impose Neumann boundary condition. However, the new schemes can still easily trigger instability at free surface because of the biased error caused by the incomplete/biased neighbor support. Therefore, the existing stable schemes based on virtual particles and conservative gradient models are applied to free surface and nearby particles to produce a stable transitional layer at free surface. The new corrective matrix schemes are only applied to the particles under the stable transitional layer for improving the wall boundary conditions. Three numerical examples of free surface flows demonstrate that the proposed method can help to reduce the pressure/velocity fluctuations and hence enhance accuracy further.     

石墨基液流电池复合双极板

液流电池是一种安全性高、使用寿命长、可扩展的大规模储能系统,可以协助电网调峰储能,提高能源利用率,发展前景广阔。双极板是液流电池的重要组成部分。功能上起到了分隔、串联电池、传导电流、为电堆提供结构支撑等作用。从成本构成角度看,双极板的价格占电堆成本的比重也较大。开发高性能、低成本的双极板对加快液流电池的商业化应用具有重要意义,也是目前业界的迫切需求。虽然文献上报道了许多针对液流电池双极板开发的工作,但是目前高性能、低成本的液流电池双极板产品仍无法充分满足市场需求。本文着重介绍了石墨基复合双极板的研究现状,介绍了材料选择、工艺流程对关键性能的影响,对相关工作进行了评述,并为液流电池双极板的开发提出了建议。     

Hybrid continuum–coarse-grained modeling of erythrocytes

The red blood cell (RBC) membrane is a composite structure, consisting of a phospholipid bilayer and an underlying membrane-associated cytoskeleton. Both continuum and particle-based coarse-grained RBC models make use of a set of vertices connected by edges to represent the RBC membrane, which can be seen as a triangular surface mesh for the former and a spring network for the latter. Here, we present a modeling approach combining an existing continuum vesicle model with a coarse-grained model for the cytoskeleton. Compared to other two-component approaches, our method relies on only one mesh, representing the cytoskeleton, whose velocity in the tangential direction of the membrane may be different from that of the lipid bilayer. The finitely extensible nonlinear elastic (FENE) spring force law in combination with a repulsive force defined as a power function (POW), called FENE–POW, is used to describe the elastic properties of the RBC membrane. The mechanical interaction between the lipid bilayer and the cytoskeleton is explicitly computed and incorporated into the vesicle model. Our model includes the fundamental mechanical properties of the RBC membrane, namely fluidity and bending rigidity of the lipid bilayer, and shear elasticity of the cytoskeleton while maintaining surface-area and volume conservation constraint. We present three simulation examples to demonstrate the effectiveness of this hybrid continuum–coarse-grained model for the study of RBCs in fluid flows.     

Pullback exponential attractors for the non-autonomous micropolar fluid flows

This paper investigates the pullback asymptotic behaviors for the non-autonomous micropolar fluid flows in 2D bounded domains. We use the energy method, combining with some important properties of the generated processes, to prove the existence of pullback exponential attractors and global pullback attractors and show that they both with finite fractal dimension. Further, we give the relationship between global pullback attractors and pullback exponential attractors.     

Particle manipulation via integration of electroosmotic flow of power-law fluids with standing surface acoustic waves (SSAW)

Standing surface acoustic wave (SSAW) based microfluidic devices have shown great promise toward fluid and particle manipulation applications in medicine, chemistry, and biotechnology. In this article, we present an analytical model for investigating continuous manipulation of particles (both synthetic and biological) within electroosmotic flow of non-Newtonian bio-fluids in a microfluidic channel under the influence of standing surface acoustic waves (SSAW). The particles are injected along the center of channel into the electroosmotically driven flow of power-law fluids, wherein their transport through the SSAW region is dictated by the hydrodynamic, electrophoretic, and acoustic forces. We first present a mathematical model to analyze the characteristics of electroosmotic flow of non-Newtonian power-law fluids in a hydrophobic slit microchannel. Next, we investigate the trajectories of particles in the flow field due to the combined effect of electroosmotic, electrophoretic, and acoustophoretic forcing mechanisms. The effect of key parameters such as particle size, their physical properties, input power, flow rate, and flow behavior index on the particle trajectories is examined while including the effect of the channel walls. The presented model delineates the methodologies of improving SSAW-based particle separation technology by considering the fluid rheology as well as the surface properties of the channel walls. Therefore, we believe that this model can serve as an efficient tool for device design and quick optimizations to explore novel applications concerning the integration of electroosmotic flows with acoustofluidic technologies.     

Interface‐fitted moving mesh method for axisymmetric two‐phase flow in microchannels

A boundary‐fitted moving mesh scheme is presented for the simulation of two‐phase flow in two‐dimensional and axisymmetric geometries. The incompressible Navier‐Stokes equations are solved using the finite element method, and the mini element is used to satisfy the inf‐sup condition. The interface between the phases is represented explicitly by an interface adapted mesh, thus allowing a sharp transition of the fluid properties. Surface tension is modelled as a volume force and is discretized in a consistent manner, thus allowing to obtain exact equilibrium (up to rounding errors) with the pressure gradient. This is demonstrated for a spherical droplet moving in a constant flow field. The curvature of the interface, required for the surface tension term, is efficiently computed with simple but very accurate geometric formulas. An adaptive moving mesh technique, where smoothing mesh velocities and remeshing are used to preserve the mesh quality, is developed and presented. Mesh refinement strategies, allowing tailoring of the refinement of the computational mesh, are also discussed. Accuracy and robustness of the present method are demonstrated on several validation test cases. The method is developed with the prospect of being applied to microfluidic flows and the simulation of microchannel evaporators used for electronics cooling. Therefore, the simulation results for the flow of a bubble in a microchannel are presented and compared to experimental data.     

Cytotoxic dimeric half-sandwich Ru(II), Os(II) and Ir(III) complexes containing the 4,4′-biphenyl-based bridging ligands

A series of dinuclear half-sandwich Ru(II), Os(II) and Ir(III) complexes [Ru₂(μ-Lⁿ)(η⁶-pcym)₂Cl₂](PF₆)₂ ( 1 , 4 ), [Os₂(μ-Lⁿ)(η⁶-pcym)₂Cl₂](PF₆)₂ ( 2 , 5 ) and [Ir₂(μ-Lⁿ)(η⁵-Cp*)₂Cl₂](PF₆)₂ ( 3 , 6 ), based on 4,4′-biphenyl-based bridging Schiff base ligands N,N′-(biphenyl-4,4′-diyldimethylidyne)bis-2-(pyridin-2-yl)methanamine (L¹; for 1 – 3 ) and N,N′-(biphenyl-4,4′-diyldimethylidyne)bis-2-(pyridin-2-yl)ethanamine (L²; for 4 – 6 ) is reported; pcym = 1-methyl-4-(propan-2-yl)benzene, Cp* = pentamethylcyclopentadienyl. The complexes were characterized by relevant analytical techniques (i.e. elemental analysis, FT-IR, NMR, ESI-MS), and their in vitro cytotoxicity was assessed at six cancerous and two non-cancerous (healthy) human cell lines. Overall, complexes 4 – 6 , containing the L² bridging ligand, revealed higher cytotoxicity as compared with 1 – 3 and, thus, they were studied in greater detail. The best-performing complex 6 exceeded at least twice the in vitro cytotoxicity of cisplatin and showed high selectivity towards the cancer cells over the normal ones, including the primary culture of human hepatocytes. In contrast to cisplatin, complexes 4 – 6 did not induce the cell cycle modification of the treated A2780 human ovarian carcinoma cells (studied by flow cytometry and Western blot analysis). High levels of superoxide anion were induced by complexes 4 – 6 at the A2780 cells. The levels of activated forms of Caspase-3 and Caspase-8 at the A2780 cells treated by Ru(II) complex 4 were comparable with cisplatin, while complexes 5 and 6 had only a minor effect on activation of these caspases.     

聚变堆相关的液态金属膜流初步实验研究

在磁约束核聚变堆的面对等离子部件设计中,液态金属锂膜流因具有带走杂质、保护面对等离子固壁等优点而被认为是优选方案之一. 然而,如何克服聚变堆中强磁场环境下产生的磁流体力学效应并形成大面积均匀铺展锂膜流动是目前亟需解决的问题.本文通过搭建室温液 态镓铟锡回路和高温液态锂回路,开展了两种不同特性的液态金属膜流实验, 并采用传统可视化方法获得了展向磁场存在时镓铟锡和锂在导电底板形成的液膜流动表面特征.实验结果 表明: 无磁场时,两种液态金属膜流流动表面波动特性与常规流体膜流均一致, 即随着流动雷诺数的增加表面波动变得更为混乱; 而展向磁场存在时,镓铟锡膜流表面波动变得更为规则, 且沿着磁场方向平行排列,表现为拟二维波动的特征; 而锂膜流却产生了明显的磁流体 力学阻力效应,表现为在流动方向局部产生锂滞留现象, 且滞留点随雷诺数增大向下游移动. 最后通过膜流受力分析,进一步阐述了锂膜流受到比镓铟锡膜流更为严重磁流体力学效应影响的原因.      

纳米阵列中气体驱替液体的流动特征

纳米尺度下气体驱动液体流动特征在纳流控芯片及页岩气开发中具有广泛的应用前景. 利用管径规格为292.8 nm,206.2 nm,89.2 nm,67.0 nm,26.1 nm的氧化铝膜为纳米阵列,进行气驱水实验和单相气体流动实验,分析纳米尺度下气驱水流动特征. 实验表明,纳米阵列中气驱水时气体流量随驱动压力变化经历三个阶段:第一阶段流量缓慢增大,且比单相气体流量降低约一个数量级;第二阶段纳米阵列中的水被大量驱替出,流量迅速增大;第三阶段纳米阵列中的水全部被驱替出,流动特征与单相气体流动保持一致. 分析表明,气驱水第一阶段存在气液界面毛细管力的“钉扎”作用及固液界面相互作用力的影响,是产生非线性流动的主要原因;而一旦“钉扎”作用破坏,气体进入管道推动界面运动,气柱与液柱之间的毛细曲面曲率变化,毛细管力减小,气体流量急剧增大,其中毛细管力随驱替压力增大急剧变化,是造成第二阶段气体流量突变的主要原因.