Structure optimization of gas-liquid combined loop reactor using a CFD-PBE coupled model

Flow characteristics, such as flow pattern, gas holdup, and bubble size distribution, in an internal loop reactor with external liquid circulation, are simulated to investigate the influence of reactor internals by using the computational fluid dynamics (CFD)-population balance equations (PBE) coupled model. Numerical results reveal that introducing a downcomer tube and a draft tube can help to improve the mass and heat transfer of the reactor through enhanced liquid circulation, increased gas holdup and reduced bubble diameter. The hydrodynamic behavior in the internal loop reactor with external liquid circulation can be managed effectively by adjusting the diameter and axial position of the draft tube.     

Numerical simulation of predicting and reducing solid particle erosion of solid-liquid two-phase flow in a choke

Chokes are one of the most important components of downhole flow-control equipment. The particle erosion mathematical model, which considers particle-particle interaction, was established and used to simulate solid particle movement as well as particle erosion characteristics of the solid-liquid two-phase flow in a choke. The corresponding erosion reduction approach by setting ribs on the inner wall of the choke was advanced. This mathematical model includes three parts： the flow field simulation of the continuous carrier fluid by an Eulerian approach, the particle interaction simulation using the discrete particle hard sphere model by a Lagrangian approach and calculation of erosion rate using semiempirical correlations. The results show that particles accumulated in a narrow region from inlet to outlet of the choke and the dominating factor affecting particle motion is the fluid drag force. As a result, the optimization of rib geometrical parameters indicates that good anti-erosion performance can be achieved by four ribs, each of them with a height （H） of 3 mm and a width （B） of 5 mm equaling the interval between ribs （L）.     

A criterion for flow mechanisms through vertical sharp-edged orifice and model for the orifice discharge coefficient

According to the experimental data of the orifice discharge coefficient for the flow through a vertical sharp-edged orifice obtained in the previous study of this work,a theoretical criterion for flow mechanisms of small orifice(viz.thick-walled orifice and nozzle) and large orifice(viz.thin-walled orifice) was proposed based on the ratio of orifice diameter to plate thickness.It can help explain the dissipation of the mechanical energy loss in the flow process for the two flow mechanisms under different operating regimes.The main parameters such as orifice diameter,plate thickness and liquid head were correlated,and a semi-empirical model for orifice coefficient and an empirical model with high precision at the stable region were developed.     

A mechanistic model of heat transfer for gas–liquid flow in vertical wellbore annuli

The most prominent aspect of multiphase flow is the variation in the physical distribution of the phases in the flow conduit known as the flow pattern. Several different flow patterns can exist under different flow conditions which have significant effects on liquid holdup, pressure gradient and heat transfer. Gas–liquid two-phase flow in an annulus can be found in a variety of practical situations. In high rate oil and gas production, it may be beneficial to flow fluids vertically through the annulus configuration between well tubing and casing. The flow patterns in annuli are different from pipe flow. There are both casing and tubing liquid films in slug flow and annular flow in the annulus. Multiphase heat transfer depends on the hydrodynamic behavior of the flow. There are very limited research results that can be found in the open literature for multiphase heat transfer in wellbore annuli. A mechanistic model of multiphase heat transfer is developed for different flow patterns of upward gas–liquid flow in vertical annuli. The required local flow parameters are predicted by use of the hydraulic model of steady-state multiphase flow in wellbore annuli recently developed by Yin et al. The modified heat-transfer model for single gas or liquid flow is verified by comparison with Manabe's experimental results. For different flow patterns, it is compared with modified unified Zhang et al. model based on representative diameters.     

CFD simulation on the gasification of asphalt water slurry in an entrained flow gasifier

Gasification technology is suggested to utilize asphalt particles, which are produced in the heavy oil deep separation process of using coupled low temperature separation of solvent and post extraction residue. In this work, the asphalt particles were first slurried with water and then gasified to produce synthesis gas. The gasification process of asphalt water slurry in an entrained flow gasifier was simulated using a three-dimensional computational fluid dynamics （CFD） model based on an Eulerian- Lagrangian method. The trajectories and residence time of asphalt particles, and the reaction rates, gas species distribution, temperature field and carbon conversion in the entrained flow gasifier were obtained. The predicted results indicated that the asphalt water slurry was a good feedstock for gasification. Moreover, the effects of particle size, oxygen equivalence ratio, and mass content of asphalt particles on the gasification performance of asphalt water slurry were investigated. These results are helpful for industrial application of asphalt water slurry gasification technology.     

Experimental study and a proposed new approach for thermodynamic modeling of wax precipitation in crude oil using a PC-SAFT model

A powerful method is necessary for thermodynamic modeling of wax phase behavior in crude oils,such as the perturbed-chain statistical associating fluid theory(PC-SAFT).In this work,a new approach based on the wax appearance temperature of crude oil was proposed to estimate PC-SAFT parameters in thermodynamic modeling of wax precipitation from crude oil.The proposed approach was verified using experimental data obtained in this work and also with those reported in the literature.In order to compare the performance of the PC-SAFT model with previous models,the wax precipitation experimental data were correlated using previous models such as the solid solution model and multi-solid phase model.The results showed that the PC-SAFT model can correlate more accurately the wax precipitation experimental data of crude oil than the previous models,with an absolute average deviation less than 0.4 %.Also,a series of dynamic experiments were carried out to determine the rheological behavior of waxy crude oil in the absence and presence of a flow improver such as ethylene–vinyl acetate copolymer.It was found that the apparent viscosity of waxy crude oil decreased with increasing shear rate.Also,the results showed that the performance of flow improver was dependent on its molecular weight.     

Physical properties of wax deposits on the walls of crude pipelines

Wax deposits on the wall of a crude oil pipeline are a solid wax network of fine crystals, filled with oil, resin, asphaltene and other impurities. In this paper, a series of experiments on wax deposition in a laboratory flow loop were performed under different conditions （flow rate, temperature differential between crude oil and pipeline wall, and dissolved wax concentration gradient）, and the wax deposits were analyzed, so quantitative relationships among wax content, wax appearance temperature （WAT）, shear stress, and radial concentration gradient of dissolved wax at the solid/liquid interface were obtained. Finally, a model was established to predict WAT and the wax content of the deposit.     

Numerical study of crude oil batch mixing in a long channel

The main objective of this work is to predict the mixing of two different miscible oils in a very long channel. The background to this problem relates to the mixing of heavy and light oil in a pipeline. As a first step, a 2D channel with an aspect ratio of 250 is considered. The batch-mixing of two miscible crude oils with different viscosities and densities is modeled using an unsteady laminar model and unsteady RANS model available in the commercial CFD solver ANSYSFluent. For a comparison, a LES model was used for a 3D version of the 2D channel. The distinguishing feature of this work is the Lagrangian coordinate system utilized to set no-slip wall boundary conditions. The global CFD model has been validated against classical analytical solutions. Excellent agreement has been achieved. Simulations were carried out for a Reynolds number of 6300(calculated using light oil properties) and a Schmidt number of 10~4. The results show that, in contrast to the unsteady RANS model, the LES and unsteady laminar models produce comparable mixing dynamics for two oils in the channel. Analysis of simulations also shows that, for a channel length of 100 m and a height of 0.4 m, the complete mixing of two oils across the channel has not been achieved. We showed that the mixing zone consists of the three different mixing sub-zones, which have been identified using the averaged mass fraction of the heavy oil along the flow direction. The first sub-zone corresponds to the main front propagation area with a length of several heights of the channel. The second and third sub-zones are characterized by so-called shear-flow-driven mixing due to the Kelvin–Helmholtz vortices occurring between oils in the axial direction. It was observed that the third sub-zone has a steeper mass fraction gradient of the heavy oil in the axial direction in comparison with the second sub-zone, which corresponds to the flow-averaged mass fraction of 0.5 for the heavy oil.     

Flow pattern identification in oil wells by electromagnetic image logging

Petroleum production logging needs to determine the interpretation models first and flow pattern identification is the foundation, but traditional flow pattern identification methods have some limitations. In this paper, a new method of flow pattern identification in oil wells by electromagnetic image logging is proposed. First, the characteristics of gas-water and oil-water flow patterns in horizontal and vertical wellbores are picked up. Then, the continuous variation of the two phase flow pattern in the vertical and horizontal pipe space is discretized into continuous fluid distribution models in the pipeline section. Second, the electromagnetic flow image measurement responses of all the eight fluid distribution models are simulated and the characteristic vector of each response is analyzed in order to distinguish the fluid distribution models. Third, the time domain changes of the fluid distribution models in the pipeline section are used to identify the flow pattern. Finally, flow simulation experiments using electromagnetic flow image logging are operated and the experimental and simulated data are compared. The results show that the method can be used for flow pattern identification of actual electromagnetic image logging data.     

New dynamic permeability upscaling method for flow simulation under depletion drive and no-crossflow conditions

The main purpose of upscaling in reservoir simulation is to capture the dynamic behavior of fine scale models at the coarse scale.Traditional static or dynamic methods use assumptions about the boundary conditions to determine the upscaled properties.In this paper,we show that the upscaled properties are strongly dependent on the flow process observed at the fine scale.We use a simple no-crossflow depletion drive process and demonstrate that an upscaled property is not a constant value.Instead,if the goal is to match the performance of the fine scale model,the upscaled permeability changes with time.We provide an analytical solution to determine the upscaled permeability and present the value of upscaled permeability under limiting conditions.Our equation suggests that it is possible that upscaled value can fall outside the range of fine scale values under certain conditions.We show that for pseudo steady state flow,using common averaging methods like arithmetic or even geometric averaging methods can lead to optimistic results.We also show that the no-crossflow solution is significantly different than crossflow solution at late times.We validate our method by comparing the results of the method with flow simulation results in two and multi-layered models.     

环烷烃混合烯烃催化裂化：2. 单事件模型评估

The developed SEMK model is used to provide an insight into the contribution of individual reactions in the cracking of methylcyclohexane as well as the site coverage by various carbenium ions. The preferred reaction pathways for the conversion of methylcyclohexane are hydride transfer reactions followed by PCP-isomerizations, deprotonation and endocyclic β-scission, accounting for 61%, 22% and 12% of its disappearance, respectively, at 693 K and 30% conversion of methylcyclohexane. Protolysis plays a minor role in the cracking of methylcyclohexane. Once cyclic diolefins are formed, all of them can be instantaneously transformed to aromatics, which are easily interconverted via disproportionation. Judging from the carbenium ion concentrations it is evident that, at the investigated operating conditions, less than 5% of the acid sites are covered by carbenium ions, less than 2% of which corresponds to cyclic type species including allylic ones.     

Innovative Production of PCMs （Phase Change Materials） Preparation by Vacuum Impregnation： Thermal and Physical Properties

Energy is essential for every human activity for more comfortable life but also consumes more natural resources. In order to control human comfort, temperature usually required when the differences in temperature swing between indoor and outdoor temperatures. PCMs （phase change materials） are the high latent heat materials which can be used in building materials for energy conservation purpose. PCMs can store thermal energy and also can prevent heat to pass through temperature control areas. Paraffin has been used as PCMs are absorbed into the pore of fly ash as paraffin/fly-ash composite and mixed into the buildings materials. Paraffin is an organic material with high melting point （-59℃）, and nonflammable materials therefore paraffin can be used as the building materials for the function of PCMs for energy saving purposes. Composite PCMs can be prepared by vacuum impregnation process. Paraffin in liquid form will be impregnated into the pore of fly ash by vacuum capillary force to form paraffin/fly ash composite PCMs. Vacuum impregnation pressures, vacuum times, impregnation times of liquid paraffin in fly ash pores and temperatures for melting the solid paraffin into the liquid form are all affect on the thermal properties of paraffin/fly ash composite PCMs. Paraffin or PCMs impregnation are also relate to the physical property including the fractal dimensions of the pores of the fly ash particles and paraffin/fly ash composite PCMs. The fractal dimensions of the pore of fly ash and paraffin/fiy ash composites PCMs are between the values of 1.0 and 2.0. Fractal dimensions of paraffin/fly-ash composite PCMs have the same trend as the thermal properties for heat capacity and latent heat of melting. These fractal dimensions technique is a novel method to measure physical property of building material related to latent heat and heat capacity.     

C-Nanostructures Cluster Models in Organic Solvents： Fullerenes,Tubes, Buds and Graphenes

The existence of nanographene in cluster form is discussed in organic solvents. Theories are developed based on the columnlet, bundlet and droplet models describing the size-distribution functions. Phenomena present a unified explanation in the columnlet model in which free energy of Cgraphene involved in cluster is combined from a volume part proportional to the number of molecules n in cluster and a constant. The columnlet model enables describing distribution function of Cgraphene clusters by size. From purely geometrical considerations the columnlet （Cgraphene）, bundlet （single-wall carbon nanotube）, CNT （carbon nanotube）, SWNT （single-wall C-nanotube）, and carbon nanobud, CNB （carbon nanobud）） and droplet （fullerene） models predict dissimilar behaviours. The interaction-energy parameters of Cgraphene are taken from C60. An CNB behaviour or further is expected. The decay of solubility with rising temperature is smaller for Cgraphene than for SWNT and CNB and, furthermore, than for C60, in agreement with lesser numbers of units in Cgraphene clusters. The discrepancy between the experimental data of the heat of solution of fullerenes, CNTs, CNBs and graphenes is ascribed to the sharp concentration dependence of the heat of solution. The diffusion coefficient drops with temperature result greater for Cgraphene than CNB and SWNT than C60 corresponding to lesser number of units in clusters. The aggregates near （C60）13, SWNT/CNB7 and （Cgraphene）3 could be representative of the droplet, bundlet and columnlet models.     

Geochemical Studies and Elemental Contaminants in the Bay of the City of Asuncion

Minor and trace elements composition of bottom sediments from the bay of Asunci6n on the Paraguay River have been investigated by XRF （X-ray fluorescence） techniques to determine their correlation as well as provenance. The analysis of complex spectra was performed by the AXIL software and the quantitative analysis by the QAES （quantitative analysis of environmental samples ） software. Analyzed trace elements were the refractory elements Rb, Sr, Y, Zr, Nb, Ba, La, Ce, Pr, Nd, Zr, Th and others with high field stabilization energy as Cr, Ni, Cu, together with Zn, As, Cd, Pb. Minor elements were Ti, Mn, Fe which are often to the above refractories related. According to their normalized spidergrams, two sets of sediments can be differentiated. Those that show LREE （light rare earth elements） enrichment, negative Nb and Ti anomalies and no spike at Zr and those that have spidergrams very alike, with strong negative anomalies at Nb, Nd and Ti whereas a of Fe versus the refractory elements except Fe-Zr in which correlation seem to be low. positive spike at Zr. In both set, there are strong correlations is negative. Potential hazards of toxic elements in sediments     

脂肪酸及其衍生物对低硫柴油润滑性影响的分子动力学模拟及实验研究

In this work, fatty acid and its derivatives were adopted as lubricity additives for low sulfur diesel. Tribological evaluation obtained from the High-Frequency Reciprocating Rig (HFRR) apparatus showe...     

Innovative Production of PCMs （Phase Change Materials） Preparation by Vacuum Impregnation： Mechanical Strength of Mortars Cement with Composite PCMs Content

An experimental investigation was conducted to study the efficiency of thermal insulation of composite PCMs （phase change materials） produced by vacuum impregnation process between paraffin （PCMs） and fly ash particles. DSC （differential scanning calorimeter） has been used to determine the thermal properties of latent heat of melting and heat capacity for composite PCMs. Vacuum impregnation pressure of 40 in.Hg, paraffin melting temperature of 90℃, vacuum time and impregnation time of paraffin of 30 min are the optimum condition of composite PCMs productions. The values of latent heat of melting and heat capacity are 74.00 J/g and 15.726 J/g.℃ for composite PCMs that produces by the optimum condition in vacuum impregnation process. Increasing the amount of composite PCMs replacing for cement in mortars causes the compressive strength, flexural strength and tensile strength reduction. Compressive strength, flexural strength and tensile strength of mortar with and without composite PCMs can be increased by the longer time of water curing for mortar specimens. Thermal conductivity （k） of mortar cement is reduced by increasing the amount of composite PCMs which replaced for cement in mortar plate compositions. Composite PCMs have the efficiency for thermal energy insulation when incorporated into the buildings. Therefore, this property of paraffin/fly ash composites PCMs can reduce the energy consumption for temperature control in the buildings.     

二次重建法制备具有大孔容高比表面的大孔-介孔γ-Al2O3

Through improving the aging process during synthesis of the support, γ-Al2O3 with large pore volume and high surface area was synthesized by a facile secondary reforming method. The synthesis parameters, such as the reaction temperature, the first aging temperature and the second aging temperature, were investigated. The textural properties of γ-Al2O3 were characterized by means of N2 adsorption-desorption isotherms, X-ray powder diffractometry(XRD), scanning electron microscopy(SEM), Fourier transform infrared(FTIR) spectroscopy and thermogravimetry(TG). The experimental results indicated that AACH and amorphous AlOOH were the precursors of alumina, which were formed via precipitation from solutions after reaction of aluminum sulphate with ammonium hydrogen carbonate. The precursor nanocrystallites grew and re-assembled during the secondary reforming process, which resulted in an increased pore size and pore volume and a decreased bulk density. The as-synthesized γ-Al2O3 materials featured meso/macroporosity, large pore volume(2.175 cm3/g), high surface area(237.8 m2/g), and low bulk density(0.284 g/mL).     

Diatomic Gaseous Sulfur Obtained at Low Temperature Catalytic Decomposition of Hydrogen Sulfide

Experimental evidences of occurrence of gaseous diatomic sulfur produced in the low temperature catalytic decomposition of hydrogen sulfide 2 H2S ←→ 2 H2 ＋ S2 （g） are summarized. The S2 molecule is suggested to be in the ground triplet state. Analysis of literature data allows concluding that the S2 metastable singlet state is realized in the thermal dissociation of hydrogen sulfide and solid sulfur. Arguments in favor of the hypothesis are been discussed.     

水滑石催化PET的合成及性能研究

Poly(ethylene terephthalate) (PET) was synthesized by the in-situ polymerization method using layered double hydrotalcite (LDH) as the catalyst, and the thermal and lfame retardation properties of PET ...     

Wear Behavior of Ductile Cast Irons with Nanoparticles Additives

The work in this study is focused on the investigation of the structure and properties of ductile cast iron with nanoparticle additives： TiN （titanium nitride）, TiN ＋ TiCN （titanium carbonitride） and cBN （cubic boron nitride）. The nanoparticles are coated with nickel prior to addition to the iron melt to improve their wetting and uniform distribution in the volume of the casting. The metallographic observation and wear test are performed to study the influence of the nanoparticle additives on the microstructure and and cast iron tribological properties.