Study of Stability of Rotational Motion of Spacecraft with Canonical Variables

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Spin-Stabilized Spacecrafts: Analytical Attitude Propagation Using Magnetic Torques

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Vibration energy harvesting using piezoelectric transducer and non-controlled rectifiers circuits

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Crossover between macroscopic and mesoscopic regimes of vortex interactions in type-II superconductors

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

On an Optimal Linear Control of a Chaotic Non-Ideal Duffing System

In this work, we use a nonlinear control based on Optimal Linear Control. We used as mathematical model a Duffing equation to model a supporting structure for an unbalanced rotating machine with limited power (non-ideal motor). Numerical simulations are performed for a set control parameter (depending on the voltage of the motor, that is, in the static and dynamic characteristic of the motor) The interaction of the non-ideal excitation with the structure may lead to the occurrence of interesting phenomena during the forward passage through the several resonance states of the system. Chaotic behavior is obtained for values of the parameters. Then, the proposed control strategy is applied in order to regulate the chaotic behavior, in order to obtain a periodic orbit and to decrease its amplitude. Both methodologies were used in complete agreement between them. The purpose of the paper is to give suggestions and recommendations to designers and engineers on how to drive this kind of system through resonance.

Dinâmica e estrutura de redes complexas no modelo de afinidade

In this work we elaborate and discuss a Complex Network model which presents connectivity scale free probability distribution (power-law degree distribution). In order to do that, we modify the rule of the preferential attachment of the Bianconi-Barabasi model, including a factor which represents the similarity of the sites. The term that corresponds to this similarity is called the affinity, and is obtained by the modulus of the difference between the fitness (or quality) of the sites. This variation in the preferential attachment generates very interesting results, by instance the time evolution of the connectivity, which follows a power-law distribution ki / ( t t0 )fi, where fi indicates the rate to the site gain connections. Certainly this depends on the affinity with other sites. Besides, we will show by numerical simulations results for the average path length and for the clustering coefficient

Efeitos de Confinamento de Vórtice Magnético em Nanodiscos de Fe e Py

We report two theoretical works, based in numerical simulations. The first study consists in the investigation of equilibrium phases and vortex formation in Ferro and Permalloy circular and square nanoelements.The another have the aim to investigate the magnetostatic interaction between pairs of nanodisks of Ferro and Permalloy and it`s impact in the vortex structure

Nucleação de vortices e paredes de domínio em nanoestruturas magnéticas

We report a theoretical investigation of the magnetic phases and hysteresis of exchange biased ferromagnetic (F) nanoelements for three di erent systems: exchange biased nanoparticles, exchange biased narrow ferromagnetic stripes and exchange biased thin ferromagnetic lms. In all cases the focus is on the new e ects produced by suitable patterns of the exchange energy coupling the ferromagnetic nanoelement with a large anisotropy antiferromagnetic (AF) substrate. We investigate the hysteresis of iron and permalloy nanoparticles with a square basis, with lateral dimensions between 45 nm and 120 nm and thickness between 12 nm and 21 nm. Interface bias is aimed at producing large domains in thin lms. Our results show that, contrary to intuition, the interface exchange coupling may generate vortex states along the hysteresis loop. Also, the threshold value of the interface eld strength for vortex nucleation is smaller for iron nanoelements. We investigate the nucleation and depinning of an array of domain walls pinned at interface defects of a vicinal stripe/AF bilayer. The interface exchange eld displays a periodic pattern corresponding to the topology of the AF vicinal substrate. The vicinal AF substrate consists of a sequence of terraces, each with spins from one AF subalattice, alternating one another. As a result the interface eld of neighboring terraces point in opposite direction, leading to the nucleation of a sequence of domain walls in the ferromagnetic stripe. We investigated iron an permalloy micrometric stripes, with width ranging from 100 nm and 300 nm and thickness of 5 nm. We focused in domain wall sequences with same chirality and alternate chirality. We have found that for 100nm terraces the same chiraility sequence is more stable, requiring a larger value of the external eld for depinning. The third system consists of an iron lm with a thickness of 5 nm, exchange coupled to an AF substrate with a periodic distribution of islands where the AF spins have the opposite direction of the spins in the background. This corresponds to a two-sublattice noncompensated AF plane (such as the surface of a (100) FeF2 lm), with monolayer-height islands containing spins of one sublattice on a surface containing spins of the opposite sublattice. The interface eld acting in the ferromagnetic spins over the islands points in the opposite direction of that in the spins over the background. This a model system for the investigation of interface roughness e ects. We have studied the coercicivity an exchange bias hysteresis shift as a function of the distance between the islands and the degree of interface roughness. We have found a relevant reduction of coercivity for nearly compensated interfaces. Also the e ective hysteresis shift is not proportional to the liquid moment of the AF plane. We also developed an analytical model which reproduces qualitatively the results of numerical simulations

Modelagem Matemática e Computacional de fenômenos eletrocinéticos em meios porosos carregados eletricamente

In this work we present a mathematical and computational modeling of electrokinetic phenomena in electrically charged porous medium. We consider the porous medium composed of three different scales (nanoscopic, microscopic and macroscopic). On the microscopic scale the domain is composed by a porous matrix and a solid phase. The pores are filled with an aqueous phase consisting of ionic solutes fully diluted, and the solid matrix consists of electrically charged particles. Initially we present the mathematical model that governs the electrical double layer in order to quantify the electric potential, electric charge density, ion adsorption and chemical adsorption in nanoscopic scale. Then, we derive the microscopic model, where the adsorption of ions due to the electric double layer and the reactions of protonation/ deprotanaç~ao and zeta potential obtained in modeling nanoscopic arise in microscopic scale through interface conditions in the problem of Stokes and Nerst-Planck equations respectively governing the movement of the aqueous solution and transport of ions. We developed the process of upscaling the problem nano/microscopic using the homogenization technique of periodic structures by deducing the macroscopic model with their respectives cell problems for effective parameters of the macroscopic equations. Considering a clayey porous medium consisting of kaolinite clay plates distributed parallel, we rewrite the macroscopic model in a one-dimensional version. Finally, using a sequential algorithm, we discretize the macroscopic model via the finite element method, along with the interactive method of Picard for the nonlinear terms. Numerical simulations on transient regime with variable pH in one-dimensional case are obtained, aiming computational modeling of the electroremediation process of clay soils contaminated

Parâmetro de regularização em problemas inversos: estudo numérico com a transformada de Radon

In general, an inverse problem corresponds to find a value of an element x in a suitable vector space, given a vector y measuring it, in some sense. When we discretize the problem, it usually boils down to solve an equation system f(x) = y, where f : U Rm ! Rn represents the step function in any domain U of the appropriate Rm. As a general rule, we arrive to an ill-posed problem. The resolution of inverse problems has been widely researched along the last decades, because many problems in science and industry consist in determining unknowns that we try to know, by observing its effects under certain indirect measures. Our general subject of this dissertation is the choice of Tykhonov´s regulaziration parameter of a poorly conditioned linear problem, as we are going to discuss on chapter 1 of this dissertation, focusing on the three most popular methods in nowadays literature of the area. Our more specific focus in this dissertation consists in the simulations reported on chapter 2, aiming to compare the performance of the three methods in the recuperation of images measured with the Radon transform, perturbed by the addition of gaussian i.i.d. noise. We choosed a difference operator as regularizer of the problem. The contribution we try to make, in this dissertation, mainly consists on the discussion of numerical simulations we execute, as is exposed in Chapter 2. We understand that the meaning of this dissertation lays much more on the questions which it raises than on saying something definitive about the subject. Partly, for beeing based on numerical experiments with no new mathematical results associated to it, partly for being about numerical experiments made with a single operator. On the other hand, we got some observations which seemed to us interesting on the simulations performed, considered the literature of the area. In special, we highlight observations we resume, at the conclusion of this work, about the different vocations of methods like GCV and L-curve and, also, about the optimal parameters tendency observed in the L-curve method of grouping themselves in a small gap, strongly correlated with the behavior of the generalized singular value decomposition curve of the involved operators, under reasonably broad regularity conditions in the images to be recovered

Desenvolvimento de dobras e falhas em ambiente distensional: Aplicação da modelagem física

The geological modeling allows, at laboratory scaling, the simulation of the geometric and kinematic evolution of geological structures. The importance of the knowledge of these structures grows when we consider their role in the creation of traps or conduits to oil and water. In the present work we simulated the formation of folds and faults in extensional environment, through physical and numerical modeling, using a sandbox apparatus and MOVE2010 software. The physical modeling of structures developed in the hangingwall of a listric fault, showed the formation of active and inactive axial zones. In consonance with the literature, we verified the formation of a rollover between these two axial zones. The crestal collapse of the anticline formed grabens, limited by secondary faults, perpendicular to the extension, with a curvilinear aspect. Adjacent to these faults we registered the formation of transversal folds, parallel to the extension, characterized by a syncline in the fault hangingwall. We also observed drag folds near the faults surfaces, these faults are parallel to the fault surface and presented an anticline in the footwall and a syncline hangingwall. To observe the influence of geometrical variations (dip and width) in the flat of a flat-ramp fault, we made two experimental series, being the first with the flat varying in dip and width and the second maintaining the flat variation in width but horizontal. These experiments developed secondary faults, perpendicular to the extension, that were grouped in three sets: i) antithetic faults with a curvilinear geometry and synthetic faults, with a more rectilinear geometry, both nucleated in the base of sedimentary pile. The normal antithetic faults can rotate, during the extension, presenting a pseudo-inverse kinematics. ii) Faults nucleated at the top of the sedimentary pile. The propagation of these faults is made through coalescence of segments, originating, sometimes, the formation of relay ramps. iii) Reverse faults, are nucleated in the flat-ramp interface. Comparing the two models we verified that the dip of the flat favors a differentiated nucleation of the faults at the two extremities of the mater fault. V These two flat-ramp models also generated an anticline-syncline pair, drag and transversal folds. The anticline was formed above the flat being sub-parallel to the master fault plane, while the syncline was formed in more distal areas of the fault. Due the geometrical variation of these two folds we can define three structural domains. Using the physical experiments as a template, we also made numerical modeling experiments, with flat-ramp faults presenting variation in the flat. Secondary antithetic, synthetic and reverse faults were generated in both models. The numerical modeling formed two folds, and anticline above the flat and a syncline further away of the master fault. The geometric variation of these two folds allowed the definition of three structural domains parallel to the extension. These data reinforce the physical models. The comparisons between natural data of a flat-ramp fault in the Potiguar basin with the data of physical and numerical simulations, showed that, in both cases, the variation of the geometry of the flat produces, variation in the hangingwall geometry

Modelling the lethargic crab disease.

The lethargic crab disease (LCD) is an emergent infirmity that has decimated native populations of the mangrove land crab (Ucides cordatus, Decapoda: Ocypodidae) along the Brazilian coast. Several potential etiological agents have been linked with LCD, but only in 2005 was it proved that it is caused by an ascomycete fungus. This is the first attempt to develop a mathematical model to describe the epidemiological dynamics of LCD. The model presents four possible scenarios, namely, the trivial equilibrium, the disease-free equilibrium, endemic equilibrium, and limit cycles arising from a Hopf bifurcation. The threshold values depend on the basic reproductive number of crabs and fungi, and on the infection rate. These scenarios depend on both the biological assumptions and the temporal evolution of the disease. Numerical simulations corroborate the analytical results and illustrate the different temporal dynamics of the crab and fungus populations.

CHARGED PARTICLE DYNAMICS IN TIME DEPENDENT MAGNETIC FIELDS

In this work we study the behavior of charged particles immersed in a peculiar configuration of magnetic fields, which has a main constant field B(0) and a superimposed, transversal perturbation field B(1) sin(omega(p)t), with B(1) << B(0). By taking Cartesian coordinates and placing B(0) along the z axis and B(1) sin (omega(p)t) on the x axis, an analytical solution for y(t) may be obtained by solving an integrodifferential equation. Besides, the solution z(t) also exhibits a very interesting dynamics, and the entire system is conditioned by resonances between the particle orbit frequencies and the frequency of the magnetic transversal perturbation, omega(p). In this work we also discuss numerical simulations for the related particle trajectories, as well as potential applications in the context of separation phenomena.

Dynamics of the deconfinement transition of quarks and gluons in a finite volume

We employ finite elements methods for the approximation of solutions of the Ginzburg-Landau equations describing the deconfinement transition in quantum chromodynamics. These methods seem appropriate for situations where the deconfining transition occurs over a finite volume as in relativistic heavy ion collisions. where in addition expansion of the system and flow of matter are important. Simulation results employing finite elements are presented for a Ginzburg-Landau equation based on a model free energy describing the deconfining transition in pure gauge SU(2) theory. Results for finite and infinite system are compared. (C) 2009 Elsevier B.V. All rights reserved.

Dissipationless shock waves in Bose-Einstein condensates with repulsive interaction between atoms

We consider formation of dissipationless shock waves in Bose-Einstein condensates with repulsive interaction between atoms. It is shown that for big enough initial inhomogeneity of density, interplay of nonlinear and dispersion effects leads to wave breaking phenomenon followed by generation of a train of dark solitons. Analytical theory is confirmed by numerical simulations.