T形微通道中互不相溶两相流数值模拟

采用摄动有限体积(PFV)算法和水平集(level set)技术对T形微通道内互不相溶两相流动进行了数值模拟研究。考察了两相界面张力和微通道壁面润湿性对流动的影响,精确地捕捉到了油水两相流动的界面。对一些典型的T形微通道油水两相流动进行了数值计算,模拟结果和实验结果吻合较好。分析总结出了微通道内两相流动过程中的一些基本规律,为微通道内的液液两相流动实验设计和工业应用提供了新的数值预测手段.

蒸汽-冷流体接触冷凝流动的数值模拟

介绍了关于蒸汽-冷流体直接接触冷凝流动与传热的数值计算模型与部分研究结果。用Level Set方法确定蒸汽-冷流体接触界面的位置和形状,建立了对蒸汽和冷流体普遍适用的动量、能量和质量守恒方程,在能量和质量寺恒方程中增加了部分项用于计算蒸汽冷凝所产生的影响。用有限差分法在交错网格上离散控制方程,用Runge-Kutta法-五阶WENO组合格式求解Level Set输运方程,用压力修正的迭代Projection方法求解动量方程,而用SIMPLE方法求解温度控制方程。对算例的计算结果表明,本文所建立的数值计算模型能反映物理现象的宏观特性。根据计算结果,分析了本文模型的优缺点,并指出了今后改进的方向。

气-液两相流界面迁移现象的数值模拟研究

借助于Level Set函数,建立了气-液两相的统一控制方程组,并在交错网格中进行离散.用两种格式,即Superbee-TVD格式和5阶WENO格式求解Level Set函数的输运方程,用SIMPLER算法的思想对主流场控制方程的求解方法进行改进.数值实验结果表明,在求解Level Set的控制方程时,5阶WENO方法比Superbee-TVD格式的结果更准确;用改进的数值算法可成功实现对密度比大于1 000/1的气-液两相流界面迁移问题的数值模拟.对几种典型大密度比气-液两相流问题的计算结果与实际问题的物理规律完全一致,验证了该方法的有效性和可靠性.

Influência do sistema de cimentação e do modo de ativação na resistência ao push-out de pinos de fibra de vidro

O objetivo deste estudo foi avaliar a influência do tipo de sistema de cimentação (condicionamento ácido total ou autoadesivo), do modo de ativação (autoativado ou dual), do terço do conduto radicular (cervical, médio ou apical) e da espessura do filme de cimento sobre a resistência de união de pinos de fibra de vidro cimentados em dentes humanos. Quarenta raízes foram incluídas em resina epóxi, submetidas a tratamento endodôntico e obturadas com guta percha e cimento endodôntico sem eugenol. Decorridos sete dias, os condutos foram preparados a uma profundidade de 10mm com brocas padronizadas do sistema dos pinos de fibra (WhitePost DC #2) e aleatoriamente divididos em 4 grupos, conforme o sistema de cimentação e o modo de ativação: (G1) RelyX ARC/Adper Scotchbond Multi-Purpose Plus (condicionamento ácido total), ativação dual, (G2) RelyX ARC/Adper Scotchbond Multi-Purpose Plus, autoativado, (G3) RelyX U100 (autoadesivo), dual e (G4) RelyX U100, autoativado. Após uma semana, cada raiz foi seccionada em máquina de corte, originando 6 fatias de 1 mm de espessura (n=60). Antes do ensaio de push-out cada fatia foi fotografada em ambas as faces, para determinação do raio dos pinos e da espessura do filme de cimento. Após o ensaio mecânico, novas imagens foram capturadas para determinação do modo de falha. Para automatizar a determinação da espessura de cimento, foi desenvolvida uma macro no software KS 400. Os dados foram estatisticamente analisados com ANOVA 3 fatores (resistência de união) e teste de Kruskall-Wallis (espessura do cimento). Comparações múltiplas foram realizadas com o teste Student-Newman-Keuls. Análise de regressão, modelo linear, foi empregada para verificar a correlação entre espessura do cimento e resistência de união. Todos os testes foram aplicados com α = 0,05. O fator cimento exerceu influência significativa para a resistência de união (p = 0,0402): o RelyX U100 apresentou a maior média. A ativação dual elevou os valores de resistência de união em comparação ao modo quimicamente ativado (p < 0,0001). Houve diferenças significantes entre os grupos, sendo G1 (22,4 4,0 MPa) > G3 (20,4 3,6 MPa) > G4 (17,8 5,2 MPa) > G2 (13,5 4,3 MPa). O terço do conduto não exerceu influência significativa sobre a resistência adesiva (p = 0,4749). As espessuras dos filmes de cimento foram estatisticamente diferentes nos diferentes terços: cervical (102 45 m) > médio (75 29 m) > apical (52 28m). Não foi observada forte correlação entre os valores de espessura e os de resistência ao push-out (r = - 0,2016, p = 0,0033). O tipo de falha predominante foi a mista, exceto para o G2, que apresentou 74% das falhas na interface cimento-pino. Dessa forma, o cimento autoadesivo apresentou melhor desempenho que o convencional, e ambos os sistemas duais, sobretudo o RelyX ARC, apresentaram dependência da fotoativação para atingirem maiores valores de resistência de união.

Paralelización mediante OpenMP de segmentación de imágenes para el análisis de materiales bidimensionales

A lo largo de este Proyecto de Fin de Carrera se ha implementado una versión paralela mediante OpenMP de una aproximación al algoritmo de seg- mentación de imágenes level set. El fin de esta segmentación es la extracción de ciertas características de las imágenes necesarias para el cliente con el que se trabaja a lo largo de este proyecto. Además, se ha realizado un estudio previo de los tipos de técnicas de segmentación existentes y se han realizado pruebas experimentales con la paralelización realizada, con el fin de analizar la mejora realizada frente al algoritmo en serie.

Computational fluid-structure interaction methods for simulation of inflatable aerodynamic decelerators

Inflatable aerodynamic decelerators have potential advantages for planetary re-entry in robotic and human exploration missions. It is theorized that volume-mass characteristics of these decelerators are superior to those of common supersonic/subsonic parachutes and after deployment they may suffer no instabilities at high Mach numbers. A high fidelity computational fluid-structure interaction model is employed to investigate the behavior of tension cone inflatable aeroshells at supersonic speeds up to Mach 2.0. The computational framework targets the large displacements regime encountered during the inflation of the decelerator using fast level set techniques to incorporate boundary conditions of the moving structure. The preliminary results indicate large but steady aeroshell displacement with rich dynamics, including buckling of the inflatable torus that maintains the decelerator open under normal operational conditions, owing to interactions with the turbulent wake. Copyright © 2009 by the American Institute of Aeronautics and Astronautics, Inc.

Towards topology-free optimization - An application to turbine internal cooling geometries

The application of automated design optimization to real-world, complex geometry problems is a significant challenge - especially if the topology is not known a priori like in turbine internal cooling. The long term goal of our work is to focus on an end-to-end integration of the whole CFD Process, from solid model through meshing, solving and post-processing to enable this type of design optimization to become viable & practical. In recent papers we have reported the integration of a Level Set based geometry kernel with an octree-based cut- Cartesian mesh generator, RANS flow solver, post-processing & geometry editing all within a single piece of software - and all implemented in parallel with commodity PC clusters as the target. The cut-cells which characterize the approach are eliminated by exporting a body-conformal mesh guided by the underpinning Level Set. This paper extends this work still further with a simple scoping study showing how the basic functionality can be scripted & automated and then used as the basis for automated optimization of a generic gas turbine cooling geometry. Copyright © 2008 by W.N.Dawes.

A practical demonstration of scalable, parallel mesh generation

Cambridge Flow Solutions Ltd, Compass House, Vision Park, Cambridge, CB4 9AD, UK Real-world simulation challenges are getting bigger: virtual aero-engines with multistage blade rows coupled with their secondary air systems & with fully featured geometry; environmental flows at meta-scales over resolved cities; synthetic battlefields. It is clear that the future of simulation is scalable, end-to-end parallelism. To address these challenges we have reported in a sequence of papers a series of inherently parallel building blocks based on the integration of a Level Set based geometry kernel with an octree-based cut-Cartesian mesh generator, RANS flow solver, post-processing and geometry management & editing. The cut-cells which characterize the approach are eliminated by exporting a body-conformal mesh driven by the underpinning Level Set and managed by mesh quality optimization algorithms; this permits third party flow solvers to be deployed. This paper continues this sequence by reporting & demonstrating two main novelties: variable depth volume mesh refinement enabling variable surface mesh refinement and a radical rework of the mesh generation into a bottom-up system based on Space Filling Curves. Also reported are the associated extensions to body-conformal mesh export. Everything is implemented in a scalable, parallel manner. As a practical demonstration, meshes of guaranteed quality are generated for a fully resolved, generic aircraft carrier geometry, a cooled disc brake assembly and a B747 in landing configuration. Copyright © 2009 by W.N.Dawes.

Finite volume distance field solution applied to medial axis transform

Accurate and efficient computation of the nearest wall distance d (or level set) is important for many areas of computational science/engineering. Differential equation-based distance/ level set algorithms, such as the hyperbolic-natured Eikonal equation, have demonstrated valuable computational efficiency. Here, in the context, as an 'auxiliary' equation to the main flow equations, the Eikonal equation is solved efficiently with two different finite volume approaches (the cell vertex and cell-centered). Application of the distance solution is studied for various geometries. Moreover, a procedure using the differential field to obtain the medial axis transform (MAT) for different geometries is presented. The latter provides a skeleton representation of geometric models that has many useful analysis properties. As an alternative approach to the pure geometric methods (e.g. the Voronoi approach), the current d-MAT procedure bypasses many difficulties that are usually encountered by pure geometric methods, especially in three dimensional space. It is also shown that the d-MAT approach provides the potential to sculpt/control the MAT form for specialized solution purposes. Copyright © 2010 by the American Institute of Aeronautics and Astronautics, Inc.

Parallel implementation of Multilevel BDDC

In application of the Balancing Domain Decomposition by Constraints (BDDC) to a case with many substructures, solving the coarse problem exactly becomes the bottleneck which spoils scalability of the solver. However, it is straightforward for BDDC to substitute the exact solution of the coarse problem by another step of BDDC method with subdomains playing the role of elements. In this way, the algorithm of three-level BDDC method is obtained. If this approach is applied recursively, multilevel BDDC method is derived. We present a detailed description of a recently developed parallel implementation of this algorithm. The implementation is applied to an engineering problem of linear elasticity and a benchmark problem of Stokes flow in a cavity. Results by the multilevel approach are compared to those by the standard (two-level) BDDC method.

Numerical simulation of a bubble rising in shear-thinning fluids

The motion of a single bubble rising freely in quiescent non-Newtonian viscous fluids was investigated experimentally and computationally. The non-Newtonian effects in the flow of viscous inelastic fluids are modeled by the Carreau theological model. An improved level set approach for computing the incompressible two-phase flow with deformable free interface is used. The control volume formulation with the SIMPLEC algorithm incorporated is used to solve the governing equations on a staggered Eulerian grid. The simulation results demonstrate that the algorithm is robust for shear-thinning liquids with large density (rho(1)/rho(g) up to 10(3)) and high viscosity (eta(1)/eta(g) up to 10(4)). The comparison of the experimental measurements of terminal bubble shape and velocity with the computational results is satisfactory. It is shown that the local change in viscosity around a bubble greatly depends on the bubble shape and the zero-shear viscosity of non-Newtonian shear-thinning liquids. The shear-rate distribution and velocity fields are used to elucidate the formation of a region of large viscosity at the rear of a bubble as a result of the rather stagnant flow behind the bubble. The numerical results provide the basis for further investigations, such as the numerical simulation of viscoelastic fluids. (C) 2010 Elsevier B.V. All rights reserved.

气泡/液滴Marangoni迁移数值模拟研究

本文利用Level Set方法,数值模拟了微重力情况下气泡/液滴的Marangoni迁移现象,分析了Marangoni数对迁移速度的影响.数值模拟结果表明,随着Marangoni数的增大,非线性热对流效应的影响会逐渐增大,导致沿相界面温度分布趋于均匀,从而降低迁移运动的驱动力,使气泡/液滴的迁移速度随Marangoni数的增加而逐渐减小.

单气泡池沸腾过程中的气泡动力学数值模拟研究

本文以单气泡池沸腾作为核态池沸腾现象的模型体系,以避免汽化核心分布的随机性及气泡间相互作用等引起的复杂性,深入研究核态池沸腾过程中的相变及气泡周围细观流动与传热规律。基于Level Set方法建立了描述单气泡池沸腾问题的数值模型,模拟了气泡的生长与脱落过程,并分析了气泡周围细观流动与传热。

液滴热毛细迁移流场和温度场拓扑结构的演变

热毛细迁移现象是流体颗粒(液滴/气泡)在非均匀温度场中由于界面温度梯度引起的非均匀界面张力驱动的运动。它不仅是流体力学中的经典问题之一,而且在诸如空间材料制备、空间流体和热管理系统等应用中也有着重要的应用。本文利用投影法求解了微重力条件下可变形液滴的轴对称热毛细迁移问题,控制方程组基于Level-set方法和连续界面张力模型。数值计算的终端迁移速度和空间实验结果相一致。计算结果表明,不同Marangoni数(Ma)情形具有相同的流场拓扑结构;但随着Ma数的增加,温度场的拓扑结构变化极大——在较小Ma数情形中,液滴内的最小温度发生在液滴尾部滞止点处。当Ma数超过某个临界值(10～20之间)后,最小温度点跳进液滴内部,并随着Ma数的增加而不断上移;液滴内部冷区最初呈球帽状,但其中部厚度随Ma数增加不断减小,同时向外扩展,形成外缘不断增厚的球壳状冷区;当Ma数超过另一临界值(约100)后,球壳状冷区在液滴轴线处破裂,冷区的拓扑结构转变为环状,最小温度点也随之远离轴线,不断向外、向下移动。温度场拓扑结构的变化反映了对流效应对液滴内部热量传输的影响不断增强,也强烈影响着液滴热毛细迁移终端速度随Ma数的变化特征。

The study on a PVC membrane electrode for gemfibrozil

In this paper, a poly(vinyl chloride) (PVC) membrane electrode is prepared for gemfibrozil, 2, 2-dimethyl-5(2,5-xylyloxy) valeric acid, based on its ion pair complexes with hexadecyltrioctyl ammonium iodide (HTOA). The membrane composition of the electrode was optimized by using the sequential level elimination method for orthogonal experimental design. The electrode has a Nernstian response range from 2.5 X 10(-5) to 0.1 mol/l with an average slope of 55.3 mV/decade. The limit of detection is 7.1 X 10(-6) mol/l. The electrode responses were not affected by pH in the range 10.0-12.3. A Na2B4O7-Na2CO3 buffer of pH = 11.0 was selected as the background electrolyte solution for potentiometric measurements. The electrode was used for determining gemfibrozil in pharmaceutical preparations with satisfactory results.