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.

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)

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

Water Hammer Simulations of MMH Propellant - New Capability Demonstration of the Generalized Fluid Flow Simulation Program

Science.gov (United States)

Burkhardt, Z.; Ramachandran, N.; Majumdar, A.

2017-01-01

Fluid Transient analysis is important for the design of spacecraft propulsion system to ensure structural stability of the system in the event of sudden closing or opening of the valve. Generalized Fluid System Simulation Program (GFSSP), a general purpose flow network code developed at NASA/MSFC is capable of simulating pressure surge due to sudden opening or closing of valve when thermodynamic properties of real fluid are available for the entire range of simulation. Specifically GFSSP needs an accurate representation of pressure-density relationship in order to predict pressure surge during a fluid transient. Unfortunately, the available thermodynamic property programs such as REFPROP, GASP or GASPAK does not provide the thermodynamic properties of Monomethylhydrazine (MMH). This paper will illustrate the process used for building a customized table of properties of state variables from available properties and speed of sound that is required by GFSSP for simulation. Good agreement was found between the simulations and measured data. This method can be adopted for modeling flow networks and systems with other fluids whose properties are not known in detail in order to obtain general technical insight. Rigorous code validation of this approach will be done and reported at a future date.

Two-Fluid Models for Simulating Dispersed Multiphase Flows-A Review

Directory of Open Access Journals (Sweden)

L.X. Zhou

2009-01-01

Full Text Available The development of two-fluid models for simulating dispersed multiphase flows (gas-particle, gas-droplet, bubble-liquid, liquid-particle flows by the present author within the last 20 years is systematically reviewed. The two-fluid models based on Reynolds expansion, time averaging and mass-weighed averaging, and also PDF transport equations are described. Different versions of two-phase turbulence models, including the unified second-order moment (USM and k-ε-kp models, the DSM-PDF model, the SOM-MC model, the nonlinear k-e-kp model, and the USM-Θ model for dense gas-particle flows and their application and experimental validation are discussed.

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

Numerical simulations of heterogeneous chemical reactions coupled to fluid flow in varying thermal fields

International Nuclear Information System (INIS)

Carnahan, C.L.

1991-11-01

A numerical simulator of reactive chemical transport with coupling from precipitation-dissolution reactions to fluid flow, via changes of porosity and permeability, is applied to precipitation-dissolution of quartz and calcite in spatially and temporally variable fields of temperature. Significant effects on fluid flow are found in the quartz-silicic acid system in the presence of persistent, strong gradient of temperature. Transient heat flow in the quartz-silicic acid system and in a calcite-calcium ion-carbonato species system produces vanishingly small effects on fluid flow

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.

Effect of fluid friction on interstitial fluid flow coupled with blood flow through solid tumor microvascular network.

Science.gov (United States)

Sefidgar, Mostafa; Soltani, M; Raahemifar, Kaamran; Bazmara, Hossein

2015-01-01

A solid tumor is investigated as porous media for fluid flow simulation. Most of the studies use Darcy model for porous media. In Darcy model, the fluid friction is neglected and a few simplified assumptions are implemented. In this study, the effect of these assumptions is studied by considering Brinkman model. A multiscale mathematical method which calculates fluid flow to a solid tumor is used in this study to investigate how neglecting fluid friction affects the solid tumor simulation. The mathematical method involves processes such as blood flow through vessels and solute and fluid diffusion, convective transport in extracellular matrix, and extravasation from blood vessels. The sprouting angiogenesis model is used for generating capillary network and then fluid flow governing equations are implemented to calculate blood flow through the tumor-induced capillary network. Finally, the two models of porous media are used for modeling fluid flow in normal and tumor tissues in three different shapes of tumors. Simulations of interstitial fluid transport in a solid tumor demonstrate that the simplifications used in Darcy model affect the interstitial velocity and Brinkman model predicts a lower value for interstitial velocity than the values that Darcy model predicts.

Development of a computational framework on fluid-solid mixture flow simulations for the COMPASS code

International Nuclear Information System (INIS)

Zhang, Shuai; Morita, Koji; Shirakawa, Noriyuki; Yamamoto, Yuichi

2010-01-01

The COMPASS code is designed based on the moving particle semi-implicit method to simulate various complex mesoscale phenomena relevant to core disruptive accidents of sodium-cooled fast reactors. In this study, a computational framework for fluid-solid mixture flow simulations was developed for the COMPASS code. The passively moving solid model was used to simulate hydrodynamic interactions between fluid and solids. Mechanical interactions between solids were modeled by the distinct element method. A multi-time-step algorithm was introduced to couple these two calculations. The proposed computational framework for fluid-solid mixture flow simulations was verified by the comparison between experimental and numerical studies on the water-dam break with multiple solid rods. (author)

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)

Numerical Simulation of Complex Multi-Fluid Flows using a Combined Immersed Boundary and Volume of Fluid Approach

NARCIS (Netherlands)

Deen, N.G.; van Sint Annaland, M.; Kuipers, J.A.M.

2007-01-01

In this paper a simulation model is presented for the Direct Numerical Simulation (DNS) of complex multi-fluid flows in which simultaneously (moving) deformable (drops or bubbles) and non-deformable (moving) elements (particles) are present, possibly with the additional presence of free surfaces.

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.

Hydrodynamically Coupled Brownian Dynamics simulations for flow on non-Newtonian fluids

NARCIS (Netherlands)

Ahuja, Vishal Raju

2018-01-01

This thesis deals with model development for particle-based flow simulations of non-Newtonian fluids such as polymer solutions. A novel computational technique called Hydrodynamically Coupled Brownian Dynamics (HCBD) is presented in this thesis. This technique essentially couples the Brownian motion

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.)

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.

Numerical simulation of fluid flow in microporous media

International Nuclear Information System (INIS)

Xu Ruina; Jiang Peixue

2008-01-01

The flow characteristics of water and air in microporous media with average diameters of 200 μm, 125 μm, 90 μm, 40 μm, 20 μm, and 10 μm were studied numerically. The calculated friction factors for water and air in the non-slip-flow regime in the microporous media agree well with the known correlation suitable for normal size porous media. The numerically predicted friction factors for air in the slip-flow regime in the microporous media with 90 μm, 40 μm, 20 μm, and 10 μm diameter particles were less than the correlation for normal size porous media but close to experimental data and a modified correlation that accounts for rarefaction. Comparisons of the numerical results with the experimental data and the modified correlations show that rarefaction effects occur in air flows in the microporous media with particle diameters less than 90 μm and that the numerical calculations with velocity slip on the boundary can properly simulate the fluid flow in microporous media

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...

3D code for simulations of fluid flows

International Nuclear Information System (INIS)

Skandera, D.

2004-01-01

In this paper, a present status in the development of the new numerical code is reported. The code is considered for simulations of fluid flows. The finite volume approach is adopted for solving standard fluid equations. They are treated in a conservative form to ensure a correct conservation of fluid quantities. Thus, a nonlinear hyperbolic system of conservation laws is numerically solved. The code uses the Eulerian description of the fluid and is designed as a high order central numerical scheme. The central approach employs no (approximate) Riemann solver and is less computational expensive. The high order WENO strategy is adopted in the reconstruction step to achieve results comparable with more accurate Riemann solvers. A combination of the central approach with an iterative solving of a local Riemann problem is tested and behaviour of such numerical flux is reported. An extension to three dimensions is implemented using a dimension by dimension approach, hence, no complicated dimensional splitting need to be introduced. The code is fully parallelized with the MPI library. Several standard hydrodynamic tests in one, two and three dimensions were performed and their results are presented. (author)

Virtual reality aided visualization of fluid flow simulations with application in medical education and diagnostics.

Science.gov (United States)

Djukic, Tijana; Mandic, Vesna; Filipovic, Nenad

2013-12-01

Medical education, training and preoperative diagnostics can be drastically improved with advanced technologies, such as virtual reality. The method proposed in this paper enables medical doctors and students to visualize and manipulate three-dimensional models created from CT or MRI scans, and also to analyze the results of fluid flow simulations. Simulation of fluid flow using the finite element method is performed, in order to compute the shear stress on the artery walls. The simulation of motion through the artery is also enabled. The virtual reality system proposed here could shorten the length of training programs and make the education process more effective. © 2013 Published by Elsevier Ltd.

Parallel Simulation of Three-Dimensional Free Surface Fluid Flow Problems

International Nuclear Information System (INIS)

BAER, THOMAS A.; SACKINGER, PHILIP A.; SUBIA, SAMUEL R.

1999-01-01

Simulation of viscous three-dimensional fluid flow typically involves a large number of unknowns. When free surfaces are included, the number of unknowns increases dramatically. Consequently, this class of problem is an obvious application of parallel high performance computing. We describe parallel computation of viscous, incompressible, free surface, Newtonian fluid flow problems that include dynamic contact fines. The Galerkin finite element method was used to discretize the fully-coupled governing conservation equations and a ''pseudo-solid'' mesh mapping approach was used to determine the shape of the free surface. In this approach, the finite element mesh is allowed to deform to satisfy quasi-static solid mechanics equations subject to geometric or kinematic constraints on the boundaries. As a result, nodal displacements must be included in the set of unknowns. Other issues discussed are the proper constraints appearing along the dynamic contact line in three dimensions. Issues affecting efficient parallel simulations include problem decomposition to equally distribute computational work among a SPMD computer and determination of robust, scalable preconditioners for the distributed matrix systems that must be solved. Solution continuation strategies important for serial simulations have an enhanced relevance in a parallel coquting environment due to the difficulty of solving large scale systems. Parallel computations will be demonstrated on an example taken from the coating flow industry: flow in the vicinity of a slot coater edge. This is a three dimensional free surface problem possessing a contact line that advances at the web speed in one region but transitions to static behavior in another region. As such, a significant fraction of the computational time is devoted to processing boundary data. Discussion focuses on parallel speed ups for fixed problem size, a class of problems of immediate practical importance

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.

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.

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.

Simulating single-phase and two-phase non-Newtonian fluid flow of a digital rock scanned at high resolution

Science.gov (United States)

Tembely, Moussa; Alsumaiti, Ali M.; Jouini, Mohamed S.; Rahimov, Khurshed; Dolatabadi, Ali

2017-11-01

Most of the digital rock physics (DRP) simulations focus on Newtonian fluids and overlook the detailed description of rock-fluid interaction. A better understanding of multiphase non-Newtonian fluid flow at pore-scale is crucial for optimizing enhanced oil recovery (EOR). The Darcy scale properties of reservoir rocks such as the capillary pressure curves and the relative permeability are controlled by the pore-scale behavior of the multiphase flow. In the present work, a volume of fluid (VOF) method coupled with an adaptive meshing technique is used to perform the pore-scale simulation on a 3D X-ray micro-tomography (CT) images of rock samples. The numerical model is based on the resolution of the Navier-Stokes equations along with a phase fraction equation incorporating the dynamics contact model. The simulations of a single phase flow for the absolute permeability showed a good agreement with the literature benchmark. Subsequently, the code is used to simulate a two-phase flow consisting of a polymer solution, displaying a shear-thinning power law viscosity. The simulations enable to access the impact of the consistency factor (K), the behavior index (n), along with the two contact angles (advancing and receding) on the relative permeability.

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...

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

Simulations of fluid flow through porous media based on cellular automata and non-linear dynamics

Energy Technology Data Exchange (ETDEWEB)

Paulson, K V

1992-05-15

A study is being carried out to apply cellular automata and non-linear dynamics in the construction of efficient and accurate computer simulations of multiphase fluid flow through porous media, with the objective of application to reservoir modelling for hydrocarbon recovery. An algorithm based on Boolean operations has been developed which transforms a PC clone into a highly efficient vector processor capable of cellular automata simulation of single fluid flow through two-dimensional rock matrix models of varying porosities. Macroscopic flow patterns have been established through spatial and temporal averaging with no floating point operations. Permeabilities of the different models have been calculated. Hardware allows the algorithm to function on dual processors on a PC platform using a video recording and editing facility. Very encouraging results have been obtained. 4 figs.

The simulation of multidimensional multiphase flows

International Nuclear Information System (INIS)

Lahey, Richard T.

2005-01-01

This paper presents an assessment of various models which can be used for the multidimensional simulation of multiphase flows, such as may occur in nuclear reactors. In particular, a model appropriate for the direct numerical simulation (DNS) of multiphase flows and a mechanistically based, three-dimensional, four-field, turbulent, two-fluid computational multiphase fluid dynamics (CMFD) model are discussed. A two-fluid bubbly flow model, which was derived using potential flow theory, can be extended to other flow regimes, but this will normally involve ensemble-averaging the results from direct numerical simulations (DNS) of various flow regimes to provide the detailed numerical data necessary for the development of flow-regime-specific interfacial and wall closure laws

Engineering Fracking Fluids with Computer Simulation

Science.gov (United States)

Shaqfeh, Eric

2015-11-01

There are no comprehensive simulation-based tools for engineering the flows of viscoelastic fluid-particle suspensions in fully three-dimensional geometries. On the other hand, the need for such a tool in engineering applications is immense. Suspensions of rigid particles in viscoelastic fluids play key roles in many energy applications. For example, in oil drilling the ``drilling mud'' is a very viscous, viscoelastic fluid designed to shear-thin during drilling, but thicken at stoppage so that the ``cuttings'' can remain suspended. In a related application known as hydraulic fracturing suspensions of solids called ``proppant'' are used to prop open the fracture by pumping them into the well. It is well-known that particle flow and settling in a viscoelastic fluid can be quite different from that which is observed in Newtonian fluids. First, it is now well known that the ``fluid particle split'' at bifurcation cracks is controlled by fluid rheology in a manner that is not understood. Second, in Newtonian fluids, the presence of an imposed shear flow in the direction perpendicular to gravity (which we term a cross or orthogonal shear flow) has no effect on the settling of a spherical particle in Stokes flow (i.e. at vanishingly small Reynolds number). By contrast, in a non-Newtonian liquid, the complex rheological properties induce a nonlinear coupling between the sedimentation and shear flow. Recent experimental data have shown both the shear thinning and the elasticity of the suspending polymeric solutions significantly affects the fluid-particle split at bifurcations, as well as the settling rate of the solids. In the present work, we use the Immersed Boundary Method to develop computer simulations of viscoelastic flow in suspensions of spheres to study these problems. These simulations allow us to understand the detailed physical mechanisms for the remarkable physical behavior seen in practice, and actually suggest design rules for creating new fluid recipes.

Network Flow Simulation of Fluid Transients in Rocket Propulsion Systems

Science.gov (United States)

Bandyopadhyay, Alak; Hamill, Brian; Ramachandran, Narayanan; Majumdar, Alok

2011-01-01

Fluid transients, also known as water hammer, can have a significant impact on the design and operation of both spacecraft and launch vehicle propulsion systems. These transients often occur at system activation and shutdown. The pressure rise due to sudden opening and closing of valves of propulsion feed lines can cause serious damage during activation and shutdown of propulsion systems. During activation (valve opening) and shutdown (valve closing), pressure surges must be predicted accurately to ensure structural integrity of the propulsion system fluid network. In the current work, a network flow simulation software (Generalized Fluid System Simulation Program) based on Finite Volume Method has been used to predict the pressure surges in the feed line due to both valve closing and valve opening using two separate geometrical configurations. The valve opening pressure surge results are compared with experimental data available in the literature and the numerical results compared very well within reasonable accuracy (< 5%) for a wide range of inlet-to-initial pressure ratios. A Fast Fourier Transform is preformed on the pressure oscillations to predict the various modal frequencies of the pressure wave. The shutdown problem, i.e. valve closing problem, the simulation results are compared with the results of Method of Characteristics. Most rocket engines experience a longitudinal acceleration, known as "pogo" during the later stage of engine burn. In the shutdown example problem, an accumulator has been used in the feed system to demonstrate the "pogo" mitigation effects in the feed system of propellant. The simulation results using GFSSP compared very well with the results of Method of Characteristics.

Verification of Eulerian-Eulerian and Eulerian-Lagrangian simulations for fluid-particle flows

Science.gov (United States)

Kong, Bo; Patel, Ravi G.; Capecelatro, Jesse; Desjardins, Olivier; Fox, Rodney O.

2017-11-01

In this work, we study the performance of three simulation techniques for fluid-particle flows: (1) a volume-filtered Euler-Lagrange approach (EL), (2) a quadrature-based moment method using the anisotropic Gaussian closure (AG), and (3) a traditional two-fluid model. By simulating two problems: particles in frozen homogeneous isotropic turbulence (HIT), and cluster-induced turbulence (CIT), the convergence of the methods under grid refinement is found to depend on the simulation method and the specific problem, with CIT simulations facing fewer difficulties than HIT. Although EL converges under refinement for both HIT and CIT, its statistical results exhibit dependence on the techniques used to extract statistics for the particle phase. For HIT, converging both EE methods (TFM and AG) poses challenges, while for CIT, AG and EL produce similar results. Overall, all three methods face challenges when trying to extract converged, parameter-independent statistics due to the presence of shocks in the particle phase. National Science Foundation and National Energy Technology Laboratory.

A ghost fluid method for sharp interface simulations of compressible multiphase flows

International Nuclear Information System (INIS)

Majidi, Sahand; Afshari, Asghar

2016-01-01

A ghost fluid based computational tool is developed to study a wide range of compressible multiphase flows involving strong shocks and contact discontinuities while accounting for surface tension, viscous stresses and gravitational forces. The solver utilizes constrained reinitialization method to predict the interface configuration at each time step. Surface tension effect is handled via an exact interface Riemann problem solver. Interfacial viscous stresses are approximated by considering continuous velocity and viscous stress across the interface. To assess the performance of the solver several benchmark problems are considered: One-dimensional gas-water shock tube problem, shock-bubble interaction, air cavity collapse in water, underwater explosion, Rayleigh-Taylor Instability, and ellipsoidal drop oscillations. Results obtained from the numerical simulations indicate that the numerical methodology performs reasonably well in predicting flow features and exhibit a very good agreement with prior experimental and numerical observations. To further examine the accuracy of the developed ghost fluid solver, the obtained results are compared to those by a conventional diffuse interface solver. The comparison shows the capability of our ghost fluid method in reproducing the experimentally observed flow characteristics while revealing more details regarding topological changes of the interface.

A ghost fluid method for sharp interface simulations of compressible multiphase flows

Energy Technology Data Exchange (ETDEWEB)

Majidi, Sahand; Afshari, Asghar [University of Tehran, Teheran (Iran, Islamic Republic of)

2016-04-15

A ghost fluid based computational tool is developed to study a wide range of compressible multiphase flows involving strong shocks and contact discontinuities while accounting for surface tension, viscous stresses and gravitational forces. The solver utilizes constrained reinitialization method to predict the interface configuration at each time step. Surface tension effect is handled via an exact interface Riemann problem solver. Interfacial viscous stresses are approximated by considering continuous velocity and viscous stress across the interface. To assess the performance of the solver several benchmark problems are considered: One-dimensional gas-water shock tube problem, shock-bubble interaction, air cavity collapse in water, underwater explosion, Rayleigh-Taylor Instability, and ellipsoidal drop oscillations. Results obtained from the numerical simulations indicate that the numerical methodology performs reasonably well in predicting flow features and exhibit a very good agreement with prior experimental and numerical observations. To further examine the accuracy of the developed ghost fluid solver, the obtained results are compared to those by a conventional diffuse interface solver. The comparison shows the capability of our ghost fluid method in reproducing the experimentally observed flow characteristics while revealing more details regarding topological changes of the interface.

Buoyant miscible displacement flow of shear-thinning fluids: Experiments and Simulations

Science.gov (United States)

Ale Etrati Khosroshahi, Seyed Ali; Frigaard, Ian

2017-11-01

We study displacement flow of two miscible fluids with density and viscosity contrast in an inclined pipe. Our focus is mainly on displacements where transverse mixing is not significant and thus a two-layer, stratified flow develops. Our experiments are carried out in a long pipe, covering a wide range of flow-rates, inclination angles and viscosity ratios. Density and viscosity contrasts are achieved by adding Glycerol and Xanthan gum to water, respectively. At each angle, flow rate and viscosity ratio are varied and density contrast is fixed. We identify and map different flow regimes, instabilities and front dynamics based on Fr , Re / Frcosβ and viscosity ratio m. The problem is also studied numerically to get a better insight into the flow structure and shear-thinning effects. Numerical simulations are completed using OpenFOAM in both pipe and channel geometries and are compared against the experiments. Schlumberger, NSERC.

Direct numerical simulations of fluid flow, heat transfer and phase changes

Science.gov (United States)

Juric, D.; Tryggvason, G.; Han, J.

1997-01-01

Direct numerical simulations of fluid flow, heat transfer, and phase changes are presented. The simulations are made possible by a recently developed finite difference/front tracking method based on the one-field formulation of the governing equations where a single set of conservation equations is written for all the phases involved. The conservation equations are solved on a fixed rectangular grid, but the phase boundaries are kept sharp by tracking them explicitly by a moving grid of lower dimension. The method is discussed and applications to boiling heat transfer and the solidification of drops colliding with a wall are shown.

Fluid simulation for computer graphics

CERN Document Server

Bridson, Robert

2008-01-01

Animating fluids like water, smoke, and fire using physics-based simulation is increasingly important in visual effects, in particular in movies, like The Day After Tomorrow, and in computer games. This book provides a practical introduction to fluid simulation for graphics. The focus is on animating fully three-dimensional incompressible flow, from understanding the math and the algorithms to the actual implementation.

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.

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.

Fluid Flow Simulation and Energetic Analysis of Anomalocarididae Locomotion

Science.gov (United States)

Mikel-Stites, Maxwell; Staples, Anne

2014-11-01

While an abundance of animal locomotion simulations have been performed modeling the motions of living arthropods and aquatic animals, little quantitative simulation and reconstruction of gait parameters has been done to model the locomotion of extinct animals, many of which bear little physical resemblance to their modern descendants. To that end, this project seeks to analyze potential swimming patterns used by the anomalocaridid family, (specifically Anomalocaris canadensis, a Cambrian Era aquatic predator), and determine the most probable modes of movement. This will serve to either verify or cast into question the current assumed movement patterns and properties of these animals and create a bridge between similar flexible-bodied swimmers and their robotic counterparts. This will be accomplished by particle-based fluid flow simulations of the flow around the fins of the animal, as well as an energy analysis of a variety of sample gaits. The energy analysis will then be compared to the extant information regarding speed/energy use curves in an attempt to determine which modes of swimming were most energy efficient for a given range of speeds. These results will provide a better understanding of how these long-extinct animals moved, possibly allowing an improved understanding of their behavioral patterns, and may also lead to a novel potential platform for bio-inspired underwater autonomous vehicles (UAVs).

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.

Using artificial intelligence to control fluid flow computations

Science.gov (United States)

Gelsey, Andrew

1992-01-01

Computational simulation is an essential tool for the prediction of fluid flow. Many powerful simulation programs exist today. However, using these programs to reliably analyze fluid flow and other physical situations requires considerable human effort and expertise to set up a simulation, determine whether the output makes sense, and repeatedly run the simulation with different inputs until a satisfactory result is achieved. Automating this process is not only of considerable practical importance but will also significantly advance basic artificial intelligence (AI) research in reasoning about the physical world.

Fluid, solid and fluid-structure interaction simulations on patient-based abdominal aortic aneurysm models.

Science.gov (United States)

Kelly, Sinead; O'Rourke, Malachy

2012-04-01

This article describes the use of fluid, solid and fluid-structure interaction simulations on three patient-based abdominal aortic aneurysm geometries. All simulations were carried out using OpenFOAM, which uses the finite volume method to solve both fluid and solid equations. Initially a fluid-only simulation was carried out on a single patient-based geometry and results from this simulation were compared with experimental results. There was good qualitative and quantitative agreement between the experimental and numerical results, suggesting that OpenFOAM is capable of predicting the main features of unsteady flow through a complex patient-based abdominal aortic aneurysm geometry. The intraluminal thrombus and arterial wall were then included, and solid stress and fluid-structure interaction simulations were performed on this, and two other patient-based abdominal aortic aneurysm geometries. It was found that the solid stress simulations resulted in an under-estimation of the maximum stress by up to 5.9% when compared with the fluid-structure interaction simulations. In the fluid-structure interaction simulations, flow induced pressure within the aneurysm was found to be up to 4.8% higher than the value of peak systolic pressure imposed in the solid stress simulations, which is likely to be the cause of the variation in the stress results. In comparing the results from the initial fluid-only simulation with results from the fluid-structure interaction simulation on the same patient, it was found that wall shear stress values varied by up to 35% between the two simulation methods. It was concluded that solid stress simulations are adequate to predict the maximum stress in an aneurysm wall, while fluid-structure interaction simulations should be performed if accurate prediction of the fluid wall shear stress is necessary. Therefore, the decision to perform fluid-structure interaction simulations should be based on the particular variables of interest in a given

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)

ECMOR 4. 4th European conference on the mathematics of oil recovery. Topic D: Simulation of fluid flow. Proceedings

Energy Technology Data Exchange (ETDEWEB)

1994-01-01

The report with collected proceedings from a conference, deals with mathematics of oil recovery with the focus on simulation of fluid flow. Topics of proceedings are as follow: Validity of macroscopic viscous fingering models for 2D and 3D-flows; pressure equation for fluid flow in a stochastic medium; predicting multicomponent, multiphase flow in heterogeneous systems using streamtubes; analytic techniques in pressure transient testing; global triangular structure in four-component conservation laws; exact solution of the problem on hydrodynamic interaction between noncommunicating layers under conditions of their joint development; fluid rate in flowing granular medium with moving boundary; complex variable boundary element method for tracking streamlines across fractures; transport equations for miscible displacements in heterogeneous porous media - a streamtube approach; mathematical modelling of condensate film flow by gravity drainage; effect of capillary forces on immiscible two-phase flow in strongly heterogeneous porous media; multidomain direct method and local time steps in reservoir simulation; adaptive methods for chemical flooding; flux continuous for the full tensor equation; discretization on non-orthogonal, curvilinear grids for multi-phase flow; blending finite elements and finite volumes for the solution of miscible incompressible flow. 16 papers are prepared. 240 refs., 122 figs., 6 tabs.

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

KAUST Repository

Calo, Victor M.

2015-04-27

Pleated filters are broadly used for various applications. In certain cases, especially in solid-liquid separation case, the filtering media may get deflected and that may change the overall performance characteristics of the filter. From the modeling point of view, this is a challenging multiphysics problem, namely the interaction of the fluid with a so-called poroelastic structure. This work focuses on the development of an algorithm for the simulation of the Fluid Porous Structure Interaction (FPSI) problem in the case of pleated filtering media. The first part of the work is concerned with the development of a robust and accurate numerical method for solving the Stokes-Brinkman system of equations on quadrilateral grids. The mathematical 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. The Stokes-Brinkman system of equations is discretized here using a sophisticated finite volume method, namely multi-point flux approximation (MPFA) O-method. MPFA is widely used, e.g., in solving scalar elliptic equations with full tensor and highly varying coefficients and/or solving on heterogeneous non-orthogonalgrids. Up to the authors’ knowledge, there was no investigation of MPFA discretization for Stokes-Brinkman problems, and this study aims to fill this gap. Some numerical experiments are presented in order to demonstrate the robustness of the proposed numerical algorithm[1]. The second part of this study focuses on the coupling of the flow model with the deflection of the filtering media. For the consideration of the FPSI problem in 3D, the classical Biot system describes coupled flow and deformations in a porous body due to difference in the upstream and downstream pressures. Solving the Biot system of equations is complicated and requires a significant amount of

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.

Molecular Dynamics Simulation of Binary Fluid in a Nanochannel

International Nuclear Information System (INIS)

Mullick, Shanta; Ahluwalia, P. K.; Pathania, Y.

2011-01-01

This paper presents the results from a molecular dynamics simulation of binary fluid (mixture of argon and krypton) in the nanochannel flow. The computational software LAMMPS is used for carrying out the molecular dynamics simulations. Binary fluids of argon and krypton with varying concentration of atom species were taken for two densities 0.65 and 0.45. The fluid flow takes place between two parallel plates and is bounded by horizontal walls in one direction and periodic boundary conditions are imposed in the other two directions. To drive the flow, a constant force is applied in one direction. Each fluid atom interacts with other fluid atoms and wall atoms through Week-Chandler-Anderson (WCA) potential. The velocity profile has been looked at for three nanochannel widths i.e for 12σ, 14σ and 16σ and also for the different concentration of two species. The velocity profile of the binary fluid predicted by the simulations agrees with the quadratic shape of the analytical solution of a Poiseuille flow in continuum theory.

Direct numerical simulation of complex multi-fluid flows using a combined front tracking and immersed boundary method

NARCIS (Netherlands)

Deen, N.G.; van Sint Annaland, M.; Kuipers, J.A.M.

2009-01-01

In this paper a simulation model is presented for the Direct Numerical Simulation (DNS) of complex multi-fluid flows in which simultaneously (moving) deformable (drops or bubbles) and non-deformable (moving) elements (particles) are present, possibly with the additional presence of free surfaces.

Direct numerical simulation of solidification microstructures affected by fluid flow

International Nuclear Information System (INIS)

Juric, D.

1997-12-01

The effects of fluid flow on the solidification morphology of pure materials and solute microsegregation patterns of binary alloys are studied using a computational methodology based on a front tracking/finite difference method. A general single field formulation is presented for the full coupling of phase change, fluid flow, heat and solute transport. This formulation accounts for interfacial rejection/absorption of latent heat and solute, interfacial anisotropies, discontinuities in material properties between the liquid and solid phases, shrinkage/expansion upon solidification and motion and deformation of the solid. Numerical results are presented for the two dimensional dendritic solidification of pure succinonitrile and the solidification of globulitic grains of a plutonium-gallium alloy. For both problems, comparisons are made between solidification without fluid flow and solidification within a shear flow

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

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

Simulation of sound waves using the Lattice Boltzmann Method for fluid flow: Benchmark cases for outdoor sound propagation

NARCIS (Netherlands)

Salomons, E.M.; Lohman, W.J.A.; Zhou, H.

2016-01-01

Propagation of sound waves in air can be considered as a special case of fluid dynamics. Consequently, the lattice Boltzmann method (LBM) for fluid flow can be used for simulating sound propagation. In this article application of the LBM to sound propagation is illustrated for various cases:

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.

Prediction of Separation Length of Turbulent Multiphase Flow Using Radiotracer and Computational Fluid Dynamics Simulation

International Nuclear Information System (INIS)

Sugiharto, S.; Kurniadi, R.; Abidin, Z.; Stegowski, Z.; Furman, L.

2013-01-01

Multiphase flow modeling presents great challenges due to its extreme importance in various industrial and environmental applications. In the present study, prediction of separation length of multiphase flow is examined experimentally by injection of two kinds of iodine-based radiotracer solutions into a hydrocarbon transport pipeline (HCT) having an inner diameter of 24 in (60,96 m). The main components of fluids in the pipeline are water 95%, crude oil 3% and gas 2%. A radiotracing experiment was carried out at the segment of pipe which is located far from branch points with assumptions that stratified flows in such segment were achieved. Two radiation detectors located at 80 and 100 m from injection point were used to generate residence time distribution (RTD) curve resulting from injection of radiotracer solutions. Multiphase computational fluid dynamics (CFD) simulations using Eulerian-Eulerian control volume and commercial CFD package Fluent 6.2 were employed to simulate separation length of multiphase flow. The results of study shows that the flow velocity of water is higher than the flow rate of crude oil in water-dominated system despite the higher density of water than the density of the crude oil. The separation length in multiphase flow predicted by Fluent mixture model is approximately 20 m, measured from injection point. This result confirms that the placement of the first radiation detector at the distance 80 m from the injection point was correct. (author)

Prediction of Separation Length of Turbulent Multiphase Flow Using Radiotracer and Computational Fluid Dynamics Simulation

Directory of Open Access Journals (Sweden)

S. Sugiharto1

2013-04-01

Full Text Available Multiphase flow modeling presents great challenges due to its extreme importance in various industrial and environmental applications. In the present study, prediction of separation length of multiphase flow is examined experimentally by injection of two kinds of iodine-based radiotracer solutions into a hydrocarbon transport pipeline (HCT having an inner diameter of 24 in (60,96 m. The main components of fluids in the pipeline are water 95%, crude oil 3% and gas 2%. A radiotracing experiment was carried out at the segment of pipe which is located far from branch points with assumptions that stratified flows in such segment were achieved. Two radiation detectors located at 80 and 100 m from injection point were used to generate residence time distribution (RTD curve resulting from injection of radiotracer solutions. Multiphase computational fluid dynamics (CFD simulations using Eulerian-Eulerian control volume and commercial CFD package Fluent 6.2 were employed to simulate separation length of multiphase flow. The results of study shows that the flow velocity of water is higher than the flow rate of crude oil in water-dominated system despite the higher density of water than the density of the crude oil. The separation length in multiphase flow predicted by Fluent mixture model is approximately 20 m, measured from injection point. This result confirms that the placement of the first radiation detector at the distance 80 m from the injection point was correct

Shaded computer graphic techniques for visualizing and interpreting analytic fluid flow models

Science.gov (United States)

Parke, F. I.

1981-01-01

Mathematical models which predict the behavior of fluid flow in different experiments are simulated using digital computers. The simulations predict values of parameters of the fluid flow (pressure, temperature and velocity vector) at many points in the fluid. Visualization of the spatial variation in the value of these parameters is important to comprehend and check the data generated, to identify the regions of interest in the flow, and for effectively communicating information about the flow to others. The state of the art imaging techniques developed in the field of three dimensional shaded computer graphics is applied to visualization of fluid flow. Use of an imaging technique known as 'SCAN' for visualizing fluid flow, is studied and the results are presented.

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

Computational Fluid Dynamics simulation of hydrothermal liquefaction of microalgae in a continuous plug-flow reactor.

Science.gov (United States)

Ranganathan, Panneerselvam; Savithri, Sivaraman

2018-06-01

Computational Fluid Dynamics (CFD) technique is used in this work to simulate the hydrothermal liquefaction of Nannochloropsis sp. microalgae in a lab-scale continuous plug-flow reactor to understand the fluid dynamics, heat transfer, and reaction kinetics in a HTL reactor under hydrothermal condition. The temperature profile in the reactor and the yield of HTL products from the present simulation are obtained and they are validated with the experimental data available in the literature. Furthermore, the parametric study is carried out to study the effect of slurry flow rate, reactor temperature, and external heat transfer coefficient on the yield of products. Though the model predictions are satisfactory in comparison with the experimental results, it still needs to be improved for better prediction of the product yields. This improved model will be considered as a baseline for design and scale-up of large-scale HTL reactor. Copyright © 2018 Elsevier Ltd. All rights reserved.

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.

Application of computational fluid dynamics to closed-loop bioreactors: I. Characterization and simulation of fluid-flow pattern and oxygen transfer.

Science.gov (United States)

Littleton, Helen X; Daigger, Glen T; Strom, Peter F

2007-06-01

A full-scale, closed-loop bioreactor (Orbal oxidation ditch, Envirex brand technologies, Siemens, Waukesha, Wisconsin), previously examined for simultaneous biological nutrient removal (SBNR), was further evaluated using computational fluid dynamics (CFD). A CFD model was developed first by imparting the known momentum (calculated by tank fluid velocity and mass flowrate) to the fluid at the aeration disc region. Oxygen source (aeration) and sink (consumption) terms were introduced, and statistical analysis was applied to the CFD simulation results. The CFD model was validated with field data obtained from a test tank and a full-scale tank. The results indicated that CFD could predict the mixing pattern in closed-loop bioreactors. This enables visualization of the flow pattern, both with regard to flow velocity and dissolved-oxygen-distribution profiles. The velocity and oxygen-distribution gradients suggested that the flow patterns produced by directional aeration in closed-loop bioreactors created a heterogeneous environment that can result in dissolved oxygen variations throughout the bioreactor. Distinct anaerobic zones on a macroenvironment scale were not observed, but it is clear that, when flow passed around curves, a secondary spiral flow was generated. This second current, along with the main recirculation flow, could create alternating anaerobic and aerobic conditions vertically and horizontally, which would allow SBNR to occur. Reliable SBNR performance in Orbal oxidation ditches may be a result, at least in part, of such a spatially varying environment.

High definition graphics application in fluid flow simulations

International Nuclear Information System (INIS)

Bancroft, G.; Merritt, F.; Buning, P.; Watson, V.

1987-01-01

NASA Ames is using high spatial/color-resolution computer graphics to interactively visualize flow fields generated by supercomputer solutions for the flow about such vehicles as the Space Shuttle. The scientist's viewing position in three-dimensional space can be interactively changed while the fluid flow is either frozen or moving in time. Animated sequences can then be made for workstation viewing with the aid of specialized software that allows easy editing and automatic tweening of the sequences. Attention is presently given to the software generating the three-dimensional flow field displays, as well as that for creating the animation sequences

Digital Rock Simulation of Flow in Carbonate Samples

Science.gov (United States)

Klemin, D.; Andersen, M.

2014-12-01

Reservoir engineering has becomes more complex to deal with current challenges, so core analysts must understand and model pore geometries and fluid behaviors at pores scales more rapidly and realistically. We introduce an industry-unique direct hydrodynamic pore flow simulator that operates on pore geometries from digital rock models obtained using microCT or 3D scanning electron microscope (SEM) images. The PVT and rheological models used in the simulator represent real reservoir fluids. Fluid-solid interactions are introduced using distributed micro-scale wetting properties. The simulator uses density functional approach applied for hydrodynamics of complex systems. This talk covers selected applications of the simulator. We performed microCT scanning of six different carbonate rock samples from homogeneous limestones to vuggy carbonates. From these, we constructed digital rock models representing pore geometries for the simulator. We simulated nonreactive tracer flow in all six digital models using a digital fluid description that included a passive tracer solution. During the simulation, we evaluated the composition of the effluent. Results of tracer flow simulations corresponded well with experimental data of nonreactive tracer floods for the same carbonate rock types. This simulation data of the non-reactive tracer flow can be used to calculate the volume of the rock accessible by the fluid, which can be further used to predict response of a porous medium to a reactive fluid. The described digital core analysis workflow provides a basis for a wide variety of activities, including input to design acidizing jobs and evaluating treatment efficiency and EOR economics. Digital rock multiphase flow simulations of a scanned carbonate rock evaluated the effect of wettability on flow properties. Various wetting properties were tested: slightly oil wet, slightly water wet, and water wet. Steady-state relative permeability simulations yielded curves for all three

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.

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.

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

Generalized Fluid System Simulation Program (GFSSP) Version 6 - General Purpose Thermo-Fluid Network Analysis Software

Science.gov (United States)

Majumdar, Alok; Leclair, Andre; Moore, Ric; Schallhorn, Paul

2011-01-01

GFSSP stands for Generalized Fluid System Simulation Program. It is a general-purpose computer program to compute pressure, temperature and flow distribution in a flow network. GFSSP calculates pressure, temperature, and concentrations at nodes and calculates flow rates through branches. It was primarily developed to analyze Internal Flow Analysis of a Turbopump Transient Flow Analysis of a Propulsion System. GFSSP development started in 1994 with an objective to provide a generalized and easy to use flow analysis tool for thermo-fluid systems.

From capture to simulation: connecting forward and inverse problems in fluids

KAUST Repository

Gregson, James; Ihrke, Ivo; Thuerey, Nils; Heidrich, Wolfgang

2014-01-01

We explore the connection between fluid capture, simulation and proximal methods, a class of algorithms commonly used for inverse problems in image processing and computer vision. Our key finding is that the proximal operator constraining fluid velocities to be divergence-free is directly equivalent to the pressure-projection methods commonly used in incompressible flow solvers. This observation lets us treat the inverse problem of fluid tracking as a constrained flow problem all while working in an efficient, modular framework. In addition it lets us tightly couple fluid simulation into flow tracking, providing a global prior that significantly increases tracking accuracy and temporal coherence as compared to previous techniques. We demonstrate how we can use these improved results for a variety of applications, such as re-simulation, detail enhancement, and domain modification. We furthermore give an outlook of the applications beyond fluid tracking that our proximal operator framework could enable by exploring the connection of deblurring and fluid guiding.

From capture to simulation: connecting forward and inverse problems in fluids

KAUST Repository

Gregson, James

2014-07-27

We explore the connection between fluid capture, simulation and proximal methods, a class of algorithms commonly used for inverse problems in image processing and computer vision. Our key finding is that the proximal operator constraining fluid velocities to be divergence-free is directly equivalent to the pressure-projection methods commonly used in incompressible flow solvers. This observation lets us treat the inverse problem of fluid tracking as a constrained flow problem all while working in an efficient, modular framework. In addition it lets us tightly couple fluid simulation into flow tracking, providing a global prior that significantly increases tracking accuracy and temporal coherence as compared to previous techniques. We demonstrate how we can use these improved results for a variety of applications, such as re-simulation, detail enhancement, and domain modification. We furthermore give an outlook of the applications beyond fluid tracking that our proximal operator framework could enable by exploring the connection of deblurring and fluid guiding.

An adaptive mesh refinement-multiphase lattice Boltzmann flux solver for simulation of complex binary fluid flows

Science.gov (United States)

Yuan, H. Z.; Wang, Y.; Shu, C.

2017-12-01

This paper presents an adaptive mesh refinement-multiphase lattice Boltzmann flux solver (AMR-MLBFS) for effective simulation of complex binary fluid flows at large density ratios. In this method, an AMR algorithm is proposed by introducing a simple indicator on the root block for grid refinement and two possible statuses for each block. Unlike available block-structured AMR methods, which refine their mesh by spawning or removing four child blocks simultaneously, the present method is able to refine its mesh locally by spawning or removing one to four child blocks independently when the refinement indicator is triggered. As a result, the AMR mesh used in this work can be more focused on the flow region near the phase interface and its size is further reduced. In each block of mesh, the recently proposed MLBFS is applied for the solution of the flow field and the level-set method is used for capturing the fluid interface. As compared with existing AMR-lattice Boltzmann models, the present method avoids both spatial and temporal interpolations of density distribution functions so that converged solutions on different AMR meshes and uniform grids can be obtained. The proposed method has been successfully validated by simulating a static bubble immersed in another fluid, a falling droplet, instabilities of two-layered fluids, a bubble rising in a box, and a droplet splashing on a thin film with large density ratios and high Reynolds numbers. Good agreement with the theoretical solution, the uniform-grid result, and/or the published data has been achieved. Numerical results also show its effectiveness in saving computational time and virtual memory as compared with computations on uniform meshes.

Lattice Boltzmann model for three-phase viscoelastic fluid flow

Science.gov (United States)

Xie, Chiyu; Lei, Wenhai; Wang, Moran

2018-02-01

A lattice Boltzmann (LB) framework is developed for simulation of three-phase viscoelastic fluid flows in complex geometries. This model is based on a Rothman-Keller type model for immiscible multiphase flows which ensures mass conservation of each component in porous media even for a high density ratio. To account for the viscoelastic effects, the Maxwell constitutive relation is correctly introduced into the momentum equation, which leads to a modified lattice Boltzmann evolution equation for Maxwell fluids by removing the normal but excess viscous term. Our simulation tests indicate that this excess viscous term may induce significant errors. After three benchmark cases, the displacement processes of oil by dispersed polymer are studied as a typical example of three-phase viscoelastic fluid flow. The results show that increasing either the polymer intrinsic viscosity or the elastic modulus will enhance the oil recovery.

Contribution to numerical methods for all Mach flow regimes and to fluid-porous coupling for the simulation of homogeneous two-phase flows in nuclear reactors

International Nuclear Information System (INIS)

Zaza, Chady

2015-01-01

The numerical simulation of steam generators of pressurized water reactors is a complex problem, involving different flow regimes and a wide range of length and time scales. An accidental scenario may be associated with very fast variations of the flow with an important Mach number. In contrast in the nominal regime the flow may be stationary, at low Mach number. Moreover whatever the regime under consideration, the array of U-tubes is modelled by a porous medium in order to avoid taking into account the complex geometry of the steam generator, which entails the issue of the coupling conditions at the interface with the free-fluid. We propose a new pressure-correction scheme for cell-centered finite volumes for solving the compressible Navier-Stokes and Euler equations at all Mach number. The existence of a discrete solution, the consistency of the scheme in the Lax sense and the positivity of the internal energy were proved. Then the scheme was extended to the homogeneous two-phase flow models of the GENEPI code developed at CEA. Lastly a multigrid-AMR algorithm was adapted for using our pressure-correction scheme on adaptive grids. Regarding the second issue addressed in this work, the numerical simulation of a fluid flow over a porous bed involves very different length scales. Macroscopic interface models - such as Ochoa-Tapia-Whitaker or Beavers-Joseph law for a viscous flow - represent the transition region between the free-fluid and the porous region by an interface of discontinuity associated with specific transmission conditions. An extension to the Beavers-Joseph law was proposed for the convective regime. By introducing a jump in the kinetic energy at the interface, we recover an interface condition close to the Beavers-Joseph law but with a non-linear slip coefficient, which depends on the free-fluid velocity at the interface and on the Darcy velocity. The validity of this new transmission condition was assessed with direct numerical simulations at

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.

Parallel Simulation of HGMS of Weakly Magnetic Nanoparticles in Irrotational Flow of Inviscid Fluid

Directory of Open Access Journals (Sweden)

Kanok Hournkumnuard

2014-01-01

Full Text Available The process of high gradient magnetic separation (HGMS using a microferromagnetic wire for capturing weakly magnetic nanoparticles in the irrotational flow of inviscid fluid is simulated by using parallel algorithm developed based on openMP. The two-dimensional problem of particle transport under the influences of magnetic force and fluid flow is considered in an annular domain surrounding the wire with inner radius equal to that of the wire and outer radius equal to various multiples of wire radius. The differential equations governing particle transport are solved numerically as an initial and boundary values problem by using the finite-difference method. Concentration distribution of the particles around the wire is investigated and compared with some previously reported results and shows the good agreement between them. The results show the feasibility of accumulating weakly magnetic nanoparticles in specific regions on the wire surface which is useful for applications in biomedical and environmental works. The speedup of parallel simulation ranges from 1.8 to 21 depending on the number of threads and the domain problem size as well as the number of iterations. With the nature of computing in the application and current multicore technology, it is observed that 4–8 threads are sufficient to obtain the optimized speedup.

Predicting Flow Reversals in a Computational Fluid Dynamics Simulated Thermosyphon Using Data Assimilation.

Science.gov (United States)

Reagan, Andrew J; Dubief, Yves; Dodds, Peter Sheridan; Danforth, Christopher M

2016-01-01

A thermal convection loop is a annular chamber filled with water, heated on the bottom half and cooled on the top half. With sufficiently large forcing of heat, the direction of fluid flow in the loop oscillates chaotically, dynamics analogous to the Earth's weather. As is the case for state-of-the-art weather models, we only observe the statistics over a small region of state space, making prediction difficult. To overcome this challenge, data assimilation (DA) methods, and specifically ensemble methods, use the computational model itself to estimate the uncertainty of the model to optimally combine these observations into an initial condition for predicting the future state. Here, we build and verify four distinct DA methods, and then, we perform a twin model experiment with the computational fluid dynamics simulation of the loop using the Ensemble Transform Kalman Filter (ETKF) to assimilate observations and predict flow reversals. We show that using adaptively shaped localized covariance outperforms static localized covariance with the ETKF, and allows for the use of less observations in predicting flow reversals. We also show that a Dynamic Mode Decomposition (DMD) of the temperature and velocity fields recovers the low dimensional system underlying reversals, finding specific modes which together are predictive of reversal direction.

Predicting Flow Reversals in a Computational Fluid Dynamics Simulated Thermosyphon Using Data Assimilation.

Directory of Open Access Journals (Sweden)

Andrew J Reagan

Full Text Available A thermal convection loop is a annular chamber filled with water, heated on the bottom half and cooled on the top half. With sufficiently large forcing of heat, the direction of fluid flow in the loop oscillates chaotically, dynamics analogous to the Earth's weather. As is the case for state-of-the-art weather models, we only observe the statistics over a small region of state space, making prediction difficult. To overcome this challenge, data assimilation (DA methods, and specifically ensemble methods, use the computational model itself to estimate the uncertainty of the model to optimally combine these observations into an initial condition for predicting the future state. Here, we build and verify four distinct DA methods, and then, we perform a twin model experiment with the computational fluid dynamics simulation of the loop using the Ensemble Transform Kalman Filter (ETKF to assimilate observations and predict flow reversals. We show that using adaptively shaped localized covariance outperforms static localized covariance with the ETKF, and allows for the use of less observations in predicting flow reversals. We also show that a Dynamic Mode Decomposition (DMD of the temperature and velocity fields recovers the low dimensional system underlying reversals, finding specific modes which together are predictive of reversal direction.

Simulation of fluid flow in fractured rock: a probabilistic approach

International Nuclear Information System (INIS)

Samaniego, J.A.; Priest, S.D.

1985-02-01

This report describes the results of a research project designed to investigate the influence of discontinuities on fluid flow through fractured rock masses. The aim has been to provide a rational basis for the assessment of prospective intermediate level radioactive waste repository sites. The results of this work are presented in the form of two groups of FORTRAN computer programs. The first of these is designed to process data obtained from exposed rock faces and thereby provide an unbiased estimate of discontinuity characteristics. The resulting data are input to the second group of programs which generate a two-dimensional random realisation of discontinuity geometry. When appropriate boundary conditions have been specified, the resulting network of channels is solved numerically to determine nodal potentials, flow quantities and equivalent permeabilities. A number of validation runs are presented, together with some parametric studies, to investigate the influence of excavation size and discontinuity geometry on fluid flow. A practical application is given in the form of a case study involving the prediction of fluid flow into a 2.8 m diameter tunnel in water bearing, fractured rock. Finally, the applications and limitations of the programs in site assessment for radioactive waste repositories are discussed. (author)

An Efficient and Robust Method for Lagrangian Magnetic Particle Tracking in Fluid Flow Simulations on Unstructured Grids

NARCIS (Netherlands)

Cohen Stuart, D.C.; Kleijn, C.R.; Kenjeres, S.

2010-01-01

In this paper we report on a newly developed particle tracking scheme for fluid flow simulations on 3D unstructured grids, aiming to provide detailed insights in the particle behaviour in complex geometries. A possible field of applications is the Magnetic Drug Targeting (MDT) technique, on which

TOUGH Simulations of the Updegraff's Set of Fluid and Heat Flow Problems

Energy Technology Data Exchange (ETDEWEB)

Moridis, G.J.; Pruess (editor), K.

1992-11-01

The TOUGH code [Pruess, 1987] for two-phase flow of water, air, and heat in penneable media has been exercised on a suite of test problems originally selected and simulated by C. D. Updegraff [1989]. These include five 'verification' problems for which analytical or numerical solutions are available, and three 'validation' problems that model laboratory fluid and heat flow experiments. All problems could be run without any code modifications (*). Good and efficient numerical performance, as well as accurate results were obtained throughout. Additional code verification and validation problems from the literature are briefly summarized, and suggestions are given for proper applications of TOUGH and related codes.

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

MicroHH 1.0: a computational fluid dynamics code for direct numerical simulation and large-eddy simulation of atmospheric boundary layer flows

Science.gov (United States)

van Heerwaarden, Chiel C.; van Stratum, Bart J. H.; Heus, Thijs; Gibbs, Jeremy A.; Fedorovich, Evgeni; Mellado, Juan Pedro

2017-08-01

This paper describes MicroHH 1.0, a new and open-source (www.microhh.org) computational fluid dynamics code for the simulation of turbulent flows in the atmosphere. It is primarily made for direct numerical simulation but also supports large-eddy simulation (LES). The paper covers the description of the governing equations, their numerical implementation, and the parameterizations included in the code. Furthermore, the paper presents the validation of the dynamical core in the form of convergence and conservation tests, and comparison of simulations of channel flows and slope flows against well-established test cases. The full numerical model, including the associated parameterizations for LES, has been tested for a set of cases under stable and unstable conditions, under the Boussinesq and anelastic approximations, and with dry and moist convection under stationary and time-varying boundary conditions. The paper presents performance tests showing good scaling from 256 to 32 768 processes. The graphical processing unit (GPU)-enabled version of the code can reach a speedup of more than an order of magnitude for simulations that fit in the memory of a single GPU.

Numerical simulation of effects of a non-ionized fluid injection and suction on the MHD flow in a circular channel

International Nuclear Information System (INIS)

Rahimi Eosboee, M.; Pourmahmoud, N.; Mirzaie, I.; Mohajeri Khameneh, P.; Majidyfar, S.

2012-01-01

Control of a fluid flow velocity profile by injection and suction of a non-ionized fluid in presence of a uniform steady magnetic field has important technical applications. In this paper, the unsteady incompressible and viscous conducting fluid flow has been investigated in a circular channel. The channel walls are assumed to be non-conducting and porous. They are subjected to a uniform steady magnetic field which is perpendicular to the axis of channel, then and suction and injection are applied at the walls. The well known equations of Magnetohydrodynamics are governed to the motion of an electrically conducting fluid flow that is subjected to magnetic field. The numerical solution is carried out by finite difference approach. The results of present numerical simulation shown that the flow injection and suction through the wall can be controlled effectively, the main flow in channel especially in industrial purposes. The results are obtained for different values of the injected and sucked non-ionized flow rate and the effect of Hartman number on the velocity profile is investigated. Finally, a good agreement is seen between the presented results and the corresponding data of finite element method.

Experimental study of MHD effects on turbulent flow of flibe simulant fluid in a circular pipe

International Nuclear Information System (INIS)

Takeuchi, Junichi; Morley, N.B.; Abdou, M.A.; Satake, Shin-ichi; Yokomine, Takehiko

2007-01-01

Experimental studies of MHD turbulent pipe flow of Flibe simulant fluid have been conducted as a part of US-Japan JUPITER-II collaboration. Flibe is considered as a promising candidate for coolant and tritium breeder in some fusion reactor design concepts because of its low electrical conductivity compared to liquid metals. This reduces the MHD pressure drop to a negligible level; however, turbulence can be significantly suppressed by MHD effects in fusion reactor magnetic field conditions. Heat transfer in the Flibe coolant is characterized by its high Prandtl number. In order to achieve sufficient heat transfer and to prevent localized heat concentration in a high Prandtl number coolant, high turbulence is essential. Even though accurate prediction of the MHD effects on heat transfer for high Prandtl number fluids in the fusion environment is very important, reliable data is not available. In these experiments, an aqueous solution of potassium hydroxide is used as a simulant fluid for Flibe. This paper presents the experimental results obtained by flow field measurement using particle image velocimetry (PIV) technique. The PIV measurements provide 2-dimensional 2-velocity component information on the MHD flow field. The test section is a circular pipe with 89 mm inner diameter and 7.0 m in length, which is 79 times pipe diameter. This relatively large diameter pipe is selected in order to maximize the MHD effects measured by Hartmann number (Ha=BL(sigma/mu)1/2), and to allow better resolution of the flow in the near-wall region. The test section is placed under maximum 2 Tesla magnetic fields for 1.4m of the axial length. The hydrodynamic developing length under the magnetic field is expected to be 1.2 m. In order to apply PIV technique in the magnetic field condition, special optical devices and visualization sections were created. PIV measurements are performed for Re = 11600 with variable Hartmann numbers. The turbulence statistics of the MHD turbulent flow

Formation of the bottom-simulating reflector and its link to vertical fluid flow

Energy Technology Data Exchange (ETDEWEB)

Haacke, R.R.; Hyndman, R.D. [Natural Resources Canada, Sidney, BC (Canada). Geological Survey of Canada, Pacific Geoscience Centre; Westbrook, G.K. [Birmingham Univ., Edgbaston (United Kingdom). Dept. of Geography, Earth and Environmental Sciences

2008-07-01

Natural gas hydrates typically occur with a bottom-simulating reflector (BSR) marking the base of its hydrate stability field. This paper outlined the 2 most important mechanisms that produce free gas beneath the gas hydrate stability zone (GHSZ), consequently producing the BSR. It discussed the importance of hydrate recycling and the solubility-curvature mechanisms in different tectonic environments. It also explained why some areas, such as the Mackenzie Delta in the Canadian Arctic or the northern Gulf of Mexico, have natural gas hydrates without an underlying free-gas zone (FGZ) and associated BSR. The BSR is created primarily by the presence of low-velocity free gas in the pore space under the stability field. This paper focused on the widespread, diffuse distribution of natural gas hydrate in relatively low concentrations that is produced by the vertical migration of gas-rich fluids. The FGZ that occurs under the BSR achieves a steady-state thickness that depends on the diffuse, vertical fluid flux in the system. The opposite is also true, notably if the steady-state thickness of the FGZ can be measured, then the diffuse vertical fluid flux can be estimated. The presence of free gas is easier to detect than gas hydrate because of its very low seismic velocity. This enables the measurement of vertical fluid flux using geophysical methods. The regional hydrate concentration can then be predicted. This study showed that if the gas-water solubility decreases downward beneath the GHSZ, low rates of upward fluid flow enable pore water to become saturated in a thick layer beneath the GHSZ. The FGZ that this produces achieves a steady-state thickness that is sensitive to the rate of upward fluid flow. Geophysical observations that constrain the thickness of sub-BSR FGZs can therefore be used to estimate the regional, diffuse, upward fluid flux through natural gas-hydrate systems. 23 refs., 6 figs.

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.

Numerical methodologies for investigation of moderate-velocity flow using a hybrid computational fluid dynamics - molecular dynamics simulation approach

International Nuclear Information System (INIS)

Ko, Soon Heum; Kim, Na Yong; Nikitopoulos, Dimitris E.; Moldovan, Dorel; Jha, Shantenu

2014-01-01

Numerical approaches are presented to minimize the statistical errors inherently present due to finite sampling and the presence of thermal fluctuations in the molecular region of a hybrid computational fluid dynamics (CFD) - molecular dynamics (MD) flow solution. Near the fluid-solid interface the hybrid CFD-MD simulation approach provides a more accurate solution, especially in the presence of significant molecular-level phenomena, than the traditional continuum-based simulation techniques. It also involves less computational cost than the pure particle-based MD. Despite these advantages the hybrid CFD-MD methodology has been applied mostly in flow studies at high velocities, mainly because of the higher statistical errors associated with low velocities. As an alternative to the costly increase of the size of the MD region to decrease statistical errors, we investigate a few numerical approaches that reduce sampling noise of the solution at moderate-velocities. These methods are based on sampling of multiple simulation replicas and linear regression of multiple spatial/temporal samples. We discuss the advantages and disadvantages of each technique in the perspective of solution accuracy and computational cost.

Sensitivity analysis of Immersed Boundary Method simulations of fluid flow in dense polydisperse random grain packings

Directory of Open Access Journals (Sweden)

Knight Chris

2017-01-01

Full Text Available Polydisperse granular materials are ubiquitous in nature and industry. Despite this, knowledge of the momentum coupling between the fluid and solid phases in dense saturated grain packings comes almost exclusively from empirical correlations [2–4, 8] with monosized media. The Immersed Boundary Method (IBM is a Computational Fluid Dynamics (CFD modelling technique capable of resolving pore scale fluid flow and fluid-particle interaction forces in polydisperse media at the grain scale. Validation of the IBM in the low Reynolds number, high concentration limit was performed by comparing simulations of flow through ordered arrays of spheres with the boundary integral results of Zick and Homsy [10]. Random grain packings were studied with linearly graded particle size distributions with a range of coefficient of uniformity values (Cu = 1.01, 1.50, and 2.00 at a range of concentrations (ϕ ∈ [0.396; 0.681] in order to investigate the influence of polydispersity on drag and permeability. The sensitivity of the IBM results to the choice of radius retraction parameter [1] was investigated and a comparison was made between the predicted forces and the widely used Ergun correlation [3].

Flow Simulation and Performance Prediction of Centrifugal Pumps ...

African Journals Online (AJOL)

With the aid of computational fluid dynamics, the complex internal flows in water pump impellers can be well predicted, thus facilitating the product development process of pumps. In this paper a commercial CFD code was used to solve the governing equations of the flow field. A 2-D simulation of turbulent fluid flow is ...

A computational fluid dynamics and effectiveness-NTU based co-simulation approach for flow mal-distribution analysis in microchannel heat exchanger headers

International Nuclear Information System (INIS)

Huang, Long; Lee, Moon Soo; Saleh, Khaled; Aute, Vikrant; Radermacher, Reinhard

2014-01-01

Refrigerant flow mal-distribution is a practical challenge in most microchannel heat exchangers (MCHXs) applications. Geometry design, uneven heat transfer and pressure drop in the different microchannel tubes are three main reasons leading to the flow mal-distribution. To efficiently and accurately account for these three effects, a new MCHX co-simulation approach is proposed in this paper. The proposed approach combines a detailed header simulation based on computational fluid dynamics (CFD) and a robust effectiveness-based finite volume tube-side heat transfer and refrigerant flow modeling tool. The co-simulation concept is demonstrated on a ten-tube MCHX case study. Gravity effect and uneven airflow effect were numerically analyzed using both water and condensing R134a as the working fluids. The approach was validated against experimental data for an automotive R134a condenser. The inlet header was cut open after the experimental data had been collected. The detailed header geometry was reproduced using the proposed CFD header model. Good prediction accuracy was achieved compared to the experimental data. The presented co-simulation approach is capable of predicting detailed refrigerant flow behavior while accurately predicts the overall heat exchanger performance. - Highlights: •MCHX header flow distribution is analyzed by a co-simulation approach. •The proposed method is capable of simulating both single-phase and two-phase flow. •An actual header geometry is reproduced in the CFD header model. •The modeling work is experimentally validated with good accuracy. •Gravity effect and air side mal-distribution are accounted for

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

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

Generalized Fluid System Simulation Program (GFSSP) - Version 6

Science.gov (United States)

Majumdar, Alok; LeClair, Andre; Moore, Ric; Schallhorn, Paul

2015-01-01

The Generalized Fluid System Simulation Program (GFSSP) is a finite-volume based general-purpose computer program for analyzing steady state and time-dependent flow rates, pressures, temperatures, and concentrations in a complex flow network. The program is capable of modeling real fluids with phase changes, compressibility, mixture thermodynamics, conjugate heat transfer between solid and fluid, fluid transients, pumps, compressors, flow control valves and external body forces such as gravity and centrifugal. The thermo-fluid system to be analyzed is discretized into nodes, branches, and conductors. The scalar properties such as pressure, temperature, and concentrations are calculated at nodes. Mass flow rates and heat transfer rates are computed in branches and conductors. The graphical user interface allows users to build their models using the 'point, drag, and click' method; the users can also run their models and post-process the results in the same environment. The integrated fluid library supplies thermodynamic and thermo-physical properties of 36 fluids, and 24 different resistance/source options are provided for modeling momentum sources or sinks in the branches. Users can introduce new physics, non-linear and time-dependent boundary conditions through user-subroutine.

A direct force model for Galilean invariant lattice Boltzmann simulation of fluid-particle flows

Science.gov (United States)

Tao, Shi; He, Qing; Chen, Baiman; Yang, Xiaoping; Huang, Simin

The lattice Boltzmann method (LBM) has been widely used in the simulation of particulate flows involving complex moving boundaries. Due to the kinetic background of LBM, the bounce-back (BB) rule and the momentum exchange (ME) method can be easily applied to the solid boundary treatment and the evaluation of fluid-solid interaction force, respectively. However, recently it has been found that both the BB and ME schemes may violate the principle of Galilean invariance (GI). Some modified BB and ME methods have been proposed to reduce the GI error. But these remedies have been recognized subsequently to be inconsistent with Newton’s Third Law. Therefore, contrary to those corrections based on the BB and ME methods, a unified iterative approach is adopted to handle the solid boundary in the present study. Furthermore, a direct force (DF) scheme is proposed to evaluate the fluid-particle interaction force. The methods preserve the efficiency of the BB and ME schemes, and the performance on the accuracy and GI is verified and validated in the test cases of particulate flows with freely moving particles.

Dissipative particle dynamics simulation of fluid motion through an unsaturated fracture and fracture junction

International Nuclear Information System (INIS)

Liu Moubin; Meakin, Paul; Huang Hai

2007-01-01

Multiphase fluid motion in unsaturated fractures and fracture networks involves complicated fluid dynamics, which is difficult to model using grid-based continuum methods. In this paper, the application of dissipative particle dynamics (DPD), a relatively new mesoscale method to simulate fluid motion in unsaturated fractures is described. Unlike the conventional DPD method that employs a purely repulsive conservative (non-dissipative) particle-particle interaction to simulate the behavior of gases, we used conservative particle-particle interactions that combine short-range repulsive and long-range attractive interactions. This new conservative particle-particle interaction allows the behavior of multiphase systems consisting of gases, liquids and solids to be simulated. Our simulation results demonstrate that, for a fracture with flat parallel walls, the DPD method with the new interaction potential function is able to reproduce the hydrodynamic behavior of fully saturated flow, and various unsaturated flow modes including thin film flow, wetting and non-wetting flow. During simulations of flow through a fracture junction, the fracture junction can be fully or partially saturated depending on the wetting property of the fluid, the injection rate and the geometry of the fracture junction. Flow mode switching from a fully saturated flow to a thin film flow can also be observed in the fracture junction

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.

User's Guide of TOUGH2-EGS. A Coupled Geomechanical and Reactive Geochemical Simulator for Fluid and Heat Flow in Enhanced Geothermal Systems Version 1.0

Energy Technology Data Exchange (ETDEWEB)

Fakcharoenphol, Perapon [Colorado School of Mines, Golden, CO (United States); Xiong, Yi [Colorado School of Mines, Golden, CO (United States); Hu, Litang [Colorado School of Mines, Golden, CO (United States); Winterfeld, Philip H. [Colorado School of Mines, Golden, CO (United States); Xu, Tianfu [Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States); Wu, Yu-Shu [Colorado School of Mines, Golden, CO (United States)

2013-05-01

TOUGH2-EGS is a numerical simulation program coupling geomechanics and chemical reactions for fluid and heat flows in porous media and fractured reservoirs of enhanced geothermal systems. The simulator includes the fully-coupled geomechanical (THM) module, the fully-coupled geochemical (THC) module, and the sequentially coupled reactive geochemistry (THMC) module. The fully-coupled flow-geomechanics model is developed from the linear elastic theory for the thermo-poro-elastic system and is formulated with the mean normal stress as well as pore pressure and temperature. The chemical reaction is sequentially coupled after solution of flow equations, which provides the flow velocity and phase saturation for the solute transport calculation at each time step. In addition, reservoir rock properties, such as porosity and permeability, are subjected to change due to rock deformation and chemical reactions. The relationships between rock properties and geomechanical and chemical effects from poro-elasticity theories and empirical correlations are incorporated into the simulator. This report provides the user with detailed information on both mathematical models and instructions for using TOUGH2-EGS for THM, THC or THMC simulations. The mathematical models include the fluid and heat flow equations, geomechanical equation, reactive geochemistry equations, and discretization methods. Although TOUGH2-EGS has the capability for simulating fluid and heat flows coupled with both geomechanical and chemical effects, it is up to the users to select the specific coupling process, such as THM, THC, or THMC in a simulation. There are several example problems illustrating the applications of this program. These example problems are described in details and their input data are presented. The results demonstrate that this program can be used for field-scale geothermal reservoir simulation with fluid and heat flow, geomechanical effect, and chemical reaction in porous and fractured media.

Numerical Simulation of Flow Features and Energy Exchange Physics in Near-Wall Region with Fluid-Structure Interaction

Science.gov (United States)

Zhang, Lixiang; Wang, Wenquan; Guo, Yakun

Large eddy simulation is used to explore flow features and energy exchange physics between turbulent flow and structure vibration in the near-wall region with fluid-structure interaction (FSI). The statistical turbulence characteristics in the near-wall region of a vibrating wall, such as the skin frictional coefficient, velocity, pressure, vortices, and the coherent structures have been studied for an aerofoil blade passage of a true three-dimensional hydroturbine. The results show that (i) FSI greatly strengthens the turbulence in the inner region of y+ < 25; and (ii) the energy exchange mechanism between the flow and the vibration depends strongly on the vibration-induced vorticity in the inner region. The structural vibration provokes a frequent action between the low- and high-speed streaks to balance the energy deficit caused by the vibration. The velocity profile in the inner layer near the vibrating wall has a significant distinctness, and the viscosity effect of the fluid in the inner region decreases due to the vibration. The flow features in the inner layer are altered by a suitable wall vibration.

CFD Simulation of Heat and Fluid Flow for Spent Fuel in a Dry Storage

International Nuclear Information System (INIS)

In, Wangkee; Kwack, Youngkyun; Kook, Donghak; Koo, Yanghyun

2014-01-01

A dry storage system is used for the interim storage of spent fuel prior to permanent depository and/or recycling. The spent fuel is initially stored in a water pool for more than 5 years at least after dispatch from the reactor core and is transported to dry storage. The dry cask contains a multiple number of spent fuel assemblies, which are cooled down in the spent fuel pool. The dry cask is usually filled up with helium gas for increasing the heat transfer to the environment outside the cask. The dry storage system has been used for more than a decade in United States of America (USA) and the European Union (EU). Korea is also developing a dry storage system since its spent fuel pool is anticipated to be full within 10 years. The spent fuel will be stored in a dry cask for more than 40 years. The integrity and safety of spent fuel are important for long-term dry storage. The long-term storage will experience the degradation of spent fuel such as the embrittlement of fuel cladding, thermal creep and hydride reorientation. High burn-up fuel may expedite the material degradation. It is known that the cladding temperature has a strong influence on the material degradation. Hence, it is necessary to accurately predict the local distribution of the cladding temperature using the Computational Fluid Dynamics (CFD) approach. The objective of this study is to apply the CFD method for predicting the three-dimensional distribution of fuel temperature in a dry cask. This CFD study simulated the dry cask for containing the 21 fuel assemblies under development in Korea. This paper presents the fluid velocity and temperature distribution as well as the fuel temperature. A two-step CFD approach was applied to simulate the heat and fluid flow in a dry storage of 21 spent fuel assemblies. The first CFD analysis predicted the helium flow and temperature in a dry cask by a assuming porous body of the spent fuel. The second CFD analysis was to simulate a spent fuel assembly in the

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.

Linking fault permeability, fluid flow, and earthquake triggering in a hydrothermally active tectonic setting: Numerical Simulations of the hydrodynamics in the Tjörnes Fracture Zone, Iceland.

Science.gov (United States)

Lupi, M.; Geiger, S.; Graham, C.; Claesson, L.; Richter, B.

2007-12-01

A good insight into the transient fluid flow evolution within a hydrothermal system is of primary importance for the understanding of several geologic processes, for example the hydrodynamic triggering of earthquakes or the formation of mineral deposits. The strong permeability contrast between different crustal layers as well as the high geothermal gradient of these areas are elements that strongly affect the flow behaviour. In addition, the sudden and transient occurrence of joints, faults and magmatic intrusions are likely to change the hydrothermal flow paths in very short time. The Tjörnes Fracture Zone (TFZ) north of Iceland, is such a hydrothermal area where a high geothermal gradient, magmatic bodies, faults, and the strong contrast between sediments and fractured lava layers govern the large-scale fluid flow. The TFZ offsets the Kolbeinsey Ridge and the Northern Rift Zone. It is characterized by km-scale faults that link sub-seafloor sediments and lava layers with deeper crystalline rocks. These structures focus fluid flow and allow for the mixing between cold seawater and deep hydrothermal fluids. A strong seismic activity is present in the TFZ: earthquakes up to magnitude 7 have been recorded over the past years. Hydrogeochemical changes before, during and after a magnitude 5.8 earthquake suggest that the evolving stress state before the earthquake leads to (remote) permeability variations, which alter the fluid flow paths. This is in agreement with recent numerical fluid flow simulations which demonstrate that fluid flow in magmatic- hydrothermal systems is often convective and very sensitive to small variations in permeability. In order to understand the transient fluid flow behaviour in this complex geological environment, we have conducted numerical simulations of heat and mass transport in two geologically realistic cross-sectional models of the TFZ. The geologic models are discretised using finite element and finite volume methods. They hence have

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.

Numerical simulation of complex multi-phase fluid of casting process and its applications

Directory of Open Access Journals (Sweden)

CHEN Li-liang

2006-05-01

Full Text Available The fluid of casting process is a typical kind of multi-phase flow. Actually, many casting phenomena have close relationship with the multi-phase flow, such as molten metal filling process, air entrapment, slag movement, venting process of die casting, gas escaping of lost foam casting and so on. Obviously, in order to analyze these phenomena accurately, numerical simulation of the multi-phase fluid is necessary. Unfortunately, so far, most of the commercial casting simulation systems do not have the ability of multi-phase flow modeling due to the difficulty in the multi-phase flow calculation. In the paper, Finite Different Method (FDM technique was adopt to solve the multi-phase fluid model. And a simple object of the muiti-phase fluid was analyzed to obtain the fluid rates of the liquid phase and the entrapped air phase.

High viscosity fluid simulation using particle-based method

KAUST Repository

Chang, Yuanzhang

2011-03-01

We present a new particle-based method for high viscosity fluid simulation. In the method, a new elastic stress term, which is derived from a modified form of the Hooke\\'s law, is included in the traditional Navier-Stokes equation to simulate the movements of the high viscosity fluids. Benefiting from the Lagrangian nature of Smoothed Particle Hydrodynamics method, large flow deformation can be well handled easily and naturally. In addition, in order to eliminate the particle deficiency problem near the boundary, ghost particles are employed to enforce the solid boundary condition. Compared with Finite Element Methods with complicated and time-consuming remeshing operations, our method is much more straightforward to implement. Moreover, our method doesn\\'t need to store and compare to an initial rest state. The experimental results show that the proposed method is effective and efficient to handle the movements of highly viscous flows, and a large variety of different kinds of fluid behaviors can be well simulated by adjusting just one parameter. © 2011 IEEE.

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.

Status of the TOUGH-FLAC simulator and recent applications related to coupled fluid flow and crustal deformations

Energy Technology Data Exchange (ETDEWEB)

Rutqvist, J.

2010-06-01

This paper presents recent advancement in and applications of TOUGH-FLAC, a simulator for multiphase fluid flow and geomechanics. The TOUGH-FLAC simulator links the TOUGH family multiphase fluid and heat transport codes with the commercial FLAC{sup 3D} geomechanical simulator. The most significant new TOUGH-FLAC development in the past few years is a revised architecture, enabling a more rigorous and tight coupling procedure with improved computational efficiency. The applications presented in this paper are related to modeling of crustal deformations caused by deep underground fluid movements and pressure changes as a result of both industrial activities (the In Salah CO{sub 2} Storage Project and the Geysers Geothermal Field) and natural events (the 1960s Matsushiro Earthquake Swarm). Finally, the paper provides some perspectives on the future of TOUGH-FLAC in light of its applicability to practical problems and the need for high-performance computing capabilities for field-scale problems, such as industrial-scale CO{sub 2} storage and enhanced geothermal systems. It is concluded that despite some limitations to fully adapting a commercial code such as FLAC{sup 3D} for some specialized research and computational needs, TOUGH-FLAC is likely to remain a pragmatic simulation approach, with an increasing number of users in both academia and industry.

Numerical Modelling of Three-Fluid Flow Using The Level-set Method

Science.gov (United States)

Li, Hongying; Lou, Jing; Shang, Zhi

2014-11-01

This work presents a numerical model for simulation of three-fluid flow involving two different moving interfaces. These interfaces are captured using the level-set method via two different level-set functions. A combined formulation with only one set of conservation equations for the whole physical domain, consisting of the three different immiscible fluids, is employed. Numerical solution is performed on a fixed mesh using the finite volume method. Surface tension effect is incorporated using the Continuum Surface Force model. Validation of the present model is made against available results for stratified flow and rising bubble in a container with a free surface. Applications of the present model are demonstrated by a variety of three-fluid flow systems including (1) three-fluid stratified flow, (2) two-fluid stratified flow carrying the third fluid in the form of drops and (3) simultaneous rising and settling of two drops in a stationary third fluid. The work is supported by a Thematic and Strategic Research from A*STAR, Singapore (Ref. #: 1021640075).

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.

Simulations of magnetorheological suspensions in Poiseuille flow

Energy Technology Data Exchange (ETDEWEB)

Pappas, Yannis; Klingenberg, Daniel J. [University of Wisconsin-Madison, Department of Chemical and Biological Engineering, Madison, WI (United States)

2006-06-15

Particle-level simulations are conducted to study magnetorheological fluids in plane Poiseuille flow. The importance of the boundary conditions for the particles at the channel walls is examined by considering two extreme cases: no friction and infinite coefficient of friction. The inclusion of friction produces Bingham fluid behavior, as commonly observed experimentally for MR suspensions. Lamellar structures, similar to those reported for electrorheological fluids in shear flow, are observed in the post-yield region for both particle boundary conditions. The formation of these lamellae is accompanied by an increase in the bulk fluid velocity. The slip boundary condition produces higher fluid velocities and thicker lamellar structures. (orig.)

Parallel computing simulation of fluid flow in the unsaturated zone of Yucca Mountain, Nevada

International Nuclear Information System (INIS)

Zhang, Keni; Wu, Yu-Shu; Bodvarsson, G.S.

2001-01-01

This paper presents the application of parallel computing techniques to large-scale modeling of fluid flow in the unsaturated zone (UZ) at Yucca Mountain, Nevada. In this study, parallel computing techniques, as implemented into the TOUGH2 code, are applied in large-scale numerical simulations on a distributed-memory parallel computer. The modeling study has been conducted using an over-one-million-cell three-dimensional numerical model, which incorporates a wide variety of field data for the highly heterogeneous fractured formation at Yucca Mountain. The objective of this study is to analyze the impact of various surface infiltration scenarios (under current and possible future climates) on flow through the UZ system, using various hydrogeological conceptual models with refined grids. The results indicate that the one-million-cell models produce better resolution results and reveal some flow patterns that cannot be obtained using coarse-grid modeling models

Numerical simulation of fluid flow in a rotational bioreactor

Science.gov (United States)

Ganimedov, V. L.; Papaeva, E. O.; Maslov, N. A.; Larionov, P. M.

2017-10-01

Application of scaffold technology for the problem of bone tissue regeneration has great prospects in modern medicine. The influence of fluid shear stress on stem cells cultivation and its differentiation into osteoblasts is the subject of intensive research. Mathematical modeling of fluid flow in bioreactor allowed us to determine the structure of flow and estimate the level of mechanical stress on cells. The series of computations for different rotation frequencies (0.083, 0.124, 0.167, 0.2 and 0.233 Hz) was performed for the laminar flow regime approximation. It was shown that the Taylor vortices in the gap between the cylinders qualitatively change the distribution of static pressure and shear stress in the region of vortices connection. It was shown that an increase in the rotation frequency leads to an increase of the unevenness in distribution of the above mentioned functions. The obtained shear stress and static pressure dependence on the rotational frequency make it possible to choose the operating mode of the reactor depending on the provided requirements. It was shown that in the range of rotation frequencies chosen in this work (0.083 < f < 0.233 Hz), the shear stress does not exceed the known literature data (0.002 - 0.1 Pa).

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.

Hydrodynamic cavitation in Stokes flow of anisotropic fluids

OpenAIRE

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

2017-01-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 domai...

Numerical solver for compressible two-fluid flow

NARCIS (Netherlands)

J. Naber (Jorick)

2005-01-01

textabstractThis report treats the development of a numerical solver for the simulation of flows of two non-mixing fluids described by the two-dimensional Euler equations. A level-set equation in conservative form describes the interface. After each time step the deformed level-set function is

Cellular-automation fluids: A model for flow in porous media

International Nuclear Information System (INIS)

Rothman, D.H.

1987-01-01

Because the intrinsic inhomogeneity of porous media makes the application of proper boundary conditions difficult, fluid flow through microgeometric models has typically been achieved with idealized arrays of geometrically simple pores, throats, and cracks. The author proposes here an attractive alternative, capable of freely and accurately modeling fluid flow in grossly irregular geometries. This new method numerically solves the Navier-Stokes equations using the cellular-automation fluid model introduced by Frisch, Hasslacher, and Pomeau. The cellular-automation fluid is extraordinarily simple - particles of unit mass traveling with unit velocity reside on a triangular lattice and obey elementary collisions rules - but capable of modeling much of the rich complexity of real fluid flow. The author shows how cellular-automation fluids are applied to the study of porous media. In particular, he discusses issues of scale on the cellular-automation lattice and present the results of 2-D simulations, including numerical estimation of permeability and verification of Darcy's law

CFD simulation of subcooled flow boiling at low pressure

International Nuclear Information System (INIS)

Koncar, B.; Mavko, B.

2001-01-01

An increased interest to numerically simulate the subcooled flow boiling at low pressures (1 to 10 bar) has been aroused in recent years, pursued by the need to perform safety analyses of research nuclear reactors and to investigate the sump cooling concept for future light water reactors. In this paper the subcooled flow boiling has been simulated with a multidimensional two-fluid model used in a CFX-4.3 computational fluid dynamics (CFD) code. The existing model was adequately modified for low pressure conditions. It was shown that interfacial forces, which are usually used for adiabatic flows, need to be modeled to simulate subcooled boiling at low pressure conditions. Simulation results are compared against published experimental data [1] and agree well with experiments.(author)

Lattice Boltzmann model for simulating immiscible two-phase flows

International Nuclear Information System (INIS)

Reis, T; Phillips, T N

2007-01-01

The lattice Boltzmann equation is often promoted as a numerical simulation tool that is particularly suitable for predicting the flow of complex fluids. This paper develops a two-dimensional 9-velocity (D2Q9) lattice Boltzmann model for immiscible binary fluids with variable viscosities and density ratio using a single relaxation time for each fluid. In the macroscopic limit, this model is shown to recover the Navier-Stokes equations for two-phase flows. This is achieved by constructing a two-phase component of the collision operator that induces the appropriate surface tension term in the macroscopic equations. A theoretical expression for surface tension is determined. The validity of this analysis is confirmed by comparing numerical and theoretical predictions of surface tension as a function of density. The model is also shown to predict Laplace's law for surface tension and Poiseuille flow of layered immiscible binary fluids. The spinodal decomposition of two fluids of equal density but different viscosity is then studied. At equilibrium, the system comprises one large low viscosity bubble enclosed by the more viscous fluid in agreement with theoretical arguments of Renardy and Joseph (1993 Fundamentals of Two-Fluid Dynamics (New York: Springer)). Two other simulations, namely the non-equilibrium rod rest and the coalescence of two bubbles, are performed to show that this model can be used to simulate two fluids with a large density ratio

Numerical simulation of heat transfer and fluid flow in a DC plasma-arc device for waste thermal treatment

International Nuclear Information System (INIS)

Deng, Jing; Li, Yaojian; Xu, Yongxiang; Sheng, Hongzhi

2010-01-01

In this work, Magnetic Fluid dynamics (MHD) model is used to stimulate the electromagnetic field, heat transfer and fluid flow in a DC non-transferred arc plasma torch. Through the coupled iterative computation about the electromagnetic equations described by magnetic vector potential format and the modified fluid dynamics equations, the electric potential, temperature and velocity distributions in the torch are obtained. The fluid-solid coupled computation method is applied to treat the electric current and heat transfer at the interface between the electrodes and fluid. The location of arc root attachment at the inside surface of anode and the arc voltage of the torch that we have predicted are very consistent with the corresponding experimental results. The calculated results of the torch are applied to the numerical simulation of the plasma jets under the laminar and turbulent condition. (author)

Modeling of Compressible Flow with Friction and Heat Transfer Using the Generalized Fluid System Simulation Program (GFSSP)

Science.gov (United States)

Bandyopadhyay, Alak; Majumdar, Alok

2007-01-01

The present paper describes the verification and validation of a quasi one-dimensional pressure based finite volume algorithm, implemented in Generalized Fluid System Simulation Program (GFSSP), for predicting compressible flow with friction, heat transfer and area change. The numerical predictions were compared with two classical solutions of compressible flow, i.e. Fanno and Rayleigh flow. Fanno flow provides an analytical solution of compressible flow in a long slender pipe where incoming subsonic flow can be choked due to friction. On the other hand, Raleigh flow provides analytical solution of frictionless compressible flow with heat transfer where incoming subsonic flow can be choked at the outlet boundary with heat addition to the control volume. Nonuniform grid distribution improves the accuracy of numerical prediction. A benchmark numerical solution of compressible flow in a converging-diverging nozzle with friction and heat transfer has been developed to verify GFSSP's numerical predictions. The numerical predictions compare favorably in all cases.

Numerical Modeling of Porous Structure of Biomaterial and Fluid Flowing Through Biomaterial

Institute of Scientific and Technical Information of China (English)

无

2005-01-01

A Cellular Automata model of simulating body fluid flowing into porous bioceramic implants generated with stochastic methods is described, of which main parameters and evolvement rule are determined in terms of flow behavior of body fluid in porous biomaterials. The model is implemented by GUI( Graphical User Interface) program in MATLAB, and the results of numerical modeling show that the body fluid percolation is related to the size of pores and porosity.

SINDA/SINFLO computer routine, volume 1, revision A. [for fluid flow system analysis

Science.gov (United States)

Oren, J. A.; Williams, D. R.

1975-01-01

The SINFLO package was developed to modify the SINDA preprocessor to accept and store the input data for fluid flow systems analysis and adding the FLOSOL user subroutine to perform the flow solution. This reduced and simplified the user input required for analysis of flow problems. A temperature calculation method, the flow-hybrid method which was developed in previous VSD thermal simulator routines, was incorporated for calculating fluid temperatures. The calculation method accuracy was improved by using fluid enthalpy rather than specific heat for the convective term of the fluid temperature equation. Subroutines and data input requirements are described along with user subroutines, flow data storage, and usage of the plot program.

Streamline topology: Patterns in fluid flows and their bifurcations

DEFF Research Database (Denmark)

Brøns, Morten

2007-01-01

Using dynamical systems theory, we consider structures such as vortices and separation in the streamline patterns of fluid flows. Bifurcation of patterns under variation of external parameters is studied using simplifying normal form transformations. Flows away from boundaries, flows close to fix...... walls, and axisymmetric flows are analyzed in detail. We show how to apply the ideas from the theory to analyze numerical simulations of the vortex breakdown in a closed cylindrical container....

User's guide of TOUGH2-EGS-MP: A Massively Parallel Simulator with Coupled Geomechanics for Fluid and Heat Flow in Enhanced Geothermal Systems VERSION 1.0

Energy Technology Data Exchange (ETDEWEB)

Xiong, Yi [Colorado School of Mines, Golden, CO (United States); Fakcharoenphol, Perapon [Colorado School of Mines, Golden, CO (United States); Wang, Shihao [Colorado School of Mines, Golden, CO (United States); Winterfeld, Philip H. [Colorado School of Mines, Golden, CO (United States); Zhang, Keni [Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States); Wu, Yu-Shu [Colorado School of Mines, Golden, CO (United States)

2013-12-01

TOUGH2-EGS-MP is a parallel numerical simulation program coupling geomechanics with fluid and heat flow in fractured and porous media, and is applicable for simulation of enhanced geothermal systems (EGS). TOUGH2-EGS-MP is based on the TOUGH2-MP code, the massively parallel version of TOUGH2. In TOUGH2-EGS-MP, the fully-coupled flow-geomechanics model is developed from linear elastic theory for thermo-poro-elastic systems and is formulated in terms of mean normal stress as well as pore pressure and temperature. Reservoir rock properties such as porosity and permeability depend on rock deformation, and the relationships between these two, obtained from poro-elasticity theories and empirical correlations, are incorporated into the simulation. This report provides the user with detailed information on the TOUGH2-EGS-MP mathematical model and instructions for using it for Thermal-Hydrological-Mechanical (THM) simulations. The mathematical model includes the fluid and heat flow equations, geomechanical equation, and discretization of those equations. In addition, the parallel aspects of the code, such as domain partitioning and communication between processors, are also included. Although TOUGH2-EGS-MP has the capability for simulating fluid and heat flows coupled with geomechanical effects, it is up to the user to select the specific coupling process, such as THM or only TH, in a simulation. There are several example problems illustrating applications of this program. These example problems are described in detail and their input data are presented. Their results demonstrate that this program can be used for field-scale geothermal reservoir simulation in porous and fractured media with fluid and heat flow coupled with geomechanical effects.

Measurement bias of fluid velocity in molecular simulations

International Nuclear Information System (INIS)

Tysanner, Martin W.; Garcia, Alejandro L.

2004-01-01

In molecular simulations of fluid flow, the measurement of mean fluid velocity is considered to be a straightforward computation, yet there is some ambiguity in its definition. We show that in systems far from equilibrium, such as those with large temperature or velocity gradients, two commonly used definitions give slightly different results. Specifically, a bias can arise when computing the mean fluid velocity by measuring the mean particle velocity in a cell and averaging this mean over samples. We show that this bias comes from the correlation of momentum and density fluctuations in non-equilibrium fluids, obtain an analytical expression for predicting it, and discuss what system characteristics (e.g., number of particles per cell, temperature gradients) reduce or magnify the error. The bias has a physical origin so although we demonstrate it by direct simulation Monte Carlo (DSMC) computations, the same effect will be observed with other particle-based simulation methods, such as molecular dynamics and lattice gases

Simulation of Tailrace Hydrodynamics Using Computational Fluid Dynamics Models

Energy Technology Data Exchange (ETDEWEB)

Cook, Christopher B.; Richmond, Marshall C.

2001-05-01

This report investigates the feasibility of using computational fluid dynamics (CFD) tools to investigate hydrodynamic flow fields surrounding the tailrace zone below large hydraulic structures. Previous and ongoing studies using CFD tools to simulate gradually varied flow with multiple constituents and forebay/intake hydrodynamics have shown that CFD tools can provide valuable information for hydraulic and biological evaluation of fish passage near hydraulic structures. These studies however are incapable of simulating the rapidly varying flow fields that involving breakup of the free-surface, such as those through and below high flow outfalls and spillways. Although the use of CFD tools for these types of flow are still an active area of research, initial applications discussed in this report show that these tools are capable of simulating the primary features of these highly transient flow fields.

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

Simulation of solid-liquid flows in a stirred bead mill based on computational fluid dynamics (CFD)

Science.gov (United States)

Winardi, S.; Widiyastuti, W.; Septiani, E. L.; Nurtono, T.

2018-05-01

The selection of simulation model is an important step in computational fluid dynamics (CFD) to obtain an agreement with experimental work. In addition, computational time and processor speed also influence the performance of the simulation results. Here, we report the simulation of solid-liquid flow in a bead mill using Eulerian model. Multiple Reference Frame (MRF) was also used to model the interaction between moving (shaft and disk) and stationary (chamber exclude shaft and disk) zones. Bead mill dimension was based on the experimental work of Yamada and Sakai (2013). The effect of shaft rotation speed of 1200 and 1800 rpm on the particle distribution and the flow field was discussed. For rotation speed of 1200 rpm, the particles spread evenly throughout the bead mill chamber. On the other hand, for the rotation speed of 1800 rpm, the particles tend to be thrown to the near wall region resulting in the dead zone and found no particle in the center region. The selected model agreed well to the experimental data with average discrepancies less than 10%. Furthermore, the simulation was run without excessive computational cost.

On fluid flow driven by topography in a librating body

Science.gov (United States)

Wu, C.; Roberts, P. H.

2009-12-01

Currently considerable effort and resources are being devoted to studies of Mercury, the Moon and Europa. Measuring the libration of these bodies can provide significant knowledge about their internal structures and physical properties; see Williams et al., 2001, Peale et al., 2002, Wu et al., 2007. To interpret such observations, it is important to understand better how libration affects the motion of the fluid in their interiors. To this end, Noir et al. (2009) investigated, via laboratory experiments and numerical simulations, the flow in a fluid filling a rotating spherical cavity driven by an axial oscillation of the container about a diameter. More realistically, the cavity is better represented by a triaxial ellipsoid. We may then distinguish between topographic and axisymmetricli libration. The latter refers to libration about a symmetry axis of the container which is therefore only viscously coupled to the fluid. In topographic libration, pressure forces on the boundary also affect the fluid motions in the cavity. We describe results from preliminary studies of topographic libration obtained through numerical simulation of incompressible fluid motion in an oblate spheroidal cavity with a libration axis perpendicular to the symmetry axis of the container. The computer code is a modification of one recently developed to study precessionally-driven flows in a spheroidal body of fluid (Wu and Roberts, 2009). It advances the flow in time using finite differences on overlapping grids; in this way the numerical difficulty known as the pole problem, is completely avoided.

Visualization of numerically simulated aerodynamic flow fields

International Nuclear Information System (INIS)

Hian, Q.L.; Damodaran, M.

1991-01-01

The focus of this paper is to describe the development and the application of an interactive integrated software to visualize numerically simulated aerodynamic flow fields so as to enable the practitioner of computational fluid dynamics to diagnose the numerical simulation and to elucidate essential flow physics from the simulation. The input to the software is the numerical database crunched by a supercomputer and typically consists of flow variables and computational grid geometry. This flow visualization system (FVS), written in C language is targetted at the Personal IRIS Workstations. In order to demonstrate the various visualization modules, the paper also describes the application of this software to visualize two- and three-dimensional flow fields past aerodynamic configurations which have been numerically simulated on the NEC-SXIA Supercomputer. 6 refs

Numerical simulation for gas-liquid two-phase flow in pipe networks

International Nuclear Information System (INIS)

Li Xiaoyan; Kuang Bo; Zhou Guoliang; Xu Jijun

1998-01-01

The complex pipe network characters can not directly presented in single phase flow, gas-liquid two phase flow pressure drop and void rate change model. Apply fluid network theory and computer numerical simulation technology to phase flow pipe networks carried out simulate and compute. Simulate result shows that flow resistance distribution is non-linear in two phase pipe network

Flow MRI simulation in complex 3D geometries: Application to the cerebral venous network.

Science.gov (United States)

Fortin, Alexandre; Salmon, Stéphanie; Baruthio, Joseph; Delbany, Maya; Durand, Emmanuel

2018-02-05

Develop and evaluate a complete tool to include 3D fluid flows in MRI simulation, leveraging from existing software. Simulation of MR spin flow motion is of high interest in the study of flow artifacts and angiography. However, at present, only a few simulators include this option and most are restricted to static tissue imaging. An extension of JEMRIS, one of the most advanced high performance open-source simulation platforms to date, was developed. The implementation of a Lagrangian description of the flow allows simulating any MR experiment, including both static tissues and complex flow data from computational fluid dynamics. Simulations of simple flow models are compared with real experiments on a physical flow phantom. A realistic simulation of 3D flow MRI on the cerebral venous network is also carried out. Simulations and real experiments are in good agreement. The generality of the framework is illustrated in 2D and 3D with some common flow artifacts (misregistration and inflow enhancement) and with the three main angiographic techniques: phase contrast velocimetry (PC), time-of-flight, and contrast-enhanced imaging MRA. The framework provides a versatile and reusable tool for the simulation of any MRI experiment including physiological fluids and arbitrarily complex flow motion. © 2018 International Society for Magnetic Resonance in Medicine.

Simulations of Bingham plastic flows with the multiple-relaxation-time lattice Boltzmann model

Science.gov (United States)

Chen, SongGui; Sun, QiCheng; Jin, Feng; Liu, JianGuo

2014-03-01

Fresh cement mortar is a type of workable paste, which can be well approximated as a Bingham plastic and whose flow behavior is of major concern in engineering. In this paper, Papanastasiou's model for Bingham fluids is solved by using the multiplerelaxation-time lattice Boltzmann model (MRT-LB). Analysis of the stress growth exponent m in Bingham fluid flow simulations shows that Papanastasiou's model provides a good approximation of realistic Bingham plastics for values of m > 108. For lower values of m, Papanastasiou's model is valid for fluids between Bingham and Newtonian fluids. The MRT-LB model is validated by two benchmark problems: 2D steady Poiseuille flows and lid-driven cavity flows. Comparing the numerical results of the velocity distributions with corresponding analytical solutions shows that the MRT-LB model is appropriate for studying Bingham fluids while also providing better numerical stability. We further apply the MRT-LB model to simulate flow through a sudden expansion channel and the flow surrounding a round particle. Besides the rich flow structures obtained in this work, the dynamics fluid force on the round particle is calculated. Results show that both the Reynolds number Re and the Bingham number Bn affect the drag coefficients C D , and a drag coefficient with Re and Bn being taken into account is proposed. The relationship of Bn and the ratio of unyielded zone thickness to particle diameter is also analyzed. Finally, the Bingham fluid flowing around a set of randomly dispersed particles is simulated to obtain the apparent viscosity and velocity fields. These results help simulation of fresh concrete flowing in porous media.

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

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 gas compressible flow by free surface water flow

International Nuclear Information System (INIS)

Altafini, C.R.; Silva Ferreira, R.T. da

1981-01-01

The analogy between the water flow with a free surface and the compressible fluid flow, commonly called hydraulic analogy, is analyzed and its limitations are identified. The water table is the equipment used for this simulation, which allows the quatitative analysis of subsonic and supersonic flow with a low cost apparatus. The hydraulic analogy is applied to subsonic flow around circular cylinders and supersonic flow around cones. The results are compared with available theoretical and experimental data and a good agreement is achieved. (Author) [pt

Methods and models for accelerating dynamic simulation of fluid power circuits

Energy Technology Data Exchange (ETDEWEB)

Aaman, R.

2011-07-01

are proposed. Generally, in hydraulic power transmission systems the orifice flow is clearly in the turbulent area. The flow becomes laminar as the pressure drop over the orifice approaches zero only in rare situations. These are e.g. when a valve is closed, or an actuator is driven against an end stopper, or external force makes actuator to switch its direction during operation. This means that in terms of accuracy, the description of laminar flow is not necessary. But, unfortunately, when a purely turbulent description of the orifice is used, numerical problems occur when the pressure drop comes close to zero since the first derivative of flow with respect to the pressure drop approaches infinity when the pressure drop approaches zero. The second derivative becomes discontinuous, which causes numerical noise and an infinitely small integration step when a variable step integrator is used. A numerically efficient model for the orifice flow is proposed using a cubic spline function to describe the flow in the laminar and transition areas. Parameters for the cubic spline function are selected such that its first derivative is equal to the first derivative of the pure turbulent orifice flow model in the boundary condition. In the dynamic simulation of fluid power circuits, a tradeoff exists between accuracy and calculation speed. This investigation is made for the two-regime flow orifice model. Especially inside of many types of valves, as well as between them, there exist very small volumes. The integration of pressures in small fluid volumes causes numerical problems in fluid power circuit simulation. Particularly in realtime simulation, these numerical problems are a great weakness. The system stiffness approaches infinity as the fluid volume approaches zero. If fixed step explicit algorithms for solving ordinary differential equations (ODE) are used, the system stability would easily be lost when integrating pressures in small volumes. To solve the problem caused by

Numerical Simulation of Steady Supercavitating Flows

OpenAIRE

Ali Jafarian; Ahmad-Reza Pishevar

2016-01-01

In this research, the Supercavitation phenomenon in compressible liquid flows is simulated. The one-fluid method based on a new exact two-phase Riemann solver is used for modeling. The cavitation is considered as an isothermal process and a consistent equation of state with the physical behavior of the water is used. High speed flow of water over a cylinder and a projectile are simulated and the results are compared with the previous numerical and experimental results. The cavitation bubble p...

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)

Measurement with microscopic MRI and simulation of flow in different aneurysm models

Energy Technology Data Exchange (ETDEWEB)

Edelhoff, Daniel, E-mail: daniel.edelhoff@tu-dortmund.de; Frank, Frauke; Heil, Marvin; Suter, Dieter [Experimental Physics III, TU Dortmund University, Otto-Hahn-Street 4, Dortmund 44227 (Germany); Walczak, Lars; Weichert, Frank [Computer Science VII, TU Dortmund University, Otto-Hahn-Street 16, Dortmund 44227 (Germany); Schmitz, Inge [Institute for Pathology, Ruhr Universität Bochum, Bürkle-de-la-Camp-Platz 1, Bochum 44789 (Germany)

2015-10-15

Purpose: The impact and the development of aneurysms depend to a significant degree on the exchange of liquid between the regular vessel and the pathological extension. A better understanding of this process will lead to improved prediction capabilities. The aim of the current study was to investigate fluid-exchange in aneurysm models of different complexities by combining microscopic magnetic resonance measurements with numerical simulations. In order to evaluate the accuracy and applicability of these methods, the fluid-exchange process between the unaltered vessel lumen and the aneurysm phantoms was analyzed quantitatively using high spatial resolution. Methods: Magnetic resonance flow imaging was used to visualize fluid-exchange in two different models produced with a 3D printer. One model of an aneurysm was based on histological findings. The flow distribution in the different models was measured on a microscopic scale using time of flight magnetic resonance imaging. The whole experiment was simulated using fast graphics processing unit-based numerical simulations. The obtained simulation results were compared qualitatively and quantitatively with the magnetic resonance imaging measurements, taking into account flow and spin–lattice relaxation. Results: The results of both presented methods compared well for the used aneurysm models and the chosen flow distributions. The results from the fluid-exchange analysis showed comparable characteristics concerning measurement and simulation. Similar symmetry behavior was observed. Based on these results, the amount of fluid-exchange was calculated. Depending on the geometry of the models, 7% to 45% of the liquid was exchanged per second. Conclusions: The result of the numerical simulations coincides well with the experimentally determined velocity field. The rate of fluid-exchange between vessel and aneurysm was well-predicted. Hence, the results obtained by simulation could be validated by the experiment. The

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

Particle hopping vs. fluid-dynamical models for traffic flow

Energy Technology Data Exchange (ETDEWEB)

Nagel, K.

1995-12-31

Although particle hopping models have been introduced into traffic science in the 19509, their systematic use has only started recently. Two reasons for this are, that they are advantageous on modem computers, and that recent theoretical developments allow analytical understanding of their properties and therefore more confidence for their use. In principle, particle hopping models fit between microscopic models for driving and fluiddynamical models for traffic flow. In this sense, they also help closing the conceptual gap between these two. This paper shows connections between particle hopping models and traffic flow theory. It shows that the hydrodynamical limits of certain particle hopping models correspond to the Lighthill-Whitham theory for traffic flow, and that only slightly more complex particle hopping models produce already the correct traffic jam dynamics, consistent with recent fluid-dynamical models for traffic flow. By doing so, this paper establishes that, on the macroscopic level, particle hopping models are at least as good as fluid-dynamical models. Yet, particle hopping models have at least two advantages over fluid-dynamical models: they straightforwardly allow microscopic simulations, and they include stochasticity.

Generalized Fluid System Simulation Program, Version 5.0-Educational

Science.gov (United States)

Majumdar, A. K.

2011-01-01

The Generalized Fluid System Simulation Program (GFSSP) is a finite-volume based general-purpose computer program for analyzing steady state and time-dependent flow rates, pressures, temperatures, and concentrations in a complex flow network. The program is capable of modeling real fluids with phase changes, compressibility, mixture thermodynamics, conjugate heat transfer between solid and fluid, fluid transients, pumps, compressors and external body forces such as gravity and centrifugal. The thermofluid system to be analyzed is discretized into nodes, branches, and conductors. The scalar properties such as pressure, temperature, and concentrations are calculated at nodes. Mass flow rates and heat transfer rates are computed in branches and conductors. The graphical user interface allows users to build their models using the point, drag and click method; the users can also run their models and post-process the results in the same environment. The integrated fluid library supplies thermodynamic and thermo-physical properties of 36 fluids and 21 different resistance/source options are provided for modeling momentum sources or sinks in the branches. This Technical Memorandum illustrates the application and verification of the code through 12 demonstrated example problems.

Numerical simulation of flow-induced vibrations in tube bundles

International Nuclear Information System (INIS)

Elisabeth Longatte; Zaky Bendjeddou; Mhamed Souli

2005-01-01

Full text of publication follows: In many industrial components mechanical structures like rod cluster control assembly, fuel assembly and heat exchanger tube bundles are submitted to complex flows causing possible vibrations and damage. Fluid forces are usually split into two parts: structure motion independent forces and fluid-elastic forces coupled with tube motion and responsible for possible dynamic instability development leading to possible short term failures through high amplitude vibrations. Most classical fluid force identification methods rely on structure response experimental measurements associated with convenient data processes. Owing to recent improvements in Computational Fluid Dynamics (C.F.D.), numerical fluid force identification is now practicable in the presence of industrial configurations. The present paper is devoted to numerical simulation of flow-induced vibrations of tube bundles submitted to single-phase cross flows by using C.F.D. codes. Direct Numerical Simulation (D.N.S.), Arbitrary Lagrange Euler formulation (A.L.E.) and code coupling process are involved to predict fluid forces responsible for tube bundle vibrations in the presence of fluid structure and fluid-elastic coupling effects. In the presence of strong multi-physics coupling, simulation of flow-induced vibrations requires a fluid structure code coupling process. The methodology consists in solving in the same time thermohydraulics and mechanics problems by using an A.L.E. formulation for the fluid computation. The purpose is to take into account coupling between flow and structure motions in order to be able to capture coupling effects. From a numerical point of view, there are three steps in the computation: the fluid problem is solved on the computational domain; fluid forces acting on the moving tube are estimated; finally they are introduced in the structure solver providing the tube displacement that is used to actualize the fluid computational domain. Specific

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.

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.

Computational fluid dynamics analysis of a mixed flow pump impeller

African Journals Online (AJOL)

ATHARVA

International Journal of Engineering, Science and Technology ... From the CFD analysis software and advanced post processing tools the complex flow inside the ... The numerical simulation can provide quite accurate information on the fluid ...

Simulation and Verificaiton of Flow in Test Methods

DEFF Research Database (Denmark)

Thrane, Lars Nyholm; Szabo, Peter; Geiker, Mette Rica

2005-01-01

Simulations and experimental results of L-box and slump flow test of a self-compacting mortar and a self-compacting concrete are compared. The simulations are based on a single fluid approach and assume an ideal Bingham behavior. It is possible to simulate the experimental results of both tests...

Pore scale simulations for the extension of the Darcy-Forchheimer law to shear thinning fluids

Science.gov (United States)

Tosco, Tiziana; Marchisio, Daniele; Lince, Federica; Boccardo, Gianluca; Sethi, Rajandrea

2014-05-01

Flow of non-Newtonian fluids through porous media at high Reynolds numbers is often encountered in chemical, pharmaceutical and food as well as petroleum and groundwater engineering and in many other industrial applications (1 - 2). In particular, the use of shear thinning polymeric solutions has been recently proposed to improve colloidal stability of micro- and nanoscale zerovalent iron particles (MZVI and NZVI) for groundwater remediation. In all abovementioned applications, it is of paramount importance to correctly predict the pressure drop resulting from non-Newtonian fluid flow through the porous medium. For small Reynolds numbers, usually up to 1, typical of laboratory column tests, the extended Darcy law is known to be applicable also to non Newtonian fluids, provided that all non-Newtonian effects are lumped together into a proper viscosity parameter (1,3). For higher Reynolds numbers (eg. close to the injection wells) non linearities between pressure drop and flow rate arise, and the Darcy-Forchheimer law holds for Newtonian fluids, while for non-Newtonian fluids, it has been demonstrated that, at least for simple rheological models (eg. power law fluids) a generalized Forchheimer law can be applied, even if the determination of the flow parameters (permeability K, inertial coefficient β, and equivalent viscosity) is not straightforward. This work (co-funded by European Union project AQUAREHAB FP7 - Grant Agreement Nr. 226565) aims at proposing an extended formulation of the Darcy-Forchheimer law also for shear-thinning fluids, and validating it against results of pore-scale simulations via computational fluid dynamics (4). Flow simulations were performed using Fluent 12.0 on four different 2D porous domains for Newtonian and non-Newtonian fluids (Cross, Ellis and Carreau models). The micro-scale flow simulation results are analyzed in terms of 'macroscale' pressure drop between inlet and outlet of the model domain as a function of flow rate. The

Leveraging Understanding of Flow of Variable Complex Fluid to Design Better Absorbent Hygiene Products

Science.gov (United States)

Krautkramer, C.; Rend, R. R.

2014-12-01

Menstrual flow, which is a result of shedding of uterus endometrium, occurs periodically in sync with a women's hormonal cycle. Management of this flow while allowing women to pursue their normal daily lives is the purpose of many commercial products. Some of these products, e.g. feminine hygiene pads and tampons, utilize porous materials in achieving their goal. In this paper we will demonstrate different phenomena that have been observed in flow of menstrual fluid through these porous materials, share some of the advances made in experimental and analytical study of these phenomena, and also present some of the unsolved challenges and difficulties encountered while studying this kind of flow. Menstrual fluid is generally composed of four main components: blood plasma, blood cells, cervical mucus, and tissue debris. This non-homogeneous, multiphase fluid displays very complex rheological behavior, e. g., yield stress, thixotropy, and visco-elasticity, that varies throughout and between menstrual cycles and among women due to various factors. Flow rates are also highly variable during menstruation and across the population and the rheological properties of the fluid change during the flow into and through the product. In addition to these phenomena, changes to the structure of the porous medium within the product can also be seen due to fouling and/or swelling of the material. This paper will, also, share how the fluid components impact the flow and the consequences for computer simulation, the creation of a simulant fluid and testing methods, and for designing products that best meet consumer needs. We hope to bring to light the challenges of managing this complex flow to meet a basic need of women all over the world. An opportunity exists to apply learnings from research in other disciplines to improve the scientific knowledge related to the flow of this complex fluid through the porous medium that is a sanitary product.

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.

Simulation of the flow inside an annular can combustor

OpenAIRE

Alqaraghuli, W; Alkhafagiy, D; Shires, A

2014-01-01

In the gas turbine combustion system, the external flows in annuli play one of the key roles in controlling pressure loss, air flow distribution around the combustor liner, and the attendant effects on performance, durability, and stability.  This paper describes a computational fluid dynamics (CFD) simulation of the flow in the outer annulus of a can combustor. Validating this simulation was done with experimental results obtained from analyzing the flow inside a can combustor annulus that w...

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.

Numerical Simulation Of Flow Through An Artificial Heart

Science.gov (United States)

Rogers, Stuart; Kutler, Paul; Kwak, Dochan; Kiris, Centin

1991-01-01

Research in both artificial hearts and fluid dynamics benefits from computational studies. Algorithm that implements Navier-Stokes equations of flow extended to simulate flow of viscous, incompressible blood through articifial heart. Ability to compute details of such flow important for two reasons: internal flows with moving boundaries of academic interest in their own right, and many of deficiencies of artificial hearts attributable to dynamics of flow.

Fluid dynamics of cryogenic two-phase flows

International Nuclear Information System (INIS)

Verfondern, K.; Jahn, W.

2004-01-01

The objective of this study was to examine the flow behavior of a methane hydrate/methane-liquid hydrogen dispersed two-phase fluid through a given design of a moderator chamber for the ESS target system. The calculations under simplified conditions, e.g., taking no account of heat input from outside, have shown that the computer code used, CFX, was able to simulate the behavior of the two-phase flow through the moderator chamber, producing reasonable results up to a certain level of the solid phase fraction, that allowed a continuous flow process through the chamber. Inlet flows with larger solid phase fractions than 40 vol% were found to be a ''problem'' for the computer code. From the computer runs based on fractions between 20 and 40 vol%, it was observed that with increasing solid phase fraction at the inlet, the resulting flow pattern revealed a strong tendency for blockage within the chamber, supported by the ''heavy weight'' of the pellets compared to the carrying liquid. Locations which are prone to the development of such uneven flow behavior are the areas around the turning points in the semispheres and near the exit of the moderator. The considered moderator chamber with horizontal inlet and outlet flow for a solid-liquid two-phase fluid does not seem to be an appropriate design. (orig.)

Two-dimensional simulation of clastic and carbonate sedimentation, consolidation, subsidence, fluid flow, heat flow and solute transport during the formation of sedimentary basins

Science.gov (United States)

Bitzer, Klaus

1999-05-01

Geological processes that create sedimentary basins or act during their formation can be simulated using the public domain computer code `BASIN'. For a given set of geological initial and boundary conditions the sedimentary basin evolution is calculated in a forward modeling approach. The basin is represented in a two-dimensional vertical cross section with individual layers. The stratigraphic, tectonic, hydrodynamic and thermal evolution is calculated beginning at an initial state, and subsequent changes of basin geometry are calculated from sedimentation rates, compaction and pore fluid mobilization, isostatic compensation, fault movement and subsidence. The sedimentologic, hydraulic and thermal parameters are stored at discrete time steps allowing the temporal evolution of the basin to be analyzed. A maximum flexibility in terms of geological conditions is achieved by using individual program modules representing geological processes which can be switched on and off depending on the data available for a specific simulation experiment. The code incorporates a module for clastic and carbonate sedimentation, taking into account the impact of clastic sediment supply on carbonate production. A maximum of four different sediment types, which may be mixed during sedimentation, can be defined. Compaction and fluid flow are coupled through the consolidation equation and the nonlinear form of the equation of state for porosity, allowing nonequilibrium compaction and overpressuring to be calculated. Instead of empirical porosity-effective stress equations, a physically consistent consolidation model is applied which incorporates a porosity dependent sediment compressibility. Transient solute transport and heat flow are calculated as well, applying calculated fluid flow rates from the hydraulic model. As a measure for hydrocarbon generation, the Time-Temperature Index (TTI) is calculated. Three postprocessing programs are available to provide graphic output in Post

Computational Fluid Dynamics and Building Energy Performance Simulation

DEFF Research Database (Denmark)

Nielsen, Peter Vilhelm; Tryggvason, T.

1998-01-01

An interconnection between a building energy performance simulation program and a Computational Fluid Dynamics program (CFD) for room air distribution will be introduced for improvement of the predictions of both the energy consumption and the indoor environment. The building energy performance...... simulation program requires a detailed description of the energy flow in the air movement which can be obtained by a CFD program. The paper describes an energy consumption calculation in a large building, where the building energy simulation program is modified by CFD predictions of the flow between three...... zones connected by open areas with pressure and buoyancy driven air flow. The two programs are interconnected in an iterative procedure. The paper shows also an evaluation of the air quality in the main area of the buildings based on CFD predictions. It is shown that an interconnection between a CFD...

A numerical model for dynamic crustal-scale fluid flow

Science.gov (United States)

Sachau, Till; Bons, Paul; Gomez-Rivas, Enrique; Koehn, Daniel

2015-04-01

Fluid flow in the crust is often envisaged and modeled as continuous, yet minimal flow, which occurs over large geological times. This is a suitable approximation for flow as long as it is solely controlled by the matrix permeability of rocks, which in turn is controlled by viscous compaction of the pore space. However, strong evidence (hydrothermal veins and ore deposits) exists that a significant part of fluid flow in the crust occurs strongly localized in both space and time, controlled by the opening and sealing of hydrofractures. We developed, tested and applied a novel computer code, which considers this dynamic behavior and couples it with steady, Darcian flow controlled by the matrix permeability. In this dual-porosity model, fractures open depending on the fluid pressure relative to the solid pressure. Fractures form when matrix permeability is insufficient to accommodate fluid flow resulting from compaction, decompression (Staude et al. 2009) or metamorphic dehydration reactions (Weisheit et al. 2013). Open fractures can close when the contained fluid either seeps into the matrix or escapes by fracture propagation: mobile hydrofractures (Bons, 2001). In the model, closing and sealing of fractures is controlled by a time-dependent viscous law, which is based on the effective stress and on either Newtonian or non-Newtonian viscosity. Our simulations indicate that the bulk of crustal fluid flow in the middle to lower upper crust is intermittent, highly self-organized, and occurs as mobile hydrofractures. This is due to the low matrix porosity and permeability, combined with a low matrix viscosity and, hence, fast sealing of fractures. Stable fracture networks, generated by fluid overpressure, are restricted to the uppermost crust. Semi-stable fracture networks can develop in an intermediate zone, if a critical overpressure is reached. Flow rates in mobile hydrofractures exceed those in the matrix porosity and fracture networks by orders of magnitude

Lattice fluid dynamics from perfect discretizations of continuum flows

International Nuclear Information System (INIS)

Katz, E.; Wiese, U.

1998-01-01

We use renormalization group methods to derive equations of motion for large scale variables in fluid dynamics. The large scale variables are averages of the underlying continuum variables over cubic volumes and naturally exist on a lattice. The resulting lattice dynamics represents a perfect discretization of continuum physics, i.e., grid artifacts are completely eliminated. Perfect equations of motion are derived for static, slow flows of incompressible, viscous fluids. For Hagen-Poiseuille flow in a channel with a square cross section the equations reduce to a perfect discretization of the Poisson equation for the velocity field with Dirichlet boundary conditions. The perfect large scale Poisson equation is used in a numerical simulation and is shown to represent the continuum flow exactly. For nonsquare cross sections one can use a numerical iterative procedure to derive flow equations that are approximately perfect. copyright 1998 The American Physical Society

A Level-set based framework for viscous simulation of particle-laden supersonic flows

Science.gov (United States)

Das, Pratik; Sen, Oishik; Jacobs, Gustaaf; Udaykumar, H. S.

2017-06-01

Particle-laden supersonic flows are important in natural and industrial processes, such as, volcanic eruptions, explosions, pneumatic conveyance of particle in material processing etc. Numerical study of such high-speed particle laden flows at the mesoscale calls for a numerical framework which allows simulation of supersonic flow around multiple moving solid objects. Only a few efforts have been made toward development of numerical frameworks for viscous simulation of particle-fluid interaction in supersonic flow regime. The current work presents a Cartesian grid based sharp-interface method for viscous simulations of interaction between supersonic flow with moving rigid particles. The no-slip boundary condition is imposed at the solid-fluid interfaces using a modified ghost fluid method (GFM). The current method is validated against the similarity solution of compressible boundary layer over flat-plate and benchmark numerical solution for steady supersonic flow over cylinder. Further validation is carried out against benchmark numerical results for shock induced lift-off of a cylinder in a shock tube. 3D simulation of steady supersonic flow over sphere is performed to compare the numerically obtained drag co-efficient with experimental results. A particle-resolved viscous simulation of shock interaction with a cloud of particles is performed to demonstrate that the current method is suitable for large-scale particle resolved simulations of particle-laden supersonic flows.

Multiple stable isotope fronts during non-isothermal fluid flow

Science.gov (United States)

Fekete, Szandra; Weis, Philipp; Scott, Samuel; Driesner, Thomas

2018-02-01

Stable isotope signatures of oxygen, hydrogen and other elements in minerals from hydrothermal veins and metasomatized host rocks are widely used to investigate fluid sources and paths. Previous theoretical studies mostly focused on analyzing stable isotope fronts developing during single-phase, isothermal fluid flow. In this study, numerical simulations were performed to assess how temperature changes, transport phenomena, kinetic vs. equilibrium isotope exchange, and isotopic source signals determine mineral oxygen isotopic compositions during fluid-rock interaction. The simulations focus on one-dimensional scenarios, with non-isothermal single- and two-phase fluid flow, and include the effects of quartz precipitation and dissolution. If isotope exchange between fluid and mineral is fast, a previously unrecognized, significant enrichment in heavy oxygen isotopes of fluids and minerals occurs at the thermal front. The maximum enrichment depends on the initial isotopic composition of fluid and mineral, the fluid-rock ratio and the maximum change in temperature, but is independent of the isotopic composition of the incoming fluid. This thermally induced isotope front propagates faster than the signal related to the initial isotopic composition of the incoming fluid, which forms a trailing front behind the zone of transient heavy oxygen isotope enrichment. Temperature-dependent kinetic rates of isotope exchange between fluid and rock strongly influence the degree of enrichment at the thermal front. In systems where initial isotope values of fluids and rocks are far from equilibrium and isotope fractionation is controlled by kinetics, the temperature increase accelerates the approach of the fluid to equilibrium conditions with the host rock. Consequently, the increase at the thermal front can be less dominant and can even generate fluid values below the initial isotopic composition of the input fluid. As kinetics limit the degree of isotope exchange, a third front may

Simulation of two-phase flow in horizontal fracture networks with numerical manifold method

Science.gov (United States)

Ma, G. W.; Wang, H. D.; Fan, L. F.; Wang, B.

2017-10-01

The paper presents simulation of two-phase flow in discrete fracture networks with numerical manifold method (NMM). Each phase of fluids is considered to be confined within the assumed discrete interfaces in the present method. The homogeneous model is modified to approach the mixed fluids. A new mathematical cover formation for fracture intersection is proposed to satisfy the mass conservation. NMM simulations of two-phase flow in a single fracture, intersection, and fracture network are illustrated graphically and validated by the analytical method or the finite element method. Results show that the motion status of discrete interface significantly depends on the ratio of mobility of two fluids rather than the value of the mobility. The variation of fluid velocity in each fracture segment and the driven fluid content are also influenced by the ratio of mobility. The advantages of NMM in the simulation of two-phase flow in a fracture network are demonstrated in the present study, which can be further developed for practical engineering applications.

Benchmark exercise for fluid flow simulations in a liquid metal fast reactor fuel assembly

Energy Technology Data Exchange (ETDEWEB)

Merzari, E., E-mail: emerzari@anl.gov [Mathematics and Computer Science Division, Argonne National Laboratory, 9700 S. Cass Avenue, Lemont, IL 60439 (United States); Fischer, P. [Mathematics and Computer Science Division, Argonne National Laboratory, 9700 S. Cass Avenue, Lemont, IL 60439 (United States); Yuan, H. [Nuclear Engineering Division, Argonne National Laboratory, Lemont, IL (United States); Van Tichelen, K.; Keijers, S. [SCK-CEN, Boeretang 200, Mol (Belgium); De Ridder, J.; Degroote, J.; Vierendeels, J. [Ghent University, Ghent (Belgium); Doolaard, H.; Gopala, V.R.; Roelofs, F. [NRG, Petten (Netherlands)

2016-03-15

Highlights: • A EUROTAM-US INERI consortium has performed a benchmark exercise related to fast reactor assembly simulations. • LES calculations for a wire-wrapped rod bundle are compared with RANS calculations. • Results show good agreement for velocity and cross flows. - Abstract: As part of a U.S. Department of Energy International Nuclear Energy Research Initiative (I-NERI), Argonne National Laboratory (Argonne) is collaborating with the Dutch Nuclear Research and consultancy Group (NRG), the Belgian Nuclear Research Centre (SCK·CEN), and Ghent University (UGent) in Belgium to perform and compare a series of fuel-pin-bundle calculations representative of a fast reactor core. A wire-wrapped fuel bundle is a complex configuration for which little data is available for verification and validation of new simulation tools. UGent and NRG performed their simulations with commercially available computational fluid dynamics (CFD) codes. The high-fidelity Argonne large-eddy simulations were performed with Nek5000, used for CFD in the Simulation-based High-efficiency Advanced Reactor Prototyping (SHARP) suite. SHARP is a versatile tool that is being developed to model the core of a wide variety of reactor types under various scenarios. It is intended both to serve as a surrogate for physical experiments and to provide insight into experimental results. Comparison of the results obtained by the different participants with the reference Nek5000 results shows good agreement, especially for the cross-flow data. The comparison also helps highlight issues with current modeling approaches. The results of the study will be valuable in the design and licensing process of MYRRHA, a flexible fast research reactor under design at SCK·CEN that features wire-wrapped fuel bundles cooled by lead-bismuth eutectic.

Displacement of one Newtonian fluid by another: density effects in axial annular flow

DEFF Research Database (Denmark)

Szabo, Peter; Hassager, Ole

1997-01-01

The arbitrary Lagrange-Euler (ALE) finite elementtechnique is used to simulate 3D displacement oftwo immiscible Newtonian fluids in vertical annular wells. For equally viscous fluids the effect of distinct fluid densities is investigated in the region of low to intermediate Reynolds numbers......, the efficiency of the displacement is analysed for various flow situations....

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

High order methods for incompressible fluid flow: Application to moving boundary problems

Energy Technology Data Exchange (ETDEWEB)

Bjoentegaard, Tormod

2008-04-15

Fluid flows with moving boundaries are encountered in a large number of real life situations, with two such types being fluid-structure interaction and free-surface flows. Fluid-structure phenomena are for instance apparent in many hydrodynamic applications; wave effects on offshore structures, sloshing and fluid induced vibrations, and aeroelasticity; flutter and dynamic response. Free-surface flows can be considered as a special case of a fluid-fluid interaction where one of the fluids are practically inviscid, such as air. This type of flows arise in many disciplines such as marine hydrodynamics, chemical engineering, material processing, and geophysics. The driving forces for free-surface flows may be of large scale such as gravity or inertial forces, or forces due to surface tension which operate on a much smaller scale. Free-surface flows with surface tension as a driving mechanism include the flow of bubbles and droplets, and the evolution of capillary waves. In this work we consider incompressible fluid flow, which are governed by the incompressible Navier-Stokes equations. There are several challenges when simulating moving boundary problems numerically, and these include - Spatial discretization - Temporal discretization - Imposition of boundary conditions - Solution strategy for the linear equations. These are some of the issues which will be addressed in this introduction. We will first formulate the problem in the arbitrary Lagrangian-Eulerian framework, and introduce the weak formulation of the problem. Next, we discuss the spatial and temporal discretization before we move to the imposition of surface tension boundary conditions. In the final section we discuss the solution of the resulting linear system of equations. (Author). refs., figs., tabs

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.

Direct numerical simulation of particulate flow with heat transfer

NARCIS (Netherlands)

Tavassoli Estahbanati, H; Kriebitzsch, S.H.L.; Hoef, van der M.A.; Peters, E.A.J.F.; Kuipers, J.A.M.

2013-01-01

The Immersed Boundary (IB) method proposed by Uhlmann for Direct Numerical Simulation (DNS) of fluid flow through dense fluid-particle systems is extended to systems with interphase heat transport. A fixed Eulerian grid is employed to solve the momentum and energy equations by traditional

Production of Natural Gas and Fluid Flow in Tight Sand Reservoirs

Energy Technology Data Exchange (ETDEWEB)

Maria Cecilia Bravo

2006-06-30

This document reports progress of this research effort in identifying relationships and defining dependencies between macroscopic reservoir parameters strongly affected by microscopic flow dynamics and production well performance in tight gas sand reservoirs. These dependencies are investigated by identifying the main transport mechanisms at the pore scale that should affect fluids flow at the reservoir scale. A critical review of commercial reservoir simulators, used to predict tight sand gas reservoir, revealed that many are poor when used to model fluid flow through tight reservoirs. Conventional simulators ignore altogether or model incorrectly certain phenomena such as, Knudsen diffusion, electro-kinetic effects, ordinary diffusion mechanisms and water vaporization. We studied the effect of Knudsen's number in Klinkenberg's equation and evaluated the effect of different flow regimes on Klinkenberg's parameter b. We developed a model capable of explaining the pressure dependence of this parameter that has been experimentally observed, but not explained in the conventional formalisms. We demonstrated the relevance of this, so far ignored effect, in tight sands reservoir modeling. A 2-D numerical simulator based on equations that capture the above mentioned phenomena was developed. Dynamic implications of new equations are comprehensively discussed in our work and their relative contribution to the flow rate is evaluated. We performed several simulation sensitivity studies that evidenced that, in general terms, our formalism should be implemented in order to get more reliable tight sands gas reservoirs' predictions.

Numerical simulation of two-phase flow with front-capturing

International Nuclear Information System (INIS)

Tzanos, C.P.; Weber, D.P.

2000-01-01

Because of the complexity of two-phase flow phenomena, two-phase flow codes rely heavily on empirical correlations. This approach has a number of serious shortcomings. Advances in parallel computing and continuing improvements in computer speed and memory have stimulated the development of numerical simulation tools that rely less on empirical correlations and more on fundamental physics. The objective of this work is to take advantage of developments in massively parallel computing, single-phase computational fluid dynamics of complex systems, and numerical methods for front capturing in two-phase flows to develop a computer code for direct numerical simulation of two-phase flow. This includes bubble/droplet transport, interface deformation and topology change, bubble-droplet interactions, interface mass, momentum, and energy transfer. In this work, the Navier-Stokes and energy equations are solved by treating both phases as a single fluid with interfaces between the two phases, and a discontinuity in material properties across the moving interfaces. The evolution of the interfaces is simulated by using the front capturing technique of the level-set methods. In these methods, the boundary of a two-fluid interface is modeled as the zero level set of a smooth function φ. The level-set function φ is defined as the signed distance from the interface (φ is negative inside a droplet/bubble and positive outside). Compared to other front-capturing or front-tracking methods, the level-set approach is relatively easy to implement even in three-dimensional flows, and it has been shown to simulate well the coalescence and breakup of droplets/bubbles

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.

Simulating the fluid-structure interaction of a flexible tube in an array of rigid tubes

International Nuclear Information System (INIS)

Warnica, D.; Maleki, M.; Hariri, A.; Feldman, H.

2011-01-01

Two important single-phase mechanisms for flow-induced vibration of heat-exchanger tube bundles were used to demonstrate the capabilities of commercial software to simulate unsteady fluid-structure interactions (FSI). Reasonable agreement was obtained between the FSI simulations and experimental data for the onset of fluid elastic instability. There was also reasonable agreement between the FSI simulations and empirical correlations for the dynamic tube response to random turbulence excitation. Additional benefits of performing FSI simulations were the ability to characterize important features of the unsteady flow fields and hydrodynamic parameters such as viscous damping coefficients, which would otherwise require elaborate experimental measurements. (author)

Stochastic Rotation Dynamics simulations of wetting multi-phase flows

Science.gov (United States)

Hiller, Thomas; Sanchez de La Lama, Marta; Brinkmann, Martin

2016-06-01

Multi-color Stochastic Rotation Dynamics (SRDmc) has been introduced by Inoue et al. [1,2] as a particle based simulation method to study the flow of emulsion droplets in non-wetting microchannels. In this work, we extend the multi-color method to also account for different wetting conditions. This is achieved by assigning the color information not only to fluid particles but also to virtual wall particles that are required to enforce proper no-slip boundary conditions. To extend the scope of the original SRDmc algorithm to e.g. immiscible two-phase flow with viscosity contrast we implement an angular momentum conserving scheme (SRD+mc). We perform extensive benchmark simulations to show that a mono-phase SRDmc fluid exhibits bulk properties identical to a standard SRD fluid and that SRDmc fluids are applicable to a wide range of immiscible two-phase flows. To quantify the adhesion of a SRD+mc fluid in contact to the walls we measure the apparent contact angle from sessile droplets in mechanical equilibrium. For a further verification of our wettability implementation we compare the dewetting of a liquid film from a wetting stripe to experimental and numerical studies of interfacial morphologies on chemically structured surfaces.

Grid studies for the simulation of resolved structures in an Eulerian two-fluid framework

Energy Technology Data Exchange (ETDEWEB)

Gauss, Friederike, E-mail: f.gauss@hzdr.de; Lucas, Dirk; Krepper, Eckhard

2016-08-15

Highlights: • Elaborated Eulerian two-fluid methods may predict multiphase flow with large differences in interfacial length scales. • A study on the grid requirements of resolved structures in such two-fluid methods is presented. • The two-fluid results are only little dependent on the grid size. • The results justify the resolved treatment of flow structures covering only few grid cells. • A grid-dependent limit between resolved an modeled structures may be established. - Abstract: The influence of the grid size on the rise velocity of a single bubble simulated with an Eulerian two-fluid method is investigated. This study is part of the development of an elaborated Eulerian two-fluid framework, which is able to predict complex flow phenomena as arising in nuclear reactor safety research issues. Such flow phenomena cover a wide range of interfacial length scales. An important aspect of the simulation method is the distinction into small flow structures, which are modeled, and large structures, which are resolved. To investigate the requirements on the numerical grid for the simulation of such resolved structures the velocity of rising gas bubbles is a good example since theoretical values are available. It is well known that the rise velocity of resolved bubbles is clearly underestimated in a one-fluid approach if they span over only few numerical cells. In the present paper it is shown that in the case of the two-fluid model the bubble rise velocity depends only slightly on the grid size. This is explained with the use of models for the gas–liquid interfacial forces. Good approximations of the rise velocity and the bubble shape are obtained with only few grid points per bubble diameter. This result justifies the resolved treatment of flow structures, which cover only few grid cells. Thus, a limit for the distinction into resolved and modeled structures in the two-fluid context may be established.

Direct numerical simulations of gas-liquid multiphase flows

CERN Document Server

Tryggvason, Grétar; Zaleski, Stéphane

2011-01-01

Accurately predicting the behaviour of multiphase flows is a problem of immense industrial and scientific interest. Modern computers can now study the dynamics in great detail and these simulations yield unprecedented insight. This book provides a comprehensive introduction to direct numerical simulations of multiphase flows for researchers and graduate students. After a brief overview of the context and history the authors review the governing equations. A particular emphasis is placed on the 'one-fluid' formulation where a single set of equations is used to describe the entire flow field and

Computational investigation of fluid flow and heat transfer of an economizer by porous medium approach

Science.gov (United States)

Babu, C. Rajesh; Kumar, P.; Rajamohan, G.

2017-07-01

Computation of fluid flow and heat transfer in an economizer is simulated by a porous medium approach, with plain tubes having a horizontal in-line arrangement and cross flow arrangement in a coal-fired thermal power plant. The economizer is a thermal mechanical device that captures waste heat from the thermal exhaust flue gasses through heat transfer surfaces to preheat boiler feed water. In order to evaluate the fluid flow and heat transfer on tubes, a numerical analysis on heat transfer performance is carried out on an 110 t/h MCR (Maximum continuous rating) boiler unit. In this study, thermal performance is investigated using the computational fluid dynamics (CFD) simulation using ANSYS FLUENT. The fouling factor ε and the overall heat transfer coefficient ψ are employed to evaluate the fluid flow and heat transfer. The model demands significant computational details for geometric modeling, grid generation, and numerical calculations to evaluate the thermal performance of an economizer. The simulation results show that the overall heat transfer coefficient 37.76 W/(m2K) and economizer coil side pressure drop of 0.2 (kg/cm2) are found to be conformity within the tolerable limits when compared with existing industrial economizer data.

Simulation of ferromagnetic nanomaterial flow of Maxwell fluid

Directory of Open Access Journals (Sweden)

T. Hayat

2018-03-01

Full Text Available Ferromagnetic flow of rate type liquid over a stretched surface is addressed in this article. Heat and mass transport are investigated with Brownian movement and thermophoresis effects. Magnetic dipole is also taken into consideration. Procedure of similarity transformation is employed. The obtained nonlinear expressions have been tackled numerically by means of Shooting method. Graphical results are shown and analyzed for the impact of different variables. Temperature and concentration gradients are numerically computed in Tables 1 and 2. The results described here demonstrate that ferromagnetic variable boosts the thermal field. It is noticed that velocity and concentration profiles are higher when elastic and thermophoresis variables are enhanced. Keywords: Rate type fluid, Brownian movement, Thermophoresis effect, Magnetic dipole

Fluid flows due to earthquakes with reference to Yucca Mountain, Nevada

International Nuclear Information System (INIS)

Davies, J.B.

1993-01-01

Yucca Mountain geohydrology is dominated by a deep water table in volcanic tuffa beds which are cut by numerous faults. Certain zones in these tuffas and most of the fault apertures are filled with a fine-grained calcitic cement. Earthquakes have occured in this region with the most recent being of magnitude 5.6 and at a distance of about 20 km. Earthquakes in western U.S.A. have been observed to cause fluid flows through and out of the crust of the Earth. These flows are concentrated along the faults with normal faulting producing the largest flows. An earthquake produces rapid pressure changes at and below the ground surface, thereby forcing flows of gas, water, slurries and dissolved salts. In order to examine the properties of flows produced by earthquakes, we simulate the phenomena using computer-based modeling. We investigate the effects of adults and high permeability zones on the pattern of flows induced by the earthquake. We demonstrate that faults act as conduits to the surface and that the higher the permeability of a zone, the more the flows will concentrate there. Numerical estimates of flow rates from these simulations compare favorably with data from observed flows due to earthquakes. Simple volumetric arguments demonstrate the ease with which fluids from the deep water table can reach the surface along fault conduits

Numerical simulations of seepage flow in rough single rock fractures

Directory of Open Access Journals (Sweden)

Qingang Zhang

2015-09-01

Full Text Available To investigate the relationship between the structural characteristics and seepage flow behavior of rough single rock fractures, a set of single fracture physical models were produced using the Weierstrass–Mandelbrot functions to test the seepage flow performance. Six single fractures, with various surface roughnesses characterized by fractal dimensions, were built using COMSOL multiphysics software. The fluid flow behavior through the rough fractures and the influences of the rough surfaces on the fluid flow behavior was then monitored. The numerical simulation indicates that there is a linear relationship between the average flow velocity over the entire flow path and the fractal dimension of the rough surface. It is shown that there is good a agreement between the numerical results and the experimental data in terms of the properties of the fluid flowing through the rough single rock fractures.

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)

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.

Fluid Flow through Porous Sandstone with Overprinting and Intersecting Geological Structures of Various Types

Science.gov (United States)

Zhou, X.; Karimi-Fard, M.; Durlofsky, L.; Aydin, A.

2010-12-01

Impact of a wide variety of structural heterogeneities on fluid flow in an aeolian sandstone in the Valley of Fire State Park (NV), such as (1) dilatant fractures (joints), (2) shear fractures (faults), and (3) contraction/compaction structures (compaction bands), are considered. Each type of these structures has its own geometry, spacing, distribution, connectivity, and hydraulic properties, which either enhance or impede subsurface fluid flow. Permeability of these structures may, on average, be a few orders of magnitude higher or lower than those of the corresponding matrix rocks. In recent years, the influence of a single type of these heterogeneities on fluid flow has been studied individually, such as joints, compaction bands or faults. However, as different types of geological structures are commonly present together in the same rock volume, their combined effect requires a more detailed assessment. In this study, fluid flow simulations are performed using a special finite-volume discretization technique that was developed by Karimi-Fard et al. (2004; 2006). Using this approach, thin features such as fractures and compaction bands are represented as linear elements in unstructured 2D models and as planar elements in 3D models, which significantly reduces the total number of cells and simplifies grid generation. The cell geometric information and the cell-to-cell transmissibility obtained from this discretization technique are input to Stanford’s General Purpose Research Simulator (GPRS) for fluid flow simulation. To account for the effects of the various geological structures on subsurface flow, we perform permeability upscaling over regions corresponding to large-scale simulation grid blocks in order to obtain equivalent permeability components in two principal directions. We will focus on the following problems: (1) compaction bands of multisets; (2) compartmentalization of compaction bands of high-angle, low-angle and horizontal; (3) joints overprinting

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.

Atomistic Simulations of Fluid Flow through Graphene Channels and Carbon Nanotubes

DEFF Research Database (Denmark)

Zambrano, Harvey A.; Walther, Jens Honore; Oyarzua, Elton E.

2015-01-01

conductivity, extremely low surface friction and superior mechanical properties, graphene channels and carbon nanotubes (CNTs) are promising candidates to be implemented as fluid conduits in nanosystems. Performing Non-equilibrium Molecular Dynamics simulations, we study the transport of water......-eletrolyte solutions inside single and multi-wall graphene channels and inside zig-zag and armchair CNTs of similar cross sectional area. In order to calibrate the force fields, we use dedicated criteria relevant to the hydrodynamics of the systems of interest. Different fluid driving mechanisms such as pressure...

Development and computational simulation of thermoelectric electromagnetic pumps for controlling the fluid flow in liquid metal cooled space nuclear reactors

International Nuclear Information System (INIS)

Borges, E.M.

1991-01-01

Thermoelectric Electromagnetic (TEEM) Pumps can be used for controlling the fluid flow in the primary and secondary circuits of liquid metal cooled space nuclear reactor. In order to simulate and to evaluate the pumps performance, in steady-state, the computer program BEMTE has been developed to study the main operational parameters and to determine the system actuation point, for a given reactor operating power. The results for each stage of the program were satisfactory, compared to experimental data. The program shows to be adequate for the design and simulating of direct current electromagnetic pumps. (author)

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

Characterization of cardiac flow in heart disease patients by computational fluid dynamics and 4D flow MRI

Science.gov (United States)

Lantz, Jonas; Gupta, Vikas; Henriksson, Lilian; Karlsson, Matts; Persson, Ander; Carhall, Carljohan; Ebbers, Tino

2017-11-01

In this study, cardiac blood flow was simulated using Computational Fluid Dynamics and compared to in vivo flow measurements by 4D Flow MRI. In total, nine patients with various heart diseases were studied. Geometry and heart wall motion for the simulations were obtained from clinical CT measurements, with 0.3x0.3x0.3 mm spatial resolution and 20 time frames covering one heartbeat. The CFD simulations included pulmonary veins, left atrium and ventricle, mitral and aortic valve, and ascending aorta. Mesh sizes were on the order of 6-16 million cells, depending on the size of the heart, in order to resolve both papillary muscles and trabeculae. The computed flow field agreed visually very well with 4D Flow MRI, with characteristic vortices and flow structures seen in both techniques. Regression analysis showed that peak flow rate as well as stroke volume had an excellent agreement for the two techniques. We demonstrated the feasibility, and more importantly, fidelity of cardiac flow simulations by comparing CFD results to in vivo measurements. Both qualitative and quantitative results agreed well with the 4D Flow MRI measurements. Also, the developed simulation methodology enables ``what if'' scenarios, such as optimization of valve replacement and other surgical procedures. Funded by the Wallenberg Foundation.

Simulation of transient fluid flow in mold region during steel continuous casting

International Nuclear Information System (INIS)

Liu, R; Thomas, B G; Sengupta, J

2012-01-01

A system of models has been developed to study transient flow during continuous casting and applied to simulate an event of multiple stopper-rod movements. It includes four sub-models to incorporate different aspects in this transient event. A three-dimensional (3-D) porous-flow model of the nozzle wall calculates the rate argon gas flow into the liquid steel, and the initial mean bubble size is estimated. Transient CFD models simulate multiphase flow of steel and gas bubbles in the Submerged Entry Nozzle (SEN) and mold and have been validated with experimental data from both nail dipping and Sub-meniscus Velocity Control (SVC) measurements. To obtain the transient inlet boundary conditions for the simulation, two semi-empirical models, a stopper-rod-position based model and a metal-level-based model, predict the liquid steel flow rate through the SEN based on recorded plant data. Finally the model system was applied to study the effects of stopper rod movements on SEN/mold flow patterns. Meniscus level fluctuations were calculated using a simple pressure method and compared well with plant measurements. Insights were gained from the simulation results to explain the cause of meniscus level fluctuations and the formation of sliver defects during stopper rod movements.

Towards numerical simulations of supersonic liquid jets using ghost fluid method

International Nuclear Information System (INIS)

Majidi, Sahand; Afshari, Asghar

2015-01-01

Highlights: • A ghost fluid method based solver is developed for numerical simulation of compressible multiphase flows. • The performance of the numerical tool is validated via several benchmark problems. • Emergence of supersonic liquid jets in quiescent gaseous environment is simulated using ghost fluid method for the first time. • Bow-shock formation ahead of the liquid jet is clearly observed in the obtained numerical results. • Radiation of mach waves from the phase-interface witnessed experimentally is evidently captured in our numerical simulations. - Abstract: A computational tool based on the ghost fluid method (GFM) is developed to study supersonic liquid jets involving strong shocks and contact discontinuities with high density ratios. The solver utilizes constrained reinitialization method and is capable of switching between the exact and approximate Riemann solvers to increase the robustness. The numerical methodology is validated through several benchmark test problems; these include one-dimensional multiphase shock tube problem, shock–bubble interaction, air cavity collapse in water, and underwater-explosion. A comparison between our results and numerical and experimental observations indicate that the developed solver performs well investigating these problems. The code is then used to simulate the emergence of a supersonic liquid jet into a quiescent gaseous medium, which is the very first time to be studied by a ghost fluid method. The results of simulations are in good agreement with the experimental investigations. Also some of the famous flow characteristics, like the propagation of pressure-waves from the liquid jet interface and dependence of the Mach cone structure on the inlet Mach number, are reproduced numerically. The numerical simulations conducted here suggest that the ghost fluid method is an affordable and reliable scheme to study complicated interfacial evolutions in complex multiphase systems such as supersonic liquid

Molecular dynamics simulations of oscillatory flows in microfluidic channels

DEFF Research Database (Denmark)

Hansen, J.S.; Ottesen, Johnny T.

2006-01-01

In this paper we apply the direct non-equilibrium molecular dynamics technique to oscillatory flows of fluids in microscopic channels. Initially, we show that the microscopic simulations resemble the macroscopic predictions based on the Navier–Stokes equation very well for large channel width, high...... density and low temperature. Further simulations for high temperature and low density show that the non-slip boundary condition traditionally used in the macroscopic equation is greatly compromised when the fluid–wall interactions are the same as the fluid–fluid interactions. Simulations of a system...

Numerical simulation of transient, adiabatic, two-dimensional two-phase flow using the two-fluid model

International Nuclear Information System (INIS)

Neves Conti, T. das.

1983-01-01

A numerical method is developed to simulate adiabatic, transient, two-dimensional two-phase flow. The two-fluid model is used to obtain the mass and momentum conservation equations. These are solved by an iterative algorithm emphoying a time-marching scheme. Based on the corrective procedure of Hirt and Harlow a poisson equation is derived for the pressure field. This equation is finite-differenced and solved by a suitable matrix inversion technique. In the absence of experiment results several numerical tests were made in order to chec accuracy, convergence and stability of the proposed method. Several tests were also performed to check whether the behavior of void fraction and phasic velocities conforms with previous observations. (Author) [pt

Mass and momentum conservation for fluid simulation

KAUST Repository

Lentine, Michael; Aanjaneya, Mridul; Fedkiw, Ronald

2011-01-01

Momentum conservation has long been used as a design principle for solid simulation (e.g. collisions between rigid bodies, mass-spring elastic and damping forces, etc.), yet it has not been widely used for fluid simulation. In fact, semi-Lagrangian advection does not conserve momentum, but is still regularly used as a bread and butter method for fluid simulation. In this paper, we propose a modification to the semi-Lagrangian method in order to make it fully conserve momentum. While methods of this type have been proposed earlier in the computational physics literature, they are not necessarily appropriate for coarse grids, large time steps or inviscid flows, all of which are common in graphics applications. In addition, we show that the commonly used vorticity confinement turbulence model can be modified to exactly conserve momentum as well. We provide a number of examples that illustrate the benefits of this new approach, both in conserving fluid momentum and passively advected scalars such as smoke density. In particular, we show that our new method is amenable to efficient smoke simulation with one time step per frame, whereas the traditional non-conservative semi-Lagrangian method experiences serious artifacts when run with these large time steps, especially when object interaction is considered. Copyright © 2011 by the Association for Computing Machinery, Inc.

Hydrologic mechanisms governing fluid flow in partially saturated, fractured, porous tuff at Yucca Mountain

International Nuclear Information System (INIS)

Wang, J.S.Y.; Narasimhan, T.N.

1984-10-01

In contrast to the saturated zone where fluid moves rapidly along fractures, the fractures (with apertures large relative to the size of matrix pores) will desaturate first during drainage process and the bulk of fluid flow would be through interconnected pores in the matrix. Within a partially drained fracture, the presence of a relatively continuous air phase will produce practically an infinite resistance to liquid flow in the direction parallel to the fracture. The residual liquid will be held by capillary force in regions around fracture contact areas where the apertures are small. Normal to the fracture surfaces, the drained portion of the fractures will reduce the effective area for liquid flow from one matrix block to another matrix block. A general statistical theory is constructed for flow along the fracture and for flow between the matrix blocks to the fractures under partially saturated conditions. Results are obtained from an aperture distribution model for fracture saturation, hydraulic conductivity, and effective matrix-fracture flow areas as functions of pressure. Drainage from a fractured tuff column is simulated. The parameters for the simulations are deduced from fracture surface characteristics, spacings and orientations based on core analyses, and from matrix characteristics curve based on laboratory measurements. From the cases simulated for the fractured, porous column with discrete vertical and horizontal fractures and porous matrix blocks explicitly taken into account, it is observed that the highly transient changes from fully saturated conditions to partially saturated conditions are extremely sensitive to the fracture properties. However, the quasi-steady changes of the fluid flow of a partially saturated, fractured, porous system could be approximately simulated without taking the fractures into account. 22 references, 16 figures

Efficient kinetic method for fluid simulation beyond the Navier-Stokes equation.

Science.gov (United States)

Zhang, Raoyang; Shan, Xiaowen; Chen, Hudong

2006-10-01

We present a further theoretical extension to the kinetic-theory-based formulation of the lattice Boltzmann method of Shan [J. Fluid Mech. 550, 413 (2006)]. In addition to the higher-order projection of the equilibrium distribution function and a sufficiently accurate Gauss-Hermite quadrature in the original formulation, a regularization procedure is introduced in this paper. This procedure ensures a consistent order of accuracy control over the nonequilibrium contributions in the Galerkin sense. Using this formulation, we construct a specific lattice Boltzmann model that accurately incorporates up to third-order hydrodynamic moments. Numerical evidence demonstrates that the extended model overcomes some major defects existing in conventionally known lattice Boltzmann models, so that fluid flows at finite Knudsen number Kn can be more quantitatively simulated. Results from force-driven Poiseuille flow simulations predict the Knudsen's minimum and the asymptotic behavior of flow flux at large Kn.

Detached-eddy simulation of flow around the NREL phase VI blade

DEFF Research Database (Denmark)

Johansen, Jeppe; Sørensen, Niels N.; Michelsen, J.A.

2002-01-01

the blade axis. Computed blade characteristics are compared with experimental data from the NREL/NASA Ames Phase VI unsteady experiment. The detached-eddy simulation model is a method for predicting turbulence in computational fluid dynamics computations, which combines a Reynolds-averaged Navier......-eddy simulation show considerably more three-dimensional flow structures compared to conventional two-equation Reynolds-averaged Navier–Stokes turbulence models, but no particular improvements are seen in the global blade characteristics. Copyright © 2002 John Wiley & Sons, Ltd.......The detached-eddy simulation model implemented in the computational fluid dynamics code EllipSys3D is used to calculate the flow around the non-rotating NREL Phase VI wind turbine blade. Results are presented for flow around a parked blade at fixed angle of attack and a blade pitching along...

Numerical Modeling of Interstitial Fluid Flow Coupled with Blood Flow through a Remodeled Solid Tumor Microvascular Network.

Science.gov (United States)

Soltani, M; Chen, P

2013-01-01

Modeling of interstitial fluid flow involves processes such as fluid diffusion, convective transport in extracellular matrix, and extravasation from blood vessels. To date, majority of microvascular flow modeling has been done at different levels and scales mostly on simple tumor shapes with their capillaries. However, with our proposed numerical model, more complex and realistic tumor shapes and capillary networks can be studied. Both blood flow through a capillary network, which is induced by a solid tumor, and fluid flow in tumor's surrounding tissue are formulated. First, governing equations of angiogenesis are implemented to specify the different domains for the network and interstitium. Then, governing equations for flow modeling are introduced for different domains. The conservation laws for mass and momentum (including continuity equation, Darcy's law for tissue, and simplified Navier-Stokes equation for blood flow through capillaries) are used for simulating interstitial and intravascular flows and Starling's law is used for closing this system of equations and coupling the intravascular and extravascular flows. This is the first study of flow modeling in solid tumors to naturalistically couple intravascular and extravascular flow through a network. This network is generated by sprouting angiogenesis and consisting of one parent vessel connected to the network while taking into account the non-continuous behavior of blood, adaptability of capillary diameter to hemodynamics and metabolic stimuli, non-Newtonian blood flow, and phase separation of blood flow in capillary bifurcation. The incorporation of the outlined components beyond the previous models provides a more realistic prediction of interstitial fluid flow pattern in solid tumors and surrounding tissues. Results predict higher interstitial pressure, almost two times, for realistic model compared to the simplified model.

Numerical Modeling of Interstitial Fluid Flow Coupled with Blood Flow through a Remodeled Solid Tumor Microvascular Network.

Directory of Open Access Journals (Sweden)

M Soltani

Full Text Available Modeling of interstitial fluid flow involves processes such as fluid diffusion, convective transport in extracellular matrix, and extravasation from blood vessels. To date, majority of microvascular flow modeling has been done at different levels and scales mostly on simple tumor shapes with their capillaries. However, with our proposed numerical model, more complex and realistic tumor shapes and capillary networks can be studied. Both blood flow through a capillary network, which is induced by a solid tumor, and fluid flow in tumor's surrounding tissue are formulated. First, governing equations of angiogenesis are implemented to specify the different domains for the network and interstitium. Then, governing equations for flow modeling are introduced for different domains. The conservation laws for mass and momentum (including continuity equation, Darcy's law for tissue, and simplified Navier-Stokes equation for blood flow through capillaries are used for simulating interstitial and intravascular flows and Starling's law is used for closing this system of equations and coupling the intravascular and extravascular flows. This is the first study of flow modeling in solid tumors to naturalistically couple intravascular and extravascular flow through a network. This network is generated by sprouting angiogenesis and consisting of one parent vessel connected to the network while taking into account the non-continuous behavior of blood, adaptability of capillary diameter to hemodynamics and metabolic stimuli, non-Newtonian blood flow, and phase separation of blood flow in capillary bifurcation. The incorporation of the outlined components beyond the previous models provides a more realistic prediction of interstitial fluid flow pattern in solid tumors and surrounding tissues. Results predict higher interstitial pressure, almost two times, for realistic model compared to the simplified model.

STUDY ON SIMULATION METHOD OF AVALANCHE : FLOW ANALYSIS OF AVALANCHE USING PARTICLE METHOD

OpenAIRE

塩澤, 孝哉

2015-01-01

In this paper, modeling for the simulation of the avalanche by a particle method is discussed. There are two kinds of the snow avalanches, one is the surface avalanche which shows a smoke-like flow, and another is the total-layer avalanche which shows a flow like Bingham fluid. In the simulation of the surface avalanche, the particle method in consideration of a rotation resistance model is used. The particle method by Bingham fluid is used in the simulation of the total-layer avalanche. At t...

Computational fluid dynamics simulations of blood flow regularized by 3D phase contrast MRI

DEFF Research Database (Denmark)

Rispoli, Vinicius C; Nielsen, Jon; Nayak, Krishna S

2015-01-01

BACKGROUND: Phase contrast magnetic resonance imaging (PC-MRI) is used clinically for quantitative assessment of cardiovascular flow and function, as it is capable of providing directly-measured 3D velocity maps. Alternatively, vascular flow can be estimated from model-based computation fluid dyn...

Numerical investigation on lateral migration and lift force of single bubble in simple shear flow in low viscosity fluid using volume of fluid method

International Nuclear Information System (INIS)

Zhongchun, Li; Xiaoming, Song; Shengyao, Jiang; Jiyang, Yu

2014-01-01

Highlights: • A VOF simulation of bubble in low viscosity fluid was conducted. • Lift force in different viscosity fluid had different lateral migration characteristics. • Bubble with different size migrated to different direction. • Shear stress triggered the bubble deformation process and the bubble deformation came along with the oscillation behaviors. - Abstract: Two phase flow systems have been widely used in industrial engineering. Phase distribution characteristics are vital to the safety operation and optimization design of two phase flow systems. Lift force has been known as perpendicular to the bubbles’ moving direction, which is one of the mechanisms of interfacial momentum transfer. While most widely used lift force correlations, such as the correlation of Tomiyama et al. (2002), were obtained by experimentally tracking single bubble trajectories in high viscosity glycerol–water mixture, the applicability of these models into low viscosity fluid, such as water in nuclear engineering system, needs to be further evaluated. In the present paper, bubble in low viscosity fluid in shear flow was investigated in a full 3D numerical simulation and the volume of fluid (VOF) method was applied to capture the interface. The fluid parameter: fluid viscosity, bubble parameter: diameter and external flow parameters: shear stress magnitude and liquid velocity were examined. Comparing with bubble in high viscosity shear flow and bubble in low viscosity still flow, relative large bubble in low viscosity shear flow keep an oscillation way towards the moving wall and experienced a shape deformation process. The oscillation amplitude increased as the viscosity of fluid decreased. Small bubble migrated to the static wall in a line with larger migration velocity than that in high viscosity fluid and no deformation occurred. The shear stress triggered the oscillation behaviors while it had no direct influence with the behavior. The liquid velocity had no effect on

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)

Thermal effects on fluid flow and hydraulic fracturing from wellbores and cavities in low-permeability formations

Energy Technology Data Exchange (ETDEWEB)

Yarlong Wang [Petro-Geotech Inc., Calgary, AB (Canada); Papamichos, Euripides [IKU Petroleum Research, Trondheim (Norway)

1999-07-01

The coupled heat-fluid-stress problem of circular wellbore or spherical cavity subjected to a constant temperature change and a constant fluid flow rate is considered. Transient analytical solutions for temperature, pore pressure and stress are developed by coupling conductive heat transfer with Darcy fluid flow in a poroelastic medium. They are applicable to lower permeability porous media suitable for liquid-waste disposal and also simulating reservoir for enhanced oil recovery, where conduction dominates the heat transfer process. A full range of solutions is presented showing separately the effects of temperature and fluid flow on pore pressure and stress development. It is shown that injection of warm fluid can be used to restrict fracture development around wellbores and cavities and generally to optimise a fluid injection operation. Both the limitations of the solutions and the convective flow effect are addressed. (Author)

Validation of model predictions of pore-scale fluid distributions during two-phase flow

Science.gov (United States)

Bultreys, Tom; Lin, Qingyang; Gao, Ying; Raeini, Ali Q.; AlRatrout, Ahmed; Bijeljic, Branko; Blunt, Martin J.

2018-05-01

Pore-scale two-phase flow modeling is an important technology to study a rock's relative permeability behavior. To investigate if these models are predictive, the calculated pore-scale fluid distributions which determine the relative permeability need to be validated. In this work, we introduce a methodology to quantitatively compare models to experimental fluid distributions in flow experiments visualized with microcomputed tomography. First, we analyzed five repeated drainage-imbibition experiments on a single sample. In these experiments, the exact fluid distributions were not fully repeatable on a pore-by-pore basis, while the global properties of the fluid distribution were. Then two fractional flow experiments were used to validate a quasistatic pore network model. The model correctly predicted the fluid present in more than 75% of pores and throats in drainage and imbibition. To quantify what this means for the relevant global properties of the fluid distribution, we compare the main flow paths and the connectivity across the different pore sizes in the modeled and experimental fluid distributions. These essential topology characteristics matched well for drainage simulations, but not for imbibition. This suggests that the pore-filling rules in the network model we used need to be improved to make reliable predictions of imbibition. The presented analysis illustrates the potential of our methodology to systematically and robustly test two-phase flow models to aid in model development and calibration.

Modelling and numerical simulation of two-phase flows using the two-fluid two-pressure approach; Modelisation et simulation numerique des ecoulements diphasiques par une approche bifluide a deux pressions

Energy Technology Data Exchange (ETDEWEB)

Guillemaud, V

2007-03-15

This thesis is devoted to the modelling and numerical simulation of liquid-vapor flows. In order to describe these phase transition flows, a two-fluid two-pressure approach is considered. This description of the liquid-vapor mixing is associated to the seven-equation model introduced by Baer and Nunziato. This work investigates the properties of this model in order to simulate the phase transition flows occurring in nuclear engineering. First, a theoretical thermodynamic framework is constructed to describe the liquid-vapor mixing. Provided with this framework, various modelling choices are suggested for the interaction terms between the phases. These closure laws comply with an entropy inequality. The mathematical properties of this model are thereafter examined. The convective part is associated to a nonconservative hyperbolic system. First, we focus on the definition of its weak solutions. Several flow regimes for the two-phase mixing derive from this analysis. Such regimes for the two-phase flows are analogous to the torrential and fluvial regimes for the shallow-water equations. Furthermore, we establish the linear and nonlinear stabilities of the liquid-vapor equilibrium. Finally, the implementation of a turbulence model and the introduction of a reconstruction process for the interfacial area are investigated in order to refine the description of the interfacial transfers. Using a fractional step approach, a Finite Volume method is at last constructed to simulate this model. First, various nonconservative adaptations of standard Riemann solvers are developed to approach the convective part. Unlike the classic nonconservative framework, these schemes converge towards the same solution. Furthermore, a new relaxation scheme is proposed to approach the interfacial transfers. Provided with these schemes, the whole numerical method preserves the liquid-vapor equilibria. Using this numerical method, a careful comparison between the one- and two-pressure two-fluid

Numerical simulation of secondary flow in bubbly turbulent flow in sub-channel

International Nuclear Information System (INIS)

Ikeno, Tsutomu; Kataoka, Isao

2009-01-01

Secondary flow in bubbly turbulent flow in sub-channel was simulated by using an algebraic turbulence stress model. The mass, momentum, turbulence energy and bubble diffusion equations were used as fundamental equation. The basis for these equations was the two-fluid model: the equation of liquid phase was picked up from the equation system theoretically derived for the gas-liquid two-fluid turbulent flow. The fundamental equation was transformed onto a generalized coordinate system fitted to the computational domain in sub-channel. It was discretized for the SIMPLE algorism using the finite-volume method. The shape of sub-channel causes a distortion of the computational mesh, and orthogonal nature of the mesh is sometimes broken. An iterative method to satisfy a requirement for the contra-variant velocity was introduced to represent accurate symmetric boundary condition. Two-phase flow at a steady state was simulated for different magnitude of secondary flow and void fraction. The secondary flow enhanced the momentum transport in sub-channel and accelerated the liquid phase in the rod gap. This effect was slightly mitigated when the void fraction increased. The acceleration can contribute to effective cooling in the rod gap. The numerical result implied a phenomenon of industrial interest. This suggested that experimental approach is necessary to validate the numerical model and to identify the phenomenon. (author)

Lattice Boltzmann Simulation of Permeability and Tortuosity for Flow through Dense Porous Media

Directory of Open Access Journals (Sweden)

Ping Wang

2014-01-01

Full Text Available Discrete element method (DEM is used to produce dense and fixed porous media with rigid mono spheres. Lattice Boltzmann method (LBM is adopted to simulate the fluid flow in interval of dense spheres. To simulating the same physical problem, the permeability is obtained with different lattice number. We verify that the permeability is irrelevant to the body force and the media length along flow direction. The relationships between permeability, tortuosity and porosity, and sphere radius are researched, and the results are compared with those reported by other authors. The obtained results indicate that LBM is suited to fluid flow simulation of porous media due to its inherent theoretical advantages. The radius of sphere should have ten lattices at least and the media length along flow direction should be more than twenty radii. The force has no effect on the coefficient of permeability with the limitation of slow fluid flow. For mono spheres porous media sample, the relationship of permeability and porosity agrees well with the K-C equation, and the tortuosity decreases linearly with increasing porosity.

Large interface simulation in an averaged two-fluid code

International Nuclear Information System (INIS)

Henriques, A.

2006-01-01

Different ranges of size of interfaces and eddies are involved in multiphase flow phenomena. Classical formalisms focus on a specific range of size. This study presents a Large Interface Simulation (LIS) two-fluid compressible formalism taking into account different sizes of interfaces. As in the single-phase Large Eddy Simulation, a filtering process is used to point out Large Interface (LI) simulation and Small interface (SI) modelization. The LI surface tension force is modelled adapting the well-known CSF method. The modelling of SI transfer terms is done calling for classical closure laws of the averaged approach. To simulate accurately LI transfer terms, we develop a LI recognition algorithm based on a dimensionless criterion. The LIS model is applied in a classical averaged two-fluid code. The LI transfer terms modelling and the LI recognition are validated on analytical and experimental tests. A square base basin excited by a horizontal periodic movement is studied with the LIS model. The capability of the model is also shown on the case of the break-up of a bubble in a turbulent liquid flow. The break-up of a large bubble at a grid impact performed regime transition between two different scales of interface from LI to SI and from PI to LI. (author) [fr

Visualization and simulation of complex flows in biomedical engineering

CERN Document Server

Imai, Yohsuke; Ishikawa, Takuji; Oliveira, Mónica

2014-01-01

This book focuses on the most recent advances in the application of visualization and simulation methods to understand the flow behavior of complex fluids used in biomedical engineering and other related fields. It shows the physiological flow behavior in large arteries, microcirculation, respiratory systems and in biomedical microdevices.

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.

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

Development of simulation programs for three dimensional movement of fluid

Energy Technology Data Exchange (ETDEWEB)

Kim, Bok Youn; Kang, Chang Hee; Jo, Young Do; Lim, Sang Taek [Korea Institute of Geology Mining and Materials, Taejon (Korea, Republic of)

1996-12-01

It is third project year on `Application of mobile diesel equipment in underground mines` for providing appropriate measures to improve underground working environment contaminated by the diesel exhaust pollutants. For simulating three dimensional movement of fluid in closed spaces like dead-end workings ; a software (3D-Flow) based on computerized fluid dynamics (CFD) has been developed. According to simulations using this software, following results are derived. 1) The heading faces where a diesel equipment is employed should be ventilated even though the extension is as short as 20 meters long. 2) Even there is no ventilation, the gas concentration of the heading face reaches no more than 15 % in ceiling and 5 % in the region where workers are in action approximately. 3) Gases are formed regular horizontal laminar flow when there is no ventilation. 4) In ventilated headings, it is same that the gas concentration of working region is about 30 % of the concentration of ceiling and gases are formed irregular turbulent flow. (author). 2 tabs., 15 figs.

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.

Generalized Fluid System Simulation Program, Version 6.0

Science.gov (United States)

Majumdar, A. K.; LeClair, A. C.; Moore, R.; Schallhorn, P. A.

2016-01-01

The Generalized Fluid System Simulation Program (GFSSP) is a general purpose computer program for analyzing steady state and time-dependent flow rates, pressures, temperatures, and concentrations in a complex flow network. The program is capable of modeling real fluids with phase changes, compressibility, mixture thermodynamics, conjugate heat transfer between solid and fluid, fluid transients, pumps, compressors, and external body forces such as gravity and centrifugal. The thermofluid system to be analyzed is discretized into nodes, branches, and conductors. The scalar properties such as pressure, temperature, and concentrations are calculated at nodes. Mass flow rates and heat transfer rates are computed in branches and conductors. The graphical user interface allows users to build their models using the 'point, drag, and click' method; the users can also run their models and post-process the results in the same environment. Two thermodynamic property programs (GASP/WASP and GASPAK) provide required thermodynamic and thermophysical properties for 36 fluids: helium, methane, neon, nitrogen, carbon monoxide, oxygen, argon, carbon dioxide, fluorine, hydrogen, parahydrogen, water, kerosene (RP-1), isobutene, butane, deuterium, ethane, ethylene, hydrogen sulfide, krypton, propane, xenon, R-11, R-12, R-22, R-32, R-123, R-124, R-125, R-134A, R-152A, nitrogen trifluoride, ammonia, hydrogen peroxide, and air. The program also provides the options of using any incompressible fluid with constant density and viscosity or ideal gas. The users can also supply property tables for fluids that are not in the library. Twenty-four different resistance/source options are provided for modeling momentum sources or sinks in the branches. These options include pipe flow, flow through a restriction, noncircular duct, pipe flow with entrance and/or exit losses, thin sharp orifice, thick orifice, square edge reduction, square edge expansion, rotating annular duct, rotating radial duct

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

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)

An Embedded 3D Fracture Modeling Approach for Simulating Fracture-Dominated Fluid Flow and Heat Transfer in Geothermal Reservoirs

Energy Technology Data Exchange (ETDEWEB)

Johnston, Henry [National Renewable Energy Laboratory (NREL), Golden, CO (United States); Wang, Cong [Colorado School of Mines; Winterfeld, Philip [Colorado School of Mines; Wu, Yu-Shu [Colorado School of Mines

2018-02-14

An efficient modeling approach is described for incorporating arbitrary 3D, discrete fractures, such as hydraulic fractures or faults, into modeling fracture-dominated fluid flow and heat transfer in fractured geothermal reservoirs. This technique allows 3D discrete fractures to be discretized independently from surrounding rock volume and inserted explicitly into a primary fracture/matrix grid, generated without including 3D discrete fractures in prior. An effective computational algorithm is developed to discretize these 3D discrete fractures and construct local connections between 3D fractures and fracture/matrix grid blocks of representing the surrounding rock volume. The constructed gridding information on 3D fractures is then added to the primary grid. This embedded fracture modeling approach can be directly implemented into a developed geothermal reservoir simulator via the integral finite difference (IFD) method or with TOUGH2 technology This embedded fracture modeling approach is very promising and computationally efficient to handle realistic 3D discrete fractures with complicated geometries, connections, and spatial distributions. Compared with other fracture modeling approaches, it avoids cumbersome 3D unstructured, local refining procedures, and increases computational efficiency by simplifying Jacobian matrix size and sparsity, while keeps sufficient accuracy. Several numeral simulations are present to demonstrate the utility and robustness of the proposed technique. Our numerical experiments show that this approach captures all the key patterns about fluid flow and heat transfer dominated by fractures in these cases. Thus, this approach is readily available to simulation of fractured geothermal reservoirs with both artificial and natural fractures.

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.

Numerical simulations of air–water cap-bubbly flows using two-group interfacial area transport equation

International Nuclear Information System (INIS)

Wang, Xia; Sun, Xiaodong

2014-01-01

Highlights: • Two-group interfacial area transport equation was implemented into a three-field two-fluid model in Fluent. • Numerical model was developed for cap-bubbly flows in a narrow rectangular flow channel. • Numerical simulations were performed for cap-bubbly flows with uniform void inlets and with central peaked void inlets. • Code simulations showed a significant improve over the conventional two-fluid model. - Abstract: Knowledge of cap-bubbly flows is of great interest due to its role in understanding of the flow regime transition from bubbly to slug or churn-turbulent flows. One of the key characteristics of such flows is the existence of bubbles in different sizes and shapes associated with their distinctive dynamic natures. This important feature is, however, generally not well captured by many available two-phase flow modeling approaches. In this study, a modified two-fluid model, namely a three-field, two-fluid model, is proposed. In this model, bubbles are categorized into two groups, i.e., spherical/distorted bubbles as Group-1 while cap/churn-turbulent bubbles as Group-2. A two-group interfacial area transport equation (IATE) is implemented to describe dynamic changes of interfacial structure in each bubble group, resulting from intra- and inter-group interactions and phase changes due to evaporation and condensation. Attention is also paid to appropriate constitutive relations of the interfacial transfers due to mechanical and thermal non-equilibrium between the different fields. The proposed three-field, two-fluid model is used to predict the phase distributions of adiabatic air–water flows in a confined rectangular duct. Good agreement between the simulation results from the proposed model and relevant experimental data indicates that the proposed model is promising as an improved computational tool for two-phase cap-bubbly flow simulations in rectangular flow ducts

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.

Cilia induced cerebrospinal fluid flow in the third ventricle of brain

Science.gov (United States)

Wang, Yong; Westendorf, Christian; Faubel, Regina; Eichele, Gregor; Bodenschatz, Eberhard

2016-11-01

Cerebrospinal fluid (CSF) conveys many physiologically important signaling factors through the ventricles of the mammalian brain. The walls of the ventricles are covered with motile cilia that were thought to generate a laminar flow purely following the curvature of walls. However, we recently discovered that cilia of the ventral third ventricle (v3V) generate a complex flow network along the wall, leading to subdivision of the v3V. The contribution of such cilia induced flow to the overall three dimensional volume flow remains to be investigated by using numerical simulation, arguably the best approach for such investigations. The lattice Boltzmann method is used to study the CFS flow in a reconstructed geometry of the v3V. Simulation of CSF flow neglecting cilia in this geometry confirmed that the previous idea about pure confined flow does not reflect the reality observed in experiment. The experimentally recorded ciliary flow network along the wall was refined with the smoothed particle hydrodynamics and then adapted as boundary condition in simulation. We study the contribution of the ciliary network to overall CSF flow and identify site-specific delivery of CSF constituents with respect to the temporal changes.

Numerical simulation of heat transfer in power law fluid flow through a stenosed artery

Science.gov (United States)

Talib, Amira Husni; Abdullah, Ilyani

2017-11-01

A numerical study of heat transfer in a power law fluid is investigated in this paper. The blood flow is treated as power law fluid with a presence of cosine shaped stenosis. This study reveals the effect of stenosis on the heat transfer and velocity of blood flowing in the constricted artery. The governing and energy equations are formulated in a cylindrical coordinate system. Hence, the set of equations and boundary conditions are solved numerically by Marker and Cell (MAC) method. The graphical result shows the profile of blood temperature is increased while the blood velocity is decreased at the critical height of stenosis.

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.

Tracer technology modeling the flow of fluids

CERN Document Server

Levenspiel, Octave

2012-01-01

A vessel’s behavior as a heat exchanger, absorber, reactor, or other process unit is dependent upon how fluid flows through the vessel.  In early engineering, the designer would assume either plug flow or mixed flow of the fluid through the vessel.  However, these assumptions were oftentimes inaccurate, sometimes being off by a volume factor of 100 or more.  The result of this unreliable figure produced ineffective products in multiple reaction systems.   Written by a pioneering researcher in the field of chemical engineering, the tracer method was introduced to provide more accurate flow data.  First, the tracer method measured the actual flow of fluid through a vessel.  Second, it developed a suitable model to represent the flow in question.  Such models are used to follow the flow of fluid in chemical reactors and other process units, like in rivers and streams, or solid and porous structures.  In medicine, the tracer method is used to study the flow of chemicals—harmful  and harmless—in the...

Numerical simulation of the effect of upstream swirling flow on swirl meter performance

Science.gov (United States)

Chen, Desheng; Cui, Baoling; Zhu, Zuchao

2018-04-01

Flow measurement is important in the fluid process and transmission system. For the need of accuracy measurement of fluid, stable flow is acquired. However, the elbows and devices as valves and rotary machines may produce swirling flow in the natural gas pipeline networks system and many other industry fields. In order to reveal the influence of upstream swirling flow on internal flow fields and the metrological characteristics, numerical simulations are carried out on the swirl meter. Using RNG k-ɛ turbulent model and SIMPLE algorithm, the flow field is numerically simulated under swirling flows generated from co-swirl and counter-swirl flow. Simulation results show fluctuation is enhanced or weakened depending on the rotating direction of swirling flow. A counter- swirl flow increases the entropy production rate at the inlet and outlet of the swirler, the junction region between throat and divergent section, and then the pressure loss is increased. The vortex precession dominates the static pressure distributions on the solid walls and in the channel, especially at the end region of the throat.

Turbulence prediction in two-dimensional bundle flows using large eddy simulation

Energy Technology Data Exchange (ETDEWEB)

Ibrahim, W.A.; Hassan, Y.A. [Texas A& M Univ., College Station, TX (United States)

1995-09-01

Turbulent flow is characterized by random fluctuations in the fluid velocity and by intense mixing of the fluid. Due to velocity fluctuations, a wide range of eddies exists in the flow field. Because these eddies carry mass, momentum, and energy, this enhanced mixing can sometimes lead to serious problems, such as tube vibrations in many engineering systems that include fluid-tube bundle combinations. Nuclear fuel bundles and PWR steam generators are existing examples in nuclear power plants. Fluid-induced vibration problems are often discovered during the operation of such systems because some of the fluid-tube interaction characteristics are not fully understood. Large Eddy Simulation, incorporated in a three dimensional computer code, became one of the promising techniques to estimate flow turbulence, predict and prevent of long-term tube fretting affecting PWR steam generators. the present turbulence investigations is a step towards more understanding of fluid-tube interaction characteristics by comparing the tube bundles with various pitch-to-diameter ratios were performed. Power spectral densities were used for comparison with experimental data. Correlations, calculations of different length scales in the flow domain and other important turbulent-related parameters were calculated. Finally, important characteristics of turbulent flow field were presented with the aid of flow visualization with tracers impeded in the flow field.

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.

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

Moving on to the modeling and simulation using computational fluid dynamics

International Nuclear Information System (INIS)

Norasalwa Zakaria; Rohyiza Baan; Muhd Noor Muhd Yunus

2006-01-01

The heat is on but not at the co-combustor plant. Using the Computational Fluid Dynamics (CFD), modeling and simulation of an incinerator has been made easy and possible from the comfort of cozy room. CFD has become an important design tool in nearly every industrial field because it provides understanding of flow patterns. CFD provide values for fluid velocity, fluid temperature, pressure and species concentrations throughout a flow domain. MINT has acquired a complete CFD software recently, consisting of GAMBIT, which is use to build geometry and meshing, and FLUENT as the processor or solver. This paper discusses on several trial runs that was carried out on several parts of the co-combustor plant namely the under fire section and the mixing chamber section

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.

Numerical simulation of an elementary Vortex-Induced-Vibration problem by using fully-coupled fluid solid system computation

Directory of Open Access Journals (Sweden)

M Pomarède

2016-09-01

Full Text Available Numerical simulation of Vortex-Induced-Vibrations (VIV of a rigid circular elastically-mounted cylinder submitted to a fluid cross-flow has been extensively studied over the past decades, both experimentally and numerically, because of its theoretical and practical interest for understanding Flow-Induced-Vibrations (FIV problems. In this context, the present article aims to expose a numerical study based on fully-coupled fluid-solid computations compared to previously published work [34], [36]. The computational procedure relies on a partitioned method ensuring the coupling between fluid and structure solvers. The fluid solver involves a moving mesh formulation for simulation of the fluid structure interface motion. Energy exchanges between fluid and solid models are ensured through convenient numerical schemes. The present study is devoted to a low Reynolds number configuration. Cylinder motion magnitude, hydrodynamic forces, oscillation frequency and fluid vortex shedding modes are investigated and the “lock-in” phenomenon is reproduced numerically. These numerical results are proposed for code validation purposes before investigating larger industrial applications such as configurations involving tube arrays under cross-flows [4].

Simulation of granular and gas-solid flows using discrete element method

Science.gov (United States)

Boyalakuntla, Dhanunjay S.

2003-10-01

In recent years there has been increased research activity in the experimental and numerical study of gas-solid flows. Flows of this type have numerous applications in the energy, pharmaceuticals, and chemicals process industries. Typical applications include pulverized coal combustion, flow and heat transfer in bubbling and circulating fluidized beds, hopper and chute flows, pneumatic transport of pharmaceutical powders and pellets, and many more. The present work addresses the study of gas-solid flows using computational fluid dynamics (CFD) techniques and discrete element simulation methods (DES) combined. Many previous studies of coupled gas-solid flows have been performed assuming the solid phase as a continuum with averaged properties and treating the gas-solid flow as constituting of interpenetrating continua. Instead, in the present work, the gas phase flow is simulated using continuum theory and the solid phase flow is simulated using DES. DES treats each solid particle individually, thus accounting for its dynamics due to particle-particle interactions, particle-wall interactions as well as fluid drag and buoyancy. The present work involves developing efficient DES methods for dense granular flow and coupling this simulation to continuum simulations of the gas phase flow. Simulations have been performed to observe pure granular behavior in vibrating beds. Benchmark cases have been simulated and the results obtained match the published literature. The dimensionless acceleration amplitude and the bed height are the parameters governing bed behavior. Various interesting behaviors such as heaping, round and cusp surface standing waves, as well as kinks, have been observed for different values of the acceleration amplitude for a given bed height. Furthermore, binary granular mixtures (granular mixtures with two particle sizes) in a vibrated bed have also been studied. Gas-solid flow simulations have been performed to study fluidized beds. Benchmark 2D

Simultaneous measurement of particle and fluid velocities in particle-laden flows

International Nuclear Information System (INIS)

Jin, D. X.; Lee, D. Y.

2009-01-01

For the velocity measurement in a particle-laden fluid flow, the fluid velocity and the inherently dispersed particle velocity can be analyzed by using PIV and PTV, respectively. Since the PIV result statistically represents the average displacement of all the particles in a PIV image, it is inevitable that the PIV result includes the influence of the dispersed particles' displacement if a single CCD camera is used to simultaneously measure the fluid velocity and the dispersed particle velocity. The influence of dispersed particles should be excluded before the PIV analysis in order to evaluate the fluid velocity accurately. In this study, the optimum replacement brightness of dispersed particles to minimize the false influence of dispersed particles on the PIV analysis was theoretically derived. Simulation results show that the modification of dispersed particle brightness can significantly reduce the PIV error caused by the dispersed particles. This modification method was also verified in the analysis of an actual experimental case of the particle-laden fluid flow in a triangular grooved channel

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...

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.

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.

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.

A Numerical Study of Mesh Adaptivity in Multiphase Flows with Non-Newtonian Fluids

Science.gov (United States)

Percival, James; Pavlidis, Dimitrios; Xie, Zhihua; Alberini, Federico; Simmons, Mark; Pain, Christopher; Matar, Omar

2014-11-01

We present an investigation into the computational efficiency benefits of dynamic mesh adaptivity in the numerical simulation of transient multiphase fluid flow problems involving Non-Newtonian fluids. Such fluids appear in a range of industrial applications, from printing inks to toothpastes and introduce new challenges for mesh adaptivity due to the additional ``memory'' of viscoelastic fluids. Nevertheless, the multiscale nature of these flows implies huge potential benefits for a successful implementation. The study is performed using the open source package Fluidity, which couples an unstructured mesh control volume finite element solver for the multiphase Navier-Stokes equations to a dynamic anisotropic mesh adaptivity algorithm, based on estimated solution interpolation error criteria, and conservative mesh-to-mesh interpolation routine. The code is applied to problems involving rheologies ranging from simple Newtonian to shear-thinning to viscoelastic materials and verified against experimental data for various industrial and microfluidic flows. This work was undertaken as part of the EPSRC MEMPHIS programme grant EP/K003976/1.

A new algorithm for extended nonequilibrium molecular dynamics simulations of mixed flow

NARCIS (Netherlands)

Hunt, T.A.; Hunt, Thomas A.; Bernardi, Stefano; Todd, B.D.

2010-01-01

In this work, we develop a new algorithm for nonequilibrium molecular dynamics of fluids under planar mixed flow, a linear combination of planar elongational flow and planar Couette flow. To date, the only way of simulating mixed flow using nonequilibrium molecular dynamics techniques was to impose

Numerical Simulation of Flow Behavior within a Venturi Scrubber

Directory of Open Access Journals (Sweden)

M. M. Toledo-Melchor

2014-01-01

Full Text Available The present work details the three-dimensional numerical simulation of single-phase and two-phase flow (air-water in a venturi scrubber with an inlet and throat diameters of 250 and 122.5 mm, respectively. The dimensions and operating parameters correspond to industrial applications. The mass flow rate conditions were 0.483 kg/s, 0.736 kg/s, 0.861 kg/s, and 0.987 kg/s for the gas only simulation; the mass flow rate for the liquid was 0.013 kg/s and 0.038 kg/s. The gas flow was simulated in five geometries with different converging and diverging angles while the two-phase flow was only simulated for one geometry. The results obtained were validated with experimental data obtained by other researchers. The results show that the pressure drop depends significantly on the gas flow rate and that water flow rate does not have significant effects neither on the pressure drop nor on the fluid maximum velocity within the scrubber.

The PDF method for Lagrangian two-phase flow simulations

International Nuclear Information System (INIS)

Minier, J.P.; Pozorski, J.

1996-04-01

A recent turbulence model put forward by Pope (1991) in the context of PDF modelling has been used. In this approach, the one-point joint velocity-dissipation pdf equation is solved by simulating the instantaneous behaviour of a large number of Lagrangian fluid particles. Closure of the evolution equations of these Lagrangian particles is based on stochastic models and more specifically on diffusion processes. Such models are of direct use for two-phase flow modelling where the so-called fluid seen by discrete inclusions has to be modelled. Full Lagrangian simulations have been performed for shear-flows. It is emphasized that this approach gives far more information than traditional turbulence closures (such as the K-ε model) and therefore can be very useful for situations involving complex physics. It is also believed that the present model represents the first step towards a complete Lagrangian-Lagrangian model for dispersed two-phase flow problems. (authors). 21 refs., 6 figs

2-D CFD time-dependent thermal-hydraulic simulations of CANDU-6 moderator flows

Energy Technology Data Exchange (ETDEWEB)

Mehdi Zadeh, Foad [Department of Engineering Physics/Polytechnique Montréal, Montréal, QC (Canada); Étienne, Stéphane [Department of Mechanical Engineering/Polytechnique Montréal, Montréal, QC (Canada); Teyssedou, Alberto, E-mail: alberto.teyssedou@polymtl.ca [Department of Engineering Physics/Polytechnique Montréal, Montréal, QC (Canada)

2016-12-01

Highlights: • 2-D time-dependent CFD simulations of CANDU-6 moderator flows are presented. • A thermal-hydraulic code using thermal physical fluid properties is used. • The numerical approach and convergence is validated against available data. • Flow configurations are correlated using Richardson’s number. • Frequency components indicate moderator flow oscillations vs. Richardson numbers. - Abstract: The distribution of the fluid temperature and mass density of the moderator flow in CANDU-6 nuclear power reactors may affect the reactivity coefficient. For this reason, any possible moderator flow configuration and consequently the corresponding temperature distributions must be studied. In particular, the variations of the reactivity may result in major safety issues. For instance, excessive temperature excursions in the vicinity of the calandria tubes nearby local flow stagnation zones, may bring about partial boiling. Moreover, steady-state simulations have shown that for operating condition, intense buoyancy forces may be dominant, which can trigger a thermal stratification. Therefore, the numerical study of the time-dependent flow transition to such a condition, is of fundamental safety concern. Within this framework, this paper presents detailed time-dependent numerical simulations of CANDU-6 moderator flow for a wide range of flow conditions. To get a better insight of the thermal-hydraulic phenomena, the simulations were performed by covering long physical-time periods using an open-source code (Code-Saturne V3) developed by Électricité de France. The results show not only a region where the flow is characterized by coherent structures of flow fluctuations but also the existence of two limit cases where fluid oscillations disappear almost completely.

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 and simulation of flow field in giant magnetostrictive pump

Science.gov (United States)

Zhao, Yapeng; Ren, Shiyong; Lu, Quanguo

2017-09-01

Recent years, there has been significant research in the design and analysis of giant magnetostrictive pump. In this paper, the flow field model of giant magnetostrictive pump was established and the relationship between pressure loss and working frequency of piston was studied by numerical simulation method. Then, the influence of different pump chamber height on pressure loss in giant magnetostrictive pump was studied by means of flow field simulation. Finally, the fluid pressure and velocity vector distribution in giant magnetostrictive pump chamber were simulated.

Modeling fluid flow and heat transfer at Basin and Range faults: preliminary results for Leach hot springs, Nevada

Science.gov (United States)

López, Dina L.; Smith, Leslie; Storey, Michael L.; Nielson, Dennis L.

1994-01-01

The hydrothermal systems of the Basin and Range Province are often located at or near major range bounding normal faults. The flow of fluid and energy at these faults is affected by the advective transfer of heat and fluid from an to the adjacent mountain ranges and valleys, This paper addresses the effect of the exchange of fluid and energy between the country rock, the valley fill sediments, and the fault zone, on the fluid and heat flow regimes at the fault plane. For comparative purposes, the conditions simulated are patterned on Leach Hot Springs in southern Grass Valley, Nevada. Our simulations indicated that convection can exist at the fault plane even when the fault is exchanging significant heat and fluid with the surrounding country rock and valley fill sediments. The temperature at the base of the fault decreased with increasing permeability of the country rock. Higher groundwater discharge from the fault and lower temperatures at the base of the fault are favored by high country rock permabilities and fault transmissivities. Preliminary results suggest that basal temperatures and flow rates for Leach Hot Springs can not be simulated with a fault 3 km deep and an average regional heat flow of 150 mW/m2 because the basal temperature and mass discharge rates are too low. A fault permeable to greater depths or a higher regional heat flow may be indicated for these springs.

Vibration energy harvesting in a small channel fluid flow using piezoelectric transducer

Energy Technology Data Exchange (ETDEWEB)

Hassan, Md. Mehedi, E-mail: buetmehedi10@gmail.com; Hossain, Md. Yeam, E-mail: yeamhossain@gmail.com; Mazumder, Rakib, E-mail: rakibmazumder46075@gmail.com; Rahman, Roussel, E-mail: roussel.rahman@gmail.com; Rahman, Md. Ashiqur, E-mail: ashiqurrahman@me.buet.ac.bd [Department of Mechanical Engineering, Bangladesh University of Engineering and Technology (BUET), Dhaka-1000 (Bangladesh)

2016-07-12

This work is aimed at developing a way to harvest energy from a fluid stream with the application of piezoelectric transducers in a small channel. In this COMSOL Multiphysics based simulation study, it is attempted to harvest energy from the abundant renewable source of energy available in the form of kinetic energy of naturally occurring flow of fluids. The strategy involves harnessing energy from a fluid-actuator through generation of couples, eddies and vortices, resulting from the stagnation and separation of flow around a semi-circular bluff-body attached to a cantilever beam containing a piezoceramic layer. Fluctuation of fluidic pressure impulse on the beam due to vortex shedding and varying lift forces causes the flexible cantilever beam to oscillate in the direction normal to the fluid flow in a periodic manner. The periodic application and release of a mechanical strain upon the beam effected a generation of electric potential within the piezoelectric layer, thus enabling extraction of electrical energy from the kinetic energy of the fluid. The piezoelectric material properties and transducer design are kept unchanged throughout the study, whereas the configuration is tested with different fluids and varying flow characteristics. The size and geometry of the obstructing entity are systematically varied to closely inspect the output from different iterations and for finding the optimum design parameters. The intermittent changes in the generated forces and subsequent variation in the strain on the beam are also monitored to find definitive relationship with the electrical energy output.

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.

Computational Fluid Dynamics and Visualisation of Coastal Flows in Tidal Channels Supporting Ocean Energy Development

Directory of Open Access Journals (Sweden)

Enayatollah Zangiabadi

2015-06-01

Full Text Available Flow characteristics in coastal regions are strongly influenced by the topography of the seabed and understanding the fluid dynamics is necessary before installation of tidal stream turbines (TST. In this paper, the bathymetry of a potential TST deployment site is used in the development of the a CFD (Computational Fluid Dynamics model. The steady state k-ϵ and transient Large Eddy Simulation (LES turbulence methods are employed and compared. The simulations are conducted with a fixed representation of the ocean surface, i.e., a rigid lid representation. In the vicinity of Horse Rock a study of the pressure difference shows that the small change in height of the water column is negligible, providing confidence in the simulation results. The stream surface method employed to visualise the results has important inherent characteristics that can enhance the visual perception of complex flow structures. The results of all cases are compared with the flow data transect gathered by an Acoustic Doppler Current Profiler (ADCP. It has been understood that the k-ϵ method can predict the flow pattern relatively well near the main features of the domain and the LES model has the ability to simulate some important flow patterns caused by the bathymetry.

Multiphase reacting flows modelling and simulation

CERN Document Server

Marchisio, Daniele L

2007-01-01

The papers in this book describe the most widely applicable modeling approaches and are organized in six groups covering from fundamentals to relevant applications. In the first part, some fundamentals of multiphase turbulent reacting flows are covered. In particular the introduction focuses on basic notions of turbulence theory in single-phase and multi-phase systems as well as on the interaction between turbulence and chemistry. In the second part, models for the physical and chemical processes involved are discussed. Among other things, particular emphasis is given to turbulence modeling strategies for multiphase flows based on the kinetic theory for granular flows. Next, the different numerical methods based on Lagrangian and/or Eulerian schemes are presented. In particular the most popular numerical approaches of computational fluid dynamics codes are described (i.e., Direct Numerical Simulation, Large Eddy Simulation, and Reynolds-Averaged Navier-Stokes approach). The book will cover particle-based meth...

Rheology of granular flows immersed in a viscous fluid

International Nuclear Information System (INIS)

Amarsid, Lhassan

2015-01-01

We investigate the behavior of granular materials immersed in a viscous fluid by means of extensive simulations based on the Discrete Element Method for particle dynamics coupled with the Lattice Boltzmann method for the fluid. We show that, for a broad range of parameters such as shear rate, confining stress and viscosity, the internal friction coefficient and packing fraction are well described by a single 'visco-inertial' dimensionless parameter combining inertial and Stokes numbers. The frictional behavior under constant confining pressure is mapped into a viscous behavior under volume-controlled conditions, leading to the divergence of the effective normal and shear viscosities in inverse square of the distance to the critical packing fraction. The results are in excellent agreement with the experimental data of Boyer et al. (2011). The evolution of the force network in terms of connectivity and anisotropy as a function of the visco-inertial number, indicates that the increase of frictional strength is a direct consequence of structural anisotropy enhanced by both fluid viscosity and grain inertia. In view of application to a potential nuclear accident, we also study the fragmentation and flow of confined porous aggregates in a fluid under the action of local overpressures and pressure gradients as well as gravity-driven flow of immersed particles in an hourglass. (author)

Fluid flow simulation and permeability computation in deformed porous carbonate grainstones

Science.gov (United States)

Zambrano, Miller; Tondi, Emanuele; Mancini, Lucia; Lanzafame, Gabriele; Trias, F. Xavier; Arzilli, Fabio; Materazzi, Marco; Torrieri, Stefano

2018-05-01

In deformed porous carbonates, the architecture of the pore network may be modified by deformation or diagenetic processes altering the permeability with respect to the pristine rock. The effects of the pore texture and morphology on permeability in porous rocks have been widely investigated due to the importance during the evaluation of geofluid reservoirs. In this study, these effects are assessed by combining synchrotron X-ray computed microtomography (SR micro-CT) and computational fluid dynamics. The studied samples pertain to deformed porous carbonate grainstones highly affected by deformation bands (DBs) exposed in Northwestern Sicily and Abruzzo regions, Italy. The high-resolution SR micro-CT images of the samples, acquired at the SYRMEP beamline of the Elettra - Sincrotrone Trieste laboratory (Italy), were used for simulating a pressure-driven flow by using the lattice-Boltzmann method (LBM). For the experiments, a multiple relaxation time (MRT) model with the D3Q19 scheme was used to avoid viscosity-dependent results of permeability. The permeability was calculated using Darcy's law once steady conditions were reached. After the simulations, the pore-network properties (effective porosity, specific surface area, and geometrical tortuosity) were calculated using 3D images of the velocity fields. These images were segmented considering a velocity threshold value higher than zero. The study showed that DBs may generate significant heterogeneity and anisotropy of the permeability of the evaluated rock samples. Cataclasis and cementation process taking place within the DBs reduce the effective porosity and therefore the permeability. Contrary to this, pressure dissolution and faulting may generate connected channels which contribute to the permeability only parallel to the DB.

Numerical simulation of the planar extrudate swell of pseudoplastic and viscoelastic fluids with the streamfunction and the VOF methods

DEFF Research Database (Denmark)

Comminal, Raphaël; Pimenta, Francisco; Hattel, Jesper H.

2018-01-01

, as well as with numerical simulations performed with the open-source rheoTool toolbox in OpenFOAM®. While the simulations of the generalized Newtonian fluids achieved mesh independence for all the methods tested, the flow simulations of the viscoelastic fluids are more sensitive to mesh refinement...

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...

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.

Fluid-structure interaction analysis of the flow through a stenotic aortic valve

Science.gov (United States)

Maleki, Hoda; Labrosse, Michel R.; Durand, Louis-Gilles; Kadem, Lyes

2009-11-01

In Europe and North America, aortic stenosis (AS) is the most frequent valvular heart disease and cardiovascular disease after systemic hypertension and coronary artery disease. Understanding blood flow through an aortic stenosis and developing new accurate non-invasive diagnostic parameters is, therefore, of primarily importance. However, simulating such flows is highly challenging. In this study, we considered the interaction between blood flow and the valve leaflets and compared the results obtained in healthy valves with stenotic ones. One effective method to model the interaction between the fluid and the structure is to use Arbitrary Lagrangian-Eulerian (ALE) approach. Our two-dimensional model includes appropriate nonlinear and anisotropic materials. It is loaded during the systolic phase by applying pressure curves to the fluid domain at the inflow. For modeling the calcified stenotic valve, calcium will be added on the aortic side of valve leaflets. Such simulations allow us to determine the effective orifice area of the valve, one of the main parameters used clinically to evaluate the severity of an AS, and to correlate it with changes in the structure of the leaflets.

A parallel offline CFD and closed-form approximation strategy for computationally efficient analysis of complex fluid flows

Science.gov (United States)

Allphin, Devin

Computational fluid dynamics (CFD) solution approximations for complex fluid flow problems have become a common and powerful engineering analysis technique. These tools, though qualitatively useful, remain limited in practice by their underlying inverse relationship between simulation accuracy and overall computational expense. While a great volume of research has focused on remedying these issues inherent to CFD, one traditionally overlooked area of resource reduction for engineering analysis concerns the basic definition and determination of functional relationships for the studied fluid flow variables. This artificial relationship-building technique, called meta-modeling or surrogate/offline approximation, uses design of experiments (DOE) theory to efficiently approximate non-physical coupling between the variables of interest in a fluid flow analysis problem. By mathematically approximating these variables, DOE methods can effectively reduce the required quantity of CFD simulations, freeing computational resources for other analytical focuses. An idealized interpretation of a fluid flow problem can also be employed to create suitably accurate approximations of fluid flow variables for the purposes of engineering analysis. When used in parallel with a meta-modeling approximation, a closed-form approximation can provide useful feedback concerning proper construction, suitability, or even necessity of an offline approximation tool. It also provides a short-circuit pathway for further reducing the overall computational demands of a fluid flow analysis, again freeing resources for otherwise unsuitable resource expenditures. To validate these inferences, a design optimization problem was presented requiring the inexpensive estimation of aerodynamic forces applied to a valve operating on a simulated piston-cylinder heat engine. The determination of these forces was to be found using parallel surrogate and exact approximation methods, thus evidencing the comparative

Approaching multiphase flows from the perspective of computational fluid dynamics

International Nuclear Information System (INIS)

Banas, A.O.

1992-01-01

Thermalhydraulic simulation methodologies based on subchannel and porous-medium concepts are briefly reviewed and contrasted with the general approach of Computational Fluid Dynamics (CFD). An outline of the advanced CFD methods for single-phase turbulent flows is followed by a short discussion of the unified formulation of averaged equations for turbulent and multiphase flows. Some of the recent applications of CFD at Chalk River Laboratories are discussed, and the complementary role of CFD with regard to the established thermalhydraulic methods of analysis is indicated. (author). 8 refs

Direct Numerical Simulation of Low Capillary Number Pore Scale Flows

Science.gov (United States)

Esmaeilzadeh, S.; Soulaine, C.; Tchelepi, H.

2017-12-01

The arrangement of void spaces and the granular structure of a porous medium determines multiple macroscopic properties of the rock such as porosity, capillary pressure, and relative permeability. Therefore, it is important to study the microscopic structure of the reservoir pores and understand the dynamics of fluid displacements through them. One approach for doing this, is direct numerical simulation of pore-scale flow that requires a robust numerical tool for prediction of fluid dynamics and a detailed understanding of the physical processes occurring at the pore-scale. In pore scale flows with a low capillary number, Eulerian multiphase methods are well-known to produce additional vorticity close to the interface. This is mainly due to discretization errors which lead to an imbalance of capillary pressure and surface tension forces that causes unphysical spurious currents. At the pore scale, these spurious currents can become significantly stronger than the average velocity in the phases, and lead to unphysical displacement of the interface. In this work, we first investigate the capability of the algebraic Volume of Fluid (VOF) method in OpenFOAM for low capillary number pore scale flow simulations. Afterward, we compare VOF results with a Coupled Level-Set Volume of Fluid (CLSVOF) method and Iso-Advector method. It has been shown that the former one reduces the VOF's unphysical spurious currents in some cases, and both are known to capture interfaces sharper than VOF. As the conclusion, we will investigate that whether the use of CLSVOF or Iso-Advector will lead to less spurious velocities and more accurate results for capillary driven pore-scale multiphase flows or not. Keywords: Pore-scale multiphase flow, Capillary driven flows, Spurious currents, OpenFOAM

Hydrothermal fluid flow within a tectonically active rift-ridge transform junction: Tjörnes Fracture Zone, Iceland

Science.gov (United States)

Lupi, M.; Geiger, S.; Graham, C. M.

2010-05-01

We investigate the regional fluid flow dynamics in a highly faulted transform area, the Tjörnes Fracture Zone in northern Iceland which is characterized by steep geothermal gradients, hydrothermal activity, and strong seismicity. We simulate fluid flow within the Tjörnes Fracture Zone using a high-resolution model that was based on the available geological and geophysical data and has the aim to represent the complex geological structures and the thermodynamical processes that drive the regional fluid flow in a physically realistic way. Our results show that convective heat flow and mixing of cold and saline seawater with deep hydrothermal fluids controls the large-scale fluid flow. The distribution of faults has a strong influence on the local hydrodynamics by focusing flow around clusters of faults. This explains the nature of isolated upflow zones of hot hydrothermal fluids which are observed in the Tjörnes Fracture Zone. An important emergent characteristic of the regional fluid flow in the Tjörnes Fracture Zone are two separate flow systems: one in the sedimentary basins, comprising more vigorous convection, and one in the crystalline basement, which is dominated by conduction. These two flow systems yield fundamental insight into the connection between regional hydrothermal fluid flow and seismicity because they form the basis of a toggle switch mechanism that is thought to have caused the hydrogeochemical anomalies recorded at Húsavik before and after the 5.8 M earthquake in September 2002.

Mapping flow distortion on oceanographic platforms using computational fluid dynamics

Directory of Open Access Journals (Sweden)

N. O'Sullivan

2013-10-01

Full Text Available Wind speed measurements over the ocean on ships or buoys are affected by flow distortion from the platform and by the anemometer itself. This can lead to errors in direct measurements and the derived parametrisations. Here we computational fluid dynamics (CFD to simulate the errors in wind speed measurements caused by flow distortion on the RV Celtic Explorer. Numerical measurements were obtained from the finite-volume CFD code OpenFOAM, which was used to simulate the velocity fields. This was done over a range of orientations in the test domain from −60 to +60° in increments of 10°. The simulation was also set up for a range of velocities, ranging from 5 to 25 m s−1 in increments of 0.5 m s−1. The numerical analysis showed close agreement to experimental measurements.

Numerical simulation of bubble growth and departure during flow boiling period by lattice Boltzmann method

International Nuclear Information System (INIS)

Sun, Tao; Li, Weizhong; Yang, Shuai

2013-01-01

Highlights: • The bubble departure diameter is proportional to g −0.425 in quiescent fluid. • The bubble release frequency is proportional to g 0.678 in quiescent fluid. • The simulation result supports the transient micro-convection model. • The bubble departure diameter has exponential relation with inlet velocity. • The bubble release frequency has linear relation with inlet velocity. -- Abstract: Nucleate boiling flows on a horizontal plate are studied in this paper by a hybrid lattice Boltzmann method, where both quiescent and slowly flowing ambient are concerned. The process of a single bubble growth on and departure from the superheated wall is simulated. The simulation result supports the transient micro-convection model. The bubble departure diameter and the release frequency are investigated from the simulation result. It is found that the bubble departure diameter and the release frequency are proportional to g −0.425 and g 0.678 in quiescent fluid, respectively, where g is the gravitational acceleration. Nucleate boiling in slowly flowing ambient is also calculated in consideration of forced convection. It is presented that the bubble departure diameter and the release frequency have exponential relationship and linear relationship with inlet velocity in slowly flowing fluid, respectively

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...

Numerical simulation for Jeffery-Hamel flow and heat transfer of micropolar fluid based on differential evolution algorithm

Science.gov (United States)

Ara, Asmat; Khan, Najeeb Alam; Naz, Farah; Raja, Muhammad Asif Zahoor; Rubbab, Qammar

2018-01-01

This article explores the Jeffery-Hamel flow of an incompressible non-Newtonian fluid inside non-parallel walls and observes the influence of heat transfer in the flow field. The fluid is considered to be micropolar fluid that flows in a convergent/divergent channel. The governing nonlinear partial differential equations (PDEs) are converted to nonlinear coupled ordinary differential equations (ODEs) with the help of a suitable similarity transformation. The resulting nonlinear analysis is determined analytically with the utilization of the Taylor optimization method based on differential evolution (DE) algorithm. In order to understand the flow field, the effects of pertinent parameters such as the coupling parameter, spin gradient viscosity parameter and the Reynolds number have been examined on velocity and temperature profiles. It concedes that the good results can be attained by an implementation of the proposed method. Ultimately, the accuracy of the method is confirmed by comparing the present results with the results obtained by Runge-Kutta method.

Numerical simulation for Jeffery-Hamel flow and heat transfer of micropolar fluid based on differential evolution algorithm

Directory of Open Access Journals (Sweden)

Asmat Ara

2018-01-01

Full Text Available This article explores the Jeffery-Hamel flow of an incompressible non-Newtonian fluid inside non-parallel walls and observes the influence of heat transfer in the flow field. The fluid is considered to be micropolar fluid that flows in a convergent/divergent channel. The governing nonlinear partial differential equations (PDEs are converted to nonlinear coupled ordinary differential equations (ODEs with the help of a suitable similarity transformation. The resulting nonlinear analysis is determined analytically with the utilization of the Taylor optimization method based on differential evolution (DE algorithm. In order to understand the flow field, the effects of pertinent parameters such as the coupling parameter, spin gradient viscosity parameter and the Reynolds number have been examined on velocity and temperature profiles. It concedes that the good results can be attained by an implementation of the proposed method. Ultimately, the accuracy of the method is confirmed by comparing the present results with the results obtained by Runge-Kutta method.

Incompressible Turbulent Flow Simulation Using the κ-ɛ Model and Upwind Schemes

Directory of Open Access Journals (Sweden)

V. G. Ferreira

2007-01-01

Full Text Available In the computation of turbulent flows via turbulence modeling, the treatment of the convective terms is a key issue. In the present work, we present a numerical technique for simulating two-dimensional incompressible turbulent flows. In particular, the performance of the high Reynolds κ-ɛ model and a new high-order upwind scheme (adaptative QUICKEST by Kaibara et al. (2005 is assessed for 2D confined and free-surface incompressible turbulent flows. The model equations are solved with the fractional-step projection method in primitive variables. Solutions are obtained by using an adaptation of the front tracking GENSMAC (Tomé and McKee (1994 methodology for calculating fluid flows at high Reynolds numbers. The calculations are performed by using the 2D version of the Freeflow simulation system (Castello et al. (2000. A specific way of implementing wall functions is also tested and assessed. The numerical procedure is tested by solving three fluid flow problems, namely, turbulent flow over a backward-facing step, turbulent boundary layer over a flat plate under zero-pressure gradients, and a turbulent free jet impinging onto a flat surface. The numerical method is then applied to solve the flow of a horizontal jet penetrating a quiescent fluid from an entry port beneath the free surface.

Numerical Simulation on Forced Convective Condensation of Steam Upward Flow in a Vertical Pipe

Directory of Open Access Journals (Sweden)

Guo-Dong Qiu

2014-05-01

Full Text Available A transient three-dimensional volume of fluid (VOF simulation on condensation of upward flow of wet steam inside a 12 mm i.d. vertical pipe is presented. The effect of gravity and surface tension are taken into account. A uniform wall temperature has been fixed as boundary conditions. The mass flux is 130~6400 kg m−2's−1 and the turbulence inside the vapor phase and liquid phase have been handled by Reynolds stress model (RSM. The vapor quality of fluid is 0~0.4. The numerical simulation results show that, in all the simulation conditions, the bubbly flow, slug flow, churn flow, wispy annular flow, and annular flow are observed; in addition, the results of flow pattern are in good agreement with the regime map from Hewitt and Roberts. The typical velocity field characteristic of each flow pattern and the effect of velocity field on heat transfer of condensation are analyzed, indicating that the slug flow and churn flow have obvious local eddy. However, no obvious eddy is observed in other flow patterns and the streamlines are almost parallel to the flow direction. The simulation results of heat transfer coefficients and frictional pressure drop show good agreement with the correlations from existing literatures.

Numerical simulation of random stresses on an annular turbulent flow

International Nuclear Information System (INIS)

Marti-Moreno, Marta

2000-01-01

The flow along a circular cylinder may induce structural vibrations. For the predictive analysis of such vibrations, the turbulent forcing spectrum needs to be characterized. The aim of this work is to study the turbulent fluid forces acting on a single tube in axial flow. More precisely we have performed numerical simulations of an annular flow. These simulations were carried out on a cylindrical staggered mesh by a finite difference method. We consider turbulent flow with Reynolds number up to 10 6 . The Large Eddy Simulation Method has been used. A survey of existent experiments showed that hydraulic diameter acts as an important parameter. We first showed the accuracy of the numerical code by reproducing the experiments of Mulcahy. The agreement between pressure spectra from computations and from experiments is good. Then, we applied this code to simulate new numerical experiments varying the hydraulic diameter and the flow velocity. (author) [fr

Continuous media theory for MR fluids in non-shearing flows

International Nuclear Information System (INIS)

Ruiz-López, J A; Hidalgo-Alvarez, R; Vicente, J de

2013-01-01

The enhanced mechanical response of magnetorheological fluids under slow compression has been investigated by means of experiments, theory and particle-level simulations. A wide range of magnetic field strengths (0–354 kA/m), dispersing medium viscosities (20–500 mPa·s) and particle concentrations (5–30 vol%) were investigated. Plastic media theory in compressive flow was in good agreement with experimental data. Slight deviations from the theory were associated to the so-called strengthening effect as the yield shear stress could increase during compression. Particle-level simulations were in good agreement with both experiments and simulations.

Validation of a two-fluid model used for the simulation of dense fluidized beds; Validation d`un modele a deux fluides applique a la simulation des lits fluidises denses

Energy Technology Data Exchange (ETDEWEB)

Boelle, A.

1997-02-17

A two-fluid model applied to the simulation of gas-solid dense fluidized beds is validated on micro scale and on macro scale. Phase coupling is carried out in the momentum and energy transport equation of both phases. The modeling is built on the kinetic theory of granular media in which the gas action has been taken into account in order to get correct expressions of transport coefficients. A description of hydrodynamic interactions between particles in high Stokes number flow is also incorporated in the model. The micro scale validation uses Lagrangian numerical simulations viewed as numerical experiments. The first validation case refers to a gas particle simple shear flow. It allows to validate the competition between two dissipation mechanisms: drag and particle collisions. The second validation case is concerted with sedimenting particles in high Stokes number flow. It allows to validate our approach of hydrodynamic interactions. This last case had led us to develop an original Lagrangian simulation with a two-way coupling between the fluid and the particles. The macro scale validation uses the results of Eulerian simulations of dense fluidized bed. Bed height, particles circulation and spontaneous created bubbles characteristics are studied and compared to experimental measurement, both looking at physical and numerical parameters. (author) 159 refs.

The Interplay Between Saline Fluid Flow and Dynamic Permeability in Magmatic-Hydrothermal Systems

Science.gov (United States)

Weis, P.

2014-12-01

Magmatic-hydrothermal ore deposits document the interplay between saline fluid flow and rock permeability. Numerical simulations of multi-phase flow of variably miscible, compressible H20-NaCl fluids in concert with a dynamic permeability model can reproduce characteristics of porphyry copper and epithermal gold systems. This dynamic permeability model incorporates depth-dependent permeability profiles characteristic for tectonically active crust as well as pressure- and temperature-dependent relationships describing hydraulic fracturing and the transition from brittle to ductile rock behavior. In response to focused expulsion of magmatic fluids from a crystallizing upper crustal magma chamber, the hydrothermal system self-organizes into a hydrological divide, separating an inner part dominated by ascending magmatic fluids under near-lithostatic pressures from a surrounding outer part dominated by convection of colder meteoric fluids under near-hydrostatic pressures. This hydrological divide also provides a mechanism to transport magmatic salt through the crust, and prevents the hydrothermal system to become "clogged" by precipitation of solid halite due to depressurization of saline, high-temperature magmatic fluids. The same physical processes at similar permeability ranges, crustal depths and flow rates are relevant for a number of active systems, including geothermal resources and excess degassing at volcanos. The simulations further suggest that the described mechanism can separate the base of free convection in high-enthalpy geothermal systems from the magma chamber as a driving heat source by several kilometers in the vertical direction in tectonic settings with hydrous magmatism. This hydrology would be in contrast to settings with anhydrous magmatism, where the base of the geothermal systems may be closer to the magma chamber.

Design and numerical simulation on an auto-cumulative flowmeter in horizontal oil-water two-phase flow.

Science.gov (United States)

Xie, Beibei; Kong, Lingfu; Kong, Deming; Kong, Weihang; Li, Lei; Liu, Xingbin; Chen, Jiliang

2017-11-01

In order to accurately measure the flow rate under the low yield horizontal well conditions, an auto-cumulative flowmeter (ACF) was proposed. Using the proposed flowmeter, the oil flow rate in horizontal oil-water two-phase segregated flow can be finely extracted. The computational fluid dynamics software Fluent was used to simulate the fluid of the ACF in oil-water two-phase flow. In order to calibrate the simulation measurement of the ACF, a novel oil flow rate measurement method was further proposed. The models of the ACF were simulated to obtain and calibrate the oil flow rate under different total flow rates and oil cuts. Using the finite-element method, the structure of the seven conductance probes in the ACF was simulated. The response values for the probes of the ACF under the conditions of oil-water segregated flow were obtained. The experiments for oil-water segregated flow under different heights of the oil accumulation in horizontal oil-water two-phase flow were carried out to calibrate the ACF. The validity of the oil flow rate measurement in horizontal oil-water two-phase flow was verified by simulation and experimental results.

Influence of fluid properties, flow rate and aspect ratios on stratification in a cylindrical cavity

International Nuclear Information System (INIS)

Bouhdjar, A.; Harhad, A.; Guerri, O.

2003-01-01

The fluid flow and temperature field in a cavity are numerically simulated using finite volume techniques. The fluid flow in the vertical cylindrical cavity is assumed to be two-dimensional. Inflow occurs at the top through a ring like entrance and outflow takes place at the bottom through an exit of the same shape. The study considers a transient mixed convection flow. The governing equations are the conservation equations for laminar natural convection flow based on the Boussinesq approximation. Forced convection flow is superimposed through the appropriate boundary conditions (inflow and outflow conditions). The influence of the mass flow rate and of the fluid is made through the Reynolds number and the Prandtl number. Stratification analysis is made qualitatively through temperature distribution. The study considers two fluids i.e. water (Pr=4.5) and ethylene glycol (Pr=51) and cavity aspect ratios of 1/0.5 and 1 /2. So the objective of the work is to get more information on the influence of flow rate on the performance of the thermal energy storage. Correlations for the storage efficiency are deduced with respect to the Reynolds number. (author)

A simple method for potential flow simulation of cascades

Indian Academy of Sciences (India)

vortex panel method to simulate potential flow in cascades is presented. The cascade ... The fluid loading on the blades, such as the normal force and pitching moment, may ... of such discrete infinite array singularities along the blade surface.

Re-injection feasibility study of fracturing flow-back fluid in shale gas mining

Science.gov (United States)

Kang, Dingyu; Xue, Chen; Chen, Xinjian; Du, Jiajia; Shi, Shengwei; Qu, Chengtun; Yu, Tao

2018-02-01

Fracturing flow-back fluid in shale gas mining is usually treated by re-injecting into formation. After treatment, the fracturing flow-back fluid is injected back into the formation. In order to ensure that it will not cause too much damage to the bottom layer, feasibility evaluations of re-injection of two kinds of fracturing fluid with different salinity were researched. The experimental research of the compatibility of mixed water samples based on the static simulation method was conducted. Through the analysis of ion concentration, the amount of scale buildup and clay swelling rate, the feasibility of re-injection of different fracturing fluid were studied. The result shows that the swelling of the clay expansion rate of treated fracturing fluid is lower than the mixed water of treated fracturing fluid and the distilled water, indicating that in terms of clay expansion rate, the treated fracturing flow-back fluid is better than that of water injection after re-injection. In the compatibility test, the maximum amount of fouling in the Yangzhou oilfield is 12mg/L, and the maximum value of calcium loss rate is 1.47%, indicating that the compatibility is good. For the fracturing fluid with high salinity in the Yanchang oilfield, the maximum amount of scaling is 72mg/L, and the maximum calcium loss rate is 3.50%, indicating that the compatibility is better.

Extension of Generalized Fluid System Simulation Program's Fluid Property Database

Science.gov (United States)

Patel, Kishan

2011-01-01

This internship focused on the development of additional capabilities for the General Fluid Systems Simulation Program (GFSSP). GFSSP is a thermo-fluid code used to evaluate system performance by a finite volume-based network analysis method. The program was developed primarily to analyze the complex internal flow of propulsion systems and is capable of solving many problems related to thermodynamics and fluid mechanics. GFSSP is integrated with thermodynamic programs that provide fluid properties for sub-cooled, superheated, and saturation states. For fluids that are not included in the thermodynamic property program, look-up property tables can be provided. The look-up property tables of the current release version can only handle sub-cooled and superheated states. The primary purpose of the internship was to extend the look-up tables to handle saturated states. This involves a) generation of a property table using REFPROP, a thermodynamic property program that is widely used, and b) modifications of the Fortran source code to read in an additional property table containing saturation data for both saturated liquid and saturated vapor states. Also, a method was implemented to calculate the thermodynamic properties of user-fluids within the saturation region, given values of pressure and enthalpy. These additions required new code to be written, and older code had to be adjusted to accommodate the new capabilities. Ultimately, the changes will lead to the incorporation of this new capability in future versions of GFSSP. This paper describes the development and validation of the new capability.

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.

Wind flow simulation over flat terrain using CFD based software

International Nuclear Information System (INIS)

Petrov, Peter; Terziev, Angel; Genovski, Ivan

2009-01-01

Velocity distribution recognition over definite place (terrain) is very important because due to that the zones with high energy potential could be defined (the fields with high velocities). This is a precondition for optimal wind turbine generators micro-sitting. In current work a simulation of the open flow over the flat terrain using the CFD based software is reviewed. The simulations are made of a real fluid flow in order to be defined the velocity fields over the terrain

TOUGH2: A general-purpose numerical simulator for multiphase nonisothermal flows

Energy Technology Data Exchange (ETDEWEB)

Pruess, K. [Lawrence Berkeley Lab., CA (United States)

1991-06-01

Numerical simulators for multiphase fluid and heat flows in permeable media have been under development at Lawrence Berkeley Laboratory for more than 10 yr. Real geofluids contain noncondensible gases and dissolved solids in addition to water, and the desire to model such `compositional` systems led to the development of a flexible multicomponent, multiphase simulation architecture known as MULKOM. The design of MULKOM was based on the recognition that the mass-and energy-balance equations for multiphase fluid and heat flows in multicomponent systems have the same mathematical form, regardless of the number and nature of fluid components and phases present. Application of MULKOM to different fluid mixtures, such as water and air, or water, oil, and gas, is possible by means of appropriate `equation-of-state` (EOS) modules, which provide all thermophysical and transport parameters of the fluid mixture and the permeable medium as a function of a suitable set of primary thermodynamic variables. Investigations of thermal and hydrologic effects from emplacement of heat-generating nuclear wastes into partially water-saturated formations prompted the development and release of a specialized version of MULKOM for nonisothermal flow of water and air, named TOUGH. TOUGH is an acronym for `transport of unsaturated groundwater and heat` and is also an allusion to the tuff formations at Yucca Mountain, Nevada. The TOUGH2 code is intended to supersede TOUGH. It offers all the capabilities of TOUGH and includes a considerably more general subset of MULKOM modules with added capabilities. The paper briefly describes the simulation methodology and user features.

Damping characteristics and flow behaviors of an ER fluid with a piston sine vibration in a viscous damper

International Nuclear Information System (INIS)

Yamaguchi, Hiroshi; Zhang, Xin-Rong; Niu, Xiao-Dong

2010-01-01

The damping characteristics and flow behaviors of ER fluids inside a piston–cylinder viscous damper subjected to external electric fields are studied based on experiment, theoretical analysis and numerical simulation. The viscous damper is a closed system with an inner piston and an outer cylinder, which is designed and constructed in our laboratory. In the experiment, the test ER fluid is enclosed in the gap of a piston–cylinder system. To examine the damping characteristics of the test ER fluid, a piston sine vibration experiment is performed with accompanying theoretical analyses. In addition, in order to investigate the ER flow behaviors inside the damper, a numerical simulation is carried out. The present study discloses the damping characteristics and the fluid mechanism of the ER fluid in the piston–cylinder damper with an applied external electric field

Multiscale Simulation Framework for Coupled Fluid Flow and Mechanical Deformation

Energy Technology Data Exchange (ETDEWEB)

Hou, Thomas [California Inst. of Technology (CalTech), Pasadena, CA (United States); Efendiev, Yalchin [Stanford Univ., CA (United States); Tchelepi, Hamdi [Texas A & M Univ., College Station, TX (United States); Durlofsky, Louis [Stanford Univ., CA (United States)

2016-05-24

Our work in this project is aimed at making fundamental advances in multiscale methods for flow and transport in highly heterogeneous porous media. The main thrust of this research is to develop a systematic multiscale analysis and efficient coarse-scale models that can capture global effects and extend existing multiscale approaches to problems with additional physics and uncertainties. A key emphasis is on problems without an apparent scale separation. Multiscale solution methods are currently under active investigation for the simulation of subsurface flow in heterogeneous formations. These procedures capture the effects of fine-scale permeability variations through the calculation of specialized coarse-scale basis functions. Most of the multiscale techniques presented to date employ localization approximations in the calculation of these basis functions. For some highly correlated (e.g., channelized) formations, however, global effects are important and these may need to be incorporated into the multiscale basis functions. Other challenging issues facing multiscale simulations are the extension of existing multiscale techniques to problems with additional physics, such as compressibility, capillary effects, etc. In our project, we explore the improvement of multiscale methods through the incorporation of additional (single-phase flow) information and the development of a general multiscale framework for flows in the presence of uncertainties, compressible flow and heterogeneous transport, and geomechanics. We have considered (1) adaptive local-global multiscale methods, (2) multiscale methods for the transport equation, (3) operator-based multiscale methods and solvers, (4) multiscale methods in the presence of uncertainties and applications, (5) multiscale finite element methods for high contrast porous media and their generalizations, and (6) multiscale methods for geomechanics.

CFD Numerical Simulation of the Complex Turbulent Flow Field in an Axial-Flow Water Pump

Directory of Open Access Journals (Sweden)

Wan-You Li

2014-09-01

Full Text Available Further optimal design of an axial-flow water pump calls for a thorough recognition of the characteristics of the complex turbulent flow field in the pump, which is however extremely difficult to be measured using the up-to-date experimental techniques. In this study, a numerical simulation procedure based on computational fluid dynamics (CFD was elaborated in order to obtain the fully three-dimensional unsteady turbulent flow field in an axial-flow water pump. The shear stress transport (SST k-ω model was employed in the CFD calculation to study the unsteady internal flow of the axial-flow pump. Upon the numerical simulation results, the characteristics of the velocity field and pressure field inside the impeller region were discussed in detail. The established model procedure in this study may provide guidance to the numerical simulations of turbomachines during the design phase or the investigation of flow and pressure field characteristics and performance. The presented information can be of reference value in further optimal design of the axial-flow pump.

Interpretation SP anomaly caused by subsurface fluid flow; Chika ryudokei ni yoru shizen den`i anomaly no kaishaku ni tsuite

Energy Technology Data Exchange (ETDEWEB)

Yasukawa, K [Geological Survey of Japan, Tsukuba (Japan); Mogi, T; Kawahira, M [Kyushu University, Fukuoka (Japan). Faculty of Engineering

1996-05-01

A PTSP simulator, developed by combining a fluid flow simulator and a simulator for computing the spontaneous potential (SP) out of flow velocity distribution, was used to investigate the impact of the distribution of terrain, permeability, and resistivity upon the ground surface SP. SP computation using simple models such as highlands or cliffs indicated that the terrain-caused SP anomaly was negatively dependent upon elevation. It was also indicated, however, that, in some types of resistivity distribution, the SP profile changes, the peak goes out of place, and even the apparent polarity may reverse in extreme cases. In the study of the SP profile for the Takeyu hot spa, Oita Prefecture, PTSP-aided modeling was carried out. It was then found that fluid flows caused by the terrain were not enough to explain the peak at the middle of the mountain, which suggested the existence of a fluid flow caused by temperature distribution. 11 refs., 5 figs.

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

Simulation of reactive polydisperse sprays strongly coupled to unsteady flows in solid rocket motors: Efficient strategy using Eulerian Multi-Fluid methods

Science.gov (United States)

Sibra, A.; Dupays, J.; Murrone, A.; Laurent, F.; Massot, M.

2017-06-01

In this paper, we tackle the issue of the accurate simulation of evaporating and reactive polydisperse sprays strongly coupled to unsteady gaseous flows. In solid propulsion, aluminum particles are included in the propellant to improve the global performances but the distributed combustion of these droplets in the chamber is suspected to be a driving mechanism of hydrodynamic and acoustic instabilities. The faithful prediction of two-phase interactions is a determining step for future solid rocket motor optimization. When looking at saving computational ressources as required for industrial applications, performing reliable simulations of two-phase flow instabilities appears as a challenge for both modeling and scientific computing. The size polydispersity, which conditions the droplet dynamics, is a key parameter that has to be accounted for. For moderately dense sprays, a kinetic approach based on a statistical point of view is particularly appropriate. The spray is described by a number density function and its evolution follows a Williams-Boltzmann transport equation. To solve it, we use Eulerian Multi-Fluid methods, based on a continuous discretization of the size phase space into sections, which offer an accurate treatment of the polydispersion. The objective of this paper is threefold: first to derive a new Two Size Moment Multi-Fluid model that is able to tackle evaporating polydisperse sprays at low cost while accurately describing the main driving mechanisms, second to develop a dedicated evaporation scheme to treat simultaneously mass, moment and energy exchanges with the gas and between the sections. Finally, to design a time splitting operator strategy respecting both reactive two-phase flow physics and cost/accuracy ratio required for industrial computations. Using a research code, we provide 0D validations of the new scheme before assessing the splitting technique's ability on a reference two-phase flow acoustic case. Implemented in the industrial

Fluid Dynamics of Magnetic Nanoparticles in Simulated Blood Vessels

Science.gov (United States)

Blue, Lauren; Sewell, Mary Kathryn; Brazel, Christopher S.

2008-11-01

Magnetic nanoparticles (MNPs) can be used to locally target therapies and offer the benefit of using an AC magnetic field to combine hyperthermia treatment with the triggered release of therapeutic agents. Here, we investigate localization of MNPs in a simulated environment to understand the relationship between magnetic field intensity and bulk fluid dynamics to determine MNP retention in a simulated blood vessel. As MNPs travel through blood vessels, they can be slowed or trapped in a specific area by applying a magnetic field. Magnetic cobalt ferrite nanoparticles were synthesized and labeled with a fluorescent rhodamine tag to visualize patterns in a flow cell, as monitored by a fluorescence microscope. Particle retention was determined as a function of flow rate, concentration, and magnetic field strength. Understanding the relationship between magnetic field intensity, flow behavior and nanoparticle characteristics will aid in the development of therapeutic systems specifically targeted to diseased tissue.

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...

Improving flow distribution in influent channels using computational fluid dynamics.

Science.gov (United States)

Park, No-Suk; Yoon, Sukmin; Jeong, Woochang; Lee, Seungjae

2016-10-01

Although the flow distribution in an influent channel where the inflow is split into each treatment process in a wastewater treatment plant greatly affects the efficiency of the process, and a weir is the typical structure for the flow distribution, to the authors' knowledge, there is a paucity of research on the flow distribution in an open channel with a weir. In this study, the influent channel of a real-scale wastewater treatment plant was used, installing a suppressed rectangular weir that has a horizontal crest to cross the full channel width. The flow distribution in the influent channel was analyzed using a validated computational fluid dynamics model to investigate (1) the comparison of single-phase and two-phase simulation, (2) the improved procedure of the prototype channel, and (3) the effect of the inflow rate on flow distribution. The results show that two-phase simulation is more reliable due to the description of the free-surface fluctuations. It should first be considered for improving flow distribution to prevent a short-circuit flow, and the difference in the kinetic energy with the inflow rate makes flow distribution trends different. The authors believe that this case study is helpful for improving flow distribution in an influent channel.

Simulation of two-phase flows in vertical tubes with the CTFD code FLUBOX

International Nuclear Information System (INIS)

Graf, Udo; Papadimitriou, Pavlos

2007-01-01

The computational two-fluid dynamics (CTFD) code FLUBOX is developed at GRS for the multidimensional simulation of two-phase flows. The single-pressure two-fluid model is used as basis of the simulation. A basic mathematical property of the two-fluid model of FLUBOX is the hyperbolic character of the advection. The numerical solution methods of FLUBOX make explicit use of the hyperbolic structure of the coefficient matrices. The simulation of two-phase flow phenomena needs, apart from the conservation equations for each phase, an additional transport equation for the interfacial area concentration. The concentration of the interfacial area is one of the key parameters for the modeling of interfacial friction forces and interfacial transfer terms. A new transport equation for the interfacial area concentration is in development. It describes the dynamic change of the interfacial area concentration due to mass exchange and a force balance at the phase boundary. Results from FLUBOX calculations for different experiments of two-phase flows in vertical tubes are presented as part of the validation

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.

Simulation of flooding waves in vertical churn flow

Energy Technology Data Exchange (ETDEWEB)

Tekavčič, Matej, E-mail: matej.tekavcic@ijs.si; Končar, Boštjan; Kljenak, Ivo

2016-04-01

Highlights: • Flooding waves in air–water churn flow in a vertical pipe were studied. • Simulations using two-fluid model with interface sharpening were performed. • Calculated wave amplitudes agree with existing experimental data. • Contributions of force terms in the liquid momentum balance equation are presented. - Abstract: A transient simulation of flooding waves in the churn flow of air and water in a vertical pipe is performed by the means of two-fluid modelling approach with interface sharpening. The gas and liquid phases are considered immiscible and incompressible with no mass transfer between them. Inter-phase coupling of momentum is realized via interface drag force which is based on the interface area density and the relative velocity between the phases. Surface tension effects are modelled with the Continuum Surface Model. The flow is assumed isothermal. Turbulence is modelled for each phase separately using the two-equation eddy viscosity approach. Results are compared with the reported experimental data for churn flow regime in a vertical pipe (Wang et al., 2011a). Reynolds numbers of the gas flow are in the range from 6000 to 10,000, while the liquid mass flow rate upwards ranges from 25 to 32 g/s. Prediction of critical and maximum amplitudes of the flooding waves show good agreement with experimental values. Results for wave frequencies indicate significant deviations, which can be attributed to the choice of the liquid inlet model.

Simulation of flooding waves in vertical churn flow

International Nuclear Information System (INIS)

Tekavčič, Matej; Končar, Boštjan; Kljenak, Ivo

2016-01-01

Highlights: • Flooding waves in air–water churn flow in a vertical pipe were studied. • Simulations using two-fluid model with interface sharpening were performed. • Calculated wave amplitudes agree with existing experimental data. • Contributions of force terms in the liquid momentum balance equation are presented. - Abstract: A transient simulation of flooding waves in the churn flow of air and water in a vertical pipe is performed by the means of two-fluid modelling approach with interface sharpening. The gas and liquid phases are considered immiscible and incompressible with no mass transfer between them. Inter-phase coupling of momentum is realized via interface drag force which is based on the interface area density and the relative velocity between the phases. Surface tension effects are modelled with the Continuum Surface Model. The flow is assumed isothermal. Turbulence is modelled for each phase separately using the two-equation eddy viscosity approach. Results are compared with the reported experimental data for churn flow regime in a vertical pipe (Wang et al., 2011a). Reynolds numbers of the gas flow are in the range from 6000 to 10,000, while the liquid mass flow rate upwards ranges from 25 to 32 g/s. Prediction of critical and maximum amplitudes of the flooding waves show good agreement with experimental values. Results for wave frequencies indicate significant deviations, which can be attributed to the choice of the liquid inlet model.

Computational Fluid Dynamics (CFD) simulations of a Heisenberg Vortex Tube

Science.gov (United States)

Bunge, Carl; Sitaraman, Hariswaran; Leachman, Jake

2017-11-01

A 3D Computational Fluid Dynamics (CFD) simulation of a Heisenberg Vortex Tube (HVT) is performed to estimate cooling potential with cryogenic hydrogen. The main mechanism driving operation of the vortex tube is the use of fluid power for enthalpy streaming in a highly turbulent swirl in a dual-outlet tube. This enthalpy streaming creates a temperature separation between the outer and inner regions of the flow. Use of a catalyst on the peripheral wall of the centrifuge enables endothermic conversion of para-ortho hydrogen to aid primary cooling. A κ- ɛ turbulence model is used with a cryogenic, non-ideal equation of state, and para-orthohydrogen species evolution. The simulations are validated with experiments and strategies for parametric optimization of this device are presented.

Transient Simulation of Accumulating Particle Deposition in Pipe Flow

Science.gov (United States)

Hewett, James; Sellier, Mathieu

2015-11-01

Colloidal particles that deposit in pipe systems can lead to fouling which is an expensive problem in both the geothermal and oil & gas industries. We investigate the gradual accumulation of deposited colloids in pipe flow using numerical simulations. An Euler-Lagrangian approach is employed for modelling the fluid and particle phases. Particle transport to the pipe wall is modelled with Brownian motion and turbulent diffusion. A two-way coupling exists between the fouled material and the pipe flow; the local mass flux of depositing particles is affected by the surrounding fluid in the near-wall region. This coupling is modelled by changing the cells from fluid to solid as the deposited particles exceed each local cell volume. A similar method has been used to model fouling in engine exhaust systems (Paz et al., Heat Transfer Eng., 34(8-9):674-682, 2013). We compare our deposition velocities and deposition profiles with an experiment on silica scaling in turbulent pipe flow (Kokhanenko et al., 19th AFMC, 2014).

Etude d'un modele de Boltzmann sur reseau pour la simulation assistee par ordinateur des fluides a plusieurs phases immiscibles

Science.gov (United States)

Leclaire, Sebastien

The computer assisted simulation of the dynamics of fluid flow has been a highly rewarding topic of research for several decades now, in terms of the number of scientific problems that have been solved as a result, both in the academic world and in industry. In the fluid dynamics field, simulating multiphase immiscible fluid flow remains a challenge, because of the complexity of the interactions at the flow phase interfaces. Various numerical methods are available to study these phenomena, and, the lattice Boltzmann method has been shown in recent years to be well adapted to solving this type of complex flow. In this thesis, a lattice Boltzmann model for the simulation of two-phase immiscible flows is studied. The main objective of the thesis is to develop this promising method further, with a view to enhancing its validity. To achieve this objective, the research is divided into five distinct themes. The first two focus on correcting some of the deficiencies of the original model. The third generalizes the model to support the simulation of N-phase immiscible fluid flows. The fourth is aimed at modifying the model itself, to enable the simulation of immiscible fluid flows in which the density of the phases varies. With the lattice Boltzmann class of models studied here, this density variation has been inadequately modeled, and, after 20 years, the issue still has not been resolved. The fifth, which complements this thesis, is connected with the lattice Boltzmann method, in that it generalizes the theory of 2D and 3D isotropic gradients for a high order of spatial precision. These themes have each been the subject of a scientific article, as listed in the appendix to this thesis, and together they constitute a synthesis that explains the links between the articles, as well as their scientific contributions, and satisfy the main objective of this research. Globally, a number of qualitative and quantitative test cases based on the theory of multiphase fluid flows

Investigation of entrance length in circular and noncircular conduits by computational fluid dynamics simulation

Directory of Open Access Journals (Sweden)

Pimpun Tongpun

2014-08-01

Full Text Available This study estimated entrance length of circular and noncircular conduits, including circle, triangle, square and hexagon cross-sectional conduit, by using computational fluid dynamics (CFD. For simulation condition, the length of noncircular conduit was 10 m and the hydraulic diameter was 0.2 m. The laminar flow with Reynolds number of 500 and turbulent flow with Reynolds number of 50,000 were applied to investigate water flow in conduits. The governing equations were solved iteratively by using ANSYS FLUENT 14.0. For turbulent flow simulation, standard k-epsilon and RNG k-epsilon model were employed to simulate turbulence. The preliminary results were validated by comparison with theoretical data. At first, grid independency was evaluated to optimize the model. Norm* was employed to investigate the entrance length, which is related to velocity. The simulated results revealed that the entrance length for laminar flow was longer than turbulent flow.

Simulation of Sound Waves Using the Lattice Boltzmann Method for Fluid Flow: Benchmark Cases for Outdoor Sound Propagation.

Science.gov (United States)

Salomons, Erik M; Lohman, Walter J A; Zhou, Han

2016-01-01

Propagation of sound waves in air can be considered as a special case of fluid dynamics. Consequently, the lattice Boltzmann method (LBM) for fluid flow can be used for simulating sound propagation. In this article application of the LBM to sound propagation is illustrated for various cases: free-field propagation, propagation over porous and non-porous ground, propagation over a noise barrier, and propagation in an atmosphere with wind. LBM results are compared with solutions of the equations of acoustics. It is found that the LBM works well for sound waves, but dissipation of sound waves with the LBM is generally much larger than real dissipation of sound waves in air. To circumvent this problem it is proposed here to use the LBM for assessing the excess sound level, i.e. the difference between the sound level and the free-field sound level. The effect of dissipation on the excess sound level is much smaller than the effect on the sound level, so the LBM can be used to estimate the excess sound level for a non-dissipative atmosphere, which is a useful quantity in atmospheric acoustics. To reduce dissipation in an LBM simulation two approaches are considered: i) reduction of the kinematic viscosity and ii) reduction of the lattice spacing.

STEADY-STATE modeling and simulation of pipeline networks for compressible fluids

Directory of Open Access Journals (Sweden)

A.L.H. Costa

1998-12-01

Full Text Available This paper presents a model and an algorithm for the simulation of pipeline networks with compressible fluids. The model can predict pressures, flow rates, temperatures and gas compositions at any point of the network. Any network configuration can be simulated; the existence of cycles is not an obstacle. Numerical results from simulated data on a proposed network are shown for illustration. The potential of the simulator is explored by the analysis of a pressure relief network, using a stochastic procedure for the evaluation of system performance.

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

SAFSIM theory manual: A computer program for the engineering simulation of flow systems

Energy Technology Data Exchange (ETDEWEB)

Dobranich, D.

1993-12-01

SAFSIM (System Analysis Flow SIMulator) is a FORTRAN computer program for simulating the integrated performance of complex flow systems. SAFSIM provides sufficient versatility to allow the engineering simulation of almost any system, from a backyard sprinkler system to a clustered nuclear reactor propulsion system. In addition to versatility, speed and robustness are primary SAFSIM development goals. SAFSIM contains three basic physics modules: (1) a fluid mechanics module with flow network capability; (2) a structure heat transfer module with multiple convection and radiation exchange surface capability; and (3) a point reactor dynamics module with reactivity feedback and decay heat capability. Any or all of the physics modules can be implemented, as the problem dictates. SAFSIM can be used for compressible and incompressible, single-phase, multicomponent flow systems. Both the fluid mechanics and structure heat transfer modules employ a one-dimensional finite element modeling approach. This document contains a description of the theory incorporated in SAFSIM, including the governing equations, the numerical methods, and the overall system solution strategies.

Numerical simulation of real-world flows

Energy Technology Data Exchange (ETDEWEB)

Hayase, Toshiyuki, E-mail: hayase@ifs.tohoku.ac.jp [Institute of Fluid Science, Tohoku University, 2-1-1 Katahira, Aoba-ku, Sendai, 980-8577 (Japan)

2015-10-15

Obtaining real flow information is important in various fields, but is a difficult issue because measurement data are usually limited in time and space, and computational results usually do not represent the exact state of real flows. Problems inherent in the realization of numerical simulation of real-world flows include the difficulty in representing exact initial and boundary conditions and the difficulty in representing unstable flow characteristics. This article reviews studies dealing with these problems. First, an overview of basic flow measurement methodologies and measurement data interpolation/approximation techniques is presented. Then, studies on methods of integrating numerical simulation and measurement, namely, four-dimensional variational data assimilation (4D-Var), Kalman filters (KFs), state observers, etc are discussed. The first problem is properly solved by these integration methodologies. The second problem can be partially solved with 4D-Var in which only initial and boundary conditions are control parameters. If an appropriate control parameter capable of modifying the dynamical structure of the model is included in the formulation of 4D-Var, unstable modes are properly suppressed and the second problem is solved. The state observer and KFs also solve the second problem by modifying mathematical models to stabilize the unstable modes of the original dynamical system by applying feedback signals. These integration methodologies are now applied in simulation of real-world flows in a wide variety of research fields. Examples are presented for basic fluid dynamics and applications in meteorology, aerospace, medicine, etc. (topical review)

Pore-scale simulation of fluid flow and solute dispersion in three-dimensional porous media

KAUST Repository

Icardi, Matteo; Boccardo, Gianluca; Marchisio, Daniele L.; Tosco, Tiziana; Sethi, Rajandrea

2014-01-01

In the present work fluid flow and solute transport through porous media are described by solving the governing equations at the pore scale with finite-volume discretization. Instead of solving the simplified Stokes equation (very often employed

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.

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.

Flow-induced vibration analysis of a helical coil steam generator experiment using large eddy simulation

Energy Technology Data Exchange (ETDEWEB)

Yuan, Haomin; Solberg, Jerome; Merzari, Elia; Kraus, Adam; Grindeanu, Iulian

2017-10-01

This paper describes a numerical study of flow-induced vibration in a helical coil steam generator experiment conducted at Argonne National Laboratory in the 1980s. In the experiment, a half-scale sector model of a steam generator helical coil tube bank was subjected to still and flowing air and water, and the vibrational characteristics were recorded. The research detailed in this document utilizes the multi-physics simulation toolkit SHARP developed at Argonne National Laboratory, in cooperation with Lawrence Livermore National Laboratory, to simulate the experiment. SHARP uses the spectral element code Nek5000 for fluid dynamics analysis and the finite element code DIABLO for structural analysis. The flow around the coil tubes is modeled in Nek5000 by using a large eddy simulation turbulence model. Transient pressure data on the tube surfaces is sampled and transferred to DIABLO for the structural simulation. The structural response is simulated in DIABLO via an implicit time-marching algorithm and a combination of continuum elements and structural shells. Tube vibration data (acceleration and frequency) are sampled and compared with the experimental data. Currently, only one-way coupling is used, which means that pressure loads from the fluid simulation are transferred to the structural simulation but the resulting structural displacements are not fed back to the fluid simulation

Lattice Boltzmann simulation for temperature-sensitive magnetic fluids in a porous square cavity

International Nuclear Information System (INIS)

Jin Licong; Zhang Xinrong; Niu Xiaodong

2012-01-01

A lattice Boltzmann method is developed to simulate temperature-sensitive magnetic fluids in a porous cavity. In the simulation, the magnetic force, efficient gravity, viscous loss term and geometric loss term in porous medium are imported to the momentum equation. To test the reliability of the method, a validation with water in porous cavity is carried out. Good agreements with the previous results verify that the present lattice Boltzmann method is promising for simulation of magnetic fluids in porous medium. In this study, we investigate the change of magnetization with external magnetic field, and we present numerical results for the streamlines, isotherms, and magnetization at vertical or horizontal mid-profiles for different values of Ram. In addition, Nusselt numbers changing with magnetic Rayleigh numbers are also investigated. - Highlights: → Developed a lattice Boltzmann method for magnetic nano-fluids in porous cavity. → Clarified flow and heat transfer for different values of (magnetic) Rayleigh numbers. → Heat transfer enhancement for magnetic fluid in porous cavity.

A Two-Phase Flow Solver for Incompressible Viscous Fluids, Using a Pure Streamfunction Formulation and the Volume of Fluid Technique

DEFF Research Database (Denmark)

Comminal, Raphaël; Spangenberg, Jon; Hattel, Jesper Henri

Accurate multi-phase flow solvers at low Reynolds number are of particular interest for the simulation of interface instabilities in the co-processing of multilayered material. We present a two-phase flow solver for incompressible viscous fluids which uses the streamfunction as the primary variable...... of the flow. Contrary to fractional step methods, the streamfunction formulation eliminates the pressure unknowns, and automatically fulfills the incompressibility constraint by construction. As a result, the method circumvents the loss of temporal accuracy at low Reynolds numbers. The interface is tracked...

A Two-Phase Flow Solver for Incompressible Viscous Fluids, Using a Pure Streamfunction Formulation and the Volume of Fluid Technique

DEFF Research Database (Denmark)

Comminal, Raphaël; Spangenberg, Jon; Hattel, Jesper Henri

2014-01-01

Accurate multi-phase flow solvers at low Reynolds number are of particular interest for the simulation of interface instabilities in the co-processing of multilayered material. We present a two-phase flow solver for incompressible viscous fluids which uses the streamfunction as the primary variable...... of the flow. Contrary to fractional step methods, the streamfunction formulation eliminates the pressure unknowns, and automatically fulfills the incompressibility constraint by construction. As a result, the method circumvents the loss of temporal accuracy at low Reynolds numbers. The interface is tracked...

Influence of fluid properties, flow rate and aspect ratios on stratification in a cylindrical cavity

International Nuclear Information System (INIS)

Bouhdjar, A.; Benyoucef, B.; Harhad, A.

2005-01-01

Fluid flow and temperature field in a cavity are numerically simulated using finite volume techniques. The fluid flow in the vertical cylindrical cavity is assumed to be two-dimensional. Inflow occurs at the top through a ring like entrance and outflow takes place at the bottom through an exit of the same shape. The study considers a transient mixed convection flow. The governing equations are the conservation equations for laminar natural convection flow based on the Boussinesq approximation. Forced convection flow is superimposed through the appropriate boundary conditions (inflow and outflow conditions). The influence of the mass flow rate and of the fluid is made through the Reynolds number and the Prandtl number. Stratification analysis is made qualitatively through temperature distribution. In a previous study, consideration was given to low Reynolds numbers i.e. Re +4 ) in considering water (Pr=3.01) as the working fluid for the thermal energy storage. Correlations for the storage efficiency are deduced with respect to the Reynolds number and cavity aspect ratios of 1/0.5, 1/1 and 1/2. So the objective of the work is to get more information on the influence of flow rate on the storage efficiency as well as on the medium mean temperature. (author)

Numerical simulation of steady cavitating flow of viscous fluid in a Francis hydroturbine

Science.gov (United States)

Panov, L. V.; Chirkov, D. V.; Cherny, S. G.; Pylev, I. M.; Sotnikov, A. A.

2012-09-01

Numerical technique was developed for simulation of cavitating flows through the flow passage of a hydraulic turbine. The technique is based on solution of steady 3D Navier—Stokes equations with a liquid phase transfer equation. The approch for setting boundary conditions meeting the requirements of cavitation testing standard was suggested. Four different models of evaporation and condensation were compared. Numerical simulations for turbines of different specific speed were compared with experiment.

Particle simulation algorithms with short-range forces in MHD and fluid flow

International Nuclear Information System (INIS)

Cable, S.; Tajima, T.; Umegaki, K.

1992-07-01

Attempts are made to develop numerical algorithms for handling fluid flows involving liquids and liquid-gas mixtures. In these types of systems, the short-range intermolecular interactions are important enough to significantly alter behavior predicted on the basis of standard fluid mechanics and magnetohydrodynamics alone. We have constructed a particle-in-cell (PIC) code for the purpose of studying the effects of these interactions. Of the algorithms considered, the one which has been successfully implemented is based on a MHD particle code developed by Brunel et al. In the version presented here, short range forces are included in particle motion by, first, calculating the forces between individual particles and then, to prevent aliasing, interpolating these forces to the computational grid points, then interpolating the forces back to the particles. The code has been used to model a simple two-fluid Rayleigh-Taylor instability. Limitations to the accuracy of the code exist at short wavelengths, where the effects of the short-range forces would be expected to be most pronounced

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)

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...

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.

Fluid Simulation in the Movies: Navier and Stokes Must Be Circulating in Their Graves

Science.gov (United States)

Tessendorf, Jerry

2010-11-01

Fluid simulations based on the Incompressible Navier-Stokes equations are commonplace computer graphics tools in the visual effects industry. These simulations mostly come from custom C++ code written by the visual effects companies. Their significant impact in films was recognized in 2008 with Academy Awards to four visual effects companies for their technical achievement. However artists are not fluid dynamicists, and fluid dynamics simulations are expensive to use in a deadline-driven production environment. As a result, the simulation algorithms are modified to limit the computational resources, adapt them to production workflow, and to respect the client's vision of the film plot. Eulerian solvers on fixed rectangular grids use a mix of momentum solvers, including Semi-Lagrangian, FLIP, and QUICK. Incompressibility is enforced with FFT, Conjugate Gradient, and Multigrid methods. For liquids, a levelset field tracks the free surface. Smooth Particle Hydrodynamics is also used, and is part of a hybrid Eulerian-SPH liquid simulator. Artists use all of them in a mix and match fashion to control the appearance of the simulation. Specially designed forces and boundary conditions control the flow. The simulation can be an input to artistically driven procedural particle simulations that enhance the flow with more detail and drama. Post-simulation processing increases the visual detail beyond the grid resolution. Ultimately, iterative simulation methods that fit naturally in the production workflow are extremely desirable but not yet successful. Results from some efforts for iterative methods are shown, and other approaches motivated by the history of production are proposed.

Numerical study on the heat transfer performance of non-Newtonian fluid flow in a manifold microchannel heat sink

International Nuclear Information System (INIS)

Li, Si-Ning; Zhang, Hong-Na; Li, Xiao-Bin; Li, Qian; Li, Feng-Chen; Qian, Shizhi; Joo, Sang Woo

2017-01-01

Highlights: • Heat transfer performance of non-Newtonian fluid flow in a MHS is studied. • Pseudo-plastic fluid flow can clearly promote the heat transfer efficiency in MMC. • Heat transfer enhancement is attributed to the emergence of secondary flow. • The heat transfer uniformity can also be improved by pseudo-plastic fluid flow. - Abstract: As the miniaturization and integration become the leading trend of the micro-electro-mechanical systems, it is of great significance to improve the microscaled heat transfer performance. This paper presents a three-dimensional (3D) numerical simulation on the flow characteristics and heat transfer performance of non-Newtonian fluid flow in a manifold microchannel (MMC) heat sink and traditional microchannel (TMC) heat sink. The non-Newtonian fluid was described by the power-law model. The analyses concentrated on the non-Newtonian fluid effect on the heat transfer performance, including the heat transfer efficiency and uniformity of temperature distribution, as well as the influence of inlet/outlet configurations on fluid flow and heat transfer. Comparing with Newtonian fluid flow, pseudo-plastic fluid could reduce the drag resistance in both MMC and TMC, while the dilatant fluid brought in quite larger drag resistance. For the heat transfer performance, the introduction of pseudo-plastic fluid flow greatly improved the heat transfer efficiency owing to the generation of secondary flow due to the shear-thinning property. Besides, the temperature distribution in MMC was more uniform by using pseudo-plastic fluid. Moreover, the inlet/outlet configuration was also important for the design and arrangement of microchannel heat sinks, since the present work showed that the maximum temperature was prone to locating in the corners near the inlet and outlet. This work provides guidance for optimal design of small-scale heat transfer devices in many cooling applications, such as biomedical chips, electronic systems, and

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)