采油工程中灰色预测模型的病态性诊断

依据条件数法和特征分析法对采油工程中的灰色模型的病态性进行了分析,指出了两者在判定模型病态性方面的优劣,并根据最小二乘解析式分析了数乘变换对灰色模型参数的影响,文中还就Gershgorin定理分析了灰色模型病态的缘由。

GLOBAL CONVERGENCE OF A CLASS OF OPTIMALLY CONDITIONED SSVM METHODS

This paper explores the convergence of a class of optimally conditioned self scaling variable metric (OCSSVM) methods for unconstrained optimization. We show that this class of methods with Wolfe line search are globally convergent for general convex functions.

LOCALIZATION ALGORITHM BASED ON MINIMUM CONDITION NUMBER FOR WIRELESS SENSOR NETWORKS

During range-based self-localization of Wireless Sensor Network (WSN) nodes, the number and placement methods of beacon nodes have a great influence on the accuracy of localization. This paper proves a theorem which describes the relationship between the placement of beacon nodes and whether the node can be located in 3D indoor environment. In fact, as the highest locating accuracy can be acquired when the beacon nodes form one or more equilateral triangles in 2D plane, we generalizes this conclusion to 3D space, and proposes a beacon nodes selection algorithm based on the minimum condition number to get the higher locating accuracy, which can minimize the influence of distance measurement error. Simulation results show that the algorithm is effective and feasible.

Some Remarks On E＇-matrices

The class of E＇-matrices introduced in [4] is a subclass of Q0-matrices whose proper principal submatrices are Q--matrices. In this paper, we investigate the relation of the class E＇ to other known matrix classes and obtain a series of necessary and sufficient conditions for E＇-matrices.

An Improved Semi-Orthogonal User Selection Algorithm Based on Condition Number for Multiuser MIMO Systems

In Multiple-Input Multiple-Out （MIMO） systems, the user selection algorithm plays an important role in the realization of multiplexing gain. In this paper, an improved Semi-orthogonal User Selection algorithm based on condition number is proposed. Besides, a new MIMO pre- coding scheme is designed. The proposed SUS- CN （SUS with condition number） algorithm outperforms the SUS algorithm for the selection of users with better matrix inversion property, thus a higher information rate for selected user pair is achieved. The designed MIMO precoding matrix brings benefits of the power equality at transmitted terminals, the limited dynamic range of the power over time, and a better power efficiency. The simulation results give the key insights into the im- pact of the different condition number value and users on the sum-rate capacity.

Coverage Dependence of Growth Mode in Heteroepitaxy of Ni on Cu(100) Surface

A realistic kinetic Monte Carlo （KMC） simulation model with physical parameters is developed, which well reproduces the heteroepitaxial growth of multilayered Ni thin film on Cu（100） surfaces at room temperature. The effects of mass transport between interlayers and edge diffusion of atoms along the islands are included in the simulation model, and the surface roughness and the layer distribution versus total coverage are calculated. Speeially, the simulation model reveals the transition of growth mode with coverage and the difference between the Ni heteroepitaxy on Cu（100） and the Ni homoepitaxy on Ni（100）. Through comparison of KMC simulation with the real scanning tunneling microscopy （STM） experiments, the Ehrlich-Schwoebel （ES） barrier Ees is estimated to be 0.18±0.02 eV for Ni/Cu（100） system while 0.28 eV for Ni/Ni（100）. The simulation also shows that the growth mode depends strongly on the thickness of thin film and the surface temperature, and the critical thickness of growth mode transition is dependent on the growth condition such as surface temperature and deposition flux as well.

Simulation of Current Structure Interactions Using Least-Squares Finite Element Method

The objective of this paper is to investigate the condition number of various formulations of LSFEM （least-squares finite element method） for SWE （shallow-water equations）, and develop a better conditioned shallow-water model to simulate current structure interactions. Various formulations of LSFEM for a two-dimensional vertically-averaged non-viscous shallow-water equations can be constructed, depending on the choice of norm, variables, interpolations, and possible treatment of boundary conditions. The condition number of the resulting system of equations is systematically examined and compared. It is found that condition number of the resulting system of equations depends on the choice of variables, interpolations, and size of element （h）. Order reduction （UW） formulations, with introducing auxiliary variables, with low-order interpolation is better conditioned and more efficient than direct （U） formulation with high-order interpolation. However, to resolve large gradients and fine structures of flow filed, high-order methods are generally preferred. The developed shallow-water model is used to simulate flow past an elliptic hump and flow past a cylinder. Computed results are compared with other numerical solutions.

The Influence of Edge Restraining Stiffness on the Transverse Vibrations of Rectangular Plate Structures

This paper presents an analytical study on the influence of edge restraining stiffness on the transverse vibrations of rectangular plate structure. An improved Fourier series method was employed to analyze the transverse vibration of plate structure with general elastically restrained boundary conditions. A linear combination of a double Fourier series and eight auxiliary terms was sought as the admissible function of the flexural displacement of the plate, each term being a combination of a polynomial function and a single cosine series expansion. The auxiliary terms were introduced to ensure and improve the smoothness of the original displacement function and its derivatives at the boundaries. Several numerical examples were given to demonstrate the validity and accuracy of the current solution. The influences of translational and rotational stiffness on the natural frequencies and mode shapes of plate were analyzed by numerical results. The results show that the translational stiffness has bigger influence on the natural frequencies than the rotational stiffness. It is generally well known that little change of the rotational stiffness has little influence on the mode shapes of plate. However, the current work shows that a very little change of rotational stiffness value may lead to a large change of the mode shapes of a square plate structure.

Mixed-Weights Least-Squares Stable Predictive Control Algorithm with Soft and Hard Constraints

Mixed-weight least-squares (MWLS) predictive control algorithm, compared with quadratic programming (QP) method, has the advantages of reducing the computer burden, quick calculation speed and dealing with the case in which the optimization is infeasible. But it can only deal with soft constraints. In order to deal with hard constraints and guarantee feasibility, an improved algorithm is proposed by recalculating the setpoint according to the hard constraints before calculating the manipulated variable and MWLS algorithm is used to satisfy the requirement of soft constraints for the system with the input constraints and output constraints. The algorithm can not only guarantee stability of the system and zero steady state error, but also satisfy the hard constraints of input and output variables. The simulation results show the improved algorithm is feasible and effective.

Dynamics of Anchor Last Deployment of Submersible Buoy System

Submersible buoy systems are widely used for oceanographic research,ocean engineering and coastal defense.Severe sea environment has obvious effects on the dynamics of submersible buoy systems.Huge tension can occur and may cause the snap of cables,especially during the deployment period.This paper studies the deployment dynamics of submersible buoy systems with numerical and experimental methods.By applying the lumped mass approach,a three-dimensional multi-body model of submersible buoy system is developed considering the hydrodynamic force,tension force and impact force between components of submersible buoy system and seabed.Numerical integration method is used to solve the differential equations.The simulation output includes tension force,trajectory,profile and dropping location and impact force of submersible buoys.In addition,the deployment experiment of a simplified submersible buoy model was carried out.The profile and different nodes' velocities of the submersible buoy are obtained.By comparing the results of the two methods,it is found that the numerical model well simulates the actual process and conditions of the experiment.The simulation results agree well with the results of the experiment such as gravity anchor's location and velocities of different nodes of the submersible buoy.The study results will help to understand the conditions of submersible buoy's deployment,operation and recovery,and can be used to guide the design and optimization of the system.

Use of Wood Fibre Compost for the Cultivation of Trichoderma sp, （Isolate Td22）

The main objective of this research was to investigate the ability of a Trichoderma sp. （Td22）, inhibitory to Sclerotinia minor Jagger, to grow and survive in mature wood fibre waste （WFW） compost of paper mill origin following nutrient amendment. The growth and survival of the fungus in the WFW compost was assessed by serial dilution plate count method followed by confirmation of the fungal identity using pectic enzyme analysis as described in Cruickshank and Pitt [1]. It was found in this study that the population densities of TdE2 achieved under non-sterile conditions in the WFW compost following nutrient amendment was approximately in the range of 7.0 lOgl0 CFU/g dw - 8.5 log10 CFU/g dw after 28 days, depending on pre-treatment. The efficacy of this WFW compost-grown TdE2 for protection of lettuce from attack by S. minor was also demonstrated in glasshouse trials. This study indicates that cellulosic paper mill waste compost could provide an abundant low-cost growth medium for the large-scale cultivation of fungal antagonists, improving prospects for cost-competitiveness with chemical treatments.

Time domain system identification of unknown initial conditions

System identification is a method for using measured data to create or improve a mathematical model of the object being tested. From the measured data however, noise is noticed at the beginning of the response. One solution to avoid this noise problem is to skip the noisy data and then use the initial conditions as active parameters, to be found by using the system identification process. This paper describes the development of the equations for setting up the initial conditions as active parameters. The simulated data and response data from actual shear buildings were used to prove the accuracy of both the algorithm and the computer program, which include the initial conditions as active parameters. The numerical and experimental model analysis showed that the value of mass, stiffness and frequency were very reasonable and that the computed acceleration and measured acceleration matched very well.

BDS/GPS载波相位定位病态性对比分析

针对BDS/GPS载波相位定位的病态性问题,该文提出了采用条件数法对比分析BDS/GPS在我国不同地区、不同时段载波相位定位病态性及其在8个地面站的定位精度因子PDOP值变化情况,其中PDOP值能够反映单点定位的病态性。研究结果表明:BDS在我国大部分地区可用性良好,且其平均PDOP值随纬度的增大而线性增加;然而,BDS导航定位的病态性程度比GPS更为严重,BDS载波相位定位病态性程度不仅随观测时长变化,还与所处的位置变化有关,而GPS几乎不受观测地点的影响。研究成果可在BDS/GPS模糊度的快速、准确固定方面提供参考。

On the Existence for Some Special Primitive Elements in Finite Fields

Let F_q be a finite field of characteristic p. In this paper, by using the index sum method the authors obtain a sufficient condition for the existence of a primitive elementα∈ F_(q^n) such that α + α^(-1)is also primitive or α + α^(-1)is primitive and α is a normal element of F_(q^n) over F_q.

Effective condition numbers and small sample statistical condition estimation for the generalized Sylvester equation

Abstract In this paper, we investigate the effective condition numbers for the generalized Sylvester equation （AX - YB, DX - YE） = （C,F）, where A,D ∈ Rm×m B,E ∈ Rn×n and C,F ∈ Rm×n. We apply the small sample statistical method for the fast condition estimation of the generalized Sylvester equation, which requires （9（m2n ＋ mn2） flops, comparing with （-O（m3 ＋ n3） flops for the generalized Schur and generalized Hessenberg- Schur methods for solving the generalized Sylvester equation. Numerical examples illustrate the sharpness of our perturbation bounds.

Scaled boundary isogeometric analysis for 2D elastostatics

A new numerical method,scaled boundary isogeometric analysis(SBIGA)combining the concept of the scaled boundary finite element method(SBFEM)and the isogeometric analysis(IGA),is proposed in this study for 2D elastostatic problems with both homogenous and inhomogeneous essential boundary conditions.Scaled boundary isogeometric transformation is established at a specified scaling center with boundary isogeometric representation identical to the design model imported from CAD system,which can be automatically refined without communication with the original system and keeping geometry invariability.The field variable,that is,displacement,is constructed by the same basis as boundary isogeometric description keeping analytical features in radial direction.A Lagrange multiplier scheme is suggested to impose the inhomogeneous essential boundary conditions.The new proposed method holds the semi-analytical feature inherited from SBFEM,that is,discretization only on boundaries rather than the entire domain,and isogeometric boundary geometry from IGA,which further increases the accuracy of the solution.Numerical examples,including circular cavity in full plane,Timoshenko beam with inhomogenous boundary conditions and infinite plate with circular hole subjected to remotely tension,demonstrate that SBIGA can be applied efficiently to elastostatic problems with various boundary conditions,and powerful in accuracy of solution and less degrees of freedom(DOF)can be achieved in SBIGA than other methods.

Analysis of thin plates by the weak form quadrature element method

The recently proposed weak form quadrature element method （QEM） is applied to flexural and vibrational analysis of thin plates The integrals involved in the variational description of a thin plate are evaluated by an efficient numerical scheme and the par- tial derivatives at the integration sampling points are then approximated using differential quadrature analogs. Neither the grid pattern nor the number of nodes is fixed, being adjustable according to convergence need. The C~ continuity conditions char- acterizing the thin plate theory are discussed and the robustness of the weak form quadrature element for thin plates against shape distortion is examined. Examples are presented and comparisons with analytical solutions and the results of the finite element method are made to demonstrate the convergence and computational efficiency of the weak form quadrature element method. It is shown that the present formulation is applicable to thin plates with varying thickness as well as uniform plates.

Assessment of alternative scale-providing variables in a Reynolds-stress model using high-order methods

This paper is set in the high-order finite-difference discretization of the Reynolds-averaged Navier-Stokes(RANS)equations,which are coupled with the turbulence model equations.Three alternative scale-providing variables for the specific dissipation rate(o)are implemented in the framework of the Reynolds stress model(RSM)for improving its robustness.Specifically,g(=1/√ω)has natural boundary conditions and reduced spatial gradients,and a new numerical constraint is imposed on itω(=lnω)can preserve positivity and also has reduced spatial gradients;the eddy viscosity v,also has natural boundary conditions and its equation is improved in this work.The solution polynomials of the mean-flow and turbulence-model equations are both reconstructed by the weighted compact nonlinear scheme(WCNS).Moreover,several numerical techniques are introduced to improve the numerical stability of the equation system.A range of canonical as well as industrial turbulent flows are simulated to assess the accuracy and robustness of the scale-transformed models.Numerical results show that the scale-transformed models have significantly improved robustness compared to the w model and still keep the characteristics of RSM.Therefore,the high-order discretization of the RANS and RSM equations,which number 12 in total,can be successfully achieved.

A NEW PRINCIPAL PIVOTING SCHEME FOR BOX LINEAR COMPLEMENTARITY PROBLEMS

Judice and Pires developed in recent years principal pivoting methods for the solving of the so called box linear complementarity problems (BLCPs) where the constraint matrices are restrictedly supposed to be of P matrices. This paper aims at presenting a new principal pivoting scheme for BLCPs where the constraint matrices are loosely supposed to be row sufficient.This scheme can be applied to the solving of convex quadratic programs subject to linear constraints and arbitrary upper and lower bound constraints on variables.