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 1,567 records · Page 87

A Poisson equation formulation for pressure calculations in penalty finite element models for viscous incompressible flows

The calculation of pressures when the penalty-function approximation is used in finite-element solutions of laminar incompressible flows is addressed. A Poisson equation for the pressure is formulated that involves third derivatives of the velocity field. The second derivatives appearing in the weak formulation of the Poisson equation are calculated from the C0 velocity approximation using a least-squares method. The present scheme is shown to be efficient, free of spurious oscillations, and accurate. Examples of applications are given and compared with results obtained using mixed formulations.

Sohn, J. L.↗

Calculation Of Pneumatic Attenuation In Pressure Sensors

Errors caused by attenuation of air-pressure waves in narrow tubes calculated by method based on fundamental equations of flow. Changes in ambient pressure transmitted along narrow tube to sensor. Attenuation of high-frequency components of pressure wave calculated from wave equation derived from Navier-Stokes equations of viscous flow in tube. Developed to understand and compensate for frictional attenuation in narrow tubes used to connect aircraft pressure sensors with pressure taps on affected surfaces.

Whitmore, Stephen A.↗

Calculation Of Aeroelastic Transients Using Euler Equations

Method for calculation of transient aeroelastic effects on airplane wings based on Euler equations of flow. Conforming grids facilitate calculations of flows about changing shapes. Capability to compute aeroelastic effects reduces cost of developing aircraft.

Guruswamy, Guru P.↗

Calculation of rotor impedance for use in design analysis of helicopter airframe vibrations

Excessive vibration is one of the most prevalent technical obstacles encountered in the development of new rotorcraft. The inability to predict these vibrations is primarily due to deficiencies in analysis and simulation tools. The Langley Rotorcraft Structural Dynamics Program was instituted in 1984 to meet long term industry needs in the area of rotorcraft vibration prediction. As a part of the Langley program, this research endeavors to develop an efficient means of coupling the rotor to the airframe for preliminary design analysis of helicopter airframe vibrations. The main effort was to modify the existing computer program for modeling the dynamic and aerodynamic behavior of rotorcraft called DYSCO (DYnamic System COupler) to calculate the rotor impedance. DYSCO was recently developed for the U.S. Army and has proven to be adaptable for the inclusion of new solution methods. The solution procedure developed to use DYSCO for the calculation of rotor impedance is presented. Verification of the procedure by comparison with a known solution for a simple wind turbine model is about 75 percent completed, and initial results are encouraging. After the wind turbine impedance is confirmed, the verification effort will continue by comparison to solutions of a more sophisticated rotorcraft model. Future work includes determination of the sensitivity of the rotorcraft airframe vibrations to helicopter flight conditions and rotor modeling assumptions. When completed, this research will ascertain the feasibility and efficiency of the impedance matching method of rotor-airframe coupling for use in the analysis of airframe vibrations during the preliminary rotorcraft design process.

Nygren, Kip P.↗

Further Improvements to Nozzle Boundary Layer Calculations in BLIMPJ

Further improvements made to advance the current Boundary Layer Integral Matrix Procedure - Version J (BLIMPJ) containing previously modeled simplified calculation methods by accounting for condensed phase, thick boundary layer and free stream turbulence effects are discussed. The condensed phase effects were included through species composition effect considered via input to the code and through particle damping effect considered via a turbulence model. The thrust loss calculation procedure for thick boundary layer effects was improved and the optimization of net thrust with respect to nozzle length was performed. The effects of free stream turbulence were approximately modeled in the turbulence model.

Praharaj, S. C.↗

Calculations of separated 3-D flows with a pressure-staggered Navier-Stokes equations solver

A Navier-Stokes equations solver based on a pressure correction method with a pressure-staggered mesh and calculations of separated three-dimensional flows are presented. It is shown that the velocity pressure decoupling, which occurs when various pressure correction algorithms are used for pressure-staggered meshes, is caused by the ill-conditioned discrete pressure correction equation. The use of a partial differential equation for the incremental pressure eliminates the velocity pressure decoupling mechanism by itself and yields accurate numerical results. Example flows considered are a three-dimensional lid driven cavity flow and a laminar flow through a 90 degree bend square duct. For the lid driven cavity flow, the present numerical results compare more favorably with the measured data than those obtained using a formally third order accurate quadratic upwind interpolation scheme. For the curved duct flow, the present numerical method yields a grid independent solution with a very small number of grid points. The calculated velocity profiles are in good agreement with the measured data.

Kim, S.-W.↗

Transonic cascade flow calculations using non-periodic C-type grids

A new kind of C-type grid is proposed for turbomachinery flow calculations. This grid is nonperiodic on the wake and results in minimum skewness for cascades with high turning and large camber. Euler and Reynolds averaged Navier-Stokes equations are discretized on this type of grid using a finite volume approach. The Baldwin-Lomax eddy-viscosity model is used for turbulence closure. Jameson's explicit Runge-Kutta scheme is adopted for the integration in time, and computational efficiency is achieved through accelerating strategies such as multigriding and residual smoothing. A detailed numerical study was performed for a turbine rotor and for a vane. A grid dependence analysis is presented and the effect of artificial dissipation is also investigated. Comparison of calculations with experiments clearly demonstrates the advantage of the proposed grid.

Arnone, Andrea↗

Numerical two-dimensional calculations of the formation of the solar nebula

Numerical two dimensional calculations of the formation of the solar nebula are presented. The following subject areas are covered: (1) observational constraints of the properties of the initial solar nebula; (2) the physical problem; (3) review if two dimensional calculations of the formation phase; (4) recent models with hydrodynamics and radiative transport; and (5) further evolution of the system.

Bodenheimer, Peter H.↗

Calculation of a circular jet in crossflow with a multiple-time-scale turbulence model

Numerical calculation of a three dimensional turbulent flow of a jet in a crossflow using a multiple time scale turbulence model is presented. The turbulence in the forward region of the jet is in a stronger inequilibrium state than that in the wake region of the jet, while the turbulence level in the wake region is higher than that in the front region. The calculated flow and the concentration fields are in very good agreement with the measured data, and it indicated that the turbulent transport of mass, concentration, and momentum is strongly governed by the inequilibrium turbulence. The capability of the multiple time scale turbulence model to resolve the inequilibrium turbulence field is also discussed.

Kim, S.-W.↗

An approximate viscous shock layer technique for calculating chemically reacting hypersonic flows about blunt-nosed bodies

An approximate axisymmetric method was developed which can reliably calculate fully viscous hypersonic flows over blunt nosed bodies. By substituting Maslen's second order pressure expression for the normal momentum equation, a simplified form of the viscous shock layer (VSL) equations is obtained. This approach can solve both the subsonic and supersonic regions of the shock layer without a starting solution for the shock shape. The approach is applicable to perfect gas, equilibrium, and nonequilibrium flowfields. Since the method is fully viscous, the problems associated with a boundary layer solution with an inviscid layer solution are avoided. This procedure is significantly faster than the parabolized Navier-Stokes (PNS) or VSL solvers and would be useful in a preliminary design environment. Problems associated with a previously developed approximate VSL technique are addressed before extending the method to nonequilibrium calculations. Perfect gas (laminar and turbulent), equilibrium, and nonequilibrium solutions were generated for airflows over several analytic body shapes. Surface heat transfer, skin friction, and pressure predictions are comparable to VSL results. In addition, computed heating rates are in good agreement with experimental data. The present technique generates its own shock shape as part of its solution, and therefore could be used to provide more accurate initial shock shapes for higher order procedures which require starting solutions.

Cheatwood, F. Mcneil↗

Calculation of unsteady linearized Euler flows in cascades using harmonically deforming grids

A method for calculating unsteady, inviscid, compressible flows in cascades is presented. Using the linearized Euler technique, the flow is decomposed into a steady or mean flow plus a harmonically varying small disturbance flow. The equations that describe the small disturbance flow are linear variable coefficient equations, and are solved using a pseudo-time time marching Lax-Wendroff technique. Unlike previous linearized methods, however, the solution is computed on a harmonically deforming computational grid that conforms to the motion of the vibrating airfoils. The mean flow and perturbation flow solutions are defined in the deforming coordinate system rather than in a coordinate system fixed in space. Hence, no extrapolation terms are required to implement the upwash boundary conditions at the airfoil surfaces, significantly improving the accuracy of the method. For transonic flow calculations, unsteady shock motions are modelled using shock capturing. The unsteady loads due to the shock motion are then seen as pressure impulses. Representative computational results are presented for transonic channel flows and subsonic and transonic cascade flows.

Hall, Kenneth C.↗

Influence of channel base current and varying return stroke speed on the calculated fields of three important return stroke models

Compared here are the calculated fields of the Traveling Current Source (TCS), Modified Transmission Line (MTL), and the Diendorfer-Uman (DU) models with a channel base current assumed in Nucci et al. on the one hand and with the channel base current assumed in Diendorfer and Uman on the other hand. The characteristics of the field wave shapes are shown to be very sensitive to the channel base current, especially the field zero crossing at 100 km for the TCS and DU models, and the magnetic hump after the initial peak at close range for the TCS models. Also, the DU model is theoretically extended to include any arbitrarily varying return stroke speed with height. A brief discussion is presented on the effects of an exponentially decreasing speed with height on the calculated fields for the TCS, MTL, and DU models.

Thottappillil, Rajeev↗

A rapid method of calculating the orbital radiation environment

This paper describes the process of integrating available environment models into a single package that calculates the radiation environment for any earth orbit quickly and accurately with the only input being initial orbit parameters, launch date and length of mission, and units of output. The radiation dose is calculated for particles penetrating varying aluminum shielding thicknesses and incident upon several materials inside. Included also are modifications that have been made to the codes that account for the decrease of the magnetic field with time, which allows the use of data taken in the past to accurately predict a future environment. The complete package has relatively short run times, 20 minutes per mission on a Micro VAX II, allowing multiple iterations for application to mission design and planning.

Gates, Michele M.↗

The relationship between plate velocity and trench viscosity in Newtonian and power-law subduction calculations

The relationship between oceanic trench viscosity and oceanic plate velocity is studied using a Newtonian rheology by varying the viscosity at the trench. The plate velocity is a function of the trench viscosity for fixed Rayleigh number and plate/slab viscosity. Slab velocities for non-Newtonian rheology calculations are significantly different from slab velocities from Newtonian rheology calculations at the same effective Rayleigh number. Both models give reasonable strain rates for the slab when compared with estimates of seismic strain rate. Non-Newtonian rheology eliminates the need for imposed weak zones and provides a self-consistent fluid dynamical mechanism for subduction in numerical convection models.

King, Scott D.↗

DSMC calculations for the delta wing

Results are reported from three-dimensional direct simulation Monte Carlo (DSMC) computations, using a variable-hard-sphere molecular model, of hypersonic flow on a delta wing. The body-fitted grid is made up of deformed hexahedral cells divided into six tetrahedral subcells with well defined triangular faces; the simulation is carried out for 9000 time steps using 150,000 molecules. The uniform freestream conditions include M = 20.2, T = 13.32 K, rho = 0.00001729 kg/cu m, and T(wall) = 620 K, corresponding to lambda = 0.00153 m and Re = 14,000. The results are presented in graphs and briefly discussed. It is found that, as the flow expands supersonically around the leading edge, an attached leeside flow develops around the wing, and the near-surface density distribution has a maximum downstream from the stagnation point. Coefficients calculated include C(H) = 0.067, C(DP) = 0.178, C(DF) = 0.110, C(L) = 0.714, and C(D) = 1.089. The calculations required 56 h of CPU time on the NASA Langley Voyager CRAY-2 supercomputer.

Celenligil, M. Cevdet↗

Benchmark calculations of thermal reaction rates. I - Quantal scattering theory

The thermal rate coefficient for the prototype reaction H + H2 yields H2 + H with zero total angular momentum is calculated by summing, averaging, and numerically integrating state-to-state reaction probabilities calculated by time-independent quantum-mechanical scattering theory. The results are very carefully converged with respect to all numerical parameters in order to provide high-precision benchmark results for confirming the accuracy of new methods and testing their efficiency.

Chatfield, David C.↗

A comparison of three time-dependent wave packet methods for calculating electron-atom elastic scattering cross sections

Three time-dependent wave packet methods for performing elastic scattering calculations from screened Coulomb potentials are compared. The three methods are the time-dependent amplitude density method (TDADM), a Cayley-transform method (CTM), and the Chebyshev propagation method of Tal-Ezer and Kosloff. Both the TDADM and the CTM are based on a time-dependent integral equation for the wave function. In the first, the time-dependent amplitude density is propagated, while in the other two, the wave function is propagated. As a numerical example, phase shifts and cross sections are calculated using a screened Coulomb, Yukawa type potential over the range 200-1000 eV. It is found that, in most cases, all three methods yield comparable accuracy and are about equally efficient computationally. However for l = 0, where the Coulomb well is not screened by the centrifugal potential, the TDADM requires smaller grid spacings to maintain accuracy.

Judson, Richard S.↗

Calculation of aerosol backscatter from airborne continuous wave focused CO2 Doppler lidar measurements. I - Algorithm description

Since 1981 the Royal Signals and Radar Establishment and the Royal Aircraft Establishment, United Kindom, have made vertical and horizontal sounding measurements of aerosol backscatter coefficients at 10.6 microns, using an airborne continuous-wave-focused CO2 Doppler lidar, the Laser True Airspeed System (LATAS). In this paper, the heterodyne signal from the LATAS detector is spectrally analyzed. Then, in conjunction with aircraft flight parameters, the data are processed in a six-stage computer algorithm: set search window, search for peak signal, test peak signal, measure total signal, calculate signal-to-noise ratio, and calculate backscatter coefficient.

Rothermel, Jeffry↗