Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numerical Oscillation”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Finite-temperature vibrational full configuration interaction

Thermodynamic functions of an ideal molecular gas due to its anharmonic vibrations are evaluated in a wide range of temperature (T) by the vibrational full-configuration-interaction (FCI) method using a quartic force field and a finite number (N) of harmonic-oscillator basis functions along each normal mode. The thermodynamic functions considered are the grand potential (Ω), internal energy (U), and entropy (S). They are compared with those obtained from the Bose–Einstein theory with or without truncation of the harmonic-oscillator basis functions after quantum number N–1. The comparison reveals that the finite-basis-set errors in Ω and U are, respectively, k B Tln(k B T/Nℏω) and k B T per mode in the high-T limit, obscuring anharmonic effects when k B T > ℏω, where ω is the lowest mode frequency. Here, the benchmark data for several low-order perturbation corrections to Ω, U, and S are also obtained as the numerical derivatives of their FCI values with respect to dimensionless perturbation strength, and the domain of T and N in which these data are reliable (for the N → ∞ limits) is discussed.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

A simple numerical method for evaluating heat dissipation from curved wires with periodic applied heating

In many situations, the dual-purpose heater/thermometer wires used in the three-omega method – one of the most precise and sensitive techniques for measuring the thermal conductivity of thin films and interfaces – must include bends and curves to avoid obstructions on the surface of a sample. Although the three-omega analysis assumes that the heating wire is infinitely long and straight, recent experimental work has demonstrated that in some cases curved-wire geometries can be used without introducing detectable systematic error. We describe a general numerical method that can be used to calculate the temperature of three-omega heating wires with arbitrary wire geometries. Furthermore, this method provides experimentalists with a simple quantitative procedure for calculating how large the systematic error caused by a particular wire asymmetry will be. We show calculations of two useful cases: a straight wire with a single bend of arbitrary angle and a wire that forms a circle. We find that the amplitude of the in-phase temperature oscillations near a wire that forms a circle differs from the prediction using the analytic straight-line source solution by <12%, provided that the thermal penetration depth is less than ten times the radius of curvature of the wire path. The in-phase temperature amplitude 1.5 wire widths away from a 90° bend in a wire is within 11% of the straight-line source prediction for all penetration depths greater than the wire width. Our calculations indicate that the straight-line source solution breaks down significantly when the wire bend angle is less than 45°.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2022

Computational study of runaway electrons in MST tokamak discharges with applied resonant magnetic perturbation

A numerical study of magnetohydrodynamics (MHD) and tracer-particle evolution investigates the effects of resonant magnetic perturbations (RMPs) on the confinement of runaway electrons (REs) in tokamak discharges conducted in the Madison Symmetric Torus. In computational results of applying RMPs having a broad toroidal spectrum but a single poloidal harmonic, m = 1 RMP does not suppress REs, whereas m = 3 RMP achieves significant deconfinement but not the complete suppression obtained in the experiment. MHD simulations with the NIMROD code produce sawtooth oscillations, and the associated magnetic reconnection can affect the trajectory of REs starting in the core region. Simulations with m = 3 RMP produce chaotic magnetic topology over the outer region, but the m = 1 RMP produces negligible changes in field topology, relative to applying no RMP. Using snapshots of the MHD simulation fields, full-orbit relativistic electron test particle computations with KORC show ≈50% loss from the m = 3 RMP compared to the 10%–15% loss from the m = 1 RMP. Here, test particle computations of the m = 3 RMP in the time-evolving MHD simulation fields show correlation between MHD activity and late-time particle losses, but total electron confinement is similar to computations using magnetic-field snapshots.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Energetic Nonthermal Electrons within the Above-the-looptop Regions in Solar Flares: Acceleration, Feedback, and Quasiperiodic Pulsations

Solar flares are among the most dramatic events in the solar system, releasing substantial magnetic energy and accelerating a large number of electrons to high energies. Notably, in certain events, the above-the-looptop region may contain a significant population of nonthermal electrons, both in number and energy. For the first time, we adopt a novel numerical method that combines magnetohydrodynamics with energetic particles incorporating feedback from nonthermal electrons to investigate electron acceleration and transport in solar flares. We find that a large fraction of energetic electrons are accelerated via the current sheet and termination shock regions. Most energetic electrons are concentrated in the above-the-looptop region, carrying a sizable amount of the released energy. We observe that greater feedback of nonthermal electrons leads to steeper energy spectra. The energy density of the nonthermal electrons oscillates due to the periodic impact of magnetic islands into the above-the-looptop region, which may help explain the observed quasiperiodic pulsations. Our simulations provide new insights into the origin of nonthermal electrons and associated emissions in the above-the-looptop region.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Interpretable Joint Event-Particle Reconstruction for Neutrino Physics at NOvA with Sparse CNNs and Transformers

The complex events observed at the NOvA long-baseline neutrino oscillation experiment contain vital information for understanding the most elusive particles in the standard model. The NOvA detectors observe interactions of neutrinos from the NuMI beam at Fermilab. Associating the particles produced in these interaction events to their source particles, a process known as reconstruction, is critical for accurately measuring key parameters of the standard model. Events may contain several particles, each producing sparse high-dimensional spatial observations, and current methods are limited to evaluating individual particles. To accurately label these numerous, high-dimensional observations, we present a novel neural network architecture that combines the spatial learning enabled by convolutions with the contextual learning enabled by attention. This joint approach, TransformerCVN, simultaneously classifies each event and reconstructs every individual particle's identity. TransformerCVN classifies events with 90% accuracy and improves the reconstruction of individual particles by 6% over baseline methods which lack the integrated architecture of TransformerCVN. In addition, this architecture enables us to perform several interpretability studies which provide insights into the network's predictions and show that TransformerCVN discovers several fundamental principles that stem from the standard model.

Shmakov, Alexander↗
DOE OSTI2022

Spectral function for 4 He using the Chebyshev expansion in coupled-cluster theory

Here, we compute spectral function for 4 He by combining coupled-cluster theory with an expansion of integral transforms into Chebyshev polynomials. Our method allows us to estimate the uncertainty of spectral reconstruction. The properties of the Chebyshev polynomials make the procedure numerically stable and considerably lower in memory usage than the typically employed Lanczos algorithm. We benchmark our predictions with other calculations in the literature and with electron-scattering data in the quasi-elastic peak. The spectral function formalism allows one to extend ab initio lepton-nucleus cross sections into the relativistic regime. This makes it a promising tool for modeling this process at higher-energy transfers. The results we present open the door for studies of heavier nuclei, important for the neutrino oscillation programs.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Iterative quantum optimization of spin glass problems with rapidly oscillating transverse fields

In this work, we introduce a new iterative quantum algorithm, called Iterative Symphonic Tunneling for Satisfiability problems (IST-SAT), which solves quantum spin glass optimization problems using high-frequency oscillating transverse fields. IST-SAT operates as a sequence of iterations, in which bitstrings returned from one iteration are used to set spin-dependent phases in oscillating transverse fields in the next iteration. Over several iterations, the novel mechanism of the algorithm steers the system toward the problem ground state. We benchmark IST-SAT on sets of hard MAX-3-XORSAT problem instances with exact state vector simulation, and report polynomial speedups over Trotterized adiabatic quantum computation and the best known semi-greedy classical algorithm. When IST-SAT is seeded with a sufficiently good initial approximation, the algorithm converges to exact solution(s) in a polynomial number of iterations. Our numerical results identify a critical Hamming radius, or quality of initial approximation, where the time-to-solution crosses from exponential to polynomial scaling in problem size. This work proposes IST-SAT a new quantum algorithm, which improves upon solutions obtained from initial classical or quantum optimization algorithms. The steering mechanism we introduce through IST-SAT presents a new path toward achieving quantum advantage in optimization.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Evolution of Bose–Einstein condensate systems beyond the Gross–Pitaevskii equation

While many phenomena in cold atoms and other Bose–Einstein condensate (BEC) systems are often described using the mean-field approaches, understanding the kinetics of BECs requires the inclusion of particle scattering via the collision integral of the quantum Boltzmann equation. A rigorous approach for many problems in the dynamics of the BEC, such as the nucleation of the condensate or the decay of the persistent current, requires, in the presence of factors making a symmetry breaking possible, considering collisions with thermal atoms via the collision integral. These collisions permit the emergence of vorticity or other signatures of long-range order in the nucleation of the BEC or the transfer of angular momentum to thermal atoms in the decay of persistent current, due to corresponding terms in system Hamiltonians. Here, we also discuss the kinetics of spin–orbit-coupled BEC. The kinetic equation for the particle spin density matrix is derived. Numerical simulations demonstrate significant effects of the collision integral on the dynamics of the spin–orbit-coupled BEC upon quenching of the Raman coupling that generates synthetic electric and magnetic fields.

Bose–Einstein condensation↗
DOE OSTI2020

Theoretical Analysis of Fractional Viscoelastic Flow in Circular Pipes: General Solutions

Fractional calculus is a relatively old yet emerging field of mathematics with the widest range of engineering and biomedical applications. Despite being an incredibly powerful tool, it, however, requires promotion in the engineering community. Rheology is undoubtedly one of the fields where fractional calculus has become an integral part of cutting-edge research. There exists extensive literature on the theoretical, experimental, and numerical treatment of various fractional viscoelastic flows in constraint geometries. However, the general theoretical approach that unites several most commonly used models is missing. Here we present exact analytical solutions for fractional viscoelastic flow in a circular pipe. We find velocity profiles and shear stresses for fractional Maxwell, Kelvin–Voigt, Zener, Poynting–Thomson, and Burgers models. The dynamics of these quantities are studied with respect to normalized pipe radius, fractional orders, and elastic moduli ratio. Three different types of behavior are identified: monotonic increase, resonant, and aperiodic oscillations. The models developed are applicable in the widest material range and allow for the alteration of the balance between viscous and elastic properties of the materials.

Gritsenko, Dmitry (ORCID:0000000248181512)↗
DOE OSTI2024

A Parametric Study of the SASI Comparing General Relativistic and Nonrelativistic Treatments

We present numerical results from a parameter study of the standing accretion shock instability (SASI), investigating the impact of general relativity (GR) on the dynamics. Using GR hydrodynamics with GR gravity, and nonrelativistic (NR) hydrodynamics with Newtonian gravity, in an idealized model setting, we vary the initial radius of the shock, and by varying its mass and radius in concert, the proto-neutron star compactness. We investigate four compactnesses expected in a post-bounce core-collapse supernova (CCSN). We find that GR leads to a longer SASI oscillation period, with ratios between the GR and NR cases as large as 1.29 for the highest-compactness suite. We also find that GR leads to a slower SASI growth rate, with ratios between the GR and NR cases as low as 0.47 for the highest-compactness suite. We discuss implications of our results for CCSN simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

The ICON-A model for direct QBO simulations on GPUs (version icon-cscs:baf28a514)

Abstract. Classical numerical models for the global atmosphere, as used for numerical weather forecasting or climate research, have been developed for conventional central processing unit (CPU) architectures. This hinders the employment of such models on current top-performing supercomputers, which achieve their computing power with hybrid architectures, mostly using graphics processing units (GPUs). Thus also scientific applications of such models are restricted to the lesser computer power of CPUs. Here we present the development of a GPU-enabled version of the ICON atmosphere model (ICON-A), motivated by a research project on the quasi-biennial oscillation (QBO), a global-scale wind oscillation in the equatorial stratosphere that depends on a broad spectrum of atmospheric waves, which originates from tropical deep convection. Resolving the relevant scales, from a few kilometers to the size of the globe, is a formidable computational problem, which can only be realized now on top-performing supercomputers. This motivated porting ICON-A, in the specific configuration needed for the research project, in a first step to the GPU architecture of the Piz Daint computer at the Swiss National Supercomputing Centre and in a second step to the JUWELS Booster computer at the Forschungszentrum Jülich. On Piz Daint, the ported code achieves a single-node GPU vs. CPU speedup factor of 6.4 and allows for global experiments at a horizontal resolution of 5 km on 1024 computing nodes with 1 GPU per node with a turnover of 48 simulated days per day. On JUWELS Booster, the more modern hardware in combination with an upgraded code base allows for simulations at the same resolution on 128 computing nodes with 4 GPUs per node and a turnover of 133 simulated days per day. Additionally, the code still remains functional on CPUs, as is demonstrated by additional experiments on the Levante compute system at the German Climate Computing Center. While the application shows good weak scaling over the tested 16-fold increase in grid size and node count, making also higher resolved global simulations possible, the strong scaling on GPUs is relatively poor, which limits the options to increase turnover with more nodes. Initial experiments demonstrate that the ICON-A model can simulate downward-propagating QBO jets, which are driven by wave–mean flow interaction.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Two-dimensional simulations of internal gravity waves in a $\mathrm{5M_⊙}$ zero-age-main-sequence model

Main-sequence intermediate-mass stars present a radiative envelope that supports internal gravity waves (IGWs). Excited at the boundary with the convective core, IGWs propagate towards the stellar surface and are suspected to impact physical processes such as rotation and chemical mixing. Using the fully compressible time-implicit code MUSIC, we study IGWs in two-dimensional simulations of a zero-age-main-sequence 5 solar mass star model up to 91 percent of the stellar radius with different luminosity and radiative diffusivity enhancements. Our results show that low-frequency waves excited by core convection are strongly impacted by radiative effects as they propagate. This impact depends on the radial profile of radiative diffusivity which increases by almost 5 orders of magnitude between the centre of the star and the top of the simulation domain. In the upper layers of the simulation domain, we observe an increase of the temperature. Our study suggests that this is due to heat added in these layers by IGWs damped by radiative diffusion. We show that non-linear effects linked to large amplitude IGWs may be relevant just above the convective core. Both these effects are intensified by the artificial enhancement of the luminosity and radiative diffusivity, with enhancement factors up to 10 4 times the realistic values. Our results also highlight that direct comparison between numerical simulations with enhanced luminosity and observations must be made with caution. Finally, our work suggests that thermal effects linked to the damping of IGWs could have a non-negligible impact on stellar structure.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Bayesian Estimation of Oscillator Parameters: Toward Anomaly Detection and Cyber-Physical System Security

Cyber-physical system security presents unique challenges to conventional measurement science and technology. Anomaly detection in software-assisted physical systems, such as those employed in additive manufacturing or in DNA synthesis, is often hampered by the limited available parameter space of the underlying mechanism that is transducing the anomaly. As a result, the formulation of anomaly detection for such systems often leads to inverse or ill-posed problems, requiring statistical treatments. Here, we present Bayesian inference of unknown parameters associated with a generic actuator considered as a representative vital element of a cyber-physical system. Via a series of experimental input-output measurements, a transfer function for the actuator is obtained numerically, which serves as our model for the proposed method. Linear, nonlinear, and delayed dynamics may be assumed for the actuator response. By devising a code-based malicious signal, we study the efficacy of Bayesian inference for its potential to produce a detection, including uncertainty quantification, with a remarkably small number of input data points. Our approach should be adaptable to a variety of real-time cyber-physical anomaly detection scenarios.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Collisionless zonal-flow dynamics in quasisymmetric stellarators

The linear collisionless plasma response to a zonal-density perturbation in quasisymmetric stellarators is studied, including the geodesic-acoustic-mode oscillations and the Rosenbluth–Hinton residual flow. While the geodesic-acoustic-mode oscillations in quasiaxisymmetric configurations are similar to tokamaks, they become non-existent in quasi-helically symmetric configurations when the effective safety factor in helical-angle coordinates is small. Compared with concentric-circular tokamaks, the Rosenbluth–Hinton residual is also found to be multiplied by a geometric factor $\mathcal {C}$ that arises from the flux-surface-averaged classical polarization. Using the near-axis-expansion framework, we derive an analytic expression for $\mathcal {C}$ , which varies significantly among different configurations. These analytic results are compared with numerical simulation results from the global gyrokinetic particle-in-cell code GTC, and good agreement with the theoretical Rosenbluth–Hinton residual level is achieved when the quasisymmetry error is small enough.

Zhu, Hongxuan (ORCID:0000000198446972)↗
DOE OSTI2022

Optimization Studies for the Long-Baseline Neutrino Facility at Fermilab

The Deep Underground Neutrino Experiment and Long-Baseline Neutrino Facility (DUNE-LBNF) are under development at Fermilab since early 2010s [1]. At present, the work is being performed towards a comprehensive review conducted by US Department of Energy (DOE)-the Critical Decision 2 (CD-2)-that is planned to take place in the middle of 2022. The primary scientific objectives of DUNE are to carry out a comprehensive investigation of neutrino oscillations to test CP violation in the lepton sector, determine the ordering of the neutrino masses, and to test the three-neutrino paradigm (electron, muon and tau neutrino). The LBNF will provide a 120-GeV proton beam on a neutrino production target utilizing a new 800-MeV superconducting Linac which is expected to be completed in 2027 [2]. The neutrino beamline, which utilizes a target and horn systems, decay pipe, hadron absorber and other systems, is a core component of the LBNF. At present-as a result of numerous iterations-there exists an optimized design with a 1.5-m graphite target and focusing system consisting of three horns. Further optimization energy deposition and radiological calculations are performed towards CD-2 and beyond. This paper describes results of the most recent MARS15 [3] optimization studies.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Experimental Investigation of Elastodynamic Nonlinear Response of Dry Intact, Fractured and Saturated Rock

Nonlinear elastodynamic response of fractured rocks carries crucial information on fracture features that can be exploited to forecast flow properties, friction constitutive behavior and poromechanical response. Well-controlled laboratory experiments are designed to measure the nonlinear elastodynamic response of Westerly granite in three states: dry intact, dry fractured and saturated fractured. We study the effect of fracturing and saturation in modifying the elastodynamic response of the rock. Each sample is tested at a normal stress level of 15 MPa. We measure the elastodynamic response of an intact L-shaped sample of Westerly Granite subjected to normal stress oscillations of prescribed amplitudes (0.2–1 MPa) and frequencies (0.1, 1, 10 Hz). Ultrasonic waves transmitted across the sample are used to monitor the evolution of wave velocity before, during and after dynamic stressing. Furthermore, the nonlinearity of the elastodynamic response is measured in terms of: (1) the offset in normalized wave velocity; (2) the amplitude of wave velocity fluctuation during the oscillations; and (3) recovery rate of the wave velocity post-oscillation. We observe that the three nonlinearity parameters show a similar trend. Irrespective of the parameter, the nonlinearity measures higher for sample in dry-intact condition than that for dry-fractured and the saturated-fractured sample exhibits smaller nonlinearity than the dry-fractured sample. As expected, the saturated sample exhibits less nonlinearity than the dry intact and fractured samples due to the presence of interstitial fluid and the resulting increased interface stiffness. Conversely, the dry intact rock shows a higher nonlinearity than the dry fractured. We use numerical simulations to show that the presence of fracture significantly alters the strain distribution across the bulk of the sample and only the contacting asperities are highly strained, thus resulting in a decrease in the measured elastodynamic nonlinearity.

58 GEOSCIENCES↗
DOE OSTI2023

Qubit-Oscillator Concatenated Codes: Decoding Formalism and Code Comparison

Concatenating bosonic error-correcting codes with qubit codes can substantially boost the errorcorrecting power of the original qubit codes. It is not clear how to concatenate optimally, given that there are several bosonic codes and concatenation schemes to choose from, including the recently discovered Gottesman-Kitaev-Preskill (GKP) – stabilizer codes [Phys. Rev. Lett. 125, 080503 (2020)] that allow protection of a logical bosonic mode from fluctuations of the conjugate variables of the mode. We develop efficient maximum-likelihood decoders for and analyze the performance of three different concatenations of codes taken from the following set: qubit stabilizer codes, analog or Gaussian stabilizer codes, GKP codes, and GKP-stabilizer codes. We benchmark decoder performance against additive Gaussian white noise, corroborating our numerics with analytical calculations. We observe that the concatenation involving GKP-stabilizer codes outperforms the more conventional concatenation of a qubit stabilizer code with a GKP code in some cases. We also propose a GKP-stabilizer code that suppresses fluctuations in both conjugate variables without extra quadrature squeezing and formulate qudit versions of GKP-stabilizer codes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Yu-Shiba-Rusinov states and ordering of magnetic impurities near the boundary of a superconducting nanowire

We theoretically study the spectrum induced by one and two magnetic impurities near the boundary of a one-dimensional nanowire in proximity to a conventional s-wave superconductor and extract the ground-state magnetic configuration. We show that the energies of the subgap states, supported by the magnetic impurities, are strongly affected by the boundary for distances less than the superconducting coherence length. In particular, when the impurity is moved towards the boundary, multiple quantum phase transitions periodically occur in which the parity of the superconducting condensate oscillates between even and odd. We find that the magnetic ground-state configuration of two magnetic impurities depends not only on the distance between them, but also explicitly on their distance away from the boundary of the nanowire. As a consequence, the magnetic ground state can switch from ferromagnetic to antiferromagnetic while keeping the interimpurity distance unaltered by simultaneously moving both impurities away from the boundary. The ground-state magnetic configuration of two impurities is found analytically in the weak coupling regime and exactly for an arbitrary impurity coupling strength using numerical tight-binding simulations.

1-dimensional systems↗