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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

A new estimate of the chondrule cooling rate deduced from an analysis of compositional zoning of relict olivine

Compositional zoning in chondrule phenocrysts records the crystallization environments in the early solar nebula. We modeled the growth of olivine phenocrysts from a silicate melt and proposed a new fractional crystallization model that provides a relation between the zoning profile and the cooling rate. In our model, we took elemental partitioning at a growing solid-liquid interface and time-dependent solute diffusion in the liquid into consideration. We assumed a local equilibrium condition, namely, that the compositions at the interface are equal to the equilibrium ones at a given temperature. We carried out numerical simulations of the fractional crystallization in one-dimensional planar geometry. The simulations revealed that under a constant cooling rate the growth velocity increases exponentially with time and a linear zoning profile forms in the solid as a result. We derived analytic formulae of the zoning profile, which reproduced the numerical results for wide ranges of crystallization conditions. The formulae provide a useful tool to estimate the cooling rate from the compositional zoning. Applying the formulae to low-FeO relict olivine grains in type II porphyritic chondrules observed by Wasson and Rubin, we estimate the cooling rate to be ∼200-2000 K s{sup –1}, which is greater than that expected from furnace-based experiments by orders of magnitude. Appropriate solar nebula environments for such rapid cooling conditions are discussed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Mixed non-Gaussianity in multiple-DBI inflation

We study a model of multiple-field DBI inflation leading to mixed form of primordial non-Gaussianity, including equilateral and local bispectrum shapes. We present a general formalism based on the Hamilton-Jacobi approach, allowing us to go beyond slow-roll, combining the three-point function for the fields at Hubble-exit with the non-linear evolution of super-Hubble scales. We are able to obtain analytic results by taking a separable Ansatz for the Hubble rate. We find general expressions for both the equilateral and local type non-Gaussianity parameter f{sub NL}. The equilateral non-Gaussianity includes the usual enhancement for small sound speeds, but multiplied by an analytic factor which can lead to a suppression. We illustrate our results with two scenarios. In the first model, previously found to have detectable local non-Gaussianity, we find that the equilateral signal is not sufficiently suppressed to evade current observational bounds. In our second scenario we construct a model which exhibits both a detectable equilateral f{sub NL} and a negative local f{sub NL}.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Describing massive neutrinos in cosmology as a collection of independent flows

A new analytical approach allowing to account for massive neutrinos in the non-linear description of the growth of the large-scale structure of the universe is proposed. Unlike the standard approach in which neutrinos are described as a unique hot fluid, it is shown that the overall neutrino fluid can be equivalently decomposed into a collection of independent flows. Starting either from elementary conservation equations or from the evolution equation of the phase-space distribution function and assuming that there is no shell-crossing, we derive the two non-linear motion equations that each of these flows satisfies. Those fluid equations describe the evolution of macroscopic fields. We explain in detail the connection between the collection of flows we defined and the standard massive neutrino fluid. Then, in the particular case of adiabatic initial conditions, we explicitly check that, at linear order, the resolution of this new system of equations reproduces the results obtained in the standard approach based on the collisionless Boltzmann hierarchy. Besides, the approach advocated in this paper allows to show how each neutrino flow settles into the cold dark matter flow depending on initial velocities. Although valid up to shell-crossing only, it is a further step towards a fully non-linear treatment of the dynamical evolution of neutrinos in the framework of large-scale structure growth.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2011

Tradeoff between energy and error in the discrimination of quantum-optical devices

We address the problem of energy-error tradeoff in the discrimination between two linear passive quantum optical devices with a single use. We provide an analytical derivation of the optimal strategy for beamsplitters and an iterative algorithm converging to the optimum in the general case. We then compare the optimal strategy with a simpler strategy using coherent input states and homodyne detection. It turns out that the former requires much less energy in order to achieve the same performances.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Magnon damping in the zigzag phase of the Kitaev-Heisenberg- Γ model on a honeycomb lattice

We calculate dispersions and damping of the magnetic excitations in the model that is relevant to a description of the magnetic properties of iridium oxides α-Li 2 IrO 3 and Na 2 IrO 3 , and Ru-based materials such as α-RuCl 3 . Focusing on the parameter regime with the zigzag spin pattern in the ground state that is consistent with experiments, we are able to obtain explicit analytical expressions for magnon energies and eigenstates and go beyond the standard linear spin-wave theory approximation by calculating magnon damping and demonstrating its role in the dynamical structure factor. Here, we show that the magnon damping effects in both Born and self-consistent approximations are very significant, underscoring the importance of nonlinear magnon coupling in interpreting broad features in the neutron-scattering spectra.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Distinct universality classes of diffusive transport from full counting statistics

The hydrodynamic transport of local conserved densities furnishes an effective coarse-grained description of the dynamics of a many-body quantum system. However, the full quantum dynamics contains much more structure beyond the simplified hydrodynamic description. Here we show that systems with the same hydrodynamics can nevertheless belong to distinct dynamical universality classes, as revealed by new classes of experimental observables accessible in synthetic quantum systems, which can, for instance, measure simultaneous site-resolved snapshots of all of the particles in a system. Specifically, we study the full counting statistics of spin transport, whose first moment is related to linear-response transport, but the higher moments go beyond. We present an analytic theory of the full counting statistics of spin transport in various integrable and nonintegrable anisotropic one-dimensional spin models, including the XXZ spin chain. We find that spin transport, while diffusive on average, is governed by a distinct non-Gaussian dynamical universality class in the models considered. We consider a setup in which the left and right half of the chain are initially created at different magnetization densities, and consider the probability distribution of the magnetization transferred between the two half-chains. We derive a closed-form expression for the probability distribution of the magnetization transfer, in terms of random walks on the half-line. We show that this distribution strongly violates the large-deviation form expected for diffusive chaotic systems, and explain the physical origin of this violation. Here, we discuss the crossovers that occur as the initial state is brought closer to global equilibrium. Our predictions can directly be tested in experiments using quantum gas microscopes or superconducting qubit arrays.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2011

Linearized f(R) gravity: Gravitational radiation and Solar System tests

We investigate the linearized form of metric f(R)-gravity, assuming that f(R) is analytic about R=0 so it may be expanded as f(R)=R+a{sub 2}R{sup 2}/2+.... Gravitational radiation is modified, admitting an extra mode of oscillation, that of the Ricci scalar. We derive an effective energy-momentum tensor for the radiation. We also present weak-field metrics for simple sources. These are distinct from the equivalent Kerr (or Schwarzschild) forms. We apply the metrics to tests that could constrain f(R). We show that light deflection experiments cannot distinguish f(R)-gravity from general relativity as both have an effective post-Newtonian parameter {gamma}=1. We find that planetary precession rates are enhanced relative to general relativity; from the orbit of Mercury we derive the bound |a{sub 2}| < or approx. 1.2x10{sup 18} m{sup 2}. Gravitational-wave astronomy may be more useful: considering the phase of a gravitational waveform we estimate deviations from general relativity could be measurable for an extreme-mass-ratio inspiral about a 10{sup 6}M{sub {center_dot}} black hole if |a{sub 2}| > or approx. 10{sup 17}m{sup 2}, assuming that the weak-field metric of the black hole coincides with that of a point mass. However Eoet-Wash experiments provide the strictest bound |a{sub 2}| < or approx. 2x10{sup -9} m{sup 2}. Although the astronomical bounds are weaker, they are still of interest in the case that the effective form of f(R) is modified in different regions, perhaps through the chameleon mechanism. Assuming the laboratory bound is universal, we conclude that the propagating Ricci scalar mode cannot be excited by astrophysical sources.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

A Flexible-blocking Based Approach for Performance Tuning of Matrix Multiplication Routines for Large Matrices with Edge Cases

Efficient and scalable matrix operations are being highly demanding in the recent era of Machine Learning, Deep Learning, and Big Data Analytics. The two commonly used matrix-matrix operations in the Basic Linear Algebra Subprograms (BLAS) specification are General Matrix-Matrix multiplication (GEMM) and Symmetric Rank-k update (SYRK). The SYRK routine is a specialization of the GEMM routine, where half of the multiplications are skipped as the resultant matrix is known to be symmetric. Fortunately, several linear algebra libraries implement these BLAS routines quite efficiently. The libraries usually partition the input matrices into blocks and place them in processor caches, thus improving performance by leveraging the caches. However, the contemporary libraries are highly optimized for squarish matrices, but the performance degrades significantly for the matrices with edge case (strictly thin or strictly fat shapes) in the multicore machine. The primary reason is that the current state-of-the-art libraries make fixed block shapes based on a processor architecture, and do not consider the shape of the input matrices. In this paper, we propose a new blocking approach, we name it Flexible-blocking, to mitigate the scalability issues. In contrast to the contemporary libraries, our approach formulates the blocks of the input matrices based on the shapes of the matrices as well as the number of threads used in the implementation. Our proposed technique shows noticeable performance improvement on multicore shared-memory machines for the edge case matrices.

Hossain, Md Mosharaf↗
DOE OSTI2025

Grey-Box System Identification of Grid-Forming Inverters

This paper demonstrates the use of grey-box system identification methods for simplifying and understanding the nonlinear power dynamics of grid-forming inverters (GFMs). The power and frequency outputs of complex high-order GFM models are fed into system identification software in order to fit them to a predetermined LTI system and learn system parameters such as (synthetic) inertia and droop constants. The same process is then run for a high-order synchronous generator model, and the outputs are fit to the same set of LTI equations. Simulation of a network of GFM inverters with diverse control architecture is also performed for the same process. The intent is threefold: first, to demonstrate the appropriateness of unified LTI models for describing the power and frequency dynamics of individual resources and connected networks, in order to facilitate analysis of larger heterogeneous networked systems; second, to discover the relationship between internal control parameters of GFMs and their externally observed values; and third, to validate that grey-box data-driven system identification techniques can be a valuable tool to discover the values of important parameters in the absence of explicit vendor models.

analytical models↗
DOE OSTI2016

Electromagnetic-field amplification in finite one-dimensional photonic crystals

The electromagnetic-field distribution in a finite one-dimensional photonic crystal is studied using the numerical solution of Maxwell’s equations by the transfer-matrix method. The dependence of the transmission coefficient T on the period d (or the wavelength λ) has the characteristic form with M–1 (M is the number of periods in the structure) maxima with T = 1 in the allowed band of an infinite crystal and zero values in the forbidden band. The field-modulus distribution E(x) in the structure for parameters that correspond to the transmission maxima closest to the boundaries of forbidden bands has maxima at the center of the structure; the value at the maximum considerably exceeds the incident-field strength. For the number of periods M ~ 50, more than an order of magnitude increase in the field amplification is observed. The numerical results are interpreted with an analytic theory constructed by representing the solution in the form of a linear combination of counterpropagating Floquet modes in a periodic structure.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Parametric dependencies of resonant layer responses across linear, two-fluid, drift-MHD regimes

Non-axisymmetric magnetic fields arising in a tokamak either by external or internal perturbations can induce complex non-ideal MHD responses in their resonant surfaces while remaining ideally evolved elsewhere. This layer response can be characterized in a linear regime by a single parameter called the inner-layer Delta, which enables outer-layer matching and the prediction of torque balance to non-linear island regimes. Here, we follow strictly one of the most comprehensive analytic treatments including two-fluid and drift MHD effects and keep the fidelity of the formulation by incorporating the numerical method based on the Riccati transformation when quantifying the inner-layer Delta. The proposed scheme reproduces not only the predicted responses in essentially all asymptotic regimes but also with continuous transitions as well as improved accuracies. In particular, the Delta variations across the inertial regimes with viscous or semi-collisional effects have been further resolved, in comparison with additional analytic solutions. The results imply greater shielding of the electromagnetic torque at the layer than what would be expected by earlier work when the viscous or semi-collisional effects can compete against the inertial effects, and also due to the intermediate regulation by kinetic Alfven wave resonances as rotation slows down. These are important features that can alter the nonaxisymmetric plasma responses including the field penetration by external fields or island seeding process in rotating tokamak plasmas.

MHD instability↗
DOE OSTI2022

kynema-fmb [SWR-23-07]

Kynema-FMB (FKA: Kynema) is an open-source performance portable flexible multibody (FMB) dynamics solver designed for time-domain simulations. While originally tailored for wind turbine structural dynamics, the formulation and implementation are those of a general flexible-multidbody dynamics solver that can readily be applied to a wide range of systems. Kynema was designed with a narrow focus, namely to provide a lightweight, fast, accurate FMD solver for coupling to computational-fluid-dynamics (CFD) codes, especially the CFD codes in the Kynema suite, for fluid-structure-interaction (FSI) simulations. Kynema-FMB is equipped to model systems that can be represented as a collection of beams and rigid bodies that are connected through constraints. Degrees of freedom are defined in the inertial/global frame of reference and include displacements and rotations (formally as rotation matrices, but stored as quaternions). The underlying formulation is built on a Lie-group time integrator designed for index-3 differential-algebraic equations, which is second-order accurate in time (Bruls et al., 2012). Beam models are based on geometrically exact beam theory and are discretized as high-order spectral finite elements similar to those in BeamDyn (Wang et al., 2017). The governing equations for a FMD system like a wind turbine constitute a highly nonlinear system of constrained partial-differential equations. Kynema-FMB uses analytical Jacobians in the nonlinear-system solves in each time step. Linear systems use sparse storage and several third-party sparse-linear-system solvers are enabled. Ill conditioning of linear systems is mitigated with preconditioning described in Bottasso et al, 2008. Kynema-FMB is integrated with a simple open-source controller (ROSCO). There is an application programming interface (API) for coupling to geometry-resolved CFD (like that in Sharma et al., 2023) and actuator-force CFD (like that in Kuhn et al., 2025). In the latter, for actuator-line models, Kynema-FMB includes an internal blade-element solver that depends on user-provided lookup tables for coefficients of lift and drag, i.e., aerodynamic polars. Kynema-FMB is written in C++ and leverages Kokkos and Kokkos-Kernels (KokkosEcosystem) as its performance portability layer enabling simulations on both CPU and GPU systems. The repository is equipped with extensive automated testing at the unit and regression/system levels. The following describes the high-level development objectives conceived for Kynema: *Kynema will follow modern software development best practices, including test-driven development (TDD), version control, hierarchical automated testing, and continuous integration (CI) for a robust development environment. *The core data structures are memory efficient and enable vectorization and parallelization at multiple levels. *Data structures are data-oriented to exploit methods for accelerated computing including high utilization of chip resources (e.g., single instruction multiple data (SIMD) instruction sets) and parallelization using GP-GPUs. *The computational algorithms incorporate robust open-source libraries for mathematical operations, resource allocation, and data management. *The API design considers multiple stakeholder needs and ensure integration with existing and future ecosystems for data science, machine learning, and AI. *Kynema-FMB is written in modern C++ and leverages Kokkos as its performance-portability library with inspiration from the kynema stack.

Sprague, MichaelA.↗
DOE OSTI2020

Collaborative Research: The Evolution of Magnetic Skeletons During 3D Reconnection. Final report

We investigated how the plasmoid instability affects the evolution of topological features of a magnetic field and the process of magnetic reconnection with three sub-projects. (i) We extensively studied the criterion for plasmoid-mediated current sheet disruption and the onset of fast reconnection. Scaling relations with respect to plasma parameters have been established, using both direct numerical simulations and analytic theories, for the disruption time, the current sheet width, the instability linear growth rate, and the dominant mode wavenumber. (ii) We applied the theories of plasmoid-mediated current sheet disruption to magnetohydrodynamic turbulence and predicted a steepening of the turbulence energy spectrum. This steepening has been confirmed by direct numerical simulations. (iii) We performed numerical simulations of ultraviolet burst events in the solar atmosphere with a spine-dome magnetic topology. The plasmoid instability leads to unsteady reconnection, which is consistent with solar observations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Observable-projected ensembles

Measurements in many-body quantum systems can generate non-trivial phenomena, such as preparation of long-range entangled states, dynamical phase transitions, or measurement-altered criticality. Here, we introduce a new measurement scheme that produces an ensemble of mixed states in a subsystem, obtained by measuring a local Hermitian observable on part of its complement. We refer to this as the observable-projected ensemble . Unlike standard projected ensembles-where pure states are generated by projective measurements on the complement-our approach involves projective partial measurements of specific observables. This setup has two main advantages: theoretically, it is amenable to analytical computations, especially within conformal field theories. Experimentally, it requires only a linear number of measurements, rather than an exponential one, to probe the properties of the ensemble. As a first step in exploring the observable-projected ensemble, we investigate its entanglement properties in conformal field theory and perform a detailed analysis of the free compact boson.

Milekhin, Alexey [California Institute of Technolo↗
NASA NTRS1972

Method of characteristics for supersonic flow of a guiding centre plasma.

A study of the compressive magneto-acousic waves in a guiding centre plasma shows that the wavefront that emerges point disturbance after a finite time is a simple oblate spheroid with the axis of revolution parallel to the field lines. Thus, in a steady three-dimensional supersonic flow of guiding centre plasma a simple analytic exprencan be obtained to represent the characteristic surfaces. From a proper linear combination of the governing macroscopic equations, the characteristic equation is obtained. It represents the propagation of disturbances on the characteristic surface. The characteristic theory can be used to study the interaction of the solar wind with the moon and possibly with other planetary bodies.

Whang, Y. C.↗
NASA NTRS1973

Research in plasma physics

Three aspects of barium ion cloud dynamics are discussed. First, the effect of the ratio of ion cloud conductivity to background ionospheric conductivity on the motion of barium ion clouds is investigated and compared with observations of barium ion clouds. This study led to the suggestion that the conjugate ionosphere participates in the dynamics of barium ion clouds. Second, analytic work on the deformation of ion clouds is presented. Third, a linearized stability theory was extended to include the effect of the finite extent of an ion cloud, as well as the effect of the ratio of ion cloud to ionospheric conductivities. The stability properties of a plasma with contra-streaming ion beams parallel to a magnetic field are investigated. The results are interpreted in terms of parameters appropriate for collisionless shock waves. It is found that this particular instability can be operative only if the up-stream Alfven Mach number exceeds 5.5.

Source record↗
NASA NTRS1974

Comparison of approximate and numerical analyses of nonlinear combustion instability

At the present time, there are three general analytical techniques available to study problems of unsteady motions in rocket motors: linear stability analysis; approximate nonlinear analysis, founded on examining the behavior of coupled normal modes; and numerical calculations based on the conservation equations for one-dimensional flows. The last two yield the linear results as a limit. It is the main purpose of this paper to check the accuracy of the approximate analysis against the numerical analysis for some special cases. The results provide some justification for using the approximate analysis to study three-dimensional problems.

Culick, F. E. C.↗
NASA NTRS1980

Acoustic propagation in rigid three-dimensional waveguides

The linear acoustic propagation in finite rigid three-dimensional waveguides is determined analytically using an eigenfunction expansion of the Helmholtz equation. The geometry considered consists of straight and circular bends of rectangular cross section with continuous interfaces (branches and sharp corners are excluded). The phenomena of resonance shift and relocation are explained for a bend-straight duct combination.

El-Raheb, M.↗