No description available.

We are interested in the onset of instability of the axisymmetric flow between two concentric spherical shells that differentially rotate about a common axis in the narrow-gap limit. The expected mode of instability takes the form of roughly square axisymmetric Taylor vortices which arise in the vicinity of the equator and are modulated on a latitudinal length scale large compared to the gap width but small compared to the shell radii. At the heart of the difficulties faced is the presence of phase mixing in the system, characterised by a non-zero frequency gradient at the equator and the tendency for vortices located off the equator to oscillate. This mechanism serves to enhance viscous dissipation in the fluid with the effect that the amplitude of any initial disturbance generated at onset is ultimately driven to zero. In this thesis we study a complex Ginzburg-Landau equation derived from the weakly nonlinear analysis of Harris, Bassom and Soward [D. Harris, A. P. Bassom, A. M. Soward, Global bifurcation to travelling waves with application to narrow gap spherical Couette flow, Physica D 177 (2003) p. 122-174] (referred to as HBS) to govern the amplitude modulation of Taylor vortex disturbances in the vicinity of the equator. This equation was developed in a regime that requires the angular velocities of the bounding spheres to be very close. When the spherical shells do not co-rotate, it has the remarkable property that the linearised form of the equation has no non-trivial neutral modes. Furthermore no steady solutions to the nonlinear equation have been found. Despite these challenges Bassom and Soward [A. P. Bassom, A. M. Soward, On finite amplitude subcritical instability in narrow-gap spherical Couette flow, J. Fluid Mech. 499 (2004) p. 277-314] (referred to as BS) identified solutions to the equation in the form of pulse-trains. These pulse-trains consist of oscillatory finite amplitude solutions expressed in terms of a single complex amplitude localised as a pulse about the origin. Each pulse oscillates at a frequency proportional to its distance from the equatorial plane and the whole pulse-train is modulated under an envelope and drifts away from the equator at a relatively slow speed. The survival of the pulse-train depends upon the nonlinear mutual-interaction of close neighbours; as the absence of steady solutions suggests, self-interaction is inadequate. Though we report new solutions to the HBS co-rotation model the primary focus in this work is the physically more interesting case when the shell velocities are far from close. More specifically we concentrate on the investigation of BS-style pulse-train solutions and, in the first part of this thesis, develop a generic framework for the identification and classification of pulse-train solutions. Motivated by relaxation oscillations identified by Cole [S. J. Cole, Nonlinear rapidly rotating spherical convection, Ph.D. thesis, University of Exeter (2004)] whilst studying the related problem of thermal convection in a rapidly rotating self-gravitating sphere, we extend the HBS equation in the second part of this work. A model system is developed which captures many of the essential features exhibited by Cole's, much more complicated, system of equations. We successfully reproduce relaxation oscillations in this extended HBS model and document the solution as it undergoes a series of interesting bifurcations.

Resumo: O interesse em conhecer o comportamento supercondutor tem sido cada vez maior nas últimas décadas. Na busca de melhores características supercondutoras, descobriu-se que amostras volumétricas apresentam características muito diferentes de amostras mesoscópicas (amostras com dimensões próximas dos comprimentos de penetração de London e coerência). Como exemplo, podemos citar a não formação de rede de Abrikosov, como consequência do efeito de confinamento (efeito associado às dimensões reduzidas da amostra) e também uma mudança considerável nos valores dos campos críticos. Neste trabalho foram resolvidas as equações de Ginzburg-Landau dependentes do tempo (TDGL), para fazer uma análise detalhada da dinâmica de vórtices em filmes finos mesoscópicos. Para revolvê-las, utilizamos o método das variáveis de ligação com invariância de calibre, adaptado para o algoritmo de diferenças finitas, utilizado para obter a densidade dos pares de Cooper e também curvas de magnetização. O estudo dessa dinâmica de vórtices, foi feito em três amostras com superfícies geométricas diferentes (côncova, convexa e rugosa). Observamos que na comparação entre as duas primeiras, há uma diferença considerável nos valores dos campos críticos, bem como no comportamento da magnetização comparado com um filme plano. Já para a amostra de superfície rugosa, observamos que existe uma competição entre o efeito de confinamento e a rugosidade em relação à configuração dos vórtices. Apresentamos também, uma tabela que mostra resumidamente os estados estacionários dos vórtices nas três amostras. / Abstract: The interest to investigate the investigate the behavior of a superconductor has grown in the last few decades. Having in mind to search for better superconducting characteristics, it has been found that bulk samples present characteristics much more different than mesoscopic samples (samples with dimensions of the same order of the same order of the London penetration length and the coherence length). As an example, we can mention the non-formation of an Abrikosov vortex lattice as a consequence of the confinement effect (effect associated with the reduced dimensions of the sample) and also considerable change in the critical field values. In the present work we solved the time dependent Ginzburg-Landau equation (TDGL), in order to make a detailed analysis of the vortex dynamics in mesoscopic thin films. To solve these equations, we have used the link variables method which is gauge invariant. From this, we obtain the Cooper pair density and the magnetization curves. The vortex dynamics was investigated for three different surfaces of the film (concave, convex, and irregular). We have observed that, with respect to the parabolic geometries, there is a considerable difference for the critical fields, as well as for the behavior of the magnetization compared to a flat film. On the other hand, for a sample with an irregular surface, we have seen that there is a competition between the confinement effect and rugosity with respect to vortex configurations. We also present a table which summarizes the vortex stationary states for the three topologies mentioned above. / Orientador: Paulo Noronha Lisboa Filho / Coorientador: Edson Sardella / Banca: Wilson Aires Ortiz / Banca: Clelio Clemente de Souza Silva / Mestre

Made available in DSpace on 2014-06-11T19:30:19Z (GMT). No. of bitstreams: 0 Previous issue date: 2010-04-28Bitstream added on 2014-06-13T19:39:47Z : No. of bitstreams: 1 pascolati_mcv_me_bauru.pdf: 4047066 bytes, checksum: cf327b8f4cd0447dc8dd5e7c4d6604a5 (MD5) / O interesse em conhecer o comportamento supercondutor tem sido cada vez maior nas últimas décadas. Na busca de melhores características supercondutoras, descobriu-se que amostras volumétricas apresentam características muito diferentes de amostras mesoscópicas (amostras com dimensões próximas dos comprimentos de penetração de London e coerência). Como exemplo, podemos citar a não formação de rede de Abrikosov, como consequência do efeito de confinamento (efeito associado às dimensões reduzidas da amostra) e também uma mudança considerável nos valores dos campos críticos. Neste trabalho foram resolvidas as equações de Ginzburg-Landau dependentes do tempo (TDGL), para fazer uma análise detalhada da dinâmica de vórtices em filmes finos mesoscópicos. Para revolvê-las, utilizamos o método das variáveis de ligação com invariância de calibre, adaptado para o algoritmo de diferenças finitas, utilizado para obter a densidade dos pares de Cooper e também curvas de magnetização. O estudo dessa dinâmica de vórtices, foi feito em três amostras com superfícies geométricas diferentes (côncova, convexa e rugosa). Observamos que na comparação entre as duas primeiras, há uma diferença considerável nos valores dos campos críticos, bem como no comportamento da magnetização comparado com um filme plano. Já para a amostra de superfície rugosa, observamos que existe uma competição entre o efeito de confinamento e a rugosidade em relação à configuração dos vórtices. Apresentamos também, uma tabela que mostra resumidamente os estados estacionários dos vórtices nas três amostras. / The interest to investigate the investigate the behavior of a superconductor has grown in the last few decades. Having in mind to search for better superconducting characteristics, it has been found that bulk samples present characteristics much more different than mesoscopic samples (samples with dimensions of the same order of the same order of the London penetration length and the coherence length). As an example, we can mention the non-formation of an Abrikosov vortex lattice as a consequence of the confinement effect (effect associated with the reduced dimensions of the sample) and also considerable change in the critical field values. In the present work we solved the time dependent Ginzburg-Landau equation (TDGL), in order to make a detailed analysis of the vortex dynamics in mesoscopic thin films. To solve these equations, we have used the link variables method which is gauge invariant. From this, we obtain the Cooper pair density and the magnetization curves. The vortex dynamics was investigated for three different surfaces of the film (concave, convex, and irregular). We have observed that, with respect to the parabolic geometries, there is a considerable difference for the critical fields, as well as for the behavior of the magnetization compared to a flat film. On the other hand, for a sample with an irregular surface, we have seen that there is a competition between the confinement effect and rugosity with respect to vortex configurations. We also present a table which summarizes the vortex stationary states for the three topologies mentioned above.

Comprehensive numerical simulations (reviewed in Dissipative Solitons, Akhmediev and Ankiewicz (Eds.), Springer, Berlin, 2005) of pulse solutions of the cubic--quintic Ginzburg--Landau equation (CGLE), a canonical equation governing the weakly nonlinear behavior of dissipative systems in a wide variety of disciplines, reveal various intriguing and entirely novel classes of solutions. In particular, there are five new classes of pulse or solitary waves solutions, viz. pulsating, creeping, snake, erupting, and chaotic solitons. In contrast to the regular solitary waves investigated in numerous integrable and non--integrable systems over the last three decades, these dissipative solitons are not stationary in time. Rather, they are spatially confined pulse--type structures whose envelopes exhibit complicated temporal dynamics. The numerical simulations also reveal very interesting bifurcations sequences of these pulses as the parameters of the CGLE are varied. In this dissertation, we develop a theoretical framework for these novel classes of solutions. In the first part, we use a traveling wave reduction or a so--called spatial approximation to comprehensively investigate the bifurcations of plane wave and periodic solutions of the CGLE. The primary tools used here are Singularity Theory and Hopf bifurcation theory respectively. Generalized and degenerate Hopf bifurcations have also been considered to track the emergence of global structure such as homoclinic orbits. However, these results appear difficult to correlate to the numerical bifurcation sequences of the dissipative solitons. In the second part of this dissertation, we shift gears to focus on the issues of central interest in the area, i.e., the conditions for the occurrence of the five categories of dissipative solitons, as well the dependence of both their shape and their stability on the various parameters of the CGLE, viz. the nonlinearity, dispersion, linear and nonlinear gain, loss and spectral filtering parameters. Our predictions on the variation of the soliton amplitudes, widths and periods with the CGLE parameters agree with simulation results. For this part, we develop and discuss a variational formalism within which to explore the various classes of dissipative solitons. Given the complex dynamics of the various dissipative solutions, this formulation is, of necessity, significantly generalized over all earlier approaches in several crucial ways. Firstly, the two alternative starting formulations for the Lagrangian are recent and not well explored. Also, after extensive discussions with David Kaup, the trial functions have been generalized considerably over conventional ones to keep the shape relatively simple (and the trial function integrable!) while allowing arbitrary temporal variation of the amplitude, width, position, speed and phase of the pulses. In addition, the resulting Euler--Lagrange equations are treated in a completely novel way. Rather than consider the stable fixed points which correspond to the well--known stationary solitons or plain pulses, we use dynamical systems theory to focus on more complex attractors viz. periodic, quasiperiodic, and chaotic ones. Periodic evolution of the trial function parameters on stable periodic attractors constructed via the method of multiple scales yield solitons whose amplitudes are non--stationary or time dependent. In particular, pulsating, snake (and, less easily, creeping) dissipative solitons may be treated in this manner. Detailed results are presented here for the pulsating solitary waves --- their regimes of occurrence, bifurcations, and the parameter dependences of the amplitudes, widths, and periods agree with simulation results. Finally, we elucidate the Hopf bifurcation mechanism responsible for the various pulsating solitary waves, as well as its absence in Hamiltonian and integrable systems where such structures are absent. Results will be presented for the pulsating and snake soliton cases. Chaotic evolution of the trial function parameters in chaotic regimes identified using dynamical systems analysis would yield chaotic solitary waves. The method also holds promise for detailed modeling of chaotic solitons as well. This overall approach fails only to address the fifth class of dissipative solitons, viz. the exploding or erupting solitons.

We study the fluctuation effects in the seminal cyclic predator-prey model in population dynamics due to Robert May and Warren Leonard both in the zero-dimensional and two-dimensional spatial version. We compute the mean time to extinction of a stable set of coexisting populations driven by large fluctuations. We see that the contribution of large fluctuations to extinction can be captured by a quasi-stationary approximation and the Wentzel–Kramers–Brillouin (WKB) eikonal ansatz. We see that near the Hopf bifurcation, extinctions are fast owing to the flat non-Gaussian distribution whereas away from the bifurcation, extinctions are dominated by large fluctuations of the fat tails of the distribution. We compare our results to Gillespie simulations and a single-species theoretical calculation. In addition, we study the spatio-temporal pattern formation of the stochastic May--Leonard model through the Doi-Peliti coherent state path integral formalism to obtain a coarse-grained Langevin description, i.e. the Complex Ginzburg Landau equation with stochastic noise in one complex field. We see that when one restricts the internal reaction noise to small amplitudes, one can obtain a simple form for the stochastic noise correlations that modify the Complex Ginzburg Landau equation. Finally, we study the effect of coupling a spatially extended May--Leonard model in two dimensions with symmetric predation rates to one with asymmetric rates that is prone to reach extinction. We show that the symmetric region induces otherwise unstable coexistence spiral patterns in the asymmetric May--Leonard lattice. We obtain the stability criterion for this pattern induction as we vary the strength of the extinction inducing asymmetry. This research was sponsored by the Army Research Office and was accomplished under Grant Number W911NF-17-1-0156. / Doctor of Philosophy / In the field of ecology, the cyclic predator-prey patterns in a food web are relevant yet independent to the hierarchical archetype. We study the paradigmatic cyclic May--Leonard model of three species, both analytically and numerically. First, we employ well--established techniques in large-deviation theory to study the extinction of populations induced by large but rare fluctuations. In the zero--dimensional version of the model, we compare the mean time to extinction computed from the theory to numerical simulations. Secondly, we study the stochastic spatial version of the May--Leonard model and show that for values close to the Hopf bifurcation, in the limit of small fluctuations, we can map the coarse-grained description of the model to the Complex Ginsburg Landau Equation, with stochastic noise corrections. Finally, we explore the induction of ecodiversity through spatio-temporal spirals in the asymmetric version of the May--Leonard model, which is otherwise inclined to reach an extinction state. This is accomplished by coupling to a symmetric May-Leonard counterpart on a two-dimensional lattice. The coupled system creates conditions for spiral formation in the asymmetric subsystem, thus precluding extinction.

In this dissertation, I study a nonlinear partial differential equation, the complex Ginzburg-Landau (CGL) equation. I first employed the perturbative field-theoretic renormalization group method to investigate the critical dynamics near the continuous non-equilibrium transition limit in this equation with additive noise. Due to the fact that time translation invariance is broken following a critical quench from a random initial configuration, an independent ``initial-slip'' exponent emerges to describe the crossover temporal window between microscopic time scales and the asymptotic long-time regime. My analytic work shows that to first order in a dimensional expansion with respect to the upper critical dimension, the extracted initial-slip exponent in the complex Ginzburg-Landau equation is identical to that of the equilibrium model A. Subsequently, I studied transient behavior in the CGL through numerical calculations. I developed my own code to numerically solve this partial differential equation on a two-dimensional square lattice with periodic boundary conditions, subject to random initial configurations. Aging phenomena are demonstrated in systems with either focusing and defocusing spiral waves, and the related aging exponents, as well as the auto-correlation exponents, are numerically determined. I also investigated nucleation processes when the system is transiting from a turbulent state to the ``frozen'' state. An extracted finite dimensionless barrier in the deep-quenched case and the exponentially decaying distribution of the nucleation times in the near-transition limit are both suggestive that the dynamical transition observed here is discontinuous. This research is supported by the U. S. Department of Energy, Office of Basic Energy Sciences, Division of Materials Science and Engineering under Award DE-FG02-SC0002308 / Doctor of Philosophy / The complex Ginzburg-Landau equation is one of the most studied nonlinear partial differential equation in the physics community. I study this equation using both analytical and numerical methods. First, I employed the field theory approach to extract the critical initial-slip exponent, which emerges due to the breaking of time translation symmetry and describes the intermediate temporal window between microscopic time scales and the asymptotic long-time regime. I also numerically solved this equation on a two-dimensional square lattice. I studied the scaling behavior in non-equilibrium relaxation processes in situations where defects are interactive but not subject to strong fluctuations. I observed nucleation processes when the system under goes a transition from a strongly fluctuating disordered state to the relatively stable “frozen” state where its dynamics cease. I extracted a finite dimensionless barrier for systems that are quenched deep into the frozen state regime. An exponentially decaying long tail in the nucleation time distribution is found, which suggests a discontinuous transition. This research is supported by the U. S. Department of Energy, Office of Basic Energy Sciences, Division of Materials Science and Engineering under Award DE-FG02-SC0002308.

This thesis concerns the wake control and flow dynamic analysis for a flow around a circular cylinder at different Reynolds numbers using reduced models. The wake control and dynamics in the reduced space were addressed using the sensitivity theory and the adjoint method. In the case of wake control, it was possible to predict the physical parameters of the active and passive controllers on the wake of the main cylinder. On the other hand, in the construction of the reduced space, a new shift mode calculated from a perturbation of the mean flow was proposed using the sensitivity to base flow modifications. The mathematical basis of the reduced space was constructed using a Fourier modal decomposition of the flow enriched by the shift mode. The proposed reduced space made possible the recomposition of the flow and the comparison with the physical parameters calculated in the physical space. Additionally, using the reduced space, it was possible to determine the transition dynamics between the equilibrium point of the Navier Stokes equation and the non-linear saturation state using the Landau coefficients obtained in the reduced model, opening the possibility of solving the flow around a 2D and 3D cylinder with low computational cost. / Esta tese trata sobre o controle de esteira assim como a análise dinâmica do escoamento em torno de um cilindro a diferentes números de Reynolds usando modelos reduzidos. O controle de esteira e a dinâmica no espaço reduzido foram abordados usando a teoria da sensibilidade e o método adjunto. No caso de controle de esteira, foi possível prever os parâmetros físicos dos controladores ativos e passivos no escoamento do cilindro principal. Por outro lado, na construção do espaço reduzido, foi proposto um novo modo de deslocamento (shift mode) calculado a partir de uma perturbação do campo médio usando a sensibilidade às modificações do campo base. A base matemática do espaço reduzido foi construída usando uma decomposição modal de Fourier do escoamento enriquecido pelo modo de deslocamento (shift mode). O espaço reduzido proposto possibilitou a recomposição do escoamento e a comparação com os parâmetros físicos calculados no espaço físico. Além disso, usando o espaço reduzido, foi possível determinar a dinâmica de transição entre o ponto de equilíbrio da equação de Navier Stokes e o estado de saturação não linear usando os coeficientes de Landau obtidos no modelo reduzido, abrindo a possibilidade de resolver o escoamento em torno de um cilindro 2D e 3D com baixo custo computacional

This dissertation focuses on the further development of creating accurate numerical models of complex dynamical systems using the holistic discretisation technique [Roberts, Appl. Num. Model., 37:371-396, 2001]. I extend the application from second to fourth order systems and from only one spatial dimension in all previous work to two dimensions (2D). We see that the holistic technique provides useful and accurate numerical discretisations on coarse grids. We explore techniques to model the evolution of spatial patterns governed by pdes such as the Kuramoto-Sivashinsky equation and the real-valued Ginzburg-Landau equation. We aim towards the simulation of fluid flow and convection in three spatial dimensions. I show that significant steps have been taken in this dissertation towards achieving this aim. Holistic discretisation is based upon centre manifold theory [Carr, Applications of centre manifold theory, 1981] so we are assured that the numerical discretisation accurately models the dynamical system and may be constructed systematically. To apply centre manifold theory the domain is divided into elements and using a homotopy in the coupling parameter, subgrid scale fields are constructed consisting of actual solutions of the governing partial differential equation(pde). These subgrid scale fields interact through the introduction of artificial internal boundary conditions. View the centre manifold (macroscale) as the union of all states of the collection of subgrid fields (microscale) over the physical domain. Here we explore how to extend holistic discretisation to the fourth order Kuramoto-Sivashinsky pde. I show that the holistic models give impressive accuracy for reproducing the steady states and time dependent phenomena of the Kuramoto-Sivashinsky equation on coarse grids. The holistic method based on local dynamics compares favourably to the global methods of approximate inertial manifolds. The excellent performance of the holistic models shown here is strong evidence in support of the holistic discretisation technique. For shear dispersion in a 2D channel a one-dimensional numerical approximation is generated directly from the two-dimensional advection-diffusion dynamics. We find that a low order holistic model contains the shear dispersion term of the Taylor model [Taylor, IMA J. Appl. Math., 225:473-477, 1954]. This new approach does not require the assumption of large x scales, formerly absolutely crucial in deriving the Taylor model. I develop holistic discretisation for two spatial dimensions by applying the technique to the real-valued Ginzburg-Landau equation as a representative example of second order pdes. The techniques will apply quite generally to second order reaction-diffusion equations in 2D. This is the first study implementing holistic discretisation in more than one spatial dimension. The previous applications of holistic discretisation have developed algebraic forms of the subgrid field and its evolution. I develop an algorithm for numerical construction of the subgrid field and its evolution for 1D and 2D pdes and explore various alternatives. This new development greatly extends the class of problems that may be discretised by the holistic technique. This is a vital step for the application of the holistic technique to higher spatial dimensions and towards discretising the Navier-Stokes equations.

This thesis concerns the wake control and flow dynamic analysis for a flow around a circular cylinder at different Reynolds numbers using reduced models. The wake control and dynamics in the reduced space were addressed using the sensitivity theory and the adjoint method. In the case of wake control, it was possible to predict the physical parameters of the active and passive controllers on the wake of the main cylinder. On the other hand, in the construction of the reduced space, a new shift mode calculated from a perturbation of the mean flow was proposed using the sensitivity to base flow modifications. The mathematical basis of the reduced space was constructed using a Fourier modal decomposition of the flow enriched by the shift mode. The proposed reduced space made possible the recomposition of the flow and the comparison with the physical parameters calculated in the physical space. Additionally, using the reduced space, it was possible to determine the transition dynamics between the equilibrium point of the Navier Stokes equation and the non-linear saturation state using the Landau coefficients obtained in the reduced model, opening the possibility of solving the flow around a 2D and 3D cylinder with low computational cost. / Esta tese trata sobre o controle de esteira assim como a análise dinâmica do escoamento em torno de um cilindro a diferentes números de Reynolds usando modelos reduzidos. O controle de esteira e a dinâmica no espaço reduzido foram abordados usando a teoria da sensibilidade e o método adjunto. No caso de controle de esteira, foi possível prever os parâmetros físicos dos controladores ativos e passivos no escoamento do cilindro principal. Por outro lado, na construção do espaço reduzido, foi proposto um novo modo de deslocamento (shift mode) calculado a partir de uma perturbação do campo médio usando a sensibilidade às modificações do campo base. A base matemática do espaço reduzido foi construída usando uma decomposição modal de Fourier do escoamento enriquecido pelo modo de deslocamento (shift mode). O espaço reduzido proposto possibilitou a recomposição do escoamento e a comparação com os parâmetros físicos calculados no espaço físico. Além disso, usando o espaço reduzido, foi possível determinar a dinâmica de transição entre o ponto de equilíbrio da equação de Navier Stokes e o estado de saturação não linear usando os coeficientes de Landau obtidos no modelo reduzido, abrindo a possibilidade de resolver o escoamento em torno de um cilindro 2D e 3D com baixo custo computacional