5,228 research outputs found
On acceleration of Krylov-subspace-based Newton and Arnoldi iterations for incompressible CFD: replacing time steppers and generation of initial guess
We propose two techniques aimed at improving the convergence rate of steady
state and eigenvalue solvers preconditioned by the inverse Stokes operator and
realized via time-stepping. First, we suggest a generalization of the Stokes
operator so that the resulting preconditioner operator depends on several
parameters and whose action preserves zero divergence and boundary conditions.
The parameters can be tuned for each problem to speed up the convergence of a
Krylov-subspace-based linear algebra solver. This operator can be inverted by
the Uzawa-like algorithm, and does not need a time-stepping. Second, we propose
to generate an initial guess of steady flow, leading eigenvalue and eigenvector
using orthogonal projection on a divergence-free basis satisfying all boundary
conditions. The approach, including the two proposed techniques, is illustrated
on the solution of the linear stability problem for laterally heated square and
cubic cavities
Viscous heating effects in fluids with temperature-dependent viscosity: triggering of secondary flows
Viscous heating can play an important role in the dynamics of fluids with
strongly temperature-dependent viscosities because of the coupling between the
energy and momentum equations. The heat generated by viscous friction produces
a local temperature increase near the tube walls with a consequent decrease of
the viscosity and a strong stratification in the viscosity profile. The problem
of viscous heating in fluids was investigated and reviewed by Costa & Macedonio
(2003) because of its important implications in the study of magma flows.
Because of the strong coupling between viscosity and temperature, the
temperature rise due to the viscous heating may trigger instabilities in the
velocity field, which cannot be predicted by a simple isothermal Newtonian
model. When viscous heating produces a pronounced peak in the temperature
profile near the walls, a triggering of instabilities and a transition to
secondary flows can occur because of the stratification in the viscosity
profile. In this paper we focus on the thermal and mechanical effects caused by
viscous heating. We will present the linear stability equations and we will
show, as in certain regimes, these effects can trigger and sustain a particular
class of secondary rotational flows which appear organised in coherent
structures similar to roller vortices. This phenomenon can play a very
important role in the dynamics of magma flows in conduits and lava flows in
channels and, to our knowledge, it is the first time that it has been
investigated by a direct numerical simulation.Comment: 18 pages manuscript, 10 figures, to be published in Journal of Fluid
Mechanics (2005
A pencil distributed finite difference code for strongly turbulent wall-bounded flows
We present a numerical scheme geared for high performance computation of
wall-bounded turbulent flows. The number of all-to-all communications is
decreased to only six instances by using a two-dimensional (pencil) domain
decomposition and utilizing the favourable scaling of the CFL time-step
constraint as compared to the diffusive time-step constraint. As the CFL
condition is more restrictive at high driving, implicit time integration of the
viscous terms in the wall-parallel directions is no longer required. This
avoids the communication of non-local information to a process for the
computation of implicit derivatives in these directions. We explain in detail
the numerical scheme used for the integration of the equations, and the
underlying parallelization. The code is shown to have very good strong and weak
scaling to at least 64K cores
Towards a solution of the closure problem for convective atmospheric boundary-layer turbulence
We consider the closure problem for turbulence in the dry convective atmospheric boundary
layer (CBL). Transport in the CBL is carried by small scale eddies near the surface and large
plumes in the well mixed middle part up to the inversion that separates the CBL from the
stably stratified air above. An analytically tractable model based on a multivariate Delta-PDF
approach is developed. It is an extension of the model of Gryanik and Hartmann [1] (GH02)
that additionally includes a term for background turbulence. Thus an exact solution is derived
and all higher order moments (HOMs) are explained by second order moments, correlation
coefficients and the skewness. The solution provides a proof of the extended universality
hypothesis of GH02 which is the refinement of the Millionshchikov hypothesis (quasi-
normality of FOM). This refined hypothesis states that CBL turbulence can be considered as
result of a linear interpolation between the Gaussian and the very skewed turbulence regimes.
Although the extended universality hypothesis was confirmed by results of field
measurements, LES and DNS simulations (see e.g. [2-4]), several questions remained
unexplained. These are now answered by the new model including the reasons of the
universality of the functional form of the HOMs, the significant scatter of the values of the
coefficients and the source of the magic of the linear interpolation. Finally, the closures
61
predicted by the model are tested against measurements and LES data. Some of the other
issues of CBL turbulence, e.g. familiar kurtosis-skewness relationships and relation of area
coverage parameters of plumes (so called filling factors) with HOM will be discussed also
Parallel numerical modeling of hybrid-dimensional compositional non-isothermal Darcy flows in fractured porous media
This paper introduces a new discrete fracture model accounting for
non-isothermal compositional multiphase Darcy flows and complex networks of
fractures with intersecting, immersed and non immersed fractures. The so called
hybrid-dimensional model using a 2D model in the fractures coupled with a 3D
model in the matrix is first derived rigorously starting from the
equi-dimensional matrix fracture model. Then, it is dis-cretized using a fully
implicit time integration combined with the Vertex Approximate Gradient (VAG)
finite volume scheme which is adapted to polyhedral meshes and anisotropic
heterogeneous media. The fully coupled systems are assembled and solved in
parallel using the Single Program Multiple Data (SPMD) paradigm with one layer
of ghost cells. This strategy allows for a local assembly of the discrete
systems. An efficient preconditioner is implemented to solve the linear systems
at each time step and each Newton type iteration of the simulation. The
numerical efficiency of our approach is assessed on different meshes, fracture
networks, and physical settings in terms of parallel scalability, nonlinear
convergence and linear convergence
- …