Increase In Flow Resistance Research Articles

Experimental and numerical study of forming force in ultrasonic vibration single point incremental forming

A device combining ultrasonic vibration plastic processing technology and single point incremental forming was developed in this work, with maximum power output of 2 kW at frequency of 20 kHz. On this basis, together with experiments and finite element simulation on forming force with aluminum alloy sheet 1060, the results showed that the forming force reduced evidently by appropriate selection of amplitude and frequency, meanwhile, the change of them had a remarkable influence on plastic behavior which can be summarized into two opposite aspects: the softening and the hardening effect. When a lower amplitude and frequency was applied to the forming process, the softening effect dominated, leading to a decrease of forming force, however, under the application of a high-vibrating amplitude and frequency, the harden effect dominated resulting in the decline of plasticity and the increase of flow resistance.

Significance of thermal radiation, Lorentz force, and non-Darcian porous medium on the dynamics of second-grade fluid subject to exponential stretching sheet

ABSTRACT This paper studies the non-Newtonian fluid flow over an exponential non-Darcian permeable stretching sheet. The heat generation and non-linear radiation effects are inspected for second-grade fluid flow by employing MHD normal to the flow. The mathematical model is developed in partial differential equations under the boundary layer approximations. These equations are reduced by employing similarity transformations in ordinary differential equations. The ordinary differential equations in dimensionless forms are solved by using the numerical procedure bvp4c method. The involved physical parameters impacts are explained via graphs and tables. The horizontal velocity function enlarges the porosity parameter because the porous media behaves as a barrier (opposes) to the flow. The vertical velocity diminishes with higher values of variable viscosity parameters because the rising viscosity is associated with increasing flow resistance, which lessens the fluid's velocity. The vertical velocity gets lower values as the values of magneto hydrodynamic parameter rise, causing a reduction in flow velocity, and as a result, graph vertical velocity falls. However, as the magnetohydrodynamic parameter increases, liquid surface penetration in the wall region rises.

Spent fuel simulation during dry storage via enhancement of FRAPCON-4.0: Comparison between PWR and SMR and discharge burnup effect

Spent fuel behavior of dry storage was simulated in a continuous state from steady-state operation by modifying FRAPCON-4.0 to incorporate spent fuel–specific fuel behavior models. Spent fuel behavior of a typical PWR was compared with that of NuScale Power Module (NPM™). Current PWR discharge burnup (60 MWd/kgU) gives a sufficient margin to the hoop stress limit of 90 MPa. Most hydrogen precipitation occurs in the first 50 years of dry storage, thereby no extra phenomenological safety factor is identified for extended dry storage up to 100 years. Regulation for spent fuel management can be significantly alleviated for LWR-based SMRs. Hydride embrittlement safety criterion is irrelevant to NuScale spent fuels; they have sufficiently lower plenum pressure and hydrogen contents compared to those of PWRs. Cladding creep out during dry storage reduces the subchannel area with burnup. The most deformed cladding outer diameter after 100 years of dry storage is found to be 9.64 mm for discharge burnup of 70 MWd/kgU. It may deteriorate heat transfer of dry storage by increasing flow resistance and decreasing the view factor of radiative heat transfer. Self-regulated by decreasing rod internal pressure with opening gap, cladding creep out closely reaches the saturated point after ∼50 years of dry storage.

Experimental Study on Matched Particle Size and Elastic Modulus of Preformed Particle Gel for Oil Reservoirs

Suitable elastic modulus and particle size of preformed particle gel are the keys to both diverting water flow and avoiding permanent impairment to reservoirs. Therefore, the paper aims at finding the best matched preformed particle gel for given reservoirs using sand-pack displacement experiments. The results show that the injection pressure of preformed particle gel with excessively small size and elastic modulus is relatively low, indicating poor capacity to increase flow resistance and reduce water channeling. On the other hand, if the particle size and elastic modulus of preformed particle gel are excessively large, the reservoir may be plugged and irreversibly damaged, affecting oil development performance. In fact, the best matched particle size and elastic modulus of preformed particle gel increase with the increase in reservoir permeability. Furthermore, the paper establishes a quantitative logarithmic model between the particle size of preformed particle gel and reservoir permeability. Finally, the established matching relationship is validated via microscopic visualization oil displacement experiments using a glass etching model. The validation experiments indicate that the preformed particle gel (60–80 mesh; 2–4 Pa) selected according to the matching relationship can effectively reduce water channeling and increase sweeping efficiency by as much as 55% compared with water flooding in the glass etching model with an average permeability of 2624 × 10−3 μm2. Therefore, the established matching relationship can provide an effective guide when selecting the best suitable preformed particle gel for a given reservoir in more future applications.

CFD of roughness effects on laminar heat transfer applied to additive manufactured minichannels

AbstractAdditive manufacturing has received significant interest in the fabrication of functional channels for heat transfer; however, the inherent rough surface finish of the additively manufactured channels can influence thermal performance. This study investigates the impact of roughness on the thermo-fluid characteristics of laminar forced convection in rough minichannels. A numerical model was developed to create 3D Gaussian roughness with specified root-mean-square height. The finite volume method was used to solve the conjugate heat transfer of developed laminar flow in square minichannels. For Reynolds numbers ranging from 200 to 1600, the simulation results indicated enhanced heat transfer and increased flow resistance as Reynolds number increases, compared to a smooth minichannel, where effects on heat transfer and flow friction were negligible. For channels with relative roughness (root-mean-square height to channel hydraulic diameter) of 0.0068, 0.0113, and 0.0167, increasing the Reynolds number led to increased friction factor by 1.56, 1.71, and 2.91%, while the Nusselt number was enhanced up to 0.03%, 32.74%, and 46.05%, respectively. Heat transfer reduced in roughness valleys due to the presence of local low-velocity fluid in these regions; however, recirculation regions can occur in deep valleys of high roughness, increasing heat transfer and flow friction. Heat transfer was enhanced over roughness peaks due to flow impingement on the windward face of roughness as well as intensified energy transfer to the channel wall from roughness. Moreover, surfaces with higher roughness have a greater number of high peaks providing a thermal-flow path of a larger area and a thermal conductivity greater than that of the fluid.

Analysis of Fire Throttling in Longitudinally Ventilated Tunnels With a One-dimensional Model

Fire throttling is the increase in flow resistance due to a large fire in a longitudinally ventilated tunnel. Although the fire throttling effect has been been known and studied for tunnels over the last 40 years, there is not yet a consistent one-dimensional (1D) model that can describe this behaviour or a framework suitable for practical application. We propose a semi-empirical model, based upon pipe flow engineering principles, to describe this effect by separating the resistance to flow, or pressure loses in three parts: upstream of the fire, locally at the fire, and downstream of the fire. The proposed 1D model called TE1D is derived from a simple steady one-dimensional momentum balance in which a semi-empirical mean temperature distribution is assumed across the tunnel. We verify the model by comparing the pressures losses it predicts with those calculated in CFD simulations based on OpenFOAM and Fire Dynamics Simulator. The comparison shows good agreement between the CFD codes for the range of fires sizes considered from 5 to 50 MW and good agreement between TE1D and the CFD results with the proposed 1D model for fire sizes below 30 MW. However, for values above there are large discrepancies between the results obtained by the TE1D and CFD. We posit as a potential explanation that these differences are due to flow and temperature stratification which is not accounted for in the 1D model. The model using pipe flow principles allows engineers to adopt this model for design, together with other pressure losses considered in tunnel ventilation.

Influence of Pore Morphology on Permeability through Digital Rock Modeling: New Insights from the Euler Number and Shape Factor

The microscopic pore shape and topology significantly affect fluid transport and occurrence in porous permeable rock. A quantitive characterization of the impact of pore morphology on permeability is currently lacking, which limits the efficient development of underground hydrocarbon resources. This work introduces the Euler number and shape factor to characterize the pore topology and shape of heterogeneous sandstone based on CT imaging. The pore morphology under different pore sizes and the correlation of the Euler number, shape factor, fractal dimension, and surface area are analyzed. Furthermore, a modified Kozeny–Carman equation is established to explain the influence of the Euler number and shape factor on permeability. The results show that with the increase of pore diameter, the Euler number decreases while the shape factor increases. In a connected pore system, the smaller Euler number corresponds to the complex pore network, which leads to the increase in the surface area, shape factor, and fractal dimension. At constant porosity, the shape factor is negatively correlated with permeability, and with increasing Euler number, the heterogeneity of the pore structure increases, resulting in an increase of flow resistance and a decrease of permeability. The results provide a new pore morphology characterization method for digital rock and help to understand the flow mechanism of hydrocarbons in complex pore networks.

Hydraulic fractures simulation in non-uniform stress field of horizontal shale gas well

Shale reservoirs have low porosity and ultra-low permeability. A significant increase in gas well production requires horizontal well fracturing. In this view, the propagation behavior of hydraulic fractures during the fracturing process is critical to the effect of shale gas reservoir stimulation. However, some shale formations have developed faults and the in-situ stress field has a non-uniform distribution. Also, the maximum horizontal principal stress direction varies with the spatial position, and non-planar deflection propagation behavior in hydraulic fractures is possible. The purpose of this article is to reveal the effect of the development of faults in the shale formation on the propagation of hydraulic fractures. A calculation model for the non-uniform tectonic stress field near the fault was established, and on this basis, we constructed a shale gas horizontal well-staged multi-cluster fracturing hydraulic fracture deflection propagation model under the influence of faults. Furthermore, the fracture propagation law was simulated and analyzed in the presence of a non-uniform in-situ stress field. The findings demonstrate that hydraulic fractures near the strike-slip fault tend to run parallel to the fault direction, whereas hydraulic fractures near the normal fault tend to run perpendicular to the fault direction. The greater the fault throw and the greater the hydraulic fracture propagation deflection angle, the closer the hydraulic fracture is to the fault. The deflection behavior of hydraulic fractures increases flow resistance in the fracture, and proppant transport does not flow smoothly, increasing the difficulty and risk of the fracturing operation. The findings of this study could provide an important theoretical foundation for the design of staged multi-cluster fracturing of horizontal wells near faults in shale formations.

Investigating the Behavior of Waste Alumina Powder and Nylon Fibers for Eco-Friendly Production of Self-Compacting Concrete.

Self-compacting concrete (SCC) incorporating secondary raw materials has been extensively used around the globe due to its improved fresh, mechanical and durability properties. This study was planned to evaluate the suitability of locally available waste alumina powder (AP) and nylon textile fibers (NF) as a partial replacement for fine and coarse aggregates with the ultimate goal to locally produce SCC with desired properties. The used AP was acquired from a local market and NF was collected from a local textile factory. Various dosages of AP (10%, 20%, 30%, 40% and 50% by volume of fine aggregates) and NF (1% and 2% by volume of coarse aggregates) were studied. Tests including slump flow, V-funnel and J-ring tests were performed for examining the fresh properties of developed SCC. Results showed that the addition of AP has an insignificant effect on the superplasticizer dosage for maintaining a constant flow of 70 cm. However, a higher dosage of superplasticizer was required for a mixture with increasing dosages of NF to sustain a constant flow. Similarly, slump flow time (for a spread of 50 cm) and V-funnel time increased for mixtures with higher dosages of AP and NF. Tested SCC mixtures incorporating 40% and 50% of AP with 1% and 2% of NF showed an extreme blocking assessment due to their increased interparticle friction, the higher water absorption capacity of used AP and NF leading to increased flow resistance and hence, showed lower passing ability. The compressive strength was 16% higher for specimens incorporating 40% of AP due to the filling effect of AP which fills the micro-pores, leading to a more dense and compact internal micro-structure, confirmed through scanning electron microscopy analysis. An ultrasonic pulse velocity test conducted on hardened specimens verified the findings of the compressive strength results. Moreover, it was observed that NF has an insignificant effect on the compressive strength; however, flexural strength was increased due to the incorporation of NF, especially at higher dosages of AP.

The effect of different twisted tape inserts configurations on fluid flow characteristics, pressure drop, thermo-hydraulic performance and heat transfer enhancement in the 3D circular tube

The thermal-hydraulic performance of the circular pipe heat exchanger is determined by numerical simulations of the flow field and pressure drop inside the pipe. Initially, a three-dimensional numerical model is validated using experimental results. Turbulent flows are solved using the governing, momentum and energy equations. With constant heat flux as the boundary condition, the (RNG) k-ω turbulent model predicts flow structures. Different twisted tape inserts (NTTI) (1, 3 and 5) are considered when analysing the influence of different twisted tape geometric parameters. A variety of important parameters are compared quantitatively and qualitatively using six different twisted turns (NTT): static pressure, velocity magnitude, vorticity and static temperature. Twisted tape inserts in pipes can increase flow resistance, which leads to an increase in pressure difference. Furthermore, thetwisted tapes inside the pipe caused more vortex motion (swirl flows), which resulted in different radial velocities. As compared to the temperature difference in a smooth pipe, the numerical results show that the temperature difference increases up to 38.1%, 46.11% and 50.52% when the NTTI is increased from 1 to 5.

Modeling and Analysis of the Flow Characteristics of Liquid Hydrogen in a Pipe Suffering from External Transient Impact

Pipes can be subjected to external transient impacts such as accidental collision, which affects the safe operation of storage and transportation systems for liquid hydrogen. Fluid–structure coupling calculation for a pipe under external transient impact is performed, and the flow characteristics of liquid hydrogen in the pipe are analyzed. The pipe deforms and vibrates when suffering from external transient impact. Liquid hydrogen pressure in a cross-section plane increases along the pipe deformation direction. Additionally, external transient impact enhances the disturbance of liquid hydrogen near the pipe wall. The increased flow resistance and the energy induced by the deformed pipe both affect the flow of liquid hydrogen, and contribute to the fluctuated characteristics of liquid pressure drop. In addition, the phase state of liquid hydrogen remains unchanged in the pipe, indicating that little of the induced energy is transformed into the internal energy of liquid hydrogen. The work provides theoretical guidance for the safe operation of liquid hydrogen storage and transportation systems.

Numerical investigation on hydraulic and thermal performances of a mini-channel heat sink with twisted ribs

Micro/mini-channel heat sink is an effective cooling device for electronic components. To improve its comprehensive performance, several mini-channels with twisted ribs are designed and their hydraulic and thermal performances are investigated using numerical method under Reynolds number of 239–954. The results are compared with the mini-channels with plain ribs and with inclined ribs, and the effects of top twisted angle (0°–45°), bottom twisted angle (0°–45°) and height (0.5–1.0 mm) of twisted rib on performances are analyzed. Results show that twisted rib or inclined rib in mini-channel enhances heat transfer but also increases flow resistance compared to plain rib. Although the mini-channel with inclined ribs shows best thermal performance, the mini-channel with special twisted rib performs best comprehensive performance. Increasing top or bottom twisted angle or decreasing twisted rib height can improve thermal and comprehensive performances, as well as the flow resistance increases with the top or bottom twisted angle or twisted rib height increasing. The mini-channel with twisted rib with bottom twisted angle of 45°, top twisted angle of 0° and rib height of 0.5 mm yields the highest performance evaluation criterion among the mini-channels with ribs at the studied range of Reynolds number. The comparison of performance evaluation criterion between present study and several other studies suggests that, twisted rib structure has great potential to improve comprehensive performance of mini-channel at low Reynolds number.

Amplification of downstream flood stage due to damming of fine-grained rivers

River dams provide many benefits, including flood control. However, due to constantly evolving channel morphology, downstream conveyance of floodwaters following dam closure is difficult to predict. Here, we test the hypothesis that the incised, enlarged channel downstream of dams provides enhanced water conveyance, using a case study from the lower Yellow River, China. We find that, although flood stage is lowered for small floods, counterintuitively, flood stage downstream of a dam can be amplified for moderate and large floods. This arises because bed incision is accompanied by sediment coarsening, which facilitates development of large dunes that increase flow resistance and reduce velocity relative to pre-dam conditions. Our findings indicate the underlying mechanism for such flood amplification may occur in >80% of fine-grained rivers, and suggest the need to reconsider flood control strategies in such rivers worldwide.

Assessment of Coal Permeability using Injection Falloff Test and Effect of Temperature and Skin Factor on It: A Case Study for Coal Seams of Jamalganj Coalfield, Bangladesh

The Jamalganj coalfield with the coal deposits of the Permian Gondwana Group in the halfgraben basin was discovered in Joypurhat District, Northwestern Bangladesh, sometime in 1962. Individual coalbed thickness ranges from 0.60 m to 42 m and seams were encountered between the depth ranges of 640 m and 1158 m. Since mining has not yet begun due to the greater depth of the coal seams, several researchers have proposed a Coalbed Methane (CBM) exploration in this region. This research focuses on the permeability of the Jamalganj coal derived from the in-situ Injection Falloff Test (IFT), which is an important reservoir parameter and one of the key factors in CBM exploration and Underground Coal Gasification (UCG). In addition to that, the relationship of temperature and skin factor with permeability is one of the key findings of this research, permeability obtained from the IFT ranges from 2.57 to 121.16 mD while the skin factor ranges from - 6.11 to 50.85 and a higher temperature gradient as of about 4°C per 100 m depths was observed. The study shows that temperature has an inverse relationship with the permeability that decreases with depth and temperature increases, which is analogous with the other CBM producing reservoirs around the world. The negative skin factor denotes flow enhancement near the wellbore and a well-stimulated reservoir, and the positive skin factor indicates increased flow resistance near the wellbore, which reduces permeability. The permeability data suggest that the analyzed coal seams of the Jamalganj Coalfield are suitable for unconventional gas production by either CBM or UCG development. The Dhaka University Journal of Earth and Environmental Sciences, Centennial Special Volume June 2022: 103-112

Computational Fluid Dynamics Study on Heat Transfer Augmentation in Tube With Various V-Cut Twisted Tape

Abstract A comparative study on the enhancement mechanism by various twisted tapes (TTs) in association with the local flow pattern is investigated in detail via computational fluid dynamics (CFD). The effect of gap between TT and tube is further examined based on prior experimental validation. Result shows that the Nusselt number and friction factor will be increased by 12.6% and 18.1%, respectively, and the comprehensive thermal performance (CTP) will be improved from 1.05–1.07 to 1.12–1.15 when the gap is reduced from 1.0 mm to zero. The obtuse V-cut design in the experimental study shows superior performance, and the related enhancement mechanisms by windward, isosceles, and leeward V-cut configurations are numerically investigated. By comparing the local streamline at V-cut regions and the velocity distribution on the cross section, the windward V-cut might notably introduce fluid into the cut region and toward the tube wall subsequently, thereby resulting in effective heat transfer augmentation near the tube wall. Compared to the isosceles V-cut TT, the CTP of windward V-cut TT can be augmented from 1.15–1.18 to 1.16–1.21. However, the leeward V-cut TT will impair the CTP due to its inconsistent turbulence direction to the main helical flow. On the other hand, it is also indicated that fabricating three windward V-cuts at each helical pitch shows the highest improvement on CTP. Further increment of V-cut number may reversely reduce the CTP because of the rapid increase of flow resistance. Moreover, fabricating windward V-cuts on bilateral edges of the TT shows limited improvement.

Effects of baffles on heat and contaminant removal for tunnel smoke extraction: Evaluation criteria regarding removal efficiency and flow resistance

An exhaust vent with a baffle and its corresponding evaluation index are proposed in this paper to increase the heat and contaminant extraction efficiencies for stratified smoke flow in tunnel fires. Numerical simulation and small-scale experiments are performed to investigate the effects of baffles beneath the exhaust vent on the flow pattern and the resulting heat and contaminant extraction efficiencies. Further, we introduce performance evaluation indices PEI H and PEI M , to comprehensively evaluate the performance of the exhaust vent considering both enhancement of extraction efficiency and increase in flow resistance. The analysis results show that when a concave baffle is adopted at a height of 0.7 m from the tunnel ceiling, maximum PEI H and PEI M (1.20 and 1.17, respectively) are achieved. Compared to the vent without a baffle, the heat and contaminant extraction efficiencies are increased by 21.1% and 18.1%, respectively; however, the flow resistance is increased by only 1.8%. The proposed PEI H and PEI M can be used as a general index for evaluating the performance of other types of exhaust vents and optimizing the vent layout as well, which is helpful to improve the performance of transverse smoke exhaust system in tunnel fire.

Characterizing elastic turbulence in the three-dimensional von Kármán swirling flow using the Oldroyd-B model

We present a comprehensive three-dimensional numerical investigation of the von Kármán swirling flow between two parallel plates using the Oldroyd-B model and characterize the onset and development of elastic turbulence. We quantify the flow state with the secondary-flow strength, a measure of the average strength of the velocity fluctuations, and then define an order parameter as the time average of the secondary-flow strength. The order parameter displays a subcritical transition from the laminar to a bistable flow that switches between weakly chaotic flow and elastic turbulence. The transition to the bistable flow occurs at the critical Weissenberg number Wic=12. In the elastic turbulent state, we observe a strong increase in velocity fluctuations and flow resistance which we define as the total work performed on the fluid. Upon starting simulations in the turbulent state and subsequently lowering Wi below its critical value, we observe hysteretic behavior in the order parameter and the flow resistance, which is a common feature of a subcritical transition. Hysteresis has also been found in experiments. Additionally, we find power-law scaling in the spatial and temporal power spectra of the velocity fluctuations, a characteristic for elastic turbulence. The maximum values of the power-law exponents in our simulations are αt=3.69 for the temporal exponent and αs=3.18 for the spatial exponent, which are remarkably close to the values obtained in experiments.

Two-Phase Lung Damage Mechanisms For COVID-19 Disease, and Driving Force and Selectivity in Leukecyte Recruitment and Migration

To understand lung damages caused by COVID-19, we deduced two phases lung damage mechanisms. After the lungs are infected with COVID-19, the affected lung tissue swells and surface properties of pulmonary capillaries change, both contributing to an increased flow resistance of the capillaries. The initial damages are mainly fluid leakage in a limited number of involved alveoli. The increased vascular resistance results in retaining more white blood cells (“WBCs”) in pulmonary capillaries. Some of the WBCs may get into interstitial spaces. When more and more WBCs are dynamically retained, the vascular resistance of pulmonary capillaries further rises; and thus the overall vascular resistance of the lungs rises and pulmonary pressure rises. The rise in the pulmonary pressure in turn results in elevated capillary pressures. When pulmonary capillary pressures around the alveoli are sufficiently high, the elevated pressure causes interstitial pressures to change from normally negative values to positive values. The positive pressures cause fluid leakage to the alvoeli and thus degrade lung function. Tissue swelling, and occupation of WBCs in interstitial spaces and occupation of alvoelar spaces by leaked water result in reduced deformable and compressible spaces, and thus causes a further rise of the vascular resistance of the lungs. When the pulmonary pressure has reached a critical point as in the second phase, the blood breaks capillary walls and squeezes through interstitial spaces to reach alveolar spaces, resulting in irreversible lung damages. Among potential influencing factors, the available space in the thorax cage, temperature, and humid are expected to have great impacts. The free space in the thorax cage, lung usable capacity, and other organ usable capacities are the major factors that determine the arrival time of last- phase irreversible damage. The mechanisms imply that the top priority for protecting lungs is maintaining pulmonary micro-circulation and preserving organ functions in the entire disease course while controlling viral reproduction should be stressed in the earliest time possible. The mechanisms also explain how leukecytes are “recruited and migrated” into inflamed tissues by dynamic retention.

Influence of Colloidal Silica and Silica Fume on the Rheology and Mechanical Properties of High-Performance Shotcrete

Wet-mix shotcrete technology is widely used worldwide. One possible approach to developing this technology is the use of various admixtures to improve its pumpability and shootability. However, the commonly used slump test is insufficient to evaluate the capabilities of fresh shotcrete, because they do not assess the influence of multiple factors on its workability and placement properties. The purpose of this paper is to investigate and determine the effect of colloidal silica with a particle size of 10 nm (with 2%, 3%, and 4% replacement) on high-performance shotcrete from a rheological and mechanical perspective. In addition, to compare these indicators, silica fume is used with a 4% and 7% mixing rate. Rheological parameters such as torque viscosity and flow resistance of fresh concrete were evaluated using an ICAR rheometer, and the measured values were correlated with pumpability and shootability. The results of the study demonstrate that torque viscosity decreases and flow resistance increases with the increase in the percentage of colloidal silica. Similarly, silica fume increases flow resistance and reduces torque viscosity. These findings demonstrate the possibility of obtaining high-performance shotcrete with improved pumpability and shootability by employing more colloidal silica and silica fume. It was also established that these supplementary cement materials improve the strength and durability of shotcrete.

Effect of water behaviour on the oil transport in illite nanopores: Insights from a molecular dynamics study

Both of the water injected by hydraulic fracturing and connate water can coexist with oil inside shale nanopores, which leads to different oil/water mixtures behaviors compared to those with only oil and water. However, the behavior of oil/water mixtures in shale pores is not clear. In this study, molecular dynamics simulations are applied to investigate the behavior of oil/water mixtures in a rectangle-shaped illite nanopore that has a hydrophilic internal surface. The behavior of oil/water mixtures is much more complex than that of gas/water mixture because water can dissolve in oil at particular pressures and temperatures. The results of the calculation show that illite confinement has a strong impact on both oil distribution and water adsorption. A structure of the oil droplets being wrapped by water can be identified in the illite nanopores which may be contributed by the hydrophilic surfaces. Besides, water bridges may form due to the dispersion of unabsorbed water into the oil and could increase flow resistance in the pore. Viscosity and self-diffusion of fluids are increased and decreased separately in the nanopores compared to those without the illite confinement which is called the bulk condition in this study. In the bulk condition, increasing the water concentration will reduce the shear viscosity of fluids, while an opposite phenomenon appears in shale nanopores. The results of this study may help to deepen the understandings of water and oil transport in nanopores and provide some guidelines in extracting shale resources in an aqueous environment.