Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “CALCULATION METHODS”

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 415 records · Page 23

Deuteron-Deuteron Elastic and Three and Four-Body Breakup Scattering Using the Faddeev-Yakubovskii Equations

The deuteron-deuteron elastic and three and four-body breakup scattering cross section have been calculated using the Faddeev-Yakubovskii (FY) chain-of-partition momentum-space equations. In this calculation the initial two-cluster potential is split into separable and non-separable components, and the effective potential is reduced to elastic two-body and three and four-body breakup open-channels and a closed-channel many-body contribution. The closed-channel contribution is determined by minimizing a variational bound. The Coulomb interaction was included by expanding the initial and final Coulomb states in a Coulomb-Sturmian basis. The three sets of chain-of-partition integral equations were solved for the elastic and three and four-body breakup scattering amplitudes. The calculations were performed for the S = 2 spin/L = 0 angular-momentum state. The elastic and double-breakup calculations were performed for energies up to E = 5.48 MeV, while the single-breakup calculations were performed for energies up to E = 4.17 MeV. In the case of elastic scattering, the calculated scattering length of 5 a dd = 7.8 ± 0.3 fm is in good agreement with a FY cluster reduction calculation. The calculated phase shift is smaller than that predicted by the resonating group model and this difference is believed to be due to the differences in the potential and calculation methods. The breakup cross sections were calculated as a function of initial deuteron momentum and fragmented-deuteron momentum. Here, the d+d → d+n+p cross sections were compared with neutron yield measurements and, while the measurements also included the L > 0 components, the general features were consistent. Estimates of the calculational uncertainties/bias are provided.

42 ENGINEERING↗

Radiation transport codes for potential applications related to radiobiology and radiotherapy using protons, neutrons, and negatively charged pions

Several Monte Carlo radiation transport computer codes are used to predict quantities of interest in the fields of radiotherapy and radiobiology. The calculational methods are described and comparisions of calculated and experimental results are presented for dose distributions produced by protons, neutrons, and negatively charged pions. Comparisons of calculated and experimental cell survival probabilities are also presented.

Armstrong, T. W.↗

Equilipy

There is high demand for high-throughput calculations of phase equilibria based on the CALPHAD method (CALculation of PHAse Diagram). High-throughput calculations are possible through HPC, currently by a commercial program called Thermo-Calc. However, the number of nodes/processors are limited to the number of purchased license (16 processors per license). This program provides a toolkit for high-throughput calculations of phase equilibria by the CALPHAD method. The program is prepared for Python environment, so that it can be easily installed and used together with other open-source programs. The program can be run in supercomputer.

Kwon, Sunyong [Oak Ridge National Laboratory (ORNL↗

Comparison of spur gear efficiency prediction methods

The predictions of five spur-gear efficiency calculation methods were compared with three sets of test data using different gear geometries. The data and the analysis methods were limited to jet lubricated, ground, spur gears. The data covered a range in pitch line velocity to 1 to 20 m/sec (200 to 4000 ft/min) and K-load factor range of 17 to 1600.

Anderson, N. E.↗

A Data-Driven Global Sensitivity Analysis Framework for Three-Phase Distribution System with PVs

Global sensitivity analysis (GSA) of distribution systems with respect to stochastic PV and load variations plays an important role in designing optimal voltage control schemes. This paper proposes a data-driven framework for GSA of distribution systems. In particular, two representative surrogate modeling-based approaches are developed, including the traditional Gaussian process-based and the analysis of variance (ANOVA) kernel ones. The key idea is to develop a surrogate model that captures the hidden global relationship between voltage and real and reactive power injections from the historical data. With the surrogate model, the Sobol indices can be conveniently calculated through either the sampling-based method or the analytical method to assess the global sensitivity of voltage to variations of PV and load power injections. The sampling-based method approximates the Sobol indices using Monte Carlo simulations while the analytical method calculates them by resorting to the ANOVA expansion framework. Comparison results with other model-based GSA methods on the unbalanced three-phase IEEE 37-bus and 123-bus distribution systems show that the proposed framework can achieve much higher computational efficiency with negligible loss of accuracy. The results on a real 240-node distribution system using actual smart meter data further validate the feasibility and scalability of the proposed framework.

14 SOLAR ENERGY↗

Nonequilibrium radiation and chemistry models for aerocapture vehicle flowfields, volume 3

The computer programs developed to calculate the shock wave precursor and the method of using them are described. This method calculated the precursor flow field in a nitrogen gas including the effects of emission and absorption of radiation on the energy and composition of gas. The radiative transfer is calculated including the effects of absorption and emission through the line as well as the continuum process in the shock layer and through the continuum processes only in the precursor. The effects of local thermodynamic nonequilibrium in the shock layer and precursor regions are also included in the radiative transfer calculations. Three computer programs utilized by this computational scheme to calculate the precursor flow field solution for a given shock layer flow field are discussed.

Carlson, Leland A.↗

BATPAC--VERSION 5.0

Argonne National Laboratory has worked on electrochemical energy storage for several decades. The focus on lithium-ion chemistries started in the early 1990s, developing new materials, synthesis methods, and performance characterizations. Sponsored by the U.S. Department of Energy, Energy Efficiency and Renewable Energy, Vehicle Technologies Office (DOE-EERE-VTO), Argonne has led with many advances. The experimental activities were complemented with multi-scale modeling that ranged from the atomic to the system (manufacturing processes and automobiles) level. The Battery Performance and Cost (BatPaC) model is a calculation method based on Microsoft¿ Office Excel spreadsheets that have been developed at Argonne for estimating the performance and manufacturing cost of lithium-ion batteries for electric-drive vehicles including hybrid-electrics (HEV), plug-in hybrids (PHEV) and pure electrics. The effort is being funded by the Vehicle Technology Office (VTO), which is part of the Energy Efficiency and Renewable Energy (EERE) office of the U.S. Department of Energy (USDOE). BatPaC was first developed in 2007, was subsequently peer-reviewed, and it has served Argonne researchers and the greater battery community in studying the impact of material properties on performance at the pack level. With further developments, the model now allows the design of cells and battery packs for automotive applications, to meet performance requirements (power, energy, recharge time), and estimates the cost of manufacturing the designed batteries. Since the cost depends on the materials he design, and the manufacturing process, this bottom-up model/tool enables the user to study their effects. Designed or the lithium-ion cell and battery researcher, BatPaC helps answer many questions by being 1.Transparent in the assumptions made and the method of calculation 2.Capable of designing a battery specifically for the requirements of an application 3.Constrained by the physical limitations that govern battery performance 4. A bottom-up calculation approach to account for every cost factors. BatPaC predicts the impact of promising materials (and their properties) on the performance metrics relevant for the different applications. Researchers can use the specific capacities and the half-cell voltages of a particular set of electrode materials to calculate the mass and volume of a cell to develop a model incorporating the properties of all the other materials in the cell and the design of the cell enclosure. These calculations not only reveal the impact of an improved material but also enable researchers to calculate the material properties that would be needed to meet the performance criteria of a full battery pack. With this information, researchers can provide he battery industry with realistic expectations that will help it more successfully advance novel battery technologies an Applications.

AHMED, SHABBIR↗

A simplified method for the calculation of airfoil pressure distribution

A method is presented for the rapid calculation of the pressure distribution over an airfoil section when the normal-force distribution and the pressure distribution over the "base profile" (i.e., the profile of the same airfoil were the camber line straight and the resulting airfoil at zero angle of attack) are known. This note is intended as a supplement to N.A.C.A. Report Nos. 631 and 634 wherein methods are presented for the calculation of the normal-force distribution over plain and flapped airfoils, respectively, but not of the pressures on the individual surfaces. Base-profile pressure-coefficient distributions for the usual N.A.C.A. family of airfoils, which are also suitable for several other commonly employed airfoils, are included in tabular form. With these tabulated base-profile pressures and the computed normal-force distributions, pressure distributions adequate for most engineering purposes can be obtained.

Allen, H Julian↗

Accuracy analysis of the semi-analytical method for shape sensitivity calculation

The semianalytical method of design sensitivity analysis that is widely used for calculating derivatives of static response with respect to design variables for structures modeled by finite elements is studied. It is shown that the method can have serious accuracy problems for shape design variables in structures modeled by beam, plate, truss, frame, and solid elements. Errors are shown to be associated with an incompatibility of the sensitivity field with the structure. An error index is developed to test the accuracy of the semianalytical method. It characterizes the difference in errors between a general finite-difference method and the semianalytical method. A method for improving the accuracy of the semianalytical method (when possible) is provided. Examples are presented to demonstrate the use of the error index.

Barthelemy, Bruno↗

Method for the calculation of spacecraft umbra and penumbra shadow terminator points

A method for calculating orbital shadow terminator points is presented. The current method employs the use of an iterative process which is used for an accurate determination of shadow points. This calculation methodology is required since orbital perturbation effects can introduce large errors when a spacecraft orbits a planet in a high altitude and/or highly elliptical orbit. To compensate for the required iteration methodology, all reference frame change definitions and calculations are performed with quaternions. Quaternion algebra significantly reduces the computational time required for the accurate determination of shadow terminator points.

Ortizlongo, Carlos R.↗

A comparison of several methods for the calculation of vibration mode shape derivatives

Four methods for the calculation of derivatives of vibration mode shapes (eigenvectors) with respect to design parameters are reviewed and compared. These methods (finite difference method, Nelson's method, modal method and a modified modal method) are implemented in a general-purpose commercial finite element program and applied to a cantilever beam and a stiffened cylinder with a cutout. A beam tip mass, a beam root height and specific dimensions of the cylinder model comprise the design variables. Data are presented showing the amount of central processor time used to compute the first four eigenvector derivatives for each example problem; errors and rapidity of convergence of the approximate derivative to the exact derivative are taken into account. Nelson's method proved to be most reliable and efficient.

Sutter, T. R.↗

A simple numerical method for the calculation of the laminar boundary layer

An iteration difference method for the calculation of the incompressible laminar boundary layer is described. The method uses Prandtl's boundary layer equation and the boundary conditions directly and permits the attainment of an arbitrary accuracy. The method has been tested successfully in the continuation of the Blasius profile on the flat plate, on the circular cylinder investigated by Heimenz and on an elliptical cylinder of fineness ratio 1:4. The method makes possible the testing of previously developed methods, all of which contain important assumptions.

BOUNDARY LAYER, LAMINAR - CALCULATIONS - SCHRODER ↗

Vibration Analysis Reduces Computer Time

Improved calculation method promises to reduce computer time by factor of 10 for vibration analysis of complex structures. Method originally developed to predict mechanical response of spacecraft to maneuvers in various cases in which same booster is used with different payloads. Method is more exact than conventional normal-mode methods; it involves no initial approximations of assumptions. Calculation cost is reduced without loss of accuracy.

Engels, R. C.↗

Flow mechanism and experimental investigation of a rotating stall in transonic compressors

The flow characteristics of the rotating stall in compressors is studied, and a flow model is developed along with a theoretical calculation method based on vortex theory. A detailed theoretical calculation is completed for a two dimensional flow field in a transonic rotor in a rotating stall, and the result is in good agreement with experimental findings. The oscillograms of time-varying stall characteristic parameters recorded for the onset, growth, and cessation processes of rotating stall are analyzed, and some new flow phenomena deserving of further investigation are discovered. These include serious separation of individual blades, often preceding the onset of rotating stall in compressors with very small blade-camber angles, and periodical variation of the circumferential width of the stall cell with time, accompanied by periodical oscillation of the width of the stall cell in the radial direction of the blade. The circumferential and radial oscillation frequencies are the same.

Yajun, L.↗

Method for Determining Optimum Injector Inlet Geometry

A method for determining the optimum inlet geometry of a liquid rocket engine swirl injector includes obtaining a throttleable level phase value, volume flow rate, chamber pressure, liquid propellant density, inlet injector pressure, desired target spray angle and desired target optimum delta pressure value between an inlet and a chamber for a plurality of engine stages. The method calculates the tangential inlet area for each throttleable stage. The method also uses correlation between the tangential inlet areas and delta pressure values to calculate the spring displacement and variable inlet geometry of a liquid rocket engine swirl injector.

Trinh, Huu P.↗

Computing Viscoplastic Behavior Of A Material

Finite-element implementation developed for Robinson's unified model of viscoplasticity. Implemented via MARC general-purpose finite-element computer program by incorporating all of nonlinearity of material into initial load vector and treating vector as pseudo body force in finite-element equilibrium equations. In technique, each finite increment split into several equal subincrements, constitutive equations integrated by explicit Euler forward-difference method. HYPELA subroutine in MARC includes adaptive integration scheme selecting optimal sizes of subincrements. Spline function smooths discontinuous boundaries in Robinson's model to facilitate finite-element numerical calculations. Method used to calculate elastic, plastic, and creep deformations in solid materials at various temperatures; also used to study thermomechanical behavior of such things as hot pressure vessels and highly stressed components in engines.

Arya, V. K.↗

Relating Land Cover Change to Runoff Distribution Using NASA Earth Observations in Riley County, Kansas

Riley County, Kansas, has observed increased levels of flooding, potentially due to changes in land use/land cover (LULC) and seasonal vegetation variation. This study contrasts two methods of generating runoff curve numbers (CN) from 2006-2020. (1) The traditional Soil Conservation Service CN calculation method uses a look-up table and tracked LULC to determine runoff changes. These tables allow for CN values specific to various land and crop cover types and account for various farming best practices but lack flexibility in calculations for various seasons or plant health. (2) A dynamic method employs normalized difference vegetation index (NDVI) compiled over the rainy season each year to calculate CN using seasonal vegetation. This method allows for a more precise analysis of runoff variability within and between rainy seasons because it can be updated with greater temporal detail and it captures higher spatial resolutions by using NDVI as a proxy for LULC. This study further uses inputs from the United States Geologic Survey (USGS) National Land Cover Database (NLCD), the United States Department of Agriculture (USDA) Cropland Data Layer, and Landsat imagery to create more precise LULC raster datasets including both urban cover and crop-specific land use and curve number maps of the area. Results can guide decision makers in the City of Manhattan, Riley County Department of Planning and Development, Riley County Conservation District, the Kansas Forest Service, and the Kansas Department of Health and Environment toward informed decisions on resiliency strategies to address future flooding.

Ella Griffith↗