Search NASA⌕ Search

SEARCH · Search NASA

Results for “pair correlation function”

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.

281 records · Page 16
All sourcesNASA NTRSDOE OSTI
NASA NTRS2006

Direct Numerical Simulation of Transitional Multicomponent-Species Gaseous and Multicomponent-Liquid Drop-Laden Mixing

A model of multicomponent-liquid (MC-liquid) drop evaporation in a three-dimensional mixing layer is here exercised at larger Reynolds numbers than in a previous study, and transitional states are obtained. The gas phase is followed in an Eulerian frame and the multitude of drops is described in a Lagrangian frame. Complete coupling between phases is included with source terms in the gas conservation equations accounting for the drop/flow interaction in terms of drop drag, drop heating and species evaporation. The liquid composition, initially specified as a single-Gamma (SG) probability distribution function (PDF) depending on the molar mass is allowed to evolve into a linear combination of two SGPDFs, called the double-Gamma PDF (DGPDF). The compositions of liquid and vapor emanating from the drops are calculated through four moments of the DGPDFs, which are drop-specific and location-specific, respectively. The mixing layer is initially excited to promote the double pairing of its four initial spanwise vortices into an ultimate vortex in which small scales proliferate. Simulations are performed for four liquids of different compositions and the effect of the initial mass loading and initial free-stream gas temperature are explored. For reference, Simulations are also performed for gaseous multicomponent mixing layers for which the effect of Reynolds number is investigated. The results encompass examination of the global layer characteristics, flow visualizations and homogeneous-plane statistics at transition. Comparisons are performed with previous pre-transitional MC-liquid simulations and with transitional single-component (SC) liquid studies. It is found that MCC flows at transition, the classical energy cascade is of similar strength, but that the smallest scales contain orders of magnitude less energy than SC flows, which is confirmed by the larger viscous dissipation in the former case. Contrasting to pre-transitional MC flows, the vorticity and drop organization depend on the initial gas temperature, this being due to the drop/turbulence coupling. The vapor-composition mean molar mass and standard deviation distributions strongly correlate with the initial liquid-composition PDF; such a correlation only exists for the magnitude of the mean but not for that of the standard deviation. Unlike in pre-transitional situations, regions of large composition standard deviation no longer necessarily coincide with regions of large mean molar mass. The kinetic energy, rotational and composition characteristics, and dissipation are liquid specific and the variation among liquids is amplified with increasing free-stream gas temperature. Eulerian and Lagrangian statistics of gas-phase quantities show that the different. Observation framework may affect the perception of the flow characteristics. The gas composition, of which the first four moments are calculated, is shown to be close to, but distinct from a SGPDF. The PDF of the scalar dissipation rate is calculated for drop-laden layers and is shown to depart more significantly from the typically assumed Gaussian in gaseous flows than experimentally measured gaseous scalar dissipation rates, this being attributed to the increased heterogeneity due to drop/flow interactions.

tubulent multicompenent-species↗
DOE OSTI2025

Exploring the Nucleon Structure via Deep Electroproduction Processes

Understanding the internal structure of the nucleon is a fundamental goal of modern physics, which aims at a comprehensive framework describing the internal dynamics of quarks and gluons. Among other structure functions, Generalized Parton Distributions (GPDs) offer a powerful framework for describing the nucleon dynamics by correlating the longitudinal momentum and the transverse position of its internal partons. Such a correlation provides a three-dimensional picture of the nucleon and enables access to fundamental properties, including the internal pressure distributions and the parton’s angular momentum contribution to the nucleon's total spin, thereby playing a central role in resolving the nucleon spin puzzle. At the Thomas Jefferson National Facility (JLab), polarized electron beam experiments allow for probing GPDs through the measurement of hard exclusive processes. Among the cleanest experimental channels, we find the electro-production of a real photon through the Deeply Virtual Compton Scattering (DVCS) mechanism. The first data-taking period of the CLAS12 program, taking place in 2018, allowed for unique DVCS Beam Spin Asymmetry (BSA) measurements in the phase space covered by a 10.6 GeV polarized electron beam impinging on an unpolarized liquid hydrogen target. Although detecting all final-state particles ensures exclusivity of the process, conservation laws indicate that it is not mandatory. I adopt an approach omitting the direct detection of the recoil proton, providing a simplified yet effective event selection strategy that boosts statistics and gives access to a larger phase space sensitive to the underlying GPD dynamics through BSA and cross section measurements. The Double DVCS (DDVCS) process promises a dedicated mapping of GPDs. Contributing to the electro-production of a lepton pair cross-section, the DDVCS reaction extends DVCS by allowing the final-state photon to be virtual, enriching the kinematic phase space and providing unique access to the internal correlations encoded by GPDs. A feasibility study is conducted to assess the potential of future DDVCS measurements at Jefferson Lab and the future Electron-Ion Collider (EIC). While Jefferson Lab will provide DDVCS measurements in the valence region through the SoLID$\mu$ and $\mu$CLAS12 experimental projects, in the long term, the EIC will provide complementary measurements in the sea region, both accessing unprecedented information about GPDs in a phase space region otherwise inaccessible. Taken together, these investigations demonstrate both the current capabilities and future opportunities for probing GPDs through exclusive processes. The experimental analysis of DVCS at CLAS12 provides precise measurements within an established framework, while the phenomenological study of DDVCS opens the door to richer and more comprehensive explorations with future detectors and facilities.

Alvarado, Juan [Université Paris-Saclay: Gif-sur-Y↗
NASA NTRS2003

Capacitive Sensors for Measuring Masses of Cryogenic Fluids

An effort is under way to develop capacitive sensors for measuring the masses of cryogenic fluids in tanks. These sensors are intended to function in both microgravitational and normal gravitational settings, and should not be confused with level sensors, including capacitive ones. A sensor of this type is conceptually simple in the sense that (1) it includes only one capacitor and (2) if properly designed, its single capacitance reading should be readily convertible to a close approximation of the mass of the cryogenic fluid in the tank. Consider a pair of electrically insulated electrodes used as a simple capacitive sensor. In general, the capacitance is proportional to the permittivity of the dielectric medium (in this case, a cryogenic fluid) between the electrodes. The success of design and operation of a sensor of the present type depends on the accuracy of the assumption that to a close approximation, the permittivity of the cryogenic fluid varies linearly with the density of the fluid. Data on liquid nitrogen, liquid oxygen, and liquid hydrogen, reported by the National Institute of Standards and Technology, indicate that the permittivities and densities of these fluids are, indeed, linearly related to within a few tenths of a percent over the pressure and temperature regions of interest. Hence, ignoring geometric effects for the moment, the capacitance between two electrodes immersed in the fluid should vary linearly with the density, and, hence, with the mass of the fluid. Of course, it is necessary to take account of the tank geometry. Because most cryogenic tanks do not have uniform cross sections, the readings of level sensors, including capacitive ones, are not linearly correlated with the masses of fluids in the tanks. In a sensor of the present type, the capacitor electrodes are shaped so that at a given height, the capacitance per unit height is approximately proportional to the cross-sectional area of the tank in the horizontal plane at that height (see figure).

Nurge, Mark↗
NASA NTRS1992

Computational studies of metal-metal and metal-ligand interactions

The geometric structure of Cr(CO)6 is optimized at the modified coupled-pair functional (MCPF), single and double excitation coupled-cluster (CCSD) and CCSD(T) levels of theory (including a perturbational estimate for connected triple excitations), and the force constants for the totally symmetric representation are determined. The geometry of Cr(CO)5 is partially optimized at the MCPF, CCSD and CCSD(T) levels of theory. Comparison with experimental data shows that the CCSD(T) method gives the best results for the structures and force constants, and that remaining errors are probably due to deficiencies in the one-particle basis sets used for CO. A detailed comparison of the properties of free CO is therefore given, at both the MCPF and CCSD/CCSD(T) levels of treatment, using a variety of basis sets. With very large one-particle basis sets, the SSCD(T) method gives excellent results for the bond distance, dipole moment and harmonic frequency of free CO. The total binding energies of Cr(CO)6 and Cr(CO)5 are also determined at the MCPF, CCSD and CCSD(T) levels of theory. The CCSD(T) method gives a much larger total binding energy than either the MCPF or CCSD methods. An analysis of the basis set superposition error (BSSE) at the MCPF level of treatment points out limitations in the one-particle basis used here and in a previous study. Calculations using larger basis sets reduced the BSSE, but the total binding energy of Cr(CO)6 is still significantly smaller than the experimental value, although the first CO bond dissociation energy of Cr(CO)6 is well described. An investigation of 3s3p correlation reveals only a small effect. The remaining discrepancy between the experimental and theoretical total binding energy of Cr(CO)6 is probably due to limitations in the one-particle basis, rather than limitations in the correlation treatment. In particular an additional d function and an f function on each C and O are needed to obtain quantitative results. This is underscored by the fact that even using a very large primitive se (1042 primitive functions contracted to 300 basis functions), the superposition error for the total binding energy of Cr(CO)6 is 22 kcal/mol at the MCPF level of treatment.

Barnes, Leslie A.↗
DOE OSTI2024

Graph network heterogeneity predicts interplant wake losses

Wind plants generate large-scale wakes, which can affect the performance of neighboring installations. Such wakes are challenging to estimate due to the inherent complexity in modeling wake interactions between large quantities of turbines at various distances. Weighted directed graph networks can inform complex models by linking turbine pairs into chains of upstream and downstream neighbors for a given wind direction. A novel interpretation of the graph network adjacency matrix is proposed where each element of the matrix represents the cumulative impact of upstream turbines on an individual. In this study, wake losses were estimated with an engineering wake model across a range of inflow conditions for nine parametric variations of a system containing two neighboring wind plants. The parametric nature of the study isolates turbine spacing within the plant, separation distance between plants, and wind direction as the main drivers of wake losses. Spatial heterogeneity is computed from the weighted average adjacency matrix of each plant arrangement. The proposed method is orders of magnitude faster than wake modeling and does not require detailed turbine information or atmospheric conditions. Furthermore, the weighted average adjacency matrix provides insight on the spatial organization of wake losses at various scales. Plant heterogeneity is correlated with wake losses within and among plants. Framing wind plant wake interaction in terms of graph network spatial heterogeneity provides an efficient approach for predicting wake losses within and among neighboring wind plants with applications to other complex systems where wake interactions are key factor.

17 WIND ENERGY↗
NASA NTRS1998

Analysis and Testing of Plates with Piezoelectric Sensors and Actuators

Piezoelectric material inherently possesses coupling between electrostatics and structural dynamics. Utilizing linear piezoelectric theory results in an intrinsically coupled pair of piezoelectric constitutive equations. One equation describes the direct piezoelectric effect where strains produce an electric field and the other describes the converse effect where an applied electrical field produces strain. The purpose of this study is to compare finite element analysis and experiments of a thin plate with bonded piezoelectric material. Since an isotropic plate in combination with a thin piezoelectric layer constitutes a special case of a laminated composite, the classical laminated plate theory is used in the formulation to accommodated generic laminated composite panels with multiple bonded and embedded piezoelectric layers. Additionally, the von Karman large deflection plate theory is incorporated. The formulation results in laminate constitutive equations that are amiable to the inclusion of the piezoelectric constitutive equations yielding in a fully electro-mechanically coupled composite laminate. Using the finite element formulation, the governing differential equations of motion of a composite laminate with embedded piezoelectric layers are derived. The finite element model not only considers structural degrees of freedom (d.o.f.) but an additional electrical d.o.f. for each piezoelectric layer. Comparison between experiment and numerical prediction is performed by first treating the piezoelectric as a sensor and then again treating it as an actuator. To assess the piezoelectric layer as a sensor, various uniformly distributed pressure loads were simulated in the analysis and the corresponding generated voltages were calculated using both linear and nonlinear finite element analyses. Experiments were carried out by applying the same uniformly distributed loads and measuring the resulting generated voltages and corresponding maximum plate deflections. It is found that a highly nonlinear relationship exists between maximum deflection and voltage versus pressure loading. In order to assess comparisons of predicted and measured piezoelectric actuation, sinusoidal excitation voltages are simulated/applied and maximum deflections are calculated/measured. The maximum deflection as a function of time was determined using the linear finite elements analysis. Good correlation between prediction and measurement was achieved in all cases.

Bevan, Jeffrey S.↗
NASA NTRS1998

Cardiogenic oscillation phase relationships during single-breath tests performed in microgravity

We studied the phase relationships of the cardiogenic oscillations in the phase III portion of single-breath washouts (SBW) in normal gravity (1 G) and in sustained microgravity (microG). The SBW consisted of a vital capacity inspiration of 5% He-1.25% sulfurhexafluoride-balance O2, preceded at residual volume by a 150-ml Ar bolus. Pairs of gas signals, all of which still showed cardiogenic oscillations, were cross-correlated, and their phase difference was expressed as an angle. Phase relationships between inspired gases (e.g., He) and resident gas (n2) showed no change from 1 G (211 +/- 9 degrees) to microG (163 +/- 7 degrees). Ar bolus and He were unaltered between 1 G (173 +/- 15 degrees) and microG (211 +/- 25 degrees), showing that airway closure in microG remains in regions of high specific ventilation and suggesting that airway closure results from lung regions reaching low regional volume near residual volume. In contrast, CO2 reversed phase with He between 1 G (332 +/- 6 degrees) and microG (263 +/- 27 degrees), strongly suggesting that, in microG, areas of high ventilation are associated with high ventilation-perfusion ratio (VA/Q). This widening of the range of VA/Q in microG may explain previous measurements (G.K. Prisk, A.R. Elliott, H.J.B. Guy, J.M. Kosonen, and J.B. West J. Appl. Physiol. 79: 1290-1298, 1995) of an overall unaltered range of VA/Q in microG, despite more homogeneous distributions of both ventilation and perfusion.

Non-NASA Center↗
DOE OSTI2024

Entanglement entropy of a color flux tube in (2+1)D Yang-Mills theory

We construct a novel flux tube entanglement entropy (FTE 2 ), defined as the excess entanglement entropy relative to the vacuum of a region of color flux stretching between a heavy quark-anti-quark pair in pure gauge Yang-Mills theory. We show that FTE 2 can be expressed in terms of correlators of Polyakov loops, is manifestly gauge-invariant, and therefore free of the ambiguities in computations of the entanglement entropy in gauge theories related to the choice of the center algebra. Employing the replica trick, we compute FTE 2 for SU(2) Yang-Mills theory in (2+1)D and demonstrate that it is finite in the continuum limit. We explore the properties of FTE 2 for a half-slab geometry, which allows us to vary the width and location of the slab, and the extent to which the slab cross-cuts the color flux tube. Following the intuition provided by computations of FTE 2 in (1+1)D, and in a thin string model, we examine the extent to which our FTE 2 results can be interpreted as the sum of an internal color entropy and a vibrational entropy corresponding to the transverse excitations of the string.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Entanglement entropy of a color flux tube in (1+1)D Yang–Mills theory

In recent work Amorosso et al. (2024), we computed a novel flux tube entanglement entropy (FTE) of the color flux tube stretched between a heavy quark-antiquark pair on a Euclidean lattice in (2+1)D Yang–Mills theory. Our numerical results suggested that FTE can be partitioned into an internal color entanglement entropy and a vibrational entropy corresponding to the transverse excitations of a QCD string, with the latter described by a thin string model. Since the color flux tube does not have transverse excitations in (1+1)D Yang–Mills theory, we use this simpler framework to perform an exact analytical computation of the contribution of the internal color degrees of freedom to FTE. For the multipartite partitioning of the color flux tube, we find the remarkable result that FTE depends only on the dimension of the color group representation and the number of times the flux tube crosses the boundary between the traced and untraced spatial regions but not on the string length. Our proof is independent of whether the replica and region boundaries on the lattice are placed on vertices or in plaquette centers.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Dynamics and Structure of Molten Neodymium Titanate

Binary rare earth titanates can be vitrified using levitation melt quenching, and their glasses are promising for optical applications. To better understand their thermophysical properties and dynamics, we have used molecular dynamics (MD) and empirical potential structure refinement (EPSR) to simulate molten neodymium titanate, 83TiO 2 –17Nd 2 O 3 . The models are compared with high-energy synchrotron X-ray diffraction of the melt structure over 1060–1960 K and thermophysical property measurements of electrostatically levitated droplets over 1500–2200 K. As temperature increases, the network periodicity increases and the correlation length decreases, both in monotonic trends that suggest continuous structural evolution through the melt and supercooled liquid regimes. At 1900 K, the mean Ti–O and Nd–O coordination numbers are 5.8(5) and 6.7(5), compared to the glass, 5.72(6) and 7.70(26). The equilibrium melting point is measured as 1860(20) K, over 100 K higher than previous estimates. Viscosity data from MD are in good agreement with experimental measurements. The dynamic propensity of Nd is similar to that of the network-forming Ti, suggesting that Nd does not participate as a conventional modifier. The viscosity of this fragile liquid is likely controlled by the breaking of cation–oxygen bonds and rearrangements of Ti–O polyhedra.

molecular dynamics↗
DOE OSTI2025

Microstructure‐Dependent Sodium Storage Mechanisms in Hard Carbon Anodes

Sustainable energy storage is essential to support the transition to renewables and meet the increasing demand for energy. Sodium‐ion batteries (NIBs) are attractive for grid‐scale energy storage due to the abundance and low cost of sodium, sustainability of other battery components, and electrochemical performance. Hard carbon (HC) is a leading anode material for NIBs, but its complex microstructure complicates the understanding of sodium storage mechanisms. Using X‐ray total scattering and density functional theory calculations, this study clarifies how HC's microstructural variations influence sodium storage across the slope (high potential) and plateau (low potential) regions of the potential capacity curve. In the slope region, sodium initially adsorbs at high‐binding energy defect sites and subsequently intercalates between graphene layers, adsorbing at low‐binding energy defect sites, correlating with different slopes observed during initial sodiation. Initial irreversibility arises from sodium trapping at surface defects and solid electrolyte interface formation. In the plateau region, sodium simultaneously intercalates and fills pores, influenced by pore size, interlayer spacing, and defect concentration. HCs with larger pore sizes form larger sodium clusters. In conclusion, the proposed mechanism underscores the role of microstructure engineering in enhancing HC performance and advancing NIBs for grid‐scale energy storage.

36 MATERIALS SCIENCE↗