Prediction of flexoelectricity in BaTiO3 using molecular dynamics simulations

Flexoelectric effect, referring to the strain gradient induced polarization, widely exists in dielectric materials, but its molecular dynamics has not been studied so much so far. In this work, the radial distribution function of BaTiO₃ and the phase transition temperatures have been investigated, and the results show that the core-shell potential model is effective and the structure of BaTiO₃ is stable in a temperature range of 10 K-150 K. Molecular dynamics simulated hysteresis loops of BaTiO₃ show that anisotropy can play an important role in the coercive field. Based on the rational simulation process, the effects of cantilever beam bent angle and fixed length on the polarization are analyzed. It is found that the small bent angle of the curved cantilever beam can give a proportional relationship with a fixed end length and a non-linear relationship is presented when the bent angle is much larger. The prediction of flexoelectric coefficient in BaTiO₃ is 18.5 nC/m. This work provides a computational framework for the study of flexoelectric effect by using molecular dynamics.     

Retention mechanism of calcium ferrite and compositions of ash on selenium during chemical looping gasification

Selenium pollution by coal utilization is of increasing concern. Calcium-iron (Ca–Fe) oxygen carriers (OCs) and alkali metal ions have strong inhibitory effects on selenium, which can reduce the emissions of selenium vapor. The retention mechanisms of selenium by Fe₂O₃, CaFe₂O₄, Ca₂Fe₂O₅ and bottom ash are investigated during chemical looping gasification (CLG). Iron-based OC can oxidize H₂Se(g) to SeO₂(g); furthermore, lattice oxygen is released by Fe₂O₃, contributing to the formation of an Fe–O–Se structure to retain selenium and form selenite. Because calcium ferrite is poorly oxidizing, it cannot oxidize H₂Se(g), but the CaO produced when OCs are reduced can react with H₂Se(g) to form CaSe(s), and this process can be promoted by H₂S(g). The best retention rates reached 32.301% when Ca₂Fe₂O₅ was used. In the cyclic experiment, the selenium retention of the bottom ash gradually increases. Alkali metal ions in bottom ash are the main factor in retaining selenium. Ca²⁺ and Mg²⁺ do not easily vaporize due to their high melting points; therefore, their selenium retention is significantly better than that of K⁺ and Na⁺. This research provided a new idea for the removal of selenium by using OCs and bottom ash particles during CLG.     

具有接种的确定性和随机SIS流行病模型的全局稳定性

We study the stability of endemic equilibriums of the deterministic and stochastic SIS epidemic models with vaccination. The deterministic SIS epidemic model with vaccination was proposed by Li and Ma(2004), for which some sufficient conditions for the global stability of the endemic equilibrium were given in some earlier works. In this paper, we first prove by Lyapunov function method that the endemic equilibrium of the deterministic model is globally asymptotically stable whenever the basic reproduction number is larger than one. For the stochastic version, we obtain some sufficient conditions for the global stability of the endemic equilibrium by constructing a class of different Lyapunov functions.     

Memory access optimization for particle operations in computational fluid dynamics-discrete element method simulations

Computational Fluid Dynamics – Discrete Element Method is used to model gas-solid systems in several applications in energy, pharmaceutical and petrochemical industries. Computational performance bottlenecks often limit the problem sizes that can be simulated at industrial scale. The data structures used to store several millions of particles in such large-scale simulations have a large memory footprint that does not fit into the processor cache hierarchies on current high-performance-computing platforms, leading to reduced computational performance. This paper specifically addresses this aspect of memory access bottlenecks in industrial scale simulations. The use of space-filling curves to improve memory access patterns is described and their impact on computational performance is quantified in both shared and distributed memory parallelization paradigms. The Morton space filling curve applied to uniform grids and k-dimensional tree partitions are used to reorder the particle data-structure thus improving spatial and temporal locality in memory. The performance impact of these techniques when applied to two benchmark problems, namely the homogeneous-cooling-system and a fluidized-bed, are presented. These optimization techniques lead to approximately two-fold performance improvement in particle focused operations such as neighbor-list creation and data-exchange, with ∼ 1.5 times overall improvement in a fluidization simulation with 1.27 million particles.     

Data-driven modeling and parameter estimation of nonlinear systems

The European Physical Journal B - The lattice parameters, formation energy and electronic structure as well as elastic property and Debye temperature of Ca-doped Ti2Ni alloy have been calculated....