Streamline topology and dilute particle dynamics in a Karman vortex street flow

Three types of streamline topology in a Karman vortex street flow are shown under the variation of spatial parameters. For the motion of dilute particles in the Karman vortex street flow, there exist a route of bifurcation to a chaotic orbit and more attractors in a bifurcation diagram for the proportion of particle density to fluid density. Along with the increase of spatial parameters in the flow field, the bifurcation process is suspended, as well as more and more attractors emerge. In the motion of dilute particles, a drag term and gravity term dominate and result in the bifurcation phenomenon.

A high order accurate difference scheme for complex flow fields

A high order accurate finite difference method for direct numerical simulation of coherent structure in the mixing layers is presented. The reason for oscillation production in numerical solutions is analyzed, It is caused by a nonuniform group velocity of wavepackets. A method of group velocity control for the improvement of the shock resolution is presented. In numerical simulation the fifth-order accurate upwind compact difference relation is used to approximate the derivatives in the convection terms of the compressible N-S equations, a sixth-order accurate symmetric compact difference relation is used to approximate the viscous terms, and a three-stage R-K method is used to advance in time. In order to improve the shock resolution the scheme is reconstructed with the method of diffusion analogy which is used to control the group velocity of wavepackets. (C) 1997 Academic Press.

可压缩双组分混合层的直接数值模拟研究

通过数值求解可压缩Navier-Stokes方程模拟了三维可压缩双组分混合层问题,着重研究了双组分密度比和可压缩效应对混合层发展演化的影响.研究表明,随着密度比的增大,混合层扰动增长率降低.在混合层发展初期,可压缩性作用尤其是斜压效应对旋涡的发展起主导作用;在混合层发展后期,可压缩性影响减弱,旋涡的拉伸扭转效应占主导作用.同时,还分析讨论了可压缩混合层中湍流的转捩及拟序结构的演化过程.

Multiscale Method for Turbulent-Flow Computations

A main method of predicting turbulent flows is to solve LES equations, which was called traditional LES method. The traditional LES method solves the motions of large eddies of size larger than filtering scale An while modeling unresolved scales less than Delta_n. Hughes et al argued that many shortcomings of the traditional LES approaches were associated with their inabilities to successfully differentiate between large and small scales. One may guess that a priori scale-separation would be better, because it can predict scale-interaction well compared with posteriori scale-separation. To this end, a multi-scale method was suggested to perform scale-separation computation. The primary contents of the multiscale method are l) A space average is used to differentiate scale. 2) The basic equations include the large scale equations and fluctuation equations. 3) The large-scale equations and fluctuation equations are coupled through turbulent stress terms. We use the multiscale equations of n=2, i.e., the large and small scale (LSS) equations, to simulate 3-D evolutions of a channel flow and a planar mixing layer flow Some interesting results are given.

可压平面混合层稳定性分析及数值模拟

可压平面混合层是包含复杂多时空尺度运动的非定常流体力学部问题，具有深刻的理论意义和广泛的应用背景。针对该问题所涉及内容的多面性，本文的目的是，基于高精度、高分辨率数值算法的构造、发展和数值行为分析，采用线性稳定性分析和直接数值模拟方法。从理论和计算两方面集中研究压缩性效应、粘性效应、初值效应以及燃烧反应放热效应等对可压平面混合层早期稳定性行为和大尺度拟序涡结构非线性演化的影响。以混合层已有研究成果的分析和综述为开端，论文主体共包括四部分：第一部分是可压平面混合层时间/空间模式数值线性稳定性分析。实现了高精度对称紧致差分格式（SCD）对可压粘性扰动线性稳定性边值问题的求解，对导出的线性和非线性离散特征值问题，提出了两个高效局部解法。研究涉及二维/三维扰动波、无粘/粘性扰动波、特征函数和特征值谱、第一/第二模态、超声速快/慢模态、速度比和密度比等。验证了对流Mach数Mc为一个合理的压缩性参数。指出压缩性效应和粘性效应对最不稳定扰动波的波数（频率）和增长率呈相拟的抑制作用，且时间模式稳定性分析结果在许多方面是可信的。从随机和线性扰动场出发，采用高精度五阶迎风紧致和六阶对称紧致混合差分算法（UCD5/SCD6）对可压平面混合层的稳定性特征进行了直接数值模拟，揭示了初始主导线性扰动与一些实际涡结构非线性作用形态间的内在关联，印证了线性稳定性分析方法的合理性和有效性。第二部分是高精度迎风紧致差分格式（UCD）时空全离散数值行为分析。导出了其一维/二维一般色散表达式。研究表明，UCD格式在高波数区具有内在的全离散耗散和色散特性；其数值群速度的快/慢特征可因CFL数不同而改变；在稳定CFL数下简单附加人工粘性可强化UCD格式在高波数区的耗散量；提高时间精度可放宽稳定CFL数限制；UCD格式的二维全离散色散介质中存在三个不同性质的数值波，其全离散稳定性由数值声波主控。第三部分实现了高精度UCD5/SCD6差分算法对空间发展可压平面混合层的直接数值模拟。通过亚谐扰动波的个数和扰动频率的控制，捕捉到了基频涡的饱和、一次和二次对并等现象，显示了大尺度涡结构与入中初始扰动方式之间的内在联系。利用参数Mc观察了压缩性效应对大尺度涡空间演化及其相互作用的影响。第四部分实现了高精度UCD5/SCD6差分算法对非预混扩散火焰化学反应平面混合层的直接数值模拟。研究指出，放热效应可抑制和延迟涡的形成，使基频涡卷拉伸甚至丧失，混合层Reynolds 应力ρu'v'和流向速度波动关联项u'v'下降，以致涡结构与外流动量交换和标量输运减少，脉动输运能力被削弱，从而混合效率、产物生成率和混合层增长率下降，放热主要通过膨胀效应和斜压效应来抑制大尺度涡的演化。

Velocity resolved - scalar modeled simulations of high Schmidt number turbulent transport

The objective of this thesis is to develop a framework to conduct velocity resolved - scalar modeled (VR-SM) simulations, which will enable accurate simulations at higher Reynolds and Schmidt (Sc) numbers than are currently feasible. The framework established will serve as a first step to enable future simulation studies for practical applications. To achieve this goal, in-depth analyses of the physical, numerical, and modeling aspects related to Sc>>1 are presented, specifically when modeling in the viscous-convective subrange. Transport characteristics are scrutinized by examining scalar-velocity Fourier mode interactions in Direct Numerical Simulation (DNS) datasets and suggest that scalar modes in the viscous-convective subrange do not directly affect large-scale transport for high Sc. Further observations confirm that discretization errors inherent in numerical schemes can be sufficiently large to wipe out any meaningful contribution from subfilter models. This provides strong incentive to develop more effective numerical schemes to support high Sc simulations. To lower numerical dissipation while maintaining physically and mathematically appropriate scalar bounds during the convection step, a novel method of enforcing bounds is formulated, specifically for use with cubic Hermite polynomials. Boundedness of the scalar being transported is effected by applying derivative limiting techniques, and physically plausible single sub-cell extrema are allowed to exist to help minimize numerical dissipation. The proposed bounding algorithm results in significant performance gain in DNS of turbulent mixing layers and of homogeneous isotropic turbulence. Next, the combined physical/mathematical behavior of the subfilter scalar-flux vector is analyzed in homogeneous isotropic turbulence, by examining vector orientation in the strain-rate eigenframe. The results indicate no discernible dependence on the modeled scalar field, and lead to the identification of the tensor-diffusivity model as a good representation of the subfilter flux. Velocity resolved - scalar modeled simulations of homogeneous isotropic turbulence are conducted to confirm the behavior theorized in these a priori analyses, and suggest that the tensor-diffusivity model is ideal for use in the viscous-convective subrange. Simulations of a turbulent mixing layer are also discussed, with the partial objective of analyzing Schmidt number dependence of a variety of scalar statistics. Large-scale statistics are confirmed to be relatively independent of the Schmidt number for Sc>>1, which is explained by the dominance of subfilter dissipation over resolved molecular dissipation in the simulations. Overall, the VR-SM framework presented is quite effective in predicting large-scale transport characteristics of high Schmidt number scalars, however, it is determined that prediction of subfilter quantities would entail additional modeling intended specifically for this purpose. The VR-SM simulations presented in this thesis provide us with the opportunity to overlap with experimental studies, while at the same time creating an assortment of baseline datasets for future validation of LES models, thereby satisfying the objectives outlined for this work.

SiO2/Ta interface reaction in magnetic multilayers and its influence on Ta buffer layers

Ta is often used as a buffer layer in magnetic multilayers. In this study, Ta/Ni81Fe19/Ta multilayers were deposited by magnetron sputtering on sing-crystal Si with a 300-nm-thick SiO2 film. The composition and chemical states at the interface region of SiO2/Ta were studied using the X-ray photoelectron spectroscopy (XPS) and peak decomposition technique. The results show that there is an 'inter-mixing layer" at the SiO2/Ta interface due to a thermodynamically favorable reaction: 15 SiO2 + 37 Ta = 6 Ta2O5 + 5 Ta5Si3. Therefore, the Ta buffer layer thickness used to induce NiFe (111) texture increases.

Alguns aspectos da modelação ecológica na costa portuguesa (Aveiro) : efeitos físicos na distribuição de nutrientes e biomassa fitoplanctónica

O presente trabalho propõe estudar alguns aspectos da modelação ecológica na Costa Portuguesa, nomeadamente o efeito físico na distribuição de nutrientes e biomassa fitoplanctónica. O principal propósito foi implementar e validar um modelo acoplado tridimensional físico e ecológico, para a costa portuguesa, e aplicá-lo numa área limitada, a norte, por Vila do Conde e a sul pela Figueira da Foz centrada na região de Aveiro, para estudar a distribuição, vertical e horizontal, de temperatura, nutrientes e biomassa fitoplanctónica. A região em estudo está situada na costa oeste da Península Ibérica e faz parte da Região de Afloramento do Atlântico Norte. É caracterizada por condições meteorológicas onde a predominância e prevalência, durante uma grande parte do ano, de ventos de norte/noroeste, constitui um dos principais elementos forçadores do transporte para o largo das águas costeiras e consequente subida das águas mais frias e profundas, ricas em nutrientes. A estas condições juntam-se as boas condições de luminosidade e temperatura essenciais ao desenvolvimento fitoplanctónico, que servirá de alimento às espécies marinhas. Este facto, torna esta região, uma zona de elevada riqueza biológica favorável ao desenvolvimento de várias espécies marinhas, transformando este local num ecossistema de forte produtividade. O modelo foi calibrado e validado para a área em estudo e simulou com sucesso, a resposta do sistema à situação de ventos favoráveis ao afloramento costeiro para a região em estudo. Quando comparada com dados observados, os resultados mostram que o modelo é capaz de prever satisfatoriamente as distribuições superficiais e na coluna de água: da temperatura, dos nutrientes, do oxigénio e da clorofila-a. Os resultados evidenciam o crucial papel desempenhado pelos processos físicos no aumento de fitoplâncton que ocorre ao longo de uma estreita área da costa norte portuguesa, mostrando a estreita ligação entre a distribuição costeira de fitoplâncton e a distribuição costeira de temperatura, à superfície. A produtividade do fitoplâncton que ocorre nas águas costeiras não só é atribuída à disponibilidade de nutrientes mas também à intensidade luminosa. A luz e a intensidade dos ventos de norte/noroeste são os factores chave no controlo dos blooms de fitoplâncton observados nesta região de afloramento, sobretudo no Verão, nomeadamente na camada de mistura pouco profunda e junto à nutriclina. Deste modo os modelos numéricos em associação com dados in situ e imagens de satélite poderão ser considerados uma excelente ferramenta para a análise e previsão de cenários, presentes e futuros, de acções praticadas sobre o meio ambiente.

Seasonality and extent of extratropical TST derived from in-situ CO measurements during SPURT

We present airborne in-situ trace gas measurements which were performed on eight campaigns between November 2001 and July 2003 during the SPURT-project (SPURenstofftransport in der Tropopausenregion, trace gas transport in the tropopause region). The measurements on a quasi regular basis allowed an overview of the seasonal variations of the trace gas distribution in the tropopause region over Europe from 35°-75°N to investigate the influence of transport and mixing across the extratropical tropopause on the lowermost stratosphere. From the correlation of CO and O3 irreversible mixing of tropospheric air into the lowermost stratosphere is identified. The CO distribution indicates that transport and subsequent mixing of tropospheric air across the extratropical tropopause predominantly affects a layer, which closely follows the shape of the local tropopause. In addition, the seasonal cycle of CO2 illustrates the strong coupling of that layer to the extratropical troposphere. Both, horizontal gradients of CO on isentropes as well as the CO-O3-distribution in the lowermost stratosphere reveal that the influence of quasi-horizontal transport and subsequent mixing weakens with distance from the local tropopause. The mixing layer extends to about 25 K in potential temperature above the local tropopause exhibiting only a weak seasonality. However, at large distances from the tropopause a significant influence of tropospheric air is still evident. The relation between N2O and CO2 indicates that a significant contribution of air originating from the tropical tropopause contributes to the background air in the extratropical lowermost stratosphere.

Daytime CO2 urban surface fluxes from airborne measurements, eddy-covariance observations and emissions inventory in Greater London

Airborne measurements within the urban mixing layer (360 m) over Greater London are used to quantify CO2 emissions at the meso-scale. Daytime CO2 fluxes, calculated by the Integrative Mass Boundary Layer (IMBL) method, ranged from 46 to 104 μmol CO2 m−2 s−1 for four days in October 2011. The day-to-day variability of IMBL fluxes is at the same order of magnitude as for surface eddy-covariance fluxes observed in central London. Compared to fluxes derived from emissions inventory, the IMBL method gives both lower (by −37%) and higher (by 19%) estimates. The sources of uncertainty of applying the IMBL method in urban areas are discussed and guidance for future studies is given.

The dynamics and relationships of precipitation, temperature and convection boundaries in the dayside auroral ionosphere

A continuous band of high ion temperature, which persisted for about 8 h and zigzagged north-south across more than five degrees in latitude in the dayside (07:00– 15:00MLT) auroral ionosphere, was observed by the EISCAT VHF radar on 23 November 1999. Latitudinal gradients in the temperature of the F-region electron and ion gases (Te and Ti , respectively) have been compared with concurrent observations of particle precipitation and field-perpendicular convection by DMSP satellites, in order to reveal a physical explanation for the persistent band of high Ti , and to test the potential role of Ti and Te gradients as possible markers for the open-closed field line boundary. The north/south movement of the equatorward Ti boundary was found to be consistent with the contraction/expansion of the polar cap due to an unbalanced dayside and nightside reconnection. Sporadic intensifications in Ti , recurring on _10-min time scales, indicate that frictional heating was modulated by time-varying reconnection, and the band of high Ti was located on open flux. However, the equatorward Ti boundary was not found to be a close proxy of the open-closed boundary. The closest definable proxy of the open-closed boundary is the magnetosheath electron edge observed by DMSP. Although Te appears to be sensitive to magnetosheath electron fluxes, it is not found to be a suitable parameter for routine tracking of the open-closed boundary, as it involves case dependent analysis of the thermal balance. Finally, we have documented a region of newly-opened sunward convecting flux. This region is situated between the convection reversal boundary and the magnetosheath electron edge defining the openclosed boundary. This is consistent with a delay of several minutes between the arrival of the first (super-Alfv´enic) magnetosheath electrons and the response in the ionospheric convection, conveyed to the ionosphere by the interior Alfv´en wave. It represents a candidate footprint of the low-latitude boundary mixing layer on sunward convecting open flux

Calibration Strategies Designed For Long Term Sediment Transport Modelling

The presented work deals with the calibration of a 2D numerical model for the simulation of long term bed load transport. A settled basin along an alpine stream was used as a case study. The focus is to parameterise the used multi fractional transport model such that a dynamically balanced behavior regarding erosion and deposition is reached. The used 2D hydrodynamic model utilizes a multi-fraction multi-layer approach to simulate morphological changes and bed load transport. The mass balancing is performed between three layers: a top mixing layer, an intermediate subsurface layer and a bottom layer. Using this approach bears computational limitations in calibration. Due to the high computational demands, the type of calibration strategy is not only crucial for the result, but as well for the time required for calibration. Brute force methods such as Monte Carlo type methods may require a too large number of model runs. All here tested calibration strategies used multiple model runs utilising the parameterization and/or results from previous run. One concept was to reset to initial bed elevations after each run, allowing the resorting process to convert to stable conditions. As an alternative or in combination, the roughness was adapted, based on resulting nodal grading curves, from the previous run. Since the adaptations are a spatial process, the whole model domain is subdivided in homogeneous sections regarding hydraulics and morphological behaviour. For a faster optimization, the adaptation of the parameters is made section wise. Additionally, a systematic variation was done, considering results from previous runs and the interaction between sections. The used approach can be considered as similar to evolutionary type calibration approaches, but using analytical links instead of random parameter changes.

Avaliação de características dos regimes de umidade na flona de Caxiuanã-PA durante o experimento COBRA-PARÁ

Procura-se investigar a validade de um método de classificação de regimes de umidade, baseado na caracterização de diferentes "estados" da Camada Limite Atmosférica Tropical (CLAT), acima de uma área de floresta, de acordo com a metodologia proposta por Mahrt (1991). Para essas análises foram utilizados dados de radiossondagens e de uma torre micrometeorológica, coletados durante o período menos chuvoso da região, obtidos durante o experimento "COBRA-PARÁ" (realizado no período de 30/10 a 15/11 de 2006). A análise dos regimes de umidade consiste na representação em espaço de fase dos dados disponíveis da razão de Bowen (β), em função do parâmetro -h/L (onde h é a altura da camada de mistura turbulenta e L é o comprimento de Obukhov). Dependendo da localização dos dados nesse espaço foi possível caracterizar as seguintes classes: classe I - ar seco e instável; classe II - vento seco predominante; classe III - vento úmido; classe IV - condição úmida e instável; classe V - condensação de vapor d'água na superfície; classe VI - condição estável dominante; e classe VII -formação de orvalho induzido por radiação noturna resfriando a superfície. Das classes mencionadas, aquelas mais freqüentemente observadas em Caxiuanã, foram as III, IV e VI.

Flugzeuggetragene Spurengasmessungen in der Tropopausenregion

ZusammenfassungIm Rahmen der EU-Projekte POLSTAR (Polar StratosphericAerosol Experiment) und STREAM (Stratosphere TroposphereExperiment by Airborne Measurements) wurden flugzeuggetragene Spurengasmessungen in verschiedenen geografischen Breiten durchgeführt. Zwei Messkampagnen fanden im Januar (POLSTAR 97) und März (STREAM 97) 1997 über Kiruna (Schweden, 67°N, 20°O) statt, eine Kampagne wurde im Juli 1998 von Timmins (Kanada, 47°N, 81°W) aus durchgeführt (STREAM 98).CO und N2O wurden mittels TDLAS (Tunable Diode Laser Absorption Spectroscopy) nachgewiesen, CO2 wurde mit einem modifizierten kommerziellen Messgerät breitbandig gemessen. Zur Untersuchung von Mischungsvorgängen in der Tropopausenregion wurden Korrelationenzwischen CO, O3, N2O, CO2 und NOy herangezogen.Dabei konnte festgestellt werden, dass im Winter in der untersten Stratosphäre Mischung mit troposphärischen Luftmassen durch isentropen Transport im Bereich der Polarfront bis zu potentiellen Temperaturen von Theta=335K auftritt. Im Sommer lässt sich Mischung mit troposphärischen Luftmassen bis mindestens Theta=360K nachweisen, die effektiver als im Winter abläuft.Exemplarisch kann an einem Flug gezeigt werden, dassbis Theta=349K die Ausbildung der Mischungsschicht durchLuftmassenaustausch an der Polarfront verursacht wird, während oberhalb von Theta=349K Signaturen subtropischertroposphärischer Luftmassen gefunden werden.Der stratosphärische Hintergrund wird mitbestimmt durch photochemisch gealterte Luftmassen, die ihren Ursprung höchstwahrscheinlich in der Vortexregion des vorangegangenen Winters haben.

Tropospheric carboxylic acids in tropical rainforest environments

Der Übergang der Grenzschicht von stark ozeanisch auf kontinental beeinflusst wurde in 2 tropischen Küstenwaldgebieten Amerikas untersucht, wo Luftmassen vom Meer kommend über den Kontinent transportiert werden.Zwei Feldkampagnen wurden durchgeführt; in Costa Rica (CR; 07.1996) und in Surinam (04.1998). In CR wurde im nordöstlichen Flachgebiet (etwa 10°25' N; 84°W) Regenwasser gesammelt, in dem später Carboxylate, anorganischen Anionen, Ca2+, K+, NH4+ und Mg2+ gemessen wurden. Die Proben wurden an 5 verschiedenen Stellen entlang der Windrichtung 1, 20, 60, 60 und 80 km von der Küste gesammelt. In Surinam (Sipaliwini, 2°02' N, 56°08' W) wurden organischen Säuren aus der Gasphase gesammelt etwa 550 km von der Küste entfernt, sowohl wie O3 und CO. Die Proben wurden mittels Ionenchromatographie und Kapillarelektrophorese analysiert. Morgendliche Einmischung der nächtlichen residualen Schicht und Luft der unteren freien Troposphäre war Hauptquelle für HCOOH und CH3COOH in der Tagesmischschicht. Es wurde gezeigt, dass lokale Produktion dieser Säuren durch chemische Reaktionen eine kleine Rolle gespielt hat und dass direkte Emission vernachlässigbar war.Aus den beiden Feldkampagnen folgt, dass die Konzentrationen der sekundären Verbindungen HCOOH, CH3COOH, Ozon und CO in der Tagesmischschicht von Importen bestimmt wurden, was gilt für Regen- und Trockenzeit bis zu Entfernungen von 550 km zur Küste.