薄夹层界面热弛豫解及其在冲击温度研究中的意义

 对样品/薄夹层/窗口组成的三层介质双界面热传导问题进行求解。在3个不同初温、不同热物性参数及边界条件下，获得了夹层介质中的温度解，给出了界面温度历史随薄层厚度、热扩散系数和介质之间热扩散系数之比的变化规律。从理论分析和实验结果发现，当夹层厚度小于1 μm，热扩散系数大于10^-4 m²/s时，在实验观测的时间内，界面温度历史将很快衰减到一平台，其幅度只由样品、窗口的初温和热扩散系数确定，而与夹层介质的初温、热传导系数无关。     

气体在微孔介质中的等效扩散系数预测

本文通过数值模拟研究了微孔介质内部结构对气体等效扩散系数的影响。采用随机生成生长方法生成多孔介质的三维微观结构,使用高效的格子-Boltzmann方法求解扩散的控制方程,可预测微孔结构内气体的等效扩散系数。模拟结果表明:微孔材料的内部结构对多孔介质中气体扩散及等效扩散系数有显著影响。对于高孔隙率情况,基于Fick定律数值预测的气体等效扩散系数和实验结果吻合较好;对于部分饱和的颗粒状微孔介质,孔隙内液滴阻碍了气体的扩散路径,多孔介质的无量纲等效扩散系数与饱和度呈非线性关系。     

分区均匀介质静电场的求解

1分区均匀介质静电场求解方法的回顾分区均匀介质静电场的求解是电磁学课程中比较有难度和有特色的内容.有两种情形的分区均匀介质静电场问题可以解析求解:一种是介质-介质界面与撤去电介质时的电场线平行;另一种是介质-介质界面与撤去电介质时的电场线垂直.它是唯一性定理应用的一个生动的实例,并且涉及到静     

运动网格上的水下爆炸数值模拟方法

 在对水下爆炸过程的数值模拟中,把水-气介质界面当作一种特殊的内部边界,并且在该边界上的网格节点都有双重定义,分别对应边界两侧的不同介质状态,通过利用HLLC（Harten,Lax,van Leer,Contact）方法求解任意ALE（Arbitrary Largrangian Eulerian Method,拉格朗日-欧拉方程）方程组,以及运用界面跟踪和动网格方法,对一维球对称水下爆炸模型进行了数值模拟。结果表明,所得到的数值计算结果与实验结果吻合得较好,并且在扩展到二维柱对称水下爆炸模型时,同样得到了满意的结果。     

多介质流体界面不稳定性程序

研制了二维多介质流体程序，主要包括单介质内高精度流体力学计算，多介质混合网格内各种介质输运过程和压力驰豫平衡过程计算、实际状态方程的黎曼解计算。流体计算分别采用高分辨两步PPM（Parabolic Piecewise Method）算法、TVD（Total Variation Diminishing）算法和FCT（Flux Corrected—Transport）算法，流体界面追踪采用VOF（Volume-of-Fluid）。数值求解可压缩多流体方程组和可压缩VOF方程。二维界面追踪分别采用一阶精度Youngs方法和二阶精度Elivira方法，三维界面追踪采用一阶精度Youngs方法，     

激波和单列水柱、气水固多界面相互作用的数值模拟

研究可压缩多介质流场的激波和多介质界面相互作用问题.在Descartes固定网格采用level-set方法追踪界面,气/气界面边界条件处理采用OGFM方法,采用修正的rGFM方法提高气/水和气/固界面处构造Riemann问题精度,将Riemann近似解得到的界面参数外推到两侧真实和虚拟流体,采用五阶WENO方法求解流场Euler方程和界面level-set方程,给出不同时刻流场数值纹影图像.结果表明：在可压缩流场嵌入固体和水、气体等目标,本文方法可较精确地分辨平面运动激波和单列水柱及包含气/气、气/水和气/固等界面作用后产生的复杂激波结构.和传统的分区与贴体变换方法不同,为Descartes网格包含多介质界面复杂流场计算提供新途径.     

在直角坐标系用NGFM模拟复杂计算域水下高压气泡膨胀问题

采用NGFM(New version of Ghost Fluid Method)处理复杂计算域的固壁边界,用RGFM(Real Ghost Fluid Method)求解气-水界面附近网格节点的状态参数,从而在直角坐标系下对复杂计算域的水下高压气泡膨胀问题进行数值模拟。流场控制方程选用Euler方程,用五阶WENO格式离散空间导数项,二阶Runge-Kutta法离散时间导数项;气-水界面追踪使用Level Set方法,对Level Set方程,用五阶HJ-WENO(Hamilton-Jacobi WENO)和三阶Runge-Kutta法求解。将计算结果与任意坐标系下的结果进行对比,验证了NGFM在笛卡尔网格中处理复杂形状固壁边界的可行性。得到了水下流场压力等值线图、高压气泡的演变过程以及特定点处的压力-时间曲线。计算结果表明,高压气泡在固壁反射激波的作用下,膨胀过程受到抑制;强激波在固壁的反射会导致固壁附近出现大范围的空化流动。     

同轴流动聚焦中射流不稳定性的理论研究

壳核结构的微胶囊在医学药学材料食品农业等领域具有广泛的应用前景, 其制备方法一直是相关领域关注的焦点.同轴流动聚焦(co-flow focusing)是一种新型制备技术, 利用复合射流的破碎制备微胶囊具有包裹率高过程量化可控参数域广产率高等诸多优势.在实验中, 复合射流的破碎受到多个过程参数的影响, 并涉及了多层界面的耦合效应.利用简化的物理模型, 在时间和时空域中分析了三相水-油-水复合射流不稳定性的发展和演化.在黏性流体线性稳定性理论中, 同轴射流和驱动液体的基本速度型分别基于管流和误差函数构造, 并通过数值方法求解满足相应边界条件下的线化小扰动控制方程.结果表明:增加内外层界面的界面张力均有利于射流的破碎; 流体的黏性对同轴射流的稳定性均有着促进作用; 越大的黏性越小的内界面张力对应着越大的射流破碎波长; 内外界面的耦合作用以及复合液滴的包裹情况均与内外射流的半径比息息相关; 绝对-对流不稳定性转换的临界Weber数随Reynolds数内层界面张力的增大而增大, 随内层和驱动流体的黏性增大而减小.这些结果将有助于提高液体驱动下同轴流动聚焦技术的过程控制, 为实际应用提供理论指导.      

一维多介质可压缩流动数值方法

应用高精度界面追踪方法计算一般状态方程的多介质可压缩流动问题;应用LevelSet技术捕捉界面位置,在界面附近采用守恒数值离散,用双波近似求解一般状态方程Riemann问题,并采用统一高阶PPM格式进行内点和交界面点的计算.一维算例表明,该方法对于光滑区域以及多介质交界面具有二阶精度,能准确地模拟交界面的位置,交界面计算无数值振荡和数值耗散,并能处理一般状态方程的多介质可压缩流动问题.     

30 keV He~(2+)在不同倾斜角度的聚碳酸酯微孔膜中的传输过程

本文测量了30 keV的He~(2+)入射倾斜角度分别为-0.5~?,-1~?,-1.5~?和-2.5~?的聚碳酸酯纳米微孔膜后,出射粒子角度分布、电荷态分布以及相对穿透率随时间的演化.当微孔膜倾斜角度在-0.5~?,-1~?和-1.5~?时,出射的He~(2+)离子始终保持在入射束流方向,出射的He~0原子出射方向由微孔孔道方向逐渐转移到入射束流方向,在实验过程中观测到明显的电荷交换,这一现象与之前发现的导向效应不同,微孔内部沉积的电荷斑和微孔内表面原子的短程集体散射作用,克服入射离子的横向动量,使入射离子在微孔内表面以上以类似镜面掠射的方式出射,并发生时间演化效应,主要传输机制为电荷斑辅助的表面以上的类似镜面掠射行为.而当倾斜角度在-2.5~?时,出射的He~(2+)离子始终保持在入射束流方向,出射的He~0原子始终保持在微孔孔道方向,沉积的电荷斑很难克服入射离子的横向动量,没有时间演化效应,主要传输机制为微孔内表面以下的多次随机非弹性碰撞过程.这一物理图像使中能离子入射不同倾斜角度的微孔膜物理认识更加深入和完整.     

An NMR protocol for determining ice crystal size distributions during freezing and pore size distributions during freeze-drying

Nuclear magnetic resonance water proton relaxometry is widely used to investigate pore size distributions and pore connectivity in brine-saturated porous rocks and construction materials. In this paper we show that, by replacing water with acetone, a similar method can be used to probe the porous structure of freeze-dried starch gels and therefore the ice crystal size distribution in frozen starch gels. The method relies on the observation that the starch surface acts as a powerful relaxation sink for acetone proton transverse magnetization so that Brownstein-Tarr theory can be used to extract the pore size distribution from the relaxation data. In addition the relaxation time distribution is found to depend on the spectrometer frequency and the Carr-Purcell-Meiboom-Gill pulse spacing, consistent with the existence of large susceptibility-induced field gradients within the pores. The potential of this approach for noninvasively measuring ice crystal size distributions during freezing and pore size distributions during freeze-drying in other food systems is discussed.     

An experimental investigation into the icing and melting process of a water droplet impinging onto a superhydrophobic surface

The freezing and melting process of a small water droplet on a superhydrophobic cold surface was investigated using the Laser Induced Fluorescence (LIF) technique. The superhydrophobic surface was prepared using a sol-gel method on a red copper test plate. From the obtained fluorescence images, the phase transition characteristics during the freezing and melting process of a water droplet were clearly observed. It was found that, at the beginning of the droplet freezing process, liquid water turned into ice at a very fast rate. Such phase transition process decreased gradually with time and the volume of frozen ice approached a constant value at the end of the icing process. In addition, the freezing time was found to reduce with the decrease of the test plate temperature. Besides, when the test plate temperature is relatively high, the effect of droplet volume on the freezing time is very significant. Over all, we provide some tentative insights into the microphysical process related to the icing and melting process of water droplets.     

气孔分布特性对含水多孔介质气体扩散的影响

本文基于多孔介质的气孔分布特性,计算了多孔介质在含水状态下的扩散性能,并且比较了采用两种方式计算相对渗透率时的相对扩散性能。其结果表明,基于气孔分布的计算结果低于与气孔分布无关的计算结果。另外,疏水性含水多孔介质的扩散性能低于亲水性含水多孔介质的扩散性能,基于气孔分布计算含水多孔介质的气体扩散性能时,Wyllie公式并不适用。     

NMR characterization of the pore structure and anisotropic self-diffusion in salt water Ice

NMR imaging and one- and two-dimensional self-diffusion propagator measurements of the liquid phase in salt water ice are presented. The properties of the network of brine-filled pores are found to depend on the growth conditions of the ice. Two types of samples are compared: (a) shock-frozen ice produced in the probe in situ and (b) ice grown over several hours under controlled conditions. By shock-freezing, an ice structure could be produced which featured streak-like porous channels of diameters of up to 300 μm allowing almost unrestricted self-diffusion along one preferential axis but reduced diffusivities in the remaining directions. In ice grown under controlled conditions, the pore sizes are near the resolution limit of the imaging experiment of typically 50 μm. For this type of samples, strongly non-Gaussian self-diffusion propagators are obtained, indicating restricted self-diffusion on rms scales of 30 μm. Common to all samples was the observation of highly anisotropic self-diffusion. One- and two-dimensional propagators are compared in order to estimate the degree of anisotropy and the size of the restrictions.     

氢化多孔石墨烯反渗透特性及机理分析

单层石墨烯凭借超薄的厚度和优异的力学化学防污性能,成为新一代纳滤膜材料的最佳选择之一.本文采用经典分子动力学方法,研究了氢化多孔石墨烯反渗透膜对盐水的反渗透特性.结果表明,水渗透量会随着驱动力、孔径和温度的增加而增加;而孔径大于水合半径的条件下,盐离子截留率会随驱动力和温度的增加而降低.当反渗透膜和盐水存在切向运动时,随着切向速度的增加可以有效提高盐离子截留率和减弱浓差极化现象,但也在一定程度上牺牲了水通量.通过分析水流沿渗透方向的能障分布、水分子的氢键分布和离子水合状态,解释了各参数变化对盐水在氢化多孔石墨烯中反渗透特性的影响机理.研究结果将提供基于单层多孔石墨烯反渗透特性的理论认识,并将为纳米级反渗透膜的设计提供帮助.     

Ordered porous TiO2 films obtained by freezing and the application in dye sensitized solar cells

TiO₂ films with ordered porous structure were prepared by freezing. By simply freezing wet TiO₂ film on a cold copper plate, ice crystals could grow from bottom of film. Removing such ice crystals then led to ordered, micro-sized pores in the films. With assistance of scanning electron microscopy and N₂ adsorption–desorption isotherms, micro-structural properties of the films were studied. Coarsening behavior of ice crystals was analyzed based on evolvement of the micro-structure. It was found that, both larger film thickness and longer freezing time were beneficial for the formation of ordered porous structure, which was caused by enhanced ice crystals growth. The films were then used to fabricate photoanode of dye sensitized solar cells, of which the photo-to-electric power conversion efficiency was evaluated. It was observed that, solar cell made of ordered porous TiO₂ film came out with 13% larger photo-current density comparing to that made of conventional doctor blading method, which is due to the enhanced light scattering by the ordered porous structure.     

Influence of Position on the Microstructure of Carbon Black/Polyvinyl Alcohol Composite Obtained by the Directional Freeze-drying Process

Carbon black (CB)/polyvinyl alcohol (PVA) composites with porous structure were obtained by a directional freeze-drying process of their solution in water. The microstructures of this composite at different positions through the thickness were investigated. The study shows that the composite formed a lamellar structure at the lower surface of the specimen because of the basic crystallographic and crystal growth characteristics of ice. The composite near the top surface of the specimen showed an alignment direction of the ice vertical to the container's wall and the direction of solution lowering into the liquid nitrogen. The alignment direction deflects gradually and finally converges with the aligned direction grown from the bottom position near the middle part of the CB/PVA composite. The ice crystals grow very fast along the direction of temperature gradients, which results in that a small fraction of PVA solute was entrapped within the ice crystals, leading to the formation of trans-bridge lamellar structure near the top surface of the specimen. This result shows that the position in the sample has a great influence on the structure of the porous CB/PVA composite obtained by the directional freeze-drying process.     

Mathematical Study of the Water-Vapor Permeability of the Surface Layer of a Homogeneous Porous Material

Permeability in relation to gases and liquids is one of the most important characteristics of porous materials. A porous material can interact with gases and liquids flowing through it in different ways that depend on a material’s permeability. Studies of the interaction of water vapor with materials with uniformly distributed pores are of considerable practical interest, since many of these materials are used in building and construction. Interest is also due to the possibility of expanding the study results obtained for an individual pore to a porous medium if this medium can be represented as a structure with uniformly distributed pores with sufficient accuracy. The dependence of the permeability of an individual cylindrical pore on its radius, length, and characteristics of the process of water-molecule interaction with the pore walls is studied.     

Wave Simulation in Frozen Porous Media

We propose a numerical algorithm for simulation of wave propagation in frozen porous media, where the pore space is filled with ice and water. The model, based on a Biot-type three-phase theory, predicts three compressional waves and two shear waves and models the attenuation level observed in rocks. Attenuation is modeled with exponential relaxation functions which allow a differential formulation based on memory variables. The wavefield is obtained using a grid method based on the Fourier differential operator and a Runge–Kutta time-integration algorithm. Since the presence of slow quasistatic modes makes the differential equations stiff, a time-splitting integration algorithm is used to solve the stiff part analytically. The modeling is second-order accurate in the time discretization and has spectral accuracy in the calculation of the spatial derivatives.     

聚乙烯醇水溶液二维定向凝固的微观组织演化

为了深入探究定向多孔聚合物材料的微观组织形成机理,利用定向凝固原位实时观察手段,研究不同浓度及不同分子量聚乙烯醇(PVA)水溶液在不同抽拉速度下的定向凝固形貌演化.PVA水溶液的定向凝固形态在低浓度(1 wt%,2.5 wt%)和小分子量(M_w=24000)情况下,一次枝晶间距随着抽拉速度的增加而减小.随着PVA浓度和分子量的增加,一次枝晶间距随抽拉速度变化不明显,枝晶主轴尺寸则随速度增加呈现减小的趋势.与传统凝固形态形成机理相比,PVA水溶液的凝固形态由PVA分子的扩散引起的凝固界面不稳定性机理和PVA高分子链交联引起的局部相分离机理竞争决定.