On the cause of a magnetospheric flux transfer event

We present a detailed investigation of a magnetospheric flux transfer event (FTE) seen by the Active Magnetospheric Tracer Explorer (AMPTE) UKS and IRM satellites around 1046 UT on October 28, 1984. This event has been discussed many times previously in the literature and has been cited as support for a variety of theories of FTE formation. We make use of a model developed to reproduce ion precipitations seen in the cusp ionosphere. The analysis confirms that the FTE is well explained as a brief excursion into an open low-latitude boundary layer (LLBL), as predicted by two theories of magnetospheric FTEs: namely, that they are bulges in the open LLBL due to reconnection rate enhancements or that they are indentations of the magnetopause by magnetosheath pressure increases (but in the presence of ongoing steady reconnection). The indentation of the inner edge of the open LLBL that these two models seek to explain is found to be shallow for this event. The ion model reproduces the continuous evolution of the ion distribution function between the sheath-like population at the event center and the surrounding magnetospheric populations; it also provides an explanation of the high-pressure core of the event as comprising field lines that were reconnected considerably earlier than those that are draped over it to give the event boundary layer. The magnetopause transition parameter is used to isolate a field rotation on the boundaries of the core, which is subjected to the tangential stress balance test. The test identifies this to be a convecting structure, which is neither a rotational discontinuity (RD) nor a contact discontinuity, but could possibly be a slow shock. In addition, evidence for ion reflection off a weak RD on the magnetospheric side of this structure is found. The event structure is consistent in many ways with features predicted for the open LLBL by analytic MHD theories and by MHD and hybrid simulations. The de Hoffman-Teller velocity of the structure is significantly different from that of the magnetosheath flow, indicating that it is not an indentation caused by a high-pressure pulse in the sheath but is consistent with the motion of newly opened field lines (different from the sheath flow because of the magnetic tension force) deduced from the best fit to the ion data. However, we cannot here rule out the possibility that the sheath flow pattern has changed in the long interval between the two satellites observing the FTE and subsequently emerging into the magnetosheath; thus this test is not conclusive in this particular case. Analysis of the fitted elapsed time since reconnection shows that the core of the event was reconnected in one pulse and the event boundary layer was reconnected in a subsequent pulse. Between these two pulses is a period of very low (but nonzero) reconnection rate, which lasts about 14 mins. Thus the analysis supports, but does not definitively verify, the concept that the FTE is a partial passage into an open LLBL caused by a traveling bulge in that layer produced by a pulse in reconnection rate.

The case for transient magnetopause reconnection

Recent studies indicate that magnetopause reconnection can, at times, occur exclusively as a series of separated, short bursts. Reconnection generates “open” magnetic flux that threads the magnetospheric boundary, the magnetopause, and so connects the magnetosphere with interplanetary space. The rate at which open flux is generated by a line of unit length in the magnetopause is called the reconnection rate. By Faraday's induction law, the reconnection rate is a boundary-tangential electric field, Et, along that line (called an X-line).

Low and middle altitude cusp particle signatures for general magnetopause reconnection rate variations: 1. Theory

We present predictions of the signatures of magnetosheath particle precipitation (in the regions classified as open low-latitude boundary layer, cusp, mantle and polar cap) for periods when the interplanetary magnetic field has a southward component. These are made using the “pulsating cusp” model of the effects of time-varying magnetic reconnection at the dayside magnetopause. Predictions are made for both low-altitude satellites in the topside ionosphere and for midaltitude spacecraft in the magnetosphere. Low-altitude cusp signatures, which show a continuous ion dispersion signature, reveal "quasi-steady reconnection" (one limit of the pulsating cusp model), which persists for a period of at least 10 min. We estimate that “quasi-steady” in this context corresponds to fluctuations in the reconnection rate of a factor of 2 or less. The other limit of the pulsating cusp model explains the instantaneous jumps in the precipitating ion spectrum that have been observed at low altitudes. Such jumps are produced by isolated pulses of reconnection: that is, they are separated by intervals when the reconnection rate is zero. These also generate convecting patches on the magnetopause in which the field lines thread the boundary via a rotational discontinuity separated by more extensive regions of tangential discontinuity. Predictions of the corresponding ion precipitation signatures seen by midaltitude spacecraft are presented. We resolve the apparent contradiction between estimates of the width of the injection region from midaltitude data and the concept of continuous entry of solar wind plasma along open field lines. In addition, we reevaluate the use of pitch angle-energy dispersion to estimate the injection distance.

Comparison of terrain following and cut cell grids using a non-hydrostatic model

Terrain following coordinates are widely used in operational models but the cut cell method has been proposed as an alternative that can more accurately represent atmospheric dynamics over steep orography. Because the type of grid is usually chosen during model implementation, it becomes necessary to use different models to compare the accuracy of different grids. In contrast, here a C-grid finite volume model enables a like-for-like comparison of terrain following and cut cell grids. A series of standard two-dimensional tests using idealised terrain are performed: tracer advection in a prescribed horizontal velocity field, a test starting from resting initial conditions, and orographically induced gravity waves described by nonhydrostatic dynamics. In addition, three new tests are formulated: a more challenging resting atmosphere case, and two new advection tests having a velocity field that is everywhere tangential to the terrain following coordinate surfaces. These new tests present a challenge on cut cell grids. The results of the advection tests demonstrate that accuracy depends primarily upon alignment of the flow with the grid rather than grid orthogonality. A resting atmosphere is well-maintained on all grids. In the gravity waves test, results on all grids are in good agreement with existing results from the literature, although terrain following velocity fields lead to errors on cut cell grids. Due to semi-implicit timestepping and an upwind-biased, explicit advection scheme, there are no timestep restrictions associated with small cut cells. We do not find the significant advantages of cut cells or smoothed coordinates that other authors find.

Recording from Two Neurons: Second-Order Stimulus Reconstruction from Spike Trains and Population Coding

We study the reconstruction of visual stimuli from spike trains, representing the reconstructed stimulus by a Volterra series up to second order. We illustrate this procedure in a prominent example of spiking neurons, recording simultaneously from the two H1 neurons located in the lobula plate of the fly Chrysomya megacephala. The fly views two types of stimuli, corresponding to rotational and translational displacements. Second-order reconstructions require the manipulation of potentially very large matrices, which obstructs the use of this approach when there are many neurons. We avoid the computation and inversion of these matrices using a convenient set of basis functions to expand our variables in. This requires approximating the spike train four-point functions by combinations of two-point functions similar to relations, which would be true for gaussian stochastic processes. In our test case, this approximation does not reduce the quality of the reconstruction. The overall contribution to stimulus reconstruction of the second-order kernels, measured by the mean squared error, is only about 5% of the first-order contribution. Yet at specific stimulus-dependent instants, the addition of second-order kernels represents up to 100% improvement, but only for rotational stimuli. We present a perturbative scheme to facilitate the application of our method to weakly correlated neurons.

Efficient solutions to robust, semi-implicit discretizations of the immersed boundary method

The immersed boundary method is a versatile tool for the investigation of flow-structure interaction. In a large number of applications, the immersed boundaries or structures are very stiff and strong tangential forces on these interfaces induce a well-known, severe time-step restriction for explicit discretizations. This excessive stability constraint can be removed with fully implicit or suitable semi-implicit schemes but at a seemingly prohibitive computational cost. While economical alternatives have been proposed recently for some special cases, there is a practical need for a computationally efficient approach that can be applied more broadly. In this context, we revisit a robust semi-implicit discretization introduced by Peskin in the late 1970s which has received renewed attention recently. This discretization, in which the spreading and interpolation operators are lagged. leads to a linear system of equations for the inter-face configuration at the future time, when the interfacial force is linear. However, this linear system is large and dense and thus it is challenging to streamline its solution. Moreover, while the same linear system or one of similar structure could potentially be used in Newton-type iterations, nonlinear and highly stiff immersed structures pose additional challenges to iterative methods. In this work, we address these problems and propose cost-effective computational strategies for solving Peskin`s lagged-operators type of discretization. We do this by first constructing a sufficiently accurate approximation to the system`s matrix and we obtain a rigorous estimate for this approximation. This matrix is expeditiously computed by using a combination of pre-calculated values and interpolation. The availability of a matrix allows for more efficient matrix-vector products and facilitates the design of effective iterative schemes. We propose efficient iterative approaches to deal with both linear and nonlinear interfacial forces and simple or complex immersed structures with tethered or untethered points. One of these iterative approaches employs a splitting in which we first solve a linear problem for the interfacial force and then we use a nonlinear iteration to find the interface configuration corresponding to this force. We demonstrate that the proposed approach is several orders of magnitude more efficient than the standard explicit method. In addition to considering the standard elliptical drop test case, we show both the robustness and efficacy of the proposed methodology with a 2D model of a heart valve. (C) 2009 Elsevier Inc. All rights reserved.

The use of polyacrylamides (PAMs) for removing natural organic matter (NOM) from Paraiba do Sul River (Brazil)

This work assesses the efficiency of polyacrylamides for natural organic matter (NOM) removal from Paraiba do Sul River (Brazil) raw water for drinking purposes. Jar tests were performed following an experimental design protocol. Three kinds of polyacrylamides (anionic, cationic, and non-ionic) at 0.2 mg L(-1) were tested. After coagulation, turbidity, DOC, UVA(254) and SCAN (UV-absorbing material) were determined. Color and pH were also measured. It was found that polyacrylamides did not reduce the amounts of alum and lime needed in the process and that the amount of alum alone for removing UV-absorbing organic matter is significantly higher. Efficiency of the coagulation process decreased as follows: non-ionic -> cationic -> anionic -> no polyacrylamide. Removal efficiencies for the best case were: 100%, 90%, 83%, and 68% for turbidity, DOC, UVA(254), and SCAN, respectively.

Identification of cause of impairment in spiral drawings, using non-stationary feature extraction approach

Parkinson’s disease is a clinical syndrome manifesting with slowness and instability. As it is a progressive disease with varying symptoms, repeated assessments are necessary to determine the outcome of treatment changes in the patient. In the recent past, a computer-based method was developed to rate impairment in spiral drawings. The downside of this method is that it cannot separate the bradykinetic and dyskinetic spiral drawings. This work intends to construct the computer method which can overcome this weakness by using the Hilbert-Huang Transform (HHT) of tangential velocity. The work is done under supervised learning, so a target class is used which is acquired from a neurologist using a web interface. After reducing the dimension of HHT features by using PCA, classification is performed. C4.5 classifier is used to perform the classification. Results of the classification are close to random guessing which shows that the computer method is unsuccessful in assessing the cause of drawing impairment in spirals when evaluated against human ratings. One promising reason is that there is no difference between the two classes of spiral drawings. Displaying patients self ratings along with the spirals in the web application is another possible reason for this, as the neurologist may have relied too much on this in his own ratings.

Modelo Estocástico para bloqueio de poros e redução de permeabilidade

Modeling transport of particulate suspensions in porous media is essential for understanding various processes of industrial and scientific interest. During these processes, particles are retained due to mechanisms like size exclusion (straining), adsorption, sedimentation and diffusion. In this thesis, a mathematical model is proposed and analytical solutions are obtained. The obtained analytic solutions for the proposed model, which takes pore and particle size distributions into account, were applied to predict the particle retention, pore blocking and permeability reduction during dead-end microfiltration in membranes. Various scenarios, considering different particle and pore size distributions were studied. The obtained results showed that pore blocking and permeability reduction are highly influenced by the initial pore and particle size distributions. This feature was observed even when different initial pore and particle size distributions with the same average pore size and injected particle size were considered. Finally, a mathematical model for predicting equivalent permeability in porous media during particle retention (and pore blocking) is proposed and the obtained solutions were applied to study permeability decline in different scenarios

Influência na qualidade final de metais retificados através da variação da velocidade de mergulho

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Influência da lubrirrefrigeração na qualidade superficial de metais retificados

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Desenvolvimento de antenas de microfita com Patch em anel utilizando materiais ferrimagnéticos e metamateriais

In general, the materials used as substrates in the project of microstrip antennas are: isotropic, anisotropic dielectrics and ferrimagnetic materials (magnetic anisotropy). The use of ferrimagnetic materials as substrates in microstrip patch antennas has been concentrated on the analysis of antennas with circular and rectangular patches. However, a new class of materials, called metamaterials, has been currently the focus of a great deal of interest. These materials exhibit bianisotropic characteristics, with permittivity and permeability tensors. The main objective of this work is to develop a theoretical and numerical analysis for the radiation characteristics of annular ring microstrip antennas, using ferrites and metamaterials as substrates. The full wave analysis is performed in the Hankel transform domain through the application of the Hertz vector potentials. Considering the definition of the Hertz potentials and imposing the boundary conditions, the dyadic Green s function components are obtained relating the surface current density components at the plane of the patch to the electric field tangential components. Then, Galerkin s method is used to obtain a system of matrix equations, whose solution gives the antenna resonant frequency. From this modeling, it is possible to obtain numerical results for the resonant frequency, radiation pattern, return loss, and antenna bandwidth as a function of the annular ring physical parameters, for different configurations and substrates. The theoretical analysis was developed for annular ring microstrip antennas on a double ferrimagnetic/isotropic dielectric substrate or metamaterial/isotropic dielectric substrate. Also, the analysis for annular ring microstrip antennas on a single ferrimagnetic or metamaterial layer and for suspended antennas can be performed as particular cases

Propriedades físicas dos lenhos juvenil e adulto de Pinus elliottii Engelm var. elliottii e de Eucalyptus grandis Hill ex Maiden

The juvenile wood presents great difference in its anatomical and physical properties in relation to the mature wood. That heterogeneity of the wood causes a series of upset for the processing industry. The present research had as objective the study of the physical properties of the juvenile and mature wood of Pinus elliottii Engelm var. elliottii, with 25 years of age and of Eucalyptus grandis Hill ex Maiden, with 30 years of age, coming from the Experimental Station of Itapetininga and State Forest of Santa Barbara, belonging to the Forest Institute of the state of São Paulo. From both species there were taken from radial boards, specimens with dimensions of 20x30x50 mm, for the determination of density, swelling and swelling coefficient of the juvenile and mature wood. The results of the species showed that the mature compared with the juvenile wood presented: (1) an increase in the values of the basic density, at 0% and at 12% of moisture content and of the volumetric, radial and tangential swelling at 12% moisture content and maximum coefficient of radial and tangential swelling; (2) a reduction in the values of longitudinal swelling at 12% and maximum and of the longitudinal swelling coefficient.

Caracterização de antenas planares com substrato metamaterial

This work presents a theoretical and numerical analysis for the radiation characteristics of rectangular microstrip antenna using metamaterial substrate. The full wave analysis is performed in the Fourier transform domain through the application of the Transverse Transmission Line - TTL method. A study on metamaterial theory was conducted to obtain the constructive parameters, which were characterized through permittivity and permeability tensors to arrive at a set of electromagnetic equations. The general equations for the electromagnetic fields of the antenna are developed using the Transverse Transmission Line - TTL method. Imposing the boundary conditions, the dyadic Green s function components are obtained relating the surface current density components at the plane of the patch to the electric field tangential components. Then, Galerkin s method is used to obtain a system of matrix equations, whose solution gives the antenna resonant frequency. From this modeling, it is possible to obtain numerical results for the resonant frequency and return loss for different configurations and substrates

Antenas de microfita com patch em anel e múltiplas camadas dielétricas anisotrópicas uniaxiais

This work presents an analysis of the annular ring microstrip antennas printed on uniaxial anisotropic substrates and with superstrate.The analysis uses the full-wave formulation by means of the Hertz vector potentials method, in the Hankel transform domain. The definition of the Hertz vector potentials and the application of the appropriate boundary conditions to the structure allow determining the dyadic Green functions, relating the current densities in the conducting patch to the transforms of the tangential electric field components. Galerkin s method is then used to obtain the matrix equation whose nontrivial solution gives the complex resonant frequency of the antenna. From the modeling, it is possible to obtain results for the resonant frequency, bandwidth and quality factor, as a function of several parameters of the antenna, for different configurations. We have considered annular ring microstrip antennas on a single dielectric layer, antennas with two anisotropic dielectric layers, and annular ring microstrip antennas on suspended substrates. Numerical results for the resonant frequency of the these structures printed on isotropic substrates are also presented and compared with those published by other authors, showing a good agreement