SEARCH · Engineering Papers
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.
Lewis number corrected for extremely high airspeeds. Algorithm calculates improved, variable value of Lewis number, factor in equation for heat-transfer coefficient at stagnation point of body in hypersonic flow mach numbers of approximately 25. New algorithm improves accuracy of calculations of heat generated by hot air acting on moving body.
Approximate field values computed rapidly. MODEL computer program developed to calculate electromagnetic field values of large loop antenna at all distances to observation point. Antenna assumed to be in x-y plane with center at origin of coordinate system. Calculates field values in both rectangular and spherical components. Also solves for wave impedance. Written in MicroSoft FORTRAN 77.
Computer program incorporates van der Waals equation and correction tables. FLUID program developed to calculate thermodynamic and transport properties of pure fluids in both liquid and gas phases. Properties calculated by use of simple gas model, empirical corrections, and efficient numerical interpolation scheme. Produces results that agree very well with measured values. Much faster than older, more complex programs developed for same purpose.
Lidar observation at Fukuoka has provided over four years the data of ElChichon aerosols in the stratosphere. Analyses of the data show that an enormous amount of volcanic aerosols has continuously decreased since the beginning of 1983 with significant fluctuations. These fluctuations reveal themselves as a seasonal variation of aerosol content with a maximum in winter-spring and a minimum in summer. The vertical structure of the aerosol layer also shows the seasonal variation. Although the height of a peak around 18 km in the vertical profile of scattering ratio show littel variation, the higher second peak appears frequently from late fall and the lower third peak from late winter to late spring just as two and more tropopauses appear in these periods. The mechanism which causes the seasonal variation will be discussed in terms of the transport by the atmoshperic circulation and the removal through the tropopause gap. Radiation budget in the atmosphere was calculated taking into account the large amount of aerosols observed in the early stages of the El Chichon event. The heating rate of the atmosphere is more than 1 K in the bottom region of the stratosphere even in the nighttime. The possible effect of the volcanic aerosols on the other geophysical phenomena will be discussed using the calculated values of the heating rate.
A Navier-Stokes analysis employing the time-dependent Linearized Block Implicit scheme (LBI) was applied to two-dimensional and three-dimensional transonic turbulent cascade flows. In general, the geometrical configuration of the turbine blade impacts both the grid construction procedure and the implementation of the numerical algorithm. Since modern turbine blades of interest are characterized by very blunt leading edges, rounded trailing edges and high stacking angles, a robust grid construction procedure is required that can accommodate the severe body shape while resolving regions of large flow gradients. A constructive O-type grid generation technique, suitable for cascades with rounded trailing edges, was developed and used to construct the C3X turbine cascade coordinate grid. Two-dimensional calculations were performed employing the Navier-Stokes procedure for the C3X turbine cascade, and the predicted pressure coefficients and heat transfer rates were compared with the experimental data. Three-dimensional Navier-Stokes calculations were also performed.
A computational method is described that includes the effects of static aeroelastic wing deflections in steady transonic aerodynamic calculations. This method, known as the Transonic Aero-elastic Program System (TAPS), interacts a 3D transonic computer code with boundary layer and a linear finite element structural analysis codes to calculate wing pressures and deflections. The nonlinear nature of the transonic flow makes it necessary to couple the aerodynamic and structures codes in an iterative manner. TAPS has been arranged in a modular fashion so that different aerodynamic or structures programs may be used with a minimum of coding changes required. Results obtained using two different aerodynamic codes in TAPS are given, and those results are correlated with experimental data.
A method for the direct calculation of the wall induced interference velocity in two dimensional flow based on Cauchy's integral formula was derived. This one-step method allows the calculation of the residual corrections and the required wall adaptation for interference-free flow starting from the wall pressure distribution without any model representation. Demonstrated applications are given.
The factors affecting optical communication link performance differ substantially from those at microwave frequencies, due to the drastically differing technologies, modulation formats, and effects of quantum noise in optical communications. In addition detailed design control table calculations for optical systems are less well developed than corresponding microwave system techniques, reflecting the relatively less mature state of development of optical communications. Described below are detailed calculations of received optical signal and background power in optical communication systems, with emphasis on analytic models for accurately predicting transmitter and receiver system losses.
An energy-dependent, perturbation expansion solution for heavy ion transport in one dimension is used to perform depth-dose calculations for 670/MeV nucleon Ne-20 beams incident upon a thick water target. Comparisons of predictions obtained by using typical energy-independent approximations and those obtained with fully energy-dependent input parameters are made. It is found that the calculated doses are underestimated when the energy-independent input approximations are used. The major source of error, however, is the lack of charge and mass conservation in the Silberberg-Tsao fragmentation parameters.
An interactive inviscid core flow-boundary layer method is presented for the calculation of turbomachine channel flows. For this method, a one-dimensional inviscid core flow is assumed. The end-wall and blade surface boundary layers are calculated using an integral entrainment method. The boundary layers are assumed to be collateral and thus are two-dimensional. The boundary layer equations are written in a streamline coordinate system. The streamwise velocity profiles are approximated by power law profiles. Compressibility is accounted for in the streamwise direction but not in the normal direction. Equations are derived for the special cases of conical and two-dimensional rectangular diffusers. For these cases, the assumptions of a one-dimensional core flow and collateral boundary layers are valid. Results using the method are compared with experiment and good quantitative agreement is obtained.
Correction of airfoil data for sidewall boundary-layer effects requires a knowledge of the boundary-layer displacement thickness and the shape factor with the tunnel empty. To facilitate calculation of these quantities under various test conditions for the Langley 0.3 m Transonic Cryogenic Tunnel, a computer program was written. This program reads the various tunnel parameters and the boundary-layer rake total head pressure measurements directly from the Engineering Unit tapes to calculate the required sidewall boundary-layer parameters. Details of the method along with the results for a sample case are presented.
A computer program written to calculate the static pressure position error of airspeed systems contains five separate methods for determining position error, of which the user may select from one to five at a time. The program uses data from both the test aircraft and the ground-based radar to calculate the error. In addition, some of the methods require rawinsonde data or an atmospheric analysis, or both. The program output lists the corrections to Mach number, altitude, and static pressure that are due to position error. Reference values such as angle of attack, angle of sideslip, indicated Mach number, indicated pressure altitude, stagnation pressure, and total temperature are also listed.
The effective medium approximation (EMA) formalism developed for scalar wave calculations in solid state physics is generalized to electromagnetic wave scattering in a dense random medium. Results are applied to compute the effective propagation constant in a dense medium involving discrete spherical scatterers. When compared with a common quasicrystalline approximation (QCA), it is found that EMA accounts for backward scattering and the effect of correlation among three scatterers which are not available in QCA. It is also found that there is not much difference in the calculated normalized phase velocity between the use of these two approximations. However, there is a significant difference in the computed effective loss tangent in a nonabsorptive random medium. The computed effective loss tangent using EMA and measurements from a snow medium are compared, showing good agreement.
A computer program based on the Quasi-Vortex-Lattice Method of Lan is presented for calculating longitudinal and lateral-directional aerodynamic characteristics of nonplanar wing-body combination. The method is based on the assumption of inviscid subsonic flow. Both attached and vortex-separated flows are treated. For the vortex-separated flow, the calculation is based on the method of suction analogy. The effect of vortex breakdown is accounted for by an empirical method. A summary of the theoretical method, program capabilities, input format, output variables and program job control set-up are described. Three test cases are presented as guides for potential users of the code.
The aim of this study is to find a reliable numerical algorithm to calculate thermal design sensitivities of a transient problem with discontinuous derivatives. The thermal system of interest is a transient heat conduction problem related to the curing process of a composite laminate. A logical function which can smoothly approximate the discontinuity is introduced to modify the system equation. Two commonly used methods, the adjoint variable method and the direct differentiation method, are then applied to find the design derivatives of the modified system. The comparisons of numerical results obtained by these two methods demonstrate that the direct differentiation method is a better choice to be used in calculating thermal design sensitivity.
An amalgamation of two published works that discuss the observation and theoretical calculations of OV (T approx. 250,000K) and Hard X-rays (30 to 100keV) emitted during flares are presented. The papers are by Poland et al (1984) and Mariska and Poland (1985). The observations of Hard X-rays and OV show that the excitation processes for each type of emission are closely coupled. Except for small differences the two types of emission rise and fall together during a flare. Model calculations are able to reproduce this behavior to a large extent, only when conductive processes do not dominate the energy transport processes.
The problem concerned in this work is that of calculating magnetic field configurations in which the Lorentz force (vector)j x (vector)B is everywhere zero, subject to specified boundary conditions. The magnetic field is represented in terms of Clebsch variables in the form (vector)B = del alpha x del beta. These variables are constant on any field line. The most appropriate choice of boundary conditions is to specify the values of alpha and beta on the bounding surface. It is proposed that the field lines move in the direction of local Lorentz force and relax towards a force-free field configuration. This concept leads to an iteration procedure for modifying the variables alpha and beta that tends asymptotically towards the force-free state. This method is first applied to a simple problem in two rectangular dimensions; the calculation shows that the convergence of magnetic field energy to a minimum state (force-free) is close to exponential. This method is then applied to study some astrophysical force-free magnetic fields, such as the structures and evolution of magnetic fields of rotating sunspots and accretion disks. The implication of the results, as related to the mechanisms of solar flares, extragalactic radio sources and radio jets, are discussed.
A method is described for calculating unsteady transonic flow with viscous interaction by coupling a steady integral boundary-layer code with an unsteady, transonic, inviscid small-disturbance computer code in a quasi-steady fashion. Explicit coupling of the equations together with viscous -inviscid iterations at each time step yield converged solutions with computer times about double those required to obtain inviscid solutions. The accuracy and range of applicability of the method are investigated by applying it to four AGARD standard airfoils. The first-harmonic components of both the unsteady pressure distributions and the lift and moment coefficients have been calculated. Comparisons with inviscid calcualtions and experimental data are presented. The results demonstrate that accurate solutions for transonic flows with viscous effects can be obtained for flows involving moderate-strength shock waves.