Analytical approximation and numerical studies of one-dimensional elliptic equation with random coefficients
Not Available
Xu, Zhijie (ORCID:0000000304594531)↗
SEARCH · Search NASA
Results for “explicit correlation”
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.
Perspective: How good is DFT for water?
Kohn-Sham density functional theory (DFT) has become established as an indispensable tool for investigating aqueous systems of all kinds, including those important in chemistry, surface science, biology, and the earth sciences. Nevertheless, many widely used approximations for the exchange-correlation (XC) functional describe the properties of pure water systems with an accuracy that is not fully satisfactory. The explicit inclusion of dispersion interactions generally improves the description, but there remain large disagreements between the predictions of different dispersion-inclusive methods. We present here a review of DFT work on water clusters, ice structures, and liquid water, with the aim of elucidating how the strengths and weaknesses of different XC approximations manifest themselves across this variety of water systems. Our review highlights the crucial role of dispersion in describing the delicate balance between compact and extended structures of many different water systems, including the liquid. By referring to a wide range of published work, we argue that the correct description of exchange-overlap interactions is also extremely important, so that the choice of semi-local or hybrid functional employed in dispersion-inclusive methods is crucial. The origins and consequences of beyond-2-body errors of approximate XC functionals are noted, and we also discuss the substantial differences between different representations of dispersion. We propose a simple numerical scoring system that rates the performance of different XC functionals in describing water systems, and we suggest possible future developments.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
The origin of the Stokes–Einstein relation in simple dense liquids
Here, we investigate the origin of the universal relation between structural relaxation and diffusion in simple dense liquids, known as the Stokes–Einstein (SE) relation. The fact that this relation, originally derived from a hydrodynamic model of a macroscopic particle in a viscous medium, can describe the microscopic-scale liquid dynamics still eludes understanding. We introduce a new universal measure of structural relaxation in a system of N identical particles based on an explicit decomposition of the configuration space into N! congruent convex polyhedra. This measure makes it possible to quantify the correlation between two distinct particle configurations in terms of their minimal Euclidean distance, optimized with respect to particle permutations. Using this measure alongside a model of independent random walkers under the single-occupancy constraint, we derive a master equation that quantifies the SE relation. It allows us to demonstrate that the universal relation between structural relaxation and diffusion in simple dense liquids is caused by two conditions: (a) the confinement of the dominant density fluctuations to the first coordination shell, manifested by de Gennes narrowing, and (b) Gaussianity of the diffusion process; the former is shown to be violated in low-density fluids, and the latter is known to be violated in supercooled liquids.
Physics - Condensed matter physics↗
Optimal superdense coding over memory channels
We study the superdense coding capacity in the presence of quantum channels with correlated noise. We investigate both the cases of unitary and nonunitary encoding. Pauli channels for arbitrary dimensions are treated explicitly. The superdense coding capacity for some special channels and resource states is derived for unitary encoding. We also provide an example of a memory channel where nonunitary encoding leads to an improvement in the superdense coding capacity.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Fluctuation dynamics near the QCD critical point
The evolution of nonhydrodynamic slow processes near the QCD critical point is explored with the novel hydro+ framework, which extends the conventional hydrodynamic description by coupling it to additional explicitly evolving slow modes describing long wavelength fluctuations. Here, their slow relaxation is controlled by critical behavior of the correlation length and is independent from gradients of matter density and pressure that control the evolution of the hydrodynamic quantities. In this exploratory study, we follow the evolution of the slow modes on top of a simplified QCD matter background, allowing us to clearly distinguish and study, both separately and in combination, the main effects controlling the dynamics of critical slow modes. In particular, we show how the evolution of the slow modes depends on their wave number, the expansion of and advection by the fluid background, and the behavior of the correlation length. Nonequilibrium contributions from the slow modes to bulk matter properties that affect the bulk dynamics (entropy, pressure, temperature, and chemical potential) are discussed and found to be small.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Observables for possible QGP signatures in central pp collisions
Proton–proton ( pp ) data show collective effects, such as long-range azimuthal correlations and strangeness enhancement, which are similar to phenomenology observed in heavy ion collisions. Using simulations with and without explicit existing models of collective effects, we explore new ways to probe pp collisions at high multiplicity, in order to suggest measurements that could help identify the similarities and differences between large- and small-scale collective effects. In particular, we focus on the properties of jets produced in ultra-central pp collisions in association with a Z boson. We consider observables such as jet energy loss and jet shapes, which could point to the possible existence of an underlying quark-gluon plasma, or other new dynamical effects related to the presence of large hadronic densities.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
A theoretical approach to sound propagation and radiation for ducts with suppressors
The several phenomena involved in theoretical prediction of the far-field sound radiation attenuation from an acoustically lined duct were studied. These include absorption by the suppressor, termination reflections, and far-field radiation. Extensive parametric studies show that the suppressor absorption performance can be correlated with mode cut-off ratio or angle of propagation. The other phenomena can be shown to depend explicitly upon mode cut-off ratio. A complete system can thus be generated which can be used to evaluate aircraft sound suppressors and which can be related to the sound source through the cut-off ratio-acoustic power distribution. Although the method is most fully developed for inlet suppressors, several aft radiated noise phenomena are also discussed. This simplified suppressor design and evaluation method is summarized, the recent improvements in the technique are presented, and areas where further refinement is necessary are discussed. Noise suppressor data from engine experiments are compared with the theoretical calculations.
Rice, E. J.↗
A theoretical approach to sound propagation and radiation for ducts with suppressors
The several phenomena involved in theoretical prediction of the far-field sound radiation attenuation from an acoustically lined duct have been studied. These include absorption by the suppressor, termination reflection, and far-field radiation. Extensive parameteric studies have shown that the suppressor absorption performance can be correlated with mode cut-off ratio or angle of propagation. The other phenomena can be shown to depend explicitly upon mode cut-off ratio. A complete system can thus be generated which can be used to evaluate aircraft sound suppressors and which can be related to the sound source through the cut-off ratio-acoustic power distribution. Although the method is most fully developed for inlet suppressors, several aft radiated noise phenomena will also be discussed. This paper summarized this simplified suppressor design and evaluation method, presents the recent improvements in the technique and discusses areas where further refinement is necessary. Noise suppressor data from engine experiments are compared with the theoretical calculations.
Rice, E. J.↗
Advanced Computational Modeling of High-Level Waste Vitrification at the Hanford Site
The U.S. Department of Energy (DOE) has selected vitrification for stabilizing legacy tank waste at the Hanford site, where radioactive waste from plutonium production was historically stored in underground tanks. This waste will be separated into low-activity waste (LAW) and high-level waste (HLW) fractions and processed at the Waste Treatment and Immobilization Plant (WTP). At WTP, glass melters are used for the vitrification of radioactive tank waste, transforming it into a stable borosilicate glass form for safe long-term storage. The melter vessel is constructed from highly durable and heat-resistant materials, where the vitrification process occurs. The main regions that are modeled are the melt pool, plenum, cold cap, riser/discharge chamber, and surrounding structure with insulation layers. Forced convection induced by air bubblers at the base of the melter ensure uniform temperature distribution and provide heat to the cold cap layer. The cold cap is a region of reacting batch feed that floats on top of the molten glass and is where the batch-to-glass reactions occur. Joule heating provided by electrodes mounted along the vertical walls of the melter and immersed directly in the glass, generates the necessary heat for the net endothermic conversion processes that occur in the cold cap. The high temperatures, radioactivity, and opaque nature of the glass prevent direct observation inside the melters. Therefore, computational models are essential for providing insight into factors that affect melter throughput. Thermocouples in the plenum provide operators with plenum temperature measurements. Operational adjustments include bubbling rate, voltage supplied to the electrodes, feed adjustments, and glass removal rate. Different computational fluid dynamics (CFD) models have been developed, each serving a specific purpose. There are CFD models of different scale melters, as well as models that capture the two-phase flow interfaces of rising bubbles in the molten glass or models with a simplified molten glass region so that the surrounding structure and plenum can be feasibly incorporated. Pilot-scale melter models have been developed to serve as validation of the methods employed in the simulation of the full-scale WTP melters. Models incorporating resolved bubbling are used to develop momentum source terms to implement into a single phase, multi-region, steady-state flow model that is being validated by measured process parameters such as glass production rate, voltage, input power, plenum temperatures, etc. The resolved bubbling model uses the multiphase volume of fluid approach to model the system with a high-resolution interface capturing scheme to maintain sharp interfaces between the molten glass and the air phase. The suite of CFD models is continually being improved to incorporate more realistic physics and achieve faster turnaround time. For example, an incremental controller is implemented to automatically adjust electrode voltage within the simulation to a molten glass set point temperature of 1150°C. Newer models feature improved meshes to ensure conformal meshes between regions and eliminate unnecessary mesh refinement in areas that are not of interest (such as boundary layers in offgas ports). Instead of explicitly modeling the structural, refractory, and insulation layers of the melter, a thermal resistance approach is used with published correlations used for boundary conditions. The development of robust and efficient CFD models will be instrumental in enabling the WTP to successfully fulfill its mission of safely stabilizing legacy nuclear waste.
12 - MGMT OF RADIOACTIVE AND NON-RADIOACTIVE WASTE↗
Density cumulant functional theory from a unitary transformation: N-representability, three-particle correlation effects, and application to O{sub 4}{sup +}
A new approach to density cumulant functional theory is developed that derives density cumulant N-representability conditions from an approximate Fock space unitary transformation. We present explicit equations for the third- and fourth-order two-particle cumulant N-representability, as well as the second-order contributions that depend on the connected three-particle density cumulant. These conditions are used to formulate the ODC-13 method and the non-iterative (λ{sub 3}) correction that employ an incomplete description of the fourth-order two-particle cumulant N-representability and the second-order three-particle correlation effects, respectively. We perform an analysis of the ODC-13 N-representability description for the dissociation of H{sub 2} and apply the ODC-13 method and the (λ{sub 3}) correction to diatomic molecules with multiple bond character and the symmetry-breaking tetraoxygen cation (O{sub 4}{sup +}). For the O{sub 4}{sup +} molecule, the vibrational frequencies of the ODC-13(λ{sub 3}) method do not exhibit spatial symmetry breaking and are in a good agreement with the recent infrared photodissociation experiment. We report the O{sub 4}{sup +} equilibrium structure, harmonic frequencies, and dissociation energy computed using ODC-13(λ{sub 3}) with a diffuse, core-correlated aug-cc-pCVTZ basis set.
74 ATOMIC AND MOLECULAR PHYSICS↗
NONLINEAR BEHAVIOR OF BARYON ACOUSTIC OSCILLATIONS FROM THE ZEL'DOVICH APPROXIMATION USING A NON-FOURIER PERTURBATION APPROACH
Baryon acoustic oscillations are an excellent technique to constrain the properties of dark energy in the universe. In order to accurately characterize the dark energy equation of state, we must understand the effects of both the nonlinearities and redshift space distortions on the location and shape of the acoustic peak. In this paper, we consider these effects using the Zel'dovich approximation and a novel approach to second-order perturbation theory. The second-order term of the Zel'dovich power spectrum is built from convolutions of the linear power spectrum with polynomial kernels in Fourier space, suggesting that the corresponding term of the Zel'dovich correlation function can be written as a sum of quadratic products of a broader class of correlation functions, expressed through simple spherical Bessel transforms of the linear power spectrum. We show how to systematically perform such a computation. We explicitly prove that our result is the Fourier transform of the Zel'dovich power spectrum and compare our expressions to numerical simulations. Finally, we highlight the advantages of writing the nonlinear expansion in configuration space, as this calculation is easily extended to redshift space, and the higher-order terms are mathematically simpler than their Fourier counterparts.
79 ASTRONOMY AND ASTROPHYSICS↗
Cluster correlations in dense matter and equation of state
Correlations in strongly interacting matter at subsaturation densities lead to the formation of clusters and the appearance of phases transitions with a change of thermodynamic properties and chemical composition. These features can be described in a generalized relativistic density functional approach using clusters as explicit degrees of freedom with medium-dependent properties. The model is constructed in order to provide equations of state for astrophysical applications. It can be adapted also to nuclear structure calculations with cluster correlations on the surface of heavy nuclei. The appearance of α-particles modifies the neutron skin thickness of neutron-rich nuclei and affects the correlation with the density dependence of the symmetry energy.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Calculation of Friction Factors and Nusselt Numbers for Twisted Elliptical Tube Heat Exchangers Using Nek5000
Twisted elliptical tube geometries for heat exchanger design are expected to offer significant enhancement in heat transfer with a marginal increase in frictional losses. The presented work focuses on computational fluid dynamics (CFD) simulations of heated molten salt flows through twisted elliptical tube geometries at low modified Froude numbers. In this work, the objectives are to evaluate the available correlations at low Froude numbers and to determine the impact of using a small gap to resolve contact points. The spectral element CFD code Nek5000 was used for all simulations, which were performed in periodic domains of triangular and square unit cells surrounding a single tube through a complete twist using an explicit filtering large eddy simulation (LES) method. The simulation methodology was validated by comparing to the experimental data of Tan et al. Excellent agreement was observed for this comparison. The Dzyubenko correlation for transitional flow was evaluated for a range of Froude numbers. This correlation demonstrated good agreement with LES results at high modified Froude number. As Froude number was decreased toward the bounds of the correlation, the agreement worsened. Cases were then simulated at low Froude number, testing the effects of tube spacing. It was determined that the laminar case for the square unit cell is the most affected by increasing gap size, which should be minimized to mitigate the effect. In the triangular unit cell, the laminar flow regime is also more significantly impacted by increasing gap size compared to the turbulent flow regime which was only marginally impacted.
42 ENGINEERING↗
Ancilla-assisted sequential approximation of nonlocal unitary operations
We consider the recently proposed (no-go) theorem of Lamata et al.[Phys. Rev. Lett. 101, 180506 (2008)] on the impossibility of sequential implementation of global unitary operations with the aid of an itinerant ancillary system and view the claim within the language of Kraus representation. By virtue of an extremely useful tool for analyzing entanglement properties of quantum operations, namely, operator-Schmidt decomposition, we provide alternative proof to the no-go theorem and also study the role of initial correlations between the qubits and ancilla in sequential preparation of unitary entanglers. Despite the negative response from the no-go theorem, we demonstrate explicitly how the matrix-product operator (MPO) formalism provides a flexible structure to develop protocols for sequential implementation of such entanglers with an optimal fidelity. The proposed numerical technique, which we call variational matrix-product operator (VMPO), offers a computationally efficient tool for characterizing the ''globalness'' and entangling capabilities of nonlocal unitary operations.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Can the ring polymer molecular dynamics method be interpreted as real time quantum dynamics?
The ring polymer molecular dynamics (RPMD) method has gained popularity in recent years as a simple approximation for calculating real time quantum correlation functions in condensed media. However, the extent to which RPMD captures real dynamical quantum effects and why it fails under certain situations have not been clearly understood. Addressing this issue has been difficult in the absence of a genuine justification for the RPMD algorithm starting from the quantum Liouville equation. To this end, a new and exact path integral formalism for the calculation of real time quantum correlation functions is presented in this work, which can serve as a rigorous foundation for the analysis of the RPMD method as well as providing an alternative derivation of the well established centroid molecular dynamics method. The new formalism utilizes the cyclic symmetry of the imaginary time path integral in the most general sense and enables the expression of Kubo-transformed quantum time correlation functions as that of physical observables pre-averaged over the imaginary time path. Upon filtering with a centroid constraint function, the formulation results in the centroid dynamics formalism. Upon filtering with the position representation of the imaginary time path integral, we obtain an exact quantum dynamics formalism involving the same variables as the RPMD method. The analysis of the RPMD approximation based on this approach clarifies that an explicit quantum dynamical justification does not exist for the use of the ring polymer harmonic potential term (imaginary time kinetic energy) as implemented in the RPMD method. It is analyzed why this can cause substantial errors in nonlinear correlation functions of harmonic oscillators. Such errors can be significant for general correlation functions of anharmonic systems. We also demonstrate that the short time accuracy of the exact path integral limit of RPMD is of lower order than those for finite discretization of path. The present quantum dynamics formulation also serves as the basis for developing new quantum dynamical methods that utilize the cyclic nature of the imaginary time path integral.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Operator constraints for twist-3 functions and Lorentz invariance properties of twist-3 observables
We investigate the behavior under Lorentz transformations of perturbative coefficient functions in a collinear twist-3 formalism relevant for high-energy observables including transverse polarization of hadrons. We argue that those perturbative coefficient functions can, a priori , acquire quite different yet Lorentz-invariant forms in various frames. This somewhat surprising difference can be traced back to a general dependence of the perturbative coefficient functions on light cone vectors which are introduced by the twist-3 factorization formulas and which are frame-dependent. One can remove this spurious frame dependence by invoking so-called Lorentz invariance relations (LIRs) between twist-3 parton correlation functions. Some of those relations for twist-3 distribution functions were discussed in the literature before. In this paper we derive the corresponding LIRs for twist-3 fragmentation functions. We explicitly demonstrate that these LIRs remove the light cone vector dependence by considering transverse spin observables in the single-inclusive production of hadrons in lepton-nucleon collisions, ℓN→hX. Furthermore, with the LIRs in hand, we also show that twist-3 observables in general can be written solely in terms of three-parton correlation functions.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Toward the analytic bootstrap of energy correlators
In this paper, we present a framework for the analytic bootstrap of three-point energy correlators, a crucial observable in $\mathcal{N}$ = 4 super Yang-Mills theory and quantum chromodynamics (QCD). Our approach combines spherical contour techniques, general physical constraints such as pole cancellations, and power correction data in the singular limits to determine its analytic expression. In contrast to previous bootstrap studies restricted to scattering amplitudes for supersymmetric theories, our framework makes use of the properties of Feynman integrals, marking a significant step toward bootstrapping realistic QCD observables. Using this method, we derive analytic expressions for leading-order three-point energy correlators in the multi-collinear limit with equal and unequal energy weights, where the latter are crucial ingredients for projected N -point energy correlators. We also apply the recently developed technique of analytic regression with lattice reduction as a way to bypass needing explicit expressions for the singular limits. Bridging theoretical advances in scattering amplitudes with the renewed interest in weighted cross-sections, our work opens the door to precision tests of QCD dynamics through analytic event-shape predictions.
jet substructure↗
Notes on the D = 11 pure spinor superparticle
The D = 11 pure spinor superparticle has been shown to describe linearized D = 11supergravityinamanifestlycovariantway. Anumberofauthorshaveproposedthat its correlation functions be used to compute amplitudes. The use of the scalar structure of the eleven-dimensional pure spinor top cohomology introduces a natural measure for computing such correlation functions. This prescription requires the construction of ghost number one and zero vertex operators. In these notes, we construct explicitly a ghost number one vertex operator but show the incompatibiliy of a ghost number zero vertex operator satisfying a standard descent equation for D = 11 supergravity.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗