Sand Production Modelled with Darcy Fluid Flow Using Discrete Element Method

In the process of recovering oil in weak sandstone formations, the strength of sandstones around the wellbore is weakened due to the increase of effective stress/load from the completion activities around the cavity. The weakened and de-bonded sandstone may be eroded away by the produced fluid, which is termed sand production. It is one of the major trending subjects in the petroleum industry because of its significant negative impacts, as well as some observed positive impacts. For efficient sand management therefore, there has been need for a reliable study tool to understand the mechanism of sanding. One method of studying sand production is the use of the widely recognized Discrete Element Method (DEM), Particle Flow Code (PFC3D) which represents sands as granular individual elements bonded together at contact points. However, there is limited knowledge of the particle-scale behavior of the weak sandstone, and the parameters that affect sanding. This paper aims to investigate the reliability of using PFC3D and a simple Darcy flow in understanding the sand production behavior of a weak sandstone. An isotropic tri-axial test on a weak oil sandstone sample was first simulated at a confining stress of 1MPa to calibrate and validate the parallel bond models of PFC3D using a 10m height and 10m diameter solid cylindrical model. The effect of the confining stress on the number of bonds failure was studied using this cylindrical model. With the calibrated data and sample material properties obtained from the tri-axial test, simulations without and with fluid flow were carried out to check on the effect of Darcy flow on bonds failure using the same model geometry. The fluid flow network comprised of every four particles connected with tetrahedral flow pipes with a central pore or flow domain. Parametric studies included the effects of confining stress, and fluid pressure; as well as validating flow rate – permeability relationship to verify Darcy’s fluid flow law. The effect of model size scaling on sanding was also investigated using 4m height, 2m diameter model. The parallel bond model successfully calibrated the sample’s strength of 4.4MPa, showing a sharp peak strength before strain-softening, similar to the behavior of real cemented sandstones. There seems to be an exponential increasing relationship for the bigger model, but a curvilinear shape for the smaller model. The presence of the Darcy flow induced tensile forces and increased the number of broken bonds. For the parametric studies, flow rate has a linear relationship with permeability at constant pressure head. The higher the fluid flow pressure, the higher the number of broken bonds/sanding. The DEM PFC3D is a promising tool to studying the micromechanical behavior of cemented sandstones.

Streamwise Vorticity in the Wake of a Sliding Bubble

In many practical situations, bubbles are dispersed in a liquid phase. Understanding these complex bubbly flows is therefore a key issue for applications such as shell and tube heat exchangers, mineral flotation and oxidation in water treatment. Although a large body of work exists for bubbles rising in an unbounded medium, that of bubbles rising in constricted geometries has received less attention. The particular case of a bubble sliding underneath an inclined surface is common to two-phase flow systems. The current study intends to expand this knowledge by performing experiments to quantify the streamwise flow structures associated with a single sliding air bubble under an inclined surface in quiescent water. This is achieved by means of two-dimensional, two-component particle image velocimetry (PIV), performed with a continuous wave laser and high-speed camera. PIV vorticity fields obtained in a plane perpendicular to the sliding surface show that there is significant bulk fluid motion away from the surface. The associated momentum of the bubble means that this wake motion persists for a significant time before viscous dissipation. The magnitude and direction of the flow structures in the streamwise measurement plane are found to depend on the point on its path through which the bubble enters the plane. This entry point, represented by a phase angle, affects the nature and strength of the vortical structures. This study reconstructs the vorticity field in the wake of the bubble, converting the field at different instances in time to slices of a large-scale wake structure. This is, in essence, Taylor’s ”frozen turbulence” hypothesis. Applying this to the vorticity fields provides a pseudo three-dimensional representation from 2-D data, allowing for a more intuitive understanding of the bubble wake. This study provides insights into the complex dynamics of a situation common to many engineering applications, particularly shell and tube heat exchangers in the nucleate boiling regime.

Polydopamine Nanoparticle as a Stable and Capacious Nano-Reservoir of Rifampicin

Application of nanoscience in biomedical field has come across as a new era. This study involves the synthesis of nano drug carrier with antibiotic loading. Based on the founding that polydopamine (PDA) nanoparticles could be formed via self-polymerization of dopamine at alkaline pH, one-step synthesis of rifampicin coupled polydopamine (PDA-R) nanoparticles was achieved by adding rifampicin into the dopamine solution. The successful yield of PDA nanoparticles with or without the presence of rifampicin during the polymerization process was characterized by scanning electron microscopy, Fourier transform infrared spectroscopy, and Raman spectroscopy. Drug loading was monitored by UV-vis spectroscopy and the loading efficiency of rifampicin was calculated to be 76%. Such highly capacious nano-reservoir was found very stable with little drug leakage at pH 3.

Modeling Bessel Beams and Their Discrete Superpositions from the Generalized Lorenz-Mie Theory to Calculate Optical Forces over Spherical Dielectric Particles

In this work, we propose an algorithm developed under Python language for the modeling of ordinary scalar Bessel beams and their discrete superpositions and subsequent calculation of optical forces exerted over dielectric spherical particles. The mathematical formalism, based on the generalized Lorenz-Mie theory, is implemented in Python for its large number of free mathematical (as SciPy and NumPy), data visualization (Matplotlib and PyJamas) and multiprocessing libraries. We also propose an approach, provided by a synchronized Software as Service (SaaS) in cloud computing, to develop a user interface embedded on a mobile application, thus providing users with the necessary means to easily introduce desired unknowns and parameters and see the graphical outcomes of the simulations right at their mobile devices. Initially proposed as a free Android-based application, such an App enables data post-processing in cloud-based architectures and visualization of results, figures and numerical tables.

Transformations of Spatial Distributions of Bio-Polymers and Nanoparticles in Water Suspensions Induced by Resonance-Like Low Frequency Electrical Fields

Water suspensions of in-organic (metals and oxides) and organic nano-objects (chitozan and collagen) were subjected to the treatment of direct and alternative electrical fields. In addition to quasi-periodical spatial patterning resonance-like performance of spatial distributions of these suspensions has been found at low frequencies of alternating electrical field. These resonances are explained as the result of creation of equilibrium states of groups of charged nano-objects with opposite signs of charges at the interparticle distances where the forces of Coulomb attraction are compensated by the repulsion forces induced by relatively negative polarization of hydrated regions surrounding the nanoparticles with respect to pure water. The low frequencies of these resonances are explained by comparatively big distances between the particles and their big masses with t\respect to masses of atoms constituting molecules with high resonance frequencies. These new resonances open a new approach to detailed modeling and understanding of mechanisms of the influence of electrical fields on the functioning of internal organs of living organisms at the level of cells and neurons.

Effects of Upstream Wall Roughness on Separated Turbulent Flow over a Forward Facing Step in an Open Channel

The effect of upstream surface roughness over a smooth forward facing step in an open channel was investigated using a particle image velocimetry technique. Three different upstream surface topographies consisting of hydraulically smooth wall, sandpaper 36 grit and sand grains were examined. Besides the wall roughness conditions, all other upstream flow characteristics were kept constant. It was also observed that upstream roughness decreased the approach velocity by 2% and 10% but increased the turbulence intensity by 14% and 35% at the wall-normal distance corresponding to the top plane of the step compared to smooth upstream. The results showed that roughness decreased the reattachment lengths by 14% and 30% compared to smooth upstream. Although the magnitudes of maximum positive and negative Reynolds shear stress in separated and reattached region were 0.02Ue for all the cases, the physical size of both the maximum and minimum contour levels were decreased by increasing upstream roughness.

Energy Saving of the Paint with Mineral Insulators: Simulation and Study on Different Climates

By using an adequate thermal barrier coating in buildings the energy saving will be happened. In this study, a range of wall paints with different absorption coefficient in different climates has been investigated. In order to study these effects, heating and cooling loads of a common building with different ordinary paints and paint with mineral coating have been calculated. The effect of building paint in different climatic condition was studied and comparison was done between ordinary paints and paint with mineral insulators in temperate climate to obtain optimized energy consumption. The results have been shown that coatings with inorganic micro particles as insulation reduce the energy consumption of buildings around 14%.

Investigation of Utilization Possibility of Fluid Gas Desulfurization Waste for Industrial Waste Water Treatment

Flue gas desulfurization gypsum (FGD) is a waste material arouse from coal power plants. Hydroxyapatite (HAP) is a biomaterial with porous structure. In this study, FGD gypsum which retrieved from coal power plant in Turkey was characterized and HAP particles which can be used as an adsorbent in wastewater treatment application were synthesized from the FGD gypsum. The raw materials are characterized by using X Ray Diffraction (XRD) and Fourier transform infrared spectroscopy (FT-IR) techniques and produced HAP are characterized by using XRD. As a result, HAP particles were synthesized at the molar ratio of 5:10, 5:15, 5:20, 5:24, at room temperature, in alkaline medium (pH=11) and in 1 hour-reaction time. Among these conditions, 5:20 had the best result.

Evaluation of PTFE Composites with Mineral Tailing Considering Friction, Wear and Cost

The tribological test with Pin-On-Disc configuration measures friction and wear properties in dry or lubricated sliding surfaces of a variety of materials and coatings. Polymeric matrix composites loaded with mineral filler were used, 1%, 3%, 10%, 30%, and 50% mass percentage of filler, to reduce the material cost by using mineral tailings. Using a pin-on-disc tribometer to quantify coefficient of friction and wear resistance of the specimens. The parameters known to performing the test were 300 rpm rotation, normal load of 16N and duration of 33.5 minutes. The composite with 10% mineral filler performed better, considering that the wear resistance was good when compared to the other compositions and an average low coefficient of friction, in the order of μ ≤ 0.15.

Optical and Double Folding Model Analysis for Alpha Particles Elastically Scattered from 9Be and 11B Nuclei at Different Energies

Elastic scattering of α-particles from 9Be and 11B nuclei at different alpha energies have been analyzed. Optical model parameters (OMPs) of α-particles elastic scattering by these nuclei at different energies have been obtained. In the present calculations, the real part of the optical potential are derived by folding of nucleonnucleon (NN) interaction into nuclear matter density distribution of the projectile and target nuclei using computer code FRESCO. A density-dependent version of the M3Y interaction (CDM3Y6), which is based on the G-matrix elements of the Paris NN potential, has been used. Volumetric integrals of the real and imaginary potential depth (JR, JW) have been calculated and found to be energy dependent. Good agreement between the experimental data and the theoretical predictions in the whole angular range. In double folding (DF) calculations, the obtained normalization coefficient Nr is in the range 0.70–1.32.

Localized and Time-Resolved Velocity Measurements of Pulsatile Flow in a Rectangular Channel

The exploitation of flow pulsation in micro- and mini-channels is a potentially useful technique for enhancing cooling of high-end photonics and electronics systems. It is thought that pulsation alters the thickness of the hydrodynamic and thermal boundary layers, and hence affects the overall thermal resistance of the heat sink. Although the fluid mechanics and heat transfer are inextricably linked, it can be useful to decouple the parameters to better understand the mechanisms underlying any heat transfer enhancement. Using two-dimensional, two-component particle image velocimetry, the current work intends to characterize the heat transfer mechanisms in pulsating flow with a mean Reynolds number of 48 by experimentally quantifying the hydrodynamics of a generic liquid-cooled channel geometry. Flows circulated through the test section by a gear pump are modulated using a controller to achieve sinusoidal flow pulsations with Womersley numbers of 7.45 and 2.36 and an amplitude ratio of 0.75. It is found that the transient characteristics of the measured velocity profiles are dependent on the speed of oscillation, in accordance with the analytical solution for flow in a rectangular channel. A large velocity overshoot is observed close to the wall at high frequencies, resulting from the interaction of near-wall viscous stresses and inertial effects of the main fluid body. The steep velocity gradients at the wall are indicative of augmented heat transfer, although the local flow reversal may reduce the upstream temperature difference in heat transfer applications. While unsteady effects remain evident at the lower frequency, the annular effect subsides and retreats from the wall. The shear rate at the wall is increased during the accelerating half-cycle and decreased during deceleration compared to steady flow, suggesting that the flow may experience both enhanced and diminished heat transfer during a single period. Hence, the thickness of the hydrodynamic boundary layer is reduced for positively moving flow during one half of the pulsation cycle at the investigated frequencies. It is expected that the size of the thermal boundary layer is similarly reduced during the cycle, leading to intervals of heat transfer enhancement.

Surface Characteristics of Bacillus megaterium and Its Adsorption Behavior onto Dolomite

Surface characteristics of Bacillus megaterium strain were investigated; zeta potential, FTIR and contact angle were measured. Surface energy components including Lifshitz-van der Waals, Hamaker constant, and acid/base components (Lewis acid/Lewis base) were calculated from the contact angle data. The results showed that the microbial cells were negatively charged over all pH regions with high values at alkaline region. A hydrophilic nature for the strain was confirmed by contact angle and free energy of adhesion between microbial cells. Adsorption affinity of the strain toward dolomite was studied at different pH values. The results showed that the cells had a high affinity to dolomite at acid pH comparing to neutral and alkaline pH. Extended DLVO theory was applied to calculate interaction energy between B. megaterium cells and dolomite particles. The adsorption results were in agreement with the results of Extended DLVO approach. Surface changes occurred on dolomite surface after the bio-treatment were monitored; contact angle decreased from 69° to 38° and the mineral’s floatability decreased from 95% to 25% after the treatment.

Metallic Coating for Carbon Fiber Reinforced Polymer Matrix Composite Substrate

This paper investigates the application of metallic coatings on high fiber volume fraction carbon/epoxy polymer matrix composites. For the grip of the metallic layer, a method of modifying the surface of the composite by introducing a mixture of copper and steel powder (filler powders) which can reduce the impact of thermal spray particles. The powder was introduced to the surface at the time of the forming. Arc spray was used to project the zinc coating layer. The substrate was grit blasted to avoid poor adherence. The porosity, microstructure, and morphology of layers are characterized by optical microscopy, SEM and image analysis. The samples were studied also in terms of hardness and erosion resistance. This investigation did not reveal any visible evidence damage to the substrates. The hardness of zinc layer was about 25.94 MPa and the porosity was around (∼6.70%). The erosion test showed that the zinc coating improves the resistance to erosion. Based on the results obtained, we can conclude that thermal spraying allows the production of protective coating on PMC. Zinc coating has been identified as a compatible material with the substrate. The filler powders layer protects the substrate from the impact of hot particles and allows avoiding the rupture of brittle carbon fibers.

Enhancement Effect of Superparamagnetic Iron Oxide Nanoparticle-Based MRI Contrast Agent at Different Concentrations and Magnetic Field Strengths

Magnetic Resonance Imaging Contrast Agents (MRI-CM) are significant in the clinical and biological imaging as they have the ability to alter the normal tissue contrast, thereby affecting the signal intensity to enhance the visibility and detectability of images. Superparamagnetic Iron Oxide (SPIO) nanoparticles, coated with dextran or carboxydextran are currently available for clinical MR imaging of the liver. Most SPIO contrast agents are T2 shortening agents and Resovist (Ferucarbotran) is one of a clinically tested, organ-specific, SPIO agent which has a low molecular carboxydextran coating. The enhancement effect of Resovist depends on its relaxivity which in turn depends on factors like magnetic field strength, concentrations, nanoparticle properties, pH and temperature. Therefore, this study was conducted to investigate the impact of field strength and different contrast concentrations on enhancement effects of Resovist. The study explored the MRI signal intensity of Resovist in the physiological range of plasma from T2-weighted spin echo sequence at three magnetic field strengths: 0.47 T (r1=15, r2=101), 1.5 T (r1=7.4, r2=95), and 3 T (r1=3.3, r2=160) and the range of contrast concentrations by a mathematical simulation. Relaxivities of r1 and r2 (L mmol-1 Sec-1) were obtained from a previous study and the selected concentrations were 0.05, 0.06, 0.07, 0.08, 0.09, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0, 2.0, and 3.0 mmol/L. T2-weighted images were simulated using TR/TE ratio as 2000 ms /100 ms. According to the reference literature, with increasing magnetic field strengths, the r1 relaxivity tends to decrease while the r2 did not show any systematic relationship with the selected field strengths. In parallel, this study results revealed that the signal intensity of Resovist at lower concentrations tends to increase than the higher concentrations. The highest reported signal intensity was observed in the low field strength of 0.47 T. The maximum signal intensities for 0.47 T, 1.5 T and 3 T were found at the concentration levels of 0.05, 0.06 and 0.05 mmol/L, respectively. Furthermore, it was revealed that, the concentrations higher than the above, the signal intensity was decreased exponentially. An inverse relationship can be found between the field strength and T2 relaxation time, whereas, the field strength was increased, T2 relaxation time was decreased accordingly. However, resulted T2 relaxation time was not significantly different between 0.47 T and 1.5 T in this study. Moreover, a linear correlation of transverse relaxation rates (1/T2, s–1) with the concentrations of Resovist can be observed. According to these results, it can conclude that the concentration of SPIO nanoparticle contrast agents and the field strengths of MRI are two important parameters which can affect the signal intensity of T2-weighted SE sequence. Therefore, when MR imaging those two parameters should be considered prudently.

Investigating the Shear Behaviour of Fouled Ballast Using Discrete Element Modelling

For several hundred years, the design of railway tracks has practically remained unchanged. Traditionally, rail tracks are placed on a ballast layer due to several reasons, including economy, rapid drainage, and high load bearing capacity. The primary function of ballast is to distributing dynamic track loads to sub-ballast and subgrade layers, while also providing lateral resistance and allowing for rapid drainage. Upon repeated trainloads, the ballast becomes fouled due to ballast degradation and the intrusion of fines which adversely affects the strength and deformation behaviour of ballast. This paper presents the use of three-dimensional discrete element method (DEM) in studying the shear behaviour of the fouled ballast subjected to direct shear loading. Irregularly shaped particles of ballast were modelled by grouping many spherical balls together in appropriate sizes to simulate representative ballast aggregates. Fouled ballast was modelled by injecting a specified number of miniature spherical particles into the void spaces. The DEM simulation highlights that the peak shear stress of the ballast assembly decreases and the dilation of fouled ballast increases with an increase level of fouling. Additionally, the distributions of contact force chain and particle displacement vectors were captured during shearing progress, explaining the formation of shear band and the evolutions of volumetric change of fouled ballast.

Optical Flow Technique for Supersonic Jet Measurements

This paper outlines the development of an experimental technique in quantifying supersonic jet flows, in an attempt to avoid seeding particle problems frequently associated with particle-image velocimetry (PIV) techniques at high Mach numbers. Based on optical flow algorithms, the idea behind the technique involves using high speed cameras to capture Schlieren images of the supersonic jet shear layers, before they are subjected to an adapted optical flow algorithm based on the Horn-Schnuck method to determine the associated flow fields. The proposed method is capable of offering full-field unsteady flow information with potentially higher accuracy and resolution than existing point-measurements or PIV techniques. Preliminary study via numerical simulations of a circular de Laval jet nozzle successfully reveals flow and shock structures typically associated with supersonic jet flows, which serve as useful data for subsequent validation of the optical flow based experimental results. For experimental technique, a Z-type Schlieren setup is proposed with supersonic jet operated in cold mode, stagnation pressure of 4 bar and exit Mach of 1.5. High-speed singleframe or double-frame cameras are used to capture successive Schlieren images. As implementation of optical flow technique to supersonic flows remains rare, the current focus revolves around methodology validation through synthetic images. The results of validation test offers valuable insight into how the optical flow algorithm can be further improved to improve robustness and accuracy. Despite these challenges however, this supersonic flow measurement technique may potentially offer a simpler way to identify and quantify the fine spatial structures within the shock shear layer.

Anticorrosive Polyurethane Clear Coat with Self-Cleaning Character

We have aimed to produce a self-cleaning transparent polymer coating with polyurethane (PU) matrix as the latter is highly solvent, chemical and weather resistant having good mechanical properties. Nano-silica modified by 1H, 1H, 2H, 2Hperflurooctyltriethoxysilane was incorporated into the PU matrix for attaining self-cleaning ability through hydrophobicity. The modification was confirmed by particle size analysis and scanning electron microscopy (SEM). Thermo-gravimetric (TGA) studies were carried to ascertain the grafting of silane onto the silica. Several coating formulations were prepared by varying the silica loading content and compared to a commercial equivalent. The effect of dispersion and the morphology of the coated films were assessed by SEM analysis. All coating standardized tests like solvent resistance, adhesion, flexibility, acid, alkali, gloss etc. have been performed as per ASTM standards. Water contact angle studies were conducted to analyze the hydrophobic character of the coating. In addition, the coatings were also subjected to salt spray and accelerated weather testing to analyze the durability of the coating.

Parametric Studies of Wood Pyrolysis Particles

In the present study, a numerical approach to describe the pyrolysis of a single solid particle of wood is used to study the influence of various conditions such as particle size, heat transfer coefficient, reactor temperature and heating rate. The influence of these parameters in the change of the duration of the pyrolysis cycle was studied. Mathematical modeling was employed to simulate the heat, mass transfer, and kinetic processes inside the reactor. The evolutions of the mass loss as well as the evolution of temperature inside the thick piece are investigated numerically. The elaborated model was also employed to study the effect of the reactor temperature and the rate of heating on the change of the temperature and the local loss of the mass inside the piece of wood. The obtained results are in good agreement with the experimental data available in the literature.

Aerodynamic Study of Vehicle Wind Tunnel and Water Tunnel for Analysis of Bodies

The simulation in wind tunnel is used thoroughly to model real situations of drainages of air. Besides the automotive industry, a great number of applications can be numbered: dispersion of pollutant, studies of pedestrians’ comfort, and dispersion of particles. This work had the objective of visualizing the characteristics aerodynamics of two automobiles in different ways. To accomplish that drainage of air a fan that generated a speed exists (measured with anemometer of hot thread) of 4,1m/s and 4,95m/s. To visualize the path of the air through the cars, in the wind tunnel, smoke was used, obtained with it burns of vegetable oil. For “to do smoke” vegetable oil was used, that was burned for a tension of 20V generated by a thread of 2,5mm. The cars were placed inside of the wind tunnel with the drainage of “air-smoke” and photographed, registering like this the path lines around them, in the 3 different speeds.

Experimental Investigation of Proton Exchange Membrane Fuel Cells Operated with Nanofiber and Nanofiber/Nanoparticle

Nanofibers are defined as fibers with diameters less than 100 nanometers. In this study, behaviours of activated carbon nanofiber (ACNF), carbon nanofiber (CNF), polyacrylonitrile/ carbon nanotube (PAN/CNT), polyvinyl alcohol/nanosilver (PVA/Ag) in proton exchange membrane (PEM) fuel cells are investigated experimentally. This material was used as gas diffusion layer (GDL) in PEM fuel cells. In this study, the electrical conductivities of nanofiber and nanofiber/nanoparticles have been studied to understand their effects on PEM fuel cell performance. According to the experimental results, the maximum electrical conductivity performance of the fuel cell with nanofiber was found to be at PVA/Ag (at UConn condition). The electrical conductivities of CNF, ACNF, PAN/CNT are lower for PEM. The resistance of cell with PVA/Ag is lower than the resistance of cell with PAN/CNT, ACNF, CNF.