Expansion and tissue infiltration of an allospecific CD4+CD25+CD45RO+IL-7Ralphahigh cell population in solid organ transplant recipients.

It has been recently shown (Seddiki, N., B. Santner-Nanan, J. Martinson, J. Zaunders, S. Sasson, A. Landay, M. Solomon, W. Selby, S.I. Alexander, R. Nanan, et al. 2006. J. Exp. Med. 203:1693-1700.) that the expression of interleukin (IL) 7 receptor (R) alpha discriminates between two distinct CD4 T cell populations, both characterized by the expression of CD25, i.e. CD4 regulatory T (T reg) cells and activated CD4 T cells. T reg cells express low levels of IL-7Ralpha, whereas activated CD4 T cells are characterized by the expression of IL-7Ralpha(high). We have investigated the distribution of these two CD4 T cell populations in 36 subjects after liver and kidney transplantation and in 45 healthy subjects. According to a previous study (Demirkiran, A., A. Kok, J. Kwekkeboom, H.J. Metselaar, H.W. Tilanus, and L.J. van der Laan. 2005. Transplant. Proc. 37:1194-1196.), we observed that the T reg CD25(+)CD45RO(+)IL-7Ralpha(low) cell population was reduced in transplant recipients (P < 0.00001). Interestingly, the CD4(+)CD25(+)CD45RO(+)IL-7Ralpha(high) cell population was significantly increased in stable transplant recipients compared with healthy subjects (P < 0.00001), and the expansion of this cell population was even greater in patients with documented humoral chronic rejection compared with stable transplant recipients (P < 0.0001). The expanded CD4(+)CD25(+)CD45RO(+)IL-7Ralpha(high) cell population contained allospecific CD4 T cells and secreted effector cytokines such as tumor necrosis factor alpha and interferon gamma, thus potentially contributing to the mechanisms of chronic rejection. More importantly, CD4(+)IL-7Ralpha(+)and CD25(+)IL-7Ralpha(+) cells were part of the T cell population infiltrating the allograft of patients with a documented diagnosis of chronic humoral rejection. These results indicate that the CD4(+)CD25(+)IL-7Ralpha(+) cell population may represent a valuable, sensitive, and specific marker to monitor allospecific CD4 T cell responses both in blood and in tissues after organ transplantation.

Natural selection in a postglacial range expansion: the case of the colour cline in the European barn owl.

Gradients of variation-or clines-have always intrigued biologists. Classically, they have been interpreted as the outcomes of antagonistic interactions between selection and gene flow. Alternatively, clines may also establish neutrally with isolation by distance (IBD) or secondary contact between previously isolated populations. The relative importance of natural selection and these two neutral processes in the establishment of clinal variation can be tested by comparing genetic differentiation at neutral genetic markers and at the studied trait. A third neutral process, surfing of a newly arisen mutation during the colonization of a new habitat, is more difficult to test. Here, we designed a spatially explicit approximate Bayesian computation (ABC) simulation framework to evaluate whether the strong cline in the genetically based reddish coloration observed in the European barn owl (Tyto alba) arose as a by-product of a range expansion or whether selection has to be invoked to explain this colour cline, for which we have previously ruled out the actions of IBD or secondary contact. Using ABC simulations and genetic data on 390 individuals from 20 locations genotyped at 22 microsatellites loci, we first determined how barn owls colonized Europe after the last glaciation. Using these results in new simulations on the evolution of the colour phenotype, and assuming various genetic architectures for the colour trait, we demonstrate that the observed colour cline cannot be due to the surfing of a neutral mutation. Taking advantage of spatially explicit ABC, which proved to be a powerful method to disentangle the respective roles of selection and drift in range expansions, we conclude that the formation of the colour cline observed in the barn owl must be due to natural selection.

Adjustable Suture Removal: Increasing Aqueous Flow in "Safe" Trabeculectomy.

Purpose: Aqueous flow through trabeculectomy blebs has been suggested to influence filtration bleb survival. We investigated the relationship between the requirement to increase aqueous flow via adjustable suture removal and surgical outcomes following "safe trabeculectomy" with mitomycin C (MMC). Methods: 62 consecutive eyes of 53 patients underwent fornix based trabeculectomy with adjustable sutures, intraoperative MMC and intensive postoperative steroids. Subconjunctival antimetabolite injections and bleb needlings were administered according to bleb vascularity and IOP trends. Main outcome measures were: success rates (definition: IOP≤21mmHg and 20% IOP reduction); number of antimetabolite injections; bleb needlings; number of of eyes recommencing glaucoma medications and complications. Results: Mean age was 70.4±16.0 years (mean± SD); mean preoperative IOP was 24.5±9.1 mmHg and decreased to 12.3±8.9mmHg postoperatively. Mean number of sutures was 2.6 ± 0.7. Eyes were divided into 2 groups in relation to the number of sutures removed. The number of subconjunctival MMC injections required for those requiring 2 suture removals was significantly greater than those requiring 1 suture removal (p<0.05) The number of needlings and 5FU injections also increased but did not reach significance (p=0.09 and p=0.34 respectively). Least-squared linear regression analysis showed the number of needlings required had a statistically significant (p=0.05) trend with respect to time elapsed between surgery and first suture removal. No other interventions had significant trends. Mean time between surgery and suture removal was: 4.2±9.2 weeks (suture #1) and 5.7±9.7 weeks (suture#2). Antiglaucoma medication was restarted in only 5 eyes. Postoperative complications were infrequent: Seidel (3.2%), peripheral choroidal effusions at any time (3.2%), and shallow anterior chamber (1.6%). Conclusion: Eyes requiring a greater number of suture removals required a significantly greater number of antifibrosis interventions. The time elapsed before suture removal was inversely related to the number of postoperative needlings, suggesting these eyes may have decreased aqueous production and therefore require aggressive post-operative management to prevent bleb failure.

Estimating traveltimes and groundwater flow patterns using 3D time-lapse crosshole ERT imaging of electrical resistivity fluctuations induced by infiltrating river water

The infiltration of river water into aquifers is of high relevance to drinking-water production and is a key driver of biogeochemical processes in the hyporheic and riparian zone, but the distribution and quantification of the infiltrating water are difficult to determine using conventional hydrological methods (e.g., borehole logging and tracer tests). By time-lapse inverting crosshole ERT (electrical resistivity tomography) monitoring data, we imaged groundwater flow patterns driven by river water infiltrating a perialpine gravel aquifer in northeastern Switzerland. This was possible because the electrical resistivity of the infiltrating water changed during rainfall-runoff events. Our time-lapse resistivity models indicated rather complex flow patterns as a result of spatially heterogeneous bank filtration and aquifer heterogeneity. The upper part of the aquifer was most affected by the river infiltrate, and the highest groundwater velocities and possible preferential flow occurred at shallow to intermediate depths. Time series of the reconstructed resistivity models matched groundwater electrical resistivity data recorded on borehole loggers in the upper and middle parts of the aquifer, whereas the resistivity models displayed smaller variations and delayed responses with respect to the logging data. in the lower part. This study demonstrated that crosshole ERT monitoring of natural electrical resistivity variations of river infiltrate could be used to image and quantify 3D bank filtration and aquifer dynamics at a high spatial resolution.

Generation and interconnection capacity expansion in cross-border electricity markets : the need for policy coordination

Electricity is a strategic service in modern societies. Thus, it is extremely important for governments to be able to guarantee an affordable and reliable supply, which depends to a great extent on an adequate expansion of the generation and transmission capacities. Cross- border integration of electricity markets creates new challenges for the regulators, since the evolution of the market is now influenced by the characteristics and policies of neighbouring countries. There is still no agreement on why and how regions should integrate their electricity markets. The aim of this thesis is to improve the understanding of integrated electricity markets and how their behaviour depends on the prevailing characteristics of the national markets and the policies implemented in each country. We developed a simulation model to analyse under what circumstances integration is desirable. This model is used to study three cases of interconnection between two countries. Several policies regarding interconnection expansion and operation, combined with different generation capacity adequacy mechanisms, are evaluated. The thesis is composed of three papers. The first paper presents a detailed description of the model and an analysis of the case of Colombia and Ecuador. It shows that market coupling can bring important benefits, but the relative size of the countries can lead to import dependency issues in the smaller country. The second paper compares the case of Colombia and Ecuador with the case of Great Britain and France. These countries are significantly different in terms of electricity sources, hydro- storage capacity, complementarity and demand growth. We show that complementarity is essential in order to obtain benefits from integration, while higher demand growth and hydro- storage capacity can lead to counterintuitive outcomes, thus complicating policy design. In the third paper, an extended version of the model presented in the first paper is used to analyse the case of Finland and its interconnection with Russia. Different trading arrangements are considered. We conclude that unless interconnection capacity is expanded, the current trading arrangement, where a single trader owns the transmission rights and limits the flow during peak hours, is beneficial for Finland. In case of interconnection expansion, market coupling would be preferable. We also show that the costs of maintaining a strategic reserve in Finland are justified in order to limit import dependency, while still reaping the benefits of interconnection. In general, we conclude that electricity market integration can bring benefits if the right policies are implemented. However, a large interconnection capacity is only desirable if the countries exhibit significant complementarity and trust each other. The outcomes of policies aimed at guaranteeing security of supply at a national level can be quite counterintuitive due to the interactions between neighbouring countries and their effects on interconnection and generation investments. Thus, it is important for regulators to understand these interactions and coordinate their decisions in order to take advantage of the interconnection without putting security of supply at risk. But it must be taken into account that even when integration brings benefits to the region, some market participants lose and might try to hinder the integration process. -- Dans les sociétés modernes, l'électricité est un service stratégique. Il est donc extrêmement important pour les gouvernements de pouvoir garantir la sécurité d'approvisionnement à des prix abordables. Ceci dépend en grande mesure d'une expansion adéquate des capacités de génération et de transmission. L'intégration des marchés électriques pose des nouveaux défis pour les régulateurs, puisque l'évolution du marché est maintenant influencée par les caractéristiques et les politiques des pays voisins. Il n'est pas encore claire pourquoi ni comment les marches électriques devraient s'intégrer. L'objectif de cette thèse est d'améliorer la compréhension des marchés intégrés d'électricité et de leur comportement en fonction des caractéristiques et politiques de chaque pays. Un modèle de simulation est proposé pour étudier les conditions dans lesquelles l'intégration est désirable. Ce modèle est utilisé pour étudier trois cas d'interconnexion entre deux pays. Plusieurs politiques concernant l'expansion et l'opération de l'interconnexion, combinées avec différents mécanismes de rémunération de la capacité, sont évalués. Cette thèse est compose de trois articles. Le premier présente une description détaillée du modèle et une analyse du cas de la Colombie et de l'Equateur. Il montre que le couplage de marchés peut amener des bénéfices importants ; cependant, la différence de taille entre pays peut créer des soucis de dépendance aux importations pour le pays le plus petit. Le second papier compare le cas de la Colombie et l'Equateur avec le cas de la Grande Bretagne et de la France. Ces pays sont très différents en termes de ressources, taille des réservoirs d'accumulation pour l'hydro, complémentarité et croissance de la demande. Nos résultats montrent que la complémentarité joue un rôle essentiel dans l'obtention des bénéfices potentiels de l'intégration, alors qu'un taux élevé de croissance de la demande, ainsi qu'une grande capacité de stockage, mènent à des résultats contre-intuitifs, ce qui complique les décisions des régulateurs. Dans le troisième article, une extension du modèle présenté dans le premier article est utilisée pour analyser le cas de la Finlande et de la Russie. Différentes règles pour les échanges internationaux d'électricité sont considérées. Nos résultats indiquent qu'à un faible niveau d'interconnexion, la situation actuelle, où un marchand unique possède les droits de transmission et limite le flux pendant les heures de pointe, est bénéfique pour la Finlande. Cependant, en cas d'expansion de la capacité d'interconnexion, «market coupling» est préférable. préférable. Dans tous les cas, la Finlande a intérêt à garder une réserve stratégique, car même si cette politique entraine des coûts, elle lui permet de profiter des avantages de l'intégration tout en limitant ca dépendance envers les importations. En général, nous concluons que si les politiques adéquates sont implémentées, l'intégration des marchés électriques peut amener des bénéfices. Cependant, une grande capacité d'interconnexion n'est désirable que si les pays ont une complémentarité importante et il existe une confiance mutuelle. Les résultats des politiques qui cherchent à préserver la sécurité d'approvisionnement au niveau national peuvent être très contre-intuitifs, étant données les interactions entre les pays voisins et leurs effets sur les investissements en génération et en interconnexion. Il est donc très important pour les régulateurs de comprendre ces interactions et de coordonner décisions à fin de pouvoir profiter de l'interconnexion sans mettre en danger la sécurité d'approvisionnement. Mais il faut être conscients que même quand l'intégration amène de bénéfices pour la région, certains participants au marché sont perdants et pourraient essayer de bloquer le processus d'intégration.

Isolation, expansion and differentiation of cellular progenitors obtained from dental pulp of agouti (Dasyprocta prymnolopha Wagler, 1831)

Abstract: The study aimed to isolate, expand, differentiate and characterize progenitor cells existent in the dental pulp of agouti. The material was washed with PBS solution and dissociated mechanically with the aid of a scalpel blade on plates containing culture medium D-MEM/F-12, and incubated at 5% CO2-37⁰C. The growth curve, CFU assay, osteogenic/adipogenic differentiation and characterization were obtained from the isolation. The cells began to be released from the explant tissue around the 7th day of culture. By day 22 of culture, cells reached 80% confluence. At the UFC test, 81 colonies were counted with 12 days of cultivation. The growth curves before and after freezing showed a regular growth with intense proliferation and clonogenic potential. The cell differentiation showed formation of osteoblasts and fat in culture, starting at 15 days of culture in a specific medium. Flow cytometry (FACs) was as follows: CD34 (positive), CD14 (negative), CD45 (negative), CD73 (positive), CD79 (negative), CD90 (positive), CD105 (positive), demonstrating high specificity and commitment of isolated cells with mesenchymal stem cells strains. These results suggest the existence of a cell population of stem cells with mesenchymal features from the isolated tissue in the explants of agouti dental pulp, a potential model for study of stem cell strains obtained from the pulp tissue.

Turbulent shallow-water model for orographic subgrid-scale perturbations

A parallel pseudo-spectral method for the simulation in distributed memory computers of the shallow-water equations in primitive form was developed and used on the study of turbulent shallow-waters LES models for orographic subgrid-scale perturbations. The main characteristics of the code are: momentum equations integrated in time using an accurate pseudo-spectral technique; Eulerian treatment of advective terms; and parallelization of the code based on a domain decomposition technique. The parallel pseudo-spectral code is efficient on various architectures. It gives high performance onvector computers and good speedup on distributed memory systems. The code is being used for the study of the interaction mechanisms in shallow-water ows with regular as well as random orography with a prescribed spectrum of elevations. Simulations show the evolution of small scale vortical motions from the interaction of the large scale flow and the small-scale orographic perturbations. These interactions transfer energy from the large-scale motions to the small (usually unresolved) scales. The possibility of including the parametrization of this effects in turbulent LES subgrid-stress models for the shallow-water equations is addressed.

Acute extracellular fluid volume changes increase ileocolonic resistance to saline flow in anesthetized dogs

We determined the effect of acute extracellular fluid volume changes on saline flow through 4 gut segments (ileocolonic, ileal, ileocolonic sphincter and proximal colon), perfused at constant pressure in anesthetized dogs. Two different experimental protocols were used: hypervolemia (iv saline infusion, 0.9% NaCl, 20 ml/min, volume up to 5% body weight) and controlled hemorrhage (up to a 50% drop in mean arterial pressure). Mean ileocolonic flow (N = 6) was gradually and significantly decreased during the expansion (17.1%, P<0.05) and expanded (44.9%, P<0.05) periods while mean ileal flow (N = 7) was significantly decreased only during the expanded period (38%, P<0.05). Mean colonic flow (N = 7) was decreased during expansion (12%, P<0.05) but returned to control levels during the expanded period. Mean ileocolonic sphincter flow (N = 6) was not significantly modified. Mean ileocolonic flow (N = 10) was also decreased after hemorrhage (retracted period) by 17% (P<0.05), but saline flow was not modified in the other separate circuits (N = 6, 5 and 4 for ileal, ileocolonic sphincter and colonic groups, respectively). The expansion effect was blocked by atropine (0.5 mg/kg, iv) both on the ileocolonic (N = 6) and ileal (N = 5) circuits. Acute extracellular fluid volume retraction and expansion increased the lower gastrointestinal resistances to saline flow. These effects, which could physiologically decrease the liquid volume being supplied to the colon, are possible mechanisms activated to acutely balance liquid volume deficit and excess.

Genesis and Evolution of Brittle Structures In Southwestern Finland and Western South Africa. Insights into Fault Reactivation, Fluid Flow and Structural Maturity in Precambrian Cratons

The bedrock of old crystalline cratons is characteristically saturated with brittle structures formed during successive superimposed episodes of deformation and under varying stress regimes. As a result, the crust effectively deforms through the reactivation of pre-existing structures rather than by through the activation, or generation, of new ones, and is said to be in a state of 'structural maturity'. By combining data from Olkiluoto Island, southwestern Finland, which has been investigated as the potential site of a deep geological repository for high-level nuclear waste, with observations from southern Sweden, it can be concluded that the southern part of the Svecofennian shield had already attained structural maturity during the Mesoproterozoic era. This indicates that the phase of activation of the crust, i.e. the time interval during which new fractures were generated, was brief in comparison to the subsequent reactivation phase. Structural maturity of the bedrock was also attained relatively rapidly in Namaqualand, western South Africa, after the formation of first brittle structures during Neoproterozoic time. Subsequent brittle deformation in Namaqualand was controlled by the reactivation of pre-existing strike-slip faults.In such settings, seismic events are likely to occur through reactivation of pre-existing zones that are favourably oriented with respect to prevailing stresses. In Namaqualand, this is shown for present day seismicity by slip tendency analysis, and at Olkiluoto, for a Neoproterozoic earthquake reactivating a Mesoproterozoic fault. By combining detailed field observations with the results of paleostress inversions and relative and absolute time constraints, seven distinctm superimposed paleostress regimes have been recognized in the Olkiluoto region. From oldest to youngest these are: (1) NW-SE to NNW-SSE transpression, which prevailed soon after 1.75 Ga, when the crust had sufficiently cooled down to allow brittle deformation to occur. During this phase conjugate NNW-SSE and NE-SW striking strike-slip faults were active simultaneous with reactivation of SE-dipping low-angle shear zones and foliation planes. This was followed by (2) N-S to NE-SW transpression, which caused partial reactivation of structures formed in the first event; (3) NW-SE extension during the Gothian orogeny and at the time of rapakivi magmatism and intrusion of diabase dikes; (4) NE-SW transtension that occurred between 1.60 and 1.30 Ga and which also formed the NW-SE-trending Satakunta graben located some 20 km north of Olkiluoto. Greisen-type veins also formed during this phase. (5) NE-SW compression that postdates both the formation of the 1.56 Ga rapakivi granites and 1.27 Ga olivine diabases of the region; (6) E-W transpression during the early stages of the Mesoproterozoic Sveconorwegian orogeny and which also predated (7) almost coaxial E-W extension attributed to the collapse of the Sveconorwegian orogeny. The kinematic analysis of fracture systems in crystalline bedrock also provides a robust framework for evaluating fluid-rock interaction in the brittle regime; this is essential in assessment of bedrock integrity for numerous geo-engineering applications, including groundwater management, transient or permanent CO2 storage and site investigations for permanent waste disposal. Investigations at Olkiluoto revealed that fluid flow along fractures is coupled with low normal tractions due to in-situ stresses and thus deviates from the generally accepted critically stressed fracture concept, where fluid flow is concentrated on fractures on the verge of failure. The difference is linked to the shallow conditions of Olkiluoto - due to the low differential stresses inherent at shallow depths, fracture activation and fluid flow is controlled by dilation due to low normal tractions. At deeper settings, however, fluid flow is controlled by fracture criticality caused by large differential stress, which drives shear deformation instead of dilation.

Réduction d'expansion volumique au cours de la gestation : conséquence sur la circulation utéro-placentaire.

La grossesse est caractérisée par une augmentation du volume circulant maternel et du débit sanguin afin de subvenir aux besoins croissants du foetus. En diminuant l'expansion volumique maternelle, nous avons développé un modèle de restriction de croissance intra-utérine (RCIU) chez la rate dans lequel on observe une diminution du diamètre des artères utérines arquées ainsi qu'une diminution de la réponse à l'angiotensine II. Nous avons émis l'hypothèse que comparativement aux rates normales, les vaisseaux utérins des rates RCIU présentent une dysfonction des mécanismes de relaxation dépendant de l'endothélium lors de la réponse à des agents vasoactifs. Notre objectif était de caractériser l'influence de l'endothélium sur la régulation de la réactivité des artères utérines radiales dans notre modèle RCIU. À l'aide d'un myographe pressurisé, des courbes de concentration-réponse à des agents vasoactifs (Phényléphrine (Pe), Carbachol (Cbc) et Nitroprussiate de sodium (SNP)) ont été effectuées, en présence ou en absence d'inhibiteurs des voies principales de la relaxation dépendante de l'endothélium (monoxyde d'azote (NO), prostacycline (PGI2) et facteurs hyperpolarisants dépendant de l'endothélium (EDHF)). Nous avons ainsi démontré que les artères utérines radiales des rates RCIU ont un diamètre plus petit que celles des rates contrôles. Une utilisation plus grande du NO dans la réponse dilatatrice au Cbc a été observée chez les rates RCIU. De plus, l'inhibition de la voie de l'EDHF par l'ajout de KCl aux inhibiteurs de la voie du NO (L-NAME) et des PGI2 (Ibuprofène) bloque presque complètement la relaxation au Cbc chez les deux groupes. En conclusion, dans les artères utérines radiales des deux groupes, il y a une participation importante de l'EDHF lors de la dilatation au Cbc. De plus, dans notre modèle RCIU, la voie du NO est activée, probablement afin de compenser la diminution du diamètre des artères et la réduction de la perfusion utérine.

A High Order Finite Volume Scheme for the 2D Shallow Water Equations Including Topography

Inhalt dieser Arbeit ist ein Verfahren zur numerischen Lösung der zweidimensionalen Flachwassergleichung, welche das Fließverhalten von Gewässern, deren Oberflächenausdehnung wesentlich größer als deren Tiefe ist, modelliert. Diese Gleichung beschreibt die gravitationsbedingte zeitliche Änderung eines gegebenen Anfangszustandes bei Gewässern mit freier Oberfläche. Diese Klasse beinhaltet Probleme wie das Verhalten von Wellen an flachen Stränden oder die Bewegung einer Flutwelle in einem Fluss. Diese Beispiele zeigen deutlich die Notwendigkeit, den Einfluss von Topographie sowie die Behandlung von Nass/Trockenübergängen im Verfahren zu berücksichtigen. In der vorliegenden Dissertation wird ein, in Gebieten mit hinreichender Wasserhöhe, hochgenaues Finite-Volumen-Verfahren zur numerischen Bestimmung des zeitlichen Verlaufs der Lösung der zweidimensionalen Flachwassergleichung aus gegebenen Anfangs- und Randbedingungen auf einem unstrukturierten Gitter vorgestellt, welches in der Lage ist, den Einfluss topographischer Quellterme auf die Strömung zu berücksichtigen, sowie in sogenannten \glqq lake at rest\grqq-stationären Zuständen diesen Einfluss mit den numerischen Flüssen exakt auszubalancieren. Basis des Verfahrens ist ein Finite-Volumen-Ansatz erster Ordnung, welcher durch eine WENO Rekonstruktion unter Verwendung der Methode der kleinsten Quadrate und eine sogenannte Space Time Expansion erweitert wird mit dem Ziel, ein Verfahren beliebig hoher Ordnung zu erhalten. Die im Verfahren auftretenden Riemannprobleme werden mit dem Riemannlöser von Chinnayya, LeRoux und Seguin von 1999 gelöst, welcher die Einflüsse der Topographie auf den Strömungsverlauf mit berücksichtigt. Es wird in der Arbeit bewiesen, dass die Koeffizienten der durch das WENO-Verfahren berechneten Rekonstruktionspolynome die räumlichen Ableitungen der zu rekonstruierenden Funktion mit einem zur Verfahrensordnung passenden Genauigkeitsgrad approximieren. Ebenso wird bewiesen, dass die Koeffizienten des aus der Space Time Expansion resultierenden Polynoms die räumlichen und zeitlichen Ableitungen der Lösung des Anfangswertproblems approximieren. Darüber hinaus wird die wohlbalanciertheit des Verfahrens für beliebig hohe numerische Ordnung bewiesen. Für die Behandlung von Nass/Trockenübergangen wird eine Methode zur Ordnungsreduktion abhängig von Wasserhöhe und Zellgröße vorgeschlagen. Dies ist notwendig, um in der Rechnung negative Werte für die Wasserhöhe, welche als Folge von Oszillationen des Raum-Zeit-Polynoms auftreten können, zu vermeiden. Numerische Ergebnisse die die theoretische Verfahrensordnung bestätigen werden ebenso präsentiert wie Beispiele, welche die hervorragenden Eigenschaften des Gesamtverfahrens in der Berechnung herausfordernder Probleme demonstrieren.

A high-order arbitrarily unstructured finite-volume model of the global atmosphere: tests solving the shallow-water equations

Simulations of the global atmosphere for weather and climate forecasting require fast and accurate solutions and so operational models use high-order finite differences on regular structured grids. This precludes the use of local refinement; techniques allowing local refinement are either expensive (eg. high-order finite element techniques) or have reduced accuracy at changes in resolution (eg. unstructured finite-volume with linear differencing). We present solutions of the shallow-water equations for westerly flow over a mid-latitude mountain from a finite-volume model written using OpenFOAM. A second/third-order accurate differencing scheme is applied on arbitrarily unstructured meshes made up of various shapes and refinement patterns. The results are as accurate as equivalent resolution spectral methods. Using lower order differencing reduces accuracy at a refinement pattern which allows errors from refinement of the mountain to accumulate and reduces the global accuracy over a 15 day simulation. We have therefore introduced a scheme which fits a 2D cubic polynomial approximately on a stencil around each cell. Using this scheme means that refinement of the mountain improves the accuracy after a 15 day simulation. This is a more severe test of local mesh refinement for global simulations than has been presented but a realistic test if these techniques are to be used operationally. These efficient, high-order schemes may make it possible for local mesh refinement to be used by weather and climate forecast models.

Evidence for global circuit current flow through water droplet layers

Abstract Foggy air and clear air have appreciably different electrical conductivities. The conductivity gradient at horizontal droplet boundaries causes droplet charging, as a result of vertical current flow in the global atmospheric electrical circuit. The charging is poorly known, as both the current flow through atmospheric water droplet layers and the air conductivity are poorly characterised experimentally. Surface measurements during three days of continuous fog using new instrument techniques show that a shallow (of order 100 m deep) fog layer still permits the vertical conduction current to pass. Further, the conductivity in the fog is estimated to be approximately 20% lower than in clear air. Assuming a fog transition thickness of one metre, this implies a vertical conductivity gradient of order 10 fS m−2 at the boundary. The actual vertical conductivity gradient at a cloud boundary would probably be greater, due to the presence of larger droplets in clouds compared to fog, and cleaner, more conductive clear air aloft.

Voronoi, Delaunay, and Block-Structured Mesh Refinement for Solution of the Shallow-Water Equations on the Sphere

Alternative meshes of the sphere and adaptive mesh refinement could be immensely beneficial for weather and climate forecasts, but it is not clear how mesh refinement should be achieved. A finite-volume model that solves the shallow-water equations on any mesh of the surface of the sphere is presented. The accuracy and cost effectiveness of four quasi-uniform meshes of the sphere are compared: a cubed sphere, reduced latitude–longitude, hexagonal–icosahedral, and triangular–icosahedral. On some standard shallow-water tests, the hexagonal–icosahedral mesh performs best and the reduced latitude–longitude mesh performs well only when the flow is aligned with the mesh. The inclusion of a refined mesh over a disc-shaped region is achieved using either gradual Delaunay, gradual Voronoi, or abrupt 2:1 block-structured refinement. These refined regions can actually degrade global accuracy, presumably because of changes in wave dispersion where the mesh is highly nonuniform. However, using gradual refinement to resolve a mountain in an otherwise coarse mesh can improve accuracy for the same cost. The model prognostic variables are height and momentum collocated at cell centers, and (to remove grid-scale oscillations of the A grid) the mass flux between cells is advanced from the old momentum using the momentum equation. Quadratic and upwind biased cubic differencing methods are used as explicit corrections to a fast implicit solution that uses linear differencing.

The acoustic signature of fluid flow in complex porous media

Effective medium approximations for the frequency-dependent and complex-valued effective stiffness tensors of cracked/ porous rocks with multiple solid constituents are developed on the basis of the T-matrix approach (based on integral equation methods for quasi-static composites), the elastic - viscoelastic correspondence principle, and a unified treatment of the local and global flow mechanisms, which is consistent with the principle of fluid mass conservation. The main advantage of using the T-matrix approach, rather than the first-order approach of Eshelby or the second-order approach of Hudson, is that it produces physically plausible results even when the volume concentrations of inclusions or cavities are no longer small. The new formulae, which operates with an arbitrary homogeneous (anisotropic) reference medium and contains terms of all order in the volume concentrations of solid particles and communicating cavities, take explicitly account of inclusion shape and spatial distribution independently. We show analytically that an expansion of the T-matrix formulae to first order in the volume concentration of cavities (in agreement with the dilute estimate of Eshelby) has the correct dependence on the properties of the saturating fluid, in the sense that it is consistent with the Brown-Korringa relation, when the frequency is sufficiently low. We present numerical results for the (anisotropic) effective viscoelastic properties of a cracked permeable medium with finite storage porosity, indicating that the complete T-matrix formulae (including the higher-order terms) are generally consistent with the Brown-Korringa relation, at least if we assume the spatial distribution of cavities to be the same for all cavity pairs. We have found an efficient way to treat statistical correlations in the shapes and orientations of the communicating cavities, and also obtained a reasonable match between theoretical predictions (based on a dual porosity model for quartz-clay mixtures, involving relatively flat clay-related pores and more rounded quartz-related pores) and laboratory results for the ultrasonic velocity and attenuation spectra of a suite of typical reservoir rocks. (C) 2003 Elsevier B.V. All rights reserved.