Accelerating CFD via local POD plus Galerkin projections

Se desarrollan varias técnicas basadas en descomposición ortogonal propia (DOP) local y proyección de tipo Galerkin para acelerar la integración numérica de problemas de evolución, de tipo parabólico, no lineales. Las ideas y métodos que se presentan conllevan un nuevo enfoque para la modelización de tipo DOP, que combina intervalos temporales cortos en que se usa un esquema numérico estándard con otros intervalos temporales en que se utilizan los sistemas de tipo Galerkin que resultan de proyectar las ecuaciones de evolución sobre la variedad lineal generada por los modos DOP, obtenidos a partir de instantáneas calculadas en los intervalos donde actúa el código numérico. La variedad DOP se construye completamente en el primer intervalo, pero solamente se actualiza en los demás intervalos según las dinámicas de la solución, aumentando de este modo la eficiencia del modelo de orden reducido resultante. Además, se aprovechan algunas propiedades asociadas a la dependencia débil de los modos DOP tanto en la variable temporal como en los posibles parámetros de que pueda depender el problema. De esta forma, se aumentan la flexibilidad y la eficiencia computacional del proceso. La aplicación de los métodos resultantes es muy prometedora, tanto en la simulación de transitorios en flujos laminares como en la construcción de diagramas de bifurcación en sistemas dependientes de parámetros. Las ideas y los algoritmos desarrollados en la tesis se ilustran en dos problemas test, la ecuación unidimensional compleja de Ginzburg-Landau y el problema bidimensional no estacionario de la cavidad. Abstract Various ideas and methods involving local proper orthogonal decomposition (POD) and Galerkin projection are presented aiming at accelerating the numerical integration of nonlinear time dependent parabolic problems. The proposed methods come from a new approach to the POD-based model reduction procedures, which combines short runs with a given numerical solver and a reduced order model constructed by expanding the solution of the problem into appropriate POD modes, which span a POD manifold, and Galerkin projecting some evolution equations onto that linear manifold. The POD manifold is completely constructed from the outset, but only updated as time proceeds according to the dynamics, which yields an adaptive and flexible procedure. In addition, some properties concerning the weak dependence of the POD modes on time and possible parameters in the problem are exploited in order to increase the flexibility and efficiency of the low dimensional model computation. Application of the developed techniques to the approximation of transients in laminar fluid flows and the simulation of attractors in bifurcation problems shows very promising results. The test problems considered to illustrate the various ideas and check the performance of the algorithms are the onedimensional complex Ginzburg-Landau equation and the two-dimensional unsteady liddriven cavity problem.

Anomalies in stiffness and damping of a 2D discrete viscoelastic system due to negative stiffness components

The recent development of using negative stiffness inclusions to achieve extreme overall stiffness and mechanical damping of composite materials reveals a new avenue for constructing high performance materials. One of the negative stiffness sources can be obtained from phase transforming materials in the vicinity of their phase transition, as suggested by the Landau theory. To understand the underlying mechanism from a microscopic viewpoint, we theoretically analyze a 2D, nested triangular lattice cell with pre-chosen elements containing negative stiffness to demonstrate anomalies in overall stiffness and damping. Combining with current knowledge from continuum models, based on the composite theory, such as the Voigt, Reuss, and Hashin-Shtrikman model, we further explore the stability of the system with Lyapunov's indirect stability theorem. The evolution of the microstructure in terms of the discrete system is discussed. A potential application of the results presented here is to develop special thin films with unusual in-plane mechanical properties. © 2006 Elsevier B.V. All rights reserved.

Dissipative Raman solitons

A new type of dissipative solitons - dissipative Raman solitons - are revealed on the basis of numerical study of the generalized complex nonlinear Ginzburg-Landau equation. The stimulated Raman scattering significantly affects the energy scalability of the dissipative solitons, causing splitting to multiple pulses. We show, that an appropriate increase of the group-delay dispersion can suppress the multipulsing instability due to formation of the dissipative Raman soliton, which is chirped, has a Stokes-shifted spectrum, and chaotic modulation on its trailing edge. The strong perturbation of a soliton envelope caused by the stimulated Raman scattering confines the energy scalability, preventing the so-called dissipative soliton resonance. We show that in practical implementations, a spectral filter can extend the stability regions of high-energy pulses.

Fast and slowly evolving vector solitons in mode-locked fibre lasers

We report on a new vector model of an erbium-doped fibre laser mode locked with carbon nanotubes. This model goes beyond the limitations of the previously used models based on either coupled nonlinear Schrödinger or Ginzburg-Landau equations. Unlike the previous models, it accounts for the vector nature of the interaction between an optical field and an erbium-doped active medium, slow relaxation dynamics of erbium ions, linear birefringence in a fibre, linear and circular birefringence of a laser cavity caused by in-cavity polarization controller and light-induced anisotropy caused by elliptically polarized pump field. Interplay of aforementioned factors changes coherent coupling of two polarization modes at a long time scale and so results in a new family of vector solitons (VSs) with fast and slowly evolving states of polarization. The observed VSs can be of interest in secure communications, trapping and manipulation of atoms and nanoparticles, control of magnetization in data storage devices and many other areas.

Stabilization of standing waves through time-delay feedback

Standing waves are studied as solutions of a complex Ginzburg-Landau equation subjected to local and global time-delay feedback terms. The onset is described as an instability of the uniform oscillations with respect to spatially periodic perturbations. The solution of the standing wave pattern is given analytically and studied through simulations. © 2013 American Physical Society.

Dissipative vector solitons with fast evolving states of polarization

We report on a new vector model of an erbium doped fiber laser mode locked with carbon nanotubes. This model goes beyond the limitations of the previously used models based on either coupled nonlinear Schrödinger or Ginzburg-Landau equations. It results in a new family of vector solitons with fast evolving states of polarization experimentally observed in our previous papers.

Comments on multiple oscillatory solutions in systems with time-delay feedback

A complex Ginzburg-Landau equation subjected to local and global time-delay feedback terms is considered. In particular, multiple oscillatory solutions and their properties are studied. We present novel results regarding the disappearance of limit cycle solutions, derive analytical criteria for frequency degeneration, amplitude degeneration, frequency extrema. Furthermore, we discuss the influence of the phase shift parameter and show analytically that the stabilization of the steady state and the decay of all oscillations (amplitude death) cannot happen for global feedback only. Finally, we explain the onset of traveling wave patterns close to the regime of amplitude death.

Strongly localized moving discrete dissipative breather-solitons in Kerr nonlinear media supported by intrinsic gain

We investigate the mobility of nonlinear localized modes in a generalized discrete Ginzburg-Landau-type model, describing a one-dimensional waveguide array in an active Kerr medium with intrinsic, saturable gain and damping. It is shown that exponentially localized, traveling discrete dissipative breather-solitons may exist as stable attractors supported only by intrinsic properties of the medium, i.e., in the absence of any external field or symmetry-breaking perturbations. Through an interplay by the gain and damping effects, the moving soliton may overcome the Peierls-Nabarro barrier, present in the corresponding conservative system, by self-induced time-periodic oscillations of its power (norm) and energy (Hamiltonian), yielding exponential decays to zero with different rates in the forward and backward directions. In certain parameter windows, bistability appears between fast modes with small oscillations and slower, large-oscillation modes. The velocities and the oscillation periods are typically related by lattice commensurability and exhibit period-doubling bifurcations to chaotically "walking" modes under parameter variations. If the model is augmented by intersite Kerr nonlinearity, thereby reducing the Peierls-Nabarro barrier of the conservative system, the existence regime for moving solitons increases considerably, and a richer scenario appears including Hopf bifurcations to incommensurately moving solutions and phase-locking intervals. Stable moving breathers also survive in the presence of weak disorder. © 2014 American Physical Society.

Énergétique des vortex et comportement magnétique de la supraconductivité SO(5)

Mémoire numérisé par la Direction des bibliothèques de l'Université de Montréal.

Dissipative parametric modulation instability and pattern formation in nonlinear optical systems

We present the essential features of the dissipative parametric instability, in the universal complex Ginzburg- Landau equation. Dissipative parametric instability is excited through a parametric modulation of frequency dependent losses in a zig-zag fashion in the spectral domain. Such damping is introduced respectively for spectral components in the +ΔF and in the -ΔF region in alternating fashion, where F can represent wavenumber or temporal frequency depending on the applications. Such a spectral modulation can destabilize the homogeneous stationary solution of the system leading to growth of spectral sidebands and to the consequent pattern formation: both stable and unstable patterns in one- and in two-dimensional systems can be excited. The dissipative parametric instability provides an useful and interesting tool for the control of pattern formation in nonlinear optical systems with potentially interesting applications in technological applications, like the design of mode- locked lasers emitting pulse trains with tunable repetition rate; but it could also find realizations in nanophotonics circuits or in dissipative polaritonic Bose-Einstein condensates.

Self-pulsating nonlinear systems via dissipative parametric instability

The recently discovered dissipative parametric instability is presented in the framework of the universal complex Ginzburg-Landau equation. The pattern formation associated with the instability is discussed in connection to the relevant applications in nonlinear photonics especially as a new tool for pulsed lasers design.

Énergétique des vortex et comportement magnétique de la supraconductivité SO(5)

Mémoire numérisé par la Direction des bibliothèques de l'Université de Montréal.

Dimensionality hyper-reduction and machine learning for dynamical systems with varying parameters

Thesis (Ph.D.)--University of Washington, 2016-08

Pattern formation with a conservation law

Pattern formation in systems with a conserved quantity is considered by studying the appropriate amplitude equations. The conservation law leads to a large-scale neutral mode that must be included in the asymptotic analysis for pattern formation near onset. Near a stationary bifurcation, the usual Ginzburg--Landau equation for the amplitude of the pattern is then coupled to an equation for the large-scale mode. These amplitude equations show that for certain parameters all roll-type solutions are unstable. This new instability differs from the Eckhaus instability in that it is amplitude-driven and is supercritical. Beyond the stability boundary, there exist stable stationary solutions in the form of strongly modulated patterns. The envelope of these modulations is calculated in terms of Jacobi elliptic functions and, away from the onset of modulation, is closely approximated by a sech profile. Numerical simulations indicate that as the modulation becomes more pronounced, the envelope broadens. A number of applications are considered, including convection with fixed-flux boundaries and convection in a magnetic field, resulting in new instabilities for these systems.

Random path sampling applied to vector models of magnetism

The recently reported Monte Carlo Random Path Sampling method (RPS) is here improved and its application is expanded to the study of the 2D and 3D Ising and discrete Heisenberg models. The methodology was implemented to allow use in both CPU-based high-performance computing infrastructures (C/MPI) and GPU-based (CUDA) parallel computation, with significant computational performance gains. Convergence is discussed, both in terms of free energy and magnetization dependence on field/temperature. From the calculated magnetization-energy joint density of states, fast calculations of field and temperature dependent thermodynamic properties are performed, including the effects of anisotropy on coercivity, and the magnetocaloric effect. The emergence of first-order magneto-volume transitions in the compressible Ising model is interpreted using the Landau theory of phase transitions. Using metallic Gadolinium as a real-world example, the possibility of using RPS as a tool for computational magnetic materials design is discussed. Experimental magnetic and structural properties of a Gadolinium single crystal are compared to RPS-based calculations using microscopic parameters obtained from Density Functional Theory.