LOCAL BENDING OF THIN FILM ON VISCOUS LAYER

Effects of deposition layer position and number/density on local bending of a thin film are systematically investigated. Because the deposition layer interacts with the thin film at the interface and there is an offset between the thin film neutral surface and the interface, the deposition layer generates not only axial stress but also bending moment. The bending moment induces an instant out-of-plane deflection of the thin film, which may or may not cause the so-called local bending. The deposition layer is modeled as a local stressor, whose location and density are demonstrated to be vital to the occurrence of local bending. The thin film rests on a viscous layer, which is governed by the Navier-Stokes equation and behaves like an elastic foundation to exert transverse forces on the thin film. The unknown feature of the axial constraint force makes the governing equation highly nonlinear even for the small deflection case. The constraint force and film transverse deflection are solved iteratively through the governing equation and the displacement constraint equation of immovable edges. This research shows that in some special cases, the deposition density increase does not necessarily reduce the local bending. By comparing the thin film deflections of different deposition numbers and positions, we also present the guideline of strengthening or suppressing the local bending.

零攻角NADTPS重构高超声速流场气动力/热研究

长时间、远距离高超声速飞行要求飞行器设计兼顾减阻和防热特点,论文提出无烧蚀自适应减阻防热新概念(NADTPS)并展开相关研究。NADTPS基于支杆-钝头体结构重构流场,利用侧向射流重点防护支杆头部。本文通过数值求解二维轴对称Navier-Stokes方程,对零攻角下NADTPS重构流场进行模拟,研究流场基本特征、气动力/热特性,分析新概念的减阻和防热机理。

超声速平板绕流的全流域阻力特性研究

超声速流动中摩擦阻力的准确预测一直是气动研究中的难点。本文以零攻角的平板绕流为对象，结合传统的CFD方法求解Navier-Stokes方程和直接模拟Monte-Carlo（DSMC）方法模拟粒子的微观运动，研究平板在超声速来流下的全流域阻力特性。研究发现，在超声速流中会产生激波-边界层相互干扰作用并随着来流马赫数的增加而加强，同时平板前缘存在非平衡效应显著的区域；努森数的增大使激波和边界层均增厚，粘性作用增强，非平衡区的范围增大，非平衡强度增强，而激波强度减弱；最后通过拟合连续流条件下的摩阻系数和修正自由分子流结果，结合桥函数，获得了全流域范围内均有效的平板阻力系数表达式。

支杆-钝头体带攻角飞行流场和“军刺”挡板作用机理研究

本文通过三维Navier-Stokes方程的数值求解，研究了支杆-钝头体结构在10°攻角M∞=6.0飞行条件下的流场结构和特点，指出其气动力特性恶化的原因，提出采用“军刺”挡板改善流场和气动力特性，并通过对比两种不同挡板作用下的流场和气动力特性变化分析其作用机理，发现“军刺”挡板结构分割流场抑制三维效应形成的周向流动，迎风面形成稳定的回流区和剪切层结构，将迎风面锥激波推离轴线，降低钝头体肩部流动结构相互作用强度，并在一定程度上缓解背风面流动干扰，明显改善支杆-钝头体带攻角飞行时的气动力特性。

Heat Transfer Problems Induced By Multi-shocks Interaction

Numerical simulations of the multi-shock interactions observable around hypersonic vehicles were carried out by solving Navier-Stokes equations with the AUSMPW scheme and the new type of the IV interaction created by two incident shock waves was investigated in detail. Numerical results show that the intersection point of the second incident shock with the bow shock plays important role on the flow pattern, peak pressures and heat fluxes. In the case of two incident shocks interacting with the bow shock at the same position, the much higher peak pressure and more severe heat transfer rate are induced than the classical IV interaction. The phenomenon is referred to as the multi-shock interaction and higher requirements will be imposed on thermal protection systems.

用双流体模型模拟超声速气固两相流动

采用了双流体模型对JPL(Jet Propulsion Laboratory)喷管中气固两相流动以及超声速射流两相流动进行了数值模拟,并研究了可压缩两相流动中气相与颗粒的相互作用规律.自主开发的一般曲线坐标系下二维轴对称可压缩双流体程序Solve2D,对气相求解Navier-Stokes(N-S)方程组,采用k-ε湍流模型,颗粒相求解Euler方程组.对JPL喷管内的两相流场和湍流两相射流流场进行了数值模拟,研究了不同颗粒质量百分数以及不同颗粒直径时的气固两相流场的流动规律.

气泡/液滴运动的Level Set方法模拟研究

利用Level Set方法，结合投影法求解了描述气泡/液滴运动的Navier-Stokes方程。对地面常重力场中不同大小的空气泡在高黏度糖浆溶液中的自由上升运动现象，数值模拟结果与实验观测结果符合甚好，表明该方法能够计算大密度比和黏度比$(>1000:1)$情况下的气液两相流动。而对等密度液滴的热毛细迁移现象的数值模拟结果同样能够与实验结果相一致，表明该方法同样适于研究具有Marangoni效应的两相流动现象，特别是在空间微重力环境中的气液两相传热现象中的局部流动与传热问题。

三维波流数值水槽的构建方法

为了模拟波浪与水流的联合作用,基于Reynolds平均的Navier-Stokes方法,构建了一个三维波流数值水槽. 使用该模型对不同出流边界条件进行了对比分析. 结果表明:引入合理的入流和出流边界条件,可以保证在计算稳定后,水槽平均水面基本保持不变. 垂向时均流速分布的计算结果与实验数据吻合良好. 这表明:该波流数值水槽具有较好的波流特性,可用于计算波浪与水流、以及与三维固定结构之间相互作用.

可压缩湍流直接数值模拟

本书是以作者及其研究小组多年来的研究工作为基础撰写的，主要是可压缩湍流的直接数值模拟。全书共分两部分。第一部分为数值模拟方法：针对可压缩湍流的多尺度非定常非稳定流动特征，所构造的高精度差分格式（高于三阶精度），分析了各种高精度差分格式的数值行为特征和模拟小尺度物理结构、分辨非定常间断面的能力，重点分析了正确模拟多尺度物理结构及其非线性干扰对数值方法的要求。第二部分是物理问题的直接数值模拟：根据可压缩湍流的流动物理特征，所建立的求解可压缩Navier-Stokes方程的高精度数值模拟方法，详细描述了不稳定流和可压缩湍流的直接数值模拟，以及在此基础上所发现的一些新的物理现象。直接数值模拟主要针对典型可压缩湍流，包括均匀各向同性湍流、槽道湍流、平面混合流、射流、平板边界层湍流和钝楔、钝锥边界层湍流等，重点是可压缩湍流从层流转捩到湍流的直接数值模拟及其流动机理的研究。书中还介绍了我课组所研制的高精度、高效并行算法的软件。

Enhancements of the simulation method on the edge effect in resin transfer molding processes

In resin transfer molding processes, small clearances exist between the fiber preform and the mold edges, which result in a preferential resin flow in the edge channel and then disrupt the flow patterns during the mold filling stage. A mathematical model including the effect of cavity thickness on resin flow was developed for flow behavior involving the interface between an edge channel and a porous medium. According to the mathematical analysis of momentum equations in a fully developed rectangular duct and formulations of the equivalent edge permeability, comparing with three-dimensional Navier-Stokes equations, the governing equations were modified in the edge channel. The volume of fluid (VOF) method was applied to track the flow front. A simple case is numerically simulated using the modified governing equations. The effects of edge channel width and cavity thickness on flow front and inlet pressure are analyzed, and the evolution characteristics of simulated results are in agreement with the experimental results. (c) 2007 Elsevier B.V. All rights reserved

Reactive Mold Filling in Resin Transfer Molding Processes with Edge Effects

Reactive mold filling is one of the important stages in resin transfer molding processes, in which resin curing and edge effects are important characteristics. On the basis of previous work, volume-averaging momentum equations involving viscous and inertia terms were adopted to describe the resin flow in fiber preform, and modified governing equations derived from the Navier-Stokes equations are introduced to describe the resin flow in the edge channel. A dual-Arrhenius viscosity model is newly introduced to describe the chemorheological behavior of a modified bismaleimide resin. The influence of the curing reaction and processing parameters on the resin flow patterns was investigated.

Rational Design, Synthesis, and Optical Properties of Film-Forming, Near-Infrared Absorbing, and Fluorescent Chiromophores with Multidonors and Large Heterocyclic Acceptors

A new series of film-forming, low-bandgap chromophores (1a,b and 2a,b) were rationally designed with aid of a computational study., and then synthesized and characterized. To realize absorption and emission above the 1000 nm wavelength, the molecular design focuses on lowering the LUMO level by fusing common heterocyclic units into a large conjugated core that acts an electron acceptor and increasing the charge transfer by attaching the multiple electron-donating groups at the appropriate positions of the acceptor core. The chromophores have bandgap levels of 1.27-0.71 eV, and accordingly absorb at 746-1003 nm and emit at 1035-1290 nm in solution. By design, the relatively high molecular weight (up to 2400 g mol(-1)) and non-coplanar structure allow these near-infrared (NIR) chromophores to be readily spin-coated as uniform thin films and doped with other organic semiconductors for potential device applications. Doping with [6,6]-phenyl-C-61 butyric acid methyl ester leads to a red shift in the absorption on]), for la and 2a. An interesting NIR electrochromism was found for 2a, with absorption being turned on at 1034 nm when electrochemically switched (at 1000 mV) from its neutral state to a radical cation state. Furthermore, a large Stokes shift (256-318 nm) is also unique for this multidonor-acceptor type of chromophore.

N-Methyl-N-Allylpyrrolidinium Based Ionic Liquids for Solvent-Free Dye-Sensitized Solar Cells

We prepared four new ionic liquids consisting of N-methyl-N-allylpyrrolidinium cation in conjunction with anions including iodide, nitrate, thiocyanate, and dicyanamide, respectively, and measured their physical properties of density, viscosity, and conductivity. Owing to the relatively lower melting point of electroactive N-methyl-N-allylpyrrolidinium iodide, in combination with three other nonelectroactive ionic liquids, we could construct solvent-free electrolytes possessing high iodide concentrations for dye-sensitized solar cells. We correlated temperature-dependent electrolyte viscosity with molar conductivity and triiodide mobility through applying an empirical Walden's rule and a modified Stokes-Einstein equation, respectively. We have further found that these anions (nitrate, thiocyanate, and dicyanamide) have different influences on surface states and electron transport in the mesoporous titania film, resulting in different photovoltages and photocurrents of dye-sensitized solar cells.

Syntheses of Ag, PbSe, and PbTe nanocrystals and their binary self-assembly exploration at low size-ratio

Nanocrystals of Ag, PbSe, and PbTe were prepared via a high-temperature organic solution approach, respectively. Using a size-selection technique, the size-distribution of each set of nanocrystals could be fine-tuned and finally monodisperse products were achieved. Superlattice structure of binary self-assemblies in low size-ratio were also explored and characterized by transmission electron microscopy. It is realized that a success of achieving binary self-assembly pattern is greatly dependent on several key factors including particle size-distributions, relative concentrations of both components, as well as the size-ratios between Ag and PbSe (or PbTe) nanocrystals.