Search NASA⌕ Search

SEARCH · Search NASA

Results for “calculation method”

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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Application of a monotonic upstream-biased transport scheme to three-dimensional constituent transport calculations

The application of van Leer's scheme, a monotonic, upstream-biased differencing scheme, to three-dimensional constituent transport calculations is shown. The major disadvantage of the scheme is shown to be a self-limiting diffusion. A major advantage of the scheme is shown to be its ability to maintain constituent correlations. The scheme is adapted for a spherical coordinate system with a hybrid sigma-pressure coordinate in the vertical. Special consideration is given to cross-polar flow. The vertical wind calculation is shown to be extremely sensitive to the method of calculating the divergence. This sensitivity implies that a vertical wind formulation consistent with the transport scheme is essential for accurate transport calculations. The computational savings of the time-splitting method used to solve this equation are shown. Finally, the capabilities of this scheme are illustrated by an ozone transport and chemistry model simulation.

Allen, Dale J.↗
NASA NTRS2012

Calculation of Sensitivity Derivatives in an MDAO Framework

During gradient-based optimization of a system, it is necessary to generate the derivatives of each objective and constraint with respect to each design parameter. If the system is multidisciplinary, it may consist of a set of smaller "components" with some arbitrary data interconnection and process work ow. Analytical derivatives in these components can be used to improve the speed and accuracy of the derivative calculation over a purely numerical calculation; however, a multidisciplinary system may include both components for which derivatives are available and components for which they are not. Three methods to calculate the sensitivity of a mixed multidisciplinary system are presented: the finite difference method, where the derivatives are calculated numerically; the chain rule method, where the derivatives are successively cascaded along the system's network graph; and the analytic method, where the derivatives come from the solution of a linear system of equations. Some improvements to these methods, to accommodate mixed multidisciplinary systems, are also presented; in particular, a new method is introduced to allow existing derivatives to be used inside of finite difference. All three methods are implemented and demonstrated in the open-source MDAO framework OpenMDAO. It was found that there are advantages to each of them depending on the system being solved.

Moore, Kenneth T.↗
NASA NTRS1995

Calculating 3-D Crack-Opening Behavior Of A Fatigue Crack

Method for calculating through-the-thickness, three-dimensional (3-D) crack-opening behavior in standard fatigue-crack-growth test specimens devised. Involves determination of crack-opening loads from crack-propagation data and finite-element analyses. In method, one considers variation in rate of crack growth along crack front, relationship between rate of crack growth and effective stress intensity factor, and variation in three-dimensional stress-intensity factor along crack front. Developed to provide simpler means for obtaining fatigue-crack-opening behavior loads.

Newman, J. C., Jr.↗
DOE OSTI2021

A Statistics-Based Threshold for the RMS-Energy Oscillation Detector

System operators and reliability coordinators currently rely on extensive baselining studies to set thresholds for their oscillation detectors. The resulting thresholds are based largely on engineering judgment and may vary significantly between organizations. In this paper, statistical distribution theory is used to derive a detection threshold for the widely deployed root-mean-square (RMS)-energy detector. This expression provides a theoretical basis for the detector's configuration, simplifies the process of selecting the threshold, and enables improved consistency among organizations that need to coordinate during system-wide events. Three methods for calculating the threshold using synchrophasor measurements are also proposed. These methods ensure that the threshold can be calculated reliably for various applications. Tests with simulated and field-measured data demonstrate that the statistics-based threshold provides consistent detection of grid disturbances while maintaining a low probability of false alarm.

Follum, James D.↗
DOE OSTI2025

Impacts of Bulk Microphysics Scheme Structural Choices on Simulations of Rain Initiation Through Drop Coalescence

This study examines how different structural choices in bulk microphysics schemes impact the simulation of warm rain initiation. A single liquid category (SLC) approach prognosing up to four moments of a single drop size distribution (DSD) is compared to the traditional two-category, two-moment approach with separate DSDs for cloud and rain (four total prognostic variables). Different methods for calculating tendencies of the prognostic variables from drop collision-coalescence are also tested: a discretized numerical-integration approach, machine learning via neural networks, lookup tables, and traditional power law fits. Relative to simulations using a bin microphysics model, SLC gives smaller error overall than the two-category approach when numerical integration is used to calculate the collision-coalescence tendencies for both. Replacing the numerical integration with a pre-computed lookup table reduces computational cost with little loss of accuracy. However, using fitted power laws with SLC to represent the collision-coalescence tendencies substantially reduces accuracy and leads to an order of magnitude increase in error. It is also demonstrated that with SLC, reasonably accurate solutions are obtained using only three prognostic moments, while a two-moment SLC scheme leads to substantial error. Overall, both the choice of prognostic moments (e.g., SLC vs. two-category) and method to calculate the collision-coalescence tendencies are important to consider for minimizing errors in bulk schemes. SLC with a sufficiently detailed calculation of the collision-coalescence tendencies provides accurate solutions for a reasonable computational cost, providing a viable alternative to the traditional two-category, two-moment approach for bulk microphysics.

320 (cloud physics and chemistry)↗
NASA NTRS1975

Calculation of arbitrary-order diffraction efficiencies of thick gratings with arbitrary grating shape

A method for calculating arbitrary-order diffraction efficiencies of thick, lossless transmission gratings with arbitrary periodic grating shapes has been developed. A Fourier-series representation of the grating is employed, along with a coupled-mode theory of diffraction. For illustration, numerical values of the diffraction efficiencies at the first three Bragg angles are calculated for sinusoidal, square-wave, triangular, and saw-tooth gratings. Numerical results for the same grating shapes with the same parameters are also calculated for comparison, by extending Burckhardt's numerical method for analyzing thick sinusoidal gratings. The comparison shows that the coupled-mode theory provides results with relative computational ease and results that are in agreement with calculations obtained by extending the more-rigorous Burckhardt theory to nonsinusoidal grating shapes and to higher-order Bragg angles.

Su, S. F.↗
DOE OSTI2020

Prediction of tissue optical properties using the Monte Carlo modeling of photon transport in turbid media and integrating spheres

Monte Carlo methods are an established technique for simulating light transport in biological tissue. Integrating spheres make experimental measurements of the reflectance and transmittance of a sample straightforward and inexpensive. This work presents an extension to existing Monte Carlo photon transport methods to simulate integrating sphere experiments. Crosstalk between spheres in dual-sphere experiments is accounted for in the method. Analytical models, previous works on Monte Carlo photon transport, and experimental measurements of a synthetic tissue phantom validate this method. We present two approaches for using this method to back-calculate the optical properties of samples. Experimental and simulation uncertainties are propagated through both methods. Both back-calculation methods find the optical properties of a sample accurately and precisely. Our model is implemented in standard Python 3 and CUDA C++ [J. Nickolls, I. Buck, M. Garland, and K. Skadron, ACM Queue 6 , 40 ( 2008 ) ] and is publicly available in Code 1.

Cook, Patrick D. (ORCID:0000000279345428)↗
NASA NTRS1990

The C2 Ballik-Ramsay transitions in carbon stars

The bands of C2 are an important source of opacity in the near-infrared spectral region of carbon stars. In the 2 micron spectral region, the Delta v = -1 and -2 transitions of the Ballik-Ramsay b 3Sigma(g)(-) - a 3Pi(u) system are especially important. A general method for calculating highly accurate transition parameters of the Ballik-Ramsay band is presented. This method allows the easy calculation of the satellite lines and the forbidden lines resulting from the perturbation of the upper b 3Sigma(g)(-) state by the X 1Sigma(0)(+) state. Included in the calculation is the effect of the isotopic molecules (C-12)2, (C-12)(C-13), and (C-13)2. Three consistency tests of the calculated strengths are presented. A comparison of synthetic spectra is made with an observed spectrum of WZ Cas in the 3940/cm spectral region.

Goorvitch, David↗
NASA NTRS2010

Lightning Reporting at 45th Weather Squadron: Recent Improvements

The 45th Weather Squadron (45 WS) provides daily lightning reports to space launch customers at CCAFS/KSC. These reports are provided to assess the need to inspect the electronics of satellite payloads, space launch vehicles, and ground support equipment for induced current damage from nearby lightning strokes. The 45 WS has made several improvements to the lightning reports during 2008-2009. The 4DLSS, implemented in April 2008, provides all lightning strokes as opposed to just one stroke per flash as done by the previous system. The 45 WS discovered that the peak current was being truncated to the nearest kilo amp in the database used to generate the daily lightning reports, which led to an up to 4% underestimate in the peak current for average lightning. This error was corrected and led to elimination of this underestimate. The 45 WS and their mission partners developed lightning location error ellipses for 99% and 95% location accuracies tailored to each individual stroke and began providing them in the spring of 2009. The new procedure provides the distance from the point of interest to the best location of the stroke (the center of the error ellipse) and the distance to the closest edge of the ellipse. This information is now included in the lightning reports, along with the peak current of the stroke. The initial method of calculating the error ellipses could only be used during normal duty hours, i.e. not during nights, weekends, or holidays. This method was improved later to provide lightning reports in near real-time, 24/7. The calculation of the distance to the closest point on the ellipse was also significantly improved later. Other improvements were also implemented. A new method to calculate the probability of any nearby lightning stroke. being within any radius of any point of interest was developed and is being implemented. This may supersede the use of location error ellipses. The 45 WS is pursuing adding data from nine NLDN sensors into 4DLSS in real-time. This will overcome the problem of 4DLSS missing some of the strong local strokes. This will also improve the location accuracy, reduce the size and eccentricity of the location error ellipses, and reduce the probability of nearby strokes being inside the areas of interest when few of the 4DLSS sensors are used in the stroke solution. This will not reduce 4DLSS performance when most of the 4DLSS sensors are used in the stroke solution. Finally, several possible future improvements were discussed, especially for improving the peak current estimate and the error estimate for peak current, and upgrading the 4DLSS. Some possible approaches for both of these goals were discussed.

Finn, Frank C.↗
DOE OSTI2024

New procedure for evaluation of U(3) coupling and recoupling coefficients

A simple method to calculate Wigner coupling coefficients and Racah recoupling coefficients for U(3) in two group–subgroup chains is presented. While the canonical U(3) coupling and recoupling coefficients are applicable to any system that respects U(3) symmetry, the U(3) coupling coefficients are more specific to nuclear structure studies. This new procedure precludes the use of binomial coefficients and alternating sums which were used in the 1973 formulation of Draayer and Akiyama, and in so doing provides a faster and more accurate determination of any and all required results. The resolution of the outer multiplicity is based on the null space concept of the U(3) generators proposed by Alex et al., whereas the inner multiplicity in the angular momentum subgroup chain is obtained from the dimension of the null space of the SO(3) raising operator. It is anticipated that a C++ library will ultimately be available for determining generic coupling and recoupling coefficients associated with both the canonical and the physical group–subgroup chains of U(3).

Cross-Coupling Reaction↗
DOE OSTI2021

Measurement of the proton spin structure at long distances

Measuring the spin structure of protons and neutrons tests our understanding of how they arise from quarks and gluons, the fundamental building blocks of nuclear matter. At long distances, the coupling constant of the strong interaction becomes large, requiring non-perturbative methods to calculate quantum chromodynamics processes, such as lattice gauge theory or effective field theories. In this work, we report proton spin structure measurements from scattering a polarized electron beam off polarized protons. The spin-dependent cross-sections were measured at large distances, corresponding to the region of low momentum transfer squared between 0.012 and 1.0GeV 2 . This kinematic range provides unique tests of chiral effective field theory predictions. Our results show that a complete description of the nucleon spin remains elusive, and call for further theoretical works, for example, in lattice quantum chromodynamics. Finally, our data extrapolated to the photon point agree with the Gerasimov–Drell–Hearn sum rule, a fundamental prediction of quantum field theory that relates the anomalous magnetic moment of the proton to its integrated spin-dependent cross-sections.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1983

Homogeneous crystal nucleation in binary metallic melts

A method for calculating the homogeneous crystal nucleation frequency in binary metallic melts is developed. The free energy of crystallization is derived from regular solution models for the liquid and solid and is used, together with model-based estimates of the interfacial tension, to calculate the nucleation frequency from the classical theory. The method can account for the composition dependence of the maximum undercooling observed in a number of experiments on small droplet dispersions. It can also be used to calculate the driving force for crystal growth and to obtain more precise estimates of the homogeneous crystal nucleation frequency in glass-forming alloys. This method, although approximate, is simple to apply, and requires only knowledge of the phase diagram and a few readily available thermodynamic quantities as input data.

Thompson, C. V.↗
NASA NTRS1983

Calculation of symmetric vortex separation affecting subsonic bodies at high incidence

A method for calculating the longitudinal aerodynamic coefficients and the pressure distributions on a body at reasonably high angles of attack is presented. The body is represented by a combination of source elements and vortex-lattice elements, including separation of the vortices at increasing angles of attack. The method is self-consistent in that the body and the separated vortex wake are treated as an integrated interacting system. The location of the separation line can be included as an arbitrary input from experimental data or can be evaluated approximately by a pressure-dependent criterion. Calculated values of the aerodynamic coefficients and pressure distributions on cone-cylinder and ogive-cylinder bodies compare well, qualitatively and quantitatively, with experimental data, including simulation of the dependence on Reynolds number.

Almosnino, D.↗
NASA NTRS2007

Asymptotic Energies and QED Shifts for Rydberg States of Helium

This paper reviews progress that has been made in obtaining essentially exact solutions to the nonrelativistic three-body problem for helium by a combination of variational and asymptotic expansion methods. The calculation of relativistic and quantum electrodynamic corrections by perturbation theory is discussed, and in particular, methods for the accurate calculation of the Bethe logarithm part of the electron self energy are presented. As an example, the results are applied to the calculation of isotope shifts for the short-lived 'halo' nucleus He-6 relative to He-4 in order to determine the nuclear charge radius of He-6 from high precision spectroscopic measurements carried out at the Argonne National Laboratory. The results demonstrate that the high precision that is now available from atomic theory is creating new opportunities to create novel measurement tools, and helium, along with hydrogen, can be regarded as a fundamental atomic system whose spectrum is well understood for all practical purposes.

Drake, G.W.F.↗
NASA NTRS1958

Comparison of Calculated and Experimental Temperatures and Coolant Pressure Losses for a Cascade of Small Air-Cooled Turbine Rotor Blades

Average spanwise blade temperatures and cooling-air pressure losses through a small (1.4-in, span, 0.7-in, chord) air-cooled turbine blade were calculated and are compared with experimental nonrotating cascade data. Two methods of calculating the blade spanwise metal temperature distributions are presented. The method which considered the effect of the length-to-diameter ratio of the coolant passage on the blade-to-coolant heat-transfer coefficient and assumed constant coolant properties based on the coolant bulk temperature gave the best agreement with experimental data. The agreement obtained was within 3 percent at the midspan and tip regions of the blade. At the root region of the blade, the agreement was within 3 percent for coolant flows within the turbulent flow regime and within 10 percent for coolant flows in the laminar regime. The calculated and measured cooling-air pressure losses through the blade agreed within 5 percent. Calculated spanwise blade temperatures for assumed turboprop engine operating conditions of 2000 F turbine-inlet gas temperature and flight conditions of 300 knots at a 30,000-foot altitude agreed well with those obtained by the extrapolation of correlated experimental data of a static cascade investigation of these blades.

Stepka, Francis S↗
NASA NTRS1993

Prediction of nonlinear optical properties of organic materials. General theoretical considerations

The prediction of nonlinear optical properties of organic materials is geared to assist materials scientists in the selection of good candidate molecules. A brief summary of the quantum mechanical methods used for estimating hyperpolarizabilities will be presented. The advantages and limitations of each technique will be discussed. Particular attention will be given to the finite-field method for calculating first and second order hyperpolarizabilities, since this method is better suited for large molecules. Corrections for dynamic fields and bulk effects will be discussed in detail, focusing on solvent effects, conformational isomerization, core effects, dispersion, and hydrogen bonding. Several results will be compared with data obtained from third-harmonic-generation (THG) and dc-induced second harmonic generation (EFISH) measurements. These comparisons will demonstrate the qualitative ability of the method to predict the relative strengths of hyperpolarizabilities of a class of compounds. The future application of molecular mechanics, as well as other techniques, in the study of bulk properties and solid state defects will be addressed. The relationship between large values for nonlinear optical properties and large conjugation lengths is well known, and is particularly important for third-order processes. For this reason, the materials with the largest observed nonresonant third-order properties are conjugated polymers. An example of this type of polymer is polydiacetylene. One of the problems in dealing with polydiacetylene is that substituents which may enhance its nonlinear properties may ultimately prevent it from polymerizing. A model which attempts to predict the likelihood of solid-state polymerization is considered, along with the implications of the assumptions that are used. Calculations of the third-order optical properties and their relationship to first-order properties and energy gaps will be discussed. The relationship between monomeric and polymeric third-order optical properties will also be considered.

Cardelino, B.↗