KLUEDO RSS FeedKLUEDO Dokumente/documents
https://kluedo.ub.uni-kl.de/index/index/
Mon, 05 Dec 2011 14:09:24 +0200Mon, 05 Dec 2011 14:09:24 +0200A one-dimensional model of the pressing section of a paper machine including dynamic capillary effects
https://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/2322
This work presents the dynamic capillary pressure model (Hassanizadeh, Gray, 1990, 1993a) adapted for the needs of paper manufacturing process simulations. The dynamic capillary pressure-saturation relation is included in a one-dimensional simulation model for the pressing section of a paper machine. The one-dimensional model is derived from a two-dimensional model by averaging with respect to the vertical direction. Then, the model is discretized by the finite volume method and solved by Newton’s method. The numerical experiments are carried out for parameters typical for the paper layer. The dynamic capillary pressure-saturation relation shows significant influence on the distribution of water pressure. The behaviour of the solution agrees with laboratory experiments (Beck, 1983).O. Iliev; G. Printsypar; S. Riefreporthttps://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/2322Thu, 12 May 2011 14:09:24 +0200Finite volume discretization of equations describing nonlinear diffusion in Li-Ion batteries
https://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/2225
Numerical modeling of electrochemical process in Li-Ion battery is an emerging topic of great practical interest. In this work we present a Finite Volume discretization of electrochemical diffusive processes occurring during the operation of Li-Ion batteries. The system of equations is a nonlinear, time-dependent diffusive system, coupling the Li concentration and the electric potential. The system is formulated at length-scale at which two different types of domains are distinguished, one for the electrolyte and one for the active solid particles in the electrode. The domains can be of highly irregular shape, with electrolyte occupying the pore space of a porous electrode. The material parameters in each domain differ by several orders of magnitude and can be non-linear functions of Li ions concentration and/or the electrical potential. Moreover, special interface conditions are imposed at the boundary separating the electrolyte from the active solid particles. The field variables are discontinuous across such an interface and the coupling is highly non- linear, rendering direct iteration methods ineffective for such problems. We formulate a Newton iteration for an purely implicit Finite Volume discretization of the coupled system. A series of numerical examples are presented for different type of electrolyte/electrode configurations and material parameters. The convergence of the Newton method is characterized both as function of nonlinear material parameters as well as the nonlinearity in the interface conditions.P. Popov; Y. Vutov; S. Margenov; O. Ilievreporthttps://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/2225Tue, 10 Aug 2010 09:23:20 +0200Modeling of species and charge transport in Li-Ion Batteries based on non-equilibrium thermodynamics
https://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/2219
In order to improve the design of Li ion batteries the complex interplay of various physical phenomena in the active particles of the electrodes and in the electrolyte has to be balanced. The separate transport phenomena in the electrolyte and in the active particle as well as their coupling due to the electrochemical reactions at the interfaces between the electrode particles and the electrolyte will inuence the performance and the lifetime of a battery. Any modeling of the complex phenomena during the usage of a battery has therefore to be based on sound physical and chemical principles in order to allow for reliable predictions for the response of the battery to changing load conditions. We will present a modeling approach for the transport processes in the electrolyte and the electrodesbased on non-equilibrium thermodynamics and transport theory. The assumption of local charge neutrality, which is known to be valid in concentrated electrolytes, is explicitly used to identify the independent thermodynamic variables and uxes. The theory guarantees strictly positive entropy production. Dierences to other theories will be discussed.A. Latz; J. Zausch; O. Ilievreporthttps://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/2219Wed, 21 Jul 2010 14:31:57 +0200On a numerical subgrid upscaling algorithm for Stokes-Brinkman equations
https://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/2218
This paper discusses a numerical subgrid resolution approach for solving the Stokes-Brinkman system of equations, which is describing coupled ow in plain and in highly porous media. Various scientic and industrial problems are described by this system, and often the geometry and/or the permeability vary on several scales. A particular target is the process of oil ltration. In many complicated lters, the lter medium or the lter element geometry are too ne to be resolved by a feasible computational grid. The subgrid approach presented in the paper is aimed at describing how these ne details are accounted for by solving auxiliary problems in appropriately chosen grid cells on a relatively coarse computational grid. This is done via a systematic and a careful procedure of modifying and updating the coecients of the Stokes-Brinkman system in chosen cells. This numerical subgrid approach is motivated from one side from homogenization theory, from which we borrow the formulations for the so called cell problem, and from the other side from the numerical upscaling approaches, such as Multiscale Finite Volume, Multiscale Finite Element, etc. Results on the algorithm's eciency, both in terms of computational time and memory usage, are presented. Comparison with solutions on full ne grid (when possible) are presented in order to evaluate the accuracy. Advantages and limitations of the considered subgrid approach are discussed.O. Iliev; Z. Lakdawala; V. Starikoviciusreporthttps://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/2218Wed, 21 Jul 2010 14:30:10 +0200Variational multiscale Finite Element Method for flows in highly porous media
https://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/2210
We present a two-scale finite element method for solving Brinkman’s and Darcy’s equations. These systems of equations model fluid flows in highly porous and porous media, respectively. The method uses a recently proposed discontinuous Galerkin FEM for Stokes’ equations byWang and Ye and the concept of subgrid approximation developed by Arbogast for Darcy’s equations. In order to reduce the “resonance error” and to ensure convergence to the global fine solution the algorithm is put in the framework of alternating Schwarz iterations using subdomains around the coarse-grid boundaries. The discussed algorithms are implemented using the Deal.II finite element library and are tested on a number of model problems.O. Iliev; R. Lazarov; J. Willemsreporthttps://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/2210Wed, 21 Jul 2010 14:08:14 +0200Design of pleated filters by computer simulations
https://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/2083
Four aspects are important in the design of hydraulic lters. We distinguish between two cost factors and two performance factors. Regarding performance, filter eciencynd lter capacity are of interest. Regarding cost, there are production considerations such as spatial restrictions, material cost and the cost of manufacturing the lter. The second type of cost is the operation cost, namely the pressure drop. Albeit simulations should and will ultimately deal with all 4 aspects, for the moment our work is focused on cost. The PleatGeo Module generates three-dimensional computer models of a single pleat of a hydraulic lter interactively. PleatDict computes the pressure drop that will result for the particular design by direct numerical simulation. The evaluation of a new pleat design takes only a few hours on a standard PC compared to days or weeks used for manufacturing and testing a new prototype of a hydraulic lter. The design parameters are the shape of the pleat, the permeabilities of one or several layers of lter media and the geometry of a supporting netting structure that is used to keep the out ow area open. Besides the underlying structure generation and CFD technology, we present some trends regarding the dependence of pressure drop on design parameters that can serve as guide lines for the design of hydraulic lters. Compared to earlier two-dimensional models, the three-dimensional models can include a support structure.A. Wiegmann; L. Cheng; E. Glatt; O. Iliev; S. Riefreporthttps://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/2083Wed, 29 Apr 2009 12:21:31 +0200A Graph-Laplacian approach for calculating the effective thermal conductivity of complicated fiber geometries
https://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/2004
Abstract. An efficient approach to the numerical upscaling of thermal conductivities of fibrous media, e.g. insulation materials, is considered. First, standard cell problems for a second order elliptic equation are formulated for a proper piece of random fibrous structure, following homogenization theory. Next, a graph formed by the fibers is considered, and a second order elliptic equation with suitable boundary conditions is solved on this graph only. Replacing the boundary value problem for the full cell with an auxiliary problem with special boundary conditions on a connected subdomain of highly conductive material is justified in a previous work of the authors. A discretization on the graph is presented here, and error estimates are provided. The efficient implementation of the algorithm is discussed. A number of numerical experiments is presented in order to illustrate the performance of the proposed method.O. Iliev; R. Lazarov; J. Willemsreporthttps://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/2004Wed, 25 Jun 2008 13:01:53 +0200An efficient approach for upscaling properties of composite materials with high contrast of coefficients
https://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1987
An efficient approach for calculating the effective heat conductivity for a class of industrial composite materials, such as metal foams, fibrous glass materials, and the like, is discussed. These materials, used in insulation or in advanced heat exchangers, are characterized by a low volume fraction of the highly conductive material (glass or metal) having a complex, network-like structure and by a large volume fraction of the insulator (air). We assume that the composite materials have constant macroscopic thermal conductivity tensors, which in principle can be obtained by standard up-scaling techniques, that use the concept of representative elementary volumes (REV), i.e. the effective heat conductivities of composite media can be computed by post-processing the solutions of some special cell problems for REVs. We propose, theoretically justify, and numerically study an efficient approach for calculating the effective conductivity for media for which the ratio of low and high conductivities satisfies 1. In this case one essentially only needs to solve the heat equation in the region occupied by the highly conductive media. For a class of problems we show, that under certain conditions on the microscale geometry, the proposed approach produces an upscaled conductivity that is O() close to the exact upscaled permeability. A number of numerical experiments are presented in order to illustrate the accuracy and the limitations of the proposed method. Applicability of the presented approach to upscaling other similar problems, e.g. flow in fractured porous media, is also discussed.R. Ewing; O. Iliev; R. Lazarov; I. Rybak; J. Willemsreporthttps://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1987Wed, 18 Jun 2008 15:28:21 +0200An analysis of one regularization approach for solution of pure Neumann problem
https://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1992
In this paper, the analysis of one approach for the regularization of pure Neumann problems for second order elliptical equations, e.g., Poisson’s equation and linear elasticity equations, is presented. The main topic under consideration is the behavior of the condition number of the regularized problem. A general framework for the analysis is presented. This allows to determine a form of regularization term which leads to the “natural” asymptotic of the condition number of the regularized problem with respect to mesh parameter. Some numerical results, which support theoretical analysis are presented as well. The main motivation for the presented research is to develop theoretical background for an efficient and robust implementation of the solver for pure Neumann problems for the linear elasticity equations. Such solvers usually are needed in a number of domain decomposition methods, e.g. FETI. Developed approaches are planed to be used in software, developing in ITWM, e.g. KneeMech simulation software.E. Savenkov; H. Andrä; O. Ilievreporthttps://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1992Wed, 18 Jun 2008 15:27:25 +0200On two-level preconditioners for flow in porous media
https://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1975
Two-level domain decomposition preconditioner for 3D flows in anisotropic highly heterogeneous porous media is presented. Accurate finite volume discretization based on multipoint flux approximation (MPFA) for 3D pressure equation is employed to account for the jump discontinuities of full permeability tensors. DD/MG type preconditioner for above mentioned problem is developed. Coarse scale operator is obtained from a homogenization type procedure. The influence of the overlapping as well as the influence of the smoother and cell problem formulation is studied. Results from numerical experiments are presented and discussed.R. Ewing; O. Iliev; R. Lazarov; I. Rybakreporthttps://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1975Wed, 28 May 2008 10:24:09 +0200On upscaling heat conductivity for a class of industrial problems
https://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1974
Calculating effective heat conductivity for a class of industrial problems is discussed. The considered composite materials are glass and metal foams, fibrous materials, and the like, used in isolation or in advanced heat exchangers. These materials are characterized by a very complex internal structure, by low volume fraction of the higher conductive material (glass or metal), and by a large volume fraction of the air. The homogenization theory (when applicable), allows to calculate the effective heat conductivity of composite media by postprocessing the solution of special cell problems for representative elementary volumes (REV). Different formulations of such cell problems are considered and compared here. Furthermore, the size of the REV is studied numerically for some typical materials. Fast algorithms for solving the cell problems for this class of problems, are presented and discussed.O. Iliev; I. Rybak; J. Willemsreporthttps://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1974Wed, 28 May 2008 10:24:01 +0200On approximation property of multipoint flux approximation method
https://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1973
Approximation property of multipoint flux approximation (MPFA) approach for elliptic equations with discontinuous full tensor coefficients is discussed here. Finite volume discretization of the above problem is presented in the case of jump discontinuities for the permeability tensor. First order approximation for the fluxes is proved. Results from numerical experiments are presented and discussed.O. Iliev; I. Rybakreporthttps://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1973Wed, 28 May 2008 10:23:55 +0200On numerical upscaling for flows in heterogeneous porous media
https://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1972
A numerical upscaling approach, NU, for solving multiscale elliptic problems is discussed. The main components of this NU are: i) local solve of auxil- iary problems in grid blocks and formal upscaling of the obtained re sults to build a coarse scale equation; ii) global solve of the upscaled coarse scale equation; and iii) reconstruction of a fine scale solution by solving local block problems on a dual coarse grid. By its structure NU is similar to other methods for solving multiscale elliptic problems, such as the multiscale finite element method, the multiscale mixed finite element method, the numerical subgrid upscaling method, heterogeneous multiscale method, and the multiscale finite volume method. The difference with those methods is in the way the coarse scale equation is build and solved, and in the way the fine scale solution is reconstructed. Essential components of the presented here NU approach are the formal homogenization in the coarse blocks and the usage of so called multipoint flux approximation method, MPFA. Unlike the usual usage as MPFA as a discretiza- tion method for single scale elliptic problems with tensor discontinuous coefficients, we consider its usage as a part of a numerical upscaling approach. The main aim of this paper is to compare NU with the MsFEM. In particular, it is shown that the resonance effect, which limits the application of the Multiscale FEM, does not appear, or it is significantly relaxed, when the presented here numerical upscaling approach is applied.O. Iliev; I. Rybakreporthttps://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1972Wed, 28 May 2008 10:23:47 +0200On parallel numerical algorithms for simulating industrial filtration problems
https://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1968
The performance of oil filters used in the automotive industry can be significantly improved, especially when computer simulation is an essential component of the design process. In this paper, we consider parallel numerical algorithms for solving mathematical models describing the process of filtration, filtering out solid particles from liquid oil. The Navier-Stokes-Brinkmann system of equations is used to describe the laminar flow of incompressible isothermal oil. The space discretization in the complicated filter geometry is based on the finite-volume method. Special care is taken for an accurate approximation of velocity and pressure on the interface between the fluid and the porous media. The time discretization used here is a proper modification of the fractional time step discretization (cf. Chorin scheme) of the Navier-Stokes equations, where the Brinkmann term is considered at both, prediction and correction substeps. A data decomposition method is used to develop a parallel algorithm, where the domain is distributed among processors by using a structured reference grid. The MPI library is used to implement the data communication part of the algorithm. A theoretical model is proposed for the estimation of the complexity of the given parallel algorithm and a scalability analysis is done on the basis of this model. Results of computational experiments are presented, and the accuracy and efficiency of the parallel algorithm is tested on real industrial geometries.R. Ciegis; O. Iliev; Z. Lakdawalareporthttps://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1968Wed, 28 May 2008 10:23:18 +0200Numerical study of two-grid preconditioners for 1d elliptic problems with highly oscillating discontinuous coefficients
https://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1964
Abstract — Various advanced two-level iterative methods are studied numerically and compared with each other in conjunction with finite volume discretizations of symmetric 1-D elliptic problems with highly oscillatory discontinuous coefficients. Some of the methods considered rely on the homogenization approach for deriving the coarse grid operator. This approach is considered here as an alternative to the well-known Galerkin approach for deriving coarse grid operators. Different intergrid transfer operators are studied, primary consideration being given to the use of the so-called problemdependent prolongation. The two-grid methods considered are used as both solvers and preconditioners for the Conjugate Gradient method. The recent approaches, such as the hybrid domain decomposition method introduced by Vassilevski and the globallocal iterative procedure proposed by Durlofsky et al. are also discussed. A two-level method converging in one iteration in the case where the right-hand side is only a function of the coarse variable is introduced and discussed. Such a fast convergence for problems with discontinuous coefficients arbitrarily varying on the fine scale is achieved by a problem-dependent selection of the coarse grid combined with problem-dependent prolongation on a dual grid. The results of the numerical experiments are presented to illustrate the performance of the studied approaches.O. Iliev; R. Lazarov; J. Willemsreporthttps://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1964Wed, 28 May 2008 10:22:58 +0200Numerical Algorithems for Solving Problems of Multiphase Flows in Porous Media
https://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1744
R. Ciegis; O. Iliev; V. Starikovicius; K. Steinerreporthttps://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1744Wed, 31 May 2006 14:03:22 +0200On 3D Numerical Simulations of Viscoelastic Fluids
https://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1743
D. Niedziela; O. Iliev; A. Latzreporthttps://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1743Wed, 31 May 2006 14:03:04 +0200On Efficent Simulation of Non-Newtonian Flow in Saturated Porous Media with a Multigrid Adaptive Refinement Solver
https://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1683
Flow of non-Newtonian fluid in saturated porous media can be described by the continuity equation and the generalized Darcy law. Efficient solution of the resulting second order nonlinear elliptic equation is discussed here. The equation is discretized by a finite volume method on a cell-centered grid. Local adaptive refinement of the grid is introduced in order to reduce the number of unknowns. A special implementation approach is used, which allows us to perform unstructured local refinement in conjunction with the finite volume discretization. Two residual based error indicators are exploited in the adaptive refinement criterion. Second order accurate discretization of the fluxes on the interfaces between refined and non-refined subdomains, as well as on the boundaries with Dirichlet boundary condition, are presented here, as an essential part of the accurate and efficient algorithm. A nonlinear full approximation storage multigrid algorithm is developed especially for the above described composite (coarse plus locally refined) grid approach. In particular, second order approximation of the fluxes around interfaces is a result of a quadratic approximation of slave nodes in the multigrid - adaptive refinement (MG-AR) algorithm. Results from numerical solution of various academic and practice-induced problems are presented and the performance of the solver is discussed.W. Dörfler; O. Iliev; D. Stoyanov; D. Vassilevareporthttps://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1683Wed, 23 Nov 2005 14:29:15 +0100On convergence of certain finite difference discretizations for 1D poroelasticity interface problems
https://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1599
Finite difference discretizations of 1D poroelasticity equations with discontinuous coefficients are analyzed. A recently suggested FD discretization of poroelasticity equations with constant coefficients on staggered grid, [5], is used as a basis. A careful treatment of the interfaces leads to harmonic averaging of the discontinuous coefficients. Here, convergence for the pressure and for the displacement is proven in certain norms for the scheme with harmonic averaging (HA). Order of convergence 1.5 is proven for arbitrary located interface, and second order convergence is proven for the case when the interface coincides with a grid node. Furthermore, following the ideas from [3], modified HA discretization are suggested for particular cases. The velocity and the stress are approximated with second order on the interface in this case. It is shown that for wide class of problems, the modified discretization provides better accuracy. Second order convergence for modified scheme is proven for the case when the interface coincides with a displacement grid node. Numerical experiments are presented in order to illustrate our considerations.R. Ewing; O. Iliev; R. Lazarov; A. Naumovichreporthttps://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1599Mon, 31 Jan 2005 16:32:30 +0100On numerical solution of 1-D poroelasticity equations in a multilayered domain
https://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1579
In soil mechanics assumption of only vertical subsidence is often invoked and this leads to the one-dimensional model of poroelasticity. The classical model of linear poroelasticity is obtained by Biot [1], detailed derivation can be found e.g., in [2]. This model is applicable also to modelling certain processes in geomechanics, hydrogeology, petroleum engineering (see, e.g., [3, 8], in biomechanics (e.g., [9, 10]), in filtration (e.g., filter cake formation, see [15, 16, 17]), in paper manufacturing (e.g., [11, 12]), in printing (e.g., [13]), etc. Finite element and finite difference methods were applied by many authors for numerical solution of the Biot system of PDEs, see e.g. [3, 4, 5] and references therein. However, as it is wellknown, the standard FEM and FDM methods are subject to numerical instabilities at the first time steps. To avoid this, discretization on staggered grid was suggested in [4, 5]. A single layer deformable porous medium was considered there. This paper can be viewed as extension of [4, 5] to the case of multilayered deformable porous media. A finite volume discretization to the interface problem for the classical one-dimensional Biot model of consolidation process is applied here. Following assumptions are supposed to be valid: each of the porous layers is composed of incompressible solid matrix, it is homogeneous and isotropic. Furthermore, one of two following assumptions is valid: porous medium is not completely saturated and ﬂuid is incompressible or porous medium is completely saturated and fluid is slightly compressible. The reminder of the paper is organised as follows. Next section presents the mathematical model. Third section is devoted to the dicsretization of the continuous problem. Fourth section contains the results from the numerical experiments.F. Gaspar; O. Iliev; F. Lisbona; A. Naumovich; P. Vabishchevichreporthttps://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1579Wed, 27 Oct 2004 17:17:21 +0200Fluid structure interaction problems in deformable porous media: Toward permeability of deformable porous media
https://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1578
In this work the problem of fluid flow in deformable porous media is studied. First, the stationary fluid-structure interaction (FSI) problem is formulated in terms of incompressible Newtonian fluid and a linearized elastic solid. The flow is assumed to be characterized by very low Reynolds number and is described by the Stokes equations. The strains in the solid are small allowing for the solid to be described by the Lame equations, but no restrictions are applied on the magnitude of the displacements leading to strongly coupled, nonlinear fluid-structure problem. The FSI problem is then solved numerically by an iterative procedure which solves sequentially fluid and solid subproblems. Each of the two subproblems is discretized by finite elements and the fluid-structure coupling is reduced to an interface boundary condition. Several numerical examples are presented and the results from the numerical computations are used to perform permeability computations for different geometries.O. Iliev; A. Mikelic; P. Popovreporthttps://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1578Wed, 27 Oct 2004 17:16:34 +0200On Modelling and Simulation of Different Regimes for Liquid Polymer Moulding
https://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1556
In this paper we consider numerical algorithms for solving a system of nonlinear PDEs arising in modeling of liquid polymer injection. We investigate the particular case when a porous preform is located within the mould, so that the liquid polymer flows through a porous medium during the filling stage. The nonlinearity of the governing system of PDEs is due to the non-Newtonian behavior of the polymer, as well as due to the moving free boundary. The latter is related to the penetration front and a Stefan type problem is formulated to account for it. A finite-volume method is used to approximate the given differential problem. Results of numerical experiments are presented. We also solve an inverse problem and present algorithms for the determination of the absolute preform permeability coefficient in the case when the velocity of the penetration front is known from measurements. In both cases (direct and inverse problems) we emphasize on the specifics related to the non-Newtonian behavior of the polymer. For completeness, we discuss also the Newtonian case. Results of some experimental measurements are presented and discussed.R. Ciegis; O. Iliev; S. Rief; K. Steinerreporthttps://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1556Mon, 26 Jul 2004 11:04:53 +0200On the Performance of Certain Iterative Solvers for Coupled Systems Arising in Discretization of Non-Newtonian Flow Equations
https://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1555
Iterative solution of large scale systems arising after discretization and linearization of the unsteady non-Newtonian Navier–Stokes equations is studied. cross WLF model is used to account for the non-Newtonian behavior of the fluid. Finite volume method is used to discretize the governing system of PDEs. Viscosity is treated explicitely (e.g., it is taken from the previous time step), while other terms are treated implicitly. Different preconditioners (block–diagonal, block–triangular, relaxed incomplete LU factorization, etc.) are used in conjunction with advanced iterative methods, namely, BiCGStab, CGS, GMRES. The action of the preconditioner in fact requires inverting different blocks. For this purpose, in addition to preconditioned BiCGStab, CGS, GMRES, we use also algebraic multigrid method (AMG). The performance of the iterative solvers is studied with respect to the number of unknowns, characteristic velocity in the basic flow, time step, deviation from Newtonian behavior, etc. Results from numerical experiments are presented and discussed.O. Iliev; J. Linn; M. Moog; D. Niedziela; V. Starikoviciusreporthttps://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1555Mon, 26 Jul 2004 11:03:55 +0200Multigrid – adaptive local refinement solver for incompressible flows
https://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1511
A non-linear multigrid solver for incompressible Navier-Stokes equations, exploiting finite volume discretization of the equations, is extended by adaptive local refinement. The multigrid is the outer iterative cycle, while the SIMPLE algorithm is used as a smoothing procedure. Error indicators are used to define the refinement subdomain. A special implementation approach is used, which allows to perform unstructured local refinement in conjunction with the finite volume discretization. The multigrid - adaptive local refinement algorithm is tested on 2D Poisson equation and further is applied to a lid-driven flows in a cavity (2D and 3D case), comparing the results with bench-mark data. The software design principles of the solver are also discussed.O. Iliev; D. Stoyanovreporthttps://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1511Tue, 10 Feb 2004 10:12:45 +0100On a Multigrid Adaptive Refinement Solver for Saturated Non-Newtonian Flow in Porous Media
https://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1509
On a Multigrid Adaptive Refinement Solver for Saturated Non-Newtonian Flow in Porous Media A multigrid adaptive refinement algorithm for non-Newtonian flow in porous media is presented. The saturated flow of a non-Newtonian fluid is described by the continuity equation and the generalized Darcy law. The resulting second order nonlinear elliptic equation is discretized by a finite volume method on a cell-centered grid. A nonlinear full-multigrid, full-approximation-storage algorithm is implemented. As a smoother, a single grid solver based on Picard linearization and Gauss-Seidel relaxation is used. Further, a local refinement multigrid algorithm on a composite grid is developed. A residual based error indicator is used in the adaptive refinement criterion. A special implementation approach is used, which allows us to perform unstructured local refinement in conjunction with the finite volume discretization. Several results from numerical experiments are presented in order to examine the performance of the solver.W. Dörfler; O. Iliev; D. Stoyanov; D. Vassilevareporthttps://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1509Tue, 10 Feb 2004 10:06:38 +0100