Turbulent Viscosity (turbulent + viscosity)

Distribution by Scientific Domains


Selected Abstracts


Coupled ghost fluid/two-phase level set method for curvilinear body-fitted grids

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 9 2007
Juntao Huang
Abstract A coupled ghost fluid/two-phase level set method to simulate air/water turbulent flow for complex geometries using curvilinear body-fitted grids is presented. The proposed method is intended to treat ship hydrodynamics problems. The original level set method for moving interface flows was based on Heaviside functions to smooth all fluid properties across the interface. We call this the Heaviside function method (HFM). The HFM requires fine grids across the interface. The ghost fluid method (GFM) has been designed to explicitly enforce the interfacial jump conditions, but the implementation of the jump conditions in curvilinear grids is intricate. To overcome these difficulties a coupled GFM/HFM method was developed in which approximate jump conditions are derived for piezometric pressure and velocity and pressure gradients based on exact continuous velocity and stress and jump in momentum conditions with the jump in density maintained but continuity of the molecular and turbulent viscosities imposed. The implementation of the ghost points is such that no duplication of memory storage is necessary. The level set method is adopted to locate the air/water interface, and a fast marching method was implemented in curvilinear grids to reinitialize the level set function. Validations are performed for three tests: super- and sub-critical flow without wave breaking and an impulsive plunging wave breaking over 2D submerged bumps, and the flow around surface combatant model DTMB 5512. Comparisons are made against experimental data, HFM and single-phase level set computations. The proposed method performed very well and shows great potential to treat complicated turbulent flows related to ship flows. Copyright © 2007 John Wiley & Sons, Ltd. [source]


Progress in the Modelling of Air Flow Patterns in Softwood Timber Kilns

ASIA-PACIFIC JOURNAL OF CHEMICAL ENGINEERING, Issue 3-4 2004
T.A.G. Langrish
Progress in modelling air flow patterns in timber kilns using computational fluid dynamics (CFD) is reviewed in this work. These simulations are intended to predict the distribution of the flow in the fillet spaces between boards in a hydraulic model of a timber kiln. Here, the flow regime between the boards is transitional between laminar and turbulent flow, with Reynolds numbers of the order of 5000. Running the simulation as a transient calculation has shown few problems with convergence issues, reaching a mass residual of 0.2% of the total inflow after 40 to 100 iterations per time step for time steps of 0.01 s. Grid sensitivity studies have shown that non-uniform grids are necessary because of the sudden changes in flow cross section, and the flow simulations are insensitive to grid refinement for non-uniform grids with more than 300,000 cells. The best agreement between the experimentally-measured flow distributions between fillet spaces and those predicted by the simulation have been achieved for (effective) bulk viscosities between the laminar viscosity for water and ten times that value. This change in viscosity is not very large (less than an order of magnitude), given that effective turbulent viscosities are typically several orders of magnitude greater than laminar ones. This result is consistent with the transitional flows here. The effect of weights above the stack can reduce the degree of non-uniformity in air velocities through the stack, especially when thick weights are used, because the stack may then be separated from the eddy at the top of the plenum chamber. [source]


On coupling the Reynolds-averaged Navier,Stokes equations with two-equation turbulence model equations

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 2 2006
Seungsoo Lee
Abstract Two methods for coupling the Reynolds-averaged Navier,Stokes equations with the q,, turbulence model equations on structured grid systems have been studied; namely a loosely coupled method and a strongly coupled method. The loosely coupled method first solves the Navier,Stokes equations with the turbulent viscosity fixed. In a subsequent step, the turbulence model equations are solved with all flow quantities fixed. On the other hand, the strongly coupled method solves the Reynolds-averaged Navier,Stokes equations and the turbulence model equations simultaneously. In this paper, numerical stabilities of both methods in conjunction with the approximated factorization-alternative direction implicit method are analysed. The effect of the turbulent kinetic energy terms in the governing equations on the convergence characteristics is also studied. The performance of the two methods is compared for several two- and three-dimensional problems. Copyright © 2005 John Wiley & Sons, Ltd. [source]


Hybrid LES-RANS modelling: a one-equation SGS model combined with a k,, model for predicting recirculating flows

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 9 2003
L. Davidson
Abstract A hybrid LES-RANS modelling approach is proposed. RANS is used in the near wall regions (y+,60), and the turbulence is modelled with a k,, model. LES is used in the remaining part of the flow, and the SGS turbulence is modelled with a one-equation ksgs model. The same continuity and momentum equations are solved throughout the domain, the only difference being that the turbulent viscosity is taken from the k,, model in the RANS region, and from the one-equation ksgs model in the LES region. The new modelling approach is applied to two incompressible flow test cases. They are fully developed flow in a plane channel and the flow over a 2D-hill in a channel. Copyright © 2003 John Wiley & Sons, Ltd. [source]


Gas dynamics of the central few parsec region of NGC 1068 fuelled by the evolving nuclear star cluster

MONTHLY NOTICES OF THE ROYAL ASTRONOMICAL SOCIETY, Issue 4 2010
M. Schartmann
ABSTRACT Recently, high-resolution observations with the help of the near-infrared adaptive optics integral field spectrograph Spectrograph for INtegral Field Observations in the Near Infrared (SINFONI) at the Very Large Telescope proved the existence of massive and young nuclear star clusters in the centres of a sample of Seyfert galaxies. With the help of three-dimensional high-resolution hydrodynamical simulations with the Pluto code, we follow the evolution of such clusters, especially focusing on stellar mass loss feeding gas into the ambient interstellar medium and driving turbulence. This leads to a vertically wide distributed clumpy or filamentary inflow of gas on large scales (tens of parsec), whereas a turbulent and very dense disc builds up on the parsec scale. In order to capture the relevant physics in the inner region, we treat this disc separately by viscously evolving the radial surface density distribution. This enables us to link the tens of parsec-scale region (accessible via SINFONI observations) to the (sub-)parsec-scale region (observable with the mid-infrared interferometer instrument and via water maser emission). Thereby, this procedure provides us with an ideal testbed for data comparison. In this work, we concentrate on the effects of a parametrized turbulent viscosity to generate angular momentum and mass transfer in the disc and additionally take star formation into account. Most of the input parameters are constrained by available observations of the nearby Seyfert 2 galaxy NGC 1068, and we discuss parameter studies for the free parameters. At the current age of its nuclear starburst of 250 Myr, our simulations yield disc sizes of the order of 0.8,0.9 pc, gas masses of 106 M, and mass transfer rates of 0.025 M, yr,1 through the inner rim of the disc. This shows that our large-scale torus model is able to approximately account for the disc size as inferred from interferometric observations in the mid-infrared and compares well to the extent and mass of a rotating disc structure as inferred from water maser observations. Several other observational constraints are discussed as well. [source]


On the orbital period modulation of RS CVn binary systems

MONTHLY NOTICES OF THE ROYAL ASTRONOMICAL SOCIETY, Issue 1 2005
Antonino F. Lanza
ABSTRACT The Applegate hypothesis proposed to explain the orbital period modulation of RS Canum Venaticorum (RS CVn) close binaries (Applegate 1992) is considered in the framework of a general model to treat the angular momentum exchanges within the convective envelope of a magnetically active star. This model assumes that the convection zone is strictly adiabatic and that the Taylor,Proudman balance holds, leading to an internal angular velocity constant over cylindrical surfaces co-axial with the rotation axis. It turns out that the angular velocity perturbations, whatever their origin, can be expressed in terms of the eigenfunctions of the equation of angular momentum conservation with stress-free boundary conditions. Moreover, a lower limit for the energy dissipation rate in a turbulent convection zone can be set, thanks to the extremal properties of the eigenfunctions. This approach allows to apply precise constraints on the amplitude and the radial profile of the angular velocity variations that are required to explain the observed orbital period changes in classical RS CVn binaries (i.e. with orbital period longer than 1,2 d and a subgiant secondary component). It is found that an angular velocity change as large as 10 per cent of the unperturbed angular velocity at the base of the stellar convection zone is needed. Such a large change is not compatible with the observations. Moreover, it would produce an energy dissipation rate much larger than the typical luminosities of the active components of RS CVn systems, except in the case that fast rotation and internal magnetic fields reduce the turbulent viscosity by at least 2 orders of magnitude with respect to the value given by the mixing-length theory. Therefore, the model proposed by Applegate should be rejected, at least in the case of classical RS CVn close binaries. Possible alternative models are briefly discussed, emphasizing the effects of intense magnetic fields (, 10 T) on the internal structure of magnetically active stars and the dynamics of close binary systems. [source]


The response of a turbulent accretion disc to an imposed epicyclic shearing motion

MONTHLY NOTICES OF THE ROYAL ASTRONOMICAL SOCIETY, Issue 1 2000
Ulf Torkelsson
We excite an epicyclic motion, the amplitude of which depends on the vertical position, z, in a simulation of a turbulent accretion disc. An epicyclic motion of this kind may be caused by a warping of the disc. By studying how the epicyclic motion decays, we can obtain information about the interaction between the warp and the disc turbulence. A high-amplitude epicyclic motion decays first by exciting inertial waves through a parametric instability, but its subsequent exponential damping may be reproduced by a turbulent viscosity. We estimate the effective viscosity parameter, ,v, pertaining to such a vertical shear. We also gain new information on the properties of the disc turbulence in general, and measure the usual viscosity parameter, ,h, pertaining to a horizontal (Keplerian) shear. We find that, as is often assumed in theoretical studies, ,v is approximately equal to ,h and both are much less than unity, for the field strengths achieved in our local box calculations of turbulence. In view of the smallness (,0.01) of ,v and ,h we conclude that for ,pgaspmag,10 the time-scale for diffusion or damping of a warp is much shorter than the usual viscous time-scale. Finally, we review the astrophysical implications. [source]


Effects of local thermodynamics and of stellar mass ratio on accretion disc stability in close binaries

ASTRONOMISCHE NACHRICHTEN, Issue 8 2009
G. Lanzafame
Abstract Inflow kinematics at the inner Lagrangian point L1, gas compressibility, and physical turbulent viscosity play a fundamental role on accretion disc dynamics and structure in a close binary (CB). Physical viscosity supports the accretion disc development inside the primary gravitational potential well, developing the gas radial transport, converting mechanical energy into heat. The Stellar-Mass-Ratio (SMR) between the compact primary and the secondary star (M1/M2) is also effective, not only in the location of the inner Lagrangian point, but also in the angular kinematics of the mass transfer and in the geometry ofthe gravitational potential wells. In this work we pay attention in particular to the role ofthe SMR, evaluating boundaries, separating theoretical domains in compressibility-viscosity graphs where physical conditions allow a well-bound disc development, as a function ofmass transfer kinematic conditions. In such domains, the lower is the gas compressibility (the higher the polytropic index ,), the higher is the physical viscosity (,) requested. In this work, we show how the boundaries of such domains vary as a function of M1/M2. Conclusions as far as dwarf novae outbursts are concerned, induced by mass transfer rate variations, are also reported. The smaller M1/M2, the shorter the duration of the active-to-quiet and vice-versa transitional phases. Time-scales are of the order of outburst duration of SU Uma, OY Car, Z Cha and SS Cyg-like objects. Moreover, conclusions as far as active-quiet-active phenomena in a CB, according to viscous-thermal instabilities, in accordance to such domains, are also reported (© 2009 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim) [source]


On the maximal value of the turbulent , -parameter in accretion discs

ASTRONOMISCHE NACHRICHTEN, Issue 7 2009
P. Abolmasov
Abstract In this short paper we show that making turbulence two-rather than three-dimensional may increase the effective turbulent viscosity by about 40 %. Dimensionless hydrodynamical viscosity parameters up to ,max = 0.25 Mt2 may be obtained in this approach, which are in better agreement with the observational data on non-stationary accretion than the values obtained in numerical simulations. However, the , -parameter values known from observations are still several times higher (© 2009 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim) [source]


Predicting the Displacement of Yoghurt by Water in a Pipe Using CFD

CHEMICAL ENGINEERING & TECHNOLOGY (CET), Issue 7 2007
M. Regner
Abstract A numerical scheme based on the volume of fluid (VOF) method for predicting the displacement of one liquid by another has been verified versus electrical resistance tomography (ERT) and ultrasonic velocity profile (UVP) measurements for the displacement of yoghurt by water. The scheme using the VOF method predicts the skewed phase distribution as measured using ERT and the global structure of the velocity profile as measured using UVP. The phase distribution using the VOF method was compared with the results using the species transport model which allows for mixing between the phases. The species transport model was found to be less suitable for predicting the displacement of yoghurt by water since the turbulence model was unable to accurately predict the turbulent viscosity in the mixing zone between yoghurt and water, which resulted in a too high rate of mixing. [source]