A CFD Model to Evaluate Near-Surface Oil Spill from a Broken Loading Pipe in Shallow Coastal Waters

Oil spills continue to generate various issues and concerns regarding their effect and behavior in the marine environment,owing to the related potential for detrimental environmental,economic and social implications.It is essential to have a solid understanding of the ways in which oil interacts with the water and the coastal ecosystems that are located nearby.This study proposes a simplified model for predicting the plume-like transport behavior of heavy Bunker C fuel oil discharging downward from an acutely-angled broken pipeline located on the water surface.The results show that the spill overall profile is articulated in three major flow areas.The first,is the source field,i.e.,a region near the origin of the initial jet,followed by the intermediate or transport field,namely,the region where the jet oil flow transitions into an underwater oil plume flow and starts to move horizontally,and finally,the far-field,where the oil re-surface and spreads onto the shore at a significant distance from the spill site.The behavior of the oil in the intermediate field is investigated using a simplified injection-type oil spill model capable of mimicking the undersea trapping and lateral migration of an oil plume originating from a negatively buoyant jet spill.A rectangular domain with proper boundary conditions is used to implement the model.The Projection approach is used to discretize a modified version of the Navier-Stokes equations in two dimensions.A benchmark fluid flow issue is used to verify the model and the results indicate a reasonable relationship between specific gravity and depth as well as agreement with the aerial data and a vertical temperature profile plot.

On numerical modeling of low-head direct chill ingot caster for magnesium alloy AZ31

A comprehensive 3D turbulent CFD study has been carried out to simulate a Low-Head(LH)vertical Direct Chill(DC)rolling ingot caster for the common magnesium alloy AZ31.The model used in this study takes into account the coupled laminar/turbulent melt flow and solidification aspects of the process and is based on the control-volume finite-difference approach.Following the aluminum/magnesium DC casting industrial practices,the LH mold is taken as 30 mm with a hot top of 60 mm.The previously verified in-house code has been modified to model the present casting process.Important quantitative results are obtained for four casting speeds,for three inlet melt pouring temperatures(superheats)and for three metal-mold contact heat transfer coefficients for the steady state operational phase of the caster.The variable cooling water temperatures reported by the industry are considered for the primary and secondary cooling zones during the simulations.Specifically,the temperature and velocity fields,sump depth and sump profiles,mushy region thickness,solid shell thickness at the exit of the mold and axial temperature profiles at the center and at three strategic locations at the surface of the slab are presented and discussed.

A coupled CFD simulation approach for investigating the pyrolysis process in industrial naphtha thermal cracking furnaces

In the steam thermal cracking of naphtha,the hydrocarbon stream flows inside tubular reactors and is exposed to flames of a series of burners in the firebox.In this paper,a full three-dimensional computational fluid dynamics(CFD)model was developed to investigate the process variables in the firebox and reactor coil of an industrial naphtha furnace.This comprehensive CFD model consists of a standard k-εturbulence model accompanied by a molecular kinetic reaction for cracking,detailed combustion model,and radiative properties.In order to improve the steam cracking performance,the model is solved using a proposed iterative algorithm.With respect to temperature,product yield and specially propylene-toethylene ratio(P/E),the simulation results agreed well with industrial data obtained from a mega olefin plant of a petrochemical complex.The deviation of P/E results from industrial data was less than 2%.The obtained velocity,temperature,and concentration profiles were used to investigate the residence time,coking rate,coke concentration,and some other findings.The coke concentration at coil exit was1.9×10^(-3)%(mass)and the residence time is calculated to be 0.29 s.The results can be used as a scientific guide for process engineers.

Simulating Urban Flow and Dispersion in Beijing by Coupling a CFD Model with the WRF Model

The airflow and dispersion of a pollutant in a complex urban area of Beijing, China, were numerically examined by coupling a Computational Fluid Dynamics （CFD） model with a mesoscale weather model. The models used were Open Source Field Operation and Manipulation （OpenFOAM） software package and Weather Research and Forecasting （WRF） model. OpenFOAM was firstly validated against wind-tunnel experiment data. Then, the WRF model was integrated for 42 h starting from 0800 LST 08 September 2009, and the coupled model was used to compute the flow fields at 1000 LST and 1400 LST 09 September 2009. During the WRF-simulated period, a high pressure system was dominant over the Beijing area. The WRF-simulated local circulations were characterized by mountain valley winds, which matched well with observations. Results from the coupled model simulation demonstrated that the airflows around actual buildings were quite different from the ambient wind on the boundary provided by the WRF model, and the pollutant dispersion pattern was complicated under the influence of buildings. A higher concentration level of the pollutant near the surface was found in both the step-down and step-up notches, but the reason for this higher level in each configurations was different： in the former, it was caused by weaker vertical flow, while in the latter it was caused by a downward-shifted vortex. Overall, the results of this study suggest that the coupled WRF-OpenFOAM model is an important tool that can be used for studying and predicting urban flow and dispersions in densely built-up areas.

CFD Model of Dense Gas-solid Systems in Jetting Fluidized Beds

A CFD code has been developed based on the conservation principles describing gas and solid flow in fluidized beds. This code is employed to simulate not only the spatiotemporal gas and solid phase velocities and voidage profiles in a two dimensional bed but also fluid dynamics in the jet region. The computational results show that gas flow direction is upward in the entire bed accompanied with random local circulations, whilst solid flow direction is upward at the center and downward near the wall. The radical reason of strong back mixing of solid particles and good transfer behavior between two phases is that the jet entrains solid particles. Numerical calculation indicates that gas velocity, solid velocity and pressure profile have a significant change when the voidage is 0 8. The simulated time averaged voidage profiles agree with the experimental results and simulated data reported by Gidaspow and Ettehadieh(1983). Therefore, CFD model can be regarded as a useful tool to study the jet characteristics in dense gas solid fluidized beds.

Design and Characterization of an Aerosol Test Chamber for Emergency Response Patient Contamination Control Simulation and Research

Contaminated or infected patients present a risk of cross-contamination for emergency responders, attending medical personnel and medical facilities as they enter a treatment facility. The controlled conditions of an aerosol test chamber are required to examine factors of contamination, decontamination, and cross-contamination. This study presents the design, construction, and a method for characterizing an aerosol test chamber for a full-sized manikin on a standard North Atlantic Treaty Organization litter. The methodology combined air velocity measurements, aerosol particle counts and size distributions, and computational fluid dynamics modeling to describe the chamber’s performance in three dimensions. This detailed characterization facilitates future experimental design by predicting chamber performance for a variety of patient-focused research.

Modelling Fire in Tunnels:A Large Scale Validated Two-Step Modelling Method

Fire is a quite common phenomenon in tunnel and being able to model its consequences with a good precision is crucial to design adapted safety measures.Modelling the fire behaviour in tunnel is quite challenging.Managing large-scale experiment for all the possible configurations is however economically unrealistic.This paper presents an experimental real scale fire test that was used not only for demonstrating the fire behaviour but also for evaluating the capabilities of the FDS(fire dynamics simulator)fire code to model fire consequences too.It enables highlighting the importance of wall and inlet boundary condition treatment.Keeping in mind that predicting fire development using a CFD(computational fluid dynamics)code is quite impossible,a two-level approach is discussed with an analytical model to predict the fire curve and a CFD model for predicting smoke propagation,temperature and toxic gases distribution inside the tunnel.The comparisons show a good agreement between experimental fire test and CFD modelling but also let appear requirements when using CFD.

Semi-continuous casting of magnesium alloy AZ91 using a filtered melt delivery system

A 3-D numerical simulation of an industrial-sized slab caster for magnesium alloy AZ91 has been carried out for the steady state operational phase of the caster.The simulated model consists of an open-top melt delivery system fitted with a porous filter near the hot-top.The melt flow through the porous filter was modeled on the basis of Brinkmann-Forchimier-Extended non-Darcy model for turbulent flow.An in-house 3-D CFD code was modified to account for the melt flow through the porous filter.Results are obtained for four casting speeds namely,40,60,80,and 100 mm/min.The metal-mold contact region as well as the convective heat transfer coefficient at the mold wall were also varied.In addition to the above,the Darcy number for the porous media was also changed.All parametric studies were performed for a fixed inlet melt superheat of 64℃.The results are presented pictorially in the form of temperature and velocity fields.The sump depth,mushy region thickness,solid shell thickness at the exit of the mold and axial temperature profiles are also presented and correlated with the casting speed through regression analysis.©2015 Production and hosting by Elsevier B.V.on behalf of Chongqing University.

Effects of Street-Bottom and Building-Roof Heating on Flow in Three-Dimensional Street Canyons

Using a computational fluid dynamics （CFD） model, the effects of street-bottom and building-roof heating on flow in three-dimensional street canyons are investigated. The building and street-canyon aspect ratios are one. In the presence of street-bottom heating, as the street-bottom heating intensity increases, the mean kinetic energy increases in the spanwise street canyon formed by the upwind and downwind buildings but decreases in the lower region of the streamwise street canyon. The increase in momentum due to buoyancy force intensifies mechanically induced flow in the spanwise street canyon. The vorticity in the spanwise street canyon strengthens. The temperature increase is not large because relatively cold above-roof-level air comes into the spanwise street canyon. In the presence of both street-bottom and building-roof heating, the mean kinetic energy rather decreases in the spanwise street canyon. This is caused by the decrease in horizontal flow speed at the roof level, which results in the weakening of the mean flow circulation in the spanwise street canyon. It is found that the vorticity in the spanwise street canyon weakens. The temperature increase is relatively large compared with that in the street-bottom heating case, because relatively warm above-roof-level air comes into the spanwise street canyon.

Numerical Study of Flow and Gas Diffusion in the Near-Wake behind an Isolated Building

To assist validation of the experimental data of urban pollution dispersion, the effect of an isolated building on the flow and gaseous diffusion in the wake region have been investigated numerically in the neutrally stratified rough-walled turbulent boundary layer. Numerical studies were carried out using Computational Fluid Dynamics （CFD） models. The CFD models used for the simulation were based on the steady-state Reynolds-Average Navier-Stoke equations （RANS） with κ-ε turbulence models; standard κ-ε and RNG κ-ε models. Inlet conditions and boundary conditions were specified numerically to the best information available for each fluid modeling simulation. A gas pollutant was emitted from a point source within the recirculation cavity behind the building model. The accuracy of these simulations was examined by comparing the predicted results with wind tunnel experimental data. It was confirmed that simulation using the model accurately reproduces the velocity and concentration diffusion fields with a fine-mish resolution in the near wake region. Results indicated that there is a good agreement between the numerical simulation and the wind tunnel experiment for both wind flow and concentration diffusion. The results of this work can help to improve the understanding of mechanisms of and simulation of pollutant transport in an urban environment.

Analysis of diffusion behavior of harmful emissions from trackless rubber-wheel diesel vehicles in underground coal mines

To define the diffusion behavior of harmful exhaust substances from diesel vehicles and support safety risk assessments of underground coal mines,we performed a multi-species coupling calculation of the emission and diffusion of harmful substances from a trackless rubber-wheel diesel vehicle.A computational fluid dynamics(CFD)model of the diffusion of harmful emissions was hence established and verified.From the perspective of risk analysis,the diffusion behavior and distribution of hazardous substances emitted by the diesel vehicle were studied under 4 different conditions;moreover,we identified areas characterized by hazardous levels of emissions.When the vehicle idled upwind in the roadway,high-risk areas formed behind and to the right of the vehicle:particularly high concentrations of pollutants were measured near the rear floor of the vehicle and within 5 m behind the vehicle.When the vehicle idled downwind,high-risk areas formed in front of it:particularly high concentrations of pollutants were measured near the floor and within 5 m from the front of the vehicle.In the above cases,the driver would not breathe highly polluted air and would be relatively safe.When the vehicle idled into the chamber,however,high-risk areas formed on both sides of the vehicle and near the upper roof.Forward entry of the vehicle caused a greater increase in the concentration of pollutants in the chamber and in the driver’s breathing zone compared with reverse entry.

Effects of a Building’s Density on Flow in Urban Areas

The effects of a building＇s density on urban flows are investigated using a CFD model with the RNG k - ε turbulence closure scheme. Twenty-seven cases with different building＇s density parameters （e.g., building and street-canyon aspect ratios） are numerically simulated. As the building＇s density parameters vary, different flow regimes appear. When the street canyon is relatively narrow and high, two counterrotating vortices in the vertical direction are generated. The wind speed along streets is mainly affected by the building＇s length. However, it is very difficult to find or generalize the characteristics of the street-canyon flows in terms of a single building＇s density parameter. This is because the complicated flow patterns appear due to the variation of the vortex structure and vortex number. Volume-averaged vorticity magnitude is a very good indicator to reflect the flow characteristics despite the strong dependency of flows on the variation of the building＇s density parameters. Multi-linear regression shows that the volume-averaged vorticity magnitude is a strong function of the building＇s length and the street-canyon width. The increase in the building＇s length decreases the vorticity of the street-canyon flow, while, the increase in the street- canyon width increases the vorticity.

CAVITATION MODELING WITH A CFD APPROACH

A CFD model of two-dimensional cavity flow is developed using the full Navier-Stokes equations. Based on pseudo-compressibillty and time marching techniques, the cavity surface evolves with the flow field during the time marching process. The cavitation boundary conditions are enforced on the cavity surface. Numerical computations are carried out for a wide range of two-dimensional cavity flows over various hydrofoils, including leading edge cavities and midchord cavities. Extensions to three dimensional nows are, in principle, straight forward.

Effects of Wind Fences on the Wind Environment around Jang Bogo Antarctic Research Station

This study investigated the flow characteristics altered by Jang Bogo Antarctic Research Station using computational fluid dynamics(CFD) modeling. The topography and buildings around Jang Bogo Station were constructed with computeraided-design data in the CFD model domain. We simulated 16 cases with different inflow directions, and compared the flow characteristics with and without Jang Bogo Station for each inflow direction. The wind data recorded by the site’s automatic weather station(AWS) were used for comparison. Wind rose analysis showed that the wind speed and direction after the construction of Jang Bogo Station were quite different from those before construction. We also investigated how virtual wind fences would modify the flow patterns, changing the distance of the fence from the station as well as the porosity of the fence. For westerly inflows, when the AWS was downwind of Jang Bogo Station, the decrease in wind speed was maximized(-81% for west-northwesterly). The wind speed reduction was also greater as the distance of the fence was closer to Jang Bogo Station. With the same distance, the fence with medium porosity(25%–33%) maximized the wind speed reduction.These results suggest that the location and material of the wind fence should be selected carefully, or AWS data should be interpreted cautiously, for particular prevailing wind directions.

Study on hemodynamics in patient-specific thoracic aortic aneurysm

The objective of this study is to investigate the hemodynamics in patient-specific thoracic aortic aneurysm and discuss the reason for formation of aortic plaque.A 3-Dimensional pulsatile blood flow in thoracic aorta with a fusiform aneurysm and 3 main branched vessels was studied numerically with the average Reynolds number of 1399 and the Womersley number of 19.2.Based on the clinical 2-Dimensional CT slice data,the patient-specific geometry model was constructed using medical image process software.Unsteady,incompressible,3-Dimensional Navier-Stokes equations were employed to solve the flow field.The temporal distributions of hemodynamic variables during the cardiac cycle such as streamlines,wall shear stresses in the arteries and aneurysm were analyzed. Growth and rupture mechanisms of thoracic aortic aneurysm in the patient can be analyzed based on patient-specific model and hemodynamics simulation.

AXIAL HEAT CONDUCTION MODEL TO PREDICT MAXIMUM HEAT REMOVE OF MINIATURE HEAT PIPE BASED ON GREY MODEL THEORY

Computer chip is always accompanied by the increase of heat dissipation and miniaturization. The miniature heat pipes are widely used in notebook computer to resolve the heat dissipation problems. Maximum heat removed model of miniature heat pipes building by grey model is presented. In order to know the foundation for modeling, the smooth grade of error examination is inquired and the accuracy of grey relational grade is verified. The model can be used to select a suitable heat pipes to solve electric heat problems in the future. Final results show that the grey model only needs four experiment data and its error value is less than 10%, further, it is better than computational fluid dynamics （CFD） model.

Using computational fluid dynamic simulation with Flow-3D to reveal the origin of the mushroom stone in the Xiqiao Mountain of Guangdong, China

Mushroom shaped rocks are not uncommon in nature, but their origin is often misunderstood because they can be formed by different natural forces. A huge mushroom stone in a stream valley of the Xiqiao Mountain of Guangdong,China is widely believed to be formed through the way of stream water erosion. However, the result of the simulation performed with Flow-3 D in this study indicates that the mushroom shape of the rock could not have been sculpted by the flowing water erosion because the simulation result contradicts that of the field investigations:(1) the sediment brought by upstream flowing water tends to be deposited on the leeward side of the flowing water preventing the rock on this side form being eroded, but the narrowest part of the rock bottom is exactly located on this side of the rock and(2) the stream flow should erode the bed sediment and produce high abrasion on the rock surface on the upstream and lateral sides of the rock,respectively, but no abrasion marks are found on the same sides of the rock. Subsequent geological analyses in this study provide evidence that the narrow bottom of the mushroom stone was resulted from chemical and physical weathering happening within the sediment that used to deposit on the lee side of the stone.

Channel Slope Effect on Energy Dissipation of Flow over Broad Crested Weirs

The main purpose of broad crested weir used in open channels is to raise and control upstream (U/S) water level. In this study, a new performance was added to this weir, by making a step at downstream (D/S) of weir. The energy dissipation, the height of the weir/the upstream water height ratio and Froude number relationships (E% – P/h – Fr) for three range of flume slop S = 0.0, 0.002 and 0.004 were simulated. The experiments were performed in a laboratory horizontal channel of 4.6 m length, 0.3 m width and 0.3 m depth for a wide range of discharge. The D/S step height of the weir was 7.5 cm. FLUENT software was used as numerical model which represent a type of Computational Fluid Dynamics (CFD) model in order to simulate flow over weirs. The Volume of Fluid (VOF) method with the Standard k – ε turbulence model was used to estimate the free surface profile and the structured mesh with high concentration near the wall regions. The experimental results of the water surface profile gave a high agreement with the results of the numerical models. The maximum value 28.78 of E% was obtained in single step broad crested weir in the experimental result and 27.35 in numerical result at S = 0.004. Finally, the range of the relative error of the energy dissipation between experimental and numerical results was achieved and the maximum was 6.76 in all runs.

Comparison of two-and three-dimensional modeling of invert trap for sewer solid management

In the present study, five different invert trap configurations （rectangular with and without lids on both sides; trapezoidal, trapezoidal with rectangular base and rectangular with trapezoidal base with lids on both sides） were simulated for both two-dimensional （2D） and three-dimensional （3D） flow conditions for three sediment types （sand, styrocell and plastic beads） at six flow rates （0.35, 0.70, 1.05, 1.35, 4.55 and 9.95 L/s） for each trap. Computational fluid dynamics （CFD）-based modeling using FLUENT software with Renormalization Group （RNG） k-e model along with discrete phase model （DPM） were used in the simulations. A hexagonal/tetrahedral and map-type non-uniform grid was chosen to discretize the entire computational domain and a control volume finite difference method was used to solve the governing equations. The flow rates selected in the present study cover the entire range of flow rate expected for dry weather and monsoon. The simulation is capable of differentiating between 2D and 3D modeling of particle trajectories, the effects of flow rate and trap geometry on flow patterns developed in the trap. The sediment retention ratio for 2D is higher than that for 3D modeling for all flow conditions, particle types and model geometry due to inclusion of lateral effects in 3D modeling. The invert trap having rectangular shape with trapezoidal base is found to be the most efficient configuration in both 2D and 3D modeling.

Effects of window opening style on inside environment of solar greenhouse based on CFD simulation

Solar greenhouse is extensively used in horticultural production in China.Natural ventilation is one of the key technological means to adjust the inside environment of greenhouse,while the effects of window opening styles on the inside environment are still not well understood.In the present study,the temperature and air flow field of five operation styles,namely back roof central opening style(G1),back roof evenly spaced opening style(G2),top roof full opening style(G3),style of G1+G3,and style of G2+G3 were simulated using the CFD method.The results indicated that:(1)the simulated and measured results exhibited favorable agreement,with relative errors within 5%;(2)In the case of the windows opening area was the same and only single ventilation style was applied,back roof full opening style exhibited the best cooling effect.The inside average temperature of the greenhouse with G1 style decreased by 0.5°C and 1.6°C respectively compared with those of greenhouses with G2 and G3.(3)The cooling effect of the style of G2+G3 was more favorable than that of the style of G1+G3.The style of G2+G3 exhibited better cooling effect than the single ventilation styles,with the lowest temperature(27.5°C)and temperature uniformity coefficient(0.36).