首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
The steady flow inside a spatially modulated channel is examined for shear‐thinning and shear‐thickening fluids. The flow is induced by the translation of the lower plate. The modulation amplitude is assumed to be small. A regular perturbation expansion of the flow field is used, coupled to a variable‐step finite‐difference scheme, to solve the problem. Convergence and accuracy assessment against finite‐volume calculations indicates that there is a significant range of validity of the perturbation approach. The influence of the wall geometry, inertia and non‐Newtonian effects are investigated systematically. In particular, the influence of the flow and fluid parameters is examined on the conditions for the onset of separation, vortex size and location. Copyright © 2005 John Wiley & Sons, Ltd.  相似文献   

2.
A new algorithm, which combines the spectral element method with elastic viscous splitting stress (EVSS) method, has been developed for viscoelastic fluid flows in a planar contraction channel. The system of spectral element approximations to the velocity, pressure, extra stress and the rate of deformation variables is solved by a preconditioned conjugate gradient method based on the Uzawa iteration procedure. The numerical approach is implemented on a planar four‐to‐one contraction channel for a fluid governed by an Oldroyd‐B constitutive equation. The behaviour of the Oldroyd‐B fluids in the contraction channel is investigated with various Weissenberg numbers. It is shown that numerical solutions obtained here agree well with experimental measurements and other numerical predictions. Copyright © 2003 John Wiley & Sons, Ltd.  相似文献   

3.
This paper reports on a numerical algorithm for the steady flow of viscoelastic fluid. The conservative and constitutive equations are solved using the finite volume method (FVM) with a hybrid scheme for the velocities and first‐order upwind approximation for the viscoelastic stress. A non‐uniform staggered grid system is used. The iterative SIMPLE algorithm is employed to relax the coupled momentum and continuity equations. The non‐linear algebraic equations over the flow domain are solved iteratively by the symmetrical coupled Gauss–Seidel (SCGS) method. In both, the full approximation storage (FAS) multigrid algorithm is used. An Oldroyd‐B fluid model was selected for the calculation. Results are reported for planar 4:1 abrupt contraction at various Weissenberg numbers. The solutions are found to be stable and smooth. The solutions show that at high Weissenberg number the domain must be long enough. The convergence of the method has been verified with grid refinement. All the calculations have been performed on a PC equipped with a Pentium III processor at 550 MHz. Copyright © 2001 John Wiley & Sons, Ltd.  相似文献   

4.
In this paper, a numerical method, which is about the coupling of continuous and discontinuous Galerkin method based on the splitting scheme, is presented for the calculation of viscoelastic flows of the Oldroyd‐B fluid. The momentum equation is discretized in time by using the Adams‐Bashforth second‐order algorithm, and then decoupled via the splitting approach. Considering the Oldroyd‐B constitutive equation, the second‐order Runge‐Kutta approach is selected to complete the temporal discretization. As for the spatial discretizations, the fundamental purpose is to make the best of finite element method (FEM) and discontinuous Galerkin (DG) method to handle different types of equations. Specifically speaking, for the subequations, FEM is chosen to treat the Poisson and Helmholtz equations, and DG is employed to deal with the nonlinear convective term. In addition, because of the hyperbolic nature, DG is also utilized to discretize the Oldroyd‐B constitutive equation spatially. This coupled method avoids resorting to extra stabilization technique occurred in standard FEM framework even for moderately high values of Weissenberg number and also reduces the complexity compared with unified DG scheme. The Oldroyd‐B model is applied to investigate several typical and challenging benchmarks, such as the 4:1 planar contraction flow and the lid‐driven cavity flow, with a wide range of Weissenberg number to illustrate the feasibility, robustness, and validity of our coupled method.  相似文献   

5.
I.Introducti0nLuiCiqunandHuangJunqiI'I(l989),ZhuWeihuiandLuiCiquri1'l(l992)sequentiallystudiedtheaxialflowofsecondorderandMaxwellfluidsandanalyzedtheflowcharactersofthesefluids.Thispaperwillstudyunsteadyrotat0ryflowofsecondorderandMaxwellfluidsinannularpi…  相似文献   

6.
A finite volume, time‐marching for solving time‐dependent viscoelastic flow in two space dimensions for Oldroyd‐B and Phan Thien–Tanner fluids, is presented. A non‐uniform staggered grid system is used. The conservation and constitutive equations are solved using the finite volume method with an upwind scheme for the viscoelastic stresses and an hybrid scheme for the velocities. To calculate the pressure field, the semi‐implicit method for the pressure linked equation revised method is used. The discretized equations are solved sequentially, using the tridiagonal matrix algorithm solver with under‐relaxation. In both, the full approximation storage multigrid algorithm is used to speed up the convergence rate. Simulations of viscoelastic flows in four‐to‐one abrupt plane contraction are carried out. We will study the behaviour at the entrance corner of the four‐to‐one planar abrupt contraction. Using this solver, we show convergence up to a Weissenberg number We of 20 for the Oldroyd‐B model. No limiting Weissenberg number is observed even though a Phan Thien–Tanner model is used. Several numerical results are presented. Smooth and stable solutions are obtained for high Weissenberg number. Copyright © 2002 John Wiley & Sons, Ltd.  相似文献   

7.
The present investigation deals with the three‐dimensional flow of an Oldroyd‐B fluid over a stretching surface. The governing equations for the three‐dimensional flow are developed. Similarity transformations are invoked for the conversion of nonlinear partial differential equations into the coupled system of ordinary differential equations. Computations for the series solution are presented through implementation of homotopy analysis method. The salient features of the involved parameters have been pointed out. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

8.
A robust finite volume method for viscoelastic flow analysis on general unstructured meshes is developed. It is built upon a general‐purpose stabilization framework for high Weissenberg number flows. The numerical framework provides full combinatorial flexibility between different kinds of rheological models on the one hand, and effective stabilization methods on the other hand. A special emphasis is put on the velocity‐stress‐coupling on colocated computational grids. Using special face interpolation techniques, a semi‐implicit stress interpolation correction is proposed to correct the cell‐face interpolation of the stress in the divergence operator of the momentum balance. Investigating the entry‐flow problem of the 4:1 contraction benchmark, we demonstrate that the numerical methods are robust over a wide range of Weissenberg numbers and significantly alleviate the high Weissenberg number problem. The accuracy of the results is evaluated in a detailed mesh convergence study.  相似文献   

9.
This paper develops and applies a linear viscoelastic model for bending and torsional modes of fluid membranes, based on the nonlinear Cosserat surface fluid model. The linearized fluid membrane model in spherical and cylindrical geometries is shown to decouple bending and torsional viscoelastic modes. It is found that solutions of the membrane viscoelastic model to small-amplitude oscillatory bending and torsion allows for the measurement of the bending and the torsion viscosity. The model and its potential in characterizing the bending and torsion viscoelasticity of membranes complements the on-going efforts to establish the role of curvature in dissipative process of biological membranes.  相似文献   

10.
Some equivalence conditions are formulated for non-linear models of polymer melts and solutions that are analogous to known conditions for three-constant linear rheological equations. The resulting model is analysed in simple shear and elongational flows. The kinematics of finite elastoviscous strains is considered in an appendix.  相似文献   

11.
An experimental method is developed to perform Hopkinson tests by means of viscoelastic bars by considering the wave propagation attenuation and dispersion due to the material rheological properties and the bar radial inertia (geometric effect). A propagation coefficient, representative of the wave dispersion and attenuation, is evaluated experimentally. Thus, the Pochhammer and Chree frequency equation is not necessary. Any bar cross-section shapes can be employed, and the knowledge of the bar mechanical properties is useless. The propagation coefficients for two PMMA bars with different diameters and for an elastic aluminum alloy bar are evaluated. These coefficients are used to determine the normal forces at the free end of a bar and at the ends of two bars held in contact. As an application, the mechanical impedance of an accelerometer is evaluated. A part of this work has been performed in the Laboratoire Matériaux Endommagement Fiabilité of the Ecole Nationale Supérieure des Arts et Métiers de Bordeaux.  相似文献   

12.
Two‐phase immiscible fluids in a two‐dimensional micro‐channels network are considered. The incompressible Stokes equations are used to describe the Newtonian fluid flow, while the Oldroyd‐B rheological model is used to capture the viscoelastic behavior. In order to perform numerical simulations in a complex geometry like a micro‐channels network, the volume penalization method is implemented. To follow the interface between the two fluids, the level‐set method is used, and the dynamics of the contact line is modeled by Cox law. Numerical results show the ability of the method to simulate two‐phase flows and to follow properly the contact line between the two immiscible fluids. Finally, simulations with realistic parameters are performed to show the difference when a Newtonian fluid is pushed by a viscoelastic fluid instead of a Newtonian one. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

13.
Two‐dimensional flows past a stationary circular cylinder near a plane boundary are numerically simulated using an immersed interface method with second‐order accuracy. Instead of a fixed wall, a moving wall with no‐slip boundary is considered to avoid the complex involvement of the boundary layer and to focus only on the shear‐free wall proximity effects for investigating the force dynamics and flow fields. To analyze the convergence and accuracy of our implementation, numerical studies have been first performed on a simple test problem of rotational flow, where the second order of convergence is confirmed through numerical experiments and an optimal range of relative grid‐match ratio of Lagrangian to Eulerian grid sizes has been recommended. By comparing the force quantities and the Strouhal number, the accuracy of this method has been demonstrated on the flow past a stationary isolated cylinder. The cylinder is then put in proximity to the wall to investigate the shear‐free wall proximity effects in the low Reynolds number regime (20≤Re≤200). The gap ratio, e/D, where e denotes the gap between the cylinder and the moving wall and D denotes the diameter of the cylinder, is taken from 0.10 to 2.00 to determine the critical gap ratio, (e/D)critical, for the alternate vortex shedding, where the fluid forces, flow fields and the streamwise velocity profiles are studied. One of the key findings is that the (e/D)critical for the alternate vortex shedding decreases as the Reynolds number increases. We also find that, in this low Reynolds number regime, the mean drag coefficient increases and peaks at e/D = 0.5 with the increase of e/D and keeps decreasing gently from e/D = 0.5 to e/D = 2.0, while the mean lift coefficient decreases monotonically with the increase of e/D. New correlations are then proposed for computing force coefficients as a function of Re and e/D for a cylinder in the vicinity of a moving plane wall. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

14.
The aim of this paper is to examine the Dufour and Soret effects on the two-dimensional magnetohydrodynamic (MHD) steady flow of an electrically conducting viscous fluid bounded by infinite sheets. An incompressible viscous fluid fills the porous space. The mathematical analysis is performed in the presence of viscous dissipation, Joule heating, and a first-order chemical reaction. With suitable transformations, the governing partial differential equations through momentum, energy, and concentration laws are transformed into ordinary differential equations. The resulting equations are solved by the homotopy analysis method (HAM). The convergence of the series solutions is ensured. The effects of the emerging parameters, the skin friction coefficient, the Nusselt number, and the Sherwood number are analyzed on the dimensionless velocities, temperature, and concentration fields.  相似文献   

15.
In this paper analytical solutions for the steady fully developed laminar fluid flow in the parallel-plate and cylindrical channels partially filled with a porous medium and partially with a clear fluid are presented. The Brinkman-extended Darcy equation is utilized to model the flow in a porous region. The solutions account for the boundary effects and for the stress jump boundary condition at the interface recently suggested by Ochoa-Tapia and Whitaker. The dependence of the velocity on the Darcy number and on the adjustable coefficient in the stress jump boundary condition is investigated. It is shown that accounting for a jump in the shear stress at the interface essentially influences velocity profiles.  相似文献   

16.
ABSTRACT

In this paper, effects of two wires magnetic field on heat transfer and biomagnetic fluid flow in an aneurysm have been investigated using the ferrohydrodynamics model. Using the finite volume method and the SIMPLE algorithm, the governing equations have been discretised. Simulations have been carried out for both conditions of wires in the same and opposite directions and different magnetic numbers of 41 and 82. Results show that the magnetic field causes a decrease in heat transfer of blood flow towards the walls. Moreover, major energy loss or pressure drop, arising from mean wall shear stress, decreases but local or minor energy loss, arising from aneurysm vortexes, increases. Furthermore, risk factors of aneurysm rupture is decreased under the effect of the magnetic field. The effective contact surface between drug-coated magnetic nanoparticles and the aneurysm tissue may increase and residence time of drug on the cells of the region would decrease.  相似文献   

17.
There have been a few recent numerical implementations of the stress‐jump condition at the interface of conjugate flows, which couple the governing equations for flows in the porous and homogenous fluid domains. These previous demonstration cases were for two‐dimensional, planar flows with simple geometries, for example, flow over a porous layer or flow through a porous plug. The present study implements the interfacial stress‐jump condition for a non‐planar flow with three velocity components, which is more realistic in terms of practical flow applications. The steady, laminar, Newtonian flow in a stirred micro‐bioreactor with a porous scaffold inside was investigated. It is shown how to implement the interfacial jump condition on the radial, axial, and swirling velocity components. To avoid a full three‐dimensional simulation, the flow is assumed to be independent of the azimuthal direction, which makes it an axisymmetric flow with a swirling velocity. The present interface treatment is suitable for non‐flat surfaces, which is achieved by applying the finite volume method based on body‐fitted and multi‐block grids. The numerical simulations show that a vortex breakdown bubble, attached to the free surface, occurs above a certain Reynolds number. The presence of the porous scaffold delays the onset of vortex breakdown and confines it to a region above the scaffold. Copyright © 2008 John Wiley & Sons, Ltd.  相似文献   

18.
This work describes numerical simulations of a hot jet in cross‐flow with applications to anti‐ice systems of aircraft engine nacelles. Numerical results are compared with experimental measurements obtained at ONERA to evaluate the performances of LES in this industrial context. The combination of complex geometries requiring unstructured meshes and high Reynolds number does not allow the resolution of boundary layers so that wall models must be employed. In this framework, the relative influence of subgrid‐scale modelling and conjugate heat transfer in LESs of aerothermal flows is evaluated. After a general overview of the transverse jet simulation results, a LES coupled with a heat transfer solver in the walls is used to show that thermal boundary conditions at the wall have more influence on the results than subgrid scale models. Coupling fluid flow and heat transfer in solids simulations is the only method to specify their respective thermal boundary conditions. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

19.
Current existing main nuclear thermal‐hydraulics (T‐H) system analysis codes, such as RALAP5, TRACE, and CATHARE, play a crucial role in the nuclear engineering field for the design and safety analysis of nuclear reactor systems. However, two‐fluid model used in these T‐H system analysis codes is ill posed, easily leading to numerical oscillations, and the classical first‐order methods for temporal and special discretization are widely employed for numerical simulations, yielding excessive numerical diffusion. Two‐fluid seven‐equation two‐pressure model is of particular interest due to the inherent well‐posed advantage. Moreover, high‐order accuracy schemes have also attracted great attention to overcome the challenge of serious numerical diffusion induced by low‐order time and space schemes for accurately simulating nuclear T‐H problems. In this paper, the semi‐implicit solution algorithm with high‐order accuracy in space and time is developed for this well‐posed two‐fluid model and the robustness and accuracy are verified and assessed against several important two‐phase flow benchmark tests in the nuclear engineering T‐H field, which include two linear advection problems, the oscillation problem of the liquid column, the Ransom water faucet problem, the reversed water faucet problem, and the two‐phase shock tube problem. The following conclusions are achieved. (1) The proposed semi‐implicit solution algorithm is robust in solving two‐phase flows, even for fast transients and discontinuous solutions. (2) High‐order schemes in both time and space could prevent excessive numerical diffusion effectively and the numerical simulation results are more accurate than those of first‐order time and space schemes, which demonstrates the advantage of using high‐order schemes.  相似文献   

20.
The modelling of mass transfer in porous media presenting capacitance effects by a dispersion-convection equation for the mobile fluid and a diffusion equation for the stagnant fluid has been shown (Piquemal, 1992) to be erroneous in the general case, because it is assumed that the mean concentration of the flowing fluid equals the point concentration at the boundary of the stagnant fluid. This boundary condition cannot be realized. This paper gives the conditions that allows the use of this kind of model with an acceptable approximation. The problem has been solved in the case of two schematic structures: the first is a cylindrical tube with stagnant pockets in its wall, the second a stratified medium. The characteristic lengths of the mobile and immobile domain must obey a criterion in which the porous medium characteristics and the flow velocity appear.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号