形变标定格值计算软件

为提高地壳形变观测仪器格值标定的真实性、准确性和可靠性，根据《地震前兆台网运行管理办法》和《地震前兆台网形变观测与运行工作细则》，针对不同测项、不同仪器标定的规范要求及特点，应用Visual Basic语言编制了格值计算软件。本文主要介绍了软件的使用、格值计算、误差检查方法，结果以全国形变台网中心推荐的Word格式保存。对水管仪、水平摆、垂直摆仪直接从原始文件读取数据，提高了标定精度，实现了格值标定的半自动化。本软件已应用于河北省地震前兆数字化台站工作。

论分段格值计算重力单程观测法公式

国产石英弹簧大测程重力仪分段标定的格值相对变化 >1/ 14 0 0时 ,文中总结分段格值计算单程观测 (混合零点位移改正、固体潮改正、静掉格改正 )的公式。

基于固体潮数采记录结果中标定的格值程序计算

针对固体潮汐监测仪器的数据采集结果,研究了其信号连续性判断的约束范围限方法,并利用它有效地区分了标定、地震和突跳等信号。用最小二乘法求多项式拟合方法得到了标定过程中固体潮背景信号的表达式,从而可计算出标定信号的偏移量大小,并得到仪器记录数采曲线的格值。

格值交替树自动机

交替(树)自动机因其本身关于取补运算的简洁性及其与非确定型(树)自动机的等价性,成为自动机与模型检测领域研究的一个新方向.在格值交替自动机与经典交替树自动机概念的基础上,引入格值交替树自动机的概念,并研究了格值交替树自动机的代数封闭性和表达能力.首先,证明了对格值交替树自动机的转移函数取对偶运算,终止权重取补之后所得自动机与原自动机接受语言互补这一结论.其次,证明了格值交替树自动机关于交、并运算的封闭性.最后,讨论了格值交替树自动机和格值树自动机、格值非确定型自动机的表达能力;证明了格值交替树自动机与格值树自动机的等价性,并给出了二者相互转化的算法及其复杂度分析;同时,提供了用格值非确定型自动机来模拟格值交替树自动机的方法.

High Resolution Ship Hydrodynamics Simulations in Open Source Environment

The numerical simulation of wake and flee-surface flow around ships is a complex topic that involves multiple tasks： the generation of an optimal computational grid and the development of numerical algorithms capable to predict the flow field around a hull. In this paper, a numerical framework is developed aimed at high-resolution CFD simulations of turbulent, free-surface flows around ship hulls. The framework consists in the concatenation of ＂tools＂, partly available in the open-source finite volume library OpenFOAM. A novel, flexible mesh-generation algorithm is presented, capable of producing high-quality computational grids for free-surface ship hydrodynamics. The numerical frame work is used to solve some benchmark problems, providing results that are in excellent agreement with the experimental measures.

Study of Defect-Layers Effect in Ferroelectric Thin Film with Transverse Ising Model

By taking into account the two-spin interaction in the transverse Ising model （TIM）, the influence of the defect layers （including JB and ΩB） on the polarization and Curie temperature are calculated numerically, within the framework of the decoupling approximation under Green＇s function. The numerical results show that the polarization and Curie temperature will both become large sensitively due to the large values of JB and the small value of ΩB of the defect layers. Meanwhile, the dependence of the crossover values of the exchange interaction JA, the transverse field ΩA of the bulk material on the exchange interaction JB and the transverse field ΩB of the defect layers are shown in 3-Dimensional （3-D） figures for the first time. Moreover, the transition features of the ferroelectric thin film with defect layers are presented.

First-principles Study of Electron Transport Through Oligoacenes

The electronic transport properties of oligoacenes sandwiched between two Au（111） surfaces with serial and parrallel configurations were investigeted by using a fully self-consistent nonequilibrium Green＇s function method combined with density functional calculations. This theoretical results show that the conductivity of oligoacenes with both sandwiched configurations at low bias voltage is mainly determined by the tail of the transmission peak from the perturbed highest occupied molecular orbital. When the molecular length increases, the zero-bias voltage conductance G（0） of oligoacenes with serial configuration neither follows Magoga＇s exponential law nor displays the even-odd oscillation effect, while the G（O） of the oligoacenes sandwiched with parallel configuration monotonically increases. The reduction of energy gaps, the alignment of the Fermi level, and the spatial distribution of the perturbed molecular orbitals are used to self-consistently explore the transport mechanism through oligoacenes.

Stochastic Lattice gas Cellular Automata Model for Epidemics

The aim of this study was to develop and explore a stochastic lattice gas cellular automata （LGCA） model for epidemics. A computer program was development in order to implement the model. An irregular grid of cells was used. A susceptible-infected-recovered （SIR） scheme was represented. Stochasticity was generated by Monte Carlo method. Dynamics of model was explored by numerical simulations. Model achieves to represent the typical SIR prevalence curve. Performed simulations also show how infection, mobility and distribution of infected individuals may influence the dynamics of propagation. This simple theoretical model might be a basis for developing more realistic designs.

Heat Conduction Analytical Solutions to be Applied in Boundary Conditions Obtained from Discrete Data

Analytical solutions have varied uses. One is to provide solutions that can be used in verification of numerical methods. Another is to provide relatively simple forms of exact solutions that can be used in estimating parameters, thus, it is possible to reduce computation time in comparison with numerical methods. In this paper, an alternative procedure is presented. Here is used a hybrid solution based on Green＇s function and real characteristics （discrete data） of the boundary conditions.

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.

Zero Energy Anomaly in One-Dimensional Anderson Lattice with Exponentially Correlated Weak Diagonal Disorder

We calculated numerically the localization length of one-dimensional Anderson model with correlated diagonal disorder. For zero energy point in the weak disorder limit, we showed that the localization length changes continuously as the correlation of the disorder increases. We found that higher order terms of the correlation must be included into the current perturbation result in order to give the correct localization length, arid to connect smoothly the anomaly at zero correlation with the perturbation result for large correlation.

Three-dimensional Numerical Simulation of the Movement of the Flexible Body under Different Constraints

For the large deformation of the flexible body may cause the fluid grid distortion,which will make the numerical calculation tedious,even to end,the numerical simulation of the flexible body coupling with the fluid is always a tough problem.In this paper,the flexible body is under two kinds of constrained conditions and the ratio of length-diameter is 1:30.The Reynolds number of the airflow is 513,belonging to the area of low Reynolds number.The control equations of the coupling of flexible body with airflow are built and the adaptive grid control method is adopted to conduct the three-dimensional numerical simulation of the movement of the flexible body.The numerical results show that it is possible to simulate the characteristics of the flexible body's movement in the low Reynolds number airflow when the appropriate control equations are modeled and suitable equation-solving method is adopted.Unconstrained flexible body would turn over forward along the airflow's diffusion direction,while constrained flexible body in the flow field will make periodic rotation motion along the axis of the flexible body,and the bending deformation is more obvious than that of unconstrained flexible body.The preliminary three-dimensional numerical simulation can provide references for further research on the characteristics of the yam movement in high Reynolds number airflow.