Journal of Nanotechnology

Volume 2015 (2015), Article ID 642014, 10 pages

http://dx.doi.org/10.1155/2015/642014

## Growth Mechanisms of Nanostructured Titania in Turbulent Reacting Flows

Department of Mechanical Engineering, University of Minnesota, 111 Church Street SE, Minneapolis, MN 55455-0111, USA

Received 5 June 2015; Accepted 22 July 2015

Academic Editor: Lin-Hua Xu

Copyright © 2015 Sean C. Garrick. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

Titanium dioxide (titania) is used in chemical sensors, pigments, and paints and holds promise as an antimicrobial agent. This is due to its photoinduced activity and, in nanostructured form, its high specific surface area. Particle size and surface area result from the interplay of fluid, chemical, and thermal dynamics as well as nucleation, condensation and coagulation. After nucleation, condensation, and coagulation are the dominant phenomena affecting the particle size distribution. Manufacture of nanostructured titania via gas-phase synthesis often occurs under turbulent flow conditions. This study examines the competition between coagulation and condensation in the growth of nanostructured titania. Direct numerical simulation is utilized in simulating the hydrolysis of titanium tetrachloride to produce titania in a turbulent, planar jet. The fluid, chemical, and particle fields are resolved as a function of space and time. As a result, knowledge of titania is available as a function of space, time, and phase (vapor or particle), facilitating the analysis of the particle dynamics by mechanism. Results show that in the proximal region of the jet nucleation and condensation are the dominant mechanisms. However once the jet potential core collapses and turbulent mixing begins, coagulation is the dominant mechanism. The data also shows that the coagulation growth-rate is as much as twice the condensation growth-rate.

#### 1. Introduction

High rate synthesis of nanoparticles from vapor requires operation in the turbulent flow regime [1]. The variety of length and time scales present in turbulent multiphase flows makes them not very amenable to physical observation or analysis. This is especially true in the early stages of the nanoparticle formation and growth processes [2]. Computational fluid dynamics (CFD) has been developed for many years and a variety of numerical techniques have been developed and utilized for studying particle dynamics [3–10]. CFD has enabled engineers to achieve their goals more rapidly and cost effectively [11]. It has also become an effective tool for understanding physicochemical dynamics.

Nakaso et al. [12] modeled the titania nanoparticle growth of both agglomerates and primary particles simultaneously by using spatial zero-dimensional fluid dynamics. Johannessen et al. [13] combined CFD with a mathematical model for the particle dynamics to compare with experimental data from the synthesis of titanium dioxide particles in diffusion flames. Tsantilis et al. [5] used a moving sectional aerosol dynamics model accounting for gas-phase chemical reactions, coagulation, surface growth, and sintering with zero-dimensional fluid dynamics to investigate flame synthesis of titania nanoparticles. Also one study by Moody and Collins [14] considered titania nanoparticle nucleation and growth in a turbulent “box” located near the center of the reactor via three-dimensional DNS coupled with moment method in particle dynamics. Researchers have also considered the effects of turbulence on particle growth. Strakey et al. [15] studied the role of turbulence on the characteristics of powder made by oxidation and they found that the increased turbulence intensity narrowed the size distribution of the product powder indicating that the particle growth may be dominated by reactant mixing rather than by particle-particle collisions. These computations are quite compute intensive as they resolve all of the appropriate length and time scales. Garrick and his group studied the effects of turbulence on particle coagulation intensively and indicated that turbulence has a positive effect on the particle growth [16, 17]. The effects of turbulence on nucleation have also been elucidated [18–20]. To obtain physically accurate data, simulations must be three-dimensional and must be model-free. That is, the results must be obtained without the use of turbulence or subgrid scale models [21]. A review of the literature reveals a lack of the detailed information, especially via high-resolution direct simulation, on the interplay between the different mechanisms affecting nanoparticle growth.

In this work, the hydrolysis of titanium tetrachloride () to produce titanium dioxide () is simulated via direct numerical simulation. The turbulent, reacting, multiphase flow is obtained by solving the Navier-Stokes equations in conjunction with transport equations for all of the relevant chemical species and a nodal approach is used to represent the particle field [22, 23]. With the chemical and particle fields available as a function of space, time, and size, particle nucleation, condensation, and coagulation are illustrated individually. Additionally the two growth mechanisms are elucidated.

#### 2. Methodology

The mass, momentum, and energy equations are solved to obtain the fluid velocity , pressure , density , and the enthalpy . These variables are governed by the following conservation equations:where is the viscous stress tensor for a Newtonian fluid, is the coefficient of thermal conduction, and is the specific heat at constant pressure.

##### 2.1. Chemical Transport

The fluid contains five chemical species, the transport of which is given by the conservation of species equations:where is the mass concentration of species and is the diffusion coefficient of species . The reactants, and water vapor , undergo an irreversible one-step, isothermal chemical reaction at a temperature of 300 K and atmospheric pressure, which produces titanium dioxide and hydrochloric acid :The source term in (2) represents the effects of chemical reaction, the rate of creation or consumption of species The system is closed with the ideal gas equation of state, and , where and are the gas constant and molecular weight of species , respectively. The fluid temperature, , is obtained using the enthalpy via The fifth chemical species, nitrogen (), does not participate in any chemical reactions and serves as the carrier gas.

##### 2.2. Particle Field

The aerosol general dynamic equation (GDE) describes particle dynamics under the influence of various physicochemical phenomena, convection, diffusion, coagulation, surface growth, nucleation, and other internal/external forces. The GDE is utilized in discrete form as a population balance on each cluster or particle size. The methodology uses the nodal/sectional method of to approximate the GDE [23–26]. This approach effectively divides the aerosol population into three classes, monomers, clusters, and particles [16]. The GDE is therefore solved as a set of transport equations, one for each bin , [27]. Monomers of size nm in diameter populate bin while bins and are populated by clusters of molecules. Molecular clusters of size nm and larger are considered “particles” [28]. The general transport equation for the concentration of monomers, clusters, and particles in bin , , is written aswhere is the diffusivity given bywhere is the Boltzmann constant, is the Cunningham correction factor, and is the mean volume particle diameter [29]. The source term represents particle formation and growth processes and is given bywhere accounts for the formation of monomers in bin via chemical reaction [30–33] and is given byThe source term in (6), , represents the effects of nucleation, condensation (via monomer-cluster, monomer-particle and cluster-particle collisions), and Brownian coagulation. In this study, all particle formation and growth processes occur at 300 K. This means that all particles considered are agglomerates, meaning growth via coagulation results in fractal-like aggregate particles, consisting of primary particles (monomers). The collision frequency function, , is well documented in both the free-molecular and continuum regimes [29, 34–37]. Collisions of all particles (monomers, clusters, etc.) are considered and is given bywherewhere is the primary particle volume, is the volume of an agglomerate in the th section, is the particle density, and are the primary particle numbers in bins and , respectively, and is the fractal dimension. The particles generated by the industrial aerosol processes are usually nonspherical and the products are composed of groups of adhering particles ranging from loosely linked particles, sometimes termed agglomerates, to strongly necked particles, called aggregates or hard agglomerates. The fractal dimension is used to represent the shape of the aggregate or particle structure. A statistical concept, the fractal (or Hausdorff) dimension , was introduced to describe the shape of the agglomerate, which is obtained after averaging over many agglomerates with the same number of primary particles. The value of the fractal dimension ranges from to depending on the details of the agglomerate formation process [29]. The bins are organized in such a manner that the volume of particles in two successive bins is doubled; that is, [23, 38].

#### 3. Results

##### 3.1. Flow Configuration

The flow under consideration is a three-dimensional, isothermal, turbulent reacting jet issuing from an orifice of diameter into a coflowing stream. The jet is composed of titanium tetrachloride diluted in nitrogen , while the coflowing stream is composed of water vapor diluted in . The initial velocities are m/s for the jet and for the coflowing stream. The fluid field is characterized by the velocity ratio . The Reynolds number based on velocity of the high-speed stream and the jet diameter is = 3,000. The simulation is performed at 300 K and atmosphere pressure (1atm). To accelerate the development of large-scale structures, random perturbations with a maximum intensity of 5% are added to the cross-stream -velocity. The chemical composition of the jet stream is 0.01% TiCl_{4} and 99.99% N_{2} by mass. The simulation utilizes stoichiometric mixtures and the molar ratio of 1 : 2 for and .

In this work, several assumptions and approximations are utilized. These are stated below for clarity.(1)Ceramic powders such as have low equilibrium vapor pressures implying that single molecules may be considered particles [31]. As is formed, it appears as 0.5 nm diameter free spherical “monomers,” populating bin ; hence nucleation is treated as an instantaneous process.(2)The nanoparticles are small enough to follow the fluid path lines. Additionally, the particle volume fraction is of order . As a result the presence of the particles does not affect the fluid field.(3)Condensation is dominated by the collision rate of monomers with other monomers, clusters, and particles. This means the condensable species in the simulations is the “vapor."(4)The clusters and particles are stable because of the high supersaturation of the monomers. As a result, at these temperatures, there is no evaporation or sublimation of particles back to the gas-phase.(5)The fractal dimension, , is based on the collision time, , and sintering time, . In this work, all processes occur at K. At this temperature, particles do not sinter as [39].(6)A fractal dimension of is used when collisions occur between monomers while is used when monomers and dimers collide. For all other particle interactions (collisions between larger particles) a fractal dimension of is used [39, 40].

##### 3.2. Numerical Specifications

Ten bins are used to discretize the particle field (). The computational domain is of size and is comprised of grid points in the -, -, and -directions, respectively. The governing transport equations representing both the fluid and particle fields are solved using a MacCormack-based finite difference scheme [41, 42]. The scheme is of second order, accurate in time, and of fourth order, accurate in space. The boundary conditions are periodic in the spanwise -direction and zero-derivative in the cross-stream -direction and nonreflecting boundary conditions are used in both inflow and outflow boundaries (-direction) [43]. The simulation is performed up to a nondimensional time of which corresponds to a physical time of 0.48 ms. Both instantaneous and mean or averaged data are presented. We average in the -direction as it is the spatially homogeneous direction in planar jets. Quantities such as contours, isosurfaces, and mean data are useful in making qualitative and quantitative assessments of the nanoparticle growth dynamics as well as the underlying fluid and particle fields.

##### 3.3. Flow Field

The vorticity is the curl of the velocity vector and is an indicator of fluid mixing. The vorticity magnitude is the local rate of rotation. In nonpremixed chemically reacting flows, vorticity has the effect of increasing the interfacial area between the reactants. An isosurface of the instantaneous vorticity magnitude, the level-set, is shown in Figure 1 at time . The image shows that the flow is initially laminar and becomes turbulent as the jet travels downstream. Near , the two boundary layers initially located at merge and the jet spreads across the domain. This is aided by the presence of vortex braids, the tubular structures oriented in the stream-wise -direction, which act to draw the surrounding fluid into contact with the fluid issuing through the nozzle. Further downstream, vortex bending and stretching acts to generate small-scale structures as the flow becomes fully turbulent. It is evident that the jet suddenly spreads/amplifies near Near , the contours reveal a high concentration of intense mixing that persists throughout the latter half of the computational domain. These small-scale structures result in an increased chemical reaction as they serve to bring the reactants into contact. High-resolution DNS facilitates the capturing of the small-scale structures. If turbulence models were used, then the effects of the small-scale interactions on the chemical reaction (and particle formation and growth) would need to be accounted for [18, 44].