Energy-entropy-momentum time integration methods for coupled smooth dissipative problems

Esta tesis aborda la formulación, análisis e implementación de métodos numéricos de integración temporal para la solución de sistemas disipativos suaves de dimensión finita o infinita de manera que su estructura continua sea conservada. Se entiende por dichos sistemas aquellos que involucran acoplamiento termo-mecánico y/o efectos disipativos internos modelados por variables internas que siguen leyes continuas, de modo que su evolución es considerada suave. La dinámica de estos sistemas está gobernada por las leyes de la termodinámica y simetrías, las cuales constituyen la estructura que se pretende conservar de forma discreta. Para ello, los sistemas disipativos se describen geométricamente mediante estructuras metriplécticas que identifican claramente las partes reversible e irreversible de la evolución del sistema. Así, usando una de estas estructuras conocida por las siglas (en inglés) de GENERIC, la estructura disipativa de los sistemas es identificada del mismo modo que lo es la Hamiltoniana para sistemas conservativos. Con esto, métodos (EEM) con precisión de segundo orden que conservan la energía, producen entropía y conservan los impulsos lineal y angular son formulados mediante el uso del operador derivada discreta introducido para asegurar la conservación de la Hamiltoniana y las simetrías de sistemas conservativos. Siguiendo estas directrices, se formulan dos tipos de métodos EEM basados en el uso de la temperatura o de la entropía como variable de estado termodinámica, lo que presenta importantes implicaciones que se discuten a lo largo de esta tesis. Entre las cuales cabe destacar que las condiciones de contorno de Dirichlet son naturalmente impuestas con la formulación basada en la temperatura. Por último, se validan dichos métodos y se comprueban sus mejores prestaciones en términos de la estabilidad y robustez en comparación con métodos estándar. This dissertation is concerned with the formulation, analysis and implementation of structure-preserving time integration methods for the solution of the initial(-boundary) value problems describing the dynamics of smooth dissipative systems, either finite- or infinite-dimensional ones. Such systems are understood as those involving thermo-mechanical coupling and/or internal dissipative effects modeled by internal state variables considered to be smooth in the sense that their evolutions follow continuos laws. The dynamics of such systems are ruled by the laws of thermodynamics and symmetries which constitutes the structure meant to be preserved in the numerical setting. For that, dissipative systems are geometrically described by metriplectic structures which clearly identify the reversible and irreversible parts of their dynamical evolution. In particular, the framework known by the acronym GENERIC is used to reveal the systems' dissipative structure in the same way as the Hamiltonian is for conserving systems. Given that, energy-preserving, entropy-producing and momentum-preserving (EEM) second-order accurate methods are formulated using the discrete derivative operator that enabled the formulation of Energy-Momentum methods ensuring the preservation of the Hamiltonian and symmetries for conservative systems. Following these guidelines, two kind of EEM methods are formulated in terms of entropy and temperature as a thermodynamical state variable, involving important implications discussed throughout the dissertation. Remarkably, the formulation in temperature becomes central to accommodate Dirichlet boundary conditions. EEM methods are finally validated and proved to exhibit enhanced numerical stability and robustness properties compared to standard ones.

On uniqueness, transfer and combination of acoustic sources in one-dimensional systems

The use of electroacoustic analogies suggests that a source of acoustical energy (such as an engine, compressor, blower, turbine, loudspeaker, etc.) can be characterized by an acoustic source pressure ps and internal source impedance Zs, analogous to the open-circuit voltage and internal impedance of an electrical source. The present paper shows analytically that the source characteristics evaluated by means of the indirect methods are independent of the loads selected; that is, the evaluated values of ps and Zs are unique, and that the results of the different methods (including the direct method) are identical. In addition, general relations have been derived here for the transfer of source characteristics from one station to another station across one or more acoustical elements, and also for combining several sources into a single equivalent source. Finally, all the conclusions are extended to the case of a uniformly moving medium, incorporating the convective as well as dissipative effects of the mean flow.

Authors reply to "On uniqueness, transfer and combination of acoustic sources in one-dimensional systems"

The use of electroacoustic analogies suggests that a source of acoustical energy (such as an engine, compressor, blower, turbine, loudspeaker, etc.) can be characterized by an acoustic source pressure ps and internal source impedance Zs, analogous to the open-circuit voltage and internal impedance of an electrical source. The present paper shows analytically that the source characteristics evaluated by means of the indirect methods are independent of the loads selected; that is, the evaluated values of ps and Zs are unique, and that the results of the different methods (including the direct method) are identical. In addition, general relations have been derived here for the transfer of source characteristics from one station to another station across one or more acoustical elements, and also for combining several sources into a single equivalent source. Finally, all the conclusions are extended to the case of a uniformly moving medium, incorporating the convective as well as dissipative effects of the mean flow.

Rotor - Stator broadband noise prediction

In this paper a semi analytic model for rotor - stator broadband noise is presented. The work can be split into two sections. The first examines the distortion of the rotor wake in mean swirling flow, downstream of the fan. Previous work by Cooper and Peake4 is extended to include dissipative effects. In the second section we consider the interaction of this gust with the downstream stator row. We examine the way in which an unsteady pressure field is generated by the interaction of this wake flow with the stator blades and obtain estimates for the radiated noise. A new method is presented to extend the well known LINSUB code to the third dimension to capture the effect of the spanwise wavenumber and stator lean and sweep. Copyright © 2008 by Adrian Lloyd and Nigel Peake.

Decoherence of two Josephson charge qubits coupled via sigma(x)sigma(x)-type coupling and exposed to sigma(x) noise

Decoherence properties of two Josephson charge qubits coupled via the sigma(x)sigma(x) type are investigated. Considering the special structure of this new design, the dissipative effects arising from the circuit impedance providing the fluxes for the qubits' superconducting quantum interference device loops coupled to the sigma(x) qubit variables are considered. The results show that the overall decoherence effects are significantly strong in this qubit design. It is found that the dissipative effects are stronger in the case of coupling to two uncorrelated baths than are found in the case of one common bath.

Ressonâncias planetárias spin-órbita com perturbação de terceiro corpo

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

Hyperbolic-like estimates for higher order equations

The main goal of this paper is to derive long time estimates of the energy for the higher order hyperbolic equations with time-dependent coefficients. in particular, we estimate the energy in the hyperbolic zone of the extended phase space by means of a function f (t) which depends on the principal part and on the coefficients of the terms of order m - 1. Then we look for sufficient conditions that guarantee the same energy estimate from above in all the extended phase space. We call this class of estimates hyperbolic-like since the energy behavior is deeply depending on the hyperbolic structure of the equation. In some cases, these estimates produce a dissipative effect on the energy. (C) 2012 Elsevier Inc. All rights reserved.

Tratamiento consistente de restricciones e impacto en dinámica de sistemas multicuerpo

Esta tesis está enmarcada en el estudio de diferentes procedimientos numéricos para resolver la dinámica de un sistema multicuerpo sometido a restricciones e impacto, que puede estar compuesto por sólidos rígidos y deformables conectados entre sí por diversos tipos de uniones. Dentro de los métodos numéricos analizados se presta un especial interés a los métodos consistentes, los cuales tienen por objetivo que la energía calculada en cada paso de tiempo, para un sistema mecánico, tenga una evolución coherente con el comportamiento teórico de la energía. En otras palabras, un método consistente mantiene constante la energía total en un problema conservativo, y en presencia de fuerzas disipativas proporciona un decremento positivo de la energía total. En esta línea se desarrolla un algoritmo numérico consistente con la energía total para resolver las ecuaciones de la dinámica de un sistema multicuerpo. Como parte de este algoritmo se formulan energéticamente consistentes las restricciones y el contacto empleando multiplicadores de Lagrange, penalización y Lagrange aumentado. Se propone también un método para el contacto con sólidos rígidos representados mediante superficies implícitas, basado en una restricción regularizada que se adaptada adecuadamente para el cumplimiento exacto de la restricción de contacto y para ser consistente con la conservación de la energía total. En este contexto se estudian dos enfoques: uno para el contacto elástico puro (sin deformación) formulado con penalización y Lagrange aumentado; y otro basado en un modelo constitutivo para el contacto con penetración. En el segundo enfoque se usa un potencial de penalización que, en ausencia de componentes disipativas, restaura la energía almacenada en el contacto y disipa energía de forma consistente con el modelo continuo cuando las componentes de amortiguamiento y fricción son consideradas. This thesis focuses on the study of several numerical procedures used to solve the dynamics of a multibody system subjected to constraints and impact. The system may be composed by rigid and deformable bodies connected by different types of joints. Within this framework, special attention is paid to consistent methods, which preserve the theoretical behavior of the energy at each time step. In other words, a consistent method keeps the total energy constant in a conservative problem, and provides a positive decrease in the total energy when dissipative forces are present. A numerical algorithm has been developed for solving the dynamical equations of multibody systems, which is energetically consistent. Energetic consistency in contacts and constraints is formulated using Lagrange multipliers, penalty and augmented Lagrange methods. A contact methodology is proposed for rigid bodies with a boundary represented by implicit surfaces. The method is based on a suitable regularized constraint formulation, adapted both to fulfill exactly the contact constraint, and to be consistent with the conservation of the total energy. In this context two different approaches are studied: the first applied to pure elastic contact (without deformation), formulated with penalty and augmented Lagrange; and a second one based on a constitutive model for contact with penetration. In this second approach, a penalty potential is used in the constitutive model, that restores the energy stored in the contact when no dissipative effects are present. On the other hand, the energy is dissipated consistently with the continuous model when friction and damping are considered.

Nondiffracting Bessel plasmons

We report on the existence of nondiffracting Bessel surface plasmon polaritons (SPPs), advancing at either superluminal or subluminal phase velocities. These wave fields feature deep subwavelength FWHM, but are supported by high-order homogeneous SPPs of a metal/dielectric (MD) superlattice. The beam axis can be relocated to any MD interface, by interfering multiple converging SPPs with controlled phase matching. Dissipative effects in metals lead to a diffraction-free regime that is limited by the energy attenuation length. However, the ultra-localization of the diffracted wave field might still be maintained by more than one order of magnitude.

Many-body quantum dynamics of polarization squeezing in optical fibers

We report new experiments that test quantum dynamical predictions of polarization squeezing for ultrashort photonic pulses in a birefringent fiber, including all relevant dissipative effects. This exponentially complex many-body problem is solved by means of a stochastic phase-space method. The squeezing is calculated and compared to experimental data, resulting in excellent quantitative agreement. From the simulations, we identify the physical limits to quantum noise reduction in optical fibers. The research represents a significant experimental test of first-principles time-domain quantum dynamics in a one-dimensional interacting Bose gas coupled to dissipative reservoirs.

Temperature effects on a network of dissipative quantum harmonic oscillators: collective damping and dispersion processes

In this paper we extend the results presented in (de Ponte, Mizrahi and Moussa 2007 Phys. Rev. A 76 032101) to treat quantitatively the effects of reservoirs at finite temperature in a bosonic dissipative network: a chain of coupled harmonic oscillators whatever its topology, i.e., whichever the way the oscillators are coupled together, the strength of their couplings and their natural frequencies. Starting with the case where distinct reservoirs are considered, each one coupled to a corresponding oscillator, we also analyze the case where a common reservoir is assigned to the whole network. Master equations are derived for both situations and both regimes of weak and strong coupling strengths between the network oscillators. Solutions of these master equations are presented through the normal ordered characteristic function. These solutions are shown to be significantly involved when temperature effects are considered, making difficult the analysis of collective decoherence and dispersion in dissipative bosonic networks. To circumvent these difficulties, we turn to the Wigner distribution function which enables us to present a technique to estimate the decoherence time of network states. Our technique proceeds by computing separately the effects of dispersion and the attenuation of the interference terms of the Wigner function. A detailed analysis of the dispersion mechanism is also presented through the evolution of the Wigner function. The interesting collective dispersion effects are discussed and applied to the analysis of decoherence of a class of network states. Finally, the entropy and the entanglement of a pure bipartite system are discussed.

Nanoparticle induced miscibility in LCST polymer blends: critically assessing the enthalpic and entropic effects

The use of copolymer and polymer blends widened the possibility of creating materials with multilayered architectures. Hierarchical polymer systems with a wide array of micro and nanostructures are generated by thermally induced phase separation (TIPS) in partially miscible polymer blends. Various parameters like the interaction between the polymers, concentration, solvent/non-solvent ratio, and quenching temperature have to be optimized to obtain these micro/nanophase structures. Alternatively, the addition of nanoparticles is another strategy to design materials with desired hetero-phase structures. The dynamics of the polymer nanocomposite depends on the statistical ordering of polymers around the nanoparticle, which is dependent on the shape of the nanoparticle. The entropic loss due to deformation of polymer chains, like the repulsive interactions due to coiling and the attractive interactions in the case of swelling has been highlighted in this perspective article. The dissipative particle dynamics has been discussed and is correlated with the molecular dynamics simulation in the case of polymer blends. The Cahn Hillard Cook model on variedly shaped immobile fillers has shown difference in the propagation of the composition wave. The nanoparticle shape has a contributing effect on the polymer particle interaction, which can change the miscibility window in the case of these phase separating polymer blends. Quantitative information on the effect of spherical particles on the demixing temperature is well established and further modified to explain the percolation of rod shaped particles in the polymer blends. These models correlate well with the experimental observations in context to the dynamics induced by the nanoparticle in the demixing behavior of the polymer blend. The miscibility of the LCST polymer blend depends on the enthalpic factors like the specific interaction between the components, and the solubility product and the entropic losses occurring due to the formation of any favorable interactions. Hence, it is essential to assess the entropic and enthalpic interactions induced by the nanoparticles independently. The addition of nanoparticles creates heterogeneity in the polymer phase it is localized. This can be observed as an alteration in the relaxation behavior of the polymer. This changes the demixing behavior and the interaction parameter between the polymers. The compositional changes induced due to the incorporation of nanoparticles are also attributed as a reason for the altered demixing temperature. The particle shape anisotropy causes a direction dependent depletion, which changes the phase behavior of the blend. The polymer-grafted nanoparticles with varying grafting density show tremendous variation in the miscibility of the blend. The stretching of the polymer chains grafted on the nanoparticles causes an entropy penalty in the polymer blend. A comparative study on the different shaped particles is not available up to date for understanding these aspects. Hence, we have juxtaposed the various computational studies on nanoparticle dynamics, the shape effect of NPs on homopolymers and also the cases of various polymer blends without nanoparticles to sketch a complete picture on the effect of various particles on the miscibility of LCST blends.

Rate and Microstructure Effects on the Dynamics of Carbon Nanotube Foams

Soft hierarchical materials often present unique functional properties that are sensitive to the geometry and organization of their micro- and nano-structural features across different lengthscales. Carbon Nanotube (CNT) foams are hierarchical materials with fibrous morphology that are known for their remarkable physical, chemical and electrical properties. Their complex microstructure has led them to exhibit intriguing mechanical responses at different length-scales and in different loading regimes. Even though these materials have been studied for mechanical behavior over the past few years, their response at high-rate finite deformations and the influence of their microstructure on bulk mechanical behavior and energy dissipative characteristics remain elusive.

In this dissertation, we study the response of aligned CNT foams at the high strain-rate regime of 10² - 10⁴ s^-1. We investigate their bulk dynamic response and the fundamental deformation mechanisms at different lengthscales, and correlate them to the microstructural characteristics of the foams. We develop an experimental platform, with which to study the mechanics of CNT foams in high-rate deformations, that includes direct measurements of the strain and transmitted forces, and allows for a full field visualization of the sample’s deformation through high-speed microscopy.

We synthesize various CNT foams (e.g., vertically aligned CNT (VACNT) foams, helical CNT foams, micro-architectured VACNT foams and VACNT foams with microscale heterogeneities) and show that the bulk functional properties of these materials are highly tunable either by tailoring their microstructure during synthesis or by designing micro-architectures that exploit the principles of structural mechanics. We also develop numerical models to describe the bulk dynamic response using multiscale mass-spring models and identify the mechanical properties at length scales that are smaller than the sample height.

The ability to control the geometry of microstructural features, and their local interactions, allows the creation of novel hierarchical materials with desired functional properties. The fundamental understanding provided by this work on the key structure-function relations that govern the bulk response of CNT foams can be extended to other fibrous, soft and hierarchical materials. The findings can be used to design materials with tailored properties for different engineering applications, like vibration damping, impact mitigation and packaging.

Size effects in the conical indentation of an elasto-plastic solid

The size effect in conical indentation of an elasto-plastic solid is predicted via the Fleck and Willis formulation of strain gradient plasticity (Fleck, N.A. and Willis, J.R., 2009, A mathematical basis for strain gradient plasticity theory. Part II: tensorial plastic multiplier, J. Mech. Phys. Solids, 57, 1045-1057). The rate-dependent formulation is implemented numerically and the full-field indentation problem is analyzed via finite element calculations, for both ideally plastic behavior and dissipative hardening. The isotropic strain-gradient theory involves three material length scales, and the relative significance of these length scales upon the degree of size effect is assessed. Indentation maps are generated to summarize the sensitivity of indentation hardness to indent size, indenter geometry and material properties (such as yield strain and strain hardening index). The finite element model is also used to evaluate the pertinence of the Johnson cavity expansion model and of the Nix-Gao model, which have been extensively used to predict size effects in indentation hardness. © 2012 Elsevier Ltd.

Dissipative scheme to approach the boundary of two-qubit entangled mixed states

We discuss the generation of states close to the boundary family of maximally entangled mixed states as defined by the use of concurrence and linear entropy. The coupling of two qubits to a dissipation-affected bosonic mode is able to produce a bipartite state having, for all practical purposes, the entanglement and mixedness properties of one of such boundary states. We thoroughly study the effects that thermal and squeezed characters of the bosonic mode have in such a process and we discuss tolerance to qubit phase-damping mechanisms. The nondemanding nature of the scheme makes it realizable in a matter-light-based physical setup, which we address in some details.