Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “calculator”

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

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

All-electron molecular Dirac-Hartree-Fock calculations - Properties of the XH4 and XH2 molecules and the reaction energy XH4 yields XH2 + H2, X = Si, Ge, Sn, Pb

Relativistic corrections to a number of properties of the Group IV hydrides are calculated using the Dirac-Hartree-Fock method. The use of first-order perturbation theory is sufficient to obtain relativistic corrections for Ge, but the effects of spin-orbit interaction and other higher-order effects begin to show for Sn and become important for Pb. The energy of the reaction XH4 yields XH2 + H2 (X = Si, Ge, Sn, and Pb) is also calculated. The results are compared with relativistic effective core potential calculations, first-order perturbation theory calculations, and limited experimental data.

Dyall, Kenneth G.↗

An efficient pressure-velocity procedure for gas-droplet two-phase flow calculations

This paper presents a noniterative numerical technique for computing time-dependent gas-droplet flows. The method is a fully interacting combination of Eulerian fluid and Lagrangian particle calculations. The interaction calculations between the two phases are formulated on a pressure-velocity-coupling procedure based on the operator-splitting technique. This procedure eliminates the global iterations required in the conventional particle-source-in-cell procedure. Turbulent dispersion calculations are treated by a stochastic procedure. Numerical calculations and comparisons with available experimental data as well as efficiency assessments are given for some sprays typical of spray combustion applications.

Chen, C. P.↗

Ice Accretion and Performance Degradation Calculations with LEWICE/NS

The LEWICE ice accretion computer code has been extended to include the solution of the two-dimensional Navier-Stokes equations. The code is modular and contains separate stand-alone program elements that create a grid, calculate the flow field parameters, calculate the droplet trajectory paths, determine the amount of ice growth, calculate aeroperformance changes, and plot results. The new elements of the code are described. Calculated results are compared to experiment for several cases, including both ice shape and drag rise.

Potapczuk, Mark G.↗

Calculations of turbulent separated flows

A numerical study of incompressible turbulent separated flows is carried out by using two-equation turbulence models of the K-epsilon type. On the basis of realizability analysis, a new formulation of the eddy-viscosity is proposed which ensures the positiveness of turbulent normal stresses - a realizability condition that most existing two-equation turbulence models are unable to satisfy. The present model is applied to calculate two backward-facing step flows. Calculations with the standard K-epsilon model and a recently developed RNG-based K-epsilon model are also made for comparison. The calculations are performed with a finite-volume method. A second-order accurate differencing scheme and sufficiently fine grids are used to ensure the numerical accuracy of solutions. The calculated results are compared with the experimental data for both mean and turbulent quantities. The comparison shows that the present model performs quite well for separated flows.

Zhu, J.↗

Ice accretion and performance degradation calculations with LEWICE/NS

The LEWICE ice accretion computer code has been extended to include the solution of the two-dimensional Navier-Stokes equations. The code is modular and contains separate stand-alone program elements that create a grid, calculate the flow field parameters, calculate the droplet trajectory paths, determine the amount of ice growth, calculate aeroperformance changes, and plot results. The new elements of the code are described. Calculated results are compared to experiment for several cases, including both ice shape and drag rise.

Potapczuk, Mark G.↗

All-electron molecular Dirac-Hartree-Fock calculations - Properties of the group IV monoxides GeO, SnO, and PbO

Dirac-Hartree-Fock calculations have been carried out on the ground states of the group IV monoxides GeO, SnO and PbO. Geometries, dipole moments and infrared data are presented. For comparison, nonrelativistic, first-order perturbation and relativistic effective core potential calculations have also been carried out. Where appropriate the results are compared with the experimental data and previous calculations. Spin-orbit effects are of great importance for PbO, where first-order perturbation theory including only the mass-velocity and Darwin terms is inadequate to predict the relativistic corrections to the properties. The relativistic effective core potential results show a larger deviation from the all-electron values than for the hydrides, and confirm the conclusions drawn on the basis of the hydride calculations.

Dyall, Kenneth G.↗

Chromospheric heating by acoustic shocks - A confrontation of GHRS observations of Alpha Tauri (K5 III) with ab initio calculations

We compare ab initio calculations of semiforbidden C II line profiles near 2325 A with recently published observations of the inactive red giant Alpha Tau (K5 III) obtained using the GHRS on board the Hubble Space Telescope. Our one-dimensional, time-dependent calculations assume that the chromosphere is heated by stochastic acoustic shocks generated by photospheric convection. We calculate various models using results from traditional (mixing length) convection zone calculations as input to hydrodynamical models. The semiforbidden C II line profiles and ratios provide sensitive diagnostics of chromospheric velocity fields, electron densities, and temperatures. We identify major differences between observed and computed line profiles which are related to basic gas dynamics and which are probably not due to technical modeling restrictions. If the GHRS observations are representative of chromospheric conditions at all epochs, then one (or more) of our model assumptions must be incorrect. Several possibilities are examined. We predict time variability of semiforbidden C II lines for comparison with observations. Based upon data from the IUE archives, we argue that photospheric motions associated with supergranulation or global pulsation modes are unimportant in heating the chromosphere of Alpha Tau.

Judge, P. G.↗

Calculating potential fields using microchannel spatial light modulators

We describe and present experimental results of the optical calculation of potential field maps suitable for mobile robot navigation. The optical computation employs two write modes of a microchannel spatial light modulator (MSLM). In one mode, written patterns expand spatially, and this characteristic is used to create an extended two dimensional function representing the influence of the goal in a robot's workspace. Distinct obstacle patterns are written in a second, non-expanding, mode. A model of the mechanisms determining MSLM write mode characteristics is developed and used to derive the optical calculation time for full potential field maps. Field calculations at a few hertz are possible with current technology, and calculation time vs. map size scales favorably in comparison to digital electronic computation.

Reid, Max B.↗

Calculations of the Effect of Wing Twist on the Air Forces Acting on a Monoplane Wing

A method is presented for calculating the aerodynamic forces on a moncylane wing, taking into account the elastic twisting of the wing due to these forces. The lift distribution along the span is calculated by the formulas of Amstutz as a function of the geometrical characteristics of the wing and of the twist at stations 60 and 90 percent of the semispan. The twist for a given lift distribution is calculated by means of influence lines. As a numerical example, the forces on a Swiss military D.2V airplane are calculated. Comparisons with the strip method and with the ordinary stress-analysis method are also given.

FORMULAS-LIFT DISTRIBUTION-AMSTUTZ↗

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↗

Method for calculating the rolling and yawing moments due to rolling for unswept wings with or without flaps or ailerons by use of nonlinear section lift data

The methods of NACA reports 865 and 1090 have been applied to the calculation of the rolling-moment and yawing-moment coefficients due to rolling for unswept wings with or without flaps or ailerons. The methods are based on lifting-line theory and allow the use of nonlinear section lift data. The method presented in this report permits calculations to be made somewhat beyond maximum lift for wings having no twist or continuous twist and employing airfoil sections which do not display large discontinuities in the lift curves. Calculations can be made up to maximum lift for wings with discontinuous twist such as that produced by partial-span flaps or ailerons, or both. Two calculated examples are presented in simplified computing forms in order to illustrate the procedures involve.

Martina, Albert P↗

Calculated spanwise lift distributions, influence functions, and influence coefficients for unswept wings in subsonic flow

Spanwise lift distributions have been calculated for nineteen unswept wings with various aspect ratios and taper ratios and with a variety of angle-of-attack or twist distributions, including flap and aileron deflections, by means of the Weissinger method with eight control points on the semispan. Also calculated were aerodynamic influence coefficients which pertain to a certain definite set of stations along the span, and several methods are presented for calculating aerodynamic influence functions and coefficients for stations other than those stipulated. The information presented in this report can be used in the analysis of untwisted wings or wings with known twist distributions, as well as in aeroelastic calculations involving initially unknown twist distributions.

Diederich, Franklin W↗

Use of Navier-Stokes methods for the calculation of high-speed nozzle flow fields

Flows through three reference nozzles have been calculated to determine the capabilities and limitations of the widely used Navier-Stokes solver, PARC. The nozzles examined have similar dominant flow characteristics as those considered for supersonic transport programs. Flows from an inverted velocity profile (IVP) nozzle, an under expanded nozzle, and an ejector nozzle were examined. PARC calculations were obtained with its standard algebraic turbulence model, Thomas, and the two-equation turbulence model, Chien k-epsilon. The Thomas model was run with the default coefficient of mixing set at both 0.09 and a larger value of 0.13 to improve the mixing prediction. Calculations using the default value substantially underpredicted the mixing for all three flows. The calculations obtained with the higher mixing coefficient better predicted mixing in the IVP and underexpanded nozzle flows but adversely affected PARC's convergence characteristics for the IVP nozzle case. The ejector nozzle case did not converge with the Thomas model and the higher mixing coefficient. The Chien k-epsilon results were in better agreement with the experimental data overall than were those of the Thomas run with the default mixing coefficient, but the default boundary conditions for k and epsilon underestimated the levels of mixing near the nozzle exits.

Georgiadis, Nicholas J.↗

Use of the sampling theorem to speed up near-field physical optics scattering calculations

Physical optics scattering calculations performed on the DSN 34-m beam-waveguide antennas at Ka-band (32 GHz) require approximately 12 hr of central processing unit time on a Cray Y-MP2 computer. This is excessive in terms of resource utilization and turnaround time. Typically, the calculations involve five surfaces, and the calculations are done two surfaces at a time. The sampling theorem is used to reduce the number of current values that must be calculated over the second surface by performing a physical optics integration over the first surface. The additional number of current values required on the second surface by subsequent physical optics integrations is obtained by interpolation over the original current values. Time improvements on the order of a factor of 2 to 4 were obtained for typical scattering pairs.

Cramer, P. W.↗

Theoretical calculations of pressure broadening coefficients for H2O perturbed by hydrogen or helium gas

Halfwidths were calculated for H2O with H2 as a broadening gas and were estimated for He as the broadening species. The calculations used the model of Robert and Bonamy with parabolic trajectories and all relevant terms in the interaction potential. The calculations investigated the dependence of the halfwidth on the order of the atom-atom expansion, the rotational states, and the temperature in the range 200 to 400K. Finally, calculations were performed for many transitions of interest in the 5 micrometer window region of the spectrum. The resulting data will be supplied to Dr. R. Freedman for extracting accurate water mixing ratios from the analysis of the thermal channels for the Net Flux experiment on the Galileo probe.

Gamache, Robert R.↗

Calculation of three-dimensional (3-D) internal flow by means of the velocity-vorticity formulation on a staggered grid

A method has been developed to accurately compute the viscous flow in three-dimensional (3-D) enclosures. This method is the 3-D extension of a two-dimensional (2-D) method developed for the calculation of flow over airfoils. The 2-D method has been tested extensively and has been shown to accurately reproduce experimental results. As in the 2-D method, the 3-D method provides for the non-iterative solution of the incompressible Navier-Stokes equations by means of a fully coupled implicit technique. The solution is calculated on a body fitted computational mesh incorporating a staggered grid methodology. In the staggered grid method, the three components of vorticity are defined at the centers of the computational cell sides, while the velocity components are defined as normal vectors at the centers of the computational cell faces. The staggered grid orientation provides for the accurate definition of the vorticity components at the vorticity locations, the divergence of vorticity at the mesh cell nodes and the conservation of mass at the mesh cell centers. The solution is obtained by utilizing a fractional step solution technique in the three coordinate directions. The boundary conditions for the vorticity and velocity are calculated implicitly as part of the solution. The method provides for the non-iterative solution of the flow field and satisfies the conservation of mass and divergence of vorticity to machine zero at each time step. To test the method, the calculation of simple driven cavity flows have been computed. The driven cavity flow is defined as the flow in an enclosure driven by a moving upper plate at the top of the enclosure. To demonstrate the ability of the method to predict the flow in arbitrary cavities, results will he shown for both cubic and curved cavities.

Stremel, Paul M.↗

Calculated coupling efficiency between an elliptical-core optical fiber and an optical waveguide over temperature

To determine the feasibility of coupling the output of a single-mode optical fiber into a single-mode rib waveguide in a temperature varying environment, a theoretical calculation of the coupling efficiency between the two was investigated. Due to the complex geometry of the rib guide, there is no analytical solution to the wave equation for the guided modes, thus, approximation and/or numerical techniques must be utilized to determine the field patterns of the guide. In this study, three solution methods were used for both the fiber and guide fields; the effective-index method (EIM), Marcatili's approximation, and a Fourier method. These methods were utilized independently to calculate the electric field profile of each component at two temperatures, 20 C and 300 C, representing a nominal and high temperature. Using the electric field profile calculated from each method, the theoretical coupling efficiency between an elliptical-core optical fiber and a rib waveguide was calculated using the overlap integral and the results were compared. It was determined that a high coupling efficiency can be achieved when the two components are aligned. The coupling efficiency was more sensitive to alignment offsets in the y direction than the x, due to the elliptical modal field profile of both components. Changes in the coupling efficiency over temperature were found to be minimal.

Tuma, Margaret L.↗