摄动有限体积算法及其在两相流中应用

采用一种全新的摄动有限体积(PFV)算法和水平集(Level Set)技术对液液两相系统中液滴坠落进行数值模拟,数值结果表明,PFV新算法具有节点少、精度高,效率高,编程方便等优点,能成功模拟液液两相流动,为两相流动数值模拟提供了一种新的途径.

Numerical simulation of condensation of bubbles under microgravity conditions by moving mesh method in the double-staggered grid

A numerical 2D method for simulation of two-phase flows including phase change under microgravity conditions is presented in this paper, with a level set method being coupled with the moving mesh method in the double-staggered grid systems. When the grid lines bend very much in a curvilinear grid, great errors may be generated by using the collocated grid or the staggered grid. So the double-staggered grid was adopted in this paper. The level set method is used to track the liquid-vapor interface. The numerical analysis is fulfilled by solving the Navier-Stokes equations using the SIMPLER method, and the surface tension force is modeled by a continuum surface force approximation. A comparison of the numerical results obtained with different numerical strategies shows that the double-staggered grid moving-mesh method presented in this paper is more accurate than that used previously in the collocated grid system. Based on the method presented in this paper, the condensation of a single bubble in the cold water under different level of gravity is simulated. The results show that the condensation process under the normal gravity condition is different from the condensation process under microgravity conditions. The whole condensation time is much longer under the normal gravity than under the microgravity conditions.

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.

Microscopic Origin of Macroscopic Strength in Granular Media: A Numerical and Analytical Approach

Constitutive modeling in granular materials has historically been based on macroscopic experimental observations that, while being usually effective at predicting the bulk behavior of these type of materials, suffer important limitations when it comes to understanding the physics behind grain-to-grain interactions that induce the material to macroscopically behave in a given way when subjected to certain boundary conditions.

The advent of the discrete element method (DEM) in the late 1970s helped scientists and engineers to gain a deeper insight into some of the most fundamental mechanisms furnishing the grain scale. However, one of the most critical limitations of classical DEM schemes has been their inability to account for complex grain morphologies. Instead, simplified geometries such as discs, spheres, and polyhedra have typically been used. Fortunately, in the last fifteen years, there has been an increasing development of new computational as well as experimental techniques, such as non-uniform rational basis splines (NURBS) and 3D X-ray Computed Tomography (3DXRCT), which are contributing to create new tools that enable the inclusion of complex grain morphologies into DEM schemes.

Yet, as the scientific community is still developing these new tools, there is still a gap in thoroughly understanding the physical relations connecting grain and continuum scales as well as in the development of discrete techniques that can predict the emergent behavior of granular materials without resorting to phenomenology, but rather can directly unravel the micro-mechanical origin of macroscopic behavior.

In order to contribute towards closing the aforementioned gap, we have developed a micro-mechanical analysis of macroscopic peak strength, critical state, and residual strength in two-dimensional non-cohesive granular media, where typical continuum constitutive quantities such as frictional strength and dilation angle are explicitly related to their corresponding grain-scale counterparts (e.g., inter-particle contact forces, fabric, particle displacements, and velocities), providing an across-the-scale basis for better understanding and modeling granular media.

In the same way, we utilize a new DEM scheme (LS-DEM) that takes advantage of a mathematical technique called level set (LS) to enable the inclusion of real grain shapes into a classical discrete element method. After calibrating LS-DEM with respect to real experimental results, we exploit part of its potential to study the dependency of critical state (CS) parameters such as the critical state line (CSL) slope, CSL intercept, and CS friction angle on the grain's morphology, i.e., sphericity, roundness, and regularity.

Finally, we introduce a first computational algorithm to ``clone'' the grain morphologies of a sample of real digital grains. This cloning algorithm allows us to generate an arbitrary number of cloned grains that satisfy the same morphological features (e.g., roundness and aspect ratio) displayed by their real parents and can be included into a DEM simulation of a given mechanical phenomenon. In turn, this will help with the development of discrete techniques that can directly predict the engineering scale behavior of granular media without resorting to phenomenology.

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.

A fixed-grid b-spline finite element technique for fluid-structure interaction

We present a fixed-grid finite element technique for fluid-structure interaction problems involving incompressible viscous flows and thin structures. The flow equations are discretised with isoparametric b-spline basis functions defined on a logically Cartesian grid. In addition, the previously proposed subdivision-stabilisation technique is used to ensure inf-sup stability. The beam equations are discretised with b-splines and the shell equations with subdivision basis functions, both leading to a rotation-free formulation. The interface conditions between the fluid and the structure are enforced with the Nitsche technique. The resulting coupled system of equations is solved with a Dirichlet-Robin partitioning scheme, and the fluid equations are solved with a pressure-correction method. Auxiliary techniques employed for improving numerical robustness include the level-set based implicit representation of the structure interface on the fluid grid, a cut-cell integration algorithm based on marching tetrahedra and the conservative data transfer between the fluid and structure discretisations. A number of verification and validation examples, primarily motivated by animal locomotion in air or water, demonstrate the robustness and efficiency of our approach. © 2013 John Wiley & Sons, Ltd.

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.

Thermocapillary Migration of Deformable Bubbles at Moderate to Large Marangoni Number in Microgravity

Using the level-set method and the continuum interface model, the axisymmetric thermocapillary migration of gas bubbles in an immiscible bulk liquid with a temperature gradient at moderate to large Marangoni number is simulated numerically. Constant material properties of the two phases are assumed. Steady state of the motion can always be reached. The terminal migration velocity decreases monotonously with the increase of the Marangoni number due to the wrapping of isotherms around the front surface of the bubble. Good agreements with space experimental data and previous theoretical and numerical studies in the literature are evident. Slight deformation of bubble is observed, but no distinct influence on the motion occurs. It is also found that the influence of the convective transport of heat inside bubbles cannot be neglected at finite Marangoni number, while the influence of the convective transport of momentum inside bubbles may be actually negligible.

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

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

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

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

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

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