Search NASA⌕ Search

SEARCH · Search NASA

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 253 records · Page 14

Calculating the Lightning Protection System Downconductors' Grounding Resistance at Launch Complex 39B, Kennedy Space Center

A new Lightning Protection System (LPS) was designed and built at Launch Complex 39B (LC39B), at the Kennedy Space Center (KSC), Florida, which consists of a catenary wire system (at a height of about 181 meters above ground level) supported by three insulators installed atop three towers in a triangular configuration. Nine downconductors (each about 250 meters long) are connected to the catenary wire system. Each downconductor is connected to a 7.62-meter-radius circular counterpoise conductor with six equally spaced, 6-meter-long vertical grounding rods. Grounding requirements at LC39B call for all underground and aboveground metallic piping, enclosures, raceways, and cable trays, within 7.62 meters of the counterpoise, to be bonded to the counterpoise, which results in a complex interconnected grounding system, given the many metallic piping, raceways, and cable trays that run in multiple directions around LC39B. The complexity of this grounding system makes the fall-of-potential method, which uses multiple metallic rods or stakes, unsuitable for measuring the grounding impedances of the downconductors. To calculate the grounding impedance of the downconductors, an Earth Ground Clamp (EGC) (a stakeless device for measuring grounding impedance) and an Alternative Transient Program (ATP) model of the LPS are used. The EGC is used to measure the loop impedance plus the grounding impedance of each downconductor, and the ATP model is used to calculate the loop impedance of each downconductor circuit. The grounding resistance of the downconductors is then calculated by subtracting the ATP calculated loop impedances from the EGC measurements.

Mata, Carlos T.↗

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.↗

Non-Equilibrium Properties from Equilibrium Free Energy Calculations

Calculating free energy in computer simulations is of central importance in statistical mechanics of condensed media and its applications to chemistry and biology not only because it is the most comprehensive and informative quantity that characterizes the eqUilibrium state, but also because it often provides an efficient route to access dynamic and kinetic properties of a system. Most of applications of equilibrium free energy calculations to non-equilibrium processes rely on a description in which a molecule or an ion diffuses in the potential of mean force. In general case this description is a simplification, but it might be satisfactorily accurate in many instances of practical interest. This hypothesis has been tested in the example of the electrodiffusion equation . Conductance of model ion channels has been calculated directly through counting the number of ion crossing events observed during long molecular dynamics simulations and has been compared with the conductance obtained from solving the generalized Nernst-Plank equation. It has been shown that under relatively modest conditions the agreement between these two approaches is excellent, thus demonstrating the assumptions underlying the diffusion equation are fulfilled. Under these conditions the electrodiffusion equation provides an efficient approach to calculating the full voltage-current dependence routinely measured in electrophysiological experiments.

Pohorille, Andrew↗

Gaseous Nitrogen Orifice Mass Flow Calculator

The Gaseous Nitrogen (GN2) Orifice Mass Flow Calculator was used to determine Space Shuttle Orbiter Water Spray Boiler (WSB) GN2 high-pressure tank source depletion rates for various leak scenarios, and the ability of the GN2 consumables to support cooling of Auxiliary Power Unit (APU) lubrication during entry. The data was used to support flight rationale concerning loss of an orbiter APU/hydraulic system and mission work-arounds. The GN2 mass flow-rate calculator standardizes a method for rapid assessment of GN2 mass flow through various orifice sizes for various discharge coefficients, delta pressures, and temperatures. The calculator utilizes a 0.9-lb (0.4 kg) GN2 source regulated to 40 psia (.276 kPa). These parameters correspond to the Space Shuttle WSB GN2 Source and Water Tank Bellows, but can be changed in the spreadsheet to accommodate any system parameters. The calculator can be used to analyze a leak source, leak rate, gas consumables depletion time, and puncture diameter that simulates the measured GN2 system pressure drop.

Ritrivi, Charles↗

Comparison of Trajectory Models in Calculations of N2-broadened Half-widths and N2-induced Line Shifts for the Rotational Band of H2O-16 and Comparison with Measurements

In this work, Complex Robert-Bonamy calculations of half-widths and line shifts were done for N2-broadening of water for 1639 transitions in the rotational band using two models for the trajectories. The first is a model correct to second order in time, the Robert-Bonamy parabolic approximation. The second is the solution of Hamilton's equations. Both models use the isotropic part of the atom-atom potential to determine the trajectories. The present calculations used an intermolecular potential expanded to 20th order to assure the convergence of the half-widths and line shifts. The aim of the study is to assess if the difference in the half-widths and line shifts determined from the two trajectory models is greater than the accuracy requirements of the spectroscopic and remote sensing communities. The results of the calculations are compared with measurements of the half-widths and line shifts. It is shown that the effects of the trajectory model greatly exceed the needs of current remote sensing measurements and that line shape parameters calculated using trajectories determined by solving Hamilton's equations agree better with measurement.

Lamouroux, J.↗

Combined Uncertainty and A-Posteriori Error Bound Estimates for General CFD Calculations: Theory and Software Implementation

This workshop presentation discusses the design and implementation of numerical methods for the quantification of statistical uncertainty, including a-posteriori error bounds, for output quantities computed using CFD methods. Hydrodynamic realizations often contain numerical error arising from finite-dimensional approximation (e.g. numerical methods using grids, basis functions, particles) and statistical uncertainty arising from incomplete information and/or statistical characterization of model parameters and random fields. The first task at hand is to derive formal error bounds for statistics given realizations containing finite-dimensional numerical error [1]. The error in computed output statistics contains contributions from both realization error and the error resulting from the calculation of statistics integrals using a numerical method. A second task is to devise computable a-posteriori error bounds by numerically approximating all terms arising in the error bound estimates. For the same reason that CFD calculations including error bounds but omitting uncertainty modeling are only of limited value, CFD calculations including uncertainty modeling but omitting error bounds are only of limited value. To gain maximum value from CFD calculations, a general software package for uncertainty quantification with quantified error bounds has been developed at NASA. The package provides implementations for a suite of numerical methods used in uncertainty quantification: Dense tensorization basis methods [3] and a subscale recovery variant [1] for non-smooth data, Sparse tensorization methods[2] utilizing node-nested hierarchies, Sampling methods[4] for high-dimensional random variable spaces.

CFD↗

Contact and Bending Durability Calculation for Spiral-Bevel Gears

The objective of this project is to extend the capabilities of the gear contact analysis solver Calyx, and associated packages Transmission3D, HypoidFaceMilled, HypoidFaceHobbed. A calculation process for the surface durability was implemented using the Dowson-Higginson correlation for fluid film thickness. Comparisons to failure data from NASA's Spiral Bevel Gear Fatigue rig were carried out. A bending fatigue calculation has been implemented that allows the use of the stress-life calculation at each individual fillet point. The gears in the NASA test rig did not exhibit any bending fatigue failure, so the bending fatigue calculations are presented in this report by using significantly lowered strength numbers.

life↗

The Effects of Consistent Chemical Kinetics Calculations on the Pressure-Temperature Profiles and Emission Spectra of Hot Jupiters

In this work we investigate the impact of calculating non-equilibrium chemical abundances consistently with the temperature structure for the atmospheres of highly-irradiated, close-in gas giant exoplanets. Chemical kinetics models have been widely used in the literature to investigate the chemical compositions of hot Jupiter atmospheres which are expected to be driven away from chemical equilibrium via processes such as vertical mixing and photochemistry. All of these models have so far used pressure-temperature (P-T) profiles as fixed model input. This results in a decoupling of the chemistry from the radiative and thermal properties of the atmosphere, despite the fact that in nature they are intricately linked. We use a one-dimensional radiative-convective equilibrium model, ATMO, which includes a sophisticated chemistry scheme to calculate P-T profiles which are fully consistent with non-equilibrium chemical abundances, including vertical mixing and photochemistry. Our primary conclusion is that, in cases of strong chemical disequilibrium, consistent calculations can lead to differences in the P-T profile of up to 100 K compared to the P-T profile derived assuming chemical equilibrium. This temperature change can, in turn, have important consequences for the chemical abundances themselves as well as for the simulated emission spectra. In particular, we find that performing the chemical kinetics calculation consistently can reduce the overall impact of non-equilibrium chemistry on the observable emission spectrum of hot Jupiters. Simulated observations derived from non-consistent models could thus yield the wrong interpretation. We show that this behaviour is due to the non-consistent models violating the energy budget balance of the atmosphere.

planets and satellites↗

Experimental Verification of Buffet Calculation Procedure Using Unsteady Pressure-Sensitive Paint

Typically, a limited number of dynamic pressure sensors is employed to determine the unsteady aerodynamic forces on large, slender aerospace structures. This paper describes a robust calculation procedure based on frequency-specific correlation lengths and validation from an experiment conducted on a flat panel coated with fast response pressure-sensitive paint. The first part of the paper describes the procedure used to analyze the pressure sensitive paint images and a calibration method using dynamic pressure transducers. Excellent comparison in spectra, coherence, and phase, measured via pressure-sensitive paint and dynamic pressure sensors, validates the pressure-sensitive paint data. The second part of the paper describes the buffet validation process, the first step of which was to use pressure histories from all pixels to determine the true force fluctuations. In the next step, only a selected number of pixels was chosen as virtual sensors, and a correlation-length-based buffet calculation procedure was applied to determine modeled force fluctuations. By progressively decreasing the number of virtual sensors, it was observed that the present calculation procedure was able to make a close estimate of the true unsteady forces only from eight sensors. It is believed that the present work provides the first validation of a buffet calculation procedure.

Buffet↗

Structure Deformation Calculation Program Based on Displacement Theory for Shape Predictions

Separated programs were written in C/C++ to validate the Displacement Transfer Functions. The Structure Deformation Calculation Program was written to combine all of the programs to calculate deformed shapes of a structure using surface strain data and structural geometrical parameters. Users do not need to know the material properties, nor the complex internal structures geometry because the Displacement Theory is purely geometrical in nature. Users only need to know the structure types as defined in this report and information such as the structure length, depth factors, number of strain sensors, and the surface strains measured at the strain-sensing stations installed on the structures. Depending on the structure type, an applicable Displacement Transfer Function will be used. This program requires two input files created by users; the recorded strain data file in comma-separated values format and the structure geometry data file in text format. The program will output the out-of-plane deflections, slopes, cross-sectional twist angles, and depth factors if applicable. All output files are created in comma-separated values format. A section in this report describes step-by-step procedures on how to use the Structure Deformation Calculation Program for structure deformed shape calculations.

Displacement theory↗

Calculating Mode Mixity in Interlaminar Fracture of Composite Laminates

To accurately predict delamination under mixed-mode loading, cohesive elements must partition the energies dissipated in Mode I and Mode II correctly. The ability of cohesive models to predict mode mixity was evaluated by performing analyses of the single leg bending (SLB) test, which exhibits a variable mode mixity with crack length. Accordingly, two numerical integration methods were developed to compute the energy partitioning: (i) a local calculation of the accumulated work of tractions at each point along the crack path, and (ii) a global J-integral “snapshot” calculation along a contour around the process zone. In the case of a composite laminate with properties characterized by short fracture process zone (FPZ), the results of these two methods correlate well with each other and with a linear elastic fracture mechanics solution. However, FPZscan be long due to fiber bridging or due to the intrinsic toughness of the material. Then, the results of mode-mixity calculations obtained with various methods vary widely from each other. It is argued that the global J-integral calculation is the most accurate interpretation of the energy release rate and mode mixity for either short or long process zones. This solution will provide a more accurate characterization of resistance curves and mode mixity in experimental results, and it will also guide the development of more accurate mixed-mode cohesive laws for tougher materials characterized by strong R-curve responses, such as thermo-plastic composites.

cohesive zone analysis↗

Computer program calculates and plots surface area and pore size distribution data

Computer program calculates surface area and pore size distribution of powders, metals, ceramics, and catalysts, and prints and plots the desired data directly. Surface area calculations are based on the gas adsorption technique of Brunauer, Emmett, and Teller, and pore size distribution calculations are based on the gas adsorption technique of Pierce.

Halpert, G.↗

Program for calculating optimum dimensions of alpha radioisotope capsules exposed to varying stress and temperature

A method and computer program were developed for calculating the creep and optimizing the dimensions of capsules filled with alpha-emitting radioisotopes. The method solves an integral equation that was developed assuming linear accumulation of partial creep lives and relating life to time-dependent stress and temperature using the Larson-Miller parameter. The computer program, CAPSUL, is written in Fortran language for the IBM 360/75 computer. The program makes a least squares fit of the creep life function using conventional constant stress, constant temperature creep data. Dimensions of capsules having maximum thermal power per unit of weight, volume, or area are calculated for a given creep life and pressure-temperature history using a numerical Lagrange Multiplier formulation. The program also calculates the life to a prescribed strain for capsules of given dimensions and pressure-temperature history. The method has been used to analyze creep data for the alloys 304 stainless steel, Hastelloy N, Cb-1% Zr, FS-85, and T-222.

J. P. Nichols↗

A fast algorithm for the calculation of junction capacitance and its application for impurity profile determination

A fast algorithm is described which calculates the space charge layer width and junction capacitance for an arbitrary impurity profile and for plane, cylindrical and spherical junctions. The algorithm is based on the abrupt space charge edge (ASCE) approximation. A method to use the algorithm for the determination of impurity profiles for two-sided junctions is presented. An expression is derived for the built-in voltage to be used for capacitance calculations with the ASCE approximation. Experimental evidence is given that the algorithm permits very accurate capacitance calculations and also predicts the exact temperature dependence of the junction capacitance.

Deman, H. J. J.↗

Deformed Brueckner-Hartree-Fock calculation for light nuclei

For the first time the Brueckner-Hartree-Fock (BHF) method was applied to nuclei whose intrinsic structure is nonspherical. One aim was to investigate whether the energy dependent reaction matrix calculated from a realistic nucleon-nucleon interaction leads to deformations similar to, or different from, those obtained from energy independent interactions in Hartree-Fock (HF) calculations. Reaction matrix elements were calculated as a function of starting energy for the Hamada-Johnston interaction, using a Pauli operator appropriate to O-16 and a shifted oscillator spectrum for virtual excited states. Binding energies, single-particle energies, radii, and shape deformations of the intrinsic state in unrenormalized as well as renormalized BHF are discussed and compared with previous HF studies. Results are presented for C-12, O-16, and Ne-20.

Braley, R. C.↗

Monte Carlo calculations of high energy nucleon meson cascades and applications to galactic cosmic ray transport

Results obtained using a recently developed calculational method for determining the nucleon-meson cascade induced in thick materials by high-energy nucleons and charged pions are presented. The calculational method uses the intranuclear-cascade-evaporation model to treat nonelastic collisions by particles with energies approximately or smaller than GeV and an extrapolation model at higher energies. The following configurations are considered: (1) 19.2-GeV/c protons incident on iron; (2) 30.3-GeV/c protons incident on iron; (3) solar and galactic protons incident on the moon, and (4) galactic protons incident on tissue. For the first three configurations, experimental results are available and comparisons between the experimental and calculated results are given.

Armstrong, T. W.↗

MORSE Monte Carlo shielding calculations for the zirconium hydride reference reactor

Verification of DOT-SPACETRAN transport calculations of a lithium hydride and tungsten shield for a SNAP reactor was performed using the MORSE (Monte Carlo) code. Transport of both neutrons and gamma rays was considered. Importance sampling was utilized in the MORSE calculations. Several quantities internal to the shield, as well as dose at several points outside of the configuration, were in satisfactory agreement with the DOT calculations of the same.

Burgart, C. E.↗