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.

SourceResearch connectionFrom yearThrough year

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

At least 235 records · Page 13

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Quasi-linear theory of forced magnetic reconnection for the transition from the linear to the Rutherford regime

Abstract Using the in-viscid two-field reduced magneto-hydrodynamic (MHD) model, a new analytical theory is developed to unify the Hahm–Kulsrud–Taylor (HKT) linear solution and the Rutherford quasi-linear regime. Adopting a quasi-linear approach, we obtain a closed system of equations for plasma response in the Taylor problem. An integral form of an analytical solution is obtained for the forced magnetic reconnection, uniformly valid throughout the entire regimes from the HKT linear solution to the Rutherford quasi-linear solution. In particular, the quasi-linear effect can be described by a single coefficient K s ∝ S 8 / 5 ψ c 2 , where S = τ R τ A and ψ c are the Lundquist number and amplitude of external magnetic perturbation, respectively. The HKT linear solution for the response can be recovered when the index K s → 0. On the other hand, the quasi-linear effect plays a key role in the island growth when K s ~ 1. Our new analytical solution has also been compared with reduced MHD simulations with agreement.

Physics↗

DOE OSTI2020

Temperature scaling behavior of the linear magnetoresistance observed in high-temperature superconductors

An analytical model invoking variations in the charge-carrier density is used to generate magnetoresistance curves that are almost indistinguishable from those produced by sophisticated numerical models. This demonstrates that, though disorder is pivotal in causing linear magnetoresistance, the form of the magnetoresistance thus generated is insensitive to details of the disorder. Additionally, taken in conjunction with the temperature (T) dependence of the zero-field resistivity, realistic levels of disorder are shown to be sufficient to explain the linear magnetoresistance and field-T resistance scaling observed in high-temperature pnictide and cuprate superconductors. Hence, though the T-linear zero-field resistance is a definite signature of the “strange metal” state of high-temperature superconductors, their linear magnetoresistance and its scaling is unlikely to be so.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

Non-linear dynamics of jet quenching

We develop a comprehensive analytic framework for jet quenching in QCD media, based on a medium-induced parton cascade sourced by collinear virtual splittings. We show that the energy flow out of the jet cone, driven by turbulent gluon cascades, is governed by a non-linear rate equation that resums gluon splittings at arbitrary angles and is enhanced by the medium length, L. The solution of this equation sets the initial condition for a non-linear DGLAP-like evolution equation, which describes the collinear early vacuum cascade resolved by the medium at angles exceeding the medium resolution angle, θ c . For asymptotic jet energies, the medium-induced cascade displays an exponential behavior that generalizes the Poisson-like distribution of parton energy loss. This formulation enables the resummation of leading contributions in α s ln(1/R), and α s ln(R/θ c ), and powers of α s L. We briefly explore the limit of strong quenching, where analytic treatments are feasible, offering insights into the impact of parton cascades on jet quenching. These results provide guidance for future numerical simulations and analytical investigations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2018

Linear fitting of multi-threshold counting data

The present disclosure provides a system and method for efficiently mining multi-threshold measurements acquired using photon counting pixel-array detectors for spectral imaging and diffraction analyses. Images of X-ray intensity as a function of X-ray energy were recorded on a 6 megapixel X-ray photon counting array detector through linear fitting of the measured counts recorded as a function of counting threshold. An analytical model is disclosed for describing the probability density of detected voltage, utilizing fractional photon counting to account for edge/corner effects from voltage plumes that spread across multiple pixels. Three-parameter fits to the model were independently performed for each pixel in the array for X-ray scattering images acquired for 13.5 keV and 15.0 keV X-ray energies. From the established pixel responses, multi-threshold composite images produced from the sum of 13.5 keV and 15.0 keV data can be analytically separated to recover the monochromatic images through simple linear fitting.

47 OTHER INSTRUMENTATION↗

DOE OSTI2012

GRAVITATIONAL INSTABILITY OF ROTATING, PRESSURE-CONFINED, POLYTROPIC GAS DISKS WITH VERTICAL STRATIFICATION

We investigate the gravitational instability (GI) of rotating, vertically stratified, pressure-confined, polytropic gas disks using a linear stability analysis as well as analytic approximations. The disks are initially in vertical hydrostatic equilibrium and bounded by a constant external pressure. We find that the GI of a pressure-confined disk is in general a mixed mode of the conventional Jeans and distortional instabilities, and is thus an unstable version of acoustic-surface-gravity waves. The Jeans mode dominates in weakly confined disks or disks with rigid boundaries. On the other hand, when the disk has free boundaries and is strongly pressure confined, the mixed GI is dominated by the distortional mode that is surface-gravity waves driven unstable under their own gravity and thus incompressible. We demonstrate that the Jeans mode is gravity-modified acoustic waves rather than inertial waves and that inertial waves are almost unaffected by self-gravity. We derive an analytic expression for the effective sound speed c{sub eff} of acoustic-surface-gravity waves. We also find expressions for the gravity reduction factors relative to a razor-thin counterpart that are appropriate for the Jeans and distortional modes. The usual razor-thin dispersion relation, after correcting for c{sub eff} and the reduction factors, closely matches the numerical results obtained by solving a full set of linearized equations. The effective sound speed generalizes the Toomre stability parameter of the Jeans mode to allow for the mixed GI of vertically stratified, pressure-confined disks.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

Application of Adomian Decomposition for Multi-Machine Power System Simulation

In this paper, the merits of Adomian Decomposition Method (ADM) have been investigated for the time domain simulation of multi-machine power systems. ADM is an analytical approximation method for the solution of linear, nonlinear, deterministic and stochastic operator equations. It approximates the nonlinearities in the equations using Adomian polynomial series and that can form a rapidly convergent sequence of analytic functions. We apply the ADM approach for simulation of differential algebraic equations representing the detailed multimachine power system. The ADM approach has been applied on the widely used IEEE 3 generator 9 bus system and IEEE 10 generator 39 bus system for 3 fault simulations. In our simulations we found that the ADM approach is faster than the standard trapezoidal time integration method for the comparable accuracy.

Gurrala, Gurunath↗

DOE OSTI2020

Detecting Magnetic Fields in Exoplanets with Spectropolarimetry of the Helium Line at 1083 nm

The magnetic fields of the solar system’s planets provide valuable insights into their interiors and can have dramatic consequences for the evolution of their atmospheres and interaction with the solar wind. However, we have little direct knowledge of magnetic fields in exoplanets. Here we present a method for detecting magnetic fields in the atmospheres of close-in exoplanets based on spectropolarimetric transit observations at the wavelength of the helium line at 1083 nm. This methodology has been successfully applied for exploring magnetic fields in solar coronal filaments. Strong absorption signatures (transit depths on the order of a few percent) in the 1083 nm line have recently been observed for several close-in exoplanets. We show that in the conditions in these escaping atmospheres, metastable helium atoms should be optically pumped by the starlight and, for field strengths more than a few × 10 -4 G, should align with the magnetic field. This results in linearly polarized absorption at 1083 nm that traces the field direction (the Hanle effect), which we explore by both analytic computation and the HAZEL numerical code. The linear polarization $\sqrt{{Q}^{2}+{U}^{2}}/I$ ranges from ~10 -3 in optimistic cases down to a few × 10 -5 for particularly unfavorable cases, with very weak dependence on field strength. The line-of-sight component of the field results in a slight circular polarization (the Zeeman effect), also reaching $V/I\sim \mathrm{few}\times {10}^{-5}({B}_{\parallel }/10\,{\rm{G}})$. We discuss the detectability of these signals with current (SPIRou) and future (Extremely Large Telescope) high-resolution infrared spectropolarimeters, and we briefly comment on possible sources of astrophysical contamination.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Analytic solutions of the DGLAP evolution and theoretical uncertainties

The energy dependence for the singlet sector of Parton Distributions Functions (PDFs) is described by an entangled pair of ordinary linear differential equations. Although there are no exact analytic solutions, it is possible to provide approximated results depending on the assumptions and the methodology adopted. These results differ in their sub-leading, neglected terms and ultimately they are associated with different treatments of the theoretical uncertainties. In this work, a novel analytic approach in Mellin space is presented and a new methodology for obtaining closed and exponentiated analytic solutions is devised. Different results for the DGLAP evolution at Next-Leading-Order are compared, discussing advantages and disadvantages for each solution. The generalizations to higher orders are addressed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2013

Exact PDF equations and closure approximations for advective-reactive transport

Mathematical models of advection–reaction phenomena rely on advective flow velocity and (bio) chemical reaction rates that are notoriously random. By using functional integral methods, we derive exact evolution equations for the probability density function (PDF) of the state variables of the advection–reaction system in the presence of random transport velocity and random reaction rates with rather arbitrary distributions. These PDF equations are solved analytically for transport with deterministic flow velocity and a linear reaction rate represented mathematically by a heterog eneous and strongly-correlated random field. Our analytical solution is then used to investigate the accuracy and robustness of the recently proposed large-eddy diffusivity (LED) closure approximation [1]. We find that the solution to the LED-based PDF equation, which is exact for uncorrelated reaction rates, is accurate even in the presence of strong correlations and it provides an upper bound of predictive uncertainty.

PDF methods, Uncertainty quantification, Heterogen↗

DOE OSTI2011

Non-linear matter power spectrum from Time Renormalisation Group: efficient computation and comparison with one-loop

We address the issue of computing the non-linear matter power spectrum on mildly non-linear scales with efficient semi-analytic methods. We implemented M. Pietroni's Time Renormalization Group (TRG) method and its Dynamical 1-Loop (D1L) limit in a numerical module for the new Boltzmann code CLASS. Our publicly released module is valid for ΛCDM models, and optimized in such a way to run in less than a minute for D1L, or in one hour (divided by number of nodes) for TRG. A careful comparison of the D1L, TRG and Standard 1-Loop approaches reveals that results depend crucially on the assumed initial bispectrum at high redshift. When starting from a common assumption, the three methods give roughly the same results, showing that the partial resumation of diagrams beyond one loop in the TRG method improves one-loop results by a negligible amount. A comparison with highly accurate simulations by M. Sato and T. Matsubara shows that all three methods tend to over-predict non-linear corrections by the same amount on small wavelengths. Percent precision is achieved until k ∼ 0.2 hMpc{sup −1} for z ≥ 2, or until k ∼ 0.14 hMpc{sup −1} at z = 1.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Precision beam diagnostics at the NuMI facility using muon monitor observations

The Neutrinos at the Main Injector (NuMI) facility at Fermilab delivers an intense neutrino beam for multiple experiments by producing pions that decay into neutrinos, muons, and other particles. Magnetic horns—the primary pion focusing elements in the NuMI beamline—exhibit predominantly linear optics, enabling a predictable relationship between the proton beam and the resulting pion and muon phase spaces. This study has two primary objectives: first, to evaluate and confirm the linearity of the horn focusing mechanism using analytical models and numerical simulations; and second, to demonstrate that key beam parameters—such as proton beam intensity, beam position on target, and horn current—can be extracted from muon monitor observations within this linear optics framework. Using a machine learning model trained on spill-by-spill muon monitor data, we infer the horn current with a precision of ±0.05%, the beam intensity with ±0.1%, and the beam position on target with ±0.018⁢ mm horizontally and ±0.013⁢ mm vertically. This approach provides a reliable cross-check of beam parameters, helping to reduce systematic uncertainties that are critical for future experiments such as the Deep Underground Neutrino Experiment, which will rely on the neutrino beam produced by the Long-Baseline Neutrino Facility.

Beam control↗

DOE OSTI2011

Multisource modeling of flattening filter free (FFF) beam and the optimization of model parameters

Purpose: With the introduction of flattening filter free (FFF) linear accelerators to radiation oncology, new analytical source models for a FFF beam applicable to current treatment planning systems is needed. In this work, a multisource model for the FFF beam and the optimization of involved model parameters were designed. Methods: The model is based on a previous three source model proposed by Yang et al. [''A three-source model for the calculation of head scatter factors,'' Med. Phys. 29, 2024-2033 (2002)]. An off axis ratio (OAR) of photon fluence was introduced to the primary source term to generate cone shaped profiles. The parameters of the source model were determined from measured head scatter factors using a line search optimization technique. The OAR of the photon fluence was determined from a measured dose profile of a 40x40 cm{sup 2} field size with the same optimization technique, but a new method to acquire gradient terms for OARs was developed to enhance the speed of the optimization process. The improved model was validated with measured dose profiles from 3x3 to 40x40 cm{sup 2} field sizes at 6 and 10 MV from a TrueBeam STx linear accelerator. Furthermore, planar dose distributions for clinically used radiation fields were also calculated and compared to measurements using a 2D array detector using the gamma index method. Results: All dose values for the calculated profiles agreed with the measured dose profiles within 0.5% at 6 and 10 MV beams, except for some low dose regions for larger field sizes. A slight overestimation was seen in the lower penumbra region near the field edge for the large field sizes by 1%-4%. The planar dose calculations showed comparable passing rates (>98%) when the criterion of the gamma index method was selected to be 3%/3 mm. Conclusions: The developed source model showed good agreements between measured and calculated dose distributions. The model is easily applicable to any other linear accelerator using FFF beams as the required data include only the measured PDD, dose profiles, and output factors for various field sizes, which are easily acquired during conventional beam commissioning process.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2013

Coronal electron confinement by double layers

In observations of flare-heated electrons in the solar corona, a longstanding problem is the unexplained prolonged lifetime of the electrons compared to their transit time across the source. This suggests confinement. Recent particle-in-cell (PIC) simulations, which explored the transport of pre-accelerated hot electrons through ambient cold plasma, showed that the formation of a highly localized electrostatic potential drop, in the form of a double layer (DL), significantly inhibited the transport of hot electrons. The effectiveness of confinement by a DL is linked to the strength of the DL as defined by its potential drop. In this work, we investigate the scaling of the DL strength with the hot electron temperature by PIC simulations and find a linear scaling. We demonstrate that the strength is limited by the formation of parallel shocks. Based on this, we analytically determine the maximum DL strength, and also find a linear scaling with the hot electron temperature. The DL strength obtained from the analytic calculation is comparable to that from the simulations. At the maximum strength, the DL is capable of confining a significant fraction of hot electrons in the source.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2016

Effect of PyC Irradiation-Induced Creep on the Stress of TRISO-Coated Fuel Particles - 18474

Irradiation-induced Pyrocarbon (PyC) creep affects the stresses of a TRISO-coated particle. An analytical solution with a 4-parameter linear viscoelastic model was used to investigate the effects of PyC creep on the particle stresses. PyC transient and steady-state creep effects on particle stresses were analyzed. PyC creep parameter effects on the particle stresses were studied, including PyC creep coefficient, creep Poisson s ratio in the steady-state component, and Young s modulus in the transient component. The results showed that PyC irradiation-induced creep decreased the particle stresses. Transient creep affected particle stresses only at the irradiation early stage. Lower Young s modulus and higher creep coefficient in the transient component resulted in a larger stresses decrease, and large Poisson s ratio and creep coefficient in the steady-state creep decreased particle stresses. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2024

Terahertz-polarizing effect based on geometric anisotropy of Ti 3 C 2 T x MXene nanosheets

Two-dimensional (2D) nanomaterials have garnered considerable attention due to their unique properties, such as thinness and excellent electronic properties. The inherent anisotropy of 2D nanomaterials, which is characterized by their high aspect ratio, can maximize the directional functionality to enhance the conductivity and utilization of edges. Here, we demonstrate how to utilize the edges of titanium carbide (Ti 3 C 2 T x ) MXene nanosheet to control the linear polarization in the terahertz (THz) frequency range. By a simple and unique vertical alignment exploiting the high surface charge and excellent colloidal dispersibility of MXene nanosheets, MXene nanosheets are vertically oriented along an applied AC electric field. The well-aligned edges justify the THz polarizing effect, revealing outstanding performances such as the broadband coverage in 0.2–1.5 THz and an exceptional extinction ratio reaching 20 dB. Furthermore, we vary the local orientation of the vertical MXene nanosheets by tailoring the applied electric field, enabling spatially controllable polarization capabilities. Our results can provide a generic tool for the utilization of 2D nanomaterials in diverse potential applications in optics and electronics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Spatiotemporal spread of perturbations in power-law models at low temperatures: Exact results for classical out-of-time-order correlators

Here, we present exact results for the classical version of the out-of-time-order commutator (OTOC) for a family of power-law models consisting of N particles in one dimension and confined by an external harmonic potential. These particles are interacting via power-law interaction of the form ∝ $Σ_{i,j=(i≠j)}^{N}$ |x i -x j | -k for every k>1 where x i is the position of the i th particle. We present numerical results for the OTOC for finite N at low temperatures and short enough times so that the system is well approximated by the linearized dynamics around the many-body ground state. In the large-N limit, we compute the ground-state dispersion relation in the absence of external harmonic potential exactly and use it to arrive at analytical results for OTOC. We find excellent agreement between our analytical results and the numerics. We further obtain analytical results in the limit where only linear and leading nonlinear (in momentum) terms in the dispersion relation are included. The resulting OTOC is in agreement with numerics in the vicinity of the edge of the “light cone.” We find remarkably distinct features in OTOC below and above k=3 in terms of going from non-Airy behavior (1 < k < 3) to an Airy universality class (k>3). We present certain additional rich features for the case k=2 that stem from the underlying integrability of the Calogero-Moser model. We present a field theory approach that also assists in understanding certain aspects of OTOC such as the sound speed. Our findings are a step forward towards a more general understanding of the spatiotemporal spread of perturbations in long-range interacting systems.

1-dimensional systems↗

DOE OSTI2014

Joint discrete universality of Hurwitz zeta functions

We obtain a joint discrete universality theorem for Hurwitz zeta functions. Here the parameters of zeta functions and the step of shifts of these functions approximating a given family of analytic functions are connected by some condition of linear independence. Nesterenko's theorem gives an example satisfying this condition. The universality theorem is applied to estimate the number of zeros of a linear combination of Hurwitz zeta functions. Bibliography: 20 titles.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Optimizing Piezoelectric Nanocomposites by High‐Throughput Phase‐Field Simulation and Machine Learning

Abstract Piezoelectric nanocomposites with oxide fillers in a polymer matrix combine the merit of high piezoelectric response of the oxides and flexibility as well as biocompatibility of the polymers. Understanding the role of the choice of materials and the filler‐matrix architecture is critical to achieving desired functionality of a composite towards applications in flexible electronics and energy harvest devices. Herein, a high‐throughput phase‐field simulation is conducted to systematically reveal the influence of morphology and spatial orientation of an oxide filler on the piezoelectric, mechanical, and dielectric properties of the piezoelectric nanocomposites. It is discovered that with a constant filler volume fraction, a composite composed of vertical pillars exhibits superior piezoelectric response and electromechanical coupling coefficient as compared to the other geometric configurations. An analytical regression is established from a linear regression‐based machine learning model, which can be employed to predict the performance of nanocomposites filled with oxides with a given set of piezoelectric coefficient, dielectric permittivity, and stiffness. This work not only sheds light on the fundamental mechanism of piezoelectric nanocomposites, but also offers a promising material design strategy for developing high‐performance polymer/inorganic oxide composite‐based wearable electronics.

42 ENGINEERING↗