Evaluation of the effect of nickel clusters on the formation of incipient soot particles from polycyclic aromatic hydrocarbons via ReaxFF molecular dynamics simulations.

In the present study, the ReaxFF reactive molecular dynamics simulation method was applied to investigate the effect of a small nickel cluster (Ni13) on the formation of nascent soot from polycyclic aromatic hydrocarbon (PAH) precursors. A series of NVT simulations was performed for systems of a Ni13 cluster and various PAH monomers, namely, naphthalene, anthracene, pyrene, coronene, ovalene, and circumcoronene, at temperatures from 400 to 2500 K. At low temperatures, the PAHs form soot particles via binding and stacking around nickel clusters. Larger soot particles are formed due to the early initiation of clustering provided by nickel compared to those observed in homogenous PAH systems. At 1200-1600 K, the PAH monomers show a chemisorption tendency onto the nickel surface, which results in incipient soot particles. Chemical nucleation was observed at 2000 K where nickel-assisted dehydrogenation and chemisorption of PAH led to the growth of stable soot particles, which did not occur in the absence of Ni-clusters. At a high temperature (2500 K), nickel significantly accelerates the ring-opening and graphitization of PAH molecules and increases the size of the fullerene-type soot as compared to that of homogenous systems.

Dynamics and kinetics of reversible homo-molecular dimerization of polycyclic aromatic hydrocarbons.

Physical dimerization of polycyclic aromatic hydrocarbons (PAHs) has been investigated via molecular dynamics (MD) simulation with the ReaxFF reactive force field that is developed to bridge the gap between the quantum mechanism and classical MD. Dynamics and kinetics of homo-molecular PAH collision under different temperatures, impact parameters, and orientations are studied at an atomic level, which is of great value to understand and model the PAH dimerization. In the collision process, the enhancement factors of homo-molecular dimerizations are quantified and found to be larger at lower temperatures or with smaller PAH instead of size independent. Within the capture radius, the lifetime of the formed PAH dimer decreases as the impact parameter increases. Temperature and PAH characteristic dependent forward and reverse rate constants of homo-molecular PAH dimerization are derived from MD simulations, on the basis of which a reversible model is developed. This model can predict the tendency of PAH dimerization as validated by pyrene dimerization experiments [H. Sabbah et al., J. Phys. Chem. Lett. 1(19), 2962 (2010)]. Results from this study indicate that the physical dimerization cannot be an important source under the typical flame temperatures and PAH concentrations, which implies a more significant role played by the chemical route.

A practical procedure for the cryopreservation of marrow cells intended for autotransplantation.

A simple and practical method of unfractionated bone marrow processing and cryopreservation was studied. The date shows that RBCs can be rapidly sedimented by methylcellulose or sodium carboxymethyl starch within 15-45 min. The cells can be cryopreserved in a mixture consisting of 5% DMSO, 6% HES and 4% Albumin prepared in a Lactated Ringer solution which is widely used, and can be simply immersed into a -80 degrees C freezer and stored in liquid nitrogen until infusion. Recovery percentages of nucleated cell, cell viability and CFU-G were similar to those cryopreserved with the conventional method. Clinical toxicity was mild in 12 infused patients. Of them 8 patients had received high dose chemotherapy +/- TBI and their peripheral WBC recovery was rapid. The recovery of WBC or Platelet (PLT) in the study group was similar to that of the control group whose marrow cells were cryopreserved in 10% DMSO. Therefore, cells cryopreserved with this method can also accelerate the hematopoietic recovery in myeloablatively treated patients.

Effect of the forcing term in the pseudopotential lattice Boltzmann modeling of thermal flows.

The pseudopotential lattice Boltzmann (LB) model is a popular model in the LB community for simulating multiphase flows. Recently, several thermal LB models, which are based on the pseudopotential LB model and constructed within the framework of the double-distribution-function LB method, were proposed to simulate thermal multiphase flows [G. Házi and A. Márkus, Phys. Rev. E 77, 026305 (2008); L. Biferale, P. Perlekar, M. Sbragaglia, and F. Toschi, Phys. Rev. Lett. 108, 104502 (2012); S. Gong and P. Cheng, Int. J. Heat Mass Transfer 55, 4923 (2012); M. R. Kamali et al., Phys. Rev. E 88, 033302 (2013)]. The objective of the present paper is to show that the effect of the forcing term on the temperature equation must be eliminated in the pseudopotential LB modeling of thermal flows. First, the effect of the forcing term on the temperature equation is shown via the Chapman-Enskog analysis. For comparison, alternative treatments that are free from the forcing-term effect are provided. Subsequently, numerical investigations are performed for two benchmark tests. The numerical results clearly show that the existence of the forcing-term effect will lead to significant numerical errors in the pseudopotential LB modeling of thermal flows.

Contact angles in the pseudopotential lattice Boltzmann modeling of wetting.

In this paper we investigate the implementation of contact angles in the pseudopotential lattice Boltzmann modeling of wetting at a large density ratio ρ_{L}/ρ_{V}=500. The pseudopotential lattice Boltzmann model [X. Shan and H. Chen, Phys. Rev. E 49, 2941 (1994)10.1103/PhysRevE.49.2941] is a popular mesoscopic model for simulating multiphase flows and interfacial dynamics. In this model the contact angle is usually realized by a fluid-solid interaction. Two widely used fluid-solid interactions, the density-based interaction and the pseudopotential-based interaction, as well as a modified pseudopotential-based interaction formulated in the present paper are numerically investigated and compared in terms of the achievable contact angles, the maximum and the minimum densities, and the spurious currents. It is found that the pseudopotential-based interaction works well for simulating small static (liquid) contact angles θ<90^{∘}, however, it is unable to reproduce static contact angles close to 180^{∘}. Meanwhile, it is found that the proposed modified pseudopotential-based interaction performs better in light of the maximum and the minimum densities and is overall more suitable for simulating large contact angles θ>90^{∘} as compared with the two other types of fluid-solid interactions. Furthermore, the spurious currents are found to be enlarged when the fluid-solid interaction force is introduced. Increasing the kinematic viscosity ratio between the vapor and liquid phases is shown to be capable of reducing the spurious currents caused by the fluid-solid interactions.

Achieving tunable surface tension in the pseudopotential lattice Boltzmann modeling of multiphase flows.

In this paper, we aim to address an important issue about the pseudopotential lattice Boltzmann (LB) model, which has attracted much attention as a mesoscopic model for simulating interfacial dynamics of complex fluids, but suffers from the problem that the surface tension cannot be tuned independently of the density ratio. In the literature, a multirange potential was devised to adjust the surface tension [Sbragaglia et al., Phys. Rev. E 75, 026702 (2007)]. However, it was recently found that the density ratio of the system will be changed when the multirange potential is employed to adjust the surface tension. An alternative approach is therefore proposed in the present work. The basic strategy is to add a source term to the LB equation so as to tune the surface tension of the pseudopotential LB model. The proposed approach can guarantee that the adjustment of the surface tension does not affect the mechanical stability condition of the pseudopotential LB model, and thus provides a separate control of the surface tension and the density ratio. Meanwhile, it still retains the mesoscopic feature and the computational simplicity of the pseudopotential LB model. Numerical simulations are carried out for stationary droplets, capillary waves, and droplet splashing on a thin liquid film. The numerical results demonstrate that the proposed approach is capable of achieving a tunable surface tension over a very wide range and can keep the density ratio unchanged when adjusting the surface tension.

Lattice Boltzmann modeling of multiphase flows at large density ratio with an improved pseudopotential model.

Owing to its conceptual simplicity and computational efficiency, the pseudopotential multiphase lattice Boltzmann (LB) model has attracted significant attention since its emergence. In this work, we aim to extend the pseudopotential LB model to simulate multiphase flows at large density ratio and relatively high Reynolds number. First, based on our recent work [Q. Li, K. H. Luo, and X. J. Li, Phys. Rev. E 86, 016709 (2012)], an improved forcing scheme is proposed for the multiple-relaxation-time pseudopotential LB model in order to achieve thermodynamic consistency and large density ratio in the model. Next, through investigating the effects of the parameter a in the Carnahan-Starling equation of state, we find that the interface thickness is approximately proportional to 1/√a. Using a smaller a will lead to a wider interface thickness, which can reduce the spurious currents and enhance the numerical stability of the pseudopotential model at large density ratio. Furthermore, it is found that a lower liquid viscosity can be gained in the pseudopotential model by increasing the kinematic viscosity ratio between the vapor and liquid phases. The improved pseudopotential LB model is numerically validated via the simulations of stationary droplet and droplet oscillation. Using the improved model as well as the above treatments, numerical simulations of droplet splashing on a thin liquid film are conducted at a density ratio in excess of 500 with Reynolds numbers ranging from 40 to 1000. The dynamics of droplet splashing is correctly reproduced and the predicted spread radius is found to obey the power law reported in the literature.

Forcing scheme in pseudopotential lattice Boltzmann model for multiphase flows.

The pseudopotential lattice Boltzmann (LB) model is a widely used multiphase model in the LB community. In this model, an interaction force, which is usually implemented via a forcing scheme, is employed to mimic the molecular interactions that cause phase segregation. The forcing scheme is therefore expected to play an important role in the pseudoepotential LB model. In this paper, we aim to address some key issues about forcing schemes in the pseudopotential LB model. First, theoretical and numerical analyses will be made for Shan-Chen's forcing scheme [Shan and Chen, Phys. Rev. E 47, 1815 (1993)] and the exact-difference-method forcing scheme [Kupershtokh et al., Comput. Math. Appl. 58, 965 (2009)]. The nature of these two schemes and their recovered macroscopic equations will be shown. Second, through a theoretical analysis, we will reveal the physics behind the phenomenon that different forcing schemes exhibit different performances in the pseudopotential LB model. Moreover, based on the analysis, we will present an improved forcing scheme and numerically demonstrate that the improved scheme can be treated as an alternative approach to achieving thermodynamic consistency in the pseudopotential LB model.

Coupling lattice Boltzmann model for simulation of thermal flows on standard lattices.

In this paper, a coupling lattice Boltzmann (LB) model for simulating thermal flows on the standard two-dimensional nine-velocity (D2Q9) lattice is developed in the framework of the double-distribution-function (DDF) approach in which the viscous heat dissipation and compression work are considered. In the model, a density distribution function is used to simulate the flow field, while a total energy distribution function is employed to simulate the temperature field. The discrete equilibrium density and total energy distribution functions are obtained from the Hermite expansions of the corresponding continuous equilibrium distribution functions. The pressure given by the equation of state of perfect gases is recovered in the macroscopic momentum and energy equations. The coupling between the momentum and energy transports makes the model applicable for general thermal flows such as non-Boussinesq flows, while the existing DDF LB models on standard lattices are usually limited to Boussinesq flows in which the temperature variation is small. Meanwhile, the simple structure and general features of the DDF LB approach are retained. The model is tested by numerical simulations of thermal Couette flow, attenuation-driven acoustic streaming, and natural convection in a square cavity with small and large temperature differences. The numerical results are found to be in good agreement with the analytical solutions and/or other numerical results reported in the literature.

Additional interfacial force in lattice Boltzmann models for incompressible multiphase flows.

The existing lattice Boltzmann models for incompressible multiphase flows are mostly constructed with two distribution functions: one is the order parameter distribution function, which is used to track the interface between different phases, and the other is the pressure distribution function for solving the velocity field. In this paper, it is shown that in these models the recovered momentum equation is inconsistent with the target one: an additional force is included in the recovered momentum equation. The additional force has the following features. First, it is proportional to the macroscopic velocity. Second, it is zero in every single-phase region but is nonzero in the interface. Therefore it can be interpreted as an interfacial force. To investigate the effects of the additional interfacial force, numerical simulations are carried out for the problem of Rayleigh-Taylor instability, droplet splashing on a thin liquid film, and the evolution of a falling droplet under gravity. Numerical results demonstrate that, with the increase of the velocity or the Reynolds number, the additional interfacial force will gradually have an important influence on the interface and affect the numerical accuracy.

The pyrolytic degradation of wood-derived lignin from pulping process.

Lignin is a key component in the biomass with a complex polymeric structure of the phenyl-C(3) alkyl units. The kraft lignin from the wood pulping process is tested in TG-FTIR and Py-GC-MS. The samples are pyrolyzed in TGA coupled with FTIR from 30 to 900 degrees C at the heating rate of 20 and 40K/min. The evolution of phenolic compounds in the initial pyrolysis stage of lignin is determined by FTIR, while the second stage is mainly attributed to the production of the low molecular weight species. A bench-scale fast pyrolysis unit is employed to investigate the effect of temperature on the product yield and composition. It is found that the guaiacol-type and syringol-type compounds as the primary products of lignin pyrolysis are predominant in bio-oil, acting as the significant precursors for the formation of the derivatives such as the phenol-, cresol- and catechol-types. A series of free-radical chain-reactions, concerning the cracking of different side-chain structures and the methoxy groups on aromatic ring, are proposed to demonstrate the formation pathways for the typical compounds in bio-oil by closely relating lignin structure to the pyrolytic mechanisms. The methoxy group (-OCH(3)) is suggested to work as an important source for the formation of the small volatile species (CO, CO(2) and CH(4)) through the relevant free radical coupling reactions.