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.     

Fixed Points,Continuous Selections,and Existence of Solution of Optimization-Related Problems

We propose a definition of a generalized type of Knaster–Kuratowski–Mazurkiewicz (KKM) mappings, called a weak T-KKM mapping, and a corresponding weak KKM property. A new extension of the Fan–Glicksberg fixed-point theorem is established. Sufficient conditions for the existence of a continuous selection, a fixed point of a composition, and a coincidence point are also provided. Then, we use the obtained results to study the existence of solutions to various optimization-related problems. Discussions and detailed examples are included as well to compare our results with existing ones and to explain their advantages in many situations.     

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.     

Runge–Kutta Methods for Systems of Differential Equation with Piecewise Continuous Arguments: Convergence and Stability

This work deals with the convergence and stability of Runge–Kutta methods for systems of differential equation with piecewise continuous arguments x^′(t) = Px(t)+Qx([t+1∕2]) under two cases for coe?cient matrix. First, when P and Q are complex matrices, the su?cient condition under which the analytic solution is asymptotically stable is given. It is proven that the Runge–Kutta methods are convergent with order p. Moreover, the su?cient condition under which the analytical stability region is contained in the numerical stability region is obtained. Second, when P and Q are commutable Hermitian matrices, using the theory of characteristic, the necessary and su?cient conditions under which the analytic solution and the numerical solution are asymptotically stable are presented, respectively. Furthermore, whether the Runge–Kutta methods preserve the stability of analytic solution are investigated by the theory of Padé approximation and order star. To demonstrate the theoretical results, some numerical experiments are adopted.     

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.