Search

The dispersion (transport and breakup) of oil droplets under breaking waves is essential for evaluating the impact of oil spills on the environment, and for designing countermeasures. In this part of the work (Part I), the hydrodynamics of a wave breaker in a 1.0 m deep tank obtained experimentally by Rapp and Melville (Philos Trans R Soc Lond A Math Phys Eng Sci 331(1622):735–800, 1990) using the dispersive focusing method was investigated numerically using Reynold-averaged Navier–Stokes (RANS) within the CFD code ANSYS Fluent. The renormalization group (RNG) k–e turbulence closure model was adopted to simulate wave turbulence, and the transient water–air interface was captured using volume of fluid (VOF) method. The simulated surface excursion and velocity fields matched closely the experimental observations during wave breaking. The energy dissipation of the wave crest showed good agreement with recent laboratory work and numerical simulations through evaluating the breaking parameter. The RANS approach was able to reproduce the turbulent kinetic energy field engendered by the breakers and simulated the residual turbulence within about two wave periods after the passage of the wave train. The VOF scheme Compressive provided a better agreement with the observation than the Geo-reconstruct scheme. The approach herein suggests that RANS method coupled with VOF in ANSYS Fluent is capable of capturing the major hydrodynamic forces and turbulence, and thus could be used to predict environmental processes within the breaking waves such as oil droplet formation and transport.

Response surface methodology is used in this study to model and optimize a solar thermal-photovoltaic vacuum membrane distillation (STPVMD) system with PVDF hollow fiber membrane. Regression models have been developed to forecast the effect of various operation factors on the permeate flux and the energy consumption. The operation factors in the model include the feed inlet temperature, the feed flow rate and the vacuum pressure. Analysis of variance is used to statistically validate the regression models. With permeate flux as the objective function of optimization, the optimal operation parameters are found to be 63 °C for the feed inlet temperature, 237 L / h for the feed flow rate, and 750 kPa for the vacuum pressure. Experimental tests show that the resulting permeate flux is 6.26 L / m 2 · h , which is slightly higher than the predicted value of 6.05 L / m 2 · h . The corresponding energy consumption is predicted to be 12.52 L / kW · h , demonstrating the effectiveness and reliability of the models.

The transport of virus-laden particles was investigated numerically in an archetypical supermarket configuration of area 1,200 m2 and ceiling height of 4.5 m. The particles were tracked u...

With the increase in transport of oil by rail, the probability of oil spills in rivers has increased. Traditionally, focus has been placed on oil slicks moving on the water surface. However...

Understanding the shear breakup in jet flows and the formation of droplets from ligaments is important to determine the final droplet size distribution (DSD). The initial droplet size, which affects the final DSD, is considered to be generated by the shear breakup. Large eddy simulation (LES) was performed to investigate the shear breakup in liquid-liquid jet flows. The explicit Volume of Fluid (VOF) model with the geometric reconstruction scheme was used to capture the oil-water interface. The estimated oil distribution including wave peaks, ligaments, droplets and water streaks were compared to the experiments with a good agreement. The estimated DSD matched with the measurements favorably well. In the simulation, the formation of droplets with a smooth and curved surface from ligaments or sheet-like structures was obtained. Different mechanisms were observed along with the shear layer including the formation of droplets from ligament through the capillary forces, breakage of a droplet into smaller ones and attachment of a droplet to a ligament. The destructive shear forces and resisting surface tension forces were quantified on stretching and retracting ligaments. The influence of internal viscous force was found to be negligible due to low oil viscosity. The critical capillary number was found to be larger than 5.0 for ligaments breaking with the shear breakup. The capillary number was below unity for retracting ligaments. The coalescence of two equal-sized droplets was obtained in the shear breakup region. The shear stress magnitude at the contact region increased more than two folds. The total surface area decreased nearly 20% after the coalescence.

A numerical study was undertaken to investigate subsurface release and fate of benzene and toluene in a tidally influenced beach. The simulations were conducted by using a numerical model BIOMARUN, which coupled a multi-Monod kinetic model BIOB to a density-dependent variably saturated groundwater flow model MARUN. The fate and transport of the contaminant plume were characterized by computing its centroid trajectory, spreading area and percentage of biodegradation in beach saturated and unsaturated zone, respectively. Key factors likely affecting this process were investigated, including tide amplitude, capillarity and hydraulic conductivity. It was found that aerobic biodegradation was the major fate of the contaminant plume in the beach. Tidal action twisted the centroid of the contaminant plume during its migration in the beach, which increased the residence time of the plume in the beach. High tidal range significantly altered the spatial distribution of the contaminant biodegradation in the beach. In contrast, the capillary fringe had impacts on the percentage of benzene biodegraded in the saturated and unsaturated zone of the beach. The increase in capillary fringe enhanced the percentage of the contaminant biodegraded in the unsaturated zone, up to 40%, which is comparable to that in the saturated zone. Hydraulic conductivity seemed to have large impacts on the biodegradation rate of the contaminant in the beach. Higher hydraulic conductivity induced faster contaminant biodegradation in the beach.

Shape memory polymers (SMPs) are soft active materials that have an ability to retain a temporary shape, and revert back to their original shape when triggered by a suitable stimulus, typically an increase in temperature. These materials are finding wide use in a variety of fields such as biomedical and aerospace engineering; hence it is important to model their mechanical behavior. Crystallizable shape memory polymers (CSMPs) is an important subclass of SMPs, and their temporary shape is fixed by a crystalline phase, while return to the original shape is due to the melting of this crystalline phase. In our earlier work, we have studied the mechanical behavior of CSMPs within a mechanical setting by considering the original amorphous network above the recovery temperature as a hyperelastic material. In this article, we extend our earlier work to incorporate the temperature-dependent viscoelasticity into the developed constitutive model to study the mechanical behavior of CSMPs. The viscoelastic behavior of the polymers at high temperature is simulated through a rate type model. Furthermore, the model of the semi-crystalline polymer after the onset of crystallization is developed based on the mixture theory and the theory of “multiple natural configurations”. In addition, we have applied the model to a specific boundary value problem, namely uniaxial extension. The shape memory cycles of the CSMPs under different stretch rates have been studied. The results are consistent with what has been observed in experiments.

Summary Shape-memory polymers (SMPs) belong to a class of smart materials that have shown promise for a wide range of applications. They are characterized by their ability to maintain a temporary deformed shape and return to an original parent permanent shape. In this paper, we consider the coupled photomechanical behavior of light activated shape-memory polymers (LASMPs), focusing on the numerical aspects for finite element simulations at the engineering scale. The photomechanical continuum framework is summarized, and some specific constitutive equations for LASMPs are described. Numerical implementation of the multiphysics governing partial differential equations takes the form of a user defined element subroutine within the commercial software package ABAQUS. We verify our two-dimensional and three-dimensional finite element procedure for multiple analytically tractable cases. To show the robustness of the numerical implementation, simulations are performed under various geometries and complex photomechanical loading. Copyright © 2016 John Wiley & Sons, Ltd.

Experimental results were presented for the release of diesel oil from a one-inch (2.5 cm) vertical pipe in a crossflow at 0.27 m/s. The ratio of jet velocity to crossflow speed was 5.0 and the Reynolds number based on jet velocity and pipe diameter was 7.1 × 10 3 . In the experiments, the plume shape was photographed, and the oil droplets were measured at two vertical locations on the center axis of the plume. Acoustic Doppler velocimetry (ADV) data was also obtained and compared to numerical predictions. The plume was simulated using large eddy simulation (LES), and the mixture multiphase model. The impact of the oil buoyancy was captured by adding a transport term to the volume fraction equation. Using the rise velocity based on d50 (volume-median) droplet size in the lower part of the plume allowed us to capture the lower boundary of the plume, but the estimated upper boundary of the plume penetrated less into the crossflow as compared to the experimental findings. However, using the rise velocity of the d50 at the upper part of the plume allowed one to estimate the upper boundary of the plume. As the droplets are too small to be resolved by the LES, we could not use a systematic approach to allow the multiphase plume to spread to mimic the observations. Based on the simulation results, the interaction between the jet and crossflow yielded small-sized flow structures near the upper boundary of the plume. The wake vortices initiated from the leeward side of the plume showed an alternating vorticity pattern in the wake. The shear layer vortices were induced by Kevin-Helmholtz instabilities mostly on the windward side of the plume. The formation of counter rotating vortex pair (CVP) altered greatly the hydrodynamics of the jet from that of a vertical jet to manifest flow reversals in all directions. The formation of CVP is likely to enhance the mixing of chemicals and droplets within the plume.

Canada's production, transport, and sale of diluted bitumen (dilbit) products are expected to increase by a million barrels per day over the next decade. The anticipated growth in oil production and transport increases the risk of oil spills in aquatic areas and places greater demands on oil spill capabilities to respond to spills, which have raised stakeholder concerns. Current oil spill models only predict the transport of bitumen blends that are used in contingency plans and oil spill response strategies, rather than changes in the oil's physical properties that are relevant to spill response. We conducted weathering studies of five oil products (two conventional oils and three bitumen blends) in the Department of Fisheries and Oceans' flume tank. We also considered two initial oil slick thicknesses, 4.0 mm and 7.0 mm. We found that there is a major difference in the time evolution of oil properties (density and viscosity), raising doubts on weathering models that do not consider the thickness of oil. We also developed empirical expressions for the evolution of the density and viscosity of these oil products. The findings from the 4.0 mm results were incorporated with data from the literature to provide an update on the factors to consider during the decision making for spills of diluted bitumen products. The matrix indicated that most response options, including chemical dispersants, work much more effectively within 48 hours of the initiation of weathering. After this window of opportunity closes, natural attenuation or in situ burning is the only option remaining, but containment of oil is a limiting factor for in situ burning.

Shape memory polymers (SMPs) are soft active materials, their special property is the ability to hold a temporary shape and when exposed to a suitable trigger, they come back to their original shape. These external stimuli can be temperature, light or electro-magnetic fields. Amorphous SMPs are a class of thermally-activated SMPs that rely on glass transition to retain their temporary shape. Above the glass transition temperature (T > Tg), (amorphous SMPs exhibit finite deformation and viscoelastic behavior. In this work we develop a model to capture the viscoelastic behavior of the amorphous SMPs at elevated temperatures. The model uses an approach that was initially developed to study non-Newtonian viscoelastic fluids. We accomplish this by developing a multi-branch model based on the theory of multiple natural configurations using the maximization of the rate dissipation to determine the evolution of the natural configurations. We apply our model to study several different deformations at an elevated temperature T = 130 °C and show that this approach is able to capture the viscoelastic behavior of these polymers. The predictions of the theory are then compared with experimental results.