Defect-Driven Efficient Selective CO2 Hydrogenation with Mo-Based Clusters

Synthetic fuels produced from CO2 show promise in combating climate change. The reverse water gas shift (RWGS) reaction is the key to opening the CO2 molecule, and CO serves as a versatile intermediate for creating various hydrocarbons. Mo-based catalysts are of great interest for RWGS reactions featured for their stability and strong metal-oxygen interactions. Our study identified Mo defects as the intrinsic origin of the high activity of cluster Mo2C for CO2-selective hydrogenation. Specifically, we found that defected Mo2C clusters supported on nitrogen-doped graphene exhibited exceptional catalytic performance, attaining a reaction rate of 6.3 gCO/gcat/h at 400 °C with over 99% CO selectivity and good stability. Such a catalyst outperformed other Mo-based catalysts and noble metal-based catalysts in terms of facile dissociation of CO2, highly selective hydrogenation, and nonbarrier liberation of CO. Our study revealed that as a potential descriptor, the atomic magnetism linearly correlates to the liberation capacity of CO, and Mo defects facilitated product desorption by reducing the magnetization of the adsorption site. On the other hand, the defects were effective in neutralizing the negative charges of surface hydrogen, which is crucial for selective hydrogenation. Finally, we have successfully demonstrated that the combination of a carbon support and the carbonization process synergistically serves as a feasible strategy for creating rich Mo defects, and biochar can be a low-cost alternative option for large-scale applications.

Hetero-Diels-Alder Reaction between Singlet Oxygen and Anthracene Drives Integrative Cage Self-Sorting

A ZnII8L6 pseudocube containing anthracene-centered ligands, a ZnII4L'4 tetrahedron with a similar side length as the cube, and a trigonal prism ZnII6L3L'2 were formed in equilibrium from a common set of subcomponents. Hetero-Diels-Alder reaction with photogenerated singlet oxygen transformed the anthracene-containing "L" ligands into endoperoxide "LO" ones and ultimately drove the integrative self-sorting to form the trigonal prismatic cage ZnII6LO3L'2 exclusively. This ZnII6LO3L'2 structure lost dioxygen in a retro-Diels-Alder reaction after heating, which resulted in reversion to the initial ZnII8L6 + ZnII4L'4 ⇌ 2 × ZnII6L3L'2 equilibrating system. Whereas the ZnII8L6 pseudocube had a cavity too small for guest encapsulation, the ZnII6L3L'2 and ZnII6LO3L'2 trigonal prisms possessed peanut-shaped internal cavities with two isolated compartments divided by bulky anthracene panels. Guest binding was also observed to drive the equilibrating system toward exclusive formation of the ZnII6L3L'2 structure, even in the absence of reaction with singlet oxygen.

Demonstration of Hong-Ou-Mandel interference in an LNOI directional coupler

Interference between single photons is key for many quantum optics experiments and applications in quantum technologies, such as quantum communication or computation. It is advantageous to operate the systems at telecommunication wavelengths and to integrate the setups for these applications in order to improve stability, compactness and scalability. A new promising material platform for integrated quantum optics is lithium niobate on insulator (LNOI). Here, we realise Hong-Ou-Mandel (HOM) interference between telecom photons from an engineered parametric down-conversion source in an LNOI directional coupler. The coupler has been designed and fabricated in house and provides close to perfect balanced beam splitting. We obtain a raw HOM visibility of (93.5 ± 0.7) %, limited mainly by the source performance and in good agreement with off-chip measurements. This lays the foundation for more sophisticated quantum experiments in LNOI.

Fabrication of Z-scheme Bi7O9I3/g-C3N4 heterojunction modified by carbon quantum dots for synchronous photocatalytic removal of Cr (Ⅵ) and organic pollutants

Chromium(VI) (Cr(VI)), a highly toxic metal ion, generally co-exists with organic pollutants in industrial effluents. The clean and effective technology for water purification is an imperative issue but still a challenging task. A series of Bi₇O₉I₃/g-C₃N₄ (BOI/CN) composites modified by lignin-derived carbon quantum dots (CQDs) were fabricated by hydrothermal method and applied for synchronous photocatalytic removal of Cr (Ⅵ) and levofloxacin (LEV). With the modification of CQDs in BOI/CN heterojunction, the 0.5-CQD/BOI/CN photocatalyst (0.5% content of CQDs) exhibited stronger light-harvesting capacity, more efficient charge separation, and faster electron transfer. Compared to those of BOI (51.2%), CN (36.8%), and BOI/CN (74.4%), the photoreduction efficiency of Cr(VI) reached up to 100% by 0.5-CQD/BOI/CN under 60 min of light irradiation, together with 94.8% degradation efficiency of LEV. The degradation of LEV was dominantly controlled by active species (•OH and •O₂^-) identified by electron paramagnetic resonance analysis and free radical trapping experiments. The intermediates of LEV were determined by LC-MS and the possible degradation pathway was speculated in combination with density functional theory calculation, involving defluorination, decarboxylation, quinolone rings opening, and piperazine moieties oxidation reactions. This work provides an advanced strategy for the fabrication of high-efficiency CQDs-based Z-scheme photocatalysts for environmental remediation.

On the determination of Lennard-Jones parameters for polyatomic molecules

Characterizing the key length and energy scales of intermolecular interactions, Lennard-Jones parameters, i.e., collision diameter and well depth, are prerequisites for predicting transport properties and rate constants of chemical species in dilute gases. Due to anisotropy in molecular structures, Lennard-Jones parameters of many polyatomic molecules are only empirically estimated or even undetermined. This study focuses on determining the effective Lennard-Jones parameters between a polyatomic molecule and a bath gas molecule from interatomic interactions. An iterative search algorithm is developed to find orientation-dependent collision diameters and well depths on intermolecular potential energy surfaces. An orientation-averaging rule based on characteristic variables is proposed to derive the effective parameters. Cross-interaction parameters for twelve hydrocarbons with varying molecular shapes, including long-chain and planar ones, interacting with four bath gases He, Ar, N₂, and O₂ are predicted and reported. Three-dimensional parametric surfaces are constructed to quantitatively depict molecular anisotropy. Algorithmic complexity analysis and numerical experiments demonstrate that the iterative search algorithm is robust and efficient. By using the latest experimental diffusion data, it is found that the proposed orientation-averaging rule improves the prediction of cross-interaction Lennard-Jones parameters for polyatomic molecules, including for long-chain molecules that challenge the consistency of previous methods. By introducing characteristic variables, the present study shows a new route to determining effective Lennard-Jones parameters for polyatomic molecules.

Atomically Dispersed Zn-Stabilized Niδ+ Enabling Tunable Selectivity for CO2 Hydrogenation

For a heterogeneous catalytic process, the performance of catalysts could be improved by modifying the active metal with a second element. Determining the enhanced mechanism of the second element is essential to the rational design of catalysts. In this work, Zn was introduced as a second element into Ni/ZrO₂ for CO₂ hydrogenation. In contrast to Ni/ZrO₂ , the selectivity of NiZn/ZrO₂ is observed to shift from CH₄ to CO. A series of structural characterization results reveals that Zn is atomically dispersed in the NiO and ZrO₂ phases as NiZnO_x and ZnZrO_x , respectively during CO₂ hydrogenation, stabilizing a higher valence state of Ni (Ni^δ+ ) under a hydrogenation atmosphere over Ni-O-Zn site and thus promoting the generation of CO. These findings shed light on the O-mediated bimetallic effect of NiZn/ZrO₂ and bring new insight into the rational design of more efficient heterogeneous catalysts.

Triple-emission nitrogen and boron co-doped carbon quantum dots from lignin: Highly fluorescent sensing platform for detection of hexavalent chromium ions

Considering that hexavalent chromium ions (Cr⁶⁺) with high toxicity poses a huge threat to human health and the ecological environment, constructing a rapid and accurate sensing platform is of great significance in detecting the toxic substance. The novel nitrogen and boron co-doped carbon quantum dots (N, B-CQDs) from lignin are synthesized as fluorescent sensors for the detection of Cr⁶⁺. The synthetic processes involve the acid hydrolysis step followed by the hydrothermal treatment step. Lignin is firstly depolymerized by cleaving ether bonds in the acidolysis, and N, B-CQDs are consequently formed by the aromatic re-fusion of lignin nanoparticles in the hydrothermal process. The lignin-derived N, B-CQDs show triple emission of purple, blue and green fluorescence under the excitation of 300, 330, and 490 nm, respectively. The triple-emission N, B-CQDs are applied for the triple-channel detection of Cr⁶⁺, which exhibit highly sensitive and selective fluorescence quenching for Cr⁶⁺ with good linearity (R² ≤ 0.996) and very low limit of detection as 0.054, 0.049, and 0.077 μM under the excitation of 300, 330 and 490 nm, respectively. The utilization of renewable lignin as CQDs-based fluorescent sensors opens a new avenue for the rapid and accurate detection of Cr⁶⁺ through a multichannel sensing platform.

Green Synthesis of Tunable Fluorescent Carbon Quantum Dots from Lignin and Their Application in Anti-Counterfeit Printing

Lignin converted to carbon quantum dots (CQDs) attracts tremendous attention for large-scale production of carbon nanomaterials and value-added disposal of biomass wastes (such as the black liquor from pulping industry and the residue from hydrolysis of biomass). The green synthesis of lignin-derived CQDs is reported via a facile two-step method with the adjustment of acid additives containing N or S. The resulting series of CQDs exhibit bright fluorescence in gradient colors from blue to yellowish green, among which the N, S co-doped CQDs with the addition of 2,4-diaminobenzene sulfonic acid show an optimal fluorescence quantum yield (QY) of 30.5%. The red-shift photoluminescence emission behaviors of these CQDs can be attributed to the increased graphitization degree and reduced optical energy band gaps (2.47 → 2.17 eV) with regard to the incorporation of various heteroatoms. The improved fluorescence QYs are consistent with the variation trend of the increased N/C content in the CQDs. The yellowish green-emissive CQDs with bright fluorescence, strong water solubility, and excellent chemical stability perform well in anti-counterfeiting printing. The promising and sustainable approach for the synthesis of tunable fluorescent CQDs exhibits the value-added utilization of lignin for the fluorescence ink production.

A critical review on VOCs adsorption by different porous materials: Species, mechanisms and modification methods

Volatile organic compounds (VOCs) have attracted world-wide attention regarding their serious hazards on ecological environment and human health. Industrial processes such as fossil fuel combustion, petrochemicals, painting, coatings, pesticides, plastics, contributed to the large proportion of anthropogenic VOCs emission. Destructive methods (catalysis oxidation and biofiltration) and recovery methods (absorption, adsorption, condensation and membrane separation) have been developed for VOCs removal. Adsorption is established as one of the most promising strategies for VOCs abatement thanks to its characteristics of cost-effectiveness, simplicity and low energy consumption. The prominent progress in VOCs adsorption by different kinds of porous materials (such as carbon-based materials, oxygen-contained materials, organic polymers and composites is carefully summarized in this work, concerning the mechanism of adsorbate-adsorbent interactions, modification methods for the mentioned porous materials, and enhancement of VOCs adsorption capacity. This overview is to provide a comprehensive understanding of VOCs adsorption mechanisms and up-to-date progress of modification technologies for different porous materials.

Effects of Acoustic Excitation on the Combustion Instability of Hydrogen-Methane Lean Premixed Swirling Flames

Lean premixed flames are useful for low nitrogen oxide (NO _x ) emissions but more prone to induce combustion instability in gas turbines. Combustion instability of a lean premixed swirling flame (LPSF) with hydrogen-methane was investigated experimentally. The effects of hydrogen addition on combustion instability with equivalence ratios 0.75-1 were investigated with acoustic frequencies (90-240 Hz) and acoustic amplitudes (the ratio of velocity fluctuation to an average velocity of 0-0.5), respectively, which are characterized by the gain and phase of the flame describing function (FDF). The evolution of vortex and the flame morphologies were observed by the particle image velocimetry (PIV), intensified charge-coupled device (ICCD), photomultiplier tube (PMT), and Cassegrain optical systems. The global and local heat release fluctuations of the LPSF were shown by CH*/OH* chemiluminescence and temperature measurements. Results show that the FDF features maximum and minimum gain values in the acoustic frequency range of 90-240 Hz and reaches local maximum peaks at 110 and 180 Hz and local minimum peaks at 160 Hz. It can also be observed that varying velocity amplitudes (0-0.5) have greater effects on the gain and phase of FDF than changing equivalence ratios (0.75-1) for lean swirling flames. Higher velocity amplitudes more effectively intensified the compression of the flame length, which enhanced the mixing of the high-burning gas and the unburned gas, and then heat release fluctuations increased. However, it is more interesting that the effects of hydrogen addition on the combustion instability of the LPSF show a completely opposite phenomenon due to acoustic frequency under all experimental conditions. The FDFs were compared at typical frequencies of 140 and 180 Hz, and it was found that combustion instability enhanced with increasing hydrogen content at 140 Hz while weakened at 180 Hz. The flow field of PIV images shows that it is related to the location and development of vortices in the flame with varying acoustic frequencies. The intensity of OH*/CH* chemiluminescence, local temperature, and heat release rate show the same changing trend with the flame morphology for two acoustic parameters with the increasing hydrogen content in the LPSF. This directly affects the compression and curvature of the LPSF and thereby changes the mixture and temperature of the combustible gas, which influence the heat release fluctuation of the LPSF.

Efficient C-band single-photon upconversion with chip-scale Ti-indiffused pp-LiNbO3 waveguides

Frequency upconversion for single photons at telecom wavelengths is important to simultaneously meet the different wavelength requirements for long-distance communications and quantum memories in a quantum nodal network. It also enables the detection for the telecom "flying qubit" photons with silicon-based efficient single-photon detectors with low dark count (DC) rates. Here, we demonstrate the frequency upconversion of attenuated single photons, using a low-loss titanium-indiffused periodically poled lithium niobate waveguide, pumped with a readily available erbium-doped fiber amplifier in the L-band. Internal and conversion efficiencies up to 84.4% and 49.9% have been achieved, respectively. The DC rates are suppressed down to 44 kHz at 13.9% end-to-end quantum efficiency (including full conversion and detection), enabled by our long-wavelength pump configuration and narrow 3.5-GHz bandpass filtering.

Energy transfer in intermolecular collisions of polycyclic aromatic hydrocarbons with bath gases He and Ar

Classical trajectory simulations of intermolecular collisions were performed for a series of polycyclic aromatic hydrocarbons interacting with the bath gases helium and argon for bath gas temperature from 300 to 2500 K. The phase-space average energy transferred per deactivating collision, ⟨∆E_down⟩, was obtained. The Buckingham pairwise intermolecular potentials were validated against high-level quantum chemistry calculations and used in the simulations. The reactive force-field was used to describe intramolecular potentials. The dependence of ⟨∆E_down⟩ on initial vibrational energy is discussed. A canonical sampling method was compared with a microcanonical sampling method for selecting initial vibrational energy at high bath gas temperatures. Uncertainties introduced by the initial angular momentum distribution were identified. The dependence of the collisional energy transfer parameters on the type of bath gas and the molecular structure of polycyclic aromatic hydrocarbons was examined.

Flame evolution in shock-accelerated flow under different reactive gas mixture gradients

The interaction between a planar shock wave and a spherical flame is studied numerically for an ethylene-oxygen-nitrogen gas mixture. Influences of different initial reactive gas mixture gradients on the shock-flame interaction are investigated by using high-resolution computational simulations. The results show that the different reactive gas mixture gradients can greatly affect the flame evolution in shock accelerated flow. A detonation only emerges in the homogenous reactive gas mixture case, but a distinct shock bifurcation can be found in the inhomogeneous cases where the leftward reflected shock wave propagates in a reverse flow with a high transverse velocity gradient in the inhomogeneous cases. Also, the flame volume and heat release rate increase when the distribution of the reactive gas mixture is uniform or with a positive gradient in this paper, but decrease when the distribution of the reactive gas mixture is with a negative gradient, however, the ratio of unburned to burned regions in the flame zone shows just the opposite trends. Furthermore, the factors affecting the vorticity generation are also analyzed. It is found that the compression term has a relatively stronger influence on the vorticity generation in all the three cases except the period before the reflected shock wave impinges on the distorted flame in the homogeneous case, wherein the baroclinic effect dominates the vorticity generation in the flame zone.

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.

Large-scale molecular dynamics simulation of coupled dynamics of flow and glycocalyx: towards understanding atomic events on an endothelial cell surface

The glycocalyx has a prominent role in orchestrating multiple biological processes occurring at the plasma membrane. In this paper, an all-atom flow/glycocalyx system is constructed with the bulk flow velocity in the physiologically relevant ranges for the first time. The system is simulated by molecular dynamics using 5.8 million atoms. Flow dynamics and statistics in the presence of the glycocalyx are presented and discussed. Complex dynamic behaviours of the glycocalyx, particularly the sugar chains, are observed in response to blood flow. In turn, the motion of the glycocalyx, including swing and swirling, disturbs the flow by altering the velocity profiles and modifying the vorticity distributions. As a result, the initially one-dimensional forcing is spread to all directions in the region near the endothelial cell surface. Furthermore, the coupled dynamics exist not only between the flow and the glycocalyx but also within the glycocalyx molecular constituents. Shear stress distributions between one-dimer and three-dimer cases are also conducted. Finally, potential force transmission pathways are discussed based on the dynamics of the glycocalyx constituents, which provides new insight into the mechanism of mechanotransduction of the glycocalyx. These findings have relevance in the pathologies of glycocalyx-related diseases, for example in renal or cardiovascular conditions.

State-of-the-art catalytic hydrogenolysis of lignin for the production of aromatic chemicals

Catalytic hydrogenolysis of lignin is overviewed, concerning the cleavage of typical inter-unit linkages and the production of aromatic chemicals.

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.

Discrete Boltzmann modeling of Rayleigh-Taylor instability in two-component compressible flows

A discrete Boltzmann model (DBM) is proposed to probe the Rayleigh-Taylor instability (RTI) in two-component compressible flows. Each species has a flexible specific-heat ratio and is described by one discrete Boltzmann equation (DBE). Independent discrete velocities are adopted for the two DBEs. The collision and force terms in the DBE account for the molecular collision and external force, respectively. Two types of force terms are exploited. In addition to recovering the modified Navier-Stokes equations in the hydrodynamic limit, the DBM has the capability of capturing detailed nonequilibrium effects. Furthermore, we use the DBM to investigate the dynamic process of the RTI. The invariants of tensors for nonequilibrium effects are presented and studied. For low Reynolds numbers, both global nonequilibrium manifestations and the growth rate of the entropy of mixing show three stages (i.e., the reducing, increasing, and then decreasing trends) in the evolution of the RTI. On the other hand, the early reducing tendency is suppressed and even eliminated for high Reynolds numbers. Relevant physical mechanisms are analyzed and discussed.

A multi-component discrete Boltzmann model for nonequilibrium reactive flows

We propose a multi-component discrete Boltzmann model (DBM) for premixed, nonpremixed, or partially premixed nonequilibrium reactive flows. This model is suitable for both subsonic and supersonic flows with or without chemical reaction and/or external force. A two-dimensional sixteen-velocity model is constructed for the DBM. In the hydrodynamic limit, the DBM recovers the modified Navier-Stokes equations for reacting species in a force field. Compared to standard lattice Boltzmann models, the DBM presents not only more accurate hydrodynamic quantities, but also detailed nonequilibrium effects that are essential yet long-neglected by traditional fluid dynamics. Apart from nonequilibrium terms (viscous stress and heat flux) in conventional models, specific hydrodynamic and thermodynamic nonequilibrium quantities (high order kinetic moments and their departure from equilibrium) are dynamically obtained from the DBM in a straightforward way. Due to its generality, the developed methodology is applicable to a wide range of phenomena across many energy technologies, emissions reduction, environmental protection, mining accident prevention, chemical and process industry.

Consistent forcing scheme in the cascaded lattice Boltzmann method

In this paper, we give an alternative derivation for the cascaded lattice Boltzmann method (CLBM) within a general multiple-relaxation-time (MRT) framework by introducing a shift matrix. When the shift matrix is a unit matrix, the CLBM degrades into an MRT LBM. Based on this, a consistent forcing scheme is developed for the CLBM. The consistency of the nonslip rule, the second-order convergence rate in space, and the property of isotropy for the consistent forcing scheme is demonstrated through numerical simulations of several canonical problems. Several existing forcing schemes previously used in the CLBM are also examined. The study clarifies the relation between MRT LBM and CLBM under a general framework.

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.

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.

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.

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.

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.

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.

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.