In diesem Beitrag wird die Wirksamkeit statischer Mischer in verschiedenen Anordnungen und Strömungskonfigurationen untersucht. Auf Grundlage umfangreicher numerischer Untersuchungen werden die Anwendungsgrenzen von spiralförmigen statischen Mischern zur Verbesserung des Wärmeübergangs in Kühlkanälen von Werkzeugmaschinen aufgezeigt. Die numerischen Simulationen wurden mit der kommerziellen Computational-Fluid-Dynamics (CFD)-Software, ANSYS Fluent 2020 R2, durchgeführt. Diese Studie zeigt, dass es einen optimalen Anwendungsbereich für statische Mischer als Wärmeaustauschverstärker in Abhängigkeit von der Strömungsgeschwindigkeit, dem übertragenen Wärmestrom und der Wärmeleitfähigkeit des Werkzeugs gibt. Die Untersuchungen in diesem Beitrag beschränken sich auf einphasige Strömungen in kreisförmigen Querschnitten und geraden Kanalgeometrien. Als repräsentatives Anwendungsbeispiel für eine Werkzeugmaschine wird die Kühlung eines einfachen Spritzgießwerkzeugs untersucht. Die durchgeführten Analysen zeigen, dass der Einsatz von statischer Mischelemente zur Verbesserung der Wärmeübertragung sehr effektiv ist, insbesondere bei Strömungen mit niedrigen bis mittleren Reynoldszahlen, konturnaher Kühlung, hohen Wärmestromwerten sowie hoher Wärmeleitfähigkeit des Werkzeugmaterials. / In this contribution, the effectiveness of helical static mixers in different arrangements and flow configurations/regimes is explored. By means of a thorough numerical analysis the application limits of helical static mixers for the heat transfer enhancement inside cooling channels of machine tools is provided. The numerical simulations were processed with the commercial finite volume Computational Fluid Dynamics (CFD) code, ANSYS Fluent 2020 R2. This study shows that there exists an optimal range of application for static mixers as heat exchange intensifier depending on the flow speed, the transmitted heat flow and the thermal conductivity of the tool. The investigations of this contribution are restricted to single-phase flow in circular cross-sections and straight channel geometries. As a representative application example for a machine tooling, the cooling of a simple injection mould is investigated. The research carried out reveals that the application of static mixing elements for enhancement of heat transfer is very effective, particularly for fluid flow with low to medium Reynolds numbers, close-contour cooling, high values of heat fluxes as well as high thermal conductivity of the tooling material.

One elementary part of a centrifugal compressor is the volute, which is located downstream the impeller. Its purpose is to collect the flow and increase the static pressure by converting kinetic energy into potential energy. Despite its significant effect onto the design point and operating range of the compressor, the number of publications regarding this component is quite small. Therefore, a numerical optimization of the volute housing is performed in order to identify important geometric parameters and find an optimal volute geometry. For this purpose, a new density-based CFD solver for all Mach numbers is developed as well as an automated geometry generation tool for the volute housing. The results show, that a volute with an inlet eccentricity of 0.9 and a slightly lower radial volute channel offers the best compressor efficiency. Moreover, the actual cross-sectional shape of the volute has only a minor influence onto the performance. As a result, the isentropic efficiency could be improved by up to 2 % compared to the reference compressor model, in particular at high off-design flow rates. These results are a novelty in the scientific community and help to design more efficient compressors. / Das Spiralgehäuse eines Radialverdichters wird im Gegensatz zum Laufrad kaum in wissenschaftlichen Arbeiten untersucht. Um wichtige Geometrieparameter und Einflussfaktoren dieses Bauteils zu identifizieren, wird daher eine Optimierung mittels genetischer Algorithmen durchgeführt. Dazu wird zunächst ein dichte-basierter CFD-Löser entwickelt und validiert, um die komplexe Strömung in einem Radialverdichter mit hoher Genauigkeit simulieren zu können. Darauf aufbauend wird das Spiralgehäuse parametrisiert und ein Programm entwickelt, welches die komplexe Geometrie automatisiert erstellt. Durch die neuartige Kombination von numerischer Optimierung, automatisierter Geometrieerstellung und CFD-Simulation des Spiralgehäuses können erstmals Aussagen zur optimalen Geometrie sowie über Verlusteffekte für eine Vielzahl an Geomtrievarianten getroffen werden. Mit Hilfe dieses Wissens können sparsamere und effizientere Radialkompressoren für viele Bereiche des Maschinenbaus entwickelt werden.

One of the most exciting predictions of Einstein's theory of gravitation that have not yet been proven experimentally by a direct detection are gravitational waves. These are tiny distortions of the spacetime itself, and a world-wide effort to directly measure them for the first time with a network of large-scale laser interferometers is currently ongoing and expected to provide positive results within this decade. One potential source of measurable gravitational waves is the inspiral and merger of two compact objects, such as binary black holes. Successfully finding their signature in the noise-dominated data of the detectors crucially relies on accurate predictions of what we are looking for. In this thesis, we present a detailed study of how the most complete waveform templates can be constructed by combining the results from (A) analytical expansions within the post-Newtonian framework and (B) numerical simulations of the full relativistic dynamics. We analyze various strategies to construct complete hybrid waveforms that consist of a post-Newtonian inspiral part matched to numerical-relativity data. We elaborate on exsisting approaches for nonspinning systems by extending the accessible parameter space and introducing an alternative scheme based in the Fourier domain. Our methods can now be readily applied to multiple spherical-harmonic modes and precessing systems. In addition to that, we analyze in detail the accuracy of hybrid waveforms with the goal to quantify how numerous sources of error in the approximation techniques affect the application of such templates in real gravitational-wave searches. This is of major importance for the future construction of improved models, but also for the correct interpretation of gravitational-wave observations that are made utilizing any complete waveform family. In particular, we comprehensively discuss how long the numerical-relativity contribution to the signal has to be in order to make the resulting hybrids accurate enough, and for currently feasible simulation lengths we assess the physics one can potentially do with template-based searches. / Eine der aufregendsten Vorhersagen aus Einsteins Gravitationstheorie, die bisher noch nicht direkt durch ein Experiment nachgewiesen werden konnten, sind Gravitationswellen. Dies sind winzige Verzerrungen der Raumzeit selbst, und es wird erwartet, dass das aktuelle Netzwerk von groß angelegten Laserinterferometern im kommenden Jahrzehnt die erste direkte Gravitationswellenmessung realisieren kann. Eine potentielle Quelle von messbaren Gravitationswellen ist das Einspiralen und Verschmelzen zweier kompakter Objekte, wie z.B. ein Binärsystem von Schwarzen Löchern. Die erfolgreiche Identifizierung ihrer charakteristischen Signatur im Rausch-dominierten Datenstrom der Detektoren hängt allerdings entscheidend von genauen Vorhersagen ab, was wir eigentlich suchen. In dieser Arbeit wird detailliert untersucht, wie die komplettesten Wellenformenmodelle konstruiert werden können, indem die Ergebnisse von (A) analytischen Entwicklungen im post-Newtonschen Verfahren und (B) numerische Simulationen der voll-relativistischen Bewegungen verknüpft werden. Es werden verschiedene Verfahren zur Erstellung solcher "hybriden Wellenformen", bei denen der post-Newtonsche Teil mit numerischen Daten vervollständigt wird, analysiert. Existierende Strategien für nicht-rotierende Systeme werden vertieft und der beschriebene Parameterraum erweitert. Des Weiteren wird eine Alternative im Fourierraum eingeführt. Die entwickelten Methoden können nun auf multiple sphärisch-harmonische Moden und präzedierende Systeme angewandt werden. Zusätzlich wird die Genauigkeit der hybriden Wellenformen mit dem Ziel analysiert, den Einfluss verschiedener Fehlerquellen in den Näherungstechniken zu quantifizieren und die resultierenden Einschränkungen bei realen Anwendungen abzuschätzen. Dies ist von größter Bedeutung für die zukünftige Entwicklung von verbesserten Modellen, aber auch für die korrekte Interpretation von Gravitationswellenbeobachtungen, die auf Grundlage solcher Familien von Wellenformen gemacht worden sind. Insbesondere wird diskutiert, wie lang der numerische Anteil des Signals sein muss, um die Hybride genau genug konstruieren zu können. Für die aktuell umsetzbaren Simulationslängen wird die Physik eingeschätzt, die mit Hilfe von Modell-basierten Suchen potentiell untersucht werden kann.

Es werden nahfeldoptische Phänomene an Raster-Sonden-Geometrien untersucht. Hierzu zählen die Feldverstärkung an laserbestrahlten Rastersondenspitzen und die Feldverteilung in nahfeldoptischen Apertur- und Koaxialspitzen. Zur Berechnung der Feldverteilung werden verschiedene numerische Verfahren wie die Methode der Randelemente (BEM) und die Methode der Finiten Integration (FIT) angewendet. Die durchgeführten Berechnungen sagen eine erhebliche Feldüberhöhung im Nahfeld von laserbestrahlten Rastersondenspitzen voraus. Es wird diskutiert, ob diese Feldüberhöhung ursächlich für die Modifikation der Oberfläche unter der Sondenspitze verantwortlich ist. Zur weiteren Klärung der Strukturierungsmechanismen wird die Stromantwort des Tunnelübergangs auf Laserbestrahlung experimentell untersucht. Die Detektion des Stromsignals erfolgt mit einem neuen Vorverstärker, der die Eigenschaften einer hohen Eingangsimpedanz, einer hohen Transimpedanz und einer hohen Bandbreite miteinander vereinigt. Weiterhin wird die Feldverteilung in SNOM-Spitzen, insbesondere in den sogenannten Koaxialspitzen, bestimmt. Es wird untersucht wie sich koaxial Moden in koaxialen SNOM-Spitzen anregen lassen und wie sie sich ausbreiten. Ausgehend von bereits existierenden Entwürfen solcher Spitzen wird diskutiert, welche Strukturen sich zur Anregung von koaxialen Moden eignen.

Die Entwicklung neuer Materialien und neuer Werkstoffe wird heute international als Schlüsseltechnologie mit Querschnittscharakter und Schrittmacherfunktion für viele industrielle Bereiche eingestuft. Die Wirtschaftskraft der hoch entwickelten Industriegesellschaften hängt zunehmend von Erfolgen in der Materialwissenschaft und der Werkstofftechnologie ab. Die Forschungsaktivitäten in der Profillinie 1 sind gekennzeichnet durch Interdisziplinarität und Vernetzung von Forschungsvorhaben. Von besonderer Bedeutung ist darüber hinaus die zusätzliche Verzahnung mit der Profillinie 6 der TU Chemnitz “Modellierung, Simulation, Hochleistungsrechnen“, um die Material- und Werkstoffforschung durch den intelligenten Einsatz leistungsstarker Rechentechnik weniger kostenintensiv gestalten zu können.

Im Projekt „Weiterentwicklung und Optimierung von Be- und Entladesystemen für Tank- und Erdbeckenspeicher“ des Forschungsprogramms „Solarthermie 2000plus“ wurden von 2004 bis 2008 unterschiedliche Konstruktionen zur Be- und Entladung thermischer Speicher experimentell und mit Hilfe von Simulationsrechnungen untersucht. Wesentliche Ziele dieses Projektes waren: - die Systematisierung und Analyse bestehender Konstruktionen, - die Weiterentwicklung und Optimierung von Be- und Entladesystemen, - die Ableitung von Richtlinien für Planung und Auslegung, - die Untersuchung der Übertragbarkeit vom Experiment auf Mittel- und Großspeicher sowie - die Einbindung von Industriepartnern (Kooperation und Praxistests). Das Projekt wurde an der Technischen Universität Chemnitz in Zusammenarbeit mit dem Projektpartner Technische Universität Ilmenau durchgeführt. Die Ergebnisse der Untersuchungen und die gewonnenen Erkenntnisse werden im Bericht detailliert erläutert.

No description available.

Tagungsband des Ohde-Kolloquiums 2014. Die Fachtagung fand am 26.03.2014 an der TU Dresden statt.

This thesis reports recent numerical investigation of the rheological behavior of fresh concrete using the Distinct Element Method (DEM). Some relevant questions of the concrete rheology e.g. the influence of the concrete composition on the rheological behavior of the fresh concrete, the experimental determination of the Bingham rheological constants as well as the use of these constants in the numerical simulation were discussed thoroughly. An important topic of the performed investigation was the development of the numerical model for fresh concrete which enables simple, fast and stable predictive simulation of different technological operations with fresh concrete. Firstly, in a literature survey, the state-of-the-art of the numerical simulation of fresh concrete was presented and critically discussed in order to show advantages and disadvantages of other methods and modeling approaches. Open (unsolved) questions were highlighted and the basis for their investigation is created within this thesis. Fundamental concepts of the rheology were then presented and the basic rheological models of viscoelastic materials were considered; the rheological behaviors of different types of concretes were presented and its influencing factors were discussed. Additionally main methods for scientific investigation and testing of the fresh concrete were shown. The test methods were critically discussed in order to select the test, which has been used as a reference experimental test for the numerical simulations. Chosen reference experimental test was the slump flow test. The slump flow test was thoroughly analyzed and an analytical solution was developed which helps to interpret the results of measurements and provides a link between rheological constants and measured quantities. In a further step an extensive experimental program was carried out in order to investigate the rheological behavior of fresh concrete and get the input data for numerical simulation. Firstly, the experiments on macrolevel were performed. Here the rheological behavior of the fresh concrete flow in different tests was investigated (slump and slump flow tests, L-Box). Further, the experiments on mesolevel with polymer on Carbopol basis and mortar were developed and performed in order to investigate the interaction between distinct particles suspended in a fluid matrix. The necessary material parameters, especially those representative of the fluid suspension micromechanical behavior, i.e. the force-displacement relationship, yield force and bond strength, were determined by these experiments. The slump flow test was used as the basic test to calibrate the model for fresh concrete (key data: slump value, slump flow diameter (for concretes with a soft consistency) and the time of spreading). Thus, the decisive phenomena of the fresh concrete flow were highlighted, control points for a contact model were selected and the initial input data for the development of the contact model was obtained. Next, the user-defined contact model was developed and implemented into the Particle Flow Code ITASCA. The contact model was completely described and its limitations discussed. Then, the set of numerical tools was developed, which enable simplified and stable numerical simulation of the fresh concrete with particular behavior, i.e. automatic generation of the concrete with given particle grading, amount of fibers and air, automatic recalculation of the micromechanical parameters of the contact model from given initial yield stress and plastic viscosity. The model was calibrated by slump flow test simulations and validated by corresponding analytical approach. Further, the role of different model parameters was investigated by simulating the slump flow test. Furthermore, for verification of the model several additional experiments were simulated, i.e. L-Box and LCPC-box test. The results of modeling were compared with experimental results and discussed in detail. All numerical simulations provide qualitatively as well as quantitatively correct results and hence adequately represent the phenomena observed in real experiments. The thesis closes with general conclusions and outlook of the work. In the future, the developed contact model and tools of the “Virtual concrete laboratory” could be modified in order to extend the potential of the laboratory to cover such properties as thixotropic behavior of fresh concrete or simulating hardening of the concrete and behavior of the hardened concrete.

Die Modellierung von komplexen Systemen, wie dem Untergrund, ist ein Hilfsmittel zur Beschreibung der in der Realität ablaufenden Prozesse. Die Durchführung von Experimenten an einem Modell, um qualitative Aussagen über ein reales System zu erhalten, wird als Simulation bezeichnet. Dabei können vielfältige Modelle, wie z.B. physikalische und mathematische, zum Einsatz kommen. Die ungesättigte Bodenzone (vadose Zone) bezeichnet den Bereich zwischen der Landoberfläche und dem Grundwasserspiegel, innerhalb dessen der Wassergehalt geringer als bei Vollsättigung, und der Druck geringer als der Atmosphärendruck ist. Dieser Bodenbereich hat für die Landwirtschaft, Geobiologie, aerobe Abbauprozesse und Grundwasserneubildung eine große Bedeutung. Für die Nachbildung von Strömungs- und Stofftransportprozessen der ungesättigten Bodenzone existieren numerische Simulationsprogramme. Ziel der vorliegenden Arbeit ist eine umfangreiche Validierung des Programms PCSiWaPro® (entwickelt an der TU-Dresden, Institut für Abfallwirtschaft und Altlasten) für unterschiedliche Anwendungsfälle. Ein weiteres Ziel der Arbeit besteht in der Untersuchung der Anwendbarkeit des aktuellen Stands des Simulationsprogramms PCSiWaPro® auf unterschiedliche Praxisfälle bei Kombination der ungesättigten Bodenzone mit technischen Anlagen. Vier Anwendungsfälle mit unterschiedlichen Zielen wurden dafür im Rahmen dieser Arbeit untersucht: die Simulation von dezentraler Abwasserversickerung (Kleinkläranlage - KKA) anhand entsprechender Säulen- und Feldversuche, die Berechnung der Grundwasserneubildung am Beispiel von Lysimetern, der Wasserhaushalt von Erddämmen und die Modellierung von Deponieabdeckungssystemen. Die Anwendungsfälle unterscheiden sich durch den Zweck der Simulation, die Geometrie, die Größe, die festgelegten Anfangs- und Randbedingungen, die Simulationszeit, die Materialien, das Koordinatensystem sowie die Ein- und Ausgabewerte. Die Simulationsergebnisse konnten eindeutig zeigen, dass das Programm PCSiWaPro® für alle im Rahmen der vorliegenden Arbeit untersuchten Fälle, mit unterschiedlichen Strömungsregimen, Stofftransport-Parametern, Randbedingungen, Koordinatensystemen sowie Raum- und Zeitdiskretisierungen anwendbar ist. Die Simulationsergebnisse der Säulenversuche am Beispiel dezentraler Abwasserversickerung zeigten eine sehr gute Übereinstimmung zwischen gemessenen und mittels PCSiWa-Pro® berechneten Werten des Wasser- und Stoffhaushaltes (Druckhöhe, Abfluss und Stoff-konzentration) der untersuchten Bodentypen B3 (schwachschluffiger Sand), B4 (Grobsand) und B5 (mittelschluffiger Sand). Die Wurzel des mittleren quadratischen Fehlers (RMSE) betrug für die Berechnung der Druckhöhe 1,84 cm bei B5, 3,61 cm bei B3 und 1,27 cm bei B4. Die relative Abweichung betrug für die Berechnung der Druckhöhe 2,19 % bei B5, 1,3 % bei B3 und ca. 5,3 % bei B4. Die Durchführung der Sensitivitätsanalyse der für die Modellierung relevanten Parameter zeigte eine sehr hohe Sensitivität der VAN GENUCHTEN-Parameter und der gesättigten hydraulischen Leitfähigkeit des Bodens. Darüber hinaus führten die Parameter nach DIN 4220 und die mithilfe von Pedotransferfunktionen aus Siebanalysen genommenen Parameter zu unterschiedlichen Ergebnissen. Im Rahmen des am Institut für Abfallwirtschaft und Altlasten durchgeführten Projektes EGSIM wurden die Programme SENSIT und ISSOP (in Zusammenarbeit mit DUALIS GmbH IT Solution) entwickelt und zur Parameteridentifikation/-kalibrierung benutzt. Die im Rahmen dieser Arbeit erzielten Ergebnisse konnten nachweisen, unter welchen Bedingungen eine Nachklärung des vollbiologisch gereinigten Abwassers innerhalb der Bodenzone möglich ist, so dass am Ort der Beurteilung (Grundwasseroberfläche) kein unzulässiger Schadstoffeintrag erfolgt. In Bezug auf die KKA-Feldmodelle ist die Anwendung des rotationssymmetrischen Koordinatensystems als Voraussetzung der Realität besser zu entsprechen und nicht als Option zu betrachten. Darüber hinaus wurden anhand der Feldmodelle verschiedene Szenarien mit kontinuierlicher und diskontinuierlicher Versickerung sowie zwei unterschiedlich großen Einleitflächen durchgeführt. Das Programm PCSiWaPro® ist sowohl für ungesättigte als auch für variabel-gesättigte porösen Medien anwendbar. Dies wurde im Rahmen der Simulation des Wasserhaushaltes eines Erddamms nachgewiesen. Die durchschnittliche relative Abweichung zwischen gemessenen und mittels PCSiWaPro® berechneten Wasserständen des entsprechenden Beobachtungspunkts im untersuchten Dammkörper lag bei 0,08 % (entspricht 5,8 cm bei einer Müchtigkeit von ca. 70 m) und das Bestimmtheitsmaß (R2) betrug 0,987. Die Simulation des Wasserhaushaltes unterschiedlicher Deponieabdichtungssystemen mittels PCSiWaPro® zeigte im Allgemeinen ein funktionierendes Ableiten des auf Deponien anfallenden Regenwassers (auch bei Starkregenereignissen). Darüber hinaus haben die durchgeführten Bewuchs-Modelle nachweisen können, dass die Vegetation der Deponieoberflächen den Wassergehalt, durch Pflanzenwurzelentzug bzw. Evapotranspiration, reduzieren können. Die Simulationsergebnisse der durchgeführten Szenarien des Wasserhaushaltes von Lysimetern zur prognostischen Berechnung der Grundwasserneubildung mittels PCSiWaPro® konnten nachweisen, dass das Programm für die Berechnung der Grundwasserneubildungsrate für diesen Zweck anwendbar ist. Die relativen Abweichungen der be-rechneten von den gemessenen Grundwasserneubildungsraten sind auf die verwendeten Materialparameter sowie auf Vernachlässigung der möglicherweise in Lysimetern sich befin-denden Makroporen (duale Porosität) zurückzuführen. / The modelling of complex systems such as the underground is a means to describe the processes occurring in the reality. The conducting of experiments on a model to obtain qualitative evidence about a real system is referred to as a simulation. Thereby, various models (e.g. physical and mathematical models) can be used. The unsaturated zone (vadose zone) is the region between the land surface and the water table, in which the water content is less than full saturation, and the pressure is lower than the atmospheric pressure. The unsaturated zone is very significant for agriculture, geobiology, aerobic degradation processes and groundwater recharge. The processes of water flow and solute transport in the unsaturated zone can be described by means of numerical simulation programs. The aim of the present work is a comprehensive validation of the simulation program PCSiWaPro® (developed at the TU-Dresden, Institute of Waste Management and Contaminated Site Treatment) for different applications. Another aim of this work is to investigate the applicability of the current version of PCSiWaPro® for different cases of a combination between the unsaturated zone and technical facilities. Four application cases with different objectives were investigated within the present work, which are: the simulation of decentralized wastewater infiltration with corresponding column and field experiments, the computation of groundwater recharge by means of lysimeters, the water balance of earth dams and the modelling of landfill covering systems. The application cases differ from each other by the objective of the simulation, the geometry, the size, the specified initial and boundary conditions, the simulation time, the applied materials, the coordinate system, the input and output data. The simulation results clearly showed that PCSiWaPro® is applicable for all investigated cases under consideration of different flow and solute transport regimes, parameters, boundary conditions, spatial and temporal discretization, and coordinate systems. The simulation results of the experimental soil columns for the decentralized treated wastewater infiltration case showed a very good agreement between measured and computed values of water and solute balance (pressure head, flow and solute concentration) of the investigated soil types B3 (slightly silty sand), B4 (coarse sand / gravel) and B5 (medium silty sand). The root of the mean squared error (RMSE) for the computation of the pressure head was 1,84 cm at B5, 3,61 cm at B3 and 1,27 cm at B4. The relative deviation in case of pressure head computation was 2,19 % at B5, 1,3 % at B3 and 5,3 % at B4. The implementation of the sensitivity analysis of the relevant parameters for the modelling showed a very high sensitivity of the VAN GENUCHTEN parameters and the saturated hydraulic conductivity of the soil. Moreover, the parameters according to DIN 4220 led to different results than the estimated ones according to pedotransfer methods based on sieve analysis. Within the project EGSIM, which was carried out at the Institute for waste management and contaminated sites treatment in collaboration with DUALIS GmbH IT Solution, the programs SENSIT and ISSOP were developed and used for parameter identification/ calibration. The results obtained in this Work showed under which conditions is a secondary treatment of full biologically treated wastewater in the soil possible, so that no unallowable pollutants entry in the groundwater occurs. With regard to the field models of this application the implementation of the rotationally symmetric coordinate system should be considered as a condition and not as an option for a better corresponding to the reality. Furthermore, different scenarios of the field models were carried out with continuous and discontinuous infiltration, as well as under different initiation areas. PCSiWaPro® could be applied for both unsaturated and variably-saturated porous media. This could be proven by the simulation of the water balance in an earth dam. The average relative deviation between measured and simulated water levels of the corresponding observation point in the investigated dam embankment was 0,08 % (corresponding to 5,8 cm at 70 m thickness) and the coefficient of determination (R2) was 0,987. In general, the simulation of the water balance using PCSiWaPro® of different landfill covering systems showed a successful draining of the falling rainwater (even under heavy rainfall). In addition, the implemented vegetation models have proven that the vegetation of the landfill surface can reduce the water content in the landfill by evapotranspiration and water uptake by roots. The water balance simulation results of the scenarios for the computation of groundwater recharge by means of lysimeters showed that the program is applicable for this case. The relative deviation of the simulated from the measured groundwater recharge rates occur due to the implemented material parameters as well as to the neglect of macro pores effects (dual porosity).