Nonequilibrium self-energies, Ng approach, and heat current of a nanodevice for small bias voltage and temperature

NASA Astrophysics Data System (ADS)

Aligia, A. A.

2014-03-01

Using nonequilibrium renormalized perturbation theory to second order in the renormalized Coulomb repulsion, we calculate the lesser Σ< and and greater Σ> self-energies of the impurity Anderson model, which describes the current through a quantum dot, in the general asymmetric case. While in general a numerical integration is required to evaluate the perturbative result, we derive an analytical approximation for small frequency ω, bias voltage V, and temperature T, which is exact to total second order in these quantities. The approximation is valid when the corresponding energies ℏω, eV, and kBT are small compared to kBTK, where TK is the Kondo temperature. The result of the numerical integration is compared with the analytical one and with Ng approximation, in which Σ< and Σ> are assumed proportional to the retarded self-energy Σr times an average Fermi function. While it fails at T =0 for ℏ |ω|≲eV, we find that the Ng approximation is excellent for kBT>eV/2 and improves for asymmetric coupling to the leads. Even at T =0, the effect of the Ng approximation on the total occupation at the dot is very small. The dependence on ω and V are discussed in comparison with a Ward identity that is fulfilled by the three approaches. We also calculate the heat currents between the dot and any of the leads at finite bias voltage. One of the heat currents changes sign with the applied bias voltage at finite temperature.

Predictive Understanding of Mountainous Watershed Hydro-Biogeochemical Function and Response to Perturbations

NASA Astrophysics Data System (ADS)

Hubbard, S. S.; Williams, K. H.; Agarwal, D.; Banfield, J. F.; Beller, H. R.; Bouskill, N.; Brodie, E.; Maxwell, R. M.; Nico, P. S.; Steefel, C. I.; Steltzer, H.; Tokunaga, T. K.; Wainwright, H. M.; Dwivedi, D.; Newcomer, M. E.

2017-12-01

Recognizing the societal importance, vulnerability and complexity of mountainous watersheds, the `Watershed Function' project is developing a predictive understanding of how mountainous watersheds retain and release downgradient water, nutrients, carbon, and metals. In particular, the project is exploring how early snowmelt, drought, floods and other disturbances will influence mountainous watershed dynamics at seasonal to decadal timescales. Located in the 300km2 East River headwater catchment of the Upper Colorado River Basin, the project is guided by several constructs. First, the project considers the integrated role of surface and subsurface flow and biogeochemical reactions - from bedrock to the top of the vegetative canopy, from terrestrial through aquatic compartments, and from summit to receiving waters. The project takes a system-of-systems perspective, focused on developing new methods to quantify the cumulative watershed hydrobiogeochemical response to perturbations based on information from select subsystems within the watershed, each having distinct vegetation-subsurface biogeochemical-hydrological characteristics. A `scale-adaptive' modeling capability, in development using adaptive mesh refinement methods, serves as the organizing framework for the SFA. The scale-adaptive approach is intended to permit simulation of system-within-systems behavior - and aggregation of that behavior - from genome through watershed scales. This presentation will describe several early project discoveries and advances made using experimental, observational and numerical approaches. Among others, examples may include:quantiying how seasonal hydrological perturbations drive biogeochemical responses across critical zone compartments, with a focus on N and C transformations; metagenomic documentation of the spatial variability in floodplain meander microbial ecology; 3D reactive transport simulations of couped hydrological-biogeochemical behavior in the hyporheic zone; and new characterization and inversion approaches to quantify co-variability between above and below ground dynamics and the susceptibility of vegetation to drought stress. More information is provided at: Watershed.lbl.gov

Singular perturbation of smoothly evolving Hele-Shaw solutions

DOE Office of Scientific and Technical Information (OSTI.GOV)

Siegel, M.; Tanveer, S.

1996-01-01

We present analytical scaling results, confirmed by accurate numerics, to show that there exists a class of smoothly evolving zero surface tension solutions to the Hele-Shaw problem that are significantly perturbed by an arbitrarily small amount of surface tension in order one time. {copyright} {ital 1996 The American Physical Society.}

Hamilton-Jacobi modelling of relative motion for formation flying.

PubMed

Kolemen, Egemen; Kasdin, N Jeremy; Gurfil, Pini

2005-12-01

A precise analytic model for the relative motion of a group of satellites in slightly elliptic orbits is introduced. With this aim, we describe the relative motion of an object relative to a circular or slightly elliptic reference orbit in the rotating Hill frame via a low-order Hamiltonian, and solve the Hamilton-Jacobi equation. This results in a first-order solution to the relative motion identical to the Clohessy-Wiltshire approach; here, however, rather than using initial conditions as our constants of the motion, we utilize the canonical momenta and coordinates. This allows us to treat perturbations in an identical manner, as in the classical Delaunay formulation of the two-body problem. A precise analytical model for the base orbit is chosen with the included effect of zonal harmonics (J(2), J(3), J(4)). A Hamiltonian describing the real relative motion is formed and by differing this from the nominal Hamiltonian, the perturbing Hamiltonian is obtained. Using the Hamilton equations, the variational equations for the new constants are found. In a manner analogous to the center manifold reduction procedure, the non-periodic part of the motion is canceled through a magnitude analysis leading to simple boundedness conditions that cancel the drift terms due to the higher order perturbations. Using this condition, the variational equations are integrated to give periodic solutions that closely approximate the results from numerical integration (1 mm/per orbit for higher order and eccentricity perturbations and 30 cm/per orbit for zonal perturbations). This procedure provides a compact and insightful analytic description of the resulting relative motion.

Validation and application of auxiliary density perturbation theory and non-iterative approximation to coupled-perturbed Kohn-Sham approach for calculation of dipole-quadrupole polarizability

NASA Astrophysics Data System (ADS)

Shedge, Sapana V.; Pal, Sourav; Köster, Andreas M.

2011-07-01

Recently, two non-iterative approaches have been proposed to calculate response properties within density functional theory (DFT). These approaches are auxiliary density perturbation theory (ADPT) and the non-iterative approach to the coupled-perturbed Kohn-Sham (NIA-CPKS) method. Though both methods are non-iterative, they use different techniques to obtain the perturbed Kohn-Sham matrix. In this Letter, for the first time, both of these two independent methods have been used for the calculation of dipole-quadrupole polarizabilities. To validate these methods, three tetrahedral molecules viz., P4,CH4 and adamantane (C10H16) have been used as examples. The comparison with MP2 and CCSD proves the reliability of the methodology.

Stability under scalar perturbations and quasinormal modes of 4D Einstein-Born-Infeld dilaton spacetime: exact spectrum

NASA Astrophysics Data System (ADS)

Destounis, Kyriakos; Panotopoulos, Grigoris; Rincón, Ángel

2018-02-01

We study the stability under scalar perturbations, and we compute the quasinormal modes of the Einstein-Born-Infeld dilaton spacetime in 1+3 dimensions. Solving the full radial equation in terms of hypergeometric functions, we provide an exact analytical expression for the spectrum. We find that the frequencies are purely imaginary, and we confirm our results by computing them numerically. Although the scalar field that perturbs the black hole is electrically neutral, an instability similar to that seen in charged scalar perturbations of the Reissner-Nordström black hole is observed.

Application of ab initio many-body perturbation theory with Gaussian basis sets to the singlet and triplet excitations of organic molecules

NASA Astrophysics Data System (ADS)

Hamed, Samia; Rangel, Tonatiuh; Bruneval, Fabien; Neaton, Jeffrey B.

Quantitative understanding of charged and neutral excitations of organic molecules is critical in diverse areas of study that include astrophysics and the development of energy technologies that are clean and efficient. The recent use of local basis sets with ab initio many-body perturbation theory in the GW approximation and the Bethe-Saltpeter equation approach (BSE), methods traditionally applied to periodic condensed phases with a plane-wave basis, has opened the door to detailed study of such excitations for molecules, as well as accurate numerical benchmarks. Here, through a series of systematic benchmarks with a Gaussian basis, we report on the extent to which the predictive power and utility of this approach depend critically on interdependent underlying approximations and choices for molecules, including the mean-field starting point (eg optimally-tuned range separated hybrids, pure DFT functionals, and untuned hybrids), the GW scheme, and the Tamm Dancoff approximation. We demonstrate the effects of these choices in the context of Thiels' set while drawing analogies to linear-response time-dependent DFT and making comparisons to best theoretical estimates from higher-order wavefunction-based theories.

A New Homotopy Perturbation Scheme for Solving Singular Boundary Value Problems Arising in Various Physical Models

NASA Astrophysics Data System (ADS)

Roul, Pradip; Warbhe, Ujwal

2017-08-01

The classical homotopy perturbation method proposed by J. H. He, Comput. Methods Appl. Mech. Eng. 178, 257 (1999) is useful for obtaining the approximate solutions for a wide class of nonlinear problems in terms of series with easily calculable components. However, in some cases, it has been found that this method results in slowly convergent series. To overcome the shortcoming, we present a new reliable algorithm called the domain decomposition homotopy perturbation method (DDHPM) to solve a class of singular two-point boundary value problems with Neumann and Robin-type boundary conditions arising in various physical models. Five numerical examples are presented to demonstrate the accuracy and applicability of our method, including thermal explosion, oxygen-diffusion in a spherical cell and heat conduction through a solid with heat generation. A comparison is made between the proposed technique and other existing seminumerical or numerical techniques. Numerical results reveal that only two or three iterations lead to high accuracy of the solution and this newly improved technique introduces a powerful improvement for solving nonlinear singular boundary value problems (SBVPs).

First-Principles Lattice Dynamics Method for Strongly Anharmonic Crystals

NASA Astrophysics Data System (ADS)

Tadano, Terumasa; Tsuneyuki, Shinji

2018-04-01

We review our recent development of a first-principles lattice dynamics method that can treat anharmonic effects nonperturbatively. The method is based on the self-consistent phonon theory, and temperature-dependent phonon frequencies can be calculated efficiently by incorporating recent numerical techniques to estimate anharmonic force constants. The validity of our approach is demonstrated through applications to cubic strontium titanate, where overall good agreement with experimental data is obtained for phonon frequencies and lattice thermal conductivity. We also show the feasibility of highly accurate calculations based on a hybrid exchange-correlation functional within the present framework. Our method provides a new way of studying lattice dynamics in severely anharmonic materials where the standard harmonic approximation and the perturbative approach break down.

A Well-Balanced Path-Integral f-Wave Method for Hyperbolic Problems with Source Terms

PubMed Central

2014-01-01

Systems of hyperbolic partial differential equations with source terms (balance laws) arise in many applications where it is important to compute accurate time-dependent solutions modeling small perturbations of equilibrium solutions in which the source terms balance the hyperbolic part. The f-wave version of the wave-propagation algorithm is one approach, but requires the use of a particular averaged value of the source terms at each cell interface in order to be “well balanced” and exactly maintain steady states. A general approach to choosing this average is developed using the theory of path conservative methods. A scalar advection equation with a decay or growth term is introduced as a model problem for numerical experiments. PMID:24563581

A comparison of Lorentz, planetary gravitational, and satellite gravitational resonances

NASA Technical Reports Server (NTRS)

Hamilton, Douglas P.

1994-01-01

We consider a charged dust grain whose orbital motion is dominated by a planet's point-source gravity, but perturbed by higher-order terms in the planet's gravity field as well as by the Lorentz force arising from an asymmetric planetary magnetic field. Perturbations to Keplerian orbits due to a nonspherical gravity field are expressed in the traditional way: in terms of a disturbing function which can be expanded in a series of spherical harmonics (W. M. Kaula, 1966). In order to calculate the electromagnetic perturbation, we first write the Lorentz force in terms of the orbital elements and then substitute it into Gauss' perturbation equations. We use our result to derive strengths of Lorentz resonances and elucidate their properties. In particular, we compare Lorentz resonances to two types of gravitational resonances: those arising from periodic tugs of a satellite and those due to the attraction of an arbitrarily shaped planet. We find that Lorentz resonances share numerous properties with their gravitational counterparts and show, using simple physical arguments, that several of these patterns are fundamental, applying not only to our expansions, but to all quantities expressed in terms of orbital elements. Some of these patterns have been previously called 'd'Alembert rules' for satellite resonances. Other similarities arise because, to first-order in the perturbing force, the three problems share an integral of the motion. Yet there are also differences; for example, first-order inclination resonances exist for perturbations arising from planetary gravity and from the Lorentz force, but not for those due to an orbiting satellite. Finally, we provide a heuristic treatment of a particle's orbital evolution under the influence of drag and resonant forces. Particles brought into mean-motion resonances experience either trapping or resonant 'jumps,' depending on the direction from which the resonance is approached. We show that this behavior does not depend on the details of the perturbing force but rather is fundamental to all mean-motion resonances.

Modeling of multi-band drift in nanowires using a full band Monte Carlo simulation

NASA Astrophysics Data System (ADS)

Hathwar, Raghuraj; Saraniti, Marco; Goodnick, Stephen M.

2016-07-01

We report on a new numerical approach for multi-band drift within the context of full band Monte Carlo (FBMC) simulation and apply this to Si and InAs nanowires. The approach is based on the solution of the Krieger and Iafrate (KI) equations [J. B. Krieger and G. J. Iafrate, Phys. Rev. B 33, 5494 (1986)], which gives the probability of carriers undergoing interband transitions subject to an applied electric field. The KI equations are based on the solution of the time-dependent Schrödinger equation, and previous solutions of these equations have used Runge-Kutta (RK) methods to numerically solve the KI equations. This approach made the solution of the KI equations numerically expensive and was therefore only applied to a small part of the Brillouin zone (BZ). Here we discuss an alternate approach to the solution of the KI equations using the Magnus expansion (also known as "exponential perturbation theory"). This method is more accurate than the RK method as the solution lies on the exponential map and shares important qualitative properties with the exact solution such as the preservation of the unitary character of the time evolution operator. The solution of the KI equations is then incorporated through a modified FBMC free-flight drift routine and applied throughout the nanowire BZ. The importance of the multi-band drift model is then demonstrated for the case of Si and InAs nanowires by simulating a uniform field FBMC and analyzing the average carrier energies and carrier populations under high electric fields. Numerical simulations show that the average energy of the carriers under high electric field is significantly higher when multi-band drift is taken into consideration, due to the interband transitions allowing carriers to achieve higher energies.

Numerical simulation and parametric analysis of selective laser melting process of AlSi10Mg powder

NASA Astrophysics Data System (ADS)

Pei, Wei; Zhengying, Wei; Zhen, Chen; Junfeng, Li; Shuzhe, Zhang; Jun, Du

2017-08-01

A three-dimensional numerical model was developed to investigate effects of laser scanning speed, laser power, and hatch spacing on the thermodynamic behaviors of the molten pool during selective laser melting of AlSi10Mg powder. A randomly distributed packed powder bed was achieved using discrete element method (DEM). The powder bed can be treated as a porous media with interconnected voids in the simulation. A good agreement between numerical results and experimental results establish the validity of adopted method. The numerical results show that the Marangoni flow within the molten pool was significantly affected by the processing parameters. An intense Marangoni flow leads to a perturbation within the molten pool. In addition, a relatively high scanning speed tends to cause melt instability. The perturbation or the instability within the molten pool results in the formation of pores during SLM, which have a direct influence on the densification level.

Nonlinear ion-acoustic cnoidal waves in a dense relativistic degenerate magnetoplasma.

PubMed

El-Shamy, E F

2015-03-01

The complex pattern and propagation characteristics of nonlinear periodic ion-acoustic waves, namely, ion-acoustic cnoidal waves, in a dense relativistic degenerate magnetoplasma consisting of relativistic degenerate electrons and nondegenerate cold ions are investigated. By means of the reductive perturbation method and appropriate boundary conditions for nonlinear periodic waves, a nonlinear modified Korteweg-de Vries (KdV) equation is derived and its cnoidal wave is analyzed. The various solutions of nonlinear ion-acoustic cnoidal and solitary waves are presented numerically with the Sagdeev potential approach. The analytical solution and numerical simulation of nonlinear ion-acoustic cnoidal waves of the nonlinear modified KdV equation are studied. Clearly, it is found that the features (amplitude and width) of nonlinear ion-acoustic cnoidal waves are proportional to plasma number density, ion cyclotron frequency, and direction cosines. The numerical results are applied to high density astrophysical situations, such as in superdense white dwarfs. This research will be helpful in understanding the properties of compact astrophysical objects containing cold ions with relativistic degenerate electrons.

Nonlinear (time domain) and linearized (time and frequency domain) solutions to the compressible Euler equations in conservation law form

NASA Technical Reports Server (NTRS)

Sreenivas, Kidambi; Whitfield, David L.

1995-01-01

Two linearized solvers (time and frequency domain) based on a high resolution numerical scheme are presented. The basic approach is to linearize the flux vector by expressing it as a sum of a mean and a perturbation. This allows the governing equations to be maintained in conservation law form. A key difference between the time and frequency domain computations is that the frequency domain computations require only one grid block irrespective of the interblade phase angle for which the flow is being computed. As a result of this and due to the fact that the governing equations for this case are steady, frequency domain computations are substantially faster than the corresponding time domain computations. The linearized equations are used to compute flows in turbomachinery blade rows (cascades) arising due to blade vibrations. Numerical solutions are compared to linear theory (where available) and to numerical solutions of the nonlinear Euler equations.

Exact calculation of the time convolutionless master equation generator: Application to the nonequilibrium resonant level model

DOE Office of Scientific and Technical Information (OSTI.GOV)

Kidon, Lyran; The Sackler Center for Computational Molecular and Materials Science, Tel Aviv University, Tel Aviv 69978; Wilner, Eli Y.

2015-12-21

The generalized quantum master equation provides a powerful tool to describe the dynamics in quantum impurity models driven away from equilibrium. Two complementary approaches, one based on Nakajima–Zwanzig–Mori time-convolution (TC) and the other on the Tokuyama–Mori time-convolutionless (TCL) formulations provide a starting point to describe the time-evolution of the reduced density matrix. A key in both approaches is to obtain the so called “memory kernel” or “generator,” going beyond second or fourth order perturbation techniques. While numerically converged techniques are available for the TC memory kernel, the canonical approach to obtain the TCL generator is based on inverting a super-operatormore » in the full Hilbert space, which is difficult to perform and thus, nearly all applications of the TCL approach rely on a perturbative scheme of some sort. Here, the TCL generator is expressed using a reduced system propagator which can be obtained from system observables alone and requires the calculation of super-operators and their inverse in the reduced Hilbert space rather than the full one. This makes the formulation amenable to quantum impurity solvers or to diagrammatic techniques, such as the nonequilibrium Green’s function. We implement the TCL approach for the resonant level model driven away from equilibrium and compare the time scales for the decay of the generator with that of the memory kernel in the TC approach. Furthermore, the effects of temperature, source-drain bias, and gate potential on the TCL/TC generators are discussed.« less

Non-perturbative String Theory from Water Waves

DOE Office of Scientific and Technical Information (OSTI.GOV)

Iyer, Ramakrishnan; Johnson, Clifford V.; /Southern California U.

2012-06-14

We use a combination of a 't Hooft limit and numerical methods to find non-perturbative solutions of exactly solvable string theories, showing that perturbative solutions in different asymptotic regimes are connected by smooth interpolating functions. Our earlier perturbative work showed that a large class of minimal string theories arise as special limits of a Painleve IV hierarchy of string equations that can be derived by a similarity reduction of the dispersive water wave hierarchy of differential equations. The hierarchy of string equations contains new perturbative solutions, some of which were conjectured to be the type IIA and IIB string theoriesmore » coupled to (4, 4k ? 2) superconformal minimal models of type (A, D). Our present paper shows that these new theories have smooth non-perturbative extensions. We also find evidence for putative new string theories that were not apparent in the perturbative analysis.« less

Using scattering theory to compute invariant manifolds and numerical results for the laser-driven Hénon-Heiles system.

PubMed

Blazevski, Daniel; Franklin, Jennifer

2012-12-01

Scattering theory is a convenient way to describe systems that are subject to time-dependent perturbations which are localized in time. Using scattering theory, one can compute time-dependent invariant objects for the perturbed system knowing the invariant objects of the unperturbed system. In this paper, we use scattering theory to give numerical computations of invariant manifolds appearing in laser-driven reactions. In this setting, invariant manifolds separate regions of phase space that lead to different outcomes of the reaction and can be used to compute reaction rates.

A numerical study of electromagnetic scattering from ocean like surfaces

NASA Technical Reports Server (NTRS)

Lentz, R. R.

1972-01-01

The integral equations describing electromagnetic scattering from one dimensional conducting surfaces are formulated and numerical results are presented. The results are compared with those obtained using approximate methods such as physical optics, geometrical optics, and perturbation theory. The integral equation solutions show that the surface radius of curvature must be greater than 2.5 wavelengths for either the physical optics or geometric optics to give satisfactory results. It has also been shown that perturbation theory agrees with the exact fields as long as the root mean square surface roughness is less than one-tenth of a wavelength.

Numerical simulation of plasma response to externally applied resonant magnetic perturbation on the J-TEXT tokamak

NASA Astrophysics Data System (ADS)

Bicheng, LI; Zhonghe, JIANG; Jian, LV; Xiang, LI; Bo, RAO; Yonghua, DING

2018-05-01

Nonlinear magnetohydrodynamic (MHD) simulations of an equilibrium on the J-TEXT tokamak with applied resonant magnetic perturbations (RMPs) are performed with NIMROD (non-ideal MHD with rotation, open discussion). Numerical simulation of plasma response to RMPs has been developed to investigate magnetic topology, plasma density and rotation profile. The results indicate that the pure applied RMPs can stimulate 2/1 mode as well as 3/1 mode by the toroidal mode coupling, and finally change density profile by particle transport. At the same time, plasma rotation plays an important role during the entire evolution process.

Integrability and structural stability of solutions to the Ginzburg-Landau equation

NASA Technical Reports Server (NTRS)

Keefe, Laurence R.

1986-01-01

The integrability of the Ginzburg-Landau equation is studied to investigate if the existence of chaotic solutions found numerically could have been predicted a priori. The equation is shown not to possess the Painleveproperty, except for a special case of the coefficients that corresponds to the integrable, nonlinear Schroedinger (NLS) equation. Regarding the Ginzburg-Landau equation as a dissipative perturbation of the NLS, numerical experiments show all but one of a family of two-tori solutions, possessed by the NLS under particular conditions, to disappear under real perturbations to the NLS coefficients of O(10 to the -6th).

Entanglement dynamics in critical random quantum Ising chain with perturbations

DOE Office of Scientific and Technical Information (OSTI.GOV)

Huang, Yichen, E-mail: ychuang@caltech.edu

We simulate the entanglement dynamics in a critical random quantum Ising chain with generic perturbations using the time-evolving block decimation algorithm. Starting from a product state, we observe super-logarithmic growth of entanglement entropy with time. The numerical result is consistent with the analytical prediction of Vosk and Altman using a real-space renormalization group technique. - Highlights: • We study the dynamical quantum phase transition between many-body localized phases. • We simulate the dynamics of a very long random spin chain with matrix product states. • We observe numerically super-logarithmic growth of entanglement entropy with time.

Input-Independent Energy Harvesting in Bistable Lattices from Transition Waves.

PubMed

Hwang, Myungwon; Arrieta, Andres F

2018-02-26

We demonstrate the utilisation of transition waves for realising input-invariant, frequency-independent energy harvesting in 1D lattices of bistable elements. We propose a metamaterial-inspired design with an integrated electromechanical transduction mechanism to the unit cell, rendering the power conversion capability an intrinsic property of the lattice. Moreover, focusing of transmitted energy to desired locations is demonstrated numerically and experimentally by introducing engineered defects in the form of perturbation in mass or inter-element forcing. We achieve further localisation of energy and numerically observe a breather-like mode for the first time in this type of lattice, improving the harvesting performance by an order of magnitude. Our approach considers generic bistable unit cells and thus provides a universal mechanism to harvest energy and realise metamaterials effectively behaving as a capacitor and power delivery system.

Low-thrust trajectory analysis for the geosynchronous mission

NASA Technical Reports Server (NTRS)

Jasper, T. P.

1973-01-01

Methodology employed in development of a computer program designed to analyze optimal low-thrust trajectories is described, and application of the program to a Solar Electric Propulsion Stage (SEPS) geosynchronous mission is discussed. To avoid the zero inclination and eccentricity singularities which plague many small-force perturbation techniques, a special set of state variables (equinoctial) is used. Adjoint equations are derived for the minimum time problem and are also free from the singularities. Solutions to the state and adjoint equations are obtained by both orbit averaging and precision numerical integration; an evaluation of these approaches is made.

Linear and nonlinear pattern selection in Rayleigh-Benard stability problems

NASA Technical Reports Server (NTRS)

Davis, Sanford S.

1993-01-01

A new algorithm is introduced to compute finite-amplitude states using primitive variables for Rayleigh-Benard convection on relatively coarse meshes. The algorithm is based on a finite-difference matrix-splitting approach that separates all physical and dimensional effects into one-dimensional subsets. The nonlinear pattern selection process for steady convection in an air-filled square cavity with insulated side walls is investigated for Rayleigh numbers up to 20,000. The internalization of disturbances that evolve into coherent patterns is investigated and transient solutions from linear perturbation theory are compared with and contrasted to the full numerical simulations.

SCISEAL: A CFD code for analysis of fluid dynamic forces in seals

NASA Technical Reports Server (NTRS)

Athavale, Mahesh; Przekwas, Andrzej

1994-01-01

A viewgraph presentation is made of the objectives, capabilities, and test results of the computer code SCISEAL. Currently, the seal code has: a finite volume, pressure-based integration scheme; colocated variables with strong conservation approach; high-order spatial differencing, up to third-order; up to second-order temporal differencing; a comprehensive set of boundary conditions; a variety of turbulence models and surface roughness treatment; moving grid formulation for arbitrary rotor whirl; rotor dynamic coefficients calculated by the circular whirl and numerical shaker methods; and small perturbation capabilities to handle centered and eccentric seals.

A Green's function formulation for a nonlinear potential flow solution applicable to transonic flow

NASA Technical Reports Server (NTRS)

Baker, A. J.; Fox, C. H., Jr.

1977-01-01

Routine determination of inviscid subsonic flow fields about wing-body-tail configurations employing a Green's function approach for numerical solution of the perturbation velocity potential equation is successfully extended into the high subsonic subcritical flow regime and into the shock-free supersonic flow regime. A modified Green's function formulation, valid throughout a range of Mach numbers including transonic, that takes an explicit accounting of the intrinsic nonlinearity in the parent governing partial differential equations is developed. Some considerations pertinent to flow field predictions in the transonic flow regime are discussed.

Analysis of airfoil leading edge separation bubbles

NASA Technical Reports Server (NTRS)

Carter, J. E.; Vatsa, V. N.

1982-01-01

A local inviscid-viscous interaction technique was developed for the analysis of low speed airfoil leading edge transitional separation bubbles. In this analysis an inverse boundary layer finite difference analysis is solved iteratively with a Cauchy integral representation of the inviscid flow which is assumed to be a linear perturbation to a known global viscous airfoil analysis. Favorable comparisons with data indicate the overall validity of the present localized interaction approach. In addition numerical tests were performed to test the sensitivity of the computed results to the mesh size, limits on the Cauchy integral, and the location of the transition region.

Trajectory optimization and guidance law development for national aerospace plane applications

NASA Technical Reports Server (NTRS)

Calise, A. J.; Flandro, G. A.; Corban, J. E.

1988-01-01

The work completed to date is comprised of the following: a simple vehicle model representative of the aerospace plane concept in the hypersonic flight regime, fuel-optimal climb profiles for the unconstrained and dynamic pressure constrained cases generated using a reduced order dynamic model, an analytic switching condition for transition to rocket powered flight as orbital velocity is approached, simple feedback guidance laws for both the unconstrained and dynamic pressure constrained cases derived via singular perturbation theory and a nonlinear transformation technique, and numerical simulation results for ascent to orbit in the dynamic pressure constrained case.

Long-lived oscillons from asymmetric bubbles: Existence and stability

NASA Astrophysics Data System (ADS)

Adib, Artur B.; Gleiser, Marcelo; Almeida, Carlos A.

2002-10-01

The possibility that extremely long-lived, time-dependent, and localized field configurations (``oscillons'') arise during the collapse of asymmetrical bubbles in (2+1)-dimensional φ4 models is investigated. It is found that oscillons can develop from a large spectrum of elliptically deformed bubbles. Moreover, we provide numerical evidence that such oscillons are (a) circularly symmetric and (b) linearly stable against small arbitrary radial and angular perturbations. The latter is based on a dynamical approach designed to investigate the stability of nonintegrable time-dependent configurations that is capable of probing slowly growing instabilities not seen through the usual ``spectral'' method.

Techniques of orbital decay and long-term ephemeris prediction for satellites in earth orbit

NASA Technical Reports Server (NTRS)

Barry, B. F.; Pimm, R. S.; Rowe, C. K.

1971-01-01

In the special perturbation method, Cowell and variation-of-parameters formulations of the motion equations are implemented and numerically integrated. Variations in the orbital elements due to drag are computed using the 1970 Jacchia atmospheric density model, which includes the effects of semiannual variations, diurnal bulge, solar activity, and geomagnetic activity. In the general perturbation method, two-variable asymptotic series and automated manipulation capabilities are used to obtain analytical solutions to the variation-of-parameters equations. Solutions are obtained considering the effect of oblateness only and the combined effects of oblateness and drag. These solutions are then numerically evaluated by means of a FORTRAN program in which an updating scheme is used to maintain accurate epoch values of the elements. The atmospheric density function is approximated by a Fourier series in true anomaly, and the 1970 Jacchia model is used to periodically update the Fourier coefficients. The accuracy of both methods is demonstrated by comparing computed orbital elements to actual elements over time spans of up to 8 days for the special perturbation method and up to 356 days for the general perturbation method.

Effect of the nonlocal exchange on the performance of the orbital-dependent correlation functionals from second-order perturbation theory.

PubMed

Schweigert, Igor V; Bartlett, Rodney J

2008-09-28

Adding a fraction of the nonlocal exchange operator to the local orbital-dependent exchange potential improves the many-body perturbation expansion based on the Kohn-Sham determinant. The effect of such a hybrid scheme on the performance of the orbital-dependent correlation functional from the second-order perturbation theory (PT2H) is investigated numerically. A small fraction of the nonlocal exchange is often sufficient to ensure the existence of the self-consistent solution for the PT2H potential. In the He and Be atoms, including 37% of the nonlocal exchange leads to the correlation energies and electronic densities that are very close to the exact ones. In molecules, varying the fraction of the nonlocal exchange may result in the PT2H energy closely reproducing the CCSD(T) value; however such a fraction depends on the system and does not always result in an accurate electronic density. We also numerically verify that the "semicanonical" perturbation series includes most of the beneficial effects of the nonlocal exchange without sacrificing the locality of the exchange potential.

Linear perturbations of a Schwarzschild blackhole by thin disc - convergence

NASA Astrophysics Data System (ADS)

Čížek, P.; Semerák, O.

2012-07-01

In order to find the perturbation of a Schwarzschild space-time due to a rotating thin disc, we try to adjust the method used by [4] in the case of perturbation by a one-dimensional ring. This involves solution of stationary axisymmetric Einstein's equations in terms of spherical-harmonic expansions whose convergence however turned out questionable in numerical examples. Here we show, analytically, that the series are almost everywhere convergent, but in some regions the convergence is not absolute.

Scalar perturbation in symmetric Lee-Wick bouncing universe

DOE Office of Scientific and Technical Information (OSTI.GOV)

Cho, Inyong; Kwon, O-Kab, E-mail: iycho@seoultech.ac.kr, E-mail: okab@skku.edu

2011-11-01

We investigate the scalar perturbation in the Lee-Wick bouncing universe driven by an ordinary scalar field plus a ghost field. We consider only a symmetric evolution of the universe and the scalar fields about the bouncing point. The gauge invariant Sasaki-Mukhanov variable is numerically solved in the spatially flat gauge. We find a new form of the initial perturbation growing during the contracting phase. After the bouncing, this growing mode stabilizes to a constant mode which is responsible for the late-time power spectrum.

N-soliton interactions: Effects of linear and nonlinear gain and loss

NASA Astrophysics Data System (ADS)

Carretero-González, R.; Gerdjikov, V. S.; Todorov, M. D.

2017-10-01

We analyze the dynamical behavior of the N-soliton train in the adiabatic approximation of the nonlinear Schrödinger equation perturbed simultaneously by linear and nonlinear gain/loss terms. We derive the corresponding perturbed complex Toda chain in the case of a combination of linear, cubic, and/or quintic terms. We show that the soliton interactions dynamics for this reduced PCTC model compares favorably to full numerical results of the original perturbed nonlinear Schrödinger equation.

Do cosmological data rule out f (R ) with w ≠-1 ?

NASA Astrophysics Data System (ADS)

Battye, Richard A.; Bolliet, Boris; Pace, Francesco

2018-05-01

We review the equation of state (EoS) approach to dark sector perturbations and apply it to f (R ) gravity models of dark energy. We show that the EoS approach is numerically stable and use it to set observational constraints on designer models. Within the EoS approach we build an analytical understanding of the dynamics of cosmological perturbations for the designer class of f (R ) gravity models, characterized by the parameter B0 and the background equation of state of dark energy w . When we use the Planck cosmic microwave background temperature anisotropy, polarization, and lensing data as well as the baryonic acoustic oscillation data from SDSS and WiggleZ, we find B0<0.006 (95% C.L.) for the designer models with w =-1 . Furthermore, we find B0<0.0045 and |w +1 |<0.002 (95% C.L.) for the designer models with w ≠-1 . Previous analyses found similar results for designer and Hu-Sawicki f (R ) gravity models using the effective field theory approach [Raveri et al., Phys. Rev. D 90, 043513 (2014), 10.1103/PhysRevD.90.043513; Hu et al., Mon. Not. R. Astron. Soc. 459, 3880 (2016), 10.1093/mnras/stw775]; therefore this hints for the fact that generic f (R ) models with w ≠-1 can be tightly constrained by current cosmological data, complementary to solar system tests [Brax et al., Phys. Rev. D 78, 104021 (2008), 10.1103/PhysRevD.78.104021; Faulkner et al., Phys. Rev. D 76, 063505 (2007), 10.1103/PhysRevD.76.063505]. When compared to a w CDM fluid with the same sound speed, we find that the equation of state for f (R ) models is better constrained to be close to -1 by about an order of magnitude, due to the strong dependence of the perturbations on w .

A laboratory experiment simulating the dynamics of topographic relief: methodology and results

NASA Astrophysics Data System (ADS)

Crave, A.; Lague, D.; Davy, P.; Bonnet, S.; Laguionie, P.

2002-12-01

Theoretical analysis and numerical models of landscape evolution have advanced several scenarios for the long-term evolution of terrestrial topography. These scenarios require quantitative evaluation. Analyses of topography, sediment fluxes, and the physical mechanisms of erosion and sediment transport can provide some constraints on the range of plausible models. But in natural systems the boundary conditions (tectonic uplift, climate, base level) are often not well constrained and the spatial heterogeneity of substrate, climate, vegetation, and prevalent processes commonly confounds attempts at extrapolation of observations to longer timescales. In the laboratory, boundary conditions are known and heterogeneity and complexity can be controlled. An experimental approach can thus provide valuable constraints on the dynamics of geomorphic systems, provided that (1) the elementary processes are well calibrated and (2) the topography and sediment fluxes are sufficiently well documented. We have built an experimental setup of decimeter scale that is designed to develop a complete drainage network by the growth and propagation of erosion instabilities in response to tectonic and climatic perturbations. Uplift and precipitation rates can be changed over an order of magnitude. Telemetric lasers and 3D stereo-photography allow the precise quantification of the topographic evolution of the experimental surface. In order to calibrate the principal processes of erosion and transport we have used three approaches: (1) theoretical derivation of erosion laws deduced from the geometrical properties of experimental surfaces at steady-state under different rates of tectonic uplift; (2) comparison of the experimental transient dynamics with a numerical simulation model to test the validity of the predicted erosion laws; and (3) detailed analysis of particle detachment and transport in a millimeter sheet flow on a two-meter long flume under precisely controlled water discharge, slope and flow width. The analogy with real geomorphic systems is limited by the imperfect downscaling in both time and space of the experiments. However, these simple experiments have allowed us to probe (1) the importance of a threshold for particle mobilization to the relationship between steady-state elevation and uplift rate, (2) the role of initial drainage network organization in the transient dynamics of tectonically perturbed systems and (3) the sediment flux dynamics of climatically perturbed systems.

Computational singular perturbation analysis of stochastic chemical systems with stiffness

NASA Astrophysics Data System (ADS)

Wang, Lijin; Han, Xiaoying; Cao, Yanzhao; Najm, Habib N.

2017-04-01

Computational singular perturbation (CSP) is a useful method for analysis, reduction, and time integration of stiff ordinary differential equation systems. It has found dominant utility, in particular, in chemical reaction systems with a large range of time scales at continuum and deterministic level. On the other hand, CSP is not directly applicable to chemical reaction systems at micro or meso-scale, where stochasticity plays an non-negligible role and thus has to be taken into account. In this work we develop a novel stochastic computational singular perturbation (SCSP) analysis and time integration framework, and associated algorithm, that can be used to not only construct accurately and efficiently the numerical solutions to stiff stochastic chemical reaction systems, but also analyze the dynamics of the reduced stochastic reaction systems. The algorithm is illustrated by an application to a benchmark stochastic differential equation model, and numerical experiments are carried out to demonstrate the effectiveness of the construction.

Investigation of the spatial structure and developmental dynamics of near-Earth plasma perturbations under the action of powerful HF radio waves

DOE Office of Scientific and Technical Information (OSTI.GOV)

Belov, A. S., E-mail: alexis-belov@yandex.ru

2015-10-15

Results of numerical simulations of the near-Earth plasma perturbations induced by powerful HF radio waves from the SURA heating facility are presented. The simulations were performed using a modified version of the SAMI2 ionospheric model for the input parameters corresponding to the series of in-situ SURA–DEMETER experiments. The spatial structure and developmental dynamics of large-scale plasma temperature and density perturbations have been investigated. The characteristic formation and relaxation times of the induced large-scale plasma perturbations at the altitudes of the Earth’s outer ionosphere have been determined.

Least Squares Moving-Window Spectral Analysis.

PubMed

Lee, Young Jong

2017-08-01

Least squares regression is proposed as a moving-windows method for analysis of a series of spectra acquired as a function of external perturbation. The least squares moving-window (LSMW) method can be considered an extended form of the Savitzky-Golay differentiation for nonuniform perturbation spacing. LSMW is characterized in terms of moving-window size, perturbation spacing type, and intensity noise. Simulation results from LSMW are compared with results from other numerical differentiation methods, such as single-interval differentiation, autocorrelation moving-window, and perturbation correlation moving-window methods. It is demonstrated that this simple LSMW method can be useful for quantitative analysis of nonuniformly spaced spectral data with high frequency noise.

Time-variant random interval natural frequency analysis of structures

NASA Astrophysics Data System (ADS)

Wu, Binhua; Wu, Di; Gao, Wei; Song, Chongmin

2018-02-01

This paper presents a new robust method namely, unified interval Chebyshev-based random perturbation method, to tackle hybrid random interval structural natural frequency problem. In the proposed approach, random perturbation method is implemented to furnish the statistical features (i.e., mean and standard deviation) and Chebyshev surrogate model strategy is incorporated to formulate the statistical information of natural frequency with regards to the interval inputs. The comprehensive analysis framework combines the superiority of both methods in a way that computational cost is dramatically reduced. This presented method is thus capable of investigating the day-to-day based time-variant natural frequency of structures accurately and efficiently under concrete intrinsic creep effect with probabilistic and interval uncertain variables. The extreme bounds of the mean and standard deviation of natural frequency are captured through the embedded optimization strategy within the analysis procedure. Three particularly motivated numerical examples with progressive relationship in perspective of both structure type and uncertainty variables are demonstrated to justify the computational applicability, accuracy and efficiency of the proposed method.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Albert, Christopher G.; Heyn, Martin F.; Kapper, Gernot

Toroidal torque generated by neoclassical viscosity caused by external non-resonant, non-axisymmetric perturbations has a significant influence on toroidal plasma rotation in tokamaks. In this article, a derivation for the expressions of toroidal torque and radial transport in resonant regimes is provided within quasilinear theory in canonical action-angle variables. The proposed approach treats all low-collisional quasilinear resonant neoclassical toroidal viscosity regimes including superbanana-plateau and drift-orbit resonances in a unified way and allows for magnetic drift in all regimes. It is valid for perturbations on toroidally symmetric flux surfaces of the unperturbed equilibrium without specific assumptions on geometry or aspect ratio. Themore » resulting expressions are shown to match the existing analytical results in the large aspect ratio limit. Numerical results from the newly developed code NEO-RT are compared to calculations by the quasilinear version of the code NEO-2 at low collisionalities. The importance of the magnetic shear term in the magnetic drift frequency and a significant effect of the magnetic drift on drift-orbit resonances are demonstrated.« less

Analytical Solution of Coupled Perturbation of Tesseral Harmonic Terms of Mars's Non-Spherical Gravitational Potential

NASA Astrophysics Data System (ADS)

Zhou, Chui-hong; Yu, Sheng-xian; Liu, Lin

2012-10-01

The non-spherical gravitational potential of the planet Mars is sig- nificantly different from that of the Earth. The magnitudes of Mars' tesseral harmonic coefficients are basically ten times larger than the corresponding val- ues of the Earth. Especially, the magnitude of its second degree and order tesseral harmonic coefficient J2,2 is nearly 40 times that of the Earth, and approaches to the one tenth of its second zonal harmonic coefficient J2. For a low-orbit Mars probe, if the required accuracy of orbit prediction of 1-day arc length is within 500 m (equivalent to the order of magnitude of 10-4 standard unit), then the coupled terms of J2 with the tesseral harmonics, and even those of the tesseral harmonics themselves, which are negligible for the Earth satellites, should be considered when the analytical perturbation solution of its orbit is built. In this paper, the analytical solutions of the coupled terms are presented. The anal- ysis and numerical verification indicate that the effect of the above-mentioned coupled perturbation on the orbit may exceed 10-4 in the along-track direc- tion. The conclusion is that the solutions of Earth satellites cannot be simply used without any modification when dealing with the analytical perturbation solutions of Mars-orbiting satellites, and that the effect of the coupled terms of Mars's non-spherical gravitational potential discussed in this paper should be taken into consideration.

Dependence of locked mode behavior on frequency and polarity of a rotating external magnetic perturbation

NASA Astrophysics Data System (ADS)

Inoue, S.; Shiraishi, J.; Takechi, M.; Matsunaga, G.; Isayama, A.; Hayashi, N.; Ide, S.

2018-02-01

Active control and stabilization of locked modes (LM) via rotating external magnetic perturbations are numerically investigated under a realistic low resistivity condition. To explore plasma responses to rotating and/or static external magnetic perturbations, we have developed a resistive magnetohydrodynamic code ‘AEOLUS-IT’. By using AEOLUS-IT, dependencies of mode behavior on frequency and polarity of the rotating magnetic perturbation are successfully clarified. Here, the rotational direction of the rotating magnetic perturbation to the equilibrium plasma rotation in the laboratory frame is referred to as ‘polarity’. The rotating magnetic perturbation acts on the background rotating plasma in the presence of a static field. Under such circumstances, there exist bifurcated states of the background rotating plasma, which should be taken into account when studying the dependence of the mode behavior on the rotating magnetic perturbation. It is found that there exist an optimum frequency and polarity of the rotating magnetic perturbation to control the LM, and that the LM is effectively stabilized by a co-polarity magnetic perturbation in comparison with a counter-polarity one.

Constrained orbital intercept-evasion

NASA Astrophysics Data System (ADS)

Zatezalo, Aleksandar; Stipanovic, Dusan M.; Mehra, Raman K.; Pham, Khanh

2014-06-01

An effective characterization of intercept-evasion confrontations in various space environments and a derivation of corresponding solutions considering a variety of real-world constraints are daunting theoretical and practical challenges. Current and future space-based platforms have to simultaneously operate as components of satellite formations and/or systems and at the same time, have a capability to evade potential collisions with other maneuver constrained space objects. In this article, we formulate and numerically approximate solutions of a Low Earth Orbit (LEO) intercept-maneuver problem in terms of game-theoretic capture-evasion guaranteed strategies. The space intercept-evasion approach is based on Liapunov methodology that has been successfully implemented in a number of air and ground based multi-player multi-goal game/control applications. The corresponding numerical algorithms are derived using computationally efficient and orbital propagator independent methods that are previously developed for Space Situational Awareness (SSA). This game theoretical but at the same time robust and practical approach is demonstrated on a realistic LEO scenario using existing Two Line Element (TLE) sets and Simplified General Perturbation-4 (SGP-4) propagator.

Aeroelastic Stability of Rotor Blades Using Finite Element Analysis

NASA Technical Reports Server (NTRS)

Chopra, I.; Sivaneri, N.

1982-01-01

The flutter stability of flap bending, lead-lag bending, and torsion of helicopter rotor blades in hover is investigated using a finite element formulation based on Hamilton's principle. The blade is divided into a number of finite elements. Quasi-steady strip theory is used to evaluate the aerodynamic loads. The nonlinear equations of motion are solved for steady-state blade deflections through an iterative procedure. The equations of motion are linearized assuming blade motion to be a small perturbation about the steady deflected shape. The normal mode method based on the coupled rotating natural modes is used to reduce the number of equations in the flutter analysis. First the formulation is applied to single-load-path blades (articulated and hingeless blades). Numerical results show very good agreement with existing results obtained using the modal approach. The second part of the application concerns multiple-load-path blades, i.e. bearingless blades. Numerical results are presented for several analytical models of the bearingless blade. Results are also obtained using an equivalent beam approach wherein a bearingless blade is modelled as a single beam with equivalent properties. Results show the equivalent beam model.

Approximate analytical solutions to the double-stance dynamics of the lossy spring-loaded inverted pendulum.

PubMed

Shahbazi, Mohammad; Saranlı, Uluç; Babuška, Robert; Lopes, Gabriel A D

2016-12-05

This paper introduces approximate time-domain solutions to the otherwise non-integrable double-stance dynamics of the 'bipedal' spring-loaded inverted pendulum (B-SLIP) in the presence of non-negligible damping. We first introduce an auxiliary system whose behavior under certain conditions is approximately equivalent to the B-SLIP in double-stance. Then, we derive approximate solutions to the dynamics of the new system following two different methods: (i) updated-momentum approach that can deal with both the lossy and lossless B-SLIP models, and (ii) perturbation-based approach following which we only derive a solution to the lossless case. The prediction performance of each method is characterized via a comprehensive numerical analysis. The derived representations are computationally very efficient compared to numerical integrations, and, hence, are suitable for online planning, increasing the autonomy of walking robots. Two application examples of walking gait control are presented. The proposed solutions can serve as instrumental tools in various fields such as control in legged robotics and human motion understanding in biomechanics.

The Harmonic Oscillator with a Gaussian Perturbation: Evaluation of the Integrals and Example Applications

ERIC Educational Resources Information Center

Earl, Boyd L.

2008-01-01

A general result for the integrals of the Gaussian function over the harmonic oscillator wavefunctions is derived using generating functions. Using this result, an example problem of a harmonic oscillator with various Gaussian perturbations is explored in order to compare the results of precise numerical solution, the variational method, and…

Dynamical evolution of space debris on high-elliptical orbits near high-order resonance zones

NASA Astrophysics Data System (ADS)

Kuznetsov, Eduard; Zakharova, Polina

Orbital evolution of objects on Molniya-type orbits is considered near high-order resonance zones. Initial conditions correspond to high-elliptical orbits with the critical inclination 63.4 degrees. High-order resonances are analyzed. Resonance orders are more than 5 and less than 50. Frequencies of perturbations caused by the effect of sectorial and tesseral harmonics of the Earth's gravitational potential are linear combinations of the mean motion of a satellite, angular velocities of motion of the pericenter and node of its orbit, and the angular velocity of the Earth. Frequencies of perturbations were calculated by taking into account secular perturbations from the Earth oblateness, the Moon, the Sun, and a solar radiation pressure. Resonance splitting effect leads to three sub-resonances. The study of dynamical evolution on long time intervals was performed on the basis of the results of numerical simulation. We used "A Numerical Model of the Motion of Artificial Earth's Satellites", developed by the Research Institute of Applied Mathematics and Mechanics of the Tomsk State University. The model of disturbing forces taken into account the main perturbing factors: the gravitational field of the Earth, the attraction of the Moon and the Sun, the tides in the Earth’s body, the solar radiation pressure, taking into account the shadow of the Earth, the Poynting-Robertson effect, and the atmospheric drag. Area-to-mass ratio varied from small values corresponding to satellites to big ones corresponding to space debris. The locations and sizes of resonance zones were refined from numerical simulation. The Poynting-Robertson effect results in a secular decrease in the semi-major axis of a spherically symmetrical satellite. In resonance regions the effect weakens slightly. Reliable estimates of secular perturbations of the semi-major axis were obtained from the numerical simulation. Under the Poynting-Robertson effect objects pass through the regions of high-order resonances. The Poynting-Robertson effect and secular perturbations of the semi-major axis lead to formation weak stochastic trajectories. The integral autocorrelation function was used to analysis stochastic properties trajectories. This work was supported by the Ministry of Education and Science of the Russian Federation and the Russian Foundation for Basic Researches (grant 13-02-00026a).

Analysis of pulsating spray flames propagating in lean two-phase mixtures with unity Lewis number

DOE Office of Scientific and Technical Information (OSTI.GOV)

Nicoli, C.; Haldenwang, P.; Suard, S.

2005-11-01

Pulsating (or oscillatory) spray flames have recently been observed in experiments on two-phase combustion. Numerical studies have pointed out that such front oscillations can be obtained even with very simple models of homogeneous two-phase mixtures, including elementary vaporization schemes. The paper presents an analytical approach within the simple framework of the thermal-diffusive model, which is complemented by a vaporization rate independent of gas temperature, as soon as the latter reaches a certain thermal threshold ({theta}{sub v} in reduced form). The study involves the Damkoehler number (Da), the ratio of chemical reaction rate to vaporization rate, and the Zeldovich number (Ze)more » as essential parameters. We use the standard asymptotic method based on matched expansions in terms of 1/Ze. Linear analysis of two-phase flame stability is performed by studying, in the absence of differential diffusive effects (unity Lewis number), the linear growth rate of 2-D perturbations added to steady plane solutions and characterized by wavenumber k in the direction transverse to spreading. A domain of existence is found for the pulsating regime. It corresponds to mixture characteristics often met in air-fuel two-phase systems: low boiling temperature ({theta}{sub v} << 1), reaction rate not higher than vaporization rate (Da < 1, i.e., small droplets), and activation temperature assumed to be high compared with flame temperature (Ze {>=} 10). Satisfactory comparison with numerical simulations confirms the validity of the analytical approach; in particular, positive growth rates have been found for planar perturbations (k = 0) and for wrinkled fronts (k {ne} 0). Finally, comparison between predicted frequencies and experimental measurements is discussed.« less

Towards accurate cosmological predictions for rapidly oscillating scalar fields as dark matter

DOE Office of Scientific and Technical Information (OSTI.GOV)

Ureña-López, L. Arturo; Gonzalez-Morales, Alma X., E-mail: lurena@ugto.mx, E-mail: alma.gonzalez@fisica.ugto.mx

2016-07-01

As we are entering the era of precision cosmology, it is necessary to count on accurate cosmological predictions from any proposed model of dark matter. In this paper we present a novel approach to the cosmological evolution of scalar fields that eases their analytic and numerical analysis at the background and at the linear order of perturbations. The new method makes use of appropriate angular variables that simplify the writing of the equations of motion, and which also show that the usual field variables play a secondary role in the cosmological dynamics. We apply the method to a scalar fieldmore » endowed with a quadratic potential and revisit its properties as dark matter. Some of the results known in the literature are recovered, and a better understanding of the physical properties of the model is provided. It is confirmed that there exists a Jeans wavenumber k {sub J} , directly related to the suppression of linear perturbations at wavenumbers k > k {sub J} , and which is verified to be k {sub J} = a √ mH . We also discuss some semi-analytical results that are well satisfied by the full numerical solutions obtained from an amended version of the CMB code CLASS. Finally we draw some of the implications that this new treatment of the equations of motion may have in the prediction of cosmological observables from scalar field dark matter models.« less

First-order analytic propagation of satellites in the exponential atmosphere of an oblate planet

NASA Astrophysics Data System (ADS)

Martinusi, Vladimir; Dell'Elce, Lamberto; Kerschen, Gaëtan

2017-04-01

The paper offers the fully analytic solution to the motion of a satellite orbiting under the influence of the two major perturbations, due to the oblateness and the atmospheric drag. The solution is presented in a time-explicit form, and takes into account an exponential distribution of the atmospheric density, an assumption that is reasonably close to reality. The approach involves two essential steps. The first one concerns a new approximate mathematical model that admits a closed-form solution with respect to a set of new variables. The second step is the determination of an infinitesimal contact transformation that allows to navigate between the new and the original variables. This contact transformation is obtained in exact form, and afterwards a Taylor series approximation is proposed in order to make all the computations explicit. The aforementioned transformation accommodates both perturbations, improving the accuracy of the orbit predictions by one order of magnitude with respect to the case when the atmospheric drag is absent from the transformation. Numerical simulations are performed for a low Earth orbit starting at an altitude of 350 km, and they show that the incorporation of drag terms into the contact transformation generates an error reduction by a factor of 7 in the position vector. The proposed method aims at improving the accuracy of analytic orbit propagation and transforming it into a viable alternative to the computationally intensive numerical methods.

SMD-based numerical stochastic perturbation theory

NASA Astrophysics Data System (ADS)

Dalla Brida, Mattia; Lüscher, Martin

2017-05-01

The viability of a variant of numerical stochastic perturbation theory, where the Langevin equation is replaced by the SMD algorithm, is examined. In particular, the convergence of the process to a unique stationary state is rigorously established and the use of higher-order symplectic integration schemes is shown to be highly profitable in this context. For illustration, the gradient-flow coupling in finite volume with Schrödinger functional boundary conditions is computed to two-loop (i.e. NNL) order in the SU(3) gauge theory. The scaling behaviour of the algorithm turns out to be rather favourable in this case, which allows the computations to be driven close to the continuum limit.

Propagation of large amplitude Alfven waves in the solar wind neutral sheet

NASA Technical Reports Server (NTRS)

Malara, F.; Primavera, L.; Veltri, P.

1995-01-01

Analysis of solar wind fluctuation data show that the correlation between velocity and magnetic field fluctuations decreases when going farther away from the Sun. This decorrelation can be attributed either to the time evolution of the fluctuations, carried away by the solar wind, or to the interaction between the solar wind neutral sheet and Alfven waves. To check this second hypothesis we have numerically studied the propagation of Alfven waves in the solar wind neutral sheet. The initial conditions have been set up in order to guarantee B(exp 2) = const, so that the following numerical evolution is only due to the inhomogeneity in the background magnetic field. The analysis of the results shows that compressive structures are formed, mainly in the neutral sheet where they have been identified as pressure balanced structures, i.e., tangential discontinuities. Fast perturbations, which are also produced, have a tendency to leave the simulation domain, propagating also perpendicularly to the mean magnetic field. For this reason the level of fast perturbations is always smaller with respect to the previously cited plasma balanced structures, which are slow mode perturbations. A comparison between the numerical results and some particular observational issues is also presented.

Curvature perturbation spectra from waterfall transition, black hole constraints and non-Gaussianity

NASA Astrophysics Data System (ADS)

Bugaev, Edgar; Klimai, Peter

2011-11-01

We carried out numerical calculations of a contribution of the waterfall field to the primordial curvature perturbation (on uniform density hypersurfaces) ζ, which is produced during waterfall transition in hybrid inflation scenario. The calculation is performed for a broad interval of values of the model parameters. We show that there is a strong growth of amplitudes of the curvature perturbation spectrum in the limit when the bare mass-squared of the waterfall field becomes comparable with the square of Hubble parameter. We show that in this limit the primordial black hole constraints on the curvature perturbations must be taken into account. It is shown that, in the same limit, peak values of the curvature perturbation spectra are far beyond horizon, and the spectra are strongly non-Gaussian.

Artificial Bee Colony Optimization of Capping Potentials for Hybrid Quantum Mechanical/Molecular Mechanical Calculations.

PubMed

Schiffmann, Christoph; Sebastiani, Daniel

2011-05-10

We present an algorithmic extension of a numerical optimization scheme for analytic capping potentials for use in mixed quantum-classical (quantum mechanical/molecular mechanical, QM/MM) ab initio calculations. Our goal is to minimize bond-cleavage-induced perturbations in the electronic structure, measured by means of a suitable penalty functional. The optimization algorithm-a variant of the artificial bee colony (ABC) algorithm, which relies on swarm intelligence-couples deterministic (downhill gradient) and stochastic elements to avoid local minimum trapping. The ABC algorithm outperforms the conventional downhill gradient approach, if the penalty hypersurface exhibits wiggles that prevent a straight minimization pathway. We characterize the optimized capping potentials by computing NMR chemical shifts. This approach will increase the accuracy of QM/MM calculations of complex biomolecules.

Polaron mobility obtained by a variational approach for lattice Fröhlich models

NASA Astrophysics Data System (ADS)

Kornjača, Milan; Vukmirović, Nenad

2018-04-01

Charge carrier mobility for a class of lattice models with long-range electron-phonon interaction was investigated. The approach for mobility calculation is based on a suitably chosen unitary transformation of the model Hamiltonian which transforms it into the form where the remaining interaction part can be treated as a perturbation. Relevant spectral functions were then obtained using Matsubara Green's functions technique and charge carrier mobility was evaluated using Kubo's linear response formula. Numerical results were presented for a wide range of electron-phonon interaction strengths and temperatures in the case of one-dimensional version of the model. The results indicate that the mobility decreases with increasing temperature for all electron-phonon interaction strengths in the investigated range, while longer interaction range leads to more mobile carriers.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Scherrer, Arne; UMR 8640 ENS-CNRS-UPMC, Département de Chimie, 24 rue Lhomond, École Normale Supérieure, 75005 Paris; UPMC Université Paris 06, 4, Place Jussieu, 75005 Paris

The nuclear velocity perturbation theory (NVPT) for vibrational circular dichroism (VCD) is derived from the exact factorization of the electron-nuclear wave function. This new formalism offers an exact starting point to include correction terms to the Born-Oppenheimer (BO) form of the molecular wave function, similar to the complete-adiabatic approximation. The corrections depend on a small parameter that, in a classical treatment of the nuclei, is identified as the nuclear velocity. Apart from proposing a rigorous basis for the NVPT, we show that the rotational strengths, related to the intensity of the VCD signal, contain a new contribution beyond-BO that canmore » be evaluated with the NVPT and that only arises when the exact factorization approach is employed. Numerical results are presented for chiral and non-chiral systems to test the validity of the approach.« less

PREditOR: a synthetic biology approach to removing heterochromatin from cells.

PubMed

Molina, Oscar; Carmena, Mar; Maudlin, Isabella E; Earnshaw, William C

2016-12-01

It is widely accepted that heterochromatin is necessary to maintain genomic stability. However, direct experimental evidence supporting this is slim. Previous studies using either enzyme inhibitors, gene knockout or knockdown studies all are subject to the caveat that drugs may have off-target effects and enzymes that modify chromatin proteins to support heterochromatin formation may also have numerous other cellular targets as well. Here, we describe PREditOR (protein reading and editing of residues), a synthetic biology approach that allows us to directly remove heterochromatin from cells without either drugs or global interference with gene function. We find that removal of heterochromatin perturbs mitotic progression and causes a dramatic increase in chromosome segregation defects, possibly as a result of interfering with the normal centromeric localization of the chromosomal passenger complex.

A comparative numerical analysis of linear and nonlinear aerodynamic sound generation by vortex disturbances in homentropic constant shear flows

DOE Office of Scientific and Technical Information (OSTI.GOV)

Hau, Jan-Niklas, E-mail: hau@fdy.tu-darmstadt.de; Oberlack, Martin; GSC CE, Technische Universität Darmstadt, Dolivostraße 15, 64293 Darmstadt

2015-12-15

Aerodynamic sound generation in shear flows is investigated in the light of the breakthrough in hydrodynamics stability theory in the 1990s, where generic phenomena of non-normal shear flow systems were understood. By applying the thereby emerged short-time/non-modal approach, the sole linear mechanism of wave generation by vortices in shear flows was captured [G. D. Chagelishvili, A. Tevzadze, G. Bodo, and S. S. Moiseev, “Linear mechanism of wave emergence from vortices in smooth shear flows,” Phys. Rev. Lett. 79, 3178-3181 (1997); B. F. Farrell and P. J. Ioannou, “Transient and asymptotic growth of two-dimensional perturbations in viscous compressible shear flow,” Phys.more » Fluids 12, 3021-3028 (2000); N. A. Bakas, “Mechanism underlying transient growth of planar perturbations in unbounded compressible shear flow,” J. Fluid Mech. 639, 479-507 (2009); and G. Favraud and V. Pagneux, “Superadiabatic evolution of acoustic and vorticity perturbations in Couette flow,” Phys. Rev. E 89, 033012 (2014)]. Its source is the non-normality induced linear mode-coupling, which becomes efficient at moderate Mach numbers that is defined for each perturbation harmonic as the ratio of the shear rate to its characteristic frequency. Based on the results by the non-modal approach, we investigate a two-dimensional homentropic constant shear flow and focus on the dynamical characteristics in the wavenumber plane. This allows to separate from each other the participants of the dynamical processes — vortex and wave modes — and to estimate the efficacy of the process of linear wave-generation. This process is analyzed and visualized on the example of a packet of vortex modes, localized in both, spectral and physical, planes. Further, by employing direct numerical simulations, the wave generation by chaotically distributed vortex modes is analyzed and the involved linear and nonlinear processes are identified. The generated acoustic field is anisotropic in the wavenumber plane, which results in highly directional linear sound radiation, whereas the nonlinearly generated waves are almost omni-directional. As part of this analysis, we compare the effectiveness of the linear and nonlinear mechanisms of wave generation within the range of validity of the rapid distortion theory and show the dominance of the linear aerodynamic sound generation. Finally, topological differences between the linear source term of the acoustic analogy equation and of the anisotropic non-normality induced linear mechanism of wave generation are found.« less

Rupture dynamics along bimaterial interfaces: a parametric study of the coupling between interfacial sliding and normal traction perturbation

NASA Astrophysics Data System (ADS)

Scala, Antonio; Festa, Gaetano; Vilotte, Jean-Pierre

2017-04-01

Earthquake ruptures often develop along faults separating materials with dissimilar elastic properties. Due to the broken symmetry, the propagation of the rupture along the bimaterial interface is driven by the coupling between interfacial sliding and normal traction perturbations. We numerically investigate in-plane rupture growth along a planar interface, under slip weakening friction, separating two dissimilar isotropic linearly elastic half-spaces. We perform a parametric study of the classical Prakash-Clifton regularisation for different material contrasts. In particular mesh-dependence and regularisation-dependence of the numerical solutions are analysed in this parameter space. When regularisation involves a slip-rate dependent relaxation time, a characteristic sliding distance is identified below which numerical solutions no longer depend on the regularisation parameter, i.e. they are consistent solutions of the same physical problem. Such regularisation provides an adaptive high-frequency filter of the slip-induced normal traction perturbations, following the dynamic shrinking of the dissipation zone during the acceleration phase. In contrast, regularisation involving a constant relaxation time leads to numerical solutions that always depend on the regularisation parameter since it fails adapting to the shrinking of the process zone. Dynamic regularisation is further investigated using a non-local regularisation based on a relaxation time that depends on the dynamic length of the dissipation zone. Such reformulation is shown to provide similar results as the dynamic time scale regularisation proposed by Prakash-Clifton when slip rate is replaced by the maximum slip rate along the sliding interface. This leads to the identification of a dissipative length scale associated with the coupling between interfacial sliding and normal traction perturbations, together with a scaling law between the maximum slip rate and the dynamic size of the process zone during the rupture propagation. Dynamic time scale regularisation is show to provide mesh-independent and physically well-posed numerical solutions during the acceleration phase toward an asymptotic speed. When generalised Rayleigh wave does not exist, numerical solutions are shown to tend toward an asymptotic velocity higher than the slowest shear wave speed. When generalised Rayleigh wave speed exists, as numerical solutions tend toward this velocity, increasing spurious oscillations develop and solutions become unstable. In this regime regularisation dependent and unstable finite-size pulses may be generated. This instability is associated with the singular behaviour of the slip-induced normal traction perturbations, and of the slip rate at the rupture front, in relation with complete shrinking of the dissipation zone. This phase requires to be modelled either by more complex interface constitutive laws involving velocity-strengthening effects that may stabilize short wavelength interfacial propagating modes or by considering non-ideal interfaces that introduce a new length scale in the problem that may promote selection and stabilization of the slip pulses.

Bayesian Analysis of the Power Spectrum of the Cosmic Microwave Background

NASA Technical Reports Server (NTRS)

Jewell, Jeffrey B.; Eriksen, H. K.; O'Dwyer, I. J.; Wandelt, B. D.

2005-01-01

There is a wealth of cosmological information encoded in the spatial power spectrum of temperature anisotropies of the cosmic microwave background. The sky, when viewed in the microwave, is very uniform, with a nearly perfect blackbody spectrum at 2.7 degrees. Very small amplitude brightness fluctuations (to one part in a million!!) trace small density perturbations in the early universe (roughly 300,000 years after the Big Bang), which later grow through gravitational instability to the large-scale structure seen in redshift surveys... In this talk, I will discuss a Bayesian formulation of this problem; discuss a Gibbs sampling approach to numerically sampling from the Bayesian posterior, and the application of this approach to the first-year data from the Wilkinson Microwave Anisotropy Probe. I will also comment on recent algorithmic developments for this approach to be tractable for the even more massive data set to be returned from the Planck satellite.

Spectral methods in general relativity and large Randall-Sundrum II black holes

NASA Astrophysics Data System (ADS)

Abdolrahimi, Shohreh; Cattoën, Céline; Page, Don N.; \\\\; Yaghoobpour-Tari, Shima

2013-06-01

Using a novel numerical spectral method, we have found solutions for large static Randall-Sundrum II (RSII) black holes by perturbing a numerical AdS5-CFT4 solution to the Einstein equation with a negative cosmological constant Λ that is asymptotically conformal to the Schwarzschild metric. We used a numerical spectral method independent of the Ricci-DeTurck-flow method used by Figueras, Lucietti, and Wiseman for a similar numerical solution. We have compared our black-hole solution to the one Figueras and Wiseman have derived by perturbing their numerical AdS5-CFT4 solution, showing that our solution agrees closely with theirs. We have obtained a closed-form approximation to the metric of the black hole on the brane. We have also deduced the new results that to first order in 1/(-ΛM2), the Hawking temperature and entropy of an RSII static black hole have the same values as the Schwarzschild metric with the same mass, but the horizon area is increased by about 4.7/(-Λ).

Extended multi-configuration quasi-degenerate perturbation theory: the new approach to multi-state multi-reference perturbation theory.

PubMed

Granovsky, Alexander A

2011-06-07

The distinctive desirable features, both mathematically and physically meaningful, for all partially contracted multi-state multi-reference perturbation theories (MS-MR-PT) are explicitly formulated. The original approach to MS-MR-PT theory, called extended multi-configuration quasi-degenerate perturbation theory (XMCQDPT), having most, if not all, of the desirable properties is introduced. The new method is applied at the second order of perturbation theory (XMCQDPT2) to the 1(1)A(')-2(1)A(') conical intersection in allene molecule, the avoided crossing in LiF molecule, and the 1(1)A(1) to 2(1)A(1) electronic transition in cis-1,3-butadiene. The new theory has several advantages compared to those of well-established approaches, such as second order multi-configuration quasi-degenerate perturbation theory and multi-state-second order complete active space perturbation theory. The analysis of the prevalent approaches to the MS-MR-PT theory performed within the framework of the XMCQDPT theory unveils the origin of their common inherent problems. We describe the efficient implementation strategy that makes XMCQDPT2 an especially useful general-purpose tool in the high-level modeling of small to large molecular systems. © 2011 American Institute of Physics

Asymptotic behaviors of a cell-to-cell HIV-1 infection model perturbed by white noise

NASA Astrophysics Data System (ADS)

Liu, Qun

2017-02-01

In this paper, we analyze a mathematical model of cell-to-cell HIV-1 infection to CD4+ T cells perturbed by stochastic perturbations. First of all, we investigate that there exists a unique global positive solution of the system for any positive initial value. Then by using Lyapunov analysis methods, we study the asymptotic property of this solution. Moreover, we discuss whether there is a stationary distribution for this system and if it owns the ergodic property. Numerical simulations are presented to illustrate the theoretical results.

A SCILAB Program for Computing General-Relativistic Models of Rotating Neutron Stars by Implementing Hartle's Perturbation Method

NASA Astrophysics Data System (ADS)

Papasotiriou, P. J.; Geroyannis, V. S.

We implement Hartle's perturbation method to the computation of relativistic rigidly rotating neutron star models. The program has been written in SCILAB (© INRIA ENPC), a matrix-oriented high-level programming language. The numerical method is described in very detail and is applied to many models in slow or fast rotation. We show that, although the method is perturbative, it gives accurate results for all practical purposes and it should prove an efficient tool for computing rapidly rotating pulsars.

Perturbation-iteration theory for analyzing microwave striplines

NASA Technical Reports Server (NTRS)

Kretch, B. E.

1985-01-01

A perturbation-iteration technique is presented for determining the propagation constant and characteristic impedance of an unshielded microstrip transmission line. The method converges to the correct solution with a few iterations at each frequency and is equivalent to a full wave analysis. The perturbation-iteration method gives a direct solution for the propagation constant without having to find the roots of a transcendental dispersion equation. The theory is presented in detail along with numerical results for the effective dielectric constant and characteristic impedance for a wide range of substrate dielectric constants, stripline dimensions, and frequencies.

Optimization of block-floating-point realizations for digital controllers with finite-word-length considerations.

PubMed

Wu, Jun; Hu, Xie-he; Chen, Sheng; Chu, Jian

2003-01-01

The closed-loop stability issue of finite-precision realizations was investigated for digital controllers implemented in block-floating-point format. The controller coefficient perturbation was analyzed resulting from using finite word length (FWL) block-floating-point representation scheme. A block-floating-point FWL closed-loop stability measure was derived which considers both the dynamic range and precision. To facilitate the design of optimal finite-precision controller realizations, a computationally tractable block-floating-point FWL closed-loop stability measure was then introduced and the method of computing the value of this measure for a given controller realization was developed. The optimal controller realization is defined as the solution that maximizes the corresponding measure, and a numerical optimization approach was adopted to solve the resulting optimal realization problem. A numerical example was used to illustrate the design procedure and to compare the optimal controller realization with the initial realization.

Accurate ω-ψ Spectral Solution of the Singular Driven Cavity Problem

NASA Astrophysics Data System (ADS)

Auteri, F.; Quartapelle, L.; Vigevano, L.

2002-08-01

This article provides accurate spectral solutions of the driven cavity problem, calculated in the vorticity-stream function representation without smoothing the corner singularities—a prima facie impossible task. As in a recent benchmark spectral calculation by primitive variables of Botella and Peyret, closed-form contributions of the singular solution for both zero and finite Reynolds numbers are subtracted from the unknown of the problem tackled here numerically in biharmonic form. The method employed is based on a split approach to the vorticity and stream function equations, a Galerkin-Legendre approximation of the problem for the perturbation, and an evaluation of the nonlinear terms by Gauss-Legendre numerical integration. Results computed for Re=0, 100, and 1000 compare well with the benchmark steady solutions provided by the aforementioned collocation-Chebyshev projection method. The validity of the proposed singularity subtraction scheme for computing time-dependent solutions is also established.

Structural robustness with suboptimal responses for linear state space model

NASA Technical Reports Server (NTRS)

Keel, L. H.; Lim, Kyong B.; Juang, Jer-Nan

1989-01-01

A relationship between the closed-loop eigenvalues and the amount of perturbations in the open-loop matrix is addressed in the context of performance robustness. If the allowable perturbation ranges of elements of the open-loop matrix A and the desired tolerance of the closed-loop eigenvalues are given such that max(j) of the absolute value of Delta-lambda(j) (A+BF) should be less than some prescribed value, what is a state feedback controller F which satisfies the closed-loop eigenvalue perturbation-tolerance requirement for a class of given perturbation in A? The paper gives an algorithm to design such a controller. Numerical examples are included for illustration.

Discrimination of time-dependent inflow properties with a cooperative dynamical system.

PubMed

Ueno, Hiroshi; Tsuruyama, Tatsuaki; Nowakowski, Bogdan; Górecki, Jerzy; Yoshikawa, Kenichi

2015-10-01

Many physical, chemical, and biological systems exhibit a cooperative or sigmoidal response with respect to the input. In biochemistry, such behavior is called an allosteric effect. Here, we demonstrate that a system with such properties can be used to discriminate the amplitude or frequency of an external periodic perturbation. Numerical simulations performed for a model sigmoidal kinetics illustrate that there exists a narrow range of frequencies and amplitudes within which the system evolves toward significantly different states. Therefore, observation of system evolution should provide information about the characteristics of the perturbation. The discrimination properties for periodic perturbation are generic. They can be observed in various dynamical systems and for different types of periodic perturbation.

Effects of Time-Dependent Inflow Perturbations on Turbulent Flow in a Street Canyon

NASA Astrophysics Data System (ADS)

Duan, G.; Ngan, K.

2017-12-01

Urban flow and turbulence are driven by atmospheric flows with larger horizontal scales. Since building-resolving computational fluid dynamics models typically employ steady Dirichlet boundary conditions or forcing, the accuracy of numerical simulations may be limited by the neglect of perturbations. We investigate the sensitivity of flow within a unit-aspect-ratio street canyon to time-dependent perturbations near the inflow boundary. Using large-eddy simulation, time-periodic perturbations to the streamwise velocity component are incorporated via the nudging technique. Spatial averages of pointwise differences between unperturbed and perturbed velocity fields (i.e., the error kinetic energy) show a clear dependence on the perturbation period, though spatial structures are largely insensitive to the time-dependent forcing. The response of the error kinetic energy is maximized for perturbation periods comparable to the time scale of the mean canyon circulation. Frequency spectra indicate that this behaviour arises from a resonance between the inflow forcing and the mean motion around closed streamlines. The robustness of the results is confirmed using perturbations derived from measurements of roof-level wind speed.

The comparison of stepping responses following perturbations applied to pelvis during overground and treadmill walking.

PubMed

Zadravec, Matjaž; Olenšek, Andrej; Matjačić, Zlatko

2017-08-09

Treadmills are used frequently in rehabilitation enabling neurologically impaired subjects to train walking while being assisted by therapists. Numerous studies compared walking on treadmill and overground for unperturbed but not also perturbed conditions. The objective of this study was to compare stepping responses (step length, step width and step time) during overground and treadmill walking in a group of healthy subjects where balance assessment robots applied perturbing pushes to the subject's pelvis in sagittal and frontal planes. During walking in both balance assessment robots (overground and treadmill-based) with applied perturbations the stepping responses of a group of seven healthy subjects were assessed with a motion tracking camera. The results show high degree of similarity of stepping responses between overground and treadmill walking for all perturbation directions. Both devices reproduced similar experimental conditions with relatively small standard deviations in the unperturbed walking as well as in perturbed walking. Based on these results we may conclude that stepping responses following perturbations can be studied on an instrumented treadmill where ground reaction forces can be readily assessed which is not the case during perturbed overground walking.

High sensitivity of tidewater outlet glacier dynamics to shape

NASA Astrophysics Data System (ADS)

Enderlin, E. M.; Howat, I. M.; Vieli, A.

2013-02-01

Variability in tidewater outlet glacier behavior under similar external forcing has been attributed to differences in outlet shape (i.e. bed elevation and width), but this dependence has not been investigated in detail. Here we use a numerical ice flow model to show that the dynamics of tidewater outlet glaciers under external forcing are highly sensitive to width and bed topography. Our sensitivity tests indicate that for glaciers with similar discharge, the trunks of wider glaciers and those grounded over deeper basal depressions tend to be closer to flotation, so that less dynamically induced thinning results in rapid, unstable retreat following a perturbation. The lag time between the onset of the perturbation and unstable retreat varies with outlet shape, which may help explain intra-regional variability in tidewater outlet glacier behavior. Further, because the perturbation response is dependent on the thickness relative to flotation, varying the bed topography within the range of observational uncertainty can result in either stable or unstable retreat due to the same perturbation. Thus, extreme care must be taken when interpreting the future behavior of actual glacier systems using numerical ice flow models that are not accompanied by comprehensive sensitivity analyses.

High sensitivity of tidewater outlet glacier dynamics to shape

NASA Astrophysics Data System (ADS)

Enderlin, E. M.; Howat, I. M.; Vieli, A.

2013-06-01

Variability in tidewater outlet glacier behavior under similar external forcing has been attributed to differences in outlet shape (i.e., bed elevation and width), but this dependence has not been investigated in detail. Here we use a numerical ice flow model to show that the dynamics of tidewater outlet glaciers under external forcing are highly sensitive to width and bed topography. Our sensitivity tests indicate that for glaciers with similar discharge, the trunks of wider glaciers and those grounded over deeper basal depressions tend to be closer to flotation, so that less dynamically induced thinning results in rapid, unstable retreat following a perturbation. The lag time between the onset of the perturbation and unstable retreat varies with outlet shape, which may help explain intra-regional variability in tidewater outlet glacier behavior. Further, because the perturbation response is dependent on the thickness relative to flotation, varying the bed topography within the range of observational uncertainty can result in either stable or unstable retreat due to the same perturbation. Thus, extreme care must be taken when interpreting the future behavior of actual glacier systems using numerical ice flow models that are not accompanied by comprehensive sensitivity analyses.

Optimal energy growth in a stably stratified shear flow

NASA Astrophysics Data System (ADS)

Jose, Sharath; Roy, Anubhab; Bale, Rahul; Iyer, Krithika; Govindarajan, Rama

2018-02-01

Transient growth of perturbations by a linear non-modal evolution is studied here in a stably stratified bounded Couette flow. The density stratification is linear. Classical inviscid stability theory states that a parallel shear flow is stable to exponentially growing disturbances if the Richardson number (Ri) is greater than 1/4 everywhere in the flow. Experiments and numerical simulations at higher Ri show however that algebraically growing disturbances can lead to transient amplification. The complexity of a stably stratified shear flow stems from its ability to combine this transient amplification with propagating internal gravity waves (IGWs). The optimal perturbations associated with maximum energy amplification are numerically obtained at intermediate Reynolds numbers. It is shown that in this wall-bounded flow, the three-dimensional optimal perturbations are oblique, unlike in unstratified flow. A partitioning of energy into kinetic and potential helps in understanding the exchange of energies and how it modifies the transient growth. We show that the apportionment between potential and kinetic energy depends, in an interesting manner, on the Richardson number, and on time, as the transient growth proceeds from an optimal perturbation. The oft-quoted stabilizing role of stratification is also probed in the non-diffusive limit in the context of disturbance energy amplification.

Estimation of metal strength at very high rates using free-surface Richtmyer–Meshkov Instabilities

DOE PAGES

Prime, Michael Bruce; Buttler, William Tillman; Buechler, Miles Allen; ...

2017-03-08

Recently, Richtmyer–Meshkov Instabilities (RMI) have been proposed for studying the average strength at strain rates up to at least 10 7/s. RMI experiments involve shocking a metal interface that has initial sinusoidal perturbations. The perturbations invert and grow subsequent to shock and may arrest because of strength effects. In this work we present new RMI experiments and data on a copper target that had five regions with different perturbation amplitudes on the free surface opposite the shock. We estimate the high-rate, low-pressure copper strength by comparing experimental data with Lagrangian numerical simulations. From a detailed computational study we find thatmore » mesh convergence must be carefully addressed to accurately compare with experiments, and numerical viscosity has a strong influence on convergence. We also find that modeling the as-built perturbation geometry rather than the nominal makes a significant difference. Because of the confounding effect of tensile damage on total spike growth, which has previously been used as the metric for estimating strength, we instead use a new strength metric: the peak velocity during spike growth. Furthermore, this new metric also allows us to analyze a broader set of experimental results that are sensitive to strength because some larger initial perturbations grow unstably to failure and so do not have a finite total spike growth.« less

Estimation of metal strength at very high rates using free-surface Richtmyer–Meshkov Instabilities

DOE Office of Scientific and Technical Information (OSTI.GOV)

Prime, Michael Bruce; Buttler, William Tillman; Buechler, Miles Allen

Recently, Richtmyer–Meshkov Instabilities (RMI) have been proposed for studying the average strength at strain rates up to at least 10 7/s. RMI experiments involve shocking a metal interface that has initial sinusoidal perturbations. The perturbations invert and grow subsequent to shock and may arrest because of strength effects. In this work we present new RMI experiments and data on a copper target that had five regions with different perturbation amplitudes on the free surface opposite the shock. We estimate the high-rate, low-pressure copper strength by comparing experimental data with Lagrangian numerical simulations. From a detailed computational study we find thatmore » mesh convergence must be carefully addressed to accurately compare with experiments, and numerical viscosity has a strong influence on convergence. We also find that modeling the as-built perturbation geometry rather than the nominal makes a significant difference. Because of the confounding effect of tensile damage on total spike growth, which has previously been used as the metric for estimating strength, we instead use a new strength metric: the peak velocity during spike growth. Furthermore, this new metric also allows us to analyze a broader set of experimental results that are sensitive to strength because some larger initial perturbations grow unstably to failure and so do not have a finite total spike growth.« less

Multi-symplectic integrators: numerical schemes for Hamiltonian PDEs that conserve symplecticity

NASA Astrophysics Data System (ADS)

Bridges, Thomas J.; Reich, Sebastian

2001-06-01

The symplectic numerical integration of finite-dimensional Hamiltonian systems is a well established subject and has led to a deeper understanding of existing methods as well as to the development of new very efficient and accurate schemes, e.g., for rigid body, constrained, and molecular dynamics. The numerical integration of infinite-dimensional Hamiltonian systems or Hamiltonian PDEs is much less explored. In this Letter, we suggest a new theoretical framework for generalizing symplectic numerical integrators for ODEs to Hamiltonian PDEs in R2: time plus one space dimension. The central idea is that symplecticity for Hamiltonian PDEs is directional: the symplectic structure of the PDE is decomposed into distinct components representing space and time independently. In this setting PDE integrators can be constructed by concatenating uni-directional ODE symplectic integrators. This suggests a natural definition of multi-symplectic integrator as a discretization that conserves a discrete version of the conservation of symplecticity for Hamiltonian PDEs. We show that this approach leads to a general framework for geometric numerical schemes for Hamiltonian PDEs, which have remarkable energy and momentum conservation properties. Generalizations, including development of higher-order methods, application to the Euler equations in fluid mechanics, application to perturbed systems, and extension to more than one space dimension are also discussed.

Fourier series expansion for nonlinear Hamiltonian oscillators.

PubMed

Méndez, Vicenç; Sans, Cristina; Campos, Daniel; Llopis, Isaac

2010-06-01

The problem of nonlinear Hamiltonian oscillators is one of the classical questions in physics. When an analytic solution is not possible, one can resort to obtaining a numerical solution or using perturbation theory around the linear problem. We apply the Fourier series expansion to find approximate solutions to the oscillator position as a function of time as well as the period-amplitude relationship. We compare our results with other recent approaches such as variational methods or heuristic approximations, in particular the Ren-He's method. Based on its application to the Duffing oscillator, the nonlinear pendulum and the eardrum equation, it is shown that the Fourier series expansion method is the most accurate.

Unsteady Blood Flow with Nanoparticles Through Stenosed Arteries in the Presence of Periodic Body Acceleration

NASA Astrophysics Data System (ADS)

Fatin Jamil, Dzuliana; Roslan, Rozaini; Abdulhameed, Mohammed; Che-Him, Norziha; Sufahani, Suliadi; Mohamad, Mahathir; Ghazali Kamardan, Muhamad

2018-04-01

The effects of nanoparticles such as Fe 3O4,TiO2, and Cu on blood flow inside a stenosed artery are studied. In this study, blood was modelled as non-Newtonian Bingham plastic fluid subjected to periodic body acceleration and slip velocity. The flow governing equations were solved analytically by using the perturbation method. By using the numerical approaches, the physiological parameters were analyzed, and the blood flow velocity distributions were generated graphically and discussed. From the flow results, the flow speed increases as slip velocity increases and decreases as the values of yield stress increases.

Development of a Perfectly Matched Layer Technique for a Discontinuous-Galerkin Spectral-Element Method

NASA Technical Reports Server (NTRS)

Garai, Anirban; Diosady, Laslo T.; Murman, Scott M.; Madavan, Nateri K.

2016-01-01

The perfectly matched layer (PML) technique is developed in the context of a high- order spectral-element Discontinuous-Galerkin (DG) method. The technique is applied to a range of test cases and is shown to be superior compared to other approaches, such as those based on using characteristic boundary conditions and sponge layers, for treating the inflow and outflow boundaries of computational domains. In general, the PML technique improves the quality of the numerical results for simulations of practical flow configurations, but it also exhibits some instabilities for large perturbations. A preliminary analysis that attempts to understand the source of these instabilities is discussed.

Local conservation laws and the structure of the many-body localized states.

PubMed

Serbyn, Maksym; Papić, Z; Abanin, Dmitry A

2013-09-20

We construct a complete set of local integrals of motion that characterize the many-body localized (MBL) phase. Our approach relies on the assumption that local perturbations act locally on the eigenstates in the MBL phase, which is supported by numerical simulations of the random-field XXZ spin chain. We describe the structure of the eigenstates in the MBL phase and discuss the implications of local conservation laws for its nonequilibrium quantum dynamics. We argue that the many-body localization can be used to protect coherence in the system by suppressing relaxation between eigenstates with different local integrals of motion.

Exact relativistic Toda chain eigenfunctions from Separation of Variables and gauge theory

NASA Astrophysics Data System (ADS)

Sciarappa, Antonio

2017-10-01

We provide a proposal, motivated by Separation of Variables and gauge theory arguments, for constructing exact solutions to the quantum Baxter equation associated to the N-particle relativistic Toda chain and test our proposal against numerical results. Quantum Mechanical non-perturbative corrections, essential in order to obtain a sensible solution, are taken into account in our gauge theory approach by considering codimension two defects on curved backgrounds (squashed S 5 and degenerate limits) rather than flat space; this setting also naturally incorporates exact quantization conditions and energy spectrum of the relativistic Toda chain as well as its modular dual structure.

Cloaking of arbitrarily shaped objects with homogeneous coatings

NASA Astrophysics Data System (ADS)

Forestiere, Carlo; Dal Negro, Luca; Miano, Giovanni

2014-05-01

We present a theory for the cloaking of arbitrarily shaped objects and demonstrate electromagnetic scattering cancellation through designed homogeneous coatings. First, in the small-particle limit, we expand the dipole moment of a coated object in terms of its resonant modes. By zeroing the numerator of the resulting rational function, we accurately predict the permittivity values of the coating layer that abates the total scattered power. Then, we extend the applicability of the method beyond the small-particle limit, deriving the radiation corrections of the scattering-cancellation permittivity within a perturbation approach. Our method permits the design of invisibility cloaks for irregularly shaped devices such as complex sensors and detectors.

Dynamic optimization of distributed biological systems using robust and efficient numerical techniques.

PubMed

Vilas, Carlos; Balsa-Canto, Eva; García, Maria-Sonia G; Banga, Julio R; Alonso, Antonio A

2012-07-02

Systems biology allows the analysis of biological systems behavior under different conditions through in silico experimentation. The possibility of perturbing biological systems in different manners calls for the design of perturbations to achieve particular goals. Examples would include, the design of a chemical stimulation to maximize the amplitude of a given cellular signal or to achieve a desired pattern in pattern formation systems, etc. Such design problems can be mathematically formulated as dynamic optimization problems which are particularly challenging when the system is described by partial differential equations.This work addresses the numerical solution of such dynamic optimization problems for spatially distributed biological systems. The usual nonlinear and large scale nature of the mathematical models related to this class of systems and the presence of constraints on the optimization problems, impose a number of difficulties, such as the presence of suboptimal solutions, which call for robust and efficient numerical techniques. Here, the use of a control vector parameterization approach combined with efficient and robust hybrid global optimization methods and a reduced order model methodology is proposed. The capabilities of this strategy are illustrated considering the solution of a two challenging problems: bacterial chemotaxis and the FitzHugh-Nagumo model. In the process of chemotaxis the objective was to efficiently compute the time-varying optimal concentration of chemotractant in one of the spatial boundaries in order to achieve predefined cell distribution profiles. Results are in agreement with those previously published in the literature. The FitzHugh-Nagumo problem is also efficiently solved and it illustrates very well how dynamic optimization may be used to force a system to evolve from an undesired to a desired pattern with a reduced number of actuators. The presented methodology can be used for the efficient dynamic optimization of generic distributed biological systems.

Elimination of secular terms from the differential equations for the elements of perturbed two-body motion

NASA Technical Reports Server (NTRS)

Bond, Victor R.; Fraietta, Michael F.

1991-01-01

In 1961, Sperling linearized and regularized the differential equations of motion of the two-body problem by changing the independent variable from time to fictitious time by Sundman's transformation (r = dt/ds) and by embedding the two-body energy integral and the Laplace vector. In 1968, Burdet developed a perturbation theory which was uniformly valid for all types of orbits using a variation of parameters approach on the elements which appeared in Sperling's equations for the two-body solution. In 1973, Bond and Hanssen improved Burdet's set of differential equations by embedding the total energy (which is a constant when the potential function is explicitly dependent upon time.) The Jacobian constant was used as an element to replace the total energy in a reformulation of the differential equations of motion. In the process, another element which is proportional to a component of the angular momentum was introduced. Recently trajectories computed during numerical studies of atmospheric entry from circular orbits and low thrust beginning in near-circular orbits exhibited numerical instability when solved by the method of Bond and Gottlieb (1989) for long time intervals. It was found that this instability was due to secular terms which appear on the righthand sides of the differential equations of some of the elements. In this paper, this instability is removed by the introduction of another vector integral called the delta integral (which replaces the Laplace Vector) and another scalar integral which removes the secular terms. The introduction of these integrals requires a new derivation of the differential equations for most of the elements. For this rederivation, the Lagrange method of variation of parameters is used, making the development more concise. Numerical examples of this improvement are presented.

Development of an efficient procedure for calculating the aerodynamic effects of planform variation

NASA Technical Reports Server (NTRS)

Mercer, J. E.; Geller, E. W.

1981-01-01

Numerical procedures to compute gradients in aerodynamic loading due to planform shape changes using panel method codes were studied. Two procedures were investigated: one computed the aerodynamic perturbation directly; the other computed the aerodynamic loading on the perturbed planform and on the base planform and then differenced these values to obtain the perturbation in loading. It is indicated that computing the perturbed values directly can not be done satisfactorily without proper aerodynamic representation of the pressure singularity at the leading edge of a thin wing. For the alternative procedure, a technique was developed which saves most of the time-consuming computations from a panel method calculation for the base planform. Using this procedure the perturbed loading can be calculated in about one-tenth the time of that for the base solution.

The spatial range of protein hydration

NASA Astrophysics Data System (ADS)

Persson, Filip; Söderhjelm, Pär; Halle, Bertil

2018-06-01

Proteins interact with their aqueous surroundings, thereby modifying the physical properties of the solvent. The extent of this perturbation has been investigated by numerous methods in the past half-century, but a consensus has still not emerged regarding the spatial range of the perturbation. To a large extent, the disparate views found in the current literature can be traced to the lack of a rigorous definition of the perturbation range. Stating that a particular solvent property differs from its bulk value at a certain distance from the protein is not particularly helpful since such findings depend on the sensitivity and precision of the technique used to probe the system. What is needed is a well-defined decay length, an intrinsic property of the protein in a dilute aqueous solution, that specifies the length scale on which a given physical property approaches its bulk-water value. Based on molecular dynamics simulations of four small globular proteins, we present such an analysis of the structural and dynamic properties of the hydrogen-bonded solvent network. The results demonstrate unequivocally that the solvent perturbation is short-ranged, with all investigated properties having exponential decay lengths of less than one hydration shell. The short range of the perturbation is a consequence of the high energy density of bulk water, rendering this solvent highly resistant to structural perturbations. The electric field from the protein, which under certain conditions can be long-ranged, induces a weak alignment of water dipoles, which, however, is merely the linear dielectric response of bulk water and, therefore, should not be thought of as a structural perturbation. By decomposing the first hydration shell into polarity-based subsets, we find that the hydration structure of the nonpolar parts of the protein surface is similar to that of small nonpolar solutes. For all four examined proteins, the mean number of water-water hydrogen bonds in the nonpolar subset is within 1% of the value in bulk water, suggesting that the fragmentation and topography of the nonpolar protein-water interface has evolved to minimize the propensity for protein aggregation by reducing the unfavorable free energy of hydrophobic hydration.

Improvement and performance evaluation of the perturbation source method for an exact Monte Carlo perturbation calculation in fixed source problems

NASA Astrophysics Data System (ADS)

Sakamoto, Hiroki; Yamamoto, Toshihiro

2017-09-01

This paper presents improvement and performance evaluation of the "perturbation source method", which is one of the Monte Carlo perturbation techniques. The formerly proposed perturbation source method was first-order accurate, although it is known that the method can be easily extended to an exact perturbation method. A transport equation for calculating an exact flux difference caused by a perturbation is solved. A perturbation particle representing a flux difference is explicitly transported in the perturbed system, instead of in the unperturbed system. The source term of the transport equation is defined by the unperturbed flux and the cross section (or optical parameter) changes. The unperturbed flux is provided by an "on-the-fly" technique during the course of the ordinary fixed source calculation for the unperturbed system. A set of perturbation particle is started at the collision point in the perturbed region and tracked until death. For a perturbation in a smaller portion of the whole domain, the efficiency of the perturbation source method can be improved by using a virtual scattering coefficient or cross section in the perturbed region, forcing collisions. Performance is evaluated by comparing the proposed method to other Monte Carlo perturbation methods. Numerical tests performed for a particle transport in a two-dimensional geometry reveal that the perturbation source method is less effective than the correlated sampling method for a perturbation in a larger portion of the whole domain. However, for a perturbation in a smaller portion, the perturbation source method outperforms the correlated sampling method. The efficiency depends strongly on the adjustment of the new virtual scattering coefficient or cross section.

A design approach for improving the performance of single-grid planar retarding potential analyzers

DOE Office of Scientific and Technical Information (OSTI.GOV)

Davidson, R. L.; Earle, G. D.

2011-01-15

Planar retarding potential analyzers (RPAs) have a long flight history and have been included on numerous spaceflight missions including Dynamics Explorer, the Defense Meteorological Satellite Program, and the Communications/Navigation Outage Forecast System. RPAs allow for simultaneous measurement of plasma composition, density, temperature, and the component of the velocity vector normal to the aperture plane. Internal conductive grids are used to approximate ideal potential planes within the instrument, but these grids introduce perturbations to the potential map inside the RPA and cause errors in the measurement of the parameters listed above. A numerical technique is presented herein for minimizing these gridmore » errors for a specific mission by varying the depth and spacing of the grid wires. The example mission selected concentrates on plasma dynamics near the sunset terminator in the equatorial region. The international reference ionosphere model is used to discern the average conditions expected for this mission, and a numerical model of the grid-particle interaction is used to choose a grid design that will best fulfill the mission goals.« less

Numerical analysis of mixing by sharp-edge-based acoustofluidic micromixer

NASA Astrophysics Data System (ADS)

Nama, Nitesh; Huang, Po-Hsun; Jun Huang, Tony; Costanzo, Francesco

2015-11-01

Recently, acoustically oscillated sharp-edges have been employed to realize rapid and homogeneous mixing at microscales (Huang, Lab on a Chip, 13, 2013). Here, we present a numerical model, qualitatively validated by experimental results, to analyze the acoustic mixing inside a sharp-edge-based micromixer. We extend our previous numerical model (Nama, Lab on a Chip, 14, 2014) to combine the Generalized Lagrangian Mean (GLM) theory with the convection-diffusion equation, while also allowing for the presence of a background flow as observed in a typical sharp-edge-based micromixer. We employ a perturbation approach to divide the flow variables into zeroth-, first- and second-order fields which are successively solved to obtain the Lagrangian mean velocity. The Langrangian mean velocity and the background flow velocity are further employed with the convection-diffusion equation to obtain the concentration profile. We characterize the effects of various operational and geometrical parameters to suggest potential design changes for improving the mixing performance of the sharp-edge-based micromixer. Lastly, we investigate the possibility of generation of a spatio-temporally controllable concentration gradient by placing sharp-edge structures inside the microchannel.

Transport and fluctuation-dissipation relations in asymptotic and preasymptotic diffusion across channels with variable section.

PubMed

Forte, Giuseppe; Cecconi, Fabio; Vulpiani, Angelo

2014-12-01

We study the asymptotic and preasymptotic diffusive properties of Brownian particles in channels whose section varies periodically in space. The effective diffusion coefficient D(eff) is numerically determined by the asymptotic behavior of the root mean square displacement in different geometries, considering even cases of steep variations of the channel boundaries. Moreover, we compared the numerical results to the predictions from the various corrections proposed in the literature to the well known Fick-Jacobs approximation. Building an effective one-dimensional equation for the longitudinal diffusion, we obtain an approximation for the effective diffusion coefficient. Such a result goes beyond a perturbation approach, and it is in good agreement with the actual values obtained by the numerical simulations. We discuss also the preasymptotic diffusion which is observed up to a crossover time whose value, in the presence of strong spatial variation of the channel cross section, can be very large. In addition, we show how the Einstein's relation between the mean drift induced by a small external field and the mean square displacement of the unperturbed system is valid in both asymptotic and preasymptotic regimes.

Development and evaluation of a hybrid averaged orbit generator

NASA Technical Reports Server (NTRS)

Mcclain, W. D.; Long, A. C.; Early, L. W.

1978-01-01

A rapid orbit generator based on a first-order application of the Generalized Method of Averaging has been developed for the Research and Development (R&D) version of the Goddard Trajectory Determination System (GTDS). The evaluation of the averaged equations of motion can use both numerically averaged and recursively evaluated, analytically averaged perturbation models. These equations are numerically integrated to obtain the secular and long-period motion. Factors affecting efficient orbit prediction are discussed and guidelines are presented for treatment of each major perturbation. Guidelines for obtaining initial mean elements compatible with the theory are presented. An overview of the orbit generator is presented and comparisons with high precision methods are given.

Coulomb double helical structure

NASA Astrophysics Data System (ADS)

Kamimura, Tetsuo; Ishihara, Osamu

2012-01-01

Structures of Coulomb clusters formed by dust particles in a plasma are studied by numerical simulation. Our study reveals the presence of various types of self-organized structures of a cluster confined in a prolate spheroidal electrostatic potential. The stable configurations depend on a prolateness parameter for the confining potential as well as on the number of dust particles in a cluster. One-dimensional string, two-dimensional zigzag structure and three-dimensional double helical structure are found as a result of the transition controlled by the prolateness parameter. The formation of stable double helical structures resulted from the transition associated with the instability of angular perturbations on double strings. Analytical perturbation study supports the findings of numerical simulations.

Numerical theory of the motion of Jupiter's Galilean satellites

NASA Astrophysics Data System (ADS)

Kosmodamianskii, G. A.

2009-12-01

A numerical theory of the motion of Jupiter’s Galilean satellites was constructed using 3767 absolute observations of the satellites. The theory was based on the numerical integration of the equations of motion of the satellites. The integration was carried out by Everhart’s method using the ERA software package developed at the Institute of Applied Astronomy (IAA). Perturbations due to the oblateness of the central planet, perturbations from Saturn and the Sun, and the mutual attraction of the satellites were taken into account in the integration. As a result, the coefficients of the Chebyshev series expansion for coordinates and velocities were found for the period from 1962 to 2010. The initial coordinates and velocities of the satellites, as well as their masses, the mass of Jupiter, and the harmonic coefficient J 2 of the potential of Jupiter, were adjusted. The resulting ephemerides were compared to those of Lieske and Lainey.

Constrained Perturbation Regularization Approach for Signal Estimation Using Random Matrix Theory

NASA Astrophysics Data System (ADS)

Suliman, Mohamed; Ballal, Tarig; Kammoun, Abla; Al-Naffouri, Tareq Y.

2016-12-01

In this supplementary appendix we provide proofs and additional extensive simulations that complement the analysis of the main paper (constrained perturbation regularization approach for signal estimation using random matrix theory).

Effect of error field correction coils on W7-X limiter loads

NASA Astrophysics Data System (ADS)

Bozhenkov, S. A.; Jakubowski, M. W.; Niemann, H.; Lazerson, S. A.; Wurden, G. A.; Biedermann, C.; Kocsis, G.; König, R.; Pisano, F.; Stephey, L.; Szepesi, T.; Wenzel, U.; Pedersen, T. S.; Wolf, R. C.; W7-X Team

2017-12-01

In the first campaign Wendelstein 7-X was operated with five poloidal graphite limiters installed stellarator symmetrically. In an ideal situation the power losses would be equally distributed between the limiters. The limiter shape was designed to smoothly distribute the heat flux over two strike lines. Vertically the strike lines are not uniform because of different connection lengths. In this paper it is demonstrated both numerically and experimentally that the heat flux distribution can be significantly changed by non-resonant n=1 perturbation field of the order of 10-4 . Numerical studies are performed with field line tracing. In experiments perturbation fields are excited with five error field trim coils. The limiters are diagnosed with infrared cameras, neutral gas pressure gauges, thermocouples and spectroscopic diagnostics. Experimental results are qualitatively consistent with the simulations. With a suitable choice of the phase and amplitude of the perturbation a more symmetric plasma-limiter interaction can be potentially achieved. These results are also of interest for the later W7-X divertor operation.

Comparison of the numerical modelling and experimental measurements of DIII-D separatrix displacements during H-modes with resonant magnetic perturbations

DOE PAGES

Orlov, Dmitry M.; Moyer, Richard A.; Evans, Todd E.; ...

2014-08-15

Numerical modeling of the plasma boundary position and its displacement due to external magnetic perturbations in DIII-D low-collisionality H-mode discharges is presented. The results of the vacuum model are compared to the experimental measurements for boundary displacements including Thomson scattering electron temperature T e, charge exchange recombination spectroscopy, beam emission spectroscopy, soft x-ray, and divertor Langmuir probe measurements. Magnetically perturbed discharges with toroidal mode number n=2 and n=3 are studied. It is shown that the vacuum model predictions agree well with the measurements above and below the midplane, and disagree at the outer midplane in discharges where significant kink amplificationmore » is present. Lastly, the role of the plasma response is studied using the two-fluid MHD code M3D-C 1, and the results are compared to the vacuum model showing that the plasma response model underestimates the boundary displacements.« less

Computational singular perturbation analysis of stochastic chemical systems with stiffness

DOE PAGES

Wang, Lijin; Han, Xiaoying; Cao, Yanzhao; ...

2017-01-25

Computational singular perturbation (CSP) is a useful method for analysis, reduction, and time integration of stiff ordinary differential equation systems. It has found dominant utility, in particular, in chemical reaction systems with a large range of time scales at continuum and deterministic level. On the other hand, CSP is not directly applicable to chemical reaction systems at micro or meso-scale, where stochasticity plays an non-negligible role and thus has to be taken into account. In this work we develop a novel stochastic computational singular perturbation (SCSP) analysis and time integration framework, and associated algorithm, that can be used to notmore » only construct accurately and efficiently the numerical solutions to stiff stochastic chemical reaction systems, but also analyze the dynamics of the reduced stochastic reaction systems. Furthermore, the algorithm is illustrated by an application to a benchmark stochastic differential equation model, and numerical experiments are carried out to demonstrate the effectiveness of the construction.« less

Gravity and large black holes in Randall-Sundrum II braneworlds.

PubMed

Figueras, Pau; Wiseman, Toby

2011-08-19

We show how to construct low energy solutions to the Randall-Sundrum II (RSII) model by using an associated five-dimensional anti-de Sitter space (AdS(5)) and/or four-dimensional conformal field theory (CFT(4)) problem. The RSII solution is given as a perturbation of the AdS(5)-CFT(4) solution, with the perturbation parameter being the radius of curvature of the brane metric compared to the AdS length ℓ. The brane metric is then a specific perturbation of the AdS(5)-CFT(4) boundary metric. For low curvatures the RSII solution reproduces 4D general relativity on the brane. Recently, AdS(5)-CFT(4) solutions with a 4D Schwarzschild boundary metric were numerically constructed. We modify the boundary conditions to numerically construct large RSII static black holes with radius up to ~20ℓ. For a large radius, the RSII solutions are indeed close to the associated AdS(5)-CFT(4) solution. © 2011 American Physical Society

The QBO and weak external forcing by solar activity: A three dimensional model study

NASA Technical Reports Server (NTRS)

Dameris, M.; Ebel, A.

1989-01-01

A better understanding is attempted of the physical mechanisms leading to significant correlations between oscillations in the lower and middle stratosphere and solar variability associated with the sun's rotation. A global 3-d mechanistic model of the middle atmosphere is employed to investigate the effects of minor artificially induced perturbations. The aim is to explore the physical mechanisms of the dynamical response especially of the stratosphere to weak external forcing as it may result from UV flux changes due to solar rotation. First results of numerical experiments dealing about the external forcing of the middle atmosphere by solar activity were presented elsewhere. Different numerical studies regarding the excitation and propagation of weak perturbations have been continued since then. The model calculations presented are made to investigate the influence of the quasi-biennial oscillation (QBO) on the dynamical response of the middle atmosphere to weak perturbations by employing different initial wind fields which represent the west and east phase of the QBO.

Impact of time-dependent nonaxisymmetric velocity perturbations on dynamo action of von Kármán-like flows.

PubMed

Giesecke, André; Stefani, Frank; Burguete, Javier

2012-12-01

We present numerical simulations of the kinematic induction equation in order to examine the dynamo efficiency of an axisymmetric von Kármán-like flow subject to time-dependent nonaxisymmetric velocity perturbations. The numerical model is based on the setup of the French von Kármán-sodium dynamo (VKS) and on the flow measurements from a water experiment conducted at the University of Navarra in Pamplona, Spain. The principal experimental observations that are modeled in our simulations are nonaxisymmetric vortexlike structures which perform an azimuthal drift motion in the equatorial plane. Our simulations show that the interactions of these periodic flow perturbations with the fundamental drift of the magnetic eigenmode (including the special case of nondrifting fields) essentially determine the temporal behavior of the dynamo state. We find two distinct regimes of dynamo action that depend on the (prescribed) drift frequency of an (m=2) vortexlike flow perturbation. For comparatively slowly drifting vortices we observe a narrow window with enhanced growth rates and a drift of the magnetic eigenmode that is synchronized with the perturbation drift. The resonance-like enhancement of the growth rates takes place when the vortex drift frequency roughly equals the drift frequency of the magnetic eigenmode in the unperturbed system. Outside of this small window, the field generation is hampered compared to the unperturbed case, and the field amplitude of the magnetic eigenmode is modulated with approximately twice the vortex drift frequency. The abrupt transition between the resonant regime and the modulated regime is identified as a spectral exceptional point where eigenvalues (growth rates and frequencies) and eigenfunctions of two previously independent modes collapse. In the actual configuration the drift frequencies of the velocity perturbations that are observed in the water experiment are much larger than the fundamental drift frequency of the magnetic eigenmode that is obtained from our numerical simulations. Hence, we conclude that the fulfillment of the resonance condition might be unlikely in present day dynamo experiments. However, a possibility to increase the dynamo efficiency in the VKS experiment might be realized by an application of holes or fingers on the outer boundary in the equatorial plane. These mechanical distortions provoke an anchorage of the vortices at fixed positions thus allowing an adjustment of the temporal behavior of the nonaxisymmetric flow perturbations.

Use of multivariable asymptotic expansions in a satellite theory

NASA Technical Reports Server (NTRS)

Dallas, S. S.

1973-01-01

Initial conditions and perturbative force of satellite are restricted to yield motion of equatorial satellite about oblate body. In this manner, exact analytic solution exists and can be used as standard of comparison in numerical accuracy comparisons. Detailed numerical accuracy studies of uniformly valid asymptotic expansions were made.

Breakdown of String Perturbation Theory for Many External Particles.

PubMed

Ghosh, Sudip; Raju, Suvrat

2017-03-31

We consider massless string scattering amplitudes in a limit where the number of external particles becomes very large, while the energy of each particle remains small. Using the growth of the volume of the relevant moduli space, and by means of independent numerical evidence, we argue that string perturbation theory breaks down in this limit. We discuss some remarkable implications for the information paradox.

On the account of gravitational perturbations in computer simulation technology of meteoroid complex formation and evolution

NASA Astrophysics Data System (ADS)

Kulikova, N. V.; Chepurova, V. M.

2009-10-01

So far we investigated the nonperturbation dynamics of meteoroid complexes. The numerical integration of the differential equations of motion in the N-body problem by the Everhart algorithm (N=2-6) and introduction of the intermediate hyperbolic orbits build on the base of the generalized problem of two fixed centers permit to take into account some gravitational perturbations.

New Approaches and Applications for Monte Carlo Perturbation Theory

DOE Office of Scientific and Technical Information (OSTI.GOV)

Aufiero, Manuele; Bidaud, Adrien; Kotlyar, Dan

2017-02-01

This paper presents some of the recent and new advancements in the extension of Monte Carlo Perturbation Theory methodologies and application. In particular, the discussed problems involve Brunup calculation, perturbation calculation based on continuous energy functions, and Monte Carlo Perturbation Theory in loosely coupled systems.

Cosmological perturbations during the Bose-Einstein condensation of dark matter

DOE Office of Scientific and Technical Information (OSTI.GOV)

Freitas, R.C.; Gonçalves, S.V.B., E-mail: rodolfo.camargo@pq.cnpq.br, E-mail: sergio.vitorino@pq.cnpq.br

In the present work, we analyze the evolution of the scalar and tensorial perturbations and the quantities relevant for the physical description of the Universe, as the density contrast of the scalar perturbations and the gravitational waves energy density during the Bose-Einstein condensation of dark matter. The behavior of these parameters during the Bose-Einstein phase transition of dark matter is analyzed in details. To study the cosmological dynamics and evolution of scalar and tensorial perturbations in a Universe with and without cosmological constant we use both analytical and numerical methods. The Bose-Einstein phase transition modifies the evolution of gravitational wavesmore » of cosmological origin, as well as the process of large-scale structure formation.« less

Manual physical therapy and perturbation exercises in knee osteoarthritis.

PubMed

Rhon, Daniel; Deyle, Gail; Gill, Norman; Rendeiro, Daniel

2013-11-01

Knee osteoarthritis (OA) causes disability among the elderly and is often associated with impaired balance and proprioception. Perturbation exercises may help improve these impairments. Although manual physical therapy is generally a well-tolerated treatment for knee OA, perturbation exercises have not been evaluated when used with a manual physical therapy approach. The purpose of this study was to observe tolerance to perturbation exercises and the effect of a manual physical therapy approach with perturbation exercises on patients with knee OA. This was a prospective observational cohort study of 15 patients with knee OA. The Western Ontario and McMaster Universities Arthritis Index (WOMAC), global rating of change (GROC), and 72-hour post-treatment tolerance were primary outcome measures. Patients received perturbation balance exercises along with a manual physical therapy approach, twice weekly for 4 weeks. Follow-up evaluation was done at 1, 3, and 6 months after beginning the program. Mean total WOMAC score significantly improved (P = 0.001) after the 4-week program (total WOMAC: initial, 105; 4 weeks, 56; 3 months, 54; 6 months, 57). Mean improvements were similar to previously published trials of manual physical therapy without perturbation exercises. The GROC score showed a minimal clinically important difference (MCID)≥+3 in 13 patients (87%) at 4 weeks, 12 patients (80%) at 3 months, and 9 patients (60%) at 6 months. No patients reported exacerbation of symptoms within 72 hours following each treatment session. A manual physical therapy approach that also included perturbation exercises was well tolerated and resulted in improved outcome scores in patients with knee OA.

Stochastic uncertainty analysis for unconfined flow systems

USGS Publications Warehouse

Liu, Gaisheng; Zhang, Dongxiao; Lu, Zhiming

2006-01-01

A new stochastic approach proposed by Zhang and Lu (2004), called the Karhunen‐Loeve decomposition‐based moment equation (KLME), has been extended to solving nonlinear, unconfined flow problems in randomly heterogeneous aquifers. This approach is on the basis of an innovative combination of Karhunen‐Loeve decomposition, polynomial expansion, and perturbation methods. The random log‐transformed hydraulic conductivity field (lnKS) is first expanded into a series in terms of orthogonal Gaussian standard random variables with their coefficients obtained as the eigenvalues and eigenfunctions of the covariance function of lnKS. Next, head h is decomposed as a perturbation expansion series Σh(m), where h(m) represents the mth‐order head term with respect to the standard deviation of lnKS. Then h(m) is further expanded into a polynomial series of m products of orthogonal Gaussian standard random variables whose coefficients hi1,i2,...,im(m) are deterministic and solved sequentially from low to high expansion orders using MODFLOW‐2000. Finally, the statistics of head and flux are computed using simple algebraic operations on hi1,i2,...,im(m). A series of numerical test results in 2‐D and 3‐D unconfined flow systems indicated that the KLME approach is effective in estimating the mean and (co)variance of both heads and fluxes and requires much less computational effort as compared to the traditional Monte Carlo simulation technique.

Improvements in mode-based waveform modeling and application to Eurasian velocity structure

NASA Astrophysics Data System (ADS)

Panning, M. P.; Marone, F.; Kim, A.; Capdeville, Y.; Cupillard, P.; Gung, Y.; Romanowicz, B.

2006-12-01

We introduce several recent improvements to mode-based 3D and asymptotic waveform modeling and examine how to integrate them with numerical approaches for an improved model of upper-mantle structure under eastern Eurasia. The first step in our approach is to create a large-scale starting model including shear anisotropy using Nonlinear Asymptotic Coupling Theory (NACT; Li and Romanowicz, 1995), which models the 2D sensitivity of the waveform to the great-circle path between source and receiver. We have recently improved this approach by implementing new crustal corrections which include a non-linear correction for the difference between the average structure of several large regions from the global model with further linear corrections to account for the local structure along the path between source and receiver (Marone and Romanowicz, 2006; Panning and Romanowicz, 2006). This model is further refined using a 3D implementation of Born scattering (Capdeville, 2005). We have made several recent improvements to this method, in particular introducing the ability to represent perturbations to discontinuities. While the approach treats all sensitivity as linear perturbations to the waveform, we have also experimented with a non-linear modification analogous to that used in the development of NACT. This allows us to treat large accumulated phase delays determined from a path-average approximation non-linearly, while still using the full 3D sensitivity of the Born approximation. Further refinement of shallow regions of the model is obtained using broadband forward finite-difference waveform modeling. We are also integrating a regional Spectral Element Method code into our tomographic modeling, allowing us to move beyond many assumptions inherent in the analytic mode-based approaches, while still taking advantage of their computational efficiency. Illustrations of the effects of these increasingly sophisticated steps will be presented.

The spatial sensitivity of Sp converted waves-kernels and their applications

NASA Astrophysics Data System (ADS)

Mancinelli, N. J.; Fischer, K. M.

2017-12-01

We have developed a framework for improved imaging of strong lateral variations in crust and upper mantle seismic discontinuity structure using teleseismic S-to-P (Sp) scattered waves. In our framework, we rapidly compute scattered wave sensitivities to velocity perturbations in a one-dimensional background model using ray-theoretical methods to account for timing, scattering, and geometrical spreading effects. The kernels accurately describe the amplitude and phase information of a scattered waveform, which we confirm by benchmarking against kernels derived from numerical solutions of the wave equation. The kernels demonstrate that the amplitude of an Sp converted wave at a given time is sensitive to structure along a quasi-hyperbolic curve, such that structure far from the direct ray path can influence the measurements. We use synthetic datasets to explore two potential applications of the scattered wave sensitivity kernels. First, we back-project scattered energy back to its origin using the kernel adjoint operator. This approach successfully images mantle interfaces at depths of 120-180 km with up to 20 km of vertical relief over lateral distances of 100 km (i.e., undulations with a maximal 20% grade) when station spacing is 10 km. Adjacent measurements sum coherently at nodes where gradients in seismic properties occur, and destructively interfere at nodes lacking gradients. In cases where the station spacing is greater than 10 km, the destructive interference can be incomplete, and smearing along the isochrons can occur. We demonstrate, however, that model smoothing can dampen these artifacts. This method is relatively fast, and accurately retrieves the positions of the interfaces, but it generally does not retrieve the strength of the velocity perturbations. Therefore, in our second approach, we attempt to invert directly for velocity perturbations from our reference model using an iterative conjugate-directions scheme.

Hybrid Evidence Theory-based Finite Element/Statistical Energy Analysis method for mid-frequency analysis of built-up systems with epistemic uncertainties

NASA Astrophysics Data System (ADS)

Yin, Shengwen; Yu, Dejie; Yin, Hui; Lü, Hui; Xia, Baizhan

2017-09-01

Considering the epistemic uncertainties within the hybrid Finite Element/Statistical Energy Analysis (FE/SEA) model when it is used for the response analysis of built-up systems in the mid-frequency range, the hybrid Evidence Theory-based Finite Element/Statistical Energy Analysis (ETFE/SEA) model is established by introducing the evidence theory. Based on the hybrid ETFE/SEA model and the sub-interval perturbation technique, the hybrid Sub-interval Perturbation and Evidence Theory-based Finite Element/Statistical Energy Analysis (SIP-ETFE/SEA) approach is proposed. In the hybrid ETFE/SEA model, the uncertainty in the SEA subsystem is modeled by a non-parametric ensemble, while the uncertainty in the FE subsystem is described by the focal element and basic probability assignment (BPA), and dealt with evidence theory. Within the hybrid SIP-ETFE/SEA approach, the mid-frequency response of interest, such as the ensemble average of the energy response and the cross-spectrum response, is calculated analytically by using the conventional hybrid FE/SEA method. Inspired by the probability theory, the intervals of the mean value, variance and cumulative distribution are used to describe the distribution characteristics of mid-frequency responses of built-up systems with epistemic uncertainties. In order to alleviate the computational burdens for the extreme value analysis, the sub-interval perturbation technique based on the first-order Taylor series expansion is used in ETFE/SEA model to acquire the lower and upper bounds of the mid-frequency responses over each focal element. Three numerical examples are given to illustrate the feasibility and effectiveness of the proposed method.

Drought on the North American High Plains: Modeling Effects of Vegetation, Temperature, and Rainfall Perturbations on Regional Hydrology

NASA Astrophysics Data System (ADS)

Hein, A. E.; Condon, L. E.; Maxwell, R. M.

2017-12-01

Large scale droughts can disrupt the water supply for agriculture, municipalities and industrial use worldwide. For example, the Dustbowl drought of the 1930s severely damaged agriculture on the North American High Plains. The Dustbowl is generally attributed to three major factors: increased temperature, decreased precipitation, and a change from native grasses that might have tolerated these climate perturbations to dryland wheat farming, which did not. This study explores the individual importance of each of these factors and the feedbacks between them. Previous modeling studies have explored how the High Plains system responds to changes in precipitation or temperature, but these models often depend on simplified or lumped parameter approaches. These approaches may not fully represent all the relevant physical processes, especially those related to energy balance changes due to increased temperature. For this study, we built a high-resolution model of the High Plains using ParFlow-CLM, an integrated hydrologic model that solves both energy and water balances from the subsurface to the top of vegetation. Model inputs including geology and climate forcing, together with representative precipitation and temperature changes for a major drought were assembled from public data. Numerical experiments were run to perturb vegetation, precipitation and temperature separately, as well as a baseline scenario with no changes and a worst-case scenario with all three simultaneously. The impact of each factor on High Plains hydrology and water resources was examined by comparing soil moisture, stream flow and water table levels between the runs. The one-factor experiments were used to show which of these outputs was the most sensitive and responded most quickly to each change. The worst-case scenario revealed interactions between the three factors.

Resonant frequency calculations using a hybrid perturbation-Galerkin technique

NASA Technical Reports Server (NTRS)

Geer, James F.; Andersen, Carl M.

1991-01-01

A two-step hybrid perturbation Galerkin technique is applied to the problem of determining the resonant frequencies of one or several degree of freedom nonlinear systems involving a parameter. In one step, the Lindstedt-Poincare method is used to determine perturbation solutions which are formally valid about one or more special values of the parameter (e.g., for large or small values of the parameter). In step two, a subset of the perturbation coordinate functions determined in step one is used in Galerkin type approximation. The technique is illustrated for several one degree of freedom systems, including the Duffing and van der Pol oscillators, as well as for the compound pendulum. For all of the examples considered, it is shown that the frequencies obtained by the hybrid technique using only a few terms from the perturbation solutions are significantly more accurate than the perturbation results on which they are based, and they compare very well with frequencies obtained by purely numerical methods.

Investigation of a Diagnostic for Perturbation Theory: Comparison to the T(sub 1) Diagnostic of Coupled-Cluster Theory

NASA Technical Reports Server (NTRS)

Lee, Timothy J.; Head-Gordon, Martin; Rendell, Alistair P.; Langhoff, Stephen R. (Technical Monitor)

1995-01-01

A diagnostic for perturbation theory calculations, S(sub 2), is defined and numerical results are compared to the established T(sub 1) diagnostic from coupled-cluster theory. S(sub 2) is the lowest order non-zero contribution to a perturbation expansion of T(sub 1). S(sub 2) is a reasonable estimate of the importance of non-dynamical electron correlation, although not as reliable as T(sub 1). S(sub 2) values less than or equal to 0.012 suggest that low orders of perturbation theory should yield reasonable results; S(sub 2) values between 0.012-0.015 suggest that caution is required in interpreting results from low orders of perturbation theory; S(sub 2) values greater than or equal to 0.015 indicate that low orders of perturbation theory are not reliable for accurate results. Although not required mathematically, S(sub 2) is always less than T(sub 1) for the examples studied here.

Tube dynamics and low energy Earth-Moon transfers in the 4-body system

NASA Astrophysics Data System (ADS)

Onozaki, Kaori; Yoshimura, Hiroaki; Ross, Shane D.

2017-11-01

In this paper, we show a low energy Earth-Moon transfer in the context of the Sun-Earth-Moon-spacecraft 4-body system. We consider the 4-body system as the coupled system of the Sun-Earth-spacecraft 3-body system perturbed by the Moon (which we call the Moon-perturbed system) and the Earth-Moon-spacecraft 3-body system perturbed by the Sun (which we call the Sun-perturbed system). In both perturbed systems, analogs of the stable and unstable manifolds are computed numerically by using the notion of Lagrangian coherent structures, wherein the stable and unstable manifolds play the role of separating orbits into transit and non-transit orbits. We obtain a family of non-transit orbits departing from a low Earth orbit in the Moon-perturbed system, and a family of transit orbits arriving into a low lunar orbit in the Sun-perturbed system. Finally, we show that we can construct a low energy transfer from the Earth to the Moon by choosing appropriate trajectories from both families and patching these trajectories with a maneuver.

Turbulence and mixing from optimal perturbations to a stratified shear layer

NASA Astrophysics Data System (ADS)

Kaminski, Alexis; Caulfield, C. P.; Taylor, John

2014-11-01

The stability and mixing of stratified shear layers is a canonical problem in fluid dynamics with relevance to flows in the ocean and atmosphere. The Miles-Howard theorem states that a necessary condition for normal-mode instability in parallel, inviscid, steady stratified shear flows is that the gradient Richardson number, Rig is less than 1/4 somewhere in the flow. However, substantial transient growth of non-normal modes may be possible at finite times even when Rig > 1 / 4 everywhere in the flow. We have calculated the ``optimal perturbations'' associated with maximum perturbation energy gain for a stably-stratified shear layer. These optimal perturbations are then used to initialize direct numerical simulations. For small but finite perturbation amplitudes, the optimal perturbations grow at the predicted linear rate initially, but then experience sufficient transient growth to become nonlinear and susceptible to secondary instabilities, which then break down into turbulence. Remarkably, this occurs even in flows for which Rig > 1 / 4 everywhere. We will describe the nonlinear evolution of the optimal perturbations and characterize the resulting turbulence and mixing.

A hybrid perturbation-Galerkin technique for partial differential equations

NASA Technical Reports Server (NTRS)

Geer, James F.; Anderson, Carl M.

1990-01-01

A two-step hybrid perturbation-Galerkin technique for improving the usefulness of perturbation solutions to partial differential equations which contain a parameter is presented and discussed. In the first step of the method, the leading terms in the asymptotic expansion(s) of the solution about one or more values of the perturbation parameter are obtained using standard perturbation methods. In the second step, the perturbation functions obtained in the first step are used as trial functions in a Bubnov-Galerkin approximation. This semi-analytical, semi-numerical hybrid technique appears to overcome some of the drawbacks of the perturbation and Galerkin methods when they are applied by themselves, while combining some of the good features of each. The technique is illustrated first by a simple example. It is then applied to the problem of determining the flow of a slightly compressible fluid past a circular cylinder and to the problem of determining the shape of a free surface due to a sink above the surface. Solutions obtained by the hybrid method are compared with other approximate solutions, and its possible application to certain problems associated with domain decomposition is discussed.

Effect of the pump rate and loss perturbations on the lasing dynamics of a Fabry-Perot laser

DOE Office of Scientific and Technical Information (OSTI.GOV)

Kumar, N; Ledenev, V I

2010-11-13

Transition from generation of the fundamental mode to generation of the fundamental and first modes is studied numerically under the action of nonstationary asymmetric perturbations of pump rate and loss distributions in the active medium layer. It is shown that emergence of perturbations directly leads to excitation of the first mode with significant amplitude. The regime of two-mode lasing in the presence of perturbations is shown to appear at a pump rate that is smaller than the threshold one for two-mode lasing in the absence of perturbations. It is found that the first-mode amplitude has a maximum at a frequencymore » of intermode beatings of an unfilled Fabry-Perot resonator. It is also determined that emergence of nonstationary asymmetric perturbations leads to an increase in the average intensity of the fundamental mode. Various transition regimes to two-mode lasing are compared in different types and periods of perturbations. The operability of the scheme controlling the mode composition of laser radiation is considered. (lasers)« less

Numerical Generation of Dense Plume Fingers in Unsaturated Homogeneous Porous Media

NASA Astrophysics Data System (ADS)

Cremer, C.; Graf, T.

2012-04-01

In nature, the migration of dense plumes typically results in the formation of vertical plume fingers. Flow direction in fingers is downwards, which is counterbalanced by upwards flow of less dense fluid between fingers. In heterogeneous media, heterogeneity itself is known to trigger the formation of fingers. In homogeneous media, however, fingers are also created even if all grains had the same diameter. The reason is that pore-scale heterogeneity leading to different flow velocities also exists in homogeneous media due to two effects: (i) Grains of identical size may randomly arrange differently, e.g. forming tetrahedrons, hexahedrons or octahedrons. Each arrangement creates pores of varying diameter, thus resulting in different average flow velocities. (ii) Random variations of solute concentration lead to varying buoyancy effects, thus also resulting in different velocities. As a continuation of previously made efforts to incorporate pore-scale heterogeneity into fully saturated soil such that dense fingers are realistically generated (Cremer and Graf, EGU Assembly, 2011), the current paper extends the research scope from saturated to unsaturated soil. Perturbation methods are evaluated by numerically re-simulating a laboratory-scale experiment of plume transport in homogeneous unsaturated sand (Simmons et al., Transp. Porous Media, 2002). The following 5 methods are being discussed: (i) homogeneous sand, (ii) initial perturbation of solute concentration, (iii) spatially random, time-constant perturbation of solute source, (iv) spatially and temporally random noise of simulated solute concentration, and (v) random K-field that introduces physically insignificant but numerically significant heterogeneity. Results demonstrate that, as opposed to saturated flow, perturbing the solute source will not result in plume fingering. This is because the location of the perturbed source (domain top) and the location of finger generation (groundwater surface) do not coincide. Alternatively, similar to saturated flow, applying either a random concentration noise (iv) or a random K-field (v) generates realistic plume fingering. Future work will focus on the generation mechanisms of plume finger splitting.

A Perturbation Analysis of Harmonics Generation from Saturated Elements in Power Systems

NASA Astrophysics Data System (ADS)

Kumano, Teruhisa

Nonlinear phenomena such as saturation in magnetic flux give considerable effects in power system analysis. It is reported that a failure in a real 500kV system triggered islanding operation, where resultant even harmonics caused malfunctions in protective relays. It is also reported that the major origin of this wave distortion is nothing but unidirectional magnetization of the transformer iron core. Time simulation is widely used today to analyze this type of phenomena, but it has basically two shortcomings. One is that the time simulation takes two much computing time in the vicinity of inflection points in the saturation characteristic curve because certain iterative procedure such as N-R (Newton-Raphson) should be used and such methods tend to be caught in an ill conditioned numerical hunting. The other is that such simulation methods sometimes do not help intuitive understanding of the studied phenomenon because the whole nonlinear equations are treated in a matrix form and not properly divided into understandable parts as done in linear systems. This paper proposes a new computation scheme which is based on so called perturbation method. Magnetic saturation in iron cores in a generator and a transformer are taken into account. The proposed method has a special feature against the first shortcoming of the N-R based time simulation method stated above. In the proposed method no iterative process is used to reduce the equation residue but uses perturbation series, which means free from the ill condition problem. Users have only to calculate each perturbation terms one by one until he reaches necessary accuracy. In a numerical example treated in the present paper the first order perturbation can make reasonably high accuracy, which means very fast computing. In numerical study three nonlinear elements are considered. Calculated results are almost identical to the conventional Newton-Raphson based time simulation, which shows the validity of the method. The proposed method would be effectively used in a screening where many case studies are needed.

Electron density increases due to Lightning activity as deduced from LWPC code and VLF signal perturbations.

NASA Astrophysics Data System (ADS)

Samir, Nait Amor; Bouderba, Yasmina

VLF signal perturbations in association with thunderstorm activity appear as changes in the signal amplitude and phase. Several papers reported on the characteristics of thus perturbations and their connection to the lightning strokes amplitude and polarity. In this contribution, we quantified the electrons density increases due to lightning activity by the use of the LWPC code and VLF signal perturbations parameters. The method is similar to what people did in studying the solar eruptions effect. the results showed that the reference height (h') decreased to lower altitudes (between 70 and 80 km). From the LWPC code results the maximum of the electron density was then deduced. Therefore, a numerical simulation of the atmospheric species times dependences was performed to study the recovery times of the electrons density at different heights. The results showed that the recovery time last for several minutes and explain the observation of long recovery Early signal perturbations.

Complete Hamiltonian analysis of cosmological perturbations at all orders

DOE Office of Scientific and Technical Information (OSTI.GOV)

Nandi, Debottam; Shankaranarayanan, S., E-mail: debottam@iisertvm.ac.in, E-mail: shanki@iisertvm.ac.in

2016-06-01

In this work, we present a consistent Hamiltonian analysis of cosmological perturbations at all orders. To make the procedure transparent, we consider a simple model and resolve the 'gauge-fixing' issues and extend the analysis to scalar field models and show that our approach can be applied to any order of perturbation for any first order derivative fields. In the case of Galilean scalar fields, our procedure can extract constrained relations at all orders in perturbations leading to the fact that there is no extra degrees of freedom due to the presence of higher time derivatives of the field in themore » Lagrangian. We compare and contrast our approach to the Lagrangian approach (Chen et al. [2006]) for extracting higher order correlations and show that our approach is efficient and robust and can be applied to any model of gravity and matter fields without invoking slow-roll approximation.« less

Adaptive Modeling Procedure Selection by Data Perturbation.

PubMed

Zhang, Yongli; Shen, Xiaotong

2015-10-01

Many procedures have been developed to deal with the high-dimensional problem that is emerging in various business and economics areas. To evaluate and compare these procedures, modeling uncertainty caused by model selection and parameter estimation has to be assessed and integrated into a modeling process. To do this, a data perturbation method estimates the modeling uncertainty inherited in a selection process by perturbing the data. Critical to data perturbation is the size of perturbation, as the perturbed data should resemble the original dataset. To account for the modeling uncertainty, we derive the optimal size of perturbation, which adapts to the data, the model space, and other relevant factors in the context of linear regression. On this basis, we develop an adaptive data-perturbation method that, unlike its nonadaptive counterpart, performs well in different situations. This leads to a data-adaptive model selection method. Both theoretical and numerical analysis suggest that the data-adaptive model selection method adapts to distinct situations in that it yields consistent model selection and optimal prediction, without knowing which situation exists a priori. The proposed method is applied to real data from the commodity market and outperforms its competitors in terms of price forecasting accuracy.

Maximally Rotating Supermassive Stars at the Onset of Collapse: The Perturbative Effects of Gas Pressure, Magnetic Fields, Dark Matter and Dark Energy

NASA Astrophysics Data System (ADS)

Butler, Satya P.; Lima, Alicia R.; Baumgarte, Thomas W.; Shapiro, Stuart L.

2018-04-01

The discovery of quasars at increasingly large cosmological redshifts may favor "direct collapse" as the most promising evolutionary route to the formation of supermassive black holes. In this scenario, supermassive black holes form when their progenitors - supermassive stars - become unstable to gravitational collapse. For uniformly rotating stars supported by pure radiation pressure and spinning at the mass-shedding limit, the critical configuration at the onset of collapse is characterized by universal values of the dimensionless spin and radius parameters J/M2 and R/M, independent of mass M. We consider perturbative effects of gas pressure, magnetic fields, dark matter and dark energy on these parameters, and thereby determine the domain of validity of this universality. We obtain leading-order corrections for the critical parameters and establish their scaling with the relevant physical parameters. We compare two different approaches to approximate the effects of gas pressure, which plays the most important role, find identical results for the above dimensionless parameters, and also find good agreement with recent numerical results.

Effect of large magnetic islands on screening of external magnetic perturbation fields at slow plasma flow

NASA Astrophysics Data System (ADS)

Li, L.; Liu, Y. Q.; Huang, X.; Luan, Q.; Zhong, F. C.

2017-02-01

A toroidal resistive magneto-hydrodynamic plasma response model, involving large magnetic islands, is proposed and numerically investigated, based on local flattening of the equilibrium pressure profile near a rational surface. It is assumed that such islands can be generated near the edge of the tokamak plasma, due to the penetration of the resonant magnetic perturbations, used for the purpose of controlling the edge localized mode. Within this model, it is found that the local flattening of the equilibrium pressure helps to mitigate the toroidal curvature induced screening effect [Glasser et al., Phys. Fluids 7, 875 (1975)]—the so called Glasser-Greene-Johnson screening, when the local toroidal flow near the mode rational surface is very slow (for example, as a result of mode locking associated with the field penetration). The saturation level of the plasma response amplitude is computed, as the plasma rotation frequency approaches zero. The local modification of the plasma resistivity inside the magnetic island is found to also affect the saturation level of the plasma response at vanishing flow.

Inflation with a smooth constant-roll to constant-roll era transition

NASA Astrophysics Data System (ADS)

Odintsov, S. D.; Oikonomou, V. K.

2017-07-01

In this paper, we study canonical scalar field models, with a varying second slow-roll parameter, that allow transitions between constant-roll eras. In the models with two constant-roll eras, it is possible to avoid fine-tunings in the initial conditions of the scalar field. We mainly focus on the stability of the resulting solutions, and we also investigate if these solutions are attractors of the cosmological system. We shall calculate the resulting scalar potential and, by using a numerical approach, we examine the stability and attractor properties of the solutions. As we show, the first constant-roll era is dynamically unstable towards linear perturbations, and the cosmological system is driven by the attractor solution to the final constant-roll era. As we demonstrate, it is possible to have a nearly scale-invariant power spectrum of primordial curvature perturbations in some cases; however, this is strongly model dependent and depends on the rate of the final constant-roll era. Finally, we present, in brief, the essential features of a model that allows oscillations between constant-roll eras.

Why does MP2 work?

PubMed

Fink, Reinhold F

2016-11-14

We show analytically and numerically that the performance of second order Møller-Plesset (MP) perturbation theory (PT), coupled-cluster (CC) theory, and other perturbation theory approaches can be rationalized by analyzing the wavefunctions of these methods. While rather large deviations for the individual contributions of configurations to the electron correlation energy are found for MP wavefunctions, they profit from an advantageous and robust error cancellation: The absolute contribution to the correlation energy is generally underestimated for the critical excitations with small energy denominators and all other doubly excited configurations where the two excited electrons are coupled to a singlet. This is balanced by an overestimation of the contribution of triplet-coupled double excitations to the correlation energy. The even better performance of spin-component-scaled-MP2 theory is explained by a similar error compensation effect. The wavefunction analysis for the lowest singlet states of H 2 O, CH 2 , CO, and Cu + shows the predicted trends for MP methods, rapid but biased convergence of CC theory as well as the substantial potential of linearized CC, or retaining the excitation-degree (RE)-PT.

TB3 - Measurement of vibrational-vibrational exchange of highly excited states of diatomic molecules where the collisional probability is approaching unity

DOE Office of Scientific and Technical Information (OSTI.GOV)

Nachshon, Y.; Coleman, P.

1975-08-01

An experimental method, employing a fast population perturbation technique, is described to measure the vibrational-vibrational (VV) collisional probability P/sub r,r-1/sup/v,v+1/ of a diatomic molecule for large vibrational quantum numbers r and v. The relaxation of the perturbed gain of a pair of vibrational levels is a function of the vibrational populations and VV rate constants k/sub r,r-1/sup v,v+1/. The numerical inversion of the VV master rate equations determining this relaxation does not give unique value for k/sub r,r-1/ sup v,v+1/ (or P/sub r,r-1/sup v,v+1), but lower bounds can be evaluated and with empirical formulas, having several adjustable constants, it canmore » be shown that probabilities of the order of unity are required to satisfy the experimental data. The method has been specifically applied to the CO molecule, but other molecules such as HX(X = F, Cl, Br), NO, etc., could also be measured.« less

Maximally rotating supermassive stars at the onset of collapse: the perturbative effects of gas pressure, magnetic fields, dark matter, and dark energy

NASA Astrophysics Data System (ADS)

Butler, Satya P.; Lima, Alicia R.; Baumgarte, Thomas W.; Shapiro, Stuart L.

2018-07-01

The discovery of quasars at increasingly large cosmological redshifts may favour `direct collapse' as the most promising evolutionary route to the formation of supermassive black holes. In this scenario, supermassive black holes form when their progenitors - supermassive stars - become unstable to gravitational collapse. For uniformly rotating stars supported by pure radiation pressure and spinning at the mass-shedding limit, the critical configuration at the onset of collapse is characterized by universal values of the dimensionless spin and radius parameters J/M2 and R/M, independent of mass M. We consider perturbative effects of gas pressure, magnetic fields, dark matter, and dark energy on these parameters, and thereby determine the domain of validity of this universality. We obtain leading-order corrections for the critical parameters and establish their scaling with the relevant physical parameters. We compare two different approaches to approximate the effects of gas pressure, which plays the most important role, find identical results for the above dimensionless parameters, and also find good agreement with recent numerical results.

Effects of Slag Ejection on Solid Rocket Motor Performance

NASA Technical Reports Server (NTRS)

Whitesides, R. Harold; Purinton, David C.; Hengel, John E.; Skelley, Stephen E.

1995-01-01

In past firings of the Reusable Solid Rocket Motor (RSRM) both static test and flight motors have shown small pressure perturbations occurring primarily between 65 and 80 seconds. A joint NASA/Thiokol team investigation concluded that the cause of the pressure perturbations was the periodic ingestion and ejection of molten aluminum oxide slag from the cavity around the submerged nozzle nose which tends to trap and collect individual aluminum oxide droplets from the approach flow. The conclusions of the team were supported by numerous data and observations from special tests including high speed photographic films, real time radiography, plume calorimeters, accelerometers, strain gauges, nozzle TVC system force gauges, and motor pressure and thrust data. A simplistic slag ballistics model was formulated to relate a given pressure perturbation to a required slag quantity. Also, a cold flow model using air and water was developed to provide data on the relationship between the slag flow rate and the chamber pressure increase. Both the motor and the cold flow model exhibited low frequency oscillations in conjunction with periods of slag ejection. Motor and model frequencies were related to scaling parameters. The data indicate that there is a periodicity to the slag entrainment and ejection phenomena which is possibly related to organized oscillations from instabilities in the dividing streamline shear layer which impinges on the underneath surface of the nozzle.

Estimate of Joule Heating in a Flat Dechirper

DOE Office of Scientific and Technical Information (OSTI.GOV)

Bane, Karl; Stupakov, Gennady; Gjonaj, Erion

2017-02-10

We have performed Joule power loss calculations for a flat dechirper. We have considered the configurations of the beam on-axis between the two plates—for chirp control—and for the beam especially close to one plate—for use as a fast kicker. Our calculations use a surface impedance approach, one that is valid when corrugation parameters are small compared to aperture (the perturbative parameter regime). In our model we ignore effects of field reflections at the sides of the dechirper plates, and thus expect the results to underestimate the Joule losses. The analytical results were also tested by numerical, time-domain simulations. We findmore » that most of the wake power lost by the beam is radiated out to the sides of the plates. For the case of the beam passing by a single plate, we derive an analytical expression for the broad-band impedance, and—in Appendix B—numerically confirm recently developed, analytical formulas for the short-range wakes. While our theory can be applied to the LCLS-II dechirper with large gaps, for the nominal apertures we are not in the perturbative regime and the reflection contribution to Joule losses is not negligible. With input from computer simulations, we estimate the Joule power loss (assuming bunch charge of 300 pC, repetition rate of 100 kHz) is 21 W/m for the case of two plates, and 24 W/m for the case of a single plate.« less

Predicting geomorphic evolution through integration of numerical-model scenarios and topographic/bathymetric-survey updates

NASA Astrophysics Data System (ADS)

Plant, N. G.; Long, J.; Dalyander, S.; Thompson, D.; Miselis, J. L.

2013-12-01

Natural resource and hazard management of barrier islands requires an understanding of geomorphic changes associated with long-term processes and storms. Uncertainty exists in understanding how long-term processes interact with the geomorphic changes caused by storms and the resulting perturbations of the long-term evolution trajectories. We use high-resolution data sets to initialize and correct high-fidelity numerical simulations of oceanographic forcing and resulting barrier island evolution. We simulate two years of observed storms to determine the individual and cumulative impacts of these events. Results are separated into cross-shore and alongshore components of sediment transport and compared with observed topographic and bathymetric changes during these time periods. The discrete island change induced by these storms is integrated with previous knowledge of long-term net alongshore sediment transport to project island evolution. The approach has been developed and tested using data collected at the Chandeleur Island chain off the coast of Louisiana (USA). The simulation time period included impacts from tropical and winter storms, as well as a human-induced perturbation associated with construction of a sand berm along the island shoreline. The predictions and observations indicated that storm and long-term processes both contribute to the migration, lowering, and disintegration of the artificial berm and natural island. Further analysis will determine the relative importance of cross-shore and alongshore sediment transport processes and the dominant time scales that drive each of these processes and subsequent island morphologic response.

A numerical study of atmospheric perturbations induced by heat from a wildland fire: sensitivity to vertical canopy structure and heat source strength

Treesearch

Michael T. Kiefer; Shiyuan Zhong; Warren E. Heilman; Joseph J. Charney; Xindi Bian

2018-01-01

An improved understanding of atmospheric perturbations within and above a forest during a wildland fire has relevance to many aspects of wildland fires including fire spread, smoke transport and dispersion, and tree mortality. In this study, the ARPS-CANOPY model, a version of the Advanced Regional Prediction System (ARPS) model with a canopy parameterization, is...

A perturbation analysis of a mechanical model for stable spatial patterning in embryology

NASA Astrophysics Data System (ADS)

Bentil, D. E.; Murray, J. D.

1992-12-01

We investigate a mechanical cell-traction mechanism that generates stationary spatial patterns. A linear analysis highlights the model's potential for these heterogeneous solutions. We use multiple-scale perturbation techniques to study the evolution of these solutions and compare our solutions with numerical simulations of the model system. We discuss some potential biological applications among which are the formation of ridge patterns, dermatoglyphs, and wound healing.

Viable tensor-to-scalar ratio in a symmetric matter bounce

NASA Astrophysics Data System (ADS)

Nath Raveendran, Rathul; Chowdhury, Debika; Sriramkumar, L.

2018-01-01

Matter bounces refer to scenarios wherein the universe contracts at early times as in a matter dominated epoch until the scale factor reaches a minimum, after which it starts expanding. While such scenarios are known to lead to scale invariant spectra of primordial perturbations after the bounce, the challenge has been to construct completely symmetric bounces that lead to a tensor-to-scalar ratio which is small enough to be consistent with the recent cosmological data. In this work, we construct a model involving two scalar fields (a canonical field and a non-canonical ghost field) to drive the symmetric matter bounce and study the evolution of the scalar perturbations in the model. We find that the model can be completely described in terms of a single parameter, viz. the ratio of the scale associated with the bounce to the value of the scale factor at the bounce. We evolve the scalar perturbations numerically across the bounce and evaluate the scalar power spectra after the bounce. We show that, while the scalar and tensor perturbation spectra are scale invariant over scales of cosmological interest, the tensor-to-scalar ratio proves to be much smaller than the current upper bound from the observations of the cosmic microwave background anisotropies by the Planck mission. We also support our numerical analysis with analytical arguments.

Resonant interaction of ϕ4 kink with PT-symmetric perturbation with spatially periodic gain/loss coefficient

NASA Astrophysics Data System (ADS)

Saadatmand, Danial; Borisov, Denis I.; Kevrekidis, Panayotis G.; Zhou, Kun; Dmitriev, Sergey V.

2018-03-01

The resonant interaction of the ϕ4 kink with a PT-symmetric perturbation is observed in the numerical study performed in the frame of the continuum model and with the help of a two degree of freedom collective variable model derived in PRA 89, 010102(R). The perturbation is in the form of first partial derivative in time term with a spatially periodic gain/loss coefficient. When the kink interacts with the perturbation, the kink's internal mode is excited with the amplitude varying in time quasiperiodically. The maximal value of the amplitude was found to grow when the kink velocity is such that it travels one period of the gain/loss prefactor in nearly one period of the kink's internal mode. It is also found that the kink's translational and vibrational modes are coupled in a way that an increase in the kink's internal mode amplitude results in a decrease in kink velocity. The results obtained with the collective variable method are in a good qualitative agreement with the numerical simulations for the continuum model. The results of the present study suggest that kink dynamics in open systems with balanced gain and loss can have new features in comparison with the case of conservative systems.

Numerical method for the solution of large systems of differential equations of the boundary layer type

NASA Technical Reports Server (NTRS)

Green, M. J.; Nachtsheim, P. R.

1972-01-01

A numerical method for the solution of large systems of nonlinear differential equations of the boundary-layer type is described. The method is a modification of the technique for satisfying asymptotic boundary conditions. The present method employs inverse interpolation instead of the Newton method to adjust the initial conditions of the related initial-value problem. This eliminates the so-called perturbation equations. The elimination of the perturbation equations not only reduces the user's preliminary work in the application of the method, but also reduces the number of time-consuming initial-value problems to be numerically solved at each iteration. For further ease of application, the solution of the overdetermined system for the unknown initial conditions is obtained automatically by applying Golub's linear least-squares algorithm. The relative ease of application of the proposed numerical method increases directly as the order of the differential-equation system increases. Hence, the method is especially attractive for the solution of large-order systems. After the method is described, it is applied to a fifth-order problem from boundary-layer theory.

Insights from Modeling the Integrated Climate, Biogeochemical Cycles, Human Activities and Their Interactions in the ACME Earth System Model

NASA Astrophysics Data System (ADS)

Leung, L. R.; Thornton, P. E.; Riley, W. J.; Calvin, K. V.

2017-12-01

Towards the goal of understanding the contributions from natural and managed systems to current and future greenhouse gas fluxes and carbon-climate and carbon-CO2 feedbacks, efforts have been underway to improve representations of the terrestrial, river, and human components of the ACME earth system model. Broadly, our efforts include implementation and comparison of approaches to represent the nutrient cycles and nutrient limitations on ecosystem production, extending the river transport model to represent sediment and riverine biogeochemistry, and coupling of human systems such as irrigation, reservoir operations, and energy and land use with the ACME land and river components. Numerical experiments have been designed to understand how terrestrial carbon, nitrogen, and phosphorus cycles regulate climate system feedbacks and the sensitivity of the feedbacks to different model treatments, examine key processes governing sediment and biogeochemistry in the rivers and their role in the carbon cycle, and exploring the impacts of human systems in perturbing the hydrological and carbon cycles and their interactions. This presentation will briefly introduce the ACME modeling approaches and discuss preliminary results and insights from numerical experiments that lay the foundation for improving understanding of the integrated climate-biogeochemistry-human system.

QmeQ 1.0: An open-source Python package for calculations of transport through quantum dot devices

NASA Astrophysics Data System (ADS)

Kiršanskas, Gediminas; Pedersen, Jonas Nyvold; Karlström, Olov; Leijnse, Martin; Wacker, Andreas

2017-12-01

QmeQ is an open-source Python package for numerical modeling of transport through quantum dot devices with strong electron-electron interactions using various approximate master equation approaches. The package provides a framework for calculating stationary particle or energy currents driven by differences in chemical potentials or temperatures between the leads which are tunnel coupled to the quantum dots. The electronic structures of the quantum dots are described by their single-particle states and the Coulomb matrix elements between the states. When transport is treated perturbatively to lowest order in the tunneling couplings, the possible approaches are Pauli (classical), first-order Redfield, and first-order von Neumann master equations, and a particular form of the Lindblad equation. When all processes involving two-particle excitations in the leads are of interest, the second-order von Neumann approach can be applied. All these approaches are implemented in QmeQ. We here give an overview of the basic structure of the package, give examples of transport calculations, and outline the range of applicability of the different approximate approaches.

Phases, periphases, and interphases equilibrium by molecular modeling. I. Mass equilibrium by the semianalytical stochastic perturbations method and application to a solution between (120) gypsum faces

NASA Astrophysics Data System (ADS)

Pedesseau, Laurent; Jouanna, Paul

2004-12-01

The SASP (semianalytical stochastic perturbations) method is an original mixed macro-nano-approach dedicated to the mass equilibrium of multispecies phases, periphases, and interphases. This general method, applied here to the reflexive relation Ck⇔μk between the concentrations Ck and the chemical potentials μk of k species within a fluid in equilibrium, leads to the distribution of the particles at the atomic scale. The macroaspects of the method, based on analytical Taylor's developments of chemical potentials, are intimately mixed with the nanoaspects of molecular mechanics computations on stochastically perturbed states. This numerical approach, directly linked to definitions, is universal by comparison with current approaches, DLVO Derjaguin-Landau-Verwey-Overbeek, grand canonical Monte Carlo, etc., without any restriction on the number of species, concentrations, or boundary conditions. The determination of the relation Ck⇔μk implies in fact two problems: a direct problem Ck⇒μk and an inverse problem μk⇒Ck. Validation of the method is demonstrated in case studies A and B which treat, respectively, a direct problem and an inverse problem within a free saturated gypsum solution. The flexibility of the method is illustrated in case study C dealing with an inverse problem within a solution interphase, confined between two (120) gypsum faces, remaining in connection with a reference solution. This last inverse problem leads to the mass equilibrium of ions and water molecules within a 3 Å thick gypsum interface. The major unexpected observation is the repulsion of SO42- ions towards the reference solution and the attraction of Ca2+ ions from the reference solution, the concentration being 50 times higher within the interphase as compared to the free solution. The SASP method is today the unique approach able to tackle the simulation of the number and distribution of ions plus water molecules in such extreme confined conditions. This result is of prime importance for all coupled chemical-mechanical problems dealing with interfaces, and more generally for a wide variety of applications such as phase changes, osmotic equilibrium, surface energy, etc., in complex chemical-physics situations.

A finite-element-based perturbation model for the rotordynamic analysis of shrouded pump impellers: Part 1: Model development and applications

NASA Technical Reports Server (NTRS)

Baskharone, Erian A.

1993-01-01

This study concerns the rotor dynamic characteristics of fluid-encompassed rotors, with special emphasis on shrouded pump impellers. The core of the study is a versatile and categorically new finite-element-based perturbation model, which is based on a rigorous flow analysis and what we have generically termed the 'virtually' deformable finite-element approach. The model is first applied to the case of a smooth annular seal for verification purposes. The rotor excitation components, in this sample problem, give rise to a purely cylindrical, purely conical, and a simultaneous cylindrical/conical rotor whirl around the housing centerline. In all cases, the computed results are compared to existing experimental and analytical data involving the same seal geometry and operating conditions. Next, two labyrinth-seal configurations, which share the same tooth-to-tooth chamber geometry but differ in the total number of chambers, were investigated. The results, in this case, are compared to experimental measurements for both seal configurations. The focus is finally shifted to the shrouded-impeller problem, where the stability effects of the leakage flow in the shroud-to-housing secondary passage are investigated. To this end, the computational model is applied to a typical shrouded-impeller pump stage, fabricated and rotor dynamically tested by Sulzer Bros., and the results compared to those of a simplified 'bulk-flow' analysis and Sulzer Bros.' test data. In addition to assessing the computed rotor dynamic coefficients, the shrouded-impeller study also covers a controversial topic, namely that of the leakage-passage inlet swirl, which was previously cited as the origin of highly unconventional (resonance-like) trends of the fluid-exerted forces. In order to validate this claim, a 'microscopic' study of the fluid/shroud interaction mechanism is conducted, with the focus being on the structure of the perturbed flow field associated with the impeller whirl. The conclusions of this study were solidified by the outcome of a numerical-certainty exercise, where the grid dependency of the numerical results is objectively examined. The final phase of the shrouded-impeller investigation involves the validation of a built-in assumption, in all other perturbation models, whereby single-harmonic tangential distributions of all the flow thermophysical properties are imposed. The last phase of the investigation course is aimed at verifying the fine details of the perturbed flow field in light of recent set of detailed LDA measurements in a smooth annular seal. Grid dependency of the fluid-induced forces is also investigated, and specific recommendations are made.

Pricing a raindrop in a process-based model: general methodology and a case study of the Upper-Zambezi

NASA Astrophysics Data System (ADS)

Albersen, Peter J.; Houba, Harold E. D.; Keyzer, Michiel A.

A general approach is presented to value the stocks and flows of water as well as the physical structure of the basin on the basis of an arbitrary process-based hydrological model. This approach adapts concepts from the economic theory of capital accumulation, which are based on Lagrange multipliers that reflect market prices in the absence of markets. This permits to derive a financial account complementing the water balance in which the value of deliveries by the hydrological system fully balances with the value of resources, including physical characteristics reflected in the shape of the functions in the model. The approach naturally suggests the use of numerical optimization software to compute the multipliers, without the need to impose an immensely large number of small perturbations on the simulation model, or to calculate all derivatives analytically. A novel procedure is proposed to circumvent numerical problems in computation and it is implemented in a numerical application using AQUA, an existing model of the Upper-Zambezi River. It appears, not unexpectedly, that most end value accrues to agriculture. Irrigated agriculture receives a remarkably large share, and is by far the most rewarding activity. Furthermore, according to the model, the economic value would be higher if temperature was lower, pointing to the detrimental effect of climate change. We also find that a significant economic value is stored in the groundwater stock because of its critical role in the dry season. As groundwater comes out as the main capital of the basin, its mining could be harmful.

Trigonometrically-fitted Scheifele two-step methods for perturbed oscillators

NASA Astrophysics Data System (ADS)

You, Xiong; Zhang, Yonghui; Zhao, Jinxi

2011-07-01

In this paper, a new family of trigonometrically-fitted Scheifele two-step (TFSTS) methods for the numerical integration of perturbed oscillators is proposed and investigated. An essential feature of TFSTS methods is that they are exact in both the internal stages and the updates when solving the unperturbed harmonic oscillator y″ = -ω2 y for known frequency ω. Based on the linear operator theory, the necessary and sufficient conditions for TFSTS methods of up to order five are derived. Two specific TFSTS methods of orders four and five respectively are constructed and their stability and phase properties are examined. In the five numerical experiments carried out the new integrators are shown to be more efficient and competent than some well-known methods in the literature.

δ M formalism: a new approach to cosmological perturbation theory in anisotropic inflation

NASA Astrophysics Data System (ADS)

Talebian-Ashkezari, A.; Ahmadi, N.; Abolhasani, A. A.

2018-03-01

We study the evolution of the metric perturbations in a Bianchi background in the long-wavelength limit. By applying the gradient expansion to the equations of motion we exhibit a generalized "Separate Universe" approach to the cosmological perturbation theory. Having found this consistent separate universe picture, we introduce the δ M formalism for calculating the evolution of the linear tensor perturbations in anisotropic inflation models in almost the same way that the so-called δ N formula is applied to the super-horizon dynamics of the curvature perturbations. Similar to her twin formula, δ N, this new method can substantially reduce the amount of calculations related to the evolution of tensor modes. However, it is not as general as δ N it is a "perturbative" formula and solves the shear only to linear order. In other words, it is restricted to weak shear limit.

A Double Perturbation Method for Reducing Dynamical Degradation of the Digital Baker Map

NASA Astrophysics Data System (ADS)

Liu, Lingfeng; Lin, Jun; Miao, Suoxia; Liu, Bocheng

2017-06-01

The digital Baker map is widely used in different kinds of cryptosystems, especially for image encryption. However, any chaotic map which is realized on the finite precision device (e.g. computer) will suffer from dynamical degradation, which refers to short cycle lengths, low complexity and strong correlations. In this paper, a novel double perturbation method is proposed for reducing the dynamical degradation of the digital Baker map. Both state variables and system parameters are perturbed by the digital logistic map. Numerical experiments show that the perturbed Baker map can achieve good statistical and cryptographic properties. Furthermore, a new image encryption algorithm is provided as a simple application. With a rather simple algorithm, the encrypted image can achieve high security, which is competitive to the recently proposed image encryption algorithms.

Thermal effects and sudden decay approximation in the curvaton scenario

DOE Office of Scientific and Technical Information (OSTI.GOV)

Kitajima, Naoya; Takesako, Tomohiro; Yokoyama, Shuichiro

2014-10-01

We study the impact of a temperature-dependent curvaton decay rate on the primordial curvature perturbation generated in the curvaton scenario. Using the familiar sudden decay approximation, we obtain an analytical expression for the curvature perturbation after the decay of the curvaton. We then investigate numerically the evolution of the background and of the perturbations during the decay. We first show that the instantaneous transfer coefficient, related to the curvaton energy fraction at the decay, can be extended into a more general parameter, which depends on the net transfer of the curvaton energy into radiation energy or, equivalently, on the totalmore » entropy ratio after the complete curvaton decay. We then compute the curvature perturbation and compare this result with the sudden decay approximation prediction.« less

A Systematic Approach for Computing Zero-Point Energy, Quantum Partition Function, and Tunneling Effect Based on Kleinert's Variational Perturbation Theory.

PubMed

Wong, Kin-Yiu; Gao, Jiali

2008-09-09

In this paper, we describe an automated integration-free path-integral (AIF-PI) method, based on Kleinert's variational perturbation (KP) theory, to treat internuclear quantum-statistical effects in molecular systems. We have developed an analytical method to obtain the centroid potential as a function of the variational parameter in the KP theory, which avoids numerical difficulties in path-integral Monte Carlo or molecular dynamics simulations, especially at the limit of zero-temperature. Consequently, the variational calculations using the KP theory can be efficiently carried out beyond the first order, i.e., the Giachetti-Tognetti-Feynman-Kleinert variational approach, for realistic chemical applications. By making use of the approximation of independent instantaneous normal modes (INM), the AIF-PI method can readily be applied to many-body systems. Previously, we have shown that in the INM approximation, the AIF-PI method is accurate for computing the quantum partition function of a water molecule (3 degrees of freedom) and the quantum correction factor for the collinear H(3) reaction rate (2 degrees of freedom). In this work, the accuracy and properties of the KP theory are further investigated by using the first three order perturbations on an asymmetric double-well potential, the bond vibrations of H(2), HF, and HCl represented by the Morse potential, and a proton-transfer barrier modeled by the Eckart potential. The zero-point energy, quantum partition function, and tunneling factor for these systems have been determined and are found to be in excellent agreement with the exact quantum results. Using our new analytical results at the zero-temperature limit, we show that the minimum value of the computed centroid potential in the KP theory is in excellent agreement with the ground state energy (zero-point energy) and the position of the centroid potential minimum is the expectation value of particle position in wave mechanics. The fast convergent property of the KP theory is further examined in comparison with results from the traditional Rayleigh-Ritz variational approach and Rayleigh-Schrödinger perturbation theory in wave mechanics. The present method can be used for thermodynamic and quantum dynamic calculations, including to systematically determine the exact value of zero-point energy and to study kinetic isotope effects for chemical reactions in solution and in enzymes.

Phases and stability of non-uniform black strings

NASA Astrophysics Data System (ADS)

Emparan, Roberto; Luna, Raimon; Martínez, Marina; Suzuki, Ryotaku; Tanabe, Kentaro

2018-05-01

We construct solutions of non-uniform black strings in dimensions from D ≈ 9 all the way up to D = ∞, and investigate their thermodynamics and dynamical stability. Our approach employs the large- D perturbative expansion beyond the leading order, including corrections up to 1 /D 4. Combining both analytical techniques and relatively simple numerical solution of ODEs, we map out the ranges of parameters in which non-uniform black strings exist in each dimension and compute their thermodynamics and quasinormal modes with accuracy. We establish with very good precision the existence of Sorkin's critical dimension and we prove that not only the thermodynamic stability, but also the dynamic stability of the solutions changes at it.

Dark soliton dynamics and interactions in continuous-wave-induced lattices.

PubMed

Tsopelas, Ilias; Kominis, Yannis; Hizanidis, Kyriakos

2007-10-01

The dynamics of dark spatial soliton beams and their interaction under the presence of a continuous wave (CW), which dynamically induces a photonic lattice, are investigated. It is shown that appropriate selection of the characteristic parameters of the CW result in controllable steering of a single soliton as well as controllable interaction between two solitons. Depending on the CW parameters, the soliton angle of propagation can be changed drastically, while two-soliton interaction can be either enhanced or reduced, suggesting a reconfigurable soliton control mechanism. Our analytical approach, based on the variational perturbation method, provides a dynamical system for the dark soliton evolution parameters. Analytical results are shown in good agreement with direct numerical simulations.

Transition of planar Couette flow at infinite Reynolds numbers.

PubMed

Itano, Tomoaki; Akinaga, Takeshi; Generalis, Sotos C; Sugihara-Seki, Masako

2013-11-01

An outline of the state space of planar Couette flow at high Reynolds numbers (Re<10^{5}) is investigated via a variety of efficient numerical techniques. It is verified from nonlinear analysis that the lower branch of the hairpin vortex state (HVS) asymptotically approaches the primary (laminar) state with increasing Re. It is also predicted that the lower branch of the HVS at high Re belongs to the stability boundary that initiates a transition to turbulence, and that one of the unstable manifolds of the lower branch of HVS lies on the boundary. These facts suggest HVS may provide a criterion to estimate a minimum perturbation arising transition to turbulent states at the infinite Re limit.

Gravitational self-force meets the post-Newtonian approximation in extreme-mass ratio inspiral of binary black holes

NASA Astrophysics Data System (ADS)

Detweiler, Steven

2010-02-01

Post-Newtonian analysis, numerical relativity and, now, perturbation-based gravitational self-force analysis are all being used to describe various aspects of black hole binary systems. Recent comparisons between self-force analysis, with m1m2, and post-Newtonian analysis, with v/c 1 show excellent agreement in their common domain of validity. This lends credence to the two very different regularization procedures which are invoked in these approximations. When self-force analysis is able to create gravitational waveforms from extreme mass-ratio inspiral, then unprecedented cross cultural comparisons of these three distinct approaches to understanding gravitational waves will reveal the strengths and weaknesses of each. )

Numerical modelling of distributed vibration sensor based on phase-sensitive OTDR

NASA Astrophysics Data System (ADS)

Masoudi, A.; Newson, T. P.

2017-04-01

A Distributed Vibration Sensor Based on Phase-Sensitive OTDR is numerically modeled. The advantage of modeling the building blocks of the sensor individually and combining the blocks to analyse the behavior of the sensing system is discussed. It is shown that the numerical model can accurately imitate the response of the experimental setup to dynamic perturbations a signal processing procedure similar to that used to extract the phase information from sensing setup.

Investigation of the local stress perturbation in Long Valley, California, by coupling seismic analyses and FEM numerical modeling

NASA Astrophysics Data System (ADS)

Lin, G.; Albino, F.; Amelung, F.

2017-12-01

Long Valley Caldera in eastern California is well known for producing numerous volcanic eruptions over the past 3 Myr. There has been a stress perturbation in the vicinity of the caldera with respect to the regional stress field. In this study, we combine seismic analyses and finite-element numerical modeling to investigate this local stress anomaly. We first compute focal mechanisms for earthquakes relocated by using a three-dimensional (3-D) seismic velocity model and waveform cross-correlation data. The final 42,000 good-quality focal solutions show that the mechanisms are dominated by approximately the same amount of normal faulting and strike-slip and much fewer reverse focal types. These focal mechanisms are then used to invert for the stress field in the study area by applying the SATSI algorithm. The orientations of the inverted minimum horizontal principal stress (ShMIN) greatly agree with those in previous studies based on analyses of focal mechanisms, borehole breakouts, and fault offsets. The NE-SW oriented ShMIN in the resurgent dome and south moat of the caldera is in contrast to the dominating E-W orientation in the western Basin and Range province and Mammoth Mountain. We then investigate which mechanism most likely causes this local stress perturbation by applying 3-D Finite Element Modeling (FEM). Mechanical properties (e.g., density, Poisson's ratio, and Young's Modulus) used in the model are derived from the latest 3-D seismic tomography model. Taking into account an initial stress field, we examine stress perturbations resulting from different sources: (1) pressurization of a magma reservoir, (2) dyking event, and (3) tectonic faulting; and compute the corresponding stress field orientation for each and compare it with the observations.

Extra-dimensional models on the lattice

DOE PAGES

Knechtli, Francesco; Rinaldi, Enrico

2016-08-05

In this paper we summarize the ongoing effort to study extra-dimensional gauge theories with lattice simulations. In these models the Higgs field is identified with extra-dimensional components of the gauge field. The Higgs potential is generated by quantum corrections and is protected from divergences by the higher dimensional gauge symmetry. Dimensional reduction to four dimensions can occur through compactification or localization. Gauge-Higgs unification models are often studied using perturbation theory. Numerical lattice simulations are used to go beyond these perturbative expectations and to include nonperturbative effects. We describe the known perturbative predictions and their fate in the strongly-coupled regime formore » various extra-dimensional models.« less

Application of the perturbation iteration method to boundary layer type problems.

PubMed

Pakdemirli, Mehmet

2016-01-01

The recently developed perturbation iteration method is applied to boundary layer type singular problems for the first time. As a preliminary work on the topic, the simplest algorithm of PIA(1,1) is employed in the calculations. Linear and nonlinear problems are solved to outline the basic ideas of the new solution technique. The inner and outer solutions are determined with the iteration algorithm and matched to construct a composite expansion valid within all parts of the domain. The solutions are contrasted with the available exact or numerical solutions. It is shown that the perturbation-iteration algorithm can be effectively used for solving boundary layer type problems.

Control optimization of a lifting body entry problem by an improved and a modified method of perturbation function. Ph.D. Thesis - Houston Univ.

NASA Technical Reports Server (NTRS)

Garcia, F., Jr.

1974-01-01

A study of the solution problem of a complex entry optimization was studied. The problem was transformed into a two-point boundary value problem by using classical calculus of variation methods. Two perturbation methods were devised. These methods attempted to desensitize the contingency of the solution of this type of problem on the required initial co-state estimates. Also numerical results are presented for the optimal solution resulting from a number of different initial co-states estimates. The perturbation methods were compared. It is found that they are an improvement over existing methods.

Effect of Body Perturbations on Hypersonic Flow Over Slender Power Law Bodies

NASA Technical Reports Server (NTRS)

Mirels, Harold; Thornton, Philip R.

1959-01-01

Hypersonic-slender-body theory, in the limit as the free-stream Mach number becomes infinite, is used to find the effect of slightly perturbing the surface of slender two-dimensional and axisymmetric power law bodies, The body perturbations are assumed to have a power law variation (with streamwise distance downstream of the nose of the body). Numerical results are presented for (1) the effect of boundary-layer development on two dimensional and axisymmetric bodies, (2) the effect of very small angles of attack (on tow[dimensional bodies), and (3) the effect of blunting the nose of very slender wedges and cones.

Numerical study of acoustophoretic motion of particles in a PDMS microchannel driven by surface acoustic waves.

PubMed

Nama, Nitesh; Barnkob, Rune; Mao, Zhangming; Kähler, Christian J; Costanzo, Francesco; Huang, Tony Jun

2015-06-21

We present a numerical study of the acoustophoretic motion of particles suspended in a liquid-filled PDMS microchannel on a lithium niobate substrate acoustically driven by surface acoustic waves. We employ a perturbation approach where the flow variables are divided into first- and second-order fields. We use impedance boundary conditions to model the PDMS microchannel walls and we model the acoustic actuation by a displacement function from the literature based on a numerical study of piezoelectric actuation. Consistent with the type of actuation, the obtained first-order field is a horizontal standing wave that travels vertically from the actuated wall towards the upper PDMS wall. This is in contrast to what is observed in bulk acoustic wave devices. The first-order fields drive the acoustic streaming, as well as the time-averaged acoustic radiation force acting on suspended particles. We analyze the motion of suspended particles driven by the acoustic streaming drag and the radiation force. We examine a range of particle diameters to demonstrate the transition from streaming-drag-dominated acoustophoresis to radiation-force-dominated acoustophoresis. Finally, as an application of our numerical model, we demonstrate the capability to tune the position of the vertical pressure node along the channel width by tuning the phase difference between two incoming surface acoustic waves.

Reduced conservatism in stability robustness bounds by state transformation

NASA Technical Reports Server (NTRS)

Yedavalli, R. K.; Liang, Z.

1986-01-01

This note addresses the issue of 'conservatism' in the time domain stability robustness bounds obtained by the Liapunov approach. A state transformation is employed to improve the upper bounds on the linear time-varying perturbation of an asymptotically stable linear time-invariant system for robust stability. This improvement is due to the variance of the conservatism of the Liapunov approach with respect to the basis of the vector space in which the Liapunov function is constructed. Improved bounds are obtained, using a transformation, on elemental and vector norms of perturbations (i.e., structured perturbations) as well as on a matrix norm of perturbations (i.e., unstructured perturbations). For the case of a diagonal transformation, an algorithm is proposed to find the 'optimal' transformation. Several examples are presented to illustrate the proposed analysis.

Desynchronization of stochastically synchronized chemical oscillators

DOE Office of Scientific and Technical Information (OSTI.GOV)

Snari, Razan; Tinsley, Mark R., E-mail: mark.tinsley@mail.wvu.edu, E-mail: kshowalt@wvu.edu; Faramarzi, Sadegh

Experimental and theoretical studies are presented on the design of perturbations that enhance desynchronization in populations of oscillators that are synchronized by periodic entrainment. A phase reduction approach is used to determine optimal perturbation timing based upon experimentally measured phase response curves. The effectiveness of the perturbation waveforms is tested experimentally in populations of periodically and stochastically synchronized chemical oscillators. The relevance of the approach to therapeutic methods for disrupting phase coherence in groups of stochastically synchronized neuronal oscillators is discussed.

Scalar Hairy Black Holes in Four Dimensions are Unstable

NASA Astrophysics Data System (ADS)

Ganchev, Bogdan; Santos, Jorge E.

2018-04-01

We present a numerical analysis of the stability properties of the black holes with scalar hair constructed by Herdeiro and Radu. We prove the existence of a novel gauge where the scalar field perturbations decouple from the metric perturbations, and analyze the resulting quasinormal mode spectrum. We find unstable modes with characteristic growth rates which for uniformly small hair are almost identical to those of a massive scalar field on a fixed Kerr background.

Scalar Hairy Black Holes in Four Dimensions are Unstable.

PubMed

Ganchev, Bogdan; Santos, Jorge E

2018-04-27

We present a numerical analysis of the stability properties of the black holes with scalar hair constructed by Herdeiro and Radu. We prove the existence of a novel gauge where the scalar field perturbations decouple from the metric perturbations, and analyze the resulting quasinormal mode spectrum. We find unstable modes with characteristic growth rates which for uniformly small hair are almost identical to those of a massive scalar field on a fixed Kerr background.

A new approach to the stability analysis of transient natural convection in porous media

NASA Astrophysics Data System (ADS)

Tilton, Nils

2016-11-01

Onset of natural convection due to transient diffusion in porous media has attracted considerable attention for its applications to CO2 sequestration. Stability analyses typically investigate onset of convection using an initial value problem approach in which a perturbation is introduced to the concentration field at an initial time t =tp . This leads to debate concerning physically appropriate perturbations, the critical time tc for linear instability, and to the counter-intuitive notion of an optimal initial time tp that maximizes perturbation growth. We propose a new approach in which transient diffusion is continuously perturbed by small variations in the porosity. With this approach, instability occurs immediately (tc = 0) without violating any physical constraints, such that the concepts of initial time tp and critical time tc have less relevance. We argue that the onset time for nonlinear convection is a more physically relevant parameter, and show that it can be predicted using a simple asymptotic expansion. Using the expansion, we consider porosity perturbations that vary sinusoidally in the horizontal and vertical directions, and show there are optimal combinations of wavelengths that minimize the onset time of nonlinear convection.

Time-dependent mass of cosmological perturbations in the hybrid and dressed metric approaches to loop quantum cosmology

NASA Astrophysics Data System (ADS)

Elizaga Navascués, Beatriz; Martín de Blas, Daniel; Mena Marugán, Guillermo A.

2018-02-01

Loop quantum cosmology has recently been applied in order to extend the analysis of primordial perturbations to the Planck era and discuss the possible effects of quantum geometry on the cosmic microwave background. Two approaches to loop quantum cosmology with admissible ultraviolet behavior leading to predictions that are compatible with observations are the so-called hybrid and dressed metric approaches. In spite of their similarities and relations, we show in this work that the effective equations that they provide for the evolution of the tensor and scalar perturbations are somewhat different. When backreaction is neglected, the discrepancy appears only in the time-dependent mass term of the corresponding field equations. We explain the origin of this difference, arising from the distinct quantization procedures. Besides, given the privileged role that the big bounce plays in loop quantum cosmology, e.g. as a natural instant of time to set initial conditions for the perturbations, we also analyze the positivity of the time-dependent mass when this bounce occurs. We prove that the mass of the tensor perturbations is positive in the hybrid approach when the kinetic contribution to the energy density of the inflaton dominates over its potential, as well as for a considerably large sector of backgrounds around that situation, while this mass is always nonpositive in the dressed metric approach. Similar results are demonstrated for the scalar perturbations in a sector of background solutions that includes the kinetically dominated ones; namely, the mass then is positive for the hybrid approach, whereas it typically becomes negative in the dressed metric case. More precisely, this last statement is strictly valid when the potential is quadratic for values of the inflaton mass that are phenomenologically favored.

Nonlocal transport in the presence of transport barriers

NASA Astrophysics Data System (ADS)

Del-Castillo-Negrete, D.

2013-10-01

There is experimental, numerical, and theoretical evidence that transport in plasmas can, under certain circumstances, depart from the standard local, diffusive description. Examples include fast pulse propagation phenomena in perturbative experiments, non-diffusive scaling in L-mode plasmas, and non-Gaussian statistics of fluctuations. From the theoretical perspective, non-diffusive transport descriptions follow from the relaxation of the restrictive assumptions (locality, scale separation, and Gaussian/Markovian statistics) at the foundation of diffusive models. We discuss an alternative class of models able to capture some of the observed non-diffusive transport phenomenology. The models are based on a class of nonlocal, integro-differential operators that provide a unifying framework to describe non- Fickian scale-free transport, and non-Markovian (memory) effects. We study the interplay between nonlocality and internal transport barriers (ITBs) in perturbative transport including cold edge pulses and power modulation. Of particular interest in the nonlocal ``tunnelling'' of perturbations through ITBs. Also, flux-gradient diagrams are discussed as diagnostics to detect nonlocal transport processes in numerical simulations and experiments. Work supported by the US Department of Energy.

Airfoil Design and Optimization by the One-Shot Method

NASA Technical Reports Server (NTRS)

Kuruvila, G.; Taasan, Shlomo; Salas, M. D.

1995-01-01

An efficient numerical approach for the design of optimal aerodynamic shapes is presented in this paper. The objective of any optimization problem is to find the optimum of a cost function subject to a certain state equation (governing equation of the flow field) and certain side constraints. As in classical optimal control methods, the present approach introduces a costate variable (Lagrange multiplier) to evaluate the gradient of the cost function. High efficiency in reaching the optimum solution is achieved by using a multigrid technique and updating the shape in a hierarchical manner such that smooth (low-frequency) changes are done separately from high-frequency changes. Thus, the design variables are changed on a grid where their changes produce nonsmooth (high-frequency) perturbations that can be damped efficiently by the multigrid. The cost of solving the optimization problem is approximately two to three times the cost of the equivalent analysis problem.

Quantification of uncertainty for fluid flow in heterogeneous petroleum reservoirs

NASA Astrophysics Data System (ADS)

Zhang, Dongxiao

Detailed description of the heterogeneity of oil/gas reservoirs is needed to make performance predictions of oil/gas recovery. However, only limited measurements at a few locations are usually available. This combination of significant spatial heterogeneity with incomplete information about it leads to uncertainty about the values of reservoir properties and thus, to uncertainty in estimates of production potential. The theory of stochastic processes provides a natural method for evaluating these uncertainties. In this study, we present a stochastic analysis of transient, single phase flow in heterogeneous reservoirs. We derive general equations governing the statistical moments of flow quantities by perturbation expansions. These moments can be used to construct confidence intervals for the flow quantities (e.g., pressure and flow rate). The moment equations are deterministic and can be solved numerically with existing solvers. The proposed moment equation approach has certain advantages over the commonly used Monte Carlo approach.

Numerical Prediction of Combustion-induced Noise using a hybrid LES/CAA approach

NASA Astrophysics Data System (ADS)

Ihme, Matthias; Pitsch, Heinz; Kaltenbacher, Manfred

2006-11-01

Noise generation in technical devices is an increasingly important problem. Jet engines in particular produce sound levels that not only are a nuisance but may also impair hearing. The noise emitted by such engines is generated by different sources such as jet exhaust, fans or turbines, and combustion. Whereas the former acoustic mechanisms are reasonably well understood, combustion-generated noise is not. A methodology for the prediction of combustion-generated noise is developed. In this hybrid approach unsteady acoustic source terms are obtained from an LES and the propagation of pressure perturbations are obtained using acoustic analogies. Lighthill's acoustic analogy and a non-linear wave equation, accounting for variable speed of sound, have been employed. Both models are applied to an open diffusion flame. The effects on the far field pressure and directivity due to the variation of speed of sound are analyzed. Results for the sound pressure level will be compared with experimental data.

Airfoil optimization by the one-shot method

NASA Technical Reports Server (NTRS)

Kuruvila, G.; Taasan, Shlomo; Salas, M. D.

1994-01-01

An efficient numerical approach for the design of optimal aerodynamic shapes is presented in this paper. The objective of any optimization problem is to find the optimum of a cost function subject to a certain state equation (Governing equation of the flow field) and certain side constraints. As in classical optimal control methods, the present approach introduces a costate variable (Language multiplier) to evaluate the gradient of the cost function. High efficiency in reaching the optimum solution is achieved by using a multigrid technique and updating the shape in a hierarchical manner such that smooth (low-frequency) changes are done separately from high-frequency changes. Thus, the design variables are changed on a grid where their changes produce nonsmooth (high-frequency) perturbations that can be damped efficiently by the multigrid. The cost of solving the optimization problem is approximately two to three times the cost of the equivalent analysis problem.

Aerodynamic design and optimization in one shot

NASA Technical Reports Server (NTRS)

Ta'asan, Shlomo; Kuruvila, G.; Salas, M. D.

1992-01-01

This paper describes an efficient numerical approach for the design and optimization of aerodynamic bodies. As in classical optimal control methods, the present approach introduces a cost function and a costate variable (Lagrange multiplier) in order to achieve a minimum. High efficiency is achieved by using a multigrid technique to solve for all the unknowns simultaneously, but restricting work on a design variable only to grids on which their changes produce nonsmooth perturbations. Thus, the effort required to evaluate design variables that have nonlocal effects on the solution is confined to the coarse grids. However, if a variable has a nonsmooth local effect on the solution in some neighborhood, it is relaxed in that neighborhood on finer grids. The cost of solving the optimal control problem is shown to be approximately two to three times the cost of the equivalent analysis problem. Examples are presented to illustrate the application of the method to aerodynamic design and constraint optimization.

New imaging algorithm in diffusion tomography

NASA Astrophysics Data System (ADS)

Klibanov, Michael V.; Lucas, Thomas R.; Frank, Robert M.

1997-08-01

A novel imaging algorithm for diffusion/optical tomography is presented for the case of the time dependent diffusion equation. Numerical tests are conducted for ranges of parameters realistic for applications to an early breast cancer diagnosis using ultrafast laser pulses. This is a perturbation-like method which works for both homogeneous a heterogeneous background media. Its main innovation lies in a new approach for a novel linearized problem (LP). Such an LP is derived and reduced to a boundary value problem for a coupled system of elliptic partial differential equations. As is well known, the solution of such a system amounts to the factorization of well conditioned, sparse matrices with few non-zero entries clustered along the diagonal, which can be done very rapidly. Thus, the main advantages of this technique are that it is fast and accurate. The authors call this approach the elliptic systems method (ESM). The ESM can be extended for other data collection schemes.

Perturbative and Ab-Initio Calculations of Electrical Susceptibilities of Atoms

NASA Astrophysics Data System (ADS)

Spott, Andrew

Perturbative nonlinear optics consists of many powerful predictive theoretical methods, including the perturbative series of observables related to the interaction of light with matter. The light intensity limits of such series have been studied in the past for highly nonlinear processes such as above threshold ionization and high harmonic generation. A more recent debate focuses on the limits of applicability of perturbation theory for the nonlinear electrical susceptibility and the nonlinear index of refraction of atoms, which are important parameters to study, for example, for filamentation of laser pulses in nonlinear media. In this thesis we analyze theoretical predictions for the electrical susceptibility of atoms for the transition from the perturbative to the nonperturbative intensity regime. To this end, we apply a numerical basis state method that allows us to perform respective calculations in the framework of perturbation theory as well as using ab-initio methods. The results let us identify the intensity at which the application of perturbation theory breaks down. Furthermore, we provide an analysis of the nonlinear susceptibility as a function of time during the interaction with the laser pulse and find that theoretical predictions are in good agreement with recent experimental data.

Nonlinear power spectrum from resummed perturbation theory: a leap beyond the BAO scale

DOE Office of Scientific and Technical Information (OSTI.GOV)

Anselmi, Stefano; Pietroni, Massimo, E-mail: anselmi@ieec.uab.es, E-mail: massimo.pietroni@pd.infn.it

2012-12-01

A new computational scheme for the nonlinear cosmological matter power spectrum (PS) is presented. Our method is based on evolution equations in time, which can be cast in a form extremely convenient for fast numerical evaluations. A nonlinear PS is obtained in a time comparable to that needed for a simple 1-loop computation, and the numerical implementation is very simple. Our results agree with N-body simulations at the percent level in the BAO range of scales, and at the few-percent level up to k ≅ 1 h/Mpc at z∼>0.5, thereby opening the possibility of applying this tool to scales interestingmore » for weak lensing. We clarify the approximations inherent to this approach as well as its relations to previous ones, such as the Time Renormalization Group, and the multi-point propagator expansion. We discuss possible lines of improvements of the method and its intrinsic limitations by multi streaming at small scales and low redshifts.« less

Analytical approximation for the Einstein-dilaton-Gauss-Bonnet black hole metric

NASA Astrophysics Data System (ADS)

Kokkotas, K. D.; Konoplya, R. A.; Zhidenko, A.

2017-09-01

We construct an analytical approximation for the numerical black hole metric of P. Kanti et al. [Phys. Rev. D 54, 5049 (1996), 10.1103/PhysRevD.54.5049] in the four-dimensional Einstein-dilaton-Gauss-Bonnet (EdGB) theory. The continued fraction expansion in terms of a compactified radial coordinate, used here, converges slowly when the dilaton coupling approaches its extremal values, but for a black hole far from the extremal state, the analytical formula has a maximal relative error of a fraction of one percent already within the third order of the continued fraction expansion. The suggested analytical representation of the numerical black hole metric is relatively compact and a good approximation in the whole space outside the black hole event horizon. Therefore, it can serve in the same way as an exact solution when analyzing particles' motion, perturbations, quasinormal modes, Hawking radiation, accreting disks, and many other problems in the vicinity of a black hole. In addition, we construct the approximate analytical expression for the dilaton field.

Non-Schwarzschild black-hole metric in four dimensional higher derivative gravity: Analytical approximation

NASA Astrophysics Data System (ADS)

Kokkotas, K. D.; Konoplya, R. A.; Zhidenko, A.

2017-09-01

Higher derivative extensions of Einstein gravity are important within the string theory approach to gravity and as alternative and effective theories of gravity. H. Lü, A. Perkins, C. Pope, and K. Stelle [Phys. Rev. Lett. 114, 171601 (2015), 10.1103/PhysRevLett.114.171601] found a numerical solution describing a spherically symmetric non-Schwarzschild asymptotically flat black hole in Einstein gravity with added higher derivative terms. Using the general and quickly convergent parametrization in terms of the continued fractions, we represent this numerical solution in the analytical form, which is accurate not only near the event horizon or far from the black hole, but in the whole space. Thereby, the obtained analytical form of the metric allows one to study easily all the further properties of the black hole, such as thermodynamics, Hawking radiation, particle motion, accretion, perturbations, stability, quasinormal spectrum, etc. Thus, the found analytical approximate representation can serve in the same way as an exact solution.

Piezo-optic tensor of crystals from quantum-mechanical calculations.

PubMed

Erba, A; Ruggiero, M T; Korter, T M; Dovesi, R

2015-10-14

An automated computational strategy is devised for the ab initio determination of the full fourth-rank piezo-optic tensor of crystals belonging to any space group of symmetry. Elastic stiffness and compliance constants are obtained as numerical first derivatives of analytical energy gradients with respect to the strain and photo-elastic constants as numerical derivatives of analytical dielectric tensor components, which are in turn computed through a Coupled-Perturbed-Hartree-Fock/Kohn-Sham approach, with respect to the strain. Both point and translation symmetries are exploited at all steps of the calculation, within the framework of periodic boundary conditions. The scheme is applied to the determination of the full set of ten symmetry-independent piezo-optic constants of calcium tungstate CaWO4, which have recently been experimentally reconstructed. Present calculations unambiguously determine the absolute sign (positive) of the π61 constant, confirm the reliability of 6 out of 10 experimentally determined constants and provide new, more accurate values for the remaining 4 constants.

Piezo-optic tensor of crystals from quantum-mechanical calculations

DOE Office of Scientific and Technical Information (OSTI.GOV)

Erba, A., E-mail: alessandro.erba@unito.it; Dovesi, R.; Ruggiero, M. T.

2015-10-14

An automated computational strategy is devised for the ab initio determination of the full fourth-rank piezo-optic tensor of crystals belonging to any space group of symmetry. Elastic stiffness and compliance constants are obtained as numerical first derivatives of analytical energy gradients with respect to the strain and photo-elastic constants as numerical derivatives of analytical dielectric tensor components, which are in turn computed through a Coupled-Perturbed-Hartree-Fock/Kohn-Sham approach, with respect to the strain. Both point and translation symmetries are exploited at all steps of the calculation, within the framework of periodic boundary conditions. The scheme is applied to the determination of themore » full set of ten symmetry-independent piezo-optic constants of calcium tungstate CaWO{sub 4}, which have recently been experimentally reconstructed. Present calculations unambiguously determine the absolute sign (positive) of the π{sub 61} constant, confirm the reliability of 6 out of 10 experimentally determined constants and provide new, more accurate values for the remaining 4 constants.« less

Atomic mean-square displacement of a solid: A Green's-function approach

NASA Astrophysics Data System (ADS)

Shukla, R. C.; Hübschle, Hermann

1989-07-01

We have presented a Green's-function method of calculating the atomic mean-square displacement (MSD) of a solid. The method effectively sums a class of all anharmonic contributions to the MSD. From the point of view of perturbation theory (PT) our expression for MSD includes the lowest-order (λ2) PT contributions (cubic and quartic) with correct numerical coefficients. The numerical results obtained by this method in the high-temperature limit for a fcc nearest-neighbor Lennard-Jones-solid model are in excellent agreement with the Monte Carlo (MC) method for the same model over the entire temperature range of the solid. Highly accurate results for the order-λ2 PT contributions to MSD are obtained by eliminating the uncertainty in the convergence of the cubic contributions in the earlier work of Heiser, Shukla, and Cowly and they are now in much better agreement with the MC results but still inferior to the Green's-function method at the highest temperature.

Protecting a quantum state from environmental noise by an incompatible finite-time measurement

DOE Office of Scientific and Technical Information (OSTI.GOV)

Brasil, Carlos Alexandre; Castro, L. A. de; Napolitano, R. d. J.

We show that measurements of finite duration performed on an open two-state system can protect the initial state from a phase-noisy environment, provided the measured observable does not commute with the perturbing interaction. When the measured observable commutes with the environmental interaction, the finite-duration measurement accelerates the rate of decoherence induced by the phase noise. For the description of the measurement of an observable that is incompatible with the interaction between system and environment, we have found an approximate analytical expression, valid at zero temperature and weak coupling with the measuring device. We have tested the validity of the analyticalmore » predictions against an exact numerical approach, based on the superoperator-splitting method, that confirms the protection of the initial state of the system. When the coupling between the system and the measuring apparatus increases beyond the range of validity of the analytical approximation, the initial state is still protected by the finite-time measurement, according with the exact numerical calculations.« less

Perturbation Experiments: Approaches for Metabolic Pathway Analysis in Bioreactors.

PubMed

Weiner, Michael; Tröndle, Julia; Albermann, Christoph; Sprenger, Georg A; Weuster-Botz, Dirk

2016-01-01

In the last decades, targeted metabolic engineering of microbial cells has become one of the major tools in bioprocess design and optimization. For successful application, a detailed knowledge is necessary about the relevant metabolic pathways and their regulation inside the cells. Since in vitro experiments cannot display process conditions and behavior properly, process data about the cells' metabolic state have to be collected in vivo. For this purpose, special techniques and methods are necessary. Therefore, most techniques enabling in vivo characterization of metabolic pathways rely on perturbation experiments, which can be divided into dynamic and steady-state approaches. To avoid any process disturbance, approaches which enable perturbation of cell metabolism in parallel to the continuing production process are reasonable. Furthermore, the fast dynamics of microbial production processes amplifies the need of parallelized data generation. These points motivate the development of a parallelized approach for multiple metabolic perturbation experiments outside the operating production reactor. An appropriate approach for in vivo characterization of metabolic pathways is presented and applied exemplarily to a microbial L-phenylalanine production process on a 15 L-scale.

Modeling of Mercury tides for recovery of gravity field and interior properties

NASA Astrophysics Data System (ADS)

Padovan, S.; Margot, J.; Hauck, S. A.; Lemoine, F. G.; Mazarico, E.; Peale, S. J.; Solomon, S. C.

2011-12-01

The radio science experiment on the MESSENGER mission allows the determination of the gravitational field of Mercury. In order to secure the best possible gravity-field recovery, it is important to model all the forces acting on the spacecraft. Here we study the perturbations induced on the spacecraft by the tides raised on Mercury by the Sun. The manner by which the tides affect the orbit of MESSENGER depends on the response of the planet to the tide-raising potential. This response is directly connected to the interior properties of Mercury, and its study can help improve our understanding of the physical and chemical properties of the planet. The standard approach of modeling the strongest tidal effect on the gravitational field is by introducing a time-varying component in the degree-two harmonic coefficients of the gravity field. The amplitude of these variations depends on known quantities (mass of the Sun and Mercury, radius of Mercury and its position and relative orientation with respect to the Sun) and on the Love number k2. The value of this parameter is sensitive (among other things) to the state of the core and to the rigidity of the mantle (which in turn depends on its chemical composition). An accurate value of k2 determined from orbit perturbations can be compared to values obtained with forward modeling of the interior of Mercury. The orbital geometry and physical environment of MESSENGER make the identification of the tidal perturbation difficult. Nevertheless, recent work has shown that in the case of Mars, careful study of the effect of tides on the spacecraft trajectory can help identify which orbital and observational geometries exhibit stronger tidal signatures and are apt to provide the best possible determination of k2. Our long-term goal is to evaluate k2 for a suite of interior models and to evaluate the sensitivity of k2 to key interior properties. We will describe the orbital geometry and the tidal perturbations acting on the spacecraft trajectory with both numerical and analytical approaches, and we will report on the status of the interior modeling efforts.

Impact of basic angle variations on the parallax zero point for a scanning astrometric satellite

NASA Astrophysics Data System (ADS)

Butkevich, Alexey G.; Klioner, Sergei A.; Lindegren, Lennart; Hobbs, David; van Leeuwen, Floor

2017-07-01

Context. Determination of absolute parallaxes by means of a scanning astrometric satellite such as Hipparcos or Gaia relies on the short-term stability of the so-called basic angle between the two viewing directions. Uncalibrated variations of the basic angle may produce systematic errors in the computed parallaxes. Aims: We examine the coupling between a global parallax shift and specific variations of the basic angle, namely those related to the satellite attitude with respect to the Sun. Methods: The changes in observables produced by small perturbations of the basic angle, attitude, and parallaxes were calculated analytically. We then looked for a combination of perturbations that had no net effect on the observables. Results: In the approximation of infinitely small fields of view, it is shown that certain perturbations of the basic angle are observationally indistinguishable from a global shift of the parallaxes. If these kinds of perturbations exist, they cannot be calibrated from the astrometric observations but will produce a global parallax bias. Numerical simulations of the astrometric solution, using both direct and iterative methods, confirm this theoretical result. For a given amplitude of the basic angle perturbation, the parallax bias is smaller for a larger basic angle and a larger solar aspect angle. In both these respects Gaia has a more favourable geometry than Hipparcos. In the case of Gaia, internal metrology is used to monitor basic angle variations. Additionally, Gaia has the advantage of detecting numerous quasars, which can be used to verify the parallax zero point.

A Singular Perturbation Approach for Time-Domain Assessment of Phase Margin

NASA Technical Reports Server (NTRS)

Zhu, J. Jim; Yang, Xiaojing; Hodel, A Scottedward

2010-01-01

This paper considers the problem of time-domain assessment of the Phase Margin (PM) of a Single Input Single Output (SISO) Linear Time-Invariant (LTI) system using a singular perturbation approach, where a SISO LTI fast loop system, whose phase lag increases monotonically with frequency, is introduced into the loop as a singular perturbation with a singular perturbation (time-scale separation) parameter Epsilon. First, a bijective relationship between the Singular Perturbation Margin (SPM) max and the PM of the nominal (slow) system is established with an approximation error on the order of Epsilon(exp 2). In proving this result, relationships between the singular perturbation parameter Epsilon, PM of the perturbed system, PM and SPM of the nominal system, and the (monotonically increasing) phase of the fast system are also revealed. These results make it possible to assess the PM of the nominal system in the time-domain for SISO LTI systems using the SPM with a standardized testing system called "PM-gauge," as demonstrated by examples. PM is a widely used stability margin for LTI control system design and certification. Unfortunately, it is not applicable to Linear Time-Varying (LTV) and Nonlinear Time-Varying (NLTV) systems. The approach developed here can be used to establish a theoretical as well as practical metric of stability margin for LTV and NLTV systems using a standardized SPM that is backward compatible with PM.

Simple analysis of scattering data with the Ornstein-Zernike equation

NASA Astrophysics Data System (ADS)

Kats, E. I.; Muratov, A. R.

2018-01-01

In this paper we propose and explore a method of analysis of the scattering experimental data for uniform liquidlike systems. In our pragmatic approach we are not trying to introduce by hands an artificial small parameter to work out a perturbation theory with respect to the known results, e.g., for hard spheres or sticky hard spheres (all the more that in the agreement with the notorious Landau statement, there is no physical small parameter for liquids). Instead of it being guided by the experimental data we are solving the Ornstein-Zernike equation with a trial (variational) form of the interparticle interaction potential. To find all needed correlation functions this variational input is iterated numerically to satisfy the Ornstein-Zernike equation supplemented by a closure relation. Our method is developed for spherically symmetric scattering objects, and our numeric code is written for such a case. However, it can be extended (at the expense of more involved computations and a larger amount of required experimental input information) for nonspherical particles. What is important for our approach is that it is sufficient to know experimental data in a relatively narrow range of the scattering wave vectors (q ) to compute the static structure factor in a much broader range of q . We illustrate by a few model and real experimental examples of the x-ray and neutron scattering data how the approach works.

A Systems Medicine Approach: Translating Emerging Science into Individualized Wellness

PubMed Central

Bland, J. S.

2017-01-01

In today's aging society, more people are living with lifestyle-related noncommunicable diseases (NCDs) such as cardiovascular disease, type 2 diabetes, obesity, and cancer. Numerous opinion-leader organizations recommend lifestyle medicine as the first-line approach in NCD prevention and treatment. However, there is a strong need for a personalized approach as “one-size-fits-all” public health recommendations have been insufficient in addressing the interindividual differences in the diverse populations. Advancement in systems biology and the “omics” technologies has allowed comprehensive analysis of how complex biological systems are impacted upon external perturbations (e.g., nutrition and exercise), and therefore is gradually pushing personalized lifestyle medicine toward reality. Clinicians and healthcare practitioners have a unique opportunity in advocating lifestyle medicine because patients see them as a reliable source of advice. However, there are still numerous technical and logistic challenges to overcome before personal “big data” can be translated into actionable and clinically relevant solutions. Clinicians are also facing various issues prior to bringing personalized lifestyle medicine to their practice. Nevertheless, emerging ground-breaking research projects have given us a glimpse of how systems thinking and computational methods may lead to personalized health advice. It is important that all stakeholders work together to create the needed paradigm shift in healthcare before the rising epidemic of NCDs overwhelm the society, the economy, and the dated health system. PMID:29164177

On physical and numerical instabilities arising in simulations of non-stationary radiatively cooling shocks

NASA Astrophysics Data System (ADS)

Badjin, D. A.; Glazyrin, S. I.; Manukovskiy, K. V.; Blinnikov, S. I.

2016-06-01

We describe our modelling of the radiatively cooling shocks and their thin shells with various numerical tools in different physical and calculational setups. We inspect structure of the dense shell, its formation and evolution, pointing out physical and numerical factors that sustain its shape and also may lead to instabilities. We have found that under certain physical conditions, the circular shaped shells show a strong bending instability and successive fragmentation on Cartesian grids soon after their formation, while remain almost unperturbed when simulated on polar meshes. We explain this by physical Rayleigh-Taylor-like instabilities triggered by corrugation of the dense shell surfaces by numerical noise. Conditions for these instabilities follow from both the shell structure itself and from episodes of transient acceleration during re-establishing of dynamical pressure balance after sudden radiative cooling onset. They are also easily excited by physical perturbations of the ambient medium. The widely mentioned non-linear thin shell instability, in contrast, in tests with physical perturbations is shown to have only limited chances to develop in real radiative shocks, as it seems to require a special spatial arrangement of fluctuations to be excited efficiently. The described phenomena also set new requirements on further simulations of the radiatively cooling shocks in order to be physically correct and free of numerical artefacts.

Modularity and the spread of perturbations in complex dynamical systems

NASA Astrophysics Data System (ADS)

Kolchinsky, Artemy; Gates, Alexander J.; Rocha, Luis M.

2015-12-01

We propose a method to decompose dynamical systems based on the idea that modules constrain the spread of perturbations. We find partitions of system variables that maximize "perturbation modularity," defined as the autocovariance of coarse-grained perturbed trajectories. The measure effectively separates the fast intramodular from the slow intermodular dynamics of perturbation spreading (in this respect, it is a generalization of the "Markov stability" method of network community detection). Our approach captures variation of modular organization across different system states, time scales, and in response to different kinds of perturbations: aspects of modularity which are all relevant to real-world dynamical systems. It offers a principled alternative to detecting communities in networks of statistical dependencies between system variables (e.g., "relevance networks" or "functional networks"). Using coupled logistic maps, we demonstrate that the method uncovers hierarchical modular organization planted in a system's coupling matrix. Additionally, in homogeneously coupled map lattices, it identifies the presence of self-organized modularity that depends on the initial state, dynamical parameters, and type of perturbations. Our approach offers a powerful tool for exploring the modular organization of complex dynamical systems.

Modularity and the spread of perturbations in complex dynamical systems.

PubMed

Kolchinsky, Artemy; Gates, Alexander J; Rocha, Luis M

2015-12-01

We propose a method to decompose dynamical systems based on the idea that modules constrain the spread of perturbations. We find partitions of system variables that maximize "perturbation modularity," defined as the autocovariance of coarse-grained perturbed trajectories. The measure effectively separates the fast intramodular from the slow intermodular dynamics of perturbation spreading (in this respect, it is a generalization of the "Markov stability" method of network community detection). Our approach captures variation of modular organization across different system states, time scales, and in response to different kinds of perturbations: aspects of modularity which are all relevant to real-world dynamical systems. It offers a principled alternative to detecting communities in networks of statistical dependencies between system variables (e.g., "relevance networks" or "functional networks"). Using coupled logistic maps, we demonstrate that the method uncovers hierarchical modular organization planted in a system's coupling matrix. Additionally, in homogeneously coupled map lattices, it identifies the presence of self-organized modularity that depends on the initial state, dynamical parameters, and type of perturbations. Our approach offers a powerful tool for exploring the modular organization of complex dynamical systems.

Representing perturbed dynamics in biological network models

NASA Astrophysics Data System (ADS)

Stoll, Gautier; Rougemont, Jacques; Naef, Felix

2007-07-01

We study the dynamics of gene activities in relatively small size biological networks (up to a few tens of nodes), e.g., the activities of cell-cycle proteins during the mitotic cell-cycle progression. Using the framework of deterministic discrete dynamical models, we characterize the dynamical modifications in response to structural perturbations in the network connectivities. In particular, we focus on how perturbations affect the set of fixed points and sizes of the basins of attraction. Our approach uses two analytical measures: the basin entropy H and the perturbation size Δ , a quantity that reflects the distance between the set of fixed points of the perturbed network and that of the unperturbed network. Applying our approach to the yeast-cell-cycle network introduced by Li [Proc. Natl. Acad. Sci. U.S.A. 101, 4781 (2004)] provides a low-dimensional and informative fingerprint of network behavior under large classes of perturbations. We identify interactions that are crucial for proper network function, and also pinpoint functionally redundant network connections. Selected perturbations exemplify the breadth of dynamical responses in this cell-cycle model.

Rapid distortion analysis of high speed homogeneous turbulence subject to periodic shear

DOE PAGES

Bertsch, Rebecca L.; Girimaji, Sharath S.

2015-12-30

The effect of unsteady shear forcing on small perturbation growth in compressible flow is investigated. In particular, flow-thermodynamic field interaction and the resulting effect on the phase-lag between applied shear and Reynolds stress are examined. Simplified linear analysis of the perturbation pressure equation reveals crucial differences between steady and unsteady shear effects. The analytical findings are validated with numerical simulations of inviscid rapid distortion theory (RDT) equations. In contrast to steadily sheared compressible flows, perturbations in the unsteady (periodic) forcing case do not experience an asymptotic growth phase. Further, the resonance growth phenomenon found in incompressible unsteady shear turbulence ismore » absent in the compressible case. Overall, the stabilizing influence of both unsteadiness and compressibility is compounded leading to suppression of all small perturbations. As a result, the underlying mechanisms are explained.« less

Three-dimensional simulation of the free shear layer using the vortex-in-cell method

NASA Technical Reports Server (NTRS)

Couet, B.; Buneman, O.; Leonard, A.

1979-01-01

We present numerical simulations of the evolution of a mixing layer from an initial state of uniform vorticity with simple two- and three-dimensional small perturbations. A new method for tracing a large number of three-dimensional vortex filaments is used in the simulations. Vortex tracing by Biot-Savart interaction originally implied ideal (non-viscous) flow, but we use a 3-d mesh, Fourier transforms and filtering for vortex tracing, which implies 'modeling' of subgrid scale motion and hence some viscosity. Streamwise perturbations lead to the usual roll-up of vortex patterns with spanwise uniformity maintained. Remarkably, spanwise perturbations generate streamwise distortions of the vortex filaments and the combination of both perturbations leads to patterns with interesting features discernable in the movies and in the records of enstrophy and energy for the three components of the flow.

Qubit dephasing due to low-frequency noise.

NASA Astrophysics Data System (ADS)

Sverdlov, Victor; Rabenstein, Kristian; Averin, Dmitri

2004-03-01

We have numerically investigated the effects of the classical low-frequency noise on the qubit dynamics beyond the standard lowest-order perturbation theory in coupling. Noise is generated as a random process with a correlation function characterized by two parameters, the amplitude v0 and the cut-off frequency 2π/τ. Time evolution of the density matrix was averaged over up to 10^7 noise realizations. Contrary to the relaxation time T_1, which for v_0<ω, where ω is the qubit oscillation frequency, is always given correctly by the ``golden-rule'' expression, the dephasing time deviates from the perturbation-theory result, when (v_0/ω)^2(ωτ) ≥1. In this regime, even for unbiased qubit for which the pure dephasing vanishes in perturbation theory, the dephasing is much larger than it's perturbation-theory value 1/(2 T_1).

Rapid distortion analysis of high speed homogeneous turbulence subject to periodic shear

DOE Office of Scientific and Technical Information (OSTI.GOV)

Bertsch, Rebecca L., E-mail: rlb@lanl.gov; Girimaji, Sharath S., E-mail: girimaji@aero.tamu.edu

2015-12-15

The effect of unsteady shear forcing on small perturbation growth in compressible flow is investigated. In particular, flow-thermodynamic field interaction and the resulting effect on the phase-lag between applied shear and Reynolds stress are examined. Simplified linear analysis of the perturbation pressure equation reveals crucial differences between steady and unsteady shear effects. The analytical findings are validated with numerical simulations of inviscid rapid distortion theory (RDT) equations. In contrast to steadily sheared compressible flows, perturbations in the unsteady (periodic) forcing case do not experience an asymptotic growth phase. Further, the resonance growth phenomenon found in incompressible unsteady shear turbulence ismore » absent in the compressible case. Overall, the stabilizing influence of both unsteadiness and compressibility is compounded leading to suppression of all small perturbations. The underlying mechanisms are explained.« less

Many-body formalism for fermions: The partition function

NASA Astrophysics Data System (ADS)

Watson, D. K.

2017-09-01

The partition function, a fundamental tenet in statistical thermodynamics, contains in principle all thermodynamic information about a system. It encapsulates both microscopic information through the quantum energy levels and statistical information from the partitioning of the particles among the available energy levels. For identical particles, this statistical accounting is complicated by the symmetry requirements of the allowed quantum states. In particular, for Fermi systems, the enforcement of the Pauli principle is typically a numerically demanding task, responsible for much of the cost of the calculations. The interplay of these three elements—the structure of the many-body spectrum, the statistical partitioning of the N particles among the available levels, and the enforcement of the Pauli principle—drives the behavior of mesoscopic and macroscopic Fermi systems. In this paper, we develop an approach for the determination of the partition function, a numerically difficult task, for systems of strongly interacting identical fermions and apply it to a model system of harmonically confined, harmonically interacting fermions. This approach uses a recently introduced many-body method that is an extension of the symmetry-invariant perturbation method (SPT) originally developed for bosons. It uses group theory and graphical techniques to avoid the heavy computational demands of conventional many-body methods which typically scale exponentially with the number of particles. The SPT application of the Pauli principle is trivial to implement since it is done "on paper" by imposing restrictions on the normal-mode quantum numbers at first order in the perturbation. The method is applied through first order and represents an extension of the SPT method to excited states. Our method of determining the partition function and various thermodynamic quantities is accurate and efficient and has the potential to yield interesting insight into the role played by the Pauli principle and the influence of large degeneracies on the emergence of the thermodynamic behavior of large-N systems.

Possibilities of Bragg filtering structures based on subwavelength grating guiding mechanism (Conference Presentation)

NASA Astrophysics Data System (ADS)

Kwiecien, Pavel; Litvik, Ján.; Richter, Ivan; Ctyroký, Jirí; Cheben, Pavel

2017-05-01

Silicon-on-insulator (SOI), as the most promising platform, for advanced photonic integrated structures, employs a high refractive index contrast between the silicon "core" and surrounding media. One of the recent new ideas within this field is based on the alternative formation of the subwavelength sized (quasi)periodic structures, manifesting as an effective medium with respect to propagating light. Such structures relay on Bloch wave propagation concept, in contrast to standard index guiding mechanism. Soon after the invention of such subwavelength grating (SWG) waveguides, the scientists concentrated on various functional elements such as couplers, crossings, mode transformers, convertors, MMI couplers, polarization converters, resonators, Bragg filters, and others. Our contribution is devoted to a detailed numerical analysis and design considerations of Bragg filtering structures based on SWG idea. Based on our previous studies where we have shown impossibility of application of various 2 and "2.5" dimensional methods for the proper numerical analysis, here we effectively use two independent but similar in-house approaches based on 3D Fourier modal methods, namely aperiodic rigorous coupled wave analysis (aRCWA) and bidirectional expansion and propagation method based on Fourier series (BEX) tools. As it was recently demonstrated, SWG Bragg filters are feasible. Based on this idea, we propose, simulate, and optimize spectral characteristics of such filters. In particular, we have investigated several possibilities of modifications of original SWG waveguides towards the Bragg filtering, including firstly - simple single-segment changes in position, thickness, and width, and secondly - several types of Si inclusions, in terms of perturbed width and thickness (and their combinations). The leading idea was to obtain required (e.g. sufficiently narrow) spectral characteristic while keeping the minimum size of Si features large enough. We have found that the second approach with the single element perturbations can provide promising designs. Furthermore, even more complex filtering SWG structures can be considered.

Dissipative exciton transfer in donor-bridge-acceptor systems: numerical renormalization group calculation of equilibrium properties.

PubMed

Tornow, Sabine; Tong, Ning-Hua; Bulla, Ralf

2006-07-05

We present a detailed model study of exciton transfer processes in donor-bridge-acceptor (DBA) systems. Using a model which includes the intermolecular Coulomb interaction and the coupling to a dissipative environment we calculate the phase diagram, the absorption spectrum as well as dynamic equilibrium properties with the numerical renormalization group. This method is non-perturbative and therefore allows one to cover the full parameter space, especially the case when the intermolecular Coulomb interaction is of the same order as the coupling to the environment and perturbation theory cannot be applied. For DBA systems with up to six sites we found a transition to the localized phase (self-trapping) depending on the coupling to the dissipative environment. We discuss various criteria which favour delocalized exciton transfer.

Numerically Integrated Orbits of the Major Saturnian Satellites fit to Earthbased Observations

NASA Technical Reports Server (NTRS)

Jacobson, R. A.; Vaughan, R. M.

1993-01-01

We have fit numerically integrated orbits of the eight major satellites of Saturn to all available astrometric and meridian circle observations for the period of 1971 to 1992. The integration was carried out in cartesian coordinates in the J2000 system. The force model included the gravitational effects of the oblate primary, the mutual perturbations of the satellites, and perturbations due to Jupiter and the Sun. Values of the gravitational parameters of the Saturnian system, e.g. planet and satellite masses, were taken from Campbell, et. al., 1989, only the epoch state vectors of the satellites were adjusted to obtain orbits which fit the observations. All astrometric data was processed in the form of satellite relative positions which were weighted according to observer and opposition to reflect the varying data quality...

Singular perturbation techniques for real time aircraft trajectory optimization and control

NASA Technical Reports Server (NTRS)

Calise, A. J.; Moerder, D. D.

1982-01-01

The usefulness of singular perturbation methods for developing real time computer algorithms to control and optimize aircraft flight trajectories is examined. A minimum time intercept problem using F-8 aerodynamic and propulsion data is used as a baseline. This provides a framework within which issues relating to problem formulation, solution methodology and real time implementation are examined. Theoretical questions relating to separability of dynamics are addressed. With respect to implementation, situations leading to numerical singularities are identified, and procedures for dealing with them are outlined. Also, particular attention is given to identifying quantities that can be precomputed and stored, thus greatly reducing the on-board computational load. Numerical results are given to illustrate the minimum time algorithm, and the resulting flight paths. An estimate is given for execution time and storage requirements.

A New Method for Computing Three-Dimensional Capture Fraction in Heterogeneous Regional Systems using the MODFLOW Adjoint Code

NASA Astrophysics Data System (ADS)

Clemo, T. M.; Ramarao, B.; Kelly, V. A.; Lavenue, M.

2011-12-01

Capture is a measure of the impact of groundwater pumping upon groundwater and surface water systems. The computation of capture through analytical or numerical methods has been the subject of articles in the literature for several decades (Bredehoeft et al., 1982). Most recently Leake et al. (2010) described a systematic way to produce capture maps in three-dimensional systems using a numerical perturbation approach in which capture from streams was computed using unit rate pumping at many locations within a MODFLOW model. The Leake et al. (2010) method advances the current state of computing capture. A limitation stems from the computational demand required by the perturbation approach wherein days or weeks of computational time might be required to obtain a robust measure of capture. In this paper, we present an efficient method to compute capture in three-dimensional systems based upon adjoint states. The efficiency of the adjoint method will enable uncertainty analysis to be conducted on capture calculations. The USGS and INTERA have collaborated to extend the MODFLOW Adjoint code (Clemo, 2007) to include stream-aquifer interaction and have applied it to one of the examples used in Leake et al. (2010), the San Pedro Basin MODFLOW model. With five layers and 140,800 grid blocks per layer, the San Pedro Basin model, provided an ideal example data set to compare the capture computed from the perturbation and the adjoint methods. The capture fraction map produced from the perturbation method for the San Pedro Basin model required significant computational time to compute and therefore the locations for the pumping wells were limited to 1530 locations in layer 4. The 1530 direct simulations of capture require approximately 76 CPU hours. Had capture been simulated in each grid block in each layer, as is done in the adjoint method, the CPU time would have been on the order of 4 years. The MODFLOW-Adjoint produced the capture fraction map of the San Pedro Basin model at 704,000 grid blocks (140,800 grid blocks x 5 layers) in just 6 minutes. The capture fraction maps from the perturbation and adjoint methods agree closely. The results of this study indicate that the adjoint capture method and its associated computational efficiency will enable scientists and engineers facing water resource management decisions to evaluate the sensitivity and uncertainty of impacts to regional water resource systems as part of groundwater supply strategies. Bredehoeft, J.D., S.S. Papadopulos, and H.H. Cooper Jr, Groundwater: The water budget myth. In Scientific Basis of Water-Resources Management, ed. National Research Council (U.S.), Geophysical Study Committee, 51-57. Washington D.C.: National Academy Press, 1982. Clemo, Tom, MODFLOW-2005 Ground-Water Model-Users Guide to Adjoint State based Sensitivity Process (ADJ), BSU CGISS 07-01, Center for the Geophysical Investigation of the Shallow Subsurface, Boise State University, 2007. Leake, S.A., H.W. Reeves, and J.E. Dickinson, A New Capture Fraction Method to Map How Pumpage Affects Surface Water Flow, Ground Water, 48(5), 670-700, 2010.

Analysis and computation of a least-squares method for consistent mesh tying

DOE PAGES

Day, David; Bochev, Pavel

2007-07-10

We report in the finite element method, a standard approach to mesh tying is to apply Lagrange multipliers. If the interface is curved, however, discretization generally leads to adjoining surfaces that do not coincide spatially. Straightforward Lagrange multiplier methods lead to discrete formulations failing a first-order patch test [T.A. Laursen, M.W. Heinstein, Consistent mesh-tying methods for topologically distinct discretized surfaces in non-linear solid mechanics, Internat. J. Numer. Methods Eng. 57 (2003) 1197–1242]. This paper presents a theoretical and computational study of a least-squares method for mesh tying [P. Bochev, D.M. Day, A least-squares method for consistent mesh tying, Internat. J.more » Numer. Anal. Modeling 4 (2007) 342–352], applied to the partial differential equation -∇ 2φ+αφ=f. We prove optimal convergence rates for domains represented as overlapping subdomains and show that the least-squares method passes a patch test of the order of the finite element space by construction. To apply the method to subdomain configurations with gaps and overlaps we use interface perturbations to eliminate the gaps. Finally, theoretical error estimates are illustrated by numerical experiments.« less

Analytical method for the effects of the asteroid belt on planetary orbits

NASA Technical Reports Server (NTRS)

Mayo, A. P.

1979-01-01

Analytic expressions are derived for the perturbation of planetary orbits due to a thick constant-density asteroid belt. The derivations include extensions and adaptations of Plakhov's (1968) analytic expressions for the perturbations in five of the orbital elements for closed orbits around Saturn's rings. The equations of Plakhov are modified to include the effect of ring thickness, and additional equations are derived for the perturbations in the sixth orbital element, the mean anomaly. The gravitational potential and orbital perturbations are derived for the asteroid belt with and without thickness, and for a hoop approximation to the belt. The procedures are also applicable to Saturn's rings and the newly discovered rings of Uranus. The effects of the asteroid belt thickness on the gravitational potential coefficients and the orbital motions are demonstrated. Comparisons between the Mars orbital perturbations obtained by using the analytic expressions and those obtained by numerical integration are discussed. The effects of the asteroid belt on earth-based ranging to Mars are also demonstrated.

Optimal Linear Responses for Markov Chains and Stochastically Perturbed Dynamical Systems

NASA Astrophysics Data System (ADS)

Antown, Fadi; Dragičević, Davor; Froyland, Gary

2018-03-01

The linear response of a dynamical system refers to changes to properties of the system when small external perturbations are applied. We consider the little-studied question of selecting an optimal perturbation so as to (i) maximise the linear response of the equilibrium distribution of the system, (ii) maximise the linear response of the expectation of a specified observable, and (iii) maximise the linear response of the rate of convergence of the system to the equilibrium distribution. We also consider the inhomogeneous, sequential, or time-dependent situation where the governing dynamics is not stationary and one wishes to select a sequence of small perturbations so as to maximise the overall linear response at some terminal time. We develop the theory for finite-state Markov chains, provide explicit solutions for some illustrative examples, and numerically apply our theory to stochastically perturbed dynamical systems, where the Markov chain is replaced by a matrix representation of an approximate annealed transfer operator for the random dynamical system.

Effects of the underlying topology on perturbation spreading in the Axelrod model for cultural dissemination

NASA Astrophysics Data System (ADS)

Kim, Yup; Cho, Minsoo; Yook, Soon-Hyung

2011-10-01

We study the effects of the underlying topologies on a single feature perturbation imposed to the Axelrod model of consensus formation. From the numerical simulations we show that there are successive updates which are similar to avalanches in many self-organized criticality systems when a perturbation is imposed. We find that the distribution of avalanche size satisfies the finite-size scaling (FSS) ansatz on two-dimensional lattices and random networks. However, on scale-free networks with the degree exponent γ≤3 we show that the avalanche size distribution does not satisfy the FSS ansatz. The results indicate that the disordered configurations on two-dimensional lattices or on random networks are still stable against the perturbation in the limit N (network size) →∞. However, on scale-free networks with γ≤3 the perturbation always drives the disordered phase into an ordered phase. The possible relationship between the properties of phase transition of the Axelrod model and the avalanche distribution is also discussed.

On the Stability of Shocks with Particle Pressure

NASA Astrophysics Data System (ADS)

Finazzi, Stefano; Vietri, Mario

2008-11-01

We perform a linear stability analysis for corrugations of a Newtonian shock, with particle pressure included, for an arbitrary diffusion coefficient. We study first the dispersion relation for homogeneous media, showing that, besides the conventional pressure waves and entropy/vorticity disturbances, two new perturbation modes exist, dominated by the particles' pressure and damped by diffusion. We show that, due to particle diffusion into the upstream region, the fluid will be perturbed also upstream; we treat these perturbation in the short-wavelength (WKBJ) regime. We then show how to construct a corrugational mode for the shock itself, one, that is, where the shock executes free oscillations (possibly damped or growing) and sheds perturbations away from itself; this global mode requires the new modes. Then, using the perturbed Rankine-Hugoniot conditions, we show that this leads to the determination of the corrugational eigenfrequency. We solve numerically the equations for the eigenfrequency in the WKBJ regime for the models of Amato & Blasi, showing that they are stable. We then discuss the differences between our treatment and previous work.

Interval stability for complex systems

NASA Astrophysics Data System (ADS)

Klinshov, Vladimir V.; Kirillov, Sergey; Kurths, Jürgen; Nekorkin, Vladimir I.

2018-04-01

Stability of dynamical systems against strong perturbations is an important problem of nonlinear dynamics relevant to many applications in various areas. Here, we develop a novel concept of interval stability, referring to the behavior of the perturbed system during a finite time interval. Based on this concept, we suggest new measures of stability, namely interval basin stability (IBS) and interval stability threshold (IST). IBS characterizes the likelihood that the perturbed system returns to the stable regime (attractor) in a given time. IST provides the minimal magnitude of the perturbation capable to disrupt the stable regime for a given interval of time. The suggested measures provide important information about the system susceptibility to external perturbations which may be useful for practical applications. Moreover, from a theoretical viewpoint the interval stability measures are shown to bridge the gap between linear and asymptotic stability. We also suggest numerical algorithms for quantification of the interval stability characteristics and demonstrate their potential for several dynamical systems of various nature, such as power grids and neural networks.

Control of spiral waves and turbulent states in a cardiac model by travelling-wave perturbations

NASA Astrophysics Data System (ADS)

Wang, Peng-Ye; Xie, Ping; Yin, Hua-Wei

2003-06-01

We propose a travelling-wave perturbation method to control the spatiotemporal dynamics in a cardiac model. It is numerically demonstrated that the method can successfully suppress the wave instability (alternans in action potential duration) in the one-dimensional case and convert spiral waves and turbulent states to the normal travelling wave states in the two-dimensional case. An experimental scheme is suggested which may provide a new design for a cardiac defibrillator.

Computation of diverging sums based on a finite number of terms

NASA Astrophysics Data System (ADS)

Lv, Q. Z.; Norris, S.; Pelphrey, R.; Su, Q.; Grobe, R.

2017-10-01

We propose a numerical method that permits us to compute the sum of a diverging series from only the first N terms by generalizing the traditional Borel technique. The method is rather robust and can be used to recover the ground state energy from the diverging perturbation theory for quantum field theoretical systems that are spatially constrained. Surprisingly, even the corresponding eigenvectors can be generated despite the intrinsic non-perturbative nature of bound state problems.

A comparison of methods to estimate seismic phase delays--Numerical examples for coda wave interferometry

USGS Publications Warehouse

Mikesell, T. Dylan; Malcolm, Alison E.; Yang, Di; Haney, Matthew M.

2015-01-01

Time-shift estimation between arrivals in two seismic traces before and after a velocity perturbation is a crucial step in many seismic methods. The accuracy of the estimated velocity perturbation location and amplitude depend on this time shift. Windowed cross correlation and trace stretching are two techniques commonly used to estimate local time shifts in seismic signals. In the work presented here, we implement Dynamic Time Warping (DTW) to estimate the warping function – a vector of local time shifts that globally minimizes the misfit between two seismic traces. We illustrate the differences of all three methods compared to one another using acoustic numerical experiments. We show that DTW is comparable to or better than the other two methods when the velocity perturbation is homogeneous and the signal-to-noise ratio is high. When the signal-to-noise ratio is low, we find that DTW and windowed cross correlation are more accurate than the stretching method. Finally, we show that the DTW algorithm has better time resolution when identifying small differences in the seismic traces for a model with an isolated velocity perturbation. These results impact current methods that utilize not only time shifts between (multiply) scattered waves, but also amplitude and decoherence measurements. DTW is a new tool that may find new applications in seismology and other geophysical methods (e.g., as a waveform inversion misfit function).

A collision history-based approach to Sensitivity/Perturbation calculations in the continuous energy Monte Carlo code SERPENT

DOE Office of Scientific and Technical Information (OSTI.GOV)

Giuseppe Palmiotti

In this work, the implementation of a collision history-based approach to sensitivity/perturbation calculations in the Monte Carlo code SERPENT is discussed. The proposed methods allow the calculation of the eects of nuclear data perturbation on several response functions: the eective multiplication factor, reaction rate ratios and bilinear ratios (e.g., eective kinetics parameters). SERPENT results are compared to ERANOS and TSUNAMI Generalized Perturbation Theory calculations for two fast metallic systems and for a PWR pin-cell benchmark. New methods for the calculation of sensitivities to angular scattering distributions are also presented, which adopts fully continuous (in energy and angle) Monte Carlo estimators.

Feynman-like rules for calculating n-point correlators of the primordial curvature perturbation

NASA Astrophysics Data System (ADS)

Valenzuela-Toledo, César A.; Rodríguez, Yeinzon; Beltrán Almeida, Juan P.

2011-10-01

A diagrammatic approach to calculate n-point correlators of the primordial curvature perturbation ζ was developed a few years ago following the spirit of the Feynman rules in Quantum Field Theory. The methodology is very useful and time-saving, as it is for the case of the Feynman rules in the particle physics context, but, unfortunately, is not very well known by the cosmology community. In the present work, we extend such an approach in order to include not only scalar field perturbations as the generators of ζ, but also vector field perturbations. The purpose is twofold: first, we would like the diagrammatic approach (which we would call the Feynman-like rules) to become widespread among the cosmology community; second, we intend to give an easy tool to formulate any correlator of ζ for those cases that involve vector field perturbations and that, therefore, may generate prolonged stages of anisotropic expansion and/or important levels of statistical anisotropy. Indeed, the usual way of formulating such correlators, using the Wick's theorem, may become very clutter and time-consuming.

Primordial power spectra for scalar perturbations in loop quantum cosmology

DOE Office of Scientific and Technical Information (OSTI.GOV)

De Blas, Daniel Martín; Olmedo, Javier, E-mail: d.martindeblas@uandresbello.edu, E-mail: jolmedo@lsu.edu

We provide the power spectrum of small scalar perturbations propagating in an inflationary scenario within loop quantum cosmology. We consider the hybrid quantization approach applied to a Friedmann-Robertson-Walker spacetime with flat spatial sections coupled to a massive scalar field. We study the quantum dynamics of scalar perturbations on an effective background within this hybrid approach. We consider in our study adiabatic states of different orders. For them, we find that the hybrid quantization is in good agreement with the predictions of the dressed metric approach. We also propose an initial vacuum state for the perturbations, and compute the primordial andmore » the anisotropy power spectrum in order to qualitatively compare with the current observations of Planck mission. We find that our vacuum state is in good agreement with them, showing a suppression of the power spectrum for large scale anisotropies. We compare with other choices already studied in the literature.« less

Self-organized behavior of modeled shoreline shapes

NASA Astrophysics Data System (ADS)

Ashton, A.; Murray, A. B.

2003-04-01

Whenever waves approach a coast and break at oblique angles, they drive a current along the shore. This current, along with wave-induced sediment suspension, transports relatively large amounts of sediment, affecting the shape and evolution of a coastline. Traditionally, researchers have assumed that alongshore sediment transport will diffuse, or smooth, bumps along a shoreline. Recent research, however, shows that when the angle between wave crests in deep water and the shoreline is sufficiently high (greater than approximately 45 degrees), a shoreline is unstable. Linear stability analysis does not predict that this instability will cause a preferred wavelength of shoreline perturbation growth or that organized patterns will emerge. However, a simple numerical model of shoreline change shows those when there is a predominance of high angle waves approaching a shoreline, finite-amplitude features will develop that interact with each other and increase in wavelength over time, translating in the direction of net alongshore sediment transport. Some of these simulated features resemble naturally occurring shoreline features, such as 'alongshore sandwaves', 'ords', 'cuspate spits', and 'cuspate forelands'. By varying two wave climate parameters, one describing the relative dominance of waves approaching at high angles and the other controlling the signs of the approach angle of incoming waves (i.e., the asymmetry of waves approaching from the right vs. the left), we investigate how the attributes of the input wave climate determine the aspect ratio and characteristic form of the simulated features. Varying these two parameters also affects the wavelength of the initially fastest growing perturbation. By tracking the average wavelength of simulated features, which increases over time for all simulations, we show that more complicated phenomena, such as rapid period doubling, can dominate simulated shoreline evolution. These rich behaviors result from large-scale emergent interactions. Although the wave distribution determines the character of shoreline features, their specific configuration and evolution is sensitively dependant on both initial conditions and the stochastic sequencing of wave approach angles.

Numerical modeling of a point-source image under relative motion of radiation receiver and atmosphere

NASA Astrophysics Data System (ADS)

Kucherov, A. N.; Makashev, N. K.; Ustinov, E. V.

1994-02-01

A procedure is proposed for numerical modeling of instantaneous and averaged (over various time intervals) distant-point-source images perturbed by a turbulent atmosphere that moves relative to the radiation receiver. Examples of image calculations under conditions of the significant effect of atmospheric turbulence in an approximation of geometrical optics are presented and analyzed.

A matrix-based approach to solving the inverse Frobenius-Perron problem using sequences of density functions of stochastically perturbed dynamical systems

NASA Astrophysics Data System (ADS)

Nie, Xiaokai; Coca, Daniel

2018-01-01

The paper introduces a matrix-based approach to estimate the unique one-dimensional discrete-time dynamical system that generated a given sequence of probability density functions whilst subjected to an additive stochastic perturbation with known density.

A matrix-based approach to solving the inverse Frobenius-Perron problem using sequences of density functions of stochastically perturbed dynamical systems.

PubMed

Nie, Xiaokai; Coca, Daniel

2018-01-01

The paper introduces a matrix-based approach to estimate the unique one-dimensional discrete-time dynamical system that generated a given sequence of probability density functions whilst subjected to an additive stochastic perturbation with known density.

Impact of seasonality upon the dynamics of a novel pathogen in a seabird colony

NASA Astrophysics Data System (ADS)

O'Regan, S. M.

2008-11-01

A seasonally perturbed variant of the basic Susceptible-Infected-Recovered (SIR) model in epidemiology is considered in this paper. The effect of seasonality on an IR system of ordinary differential equations describing the dynamics of a novel pathogen, e.g., highly pathogenic avian influenza, in a seabird colony is investigated. The method of Lyapunov functions is used to determine the long-term behaviour of this system. Numerical simulations of the seasonally perturbed IR system indicate that the system exhibits complex dynamics as the amplitude of the seasonal perturbation term is increased. These findings suggest that seasonality may exert a considerable effect on the dynamics of epidemics in a seabird colony.

Perturbations of the Kerr spacetime in horizon-penetrating coordinates

NASA Astrophysics Data System (ADS)

Campanelli, Manuela; Khanna, Gaurav; Laguna, Pablo; Pullin, Jorge; Ryan, Michael P.

2001-04-01

We derive the Teukolsky equation for perturbations of a Kerr spacetime when the spacetime metric is written in either ingoing or outgoing Kerr-Schild form. We also write explicit formulae for setting up the initial data for the Teukolsky equation in the time domain in terms of a 3-metric and an extrinsic curvature. The motivation of this work is to have in place a formalism to study the evolution in the `close limit' of two recently proposed solutions to the initial-value problem in general relativity that are based on Kerr-Schild slicings. A perturbative formalism in horizon-penetrating coordinates is also very desirable in connection with numerical relativity simulations using black hole `excision'.

Adiabatic perturbation theory for atoms and molecules in the low-frequency regime

NASA Astrophysics Data System (ADS)

Martiskainen, Hanna; Moiseyev, Nimrod

2017-12-01

There is an increasing interest in the photoinduced dynamics in the low frequency, ω, regime. The multiphoton absorptions by molecules in strong laser fields depend on the polarization of the laser and on the molecular structure. The unique properties of the interaction of atoms and molecules with lasers in the low-frequency regime imply new concepts and directions in strong-field light-matter interactions. Here we represent a perturbational approach for the calculations of the quasi-energy spectrum in the low-frequency regime, which avoids the construction of the Floquet operator with extremely large number of Floquet channels. The zero-order Hamiltonian in our perturbational approach is the adiabatic Hamiltonian where the atoms/molecules are exposed to a dc electric field rather than to ac-field. This is in the spirit of the first step in the Corkum three-step model. The second-order perturbation correction terms are obtained when i ℏ ω ∂/∂ τ serves as a perturbation and τ is a dimensionless variable. The second-order adiabatic perturbation scheme is found to be an excellent approach for calculating the ac-field Floquet solutions in our test case studies of a simple one-dimensional time-periodic model Hamiltonian. It is straightforward to implement the perturbation approach presented here for calculating atomic and molecular energy shifts (positions) due to the interaction with low-frequency ac-fields using high-level electronic structure methods. This is enabled since standard quantum chemistry packages allow the calculations of atomic and molecular energy shifts due to the interaction with dc-fields. In addition to the shift of the energy positions, the energy widths (inverse lifetimes) can be obtained at the same level of theory. These energy shifts are functions of the laser parameters (low frequency, intensity, and polarization).

Numerical simulation for solution of space-time fractional telegraphs equations with local fractional derivatives via HAFSTM

NASA Astrophysics Data System (ADS)

Pandey, Rishi Kumar; Mishra, Hradyesh Kumar

2017-11-01

In this paper, the semi-analytic numerical technique for the solution of time-space fractional telegraph equation is applied. This numerical technique is based on coupling of the homotopy analysis method and sumudu transform. It shows the clear advantage with mess methods like finite difference method and also with polynomial methods similar to perturbation and Adomian decomposition methods. It is easily transform the complex fractional order derivatives in simple time domain and interpret the results in same meaning.

Development of a Perfectly Matched Layer Technique for a Discontinuous-Galerkin Spectral-Element Method

NASA Technical Reports Server (NTRS)

Garai, Anirban; Murman, Scott M.; Madavan, Nateri K.

2016-01-01

The numerical simulation of many aerodynamic non-periodic flows of practical interest involves discretized computational domains that often must be artificially truncated. Appropriate boundary conditions are required at these truncated domain boundaries, and ideally, these boundary conditions should be perfectly "absorbing" or "nonreflecting" so that they do not contaminate the flow field in the interior of the domain. The proper specification of these boundaries is critical to the stability, accuracy, convergence, and quality of the numerical solution, and has been the topic of considerable research. The need for accurate boundary specification has been underscored in recent years with efforts to apply higher-fidelity methods (DNS, LES) in conjunction with high-order low-dissipation numerical schemes to realistic flow configurations. One of the most popular choices for specifying these boundaries is the characteristics-based boundary condition where the linearized flow field at the boundaries are decomposed into characteristic waves using either one-dimensional Riemann or other multi-dimensional Riemann approximations. The values of incoming characteristics are then suitably modified. The incoming characteristics are specified at the in flow boundaries, and at the out flow boundaries the variation of the incoming characteristic is zeroed out to ensure no reflection. This, however, makes the problem ill-posed requiring the use of an ad-hoc parameter to allow small reflections that make the solution stable. Generally speaking, such boundary conditions work reasonably well when the characteristic flow direction is normal to the boundary, but reflects spurious energy otherwise. An alternative to the characteristic-based boundary condition is to add additional "buffer" regions to the main computational domain near the artificial boundaries, and solve a different set of equations in the buffer region in order to minimize acoustic reflections. One approach that has been used involves modeling the pressure fluctuations as acoustic waves propagating in the far-field relative to a single noise-source inside the buffer region. This approach treats vorticity-induced pressure fluctuations the same as acoustic waves. Another popular approach, often referred to as the "sponge layer," attempts to dampen the flow perturbations by introducing artificial dissipation in the buffer region. Although the artificial dissipation removes all perturbations inside the sponge layer, incoming waves are still reflected from the interface boundary between the computational domain and the sponge layer. The effect of these refkections can be somewhat mitigated by appropriately selecting the artificial dissipation strength and the extent of the sponge layer. One of the most promising variants on the buffer region approach is the Perfectly Matched Layer (PML) technique. The PML technique mitigates spurious reflections from boundaries and interfaces by dampening the perturbation modes inside the buffer region such that their eigenfunctions remain unchanged. The technique was first developed by Berenger for application to problems involving electromagnetic wave propagation. It was later extended to the linearized Euler, Euler and Navier-Stokes equations by Hu and his coauthors. The PML technique ensures the no-reflection property for all waves, irrespective of incidence angle, wavelength, and propagation direction. Although the technique requires the solution of a set of auxiliary equations, the computational overhead is easily justified since it allows smaller domain sizes and can provide better accuracy, stability, and convergence of the numerical solution. In this paper, the PML technique is developed in the context of a high-order spectral-element Discontinuous Galerkin (DG) method. The technique is compared to other approaches to treating the in flow and out flow boundary, such as those based on using characteristic boundary conditions and sponge layers. The superiority of the current PML technique over other approaches is demonstrated for a range of test cases, viz., acoustic pulse propagation, convective vortex, shear layer flow, and low-pressure turbine cascade flow. The paper is structured as follows. We first derive the PML equations from the non{linear Euler equations. A short description of the higher-order DG method used is then described. Preliminary results for the four test cases considered are then presented and discussed. Details regarding current work that will be included in the final paper are also provided.

Large Black Holes in the Randall-Sundrum II Model

NASA Astrophysics Data System (ADS)

Yaghoobpour Tari, Shima

The Einstein equation with a negative cosmological constant ! in the five dimensions for the Randall-Sundrum II model, which includes a black hole, has been solved numerically. We have constructed an AdS5-CFT 4 solution numerically, using a spectral method to minimize the integral of the square of the error of the Einstein equation, with 210 parameters to be determined by optimization. This metric is conformal to the Schwarzschild metric at an AdS5 boundary with an infinite scale factor. So, we consider this solution as an infinite-mass black hole solution. We have rewritten the infinite-mass black hole in the Fefferman-Graham form and obtained the numerical components of the CFT energy-momentum tensor. Using them, we have perturbed the metric to relocate the brane from infinity and obtained a large static black hole solution for the Randall- Sundrum II model. The changes of mass, entropy, temperature and area of the large black hole from the Schwarzschild metric are studied up to the first order for the perturbation parameter 1/(-Λ5M 2). The Hawking temperature and entropy for our large black hole have the same values as the Schwarzschild metric with the same mass, but the horizon area is increased by about 4.7/(-Λ5). Figueras, Lucietti, and Wiseman found an AdS5-CFT4 solution using an independent and different method from us, called the Ricci-DeTurck-flow method. Then, Figueras and Wiseman perturbed this solution in a same way as we have done and obtained the solution for the large black hole in the Randall-Sundrum II model. These two numerical solutions are the first mathematical proofs for having a large black hole in the Randall-Sundrum II. We have compared their results with ours for the CFT energy-momentum tensor components and the perturbed metric. We have shown that the results are closely in agreement, which can be considered as evidence that the solution for the large black hole in the Randall-Sundrum II model exists.

An Overview of A Perturbation Analysis for Uni-directionally Coupled Vibratory Gyroscopes

NASA Astrophysics Data System (ADS)

Vu, Huy; Palacios, Antonio; In, Visarath; Longhini, Patrick; Neff, Joseph

2011-04-01

The complex behaviours of gyroscope systems have been scientifically researched and thoroughly studied for decades. Most of scientific research involving gyroscopes specifically concentrates on studying the designs and fabrications at the circuitry level. Although gaining a recent popularity with the low cost of MEMS device that offers an attractive approach for gyroscope fabrications, its performance is far from meeting the requirements for an inertial grade guidance system. To improve the performance, our current research is theoretically focusing upon investigating the dynamics of vibratory gyroscopes coupled in a ring configuration. Particularly, a certain topology of arrangements among coupled gyroscopes can be designed and studied to enhance robustness. The main operation depends mostly on an external source for a stable oscillation in the drive axis, while an oscillatory motion in the sense axis, which is used to detect an angular rate of rotation, is enabled through the transfers of energy from the drive via the Coriolis force. With the mathematical model depicted as Duffing oscillators, however, by adding a certain coupling among gyroscopes, a similar behavior to a Duffing oscillator is expected, only with more complicated dynamics at a higher dimension. A number of Perturbation methods have popularly been carried out, to seek for a general asymptotic solution of typical Duffing oscillators. In this work as an overview, the two-time scale Perturbation expansion is asymptotically applied on the uni-directionally coupled vibratory gyroscopes to find an analytical solution which is then compared to the numerical one.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Stottmeister, Alexander, E-mail: alexander.stottmeister@gravity.fau.de; Thiemann, Thomas, E-mail: thomas.thiemann@gravity.fau.de

In this article, the second of three, we discuss and develop the basis of a Weyl quantisation for compact Lie groups aiming at loop quantum gravity-type models. This Weyl quantisation may serve as the main mathematical tool to implement the program of space adiabatic perturbation theory in such models. As we already argued in our first article, space adiabatic perturbation theory offers an ideal framework to overcome the obstacles that hinder the direct implementation of the conventional Born-Oppenheimer approach in the canonical formulation of loop quantum gravity. Additionally, we conjecture the existence of a new form of the Segal-Bargmann-Hall “coherentmore » state” transform for compact Lie groups G, which we prove for G = U(1){sup n} and support by numerical evidence for G = SU(2). The reason for conjoining this conjecture with the main topic of this article originates in the observation that the coherent state transform can be used as a basic building block of a coherent state quantisation (Berezin quantisation) for compact Lie groups G. But, as Weyl and Berezin quantisation for ℝ{sup 2d} are intimately related by heat kernel evolution, it is natural to ask whether a similar connection exists for compact Lie groups as well. Moreover, since the formulation of space adiabatic perturbation theory requires a (deformation) quantisation as minimal input, we analyse the question to what extent the coherent state quantisation, defined by the Segal-Bargmann-Hall transform, can serve as basis of the former.« less

A mathematical approach to HIV infection dynamics

NASA Astrophysics Data System (ADS)

Ida, A.; Oharu, S.; Oharu, Y.

2007-07-01

In order to obtain a comprehensive form of mathematical models describing nonlinear phenomena such as HIV infection process and AIDS disease progression, it is efficient to introduce a general class of time-dependent evolution equations in such a way that the associated nonlinear operator is decomposed into the sum of a differential operator and a perturbation which is nonlinear in general and also satisfies no global continuity condition. An attempt is then made to combine the implicit approach (usually adapted for convective diffusion operators) and explicit approach (more suited to treat continuous-type operators representing various physiological interactions), resulting in a semi-implicit product formula. Decomposing the operators in this way and considering their individual properties, it is seen that approximation-solvability of the original model is verified under suitable conditions. Once appropriate terms are formulated to describe treatment by antiretroviral therapy, the time-dependence of the reaction terms appears, and such product formula is useful for generating approximate numerical solutions to the governing equations. With this knowledge, a continuous model for HIV disease progression is formulated and physiological interpretations are provided. The abstract theory is then applied to show existence of unique solutions to the continuous model describing the behavior of the HIV virus in the human body and its reaction to treatment by antiretroviral therapy. The product formula suggests appropriate discrete models describing the dynamics of host pathogen interactions with HIV1 and is applied to perform numerical simulations based on the model of the HIV infection process and disease progression. Finally, the results of our numerical simulations are visualized and it is observed that our results agree with medical and physiological aspects.

Non-perturbative methodologies for low-dimensional strongly-correlated systems: From non-Abelian bosonization to truncated spectrum methods

DOE PAGES

James, Andrew J. A.; Konik, Robert M.; Lecheminant, Philippe; ...

2018-02-26

We review two important non-perturbative approaches for extracting the physics of low-dimensional strongly correlated quantum systems. Firstly, we start by providing a comprehensive review of non-Abelian bosonization. This includes an introduction to the basic elements of conformal field theory as applied to systems with a current algebra, and we orient the reader by presenting a number of applications of non-Abelian bosonization to models with large symmetries. We then tie this technique into recent advances in the ability of cold atomic systems to realize complex symme-tries. Secondly, we discuss truncated spectrum methods for the numerical study of systems in one andmore » two dimensions. For one-dimensional systems we provide the reader with considerable insight into the methodology by reviewing canonical applications of the technique to the Ising model (and its variants) and the sine-Gordon model. Following this we review recent work on the development of renormalization groups, both numerical and analytical, that alleviate the effects of truncating the spectrum. Using these technologies, we consider a number of applications to one-dimensional systems: properties of carbon nanotubes, quenches in the Lieb-Liniger model, 1+1D quantum chro-modynamics, as well as Landau-Ginzburg theories. In the final part we move our attention to consider truncated spectrum methods applied to two-dimensional systems. This involves combining truncated spectrum methods with matrix product state algorithms. Lastly, we describe applications of this method to two-dimensional systems of free fermions and the quantum Ising model, including their non-equilibrium dynamics.« less

Non-perturbative methodologies for low-dimensional strongly-correlated systems: From non-Abelian bosonization to truncated spectrum methods

DOE Office of Scientific and Technical Information (OSTI.GOV)

James, Andrew J. A.; Konik, Robert M.; Lecheminant, Philippe

We review two important non-perturbative approaches for extracting the physics of low-dimensional strongly correlated quantum systems. Firstly, we start by providing a comprehensive review of non-Abelian bosonization. This includes an introduction to the basic elements of conformal field theory as applied to systems with a current algebra, and we orient the reader by presenting a number of applications of non-Abelian bosonization to models with large symmetries. We then tie this technique into recent advances in the ability of cold atomic systems to realize complex symme-tries. Secondly, we discuss truncated spectrum methods for the numerical study of systems in one andmore » two dimensions. For one-dimensional systems we provide the reader with considerable insight into the methodology by reviewing canonical applications of the technique to the Ising model (and its variants) and the sine-Gordon model. Following this we review recent work on the development of renormalization groups, both numerical and analytical, that alleviate the effects of truncating the spectrum. Using these technologies, we consider a number of applications to one-dimensional systems: properties of carbon nanotubes, quenches in the Lieb-Liniger model, 1+1D quantum chro-modynamics, as well as Landau-Ginzburg theories. In the final part we move our attention to consider truncated spectrum methods applied to two-dimensional systems. This involves combining truncated spectrum methods with matrix product state algorithms. Lastly, we describe applications of this method to two-dimensional systems of free fermions and the quantum Ising model, including their non-equilibrium dynamics.« less

Wigner phase space distribution via classical adiabatic switching

DOE Office of Scientific and Technical Information (OSTI.GOV)

Bose, Amartya; Makri, Nancy; Department of Physics, University of Illinois, 1110 W. Green Street, Urbana, Illinois 61801

2015-09-21

Evaluation of the Wigner phase space density for systems of many degrees of freedom presents an extremely demanding task because of the oscillatory nature of the Fourier-type integral. We propose a simple and efficient, approximate procedure for generating the Wigner distribution that avoids the computational difficulties associated with the Wigner transform. Starting from a suitable zeroth-order Hamiltonian, for which the Wigner density is available (either analytically or numerically), the phase space distribution is propagated in time via classical trajectories, while the perturbation is gradually switched on. According to the classical adiabatic theorem, each trajectory maintains a constant action if themore » perturbation is switched on infinitely slowly. We show that the adiabatic switching procedure produces the exact Wigner density for harmonic oscillator eigenstates and also for eigenstates of anharmonic Hamiltonians within the Wentzel-Kramers-Brillouin (WKB) approximation. We generalize the approach to finite temperature by introducing a density rescaling factor that depends on the energy of each trajectory. Time-dependent properties are obtained simply by continuing the integration of each trajectory under the full target Hamiltonian. Further, by construction, the generated approximate Wigner distribution is invariant under classical propagation, and thus, thermodynamic properties are strictly preserved. Numerical tests on one-dimensional and dissipative systems indicate that the method produces results in very good agreement with those obtained by full quantum mechanical methods over a wide temperature range. The method is simple and efficient, as it requires no input besides the force fields required for classical trajectory integration, and is ideal for use in quasiclassical trajectory calculations.« less

Non-perturbative methodologies for low-dimensional strongly-correlated systems: From non-Abelian bosonization to truncated spectrum methods

NASA Astrophysics Data System (ADS)

James, Andrew J. A.; Konik, Robert M.; Lecheminant, Philippe; Robinson, Neil J.; Tsvelik, Alexei M.

2018-04-01

We review two important non-perturbative approaches for extracting the physics of low-dimensional strongly correlated quantum systems. Firstly, we start by providing a comprehensive review of non-Abelian bosonization. This includes an introduction to the basic elements of conformal field theory as applied to systems with a current algebra, and we orient the reader by presenting a number of applications of non-Abelian bosonization to models with large symmetries. We then tie this technique into recent advances in the ability of cold atomic systems to realize complex symmetries. Secondly, we discuss truncated spectrum methods for the numerical study of systems in one and two dimensions. For one-dimensional systems we provide the reader with considerable insight into the methodology by reviewing canonical applications of the technique to the Ising model (and its variants) and the sine-Gordon model. Following this we review recent work on the development of renormalization groups, both numerical and analytical, that alleviate the effects of truncating the spectrum. Using these technologies, we consider a number of applications to one-dimensional systems: properties of carbon nanotubes, quenches in the Lieb–Liniger model, 1  +  1D quantum chromodynamics, as well as Landau–Ginzburg theories. In the final part we move our attention to consider truncated spectrum methods applied to two-dimensional systems. This involves combining truncated spectrum methods with matrix product state algorithms. We describe applications of this method to two-dimensional systems of free fermions and the quantum Ising model, including their non-equilibrium dynamics.

Non-perturbative methodologies for low-dimensional strongly-correlated systems: From non-Abelian bosonization to truncated spectrum methods.

PubMed

James, Andrew J A; Konik, Robert M; Lecheminant, Philippe; Robinson, Neil J; Tsvelik, Alexei M

2018-02-26

We review two important non-perturbative approaches for extracting the physics of low-dimensional strongly correlated quantum systems. Firstly, we start by providing a comprehensive review of non-Abelian bosonization. This includes an introduction to the basic elements of conformal field theory as applied to systems with a current algebra, and we orient the reader by presenting a number of applications of non-Abelian bosonization to models with large symmetries. We then tie this technique into recent advances in the ability of cold atomic systems to realize complex symmetries. Secondly, we discuss truncated spectrum methods for the numerical study of systems in one and two dimensions. For one-dimensional systems we provide the reader with considerable insight into the methodology by reviewing canonical applications of the technique to the Ising model (and its variants) and the sine-Gordon model. Following this we review recent work on the development of renormalization groups, both numerical and analytical, that alleviate the effects of truncating the spectrum. Using these technologies, we consider a number of applications to one-dimensional systems: properties of carbon nanotubes, quenches in the Lieb-Liniger model, 1  +  1D quantum chromodynamics, as well as Landau-Ginzburg theories. In the final part we move our attention to consider truncated spectrum methods applied to two-dimensional systems. This involves combining truncated spectrum methods with matrix product state algorithms. We describe applications of this method to two-dimensional systems of free fermions and the quantum Ising model, including their non-equilibrium dynamics.

A Unified Approach for Solving Nonlinear Regular Perturbation Problems

ERIC Educational Resources Information Center

Khuri, S. A.

2008-01-01

This article describes a simple alternative unified method of solving nonlinear regular perturbation problems. The procedure is based upon the manipulation of Taylor's approximation for the expansion of the nonlinear term in the perturbed equation. An essential feature of this technique is the relative simplicity used and the associated unified…

Design Optimization Method for Composite Components Based on Moment Reliability-Sensitivity Criteria

NASA Astrophysics Data System (ADS)

Sun, Zhigang; Wang, Changxi; Niu, Xuming; Song, Yingdong

2017-08-01

In this paper, a Reliability-Sensitivity Based Design Optimization (RSBDO) methodology for the design of the ceramic matrix composites (CMCs) components has been proposed. A practical and efficient method for reliability analysis and sensitivity analysis of complex components with arbitrary distribution parameters are investigated by using the perturbation method, the respond surface method, the Edgeworth series and the sensitivity analysis approach. The RSBDO methodology is then established by incorporating sensitivity calculation model into RBDO methodology. Finally, the proposed RSBDO methodology is applied to the design of the CMCs components. By comparing with Monte Carlo simulation, the numerical results demonstrate that the proposed methodology provides an accurate, convergent and computationally efficient method for reliability-analysis based finite element modeling engineering practice.

Optimal motion planning for collision avoidance of mobile robots in non-stationary environments

NASA Technical Reports Server (NTRS)

Kyriakopoulos, K. J.; Saridis, G. N.

1992-01-01

An optimal control formulation of the problem of collision avoidance of mobile robots moving in general terrains containing moving obstacles is presented. A dynamic model of the mobile robot and the dynamic constraints are derived. Collision avoidance is guaranteed if the minimum distance between the robot and the object is nonzero. A nominal trajectory is assumed to be known from off-line planning. The main idea is to change the velocity along the nominal trajectory so that collisions are avoided. Time consistency with the nominal plan is desirable. A numerical solution of the optimization problem is obtained. A perturbation control type of approach is used to update the optimal plan. Simulation results verify the value of the proposed strategy.

DNA replication stress and cancer chemotherapy.

PubMed

Kitao, Hiroyuki; Iimori, Makoto; Kataoka, Yuki; Wakasa, Takeshi; Tokunaga, Eriko; Saeki, Hiroshi; Oki, Eiji; Maehara, Yoshihiko

2018-02-01

DNA replication is one of the fundamental biological processes in which dysregulation can cause genome instability. This instability is one of the hallmarks of cancer and confers genetic diversity during tumorigenesis. Numerous experimental and clinical studies have indicated that most tumors have experienced and overcome the stresses caused by the perturbation of DNA replication, which is also referred to as DNA replication stress (DRS). When we consider therapeutic approaches for tumors, it is important to exploit the differences in DRS between tumor and normal cells. In this review, we introduce the current understanding of DRS in tumors and discuss the underlying mechanism of cancer therapy from the aspect of DRS. © 2017 The Authors. Cancer Science published by John Wiley & Sons Australia, Ltd on behalf of Japanese Cancer Association.

Self-consistent approach to many-body localization and subdiffusion

NASA Astrophysics Data System (ADS)

Prelovšek, P.; Herbrych, J.

2017-07-01

An analytical theory, based on the perturbative treatment of the disorder and extended into a self-consistent set of equations for the dynamical density correlations, is developed and applied to the prototype one-dimensional model of many-body localization. Results show a qualitative agreement with the numerically obtained dynamical structure factor in the whole range of frequencies and wave vectors, as well as across the transition to nonergodic behavior. The theory reveals the singular nature of the one-dimensional problem, whereby on the ergodic side the dynamics is subdiffusive with dynamical conductivity σ (ω ) ∝|ω| α , i.e., with vanishing dc limit σ0=0 and α <1 varying with disorder, while we get α >1 in the localized phase.

Integral sliding mode-based attitude coordinated tracking for spacecraft formation with communication delays

NASA Astrophysics Data System (ADS)

Zhang, Jian; Hu, Qinglei; Xie, Wenbo

2017-11-01

This paper investigates the attitude coordinated tracking control for a group of rigid spacecraft under directed communication topology, in which inertia uncertainties, external disturbances, input saturation and constant time-delays between the formation members are handled. Initially, the nominal system with communication delays is studied. A delay-dependent controller is proposed by using Lyapunov-Krasovskii function and sufficient condition for system stability is derived. Then, an integral sliding manifold is designed and adaptive control approach is employed to deal with the total perturbation. Meanwhile, the boundary layer method is introduced to alleviate the unexpected chattering as system trajectories cross the switching surface. Finally, numerical simulation results are presented to validate the effectiveness and robustness of the proposed control strategy.

A Dynamical Systems Approach to the Design of the Science Orbit Around Europa

NASA Technical Reports Server (NTRS)

Gomez, Gerard; Lara, Martin; Russell, Ryan P.

2006-01-01

The science orbit for a future mission to Europa requires low eccentricity, low altitude, and high inclination. However, high inclination orbits around planetary satellites are unstable due to third-body perturbations. Without control, the orbiter impacts Europa after few weeks. To minimize control, a tour over the stable-unstable, averaged manifolds of unstable frozen orbits has been suggested. We proceed with the unaveraged equations and study the manifolds of unstable orbits that are periodic in a rotating frame attached to Europa. Massive numerical computation helps in understanding the unstable dynamics close to Europa, and, thus, in selecting long lifetime high inclination orbits. A final test of a selected set of initial conditions on a high fidelity, ephemeris model, validate the results.

Statistics of the stochastically forced Lorenz attractor by the Fokker-Planck equation and cumulant expansions.

PubMed

Allawala, Altan; Marston, J B

2016-11-01

We investigate the Fokker-Planck description of the equal-time statistics of the three-dimensional Lorenz attractor with additive white noise. The invariant measure is found by computing the zero (or null) mode of the linear Fokker-Planck operator as a problem of sparse linear algebra. Two variants are studied: a self-adjoint construction of the linear operator and the replacement of diffusion with hyperdiffusion. We also access the low-order statistics of the system by a perturbative expansion in equal-time cumulants. A comparison is made to statistics obtained by the standard approach of accumulation via direct numerical simulation. Theoretical and computational aspects of the Fokker-Planck and cumulant expansion methods are discussed.

A semi-analytical method of computation of oceanic tidal perturbations in the motion of artificial satellites

NASA Technical Reports Server (NTRS)

Musen, P.

1973-01-01

The method of expansion of the satellite's perturbations, as caused by the oceanic tides, into Fourier series is discussed. The coefficients of the expansion are purely numerical and peculiar to each particular satellite. Such a method is termed as semi-analytical in celestial mechanics. Gaussian form of the differential equations for variation of elements, with the right hand sides averaged over the orbit of the satellite, is convenient to use with the semi-analytical expansion.

Application of the enhanced homotopy perturbation method to solve the fractional-order Bagley-Torvik differential equation

NASA Astrophysics Data System (ADS)

Zolfaghari, M.; Ghaderi, R.; Sheikhol Eslami, A.; Ranjbar, A.; Hosseinnia, S. H.; Momani, S.; Sadati, J.

2009-10-01

The enhanced homotopy perturbation method (EHPM) is applied for finding improved approximate solutions of the well-known Bagley-Torvik equation for three different cases. The main characteristic of the EHPM is using a stabilized linear part, which guarantees the stability and convergence of the overall solution. The results are finally compared with the Adams-Bashforth-Moulton numerical method, the Adomian decomposition method (ADM) and the fractional differential transform method (FDTM) to verify the performance of the EHPM.

The hydrogen atom in D = 3 - 2ɛ dimensions

NASA Astrophysics Data System (ADS)

Adkins, Gregory S.

2018-06-01

The nonrelativistic hydrogen atom in D = 3 - 2 ɛ dimensions is the reference system for perturbative schemes used in dimensionally regularized nonrelativistic effective field theories to describe hydrogen-like atoms. Solutions to the D-dimensional Schrödinger-Coulomb equation are given in the form of a double power series. Energies and normalization integrals are obtained numerically and also perturbatively in terms of ɛ. The utility of the series expansion is demonstrated by the calculation of the divergent expectation value <(V‧)2 >.

THE FIRST FERMI IN A HIGH ENERGY NUCLEAR COLLISION.

DOE Office of Scientific and Technical Information (OSTI.GOV)

KRASNITZ,A.

1999-08-09

At very high energies, weak coupling, non-perturbative methods can be used to study classical gluon production in nuclear collisions. One observes in numerical simulations that after an initial formation time, the produced partons are on shell, and their subsequent evolution can be studied using transport theory. At the initial formation time, a simple non-perturbative relation exists between the energy and number densities of the produced partons, and a scale determined by the saturated parton density in the nucleus.

Planar hydrodynamic instability computations and experiments with rugby-shaped hohlraums at the Omega laser

NASA Astrophysics Data System (ADS)

Vandenboomgaerde, M.; Liberatore, S.; Galmiche, D.; Casner, A.; Huser, G.; Jadaud, J. P.; Villette, B.

2008-05-01

Implosion of inertial confinement fusion (ICF) capsule is very sensitive to the growth of sphericity perturbations. The control of the feeding of such perturbations and their transport ('feedthrough') through the ablator is a key point to reach ignition. Since 2002 [1, 2], experiments have been designed and performed on the Omega laser facility in order to study these phenomena in planar geometry. A new 'rugby shaped' hohlraum was used [3, 4]. We present experimental results and comparisons with numerical simulations.

A numerical study on weak-dissipative two-mode perturbed Burgers' and Ostrovsky models: right-left moving waves

NASA Astrophysics Data System (ADS)

Jaradat, Imad; Alquran, Marwan; Ali, Mohammed

2018-04-01

The purpose of this study is threefold. First, it derives newly developed two-mode nonlinear equations, two-mode perturbed Burgers' and two-mode Ostrovsky models. Second, it investigates the values of the nonlinearity and dispersion parameters that support the existence of two right-left (R-L) moving wave solutions to these models. Finally, it provides a graphical analysis of the "two-mode" concept and the impact of its phase velocity on the field function.

Numerical investigation of non-perturbative kinetic effects of energetic particles on toroidicity-induced Alfvén eigenmodes in tokamaks and stellarators

DOE Office of Scientific and Technical Information (OSTI.GOV)

Slaby, Christoph; Könies, Axel; Kleiber, Ralf

2016-09-15

The resonant interaction of shear Alfvén waves with energetic particles is investigated numerically in tokamak and stellarator geometry using a non-perturbative MHD-kinetic hybrid approach. The focus lies on toroidicity-induced Alfvén eigenmodes (TAEs), which are most easily destabilized by a fast-particle population in fusion plasmas. While the background plasma is treated within the framework of an ideal-MHD theory, the drive of the fast particles, as well as Landau damping of the background plasma, is modelled using the drift-kinetic Vlasov equation without collisions. Building on analytical theory, a fast numerical tool, STAE-K, has been developed to solve the resulting eigenvalue problem usingmore » a Riccati shooting method. The code, which can be used for parameter scans, is applied to tokamaks and the stellarator Wendelstein 7-X. High energetic-ion pressure leads to large growth rates of the TAEs and to their conversion into kinetically modified TAEs and kinetic Alfvén waves via continuum interaction. To better understand the physics of this conversion mechanism, the connections between TAEs and the shear Alfvén wave continuum are examined. It is shown that, when energetic particles are present, the continuum deforms substantially and the TAE frequency can leave the continuum gap. The interaction of the TAE with the continuum leads to singularities in the eigenfunctions. To further advance the physical model and also to eliminate the MHD continuum together with the singularities in the eigenfunctions, a fourth-order term connected to radiative damping has been included. The radiative damping term is connected to non-ideal effects of the bulk plasma and introduces higher-order derivatives to the model. Thus, it has the potential to substantially change the nature of the solution. For the first time, the fast-particle drive, Landau damping, continuum damping, and radiative damping have been modelled together in tokamak- as well as in stellarator geometry.« less

A normal mode treatment of semi-diurnal body tides on an aspherical, rotating and anelastic Earth

NASA Astrophysics Data System (ADS)

Lau, Harriet C. P.; Yang, Hsin-Ying; Tromp, Jeroen; Mitrovica, Jerry X.; Latychev, Konstantin; Al-Attar, David

2015-08-01

Normal mode treatments of the Earth's body tide response were developed in the 1980s to account for the effects of Earth rotation, ellipticity, anelasticity and resonant excitation within the diurnal band. Recent space-geodetic measurements of the Earth's crustal displacement in response to luni-solar tidal forcings have revealed geographical variations that are indicative of aspherical deep mantle structure, thus providing a novel data set for constraining deep mantle elastic and density structure. In light of this, we make use of advances in seismic free oscillation literature to develop a new, generalized normal mode theory for the tidal response within the semi-diurnal and long-period tidal band. Our theory involves a perturbation method that permits an efficient calculation of the impact of aspherical structure on the tidal response. In addition, we introduce a normal mode treatment of anelasticity that is distinct from both earlier work in body tides and the approach adopted in free oscillation seismology. We present several simple numerical applications of the new theory. First, we compute the tidal response of a spherically symmetric, non-rotating, elastic and isotropic Earth model and demonstrate that our predictions match those based on standard Love number theory. Second, we compute perturbations to this response associated with mantle anelasticity and demonstrate that the usual set of seismic modes adopted for this purpose must be augmented by a family of relaxation modes to accurately capture the full effect of anelasticity on the body tide response. Finally, we explore aspherical effects including rotation and we benchmark results from several illustrative case studies of aspherical Earth structure against independent finite-volume numerical calculations of the semi-diurnal body tide response. These tests confirm the accuracy of the normal mode methodology to at least the level of numerical error in the finite-volume predictions. They also demonstrate that full coupling of normal modes, rather than group coupling, is necessary for accurate predictions of the body tide response.

Local convertibility of the ground state of the perturbed toric code

NASA Astrophysics Data System (ADS)

Santra, Siddhartha; Hamma, Alioscia; Cincio, Lukasz; Subasi, Yigit; Zanardi, Paolo; Amico, Luigi

2014-12-01

We present analytical and numerical studies of the behavior of the α -Renyi entropies in the toric code in presence of several types of perturbations aimed at studying the simulability of these perturbations to the parent Hamiltonian using local operations and classical communications (LOCC)—a property called local convertibility. In particular, the derivatives, with respect to the perturbation parameter, present different signs for different values of α within the topological phase. From the information-theoretic point of view, this means that such ground states cannot be continuously deformed within the topological phase by means of catalyst assisted local operations and classical communications (LOCC). Such LOCC differential convertibility is on the other hand always possible in the trivial disordered phase. The non-LOCC convertibility is remarkable because it can be computed on a system whose size is independent of correlation length. This method can therefore constitute an experimentally feasible witness of topological order.

Perturbative stability of SFT-based cosmological models

DOE Office of Scientific and Technical Information (OSTI.GOV)

Galli, Federico; Koshelev, Alexey S., E-mail: fgalli@tena4.vub.ac.be, E-mail: alexey.koshelev@vub.ac.be

2011-05-01

We review the appearance of multiple scalar fields in linearized SFT based cosmological models with a single non-local scalar field. Some of these local fields are canonical real scalar fields and some are complex fields with unusual coupling. These systems only admit numerical or approximate analysis. We introduce a modified potential for multiple scalar fields that makes the system exactly solvable in the cosmological context of Friedmann equations and at the same time preserves the asymptotic behavior expected from SFT. The main part of the paper consists of the analysis of inhomogeneous cosmological perturbations in this system. We show numericallymore » that perturbations corresponding to the new type of complex fields always vanish. As an example of application of this model we consider an explicit construction of the phantom divide crossing and prove the perturbative stability of this process at the linear order. The issue of ghosts and ways to resolve it are briefly discussed.« less

The equilibrium and stability of the gaseous component of the galaxy, 2

NASA Technical Reports Server (NTRS)

Kellman, S. A.

1971-01-01

A time-independent, linear, plane and axially-symmetric stability analysis was performed on a self-gravitating, plane-parallel, isothermal layer of nonmagnetic, nonrotating gas. The gas layer was immersed in a plane-stratified field isothermal layer of stars which supply a self-consistent gravitational field. Only the gaseous component was perturbed. Expressions were derived for the perturbed gas potential and perturbed gas density that satisfied both the Poisson and hydrostatic equilibrium equations. The equation governing the size of the perturbations in the mid-plane was found to be analogous to the one-dimensional time-independent Schrodinger equation for a particle bound by a potential well, and with similar boundary conditions. The radius of the neutral state was computed numerically and compared with the Jeans' and Ledoux radius. The inclusion of a rigid stellar component increased the Ledoux radius, though only slightly. Isodensity contours of the neutrual or marginally unstable state were constructed.

Double power series method for approximating cosmological perturbations

NASA Astrophysics Data System (ADS)

Wren, Andrew J.; Malik, Karim A.

2017-04-01

We introduce a double power series method for finding approximate analytical solutions for systems of differential equations commonly found in cosmological perturbation theory. The method was set out, in a noncosmological context, by Feshchenko, Shkil' and Nikolenko (FSN) in 1966, and is applicable to cases where perturbations are on subhorizon scales. The FSN method is essentially an extension of the well known Wentzel-Kramers-Brillouin (WKB) method for finding approximate analytical solutions for ordinary differential equations. The FSN method we use is applicable well beyond perturbation theory to solve systems of ordinary differential equations, linear in the derivatives, that also depend on a small parameter, which here we take to be related to the inverse wave-number. We use the FSN method to find new approximate oscillating solutions in linear order cosmological perturbation theory for a flat radiation-matter universe. Together with this model's well-known growing and decaying Mészáros solutions, these oscillating modes provide a complete set of subhorizon approximations for the metric potential, radiation and matter perturbations. Comparison with numerical solutions of the perturbation equations shows that our approximations can be made accurate to within a typical error of 1%, or better. We also set out a heuristic method for error estimation. A Mathematica notebook which implements the double power series method is made available online.

Feature-based data assimilation in geophysics

NASA Astrophysics Data System (ADS)

Morzfeld, Matthias; Adams, Jesse; Lunderman, Spencer; Orozco, Rafael

2018-05-01

Many applications in science require that computational models and data be combined. In a Bayesian framework, this is usually done by defining likelihoods based on the mismatch of model outputs and data. However, matching model outputs and data in this way can be unnecessary or impossible. For example, using large amounts of steady state data is unnecessary because these data are redundant. It is numerically difficult to assimilate data in chaotic systems. It is often impossible to assimilate data of a complex system into a low-dimensional model. As a specific example, consider a low-dimensional stochastic model for the dipole of the Earth's magnetic field, while other field components are ignored in the model. The above issues can be addressed by selecting features of the data, and defining likelihoods based on the features, rather than by the usual mismatch of model output and data. Our goal is to contribute to a fundamental understanding of such a feature-based approach that allows us to assimilate selected aspects of data into models. We also explain how the feature-based approach can be interpreted as a method for reducing an effective dimension and derive new noise models, based on perturbed observations, that lead to computationally efficient solutions. Numerical implementations of our ideas are illustrated in four examples.

SYNCHRONIZATION OF HETEROGENEOUS OSCILLATORS UNDER NETWORK MODIFICATIONS: PERTURBATION AND OPTIMIZATION OF THE SYNCHRONY ALIGNMENT FUNCTION

PubMed Central

Taylor, Dane; Skardal, Per Sebastian; Sun, Jie

2016-01-01

Synchronization is central to many complex systems in engineering physics (e.g., the power-grid, Josephson junction circuits, and electro-chemical oscillators) and biology (e.g., neuronal, circadian, and cardiac rhythms). Despite these widespread applications—for which proper functionality depends sensitively on the extent of synchronization—there remains a lack of understanding for how systems can best evolve and adapt to enhance or inhibit synchronization. We study how network modifications affect the synchronization properties of network-coupled dynamical systems that have heterogeneous node dynamics (e.g., phase oscillators with non-identical frequencies), which is often the case for real-world systems. Our approach relies on a synchrony alignment function (SAF) that quantifies the interplay between heterogeneity of the network and of the oscillators and provides an objective measure for a system’s ability to synchronize. We conduct a spectral perturbation analysis of the SAF for structural network modifications including the addition and removal of edges, which subsequently ranks the edges according to their importance to synchronization. Based on this analysis, we develop gradient-descent algorithms to efficiently solve optimization problems that aim to maximize phase synchronization via network modifications. We support these and other results with numerical experiments. PMID:27872501

Simple, explicitly time-dependent, and regular solutions of the linearized vacuum Einstein equations in Bondi-Sachs coordinates

NASA Astrophysics Data System (ADS)

Mädler, Thomas

2013-05-01

Perturbations of the linearized vacuum Einstein equations in the Bondi-Sachs formulation of general relativity can be derived from a single master function with spin weight two, which is related to the Weyl scalar Ψ0, and which is determined by a simple wave equation. By utilizing a standard spin representation of tensors on a sphere and two different approaches to solve the master equation, we are able to determine two simple and explicitly time-dependent solutions. Both solutions, of which one is asymptotically flat, comply with the regularity conditions at the vertex of the null cone. For the asymptotically flat solution we calculate the corresponding linearized perturbations, describing all multipoles of spin-2 waves that propagate on a Minkowskian background spacetime. We also analyze the asymptotic behavior of this solution at null infinity using a Penrose compactification and calculate the Weyl scalar Ψ4. Because of its simplicity, the asymptotically flat solution presented here is ideally suited for test bed calculations in the Bondi-Sachs formulation of numerical relativity. It may be considered as a sibling of the Bergmann-Sachs or Teukolsky-Rinne solutions, on spacelike hypersurfaces, for a metric adapted to null hypersurfaces.

CMB anisotropies from patchy reionisation and diffuse Sunyaev-Zel'dovich effects

DOE Office of Scientific and Technical Information (OSTI.GOV)

Fidler, Christian; Ringeval, Christophe, E-mail: christophe.ringeval@uclouvain.be, E-mail: christian.fidler@uclouvain.be

Anisotropies in the Cosmic Microwave Background (CMB) can be induced during the later stages of cosmic evolution, and in particular during and after the Epoch of Reionisation. Inhomogeneities in the ionised fraction, but also in the baryon density, in the velocity fields and in the gravitational potentials are expected to generate correlated CMB perturbations. We present a complete relativistic treatment of all these effects, up to second order in perturbation theory, that we solve using the numerical Boltzmann code (\\SONG). The physical origin and relevance of all second order terms are carefully discussed. In addition to collisional and gravitational contributions,more » we identify the diffuse analogue of the blurring and kinetic Sunyaev-Zel'dovich (SZ) effects. Our approach naturally includes the correlations between the imprint from patchy reionisation and the diffuse SZ effects thereby allowing us to derive reliable estimates of the induced temperature and polarisation CMB angular power spectra. In particular, we show that the B -modes generated at intermediate length-scales (ℓ ≅ 100) have the same amplitude as the B -modes coming from primordial gravitational waves with a tensor-to-scalar ratio r =10{sup −4}.« less

Full counting statistics of conductance for disordered systems

NASA Astrophysics Data System (ADS)

Fu, Bin; Zhang, Lei; Wei, Yadong; Wang, Jian

2017-09-01

Quantum transport is a stochastic process in nature. As a result, the conductance is fully characterized by its average value and fluctuations, i.e., characterized by full counting statistics (FCS). Since disorders are inevitable in nanoelectronic devices, it is important to understand how FCS behaves in disordered systems. The traditional approach dealing with fluctuations or cumulants of conductance uses diagrammatic perturbation expansion of the Green's function within coherent potential approximation (CPA), which is extremely complicated especially for high order cumulants. In this paper, we develop a theoretical formalism based on nonequilibrium Green's function by directly taking the disorder average on the generating function of FCS of conductance within CPA. This is done by mapping the problem into higher dimensions so that the functional dependence of generating a function on the Green's function becomes linear and the diagrammatic perturbation expansion is not needed anymore. Our theory is very simple and allows us to calculate cumulants of conductance at any desired order efficiently. As an application of our theory, we calculate the cumulants of conductance up to fifth order for disordered systems in the presence of Anderson and binary disorders. Our numerical results of cumulants of conductance show remarkable agreement with that obtained by the brute force calculation.

Straightening the Hierarchical Staircase for Basis Set Extrapolations: A Low-Cost Approach to High-Accuracy Computational Chemistry

NASA Astrophysics Data System (ADS)

Varandas, António J. C.

2018-04-01

Because the one-electron basis set limit is difficult to reach in correlated post-Hartree-Fock ab initio calculations, the low-cost route of using methods that extrapolate to the estimated basis set limit attracts immediate interest. The situation is somewhat more satisfactory at the Hartree-Fock level because numerical calculation of the energy is often affordable at nearly converged basis set levels. Still, extrapolation schemes for the Hartree-Fock energy are addressed here, although the focus is on the more slowly convergent and computationally demanding correlation energy. Because they are frequently based on the gold-standard coupled-cluster theory with single, double, and perturbative triple excitations [CCSD(T)], correlated calculations are often affordable only with the smallest basis sets, and hence single-level extrapolations from one raw energy could attain maximum usefulness. This possibility is examined. Whenever possible, this review uses raw data from second-order Møller-Plesset perturbation theory, as well as CCSD, CCSD(T), and multireference configuration interaction methods. Inescapably, the emphasis is on work done by the author's research group. Certain issues in need of further research or review are pinpointed.

Spatio-temporal correlations in models of collective motion ruled by different dynamical laws.

PubMed

Cavagna, Andrea; Conti, Daniele; Giardina, Irene; Grigera, Tomas S; Melillo, Stefania; Viale, Massimiliano

2016-11-15

Information transfer is an essential factor in determining the robustness of biological systems with distributed control. The most direct way to study the mechanisms ruling information transfer is to experimentally observe the propagation across the system of a signal triggered by some perturbation. However, this method may be inefficient for experiments in the field, as the possibilities to perturb the system are limited and empirical observations must rely on natural events. An alternative approach is to use spatio-temporal correlations to probe the information transfer mechanism directly from the spontaneous fluctuations of the system, without the need to have an actual propagating signal on record. Here we test this method on models of collective behaviour in their deeply ordered phase by using ground truth data provided by numerical simulations in three dimensions. We compare two models characterized by very different dynamical equations and information transfer mechanisms: the classic Vicsek model, describing an overdamped noninertial dynamics and the inertial spin model, characterized by an underdamped inertial dynamics. By using dynamic finite-size scaling, we show that spatio-temporal correlations are able to distinguish unambiguously the diffusive information transfer mechanism of the Vicsek model from the linear mechanism of the inertial spin model.

Energy management of three-dimensional minimum-time intercept. [for aircraft flight optimization

NASA Technical Reports Server (NTRS)

Kelley, H. J.; Cliff, E. M.; Visser, H. G.

1985-01-01

A real-time computer algorithm to control and optimize aircraft flight profiles is described and applied to a three-dimensional minimum-time intercept mission. The proposed scheme has roots in two well known techniques: singular perturbations and neighboring-optimal guidance. Use of singular-perturbation ideas is made in terms of the assumed trajectory-family structure. A heading/energy family of prestored point-mass-model state-Euler solutions is used as the baseline in this scheme. The next step is to generate a near-optimal guidance law that will transfer the aircraft to the vicinity of this reference family. The control commands fed to the autopilot (bank angle and load factor) consist of the reference controls plus correction terms which are linear combinations of the altitude and path-angle deviations from reference values, weighted by a set of precalculated gains. In this respect the proposed scheme resembles neighboring-optimal guidance. However, in contrast to the neighboring-optimal guidance scheme, the reference control and state variables as well as the feedback gains are stored as functions of energy and heading in the present approach. Some numerical results comparing open-loop optimal and approximate feedback solutions are presented.

Rethinking the logistic approach for population dynamics of mutualistic interactions.

PubMed

García-Algarra, Javier; Galeano, Javier; Pastor, Juan Manuel; Iriondo, José María; Ramasco, José J

2014-12-21

Mutualistic communities have an internal structure that makes them resilient to external perturbations. Late research has focused on their stability and the topology of the relations between the different organisms to explain the reasons of the system robustness. Much less attention has been invested in analyzing the systems dynamics. The main population models in use are modifications of the r-K formulation of logistic equation with additional terms to account for the benefits produced by the interspecific interactions. These models have shortcomings as the so-called r-K formulation diverges under some conditions. In this work, we introduce a model for population dynamics under mutualism that preserves the original logistic formulation. It is mathematically simpler than the widely used type II models, although it shows similar complexity in terms of fixed points and stability of the dynamics. We perform an analytical stability analysis and numerical simulations to study the model behavior in general interaction scenarios including tests of the resilience of its dynamics under external perturbations. Despite its simplicity, our results indicate that the model dynamics shows an important richness that can be used to gain further insights in the dynamics of mutualistic communities. Copyright © 2014 Elsevier Ltd. All rights reserved.