应用气液两相流模型模拟湖泊风生流

采用了基于两相流理论的三维相间界面计算模型模拟湖泊风生流,并采用交错网格的有限体积法求解。该模型应用于昆明市滇池的风生流模拟,实际计算表明该模型算法程序简单,计算结果合理,能较真实的反映湖泊风生流的流态特征。

微重力条件下火箭液氢贮箱重定位过程的数值模拟

采用流体体积法（VOF） 对微重力条件下火箭液氢贮箱的重定位过程进行了三维数值模拟研究。首先结合美国 半人马座液氢贮箱缩比模型的落塔微重力实验，针对无挡板情形，得到了推进剂触底、反弹、碰撞、触顶、顶部清 空和完成重定位六个特征时间，与实验结果符合较好。进一步考虑了环形防晃挡板的影响，得到了重定位时间和邦 德数（Bo）的依赖关系。发现在可靠沉底的前提下，较小的沉底邦德数时完成重定位所需要的总冲量较小，可以节 省推进剂用量。

重定位过程的数值模拟

<结合火箭液氢推进剂贮箱，采用流体体积法(VOF方法)对微重力条件下火箭液氢贮箱中推进剂的重定位过程进行数值模拟研究。针对自由面变形问题建立了完整的数学模型，考虑表面张力效应，研究在不同邦德数(Bo)下自由面变形过程的特点。结合美国半人马座液氢贮箱缩比模型的微重力实验结果，首先考虑无挡板情形，对推进剂重定位过程进行了三维直接数值模拟，得到了推进剂触底、反弹、碰撞、触顶和清空五个特征时间，与美国落塔微重力实验结果符合很好。进一步考虑了环形防晃挡板的影响，完成了相应的数值模拟研究。

分岔管路内的油水相分布特性研究

分岔管路作为一种新型管式分离部件已经在油田上得到了初步应用,与沉降罐和卧式分离器等传统分离设备相比具有分离效率高、体积紧凑、压降损失小等优点。不同入口条件和运行工况下,多相介质在分岔管路内的相分布特性对其分离性能和流动特性有显著的影响,目前这方面的研究还需深入和完善。本研究通过室内实验观察了分岔接头和分支管路内的相分布情况,采用VOF模型中的PLIC方法对分岔管路内的油水相界面进行了数值重构,计算中考虑了混合流速、相含率和流量配比等参数的影响,模拟结果与室内实验符合较好。计算发现,油相在分岔接头和分支管路内发生了明显的拉伸、破碎和变形等现象,在较低的混合流速和相含率情况下油水两相可以得到更为有效的分离,流量配比是影响相分布的最重要的运行参数,这些结论对于指导分岔管路这种新型管式油水分离设备的现场试验具有重要的意义。

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

Mesoscopic simulation of the impregnating process of unidirectional fibrous preform in resin transfer molding

The resin transfer molding has gained popularity in the preparation of fiber-reinforced polymer-matrix composites because of its high efficiency and low pollution. The non-uniform inter-tow and intra-tow flows are regarded as the reason of void formation in RTM. According to the process characteristics, the axisymmetric model was developed to study the interaction between the flow in the inter-tow space and that in the intra-tow space. The flow behavior inside the fiber tows was formulated using Brinkman's equation, while that in the open space around the fiber tows was formulated by Stokes' equation. The volume of fluid (VOF) method was applied to track the flow front, and the effects of filling velocity, resin viscosity, inter-tow dimension and intra-tow permeability on fluid pressure and flow front were analyzed. The results show that the flow front difference between the inter-tow and intra-tow becomes larger with the decrease of intra-tow permeability, as well as the increase of filling velocity and inter-tow dimension.

Numerical wave flume study on wave motion around submerged plates

Nonlinear interaction between surface waves and a submerged horizontal plate is investigated in the absorbed numerical wave flume developed based on the volume of fluid (VOF) method. The governing equations of the numerical model are the continuity equation and the Reynolds-Averaged Navier-Stokes (RANS) equations with the k-epsilon turbulence equations. Incident waves are generated by an absorbing wave-maker that eliminates the waves reflected from structures. Results are obtained for a range of parameters, with consideration of the condition under which the reflection coefficient becomes maximal and the transmission coefficient minimal. Wave breaking over the plate, vortex shedding downwave, and pulsating flow below the plate are observed. Time-averaged hydrodynamic force reveals a negative drift force. All these characteristics provide a reference for construction of submerged plate breakwaters.

A volume of fluid numerical model for wave impacts at coastal structures

The first stages in the development of a new design tool, to be used by coastal engineers to improve the efficiency, analysis, design, management and operation of a wide range of coastal and harbour structures, are described. The tool is based on a two-dimensional numerical model, NEWMOTICS-2D, using the volume of fluid (VOF) method, which permits the rapid calculation of wave hydrodynamics at impermeable natural and man-made structures. The critical hydrodynamic flow processes and forces are identified together with the equations that describe these key processes. The different possible numerical approaches for the solution of these equations, and the types of numerical models currently available, are examined and assessed. Preliminary tests of the model, using comparisons with results from a series of hydraulic model test cases, are described. The results of these tests demonstrate that the VOF approach is particularly appropriate for the simulation of the dynamics of waves at coastal structures because of its flexibility in representing the complex free surfaces encountered during wave impact and breaking. The further programme of work, required to develop the existing model into a tool for use in routine engineering design, is outlined.

Numerical Study of Wave Slamming on an Oscillating Flap

Bottom hinged Oscillating Wave Surge Converters (OWSCs) are efficient devices for extracting power from ocean waves. There is limited knowledge about wave slamming on such devices. This paper deals with numerical studies of wave slamming on an oscillating flap to investigate the mechanism of slamming events. In our model, the Navier–Stokes equations are discretized using the Finite Volume method with the Volume of Fluid (VOF) approach for interface capturing. Waves are generated by a flaptype wave maker in the numerical wave tank, and the dynamic mesh method is applied to model the motion of the oscillating flap. Basic mesh and time step refinement studies are performed. The flow characteristics in a slamming event are analysed based on numerical results. Various simulations with different flap densities, water depths and wave amplitudes are performed for a better understanding of the slamming.

CFD Simulation of Multiphase Reactors

Regional Research Laboratory

Numerical Experimental Comparison of Dam Break Flows with non-Newtonian Fluids

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Different prescribed doses of high-volume peritoneal dialysis and outcome of patients with acute kidney injury.

The optimal dialysis dose for the treatment of acute kidney injury (AKI) is controversial. No studies have directly examined the effects of peritoneal dialysis (PD) dose on outcomes in AKI. From January 2005 to January 2007, we randomly assigned critically ill patients with AKI to receive higher- or lower-intensity PD therapy (prescribed Kt/Vof 0.8 and 0.5 per session respectively). The main outcome measure was death within 30 days. Of the 61 enrolled patients, 30 were randomly assigned to higher-intensity therapy, and 31, to a lower-intensity PD dose. The two study groups had similar baseline characteristics and received treatment for 6.1 days and 5.7 days respectively (p = 0.42). At 30 days after randomization, 17 deaths had occurred in the higher-intensity group (55%), and 16 deaths, in the lower-intensity group (53%, p = 0.83). There was a significant difference between the groups in the PD dose prescribed compared with the dose delivered (higher-intensity group: 0.8 vs. 0.59, p = 0.04; lower-intensity group: 0.5 vs. 0.49, p = 0.89). The groups had similar metabolic control after 4 PD sessions (blood urea nitrogen: 69.3 +/- 14.4 mg/dL and 60.3 +/- 11.1 mg/dL respectively, p = 0. 71). In critically ill patients with AKI, an intensive PD dose did not lower the mortality or improve the recovery of kidney function or metabolic control. The PD dose is limited by dialysate flow and membrane permeability, and clearance per exchange can decrease if a shorter dwell time is applied.

Rompimento de barreiras: análise experimental e numérica na previsão de velocidades de propagação de frentes de material hiperconcentrado

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Hydrodynamic forces on heave plates for offshore systems oscillating close to the seabed or the free surface

La energía eólica marina es uno de los recursos energéticos con mayor proyección pudiendo contribuir a reducir el consumo de combustibles fósiles y a cubrir la demanda de energía en todo el mundo. El concepto de aerogenerador marino está basado en estructuras fijas como jackets o en plataformas flotantes, ya sea una semisumergible o una TLP. Se espera que la energía eólica offshore juegue un papel importante en el perfil de producción energética de los próximos años; por tanto, las turbinas eólicas deben hacerse más fables y rentables para ser competitivas frente a otras fuentes de energía. Las estructuras flotantes pueden experimentar movimientos resonantes en estados de la mar con largos períodos de oleaje. Estos movimientos disminuyen su operatividad y pueden causar daños en los componentes eléctricos de las turbinas y en las palas, también en los risers y moorings. La respuesta de la componente vertical del movimiento puede reducirse mediante diferentes actuaciones: (1) aumentando la amortiguación del sistema, (2) manteniendo el período del movimiento vertical fuera del rango de la energía de la ola, y (3) reduciendo las fuerzas de excitación verticales. Un ejemplo típico para llevar a cabo esta reducción son las "Heave Plates". Las heave plates son placas que se utilizan en la industria offshore debido a sus características hidrodinámicas, ya que aumentan la masa añadida y la amortiguación del sistema. En un análisis hidrodinámico convencional, se considera una estructura sometida a un oleaje con determinadas características y se evalúan las cargas lineales usando la teoría potencial. El amortiguamiento viscoso, que juega un papel crucial en la respuesta en resonancia del sistema, es un dato de entrada para el análisis. La tesis se centra principalmente en la predicción del amortiguamiento viscoso y de la masa añadida de las heave plates usadas en las turbinas eólicas flotantes. En los cálculos, las fuerzas hidrodinámicas se han obtenido con el f n de estudiar cómo los coeficientes hidrodinámicos de masa añadida5 y amortiguamiento varían con el número de KC, que caracteriza la amplitud del movimiento respecto al diámetro del disco. Por otra parte, se ha investigado la influencia de la distancia media de la ‘heave plate’ a la superficie libre o al fondo del mar, sobre los coeficientes hidrodinámicos. En este proceso, un nuevo modelo que describe el trabajo realizado por la amortiguación en función de la enstrofía, es descrito en el presente documento. Este nuevo enfoque es capaz de proporcionar una correlación directa entre el desprendimiento local de vorticidad y la fuerza de amortiguación global. El análisis también incluye el estudio de los efectos de la geometría de la heave plate, y examina la sensibilidad de los coeficientes hidrodinámicos al incluir porosidad en ésta. Un diseño novedoso de una heave plate, basado en la teoría fractal, también fue analizado experimentalmente y comparado con datos experimentales obtenidos por otros autores. Para la resolución de las ecuaciones de Navier Stokes se ha usado un solver basado en el método de volúmenes finitos. El solver usa las librerías de OpenFOAM (Open source Field Operation And Manipulation), para resolver un problema multifásico e incompresible, usando la técnica VOF (volume of fluid) que permite capturar el movimiento de la superficie libre. Los resultados numéricos han sido comparados con resultados experimentales llevados a cabo en el Canal del Ensayos Hidrodinámicos (CEHINAV) de la Universidad Politécnica de Madrid y en el Canal de Experiencias Hidrodinámicas (CEHIPAR) en Madrid, al igual que con otros experimentos realizados en la Escuela de Ingeniería Mecánica de la Universidad de Western Australia. Los principales resultados se presentan a continuación: 1. Para pequeños valores de KC, los coeficientes hidrodinámicos de masa añadida y amortiguamiento incrementan su valor a medida que el disco se aproxima al fondo marino. Para los casos cuando el disco oscila cerca de la superficie libre, la dependencia de los coeficientes hidrodinámicos es más fuerte por la influencia del movimiento de la superficie libre. 2. Los casos analizados muestran la existencia de un valor crítico de KC, donde la tendencia de los coeficientes hidrodinámicos se ve alterada. Dicho valor crítico depende de la distancia al fondo marino o a la superficie libre. 3. El comportamiento físico del flujo, para valores de KC cercanos a su valor crítico ha sido estudiado mediante el análisis del campo de vorticidad. 4. Introducir porosidad al disco, reduce la masa añadida para los valores de KC estudiados, pero se ha encontrado que la porosidad incrementa el valor del coeficiente de amortiguamiento cuando se incrementa la amplitud del movimiento, logrando un máximo de damping para un disco con 10% de porosidad. 5. Los resultados numéricos y experimentales para los discos con faldón, muestran que usar este tipo de geometrías incrementa la masa añadida cuando se compara con el disco sólido, pero reduce considerablemente el coeficiente de amortiguamiento. 6. Un diseño novedoso de heave plate basado en la teoría fractal ha sido experimentalmente estudiado a diferentes calados y comparado con datos experimentales obtenidos por otro autores. Los resultados muestran un comportamiento incierto de los coeficientes y por tanto este diseño debería ser estudiado más a fondo. ABSTRACT Offshore wind energy is one of the promising resources which can reduce the fossil fuel energy consumption and cover worldwide energy demands. Offshore wind turbine concepts are based on either a fixed structure as a jacket or a floating offshore platform like a semisubmersible, spar or tension leg platform. Floating offshore wind turbines have the potential to be an important part of the energy production profile in the coming years. In order to accomplish this wind integration, these wind turbines need to be made more reliable and cost efficient to be competitive with other sources of energy. Floating offshore artifacts, such oil rings and wind turbines, may experience resonant heave motions in sea states with long peak periods. These heave resonances may increase the system downtime and cause damage on the system components and as well as on risers and mooring systems. The heave resonant response may be reduced by different means: (1) increasing the damping of the system, (2) keeping the natural heave period outside the range of the wave energy, and (3) reducing the heave excitation forces. A typical example to accomplish this reduction are “Heave Plates”. Heave plates are used in the offshore industry due to their hydrodynamic characteristics, i.e., increased added mass and damping. Conventional offshore hydrodynamic analysis considers a structure in waves, and evaluates the linear and nonlinear loads using potential theory. Viscous damping, which is expected to play a crucial role in the resonant response, is an empirical input to the analysis, and is not explicitly calculated. The present research has been mainly focused on the prediction of viscous damping and added mass of floating offshore wind turbine heave plates. In the calculations, the hydrodynamic forces have been measured in order to compute how the hydrodynamic coefficients of added mass1 and damping vary with the KC number, which characterises the amplitude of heave motion relative to the diameter of the disc. In addition, the influence on the hydrodynamic coefficients when the heave plate is oscillating close to the free surface or the seabed has been investigated. In this process, a new model describing the work done by damping in terms of the flow enstrophy, is described herein. This new approach is able to provide a direct correlation between the local vortex shedding processes and the global damping force. The analysis also includes the study of different edges geometry, and examines the sensitivity of the damping and added mass coefficients to the porosity of the plate. A novel porous heave plate based on fractal theory has also been proposed, tested experimentally and compared with experimental data obtained by other authors for plates with similar porosity. A numerical solver of Navier Stokes equations, based on the finite volume technique has been applied. It uses the open-source libraries of OpenFOAM (Open source Field Operation And Manipulation), to solve 2 incompressible, isothermal immiscible fluids using a VOF (volume of fluid) phase-fraction based interface capturing approach, with optional mesh motion and mesh topology changes including adaptive re-meshing. Numerical results have been compared with experiments conducted at Technical University of Madrid (CEHINAV) and CEHIPAR model basins in Madrid and with others performed at School of Mechanical Engineering in The University of Western Australia. A brief summary of main results are presented below: 1. At low KC numbers, a systematic increase in added mass and damping, corresponding to an increase in the seabed proximity, is observed. Specifically, for the cases when the heave plate is oscillating closer to the free surface, the dependence of the hydrodynamic coefficients is strongly influenced by the free surface. 2. As seen in experiments, a critical KC, where the linear trend of the hydrodynamic coefficients with KC is disrupted and that depends on the seabed or free surface distance, has been found. 3. The physical behavior of the flow around the critical KC has been explained through an analysis of the flow vorticity field. 4. The porosity of the heave plates reduces the added mass for the studied porosity at all KC numbers, but the porous heave plates are found to increase the damping coefficient with increasing amplitude of oscillation, achieving a maximum damping coefficient for the heave plate with 10% porosity in the entire KC range. 5. Another concept taken into account in this work has been the heave plates with flaps. Numerical and experimental results show that using discs with flaps will increase added mass when compared to the plain plate but may also significantly reduce damping. 6. A novel heave plate design based on fractal theory has tested experimentally for different submergences and compared with experimental data obtained by other authors for porous plates. Results show an unclear behavior in the coefficients and should be studied further. Future work is necessary in order to address a series of open questions focusing on 3D effects, optimization of the heave plates shapes, etc.