Numerical simulation of droplet ejection of thermal inkjet printheads

In this study, we present a method to predict the droplet ejection in thermal inkjet printheads including the growth and collapse of a vapor bubble and refill of the firing chamber. The three‐dimensional Navier–Stokes equations are solved using a finite‐volume approach with a fixed Cartesian mesh. The piecewise‐linear interface calculation‐based volume‐of‐fluid method is employed to track and reconstruct the ink–air interface. A geometrical computation based on Lagrangian advection is used to compute the mass flux and advance the interface. A simple and efficient model for the bubble dynamics is employed to model the effect of ink vapor on the adjacent ink liquid. To solve the surface tension‐dominated flow accurately, a hierarchical curvature‐estimation method is proposed to adapt to the local grid resolution. The numerical methods mentioned earlier have been implemented in an internal simulation code, CFD3. The numerical examples presented in the study show good performance of CFD3 in prediction of surface tension‐dominated free‐surface flows, for example, droplet ejection in thermal inkjet printing. Currently, CFD3 is used extensively for printhead development within Hewlett‐Packard. Copyright © 2015 John Wiley & Sons, Ltd.     

Assessment of a high‐order MUSCL method for rotor flows

This work presents the implementation of a high‐order, finite‐volume scheme suitable for rotor flows. The formulation is based on the variable extrapolation MUSCL‐scheme, where high‐order spatial accuracy (up to fourth‐order) is achieved using correction terms obtained through successive differentiation. A variety of results are presented, including 2‐ and 3‐dimensional test cases. Results with the proposed scheme, showed better wake and higher resolution of vortical structures compared with the standard MUSCL, even when coarse meshes were employed. The method was also demonstrated for 3‐dimensional unsteady flows using overset and moving grids for the UH‐60A rotor in forward flight and the Enhanced Rotorcraft Innovative Concept Achievement tiltrotor in aeroplane mode. For medium grids, the present method adds reasonable CPU and memory overheads and offers good accuracy on relatively coarse grids.     

Design and optimization of a cavitating device for Congo red decolorization: Experimental investigation and CFD simulation

The aim of this work is to perform design and optimization of a cavitating device based on CFD simulation. A set of operational and geometrical parameters such as convergence angle, divergence angle, length of throat, and inlet pressure that can affect the hydrodynamic cavitation phenomenon generating in a Venturi are evaluated through CFD simulation and experimental approaches. Response surface methodology (RSM) was employed to achieve the optimum geometrical configuration. The CFD results show that the maximum cavitation zone in the Venturi can be obtained when half angle of the convergence section, throat length and half angle of the divergence section are 22.7°, 4 mm, and 6.5°, respectively. A maximum decolorization of 38.8% has been obtained using the designed Venturi at cavitation number (Cv) of 0.12. Additionally, the results were compared to that of various orifice plates. A decolorization of 26.2% using 33 holes orifice plate and 11.55% in one hole orifice plate approved the superiority of the designed Venturi.     

Analysis of droplet behavior in a de-oiling hydrocyclone

A mechanical separation process in a de-oiling hydrocyclone is described in which disperse oil droplets are separated from a continuous water phase. This separation process is influenced by droplet breakage and coalescence. Based on experimental data and simulation results in a stirred tank, a modified breakage model, which can be applied to droplet breakage in the de-oiling hydrocyclone, is developed. Then, a simulation model is developed coupling the numerical solution of the flow field in the hydrocyclone based on computational fluid dynamics (CFD) with population balances. The homogenous discrete method and the inhomogeneous discrete method are applied for solving the population balance model (PBM). The investigations show that the numerical results obtained by the simulation model coupled with the modified PBM using the inhomogeneous discrete method are in good accordance with experimental data under a high flow rate. According to this simulation model, the effect of three different inlet designs on the separation efficiency of the de-oiling hydrocyclone has been discussed. The results indicate that the separation efficiency of the de-oiling hydrocyclone can be improved with an appropriate inlet design.     

Coupling of CFD Simulations and Population Balance Modeling to Predict Brownian Coagulation in an Emulsion Polymerization Reactor

In this study a framework consisting of a computational fluid dynamics simulation coupled to a population balance model for the modeling of emulsion polymerizations is proposed. The combined approach is used to understand the impact of changing length and time scales, as well as mixing conditions on the particle size distribution (PSD) of a polymer latex under different conditions. It is shown that the effect of agitation rate can have a profound impact on the distribution of ionic species in the reactor, and thus on the evolution of the PSD.

     

建筑风压场数值模拟中压力波动的平抑

在以同位网格为基础的简单流场压力计算中,通常采用动量插值方法来平抑流场中的压力波动现象;但是对于建筑风场等复杂的钝体绕流问题,由该平抑方法得到的收敛风压场仍可能存在小幅波动。为彻底解决同位网格格式下的压力波动,除采用动量插值方法外,本文提出了在压力校正方程的界面流速中添加压力梯度差值项的方法。算例分析表明,该方法计算得到的建筑风压场完全避免了压力波动现象,风压解与试验结果吻合良好。     

Fluid Flows Through Two-Dimensional Channels of Composite Materials

The present analysis relates to the study of the full two-dimensional Brinkman equation representing the fluid flow through porous medium. The steady, incompressible fluid flow, with a negligible gravitational force, is constrained to flow in an infinitely long channel in which the height assumes a series of piecewise constant values. The control volume method is used to solve the Brinkman equation which involves the parameter, =/Da, where Da is the Darcy number and is the ratio of the fluid viscosity _f to the effective viscosity . An analytical study in the fully developed section of the composite channel is presented when the channel is of constant height and composed of several layers of porous media, each of uniform porosity. In the fully developed flow regime the analytical and numerical solutions are graphically indistinguishable. A geometrical configuration involving several discontinuities of channel height, and where the entry and exit sections are layered, is considered and the effect of different permeabilities is demonstrated. Further, numerical investigations are performed to evaluate the behaviour of fluid flow through regions which mathematically model some geological structures of various sizes, positions and permeability, for example a fault or a fracture, where the outlet channel is offset at different levels. The effect on the overall pressure gradient is also considered.     

CFD SIMULATION OF FLUID CATALYTIC CRACKING IN DOWNER REACTORS

A mathematical model has been developed for the simulation of gas-particle flow and fluid catalytic cracking in downer reactors. The model takes into account both cracking reaction and flow behavior through a four-lump reaction kinetics coupled with two-phase turbulent flow. The prediction results show that the relatively large change of gas velocity affects directly the axial distribution of solids velocity and void fraction, which significantly interact with the chemical reaction. Furthermore, model simulations are carried out to determine the effects of such parameters on product yields, as bed diameter, reaction temperature and the ratio of catalyst to oil, which are helpful for optimizing the yields of desired products. The model equations are coded and solved on CFX4.4.     

CFD SIMULATION OF FLUID CATALYTIC CRACKING IN DOWNER REACTORS

A mathematical model has been developed for the simulation of gas-particle flow and fluid catalytic cracking in downer reactors. The model takes into account both cracking reaction and flow behavior through a four-lump reaction kinetics coupled with two-phase turbulent flow. The prediction results show that the relatively large change of gas velocity affects directly the axial distribution of solids velocity and void fraction, which significantly interact with the chemical reaction. Furthermore, model simulations are carried out to determine the effects of such parameters on product yields, as bed diameter, reaction temperature and the ratio of catalyst to oil, which are helpful for optimizing the yields of desired products. The model equations are coded and solved on CFX4.4.     

飞行器跨声速气动弹性数值分析

将流体和结构运动方程分别构造为含子迭代的计算格式,发展了一种紧耦合气动弹性分析方法.其中流体计算的空间离散采用改进的HLLEW(Harten—Lax-van Leer-Einfeldt-Wada)格式. TFI(transfinite inter- polation)方法用于生成随结构变形的自适应多块动网格.利用所发展的方法,对-翼-身-尾气动外形,数值预测了马赫数在0.3-1.3范围内的气动颤振边界.并详细研究了时间步长、子迭代步数、初始流场、耦合方法、疏密网格对颤振计算结果的影响.