Search NASA⌕ Search

SEARCH · Search NASA

Results for “Higher order calculations”

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 217 records · Page 12

Exciton energy transfer reveals spectral signatures of excited states in clusters

Electronic excitation and concomitant energy transfer leading to Penning ionization in argon–acetylene clusters generated in a supersonic expansion are investigated with synchrotron-based photoionization mass spectrometry and electronic structure calculations. Spectral features in the photoionization efficiency of the mixed argon–acetylene clusters reveal a blue shift from the 2 P 1/2 and 2 P 3/2 excited states of atomic argon. Analysis of this feature suggests that excited states of argon clusters transfer energy to acetylene, resulting in its ionization and successive evaporation of argon. Theoretically calculated Ar n (n = 2–6) cluster spectra are in excellent agreement with experimental observations, and provide insight into the structure and ionization dynamics of the clusters. A comparison between argon–acetylene and argon–water clusters reveals that argon solvates water better, allowing for higher-order excitons and Rydberg states to be populated. Furthermore, these results are explained by theoretical calculations of respective binding energies and structures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Increasing Accuracy in Computed Inviscid Boundary Conditions

A technique has been devised to increase the accuracy of computational simulations of flows of inviscid fluids by increasing the accuracy with which surface boundary conditions are represented. This technique is expected to be especially beneficial for computational aeroacoustics, wherein it enables proper accounting, not only for acoustic waves, but also for vorticity and entropy waves, at surfaces. Heretofore, inviscid nonlinear surface boundary conditions have been limited to third-order accuracy in time for stationary surfaces and to first-order accuracy in time for moving surfaces. For steady-state calculations, it may be possible to achieve higher accuracy in space, but high accuracy in time is needed for efficient simulation of multiscale unsteady flow phenomena. The present technique is the first surface treatment that provides the needed high accuracy through proper accounting of higher-order time derivatives. The present technique is founded on a method known in art as the Hermitian modified solution approximation (MESA) scheme. This is because high time accuracy at a surface depends upon, among other things, correction of the spatial cross-derivatives of flow variables, and many of these cross-derivatives are included explicitly on the computational grid in the MESA scheme. (Alternatively, a related method other than the MESA scheme could be used, as long as the method involves consistent application of the effects of the cross-derivatives.) While the mathematical derivation of the present technique is too lengthy and complex to fit within the space available for this article, the technique itself can be characterized in relatively simple terms: The technique involves correction of surface-normal spatial pressure derivatives at a boundary surface to satisfy the governing equations and the boundary conditions and thereby achieve arbitrarily high orders of time accuracy in special cases. The boundary conditions can now include a potentially infinite number of time derivatives of surface-normal velocity (consistent with no flow through the boundary) up to arbitrarily high order. The corrections for the first-order spatial derivatives of pressure are calculated by use of the first-order time derivative velocity. The corrected first-order spatial derivatives are used to calculate the second- order time derivatives of velocity, which, in turn, are used to calculate the corrections for the second-order pressure derivatives. The process as described is repeated, progressing through increasing orders of derivatives, until the desired accuracy is attained.

Dyson, Roger↗

Calculation of transonic flow in a linear cascade

Turbomachinery blade designs are becoming more aggressive in order to achieve higher loading and greater range. New analysis tools are required to cope with these heavily loaded blades that may operate with a thin separated region near the trailing edge on the suction surface. An existing, viscous airfoil code was adapted to cascade conditions in an attempt to provide this capability. Comparisons with recently obtained data show that calculated and experimental surface Mach numbers were in good agreement but loss coefficients and outlet air angles were not.

Donovan, L. F.↗

Calculation of transonic flow in a linear cascade

Turbomachinery blade designs are becoming more aggressive in order to achieve higher loading and greater range. New analysis tools are required to cope with these heavily loaded blades that may operate with a thin separated region near the trailing edge on the suction surface. An existing, viscous airfoil code was adapted to cascade conditions in an attempt to provide this capability. Comparisons with recently obtained data show that calculated and experimental surface Mach numbers were in good agreement but loss coefficients and outlet air angles were not. Previously announced in STAR as N84-24539

Donovan, L. F.↗

Rotational relaxation in parahydrogen and its mixtures with helium, neon, and argon at 300 deg K.

Ultrasonic‐velocity dispersion measurements have been performed in parahydrogen and its mixtures with helium, neon, and argon, all at 300°K. In each case the experimental dispersion curves can be matched successfully to those calculated under the assumption that the 0–2 rotational transition relaxes separately from the 2–4 and higher‐order terms. For pure pH 2 we find a relaxation time 𝜏 20 of 1.30 × 10 −8 sec for the 2 → 0 transition and a 𝜏 42 of 3.90 × 10 −8 sec for the 4 → 2 transition. Comparison with the quantum‐mechanical theories of Roberts and of Davison for H 2 –H 2 collisions using Morse potentials shows good agreement for 𝜏 20 over the temperature range of 75°–300°K. The Morse‐potential asymmetry parameter yielding the best fit is β = 0.113 for Roberts' calculation and 0.108 for Davison's. It is found that He–pH 2 collisions are more effective than pH 2 –pH 2 in producing the J = 0 to J = 2 transition, but less effective for higher‐order transitions. Collisions of neon with pH 2 are found to be more effective at room temperature for inducing the 0 → 2 and 2 → 4 transitions than either helium or argon.

Leonard M Valley↗

Field redefinitions can be nonlocal

We revisit the lore establishing the allowed space of field redefinitions and show that there are essentially no restrictions. Our conclusions hold to all orders in perturbation theory and for any dispersion relation. Field redefinitions can be nonlocal, symmetry breaking, or in certain cases have explicit dependence on spacetime. We address field redefinitions that can be resummed into the propagator, which demonstrates how to perform perturbative calculations away from the minimum in field space. Field redefinitions are used to derive higher-order Schwinger-Dyson equations, which imply multiparticle soft theorems. Non-standard field redefinitions are showcased using both relativistic and nonrelativistic examples.

effective field theories↗

Measurements of Lund subjet multiplicities in 13 TeV proton-proton collisions with the ATLAS detector

This Letter presents a differential cross-section measurement of Lund subjet multiplicities, suitable for testing current and future parton shower Monte Carlo algorithms. This measurement is made in dijet events in 140 fb -1 of $\sqrt{s}$ =13 TeV proton–proton collision data collected with the ATLAS detector at CERN's Large Hadron Collider. The data are unfolded to account for acceptance and detector-related effects, and are then compared with several Monte Carlo models and to recent resummed analytical calculations. The experimental precision achieved in the measurement allows tests of higher-order effects in QCD predictions. Most predictions fail to accurately describe the measured data, particularly at large values of jet transverse momentum accessible at the Large Hadron Collider, indicating the measurement's utility as an input to future parton shower developments and other studies probing fundamental properties of QCD and the production of hadronic final states up to the TeV-scale.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Further Validation of a CFD Code for Calculating the Performance of Two-Stage Light Gas Guns

Earlier validations of a higher-order Godunov code for modeling the performance of two-stage light gas guns are reviewed. These validation comparisons were made between code predictions and experimental data from the NASA Ames 1.5" and 0.28" guns and covered muzzle velocities of 6.5 to 7.2 km/s. In the present report, five more series of code validation comparisons involving experimental data from the Ames 0.22" (1.28" pump tube diameter), 0.28", 0.50", 1.00" and 1.50" guns are presented. The total muzzle velocity range of the validation data presented herein is 3 to 11.3 km/s. The agreement between the experimental data and CFD results is judged to be very good. Muzzle velocities were predicted within 0.35 km/s for 74% of the cases studied with maximum differences being 0.5 km/s and for 4 out of 50 cases, 0.5 - 0.7 km/s.

Ballistic Range↗

Radiative Susceptibility of Cloudy Atmospheres to Droplet Number Perturbations: 1. Theoretical Analysis and Examples from MODIS

Theoretical and satellite-based assessments of the sensitivity of broadband shortwave radiative fluxes in cloudy atmospheres to small perturbations in the cloud droplet number concentration (N) of liquid water clouds under constant water conditions are performed. Two approaches to study this sensitivity are adopted: absolute increases in N, for which the radiative response is referred to as absolute cloud susceptibility, and relative increases in N or relative cloud susceptibility. Estimating the former is more challenging as it requires an assumed value for either cloud liquid water content or geometrical thickness; both susceptibilities require an assumed relationship between the droplet volume and effective radius. Expanding upon previous susceptibility studies, present radiative calculations include the effect of AN perturbations on droplet asymmetry parameter and single-scattering albedo, in addition to extinction. Absolute cloud susceptibility has a strong nonlinear dependence on the droplet effective radius as expected, while relative cloud susceptibility is primarily dependent on optical thickness. Molecular absorption and reflecting surfaces both reduce the relative contribution of the cloud to the top-of-atmosphere (TOA) flux and therefore also reduce the TOA albedo susceptibility. Transmittance susceptibilities are negative with absolute values similar to albedo susceptibility, while atmospheric absorptance susceptibilities are about an order of magnitude smaller than albedo susceptibilities and can be either positive or negative. Observation-based susceptibility calculations are derived from MODIS pixel-level retrievals of liquid water cloud optical thickness, effective radius, and cloud top temperature; two data granule examples are shown. Susceptibility quantifies the aerosol indirect effect sensitivity in a way that can be easily computed from model fields. As such, susceptibilities derived from MODIS observations provide a higher-order test of model cloud properties used for indirect effect studies. MODIS-derived global distributions of cloud susceptibility and radiative forcing calculations are presented in a companion paper.

Platnick, Steven↗

Towards Arbitrary Accuracy Inviscid Surface Boundary Conditions

Inviscid nonlinear surface boundary conditions are currently limited to third order accuracy in time for non-moving surfaces and actually reduce to first order in time when the surfaces move. For steady-state calculations it may be possible to achieve higher accuracy in space, but high accuracy in time is required for efficient simulation of multiscale unsteady phenomena. A surprisingly simple technique is shown here that can be used to correct the normal pressure derivatives of the flow at a surface on a Cartesian grid so that arbitrarily high order time accuracy is achieved in idealized cases. This work demonstrates that nonlinear high order time accuracy at a solid surface is possible and desirable, but it also shows that the current practice of only correcting the pressure is inadequate.

Dyson, Rodger W.↗

Prospects for strong coupling measurement at hadron colliders using soft-drop jet mass

We compute the soft-drop jet-mass distribution from pp collisions to NNLL accuracy while including nonperturbative corrections through a field-theory based formalism. Using these calculations, we assess the theoretical uncertainties on an αs precision measurement due to higher order perturbative effects, nonperturbative corrections, and PDF uncertainty. We identify which soft-drop parameters are well-suited for measuring αs, and find that higher-logarithmic resummation has a qualitatively important effect on the shape of the jet-mass distribution. We find that quark jets and gluon jets have similar sensitivity to αs, and emphasize that experimentally distinguishing quark and gluon jets is not required for an αs measurement. We conclude that measuring αs to the 10% level is feasible now, and with improvements in theory a 5% level measurement is possible. Getting down to the 1% level to be competitive with other state-of-the-art measurements will be challenging.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

The influence of thermochemistry on the reactivity of propane, the pentane isomers and $n$-heptane in the low temperature regime

In this study, the influence of thermochemistry on the reactivity of fuels at low temperatures (600–1000 K) is studied here. Specifically, the effect of different sets of thermochemistry on chemical model predictions is explored, where various sets are calculated at different levels of theory in addition to recently updated group additivity values. Experimentally measured ignition delay times for propane, the pentane isomers and n-heptane are simulated using NUIGMech1.2 and replacing the thermochemistry of the low-temperature species with the calculated values. For propane, three different thermochemistry sets were calculated, namely CCSD(T)-F12/TZ-F12//B2PLYP-D3/TZ//B2PLYP-D3/TZ (QM1), CCSD(T)-F12/TZ-F12//B2PLYP-D3/TZ//ωB97X-D/TZ (QM2) and B2PLYP-D3/TZ/ωB97X-D/6-31G*//ωB97X-D/6-31G* (QM3). The QM2 results provide parameters to optimize new group additivity (NGA) values which are used to calculate the fourth set of thermochemistry. The model predictions using these four sets are compared to those using NUIGMech1.2 for propane. As the QM1 and QM2 calculations are expensive, the thermochemistry calculated from the QM3 and NGA calculations are used in the pentane isomer and n-heptane models. For all of the models, it is found that the thermochemistry of the species involved in the low-temperature reaction sequence ($RH$, $\dot{R}$, $RO_2H$, $R\dot{O}_2$, $\dot{Q}OOH$ and $\dot{O}_2QOOH$ species) significantly affect fuel reactivity. The NGA values were developed based on all of these species except $\dot{O}_2QOOH$ radicals. The thermochemistry of $\dot{O}_2QOOH$ species cannot be accurately calculated with the NGA representations due to the importance of non-next-nearest neighbor interactions of –OOH substitution. Further development of the NGA method to capture such interactions is in progress. Overall, the model developed using the NGA thermochemistry shows better agreement with experimental data than the model using thermochemistry from affordable and prominent QM methods, such as QM3. Based on the results presented for propane, the pentane isomers and n-heptane, the thermochemistry calculated using the NGA method can be used to model the oxidation of higher order hydrocarbons at low temperatures.

33 ADVANCED PROPULSION SYSTEMS↗

K-epsilon turbulence model assessment with reduced numerical diffusion for coaxial jets

A numerical study has been conducted to evaluate the performance of the k-epsilon turbulence model for axisymmetric, unconfined, swirling, and nonswirling coaxial jet flows. Two lower order schemes, hybrid and power-law, and two higher order schemes, flux-spline and bounded skew upwind differencing, were employed in this investigation. The predicted results indicate that the higher order numerical schemes have greater potential for future model improvements and complex flow calculations in terms of storage, accuracy, and execution time. For the nonswirling flow, computations using an algebraic stress model have also been presented.

Nikjooy, M.↗

Aerosol characterization in smoke plumes from a wetlands fire

Results are presented from airborne measurements of aerosol mass loading, size distribution, and elemental composition obtained in a smoke plume from a fire at a Florida wildlife refuge. It is found that there was a high concentration of carbon-containing aerosols and salt crystals in the 0.1-0.2 micron size range; the composition and morphology of the aerosols differed with size. Aerosol mass concentrations are used in conjunction with those obtained for CO2 in order to calculate ratios of aerosol and CO2. The ratios are noted to be higher for the smoldering phase of the fire than for its flaming phase.

Woods, David C.↗

Utilization of Low Gravity Environment for Measuring Liquid Viscosity

The method of drop coalescence is used for determining the viscosity of highly viscous undercooled liquids. Low gravity environment is necessary in order to allow for examining large volumes affording much higher accuracy for the viscosity calculations than possible for smaller volumes available under 1 - g conditions. The drop coalescence method is preferred over the drop oscillation technique since the latter method can only be applied for liquids with vanishingly small viscosities. The technique developed relies on both the highly accurate solution of the Navier-Stokes equations as well as on data from experiments conducted in near zero gravity environment. Results are presented for method validation experiments recently performed on board the NASA/KC-135 aircraft. While the numerical solution was produced using the Boundary Element Method. In these tests the viscosity of a highly viscous liquid, glycerine at room temperature, was determined using the liquid coalescence method. The results from these experiments will be discussed.

Antar, Basil N.↗

Reducing negative weights in Monte Carlo event generation with Sherpa

An increase in theoretical precision of Monte Carlo event generators is typically accompanied by an increased need for computational resources. One major obstacle are negative weighted events, which appear in Monte Carlo simulations with higher perturbative accuracy. While they can be handled somewhat easily in fixed-order calculations, they are a major concern for particle level event simulations. In this article, the origin of negative weights in the S-MC@NLO method is reviewed and mechanisms to reduce the negative weight fraction in simulations with the Sherpa event generator are presented, with a focus on V+jets and tt+jets simulations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Measurement of the Sixth-Order Cumulant of Net-Proton Multiplicity Distributions in Au + Au Collisions at s NN = 27 , 54.4, and 200 GeV at RHIC

According to first principle Lattice QCD calculations, the transition from quark-gluon plasma to hadronic matter is a smooth crossover in the region µB ≤ Tc. As a result, higher-order cumulants and their ratios are predicted to be negative, C 6 /C 2 < 0, for example. In this paper, we report the first measurement of the midrapidity net-proton C 6 /C 2 from 27, 54.4 and 200 GeV Au+Au collisions at RHIC. The dependence on collision centrality and kinematic acceptance in (p T, y ) are analyzed. While for 27 and 54.4 GeV collisions the C 6 /C 2 values are close to zero within uncertainties, it is observed that for 200 GeV collisions, the C 6 /C 2 ratio becomes progressively negative from peripheral to central collisions. Transport model calculations without critical dynamics predict values around zero. These observations seem to favor a smooth crossover in the high energy nuclear collisions at RHIC.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Understanding Lithium Local Environments in LiMn 0.5 Ni 0.5 O 2 Cathodes: A DFT-Supported 6 Li Solid-State NMR Study

LiMn 0.5 Ni 0.5 O 2 has been plagued by the structural failure driven by the irreversible phase change upon subsequent cycling and the structural degradation caused by the seemingly unavoidable Li/Ni exchange. This begs for a better understanding of the correlation between lithium local structures, synthesis conditions, and overall electrochemical properties of materials. In this study, we use density functional theory (DFT) to assist experimental 6 Li solid-state magic angle spinning (MAS) nuclear magnetic resonance (NMR) spectroscopy to understand the nature and change in lithium local structures with extended annealing times (168 vs 15 h) synthesized at 900 °C. To correlate the NMR peak changes after different synthesis conditions, we calculated the NMR spectra of representative low-energy model configurations of LiMn 0.5 Ni 0.5 O 2 composition–zigzag, row, and flower structures–as well as the influence of Ni/Mn ordering and Li/Ni exchange on each spectrum using a supercell that is large enough to allow us to consider the effect of different relative configurations of adjacent transition metal (TM) layers. The analysis based on a combination of the calculated NMR and the experimental spectra suggests that most configurations contributing to the experimental spectra are in fact composed of a blend of configurations higher in energy than the classical well-ordered ground-state structures. Further, the extended annealing of LiMn 0.5 Ni 0.5 O 2 slightly enhances Ni/Mn and Li/Mn orderings while reducing the Li/Ni mixing ratio. This study shows how DFT calculations are crucial in providing a better understanding of the experimental characterization data and therefore local environments and domain structures in oxide materials, such as Li-ion battery cathodes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗