Architectures for Parafermionic Topological Matter in Two Dimensions

NASA Astrophysics Data System (ADS)

Burrello, Michele; van Heck, Bernard; Cobanera, Emilio

2013-03-01

Recent proposals exploit edge modes of fractional topological insulators (FTIs), induced superconducting pairing, and back-scattering to realize one-dimensional systems of parafermions. We extend these proposals to two dimensions and analyze the effect of the superconducting islands' charging energy on the parafermions they host. We focus on two two-dimensional architectures, the tile and stripe configurations, characterized by different distributions of FTIs and derive the associated parafermionic effective Hamiltonians. The tile model realizes the Z2 m toric code in low-order perturbation theory and hence possesses full topological quantum order. By exploiting dualities, we obtain the phase diagram and generalized order parameters for both the tile and stripe models of parafermions. This work was supported by the Dutch Science Foundation NWO/FOM and an ERC Advanced Investigator grant.

Fast large scale structure perturbation theory using one-dimensional fast Fourier transforms

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

Schmittfull, Marcel; Vlah, Zvonimir; McDonald, Patrick

The usual fluid equations describing the large-scale evolution of mass density in the universe can be written as local in the density, velocity divergence, and velocity potential fields. As a result, the perturbative expansion in small density fluctuations, usually written in terms of convolutions in Fourier space, can be written as a series of products of these fields evaluated at the same location in configuration space. Based on this, we establish a new method to numerically evaluate the 1-loop power spectrum (i.e., Fourier transform of the 2-point correlation function) with one-dimensional fast Fourier transforms. This is exact and a fewmore » orders of magnitude faster than previously used numerical approaches. Numerical results of the new method are in excellent agreement with the standard quadrature integration method. This fast model evaluation can in principle be extended to higher loop order where existing codes become painfully slow. Our approach follows by writing higher order corrections to the 2-point correlation function as, e.g., the correlation between two second-order fields or the correlation between a linear and a third-order field. These are then decomposed into products of correlations of linear fields and derivatives of linear fields. In conclusion, the method can also be viewed as evaluating three-dimensional Fourier space convolutions using products in configuration space, which may also be useful in other contexts where similar integrals appear.« less

Fast large scale structure perturbation theory using one-dimensional fast Fourier transforms

DOE PAGES

Schmittfull, Marcel; Vlah, Zvonimir; McDonald, Patrick

2016-05-01

The usual fluid equations describing the large-scale evolution of mass density in the universe can be written as local in the density, velocity divergence, and velocity potential fields. As a result, the perturbative expansion in small density fluctuations, usually written in terms of convolutions in Fourier space, can be written as a series of products of these fields evaluated at the same location in configuration space. Based on this, we establish a new method to numerically evaluate the 1-loop power spectrum (i.e., Fourier transform of the 2-point correlation function) with one-dimensional fast Fourier transforms. This is exact and a fewmore » orders of magnitude faster than previously used numerical approaches. Numerical results of the new method are in excellent agreement with the standard quadrature integration method. This fast model evaluation can in principle be extended to higher loop order where existing codes become painfully slow. Our approach follows by writing higher order corrections to the 2-point correlation function as, e.g., the correlation between two second-order fields or the correlation between a linear and a third-order field. These are then decomposed into products of correlations of linear fields and derivatives of linear fields. In conclusion, the method can also be viewed as evaluating three-dimensional Fourier space convolutions using products in configuration space, which may also be useful in other contexts where similar integrals appear.« less

1D Resonance line Broadened Quasilinear (RBQ1D) code for fast ion Alfvenic relaxations and its validations

NASA Astrophysics Data System (ADS)

Gorelenkov, Nikolai; Duarte, Vinicius; Podesta, Mario

2017-10-01

The performance of the burning plasma can be limited by the requirements to confine the superalfvenic fusion products which are capable of resonating with the Alfvénic eigenmodes (AEs). The effect of AEs on fast ions is evaluated using the quasi-linear approach [Berk et al., Ph.Plasmas'96] generalized for this problem recently [Duarte et al., Ph.D.'17]. The generalization involves the resonance line broadened interaction regions with the diffusion coefficient prescribed to find the evolution of the velocity distribution function. The baseline eigenmode structures are found using the NOVA-K code perturbatively [Gorelenkov et al., Ph.Plasmas'99]. A RBQ1D code allowing the diffusion in radial direction is presented here. The wave particle interaction can be reduced to one-dimensional dynamics where for the Alfvénic modes typically the particle kinetic energy is nearly constant. Hence to a good approximation the Quasi-Linear (QL) diffusion equation only contains derivatives in the angular momentum. The diffusion equation is then one dimensional that is efficiently solved simultaneously for all particles with the equation for the evolution of the wave angular momentum. The RBQ1D is validated against recent DIIID results [Collins et al., PRL'16]. Supported by the US Department of Energy under DE-AC02-09CH11466.

Study of the hydrodynamics of the formation of flows caused by the interaction of a shock wave with two-dimensional density perturbations on the Iskra-5 laser facility

NASA Astrophysics Data System (ADS)

Babanov, A. V.; Barinov, M. A.; Barinov, S. P.; Garanin, R. V.; Zhidkov, N. V.; Kalmykov, N. A.; Kovalenko, V. P.; Kokorin, S. N.; Pinegin, A. V.; Solomatina, E. Yu.; Solomatin, I. I.; Suslov, N. A.

2017-03-01

The hydrodynamics of the flow formation due to the interaction of a shock wave with two-dimensional density perturbations is experimentally investigated on the Iskra-5 laser facility. Shadow images of a jet arising as a result of the impact of a shock wave (formed by a soft X-ray pulse from a target-illuminator) on a flat aluminium target with a blind cylindrical cavity are recorded in experiments with point-like X-ray backlighting having a photon energy of ~4.5 keV. The sizes and mass of the jet ejected from the aluminium cavity by this shock wave are estimated. The experimental data are compared with the results of numerical simulation of the jet formation and dynamics according to the two-dimensional MID-ND2D code.

Shape design sensitivity analysis and optimization of three dimensional elastic solids using geometric modeling and automatic regridding. Ph.D. Thesis

NASA Technical Reports Server (NTRS)

Yao, Tse-Min; Choi, Kyung K.

1987-01-01

An automatic regridding method and a three dimensional shape design parameterization technique were constructed and integrated into a unified theory of shape design sensitivity analysis. An algorithm was developed for general shape design sensitivity analysis of three dimensional eleastic solids. Numerical implementation of this shape design sensitivity analysis method was carried out using the finite element code ANSYS. The unified theory of shape design sensitivity analysis uses the material derivative of continuum mechanics with a design velocity field that represents shape change effects over the structural design. Automatic regridding methods were developed by generating a domain velocity field with boundary displacement method. Shape design parameterization for three dimensional surface design problems was illustrated using a Bezier surface with boundary perturbations that depend linearly on the perturbation of design parameters. A linearization method of optimization, LINRM, was used to obtain optimum shapes. Three examples from different engineering disciplines were investigated to demonstrate the accuracy and versatility of this shape design sensitivity analysis method.

Nonlinear oscillations of compact stars in the vicinity of the maximum mass configuration

NASA Astrophysics Data System (ADS)

Brillante, A.; Mishustin, I. N.

2015-07-01

We solve the dynamical GR equations for the spherically symmetric evolution of compact stars in the vicinity of the maximum mass, for which instability sets in according to linear perturbation theory. The calculations are done with the analytical Zeldovich-like EOS P=a≤ft(ρ-ρ_0\\right) and with the TM1 parametrisation of the RMF model. The initial configurations for the dynamical calculations are represented by spherical stars with equilibrium density profile, which are perturbed by either i) an artificially added inward velocity field proportional to the radial coordinate, or ii) a rarefaction corresponding to a static and expanded star. These configurations are evolved using a one-dimensional GR hydro code for ideal and barotropic fluids. Depending on the initial conditions we obtain either stable oscillations or the collapse to a black hole. The minimal amplitude of the perturbation, needed to trigger gravitational collapse is evaluated. The approximate independence of this energy on the type of perturbation is pointed out. At the threshold we find type-I critical behaviour for all stellar models considered and discuss the dependence of the time scaling exponent on the baryon mass and EOS.

Non-linear hydrodynamical evolution of rotating relativistic stars: numerical methods and code tests

NASA Astrophysics Data System (ADS)

Font, José A.; Stergioulas, Nikolaos; Kokkotas, Kostas D.

2000-04-01

We present numerical hydrodynamical evolutions of rapidly rotating relativistic stars, using an axisymmetric, non-linear relativistic hydrodynamics code. We use four different high-resolution shock-capturing (HRSC) finite-difference schemes (based on approximate Riemann solvers) and compare their accuracy in preserving uniformly rotating stationary initial configurations in long-term evolutions. Among these four schemes, we find that the third-order piecewise parabolic method scheme is superior in maintaining the initial rotation law in long-term evolutions, especially near the surface of the star. It is further shown that HRSC schemes are suitable for the evolution of perturbed neutron stars and for the accurate identification (via Fourier transforms) of normal modes of oscillation. This is demonstrated for radial and quadrupolar pulsations in the non-rotating limit, where we find good agreement with frequencies obtained with a linear perturbation code. The code can be used for studying small-amplitude or non-linear pulsations of differentially rotating neutron stars, while our present results serve as testbed computations for three-dimensional general-relativistic evolution codes.

Depletion Calculations Based on Perturbations. Application to the Study of a Rep-Like Assembly at Beginning of Cycle with TRIPOLI-4®.

NASA Astrophysics Data System (ADS)

Dieudonne, Cyril; Dumonteil, Eric; Malvagi, Fausto; M'Backé Diop, Cheikh

2014-06-01

For several years, Monte Carlo burnup/depletion codes have appeared, which couple Monte Carlo codes to simulate the neutron transport to deterministic methods, which handle the medium depletion due to the neutron flux. Solving Boltzmann and Bateman equations in such a way allows to track fine 3-dimensional effects and to get rid of multi-group hypotheses done by deterministic solvers. The counterpart is the prohibitive calculation time due to the Monte Carlo solver called at each time step. In this paper we present a methodology to avoid the repetitive and time-expensive Monte Carlo simulations, and to replace them by perturbation calculations: indeed the different burnup steps may be seen as perturbations of the isotopic concentration of an initial Monte Carlo simulation. In a first time we will present this method, and provide details on the perturbative technique used, namely the correlated sampling. In a second time the implementation of this method in the TRIPOLI-4® code will be discussed, as well as the precise calculation scheme a meme to bring important speed-up of the depletion calculation. Finally, this technique will be used to calculate the depletion of a REP-like assembly, studied at beginning of its cycle. After having validated the method with a reference calculation we will show that it can speed-up by nearly an order of magnitude standard Monte-Carlo depletion codes.

Two-Dimensional Description of Potential Perturbations to the Ozone Layer Due to NOx and H2O Aircraft Emissions,

DTIC Science & Technology

1984-01-01

Perturbations to the Ozone 01aparUNM fe"MAWLayer Due to NOX and H20Arrf Emissions20Araf Office of Environment and Energy Washington, D.C. 20591 uTI LEOU George F...and Energy 14. Sponsorng Agency Code Washington, D.C. 20591 15. Supplementary Notes 1 . Abstroct . This report is divided into two parts. Part I is an...obtained from the work of Louis, et al. (1974), who calculated the circulation patterns by solving the continuity and energy equations using compiled

Analysis of the Effect of Electron Density Perturbations Generated by Gravity Waves on HF Communication Links

NASA Astrophysics Data System (ADS)

Fagre, M.; Elias, A. G.; Chum, J.; Cabrera, M. A.

2017-12-01

In the present work, ray tracing of high frequency (HF) signals in ionospheric disturbed conditions is analyzed, particularly in the presence of electron density perturbations generated by gravity waves (GWs). The three-dimensional numerical ray tracing code by Jones and Stephenson, based on Hamilton's equations, which is commonly used to study radio propagation through the ionosphere, is used. An electron density perturbation model is implemented to this code based upon the consideration of atmospheric GWs generated at a height of 150 km in the thermosphere and propagating up into the ionosphere. The motion of the neutral gas at these altitudes induces disturbances in the background plasma which affects HF signals propagation. To obtain a realistic model of GWs in order to analyze the propagation and dispersion characteristics, a GW ray tracing method with kinematic viscosity and thermal diffusivity was applied. The IRI-2012, HWM14 and NRLMSISE-00 models were incorporated to assess electron density, wind velocities, neutral temperature and total mass density needed for the ray tracing codes. Preliminary results of gravity wave effects on ground range and reflection height are presented for low-mid latitude ionosphere.

Solar Filament Longitudinal Oscillations along a Magnetic Field Tube with Two Dips

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

Zhou Yu-Hao; Zhang Li-Yue; Ouyang, Y.

Large-amplitude longitudinal oscillations of solar filaments have been observed and explored for more than ten years. Previous studies are mainly based on the one-dimensional rigid flux tube model with a single magnetic dip. However, it has been noted that there might be two magnetic dips, and hence two threads, along one magnetic field line. Following previous work, we intend to investigate the kinematics of the filament longitudinal oscillations when two threads are magnetically connected, which is done by solving one-dimensional radiative hydrodynamic equations with the numerical code MPI-AMRVAC. Two different types of perturbations are considered, and the difference from previousmore » works resulting from the interaction of the two filament threads is investigated. We find that even with the inclusion of the thread–thread interaction, the oscillation period is modified weakly, by at most 20% compared to the traditional pendulum model with one thread. However, the damping timescale is significantly affected by the thread–thread interaction. Hence, we should take it into account when applying the consistent seismology to the filaments where two threads are magnetically connected.« less

2D Implosion Simulations with a Kinetic Particle Code

NASA Astrophysics Data System (ADS)

Sagert, Irina; Even, Wesley; Strother, Terrance

2017-10-01

Many problems in laboratory and plasma physics are subject to flows that move between the continuum and the kinetic regime. We discuss two-dimensional (2D) implosion simulations that were performed using a Monte Carlo kinetic particle code. The application of kinetic transport theory is motivated, in part, by the occurrence of non-equilibrium effects in inertial confinement fusion (ICF) capsule implosions, which cannot be fully captured by hydrodynamics simulations. Kinetic methods, on the other hand, are able to describe both, continuum and rarefied flows. We perform simple 2D disk implosion simulations using one particle species and compare the results to simulations with the hydrodynamics code RAGE. The impact of the particle mean-free-path on the implosion is also explored. In a second study, we focus on the formation of fluid instabilities from induced perturbations. I.S. acknowledges support through the Director's fellowship from Los Alamos National Laboratory. This research used resources provided by the LANL Institutional Computing Program.

Multidimensional Modeling of Atmospheric Effects and Surface Heterogeneities on Remote Sensing

NASA Technical Reports Server (NTRS)

Gerstl, S. A. W.; Simmer, C.; Zardecki, A. (Principal Investigator)

1985-01-01

The overall goal of this project is to establish a modeling capability that allows a quantitative determination of atmospheric effects on remote sensing including the effects of surface heterogeneities. This includes an improved understanding of aerosol and haze effects in connection with structural, angular, and spatial surface heterogeneities. One important objective of the research is the possible identification of intrinsic surface or canopy characteristics that might be invariant to atmospheric perturbations so that they could be used for scene identification. Conversely, an equally important objective is to find a correction algorithm for atmospheric effects in satellite-sensed surface reflectances. The technical approach is centered around a systematic model and code development effort based on existing, highly advanced computer codes that were originally developed for nuclear radiation shielding applications. Computational techniques for the numerical solution of the radiative transfer equation are adapted on the basis of the discrete-ordinates finite-element method which proved highly successful for one and two-dimensional radiative transfer problems with fully resolved angular representation of the radiation field.

Reconstruction of a three-dimensional, transonic rotor flow field from holographic interferogram data

NASA Technical Reports Server (NTRS)

Kittleson, John K.; Yu, Yung H.

1987-01-01

Holographic interferometry and computerized aided tomography (CAT) are used to determine the transonic velocity field of a model rotor blade in hover. A pulsed ruby laser recorded 40 interferograms with a 2 ft dia view field near the model rotor blade tip operating at a tip Mach number of 0.90. After digitizing the interferograms and extracting the fringe order functions, the data are transferred to a CAT code. The CAT code then calculates the perturbation velocity in several planes above the blade surface. The values from the holography-CAT method compare favorably with previously obtained numerical computations in most locations near the blade tip. The results demonstrate the technique's potential for three dimensional transonic rotor flow studies.

Three-dimensional dynamical and chemical modelling of the upper atmosphere

NASA Technical Reports Server (NTRS)

Prinn, R. G.; Alyea, F. N.; Cunnold, D. M.

1976-01-01

Progress in coding a 3-D upper atmospheric model and in modeling the ozone perturbation resulting from the shuttle booster exhaust is reported. A time-dependent version of a 2-D model was studied and the sulfur cycle in the stratosphere was investigated. The role of meteorology in influencing stratospheric composition measurements was also studied.

Solwnd: A 3D Compressible MHD Code for Solar Wind Studies. Version 1.0: Cartesian Coordinates

NASA Technical Reports Server (NTRS)

Deane, Anil E.

1996-01-01

Solwnd 1.0 is a three-dimensional compressible MHD code written in Fortran for studying the solar wind. Time-dependent boundary conditions are available. The computational algorithm is based on Flux Corrected Transport and the code is based on the existing code of Zalesak and Spicer. The flow considered is that of shear flow with incoming flow that perturbs this base flow. Several test cases corresponding to pressure balanced magnetic structures with velocity shear flow and various inflows including Alfven waves are presented. Version 1.0 of solwnd considers a rectangular Cartesian geometry. Future versions of solwnd will consider a spherical geometry. Some discussions of this issue is presented.

An initial investigation into methods of computing transonic aerodynamic sensitivity coefficients

NASA Technical Reports Server (NTRS)

Carlson, Leland A.

1994-01-01

The primary accomplishments of the project are as follows: (1) Using the transonic small perturbation equation as a flowfield model, the project demonstrated that the quasi-analytical method could be used to obtain aerodynamic sensitivity coefficients for airfoils at subsonic, transonic, and supersonic conditions for design variables such as Mach number, airfoil thickness, maximum camber, angle of attack, and location of maximum camber. It was established that the quasi-analytical approach was an accurate method for obtaining aerodynamic sensitivity derivatives for airfoils at transonic conditions and usually more efficient than the finite difference approach. (2) The usage of symbolic manipulation software to determine the appropriate expressions and computer coding associated with the quasi-analytical method for sensitivity derivatives was investigated. Using the three dimensional fully conservative full potential flowfield model, it was determined that symbolic manipulation along with a chain rule approach was extremely useful in developing a combined flowfield and quasi-analytical sensitivity derivative code capable of considering a large number of realistic design variables. (3) Using the three dimensional fully conservative full potential flowfield model, the quasi-analytical method was applied to swept wings (i.e. three dimensional) at transonic flow conditions. (4) The incremental iterative technique has been applied to the three dimensional transonic nonlinear small perturbation flowfield formulation, an equivalent plate deflection model, and the associated aerodynamic and structural discipline sensitivity equations; and coupled aeroelastic results for an aspect ratio three wing in transonic flow have been obtained.

Topological entanglement entropy of fracton stabilizer codes

NASA Astrophysics Data System (ADS)

Ma, Han; Schmitz, A. T.; Parameswaran, S. A.; Hermele, Michael; Nandkishore, Rahul M.

2018-03-01

Entanglement entropy provides a powerful characterization of two-dimensional gapped topological phases of quantum matter, intimately tied to their description by topological quantum field theories (TQFTs). Fracton topological orders are three-dimensional gapped topologically ordered states of matter that lack a TQFT description. We show that three-dimensional fracton phases are nevertheless characterized, at least partially, by universal structure in the entanglement entropy of their ground-state wave functions. We explicitly compute the entanglement entropy for two archetypal fracton models, the "X-cube model" and "Haah's code," and demonstrate the existence of a nonlocal contribution that scales linearly in subsystem size. We show via Schrieffer-Wolff transformations that this piece of the entanglement entropy of fracton models is robust against arbitrary local perturbations of the Hamiltonian. Finally, we argue that these results may be extended to characterize localization-protected fracton topological order in excited states of disordered fracton models.

Numerical optimization of three-dimensional coils for NSTX-U

NASA Astrophysics Data System (ADS)

Lazerson, S. A.; Park, J.-K.; Logan, N.; Boozer, A.

2015-10-01

A tool for the calculation of optimal three-dimensional (3D) perturbative magnetic fields in tokamaks has been developed. The IPECOPT code builds upon the stellarator optimization code STELLOPT to allow for optimization of linear ideal magnetohydrodynamic perturbed equilibrium (IPEC). This tool has been applied to NSTX-U equilibria, addressing which fields are the most effective at driving NTV torques. The NTV torque calculation is performed by the PENT code. Optimization of the normal field spectrum shows that fields with n  =  1 character can drive a large core torque. It is also shown that fields with n  =  3 features are capable of driving edge torque and some core torque. Coil current optimization (using the planned in-vessel and existing RWM coils) on NSTX-U suggest the planned coils set is adequate for core and edge torque control. Comparison between error field correction experiments on DIII-D and the optimizer show good agreement. Notice: This manuscript has been authored by Princeton University under Contract Number DE-AC02-09CH11466 with the U.S. Department of Energy. The publisher, by accepting the article for publication acknowledges, that the United States Government retains a non-exclusive,paid-up, irrevocable, world-wide license to publish or reproduce the published form of this manuscript, or allow others to do so, for United States Government purposes.

Signatures of extra dimensions in gravitational waves from black hole quasinormal modes

NASA Astrophysics Data System (ADS)

Chakraborty, Sumanta; Chakravarti, Kabir; Bose, Sukanta; SenGupta, Soumitra

2018-05-01

In this work, we have derived the evolution equation for gravitational perturbation in four-dimensional spacetime in the presence of a spatial extra dimension. The evolution equation is derived by perturbing the effective gravitational field equations on the four-dimensional spacetime, which inherits nontrivial higher-dimensional effects. Note that this is different from the perturbation of the five-dimensional gravitational field equations that exist in the literature and possess quantitatively new features. The gravitational perturbation has further been decomposed into a purely four-dimensional part and another piece that depends on extra dimensions. The four-dimensional gravitational perturbation now admits massive propagating degrees of freedom, owing to the existence of higher dimensions. We have also studied the influence of these massive propagating modes on the quasinormal mode frequencies, signaling the higher-dimensional nature of the spacetime, and have contrasted these massive modes with the massless modes in general relativity. Surprisingly, it turns out that the massive modes experience damping much smaller than that of the massless modes in general relativity and may even dominate over and above the general relativity contribution if one observes the ringdown phase of a black hole merger event at sufficiently late times. Furthermore, the whole analytical framework has been supplemented by the fully numerical Cauchy evolution problem, as well. In this context, we have shown that, except for minute details, the overall features of the gravitational perturbations are captured both in the Cauchy evolution as well as in the analysis of quasinormal modes. The implications on observations of black holes with LIGO and proposed space missions such as LISA are also discussed.

A finsler perturbation of the Poincaré metric

NASA Astrophysics Data System (ADS)

Rutz, Solange F.; McCarthy, Patrick J.

1993-02-01

One method of gaining some insight into Finsler geomety is that of studying small Finsler perturbations of Riemannian metrics. We consider here the the standard two-dimensional upper half plane Poincaré metric, for which the geodesics are semi-circles and vertical lines. The effect of a simple Finsler perturbation on these geodesics is given by an explicit computation of the perturbed geodesics.

Classifying the Quantum Phases of Matter

DTIC Science & Technology

2015-01-01

Kim related entanglement entropy to topological storage of quantum information [8]. Michalakis et al. showed that a particle-like excitation spectrum...Perturbative analysis of topological entanglement entropy from conditional independence, Phys. Rev. B 86, 254116 (2012), arXiv:1210.2360. [3] I. Kim...symmetries or long-range entanglement ), (2) elucidating the properties of three-dimensional quantum codes (in particular those which admit no string-like

Linus cycle calculations including plasma transport and resistive flux loss in an incompressible liner

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

Quimby, D.C.; Hoffman, A.L.; Vlases, G.C.

1980-08-01

In the LINUS fusion reactor concept, a rotating liquid metal liner is used for reversible mechanical compression of thermonuclear plasmas, where a vacuum field buffer zone is used between the plasma and wall to reduce transport losses. A one-dimensional plasma transport and burn code, including incompressible liner dynamics with heat transfer and temperature dependent flux diffusion in the liquid metal, is used to model LINUS cycles. The effects of compressibility are treated as a perturbation. Numerical coefficients are derived for simple LINUS scaling laws. The particular case of plasma contact with the liquid metal is studied to determine the effectmore » on LINUS performance.« less

Source Distribution Method for Unsteady One-Dimensional Flows With Small Mass, Momentum, and Heat Addition and Small Area Variation

NASA Technical Reports Server (NTRS)

Mirels, Harold

1959-01-01

A source distribution method is presented for obtaining flow perturbations due to small unsteady area variations, mass, momentum, and heat additions in a basic uniform (or piecewise uniform) one-dimensional flow. First, the perturbations due to an elemental area variation, mass, momentum, and heat addition are found. The general solution is then represented by a spatial and temporal distribution of these elemental (source) solutions. Emphasis is placed on discussing the physical nature of the flow phenomena. The method is illustrated by several examples. These include the determination of perturbations in basic flows consisting of (1) a shock propagating through a nonuniform tube, (2) a constant-velocity piston driving a shock, (3) ideal shock-tube flows, and (4) deflagrations initiated at a closed end. The method is particularly applicable for finding the perturbations due to relatively thin wall boundary layers.

Chaotic attractors of relaxation oscillators

NASA Astrophysics Data System (ADS)

Guckenheimer, John; Wechselberger, Martin; Young, Lai-Sang

2006-03-01

We develop a general technique for proving the existence of chaotic attractors for three-dimensional vector fields with two time scales. Our results connect two important areas of dynamical systems: the theory of chaotic attractors for discrete two-dimensional Henon-like maps and geometric singular perturbation theory. Two-dimensional Henon-like maps are diffeomorphisms that limit on non-invertible one-dimensional maps. Wang and Young formulated hypotheses that suffice to prove the existence of chaotic attractors in these families. Three-dimensional singularly perturbed vector fields have return maps that are also two-dimensional diffeomorphisms limiting on one-dimensional maps. We describe a generic mechanism that produces folds in these return maps and demonstrate that the Wang-Young hypotheses are satisfied. Our analysis requires a careful study of the convergence of the return maps to their singular limits in the Ck topology for k >= 3. The theoretical results are illustrated with a numerical study of a variant of the forced van der Pol oscillator.

Extending SIESTA capabilities: removing field-periodic and stellarator symmetric limitations

NASA Astrophysics Data System (ADS)

Cook, C. R.; Hirshman, S. P.; Sanchez, R.; Anderson, D. T.

2011-10-01

SIESTA is a three-dimensional magnetohydrodynamics equilibrium code capable of resolving magnetic islands in toroidal plasma confinement devices. Currently SIESTA assumes that plasma perturbations, and thus also magnetic islands, are field-periodic. This limitation is being removed from the code by allowing the displacement toroidal mode number to not be restricted to multiples of the number of field periods. Extending SIESTA in this manner will allow larger, lower-order resonant islands to form in devices such as CTH. An example of a non-field-periodic perturbation in CTH will be demonstrated. Currently the code also operates in a stellarator-symmetric fashion in which an ``up-down'' symmetry is present at some toroidal angle. Nearly all of the current tokamaks (and ITER in the future) operate with a divertor and as such do not possess stellarator symmetry. Removal of this symmetry restriction requires including both sine and cosine terms in the Fourier expansion for the geometry of the device and the fields contained within. The current status of this extension of the code will be discussed, along with the method of implementation. U.S. DOE Contract No. DE-AC05-00OR22725 with UT-Battelle, LLC.

Rarefaction-driven Rayleigh–Taylor instability. Part 2. Experiments and simulations in the nonlinear regime

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

Morgan, R. V.; Cabot, W. H.; Greenough, J. A.

Experiments and large eddy simulation (LES) were performed to study the development of the Rayleigh–Taylor instability into the saturated, nonlinear regime, produced between two gases accelerated by a rarefaction wave. Single-mode two-dimensional, and single-mode three-dimensional initial perturbations were introduced on the diffuse interface between the two gases prior to acceleration. The rarefaction wave imparts a non-constant acceleration, and a time decreasing Atwood number,more » $$A=(\\unicode[STIX]{x1D70C}_{2}-\\unicode[STIX]{x1D70C}_{1})/(\\unicode[STIX]{x1D70C}_{2}+\\unicode[STIX]{x1D70C}_{1})$$, where$$\\unicode[STIX]{x1D70C}_{2}$$and$$\\unicode[STIX]{x1D70C}_{1}$$are the densities of the heavy and light gas, respectively. Experiments and simulations are presented for initial Atwood numbers of$A=0.49$$,$$A=0.63$$,$$A=0.82$$and$$A=0.94$$. Nominally two-dimensional (2-D) experiments (initiated with nearly 2-D perturbations) and 2-D simulations are observed to approach an intermediate-time velocity plateau that is in disagreement with the late-time velocity obtained from the incompressible model of Goncharov (Phys. Rev. Lett., vol. 88, 2002, 134502). Reacceleration from an intermediate velocity is observed for 2-D bubbles in large wavenumber,$$k=2\\unicode[STIX]{x03C0}/\\unicode[STIX]{x1D706}=0.247~\\text{mm}^{-1}$$, experiments and simulations, where$$\\unicode[STIX]{x1D706}$is the wavelength of the initial perturbation. At moderate Atwood numbers, the bubble and spike velocities approach larger values than those predicted by Goncharov’s model. These late-time velocity trends are predicted well by numerical simulations using the LLNL Miranda code, and by the 2009 model of Mikaelian (Phys. Fluids., vol. 21, 2009, 024103) that extends Layzer type models to variable acceleration and density. Large Atwood number experiments show a delayed roll up, and exhibit a free-fall like behaviour. Finally, experiments initiated with three-dimensional perturbations tend to agree better with models and a simulation using the LLNL Ares code initiated with an axisymmetric rather than Cartesian symmetry.« less

Rarefaction-driven Rayleigh–Taylor instability. Part 2. Experiments and simulations in the nonlinear regime

DOE PAGES

Morgan, R. V.; Cabot, W. H.; Greenough, J. A.; ...

2018-01-12

Experiments and large eddy simulation (LES) were performed to study the development of the Rayleigh–Taylor instability into the saturated, nonlinear regime, produced between two gases accelerated by a rarefaction wave. Single-mode two-dimensional, and single-mode three-dimensional initial perturbations were introduced on the diffuse interface between the two gases prior to acceleration. The rarefaction wave imparts a non-constant acceleration, and a time decreasing Atwood number,more » $$A=(\\unicode[STIX]{x1D70C}_{2}-\\unicode[STIX]{x1D70C}_{1})/(\\unicode[STIX]{x1D70C}_{2}+\\unicode[STIX]{x1D70C}_{1})$$, where$$\\unicode[STIX]{x1D70C}_{2}$$and$$\\unicode[STIX]{x1D70C}_{1}$$are the densities of the heavy and light gas, respectively. Experiments and simulations are presented for initial Atwood numbers of$A=0.49$$,$$A=0.63$$,$$A=0.82$$and$$A=0.94$$. Nominally two-dimensional (2-D) experiments (initiated with nearly 2-D perturbations) and 2-D simulations are observed to approach an intermediate-time velocity plateau that is in disagreement with the late-time velocity obtained from the incompressible model of Goncharov (Phys. Rev. Lett., vol. 88, 2002, 134502). Reacceleration from an intermediate velocity is observed for 2-D bubbles in large wavenumber,$$k=2\\unicode[STIX]{x03C0}/\\unicode[STIX]{x1D706}=0.247~\\text{mm}^{-1}$$, experiments and simulations, where$$\\unicode[STIX]{x1D706}$is the wavelength of the initial perturbation. At moderate Atwood numbers, the bubble and spike velocities approach larger values than those predicted by Goncharov’s model. These late-time velocity trends are predicted well by numerical simulations using the LLNL Miranda code, and by the 2009 model of Mikaelian (Phys. Fluids., vol. 21, 2009, 024103) that extends Layzer type models to variable acceleration and density. Large Atwood number experiments show a delayed roll up, and exhibit a free-fall like behaviour. Finally, experiments initiated with three-dimensional perturbations tend to agree better with models and a simulation using the LLNL Ares code initiated with an axisymmetric rather than Cartesian symmetry.« less

On the Absolutely Continuous Spectrum of Stark Operators

NASA Astrophysics Data System (ADS)

Perelman, Galina

The stability of the absolutely continuous spectrum of the one-dimensional Stark operator under perturbations of the potential is discussed. The focus is on proving this stability under minimal assumptions on smoothness of the perturbation. A general criterion is presented together with some applications. These include the case of periodic perturbations where we show that any perturbation vL1()∩H-1/2() preserves the a.c. spectrum.

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.

Controlling spiral waves and turbulent states in cardiac tissue by traveling wave perturbations

NASA Astrophysics Data System (ADS)

Wang, Peng-Ye; Xie, Ping

2000-03-01

We propose a traveling wave perturbation method to control the spatiotemporal dynamics in cardiac tissue. With a two-variable model we demonstrate 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 traveling wave state in the two-dimensional case. An experimental scheme is suggested which may provide a new design for a cardiac defibrillator.

Topological phases in two-dimensional arrays of parafermionic zero modes

NASA Astrophysics Data System (ADS)

Burrello, M.; van Heck, B.; Cobanera, E.

2013-05-01

It has recently been realized that zero modes with projective non-Abelian statistics, generalizing the notion of Majorana bound states, may exist at the interface between a superconductor and a ferromagnet along the edge of a fractional topological insulator (FTI). Here, we study two-dimensional architectures of these non-Abelian zero modes, whose interactions are generated by the charging and Josephson energies of the superconductors. We derive low-energy Hamiltonians for two different arrays of FTIs on the plane, revealing an interesting interplay between the real-space geometry of the system and its topological properties. On the one hand, in a geometry where the length of the FTI edges is independent on the system size, the array has a topologically ordered phase, giving rise to a qudit toric code Hamiltonian in perturbation theory. On the other hand, in a geometry where the length of the edges scales with system size, we find an exact duality to an Abelian lattice gauge theory and no topological order.

Three- and Two- Dimensional Simulations of Re-shock Experiments at High Energy Densities at the National Ignition Facility

NASA Astrophysics Data System (ADS)

Wang, Ping; Raman, Kumar; MacLaren, Stephan; Huntington, Channing; Nagel, Sabrina

2016-10-01

We present simulations of recent high-energy-density (HED) re-shock experiments on the National Ignition Facility (NIF). The experiments study the Rayleigh-Taylor (RT) and Richtmyer-Meshkov (RM) instability growth that occurs after successive shocks transit a sinusoidally-perturbed interface between materials of different densities. The shock tube is driven at one or both ends using indirect-drive laser cavities or hohlraums. X-ray area-backlit imaging is used to visualize the growth at different times. Our simulations are done with the three-dimensional, radiation hydrodynamics code ARES, developed at LLNL. We show the instabilitygrowth rate, inferred from the experimental radiographs, agrees well with our 2D and 3D simulations. We also discuss some 3D geometrical effects, suggested by our simulations, which could deteriorate the images at late times, unless properly accounted for in the experiment design. Work supported by U.S. Department of Energy under Contract DE- AC52-06NA27279. LLNL-ABS-680789.

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

Favorite, Jeffrey A.

SENSMG is a tool for computing first-order sensitivities of neutron reaction rates, reaction-rate ratios, leakage, k eff, and α using the PARTISN multigroup discrete-ordinates code. SENSMG computes sensitivities to all of the transport cross sections and data (total, fission, nu, chi, and all scattering moments), two edit cross sections (absorption and capture), and the density for every isotope and energy group. It also computes sensitivities to the mass density for every material and derivatives with respect to all interface locations. The tool can be used for one-dimensional spherical (r) and two-dimensional cylindrical (r-z) geometries. The tool can be used formore » fixed-source and eigenvalue problems. The tool implements Generalized Perturbation Theory (GPT) as discussed by Williams and Stacey. Section II of this report describes the theory behind adjoint-based sensitivities, gives the equations that SENSMG solves, and defines the sensitivities that are output. Section III describes the user interface, including the input file and command line options. Section IV describes the output. Section V gives some notes about the coding that may be of interest. Section VI discusses verification, which is ongoing. Section VII lists needs and ideas for future work. Appendix A lists all of the input files whose results are presented in Sec. VI.« less

Comparison of two-dimensional and quasi-one-dimensional scramjet models by the example of VAG experiment

NASA Astrophysics Data System (ADS)

Seleznev, R. K.

2017-02-01

In the paper two-dimensional and quasi-one dimensional models for scramjet combustion chamber are described. Comparison of the results of calculations for the two-dimensional and quasi-one dimensional code by the example of VAG experiment are presented.

Effect of rotation zero-crossing on single-fluid plasma response to three-dimensional magnetic perturbations

NASA Astrophysics Data System (ADS)

Lyons, B. C.; Ferraro, N. M.; Paz-Soldan, C.; Nazikian, R.; Wingen, A.

2017-04-01

In order to understand the effect of rotation on the response of a plasma to three-dimensional magnetic perturbations, we perform a systematic scan of the zero-crossing of the rotation profile in a DIII-D ITER-similar shape equilibrium using linear, time-independent modeling with the M3D-C1 extended magnetohydrodynamics code. We confirm that the local resonant magnetic field generally increases as the rotation decreases at a rational surface. Multiple peaks in the resonant field are observed near rational surfaces, however, and the maximum resonant field does not always correspond to zero rotation at the surface. Furthermore, we show that non-resonant current can be driven at zero-crossings not aligned with rational surfaces if there is sufficient shear in the rotation profile there, leading to amplification of near-resonant Fourier harmonics of the perturbed magnetic field and a decrease in the far-off-resonant harmonics. The quasilinear electromagnetic torque induced by this non-resonant plasma response provides drive to flatten the rotation, possibly allowing for increased transport in the pedestal by the destabilization of turbulent modes. In addition, this torque acts to drive the rotation zero-crossing to dynamically stable points near rational surfaces, which would allow for increased resonant penetration. By one or both of these mechanisms, this torque may play an important role in bifurcations into suppression of edge-localized modes. Finally, we discuss how these changes to the plasma response could be detected by tokamak diagnostics. In particular, we show that the changes to the resonant field discussed here have a significant impact on the external perturbed magnetic field, which should be observable by magnetic sensors on the high-field side of tokamaks but not on the low-field side. In addition, TRIP3D-MAFOT simulations show that none of the changes to the plasma response described here substantially affects the divertor footprint structure.

Effect of rotation zero-crossing on single-fluid plasma response to three-dimensional magnetic perturbations

DOE PAGES

Lyons, Brendan C.; Ferraro, Nathaniel M.; Paz-Soldan, Carlos A.; ...

2017-02-14

In order to understand the effect of rotation on the plasma's response to three-dimensional magnetic perturbations, we perform a systematic scan of the zero-crossing of the rotation profile in a DIII-D ITER-similar shape equilibrium using linear, time-independent modeling with the M3D-C1 extended magnetohydrodynamics code. We confirm that the local resonant magnetic field generally increases as the rotation decreases at a rational surface. Multiple peaks in the resonant field are observed near rational surfaces, however, and the maximum resonant field does not always correspond to zero rotation at the surface. Furthermore, we show that non-resonant current can be driven at zero-more » crossings not aligned with rational surfaces if there is sufficient shear in the rotation profile there, leading to an amplification of near-resonant Fourier harmonics of the perturbed magnetic field and a decrease in the far-off -resonant harmonics. The quasilinear electromagnetic torque induced by this non-resonant plasma response provides drive to flatten the rotation, possibly allowing for increased transport in the pedestal by the destabilization of turbulent modes. In addition, this torque acts to drive the rotation zero-crossing to dynamically stable points near rational surfaces, which would allow for increased resonant penetration. By one or both of these mechanisms, this torque may play an important role in bifurcations into ELM suppression. Finally, we discuss how these changes to the plasma response could be detected by tokamak diagnostics. In particular, we show that the changes to the resonant field discussed here have a significant impact on the external perturbed magnetic field, which should be observable by magnetic sensors on the high-field side of tokamaks, but not on the low-field side. In addition, TRIP3D-MAFOT simulations show that none of the changes to the plasma response described here substantially affects the divertor footprint structure.« less

Synthetic Microwave Imaging Reflectometry diagnostic using 3D FDTD Simulations

NASA Astrophysics Data System (ADS)

Kruger, Scott; Jenkins, Thomas; Smithe, David; King, Jacob; Nimrod Team Team

2017-10-01

Microwave Imaging Reflectometry (MIR) has become a standard diagnostic for understanding tokamak edge perturbations, including the edge harmonic oscillations in QH mode operation. These long-wavelength perturbations are larger than the normal turbulent fluctuation levels and thus normal analysis of synthetic signals become more difficult. To investigate, we construct a synthetic MIR diagnostic for exploring density fluctuation amplitudes in the tokamak plasma edge by using the three-dimensional, full-wave FDTD code Vorpal. The source microwave beam for the diagnostic is generated and refelected at the cutoff surface that is distorted by 2D density fluctuations in the edge plasma. Synthetic imaging optics at the detector can be used to understand the fluctuation and background density profiles. We apply the diagnostic to understand the fluctuations in edge plasma density during QH-mode activity in the DIII-D tokamak, as modeled by the NIMROD code. This work was funded under DOE Grant Number DE-FC02-08ER54972.

VENTURE: a code block for solving multigroup neutronics problems applying the finite-difference diffusion-theory approximation to neutron transport

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

Vondy, D.R.; Fowler, T.B.; Cunningham, G.W.

1975-10-01

The computer code block VENTURE, designed to solve multigroup neutronics problems with application of the finite-difference diffusion-theory approximation to neutron transport (or alternatively simple P$sub 1$) in up to three- dimensional geometry is described. A variety of types of problems may be solved: the usual eigenvalue problem, a direct criticality search on the buckling, on a reciprocal velocity absorber (prompt mode), or on nuclide concentrations, or an indirect criticality search on nuclide concentrations, or on dimensions. First- order perturbation analysis capability is available at the macroscopic cross section level. (auth)

A New Interpretation of Effects of the Probabilistic Delayed Start on the One-Dimensional Traffic Flow

NASA Astrophysics Data System (ADS)

Ishibashi, Yoshihiro; Fukui, Minoru

2018-03-01

The effect of the probabilistic delayed start on the one-dimensional traffic flow is investigated on the basis of several models. Analogy with the degeneracy of the states and its resolution, as well as that with the mathematical procedures adopted for them, is utilized. The perturbation is assumed to be proportional to the probability of the delayed start, and the perturbation function is determined so that imposed conditions are fulfilled. The obtained formulas coincide with those previously derived on the basis of the mean-field analyses of the Nagel-Schreckenberg and Fukui-Ishibashi models, and reproduce the cellular automaton simulation results.

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.

Identification of Stochastically Perturbed Autonomous Systems from Temporal Sequences of Probability Density Functions

NASA Astrophysics Data System (ADS)

Nie, Xiaokai; Luo, Jingjing; Coca, Daniel; Birkin, Mark; Chen, Jing

2018-03-01

The paper introduces a method for reconstructing one-dimensional iterated maps that are driven by an external control input and subjected to an additive stochastic perturbation, from sequences of probability density functions that are generated by the stochastic dynamical systems and observed experimentally.

Black hole perturbation under a 2 +2 decomposition in the action

NASA Astrophysics Data System (ADS)

Ripley, Justin L.; Yagi, Kent

2018-01-01

Black hole perturbation theory is useful for studying the stability of black holes and calculating ringdown gravitational waves after the collision of two black holes. Most previous calculations were carried out at the level of the field equations instead of the action. In this work, we compute the Einstein-Hilbert action to quadratic order in linear metric perturbations about a spherically symmetric vacuum background in Regge-Wheeler gauge. Using a 2 +2 splitting of spacetime, we expand the metric perturbations into a sum over scalar, vector, and tensor spherical harmonics, and dimensionally reduce the action to two dimensions by integrating over the two sphere. We find that the axial perturbation degree of freedom is described by a two-dimensional massive vector action, and that the polar perturbation degree of freedom is described by a two-dimensional dilaton massive gravity action. Varying the dimensionally reduced actions, we rederive covariant and gauge-invariant master equations for the axial and polar degrees of freedom. Thus, the two-dimensional massive vector and massive gravity actions we derive by dimensionally reducing the perturbed Einstein-Hilbert action describe the dynamics of a well-studied physical system: the metric perturbations of a static black hole. The 2 +2 formalism we present can be generalized to m +n -dimensional spacetime splittings, which may be useful in more generic situations, such as expanding metric perturbations in higher dimensional gravity. We provide a self-contained presentation of m +n formalism for vacuum spacetime splittings.

Closed coronal structures. V - Gasdynamic models of flaring loops and comparison with SMM observations

NASA Technical Reports Server (NTRS)

Peres, G.; Serio, S.; Vaiana, G.; Acton, L.; Leibacher, J.; Rosner, R.; Pallavicini, R.

1983-01-01

A time-dependent one-dimensional code incorporating energy, momentum and mass conservation equations, and taking the entire solar atmospheric structure into account, is used to investigate the hydrodynamic response of confined magnetic structures to strong heating perturbations. Model calculation results are compared with flare observations which include the light curves of spectral lines formed over a wide range of coronal flare temperatures, as well as determinations of Doppler shifts for the high temperature plasma. It is shown that the numerical simulation predictions are in good overall agreement with the observed flare coronal plasma evolution, correctly reproducing the temporal profile of X-ray spectral lines and their relative intensities. The predicted upflow velocities support the interpretation of the blueshifts as due to evaporation of chromospheric material.

Numerical model for learning concepts of streamflow simulation

USGS Publications Warehouse

DeLong, L.L.; ,

1993-01-01

Numerical models are useful for demonstrating principles of open-channel flow. Such models can allow experimentation with cause-and-effect relations, testing concepts of physics and numerical techniques. Four PT is a numerical model written primarily as a teaching supplement for a course in one-dimensional stream-flow modeling. Four PT options particularly useful in training include selection of governing equations, boundary-value perturbation, and user-programmable constraint equations. The model can simulate non-trivial concepts such as flow in complex interconnected channel networks, meandering channels with variable effective flow lengths, hydraulic structures defined by unique three-parameter relations, and density-driven flow.The model is coded in FORTRAN 77, and data encapsulation is used extensively to simplify maintenance and modification and to enhance the use of Four PT modules by other programs and programmers.

Perturbations of linear delay differential equations at the verge of instability.

PubMed

Lingala, N; Namachchivaya, N Sri

2016-06-01

The characteristic equation for a linear delay differential equation (DDE) has countably infinite roots on the complex plane. This paper considers linear DDEs that are on the verge of instability, i.e., a pair of roots of the characteristic equation lies on the imaginary axis of the complex plane and all other roots have negative real parts. It is shown that when small noise perturbations are present, the probability distribution of the dynamics can be approximated by the probability distribution of a certain one-dimensional stochastic differential equation (SDE) without delay. This is advantageous because equations without delay are easier to simulate and one-dimensional SDEs are analytically tractable. When the perturbations are also linear, it is shown that the stability depends on a specific complex number. The theory is applied to study oscillators with delayed feedback. Some errors in other articles that use multiscale approach are pointed out.

Three-dimensional modeling of plasma edge transport and divertor fluxes during application of resonant magnetic perturbations on ITER

NASA Astrophysics Data System (ADS)

Schmitz, O.; Becoulet, M.; Cahyna, P.; Evans, T. E.; Feng, Y.; Frerichs, H.; Loarte, A.; Pitts, R. A.; Reiser, D.; Fenstermacher, M. E.; Harting, D.; Kirschner, A.; Kukushkin, A.; Lunt, T.; Saibene, G.; Reiter, D.; Samm, U.; Wiesen, S.

2016-06-01

Results from three-dimensional modeling of plasma edge transport and plasma-wall interactions during application of resonant magnetic perturbation (RMP) fields for control of edge-localized modes in the ITER standard 15 MA Q  =  10 H-mode are presented. The full 3D plasma fluid and kinetic neutral transport code EMC3-EIRENE is used for the modeling. Four characteristic perturbed magnetic topologies are considered and discussed with reference to the axisymmetric case without RMP fields. Two perturbation field amplitudes at full and half of the ITER ELM control coil current capability using the vacuum approximation are compared to a case including a strongly screening plasma response. In addition, a vacuum field case at high q 95  =  4.2 featuring increased magnetic shear has been modeled. Formation of a three-dimensional plasma boundary is seen for all four perturbed magnetic topologies. The resonant field amplitudes and the effective radial magnetic field at the separatrix define the shape and extension of the 3D plasma boundary. Opening of the magnetic field lines from inside the separatrix establishes scrape-off layer-like channels of direct parallel particle and heat flux towards the divertor yielding a reduction of the main plasma thermal and particle confinement. This impact on confinement is most accentuated at full RMP current and is strongly reduced when screened RMP fields are considered, as well as for the reduced coil current cases. The divertor fluxes are redirected into a three-dimensional pattern of helical magnetic footprints on the divertor target tiles. At maximum perturbation strength, these fingers stretch out as far as 60 cm across the divertor targets, yielding heat flux spreading and the reduction of peak heat fluxes by 30%. However, at the same time substantial and highly localized heat fluxes reach divertor areas well outside of the axisymmetric heat flux decay profile. Reduced RMP amplitudes due to screening or reduced RMP coil current yield a reduction of the width of the divertor flux spreading to about 20-25 cm and cause increased peak heat fluxes back to values similar to those in the axisymmetric case. The dependencies of these features on the divertor recycling regime and the perpendicular transport assumptions, as well as toroidal averaged effects mimicking rotation of the RMP field, are discussed in the paper.

Two-dimensional analytic weighting functions for limb scattering

NASA Astrophysics Data System (ADS)

Zawada, D. J.; Bourassa, A. E.; Degenstein, D. A.

2017-10-01

Through the inversion of limb scatter measurements it is possible to obtain vertical profiles of trace species in the atmosphere. Many of these inversion methods require what is often referred to as weighting functions, or derivatives of the radiance with respect to concentrations of trace species in the atmosphere. Several radiative transfer models have implemented analytic methods to calculate weighting functions, alleviating the computational burden of traditional numerical perturbation methods. Here we describe the implementation of analytic two-dimensional weighting functions, where derivatives are calculated relative to atmospheric constituents in a two-dimensional grid of altitude and angle along the line of sight direction, in the SASKTRAN-HR radiative transfer model. Two-dimensional weighting functions are required for two-dimensional inversions of limb scatter measurements. Examples are presented where the analytic two-dimensional weighting functions are calculated with an underlying one-dimensional atmosphere. It is shown that the analytic weighting functions are more accurate than ones calculated with a single scatter approximation, and are orders of magnitude faster than a typical perturbation method. Evidence is presented that weighting functions for stratospheric aerosols calculated under a single scatter approximation may not be suitable for use in retrieval algorithms under solar backscatter conditions.

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 applications of the advective-diffusive codes for the inner magnetosphere

NASA Astrophysics Data System (ADS)

Aseev, N. A.; Shprits, Y. Y.; Drozdov, A. Y.; Kellerman, A. C.

2016-11-01

In this study we present analytical solutions for convection and diffusion equations. We gather here the analytical solutions for the one-dimensional convection equation, the two-dimensional convection problem, and the one- and two-dimensional diffusion equations. Using obtained analytical solutions, we test the four-dimensional Versatile Electron Radiation Belt code (the VERB-4D code), which solves the modified Fokker-Planck equation with additional convection terms. The ninth-order upwind numerical scheme for the one-dimensional convection equation shows much more accurate results than the results obtained with the third-order scheme. The universal limiter eliminates unphysical oscillations generated by high-order linear upwind schemes. Decrease in the space step leads to convergence of a numerical solution of the two-dimensional diffusion equation with mixed terms to the analytical solution. We compare the results of the third- and ninth-order schemes applied to magnetospheric convection modeling. The results show significant differences in electron fluxes near geostationary orbit when different numerical schemes are used.

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.

Three dimensional nonlinear simulations of edge localized modes on the EAST tokamak using BOUT++ code

NASA Astrophysics Data System (ADS)

Liu, Z. X.; Xu, X. Q.; Gao, X.; Xia, T. Y.; Joseph, I.; Meyer, W. H.; Liu, S. C.; Xu, G. S.; Shao, L. M.; Ding, S. Y.; Li, G. Q.; Li, J. G.

2014-09-01

Experimental measurements of edge localized modes (ELMs) observed on the EAST experiment are compared to linear and nonlinear theoretical simulations of peeling-ballooning modes using the BOUT++ code. Simulations predict that the dominant toroidal mode number of the ELM instability becomes larger for lower current, which is consistent with the mode structure captured with visible light using an optical CCD camera. The poloidal mode number of the simulated pressure perturbation shows good agreement with the filamentary structure observed by the camera. The nonlinear simulation is also consistent with the experimentally measured energy loss during an ELM crash and with the radial speed of ELM effluxes measured using a gas puffing imaging diagnostic.

A new feedback image encryption scheme based on perturbation with dynamical compound chaotic sequence cipher generator

NASA Astrophysics Data System (ADS)

Tong, Xiaojun; Cui, Minggen; Wang, Zhu

2009-07-01

The design of the new compound two-dimensional chaotic function is presented by exploiting two one-dimensional chaotic functions which switch randomly, and the design is used as a chaotic sequence generator which is proved by Devaney's definition proof of chaos. The properties of compound chaotic functions are also proved rigorously. In order to improve the robustness against difference cryptanalysis and produce avalanche effect, a new feedback image encryption scheme is proposed using the new compound chaos by selecting one of the two one-dimensional chaotic functions randomly and a new image pixels method of permutation and substitution is designed in detail by array row and column random controlling based on the compound chaos. The results from entropy analysis, difference analysis, statistical analysis, sequence randomness analysis, cipher sensitivity analysis depending on key and plaintext have proven that the compound chaotic sequence cipher can resist cryptanalytic, statistical and brute-force attacks, and especially it accelerates encryption speed, and achieves higher level of security. By the dynamical compound chaos and perturbation technology, the paper solves the problem of computer low precision of one-dimensional chaotic function.

Electronic field emission models beyond the Fowler-Nordheim one

NASA Astrophysics Data System (ADS)

Lepetit, Bruno

2017-12-01

We propose several quantum mechanical models to describe electronic field emission from first principles. These models allow us to correlate quantitatively the electronic emission current with the electrode surface details at the atomic scale. They all rely on electronic potential energy surfaces obtained from three dimensional density functional theory calculations. They differ by the various quantum mechanical methods (exact or perturbative, time dependent or time independent), which are used to describe tunneling through the electronic potential energy barrier. Comparison of these models between them and with the standard Fowler-Nordheim one in the context of one dimensional tunneling allows us to assess the impact on the accuracy of the computed current of the approximations made in each model. Among these methods, the time dependent perturbative one provides a well-balanced trade-off between accuracy and computational cost.

Application of advanced computational codes in the design of an experiment for a supersonic throughflow fan rotor

NASA Technical Reports Server (NTRS)

Wood, Jerry R.; Schmidt, James F.; Steinke, Ronald J.; Chima, Rodrick V.; Kunik, William G.

1987-01-01

Increased emphasis on sustained supersonic or hypersonic cruise has revived interest in the supersonic throughflow fan as a possible component in advanced propulsion systems. Use of a fan that can operate with a supersonic inlet axial Mach number is attractive from the standpoint of reducing the inlet losses incurred in diffusing the flow from a supersonic flight Mach number to a subsonic one at the fan face. The design of the experiment using advanced computational codes to calculate the components required is described. The rotor was designed using existing turbomachinery design and analysis codes modified to handle fully supersonic axial flow through the rotor. A two-dimensional axisymmetric throughflow design code plus a blade element code were used to generate fan rotor velocity diagrams and blade shapes. A quasi-three-dimensional, thin shear layer Navier-Stokes code was used to assess the performance of the fan rotor blade shapes. The final design was stacked and checked for three-dimensional effects using a three-dimensional Euler code interactively coupled with a two-dimensional boundary layer code. The nozzle design in the expansion region was analyzed with a three-dimensional parabolized viscous code which corroborated the results from the Euler code. A translating supersonic diffuser was designed using these same codes.

3D-PDR: Three-dimensional photodissociation region code

NASA Astrophysics Data System (ADS)

Bisbas, T. G.; Bell, T. A.; Viti, S.; Yates, J.; Barlow, M. J.

2018-03-01

3D-PDR is a three-dimensional photodissociation region code written in Fortran. It uses the Sundials package (written in C) to solve the set of ordinary differential equations and it is the successor of the one-dimensional PDR code UCL_PDR (ascl:1303.004). Using the HEALpix ray-tracing scheme (ascl:1107.018), 3D-PDR solves a three-dimensional escape probability routine and evaluates the attenuation of the far-ultraviolet radiation in the PDR and the propagation of FIR/submm emission lines out of the PDR. The code is parallelized (OpenMP) and can be applied to 1D and 3D problems.

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.

A generalized one-dimensional computer code for turbomachinery cooling passage flow calculations

NASA Technical Reports Server (NTRS)

Kumar, Ganesh N.; Roelke, Richard J.; Meitner, Peter L.

1989-01-01

A generalized one-dimensional computer code for analyzing the flow and heat transfer in the turbomachinery cooling passages was developed. This code is capable of handling rotating cooling passages with turbulators, 180 degree turns, pin fins, finned passages, by-pass flows, tip cap impingement flows, and flow branching. The code is an extension of a one-dimensional code developed by P. Meitner. In the subject code, correlations for both heat transfer coefficient and pressure loss computations were developed to model each of the above mentioned type of coolant passages. The code has the capability of independently computing the friction factor and heat transfer coefficient on each side of a rectangular passage. Either the mass flow at the inlet to the channel or the exit plane pressure can be specified. For a specified inlet total temperature, inlet total pressure, and exit static pressure, the code computers the flow rates through the main branch and the subbranches, flow through tip cap for impingement cooling, in addition to computing the coolant pressure, temperature, and heat transfer coefficient distribution in each coolant flow branch. Predictions from the subject code for both nonrotating and rotating passages agree well with experimental data. The code was used to analyze the cooling passage of a research cooled radial rotor.

Restoration of dimensional reduction in the random-field Ising model at five dimensions

NASA Astrophysics Data System (ADS)

Fytas, Nikolaos G.; Martín-Mayor, Víctor; Picco, Marco; Sourlas, Nicolas

2017-04-01

The random-field Ising model is one of the few disordered systems where the perturbative renormalization group can be carried out to all orders of perturbation theory. This analysis predicts dimensional reduction, i.e., that the critical properties of the random-field Ising model in D dimensions are identical to those of the pure Ising ferromagnet in D -2 dimensions. It is well known that dimensional reduction is not true in three dimensions, thus invalidating the perturbative renormalization group prediction. Here, we report high-precision numerical simulations of the 5D random-field Ising model at zero temperature. We illustrate universality by comparing different probability distributions for the random fields. We compute all the relevant critical exponents (including the critical slowing down exponent for the ground-state finding algorithm), as well as several other renormalization-group invariants. The estimated values of the critical exponents of the 5D random-field Ising model are statistically compatible to those of the pure 3D Ising ferromagnet. These results support the restoration of dimensional reduction at D =5 . We thus conclude that the failure of the perturbative renormalization group is a low-dimensional phenomenon. We close our contribution by comparing universal quantities for the random-field problem at dimensions 3 ≤D <6 to their values in the pure Ising model at D -2 dimensions, and we provide a clear verification of the Rushbrooke equality at all studied dimensions.

Restoration of dimensional reduction in the random-field Ising model at five dimensions.

PubMed

Fytas, Nikolaos G; Martín-Mayor, Víctor; Picco, Marco; Sourlas, Nicolas

2017-04-01

The random-field Ising model is one of the few disordered systems where the perturbative renormalization group can be carried out to all orders of perturbation theory. This analysis predicts dimensional reduction, i.e., that the critical properties of the random-field Ising model in D dimensions are identical to those of the pure Ising ferromagnet in D-2 dimensions. It is well known that dimensional reduction is not true in three dimensions, thus invalidating the perturbative renormalization group prediction. Here, we report high-precision numerical simulations of the 5D random-field Ising model at zero temperature. We illustrate universality by comparing different probability distributions for the random fields. We compute all the relevant critical exponents (including the critical slowing down exponent for the ground-state finding algorithm), as well as several other renormalization-group invariants. The estimated values of the critical exponents of the 5D random-field Ising model are statistically compatible to those of the pure 3D Ising ferromagnet. These results support the restoration of dimensional reduction at D=5. We thus conclude that the failure of the perturbative renormalization group is a low-dimensional phenomenon. We close our contribution by comparing universal quantities for the random-field problem at dimensions 3≤D<6 to their values in the pure Ising model at D-2 dimensions, and we provide a clear verification of the Rushbrooke equality at all studied dimensions.

Instability of the sliding Luttinger liquid

NASA Astrophysics Data System (ADS)

Fleurov, V.; Kagalovsky, V.; Lerner, I. V.; Yurkevich, I. V.

2018-05-01

We revise a phase diagram for the sliding Luttinger liquid (SLL) of coupled one-dimensional quantum wires packed in two- or three-dimensional arrays in the absence of a magnetic field. We analyse whether physically justifiable (reasonable) inter-wire interactions, i.e. either the screened Coulomb or ‘Coulomb-blockade’ type interactions, stabilise the SLL phase. Calculating the scaling dimensions of the most relevant perturbations (the inter-wire single-particle hybridisation, charge-density wave, and superconducting inter-wire couplings), we find that their combination always destroys the SLL phase for the repulsive intra-wire interaction. However, suppressing the inter-wire tunnelling of repulsive fermions (when the charge-density wave is the only remaining perturbation), one can observe a stability region emerging due to the inter-wire forward scattering interaction.

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.

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

Steiner, J.L.; Lime, J.F.; Elson, J.S.

One dimensional TRAC transient calculations of the process inherent ultimate safety (PIUS) advanced reactor design were performed for a pump-trip SCRAM. The TRAC calculations showed that the reactor power response and shutdown were in qualitative agreement with the one-dimensional analyses presented in the PIUS Preliminary Safety Information Document (PSID) submitted by Asea Brown Boveri (ABB) to the US Nuclear Regulatory Commission for preapplication safety review. The PSID analyses were performed with the ABB-developed RIGEL code. The TRAC-calculated phenomena and trends were also similar to those calculated with another one-dimensional PIUS model, the Brookhaven National Laboratory developed PIPA code. A TRACmore » pump-trip SCRAM transient has also been calculated with a TRAC model containing a multi-dimensional representation of the PIUS intemal flow structures and core region. The results obtained using the TRAC fully one-dimensional PIUS model are compared to the RIGEL, PIPA, and TRAC multi-dimensional results.« less

Linearized compressible-flow theory for sonic flight speeds

NASA Technical Reports Server (NTRS)

Heaslet, Max A; Lomax, Harvard; Spreiter, John R

1950-01-01

The partial differential equation for the perturbation velocity potential is examined for free-stream Mach numbers close to and equal to one. It is found that, under the assumptions of linearized theory, solutions can be found consistent with the theory for lifting-surface problems both in stationary three-dimensional flow and in unsteady two-dimensional flow. Several examples are solved including a three dimensional swept-back wing and two dimensional harmonically-oscillating wing, both for a free stream Mach number equal to one. Momentum relations for the evaluation of wave and vortex drag are also discussed. (author)

Solitonic and chaotic behaviors of a (3+1)-dimensional nonlinear Schrödinger equation in a magnetized electron-positron-ion plasma

NASA Astrophysics Data System (ADS)

Zhen, Hui-Ling; Tian, Bo; Xie, Xi-Yang; Chai, Jun

2015-05-01

In a magnetized e-p-i plasma with two-electron temperatures, the (3+1)-dimensional nonlinear Schrödinger equation for the ion-acoustic envelope solitons is hereby investigated. Bright and dark soliton solutions are both obtained. It is found that the soliton amplitude is inversely related to P, Q and R, with P, Q and R respectively as the coefficients of the dispersive term, nonlinear term and combined effect of the transverse perturbation and magnetic field. Head-on interactions are displayed, and with Q and R decreasing, can be transferred into the overtaking ones. Bright bound-state solitons are obtained, and the interaction period decreases with Q increasing. Both the head-on and overtaking interactions between the two dark solitons are displayed, and such a head-on interaction can be transformed into the overtaking one with R increasing. Upon the introduction of the external perturbation, the developed and weak chaos are observed. Phase projections and power spectra are presented. Difference between the two chaotic motions roots in the inequality between the nonlinear and perturbed terms. Developed chaos can be weakened with Q decreasing, or with the frequency of external perturbation increasing.

Parallel processing a three-dimensional free-lagrange code

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

Mandell, D.A.; Trease, H.E.

1989-01-01

A three-dimensional, time-dependent free-Lagrange hydrodynamics code has been multitasked and autotasked on a CRAY X-MP/416. The multitasking was done by using the Los Alamos Multitasking Control Library, which is a superset of the CRAY multitasking library. Autotasking is done by using constructs which are only comment cards if the source code is not run through a preprocessor. The three-dimensional algorithm has presented a number of problems that simpler algorithms, such as those for one-dimensional hydrodynamics, did not exhibit. Problems in converting the serial code, originally written for a CRAY-1, to a multitasking code are discussed. Autotasking of a rewritten versionmore » of the code is discussed. Timing results for subroutines and hot spots in the serial code are presented and suggestions for additional tools and debugging aids are given. Theoretical speedup results obtained from Amdahl's law and actual speedup results obtained on a dedicated machine are presented. Suggestions for designing large parallel codes are given.« less

Reducing the two-loop large-scale structure power spectrum to low-dimensional, radial integrals

DOE PAGES

Schmittfull, Marcel; Vlah, Zvonimir

2016-11-28

Modeling the large-scale structure of the universe on nonlinear scales has the potential to substantially increase the science return of upcoming surveys by increasing the number of modes available for model comparisons. One way to achieve this is to model nonlinear scales perturbatively. Unfortunately, this involves high-dimensional loop integrals that are cumbersome to evaluate. Here, trying to simplify this, we show how two-loop (next-to-next-to-leading order) corrections to the density power spectrum can be reduced to low-dimensional, radial integrals. Many of those can be evaluated with a one-dimensional fast Fourier transform, which is significantly faster than the five-dimensional Monte-Carlo integrals thatmore » are needed otherwise. The general idea of this fast fourier transform perturbation theory method is to switch between Fourier and position space to avoid convolutions and integrate over orientations, leaving only radial integrals. This reformulation is independent of the underlying shape of the initial linear density power spectrum and should easily accommodate features such as those from baryonic acoustic oscillations. We also discuss how to account for halo bias and redshift space distortions.« less

Three-dimensional drift kinetic response of high- β plasmas in the DIII-D tokamak

DOE PAGES

Wang, Zhirui R.; Lanctot, Matthew J.; Liu, Y. Q.; ...

2015-04-07

A quantitative interpretation of the experimentally measured high pressure plasma response to externally applied three-dimensional (3D) magnetic field perturbations, across the no-wall Troyon limit, is achieved. The key to success is the self-consistent inclusion of the drift kinetic resonance effects in numerical modeling using the MARS-K code. This resolves an outstanding issue of ideal magneto-hydrodynamic model, which signi cantly over-predicts the plasma induced field ampli fication near the no-wall limit, as compared to experiments. The self-consistent drift kinetic model leads to quantitative agreement not only for the measured 3D field amplitude and toroidal phase, but also for the measured internalmore » 3D displacement of the plasma.« less

User's manual for three dimensional FDTD version B code for scattering from frequency-dependent dielectric materials

NASA Technical Reports Server (NTRS)

Beggs, John H.; Luebbers, Raymond J.; Kunz, Karl S.

1992-01-01

The Penn State Finite Difference Time Domain Electromagnetic Code Version B is a three dimensional numerical electromagnetic scattering code based upon the Finite Difference Time Domain Technique (FDTD). The supplied version of the code is one version of our current three dimensional FDTD code set. This manual provides a description of the code and corresponding results for several scattering problems. The manual is organized into 14 sections: introduction, description of the FDTD method, operation, resource requirements, Version B code capabilities, a brief description of the default scattering geometry, a brief description of each subroutine, a description of the include file, a discussion of radar cross section computations, a discussion of some scattering results, a sample problem setup section, a new problem checklist, references and figure titles.

Chemically Reacting One-Dimensional Gas-Particle Flows

NASA Technical Reports Server (NTRS)

Tevepaugh, J. A.; Penny, M. M.

1975-01-01

The governing equations for the one-dimensional flow of a gas-particle system are discussed. Gas-particle effects are coupled via the system momentum and energy equations with the gas assumed to be chemically frozen or in chemical equilibrium. A computer code for calculating the one-dimensional flow of a gas-particle system is discussed and a user's input guide presented. The computer code provides for the expansion of the gas-particle system from a specified starting velocity and nozzle inlet geometry. Though general in nature, the final output of the code is a startline for initiating the solution of a supersonic gas-particle system in rocket nozzles. The startline includes gasdynamic data defining gaseous startline points from the nozzle centerline to the nozzle wall and particle properties at points along the gaseous startline.

2d affine XY-spin model/4d gauge theory duality and deconfinement

NASA Astrophysics Data System (ADS)

Anber, Mohamed M.; Poppitz, Erich; Ünsal, Mithat

2012-04-01

We introduce a duality between two-dimensional XY-spin models with symmetry-breaking perturbations and certain four-dimensional SU(2) and SU(2)/ {{Z}_2} gauge theories, compactified on a small spatial circle {{R}^{{^{{{1},{2}}}}}} × {{S}^{{^{{1}}}}} , and considered at temperatures near the deconfinement transition. In a Euclidean set up, the theory is defined on {{R}^{{^{{2}}}}} × {{T}^{{^{{2}}}}} . Similarly, thermal gauge theories of higher rank are dual to new families of "affine" XY-spin models with perturbations. For rank two, these are related to models used to describe the melting of a 2d crystal with a triangular lattice. The connection is made through a multi-component electric-magnetic Coulomb gas representation for both systems. Perturbations in the spin system map to topological defects in the gauge theory, such as monopole-instantons or magnetic bions, and the vortices in the spin system map to the electrically charged W-bosons in field theory (or vice versa, depending on the duality frame). The duality permits one to use the two-dimensional technology of spin systems to study the thermal deconfinement and discrete chiral transitions in four-dimensional SU( N c ) gauge theories with n f ≥1 adjoint Weyl fermions.

Accelerating execution of the integrated TIGER series Monte Carlo radiation transport codes

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

Smith, L.M.; Hochstedler, R.D.

1997-02-01

Execution of the integrated TIGER series (ITS) of coupled electron/photon Monte Carlo radiation transport codes has been accelerated by modifying the FORTRAN source code for more efficient computation. Each member code of ITS was benchmarked and profiled with a specific test case that directed the acceleration effort toward the most computationally intensive subroutines. Techniques for accelerating these subroutines included replacing linear search algorithms with binary versions, replacing the pseudo-random number generator, reducing program memory allocation, and proofing the input files for geometrical redundancies. All techniques produced identical or statistically similar results to the original code. Final benchmark timing of themore » accelerated code resulted in speed-up factors of 2.00 for TIGER (the one-dimensional slab geometry code), 1.74 for CYLTRAN (the two-dimensional cylindrical geometry code), and 1.90 for ACCEPT (the arbitrary three-dimensional geometry code).« less

Nonambipolar Transport and Torque in Perturbed Equilibria

NASA Astrophysics Data System (ADS)

Logan, N. C.; Park, J.-K.; Wang, Z. R.; Berkery, J. W.; Kim, K.; Menard, J. E.

2013-10-01

A new Perturbed Equilibrium Nonambipolar Transport (PENT) code has been developed to calculate the neoclassical toroidal torque from radial current composed of both passing and trapped particles in perturbed equilibria. This presentation outlines the physics approach used in the development of the PENT code, with emphasis on the effects of retaining general aspect-ratio geometric effects. First, nonambipolar transport coefficients and corresponding neoclassical toroidal viscous (NTV) torque in perturbed equilibria are re-derived from the first order gyro-drift-kinetic equation in the ``combined-NTV'' PENT formalism. The equivalence of NTV torque and change in potential energy due to kinetic effects [J-K. Park, Phys. Plas., 2011] is then used to showcase computational challenges shared between PENT and stability codes MISK and MARS-K. Extensive comparisons to a reduced model, which makes numerous large aspect ratio approximations, are used throughout to emphasize geometry dependent physics such as pitch angle resonances. These applications make extensive use of the PENT code's native interfacing with the Ideal Perturbed Equilibrium Code (IPEC), and the combination of these codes is a key step towards an iterative solver for self-consistent perturbed equilibrium torque. Supported by US DOE contract #DE-AC02-09CH11466 and the DOE Office of Science Graduate Fellowship administered by the Oak Ridge Institute for Science & Education under contract #DE-AC05-06OR23100.

Self-similarity of a Rayleigh–Taylor mixing layer at low Atwood number with a multimode initial perturbation

DOE PAGES

Morgan, B. E.; Olson, B. J.; White, J. E.; ...

2017-06-29

High-fidelity large eddy simulation (LES) of a low-Atwood number (A = 0.05) Rayleigh-Taylor mixing layer is performed using the tenth-order compact difference code Miranda. An initial multimode perturbation spectrum is specified in Fourier space as a function of mesh resolution such that a database of results is obtained in which each successive level of increased grid resolution corresponds approximately to one additional doubling of the mixing layer width, or generation. The database is then analyzed to determine approximate requirements for self-similarity, and a new metric is proposed to quantify how far a given simulation is from the limit of self-similarity.more » It is determined that mixing layer growth reaches a high degree of self-similarity after approximately 4.5 generations. Statistical convergence errors and boundary effects at late time, however, make it impossible to draw similar conclusions regarding the self-similar growth of more sensitive turbulence parameters. Finally, self-similar turbulence profiles from the LES database are compared with one-dimensional simulations using the k-L-a and BHR-2 Reynolds-averaged Navier-Stokes (RANS) models. The k-L-a model, which is calibrated to reproduce a quadratic turbulence kinetic energy profile for a self-similar mixing layer, is found to be in better agreement with the LES than BHR-2 results.« less

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

Floating phase in the one-dimensional transverse axial next-nearest-neighbor Ising model.

PubMed

Chandra, Anjan Kumar; Dasgupta, Subinay

2007-02-01

To study the ground state of an axial next-nearest-neighbor Ising chain under transverse field as a function of frustration parameter kappa and field strength Gamma, we present here two different perturbative analyses. In one, we consider the (known) ground state at kappa=0.5 and Gamma=0 as the unperturbed state and treat an increase of the field from 0 to Gamma coupled with an increase of kappa from 0.5 to 0.5+rGamma/J as perturbation. The first-order perturbation correction to eigenvalue can be calculated exactly and we could conclude that there are only two phase-transition lines emanating from the point kappa=0.5, Gamma=0. In the second perturbation scheme, we consider the number of domains of length 1 as the perturbation and obtain the zeroth-order eigenfunction for the perturbed ground state. From the longitudinal spin-spin correlation, we conclude that floating phase exists for small values of transverse field over the entire region intermediate between the ferromagnetic phase and antiphase.

User's manual for the one-dimensional hypersonic experimental aero-thermodynamic (1DHEAT) data reduction code

NASA Technical Reports Server (NTRS)

Hollis, Brian R.

1995-01-01

A FORTRAN computer code for the reduction and analysis of experimental heat transfer data has been developed. This code can be utilized to determine heat transfer rates from surface temperature measurements made using either thin-film resistance gages or coaxial surface thermocouples. Both an analytical and a numerical finite-volume heat transfer model are implemented in this code. The analytical solution is based on a one-dimensional, semi-infinite wall thickness model with the approximation of constant substrate thermal properties, which is empirically corrected for the effects of variable thermal properties. The finite-volume solution is based on a one-dimensional, implicit discretization. The finite-volume model directly incorporates the effects of variable substrate thermal properties and does not require the semi-finite wall thickness approximation used in the analytical model. This model also includes the option of a multiple-layer substrate. Fast, accurate results can be obtained using either method. This code has been used to reduce several sets of aerodynamic heating data, of which samples are included in this report.

User's manual for three dimensional FDTD version C code for scattering from frequency-independent dielectric and magnetic materials

NASA Technical Reports Server (NTRS)

Beggs, John H.; Luebbers, Raymond J.; Kunz, Karl S.

1991-01-01

The Penn State Finite Difference Time Domain Electromagnetic Scattering Code Version C is a three dimensional numerical electromagnetic scattering code based upon the Finite Difference Time Domain Technique (FDTD). The supplied version of the code is one version of our current three dimensional FDTD code set. This manual provides a description of the code and corresponding results for several scattering problems. The manual is organized into fourteen sections: introduction, description of the FDTD method, operation, resource requirements, Version C code capabilities, a brief description of the default scattering geometry, a brief description of each subroutine, a description of the include file (COMMONC.FOR), a section briefly discussing Radar Cross Section (RCS) computations, a section discussing some scattering results, a sample problem setup section, a new problem checklist, references and figure titles.

User's manual for three dimensional FDTD version D code for scattering from frequency-dependent dielectric and magnetic materials

NASA Technical Reports Server (NTRS)

Beggs, John H.; Luebbers, Raymond J.; Kunz, Karl S.

1991-01-01

The Penn State Finite Difference Time Domain Electromagnetic Scattering Code Version D is a three dimensional numerical electromagnetic scattering code based upon the Finite Difference Time Domain Technique (FDTD). The supplied version of the code is one version of our current three dimensional FDTD code set. This manual provides a description of the code and corresponding results for several scattering problems. The manual is organized into fourteen sections: introduction, description of the FDTD method, operation, resource requirements, Version D code capabilities, a brief description of the default scattering geometry, a brief description of each subroutine, a description of the include file (COMMOND.FOR), a section briefly discussing Radar Cross Section (RCS) computations, a section discussing some scattering results, a sample problem setup section, a new problem checklist, references and figure titles.

User's manual for three dimensional FDTD version A code for scattering from frequency-independent dielectric materials

NASA Technical Reports Server (NTRS)

Beggs, John H.; Luebbers, Raymond J.; Kunz, Karl S.

1992-01-01

The Penn State Finite Difference Time Domain (FDTD) Electromagnetic Scattering Code Version A is a three dimensional numerical electromagnetic scattering code based on the Finite Difference Time Domain technique. The supplied version of the code is one version of our current three dimensional FDTD code set. The manual provides a description of the code and the corresponding results for the default scattering problem. The manual is organized into 14 sections: introduction, description of the FDTD method, operation, resource requirements, Version A code capabilities, a brief description of the default scattering geometry, a brief description of each subroutine, a description of the include file (COMMONA.FOR), a section briefly discussing radar cross section (RCS) computations, a section discussing the scattering results, a sample problem setup section, a new problem checklist, references, and figure titles.

User's manual for three dimensional FDTD version C code for scattering from frequency-independent dielectric and magnetic materials

NASA Technical Reports Server (NTRS)

Beggs, John H.; Luebbers, Raymond J.; Kunz, Karl S.

1992-01-01

The Penn State Finite Difference Time Domain Electromagnetic Scattering Code Version C is a three-dimensional numerical electromagnetic scattering code based on the Finite Difference Time Domain (FDTD) technique. The supplied version of the code is one version of our current three-dimensional FDTD code set. The manual given here provides a description of the code and corresponding results for several scattering problems. The manual is organized into 14 sections: introduction, description of the FDTD method, operation, resource requirements, Version C code capabilities, a brief description of the default scattering geometry, a brief description of each subroutine, a description of the include file (COMMONC.FOR), a section briefly discussing radar cross section computations, a section discussing some scattering results, a new problem checklist, references, and figure titles.

User's manual for three dimensional FDTD version B code for scattering from frequency-dependent dielectric materials

NASA Technical Reports Server (NTRS)

Beggs, John H.; Luebbers, Raymond J.; Kunz, Karl S.

1991-01-01

The Penn State Finite Difference Time Domain Electromagnetic Scattering Code Version B is a three dimensional numerical electromagnetic scattering code based upon the Finite Difference Time Domain Technique (FDTD). The supplied version of the code is one version of our current three dimensional FDTD code set. This manual provides a description of the code and corresponding results for several scattering problems. The manual is organized into fourteen sections: introduction, description of the FDTD method, operation, resource requirements, Version B code capabilities, a brief description of the default scattering geometry, a brief description of each subroutine, a description of the include file (COMMONB.FOR), a section briefly discussing Radar Cross Section (RCS) computations, a section discussing some scattering results, a sample problem setup section, a new problem checklist, references and figure titles.

FAST-PT: a novel algorithm to calculate convolution integrals in cosmological perturbation theory

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

McEwen, Joseph E.; Fang, Xiao; Hirata, Christopher M.

2016-09-01

We present a novel algorithm, FAST-PT, for performing convolution or mode-coupling integrals that appear in nonlinear cosmological perturbation theory. The algorithm uses several properties of gravitational structure formation—the locality of the dark matter equations and the scale invariance of the problem—as well as Fast Fourier Transforms to describe the input power spectrum as a superposition of power laws. This yields extremely fast performance, enabling mode-coupling integral computations fast enough to embed in Monte Carlo Markov Chain parameter estimation. We describe the algorithm and demonstrate its application to calculating nonlinear corrections to the matter power spectrum, including one-loop standard perturbation theorymore » and the renormalization group approach. We also describe our public code (in Python) to implement this algorithm. The code, along with a user manual and example implementations, is available at https://github.com/JoeMcEwen/FAST-PT.« less

Nonlinear evolution of energetic-particles-driven waves in collisionless plasmas

NASA Astrophysics Data System (ADS)

Li, Shuhan; Liu, Jinyuan; Wang, Feng; Shen, Wei; Li, Dong

2018-06-01

A one-dimensional electrostatic collisionless particle-in-cell code has been developed to study the nonlinear interaction between electrostatic waves and energetic particles (EPs). For a single wave, the results are clear and agree well with the existing theories. For coexisting two waves, although the mode nonlinear coupling between two wave fields is ignored, the second-order phase space islands can still exist between first-order islands generated by the two waves. However, the second-order phase islands are not formed by the superposed wave fields and the perturbed motions of EPs induced by the combined effect of two main resonances make these structures in phase space. Owing to these second-order islands, energy can be transferred between waves, even if the overlap of two main resonances never occurs. Depending on the distance between the main resonance islands in velocity space, the second-order island can affect the nonlinear dynamics and saturations of waves.

One-loop tests of supersymmetric gauge theories on spheres

DOE PAGES

Minahan, Joseph A.; Naseer, Usman

2017-07-14

Here, we show that a recently conjectured form for perturbative supersymmetric partition functions on spheres of general dimension d is consistent with the at space limit of 6-dimensional N = 1 super Yang-Mills. We also show that the partition functions for N = 1 8- and 9-dimensional theories are consistent with their known at space limits.

An Integrated Magnetic Circuit Model and Finite Element Model Approach to Magnetic Bearing Design

NASA Technical Reports Server (NTRS)

Provenza, Andrew J.; Kenny, Andrew; Palazzolo, Alan B.

2003-01-01

A code for designing magnetic bearings is described. The code generates curves from magnetic circuit equations relating important bearing performance parameters. Bearing parameters selected from the curves by a designer to meet the requirements of a particular application are input directly by the code into a three-dimensional finite element analysis preprocessor. This means that a three-dimensional computer model of the bearing being developed is immediately available for viewing. The finite element model solution can be used to show areas of magnetic saturation and make more accurate predictions of the bearing load capacity, current stiffness, position stiffness, and inductance than the magnetic circuit equations did at the start of the design process. In summary, the code combines one-dimensional and three-dimensional modeling methods for designing magnetic bearings.

Validation of the NCC Code for Staged Transverse Injection and Computations for a RBCC Combustor

NASA Technical Reports Server (NTRS)

Ajmani, Kumud; Liu, Nan-Suey

2005-01-01

The NCC code was validated for a case involving staged transverse injection into Mach 2 flow behind a rearward facing step. Comparisons with experimental data and with solutions from the FPVortex code was then used to perform computations to study fuel-air mixing for the combustor of a candidate rocket based combined cycle engine geometry. Comparisons with a one-dimensional analysis and a three-dimensional code (VULCAN) were performed to assess the qualitative and quantitative performance of the NCC solver.

Numerical investigations in three-dimensional internal flows

NASA Technical Reports Server (NTRS)

Rose, William C.

1991-01-01

In previous efforts, a two-dimensional full Navier-Stokes (FNS) code (SCRAM2D) was used in a design process that involved parametric modifications of the inlet geometry to arrive at what appeared to be an optimum inlet flowfield that produced a uniform flow at the exit in a very short distance. In these previous studies, the technologies for determining the contours with a 'man-in-the-loop' approach for both the ramp and cowl of the inlet were demonstrated, and nearly shock-free exiting flowfields were shown to be obtainable. The resulting two-dimensional compression contours were then used with swept sidewalls to form a three-dimensional inlet. Then the three-dimensional Navier-Stokes code (SCRAM3D) was used to investigate the inlet's three-dimensional flow. One of the major difficulties encountered in the previous studies was that associated with the relatively long time required to obtain a solution using even the 2D FNS code in the design process. Since one of the goals of high-speed inlet design is to produce inputs to the overall aircraft design in a timely manner, it was proposed for this year's research to examine 2D and 3D viscous flow solver techniques alternative to the NFS codes used to date. Areas of the inlet particularly identified for code speed up are those associated with the forebody and external flow ramp systems of the inlet. In these areas, parabolized, or space-marched, Navier-Stokes codes were proposed to be investigated for their applicability in the design process developed previously. This report describes the results of an investigation into the use of two other codes for analyzing the forebody and inlet ramp systems of high-speed inlets.

Calculation of Multistage Turbomachinery Using Steady Characteristic Boundary Conditions

NASA Technical Reports Server (NTRS)

Chima, Rodrick V.

1998-01-01

A multiblock Navier-Stokes analysis code for turbomachinery has been modified to allow analysis of multistage turbomachines. A steady averaging-plane approach was used to pass information between blade rows. Characteristic boundary conditions written in terms of perturbations about the mean flow from the neighboring blade row were used to allow close spacing between the blade rows without forcing the flow to be axisymmetric. In this report the multiblock code is described briefly and the characteristic boundary conditions and the averaging-plane implementation are described in detail. Two approaches for averaging the flow properties are also described. A two-dimensional turbine stator case was used to compare the characteristic boundary conditions with standard axisymmetric boundary conditions. Differences were apparent but small in this low-speed case. The two-stage fuel turbine used on the space shuttle main engines was then analyzed using a three-dimensional averaging-plane approach. Computed surface pressure distributions on the stator blades and endwalls and computed distributions of blade surface heat transfer coefficient on three blades showed very good agreement with experimental data from two tests.

Hydrodynamic instability growth of three-dimensional modulations in radiation-driven implosions with “low-foot” and “high-foot” drives at the National Ignition Facility

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

Smalyuk, V. A.; Weber, C. R.; Robey, H. F.

Hydrodynamic instability growth has been studied using three-dimensional (3-D) broadband modulations by comparing “high-foot” and “low-foot” spherical plastic (CH) capsule implosions at the National Ignition Facility (NIF). The initial perturbations included capsule outer-surface roughness and capsule-mounting membranes (“tents”) that were similar to those used in a majority of implosions on NIF. The tents with thicknesses of 31-nm, 46-nm, and 109-nm were used in the experiments. The outer-surface roughness in the “low-foot” experiment was similar to the standard specification, while it was increased by ~4 times in the “high-foot” experiment to compensate for the reduced growth. The ablation-front instability growth wasmore » measured using a Hydrodynamic Growth Radiography platform at a convergence ratio of 3. The dominant capsule perturbations, generated by the tent mountings, had measured perturbation amplitudes comparable to the capsule thickness with the “low-foot” drive. These tent perturbations were reduced by ~3 to 10 times in implosions with the “high-foot” drive. Unexpectedly, the measured perturbations with initially thinner tents were either larger or similar to the measured perturbations with thicker tents for both “high-foot” and “low-foot” drives. While the measured instability growth of 3-D broadband perturbations was also significantly reduced by ~5 to 10 times with the “high-foot” drive, compared to the “low-foot” drive, the growth mitigation was stronger than expected based on previous “growth-factor” results measured with two-dimensional modulations. Lastly, one of the hypotheses to explain the results is based on the 3-D modulations of the oxygen content in the bulk of the capsule having a stronger effect on the overall growth of capsule perturbations than the outer-surface capsule roughness.« less

Hydrodynamic instability growth of three-dimensional modulations in radiation-driven implosions with “low-foot” and “high-foot” drives at the National Ignition Facility

DOE PAGES

Smalyuk, V. A.; Weber, C. R.; Robey, H. F.; ...

2017-04-11

Hydrodynamic instability growth has been studied using three-dimensional (3-D) broadband modulations by comparing “high-foot” and “low-foot” spherical plastic (CH) capsule implosions at the National Ignition Facility (NIF). The initial perturbations included capsule outer-surface roughness and capsule-mounting membranes (“tents”) that were similar to those used in a majority of implosions on NIF. The tents with thicknesses of 31-nm, 46-nm, and 109-nm were used in the experiments. The outer-surface roughness in the “low-foot” experiment was similar to the standard specification, while it was increased by ~4 times in the “high-foot” experiment to compensate for the reduced growth. The ablation-front instability growth wasmore » measured using a Hydrodynamic Growth Radiography platform at a convergence ratio of 3. The dominant capsule perturbations, generated by the tent mountings, had measured perturbation amplitudes comparable to the capsule thickness with the “low-foot” drive. These tent perturbations were reduced by ~3 to 10 times in implosions with the “high-foot” drive. Unexpectedly, the measured perturbations with initially thinner tents were either larger or similar to the measured perturbations with thicker tents for both “high-foot” and “low-foot” drives. While the measured instability growth of 3-D broadband perturbations was also significantly reduced by ~5 to 10 times with the “high-foot” drive, compared to the “low-foot” drive, the growth mitigation was stronger than expected based on previous “growth-factor” results measured with two-dimensional modulations. Lastly, one of the hypotheses to explain the results is based on the 3-D modulations of the oxygen content in the bulk of the capsule having a stronger effect on the overall growth of capsule perturbations than the outer-surface capsule roughness.« less

A User's Guide to AMR1D: An Instructional Adaptive Mesh Refinement Code for Unstructured Grids

NASA Technical Reports Server (NTRS)

deFainchtein, Rosalinda

1996-01-01

This report documents the code AMR1D, which is currently posted on the World Wide Web (http://sdcd.gsfc.nasa.gov/ESS/exchange/contrib/de-fainchtein/adaptive _mesh_refinement.html). AMR1D is a one-dimensional finite element fluid-dynamics solver, capable of adaptive mesh refinement (AMR). It was written as an instructional tool for AMR on unstructured mesh codes. It is meant to illustrate the minimum requirements for AMR on more than one dimension. For that purpose, it uses the same type of data structure that would be necessary on a two-dimensional AMR code (loosely following the algorithm described by Lohner).

Vibroacoustic response of structures and perturbation Reynolds stress near structure-turbulence interface

NASA Technical Reports Server (NTRS)

Maekawa, S.; Lin, Y. K.

1977-01-01

The interaction between a turbulent flow and certain types of structures which respond to its excitation is investigated. One-dimensional models were used to develop the basic ideas applied to a second model resembling the fuselage construction of an aircraft. In the two-dimensional case a simple membrane, with a small random variation in the membrane tension, was used. A decaying turbulence was constructed by superposing infinitely many components, each of which is convected as a frozen pattern at a different velocity. Structure-turbulence interaction results are presented in terms of the spectral densities of the structural response and the perturbation Reynolds stress in the fluid at the vicinity of the interface.

One-Dimensional Czedli-Type Islands

ERIC Educational Resources Information Center

Horvath, Eszter K.; Mader, Attila; Tepavcevic, Andreja

2011-01-01

The notion of an island has surfaced in recent algebra and coding theory research. Discrete versions provide interesting combinatorial problems. This paper presents the one-dimensional case with finitely many heights, a topic convenient for student research.

Numerical simulations of the convective flame in white dwarfs

NASA Technical Reports Server (NTRS)

Livne, Eli

1993-01-01

A first step toward better understanding of the mechanism driving convective flames in exploding white dwarfs is presented. The propagation of the convective flame is examined using a two-dimensional implicit hydrodynamical code. The large scales of the instability are captured by the grid while the scales that are smaller than the grid resolution are approximated by a mixing-length approximation. It is found that largescale perturbations (of order of the pressure scale height) do grow significantly during the expansion, leading to a very nonspherical burning front. The combustion rate is strongly enhanced (compared to the unperturbed case) during the first second, but later the expansion of the star suppresses the flame speed, leading to only partial incineration of the nuclear fuel. Our results imply that large-scale perturbations by themselves are not enough to explain the mechanism by which convective flames are driven, and a study of the whole spectrum of relevant perturbations is needed. The implications of these preliminary results on future simulations, in the context of current models for Type Ia supernovae, are discussed.

Solution of the Burnett equations for hypersonic flows near the continuum limit

NASA Technical Reports Server (NTRS)

Imlay, Scott T.

1992-01-01

The INCA code, a three-dimensional Navier-Stokes code for analysis of hypersonic flowfields, was modified to analyze the lower reaches of the continuum transition regime, where the Navier-Stokes equations become inaccurate and Monte Carlo methods become too computationally expensive. The two-dimensional Burnett equations and the three-dimensional rotational energy transport equation were added to the code and one- and two-dimensional calculations were performed. For the structure of normal shock waves, the Burnett equations give consistently better results than Navier-Stokes equations and compare reasonably well with Monte Carlo methods. For two-dimensional flow of Nitrogen past a circular cylinder the Burnett equations predict the total drag reasonably well. Care must be taken, however, not to exceed the range of validity of the Burnett equations.

A stratospheric aerosol model with perturbations induced by the space shuttle particulate effluents

NASA Technical Reports Server (NTRS)

Rosen, J. M.; Hofmann, D. J.

1977-01-01

A one dimensional steady state stratospheric aerosol model is developed that considers the subsequent perturbations caused by including the expected space shuttle particulate effluents. Two approaches to the basic modeling effort were made: in one, enough simplifying assumptions were introduced so that a more or less exact solution to the descriptive equations could be obtained; in the other approach very few simplifications were made and a computer technique was used to solve the equations. The most complex form of the model contains the effects of sedimentation, diffusion, particle growth and coagulation. Results of the perturbation calculations show that there will probably be an immeasurably small increase in the stratospheric aerosol concentration for particles larger than about 0.15 micrometer radius.

Quantum vacuum interaction between two cosmic strings revisited

NASA Astrophysics Data System (ADS)

Muñoz-Castañeda, J. M.; Bordag, M.

2014-03-01

We reconsider the quantum vacuum interaction energy between two straight parallel cosmic strings. This problem was discussed several times in an approach treating both strings perturbatively and treating only one perturbatively. Here we point out that a simplifying assumption made by Bordag [Ann. Phys. (Berlin) 47, 93 (1990).] can be justified and show that, despite the global character of the background, the perturbative approach delivers a correct result. We consider the applicability of the scattering methods, developed in the past decade for the Casimir effect, for the cosmic string and find it not applicable. We calculate the scattering T-operator on one string. Finally, we consider the vacuum interaction of two strings when each carries a two-dimensional delta function potential.

Chirping and Sudden Excitation of Energetic-Particle-Driven Geodesic Acoustic Modes in a Large Helical Device Experiment

NASA Astrophysics Data System (ADS)

Wang, Hao; Todo, Yasushi; Ido, Takeshi; Suzuki, Yasuhiro

2018-04-01

Energetic-particle-driven geodesic acoustic modes (EGAMs) observed in a Large Helical Device experiment are investigated using a hybrid simulation code for energetic particles interacting with a magnetohydrodynamic (MHD) fluid. The frequency chirping of the primary mode and the sudden excitation of the half-frequency secondary mode are reproduced for the first time with the hybrid simulation using the realistic physical condition and the three-dimensional equilibrium. Both EGAMs have global spatial profiles which are consistent with the experimental measurements. For the secondary mode, the bulk pressure perturbation and the energetic particle pressure perturbation cancel each other out, and thus the frequency is lower than the primary mode. It is found that the excitation of the secondary mode does not depend on the nonlinear MHD coupling. The secondary mode is excited by energetic particles that satisfy the linear and nonlinear resonance conditions, respectively, for the primary and secondary modes.

Chirping and Sudden Excitation of Energetic-Particle-Driven Geodesic Acoustic Modes in a Large Helical Device Experiment.

PubMed

Wang, Hao; Todo, Yasushi; Ido, Takeshi; Suzuki, Yasuhiro

2018-04-27

Energetic-particle-driven geodesic acoustic modes (EGAMs) observed in a Large Helical Device experiment are investigated using a hybrid simulation code for energetic particles interacting with a magnetohydrodynamic (MHD) fluid. The frequency chirping of the primary mode and the sudden excitation of the half-frequency secondary mode are reproduced for the first time with the hybrid simulation using the realistic physical condition and the three-dimensional equilibrium. Both EGAMs have global spatial profiles which are consistent with the experimental measurements. For the secondary mode, the bulk pressure perturbation and the energetic particle pressure perturbation cancel each other out, and thus the frequency is lower than the primary mode. It is found that the excitation of the secondary mode does not depend on the nonlinear MHD coupling. The secondary mode is excited by energetic particles that satisfy the linear and nonlinear resonance conditions, respectively, for the primary and secondary modes.

Twisting perturbed parafermions

NASA Astrophysics Data System (ADS)

Belitsky, A. V.

2017-07-01

The near-collinear expansion of scattering amplitudes in maximally supersymmetric Yang-Mills theory at strong coupling is governed by the dynamics of stings propagating on the five sphere. The pentagon transitions in the operator product expansion which systematize the series get reformulated in terms of matrix elements of branch-point twist operators in the two-dimensional O(6) nonlinear sigma model. The facts that the latter is an asymptotically free field theory and that there exists no local realization of twist fields prevents one from explicit calculation of their scaling dimensions and operator product expansion coefficients. This complication is bypassed making use of the equivalence of the sigma model to the infinite-level limit of WZNW models perturbed by current-current interactions, such that one can use conformal symmetry and conformal perturbation theory for systematic calculations. Presently, to set up the formalism, we consider the O(3) sigma model which is reformulated as perturbed parafermions.

Homotopy decomposition method for solving one-dimensional time-fractional diffusion equation

NASA Astrophysics Data System (ADS)

Abuasad, Salah; Hashim, Ishak

2018-04-01

In this paper, we present the homotopy decomposition method with a modified definition of beta fractional derivative for the first time to find exact solution of one-dimensional time-fractional diffusion equation. In this method, the solution takes the form of a convergent series with easily computable terms. The exact solution obtained by the proposed method is compared with the exact solution obtained by using fractional variational homotopy perturbation iteration method via a modified Riemann-Liouville derivative.

3D Multispecies Nonlinear Perturbative Particle Simulation of Intense Nonneutral Particle Beams (Research supported by the Department of Energy and the Short Pulse Spallation Source Project and LANSCE Division of LANL.)

NASA Astrophysics Data System (ADS)

Qin, Hong; Davidson, Ronald C.; Lee, W. Wei-Li

1999-11-01

The Beam Equilibrium Stability and Transport (BEST) code, a 3D multispecies nonlinear perturbative particle simulation code, has been developed to study collective effects in intense charged particle beams described self-consistently by the Vlasov-Maxwell equations. A Darwin model is adopted for transverse electromagnetic effects. As a 3D multispecies perturbative particle simulation code, it provides several unique capabilities. Since the simulation particles are used to simulate only the perturbed distribution function and self-fields, the simulation noise is reduced significantly. The perturbative approach also enables the code to investigate different physics effects separately, as well as simultaneously. The code can be easily switched between linear and nonlinear operation, and used to study both linear stability properties and nonlinear beam dynamics. These features, combined with 3D and multispecies capabilities, provides an effective tool to investigate the electron-ion two-stream instability, periodically focused solutions in alternating focusing fields, and many other important problems in nonlinear beam dynamics and accelerator physics. Applications to the two-stream instability are presented.

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

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.

Casimir forces between defects in one-dimensional quantum liquids

NASA Astrophysics Data System (ADS)

Recati, A.; Fuchs, J. N.; Peça, C. S.; Zwerger, W.

2005-08-01

We discuss the effective interactions between two localized perturbations in one-dimensional quantum liquids. For noninteracting fermions, the interactions exhibit Friedel oscillations, giving rise to a Ruderman-Kittel-Kasuya-Yosida-type interaction familiar from impurity spins in metals. In the interacting case, at low energies, a Luttinger-liquid description applies. In the case of repulsive fermions, the Friedel oscillations of the interacting system are replaced, at long distances, by a universal Casimir-type interaction which depends only on the sound velocity and decays inversely with the separation. The Casimir-type interaction between localized perturbations embedded in a fermionic environment gives rise to a long-range coupling between quantum dots in ultracold Fermi gases, opening an alternative to couple qubits with neutral atoms. We also briefly discuss the case of bosonic quantum liquids in which the interaction between weak impurities turns out to be short ranged, decaying exponentially on the scale of the healing length.

Using the NASA GRC Sectored-One-Dimensional Combustor Simulation

NASA Technical Reports Server (NTRS)

Paxson, Daniel E.; Mehta, Vishal R.

2014-01-01

The document is a user manual for the NASA GRC Sectored-One-Dimensional (S-1-D) Combustor Simulation. It consists of three sections. The first is a very brief outline of the mathematical and numerical background of the code along with a description of the non-dimensional variables on which it operates. The second section describes how to run the code and includes an explanation of the input file. The input file contains the parameters necessary to establish an operating point as well as the associated boundary conditions (i.e. how it is fed and terminated) of a geometrically configured combustor. It also describes the code output. The third section describes the configuration process and utilizes a specific example combustor to do so. Configuration consists of geometrically describing the combustor (section lengths, axial locations, and cross sectional areas) and locating the fuel injection point and flame region. Configuration requires modifying the source code and recompiling. As such, an executable utility is included with the code which will guide the requisite modifications and insure that they are done correctly.

Three-Dimensional Hydrodynamic Simulations of the Effects of Laser Imprint in OMEGA Implosions

NASA Astrophysics Data System (ADS)

Igumenshchev, I. V.; Campbell, E. M.; Goncharov, V. N.; Regan, S. P.; Shvydky, A.; Schmitt, A. J.

2017-10-01

Illumination of direct-drive implosion targets by the OMEGA laser introduces large-amplitude broadband modulations in the absorbed energy from the largest (target size 900- μm) to smallest (speckle size 2- μm) spatial scales. These modulations ``imprint'' perturbations into a target that are amplified because of the secular and Rayleigh-Taylor growths during acceleration and deceleration of the target. The degradation of performance of room-temperature and cryogenic OMEGA implosions caused by these perturbations were simulated in three dimensions using the code ASTER. The highest-resolution simulations resolve perturbation modes as high as l 200 . The high modes l 50to 100 dominate in the perturbation spectrum during the linear growth, while the late-time nonlinear evolution results in domination of modes with l 30to 50 . Smoothing by spectral dispersion reduces the linear-phase mode amplitudes by a factor of 4 and results in substantial improvements in implosion performance that is in good agreement with measurements. The effects of imprint on implosion performance are compared with the effects of other implosion asymmetries, such as those induced because of laser beam imbalance, mistiming and mispointing, and target offset. This material is based upon work supported by the Department of Energy National Nuclear Security Administration under Award Number DE-NA0001944.

Eigenstates and dynamics of Hooke's atom: Exact results and path integral simulations

NASA Astrophysics Data System (ADS)

Gholizadehkalkhoran, Hossein; Ruokosenmäki, Ilkka; Rantala, Tapio T.

2018-05-01

The system of two interacting electrons in one-dimensional harmonic potential or Hooke's atom is considered, again. On one hand, it appears as a model for quantum dots in a strong confinement regime, and on the other hand, it provides us with a hard test bench for new methods with the "space splitting" arising from the one-dimensional Coulomb potential. Here, we complete the numerous previous studies of the ground state of Hooke's atom by including the excited states and dynamics, not considered earlier. With the perturbation theory, we reach essentially exact eigenstate energies and wave functions for the strong confinement regime as novel results. We also consider external perturbation induced quantum dynamics in a simple separable case. Finally, we test our novel numerical approach based on real-time path integrals (RTPIs) in reproducing the above. The RTPI turns out to be a straightforward approach with exact account of electronic correlations for solving the eigenstates and dynamics without the conventional restrictions of electronic structure methods.

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.

Asteroid mass estimation using Markov-Chain Monte Carlo techniques

NASA Astrophysics Data System (ADS)

Siltala, Lauri; Granvik, Mikael

2016-10-01

Estimates for asteroid masses are based on their gravitational perturbations on the orbits of other objects such as Mars, spacecraft, or other asteroids and/or their satellites. In the case of asteroid-asteroid perturbations, this leads to a 13-dimensional inverse problem where the aim is to derive the mass of the perturbing asteroid and six orbital elements for both the perturbing asteroid and the test asteroid using astrometric observations. We have developed and implemented three different mass estimation algorithms utilizing asteroid-asteroid perturbations into the OpenOrb asteroid-orbit-computation software: the very rough 'marching' approximation, in which the asteroid orbits are fixed at a given epoch, reducing the problem to a one-dimensional estimation of the mass, an implementation of the Nelder-Mead simplex method, and most significantly, a Markov-Chain Monte Carlo (MCMC) approach. We will introduce each of these algorithms with particular focus on the MCMC algorithm, and present example results for both synthetic and real data. Our results agree with the published mass estimates, but suggest that the published uncertainties may be misleading as a consequence of using linearized mass-estimation methods. Finally, we discuss remaining challenges with the algorithms as well as future plans, particularly in connection with ESA's Gaia mission.

Algorithm for loading shot noise microbunching in multi-dimensional, free-electron laser simulation codes

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

Fawley, William M.

We discuss the underlying reasoning behind and the details of the numerical algorithm used in the GINGER free-electron laser(FEL) simulation code to load the initial shot noise microbunching on the electron beam. In particular, we point out that there are some additional subtleties which must be followed for multi-dimensional codes which are not necessary for one-dimensional formulations. Moreover, requiring that the higher harmonics of the microbunching also be properly initialized with the correct statistics leads to additional complexities. We present some numerical results including the predicted incoherent, spontaneous emission as tests of the shot noise algorithm's correctness.

streamgap-pepper: Effects of peppering streams with many small impacts

NASA Astrophysics Data System (ADS)

Bovy, Jo; Erkal, Denis; Sanders, Jason

2017-02-01

streamgap-pepper computes the effect of subhalo fly-bys on cold tidal streams based on the action-angle representation of streams. A line-of-parallel-angle approach is used to calculate the perturbed distribution function of a given stream segment by undoing the effect of all impacts. This approach allows one to compute the perturbed stream density and track in any coordinate system in minutes for realizations of the subhalo distribution down to 10^5 Msun, accounting for the stream's internal dispersion and overlapping impacts. This code uses galpy (ascl:1411.008) and the streampepperdf.py galpy extension, which implements the fast calculation of the perturbed stream structure.

Global Culture: A Noise Induced Transition in Finite Systems

NASA Astrophysics Data System (ADS)

Klemm, Konstantin; Eguíluz, Victor M.; Toral, Raúl; San Miguel, Maxi

2003-04-01

We analyze Axelrod's model for the unbiased transmission of culture in the presence of noise. In a one-dimensional lattice, the dynamics is described in terms of a Lyapunov potential, where the disordered configurations are metastable states of the dynamics. In a two-dimensional lattice the dynamics is governed by the average relaxation time T for perturbations to the homogeneous configuration. If the noise rate is smaller than 1/T, the perturbations drive the system to a completely ordered configuration, whereas the system remains disordered for larger noise rates. Based on a mean-field approximation we obtain the average relaxation time T(N) = Nln(N) for system size N. Thus in the limit of infinite system size the system is disordered for any finite noise rate.

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.

Computer model of one-dimensional equilibrium controlled sorption processes

USGS Publications Warehouse

Grove, D.B.; Stollenwerk, K.G.

1984-01-01

A numerical solution to the one-dimensional solute-transport equation with equilibrium-controlled sorption and a first-order irreversible-rate reaction is presented. The computer code is written in FORTRAN language, with a variety of options for input and output for user ease. Sorption reactions include Langmuir, Freundlich, and ion-exchange, with or without equal valance. General equations describing transport and reaction processes are solved by finite-difference methods, with nonlinearities accounted for by iteration. Complete documentation of the code, with examples, is included. (USGS)

HUFF, a One-Dimensional Hydrodynamics Code for Strong Shocks

DTIC Science & Technology

1978-12-01

results for two sample problems. The first problem discussed is a one-kiloton nuclear burst in infinite sea level air. The second problem is the one...of HUFF as an effective first order hydro- dynamic computer code. 1 KT Explosion The one-kiloton nuclear explosion in infinite sea level air was

Re-forming supercritical quasi-parallel shocks. I - One- and two-dimensional simulations

NASA Technical Reports Server (NTRS)

Thomas, V. A.; Winske, D.; Omidi, N.

1990-01-01

The process of reforming supercritical quasi-parallel shocks is investigated using one-dimensional and two-dimensional hybrid (particle ion, massless fluid electron) simulations both of shocks and of simpler two-stream interactions. It is found that the supercritical quasi-parallel shock is not steady. Instread of a well-defined shock ramp between upstream and downstream states that remains at a fixed position in the flow, the ramp periodically steepens, broadens, and then reforms upstream of its former position. It is concluded that the wave generation process is localized at the shock ramp and that the reformation process proceeds in the absence of upstream perturbations intersecting the shock.

DSMC Studies of the Richtmyer-Meshkov Instability

NASA Astrophysics Data System (ADS)

Gallis, M. A.; Koehler, T. P.; Torczynski, J. R.

2014-11-01

A new exascale-capable Direct Simulation Monte Carlo (DSMC) code, SPARTA, developed to be highly efficient on massively parallel computers, has extended the applicability of DSMC to challenging, transient three-dimensional problems in the continuum regime. Because DSMC inherently accounts for compressibility, viscosity, and diffusivity, it has the potential to improve the understanding of the mechanisms responsible for hydrodynamic instabilities. Here, the Richtmyer-Meshkov instability at the interface between two gases was studied parametrically using SPARTA. Simulations performed on Sequoia, an IBM Blue Gene/Q supercomputer at Lawrence Livermore National Laboratory, are used to investigate various Atwood numbers (0.33-0.94) and Mach numbers (1.2-12.0) for two-dimensional and three-dimensional perturbations. Comparisons with theoretical predictions demonstrate that DSMC accurately predicts the early-time growth of the instability. Sandia National Laboratories is a multi-program laboratory managed and operated by Sandia Corporation, a wholly owned subsidiary of Lockheed Martin Corporation, for the U.S. Department of Energy's National Nuclear Security Administration under contract DE-AC04-94AL85000.

Response solutions and quasi-periodic degenerate bifurcations for quasi-periodically forced systems

NASA Astrophysics Data System (ADS)

Si, Wen; Si, Jianguo

2018-06-01

This paper includes two parts. In the first part, we first focus on quasi-periodic time dependent perturbations of one-dimensional quasi-periodically forced systems with degenerate equilibrium. We study the system in two cases, for one of which system admits a response solution under a non-resonant condition on the frequency vector weaker than Brjuno–Rüssmann’s and for another of which system also admits a response solution without any non-resonant conditions. Next, we investigate the existence of response solutions of a quasi-periodic perturbed system with degenerate (including completely degenerate) equilibrium under Brjuno–Rüssmann’s non-resonant condition by using the Herman method. In the second part, we consider, firstly, the quasi-periodic perturbation of a universal unfolding of one-dimensional degenerate vector field . Secondly, we consider the perturbation of a universal unfolding of normal two-dimensional Hamiltonian system with completely degenerate equilibrium. With KAM theory and singularity theory, we show that these two classes of universal unfolding can persist on large Cantor sets under Brjuno–Rüssmann’s non-resonant condition, which implies all the invariant tori in the integrable part and all the bifurcation scenario can survive on large Cantor sets. The result for Hamiltonian system can apply directly to the response context for quasi-periodically forced systems. Our results in this paper can be regarded as an improvement with respect to several results in various literature (Broer et al 2005 Nonlinearity 18 1735–69 Broer et al 2006 J. Differ. Equ. 222 233–62 Wagener 2005 J. Differ. Equ. 216 216–81 Xu 2010 J. Differ. Equ. 250 551–71 Xu and Jiang 2010 Ergod. Theor. Dynam. Syst. 31 599–611 Lu and Xu 2014 Nonlinear Differ. Equ. Appl. 21 361–70). This work was partially supported by the National Natural Science Foundation of China (Grant Nos. 11171185, 11571201).

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.

SRGULL - AN ADVANCED ENGINEERING MODEL FOR THE PREDICTION OF AIRFRAME INTEGRATED SCRAMJET CYCLE PERFORMANCE

NASA Technical Reports Server (NTRS)

Walton, J. T.

1994-01-01

The development of a single-stage-to-orbit aerospace vehicle intended to be launched horizontally into low Earth orbit, such as the National Aero-Space Plane (NASP), has concentrated on the use of the supersonic combustion ramjet (scramjet) propulsion cycle. SRGULL, a scramjet cycle analysis code, is an engineer's tool capable of nose-to-tail, hydrogen-fueled, airframe-integrated scramjet simulation in a real gas flow with equilibrium thermodynamic properties. This program facilitates initial estimates of scramjet cycle performance by linking a two-dimensional forebody, inlet and nozzle code with a one-dimensional combustor code. Five computer codes (SCRAM, SEAGUL, INLET, Progam HUD, and GASH) originally developed at NASA Langley Research Center in support of hypersonic technology are integrated in this program to analyze changing flow conditions. The one-dimensional combustor code is based on the combustor subroutine from SCRAM and the two-dimensional coding is based on an inviscid Euler program (SEAGUL). Kinetic energy efficiency input for sidewall area variation modeling can be calculated by the INLET program code. At the completion of inviscid component analysis, Program HUD, an integral boundary layer code based on the Spaulding-Chi method, is applied to determine the friction coefficient which is then used in a modified Reynolds Analogy to calculate heat transfer. Real gas flow properties such as flow composition, enthalpy, entropy, and density are calculated by the subroutine GASH. Combustor input conditions are taken from one-dimensionalizing the two-dimensional inlet exit flow. The SEAGUL portions of this program are limited to supersonic flows, but the combustor (SCRAM) section can handle supersonic and dual-mode operation. SRGULL has been compared to scramjet engine tests with excellent results. SRGULL was written in FORTRAN 77 on an IBM PC compatible using IBM's FORTRAN/2 or Microway's NDP386 F77 compiler. The program is fully user interactive, but can also run in batch mode. It operates under the UNIX, VMS, NOS, and DOS operating systems. The source code is not directly compatible with all PC compilers (e.g., Lahey or Microsoft FORTRAN) due to block and segment size requirements. SRGULL executable code requires about 490K RAM and a math coprocessor on PC's. The SRGULL program was developed in 1989, although the component programs originated in the 1960's and 1970's. IBM, IBM PC, and DOS are registered trademarks of International Business Machines. VMS is a registered trademark of Digital Equipment Corporation. UNIX is a registered trademark of Bell Laboratories. NOS is a registered trademark of Control Data Corporation.

A Non Local Electron Heat Transport Model for Multi-Dimensional Fluid Codes

NASA Astrophysics Data System (ADS)

Schurtz, Guy

2000-10-01

Apparent inhibition of thermal heat flow is one of the most ancient problems in computational Inertial Fusion and flux-limited Spitzer-Harm conduction has been a mainstay in multi-dimensional hydrodynamic codes for more than 25 years. Theoretical investigation of the problem indicates that heat transport in laser produced plasmas has to be considered as a non local process. Various authors contributed to the non local theory and proposed convolution formulas designed for practical implementation in one-dimensional fluid codes. Though the theory, confirmed by kinetic calculations, actually predicts a reduced heat flux, it fails to explain the very small limiters required in two-dimensional simulations. Fokker-Planck simulations by Epperlein, Rickard and Bell [PRL 61, 2453 (1988)] demonstrated that non local effects could lead to a strong reduction of heat flow in two dimensions, even in situations where a one-dimensional analysis suggests that the heat flow is nearly classical. We developed at CEA/DAM a non local electron heat transport model suitable for implementation in our two-dimensional radiation hydrodynamic code FCI2. This model may be envisionned as the first step of an iterative solution of the Fokker-Planck equations; it takes the mathematical form of multigroup diffusion equations, the solution of which yields both the heat flux and the departure of the electron distribution function to the Maxwellian. Although direct implementation of the model is straightforward, formal solutions of it can be expressed in convolution form, exhibiting a three-dimensional tensor propagator. Reduction to one dimension retrieves the original formula of Luciani, Mora and Virmont [PRL 51, 1664 (1983)]. Intense magnetic fields may be generated by thermal effects in laser targets; these fields, as well as non local effects, will inhibit electron conduction. We present simulations where both effects are taken into account and shortly discuss the coupling strategy between them.

TWANG-PIC, a novel gyro-averaged one-dimensional particle-in-cell code for interpretation of gyrotron experiments

NASA Astrophysics Data System (ADS)

Braunmueller, F.; Tran, T. M.; Vuillemin, Q.; Alberti, S.; Genoud, J.; Hogge, J.-Ph.; Tran, M. Q.

2015-06-01

A new gyrotron simulation code for simulating the beam-wave interaction using a monomode time-dependent self-consistent model is presented. The new code TWANG-PIC is derived from the trajectory-based code TWANG by describing the electron motion in a gyro-averaged one-dimensional Particle-In-Cell (PIC) approach. In comparison to common PIC-codes, it is distinguished by its computation speed, which makes its use in parameter scans and in experiment interpretation possible. A benchmark of the new code is presented as well as a comparative study between the two codes. This study shows that the inclusion of a time-dependence in the electron equations, as it is the case in the PIC-approach, is mandatory for simulating any kind of non-stationary oscillations in gyrotrons. Finally, the new code is compared with experimental results and some implications of the violated model assumptions in the TWANG code are disclosed for a gyrotron experiment in which non-stationary regimes have been observed and for a critical case that is of interest in high power gyrotron development.

TWANG-PIC, a novel gyro-averaged one-dimensional particle-in-cell code for interpretation of gyrotron experiments

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

Braunmueller, F., E-mail: falk.braunmueller@epfl.ch; Tran, T. M.; Alberti, S.

A new gyrotron simulation code for simulating the beam-wave interaction using a monomode time-dependent self-consistent model is presented. The new code TWANG-PIC is derived from the trajectory-based code TWANG by describing the electron motion in a gyro-averaged one-dimensional Particle-In-Cell (PIC) approach. In comparison to common PIC-codes, it is distinguished by its computation speed, which makes its use in parameter scans and in experiment interpretation possible. A benchmark of the new code is presented as well as a comparative study between the two codes. This study shows that the inclusion of a time-dependence in the electron equations, as it is themore » case in the PIC-approach, is mandatory for simulating any kind of non-stationary oscillations in gyrotrons. Finally, the new code is compared with experimental results and some implications of the violated model assumptions in the TWANG code are disclosed for a gyrotron experiment in which non-stationary regimes have been observed and for a critical case that is of interest in high power gyrotron development.« less

QR Codes 101

ERIC Educational Resources Information Center

Crompton, Helen; LaFrance, Jason; van 't Hooft, Mark

2012-01-01

A QR (quick-response) code is a two-dimensional scannable code, similar in function to a traditional bar code that one might find on a product at the supermarket. The main difference between the two is that, while a traditional bar code can hold a maximum of only 20 digits, a QR code can hold up to 7,089 characters, so it can contain much more…

Stirling Analysis Comparison of Commercial vs. High-Order Methods

NASA Technical Reports Server (NTRS)

Dyson, Rodger W.; Wilson, Scott D.; Tew, Roy C.; Demko, Rikako

2007-01-01

Recently, three-dimensional Stirling engine simulations have been accomplished utilizing commercial Computational Fluid Dynamics software. The validations reported can be somewhat inconclusive due to the lack of precise time accurate experimental results from engines, export control/ proprietary concerns, and the lack of variation in the methods utilized. The last issue may be addressed by solving the same flow problem with alternate methods. In this work, a comprehensive examination of the methods utilized in the commercial codes is compared with more recently developed high-order methods. Specifically, Lele's Compact scheme and Dyson s Ultra Hi-Fi method will be compared with the SIMPLE and PISO methods currently employed in CFD-ACE, FLUENT, CFX, and STAR-CD (all commercial codes which can in theory solve a three-dimensional Stirling model although sliding interfaces and their moving grids limit the effective time accuracy). We will initially look at one-dimensional flows since the current standard practice is to design and optimize Stirling engines with empirically corrected friction and heat transfer coefficients in an overall one-dimensional model. This comparison provides an idea of the range in which commercial CFD software for modeling Stirling engines may be expected to provide accurate results. In addition, this work provides a framework for improving current one-dimensional analysis codes.

Stirling Analysis Comparison of Commercial Versus High-Order Methods

NASA Technical Reports Server (NTRS)

Dyson, Rodger W.; Wilson, Scott D.; Tew, Roy C.; Demko, Rikako

2005-01-01

Recently, three-dimensional Stirling engine simulations have been accomplished utilizing commercial Computational Fluid Dynamics software. The validations reported can be somewhat inconclusive due to the lack of precise time accurate experimental results from engines, export control/proprietary concerns, and the lack of variation in the methods utilized. The last issue may be addressed by solving the same flow problem with alternate methods. In this work, a comprehensive examination of the methods utilized in the commercial codes is compared with more recently developed high-order methods. Specifically, Lele's compact scheme and Dyson's Ultra Hi-Fi method will be compared with the SIMPLE and PISO methods currently employed in CFD-ACE, FLUENT, CFX, and STAR-CD (all commercial codes which can in theory solve a three-dimensional Stirling model with sliding interfaces and their moving grids limit the effective time accuracy). We will initially look at one-dimensional flows since the current standard practice is to design and optimize Stirling engines with empirically corrected friction and heat transfer coefficients in an overall one-dimensional model. This comparison provides an idea of the range in which commercial CFD software for modeling Stirling engines may be expected to provide accurate results. In addition, this work provides a framework for improving current one-dimensional analysis codes.

Verification of low-Mach number combustion codes using the method of manufactured solutions

NASA Astrophysics Data System (ADS)

Shunn, Lee; Ham, Frank; Knupp, Patrick; Moin, Parviz

2007-11-01

Many computational combustion models rely on tabulated constitutive relations to close the system of equations. As these reactive state-equations are typically multi-dimensional and highly non-linear, their implications on the convergence and accuracy of simulation codes are not well understood. In this presentation, the effects of tabulated state-relationships on the computational performance of low-Mach number combustion codes are explored using the method of manufactured solutions (MMS). Several MMS examples are developed and applied, progressing from simple one-dimensional configurations to problems involving higher dimensionality and solution-complexity. The manufactured solutions are implemented in two multi-physics hydrodynamics codes: CDP developed at Stanford University and FUEGO developed at Sandia National Laboratories. In addition to verifying the order-of-accuracy of the codes, the MMS problems help highlight certain robustness issues in existing variable-density flow-solvers. Strategies to overcome these issues are briefly discussed.

Three dimensional δf simulations of beams in the SSC

NASA Astrophysics Data System (ADS)

Koga, J.; Tajima, T.; Machida, S.

1993-12-01

A three dimensional δf strong-strong algorithm has been developed to apply to the study of such effects as space charge and beam-beam interaction phenomena in the Superconducting Super Collider (SSC). The algorithm is obtained from the merging of the particle tracking code Simpsons used for 3 dimensional space charge effects and a δf code. The δf method is used to follow the evolution of the non-gaussian part of the beam distribution. The advantages of this method are twofold. First, the Simpsons code utilizes a realistic accelerator model including synchrotron oscillations and energy ramping in 6 dimensional phase space with electromagnetic fields of the beams calculated using a realistic 3 dimensional field solver. Second, the beams are evolving in the fully self-consistent strong-strong sense with finite particle fluctuation noise is greatly reduced as opposed to the weak-strong models where one beam is fixed.

A generic efficient adaptive grid scheme for rocket propulsion modeling

NASA Technical Reports Server (NTRS)

Mo, J. D.; Chow, Alan S.

1993-01-01

The objective of this research is to develop an efficient, time-accurate numerical algorithm to discretize the Navier-Stokes equations for the predictions of internal one-, two-dimensional and axisymmetric flows. A generic, efficient, elliptic adaptive grid generator is implicitly coupled with the Lower-Upper factorization scheme in the development of ALUNS computer code. The calculations of one-dimensional shock tube wave propagation and two-dimensional shock wave capture, wave-wave interactions, shock wave-boundary interactions show that the developed scheme is stable, accurate and extremely robust. The adaptive grid generator produced a very favorable grid network by a grid speed technique. This generic adaptive grid generator is also applied in the PARC and FDNS codes and the computational results for solid rocket nozzle flowfield and crystal growth modeling by those codes will be presented in the conference, too. This research work is being supported by NASA/MSFC.

Development of a turbomachinery design optimization procedure using a multiple-parameter nonlinear perturbation method

NASA Technical Reports Server (NTRS)

Stahara, S. S.

1984-01-01

An investigation was carried out to complete the preliminary development of a combined perturbation/optimization procedure and associated computational code for designing optimized blade-to-blade profiles of turbomachinery blades. The overall purpose of the procedures developed is to provide demonstration of a rapid nonlinear perturbation method for minimizing the computational requirements associated with parametric design studies of turbomachinery flows. The method combines the multiple parameter nonlinear perturbation method, successfully developed in previous phases of this study, with the NASA TSONIC blade-to-blade turbomachinery flow solver, and the COPES-CONMIN optimization procedure into a user's code for designing optimized blade-to-blade surface profiles of turbomachinery blades. Results of several design applications and a documented version of the code together with a user's manual are provided.

NASA One-Dimensional Combustor Simulation--User Manual for S1D_ML

NASA Technical Reports Server (NTRS)

Stueber, Thomas J.; Paxson, Daniel E.

2014-01-01

The work presented in this paper is to promote research leading to a closed-loop control system to actively suppress thermo-acoustic instabilities. To serve as a model for such a closed-loop control system, a one-dimensional combustor simulation composed using MATLAB software tools has been written. This MATLAB based process is similar to a precursor one-dimensional combustor simulation that was formatted as FORTRAN 77 source code. The previous simulation process requires modification to the FORTRAN 77 source code, compiling, and linking when creating a new combustor simulation executable file. The MATLAB based simulation does not require making changes to the source code, recompiling, or linking. Furthermore, the MATLAB based simulation can be run from script files within the MATLAB environment or with a compiled copy of the executable file running in the Command Prompt window without requiring a licensed copy of MATLAB. This report presents a general simulation overview. Details regarding how to setup and initiate a simulation are also presented. Finally, the post-processing section describes the two types of files created while running the simulation and it also includes simulation results for a default simulation included with the source code.

Continuum modeling of catastrophic collisions

NASA Technical Reports Server (NTRS)

Ryan, Eileen V.; Aspaug, Erik; Melosh, H. J.

1991-01-01

A two dimensional hydrocode based on 2-D SALE was modified to include strength effects and fragmentation equations for fracture resulting from tensile stress in one dimension. Output from this code includes a complete fragmentation summary for each cell of the modeled object: fragment size (mass) distribution, vector velocities of particles, peak values of pressure and tensile stress, and peak strain rates associated with fragmentation. Contour plots showing pressure and temperature at given times within the object are also produced. By invoking axial symmetry, three dimensional events can be modeled such as zero impact parameter collisions between asteroids. The code was tested against the one dimensional model and the analytical solution for a linearly increasing tensile stress under constant strain rate.

One-Dimensional Modelling of Internal Ballistics

NASA Astrophysics Data System (ADS)

Monreal-González, G.; Otón-Martínez, R. A.; Velasco, F. J. S.; García-Cascáles, J. R.; Ramírez-Fernández, F. J.

2017-10-01

A one-dimensional model is introduced in this paper for problems of internal ballistics involving solid propellant combustion. First, the work presents the physical approach and equations adopted. Closure relationships accounting for the physical phenomena taking place during combustion (interfacial friction, interfacial heat transfer, combustion) are deeply discussed. Secondly, the numerical method proposed is presented. Finally, numerical results provided by this code (UXGun) are compared with results of experimental tests and with the outcome from a well-known zero-dimensional code. The model provides successful results in firing tests of artillery guns, predicting with good accuracy the maximum pressure in the chamber and muzzle velocity what highlights its capabilities as prediction/design tool for internal ballistics.

Analysis of a two-dimensional type 6 shock-interference pattern using a perfect-gas code and a real-gas code

NASA Technical Reports Server (NTRS)

Bertin, J. J.; Graumann, B. W.

1973-01-01

Numerical codes were developed to calculate the two dimensional flow field which results when supersonic flow encounters double wedge configurations whose angles are such that a type 4 pattern occurs. The flow field model included the shock interaction phenomena for a delta wing orbiter. Two numerical codes were developed, one which used the perfect gas relations and a second which incorporated a Mollier table to define equilibrium air properties. The two codes were used to generate theoretical surface pressure and heat transfer distributions for velocities from 3,821 feet per second to an entry condition of 25,000 feet per second.

Macroscopic Lagrangian description of warm plasmas. II Nonlinear wave interactions

NASA Technical Reports Server (NTRS)

Kim, H.; Crawford, F. W.

1983-01-01

A macroscopic Lagrangian is simplified to the adiabatic limit and expanded about equilibrium, to third order in perturbation, for three illustrative cases: one-dimensional compression parallel to the static magnetic field, two-dimensional compression perpendicular to the static magnetic field, and three-dimensional compression. As examples of the averaged-Lagrangian method applied to nonlinear wave interactions, coupling coefficients are derived for interactions between two electron plasma waves and an ion acoustic wave, and between an ordinary wave, an electron plasma wave, and an ion acoustic wave.

Study of the early signal perturbations due to GJ and Elves using the LWPC code

NASA Astrophysics Data System (ADS)

Nait Amor, Samir; Ghalila, Hassen; Bouderba, Yasmina

2015-04-01

Early events are a Very Low Frequencies (VLF) signal perturbations recorded during a lightning activity. The properties of these signal perturbations and their association to the lightning peak current and/or Transient Luminous Events (TLEs) were widely studied. In a recently analysis a new early signal perturbations whose recovery time persists for several minutes were discovered. The underlying cause of these events is still unclear. In a recently published work, these events were attributed to the lightning peak current and the type of associated TLE. In others, and newly published papers, analyzes were done where all kind of early events were considered. Statistical results showed that the occurrence of long recovery events is independent of the lightning current amplitude and/or TLEs type. To understand which is the main cause of these events, we analyzed two types of early signal perturbations: One was a typical event (~200s time duration) in association with a Gigantic Jet and the second was a long recovery event in association with an elve recorded on December 12 2009 during the EuroSprite campaign. In addition to the VLF signal analysis, we used the Long Wave Propagation Capability (LWPC) code to simulate the unperturbed and perturbed signal parameters (amplitude and phase), to determine the signal modes attenuation coefficient and then to infer the electron density increases in the disturbed region. The results showed that the reference height was reduced from its ambient value (87km) to 66.4 km in the case of the GJ and 74.3 km for the elve. These reference heights decreases affected the propagating signal at the disturbed region by increasing the modes attenuation coefficient. Effectively, the number of modes was reduced from 28 at ambient condition to 9 modes (in the case of GJ) and 17 (in the case of elve). This high attenuation of modes leads to the appearance of null signal perturbations positions due to the interferences. Between two null positions the signal perturbation was negative (or decreasing) and sometimes positive (or increasing). It was also observed from the LWPC code results that the perturbation amplitude was maximum when the perturbed and unperturbed signals were in phase. Thus the main reason of these observations is the modal structure of the signal at the disturbed region and the receiver location. The electron density increases reached 104 cm-3 at 85 km independently on TLE kind. By the use of the signal perturbation parameters due to the long recovery event and the LWPC code, a recovery time profile of the electron density at each height below 87 km is obtained. The first order exponential decay fit gives different recovery constants depending on the height. This is in good agreement with the atmospheric model where the loss terms rates vary with altitude.

Investigation of surface evolution for stainless steel electrode under pulsed megagauss magnetic field

NASA Astrophysics Data System (ADS)

Zou, Wenkang; Dan, Jiakun; Wang, Guilin; Duan, Shuchao; Wei, Bing; Zhang, Hengdi; Huang, Xianbin; Zhang, Zhaohui; Guo, Fan; Gong, Boyi; Chen, Lin; Wang, Meng; Feng, Shuping; Xie, Weiping; Deng, Jianjun

2018-02-01

Surface evolution for a conductor electrode under pulsed megagauss (MG) magnetic field was investigated. Stainless steel rods with 3 mm diameter were driven by 8 MA, 130 ns (10%-90%) current pulse in a series of shots on the Primary Test Stand. Experimental data from two complementary diagnostic systems and simulation results from one-dimensional magneto-hydrodynamics code reveal a transition phase for instability development. The transition, which begins as the conductor surface starts to expand, lasts about 40 ns in the pulse. It ends after the thermal plasma is formed, and striation electrothermal instability growth stops but magneto-Rayleigh-Taylor instability (MRTI) starts to develop. An expanding velocity which grows to about 2.0 km/s during the transition phase was directly measured for the first time. The threshold magnetic field for thermal plasma formation on the stainless steel surface was inferred to be 3.3 MG under a rising rate of about 66 MG/μs, and after that MRTI becomes predominant for amplitude growth in surface perturbation.

A comparative study of space radiation organ doses and associated cancer risks using PHITS and HZETRN.

PubMed

Bahadori, Amir A; Sato, Tatsuhiko; Slaba, Tony C; Shavers, Mark R; Semones, Edward J; Van Baalen, Mary; Bolch, Wesley E

2013-10-21

NASA currently uses one-dimensional deterministic transport to generate values of the organ dose equivalent needed to calculate stochastic radiation risk following crew space exposures. In this study, organ absorbed doses and dose equivalents are calculated for 50th percentile male and female astronaut phantoms using both the NASA High Charge and Energy Transport Code to perform one-dimensional deterministic transport and the Particle and Heavy Ion Transport Code System to perform three-dimensional Monte Carlo transport. Two measures of radiation risk, effective dose and risk of exposure-induced death (REID) are calculated using the organ dose equivalents resulting from the two methods of radiation transport. For the space radiation environments and simplified shielding configurations considered, small differences (<8%) in the effective dose and REID are found. However, for the galactic cosmic ray (GCR) boundary condition, compensating errors are observed, indicating that comparisons between the integral measurements of complex radiation environments and code calculations can be misleading. Code-to-code benchmarks allow for the comparison of differential quantities, such as secondary particle differential fluence, to provide insight into differences observed in integral quantities for particular components of the GCR spectrum.

A comparative study of space radiation organ doses and associated cancer risks using PHITS and HZETRN

NASA Astrophysics Data System (ADS)

Bahadori, Amir A.; Sato, Tatsuhiko; Slaba, Tony C.; Shavers, Mark R.; Semones, Edward J.; Van Baalen, Mary; Bolch, Wesley E.

2013-10-01

NASA currently uses one-dimensional deterministic transport to generate values of the organ dose equivalent needed to calculate stochastic radiation risk following crew space exposures. In this study, organ absorbed doses and dose equivalents are calculated for 50th percentile male and female astronaut phantoms using both the NASA High Charge and Energy Transport Code to perform one-dimensional deterministic transport and the Particle and Heavy Ion Transport Code System to perform three-dimensional Monte Carlo transport. Two measures of radiation risk, effective dose and risk of exposure-induced death (REID) are calculated using the organ dose equivalents resulting from the two methods of radiation transport. For the space radiation environments and simplified shielding configurations considered, small differences (<8%) in the effective dose and REID are found. However, for the galactic cosmic ray (GCR) boundary condition, compensating errors are observed, indicating that comparisons between the integral measurements of complex radiation environments and code calculations can be misleading. Code-to-code benchmarks allow for the comparison of differential quantities, such as secondary particle differential fluence, to provide insight into differences observed in integral quantities for particular components of the GCR spectrum.

Perturbative Quantum Gravity and its Relation to Gauge Theory.

PubMed

Bern, Zvi

2002-01-01

In this review we describe a non-trivial relationship between perturbative gauge theory and gravity scattering amplitudes. At the semi-classical or tree-level, the scattering amplitudes of gravity theories in flat space can be expressed as a sum of products of well defined pieces of gauge theory amplitudes. These relationships were first discovered by Kawai, Lewellen, and Tye in the context of string theory, but hold more generally. In particular, they hold for standard Einstein gravity. A method based on D -dimensional unitarity can then be used to systematically construct all quantum loop corrections order-by-order in perturbation theory using as input the gravity tree amplitudes expressed in terms of gauge theory ones. More generally, the unitarity method provides a means for perturbatively quantizing massless gravity theories without the usual formal apparatus associated with the quantization of constrained systems. As one application, this method was used to demonstrate that maximally supersymmetric gravity is less divergent in the ultraviolet than previously thought.

Programmers manual for a one-dimensional Lagrangian transport model

USGS Publications Warehouse

Schoellhamer, D.H.; Jobson, H.E.

1986-01-01

A one-dimensional Lagrangian transport model for simulating water-quality constituents such as temperature, dissolved oxygen , and suspended sediment in rivers is presented in this Programmers Manual. Lagrangian transport modeling techniques, the model 's subroutines, and the user-written decay-coefficient subroutine are discussed in detail. Appendices list the program codes. The Programmers Manual is intended for the model user who needs to modify code either to adapt the model to a particular need or to use reaction kinetics not provided with the model. (Author 's abstract)

Radiation hydrodynamic effects in two beryllium plates with an idealized aluminum joint

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

Belkov, S.A.; Mkhitarian, L.S.; Vinokurov, O.A.

A beryllium capsule formed from two hemispherical shells with a thin bond is one possible ignition target for the National Ignition Facility [J. A. Paisner {ital et al.}, Laser Focus World {bold 30}, 75 (1994)] Nonuniformities in density, opacity, and interface position at the joint between these hemishells will initiate two-dimensional (2-D) perturbations of the shock wave and material behind the shock as the shock passes through the shell perpendicular to the joint width. Rarefaction of material flow behind the shock front can cause the interface between the shell and joint material to oscillate in position. The amplitude of thesemore » oscillations may be comparable to the joint width. The evolution of these perturbations is studied by numerically simulating shock passage through flat beryllium plates containing aluminum joints. Using the MIMOSA-ND code [D. Sofronov {ital et al.}, Vopr. At. Nauki Tekh., Ser: Mat. modelirovanie fizicheskih processov {bold 2}, 3 (1990)] two different cases are calculated{emdash}a wide (10 {mu}m) and a narrow (1 {mu}m) joint of aluminum between two 150 {mu}m long semiinfinite beryllium plates. Both cases showed good agreement with an analytic representation of the oscillation behavior. For the narrow joint, a special technique allows the calculation of mixing between the joint and surrounding material caused by the Kelvin{endash}Helmholtz instability. {copyright} {ital 1999 American Institute of Physics.}« less

Cyclic behavior at quasi-parallel collisionless shocks

NASA Technical Reports Server (NTRS)

Burgess, D.

1989-01-01

Large scale one-dimensional hybrid simulations with resistive electrons have been carried out of a quasi-parallel high-Mach-number collisionless shock. The shock initially appears stable, but then exhibits cyclic behavior. For the magnetic field, the cycle consists of a period when the transition from upstream to downstream is steep and well defined, followed by a period when the shock transition is extended and perturbed. This cyclic shock solution results from upstream perturbations caused by backstreaming gyrating ions convecting into the shock. The cyclic reformation of a sharp shock transition can allow ions, at one time upstream because of reflection or leakage, to contribute to the shock thermalization.

Extension of the SIESTA MHD equilibrium code to free-plasma-boundary problems

DOE PAGES

Peraza-Rodriguez, Hugo; Reynolds-Barredo, J. M.; Sanchez, Raul; ...

2017-08-28

Here, SIESTA is a recently developed MHD equilibrium code designed to perform fast and accurate calculations of ideal MHD equilibria for three-dimensional magnetic configurations. Since SIESTA does not assume closed magnetic surfaces, the solution can exhibit magnetic islands and stochastic regions. In its original implementation SIESTA addressed only fixed-boundary problems. That is, the shape of the plasma edge, assumed to be a magnetic surface, was kept fixed as the solution iteratively converges to equilibrium. This condition somewhat restricts the possible applications of SIESTA. In this paper we discuss an extension that will enable SIESTA to address free-plasma-boundary problems, opening upmore » the possibility of investigating problems in which the plasma boundary is perturbed either externally or internally. As an illustration, SIESTA is applied to a configuration of the W7-X stellarator.« less

Extension of the SIESTA MHD equilibrium code to free-plasma-boundary problems

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

Peraza-Rodriguez, Hugo; Reynolds-Barredo, J. M.; Sanchez, Raul

Here, SIESTA is a recently developed MHD equilibrium code designed to perform fast and accurate calculations of ideal MHD equilibria for three-dimensional magnetic configurations. Since SIESTA does not assume closed magnetic surfaces, the solution can exhibit magnetic islands and stochastic regions. In its original implementation SIESTA addressed only fixed-boundary problems. That is, the shape of the plasma edge, assumed to be a magnetic surface, was kept fixed as the solution iteratively converges to equilibrium. This condition somewhat restricts the possible applications of SIESTA. In this paper we discuss an extension that will enable SIESTA to address free-plasma-boundary problems, opening upmore » the possibility of investigating problems in which the plasma boundary is perturbed either externally or internally. As an illustration, SIESTA is applied to a configuration of the W7-X stellarator.« less

Highly Parallel Alternating Directions Algorithm for Time Dependent Problems

NASA Astrophysics Data System (ADS)

Ganzha, M.; Georgiev, K.; Lirkov, I.; Margenov, S.; Paprzycki, M.

2011-11-01

In our work, we consider the time dependent Stokes equation on a finite time interval and on a uniform rectangular mesh, written in terms of velocity and pressure. For this problem, a parallel algorithm based on a novel direction splitting approach is developed. Here, the pressure equation is derived from a perturbed form of the continuity equation, in which the incompressibility constraint is penalized in a negative norm induced by the direction splitting. The scheme used in the algorithm is composed of two parts: (i) velocity prediction, and (ii) pressure correction. This is a Crank-Nicolson-type two-stage time integration scheme for two and three dimensional parabolic problems in which the second-order derivative, with respect to each space variable, is treated implicitly while the other variable is made explicit at each time sub-step. In order to achieve a good parallel performance the solution of the Poison problem for the pressure correction is replaced by solving a sequence of one-dimensional second order elliptic boundary value problems in each spatial direction. The parallel code is implemented using the standard MPI functions and tested on two modern parallel computer systems. The performed numerical tests demonstrate good level of parallel efficiency and scalability of the studied direction-splitting-based algorithm.

Network perturbation by recurrent regulatory variants in cancer

PubMed Central

Cho, Ara; Lee, Insuk; Choi, Jung Kyoon

2017-01-01

Cancer driving genes have been identified as recurrently affected by variants that alter protein-coding sequences. However, a majority of cancer variants arise in noncoding regions, and some of them are thought to play a critical role through transcriptional perturbation. Here we identified putative transcriptional driver genes based on combinatorial variant recurrence in cis-regulatory regions. The identified genes showed high connectivity in the cancer type-specific transcription regulatory network, with high outdegree and many downstream genes, highlighting their causative role during tumorigenesis. In the protein interactome, the identified transcriptional drivers were not as highly connected as coding driver genes but appeared to form a network module centered on the coding drivers. The coding and regulatory variants associated via these interactions between the coding and transcriptional drivers showed exclusive and complementary occurrence patterns across tumor samples. Transcriptional cancer drivers may act through an extensive perturbation of the regulatory network and by altering protein network modules through interactions with coding driver genes. PMID:28333928

Travel-time sensitivity kernels in long-range propagation.

PubMed

Skarsoulis, E K; Cornuelle, B D; Dzieciuch, M A

2009-11-01

Wave-theoretic travel-time sensitivity kernels (TSKs) are calculated in two-dimensional (2D) and three-dimensional (3D) environments and their behavior with increasing propagation range is studied and compared to that of ray-theoretic TSKs and corresponding Fresnel-volumes. The differences between the 2D and 3D TSKs average out when horizontal or cross-range marginals are considered, which indicates that they are not important in the case of range-independent sound-speed perturbations or perturbations of large scale compared to the lateral TSK extent. With increasing range, the wave-theoretic TSKs expand in the horizontal cross-range direction, their cross-range extent being comparable to that of the corresponding free-space Fresnel zone, whereas they remain bounded in the vertical. Vertical travel-time sensitivity kernels (VTSKs)-one-dimensional kernels describing the effect of horizontally uniform sound-speed changes on travel-times-are calculated analytically using a perturbation approach, and also numerically, as horizontal marginals of the corresponding TSKs. Good agreement between analytical and numerical VTSKs, as well as between 2D and 3D VTSKs, is found. As an alternative method to obtain wave-theoretic sensitivity kernels, the parabolic approximation is used; the resulting TSKs and VTSKs are in good agreement with normal-mode results. With increasing range, the wave-theoretic VTSKs approach the corresponding ray-theoretic sensitivity kernels.

Structures and Intermittency in a Passive Scalar Model

NASA Astrophysics Data System (ADS)

Vergassola, M.; Mazzino, A.

1997-09-01

Perturbative expansions for intermittency scaling exponents in the Kraichnan passive scalar model [Phys. Rev. Lett. 72, 1016 (1994)] are investigated. A one-dimensional compressible model is considered for this purpose. High resolution Monte Carlo simulations using an Ito approach adapted to an advecting velocity field with a very short correlation time are performed and lead to clean scaling behavior for passive scalar structure functions. Perturbative predictions for the scaling exponents around the Gaussian limit of the model are derived as in the Kraichnan model. Their comparison with the simulations indicates that the scale-invariant perturbative scheme correctly captures the inertial range intermittency corrections associated with the intense localized structures observed in the dynamics.

Improving performance of DS-CDMA systems using chaotic complex Bernoulli spreading codes

NASA Astrophysics Data System (ADS)

Farzan Sabahi, Mohammad; Dehghanfard, Ali

2014-12-01

The most important goal of spreading spectrum communication system is to protect communication signals against interference and exploitation of information by unintended listeners. In fact, low probability of detection and low probability of intercept are two important parameters to increase the performance of the system. In Direct Sequence Code Division Multiple Access (DS-CDMA) systems, these properties are achieved by multiplying the data information in spreading sequences. Chaotic sequences, with their particular properties, have numerous applications in constructing spreading codes. Using one-dimensional Bernoulli chaotic sequence as spreading code is proposed in literature previously. The main feature of this sequence is its negative auto-correlation at lag of 1, which with proper design, leads to increase in efficiency of the communication system based on these codes. On the other hand, employing the complex chaotic sequences as spreading sequence also has been discussed in several papers. In this paper, use of two-dimensional Bernoulli chaotic sequences is proposed as spreading codes. The performance of a multi-user synchronous and asynchronous DS-CDMA system will be evaluated by applying these sequences under Additive White Gaussian Noise (AWGN) and fading channel. Simulation results indicate improvement of the performance in comparison with conventional spreading codes like Gold codes as well as similar complex chaotic spreading sequences. Similar to one-dimensional Bernoulli chaotic sequences, the proposed sequences also have negative auto-correlation. Besides, construction of complex sequences with lower average cross-correlation is possible with the proposed method.

Control of Combustion-Instabilities Through Various Passive Devices

NASA Technical Reports Server (NTRS)

Frendi, Kader

2005-01-01

It is well known that under some operating conditions, rocket engines (using solid or liquid fuels) exhibit unstable modes of operation that can lead to engine malfunction and shutdown. The sources of these instabilities are diverse and are dependent on fuel, chamber geometry and various upstream sources such as pumps, valves and injection mechanism. It is believed that combustion-acoustic instabilities occur when the acoustic energy increase due to the unsteady heat release of the flame is greater than the losses of acoustic energy from the system [1, 2]. Giammar and Putnam [3] performed a comprehensive study of noise generated by gasfired industrial burners and made several key observations; flow noise was sometimes more intense than combustion roar, which tended to have a characteristic frequency spectrum. Turbulence was amplified by the flame. The noise power varied directly with combustion intensity and also with the product of pressure drop and heat release rate. Karchmer [4] correlated the noise emitted from a turbofan jet engine with that in the combustion chamber. This is important, since it quantified how much of the noise from an engine originates in the combustor. A physical interpretation of the interchange of energy between sound waves and unsteady heat release rates was given by Rayleigh [5] for inviscid, linear perturbations. Bloxidge et al [6] extended Rayleigh s criterion to describe the interaction of unsteady combustion with one-dimensional acoustic waves in a duct. Solutions to the mass, momentum and energy conservation equations in the pre- and post-flame zones were matched by making several assumptions about the combustion process. They concluded that changes in boundary conditions affect the energy balance of acoustic waves in the combustor. Abouseif et al [7] also solved the one-dimensional flow equations, but they used a onestep reaction to evaluate the unsteady heat release rate by relating it to temperature and velocity perturbations. Their analysis showed that oscillations arise from coupling between entropy waves produced at the flame and pressure waves originating from the nozzle. Yang and Culick [8] assumed a thin flame sheet, which is distorted by velocity and pressure oscillations. Conservation equations were expressed in integral form and solutions for the acoustic wave equations and complex frequencies were obtained. The imaginary part of the frequency indicated stability regions of the flame. Activation energy asymptotics together with a one-step reaction were used by McIntosh [9] to study the effects of acoustic forcing and feedback on unsteady, one-dimensional flames. He found that the flame stability was altered by the upstream acoustic feedback. Shyy et al [10] used a high-accuracy TVD scheme to simulate unsteady, one-dimensional longitudinal, combustion instabilities. However, numerical diffusion was not completely eliminated. Recently, Prasad [11] investigated numerically the interactions of pressure perturbations with premixed flames. He used complex chemistry to study responses of pressure perturbations in one-dimensional combustors. His results indicated that reflected and transmitted waves differed significantly from incident waves.

Decomposition of Magnetic Field Boundary Conditions into Parts Produced by Internal and External Sources

NASA Astrophysics Data System (ADS)

Lazanja, David; Boozer, Allen

2006-10-01

Given the total magnetic field on a toroidal plasma surface, a method for decomposing the field into a part due to internal currents (often the plasma) and a part due to external currents is presented. The method exploits Laplace theory which is valid in the vacuum region between the plasma surface and the chamber walls. The method is developed for the full three dimensional case which is necessary for studying stellarator plasma configurations. A change in the plasma shape is produced by the total normal field perturbation on the plasma surface. This method allows a separation of the total normal field perturbation into a part produced by external currents and a part produced by the plasma response. There are immediate applications to coil design. The computational procedure is based on Merkel's 1986 work on vacuum field computations. Several test cases are presented for toroidal surfaces which verify the method and computational robustness of the code.

A computer program to obtain time-correlated gust loads for nonlinear aircraft using the matched-filter-based method

NASA Technical Reports Server (NTRS)

Scott, Robert C.; Pototzky, Anthony S.; Perry, Boyd, III

1994-01-01

NASA Langley Research Center has, for several years, conducted research in the area of time-correlated gust loads for linear and nonlinear aircraft. The results of this work led NASA to recommend that the Matched-Filter-Based One-Dimensional Search Method be used for gust load analyses of nonlinear aircraft. This manual describes this method, describes a FORTRAN code which performs this method, and presents example calculations for a sample nonlinear aircraft model. The name of the code is MFD1DS (Matched-Filter-Based One-Dimensional Search). The program source code, the example aircraft equations of motion, a sample input file, and a sample program output are all listed in the appendices.

Perturbative Aspects of Low-Dimensional Quantum Field Theory

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

Wardaya, Asep Y.; Theoretical Physics Laboratory, Theoretical High Energy Physics and Instrumentation Research Group, FMIPA, Institut Teknologi Bandung, Jl. Ganesha 10 Bandung 40132; Zen, Freddy P.

We investigate the low-dimensional applications of Quantum Field Theory (QFT), namely Chern-Simons-Witten Theory (CSWT) and Affine Toda Field Theory (ATFT) in 3- and 2- dimensions. We discuss the perturbative aspects of both theories and compare the results to the exact solutions obtained nonperturbatively. For the three dimensions CSWT case, the perturbative term agree with the nonperturbative polynomial invariants up to third order of the coupling constant 1/k. In the two dimensions ATFT, we investigate the perturbative aspect of S-matrices for A{sub 1}{sup (1)} case in eighth order of the coupling constant {beta}.

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.

Morse Monte Carlo Radiation Transport Code System

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

Emmett, M.B.

1975-02-01

The report contains sections containing descriptions of the MORSE and PICTURE codes, input descriptions, sample problems, deviations of the physical equations and explanations of the various error messages. The MORSE code is a multipurpose neutron and gamma-ray transport Monte Carlo code. Time dependence for both shielding and criticality problems is provided. General three-dimensional geometry may be used with an albedo option available at any material surface. The PICTURE code provide aid in preparing correct input data for the combinatorial geometry package CG. It provides a printed view of arbitrary two-dimensional slices through the geometry. By inspecting these pictures one maymore » determine if the geometry specified by the input cards is indeed the desired geometry. 23 refs. (WRF)« less

Linear chirp phase perturbing approach for finding binary phased codes

NASA Astrophysics Data System (ADS)

Li, Bing C.

2017-05-01

Binary phased codes have many applications in communication and radar systems. These applications require binary phased codes to have low sidelobes in order to reduce interferences and false detection. Barker codes are the ones that satisfy these requirements and they have lowest maximum sidelobes. However, Barker codes have very limited code lengths (equal or less than 13) while many applications including low probability of intercept radar, and spread spectrum communication, require much higher code lengths. The conventional techniques of finding binary phased codes in literatures include exhaust search, neural network, and evolutionary methods, and they all require very expensive computation for large code lengths. Therefore these techniques are limited to find binary phased codes with small code lengths (less than 100). In this paper, by analyzing Barker code, linear chirp, and P3 phases, we propose a new approach to find binary codes. Experiments show that the proposed method is able to find long low sidelobe binary phased codes (code length >500) with reasonable computational cost.

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.

One-loop perturbative coupling of A and A? through the chiral overlap operator

NASA Astrophysics Data System (ADS)

Makino, Hiroki; Morikawa, Okuto; Suzuki, Hiroshi

2018-03-01

Recently, Grabowska and Kaplan constructed a four-dimensional lattice formulation of chiral gauge theories on the basis of the chiral overlap operator. At least in the tree-level approximation, the left-handed fermion is coupled only to the original gauge field A, while the right-handed one is coupled only to the gauge field A*, a deformation of A by the gradient flow with infinite flow time. In this paper, we study the fermion one-loop effective action in their formulation. We show that the continuum limit of this effective action contains local interaction terms between A and A*, even if the anomaly cancellation condition is met. These non-vanishing terms would lead an undesired perturbative spectrum in the formulation.

Extended resolvent and inverse scattering with an application to KPI

NASA Astrophysics Data System (ADS)

Boiti, M.; Pempinelli, F.; Pogrebkov, A. K.; Prinari, B.

2003-08-01

We present in detail an extended resolvent approach for investigating linear problems associated to 2+1 dimensional integrable equations. Our presentation is based as an example on the nonstationary Schrödinger equation with potential being a perturbation of the one-soliton potential by means of a decaying two-dimensional function. Modification of the inverse scattering theory as well as properties of the Jost solutions and spectral data as follows from the resolvent approach are given.

XGC developments for a more efficient XGC-GENE code coupling

NASA Astrophysics Data System (ADS)

Dominski, Julien; Hager, Robert; Ku, Seung-Hoe; Chang, Cs

2017-10-01

In the Exascale Computing Program, the High-Fidelity Whole Device Modeling project initially aims at delivering a tightly-coupled simulation of plasma neoclassical and turbulence dynamics from the core to the edge of the tokamak. To permit such simulations, the gyrokinetic codes GENE and XGC will be coupled together. Numerical efforts are made to improve the numerical schemes agreement in the coupling region. One of the difficulties of coupling those codes together is the incompatibility of their grids. GENE is a continuum grid-based code and XGC is a Particle-In-Cell code using unstructured triangular mesh. A field-aligned filter is thus implemented in XGC. Even if XGC originally had an approximately field-following mesh, this field-aligned filter permits to have a perturbation discretization closer to the one solved in the field-aligned code GENE. Additionally, new XGC gyro-averaging matrices are implemented on a velocity grid adapted to the plasma properties, thus ensuring same accuracy from the core to the edge regions.

Development of a thermal and structural analysis procedure for cooled radial turbines

NASA Technical Reports Server (NTRS)

Kumar, Ganesh N.; Deanna, Russell G.

1988-01-01

A procedure for computing the rotor temperature and stress distributions in a cooled radial turbine is considered. Existing codes for modeling the external mainstream flow and the internal cooling flow are used to compute boundary conditions for the heat transfer and stress analyses. An inviscid, quasi three-dimensional code computes the external free stream velocity. The external velocity is then used in a boundary layer analysis to compute the external heat transfer coefficients. Coolant temperatures are computed by a viscous one-dimensional internal flow code for the momentum and energy equation. These boundary conditions are input to a three-dimensional heat conduction code for calculation of rotor temperatures. The rotor stress distribution may be determined for the given thermal, pressure and centrifugal loading. The procedure is applied to a cooled radial turbine which will be tested at the NASA Lewis Research Center. Representative results from this case are included.

Benchmark problems in computational aeroacoustics

NASA Technical Reports Server (NTRS)

Porter-Locklear, Freda

1994-01-01

A recent directive at NASA Langley is aimed at numerically predicting principal noise sources. During my summer stay, I worked with high-order ENO code, developed by Dr. Harold Atkins, for solving the unsteady compressible Navier-Stokes equations, as it applies to computational aeroacoustics (CAA). A CAA workshop, composed of six categories of benchmark problems, has been organized to test various numerical properties of code. My task was to determine the robustness of Atkins' code for these test problems. In one category, we tested the nonlinear wave propagation of the code for the one-dimensional Euler equations, with initial pressure, density, and velocity conditions. Using freestream boundary conditions, our results were plausible. In another category, we solved the linearized two-dimensional Euler equations to test the effectiveness of radiation boundary conditions. Here we utilized MAPLE to compute eigenvalues and eigenvectors of the Jacobian given variable and flux vectors. We experienced a minor problem with inflow and outflow boundary conditions. Next, we solved the quasi one dimensional unsteady flow equations with an incoming acoustic wave of amplitude 10(exp -6). The small amplitude sound wave was incident on a convergent-divergent nozzle. After finding a steady-state solution and then marching forward, our solution indicated that after 30 periods the acoustic wave had dissipated (a period is time required for sound wave to traverse one end of nozzle to other end).

Large perturbation flow field analysis and simulation for supersonic inlets

NASA Technical Reports Server (NTRS)

Varner, M. O.; Martindale, W. R.; Phares, W. J.; Kneile, K. R.; Adams, J. C., Jr.

1984-01-01

An analysis technique for simulation of supersonic mixed compression inlets with large flow field perturbations is presented. The approach is based upon a quasi-one-dimensional inviscid unsteady formulation which includes engineering models of unstart/restart, bleed, bypass, and geometry effects. Numerical solution of the governing time dependent equations of motion is accomplished through a shock capturing finite difference algorithm, of which five separate approaches are evaluated. Comparison with experimental supersonic wind tunnel data is presented to verify the present approach for a wide range of transient inlet flow conditions.

Use Them ... or Lose Them? The Case for and against Using QR Codes

ERIC Educational Resources Information Center

Cunningham, Chuck; Dull, Cassie

2011-01-01

A quick-response (QR) code is a two-dimensional, black-and-white square barcode and links directly to a URL of one's choice. When the code is scanned with a smartphone, it will automatically redirect the user to the designated URL. QR codes are popping up everywhere--billboards, magazines, posters, shop windows, TVs, computer screens, and more.…

Stochastic Stabilityfor Contracting Lorenz Maps and Flows

NASA Astrophysics Data System (ADS)

Metzger, R. J.

In a previous work [M], we proved the existence of absolutely continuous invariant measures for contracting Lorenz-like maps, and constructed Sinai-Ruelle-Bowen measures f or the flows that generate them. Here, we prove stochastic stability for such one-dimensional maps and use this result to prove that the corresponding flows generating these maps are stochastically stable under small diffusion-type perturbations, even though, as shown by Rovella [Ro], they are persistent only in a measure theoretical sense in a parameter space. For the one-dimensional maps we also prove strong stochastic stability in the sense of Baladi and Viana[BV].

Modelling and calculation of flotation process in one-dimensional formulation

NASA Astrophysics Data System (ADS)

Amanbaev, Tulegen; Tilleuov, Gamidulla; Tulegenova, Bibigul

2016-08-01

In the framework of the assumptions of the mechanics of the multiphase media is constructed a mathematical model of the flotation process in the dispersed mixture of liquid, solid and gas phases, taking into account the degree of mineralization of the surface of the bubbles. Application of the constructed model is demonstrated on the example of one-dimensional stationary flotation and it is shown that the equations describing the process of ascent of the bubbles are singularly perturbed ("rigid"). The effect of size and concentration of bubbles and the volumetric content of dispersed particles on the flotation process are analyzed.

HITEMP Material and Structural Optimization Technology Transfer

NASA Technical Reports Server (NTRS)

Collier, Craig S.; Arnold, Steve (Technical Monitor)

2001-01-01

The feasibility of adding viscoelasticity and the Generalized Method of Cells (GMC) for micromechanical viscoelastic behavior into the commercial HyperSizer structural analysis and optimization code was investigated. The viscoelasticity methodology was developed in four steps. First, a simplified algorithm was devised to test the iterative time stepping method for simple one-dimensional multiple ply structures. Second, GMC code was made into a callable subroutine and incorporated into the one-dimensional code to test the accuracy and usability of the code. Third, the viscoelastic time-stepping and iterative scheme was incorporated into HyperSizer for homogeneous, isotropic viscoelastic materials. Finally, the GMC was included in a version of HyperSizer. MS Windows executable files implementing each of these steps is delivered with this report, as well as source code. The findings of this research are that both viscoelasticity and GMC are feasible and valuable additions to HyperSizer and that the door is open for more advanced nonlinear capability, such as viscoplasticity.

Confined One Dimensional Harmonic Oscillator as a Two-Mode System

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

Gueorguiev, V G; Rau, A P; Draayer, J P

2005-07-11

The one-dimensional harmonic oscillator in a box problem is possibly the simplest example of a two-mode system. This system has two exactly solvable limits, the harmonic oscillator and a particle in a (one-dimensional) box. Each of the two limits has a characteristic spectral structure describing the two different excitation modes of the system. Near each of these limits, one can use perturbation theory to achieve an accurate description of the eigenstates. Away from the exact limits, however, one has to carry out a matrix diagonalization because the basis-state mixing that occurs is typically too large to be reproduced in anymore » other way. An alternative to casting the problem in terms of one or the other basis set consists of using an ''oblique'' basis that uses both sets. Through a study of this alternative in this one-dimensional problem, we are able to illustrate practical solutions and infer the applicability of the concept for more complex systems, such as in the study of complex nuclei where oblique-basis calculations have been successful.« less

Pumping approximately integrable systems

PubMed Central

Lange, Florian; Lenarčič, Zala; Rosch, Achim

2017-01-01

Weak perturbations can drive an interacting many-particle system far from its initial equilibrium state if one is able to pump into degrees of freedom approximately protected by conservation laws. This concept has for example been used to realize Bose–Einstein condensates of photons, magnons and excitons. Integrable quantum systems, like the one-dimensional Heisenberg model, are characterized by an infinite set of conservation laws. Here, we develop a theory of weakly driven integrable systems and show that pumping can induce large spin or heat currents even in the presence of integrability breaking perturbations, since it activates local and quasi-local approximate conserved quantities. The resulting steady state is qualitatively captured by a truncated generalized Gibbs ensemble with Lagrange parameters that depend on the structure but not on the overall amplitude of perturbations nor the initial state. We suggest to use spin-chain materials driven by terahertz radiation to realize integrability-based spin and heat pumps. PMID:28598444

Zero Dimensional Field Theory of Tachyon Matter

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

Dimitrijevic, D. D.; Djordjevic, G. S.

2007-04-23

The first issue about the object (now) called tachyons was published almost one century ago. Even though there is no experimental evidence of tachyons there are several reasons why tachyons are still of interest today, in fact interest in tachyons is increasing. Many string theories have tachyons occurring as some of the particles in the theory. In this paper we consider the zero dimensional version of the field theory of tachyon matter proposed by A. Sen. Using perturbation theory and ideas of S. Kar, we demonstrate how this tachyon field theory can be connected with a classical mechanical system, suchmore » as a massive particle moving in a constant field with quadratic friction. The corresponding Feynman path integral form is proposed using a perturbative method. A few promising lines for further applications and investigations are noted.« less

Non-perturbative background field calculations

NASA Astrophysics Data System (ADS)

Stephens, C. R.

1988-01-01

New methods are developed for calculating one loop functional determinants in quantum field theory. Instead of relying on a calculation of all the eigenvalues of the small fluctuation equation, these techniques exploit the ability of the proper time formalism to reformulate an infinite dimensional field theoretic problem into a finite dimensional covariant quantum mechanical analog, thereby allowing powerful tools such as the method of Jacobi fields to be used advantageously in a field theory setting. More generally the methods developed herein should be extremely valuable when calculating quantum processes in non-constant background fields, offering a utilitarian alternative to the two standard methods of calculation—perturbation theory in the background field or taking the background field into account exactly. The formalism developed also allows for the approximate calculation of covariances of partial differential equations from a knowledge of the solutions of a homogeneous ordinary differential equation.

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

Kung, Y. F.; Bazin, C.; Wohlfeld, K.

Using determinant quantum Monte Carlo (DQMC) simulations, we systematically study the doping dependence of the crossover from one to two dimensions and its impact on the magnetic properties of the Hubbard model. A square lattice of chains is used, in which the dimensionality can be tuned by varying the interchain coupling t ⊥. The dynamical spin structure factor and static quantities, such as the static spin susceptibility and nearest-neighbor spin correlation function, are characterized in the one- and two-dimensional limits as a benchmark. When the dimensionality is tuned between these limits, the magnetic properties, while evolving smoothly from one tomore » two dimensions, drastically change regardless of the doping level. This suggests that the spin excitations in the two-dimensional Hubbard model, even in the heavily doped case, cannot be explained using the spinon picture known from one dimension. In conclusion, the DQMC calculations are complemented by cluster perturbation theory studies to form a more complete picture of how the crossover occurs as a function of doping and how doped holes impact magnetic order.« less

THR-TH: a high-temperature gas-cooled nuclear reactor core thermal hydraulics code

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

Vondy, D.R.

1984-07-01

The ORNL version of PEBBLE, the (RZ) pebble bed thermal hydraulics code, has been extended for application to a prismatic gas cooled reactor core. The supplemental treatment is of one-dimensional coolant flow in up to a three-dimensional core description. Power density data from a neutronics and exposure calculation are used as the basic information for the thermal hydraulics calculation of heat removal. Two-dimensional neutronics results may be expanded for a three-dimensional hydraulics calculation. The geometric description for the hydraulics problem is the same as used by the neutronics code. A two-dimensional thermal cell model is used to predict temperatures inmore » the fuel channel. The capability is available in the local BOLD VENTURE computation system for reactor core analysis with capability to account for the effect of temperature feedback by nuclear cross section correlation. Some enhancements have also been added to the original code to add pebble bed modeling flexibility and to generate useful auxiliary results. For example, an estimate is made of the distribution of fuel temperatures based on average and extreme conditions regularly calculated at a number of locations.« less

Numerically exploring the 1D-2D dimensional crossover on spin dynamics in the doped Hubbard model

DOE PAGES

Kung, Y. F.; Bazin, C.; Wohlfeld, K.; ...

2017-11-02

Using determinant quantum Monte Carlo (DQMC) simulations, we systematically study the doping dependence of the crossover from one to two dimensions and its impact on the magnetic properties of the Hubbard model. A square lattice of chains is used, in which the dimensionality can be tuned by varying the interchain coupling t ⊥. The dynamical spin structure factor and static quantities, such as the static spin susceptibility and nearest-neighbor spin correlation function, are characterized in the one- and two-dimensional limits as a benchmark. When the dimensionality is tuned between these limits, the magnetic properties, while evolving smoothly from one tomore » two dimensions, drastically change regardless of the doping level. This suggests that the spin excitations in the two-dimensional Hubbard model, even in the heavily doped case, cannot be explained using the spinon picture known from one dimension. In conclusion, the DQMC calculations are complemented by cluster perturbation theory studies to form a more complete picture of how the crossover occurs as a function of doping and how doped holes impact magnetic order.« less

One-loop corrections to light cone wave functions: The dipole picture DIS cross section

NASA Astrophysics Data System (ADS)

Hänninen, H.; Lappi, T.; Paatelainen, R.

2018-06-01

We develop methods to perform loop calculations in light cone perturbation theory using a helicity basis, refining the method introduced in our earlier work. In particular this includes implementing a consistent way to contract the four-dimensional tensor structures from the helicity vectors with d-dimensional tensors arising from loop integrals, in a way that can be fully automatized. We demonstrate this explicitly by calculating the one-loop correction to the virtual photon to quark-antiquark dipole light cone wave function. This allows us to calculate the deep inelastic scattering cross section in the dipole formalism to next-to-leading order accuracy. Our results, obtained using the four dimensional helicity scheme, agree with the recent calculation by Beuf using conventional dimensional regularization, confirming the regularization scheme independence of this cross section.

Deformed quantum double realization of the toric code and beyond

NASA Astrophysics Data System (ADS)

Padmanabhan, Pramod; Ibieta-Jimenez, Juan Pablo; Bernabe Ferreira, Miguel Jorge; Teotonio-Sobrinho, Paulo

2016-09-01

Quantum double models, such as the toric code, can be constructed from transfer matrices of lattice gauge theories with discrete gauge groups and parametrized by the center of the gauge group algebra and its dual. For general choices of these parameters the transfer matrix contains operators acting on links which can also be thought of as perturbations to the quantum double model driving it out of its topological phase and destroying the exact solvability of the quantum double model. We modify these transfer matrices with perturbations and extract exactly solvable models which remain in a quantum phase, thus nullifying the effect of the perturbation. The algebra of the modified vertex and plaquette operators now obey a deformed version of the quantum double algebra. The Abelian cases are shown to be in the quantum double phase whereas the non-Abelian phases are shown to be in a modified phase of the corresponding quantum double phase. These are illustrated with the groups Zn and S3. The quantum phases are determined by studying the excitations of these systems namely their fusion rules and the statistics. We then go further to construct a transfer matrix which contains the other Z2 phase namely the double semion phase. More generally for other discrete groups these transfer matrices contain the twisted quantum double models. These transfer matrices can be thought of as being obtained by introducing extra parameters into the transfer matrix of lattice gauge theories. These parameters are central elements belonging to the tensor products of the algebra and its dual and are associated to vertices and volumes of the three dimensional lattice. As in the case of the lattice gauge theories we construct the operators creating the excitations in this case and study their braiding and fusion properties.

Verification and benchmark testing of the NUFT computer code

NASA Astrophysics Data System (ADS)

Lee, K. H.; Nitao, J. J.; Kulshrestha, A.

1993-10-01

This interim report presents results of work completed in the ongoing verification and benchmark testing of the NUFT (Nonisothermal Unsaturated-saturated Flow and Transport) computer code. NUFT is a suite of multiphase, multicomponent models for numerical solution of thermal and isothermal flow and transport in porous media, with application to subsurface contaminant transport problems. The code simulates the coupled transport of heat, fluids, and chemical components, including volatile organic compounds. Grid systems may be cartesian or cylindrical, with one-, two-, or fully three-dimensional configurations possible. In this initial phase of testing, the NUFT code was used to solve seven one-dimensional unsaturated flow and heat transfer problems. Three verification and four benchmarking problems were solved. In the verification testing, excellent agreement was observed between NUFT results and the analytical or quasianalytical solutions. In the benchmark testing, results of code intercomparison were very satisfactory. From these testing results, it is concluded that the NUFT code is ready for application to field and laboratory problems similar to those addressed here. Multidimensional problems, including those dealing with chemical transport, will be addressed in a subsequent report.

Effects of Radiation Damping in Extreme Ultra-intense Laser-Plasma Interaction

NASA Astrophysics Data System (ADS)

Pandit, Rishi R.

Recent advances in the development of intense short pulse lasers are significant. Now it is available to access a laser with intensity 1021W/cm2 by focusing a petawatt class laser. In a few years, the intensity will exceed 1022W/cm2 , at which intensity electrons accelerated by the laser get energy more than 100 MeV and start to emit radiation strongly. Resultingly, the damping of electron motion can become large. In order to study this problem, we developed a code to solve a set of equations describing the evolution of a strong electromagnetic wave interacting with a single electron. Usually the equation of motion of an electron including radiation damping under the influence of electromagnetic fields is derived from the Lorentz-Dirac equation treating the damping as a perturbation. So far people had used the first order damping equation. This is because the second order term seems to be small and actually it is negligible under 1022W/cm2 intensity. The derivation of 2nd order equation is also complicated and challenging. We derived the second order damping equations for the first time and implemented in the code. The code was then tested via single particle motion in the extreme intensity laser. It was found that the 1st order damping term is reasonable up to the intensity 1022W/cm2, but the 2nd oder term becomes not negligible and comparable in magnitude to the first order term beyond 1023W/cm2. The radiation damping model was introduced using a one-dimensional particle-in-cell code (PIC), and tested in the laser-plasma interaction at extreme intensity. The strong damping of hot electrons in high energy tail was demonstrated in PIC simulations.

The development of global GRAPES 4DVAR

NASA Astrophysics Data System (ADS)

Liu, Yongzhu

2017-04-01

Four-dimensional variation data assimilation (4DVAR) has given a great contribution to the improvement of NWP system over the past twenty years. Therefore, our strategy is to develop an operational global 4D-Var system from the outset. The aim at the paper is to introduce the development of the global GRAPES four-dimensional variation data assimilation (4DVAR) using incremental analysis schemes and to presents results of a comparison between 4DVAR using 6-hour assimilation window and simplified physics during the minimization with three-dimensional variation data assimilation (3DVAR). The dynamical cores of the tangent-linear and adjoint models are developed directly based on the non-hydrostatic forecast model. In addition, the standard correctness checks have been performed. As well as the development adjoint codes, most of our work is focused on improving the computational efficiency since the bulk of the computational cost of 4D-Var is in the integration of the tangent-linear and adjoint models. In terms of tangent-linear model, the wall-clock time is reduced to about 1.2 times as much as one of nonlinear model through the optimizing of the software framework. The significant computational cost savings on adjoint model result from the removing the redundant recompilations of model trajectories. It is encouraging that the wall-clock time of adjoint model is less than 1.5 times as much as one of nonlinear model. The current difficulty is that the numerical scheme used within the linear model is based on strategically on the numeric of the corresponding nonlinear model. Further computational acceleration should be expected from the improvement on nonlinear numerical algorithm. A series of linearized physical parameterization schemes has been developed to improve the representation of perturbed fields in the linear model. It consists of horizontal and vertical diffusion, sub-grid scale orographic gravity wave drag, large-scale condensation and cumulus convection schemes. We also found the straightforward linearization based on the nonlinear physical scheme might lead to significant growing of spurious unstable perturbations. It is essential to simplify the linear physics with respect to the non-linear schemes. The improvement on the perturbed fields in the tangent-linear model is visible with the linear physics included, especially at the low level. GRAPES variation data assimilation system adopts the incremental approach. The work is ongoing to develop a pre-operational 4DVAR suite with 0.25° outer loop resolution and multiple outer-loops configurations. One 4DVAR analysis using 6-hour assimilation windows can be finished within 40-minutes when using the available conventional and satellite data. In summary, it was found that the analysis over the northern, southern hemispheres, tropical region and East Asian area of GRAPES 4DVAR performed better than GRAPES 3DVAR for one month experiments. Moreover, the forecast results show that northern and southern extra-tropical scores for GRAPES 4DVAR are already better than GRAPES 3DVAR, but the tropical performance needs further investigations. Therefore, the subsequent main improvements will aim to enhance its computational efficiency and accuracy in 2017. The global GRAPES 4DVAR is planned for operation in 2018.

Modelling of three dimensional equilibrium and stability of MAST plasmas with magnetic perturbations using VMEC and COBRA

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

Ham, C. J., E-mail: christopher.ham@ccfe.ac.uk; Chapman, I. T.; Kirk, A.

2014-10-15

It is known that magnetic perturbations can mitigate edge localized modes (ELMs) in experiments, for example, MAST [Kirk et al., Nucl. Fusion 53, 043007 (2013)]. One hypothesis is that the magnetic perturbations cause a three dimensional corrugation of the plasma and this corrugated plasma has different stability properties to peeling-ballooning modes compared to an axisymmetric plasma. It has been shown in an up-down symmetric plasma that magnetic perturbations in tokamaks will break the usual axisymmetry of the plasma causing three dimensional displacements [Chapman et al., Plasma Phys. Controlled Fusion 54, 105013 (2012)]. We produce a free boundary three-dimensional equilibrium ofmore » a lower single null MAST relevant plasma using VMEC [S. P. Hirshman and J. C. Whitson, Phys. Fluids 26, 3553 (1983)]. The safety factor and pressure profiles used for the modelling are similar to those deduced from axisymmetric analysis of experimental data with ELMs. We focus on the effect of applying n = 3 and n = 6 magnetic perturbations using the resonant magnetic perturbation (RMP) coils. A midplane displacement of over ±1 cm is seen when the full current is applied. The current in the coils is scanned and a linear relationship between coil current and midplane displacement is found. The pressure gradient in real space in different toroidal locations is shown to change when RMPs are applied. This effect should be taken into account when diagnosing plasmas with RMPs applied. The helical Pfirsch-Schlüter currents which arise as a result of the assumption of nested flux surfaces are estimated for this equilibrium. The effect of this non-axisymmetric equilibrium on infinite n ballooning stability is investigated using COBRA [Sanchez et al., J. Comput. Phys. 161, 576–588 (2000)]. The infinite n ballooning stability is analysed for two reasons; it may give an indication of the effect of non-axisymmetry on finite n peeling-ballooning modes, responsible for ELMs; and infinite n ballooning modes are correlated to kinetic ballooning modes which are thought to limit the pressure gradient of the pedestal [Snyder et al., Phys. Plasmas 16, 056118 (2009)]. The ballooning mode growth rate gains a variation in toroidal angle. The equilibria with midplane displacements due to RMP coils have a higher ballooning mode growth rate than the axisymmetric case and the possible implications are discussed.« less

Exploration of magnetic perturbation effects on advanced divertor configurations in NSTX-U

DOE Data Explorer

Frerichs, H. [Princeton Plasma Physics Lab. (PPPL), Princeton, NJ (United States); Waters, I. [Princeton Plasma Physics Lab. (PPPL), Princeton, NJ (United States); Schmitz, O. [Princeton Plasma Physics Lab. (PPPL), Princeton, NJ (United States); Canal, G. P. [Princeton Plasma Physics Lab. (PPPL), Princeton, NJ (United States); Evans, T. E. [Princeton Plasma Physics Lab. (PPPL), Princeton, NJ (United States); Feng, Y. [Princeton Plasma Physics Lab. (PPPL), Princeton, NJ (United States); Soukhanovskii, V. A. [Princeton Plasma Physics Lab. (PPPL), Princeton, NJ (United States)

2016-06-01

The control of divertor heat loads - both steady state and transient - remains a key challenge for the successful operation of ITER and FNSF. Magnetic perturbations provide a promising technique to control ELMs (transients), but understanding their detailed impact is difficult due to their symmetry breaking nature. One approach for reducing steady state heat loads are so called 'advanced divertors' which aim at optimizing the magnetic field configuration: the snowflake and the (super-)X-divertor. It is likely that both concepts - magnetic perturbations and advanced divertors - will have to work together, and we explore their interaction based on the NSTX-U setup. An overview of different divertor configurations under the impact of magnetic perturbations is presented, and the resulting impact on plasma edge transport is investigated with the EMC3-EIRENE code. Variations in size of the magnetic footprint of the perturbed separatrix are found, which is related to the level of flux expansion on the divertor target. Non-axisymmetric peaking of the heat flux related to the perturbed separatrix is found at the outer strike point, but only in locations where flux expansion is not too large.

Pinning by rare defects and effective mobility for elastic interfaces in high dimensions

NASA Astrophysics Data System (ADS)

Cao, Xiangyu; Démery, Vincent; Rosso, Alberto

2018-06-01

The existence of a depinning transition for a high dimensional interface in a weakly disordered medium is controversial. Following Larkin arguments and a perturbative expansion, one expects a linear response with a renormalized mobility . In this paper, we compare these predictions with the exact solution of a fully connected model, which displays a finite critical force . At small disorder, we unveil an intermediary linear regime for characterized by the renormalized mobility . Our results suggest that in high dimension the critical force is always finite and determined by the effect of rare impurities that is missed by the perturbative expansion. However, the perturbative expansion correctly describes an intermediate regime that should be visible at small disorder.

Apparently noninvariant terms of nonlinear sigma models in lattice perturbation theory

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

Harada, Koji; Hattori, Nozomu; Kubo, Hirofumi

2009-03-15

Apparently noninvariant terms (ANTs) that appear in loop diagrams for nonlinear sigma models are revisited in lattice perturbation theory. The calculations have been done mostly with dimensional regularization so far. In order to establish that the existence of ANTs is independent of the regularization scheme, and of the potential ambiguities in the definition of the Jacobian of the change of integration variables from group elements to 'pion' fields, we employ lattice regularization, in which everything (including the Jacobian) is well defined. We show explicitly that lattice perturbation theory produces ANTs in the four-point functions of the pion fields at one-loopmore » and the Jacobian does not play an important role in generating ANTs.« less

Effects of two-dimensional magnetic uncertainties and three-dimensional error and perturbation fields on the Small Angle Slot divertor geometry and topology [Effects of two- and three-dimensional magnetic fields on the Small Angle Slot divertor magnetic topology

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

Trevisan, Gregorio L.; Lao, Lang L.; Evans, Todd E.

The Small Angle Slot (SAS) was recently installed on DIII-D as an advanced divertor, promising easier plasma detachment and lower temperatures across the whole target. A twofold study of the SAS magnetic topology is presented in this paper. On one hand, a twodimensional uncertainty quantification analysis is carried out through a Monte Carlo approach in order to understand the level of accuracy of two-dimensional equilibrium computations in reconstructing the strike point and angle onto the divertor. Under typical experimental conditions, the uncertainties are found to be roughly 6.8 mm and 0.56 deg, respectively. On the other hand, a three-dimensional “vacuum”more » analysis is carried out to understand the effects of typical external perturbation fields on the scrape-off layer topology. When the threedimensional I-coils are switched on, poloidally-localized lobes are found to appear, grow, and hit the SAS target, although barely, even for 5 kA; at the same time, the strike point modulation is found to be roughly 1.8 mm and thus negligible for most purposes. Furthermore, such results complement previous two-dimensional analyses in characterizing typical SAS equilibria and provide useful background information for planning and interpreting SAS experiments.« less

Effects of two-dimensional magnetic uncertainties and three-dimensional error and perturbation fields on the Small Angle Slot divertor geometry and topology [Effects of two- and three-dimensional magnetic fields on the Small Angle Slot divertor magnetic topology

DOE PAGES

Trevisan, Gregorio L.; Lao, Lang L.; Evans, Todd E.; ...

2018-01-04

The Small Angle Slot (SAS) was recently installed on DIII-D as an advanced divertor, promising easier plasma detachment and lower temperatures across the whole target. A twofold study of the SAS magnetic topology is presented in this paper. On one hand, a twodimensional uncertainty quantification analysis is carried out through a Monte Carlo approach in order to understand the level of accuracy of two-dimensional equilibrium computations in reconstructing the strike point and angle onto the divertor. Under typical experimental conditions, the uncertainties are found to be roughly 6.8 mm and 0.56 deg, respectively. On the other hand, a three-dimensional “vacuum”more » analysis is carried out to understand the effects of typical external perturbation fields on the scrape-off layer topology. When the threedimensional I-coils are switched on, poloidally-localized lobes are found to appear, grow, and hit the SAS target, although barely, even for 5 kA; at the same time, the strike point modulation is found to be roughly 1.8 mm and thus negligible for most purposes. Furthermore, such results complement previous two-dimensional analyses in characterizing typical SAS equilibria and provide useful background information for planning and interpreting SAS experiments.« less

Computational Optimization of a Natural Laminar Flow Experimental Wing Glove

NASA Technical Reports Server (NTRS)

Hartshom, Fletcher

2012-01-01

Computational optimization of a natural laminar flow experimental wing glove that is mounted on a business jet is presented and discussed. The process of designing a laminar flow wing glove starts with creating a two-dimensional optimized airfoil and then lofting it into a three-dimensional wing glove section. The airfoil design process does not consider the three dimensional flow effects such as cross flow due wing sweep as well as engine and body interference. Therefore, once an initial glove geometry is created from the airfoil, the three dimensional wing glove has to be optimized to ensure that the desired extent of laminar flow is maintained over the entire glove. TRANAIR, a non-linear full potential solver with a coupled boundary layer code was used as the main tool in the design and optimization process of the three-dimensional glove shape. The optimization process uses the Class-Shape-Transformation method to perturb the geometry with geometric constraints that allow for a 2-in clearance from the main wing. The three-dimensional glove shape was optimized with the objective of having a spanwise uniform pressure distribution that matches the optimized two-dimensional pressure distribution as closely as possible. Results show that with the appropriate inputs, the optimizer is able to match the two dimensional pressure distributions practically across the entire span of the wing glove. This allows for the experiment to have a much higher probability of having a large extent of natural laminar flow in flight.

Cosmological perturbations in the (1 + 3 + 6)-dimensional space-times

NASA Astrophysics Data System (ADS)

Tomita, K.

2014-12-01

Cosmological perturbations in the (1+3+6)-dimensional space-times including photon gas without viscous processes are studied on the basis of Abbott et al.'s formalism [R. B. Abbott, B. Bednarz, and S. D. Ellis, Phys. Rev. D 33, 2147 (1986)]. Space-times consist of outer space (the 3-dimensional expanding section) and inner space (the 6-dimensional section). The inner space expands initially and later contracts. Abbott et al. derived only power-type solutions, which appear at the final stage of the space-times, in the small wave-number limit. In this paper, we derive not only small wave-number solutions, but also large wave-number solutions. It is found that the latter solutions depend on the two wave-numbers k_r and k_R (which are defined in the outer and inner spaces, respectively), and that the k_r-dependent and k_R-dependent parts dominate the total perturbations when (k_r/r(t))/(k_R/R(t)) ≫ 1 or ≪ 1, respectively, where r(t) and R(t) are the scale-factors in the outer and inner spaces. By comparing the behaviors of these perturbations, moreover, changes in the spectrum of perturbations in the outer space with time are discussed.

Users manual for the NASA Lewis three-dimensional ice accretion code (LEWICE 3D)

NASA Technical Reports Server (NTRS)

Bidwell, Colin S.; Potapczuk, Mark G.

1993-01-01

A description of the methodology, the algorithms, and the input and output data along with an example case for the NASA Lewis 3D ice accretion code (LEWICE3D) has been produced. The manual has been designed to help the user understand the capabilities, the methodologies, and the use of the code. The LEWICE3D code is a conglomeration of several codes for the purpose of calculating ice shapes on three-dimensional external surfaces. A three-dimensional external flow panel code is incorporated which has the capability of calculating flow about arbitrary 3D lifting and nonlifting bodies with external flow. A fourth order Runge-Kutta integration scheme is used to calculate arbitrary streamlines. An Adams type predictor-corrector trajectory integration scheme has been included to calculate arbitrary trajectories. Schemes for calculating tangent trajectories, collection efficiencies, and concentration factors for arbitrary regions of interest for single droplets or droplet distributions have been incorporated. A LEWICE 2D based heat transfer algorithm can be used to calculate ice accretions along surface streamlines. A geometry modification scheme is incorporated which calculates the new geometry based on the ice accretions generated at each section of interest. The three-dimensional ice accretion calculation is based on the LEWICE 2D calculation. Both codes calculate the flow, pressure distribution, and collection efficiency distribution along surface streamlines. For both codes the heat transfer calculation is divided into two regions, one above the stagnation point and one below the stagnation point, and solved for each region assuming a flat plate with pressure distribution. Water is assumed to follow the surface streamlines, hence starting at the stagnation zone any water that is not frozen out at a control volume is assumed to run back into the next control volume. After the amount of frozen water at each control volume has been calculated the geometry is modified by adding the ice at each control volume in the surface normal direction.

Superimposed Code Theorectic Analysis of DNA Codes and DNA Computing

DTIC Science & Technology

2010-03-01

because only certain collections (partitioned by font type) of sequences are allowed to be in each position (e.g., Arial = position 0, Comic ...rigidity of short oligos and the shape of the polar charge. Oligo movement was modeled by a Brownian motion 3 dimensional random walk. The one...temperature, kB is Boltz he viscosity of the medium. The random walk motion is modeled by assuming the oligo is on a three dimensional lattice and may

COMMIX-PPC: A three-dimensional transient multicomponent computer program for analyzing performance of power plant condensers

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

Chien, T.H.; Domanus, H.M.; Sha, W.T.

1993-02-01

The COMMIX-PPC computer pregrain is an extended and improved version of earlier COMMIX codes and is specifically designed for evaluating the thermal performance of power plant condensers. The COMMIX codes are general-purpose computer programs for the analysis of fluid flow and heat transfer in complex Industrial systems. In COMMIX-PPC, two major features have been added to previously published COMMIX codes. One feature is the incorporation of one-dimensional equations of conservation of mass, momentum, and energy on the tube stile and the proper accounting for the thermal interaction between shell and tube side through the porous-medium approach. The other added featuremore » is the extension of the three-dimensional conservation equations for shell-side flow to treat the flow of a multicomponent medium. COMMIX-PPC is designed to perform steady-state and transient. Three-dimensional analysis of fluid flow with heat transfer tn a power plant condenser. However, the code is designed in a generalized fashion so that, with some modification, it can be used to analyze processes in any heat exchanger or other single-phase engineering applications. Volume I (Equations and Numerics) of this report describes in detail the basic equations, formulation, solution procedures, and models for a phenomena. Volume II (User's Guide and Manual) contains the input instruction, flow charts, sample problems, and descriptions of available options and boundary conditions.« less

An exact solution of solute transport by one-dimensional random velocity fields

USGS Publications Warehouse

Cvetkovic, V.D.; Dagan, G.; Shapiro, A.M.

1991-01-01

The problem of one-dimensional transport of passive solute by a random steady velocity field is investigated. This problem is representative of solute movement in porous media, for example, in vertical flow through a horizontally stratified formation of variable porosity with a constant flux at the soil surface. Relating moments of particle travel time and displacement, exact expressions for the advection and dispersion coefficients in the Focker-Planck equation are compared with the perturbation results for large distances. The first- and second-order approximations for the dispersion coefficient are robust for a lognormal velocity field. The mean Lagrangian velocity is the harmonic mean of the Eulerian velocity for large distances. This is an artifact of one-dimensional flow where the continuity equation provides for a divergence free fluid flux, rather than a divergence free fluid velocity. ?? 1991 Springer-Verlag.

Deep linear autoencoder and patch clustering-based unified one-dimensional coding of image and video

NASA Astrophysics Data System (ADS)

Li, Honggui

2017-09-01

This paper proposes a unified one-dimensional (1-D) coding framework of image and video, which depends on deep learning neural network and image patch clustering. First, an improved K-means clustering algorithm for image patches is employed to obtain the compact inputs of deep artificial neural network. Second, for the purpose of best reconstructing original image patches, deep linear autoencoder (DLA), a linear version of the classical deep nonlinear autoencoder, is introduced to achieve the 1-D representation of image blocks. Under the circumstances of 1-D representation, DLA is capable of attaining zero reconstruction error, which is impossible for the classical nonlinear dimensionality reduction methods. Third, a unified 1-D coding infrastructure for image, intraframe, interframe, multiview video, three-dimensional (3-D) video, and multiview 3-D video is built by incorporating different categories of videos into the inputs of patch clustering algorithm. Finally, it is shown in the results of simulation experiments that the proposed methods can simultaneously gain higher compression ratio and peak signal-to-noise ratio than those of the state-of-the-art methods in the situation of low bitrate transmission.

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.

Effects of Nonlinear Inhomogeneity on the Cosmic Expansion with Numerical Relativity.

PubMed

Bentivegna, Eloisa; Bruni, Marco

2016-06-24

We construct a three-dimensional, fully relativistic numerical model of a universe filled with an inhomogeneous pressureless fluid, starting from initial data that represent a perturbation of the Einstein-de Sitter model. We then measure the departure of the average expansion rate with respect to this homogeneous and isotropic reference model, comparing local quantities to the predictions of linear perturbation theory. We find that collapsing perturbations reach the turnaround point much earlier than expected from the reference spherical top-hat collapse model and that the local deviation of the expansion rate from the homogeneous one can be as high as 28% at an underdensity, for an initial density contrast of 10^{-2}. We then study, for the first time, the exact behavior of the backreaction term Q_{D}. We find that, for small values of the initial perturbations, this term exhibits a 1/a scaling, and that it is negative with a linearly growing absolute value for larger perturbation amplitudes, thereby contributing to an overall deceleration of the expansion. Its magnitude, on the other hand, remains very small even for relatively large perturbations.

Behaviour of a Premixed Flame Subjected to Acoustic Oscillations

PubMed Central

Qureshi, Shafiq R.; Khan, Waqar A.; Prosser, Robert

2013-01-01

In this paper, a one dimensional premixed laminar methane flame is subjected to acoustic oscillations and studied. The purpose of this analysis is to investigate the effects of acoustic perturbations on the reaction rates of different species, with a view to their respective contribution to thermoacoustic instabilities. Acoustically transparent non reflecting boundary conditions are employed. The flame response has been studied with acoustic waves of different frequencies and amplitudes. The integral values of the reaction rates, the burning velocities and the heat release of the acoustically perturbed flame are compared with the unperturbed case. We found that the flame's sensitivity to acoustic perturbations is greatest when the wavelength is comparable to the flame thickness. Even in this case, the perturbations are stable with time. We conclude that acoustic fields acting on the chemistry do not contribute significantly to the emergence of large amplitude pressure oscillations. PMID:24376501

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.

A new method for computing the gyrocenter orbit in the tokamak configuration

NASA Astrophysics Data System (ADS)

Xu, Yingfeng

2013-10-01

Gyrokinetic theory is an important tool for studying the long-time behavior of magnetized plasmas in Tokamaks. The gyrocenter trajectory determined by the gyrocenter equations of motion can be computed by using a special kind of the Lie-transform perturbation method. The corresponding Lie-transform called I-transform makes that the transformed equations of motion have the same form as the unperturbed ones. The gyrocenter trajectory in short time is divided into two parts. One is along the unperturbed orbit. The other one, which is related to perturbation, is determined by the I-transform generating vector. The numerical gyrocenter orbit code based on this new method has been developed in the tokamak configuration and benchmarked with the other orbit code in some simple cases. Furthermore, it is clearly demonstrated that this new method for computing gyrocenter orbit is equivalent to the gyrocenter Hamilton equations of motion up to the second order in timestep. The new method can be applied to the gyrokinetic simulation. The gyrocenter orbit of the unperturbed part determined by the equilibrium fields can be computed previously in the gyrokinetic simulation, and the corresponding time consumption is neglectable.

Numerical simulations of Richtmyer{endash}Meshkov instabilities in finite-thickness fluid layers

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

Mikaelian, K.O.

1996-05-01

Direct numerical simulations of Richtmyer{endash}Meshkov instabilities in shocked fluid layers are reported and compared with analytic theory. To investigate new phenomena such as freeze-out, interface coupling, and feedthrough, several new configurations are simulated on a two-dimensional hydrocode. The basic system is an {ital A}/{ital B}/{ital A} combination, where {ital A} is air and {ital B} is a finite-thickness layer of freon, SF{sub 6}, or helium. The middle layer {ital B} has perturbations either on its upstream or downstream side, or on both sides, in which case the perturbations may be in phase (sinuous) or out of phase (varicose). The evolutionmore » of such perturbations under a Mach 1.5 shock is calculated, including the effect of a reshock. Recently reported gas curtain experiments [J. M. Budzinski {ital et} {ital al}., Phys. Fluids {bold 6}, 3510 (1994)] are also simulated and the code results are found to agree very well with the experiments. A new gas curtain configuration is also considered, involving an initially sinuous SF{sub 6} or helium layer and a new pattern, opposite mushrooms, is predicted to emerge. Upon reshock a relatively simple sinuous gas curtain is found to evolve into a highly complex pattern of nested mushrooms. {copyright} {ital 1996 American Institute of Physics.}« less

QSL Squasher: A Fast Quasi-separatrix Layer Map Calculator

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

Tassev, Svetlin; Savcheva, Antonia, E-mail: svetlin.tassev@cfa.harvard.edu

Quasi-Separatrix Layers (QSLs) are a useful proxy for the locations where current sheets can develop in the solar corona, and give valuable information about the connectivity in complicated magnetic field configurations. However, calculating QSL maps, even for two-dimensional slices through three-dimensional models of coronal magnetic fields, is a non-trivial task, as it usually involves tracing out millions of magnetic field lines with immense precision. Thus, extending QSL calculations to three dimensions has rarely been done until now. In order to address this challenge, we present QSL Squasher—a public, open-source code, which is optimized for calculating QSL maps in both twomore » and three dimensions on graphics processing units. The code achieves large processing speeds for three reasons, each of which results in an order-of-magnitude speed-up. (1) The code is parallelized using OpenCL. (2) The precision requirements for the QSL calculation are drastically reduced by using perturbation theory. (3) A new boundary detection criterion between quasi-connectivity domains is used, which quickly identifies possible QSL locations that need to be finely sampled by the code. That boundary detection criterion relies on finding the locations of abrupt field-line length changes, which we do by introducing a new Field-line Length Edge (FLEDGE) map. We find FLEDGE maps useful on their own as a quick-and-dirty substitute for QSL maps. QSL Squasher allows construction of high-resolution 3D FLEDGE maps in a matter of minutes, which is two orders of magnitude faster than calculating the corresponding 3D QSL maps. We include a sample of calculations done using QSL Squasher to demonstrate its capabilities as a QSL calculator, as well as to compare QSL and FLEDGE maps.« less

Quasi 1D Modeling of Mixed Compression Supersonic Inlets

NASA Technical Reports Server (NTRS)

Kopasakis, George; Connolly, Joseph W.; Paxson, Daniel E.; Woolwine, Kyle J.

2012-01-01

The AeroServoElasticity task under the NASA Supersonics Project is developing dynamic models of the propulsion system and the vehicle in order to conduct research for integrated vehicle dynamic performance. As part of this effort, a nonlinear quasi 1-dimensional model of the 2-dimensional bifurcated mixed compression supersonic inlet is being developed. The model utilizes computational fluid dynamics for both the supersonic and subsonic diffusers. The oblique shocks are modeled utilizing compressible flow equations. This model also implements variable geometry required to control the normal shock position. The model is flexible and can also be utilized to simulate other mixed compression supersonic inlet designs. The model was validated both in time and in the frequency domain against the legacy LArge Perturbation INlet code, which has been previously verified using test data. This legacy code written in FORTRAN is quite extensive and complex in terms of the amount of software and number of subroutines. Further, the legacy code is not suitable for closed loop feedback controls design, and the simulation environment is not amenable to systems integration. Therefore, a solution is to develop an innovative, more simplified, mixed compression inlet model with the same steady state and dynamic performance as the legacy code that also can be used for controls design. The new nonlinear dynamic model is implemented in MATLAB Simulink. This environment allows easier development of linear models for controls design for shock positioning. The new model is also well suited for integration with a propulsion system model to study inlet/propulsion system performance, and integration with an aero-servo-elastic system model to study integrated vehicle ride quality, vehicle stability, and efficiency.

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

Bovy, Jo, E-mail: bovy@ias.edu

I describe the design, implementation, and usage of galpy, a python package for galactic-dynamics calculations. At its core, galpy consists of a general framework for representing galactic potentials both in python and in C (for accelerated computations); galpy functions, objects, and methods can generally take arbitrary combinations of these as arguments. Numerical orbit integration is supported with a variety of Runge-Kutta-type and symplectic integrators. For planar orbits, integration of the phase-space volume is also possible. galpy supports the calculation of action-angle coordinates and orbital frequencies for a given phase-space point for general spherical potentials, using state-of-the-art numerical approximations for axisymmetricmore » potentials, and making use of a recent general approximation for any static potential. A number of different distribution functions (DFs) are also included in the current release; currently, these consist of two-dimensional axisymmetric and non-axisymmetric disk DFs, a three-dimensional disk DF, and a DF framework for tidal streams. I provide several examples to illustrate the use of the code. I present a simple model for the Milky Way's gravitational potential consistent with the latest observations. I also numerically calculate the Oort functions for different tracer populations of stars and compare them to a new analytical approximation. Additionally, I characterize the response of a kinematically warm disk to an elliptical m = 2 perturbation in detail. Overall, galpy consists of about 54,000 lines, including 23,000 lines of code in the module, 11,000 lines of test code, and about 20,000 lines of documentation. The test suite covers 99.6% of the code. galpy is available at http://github.com/jobovy/galpy with extensive documentation available at http://galpy.readthedocs.org/en/latest.« less

QSL Squasher: A Fast Quasi-separatrix Layer Map Calculator

NASA Astrophysics Data System (ADS)

Tassev, Svetlin; Savcheva, Antonia

2017-05-01

Quasi-Separatrix Layers (QSLs) are a useful proxy for the locations where current sheets can develop in the solar corona, and give valuable information about the connectivity in complicated magnetic field configurations. However, calculating QSL maps, even for two-dimensional slices through three-dimensional models of coronal magnetic fields, is a non-trivial task, as it usually involves tracing out millions of magnetic field lines with immense precision. Thus, extending QSL calculations to three dimensions has rarely been done until now. In order to address this challenge, we present QSL Squasher—a public, open-source code, which is optimized for calculating QSL maps in both two and three dimensions on graphics processing units. The code achieves large processing speeds for three reasons, each of which results in an order-of-magnitude speed-up. (1) The code is parallelized using OpenCL. (2) The precision requirements for the QSL calculation are drastically reduced by using perturbation theory. (3) A new boundary detection criterion between quasi-connectivity domains is used, which quickly identifies possible QSL locations that need to be finely sampled by the code. That boundary detection criterion relies on finding the locations of abrupt field-line length changes, which we do by introducing a new Field-line Length Edge (FLEDGE) map. We find FLEDGE maps useful on their own as a quick-and-dirty substitute for QSL maps. QSL Squasher allows construction of high-resolution 3D FLEDGE maps in a matter of minutes, which is two orders of magnitude faster than calculating the corresponding 3D QSL maps. We include a sample of calculations done using QSL Squasher to demonstrate its capabilities as a QSL calculator, as well as to compare QSL and FLEDGE maps.

Spectral calculations for pressure-velocity and pressure-strain correlations in homogeneous shear turbulence

NASA Astrophysics Data System (ADS)

Dutta, Kishore

2018-02-01

Theoretical analyses of pressure related turbulent statistics are vital for a reliable and accurate modeling of turbulence. In the inertial subrange of turbulent shear flow, pressure-velocity and pressure-strain correlations are affected by anisotropy imposed at large scales. Recently, Tsuji and Kaneda (2012 J. Fluid Mech. 694 50) performed a set of experiments on homogeneous shear flow, and estimated various one-dimensional pressure related spectra and the associated non-dimensional universal numbers. Here, starting from the governing Navier-Stokes dynamics for the fluctuating velocity field and assuming the anisotropy at inertial scales as a weak perturbation of an otherwise isotropic dynamics, we analytically derive the form of the pressure-velocity and pressure-strain correlations. The associated universal numbers are calculated using the well-known renormalization-group results, and are compared with the experimental estimates of Tsuji and Kaneda. Approximations involved in the perturbative calculations are discussed.

Spin(7) compactifications and 1/4-BPS vacua in heterotic supergravity

NASA Astrophysics Data System (ADS)

Angus, Stephen; Matti, Cyril; Svanes, Eirik E.

2016-03-01

We continue the investigation into non-maximally symmetric compactifications of the heterotic string. In particular, we consider compactifications where the internal space is allowed to depend on two or more external directions. For preservation of supersymmetry, this implies that the internal space must in general be that of a Spin(7) manifold, which leads to a 1/4-BPS four-dimensional supersymmetric perturbative vacuum breaking all but one supercharge. We find that these solutions allow for internal geometries previously excluded by the domain-wall-type solutions, and hence the resulting four-dimensional superpotential is more generic. In particular, we find an interesting resemblance to the superpotentials that appear in non-geometric flux compactifications of type II string theory. If the vacua are to be used for phenomenological applications, they must be lifted to maximal symmetry by some non-perturbative or higher-order effect.

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.

Jet-A fuel evaporation analysis in conical tube injectors

NASA Technical Reports Server (NTRS)

Lai, M.-C.; Chue, T.-H.; Zhu, G.; Sun, H.; Tacina, R.; Chun, K.; Hicks, Y.

1991-01-01

A simple one-dimensional drop-life-history analysis and a multidimensional spray calculation using KIVA-II code are applied to the vaporization of Jet-A fuel in multiple tube injectors. Within the assumptions of the analysis, the one-dimensional results are useful for design purposes. The pressure-atomizer breakup models do not accurately predict the dropsize measured experimentally or deduced from the one-dimensional analysis. Cold flow visualization and dropsize measurements show that capillary wave breakup mechanism plays an important role in the spray angle and droplet impingement on the tube wall.

Superimposed Code Theoretic Analysis of Deoxyribonucleic Acid (DNA) Codes and DNA Computing

DTIC Science & Technology

2010-01-01

partitioned by font type) of sequences are allowed to be in each position (e.g., Arial = position 0, Comic = position 1, etc. ) and within each collection...movement was modeled by a Brownian motion 3 dimensional random walk. The one dimensional diffusion coefficient D for the ellipsoid shape with 3...temperature, kB is Boltzmann’s constant, and η is the viscosity of the medium. The random walk motion is modeled by assuming the oligo is on a three

Subscale Fast Cookoff Testing and Modeling for the Hazard Assessment of Large Rocket Motors

DTIC Science & Technology

2001-03-01

41 LIST OF TABLES Table 1 Heats of Vaporization Parameter for Two-liner Phase Transformation - Complete Liner Sublimation and/or Combined Liner...One-dimensional 2-D Two-dimensional ALE3D Arbitrary-Lagrange-Eulerian (3-D) Computer Code ALEGRA 3-D Arbitrary-Lagrange-Eulerian Computer Code for...case-liner bond areas and in the grain inner bore to explore the pre-ignition and ignition phases , as well as burning evolution in rocket motor fast

Effects of Normal and Perturbed Social Play on the Duration and Amplitude of Different Types of Infant Smiles

ERIC Educational Resources Information Center

Fogel, Alan; Hsu, Hui-Chin; Shapiro, Alyson F.; Nelson-Goens, G. Christina; Secrist, Cory

2006-01-01

Different types of smiling varying in amplitude of lip corner retraction were investigated during 2 mother-infant games--peekaboo and tickle--at 6 and 12 months and during normally occurring and perturbed games. Using Facial Action Coding System (FACS), infant smiles were coded as simple (lip corner retraction only), Duchenne (simple plus cheek…

Inflation in Flatland

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

Hinterbichler, Kurt; Joyce, Austin; Khoury, Justin, E-mail: kurt.hinterbichler@case.edu, E-mail: austin.joyce@columbia.edu, E-mail: jkhoury@sas.upenn.edu

We investigate the symmetry structure of inflation in 2+1 dimensions. In particular, we show that the asymptotic symmetries of three-dimensional de Sitter space are in one-to-one correspondence with cosmological adiabatic modes for the curvature perturbation. In 2+1 dimensions, the asymptotic symmetry algebra is infinite-dimensional, given by two copies of the Virasoro algebra, and can be traced to the conformal symmetries of the two-dimensional spatial slices of de Sitter. We study the consequences of this infinite-dimensional symmetry for inflationary correlation functions, finding new soft theorems that hold only in 2+1 dimensions. Expanding the correlation functions as a power series in themore » soft momentum q , these relations constrain the traceless part of the tensorial coefficient at each order in q in terms of a lower-point function. As a check, we verify that the O( q {sup 2}) identity is satisfied by inflationary correlation functions in the limit of small sound speed.« less

Calculation of three-dimensional compressible laminar and turbulent boundary flows. Three-dimensional compressible boundary layers of reacting gases over realistic configurations

NASA Technical Reports Server (NTRS)

Kendall, R. M.; Bonnett, W. S.; Nardo, C. T.; Abbett, M. J.

1975-01-01

A three-dimensional boundary-layer code was developed for particular application to realistic hypersonic aircraft. It is very general and can be applied to a wide variety of boundary-layer flows. Laminar, transitional, and fully turbulent flows of compressible, reacting gases are efficiently calculated by use of the code. A body-oriented orthogonal coordinate system is used for the calculation and the user has complete freedom in specifying the coordinate system within the restrictions that one coordinate must be normal to the surface and the three coordinates must be mutually orthogonal.

Hypergeometric continuation of divergent perturbation series: I. Critical exponents of the Bose-Hubbard model

NASA Astrophysics Data System (ADS)

Sanders, Sören; Holthaus, Martin

2017-10-01

We study the connection between the exponent of the order parameter of the Mott insulator-to-superfluid transition occurring in the two-dimensional Bose-Hubbard model, and the divergence exponents of its one- and two-particle correlation functions. We find that at the multicritical points all divergence exponents are related to each other, allowing us to express the critical exponent in terms of one single divergence exponent. This approach correctly reproduces the critical exponent of the three-dimensional XY universality class. Because divergence exponents can be computed in an efficient manner by hypergeometric analytic continuation, our strategy is applicable to a wide class of systems.

On the relationship between the greenhouse effect, atmospheric photochemistry, and species distribution

NASA Technical Reports Server (NTRS)

Callis, L. B.; Boughner, R. E.; Natarajan, M.

1983-01-01

The coupling that exists between infrared opacity changes and tropospheric (and to a lesser extent stratospheric) chemistry is explored in considerable detail, and the effects arising from various perturbations are examined. The studies are carried out with a fully coupled one-dimensional radiative-convective-photochemical model (RCP) that extends from the surface to 53.5 km and has the capability of calculating surface temperature changes due to both chemical and radiative perturbations. The model encompasses contemporary atmospheric chemistry and photochemistry involving the O(x), HO(x), NO(x), and Cl(x) species.

The Dalgarno-Lewis summation technique: Some comments and examples

NASA Astrophysics Data System (ADS)

Mavromatis, Harry A.

1991-08-01

The Dalgarno-Lewis technique [A. Dalgarno and J. T. Lewis, ``The exact calculation of long-range forces between atoms by perturbation theory,'' Proc. R. Soc. London Ser. A 233, 70-74 (1955)] provides an elegant method to obtain exact results for various orders in perturbation theory, while avoiding the infinite sums which arise in each order. In the present paper this technique, which perhaps has not been exploited as much as it could be, is first reviewed with attention to some of its not-so-straightforward details, and then six examples of the method are given using three different one-dimensional bases.

Effects of non-Hermitian perturbations on Weyl Hamiltonians with arbitrary topological charges

NASA Astrophysics Data System (ADS)

Cerjan, Alexander; Xiao, Meng; Yuan, Luqi; Fan, Shanhui

2018-02-01

We provide a systematic study of non-Hermitian topologically charged systems. Starting from a Hermitian Hamiltonian supporting Weyl points with arbitrary topological charge, adding a non-Hermitian perturbation transforms the Weyl points to one-dimensional exceptional contours. We analytically prove that the topological charge is preserved on the exceptional contours. In contrast to Hermitian systems, the addition of gain and loss allows for a new class of topological phase transition: when two oppositely charged exceptional contours touch, the topological charge can dissipate without opening a gap. These effects can be demonstrated in realistic photonics and acoustics systems.

First moments of nucleon generalized parton distributions

DOE PAGES

Wang, P.; Thomas, A. W.

2010-06-01

We extrapolate the first moments of the generalized parton distributions using heavy baryon chiral perturbation theory. The calculation is based on the one loop level with the finite range regularization. The description of the lattice data is satisfactory, and the extrapolated moments at physical pion mass are consistent with the results obtained with dimensional regularization, although the extrapolation in the momentum transfer to t=0 does show sensitivity to form factor effects, which lie outside the realm of chiral perturbation theory. We discuss the significance of the results in the light of modern experiments as well as QCD inspired models.

Brane-World Gravity.

PubMed

Maartens, Roy; Koyama, Kazuya

2010-01-01

The observable universe could be a 1+3-surface (the "brane") embedded in a 1+3+ d -dimensional spacetime (the "bulk"), with Standard Model particles and fields trapped on the brane while gravity is free to access the bulk. At least one of the d extra spatial dimensions could be very large relative to the Planck scale, which lowers the fundamental gravity scale, possibly even down to the electroweak (∼ TeV) level. This revolutionary picture arises in the framework of recent developments in M theory. The 1+10-dimensional M theory encompasses the known 1+9-dimensional superstring theories, and is widely considered to be a promising potential route to quantum gravity. At low energies, gravity is localized at the brane and general relativity is recovered, but at high energies gravity "leaks" into the bulk, behaving in a truly higher-dimensional way. This introduces significant changes to gravitational dynamics and perturbations, with interesting and potentially testable implications for high-energy astrophysics, black holes, and cosmology. Brane-world models offer a phenomenological way to test some of the novel predictions and corrections to general relativity that are implied by M theory. This review analyzes the geometry, dynamics and perturbations of simple brane-world models for cosmology and astrophysics, mainly focusing on warped 5-dimensional brane-worlds based on the Randall-Sundrum models. We also cover the simplest brane-world models in which 4-dimensional gravity on the brane is modified at low energies - the 5-dimensional Dvali-Gabadadze-Porrati models. Then we discuss co-dimension two branes in 6-dimensional models.

Comparison of theoretical and flight-measured local flow aerodynamics for a low-aspect-ratio fin

NASA Technical Reports Server (NTRS)

Johnson, J. B.; Sandlin, D. R.

1984-01-01

Flight test and theoretical aerodynamic data were obtained for a flight test fixture mounted on the underside of an F-104G aircraft. The theoretical data were generated using two codes, a two dimensional transonic code called Code H, and a three dimensional subsonic and supersonic code call wing-body. Pressure distributions generated by the codes for the flight test fixture as well as boundary layer displacement thickness generated by the two dimensional code were compared to the flight test data. The two dimensional code pressure distributions compared well except at the minimum pressure point and trailing edge. Shock locations compared well except at high transonic speeds. The three dimensional code pressure distributions compared well except at the trailing edge of the flight test fixture. The two dimensional code does not predict displacement thickness of the flight test fixture well.

A mechanism for hot-spot generation in a reactive two-dimensional sheared viscous layer

NASA Astrophysics Data System (ADS)

Timms, Robert; Purvis, Richard; Curtis, John P.

2018-05-01

A two-dimensional model for the non-uniform melting of a thin sheared viscous layer is developed. An asymptotic solution is presented for both a non-reactive and a reactive material. It is shown that the melt front is linearly stable to small perturbations in the non-reactive case, but becomes linearly unstable upon introduction of an Arrhenius source term to model the chemical reaction. Results demonstrate that non-uniform melting acts as a mechanism to generate hot spots that are found to be sufficient to reduce the time to ignition when compared with the corresponding one-dimensional model of melting.

Validation of extended magnetohydrodynamic simulations of the HIT-SI3 experiment using the NIMROD code

NASA Astrophysics Data System (ADS)

Morgan, K. D.; Jarboe, T. R.; Hossack, A. C.; Chandra, R. N.; Everson, C. J.

2017-12-01

The HIT-SI3 experiment uses a set of inductively driven helicity injectors to apply a non-axisymmetric current drive on the edge of the plasma, driving an axisymmetric spheromak equilibrium in a central confinement volume. These helicity injectors drive a non-axisymmetric perturbation that oscillates in time, with relative temporal phasing of the injectors modifying the mode structure of the applied perturbation. A set of three experimental discharges with different perturbation spectra are modelled using the NIMROD extended magnetohydrodynamics code, and comparisons are made to both magnetic and fluid measurements. These models successfully capture the bulk dynamics of both the perturbation and the equilibrium, though disagreements related to the pressure gradients experimentally measured exist.

EFDC1D - A ONE DIMENSIONAL HYDRODYNAMIC AND SEDIMENT TRANSPORT MODEL FOR RIVER AND STREAM NETWORKS: MODEL THEORY AND USERS GUIDE

EPA Science Inventory

This technical report describes the new one-dimensional (1D) hydrodynamic and sediment transport model EFDC1D. This model that can be applied to stream networks. The model code and two sample data sets are included on the distribution CD. EFDC1D can simulate bi-directional unstea...

Combination of large and small basis sets in electronic structure calculations on large systems

NASA Astrophysics Data System (ADS)

Røeggen, Inge; Gao, Bin

2018-04-01

Two basis sets—a large and a small one—are associated with each nucleus of the system. Each atom has its own separate one-electron basis comprising the large basis set of the atom in question and the small basis sets for the partner atoms in the complex. The perturbed atoms in molecules and solids model is at core of the approach since it allows for the definition of perturbed atoms in a system. It is argued that this basis set approach should be particularly useful for periodic systems. Test calculations are performed on one-dimensional arrays of H and Li atoms. The ground-state energy per atom in the linear H array is determined versus bond length.

Coupled Kinetic-MHD Simulations of Divertor Heat Load with ELM Perturbations

NASA Astrophysics Data System (ADS)

Cummings, Julian; Chang, C. S.; Park, Gunyoung; Sugiyama, Linda; Pankin, Alexei; Klasky, Scott; Podhorszki, Norbert; Docan, Ciprian; Parashar, Manish

2010-11-01

The effect of Type-I ELM activity on divertor plate heat load is a key component of the DOE OFES Joint Research Target milestones for this year. In this talk, we present simulations of kinetic edge physics, ELM activity, and the associated divertor heat loads in which we couple the discrete guiding-center neoclassical transport code XGC0 with the nonlinear extended MHD code M3D using the End-to-end Framework for Fusion Integrated Simulations, or EFFIS. In these coupled simulations, the kinetic code and the MHD code run concurrently on the same massively parallel platform and periodic data exchanges are performed using a memory-to-memory coupling technology provided by EFFIS. The M3D code models the fast ELM event and sends frequent updates of the magnetic field perturbations and electrostatic potential to XGC0, which in turn tracks particle dynamics under the influence of these perturbations and collects divertor particle and energy flux statistics. We describe here how EFFIS technologies facilitate these coupled simulations and discuss results for DIII-D, NSTX and Alcator C-Mod tokamak discharges.

Spin-triplet superconductivity in a weak-coupling Hubbard model for the quasi-one-dimensional compound Li0.9Mo6O17

NASA Astrophysics Data System (ADS)

Cho, Weejee; Platt, Christian; McKenzie, Ross H.; Raghu, Srinivas

2015-10-01

The purple bronze Li0.9Mo6O17 is of interest due to its quasi-one-dimensional electronic structure and the possible Luttinger liquid behavior resulting from it. For sufficiently low temperatures, it is a superconductor with a pairing symmetry that is still to be determined. To shed light on this issue, we analyze a minimal Hubbard model for this material involving four molybdenum orbitals per unit cell near quarter filling, using asymptotically exact perturbative renormalization group methods. We find that spin-triplet odd-parity superconductivity is the dominant instability. Approximate nesting properties of the two quasi-one-dimensional Fermi surfaces enhance certain second-order processes, which play crucial roles in determining the structure of the pairing gap. Notably, we find that the gap has more sign changes than required by the point-group symmetry.

Spin triplet superconductivity in a weak-coupling Hubbard model for the quasi-one-dimensional compound Li0.9 Mo6 O17

NASA Astrophysics Data System (ADS)

Platt, Christian; Cho, Weejee; McKenzie, Ross H.; Raghu, Sri

The purple bronze Li0.9Mo6O17 is of interest due to its quasi-one-dimensional electronic structure and the possible Luttinger liquid behavior resulting from it. For sufficiently low temperatures, it is a superconductor with a pairing symmetry that is still to be determined. To shed light on this issue, we analyze a minimal Hubbard model for this material involving four Molybdenum orbitals per unit cell near quarter filling, using asymptotically exact perturbative renormalization group methods. We find that spin triplet odd-parity superconductivity is the dominant instability. Approximate nesting properties of the two quasi-one-dimensional Fermi surfaces enhance certain second-order processes, which play crucial roles in determining the structure of the pairing gap. Notably, we find that the gap has accidental nodes, i.e. it has more sign changes than required by the point-group symmetry.

Dexter - A one-dimensional code for calculating thermionic performance of long converters.

NASA Technical Reports Server (NTRS)

Sawyer, C. D.

1971-01-01

This paper describes a versatile code for computing the coupled thermionic electric-thermal performance of long thermionic converters in which the temperature and voltage variations cannot be neglected. The code is capable of accounting for a variety of external electrical connection schemes, coolant flow paths and converter failures by partial shorting. Example problem solutions are given.

Structural-change localization and monitoring through a perturbation-based inverse problem.

PubMed

Roux, Philippe; Guéguen, Philippe; Baillet, Laurent; Hamze, Alaa

2014-11-01

Structural-change detection and characterization, or structural-health monitoring, is generally based on modal analysis, for detection, localization, and quantification of changes in structure. Classical methods combine both variations in frequencies and mode shapes, which require accurate and spatially distributed measurements. In this study, the detection and localization of a local perturbation are assessed by analysis of frequency changes (in the fundamental mode and overtones) that are combined with a perturbation-based linear inverse method and a deconvolution process. This perturbation method is applied first to a bending beam with the change considered as a local perturbation of the Young's modulus, using a one-dimensional finite-element model for modal analysis. Localization is successful, even for extended and multiple changes. In a second step, the method is numerically tested under ambient-noise vibration from the beam support with local changes that are shifted step by step along the beam. The frequency values are revealed using the random decrement technique that is applied to the time-evolving vibrations recorded by one sensor at the free extremity of the beam. Finally, the inversion method is experimentally demonstrated at the laboratory scale with data recorded at the free end of a Plexiglas beam attached to a metallic support.

Discrete Cosine Transform Image Coding With Sliding Block Codes

NASA Astrophysics Data System (ADS)

Divakaran, Ajay; Pearlman, William A.

1989-11-01

A transform trellis coding scheme for images is presented. A two dimensional discrete cosine transform is applied to the image followed by a search on a trellis structured code. This code is a sliding block code that utilizes a constrained size reproduction alphabet. The image is divided into blocks by the transform coding. The non-stationarity of the image is counteracted by grouping these blocks in clusters through a clustering algorithm, and then encoding the clusters separately. Mandela ordered sequences are formed from each cluster i.e identically indexed coefficients from each block are grouped together to form one dimensional sequences. A separate search ensues on each of these Mandela ordered sequences. Padding sequences are used to improve the trellis search fidelity. The padding sequences absorb the error caused by the building up of the trellis to full size. The simulations were carried out on a 256x256 image ('LENA'). The results are comparable to any existing scheme. The visual quality of the image is enhanced considerably by the padding and clustering.

Turbulence closure for mixing length theories

NASA Astrophysics Data System (ADS)

Jermyn, Adam S.; Lesaffre, Pierre; Tout, Christopher A.; Chitre, Shashikumar M.

2018-05-01

We present an approach to turbulence closure based on mixing length theory with three-dimensional fluctuations against a two-dimensional background. This model is intended to be rapidly computable for implementation in stellar evolution software and to capture a wide range of relevant phenomena with just a single free parameter, namely the mixing length. We incorporate magnetic, rotational, baroclinic, and buoyancy effects exactly within the formalism of linear growth theories with non-linear decay. We treat differential rotation effects perturbatively in the corotating frame using a novel controlled approximation, which matches the time evolution of the reference frame to arbitrary order. We then implement this model in an efficient open source code and discuss the resulting turbulent stresses and transport coefficients. We demonstrate that this model exhibits convective, baroclinic, and shear instabilities as well as the magnetorotational instability. It also exhibits non-linear saturation behaviour, and we use this to extract the asymptotic scaling of various transport coefficients in physically interesting limits.

COMMIX-PPC: A three-dimensional transient multicomponent computer program for analyzing performance of power plant condensers. Volume 1, Equations and numerics

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

Chien, T.H.; Domanus, H.M.; Sha, W.T.

1993-02-01

The COMMIX-PPC computer pregrain is an extended and improved version of earlier COMMIX codes and is specifically designed for evaluating the thermal performance of power plant condensers. The COMMIX codes are general-purpose computer programs for the analysis of fluid flow and heat transfer in complex Industrial systems. In COMMIX-PPC, two major features have been added to previously published COMMIX codes. One feature is the incorporation of one-dimensional equations of conservation of mass, momentum, and energy on the tube stile and the proper accounting for the thermal interaction between shell and tube side through the porous-medium approach. The other added featuremore » is the extension of the three-dimensional conservation equations for shell-side flow to treat the flow of a multicomponent medium. COMMIX-PPC is designed to perform steady-state and transient. Three-dimensional analysis of fluid flow with heat transfer tn a power plant condenser. However, the code is designed in a generalized fashion so that, with some modification, it can be used to analyze processes in any heat exchanger or other single-phase engineering applications. Volume I (Equations and Numerics) of this report describes in detail the basic equations, formulation, solution procedures, and models for a phenomena. Volume II (User`s Guide and Manual) contains the input instruction, flow charts, sample problems, and descriptions of available options and boundary conditions.« less

Influence of vibrational relaxation on perturbations in a shock layer on a plate

NASA Astrophysics Data System (ADS)

Kirilovskiy, S. V.; Maslov, A. A.; Poplavskaya, T. V.; Tsyryul'nikov, I. S.

2015-05-01

The influence of excitation of molecular vibrational degrees of freedom on the mean flow and perturbation development in a hypersonic (M = 6-14) viscous shock layer is studied. The layer originates on a plate placed in a flow of air, carbon dioxide, or their mixture at high stagnation temperatures (2000-3000 K). The mean flow and pressure pulsation on the surface of the plate are measured in an IT-302M pulsed wind tunnel (Khristianovich Institute of Theoretical and Applied Mechanics, Siberian Branch, Russian Academy of Sciences). Numerical simulation is carried out in terms of a model of a thermally perfect gas using the ANSYS Fluent program package based on solving nonstationary two-dimensional Navier-Stokes equations. External flow perturbations are introduced into the computational domain in the form of plane monochromatic acoustic waves using UDF modules built in the computational code. It is shown that the excitation of vibrational degrees of freedom in carbon dioxide molecules considerably influences the position of the head wave and intensifies perturbations in contrast to air in which the fraction of vibrationally excited molecules is low at the same parameters of the oncoming low. The influence of the excitation of vibrational degrees of freedom is studied both for equilibrium gas and for a vibrationally nonequilibrium gas. Nonequilibrium vibrational degrees of freedom are simulated using a two-temperature model of relaxation flows in which the time variation of the vibrational energy is described by the Landau-Teller equation with regard to a finite time of energy exchange between vibrational and translational-rotational degrees of freedom of molecules. It is found that the vibrational nonequilibrium has a damping effect on perturbations.

Pickup ion processes associated with spacecraft thrusters: Implications for solar probe plus

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

Clemens, Adam, E-mail: a.j.clemens@qmul.ac.uk; Burgess, David

2016-03-15

Chemical thrusters are widely used in spacecraft for attitude control and orbital manoeuvres. They create an exhaust plume of neutral gas which produces ions via photoionization and charge exchange. Measurements of local plasma properties will be affected by perturbations caused by the coupling between the newborn ions and the plasma. A model of neutral expansion has been used in conjunction with a fully three-dimensional hybrid code to study the evolution and ionization over time of the neutral cloud produced by the firing of a mono-propellant hydrazine thruster as well as the interactions of the resulting ion cloud with the ambientmore » solar wind. Results are presented which show that the plasma in the region near to the spacecraft will be perturbed for an extended period of time with the formation of an interaction region around the spacecraft, a moderate amplitude density bow wave bounding the interaction region and evidence of an instability at the forefront of the interaction region which causes clumps of ions to be ejected from the main ion cloud quasi-periodically.« less

Images multiplexing by code division technique

NASA Astrophysics Data System (ADS)

Kuo, Chung J.; Rigas, Harriett

Spread Spectrum System (SSS) or Code Division Multiple Access System (CDMAS) has been studied for a long time, but most of the attention was focused on the transmission problems. In this paper, we study the results when the code division technique is applied to the image at the source stage. The idea is to convolve the N different images with the corresponding m-sequence to obtain the encrypted image. The superimposed image (summation of the encrypted images) is then stored or transmitted. The benefit of this is that no one knows what is stored or transmitted unless the m-sequence is known. The recovery of the original image is recovered by correlating the superimposed image with corresponding m-sequence. Two cases are studied in this paper. First, the two-dimensional image is treated as a long one-dimensional vector and the m-sequence is employed to obtain the results. Secondly, the two-dimensional quasi m-array is proposed and used for the code division multiplexing. It is shown that quasi m-array is faster when the image size is 256 x 256. The important features of the proposed technique are not only the image security but also the data compactness. The compression ratio depends on how many images are superimposed.

CAFE: A New Relativistic MHD Code

NASA Astrophysics Data System (ADS)

Lora-Clavijo, F. D.; Cruz-Osorio, A.; Guzmán, F. S.

2015-06-01

We introduce CAFE, a new independent code designed to solve the equations of relativistic ideal magnetohydrodynamics (RMHD) in three dimensions. We present the standard tests for an RMHD code and for the relativistic hydrodynamics regime because we have not reported them before. The tests include the one-dimensional Riemann problems related to blast waves, head-on collisions of streams, and states with transverse velocities, with and without magnetic field, which is aligned or transverse, constant or discontinuous across the initial discontinuity. Among the two-dimensional (2D) and 3D tests without magnetic field, we include the 2D Riemann problem, a one-dimensional shock tube along a diagonal, the high-speed Emery wind tunnel, the Kelvin-Helmholtz (KH) instability, a set of jets, and a 3D spherical blast wave, whereas in the presence of a magnetic field we show the magnetic rotor, the cylindrical explosion, a case of Kelvin-Helmholtz instability, and a 3D magnetic field advection loop. The code uses high-resolution shock-capturing methods, and we present the error analysis for a combination that uses the Harten, Lax, van Leer, and Einfeldt (HLLE) flux formula combined with a linear, piecewise parabolic method and fifth-order weighted essentially nonoscillatory reconstructors. We use the flux-constrained transport and the divergence cleaning methods to control the divergence-free magnetic field constraint.

A More Accurate and Efficient Technique Developed for Using Computational Methods to Obtain Helical Traveling-Wave Tube Interaction Impedance

NASA Technical Reports Server (NTRS)

Kory, Carol L.

1999-01-01

The phenomenal growth of commercial communications has created a great demand for traveling-wave tube (TWT) amplifiers. Although the helix slow-wave circuit remains the mainstay of the TWT industry because of its exceptionally wide bandwidth, until recently it has been impossible to accurately analyze a helical TWT using its exact dimensions because of the complexity of its geometrical structure. For the first time, an accurate three-dimensional helical model was developed that allows accurate prediction of TWT cold-test characteristics including operating frequency, interaction impedance, and attenuation. This computational model, which was developed at the NASA Lewis Research Center, allows TWT designers to obtain a more accurate value of interaction impedance than is possible using experimental methods. Obtaining helical slow-wave circuit interaction impedance is an important part of the design process for a TWT because it is related to the gain and efficiency of the tube. This impedance cannot be measured directly; thus, conventional methods involve perturbing a helical circuit with a cylindrical dielectric rod placed on the central axis of the circuit and obtaining the difference in resonant frequency between the perturbed and unperturbed circuits. A mathematical relationship has been derived between this frequency difference and the interaction impedance (ref. 1). However, because of the complex configuration of the helical circuit, deriving this relationship involves several approximations. In addition, this experimental procedure is time-consuming and expensive, but until recently it was widely accepted as the most accurate means of determining interaction impedance. The advent of an accurate three-dimensional helical circuit model (ref. 2) made it possible for Lewis researchers to fully investigate standard approximations made in deriving the relationship between measured perturbation data and interaction impedance. The most prominent approximations made in the analysis were addressed and fully investigated for their accuracy by using the three-dimensional electromagnetic simulation code MAFIA (Solution of Maxwell's Equations by the Finite Integration Algorithm) (refs. 3 and 4). We found that several approximations introduced significant error (ref. 5).

A study of the application of singular perturbation theory. [development of a real time algorithm for optimal three dimensional aircraft maneuvers

NASA Technical Reports Server (NTRS)

Mehra, R. K.; Washburn, R. B.; Sajan, S.; Carroll, J. V.

1979-01-01

A hierarchical real time algorithm for optimal three dimensional control of aircraft is described. Systematic methods are developed for real time computation of nonlinear feedback controls by means of singular perturbation theory. The results are applied to a six state, three control variable, point mass model of an F-4 aircraft. Nonlinear feedback laws are presented for computing the optimal control of throttle, bank angle, and angle of attack. Real Time capability is assessed on a TI 9900 microcomputer. The breakdown of the singular perturbation approximation near the terminal point is examined Continuation methods are examined to obtain exact optimal trajectories starting from the singular perturbation solutions.

Transfer function analysis of thermospheric perturbations

NASA Technical Reports Server (NTRS)

Mayr, H. G.; Harris, I.; Varosi, F.; Herrero, F. A.; Spencer, N. W.

1986-01-01

Applying perturbation theory, a spectral model in terms of vectors spherical harmonics (Legendre polynomials) is used to describe the short term thermospheric perturbations originating in the auroral regions. The source may be Joule heating, particle precipitation or ExB ion drift-momentum coupling. A multiconstituent atmosphere is considered, allowing for the collisional momentum exchange between species including Ar, O2, N2, O, He and H. The coupled equations of energy, mass and momentum conservation are solved simultaneously for the major species N2 and O. Applying homogeneous boundary conditions, the integration is carred out from the Earth's surface up to 700 km. In the analysis, the spherical harmonics are treated as eigenfunctions, assuming that the Earth's rotation (and prevailing circulation) do not significantly affect perturbations with periods which are typically much less than one day. Under these simplifying assumptions, and given a particular source distribution in the vertical, a two dimensional transfer function is constructed to describe the three dimensional response of the atmosphere. In the order of increasing horizontal wave numbers (order of polynomials), this transfer function reveals five components. To compile the transfer function, the numerical computations are very time consuming (about 100 hours on a VAX for one particular vertical source distribution). However, given the transfer function, the atmospheric response in space and time (using Fourier integral representation) can be constructed with a few seconds of a central processing unit. This model is applied in a case study of wind and temperature measurements on the Dynamics Explorer B, which show features characteristic of a ringlike excitation source in the auroral oval. The data can be interpreted as gravity waves which are focused (and amplified) in the polar region and then are reflected to propagate toward lower latitudes.

Transfer function analysis of thermospheric perturbations

NASA Astrophysics Data System (ADS)

Mayr, H. G.; Harris, I.; Varosi, F.; Herrero, F. A.; Spencer, N. W.

1986-06-01

Applying perturbation theory, a spectral model in terms of vectors spherical harmonics (Legendre polynomials) is used to describe the short term thermospheric perturbations originating in the auroral regions. The source may be Joule heating, particle precipitation or ExB ion drift-momentum coupling. A multiconstituent atmosphere is considered, allowing for the collisional momentum exchange between species including Ar, O2, N2, O, He and H. The coupled equations of energy, mass and momentum conservation are solved simultaneously for the major species N2 and O. Applying homogeneous boundary conditions, the integration is carred out from the Earth's surface up to 700 km. In the analysis, the spherical harmonics are treated as eigenfunctions, assuming that the Earth's rotation (and prevailing circulation) do not significantly affect perturbations with periods which are typically much less than one day. Under these simplifying assumptions, and given a particular source distribution in the vertical, a two dimensional transfer function is constructed to describe the three dimensional response of the atmosphere. In the order of increasing horizontal wave numbers (order of polynomials), this transfer function reveals five components. To compile the transfer function, the numerical computations are very time consuming (about 100 hours on a VAX for one particular vertical source distribution). However, given the transfer function, the atmospheric response in space and time (using Fourier integral representation) can be constructed with a few seconds of a central processing unit. This model is applied in a case study of wind and temperature measurements on the Dynamics Explorer B, which show features characteristic of a ringlike excitation source in the auroral oval. The data can be interpreted as gravity waves which are focused (and amplified) in the polar region and then are reflected to propagate toward lower latitudes.

Exploration of magnetic perturbation effects on advanced divertor configurations in NSTX-U

DOE PAGES

Frerichs, H.; Schmitz, O.; Waters, I.; ...

2016-06-01

The control of divertor heat loads - both steady state and transient - remains a key challenge for the successful operation of ITER and FNSF. Magnetic perturbations provide a promising technique to control ELMs (transients), but understanding their detailed impact is difficult due to their symmetry breaking nature. One approach for reducing steady state heat loads are so called 'advanced divertors' which aim at optimizing the magnetic field configuration: the snowflake and the (super-)X-divertor. It is likely that both concepts - magnetic perturbations and advanced divertors - will have to work together, and we explore their inter- action based onmore » the NSTX-U setup. An overview of different divertor con gurations under the impact of magnetic perturbations is presented, and the resulting impact on plasma edge transport is investigated with the EMC3-EIRENE code. Variations in size of the magnetic footprint of the perturbed separatrix are found, which is related to the level of flux expansion on the divertor target. Non-axisymmetric peaking of the heat flux related to the perturbed separatrix is found at the outer strike point, but only in locations where flux expansion is not too large.« less

The Effects of Aperiodic Perturbation on the Last Good Surface of a Single-Null Divertor Tokamak

NASA Astrophysics Data System (ADS)

Jordan, Joseph; Punjabi, Alkesh; Ali, Halima

2003-10-01

An area-preserving discrete mapping, called the Simple Map is used to analyze the effects of aperiodic perturbation on the last good magnetic surface of a single-null divertor tokamak. The SM was developed by Punjabi and Boozer /1/. The SM equations are modified to reflect the effect of aperiodic perturbation: Xn + 1 = Xn - (k+δ_aperiodics(n))Y_n(1-Y_n), Y_n+1 = Yn + kX_n+1 where δ_aperiodic is the amplitude of the aperiodic perturbation. The equations are coded into FORTRAN and used in conjunction with a graphing program to show the trajectories for different amplitudes of aperiodic perturbation. As the aperiodic perturbation is increased, the trajectories break from a clean curve and develop strange behavior as they approach chaos. The progression of the trajectories from a clean curve to a chaotic one is observed, and the results are presented in this paper. This project is supported by NASA SHARP and US DOE Grant Number DE-F02-02ER54673. The research was done under the mentorship of Profs. Punjabi and Ali. 1. Punjabi A., Verma A., and Boozer A., Phys. Rev. Lett, 69, 3322 (1992)

Many-body perturbation theory and non-perturbative approaches: screened interaction as the key ingredient

NASA Astrophysics Data System (ADS)

Tarantino, Walter; Mendoza, Bernardo S.; Romaniello, Pina; Berger, J. A.; Reining, Lucia

2018-04-01

Many-body perturbation theory is often formulated in terms of an expansion in the dressed instead of the bare Green’s function, and in the screened instead of the bare Coulomb interaction. However, screening can be calculated on different levels of approximation, and it is important to define what is the most appropriate choice. We explore this question by studying a zero-dimensional model (so called ‘one-point model’) that retains the structure of the full equations. We study both linear and non-linear response approximations to the screening. We find that an expansion in terms of the screening in the random phase approximation is the most promising way for an application in real systems. Moreover, by making use of the nonperturbative features of the Kadanoff-Baym equation for the one-body Green’s function, we obtain an approximate solution in our model that is very promising, although its applicability to real systems has still to be explored.

Shell-crossing in quasi-one-dimensional flow

NASA Astrophysics Data System (ADS)

Rampf, Cornelius; Frisch, Uriel

2017-10-01

Blow-up of solutions for the cosmological fluid equations, often dubbed shell-crossing or orbit crossing, denotes the breakdown of the single-stream regime of the cold-dark-matter fluid. At this instant, the velocity becomes multi-valued and the density singular. Shell-crossing is well understood in one dimension (1D), but not in higher dimensions. This paper is about quasi-one-dimensional (Q1D) flow that depends on all three coordinates but differs only slightly from a strictly 1D flow, thereby allowing a perturbative treatment of shell-crossing using the Euler-Poisson equations written in Lagrangian coordinates. The signature of shell-crossing is then just the vanishing of the Jacobian of the Lagrangian map, a regular perturbation problem. In essence, the problem of the first shell-crossing, which is highly singular in Eulerian coordinates, has been desingularized by switching to Lagrangian coordinates, and can then be handled by perturbation theory. Here, all-order recursion relations are obtained for the time-Taylor coefficients of the displacement field, and it is shown that the Taylor series has an infinite radius of convergence. This allows the determination of the time and location of the first shell-crossing, which is generically shown to be taking place earlier than for the unperturbed 1D flow. The time variable used for these statements is not the cosmic time t but the linear growth time τ ˜ t2/3. For simplicity, calculations are restricted to an Einstein-de Sitter universe in the Newtonian approximation, and tailored initial data are used. However it is straightforward to relax these limitations, if needed.

SCORE-EVET: a computer code for the multidimensional transient thermal-hydraulic analysis of nuclear fuel rod arrays. [BWR; PWR

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

Benedetti, R. L.; Lords, L. V.; Kiser, D. M.

1978-02-01

The SCORE-EVET code was developed to study multidimensional transient fluid flow in nuclear reactor fuel rod arrays. The conservation equations used were derived by volume averaging the transient compressible three-dimensional local continuum equations in Cartesian coordinates. No assumptions associated with subchannel flow have been incorporated into the derivation of the conservation equations. In addition to the three-dimensional fluid flow equations, the SCORE-EVET code ocntains: (a) a one-dimensional steady state solution scheme to initialize the flow field, (b) steady state and transient fuel rod conduction models, and (c) comprehensive correlation packages to describe fluid-to-fuel rod interfacial energy and momentum exchange. Velocitymore » and pressure boundary conditions can be specified as a function of time and space to model reactor transient conditions such as a hypothesized loss-of-coolant accident (LOCA) or flow blockage.« less

The effect of dissipative inhomogeneous medium on the statistics of the wave intensity

NASA Technical Reports Server (NTRS)

Saatchi, Sasan S.

1993-01-01

One of the main theoretical points in the theory of wave propagation in random medium is the derivation of closed form equations to describe the statistics of the propagating waves. In particular, in one dimensional problems, the closed form representation of the multiple scattering effects is important since it contributes in understanding such problems like wave localization, backscattering enhancement, and intensity fluctuations. In this the propagation of plane waves in a layer of one-dimensional dissipative random medium is considered. The medium is modeled by a complex permittivity whose real part is a constant representing the absorption. The one dimensional problem is mathematically equivalent to the analysis of a transmission line with randomly perturbed distributed parameters and a single mode lossy waveguide and the results can be used to study the propagation of radio waves through atmosphere and the remote sensing of geophysical media. It is assumed the scattering medium consists of an ensemble of one-dimensional point scatterers randomly positioned in a layer of thickness L with diffuse boundaries. A Poisson impulse process with density lambda is used to model the position of scatterers in the medium. By employing the Markov properties of this process an exact closed form equation of Kolmogorov-Feller type was obtained for the probability density of the reflection coefficient. This equation was solved by combining two limiting cases: (1) when the density of scatterers is small; and (2) when the medium is weakly dissipative. A two variable perturbation method for small lambda was used to obtain solutions valid for thick layers. These solutions are then asymptotically evaluated for small dissipation. To show the effect of dissipation, the mean and fluctuations of the reflected power are obtained. The results were compared with a lossy homogeneous medium and with a lossless inhomogeneous medium and the regions where the effect of absorption is not essential were discussed.

Magnetic Field Line Random Walk in Arbitrarily Stretched Isotropic Turbulence

NASA Astrophysics Data System (ADS)

Wongpan, P.; Ruffolo, D.; Matthaeus, W. H.; Rowlands, G.

2006-12-01

Many types of space and laboratory plasmas involve turbulent fluctuations with an approximately uniform mean magnetic field B_0, and the field line random walk plays an important role in guiding particle motions. Much of the relevant literature concerns isotropic turbulence, and has mostly been perturbative, i.e., for small fluctuations, or based on numerical simulations for specific conditions. On the other hand, solar wind turbulence is apparently anisotropic, and has been modeled as a sum of idealized two-dimensional and one dimensional (slab) components, but with the deficiency of containing no oblique wave vectors. In the present work, we address the above issues with non-perturbative analytic calculations of diffusive field line random walks for unpolarized, arbitrarily stretched isotropic turbulence, including the limits of nearly one-dimensional (highly stretched) and nearly two-dimensional (highly squashed) turbulence. We develop implicit analytic formulae for the diffusion coefficients D_x and D_z, two coupled integral equations in which D_x and D_z appear inside 3-dimensional integrals over all k-space, are solved numerically with the aid of Mathematica routines for specific cases. We can vary the parameters B0 and β, the stretching along z for constant turbulent energy. Furthermore, we obtain analytic closed-form solutions in all extreme cases. We obtain 0.54 < D_z/D_x < 2, indicating an approximately isotropic random walk even for very anisotropic (unpolarized) turbulence, a surprising result. For a given β, the diffusion coefficient vs. B0 can be described by a Padé approximant. We find quasilinear behavior at high B0 and percolative behavior at low B_0. Partially supported by a Sritrangthong Scholarship from the Faculty of Science, Mahidol University; the Thailand Research Fund; NASA Grant NNG05GG83G; and Thailand's Commission for Higher Education.

Ground state energies from converging and diverging power series expansions

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

Lisowski, C.; Norris, S.; Pelphrey, R.

2016-10-15

It is often assumed that bound states of quantum mechanical systems are intrinsically non-perturbative in nature and therefore any power series expansion methods should be inapplicable to predict the energies for attractive potentials. However, if the spatial domain of the Schrödinger Hamiltonian for attractive one-dimensional potentials is confined to a finite length L, the usual Rayleigh–Schrödinger perturbation theory can converge rapidly and is perfectly accurate in the weak-binding region where the ground state’s spatial extension is comparable to L. Once the binding strength is so strong that the ground state’s extension is less than L, the power expansion becomes divergent,more » consistent with the expectation that bound states are non-perturbative. However, we propose a new truncated Borel-like summation technique that can recover the bound state energy from the diverging sum. We also show that perturbation theory becomes divergent in the vicinity of an avoided-level crossing. Here the same numerical summation technique can be applied to reproduce the energies from the diverging perturbative sums.« less

Type IIB Colliding Plane Waves

NASA Astrophysics Data System (ADS)

Gutperle, M.; Pioline, B.

2003-09-01

Four-dimensional colliding plane wave (CPW) solutions have played an important role in understanding the classical non-linearities of Einstein's equations. In this note, we investigate CPW solutions in 2n+2-dimensional Einstein gravity with a n+1-form flux. By using an isomorphism with the four-dimensional problem, we construct exact solutions analogous to the Szekeres vacuum solution in four dimensions. The higher-dimensional versions of the Khan-Penrose and Bell-Szekeres CPW solutions are studied perturbatively in the vicinity of the light-cone. We find that under small perturbations, a curvature singularity is generically produced, leading to both space-like and time-like singularities. For n = 4, our results pertain to the collision of two ten-dimensional type-IIB Blau-Figueroa o'Farrill-Hull-Papadopoulos plane waves.

DEXTER: A one-dimensional code for calculating thermionic performance of long converters

NASA Technical Reports Server (NTRS)

Sawyer, C. D.

1971-01-01

A versatile code is described for computing the coupled thermionic electric-thermal performance of long thermionic converters in which the temperature and voltage variations cannot be neglected. The code is capable of accounting for a variety of external electrical connection schemes, coolant flow paths and converter failures by partial shorting. Example problem solutions are included along with a user's manual.

Comparison of analytical and experimental performance of a wind-tunnel diffuser section

NASA Technical Reports Server (NTRS)

Shyne, R. J.; Moore, R. D.; Boldman, D. R.

1986-01-01

Wind tunnel diffuser performance is evaluated by comparing experimental data with analytical results predicted by an one-dimensional integration procedure with skin friction coefficient, a two-dimensional interactive boundary layer procedure for analyzing conical diffusers, and a two-dimensional, integral, compressible laminar and turbulent boundary layer code. Pressure, temperature, and velocity data for a 3.25 deg equivalent cone half-angle diffuser (37.3 in., 94.742 cm outlet diameter) was obtained from the one-tenth scale Altitude Wind Tunnel modeling program at the NASA Lewis Research Center. The comparison is performed at Mach numbers of 0.162 (Re = 3.097x19(6)), 0.326 (Re = 6.2737x19(6)), and 0.363 (Re = 7.0129x10(6)). The Reynolds numbers are all based on an inlet diffuser diameter of 32.4 in., 82.296 cm, and reasonable quantitative agreement was obtained between the experimental data and computational codes.

From nonfinite to finite 1D arrays of origami tiles.

PubMed

Wu, Tsai Chin; Rahman, Masudur; Norton, Michael L

2014-06-17

CONSPECTUS: DNA based nanotechnology provides a basis for high-resolution fabrication of objects almost without physical size limitations. However, the pathway to large-scale production of large objects is currently unclear. Operationally, one method forward is to use high information content, large building blocks, which can be generated with high yield and reproducibility. Although flat DNA origami naturally invites comparison to pixels in zero, one, and two dimensions and voxels in three dimensions and has provided an excellent mechanism for generating blocks of significant size and complexity and a multitude of shapes, the field is young enough that a single "brick" has not become the standard platform used by the majority of researchers in the field. In this Account, we highlight factors we considered that led to our adoption of a cross-shaped, non-space-filling origami species, designed by Dr. Liu of the Seeman laboratory, as the building block ideal for use in the fabrication of finite one-dimensional arrays. Three approaches that can be employed for uniquely coding origami-origami linkages are presented. Such coding not only provides the energetics for tethering the species but also uniquely designates the relative orientation of the origami building blocks. The strength of the coding approach implemented in our laboratory is demonstrated using examples of oligomers ranging from finite multimers composed of four, six, and eight origami structures to semi-infinite polymers (100mers). Two approaches to finite array design and the series of assembly steps that each requires are discussed. The process of AFM observation for array characterization is presented as a critical case study. For these soft species, the array images do not simply present the solution phase geometry projected onto a two-dimensional surface. There are additional perturbations associated with fluidic forces associated with sample preparation. At this time, reconstruction of the "true" or average solution structures for blocks is more readily achieved using computer models than using direct imaging methods. The development of scalable 1D-origami arrays composed of uniquely addressable components is a logical, if not necessary, step in the evolution of higher order fully addressable structures. Our research into the fabrication of arrays has led us to generate a listing of several important areas of future endeavor. Of high importance is the re-enforcement of the mechanical properties of the building blocks and the organization of multiple arrays on a surface of technological importance. While addressing this short list of barriers to progress will prove challenging, coherent development along each of these lines of inquiry will accelerate the appearance of commercial scale molecular manufacturing.

Capillary instability of elliptic liquid jets

NASA Astrophysics Data System (ADS)

Amini, Ghobad; Dolatabadi, Ali

2011-08-01

Instability of a liquid jet issuing from an elliptic nozzle in Rayleigh mode is investigated and its behavior is compared with a circular jet. Mathematical solution of viscous free-surface flow for asymmetric geometry is complicated if 3-D analytical solutions are to be obtained. Hence, one-dimensional Cosserat (directed curve) equations are used which can be assumed as a low order form of Navier-Stokes equations for slender jets. Linear solution is performed using perturbation method. Temporal dispersion equation is derived to find the most unstable wavelength responsible for the jet breakup. The obtained results for a circular jet (i.e., an ellipse with an aspect ratio of one) are compared with the classical results of Rayleigh and Weber for inviscid and viscous cases, respectively. It is shown that in the Rayleigh regime, which is the subject of this research, symmetric perturbations are unstable while asymmetric perturbations are stable. Consequently, spatial analysis is performed and the variation of growth rate under the effect of perturbation frequencies for various jet velocities is demonstrated. Results reveal that in comparison with a circular jet, the elliptic jet is more unstable. Furthermore, among liquid jets with elliptical cross sections, those with larger ellipticities have a larger instability growth rate.

Extension of a System Level Tool for Component Level Analysis

NASA Technical Reports Server (NTRS)

Majumdar, Alok; Schallhorn, Paul

2002-01-01

This paper presents an extension of a numerical algorithm for network flow analysis code to perform multi-dimensional flow calculation. The one dimensional momentum equation in network flow analysis code has been extended to include momentum transport due to shear stress and transverse component of velocity. Both laminar and turbulent flows are considered. Turbulence is represented by Prandtl's mixing length hypothesis. Three classical examples (Poiseuille flow, Couette flow and shear driven flow in a rectangular cavity) are presented as benchmark for the verification of the numerical scheme.

Extension of a System Level Tool for Component Level Analysis

NASA Technical Reports Server (NTRS)

Majumdar, Alok; Schallhorn, Paul; McConnaughey, Paul K. (Technical Monitor)

2001-01-01

This paper presents an extension of a numerical algorithm for network flow analysis code to perform multi-dimensional flow calculation. The one dimensional momentum equation in network flow analysis code has been extended to include momentum transport due to shear stress and transverse component of velocity. Both laminar and turbulent flows are considered. Turbulence is represented by Prandtl's mixing length hypothesis. Three classical examples (Poiseuille flow, Couette flow, and shear driven flow in a rectangular cavity) are presented as benchmark for the verification of the numerical scheme.

Inverse Heat Conduction Methods in the CHAR Code for Aerothermal Flight Data Reconstruction

NASA Technical Reports Server (NTRS)

Oliver, A Brandon; Amar, Adam J.

2016-01-01

Reconstruction of flight aerothermal environments often requires the solution of an inverse heat transfer problem, which is an ill-posed problem of specifying boundary conditions from discrete measurements in the interior of the domain. This paper will present the algorithms implemented in the CHAR code for use in reconstruction of EFT-1 flight data and future testing activities. Implementation nuances will be discussed, and alternative hybrid-methods that are permitted by the implementation will be described. Results will be presented for a number of one-dimensional and multi-dimensional problems

Skyshine radiation from a pressurized water reactor containment dome

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

Peng, W.H.

1986-06-01

The radiation dose rates resulting from airborne activities inside a postaccident pressurized water reactor containment are calculated by a discrete ordinates/Monte Carlo combined method. The calculated total dose rates and the skyshine component are presented as a function of distance from the containment at three different elevations for various gamma-ray source energies. The one-dimensional (ANISN code) is used to approximate the skyshine dose rates from the hemisphere dome, and the results are compared favorably to more rigorous results calculated by a three-dimensional Monte Carlo code.

Asymptotic methods for internal transonic flows

NASA Technical Reports Server (NTRS)

Adamson, T. C., Jr.; Messiter, A. F.

1989-01-01

For many internal transonic flows of practical interest, some of the relevant nondimensional parameters typically are small enough that a perturbation scheme can be expected to give a useful level of numerical accuracy. A variety of steady and unsteady transonic channel and cascade flows is studied with the help of systematic perturbation methods which take advantage of this fact. Asymptotic representations are constructed for small changes in channel cross-section area, small flow deflection angles, small differences between the flow velocity and the sound speed, small amplitudes of imposed oscillations, and small reduced frequencies. Inside a channel the flow is nearly one-dimensional except in thin regions immediately downstream of a shock wave, at the channel entrance and exit, and near the channel throat. A study of two-dimensional cascade flow is extended to include a description of three-dimensional compressor-rotor flow which leads to analytical results except in thin edge regions which require numerical solution. For unsteady flow the qualitative nature of the shock-wave motion in a channel depends strongly on the orders of magnitude of the frequency and amplitude of impressed wall oscillations or fluctuations in back pressure. One example of supersonic flow is considered, for a channel with length large compared to its width, including the effect of separation bubbles and the possibility of self-sustained oscillations. The effect of viscosity on a weak shock wave in a channel is discussed.

Streamer properties and associated x-rays in perturbed air

NASA Astrophysics Data System (ADS)

Köhn, C.; Chanrion, O.; Babich, L. P.; Neubert, T.

2018-01-01

Streamers are ionization waves in electric discharges. One of the key ingredients of streamer propagation is an ambient gas that serves as a source of free electrons. Here, we explore the dependence of streamer dynamics on different spatial distributions of ambient air molecules. We vary the spatial profile of air parallel and perpendicular to the ambient electric field. We consider local sinusoidal perturbations of 5%-100%, as induced from discharge shock waves. We use a cylindrically symmetric particle-in-cell code to simulate the evolution of bidirectional streamers and compare the electron density, electric field, streamer velocity and electron energy of streamers in uniform air and in perturbed air. In all considered cases, the motion is driven along in decreasing air density and damped along increasing air density. Perturbations of at most 5%-10% change the velocity differences by up to approximately 40%. Perturbations perpendicular to the electric field additionally squeeze or branch streamers. Air variations can thus partly explain the difference of velocities and morphologies of streamer discharges. In cases with large perturbations, electrons gain energies of up to 30 keV compared to 100 eV in uniformly distributed air. For such perturbations parallel to the ambient electric field, we see the spontaneous initiation of a negative streamer; for perpendicular perturbations, x-rays with energies of up to 20 keV are emitted within 0.17 ns.

Accuracy and Resolution in Micro-earthquake Tomographic Inversion Studies

NASA Astrophysics Data System (ADS)

Hutchings, L. J.; Ryan, J.

2010-12-01

Accuracy and resolution are complimentary properties necessary to interpret the results of earthquake location and tomography studies. Accuracy is the how close an answer is to the “real world”, and resolution is who small of node spacing or earthquake error ellipse one can achieve. We have modified SimulPS (Thurber, 1986) in several ways to provide a tool for evaluating accuracy and resolution of potential micro-earthquake networks. First, we provide synthetic travel times from synthetic three-dimensional geologic models and earthquake locations. We use this to calculate errors in earthquake location and velocity inversion results when we perturb these models and try to invert to obtain these models. We create as many stations as desired and can create a synthetic velocity model with any desired node spacing. We apply this study to SimulPS and TomoDD inversion studies. “Real” travel times are perturbed with noise and hypocenters are perturbed to replicate a starting location away from the “true” location, and inversion is performed by each program. We establish travel times with the pseudo-bending ray tracer and use the same ray tracer in the inversion codes. This, of course, limits our ability to test the accuracy of the ray tracer. We developed relationships for the accuracy and resolution expected as a function of the number of earthquakes and recording stations for typical tomographic inversion studies. Velocity grid spacing started at 1km, then was decreased to 500m, 100m, 50m and finally 10m to see if resolution with decent accuracy at that scale was possible. We considered accuracy to be good when we could invert a velocity model perturbed by 50% back to within 5% of the original model, and resolution to be the size of the grid spacing. We found that 100 m resolution could obtained by using 120 stations with 500 events, bu this is our current limit. The limiting factors are the size of computers needed for the large arrays in the inversion and a realistic number of stations and events needed to provide the data.

Statistics of Smoothed Cosmic Fields in Perturbation Theory. I. Formulation and Useful Formulae in Second-Order Perturbation Theory

NASA Astrophysics Data System (ADS)

Matsubara, Takahiko

2003-02-01

We formulate a general method for perturbative evaluations of statistics of smoothed cosmic fields and provide useful formulae for application of the perturbation theory to various statistics. This formalism is an extensive generalization of the method used by Matsubara, who derived a weakly nonlinear formula of the genus statistic in a three-dimensional density field. After describing the general method, we apply the formalism to a series of statistics, including genus statistics, level-crossing statistics, Minkowski functionals, and a density extrema statistic, regardless of the dimensions in which each statistic is defined. The relation between the Minkowski functionals and other geometrical statistics is clarified. These statistics can be applied to several cosmic fields, including three-dimensional density field, three-dimensional velocity field, two-dimensional projected density field, and so forth. The results are detailed for second-order theory of the formalism. The effect of the bias is discussed. The statistics of smoothed cosmic fields as functions of rescaled threshold by volume fraction are discussed in the framework of second-order perturbation theory. In CDM-like models, their functional deviations from linear predictions plotted against the rescaled threshold are generally much smaller than that plotted against the direct threshold. There is still a slight meatball shift against rescaled threshold, which is characterized by asymmetry in depths of troughs in the genus curve. A theory-motivated asymmetry factor in the genus curve is proposed.

Transport and equilibrium in field-reversed mirrors

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

Boyd, J.K.

Two plasma models relevant to compact torus research have been developed to study transport and equilibrium in field reversed mirrors. In the first model for small Larmor radius and large collision frequency, the plasma is described as an adiabatic hydromagnetic fluid. In the second model for large Larmor radius and small collision frequency, a kinetic theory description has been developed. Various aspects of the two models have been studied in five computer codes ADB, AV, NEO, OHK, RES. The ADB code computes two dimensional equilibrium and one dimensional transport in a flux coordinate. The AV code calculates orbit average integralsmore » in a harmonic oscillator potential. The NEO code follows particle trajectories in a Hill's vortex magnetic field to study stochasticity, invariants of the motion, and orbit average formulas. The OHK code displays analytic psi(r), B/sub Z/(r), phi(r), E/sub r/(r) formulas developed for the kinetic theory description. The RES code calculates resonance curves to consider overlap regions relevant to stochastic orbit behavior.« less

Global Magnetohydrodynamic Simulation Using High Performance FORTRAN on Parallel Computers

NASA Astrophysics Data System (ADS)

Ogino, T.

High Performance Fortran (HPF) is one of modern and common techniques to achieve high performance parallel computation. We have translated a 3-dimensional magnetohydrodynamic (MHD) simulation code of the Earth's magnetosphere from VPP Fortran to HPF/JA on the Fujitsu VPP5000/56 vector-parallel supercomputer and the MHD code was fully vectorized and fully parallelized in VPP Fortran. The entire performance and capability of the HPF MHD code could be shown to be almost comparable to that of VPP Fortran. A 3-dimensional global MHD simulation of the earth's magnetosphere was performed at a speed of over 400 Gflops with an efficiency of 76.5 VPP5000/56 in vector and parallel computation that permitted comparison with catalog values. We have concluded that fluid and MHD codes that are fully vectorized and fully parallelized in VPP Fortran can be translated with relative ease to HPF/JA, and a code in HPF/JA may be expected to perform comparably to the same code written in VPP Fortran.

Shock-turbulence interaction in core-collapse supernovae

NASA Astrophysics Data System (ADS)

Abdikamalov, Ernazar; Zhaksylykov, Azamat; Radice, David; Berdibek, Shapagat

2016-10-01

Nuclear shell burning in the final stages of the lives of massive stars is accompanied by strong turbulent convection. The resulting fluctuations aid supernova explosion by amplifying the non-radial flow in the post-shock region. In this work, we investigate the physical mechanism behind this amplification using a linear perturbation theory. We model the shock wave as a one-dimensional planar discontinuity and consider its interaction with vorticity and entropy perturbations in the upstream flow. We find that, as the perturbations cross the shock, their total turbulent kinetic energy is amplified by a factor of ˜2, while the average linear size of turbulent eddies decreases by about the same factor. These values are not sensitive to the parameters of the upstream turbulence and the nuclear dissociation efficiency at the shock. Finally, we discuss the implication of our results for the supernova explosion mechanism. We show that the upstream perturbations can decrease the critical neutrino luminosity for producing explosion by several per cent.

Shock wave viscosity measurements

NASA Astrophysics Data System (ADS)

Celliers, Peter

2013-06-01

Several decades ago a method was proposed and demonstrated to measure the viscosity of fluids at high pressure by observing the oscillatory damping of sinusoidal perturbations on a shock front. A detailed mathematical analysis of the technique carried out subsequently by Miller and Ahrens revealed its potential, as well as a deep level of complexity in the analysis. We revisit the ideas behind this technique in the context of a recent experimental development: two-dimensional imaging velocimetry. The new technique allows one to capture a broad spectrum of perturbations down to few micron scale-lengths imposed on a shock front from an initial perturbation. The detailed evolution of the perturbation spectrum is sensitive to the viscosity in the fluid behind the shock front. Initial experiments are aimed at examining the viscosity of shock compressed SiO2 just above the shock melting transition. This work performed under the auspices of the U.S. Department of Energy by Lawrence Livermore National Laboratory under Contract DE-AC52-07NA27344.

Pechukas-Yukawa approach to the evolution of the quantum state of a parametrically perturbed system

NASA Astrophysics Data System (ADS)

Qureshi, Mumnuna A.; Zhong, Johnny; Qureshi, Zihad; Mason, Peter; Betouras, Joseph J.; Zagoskin, Alexandre M.

2018-03-01

We consider the evolution of the quantum states of a Hamiltonian that is parametrically perturbed via a term proportional to the adiabatic parameter λ (t ) . Starting with the Pechukas-Yukawa mapping of the energy eigenvalue evolution in a generalized Calogero-Sutherland model of a one-dimensional classical gas, we consider the adiabatic approximation with two different expansions of the quantum state in powers of d λ /d t and compare them with a direct numerical simulation. We show that one of these expansions (Magnus series) is especially convenient for the description of nonadiabatic evolution of the system. Applying the expansion to the exact cover 3-satisfiability problem, we obtain the occupation dynamics, which provides insight into the population of states and sources of decoherence in a quantum system.

Motion makes sense: an adaptive motor-sensory strategy underlies the perception of object location in rats.

PubMed

Saraf-Sinik, Inbar; Assa, Eldad; Ahissar, Ehud

2015-06-10

Tactile perception is obtained by coordinated motor-sensory processes. We studied the processes underlying the perception of object location in freely moving rats. We trained rats to identify the relative location of two vertical poles placed in front of them and measured at high resolution the motor and sensory variables (19 and 2 variables, respectively) associated with this whiskers-based perceptual process. We found that the rats developed stereotypic head and whisker movements to solve this task, in a manner that can be described by several distinct behavioral phases. During two of these phases, the rats' whiskers coded object position by first temporal and then angular coding schemes. We then introduced wind (in two opposite directions) and remeasured their perceptual performance and motor-sensory variables. Our rats continued to perceive object location in a consistent manner under wind perturbations while maintaining all behavioral phases and relatively constant sensory coding. Constant sensory coding was achieved by keeping one group of motor variables (the "controlled variables") constant, despite the perturbing wind, at the cost of strongly modulating another group of motor variables (the "modulated variables"). The controlled variables included coding-relevant variables, such as head azimuth and whisker velocity. These results indicate that consistent perception of location in the rat is obtained actively, via a selective control of perception-relevant motor variables. Copyright © 2015 the authors 0270-6474/15/358777-13$15.00/0.

Using the theory of small perturbations in performance calculations of the RBMK

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

Isaev, N.V.; Druzhinin, V.E.; Pogosbekyan, L.R.

The theory of small perturbations in reactor physics is discussed and applied to two-dimensional calculations of the RBMK. The classical theory of small perturbations implies considerable errors in calculations because the perturbations cannot be considered small. The modified theory of small perturbations presented here can be used in atomic power stations for determining reactivity effects and reloading rates of channels in reactors and also for assessing the reactivity storage in control rods.

An efficient HZETRN (a galactic cosmic ray transport code)

NASA Technical Reports Server (NTRS)

Shinn, Judy L.; Wilson, John W.

1992-01-01

An accurate and efficient engineering code for analyzing the shielding requirements against the high-energy galactic heavy ions is needed. The HZETRN is a deterministic code developed at Langley Research Center that is constantly under improvement both in physics and numerical computation and is targeted for such use. One problem area connected with the space-marching technique used in this code is the propagation of the local truncation error. By improving the numerical algorithms for interpolation, integration, and grid distribution formula, the efficiency of the code is increased by a factor of eight as the number of energy grid points is reduced. The numerical accuracy of better than 2 percent for a shield thickness of 150 g/cm(exp 2) is found when a 45 point energy grid is used. The propagating step size, which is related to the perturbation theory, is also reevaluated.

Verification of the helioseismology travel-time measurement technique and the inversion procedure for sound speed using artificial data

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

Parchevsky, K. V.; Zhao, J.; Hartlep, T.

We performed three-dimensional numerical simulations of the solar surface acoustic wave field for the quiet Sun and for three models with different localized sound-speed perturbations in the interior with deep, shallow, and two-layer structures. We used the simulated data generated by two solar acoustics codes that employ the same standard solar model as a background model, but utilize different integration techniques and different models of stochastic wave excitation. Acoustic travel times were measured using a time-distance helioseismology technique, and compared with predictions from ray theory frequently used for helioseismic travel-time inversions. It is found that the measured travel-time shifts agreemore » well with the helioseismic theory for sound-speed perturbations, and for the measurement procedure with and without phase-speed filtering of the oscillation signals. This testing verifies the whole measuring-filtering-inversion procedure for static sound-speed anomalies with small amplitude inside the Sun outside regions of strong magnetic field. It is shown that the phase-speed filtering, frequently used to extract specific wave packets and improve the signal-to-noise ratio, does not introduce significant systematic errors. Results of the sound-speed inversion procedure show good agreement with the perturbation models in all cases. Due to its smoothing nature, the inversion procedure may overestimate sound-speed variations in regions with sharp gradients of the sound-speed profile.« less

Differentially Private Synthesization of Multi-Dimensional Data using Copula Functions

PubMed Central

Li, Haoran; Xiong, Li; Jiang, Xiaoqian

2014-01-01

Differential privacy has recently emerged in private statistical data release as one of the strongest privacy guarantees. Most of the existing techniques that generate differentially private histograms or synthetic data only work well for single dimensional or low-dimensional histograms. They become problematic for high dimensional and large domain data due to increased perturbation error and computation complexity. In this paper, we propose DPCopula, a differentially private data synthesization technique using Copula functions for multi-dimensional data. The core of our method is to compute a differentially private copula function from which we can sample synthetic data. Copula functions are used to describe the dependence between multivariate random vectors and allow us to build the multivariate joint distribution using one-dimensional marginal distributions. We present two methods for estimating the parameters of the copula functions with differential privacy: maximum likelihood estimation and Kendall’s τ estimation. We present formal proofs for the privacy guarantee as well as the convergence property of our methods. Extensive experiments using both real datasets and synthetic datasets demonstrate that DPCopula generates highly accurate synthetic multi-dimensional data with significantly better utility than state-of-the-art techniques. PMID:25405241

Multitasking the INS3D-LU code on the Cray Y-MP

NASA Technical Reports Server (NTRS)

Fatoohi, Rod; Yoon, Seokkwan

1991-01-01

This paper presents the results of multitasking the INS3D-LU code on eight processors. The code is a full Navier-Stokes solver for incompressible fluid in three dimensional generalized coordinates using a lower-upper symmetric-Gauss-Seidel implicit scheme. This code has been fully vectorized on oblique planes of sweep and parallelized using autotasking with some directives and minor modifications. The timing results for five grid sizes are presented and analyzed. The code has achieved a processing rate of over one Gflops.

Non-Abelian supertubes

NASA Astrophysics Data System (ADS)

Fernández-Melgarejo, José J.; Park, Minkyu; Shigemori, Masaki

2017-12-01

A supertube is a supersymmetric configuration in string theory which occurs when a pair of branes spontaneously polarizes and generates a new dipole charge extended along a closed curve. The dipole charge of a codimension-2 supertube is characterized by the U-duality monodromy as one goes around the supertube. For multiple codimension-2 supertubes, their monodromies do not commute in general. In this paper, we construct a supersymmetric solution of five-dimensional supergravity that describes two supertubes with such non-Abelian monodromies, in a certain perturbative expansion. In supergravity, the monodromies are realized as the multi-valuedness of the scalar fields, while in higher dimensions they correspond to non-geometric duality twists of the internal space. The supertubes in our solution carry NS5 and 5 2 2 dipole charges and exhibit the same monodromy structure as the SU(2) Seiberg-Witten geometry. The perturbative solution has AdS2 × S 2 asymptotics and vanishing four-dimensional angular momentum. We argue that this solution represents a microstate of four-dimensional black holes with a finite horizon and that it provides a clue for the gravity realization of a pure-Higgs branch state in the dual quiver quantum mechanics.

Application of Linear and Non-Linear Harmonic Methods for Unsteady Transonic Flow

NASA Astrophysics Data System (ADS)

Gundevia, Rayomand

This thesis explores linear and non-linear computational methods for solving unsteady flow. The eventual goal is to apply these methods to two-dimensional and three-dimensional flutter predictions. In this study the quasi-one-dimensional nozzle is used as a framework for understanding these methods and their limitations. Subsonic and transonic cases are explored as the back-pressure is forced to oscillate with known amplitude and frequency. A steady harmonic approach is used to solve this unsteady problem for which perturbations are said to be small in comparison to the mean flow. The use of a linearized Euler equations (LEE) scheme is good at capturing the flow characteristics but is limited by accuracy to relatively small amplitude perturbations. The introduction of time-averaged second-order terms in the Non-Linear Harmonic (NLH) method means that a better approximation of the mean-valued solution, upon which the linearization is based, can be made. The nonlinear time-accurate Euler solutions are used for comparison and to establish the regimes of unsteadiness for which these schemes fails. The usefulness of the LEE and NLH methods lie in the gains in computational efficiency over the full equations.

Revisiting Yasinsky and Henry`s benchmark using modern nodal codes

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

Feltus, M.A.; Becker, M.W.

1995-12-31

The numerical experiments analyzed by Yasinsky and Henry are quite trivial by comparison with today`s standards because they used the finite difference code WIGLE for their benchmark. Also, this problem is a simple slab (one-dimensional) case with no feedback mechanisms. This research attempts to obtain STAR (Ref. 2) and NEM (Ref. 3) code results in order to produce a more modern kinetics benchmark with results comparable WIGLE.

Anderson localized state as a predissipative state: irreversible emission of thermalized quanta from a dynamically delocalized state.

PubMed

Yamada, Hiroaki; Ikeda, Kensuke S

2002-04-01

It was shown that localization in one-dimensional disordered (quantum) electronic system is destroyed against coherent harmonic perturbations and the delocalized electron exhibits an unlimited diffusive motion [Yamada and Ikeda, Phys. Rev. E 59, 5214 (1999)]. The appearance of diffusion implies that the system has potential for irreversibility and dissipation. In the present paper, we investigate dissipative property of the dynamically delocalized state, and we show that an irreversible quasistationary energy flow indeed appears in the form of a "heat" flow when we couple the system with another dynamical degree of freedom. In the concrete we numerically investigate dissipative properties of a one-dimensional tight-binding electronic system perturbed by time-dependent harmonic forces, by coupling it with a quantum harmonic oscillator or a quantum anharmonic oscillator. It is demonstrated that if the on-site potential is spatially irregular an irreversible energy transfer from the scattered electron to the test oscillator occurs. Moreover, the test oscillator promptly approaches a thermalized state characterized by a well-defined time-dependent temperature. On the contrary, such a relaxation process cannot be observed at all for periodic potential systems. Our system is one of the minimal quantum systems in which a distinct nonequilibrium statistical behavior is self-induced.

Selection rule engineering of forbidden transitions of a hydrogen atom near a nanogap

NASA Astrophysics Data System (ADS)

Kim, Hyunyoung Y.; Kim, Daisik S.

2018-01-01

We perform an analytical study on the allowance of forbidden transitions for a hydrogen atom placed near line dipole sources, mimicking light emanating from a one-dimensional metallic nanogap. It is shown that the rapid variation of the electric field vector, inevitable in the near zone, completely breaks the selection rule of Δl=±1. While the forbidden transitions between spherically symmetric S states, such as 2S to 1S or 3S to 1S (Δl=0), are rather robust against selection rule breakage, Δl=±2 transitions such as between 3D and 1S or 3D and 2S states are very vulnerable to the spatial variation of the perturbing electric field. Transitions between 2S and 3D states are enhanced by many orders of magnitude, aided by the quadratic nature of both the perturbing Hamiltonian and D wavefunctions. The forbidden dipole moment, which approaches one Bohr radius times the electric charge in the vicinity of the gap, can be written in a simple closed form owing to the one-dimensional nature of our gap. With large enough effective volume together with the symmetric nature of the excited state wavefunctions, our work paves way towards atomic physics application of infinitely long nanogaps.

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

Morgan, B. E.; Olson, B. J.; White, J. E.

High-fidelity large eddy simulation (LES) of a low-Atwood number (A = 0.05) Rayleigh-Taylor mixing layer is performed using the tenth-order compact difference code Miranda. An initial multimode perturbation spectrum is specified in Fourier space as a function of mesh resolution such that a database of results is obtained in which each successive level of increased grid resolution corresponds approximately to one additional doubling of the mixing layer width, or generation. The database is then analyzed to determine approximate requirements for self-similarity, and a new metric is proposed to quantify how far a given simulation is from the limit of self-similarity.more » It is determined that mixing layer growth reaches a high degree of self-similarity after approximately 4.5 generations. Statistical convergence errors and boundary effects at late time, however, make it impossible to draw similar conclusions regarding the self-similar growth of more sensitive turbulence parameters. Finally, self-similar turbulence profiles from the LES database are compared with one-dimensional simulations using the k-L-a and BHR-2 Reynolds-averaged Navier-Stokes (RANS) models. The k-L-a model, which is calibrated to reproduce a quadratic turbulence kinetic energy profile for a self-similar mixing layer, is found to be in better agreement with the LES than BHR-2 results.« less

Vapor-liquid equilibrium and equation of state of two-dimensional fluids from a discrete perturbation theory

NASA Astrophysics Data System (ADS)

Trejos, Víctor M.; Santos, Andrés; Gámez, Francisco

2018-05-01

The interest in the description of the properties of fluids of restricted dimensionality is growing for theoretical and practical reasons. In this work, we have firstly developed an analytical expression for the Helmholtz free energy of the two-dimensional square-well fluid in the Barker-Henderson framework. This equation of state is based on an approximate analytical radial distribution function for d-dimensional hard-sphere fluids (1 ≤ d ≤ 3) and is validated against existing and new simulation results. The so-obtained equation of state is implemented in a discrete perturbation theory able to account for general potential shapes. The prototypical Lennard-Jones and Yukawa fluids are tested in its two-dimensional version against available and new simulation data with semiquantitative agreement.

Density matrix perturbation theory for magneto-optical response of periodic insulators

NASA Astrophysics Data System (ADS)

Lebedeva, Irina; Tokatly, Ilya; Rubio, Angel

2015-03-01

Density matrix perturbation theory offers an ideal theoretical framework for the description of response of solids to arbitrary electromagnetic fields. In particular, it allows to consider perturbations introduced by uniform electric and magnetic fields under periodic boundary conditions, though the corresponding potentials break the translational invariance of the Hamiltonian. We have implemented the density matrix perturbation theory in the open-source Octopus code on the basis of the efficient Sternheimer approach. The procedures for responses of different order to electromagnetic fields, including electric polarizability, orbital magnetic susceptibility and magneto-optical response, have been developed and tested by comparison with the results for finite systems and for wavefunction-based perturbation theory, which is already available in the code. Additional analysis of the orbital magneto-optical response is performed on the basis of analytical models. Symmetry limitations to observation of the magneto-optical response are discussed. The financial support from the Marie Curie Fellowship PIIF-GA-2012-326435 (RespSpatDisp) is gratefully acknowledged.

A new line-of-sight approach to the non-linear Cosmic Microwave Background

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

Fidler, Christian; Koyama, Kazuya; Pettinari, Guido W., E-mail: christian.fidler@port.ac.uk, E-mail: kazuya.koyama@port.ac.uk, E-mail: guido.pettinari@gmail.com

2015-04-01

We develop the transport operator formalism, a new line-of-sight integration framework to calculate the anisotropies of the Cosmic Microwave Background (CMB) at the linear and non-linear level. This formalism utilises a transformation operator that removes all inhomogeneous propagation effects acting on the photon distribution function, thus achieving a split between perturbative collisional effects at recombination and non-perturbative line-of-sight effects at later times. The former can be computed in the framework of standard cosmological perturbation theory with a second-order Boltzmann code such as SONG, while the latter can be treated within a separate perturbative scheme allowing the use of non-linear Newtonianmore » potentials. We thus provide a consistent framework to compute all physical effects contained in the Boltzmann equation and to combine the standard remapping approach with Boltzmann codes at any order in perturbation theory, without assuming that all sources are localised at recombination.« less

Three ways to solve critical ϕ4 theory on 4 ‑ 𝜖 dimensional real projective space: Perturbation, bootstrap, and Schwinger-Dyson equation

NASA Astrophysics Data System (ADS)

Hasegawa, Chika; Nakayama, Yu

2018-03-01

In this paper, we solve the two-point function of the lowest dimensional scalar operator in the critical ϕ4 theory on 4 ‑ 𝜖 dimensional real projective space in three different methods. The first is to use the conventional perturbation theory, and the second is to impose the cross-cap bootstrap equation, and the third is to solve the Schwinger-Dyson equation under the assumption of conformal invariance. We find that the three methods lead to mutually consistent results but each has its own advantage.

Three-dimensional modeling of electron quasiviscous dissipation in guide-field magnetic reconnection

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

Hesse, Michael; Kuznetsova, Masha; Schindler, Karl

2005-10-01

A numerical study of guide-field magnetic reconnection in a three-dimensional model is presented. Starting from an initial, perturbed, force-free current sheet, it is shown that reconnection develops to an almost translationally invariant state, where magnetic perturbations are aligned primarily along the main current flow direction. An analysis of guide-field and electron flow signatures indicates behavior that is very similar to earlier, albeit not three-dimensional, simulations. Furthermore, a detailed investigation of electron pressure nongyrotropies in the central diffusion region confirms the major role the associated dissipation process plays in establishing the reconnection electric field.

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

Lyons, Brendan C.; Ferraro, Nathaniel M.; Paz-Soldan, Carlos A.

In order to understand the effect of rotation on the plasma's response to three-dimensional magnetic perturbations, we perform a systematic scan of the zero-crossing of the rotation profile in a DIII-D ITER-similar shape equilibrium using linear, time-independent modeling with the M3D-C1 extended magnetohydrodynamics code. We confirm that the local resonant magnetic field generally increases as the rotation decreases at a rational surface. Multiple peaks in the resonant field are observed near rational surfaces, however, and the maximum resonant field does not always correspond to zero rotation at the surface. Furthermore, we show that non-resonant current can be driven at zero-more » crossings not aligned with rational surfaces if there is sufficient shear in the rotation profile there, leading to an amplification of near-resonant Fourier harmonics of the perturbed magnetic field and a decrease in the far-off -resonant harmonics. The quasilinear electromagnetic torque induced by this non-resonant plasma response provides drive to flatten the rotation, possibly allowing for increased transport in the pedestal by the destabilization of turbulent modes. In addition, this torque acts to drive the rotation zero-crossing to dynamically stable points near rational surfaces, which would allow for increased resonant penetration. By one or both of these mechanisms, this torque may play an important role in bifurcations into ELM suppression. Finally, we discuss how these changes to the plasma response could be detected by tokamak diagnostics. In particular, we show that the changes to the resonant field discussed here have a significant impact on the external perturbed magnetic field, which should be observable by magnetic sensors on the high-field side of tokamaks, but not on the low-field side. In addition, TRIP3D-MAFOT simulations show that none of the changes to the plasma response described here substantially affects the divertor footprint structure.« less

Self-consistent perturbation theory for two dimensional twisted bilayers

NASA Astrophysics Data System (ADS)

Shirodkar, Sharmila N.; Tritsaris, Georgios A.; Kaxiras, Efthimios

Theoretical modeling and ab-initio simulations of two dimensional heterostructures with arbitrary angles of rotation between layers involve unrealistically large and expensive calculations. To overcome this shortcoming, we develop a methodology for weakly interacting heterostructures that treats the effect of one layer on the other as perturbation, and restricts the calculations to their primitive cells. Thus, avoiding computationally expensive supercells. We start by approximating the interaction potential between the twisted bilayers to that of a hypothetical configuration (viz. ideally stacked untwisted layers), which produces band structures in reasonable agreement with full-scale ab-initio calculations for commensurate and twisted bilayers of graphene (Gr) and Gr/hexagonal boron nitride (h-BN) heterostructures. We then self-consistently calculate the charge density and hence, interaction potential of the heterostructures. In this work, we test our model for bilayers of various combinations of Gr, h-BN and transition metal dichalcogenides, and discuss the advantages and shortcomings of the self-consistently calculated interaction potential. Department of Physics, Harvard University, Cambridge, Massachusetts 02138, USA.

Two-dimensional modulated ion-acoustic excitations in electronegative plasmas

NASA Astrophysics Data System (ADS)

Panguetna, Chérif S.; Tabi, Conrad B.; Kofané, Timoléon C.

2017-09-01

Two-dimensional modulated ion-acoustic waves are investigated in an electronegative plasma. Through the reductive perturbation expansion, the governing hydrodynamic equations are reduced to a Davey-Stewartson system with two-space variables. The latter is used to study the modulational instability of ion-acoustic waves along with the effect of plasma parameters, namely, the negative ion concentration ratio (α) and the electron-to-negative ion temperature ratio (σn). A parametric analysis of modulational instability is carried out, where regions of plasma parameters responsible for the emergence of modulated ion-acoustic waves are discussed, with emphasis on the behavior of the instability growth rate. Numerically, using perturbed plane waves as initial conditions, parameters from the instability regions give rise to series of dromion solitons under the activation of modulational instability. The sensitivity of the numerical solutions to plasma parameters is discussed. Some exact solutions in the form one- and two-dromion solutions are derived and their response to the effect of varying α and σn is discussed as well.

Moduli stabilising in heterotic nearly Kähler compactifications

NASA Astrophysics Data System (ADS)

Klaput, Michael; Lukas, Andre; Matti, Cyril; Svanes, Eirik E.

2013-01-01

We study heterotic string compactifications on nearly Kähler homogeneous spaces, including the gauge field effects which arise at order α'. Using Abelian gauge fields, we are able to solve the Bianchi identity and supersymmetry conditions to this order. The four-dimensional external space-time consists of a domain wall solution with moduli fields varying along the transverse direction. We find that the inclusion of α' corrections improves the moduli stabilization features of this solution. In this case, one of the dilaton and the volume modulus asymptotes to a constant value away from the domain wall. It is further shown that the inclusion of non-perturbative effects can stabilize the remaining modulus and "lift" the domain wall to an AdS vacuum. The coset SU(3)/U(1)2 is used as an explicit example to demonstrate the validity of this AdS vacuum. Our results show that heterotic nearly Kähler compactifications can lead to maximally symmetric four-dimensional space-times at the non-perturbative level.

The development of an explicit thermochemical nonequilibrium algorithm and its application to compute three dimensional AFE flowfields

NASA Technical Reports Server (NTRS)

Palmer, Grant

1989-01-01

This study presents a three-dimensional explicit, finite-difference, shock-capturing numerical algorithm applied to viscous hypersonic flows in thermochemical nonequilibrium. The algorithm employs a two-temperature physical model. Equations governing the finite-rate chemical reactions are fully-coupled to the gas dynamic equations using a novel coupling technique. The new coupling method maintains stability in the explicit, finite-rate formulation while allowing relatively large global time steps. The code uses flux-vector accuracy. Comparisons with experimental data and other numerical computations verify the accuracy of the present method. The code is used to compute the three-dimensional flowfield over the Aeroassist Flight Experiment (AFE) vehicle at one of its trajectory points.

Perturbation Effects on a Supercritical C7H16/N2 Mixing Layer

NASA Technical Reports Server (NTRS)

Okongo'o, Nora; Bellan, Josette

2008-01-01

A computational-simulation study has been presented of effects of perturbation wavelengths and initial Reynolds numbers on the transition to turbulence of a heptane/nitrogen mixing layer at supercritical pressure. The governing equations for the simulations were the same as those of related prior studies reported in NASA Tech Briefs. Two-dimensional (2D) simulations were performed with initially im posed span wise perturbations whereas three-dimensional (3D) simulations had both streamwise and spanwise initial perturbations. The 2D simulations were undertaken to ascertain whether perturbations having the shortest unstable wavelength obtained from a linear stability analysis for inviscid flow are unstable in viscous nonlinear flows. The goal of the 3D simulations was to ascertain whether perturbing the mixing layer at different wavelengths affects the transition to turbulence. It was found that transitions to turbulence can be obtained at different perturbation wavelengths, provided that they are longer than the shortest unstable wavelength as determined by 2D linear stability analysis for the inviscid case and that the initial Reynolds number is proportionally increased as the wavelength is decreased. The transitional states thus obtained display different dynamic and mixture characteristics, departing strongly from the behaviors of perfect gases and ideal mixtures.

Bar coded retroreflective target

DOEpatents

Vann, Charles S.

2000-01-01

This small, inexpensive, non-contact laser sensor can detect the location of a retroreflective target in a relatively large volume and up to six degrees of position. The tracker's laser beam is formed into a plane of light which is swept across the space of interest. When the beam illuminates the retroreflector, some of the light returns to the tracker. The intensity, angle, and time of the return beam is measured to calculate the three dimensional location of the target. With three retroreflectors on the target, the locations of three points on the target are measured, enabling the calculation of all six degrees of target position. Until now, devices for three-dimensional tracking of objects in a large volume have been heavy, large, and very expensive. Because of the simplicity and unique characteristics of this tracker, it is capable of three-dimensional tracking of one to several objects in a large volume, yet it is compact, light-weight, and relatively inexpensive. Alternatively, a tracker produces a diverging laser beam which is directed towards a fixed position, and senses when a retroreflective target enters the fixed field of view. An optically bar coded target can be read by the tracker to provide information about the target. The target can be formed of a ball lens with a bar code on one end. As the target moves through the field, the ball lens causes the laser beam to scan across the bar code.

An improved lambda-scheme for one-dimensional flows

NASA Technical Reports Server (NTRS)

Moretti, G.; Dipiano, M. T.

1983-01-01

A code for the calculation of one-dimensional flows is presented, which combines a simple and efficient version of the lambda-scheme with tracking of discontinuities. The latter is needed to identify points where minor departures from the basic integration scheme are applied to prevent infiltration of numerical errors. Such a tracking is obtained via a systematic application of Boolean algebra. It is, therefore, very efficient. Fifteen examples are presented and discussed in detail. The results are exceptionally good. All discontinuites are captured within one mesh interval.

Analytical three-dimensional neutron transport benchmarks for verification of nuclear engineering codes. Final report

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

Ganapol, B.D.; Kornreich, D.E.

Because of the requirement of accountability and quality control in the scientific world, a demand for high-quality analytical benchmark calculations has arisen in the neutron transport community. The intent of these benchmarks is to provide a numerical standard to which production neutron transport codes may be compared in order to verify proper operation. The overall investigation as modified in the second year renewal application includes the following three primary tasks. Task 1 on two dimensional neutron transport is divided into (a) single medium searchlight problem (SLP) and (b) two-adjacent half-space SLP. Task 2 on three-dimensional neutron transport covers (a) pointmore » source in arbitrary geometry, (b) single medium SLP, and (c) two-adjacent half-space SLP. Task 3 on code verification, includes deterministic and probabilistic codes. The primary aim of the proposed investigation was to provide a suite of comprehensive two- and three-dimensional analytical benchmarks for neutron transport theory applications. This objective has been achieved. The suite of benchmarks in infinite media and the three-dimensional SLP are a relatively comprehensive set of one-group benchmarks for isotropically scattering media. Because of time and resource limitations, the extensions of the benchmarks to include multi-group and anisotropic scattering are not included here. Presently, however, enormous advances in the solution for the planar Green`s function in an anisotropically scattering medium have been made and will eventually be implemented in the two- and three-dimensional solutions considered under this grant. Of particular note in this work are the numerical results for the three-dimensional SLP, which have never before been presented. The results presented were made possible only because of the tremendous advances in computing power that have occurred during the past decade.« less

Solitary wave solutions of two-dimensional nonlinear Kadomtsev-Petviashvili dynamic equation in dust-acoustic plasmas

NASA Astrophysics Data System (ADS)

Seadawy, Aly R.

2017-09-01

Nonlinear two-dimensional Kadomtsev-Petviashvili (KP) equation governs the behaviour of nonlinear waves in dusty plasmas with variable dust charge and two temperature ions. By using the reductive perturbation method, the two-dimensional dust-acoustic solitary waves (DASWs) in unmagnetized cold plasma consisting of dust fluid, ions and electrons lead to a KP equation. We derived the solitary travelling wave solutions of the two-dimensional nonlinear KP equation by implementing sech-tanh, sinh-cosh, extended direct algebraic and fraction direct algebraic methods. We found the electrostatic field potential and electric field in the form travelling wave solutions for two-dimensional nonlinear KP equation. The solutions for the KP equation obtained by using these methods can be demonstrated precisely and efficiency. As an illustration, we used the readymade package of Mathematica program 10.1 to solve the original problem. These solutions are in good agreement with the analytical one.

Linearization instability for generic gravity in AdS spacetime

NASA Astrophysics Data System (ADS)

Altas, Emel; Tekin, Bayram

2018-01-01

In general relativity, perturbation theory about a background solution fails if the background spacetime has a Killing symmetry and a compact spacelike Cauchy surface. This failure, dubbed as linearization instability, shows itself as non-integrability of the perturbative infinitesimal deformation to a finite deformation of the background. Namely, the linearized field equations have spurious solutions which cannot be obtained from the linearization of exact solutions. In practice, one can show the failure of the linear perturbation theory by showing that a certain quadratic (integral) constraint on the linearized solutions is not satisfied. For non-compact Cauchy surfaces, the situation is different and for example, Minkowski space having a non-compact Cauchy surface, is linearization stable. Here we study, the linearization instability in generic metric theories of gravity where Einstein's theory is modified with additional curvature terms. We show that, unlike the case of general relativity, for modified theories even in the non-compact Cauchy surface cases, there are some theories which show linearization instability about their anti-de Sitter backgrounds. Recent D dimensional critical and three dimensional chiral gravity theories are two such examples. This observation sheds light on the paradoxical behavior of vanishing conserved charges (mass, angular momenta) for non-vacuum solutions, such as black holes, in these theories.

Effects of local defect growth in direct-drive cryogenic implosions on OMEGA

NASA Astrophysics Data System (ADS)

Igumenshchev, I. V.; Goncharov, V. N.; Shmayda, W. T.; Harding, D. R.; Sangster, T. C.; Meyerhofer, D. D.

2013-08-01

Spherically symmetric, low-adiabat (adiabat α ≲ 3) cryogenic direct-drive-implosion experiments on the OMEGA laser [T. R. Boehly et al., Opt. Commun. 133, 495 (1995)] yield less than 10% of the neutrons predicted in one-dimensional hydrodynamic simulations. Two-dimensional hydrodynamic simulations suggest that this performance degradation can be explained assuming perturbations from isolated defects of submicron to tens-of-micron scale on the outer surface or inside the shell of implosion targets. These defects develop during the cryogenic filling process and typically number from several tens up to hundreds for each target covering from about 0.2% to 1% of its surface. The simulations predict that such defects can significantly perturb the implosion and result in the injection of about 1 to 2 μg of the hot ablator (carbon-deuterium) and fuel (deuterium-tritium) materials from the ablation surface into the targets. Both the hot mass injection and perturbations of the shell reduce the final shell convergence ratio and implosion performance. The injected carbon ions radiatively cool the hot spot, reducing the fuel temperature, and further reducing the neutron yield. The negative effect of local defects can be minimized by decreasing the number and size of these defects and/or using more hydrodynamically stable implosion designs with higher shell adiabat.

Development of new two-dimensional spectral/spatial code based on dynamic cyclic shift code for OCDMA system

NASA Astrophysics Data System (ADS)

Jellali, Nabiha; Najjar, Monia; Ferchichi, Moez; Rezig, Houria

2017-07-01

In this paper, a new two-dimensional spectral/spatial codes family, named two dimensional dynamic cyclic shift codes (2D-DCS) is introduced. The 2D-DCS codes are derived from the dynamic cyclic shift code for the spectral and spatial coding. The proposed system can fully eliminate the multiple access interference (MAI) by using the MAI cancellation property. The effect of shot noise, phase-induced intensity noise and thermal noise are used to analyze the code performance. In comparison with existing two dimensional (2D) codes, such as 2D perfect difference (2D-PD), 2D Extended Enhanced Double Weight (2D-Extended-EDW) and 2D hybrid (2D-FCC/MDW) codes, the numerical results show that our proposed codes have the best performance. By keeping the same code length and increasing the spatial code, the performance of our 2D-DCS system is enhanced: it provides higher data rates while using lower transmitted power and a smaller spectral width.

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.

TWINTAN: A program for transonic wall interference assessment in two-dimensional wind tunnels

NASA Technical Reports Server (NTRS)

Kemp, W. B., Jr.

1980-01-01

A method for assessing the wall interference in transonic two dimensional wind tunnel test was developed and implemented in a computer program. The method involves three successive solutions of the transonic small disturbance potential equation to define the wind tunnel flow, the perturbation attriburable to the model, and the equivalent free air flow around the model. Input includes pressure distributions on the model and along the top and bottom tunnel walls which are used as boundary conditions for the wind tunnel flow. The wall induced perturbation fields is determined as the difference between the perturbation in the tunnel flow solution and the perturbation attributable to the model. The methodology used in the program is described and detailed descriptions of the computer program input and output are presented. Input and output for a sample case are given.

Superradiant instability of near extremal and extremal four-dimensional charged hairy black holes in anti-de Sitter spacetime

NASA Astrophysics Data System (ADS)

González, P. A.; Papantonopoulos, Eleftherios; Saavedra, Joel; Vásquez, Yerko

2017-03-01

We study the instability of near extremal and extremal four-dimensional anti-de Sitter charged hairy black holes to radial neutral massive and charged massless scalar field perturbations. We solve the scalar field equation by using the improved asymptotic iteration method and the time domain analysis, and we find the quasinormal frequencies. For the charged scalar perturbations, we find the superradiance condition by computing the reflection coefficient in the low-frequency limit, and we show that in the superradiance regime, which depends on the scalar hair charge, all modes of radial charged massless perturbations are unstable, indicating that the charged hairy black hole is superradiantly unstable. On the other hand, calculating the quasinormal frequencies of radial neutral scalar perturbations in this background, we find stability of the charged hairy black hole.

Impacts and environmental catastrophes: A study of the effects of impact events on the climate system

NASA Technical Reports Server (NTRS)

Pierazzo, E.

2005-01-01

The goal of this work is to investigate the perturbation of the climate system due to large impact events. Impacts are among the most important mechanisms for the evolution, distribution, and destruction of life in the universe. However, the possible climatic effects of an impact were not seriously considered until 1980, when Louis and Walter Alvarez suggested that the profound end-Cretaceous extinction might have been caused by the impact of an asteroid or comet about 10 km in diameter. Since then, the climatic change associated with the end-Cretaceous impact has become one of the most interesting and still unresolved questions in linking the well-known Chicxulub impact event and the end- Cretaceous mass extinction. While the end-Cretaceous impact offers the best-documented case of an impact affecting the Earth's climate and biota, even smaller (and more frequent in time) impacts could introduce significant perturbations of the climate comparable, if not larger, to the largest known volcanic perturbations. We propose to study the mechanical and thermal state of the atmosphere following an impact event. This will be done by using both one-dimensional and three-dimensional climate models. When necessary, modifications of the state-of-the-art general circulation models will b e carried out. We want to use the end-Cretaceous impact event as a case study. This allows us to take advantage of the extensive modeling of this impact event that has already been carried out through a previous Exobiology grant. Furthermore, a large experimental dataset, that can be used to constrain and test our models, is associated with the end-Cretaceous mass extinction (one of the largest of the Phanerozoic) and impact event.

Development of the PARVMEC Code for Rapid Analysis of 3D MHD Equilibrium

NASA Astrophysics Data System (ADS)

Seal, Sudip; Hirshman, Steven; Cianciosa, Mark; Wingen, Andreas; Unterberg, Ezekiel; Wilcox, Robert; ORNL Collaboration

2015-11-01

The VMEC three-dimensional (3D) MHD equilibrium has been used extensively for designing stellarator experiments and analyzing experimental data in such strongly 3D systems. Recent applications of VMEC include 2D systems such as tokamaks (in particular, the D3D experiment), where application of very small (delB/B ~ 10-3) 3D resonant magnetic field perturbations render the underlying assumption of axisymmetry invalid. In order to facilitate the rapid analysis of such equilibria (for example, for reconstruction purposes), we have undertaken the task of parallelizing the VMEC code (PARVMEC) to produce a scalable and temporally rapidly convergent equilibrium code for use on parallel distributed memory platforms. The parallelization task naturally splits into three distinct parts 1) radial surfaces in the fixed-boundary part of the calculation; 2) two 2D angular meshes needed to compute the Green's function integrals over the plasma boundary for the free-boundary part of the code; and 3) block tridiagonal matrix needed to compute the full (3D) pre-conditioner near the final equilibrium state. Preliminary results show that scalability is achieved for tasks 1 and 3, with task 2 still nearing completion. The impact of this work on the rapid reconstruction of D3D plasmas using PARVMEC in the V3FIT code will be discussed. Work supported by U.S. DOE under Contract DE-AC05-00OR22725 with UT-Battelle, LLC.

Perturbed redshifts from N -body simulations

NASA Astrophysics Data System (ADS)

Adamek, Julian

2018-01-01

In order to keep pace with the increasing data quality of astronomical surveys the observed source redshift has to be modeled beyond the well-known Doppler contribution. In this article I want to examine the gauge issue that is often glossed over when one assigns a perturbed redshift to simulated data generated with a Newtonian N -body code. A careful analysis reveals the presence of a correction term that has so far been neglected. It is roughly proportional to the observed length scale divided by the Hubble scale and therefore suppressed inside the horizon. However, on gigaparsec scales it can be comparable to the gravitational redshift and hence amounts to an important relativistic effect.

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.

A comparative study of Rayleigh-Taylor and Richtmyer-Meshkov instabilities in 2D and 3D in tantalum

NASA Astrophysics Data System (ADS)

Sternberger, Z.; Maddox, B. R.; Opachich, Y. P.; Wehrenberg, C. E.; Kraus, R. G.; Remington, B. A.; Randall, G. C.; Farrell, M.; Ravichandran, G.

2017-01-01

Driving a shock wave through the interface between two materials with different densities can result in the Richtmyer-Meshkov or Rayleigh-Taylor instability and initial perturbations at the interface will grow. If the shock wave is sufficiently strong, the instability will lead to plastic flow at the interface. Material strength will reduce the amount of plastic flow and suppress growth. While such instabilities have been investigated in 2D, no studies of this phenomena have been performed in 3D on materials with strength. Initial perturbations to seed the hydrodynamic instability were coined into tantalum recovery targets. Two types of perturbations were used, two dimensional (2D) perturbations (hill and valley) and three-dimensional (3D) perturbations (egg crate pattern). The targets were subjected to dynamic loading using the Janus laser at the Jupiter Laser Facility. Shock pressures ranged from 50 GPa up to 150 GPa and were calibrated using VISAR drive targets.

Two dimensional numerical prediction of deflagration-to-detonation transition in porous energetic materials.

PubMed

Narin, B; Ozyörük, Y; Ulas, A

2014-05-30

This paper describes a two-dimensional code developed for analyzing two-phase deflagration-to-detonation transition (DDT) phenomenon in granular, energetic, solid, explosive ingredients. The two-dimensional model is constructed in full two-phase, and based on a highly coupled system of partial differential equations involving basic flow conservation equations and some constitutive relations borrowed from some one-dimensional studies that appeared in open literature. The whole system is solved using an optimized high-order accurate, explicit, central-difference scheme with selective-filtering/shock capturing (SF-SC) technique, to augment central-diffencing and prevent excessive dispersion. The sources of the equations describing particle-gas interactions in terms of momentum and energy transfers make the equation system quite stiff, and hence its explicit integration difficult. To ease the difficulties, a time-split approach is used allowing higher time steps. In the paper, the physical model for the sources of the equation system is given for a typical explosive, and several numerical calculations are carried out to assess the developed code. Microscale intergranular and/or intragranular effects including pore collapse, sublimation, pyrolysis, etc. are not taken into account for ignition and growth, and a basic temperature switch is applied in calculations to control ignition in the explosive domain. Results for one-dimensional DDT phenomenon are in good agreement with experimental and computational results available in literature. A typical shaped-charge wave-shaper case study is also performed to test the two-dimensional features of the code and it is observed that results are in good agreement with those of commercial software. Copyright © 2014 Elsevier B.V. All rights reserved.

Impacts of DNAPL Source Treatment: Experimental and Modeling Assessment of the Benefits of Partial DNAPL Source Removal

DTIC Science & Technology

2009-09-01

nuclear industry for conducting performance assessment calculations. The analytical FORTRAN code for the DNAPL source function, REMChlor, was...project. The first was to apply existing deterministic codes , such as T2VOC and UTCHEM, to the DNAPL source zone to simulate the remediation processes...but describe the spatial variability of source zones unlike one-dimensional flow and transport codes that assume homogeneity. The Lagrangian models

Numerical investigation of stability of breather-type solutions of the nonlinear Schrödinger equation

NASA Astrophysics Data System (ADS)

Calini, A.; Schober, C. M.

2013-09-01

In this article we present the results of a broad numerical investigation on the stability of breather-type solutions of the nonlinear Schrödinger (NLS) equation, specifically the one- and two-mode breathers for an unstable plane wave, which are frequently used to model rogue waves. The numerical experiments involve large ensembles of perturbed initial data for six typical random perturbations. Ensemble estimates of the "closeness", A(t), of the perturbed solution to an element of the respective unperturbed family indicate that the only neutrally stable breathers are the ones of maximal dimension, that is: given an unstable background with N unstable modes, the only neutrally stable breathers are the N-dimensional ones (obtained as a superimposition of N simple breathers via iterated Backlund transformations). Conversely, breathers which are not fully saturated are sensitive to noisy environments and are unstable. Interestingly, A(t) is smallest for the coalesced two-mode breather indicating the coalesced case may be the most robust two-mode breather in a laboratory setting. The numerical simulations confirm and provide a realistic realization of the stability behavior established analytically by the authors.