非均质材料损伤破坏的数值模拟软件开发及强韧化机理研究

本文针对平面链网模型与空间链网模型，发展了模拟非均质材料的数值软件，旨在分析材料的损伤破坏特征，寻找材料细观结构同宏观力学行为之间的相互联系，探讨材料强韧化机理以及评估材料的力学性能。采用可视化技术设计的具有交互式功能的图象用户界面，能根据具体材料对象构造特定材料模型来分析其力学行为与服役能力。初步结果表明，无论是材料弹性阶段的性质还是损伤演化阶段材料中裂纹的萌生、扩展、汇集和材料的最终失效等现象都与已有的理论结果及试验结果符合得很好。另外，本文还对材料的强韧化机理作了探索。

Color to Gray: Visual Cue Preservation

Both commercial and scientific applications often need to transform color images into gray-scale images, e. g., to reduce the publication cost in printing color images or to help color blind people see visual cues of color images. However, conventional color to gray algorithms are not ready for practical applications because they encounter the following problems: 1) Visual cues are not well defined so it is unclear how to preserve important cues in the transformed gray-scale images; 2) some algorithms have extremely high time cost for computation; and 3) some require human-computer interactions to have a reasonable transformation. To solve or at least reduce these problems, we propose a new algorithm based on a probabilistic graphical model with the assumption that the image is defined over a Markov random field. Thus, color to gray procedure can be regarded as a labeling process to preserve the newly well-defined visual cues of a color image in the transformed gray-scale image. Visual cues are measurements that can be extracted from a color image by a perceiver. They indicate the state of some properties of the image that the perceiver is interested in perceiving. Different people may perceive different cues from the same color image and three cues are defined in this paper, namely, color spatial consistency, image structure information, and color channel perception priority. We cast color to gray as a visual cue preservation procedure based on a probabilistic graphical model and optimize the model based on an integral minimization problem. We apply the new algorithm to both natural color images and artificial pictures, and demonstrate that the proposed approach outperforms representative conventional algorithms in terms of effectiveness and efficiency. In addition, it requires no human-computer interactions.

Computation of Viscous Incompressible Flow using Pressure Correction Method on Unstructured Chimera Grid

In this paper, a pressure correction algorithm for computing incompressible flows is modified and implemented on unstructured Chimera grid. Schwarz method is used to couple the solutions of different sub-domains. A new interpolation to ensure consistency between primary variables and auxiliary variables is proposed. Other important issues such as global mass conservation and order of accuracy in the interpolations are also discussed. Two numerical simulations are successfully performed. They include one steady case, the lid-driven cavity and one unsteady case, the flow around a circular cylinder. The results demonstrate a very good performance of the proposed scheme on unstructured Chimera grids. It prevents the decoupling of pressure field in the overlapping region and requires only little modification to the existing unstructured Navier–Stokes (NS) solver. The numerical experiments show the reliability and potential of this method in applying to practical problems.

A time domain computation method for dynamic behavior of mooring system

A quasi-steady time domain method is developed for the prediction of dynamic behavior of a mooring system under the environmental disturbances, such as regular or irregular waves, winds and currents. The mooring forces are obtained in a static sense at each instant. The dynamic feature of the mooring cables can be obtained by incorporating the extended 3-D lumped-mass method with the known ship motion history. Some nonlinear effects, such as the influence of the instantaneous change of the wetted hull surface on the hydrostatic restoring forces and Froude-Krylov forces, are included. The computational results show a satisfactory agreement with the experimental ones.

Semi-weight function method on computation of stress intensity factors in dissimilar materials

Semi-weight function method is developed to solve the plane problem of two bonded dissimilar materials containing a crack along the bond. From equilibrium equation, stress and strain relationship, conditions of continuity across interface and free crack surface, the stress and displacement fields were obtained. The eigenvalue of these fields is lambda. Semi-weight functions were obtained as virtual displacement and stress fields with eigenvalue-lambda. Integral expression of fracture parameters, K-I and K-II, were obtained from reciprocal work theorem with semi-weight functions and approximate displacement and stress values on any integral path around crack tip. The calculation results of applications show that the semi-weight function method is a simple, convenient and high precision calculation method.

On the computation of space-time correlations by large-eddy simulation

The effect of subgrid-scale (SGS) modeling on velocity (space-) time correlations is investigated in decaying isotropic turbulence. The performance of several SGS models is evaluated, which shows superiority of the dynamic Smagorinsky model used in conjunction with the multiscale large-eddy simulation (LES) procedure. Compared to the results of direct numerical simulation, LES is shown to underpredict the (un-normalized) correlation magnitude and slightly overpredict the decorrelation time scales. This can lead to inaccurate solutions in applications such as aeroacoustics. The underprediction of correlation functions is particularly severe for higher wavenumber modes which are swept by the most energetic modes. The classic sweeping hypothesis for stationary turbulence is generalized for decaying turbulence and used to analyze the observed discrepancies. Based on this analysis, the time correlations are determined by the wavenumber energy spectra and the sweeping velocity, which is the square root of the total energy. Hence, an accurate prediction of the instantaneous energy spectra is most critical to the accurate computation of time correlations. (C) 2004 American Institute of Physics.

Computation of stress intensity factors by the sub-region mixed finite element method of lines

Based on the sub-region generalized variational principle, a sub-region mixed version of the newly-developed semi-analytical 'finite element method of lines' (FEMOL) is proposed in this paper for accurate and efficient computation of stress intensity factors (SIFs) of two-dimensional notches/cracks. The circular regions surrounding notch/crack tips are taken as the complementary energy region in which a number of leading terms of singular solutions for stresses are used, with the sought SIFs being among the unknown coefficients. The rest of the arbitrary domain is taken as the potential energy region in which FEMOL is applied to obtain approximate displacements. A mixed system of ordinary differential equations (ODEs) and algebraic equations is derived via the sub-region generalized variational principle. A singularity removal technique that eliminates the stress parameters from the mixed equation system eventually yields a standard FEMOL ODE system, the solution of which is no longer singular and is simply and efficiently obtained using a standard general-purpose ODE solver. A number of numerical examples, including bi-material notches/cracks in anti-plane and plane elasticity, are given to show the generally excellent performance of the proposed method.

用Symbolic Computation计算毛细重力波六阶解

本文首先运用Symbolic Computation在半物理平面(x,)上计算了毛细重力波的六阶解,得到了波形与色散关系,低阶解与 Hogan 结果一致。

Computation of turbulent-flow in a square duct - aspects of the secondary flow and its origin

A numerical study of turbulent flow in a straight duct of square cross-section is made. An order-of-magnitude analysis of the 3-D, time-averaged Navier-Stokes equations resulted in a parabolic form of the Navier-Stokes equations. The governing equations, expressed in terms of a new vector-potential formulation, are expanded as a multi-deck structure with each deck characterized by its dominant physical forces. The resulting equations are solved using a finite-element approach with a bicubic element representation on each cross-sectional plane. The numerical integration along the streamwise direction is carried out with finite-difference approximations until a fully-developed state is reached. The computed results agree well with other numerical studies and compare very favorably with the available experimental data. One important outcome of the current investigation is the interpretation analytically that the driving force of the secondary flow in a square duct comes mainly from the second-order terms of the difference in the gradients of the normal and transverse Reynolds stresses in the axial vorticity equation.

Visual Basic在油藏数值模拟中的应用

随着油藏数值模拟技术的发展以及油藏数值模拟软件的不断改进和完善,油藏数值模拟软件在油田开发中的应用越来越广泛。对油藏数值模拟软件计算出的数据进行整理不仅枯燥而且花费了大量时间。本文利用Visual Basic语言编制了处理ECLIPSE软件生成的油气田开发指标数据软件(RSMAN)。该软件方便了油气田开发指标的整理和汇总,软件界面友好,操作简单。

Computation of turbulence-generated noise by large-eddy simulation

The hybrid method of large eddy simulation (LES) and the Lighthill analogy is being developed to compute the sound radiated from turbulent flows. The results obtained from the hybrid method are often contaminated by the absence of small scales in LES, since the energy level of sound is much smaller than that of turbulent flows. Previous researches investigate the effects of subgrid sacle (SGS) eddies on the frequency spectra of sound radiated by isotropic turbulence and suggest a SGS noise model to represent the SGS contributions to the frequency spectra. Their investigations are conducted in physical space and are unavoidably influenced by boundary conditions. In this paper, we propose to perform such calculations in Fourier space so that the effects of boundary conditions can be correctly treated. Posteriori tests are carried out to investigate the SGS contribution to the sound. The results obtained recover the -7/2 law within certain wave-number ranges, but under-estimate the amplitudes of the frequency spectra. The reason for the underestimation is also discussed.

Computation of flow over a rotating body on unstructured Chimera mesh

Flow around moving boundary is ubiquitous in engineering applications. To increse the efficienly of the algorithm to handle moving boundaries is still a major challenge in Computational Fluid Dynamics (CFD). The Chimera grid method is one type of method to handle moving boundaries. A concept of domain de-composition has been proposed in this paper. In this method, sub-domains are meshed independently and governing equations are also solved separately on them. The Chimera grid method was originally used only on structured (curvilinear) meshes. However, in a problem which involves both moving boundary and complex geometry, the number of sub-domains required in a traditional (structured) Chimera method becomes fairly large. Thus the time required in the interior boundary locating, link-building and data exchanging also increases. The use of unstructured Chimera grid can reduce the time consumption significantly by the reduction of domain(block) number. Generally speaking, unstructured Chimera grid method has not been developed. In this paper, a well-known pressure correction scheme - SIMPLEC is modified and implemented on unstructured Chimera mesh. A new interpolation scheme regarding the pressure correction is proposed to prevent the possible decoupling of pressure. A moving-mesh finite volume approach is implemented in an inertial reference frame. This approach is then used to compute incompressible flow around a rotating circular and elliptic cylinder. These numerical examples demonstrate the capability of the proposed scheme in handling moving boundaries. The numerical results are in good agreement with other experimental and computational data in literature. The method proposed in this paper can be efficiently applied to more challenge cases such as free-falling objects or heavy particles in fluid.

Visual C++ 6.0学习教程

MicrosoftVisualC十十6.0作为Microsoft Visual Studio的重要组成部分，包含了迄今为止功能最为强大的基于Windows的应用框架，在同类产品中处于领先地位。VisualC十十6.0是Microsoft迄今为止最全面、最完善的程序开发工具，为了适应各种编程风格，该软件提供了各种各样的辅助工具，在发挥编程能力和提高灵活性方面达到了空前的水平。与以往VisualC十十的各种版本相比较，VisualC十十6.0在编程环境、程序语言技术等方面做了许多改进，从而使VisualC十十更加适合专业程序员快速进行应用程序的开发。

本书内容丰富、图文并茂，是一本适合各种读者学习VisualC十十6.0的优秀参考书。

目 录
第一章 VisualC十十6.0简介及安装
1.1VisualC十十6.0新特性
1.2viSualC十十6.0开发环境简介
1.3如何学习使用VisualC十十6.0
1.4VisualC十十6.0的安装
第二章 走进C十十的世界
2.1类和对象的简介
2.2继承和多态性――一个具体的例子
2.3内嵌对象
2.4在栈中申请对象
2.5全程对象的申请
2.6对象之间的相互关系――指针数据成员
2.7this指针的使用
2.8对指针的引用
2.9友元类和友元函数
2.10静态类成员
2.11重载运算符
2.12从代码中分离出类定义
2.13匈牙利表示法
第三章 VisualC十十6.0的编程环境
3.1VisualC十十6.0主窗口
3.2VisualC十十6.0工具栏
3.3VisualC十十6.0菜单栏
3.4项目与项目工作区
3.5资源与资源编辑器
第四章 编一个最简单的VC十十程序
4.1什么是AppWizard？
4.2迎接你的第一个AppWizard程序
4.3“Iamaprogrammer.”在哪儿？
第五章 程序框架入门
5.1一个简化过的程序框架
5.2WinMain（）：第一个动作
5.3登记窗口类
5.4创建一个窗口
5.5显示窗口
5.6显示出那条消息
5.7窗口类与窗口对象
第六章 消息循环
6.1在消息循环中兜圈子
6.2对事件做出响应：WindowFun（）
6.3响应不同的消息
6.4现在你还跟得上吗？
6.5设备界面进行交互
第七章 精通程序框架
7.1WinMain（）函数在哪儿？
7.2应用程序框架和源文件
7.3工具条、状态条和打印等选项
7.4程序的控制流程
第八章 使用classWizard编程
8.1使用ClassWizard添加消息处理函数
8.2classWizard功能介绍
8.3传送鼠标消息
8.4保存鼠标绘图的信息
第九章 视图与文档
9.1Document－View模式
9.2从视图中分离出文档
9.3保存文档
9.4再访MyProg2.cpp
第十章 对象连接与嵌入（OLE）及其自动化
10.1公共对象模式（COM）
10.2类厂（classfactory）
10.3OLE自动化
10.4IDispatch接口
第十一章 动态连接库（DLLs）
11.1为什么使用DLL
11.2传统的DLL
11.3MFC库DLL
11.4MyProg4A――编写自己的类库扩展DLL
11.5MyProg4B――使用MFC库扩展DLL
11.6资源访问
第十二章 图形设备接口
12.1设备环境类
12.2GDI对象
12.3Windows的颜色映射
12.4映射方式
12.5字体
12.6MyProg3例程序
12.7MyProg3B程序
12.8MyPr0g3C例程序――使用CScrollView
第十三章 对话框
13.1在状态条上显示对话控件的帮助信息
13.2利用Fi1eOpen通用对话框打开多个文件
13.3定制通用文件对话框
13.4扩展和缩减一个对话框
13.5显示一个模式或无模式对话框
13.6编写定制的DDX/DDV例程
第十四章 剖析工具Spy十＋
14.1窗体
14.2消息
14.3进程与线程
第十五章 代码调试
15.1TRACE
15.2调试框架
15.3自我诊断
15.4调试代码的作用
15.5用Dump（）显示对象的信息
15.6检查内存

Parallel Fully-Implicit Computation of MHD Acceleration Experiments

A three-dimensional MHD solver is described in the paper. The solver simulates reacting flows with nonequilibrium between translational-rotational, vibrational and electron translational modes. The conservation equations are discretized with implicit time marching and the second-order modified Steger-Warming scheme, and the resulted linear system is solved iteratively with Newton-Krylov-Schwarz method that is implemented by PETSc package. The results of convergence tests are plotted, which show good scalability and convergence around twice faster when compared with the DPLR method. Then five test runs are conducted simulating the experiments done at the NASA Ames MHD channel, and the calculated pressures, temperatures, electrical conductivity, back EMF, load factors and flow accelerations are shown to agree with the experimental data. Our computation shows that the electrical conductivity distribution is not uniform in the powered section of the MHD channel, and that it is important to include Joule heating in order to calculate the correct conductivity and the MHD acceleration.

Universal quantum computation with trapped ions in thermal motion by adiabatic passage

We propose a universal quantum computation scheme for trapped ions in thermal motion via the technique of adiabatic passage, which incorporates the advantages of both the adiabatic passage and the model of trapped ions in thermal motion. Our scheme is immune from the decoherence due to spontaneous emission from excited states as the system in our scheme evolves along a dark state. In our scheme the vibrational degrees of freedom are not required to be cooled to their ground states because they are only virtually excited. It is shown that the fidelity of the resultant gate operation is still high even when the magnitude of the effective Rabi frequency moderately deviates from the desired value.