Abstract: Shorelines are vulnerable to significant environmental impacts from oil spills. Stranded oil can cause potential short- to long-term detrimental effects along beaches that include injuries to ecosystem, socio-economic and cultural resources. In this study, a three-dimensional (3D) numerical modeling approach is developed to evaluate the fate and transport of spilled oil for hypothetical oiled shoreline cases under various combinations of beach geomorphology and environmental conditions. The developed model estimates the spatial and temporal distribution of spilled oil for the various test conditions, using the finite volume method and considering the physical transport (dispersion and advection), sinks, and sorption processes. The model includes a user-friendly interface for data input on variables such as beach properties, environmental conditions, and physical-chemical properties of spilled oil. An experimental meso-scale tank design was used to test the developed model for dissolved petroleum hydrocarbon within shorelines. The simulated results for effects of different sediment substrates, oil types, and shoreline features for the transport of spilled oil are comparable to that obtained with a commercially available model. Results show that the properties of substrates and the oil removal by shoreline effects have significant impacts on oil transport in the beach area. Sensitivity analysis, through the application of the one-step-at-a-time method (OAT), for the 3D model identified hydraulic conductivity as the most sensitive parameter. The 3D numerical model allows users to examine the behavior of oil on and within beaches, assess potential environmental impacts, and provide technical support for decisions related to shoreline clean-up operations.
Abstract: The efforts to understand the heat transfer behavior of supercritical water in supercritical water cooled reactor (SCWR) are ongoing worldwide to fulfill the future energy demand. The higher thermal efficiency of these reactors compared to a conventional nuclear reactor is one of the driving forces for attracting the attention of nuclear scientists. In this work, a solution procedure has been described for solving supercritical fluid flow problems in complex geometries. The solution procedure is based on non-staggered grid. All governing equations are discretized by finite volume method (FVM) in curvilinear coordinate system. Convective terms are discretized by first-order upwind scheme and central difference approximation has been used to discretize the diffusive parts. k-ε turbulence model with standard wall function has been employed. SIMPLE solution procedure has been implemented for the curvilinear coordinate system. Based on this solution method, 3-D Computational Fluid Dynamics (CFD) code has been developed. In order to demonstrate the capability of this CFD code in supercritical fluid flows, heat transfer to supercritical water in circular tubes has been considered as a test problem. Results obtained by code have been compared with experimental results reported in literature.
Abstract: Micro-mixers play an important role in the lab-on-a-chip applications and micro total analysis systems to acquire the correct level of mixing for any given process. The mixing process can be classified as active or passive according to the use of external energy. Literature of microfluidics reports that most of the work is done on the models of steady laminar flow; however, the study of unsteady laminar flow is an active area of research at present. There are wide applications of this, out of which, we consider nanoparticle synthesis in micro-mixers. In this work, we have developed a model for unsteady flow to study the mixing performance of a passive micro mixer for reactants used for such synthesis. The model is developed in Finite Volume Method (FVM)-based software, OpenFOAM. The model is tested by carrying out the simulations at Re of 0.5. Mixing performance of the micro-mixer is investigated using simulated concentration values of mixed species across the width of the micro-mixer and calculating the variance across a line profile. Experimental validation is done by passing dyes through a Y shape micro-mixer fabricated using polydimethylsiloxane (PDMS) polymer and comparing variances with the simulated ones. Gold nanoparticles are later synthesized through the micro-mixer and collected at two different times leading to significantly different size distributions. These times match with the time scales over which reactant concentrations vary as obtained from simulations. Our simulations could thus be used to create design aids for passive micro-mixers used in nanoparticle synthesis.
Abstract: An essential component of a finite volume method (FVM) is the advection scheme that estimates values on the cell faces based on the calculated values on the nodes or cell centers. The most widely used advection schemes are upwind schemes. These schemes have been developed in FVM on different kinds of structured and unstructured grids. In this research, the physical influence scheme (PIS) is developed for a cell-centered FVM that uses an implicit coupled solver. Results are compared with the exponential differencing scheme (EDS) and the skew upwind differencing scheme (SUDS). Accuracy of these schemes is evaluated for a lid-driven cavity flow at Re = 1000, 3200, and 5000 and a backward-facing step flow at Re = 800. Simulations show considerable differences between the results of EDS scheme with benchmarks, especially for the lid-driven cavity flow at high Reynolds numbers. These differences occur due to false diffusion. Comparing SUDS and PIS schemes shows relatively close results for the backward-facing step flow and different results in lid-driven cavity flow. The poor results of SUDS in the lid-driven cavity flow can be related to its lack of sensitivity to the pressure difference between cell face and upwind points, which is critical for the prediction of such vortex dominant flows.
Abstract: In this work, numerical simulations were carried out using a specific CFD code in order to study the performance of an innovative Scraped Surface Heat Exchanger (SSHE) with helical ribbons for Bingham fluids (threshold fluids). The resolution of three-dimensional form of the conservation equations (continuity, momentum and energy equations) was carried out basing on the finite volume method (FVM). After studying the effect of dimensionless numbers (axial Reynolds, rotational Reynolds and Oldroyd numbers) on the hydrodynamic and thermal behaviors within SSHE, a parametric study was developed, by varying the width of the helical ribbon, the clearance between the stator wall and the tip of the ribbon and the number of turns of the helical ribbon, in order to improve the heat transfer inside the exchanger. The effect of these geometrical numbers on the hydrodynamic and thermal behaviors was discussed.
Abstract: The objective of this research work is to investigate
for one dimensional transient radiative transfer equations with
conduction using finite volume method. Within the infrastructure of
finite-volume, we obtain the conservative discretization of the terms
in order to preserve the overall conservative property of finitevolume
schemes. Coupling of conductive and radiative equation
resulting in fluxes is governed by the magnitude of emissivity,
extinction coefficient, and temperature of the medium as well as
geometry of the problem.
The problem under consideration has been solved, for a slab
dominating radiation coupled with transient conduction based on
finite volume method. The boundary conditions are also chosen so as
to give a good model of the discretized form of radiation transfer
equation. The important feature of the present method is flexibility in
specifying the control angles in the FVM, while keeping the
simplicity in the solution procedure.
Effects of various model parameters are examined on the
distributions of temperature, radiative and conductive heat fluxes and
incident radiation energy etc. The finite volume method is considered
to effectively evaluate the propagation of radiation intensity through
a participating medium.
Abstract: This paper discusses the applicability of the numerical model for a damage prediction method of the accidental hydrogen explosion occurring in a hydrogen facility. The numerical model was based on an unstructured finite volume method (FVM) code “NuFD/FrontFlowRed”. For simulating unsteady turbulent combustion of leaked hydrogen gas, a combination of Large Eddy Simulation (LES) and a combustion model were used. The combustion model was based on a two scalar flamelet approach, where a G-equation model and a conserved scalar model expressed a propagation of premixed flame surface and a diffusion combustion process, respectively. For validation of this numerical model, we have simulated the previous two types of hydrogen explosion tests. One is open-space explosion test, and the source was a prismatic 5.27 m3 volume with 30% of hydrogen-air mixture. A reinforced concrete wall was set 4 m away from the front surface of the source. The source was ignited at the bottom center by a spark. The other is vented enclosure explosion test, and the chamber was 4.6 m × 4.6 m × 3.0 m with a vent opening on one side. Vent area of 5.4 m2 was used. Test was performed with ignition at the center of the wall opposite the vent. Hydrogen-air mixtures with hydrogen concentrations close to 18% vol. were used in the tests. The results from the numerical simulations are compared with the previous experimental data for the accuracy of the numerical model, and we have verified that the simulated overpressures and flame time-of-arrival data were in good agreement with the results of the previous two explosion tests.
Abstract: In this study, the three-dimensional cavitating
turbulent flow in a complete Francis turbine is simulated using
mixture model for cavity/liquid two-phase flows. Numerical analysis
is carried out using ANSYS CFX software release 12, and standard k-ε
turbulence model is adopted for this analysis. The computational
fluid domain consist of spiral casing, stay vanes, guide vanes, runner
and draft tube. The computational domain is discretized with a threedimensional
mesh system of unstructured tetrahedron mesh. The
finite volume method (FVM) is used to solve the governing equations
of the mixture model. Results of cavitation on the runner’s blades
under three different boundary conditions are presented and
discussed. From the numerical results it has been found that the
numerical method was successfully applied to simulate the cavitating
two-phase turbulent flow through a Francis turbine, and also
cavitation is clearly predicted in the form of water vapor formation
inside the turbine. By comparison the numerical prediction results
with a real runner; it’s shown that the region of higher volume
fraction obtained by simulation is consistent with the region of runner
cavitation damage.
Abstract: The radiative heat transfer problem is investigated numerically for 2D complex geometry biomass pyrolysis reactor composed of two pyrolysis chambers and a heat recuperator. The fumes are a mixture of carbon dioxide and water vapor charged with absorbing and scattering particles and soot. In order to increase gases residence time and heat transfer, the heat recuperator is provided with many inclined, vertical, horizontal, diffuse and grey baffles of finite thickness and has a complex geometry. The Finite Volume Method (FVM) is applied to study radiative heat transfer. The blocked-off region procedure is used to treat the geometrical irregularities. Eight cases are considered in order to demonstrate the effect of adding baffles on the walls of the heat recuperator and on the walls of the pyrolysis rooms then choose the best case giving the maximum heat flux transferred to the biomass in the pyrolysis chambers. Ray effect due to the presence of baffles is studied and demonstrated to have a crucial effect on radiative heat flux on the walls of the pyrolysis rooms. Shadow effect caused by the presence of the baffles is also studied. The non grey radiative heat transfer is studied for the real existent configuration. The Weighted Sum of The Grey Gases (WSGG) Model of Kim and Song is used as non grey model. The effect of soot volumetric fraction on the non grey radiative heat flux is investigated and discussed.
Abstract: The practical application of the Computational Fluid Dynamics (CFD), for predicting the flow pattern around Multipurpose Amphibious Vehicle (MAV) hull has made much progress over the last decade. Today, several of the CFD tools play an important role in the land and water going vehicle hull form design. CFD has been used for analysis of MAV hull resistance, sea-keeping, maneuvering and investigating its variation when changing the hull form due to varying its parameters, which represents a very important task in the principal and final design stages. Resistance analysis based on CFD (Computational Fluid Dynamics) simulation has become a decisive factor in the development of new, economically efficient and environmentally friendly hull forms. Three-dimensional finite volume method (FVM) based on Reynolds Averaged Navier-Stokes equations (RANS) has been used to simulate incompressible flow around three types of MAV hull bow models in steady-state condition. Finally, the flow structure and streamlines, friction and pressure resistance and velocity contours of each type of hull bow will be compared and discussed.
Abstract: Numerical investigation on the generality of nanoparticle velocity equation had been done on the previous published work. The three dimensional governing equations (continuity, momentum and energy) were solved using finite volume method (FVM). Parametric study of thermal performance between pure water-cooled and nanofluid-cooled are evaluated for volume fraction in the range of 1% to 4%, and nanofluid type of gamma-Al2O3 at Reynolds number range of 67.41 to 286.77. The nanofluid is modeled using single and two phase approach. Three different existing Brownian motion velocities are applied in comparing the generality of the equation for a wide parametric condition. Deviation in between the Brownian motion velocity is identified to be due to the different means of mean free path and constant value used in diffusion equation.
Abstract: In this research, a 2-D computational analysis of
steady state free convection in a rectangular enclosure filled with an
electrically conducting fluid under Effect of Magnetic Field has been
performed. The governing equations (mass, momentum, and energy)
are formulated and solved by a finite volume method (FVM)
subjected to different boundary conditions. A parametric study has
been conducted to consider the influence of Grashof number (Gr),
Prantdl number (Pr) and the orientation of magnetic field on the flow
and heat transfer characteristics. It is observed that Nusselt number
(Nu) and heat flux will increase with increasing Grashof and Prandtl
numbers and decreasing the slope of the orientation of magnetic field.
Abstract: Calcium [Ca2+] is an important second messenger
which plays an important role in signal transduction. There are
several parameters that affect its concentration profile like buffer
source etc. The effect of stationary immobile buffer on Ca2+
concentration has been incorporated which is a very important
parameter needed to be taken into account in order to make the
model more realistic. Interdependence of all the important parameters
like diffusion coefficient and influx over [Ca2+] profile has been
studied. Model is developed in the form of advection diffusion
equation together with buffer concentration. A program has been
developed using finite volume method for the entire problem and
simulated on an AMD-Turion 32-bit machine to compute the
numerical results.
Abstract: To calculate the temperature distribution of the slab in
a hot rolled reheating furnace a mathematical model has been
developed by considering the thermal radiation in the furnace and
transient conduction in the slab. The furnace is modeled as radiating
medium with spatially varying temperature. Radiative heat flux within
the furnace including the effect of furnace walls, combustion gases,
skid beams and buttons is calculated using the FVM and is applied as
the boundary condition of the transient conduction equation of the
slab. After determining the slab emissivity by comparison between
simulation and experimental work, variation of heating characteristics
in the slab is investigated in the case of changing furnace temperature
with various time and the slab residence time is optimized with this
evaluation.
Abstract: Channel junctions can be analyzed in two ways of
division (lateral intake) and combined flows (confluence). The
present paper investigates 3D flow pattern at lateral intake using
Navier-Stokes equation and κ -ε (RNG) turbulent model. The
equations are solved by Finite-Volume Method (FVM) and results
are compared with the experimental data of (Barkdoll, B.D., 1997)
to test the validity of the findings. Comparison of the results with
the experimental data indicated a close proximity between the two
sets of data which suggest a very close simulation. Results further
indicated an inverse relation between the effects of discharge ratio
( r Q ) on the length and width of the separation zone. In other words,
as the discharge ration increases, the length and width of separation
zone decreases.
Abstract: Mathematical and computational modeling of calcium
signalling in nerve cells has produced considerable insights into how
the cells contracts with other cells under the variation of biophysical
and physiological parameters. The modeling of calcium signaling in
astrocytes has become more sophisticated. The modeling effort has
provided insight to understand the cell contraction. Main objective
of this work is to study the effect of voltage gated (Operated)
calcium channel (VOC) on calcium profile in the form of advection
diffusion equation. A mathematical model is developed in the form
of advection diffusion equation for the calcium profile. The model
incorporates the important physiological parameter like diffusion
coefficient etc. Appropriate boundary conditions have been framed.
Finite volume method is employed to solve the problem. A program
has been developed using in MATLAB 7.5 for the entire problem
and simulated on an AMD-Turion 32-bite machine to compute the
numerical results.
Abstract: This paper analytically investigates the 3D flow
pattern at the confluences of two rectangular channels having 900
angles using Navier-Stokes equations based on Reynolds Stress
Turbulence Model (RSM). The equations are solved by the Finite-
Volume Method (FVM) and the flow is analyzed in terms of steadystate
(single-phased) conditions. The Shumate experimental findings
were used to test the validity of data. Comparison of the simulation
model with the experimental ones indicated a close proximity
between the flow patterns of the two sets. Effects of the discharge
ratio on separation zone dimensions created in the main-channel
downstream of the confluence indicated an inverse relation, where a
decrease in discharge ratio, will entail an increase in the length and
width of the separation zone. The study also found the model as a
powerful analytical tool in the feasibility study of hydraulic
engineering projects.
Abstract: A cell-centered finite volume scheme for discretizing diffusion operators on distorted quadrilateral meshes has recently been designed and added to APMFCG to enable that code to be used as a tool for studying explosive magnetic flux compression generators. This paper describes this scheme. Comparisons with analytic results for 2-D test cases are presented, as well as 2-D results from a test of a "realistic" generator configuration.
Abstract: Pressure driven microscale gas flow-separation has
been investigated by solving the compressible Navier-Stokes (NS)
system of equations. A two dimensional explicit finite volume (FV)
compressible flow solver has been developed using modified
advection upwind splitting methods (AUSM+) with no-slip/first
order Maxwell-s velocity slip conditions to predict the flowseparation
behavior in microdimensions. The effects of scale-factor
of the flow geometry and gas species on the microscale gas flowseparation
have been studied in this work. The intensity of flowseparation
gets reduced with the decrease in scale of the flow
geometry. In reduced dimension, flow-separation may not at all be
present under similar flow conditions compared to the larger flow
geometry. The flow-separation patterns greatly depend on the
properties of the medium under similar flow conditions.