A two-grid finite element method for the Allen-Cahn equation with the logarithmic potential

In this paper, we present a two-grid finite element method for the Allen-Cahn equation with the logarithmic potential. This method consists of two steps. In the first step, based on a fully implicit finite element method, the Allen-Cahn equation is solved on a coarse grid with mesh size H. In the second step, a linearized system whose nonlinear term is replaced by the value of the first step is solved on a fine grid with mesh size h. We give the energy stabilities of the traditional finite element method and the two-grid finite element method. The optimal convergence order of the two-grid finite element method in H¹ norm is achieved when the mesh sizes satisfy h = O(H²). Numerical examples are given to demonstrate the validity of the proposed scheme. The results show that the two-grid method can save the CPU time while keeping the same convergence rate.     

Adaptive Evolution Analysis of a Predator-Prey Model With Group Defense北大核心CSCD

Based on the theoretical framework of adaptive dynamics, the evolution of the predator-prey model with functional response of group defense effect on the predator handling time, was investigated. Firstly, in view of the interaction of predator populations with interspecific competition, the evolutionary conditions for a single predator population to split into 2 populations with different strategies through evolutionary branching were given. Secondly, when the ecological equilibrium of the model is unstable and the system has a limit cycle, the population will have strong coexistence under large mutation, but this coexistence will be evolutionarily unstable. Finally, the conclusions for the model with Holling-Ⅱ type functional response were compared. The results indicate that, with a sufficiently large prey carrying capacity, group defense effects can evolutionarily lead to the extinction of predators. © 2023 Editorial Office of Applied Mathematics and Mechanics. All rights reserved.     

Absorbing Interface Conditions for the Simulation of Wave Propagation on Non-Uniform Meshes

We proposed absorbing interface conditions for the simulation of linear wave propagation on non-uniform meshes. Based on the superposition principle of second-order linear wave equations, we decompose the interface condition problem into two subproblems around the interface: for the first one the conventional artificial absorbing boundary conditions is applied, while for the second one, the local analytic solutions can be derived. The proposed interface conditions permit a two-way transmission of low-frequency waves across mesh interfaces which can be supported by both coarse and fine meshes, and perform a one-way absorption of high-frequency waves which can only be supported by fine meshes when they travel from fine mesh regions to coarse ones. Numerical examples are presented to illustrate the efficiency of the proposed absorbing interface conditions.     

Input-to-state stability and stabilization for switched nonlinear positive systems

Input-to-state stability (ISS) analysis and stabilization are concerned in this paper for switched nonlinear positive systems (SNPS), where the deterministic and random switching are both included. For general SNPS, switched affine nonlinear positive systems (SANPS) and switched linear positive systems (SLPS) with deterministic and some kinds of random ”slow” switching, some criterions on ISS are provided. From the criterions for SANPS and SLPS, the ISS properties can be judged just by the differential, algebraic and switching characteristics of the systems. Further, based on the criterions for SANPS and SLPS, some state feedback controllers are designed such that the closed-loop systems be positive, ISS or ISS in some stochastic senses. Four simulation examples verify the validity of our results.     

Effects of the formation time of parton shower on jet quenching in heavy-ion collisions

Jet quenching has successfully served as a hard probe to study the properties of Quark-Gluon Plasma (QGP). As a multi-particle system, jets require time to develop from a highly virtual parton to a group of partons close to mass shells. In this study, we present a systematical analysis on the effects of this formation time on jet quenching in relativistic nuclear collisions. Jets from initial hard scatterings were simulated with Pythia, and their interactions with QGP were described using a Linear Boltzmann Transport (LBT) model that incorporates both elastic and inelastic scatterings between jet partons and the thermal medium. Three different estimations of the jet formation time were implemented and compared, including instantaneous formation, formation from single splitting, and formation from sequential splittings, before which no jet-medium interaction was assumed. We found that deferring the jet-medium interaction with a longer formation time not only affects the overall magnitude of the nuclear modification factor of jets but also its dependence on the jet transverse momentum.     

A lattice Boltzmann method for KDV equation

We prepose a 5-bit lattice Boltzmann model for KdV equation. Using Chapman-Enskog expansion and multiscale technique, we obtained high order moments of equilibrium distribution function, and the 3rd dispersion coefficient and 4th order viscosity. The parameters of this scheme can be determined by analysing the energy dissipation. The project supported by the Foundation of the Laboratory for Nonlinear Mechanics of Continuous Media, Institute of Mechanics, Chinese Academy of Sciences     

Mechanical model of lung parenchyma as a two-phase porous medium

The anatomy and geometry of the lung at the micro- and macroscopic level have been described briefly. A notion of lung parenchyma — a macroscopically continuous medium whose mechanical properties result from those of microstructural components — has been adapted. Simplifying assumptions propounded in the constitutive model have been discussed. Two phases have been distinguished in the medium: the solid phase — a highly deformable, nonlinearly elastic skeleton in the form of a thin-walled tissue structure on the micro-scale — and the fluid phase — perfect gas (air) filterating through the structure. General constitutive relations for both phases and their mechanical interactions have ben formulated. Further, the fundamental set of differential equations of the quasi-static coupled problem has been developed. Large deformations, material nonlinearities, and dependence of permeability on skeleton deformation have been included. Matrix formulation of the problem has been presented from the point of view of the finite element method. An implicit iterative time integration scheme has been proposed. The algorithm has been illustrated with results of simple numerical tests.     

停歇对疲劳裂纹扩展超载迟滞效应影响的研究

本文用 16 Mn和 A 537 两种普通低合金钢试件,较系统地研究了停歇对疲劳超载迟滞效应的影响。研究表明,疲劳载荷中的停歇对超载迟滞效应有明显影响,其影响的程度大小与停歇时的载荷 P_H,停歇时间 t_H 等因素密切相关;不但可能使试件的疲劳寿命减少,也可能使试件的疲劳寿命增加。在工程结构疲劳设计、分析及应用超载迟滞效应时应予以注意。     

Investigation on time scales in a low Reynolds number jet flow using particle-tracking velocimetry

The near flow field of an axisymmetric water jet at Reynolds numbers between 2000 and 5000 is investigated using Particle-Tracking Velocimetry. Measurements are taken in the longitudinal section (along the mean flow) and in cross-sections (orthogonal to the mean flow). From the former, correlation coefficients of the two in-plane velocity components in a Lagrangian framework are obtained: thus Lagrangian integral scales can be computed. Those of the streamwise velocity (axial) component increase on moving away from the centreline, whereas the opposite happens for the vertical velocity (radial) component: integral time scales of the two components are almost equal at the interface between jet and ambient fluids. On the other hand, integral scales are almost constant or increase slightly with the axial direction. In cross-sections, fluid ejection and injection from the jet centreline are observed to be connected to counter-rotating vortices (mushroom): their number and size change with Reynolds number in agreement with results from other authors. The maximum ejection velocity (orthogonal to the mean jet flow), at 3 nozzle diameters downstream of the outlet, is found to be one half of the mean outlet velocity.     

Numerical study of steady/unsteady flow and heat transfer in porous media using a characteristics-based matrix-free implicit FV method on unstructured grids

In this study, a high-resolution characteristic-based finite-volume (FV) method on unstructured grids [Int. J. Numer. Method Eng. 50 (2001) 11; Int. J. Heat Fluid Flow 21 (2000) 432] is extended by a matrix-free implicit dual-time stepping scheme for the numerical simulation of steady and unsteady flow and heat transfer with porous media. The method has been used to study the characteristics of a complex problem: flow and heat transfer in a channel with multiple discrete porous blocks, which was originally proposed by Huang and Vafai [J. Thermophys. Heat Transfer 8 (3) (1994) 563]. In addition, flow and heat transfer in a channel partially or fully filled with porous layers and containing solid protruding blocks with constant heat flux on its lower surface are also investigated in details. Hydrodynamic and heat transfer results are reported for both steady and transient flow cases. In particular, the effects of Darcy and Reynolds numbers on heat transfer augmentation and pressure loss are studied. An in-depth discussion of the formation and variation of recirculation is presented and the existence of optimum porous insert is demonstrated. At high Reynolds numbers the flow in the porous channel exhibits a cyclic characteristics although unlike the non-porous channel flow, the cyclic vortex development is only restricted to a small area behind the last solid block, while temperature changes more slowly and does not exhibit cyclic variations over a long period of time. It is shown that for all the cases studied altering some parametric values can have significant and interesting effects on both flow pattern as well as heat transfer characteristics.