High‐order accurate numerical solutions of incompressible flows with the artificial compressibility method

A high‐order accurate, finite‐difference method for the numerical solution of incompressible flows is presented. This method is based on the artificial compressibility formulation of the incompressible Navier–Stokes equations. Fourth‐ or sixth‐order accurate discretizations of the metric terms and the convective fluxes are obtained using compact, centred schemes. The viscous terms are also discretized using fourth‐order accurate, centred finite differences. Implicit time marching is performed for both steady‐state and time‐accurate numerical solutions. High‐order, spectral‐type, low‐pass, compact filters are used to regularize the numerical solution and remove spurious modes arising from unresolved scales, non‐linearities, and inaccuracies in the application of boundary conditions. The accuracy and efficiency of the proposed method is demonstrated for test problems. Copyright © 2004 John Wiley & Sons, Ltd.     

Modelling of a single drop impact onto liquid film using particle method

The process of single liquid drop impact on thin liquid surface is numerically simulated with moving particle semi‐implicit method. The mathematical model involves gravity, viscosity and surface tension. The model is validated by the simulation of the experimental cases. It is found that the dynamic processes after impact are sensitive to the liquid pool depth and the initial drop velocity. In the cases that the initial drop velocity is low, the drop will be merged with the liquid pool and no big splash is seen. If the initial drop velocity is high enough, the dynamic process depends on the liquid depth. If the liquid film is very thin, a bowl‐shaped thin crown is formed immediately after the impact. The total crown subsequently expands outward and breaks into many tiny droplets. When the thickness of the liquid film increases, the direction of the liquid crown becomes normal to the surface and the crown propagates outward. It is also found that the radius of the crown is described by a square function of time: r_C = [c(t ? t₀)]^0.5. When the liquid film is thick enough, a crown and a deep cavity inside it are formed shortly after the impact. The bottom of the cavity is initially oblate and then the base grows downward to form a sharp corner and subsequently the corner moves downward. Copyright © 2004 John Wiley & Sons, Ltd.     

Highly strained InxGa(1−x)As−InyAl(1−y)As (x>0.8,y<0.3) layers for short wavelength QWIP and QCL structures grown by MBE

Highly strained quantum cascade laser (QCL) and quantum well infrared photodetector (QWIPs) structures based on In_xGa_(1−x)As−In_yAl_(1−y)As (x>0.8,y<0.3) layers have been grown by molecular beam epitaxy. Conditions of exact stoichiometric growth were used at a temperature of 420°C to produce structures that are suitable for both emission and detection in the 2–5 μm mid-infrared regime. High structural integrity, as assessed by double crystal X-ray diffraction, room temperature photoluminescence and electrical characteristics were observed. Strong room temperature intersubband absorption in highly tensile strained and strain-compensated In_0.84Ga_0.16As/AlAs/In_0.52Al_0.48As double barrier quantum wells grown on InP substrates is demonstrated. Γ–Γ intersubband transitions have been observed across a wide range of the mid-infrared spectrum (2–7 μm) in three structures of differing In_0.84Ga_0.16As well width (30, 45, and 80 Å). We demonstrate short-wavelength IR, intersubband operation in both detection and emission for application in QC and QWIP structures. By pushing the InGaAs–InAlAs system to its ultimate limit, we have obtained the highest band offsets that are theoretically possible in this system both for the Γ–Γ bands and the Γ–X bands, thereby opening up the way for both high power and high efficiency coupled with short-wavelength operation at room temperature. The versatility of this material system and technique in covering a wide range of the infrared spectrum is thus demonstrated.     

Self-bonding in an amorphous polymer below the glass transition: A T-peel test investigation

Films of amorphous polystyrene (PS) with a weight-average molecular weight (M_w) of 225 × 10³ g/mol were bonded in a T-peel test geometry, and the fracture energy (G) of a PS/PS interface was measured at the ambient temperature as a function of the healing time (t_h) and healing temperature (T_h). G was found to develop with (t_h)^1/2 at T_h = T_g-bulk − 33 °C (where T_g-bulk is the glass-transition temperature of the bulk sample), and log G was found to develop with 1/T_h at T_g-bulk − 43 °C ≤ T_h ≤ T_g-bulk − 23 °C. The smallest measured value of G = 1.4 J/m² was at least one order of magnitude larger than the work of adhesion required to reversibly separate the PS surfaces. These three observations indicated that the development of G at the PS/PS interface in the temperature range investigated (<T_g-bulk) was controlled by the diffusion of chain segments feasible above the glass-transition temperature of the interfacial layer, in agreement with our previous findings for fracture stress development at several polymer/polymer interfaces well below T_g-bulk. Close values of G = 8–9 J/m² were measured for the symmetric interfaces of polydisperse PS [M_w = 225 × 10³, weight-average molecular weight/number-average molecular weight (M_w/M_n) = 3] and monodisperse PS (M_w = 200 × 10³, M_w/M_n = 1.04) after healing at T_h = T_g-bulk − 33 °C for 24 h. This implies that the self-bonding of high-molecular-weight PS at such relatively low temperatures is not governed by polydispersity. © 2004 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 42: 1861–1867, 2004     

亚像元动态成像技术中系统调制传递函数的数值分析

亚像元动态成像技术是目前实现遥感器小型化非常有效的新方法,因而研究亚像元动态成像系统的调制传递函数MTF显得非常必要和重要.介绍了亚像元动态成像的基本原理及成像过程,指出了光程差对亚像元动态成像系统MTF的影响.对亚像元动态成像系统的MTF进行了数值分析,并指出了亚像元动态成像系统的成像质量及空间分辨率与一般成像系统相比较的优点.     

The long‐time behaviour of the thermoconvective flow in a porous medium

For the Boussinesq approximation of the equations of coupled heat and fluid flow in a porous medium we show that the corresponding system of partial differential equations possesses a global attractor. We give lower and upper bounds of the Hausdorff dimension of the attractor depending on a physical parameter of the system, namely the Rayleigh number of the flow. Numerical experiments confirm the theoretical findings and raise new questions on the structure of the solutions of the system. Copyright © 2004 John Wiley & Sons, Ltd.     

晶体微观结构的数值计算

晶体微观结构是晶体材料在特定物理条件下其多个能量极小平衔态在空间形成的某种微尺度的规则分布．几何非线性的连续介质力学理论可以用能量极小化原理来解释晶体微观结构的形成，并用Young测度来刻画平衡态各变体在空间的概率分布．定性的理解与定量地分析和计算晶体材料的微观结构对于发展和改进高级晶体功能材料，如形状记忆合金、铁电体、磁至伸缩材料等，有重要的意义．本文回顾了近年来晶体微观结构数值计算方面的最新进展．介绍了计算晶体微观结构的几种数值方法及有关的数值分析结果。     

光辐射计量测试技术

简要介绍国外光辐射计量采用的主要方法、测试手段和发展趋势。为了满足国防系统对光辐射计量的需要建立了一批计量标难，解决了光辐射源辐射特性和探测器特性计量测试问题。     

Continuation of travelling‐wave solutions of the Navier–Stokes equations

An efficient way of obtaining travelling waves in a periodic fluid system is described and tested. We search for steady states in a reference frame travelling at the wave phase velocity using a first‐order pseudospectral semi‐implicit time scheme adapted to carry out the Newton's iterations. The method is compared to a standard Newton–Raphson solver and is shown to be highly efficient in performing this task, even when high‐resolution grids are used. This method is well suited to three‐dimensional calculations in cylindrical or spherical geometries. Copyright © 2006 John Wiley & Sons, Ltd.     

Electronic molecule-like spectrum of quantum wells in crossed fields

We show that an electron confined to a single finite parabolic quantum well in crossed electric and magnetic fields can behave as a double quantum well system. The magnetic field is parallel to the heterostructure layers and the electric field is perpendicular to those. For a suitable choice of both fields and quantum well width, the electron can be confined to a double quantum well effective potential that is very similar to the electronic potential model for diatomic molecules. The double quantum well spectrum is calculated using a numerical algorithm based on semiclassical methods. A physical interpretation of this quantum system is given based on the analogy to the electrons bound to diatomic molecules.