Time to most recent common ancestor for stationary continuous state branching processes with immigration

Motivated by sample path decomposition of the stationary continuous state branching process with immigration, a general population model is considered using the idea of immortal individual. We compute the joint distribution of the random variables： the time to the most recent common ancestor （MRCA）, the size of the current population, and the size of the population just before MRCA. We obtain the bottleneck effect as well. The distribution of the number of the oldest families is also established. These generalize the results obtained by Y. T. Chen and J. F. Delmas.     

Mesoscopic Numerical Study on Flow Boiling Heat Transfer Performance in Channels With Multiple Rectangular Heaters北大核心CSCD

The flow boiling phenomenon in a channel with multiple rectangular heaters under a constant wall temperature was numerically studied with the lattice Boltzmann method. The effects of spacings between heaters, heater lengths and heater surface wettabilities on the bubble morphology, the bubble area and the heat flux on the heater surface, were studied. The results show that, the bubble growth rate increases with the spacing between heaters. The larger the bubble area is, the earlier the nucleated bubbles will leave the heater surface. The corresponding boiling heat transfer performance increases by 12% with the spacing between heaters growing from 250 lattices to 1 000 lattices. On the other hand, the longer the heater length is, the earlier the bubble will nucleate and leave the heater surface, and the better the boiling heat transfer performance will be. The boiling heat transfer performance increases by 13% with the heater length rising from 16 lattices to 22 lattices. In addition, the bubble nucleates later on the hydrophilic surface than on the hydrophobic surface. Compared with the hydrophilic surface, the hydrophobic surface retains residual bubbles after the leaving of bubbles from the heater. The average heat flux and the bubble area of the hydrophilic surface are less than those of the hydrophobic surface. With the contact angle changing from 77° to 120°, the heat transfer performance increases by 26%. Finally, the orthogonal test results indicate that, the wettability of the heat exchanger surface has the greatest influence on the flow boiling heat transfer performance, while the heater length has the least influence. © 2022 Editorial Office of Applied Mathematics and Mechanics. All rights reserved.     

Uzawa iteration method for stokes type variational inequality of the second kind

In this paper,the Uzawa iteration algorithm is applied to the Stokes problem with nonlinear slip boundary conditions whose variational formulation is the variational inequality of the second kind.Firstly, the multiplier in a convex set is introduced such that the variational inequality is equivalent to the variational identity.Moreover,the solution of the variational identity satisfies the saddle-point problem of the Lagrangian functional ?.Subsequently,the Uzawa algorithm is proposed to solve the solution of the saddle-point problem. We show the convergence of the algorithm and obtain the convergence rate.Finally,we give the numerical results to verify the feasibility of the Uzawa algorithm.     

Bregman iterative model using the G-norm

In this paper, we analyze the Bregman iterative model using the G-norm. Firstly, we show the convergence of the iterative model. Secondly, using the source condition and the symmetric Bregman distance, we consider the error estimations between the iterates and the exact image both in the case of clean and noisy data. The results show that the Bregman iterative model using the G-norm has the similar good properties as the Bregman iterative model using the L2-norm.     

Inverse Fluid-solid Interaction Scattering Problem Using Phased and Phaseless Far Field Data

We consider the inverse fluid-solid interaction scattering of incident plane wave from the knowledge of the phased and phaseless far field patterns.For the phased data,one direct sampling method for location and shape reconstruction is proposed.Only inner product is involved in the computation,which makes it very simple and fast to be implemented.With the help of the factorization of the far field operator,we give a lower bound of the proposed indicator functional for the sampling points inside the elastic body.While for the sampling points outside,we show that the indicator functional decays like the Bessel function as the points go away from the boundaries of the elastic body.We also show that the proposed indicator functional continuously dependents on the far field patterns,which further implies that the novel sampling method is extremely stable with respect to data error.For the phaseless data,to overcome the translation invariance,we consider the scattering of point sources simultaneously.By adding a reference sound-soft obstacle into the scattering system,we show some uniqueness results with phaseless far field data.Numerically,we introduce a phase retrieval algorithm to retrieve the phased data without the additional obstacle.The novel phase retrieval algorithm can also be combined with the sampling method for phased data.We also design two novel direct sampling methods using the phaseless data directly.Finally,some numerical simulations in two dimensions are conducted with noisy data,and the results further verify the effectiveness and robustness of the proposed numerical methods.     

The Mathematical Model and Research Progress of the Boundary Layer Flashback in Premixed Combustion北大核心CSCD

Flashback is a key problem influencing the normal operation of power equipment such as gas turbines. As one of the main mechanisms that cause flashback, the boundary layer flashback has an important effect on the design and operation of gas turbine combustors and other combustion devices. Since the critical gradient model for the boundary layer flashback was put forward by Lewis et al. in 1945, the theoretical models for the boundary layer flashback, such as the Peclet number model, the Damköhler number model and the flame angle theory, were developed one after another. However, these theoretical models still need improvements. Until now, the theoretical models for the boundary layer flashback are still in continuous development and modification. The history of the boundary layer flashback was reviewed, and the background, pertinence and shortcomings of the theoretical models were elucidated in the order of the model establishment time. In addition, the development status and research progress of the theoretical models for the boundary layer flashback in recent years were summarized, especially the progress made with new methods such as numerical simulation and statistical analysis. Further, the theoretical research direction and breakthrough points of the combustion boundary layer flashback at present and in the future were put forward. © 2023 Editorial Office of Applied Mathematics and Mechanics. All rights reserved.     

Interlaminar Crack Propagation Analysis of ENF Specimens Based on the Cohesive Zone Model北大核心CSCD

Based on the classical laminated plate theory and the cohesive zone model, a theoretical model for general delamination cracked laminates was established for crack propagation of pure mode Ⅱ ENF specimens. Compared with the conventional beam theory, the proposed model fully considered the softening process of the cohesive zone and introduced the nonlinear behavior of ENF specimens before failure. The predicted failure load is smaller than that under the beam theory and closer to the experimental data in literatures. Compared with the beam theory with only fracture toughness considered, the proposed model can simultaneously analyze the influences of the interface strength, the fracture toughness and the initial interface stiffness on the load-displacement curves in ENF tests. The results show that, the interface strength mainly affects the mechanical behavior of specimens before failure, but has no influence on crack propagation. The fracture toughness is the main parameter affecting crack propagation, and the initial interface stiffness only affects the linear elastic loading stage. The cohesive zone length increases with the fracture toughness and decreases with the interface strength. The effect of the interface strength on the cohesive zone length is more obvious than that of the fracture toughness. When the adhesive zone tip reaches the half length of the specimen, the adhesive zone length will decrease to a certain extent. Copyright ©2022 Applied Mathematics and Mechanics. All rights reserved.     

ON CONTRACTION AND SEMI-CONTRACTION FACTORS OF GSOR METHOD FOR AUGMENTED LINEAR SYSTEMS

The generalized successive overrelaxation （GSOR） method was presented and studied by Bai, Parlett and Wang [Numer. Math. 102（2005）, pp.1-38] for solving the augmented system of linear equations, and the optimal iteration parameters and the corresponding optimal convergence factor were exactly obtained. In this paper, we further estimate the contraction and the semi-contraction factors of the GSOR method. The motivation of the study is that the convergence speed of an iteration method is actually decided by the contraction factor but not by the spectral radius in finite-step iteration computations. For the nonsingular augmented linear system, under some restrictions we obtain the contraction domain of the parameters involved, which guarantees that the contraction factor of the GSOR method is less than one. For the singular but consistent augmented linear system, we also obtain the semi-contraction domain of the parameters in a similar fashion. Finally, we use two numerical examples to verify the theoretical results and the effectiveness of the GSOR method.     

ON THE FINITE ELEMENT APPROXIMATION OF SYSTEMS OF REACTION-DIFFUSION EQUATIONS BY p/hp METHODS

We consider the approximation of systems of reaction-diffusion equations, with the finite element method. The highest derivative in each equation is multiplied by a parameter ε∈ （0, 1], and as ε → 0 the solution of the system will contain boundary layers. We extend the analysis of the corresponding scalar problem from [Melenk, IMA J. Numer. Anal. 17（1997）, pp. 577-601], to construct a finite element scheme which includes elements of size O（εp） near the boundary, where p is the degree of the approximating polynomials. We show that, under the assumption of analytic input data, the method yields exponential rates of convergence, independently of ε, when the error is measured in the energy norm associated with the problem. Numerical computations supporting the theory are also presented, which also show that the method yields robust exponential convergence rates when the error in the maximum norm is used.     

THE AVERAGE ABUNDANCE FUNCTION WITH MUTATION OF THE MULTI-PLAYER SNOWDRIFT EVOLUTIONARY GAME MODEL

This article explores the characteristics of the average abundance function with mutation on the basis of the multi-player snowdrift evolutionary game model by analytical analysis and numerical simulation.The specific field of this research concerns the approximate expressions of the average abundance function with mutation on the basis of different levels of selection intensity and an analysis of the results of numerical simulation on the basis of the intuitive expression of the average abundance function.In addition,the biological background of this research lies in research on the effects of mutation,which is regarded as a biological concept and a disturbance to game behavior on the average abundance function.The mutation will make the evolutionary result get closer to the neutral drift state.It can be deduced that this affection is not only related to mutation,but also related to selection intensity and the gap between payoff and aspiration level.The main research findings contain four aspects.First,we have deduced the concrete expression of the expected payoff function.The asymptotic property and change trend of the expected payoff function has been basically obtained.In addition,the intuitive expression of the average abundance function with mutation has been obtained by taking the detailed balance condition as the point of penetration.It can be deduced that the effect of mutation is to make the average abundance function get close to 1/2.In addition,this affection is related to selection intensity and the gap.Secondly,the first-order Taylor expansion of the average abundance function has been deduced for when selection intensity is sufficiently small.The expression of the average abundance function with mutation can be simplified from a composite function to a linear function because of this Taylor expansion.This finding will play a significant role when analyzing the results of the numerical simulation.Thirdly,we have obtained the approximate expressions of the average abundance function corresponding to small and large selection intensity.The significance of the above approximate analysis lies in that we have grasped the basic characteristics of the effect of mutation.The effect is slight and can be neglected when mutation is very small.In addition,the effect begins to increase when mutation rises,and this effect will become more remarkable with the increase of selection intensity.Fourthly,we have explored the influences of parameters on the average abundance function with mutation through numerical simulation.In addition,the corresponding results have been explained on the basis of the expected payoff function.It can be deduced that the influences of parameters on the average abundance function with mutation will be slim when selection intensity is small.Moreover,the corresponding explanation is related to the first-order Taylor expansion.Furthermore,the influences will become notable when selection intensity is large.     

The generalized hierarchical product of graphs

A generalization of both the hierarchical product and the Cartesian product of graphs is introduced and some of its properties are studied. We call it the generalized hierarchical product. In fact, the obtained graphs turn out to be subgraphs of the Cartesian product of the corresponding factors. Thus, some well-known properties of this product, such as a good connectivity, reduced mean distance, radius and diameter, simple routing algorithms and some optimal communication protocols, are inherited by the generalized hierarchical product. Besides some of these properties, in this paper we study the spectrum, the existence of Hamiltonian cycles, the chromatic number and index, and the connectivity of the generalized hierarchical product.     

The q-deformed harmonic oscillator, coherent states, and the uncertainty relation

For a q-deformed harmonic oscillator, we find explicit coordinate representations of the creation and annihilation operators, eigenfunctions, and coherent states (the last being defined as eigenstates of the annihilation operator). We calculate the product of the “coordinate-momentum” uncertainties in q-oscillator eigenstates and in coherent states. For the oscillator, this product is minimum in the ground state and equals 1/2, as in the standard quantum mechanics. For coherent states, the q-deformation results in a violation of the standard uncertainty relation; the product of the coordinate-and momentum-operator uncertainties is always less than 1/2. States with the minimum uncertainty, which tends to zero, correspond to the values of λ near the convergence radius of the q-exponential. __________ Translated from Teoreticheskaya i Matematicheskaya Fizika, Vol. 147, No. 2, pp. 315–322, May, 2006.     

Gaussian Versus Optimal Integration of Analytic Functions

We consider error estimates for optimal and Gaussian quadrature formulas if the integrand is analytic and bounded in a certain complex region. First, a simple technique for the derivation of lower bounds for the optimal error constants is presented. This method is applied to Szeg?-type weight functions and ellipses as regions of analyticity. In this situation, the error constants for the Gaussian formulas are close to the obtained lower bounds, which proves the quality of the Gaussian formulas and also of the lower bounds. In the sequel, different regions of analyticity are investigated. It turns out that almost exclusively for ellipses, the Gaussian formulas are near-optimal. For classes of simply connected regions of analyticity, which are additionally symmetric to the real axis, the asymptotic of the worst ratio between the error constants of the Gaussian formulas and the optimal error constants is calculated. As a by-product, we prove explicit lower bounds for the Christoffel-function for the constant weight function and arguments outside the interval of integration. September 7, 1995. Date revised: October 25, 1996.     

The permanence and extinction of a nonlinear growth rate single-species non-autonomous dispersal models with time delays

In this paper, we consider the effect of diffusion on the permanence and extinction of a non-autonomous nonlinear growth rate single-species dispersal model with time delays. Firstly, the sufficient conditions of the permanence and extinction of the species are established, which shows if the growth rate and dispersal coefficients is suitable, the species is permanent, on the contrary, it is extinction. Secondly, an interesting result is established, that is, if only the species in some patches even in one patch is permanent, then it is also permanent in other patches. Finally, some examples together with their numerical simulations show the feasibility of our main results.     

Rolling-ball method for estimating the boundary of the support of a point-process intensity

We suggest a generalisation of the convex-hull method, or ‘DEA’ approach, for estimating the boundary or frontier of the support of a point cloud. Figuratively, our method involves rolling a ball around the cloud, and using the equilibrium positions of the ball to define an estimator of the envelope of the point cloud. Constructively, we use these ideas to remove lines from a triangulation of the points, and thereby compute a generalised form of a convex hull. The radius of the ball acts as a smoothing parameter, with the convex-hull estimator being obtained by taking the radius to be infinite. Unlike the convex-hull approach, however, our method applies to quite general frontiers, which may be neither convex nor concave. It brings to these contexts the attractive features of the convex hull: simplicity of concept, rotation-invariance, and ready extension to higher dimensions. It admits bias corrections, which we describe and illustrate through implementation.     

Moment information for probability distributions,without solving the moment problem,II: Main-mass,tails and shape approximation

How much information does a small number of moments carry about the unknown distribution function? Is it possible to explicitly obtain from these moments some useful information, e.g., about the support, the modality, the general shape, or the tails of a distribution, without going into a detailed numerical solution of the moment problem? In this, previous and subsequent papers, clear and easy to implement answers will be given to some questions of this type. First, the question of how to distinguish between the main-mass interval and the tail regions, in the case we know only a number of moments of the target distribution function, will be addressed. The answer to this question is based on a version of the Chebyshev–Stieltjes–Markov inequality, which provides us with upper and lower, moment-based, bounds for the target distribution. Then, exploiting existing asymptotic results in the main-mass region, an explicit, moment-based approximation of the target probability density function is provided. Although the latter cannot be considered, in general, as a satisfactory solution, it can always serve as an initial approximation in any iterative scheme for the numerical solution of the moment problem. Numerical results illustrating all the theoretical statements are also presented.     

Is prediction possible? Chaotic behavior of Multiple Equilibria Regulation Model in cellular automata topology

In this article, we present the Multiple Equilibria Regulation (MER) Model in cellular automata topology. As argued in previous explorations of the model, for certain parameter values, the behavior of the system exhibits transient chaos (namely, the system is unpredictable but ends in a final steady state). In order to approach empirical reality, we introduce a cellular automata topology. Examining the outcome of the simulations leads us to conclude that for certain parameter values tested, the system yields chaotic behavior. Thus, cellular automata contribution has proven crucial, because the introduced topology converts the behavior of the system from transient chaos to “pure” chaos, i.e., the system is not only unpredictable on the long run but, in addition, it will never rest in a final steady state. According to these findings, authors argue the theoretical hypothesis that the urge for “prediction” in social sciences should be reconsidered in terms of “predictability horizon”. © 2004 Wiley Periodicals, Inc. Complexity 10: 23–36, 2004     

The limiting state of a rigidly fixed nonlinearly elastic multilayer rod

The loss of the load-carrying capacity of a nonlinearly elastic multilayer rod is investigated. The rod, whose layers have various thickness and are made of different materials, is rigidly fixed at both its ends. Rigid contact conditions between the layers are assumed. The problem posed is solved by using the variational method of mixed type in combination with the Rayleigh-Ritz method. The initial analysis is reduced to the solution of the Cauchy problem for a nonlinear ordinary differential equation solved for the first derivative. As the initial condition, the maximum initial eccentricity of the rod is assumed. In the case of zero eccentricity, the Shanley critical force for an axially compressed rod is determined. For a three-layer rod whose outer layers have equal thickness and are made of the same material, numerically, for various degrees of nonlinearity, the effect of physicomechanical and geometric parameters on the critical load of buckling instability is determined. It is found that, by matching the heterogeneity of the rod, it is possible to raise its load-carrying capacity. __________ Translated from Mekhanika Kompozitnykh Materialov, Vol. 42, No. 3, pp. 347–360, May–June, 2006.