Development of a Cell Size Relaxed Scheme for the Direct Simulation Monte Carlo Method

The conventional direct simulation Monte Carlo (DSMC) method has a strong restriction on the cell size because simulated particles are selected randomly within the cell for collisions. Cells with size larger than the molecular mean free path are generally not allowed in correct DSMC simulations. However, the cell-size induced numerical error can be controlled if the gradients of flow properties are properly involved during collisions. In this study, a large cell DSMC scheme is proposed to relax the cell size restriction. The scheme is applied to simulate several test problems and promising results are obtained even when the cell size is greater than 10 mean free paths of gas molecules. However, it is still necessary, of course, that the cell size be small with respect to the flow field structures that must be resolved.

Orientation-Dependent Adhesion Strength Of A Rigid Cylinder In Non-Slipping Contact With A Transversely Isotropic Half-Space

Recently, Chen and Gao [Chen, S., Gao, H., 2007. Bio-inspired mechanics of reversible adhesion: orientation-dependent adhesion strength for non-slipping adhesive contact with transversely isotropic elastic materials. J. Mech. Phys. solids 55, 1001-1015] studied the problem of a rigid cylinder in non-slipping adhesive contact with a transversely isotropic solid subjected to an inclined pulling force. An implicit assumption made in their study was that the contact region remains symmetric with respect to the center of the cylinder. This assumption is, however, not self-consistent because the resulting energy release rates at two contact edges, which are supposed to be identical, actually differ from each other. Here we revisit the original problem of Chen and Gao and derive the correct solution by removing this problematic assumption. The corrected solution provides a proper insight into the concept of orientation-dependent adhesion strength in anisotropic elastic solids. (c) 2008 Elsevier Ltd. All rights reserved.

Subgrid-scale contributions to Lagrangian time correlations in isotropic turbulence

The application of large-eddy simulation (LES) to particle-laden turbulence raises such a fundamental question as whether the LES with a subgrid scale (SGS) model can correctly predict Lagrangian time correlations (LTCs). Most of the currently existing SGS models are constructed based on the energy budget equations. Therefore, they are able to correctly predict energy spectra, but they may not ensure the correct prediction on the LTCs. Previous researches investigated the effect of the SGS modeling on the Eulerian time correlations. This paper is devoted to study the LTCs in LES. A direct numerical simulation (DNS) and the LES with a spectral eddy viscosity model are performed for isotropic turbulence and the LTCs are calculated using the passive vector method. Both a priori and a posteriori tests are carried out. It is observed that the subgrid-scale contributions to the LTCs cannot be simply ignored and the LES overpredicts the LTCs than the DNS. It is concluded from the straining hypothesis that an accurate prediction of enstrophy spectra is most critical to the prediction of the LTCs.

A Hybrid Scheme Based on Finite Element/Volume Methods for Two Immiscible Fluid Flows

We have successfully extended our implicit hybrid finite element/volume (FE/FV) solver to flows involving two immiscible fluids. The solver is based on the segregated pressure correction or projection method on staggered unstructured hybrid meshes. An intermediate velocity field is first obtained by solving the momentum equations with the matrix-free implicit cell-centered FV method. The pressure Poisson equation is solved by the node-based Galerkin FE method for an auxiliary variable. The auxiliary variable is used to update the velocity field and the pressure field. The pressure field is carefully updated by taking into account the velocity divergence field. This updating strategy can be rigorously proven to be able to eliminate the unphysical pressure boundary layer and is crucial for the correct temporal convergence rate. Our current staggered-mesh scheme is distinct from other conventional ones in that we store the velocity components at cell centers and the auxiliary variable at vertices. The fluid interface is captured by solving an advection equation for the volume fraction of one of the fluids. The same matrix-free FV method, as the one used for momentum equations, is used to solve the advection equation. We will focus on the interface sharpening strategy to minimize the smearing of the interface over time. We have developed and implemented a global mass conservation algorithm that enforces the conservation of the mass for each fluid.

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.

Numerical Study On Transient Flow In The Deep Naturally Fractured Reservoir With High Pressure

According to the experimental results and the characteristics of the pressure-sensitive fractured formation, a transient flow model is developed for the deep naturally-fractured reservoirs with different outer boundary conditions. The finite element equations for the model are derived. After generating the unstructured grids in the solution regions, the finite element method is used to calculate the pressure type curves for the pressure-sensitive fractured reservoir with different outer boundaries, such as the infinite boundary, circle boundary and combined linear boundaries, and the characteristics of the type curves are comparatively analyzed. The effects on the pressure curves caused by pressure sensitivity module and the effective radius combined parameter are determined, and the method for calculating the pressure-sensitive reservoir parameters is introduced. By analyzing the real field case in the high temperature and pressure reservoir, the perfect results show that the transient flow model for the pressure-sensitive fractured reservoir in this paper is correct.

激光光束实时监测与自动准直系统设计

设计了一个激光光斑实时监测与光路自动准直装置,能够实时监测激光光斑并自动准直激光输出方向。基于透镜成像原理,使用CCD探测器获得光斑的二维成像,并根据两点确定一条直线原理和使用压电陶瓷电动调整架实现光路自动准直;监测控制程序采用虚拟仪器开发软件Lab View编写,可以实时监测激光光斑模式与光斑位置抖动情况,并进行反馈控制。经测试,设计装置的调整精度达0.5μrad,反馈控制频率约1 Hz,完全可降低或消除抖动周期在1 s以上的光斑飘移。

Freestanding optical fibers fabricated in a glass chip by femtosecond laser micromachining for lab-on-a-chip application (vol 13, pg 7225, 2005)

An erratum is presented to correct the propagation loss of the freestanding optical fibers fabricated in glass chip. (c) 2006 Optical Society of America.

光学干涉图像处理中基于质量权值的离散余弦变换解包裹相位

为减少噪声对相位恢复过程的影响,快速得到正确的解包裹相位,提出了一种改进的相位解包裹方法——加权离散余弦变换解包裹算法。该方法把离散余弦变换和标识相位数据好坏的质量权值结合起来,兼有速度快和可靠度高的优势。为验证此算法,对模拟和实验得到的包裹相位图添加随机噪声和散粒噪声,同时采用加权与非加权离散余弦变换算法进行处理,所得到的解包裹结果与未加噪声的解包裹相位值进行比较,结果表明,通过加权离散余弦变换算法恢复的相位图比非加权离散余弦变换算法所恢复的相位图更接近于理想值,而且两种算法的运行速度基本相同,这证明提

利用晶体串接实现高效宽带三倍频

从耦合波方程出发,分别在小信号、高功率(1.5 GW/cm2)条件下研究KDP晶体串接三次谐波转换。当两块混频晶体的长度选择为8 mm和6 mm,晶体分别偏离原混频匹配角0.35 mrad和-0.25 mrad时可以有0.3 nm的谐波转换带宽,同时系统的三次谐波转换效率与两块混频晶体之间的距离有密切关系,当两块晶体之间的距离使从第一块混频晶体出射的光波之间的相位差改变π时,会使第一块混频晶体产生的三次谐波大部分回流到基频和倍频光,从而使转换效率大幅度下降,最合适的距离应当使光波之间的相位差改变为2π。

Large three-photon absorption cross-section in a novel class of bis-(N-carbazolyl) fluorene derivatives

Two novel symmetrical charge transfer fluorene derivatives (abbreviated as BCZF and BVCZF) with carbazole end-group as the donor moieties have been synthesized. Three-photon absorption cross-sections of these two compounds have been determined by using a Q-switched Nd:YAG laser pumped with 38 ps pulses at 1064 nm in DMF. The measured 3PA cross-sections are 140 x 10(-78) and 400 x 10(-78) cm(6) s(2) for BCZF and BVCZF, respectively. The geometries, electronic structures and electronic spectra of these two compounds are systematically studied by AM1 and ZINDO/S methods. On the basis of correct UV-vis spectra, the influence of different molecular structure on three-photon absorption cross-sections is discussed micromechanically. (C) 2005 Elsevier B.V. All rights reserved.

河南南阳独山玉的岩相结构和无损分析(英文)

河南南阳独山玉的开采历史可以追溯到新石器时代,在我国玉文化中占有重要地位。鉴于当前对独山玉进行无损鉴别的方法较少,利用质子激发X荧光技术(proton induced X-ray emission,PIXE)、X射线衍射(X-ray diffraction,XRD)、激光Raman光谱(laser Raman spectroscopy,LRS)和扫描电子显微镜(scanning electron microscope,SEM)等技术对河南南阳独山玉料进行岩石矿物学分析。结果表明:独山玉主要由钙长石矿物构成,晶粒细小且结合紧密的显微结构与独山玉具有极高的稳定性有较大关系。PIXE,XRD和Raman技术作为无损分析方法为鉴定独山玉提供了准确有效的方法,为研究贵重的古玉样品提供了技术支持。

光栅位移测量系统的误差自修正方法研究

文中介绍的误差自修正方法是通过光栅位移测量系统中单片机对光栅传感器的多个零位信号进行计数,并根据测量值和系统设定值得到的误差函数自动进行误差修正。实验结果表明,该方法对光栅位移测量系统的误差既可自动进行有效的修正,又可提高系统的测量精度。

伦纳德琼斯薄膜材料导热系数的分子动力学研究

提出从微观的角度,借助计算机工具,将薄膜破坏发展的细节展现出来的分子动力学研究的思想。使得实验上难以观察的现象变得形象而便于理解。应用分子动力学理论,使用伦纳德琼斯势函数,采用预校正积分法和虚拟外力约束标定方法,模拟薄膜体系的传热系数受体系的密度、温度的影响,同时结合体系粒子的径向分布函数和长程分布函数分析了相应的系统结构特性。另外,采用不同的模拟尺寸获得了低维材料所特有的“高温尺寸效应”。结果显示,导热系数随密度的增加变大,随温度的上升而变大。这些数据现有测量手段是难以得到的,这类模拟可以为研究提供一些