Smoothing estimation of rate function for recurrent event data with informative censoring

This paper proposes kernel estimation of the occurrence rate function for recurrent event data with informative censoring. An informative censoring model is considered with assumptions made on the joint distribution of the recurrent event process and the censoring time without modeling the censoring distribution. Under the validity of the informative censoring model, we also show that an estimator based on the assumption of independent censoring becomes inappropriate and is generally asymptotically biased. To investigate the asymptotic properties of the proposed estimator, the explicit form of its asymptotic mean squared risk and the asymptotic normality are derived. Meanwhile, the empirical consistent smoothing estimator for the variance function of the estimator is suggested. The performance of the estimators are also studied through Monte Carlo simulations. An epidemiological example of intravenous drug user data is used to show the influence of informative censoring in the estimation of the occurrence rate functions for inpatient cares over time.     

Large scale petroleum reservoir simulation and parallel preconditioning algorithms research

Solving large scale linear systems efficiently plays an important role in a petroleum reservoir simulator, and the key part is how to choose an effective parallel preconditioner. Properly choosing a good preconditioner has been beyond the pure algebraic field. An integrated preconditioner should include such components as physical background, characteristics of PDE mathematical model, nonlinear solving method, linear     

非线性Schrdinger方程的高精度守恒差分格式

本文首先分析线性Schrodinger方程一种高阶差分格式的构造方法,得到方程的耗散项.在此基础上对三次非线性Schrodinger方程,提出了一种精度为O(r2 h2)的差分格式,证明了该格式保持了连续方程的两个守恒量,且是收敛的与稳定的.并通过数值例子与已有隐格式进行了比较,结果表明,本文格式在计算量类似的情况下,提高了数值精度.     

Bounds to optimal burn‐in and optimal work size

Burn‐in is a widely used method to improve the quality of products or systems after they have been produced. In this paper, we consider the problem of determining the optimal burn‐in time and optimal work size maximizing the long‐run average amount of work saved per time unit in the computer applications. Assuming that the underlying lifetime distribution of the computer has an initially decreasing or/and eventually increasing failure rate function, an upper bound for the optimal burn‐in time is derived for each fixed work size and a uniform (with respect to the burn‐in time) upper bound for the optimal work size is also obtained. Furthermore, it is shown that a non‐trivial lower bound for the optimal burn‐in time can be derived if the underlying lifetime distribution has a large initial failure rate. Copyright © 2005 John Wiley & Sons, Ltd.     

190区超形变核中转动带的组态结构

通过系统研究A～190区超形变核中转动带的转动惯量、角动量顺排、旋称分裂随转动频率的变化规律, 结合我们用处理对力的粒子数守恒方法的计算结果, 对A～190区所有转动带的组态结构给出了一个整体的描述. 绝大多数超形变带都建立在强耦合轨道上, 例如中子[512]5/2, [624]9/2. 少数超形变带则建立在高j闯入轨道上, 即中子[761]3/2, [752]5/2. 根据我们提出的组态结构所进行的理论计算结果表明, A～190区所有转动带的一般行为、反常变化和带交叉都得到了满意的解释.     

从像上提取非旋转对称系统的传递函数制作解卷空间滤波器

本文提出从高分辨电子显微像中,利用光学、全息和部分相干等方法,提取传递函数,制作相位和振幅滤波器,实现解卷,提高高分辨电子显微像的分辨率.     

Direct numerical simulation of settling of a large solid particle during bioconvection

Settling of a large solid particle in bioconvection flow caused by gyrotactic microorganisms is investigated. The particle is released from the top of the bioconvection chamber; its settling pattern depends on whether it is released in the centre of the bioconvection plume or at its periphery. The Chimera method is utilized; a subgrid is generated around a moving particle. The method suggested by Liu and Wang (Comput. Fluid 2004; 33 :223–255) is further developed to account for the presence of a moving boundary in the streamfunction‐vorticity formulation using the finite‐difference method. A number of cases for different release positions of the particle are computed. It is demonstrated that bioconvection can either accelerate or decelerate settling of the particle depending on the initial position of the particle relative to the plume centre. It is also shown that the particle impacts bioconvection plume by changing its shape and location in the chamber. Copyright © 2005 John Wiley & Sons, Ltd.     

Computations of two passing‐by high‐speed trains by a relaxation overset‐grid algorithm

This paper presents a relaxation algorithm, which is based on the overset grid technology, an unsteady three‐dimensional Navier–Stokes flow solver, and an inner‐ and outer‐relaxation method, for simulation of the unsteady flows of moving high‐speed trains. The flow solutions on the overlapped grids can be accurately updated by introducing a grid tracking technique and the inner‐ and outer‐relaxation method. To evaluate the capability and solution accuracy of the present algorithm, the computational static pressure distribution of a single stationary TGV high‐speed train inside a long tunnel is investigated numerically, and is compared with the experimental data from low‐speed wind tunnel test. Further, the unsteady flows of two TGV high‐speed trains passing by each other inside a long tunnel and at the tunnel entrance are simulated. A series of time histories of pressure distributions and aerodynamic loads acting on the train and tunnel surfaces are depicted for detailed discussions. Copyright © 2004 John Wiley & Sons, Ltd.     

Assessment of two‐equation turbulence modelling for high Reynolds number hydrofoil flows

This paper presents an evaluation of the capability of turbulence models available in the commercial CFD code FLUENT 6.0 for their application to hydrofoil turbulent boundary layer separation flow at high Reynolds numbers. Four widely applied two‐equation RANS turbulence models were assessed through comparison with experimental data at Reynolds numbers of 8.284×106 and 1.657×107. They were the standard k–εmodel, the realizable k–εmodel, the standard k–ωmodel and the shear‐stress‐transport (SST) k–ωmodel. It has found that the realizable k–εturbulence model used with enhanced wall functions and near‐wall modelling techniques, consistently provides superior performance in predicting the flow characteristics around the hydrofoil. Copyright © 2004 John Wiley & Sons, Ltd.     

A virtual boundary method with improved computational efficiency using a multi‐grid method

The flow around spherical, solid objects is considered. The boundary conditions on the solid boundaries have been applied by replacing the boundary with a surface force distribution on the surface, such that the required boundary conditions are satisfied. The velocity on the boundary is determined by extrapolation from the flow field. The source terms are determined iteratively, as part of the solution. They are then averaged and are smoothed out to nearby computational grid points. A multi‐grid scheme has been used to enhance the computational efficiency of the solution of the force equations. The method has been evaluated for flow around both moving and stationary spherical objects at very low and intermediate Reynolds numbers. The results shows a second order accuracy of the method both at creeping flow and at Re=100. The multi‐grid scheme is shown to enhance the convergence rate up to a factor 10 as compared to single grid approach. Copyright © 2004 John Wiley & Sons, Ltd.