Transient non-Darcy forced convective heat transfer from a flat plate embedded in a fluid-saturated porous medium

Transient non-Darcy forced convection on a flat plate embedded in a porous medium is investigated using the Forchheimer-extended Darcy law. A sudden uniform pressure gradient is applied along the flat plate, and at the same time, its wall temperature is suddenly raised to a high temperature. Both the momentum and energy equations are solved by retaining the unsteady terms. An exact velocity solution is obtained and substituted into the energy equation, which then is solved by means of a quasi-similarity transformation. The temperature field can be divided into the one-dimensional transient (downstream) region and the quasi-steady-state (upstream) region. Thus the transient local heat transfer coefficient can be described by connecting the quasi-steady-state solution and the one-dimensional transient solution. The non-Darcy porous inertia works to decrease the velocity level and the time required for reaching the steady-state velocity level. The porous-medium inertia delays covering of the plate by the steady-state thermal boundary layer. © 1990.

Estudo experimental do escoamento sobre placas planas munidas de protuberâncias quadradas e onduladas

Pós-graduação em Engenharia Mecânica - FEIS

Análise computacional dos campos de velocidade e temperatura do ambiente interno da usina termelétrica Santana - Amapá

A simulação numérica do escoamento de ar em ambientes internos é na atualidade o método mais apropriado para análise de conforto térmico em ambientes internos. O escoamento de ar nesses ambientes configura-se como um escoamento complexo, pois, em regra geral, é uma combinação de escoamentos cisalhantes livres (jatos) e cisalhantes de parede, além disso, esses escoamentos são governados por forças de inércia e forças de empuxo, caracterizando-o como de convecção mista. A combinação desses mecanismos cria um escoamento com características complexas, como zonas de recirculação, vórtices, descolamento e recolamento de camada-limite dentre outras. Portanto, a precisão da solução estará diretamente ligada, principalmente, na habilidade do modelo de turbulência adotado de reproduzir as características turbulentas do escoamento de ar e da transferência térmica. O objetivo principal do presente trabalho foi a simulação computacional do ambiente térmico interno do galpão que abriga os geradores e motores Wärtzilä da Usina Termelétrica Santana no estado do Amapá. A formulação matemática baseada na solução das equações gerais de conservação inclui uma análise dos principais modelos de turbulência aplicados ao escoamento de ar em ambientes internos, assim como os processos de transferência de calor associados. Na modelagem numérica o método de volumes finitos é usado na discretização das equações de conservação, através do código comercial Fluent-Airpak, que foi usado nas simulações computacionais para a análise dos campos de velocidade e temperatura do ar. A utilização correta do programa computacional foi testada e validada para o problema através da simulação precisa de casos retirados da literatura. Os resultados numéricos foram comparados a dados obtidos de medições experimentais realizados no galpão e apresentou boa concordância, considerando a complexidade do problema simulado, o objetivo da simulação em face da diminuição da temperatura no interior do galpão e, também, em função das limitações encontradas quando da tomada das medições experimentais. Além disso, foram feitas simulações de estratégias de melhoria do ambiente térmico da Usina, baseadas na realidade levantada e nos resultados da simulação numérica. Finalmente, foram realizadas simulações do protótipo de solução proposto para a diminuição da temperatura interna do galpão o que possibilitará um aumento, na faixa de 20 a 30%, do tempo de permanência no interior do galpão.

Influence of Goertler vortices spanwise wavelength on heat transfer rates

The boundary layer over concave surfaces can be unstable due to centrifugal forces, giving rise to Goertler vortices. These vortices create two regions in the spanwise direction—the upwash and downwash regions. The downwash region is responsible for compressing the boundary layer toward the wall, increasing the heat transfer rate. The upwash region does the opposite. In the nonlinear development of the Goertler vortices, it can be observed that the upwash region becomes narrow and the spanwise–average heat transfer rate is higher than that for a Blasius boundary layer. This paper analyzes the influence of the spanwise wavelength of the Goertler the heat transfer. The equation is written in vorticity-velocity formulation. The time integration is done via a classical fourth-order Runge-Kutta method. The spatial derivatives are calculated using high-order compact finite difference and spectral methods. Three different wavelengths are analyzed. The results show that steady Goertler flow can increase the heat transfer rates to values close to the values of turbulence, without the existence of a secondary instability. The geometry (and computation domain) are presented

Simulazione numerica di fenomeni aeroelastici con applicazione alle strutture civili

Il presente lavoro tratta lo studio dei fenomeni aeroelastici di interazione fra fluido e struttura, con il fine di provare a simularli mediante l’ausilio di un codice agli elementi finiti. Nel primo capitolo sono fornite alcune nozioni di fluidodinamica, in modo da rendere chiari i passaggi teorici fondamentali che portano alle equazioni di Navier-Stokes governanti il moto dei fluidi viscosi. Inoltre è illustrato il fenomeno della formazione di vortici a valle dei corpi tozzi dovuto alla separazione dello strato limite laminare, con descrizione anche di alcuni risultati ottenuti dalle simulazioni numeriche. Nel secondo capitolo vengono presi in rassegna i principali fenomeni di interazione fra fluido e struttura, cercando di metterne in luce le fondamenta della trattazione analitica e le ipotesi sotto le quali tale trattazione è valida. Chiaramente si tratta solo di una panoramica che non entra in merito degli sviluppi della ricerca più recente ma fornisce le basi per affrontare i vari problemi di instabilità strutturale dovuti a un particolare fenomeno di interazione con il vento. Il terzo capitolo contiene una trattazione più approfondita del fenomeno di instabilità per flutter. Tra tutti i fenomeni di instabilità aeroelastica delle strutture il flutter risulta il più temibile, soprattutto per i ponti di grande luce. Per questo si è ritenuto opportuno dedicargli un capitolo, in modo da illustrare i vari procedimenti con cui si riesce a determinare analiticamente la velocità critica di flutter di un impalcato da ponte, a partire dalle funzioni sperimentali denominate derivate di flutter. Al termine del capitolo è illustrato il procedimento con cui si ricavano sperimentalmente le derivate di flutter di un impalcato da ponte. Nel quarto capitolo è presentato l’esempio di studio dell’impalcato del ponte Tsing Ma ad Hong Kong. Sono riportati i risultati analitici dei calcoli della velocità di flutter e di divergenza torsionale dell’impalcato e i risultati delle simulazioni numeriche effettuate per stimare i coefficienti aerodinamici statici e il comportamento dinamico della struttura soggetta all’azione del vento. Considerazioni e commenti sui risultati ottenuti e sui metodi di modellazione numerica adottati completano l’elaborato.

Diagnostic techniques for EHD and MHD interaction

The impact of plasma technologies is growing both in the academic and in the industrial fields. Nowadays, a great interest is focused in plasma applications in aeronautics and astronautics domains. Plasma actuators based on the Magneto-Hydro-Dynamic (MHD) and Electro- Hydro-Dynamic (EHD) interactions are potentially able to suitably modify the fluid-dynamics characteristics around a flying body without utilizing moving parts. This could lead to the control of an aircraft with negligible response time, more reliability and improvements of the performance. In order to study the aforementioned interactions, a series of experiments and a wide number of diagnostic techniques have been utilized. The EHD interaction, realized by means of a Dielectric Barrier Discharge (DBD) actuator, and its impact on the boundary layer have been evaluated by means of two different experiments. In the first one a three phase multi-electrode flat panel actuator is used. Different external flow velocities (from 1 to 20m/s) and different values of the supplied voltage and frequency have been considered. Moreover a change of the phase sequence has been done to verify the influence of the electric field existing between successive phases. Measurements of the induced speed had shown the effect of the supply voltage and the frequency, and the phase order in the momentum transfer phenomenon. Gains in velocity, inside the boundary layer, of about 5m/s have been obtained. Spectroscopic measurements allowed to determine the rotational and the vibrational temperature of the plasma which lie in the range of 320 ÷ 440°K and of 3000 ÷ 3900°K respectively. A deviation from thermodynamic equilibrium had been found. The second EHD experiment is realized on a single electrode pair DBD actuator driven by nano-pulses superimposed to a DC or an AC bias. This new supply system separates the plasma formation mechanism from the acceleration action on the fluid, leading to an higher degree of the control of the process. Both the voltage and the frequency of the nano-pulses and the amplitude and the waveform of the bias have been varied during the experiment. Plasma jets and vortex behavior had been observed by means of fast Schlieren imaging. This allowed a deeper understanding of the EHD interaction process. A velocity increase in the boundary layer of about 2m/s had been measured. Thrust measurements have been performed by means of a scales and compared with experimental data reported in the literature. For similar voltage amplitudes thrust larger than those of the literature, had been observed. Surface charge measurements led to realize a modified DBD actuator able to obtain similar performances when compared with that of other experiments. However in this case a DC bias replacing the AC bias had been used. MHD interaction experiments had been carried out in a hypersonic wind tunnel in argon with a flow of Mach 6. Before the MHD experiments a thermal, fluid-dynamic and plasma characterization of the hypersonic argon plasma flow have been done. The electron temperature and the electron number density had been determined by means of emission spectroscopy and microwave absorption measurements. A deviation from thermodynamic equilibrium had been observed. The electron number density showed to be frozen at the stagnation region condition in the expansion through the nozzle. MHD experiments have been performed using two axial symmetric test bodies. Similar magnetic configurations were used. Permanent magnets inserted into the test body allowed to generate inside the plasma azimuthal currents around the conical shape of the body. These Faraday currents are responsible of the MHD body force which acts against the flow. The MHD interaction process has been observed by means of fast imaging, pressure and electrical measurements. Images showed bright rings due to the Faraday currents heating and exciting the plasma particles. Pressure measurements showed increases of the pressure in the regions where the MHD interaction is large. The pressure is 10 to 15% larger than when the MHD interaction process is silent. Finally by means of electrostatic probes mounted flush on the test body lateral surface Hall fields of about 500V/m had been measured. These results have been used for the validation of a numerical MHD code.

Stability, viscous dissipation and local thermal non-equilibrium in fluid saturated porous media

In this thesis, the field of study related to the stability analysis of fluid saturated porous media is investigated. In particular the contribution of the viscous heating to the onset of convective instability in the flow through ducts is analysed. In order to evaluate the contribution of the viscous dissipation, different geometries, different models describing the balance equations and different boundary conditions are used. Moreover, the local thermal non-equilibrium model is used to study the evolution of the temperature differences between the fluid and the solid matrix in a thermal boundary layer problem. On studying the onset of instability, different techniques for eigenvalue problems has been used. Analytical solutions, asymptotic analyses and numerical solutions by means of original and commercial codes are carried out.

Influences of deep convective cloud systems on tropospheric trace gases and photochemistry over the tropical West Pacific

A numerical model for studying the influences of deep convective cloud systems on photochemistry was developed based on a non-hydrostatic meteorological model and chemistry from a global chemistry transport model. The transport of trace gases, the scavenging of soluble trace gases, and the influences of lightning produced nitrogen oxides (NOx=NO+NO2) on the local ozone-related photochemistry were investigated in a multi-day case study for an oceanic region located in the tropical western Pacific. Model runs considering influences of large scale flows, previously neglected in multi-day cloud resolving and single column model studies of tracer transport, yielded that the influence of the mesoscale subsidence (between clouds) on trace gas transport was considerably overestimated in these studies. The simulated vertical transport and scavenging of highly soluble tracers were found to depend on the initial profiles, reconciling contrasting results from two previous studies. Influences of the modeled uptake of trace gases by hydrometeors in the liquid and the ice phase were studied in some detail for a small number of atmospheric trace gases and novel aspects concerning the role of the retention coefficient (i.e. the fraction of a dissolved trace gas that is retained in the ice phase upon freezing) on the vertical transport of highly soluble gases were illuminated. Including lightning NOx production inside a 500 km 2-D model domain was found to be important for the NOx budget and caused small to moderate changes in the domain averaged ozone concentrations. A number of sensitivity studies yielded that the fraction of lightning associated NOx which was lost through photochemical reactions in the vicinity of the lightning source was considerable, but strongly depended on assumptions about the magnitude and the altitude of the lightning NOx source. In contrast to a suggestion from an earlier study, it was argued that the near zero upper tropospheric ozone mixing ratios which were observed close to the study region were most probably not caused by the formation of NO associated with lightning. Instead, it was argued in agreement with suggestions from other studies that the deep convective transport of ozone-poor air masses from the relatively unpolluted marine boundary layer, which have most likely been advected horizontally over relatively large distances (both before and after encountering deep convection) probably played a role. In particular, it was suggested that the ozone profiles observed during CEPEX (Central Equatorial Pacific Experiment) were strongly influenced by the deep convection and the larger scale flow which are associated with the intra-seasonal oscillation.

Untersuchungen zur troposphärischen Photochemie am GAW-Observatorium Izana

Im Rahmen dieser Promotionsarbeit wurden Fragestellungen der troposphärischen Ozonproduktion, des photostationären Gleichgewicht des NO2-NO-O3-ROx-Systems, des Einflusses von heterogenen Reaktionen an Mineralstauboberflächen und der Bedeutung lokaler anthropogener Emissionen für Spurengasmessungen an einem Gebirgsobservatorium untersucht. Dazu wurde insbesondere ein atmosphärenchemisch-meteorologischer Datensatz ausgewertet, welcher während der Feldmesskampagne MINATROC II im Juli/August 2002 am Global Atmosphere Watch (GAW) Observatorium Izana (2360m über NN) auf Teneriffa gewonnen wurde. Der Datensatz enthält Messdaten von Aerosolen, Spurengasen, Strahlungsgrößen und meteorologischen Parametern. Eine dichte Sahara-Staubwolke passierte die Messstation zwischen dem 28. und 31.07.2002. Die Mischungsverhältnisse von O3, H2O2, NOx und ROx zeigten während dieses Ereignisses Abnahmen zwischen 10% und 50%. Insbesondere infolge der Reduktion von NO und HO2 wurde eine deutliche Verringerung der Netto-Ozon-Produktion in der Staubwolke berechnet. Die Bewertung der verschiedenen Terme der Ozonproduktion bzw. –vernichtung zeigte relativ geringe Beiträge der einzelnen photochemischen Verlustreaktionen am Ozonbudget. Im Rahmen einer Boxmodellstudie, welche die Eingangsdaten von MINATROC II verwendete, wurde ein Vergleich zwischen Simulation und Messung kurzlebiger Spurengase angestellt. Die Modellrechnungen weisen darauf hin, dass die reduzierten H2O2-Mischungsverhältnisse in der Sahara-Staubwolke unter Berücksichtigung heterogener Reaktionen von HO2 bzw. H2O2 an Stauboberflächen erklärt werden können. Die Betrachtungen zum photostationären Gleichgewicht des NO2-NO-O3-ROx-Systems ergeben hauptsächlich stark überhöhte Leighton-Verhältnisse. Die Abweichungen der Messdaten vom photostationären Gleichgewichtszustand sind dabei zum Verschmutzungsgrad der Atmosphäre antikorreliert. Eine Reihe potentieller Einflüsse auf das NO2-NO-O3-ROx-System wurde diskutiert, wobei gezeigt werden konnte, dass das bisherige Betrachtungsmodell nicht ausreicht, um die beobachteten Abweichungen vom photostationären Gleichgewicht zu erklären. Während für die Nachtmessungen von MINATROC II Bedingungen der freien Troposphäre vorlagen, waren die Messungen am Tage von der Kopplung zwischen photochemischer Aktivität und thermisch getriebener Hangströmung betroffen. Die zum Messort gelangte anthropogen beeinflusste Grenzschichtluft konnte dabei auf maximalen Mischungsanteilen zwischen 15% und 85% quantifiziert werden. Abschließend wurden verschiedene Aspekte der Quellen und der Verteilung anthropogener Verschmutzungen auf Teneriffa im Zusammenhang der zuvor behandelten Fragen diskutiert.

Measurement and modeling of cloud condensation nuclei in continental air

Atmospheric aerosol particles serving as cloud condensation nuclei (CCN) are key elements of the hydrological cycle and climate. Knowledge of the spatial and temporal distribution of CCN in the atmosphere is essential to understand and describe the effects of aerosols in meteorological models. In this study, CCN properties were measured in polluted and pristine air of different continental regions, and the results were parameterized for efficient prediction of CCN concentrations.The continuous-flow CCN counter used for size-resolved measurements of CCN efficiency spectra (activation curves) was calibrated with ammonium sulfate and sodium chloride aerosols for a wide range of water vapor supersaturations (S=0.068% to 1.27%). A comprehensive uncertainty analysis showed that the instrument calibration depends strongly on the applied particle generation techniques, Köhler model calculations, and water activity parameterizations (relative deviations in S up to 25%). Laboratory experiments and a comparison with other CCN instruments confirmed the high accuracy and precision of the calibration and measurement procedures developed and applied in this study.The mean CCN number concentrations (NCCN,S) observed in polluted mega-city air and biomass burning smoke (Beijing and Pearl River Delta, China) ranged from 1000 cm−3 at S=0.068% to 16 000 cm−3 at S=1.27%, which is about two orders of magnitude higher than in pristine air at remote continental sites (Swiss Alps, Amazonian rainforest). Effective average hygroscopicity parameters, κ, describing the influence of chemical composition on the CCN activity of aerosol particles were derived from the measurement data. They varied in the range of 0.3±0.2, were size-dependent, and could be parameterized as a function of organic and inorganic aerosol mass fraction. At low S (≤0.27%), substantial portions of externally mixed CCN-inactive particles with much lower hygroscopicity were observed in polluted air (fresh soot particles with κ≈0.01). Thus, the aerosol particle mixing state needs to be known for highly accurate predictions of NCCN,S. Nevertheless, the observed CCN number concentrations could be efficiently approximated using measured aerosol particle number size distributions and a simple κ-Köhler model with a single proxy for the effective average particle hygroscopicity. The relative deviations between observations and model predictions were on average less than 20% when a constant average value of κ=0.3 was used in conjunction with variable size distribution data. With a constant average size distribution, however, the deviations increased up to 100% and more. The measurement and model results demonstrate that the aerosol particle number and size are the major predictors for the variability of the CCN concentration in continental boundary layer air, followed by particle composition and hygroscopicity as relatively minor modulators. Depending on the required and applicable level of detail, the measurement results and parameterizations presented in this study can be directly implemented in detailed process models as well as in large-scale atmospheric and climate models for efficient description of the CCN activity of atmospheric aerosols.

Contributo sperimentale alla ridefinizione della diffusione in parete della camera di prova della galleria del vento Dallara

Joseph Nicolas Cugnot built the first primitive car in 1769 and approximately one hundred year later the first automotive race took place. Thanks to this, for the first time the aerodynamics principles began to be applied to cars. The aerodynamic study of a car is important to improve the performance on the road, or on the track. It purposely enhances the stability in the turns and increases the maximum velocity. However, it is also useful, decrease the fuel consumption, in order to reduce the pollution. Given that cars are a very complex body, the aerodynamic study cannot be conducted following an analytical method, but it is possible, in general, to choose between two different approaches: the numerical or the experimental one. The results of numerical studies depend on the computers’ potential and on the method use to implement the mathematical model. Today, the best way to perform an aerodynamic study is still experimental, which means that in the first phase of the design process the study is performed in a wind tunnel and in later phases directly on track. The automotive wind tunnels are singular mainly due to the test chamber, which typically contains a ground simulation system. The test chamber can have different types of walls: open walls, closed walls, adaptive walls or slotted walls. The best solution is to use the slotted walls because they minimize the interference between the walls and the streamlines, the interaction between the flow and the environment, and also to contain the overall costs. Furthermore, is necessary minimize the boundary layer at the walls, without accelerating the flow, in order to provide the maximum section of homogeneous flow. This thesis aims at redefining the divergent angle of the Dallara Automobili S.P.A. wind tunnel’s walls, in order to improve the overall homogeneity. To perform this study it was necessary to acquire the pressure data of the boundary layer, than it was created the profile of the boundary layer velocity and, to minimize the experimental errors, it was calculated the displacement thickness. The results obtained shows, even if the instrument used to the experiment was not the best one, that the boundary layer thickness could be minor in case of a low diffusion angle. So it is convenient to perform another experiment with a most sensitive instrument to verified what is the better wall configuration.

Dynamics of Bora wind over the Adriatic sea: atmospheric water balance and role of air-sea fluxes and orography

The Bora wind is a mesoscale phenomenon which typically affects the Adriatic Sea basin for several days each year, especially during winter. The Bora wind has been studied for its intense outbreak across the Dinaric Alps. The properties of the Bora wind are widely discussed in the literature and scientific papers usually focus on the eastern Adriatic coast where strong turbulence and severe gust intensity are more pronounced. However, the impact of the Bora wind can be significant also over Italy, not only in terms of wind speed instensity. Depending on the synoptic pressure pattern (cyclonic or anticyclonic Bora) and on the season, heavy snowfall, severe storms, storm surges and floods can occur along the Adriatic coast and on the windward flanks of the Apennines. In the present work five Bora cases that occurred in recent years have been selected and their evolution has been simulated with the BOLAM-MOLOCH model set, developed at ISAC-CNR in Bologna. Each case study has been addressed by a control run and by several sensitivity tests, performed with the purpose of better understanding the role played by air-sea latent and sensible heat fluxes. The tests show that the removal of the fluxes induces modifications in the wind approching the coast and a decrease of the total precipitation amount predicted over Italy. In order to assess the role of heat fluxes, further analysis has been carried out: column integrated water vapour fluxes have been computed along the Italian coastline and an atmospheric water balance has been evaluated inside a box volume over the Adriatic Sea. The balance computation shows that, although latent heat flux produces a significant impact on the precipitation field, its contribution to the balance is relatively minor. The most significant and lasting case study, that of February 2012, has been studied in more detail in order to explain the impressive drop in the total precipitation amount simulated in the sensitivity tests with removed heat fluxes with respect to the CNTRL run. In these experiments relative humidity and potential temperature distribution over different cross-sections have been examined. With respect to the CNTRL run a drier and more stable boundary layer, characterised by a more pronounced wind shear at the lower levels, has been observed to establish above the Adriatic Sea. Finally, in order to demonstrate that also the interaction of the Bora flow with the Apennines plays a crucial role, sensitivity tests varying the orography height have been considered. The results of such sensitivity tests indicate that the propagation of the Bora wind over the Adriatic Sea, and in turn its meteorological impact over Italy, is influenced by both the large air-sea heat fluxes and the interaction with the Apennines that decelerate the upstream flow.

A class of exact Navier-Stokes solutions for homogeneous flat-plate boundary layers and their linear stability

We introduce a new boundary layer formalism on the basis of which a class of exact solutions to the Navier–Stokes equations is derived. These solutions describe laminar boundary layer flows past a flat plate under the assumption of one homogeneous direction, such as the classical swept Hiemenz boundary layer (SHBL), the asymptotic suction boundary layer (ASBL) and the oblique impingement boundary layer. The linear stability of these new solutions is investigated, uncovering new results for the SHBL and the ASBL. Previously, each of these flows had been described with its own formalism and coordinate system, such that the solutions could not be transformed into each other. Using a new compound formalism, we are able to show that the ASBL is the physical limit of the SHBL with wall suction when the chordwise velocity component vanishes while the homogeneous sweep velocity is maintained. A corresponding non-dimensionalization is proposed, which allows conversion of the new Reynolds number definition to the classical ones. Linear stability analysis for the new class of solutions reveals a compound neutral surface which contains the classical neutral curves of the SHBL and the ASBL. It is shown that the linearly most unstable Görtler–Hämmerlin modes of the SHBL smoothly transform into Tollmien–Schlichting modes as the chordwise velocity vanishes. These results are useful for transition prediction of the attachment-line instability, especially concerning the use of suction to stabilize boundary layers of swept-wing aircraft.

Laminar and turbulent nozzle-jet flows and their acoustic near-field

We investigate numerically the effects of nozzle-exit flow conditions on the jet-flow development and the near-field sound at a diameter-based Reynolds number of Re D = 18 100 and Mach number Ma = 0.9. Our computational setup features the inclusion of a cylindrical nozzle which allows to establish a physical nozzle-exit flow and therefore well-defined initial jet-flow conditions. Within the nozzle, the flow is modeled by a potential flow core and a laminar, transitional, or developing turbulent boundary layer. The goal is to document and to compare the effects of the different jet inflows on the jet flow development and the sound radiation. For laminar and transitional boundary layers, transition to turbulence in the jet shear layer is governed by the development of Kelvin-Helmholtz instabilities. With the turbulent nozzle boundary layer, the jet flow development is characterized by a rapid changeover to a turbulent free shear layer within about one nozzle diameter. Sound pressure levels are strongly enhanced for laminar and transitional exit conditions compared to the turbulent case. However, a frequency and frequency-wavenumber analysis of the near-field pressure indicates that the dominant sound radiation characteristics remain largely unaffected. By applying a recently developed scaling procedure, we obtain a close match of the scaled near-field sound spectra for all nozzle-exit turbulence levels and also a reasonable agreement with experimental far-field data.

Pollen record from Golkow, Poland

Janczyk-Kopikowa (1966): The series of the organic deposits, developed in the vicinity of Golkow near Warsaw as oil shales and peats, was laid down in a grough valley and now rests on the deposits of the Middle Polish Glaciation (Riss). The organic deposits are overlain by the fluviale deposits of the North Polish Glaciation (Würm). The locality Golkow occurs beyond the extent of the continental glacier of this glaciation. Polen analysis completed by microfloristic examinations allows to determine the age of the organic series that is thought to be Eemian. The pollen diagram from Golkow does not call in question the stratigraphical position of the deposits investigated mainly due to its characteristic features such as minimum content of coniferous trees in the climatic optimum - about 5%, high percentage of Corylus - 77.5% and well developed phase of hornbeam. It may be well compared with other Eemian diagrams from the area of Poland and reveals much similar features. The development of vegetation at Golkow has depended upon the prevailing climate. At first, the cool climate brings about the development of plants having small thermal requirements. Here belong thin, park-like forests with pine and birch (Pinus, Betula) accompanied by the heliophilic plants such as Hippohäe and Ephedra. Improvement of climate that becomes warm and humid provides for development of deciduous forests prevailing in the climatic optimum, of the interglacial. Decrease of temperature causes a repeated change in the type of forest. This latter changes into coniferous forest with prevailing spruce (Picea) and fir (Abies) at the beginning, and then with pine (Pinus) and birch (Betula). During the Eemian Interglacial, the development of plants at Golkow terminates with a new and long-lasting predominance of pine-birch forests. However, such a longevity may be apparent only. Apparent character of this phenomenon is proved by a fact that the pollen spectra of the warm climatic periods have found their reflex in the oil shale that increased considerably slower than the layers off feebly decomposed peat evidencing the existence of cool pine-birch forests from the decline of the Interglacial. The water basin, in which the polen grains were laid down from surrounding plants is characterized by a calm sedimentation as proved by the occurrence of the oil shale. An insignificant water flow left behind some thin sand laminae. The not too deep basin becomes shallower owing to the growing water vegetation, and marshy vegetation. The growing of the plants causes a complete shallowing of the basin and formation of peat bog in situ, as proved by the peat beds occurring in the section. ---- Gadomska (1966): In the vicinity of Golków a series of organic deposits occurs amounting to 6.5-9.3 m in thickness, and consisting of oil shales, lacustrine silts and sands, as well as peats and peaty silts. The organic deposits fill up an old, small, but fairly deep lake basin, probably of finger-lake origin. It may be seen to-day as a slight lowering of the relief, filled up with soaked ground, stretching from north to south. On the basis of palaeobotanical examinations the organic deposits considered are of Eemian Interglacial age (Z. Janczyk-Kopikowa, 1063). The lower part of the organic series consists of a compact oil shale horizon, the maximum thickness of which may attain up to 8 m. The oil shales contain particularly in their upper part, numerous intercalations of arenaceous silts, dark grey or black in colour, or of sands mainly of lacustrine provenance. At the top of the oil shales are found peats, up to 2.5 m in thickness, covered by black, humus silts with numerous plant remains. The Eemian Interglacial deposits are covered by a series of fluviatile sands belonging partly to the Baltic Glaciation (bottom part of the series), partly to the Holocene (top part of the series). The thickness of the sands is 0.5-3.7 m. Higher up, there are found the Holocene and present-day deposits developed as clayey alluvion, or arenaceous slide rocks, or arenaceous-silty soil.