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 181 records · Page 10

Experiments and Simulations on the Turbulent, Rarefaction Wave Driven Rayleigh–Taylor Instability

Experiments were performed to observe the growth of the turbulent, Rayleigh–Taylor unstable mixing layer generated between air and SF 6 , with an Atwood number of A=(ρ 2 –ρ 1 )/(ρ 2 +ρ 1 ) = 0.64, where ρ 1 and ρ 2 are the densities of air and SF 6 , respectively. A nonconstant acceleration with an average value of 2300g 0 , where g 0 is the acceleration due to gravity, was generated by interaction of the interface between the two gases with a rarefaction wave. Three-dimensional, multimode perturbations were generated on the diffuse interface, with a diffusion layer thickness of δ = 3.6 mm, using a membraneless vertical oscillation technique, and 20 experiments were performed to establish a statistical ensemble. The average perturbation from this ensemble was extracted and used as input for a numerical simulation using the Lawrence Livermore National Laboratory (LLNL) Miranda code. Good qualitative agreement between the experiment and simulation was observed, while quantitative agreement was best at early to intermediate times. Several methods were used to extract the turbulent growth constant α from experiments and simulations while accounting for time varying acceleration. Experimental, average bubble and spike asymptotic self-similar growth rate values range from α = 0.022 to α = 0.032 depending on the method used, and accounting for variable acceleration. Values found from the simulations range from α = 0.024 to α = 0.041. Here, the values of α measured in the experiments are lower than what are typically measured in the literature but are more in line with those found in recent simulations.

42 ENGINEERING↗

NOvA 2020 official data release (13.6E20 neutrino + 12.5E20 antineutrino)

NOvA results of the 2020 joint analysis of nue appearance and numu disappearance data presented in Phys. Rev. D 106, 032004:https://doi.org/10.1103/PhysRevD.106.032004Note: these contained plots are slightly different from the ones shown at NEUTRINO 2020.Exposure: 13.6E20 protons on target, neutrino-enhanced beam12.5E20 protons on target, antineutrino-enhanced beamParameters not constrained by NOvA:- th13 is assumed to be Gaussian, ss2th13 = 0.085 +/- 0.003 (PDG 2019, Phys. Rev. D98, 030001 (2018) and 2019 update, https://pdg.lbl.gov/2019/listings/rpp2019-list-neutrino-mixing.pdf)- ss2th12 = 0.851, dm21=7.53e-5 are fixed at 2019 PDG valuesSignificances are obtained with the Feldman-Cousins procedure.Files contain TGraphs for the following:- ssth23 vs deltaCP 1,2,3 sigma and 90% CL contours, for NO or IO- dmsq32 vs ssth23 1,2,3 sigma and 90% CL contours, for NO or IO- deltaCP profile for NO UO, NO LO, IO LO, IO UO- ssth23 profile for NO, IO- dmsq32 profile for NO UO, NO LO, IO LO, IO UOContour files also include a TMarker with the best-fit parameters (in the NO)An additional file containing the predictions for all channels’ signal and background components computed at the NOvA best-fit oscillation parameters and systematic pull terms. Numu samples include the no-oscillation case as well.Please note that any fits performed with these histograms are not expected to exactly reproduce the official NOvA results as parameterizations of the numerous systematic uncertainties considered in the official fits are not included in this release.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Initialized Earth system prediction from subseasonal to decadal timescales

Initialized Earth system predictions are made by starting a numerical prediction model in a state as consistent as possible to observations, and running it forward in time for up to ten years. Skillful predictions at time slices from subseasonal to seasonal (S2S), seasonal to interannual (S2I) and seasonal to decadal (S2D) offer information useful for various stakeholders, from agriculture to water resource management, and human and infrastructure safety. In this Review, we examine the processes influencing predictability, and discuss estimates of skill across S2S, S2I and S2D timescales. There are encouraging signs that skillful predictions can be made: at S2S timescales, there has been some skill in predicting the Madden-Julian Oscillation and North Atlantic Oscillation; at S2I in predicting the El Niño-Southern Oscillation; and at S2D, in predicting variability in North Atlantic sea surface temperatures. However, challenges remain, and future work must prioritise reducing model error, more effectively communicating forecasts to users, and increasing process and mechanistic understanding that could increase predictive skill and, in turn, confidence. As numerical models progress towards Earth system models, initialized predictions are expanding to include prediction of sea-ice, air pollution, terrestrial and ocean biochemistry which can bring clear benefit to society and various stakeholders.

climate prediction↗

General framework for quantifying dissipation pathways in open quantum systems. III. Off-diagonal subsystem–bath couplings

This paper extends the previously reported theory of dissipation pathways [C. W. Kim and I. Franco, J. Chem. Phys. 160, 214111 (2024)] to incorporate off-diagonal subsystem–bath coupling, which is often required to model molecular systems where the environment directly influences transitions and couplings between subsystem states. We systematically derive master equations for both population transfer and dissipation into individual bath components, for which we also rigorously prove energy conservation and detailed balance. The approach is based on second-order perturbation theory with respect to the subsystem–bath couplings, whose form is not limited to any specific model. The accuracy of the developed method is tested by applying it to diverse model Hamiltonians involving linearly coupled harmonic oscillator baths and comparing the outcomes against the hierarchical equations of motion (HEOM) method. Overall, our method accurately quantifies the contributions of specific bath components to the overall dissipation while significantly reducing the computational cost compared to numerically exact methods such as HEOM, thus offering a path to examine how vibronic interactions steer non-adiabatic processes in realistic chemical systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An LES study of secondary motion and wall shear stresses in a pipe bend

In this study, large-eddy simulations (LES) of a single-phase, turbulent flow in a 90 ° pipe bend are performed at three Reynolds numbers (5300, 27 000, and 45 000) to investigate the correlation between secondary flow motion and wall shear stresses, which is suspected to be a potential mechanism responsible for material erosion. The isothermal flows are validated against available experimental and numerical data first. The snapshot proper orthogonal decomposition (POD) is applied for the medium and high Reynolds number flows to identify the secondary flow motions and the oscillation of the Dean vortices that are found to cause swirl-switching. Distinguished frequencies of the POD time coefficients at Strouhal numbers of 0.25 and 0.28 are identified for Reynolds numbers at 27 000 and 45 000, respectively. Moreover, shear stress on the pipe wall and the associated power spectral density are obtained and shown to have the same oscillating frequency as the swirl-switching.

42 ENGINEERING↗

Interference effects in the formation of the light bullet spectrum under axicon focusing

Specific features of the anti-Stokes band of the light bullet supercontinuum spectrum under filamentation of 1900-nm femtosecond pulses, focused by an axicon in fused silica and calcium fluoride, have been experimentally, numerically, and analytically investigated. It is found that the formation of two or more light bullets in a filament leads to modulation of spectral components’ intensity. The oscillation period of spectral components, defined as the wavelength interval between adjacent intensity minima, decreases with decreasing their wavelength. On a logarithmic scale, the relative change in the oscillation period in the spectrum is proportional to the spectral component wavelength and is independent of the pulse energy and the position of the light bullet formation region in the filament. It is established that the modulation in the spectrum is due to the interference of the broadband supercontinuum emitted by light bullets, and its period decreases with increasing interval between them. (paper)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

An analytic implementation of the IR-resummation for the BAO peak

We develop an analytic method for implementing the IR-resummation of [1], which allows one to correctly and consistently describe the imprint of baryon acoustic oscillations (BAO) on statistical observables in large-scale structure. We show that the final IR-resummed correlation function can be computed analytically without relying on numerical integration, thus allowing for an efficient and accurate use of these predictions on real data in cosmological parameter fitting. In this work we focus on the one-loop correlation function and the BAO peak. We show that, compared with the standard numerical integration method of IR-resummation, the new method is accurate to better than 0.2%, and is quite easily improvable. We also give an approximate resummation scheme which is based on using the linear displacements of a fixed fiducial cosmology, which when combined with the method described above, is about six times faster than the standard numerical integration. Finally, we show that this analytic method is generalizable to higher loop computations.

79 ASTRONOMY AND ASTROPHYSICS↗

Tapering-enhanced high-efficiency THz waveguide oscillator

Using a waveguide in a THz FEL has been shown to maximize the coupling and enable efficient extraction of relativistic electron beam energy from a single passage through a tapered helical undulator. An oscillator configuration can further boost energy extraction above the single-pass limits and open the door toward very high average power THz sources. Embedding the undulator in an oscillator cavity is particularly useful in combination with high repetition rate electron sources, even if at reduced peak brightness, since recycling a fraction of the radiation as an intense seed can compensate for lower single-pass gain. In this paper, we investigate the efficiency scaling of a tapering-enhanced waveguide oscillator, showcasing its capability for frequency-tuning operation and high-efficiency operation at different wavelengths. Using a thermionic-driven beamline equipped with compression elements, numerical start-to-end simulation results indicate a 16% efficiency at 200 GHz and a 2.1% efficiency at 1.5 THz, resulting in kW-level average power out-coupled in radiation pulses with few hundred μ J energy and tens of MW peak power. Published by the American Physical Society 2024

Yang, Y. N. (ORCID:0009000849540151)↗

Where, when, and why: Occurrence of fast-pairwise collective neutrino oscillation in three-dimensional core-collapse supernova models

Fast-pairwise collective neutrino oscillation represents a key uncertainty in the theory of core-collapse supernova (CCSN). Despite the potentially significant impact on CCSN dynamics, it is usually neglected in numerical models of CCSN because of the formidable technical difficulties of self-consistently incorporating this physics. In this paper, we investigate the prospects for the occurrence of fast flavor conversion by diagnosing electron neutrino lepton number (ELN) crossing in more than a dozen state-of-the-art three-dimensional CCSN models. ELN crossings provide a necessary condition for triggering flavor conversion. Although only zeroth and first angular moments are available from the simulations, our new method enables us to look into the angular distributions of neutrinos in momentum space and provide accurate insight into ELN crossings. Our analysis suggests that fast flavor conversion generally occurs in the postshock region of CCSNe, and that explosive models provide more favorable conditions for the flavor conversion than failed CCSNe. Here we also find that there are both common and progenitor-dependent characteristics. Classifying ELN crossings into two types, we analyze the generation mechanism of each case by scrutinizing the neutrino radiation field and matter interactions. We find key ingredients of CCSN dynamics driving the ELN crossings: proto-neutron star convection, asymmetric neutrino emission, neutrino absorptions, and scatterings. This study suggests that we need to accommodate fast flavor conversions in realistic CCSN models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Beehive haloscope for high-mass axion dark matter

We propose a new haloscope geometry that can arbitrarily increase the resonator volume for a given target axion mass. This geometry consists of closely packed, overlapping coaxial cavities operating as a single resonator. While the resonant frequency is still determined by the dimensions of the individual “cells,” the strong interactions between the cells encourage the entire “beehive” to oscillate in phase, a phenomenon expected of tightly coupled harmonic oscillators. This synchronization behavior allows the construction of a singly connected large volume resonator at high frequency by simply increasing the number of the cells. Using direct numerical simulations, we verify the existence of a global eigenmode that has a high (40%) form factor in a 169-element beehive resonator. The resonant frequency of the eigenmode is tunable by moving the center rods laterally in unison. The form factor is very tolerant to dimensional deviations and misalignment as a result of mode hybridization due to strong coupling. The beehive haloscope inherits many appealing properties from the conventional coaxial cavity: a high quality factor, compatibility with a solenoid magnet, and ease of fabrication, tuning, and coupling. We argue that this geometry is an excellent candidate for high-mass axion searches covering the post-inflationary parameter space (>5 GHz).

47 OTHER INSTRUMENTATION↗

A Study on Particle Trajectory Error in Finite-Element Particle-in-Cell Algorithms

Particle-in-cell (PIC) algorithms are widely used for the simulation of kinetic plasmas. PIC algorithms account for the interaction between charged particles in a plasma and the electromagnetic field in ambient space, including self-field interactions. The objective of this article is to study the error in charged particle trajectories present in finite-element (FE)-based PIC algorithms on unstructured meshes. We study how the trajectory error behaves according to the FE mesh resolution and the matrix solver employed in the FE algorithm. The study is performed by considering a trajectory established by a parabolic electric potential and an axial magnetic force acting on the charged particle. Under a proper choice combination of electric and magnetic field parameters, the 2-D particle trajectories comprise closed orbits. Numerical errors cause small perturbations on the trajectories, with cumulative effects. As a result, the resulting orbital trajectories exhibit distortions including spurious apsidal precession. These distortions provide a clear imprint of the numerical error. Here, we also study the numerical error in a quantitative fashion by computing the distance norm between the trajectories obtained by the exact fields and numerical fields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NOvA 2024 official data release (26.61E20 neutrino + 12.5E20 antineutrino)

This data release corresponds to the Bayesian 2024 analysis of NOvA $\nu_e$ appearance and $\nu_{\mu}$ disappearance data, corresponding to analysis described in https://arxiv.org/abs/2509.04361. File `NOvA_2024_data_histograms.root` contains data histograms for all the NOvA data samples. Exposure: * neutrino-enhanced beam: 26.61E20 protons on target * antineutrino-enhanced beam: 12.5E20 protons on target External constraints: * ss2th12=0.851, dm21=7.53e-5 are fixed at 2019 PDG values, with negligible effect on NOvA predictions. * ss2th13 & dm32: * RCDB1D: 1D constraint from Daya Bay for ss2th13, 0.0851+/-0.0024 * RCDB2D: Correlated 2D constraint from Daya Bay on ss2th13 & dm32, available in their official 2023 data release: https://journals.aps.org/prl/abstract/10.1103/PhysRevLett.130.161802#supplemental An additional file containing the predictions (`NOvA_2024_prediction_with_systs_histograms.root`) for all channels’ signal and background components computed at the NOvA best-fit oscillation parameters and systematic pull terms. Numu samples include the no-oscillation case as well. Please note that any fits performed with these histograms are not expected to exactly reproduce the official NOvA results as parameterizations of the numerous systematic uncertainties considered in the official fits are not included in this release. The zip file also contains the 2D credible interval contours, with details in a README.md

Sztuc, Artur [University Coll. London] (ORCID:0000↗

Imaging Bias‐Driven Domain Wall Motion With Scanning Oscillator Piezoresponse Force Microscopy

Understanding ferroelectric domain wall dynamics at the nanoscale across a broad range of timescales requires measuring domain wall position under different applied electric fields. The success of piezoresponse force microscopy (PFM) as a tool to apply local electric fields at different positions and imaging their changing position, together with the information obtained from associated switching spectroscopies has fueled numerous studies of the dynamics of ferroelectric domains to determine the impact of intrinsic parameters such as crystalline order, defects and pinning centers, as well as boundary conditions such as environment. However, the investigation of sub-coercive reversible domain wall vibrational modes requires the development of new tools that enable visualizing domain wall motion under varying applied fields with high temporal and spatial resolution while also accounting for spurious electrostatic effects. Here, scanning oscillator piezoresponse force microscopy extends the investigation of domain wall dynamics to new regimes, providing direct visualization of domain wall position as a function of an external electric field that varies in time and location. In conclusion, this enables studying the energetics of field-driven ferroelectric domain wall motion, which is shown to obey a thermally activated flow regime in the millisecond timescale.

Piezoresponse Force Microscopy↗

Observations and Modeling of Ocean Circulation in the Seychelles Plateau Region

The ocean circulation around and over the Seychelles Plateau is characterized using 35 months of temperature and velocity measurements along with a numerical model. The results here provide the first documented description of the ocean circulation atop the Seychelles Plateau. The Seychelles Plateau is an unusually broad (~200 km), shallow (~50 m) plateau, dropping off steeply to the abyss. It is situated in a dynamic location (3.5-5.5S, 54-57E) in the south-western tropical Indian Ocean where northwesterly winds are present during austral summer and become southeasterly in austral winter, following the reversal of the western Indian Ocean monsoon winds. Measurements around the Inner Islands, on the Seychelles Plateau, have been carried out since 2015. Velocity measurements show that most of the depth-averaged current variance on the Seychelles Plateau arises from near-inertial oscillations and lower-frequency variability. Lower-frequency variability encompasses seasonal and intraseasonal variability, the latter of which includes the effects of mixed Rossby-gravity waves and mesoscale eddies. Here, a global 0.1-degree numerical ocean simulation is used in conjunction with these observations to describe the regional circulation around the Seychelles Plateau. Atop the Seychelles Plateau, circulation is dominated by ageostrophic processes consistent with Ekman dynamics, while around the Seychelles Plateau, both geostrophic and ageostrophic processes are important and vary seasonally. Stratification responds to the sea surface height semi-annual signal which is due to Ekman pumping-driven upwelling (related to the Seychelles-Chagos Thermocline Ridge) and the arrival of an annual downwelling Rossby wave.

58 GEOSCIENCES↗

A Geometric Volume of Fluid-Based Multiphase Flow Solver Extension to the Reacting Flow Solver, PeleLM

A new algorithm is presented to simulate multiphase flows with surface tension in a pathway for spray combustion simulation. The algorithm combines capabilities from two open-source packages, including the interface reconstruction library (IRL), a library of computational geometry routines to enable the volume of fluid (VOF) method, and PeleLM, a solver for the reacting Navier-Stokes equations. Additionally, surface tension is implemented using the continuum surface force (CSF) model with an improved height function technique in the volume fraction field. Spurious errors in volume fraction arising from our combined strategy are corrected through a topology-based method that improves both numerical stability and accuracy. Multiple validation simulations are conducted, including (i) translations and rotations of Zalesak's disk, (ii) a stationary circular droplet with surface tension, (iii) an oscillating elliptical droplet, and (iv) three-dimensional deformation of a spherical droplet. Results indicate that the combined scheme retains the favorable properties of each of the component algorithms.

42 ENGINEERING↗

Large-eddy simulation study on cycle-to-cycle variation of knocking combustion in a spark-ignition engine

The cycle-to-cycle variation in the knock intensity is commonly encountered under abnormal combustion conditions. The severity of these abnormal combustion events can vary significantly, and the efficiency of engines at high loads is limited in practice by heavy knocking phenomena. Since, a thorough analysis of such recurrent but non-cyclic phenomena via experiments alone becomes highly cumbersome, in the present work, a multi-cycle large-eddy simulation study was performed to quantitatively predict cyclic variability in the combustion process and cyclic knock intensity variability in a direct injection spark-ignition engine. To account for the turbulence-chemistry interaction effects on flame propagation, the G-equation combustion model was used. Detailed chemistry was solved outside the flame front with a toluene primary reference fuel skeletal kinetic mechanism. For both the mild knock and heavy knock conditions, the numerical results were validated against experimental measurements. Based on the simulation results, a correlation analysis was performed considering combustion phasing, peak cylinder pressure and maximum amplitude of pressure oscillation. Furthermore, a detailed three-dimensional spatial analysis illustrated the evolution of auto-ignition kernel development and propagation of pressure waves during knocking combustion for three typical cycles with different knock intensities. In this process, it was found that an early occurrence of auto-ignition in the end gas was prone to high knock intensity. Although multiple auto-ignition kernels were observed in different cycles, the degree of coupling between chemical heat release and pressure waves varied, thereby leading to different maximum amplitude of pressure oscillation values.

42 ENGINEERING↗

An asymptotic approach for the statistical thermodynamics of certain model systems

In classical statistical thermodynamics, calculating the configuration integral is both vital and elusive. Analytic relations for configuration integrals are desirable for modeling purposes, but it is typically impossible to obtain them. Certain systems become analytically tractable after replacing steep potential energies with harmonic potentials or athermal rigid constraints, but these approximations are often inadequate, especially when modeling the stretching of molecules. It is therefore necessary to develop a systematic approach to improve upon the approximations provided by these reference systems. Here, a general asymptotic approach is introduced, where the configuration integral for the full system is obtained in terms of that of the reference system and several corrections. This asymptotic approach is first demonstrated using the simple example of a classical three-dimensional oscillator. Next, the approach is applied to modeling the stretching of single polymer chains and to modeling thermally assisted crack growth, where results are verified with respect to numerical calculations. Overall, this asymptotic approach is a valid and effective tool for statistical thermodynamics in general.

Buche, Michael Robert↗

Simulations of vertical displacement oscillatory modes and global Alfvén Eigenmodes in JET geometry

Vertical Displacement Oscillatory Modes (VDOM), with frequency in the Alfvén range, are natural modes of oscillation of magnetically confined laboratory plasmas with elongated cross-section. These axisymmetric modes arise from the interaction between the plasma current, which is in equilibrium with currents flowing in external coils, and perturbed currents induced on a nearby conducting wall. The restoring force exerted by these perturbed currents on the vertical motion of the plasma column leads to its oscillatory behavior. An analytic model for VDOM was proposed based on an idealized 'straight tokamak' equilibrium with uniform equilibrium current density. This article introduces the first numerical simulations of VDOM in a realistic JET tokamak configuration, using the extended-MHD code NIMROD and drawing comparisons with Global Alfvén Eigenmodes (GAE). The results show qualitative agreement with analytic predictions regarding mode frequency and radial structure, supporting the identification of VDOM as a fundamental oscillation mode in tokamak plasmas. VDOM and GAE are modeled in a representative JET discharge, where axisymmetric perturbations with toroidal mode number n = 0 driven unstable by fast ions were observed. The two modes are examined separately using a forced oscillator within the NIMROD code, which enables a comparison of their characteristics and helps identify the experimentally observed mode possibly as a GAE.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗