NUMERICAL SIMULATION OF FIBROBLAST DUROTAXIS MIGRATION INDUCED BY STIFFNESS GRADIENT OF SUBSTRATE1)

建立了一种细胞趋硬性迁移的理论模型和有限元分析框架,为连续变刚度人工基质的试验设计提供理论依据。考虑了细胞体的黏弹性属性,以及细胞与基质间的配受体动态反应过程,并以配受体合成时间为时间步长,将细胞运动方程化为静力学形式进行求解。对有限元过程提出一种动约束,便于消除其结构矩阵的奇异性。结果表明,模型能够模拟黏着斑内部力的快速波动现象,细胞的运动速度与观测数据一致,可有效模拟20,h以上的长时程问题。     

A new endocytic model based on the co-rotational grid method

Endocytosis plays important roles in many cellular physiological processes, such as metabolism and immune. Many theoretical models have been proposed to study the endocytic process, but little has considered the tensile deformation of the membrane and the actin forces. In this paper, a new endocytic model is proposed based on the co-rotational grid method. Our model gives a direct estimation of the in-plane strain of the plasma membrane and provides a basis for the calculation of further scission process of the vesicle. The results fit well with experimental data in the literature. Moreover, it is suggested that the active forces of actin at the endocytic site is the main mechanism driving the invagination of the plasma membrane.     

On U(1) Gauge Theory Transfer-Matrix in Fourier Basis*

The properties of the transfer-matrix of U(1) lattice gauge theory in the Fourier basis are explored. Among other statements it is shown: (i) the transfer-matrix is block-diagonal, (ii) all consisting vectors of a block are known based on an arbitrary block vector, (iii) the ground-state belongs to the zero-mode's block. The emergence of maximum-points in matrix-elements as functions of the gauge coupling is clarified. Based on explicit expressions for the matrix-elements we present numerical results as tests of our statements.     

Detecting nonclassicality of light via Lieb's concavity

Nonclassical light states are important for both conceptual and practical reasons: they are basic ingredients in testing and exploring quantum foundations, and are crucial resources in quantum technologies. Various useful criteria have been developed to detect nonclassicality in the literature, and several meaningful measures of nonclassicality have been introduced and measured experimentally. In this work, by use of a non-Hermitian generalization of the Wigner-Yanase-Dyson skew information and playing with operator ordering in evaluating average photon number, we develop a novel family of criteria for detecting nonclassicality of light based on Lieb's concavity, which is a deep and powerful result concerning interaction between quantum states and observables. We elucidate the information-theoretic as well as the physical meaning of the criteria, and illustrate their effectiveness in capturing and quantifying nonclassicality of various important light states.     

High sensitive and large dynamic range quasi-distributed sensing system based on slow-light π-phase-shifted fiber Bragg gratings

In this paper, we theoretically analyze the slow-light π-phase-shifted fiber Bragg grating (π-FBG) and its applications for single and multipoint/quasi-distributed sensing. Coupled-mode theory (CMT) and transfer matrix method (TMM) are used to establish the numerical modeling of slow-light π-FBG. The impact of slow-light FBG parameters, such as grating length (L), index change (Δn), and loss coefficient (α) on the spectral properties of π-FBG along with strain and thermal sensitivities are presented. Simulation results show that for the optimum grating parameters L = 50 mm, Δn = 1.5×10⁻⁴, and α = 0.10 m^-1, the proposed slow-light π-FBG is characterized with a peak transmissivity of 0.424, the maximum delay of 31.95 ns, strain sensitivity of 8.380 με^-1, and temperature sensitivity of 91.064 °C^-1. The strain and temperature sensitivity of proposed slow-light π-FBG is the highest as compared to the slow-light sensitivity of apodized FBGs reported in the literature. The proposed grating have the overall full-width at half maximum (FWHM) of 0.2245 nm, and the FWHM of the Bragg wavelength peak transmissivity is of 0.0798 pm. The optimized slow-light π-FBG is used for quasi-distributed sensing applications. For the five-stage strain quasi-distributed sensing network, a high strain dynamic range of value 1469 με is obtained for sensors wavelength spacing as small as 2 nm. In the case of temperature of quasi-distributed sensing network, the obtained dynamic range is of 133 °C. For measurement system with a sufficiently wide spectral range, the π-FBGs wavelength grid can be broadened which results in substantial increase of dynamic range of the system.     

Determining matchdays in sports league schedules to minimize rest differences

Many sports leagues first announce the games to be played in each round and then determine their matchdays as the season progresses. This study focuses on the fairness criterion of minimizing the total rest difference among opposing teams to find the matchdays for an announced schedule. We show that the problem is decomposable into optimizing the rounds separately. We also provide a polynomial-time exact algorithm for canonical schedules.     

Mathematical studies of the solution of Burgers' equations by Adomian decomposition method

In this paper, a novel Adomian decomposition method (ADM) is developed for the solution of Burgers' equation. While high level of this method for differential equations are found in the literature, this work covers most of the necessary details required to apply ADM for partial differential equations. The present ADM has the capability to produce three different types of solutions, namely, explicit exact solution, analytic solution, and semi-analytic solution. In the best cases, when a closed-form solution exists, ADM is able to capture this exact solution, while most of the numerical methods can only provide an approximation solution. The proposed ADM is validated using different test cases dealing with inviscid and viscous Burgers' equations. Satisfactory results are obtained for all test cases, and, particularly, results reported in this paper agree well with those reported by other researchers.     

Optimal control and stability analysis of an online game addiction model with two stages

In this paper, we establish a mathematical model of online game addiction with two stages to research the dynamic properties of it. The existence of all equilibria is obtained, and the basic reproduction number is calculated by the method of next-generation matrix. The global asymptotic stability of disease-free equilibrium (DFE) is proved by comparison principle, and the global asymptotic stability of endemic equilibrium (EE) is proved by constructing an appropriate Lyapunov function. Then we use the Pontryagin's maximum principle to find the optimal solution of the model, so that the number of infected people can be minimized. In numerical simulation, firstly, we validate the global stability of DFE and EE. Secondly, we consider three kind of control measures (treatment, isolation, and education) and divide them into four cases. The models with control and without control are solved numerically by using forward and backward sweep Runge-Kutta method. In order to achieve the best control effect, we suggest that three kind of measures should be used simultaneously according to the optimal control strategy.     

Chebyshev spectral variational integrator and applications

The Chebyshev spectral variational integrator(CSVI) is presented in this paper. Spectral methods have aroused great interest in approximating numerically a smooth problem for their attractive geometric convergence rates. The geometric numerical methods are praised for their excellent long-time geometric structure-preserving properties.According to the generalized Galerkin framework, we combine two methods together to construct a variational integrator, which captures the merits of both methods. Since the interpolating points of the variational integrator are chosen as the Chebyshev points,the integration of Lagrangian can be approximated by the Clenshaw-Curtis quadrature rule, and the barycentric Lagrange interpolation is presented to substitute for the classic Lagrange interpolation in the approximation of configuration variables and the corresponding derivatives. The numerical float errors of the first-order spectral differentiation matrix can be alleviated by using a trigonometric identity especially when the number of Chebyshev points is large. Furthermore, the spectral variational integrator(SVI) constructed by the Gauss-Legendre quadrature rule and the multi-interval spectral method are carried out to compare with the CSVI, and the interesting kink phenomena for the Clenshaw-Curtis quadrature rule are discovered. The numerical results reveal that the CSVI has an advantage on the computing time over the whole progress and a higher accuracy than the SVI before the kink position. The effectiveness of the proposed method is demonstrated and verified perfectly through the numerical simulations for several classical mechanics examples and the orbital propagation for the planet systems and the Solar system.