Period adding cascades: Experiment and modeling in air bubbling

Period adding cascades have been observed experimentally/numerically in the dynamics of neurons and pancreatic cells, lasers, electric circuits, chemical reactions, oceanic internal waves, and also in air bubbling. We show that the period adding cascades appearing in bubbling from a nozzle submerged in a viscous liquid can be reproduced by a simple model, based on some hydrodynamical principles, dealing with the time evolution of two variables, bubble position and pressure of the air chamber, through a system of differential equations with a rule of detachment based on force balance. The model further reduces to an iterating one-dimensional map giving the pressures at the detachments, where time between bubbles come out as an observable of the dynamics. The model has not only good agreement with experimental data, but is also able to predict the influence of the main parameters involved, like the length of the hose connecting the air supplier with the needle, the needle radius and the needle length. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.3695345]

Use of Chlorine to Remove Magnesium from Molten Aluminum

Removal of Mg from aluminum scraps, known as demagging, has been widely applied in the,aluminum industry. This work discusses bubble-formation theories and magnesium kinetic removal from aluminum scraps using chlorine and inert gas fluxing. The interfacial area of the bubbles and residence time were estimated using a mathematical model. To inject gaseous chlorine, three types of nozzles were used with varying internal diameter. In addition, a porous plug, as well as varying input chlorine flow and concentration were used. The use of lower chlorine concentration improves efficiency because the interfacial tension is reduced therefore, more and smaller bubbles are formed. The model proposed herein is consistent with the experimental data. [doi:10.2320/matertrans.M2011256]

Parâmetros de dispersão de ar em células mecânicas de flotação em usinas brasileiras

O objetivo desse trabalho é apresentar e discutir a determinação dos parâmetros de dispersão de ar em células mecânicas de flotação, em circuitos industriais brasileiros. Foram investigados os circuitos rougher e scavenger das usinas RPM, Taquari-Vassouras, Copebras e Fosfértil-Catalão (atual Vale Fertilizantes). Nas usinas RPM e Taquari-Vassouras, foram determinados a velocidade superficial do ar (Jg), o hold-up do ar (εg), o diâmetro médio de bolhas (d32) e o fluxo de área superficial de bolhas (Sb). Para as outras usinas, foram determinados somente o Jg e εg. Os resultados mostraram que, com exceção da RPM, todos os circuitos investigados apresentaram baixos valores para os parâmetros de dispersão do ar. Taquari-Vassouras apresentou os mais baixos valores médios para Jg e εg, respectivamente, 0,05 cm/s e 5%. Os maiores valores de Jg foram determinados no circuito rougher da Copebras, 1,17 cm/s, ao passo que a RPM apresentou os maiores valores para o hold-up e tamanho de bolha, respectivamente, 31% e 2,61mm. Os valores obtidos, para Sb, situaram-se abaixo daqueles reportados na literatura.

THE INTERSTELLAR MAGNETIC FIELD CLOSE TO THE SUN. II.

The magnetic field in the local interstellar medium (ISM) provides a key indicator of the galactic environment of the Sun and influences the shape of the heliosphere. We have studied the interstellar magnetic field (ISMF) in the solar vicinity using polarized starlight for stars within 40 pc of the Sun and 90 degrees of the heliosphere nose. In Frisch et al. (Paper I), we developed a method for determining the local ISMF direction by finding the best match to a group of interstellar polarization position angles obtained toward nearby stars, based on the assumption that the polarization is parallel to the ISMF. In this paper, we extend the analysis by utilizing weighted fits to the position angles and by including new observations acquired for this study. We find that the local ISMF is pointed toward the galactic coordinates l, b = 47 degrees +/- 20 degrees, 25 degrees +/- 20 degrees. This direction is close to the direction of the ISMF that shapes the heliosphere, l, b = 33 degrees +/- 4 degrees, 55 degrees +/- 4 degrees, as traced by the center of the "Ribbon" of energetic neutral atoms discovered by the Interstellar Boundary Explorer (IBEX) mission. Both the magnetic field direction and the kinematics of the local ISM are consistent with a scenario where the local ISM is a fragment of the Loop I superbubble. A nearby ordered component of the local ISMF has been identified in the region l approximate to 0 degrees -> 80 degrees and b approximate to 0 degrees -> 30 degrees, where PlanetPol data show a distance-dependent increase of polarization strength. The ordered component extends to within 8 pc of the Sun and implies a weak curvature in the nearby ISMF of +/- 0 degrees.25 pc(-1). This conclusion is conditioned on the small sample of stars available for defining this rotation. Variations from the ordered component suggest a turbulent component of +/- 23 degrees. The ordered component and standard relations between polarization, color excess, and H-o column density predict a reasonable increase of N(H) with distance in the local ISM. The similarity of the ISMF directions traced by the polarizations, the IBEX Ribbon, and pulsars inside the Local Bubble in the third galactic quadrant suggest that the ISMF is relatively uniform over spatial scales of 8-200 pc and is more similar to interarm than spiral-arm magnetic fields. The ISMF direction from the polarization data is also consistent with small-scale spatial asymmetries detected in GeV-TeV cosmic rays with a galactic origin. The peculiar geometrical relation found earlier between the cosmic microwave background dipole moment, the heliosphere nose, and the ISMF direction is supported by this study. The interstellar radiation field at +/- 975 angstrom does not appear to play a role in grain alignment for the low-density ISM studied here.

Experimental investigation of gravitational gas separation in an inclined annular channel

The oil industry uses gas separators in production wells as the free gas present in the suction of the pump reduces the pumping efficiency and pump lifetime. Therefore, free gas is one of the most important variables in the design of pumping systems. However, in the literature there is little information on these separators. It is the case of the inverted-shroud gravitational gas separator. It has an annular geometry due to the installation of a cylindrical container in between the well casing and pioduction pipe (tubing). The purpose of the present study is to understand the phenomenology and behavior of inverted-shroud separator. Experimental tests were performed in a 10.5-m-length inclinable glass tube with air and water as working fluids. The water flow rate was in the range of 8.265-26.117 l/min and the average inlet air mass flow rate was 1.1041 kg/h, with inclination angles of 15 degrees, 30 degrees, 45 degrees, 60 degrees, 75 degrees, 80 degrees and 85 degrees. One of the findings is that the length between the inner annular level and production pipe inlet is one of the most important design parameters and based on that a new criterion for total gas separation is proposed. We also found that the phenomenology of the studied separator is not directly dependent on the gas flow rate, but on the average velocity of the free surface flow generated inside the separator. Maps of efficiency of gas separation were plotted and showed that liquid flow rate, inclination angle and pressure difference between casing and production pipe outlet are the main variables related to the gas separation phenomenon. The new data can be used for the development of design tools aiming to the optimized project of the pumping system for oil production in directional wells. (C) 2012 Elsevier Inc. All rights reserved.

Ajuste e balanceamento do portfólio de projetos: o caso de uma empresa do setor químico

O objetivo deste trabalho é compreender a etapa de ajuste no contexto da gestão do portfólio de projetos, destacando sua relação com os processos de categorização e balanceamento. A pesquisa realizada tem caráter qualitativo, sendo a abordagem adotada o estudo de caso longitudinal. A pesquisa foi desenvolvida em uma empresa do setor químico brasileiro. As evidências, de várias fontes, foram coletadas através de entrevistas, documentos e dados dos sistemas corporativos. Para a compreensão do portfólio de projetos da empresa foram coletados e analisados dados de mil projetos realizados entre 2001 e 2005. Os resultados indicam que maior atenção é dada à etapa de seleção, negligenciando a etapa de ajuste. A adoção de ferramentas de balanceamento permitiu evidenciar lacunas e fontes de desbalanceamento no portfólio de projetos, promovendo o debate entre os tomadores de decisão no que concerne ao viés introduzido pelos critérios adotados na etapa de seleção e levantando a necessidade de introdução de uma sistemática de ajuste e balanceamento. Observou-se que sem uma adequada categorização dos projetos da empresa seria difícil promover a análise de balanceamento.

Usable work of macro-scale cavities in liquids

It is shown that the generation of cavities in a liquid can produce usable work, which is illustrated by the stretching of a string. This work is done during the expansion of the cavity, and not with its collapse. Basic equations are presented for the movement of a device moved by the so called cavity events. A theoretical solution is also proposed, which uses polynomial functions relating the so called "excess of pressure" in the cavity and time. Evaluations of the force generated during the expansion of the cavity showed a mean peak value of about 58 N for the moving container, while measurements with the container fixed to a support showed a peak value of 476 N, considered somewhat overestimated, because high frequency oscillations seem to superpose the mean behavior. Simultaneous phenomena occurring during the cavity events are also described. Series of pictures of the experiments are presented.

Two fluid numerical model for coastal applications

Wave breaking is an important coastal process, influencing hydro-morphodynamic processes such as turbulence generation and wave energy dissipation, run-up on the beach and overtopping of coastal defence structures. During breaking, waves are complex mixtures of air and water (“white water”) whose properties affect velocity and pressure fields in the vicinity of the free surface and, depending on the breaker characteristics, different mechanisms for air entrainment are usually observed. Several laboratory experiments have been performed to investigate the role of air bubbles in the wave breaking process (Chanson & Cummings, 1994, among others) and in wave loading on vertical wall (Oumeraci et al., 2001; Peregrine et al., 2006, among others), showing that the air phase is not negligible since the turbulent energy dissipation involves air-water mixture. The recent advancement of numerical models has given valuable insights in the knowledge of wave transformation and interaction with coastal structures. Among these models, some solve the RANS equations coupled with a free-surface tracking algorithm and describe velocity, pressure, turbulence and vorticity fields (Lara et al. 2006 a-b, Clementi et al., 2007). The single-phase numerical model, in which the constitutive equations are solved only for the liquid phase, neglects effects induced by air movement and trapped air bubbles in water. Numerical approximations at the free surface may induce errors in predicting breaking point and wave height and moreover, entrapped air bubbles and water splash in air are not properly represented. The aim of the present thesis is to develop a new two-phase model called COBRAS2 (stands for Cornell Breaking waves And Structures 2 phases), that is the enhancement of the single-phase code COBRAS0, originally developed at Cornell University (Lin & Liu, 1998). In the first part of the work, both fluids are considered as incompressible, while the second part will treat air compressibility modelling. The mathematical formulation and the numerical resolution of the governing equations of COBRAS2 are derived and some model-experiment comparisons are shown. In particular, validation tests are performed in order to prove model stability and accuracy. The simulation of the rising of a large air bubble in an otherwise quiescent water pool reveals the model capability to reproduce the process physics in a realistic way. Analytical solutions for stationary and internal waves are compared with corresponding numerical results, in order to test processes involving wide range of density difference. Waves induced by dam-break in different scenarios (on dry and wet beds, as well as on a ramp) are studied, focusing on the role of air as the medium in which the water wave propagates and on the numerical representation of bubble dynamics. Simulations of solitary and regular waves, characterized by both spilling and plunging breakers, are analyzed with comparisons with experimental data and other numerical model in order to investigate air influence on wave breaking mechanisms and underline model capability and accuracy. Finally, modelling of air compressibility is included in the new developed model and is validated, revealing an accurate reproduction of processes. Some preliminary tests on wave impact on vertical walls are performed: since air flow modelling allows to have a more realistic reproduction of breaking wave propagation, the dependence of wave breaker shapes and aeration characteristics on impact pressure values is studied and, on the basis of a qualitative comparison with experimental observations, the numerical simulations achieve good results.

Correlazioni tra proprietá reologiche, struttura e processabilitá di blend di LLDPE/LDPE

Polymer blends constitute a valuable way to produce relatively low cost new materials. A still open question concerns the miscibility of polyethylene blends. Deviations from the log-additivity rule of the newtonian viscosity are often taken as a signature of immiscibility of the two components. The aim of this thesis is to characterize the rheological behavior in shear and elongation of five series of LLDPE/LDPE blends whose parent polymers have been chosen with different viscosity and SCB content and length. Synergistic effects have been measured for both zero shear viscosity and melt strength. Both SCB length and viscosity ratio between the components have been found to be key parameters for the miscibility of the pure polymers. In particular the miscibility increases with increasing SCB length and with decreasing the LDPE molecular weight and viscosity. This rheological behavior has significant effects on the processability window of these blends when the uni or biaxial elongational flows are involved. The film blowing is one of the processes for which the synergistic effects above mentioned can be crucial. Small scale experiments of film blowing performed for one of the series of blends has demonstrated that the positive deviation of the melt strength enlarges the processability window. In particular, the bubble stability was found to improve or disappear when the melt strength of the samples increased. The blending of LDPE and LLDPE can even reduce undesired melt flow instability phenomena widening, as a consequence, the processability window in extrusion. One of the series of blends has been characterized by means of capillary rheometry in order to allow a careful morphological analysis of the surface of the extruded polymer jets by means of Scanning Electron Microscopy (SEM) with the aim to detect the very early stages of the small scale melt instabilty at low shear rates (sharksin) and to follow its subsequent evolution as long as the shear rate was increased. With this experimental procedure it was possible to evaluate the shear rate ranges corresponding to different flow regions: smooth extrudate surface (absence of instability), sharkskin (small scale instability produced at the capillary exit), stick-slip transition (instability involving the whole capillary wall) and gross melt fracture (i.e. a large scale "upstream" instability originating from the entrance region of the capillary). A quantitative map was finally worked out using which an assessment of the flow type for a given shear rate and blend composition can be predicted.

Astrophysikalische Bedingungen für einen r-Prozess im Hoch-Entropie-Wind von Typ II Supernovae

Sterne mit einer Anfangsmasse zwischen etwa 8 und 25 Sonnenmassen enden ihre Existenz mit einer gewaltigen Explosion, einer Typ II Supernova. Die hierbei entstehende Hoch-Entropie-Blase ist ein Bereich am Rande des sich bildenden Neutronensterns und gilt als möglicher Ort für den r-Prozess. Wegen der hohen Temperatur T innerhalb der Blase ist die Materie dort vollkommen photodesintegriert. Das Verhältnis von Neutronen zu Protonen wird durch die Elektronenhäufigkeit Ye beschrieben. Die thermodynamische Entwicklung des Systems wird durch die Entropie S gegeben. Da die Expansion der Blase schnell vonstatten geht, kann sie als adiabatisch betrachtet werden. Die Entropie S ist dann proportional zu T^3/rho, wobei rho die Dichte darstellt. Die explizite Zeitentwicklung von T und rho sowie die Prozessdauer hängen von Vexp, der Expansionsgeschwindigkeit der Blase, ab. Der erste Teil dieser Dissertation beschäftigt sich mit dem Prozess der Reaktionen mit geladenen Teilchen, dem alpha-Prozess. Dieser Prozess endet bei Temperaturen von etwa 3 mal 10^9 K, dem sogenannten "alpha-reichen" Freezeout, wobei überwiegend alpha-Teilchen, freie Neutronen sowie ein kleiner Anteil von mittelschweren "Saat"-Kernen im Massenbereich um A=100 gebildet werden. Das Verhältnis von freien Neutronen zu Saatkernen Yn/Yseed ist entscheidend für den möglichen Ablauf eines r-Prozesses. Der zweite Teil dieser Arbeit beschäftigt sich mit dem eigentlichen r-Prozess, der bei Neutronenanzahldichten von bis zu 10^27 Neutronen pro cm^3 stattfindet, und innerhalb von maximal 400 ms sehr neutronenreiche "Progenitor"-Isotope von Elementen bis zum Thorium und Uran bildet. Bei dem sich anschliessendem Ausfrieren der Neutroneneinfangreaktionen bei 10^9 K und 10^20 Neutronen pro cm^3 erfolgt dann der beta-Rückzerfall der ursprünglichen r-Prozesskerne zum Tal der Stabilität. Diese Nicht-Gleichgewichts-Phase wird in der vorliegenden Arbeit in einer Parameterstudie eingehend untersucht. Abschliessend werden astrophysikalische Bedingungen definiert, unter denen die gesamte Verteilung der solaren r-Prozess-Isotopenhäufigkeiten reproduziert werden können.

CFD modeling of two-phase boiling flows in the slug flow regime with an interface capturing technique

The objective of this thesis was to improve the commercial CFD software Ansys Fluent to obtain a tool able to perform accurate simulations of flow boiling in the slug flow regime. The achievement of a reliable numerical framework allows a better understanding of the bubble and flow dynamics induced by the evaporation and makes possible the prediction of the wall heat transfer trends. In order to save computational time, the flow is modeled with an axisymmetrical formulation. Vapor and liquid phases are treated as incompressible and in laminar flow. By means of a single fluid approach, the flow equations are written as for a single phase flow, but discontinuities at the interface and interfacial effects need to be accounted for and discretized properly. Ansys Fluent provides a Volume Of Fluid technique to advect the interface and to map the discontinuous fluid properties throughout the flow domain. The interfacial effects are dominant in the boiling slug flow and the accuracy of their estimation is fundamental for the reliability of the solver. Self-implemented functions, developed ad-hoc, are introduced within the numerical code to compute the surface tension force and the rates of mass and energy exchange at the interface related to the evaporation. Several validation benchmarks assess the better performances of the improved software. Various adiabatic configurations are simulated in order to test the capability of the numerical framework in modeling actual flows and the comparison with experimental results is very positive. The simulation of a single evaporating bubble underlines the dominant effect on the global heat transfer rate of the local transient heat convection in the liquid after the bubble transit. The simulation of multiple evaporating bubbles flowing in sequence shows that their mutual influence can strongly enhance the heat transfer coefficient, up to twice the single phase flow value.

Identifizierung und Charakterisierung von Intermediärfilament-Organisatoren in Caenorhabditis elegans

Intermediärfilamente (IFs) sind neben Mikrotubuli und Aktinfilamenten die dritte filamentäre Komponente des Zytoskeletts. Sie wirken als mechanische Stabilisatoren, sind außerdem an Zelldifferenzierung, Proliferation und Apoptose beteiligt und tragen zu Zellpolarität bei. IFs sind dynamische Strukturen, die zelltypspezifisch in unterschiedlichen Anordnungen und Abundanzen vorkommen und von Signalkaskaden beeinflusst werden. Die zugrundeliegenden molekularen Mechanismen dieser fein abgestimmten Prozesse sind weitgehend unbekannt. In dieser Arbeit sollte deswegen ein Tiermodell entwickelt werden, um Regulatoren der IF-(Netzwerk)-Organisation in vivo zu untersuchen und zu identifizieren. Dazu wurde C. elegans ausgewählt, da es sich hierbei um einen genetisch gut charakterisierten und leicht manipulierbaren Organismus handelt, in dessen Genom elf Gene für zytoplasmatische IFs kodieren. Zunächst wurden stabil transgene C. elegans-Linien generiert, die fluoreszierende IFs exprimieren. Es konnte gezeigt werden, dass das darmspezifische IFB-2::CFP im Bereich des apikalen Junktionskomplex verankert ist und nahezu vollständig im subapikalen Terminalgeflecht der Enterozyten lokalisiert, das als Teil der endotube besonders stabil und widerstandsfähig ist. Wenn diese Tiere mit dsRNA gegen das ebenfalls im Terminalgeflecht exprimierte IF ifc-2 behandelt wurden, entwickelten sich blasenförmige Ausstülpungen des Darmlumens, die auf eine Schwächung der rigiden und formgebenden endotube hinwiesen und damit einen direkten in vivo-Beweis für die stressprotektive Funktion des intestinalen IF-Netzwerks lieferten. Die leichte Detektierbarkeit des IFB-2::CFP-Musters wurde in einem optischen Screen ausgenutzt, bei dem nach chemischer Mutagenese nach Veränderungen im IF-Muster gefahndet wurde. Hierbei wurden drei Mutanten isoliert. In Komplementationsanalysen stellte sich heraus, dass es sich in zwei Fällen um Allele desselben Gens handelt. Die Identifizierung der betroffenen Gene gelang durch eine PCR-basierte Kartierung von single nucleotide polymorphisms nach Verpaarung mit dem Hawaii-Stamm (snp-mapping) und anschließender RNAi-Analyse der Einzelgene in den identifizierten Chromosomenabschnitten. Im einen Fall handelte es sich um das sma-5-Gen, einer Serin/Threonin-Kinase mit Homologie zu den MAP-Kinasen MAPK7/ERK5 der Säuger. Hier wurden, ebenso wie beim ifc-2 (RNAi)-Phänotyp, progressive blasenförmige Ausstülpungen des Darmlumens beobachtet. Die beiden anderen Allele tragen Mutationen in einem bisher nicht näher charakterisierten Gen. In diesen Würmern kommt es zu einem vollständigen Auflösung des IFB-2::CFP-Netzwerks mit prominenten Akkumulationen um die apikalen Junktionen. Das Darmlumen ist stellenweise geweitet und das elektronendichte Terminalgeflecht fehlt fast vollständig, die Integrität des Darmepithels ist jedoch nicht kompromittiert. Die anderen IFs des Terminalgeflechts sind ebenfalls fehlverteilt, und die intestinale Expression von Aktin ist stark reduziert. Expressionskonstrukte des Gens zeigten weiterhin, dass es darmspezifisch synthetisiert wird und mit den IFs im Terminalgeflecht kolokalisiert. Das Protein ist, ähnlich wie das IF-assoziierte Filaggrin der Säuger ausgesprochen histidinreich. Es enthält außerdem eine Prolin-reiche Domäne, die Teil einer potentiellen Aktin-Bindedomäne ist. Auf Grund all dieser Eigenschaften wird die Bezeichnung IFO-1 (intermediate filament organizer) für das neue Protein vorgeschlagen, das möglicherweise als struktureller Zytoskelett-Linker wirkt. Die vorgestellten Ergebnisse untermauern die Bedeutung von C. elegans für die Identifizierung von Faktoren, die IF-Netzwerke regulieren, und die Möglichkeit, Defekte im lebenden Gesamtorganismus zu bestimmen.

Sintesi, modifica e caratterizzazione di polimeri da fonti rinnovabili

Gli argomenti trattati all’interno di questa tesi di dottorato riguardano la sintesi e la modifica di polimeri ottenuti a partire da fonti rinnovabili. L’acido polilattico (PLA) è stato modificato per ottenere film estensibili per uso alimentare. La scelta del materiale si è basata sull’analisi del suo ciclo di vita e perché è riconosciuto come sicuro per l’utilizzo nel campo alimentare. Le formulazioni testate, a base di PLA, sono state preparate con l’aggiunta di una serie di additivi utilizzati per migliorare le proprietà meccaniche del materiale. La lavorazione è stata eseguita mediante estrusione, ottenendo dei pellet a composizione omogenea successivamente lavorati nell’estrusore a bolla, modalità industriale di produzione dei film estensibili. È stata poi effettuata la sintesi diretta di nuovi poliesteri insaturi a base di dimetil succinato e 1,6-esandiolo. L’insaturazione della catena è stata ottenuta mediante l’uso, durante la sintesi, di derivati dell’olio di ricino, l’acido ricinoleico e il suo derivato insaturo metil undecenoato. Un’altra molecola insatura utilizzata è stata il citronellolo, scelto tra i terpeni. I polimeri così ottenuti sono stati modificati tramite reazioni radicaliche indotte con radiazioni UV utilizzando sostanze antibatteriche (sale 3-esadecil-1-vinilimidazolo bromuro) al fine di ottenere materiali con attività biocida a lungo termine e senza rilascio. Si è proceduto inoltre alla polimerizzazione reversibile di monomeri furanici con oli vegetali utilizzando una strategia di tipo double click. Si è partiti dalla sintesi di monomeri derivanti da oli vegetali contenenti eterocicli furanici attaccati mediante addizione tiol-enica (prima reazione click chemistry) e si è proseguito con la loro successiva polimerizzazione attraverso una reazione di tipo Diels-Alder con molecole con gruppi maleimmidici (seconda reazione click chemistry). I polimeri così ottenuti sono materiali potenzialmente auto-riparanti, grazie alla possibilità di spostare l’equilibrio verso i prodotti o i reagenti semplicemente variando le condizioni di temperatura.

Numerical and Analytical Methods for Laser-Plasma Acceleration Physics

Theories and numerical modeling are fundamental tools for understanding, optimizing and designing present and future laser-plasma accelerators (LPAs). Laser evolution and plasma wave excitation in a LPA driven by a weakly relativistically intense, short-pulse laser propagating in a preformed parabolic plasma channel, is studied analytically in 3D including the effects of pulse steepening and energy depletion. At higher laser intensities, the process of electron self-injection in the nonlinear bubble wake regime is studied by means of fully self-consistent Particle-in-Cell simulations. Considering a non-evolving laser driver propagating with a prescribed velocity, the geometrical properties of the non-evolving bubble wake are studied. For a range of parameters of interest for laser plasma acceleration, The dependence of the threshold for self-injection in the non-evolving wake on laser intensity and wake velocity is characterized. Due to the nonlinear and complex nature of the Physics involved, computationally challenging numerical simulations are required to model laser-plasma accelerators operating at relativistic laser intensities. The numerical and computational optimizations, that combined in the codes INF&RNO and INF&RNO/quasi-static give the possibility to accurately model multi-GeV laser wakefield acceleration stages with present supercomputing architectures, are discussed. The PIC code jasmine, capable of efficiently running laser-plasma simulations on Graphics Processing Units (GPUs) clusters, is presented. GPUs deliver exceptional performance to PIC codes, but the core algorithms had to be redesigned for satisfying the constraints imposed by the intrinsic parallelism of the architecture. The simulation campaigns, run with the code jasmine for modeling the recent LPA experiments with the INFN-FLAME and CNR-ILIL laser systems, are also presented.

Untersuchungen zur Genauigkeit adaptiver unstetiger Galerkin-Simulationen mit Hilfe von Luftblasen-Testfällen

Allgemein erlaubt adaptive Gitterverfeinerung eine Steigerung der Effizienz numerischer Simulationen ohne dabei die Genauigkeit des Ergebnisses signifikant zu verschlechtern. Es ist jedoch noch nicht erforscht, in welchen Bereichen des Rechengebietes die räumliche Auflösung tatsächlich vergröbert werden kann, ohne die Genauigkeit des Ergebnisses signifikant zu beeinflussen. Diese Frage wird hier für ein konkretes Beispiel von trockener atmosphärischer Konvektion untersucht, nämlich der Simulation von warmen Luftblasen. Zu diesem Zweck wird ein neuartiges numerisches Modell entwickelt, das auf diese spezielle Anwendung ausgerichtet ist. Die kompressiblen Euler-Gleichungen werden mit einer unstetigen Galerkin Methode gelöst. Die Zeitintegration geschieht mit einer semi-implizite Methode und die dynamische Adaptivität verwendet raumfüllende Kurven mit Hilfe der Funktionsbibliothek AMATOS. Das numerische Modell wird validiert mit Hilfe einer Konvergenzstudie und fünf Standard-Testfällen. Eine Methode zum Vergleich der Genauigkeit von Simulationen mit verschiedenen Verfeinerungsgebieten wird eingeführt, die ohne das Vorhandensein einer exakten Lösung auskommt. Im Wesentlichen geschieht dies durch den Vergleich von Eigenschaften der Lösung, die stark von der verwendeten räumlichen Auflösung abhängen. Im Fall einer aufsteigenden Warmluftblase ist der zusätzliche numerische Fehler durch die Verwendung der Adaptivität kleiner als 1% des gesamten numerischen Fehlers, wenn die adaptive Simulation mehr als 50% der Elemente einer uniformen hoch-aufgelösten Simulation verwendet. Entsprechend ist die adaptive Simulation fast doppelt so schnell wie die uniforme Simulation.