Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Calculations”

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 901 records · Page 50

Heats of Segregation of BCC Binaries from ab Initio and Quantum Approximate Calculations

We compare dilute-limit heats of segregation for selected BCC transition metal binaries computed using ab initio and quantum approximate energy methods. Ab initio calculations are carried out using the CASTEP plane-wave pseudopotential computer code, while quantum approximate results are computed using the Bozzolo-Ferrante-Smith (BFS) method with the most recent LMTO-based parameters. Quantum approximate segregation energies are computed with and without atomistic relaxation, while the ab initio calculations are performed without relaxation. Results are discussed within the context of a segregation model driven by strain and bond-breaking effects. We compare our results with full-potential quantum calculations and with available experimental results.

Good, Brian S.↗

Surface Segregation Energies of BCC Binaries from Ab Initio and Quantum Approximate Calculations

We compare dilute-limit segregation energies for selected BCC transition metal binaries computed using ab initio and quantum approximate energy method. Ab initio calculations are carried out using the CASTEP plane-wave pseudopotential computer code, while quantum approximate results are computed using the Bozzolo-Ferrante-Smith (BFS) method with the most recent parameterization. Quantum approximate segregation energies are computed with and without atomistic relaxation. The ab initio calculations are performed without relaxation for the most part, but predicted relaxations from quantum approximate calculations are used in selected cases to compute approximate relaxed ab initio segregation energies. Results are discussed within the context of segregation models driven by strain and bond-breaking effects. We compare our results with other quantum approximate and ab initio theoretical work, and available experimental results.

Good, Brian S.↗

Calculated X-ray Intensities Using Monte Carlo Algorithms: A Comparison to Experimental EPMA Data

Monte Carlo (MC) modeling has been used extensively to simulate electron scattering and x-ray emission from complex geometries. Here are presented comparisons between MC results and experimental electron-probe microanalysis (EPMA) measurements as well as phi(rhoz) correction algorithms. Experimental EPMA measurements made on NIST SRM 481 (AgAu) and 482 (CuAu) alloys, at a range of accelerating potential and instrument take-off angles, represent a formal microanalysis data set that has been widely used to develop phi(rhoz) correction algorithms. X-ray intensity data produced by MC simulations represents an independent test of both experimental and phi(rhoz) correction algorithms. The alpha-factor method has previously been used to evaluate systematic errors in the analysis of semiconductor and silicate minerals, and is used here to compare the accuracy of experimental and MC-calculated x-ray data. X-ray intensities calculated by MC are used to generate a-factors using the certificated compositions in the CuAu binary relative to pure Cu and Au standards. MC simulations are obtained using the NIST, WinCasino, and WinXray algorithms; derived x-ray intensities have a built-in atomic number correction, and are further corrected for absorption and characteristic fluorescence using the PAP phi(rhoz) correction algorithm. The Penelope code additionally simulates both characteristic and continuum x-ray fluorescence and thus requires no further correction for use in calculating alpha-factors.

Carpenter, P. K.↗

Calculation of Atomic Data for NASA Missions

The interpretation of cosmic spectra relies on a vast sea of atomic data which are not readily obtainable from analytic expressions or simple calculations. Rather, their evaluation typically requires state-of-the-art atomic physics calculations, with the inclusion of weaker effects (spin-orbit and configuration interactions, relaxation, Auger broadening, etc.), to achieve the level of accuracy needed for use by astrophysicists. Our NASA-supported research program is focused on calculating data for three important atomic processes, 1) dielectronic recombination (DR), 2) inner-shell photoabsorption, and 3) fluorescence and Auger decay of inner-shell vacancy states. Some additional details and examples of our recent findings are given.

Gorczyca, T. W.↗

Numerical Strip-Yield Calculation of CTOD

Many existing methods to calculate CTOD can be costly and complicated, or apply only to particular configurations. A new numerical method for calculating CTOD was investigated. NASGRO's Boundary Element module NASBEM was adapted to calculate displacements at any point on the crack. Demonstrated for a number of crack configurations: a) finite and infinite domains; b) center and edge cracks; and c) complex cases with several cracks and holes. Great accuracy at minimal computational cost.

Beek, Joachim↗

Blocked Force and Loading Calculations for LaRC THUNDER Actuators

An analytic approach is developed to predict the performance of LaRC Thunder actuators under load and under blocked conditions. The problem is treated with the Von Karman non-linear analysis combined with a simple Raleigh-Ritz calculation. From this, shape and displacement under load combined with voltage are calculated. A method is found to calculate the blocked force vs voltage and spring force vs distance. It is found that under certain conditions, the blocked force and displacement is almost linear with voltage. It is also found that the spring force is multivalued and has at least one bifurcation point. This bifurcation point is where the device collapses under load and locks to a different bending solution. This occurs at a particular critical load. It is shown this other bending solution has a reduced amplitude and is proportional to the original amplitude times the square of the aspect ratio.

Campbell, Joel F.↗

Rapid Calculations of Three-Dimensional Inlet/Fan Interaction

Two computational fluid dynamics codes have been merged to permit rapid calculations of inlet/fan interaction. Inlets are modeled using the WIND-US Navier-Stokes code. Fans are modeled using a new three-dimensional Euler code called CSTALL that solves the flow through the entire compression system but models blade rows using body forces for turning and loss. The body force model is described and it is shown how unknown terms in the model can be estimated from other Navier-Stokes solutions of the blade rows run separately. The inlet and fan calculations are run simultaneously and are coupled at an interface plane using a third code called SYNCEX that is described briefly. Results are shown for an axisymmetric nacelle at high angle of attack modeled both as an isolated inlet and coupled to a single stage fan. The isolated inlet calculations are unrealistic after the flow separates but the coupled codes can model large regions of separated flow extending from the lower lip of the nacelle into the fan rotor.

Chima, Rodrick V.↗

Assessment of Comprehensive Analysis Calculation of Airloads on Helicopter Rotors

Blade section normal force and pitching moment were investigated for six rotors operating at transition and high speeds: H-34 in flight and wind tunnel, SA 330 (research Puma), SA 349/2, UH-60A full-scale and BO-105 model (HART-I). The measured data from flight and wind tunnel tests were compared with calculations obtained using the comprehensive analysis CAMRAD II. The calculations were made using two free wake models: rolled-up and multiple-trailer with consolidation models. At transition speed, there is fair to good agreement for the blade section normal force between the test data and analysis for the H-34, research Puma, and SA 349/2 with the rolled-up wake. The calculated airloads differ significantly from the measurements for the UH-60A and BO-105. Better correlation is obtained for the UH-60A and BO-105 by using the multiple-trailer with consolidation wake model. In the high speed condition, the analysis shows generally good agreement with the research Puma flight data in both magnitude and phase. However, poor agreement is obtained for the other rotors examined. The analysis shows that the aerodynamic tip design (chord length and quarter chord location) of the Puma has an important influence on the phase correlation.

Yeo, Hyeonsoo↗

Calculational and Experimental Investigations of the Pressure Effects on Radical - Radical Cross Combinations Reactions: C2H5 + C2H3

Pressure-dependent product yields have been experimentally determined for the cross-radical reaction C2H5 + C2H3. These results have been extended by calculations. It is shown that the chemically activated combination adduct, 1-C4H8*, is either stabilized by bimolecular collisions or subject to a variety of unimolecular reactions including cyclizations and decompositions. Therefore the "apparent" combination/disproportionation ratio exhibits a complex pressure dependence. The experimental studies were performed at 298 K and at selected pressures between about 4 Torr (0.5 kPa) and 760 Torr (101 kPa). Ethyl and vinyl radicals were simultaneously produced by 193 nm excimer laser photolysis of C2H5COC2H3 or photolysis of C2H3Br and C2H5COC2H5. Gas chromatograph/mass spectrometry/flame ionization detection (GC/MS/FID) were used to identify and quantify the final reaction products. The major combination reactions at pressures between 500 (66.5 kPa) and 760 Torr are (1c) C2H5 + C2H3 yields 1-butene, (2c) C2H5 + C2H5 yields n-butane, and (3c) C2H3 + C2H3 yields 1,3-butadiene. The major products of the disproportionation reactions are ethane, ethylene, and acetylene. At moderate and lower pressures, secondary products, including propene, propane, isobutene, 2-butene (cis and trans), 1-pentene, 1,4-pentadiene, and 1,5-hexadiene are also observed. Two isomers of C4H6, cyclobutene and/or 1,2-butadiene, were also among the likely products. The pressure-dependent yield of the cross-combination product, 1-butene, was compared to the yield of n-butane, the combination product of reaction (2c), which was found to be independent of pressure over the range of this study. The [ 1-C4H8]/[C4H10] ratio was reduced from approx.1.2 at 760 Torr (101 kPa) to approx.0.5 at 100 Torr (13.3 kPa) and approx.0.1 at pressures lower than about 5 Torr (approx.0.7 kPa). Electronic structure and RRKM calculations were used to simulate both unimolecular and bimolecular processes. The relative importance of C-C and C-H bond ruptures, cyclization, decyclization, and complex decompositions are discussed in terms of energetics and structural properties. The pressure dependence of the product yields were computed and dominant reaction paths in this chemically activated system were determined. Both modeling and experiment suggest that the observed pressure dependence of [1-C4H8]/[C4H10] is due to decomposition of the chemically activated combination adduct 1-C4H8* in which the weaker allylic C-C bond is broken: H2C=CHCH2CH3 yields C3H5 + CH3. This reaction occurs even at moderate pressures of approx.200 Torr (26 kPa) and becomes more significant at lower pressures. The additional products detected at lower pressures are formed from secondary radical-radical reactions involving allyl, methyl, ethyl, and vinyl radicals. The modeling studies have extended the predictions of product distributions to different temperatures (200-700 K) and a wider range of pressures (10(exp -3) - 10(exp 5) Torr). These calculations indicate that the high-pressure [1-C4H8]/[C4H10] yield ratio is 1.3 +/- 0.1.

Fahr, Askar↗

The Calculated Effect of Various Hydrodynamic and Aerodynamic Factors on the Take-Off of a Large Flying Boat

Present designs for large flying boats are characterized by high wing loading, high aspect ratio, and low parasite drag. The high wing loading results in the universal use of flaps for reducing the takeoff and landing speeds. These factors have an effect on takeoff performance and influence to a certain extent the design of the hull. An investigation was made of the influence of various factors and design parameters on the takeoff performance of a hypothetical large flying boat by means of takeoff calculations. The parameters varied in the calculations were size of hull (load coefficient), wing setting, trim, deflection of flap, wing loading, aspect ratio, and parasite drag. The takeoff times and distances were calculated to the stalling speeds and the performance above these speeds was studied separately to determine piloting technique for optimum takeoff. The advantage of quick deflection of the flap at high water speeds is shown.

Olson, R.E.↗

Calculation of Oxygen Fugacity in High Pressure Metal-Silicate Experiments and Comparison to Standard Approaches

Calculation of oxygen fugacity in high pressure and temperature experiments in metal-silicate systems is usually approximated by the ratio of Fe in the metal and FeO in the silicate melt: (Delta)IW=2*log(X(sub Fe)/X(sub FeO)), where IW is the iron-wustite reference oxygen buffer. Although this is a quick and easy calculation to make, it has been applied to a huge variety of metallic (Fe- Ni-S-C-O-Si systems) and silicate liquids (SiO2, Al2O3, TiO2, FeO, MgO, CaO, Na2O, K2O systems). This approach has surely led to values that have little meaning, yet are applied with great confidence, for example, to a terrestrial mantle at "IW-2". Although fO2 can be circumvented in some cases by consideration of Fe-M distribution coefficient, these do not eliminate the effects of alloy or silicate liquid compositional variation, or the specific chemical effects of S in the silicate liquid, for example. In order to address the issue of what the actual value of fO2 is in any given experiment, we have calculated fO2 from the equilibria 2Fe (metal) + SiO2 (liq) + O2 = Fe2SiO4 (liq).

Righter, K.↗

SPENVIS Implementation of End-of-Life Solar Cell Calculations Using the Displacement Damage Dose Methodology

This paper presents a method for using the SPENVIS on-line computational suite to implement the displacement damage dose (D(sub d)) methodology for calculating end-of-life (EOL) solar cell performance for a specific space mission. This paper builds on our previous work that has validated the D(sub d) methodology against both measured space data [1,2] and calculations performed using the equivalent fluence methodology developed by NASA JPL [3]. For several years, the space solar community has considered general implementation of the D(sub d) method, but no computer program exists to enable this implementation. In a collaborative effort, NRL, NASA and OAI have produced the Solar Array Verification and Analysis Tool (SAVANT) under NASA funding, but this program has not progressed beyond the beta-stage [4]. The SPENVIS suite with the Multi Layered Shielding Simulation Software (MULASSIS) contains all of the necessary components to implement the Dd methodology in a format complementary to that of SAVANT [5]. NRL is currently working with ESA and BIRA to include the Dd method of solar cell EOL calculations as an integral part of SPENVIS. This paper describes how this can be accomplished.

Walters, Robert↗

Using FLUKA to Calculate Spacecraft: Single Event Environments: A Practical Approach

The FLUKA nuclear transport and reaction code can be developed into a practical tool for calculation of spacecraft and planetary surface asset SEE and TID environments. Nuclear reactions and secondary particle shower effects can be estimated with acceptable accuracy both in-flight and in test. More detailed electronic device and/or spacecraft geometries than are reported here are possible using standard FLUKA geometry utilities. Spacecraft structure and shielding mass. Effects of high Z elements in microelectronic structure as reported previously. Median shielding mass in a generic slab or concentric sphere target geometry are at least approximately applicable to more complex spacecraft shapes. Need the spacecraft shielding mass distribution function applicable to the microelectronic system of interest. SEE environment effects can be calculated for a wide range of spacecraft and microelectronic materials with complete nuclear physics. Evaluate benefits of low Z shielding mass can be evaluated relative to aluminum. Evaluate effects of high Z elements as constituents of microelectronic devices. The principal limitation on the accuracy of the FLUKA based method reported here are found in the limited accuracy and incomplete character of affordable heavy ion test data. To support accurate rate estimates with any calculation method, the aspect ratio of the sensitive volume(s) and the dependence must be better characterized.

Koontz, Steve↗

Absorbed Dose and Dose Equivalent Calculations for Modeling Effective Dose

While in orbit, Astronauts are exposed to a much higher dose of ionizing radiation than when on the ground. It is important to model how shielding designs on spacecraft reduce radiation effective dose pre-flight, and determine whether or not a danger to humans is presented. However, in order to calculate effective dose, dose equivalent calculations are needed. Dose equivalent takes into account an absorbed dose of radiation and the biological effectiveness of ionizing radiation. This is important in preventing long-term, stochastic radiation effects in humans spending time in space. Monte carlo simulations run with the particle transport code FLUKA, give absorbed and equivalent dose data for relevant shielding. The shielding geometry used in the dose calculations is a layered slab design, consisting of aluminum, polyethylene, and water. Water is used to simulate the soft tissues that compose the human body. The results obtained will provide information on how the shielding performs with many thicknesses of each material in the slab. This allows them to be directly applicable to modern spacecraft shielding geometries.

Welton, Andrew↗

Initial Assessment of a Rapid Method of Calculating CEV Environmental Heating

An innovative method for rapidly calculating spacecraft environmental absorbed heats in planetary orbit is described. The method employs reading a database of pre-calculated orbital absorbed heats and adjusting those heats for desired orbit parameters. The approach differs from traditional Monte Carlo methods that are orbit based with a planet centered coordinate system. The database is based on a spacecraft centered coordinated system where the range of all possible sun and planet look angles are evaluated. In an example case 37,044 orbit configurations were analyzed for average orbital heats on selected spacecraft surfaces. Calculation time was under 2 minutes while a comparable Monte Carlo evaluation would have taken an estimated 26 hours

Pickney, John T.↗

Comparison of Iterative and Non-Iterative Strain-Gage Balance Load Calculation Methods

The accuracy of iterative and non-iterative strain-gage balance load calculation methods was compared using data from the calibration of a force balance. Two iterative and one non-iterative method were investigated. In addition, transformations were applied to balance loads in order to process the calibration data in both direct read and force balance format. NASA's regression model optimization tool BALFIT was used to generate optimized regression models of the calibration data for each of the three load calculation methods. This approach made sure that the selected regression models met strict statistical quality requirements. The comparison of the standard deviation of the load residuals showed that the first iterative method may be applied to data in both the direct read and force balance format. The second iterative method, on the other hand, implicitly assumes that the primary gage sensitivities of all balance gages exist. Therefore, the second iterative method only works if the given balance data is processed in force balance format. The calibration data set was also processed using the non-iterative method. Standard deviations of the load residuals for the three load calculation methods were compared. Overall, the standard deviations show very good agreement. The load prediction accuracies of the three methods appear to be compatible as long as regression models used to analyze the calibration data meet strict statistical quality requirements. Recent improvements of the regression model optimization tool BALFIT are also discussed in the paper.

Ulbrich, N.↗

Program Calculates Forces in Bolted Structural Joints

FORTRAN 77 computer program calculates forces in bolts in the joints of structures. This program is used in conjunction with the NASTRAN finite-element structural-analysis program. A mathematical model of a structure is first created by approximating its load-bearing members with representative finite elements, then NASTRAN calculates the forces and moments that each finite element contributes to grid points located throughout the structure. The user selects the finite elements that correspond to structural members that contribute loads to the joints of interest, and identifies the grid point nearest to each such joint. This program reads the pertinent NASTRAN output, combines the forces and moments from the contributing elements to determine the resultant force and moment acting at each proximate grid point, then transforms the forces and moments from these grid points to the centroids of the affected joints. Then the program uses these joint loads to obtain the axial and shear forces in the individual bolts. The program identifies which bolts bear the greatest axial and/or shear loads. The program also performs a fail-safe analysis in which the foregoing calculations are repeated for a sequence of cases in which each fastener, in turn, is assumed not to transmit an axial force.

Buder, Daniel A.↗

Bulk Enthalpy Calculations in the Arc Jet Facility at NASA ARC

The Arc Jet Facilities at NASA Ames Research Center generate test streams with enthalpies ranging from 5 MJ/kg to 25 MJ/kg. The present work describes a rigorous method, based on equilibrium thermodynamics, for calculating the bulk enthalpy of the flow produced in two of these facilities. The motivation for this work is to determine a dimensionally-correct formula for calculating the bulk enthalpy that is at least as accurate as the conventional formulas that are currently used. Unlike previous methods, the new method accounts for the amount of argon that is present in the flow. Comparisons are made with bulk enthalpies computed from an energy balance method. An analysis of primary facility operating parameters and their associated uncertainties is presented in order to further validate the enthalpy calculations reported herein.

Thompson, Corinna S.↗