Please wait a minute...
SCI和EI收录∣中国化工学会会刊
[an error occurred while processing this directive]
Table of Content
28 July 2017, Volume 25 Issue 7
    Fluid Dynamics and Transport Phenomena
    Modeling the pressure drop of wet gas in horizontal pipe
    Peining Yu, Yi Li, Jing Wei, Ying Xu, Tao Zhang
    2017, 25(7):  829-837.  doi:10.1016/j.cjche.2016.10.024
    Abstract ( )   PDF (641KB) ( )  
    References | Related Articles | Metrics
    A model for gas-liquid annular and stratified flow through a horizontal pipe is investigated, using the two-phase hydrokinetics theory. Taking into consideration the flow factors including the void fraction, the friction between the two phases and the entrainment in the gas core, the one-dimensional momentum equation for gas has been solved. The differential pressure of the wet gas between the two tapings in the straight pipe has been modeled in the pressure range of 0.1-0.8 MPa. In addition a more objective iteration approach to determine the local void fraction is proposed. Compared with the experimental data, more than 83% deviation of the test data distributed evenly within the band of ±10%. Since the model is less dependent on the specific empirical apparatus and data, it forms the foundation for further establishing a flow measurement model of wet gas which will produce fewer biases in results when it is extrapolated.
    Experimental detection of bubble-wall interactions in a vertical gas-liquid flow
    Xing Wang, Jiao Sun, Jie Zhao, Wenyi Chen
    2017, 25(7):  838-847.  doi:10.1016/j.cjche.2016.11.013
    Abstract ( )   PDF (5750KB) ( )  
    References | Related Articles | Metrics
    Bubble motions and bubble-wall interactions in stagnant liquid were experimentally investigated by high-speed CCD and PIV technique with the main feature parameters such as Eötvös numbers Eo=0.98-1.10, Morton number Mo=3.21×10-9 and Reynolds numbers Re=180~190. The effect of bubble injecting frequency and the distance S between the gas injection nozzle and the wall on the statistical trajectory of bubbles, average velocity distribution of flow field and Reynolds shear stress were studied in detail. It was shown that the combination of bubble injecting frequency and the distance S caused different bubble motion forms and hydrodynamic characteristics. When the normalized initial distance was very little, like S* ≈ 1.2 (here S*=2S/de, and de is the bubble equivalent diameter), bubbles ascended in a zigzag trajectory with alternant structure of high and low speed flow field around the bubbles, and the distribution of positive and negative Reynolds shear stress looked like a blob. With the increase of distance S*, bubbles' trajectory would tend to be smooth and straight from the zigzag curve. Meanwhile, with the increase of bubble injecting frequency, the camber of bubble trajectory at 20 < y < 60 mm had a slight increase due to the inhibitory effect from the vertical wall. Under larger spacing, such as S* ≈ 3.6, the low-frequency bubbles gradually moved away from the vertical plane wall in a straight trajectory and the high-frequency bubbles gradually moved close to the vertical wall in a similar straight trajectory after an unstable camber motion. Under the circumstances, high-speed fluid was mainly distributed in the region between the wall and the bubbles, while the relative large Reynolds shear stress mainly existed in the region far away from the wall.
    Horizontal gas mixing in rectangular fluidized bed:A novel method for gas dispersion coefficients in various conditions and distributor designs
    Asheesh Nautiyal, Chien-Song Chyang, Pin-Wei Li, Hsin-Yung Hou
    2017, 25(7):  848-861.  doi:10.1016/j.cjche.2016.12.009
    Abstract ( )   PDF (1987KB) ( )  
    References | Related Articles | Metrics
    In a rectangular fluidized bed combustor, the tracer gas is injected continuously into the bed from a point source at the center of the distributor plate. In this study, a general governing equation is formulated for tracer gas dispersion in the bed. An analytical solution is derived to estimate the dispersion coefficients, Dx and Dy, in a horizontal plane. The concentration profiles at different sampling heights with various gas velocities are plotted. Subsequently, to estimate the dispersion coefficients, surface fitting of the obtained analytical solution to the experimental data is performed. The dispersion coefficients obtained from this model are compared with those of a conventional model. Additionally, the effect of walls, bed height and gas injection rate on the dispersion coefficients in a horizontal plane is investigated, and the effect of distributor design on the dispersion coefficients in a horizontal plane is investigated with different tracer positions. It is found that Dx and Dy are nearly equivalent at a lower tracer gas ratio of the injected gas to the total gas flow rate. It is also demonstrated that the effect of bed height on Dx is minor. This model is also able to estimate the dispersion coefficients in the case of a multihorizontal nozzle distributor.
    Hydrodynamic dispersion of reactive solute in a Hagen-Poiseuille flow of a layered liquid
    Sudip Debnath, Apu Kumar Saha, B. S. Mazumder, Ashis Kumar Roy
    2017, 25(7):  862-873.  doi:10.1016/j.cjche.2017.03.005
    Abstract ( )   PDF (3821KB) ( )  
    References | Related Articles | Metrics
    An analysis of the solute dispersion in the liquid flowing through a pipe by means of Aris-Barton's ‘method of moments’, under the joint effect of some finite yield stress and irreversible absorption into the wall is presented in this paper. The liquid is considered as a three-layer liquid where the center region is Casson liquid surrounded by Newtonian liquid layer. A significant change from previous modelling exercises in the study of hydrodynamic dispersion, different molecular diffusivity has been considered for the different region yet to be constant. For all time period, finite difference implicit scheme has been adopted to solve the integral moment equation arising from the unsteady convective diffusion equation. The purpose of the study is to find the dependency of solute transport coefficients on absorption parameter, yield stress, viscosity ratio, peripheral layer variation and in addition with various diffusivity coefficients in different liquid layers. This kind of study may be useful for understanding the dispersion process in the blood flow analysis.
    Gas-liquid two-phase flow in serpentine microchannel with different wall wettability
    Yunlong Zhou, He Chang, Tianyu Qi
    2017, 25(7):  874-881.  doi:10.1016/j.cjche.2016.10.006
    Abstract ( )   PDF (905KB) ( )  
    References | Related Articles | Metrics
    Gas-liquid flow in serpentine microchannel with different surface properties exhibits drastically different flow behavior. With water and air as working fluids, the method of numerical simulation was adopted in this paper based on CLSVOF (coupled level set and volume of fluid method) multiphase model. After verifying the reasonability of the model through experiment, by changing wall properties and Re number (Re < 1500), the influences of contact angle and surface roughness on flow regime and Po number were discussed. Moreover, the difference of pressure drop between curve and straight microchannel was also calculated. Beyond that, the combined effect of curve channel and wall properties on flow resistance was analyzed. This paper finds that wall properties have great influence on gas-liquid flow in microchannels not only on flow regime but also flow characteristics. Meanwhile, the pressure drop in curve microchannels is larger than straight. It is more beneficial for fluid flowing when the straight part of microchannel is hydrophilic smooth wall and curve part is hydrophobic with large roughness.
    Separation Science and Engineering
    Enhanced permeation performance of polyether-polyamide block copolymer membranes through incorporating ZIF-8 nanocrystals
    Longwei Xu, Long Xiang, Chongqing Wang, Jian Yu, Lixiong Zhang, Yichang Pan
    2017, 25(7):  882-891.  doi:10.1016/j.cjche.2016.11.007
    Abstract ( )   PDF (1067KB) ( )  
    References | Related Articles | Metrics
    Membrane-based CO2 separation is a promising alternative in terms of energy and environmental issues to other conventional techniques. Polyether-polyamide block copolymer (Pebax) membranes are promising for CO2 separation because of their excellent selectivity, but limited by their moderate gas permeability. In this study, fresh-prepared zeolitic imidazolate framework-8 (ZIF-8) nanocrystals were integrated into the Pebax®1657 matrices to form mixed matrix membranes. The resulting membrane exhibits significantly improved CO2 permeability (as high as 300% increase), without the sacrifice of the selectivity, to the pristine polymer membrane. Several physical characterization techniques were employed to confirm the good interfacial interaction between ZIF-8 fillers and Pebax matrices. The effect of added ZIF-8 fillers on the transport mechanism through MMMs is also explored. Mixed-gas permeation for both CO2/N2 and CO2/CH4 was also evaluated. The separation performance for CO2/CH4 mixtures on the ZIF-8/Pebax MMMs is very close to the Roberson upper bound, and thus is technologically attractive for purification of natural gas.
    Modeling of U-shaped Ba0.5Sr0.5Co0.8Fe0.2O3-δ hollow-fiber membrane for oxygen permeation
    Huiqi Xie, Yanying Wei, Haihui Wang
    2017, 25(7):  892-897.  doi:10.1016/j.cjche.2017.02.002
    Abstract ( )   PDF (693KB) ( )  
    References | Related Articles | Metrics
    A mathematic model is developed for the perovskite-type mixed ionic-electronic conducting (MIEC) membrane, which makes it possible to simulate the process of oxygen separation in the U-shaped Ba0.5Sr0.5Co0.8Fe0.2O3-δ hollow-fiber membrane. The model correlates the oxygen permeation flux to the measurable variables. The trends of calculated results for purge operation coincide well with the experimental data, therefore the model is considerable for flux prediction under vacuum operation. Higher oxygen separation efficiency can be achieved with vacuum operation than purge operation. Parameter study with vacuum operation reveals that oxygen permeation flux increases with higher vacuum levels, and vacuum pressure of around 1.013×103 Pa is the optimal. Also, vacuum operation on the lumen side is much more efficient to achieve higher oxygen permeation flux compared with compression mode on the shell side.
    Preparation and evaluation of mixed-mode resins with tryptophan analogues as functional ligands for human serum albumin separation
    Qici Wu, Qilei Zhang, Shiwen Xu, Cheng ong Ge, Shanjing Yao, Dongqiang Lin
    2017, 25(7):  898-905.  doi:10.1016/j.cjche.2016.12.010
    Abstract ( )   PDF (793KB) ( )  
    References | Related Articles | Metrics
    Five tryptophan analogues with a hydrophobic indole ring and an amino group on each molecule were used as functional ligands of mixed-mode resins for human serum albumin (HSA) purification. Their adsorption performance was evaluated and the effects of pH and salt addition on HSA adsorption were studied. The resins prepared showed typical pH-dependent adsorption and the highest adsorption capacity and affinity were found at pH 5.0 for all the resins tested. The saturated adsorption capacity was 138.02 mg·g-1 with the tryptaminefunctionalized resin, which significantly decreased at pH below 4.0 due to electrostatic repulsion between ligands and HSA. Moreover, the addition of NaCl or (NH4)2SO4 in media reduced HSA adsorption capacity, although the two salts showed different affecting profiles. The tryptamine-functionalized resin showed the best salt-tolerant performance, and its high adsorption capacity was maintained under high salt concentrations. In addition, the five resins prepared showed good adsorption selectivity for recombinant HSA from Pichia pastoris broth. Molecular docking results between tryptamine and HSA indicated that tryptamine was favorable to bind on Site Ⅱ (indole-binding site) of HSA.
    Implications from protein uptake kinetics onto dextran-grafted Sepharose FF coupled with ion exchange and affinity ligands
    Aiying Xue, Linling Yu, Yan Sun
    2017, 25(7):  906-910.  doi:10.1016/j.cjche.2017.02.003
    Abstract ( )   PDF (346KB) ( )  
    References | Related Articles | Metrics
    Our previous studies have reported the presence of "chain delivery" effects of protein adsorption onto ion exchangers with polymer-grafted ion-exchange groups, such as dextran-grafted and poly(ethylenimine)-modified Sepharose gels. However, it is unclear if the "chain delivery" occurs on affinity adsorption with specific interactions. This work is designed to address this issue. A dextran-grafted Sepharose gel was prepared, and then the matrix was modified using diethylaminoethyl, a typical ion-exchange group, or octapeptide (FYCHWQDE), an affinity ligand for human immunoglobulin G (hIgG) to prepare ion-exchange or affinity adsorbents, respectively. Results of hIgG adsorption showed that the uptake rate represented by the effective diffusivity of hIgG onto the dextran-grafted ion exchangers was obviously enhanced by the dextran grafting, indicating the presence of "chain delivery" of the bound proteins on the charged groups on the dextran chains. By contrast, the effective diffusivity of hIgG changed little as ligand density increased on the dextran-grafted FYCHWQDE adsorbents. Their adsorption capacities decreased and effective diffusivities were not accelerated by the dextran grafting. Thus, this work clarified that grafted dextran could not accelerate hIgG uptake rate on the affinity resins, or in other words, chain delivery did not occur on the specific interaction-based affinity adsorption.
    Effects of solubility parameter differences among PEG,PVP and CA on the preparation of ultrafiltration membranes:Impacts of solvents and additives on morphology,permeability and fouling performances
    Kibrom Alebel Gebru, Chandan Das
    2017, 25(7):  911-923.  doi:10.1016/j.cjche.2016.11.017
    Abstract ( )   PDF (1096KB) ( )  
    References | Related Articles | Metrics
    The effects of two different hydrophilic additives and two solvents on the membrane morphological structure, permeability property and anti-fouling performances of cellulose acetate (CA) ultrafiltration membranes were investigated. During the phase-inversion process, cellulose acetate was selected as a membrane forming polymer; polyethylene glycol (PEG) and polyvinyl pyrrolidone (PVP) were used as additives; acetone (Ac):N, N-Dimethylacetamide (DMAc) and N, N-Dimethylformamide (DMF) were used as solvents; and deionized (DI) water was used in the coagulation bath. All the prepared membranes were characterized in terms of hydraulic permeability (Pm), membrane resistance, average pore radius, and hydrophilicity. The top surface and crosssectional view of the prepared membranes were also observed by using field emission scanning electron microscopy. Membrane fouling and rejection experimentations were done using a stirred batch-cell filtration set-up. The experimental studies of fouling/rinsing cycles, rejection, and permeate fluxes were used to investigate the effect of PEG and PVP additives and effect of the two solvents on the fabricated membranes using bovine serum albumin (BSA) as a model protein.
    Synthesis of clay-supported nanoscale zero-valent iron using green tea extract for the removal of phosphorus from aqueous solutions
    Akbar Soliemanzadeh, Majid Fekri
    2017, 25(7):  924-930.  doi:10.1016/j.cjche.2016.12.006
    Abstract ( )   PDF (508KB) ( )  
    References | Related Articles | Metrics
    This study addresses the synthesis of nanoscale zero-valent iron (nZVI) in the presence of natural bentonite (B-nZVI) using green tea extract. The natural bentonite and B-nZVI were then applied for the removal of phosphorus from aqueous solutions at various concentrations, pH levels and contact time. The desorption of phosphorus (P) from adsorbents was done immediately after sorption at the maximum initial concentration using the successive dilution method. The characterization of FTIR, SEM, and XRD indicated that nZVI was successfully loaded to the surface of natural bentonite. The sorption of phosphorus on B-nZVI was observed to be pH-dependent, with maximum phosphorus removal occurring at the pH range of 2 to 5. The results demonstrate that the maximum sorption capacities of natural bentonite and B-nZVI were 4.61 and 27.63 mg·g-1, respectively. Langmuir, Freundlich, and Redlich-Peterson models properly described the sorption isotherm data. For either adsorbent, desorption isotherms did not coincide with their corresponding sorption isotherms, suggesting the occurrence of irreversibility and hysteresis. The average percentages of retained phosphorus released from natural bentonite and B-nZVI were 80% and 9%, respectively. The results indicated that sorption kinetics was best described by the pseudo-second-order model. The present study suggests that B-nZVI could be used as a suitable adsorbent for the removal of phosphorus from aqueous solutions.
    Suppression of gold nanoparticle agglomeration and its separation via nylon membranes
    Ayyavoo Jayalakshmi, In-Chul Kim, Young-Nam Kwon
    2017, 25(7):  931-937.  doi:10.1016/j.cjche.2017.01.009
    Abstract ( )   PDF (716KB) ( )  
    References | Related Articles | Metrics
    Use of ultraporous nylon membrane is one of the most widely employed techniques for removal of hard and soft nanoparticles in the semiconductor industry, and the accurate determination of membrane pore size is necessary in order to avoid manufacturing defects caused by contamination. The gold nanoparticle has several benefits for the evaluation of polymeric membranes; however, the nanoparticles agglomerate easily on the nylon membrane and make it difficult to evaluate the membrane precisely. The properties of 2-amino-2-hydroxymethyl-1,3-propanediol (ADP) ligand in gold nanoparticle solution were systematically investigated, and ADP was utilized for improved evaluation of the nylon membranes. Nylon membrane used in this study was prepared by phase inversion techniques. Ultrathin dense layer on top of the membrane surface and Darcy structures in the microporous membrane support were observed. The gold particle rejection was carried out at various pH values from 4 to 14 and higher rejection was observed at pH 4 and 8. The suppression of gold colloid agglomeration using ADP and monodispersity of gold colloids was also analyzed by confocal laser scanning microscopy (CLSM), transmission electron microscopy (TEM), and scanning electron microscopy (SEM). van der Waals interaction energy of the particles was reduced in the addition of ADP. The presence of ADP ligand in the gold solutions prevented the agglomeration of gold nanoparticles and reduced the adsorption of the particles on the nylon membrane surface, leading to precise evaluation of membrane pore sizes.
    Interpenetrating polymers supported on microporous polypropylene membranes for the transport of chromium ions
    Yesid Tapiero, Julio Sánchez, Bernabé L. Rivas
    2017, 25(7):  938-946.  doi:10.1016/j.cjche.2016.12.011
    Abstract ( )   PDF (684KB) ( )  
    References | Related Articles | Metrics
    Modifying polypropylene membranes with interpenetrating polymer networks (IPNs) through the incorporation of poly(glycidyl methacrylate-N-methyl-D-glucamine) (P(GMA-NMG)) was performed by in situ synthesis via radical polymerization. The surface of the polypropylene membrane was activated by hydrophilic grafted polyelectrolyte, and then, pressure injection was used for the impregnation of the reactive solution in the membrane. Two types of pore-filled membranes were synthesized, chelating interpenetrating homopolymer networks of P(GMA-NMG), and chelating-ion exchange interpenetrating polymer networks (e.g., P(GMA-NMG)/P(AA), P(GMA-NMG)/P(AMPSA), and P(GMA-NMG)/P(ClVBTA)). After their synthesis, the modified polypropylene membranes were characterized using techniques such as the electrokinetic potential, SEM, FT-IR, and Donnan dialysis to corroborate the chromium ion transport. The P(GMA-NMG) and complex network membranes exhibited a hydrophilic character with a water-uptake capacity between 20% and 35% and a percentage of modification between 4.0% and 7.0% in comparison with the behavior of the unmodified polypropylene membrane. Hexavalent chromium ions were efficiently transported from the food chamber at pH 9.0 when the 65.2% MTA1 P(ClVBTA) homopolymer IPN membrane and 48.5% MTAG P(GMA-NMG)/P(ClVBTA) IPN membrane were used. Similarly, hexavalent chromium ions were removed from the food chamber at pH 3.0 when MTAG (63.30%) and MTA1 (35.68%) were used in 1 mol·L-1 NaCl solution as the extraction reagent.
    Elevating the flexibility and operability of dividing-wall distillation columns via feed thermal condition adjustment
    Yang Yuan, Kejin Huang, Liang Zhang, Haisheng Chen, Shaofeng Wang, Yingjie Jiao
    2017, 25(7):  947-954.  doi:10.1016/j.cjche.2017.03.011
    Abstract ( )   PDF (716KB) ( )  
    References | Related Articles | Metrics
    Because of the complicated interplay between the prefractionator and main distillation column involved, the black-hole problem might occur and prohibit the assignment of four specifications to dividing-wall distillation columns (DWDCs) (e.g., the three main product compositions plus an impurity ratio in the intermediate product), which lowers terribly process flexibility and operability. In this paper, a feed thermal condition adjustment strategy, achieved by the installation of a pre-heater in feed pipeline, is employed to eliminate the black-hole problem and serve to enhance process flexibility and operability. Through the strong influence to the overall mass and energy balance of the DWDC, the feed thermal condition adjustment can alter the interlinking flows between the thermally coupled prefractionator and main distillation column and work effectively to coordinate their relationship. A DWDC separating a benzene, toluene, and o-xylene mixture is chosen to ascertain the feasibility of the philosophy proposed. The static and dynamic studies demonstrate that the feed thermal condition adjustment is an effective way to refine process design and can completely eliminate the black-hole problem and elevate consequently process flexibility and operability.
    Chemical Engineering Thermodynamics
    Phase equilibria for the pseudo-ternary system (NaCl+Na2SO4+H2O) of coal gasification wastewater at T=(268.15 to 373.15) K
    Haijiao Lu, Jingkang Wang, Jun Yu, Yuefeng Wu, Ting Wang, Ying Bao, Dou Ma, Hongxun Hao
    2017, 25(7):  955-962.  doi:10.1016/j.cjche.2016.08.016
    Abstract ( )   PDF (1778KB) ( )  
    References | Related Articles | Metrics
    The pseudo-ternary system (NaCl + Na2SO4 + H2O) of coal gasification wastewater was studied at T=(268.15 to 373.15) K. The solubility and density of the equilibrium liquid phase were determined by the isothermal solution saturation method. The equilibrium solids were also investigated by the Schreinemaker's method of wet residues and X-ray powder diffraction (XRD). According to the experimental data, the phase diagrams were determined. It was found that there was no significant solubility difference on the NaCl-rich side between the ternary system (NaCl + Na2SO4 + H2O) in coal gasification wastewater and in pure water. However, the solubility on the Na2SO4-rich side of coal gasification wastewater was apparently higher than that of pure water. The increase in the solubility of Na2SO4 was most likely caused by the effects of other impurities apart from NaCl and Na2SO4 in coal gasification wastewater. The measured data and phase equilibrium diagrams can provide fundamental basis for salt recovery in coal gasification wastewater.
    Non-isothermal crystallization kinetics of Nylon 10T and Nylon 10T/1010 copolymers:Effect of sebacic acid as a third comonomer
    Zhongqiang Wang, Guosheng Hu, Jingting Zhang, Jiusheng Xu, Wenbo Shi
    2017, 25(7):  963-970.  doi:10.1016/j.cjche.2016.06.006
    Abstract ( )   PDF (702KB) ( )  
    References | Related Articles | Metrics
    Nylon 10T and Nylon 10T/1010 samples were synthesized by direct melt polymerization. The non-isothermal crystallization kinetics of Nylon 10T and Nylon 10T/1010 was investigated by means of differential scanning calorimetry (DSC). Jeziorny equation and Mo equation were applied to describe the non-isothermal crystallization kinetics of the Nylon 10T and the Nylon 10T/1010. The activation energies for non-isothermal crystallization were obtained by Vyazovkin's method and Friedman's method, respectively. These results showed that Jeziorny equation and Mo equation well described the non-isothermal crystallization kinetics of the Nylon 10T and the Nylon 10T/1010. It was found that the values of the activation energy for non-isothermal crystallization of the Nylon 10T/1010 were lower than those of the Nylon 10T at a given temperature or relative crystallinity degree, which revealed that crystallization ability of the Nylon 10T/1010 was higher. The crystal morphology was observed by means of a polarized optical microscope (POM) and X-ray diffraction (XRD). It was found that the addition of sebacic acid comonomer not only did not change the crystal form of the Nylon 10T, but also significantly increased the number and decreased the size of spherulites. Comparing with the Nylon 10T, the crystallization rate was increased with the addition of the sebacic acid comonomer.
    Enthalpy of phase transition of isonicotinic acid
    Dongfang Zhao, Guanghui Liu, Jian Sun, Lisheng Wang
    2017, 25(7):  971-975.  doi:10.1016/j.cjche.2016.08.014
    Abstract ( )   PDF (431KB) ( )  
    References | Related Articles | Metrics
    In this work, the group contribution method of Chickos et al. was applied to estimate the fusion enthalpy of isonicotinic acid, and the obtained result (29.2 kJ·mol-1) showed a large difference with the value (135 kJ·mol-1) as reported from literatures and as determined by differential scanning calorimetry (DSC). The results of DSC/TG measurement showed that the phase transition of isonicotinic acid from 187.27℃ to 277.47℃ underwent a sublimation process, with a sublimation enthalpy of 128.03 kJ·mol-1. An efficient analytical technique combining pyrolysis and gas chromatography/mass spectrometry (Py-GC/MS) was used to prove this conclusion.
    Viscosity profile prediction of a heavy crude oil during lifting in two deep artesian wells
    Yibo Li, Hao Gao, Wanfen Pu, Bing Wei, Yafei Chen, Dong Li, Qiang Luo
    2017, 25(7):  976-982.  doi:10.1016/j.cjche.2016.11.014
    Abstract ( )   PDF (1138KB) ( )  
    References | Related Articles | Metrics
    It has been known that the productivity of artesian wells is strongly dependent on the rheological properties of crude oils. This work targets two deep artesian wells (>5000 m) that are producing heavy crude oil. The impacts of well conditions including temperature, pressure and shear rate, on the crude oil rheology were comprehensively investigated and correlated using several empirical rheological models. The experimental data indicate that this heavy oil is very sensitive to temperature as result of microstructure change caused by hydrogen bonding. The rheological behavior of the heavy oil is also significantly impacted by the imposed pressure, i.e., the viscosity flow activation energy (Eμ) gently increases with the increasing pressure. The viscosity-shear rate data are well fitted to the power law model at low temperature. However, due to the transition of fluid feature at high temperature (Newtonian fluid), the measured viscosity was found to slightly deviate from the fitting data. Combining the evaluated correlations, the viscosity profile of the heavy crude oil in these two deep artesian wells as a function of well depth was predicted using the oilfield producing data.