Search NASA⌕ Search

SEARCH · Search NASA

Results for “TAYLOR THEOREM”

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.

Task-based detectability in CT image reconstruction by filtered backprojection and penalized likelihood estimation

Purpose: Nonstationarity is an important aspect of imaging performance in CT and cone-beam CT (CBCT), especially for systems employing iterative reconstruction. This work presents a theoretical framework for both filtered-backprojection (FBP) and penalized-likelihood (PL) reconstruction that includes explicit descriptions of nonstationary noise, spatial resolution, and task-based detectability index. Potential utility of the model was demonstrated in the optimal selection of regularization parameters in PL reconstruction. Methods: Analytical models for local modulation transfer function (MTF) and noise-power spectrum (NPS) were investigated for both FBP and PL reconstruction, including explicit dependence on the object and spatial location. For FBP, a cascaded systems analysis framework was adapted to account for nonstationarity by separately calculating fluence and system gains for each ray passing through any given voxel. For PL, the point-spread function and covariance were derived using the implicit function theorem and first-order Taylor expansion according toFessler [“Mean and variance of implicitly defined biased estimators (such as penalized maximum likelihood): Applications to tomography,” IEEE Trans. Image Process. 5(3), 493–506 (1996)]. Detectability index was calculated for a variety of simple tasks. The model for PL was used in selecting the regularization strength parameter to optimize task-based performance, with both a constant and a spatially varying regularization map. Results: Theoretical models of FBP and PL were validated in 2D simulated fan-beam data and found to yield accurate predictions of local MTF and NPS as a function of the object and the spatial location. The NPS for both FBP and PL exhibit similar anisotropic nature depending on the pathlength (and therefore, the object and spatial location within the object) traversed by each ray, with the PL NPS experiencing greater smoothing along directions with higher noise. The MTF of FBP is isotropic and independent of location to a first order approximation, whereas the MTF of PL is anisotropic in a manner complementary to the NPS. Task-based detectability demonstrates dependence on the task, object, spatial location, and smoothing parameters. A spatially varying regularization “map” designed from locally optimal regularization can improve overall detectability beyond that achievable with the commonly used constant regularization parameter. Conclusions: Analytical models for task-based FBP and PL reconstruction are predictive of nonstationary noise and resolution characteristics, providing a valuable framework for understanding and optimizing system performance in CT and CBCT.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Developments in TEM Nanotomography of Calcium Silicate Hydrate

This investigation was designed to explore the possibility of using transmission electron microscope ( TEM ) tomography on cement‐based systems gain a greater understanding of their nanostructure and pore network. The preliminary results show a clearly a well‐defined pore network at the nanoscale, with pore size approximately 1.7–2.4 nm in diameter and spaced around 5–8 nm apart. A comparison of small angle X‐ray scattering data with 2‐D TEM images analyzed with the Fourier slice theorem documents an excellent structural correlation.

Taylor, Rae↗

Schrieffer-Wolff transformation for quantum many-body systems

The Schrieffer-Wolff (SW) method is a version of degenerate perturbation theory in which the low-energy effective Hamiltonian H{sub eff} is obtained from the exact Hamiltonian by a unitary transformation decoupling the low-energy and high-energy subspaces. We give a self-contained summary of the SW method with a focus on rigorous results. We begin with an exact definition of the SW transformation in terms of the so-called direct rotation between linear subspaces. From this we obtain elementary proofs of several important properties of H{sub eff} such as the linked cluster theorem. We then study the perturbative version of the SW transformation obtained from a Taylor series representation of the direct rotation. Our perturbative approach provides a systematic diagram technique for computing high-order corrections to H{sub eff}. We then specialize the SW method to quantum spin lattices with short-range interactions. We establish unitary equivalence between effective low-energy Hamiltonians obtained using two different versions of the SW method studied in the literature. Finally, we derive an upper bound on the precision up to which the ground state energy of the nth-order effective Hamiltonian approximates the exact ground state energy. - Highlights: > The Schrieffer-Wolff transformation is specialized to quantum spin lattices with short-range interactions. > We provide a diagram technique for computing high-order corrections to the effective low-energy Hamiltonian. > We derive a rigorous bound on the error up to which the nth-order effective low-energy dynamics approximates the exact dynamics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A high order Cartesian grid, finite volume method for elliptic interface problems

We present a higher-order finite volume method for solving elliptic PDEs with jump conditions on interfaces embedded in a 2D Cartesian grid. Second, fourth, and sixth order accuracy is demonstrated on a variety of tests including problems with high-contrast and spatially varying coefficients, large discontinuities in the source term, and complex interface geometries. We include a generalized truncation error analysis based on cell-centered Taylor series expansions, which then define stencils in terms of local discrete solution data and geometric information. In the process, we develop a simple method based on Green's theorem for computing exact geometric moments directly from an implicit function definition of the embedded interface. This approach produces stencils with a simple bilinear representation, where spatially-varying coefficients and jump conditions can be easily included and finite volume conservation can be enforced.

97 MATHEMATICS AND COMPUTING↗

Diffusive magnetohydrodynamic instabilities beyond the Chandrasekhar theorem

We consider the stability of axially unbounded cylindrical flows that contain a toroidal magnetic background field with the same radial profile as their azimuthal velocity. For ideal fluids, Chandrasekhar had shown the stability of this configuration if the Alfvén velocity of the field equals the velocity of the background flow, i.e., if the magnetic Mach number Mm=1. We demonstrate that magnetized Taylor–Couette flows with such profiles become unstable against non-axisymmetric perturbations if at least one of the diffusivities is finite. We also find that for small magnetic Prandtl numbers Pm the lines of marginal instability scale with the Reynolds number and the Hartmann number. In the limit Pm→0 the lines of marginal instability completely lie below the line for Mm=1 and for Pm→∞ they completely lie above this line. For any finite value of Pm, however, the lines of marginal instability cross the line Mm=1, which separates slow from fast rotation. The minimum values of the field strength and the rotation rate that are needed for the instability (slightly) grow if the rotation law becomes flat. In this case, the electric current of the background field becomes so strong that the current-driven Tayler instability (which also exists without rotation) appears in the bifurcation map at low Hartmann numbers.

79 ASTRONOMY AND ASTROPHYSICS↗

A new data processing technique for Rayleigh-Taylor instability growth experiments

Typical face-on experiments for Rayleigh-Taylor instability study involve the time-resolved radiography of an accelerated foil with line-of-sight of the radiography along the direction of motion. The usual method which derives perturbation amplitudes from the face-on images reverses the actual image transmission procedure, so the obtained results will have a large error in the case of large optical depth. In order to improve the accuracy of data processing, a new data processing technique has been developed to process the face-on images. This technique based on convolution theorem, refined solutions of optical depth can be achieved by solving equations. Furthermore, we discuss both techniques for image processing, including the influence of modulation transfer function of imaging system and the backlighter spatial profile. Besides, we use the two methods to the process the experimental results in Shenguang-II laser facility and the comparison shows that the new method effectively improve the accuracy of data processing.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Analytical examples of reversal current, zero core current, and surface current, toroidal magnetostatic equilibria with nested flux surfaces

We present exact analytical examples of three types of axisymmetric toroidal magnetostatic equilibria with nested flux surfaces: (1) current reversal equilibria, for which the net toroidal current switches from a negative to a positive value when moving away from the magnetic axis; these equilibria have a non-monotonic pressure profile, in accordance with Hammett et al.'s theorem stating that the pressure on the current reversal surface has to exceed the volume-averaged pressure within that surface; (2) zero core current equilibria, in which the toroidal current density vanishes inside some flux surface; and (3) surface current equilibria, constituted of an arbitrary number of nested layers inside which the plasma pressure is constant and the magnetic field force-free, with two adjacent layers being separated by a current sheet. All these configurations are obtained by shaping in an adequate way the arbitrary function which intervenes in the class of generalized isodynamic equilibria first constructed by Palumbo and recovered later on by Bishop and Taylor. A derivation of these equilibria by a method slightly different from Palumbo's is given in an Appendix.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Hall interchange instability as a seed for helical magneto-Rayleigh–Taylor instabilities in magnetized liner inertial fusion Z-Pinches scaled from Z-Machine parameters to a next generation pulsed power facility

Magnetized liner inertial fusion (MagLIF) is a magneto-inertial-fusion concept that is studied on the 20-MA, 100-ns rise time Z Pulsed Power Facility at Sandia National Laboratories. Given the relative success of the platform, there is a wide interest in studying the scaled performance of this concept at a next-generation pulsed-power facility that may produce peak currents upward of 60 MA. An important aspect that requires more research is the instability dynamics of the imploding MagLIF liner, specifically how instabilities are initially seeded. It has been shown in magnetized 1-MA thin-foil liner Z-pinch implosion simulations that a Hall interchange instability (HII) effect can provide an independent seeding mechanism for helical magneto-Rayleigh–Taylor instabilities. Here in this paper, we explore this instability at higher peak currents for MagLIF using 2D discontinuous Galerkin PERSEUS simulations, an extended magneto-hydrodynamics code, which includes Hall physics. Our simulations of scaled MagLIF loads show that the growth rate of the HII is invariant to the peak current, suggesting that studies at 20-MA are directly relevant to 60-MA class machines.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Towards Provable Security in Industrial Control Systems Via Dynamic Protocol Attestation

Industrial control systems (ICSs) increasingly rely on digital technologies vulnerable to cyber attacks. Cyber attackers can infiltrate ICSs and execute malicious actions. Individually, each action seems innocuous. But taken together, they cause the system to enter an unsafe state. These attacks have resulted in dramatic consequences such as physical damage, economic loss, and environmental catastrophes. This paper introduces a methodology that restricts actions using protocols. These protocols only allow safe actions to execute. Protocols are written in a domain specific language we have embedded in an interactive theorem prover (ITP). The ITP enables formal, machine-checked proofs to ensure protocols maintain safety properties. We use dynamic attestation to ensure ICSs conform to their protocol even if an adversary compromises a component. Since protocol conformance prevents unsafe actions, the previously mentioned cyber attacks become impossible. We demonstrate the effectiveness of our methodology using an example from the Fischertechnik Industry 4.0 platform. We measure dynamic attestation's impact on latency and throughput. Our approach is a starting point for studying how to combine formal methods and protocol design to thwart attacks intended to cripple ICSs.

97 MATHEMATICS AND COMPUTING↗

Magnetized liner inertial fusion platform development to assess performance scaling with drive parameters

Magnetized liner inertial fusion (MagLIF) experiments have demonstrated fusion-relevant ion temperatures up to 3.1 keV and thermonuclear production of up to 1.1 × 1013 deuterium–deuterium neutrons. This performance was enabled through platform development that provided increases in applied magnetic field, coupled preheat energy, and drive current. Advanced coil designs with internal reinforcement enabled an increase from 10 to 20 T. An improved laser pulse shape, beam smoothing, and thinner laser entrance foils increased preheat energy coupling from less than 1 to 2.3 kJ. A redesign of the final transmission line and load region increased peak load current from 16 to 20 MA. The wider range of input parameters was leveraged to study target performance trends with preheat energy, applied magnetic field, and peak load current. Ion temperature and neutron yield generally followed trends in two-dimensional clean Lasnex calculations. Stagnation performance improved with peak load current when other input parameters were also increased such that convergence was maintained. This dataset suggests that reducing convergence to less than 30 would improve predictability of target performance. Lasnex was used to identify a simulation-optimized scaling path, which suggests 10+ kJ of fusion yield is possible on the Z facility with achievable input parameters. This path also indicates >10 MJ could be generated through volume burn on a future facility with a path to high yield (>200 MJ) using cryogenic dense fuel layers. The newly developed MagLIF platform enables exploration of both this simulation optimized scaling path and a recently developed similarity-scaling path.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

STOCHASTIC ACCELERATION AND THE EVOLUTION OF SPECTRAL DISTRIBUTIONS IN SYNCHRO-SELF-COMPTON SOURCES: A SELF-CONSISTENT MODELING OF BLAZARS' FLARES

The broadband spectral distributions of non-thermal sources, such as those of several known blazars, are well described by a log-parabolic fit. The second-degree term in these fits measures the curvature in the spectrum. In this paper, we investigate whether the curvature parameter observed in the spectra of the synchrotron emission can be used as a fingerprint of stochastic acceleration. As a first approach, we use the multiplicative central limit theorem to show how fluctuations in the energy gain result in the broadening of the spectral shape, introducing a curvature into the energy distribution. Then, by means of a Monte Carlo description, we investigate how the curvature produced in the electron distribution is linked to the diffusion in momentum space. To get a more generic description of the problem we turn to the diffusion equation in momentum space. We first study some 'standard' scenarios, in order to understand the conditions that make the curvature in the spectra significant, and the relevance of cooling during the acceleration process. We try to quantify the correlation between the curvature and the diffusive process in the pre-equilibrium stage, and investigate how the transition between the Klein-Nishina and the Thomson regimes, in inverse Compton cooling, determine the curvature in the distribution at equilibrium. We apply these results to some observed trends, such as the anticorrelation between the peak energy and the curvature term observed in the spectra of Mrk 421, and a sample of BL Lac objects whose synchrotron emission peaks at X-ray energies.

79 ASTRONOMY AND ASTROPHYSICS↗

Hall instability driven seeding of helical magneto-Rayleigh–Taylor instabilities in axially premagnetized thin-foil liner Z-pinch implosions

Helical magneto-Rayleigh–Taylor instability (MRTI) structures have been observed in z-pinch-driven liner implosion experiments with a pre-imposed axial magnetic field. We show that the formation of these helical structures can be described by a Hall magnetohydrodynamical (HMHD) model. We used the 3D extended magnetohydrodynamics simulation code PERSEUS (which includes Hall physics) [Seyler and Martin, Phys. Plasmas 18, 012703 (2011)] to study these helical instabilities and show that a Hall interchange instability in low-density coronal plasma immediately surrounding the dense liner is responsible for producing helically oriented effects in the magnetic field and current density within the coronal layer. Here, this seeds the helical pitch angle of the MRTI even when other proposed helical seeding mechanisms are either not present in the experiments or not accounted for in the simulations. For example, this mechanism does not require low-density power-feed plasmas to be swept in from large radius or the development of electrothermal instabilities. The Hall Instability is, thus, a new, independent explanation for the origin of the helical instabilities observed in axially premagnetized liner experiments. Simulation results supporting this mechanism are presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Resilience Measurement Framework For Post-deployment Artificial Intelligence (ai) Integrated Systems

Resilience is largely defined as the ability to adapt or recover from adverse conditions, stresses, attacks, or compromises on systems that use or are enabled by digital resources. In Artificial Intelligence Management and Research for Advanced Networked Testbed Hub (AMARANTH), resilience is measured in the amount of time it took from the beginning of a testing period for the model to reach predictions outside of the original 95% confidence interval or using the Kullback-Leibler (KL) divergence theorem, the Population Stability Index (PSI), and traditional methods such as root mean squared error (RMSE) threshold. Artificial Intelligence (AI) model drift is of significant concern when deploying AI-integrated systems into critical and/or secure environments. Drift can impact resilience of the AI-integrated system post-deployment and requires consistent maintenance and upkeep to ensure the model is accurate and precise. To quantify model drift and predict the point when a model's drift becomes unacceptable, we describe using Kullback-Leibler (KL) divergence, Population Stability Index (PSI) and/or confidence interval width estimations to determine the point of failure and time to failure of a model post-deployment. Through simple code functions, the KL-divergence, PSI, confidence interval, and root mean squared (RMSE) point of failures can be used to derive when a model needs to be maintained as well as the impact of adversarial action through statistical means.

Yockey, Patience [Idaho National Laboratory (INL),↗

First measurement of elastic, inelastic and total cross-section at $$\sqrt{s}=13$$ s = 13 TeV by TOTEM and overview of cross-section data at LHC energies: TOTEM Collaboration

The TOTEM collaboration has measured the proton–proton total cross section at s√=13 TeV with a luminosity-independent method. Using dedicated β*=90 m beam optics, the Roman Pots were inserted very close to the beam. The inelastic scattering rate has been measured by the T1 and T2 telescopes during the same LHC fill. After applying the optical theorem the total proton–proton cross section is σ tot =(110.6 ± 3.4 ) mb, well in agreement with the extrapolation from lower energies. This method also allows one to derive the luminosity-independent elastic and inelastic cross sections σ el =(31.0 ± 1.7) mb and σ inel =(79.5 ± 1.8) mb .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗