Understanding Optical Trapping Phenomena: a Simulation for Undergraduates

Optical trapping is an attractive and multidisciplinary topic that has become the center of attention to a large number of researchers. Moreover, it is a suitable subject for advanced students that requires a knowledge of a wide range of topics. As a result, it has been incorporated into some syllabuses of both undergraduate and graduate programs. In this paper, basic concepts in laser trapping theory are reviewed. To provide a better understanding of the underlying concepts for students, a Java application for simulating the behavior of a dielectric particle trapped in a highly focused beam has been developed. The program illustrates a wide range of theoretical results and features, such as the calculation of the force exerted by a beam in the Mie and Rayleigh regimes or the calibration of the trap stiffness. Some examples that are ready to be used in the classroom or in the computer lab are also supplied.

Optical limiting and thermal lensing studies in C60

Optical limiting and thermo-optic properties of C60 in toluene are studied using 532 nm, 9 ns pulses from a frequency-doubled Nd:YAG laser. Optical limiting studies in these fullerene molecules lead to the conclusion that reverse saturable absorption is the major mechanism for limiting properties in these molecules. Thermal lensing measurements are also performed in fullerene solutions. The quadratic dependence of thermal lens signal on incident energy confirms that enhanced optical absorption by the sample via excited triplet state absorption may play a leading role in the limiting property.

In situ fast ellipsometric analysis of repetitive surface phenomena

We present an ellipsometric technique and ellipsometric analysis of repetitive phenomena, based on the experimental arrangement of conventional phase modulated ellipsometers (PME) c onceived to study fast surface phenomena in repetitive processes such as periodic and triggered experiments. Phase modulated ellipsometry is a highly sensitive surface characterization technique that is widely used in the real-time study of several processes such as thin film deposition and etching. However, fast transient phenomena cannot be analyzed with this technique because precision requirements limit the data acquisition rate to about 25 Hz. The presented new ellipsometric method allows the study of fast transient phenomena in repetitive processes with a time resolution that is mainly limited by the data acquisition system. As an example, we apply this new method to the study of surface changes during plasma enhanced chemical vapor deposition of amorphous silicon in a modulated radio frequency discharge of SiH4. This study has revealed the evolution of the optical parameters of the film on the millisecond scale during the plasma on and off periods. The presented ellipsometric method extends the capabilities of PME arrangements and permits the analysis of fast surface phenomena that conventional PME cannot achieve.

Simple holographic polarogram

A new technique capable of obtaining quantitative values of the rotation angle of the polarization vector by using holography is presented. This is a two-stage holographic process; during the recording stage a hologram of the object of interest is obtained. The reference beam is composed of two beams that form a small angle between them and keep their polarization states at right angles to each other. In the reconstruction stage of the hologram, two images from the hologram are obtained along two different angles. As a result of the interference between these two images, a set of parallel fringes is formed at the image plane. The fringe contrast on the reconstruction is related to the angle of the polarization vector of the light at each position on the image plane. Measurements of the rotation of the polarization angle of a fraction of a degree were obtained. The main application of this technique is in the study of transient phenomena, where single-shot measurements are the only means of obtaining reliable data.

SN 2011A: A low-luminosity interacting transient with a double plateau and strong sodium absorpton

We present optical photometry and spectroscopy of the optical transient SN 2011A. Our data span 140 days after discovery including BVRI u′g′r′i′z′ photometry and 11 epochs of optical spectroscopy. Originally classified as a type IIn supernova (SN IIn) due to the presence of narrow Hα emission, this object shows exceptional characteristics. First, the light curve shows a double plateau, a property only observed before in the impostor SN 1997bs. Second, SN 2011A has a very low luminosity (M_V=-15.72), placing it between normal luminous SNe IIn and SN impostors. Third, SN 2011A shows low velocity and high equivalent width absorption close to the sodium doublet, which increases with time and is most likely of circumstellar origin. This evolution is also accompanied by a change in line profile; when the absorption becomes stronger, a P Cygni profile appears. We discuss SN 2011A in the context of interacting SNe IIn and SN impostors, which appears to confirm the uniqueness of this transient. While we favor an impostor origin for SN 2011A, we highlight the difficulty in differentiating between terminal and non-terminal interacting transients.

Alterações da Actividade Autonómica Durante o Teste de Inclinação em Doentes com Fibrilhação Auricular Paroxística: Análise com Wavelets

The autonomic nervous system (ANS) is known to be an important modulator in the pathogenesis of paroxysmal atrial fibrillation (PAF). Changes in ANS control of heart rate variability (HRV) occur during orthostatism to maintain cardiovascular homeostasis. Wavelet transform has emerged as a useful tool that provides time-frequency decomposition of the signal under investigation, enabling intermittent components of transient phenomena to be analyzed. AIM: To study HRV during head-up tilt (HUT) with wavelet transform analysis in PAF patients and healthy individuals (normals). METHODS: Twenty-one patients with PAF (8 men; age 58 +/- 14 yrs) were examined and compared with 21 normals (7 men, age 48 +/- 12 yrs). After a supine resting period, all subjects underwent passive HUT (60 degrees) while in sinus rhythm. Continuous monitoring of ECG and blood pressure was carried out (Task Force Monitor, CNSystems). Acute changes in RR-intervals were assessed by wavelet analysis and low-frequency power (LF: 0.04-0.15 Hz), high-frequency power (HF: 0.15-0.60 Hz) and LF/HF (sympathovagal) were calculated for 1) the last 2 min of the supine period; 2) the 15 sec of tilting movement (TM); and 3) the 1st (TT1) and 2nd (TT2) min of HUT. Data are expressed as means +/- SEM. RESULTS: Baseline and HUT RR-intervals were similar for the two groups. Supine basal blood pressure was also similar for the two groups, with a sustained increase in PAF patients, and a decrease followed by an increase and then recovery in normals. Basal LF, HF and LF/ HF values in PAF patients were 632 +/- 162 ms2, 534 +/- 231 ms2 and 1.95 +/- 0.39 respectively, and 1058 +/- 223 ms2, 789 +/- 244 ms2 and 2.4 +/- 0.36 respectively in normals (p = NS). During TM, LF, HF and LF/HF values for PAF patients were 747 +/- 277 ms2, 387 +/- 94 ms2 and 2.9 +/- 0.6 respectively, and 1316 +/- 315 ms2, 698 +/- 148 ms2 and 2.8 +/- 0.6 respectively in normals (p < 0.05 for LF and HF). During TF1, LF, HF and LF/ HF values for PAF patients were 1243 +/- 432 ms2, 302 +/- 88 ms2 and 7.7 +/- 2.4 respectively, and 1992 +/- 398 ms2, 333 +/- 76 ms2 and 7.8 +/- 0.98 respectively for normals (p < 0.05 for LF). During TF2, LF, HF and LF/HF values for PAF patients were 871 +/- 256 ms2, 242 +/- 51 ms2 and 4.7 +/- 0.9 respectively, and 1263 +/- 335 ms2, 317 +/- 108 ms2 and 8.6 +/- 0.68 respectively for normals (p < 0.05 for LF/HF). The dynamic profile of HRV showed that LF and HF values in PAF patients did not change significantly during TM or TT2, and LF/HF did not change during TM but increased in TT1 and TT2. CONCLUSION: Patients with PAF present alterations in HRV during orthostatism, with decreased LF and HF power during TM, without significant variations during the first minutes of HUT. These findings suggest that wavelet transform analysis may provide new insights when assessing autonomic heart regulation and highlight the presence of ANS disturbances in PAF.

Effects of modern variable frequency drive on the reliability of an old crane induction motor

The work aims to analyze the possibilities of utilizing old crane driving AC induction motors in modern pulse-width-modulated variable frequency drives. Bearing currents and voltage stresses are the two main problems associated with modern IGBT inverters, and they may cause premature failure of an old induction motor. The origins of these two problems are studied. An analysis of the mechanism of bearing failure is proposed. Certain types of bearing currents are considered in detail. The most effective and economical means are chosen for bearing currents mitigation. Transient phenomena of cables and mechanism of over voltages occurring at motor terminals are studied in the work. The weakest places of the stator winding insulation system are shown and recommendations are given considering the mitigation of voltage stresses. Only the most appropriate and cost effective preventative methods are chosen for old motor drives. Rewinding of old motors is also considered.

Auroral and plasma flow transients at magnetic noon

We present observations of a transient event in the dayside auroral ionosphere at magnetic noon. F-region plasma convection measurements were made by the EISCAT radar, operating in the beamswinging “Polar” experiment mode, and simultaneous observations of the dayside auroral emissions were made by optical meridian-scanning photometers and all-sky TV cameras at Ny Ålesund, Spitzbergen. The data were recorded on 9 January 1989, and a sequence of bursts of flow, with associated transient aurora, were observed between 08:45 and 11:00 U.T. In this paper we concentrate on an event around 09:05 U.T. because that is very close to local magnetic noon. The optical data show a transient intensification and widening (in latitude) of the cusp/cleft region, as seen in red line auroral emissions. Over an interval of about 10 min, the band of 630 nm aurora widened from about 1.5° of invariant latitude to over 5° and returned to its original width. Embedded within the widening band of 630 nm emissions were two intense, active 557.7 nm arc fragments with rays which persisted for about 2 min each. The flow data before and after the optical transient show eastward flows, with speeds increasing markedly with latitude across the band of 630 nm aurora. Strong, apparently westward, flows appeared inside the band while it was widening, but these rotated round to eastward, through northward, as the band shrunk to its original width. The observed ion temperatures verify that the flow speeds during the transient were, to a large extent, as derived using the beamswinging technique; but they also show that the flow increase initially occurred in the western azimuth only. This spatial gradient in the flow introduces ambiguity in the direction of these initial flows and they could have been north-eastward rather than westward. However, the westward direction derived by the beamswinging is consistent with the motion of the colocated and coincident active 557.7 nm arc fragment, A more stable transient 557.7 nm aurora was found close to the shear between the inferred westward flows and the persisting eastward flows to the North. Throughout the transient, northward flow was observed across the equatorward boundary of the 630 nm aurora. Interpretation of the data is made difficult by lack of IMF data, problems in distinguishing the cusp and cleft aurora and uncertainty over which field lines are open and which are closed. However, at magnetic noon there is a 50% probability that we were observing the cusp, in which case from its southerly location we infer that the IMF was southward and many features are suggestive of time-varying reconnection at a single X-line on the dayside magnetopause. This IMF orientation is also consistent with the polar rain precipitation observed simultaneously by the DMSP-F9 satellite in the southern polar cap. There is also a 25% chance that we were observing the cleft (or the mantle poleward of the cleft). In this case we infer that the IMF was northward and the transient is well explained by reconnection which is not only transient in time but occurs at various sites located randomly on the dayside magnetopause (i.e. patchy in space). Lastly, there is a 25% chance that we were observing the cusp poleward of the cleft, in which case we infer that IMF Bz was near zero and the transient is explained by a mixture of the previous two interpretations.

The role of an insulating surface layer on the relaxation time of the ionic redistribution in an electrolytic cell

We analyze the influence of a surface dielectric layer on the transient phenomena related to the ionic redistribution in an electrolytic cell submitted to a step-like external voltage. The adsorption-desorption phenomenon is taken into account in the famework of the Gouy-Chapman approximation, where the ions are assumed dimensionless. In the limit of small amplitude of the applied voltage, where the equations of the problem can be linearized, we obtain an analytical solution for the surface densities of ions, for the electrical potential and for the relaxation time for the transient phenomena. In the general case, when the linearized analysis is no longer valid, the solution of the problem is obtained numerically. The role of the thickness of the dielectric layer on the relaxation time is also discussed.

Coloring ionic trapping states in WO3 and Nb2O5 electrochromic materials

Complex electro-optical analysis is a very useful approach to separate different kinetic processes that occur during ionic insertion reactions in electrochromic oxide materials. In this paper, we use this type of combined technique to investigate ionic and optical changes in different oxide host systems, i.e., in two oxide hosts, specifically WO3 and Nb2O5. A comparison of their electro-optical responses revealed the presence of an ionic trapping contribution to the kinetics of the coloring sites, which was named here as coloring ionic trapping state. As expected, this coloring trapping process is slower in Nb2O5 since the reduction potential of Nb2O5 is more negative (more energy is needed for a higher degree of coloration). A phenomenological solid-state model that encompasses homogeneous charge transfer and valence trapping was proposed to explain the coloring ionic trapping process. Basically the model is able to explain how ionic dynamics at low frequency region, i.e., the slower kinetic step, controls the coloring kinetics, i.e., how it is capable to regulate the coloring rates.Optical transient analyses demonstrated the possibility of the presence of more than one coloring ionic trap, indicating the complexity of the processes involved in coloration phenomenon in metal oxide host systems. (C) 2008 Published by Elsevier Ltd.

Statistical properties of Radio Halos and the re-acceleration model

Galaxy clusters occupy a special position in the cosmic hierarchy as they are the largest bound structures in the Universe. There is now general agreement on a hierarchical picture for the formation of cosmic structures, in which galaxy clusters are supposed to form by accretion of matter and merging between smaller units. During merger events, shocks are driven by the gravity of the dark matter in the diffuse barionic component, which is heated up to the observed temperature. Radio and hard-X ray observations have discovered non-thermal components mixed with the thermal Intra Cluster Medium (ICM) and this is of great importance as it calls for a “revision” of the physics of the ICM. The bulk of present information comes from the radio observations which discovered an increasing number of Mpcsized emissions from the ICM, Radio Halos (at the cluster center) and Radio Relics (at the cluster periphery). These sources are due to synchrotron emission from ultra relativistic electrons diffusing through µG turbulent magnetic fields. Radio Halos are the most spectacular evidence of non-thermal components in the ICM and understanding the origin and evolution of these sources represents one of the most challenging goal of the theory of the ICM. Cluster mergers are the most energetic events in the Universe and a fraction of the energy dissipated during these mergers could be channelled into the amplification of the magnetic fields and into the acceleration of high energy particles via shocks and turbulence driven by these mergers. Present observations of Radio Halos (and possibly of hard X-rays) can be best interpreted in terms of the reacceleration scenario in which MHD turbulence injected during these cluster mergers re-accelerates high energy particles in the ICM. The physics involved in this scenario is very complex and model details are difficult to test, however this model clearly predicts some simple properties of Radio Halos (and resulting IC emission in the hard X-ray band) which are almost independent of the details of the adopted physics. In particular in the re-acceleration scenario MHD turbulence is injected and dissipated during cluster mergers and thus Radio Halos (and also the resulting hard X-ray IC emission) should be transient phenomena (with a typical lifetime <» 1 Gyr) associated with dynamically disturbed clusters. The physics of the re-acceleration scenario should produce an unavoidable cut-off in the spectrum of the re-accelerated electrons, which is due to the balance between turbulent acceleration and radiative losses. The energy at which this cut-off occurs, and thus the maximum frequency at which synchrotron radiation is produced, depends essentially on the efficiency of the acceleration mechanism so that observations at high frequencies are expected to catch only the most efficient phenomena while, in principle, low frequency radio surveys may found these phenomena much common in the Universe. These basic properties should leave an important imprint in the statistical properties of Radio Halos (and of non-thermal phenomena in general) which, however, have not been addressed yet by present modellings. The main focus of this PhD thesis is to calculate, for the first time, the expected statistics of Radio Halos in the context of the re-acceleration scenario. In particular, we shall address the following main questions: • Is it possible to model “self-consistently” the evolution of these sources together with that of the parent clusters? • How the occurrence of Radio Halos is expected to change with cluster mass and to evolve with redshift? How the efficiency to catch Radio Halos in galaxy clusters changes with the observing radio frequency? • How many Radio Halos are expected to form in the Universe? At which redshift is expected the bulk of these sources? • Is it possible to reproduce in the re-acceleration scenario the observed occurrence and number of Radio Halos in the Universe and the observed correlations between thermal and non-thermal properties of galaxy clusters? • Is it possible to constrain the magnetic field intensity and profile in galaxy clusters and the energetic of turbulence in the ICM from the comparison between model expectations and observations? Several astrophysical ingredients are necessary to model the evolution and statistical properties of Radio Halos in the context of re-acceleration model and to address the points given above. For these reason we deserve some space in this PhD thesis to review the important aspects of the physics of the ICM which are of interest to catch our goals. In Chapt. 1 we discuss the physics of galaxy clusters, and in particular, the clusters formation process; in Chapt. 2 we review the main observational properties of non-thermal components in the ICM; and in Chapt. 3 we focus on the physics of magnetic field and of particle acceleration in galaxy clusters. As a relevant application, the theory of Alfv´enic particle acceleration is applied in Chapt. 4 where we report the most important results from calculations we have done in the framework of the re-acceleration scenario. In this Chapter we show that a fraction of the energy of fluid turbulence driven in the ICM by the cluster mergers can be channelled into the injection of Alfv´en waves at small scales and that these waves can efficiently re-accelerate particles and trigger Radio Halos and hard X-ray emission. The main part of this PhD work, the calculation of the statistical properties of Radio Halos and non-thermal phenomena as expected in the context of the re-acceleration model and their comparison with observations, is presented in Chapts.5, 6, 7 and 8. In Chapt.5 we present a first approach to semi-analytical calculations of statistical properties of giant Radio Halos. The main goal of this Chapter is to model cluster formation, the injection of turbulence in the ICM and the resulting particle acceleration process. We adopt the semi–analytic extended Press & Schechter (PS) theory to follow the formation of a large synthetic population of galaxy clusters and assume that during a merger a fraction of the PdV work done by the infalling subclusters in passing through the most massive one is injected in the form of magnetosonic waves. Then the processes of stochastic acceleration of the relativistic electrons by these waves and the properties of the ensuing synchrotron (Radio Halos) and inverse Compton (IC, hard X-ray) emission of merging clusters are computed under the assumption of a constant rms average magnetic field strength in emitting volume. The main finding of these calculations is that giant Radio Halos are naturally expected only in the more massive clusters, and that the expected fraction of clusters with Radio Halos is consistent with the observed one. In Chapt. 6 we extend the previous calculations by including a scaling of the magnetic field strength with cluster mass. The inclusion of this scaling allows us to derive the expected correlations between the synchrotron radio power of Radio Halos and the X-ray properties (T, LX) and mass of the hosting clusters. For the first time, we show that these correlations, calculated in the context of the re-acceleration model, are consistent with the observed ones for typical µG strengths of the average B intensity in massive clusters. The calculations presented in this Chapter allow us to derive the evolution of the probability to form Radio Halos as a function of the cluster mass and redshift. The most relevant finding presented in this Chapter is that the luminosity functions of giant Radio Halos at 1.4 GHz are expected to peak around a radio power » 1024 W/Hz and to flatten (or cut-off) at lower radio powers because of the decrease of the electron re-acceleration efficiency in smaller galaxy clusters. In Chapt. 6 we also derive the expected number counts of Radio Halos and compare them with available observations: we claim that » 100 Radio Halos in the Universe can be observed at 1.4 GHz with deep surveys, while more than 1000 Radio Halos are expected to be discovered in the next future by LOFAR at 150 MHz. This is the first (and so far unique) model expectation for the number counts of Radio Halos at lower frequency and allows to design future radio surveys. Based on the results of Chapt. 6, in Chapt.7 we present a work in progress on a “revision” of the occurrence of Radio Halos. We combine past results from the NVSS radio survey (z » 0.05 − 0.2) with our ongoing GMRT Radio Halos Pointed Observations of 50 X-ray luminous galaxy clusters (at z » 0.2−0.4) and discuss the possibility to test our model expectations with the number counts of Radio Halos at z » 0.05 − 0.4. The most relevant limitation in the calculations presented in Chapt. 5 and 6 is the assumption of an “averaged” size of Radio Halos independently of their radio luminosity and of the mass of the parent clusters. This assumption cannot be released in the context of the PS formalism used to describe the formation process of clusters, while a more detailed analysis of the physics of cluster mergers and of the injection process of turbulence in the ICM would require an approach based on numerical (possible MHD) simulations of a very large volume of the Universe which is however well beyond the aim of this PhD thesis. On the other hand, in Chapt.8 we report our discovery of novel correlations between the size (RH) of Radio Halos and their radio power and between RH and the cluster mass within the Radio Halo region, MH. In particular this last “geometrical” MH − RH correlation allows us to “observationally” overcome the limitation of the “average” size of Radio Halos. Thus in this Chapter, by making use of this “geometrical” correlation and of a simplified form of the re-acceleration model based on the results of Chapt. 5 and 6 we are able to discuss expected correlations between the synchrotron power and the thermal cluster quantities relative to the radio emitting region. This is a new powerful tool of investigation and we show that all the observed correlations (PR − RH, PR − MH, PR − T, PR − LX, . . . ) now become well understood in the context of the re-acceleration model. In addition, we find that observationally the size of Radio Halos scales non-linearly with the virial radius of the parent cluster, and this immediately means that the fraction of the cluster volume which is radio emitting increases with cluster mass and thus that the non-thermal component in clusters is not self-similar.

Wavelet analysis of heart rate variability related to nocturnal frontal lobe epilepsy seizures

Introduction: Nocturnal frontal lobe epilepsy (NFLE) is a distinct syndrome of partial epilepsy whose clinical features comprise a spectrum of paroxysmal motor manifestations of variable duration and complexity, arising from sleep. Cardiovascular changes during NFLE seizures have previously been observed, however the extent of these modifications and their relationship with seizure onset has not been analyzed in detail. Objective: Aim of present study is to evaluate NFLE seizure related changes in heart rate (HR) and in sympathetic/parasympathetic balance through wavelet analysis of HR variability (HRV). Methods: We evaluated the whole night digitally recorded video-polysomnography (VPSG) of 9 patients diagnosed with NFLE with no history of cardiac disorders and normal cardiac examinations. Events with features of NFLE seizures were selected independently by three examiners and included in the study only if a consensus was reached. Heart rate was evaluated by measuring the interval between two consecutive R-waves of QRS complexes (RRi). RRi series were digitally calculated for a period of 20 minutes, including the seizures and resampled at 10 Hz using cubic spline interpolation. A multiresolution analysis was performed (Daubechies-16 form), and the squared level specific amplitude coefficients were summed across appropriate decomposition levels in order to compute total band powers in bands of interest (LF: 0.039062 - 0.156248, HF: 0.156248 - 0.624992). A general linear model was then applied to estimate changes in RRi, LF and HF powers during three different period (Basal) (30 sec, at least 30 sec before seizure onset, during which no movements occurred and autonomic conditions resulted stationary); pre-seizure period (preSP) (10 sec preceding seizure onset) and seizure period (SP) corresponding to the clinical manifestations. For one of the patients (patient 9) three seizures associated with ictal asystole were recorded, hence he was treated separately. Results: Group analysis performed on 8 patients (41 seizures) showed that RRi remained unchanged during the preSP, while a significant tachycardia was observed in the SP. A significant increase in the LF component was instead observed during both the preSP and the SP (p<0.001) while HF component decreased only in the SP (p<0.001). For patient 9 during the preSP and in the first part of SP a significant tachycardia was observed associated with an increased sympathetic activity (increased LF absolute values and LF%). In the second part of the SP a progressive decrease in HR that gradually exceeded basal values occurred before IA. Bradycardia was associated with an increase in parasympathetic activity (increased HF absolute values and HF%) contrasted by a further increase in LF until the occurrence of IA. Conclusions: These data suggest that changes in autonomic balance toward a sympathetic prevalence always preceded clinical seizure onset in NFLE, even when HR changes were not yet evident, confirming that wavelet analysis is a sensitive technique to detect sudden variations of autonomic balance occurring during transient phenomena. Finally we demonstrated that epileptic asystole is associated with a parasympathetic hypertonus counteracted by a marked sympathetic activation.

Método para análise da interação fluido-estrutura em travessas do pré-distribuidor de turbinas hidráulicas.

Um dos grandes desafios enfrentados pelos fabricantes de turbinas hidráulicas é prevenir o aparecimento de vibrações induzidas pelo escoamento nas travessas do pré-distribuidor e pás do rotor. Considerando apenas as travessas, e atribuídos a tais vibrações, foram relatados 28 casos de trincas ou ruídos anormais nas últimas décadas, que acarretaram enormes prejuízos associados a reparos, atrasos e perda de geração. O estado da arte na prevenção destes problemas baseia-se na utilização de sofisticados, e caros, programas comerciais de dinâmica dos fluidos computacional para o cálculo transiente do fenômeno. Este trabalho faz uma ampla revisão bibliográfica e levantamento de eventos de trincas ou ruídos ocorridos em travessas nos últimos 50 anos. Propõe, então, um enfoque alternativo, baseado exclusivamente em ferramentas de código aberto. A partir de hipóteses simplificadoras devidamente justificadas, o problema é formulado matematicamente de forma bidimensional, no plano da seção transversal da travessa, levando em conta a interação fluido-estrutura. Nesta estratégia, as equações de Navier-Stokes são resolvidas pelo método dos elementos finitos por meio da biblioteca gratuita oomph-lib. Um código especial em C++ é desenvolvido para o problema de interação fluido-estrutura, no qual o fenômeno de turbulência é levado em consideração por meio de um algoritmo baseado no modelo de Baldwin-Lomax. O método proposto é validado por meio da comparação dos resultados obtidos com referências e medições disponíveis na literatura, que tratam de problemas de barras retangulares suportadas elasticamente. O trabalho finaliza com a aplicação do método a um estudo de caso envolvendo uma travessa particular.

Solar magnetic carpet III: coronal modelling of synthetic magnetograms

This article is the third in a series working towards the construction of a realistic, evolving, non-linear force-free coronal-field model for the solar magnetic carpet. Here, we present preliminary results of 3D time-dependent simulations of the small-scale coronal field of the magnetic carpet. Four simulations are considered, each with the same evolving photospheric boundary condition: a 48-hour time series of synthetic magnetograms produced from the model of Meyer et al. ( Solar Phys. 272, 29, 2011). Three simulations include a uniform, overlying coronal magnetic field of differing strength, the fourth simulation includes no overlying field. The build-up, storage, and dissipation of magnetic energy within the simulations is studied. In particular, we study their dependence upon the evolution of the photospheric magnetic field and the strength of the overlying coronal field. We also consider where energy is stored and dissipated within the coronal field. The free magnetic energy built up is found to be more than sufficient to power small-scale, transient phenomena such as nanoflares and X-ray bright points, with the bulk of the free energy found to be stored low down, between 0.5 - 0.8 Mm. The energy dissipated is currently found to be too small to account for the heating of the entire quiet-Sun corona. However, the form and location of energy-dissipation regions qualitatively agree with what is observed on small scales on the Sun. Future MHD modelling using the same synthetic magnetograms may lead to a higher energy release.

Numerical methods for solar tomography in the STEREO era

In this thesis, we propose several advances in the numerical and computational algorithms that are used to determine tomographic estimates of physical parameters in the solar corona. We focus on methods for both global dynamic estimation of the coronal electron density and estimation of local transient phenomena, such as coronal mass ejections, from empirical observations acquired by instruments onboard the STEREO spacecraft. We present a first look at tomographic reconstructions of the solar corona from multiple points-of-view, which motivates the developments in this thesis. In particular, we propose a method for linear equality constrained state estimation that leads toward more physical global dynamic solar tomography estimates. We also present a formulation of the local static estimation problem, i.e., the tomographic estimation of local events and structures like coronal mass ejections, that couples the tomographic imaging problem to a phase field based level set method. This formulation will render feasible the 3D tomography of coronal mass ejections from limited observations. Finally, we develop a scalable algorithm for ray tracing dense meshes, which allows efficient computation of many of the tomographic projection matrices needed for the applications in this thesis.