Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Calculation”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Lattice QCD calculation of 𝑥-dependent meson distribution amplitudes at physical pion mass with threshold logarithm resummation

We present a lattice quantum chromodynamics (QCD) calculation of the 𝑥-dependent pion and kaon distribution amplitudes (DA) in the framework of large momentum effective theory. This calculation is performed on a fine lattice of 𝑎 = 0.076 fm at physical pion mass, with the pion boosted to 1.8 GeV and kaon boosted to 2.3 GeV. We renormalize the matrix elements in the hybrid scheme and match to Math output error with a subtraction of the leading renormalon in the Wilson-line mass. The perturbative matching is improved by resumming the large logarithms related to the small quark and gluon momenta in the soft-gluon limit. After resummation, we demonstrate that we are able to calculate a range of 𝑥 ∈[𝑥 0 ,1 − 𝑥 0 ] with 𝑥 0 = 0.25 for pion and 𝑥 0 = 0.2 for kaon with theoretical systematic errors under control. The kaon DA is shown to be slighted skewed, and narrower than pion DA. Although the 𝑥-dependence cannot be direct calculated beyond these ranges, we estimate higher moments of the pion and kaon DAs by complementing our calculation with short-distance factorization.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Isomeric yield ratios of fission products: A missing piece in reactor antineutrino summation calculations

The calculation of the spectrum of antineutrinos ($\bar{v}_e$) from a reactor is a complicated problem requiring several nuclear data and physics inputs. Many of these have been investigated thoroughly to improve calculations and properly account for uncertainties. One input which has heretofore escaped consideration is the fission-yield distribution between ground and isomeric states. Here, in this work, we explore the effect of incorporating newly evaluated isomeric yield ratios (IYR) for 43 fission products into summation calculations and identify the disproportionate importance of certain isotopes, particularly at higher energies. Our analysis indicates that updated IYRs contribute to a significant increase in the $\bar{v}_e$ spectrum around and above 7 MeV, with increases of more than 50% at higher energies. Through a detailed sensitivity study, we highlight a number of isotopes for which the IYR has a substantial effect on the $\bar{v}_e$ spectrum. This work stresses the critical role of isomeric yields in calculations of reactor $\bar{v}_e$ spectra and points to the necessity for their accurate experimental determination, especially for key fission products, in order to refine our understanding and address the observed anomalies between measured and calculated $\bar{v}_e$ spectra.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Adapting CLUTCH methodology to multigroup TSUNAMI-3D for eigenvalue sensitivity calculations

The sensitivity of the eigenvalue to uncertainties in nuclear data and its evaluation are important for nuclear criticality safety. TSUNAMI-3D sequences within the SCALE code system offer several options to the user community for calculating eigenvalue sensitivity coefficients with multigroup (MG) and continuous energy (CE) 3D transport capabilities. TSUNAMI-3D sequences implement the adjoint-based perturbation theory with MG KENO code, the Contributon Linked eigenvalue sensitivity/Uncertainty estimation via Track length importance CHaracterization (CLUTCH) method with CE KENO code, and the Iterated Fission Probability (IFP) method with CE KENO and Shift codes. Each method has benefits and limitations depending on the problem that is run. The work presented here aims to adapt the CLUTCH method, which enables the Contributon method's mesh-free, memory-efficient approach for calculating adjoint-weighted tallies for sensitivity calculations, to the MG TSUNAMI-3D sequence. This application would eliminate the explicit adjoint KENO calculation, as well as the memory-consuming mesh flux moment tallies required by the conventional MG TSUNAMI-3D. Smaller memory footprints in the CLUTCH methodology and relatively shorter runtimes in MG KENO transport can make MG TSUNAMI-3D a viable method for some complex problems. Moreover, this adaptation allows MG sensitivity calculations with Shift, ORNL's next-generation high-performance Monte Carlo transport code, which currently does not offer any sensitivity capabilities with MG particle transport simulations. Initial implementation of the new MG TSUNAMI-3D sequence and its preliminary results with a selected critical benchmark experiment in the Verified, Archived Library of Inputs and Data (VALID) are presented in this study.

KENO↗
DOE OSTI2025

Isomeric Yield Ratios of fission products: a missing piece in reactor antineutrino summation calculations

The calculation of the spectrum of antineutrinos ( $\overline{V}$ e ) from a reactor is a complicated problem requiring several nuclear data and physics inputs. Many of these have been investigated thoroughly to improve calculations and properly account for uncertainties. One input which has heretofore escaped consideration is the fission yield distribution between ground and isomeric states. In this work, we explore the effect of incorporating newly evaluated isomeric yield ratios (IYR) for 43 fission products into summation calculations and identify the disproportionate importance of certain isotopes, particularly at higher energies. Our analysis indicates that updated IYRs contribute to a significant increase in the $\overline{V}$ e spectrum around and above 7 MeV, with increases of more than 50% at higher energies. Through a detailed sensitivity study, we highlight a number of isotopes for which the IYR has a substantial effect on the $\overline{V}$ e spectrum. This work stresses the critical role of isomeric yields in calculations of reactor $\overline{V}$ e spectra and points to the necessity for their accurate experimental determination, especially for key fission products, in order to refine our understanding and address the observed anomalies between measured and calculated $\overline{V}$ e spectra.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

A Novel Approach to Fast Calculation of High-Order Q-Cumulants

The method of Q-cumulants has been shown as a powerful tool to study the fine details of the azimuthal anisotropies in high-energy nucleus-nucleus collisions. A new method for the fast calculation of arbitrary order Q-cumulant $v_{n}\{2k\}$ values, based on the partition of a non-negative integer l $\le$ m for calculation of the 2m-particle azimuthal correlations is presented in this paper. Unlike the standard Q-cumulants method in which the calculation of high-order multi-particle calculations is impractical, the newly proposed method enables easy calculation. The validity of the method is proven via a toy model that uses the elliptic power distribution to simulate anisotropic emission of particles. The method enables the study of fine details of the $v_{2}$ distribution, such as higher-order central moments of the $v_{2}$ distribution, as well as the hydrodynamic behavior of the Quark-Gluon Plasma.

FOS: Physical sciences↗
NASA NTRS1974

Analysis and calculation of lightning-induced voltages in aircraft electrical circuits

Techniques to calculate the transfer functions relating lightning-induced voltages in aircraft electrical circuits to aircraft physical characteristics and lightning current parameters are discussed. The analytical work was carried out concurrently with an experimental program of measurements of lightning-induced voltages in the electrical circuits of an F89-J aircraft. A computer program, ETCAL, developed earlier to calculate resistive and inductive transfer functions is refined to account for skin effect, providing results more valid over a wider range of lightning waveshapes than formerly possible. A computer program, WING, is derived to calculate the resistive and inductive transfer functions between a basic aircraft wing and a circuit conductor inside it. Good agreement is obtained between transfer inductances calculated by WING and those reduced from measured data by ETCAL. This computer program shows promise of expansion to permit eventual calculation of potential lightning-induced voltages in electrical circuits of complete aircraft in the design stage.

Plumer, J. A.↗
NASA NTRS1973

Polarizability calculations on water, hydrogen, oxygen, and carbon dioxide

A semiclassical model of damped oscillators is used as a basis for the calculation of the dispersion of the refractive index, polarizability, and dielectric permeability in water, hydrogen, and oxygen in liquid and gaseous states, and in gaseous carbon dioxide. The absorption coefficient and the imaginary part of the refractive index are also calculated at corresponding wavelengths. A good agreement is obtained between the observed and calculated values of refractive indices, and between those of absorption coefficients in the region of absorption bands. The calculated values of oscillator strengths and damping factors are also discussed. The value of the polarizability of liquid water was about 2.8 times that of previous calculations.

Nir, S.↗
NASA NTRS1976

A submerged singularity method for calculating potential flow velocities at arbitrary near-field points

A discrete singularity method has been developed for calculating the potential flow around two-dimensional airfoils. The objective was to calculate velocities at any arbitrary point in the flow field, including points that approach the airfoil surface. That objective was achieved and is demonstrated here on a Joukowski airfoil. The method used combined vortices and sources ''submerged'' a small distance below the airfoil surface and incorporated a near-field subvortex technique developed earlier. When a velocity calculation point approached the airfoil surface, the number of discrete singularities effectively increased (but only locally) to keep the point just outside the error region of the submerged singularity discretization. The method could be extended to three dimensions, and should improve nonlinear methods, which calculate interference effects between multiple wings, and which include the effects of force-free trailing vortex sheets. The capability demonstrated here would extend the scope of such calculations to allow the close approach of wings and vortex sheets (or vortices).

Maskew, B.↗
NASA NTRS1977

Method of fan sound mode structure determination computer program user's manual: Modal calculation program

A computer user's manual describing the operation and the essential features of the Modal Calculation Program is presented. The modal Calculation Program calculates the amplitude and phase of modal structures by means of acoustic pressure measurements obtained from microphones placed at selected locations within the fan inlet duct. In addition, the Modal Calculation Program also calculates the first-order errors in the modal coefficients that are due to tolerances in microphone location coordinates and inaccuracies in the acoustic pressure measurements.

Pickett, G. F.↗
NASA NTRS1978

H-H, C-H, and C-C NMR spin-spin coupling constants calculated by the FP-INDO method for aromatic hydrocarbons

The FP-INDO (finite perturbation-intermediate neglect of differential overlap) method is used to calculate the H-H, C-H, and C-C coupling constants in hertz for molecules of six different benzenoid hydrocarbons: benzene, naphthalene, biphenyl, anthracene, phenanthrene, and pyrene. The calculations are based on both the actual and the average molecular geometries. It is found that only the actual molecular geometries can always yield the correct relative order of values for the H-H coupling constants. For the calculated C-C coupling constants, as for the calculated C-H coupling constants, the signs are positive (negative) for an odd (even) number of bonds connecting the two nuclei. Agreements between the calculated and experimental values of the coupling constants for all six molecules are comparable to those reported previously for other molecules.

Long, S. A. T.↗
NASA NTRS1980

Numerical calculations of the collapse of nonrotating, magnetic gas clouds

Results of the first self-consistent numerical calculations of the dynamic collapse of a magnetized protostellar gas cloud are presented. Symmetry about an axis parallel to the initial magnetic field direction has been assumed, so that the calculations could be performed on a two-dimensional grid. Also, the cloud was taken to be nonrotating and isothermal, and the magnetic field was assumed to remain frozen in to the gas. As starting models for the calculations, gas spheres with uniform density and magnetic field were used. The time evolution of the clouds has been calculated for roughly two initial free-fall times, at which point the central density has increased by a factor of approximately 10,000 to 1,000,000. Several such calculations have been performed for different values of the cloud's initial thermal, magnetic, and gravitational energies. In virtually all cases it is found that, once a flattened core forms in the cloud, the central magnetic field strength, B, varies with gas density, rho, according to (d log B/d log rho) = 1/2. This behavior is independent of the initial energy ratios mentioned above. It is also found that the magnetic field is able to prevent completely the collapse of part of the outer envelope of the cloud.

Scott, E. H.↗
NASA NTRS1980

Quasi-projection-operator calculation of autoionization states of O5/+/ /O VI/ and Si11/+/ /Si XII/

A practical approach to calculating diagnostic parameters includes both scattering and statelike (autoionization) elements. The paper presents inelastic quasi-projection-operator calculations of autoionization states of the three-electron ions of oxygen and silicon for states below the n = 3 manifod of the respective two-electron parent (O VI and Si XII). Wave functions containing up to 40 configurations are calculated. The O VI calculated results agree with experimental results within the experimental error (about + or - 2 eV) while for the Si XII case, it was necessary to include relativistic effects. Where formulas are available (states for which p + q = 2), the Si XII results are also in good agreement with experimental results for 2Se, 2Po, and 2De states. Since the calculated results for some of the higher states did not agree with identifications proposed by Trabert et al. (1979), alternative identifications are suggested.

Wakid, S.↗
NASA NTRS1982

Improved Venus ionopause altitude calculation and comparison with measurement

Improved approximations are incorporated into the inviscid fluid method originated by Spreiter et al (1970) in calculations of the altitude of the Venus ionopause. The altitude calculations are then compared with median altitudes measured by the Pioneer Venus retarding potential analyzer. The calculated ionopause shape is found to closely approximate the measured meridian shape in the solar zenith angle (SZA) range of zero-135 deg. The use of improved ionospheric pressure field ionosheath pressure and Prandtl-Meyer expansion approximations lowers the terminator ionopause altitude to approximately half that obtained with the usual Spreiter approximations. It is also determined that the calculated dawn ionopause is about 300 km higher than the dusk ionopause, and that both ionopauses are close to their respective measured meridian values. It is concluded that median ionopause altitude within the SZA angle range may be calculated without inclusion of a viscous interaction in the theory.

Knudsen, W. C.↗
NASA NTRS1982

Numerical calculation of the parameters of the efflux from a helium dewar used for cooling of heat shields in a satellite

The parameters of the efflux from a helium dewar in space were numerically calculated. The flow was modeled as a one dimensional compressible ideal gas with variable properties. The primary boundary conditions are flow with friction and flow with heat transfer and friction. Two PASCAL programs were developed to calculate the efflux parameters: EFFLUZD and EFFLUXM. EFFLUXD calculates the minimum mass flow for the given shield temperatures and shield heat inputs. It then calculates the pipe lengths, diameter, and fluid parameters which satisfy all boundary conditions. Since the diameter returned by EFFLUXD is only rarely of nominal size, EFFLUXM calculates the mass flow and shield heat exchange for given pipe lengths, diameter, and shield temperatures.

Brendley, K.↗
NASA NTRS1982

Analytical calculation of a single jet in crossflow and comparison with experiment

A series of calculations of a jet in crossflow using a three-dimensional finite difference model of the Navier-Stokes equations were made on a series of progressively finer grids. With a reasonable number of grid point (40x30x20), calculated jet penetration and mixing characteristics were found to compare favorably with the experimental measurements of Crabb, Durao and Whitelaw. The calculated results were grid dependent (i.e., numerical or false diffusion was present in the result). A finer grid calculation indicated that turbulence model deficiencies may become more noticeable in the calculated results as the magnitude of numerical diffusion is reduced.

Clauss, R. W.↗
NASA NTRS1985

Transonic calculations for a flexible supercritical wing and comparison with experiment

Pressure data measured on the flexible DAST ARW-2 wing are compared with results calculated using the transonic small perturbation code XTRAN3S. A brief description of the analysis is given and a recently-developed grid coordinate transformation is described. Calculations are presented for the rigid and flexible wing for Mach numbers from 0.60 to 0.90 and dynamic pressures from 0 to 1000 psf. Calculated and measured static pressures and wing deflections are compared, and calculated static aeroelastic trends are given. Attempts to calculate the transonic instability boundary of the wing are described.

Bennett, R. M.↗
NASA NTRS1985

Comparison of big event with calculations of the air shower development

The incidence of high energy hadrons and electron-photons in air showers at various stages of development is calculated. Numerical calculation is used to solve the diffusion equation for a nuclear cascade and analytical calculation for cascade shower induced gamma rays. From these calculations, one can get the longitudinal development of the high energy hadron and electron-photon components, and the energy spectra of these components at various depths of air shower development. The total number of hadrons (N sub H) and electron-photon components (N sub gamma) are related according to stages of the air shower development and primary energy. The relation of the total energy of hadron and electron-photon component above the threshold energy is given. The energy balance between both components is also a useful parameter to study high energy events accompanying air showers. The relation of N sub H and fractional hadronic energy E (sum E sub H sup gamma/sum E sub H sup gamma + Sum E sub gamma) is calculated. This relation is helpful to understand the stage of air shower development(t) and primary energy (E sub p).

Niwa, M.↗
NASA NTRS1985

Transonic calculations for a flexible supercritical wing and comparison with experiment

Pressure data measured on the flexible DAST ARW-2 wing are compared with results calculated using the transonic small perturbation code XTRAN3S. A brief description of the analysis is given and a recently developed grid coordinate transformation is described. Calculations are presented for the rigid and flexible wing for Mach numbers from 0.60 to 0.90 and dynamic pressures from 0 to 1000 psf. Calculated and measured static pressures and wing deflections are compared, and calculated static aeroelastic trends are given. Attempts to calculate the transonic instability boundary of the wing are described.

Bennett, R. M.↗