Loading...
Association Journal of CSIAM
Supervised by Ministry of Education of PRC
Sponsored by Xi'an Jiaotong University
ISSN 1005-3085  CN 61-1269/O1

Table of Content

    15 June 2022, Volume 39 Issue 3 Previous Issue   
    Analysis of Solid Fluidized Hydrate Production in Double-layer Continuous Tube
    YANG Yan, BAI Jianchao, WANG Guorong, CAI Mingjie, MAO Liangjie
    2022, 39 (3):  341-356.  doi: 10.3969/j.issn.1005-3085.2022.03.001
    Abstract ( 74 )   PDF (1556KB) ( 177 )   Save
    In order to study the multi-phase wellbore flow law of the two-layer continuous tube lift hydrate, a multi-phase wellbore flow model using the two-layer continuous tube solid fluidization was established in this paper, and the model was solved by the finite difference method. Based on the above model, the multi-phase wellbore flow behavior of gas hydrate production in a case well in the South China Sea using solid fluidization of two-layer continuous pipe was simulated. The results show that the solid phase content and liquid phase content are relatively stable and the gas content is 0 when the fluid rises to a height of more than 700 meters. When the fluid rises to a height more than 700 meters, the solid phase hydrate begins to decompose gradually. When the density and displacement of drilling fluid are too high or too low, the bottomhole pressure will escape from the safety window and threaten the safety of well control. Therefore, the reasonable selection of drilling parameters is particularly important for ensuring well control safety.
    Related Articles | Metrics
    Temperature Modeling and Weather Derivative Pricing with Forward-looking Information---A Case Study of Zhengzhou
    LI Peng, NIU Zhikang, ZHONG Weizhou
    2022, 39 (3):  357-378.  doi: 10.3969/j.issn.1005-3085.2022.03.002
    Abstract ( 127 )   PDF (5466KB) ( 185 )   Save
    Through incorporating forward-looking information, this paper constructs a consistent two-factor temperature model for Zhengzhou, and prices weather derivatives based on the model. First, we analyze the historical weather forecast data, and verify the rationality of the consistent two-factor model. Using the historical temperature data and historical temperature forecast data of Zhengzhou, we then construct a consistent two-factor model and derive the approximate pricing formulas with the market price of risk extracted from utility indifference valuation. In the end, we apply Monte Carlo simulations to compute the price surface for temperature options. This paper presents a complete pricing procedure based on a consistent two-factor temperature model, and is of both theoretical and practical significance for weather risk management and the development of China's weather derivative markets.
    Related Articles | Metrics
    Optimization Model of Group Decision-making for Subjective Scoring Competition Ranking
    GUO Dongwei, DING Genhong
    2022, 39 (3):  379-388.  doi: 10.3969/j.issn.1005-3085.2022.03.003
    Abstract ( 173 )   PDF (256KB) ( 246 )   Save
    Through examining the characteristics of the judges' scoring and the group decision-making mechanism, it is found that in the incomplete scoring system, especially in the subjective scoring competition, there exists a large error in the ranking results of the traditional method or the T score method. In order to improve the scientificity of ranking, a mathematical model is established to distribute the papers evenly among the judges, and then the weight of each judge is determined by the minimum sum of the squares of ranges. Finally, the ranking is determined according to the weighted average value of each paper's scores. Simulation results show that, compared with the traditional and standard score methods, the new method can improve the coincidence of subjective and objective ranking and reduces the disorder degree. By calculating the Euclidean distance, it is proved that the new method can reduce the experts' scoring error and the resulting controversy.
    Related Articles | Metrics
    Decision Making Method Based on Intuitionistic Fuzzy Tensor and Its Application
    DENG Shengyue, TAN Jintao, SONG Yuqin, SHEN Jing
    2022, 39 (3):  389-400.  doi: 10.3969/j.issn.1005-3085.2022.03.004
    Abstract ( 69 )   PDF (195KB) ( 116 )   Save
    The intuitionistic fuzzy decision-making problem with high-dimensional data is stu-died in this paper. The general form and operational laws of the intuitionistic fuzzy tensor are defined, the generalized intuitionistic fuzzy weighted averaging (GIFWA) operator based on the intuitionistic fuzzy tensor is subsequently established. By exploring the basic properties of the GIFWA operator, it is proved that the GIFWA operator is a general form of the intuition-istic fuzzy weighted average operator. A novel decision-making method based on the GIFWA operator is presented to solve the intuitionistic fuzzy decision problems with high-dimensional data. One typical example is also provided to demonstrate the efficiency and applicability of our proposed method.
    Related Articles | Metrics
    Multiple Mean Change-points Estimation Based on Screening and Ranking Algorithm
    LI Yang, WU Mixia, HU Yao, YANG Chao
    2022, 39 (3):  401-412.  doi: 10.3969/j.issn.1005-3085.2022.03.005
    Abstract ( 87 )   PDF (431KB) ( 228 )   Save
    The multiple change-points estimation problem is a hot issue in current statistics, and there are many algorithms in the literature. Among them, Screening and Ranking algorithm (SaRa) has attracted wide attention due to its fast detection and high precision characteristics. However, this algorithm tends to be conservative in the threshold selection of the screening procedure. The reason is that the variance in SaRa is separately estimated in each segment process. The main purpose of this paper is to improve SaRa. Firstly, a global estimate of the variance is calculated through local polynomial approximation with the bandwith selected by the cross validation method. The initial change-points are obtained from screening based on the improved threshold. Then, to order those points in terms of the local diagnostic function values, the number of final change-points is determined by maximizing the MBIC. Numerical results show that the proposed algorithm has high accuracies in the estimation of the number and locations of change points in comparison to existing methods. Finally, this method is app-lied to the actual traffic flow data of Shenzhen city. The distribution characteristics of change points on working days and non-working days in this area are analyzed, which can provide some guidance for traffic control departments and travelers.
    Related Articles | Metrics
    Improved Full-Newton Step Infeasible Interior Point Algorithm for Weighted Linear Complementarity Problem
    CHI Xiaoni, LIU Sanyang, WANG Boda
    2022, 39 (3):  413-427.  doi: 10.3969/j.issn.1005-3085.2022.03.006
    Abstract ( 93 )   PDF (648KB) ( 211 )   Save
    The weighted complementarity problem (WCP) aims at finding a pair of vectors belonging to the intersection of a manifold and a cone, such that the product of the vectors under a certain algebra equals a given weighted vector. As a nontrivial generalization of complementarity problems, WCP can be used to solve various equilibrium problems in science, economics and engineering, which in some cases may lead to highly efficient algorithms. Considering a weighted linear complementarity problem (WLCP) over the nonnegative orthant, an improved full-Newton step infeasible interior point algorithm is presented for its numerical solution. By extending a full-Newton infeasible interior-point algorithm for linear optimization, the perturbed problem of WLCP, its central path, and the induced Newton direction are introduced. The algorithm constructs strictly feasible iterates for a sequence of perturbed problems of WLCP. Each main iteration of the algorithm consists of one feasibility step and several central steps, which uses only full-Newton steps, and therefore it is not necessary to calculate the steplength; at each iteration, the algorithm reduces the feasibility residuals and the weight vector residuals at the same rate; based on the quadratic convergence result of the central step, a slightly wider neighborhood is provided for the feasibility step. The feasibility step, centering step and convergence are anaylzed. Then the algorithm is shown to possess global convergence and polynomial-time complexity. Finally, some numerical examples illustrate the efficiency of the proposed algorithm for solving WLCP.
    Related Articles | Metrics
    Maximum Norm a Posteriori Error Estimation for Nonlinear Singularly Perturbed Problems Based on a Hybrid Finite Difference Scheme
    BAO Xiaobing, LIU Libin, LIANG Zhifang
    2022, 39 (3):  428-438.  doi: 10.3969/j.issn.1005-3085.2022.03.007
    Abstract ( 69 )   PDF (249KB) ( 274 )   Save
    The adaptive moving mesh algorithm plays a very important role in the numerical solution of singularly perturbed differential equations. The key technology here is the construction of an effective discrete scheme and the corresponding a posteriori error estimation. Based on this, for a class of nonlinear singularly parameterized problems, the stability estimates of continuous solutions and related corollaries are given. Then, a hybrid finite difference scheme is established by using the backward Euler formula and the first-order central finite difference scheme on an arbitrary nonuniform grid, and the stability of the discrete solution is analyzed. Based on this stability estimation and the piecewise linear interpolation technique, an a posterior error estimation of the maximum norm of the mixed finite difference scheme is given. Using the a posterior error estimation, an optimal grid monitor function is selected, and an adaptive grid generation algorithm is designed based on the mesh equidistribution principle. Finally, numerical experiments verify the effectiveness of the adaptive moving mesh algorithm, and the average convergence order of the algorithm can reach the second order. Furthermore, it is shown from the numerical results that the error of the adaptive moving mesh is obviously smaller than that of the Shishkin mesh, and its convergence order is higher than that of the Shishkin mesh.
    Related Articles | Metrics
    The Control and Eigenvalue Analysis for a Pendulum System with Multiple Time-delay
    ZHAO Dongxia, FAN Dongxia, WANG Tingting, MAO Li
    2022, 39 (3):  439-450.  doi: 10.3969/j.issn.1005-3085.2022.03.008
    Abstract ( 87 )   PDF (287KB) ( 291 )   Save
    The stabilization of a pendulum system with position feedback and delayed position feedback controller is studied. Considering the fact that the controller itself has delay, a second-order differential equation with  two delays is established. We present the characteristic equation, obtain the relation between the system parameters and stability, and give the results for delay-dependent and delay-independent stability. By the technique for finding roots of algebra equation, the multiplicity of characteristic roots is analyzed, the parameter condition that the multiplicity is at most 4 is obtained, and it is proved that there are at most two triple characteristic roots. When the modulus of the characteristic root tends to infinity, we prove that the real part must tend to negative infinity. By using the asymptotic analysis technique, the asymptotic expression of the characteristic root is calculated. The validity of the conclusion is verified by Matlab numerical simulations. This method can be extended to $n$-order differential equations with multiple time-delays.
    Related Articles | Metrics
    Global Attractivity of a Competition and Cooperation Model of Enterprises with Two Time Delays
    MUHAMMADHAJI Ahmadjan
    2022, 39 (3):  451-462.  doi: 10.3969/j.issn.1005-3085.2022.03.009
    Abstract ( 71 )   PDF (2273KB) ( 212 )   Save
    This paper studies a class of non-autonomous competition and cooperation model of enterprises with two time delays. Sufficient conditions for the ultimately boundedness and global attractivity of the model are derived by using inequality techniques and the construction of multiple Lyapunov functions. As an application of the main results, also studied are three special cases of the considered model, and obtained the sufficient conditions on the ultimately boundedness and global attractivity of the considered models. Finally, two numerical examples further verified the validity and effectiveness of the obtained results. From the main results and two numerical examples, it can be concluded that the two discrete delays affect the global attractiveness of the model. The obtained results improvemed and extended the previously related studies.
    Related Articles | Metrics
    Bifurcation Analysis of a SIRS Age-structured Model with General Incidence Function
    ZHANG Suxia, LIU Yanna, XU Xiaxia
    2022, 39 (3):  463-476.  doi: 10.3969/j.issn.1005-3085.2022.03.010
    Abstract ( 84 )   PDF (1068KB) ( 140 )   Save
    Considering the impact of age in the transmission of infectious diseases, we present an age-structured model with SIRS type and general incidence function. By reformulating the model as an abstract Cauchy problem and applying theorems related with the Hille-Yosida operator, we investigated the dynamic properties, including stability of equilibria and the condition for Hopf bifurcation due to the destabilization of endemic equilibrium. The results reveal that, if the basic reproductive number is less than 1, the infection-free equilibrium is globally stable, without being influenced by the immune age. Conversely, if the basic reproductive number is larger than 1, the endemic equilibrium may be destabilized by the perturbation of the immune age and a Hopf bifurcation can occur. Meanwhile, numerical simulations are conducted to ill-ustrate the theoretical results and to show the influence of the immune age on the dynamical behaviors of the model.
    Related Articles | Metrics
    On the 2-restricted Edge Connectivity of Folded Crossed Cubes
    CAI Xuepeng, FAN Dandan, XU Ganggang
    2022, 39 (3):  477-486.  doi: 10.3969/j.issn.1005-3085.2022.03.011
    Abstract ( 76 )   PDF (182KB) ( 109 )   Save
    The $h$-restricted edge connectivity is an important parameter in measuring the reliability and fault tolerance of large interconnection networks. Let $G$ be a connected graph and $h$ be a non-negative integer. The $h$-restricted edge connectivity of $G$ is the minimum cardinality of a set of edges, if it exists, whose deletion disconnects $G$ and the degree of each vertex in every remaining component is at least $h$. The $n$-dimensional folded crossed cube is obtained from the $n$-dimensional crossed cube by adding extra edges. The $h$-restricted edge connectivity, which is an important measure in evaluating the reliability, is utilized to analyze the reliability of folded crossed cube. Then the $h$-restricted edge connectivity of folded crossed cubes is obtained. Finally, it is proved that the $2$-restricted edge connectivity of a folded crossed cube is equal to $4n-4 (n\geq 4)$. It means that at least $4n-4$ edges must be removed to disconnect a $n$-dimensional folded crossed cube, provided that the removal of these vertices does not leave a vertex that has degree less than two.
    Related Articles | Metrics
    Pricing of Equity Warrant Bond Based on the Theory of Backward Stochastic Differential Equations
    MIAO Jie
    2022, 39 (3):  487-494.  doi: 10.3969/j.issn.1005-3085.2022.03.012
    Abstract ( 76 )   PDF (195KB) ( 442 )   Save
    The backward stochastic differential equation theory builds a bridge between randomness and certainty, which makes it possible for using deterministic strategies to solve random and uncertain problems, and opens up a new way for the pricing of financial products. Thus, the pricing problem of equity warrant bond is studied by using the theory of backward stochastic differential equations. Firstly, under the no arbitrage assumption, the portfolio is reasonably established. Through the self-financing strategy, we derive a backward stochastic differential equation that the price of equity warrant bond satisfies by using the backward sto-chastic differential equation theory. Then, using the nonlinear Feynman-Kac formula, we get a partial differential equation that the price of equity warrant bond follows, and it is proved that the price of equity warrant bond at time 0 is equal to the conditional expectation of the maturity cash flow. The explicit formula for the price of an equity warrant bond is obtained by the martingale method. Finally, taking the equity warrant bond of Masteel as an example, the empirical analysis verifies that the pricing model obtained in this paper is reasonable.
    Related Articles | Metrics
    Stability and Hopf Bifurcation of a Reaction Diffusion Predator-prey Model with Prey Fear and Refuge
    CHEN Qingwan, LIU Wenqing
    2022, 39 (3):  495-501.  doi: 10.3969/j.issn.1005-3085.2022.03.013
    Abstract ( 81 )   PDF (2371KB) ( 231 )   Save
    Both fear factors and prey refuge have important effects in predation on the ecos-ystems. A class of reaction diffusion predator-prey models with fear effect and prey refuge is studied. We first provide the local asymptotic stability of the equilibrium point by using the linearization method and local bifurcation theory. Next, the existence of Hopf bifurcation and limit cycle is examined by choosing the ratio of un-protected prey as the bifurcation parameter. The results show that the existence of the refuge leads to Hopf bifurcation and spatially homogeneous periodic solution, and the addition of diffusion leads to new Hopf bifurcation points and inhomogeneous periodic solutions. This shows that the biological population can coexist by setting up an appropriate prey refuge or reducing the diffusion of predators. Finally, the conclusions are verified through numerical simulation.
    Related Articles | Metrics
    Solvability of Three-point Boundary Value Problems for $p$-Laplacian Equation under Resonant Conditions
    DUAN Lei, CHEN Tianlan
    2022, 39 (3):  502-510.  doi: 10.3969/j.issn.1005-3085.2022.03.014
    Abstract ( 103 )   PDF (170KB) ( 88 )   Save
    By using the solution set connectivity theory of the compact vector field equation and the method of well-ordered upper and lower solutions, studied in this paper is the three-point boundary value problem of $p$-Laplacian equation with the nonlinear term being continuous under resonant conditions. The problem is rewritten as a functional equation defined on Banach space, and the existence of solutions is obtained.
    Related Articles | Metrics