A Cip/Multi-Moment Finite Volume Method For Shallow Water Equations With Source Terms

A novel finite volume method has been presented to solve the shallow water equations. In addition to the volume-integrated average (VIA) for each mesh cell, the surface-integrated average (SIA) is also treated as the model variable and is independently predicted. The numerical reconstruction is conducted based on both the VIA and the SIA. Different approaches are used to update VIA and SIA separately. The SIA is updated by a semi-Lagrangian scheme in terms of the Riemann invariants of the shallow water equations, while the VIA is computed by a flux-based finite volume formulation and is thus exactly conserved. Numerical oscillation can be effectively avoided through the use of a non-oscillatory interpolation function. The numerical formulations for both SIA and VIA moments maintain exactly the balance between the fluxes and the source terms. 1D and 2D numerical formulations are validated with numerical experiments. Copyright (c) 2007 John Wiley & Sons, Ltd.

Reflection of shock-wave in pur foam from rigid wall

Keller proposed that a building, a mechanical installation or a body wrapped bya layer of foam plastics may be an efficient means for protection from damage ofblast wave. However, the practical effect was beyond expectation. For example, agunner wearing the foam plastics-padded waistcoat was injured more seriously by theblast wave from a muzzle. Monti took the foam plastics as homogeneous two-phasemedium and analyzed it with the theory of dusty flow. The obtained results showthat the peak pressure behind the reflected shock wave from rigid wall with foamcoat exceeds obviously that without foam coat under the same condition. Gel'fand,Patz and Weaver made experimental observations by means of shock tubes and veri-

Approximate theory of natural-convection with small grashof number

To gain some insight into the behaviour of low-gravity flows in the material processing in space, an approximate theory has been developed for the convective motion of fluids with a small Grashof number Gr. The expansion of the variables into a series of Gr reduces the Boussinesq equation to a system of weakly coupled linearly inhomogeneous equations. Moreover, the analogy concept is proposed and utilized in the study of the plate bending problems in solid mechanics. Two examples are investigated in detail, i. e. the 2-dimensional steady flows in either circular or square infinite closed cylinder, which is horizontally imposed at a specified temperature of linear distribution on the boundaries. The results for stream function ψ, velocity u and temperature T are provided. The analysis of the influences of some parameters such as the Grashof number Gr and the Prandtl number Pr, on motions will lead to several interesting conclusions. The theory seems to be useful for seeking for an analytical solutions. At least, it will greatly simplify the complicated problems originally governed by the Navier-Stokes equation including buoyancy. It is our hope that the theory might be applicable to unsteady or 3-dimensional cases in future.

基于多矩的有限体积浸入边界法研究

浸入边界法（Immersed Boundary Method）是计算流体力学中求解具有复杂、移动边界流动问题的一类有效途径，该方法在笛卡尔坐标系上离散求解流体控制方程，并通过在控制方程中添加相应源相来代表浸入边界。尽管浸入边界法借助其简单、高效的显著特点在计算流体力学应用中显示出极强的生命力，特别是针对复杂的实际流动及动边界流动问题有着无可比拟的优势，但仍有许多问题需要进一步的研究。 本论文基于浸入边界方法及多矩VSIAM3（Volume/Surface Integrated Average Multi-Moment Method）格式提出了一种不可压缩流体求解数值格式。不可压N-S方程使用VSIAM3格式进行法进行离散，引入浸入边界法处理复杂、移动流动边界条件，使用虚拟网格方法计算动量方程修正项，同时还考虑了对连续方程的修正。VSIAM3格式是一种基于多矩的有限体积法，在方程的离散中总是使用两种或两种以上的矩，如：VIA（Volume Integrated Average）和SIA（Surface Integrated Average）。而不同的矩在求解过程中依据不同形式的控制方程使用不同的离散方法进行更新。VSIAM3格式更多的局地自由度及同时使用交错网格和同位网格的特点使浸入边界法的实施更加便利、高效。研究中，浸入边界法不仅应用于处理动力边界条件，同样可以处理热动力边界条件。 研究中对大量经典算例进行了数值实验，包括一维线性初始问题、方腔流问题、二维绕静止及振荡圆柱流动、三维绕球流动及热对流问题等。数值结果同实验值及其它计算结果保持一致，该算法可准确、高效处理具有复杂、移动边界及存在热对流的不可压流动问题，为实际应用打下了基础。

Numerical investigation on a transition prediction model

A new transition prediction model is introduced, which couples the intermittency effect into the turbulence transport equations and takes the characteristics of fluid transition into consideration to mimic the exact process of transition. Test cases include a two-dimensional incompressible plate and a two-dimensional NACA0012 airfoil. Performance of this transition model for incompressible flows is studied, with numerical results consistent to experimental data. The requirement of grid resolution for this transition model is also studied.

光滑粒子动力学方法中粒子分布与数值稳定性分析

光滑粒子动力学(SPH)作为一种拉格朗日型无网格粒子方法,已经成功地应用于包括含多相流动界面以及移动边界的可压缩和不可压缩流体运动的研究中.通过对Poiseuille流动的深入研究,探索了SPH方法中粒子分布对计算精度的影响,揭示了一种因为粒子不规则分布而导致的数值不稳定现象.研究显示,这种数值不稳定性起源于SPH方法粒子近似过程中的不连续性.使用了一种新的粒子近似格式以确保SPH方法中粒子近似的连续性.计算结果表明,这种新的粒子近似格式对于规则和不规则的粒子分布都能得到稳定精度的结果.

High Order Multi-Moment Constrained Finite Volume Method. Part I: Basic Formulation

A new high-order finite volume method based on local reconstruction is presented in this paper. The method, so-called the multi-moment constrained finite volume (MCV) method, uses the point values defined within single cell at equally spaced points as the model variables (or unknowns). The time evolution equations used to update the unknowns are derived from a set of constraint conditions imposed on multi kinds of moments, i.e. the cell-averaged value and the point-wise value of the state variable and its derivatives. The finite volume constraint on the cell-average guarantees the numerical conservativeness of the method. Most constraint conditions are imposed on the cell boundaries, where the numerical flux and its derivatives are solved as general Riemann problems. A multi-moment constrained Lagrange interpolation reconstruction for the demanded order of accuracy is constructed over single cell and converts the evolution equations of the moments to those of the unknowns. The presented method provides a general framework to construct efficient schemes of high orders. The basic formulations for hyperbolic conservation laws in 1- and 2D structured grids are detailed with the numerical results of widely used benchmark tests. (C) 2009 Elsevier Inc. All rights reserved.

一种制作掺Yb相移光纤光栅激光器的实验方案

采用遮挡法引入相移制作了掺Yb相移光纤光栅(PS-FBG)。在制作光栅的过程中,将其作为激光器的谐振腔,通过监测激光器的输出功率来确定相移大小。当激光器的输出功率开始下降时,停止曝光,此时引入的相移为π/2。为了使光栅的特性尽快稳定下来需要对光栅进行退火,这将导致引入的相移小于π/2。为了弥补退火过程中引起的相移降低,需要对退火后的光栅进行二次曝光,以使光栅的相移恢复π/2。利用该方法制作了一只光纤光栅激光器。当抽运功率为100 mW时,获得了25 mW的输出功率,信噪比(SNR)为60 dB。在1 h内

Desenvolvimento de um simulador numérico empregando o método Smoothed Particle Hydrodynamics para a resolução de escoamentos incompressíveis. Implementação computacional em paralelo (CUDA)

Neste trabalho, foi desenvolvido um simulador numérico baseado no método livre de malhas Smoothed Particle Hydrodynamics (SPH) para a resolução de escoamentos de fluidos newtonianos incompressíveis. Diferentemente da maioria das versões existentes deste método, o código numérico faz uso de uma técnica iterativa na determinação do campo de pressões. Este procedimento emprega a forma diferencial de uma equação de estado para um fluido compressível e a equação da continuidade a fim de que a correção da pressão seja determinada. Uma versão paralelizada do simulador numérico foi implementada usando a linguagem de programação C/C++ e a Compute Unified Device Architecture (CUDA) da NVIDIA Corporation. Foram simulados três problemas, o problema unidimensional do escoamento de Couette e os problemas bidimensionais do escoamento no interior de uma Cavidade (Shear Driven Cavity Problem) e da Quebra de Barragem (Dambreak).

Simulação de escoamentos incompressíveis empregando o método Smoothed Particle Hydrodynamics utilizando algoritmos iterativos na determinação do campo de pressões

Nesse trabalho, foi desenvolvido um simulador numérico (C/C++) para a resolução de escoamentos de fluidos newtonianos incompressíveis, baseado no método de partículas Lagrangiano, livre de malhas, Smoothed Particle Hydrodynamics (SPH). Tradicionalmente, duas estratégias são utilizadas na determinação do campo de pressões de forma a garantir-se a condição de incompressibilidade do fluido. A primeira delas é a formulação chamada Weak Compressible Smoothed Particle Hydrodynamics (WCSPH), onde uma equação de estado para um fluido quase-incompressível é utilizada na determinação do campo de pressões. A segunda, emprega o Método da Projeção e o campo de pressões é obtido mediante a resolução de uma equação de Poisson. No estudo aqui desenvolvido, propõe-se três métodos iterativos, baseados noMétodo da Projeção, para o cálculo do campo de pressões, Incompressible Smoothed Particle Hydrodynamics (ISPH). A fim de validar os métodos iterativos e o código computacional, foram simulados dois problemas unidimensionais: os escoamentos de Couette entre duas placas planas paralelas infinitas e de Poiseuille em um duto infinito e foram usadas condições de contorno do tipo periódicas e partículas fantasmas. Um problema bidimensional, o escoamento no interior de uma cavidade com a parede superior posta em movimento, também foi considerado. Na resolução deste problema foi utilizado o reposicionamento periódico de partículas e partículas fantasmas.

CFD-DEM simulation of propagation of sound waves in fluid particles fluidised medium

In this work, speed of sound in 2 phase mixture has been explored using CFD-DEM (Computational Fluid Dynamcis - Discrete Element Modelling). In this method volume averaged Navier Stokes, continuity and energy equations are solved for fluid. Particles are simulated as individual entities; their behaviour is captured by Newton's laws of motion and classical contact mechanics. Particle-fluid interaction is captured using drag laws given in literature.The speed of sound in a medium depends on physical properties. It has been found experimentally that speed of sound drops significantly in 2 phase mixture of fluidised particles because of its increased density relative to gas while maintaining its compressibility. Due to the high rate of heat transfer within 2 phase medium as given in Roy et al. (1990), it has been assumed that the fluidised gas-particle medium is isothermal.The similar phenomenon has been tried to be captured using CFD-DEM numerical simulation. The disturbance is introduced and fundamental frequency in the medium is noted to measure the speed of sound for e.g. organ pipe. It has been found that speed of sound is in agreement with the relationship given in Roy et al. (1990). Their assumption that the system is isothermal also appears to be valid.

Transient flow patterns and bubble slug lengths in parallel microchannels with oxygen gas bubbles produced by catalytic chemical reactions

Transient flow patterns and bubble slug lengths were investigated with oxygen gas (O-2) bubbles produced by catalytic chemical reactions using a high speed camera bonded with a microscope. The microreactor consists of an inlet liquid plenum, nine parallel rectangular microchannels followed by a micronozzle, using the MEMS fabrication technique. The etched surface was deposited by the thin platinum film, which is acted as the catalyst. Experiments were performed with the inlet mass concentration of the hydrogen peroxide from 50% to 90% and the pressure drop across the silicon chip from 2.5 to 20.0 kPa. The silicon chip is directly exposed in the environment thus the heat released via the catalytic chemical reactions is dissipated into the environment and the experiment was performed at the room temperature level. It is found that the two-phase flow with the catalytic chemical reactions display the cyclic behavior. A full cycle consists of a short fresh liquid refilling stage, a liquid decomposition stage followed by the bubble slug flow stage. At the beginning of the bubble slug flow stage, the liquid slug number reaches maximum, while at the end of the bubble slug flow stage the liquid slugs are quickly flushed out of the microchannels. Two or three large bubbles are observed in the inlet liquid plenum, affecting the two-phase distributions in microchannels. The bubble slug lengths, cycle periods as well as the mass flow rates are analyzed with different mass concentrations of hydrogen peroxide and pressure drops. The bubble slug length is helpful for the selection of the future microreactor length ensuring the complete hydrogen peroxide decomposition. Future studies on the temperature effect on the transient two-phase flow with chemical reactions are recommended.

六联苯超薄膜与弱取向外延生长

有机半导体薄膜的形态结构对器件性能具有非常重要的影响，近年来高质量有机半导体薄膜的制备成为有机半导体薄膜器件研究的核心内容之一。弱取向外延生长能够获得类单晶迁移率的高质量有机半导体薄膜。而制备出高取向、大尺寸、低缺陷、连续的高质量外延生长基底是获得类单晶外延薄膜的前提条件，这样就需要深入了解它们的生长行为、生长机理及薄膜相态特性。因此，高质量外延基底生长和弱取向外延生长行为及机理的研究就具有十分重要的理论价值和实用价值，是选择和扩展弱取向外延生长材料体系的基础。由此，本论文中工作主要分为两部分，第一部分工作是对外延基底六联苯超薄膜生长的研究：一方面生长出大尺寸、高取向、连续的超薄膜为弱取向外延生长提供高质量的外延基底，另一方面丰富和发展了有机半导体薄膜生长理论。第二部分工作是对酞菁化合物在六联苯超薄膜上弱取向外延生长行为和机理的研究，为弱取向外延材料体系的选择与扩展提供实验依据和理论指导。 首先，研究了六联苯超薄膜生长行为。研究结果表明：(1) 六联苯超薄膜在高温(＞60 oC)和低温(≤60 oC)的二氧化硅(SiO2)基底上生长具有不同的生长机理：高温时薄膜生长符合扩散受限凝聚生长(DLA)机理，低温时薄膜生长是由低有序薄膜经过亚稳分解后重组向团状岛转变。(2) 六联苯超薄膜的生长行为和相结构表明单层和双层薄膜是两种不同的相态：单层薄膜是高取向的具有液晶特性的薄膜相态，双层薄膜是长程有序的近似体相β-phase结构的结晶相。(3) 通过优化基底温度和生长参数，可以制备出高取向、大尺寸、连续的六联苯超薄膜，即可以为酞菁化合物的弱取向外延生长提供高质量的外延基底。 然后，以平面型自由酞菁(H2Pc)及酞菁锌(ZnPc)和非平面型酞菁氧钒(VOPc)为例，深入研究了酞菁化合物在六联苯超薄膜上的弱取向外延生长行为及机理。研究结果表明：(1) 弱取向外延生长的酞菁分子在六联苯超薄膜上立着生长，p-p共轭的方向平行于基底，同时酞菁分子在薄膜平面内具有规则的取向织构。这种高取向的酞菁化合物薄膜有利于载流子在薄膜平面内的传输，其迁移率达到了相应的单晶水平。(2) 由于六联苯双层及单层薄膜结构和相态的差别，平面型酞菁化合物表现出不同的外延生长行为：在六联苯双层薄膜上生长的酞菁化合物薄膜在薄膜平面内只有一种取向，对应于有公度外延生长(Commensurate Epitaxy)；在六联苯单层薄膜上生长的酞菁化合物薄膜在薄膜平面内有三种取向，同时兼具有公度外延生长和无公度外延生长(Incommensurate Epitaxy)。但非平面型VOPc由于分子排列方式及三斜晶体结构的本质，在六联苯单层及双层薄膜上都只表现出无公度外延生长。(3) 六联苯(001)晶面上突起的氢原子所形成的[110]、[1-10]和[010]沟道对酞菁分子具有强烈的预取向作用，从而形成取向的单分子柱晶核。然后在晶格匹配效应和基底沟道效应下，分别形成有公度外延生长和无公度外延生长。即晶格匹配关系和基底沟道效应为寻找弱取向外延生长有机半导体材料体系提供了理论指导依据。

Experimental Investigation of Performance of a Miniature Direct Methanol Fuel Cell in Short-Term Microgravity

Experimental study of a liquid fed direct methanol fuel cell has been conducted in different gravity environments. A small single cell with 5 cm x 5 cm active area has single serpentine channel on the graphite cathode polar plate and 11 parallel straight channels on the graphite anode flow bed. Cell voltage and current have been measured and two-phase flow in anode channels has been in situ visually observed. The experimental results indicate that the effect of gravity on power performance of the direct methanol fuel cell is large when the concentration polarization governs fuel cells operation. Gravitational effect becomes larger at higher current density. Increasing methanol feeding molarity is conducive to weaken the influence of gravity on performance of liquid fed direct methanol fuel cells. Increasing feeding flow rate of methanol solution from 6 to 15 ml/min could reduce the size of carbon dioxide bubbles, while the influence of gravity still exist. Transport phenomena inside direct methanol fuel cells in microgravity is also analyzed and discussed.

Cell cycle and apoptosis alteration of human hepatocarcinoma cells by subclinical-dose C-12(6+)-beam irradiation

Objective The purpose of this study is to investigate the effect of subdinical-dose C-12(6+)-beam irradiation on cell cycle and cell apoptosis in hepatocarcinoma cells. Materials and methods The HepG(2) cells were exposed to 0-2.0 Gy of either the C-12(6+) beam or a gamma-ray. Cell survival was detected by clonogenic assay. Cell cycle was determined by flow-cytometry analysis. The apoptosis was monitored by fluorescence microscope with DAPI staining. p53 and p21 expression were detected by Western blot. Results The G(0)/G(1) cells in the irradiated groups were significantly more than those in the control (P<0.05). The C-12(6+)-ion irradiation had a greater effect on the cell cycle of HepG(2) cells (including promoting G(1)-phase and G(2)-phase arrest) than gamma-ray irradiation. The apoptotic cells induced by C-12(6+) beam were significantly more numerous than those induced by gamma-ray (P<0.05). The carbon ions had a stronger effect on p53 and p21 expression than the gamma-ray irradiation. The survival fractions for cells irradiated by C-12(6+) beam were significantly smaller than those irradiated by gamma-ray (P<0.05).