Search NASA⌕ Search

SEARCH · Search NASA

Results for “Analytic linearization”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2012

PHURBAS: AN ADAPTIVE, LAGRANGIAN, MESHLESS, MAGNETOHYDRODYNAMICS CODE. II. IMPLEMENTATION AND TESTS

We present an algorithm for simulating the equations of ideal magnetohydrodynamics and other systems of differential equations on an unstructured set of points represented by sample particles. The particles move with the fluid, so the time step is not limited by the Eulerian Courant-Friedrichs-Lewy condition. Full spatial adaptivity is required to ensure the particles fill the computational volume and gives the algorithm substantial flexibility and power. A target resolution is specified for each point in space, with particles being added and deleted as needed to meet this target. We have parallelized the code by adapting the framework provided by GADGET-2. A set of standard test problems, including 10{sup -6} amplitude linear magnetohydrodynamics waves, magnetized shock tubes, and Kelvin-Helmholtz instabilities is presented. Finally, we demonstrate good agreement with analytic predictions of linear growth rates for magnetorotational instability in a cylindrical geometry. This paper documents the Phurbas algorithm as implemented in Phurbas version 1.1.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2011

Linearized T-Matrix and Mie Scattering Computations

We present a new linearization of T-Matrix and Mie computations for light scattering by non-spherical and spherical particles, respectively. In addition to the usual extinction and scattering cross-sections and the scattering matrix outputs, the linearized models will generate analytical derivatives of these optical properties with respect to the real and imaginary parts of the particle refractive index, and (for non-spherical scatterers) with respect to the ''shape'' parameter (the spheroid aspect ratio, cylinder diameter/height ratio, Chebyshev particle deformation factor). These derivatives are based on the essential linearity of Maxwell's theory. Analytical derivatives are also available for polydisperse particle size distribution parameters such as the mode radius. The T-matrix formulation is based on the NASA Goddard Institute for Space Studies FORTRAN 77 code developed in the 1990s. The linearized scattering codes presented here are in FORTRAN 90 and will be made publicly available.

Analytic linearization↗
NASA NTRS1982

Non-linear inplane vibrations of wrinkling films

This paper presents analytical investigations of non-linear effects caused by wrinkling of films in the course of contraction and shearing. As shown, such non-linearities are important even for very small amplitudes of vibrations. Approximate analytical solutions for inplane vibrations are derived.

Zak, M.↗
NASA NTRS1974

Methods studies toward simplified rotor-body dynamics, part 1

This report is directed to the problem of developing an adequate but not overly complex linear flight dynamics analytical model of a rotorcraft to study stability, control, gust and random turbulence responses. Since the conventional flight dynamics analysis using quasisteady rotor derivatives is adequate for the long period modes like the phugoid mode, only short time responses are considered here, where rotor-body coupling is of importance. Thus the body motion consists of pitch, roll and vertical motion, omitting linear longitudinal and lateral and yaw perturbations. Five analytical models of varying degree of sophistication are applied to a hypothetical hingeless compound helicopter operating up to .8 rotor advance ratio. Stability and response data are obtained for the basic helicopter and for the vehicle with two simple control feedback systems.

Hohenemser, K. H.↗
DOE OSTI2011

Commissioning of Photon Beams of a Flattening Filter-Free Linear Accelerator and the Accuracy of Beam Modeling Using an Anisotropic Analytical Algorithm;TrueBeam; Flattening filter free; Commissioning; Anisotropic analytical algorithm

Purpose: To investigate dosimetric characteristics of a new linear accelerator designed to deliver flattened, as well as flattening filter-free (FFF), beams. To evaluate the accuracy of beam modeling under physical conditions using an anisotropic analytical algorithm. Methods and Materials: Dosimetric data including depth dose curves, profiles, surface dose, penumbra, out-of-field dose, output, total and scatter factors were examined for four beams (X6, X6FFF, X10, and X10FFF) of Varian's TrueBeam machine. Beams modeled by anisotropic analytical algorithm were compared with measured dataset. Results: FFF beams have lower mean energy (tissue-phantom ratio at the depths of 20 and 10 cm (TPR 20/10): X6, 0.667; X6FFF, 0.631; X10, 0.738; X10FFF, 0.692); maximum dose is located closer to the surface; and surface dose increases by 10%. FFF profiles have sharper but faster diverging penumbra. For small fields and shallow depths, dose outside a field is lower for FFF beams; however, the advantage fades with increasing phantom scatter. Output increases 2.26 times for X6FFF and 4.03 times for X10FFF and is less variable with field size; collimator exchange effect is reduced. A good agreement between modeled and measured data is observed. Criteria of 2% depth-dose and 2-mm distance-to-agreement are always met. Conclusion: Reference dosimetric characteristics of TrueBeam photon bundles were obtained, and successful modeling of the beams was achieved.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
NASA NTRS1986

The 15 meter hoop-column antenna dynamics: Test and analysis

A 15 meter model of the hoop-column antenna concept has been vibration tested for model characterization and analytical model verification. Linear finite element analysis predicted the global vibration frequencies accurately. Good agreement between analysis and test data was obtained only after the analytical model was refined using static test data. As structures become more flexible, structural properties determined from static data become more accurate and should be used to update analytical models. Global vibration modes are not significantly affected by the surface mesh which permits simplified analytical models to be used for prediction of global behavior. These reduced models are believed sufficient for preliminary design and controls simulations where only global behavior is desired. The mesh modes were highly damped due to the knit mesh used for the reflector surface. These modes were also highly coupled and very difficult to measure in the laboratory. The inability to fully characterize the antenna mesh modes in the laboratory indicates robust methods for active surface vibration suppression will be needed. Fortunately, the surface mesh exhibits high passive damping which should be beneficial to active control systems.

Belvin, W. Keith↗
NASA NTRS1991

Generalized hydrodynamic transport in lattice-gas automata

The generalized hydrodynamics of two-dimensional lattice-gas automata is solved analytically in the linearized Boltzmann approximation. The dependence of the transport coefficients (kinematic viscosity, bulk viscosity, and sound speed) upon wave number k is obtained analytically. Anisotropy of these coefficients due to the lattice symmetry is studied for the entire range of wave number, k. Boundary effects due to a finite mean free path (Knudsen layer) are analyzed, and accurate comparisons are made with lattice-gas simulations.

Luo, Li-Shi↗
NASA NTRS1991

A methodology for using nonlinear aerodynamics in aeroservoelastic analysis and design

A methodology is presented for using the Volterra-Wiener theory of nonlinear systems in aeroservoelastic (ASE) analyses and design. The theory is applied to the development of nonlinear aerodynamic response models that can be defined in state-space form and are, therefore, appropriate for use in modern control theory. The theory relies on the identification of nonlinear kernels that can be used to predict the response of a nonlinear system due to an arbitrary input. A numerical kernel identification technique, based on unit impulse responses, is presented and applied to a simple bilinear, single-input single-output (SISO) system. The linear kernel (unit impulse response) and the nonlinear second-order kernel of the system are numerically-identified and compared with the exact, analytically-defined and linear and second-order kernels. This kernel identification technique is then applied to the CAP-TSD (Computational Aeroelasticity Program-Transonic Small Disturbance) code for identification of the linear and second-order kernels of a NACA64A010 rectangular wing undergoing pitch at M = 0.5, M = 8.5 (transonic), and M = 0.93 (transonic). Results presented demonstrate the feasibility of this approach for use with nonlinear, unsteady aerodynamic responses.

Silva, Walter A.↗
NASA NTRS1980

Practical implementation of the double linear damage rule and damage curve approach for treating cumulative fatigue damage

Simple procedures are presented for treating cumulative fatigue damage under complex loading history using either the damage curve concept or the double linear damage rule. A single equation is provided for use with the damage curve approach; each loading event providing a fraction of damage until failure is presumed to occur when the damage sum becomes unity. For the double linear damage rule, analytical expressions are provided for determining the two phases of life. The procedure involves two steps, each similar to the conventional application of the commonly used linear damage rule. When the sum of cycle ratios based on phase 1 lives reaches unity, phase 1 is presumed complete, and further loadings are summed as cycle ratios on phase 2 lives. When the phase 2 sum reaches unity, failure is presumed to occur. No other physical properties or material constants than those normally used in a conventional linear damage rule analysis are required for application of either of the two cumulative damage methods described. Illustrations and comparisons of both methods are discussed.

Manson, S. S.↗
NASA NTRS1981

Practical implementation of the double linear damage rule and damage curve approach for treating cumulative fatigue damage

Simple procedures are given for treating cumulative fatigue damage under complex loading history using either the damage curve concept or the double linear damage rule. A single equation is given for use with the damage curve approach; each loading event providing a fraction of damage until failure is presumed to occur when the damage sum becomes unity. For the double linear damage rule, analytical expressions are given for determining the two phases of life. The procedure comprises two steps, each similar to the conventional application of the commonly used linear damage rule. Once the sum of cycle ratios based on Phase I lives reaches unity, Phase I is presumed complete, and further loadings are summed as cycle ratios based on Phase II lives. When the Phase II sum attains unity, failure is presumed to occur. It is noted that no physical properties or material constants other than those normally used in a conventional linear damage rule analysis are required for application of either of the two cumulative damage methods described. Illustrations and comparisons are discussed for both methods.

Manson, S. S.↗
DOE OSTI2024

Analytic Nuclear Gradients for Complete Active Space Linearized Pair-Density Functional Theory

Accurately modeling photochemical reactions is difficult due to the presence of conical intersections and locally avoided crossings, as well as the inherently multiconfigurational character of excited states. As such, one needs a multistate method that incorporates state interaction in order to accurately model the potential energy surface at all nuclear coordinates. The recently developed linearized pair-density functional theory (L-PDFT) is a multistate extension of multiconfiguration PDFT, and it has been shown to be a cost-effective post-MCSCF method (as compared to more traditional and expensive multireference many-body perturbation methods or multireference configuration interaction methods) that can accurately model potential energy surfaces in regions of strong nuclear–electronic coupling in addition to accurately predicting Franck–Condon vertical excitations. Here, in this paper, we report the derivation of analytic gradients for L-PDFT and their implementation in the PySCF-forge software, and we illustrate the utility of these gradients for predicting ground- and excited-state equilibrium geometries and adiabatic excitation energies for formaldehyde, s-trans-butadiene, phenol, and cytosine.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1985

Construction of approximate analytical solutions to a new class of non-linear oscillator equation

The principle of harmonic balance is invoked in the development of an approximate analytic model for a class of nonlinear oscillators typified by a mass attached to a stretched wire. By assuming that harmonic balance will hold, solutions are devised for a steady state limit cycle and/or limit point motion. A method of slowly varying amplitudes then allows derivation of approximate solutions by determining the form of the exact solutions and substituting into them the lowest order terms of their respective Fourier expansions. The latter technique is actually a generalization of the method proposed by Kryloff and Bogoliuboff (1943).

Mickens, R. E.↗
NASA NTRS1994

A note on the accuracy of spectral method applied to nonlinear conservation laws

Fourier spectral method can achieve exponential accuracy both on the approximation level and for solving partial differential equations if the solutions are analytic. For a linear partial differential equation with a discontinuous solution, Fourier spectral method produces poor point-wise accuracy without post-processing, but still maintains exponential accuracy for all moments against analytic functions. In this note we assess the accuracy of Fourier spectral method applied to nonlinear conservation laws through a numerical case study. We find that the moments with respect to analytic functions are no longer very accurate. However the numerical solution does contain accurate information which can be extracted by a post-processing based on Gegenbauer polynomials.

Shu, Chi-Wang↗
DOE OSTI2010

Discrimination of the effects of saturation and optical pumping in velocity-dependent pump-probe spectroscopy of rubidium: A simple analytical study

This paper presents a simple analytical theory for the velocity-dependent pump-probe laser spectroscopy of {sup 87}Rb and {sup 85}Rb atoms where the pump and the probe beams are circularly or linearly polarized. The analytical solutions of the line shapes of the velocity-selective optical pumping spectroscopy [G. Moon and H. R. Noh, Phys. Rev. A 78, 032506 (2008)] and saturated absorption spectroscopy [G. Moon and H. R. Noh, J. Opt. Soc. Am. B 25, 701 (2008); 27, 1741 (2010)] obtained in the previous reports, expressed as a sum of several Lorentzian functions, could be approximated as one (or in some cases, two) Lorentzian function(s). In particular, the contributions of the saturation and optical pumping effects could be discriminated explicitly in these simple analytical solutions, which is not possible in existing theories such as Nakayama's model. The simple analytical results for the saturation spectroscopy were compared with experimental results, and good agreement between them was observed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2011

Quantum interference in laser-assisted photoionization and analytical methods for the measurement of an attosecond xuv pulse

Investigations of the quantum interference in laser-assisted photoionization by an attosecond extreme ultraviolet (xuv) pulse shows an approximately constant value for the total photoionizations for different laser intensities. The square of the full width at half maximum of a photoelectron energy spectrum (PES) linearly depends on the laser intensity. By determining the laser-related phase of each streaked electron and using a transfer equation with linear corrections, an analytically quick method is proposed for precisely reconstructing the xuv pulse intensity (chirp) from one (two) measured PES(s) with a theoretical root-mean-square temporal (energy) difference of less than 1 attosecond (0.1 eV).

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS1988

Light propagation in inhomogeneous universes

Refsdal's (1970) method is generalized to study the propagation of light rays through an inhomogeneous universe. The probability distribution for the linear component of the cumulative shear (CS) along light rays is derived, and it is shown that the CS can be dominated by nonlinear components, espcially for light rays in empty cones. The amplification tail of the amplification probability distribution is compared with analytic results; these linear investigations are shown to underestimate the high-amplification probability and hence the importance of the amplification bias in source counts. The distribution of the ellipticity of images of infinitesimal circular sources is derived, and it is shown that this can be dominated by the nonlinear contributions to the CS.

Schneider, Peter↗
NASA NTRS1993

Linear theory of fast reconnection at an X-type neutral point

A complete analytic theory of dynamic linear reconnection at an X-type neutral point is developed. An eigenmode analysis, using cylindrical coordinates centered on the neutral point, extends the work of Craig & McClymont (1991) to include nonazimuthally symmetric perturbations on the flux function. It is demonstrated that all physically significant disturbances, both reconnective and nonreconnective, decay resistively on a 'fast' time scale. The significance of the linear theory is discussed within the general context of steady state and dynamic reconnection studies. Disturbances can be expected to focus explosively in the vicinity of the neutral point. This suggests the formation of a 'flux pile-up' current layer in which the bulk of the magnetic energy is released as heat rather than kinetic energy of mass motion.

Craig, I. J. D.↗
DOE OSTI2020

Evaluation of Graph Analytics Frameworks Using the GAP Benchmark Suite

The analysis of connected data is an increasingly important application in high-performance computing. Such analyses can reveal fraudulent patterns in financial transactions, optimize telecommunications networks, predict information flow in social networks, etc. However, the landscape of graph analytics is highly diverse. Graph algorithms stress processor architectures differently, and no one graph can represent all topologies. Consequently, no single approach or framework is expected to be optimal for all graph analytics problems. To help make sense of this diverse landscape, we evaluated four approaches to graph analytics: GraphBLAS, Galois, BGL17, GraphIt; and compare them against hand-tuned implementations that take advantage of hardware features on our test platform. Graph- BLAS formulates graph analytics as sparse linear algebra. Galois provides syntactic constructs for data parallelism over irregular data structures. BGL17 is a generic C++ template library for implementing graph algorithms. GraphIt provides a domain- specific language to describe and optimize graph algorithms. We use the GAP Benchmark Suite to establish baseline performance and guide the side-by-side evaluation of each framework. GAP consists of 30 tests: six graph analytics algorithms (breadth- first search, single-source shortest path, PageRank, betweenness centrality, connected components, and triangle counting) run on five graphs, each with different topological characteristics (e.g., high diameter, skewed degree distribution, high average degree). High-performance reference implementations are included for each benchmark algorithm. Because a graph can be loaded into memory a number of ways (e.g., flat file on disk, compressed sparse format, data frames, retrieved from SQL or NoSQL databases), our evaluation focused on computational performance rather than I/O. Our results show the relative strengths of each framework.

Graph algorithms, Benchmarking, shared-memory prog↗