An Investigation towards Coupling Molecular Dynamics with Computational Fluid Dynamics for Modelling Polymer Pyrolysis

Building polymers implemented into building panels and exterior façades have been determined as the major contributor to severe fire incidents, including the 2017 Grenfell Tower fire incident. To gain a deeper understanding of the pyrolysis process of these polymer composites, this work proposes a multi-scale modelling framework comprising of applying the kinetics parameters and detailed pyrolysis gas volatiles (parent combustion fuel and key precursor species) extracted from Molecular Dynamics models to a macro-scale Computational Fluid Dynamics fire model. The modelling framework was tested for pure and flame-retardant polyethylene systems. Based on the modelling results, the chemical distribution of the fully decomposed chemical compounds was realised for the selected polymers. Subsequently, the identified gas volatiles from solid to gas phases were applied as the parent fuel in the detailed chemical kinetics combustion model for enhanced predictions of toxic gas, charring, and smoke particulate predictions. The results demonstrate the potential application of the developed model in the simulation of different polymer materials without substantial prior knowledge of the thermal degradation properties from costly experiments.     

Computational Analysis of Mixing and Scale-Up in Emulsion Polymerization Reactors

Summary: A hybrid multi-zonal/computational fluid dynamics (CFD) framework is currently being developed to aid in the scale-up of high solid content latex production and processing. Poly3D, a commercial laminar CFD code tailored to modelling the mixing of non-Newtonian fluids, has been coupled to a population balance model via a customized interface. CFD is used to generate flow fields inside a series of reactors; this information is then transferred to a multi-zone population balance model to assess the impact of non-homogenous mixing on the evolution of the latex particle size distribution (PSD) when concentrated latex suspension is altered via the addition of a coagulant. The rheological properties of high solid content latexes are sensitive to changes in the PSD, so the flow field is periodically updated if significant changes in the rheological properties of the latex are detected in any of the zones. The details of the models comprising the framework are presented and the utility of the framework is demonstrated.     

On the Importance of Mixing for the Production of Nanoparticles

The precipitation process of nanoparticles of barium‐sulphate in a confined impinging jet reactor (CIJR) is here modelled by means of a precipitation model based on computational fluid dynamics (CFD). The micromixing was described with the interchange by exchange with the mean model coupled with the direct quadrature method of moments algorithm (DQMOM‐IEM), and the population balance equation (PBE), which describes the nucleation, growth and aggregation steps, was solved with the approach of the quadrature method of moments (QMOM). In this work preliminary results for the validation of the model are presented.     

Hemodynamics and Wall Shear Stress of Blood Vessels in Aortic Coarctation with Computational Fluid Dynamics Simulation

The purpose of this study was to identify the characteristics of blood flow in aortic coarctation based on stenotic shape structure, stenosis rate, and the distribution of the wall load delivered into the blood vessels and to predict the impact on aneurysm formation and rupture of blood vessels by using a computational fluid dynamics modeling method. It was applied on the blood flow in abdominal aortic blood vessels in which stenosis occurred by using the commercial finite element software ADINA on fluid-solid interactions. The results of modeling, with an increasing stenosis rate and Reynolds number, showed the pressure drop was increased and the velocity was greatly changed. When the stenosis rate was the same, the pressure drop and the velocity change were larger in the stenosis with a symmetric structure than in the stenosis with an asymmetric one. Maximal changes in wall shear stress were observed in the area before stenosis and minimal changes were shown in stenosis areas. The minimal shear stress occurred at different locations depending on the stenosis shape models. With an increasing stenosis rate and Reynolds number, the maximal wall shear stress was increased and the minimal wall shear stress was decreased. Through such studies, it is thought that the characteristics of blood flow in the abdominal aorta where a stenosis is formed will be helpful in understanding the mechanism of growth of atherosclerosis and the occurrence and rupture of the abdominal aortic flow.     

制备超细微粒的超临界流体沉淀技术新进展

超临界流体沉淀(SFP)技术以其特有的优点成为具有广阔应用前景的超细微粒制备方法.本文着重综述了气体饱和溶液沉析(PGSS)、超临界辅助雾化(SAA)和膨胀液体有机溶液减压(DELOS)等最新SFP技术的原理、特点及相关的应用,指出了目前SFP技术存在的问题及今后的研究方向.     

Characterization of the Engineering Flow in a Miniaturised Bioreactor by Computational Fluid Dynamics

Three approaches based on computational fluid dynamics(CFD) techniques have been assessed for their ability to describe the engineering flow environment in a miniaturized mechanically agitated bioreactor. The three approaches tested were the source-sink(SS), the multiple reference frames (MRF) and the sliding grids(SG). In all the cases, the predictions of the velocity components agree with reported experimental data. However, the analysis of the results of the turbulent intensities predicted by the three approaches indicates the MRF and the SG techniques under predicted turbulent intensities are comparable to both experimental measurements and the SS method. The predicted power number and pumping number based on the SS ap-proach are closer to typical reported experimental values compared to those obtained from the MRF and SG methods.     

Flow field in a shrinking-bed reactor for pretreatment of cellulosic biomass

A shrinking-bed reactor was designed by the National Renewable Energy Laboratory to maintain a constant bulk packing density of cellulosic biomass. The high solid-to-liquid ratio in the pretreatment process allows a high sugar yield and avoids the need to flush large volumes of solution through the reactor. The shrinking-bed reactor is a promising pretreatment reactor with the potential for scale-up for commercial applications. To scale up the shrinking-bed reactor, it is necessary to understand the flow pattern in the reactor. In this study, flow field is simulated with computational fluid dynamics using a porous medium model. Different discrete “snapshots” and multiple steady states are utilized. The bulk flow pattern, velocity distribution, and pressure drop are determined from the simulation and can be used to guide reactor design and scale-up.     

Free Energy Change of Micelle Formation for Sodium Dodecyl Sulfate from a Dispersed State in Solution to Complete Micelles along Its Aggregation Pathways Evaluated by Chemical Species Model Combined with Molecular Dynamics Calculations

Surfactant molecules, when dispersed in solution, have been shown to spontaneously form aggregates. Our previous studies on molecular dynamics(MD) calculations have shown that ionic sodium dodecyl sulfate molecules quickly aggregated even when the aggregation number is small. The aggregation rate, however, decreased for larger aggregation numbers. In addition, studies have shown that micelle formation was not completed even after a 100 ns-long MD run(Chem. Phys. Lett. 2016, 646, 36). Herein, we analyze the free energy change of micelle formation based on chemical species model combined with molecular dynamics calculations. First, the free energy landscape of the aggregation, ?G_(i+j)~+, where two aggregates with sizes i and j associate to form the(i + j)-mer, was investigated using the free energy of micelle formation of the i-mer, G_i~+, which was obtained through MD calculations. The calculated ?G_(i+j)~+ was negative for all the aggregations where the sum of DS ions in the two aggregates was 60 or less. From the viewpoint of chemical equilibrium, aggregation to the stable micelle is desired. Further, the free energy profile along possible aggregation pathways was investigated, starting from small aggregates and ending with the complete thermodynamically stable micelles in solution. The free energy profiles, G(l, k), of the aggregates at l-th aggregation path and k-th state were evaluated by the formation free energy ∑_in_i( l,k)G_i~+ and the free energy of mixing ∑_in_i( l,k)k_BTln( n_i( l,k)/n( l,k)), where ni(l, k) is the number of i-mer in the system at the l-th i aggregation path and k-th state, with n(l,k)= ∑_n_i( l,k). All the aggregation pathways were obtained from the initial i state of 12 pentamers to the stable micelle with i = 60. All the calculated G(l, k) values monotonically decreased with increasing k. This indicates that there are no free energy barriers along the pathways. Hence, the slowdown is not due to the thermodynamic stability of the aggregates, but rather the kinetics that inhibit the association of the fragments. The time required for a collision between aggregates, one of the kinetic factors, was evaluated using the fast passage time, t_(FPT). The calculated t_(FPT) was about 20 ns for the aggregates with N = 31. Therefore, if aggregation is a diffusion-controlled process, it should be completed within the 100 ns-simulation. However, aggregation does not occur due to the free energy barrier between the aggregates, that is, the repulsive force acting on them. This may be caused by electrostatic repulsions produced by the overlap of the electric double layers, which are formed by the negative charge of the hydrophilic groups and counter sodium ions on the surface of the aggregates.     

Model reactions for abiotic mercury(II) methylation: Kinetics of methylation of mercury(II) by mono-, di-, and tri-methyltin in seawater

The usual presence of mercury(II) with monodi-, and tri-methyltin in water, sediments, and plants in estuarine environments suggests possible abiotic formation of methylmercury via methyl transfer from methyltin compounds. Kinetics studies of reactions between mercury(II) and methyltin compounds under pseudo-first-order conditions in seawater show that relative rate of methylmercury formation under the same conditions are: monomethyltin <trimethyltin> dimethyltin. This order is explainable mainly by the speciation and charge of methyltin compounds in seawater and by the existence of mercury(II) as a tetrachloro anion. A factorial experiment with the variables pH and salinity (seawater diluted with deionized water) showed that pH, but not salinity, is significant at the 95% confidence level; and that reaction rates increase as pH increases. These results suggest the possibility of abiotic methylation of mercury(II) in seawater. Additional experiments in seawater demonstrated an absence of methylation of mercury(II) (14 days) and mercury(0) (35 days) by methyl iodide.