Chapman-Enskog analysis for finite-volume formulation of lattice Boltzmann equation

The classical Chapman-Enskog expansion is performed for the recently proposed finite-volume formulation of lattice Boltzmann equation (LBE) method D.V. Patil, K.N. Lakshmisha, Finite volume TVD formulation of lattice Boltzmann simulation on unstructured mesh, J. Comput. Phys. 228 (2009) 5262-5279]. First, a modified partial differential equation is derived from a numerical approximation of the discrete Boltzmann equation. Then, the multi-scale, small parameter expansion is followed to recover the continuity and the Navier-Stokes (NS) equations with additional error terms. The expression for apparent value of the kinematic viscosity is derived for finite-volume formulation under certain assumptions. The attenuation of a shear wave, Taylor-Green vortex flow and driven channel flow are studied to analyze the apparent viscosity relation.

Heat and mass transfer and chemical transformation in a cerium nitrate droplet

This paper deals with the thermo-physical changes that a droplet undergoes when it is radiatively heated in a levitated environment. The heat and mass transport model has been developed along with chemical kinetics within a cerium nitrate droplet. The chemical transformation of cerium nitrate to ceria during the process is predicted using Kramers' reaction mechanism which justifies the formation of ceria at a very low temperature as observed in experiments. The rate equation modeled by Kramers is modified suitably to be applicable within the framework of a droplet, and predicts experimental results well in both bulk form of cerium nitrate and in aqueous cerium nitrate droplet. The dependence of dissociation reaction rate on droplet size is determined and the transient mass concentration of unreacted cerium nitrate is reported. The model is validated with experiments both for liquid phase vaporization and chemical reaction. Vaporization and chemical conversion are simulated for different ambient conditions. The competitive effects of sensible heating rate and the rate of vaporization with diffusion of cerium nitrate is seen to play a key role in determining the mass fraction of ceria formed within the droplet. Spatially resolved modeling of the droplet enables the understanding of the conversion of chemical species in more detail.

EFFECT OF CORE-SHELL STRUCTURE OF HYDROGEL BEADS ON THE THRESHOLD CONCENTRATION OF WATER FOR SWELLING AND ITS PH SENSITIVITY

In this paper we investigate the effect of core-shell structure of Sodium Alginate based hydrogel beads and their size on certain activation threshold concentration of water for applications in swelling and pH sensing. This type of hydrogel experiences diffusive pressure due to transport of certain free charges across its interface with a solvent or electrolyte. This process is essentially a dynamic equilibrium of the electric force field, stress in the polymeric network with cage like structure and molecular diffusion including phase transformation due to pressure imbalance between the hydrogel and its surroundings. The effect of pH of the solvant on the swelling rate of these beads has been studied experimentally. A mathematical model of the swelling process has been developed by considering Nernst-Planck equation representing the migration of mobile ions and Er ions, Poisson equation representing the equilibrium of the electric field and mechanical field equation representing swelling of the gel. An attempt has been made to predict the experimentally observed phenomena using these numerical simulations. It is observed experimentally that certain minimum concentration called activation threshold concentration of the water molecules must be present in the hydrogel in order to activate the swelling process. For the required activation threshold concentration of water in the beads, the pH induced change in the rate of swelling is also investigated. This effect is analyzed for various different core-shell structures of the beads.

The diffusion and relaxation of Gaussian chains in narrow rectangular slits

The confinement of a polymer to volumes whose characteristic linear dimensions are comparable to or smaller than its bulk radius of gyration R-G,R-bulk can produce significant changes in its static and dynamic properties, with important implications for the understanding of single-molecule processes in biology and chemistry. In this paper, we present calculations of the effects of a narrow rectangular slit of thickness d on the scaling behavior of the diffusivity D and relaxation time tau(r) of a Gaussian chain of polymerization index N and persistence length l(0). The calculations are based on the Rouse-Zimm model of chain dynamics, with the pre-averaged hydrodynamic interaction being obtained from the solutions to Stokes equations for an incompressible fluid in a parallel plate geometry in the limit of small d. They go beyond de Gennes' purely phenomenological analysis of the problem based on blobs, which has so far been the only analytical route to the determination of chain scaling behavior for this particular geometry. The present model predicts that D similar to dN(-1) ln(N/d(2)) and tau(r) similar to N(2)d(-1) ln(N/d(2))(-1) in the regime of moderate confinement, where l(0) << d < R-G,R-bulk. The corresponding results for the blob model have exactly the same power law behavior, but contain no logarithmic corrections; the difference suggests that segments within a blob may actually be partially draining and not non-draining as generally assumed.

Study on important diffusion parameters of binary Ni3Sn4 phase

Important diffusion parameters, such as-parabolic growth constant, integrated diffusivity, ratio of intrinsic diffusivities of species Ni and Sn, Kirkendall marker velocity and the activation energy for diffusion kinetics of binary Ni3Sn4 phase have been investigated with the help of incremental diffusion couple technique (Sn/Ni0.57Sn0.43) in the temperature range 200-150 degrees C. Low activation energy extracted from Arrhenius plot indicates grain boundary controlled diffusion process. The species Sn is three times faster than Ni at 200 degrees C. Further, the activation energy of Sn tracer diffusivity is greater than that of Ni.

Li diffusion through doped and defected graphene

We investigate the effect of nitrogen and boron doping on Li diffusion through defected graphene using first principles based density functional theory. While a high energy barrier rules out the possibility of Li-diffusion through the pristine graphene, the barrier reduces with the incorporation of defects. Among the most common defects in pristine graphene, Li diffusion through the divacancy encounters the lowest energy barrier of 1.34 eV. The effect of nitrogen and boron doping on the Li diffusion through doped defected-graphene sheets has been studied. N-doping in graphene with a monovacancy reduces the energy barrier significantly. The barrier reduces with the increasing number of N atoms. On the other hand, for N doped graphene with a divacancy, Li binds in the plane of the sheet, with an enhanced binding energy. The B doping in graphene with a monovacancy leads to the enhancement of the barrier. However, in the case of B-doped graphene with a divacancy, the barrier reduces to 1.54 eV, which could lead to good kinetics. The barriers do not change significantly with B concentration. Therefore, divacancy, B and N doped defected graphene has emerged as a better alternative to pristine graphene as an anode material for Li ion battery.

Combining indentation and diffusion couple techniques for combinatorial discovery of high temperature shape memory alloys

We demonstrate the possibility of accelerated identification of potential compositions for high-temperature shape memory alloys (SMAs) through a combinatorial material synthesis and analysis approach, wherein we employ the combination of diffusion couple and indentation techniques. The former was utilized to generate smooth and compositionally graded inter-diffusion zones (IDZs) in the Ni-Ti-Pd ternary alloy system of varying IDZ thickness, depending on the annealing time at high temperature. The IDZs thus produced were then impressed with an indenter with a spherical tip so as to inscribe a predetermined indentation strain. Subsequent annealing of the indented samples at various elevated temperatures, T-a, ranging between 150 and 550 degrees C allows for partial to full relaxation of the strain imposed due to the shape memory effect. If T-a is above the austenite finish temperature, A(f), the relaxation will be complete. By measuring the depth recovery, which serves as a proxy for the shape recovery characteristic of the SMA, a three-dimensional map in the recovery temperature composition space is constructed. A comparison of the published Af data for different compositions with the Ta data shows good agreement when the depth recovery is between 70% and 80%, indicating that the methodology proposed in this paper can be utilized for the identification of promising compositions. Advantages and further possibilities of this methodology are discussed.

Scalar transport in diffusion flame wrapped up by an air and fuel side vortex

The present work involves a computational study of soot (chosen as a scalar which is a primary pollutant source) formation and transport in a laminar acetylene diffusion flame perturbed by a convecting line vortex. The topology of soot contours resulting from flame vortex interactions has been investigated. More soot was produced when vortex was introduced from the air side in comparison to the fuel side. Also, the soot topography was spatially more diffuse in the case of air side vortex. The computational model was found to be in good agreement with the experimental work previously reported in the literature. The computational simulation enabled a study of various parameters like temperature, equivalence ratio and temperature gradient affecting the soot production and transport. Temperatures were found to be higher in the case of air side vortex in contrast to the fuel side one. In case of fuel side vortex, abundance of fuel in the vortex core resulted in fuel-rich combustion zone in the core and a more discrete soot topography. Besides, the overall soot production was observed to be low in the fuel side vortex. However, for the air side vortex, air abundance in the core resulted in higher temperatures and greater soot production. Probability density functions (PDFs) have been introduced to investigate the spatiotemporal variation of soot yield and transport and their dependence on temperature and acetylene concentration from statistical view point. In addition, the effect of flame curvature on soot production is also studied. The regions convex to fuel stream side witnessed thicker soot layer. All numerical simulations have been carried out on Fluent 6.3.26. (C) 2013 Elsevier Ltd. All rights reserved.

The diffusion and relaxation of Gaussian chains in narrow rectangular slits

The confinement of a polymer to volumes whose characteristic linear dimensions are comparable to or smaller than its bulk radius of gyration R-G,R-bulk can produce significant changes in its static and dynamic properties, with important implications for the understanding of single-molecule processes in biology and chemistry. In this paper, we present calculations of the effects of a narrow rectangular slit of thickness d on the scaling behavior of the diffusivity D and relaxation time tau(r) of a Gaussian chain of polymerization index N and persistence length l(0). The calculations are based on the Rouse-Zimm model of chain dynamics, with the pre-averaged hydrodynamic interaction being obtained from the solutions to Stokes equations for an incompressible fluid in a parallel plate geometry in the limit of small d. They go beyond de Gennes' purely phenomenological analysis of the problem based on blobs, which has so far been the only analytical route to the determination of chain scaling behavior for this particular geometry. The present model predicts that D similar to dN(-1) ln(N/d(2)) and tau(r) similar to N(2)d(-1) ln(N/d(2))(-1) in the regime of moderate confinement, where l(0) << d < R-G,R-bulk. The corresponding results for the blob model have exactly the same power law behavior, but contain no logarithmic corrections; the difference suggests that segments within a blob may actually be partially draining and not non-draining as generally assumed. (C) 2013 AIP Publishing LLC.

Solution of Time Dependent Joule Heat Equation for a Graphene Sheet Under Thomson Effect

We address a physics-based solution of joule heating phenomenon in a single-layer graphene (SLG) sheet under the presence of Thomson effect. We demonstrate that the temperature in an isotopically pure (containing only C-12) SLG sheet attains its saturation level quicker than when doped with its isotopes (C-13). From the solution of the joule heating equation, we find that the thermal time constant of the SLG sheet is in the order of tenths of a nanosecond for SLG dimensions of a few micrometers. These results have been formulated using the electron interactions with the inplane and flexural phonons to demonstrate a field-dependent Landauer transmission coefficient. We further develop an analytical model of the SLG specific heat using the quadratic (out of plane) phonon band structure over the room temperature. Additionally, we show that a cooling effect in the SLG sheet can be substantially enhanced with the addition of C-13. The methodologies as discussed in this paper can be put forward to analyze the graphene heat spreader theory.

Ground motion prediction equation considering combined dataset of recorded and simulated ground motions

Himalayan region is one of the most active seismic regions in the world and many researchers have highlighted the possibility of great seismic event in the near future due to seismic gap. Seismic hazard analysis and microzonation of highly populated places in the region are mandatory in a regional scale. Region specific Ground Motion Predictive Equation (GMPE) is an important input in the seismic hazard analysis for macro- and micro-zonation studies. Few GMPEs developed in India are based on the recorded data and are applicable for a particular range of magnitudes and distances. This paper focuses on the development of a new GMPE for the Himalayan region considering both the recorded and simulated earthquakes of moment magnitude 5.3-8.7. The Finite Fault simulation model has been used for the ground motion simulation considering region specific seismotectonic parameters from the past earthquakes and source models. Simulated acceleration time histories and response spectra are compared with available records. In the absence of a large number of recorded data, simulations have been performed at unavailable locations by adopting Apparent Stations concept. Earthquakes recorded up to 2007 have been used for the development of new GMPE and earthquakes records after 2007 are used to validate new GMPE. Proposed GMPE matched very well with recorded data and also with other highly ranked GMPEs developed elsewhere and applicable for the region. Comparison of response spectra also have shown good agreement with recorded earthquake data. Quantitative analysis of residuals for the proposed GMPE and region specific GMPEs to predict Nepal-India 2011 earthquake of Mw of 5.7 records values shows that the proposed GMPE predicts Peak ground acceleration and spectral acceleration for entire distance and period range with lower percent residual when compared to exiting region specific GMPEs. Crown Copyright (C) 2013 Published by Elsevier Ltd. All rights reserved.

Scaling of turbulent flame speed for expanding flames with Markstein diffusion considerations

In this paper we clarify the role of Markstein diffusivity, which is the product of the planar laminar flame speed and the Markstein length, on the turbulent flame speed and its scaling, based on experimental measurements on constant-pressure expanding turbulent flames. Turbulent flame propagation data are presented for premixed flames of mixtures of hydrogen, methane, ethylene, n-butane, and dimethyl ether with air, in near-isotropic turbulence in a dual-chamber, fan-stirred vessel. For each individual fuel-air mixture presented in this work and the recently published iso-octane data from Leeds, normalized turbulent flame speed data of individual fuel-air mixtures approximately follow a Re-T,f(0.5) scaling, for which the average radius is the length scale and thermal diffusivity is the transport property of the turbulence Reynolds number. At a given Re-T,Re-f, it is experimentally observed that the normalized turbulent flame speed decreases with increasing Markstein number, which could be explained by considering Markstein diffusivity as the leading dissipation mechanism for the large wave number flame surface fluctuations. Consequently, by replacing thermal diffusivity with the Markstein diffusivity in the turbulence Reynolds number definition above, it is found that normalized turbulent flame speeds could be scaled by Re-T,M(0.5) irrespective of the fuel, equivalence ratio, pressure, and turbulence intensity for positive Markstein number flames.

Turbulence in the two-dimensional Fourier-truncated Gross-Pitaevskii equation

We undertake a systematic, direct numerical simulation of the twodimensional, Fourier-truncated, Gross-Pitaevskii equation to study the turbulent evolutions of its solutions for a variety of initial conditions and a wide range of parameters. We find that the time evolution of this system can be classified into four regimes with qualitatively different statistical properties. Firstly, there are transients that depend on the initial conditions. In the second regime, powerlaw scaling regions, in the energy and the occupation-number spectra, appear and start to develop; the exponents of these power laws and the extents of the scaling regions change with time and depend on the initial condition. In the third regime, the spectra drop rapidly for modes with wave numbers k > kc and partial thermalization takes place for modes with k < kc; the self-truncation wave number kc(t) depends on the initial conditions and it grows either as a power of t or as log t. Finally, in the fourth regime, complete thermalization is achieved and, if we account for finite-size effects carefully, correlation functions and spectra are consistent with their nontrivial Berezinskii-Kosterlitz-Thouless forms. Our work is a natural generalization of recent studies of thermalization in the Euler and other hydrodynamical equations; it combines ideas from fluid dynamics and turbulence, on the one hand, and equilibrium and nonequilibrium statistical mechanics on the other.

Mechanistic Insights into a Non-Classical Diffusion Pathway for the Formation of Hollow Intermetallics: A Route to Multicomponent Hollow Structures

Hollow nanostructures are used for various applications including catalysis, sensing, and drug delivery. Methods based on the Kirkendall effect have been the most successful for obtaining hollow nanostructures of various multicomponent systems. The classical Kirkendall effect relies on the presence of a faster diffusing species in the core; the resultant imbalance in flux results in the formation of hollow structures. Here, an alternate non-Kirkendall mechanism that is operative for the formation of hollow single crystalline particles of intermetallic PtBi is demonstrated. The synthesis method involves sequential reduction of Pt and Bi salts in ethylene glycol under microwave irradiation. Detailed analysis of the reaction at various stages indicates that the formation of the intermetallic PtBi hollow nanoparticles occurs in steps. The mechanistic details are elucidated using control experiments. The use of microwave results in a very rapid synthesis of intermetallics PtBi that exhibits excellent electrocatalytic activity for formic acid oxidation reaction. The method presented can be extended to various multicomponent systems and is independent of the intrinsic diffusivities of the species involved.