Computational fluid dynamics (CFD) based approach to consequence assessment of accidental release of hydrocarbon during storage and transportation

This thesis investigated the risk of accidental release of hydrocarbons during transportation and storage. Transportation of hydrocarbons from an offshore platform to processing units through subsea pipelines involves risk of release due to pipeline leakage resulting from corrosion, plastic deformation caused by seabed shakedown or damaged by contact with drifting iceberg. The environmental impacts of hydrocarbon dispersion can be severe. Overall safety and economic concerns of pipeline leakage at subsea environment are immense. A large leak can be detected by employing conventional technology such as, radar, intelligent pigging or chemical tracer but in a remote location like subsea or arctic, a small chronic leak may be undetected for a period of time. In case of storage, an accidental release of hydrocarbon from the storage tank could lead pool fire; further it could escalate to domino effects. This chain of accidents may lead to extremely severe consequences. Analyzing past accident scenarios it is observed that more than half of the industrial domino accidents involved fire as a primary event, and some other factors for instance, wind speed and direction, fuel type and engulfment of the compound. In this thesis, a computational fluid dynamics (CFD) approach is taken to model the subsea pipeline leak and the pool fire from a storage tank. A commercial software package ANSYS FLUENT Workbench 15 is used to model the subsea pipeline leakage. The CFD simulation results of four different types of fluids showed that the static pressure and pressure gradient along the axial length of the pipeline have a sharp signature variation near the leak orifice at steady state condition. Transient simulation is performed to obtain the acoustic signature of the pipe near leak orifice. The power spectral density (PSD) of acoustic signal is strong near the leak orifice and it dissipates as the distance and orientation from the leak orifice increase. The high-pressure fluid flow generates more noise than the low-pressure fluid flow. In order to model the pool fire from the storage tank, ANSYS CFX Workbench 14 is used. The CFD results show that the wind speed has significant contribution on the behavior of pool fire and its domino effects. The radiation contours are also obtained from CFD post processing, which can be applied for risk analysis. The outcome of this study will be helpful for better understanding of the domino effects of pool fire in complex geometrical settings of process industries. The attempt to reduce and prevent risks is discussed based on the results obtained from the numerical simulations of the numerical models.

Structure and dynamics of thin colloidal films

This thesis begins by studying the thickness of evaporative spin coated colloidal crystals and demonstrates the variation of the thickness as a function of suspension concentration and spin rate. Particularly, the films are thicker with higher suspension concentration and lower spin rate. This study also provides evidence for the reproducibility of spin coating in terms of the thickness of the resulting colloidal films. These colloidal films, as well as the ones obtained from various other methods such as convective assembly and dip coating, usually possess a crystalline structure. Due to the lack of a comprehensive method for characterization of order in colloidal structures, a procedure is developed for such a characterization in terms of local and longer range translational and orientational order. Translational measures turn out to be adequate for characterizing small deviations from perfect order, while orientational measures are more informative for polycrystalline and highly disordered crystals. Finally, to obtain an understanding of the relationship between dynamics and structure, the dynamics of colloids in a quasi-2D suspension as a function of packing fraction is studied. The tools that are used are mean square displacement (MSD) and the self part of the van Hove function. The slow down of dynamics is observed as the packing fraction increases, accompanied with the emergence of 6-fold symmetry within the system. The dynamics turns out to be non-Gaussian at early times and Gaussian at later times for packing fractions below 0.6. Above this packing fraction, the dynamics is non-Gaussian at all times. Also the diffusion coefficient is calculated from MSD and the van Hove function. It goes down as the packing fraction is increased.