Novel damage constitutive models and new quantitative identification method for stress thresholds of rocks under uniaxial compression

Four key stress thresholds exist in the compression process of rocks,i.e.,crack closure stress(σ_(cc)),crack initiation stress(σ_(ci)),crack damage stress(σ_(cd))and compressive strength(σ_(c)).The quantitative identifications of the first three stress thresholds are of great significance for characterizing the microcrack growth and damage evolution of rocks under compression.In this paper,a new method based on damage constitutive model is proposed to quantitatively measure the stress thresholds of rocks.Firstly,two different damage constitutive models were constructed based on acoustic emission(AE)counts and Weibull distribution function considering the compaction stages of the rock and the bearing capacity of the damage element.Then,the accumulative AE counts method(ACLM),AE count rate method(CRM)and constitutive model method(CMM)were introduced to determine the stress thresholds of rocks.Finally,the stress thresholds of 9 different rocks were identified by ACLM,CRM,and CMM.The results show that the theoretical stress−strain curves obtained from the two damage constitutive models are in good agreement with that of the experimental data,and the differences between the two damage constitutive models mainly come from the evolutionary differences of the damage variables.The results of the stress thresholds identified by the CMM are in good agreement with those identified by the AE methods,i.e.,ACLM and CRM.Therefore,the proposed CMM can be used to determine the stress thresholds of rocks.

AUTOMATIC MULTILEVEL THRESHOLDING METHOD BASED ON MAXIMUM ENTROPY

In the multilevel thresholding segmentation of the image, the classification number is always given by the supervisor. To solve this problem, a fast multilevel thresholding algorithm considering both the threshold value and the classification number is proposed based on the maximum entropy, and the self-adaptive criterion of the classification number is given. The algorithm can obtain thresholds and automatically decide the classification number. Experimental results show that the algorithm is effective.

NOVEL IMMERSED BOUNDARY-LATTICE BOLTZMANN METHOD BASED ON FEEDBACK LAW

The lattice Boltzmann method （LBM） and the immersed boundary method （IBM） are alternative, com- putational techniques for solving complex fluid dynamics systems, and can take the place of the Navier-Stokes（N- S） equation. This paper proposes a novel immersed boundary-lattice Boltzmann method （IB-LBM） based on the feedback law. The method uses the immersed boundary concept in the LBM framework to capture the coupling between a body with complex geometry and a uniform fluid, Then, the flows around a stationary circular cylinder and two circular cylinders in a side by side arrangement are simulated by using the method. Results are agreed well with the benchmark data, so, the capability of the method for complex geometry is demonstrated. Different from the conventional IB-LBM, which uses the Hook＇s law or the direct forcing method to compute the interae- tion force, the method uses the feedback law--the feedback of velocity field and displacement information to calculate the force, thus ensuring the method has advantages of easy implementation and full parallelism.

A local pseudo arc-length method for hyperbolic conservation laws

A local pseudo arc-length method（LPALM）for solving hyperbolic conservation laws is presented in this paper.The key idea of this method comes from the original arc-length method,through which the critical points are bypassed by transforming the computational space.The method is based on local changes of physical variables to choose the discontinuous stencil and introduce the pseudo arc-length parameter,and then transform the governing equations from physical space to arc-length space.In order to solve these equations in arc-length coordinate,it is necessary to combine the velocity of mesh points in the moving mesh method,and then convert the physical variable in arclength space back to physical space.Numerical examples have proved the effectiveness and generality of the new approach for linear equation,nonlinear equation and system of equations with discontinuous initial values.Non-oscillation solution can be obtained by adjusting the parameter and the mesh refinement number for problems containing both shock and rarefaction waves.

Identifying influential spreaders in complex networks based on entropy weight method and gravity law

In complex networks,identifying influential spreader is of great significance for improving the reliability of networks and ensuring the safe and effective operation of networks.Nowadays,it is widely used in power networks,aviation networks,computer networks,and social networks,and so on.Traditional centrality methods mainly include degree centrality,closeness centrality,betweenness centrality,eigenvector centrality,k-shell,etc.However,single centrality method is onesided and inaccurate,and sometimes many nodes have the same centrality value,namely the same ranking result,which makes it difficult to distinguish between nodes.According to several classical methods of identifying influential nodes,in this paper we propose a novel method that is more full-scaled and universally applicable.Taken into account in this method are several aspects of node’s properties,including local topological characteristics,central location of nodes,propagation characteristics,and properties of neighbor nodes.In view of the idea of the multi-attribute decision-making,we regard the basic centrality method as node’s attribute and use the entropy weight method to weigh different attributes,and obtain node’s combined centrality.Then,the combined centrality is applied to the gravity law to comprehensively identify influential nodes in networks.Finally,the classical susceptible-infected-recovered(SIR)model is used to simulate the epidemic spreading in six real-society networks.Our proposed method not only considers the four topological properties of nodes,but also emphasizes the influence of neighbor nodes from the aspect of gravity.It is proved that the new method can effectively overcome the disadvantages of single centrality method and increase the accuracy of identifying influential nodes,which is of great significance for monitoring and controlling the complex networks.

An Improved Double-Threshold Method Based on Gradient Histogram

This paper analyzes the characteristics of the output gradient histogram and shortages of several traditional automatic threshold methods in order to segment the gradient image better. Then an improved double-threshold method is proposed, which is combined with the method of maximum classes variance, estimating-area method and double-threshold method. This method can automatically select two different thresholds to segment gradient images. The computer simulation is performed on the traditional methods and this algorithm and proves that this method can get satisfying result. Key words gradient histogram image - threshold selection - double-threshold method - maximum classes variance method CLC number TP 391. 41 Foundation item: Supported by the National Nature Science Foundation of China (50099620) and the Project of Chenguang Plan in Wuhan (985003062)Biography: YANG Shen (1977-), female, Ph. D. candidate, research direction: multimedia information processing and network technology.

Study on the critical stress threshold of weakly cemented sandstone damage based on the renormalization group method

During the microstructural analysis of weakly cemented sandstone,the granule components and ductile structural parts of the sandstone are typically generalized.Considering the contact between granules in the microstructure of weakly cemented sandstone,three basic units can be determined:regular tetrahedra,regular hexahedra,and regular octahedra.Renormalization group models with granule-and pore-centered weakly cemented sandstone were established,and,according to the renormalization group transformation rule,the critical stress threshold of damage was calculated.The results show that the renormalization model using regular octahedra as the basic units has the highest critical stress threshold.The threshold obtained by iterative calculations of the granule-centered model is smaller than that obtained by the pore-centered model.The granule-centered calculation provides the lower limit(18.12%),and the pore-centered model provides the upper limit(36.36%).Within this range,the weakly cemented sandstone is in a phase-like critical state.That is,the state of granule aggregation transforms from continuous to discrete.In the relative stress range of 18.12%-36.36%,the weakly cemented sandstone exhibits an increased proportion of high-frequency signals(by 83.3%)and a decreased proportion of low-frequency signals(by 23.6%).The renormalization calculation results for weakly cemented sandstone explain the high-low frequency conversion of acoustic emission signals during loading.The research reported in this paper has important significance for elucidating the damage mechanism of weakly cemented sandstone.

Mammogram Images Thresholding for Breast Cancer Detection Using Different Thresholding Methods

The purpose of this study is to apply different thresholding in mammogram images, and then we will determine which technique is the best in thresholding (extraction) malignant and benign tumors from the rest breast tissues. The used technique is Otsu method, because it is one of the most effective methods for most real world views with regard to uniformity and shape measures. Also, we present all the thresholding methods that used the concept of between class variance. We found from the experimental results that all the used thresholding techniques work well in detection normal breast tissues. But in abnormal tissues (breast tumors), we found that only neighborhood valley emphasis method gave best detection of malignant tumors. Also, the results demonstrate that variance and intensity contrast technique is the best in extraction the micro calcifications which represent the first signs of breast cancer.

Dual threshold search method for asperity boundary determination based on geodetic and seismic catalog data

As an important model for explaining the seismic rupture mode,the asperity model plays an important role in studying the stress accumulation of faults and the location of earthquake initiation.Taking Qilian-Haiyuan fault as an example,this paper combines geodetic method and b-value method to propose a multi-source observation data fusion detection method that accurately determines the asperity boundary named dual threshold search method.The method is based on the criterion that the b-value asperity boundary should be most consistent with the slip deficit rate asperity boundary.Then the optimal threshold combination of slip deficit rate and b-value is obtained through threshold search,which can be used to determine the boundary of the asperity.Based on this method,the study finds that there are four potential asperities on the Qilian-Haiyuan fault:two asperities(A1 and A2)are on the Tuolaishan segment and the other two asperities(B and C)are on Lenglongling segment and Jinqianghe segment,respectively.Among them,the lengths of asperities A1 and A2 on Tuolaishan segment are 17.0 km and 64.8 km,respectively.And the lower boundaries are 5.5 km and 15.5 km,respectively;The length of asperity B on Lenglongling segment is 70.7 km,and the lower boundary is 10.2 km.The length of asperity C on Jinqianghe segment is 42.3 km,and the lower boundary is 8.3 km.

A New Regularized Minimum Error Thresholding Method

To overcome the shortcoming that the traditional minimum error threshold method can obtain satisfactory image segmentation results only when the object and background of the image strictly obey a certain type of probability distribution,one proposes the regularized minimum error threshold method and treats the traditional minimum error threshold method as its special case.Then one constructs the discrete probability distribution by using the separation between segmentation threshold and the average gray-scale values of the object and background of the image so as to compute the information energy of the probability distribution.The impact of the regularized parameter selection on the optimal segmentation threshold of the regularized minimum error threshold method is investigated.To verify the effectiveness of the proposed regularized minimum error threshold method,one selects typical grey-scale images and performs segmentation tests.The segmentation results obtained by the regularized minimum error threshold method are compared with those obtained with the traditional minimum error threshold method.The segmentation results and their analysis show that the regularized minimum error threshold method is feasible and produces more satisfactory segmentation results than the minimum error threshold method.It does not exert much impact on object acquisition in case of the addition of a certain noise to an image.Therefore,the method can meet the requirements for extracting a real object in the noisy environment.

Extraction of LUCC with different methods and threshold value

The research of land use and land cover (LUCC) is an important aspect in the global change research. The goal of this study is to find methods of extraction of LUCC’s change outlined and change type from remotely sensed data. Take the country of Fengxian in Shanghai as an example, it was supposed two steps to finish extraction of LUCC information: the first step was to use different methods, which is used to outline change areas; the second step include methods of false composing of two-temporal and threshold value. Through combining two methods, a model rule is built and the LUCC product is obtained, four kinds of change type within the study area are given, and the results are obvious. Finally, the results support the application of the high resolution image and tasseled cap composition (greenness and wetness) in the specific regional too.

SOLUTION OF THE RAYLEIGH PROBLEM FOR A POWER-LAW NON-NEWTONIAN CONDUCTING FLUID VIA GROUP METHOD

An investigation is made of the magnetic Rayleigh problem where a semi_infinite plate is given an impulsive motion and thereafter moves with constant velocity in a non_Newtonian power law fluid of infinite extent. The solution of this highly non_linear problem is obtained by means of the transformation group theoretic approach. The one_parameter group transformation reduces the number of independent variables by one and the governing partial differential equation with the boundary conditions reduce to an ordinary differential equation with the appropriate boundary conditions. Effect of the some parameters on the velocity u(y,t) has been studied and the results are plotted.

SPECTRAL/HP ELEMENT METHOD WITH HIERARCHICAL RECONSTRUCTION FOR SOLVING NONLINEAR HYPERBOLIC CONSERVATION LAWS

The hierarchical reconstruction （HR） [Liu, Shu, Tadmor and Zhang, SINUM ＇07] has been successfully applied to prevent oscillations in solutions computed by finite volume, Runge-Kutta discontinuous Galerkin, spectral volume schemes for solving hyperbolic conservation laws. In this paper, we demonstrate that HR can also be combined with spectral/hp element method for solving hyperbolic conservation laws. An orthogonal spectral basis written in terms of Jacobi polynomials is applied. High computational efficiency is obtained due to such matrix-free algorithm. The formulation is conservative, and essential nomoscillation is enforced by the HR limiter. We show that HR preserves the order of accuracy of the spectral/hp element method for smooth solution problems and generate essentially non-oscillatory solutions profiles for capturing discontinuous solutions without local characteristic decomposition. In addition, we introduce a postprocessing technique to improve HR for limiting high degree numerical solutions.

An Information Entropy-Based Methodology to Construct the Avulsion Threshold in the Tail Reach of the Estuarine Alluvial Plain

Channel avulsion is a natural phenomenon that occurs abruptly on alluvial river deltas,which can affect the channel stability.The causes for avulsion could be generally categorized as topography-and flood-driven factors.However,previous studies on avulsion thresholds usually focused on topography-driven factors due to the centurial or millennial avulsion timescales of the world’s most deltas,but neglected the impacts of flood-driven factors.In the current study,a novel demarcation equation including the two driven factors was proposed,with the decadal timescale of avulsion being considered in the Yellow River Estuary(YRE).In order to quantify the contributions of different factors in each category,an entropy-based methodology was used to calculate the contributing weights of these factors.The factor with the highest weight in each category was then used to construct the demarcation equation,based on avulsion datasets associated with the YRE.An avulsion threshold was deduced according to the demarcation equation.This avulsion threshold was then applied to conduct the risk assessment of avulsion in the YRE.The results show that:two dominant factors cover respectively geomorphic coefficient representing the topography-driven factor and fluvial erosion intensity representing the flood-driven factor,which were thus employed to define a two dimensional mathematical space in which the demarcation equation can be obtained;the avulsion threshold derived from the equation was also applied in the risk assessment of avulsion;and the avulsion threshold proposed in this study is more accurate,as compared with the existing thresholds.

Godunov＇s method for initial-boundary value problem of scalar conservation laws

This paper is concerned with Godunov＇s scheme for the initial-boundary value problem of scalar conservation laws. A kind of Godunov＇s scheme, which satisfies the boundary entropy condition, was given. By use of the scheme, numerical simulation for the weak entropy solution to the initial-boundary value problem of scalar conservation laws is conducted.

Experimental determination of threshold pressure and permeability based on equation-solving method for liquid metal infiltration processes

A new method for determining two key parameters(threshold pressure and permeability)for fabricating metal matrix composites was proposed based on the equation-solving method.An infiltration experimental device was devised to measure the infiltration behavior precisely with controllable infiltration velocity.Two experiments with alloy Pb-Sn infiltrating into Al2O3 preform were conducted independently under two different pressures so as to get two different infiltration curves.Two sets of coefficients which are functions of threshold pressure and permeability can be obtained through curve fitting method.By solving the two-variable equation set,two unknown variables were determined.It is shown that the determined threshold pressure and permeability are very close to the calculated ones and are also verified by another independent infiltration experiment.The proposed method is also feasible to determine the key infiltration parameters for other metal matrix composite systems.

Conservative and Easily Implemented Finite Volume Semi-Lagrangian WENO Methods for 1D and 2D Hyperbolic Conservation Laws

The paper is devised to propose finite volume semi-Lagrange scheme for approximating linear and nonlinear hyperbolic conservation laws. Based on the idea of semi-Lagrangian scheme, we transform the integration of flux in time into the integration in space. Compared with the traditional semi-Lagrange scheme, the scheme devised here tries to directly evaluate the average fluxes along cell edges. It is this difference that makes the scheme in this paper simple to implement and easily extend to nonlinear cases. The procedure of evaluation of the average fluxes only depends on the high-order spatial interpolation. Hence the scheme can be implemented as long as the spatial interpolation is available, and no additional temporal discretization is needed. In this paper, the high-order spatial discretization is chosen to be the classical 5th-order weighted essentially non-oscillatory spatial interpolation. In the end, 1D and 2D numerical results show that this method is rather robust. In addition, to exhibit the numerical resolution and efficiency of the proposed scheme, the numerical solutions of the classical 5th-order WENO scheme combined with the 3rd-order Runge-Kutta temporal discretization (WENOJS) are chosen as the reference. We find that the scheme proposed in the paper generates comparable solutions with that of WENOJS, but with less CPU time.

Numerical Method for Non-Linear Conservation Laws: Inviscid Burgers Equation

This paper deals with the Burgers equation which is the most common model used in the nonlinear conservation laws. Here the theoretical aspect of conservation law is discussed by using inviscid Burgers equation. At first, we introduce the general non-linear conservation law as a partial differential equation and its solution procedure by the method of characteristic. Next, we present the weak solution of the problem with entropy condition. Taking into account shock wave and rarefaction wave, the Riemann problem has also been discussed. Finally, the finite volume method is considered to approximate the numerical solution of the inviscid Burgers equation with continuous and discontinuous initial data. An illustration of the problem is provided by some examples. Moreover, the Godunov method provides a good approximation for the problem.

A Riemann-Solver Free Spacetime Discontinuous Galerkin Method for General Conservation Laws

This paper summarizes a Riemann-solver-free spacetime discontinuous Galerkin method developed for general conservation laws. The method integrates the best features of the spacetime Conservation Element/Solution Element (CE/SE) method and the discontinuous Galerkin (DG) method. The core idea is to construct a staggered spacetime mesh through alternate cell-centered CEs and vertex-centered CEs within each time step. Inside each SE, the solution is approximated using high-order spacetime DG basis polynomials. The spacetime flux conservation is enforced inside each CE using the DG concept. The unknowns are stored at both vertices and cell centroids of the spatial mesh. However, the solutions at vertices and cell centroids are updated at different time levels within each time step in an alternate fashion. Thanks to the staggered spacetime formulation, there are no left and right states for the solution at the spacetime interface. Instead, the solution available to evaluate the flux is continuous across the interface. Therefore, no (approximate) Riemann solvers are needed to provide a unique numerical flux. The current method can be used to solve arbitrary conservation laws including the compressible Euler equations, shallow water equations and magnetohydrodynamics (MHD) equations without the need of any form of Riemann solvers. A set of benchmark problems of various conservation laws are presented to demonstrate the accuracy of the method.

Research on the Method of Balance Sheet Quality Inspection based on the Ford Law

This paper based on Benford＇s Law （Benford Law）, the Shanghai stock exchange 15 listed bank balance sheets for the past 5 years data for the study sample, results show that the model in the inspection report can quantify the quality in the process of accounting information quality comparison and has the outstanding advantages of simple, clear, this paper focuses on the Benford＇s law of balance sheet quality inspection method based on in our country and how to test the quality of balance sheets in the introduction and promotion.