High-order curvilinear mesh generation technique based on an improved radius basic function approach

A high-order curvilinear hybrid mesh generation technique is developed for high-order numerical method (eg, discontinuous Galerkin method) applications to improve the accuracy for problems with curve boundary. The grid generation technique is based on an improved radius basic function (RBF) approach by which the straight-edge mesh is converted into high-order curve mesh. Firstly, an initial straight-edge mesh is prepared by traditional grid generation software. Then, high-order interpolation points are inserted into the mesh entities such as edges, faces, and cells according to the final demand of mesh order. To preserve the original geometry, the inserted points on solid wall are then projected onto the CAD model using an open source tool “Open Cascade.” Finally, other inserted points in the field near the solid wall are moved to appropriate positions by the improved RBF approach to avoid tangled cells. If we use the original RBF approach, then the inserted points on the edge and face entities normal to the solid boundary in the region of boundary layer will move to improper positions. To overcome this problem, a weighting based on the local grid aspect ratio between normal direction and tangential direction is introduced into the baseline RBF approach. Three typical configurations are tested to validate the mesh generator. Meanwhile, a third-order solution of subsonic flow over an analytical 3D body of revolution in the second International Workshop on High-Order CFD Methods is supplied by a discontinuous Galerkin solver. These numerical tests demonstrate the potential capability of present technique for high-order simulations of complex geometries.     

Cover Image

This paper presents a framework for incorporating arbitrary implicit multistep schemes into the lattice Boltzmann method. While the temporal discretization of the lattice Boltzmann equation is usually derived using a second-order trapezoidal rule, it appears natural to augment the time discretization by using multistep methods. The effect of incorporating multistep methods into the lattice Boltzmann method is studied in terms of accuracy and stability. Numerical tests for the third-order accurate Adams-Moulton method and the second-order backward differentiation formula show that the temporal order of the method can be increased when the stability properties of multistep methods are considered in accordance with the second Dahlquist barrier.     

Nonlinear corrector for Reynolds-averaged Navier-Stokes equations

The scope of this paper is to present a nonlinear error estimation and correction for Navier-Stokes and Reynolds-averaged Navier-Stokes equations. This nonlinear corrector enables better solution or functional output predictions at fixed mesh complexity and can be considered in a mesh adaptation process. After solving the problem at hand, a corrected solution is obtained by solving again the problem with an added source term. This source term is deduced from the evaluation of the residual of the numerical solution interpolated on the h/2 mesh. To avoid the generation of the h/2 mesh (which is prohibitive for realistic applications), the residual at each vertex is computed by local refinement only in the neighborhood of the considered vertex. One of the main feature of this approach is that it automatically takes into account all the properties of the considered numerical method. The numerical examples point out that it successfully improves solution predictions and yields a sharp estimate of the numerical error. Moreover, we demonstrate the superiority of the nonlinear corrector with respect to linear corrector that can be found in the literature.     

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.