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 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
NASA NTRS2003

Numerical Modeling of Fluid Transient in Cryogenic Fluid Network of Rocket Propulsion System

Fluid transients, also known as water hammer, can have a significant impact on the design and operation of both spacecraft and launch vehicles propulsion systems. These transients often occur at system activation and shut down. For ground safety reasons, many spacecrafts are launched with the propellant lines dry. These lines are often evacuated by the time the spacecraft reaches orbit. When the propellant isolation valve opens during propulsion system activation, propellant rushes into lines creating a pressure surge. During propellant system shutdown, a pressure surge is created due to sudden closure of a valve. During both activation and shutdown, pressure surges must be predicted accurately to ensure structural integrity of the propulsion system fluid network. The method of characteristics is the most widely used method of calculating fluid transients in pipeline [ 1,2]. The method of characteristics, however, has limited applications in calculating flow distribution in complex flow circuits with phase change, heat transfer and rotational effects. A robust cryogenic propulsion system analyzer must have the capability to handle phase change, heat transfer, chemical reaction, rotational effects and fluid transients in conjunction with subsystem flow model for pumps, valves and various pipe fittings. In recent years, such a task has been undertaken at Marshall Space Flight Center with the development of the Generalized Fluid System Simulation Program (GFSSP), which is based on finite volume method in fluid network [3]. GFSSP has been extensively verified and validated by comparing its predictions with test data and other numerical methods for various applications such as internal flow of turbo-pump [4], propellant tank pressurization [5,6], chilldown of cryogenic transfer line [7] and squeeze film damper rotordynamics [8]. The purpose of the present paper is to investigate the applicability of the finite volume method to predict fluid transient in cryogenic flow circuits.

Majumdar, Alok↗
NASA NTRS1934

The calculation of take-off run

A comparatively simple method of calculating length of take-off run is developed from the assumption of a linear variation in net accelerating force with air speed and it is shown that the error involved is negligible.

Diehl, Walter S↗
NASA NTRS1962

Estimation of Critical Temperature for Surface Ion Currents from Electron Emission Data

A method of calculating a relation for the critical surface temperature for ion current emission is presented. The method is based on the S-shaped electron emission curves for surfaces in the presence of ionizable vapors and upon the assumptions of thermodynamic equilibrium involved in the Saha-Langmuir relation. A comparison of the critical temperatures so calculated with the relation obtained from actual ion emission data on the cesium-tungsten system shows good agreement over a wide current density range. Critical temperatures for cesium ion current densities of 0.21 and 1.9 amperes per square centimeter calculated by this method are presented for surfaces of rhenium, molybdenum, tantalum, and niobium. Tentative relations of the form log j equals A (sub p) plus (B (sub p) over T) are presented for these same systems, where j is the ion current density, A (sub p) and B (sub p) are constants, and T is the temperature.

ELECTRON EMISSION↗
NASA NTRS1961

Calculation of Wind Compensation for Launching of Unguided Rockets

A method for calculating wind compensation for unguided missiles s derived which has a greater degree of flexibility than the previously proposed methods. Utilization of the wind-compensation technique is demonstrated by using the Shotput vehicle as a model. Postflight simulations of four of these missiles with the use of measured winds show that if the winds are known, very good accuracy can be obtained by using the proposed method.

James, Robert L., Jr.↗
NASA NTRS1992

Three-dimensional boundary layer calculation by a characteristic method

A numerical method for solving the three-dimensional boundary layer equations for bodies of arbitrary shape is presented. In laminar flows, the application domain extends from incompressible to hypersonic flows with the assumption of chemical equilibrium. For turbulent boundary layers, the application domain is limited by the validity of the mixing length model used. In order to respect the hyperbolic nature of the equations reduced to first order partial derivative terms, the momentum equations are discretized along the local streamlines using of the osculator tangent plane at each node of the body fitted coordinate system. With this original approach, it is possible to overcome the use of the generalized coordinates, and therefore, it is not necessary to impose an extra hypothesis about the regularity of the mesh in which the boundary conditions are given. By doing so, it is possible to limit, and sometimes to suppress, the pre-treatment of the data coming from an inviscid calculation. Although the proposed scheme is only semi-implicit, the method remains numerically very efficient.

Houdeville, R.↗
DOE OSTI2020

Unsaturated and Saturated Hydraulic Parameters of Soils at the East River Watershed, Colorado Calculated using the Pedotransfer Function Method.

The purpose of the data package is to present the results of calculations of the following unsaturated-saturated hydraulic parameters of soils of the East River Watershed, Colorado: fully saturated water content (Qs), residual saturation (Qr), van genuchten parameters alpha and n, field capacity (FC), and saturated hydraulic conductivity (Ks). Calculations were conducted for the following soil texture types: sandy_loam, loam, silty_clay_loam, clay_loam, silt loam; the particle size distribution included the percentage of sand, silt, clay, and the bulk density (BD) of soils. Calculations were conducted for the following areas- Lower_Subalpine (6 samples),- Lower_Montane_lower_hillslope (7 samples),- Floodplain_at_base_of_lower_montane_hillslope, (11 samples),- PLM-0 (10 samples),- PLM-1 (10 samples), and - PLM-3 (13 samples)Calculations were conducted using the Pedotransfer Function Method developed by Tóth et al. (2014) with the application of the code ”euptf” in Rstudio. The data are being summarized in the Excel file ”Unsaturated_Saturated_Hydraulic_Parameters_PTF_calculated_EASt_EIver.xls" including two worksheets:“full_data_EAST_RIVER” which includes (a) the input data used for calculations, the soil texture triangle, and the plots of the water retention curves and unsaturated hydraulic conductivity for several types of soil texture, and (b) the summary table with the averaged parameters for each area. The the texture triangle was plotted using the Rstudio "soiltexture." Microsoft Excel software is required to open the Excel file. The metadata and a list of references to the Rstudio libraries used for calculations and data visualization are giving in the Word Readme file.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1959

Supersonic and Moment-of-Area Rules Combined for Rapid Zero-Lift Wave-Drag Calculations

The concepts of the supersonic area rule and the moment-of-area rule are combined to develop a new method for calculating zero-lift wave drag which is amenable to the use of ordinary desk calculators. The total zero-lift wave drag of a configuration is calculated by the new method as the sum of the wave drag of each component alone plus the interference between components. In calculating the separate contributions each component or pair of components is analyzed over the smallest allowable length in order to improve the convergence of the series expression for the wave drag. The accuracy of the present method is evaluated by comparing the total zero-lift wave-drag solutions for several simplified configurations obtained by the present method with solutions given by slender-body and linearized theory. The accuracy and computational time required by the present method are also evaluated relative to the supersonic area rule and the moment-of-area rule. The results of the evaluation indicate that total zero-lift wave-drag solutions for simplified configurations can be obtained by the present method which differ from solutions given by slender-body and linearized theory by less than 6 percent. This accuracy for simplified configurations was obtained from only nine terms of the series expression for the wave drag as a result of calculating the total zero-lift wave drag by parts. For the same number of terms these results represent an accuracy greater than that for solutions obtained by either of the two methods upon which the present method is based, except in a few isolated cases. For the excepted cases, solutions by the present method and the supersonic area rule are identical. Solutions by the present method are obtained in one fifth the computing time required by the supersonic area rule. This difference in computing time of course would be substantially reduced if the complete procedures for both methods were programmed on electronic computing machines.

Levy, Lionel L., Jr.↗
NASA NTRS1989

Application of a full potential method to AGARD standard airfoils

One of the most important uses of method that calculate unsteady aerodynamic loads is to predict and analyze the aeroelastic responses of flight vehicles. Currently, methods based on transonic small disturbance potential aerodynamics are the primary tools for aeroelastic analysis. Flow solutions obtained using isentropic potential theory can be highly inaccurate and even multivalued, because they do not model the effects of entropy that is produced when shock waves are in the flow field. From the results that are presented, it is concluded that nonisentropic potential methods more accurately model Euler solutions than do isentropic methods. The primary effects of modeling shock generated entropy are: (1) to eliminate mulitple flow solutions when strong shock waves are in the flow field; and (2) to bring the strengths and locations of computed shock waves into better agreement with those calculated using Euler method and those measured during experiments.

Whitlow, Woodrow, Jr.↗
NASA NTRS1985

Preconditioned minimal residual methods for Chebyshev spectral calculations

The problem of preconditioning the pseudospectral Chebyshev approximation of an elliptic operator is considered. The numerical sensitiveness to variations of the coefficients of the operator are investigated for two classes of preconditioning matrices: one arising from finite differences, the other from finite elements. The preconditioned system is solved by a conjugate gradient type method, and by a Dufort-Frankel method with dynamical parameters. The methods are compared on some test problems with the Richardson method and with the minimal residual Richardson method.

Canuto, C.↗
NASA NTRS1980

Calculation of two dimensional vortex/surface interference using panel methods

The application of panel methods to the calculation of vortex/surface interference characteristics in two dimensional flow was studied over a range of situations starting with the simple case of a vortex above a plane and proceeding to the case of vortex separation from a prescribed point on a thick section. Low order and high order panel methods were examined, but the main factor influencing the accuracy of the solution was the distance between control stations in relation to the height of the vortex above the surface. Improvements over the basic solutions were demonstrated using a technique based on subpanels and an applied doublet distribution.

Maskew, B.↗
DOE OSTI2020

Lattice Physics Calculations Using the Embedded Self-Shielding Method in Polaris, Part I: Methods and Implementation

Polaris is a 2-dimensional multigroup lattice physics capability in the SCALE code system for the analysis of light water reactor fuel designs. The goal of light water reactor lattice physics codes is to generate few-group homogenized cross sections for downstream full-core nodal diffusion calculations. Additionally, lattice physics calculations contain three primary components: the cross section processing calculation, the 2D transport calculation, and the depletion calculation. This paper summarizes the calculational methods and their implementation into Polaris, with an emphasis on implementation of the embedded self-shielding method. The accuracy of the embedded self-shielding method depends on the procedure used to generate self-shielding factors on the multigroup library. Numerical benchmarks calculations reveal that the accuracy of Polaris eigenvalue predictions is enhanced by (1) using heterogeneous unit cell models to generate the self-shielding factors on the library and (2) using self-shielding factors for within-group scattering cross section.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1981

Calculation of guaranteed mean power from wind turbine generators

A method for calculating the 'guaranteed mean' power output of a wind turbine generator is proposed. The term 'mean power' refers to the average power generated at specified wind speeds during short-term tests. Correlation of anemometers, the method of bins for analyzing non-steady data, the PROP Code for predicting turbine power, and statistical analysis of deviations in test data from theory are discussed. Guaranteed mean power density for the Clayton Mod-OA system was found to be 8 watts per square meter less than theoretical power density at all power levels, with a confidence level of 0.999. This amounts to 4 percent of rated power.

Spera, D. A.↗
NASA NTRS1987

Calculation of aerodynamic heat transfer for blunt-nosed thin wings at different angles of attack at supersonic speeds

A method for calculating aerodynamic heat transfer on thin wings at angles of attack is provided, based on the assumption of small cross flow and the calculation of an infinite cylinder. It is valid in the range of supersonic and hypersonic speeds. The method does not require calculation of the details of the stream lines and can be used for various swept wings at moderate angle of attack or yaw angle. A comparison between theoretical and experimental results is given.

Xu, Shuxuan↗
NASA NTRS1977

Calculation of transonic flow over supercritical airfoil sections

Two numerical methods for calculating the transonic flow, including viscous effects over lifting airfoil sections and experimental data are compared for turbulent flow over a supercritical airfoil. In addition, results for a NACA 64A010 airfoil at nonzero angle of attack are compared to demonstrate the applicability of the numerical methods to classical, lifting airfoils. One numerical method is a solution to the time-averaged Navier-Stokes equations throughout the entire flow field. The other is a hybrid method that combines inviscid, boundary-layer, and Navier-Stokes equations in appropriate regions of the flow field. Both methods adequately predict the surface pressures and flow field about the 64A010 airfoil at M = 0.8 and alpha = 2 deg when an appropriate turbulence model is used. The methods are not as successful for the supercritical airfoil.

Rose, W. C.↗
NASA NTRS1974

Calculations of the pressure distribution on axisymmetric boattails including effects of viscous interactions and exhaust jets in subsonic flow

A method of calculating the pressure distributions on boattails is proposed. This method accounts for viscous effects including the presence of a separated region for base flows by combining an inviscid analysis with a boundary layer analysis in an iterative calculation. Details of the reversed flow region are not considered. Some preliminary results have been obtained for boattails at subsonic free stream Mach number with turbulent boundary layers separating at the boattail base. In some cases convergence could not be obtained using the present computer program. It is possible, in principle, to extend this method to the calculation of boattail flows with pressure gradient induced separation on the boattail.

Rom, J.↗
NASA NTRS1996

Study of the Radiative Properties of Inhomogeneous Stratocumulus Clouds

Clouds play an important role in the radiation budget of the atmosphere. A good understanding of how clouds interact with solar radiation is necessary when considering their effects in both general circulation models and climate models. This study examined the radiative properties of clouds in both an inhomogeneous cloud system, and a simplified cloud system through the use of a Monte Carlo model. The purpose was to become more familiar with the radiative properties of clouds, especially absorption, and to investigate the excess absorption of solar radiation from observations over that calculated from theory. The first cloud system indicated that the absorptance actually decreased as the cloud's inhomogeneity increased, and that cloud forcing does not indicate any changes. The simplified cloud system looked at two different cases of absorption of solar radiation in the cloud. The absorptances calculated from the Monte Carlo is compared to a correction method for calculating absorptances and found that the method can over or underestimate absorptances at cloud edges. Also the cloud edge effects due to solar radiation points to a possibility of overestimating the retrieved optical depth at the edge, and indicates a possible way to correct for it. The effective cloud fraction (Ne) for a long time has been calculated from a cloud's reflectance. From the reflectance it has been observed that the N, for most cloud geometries is greater than the actual cloud fraction (Nc) making a cloud appear wider than it is optically. Recent studies we have performed used a Monte Carlo model to calculate the N, of a cloud using not only the reflectance but also the absorptance. The derived Ne's from the absorptance in some of the Monte Carlo runs did not give the same results as derived from the reflectance. This study also examined the inhomogeneity of clouds to find a relationship between larger and smaller scales, or wavelengths, of the cloud. Both Fourier transforms and wavelet transforms were used to analyze the liquid water content of marine stratocumulus clouds taken during the ASTEX project. From the analysis it was found that the energy in the cloud is not uniformly distributed but is greater at the larger scales than at the smaller scales. This was determined by examining the slope of the power spectrum, and by comparing the variability at two scales from a wavelet analysis.

Batey, Michael↗
DOE OSTI2024

The Effective Fragment Molecular Orbital Method: Achieving High Scalability and Accuracy for Large Systems

The effective fragment molecular orbital (EFMO) method has been developed to predict the total energy of a very large molecular system accurately (with respect to the underlying quantum mechanical method) and efficiently by taking advantage of the locality of strong chemical interactions and employing a two-level hierarchical parallelism. The accuracy of the EFMO method is partly attributed to the accurate and robust intermolecular interaction prediction between distant fragments, in particular, the many-body polarization and dispersion effects, which require the generation of static and dynamic polarizability tensors by solving the coupled perturbed Hartree–Fock (CPHF) and time-dependent HF (TDHF) equations, respectively. Solving the CPHF and TDHF equations is the main EFMO computational bottleneck due to the inefficient (serial) and I/O-intensive implementation of the CPHF and TDHF solvers. In this work, the efficiency and scalability of the EFMO method are significantly improved with a new CPU memory-based implementation for solving the CPHF and TDHF equations that are parallelized by either message passing interface (MPI) or hybrid MPI/OpenMP. Here, the accuracy of the EFMO method is demonstrated for both covalently bonded systems and noncovalently bound molecular clusters by systematically examining the effects of basis sets and a key distance-related cutoff parameter, R cut . R cut determines whether a fragment pair (dimer) is treated by the chosen ab initio method or calculated using the effective fragment potential (EFP) method (separated dimers). Decreasing the value of Rcut increases the number of separated (EFP) dimers, thereby decreasing the computational effort. It is demonstrated that excellent accuracy (<1 kcal/mol error per fragment) can be achieved when using a sufficiently large basis set with diffuse functions coupled with a small R cut value. With the new parallel implementation, the total EFMO wall time is substantially reduced, especially with a high number of MPI ranks. Given a sufficient workload, nearly ideal strong scaling is achieved for the CPHF and TDHF parts of the calculation. For the first time, EFMO calculations with the inclusion of long-range polarization and dispersion interactions on a hydrated mesoporous silica nanoparticle with explicit water solvent molecules (more than 15k atoms) are achieved on a massively parallel supercomputer using nearly 1000 physical nodes. In addition, EFMO calculations on the carbinolamine formation step of an amine-catalyzed aldol reaction at the nanoscale with explicit solvent effects are presented.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗