Numerical simulation of countercurrent flow based on two-fluid model

Energy Technology Data Exchange (ETDEWEB)

Chen, H.D. [Institute of Nuclear Engineering and Technology, Sun Yat-sen University, Zhuhai 519082 (China); School of Electric Power, South China University of Technology, Guangzhou 510640 (China); Zhang, X.Y., E-mail: zxiaoying@mail.sysu.edu.cn [Institute of Nuclear Engineering and Technology, Sun Yat-sen University, Zhuhai 519082 (China)

2017-03-15

Highlights: • Using one-dimensional two-fluid model to help understanding counter-current flow two-phase flows. • Using surface tension model to make the one-dimensional two-fluid flow model well-posed. • Solving the governing equations with a modified SIMPLE algorithm. • Validating code with experimental data and applying it to vertical air/steam countercurrent flow condition - Abstract: In order to improve the understanding of counter-current two-phase flows, a transient analysis code is developed based on one-dimensional two-fluid model. A six equation model has been established and a two phase pressure model with surface tension term, wall drag force and interface shear terms have been used. Taking account of transport phenomenon, heat and mass transfer models of interface were incorporated. The staggered grids have been used in discretization of equations. For validation of the model and code, a countercurrent air-water problem in one experimental horizontal stratified flow has been considered firstly. Comparison of the computed results and the experimental one shows satisfactory agreement. As the full problem for investigation, one vertical pipe with countercurrent flow of steam-water and air-water at same boundary condition has been taken for study. The transient distribution of liquid fraction, liquid velocity and gas velocity for selected positions of steam-water and air-water problem were presented and discussed. The results show that these two simulations have similar transient behavior except that the distribution of gas velocity for steam-water problem have larger oscillation than the one for air-water. The effect of mesh size on wavy characteristics of interface surface was also investigated. The mesh size has significant influence on the simulated results. With the increased refinement, the oscillation gets stronger.

Flow modelling of a newtonian fluid by two regions- the region of pure fluid and porous region

International Nuclear Information System (INIS)

Sampaio, R.; Gama, R.M.S. da

1983-01-01

A model of flow with two regions is presented using mixture theory. One region contains only pure fluid and the other a mixture of fluid and porous rigid solid. Compatibility conditons on the pure fluid-mixture interface are carefully discussed. The theory is used to solve a problem of a flow induced by pressure gradient and helicoidal motion of an impermeable cylinder on two rings one of pure fluid and another of mixture. (Author) [pt

Mathematical modeling of fluid flow in aluminum ladles for degasification with impeller - injector

Science.gov (United States)

Ramos-Gómez, E.; González-Rivera, C.; Ramírez-Argáez, M. A.

2012-09-01

In this work a fundamental Eulerian mathematical model was developed to simulate fluid flow in a water physical model of an aluminum ladle equipped with impeller for degassing treatment. The effect of critical process parameters such as rotor speed, gas flow rate on the fluid flow and vortex formation was analyzed with this model. Commercial CFD code PHOENICS 3.4 was used to solve all conservation equations governing the process for this twophase fluid flow system. The mathematical model was successfully validated against experimentally measured liquid velocity and turbulent profiles in a physical model. From the results it was concluded that the angular speed of the impeller is the most important parameter promoting better stirred baths. Pumping effect of the impeller is increased as impeller rotation speed increases. Gas flow rate is detrimental on bath stirring and diminishes pumping effect of impeller.

Fully coupled thermal-mechanical-fluid flow model for nonliner geologic systems

International Nuclear Information System (INIS)

Hart, R.D.

1981-01-01

A single model is presented which describes fully coupled thermal-mechanical-fluid flow behavior of highly nonlinear, dynamic or quasistatic, porous geologic systems. The mathematical formulation for the model utilizes the continuum theory of mixtures to describe the multiphase nature of the system, and incremental linear constitutive theory to describe the path dependency of nonlinear material behavior. The model, incorporated in an explicit finite difference numerical procedure, was implemented in two different computer codes. A special-purpose one-dimensional code, SNEAKY, was written for initial validation of the coupling mechanisms and testing of the coupled model logic. A general purpose commercially available code, STEALTH, developed for modeling dynamic nonlinear thermomechanical processes, was modified to include fluid flow behavior and the coupling constitutive model. The fully explicit approach in the coupled calculation facilitated the inclusion of the coupling mechanisms and complex constitutive behavior. Analytical solutions pertaining to consolidation theory for soils, thermoelasticity for solids, and hydrothermal convection theory provided verification of stress and fluid flow, stress and conductive heat transfer, and heat transfer and fluid flow couplings, respectively, in the coupled model. A limited validation of the adequacy of the coupling constitutive assumptions was also performed by comparison with the physical response from two laboratory tests. Finally, the full potential of the coupled model is illustrated for geotechnical applications in energy-resource related areas. Examples in the areas of nuclear waste isolation and cut-and-fill mining are cited

A computer model for dispersed fluid-solid turbulent flows

International Nuclear Information System (INIS)

Liu, C.H.; Tulig, T.J.

1985-01-01

A computer model is being developed to simulate two-phase turbulent flow phenomena in fluids containing finely dispersed solids. The model is based on a dual-continuum picture of the individual phases and an extension of a two-equation turbulence closure theory. The resulting set of nonlinear partial differential equations are solved using a finite difference procedure with special treatment to promote convergence. The model has been checked against a number of idealized flow problems with known solutions. The authors are currently comparing model predictions with measurements to determine a proper set of turbulence parameters needed for simulating two-phase turbulent flows

Modeling Fluid Flow in Faulted Basins

Directory of Open Access Journals (Sweden)

Faille I.

2014-07-01

Full Text Available This paper presents a basin simulator designed to better take faults into account, either as conduits or as barriers to fluid flow. It computes hydrocarbon generation, fluid flow and heat transfer on the 4D (space and time geometry obtained by 3D volume restoration. Contrary to classical basin simulators, this calculator does not require a structured mesh based on vertical pillars nor a multi-block structure associated to the fault network. The mesh follows the sediments during the evolution of the basin. It deforms continuously with respect to time to account for sedimentation, erosion, compaction and kinematic displacements. The simulation domain is structured in layers, in order to handle properly the corresponding heterogeneities and to follow the sedimentation processes (thickening of the layers. In each layer, the mesh is unstructured: it may include several types of cells such as tetrahedra, hexahedra, pyramid, prism, etc. However, a mesh composed mainly of hexahedra is preferred as they are well suited to the layered structure of the basin. Faults are handled as internal boundaries across which the mesh is non-matching. Different models are proposed for fault behavior such as impervious fault, flow across fault or conductive fault. The calculator is based on a cell centered Finite Volume discretisation, which ensures conservation of physical quantities (mass of fluid, heat at a discrete level and which accounts properly for heterogeneities. The numerical scheme handles the non matching meshes and guaranties appropriate connection of cells across faults. Results on a synthetic basin demonstrate the capabilities of this new simulator.

Collisionless kinetic-fluid model of zonal flows in toroidal plasmas

International Nuclear Information System (INIS)

Sugama, H.; Watanabe, T.-H.; Horton, W.

2006-12-01

A novel kinetic-fluid model is presented, which describes collisionless time evolution of zonal flows in tokamaks. In the new zonal-flow closure relations, the parallel heat fluxes are written by the sum of short- and long-time-evolution parts. The former part is given in the dissipative form of the parallel heat diffusion and relates to collisionless damping processes. The latter is derived from the long-time-averaged gyrocenter distribution and plays a major role in describing low-frequency or stationary zonal flows, for which the parallel heat fluxes are expressed in terms of the parallel flow as well as the nonlinear-source and initial-condition terms. It is shown analytically and numerically that, when applied to the zonal flow driven by either ion or electron temperature gradient turbulence, the kinetic-fluid equations including the new closure relations can reproduce the same long-time zonal-flow responses to the initial condition and to the turbulence source as those obtained from the gyrokinetic model. (author)

Beyond Poiseuille: Preservation Fluid Flow in an Experimental Model

Directory of Open Access Journals (Sweden)

Saurabh Singh

2013-01-01

Full Text Available Poiseuille’s equation describes the relationship between fluid viscosity, pressure, tubing diameter, and flow, yet it is not known if cold organ perfusion systems follow this equation. We investigated these relationships in an ex vivo model and aimed to offer some rationale for equipment selection. Increasing the cannula size from 14 to 20 Fr increased flow rate by a mean (SD of 13 (12%. Marshall’s hyperosmolar citrate was three times less viscous than UW solution, but flows were only 45% faster. Doubling the bag pressure led to a mean (SD flow rate increase of only 19 (13%, not twice the rate. When external pressure devices were used, 100 mmHg of continuous pressure increased flow by a mean (SD of 43 (17% when compared to the same pressure applied initially only. Poiseuille’s equation was not followed; this is most likely due to “slipping” of preservation fluid within the plastic tubing. Cannula size made little difference over the ranges examined; flows are primarily determined by bag pressure and fluid viscosity. External infusor devices require continuous pressurisation to deliver high flow. Future studies examining the impact of perfusion variables on graft outcomes should include detailed equipment descriptions.

Computational fluid dynamics incompressible turbulent flows

CERN Document Server

Kajishima, Takeo

2017-01-01

This textbook presents numerical solution techniques for incompressible turbulent flows that occur in a variety of scientific and engineering settings including aerodynamics of ground-based vehicles and low-speed aircraft, fluid flows in energy systems, atmospheric flows, and biological flows. This book encompasses fluid mechanics, partial differential equations, numerical methods, and turbulence models, and emphasizes the foundation on how the governing partial differential equations for incompressible fluid flow can be solved numerically in an accurate and efficient manner. Extensive discussions on incompressible flow solvers and turbulence modeling are also offered. This text is an ideal instructional resource and reference for students, research scientists, and professional engineers interested in analyzing fluid flows using numerical simulations for fundamental research and industrial applications. • Introduces CFD techniques for incompressible flow and turbulence with a comprehensive approach; • Enr...

Modeling of dilute and dense dispersed fluid-particle flow

Energy Technology Data Exchange (ETDEWEB)

Laux, Harald

1998-08-01

A general two-fluid model is derived and applied in CFD computations to various test cases of important industrial multiphase flows. It is general in the sense of its applicability to dilute and dense dispersed fluid-particle flows. The model is limited to isothermal flow without mass transfer and only one particle phase is described. The instantaneous fluid phase equations, including the phase interaction terms, are derived from a volume averaging technique, and the instantaneous particle phase equations are derived from the kinetic theory of granular material. Whereas the averaging procedure, the treatment of the interaction terms, and the kinetic theory approach have been reported in literature prior to this work the combination of the approaches is new. The resulting equations are derived without ambiguity in the interpretation of the particle phase pressure (equation-of-state of particle phase). The basic modeling for the particle phase is improved in two steps. Because in the basic modeling only stresses due to kinetic and collisional interactions are included, a simple model for an effective viscosity is developed in order to allow also frictional stresses within the particle phase. Moreover, turbulent stresses and turbulent dispersion of particles play often an important role for the transport processes. Therefore in a second step, a two-equation turbulence model for both fluid and particle phase turbulence is derived by applying the phasic average to the instantaneous equations. The resulting k-{epsilon}-k{sup d}-{epsilon}{sup d} model is new. Mathematical closure is attempted such that the resulting set of equations is valid for both dilute arid dense flows. During the development of the closure relations a clear distinction is made between granular or ''viscous'' microscale fluctuations and turbulent macro scale fluctuations (true particle turbulence) within the particle phase. The set of governing equations is discretized by using a finite volume method

Modeling of dilute and dense dispersed fluid-particle flow

Energy Technology Data Exchange (ETDEWEB)

Laux, Harald

1998-08-01

A general two-fluid model is derived and applied in CFD computations to various test cases of important industrial multiphase flows. It is general in the sense of its applicability to dilute and dense dispersed fluid-particle flows. The model is limited to isothermal flow without mass transfer and only one particle phase is described. The instantaneous fluid phase equations, including the phase interaction terms, are derived from a volume averaging technique, and the instantaneous particle phase equations are derived from the kinetic theory of granular material. Whereas the averaging procedure, the treatment of the interaction terms, and the kinetic theory approach have been reported in literature prior to this work the combination of the approaches is new. The resulting equations are derived without ambiguity in the interpretation of the particle phase pressure (equation-of-state of particle phase). The basic modeling for the particle phase is improved in two steps. Because in the basic modeling only stresses due to kinetic and collisional interactions are included, a simple model for an effective viscosity is developed in order to allow also frictional stresses within the particle phase. Moreover, turbulent stresses and turbulent dispersion of particles play often an important role for the transport processes. Therefore in a second step, a two-equation turbulence model for both fluid and particle phase turbulence is derived by applying the phasic average to the instantaneous equations. The resulting k-{epsilon}-k{sup d}-{epsilon}{sup d} model is new. Mathematical closure is attempted such that the resulting set of equations is valid for both dilute arid dense flows. During the development of the closure relations a clear distinction is made between granular or ''viscous'' microscale fluctuations and turbulent macro scale fluctuations (true particle turbulence) within the particle phase. The set of governing equations is discretized by using a

Mathematical Modelling of Fluid Flow in Cone and Cavitation Formation

Directory of Open Access Journals (Sweden)

Milada KOZUBKOVÁ

2011-06-01

Full Text Available Problem of cavitation is the undesirable phenomena occuring in the fluid flow in many hydraulic application (pumps, turbines, valves, etc.. Therefore this is in the focus of interest using experimental and mathematical methods. Based on cavitation modelling in Laval nozzle results and experience [1], [2], [4], following problem described as the water flow at the outlet from turbine blade wheel was solved. Primarily the problem is simplified into modelling of water flow in cone. Profiles of axial, radial and tangential velocity are defined on inlet zone. The value of pressure is defined on the outlet. Boundary conditions were defined by main investigator of the grant project – Energy Institute, Victor Kaplan’s Department of Fluid Engineering, Faculty of Mechanical Engineering, Brno University of Technology. The value of air volume was insignificant. Cavitation was solved by Singhal model of cavitation.

Simulation of horizontal pipe two-phase slug flows using the two-fluid model

Energy Technology Data Exchange (ETDEWEB)

Ortega Malca, Arturo J. [Pontificia Univ. Catolica do Rio de Janeiro, RJ (Brazil). Dept. de Engenharia Mecanica. Nucleo de Simulacao Termohidraulica de Dutos (SIMDUT); Nieckele, Angela O. [Pontificia Univ. Catolica do Rio de Janeiro, RJ (Brazil). Dept. de Engenharia Mecanica

2005-07-01

Slug flow occurs in many engineering applications, mainly in the transport of hydrocarbon fluids in pipelines. The intermittency of slug flow causes severe unsteady loading on the pipelines carrying the fluids, which gives rise to design problems. Therefore, it is important to be able to predict the onset and development of slug flow as well as slug characteristics. The present work consists in the simulation of two-phase flow in slug pattern through horizontal pipes using the two-fluid model in its transient and one-dimensional form. The advantage of this model is that the flow field is allowed to develop naturally from a given initial conditions as part of the transient calculation; the slug evolves automatically as a product of the computed flow development. Simulations are then carried out for a large number of flow conditions that lead a slug flow. (author)

The model coupling fluid flow in reservoir with flow in horizontal wellbore

Energy Technology Data Exchange (ETDEWEB)

Liu, Xiangping; Jiang, Zhixiang [RIPED-TEXACO Horizontal Well Technology Laboratory (United States)

1998-12-31

Three-dimensional pressure distributions of oil flow in a reservoir with horizontal well were derived, and a new formula to calculate pressure drop along the horizontal wellbore was developed based on the principle of conservation of matter and momentum. The formula considers the effect of influx into the horizontal wellbore from the reservoir on pressure drop in the wellbore. A mathematical model to couple fluid flow in the reservoir with flow in the horizontal wellbore is presented. Model results and experimental data showed good correspondence. Results showed the influence of pressure drop on well performance. 13 refs., 2 tabs., 7 figs.

Predictions of bubbly flows in vertical pipes using two-fluid models in CFDS-FLOW3D code

International Nuclear Information System (INIS)

Banas, A.O.; Carver, M.B.; Unrau, D.

1995-01-01

This paper reports the results of a preliminary study exploring the performance of two sets of two-fluid closure relationships applied to the simulation of turbulent air-water bubbly upflows through vertical pipes. Predictions obtained with the default CFDS-FLOW3D model for dispersed flows were compared with the predictions of a new model (based on the work of Lee), and with the experimental data of Liu. The new model, implemented in the CFDS-FLOW3D code, included additional source terms in the open-quotes standardclose quotes κ-ε transport equations for the liquid phase, as well as modified model coefficients and wall functions. All simulations were carried out in a 2-D axisymmetric format, collapsing the general multifluid framework of CFDS-FLOW3D to the two-fluid (air-water) case. The newly implemented model consistently improved predictions of radial-velocity profiles of both phases, but failed to accurately reproduce the experimental phase-distribution data. This shortcoming was traced to the neglect of anisotropic effects in the modelling of liquid-phase turbulence. In this sense, the present investigation should be considered as the first step toward the ultimate goal of developing a theoretically sound and universal CFD-type two-fluid model for bubbly flows in channels

Predictions of bubbly flows in vertical pipes using two-fluid models in CFDS-FLOW3D code

Energy Technology Data Exchange (ETDEWEB)

Banas, A.O.; Carver, M.B. [Chalk River Laboratories (Canada); Unrau, D. [Univ. of Toronto (Canada)

1995-09-01

This paper reports the results of a preliminary study exploring the performance of two sets of two-fluid closure relationships applied to the simulation of turbulent air-water bubbly upflows through vertical pipes. Predictions obtained with the default CFDS-FLOW3D model for dispersed flows were compared with the predictions of a new model (based on the work of Lee), and with the experimental data of Liu. The new model, implemented in the CFDS-FLOW3D code, included additional source terms in the {open_quotes}standard{close_quotes} {kappa}-{epsilon} transport equations for the liquid phase, as well as modified model coefficients and wall functions. All simulations were carried out in a 2-D axisymmetric format, collapsing the general multifluid framework of CFDS-FLOW3D to the two-fluid (air-water) case. The newly implemented model consistently improved predictions of radial-velocity profiles of both phases, but failed to accurately reproduce the experimental phase-distribution data. This shortcoming was traced to the neglect of anisotropic effects in the modelling of liquid-phase turbulence. In this sense, the present investigation should be considered as the first step toward the ultimate goal of developing a theoretically sound and universal CFD-type two-fluid model for bubbly flows in channels.

Modeling of magnetorheological fluid in quasi-static squeeze flow mode

Science.gov (United States)

Horak, Wojciech

2018-06-01

This work presents a new nonlinear model to describe MR fluid behavior in the squeeze flow mode. The basis for deriving the model were the principles of continuum mechanics and the theory of tensor transformation. The analyzed case concerned quasi-static squeeze with a constant area, between two parallel plates with non-slip boundary conditions. The developed model takes into account the rheological properties or MR fluids as a viscoplastic material for which yield stress increases due to compression. The model also takes into account the formation of normal force in the MR fluid as a result of the magnetic field impact. Moreover, a new parameter has been introduced which characterizes the behavior of MR fluid subjected to compression. The proposed model has been experimentally validated and the obtained results suggest that the assumptions made in the model development are reasonable, as good model compatibility with the experiments was obtained.

A critical review of the data requirements for fluid flow models through fractured rock

International Nuclear Information System (INIS)

Priest, S.D.

1986-01-01

The report is a comprehensive critical review of the data requirements for ten models of fluid flow through fractured rock, developed in Europe and North America. The first part of the report contains a detailed review of rock discontinuities and how their important geometrical properties can be quantified. This is followed by a brief summary of the fundamental principles in the analysis of fluid flow through two-dimensional discontinuity networks and an explanation of a new approach to the incorporation of variability and uncertainty into geotechnical models. The report also contains a review of the geological and geotechnical properties of anhydrite and granite. Of the ten fluid flow models reviewed, only three offer a realistic fracture network model for which it is feasible to obtain the input data. Although some of the other models have some valuable or novel features, there is a tendency to concentrate on the simulation of contaminant transport processes, at the expense of providing a realistic fracture network model. Only two of the models reviewed, neither of them developed in Europe, have seriously addressed the problem of analysing fluid flow in three-dimensional networks. (author)

Mathematical modeling for laminar flow of power law fluid in porous media

Energy Technology Data Exchange (ETDEWEB)

Silva, Renato A.; Mesquita, Maximilian S. [Universidade Federal do Espirito Santo (UFES), Sao Mateus, ES (Brazil). Centro Universitario Norte do Espirito Santo. Dept. de Engenharias e Computacao

2010-07-01

In this paper, the macroscopic equations for laminar power-law fluid flow is obtained for a porous medium starting from traditional equations (Navier-Stokes). Then, the volume averaging is applied in traditional transport equations with the power-law fluid model. This procedure leads to macroscopic transport equations set for non-Newtonian fluid. (author)

Computational modelling of a non-viscous fluid flow in a multi-walled carbon nanotube modelled as a Timoshenko beam

International Nuclear Information System (INIS)

Khosravian, N; Rafii-Tabar, H

2008-01-01

In the design of nanotube-based fluidic devices, a critical issue is the effect of the induced vibrations in the nanotube arising from the fluid flow, since these vibrations can promote structural instabilities, such as buckling transitions. It is known that the induced resonant frequencies depend on the fluid flow velocity in a significant manner. We have studied, for the first time, the flow of a non-viscous fluid in stubby multi-walled carbon nanotubes, using the Timoshenko classical beam theory to model the nanotubes as a continuum structure. We have obtained the variations of the resonant frequencies with the fluid flow velocity under several experimentally interesting boundary conditions and aspect ratios of the nanotube. The main finding from our work is that, compared to an Euler-Bernoulli classical beam model of a nanotube, the Timoshenko beam predicts the loss of stability at lower fluid flow velocities

Computational modelling of a non-viscous fluid flow in a multi-walled carbon nanotube modelled as a Timoshenko beam

Energy Technology Data Exchange (ETDEWEB)

Khosravian, N; Rafii-Tabar, H [Computational Physical Sciences Research Laboratory, Department of Nano-Science, Institute for Research in Fundamental Sciences (IPM), PO Box 19395-5531, Tehran (Iran, Islamic Republic of)], E-mail: rafii-tabar@nano.ipm.ac.ir

2008-07-09

In the design of nanotube-based fluidic devices, a critical issue is the effect of the induced vibrations in the nanotube arising from the fluid flow, since these vibrations can promote structural instabilities, such as buckling transitions. It is known that the induced resonant frequencies depend on the fluid flow velocity in a significant manner. We have studied, for the first time, the flow of a non-viscous fluid in stubby multi-walled carbon nanotubes, using the Timoshenko classical beam theory to model the nanotubes as a continuum structure. We have obtained the variations of the resonant frequencies with the fluid flow velocity under several experimentally interesting boundary conditions and aspect ratios of the nanotube. The main finding from our work is that, compared to an Euler-Bernoulli classical beam model of a nanotube, the Timoshenko beam predicts the loss of stability at lower fluid flow velocities.

Multiscale modeling of fluid flow and mass transport

Science.gov (United States)

Masuoka, K.; Yamamoto, H.; Bijeljic, B.; Lin, Q.; Blunt, M. J.

2017-12-01

In recent years, there are some reports on a simulation of fluid flow in pore spaces of rocks using Navier-Stokes equations. These studies mostly adopt a X-ray CT to create 3-D numerical grids of the pores in micro-scale. However, results may be of low accuracy when the rock has a large pore size distribution, because pores, whose size is smaller than resolution of the X-ray CT may be neglected. We recently found out by tracer tests in a laboratory using a brine saturated Ryukyu limestone and inject fresh water that a decrease of chloride concentration took longer time. This phenomenon can be explained due to weak connectivity of the porous networks. Therefore, it is important to simulate entire pore spaces even those of very small sizes in which diffusion is dominant. We have developed a new methodology for multi-level modeling for pore scale fluid flow in porous media. The approach is to combine pore-scale analysis with Darcy-flow analysis using two types of X-ray CT images in different resolutions. Results of the numerical simulations showed a close match with the experimental results. The proposed methodology is an enhancement for analyzing mass transport and flow phenomena in rocks with complicated pore structure.

Fluid friction and wall viscosity of the 1D blood flow model.

Science.gov (United States)

Wang, Xiao-Fei; Nishi, Shohei; Matsukawa, Mami; Ghigo, Arthur; Lagrée, Pierre-Yves; Fullana, Jose-Maria

2016-02-29

We study the behavior of the pulse waves of water into a flexible tube for application to blood flow simulations. In pulse waves both fluid friction and wall viscosity are damping factors, and difficult to evaluate separately. In this paper, the coefficients of fluid friction and wall viscosity are estimated by fitting a nonlinear 1D flow model to experimental data. In the experimental setup, a distensible tube is connected to a piston pump at one end and closed at another end. The pressure and wall displacements are measured simultaneously. A good agreement between model predictions and experiments was achieved. For amplitude decrease, the effect of wall viscosity on the pulse wave has been shown as important as that of fluid viscosity. Copyright © 2016 Elsevier Ltd. All rights reserved.

Bifurcated SEN with Fluid Flow Conditioners

Directory of Open Access Journals (Sweden)

F. Rivera-Perez

2014-01-01

Full Text Available This work evaluates the performance of a novel design for a bifurcated submerged entry nozzle (SEN used for the continuous casting of steel slabs. The proposed design incorporates fluid flow conditioners attached on SEN external wall. The fluid flow conditioners impose a pseudosymmetric pattern in the upper zone of the mold by inhibiting the fluid exchange between the zones created by conditioners. The performance of the SEN with fluid flow conditioners is analyzed through numerical simulations using the CFD technique. Numerical results were validated by means of physical simulations conducted on a scaled cold water model. Numerical and physical simulations confirmed that the performance of the proposed SEN is superior to a traditional one. Fluid flow conditioners reduce the liquid free surface fluctuations and minimize the occurrence of vortexes at the free surface.

Analysis and Modeling of Structure Formation in Granular and Fluid-Solid Flows

Science.gov (United States)

Murphy, Eric

Granular and multiphase flows are encountered in a number of industrial processes with particular emphasis in this manuscript given to the particular applications in cement pumping, pneumatic conveying, fluid catalytic cracking, CO2 capture, and fast pyrolysis of bio-materials. These processes are often modeled using averaged equations that may be simulated using computational fluid dynamics. Closure models are then required that describe the average forces that arise from both interparticle interactions, e.g. shear stress, and interphase interactions, such as mean drag. One of the biggest hurdles to this approach is the emergence of non-trivial spatio-temporal structures in the particulate phase, which can significantly modify the qualitative behavior of these forces and the resultant flow phenomenology. For example, the formation of large clusters in cohesive granular flows is responsible for a transition from solid-like to fluid-like rheology. Another example is found in gas-solid systems, where clustering at small scales is observed to significantly lower in the observed drag. Moreover, there remains the possibility that structure formation may occur at all scales, leading to a lack of scale separation required for traditional averaging approaches. In this context, several modeling problems are treated 1) first-principles based modeling of the rheology of cement slurries, 2) modeling the mean solid-solid drag experienced by polydisperse particles undergoing segregation, and 3) modeling clustering in homogeneous gas-solid flows. The first and third components are described in greater detail. In the study on the rheology of cements, several sub-problems are introduced, which systematically increase in the number and complexity of interparticle interactions. These interparticle interactions include inelasticity, friction, cohesion, and fluid interactions. In the first study, the interactions between cohesive inelastic particles was fully characterized for the

Analysis of the laminar Newtonian fluid flow through a thin fracture modelled as a fluid-saturated sparsely packed porous medium

Energy Technology Data Exchange (ETDEWEB)

Pazanin, Igor [Zagreb Univ. (Croatia). Dept. of Mathematics; Siddheshwar, Pradeep G. [Bangalore Univ., Bengaluru (India). Dept. of Mathematics

2017-06-01

In this article we investigate the fluid flow through a thin fracture modelled as a fluid-saturated porous medium. We assume that the fracture has constrictions and that the flow is governed by the prescribed pressure drop between the edges of the fracture. The problem is described by the Darcy-Lapwood-Brinkman model acknowledging the Brinkman extension of the Darcy law as well as the flow inertia. Using asymptotic analysis with respect to the thickness of the fracture, we derive the explicit higher-order approximation for the velocity distribution. We make an error analysis to comment on the order of accuracy of the method used and also to provide rigorous justification for the model.

Non-Newtonian fluid flow in 2D fracture networks

Science.gov (United States)

Zou, L.; Håkansson, U.; Cvetkovic, V.

2017-12-01

Modeling of non-Newtonian fluid (e.g., drilling fluids and cement grouts) flow in fractured rocks is of interest in many geophysical and industrial practices, such as drilling operations, enhanced oil recovery and rock grouting. In fractured rock masses, the flow paths are dominated by fractures, which are often represented as discrete fracture networks (DFN). In the literature, many studies have been devoted to Newtonian fluid (e.g., groundwater) flow in fractured rock using the DFN concept, but few works are dedicated to non-Newtonian fluids.In this study, a generalized flow equation for common non-Newtonian fluids (such as Bingham, power-law and Herschel-Bulkley) in a single fracture is obtained from the analytical solutions for non-Newtonian fluid discharge between smooth parallel plates. Using Monte Carlo sampling based on site characterization data for the distribution of geometrical features (e.g., density, length, aperture and orientations) in crystalline fractured rock, a two dimensional (2D) DFN model is constructed for generic flow simulations. Due to complex properties of non-Newtonian fluids, the relationship between fluid discharge and the pressure gradient is nonlinear. A Galerkin finite element method solver is developed to iteratively solve the obtained nonlinear governing equations for the 2D DFN model. Using DFN realizations, simulation results for different geometrical distributions of the fracture network and different non-Newtonian fluid properties are presented to illustrate the spatial discharge distributions. The impact of geometrical structures and the fluid properties on the non-Newtonian fluid flow in 2D DFN is examined statistically. The results generally show that modeling non-Newtonian fluid flow in fractured rock as a DFN is feasible, and that the discharge distribution may be significantly affected by the geometrical structures as well as by the fluid constitutive properties.

Two-phase flow modeling for low concentration spherical particle motion through a Newtonian fluid

CSIR Research Space (South Africa)

Smit GJF

2010-11-01

Full Text Available the necessity to model the discrete nature of sep- cite this article in press as: G.J.F. Smit et al., Two-phase flow modeling for low concentration spherical particle motion through a ian fluid, Appl. Math. Comput. (2010), doi:10.1016/j.amc.2010.07.055 2... and Ribberin large-scale and long term morphologica Please cite this article in press as: G.J.F. Smit Newtonian fluid, Appl. Math. Comput. (2010), � 2010 Elsevier Inc. All rights reserved. modeling of multiphase flow has increasingly become the subject...

Experimental and modeling hydraulic studies of foam drilling fluid flowing through vertical smooth pipes

Directory of Open Access Journals (Sweden)

Amit Saxena

2017-06-01

Full Text Available Foam has emerged as an efficient drilling fluid for the drilling of low pressure, fractured and matured reservoirs because of its the ability to reduce formation damage, fluid loss, differential sticking etc. However the compressible nature along with its complicated rheology has made its implementation a multifaceted task. Knowledge of the hydrodynamic behavior of drilling fluid within the borehole is the key behind successful implementation of drilling job. However, little effort has been made to develop the hydrodynamic models for the foam flowing with cuttings through pipes of variable diameter. In the present study, hydrodynamics of the foam fluid was investigated through the vertical smooth pipes of different pipe diameters, with variable foam properties in a flow loop system. Effect of cutting loading on pressure drop was also studied. Thus, the present investigation estimates the differential pressure loss across the pipe. The flow loop permits foam flow through 25.4 mm, 38.1 mm and 50.8 mm diameter pipes. The smaller diameter pipes are used to replicate the annular spaces between the drill string and wellbore. The developed model determines the pressure loss along the pipe and the results are compared with a number of existing models. The developed model is able to predict the experimental results more accurately.

Unified fluid flow model for pressure transient analysis in naturally fractured media

International Nuclear Information System (INIS)

Babak, Petro; Azaiez, Jalel

2015-01-01

Naturally fractured reservoirs present special challenges for flow modeling with regards to their internal geometrical structure. The shape and distribution of matrix porous blocks and the geometry of fractures play key roles in the formulation of transient interporosity flow models. Although these models have been formulated for several typical geometries of the fracture networks, they appeared to be very dissimilar for different shapes of matrix blocks, and their analysis presents many technical challenges. The aim of this paper is to derive and analyze a unified approach to transient interporosity flow models for slightly compressible fluids that can be used for any matrix geometry and fracture network. A unified fractional differential transient interporosity flow model is derived using asymptotic analysis for singularly perturbed problems with small parameters arising from the assumption of a much smaller permeability of the matrix blocks compared to that of the fractures. This methodology allowed us to unify existing transient interporosity flow models formulated for different shapes of matrix blocks including bounded matrix blocks, unbounded matrix cylinders with any orthogonal crossection, and matrix slabs. The model is formulated using a fractional order diffusion equation for fluid pressure that involves Caputo derivative of order 1/2 with respect to time. Analysis of the unified fractional derivative model revealed that the surface area-to-volume ratio is the key parameter in the description of the flow through naturally fractured media. Expressions of this parameter are presented for matrix blocks of the same geometrical shape as well as combinations of different shapes with constant and random sizes. Numerical comparisons between the predictions of the unified model and those obtained from existing transient interporosity ones for matrix blocks in the form of slabs, spheres and cylinders are presented for linear, radial and spherical flow types for

Experiments and Simulations of Fluid Flow in Heterogeneous Reservoir Models - Emphasis on Impacts from Crossbeds and Fractures

Energy Technology Data Exchange (ETDEWEB)

Boerresen, Knut Arne

1996-12-31

Hydrocarbon recovery from subsurface reservoirs has become increasingly dependent on advanced recovery techniques that require improved understanding of the physics of fluid flow within and across geological units including small-scale heterogeneities and fractures. In this thesis, impacts from heterogeneities on local fluid flow are studied experimentally by means of imaging techniques to visualize fluid flow in two dimensions during flooding of larger reservoir models. Part 1 reflects the multi-disciplinary collaboration, by briefly introducing the relevant geology, the literature on experiments on fluid flow in bedded structures, and outlining the applied numerical simulator and imaging techniques applied to visualize fluid flow. The second part contains a synopsis of displacement experiments in naturally laminated sandstones and in crossbed laboratory models, and of the impact from incipient shear fractures on oil recovery. The detailed results obtained from the experiments and simulations are described in six papers, all included. 215 refs., 108 figs., 16 tabs.

Experiments and Simulations of Fluid Flow in Heterogeneous Reservoir Models - Emphasis on Impacts from Crossbeds and Fractures

Energy Technology Data Exchange (ETDEWEB)

Boerresen, Knut Arne

1997-12-31

Hydrocarbon recovery from subsurface reservoirs has become increasingly dependent on advanced recovery techniques that require improved understanding of the physics of fluid flow within and across geological units including small-scale heterogeneities and fractures. In this thesis, impacts from heterogeneities on local fluid flow are studied experimentally by means of imaging techniques to visualize fluid flow in two dimensions during flooding of larger reservoir models. Part 1 reflects the multi-disciplinary collaboration, by briefly introducing the relevant geology, the literature on experiments on fluid flow in bedded structures, and outlining the applied numerical simulator and imaging techniques applied to visualize fluid flow. The second part contains a synopsis of displacement experiments in naturally laminated sandstones and in crossbed laboratory models, and of the impact from incipient shear fractures on oil recovery. The detailed results obtained from the experiments and simulations are described in six papers, all included. 215 refs., 108 figs., 16 tabs.

Computational Fluid Dynamics Analysis of Pulsatile Blood Flow Behavior in Modelled Stenosed Vessels with Different Severities

Directory of Open Access Journals (Sweden)

Mohsen Mehrabi

2012-01-01

Full Text Available This study focuses on the behavior of blood flow in the stenosed vessels. Blood is modelled as an incompressible non-Newtonian fluid which is based on the power law viscosity model. A numerical technique based on the finite difference method is developed to simulate the blood flow taking into account the transient periodic behaviour of the blood flow in cardiac cycles. Also, pulsatile blood flow in the stenosed vessel is based on the Womersley model, and fluid flow in the lumen region is governed by the continuity equation and the Navier-Stokes equations. In this study, the stenosis shape is cosine by using Tu and Devil model. Comparing the results obtained from three stenosed vessels with 30%, 50%, and 75% area severity, we find that higher percent-area severity of stenosis leads to higher extrapressure jumps and higher blood speeds around the stenosis site. Also, we observe that the size of the stenosis in stenosed vessels does influence the blood flow. A little change on the cross-sectional value makes vast change on the blood flow rate. This simulation helps the people working in the field of physiological fluid dynamics as well as the medical practitioners.

Advances in fluid modeling and turbulence measurements

International Nuclear Information System (INIS)

Wada, Akira; Ninokata, Hisashi; Tanaka, Nobukazu

2002-01-01

The context of this book consists of four fields: Environmental Fluid Mechanics; Industrial Fluid Mechanics; Fundamentals of Fluid Mechanics; and Turbulence Measurements. Environmental Fluid Mechanics includes free surface flows in channels, rivers, seas, and estuaries. It also discusses wind engineering issues, ocean circulation model and dispersion problems in atmospheric, water and ground water environments. In Industrial Fluid Mechanics, fluid phenomena in energy exchanges, modeling of turbulent two- or multi-phase flows, swirling flows, flows in combustors, variable density flows and reacting flows, flows in turbo-machines, pumps and piping systems, and fluid-structure interaction are discussed. In Fundamentals of Fluid Mechanics, progress in modeling turbulent flows and heat/mass transfers, computational fluid dynamics/numerical techniques, parallel computing algorithms, applications of chaos/fractal theory in turbulence are reported. In Turbulence Measurements, experimental studies of turbulent flows, experimental and post-processing techniques, quantitative and qualitative flow visualization techniques are discussed. Separate abstracts were presented for 15 of the papers in this issue. The remaining 89 were considered outside the subject scope of INIS. (J.P.N.)

Leaback of Pulsatile Flow of Particle Fluid Suspension Model of ...

African Journals Online (AJOL)

Leaback of Pulsatile Flow of Particle Fluid Suspension Model of Blood Under Periodic Body Acceleration. ... The variation in body acceleration amplitude though affects the velocity profile in the capillary tubes, it has no effect on the leakback in the tubes. Leakback is mainly determined by the balance of the viscous drag and ...

Fluid simulation of tokamak ion temperature gradient turbulence with zonal flow closure model

Energy Technology Data Exchange (ETDEWEB)

Yamagishi, Osamu, E-mail: yamagisi@nifs.ac.jp; Sugama, Hideo [National Institute for Fusion Science, Toki, Gifu 509-5292 (Japan)

2016-03-15

Nonlinear fluid simulation of turbulence driven by ion temperature gradient modes in the tokamak fluxtube configuration is performed by combining two different closure models. One model is a gyrofluid model by Beer and Hammett [Phys. Plasmas 3, 4046 (1996)], and the other is a closure model to reproduce the kinetic zonal flow response [Sugama et al., Phys. Plasmas 14, 022502 (2007)]. By including the zonal flow closure, generation of zonal flows, significant reduction in energy transport, reproduction of the gyrokinetic transport level, and nonlinear upshift on the critical value of gradient scale length are observed.

Fluid simulation of tokamak ion temperature gradient turbulence with zonal flow closure model

Science.gov (United States)

Yamagishi, Osamu; Sugama, Hideo

2016-03-01

Nonlinear fluid simulation of turbulence driven by ion temperature gradient modes in the tokamak fluxtube configuration is performed by combining two different closure models. One model is a gyrofluid model by Beer and Hammett [Phys. Plasmas 3, 4046 (1996)], and the other is a closure model to reproduce the kinetic zonal flow response [Sugama et al., Phys. Plasmas 14, 022502 (2007)]. By including the zonal flow closure, generation of zonal flows, significant reduction in energy transport, reproduction of the gyrokinetic transport level, and nonlinear upshift on the critical value of gradient scale length are observed.

Multi-fluid modelling of pulsed discharges for flow control applications

Science.gov (United States)

Poggie, J.

2015-02-01

Experimental evidence suggests that short-pulse dielectric barrier discharge actuators are effective for speeds corresponding to take-off and approach of large aircraft, and thus are a fruitful direction for flow control technology development. Large-eddy simulations have reproduced some of the main fluid dynamic effects. The plasma models used in such simulations are semi-empirical, however, and need to be tuned for each flowfield under consideration. In this paper, the discharge physics is examined in more detail with multi-fluid modelling, comparing a five-moment model (continuity, momentum, and energy equations) to a two-moment model (continuity and energy equations). A steady-state, one-dimensional discharge was considered first, and the five-moment model was found to predict significantly lower ionisation rates and number densities than the two-moment model. A two-dimensional, transient discharge problem with an elliptical cathode was studied next. Relative to the two-moment model, the five-moment model predicted a slower response to the activation of the cathode, and lower electron velocities and temperatures as the simulation approached steady-state. The primary reason for the differences in the predictions of the two models can be attributed to the effects of particle inertia, particularly electron inertia in the cathode layer. The computational cost of the five-moment model is only about twice that of the simpler variant, suggesting that it may be feasible to use the more sophisticated model in practical calculations for flow control actuator design.

Contribution to the modeling of particulate hypersonic flows. Study and validation of a discrete two-fluid model

International Nuclear Information System (INIS)

Papin, M.

2005-06-01

This work dedicated to the study of the hypersonic re-entry of vehicles in the atmosphere crossing clouds of particles implies the study of two-fluid flow and it is shown that some developments can be applied to the two-fluid models used to describe the phase transformation occurring in a target irradiated by laser beams. The calculation of wall fluxes on hypersonic re-entry vehicles requires the modeling of the interactions with clouds. Two-fluid flows posing many physical and mathematical problems, one studies an alternative model due to Abgrall and Saurel: the discrete equation method (DEM). Three axis are chosen. The first proposes a finite volume discretization of the Navier-Stokes equations on hybrid grids adapted to the context. The second extends the DEM within a multi-fluid not-structured N-D framework. A limit study associates an original continuous model to him: it allows to modify usual two-fluid seven equations models to obtain a phasic entropy principle. In spite of good properties, the continuous description of the particles is unsuited to the problem. The last axis is a study of the follow-up of pointwise particles which does not allow realistic calculation of parietal fluxes. An original model, extending the usual hydro-erosion models, however makes it possible to evaluate rebounds, erosion of the body and wall fluxes. The appendices expose approximate and exact Riemann solvers between pure fluids, discretization of the Baer and Nunziato model, and relations describing the atmosphere, water and heat fluxes

Modeling of flow of particles in a non-Newtonian fluid using lattice Boltzmann method

DEFF Research Database (Denmark)

Skocek, Jan; Svec, Oldrich; Spangenberg, Jon

2011-01-01

is necessary. In this contribution, the model at the scale of aggregates is introduced. The conventional lattice Boltzmann method for fluid flow is enriched with the immersed boundary method with direct forcing to simulate the flow of rigid particles in a non- Newtonian liquid. Basic ingredients of the model...

Validation of a multidimensional computational fluid dynamics model for subcooled flow boiling analysis

Energy Technology Data Exchange (ETDEWEB)

Braz Filho, Francisco A.; Caldeira, Alexandre D.; Borges, Eduardo M., E-mail: fbraz@ieav.cta.b, E-mail: alexdc@ieav.cta.b, E-mail: eduardo@ieav.cta.b [Instituto de Estudos Avancados (IEAv/CTA), Sao Jose dos Campos, SP (Brazil). Div. de Energia Nuclear

2011-07-01

In a heated vertical channel, the subcooled flow boiling regime occurs when the bulk fluid temperature is lower than the saturation temperature, but the fluid temperature reaches the saturation point near the channel wall. This phenomenon produces a significant increase in heat flux, limited by the critical heat flux. This study is particularly important to the thermal-hydraulics analysis of pressurized water reactors. The purpose of this work is the validation of a multidimensional model to analyze the subcooled flow boiling comparing the results with experimental data found in literature. The computational fluid dynamics code FLUENT was used with Eulerian multiphase model option. The calculated values of wall temperature in the liquid-solid interface presented an excellent agreement when compared to the experimental data. Void fraction calculations presented satisfactory results in relation to the experimental data in pressures of 15, 30 and 45 bars. (author)

Validation of a multidimensional computational fluid dynamics model for subcooled flow boiling analysis

International Nuclear Information System (INIS)

Braz Filho, Francisco A.; Caldeira, Alexandre D.; Borges, Eduardo M.

2011-01-01

In a heated vertical channel, the subcooled flow boiling regime occurs when the bulk fluid temperature is lower than the saturation temperature, but the fluid temperature reaches the saturation point near the channel wall. This phenomenon produces a significant increase in heat flux, limited by the critical heat flux. This study is particularly important to the thermal-hydraulics analysis of pressurized water reactors. The purpose of this work is the validation of a multidimensional model to analyze the subcooled flow boiling comparing the results with experimental data found in literature. The computational fluid dynamics code FLUENT was used with Eulerian multiphase model option. The calculated values of wall temperature in the liquid-solid interface presented an excellent agreement when compared to the experimental data. Void fraction calculations presented satisfactory results in relation to the experimental data in pressures of 15, 30 and 45 bars. (author)

Map of fluid flow in fractal porous medium into fractal continuum flow.

Science.gov (United States)

Balankin, Alexander S; Elizarraraz, Benjamin Espinoza

2012-05-01

This paper is devoted to fractal continuum hydrodynamics and its application to model fluid flows in fractally permeable reservoirs. Hydrodynamics of fractal continuum flow is developed on the basis of a self-consistent model of fractal continuum employing vector local fractional differential operators allied with the Hausdorff derivative. The generalized forms of Green-Gauss and Kelvin-Stokes theorems for fractional calculus are proved. The Hausdorff material derivative is defined and the form of Reynolds transport theorem for fractal continuum flow is obtained. The fundamental conservation laws for a fractal continuum flow are established. The Stokes law and the analog of Darcy's law for fractal continuum flow are suggested. The pressure-transient equation accounting the fractal metric of fractal continuum flow is derived. The generalization of the pressure-transient equation accounting the fractal topology of fractal continuum flow is proposed. The mapping of fluid flow in a fractally permeable medium into a fractal continuum flow is discussed. It is stated that the spectral dimension of the fractal continuum flow d(s) is equal to its mass fractal dimension D, even when the spectral dimension of the fractally porous or fissured medium is less than D. A comparison of the fractal continuum flow approach with other models of fluid flow in fractally permeable media and the experimental field data for reservoir tests are provided.

A Constitutive Model for Flow-Induced Anisotropic Behavior of Viscoelastic Complex Fluids

International Nuclear Information System (INIS)

Zhu, H.; De Kee, D.

2008-01-01

Flow-induced structural anisotropy could result when a complex fluid system is removed from equilibrium by means of hydrodynamic forces. In this paper, a general theory is developed to model flow induced anisotropic behavior of complex viscoelastic systems, e.g. polymer solutions/melts and suspensions. The rheological properties are characterized by viscosity and relaxation time tensors. We consider a second-rank tensor as a measure of the microstructure. We consider the effect of the flow on the structural changes: i.e. the evolution of the microstructure tensor is governed by a relaxation-type differential equation. We also propose that the viscosity and the relaxation time tensors depend on the second-rank microstructure tensor. That is as the microstructure tensor changes with the applied rate of deformation, the viscosity and relaxation time tensors evolve accordingly. As an example we consider elongational flow of two complex fluids

Optimization of a Two-Fluid Hydrodynamic Model of Churn-Turbulent Flow

Energy Technology Data Exchange (ETDEWEB)

Donna Post Guillen

2009-07-01

A hydrodynamic model of two-phase, churn-turbulent flows is being developed using the computational multiphase fluid dynamics (CMFD) code, NPHASE-CMFD. The numerical solutions obtained by this model are compared with experimental data obtained at the TOPFLOW facility of the Institute of Safety Research at the Forschungszentrum Dresden-Rossendorf. The TOPFLOW data is a high quality experimental database of upward, co-current air-water flows in a vertical pipe suitable for validation of computational fluid dynamics (CFD) codes. A five-field CMFD model was developed for the continuous liquid phase and four bubble size groups using mechanistic closure models for the ensemble-averaged Navier-Stokes equations. Mechanistic models for the drag and non-drag interfacial forces are implemented to include the governing physics to describe the hydrodynamic forces controlling the gas distribution. The closure models provide the functional form of the interfacial forces, with user defined coefficients to adjust the force magnitude. An optimization strategy was devised for these coefficients using commercial design optimization software. This paper demonstrates an approach to optimizing CMFD model parameters using a design optimization approach. Computed radial void fraction profiles predicted by the NPHASE-CMFD code are compared to experimental data for four bubble size groups.

Oscillating flow of a Burgers' fluid in a pipe

International Nuclear Information System (INIS)

Khan, M.; Asghar, S.; Hayat, T.

2005-12-01

An analysis is made to see the influences of Hall current on the flow of a Burgers' fluid. The velocity field corresponding to flow in a pipe is determined. The closed form analytical solutions for several Newtonian and non-Newtonian fluid models can be obtained from the present analysis as the limiting cases. The purpose of this work is twofold. Firstly, to investigate the oscillating flow in a pipe using Burgers? fluid model. Secondly, to see the effects of Hall current on the velocity field. The flow in a pipe is induced due to imposition of an oscillating pressure gradient. An exact analytical solution to the governing problem is given using the Fourier transform technique. The obtained expression for the velocity field shows that there are pronounced effects of Hall and rheological parameters. The considered fluid model is a viscoelastic model and has been used to characterize food products such as cheese, soil, asphalt and asphalt mixes etc. (author)

Steady laminar flow of fractal fluids

Energy Technology Data Exchange (ETDEWEB)

Balankin, Alexander S., E-mail: abalankin@ipn.mx [Grupo Mecánica Fractal, ESIME, Instituto Politécnico Nacional, México D.F., 07738 (Mexico); Mena, Baltasar [Laboratorio de Ingeniería y Procesos Costeros, Instituto de Ingeniería, Universidad Nacional Autónoma de México, Sisal, Yucatán, 97355 (Mexico); Susarrey, Orlando; Samayoa, Didier [Grupo Mecánica Fractal, ESIME, Instituto Politécnico Nacional, México D.F., 07738 (Mexico)

2017-02-12

We study laminar flow of a fractal fluid in a cylindrical tube. A flow of the fractal fluid is mapped into a homogeneous flow in a fractional dimensional space with metric induced by the fractal topology. The equations of motion for an incompressible Stokes flow of the Newtonian fractal fluid are derived. It is found that the radial distribution for the velocity in a steady Poiseuille flow of a fractal fluid is governed by the fractal metric of the flow, whereas the pressure distribution along the flow direction depends on the fractal topology of flow, as well as on the fractal metric. The radial distribution of the fractal fluid velocity in a steady Couette flow between two concentric cylinders is also derived. - Highlights: • Equations of Stokes flow of Newtonian fractal fluid are derived. • Pressure distribution in the Newtonian fractal fluid is derived. • Velocity distribution in Poiseuille flow of fractal fluid is found. • Velocity distribution in a steady Couette flow is established.

Analysis of two-phase flow inter-subchannel mass and momentum exchanges by the two-fluid model approach

Energy Technology Data Exchange (ETDEWEB)

Ninokata, H. [Tokyo Institute of Technology (Japan); Deguchi, A. [ENO Mathematical Analysis, Tokyo (Japan); Kawahara, A. [Kumamoto Univ., Kumamoto (Japan)

1995-09-01

A new void drift model for the subchannel analysis method is presented for the thermohydraulics calculation of two-phase flows in rod bundles where the flow model uses a two-fluid formulation for the conservation of mass, momentum and energy. A void drift model is constructed based on the experimental data obtained in a geometrically simple inter-connected two circular channel test sections using air-water as working fluids. The void drift force is assumed to be an origin of void drift velocity components of the two-phase cross-flow in a gap area between two adjacent rods and to overcome the momentum exchanges at the phase interface and wall-fluid interface. This void drift force is implemented in the cross flow momentum equations. Computational results have been successfully compared to experimental data available including 3x3 rod bundle data.

Williamson Fluid Model for the Peristaltic Flow of Chyme in Small Intestine

Directory of Open Access Journals (Sweden)

Sohail Nadeem

2012-01-01

Full Text Available Mathematical model for the peristaltic flow of chyme in small intestine along with inserted endoscope is considered. Here, chyme is treated as Williamson fluid, and the flow is considered between the annular region formed by two concentric tubes (i.e., outer tube as small intestine and inner tube as endoscope. Flow is induced by two sinusoidal peristaltic waves of different wave lengths, traveling down the intestinal wall with the same speed. The governing equations of Williamson fluid in cylindrical coordinates have been modeled. The resulting nonlinear momentum equations are simplified using long wavelength and low Reynolds number approximations. The resulting problem is solved using regular perturbation method in terms of a variant of Weissenberg number We. The numerical solution of the problem is also computed by using shooting method, and comparison of results of both solutions for velocity field is presented. The expressions for axial velocity, frictional force, pressure rise, stream function, and axial pressure gradient are obtained, and the effects of various emerging parameters on the flow characteristics are illustrated graphically. Furthermore, the streamlines pattern is plotted, and it is observed that trapping occurs, and the size of the trapped bolus varies with varying embedded flow parameters.

Simulation of dendritic growth of magnesium alloys with fluid flow

Directory of Open Access Journals (Sweden)

Meng-wu Wu

2017-11-01

Full Text Available Fluid flow has a significant impact on the microstructure evolution of alloys during solidification. Based on the previous work relating simulation of the dendritic growth of magnesium alloys with hcp (hexagonal close-packed structure, an extension was made to the formerly established CA (cellular automaton model with the purpose of studying the effect of fluid flow on the dendritic growth of magnesium alloys. The modified projection method was used to solve the transport equations of flow field. By coupling the flow field with the solute field, simulation results of equiaxed and columnar dendritic growth of magnesium alloys with fluid flow were achieved. The simulated results were quantitatively compared with those without fluid flow. Moreover, a comparison was also made between the present work and previous works conducted by others. It can be concluded that a deep understanding of the dendritic growth of magnesium alloys with fluid flow can be obtained by applying the present numerical model.

Insertable fluid flow passage bridgepiece and method

Science.gov (United States)

Jones, Daniel O.

2000-01-01

A fluid flow passage bridgepiece for insertion into an open-face fluid flow channel of a fluid flow plate is provided. The bridgepiece provides a sealed passage from a columnar fluid flow manifold to the flow channel, thereby preventing undesirable leakage into and out of the columnar fluid flow manifold. When deployed in the various fluid flow plates that are used in a Proton Exchange Membrane (PEM) fuel cell, bridgepieces of this invention prevent mixing of reactant gases, leakage of coolant or humidification water, and occlusion of the fluid flow channel by gasket material. The invention also provides a fluid flow plate assembly including an insertable bridgepiece, a fluid flow plate adapted for use with an insertable bridgepiece, and a method of manufacturing a fluid flow plate with an insertable fluid flow passage bridgepiece.

Poiseuille equation for steady flow of fractal fluid

Science.gov (United States)

Tarasov, Vasily E.

2016-07-01

Fractal fluid is considered in the framework of continuous models with noninteger dimensional spaces (NIDS). A recently proposed vector calculus in NIDS is used to get a description of fractal fluid flow in pipes with circular cross-sections. The Navier-Stokes equations of fractal incompressible viscous fluids are used to derive a generalization of the Poiseuille equation of steady flow of fractal media in pipe.

Thermo-Fluid Dynamics of Two-Phase Flow

CERN Document Server

Ishii, Mamrou

2011-01-01

"Thermo-fluid Dynamics of Two-Phase Flow, Second Edition" is focused on the fundamental physics of two-phase flow. The authors present the detailed theoretical foundation of multi-phase flow thermo-fluid dynamics as they apply to: Nuclear reactor transient and accident analysis; Energy systems; Power generation systems; Chemical reactors and process systems; Space propulsion; Transport processes. This edition features updates on two-phase flow formulation and constitutive equations and CFD simulation codes such as FLUENT and CFX, new coverage of the lift force model, which is of part

3D Printing of Fluid Flow Structures

OpenAIRE

Taira, Kunihiko; Sun, Yiyang; Canuto, Daniel

2017-01-01

We discuss the use of 3D printing to physically visualize (materialize) fluid flow structures. Such 3D models can serve as a refreshing hands-on means to gain deeper physical insights into the formation of complex coherent structures in fluid flows. In this short paper, we present a general procedure for taking 3D flow field data and producing a file format that can be supplied to a 3D printer, with two examples of 3D printed flow structures. A sample code to perform this process is also prov...

a Matlab Toolbox for Basin Scale Fluid Flow Modeling Applied to Hydrology and Geothermal Energy

Science.gov (United States)

Alcanie, M.; Lupi, M.; Carrier, A.

2017-12-01

Recent boosts in the development of geothermal energy were fostered by the latest oil crises and by the need of reducing CO2 emissions generated by the combustion of fossil fuels. Various numerical codes (e.g. FEHM, CSMP++, HYDROTHERM, TOUGH) have thus been implemented for the simulation and quantification of fluid flow in the upper crust. One possible limitation of such codes is the limited accessibility and the complex structure of the simulators. For this reason, we began to develop a Hydrothermal Fluid Flow Matlab library as part of MRST (Matlab Reservoir Simulation Toolbox). MRST is designed for the simulation of oil and gas problems including carbon capture storage. However, a geothermal module is still missing. We selected the Geneva Basin as a natural laboratory because of the large amount of data available in the region. The Geneva Basin has been intensely investigated in the past with exploration wells, active seismic and gravity surveys. In addition, the energy strategy of Switzerland promotes the development of geothermal energy that lead to recent geophysical prospections. Previous and ongoing projects have shown the geothermal potential of the Geneva Basin but a consistent fluid flow model assessing the deep circulation in the region is yet to be defined. The first step of the study was to create the basin-scale static model. We integrated available active seismic, gravity inversions and borehole data to describe the principal geologic and tectonic features of the Geneva Basin. Petrophysical parameters were obtained from available and widespread well logs. This required adapting MRST to standard text format file imports and outline a new methodology for quick static model creation in an open source environment. We implemented several basin-scale fluid flow models to test the effects of petrophysical properties on the circulation dynamics of deep fluids in the Geneva Basin. Preliminary results allow the identification of preferential fluid flow

Evaluation of Interfacial Heat Transfer Models for Flashing Flow with Two-Fluid CFD

Directory of Open Access Journals (Sweden)

Yixiang Liao

2018-06-01

Full Text Available The complexity of flashing flows is increased vastly by the interphase heat transfer as well as its coupling with mass and momentum transfers. A reliable heat transfer coefficient is the key in the modelling of such kinds of flows with the two-fluid model. An extensive literature survey on computational modelling of flashing flows has been given in previous work. The present work is aimed at giving a brief review on available theories and correlations for the estimation of interphase heat transfer coefficient, and evaluating them quantitatively based on computational fluid dynamics simulations of bubble growth in superheated liquid. The comparison of predictions for bubble growth rate obtained by using different correlations with the experimental as well as direct numerical simulation data reveals that the performance of the correlations is dependent on the Jakob number and Reynolds number. No generally applicable correlations are available. Both conduction and convection are important in cases of bubble rising and translating in stagnant liquid at high Jakob numbers. The correlations combining the analytical solution for heat diffusion and the theoretical relation for potential flow give the best agreement.

Modeling studies of multiphase fluid and heat flow processes in nuclear waste isolation

International Nuclear Information System (INIS)

Pruess, K.

1989-01-01

Multiphase fluid and heat flow plays an important role in many problems relating to the disposal of nuclear wastes in geologic media. Examples include boiling and condensation processes near heat-generating wastes, flow of water and formation gas in partially saturated formations, evolution of a free gas phase from waste package corrosion in initially water-saturated environments, and redistribution (dissolution, transport and precipitation) of rock minerals in non-isothermal flow fields. Such processes may strongly impact upon waste package and repository design considerations and performance. This paper summarizes important physical phenomena occurring in multiphase and nonisothermal flows, as well as techniques for their mathematical modeling and numerical simulation. Illustrative applications are given for a number of specific fluid and heat flow problems, including: thermohydrologic conditions near heat-generating waste packages in the unsaturated zone; repositorywide convection effects in the unsaturated zone; effects of quartz dissolution and precipitation for disposal in the saturated zone; and gas pressurization and flow effects from corrosion of low-level waste packages

Acoustic concentration of particles in fluid flow

Science.gov (United States)

Ward, Michael D.; Kaduchak, Gregory

2010-11-23

An apparatus for acoustic concentration of particles in a fluid flow includes a substantially acoustically transparent membrane and a vibration generator that define a fluid flow path therebetween. The fluid flow path is in fluid communication with a fluid source and a fluid outlet and the vibration generator is disposed adjacent the fluid flow path and is capable of producing an acoustic field in the fluid flow path. The acoustic field produces at least one pressure minima in the fluid flow path at a predetermined location within the fluid flow path and forces predetermined particles in the fluid flow path to the at least one pressure minima.

Acoustic concentration of particles in fluid flow

Science.gov (United States)

Ward, Michael W.; Kaduchak, Gregory

2017-08-15

Disclosed herein is a acoustic concentration of particles in a fluid flow that includes a substantially acoustically transparent membrane and a vibration generator that define a fluid flow path therebetween. The fluid flow path is in fluid communication with a fluid source and a fluid outlet and the vibration generator is disposed adjacent the fluid flow path and is capable of producing an acoustic field in the fluid flow path. The acoustic field produces at least one pressure minima in the fluid flow path at a predetermined location within the fluid flow path and forces predetermined particles in the fluid flow path to the at least one pressure minima.

A numerical model of two-phase flow at the micro-scale using the volume-of-fluid method

Science.gov (United States)

Shams, Mosayeb; Raeini, Ali Q.; Blunt, Martin J.; Bijeljic, Branko

2018-03-01

This study presents a simple and robust numerical scheme to model two-phase flow in porous media where capillary forces dominate over viscous effects. The volume-of-fluid method is employed to capture the fluid-fluid interface whose dynamics is explicitly described based on a finite volume discretization of the Navier-Stokes equations. Interfacial forces are calculated directly on reconstructed interface elements such that the total curvature is preserved. The computed interfacial forces are explicitly added to the Navier-Stokes equations using a sharp formulation which effectively eliminates spurious currents. The stability and accuracy of the implemented scheme is validated on several two- and three-dimensional test cases, which indicate the capability of the method to model two-phase flow processes at the micro-scale. In particular we show how the co-current flow of two viscous fluids leads to greatly enhanced flow conductance for the wetting phase in corners of the pore space, compared to a case where the non-wetting phase is an inviscid gas.

Simulation of uncompressible fluid flow through a porous media

International Nuclear Information System (INIS)

Ramirez, A.; Gonzalez, J.L.; Carrillo, F.; Lopez, S.

2009-01-01

Recently, a great interest has been focused for investigations about transport phenomena in disordered systems. One of the most treated topics is fluid flow through anisotropic materials due to the importance in many industrial processes like fluid flow in filters, membranes, walls, oil reservoirs, etc. In this work is described the formulation of a 2D mathematical model to simulate the fluid flow behavior through a porous media (PM) based on the solution of the continuity equation as a function of the Darcy's law for a percolation system; which was reproduced using computational techniques reproduced using a random distribution of the porous media properties (porosity, permeability and saturation). The model displays the filling of a partially saturated porous media with a new injected fluid showing the non-defined advance front and dispersion of fluids phenomena.

Dynamic analysis of electro- and magneto-rheological fluid dampers using duct flow models

International Nuclear Information System (INIS)

Esteki, Kambiz; Bagchi, Ashutosh; Sedaghati, Ramin

2014-01-01

Magneto-rheological (MR) and electro-rheological (ER) fluid dampers provide a semi-active control mechanism for suppressing vibration responses of a structure. MR and ER fluids change their viscosity under the influence of magnetic and electrical fields, respectively, which facilitates automatic control when these fluids are used in damping devices. The existing models, namely the phenomenological models for simulating the behavior of MR and ER dampers, rely on various parameters determined experimentally by the manufacturers for each damper configuration. It is of interest to develop mechanistic models of these dampers which can be applied to various configurations so that their fundamental characteristics can be studied to develop flexible design solutions for smart structures. This paper presents a formulation for dynamic analysis of electro-rheological (ER) and magneto-rheological (MR) fluid dampers in flow and mix mode configurations under harmonic and random excitations. The procedure employs the vorticity transport equation and the regularization function to deal with the unsteady flow and nonlinear behavior of ER/MR fluid in general motion. The finite difference method has been used to solve the governing differential equations. Using the developed approach, the damping force of ER/MR dampers can be calculated under any type of excitation. (paper)

Conjugate Compressible Fluid Flow and Heat Transfer in Ducts

Science.gov (United States)

Cross, M. F.

2011-01-01

A computational approach to modeling transient, compressible fluid flow with heat transfer in long, narrow ducts is presented. The primary application of the model is for analyzing fluid flow and heat transfer in solid propellant rocket motor nozzle joints during motor start-up, but the approach is relevant to a wide range of analyses involving rapid pressurization and filling of ducts. Fluid flow is modeled through solution of the spatially one-dimensional, transient Euler equations. Source terms are included in the governing equations to account for the effects of wall friction and heat transfer. The equation solver is fully-implicit, thus providing greater flexibility than an explicit solver. This approach allows for resolution of pressure wave effects on the flow as well as for fast calculation of the steady-state solution when a quasi-steady approach is sufficient. Solution of the one-dimensional Euler equations with source terms significantly reduces computational run times compared to general purpose computational fluid dynamics packages solving the Navier-Stokes equations with resolved boundary layers. In addition, conjugate heat transfer is more readily implemented using the approach described in this paper than with most general purpose computational fluid dynamics packages. The compressible flow code has been integrated with a transient heat transfer solver to analyze heat transfer between the fluid and surrounding structure. Conjugate fluid flow and heat transfer solutions are presented. The author is unaware of any previous work available in the open literature which uses the same approach described in this paper.

Robust and general method for determining surface fluid flow boundary conditions in articular cartilage contact mechanics modeling.

Science.gov (United States)

Pawaskar, Sainath Shrikant; Fisher, John; Jin, Zhongmin

2010-03-01

Contact detection in cartilage contact mechanics is an important feature of any analytical or computational modeling investigation when the biphasic nature of cartilage and the corresponding tribology are taken into account. The fluid flow boundary conditions will change based on whether the surface is in contact or not, which will affect the interstitial fluid pressurization. This in turn will increase or decrease the load sustained by the fluid phase, with a direct effect on friction, wear, and lubrication. In laboratory experiments or clinical hemiarthroplasty, when a rigid indenter or metallic prosthesis is used to apply load to the cartilage, there will not be any fluid flow normal to the surface in the contact region due to the impermeable nature of the indenter/prosthesis. In the natural joint, on the other hand, where two cartilage surfaces interact, flow will depend on the pressure difference across the interface. Furthermore, in both these cases, the fluid would flow freely in non-contacting regions. However, it should be pointed out that the contact area is generally unknown in advance in both cases and can only be determined as part of the solution. In the present finite element study, a general and robust algorithm was proposed to decide nodes in contact on the cartilage surface and, accordingly, impose the fluid flow boundary conditions. The algorithm was first tested for a rigid indenter against cartilage model. The algorithm worked well for two-dimensional four-noded and eight-noded axisymmetric element models as well as three-dimensional models. It was then extended to include two cartilages in contact. The results were in excellent agreement with the previous studies reported in the literature.

Fluid dynamics of bubbly flows

International Nuclear Information System (INIS)

Ziegenhein, Thomas

2016-01-01

Bubbly flows can be found in many applications in chemical, biological and power engineering. Reliable simulation tools of such flows that allow the design of new processes and optimization of existing one are therefore highly desirable. CFD-simulations applying the multi-fluid approach are very promising to provide such a design tool for complete facilities. In the multi-fluid approach, however, closure models have to be formulated to model the interaction between the continuous and dispersed phase. Due to the complex nature of bubbly flows, different phenomena have to be taken into account and for every phenomenon different closure models exist. Therefore, reliable predictions of unknown bubbly flows are not yet possible with the multi-fluid approach. A strategy to overcome this problem is to define a baseline model in which the closure models including the model constants are fixed so that the limitations of the modeling can be evaluated by validating it on different experiments. Afterwards, the shortcomings are identified so that the baseline model can be stepwise improved without losing the validity for the already validated cases. This development of a baseline model is done in the present work by validating the baseline model developed at the Helmholtz-Zentrum Dresden-Rossendorf mainly basing on experimental data for bubbly pipe flows to bubble columns, bubble plumes and air-lift reactors that are relevant in chemical and biological engineering applications. In the present work, a large variety of such setups is used for validation. The buoyancy driven bubbly flows showed thereby a transient behavior on the scale of the facility. Since such large scales are characterized by the geometry of the facility, turbulence models cannot describe them. Therefore, the transient simulation of bubbly flows with two equation models based on the unsteady Reynolds-averaged Navier-Stokes equations is investigated. In combination with the before mentioned baseline model these

Fluid dynamics of bubbly flows

Energy Technology Data Exchange (ETDEWEB)

Ziegenhein, Thomas

2016-07-08

Bubbly flows can be found in many applications in chemical, biological and power engineering. Reliable simulation tools of such flows that allow the design of new processes and optimization of existing one are therefore highly desirable. CFD-simulations applying the multi-fluid approach are very promising to provide such a design tool for complete facilities. In the multi-fluid approach, however, closure models have to be formulated to model the interaction between the continuous and dispersed phase. Due to the complex nature of bubbly flows, different phenomena have to be taken into account and for every phenomenon different closure models exist. Therefore, reliable predictions of unknown bubbly flows are not yet possible with the multi-fluid approach. A strategy to overcome this problem is to define a baseline model in which the closure models including the model constants are fixed so that the limitations of the modeling can be evaluated by validating it on different experiments. Afterwards, the shortcomings are identified so that the baseline model can be stepwise improved without losing the validity for the already validated cases. This development of a baseline model is done in the present work by validating the baseline model developed at the Helmholtz-Zentrum Dresden-Rossendorf mainly basing on experimental data for bubbly pipe flows to bubble columns, bubble plumes and air-lift reactors that are relevant in chemical and biological engineering applications. In the present work, a large variety of such setups is used for validation. The buoyancy driven bubbly flows showed thereby a transient behavior on the scale of the facility. Since such large scales are characterized by the geometry of the facility, turbulence models cannot describe them. Therefore, the transient simulation of bubbly flows with two equation models based on the unsteady Reynolds-averaged Navier-Stokes equations is investigated. In combination with the before mentioned baseline model these

Viscous Flow with Large Fluid-Fluid Interface Displacement

DEFF Research Database (Denmark)

Rasmussen, Henrik Koblitz; Hassager, Ole; Saasen, Arild

1998-01-01

The arbitrary Lagrange-Euler (ALE) kinematic description has been implemented in a 3D transient finite element program to simulate multiple fluid flows with fluid-fluid interface or surface displacements. The description of fluid interfaces includes variable interfacial tension, and the formulation...... is useful in the simulation of low and intermediate Reynolds number viscous flow. The displacement of two immiscible Newtonian fluids in a vertical (concentric and eccentric) annulus and a (vertical and inclined)tube is simulated....

Partitioned fluid-solid coupling for cardiovascular blood flow: left-ventricular fluid mechanics.

Science.gov (United States)

Krittian, Sebastian; Janoske, Uwe; Oertel, Herbert; Böhlke, Thomas

2010-04-01

We present a 3D code-coupling approach which has been specialized towards cardiovascular blood flow. For the first time, the prescribed geometry movement of the cardiovascular flow model KaHMo (Karlsruhe Heart Model) has been replaced by a myocardial composite model. Deformation is driven by fluid forces and myocardial response, i.e., both its contractile and constitutive behavior. Whereas the arbitrary Lagrangian-Eulerian formulation (ALE) of the Navier-Stokes equations is discretized by finite volumes (FVM), the solid mechanical finite elasticity equations are discretized by a finite element (FEM) approach. Taking advantage of specialized numerical solution strategies for non-matching fluid and solid domain meshes, an iterative data-exchange guarantees the interface equilibrium of the underlying governing equations. The focus of this work is on left-ventricular fluid-structure interaction based on patient-specific magnetic resonance imaging datasets. Multi-physical phenomena are described by temporal visualization and characteristic FSI numbers. The results gained show flow patterns that are in good agreement with previous observations. A deeper understanding of cavity deformation, blood flow, and their vital interaction can help to improve surgical treatment and clinical therapy planning.

Three-dimensional granular model of semi-solid metallic alloys undergoing solidification: Fluid flow and localization of feeding

International Nuclear Information System (INIS)

Sistaninia, M.; Phillion, A.B.; Drezet, J.-M.; Rappaz, M.

2012-01-01

A three-dimensional (3-D) granular model which simulates fluid flow within solidifying alloys with a globular microstructure, such as that found in grain refined Al alloys, is presented. The model geometry within a representative volume element (RVE) consists of a set of prismatic triangular elements representing the intergranular liquid channels. The pressure field within the liquid channels is calculated using a finite elements (FEs) method assuming a Poiseuille flow within each channel and flow conservation at triple lines. The fluid flow is induced by solidification shrinkage and openings at grain boundaries due to deformation of the coherent solid. The granular model predictions are validated against bulk data calculated with averaging techniques. The results show that a fluid flow simulation of globular semi-solid materials is able to reproduce both a map of the 3-D intergranular pressure and the localization of feeding within the mushy zone. A new hot cracking sensitivity coefficient is then proposed. Based on a mass balance performed over a solidifying isothermal volume element, this coefficient accounts for tensile deformation of the semi-solid domain and for the induced intergranular liquid feeding. The fluid flow model is then used to calculate the pressure drop in the mushy zone during the direct chill casting of aluminum alloy billets. The predicted pressure demonstrates that deep in the mushy zone where the permeability is low the local pressure can be significantly lower than the pressure predicted by averaging techniques.

Nonlinear modeling and testing of magneto-rheological fluids in low shear rate squeezing flows

International Nuclear Information System (INIS)

Farjoud, Alireza; Ahmadian, Mehdi; Craft, Michael; Mahmoodi, Nima; Zhang, Xinjie

2011-01-01

A novel analytical investigation of magneto-rheological (MR) fluids in squeezing flows is performed and the results are validated with experimental test data. The squeeze flow of MR fluids has recently been of great interest to researchers. This is due to the large force capacity of MR fluids in squeeze mode compared to other modes (valve and shear modes), which makes the squeeze mode appropriate for a wide variety of applications such as impact dampers and engine mounts. Tested MR fluids were capable of providing a large range of controllable force along a short stroke in squeeze mode. A mathematical model was developed using perturbation techniques to predict closed-form solutions for velocity field, shear rate distribution, pressure distribution and squeeze force. Therefore, the obtained solutions greatly help with the design process of intelligent devices that use MR fluids in squeeze mode. The mathematical model also reduces the need for complicated and computationally expensive numerical simulations. The analytical results are validated by performing experimental tests on a novel MR device called an 'MR pouch' in an MR squeeze mode rheometer, both designed and built at CVeSS

Simulation of uncompressible fluid flow through a porous media

Energy Technology Data Exchange (ETDEWEB)

Ramirez, A. [Instituto Politecnico Nacional (SEPI-ESIQIE-IPN), Unidad Profesional Zacatenco, Laboratorio de Analisis Met. (Edif. ' Z' y Edif. ' 6' P.B.), Mexico City (Mexico)], E-mail: adaramil@yahoo.com.mx; Gonzalez, J.L. [Instituto Politecnico Nacional (SEPI-ESIQIE-IPN), Unidad Profesional Zacatenco, Laboratorio de Analisis Met. (Edif. ' Z' y Edif. ' 6' P.B.), Mexico City (Mexico); Carrillo, F. [Instituto Politecnico Nacional (SEPI-CICATA-IPN), Unidad Altamira Tamaulipas, Mexico (Mexico); Lopez, S. [Instituto Mexicano del Petroleo (I.M.P.-D.F.), Mexico (Mexico)

2009-02-28

Recently, a great interest has been focused for investigations about transport phenomena in disordered systems. One of the most treated topics is fluid flow through anisotropic materials due to the importance in many industrial processes like fluid flow in filters, membranes, walls, oil reservoirs, etc. In this work is described the formulation of a 2D mathematical model to simulate the fluid flow behavior through a porous media (PM) based on the solution of the continuity equation as a function of the Darcy's law for a percolation system; which was reproduced using computational techniques reproduced using a random distribution of the porous media properties (porosity, permeability and saturation). The model displays the filling of a partially saturated porous media with a new injected fluid showing the non-defined advance front and dispersion of fluids phenomena.

Mathematical Modeling of Fluid Flow in a Water Physical Model of an Aluminum Degassing Ladle Equipped with an Impeller-Injector

Science.gov (United States)

Gómez, Eudoxio Ramos; Zenit, Roberto; Rivera, Carlos González; Trápaga, Gerardo; Ramírez-Argáez, Marco A.

2013-04-01

In this work, a 3D numerical simulation using a Euler-Euler-based model implemented into a commercial CFD code was used to simulate fluid flow and turbulence structure in a water physical model of an aluminum ladle equipped with an impeller for degassing treatment. The effect of critical process parameters such as rotor speed, gas flow rate, and the point of gas injection (conventional injection through the shaft vs a novel injection through the bottom of the ladle) on the fluid flow and vortex formation was analyzed with this model. The commercial CFD code PHOENICS 3.4 was used to solve all conservation equations governing the process for this two-phase fluid flow system. The mathematical model was reasonably well validated against experimentally measured liquid velocity and vortex sizes in a water physical model built specifically for this investigation. From the results, it was concluded that the angular speed of the impeller is the most important parameter in promoting better stirred baths and creating smaller and better distributed bubbles in the liquid. The pumping effect of the impeller is increased as the impeller rotation speed increases. Gas flow rate is detrimental to bath stirring and diminishes the pumping effect of the impeller. Finally, although the injection point was the least significant variable, it was found that the "novel" injection improves stirring in the ladle.

A coupled model on fluid flow, heat transfer and solidification in continuous casting mold

Directory of Open Access Journals (Sweden)

Xu-bin Zhang

2017-11-01

Full Text Available Fluid flow, heat transfer and solidification of steel in the mold are so complex but crucial, determining the surface quality of the continuous casting slab. In the current study, a 2D numerical model was established by Fluent software to simulate the fluid flow, heat transfer and solidification of the steel in the mold. The VOF model and k-ε model were applied to simulate the flow field of the three phases (steel, slag and air, and solidification model was used to simulate the solidification process. The phenomena at the meniscus were also explored through interfacial tension between the liquid steel and slag as well as the mold oscillation. The model included a 20 mm thick mold to clarify the heat transfer and the temperature distribution of the mold. The simulation results show that the liquid steel flows as upper backflow and lower backflow in the mold, and that a small circulation forms at the meniscus. The liquid slag flows away from the corner at the meniscus or infiltrates into the gap between the mold and the shell with the mold oscillating at the negative strip stage or at the positive strip stage. The simulated pitch and the depth of oscillation marks approximate to the theoretical pitch and measured depth on the slab.

Computational fluid dynamics modeling of mixed convection flows in buildings enclosures

Energy Technology Data Exchange (ETDEWEB)

Kayne, Alexander; Agarwal, Ramesh K. [Department of Mechanical Engineering and Materials Science, Washington University, St. Louis, MO 63130 (United States)

2013-07-01

In recent years Computational Fluid Dynamics (CFD) simulations are increasingly used to model the air circulation and temperature environment inside the rooms of residential and office buildings to gain insight into the relative energy consumptions of various HVAC systems for cooling/heating for climate control and thermal comfort. This requires accurate simulation of turbulent flow and heat transfer for various types of ventilation systems using the Reynolds-Averaged Navier-Stokes (RANS) equations of fluid dynamics. Large Eddy Simulation (LES) or Direct Numerical Simulation (DNS) of Navier-Stokes equations is computationally intensive and expensive for simulations of this kind. As a result, vast majority of CFD simulations employ RANS equations in conjunction with a turbulence model. In order to assess the modeling requirements (mesh, numerical algorithm, turbulence model etc.) for accurate simulations, it is critical to validate the calculations against the experimental data. For this purpose, we use three well known benchmark validation cases, one for natural convection in 2D closed vertical cavity, second for forced convection in a 2D rectangular cavity and the third for mixed convection in a 2D square cavity. The simulations are performed on a number of meshes of different density using a number of turbulence models. It is found that k-epsilon two-equation turbulence model with a second-order algorithm on a reasonable mesh gives the best results. This information is then used to determine the modeling requirements (mesh, numerical algorithm, turbulence model etc.) for flows in 3D enclosures with different ventilation systems. In particular two cases are considered for which the experimental data is available. These cases are (1) air flow and heat transfer in a naturally ventilated room and (2) airflow and temperature distribution in an atrium. Good agreement with the experimental data and computations of other investigators is obtained.

Modeling studies for multiphase fluid and heat flow processes in nuclear waste isolation

International Nuclear Information System (INIS)

Pruess, K.

1988-07-01

Multiphase fluid and heat flow plays an important role in many problems relating to the disposal of nuclear wastes in geologic media. Examples include boiling and condensation processes near heat-generating wastes, flow of water and formation gas in partially saturated formations, evolution of a free gas phase from waste package corrosion in initially water-saturated environments, and redistribution (dissolution, transport, and precipitation) of rock minerals in non-isothermal flow fields. Such processes may strongly impact upon waste package and repository design considerations and performance. This paper summarizes important physical phenomena occurring in multiphase and nonisothermal flows, as well as techniques for their mathematical modeling and numerical simulation. Illustrative applications are given for a number of specific fluid and heat flow problems, including: thermohydrologic conditions near heat-generating waste packages in the unsaturated zone; repository-wide convection effects in the unsaturated zone; effects of quartz dissolution and precipitation for disposal in the saturated zone; and gas pressurization and flow corrosion of low-level waste packages. 34 refs; 7 figs; 2 tabs

Lattice Boltzmann Simulations of Fluid Flow in Continental Carbonate Reservoir Rocks and in Upscaled Rock Models Generated with Multiple-Point Geostatistics

Directory of Open Access Journals (Sweden)

J. Soete

2017-01-01

Full Text Available Microcomputed tomography (μCT and Lattice Boltzmann Method (LBM simulations were applied to continental carbonates to quantify fluid flow. Fluid flow characteristics in these complex carbonates with multiscale pore networks are unique and the applied method allows studying their heterogeneity and anisotropy. 3D pore network models were introduced to single-phase flow simulations in Palabos, a software tool for particle-based modelling of classic computational fluid dynamics. In addition, permeability simulations were also performed on rock models generated with multiple-point geostatistics (MPS. This allowed assessing the applicability of MPS in upscaling high-resolution porosity patterns into large rock models that exceed the volume limitations of the μCT. Porosity and tortuosity control fluid flow in these porous media. Micro- and mesopores influence flow properties at larger scales in continental carbonates. Upscaling with MPS is therefore necessary to overcome volume-resolution problems of CT scanning equipment. The presented LBM-MPS workflow is applicable to other lithologies, comprising different pore types, shapes, and pore networks altogether. The lack of straightforward porosity-permeability relationships in complex carbonates highlights the necessity for a 3D approach. 3D fluid flow studies provide the best understanding of flow through porous media, which is of crucial importance in reservoir modelling.

Numerical modelling of thermal and fluid flow phenomena in the mould channel

Directory of Open Access Journals (Sweden)

L. Sowa

2007-12-01

Full Text Available In the paper, a mathematical and a numerical model of the solidification of a cylindrical slender shaped casting, which take into account the process of filling the mould cavity with molten metal, has been proposed. Pressure and velocity fields were obtained by solving the momentum equations and the continuity equation, while the thermal fields were obtained by solving the heat conduction equation containing the convection term. Next, the numerical analysis of the solidification process of metals alloy in a cylindrical mould channel has been made. In the model one takes into account interdependence the heat transfer and fluid flow phenomena. Coupling of the thermal and fluid flow phenomena has been taken into consideration by the changes of the fluidity function and thermophysical parameters of alloy with respect to the temperature. The influence of the pressure and the temperature of metal pouring on the solid phase growth kinetics were estimated. The problem has been solved by the finite element method.

Experimental and CFD flow studies in an intracranial aneurysm model with Newtonian and non-Newtonian fluids.

Science.gov (United States)

Frolov, S V; Sindeev, S V; Liepsch, D; Balasso, A

2016-05-18

According to the clinical data, flow conditions play a major role in the genesis of intracranial aneurysms. The disorder of the flow structure is the cause of damage of the inner layer of the vessel wall, which leads to the development of cerebral aneurysms. Knowledge of the alteration of the flow field in the aneurysm region is important for treatment. The aim is to study quantitatively the flow structure in an patient-specific aneurysm model of the internal carotid artery using both experimental and computational fluid dynamics (CFD) methods with Newtonian and non-Newtonian fluids. A patient-specific geometry of aneurysm of the internal carotid artery was used. Patient data was segmented and smoothed to obtain geometrical model. An elastic true-to-scale silicone model was created with stereolithography. For initial investigation of the blood flow, the flow was visualized by adding particles into the silicone model. The precise flow velocity measurements were done using 1D Laser Doppler Anemometer with a spatial resolution of 50 μ m and a temporal resolution of 1 ms. The local velocity measurements were done at a distance of 4 mm to each other. A fluid with non-Newtonian properties was used in the experiment. The CFD simulations for unsteady-state problem were done using constructed hexahedral mesh for Newtonian and non-Newtonian fluids. Using 1D laser Doppler Anemometer the minimum velocity magnitude at the end of systole -0.01 m/s was obtained in the aneurysm dome while the maximum velocity 1 m/s was at the center of the outlet segment. On central cross section of the aneurysm the maximum velocity value is only 20% of the average inlet velocity. The average velocity on the cross-section is only 11% of the inlet axial velocity. Using the CFD simulation the wall shear stresses for Newtonian and non-Newtonian fluid at the end of systolic phase (t= 0.25 s) were computed. The wall shear stress varies from 3.52 mPa (minimum value) to 10.21 Pa (maximum value) for the

A Well-Posed Two Phase Flow Model and its Numerical Solutions for Reactor Thermal-Fluids Analysis

Energy Technology Data Exchange (ETDEWEB)

Kadioglu, Samet Y. [Idaho National Lab. (INL), Idaho Falls, ID (United States); Berry, Ray [Idaho National Lab. (INL), Idaho Falls, ID (United States); Martineau, Richard [Idaho National Lab. (INL), Idaho Falls, ID (United States)

2016-08-01

A 7-equation two-phase flow model and its numerical implementation is presented for reactor thermal-fluids applications. The equation system is well-posed and treats both phases as compressible flows. The numerical discretization of the equation system is based on the finite element formalism. The numerical algorithm is implemented in the next generation RELAP-7 code (Idaho National Laboratory (INL)’s thermal-fluids code) built on top of an other INL’s product, the massively parallel multi-implicit multi-physics object oriented code environment (MOOSE). Some preliminary thermal-fluids computations are presented.

A Well-Posed Two Phase Flow Model and its Numerical Solutions for Reactor Thermal-Fluids Analysis

International Nuclear Information System (INIS)

Kadioglu, Samet Y.; Berry, Ray; Martineau, Richard

2016-01-01

A 7-equation two-phase flow model and its numerical implementation is presented for reactor thermal-fluids applications. The equation system is well-posed and treats both phases as compressible flows. The numerical discretization of the equation system is based on the finite element formalism. The numerical algorithm is implemented in the next generation RELAP-7 code (Idaho National Laboratory (INL)'s thermal-fluids code) built on top of an other INL's product, the massively parallel multi-implicit multi-physics object oriented code environment (MOOSE). Some preliminary thermal-fluids computations are presented.

A discontinuous finite element approach to cracking in coupled poro-elastic fluid flow models

Science.gov (United States)

Wilson, C. R.; Spiegelman, M. W.; Evans, O.; Ulven, O. I.; Sun, W.

2016-12-01

Reaction-driven cracking is a coupled process whereby fluid-induced reactions drive large volume changes in the host rock which produce stresses leading to crack propagation and failure. This in turn generates new surface area and fluid-flow pathways for subsequent reaction in a potentially self-sustaining system. This mechanism has has been proposed for the pervasive serpentinization and carbonation of peridotite, as well as applications to mineral carbon sequestration and hydrocarbon extraction. The key computational issue in this problem is implementing algorithms that adequately model the formation of discrete fractures. Here we present models using a discontinuous finite element method for modeling fracture formation (Radovitsky et al., 2011). Cracks are introduced along facets of the mesh by the relaxation of penalty parameters once a failure criterion is met. It is fully described in the weak form of the equations, requiring no modification of the underlying mesh structure and allowing fluid properties to be easily adjusted along cracked facets. To develop and test the method, we start by implementing the algorithm for the simplified Biot equations for poro-elasticity using the finite element model assembler TerraFERMA. We consider hydro-fracking around a borehole (Grassl et al., 2015), where elevated fluid pressure in the poro-elastic solid causes it to fail radially in tension. We investigate the effects of varying the Biot coefficient and adjusting the fluid transport properties in the vicinity of the crack and compare our results to related dual-graph models (Ulven & Sun, submitted). We discuss issues arising from this method, including the formation of null spaces and appropriate preconditioning and solution strategies. Initial results suggest that this method provides a promising way to incorporate cracking into our reactive fluid flow models and future work aims to integrate the mechanical and chemical aspects of this process.

A mathematical model of fluid and gas flow in nanoporous media.

Science.gov (United States)

Monteiro, Paulo J M; Rycroft, Chris H; Barenblatt, Grigory Isaakovich

2012-12-11

The mathematical modeling of the flow in nanoporous rocks (e.g., shales) becomes an important new branch of subterranean fluid mechanics. The classic approach that was successfully used in the construction of the technology to develop oil and gas deposits in the United States, Canada, and the Union of Soviet Socialist Republics becomes insufficient for deposits in shales. In the present article a mathematical model of the flow in nanoporous rocks is proposed. The model assumes the rock consists of two components: (i) a matrix, which is more or less an ordinary porous or fissurized-porous medium, and (ii) specific organic inclusions composed of kerogen. These inclusions may have substantial porosity but, due to the nanoscale of pores, tubes, and channels, have extremely low permeability on the order of a nanodarcy (~109-²¹ m² ) or less. These inclusions contain the majority of fluid: oil and gas. Our model is based on the hypothesis that the permeability of the inclusions substantially depends on the pressure gradient. At the beginning of the development of the deposit, boundary layers are formed at the boundaries of the low-permeable inclusions, where the permeability is strongly increased and intensive flow from inclusions to the matrix occurs. The resulting formulae for the production rate of the deposit are presented in explicit form. The formulae demonstrate that the production rate of deposits decays with time following a power law whose exponent lies between -1/2 and -1. Processing of experimental data obtained from various oil and gas deposits in shales demonstrated an instructive agreement with the prediction of the model.

Hydraulic Fracturing and Production Optimization in Eagle Ford Shale Using Coupled Geomechanics and Fluid Flow Model

Science.gov (United States)

Suppachoknirun, Theerapat; Tutuncu, Azra N.

2017-12-01

With increasing production from shale gas and tight oil reservoirs, horizontal drilling and multistage hydraulic fracturing processes have become a routine procedure in unconventional field development efforts. Natural fractures play a critical role in hydraulic fracture growth, subsequently affecting stimulated reservoir volume and the production efficiency. Moreover, the existing fractures can also contribute to the pressure-dependent fluid leak-off during the operations. Hence, a reliable identification of the discrete fracture network covering the zone of interest prior to the hydraulic fracturing design needs to be incorporated into the hydraulic fracturing and reservoir simulations for realistic representation of the in situ reservoir conditions. In this research study, an integrated 3-D fracture and fluid flow model have been developed using a new approach to simulate the fluid flow and deliver reliable production forecasting in naturally fractured and hydraulically stimulated tight reservoirs. The model was created with three key modules. A complex 3-D discrete fracture network model introduces realistic natural fracture geometry with the associated fractured reservoir characteristics. A hydraulic fracturing model is created utilizing the discrete fracture network for simulation of the hydraulic fracture and flow in the complex discrete fracture network. Finally, a reservoir model with the production grid system is used allowing the user to efficiently perform the fluid flow simulation in tight formations with complex fracture networks. The complex discrete natural fracture model, the integrated discrete fracture model for the hydraulic fracturing, the fluid flow model, and the input dataset have been validated against microseismic fracture mapping and commingled production data obtained from a well pad with three horizontal production wells located in the Eagle Ford oil window in south Texas. Two other fracturing geometries were also evaluated to optimize

A two-dimensional continuum model of biofilm growth incorporating fluid flow and shear stress based detachment

KAUST Repository

Duddu, Ravindra; Chopp, David L.; Moran, Brian

2009-01-01

of the biofilm. The model considers fluid flow around the biofilm surface, the advection-diffusion and reaction of substrate, variable biomass volume fraction and erosion due to the interfacial shear stress at the biofilm-fluid interface. The key assumptions

Standardization of Thermo-Fluid Modeling in Modelica.Fluid

Energy Technology Data Exchange (ETDEWEB)

Franke, Rudiger; Casella, Francesco; Sielemann, Michael; Proelss, Katrin; Otter, Martin; Wetter, Michael

2009-09-01

This article discusses the Modelica.Fluid library that has been included in the Modelica Standard Library 3.1. Modelica.Fluid provides interfaces and basic components for the device-oriented modeling of onedimensional thermo-fluid flow in networks containing vessels, pipes, fluid machines, valves and fittings. A unique feature of Modelica.Fluid is that the component equations and the media models as well as pressure loss and heat transfer correlations are decoupled from each other. All components are implemented such that they can be used for media from the Modelica.Media library. This means that an incompressible or compressible medium, a single or a multiple substance medium with one or more phases might be used with one and the same model as long as the modeling assumptions made hold. Furthermore, trace substances are supported. Modeling assumptions can be configured globally in an outer System object. This covers in particular the initialization, uni- or bi-directional flow, and dynamic or steady-state formulation of mass, energy, and momentum balance. All assumptions can be locally refined for every component. While Modelica.Fluid contains a reasonable set of component models, the goal of the library is not to provide a comprehensive set of models, but rather to provide interfaces and best practices for the treatment of issues such as connector design and implementation of energy, mass and momentum balances. Applications from various domains are presented.

Using FlowLab, an educational computational fluid dynamics tool, to perform a comparative study of turbulence models

International Nuclear Information System (INIS)

Parihar, A.; Kulkarni, A.; Stern, F.; Xing, T.; Moeykens, S.

2005-01-01

Flow over an Ahmed body is a key benchmark case for validating the complex turbulent flow field around vehicles. In spite of the simple geometry, the flow field around an Ahmed body retains critical features of real, external vehicular flow. The present study is an attempt to implement such a real life example into the course curriculum for undergraduate engineers. FlowLab, which is a Computational Fluid Dynamics (CFD) tool developed by Fluent Inc. for use in engineering education, allows students to conduct interactive application studies. This paper presents a synopsis of FlowLab, a description of one FlowLab exercise, and an overview of the educational experience gained by students through using FlowLab, which is understood through student surveys and examinations. FlowLab-based CFD exercises were implemented into 57:020 Mechanics of Fluids and Transport Processes and 58:160 Intermediate Mechanics of Fluids courses at the University of Iowa in the fall of 2004, although this report focuses only on experiences with the Ahmed body exercise, which was used only in the intermediate-level fluids class, 58:160. This exercise was developed under National Science Foundation funding by the authors of this paper. The focus of this study does not include validating the various turbulence models used for the Ahmed body simulation, because a two-dimensional simplification was applied. With the two-dimensional simplification, students may setup, run, and post process this model in a 50 minute class period using a single-CPU PC, as required for the 58:160 class at the University of Iowa. It is educational for students to understand the implication of a two- dimensional approximation for essentially a three-dimensional flow field, along with the consequent variation in both qualitative and quantitative results. Additionally, through this exercise, students may realize that the choice of the respective turbulence model will affect simulation prediction. (author)

Fluid flow control system

International Nuclear Information System (INIS)

Rion, Jacky.

1982-01-01

Fluid flow control system featuring a series of grids placed perpendicular to the fluid flow direction, characterized by the fact that it is formed of a stack of identical and continuous grids, each of which consists of identical meshes forming a flat lattice. The said meshes are offset from one grid to the next. This system applies in particular to flow control of the coolant flowing at the foot of an assembly of a liquid metal cooled nuclear reactor [fr

Resolution of through tubing fluid flow and behind casing fluid flow in multiple completion wells

International Nuclear Information System (INIS)

Arnold, D.M.

1977-01-01

A method is provided for resolving undesired fluid flow in cement channels behind casing in one producing zone of a multi zone completion well operating on gas lift from the fluid flow from lower producing zones in the same well which is contained in production tubing passing through the producing zone being investigated. Gamma rays which are characteristic of the decay of the unstable isotope nitrogen 16 produced by activation of elemental oxygen nuclei comprising the molecular structure of both the tubing fluid flow and the undesired fluid flow are detected in at least two energy bonds at two longitudinally spaced detectors in a well borehole. By appropriately combining the four count rate signals so producing according to predetermined relationships the two fluid flow components in the same direction may be uniquely distinguished on the basis of their differing distances from the gamma ray detectors. 9 claims, 17 figures

Simulation of the two-fluid model on incompressible flow with Fractional Step method for both resolved and unresolved scale interfaces

International Nuclear Information System (INIS)

Hou, Xiaofei; Rigola, Joaquim; Lehmkuhl, Oriol; Oliva, Assensi

2015-01-01

Highlights: • Two phase flow with free surface is solved by means of two-fluid model (TFM). • Fractional Step method and finite volume technique is used to solve TFM. • Conservative Level Set method reduces interface sharpening diffusion problem. • Cases including high density ratios and high viscosities validate the models. - Abstract: In the present paper, the Fractional Step method usually used in single fluid flow is here extended and applied for the two-fluid model resolution using the finite volume discretization. The use of a projection method resolution instead of the usual pressure-correction method for multi-fluid flow, successfully avoids iteration processes. On the other hand, the main weakness of the two fluid model used for simulations of free surface flows, which is the numerical diffusion of the interface, is also solved by means of the conservative Level Set method (interface sharpening) (Strubelj et al., 2009). Moreover, the use of the algorithm proposed has allowed presenting different free-surface cases with or without Level Set implementation even under coarse meshes under a wide range of density ratios. Thus, the numerical results presented, numerically verified, experimentally validated and converged under high density ratios, shows the capability and reliability of this resolution method for both mixed and unmixed flows

Local invariants in non-ideal flows of neutral fluids and two-fluid plasmas

Science.gov (United States)

Zhu, Jian-Zhou

2018-03-01

The main objective is the locally invariant geometric object of any (magneto-)fluid dynamics with forcing and damping (nonideal), while more attention is paid to the untouched dynamical properties of two-fluid fashion. Specifically, local structures, beyond the well-known "frozen-in" to the barotropic flows of the generalized vorticities, of the two-fluid model of plasma flows are presented. More general non-barotropic situations are also considered. A modified Euler equation [T. Tao, "Finite time blowup for Lagrangian modifications of the three-dimensional Euler equation," Ann. PDE 2, 9 (2016)] is also accordingly analyzed and remarked from the angle of view of the two-fluid model, with emphasis on the local structures. The local constraints of high-order differential forms such as helicity, among others, find simple formulation for possible practices in modeling the dynamics. Thus, the Cauchy invariants equation [N. Besse and U. Frisch, "Geometric formulation of the Cauchy invariants for incompressible Euler flow in flat and curved spaces," J. Fluid Mech. 825, 412 (2017)] may be enabled to find applications in non-ideal flows. Some formal examples are offered to demonstrate the calculations, and particularly interestingly the two-dimensional-three-component (2D3C) or the 2D passive scalar problem presents that a locally invariant Θ = 2θζ, with θ and ζ being, respectively, the scalar value of the "vertical velocity" (or the passive scalar) and the "vertical vorticity," may be used as if it were the spatial density of the globally invariant helicity, providing a Lagrangian prescription to control the latter in some situations of studying its physical effects in rapidly rotating flows (ubiquitous in atmosphere of astrophysical objects) with marked 2D3C vortical modes or in purely 2D passive scalars.

Simulation based engineering in fluid flow design

CERN Document Server

Rao, J S

2017-01-01

This volume offers a tool for High Performance Computing (HPC). A brief historical background on the subject is first given. Fluid Statics dealing with Pressure in fluids at rest, Buoyancy and Basics of Thermodynamics are next presented. The Finite Volume Method, the most convenient process for HPC, is explained in one-dimensional approach to diffusion with convection and pressure velocity coupling. Adiabatic, isentropic and supersonic flows in quasi-one dimensional flows in axisymmetric nozzles is considered before applying CFD solutions. Though the theory is restricted to one-dimensional cases, three-dimensional CFD examples are also given. Lastly, nozzle flows with normal shocks are presented using turbulence models. Worked examples and exercises are given in each chapter. Fluids transport thermal energy for its conversion to kinetic energy, thus playing a major role that is central to all heat engines. With the advent of rotating machinery in the 20th century, Fluid Engineering was developed in the form o...

Modeling fluid forces and response of a tube bundle in cross-flow induced vibrations

International Nuclear Information System (INIS)

Khushnood, Shahab; Khan, Zaffar M.; Malik, M. Afzaal; Koreshi, Zafarullah; Khan, Mahmood Anwar

2003-01-01

Flow induced vibrations occur in process heat exchangers, condensers, boilers and nuclear steam generators. Under certain flow conditions and fluid velocities, the fluid forces result in tube vibrations and possible damage of tube, tube sheet or baffle due to fretting and fatigue. Prediction of these forces is an important consideration. The characteristics of vibration depend greatly on the fluid dynamic forces and structure of the tube bundle. It is undesirable for the tube bundles to vibrate excessively under normal operating conditions because tubes wear and eventual leakage can occur leading to costly shutdowns. In this paper modeling of fluid forces and vibration response of a tube in a heat exchanger bundle has been carried out. Experimental validation has been performed on an existing refinery heat exchanger tube bundle. The target tube has been instrumented with an accelerometer and strain gages. The bundle has been studied for pulse, sinusoidal and random excitations. Natural frequencies and damping of the tubes have also been computed. Experimental fluid forces and response shows a reasonable agreement with the predictions. (author)

Implicit approximate Riemann solver for two fluid two phase flow models

International Nuclear Information System (INIS)

Raymond, P.; Toumi, I.; Kumbaro, A.

1993-01-01

This paper is devoted to the description of new numerical methods developed for the numerical treatment of two phase flow models with two velocity fields which are now widely used in nuclear engineering for design or safety calculations. These methods are finite volumes numerical methods and are based on the use of Approximate Riemann Solver's concepts in order to define convective flux versus mean cell quantities. The first part of the communication will describe the numerical method for a three dimensional drift flux model and the extensions which were performed to make the numerical scheme implicit and to have fast running calculations of steady states. Such a scheme is now implemented in the FLICA-4 computer code devoted to 3-D steady state and transient core computations. We will present results obtained for a steady state flow with rod bow effect evaluation and for a Steam Line Break calculation were the 3-D core thermal computation was coupled with a 3-D kinetic calculation and a thermal-hydraulic transient calculation for the four loops of a Pressurized Water Reactor. The second part of the paper will detail the development of an equivalent numerical method based on an approximate Riemann Solver for a two fluid model with two momentum balance equations for the liquid and the gas phases. The main difficulty for these models is due to the existence of differential modelling terms such as added mass effects or interfacial pressure terms which make hyperbolic the model. These terms does not permit to write the balance equations system in a conservative form, and the classical theory for discontinuity propagation for non-linear systems cannot be applied. Meanwhile, the use of non-conservative products theory allows the study of discontinuity propagation for a non conservative model and this will permit the construction of a numerical scheme for two fluid two phase flow model. These different points will be detailed in that section which will be illustrated by

A two-fluid model for vertical flow applied to CO₂ injection wells

DEFF Research Database (Denmark)

Linga, Gaute; Lund, Halvor

2016-01-01

Flow of CO2 in wells is associated with substantial variations in thermophysical properties downhole, due to the coupled transient processes involved: complex flow patterns, density changes, phase transitions, and heat transfer to and from surroundings. Large temperature variations can lead...... the well, including tubing, packer fluid, casing, cement or drilling mud, and rock formation. This enables prediction of the temperature in the well fluid and in each layer of the well. The model is applied to sudden shut-in and blowout cases of a CO2 injection well, where we employ the highly accurate...

Modelling of hydrothermal fluid flow and structural architecture in an extensional basin, Ngakuru Graben, Taupo Rift, New Zealand

Science.gov (United States)

Kissling, W. M.; Villamor, P.; Ellis, S. M.; Rae, A.

2018-05-01

Present-day geothermal activity on the margins of the Ngakuru graben and evidence of fossil hydrothermal activity in the central graben suggest that a graben-wide system of permeable intersecting faults acts as the principal conduit for fluid flow to the surface. We have developed numerical models of fluid and heat flow in a regional-scale 2-D cross-section of the Ngakuru Graben. The models incorporate simplified representations of two 'end-member' fault architectures (one symmetric at depth, the other highly asymmetric) which are consistent with the surface locations and dips of the Ngakuru graben faults. The models are used to explore controls on buoyancy-driven convective fluid flow which could explain the differences between the past and present hydrothermal systems associated with these faults. The models show that the surface flows from the faults are strongly controlled by the fault permeability, the fault system architecture and the location of the heat source with respect to the faults in the graben. In particular, fault intersections at depth allow exchange of fluid between faults, and the location of the heat source on the footwall of normal faults can facilitate upflow along those faults. These controls give rise to two distinct fluid flow regimes in the fault network. The first, a regular flow regime, is characterised by a nearly unchanging pattern of fluid flow vectors within the fault network as the fault permeability evolves. In the second, complex flow regime, the surface flows depend strongly on fault permeability, and can fluctuate in an erratic manner. The direction of flow within faults can reverse in both regimes as fault permeability changes. Both flow regimes provide insights into the differences between the present-day and fossil geothermal systems in the Ngakuru graben. Hydrothermal upflow along the Paeroa fault seems to have occurred, possibly continuously, for tens of thousands of years, while upflow in other faults in the graben has

Prediction of subcooled flow boiling characteristics using two-fluid Eulerian CFD model

Energy Technology Data Exchange (ETDEWEB)

Braz Filho, Francisco A.; Ribeiro, Guilherme B., E-mail: gbribeiro@ieav.cta.br; Caldeira, Alexandre D.

2016-11-15

Highlights: • CFD multiphase model is used to predict subcooled flow boiling characteristics. • Better agreement is achieved for higher saturation pressures. • Onset of nucleate boiling and saturated boiling are well predicted. • CFD multiphase model tends to underestimate the void fraction. • Factors were adjusted in order to improve the void fraction results. - Abstract: The present study concerns a detailed analysis of flow boiling phenomena under high pressure systems using a two-fluid Eulerian approach provided by a Computational Fluid Dynamics (CFD) solver. For this purpose, a vertical heated pipe made of stainless steel with an internal diameter of 15.4 mm was considered as the modeled domain. Two different uniform heat fluxes and three saturation pressures were applied to the channel wall, whereas water mass flux of 900 kg/m{sup 2} s was considered for all simulation cases. The model was validated against a set of experimental data and results have indicated a promising use of the CFD technique for estimation of the wall temperature, the liquid bulk temperature and the location of the departure of nucleate boiling. Changes in factors applied in the modeling of the interfacial heat transfer coefficient and bubble departure frequency were suggested, allowing a better prediction of the void fraction along the heated channel. The commercial CFD solver FLUENT 14.5 was used for the model implementation.

Prediction of subcooled flow boiling characteristics using two-fluid Eulerian CFD model

International Nuclear Information System (INIS)

Braz Filho, Francisco A.; Ribeiro, Guilherme B.; Caldeira, Alexandre D.

2016-01-01

Highlights: • CFD multiphase model is used to predict subcooled flow boiling characteristics. • Better agreement is achieved for higher saturation pressures. • Onset of nucleate boiling and saturated boiling are well predicted. • CFD multiphase model tends to underestimate the void fraction. • Factors were adjusted in order to improve the void fraction results. - Abstract: The present study concerns a detailed analysis of flow boiling phenomena under high pressure systems using a two-fluid Eulerian approach provided by a Computational Fluid Dynamics (CFD) solver. For this purpose, a vertical heated pipe made of stainless steel with an internal diameter of 15.4 mm was considered as the modeled domain. Two different uniform heat fluxes and three saturation pressures were applied to the channel wall, whereas water mass flux of 900 kg/m"2 s was considered for all simulation cases. The model was validated against a set of experimental data and results have indicated a promising use of the CFD technique for estimation of the wall temperature, the liquid bulk temperature and the location of the departure of nucleate boiling. Changes in factors applied in the modeling of the interfacial heat transfer coefficient and bubble departure frequency were suggested, allowing a better prediction of the void fraction along the heated channel. The commercial CFD solver FLUENT 14.5 was used for the model implementation.

Simultaneous effects of slip and MHD on peristaltic blood flow of Jeffrey fluid model through a porous medium

Directory of Open Access Journals (Sweden)

M.M. Bhatti

2016-06-01

Full Text Available In this article, the simultaneous effects of slip and Magnetohydrodynamics (MHD on peristaltic blood flow of Jeffrey fluid model have been investigated in a non-uniform porous channel. The governing equation of blood flow for Jeffrey fluid model is solved with the help of long wavelength and creeping flow regime. The solution of the resulting differential equation is solved analytically and a closed form solution is presented. The impact of all the physical parameters is plotted for velocity profile and pressure rise. Nowadays, Magnetohydrodynamics is applicable in various magnetic drug targeting for cancer diseases and also very helpful to control the flow. The present analysis is also described for Newtonian fluid (λ1→0 as a special case of our study. It is observed that magnitude of the velocity is opposite near the walls due to slip effects whereas similar behavior has been observed for magnetic field.

A PISO-like algorithm to simulate superfluid helium flow with the two-fluid model

CERN Document Server

Soulaine, Cyprien; Allain, Hervé; Baudouy, Bertrand; Van Weelderen, Rob

2015-01-01

This paper presents a segregated algorithm to solve numerically the superfluid helium (He II) equations using the two-fluid model. In order to validate the resulting code and illustrate its potential, different simulations have been performed. First, the flow through a capillary filled with He II with a heated area on one side is simulated and results are compared to analytical solutions in both Landau and Gorter–Mellink flow regimes. Then, transient heat transfer of a forced flow of He II is investigated. Finally, some two-dimensional simulations in a porous medium model are carried out.

Unsteady magnetohydrodynamics micropolar fluid in boundary layer flow past a sphere influenced by magnetic fluid

Science.gov (United States)

Pratomo, Rizky Verdyanto; Widodo, Basuki; Adzkiya, Dieky

2017-12-01

Research about fluid flow was very interesting because have a lot of advantages and it can be applied in many aspects of life. The study on fluid flow which is now widely studied is on magnetohydrodynamic (MHD). Magnetohydrodynamic is a conductive and electrical in a magnetic field. This paper considers the effect of unsteady magnetic fields on the flow of magneto-hydrodynamic fluid on the boundary layer that flows past a sphere in micropolar fluid influenced by magnetic field. Our approach is as follows. First, we construct a mathematical model and then the system of equations obtained will be solved numerically using the Keller-Box scheme. Then the system is simulated to assess its effect on the fluid flow velocity profile and the profile of microrotation particles. The result of this research indicates, that when the magnetic parameters increase, then velocity profile increases. If material parameters increase, then velocity profile decreases and magnetic parameters increase for n = 0. For n = 0.5, if magnetic parameters increase, then microrotation profile decreases.

The analysis of two-phase flow and heat transfer using a multidimensional, four field, two-fluid model

International Nuclear Information System (INIS)

Lahey, Richard T.; Drew, Donald A.

2001-01-01

This paper reviews the state-of-the-art in the prediction of multidimensional multiphase flow and heat transfer phenomena using a four field, two-fluid model. It is shown that accurate mechanistic computational fluid dynamic (CFD) predictions are possible for a wide variety of adiabatic and diabatic flows using this computational model. In particular, the model is able to predict the bubbly air/water upflow data of Serizawa (Serizawa, A., 1974. Fluid dynamic characteristics of two-phase flow. Ph.D. thesis, (Nuclear Engineering), Kyoto University, Japan), the downflow data of Wang et al. (Wang, S.K., Lee, S.J., Lahey Jr., R.T., Jones, O.C., 1987. 3-D turbulence structure and phase distribution measurements in bubbly two-phase flows. Int. J. Multiphase Flow 13 (3), 327-343), the isosceles triangle upflow data of Lopez de Bertodano et al. (Lopez de Bertodano, M., Lahey Jr., R.T., Jones, O.C., 1994b. Phase distribution in bubbly two-phase flow in vertical ducts. Int. J. Multiphase Flow 20 (5), 805-818), the heated annular R-113 subcooled boiling data of Velidandala, et al. (Velidandla, V., Pulta, S., Roy, P., Kaira, S.P., 1995. Velocity field in turbulent subcooled boiling flow. ASME Preprint HTD-314, 107-123) and the R-113 CHF data of Hino and Ueda (Hino, R., Ueda, T., 1985. Studies on heat transfer and flow characteristics in subcooled boiling-part 2, flow characteristics. Int. J. Multiphase Flow 11, 283-297). It can also predict external two-phase flows, such as those for spreading two-phase jets (Bonetto, F., Lahey Jr., R.T., 1993. An experimental study on air carryunder due to a plunging liquid jet. Int. J. Multiphase Flow 19 (2), 281-294) and multiphase flows around the hull of naval surface ships (Carrica, P.M., Bonetto, F., Drew, D.A., Lahey, R.T., 1999. A polydispersed model for bubbly two-phase flow around a surface ship. Int. J. Multiphase Flow 25 (2), 257-305)

Flow Diode and Method for Controlling Fluid Flow Origin of the Invention

Science.gov (United States)

Dyson, Rodger W (Inventor)

2015-01-01

A flow diode configured to permit fluid flow in a first direction while preventing fluid flow in a second direction opposite the first direction is disclosed. The flow diode prevents fluid flow without use of mechanical closures or moving parts. The flow diode utilizes a bypass flowline whereby all fluid flow in the second direction moves into the bypass flowline having a plurality of tortuous portions providing high fluidic resistance. The portions decrease in diameter such that debris in the fluid is trapped. As fluid only travels in one direction through the portions, the debris remains trapped in the portions.

Coupling two-phase fluid flow with two-phase darcy flow in anisotropic porous media

KAUST Repository

Chen, J.

2014-06-03

This paper reports a numerical study of coupling two-phase fluid flow in a free fluid region with two-phase Darcy flow in a homogeneous and anisotropic porous medium region. The model consists of coupled Cahn-Hilliard and Navier-Stokes equations in the free fluid region and the two-phase Darcy law in the anisotropic porous medium region. A Robin-Robin domain decomposition method is used for the coupled Navier-Stokes and Darcy system with the generalized Beavers-Joseph-Saffman condition on the interface between the free flow and the porous media regions. Obtained results have shown the anisotropic properties effect on the velocity and pressure of the two-phase flow. 2014 Jie Chen et al.

Coupling Two-Phase Fluid Flow with Two-Phase Darcy Flow in Anisotropic Porous Media

Directory of Open Access Journals (Sweden)

Jie Chen

2014-06-01

Full Text Available This paper reports a numerical study of coupling two-phase fluid flow in a free fluid region with two-phase Darcy flow in a homogeneous and anisotropic porous medium region. The model consists of coupled Cahn-Hilliard and Navier-Stokes equations in the free fluid region and the two-phase Darcy law in the anisotropic porous medium region. A Robin-Robin domain decomposition method is used for the coupled Navier-Stokes and Darcy system with the generalized Beavers-Joseph-Saffman condition on the interface between the free flow and the porous media regions. Obtained results have shown the anisotropic properties effect on the velocity and pressure of the two-phase flow.

A Two-Phase Solid/Fluid Model for Dense Granular Flows Including Dilatancy Effects

Science.gov (United States)

Mangeney, Anne; Bouchut, Francois; Fernandez-Nieto, Enrique; Narbona-Reina, Gladys

2015-04-01

We propose a thin layer depth-averaged two-phase model to describe solid-fluid mixtures such as debris flows. It describes the velocity of the two phases, the compression/dilatation of the granular media and its interaction with the pore fluid pressure, that itself modifies the friction within the granular phase (Iverson et al., 2010). The model is derived from a 3D two-phase model proposed by Jackson (2000) based on the 4 equations of mass and momentum conservation within the two phases. This system has 5 unknowns: the solid and fluid velocities, the solid and fluid pressures and the solid volume fraction. As a result, an additional equation inside the mixture is necessary to close the system. Surprisingly, this issue is inadequately accounted for in the models that have been developed on the basis of Jackson's work (Bouchut et al., 2014). In particular, Pitman and Le replaced this closure simply by imposing an extra boundary condition at the surface of the flow. When making a shallow expansion, this condition can be considered as a closure condition. However, the corresponding model cannot account for a dissipative energy balance. We propose here an approach to correctly deal with the thermodynamics of Jackson's equations. We close the mixture equations by a weak compressibility relation involving a critical density, or equivalently a critical pressure. Moreover, we relax one boundary condition, making it possible for the fluid to escape the granular media when compression of the granular mass occurs. Furthermore, we introduce second order terms in the equations making it possible to describe the evolution of the pore fluid pressure in response to the compression/dilatation of the granular mass without prescribing an extra ad-hoc equation for the pore pressure. We prove that the energy balance associated with this Jackson closure is dissipative, as well as its thin layer associated model. We present several numerical tests for the 1D case that are compared to the

Advanced Fluid Reduced Order Models for Compressible Flow.

Energy Technology Data Exchange (ETDEWEB)

Tezaur, Irina Kalashnikova [Sandia National Lab. (SNL-CA), Livermore, CA (United States); Fike, Jeffrey A. [Sandia National Lab. (SNL-NM), Albuquerque, NM (United States); Carlberg, Kevin Thomas [Sandia National Lab. (SNL-CA), Livermore, CA (United States); Barone, Matthew F. [Sandia National Lab. (SNL-NM), Albuquerque, NM (United States); Maddix, Danielle [Stanford Univ., CA (United States); Mussoni, Erin E. [Sandia National Lab. (SNL-CA), Livermore, CA (United States); Balajewicz, Maciej [Univ. of Illinois, Urbana-Champaign, IL (United States)

2017-09-01

This report summarizes fiscal year (FY) 2017 progress towards developing and implementing within the SPARC in-house finite volume flow solver advanced fluid reduced order models (ROMs) for compressible captive-carriage flow problems of interest to Sandia National Laboratories for the design and qualification of nuclear weapons components. The proposed projection-based model order reduction (MOR) approach, known as the Proper Orthogonal Decomposition (POD)/Least- Squares Petrov-Galerkin (LSPG) method, can substantially reduce the CPU-time requirement for these simulations, thereby enabling advanced analyses such as uncertainty quantification and de- sign optimization. Following a description of the project objectives and FY17 targets, we overview briefly the POD/LSPG approach to model reduction implemented within SPARC . We then study the viability of these ROMs for long-time predictive simulations in the context of a two-dimensional viscous laminar cavity problem, and describe some FY17 enhancements to the proposed model reduction methodology that led to ROMs with improved predictive capabilities. Also described in this report are some FY17 efforts pursued in parallel to the primary objective of determining whether the ROMs in SPARC are viable for the targeted application. These include the implemen- tation and verification of some higher-order finite volume discretization methods within SPARC (towards using the code to study the viability of ROMs on three-dimensional cavity problems) and a novel structure-preserving constrained POD/LSPG formulation that can improve the accuracy of projection-based reduced order models. We conclude the report by summarizing the key takeaways from our FY17 findings, and providing some perspectives for future work.

Optimization of a new flow design for solid oxide cells using computational fluid dynamics modelling

DEFF Research Database (Denmark)

Duhn, Jakob Dragsbæk; Jensen, Anker Degn; Wedel, Stig

2016-01-01

Design of a gas distributor to distribute gas flow into parallel channels for Solid Oxide Cells (SOC) is optimized, with respect to flow distribution, using Computational Fluid Dynamics (CFD) modelling. The CFD model is based on a 3d geometric model and the optimized structural parameters include...... the width of the channels in the gas distributor and the area in front of the parallel channels. The flow of the optimized design is found to have a flow uniformity index value of 0.978. The effects of deviations from the assumptions used in the modelling (isothermal and non-reacting flow) are evaluated...... and it is found that a temperature gradient along the parallel channels does not affect the flow uniformity, whereas a temperature difference between the channels does. The impact of the flow distribution on the maximum obtainable conversion during operation is also investigated and the obtainable overall...

A Comprehensive Prediction Model of Hydraulic Extended-Reach Limit Considering the Allowable Range of Drilling Fluid Flow Rate in Horizontal Drilling.

Science.gov (United States)

Li, Xin; Gao, Deli; Chen, Xuyue

2017-06-08

Hydraulic extended-reach limit (HERL) model of horizontal extended-reach well (ERW) can predict the maximum measured depth (MMD) of the horizontal ERW. The HERL refers to the well's MMD when drilling fluid cannot be normally circulated by drilling pump. Previous model analyzed the following two constraint conditions, drilling pump rated pressure and rated power. However, effects of the allowable range of drilling fluid flow rate (Q min  ≤ Q ≤ Q max ) were not considered. In this study, three cases of HERL model are proposed according to the relationship between allowable range of drilling fluid flow rate and rated flow rate of drilling pump (Q r ). A horizontal ERW is analyzed to predict its HERL, especially its horizontal-section limit (L h ). Results show that when Q min  ≤ Q r  ≤ Q max (Case I), L h depends both on horizontal-section limit based on rated pump pressure (L h1 ) and horizontal-section limit based on rated pump power (L h2 ); when Q min  drilling fluid flow rate, while L h2 keeps decreasing as the drilling fluid flow rate increases. The comprehensive model provides a more accurate prediction on HERL.

Squeeze flow of a Carreau fluid during sphere impact

KAUST Repository

Uddin, J.

2012-07-19

We present results from a combined numerical and experimental investigation into the squeeze flow induced when a solid sphere impacts onto a thin, ultra-viscous film of non-Newtonian fluid. We examine both the sphere motion through the liquid as well as the fluid flow field in the region directly beneath the sphere during approach to a solid plate. In the experiments we use silicone oil as the model fluid, which is well-described by the Carreau model. We use high-speed imaging and particle tracking to achieve flow visualisation within the film itself and derive the corresponding velocity fields. We show that the radial velocity either diverges as the gap between the sphere and the wall diminishes (Z tip → 0) or that it reaches a maximum value and then decays rapidly to zero as the sphere comes to rest at a non-zero distance (Z tip = Z min ) away from the wall. The horizontal shear rate is calculated and is responsible for significant viscosity reduction during the approach of the sphere. Our model of this flow, based on lubrication theory, is solved numerically and compared to experimental trials. We show that our model is able to correctly describe the physical features of the flow observed in the experiments.

Squeeze flow of a Carreau fluid during sphere impact

KAUST Repository

Uddin, J.; Marston, J. O.; Thoroddsen, Sigurdur T

2012-01-01

We present results from a combined numerical and experimental investigation into the squeeze flow induced when a solid sphere impacts onto a thin, ultra-viscous film of non-Newtonian fluid. We examine both the sphere motion through the liquid as well as the fluid flow field in the region directly beneath the sphere during approach to a solid plate. In the experiments we use silicone oil as the model fluid, which is well-described by the Carreau model. We use high-speed imaging and particle tracking to achieve flow visualisation within the film itself and derive the corresponding velocity fields. We show that the radial velocity either diverges as the gap between the sphere and the wall diminishes (Z tip → 0) or that it reaches a maximum value and then decays rapidly to zero as the sphere comes to rest at a non-zero distance (Z tip = Z min ) away from the wall. The horizontal shear rate is calculated and is responsible for significant viscosity reduction during the approach of the sphere. Our model of this flow, based on lubrication theory, is solved numerically and compared to experimental trials. We show that our model is able to correctly describe the physical features of the flow observed in the experiments.

Computational fluid dynamics modeling of two-phase flow in a BWR fuel assembly

International Nuclear Information System (INIS)

Andrey Ioilev; Maskhud Samigulin; Vasily Ustinenko; Simon Lo; Adrian Tentner

2005-01-01

Full text of publication follows: The goal of this project is to develop an advanced Computational Fluid Dynamics (CFD) computer code (CFD-BWR) that allows the detailed analysis of the two-phase flow and heat transfer phenomena in a Boiling Water Reactor (BWR) fuel bundle under various operating conditions. This code will include more fundamental physical models than the current generation of sub-channel codes and advanced numerical algorithms for improved computational accuracy, robustness, and speed. It is highly desirable to understand the detailed two-phase flow phenomena inside a BWR fuel bundle. These phenomena include coolant phase changes and multiple flow regimes which directly influence the coolant interaction with fuel assembly and, ultimately, the reactor performance. Traditionally, the best analysis tools for the analysis of two-phase flow phenomena inside the BWR fuel assembly have been the sub-channel codes. However, the resolution of these codes is still too coarse for analyzing the detailed intra-assembly flow patterns, such as flow around a spacer element. Recent progress in Computational Fluid Dynamics (CFD), coupled with the rapidly increasing computational power of massively parallel computers, shows promising potential for the fine-mesh, detailed simulation of fuel assembly two-phase flow phenomena. However, the phenomenological models available in the commercial CFD programs are not as advanced as those currently being used in the sub-channel codes used in the nuclear industry. In particular, there are no models currently available which are able to reliably predict the nature of the flow regimes, and use the appropriate sub-models for those flow regimes. The CFD-BWR code is being developed as a customized module built on the foundation of the commercial CFD Code STAR-CD which provides general two-phase flow modeling capabilities. The paper describes the model development strategy which has been adopted by the development team for the

Investigation of the fluid flow dynamic parameters for Newtonian and non-Newtonian materials: an approach to understanding the fluid flow-like structures within fault zones

Science.gov (United States)

Tanaka, H.; Shiomi, Y.; Ma, K.-F.

2017-11-01

To understand the fault zone fluid flow-like structure, namely the ductile deformation structure, often observed in the geological field (e.g., Ramsay and Huber The techniques of modern structure geology, vol. 1: strain analysis, Academia Press, London, 1983; Hobbs and Ord Structure geology: the mechanics of deforming metamorphic rocks, Vol. I: principles, Elsevier, Amsterdam, 2015), we applied a theoretical approach to estimate the rate of deformation, the shear stress and the time to form a streak-line pattern in the boundary layer of viscous fluids. We model the dynamics of streak lines in laminar boundary layers for Newtonian and pseudoplastic fluids and compare the results to those obtained via laboratory experiments. The structure of deformed streak lines obtained using our model is consistent with experimental observations, indicating that our model is appropriate for understanding the shear rate, flow time and shear stress based on the profile of deformed streak lines in the boundary layer in Newtonian and pseudoplastic viscous materials. This study improves our understanding of the transportation processes in fluids and of the transformation processes in fluid-like materials. Further application of this model could facilitate understanding the shear stress and time history of the fluid flow-like structure of fault zones observed in the field.[Figure not available: see fulltext.

Modeling Asymmetric Flow of Viscoelastic Fluid in Symmetric Planar Sudden Expansion Geometry Based on User-Defined Function in FLUENT CFD Package

Directory of Open Access Journals (Sweden)

Zhi-Ying Zheng

2013-01-01

Full Text Available Through embedding an in-house subroutine into FLUENT code by utilizing the functionalization of user-defined function provided by the software, a new numerical simulation methodology on viscoelastic fluid flows has been established. In order to benchmark this methodology, numerical simulations under different viscoelastic fluid solution concentrations (with solvent viscosity ratio varied from 0.2 to 0.9, extensibility parameters (100≤L2≤500, Reynolds numbers (0.1 ≤ Re ≤ 100, and Weissenberg numbers (0 ≤ Wi ≤ 20 are conducted on unsteady laminar flows through a symmetric planar sudden expansion with expansion ratio of 1: 3 for viscoelastic fluid flows. The constitutive model used to describe the viscoelastic effect of viscoelastic fluid flow is FENE-P (finitely extensive nonlinear elastic-Peterlin model. The numerical simulation results show that the influences of elasticity, inertia, and concentration on the flow bifurcation characteristics are more significant than those of extensibility. The present simulation results including the critical Reynolds number for which the flow becomes asymmetric, vortex size, bifurcation diagram, velocity distribution, streamline, and pressure loss show good agreements with some published results. That means the newly established method based on FLUENT software platform for simulating peculiar flow behaviors of viscoelastic fluid is credible and suitable for the study of viscoelastic fluid flows.

Computational modeling of unsteady third-grade fluid flow over a vertical cylinder: A study of heat transfer visualization

Science.gov (United States)

Reddy, G. Janardhana; Hiremath, Ashwini; Kumar, Mahesh

2018-03-01

The present paper aims to investigate the effect of Prandtl number for unsteady third-grade fluid flow over a uniformly heated vertical cylinder using Bejan's heat function concept. The mathematical model of this problem is given by highly time-dependent non-linear coupled equations and are resolved by an efficient unconditionally stable implicit scheme. The time histories of average values of momentum and heat transport coefficients as well as the steady-state flow variables are displayed graphically for distinct values of non-dimensional control parameters arising in the system. As the non-dimensional parameter value gets amplified, the time taken for the fluid flow variables to attain the time-independent state is decreasing. The dimensionless heat function values are closely associated with an overall rate of heat transfer. Thermal energy transfer visualization implies that the heat function contours are compact in the neighborhood of the leading edge of the hot cylindrical wall. It is noticed that the deviations of flow-field variables from the hot wall for a non-Newtonian third-grade fluid flow are significant compared to the usual Newtonian fluid flow.

Modeling and analysis of hydrodynamic instabilities in two-phase flow using two-fluid model

International Nuclear Information System (INIS)

Zhou, J.; Podowski, M.Z.

2001-01-01

Because of the practical importance of two-phase flow instabilities, especially in boiling water nuclear reactor technology, substantial efforts have been made to date to understand the physical phenomena governing such instabilities and to develop computational tools to model the dynamics of marginally-stable/unstable boiling systems. The purpose of this paper is to present an integrated methodology for the analysis of flow-induced instabilities in boiling channels and systems. The major novel aspects of the proposed approach are: (a) it is based on the combined frequency-domain and time-domain methods, the former used to quantify stability margins and to determine the onset of instability conditions, the latter to study the nonlinear system response outside the stability boundaries identified using the nearly-exact results of the frequency-domain analysis; (b) the two-fluid model of two-phase flow has been used for the first time to analytically derive the boiling channel transfer functions for the parallel-channel and channel-to-channel instability modes. In this way, the major characteristics of a boiling system, including the onset-of-instability conditions, can be readily evaluated by using the qualitative frequency-domain approach, whereas the explicit time-domain integration is performed, if necessary, only for the operating conditions that have already been identified as unstable. Both methods use the same physical two-fluid model that, in one case, is linearized and used to derive a rigorous analytical solution in the complex domain, and, in the other case, is solved numerically using an algorithm developed especially for this purpose. The results using both methods have been compared against each other and extensively tested. The testing and validation of the new model included comparisons of the predicted steady-state distributions of major parameters and of the transient channel response against experimental data

Slip analysis of squeezing flow using doubly stratified fluid

Science.gov (United States)

Ahmad, S.; Farooq, M.; Javed, M.; Anjum, Aisha

2018-06-01

The non-isothermal flow is modeled and explored for squeezed fluid. The influence of velocity, thermal and solutal slip effects on transport features of squeezed fluid are analyzed through Darcy porous channel when fluid is moving due to squeezing of upper plate towards the stretchable lower plate. Dual stratification effects are illustrated in transport equations. A similarity analysis is performed and reduced governing flow equations are solved using moderated and an efficient convergent approach i.e. Homotopic technique. The significant effects of physical emerging parameters on flow velocity, temperature and fluid concentration are reporting through various plots. Graphical explanations for drag force, Nusselt and Sherwood numbers are stated and examined. The results reveal that minimum velocity field occurs near the plate, whereas it increases far away from the plate for strong velocity slip parameter. Furthermore, temperature and fluid concentration significantly decreases with increased slip effects. The current analysis is applicable in some advanced technological processes and industrial fluid mechanics.

On Laminar Flow of Non-Newtonian Fluids in Porous Media

KAUST Repository

Fayed, Hassan E.

2015-10-20

Flow of generalized Newtonian fluids in porous media can be modeled as a bundle of capillary tubes or a pore-scale network. In general, both approaches rely on the solution of Hagen–Poiseuille equation using power law to estimate the variations in the fluid viscosity due to the applied shear rate. Despite the effectiveness and simplicity, power law tends to provide unrealistic values for the effective viscosity especially in the limits of zero and infinite shear rates. Here, instead of using power law, Carreau model (bubbles, drops, and particles in non-Newtonian fluids. Taylor & Francis Group, New York, 2007) is used to determine the effective viscosity as a function of the shear strain rate. Carreau model can predict accurately the variation in the viscosity at all shear rates and provide more accurate solution for the flow physics in a single pore. Using the results for a single pore, normalized Fanning friction coefficient has been calculated and plotted as a function of the newly defined Reynolds number based on pressure gradient. For laminar flow, the variation in the friction coefficient with Reynolds number has been plotted and scaled. It is observed that generalized Newtonian fluid flows show Newtonian nature up to a certain Reynolds number. At high Reynolds number, deviation from the Newtonian behavior is observed. The main contribution of this paper is to present a closed-form solution for the flow in a single pore using Carreau model, which allows for fast evaluation of the relationship between flux and pressure gradient in an arbitrary pore diameter. In this way, we believe that our development will open the perspectives for using Carreau models in pore-network simulations at low computational costs to obtain more accurate prediction for generalized Newtonian fluid flows in porous media.

On Laminar Flow of Non-Newtonian Fluids in Porous Media

KAUST Repository

Fayed, Hassan E.; Sheikh, Nadeem A.; Iliev, Oleg

2015-01-01

Flow of generalized Newtonian fluids in porous media can be modeled as a bundle of capillary tubes or a pore-scale network. In general, both approaches rely on the solution of Hagen–Poiseuille equation using power law to estimate the variations in the fluid viscosity due to the applied shear rate. Despite the effectiveness and simplicity, power law tends to provide unrealistic values for the effective viscosity especially in the limits of zero and infinite shear rates. Here, instead of using power law, Carreau model (bubbles, drops, and particles in non-Newtonian fluids. Taylor & Francis Group, New York, 2007) is used to determine the effective viscosity as a function of the shear strain rate. Carreau model can predict accurately the variation in the viscosity at all shear rates and provide more accurate solution for the flow physics in a single pore. Using the results for a single pore, normalized Fanning friction coefficient has been calculated and plotted as a function of the newly defined Reynolds number based on pressure gradient. For laminar flow, the variation in the friction coefficient with Reynolds number has been plotted and scaled. It is observed that generalized Newtonian fluid flows show Newtonian nature up to a certain Reynolds number. At high Reynolds number, deviation from the Newtonian behavior is observed. The main contribution of this paper is to present a closed-form solution for the flow in a single pore using Carreau model, which allows for fast evaluation of the relationship between flux and pressure gradient in an arbitrary pore diameter. In this way, we believe that our development will open the perspectives for using Carreau models in pore-network simulations at low computational costs to obtain more accurate prediction for generalized Newtonian fluid flows in porous media.

Thermohydrodynamic analysis of cryogenic liquid turbulent flow fluid film bearings

Science.gov (United States)

Andres, Luis San

1993-01-01

A thermohydrodynamic analysis is presented and a computer code developed for prediction of the static and dynamic force response of hydrostatic journal bearings (HJB's), annular seals or damper bearing seals, and fixed arc pad bearings for cryogenic liquid applications. The study includes the most important flow characteristics found in cryogenic fluid film bearings such as flow turbulence, fluid inertia, liquid compressibility and thermal effects. The analysis and computational model devised allow the determination of the flow field in cryogenic fluid film bearings along with the dynamic force coefficients for rotor-bearing stability analysis.

Fluid Flow in a Porous Tree-Shaped Network

OpenAIRE

Miguel, A. F.

2014-01-01

Tree-shaped flow networks connect one point to an infinity of points and are everywhere in Nature. These networks often own minimal flow resistance and vessel sizes obey to scaling power-laws. In this paper presents a model for fluid flow through a tree-shaped network with porous tubes. Hagen–Poiseuille flow is assumed for tubes and Darcy flow for the porous wall.

Computational fluid dynamics modeling of a lithium/thionyl chloride battery with electrolyte flow

Energy Technology Data Exchange (ETDEWEB)

Gu, W.B.; Wang, C.Y.; Weidner, J.W.; Jungst, R.G.; Nagasubramanian, G.

2000-02-01

A two-dimensional model is developed to simulate discharge of a lithium/thionyl chloride primary battery. As in earlier one-dimensional models, the model accounts for transport of species and charge, and electrode porosity variations and electrolyte flow induced by the volume reduction caused by electrochemical reactions. Numerical simulations are performed using a finite volume method of computational fluid dynamics. The predicted discharge curves for various temperatures show good agreement with published experimental data, and are essentially identical to results published for one-dimensional models. The detailed two-dimensional flow simulations show that the electrolyte is replenished from the cell head space predominantly through the separator into the front of the cathode during most parts of the discharge, especially for higher cell temperatures.

A two-phase debris-flow model that includes coupled evolution of volume fractions, granular dilatancy, and pore-fluid pressure

Science.gov (United States)

George, David L.; Iverson, Richard M.

2011-01-01

Pore-fluid pressure plays a crucial role in debris flows because it counteracts normal stresses at grain contacts and thereby reduces intergranular friction. Pore-pressure feedback accompanying debris deformation is particularly important during the onset of debrisflow motion, when it can dramatically influence the balance of forces governing downslope acceleration. We consider further effects of this feedback by formulating a new, depth-averaged mathematical model that simulates coupled evolution of granular dilatancy, solid and fluid volume fractions, pore-fluid pressure, and flow depth and velocity during all stages of debris-flow motion. To illustrate implications of the model, we use a finite-volume method to compute one-dimensional motion of a debris flow descending a rigid, uniformly inclined slope, and we compare model predictions with data obtained in large-scale experiments at the USGS debris-flow flume. Predictions for the first 1 s of motion show that increasing pore pressures (due to debris contraction) cause liquefaction that enhances flow acceleration. As acceleration continues, however, debris dilation causes dissipation of pore pressures, and this dissipation helps stabilize debris-flow motion. Our numerical predictions of this process match experimental data reasonably well, but predictions might be improved by accounting for the effects of grain-size segregation.

Integrated geophysical and hydrothermal models of flank degassing and fluid flow at Masaya Volcano, Nicaragua

Science.gov (United States)

Sanford, Ward E.; Pearson, S.C.P.; Kiyosugi, K.; Lehto, H.L.; Saballos, J.A.; Connor, C.B.

2012-01-01

We investigate geologic controls on circulation in the shallow hydrothermal system of Masaya volcano, Nicaragua, and their relationship to surface diffuse degassing. On a local scale (~250 m), relatively impermeable normal faults dipping at ~60° control the flowpath of water vapor and other gases in the vadose zone. These shallow normal faults are identified by modeling of a NE-SW trending magnetic anomaly of up to 2300 nT that corresponds to a topographic offset. Elevated SP and CO2 to the NW of the faults and an absence of CO2 to the SE suggest that these faults are barriers to flow. TOUGH2 numerical models of fluid circulation show enhanced flow through the footwalls of the faults, and corresponding increased mass flow and temperature at the surface (diffuse degassing zones). On a larger scale, TOUGH2 modeling suggests that groundwater convection may be occurring in a 3-4 km radial fracture zone transecting the entire flank of the volcano. Hot water rising uniformly into the base of the model at 1 x 10-5 kg/m2s results in convection that focuses heat and fluid and can explain the three distinct diffuse degassing zones distributed along the fracture. Our data and models suggest that the unusually active surface degassing zones at Masaya volcano can result purely from uniform heat and fluid flux at depth that is complicated by groundwater convection and permeability variations in the upper few km. Therefore isolating the effects of subsurface geology is vital when trying to interpret diffuse degassing in light of volcanic activity.

Fluid flow solidification simulation of molten alloys

International Nuclear Information System (INIS)

Kaschnitz, E.

1997-01-01

In an effort to minimize costs and to obtain optimum designs, computer simulation of shape casting processes is more and more used as a development tool. Accurate predictions are possible by means of three dimensional fluid flow and solidification modelling. The bases of the model are the transient laminar Navier-Stokes-equations for a Newtonian fluid including the tracking of the free surface. They are describing the melt flow pattern during the mold filling sequence. Simultaneously, the temperature development in the alloy and mold is calculated using Fourier's heat transfer equation. At OEGI, a commercial software package (MAGMAsoft) with a finite difference equation solver is used for improvement of casting processes. Different examples of industrial applications will be shown. (author)

A two-dimensional continuum model of biofilm growth incorporating fluid flow and shear stress based detachment

KAUST Repository

Duddu, Ravindra

2009-05-01

We present a two-dimensional biofilm growth model in a continuum framework using an Eulerian description. A computational technique based on the eXtended Finite Element Method (XFEM) and the level set method is used to simulate the growth of the biofilm. The model considers fluid flow around the biofilm surface, the advection-diffusion and reaction of substrate, variable biomass volume fraction and erosion due to the interfacial shear stress at the biofilm-fluid interface. The key assumptions of the model and the governing equations of transport, biofilm kinetics and biofilm mechanics are presented. Our 2D biofilm growth results are in good agreement with those obtained by Picioreanu et al. (Biotechnol Bioeng 69(5):504-515, 2000). Detachment due to erosion is modeled using two continuous speed functions based on: (a) interfacial shear stress and (b) biofilm height. A relation between the two detachment models in the case of a 1D biofilm is established and simulated biofilm results with detachment in 2D are presented. The stress in the biofilm due to fluid flow is evaluated and higher stresses are observed close to the substratum where the biofilm is attached. © 2008 Wiley Periodicals, Inc.

Development of bubble-induced turbulence model for advanced two-fluid model

International Nuclear Information System (INIS)

Hosoi, Hideaki; Yoshida, Hiroyuki

2011-01-01

A two-fluid model can simulate two-phase flow by computational cost less than detailed two-phase flow simulation method such as interface tracking method. The two-fluid model is therefore useful for thermal hydraulic analysis in the large-scale domain such as rod bundles. However, since the two-fluid model includes a lot of constitutive equations verified by use of experimental results, it has problems that the result of analyses depends on accuracy of the constitutive equations. To solve these problems, an advanced two-fluid model has been developed by Japan Atomic Energy Agency. In this model, interface tracking method is combined with two-fluid model to accurately predict large interface structure behavior. Liquid clusters and bubbles larger than a computational cell are calculated using the interface tracking method, and those smaller than the cell are simulated by the two-fluid model. The constitutive equations to evaluate the effects of small bubbles or droplets on two-phase flow are also required in the advanced two-fluid model, just as with the conventional two-fluid model. However, the dependency of small bubbles and droplets on two-phase flow characteristics is relatively small, and fewer experimental results are required to verify the characteristics of large interface structures. Turbulent dispersion force model is one of the most important constitutive equations for the advanced two-fluid model. The turbulent dispersion force model has been developed by many researchers for the conventional two-fluid model. However, existing models implicitly include the effects of large bubbles and the deformation of bubbles, and are unfortunately not applicable to the advanced two-fluid model. In the previous study, the authors suggested the turbulent dispersion force model based on the analogy of Brownian motion. And the authors improved the turbulent dispersion force model in consideration of bubble-induced turbulence to improve the analysis results for small

Numerical simulation of travelling wave induced electrothermal fluid flow

International Nuclear Information System (INIS)

Perch-Nielsen, Ivan R; Green, Nicolas G; Wolff, Anders

2004-01-01

Many microdevices for manipulating particles and cells use electric fields to produce a motive force on the particles. The movement of particles in non-uniform electric fields is called dielectrophoresis, and the usual method of applying this effect is to pass the particle suspension over a microelectrode structure. If the suspension has a noticeable conductivity, one important side effect is that the electric field drives a substantial conduction current through the fluid, causing localized Joule-heating. The resulting thermal gradient produces local conductivity and permittivity changes in the fluid. Dielectrophoretic forces acting upon these pockets of fluid will then produce motion of both the fluid and the particles. This paper presents a numerical solution of the electrical force and the resulting electrothermal driven fluid flow on a travelling wave structure. This common electrode geometry consists of interdigitated electrodes laid down in a long array, with the phase of the applied potential shifted by 90 0 on each subsequent electrode. The resulting travelling electric field was simulated and the thermal field and electrical body force on the fluid calculated, for devices constructed from two typical materials: silicon and glass. The electrothermal fluid flow in the electrolyte over the electrode array was then numerically simulated. The model predicts that the thermal field depends on the conductivity and applied voltage, but more importantly on the geometry of the system and the material used in the construction of the device. The velocity of the fluid flow depends critically on the same parameters, with slight differences in the thermal field for glass and silicon leading to diametrically opposite flow direction with respect to the travelling field for the two materials. In addition, the imposition of slight external temperature gradients is shown to have a large effect on the fluid flow in the device, under certain conditions leading to a reversal of

Acoustic geometry for general relativistic barotropic irrotational fluid flow

International Nuclear Information System (INIS)

Visser, Matt; Molina-ParIs, Carmen

2010-01-01

'Acoustic spacetimes', in which techniques of differential geometry are used to investigate sound propagation in moving fluids, have attracted considerable attention over the last few decades. Most of the models currently considered in the literature are based on non-relativistic barotropic irrotational fluids, defined in a flat Newtonian background. The extension, first to special relativistic barotropic fluid flow and then to general relativistic barotropic fluid flow in an arbitrary background, is less straightforward than it might at first appear. In this paper, we provide a pedagogical and simple derivation of the general relativistic 'acoustic spacetime' in an arbitrary (d+1)-dimensional curved-space background.

Gas-liquid Two Phase Flow Modelling of Incompressible Fluid and Experimental Validation Studies in Vertical Centrifugal Casting

International Nuclear Information System (INIS)

Zhou, J X; Shen, X; Yin, Y J; Guo, Z; Wang, H

2015-01-01

In this paper, Gas-liquid two phase flow mathematic models of incompressible fluid were proposed to explore the feature of fluid under certain centrifugal force in vertical centrifugal casting (VCC). Modified projection-level-set method was introduced to solve the mathematic models. To validate the simulation results, two methods were used in this study. In the first method, the simulation result of basic VCC flow process was compared with its analytic solution. The relationship between the numerical solution and deterministic analytic solution was presented to verify the correctness of numerical algorithms. In the second method, systematic water simulation experiments were developed. In this initial experiment, special experimental vertical centrifugal device and casting shapes were designed to describe typical mold-filling processes in VCC. High speed camera system and data collection devices were used to capture flow shape during the mold-filling process. Moreover, fluid characteristic at different rotation speed (from 40rpm, 60rpmand 80rpm) was discussed to provide comparative resource for simulation results. As compared with the simulation results, the proposed mathematical models could be proven and the experimental design could help us advance the accuracy of simulation and further studies for VCC. (paper)

Flask fluid flow simulation using CFD

International Nuclear Information System (INIS)

Swindlehurst, W.E.; Livesey, E.; Worthington, D.

1989-01-01

BNFL and its subsidiary Company, PNTL, design and operate waterfilled LWR fuel transport flasks for the international transport of irradiated fuel. Although some 150 flasks are currently in operation, new flask designs are being developed. As part of the supporting R and D program, Computational Fluid Dynamics (CFD) codes are being investigated as a means of predicting fluid movements and temperatures within the complex internal geometry of flasks. The ability to simulate fluid flow is particularly important when convection heat transfer is significant. Although obviously relevant to water filled flasks, the technique is applicable to dry flask thermal assessments (where experience shows that convection heat transfer is often underestimated). Computational Fluid Dynamics has emerged in recent years as an important technique in engineering design and safety assessments. Cheaper computing and the development of general CFD codes allows complex engineering structures to be analyzed. However, because of this complexity, it is essential that the application and associated modeling assumptions are critically reviewed. To assess the ability of a CFD code to model flask internals, the code PHOENICS has been used to model the fluid movements in a BNFL Excellox-type flask and the results compared with test data

Model validation and parametric study of fluid flows and heat transfer of aviation kerosene with endothermic pyrolysis at supercritical pressure

Directory of Open Access Journals (Sweden)

Keke Xu

2015-12-01

Full Text Available The regenerative cooling technology is a promising approach for effective thermal protection of propulsion and power-generation systems. A mathematical model has been used to examine fluid flows and heat transfer of the aviation kerosene RP-3 with endothermic fuel pyrolysis at a supercritical pressure of 5 MPa. A pyrolytic reaction mechanism, which consists of 18 species and 24 elementary reactions, is incorporated to account for fuel pyrolysis. Detailed model validations are conducted against a series of experimental data, including fluid temperature, fuel conversion rate, various product yields, and chemical heat sink, fully verifying the accuracy and reliability of the model. Effects of fuel pyrolysis and inlet flow velocity on flow dynamics and heat transfer characteristics of RP-3 are investigated. Results reveal that the endothermic fuel pyrolysis significantly improves the heat transfer process in the high fluid temperature region. During the supercritical-pressure heat transfer process, the flow velocity significantly increases, caused by the drastic variations of thermophysical properties. Under all the tested conditions, the Nusselt number initially increases, consistent with the increased flow velocity, and then slightly decreases in the high fluid temperature region, mainly owing to the decreased heat absorption rate from the endothermic pyrolytic chemical reactions.

Fluid flow in a porous medium with transverse permeability discontinuity

Science.gov (United States)

Pavlovskaya, Galina E.; Meersmann, Thomas; Jin, Chunyu; Rigby, Sean P.

2018-04-01

Magnetic resonance imaging (MRI) velocimetry methods are used to study fully developed axially symmetric fluid flow in a model porous medium of cylindrical symmetry with a transverse permeability discontinuity. Spatial mapping of fluid flow results in radial velocity profiles. High spatial resolution of these profiles allows estimating the slip in velocities at the boundary with a permeability discontinuity zone in a sample. The profiles are compared to theoretical velocity fields for a fully developed axially symmetric flow in a cylinder derived from the Beavers-Joseph [G. S. Beavers and D. D. Joseph, J. Fluid Mech. 30, 197 (1967), 10.1017/S0022112067001375] and Brinkman [H. C. Brinkman, Appl. Sci. Res. A 1, 27 (1947), 10.1007/BF02120313] models. Velocity fields are also computed using pore-scale lattice Boltzmann modeling (LBM) where the assumption about the boundary could be omitted. Both approaches give good agreement between theory and experiment, though LBM velocity fields follow the experiment more closely. This work shows great promise for MRI velocimetry methods in addressing the boundary behavior of fluids in opaque heterogeneous porous media.

Geometric analysis of the solutions of two-phase flows: two-fluid model

International Nuclear Information System (INIS)

Kestin, J.; Zeng, D.L.

1984-01-01

This report contains a lightly edited draft of a study of the two-fluid model in two-phase flow. The motivation for the study stems from the authors' conviction that the construction of a computer code for any model should be preceded by a geometrical analysis of the pattern of trajectories in the phase space appropriate for the model. Such a study greatly facilitates the understanding of the phenomenon of choking and anticipates the computational difficulties which arise from the existence of singularities. The report contains a derivation of the six conservation equations of the model which includes a consideration of the simplifications imposed on a one-dimensional treatment by the presence of boundary layers at the wall and between the phases. The model is restricted to one-dimensional adiabatic flows of a single substance present in two phases, but thermodynamic equilibrium between the phases is not assumed. The role of closure conditions is defined but no specific closure conditions, or explicit equations of state, are introduced

Interfacial Fluid Mechanics A Mathematical Modeling Approach

CERN Document Server

Ajaev, Vladimir S

2012-01-01

Interfacial Fluid Mechanics: A Mathematical Modeling Approach provides an introduction to mathematical models of viscous flow used in rapidly developing fields of microfluidics and microscale heat transfer. The basic physical effects are first introduced in the context of simple configurations and their relative importance in typical microscale applications is discussed. Then,several configurations of importance to microfluidics, most notably thin films/droplets on substrates and confined bubbles, are discussed in detail.  Topics from current research on electrokinetic phenomena, liquid flow near structured solid surfaces, evaporation/condensation, and surfactant phenomena are discussed in the later chapters. This book also:  Discusses mathematical models in the context of actual applications such as electrowetting Includes unique material on fluid flow near structured surfaces and phase change phenomena Shows readers how to solve modeling problems related to microscale multiphase flows Interfacial Fluid Me...

Heat and fluid flow during rapid solidification of non-equilibrium materials

International Nuclear Information System (INIS)

Negli, S.C.; Eddingfield, D.L.; Brower, W.E. Jr.

1990-01-01

Rapid solidification technology (RST) is an advanced solidification process which is being utilized to produce non-equilibrium structures with properties not previously available with conventionally cast materials. An iron based alloy rapidly quenched to form a metallic glass is being installed on a large scale in electric power transformers where it cuts heat losses dramatically. The formation of a non-equilibrium structure usually requires a cooling rate of at least a million degrees per second. Achieving this high a cooling rate depends not only on the heat transfer conditions during the quenching process, but also on the fluid flow conditions in the molten metal before and during solidification. This paper presents a model of both heat and fluid flow during RST by the hammer and anvil method. The symmetry of two sided cooling permits analysis which is still applicable to the one sided cooling that occurs during melt spinning, the prevalent method of RST. The heat flow is modeled as one dimensional, normal to the quench surface. Previous models have shown the heat flow in the plane of the quench surface not to be significant. The fluid flow portion of the model utilizes the squeeze film solution for flow between two parallel flat plates. The model predicts the effects of superheat of the melt and of the quench hammer speed upon cooling rate during the formation of nonequilibrium phases. An unexpected result is that increased superheat results in much higher cooling rates, due to fluid flow before a potential transformation would take place; this enhanced liquid metal flow results in a thinner section casting which in turn has a dominant effect on the cooling rate. The model also predicts an expanded regime of Newtonian (interface controlled) cooling by about a factor of ten as compared to previous model of RST

Heat transfer and fluid flow in biological processes advances and applications

CERN Document Server

Becker, Sid

2015-01-01

Heat Transfer and Fluid Flow in Biological Processes covers emerging areas in fluid flow and heat transfer relevant to biosystems and medical technology. This book uses an interdisciplinary approach to provide a comprehensive prospective on biofluid mechanics and heat transfer advances and includes reviews of the most recent methods in modeling of flows in biological media, such as CFD. Written by internationally recognized researchers in the field, each chapter provides a strong introductory section that is useful to both readers currently in the field and readers interested in learning more about these areas. Heat Transfer and Fluid Flow in Biological Processes is an indispensable reference for professors, graduate students, professionals, and clinical researchers in the fields of biology, biomedical engineering, chemistry and medicine working on applications of fluid flow, heat transfer, and transport phenomena in biomedical technology. Provides a wide range of biological and clinical applications of fluid...

Fluid flow behaviour of gas-condensate and near-miscible fluids at the pore scale

Energy Technology Data Exchange (ETDEWEB)

Dawe, Richard A. [Department of Chemical Engineering, University of West Indies, St. Augustine (Trinidad and Tobago); Grattoni, Carlos A. [Department of Earth Science and Engineering, Imperial College, London, SW7 2BP (United Kingdom)

2007-02-15

Retrograde condensate reservoir behaviour is complex with much of the detailed mechanisms of the multiphase fluid transport and mass transfer between the phases within the porous matrix still speculative. Visual modelling of selected processes occurring at the pore level under known and controlled boundary conditions can give an insight to fluid displacements at the core scale and help the interpretation of production behaviour at reservoir scale. Visualisation of the pore scale two-phase flow mechanisms has been studied experimentally at low interfacial tensions, < 0.5 mN/m, using a partially miscible fluid system in glass visual micro models. As the interfacial tension decreases the balance between fluid-fluid forces (interfacial, spreading and viscous) and fluid-solid interactions (wettability and viscous interactions) changes. Data measurements in the laboratory, particularly relative permeability, will therefore always be difficult especially for condensate fluids just below their dew point. What is certain is that gas production from a gas-condensate leads to condensate dropout when pressure falls below the dew point, either within the wellbore or, more importantly, in the reservoir. This paper illustrates some pore scale physics, particularly interfacial phenomena at low interfacial tension, which has relevance to appreciating the flow of condensate fluids close to their dew point either near the wellbore (which affects well productivity) or deep inside the reservoir (which affects condensate recovery). (author)

The influence of tip clearance on performance and internal flow condition of fluid food pump using low viscous fluid

International Nuclear Information System (INIS)

Kubo, S; Ishioka, T; Fukutomi, J; Shigemitsu, T

2012-01-01

Fluid machines for fluid food have been used in wide variety of fields i.e. transportation, the filling, and for the improvement of quality of fluid foods. However, flow conditions of it are quite complicated because fluid foods are different from water. Therefore, design methods based on internal flow conditions have not been conducted. In this research, turbo-pumps having a small number of blades were used to decrease shear loss and keep wide flow passage. The influence of the tip clearance was investigated by the numerical analysis using the model with and without the tip clearance. In this paper, the influence of tip clearance on performances and internal flow conditions of turbo-pump using low viscous fluid were clarified by experimental and numerical analysis results. In addition, design methods based on the internal flow were considered. Further, the influences of viscosity on the performance characteristic and internal flow were investigated.

Numerical study on modeling of liquid film flow under countercurrent flow limitation in volume of fluid method

Energy Technology Data Exchange (ETDEWEB)

Watanabe, Taro, E-mail: watanabe_t@qe.see.eng.osaka-u.ac.jp [Graduate School of Engineering, Osaka University, 2-1, Yamadaoka, Suita-shi, Osaka 565-7895 (Japan); Takata, Takashi, E-mail: takata.takashi@jaea.go.jp [Japan Atomic Energy Agency, 4002 Narita-chou, Oarai-machi, Higashi-Ibaraki-gun, Ibaraki 331-1393 (Japan); Yamaguchi, Akira, E-mail: yamaguchi@n.t.u-tokyo.ac.jp [Graduate School of Engineering, The University of Tokyo, 2-22 Shirakata-Shirane, Tokai-mura, Naka-gun, Ibaraki 319-1188 (Japan)

2017-03-15

Highlights: • Thin liquid film flow under CCFL was modeled and coupled with the VOF method. • The difference of the liquid flow rate in experiments of CCFL was evaluated. • The proposed VOF method can quantitatively predict CCFL with low computational cost. - Abstract: Countercurrent flow limitation (CCFL) in a heat transfer tube at a steam generator (SG) of pressurized water reactor (PWR) is one of the important issues on the core cooling under a loss of coolant accident (LOCA). In order to improve the prediction accuracy of the CCFL characteristics in numerical simulations using the volume of fluid (VOF) method with less computational cost, a thin liquid film flow in a countercurrent flow is modeled independently and is coupled with the VOF method. The CCFL characteristics is evaluated analytically in condition of a maximizing down-flow rate as a function of a void fraction or a liquid film thickness considering a critical thickness. Then, we have carried out numerical simulations of a countercurrent flow in a vertical tube so as to investigate the CCFL characteristics and compare them with the previous experimental results. As a result, it has been concluded that the effect of liquid film entrainment by upward gas flux will cause the difference in the experiments.

Mathematical modeling and exact solutions to rotating flows of a Burgers' fluid

International Nuclear Information System (INIS)

Hayat, T.

2005-12-01

The aim of this study is to provide the modeling and exact analytic solutions for hydromagnetic oscillatory rotating flows of an incompressible Burgers' fluid bounded by a plate. The governing time-dependent equation for the Burgers' fluid is different than those from the Navier-Stokes' equation. The entire system is assumed to rotate around an axis normal to the plate. The governing equations for this investigation are solved analytically for two physical problems. The solutions for the three cases, when the two times angular velocity is greater than the frequency of oscillation or it is smaller than the frequency or it is equal to the frequency (resonant case), are discussed in second problem. In Burgers' fluid, it is also found that hydromagnetic solution in the resonant case satisfies the boundary condition at infinity. Moreover, the obtained analytical results reduce to several previously published results as the special cases. (author)

Efficient simulations of fluid flow coupled with poroelastic deformations in pleated filters

KAUST Repository

Calo, Victor M.; Iliev, Dimitar; Iliev, Oleg; Kirsch, Ralf; Lakdawala, Zahra; Printsypar, Galina

2015-01-01

model describes a free fluid flow coupled with a flow in porous media in a domain that contains the filtering media. To discretize the complex computational domain we use quadrilateral boundary fitted grids which resolve porous-fluid interfaces

Fluids in crustal deformation: Fluid flow, fluid-rock interactions, rheology, melting and resources

Science.gov (United States)

Lacombe, Olivier; Rolland, Yann

2016-11-01

Fluids exert a first-order control on the structural, petrological and rheological evolution of the continental crust. Fluids interact with rocks from the earliest stages of sedimentation and diagenesis in basins until these rocks are deformed and/or buried and metamorphosed in orogens, then possibly exhumed. Fluid-rock interactions lead to the evolution of rock physical properties and rock strength. Fractures and faults are preferred pathways for fluids, and in turn physical and chemical interactions between fluid flow and tectonic structures, such as fault zones, strongly influence the mechanical behaviour of the crust at different space and time scales. Fluid (over)pressure is associated with a variety of geological phenomena, such as seismic cycle in various P-T conditions, hydrofracturing (including formation of sub-horizontal, bedding-parallel veins), fault (re)activation or gravitational sliding of rocks, among others. Fluid (over)pressure is a governing factor for the evolution of permeability and porosity of rocks and controls the generation, maturation and migration of economic fluids like hydrocarbons or ore forming hydrothermal fluids, and is therefore a key parameter in reservoir studies and basin modeling. Fluids may also help the crust partially melt, and in turn the resulting melt may dramatically change the rheology of the crust.

Vertically aligned carbon nanotubes for sensing unidirectional fluid flow

Energy Technology Data Exchange (ETDEWEB)

Kiani, Keivan, E-mail: k_kiani@kntu.ac.ir

2015-05-15

From applied mechanics points of view, potential application of ensembles of single-walled carbon nanotubes (SWCNTs) as fluid flow sensors is aimed to be examined. To this end, useful nonlocal analytical and numerical models are developed. The deflection of the ensemble of SWCNTs at the tip is introduced as a measure of its sensitivity. The influences of the length and radius of the SWCNT, intertube distance, fluid flow velocity, and distance of the ensemble from the leading edge of the rigid base on the deflection field of the ensemble are comprehensively examined. The obtained results display how calibration of an ensemble of SWCNTs can be methodically carried out in accordance with the characteristics of the ensemble and the external fluid flow.

Numerical computation of fluid flow in different nonferrous metallurgical reactors

International Nuclear Information System (INIS)

Lackner, A.

1996-10-01

Heat, mass and fluid flow phenomena in metallurgical reactor systems such as smelting cyclones or electrolytic cells are complex and intricately linked through the governing equations of fluid flow, chemical reaction kinetics and chemical thermodynamics. The challenges for the representation of flow phenomena in such reactors as well as the transfers of these concepts to non-specialist modelers (e.g. plant operators and management personnel) can be met through scientific flow visualization techniques. In the first example the fluid flow of the gas phase and of concentrate particles in a smelting cyclone for copper production are calculated three dimensionally. The effect of design parameters (length and diameter of reactor, concentrate feeding tangentially or from the top, ..) and operating conditions are investigated. Single particle traces show, how to increase particle retention time before the particles reach the liquid film flowing down the cyclone wall. Cyclone separators are widely used in the metallurgical and chemical industry for collection of large quantities of dust. Most of the empirical models, which today are applied for the design, are lacking in being valid in the high temperature region. Therefore the numerical prediction of the collection efficiency of dust particles is done. The particle behavior close to the wall is considered by applying a particle restitution model, which calculates individual particle restitution coefficients as functions of impact velocity and impact angle. The effect of design parameters and operating are studied. Moreover, the fluid flow inside a copper refining electrolysis cell is modeled. The simulation is based on density variations in the boundary layer at the electrode surface. Density and thickness of the boundary layer are compared to measurements in a parametric study. The actual inhibitor concentration in the cell is calculated, too. Moreover, a two-phase flow approach is developed to simulate the behavior of

The use of paleo-thermo-barometers and coupled thermal, fluid flow and pore fluid pressure modelling for hydrocarbon and reservoir prediction in fold and thrust belts

NARCIS (Netherlands)

Roure, F.; Andriessen, P.A.M.; Callot, J.P.; Ferket, H.; Gonzales, E.; Guilhaumou, N.; Hardebol, N.J.; Lacombe, O.; Malandain, J.; Mougin, P.; Muska, K.; Ortuno, S.; Sassi, W.; Swennen, R.; Vilasi, N.

2010-01-01

Basin modelling tools are now more efficient to reconstruct palinspastic structural cross sections and compute the history of temperature, pore-fluid pressure and fluid flow circulations in complex structural settings. In many cases and especially in areas where limited erosion occurred, the use of

Fluid dynamics characterization of riser in a FCC cold flow model using gas radiotracer

International Nuclear Information System (INIS)

Santos, Valdemir A. dos; Lima, Emerson A.O.

2013-01-01

Was carried out the characterization of a diameter small riser of a cold flow model of a circulating fluidized bed (CFB), with aid of a radioactive tracer. Compressed air and catalytic cracking of petroleum flow through solids pneumatic transport regime, made of transparent material (glass, acrylic, PVC, polycarbonate) for study of problems in Fluid Catalytic Cracking (FCC) unit and development of methods of measurement of fluid dynamic parameters. The CFB model consisted of a mixer component solid-gas (compressed air at 25 deg C and 200 kN/m 2 ; cracking catalyst with an average diameter of 72μm and specific mass of 1,500 kg/m 3 ), comprising a riser pipe glass 0.02m internal diameter and 1.8m height, a gas solid separation vessel by flash effect, with the filter in the gas outlet, and a return column (a glass tube with an internal diameter of 0.0254m) to redirect the catalyst for the riser base. Recorded data allowed studies on residence time distribution of the gaseous phase in the riser, with the identification and characterization of the flow of gas-solid components in the CFB riser of small diameter. A plug flow type with deviations due to back mixing of catalyst close to the walls, associated with the density difference between this component was observed. (author)

Numerical Modeling and Investigation of Fluid-Driven Fracture Propagation in Reservoirs Based on a Modified Fluid-Mechanically Coupled Model in Two-Dimensional Particle Flow Code

Directory of Open Access Journals (Sweden)

Jian Zhou

2016-09-01

Full Text Available Hydraulic fracturing is a useful tool for enhancing rock mass permeability for shale gas development, enhanced geothermal systems, and geological carbon sequestration by the high-pressure injection of a fracturing fluid into tight reservoir rocks. Although significant advances have been made in hydraulic fracturing theory, experiments, and numerical modeling, when it comes to the complexity of geological conditions knowledge is still limited. Mechanisms of fluid injection-induced fracture initiation and propagation should be better understood to take full advantage of hydraulic fracturing. This paper presents the development and application of discrete particle modeling based on two-dimensional particle flow code (PFC2D. Firstly, it is shown that the modeled value of the breakdown pressure for the hydraulic fracturing process is approximately equal to analytically calculated values under varied in situ stress conditions. Furthermore, a series of simulations for hydraulic fracturing in competent rock was performed to examine the influence of the in situ stress ratio, fluid injection rate, and fluid viscosity on the borehole pressure history, the geometry of hydraulic fractures, and the pore-pressure field, respectively. It was found that the hydraulic fractures in an isotropic medium always propagate parallel to the orientation of the maximum principal stress. When a high fluid injection rate is used, higher breakdown pressure is needed for fracture propagation and complex geometries of fractures can develop. When a low viscosity fluid is used, fluid can more easily penetrate from the borehole into the surrounding rock, which causes a reduction of the effective stress and leads to a lower breakdown pressure. Moreover, the geometry of the fractures is not particularly sensitive to the fluid viscosity in the approximate isotropic model.

A general coarse and fine mesh solution scheme for fluid flow modeling in VHTRS

International Nuclear Information System (INIS)

Clifford, I; Ivanov, K; Avramova, M.

2011-01-01

Coarse mesh Computational Fluid Dynamics (CFD) methods offer several advantages over traditional coarse mesh methods for the safety analysis of helium-cooled graphite-moderated Very High Temperature Reactors (VHTRs). This relatively new approach opens up the possibility for system-wide calculations to be carried out using a consistent set of field equations throughout the calculation, and subsequently the possibility for hybrid coarse/fine mesh or hierarchical multi scale CFD simulations. To date, a consistent methodology for hierarchical multi-scale CFD has not been developed. This paper describes work carried out in the initial development of a multi scale CFD solver intended to be used for the safety analysis of VHTRs. The VHTR is considered on any scale to consist of a homogenized two-phase mixture of fluid and stationary solid material of varying void fraction. A consistent set of conservation equations was selected such that they reduce to the single-phase conservation equations for the case where void fraction is unity. The discretization of the conservation equations uses a new pressure interpolation scheme capable of capturing the discontinuity in pressure across relatively large changes in void fraction. Based on this, a test solver was developed which supports fully unstructured meshes for three-dimensional time-dependent compressible flow problems, including buoyancy effects. For typical VHTR flow phenomena the new solver shows promise as an effective candidate for predicting the flow behavior on multiple scales, as it is capable of modeling both fine mesh single phase flows as well as coarse mesh flows in homogenized regions containing both fluid and solid materials. (author)

Modeling of turbulent flows in porous media and at the interface with a free fluid medium

International Nuclear Information System (INIS)

Chandesris, M.

2006-12-01

This work deals with the numerical simulation of turbulent flows in the whole nuclear reactor core, using multi-scale approaches. First, a macroscopic turbulence model is built, based on a porous media approach, to describe the flow in the fuel assemblies part of the nuclear core. Then, we study the jump conditions that have to be applied at a free fluid/porous interface. A thorough analytical study is carried out for laminar flows. This study allows to answer some fundamental questions about the physical meaning of the jump conditions, the values of the jump parameters and the location of the interface. Using these results, jump conditions for turbulent flows are proposed. The model is then applied to the simulation of a turbulent flow in a simplified model of a reactor core. (author)

Transient flows of a Burgers' fluid

International Nuclear Information System (INIS)

Khan, M.

2005-12-01

An analysis is performed to develop the analytical solutions for some unsteady magnetohydrodynamic (MHD) flows of a Burgers' fluid between two plates. A uniform magnetic field is applied transversely to the fluid motion. The exact solutions are given for three problems. Results for the velocity fields are discussed and compared with the flows of Oldroyd-B, Maxwell, second grade and Newtonian fluids. (author)

Statistically derived conservation equations for fluid particle flows

International Nuclear Information System (INIS)

Reyes, J.N. Jr.

1989-01-01

The behavior of water droplets in a heated nuclear fuel channel is of significant interest to nuclear reactor safety studies pertaining to loss-of-coolant accidents. This paper presents the derivation of the mass, momentum, and energy conservation equations for a distribution of fluid particles (bubbles or droplets) transported by a continuous fluid medium. When coupled with the appropriate closure equations, the conservation equations can be used to model nonequilibrium, two-phase, dispersed, fluid flow behavior

An investigation of dentinal fluid flow in dental pulp during food mastication: simulation of fluid-structure interaction.

Science.gov (United States)

Su, Kuo-Chih; Chuang, Shu-Fen; Ng, Eddie Yin-Kwee; Chang, Chih-Han

2014-06-01

This study uses fluid-structure interaction (FSI) simulation to investigate the relationship between the dentinal fluid flow in the dental pulp of a tooth and the elastic modulus of masticated food particles and to investigate the effects of chewing rate on fluid flow in the dental pulp. Three-dimensional simulation models of a premolar tooth (enamel, dentine, pulp, periodontal ligament, cortical bone, and cancellous bone) and food particle were created. Food particles with elastic modulus of 2,000 and 10,000 MPa were used, respectively. The external displacement loading (5 μm) was gradually directed to the food particle surface for 1 and 0.1 s, respectively, to simulate the chewing of food particles. The displacement and stress on tooth structure and fluid flow in the dental pulp were selected as evaluation indices. The results show that masticating food with a high elastic modulus results in high stress and deformation in the tooth structure, causing faster dentinal fluid flow in the pulp in comparison with that obtained with soft food. In addition, fast chewing of hard food particles can induce faster fluid flow in the pulp, which may result in dental pain. FSI analysis is shown to be a useful tool for investigating dental biomechanics during food mastication. FSI simulation can be used to predict intrapulpal fluid flow in dental pulp; this information may provide the clinician with important concept in dental biomechanics during food mastication.

Mechanics of fluid flow

CERN Document Server

Basniev, Kaplan S; Chilingar, George V 0

2012-01-01

The mechanics of fluid flow is a fundamental engineering discipline explaining both natural phenomena and human-induced processes, and a thorough understanding of it is central to the operations of the oil and gas industry.  This book, written by some of the world's best-known and respected petroleum engineers, covers the concepts, theories, and applications of the mechanics of fluid flow for the veteran engineer working in the field and the student, alike.  It is a must-have for any engineer working in the oil and gas industry.

Modelling of fluid flow phenomenon in laser+GMAW hybrid welding of aluminum alloy considering three phase coupling and arc plasma shear stress

Science.gov (United States)

Xu, Guoxiang; Li, Pengfei; Cao, Qingnan; Hu, Qingxian; Gu, Xiaoyan; Du, Baoshuai

2018-03-01

The present study aims to develop a unified three dimensional numerical model for fiber laser+GMAW hybrid welding, which is used to study the fluid flow phenomena in hybrid welding of aluminum alloy and the influence of laser power on weld pool dynamic behavior. This model takes into account the coupling of gas, liquid and metal phases. Laser heat input is described using a cone heat source model with changing peak power density, its height being determined based on the keyhole size. Arc heat input is modeled as a double ellipsoid heat source. The arc plasma flow and droplet transfer are simulated through the two simplified models. The temperature and velocity fields for different laser powers are calculated. The computed results are in general agreement with the experimental data. Both the peak and average values of fluid flow velocity during hybrid welding are much higher than those of GMAW. At a low level of laser power, both the arc force and droplet impingement force play a relatively large role on fluid flow in the hybrid welding. Keyhole depth always oscillates within a range. With an increase in laser power, the weld pool behavior becomes more complex. An anti-clockwise vortex is generated and the stability of keyhole depth is improved. Besides, the effects of laser power on different driving forces of fluid flow in weld pool are also discussed.

Unsteady flow of fractional Oldroyd-B fluids through rotating annulus

Science.gov (United States)

Tahir, Madeeha; Naeem, Muhammad Nawaz; Javaid, Maria; Younas, Muhammad; Imran, Muhammad; Sadiq, Naeem; Safdar, Rabia

2018-04-01

In this paper exact solutions corresponding to the rotational flow of a fractional Oldroyd-B fluid, in an annulus, are determined by applying integral transforms. The fluid starts moving after t = 0+ when pipes start rotating about their axis. The final solutions are presented in the form of usual Bessel and hypergeometric functions, true for initial and boundary conditions. The limiting cases for the solutions for ordinary Oldroyd-B, fractional Maxwell and Maxwell and Newtonian fluids are obtained. Moreover, the solution is obtained for the fluid when one pipe is rotating and the other one is at rest. At the end of this paper some characteristics of fluid motion, the effect of the physical parameters on the flow and a correlation between different fluid models are discussed. Finally, graphical representations confirm the above affirmation.

Numerical Modeling of Conjugate Heat Transfer in Fluid Network

Science.gov (United States)

Majumdar, Alok

2004-01-01

Fluid network modeling with conjugate heat transfer has many applications in Aerospace engineering. In modeling unsteady flow with heat transfer, it is important to know the variation of wall temperature in time and space to calculate heat transfer between solid to fluid. Since wall temperature is a function of flow, a coupled analysis of temperature of solid and fluid is necessary. In cryogenic applications, modeling of conjugate heat transfer is of great importance to correctly predict boil-off rate in propellant tanks and chill down of transfer lines. In TFAWS 2003, the present author delivered a paper to describe a general-purpose computer program, GFSSP (Generalized Fluid System Simulation Program). GFSSP calculates flow distribution in complex flow circuit for compressible/incompressible, with or without heat transfer or phase change in all real fluids or mixtures. The flow circuit constitutes of fluid nodes and branches. The mass, energy and specie conservation equations are solved at the nodes where as momentum conservation equations are solved at the branches. The proposed paper describes the extension of GFSSP to model conjugate heat transfer. The network also includes solid nodes and conductors in addition to fluid nodes and branches. The energy conservation equations for solid nodes solves to determine the temperatures of the solid nodes simultaneously with all conservation equations governing fluid flow. The numerical scheme accounts for conduction, convection and radiation heat transfer. The paper will also describe the applications of the code to predict chill down of cryogenic transfer line and boil-off rate of cryogenic propellant storage tank.

Analysis of time integration methods for the compressible two-fluid model for pipe flow simulations

NARCIS (Netherlands)

B. Sanderse (Benjamin); I. Eskerud Smith (Ivar); M.H.W. Hendrix (Maurice)

2017-01-01

textabstractIn this paper we analyse different time integration methods for the two-fluid model and propose the BDF2 method as the preferred choice to simulate transient compressible multiphase flow in pipelines. Compared to the prevailing Backward Euler method, the BDF2 scheme has a significantly

Radiotracer techniques for measuring fluid flow and calibrating flow meters

International Nuclear Information System (INIS)

Cooper, E.L.

1987-08-01

Radiotracer techniques can be used to measure accurately both gas and liquid flow rates under operating conditions in a wide range of flow systems. They are ideally suited for calibrating flow meters as well as for measuring unmetered flows in industrial plants. Applications of these techniques range from measuring the flows of fuels and process fluids for energy and mass balance studies to measuring the flows of liquid and airborne effluents for pollution control. This report describes the various radiotracer techniques which can be used to measure fluid flows. The range of application and inherent accuracy of each technique is discussed

Fluid Mechanics An Introduction to the Theory of Fluid Flows

CERN Document Server

Durst, Franz

2008-01-01

Advancements of fluid flow measuring techniques and of computational methods have led to new ways to treat laminar and turbulent flows. These methods are extensively used these days in research and engineering practise. This also requires new ways to teach the subject to students at higher educational institutions in an introductory manner. The book provides the knowledge to students in engineering and natural science needed to enter fluid mechanics applications in various fields. Analytical treatments are provided, based on the Navier-Stokes equations. Introductions are also given into numerical and experimental methods applied to flows. The main benefit the reader will derive from the book is a sound introduction into all aspects of fluid mechanics covering all relevant subfields.

An artificial compressibility CBS method for modelling heat transfer and fluid flow in heterogeneous porous materials

CSIR Research Space (South Africa)

Malan, AG

2011-08-01

Full Text Available to modelling both forced convection as well as heat transfer and fluid flow through heterogeneous saturated porous materials via an edge-based finite volume discretization scheme. A volume-averaged set of local thermal disequilibrium governing equations...

Benchmarking computational fluid dynamics models of lava flow simulation for hazard assessment, forecasting, and risk management

Science.gov (United States)

Dietterich, Hannah; Lev, Einat; Chen, Jiangzhi; Richardson, Jacob A.; Cashman, Katharine V.

2017-01-01

Numerical simulations of lava flow emplacement are valuable for assessing lava flow hazards, forecasting active flows, designing flow mitigation measures, interpreting past eruptions, and understanding the controls on lava flow behavior. Existing lava flow models vary in simplifying assumptions, physics, dimensionality, and the degree to which they have been validated against analytical solutions, experiments, and natural observations. In order to assess existing models and guide the development of new codes, we conduct a benchmarking study of computational fluid dynamics (CFD) models for lava flow emplacement, including VolcFlow, OpenFOAM, FLOW-3D, COMSOL, and MOLASSES. We model viscous, cooling, and solidifying flows over horizontal planes, sloping surfaces, and into topographic obstacles. We compare model results to physical observations made during well-controlled analogue and molten basalt experiments, and to analytical theory when available. Overall, the models accurately simulate viscous flow with some variability in flow thickness where flows intersect obstacles. OpenFOAM, COMSOL, and FLOW-3D can each reproduce experimental measurements of cooling viscous flows, and OpenFOAM and FLOW-3D simulations with temperature-dependent rheology match results from molten basalt experiments. We assess the goodness-of-fit of the simulation results and the computational cost. Our results guide the selection of numerical simulation codes for different applications, including inferring emplacement conditions of past lava flows, modeling the temporal evolution of ongoing flows during eruption, and probabilistic assessment of lava flow hazard prior to eruption. Finally, we outline potential experiments and desired key observational data from future flows that would extend existing benchmarking data sets.

Flow of a non-Newtonian fluid through channels with permeable wall

Energy Technology Data Exchange (ETDEWEB)

Martins-Costa, Maria Laura [Universidade Federal Fluminense, Niteroi, RJ (Brazil). Dept. de Engenharia Mecanica. Lab. de Matematica Teorica e Aplicada]. E-mail: laura@mec.uff.br; Gama, Rogerio M. Saldanha da [Laboratorio Nacional de Computacao Cientifica (LNCC), Petropolis, RJ (Brazil)]. E-mail: rsgama@domain.com.br; Frey, Sergio [Rio Grande do Sul Univ., Porto Alegre, RS (Brazil). Dept. de Engenharia Mecanica. Grupo de Estudos Termicos e Energeticos

2000-07-01

In the present work the momentum transport in two adjacent flow regions is described by means of a continuum theory of mixtures, specially developed to model multiphase phenomena. A generalized Newtonian fluid flows through the permeable wall channel, originating a pure fluid region and a mixture region - where the fluid saturates the porous matrix. The fluid and the porous matrix are treated as continuous constituents of a binary mixture coexisting superposed, each of them occupying simultaneously the whole volume of the mixture. An Ostwald-de Waele behavior is assumed for both the fluid constituent (in the mixture region) and the fluid (in the so-called pure fluid region), while the porous matrix, represented by the solid constituent, is assumed rigid, homogeneous, isotropic and at rest. Compatibility conditions at the interface (pure fluid-mixture) for momentum transfer are proposed and discussed. Assuming no flow across the interface, the velocity should be zero on the solid parts of the boundary and should match the fluid diffusing velocity on the fluid parts of the boundary. Also the shear stress at the pure fluid region is to be balanced by a multiple of the partial shear stress at the mixture region. A minimum principle for the above-described problem, assuming fully developed flow in both regions, is presented, providing an easy and reliable way for carrying out numerical simulations. (author)

Basic study on an energy conversion system using boiling two-phase flows of temperature-sensitive magnetic fluid. Theoretical analysis based on thermal nonequilibrium model and flow visualization using ultrasonic echo

International Nuclear Information System (INIS)

Ishimoto, Jun; Kamiyama, Shinichi; Okubo, Masaaki.

1995-01-01

Effects of magnetic field on the characteristics of boiling two-phase pipe flow of temperature-sensitive magnetic fluid are clarified in detail both theoretically and experimentally. Firstly, governing equations of two-phase magnetic fluid flow based on the thermal nonequilibrium two-fluid model are presented and numerically solved considering evaporation and condensation between gas- and liquid-phases. Next, behaviour of vapor bubbles is visualized with ultrasonic echo in the region of nonuniform magnetic field. This is recorded and processed with an image processor. As a result, the distributions of void fraction in the two-phase flow are obtained. Furthermore, detailed characteristics of the two-phase magnetic fluid flow are investigated using a small test loop of the new energy conversion system. From the numerical and experimental results, it is known that the precise control of the boiling two-phase flow and bubble generation is possible by using the nonuniform magnetic field effectively. These fundamental studies on the characteristics of two-phase magnetic fluid flow will contribute to the development of the new energy conversion system using a gas-liquid boiling two-phase flow of magnetic fluid. (author)

Lattice Boltzmann computation of creeping fluid flow in roll-coating applications

Science.gov (United States)

Rajan, Isac; Kesana, Balashanker; Perumal, D. Arumuga

2018-04-01

Lattice Boltzmann Method (LBM) has advanced as a class of Computational Fluid Dynamics (CFD) methods used to solve complex fluid systems and heat transfer problems. It has ever-increasingly attracted the interest of researchers in computational physics to solve challenging problems of industrial and academic importance. In this current study, LBM is applied to simulate the creeping fluid flow phenomena commonly encountered in manufacturing technologies. In particular, we apply this novel method to simulate the fluid flow phenomena associated with the "meniscus roll coating" application. This prevalent industrial problem encountered in polymer processing and thin film coating applications is modelled as standard lid-driven cavity problem to which creeping flow analysis is applied. This incompressible viscous flow problem is studied in various speed ratios, the ratio of upper to lower lid speed in two different configurations of lid movement - parallel and anti-parallel wall motion. The flow exhibits interesting patterns which will help in design of roll coaters.

The flow and spray characteristics of gelled fluids; Die Stroemungs- und Verspruehungseigenschaften gelfoermiger Fluide

Energy Technology Data Exchange (ETDEWEB)

Madlener, K.

2008-07-01

In the present study gelled fluids are investigated concerning their application as propellants in storable and thrust controllable rocket propulsion systems. The correlations between the non-Newtonian viscosity properties and the flow and spray characteristics are discussed. Based on the proposed viscosity model Herschel-Bulkley-Extended (HBE) the laminar pipe flow is calculated for the investigated propellants. With the introduction of a generalized form of the Reynolds number and the presentation of a possibility to determine the critical values of this number it is possible to calculate the laminar-turbulent transition in a pipe flow. The theoretical results are evaluated with experimental data. The spray characteristics of various gelled fluids are examined using an experimental setup with impinging-jet-injectors. (orig.)

Centrifuge modelling of granular flows

Science.gov (United States)

Cabrera, Miguel Angel; Wu, Wei

2015-04-01

A common characteristic of mass flows like debris flows, rock avalanches and mudflows is that gravity is their main driving force. Gravity defines the intensity and duration of the main interactions between particles and their surrounding media (particle-particle, particle-fluid, fluid-fluid). At the same time, gravity delimits the occurrence of phase separation, inverse segregation, and mass consolidation, among other phenomena. Therefore, in the understanding of the flow physics it is important to account for the scaling of gravity in scaled models. In this research, a centrifuge model is developed to model free surface granular flows down an incline at controlled gravity conditions. Gravity is controlled by the action of an induced inertial acceleration field resulting from the rotation of the model in a geotechnical centrifuge. The characteristics of the induced inertial acceleration field during flow are discussed and validated via experimental data. Flow heights, velocity fields, basal pressure and impact forces are measured for a range of channel inclinations and gravity conditions. Preliminary results enlighten the flow characteristics at variable gravity conditions and open a discussion on the simulation of large scale processes at a laboratory scale. Further analysis on the flow physics brings valuable information for the validation of granular flows rheology.

Simulation of swimming strings immersed in a viscous fluid flow

Science.gov (United States)

Huang, Wei-Xi; Sung, Hyung Jin

2006-11-01

In nature, many phenomena involve interactions between flexible bodies and their surrounding viscous fluid, such as a swimming fish or a flapping flag. The intrinsic dynamics is complicate and not well understood. A flexible string can be regarded as a one-dimensional flag model. Many similarities can be found between the flapping string and swimming fish, although different wake speed results in a drag force for the flapping string and a propulsion force for the swimming fish. In the present study, we propose a mathematical formulation for swimming strings immersed in a viscous fluid flow. Fluid motion is governed by the Navier-Stokes equations and a momentum forcing is added in order to bring the fluid to move at the same velocity with the immersed surface. A flexible inextensible string model is described by another set of equations with an additional momentum forcing which is a result of the fluid viscosity and the pressure difference across the string. The momentum forcing is calculated by a feedback loop. Simulations of several numerical examples are carried out, including a hanging string which starts moving under gravity without ambient fluid, a swinging string immersed in a quiescent viscous fluid, a string swimming within a uniform surrounding flow, and flow over two side-by-side strings. The numerical results agree well with the theoretical analysis and previous experimental observations. Further simulation of a swimming fish is under consideration.

Newtonian heating effects in three-dimensional flow of viscoelastic fluid

International Nuclear Information System (INIS)

Qayyum, A.; Hayat, T.; Alhuthali, M. S.; Malaikah, H. M.

2014-01-01

A mathematical model is constructed to investigate the three-dimensional flow of a non-Newtonian fluid. An incompressible viscoelastic fluid is used in mathematical formulation. The conjugate convective process (in which heat the transfer rate from the bounding surface with a finite capacity is proportional to the local surface temperature) in three-dimensional flow of a differential type of non-Newtonian fluid is analyzed for the first time. Series solutions for the nonlinear differential system are computed. Plots are presented for the description of emerging parameters entering into the problem. It is observed that the conjugate heating phenomenon causes an appreciable increase in the temperature at the stretching wall. (electromagnetism, optics, acoustics, heat transfer, classical mechanics, and fluid dynamics)

Enhanced Geothermal Systems Research and Development: Models of Subsurface Chemical Processes Affecting Fluid Flow

Energy Technology Data Exchange (ETDEWEB)

Moller, Nancy; Weare J. H.

2008-05-29

Successful exploitation of the vast amount of heat stored beneath the earth’s surface in hydrothermal and fluid-limited, low permeability geothermal resources would greatly expand the Nation’s domestic energy inventory and thereby promote a more secure energy supply, a stronger economy and a cleaner environment. However, a major factor limiting the expanded development of current hydrothermal resources as well as the production of enhanced geothermal systems (EGS) is insufficient knowledge about the chemical processes controlling subsurface fluid flow. With funding from past grants from the DOE geothermal program and other agencies, we successfully developed advanced equation of state (EOS) and simulation technologies that accurately describe the chemistry of geothermal reservoirs and energy production processes via their free energies for wide XTP ranges. Using the specific interaction equations of Pitzer, we showed that our TEQUIL chemical models can correctly simulate behavior (e.g., mineral scaling and saturation ratios, gas break out, brine mixing effects, down hole temperatures and fluid chemical composition, spent brine incompatibilities) within the compositional range (Na-K-Ca-Cl-SO4-CO3-H2O-SiO2-CO2(g)) and temperature range (T < 350°C) associated with many current geothermal energy production sites that produce brines with temperatures below the critical point of water. The goal of research carried out under DOE grant DE-FG36-04GO14300 (10/1/2004-12/31/2007) was to expand the compositional range of our Pitzer-based TEQUIL fluid/rock interaction models to include the important aluminum and silica interactions (T < 350°C). Aluminum is the third most abundant element in the earth’s crust; and, as a constituent of aluminosilicate minerals, it is found in two thirds of the minerals in the earth’s crust. The ability to accurately characterize effects of temperature, fluid mixing and interactions between major rock-forming minerals and hydrothermal and

Stress dependent fluid flow in porous rock: experiments and network modelling

Energy Technology Data Exchange (ETDEWEB)

Flornes, Olav

2005-07-01

During the lifetime of a hydrocarbon reservoir, the pore pressure decreases because fluids are drained. Changed pore pressure causes a deformation of the reservoir rock, and the flow channels may be narrowed by the increased weight carried by the rock matrix. Knowledge of how the rocks ability to transport fluids, the permeability, is changed by increased stress can be important for effective reservoir management. In this work, we present experimental results for how permeability changes with applied stress. The materials tested are several different sandstones and one limestone, all having porosities higher than 19 percent. Application of stress is done in a number of different ways. We subject the sample to an isotropic stress, and see how changing this applied stress affects permeability as opposed to changing the pore fluid pressure. This allows for investigating the effective stress law for permeability. Permeability decreased by 10 to 20 percent, when we deformed the materials hydro statically within the elastic regime. For all of our samples, we observed a higher permeability change than predicted by a conventional model for relating porosity and permeability, the Kozeny Carman model. For Red Wildmoor, a sandstone having some clay content, we observed that a change in pore pressure was slightly more important for permeability than a change in the applied stress with the same amount. A sandstone with no clay content, Bad Durckheim, showed the opposite behavior, with applied stress slightly more important than pore pressure. We present a new method for measuring permeability in two directions in the same experiment. We apply different anisotropic stresses, and see if a high stress in one direction causes a difference in permeability changes parallel and perpendicular to maximum stress. We observe that deforming the sample axially, causes a larger decrease in axial permeability than in the radial at low confining pressure. At high confining pressure, the

Time-Dependent Thermally-Driven Interfacial Flows in Multilayered Fluid Structures

Science.gov (United States)

Haj-Hariri, Hossein; Borhan, A.

1996-01-01

A computational study of thermally-driven convection in multilayered fluid structures will be performed to examine the effect of interactions among deformable fluid-fluid interfaces on the structure of time-dependent flow in these systems. Multilayered fluid structures in two models configurations will be considered: the differentially heated rectangular cavity with a free surface, and the encapsulated cylindrical liquid bridge. An extension of a numerical method developed as part of our recent NASA Fluid Physics grant will be used to account for finite deformations of fluid-fluid interfaces.

Modelling of fluid flow in tape casting of thin ceramics: Analytical approaches and numerical investigations

DEFF Research Database (Denmark)

Jabbaribehnam, Mirmasoud; Hattel, Jesper Henri

2016-01-01

Tape casting has been used to produce thin layers of ceramics that can be used as single layers or can be stacked and laminated into multilayered structures. Many startup products such as multilayered inductors, multilayered varistors, piezoelectrics, ceramic fuel cells and lithium ion battery...... components are dependent upon tape casting technology. One of the growing sciences in the processing of ceramics by tape casting is the use of fluid flow analysis to control and enhance the final tapes. The fluid dynamics analysis of the ceramic slurries during tape casting is an efficient mean to elucidate...... the physical parameters crucial to the process. A review of the development of the tape casting process with particular focus on modelling the material flow is presented and in this context the current status is examined and future potential discussed....

Magnetic field effect on blood flow of Casson fluid in axisymmetric cylindrical tube: A fractional model

Energy Technology Data Exchange (ETDEWEB)

Ali, Farhad, E-mail: farhadaliecomaths@yahoo.com [Department of Mathematics, City University of Science and Information Technology, Peshawar 25000 (Pakistan); Sheikh, Nadeem Ahmad [Department of Mathematics, City University of Science and Information Technology, Peshawar 25000 (Pakistan); Khan, Ilyas [Basic Engineering Sciences Department, College of Engineering Majmaah University, Majmaah 11952 (Saudi Arabia); Saqib, Muhammad [Department of Mathematics, City University of Science and Information Technology, Peshawar 25000 (Pakistan)

2017-02-01

The effects of magnetohydrodynamics on the blood flow when blood is represented as a Casson fluid, along with magnetic particles in a horizontal cylinder is studied. The flow is due to an oscillating pressure gradient. The Laplace and finite Hankel transforms are used to obtain the closed form solutions of the fractional partial differential equations. Effects of various parameters on the flow of both blood and magnetic particles are shown graphically. The analysis shows that, the model with fractional order derivatives bring a remarkable changes as compared to the ordinary model. The study highlights that applied magnetic field reduces the velocities of both the blood and magnetic particles.

Magnetic field effect on blood flow of Casson fluid in axisymmetric cylindrical tube: A fractional model

International Nuclear Information System (INIS)

Ali, Farhad; Sheikh, Nadeem Ahmad; Khan, Ilyas; Saqib, Muhammad

2017-01-01

The effects of magnetohydrodynamics on the blood flow when blood is represented as a Casson fluid, along with magnetic particles in a horizontal cylinder is studied. The flow is due to an oscillating pressure gradient. The Laplace and finite Hankel transforms are used to obtain the closed form solutions of the fractional partial differential equations. Effects of various parameters on the flow of both blood and magnetic particles are shown graphically. The analysis shows that, the model with fractional order derivatives bring a remarkable changes as compared to the ordinary model. The study highlights that applied magnetic field reduces the velocities of both the blood and magnetic particles.

A flow meter for ultrasonically measuring the flow velocity of fluids

DEFF Research Database (Denmark)

2015-01-01

The invention regards a flow meter for ultrasonically measuring the flow velocity of fluids comprising a duct having a flow channel with an internal cross section comprising variation configured to generate at least one acoustic resonance within the flow channel for a specific ultrasonic frequency......, and at least two transducers for generating and sensing ultrasonic pulses, configured to transmit ultrasonic pulses at least at said specific ultrasonic frequency into the flow channel such that the ultrasonic pulses propagate through a fluid flowing in the flow channel, wherein the flow meter is configured...

Persistent Homology to describe Solid and Fluid Structures during Multiphase Flow

Science.gov (United States)

Herring, A. L.; Robins, V.; Liu, Z.; Armstrong, R. T.; Sheppard, A.

2017-12-01

The question of how to accurately and effectively characterize essential fluid and solid distributions and structures is a long-standing topic within the field of porous media and fluid transport. For multiphase flow applications, considerable research effort has been made to describe fluid distributions under a range of conditions; including quantification of saturation levels, fluid-fluid pressure differences and interfacial areas, and fluid connectivity. Recent research has effectively used topological metrics to describe pore space and fluid connectivity, with researchers demonstrating links between pore-scale nonwetting phase topology to fluid mobilization and displacement mechanisms, relative permeability, fluid flow regimes, and thermodynamic models of multiphase flow. While topology is clearly a powerful tool to describe fluid distribution, topological metrics by definition provide information only on the connectivity of a phase, not its geometry (shape or size). Physical flow characteristics, e.g. the permeability of a fluid phase within a porous medium, are dependent on the connectivity of the pore space or fluid phase as well as the size of connections. Persistent homology is a technique which provides a direct link between topology and geometry via measurement of topological features and their persistence from the signed Euclidean distance transform of a segmented digital image (Figure 1). We apply persistent homology analysis to measure the occurrence and size of pore-scale topological features in a variety of sandstones, for both the dry state and the nonwetting phase fluid during two-phase fluid flow (drainage and imbibition) experiments, visualized with 3D X-ray microtomography. The results provide key insights into the dominant topological features and length scales of a media which control relevant field-scale engineering properties such as fluid trapping, absolute permeability, and relative permeability.

Fluid flow through replicated microcellular materials in the Darcy-Forchheimer regime

International Nuclear Information System (INIS)

Weber, L.; Ingram, D.; Guardia, S.; Athanasiou-Ioannou, A.; Mortensen, A.

2017-01-01

We extend here a “bottleneck” flow model derived earlier for incompressible fluids flowing under creeping flow conditions [Despois, J. and Mortensen, A: Acta Materialia 53 (2005) 1381] to flow regimes where inertial losses are no longer negligible, causing the governing flow law to deviate from Darcy's law and become the Darcy-Forchheimer law. The proposed law is compared with measurements of the Darcian permeability K_D and of the Forchheimer coefficient C in forced-flow of air through microcellular aluminium made by the replication process. The geometrical features of the cellular medium are varied in terms of volume fraction of porosity (in the range of 0.66–0.86) and the average cell diameter from (108–425 μm). As found previously in measurements with water, the Darcy permeability of the foams for airflow is also reasonably well captured by the model. In the Forchheimer-regime the model gives good quantitative agreement with data if one assumes that the amount of air kinetic energy that is dissipated when passing across each bottleneck linking one pore to its neighbour along the fluid flow path corresponds to the difference, in a stream of constant cross-sectional area, between a uniform fluid velocity profile and the non-uniform profile that is created by the no-slip condition along the window boundary.

Wave Interactions and Fluid Flows

Science.gov (United States)

Craik, Alex D. D.

1988-07-01

This up-to-date and comprehensive account of theory and experiment on wave-interaction phenomena covers fluids both at rest and in their shear flows. It includes, on the one hand, water waves, internal waves, and their evolution, interaction, and associated wave-driven means flow and, on the other hand, phenomena on nonlinear hydrodynamic stability, especially those leading to the onset of turbulence. This study provide a particularly valuable bridge between these two similar, yet different, classes of phenomena. It will be of value to oceanographers, meteorologists, and those working in fluid mechanics, atmospheric and planetary physics, plasma physics, aeronautics, and geophysical and astrophysical fluid dynamics.

Interstitial fluid flow and drug delivery in vascularized tumors: a computational model.

Directory of Open Access Journals (Sweden)

Michael Welter

Full Text Available Interstitial fluid is a solution that bathes and surrounds the human cells and provides them with nutrients and a way of waste removal. It is generally believed that elevated tumor interstitial fluid pressure (IFP is partly responsible for the poor penetration and distribution of therapeutic agents in solid tumors, but the complex interplay of extravasation, permeabilities, vascular heterogeneities and diffusive and convective drug transport remains poorly understood. Here we consider-with the help of a theoretical model-the tumor IFP, interstitial fluid flow (IFF and its impact upon drug delivery within tumor depending on biophysical determinants such as vessel network morphology, permeabilities and diffusive vs. convective transport. We developed a vascular tumor growth model, including vessel co-option, regression, and angiogenesis, that we extend here by the interstitium (represented by a porous medium obeying Darcy's law and sources (vessels and sinks (lymphatics for IFF. With it we compute the spatial variation of the IFP and IFF and determine its correlation with the vascular network morphology and physiological parameters like vessel wall permeability, tissue conductivity, distribution of lymphatics etc. We find that an increased vascular wall conductivity together with a reduction of lymph function leads to increased tumor IFP, but also that the latter does not necessarily imply a decreased extravasation rate: Generally the IF flow rate is positively correlated with the various conductivities in the system. The IFF field is then used to determine the drug distribution after an injection via a convection diffusion reaction equation for intra- and extracellular concentrations with parameters guided by experimental data for the drug Doxorubicin. We observe that the interplay of convective and diffusive drug transport can lead to quite unexpected effects in the presence of a heterogeneous, compartmentalized vasculature. Finally we discuss

Interstitial fluid flow and drug delivery in vascularized tumors: a computational model.

Science.gov (United States)

Welter, Michael; Rieger, Heiko

2013-01-01

Interstitial fluid is a solution that bathes and surrounds the human cells and provides them with nutrients and a way of waste removal. It is generally believed that elevated tumor interstitial fluid pressure (IFP) is partly responsible for the poor penetration and distribution of therapeutic agents in solid tumors, but the complex interplay of extravasation, permeabilities, vascular heterogeneities and diffusive and convective drug transport remains poorly understood. Here we consider-with the help of a theoretical model-the tumor IFP, interstitial fluid flow (IFF) and its impact upon drug delivery within tumor depending on biophysical determinants such as vessel network morphology, permeabilities and diffusive vs. convective transport. We developed a vascular tumor growth model, including vessel co-option, regression, and angiogenesis, that we extend here by the interstitium (represented by a porous medium obeying Darcy's law) and sources (vessels) and sinks (lymphatics) for IFF. With it we compute the spatial variation of the IFP and IFF and determine its correlation with the vascular network morphology and physiological parameters like vessel wall permeability, tissue conductivity, distribution of lymphatics etc. We find that an increased vascular wall conductivity together with a reduction of lymph function leads to increased tumor IFP, but also that the latter does not necessarily imply a decreased extravasation rate: Generally the IF flow rate is positively correlated with the various conductivities in the system. The IFF field is then used to determine the drug distribution after an injection via a convection diffusion reaction equation for intra- and extracellular concentrations with parameters guided by experimental data for the drug Doxorubicin. We observe that the interplay of convective and diffusive drug transport can lead to quite unexpected effects in the presence of a heterogeneous, compartmentalized vasculature. Finally we discuss various

Stochastic Geometric Models with Non-stationary Spatial Correlations in Lagrangian Fluid Flows

Science.gov (United States)

Gay-Balmaz, François; Holm, Darryl D.

2018-01-01

Inspired by spatiotemporal observations from satellites of the trajectories of objects drifting near the surface of the ocean in the National Oceanic and Atmospheric Administration's "Global Drifter Program", this paper develops data-driven stochastic models of geophysical fluid dynamics (GFD) with non-stationary spatial correlations representing the dynamical behaviour of oceanic currents. Three models are considered. Model 1 from Holm (Proc R Soc A 471:20140963, 2015) is reviewed, in which the spatial correlations are time independent. Two new models, called Model 2 and Model 3, introduce two different symmetry breaking mechanisms by which the spatial correlations may be advected by the flow. These models are derived using reduction by symmetry of stochastic variational principles, leading to stochastic Hamiltonian systems, whose momentum maps, conservation laws and Lie-Poisson bracket structures are used in developing the new stochastic Hamiltonian models of GFD.

Stochastic Geometric Models with Non-stationary Spatial Correlations in Lagrangian Fluid Flows

Science.gov (United States)

Gay-Balmaz, François; Holm, Darryl D.

2018-06-01

Inspired by spatiotemporal observations from satellites of the trajectories of objects drifting near the surface of the ocean in the National Oceanic and Atmospheric Administration's "Global Drifter Program", this paper develops data-driven stochastic models of geophysical fluid dynamics (GFD) with non-stationary spatial correlations representing the dynamical behaviour of oceanic currents. Three models are considered. Model 1 from Holm (Proc R Soc A 471:20140963, 2015) is reviewed, in which the spatial correlations are time independent. Two new models, called Model 2 and Model 3, introduce two different symmetry breaking mechanisms by which the spatial correlations may be advected by the flow. These models are derived using reduction by symmetry of stochastic variational principles, leading to stochastic Hamiltonian systems, whose momentum maps, conservation laws and Lie-Poisson bracket structures are used in developing the new stochastic Hamiltonian models of GFD.

Experimental and Numerical Modeling of Fluid Flow Processes in Continuous Casting: Results from the LIMMCAST-Project

Science.gov (United States)

Timmel, K.; Kratzsch, C.; Asad, A.; Schurmann, D.; Schwarze, R.; Eckert, S.

2017-07-01

The present paper reports about numerical simulations and model experiments concerned with the fluid flow in the continuous casting process of steel. This work was carried out in the LIMMCAST project in the framework of the Helmholtz alliance LIMTECH. A brief description of the LIMMCAST facilities used for the experimental modeling at HZDR is given here. Ultrasonic and inductive techniques and the X-ray radioscopy were employed for flow measurements or visualizations of two-phase flow regimes occurring in the submerged entry nozzle and the mold. Corresponding numerical simulations were performed at TUBAF taking into account the dimensions and properties of the model experiments. Numerical models were successfully validated using the experimental data base. The reasonable and in many cases excellent agreement of numerical with experimental data allows to extrapolate the models to real casting configurations. Exemplary results will be presented here showing the effect of electromagnetic brakes or electromagnetic stirrers on the flow in the mold or illustrating the properties of two-phase flows resulting from an Ar injection through the stopper rod.

STAFAN, Fluid Flow, Mechanical Stress in Fractured Rock of Nuclear Waste Repository

International Nuclear Information System (INIS)

Huyakorn, P.; Golis, M.J.

1989-01-01

1 - Description of program or function: STAFAN (Stress And Flow Analysis) is a two-dimensional, finite-element code designed to model fluid flow and the interaction of fluid pressure and mechanical stresses in a fractured rock surrounding a nuclear waste repository. STAFAN considers flow behavior of a deformable fractured system with fracture-porous matrix interactions, the coupling effects of fluid pressure and mechanical stresses in a medium containing discrete joints, and the inelastic response of the individual joints of the rock mass subject to the combined fluid pressure and mechanical loading. 2 - Restrictions on the complexity of the problem: STAFAN does not presently contain thermal coupling, and it is unable to simulate inelastic deformation of the rock mass and variably saturated or two-phase flow in the fractured porous medium system

Bifurcation and stability of an improved time-delayed fluid flow model in internet congestion control

International Nuclear Information System (INIS)

Yu-Liang, Liu; Jie, Zhu; Xiao-Shu, Luo

2009-01-01

Based on the fluid flow time-delayed model proposed by Misra et al in internet congestion control, one modified time-delayed model is presented, where the influence of the communication delay on the router queue length is investigated in detail. The main advantage of the new model is that its stability domain is larger even without an extra controller. By linear stability analysis and numerical simulation, the effectiveness and feasibility of the novel model in internet congestion control are verified

Bifurcation and stability of an improved time-delayed fluid flow model in internet congestion control

Science.gov (United States)

Liu, Yu-Liang; Zhu, Jie; Luo, Xiao-Shu

2009-09-01

Based on the fluid flow time-delayed model proposed by Misra et al in internet congestion control, one modified time-delayed model is presented, where the influence of the communication delay on the router queue length is investigated in detail. The main advantage of the new model is that its stability domain is larger even without an extra controller. By linear stability analysis and numerical simulation, the effectiveness and feasibility of the novel model in internet congestion control are verified.

Assessment of fluid-to-fluid modelling of critical heat flux in horizontal 37-element bundle flows

International Nuclear Information System (INIS)

Yang, S.K.

2006-01-01

Fluid-to-fluid modelling laws of critical heat flux (CHF) available in the literature were reviewed. The applicability of the fluid-to-fluid modelling laws was assessed using available data ranging from low to high mass fluxes in horizontal 37-element bundles simulating a CANDU fuel string. Correlations consisting of dimensionless similarity groups were derived using modelling fluid data (Freon-12) to predict water CHF data in horizontal 37-element bundles with uniform and non-uniform axial-heat flux distribution (AFD). The results showed that at mass fluxes higher than ∼4,000 kg/m 2 s (water equivalent value), the vertical fluid-to-fluid modelling laws of Ahmad (1973) and Katto (1979) predict water CHF in horizontal 37-element bundles with non-uniform AFD with average errors of 1.4% and 3.0% and RMS errors of 5.9% and 6.1%, respectively. The Francois and Berthoud (2003) fluid-to-fluid modelling law predicts CHF in non-uniformly heated 37-element bundles in the horizontal orientation with an average error of 0.6% and an RMS error of 10.4% over the available range of 2,000 to 6,200 kg/m 2 s. (author)

Fluid Flow Behaviour under Different Gases and Flow Rate during Gas Metal Arc Welding

OpenAIRE

Jaison Peter

2013-01-01

Gas metal arc welding (GMAW) is a highly efficient and fast process for fabricating high quality weld. High quality welds are fabricated by proper selection of consumable includes gas and filler metals. The optimum flow rate of gas will ensure the proper quality of weld. In this project, a fluid flow behavior of different flow rate is modeled and the change quality will be studied.

Computational fluid dynamic modelling of cavitation

Science.gov (United States)

Deshpande, Manish; Feng, Jinzhang; Merkle, Charles L.

1993-01-01

Models in sheet cavitation in cryogenic fluids are developed for use in Euler and Navier-Stokes codes. The models are based upon earlier potential-flow models but enable the cavity inception point, length, and shape to be determined as part of the computation. In the present paper, numerical solutions are compared with experimental measurements for both pressure distribution and cavity length. Comparisons between models are also presented. The CFD model provides a relatively simple modification to an existing code to enable cavitation performance predictions to be included. The analysis also has the added ability of incorporating thermodynamic effects of cryogenic fluids into the analysis. Extensions of the current two-dimensional steady state analysis to three-dimensions and/or time-dependent flows are, in principle, straightforward although geometrical issues become more complicated. Linearized models, however offer promise of providing effective cavitation modeling in three-dimensions. This analysis presents good potential for improved understanding of many phenomena associated with cavity flows.

Principal component analysis in an experimental cold flow model of a fluid catalytic cracking unit by gammametry

International Nuclear Information System (INIS)

Araujo, Janeo Severino C. de; Dantas, Carlos Costa; Santos, Valdemir A. dos; Souza, Jose Edson G. de; Luna-Finkler, Christine L.

2009-01-01

The fluid dynamic behavior of riser of a cold flow model of a Fluid Catalytic Cracking Unit (FCCU) was investigated. The experimental data were obtained by the nuclear technique of gamma transmission. A gamma source was placed diametrically opposite to a detector in any straight section of the riser. The gas-solid flow through riser was monitored with a source of Americium-241 what allowed obtaining information of the axial solid concentration without flow disturbance and also identifying the dependence of this concentration profile with several independent variables. The MatLab R and Statistica R software were used. Statistica tool employed was the Principal Components Analysis (PCA), that consisted of the job of the data organization, through two-dimensional head offices to allow extract relevant information about the importance of the independent variables on axial solid concentration in a cold flow riser. The variables investigated were mass flow rate of solid, mass flow rate of gas, pressure in the riser base and the relative height in the riser. The first two components reached about 98 % of accumulated percentage of explained variance. (author)

Fluid flow near the surface of earth's outer core

Science.gov (United States)

Bloxham, Jeremy; Jackson, Andrew

1991-01-01

This review examines the recent attempts at extracting information on the pattern of fluid flow near the surface of the outer core from the geomagnetic secular variation. Maps of the fluid flow at the core surface are important as they may provide some insight into the process of the geodynamo and may place useful constraints on geodynamo models. In contrast to the case of mantle convection, only very small lateral variations in core density are necessary to drive the flow; these density variations are, by several orders of magnitude, too small to be imaged seismically; therefore, the geomagnetic secular variation is utilized to infer the flow. As substantial differences exist between maps developed by different researchers, the possible underlying reasons for these differences are examined with particular attention given to the inherent problems of nonuniqueness.

Dynamic Modeling Strategy for Flow Regime Transition in Gas-Liquid Two-Phase Flows

Directory of Open Access Journals (Sweden)

Xia Wang

2012-12-01

Full Text Available In modeling gas-liquid two-phase flows, the concept of flow regimes has been widely used to characterize the global interfacial structure of the flows. Nearly all constitutive relations that provide closures to the interfacial transfers in two-phase flow models, such as the two-fluid model, are flow regime dependent. Current nuclear reactor safety analysis codes, such as RELAP5, classify flow regimes using flow regime maps or transition criteria that were developed for steady-state, fully-developed flows. As two-phase flows are dynamic in nature, it is important to model the flow regime transitions dynamically to more accurately predict the two-phase flows. The present work aims to develop a dynamic modeling strategy to determine flow regimes in gas-liquid two-phase flows through introduction of interfacial area transport equations (IATEs within the framework of a two-fluid model. The IATE is a transport equation that models the interfacial area concentration by considering the creation of the interfacial area, fluid particle (bubble or liquid droplet disintegration, boiling and evaporation, and the destruction of the interfacial area, fluid particle coalescence and condensation. For flow regimes beyond bubbly flows, a two-group IATE has been proposed, in which bubbles are divided into two groups based on their size and shapes, namely group-1 and group-2 bubbles. A preliminary approach to dynamically identify the flow regimes is discussed, in which discriminators are based on the predicted information, such as the void fraction and interfacial area concentration. The flow regime predicted with this method shows good agreement with the experimental observations.

Flow regime classification in air-magnetic fluid two-phase flow.

Science.gov (United States)

Kuwahara, T; De Vuyst, F; Yamaguchi, H

2008-05-21

A new experimental/numerical technique of classification of flow regimes (flow patterns) in air-magnetic fluid two-phase flow is proposed in the present paper. The proposed technique utilizes the electromagnetic induction to obtain time-series signals of the electromotive force, allowing us to make a non-contact measurement. Firstly, an experiment is carried out to obtain the time-series signals in a vertical upward air-magnetic fluid two-phase flow. The signals obtained are first treated using two kinds of wavelet transforms. The data sets treated are then used as input vectors for an artificial neural network (ANN) with supervised training. In the present study, flow regimes are classified into bubbly, slug, churn and annular flows, which are generally the main flow regimes. To validate the flow regimes, a visualization experiment is also performed with a glycerin solution that has roughly the same physical properties, i.e., kinetic viscosity and surface tension, as a magnetic fluid used in the present study. The flow regimes from the visualization are used as targets in an ANN and also used in the estimation of the accuracy of the present method. As a result, ANNs using radial basis functions are shown to be the most appropriate for the present classification of flow regimes, leading to small classification errors.

A non-traditional fluid problem: transition between theoretical models from Stokes’ to turbulent flow

Science.gov (United States)

Salomone, Horacio D.; Olivieri, Néstor A.; Véliz, Maximiliano E.; Raviola, Lisandro A.

2018-05-01

In the context of fluid mechanics courses, it is customary to consider the problem of a sphere falling under the action of gravity inside a viscous fluid. Under suitable assumptions, this phenomenon can be modelled using Stokes’ law and is routinely reproduced in teaching laboratories to determine terminal velocities and fluid viscosities. In many cases, however, the measured physical quantities show important deviations with respect to the predictions deduced from the simple Stokes’ model, and the causes of these apparent ‘anomalies’ (for example, whether the flow is laminar or turbulent) are seldom discussed in the classroom. On the other hand, there are various variable-mass problems that students tackle during elementary mechanics courses and which are discussed in many textbooks. In this work, we combine both kinds of problems and analyse—both theoretically and experimentally—the evolution of a system composed of a sphere pulled by a chain of variable length inside a tube filled with water. We investigate the effects of different forces acting on the system such as weight, buoyancy, viscous friction and drag force. By means of a sequence of mathematical models of increasing complexity, we obtain a progressive fit that accounts for the experimental data. The contrast between the various models exposes the strengths and weaknessess of each one. The proposed experience can be useful for integrating concepts of elementary mechanics and fluids, and is suitable as laboratory practice, stressing the importance of the experimental validation of theoretical models and showing the model-building processes in a didactic framework.

Computational fluid dynamics (CFD) simulation of hot air flow ...

African Journals Online (AJOL)

Computational Fluid Dynamics simulation of air flow distribution, air velocity and pressure field pattern as it will affect moisture transient in a cabinet tray dryer is performed using SolidWorks Flow Simulation (SWFS) 2014 SP 4.0 program. The model used for the drying process in this experiment was designed with Solid ...

Fast intraslab fluid-flow events linked to pulses of high pore fluid pressure at the subducted plate interface

Science.gov (United States)

Taetz, Stephan; John, Timm; Bröcker, Michael; Spandler, Carl; Stracke, Andreas

2018-01-01

A better understanding of the subduction zone fluid cycle and its chemical-mechanical feedback requires in-depth knowledge about how fluids flow within and out of descending slabs. Relicts of fluid-flow systems in exhumed rocks of fossil subduction zones allow for identification of the general relationships between dehydration reactions, fluid pathway formation, the dimensions and timescales of distinct fluid flow events; all of which are required for quantitative models for fluid-induced subduction zone processes. Two types of garnet-quartz-phengite veins can be distinguished in an eclogite-facies mélange block from the Pouébo Eclogite Mélange, New Caledonia. These veins record synmetamorphic internal fluid release by mineral breakdown reactions (type I veins), and infiltration of an external fluid (type II veins) with the associated formation of a reaction selvage. The dehydration and fluid migration documented by the type I veins likely occurred on a timescale of 105-106 years, based on average subduction rates and metamorphic conditions required for mineral dehydration and fluid flow. The timeframe of fluid-rock interaction between the external fluid and the wall-rock of the type II veins is quantified using a continuous bulk-rock Li-diffusion profile perpendicular to a vein and its metasomatic selvage. Differences in Li concentration between the internal and external fluid reservoirs resulted in a distinct diffusion profile (decreasing Li concentration and increasing δ7 Li) as the reaction front propagated into the host rock. Li-chronometric constraints indicate that the timescales of fluid-rock interaction associated with type II vein formation are on the order of 1 to 4 months (0.150-0.08+0.14 years). The short-lived, pulse-like character of this process is consistent with the notion that fluid flow caused by oceanic crust dehydration at the blueschist-to-eclogite transition contributes to or even dominates episodic pore fluid pressure increases at the

A Novel Plasma-Based Fluid for Particle Image Velocimetry (PIV): In-Vitro Feasibility Study of Flow Diverter Effects in Aneurysm Model.

Science.gov (United States)

Clauser, Johanna; Knieps, Marius S; Büsen, Martin; Ding, Andreas; Schmitz-Rode, Thomas; Steinseifer, Ulrich; Arens, Jutta; Cattaneo, Giorgio

2018-02-27

Particle image velocimetry (PIV) is a commonly used method for in vitro investigation of fluid dynamics in biomedical devices, such as flow diverters for intracranial aneurysm treatment. Since it is limited to transparent blood substituting fluids like water-glycerol mixture, the influence of coagulation and platelet aggregation is neglected. We aimed at the development and the application of a modified platelet rich plasma as a new PIV fluid with blood-like rheological and coagulation properties. In standardized intracranial aneurysm silicone models, the effect of this new PIV plasma on the fluid dynamics before and after flow diverter implantation was evaluated and compared with water-glycerol measurements. The flow diverting effect was strongly dependent on the used fluid, with considerably lower velocities achieved using PIV plasma, despite the same starting viscosity of both fluids. Moreover, triggering coagulation of PIV plasma allowed for intra-aneurysmal clot formation. We presented the first in vitro PIV investigation using a non-Newtonian, clottable PIV plasma, demonstrating a mismatch to a standard PIV fluid and allowing for thrombus formation.

Considerations of a nonhomogeneous fluid in the deep groundwater flow system at Hanford

International Nuclear Information System (INIS)

Nelson, R.W.

1988-11-01

This report presents such a general theory capable of describing the flow on nonhomogeneous fluids in porous media, theory that is a composite from several disciplines including groundwater hydrology, soil physics, civil engineering, petroleum reservoir engineering, mechanics, and mathematical physics. The report discussed the conceptual basis for considering the flow of nonhomogeneous fluids. From this conceptual basis emphasis shifts to providing complete definitions and then appropriately describing those definitions in mathematical terms. Throughout the report, the necessary assumptions are stated in detail because the limitations of any theory are best assessed through careful scrutiny of the assumptions. From the mathematical definitions with appropriate functional dependence the results and constraints needed are derived to provide the general theory necessary to describe the flow of nonhomogeneous fluids in porous media. Particular attention is given to comparing the general theory with the classical theory of flow for a homogeneous fluid. Such comparison provides significant insight to the effects of variable fluid properties on subsurface flow systems. The comparisons also indicate the importance of carefully formulating subsurface flow models within the more general theoretical framework describing the flow of nonhomogeneous fluids in porous media. 29 refs.; 6 figs.; 1 tab

Mathematical well-posedness of a two-fluid equations for bubbly two-phase flows

International Nuclear Information System (INIS)

Okawa, Tomio; Kataoka, Isao

2000-01-01

It is widely known that two-fluid equations used in most engineering applications do not satisfy the necessary condition for being mathematical well-posed as initial-value problems. In the case of stratified two-phase flows, several researchers have revealed that differential models satisfying the necessary condition are to be derived if the pressure difference between the phases is related to the spatial gradient of the void fraction through the effects of gravity or surface tension. While, in the case of dispersed two-phase flows, no physically reasonable method to derive mathematically well-posed two-fluid model has been proposed. In the present study, particularly focusing on the effect of interfacial pressure terms, we derived the mathematically closed form of the volume-averaged two-fluid model for bubbly two-phase flows. As a result of characteristic analyses, it was shown that the proposed two-fluid equations satisfy the necessary condition of mathematical well-posedness if the void fraction is sufficiently small. (author)

Numerical Modeling of Fluid Flow in the Tape Casting Process

DEFF Research Database (Denmark)

Jabbari, Masoud; Hattel, Jesper Henri

2011-01-01

The flow behavior of the fluid in the tape casting process is analyzed. A simple geometry is assumed for running the numerical calculations in ANSYS Fluent and the main parameters are expressed in non-dimensional form. The effect of different values for substrate velocity and pressure force...

Boundary control of fluid flow through porous media

DEFF Research Database (Denmark)

Hasan, Agus; Foss, Bjarne; Sagatun, Svein Ivar

2010-01-01

The flow of fluids through porous media can be described by the Boussinesq’s equation with mixed boundary conditions; a Neumann’s boundary condition and a nonlinear boundary condition. The nonlinear boundary condition provides a means to control the fluid flow through porous media. In this paper,......, some stabilizing controllers are constructed for various cases using Lyapunov design.......The flow of fluids through porous media can be described by the Boussinesq’s equation with mixed boundary conditions; a Neumann’s boundary condition and a nonlinear boundary condition. The nonlinear boundary condition provides a means to control the fluid flow through porous media. In this paper...

Development of a numerical model for fluid-structure interaction analysis of flow through and around an aquaculture net cage

DEFF Research Database (Denmark)

Chen, Hao; Christensen, Erik Damgaard

2017-01-01

In the present work, we developed a numerical model for fluid-structure interaction analysis of flow through and around an aquaculture net cage. The numerical model is based on the coupling between the porous media model and the lumped mass structural model. A novel interface was implemented...

Modeling fluid transport in 2d paper networks

Science.gov (United States)

Tirapu Azpiroz, Jaione; Fereira Silva, Ademir; Esteves Ferreira, Matheus; Lopez Candela, William Fernando; Bryant, Peter William; Ohta, Ricardo Luis; Engel, Michael; Steiner, Mathias Bernhard

2018-02-01

Paper-based microfluidic devices offer great potential as a low-cost platform to perform chemical and biochemical tests. Commercially available formats such as dipsticks and lateral-flow test devices are widely popular as they are easy to handle and produce fast and unambiguous results. While these simple devices lack precise control over the flow to enable integration of complex functionality for multi-step processes or the ability to multiplex several tests, intense research in this area is rapidly expanding the possibilities. Modeling and simulation is increasingly more instrumental in gaining insight into the underlying physics driving the processes inside the channels, however simulation of flow in paper-based microfluidic devices has barely been explored to aid in the optimum design and prototyping of these devices for precise control of the flow. In this paper, we implement a multiphase fluid flow model through porous media for the simulation of paper imbibition of an incompressible, Newtonian fluid such as when water, urine or serum is employed. The formulation incorporates mass and momentum conservation equations under Stokes flow conditions and results in two coupled Darcy's law equations for the pressures and saturations of the wetting and non-wetting phases, further simplified to the Richard's equation for the saturation of the wetting fluid, which is then solved using a Finite Element solver. The model tracks the wetting fluid front as it displaces the non-wetting fluid by computing the time-dependent saturation of the wetting fluid. We apply this to the study of liquid transport in two-dimensional paper networks and validate against experimental data concerning the wetting dynamics of paper layouts of varying geometries.

The correct Q1D electrodynamical part of a model for the fluid flow in a faraday segmented magnetohydrodynamic generator channel

International Nuclear Information System (INIS)

Bajovic, V.S.

1995-01-01

A Faraday ideally segmented (the absence of the Hall electric current I x =0) MHD generator channel, with a weakly ionized plasma as a working fluid, is considered. The magnetic field is applied along the z axis and the working fluid flows along the x axis. The stationary state of the flow and the stationary electric current are considered. The new quasi-one-dimensional (Q1D) electrodynamical part of a model is developed. The main assumption (besides j z = 0, E z = 0) taken in the whole working fluid flow, is j x = 0. It means that the bending of the electric current pattern in the working fluid, due to the nonmassive electrodes and the presence of the magnetic field, is not explicitly considered. In the frame of the assumption taken, the legitimacy regarding the equation of conservation of charges: div rvec j = 0 suggests that the straight electric current pattern in the model should be narrow. In other words, it suggests replacing of a rather complicated electric current pattern in a working fluid by an artificially straight and narrow one in the modeling, promising that it would still be possible to describe the global channel parameters - its electrical output, the influence of the shape and the size of the channel on the generator performance etc

Fluid flow and convective transport of solutes within the intervertebral disc.

Science.gov (United States)

Ferguson, Stephen J; Ito, Keita; Nolte, Lutz P

2004-02-01

Previous experimental and analytical studies of solute transport in the intervertebral disc have demonstrated that for small molecules diffusive transport alone fulfils the nutritional needs of disc cells. It has been often suggested that fluid flow into and within the disc may enhance the transport of larger molecules. The goal of the study was to predict the influence of load-induced interstitial fluid flow on mass transport in the intervertebral disc. An iterative procedure was used to predict the convective transport of physiologically relevant molecules within the disc. An axisymmetric, poroelastic finite-element structural model of the disc was developed. The diurnal loading was divided into discrete time steps. At each time step, the fluid flow within the disc due to compression or swelling was calculated. A sequentially coupled diffusion/convection model was then employed to calculate solute transport, with a constant concentration of solute being provided at the vascularised endplates and outer annulus. Loading was simulated for a complete diurnal cycle, and the relative convective and diffusive transport was compared for solutes with molecular weights ranging from 400 Da to 40 kDa. Consistent with previous studies, fluid flow did not enhance the transport of low-weight solutes. During swelling, interstitial fluid flow increased the unidirectional penetration of large solutes by approximately 100%. Due to the bi-directional temporal nature of disc loading, however, the net effect of convective transport over a full diurnal cycle was more limited (30% increase). Further study is required to determine the significance of large solutes and the timing of their delivery for disc physiology.

Performance of intact and partially degraded concrete barriers in limiting fluid flow

International Nuclear Information System (INIS)

Walton, J.C.; Seitz, R.R.

1991-07-01

Concrete barriers will play a critical role in the long-term isolation of low-level radioactive wastes. Over time the barriers will degrade, and in many cases, the fundamental processes controlling performance of the barriers will be different for intact and degraded conditions. This document examines factors controlling fluid flow through intact and degraded concrete disposal facilities. Simplified models are presented fro predicting build up of fluid above a vault; fluid flow through and around intact vaults, through flaws in coatings/liners applied to a vault, and through cracks in a concrete vault; and the influence of different backfill materials around the outside of the vault. Example calculations are presented to illustrate the parameters and processes that influence fluid flow. 46 refs., 49 figs., 2 tabs

Numerical solution of chemically reactive non-Newtonian fluid flow: Dual stratification

Science.gov (United States)

Rehman, Khalil Ur; Malik, M. Y.; Khan, Abid Ali; Zehra, Iffat; Zahri, Mostafa; Tahir, M.

2017-12-01

We have found that only a few attempts are available in the literature relatively to the tangent hyperbolic fluid flow induced by stretching cylindrical surfaces. In particular, temperature and concentration stratification effects have not been investigated until now with respect to the tangent hyperbolic fluid model. Therefore, we have considered the tangent hyperbolic fluid flow induced by an acutely inclined cylindrical surface in the presence of both temperature and concentration stratification effects. To be more specific, the fluid flow is attained with the no slip condition, which implies that the bulk motion of the fluid particles is the same as the stretching velocity of a cylindrical surface. Additionally, the flow field situation is manifested with heat generation, mixed convection and chemical reaction effects. The flow partial differential equations give a complete description of the present problem. Therefore, to trace out the solution, a set of suitable transformations is introduced to convert these equations into ordinary differential equations. In addition, a self-coded computational algorithm is executed to inspect the numerical solution of these reduced equations. The effect logs of the involved parameters are provided graphically. Furthermore, the variations of the physical quantities are examined and given with the aid of tables. It is observed that the fluid temperature is a decreasing function of the thermal stratification parameter and a similar trend is noticed for the concentration via the solutal stratification parameter.

Fluid-elastic force measurements acting on a tube bundle in two-phase cross flow

International Nuclear Information System (INIS)

Inada, Fumio; Kawamura, Koji; Yasuo, Akira

1996-01-01

Fluid-elastic force acting on a square tube bundle of P/D = 1.47 in air-water two-phase cross flow was measured to investigate the characteristics and to clarify whether the fluid elastic vibration characteristics could be expressed using two-phase mixture characteristics. Measured fluid elastic forces were separated into fluid-elastic force coefficients such as added mass, added stiffness, and added damping coefficient. The added damping coefficient was separated into a two-phase damping and a flow-dependent component as in previous research (Carlucci, 1981 and 1983; Pettigrew, 1994). These coefficients were nondimensionalized with two-phase mixture characteristics such as void fraction, mixture density and mixture velocity, which were obtained using the drift-flux model with consideration given to the model. The result was compared with the result obtained with the homogeneous model. It was found that fluid-elastic force coefficients could be expressed with two-phase flow mixture characteristics very well in the experimental result, and that better result can be derived using the slip model as compared to the homogeneous model. Added two-phase flow, which could be expressed as a function of void fraction, where two-phase damping was nondimensionalized with the relative velocity between the gas and liquid phases used as a reference velocity. Using these, the added stiffness coefficient and flow-dependent component of damping could be expressed very well as a function of nondimensional mixture velocity

Computational fluid dynamics simulations of single-phase flow in a filter-press flow reactor having a stack of three cells

International Nuclear Information System (INIS)

Sandoval, Miguel A.; Fuentes, Rosalba; Walsh, Frank C.; Nava, José L.; Ponce de León, Carlos

2016-01-01

Highlights: • Computational fluid dynamic simulations in a filter-press stack of three cells. • The fluid velocity was different in each cell due to local turbulence. • The upper cell link pipe of the filter press cell acts as a fluid mixer. • The fluid behaviour tends towards a continuous mixing flow pattern. • Close agreement between simulations and experimental data was achieved. - Abstract: Computational fluid dynamics (CFD) simulations were carried out for single-phase flow in a pre-pilot filter press flow reactor with a stack of three cells. Velocity profiles and streamlines were obtained by solving the Reynolds-Averaged Navier-Stokes (RANS) equations with a standard k − ε turbulence model. The flow behaviour shows the appearance of jet flow at the entrance to each cell. At lengths from 12 to 15 cm along the cells channels, a plug flow pattern is developed at all mean linear flow rates studied here, 1.2 ≤ u ≤ 2.1 cm s −1 . The magnitude of the velocity profiles in each cell was different, due to the turbulence generated by the change of flow direction in the last fluid manifold. Residence time distribution (RTD) simulations indicated that the fluid behaviour tends towards a continuous mixing flow pattern, owing to flow at the output of each cell across the upper cell link pipe, which acts as a mixer. Close agreement between simulations and experimental RTD was obtained.

Simulation of biological flow and transport in complex geometries using embedded boundary/volume-of-fluid methods

International Nuclear Information System (INIS)

Trebotich, David

2007-01-01

We have developed a simulation capability to model multiscale flow and transport in complex biological systems based on algorithms and software infrastructure developed under the SciDAC APDEC CET. The foundation of this work is a new hybrid fluid-particle method for modeling polymer fluids in irregular microscale geometries that enables long-time simulation of validation experiments. Both continuum viscoelastic and discrete particle representations have been used to model the constitutive behavior of polymer fluids. Complex flow environment geometries are represented on Cartesian grids using an implicit function. Direct simulation of flow in the irregular geometry is then possible using embedded boundary/volume-of-fluid methods without loss of geometric detail. This capability has been used to simulate biological flows in a variety of application geometries including biomedical microdevices, anatomical structures and porous media

Interfacing a General Purpose Fluid Network Flow Program with the SINDA/G Thermal Analysis Program

Science.gov (United States)

Schallhorn, Paul; Popok, Daniel

1999-01-01

A general purpose, one dimensional fluid flow code is currently being interfaced with the thermal analysis program Systems Improved Numerical Differencing Analyzer/Gaski (SINDA/G). The flow code, Generalized Fluid System Simulation Program (GFSSP), is capable of analyzing steady state and transient flow in a complex network. The flow code is capable of modeling several physical phenomena including compressibility effects, phase changes, body forces (such as gravity and centrifugal) and mixture thermodynamics for multiple species. The addition of GFSSP to SINDA/G provides a significant improvement in convective heat transfer modeling for SINDA/G. The interface development is conducted in multiple phases. This paper describes the first phase of the interface which allows for steady and quasi-steady (unsteady solid, steady fluid) conjugate heat transfer modeling.

Topological fluid mechanics of Axisymmetric Flow

DEFF Research Database (Denmark)

Brøns, Morten

1998-01-01

Topological fluid mechanics in the sense of the present paper is the study and classification of flow patterns close to a critical point. Here we discuss the topology of steady viscous incompressible axisymmetric flows in the vicinity of the axis. Following previous studies the velocity field v...... to the authors knowledge has not been used systematically to high orders in topological fluid mechanics. We compare the general results with experimental and computational results on the Vogel-Ronneberg flow. We show that the topology changes observed when recirculating bubbles on the vortex axis are created...

Fluid Mechanics of Taylor Bubbles and Slug Flows in Vertical Channels

International Nuclear Information System (INIS)

Anglart, Henryk; Podowski, Michael Z.

2002-01-01

Fluid mechanics of Taylor bubbles and slug flows is investigated in vertical, circular channels using detailed, three-dimensional computational fluid dynamics simulations. The Volume of Fluid model with the interface-sharpening algorithm, implemented in the commercial CFX4 code, is used to predict the shape and velocity of Taylor bubbles moving along a vertical channel. Several cases are investigated, including both a single Taylor bubble and a train of bubbles rising in water. It is shown that the potential flow solution underpredicts the water film thickness around Taylor bubbles. Furthermore, the computer simulations that are performed reveal the importance of properly modeling the three-dimensional nature of phenomena governing the motion of Taylor bubbles. Based on the present results, a new formula for the evaluation of bubble shape is derived. Both the shape of Taylor bubbles and the bubble rise velocity predicted by the proposed model agree well with experimental observations. Furthermore, the present model shows good promise in predicting the coalescence of Taylor bubbles

Application of stereolithographic custom models for studying the impact of biofilms and mineral precipitation on fluid flow.

Science.gov (United States)

Stoner, D L; Watson, S M; Stedtfeld, R D; Meakin, P; Griffel, L K; Tyler, T L; Pegram, L M; Barnes, J M; Deason, V A

2005-12-01

Here we introduce the use of transparent experimental models fabricated by stereolithography for studying the impacts of biomass accumulation, minerals precipitation, and physical configuration of flow paths on liquid flow in fracture apertures. The internal configuration of the models ranged in complexity from simple geometric shapes to those that incorporate replicated surfaces of natural fractures and computationally derived fracture surfaces. High-resolution digital time-lapse imaging was employed to qualitatively observe the migration of colloidal and soluble dyes through the flow models. In this study, a Sphingomonas sp. and Sporosarcina (Bacillus) pasteurii influenced the fluid dynamics by physically altering flow paths. Microbial colonization and calcite deposition enhanced the stagnant regions adjacent to solid boundaries. Microbial growth and calcite precipitation occurred to a greater extent in areas behind the fabricated obstacles and less in high-velocity orifices.

Analysis of fluid flow around a beating artificial cilium

Directory of Open Access Journals (Sweden)

Mojca Vilfan

2012-02-01

Full Text Available Biological cilia are found on surfaces of some microorganisms and on surfaces of many eukaryotic cells where they interact with the surrounding fluid. The periodic beating of the cilia is asymmetric, resulting in directed swimming of unicellular organisms or in generation of a fluid flow above a ciliated surface in multicellular ones. Following the biological example, externally driven artificial cilia have recently been successfully implemented as micropumps and mixers. However, biomimetic systems are useful not only in microfluidic applications, but can also serve as model systems for the study of fundamental hydrodynamic phenomena in biological samples. To gain insight into the basic principles governing propulsion and fluid pumping on a micron level, we investigated hydrodynamics around one beating artificial cilium. The cilium was composed of superparamagnetic particles and driven along a tilted cone by a varying external magnetic field. Nonmagnetic tracer particles were used for monitoring the fluid flow generated by the cilium. The average flow velocity in the pumping direction was obtained as a function of different parameters, such as the rotation frequency, the asymmetry of the beat pattern, and the cilium length. We also calculated the velocity field around the beating cilium by using the analytical far-field expansion. The measured average flow velocity and the theoretical prediction show an excellent agreement.

CFD simulations of flow erosion and flow-induced deformation of needle valve: Effects of operation, structure and fluid parameters

Energy Technology Data Exchange (ETDEWEB)

Zhu, Hongjun, E-mail: ticky863@126.com [State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation, Southwest Petroleum University, Chengdu 610500, Sichuan (China); State Key Laboratory of Hydraulics and Mountain River Engineering, Sichuan University, Chengdu 610065, Sichuan (China); Pan, Qian; Zhang, Wenli; Feng, Guang; Li, Xue [State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation, Southwest Petroleum University, Chengdu 610500, Sichuan (China)

2014-07-01

Highlights: • A combined FSI–CFD and DPM computational method is used to investigate flow erosion and deformation of needle valve. • The numerical model is validated with the comparison of measured and predicted erosion rate. • Effects of operation, structure and fluid parameters on flow erosion and flow-induced deformation are discussed. • Particle diameter has the most significant effect on flow erosion. • Inlet rate has the most obvious effect on flow-induced deformation. - Abstract: A three-dimensional fluid–structure interaction (FSI) computational model coupling with a combined continuum and discrete model has been used to predict the flow erosion rate and flow-induced deformation of needle valve. Comparisons with measured data demonstrate good agreement with the predictions of erosion rate. The flow field distribution of gas-particle flow and the erosion rate and deformation of valve core are captured under different operating and structural conditions with different fluid parameters. The effects of inlet velocity, valve opening and inlet valve channel size, particle concentration, particle diameter and particle phase components are discussed in detail. The results indicate that valve tip has the most severe erosion and deformation, and flow field, erosion rate and deformation of valve are all sensitive to inlet condition changes, structural changes and fluid properties changes. The effect of particle diameter on erosion is the most significant, while the influence of inlet rate on deformation is the greatest one.

Numerical solution of viscous and viscoelastic fluids flow through the branching channel by finite volume scheme

Science.gov (United States)

Keslerová, Radka; Trdlička, David

2015-09-01

This work deals with the numerical modelling of steady flows of incompressible viscous and viscoelastic fluids through the three dimensional channel with T-junction. The fundamental system of equations is the system of generalized Navier-Stokes equations for incompressible fluids. This system is based on the system of balance laws of mass and momentum for incompressible fluids. Two different mathematical models for the stress tensor are used for simulation of Newtonian and Oldroyd-B fluids flow. Numerical solution of the described models is based on cetral finite volume method using explicit Runge-Kutta time integration.

AFDM: An Advanced Fluid-Dynamics Model

International Nuclear Information System (INIS)

Wilhelm, D.

1990-09-01

This volume describes the Advanced Fluid-Dynamics Model (AFDM) for topologies, flow regimes, and interfacial areas. The objective of these models is to provide values for the interfacial areas between all components existing in a computational cell. The interfacial areas are then used to evaluate the mass, energy, and momentum transfer between the components. A new approach has been undertaken in the development of a model to convect the interfacial areas of the discontinuous velocity fields in the three-velocity-field environment of AFDM. These interfacial areas are called convectible surface areas. The continuous and discontinuous components are chosen using volume fraction and levitation criteria. This establishes so-called topologies for which the convectible surface areas can be determined. These areas are functions of space and time. Solid particulates that are limited to being discontinuous within the bulk fluid are assumed to have a constant size. The convectible surface areas are subdivided to model contacts between two discontinuous components or discontinuous components and the structure. The models have been written for the flow inside of large pools. Therefore, the structure is tracked only as a boundary to the fluid volume without having a direct influence on velocity or volume fraction distribution by means of flow regimes or boundary layer models. 17 refs., 7 tabs., 18 figs

Fluid flow dynamics in MAS systems

Science.gov (United States)

Wilhelm, Dirk; Purea, Armin; Engelke, Frank

2015-08-01

The turbine system and the radial bearing of a high performance magic angle spinning (MAS) probe with 1.3 mm-rotor diameter has been analyzed for spinning rates up to 67 kHz. We focused mainly on the fluid flow properties of the MAS system. Therefore, computational fluid dynamics (CFD) simulations and fluid measurements of the turbine and the radial bearings have been performed. CFD simulation and measurement results of the 1.3 mm-MAS rotor system show relatively low efficiency (about 25%) compared to standard turbo machines outside the realm of MAS. However, in particular, MAS turbines are mainly optimized for speed and stability instead of efficiency. We have compared MAS systems for rotor diameter of 1.3-7 mm converted to dimensionless values with classical turbomachinery systems showing that the operation parameters (rotor diameter, inlet mass flow, spinning rate) are in the favorable range. This dimensionless analysis also supports radial turbines for low speed MAS probes and diagonal turbines for high speed MAS probes. Consequently, a change from Pelton type MAS turbines to diagonal turbines might be worth considering for high speed applications. CFD simulations of the radial bearings have been compared with basic theoretical values proposing considerably smaller frictional loss values. The discrepancies might be due to the simple linear flow profile employed for the theoretical model. Frictional losses generated inside the radial bearings result in undesired heat-up of the rotor. The rotor surface temperature distribution computed by CFD simulations show a large temperature gradient over the rotor.

Numerical investigation of fluid mud motion using a three-dimensional hydrodynamic and two-dimensional fluid mud coupling model

Science.gov (United States)

Yang, Xiaochen; Zhang, Qinghe; Hao, Linnan

2015-03-01

A water-fluid mud coupling model is developed based on the unstructured grid finite volume coastal ocean model (FVCOM) to investigate the fluid mud motion. The hydrodynamics and sediment transport of the overlying water column are solved using the original three-dimensional ocean model. A horizontal two-dimensional fluid mud model is integrated into the FVCOM model to simulate the underlying fluid mud flow. The fluid mud interacts with the water column through the sediment flux, current, and shear stress. The friction factor between the fluid mud and the bed, which is traditionally determined empirically, is derived with the assumption that the vertical distribution of shear stress below the yield surface of fluid mud is identical to that of uniform laminar flow of Newtonian fluid in the open channel. The model is validated by experimental data and reasonable agreement is found. Compared with numerical cases with fixed friction factors, the results simulated with the derived friction factor exhibit the best agreement with the experiment, which demonstrates the necessity of the derivation of the friction factor.

Lattice gas automaton scheme with stochastic particle movement for a rotated fluid flow

International Nuclear Information System (INIS)

Ishiguro, Misako

2002-01-01

Lattice gas automaton (LGA) models developed so far are just for Cartesian geometries, and no direct approach to rotated fluid flows is found. In this paper, LGA method is applied to model a two-dimensional rotated flow. Several problems specific to the rotated flow are to be solved: hexagonal lattice geometry to effectively identify the neighbors, boundary condition for irregular walls, multi-speed scheme to represent angular-oriented fluid velocity υ θ ≅γω, shape of macroscopic domain for statistics, formula to obtain macroscopic quantities such as density and mean fluid velocities, application method of Fermi-Dirac function to the initial particle arrangement. For this purpose, FHP-I type hexagonal lattice model is revised and a new LGA model with stochastic particle movement is proposed. The results of the trial calculation are shown. It is also investigated whether or not the underlying microscopic Boolean equations newly introduced leads to Navier-Stokes equation. (author)

Interface model coupling in fluid dynamics: application to two-phase flows

International Nuclear Information System (INIS)

Galie, Th.

2009-03-01

This thesis is devoted to the study of interface model coupling problems in space between different models of compressible flows. We consider one-dimensional problems where the interface is sharp, fixed and separating two regions of space corresponding to the two coupled models. Our goal is to define a coupling condition at the interface and to solve numerically the coupling problem with this condition. After a state of art on the interface model coupling of hyperbolic systems of conservation laws, we propose a new coupling condition by adding in the equations of the coupled problem a measure source term at the interface. We first suppose a given constant weight associated to this source term. Two Riemann solvers are developed and one of them is based on a relaxation approach preserving equilibrium solutions of the coupled problem. This relaxation method is then used in an optimization problem, defined by several motivations at the interface, which permits to calculate a time dynamical weight. In a second part, we develop an approached Riemann solver for a two-phase two-pressure model in the particular case of a two-phase isentropic flow. Such a model contains non conservative terms that we write under the form of measure source terms. The previous relaxation method is thus extended to the case of the two-phase two-pressure model with an a priori estimation of the non conservative term contributions. The method allows us to solve, in the next and last chapter, the coupling problem of a two-fluid two-pressure model with a drift-flux model thanks to the father model approach. (authors)

The Challenge of Fluid Flow

Indian Academy of Sciences (India)

makes fluid flows so rich, so complex - some times so highly ordered ..... to other frequencies, which again can grow in amplitude before they also eventually decay again. On the ..... think of it as a slice of flow issuing towards this sheet of paper.

Nonlinear transport processes and fluid dynamics: Cylindrical Couette flow of Lennard-Jones fluids

International Nuclear Information System (INIS)

Khayat, R.E.; Eu, B.C.

1988-01-01

In this paper we report on calculations of flow profiles for cylindrical Couette flow of a Lennard-Jones fluid. The flow is subjected to a temperature gradient and thermoviscous effects are taken into consideration. We apply the generalized fluid dynamic equations which are provided by the modified moment method for the Boltzmann equation reported previously. The results of calculations are in good agreement with the Monte Carlo direct simulation method by K. Nanbu [Phys. Fluids 27, 2632 (1984)] for most of Knudsen numbers for which the simulation data are available

Second order numerical method of two-fluid model of air-water flow

International Nuclear Information System (INIS)

Tiselj, I.; Petelin, S.

1995-01-01

Model considered in this paper is six-equation two-fluid model used in computer code RELAP5. Air-water equations were taken in a code named PDE to avoid additional problems caused by condensation or vaporization. Terms with space derivatives were added in virtual mass term in momentum equations to ensure the hyperbolicity of the equations. Numerical method in PDE code is based on approximate Riemann solvers. Equations are solved on non-staggered grid with explicit time advancement and with upwind discretization of the convective terms in characteristic form of the equations. Flux limiters are used to find suitable combinations of the first (upwind) and the second order (Lax-Wendroff) discretization s which ensure second order accuracy on smooth solutions and damp oscillations around the discontinuities. Because of the small time steps required and because of its non-dissipative nature the scheme is suitable for the prediction of the fast transients: pressure waves, shock and rarefaction waves, water hammer or critical flow. Some preliminary results are presented for a shock tube problem and for Water Faucet problem - problems usually used as benchmarks for two-fluid computer codes. (author)

A solution algorithm for fluid-particle flows across all flow regimes

Science.gov (United States)

Kong, Bo; Fox, Rodney O.

2017-09-01

Many fluid-particle flows occurring in nature and in technological applications exhibit large variations in the local particle volume fraction. For example, in circulating fluidized beds there are regions where the particles are close-packed as well as very dilute regions where particle-particle collisions are rare. Thus, in order to simulate such fluid-particle systems, it is necessary to design a flow solver that can accurately treat all flow regimes occurring simultaneously in the same flow domain. In this work, a solution algorithm is proposed for this purpose. The algorithm is based on splitting the free-transport flux solver dynamically and locally in the flow. In close-packed to moderately dense regions, a hydrodynamic solver is employed, while in dilute to very dilute regions a kinetic-based finite-volume solver is used in conjunction with quadrature-based moment methods. To illustrate the accuracy and robustness of the proposed solution algorithm, it is implemented in OpenFOAM for particle velocity moments up to second order, and applied to simulate gravity-driven, gas-particle flows exhibiting cluster-induced turbulence. By varying the average particle volume fraction in the flow domain, it is demonstrated that the flow solver can handle seamlessly all flow regimes present in fluid-particle flows.

Simulation of Two-Fluid Flows by the Least-Squares Finite Element Method Using a Continuum Surface Tension Model

Science.gov (United States)

Wu, Jie; Yu, Sheng-Tao; Jiang, Bo-nan

1996-01-01

In this paper a numerical procedure for simulating two-fluid flows is presented. This procedure is based on the Volume of Fluid (VOF) method proposed by Hirt and Nichols and the continuum surface force (CSF) model developed by Brackbill, et al. In the VOF method fluids of different properties are identified through the use of a continuous field variable (color function). The color function assigns a unique constant (color) to each fluid. The interfaces between different fluids are distinct due to sharp gradients of the color function. The evolution of the interfaces is captured by solving the convective equation of the color function. The CSF model is used as a means to treat surface tension effect at the interfaces. Here a modified version of the CSF model, proposed by Jacqmin, is used to calculate the tension force. In the modified version, the force term is obtained by calculating the divergence of a stress tensor defined by the gradient of the color function. In its analytical form, this stress formulation is equivalent to the original CSF model. Numerically, however, the use of the stress formulation has some advantages over the original CSF model, as it bypasses the difficulty in approximating the curvatures of the interfaces. The least-squares finite element method (LSFEM) is used to discretize the governing equation systems. The LSFEM has proven to be effective in solving incompressible Navier-Stokes equations and pure convection equations, making it an ideal candidate for the present applications. The LSFEM handles all the equations in a unified manner without any additional special treatment such as upwinding or artificial dissipation. Various bench mark tests have been carried out for both two dimensional planar and axisymmetric flows, including a dam breaking, oscillating and stationary bubbles and a conical liquid sheet in a pressure swirl atomizer.

Flow Characterization of Vapor Phase of Geothermal Fluid in Pipe Using Isotope 85Kr and Residence Time Distribution Modeling

Directory of Open Access Journals (Sweden)

S. Sugiharto

2014-08-01

Full Text Available Measurement of vapor flow in geothermal pipe faces great challenges due to fast fluids flow in high-temperature and high-pressure environment. In present study the flow rate measurement has been performed to characterization the geothermal vapor flow in a pipe. The experiment was carried out in a pipe which is connected to a geothermal production well, KMJ-14. The pipe has a 10” outside diameter and contains dry vapor at a pressure of 8 kg/cm2 and a temperature of 170 oC. Krypton-85 gas isotope (85Kr has been injected into the pipe. Three collimated radiation detectors positioned respectively at 127, 177 and 227m from injection point were used to obtain experimental data which represent radiotracer residence time distribution (RTD in the pipe. The last detector at the position of 227 m did not respond, which might be due to problems in cable connections. Flow properties calculated using mean residence time (MRT shows that the flow rate of the vapor in pipe is 10.98 m/s, much faster than fluid flow commonly found in various industrial process plants. Best fitting evaluated using dedicated software developed by IAEA expert obtained the Péclet number Pe as 223. This means that the flow of vapor of geothermal fluids in pipe is plug flow in character. The molecular diffusion coefficient is 0.45 m2/s, calculated from the axial dispersion model.

Computational Fluid Dynamics Modeling of a Lithium/Thionyl Chloride Battery with Electrolyte Flow

Energy Technology Data Exchange (ETDEWEB)

Gu, W.B.; Jungst, Rudolph G.; Nagasubramanian, Ganesan; Wang, C.Y.; Weidner, John.

1999-06-11

A two-dimensional model is developed to simulate discharge of a lithium/thionyl chloride primary battery. The model accounts for not only transport of species and charge, but also the electrode porosity variations and the electrolyte flow induced by the volume reduction caused by electrochemical reactions. Numerical simulations are performed using a finite volume method of computational fluid dynamics. The predicted discharge curves for various temperatures are compared to the experimental data with excellent agreement. Moreover, the simulation results. in conjunction with computer visualization and animation techniques, confirm that cell utilization in the temperature and current range of interest is limited by pore plugging or clogging of the front side of the cathode as a result of LiCl precipitation. The detailed two-dimensional flow simulation also shows that the electrolyte is replenished from the cell header predominantly through the separator into the front of the cathode during most parts of the discharge, especially for higher cell temperatures.

Magnetic resonance imaging of cerebrospinal fluid flow in pediatrics

Energy Technology Data Exchange (ETDEWEB)

Heroux, R. [Children' s Hospital of Eastern Ontario, Magnetic Resonance Imaging Dept., Ottawa, Ontario (Canada)

2000-06-30

Magnetic Resonance Imaging of flowing protons in cerebrospinal fluid is useful for demonstrating areas of obstruction or stenosis of the ventricular system causing hydrocephalus. This is used in pediatric patients to assess the circulation of the cerebrospinal fluid. This article discusses two studies. In the first, the cerebrospinal fluid flow study helped the neurosurgeon assess the patency after a third ventriculocisternostomy. The second study evaluated the cerebrospinal fluid flowing through the foramen magnum in a patient with cerebellar tonsilar descent (Chiari malformation) and a syringomyelia. Different techniques to evaluate the flow studies are also discussed. (author)

Free vibration of a single-walled carbon nanotube containing a fluid flow using the Timoshenko beam model

International Nuclear Information System (INIS)

Chang, W.-J.; Lee, H.-L.

2009-01-01

The flexural vibration of the fluid-conveying single-walled carbon nanotube (SWCNT) is derived by the Timoshenko beam model, including rotary inertia and transverse shear deformation. The effects of the flow velocity and the aspect ratio of length to diameter on the vibration frequency and mode shape of the SWCNT are analyzed. Results show that the effects of rotary inertia and transverse shear deformation result in a reduction of the vibration frequencies, especially for higher modes of vibration and short nanotubes. The frequency is also compared with the previous study based on Euler beam model. In addition, if the ratio of length to diameter increased to 60, the influence of the shear deformation and rotary inertia on the mode shape and the resonant frequencies can be neglected. However, the influence is very obvious when the ratio decreased to 20. As the flow velocity of the fluid increases in the vicinity of 2π, the SWCNT reveals the divergence instability. It regains stability when the flow velocity reaches about 9. As the velocity increases further, the SWCNT undergoes a coupled-mode flutter and results in a larger amplitude

Fluid Flow Programming in Paper-Derived Silica-Polymer Hybrids.

Science.gov (United States)

Dubois, Christelle; Herzog, Nicole; Rüttiger, Christian; Geißler, Andreas; Grange, Eléonor; Kunz, Ulrike; Kleebe, Hans-Joachim; Biesalski, Markus; Meckel, Tobias; Gutmann, Torsten; Gallei, Markus; Andrieu-Brunsen, Annette

2017-01-10

In paper-based devices, capillary fluid flow is based on length-scale selective functional control within a hierarchical porous system. The fluid flow can be tuned by altering the paper preparation process, which controls parameters such as the paper grammage. Interestingly, the fiber morphology and nanoporosity are often neglected. In this work, porous voids are incorporated into paper by the combination of dense or mesoporous ceramic silica coatings with hierarchically porous cotton linter paper. Varying the silica coating leads to significant changes in the fluid flow characteristics, up to the complete water exclusion without any further fiber surface hydrophobization, providing new approaches to control fluid flow. Additionally, functionalization with redox-responsive polymers leads to reversible, dynamic gating of fluid flow in these hybrid paper materials, demonstrating the potential of length scale specific, dynamic, and external transport control.

Nonlinear state-space modelling of the kinematics of an oscillating circular cylinder in a fluid flow

Science.gov (United States)

Decuyper, J.; De Troyer, T.; Runacres, M. C.; Tiels, K.; Schoukens, J.

2018-01-01

The flow-induced vibration of bluff bodies is an important problem of many marine, civil, or mechanical engineers. In the design phase of such structures, it is vital to obtain good predictions of the fluid forces acting on the structure. Current methods rely on computational fluid dynamic simulations (CFD), with a too high computational cost to be effectively used in the design phase or for control applications. Alternative methods use heuristic mathematical models of the fluid forces, but these lack the accuracy (they often assume the system to be linear) or flexibility to be useful over a wide operating range. In this work we show that it is possible to build an accurate, flexible and low-computational-cost mathematical model using nonlinear system identification techniques. This model is data driven: it is trained over a user-defined region of interest using data obtained from experiments or simulations, or both. Here we use a Van der Pol oscillator as well as CFD simulations of an oscillating circular cylinder to generate the training data. Then a discrete-time polynomial nonlinear state-space model is fit to the data. This model relates the oscillation of the cylinder to the force that the fluid exerts on the cylinder. The model is finally validated over a wide range of oscillation frequencies and amplitudes, both inside and outside the so-called lock-in region. We show that forces simulated by the model are in good agreement with the data obtained from CFD.

Stretched flow of Oldroyd-B fluid with Cattaneo-Christov heat flux

Directory of Open Access Journals (Sweden)

T. Hayat

Full Text Available The objective of present attempt is to analyse the flow and heat transfer in the flow of an Oldroyd-B fluid over a non-linear stretching sheet having variable thickness. Characteristics of heat transfer are analyzed with temperature dependent thermal conductivity and heat source/sink. Cattaneo-Christov heat flux model is considered rather than Fourier’s law of heat conduction in the present flow analysis. Thermal conductivity varies with temperature. Resulting partial differential equations through laws of conservation of mass, linear momentum and energy are converted into ordinary differential equations by suitable transformations. Convergent series solutions for the velocity and temperature distributions are developed and discussed. Keywords: Oldroyd-B fluid, Variable sheet thickness, Cattaneo-Christov heat flux model, Heat source/sink, Temperature dependent thermal conductivity

A novel investigation of a micropolar fluid characterized by nonlinear constitutive diffusion model in boundary layer flow and heat transfer

Science.gov (United States)

Sui, Jize; Zhao, Peng; Cheng, Zhengdong; Zheng, Liancun; Zhang, Xinxin

2017-02-01

The rheological and heat-conduction constitutive models of micropolar fluids (MFs), which are important non-Newtonian fluids, have been, until now, characterized by simple linear expressions, and as a consequence, the non-Newtonian performance of such fluids could not be effectively captured. Here, we establish the novel nonlinear constitutive models of a micropolar fluid and apply them to boundary layer flow and heat transfer problems. The nonlinear power law function of angular velocity is represented in the new models by employing generalized "n-diffusion theory," which has successfully described the characteristics of non-Newtonian fluids, such as shear-thinning and shear-thickening fluids. These novel models may offer a new approach to the theoretical understanding of shear-thinning behavior and anomalous heat transfer caused by the collective micro-rotation effects in a MF with shear flow according to recent experiments. The nonlinear similarity equations with a power law form are derived and the approximate analytical solutions are obtained by the homotopy analysis method, which is in good agreement with the numerical solutions. The results indicate that non-Newtonian behaviors involving a MF depend substantially on the power exponent n and the modified material parameter K 0 introduced by us. Furthermore, the relations of the engineering interest parameters, including local boundary layer thickness, local skin friction, and Nusselt number are found to be fitted by a quadratic polynomial to n with high precision, which enables the extraction of the rapid predictions from a complex nonlinear boundary-layer transport system.

Stochastic porous media modeling and high-resolution schemes for numerical simulation of subsurface immiscible fluid flow transport

Science.gov (United States)

Brantson, Eric Thompson; Ju, Binshan; Wu, Dan; Gyan, Patricia Semwaah

2018-04-01

This paper proposes stochastic petroleum porous media modeling for immiscible fluid flow simulation using Dykstra-Parson coefficient (V DP) and autocorrelation lengths to generate 2D stochastic permeability values which were also used to generate porosity fields through a linear interpolation technique based on Carman-Kozeny equation. The proposed method of permeability field generation in this study was compared to turning bands method (TBM) and uniform sampling randomization method (USRM). On the other hand, many studies have also reported that, upstream mobility weighting schemes, commonly used in conventional numerical reservoir simulators do not accurately capture immiscible displacement shocks and discontinuities through stochastically generated porous media. This can be attributed to high level of numerical smearing in first-order schemes, oftentimes misinterpreted as subsurface geological features. Therefore, this work employs high-resolution schemes of SUPERBEE flux limiter, weighted essentially non-oscillatory scheme (WENO), and monotone upstream-centered schemes for conservation laws (MUSCL) to accurately capture immiscible fluid flow transport in stochastic porous media. The high-order schemes results match well with Buckley Leverett (BL) analytical solution without any non-oscillatory solutions. The governing fluid flow equations were solved numerically using simultaneous solution (SS) technique, sequential solution (SEQ) technique and iterative implicit pressure and explicit saturation (IMPES) technique which produce acceptable numerical stability and convergence rate. A comparative and numerical examples study of flow transport through the proposed method, TBM and USRM permeability fields revealed detailed subsurface instabilities with their corresponding ultimate recovery factors. Also, the impact of autocorrelation lengths on immiscible fluid flow transport were analyzed and quantified. A finite number of lines used in the TBM resulted into visual

Optimization of micropillar sequences for fluid flow sculpting

Energy Technology Data Exchange (ETDEWEB)

Stoecklein, Daniel; Ganapathysubramanian, Baskar [Department of Mechanical Engineering, Iowa State University, Ames, Iowa 50011 (United States); Wu, Chueh-Yu; Kim, Donghyuk; Di Carlo, Dino [Department of Bioengineering, University of California at Los Angeles, Los Angeles, California 90095 (United States)

2016-01-15

Inertial fluid flow deformation around pillars in a microchannel is a new method for controlling fluid flow. Sequences of pillars have been shown to produce a rich phase space with a wide variety of flow transformations. Previous work has successfully demonstrated manual design of pillar sequences to achieve desired transformations of the flow cross section, with experimental validation. However, such a method is not ideal for seeking out complex sculpted shapes as the search space quickly becomes too large for efficient manual discovery. We explore fast, automated optimization methods to solve this problem. We formulate the inertial flow physics in microchannels with different micropillar configurations as a set of state transition matrix operations. These state transition matrices are constructed from experimentally validated streamtraces for a fixed channel length per pillar. This facilitates modeling the effect of a sequence of micropillars as nested matrix-matrix products, which have very efficient numerical implementations. With this new forward model, arbitrary micropillar sequences can be rapidly simulated with various inlet configurations, allowing optimization routines quick access to a large search space. We integrate this framework with the genetic algorithm and showcase its applicability by designing micropillar sequences for various useful transformations. We computationally discover micropillar sequences for complex transformations that are substantially shorter than manually designed sequences. We also determine sequences for novel transformations that were difficult to manually design. Finally, we experimentally validate these computational designs by fabricating devices and comparing predictions with the results from confocal microscopy.

Theory and computer simulation of structure, transport, and flow of fluid in micropores

International Nuclear Information System (INIS)

Davis, H.T.; Bitsanis, I.; Vanderlick, T.K.; Tirrell, M.V.

1987-01-01

An overview is given of recent progress made in our laboratory on this topic. The density profiles of fluid in micropores are found by solving numerically an approximate Yvon-Born-Green equation. A related local average density model (LADM) allows prediction of transport and flow in inhomogeneous fluids from density profiles. A rigorous extension of the Enskog theory of transport is also outlined. Simple results of this general approach for the tracer diffusion and Couette flow between planar micropore walls are presented. Equilibrium and flow (molecular dynamics) simulations are compared with the theoretical predictions. Simulated density profiles of the micropore fluid exhibit substantial fluid layering. The number and sharpness of fluid layers depend sensitively on the pore width. The solvation force and the pore average density and diffusivity are oscillating functions of the pore width. The theoretical predictions for these quantities agree qualitatively with the simulation results. The flow simulations indicate that the flow does not affect the fluid structure and diffusivity even at extremely high shear rates (10/sup 10/s/sup -1/). The fluid structure induces large deviations of the shear stress and the effective viscosity from the bulk fluid values. The flow velocity profiles are correlated with the density profiles and differ from those of a bulk fluid. The LADM and extended Enskog theory predictions for the velocity profiles and the pore average diffusivity agree very well with each other and with the simulation results. The LADM predictions for the shear stress and the effective viscosity agrees fairly well with the simulation results

Advances in Computational Fluid-Structure Interaction and Flow Simulation Conference

CERN Document Server

Takizawa, Kenji

2016-01-01

This contributed volume celebrates the work of Tayfun E. Tezduyar on the occasion of his 60th birthday. The articles it contains were born out of the Advances in Computational Fluid-Structure Interaction and Flow Simulation (AFSI 2014) conference, also dedicated to Prof. Tezduyar and held at Waseda University in Tokyo, Japan on March 19-21, 2014. The contributing authors represent a group of international experts in the field who discuss recent trends and new directions in computational fluid dynamics (CFD) and fluid-structure interaction (FSI). Organized into seven distinct parts arranged by thematic topics, the papers included cover basic methods and applications of CFD, flows with moving boundaries and interfaces, phase-field modeling, computer science and high-performance computing (HPC) aspects of flow simulation, mathematical methods, biomedical applications, and FSI. Researchers, practitioners, and advanced graduate students working on CFD, FSI, and related topics will find this collection to be a defi...

Numerical Investigation on Fluid Flow in a 90-Degree Curved Pipe with Large Curvature Ratio

Directory of Open Access Journals (Sweden)

Yan Wang

2015-01-01

Full Text Available In order to understand the mechanism of fluid flows in curved pipes, a large number of theoretical and experimental researches have been performed. As a critical parameter of curved pipe, the curvature ratio δ has received much attention, but most of the values of δ are very small (δ<0.1 or relatively small (δ≤0.5. As a preliminary study and simulation this research studied the fluid flow in a 90-degree curved pipe of large curvature ratio. The Detached Eddy Simulation (DES turbulence model was employed to investigate the fluid flows at the Reynolds number range from 5000 to 20000. After validation of the numerical strategy, the pressure and velocity distribution, pressure drop, fluid flow, and secondary flow along the curved pipe were illustrated. The results show that the fluid flow in a curved pipe with large curvature ratio seems to be unlike that in a curved pipe with small curvature ratio. Large curvature ratio makes the internal flow more complicated; thus, the flow patterns, the separation region, and the oscillatory flow are different.

Probability model of solid to liquid-like transition of a fluid suspension after a shear flow onset

Czech Academy of Sciences Publication Activity Database

Nouar, C.; Říha, Pavel

2008-01-01

Roč. 34, č. 5 (2008), s. 477-483 ISSN 0301-9322 R&D Projects: GA AV ČR IAA200600803 Institutional research plan: CEZ:AV0Z20600510 Keywords : laminar suspension flow * liquid-liquid interface * probability model Subject RIV: BK - Fluid Dynamics Impact factor: 1.497, year: 2008

Hydrodynamic cavitation in Stokes flow of anisotropic fluids

Science.gov (United States)

Stieger, Tillmann; Agha, Hakam; Schoen, Martin; Mazza, Marco G.; Sengupta, Anupam

2017-05-01

Cavitation, the nucleation of vapour in liquids, is ubiquitous in fluid dynamics, and is often implicated in a myriad of industrial and biomedical applications. Although extensively studied in isotropic liquids, corresponding investigations in anisotropic liquids are largely lacking. Here, by combining liquid crystal microfluidic experiments, nonequilibrium molecular dynamics simulations and theoretical arguments, we report flow-induced cavitation in an anisotropic fluid. The cavitation domain nucleates due to sudden pressure drop upon flow past a cylindrical obstacle within a microchannel. For an anisotropic fluid, the inception and growth of the cavitation domain ensued in the Stokes regime, while no cavitation was observed in isotropic liquids flowing under similar hydrodynamic parameters. Using simulations we identify a critical value of the Reynolds number for cavitation inception that scales inversely with the order parameter of the fluid. Strikingly, the critical Reynolds number for anisotropic fluids can be 50% lower than that of isotropic fluids.

Biomimetic structures for fluid drag reduction in laminar and turbulent flows

International Nuclear Information System (INIS)

Jung, Yong Chae; Bhushan, Bharat

2010-01-01

Biomimetics allows one to mimic nature to develop materials and devices of commercial interest for engineers. Drag reduction in fluid flow is one of the examples found in nature. In this study, nano, micro, and hierarchical structures found in lotus plant surfaces, as well as shark skin replica and a rib patterned surface to simulate shark skin structure were fabricated. Drag reduction efficiency studies on the surfaces were systematically carried out using water flow. An experimental flow channel was used to measure the pressure drop in laminar and turbulent flows, and the trends were explained in terms of the measured and predicted values by using fluid dynamics models. The slip length for various surfaces in laminar flow was also investigated based on the measured pressure drop. For comparison, the pressure drop for various surfaces was also measured using air flow.

Mathematical simulation of fluid flow and analysis of flow pattern in the flow path of low-head Kaplan turbine

Directory of Open Access Journals (Sweden)

A. V. Rusanov

2016-12-01

Full Text Available The results of numerical investigation of spatial flow of viscous incompressible fluid in flow part of Kaplan turbine PL20 Kremenchug HPP at optimum setting angle of runner blade φb = 15° and at maximum setting angle φb = 35° are shown. The flow simulation has been carried out on basis of numerical integration of the Reynolds equations with an additional term containing artificial compressibility. The differential two-parameter model of Menter (SST has been applied to take into account turbulent effects. Numerical integration of the equations is carried out using an implicit quasi-monotone Godunov type scheme of second - order accuracy in space and time. The calculations have been conducted with the help of the software system IPMFlow. The analysis of fluid flow in the flow part elements is shown and the values of hydraulic losses and local cavitation coefficient have been obtained. Comparison of calculated and experimental results has been carried out.

Multiphase flow models for hydraulic fracturing technology

Science.gov (United States)

Osiptsov, Andrei A.

2017-10-01

The technology of hydraulic fracturing of a hydrocarbon-bearing formation is based on pumping a fluid with particles into a well to create fractures in porous medium. After the end of pumping, the fractures filled with closely packed proppant particles create highly conductive channels for hydrocarbon flow from far-field reservoir to the well to surface. The design of the hydraulic fracturing treatment is carried out with a simulator. Those simulators are based on mathematical models, which need to be accurate and close to physical reality. The entire process of fracture placement and flowback/cleanup can be conventionally split into the following four stages: (i) quasi-steady state effectively single-phase suspension flow down the wellbore, (ii) particle transport in an open vertical fracture, (iii) displacement of fracturing fluid by hydrocarbons from the closed fracture filled with a random close pack of proppant particles, and, finally, (iv) highly transient gas-liquid flow in a well during cleanup. The stage (i) is relatively well described by the existing hydralics models, while the models for the other three stages of the process need revisiting and considerable improvement, which was the focus of the author’s research presented in this review paper. For stage (ii), we consider the derivation of a multi-fluid model for suspension flow in a narrow vertical hydraulic fracture at moderate Re on the scale of fracture height and length and also the migration of particles across the flow on the scale of fracture width. At the stage of fracture cleanaup (iii), a novel multi-continua model for suspension filtration is developed. To provide closure relationships for permeability of proppant packings to be used in this model, a 3D direct numerical simulation of single phase flow is carried out using the lattice-Boltzmann method. For wellbore cleanup (iv), we present a combined 1D model for highly-transient gas-liquid flow based on the combination of multi-fluid and

On Cattaneo-Christov heat flux model for Carreau fluid flow over a slendering sheet

Science.gov (United States)

Hashim; Khan, Masood

The underlying intentions of this article are to investigate the impact of non-Fourier heat flux model on the stagnation-point flow of non-Newtonian Carreau fluid. In this study, the innovative Cattaneo-Christov constitutive model is introduced to study the characteristics of thermal relaxation time. The flow is impelled by a slendering surface which is of the variable thickness. In the model, the physical mechanism responsible for homogeneous-heterogeneous reactions are further taken into account. Also, the diffusion coefficients of the reactant and auto catalyst are considered to be equal. The governing non-linear partial differential equations consisting of the momentum, energy and concentration equations are reduced to the coupled ordinary differential equations by means of local similarity transformations. The transformed ODEs are tackled numerically by employing an effective shooting algorithm along with the Runge-Kutta Fehlberg scheme. The physical characteristics of the fluid velocity, temperature and concentration profiles are illuminated with the variation of numerous governing factors and are presented graphically. For instance, our result indicates that the temperature and thermal boundary layer thickness are lower in case of Cattaneo-Christov heat flux model when compared to classical Fourier's heat model. Meanwhile, the rate of heat transfer is significantly improved by a high wall thickness parameter and an opposite influence is found due to the thermal relaxation parameter. We further noticed that a higher value of homogeneous and heterogeneous reaction parameter corresponds to a deceleration in the concentration field and it shows an inverse relation for the Schmidt number. A correlation with accessible results for specific cases is found with fabulous consent.

Two-Phase Acto-Cytosolic Fluid Flow in a Moving Keratocyte: A 2D Continuum Model.

Science.gov (United States)

Nikmaneshi, M R; Firoozabadi, B; Saidi, M S

2015-09-01

The F-actin network and cytosol in the lamellipodia of crawling cells flow in a centripetal pattern and spout-like form, respectively. We have numerically studied this two-phase flow in the realistic geometry of a moving keratocyte. Cytosol has been treated as a low viscosity Newtonian fluid flowing through the high viscosity porous medium of F-actin network. Other involved phenomena including myosin activity, adhesion friction, and interphase interaction are also discussed to provide an overall view of this problem. Adopting a two-phase coupled model by myosin concentration, we have found new accurate perspectives of acto-cytosolic flow and pressure fields, myosin distribution, as well as the distribution of effective forces across the lamellipodia of a keratocyte with stationary shape. The order of magnitude method is also used to determine the contribution of forces in the internal dynamics of lamellipodia.

Two-fluid equilibria with flow

International Nuclear Information System (INIS)

Steinhauer, L.

1999-01-01

The formalism is developed for flowing two-fluid equilibria. The equilibrium system is governed by a pair of second order partial differential equations for the magnetic stream function and the ion stream function plus a Bernoulli-like equation for the density. There are six arbitrary surface function. There are separate characteristic surfaces for each species, which are the guiding-center surfaces. This system is a generalization of the familiar Grad-Shafranov system for a single-fluid equilibrium without flow, which has only one equation and two arbitrary surface functions. In the case of minimum energy equilibria, the six surface functions take on particular forms. (author)

Natural convection flow and heat transfer between a fluid layer and a porous layer inside a rectangular enclosure

International Nuclear Information System (INIS)

Beckermann, C.; Ramadhyani, S.; Viskanta, R.

1986-01-01

A numerical and experimental study is performed to analyze the steady-state natural convection fluid flow and heat transfer in a vertical rectangular enclosure that is partially filled with a vertical layer of a fluid-saturated porous medium. The flow in the porous layer is modeled utilizing the Brinkman-Forchheimer-extended Darcy equations. The numerical model is verified by conducting a number of experiments with spherical glass beads as the porous medium and water and glycerin as the fluids in rectangular test-cells. The agreement between the flow visualization results and temperature measurements and the numerical model is, in general, good. It is found that the amount of fluid penetrating from the fluid region into the porous layer depends strongly on the Darcy (Da) and Rayleigh (Ra) numbers. For a relatively low product of Ra x Da, the flow takes place primarily in the fluid layer, and heat transfer in the porous layer is by conduction only. On the other hand, fluid penetrating into a relatively highly permeable porous layer has a significant impact on the natural convection flow patterns in the entire enclosure

CFD simulation of IPR-R1 Triga subchannels fluid flow

International Nuclear Information System (INIS)

Silva, Vitor V.; Santos, A.; Mesquita, Amir Z.; Silva, P.S. da; Pereira, C.

2013-01-01

Computational fluid dynamics (CFD) codes have been extensively used in engineering problems, with increasing use in nuclear engineering. One of these computer codes is OpenFOAM. It is freely distributed with source code and offers a great flexibility in simulating particular conditions like those found in many problems in nuclear reactor analysis. The aim of this work is to simulate fluid flow and heat flux in three different configurations of subchannels of IPR-R1 TRIGA reactor using OpenFOAM. The data will be then validated against real experimental data obtained during the operation of the reactor at 100kW. This validation process is fundamental to allow the use of the software and associated model to simulate reactor's operation at different conditions, namely different power e fluid flow velocities. (author)

CFD simulation of IPR-R1 Triga subchannels fluid flow

Energy Technology Data Exchange (ETDEWEB)

Silva, Vitor V.; Santos, A.; Mesquita, Amir Z.; Silva, P.S. da, E-mail: vitors@cdtn.br, E-mail: aacs@cdtn.br, E-mail: amir@cdtn.br, E-mail: psblsg@cdtn.br [Centro de Desenvolvimento da Tecnologia Nuclear (CDTN/CNEN - MG), Belo Horizonte, MG (Brazil); Pereira, C., E-mail: claubia@nuclear.ufmg.br [Universidade Federal de Minas Gerais (UFMG), Belo Horizonte, MG (Brazil). Dept. de Engenharia Nuclear

2013-07-01

Computational fluid dynamics (CFD) codes have been extensively used in engineering problems, with increasing use in nuclear engineering. One of these computer codes is OpenFOAM. It is freely distributed with source code and offers a great flexibility in simulating particular conditions like those found in many problems in nuclear reactor analysis. The aim of this work is to simulate fluid flow and heat flux in three different configurations of subchannels of IPR-R1 TRIGA reactor using OpenFOAM. The data will be then validated against real experimental data obtained during the operation of the reactor at 100kW. This validation process is fundamental to allow the use of the software and associated model to simulate reactor's operation at different conditions, namely different power e fluid flow velocities. (author)

Application of x-ray microtomography to environmental fluid flow problems

International Nuclear Information System (INIS)

Wildenschild, D.; Culligan, K.A.; Christensen, B.S.B.

2005-01-01

Many environmental processes are controlled by the micro-scale interaction of water and air with the solid phase (soils, sediments, rock) in pore spaces within the subsurface. The distribution in time and space of fluids in pores ultimately controls subsurface flow and contaminant transport relevant to groundwater resource management, contaminant remediation, and agriculture. Many of these physical processes operative at the pore-scale cannot be directly investigated using conventional hydrologic techniques, however recent developments in synchrotron-based micro-imaging have made it possible to observe and quantify pore-scale processes non-invasively. Micron-scale resolution makes it possible to track fluid flow within individual pores and therefore facilitates previously unattainable measurements. We report on experiments performed at the GSECARS** (Advanced Photon Source) microtomography facility and have measured properties such as porosity, fluid saturation and distribution within the pore space, as well as interfacial characteristics of the fluids involved (air, water, contaminant). Different image processing techniques were applied following mathematical reconstruction to produce accurate measurements of the physical flow properties. These new micron-scale measurements make it possible to test existing and new theory, as well as emerging numerical modeling schemes aimed at the pore scale.

Bone tissue engineering: the role of interstitial fluid flow

Science.gov (United States)

Hillsley, M. V.; Frangos, J. A.

1994-01-01

It is well established that vascularization is required for effective bone healing. This implies that blood flow and interstitial fluid (ISF) flow are required for healing and maintenance of bone. The fact that changes in bone blood flow and ISF flow are associated with changes in bone remodeling and formation support this theory. ISF flow in bone results from transcortical pressure gradients produced by vascular and hydrostatic pressure, and mechanical loading. Conditions observed to alter flow rates include increases in venous pressure in hypertension, fluid shifts occurring in bedrest and microgravity, increases in vascularization during the injury-healing response, and mechanical compression and bending of bone during exercise. These conditions also induce changes in bone remodeling. Previously, we hypothesized that interstitial fluid flow in bone, and in particular fluid shear stress, serves to mediate signal transduction in mechanical loading- and injury-induced remodeling. In addition, we proposed that a lack or decrease of ISF flow results in the bone loss observed in disuse and microgravity. The purpose of this article is to review ISF flow in bone and its role in osteogenesis.

Fluid flow in 0.5-m scale blocks of Topopah Spring tuff

International Nuclear Information System (INIS)

Blair, S. C.; Carlson, S. R.; Constantino, M. S.

1999-01-01

A laboratory experiment was conducted on a 0.5-m scale block of Topopah Spring tuff, to measure fluid flow and mechanical deformation properties under conditions that approximate the near-field environment of a potential nuclear waste repository, and to provide an intermediate-scale test case for numerical model validation. The test specimen is a 0.25 x 0.25 x 0.50 m rectangular prism bisected by an artificial (saw-cut) fracture orthogonal to the tuff fabric. Water was supplied by a point source at the center of the fracture under various pressures of up to 0.04 MPa. Both fluid flow and mechanical properties were found to be anisotropic and strongly correlated with the ash flow fabric. Fluid mass-balance measurements revealed that only minor imbibition of water occurred through the fracture surfaces and that flow rates were independent of normal stress to 14.0 MPa and temperature to 140 C. Flow through the fracture occurred largely through uncorrelated porosity that intersected the fracture plane

Computer modeling of fluid flow and combustion in the ISV [In Situ Vitrification] confinement hood

International Nuclear Information System (INIS)

Johnson, R.W.; Paik, S.

1990-09-01

Safety and suitability objectives for the application of the In Situ Vitrification (ISV) technology at the INEL require that the physical processes involved in ISVV be modeled to determine their operational behavior. The mathematical models that have been determined to address the modeling needs adequately for the ISV analysis package are detailed elsewhere. The present report is concerned with the models required for simulating the reacting flow that occurs in the ISV confinement hood. An experimental code named COYOTE has been secured that appears adequate to model the combustion in the confinement hood. The COYOTE code is a two-dimensional, transient, compressible, Eulerian, gas dynamics code for modeling reactive flows. It recognizes nonuniform Cartesian and cylindrical geometry and is based on the ICE (Implicit Continuous-fluid Eulerian) family of solution methods. It includes models for chemical reactions based on chemical kinetics as well as equilibrium chemistry. The mathematical models contained in COYOTE, their discrete analogs, the solution procedure, code structure and some test problems are presented in the report. 12 refs., 17 figs., 6 tabs

Scaling of two-phase flow transients using reduced pressure system and simulant fluid

International Nuclear Information System (INIS)

Kocamustafaogullari, G.; Ishii, M.

1987-01-01

Scaling criteria for a natural circulation loop under single-phase flow conditions are derived. Based on these criteria, practical applications for designing a scaled-down model are considered. Particular emphasis is placed on scaling a test model at reduced pressure levels compared to a prototype and on fluid-to-fluid scaling. The large number of similarty groups which are to be matched between modell and prototype makes the design of a scale model a challenging tasks. The present study demonstrates a new approach to this clasical problen using two-phase flow scaling parameters. It indicates that a real time scaling is not a practical solution and a scaled-down model should have an accelerated (shortened) time scale. An important result is the proposed new scaling methodology for simulating pressure transients. It is obtained by considerung the changes of the fluid property groups which appear within the two-phase similarity parameters and the single-phase to two-phase flow transition prameters. Sample calculations are performed for modeling two-phase flow transients of a high pressure water system by a low-pressure water system or a Freon system. It is shown that modeling is possible for both cases for simulation pressure transients. However, simulation of phase change transitions is not possible by a reduced pressure water system without distortion in either power or time. (orig.)

Electrification of particulate entrained fluid flows-Mechanisms, applications, and numerical methodology

Science.gov (United States)

Wei, Wei; Gu, Zhaolin

2015-10-01

Particulates in natural and industrial flows have two basic forms: liquid (droplet) and solid (particle). Droplets would be charged in the presence of the applied electric field (e.g. electrospray). Similar to the droplet charging, particles can also be charged under the external electric field (e.g. electrostatic precipitator), while in the absence of external electric field, tribo-electrostatic charging is almost unavoidable in gas-solid two-phase flows due to the consecutive particle contacts (e.g. electrostatic in fluidized bed or wind-blown sand). The particle charging may be beneficial, or detrimental. Although electrostatics in particulate entrained fluid flow systems have been so widely used and concerned, the mechanisms of particulate charging are still lack of a thorough understanding. The motivation of this review is to explore a clear understanding of particulate charging and movement of charged particulate in two-phase flows, by summarizing the electrification mechanisms, physical models of particulate charging, and methods of charging/charged particulate entrained fluid flow simulations. Two effective methods can make droplets charged in industrial applications: corona charging and induction charging. The droplet charge to mass ratio by corona charging is more than induction discharge. The particle charging through collisions could be attributed to electron transfer, ion transfer, material transfer, and/or aqueous ion shift on particle surfaces. The charges on charged particulate surface can be measured, nevertheless, the charging process in nature or industry is difficult to monitor. The simulation method might build a bridge of investigating from the charging process to finally charged state on particulate surface in particulate entrained fluid flows. The methodology combining the interface tracking under the action of the applied electric with the fluid flow governing equations is applicable to the study of electrohydrodynamics problems. The charge

Topological Fluid Mechanics with Applications to Free Surfaces and Axisymmetric Flows

DEFF Research Database (Denmark)

Brøns, Morten

1996-01-01

Topological fluid mechanics is the study of qualitative features of fluid patterns. We discuss applications to the flow beneath a stagnant surface film, and to patterns in axisymmetric flow.......Topological fluid mechanics is the study of qualitative features of fluid patterns. We discuss applications to the flow beneath a stagnant surface film, and to patterns in axisymmetric flow....

Acceleration of coupled granular flow and fluid flow simulations in pebble bed energy systems

Energy Technology Data Exchange (ETDEWEB)

Li, Yanheng, E-mail: liy19@rpi.edu [Department of Mechanical, Aerospace, and Nuclear Engineering, Rensselaer Polytechnic Institute, 110 8th Street, Troy, NY (United States); Ji, Wei, E-mail: jiw2@rpi.edu [Department of Mechanical, Aerospace, and Nuclear Engineering, Rensselaer Polytechnic Institute, 110 8th Street, Troy, NY (United States)

2013-05-15

Highlights: ► Fast simulation of coupled pebble flow and coolant flow in PBR systems is studied. ► Dimension reduction based on axisymmetric geometry shows significant speedup. ► Relaxation of coupling frequency is investigated and an optimal range is determined. ► A total of 80% efficiency increase is achieved by the two fast strategies. ► Fast strategies can be applied to simulating other general fluidized bed systems. -- Abstract: Fast and accurate approaches to simulating the coupled particle flow and fluid flow are of importance to the analysis of large particle-fluid systems. This is especially needed when one tries to simulate pebble flow and coolant flow in Pebble Bed Reactor (PBR) energy systems on a routine basis. As one of the Generation IV designs, the PBR design is a promising nuclear energy system with high fuel performance and inherent safety. A typical PBR core can be modeled as a particle-fluid system with strong interactions among pebbles, coolants and reactor walls. In previous works, the coupled Discrete Element Method (DEM)-Computational Fluid Dynamics (CFD) approach has been investigated and applied to modeling PBR systems. However, the DEM-CFD approach is computationally expensive due to large amounts of pebbles in PBR systems. This greatly restricts the PBR analysis for the real time prediction and inclusion of more physics. In this work, based on the symmetry of the PBR geometry and the slow motion characteristics of the pebble flow, two acceleration strategies are proposed. First, a simplified 3D-DEM/2D-CFD approach is proposed to speed up the DEM-CFD simulation without loss of accuracy. Pebble flow is simulated by a full 3D DEM, while the coolant flow field is calculated with a 2D CFD simulation by averaging variables along the annular direction in the cylindrical and annular geometries. Second, based on the slow motion of pebble flow, the impact of the coupling frequency on the computation accuracy and efficiency is

Acceleration of coupled granular flow and fluid flow simulations in pebble bed energy systems

International Nuclear Information System (INIS)

Li, Yanheng; Ji, Wei

2013-01-01

Highlights: ► Fast simulation of coupled pebble flow and coolant flow in PBR systems is studied. ► Dimension reduction based on axisymmetric geometry shows significant speedup. ► Relaxation of coupling frequency is investigated and an optimal range is determined. ► A total of 80% efficiency increase is achieved by the two fast strategies. ► Fast strategies can be applied to simulating other general fluidized bed systems. -- Abstract: Fast and accurate approaches to simulating the coupled particle flow and fluid flow are of importance to the analysis of large particle-fluid systems. This is especially needed when one tries to simulate pebble flow and coolant flow in Pebble Bed Reactor (PBR) energy systems on a routine basis. As one of the Generation IV designs, the PBR design is a promising nuclear energy system with high fuel performance and inherent safety. A typical PBR core can be modeled as a particle-fluid system with strong interactions among pebbles, coolants and reactor walls. In previous works, the coupled Discrete Element Method (DEM)-Computational Fluid Dynamics (CFD) approach has been investigated and applied to modeling PBR systems. However, the DEM-CFD approach is computationally expensive due to large amounts of pebbles in PBR systems. This greatly restricts the PBR analysis for the real time prediction and inclusion of more physics. In this work, based on the symmetry of the PBR geometry and the slow motion characteristics of the pebble flow, two acceleration strategies are proposed. First, a simplified 3D-DEM/2D-CFD approach is proposed to speed up the DEM-CFD simulation without loss of accuracy. Pebble flow is simulated by a full 3D DEM, while the coolant flow field is calculated with a 2D CFD simulation by averaging variables along the annular direction in the cylindrical and annular geometries. Second, based on the slow motion of pebble flow, the impact of the coupling frequency on the computation accuracy and efficiency is

Computational fluid dynamics modeling of two-phase flow in a BWR fuel assembly. Final CRADA Report

International Nuclear Information System (INIS)

Tentner, A.

2009-01-01

A direct numerical simulation capability for two-phase flows with heat transfer in complex geometries can considerably reduce the hardware development cycle, facilitate the optimization and reduce the costs of testing of various industrial facilities, such as nuclear power plants, steam generators, steam condensers, liquid cooling systems, heat exchangers, distillers, and boilers. Specifically, the phenomena occurring in a two-phase coolant flow in a BWR (Boiling Water Reactor) fuel assembly include coolant phase changes and multiple flow regimes which directly influence the coolant interaction with fuel assembly and, ultimately, the reactor performance. Traditionally, the best analysis tools for this purpose of two-phase flow phenomena inside the BWR fuel assembly have been the sub-channel codes. However, the resolution of these codes is too coarse for analyzing the detailed intra-assembly flow patterns, such as flow around a spacer element. Advanced CFD (Computational Fluid Dynamics) codes provide a potential for detailed 3D simulations of coolant flow inside a fuel assembly, including flow around a spacer element using more fundamental physical models of flow regimes and phase interactions than sub-channel codes. Such models can extend the code applicability to a wider range of situations, which is highly important for increasing the efficiency and to prevent accidents.

Finite element modeling of melting and fluid flow in the laser-heated diamond-anvil cell

Science.gov (United States)

Gomez-Perez, N.; Rodriguez, J. F.; McWilliams, R. S.

2017-04-01

The laser-heated diamond anvil cell is widely used in the laboratory study of materials behavior at high-pressure and high-temperature, including melting curves and liquid properties at extreme conditions. Laser heating in the diamond cell has long been associated with fluid-like motion in samples, which is routinely used to determine melting points and is often described as convective in appearance. However, the flow behavior of this system is poorly understood. A quantitative treatment of melting and flow in the laser-heated diamond anvil cell is developed here to physically relate experimental motion to properties of interest, including melting points and viscosity. Numerical finite-element models are used to characterize the temperature distribution, melting, buoyancy, and resulting natural convection in samples. We find that continuous fluid motion in experiments can be explained most readily by natural convection. Fluid velocities, peaking near values of microns per second for plausible viscosities, are sufficiently fast to be detected experimentally, lending support to the use of convective motion as a criterion for melting. Convection depends on the physical properties of the melt and the sample geometry and is too sluggish to detect for viscosities significantly above that of water at ambient conditions, implying an upper bound on the melt viscosity of about 1 mPa s when convective motion is detected. A simple analytical relationship between melt viscosity and velocity suggests that direct viscosity measurements can be made from flow speeds, given the basic thermodynamic and geometric parameters of samples are known.

A computational analysis on homogeneous-heterogeneous mechanism in Carreau fluid flow

Science.gov (United States)

Khan, Imad; Rehman, Khalil Ur; Malik, M. Y.; Shafquatullah

2018-03-01

In this article magnetohydrodynamic Carreau fluid flow towards stretching cylinder is considered in the presence of homogeneous-heterogeneous reactions effect. The flow model is structured by utilizing theoretical grounds. For the numerical solution a shooting method along with Runge-Kutta algorithm is executed. The outcomes are provided through graphs. It is observed that the Carreau fluid concentration shows decline values via positive iterations of homogeneous-heterogeneous reaction parameters towards both shear thinning and thickening case. The present work is certified through comparison with already existing literature in a limiting sense.

Approximate Solution of Dam-break Flow of Low Viscosity Bingham Fluid

Science.gov (United States)

Puay, How Tion; Hosoda, Takashi

In this study, we investigate the characteristics of dam-break flow of low viscosity Bingham fluid by deriving an approximate solution for the time development of the front position and depth at the origin of the flow. The asymptotic solutions representing the characteristic of Bingham fluid in the limit of low plastic viscosity are verified with a depth-averaged numerical model. Numerical simulations showed that with the decrease of plastic viscosity, the time development of the front position and depth at the origin approach to the theoretical asymptotic solution.

Modelling of fluid flow in fractured porous media by the singular integral equations method

International Nuclear Information System (INIS)

Vu, M.N.

2012-01-01

This thesis aims to develop a method for numerical modelling of fluid flow through fractured porous media and for determination of their effective permeability by taking advantage of recent results based on formulation of the problem by Singular Integral Equations. In parallel, it was also an occasion to continue on the theoretical development and to obtain new results in this area. The governing equations for flow in such materials are reviewed first and mass conservation at the fracture intersections is expressed explicitly. Using the theory of potential, the general potential solutions are proposed in the form of a singular integral equation that describes the steady-state flow in and around several fractures embedded in an infinite porous matrix under a far-field pressure condition. These solutions represent the pressure field in the whole body as functions of the infiltration in the fractures, which fully take into account the fracture interaction and intersections. Closed-form solutions for the fundamental problem of fluid flow around a single fracture are derived, which are considered as the benchmark problems to validate the numerical solutions. In particular, the solution obtained for the case of an elliptical disc-shaped crack obeying to the Poiseuille law has been compared to that obtained for ellipsoidal inclusions with Darcy law.The numerical programs have been developed based on the singular integral equations method to resolve the general potential equations. These allow modeling the fluid flow through a porous medium containing a great number of fractures. Besides, this formulation of the problem also allows obtaining a semi-analytical infiltration solution over a single fracture depending on the matrice permeability, the fracture conductivity and the fracture geometry. This result is the important key to up-scaling the effective permeability of a fractured porous medium by using different homogenisation schemes. The results obtained by the self

Linear analysis of sheared flow stabilization of global magnetohydrodynamic instabilities based on the Hall fluid model

International Nuclear Information System (INIS)

Sotnikov, V.I.; Paraschiv, I.; Makhin, V.; Bauer, B.S.; Leboeuf, J.N.; Dawson, J.M.

2002-01-01

A systematic study of the linear stage of sheared flow stabilization of Z-pinch plasmas based on the Hall fluid model with equilibrium that contains sheared flow and an axial magnetic field is presented. In the study we begin with the derivation of a general set of equations that permits the evaluation of the combined effect of sheared flow and axial magnetic field on the development of the azimuthal mode number m=0 sausage and m=1 kink magnetohydrodynamic (MHD) instabilities, with the Hall term included in the model. The incorporation of sheared flow, axial magnetic field, and the Hall term allows the Z-pinch system to be taken away from the region in parameter space where ideal MHD is applicable to a regime where nonideal effects tend to govern stability. The problem is then treated numerically by following the linear development in time of an initial perturbation. The numerical results for linear growth rates as a function of axial sheared flow, an axial magnetic field, and the Hall term are reported

Biorheological Model on Flow of Herschel-Bulkley Fluid through a Tapered Arterial Stenosis with Dilatation.

Science.gov (United States)

Priyadharshini, S; Ponalagusamy, R

2015-01-01

An analysis of blood flow through a tapered artery with stenosis and dilatation has been carried out where the blood is treated as incompressible Herschel-Bulkley fluid. A comparison between numerical values and analytical values of pressure gradient at the midpoint of stenotic region shows that the analytical expression for pressure gradient works well for the values of yield stress till 2.4. The wall shear stress and flow resistance increase significantly with axial distance and the increase is more in the case of converging tapered artery. A comparison study of velocity profiles, wall shear stress, and flow resistance for Newtonian, power law, Bingham-plastic, and Herschel-Bulkley fluids shows that the variation is greater for Herschel-Bulkley fluid than the other fluids. The obtained velocity profiles have been compared with the experimental data and it is observed that blood behaves like a Herschel-Bulkley fluid rather than power law, Bingham, and Newtonian fluids. It is observed that, in the case of a tapered stenosed tube, the streamline pattern follows a convex pattern when we move from r/R = 0 to r/R = 1 and it follows a concave pattern when we move from r/R = 0 to r/R = -1. Further, it is of opposite behaviour in the case of a tapered dilatation tube which forms new information that is, for the first time, added to the literature.

Fluid-elastic vibration in two-phase cross flow

International Nuclear Information System (INIS)

Sasakawa, T.; Serizawa, A.; Kawara, Z.

2003-01-01

The present work aims at clarifying the mechanisms of fluid elastic vibration of tube bundles in two-phase cross flow. The experiment is conducted using air-water two-phase flow under atmospheric pressure. The test section is a 1.03m long transparent acrylic square duct with 128 x 128 mm 2 cross section, which consists of 3 rod-rows with 5 rods in each row. The rods are 125mm long aluminum rods with 22 mm in diameter (p/D=1.45). The natural frequency of rod vibration is about 30Hz. The result indicated a diversion of observed trend in vibration behavior depending on two-phase flow patterns either bubbly flow or churn flow. Specifically, in churn flow, the fluid elastic vibration has been observed to occur when the frequency in void fraction fluctuation approached to the natural frequency of the rods, but this was not the case in fluid elastic vibration in bubbly flow. This fact suggests the existence of mechanisms closely coupled with two-phase flow structures depending on the flow patterns, that is, static two-phase character-controlled mechanism in bubbly flow and dynamic character- controlled in churn flow

Newtonian heating effect on unsteady hydromagnetic Casson fluid flow past a flat plate with heat and mass transfer

Directory of Open Access Journals (Sweden)

M. Das

2015-12-01

Full Text Available The influence of Newtonian heating on heat and mass transfer in unsteady hydromagnetic flow of a Casson fluid past a vertical plate in the presence of thermal radiation and chemical reaction is studied. The Casson fluid model is used to distinguish the non-Newtonian fluid behavior. The fluid flow is induced due to periodic oscillations of the plate along its length and a uniform transverse magnetic field is applied in a direction which is normal to the direction of fluid flow. The partial differential equations governing the flow, heat, and mass transfer are transformed to non-dimensional form using suitable non-dimensional variables which are then solved analytically by using Laplace transform technique. The numerical values of the fluid velocity, fluid temperature, and species concentration are depicted graphically whereas the values of skin-friction, Nusselt number, and Sherwood number are presented in tabular form. It is noticed that the fluid velocity and temperature decrease with increasing values of Casson parameter while concentration decreases with increasing values of chemical reaction parameter and Schmidt number. Such a fluid flow model has several industrial and medical applications such as in glass manufacturing, paper production, purification of crude oil and study of blood flow in the cardiovascular system.

Modelling and simulation of compressible fluid flow in oil reservoir: a case study of the Jubilee Field, Tano Basin (Ghana)

International Nuclear Information System (INIS)

Gawusu, S.

2015-07-01

Oil extraction represents an important investment and the control of a rational exploitation of a field means mastering various scientific techniques including the understanding of the dynamics of fluids in place. This thesis presents a theoretical investigation of the dynamic behaviour of an oil reservoir during its exploitation. The study investigated the dynamics of fluid flow patterns in a homogeneous oil reservoir using the Radial Diffusivity Equation (RDE) as well as two phase oil-water flow equations. The RDE model was solved analytically and numerically for pressure using the Constant Terminal Rate Solution (CTRS) and the fully implicit Finite Difference Method (FDM) respectively. The mathematical derivations of the models and their solution procedures were presented to allow for easy utilization of the techniques for reservoir and engineering applications. The study predicted that the initial oil reservoir pressure will be able to do the extraction for a very long time before any other recovery method will be used to aid in the extraction process depending on the rate of production. Reservoir simulation describing a one dimensional radial flow of a compressible fluid in porous media may be adequately performed using ordinary laptop computers as revealed by the study. For the simulation of MATLAB, the case of the Jubilee Fields, Tano Basin was studied, an algorithm was developed for the simulation of pressure in the reservoir. It ensues from the analysis of the plots of pressure vrs time and space that the Pressure Transient Analysis (PTA) was duly followed. The approximate solutions of the analytical and numerical solutions to the Radial Diffusivity Equation (RDE) were in excellent agreement, thus the reservoir simulation model developed can be used to describe typical pressure-time relationships that are used in conventional Pressure Transient Analysis (PTA). The study was extended to two phase oil-water flow in reservoirs. The flow of fluids in multi

Fluid flow and heat transfer investigation of pebble bed reactors using mesh adaptive large-eddy simulation

International Nuclear Information System (INIS)

Pavlidis, D.; Lathouwers, D.

2011-01-01

A computational fluid dynamics model with anisotropic mesh adaptivity is used to investigate coolant flow and heat transfer in pebble bed reactors. A novel method for implicitly incorporating solid boundaries based on multi-fluid flow modelling is adopted. The resulting model is able to resolve and simulate flow and heat transfer in randomly packed beds, regardless of the actual geometry, starting off with arbitrarily coarse meshes. The model is initially evaluated using an orderly stacked square channel of channel-height-to-particle diameter ratio of unity for a range of Reynolds numbers. The model is then applied to the face-centred cubical geometry. Coolant flow and heat transfer patterns are investigated. (author)

Fluid Damping Variation of a Slender Rod in Axial Flow Field

Energy Technology Data Exchange (ETDEWEB)

Park, Nam-Gyu; Yoo, Jong-Sung; Jung, Yil-Sup [KEPCO Nuclear Fuel Co., Daejeon (Korea, Republic of)

2016-10-15

This study proposed an analytic damping model considering the axial flow condition. In addition, the specific damping values with respect to the flow speeds are calculated. The flow induced damping is beneficial to fuel integrity in that impact energy due to severe accidents such as earthquake dissipates rapidly. A nuclear fuel bundle is composed of many slender fuel rods which contain fission material. The slender rod is typical structure in the fuel, therefore fluid damping estimation on the rod should be an important clue leading to fuel bundle damping identification. Severe accidents could cause fuel assembly vibration in the core, but large motion could be damped out rapidly when a strong damping mechanism is involved. This paper suggested a mathematical model of the slender structure. The physical meaning of the model is described, and the simulation results with the model are also provided. Actual damping due to the fluid is nonlinear, therefore further works are required to explain the detail behavior with the nonlinearity. The model validation test is on-going in KEPCO Nuclear Fuel, but it is believed that performance of the model is well correlated to the published work.

Simulations of flow induced ordering in viscoelastic fluids

NARCIS (Netherlands)

Santos de Oliveira, I.S.

2012-01-01

In this thesis we report on simulations of colloidal ordering phenomena in shearthinning viscoelastic fluids under shear flow. Depending on the characteristics of the fluid, the colloids are observed to align in the direction of the flow. These string-like structures remain stable as long as the

Analysis of the flow dynamics characteristics of an axial piston pump based on the computational fluid dynamics method

Directory of Open Access Journals (Sweden)

Bin Zhang

2017-01-01

Full Text Available To improve its working performance, the flow ripple characteristics of an axial piston pump were investigated with software which uses computational fluid dynamics (CFD technology. The simulation accuracy was significantly optimized through the use of the improved compressible fluid model. Flow conditions of the pump were tested using a pump flow ripple test rig, and the simulation results of the CFD model showed good agreement with the experimental data. Additionally, the composition of the flow ripple was analyzed using the improved CFD model, and the results showed that the compression ripple makes up 88% of the flow ripple. The flow dynamics of the piston pump is mainly caused by the pressure difference between the intake and discharge ports of the valve plates and the fluid oil compressibility.

The influence of fluid - flexible particle interaction on fluid flow optical non-homogeneity in channel bifurcation

Science.gov (United States)

Tazyukov, F. Kh; Kutuzova, E. R.; Garifullin, F. A.

2017-01-01

In the present paper the peculiar properties of convergent fluid flow in T-junction channel is considered. There is no interaction between flexible particles in the flow. Such kind of situation is described by rheological FENE-P and Oldroyd-B models. The first one predicts viscosity anomaly, dependence of longitudinal viscosity on longitudinal strain rate and elastic properties; the last one - existence of longitudinal viscosity depending on longitudinal strain rate and having a physical sense only for and elastic properties. The model’s governing parameters are the Weissenberg number (We), the Reynolds number (Re), the ability of flexible particle to change its orientation and stretching degree (L2 ) in the main flow. The bifurcation area is of great importance due to possibility of high stresses and velocities existence not only in central area, but also on the walls and near the corners. The symmetry-loss effect at creeping flows regime (Re≪1) is investigated. It has been showed that at certain set of We and L2 values the symmetrical shape of fluid flow turns to asymmetrical shape.

The influence of fluid - flexible particle interaction on fluid flow optical non-homogeneity in channel bifurcation

International Nuclear Information System (INIS)

Tazyukov, F Kh; Kutuzova, E R; Garifullin, F A

2017-01-01

In the present paper the peculiar properties of convergent fluid flow in T-junction channel is considered. There is no interaction between flexible particles in the flow. Such kind of situation is described by rheological FENE-P and Oldroyd-B models. The first one predicts viscosity anomaly, dependence of longitudinal viscosity on longitudinal strain rate and elastic properties; the last one – existence of longitudinal viscosity depending on longitudinal strain rate and having a physical sense only for and elastic properties. The model’s governing parameters are the Weissenberg number ( We ), the Reynolds number ( Re ), the ability of flexible particle to change its orientation and stretching degree ( L 2 ) in the main flow. The bifurcation area is of great importance due to possibility of high stresses and velocities existence not only in central area, but also on the walls and near the corners. The symmetry-loss effect at creeping flows regime ( Re ≪1) is investigated. It has been showed that at certain set of We and L 2 values the symmetrical shape of fluid flow turns to asymmetrical shape. (paper)

Fluid flow and heat transfer investigation of pebble bed reactors using mesh-adaptive LES

International Nuclear Information System (INIS)

Pavlidis, Dimitrios; Lathouwers, Danny

2013-01-01

The very high temperature reactor is one of the designs currently being considered for nuclear power generation. One its variants is the pebble bed reactor in which the coolant passes through complex geometries (pores) at high Reynolds numbers. A computational fluid dynamics model with anisotropic mesh adaptivity is used to investigate coolant flow and heat transfer in such reactors. A novel method for implicitly incorporating solid boundaries based on multi-fluid flow modelling is adopted. The resulting model is able to resolve and simulate flow and heat transfer in randomly packed beds, regardless of the actual geometry, starting off with arbitrarily coarse meshes. The model is initially evaluated using an orderly stacked square channel of channel-height-to-particle diameter ratio of unity for a range of Reynolds numbers. The model is then applied to the face-centred cubical geometry. coolant flow and heat transfer patterns are investigated

Development of the tube bundle structure for fluid-structure interaction analysis model

International Nuclear Information System (INIS)

Yoon, Kyung Ho; Kim, Jae Yong

2010-02-01

Tube bundle structures within a Boiler or heat exchanger are laid the fluid-structure, thermal-structure and fluid-thermal-structure coupled boundary condition. In these complicated boundary conditions, Fluid-structure interaction (FSI) occurs when fluid flow causes deformation of the structure. This deformation, in turn, changes the boundary conditions for the fluid flow. The structural analysis discipline, and then independently analyzed each other. However, the fluid dynamic force effect the behavior of the structure, and the vibration amplitude of the structure to fluid. FSI analysis model was separately created fluid and structure model, and then defined the fsi boundary condition, and simultaneously analyzed in one domain. The analysis results were compared with those of the experimental method for validating the analysis model. Flow-induced vibration test was executed with single rod configuration. The vibration amplitudes of a fuel rod were measured by the laser vibro-meter system in x and y-direction. The analyses results were not closely with the test data, but the trend was very similar with the test result. In fsi coupled analysis case, the turbulent model was very important with the reliability of the accuracy of the analysis model. Therefore, the analysis model will be needed to further study

Topological fluid dynamics of interfacial flows

DEFF Research Database (Denmark)

Brøns, Morten

1994-01-01

The topological description of flows in the vicinity of a solid boundary, that is familiar from the aerodynamics literature, has recently been extended to the case of flow at a liquid–gas interface or a free surface by Lugt [Phys. Fluids 30, 3647 (1987)]. Lugt's work is revisited in a more general...... setting, including nonconstant curvature of the interface and gradients of surface tension, using tools of modern nonlinear dynamics. Bifurcations of the flow pattern occur at degenerate configurations. Using the theory of unfolding, this paper gives a complete description of the bifurcations that depend...... on terms up to the second order. The general theory of this paper is applied to the topology of streamlines during the breaking of a wave and to the flow below a stagnant surface film. Physics of Fluids is copyrighted by The American Institute of Physics....

Investigation of Swirling Flow in Rod Bundle Subchannels Using Computational Fluid Dynamics

International Nuclear Information System (INIS)

Holloway, Mary V.; Beasley, Donald E.; Conner, Michael E.

2006-01-01

The fluid dynamics for turbulent flow through rod bundles representative of those used in pressurized water reactors is examined using computational fluid dynamics (CFD). The rod bundles of the pressurized water reactor examined in this study consist of a square array of parallel rods that are held on a constant pitch by support grids spaced axially along the rod bundle. Split-vane pair support grids are often used to create swirling flow in the rod bundle in an effort to improve the heat transfer characteristics for the rod bundle during both normal operating conditions and in accident condition scenarios. Computational fluid dynamics simulations for a two subchannel portion of the rod bundle were used to model the flow downstream of a split-vane pair support grid. A high quality computational mesh was used to investigate the choice of turbulence model appropriate for the complex swirling flow in the rod bundle subchannels. Results document a central swirling flow structure in each of the subchannels downstream of the split-vane pairs. Strong lateral flows along the surface of the rods, as well as impingement regions of lateral flow on the rods are documented. In addition, regions of lateral flow separation and low axial velocity are documented next to the rods. Results of the CFD are compared to experimental particle image velocimetry (PIV) measurements documenting the lateral flow structures downstream of the split-vane pairs. Good agreement is found between the computational simulation and experimental measurements for locations close to the support grid. (authors)

Modal analysis of fluid flows using variants of proper orthogonal decomposition

Science.gov (United States)

Rowley, Clarence; Dawson, Scott

2017-11-01

This talk gives an overview of several methods for analyzing fluid flows, based on variants of proper orthogonal decomposition. These methods may be used to determine simplified, approximate models that capture the essential features of these flows, in order to better understand the dominant physical mechanisms, and potentially to develop appropriate strategies for model-based flow control. We discuss balanced proper orthogonal decomposition as an approximation of balanced truncation, and explain connections with system identification methods such as the eigensystem realization algorithm. We demonstrate the methods on several canonical examples, including a linearized channel flow and the flow past a circular cylinder. Supported by AFOSR, Grant FA9550-14-1-0289.

Free surface modelling with two-fluid model and reduced numerical diffusion of the interface

International Nuclear Information System (INIS)

Strubelj, Luka; Tiselj, Izrok

2008-01-01

Full text of publication follows: The free surface flows are successfully modelled with one of existing free surface models, such as: level set method, volume of fluid method (with/without surface reconstruction), front tracking, two-fluid model (two momentum equations) with modified interphase force and others. The main disadvantage of two-fluid model used for simulations of free surface flows is numerical diffusion of the interface, which can be significantly reduced using the method presented in this paper. Several techniques for reduction of numerical diffusion of the interface have been implemented in the volume of fluid model and are based on modified numerical schemes for advection of volume fraction near the interface. The same approach could be used also for two-fluid method, but according to our experience more successful reduction of numerical diffusion of the interface can be achieved with conservative level set method. Within the conservative level set method, continuity equation for volume fraction is solved and after that the numerical diffusion of the interface is reduced in such a way that the thickness of the interface is kept constant during the simulation. Reduction of the interface diffusion can be also called interface sharpening. In present paper the two-fluid model with interface sharpening is validated on Rayleigh-Taylor instability. Under assumptions of isothermal and incompressible flow of two immiscible fluids, we simulated a system with the fluid of higher density located above the fluid of smaller density in two dimensions. Due to gravity in the system, fluid with higher density moves below the fluid with smaller density. Initial condition is not a flat interface between the fluids, but a sine wave with small amplitude, which develops into a mushroom-like structure. Mushroom-like structure in simulation of Rayleigh-Taylor instability later develops to small droplets as result of numerical dispersion of interface (interface sharpening

Electroosmotic flow of Phan-Thien-Tanner fluids at high zeta potentials: An exact analytical solution

Science.gov (United States)

Sarma, Rajkumar; Deka, Nabajit; Sarma, Kuldeep; Mondal, Pranab Kumar

2018-06-01

We present a mathematical model to study the electroosmotic flow of a viscoelastic fluid in a parallel plate microchannel with a high zeta potential, taking hydrodynamic slippage at the walls into account in the underlying analysis. We use the simplified Phan-Thien-Tanner (s-PTT) constitutive relationships to describe the rheological behavior of the viscoelastic fluid, while Navier's slip law is employed to model the interfacial hydrodynamic slip. Here, we derive analytical solutions for the potential distribution, flow velocity, and volumetric flow rate based on the complete Poisson-Boltzmann equation (without considering the frequently used Debye-Hückel linear approximation). For the underlying electrokinetic transport, this investigation primarily reveals the influence of fluid rheology, wall zeta potential as modulated by the interfacial electrochemistry and interfacial slip on the velocity distribution, volumetric flow rate, and fluid stress, as well as the apparent viscosity. We show that combined with the viscoelasticity of the fluid, a higher wall zeta potential and slip coefficient lead to a phenomenal enhancement in the volumetric flow rate. We believe that this analysis, besides providing a deep theoretical insight to interpret the transport process, will also serve as a fundamental design tool for microfluidic devices/systems under electrokinetic influence.

On exact solutions for oscillatory flows in a generalized Burgers fluid with slip condition

Energy Technology Data Exchange (ETDEWEB)

Hayat, Tasawar [Dept. of Mathematics, Quaid-i-Azam Univ., Islamabad (Pakistan); Dept. of Mathematics, Coll. of Sciences, KS Univ., Riyadh (Saudi Arabia); Najam, Saher [Theoretical Plasma Physics Div., PINSTECH, P.O. Nilore, Islamabad (Pakistan); Sajid, Muhammad; Mesloub, Said [Dept. of Mathematics, Coll. of Sciences, KS Univ., Riyadh (Saudi Arabia); Ayub, Muhammad [Dept. of Mathematics, Quaid-i-Azam Univ., Islamabad (Pakistan)

2010-05-15

An analysis is performed for the slip effects on the exact solutions of flows in a generalized Burgers fluid. The flow modelling is based upon the magnetohydrodynamic (MHD) nature of the fluid and modified Darcy law in a porous space. Two illustrative examples of oscillatory flows are considered. The results obtained are compared with several limiting cases. It has been shown here that the derived results hold for all values of frequencies including the resonant frequency. (orig.)

Biorheological Model on Flow of Herschel-Bulkley Fluid through a Tapered Arterial Stenosis with Dilatation

Directory of Open Access Journals (Sweden)

S. Priyadharshini

2015-01-01

Full Text Available An analysis of blood flow through a tapered artery with stenosis and dilatation has been carried out where the blood is treated as incompressible Herschel-Bulkley fluid. A comparison between numerical values and analytical values of pressure gradient at the midpoint of stenotic region shows that the analytical expression for pressure gradient works well for the values of yield stress till 2.4. The wall shear stress and flow resistance increase significantly with axial distance and the increase is more in the case of converging tapered artery. A comparison study of velocity profiles, wall shear stress, and flow resistance for Newtonian, power law, Bingham-plastic, and Herschel-Bulkley fluids shows that the variation is greater for Herschel-Bulkley fluid than the other fluids. The obtained velocity profiles have been compared with the experimental data and it is observed that blood behaves like a Herschel-Bulkley fluid rather than power law, Bingham, and Newtonian fluids. It is observed that, in the case of a tapered stenosed tube, the streamline pattern follows a convex pattern when we move from r/R=0 to r/R=1 and it follows a concave pattern when we move from r/R=0 to r/R=-1. Further, it is of opposite behaviour in the case of a tapered dilatation tube which forms new information that is, for the first time, added to the literature.

An analytical model for the prediction of fluid-elastic forces in a rod bundle subjected to axial flow: theory, experimental validation and application to PWR fuel assemblies

International Nuclear Information System (INIS)

Beaud, F.

1997-01-01

A model predicting the fluid-elastic forces in a bundle of circular cylinders subjected to axial flow is presented in this paper. Whereas previously published models were limited to circular flow channel, the present one allows to take a rectangular flow external boundary into account. For that purpose, an original approach is derived from the standard method of images. This model will eventually be used to predict the fluid-structure coupling between the flow of primary coolant and a fuel assemblies in PWR nuclear reactors. It is indeed of major importance since the flow is shown to induce quite high damping and could therefore mitigate the incidence of an external load like a seismic excitation on the dynamics of the assemblies. The proposed model is validated on two cases from the literature but still needs further comparisons with the experiments being currently carried out on the EDF set-up. The flow has been shown to induce an approximate 12% damping on a PWR fuel assembly, at nominal reactor conditions. The possible grid effect on the fluid-structure coupling has been neglected so far but will soon be investigated at EDF. (author)

A critical flow model for the Cathena thermalhydraulic code

International Nuclear Information System (INIS)

Popov, N.K.; Hanna, B.N.

1990-01-01

The calculation of critical flow rate, e.g., of choked flow through a break, is required for simulating a loss of coolant transient in a reactor or reactor-like experimental facility. A model was developed to calculate the flow rate through the break for given geometrical parameters near the break and fluid parameters upstream of the break for ordinary water, as well as heavy water, with or without non- condensible gases. This model has been incorporated in the CATHENA, one-dimensional, two-fluid thermalhydraulic code. In the CATHENA code a standard staggered-mesh, finite-difference representation is used to solve the thermalhydraulic equations. This model compares the fluid mixture velocity, calculated using the CATHENA momentum equations, with a critical velocity. When the mixture velocity is smaller than the critical velocity, the flow is assumed to be subcritical, and the model remains passive. When the fluid mixture velocity is higher than the critical velocity, the model sets the fluid mixture velocity equal to the critical velocity. In this paper the critical velocity at a link (momentum cell) is first estimated separately for single-phase liquid, two- phase, or single-phase gas flow condition at the upstream node (mass/energy cell). In all three regimes non-condensible gas can be present in the flow. For single-phase liquid flow, the critical velocity is estimated using a Bernoulli- type of equation, the pressure at the link is estimated by the pressure undershoot method

On Cattaneo–Christov heat flux model for Carreau fluid flow over a slendering sheet

Directory of Open Access Journals (Sweden)

Hashim

Full Text Available The underlying intentions of this article are to investigate the impact of non-Fourier heat flux model on the stagnation-point flow of non-Newtonian Carreau fluid. In this study, the innovative Cattaneo–Christov constitutive model is introduced to study the characteristics of thermal relaxation time. The flow is impelled by a slendering surface which is of the variable thickness. In the model, the physical mechanism responsible for homogeneous–heterogeneous reactions are further taken into account. Also, the diffusion coefficients of the reactant and auto catalyst are considered to be equal. The governing non-linear partial differential equations consisting of the momentum, energy and concentration equations are reduced to the coupled ordinary differential equations by means of local similarity transformations. The transformed ODEs are tackled numerically by employing an effective shooting algorithm along with the Runge–Kutta Fehlberg scheme. The physical characteristics of the fluid velocity, temperature and concentration profiles are illuminated with the variation of numerous governing factors and are presented graphically. For instance, our result indicates that the temperature and thermal boundary layer thickness are lower in case of Cattaneo–Christov heat flux model when compared to classical Fourier’s heat model. Meanwhile, the rate of heat transfer is significantly improved by a high wall thickness parameter and an opposite influence is found due to the thermal relaxation parameter. We further noticed that a higher value of homogeneous and heterogeneous reaction parameter corresponds to a deceleration in the concentration field and it shows an inverse relation for the Schmidt number. A correlation with accessible results for specific cases is found with fabulous consent. Keywords: Cattaneo–Christov model, Carreau fluid, Slendering sheet, Homogeneous–heterogeneous reactions, Runge–Kutta method

Interstitial Fluid Flow: The Mechanical Environment of Cells and Foundation of Meridians

Directory of Open Access Journals (Sweden)

Wei Yao

2012-01-01

Full Text Available Using information from the deep dissection, microobservation, and measurement of acupoints in the upper and lower limbs of the human body, we developed a three-dimensional porous medium model to simulate the flow field using FLUENT software and to study the shear stress on the surface of interstitial cells (mast cells caused by interstitial fluid flow. The numerical simulation results show the following: (i the parallel nature of capillaries will lead to directional interstitial fluid flow, which may explain the long interstitial tissue channels or meridians observed in some experiments; (ii when the distribution of capillaries is staggered, increases in the velocity alternate, and the velocity tends to be uniform, which is beneficial for substance exchange; (iii interstitial fluid flow induces a shear stress, with magnitude of several Pa, on interstitial cell membranes, which will activate cells and lead to a biological response; (iv capillary and interstitial parameters, such as capillary density, blood pressure, capillary permeability, interstitial pressure, and interstitial porosity, affect the shear stress on cell surfaces. The numerical simulation results suggest that in vivo interstitial fluid flow constitutes the mechanical environment of cells and plays a key role in guiding cell activities, which may explain the meridian phenomena and the acupuncture effects observed in experiments.

Calibration of a Numerical Model for Heat Transfer and Fluid Flow in an Extruder

DEFF Research Database (Denmark)

Hofstätter, Thomas; Pedersen, David Bue; Nielsen, Jakob Skov

2016-01-01

This paper discusses experiments performed in order to validate simulations on a fused deposition modelling (FDM) extruder. The nozzle has been simulated in terms of heat transfer and fluid flow. In order to calibrate and validate these simulations, experiments were performed giving a significant...... look into the physical behaviour of the nozzle, heating and cooling systems. Experiments on the model were performed at different sub-mm diameters of the extruder. Physical parameters of the model – especially temperature dependent parameters – were set into analytical relationships in order to receive...... dynamical parameters. This research sets the foundation for further research within melted extrusion based additive manufacturing. The heating process of the extruder will be described and a note on the material feeding will be given....

A model to analyse the flow of an incompressible Newtonian fluid through a rigid, homogeneous, isotropic and infinite porous medium

International Nuclear Information System (INIS)

Gama, R.M.S. da; Sampaio, R.

1985-01-01

The flow of an incompressible Newtonian fluid through a rigid, homogeneous, isotropic and infinite porous medium which has a given inicial distribuition of the mentioned fluid, is analyzed. It is proposed a model that assumes that the motion is caused by concentration gradient, but it does not consider the friction between the porous medium and the fluid. We solve an onedimensional case where the mathematical problem is reduced to the solution of a non-linear hyperbolic system of differential equations, subjected to an inicial condition given by a step function, called 'Riemann Problem'. (Author) [pt

Regulation of tumor invasion by interstitial fluid flow

International Nuclear Information System (INIS)

Shieh, Adrian C; Swartz, Melody A

2011-01-01

The importance of the tumor microenvironment in cancer progression is undisputed, yet the significance of biophysical forces in the microenvironment remains poorly understood. Interstitial fluid flow is a nearly ubiquitous and physiologically relevant biophysical force that is elevated in tumors because of tumor-associated angiogenesis and lymphangiogenesis, as well as changes in the tumor stroma. Not only does it apply physical forces to cells directly, but interstitial flow also creates gradients of soluble signals in the tumor microenvironment, thus influencing cell behavior and modulating cell–cell interactions. In this paper, we highlight our current understanding of interstitial fluid flow in the context of the tumor, focusing on the physical changes that lead to elevated interstitial flow, how cells sense flow and how they respond to changes in interstitial flow. In particular, we emphasize that interstitial flow can directly promote tumor cell invasion through a mechanism known as autologous chemotaxis, and indirectly support tumor invasion via both biophysical and biochemical cues generated by stromal cells. Thus, interstitial fluid flow demonstrates how important biophysical factors are in cancer, both by modulating cell behavior and coupling biophysical and biochemical signals

Nonisothermal flow of a non-Newtonian fluid with viscous heating between two parallel plates

International Nuclear Information System (INIS)

Imal, M.; Pinarbasi, A.

2004-01-01

In this study the pressure gradient-flow rate relationship for steady-state nonisothermal pressure-driven flow of a non-Newtonian fluid in a channel is investigated including the effect of viscous heating is taken into account. The viscosity of the fluid depends on both temperature and shear-rate. Exponential dependence of viscosity on temperature is modelled through Arrhenius law. Non-Newtonian behaviour of the fluid is modelled according to the Carreau rheological equation, which reflects the characteristics of most polymers adequately with an exponential temperature dependence of viscosity. Flow governing motion and energy balance equations are coupled and solution of this non-linear boundary value problem is found iteratively using a pseudo spectral method based on Chebyshev polynomials. The effect of activation energy parameter and Brinkman number, as well as the power-law index and material time constant on the flow is studied. It is found that while the pressure gradient-flow rate graph is monotonic for certain ranges of flow controlling parameters, there is a large jump in the graph under certain values of these parameters.(1 table and 5 figures are included.)

Modelling of Emulsion Flow in Porous Media

Energy Technology Data Exchange (ETDEWEB)

Abou-Kassem, J.H. [UAE University (United Arab Emirates); Farouq Ali, S.M. [UAE University (United Arab Emirates)

1995-06-01

Oil recovery methods predominantly involve emulsion formation. Oil recovery simulation requires the incorporation of emulsion characteristics and flow in porous media, in order to optimize oil recovery from petroleum reservoirs. This paper explored the nature and rheology of emulsions, and evaluated several models of flow of Newtonian and non-Newtonian fluids in porous media. It also summarized in situ emulsion formation in porous media. A model for both Newtonian and non-Newtonian emulsion fluid flow was proposed, with special emphasis on pore size, and tortuosity in the porous media.

CFD simulation of gas and non-Newtonian fluid two-phase flow in anaerobic digesters.

Science.gov (United States)

Wu, Binxin

2010-07-01

This paper presents an Eulerian multiphase flow model that characterizes gas mixing in anaerobic digesters. In the model development, liquid manure is assumed to be water or a non-Newtonian fluid that is dependent on total solids (TS) concentration. To establish the appropriate models for different TS levels, twelve turbulence models are evaluated by comparing the frictional pressure drops of gas and non-Newtonian fluid two-phase flow in a horizontal pipe obtained from computational fluid dynamics (CFD) with those from a correlation analysis. The commercial CFD software, Fluent12.0, is employed to simulate the multiphase flow in the digesters. The simulation results in a small-sized digester are validated against the experimental data from literature. Comparison of two gas mixing designs in a medium-sized digester demonstrates that mixing intensity is insensitive to the TS in confined gas mixing, whereas there are significant decreases with increases of TS in unconfined gas mixing. Moreover, comparison of three mixing methods indicates that gas mixing is more efficient than mixing by pumped circulation while it is less efficient than mechanical mixing.

Double stratification effects in chemically reactive squeezed Sutterby fluid flow with thermal radiation and mixed convection

Directory of Open Access Journals (Sweden)

S. Ahmad

2018-03-01

Full Text Available A current analysis is carried out to study theoretically the mixed convection characteristics in squeezing flow of Sutterby fluid in squeezed channel. The constitutive equation of Sutterby model is utilized to characterize the rheology of squeezing phenomenon. Flow characteristics are explored with dual stratification. In flowing fluid which contains heat and mass transport, the first order chemical reaction and radiative heat flux affect the transport phenomenon. The systems of non-linear governing equations have been modulating which then solved by mean of convergent approach (Homotopy Analysis Method. The graphs are reported and illustrated for emerging parameters. Through graphical explanations, drag force, rate of heat and mass transport are conversed for different pertinent parameters. It is found that heat and mass transport rate decays with dominant double stratified parameters and chemical reaction parameter. The present two-dimensional examination is applicable in some of the engineering processes and industrial fluid mechanics. Keywords: Squeezing flow, Sutterby fluid model, Mixed convection, Double stratification, Thermal radiation, Chemical reaction

Modelling of reactive fluid transport in deformable porous rocks

Science.gov (United States)

Yarushina, V. M.; Podladchikov, Y. Y.

2009-04-01

One outstanding challenge in geology today is the formulation of an understanding of the interaction between rocks and fluids. Advances in such knowledge are important for a broad range of geologic settings including partial melting and subsequent migration and emplacement of a melt into upper levels of the crust, or fluid flow during regional metamorphism and metasomatism. Rock-fluid interaction involves heat and mass transfer, deformation, hydrodynamic flow, and chemical reactions, thereby necessitating its consideration as a complex process coupling several simultaneous mechanisms. Deformation, chemical reactions, and fluid flow are coupled processes. Each affects the others. Special effort is required for accurate modelling of the porosity field through time. Mechanical compaction of porous rocks is usually treated under isothermal or isoentropic simplifying assumptions. However, joint consideration of both mechanical compaction and reactive porosity alteration requires somewhat greater than usual care about thermodynamic consistency. Here we consider the modelling of multi-component, multi-phase systems, which is fundamental to the study of fluid-rock interaction. Based on the conservation laws for mass, momentum, and energy in the form adopted in the theory of mixtures, we derive a thermodynamically admissible closed system of equations describing the coupling of heat and mass transfer, chemical reactions, and fluid flow in a deformable solid matrix. Geological environments where reactive transport is important are located at different depths and accordingly have different rheologies. In the near surface, elastic or elastoplastic properties would dominate, whereas viscoplasticity would have a profound effect deeper in the lithosphere. Poorly understood rheologies of heterogeneous porous rocks are derived from well understood processes (i.e., elasticity, viscosity, plastic flow, fracturing, and their combinations) on the microscale by considering a

Direct numerical simulation of granular flows with fluid; Simulation numerique directe d'ecoulements granulaires en presence de fluide

Energy Technology Data Exchange (ETDEWEB)

Komiwes, V.

1999-09-01

Numerical models applied to simulation of granular flow with fluid are developed. The physical model selected to describe particles flow is a discrete approach. Particle trajectories are calculated by the Newton law and collision is describe by a soft-sphere approach. The fluid flow is modelled by Navier-Stokes equations. The modelling of the momentum transfer depends on the resolution scale: for a scale of the order of the particle diameter, it is modelled by a drag-law and for a scale smaller than the particle diameter, it is directly calculated by stress tensor computation around particles. The direct model is used to find representative elementary volume and prove the local character of the Ergun's law. This application shows the numerical (mesh size), physical (Reynolds number) and computational (CPU time and memory consumptions) limitations. The drag law model and the direct model are validated with analytical and empirical solutions and compared. For the two models, the CPU time and the memory consumptions are discussed. The drag law model is applied to the simulation of gas-solid dense fluidized-beds. In the case of uniform gas distribution, the fluidized-bed simulation heights are compared to experimental data for particle of group A and B of the Geldart classification. (author)

Focused Fluid Flow along Convergent Plate Boundaries - Deriving Flow Rates along Faults from Local Upwarping of the Base of the Gas Hydrate Stability Zone

Science.gov (United States)

Kunath, P.; Chi, W. C.; Liu, C. S.

2017-12-01

Convergent plate boundaries provide the ideal opportunity to examine the interactions of deformation, fluid flow and gas hydrate stability. However, there are still processes and parameters that remain unclear or scarce. This may be in part due to the fact that in situ determination of fluid flow rate is very difficult. Here, we present a newly developed 2-D hydrothermal model for (1) simulating the steady state, thermal effect of forced heat advection along a thin and shallow dipping fault and (2) quantifying fluid velocities required to deliver a thermal anomalies manifested at the bottom-simulating reflector (BSR) at its intersection with the fault zone. Assuming the horizontal thermal conduction is negligible, we derive our model using only a few crucial parameters: (a) the thermal conductivity structure between seafloor and fault; (b) the temperature at BSR depth and the seafloor; (c) fluid flow rate; (d) geometry of the fault conduit, including depth and thickness. Temperature disturbance is then described as a function of Peclet number and of the dip of the fault. Application of our model to Site 892 at Cascadia accretionary wedge (ODP Leg 146), where borehole data provide excellent thermal constraints on the hydrology, shows consistent results. By comparing the temperatures derived at the BSRs with the temperature field of our model, the results demonstrate that the temperature discrepancy is about 0 - 0.5 oC. We propose that this simple approach can provide, on the basis of a few parameters, rough estimate of the disturbance of the temperature caused by advecting fluid. Localized lateral BSR-based heat flow variations have been observed near thrust faults along many convergent plate boundaries around the world and are associated with strong localized fluid flow. We wish to further testing this approach using other seismic datasets to estimate first order of magnitude fluid migration patterns in other convergent boundaries.

Fluid-structure coupling between a vibrating cylinder and a narrow annular flow

International Nuclear Information System (INIS)

Perotin, L.

1994-01-01

This paper presents an analytical investigation of the fluidelastic coupling between an axial annular flow and a flexible vibrating axisymmetrical structure. The model presented is suited to single-phase, incompressible, viscous fluids and to annular flows of variable cross-section, axially symmetrical when the structure is motionless.An experimental validation of this model is presented at the end of the paper: the results obtained with the numerical model are compared with experimental data for an oscillating cylinder free to vibrate under the effect of a variable-cross-section annular flow. ((orig.))

Computational electrochemo-fluid dynamics modeling in a uranium electrowinning cell

International Nuclear Information System (INIS)

Kim, K.R.; Choi, S.Y.; Kim, S.H.; Shim, J.B.; Paek, S.; Kim, I.T.

2014-01-01

A computational electrochemo-fluid dynamics model has been developed to describe the electrowinning behavior in an electrolyte stream through a planar electrode cell system. Electrode reaction of the uranium electrowinning process from a molten-salt electrolyte stream was modeled to illustrate the details of the flow-assisted mass transport of ions to the cathode. This modeling approach makes it possible to represent variations of the convective diffusion limited current density by taking into account the concentration profile at the electrode surface as a function of the flow characteristics and applied current density in a commercially available computational fluid dynamics platform. It was possible to predict the conventional current-voltage relation in addition to details of electrolyte fluid dynamics and electrochemical variables, such as the flow field, species concentrations, potential, and current distributions throughout the galvanostatic electrolysis cell. (author)

Some numerical methods for two-fluid two-phase flows in oil pipes; Quelques methodes numeriques pour les ecoulements diphasiques bi-fluide en conduites petrolieres