Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “characteristic equations”

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 865 records · Page 48

Metal-insulator-semi-conductor studies of lead telluride

The capacitance and conductance-voltage characteristics were measured on metal-insulator-semiconductor capacitors fabricated with zirconium dioxide films on single-crystal lead telluride. At 77 K, on both n- and p-type substrates, evidence of surface potential control was obtained. Comparison of the measured capacitance-voltage characteristics with those calculated from the equilibrium solution of the one-dimensional Poisson equation indicated qualitative agreement, although the slope of the measured capacitance in the region near the capacitance minimum was less steep than calculated.

Lilly, D. A.↗

Mercury's helium exosphere after Mariner 10's third encounter

From a comparison of the Mariner 10 third encounter UV spectrometer data with intensities generated from a newly constructed model exosphere, a new value of 4.5 x 10 to the -4th power for the fraction of the solar wind He(++) flux to be intercepted and captured by Mercury's magnetosphere, if the observed He atmosphere is maintained by the solar wind, has been derived. If an internal source for He prevails, the corresponding upper bound for the global outgassing rate is estimated to be 4.5 x 10 to the 22nd power/s. These values differ from those given earlier because of the present use of a surface temperature distribution satisfying the heat equation over Mercury's entire surface which employs Mariner-10-determined mean surface thermal characteristics. The mean standoff distance of Mercury's magnetopause averaged over Mercury's orbit is also used. Agreement between the observed and calculated intensities is found to be good. Minor discrepancies on the nightside of the terminator are explicable in terms of differences between actual and computed temperatures and scale height structure changes.

Curtis, S. A.↗

IMP: Interactive mass properties program. Volume 1: Program description

A method of computing a weights and center of gravity analysis of a flight vehicle using interactive graphical capabilities of the Adage 340 computer is described. The equations used to calculate area, volume, and mass properties are based on elemental surface characteristics. The input/output methods employ the graphic support of the Adage computer. Several interactive program options are available for analyzing the mass properties of a vehicle. These options are explained.

Stewart, W. A.↗

Finite element analysis of transonic flows in highly staggered cascades

Finite element method is employed in the solution of full-potential equation for analyzing steady, transonic flows through a series of highly staggered cascades. Convergence characteristics and efficiency of the developed numerical procedure are discussed. Accuracy of the mathematical model is investigated in terms of the employed computational grid. In particular, the importance of accurate representation of the trailing edge together with the shock on the obtained flow results are studied. Numerical results include choked flows, for a cascade of NACA 0012 airfoils with pitch/chord length = 1, and stagger angle of 45 deg.

Akay, H. U.↗

One-dimensional unsteady modeling of supersonic inlet unstart/restart

A quasi-one-dimensional unsteady inviscid analysis of mixed-compression supersonic inlet flow is presented with emphasis on modeling of inlet unstart/restart phenomena. Numerical solution of the governing equations of motion is performed using a computationally efficient shock-capturing split-characteristics algorithm. Inlet unstart is modeled using a mass balance method which relates the expelled normal shock position ahead of the inlet cowl to the amount of spilled mass flow over the inlet housing. Comparison of computed results with experimental data for an axisymmetric inlet at a free-stream Mach number of 2.50 shows quite reasonable agreement over an entire unstart/restart transient which includes centerbody translation and retraction as well as bypass mass flow variations.

Adams, J. C., Jr.↗

Electromagnetic wave scattering in a two-layer anisotropic random medium

For electromagnetic wave propagation and scattering in an anisotropic random medium, the Dyson equation for the mean field and the Bethe-Salpeter equation for the correlation or the covariance of the field were derived. With the random permittivity expressed in a general anisotropic form, the bilocal and the nonlinear approximations are employed to solve the Dyson equation, and the ladder approximation to solve the Bethe-Salpeter equation. The mean dyadic Green's function for a two-layer anisotropic random medium with arbitrary three-dimensional correlation functions has been investigated with the zeroth-order solutions to the Dyson equation under the nonlinear approximation. The effective propagation constants are calculated for the four characteristic waves associated with the coherent vector fields, propagating in an anisotropic random-medium layer, which are the ordinary and extraordinary waves with upward- and downward-propagating vectors.

Lee, J. K.↗

Acoustic propagation using computational fluid dynamics

The propagation characteristics of several helicopter airfoil profiles have been investigated using the transonic small disturbance equation. A test case was performed to generate a moving shock that propagated off the airfoil. Various grids were then examined to determine their ability to accurately capture these propagating shock waves. Finally, the case of airfoil-vortex interactions was thoroughly studied over a wide range of Mach numbers and airfoil shapes with particular emphasis on the transonic regime; this results in a highly conplicated fluctuation of lift, drag, and pitching moment. The calculated acoustic intensity levels, along with the details of the computational flow field, provide new insights into the understanding of transonic airfoil-vortex interactions.

Baeder, J. D.↗

Numerical optimization of conical flow waveriders including detailed viscous effects

A family of optimized hypersonic waveriders is generated and studied wherein detailed viscous effects are included within the optimization process itself. This is in contrast to previous optimized waverider work, wherein purely inviscid flow is used to obtain the waverider shapes. For the present waveriders, the undersurface is a streamsurface of an inviscid conical flowfield, the upper surface is a streamsurface of the inviscid flow over a tapered cylinder (calculated by the axisymmetric method of characteristics), and the viscous effects are treated by integral solutions of the boundary layer equations. Transition from laminar to turbulent flow is included within the viscous calculations. The optimization is carried out using a nonlinear simplex method. The resulting family of viscous hypersonic waveriders yields predicted high values of lift/drag, high enough to break the L/D barrier based on experience with other hypersonic configurations. Moreover, the numerical optimization process for the viscous waveriders results in distinctly different shapes compared to previous work with inviscid-designed waveriders. Also, the fine details of the viscous solution, such as how the shear stress is distributed over the surface, and the location of transition, are crucial to the details of the resulting waverider geometry. Finally, the moment coefficient variations and heat transfer distributions associated with the viscous optimized waveriders are studied.

Bowcutt, Kevin G.↗

Distributed parameter estimation for large beam-like structures for general inputs

In this paper, maximum likelihood estimation for distributed parameter models of large flexible structures has been formulated. Distributed parameter models involve far fewer unknown parameters than independent modal characteristics of finite element models. The closed-form solutions for the partial differential equations with general forced inputs have been derived. The closed-form expressions of the sensitivity functions lead to highly efficient algorithms for analyzing ground or on-orbit test results. Numerical simulations with impulse and step inputs to the NASA Mini-MAST trust have been demonstrated. The estimations of modal properties involve its longitudinal elongation modes, lateral bending modes and torsional modes. The results show that distributed parameter models are promising in the parameter estimation of large flexible structures.

Shen, Ji-Yao↗

Toward the routine aerodynamic analysis of complex configurations

An evaluation is made of the ability of the TranAir CFD code to routinely compute the aerodynamic characteristics of complex subsonic and supersonic aircraft configurations. TranAir solves the full-potential equation for transonic flow about completely arbitrary geometries, using the surface-paneling PanAir technique in geometry definition. The uniform global grid may be locally refined in regions where flow properties are rapidly changing, such as regions where shocks arise, and around wing leading edges. Unlike panel method codes, TranAir solutions are not undermined by small-perturbation assumptions. Illustrative results are presented for such configurations as the F-16A with wingtip-mounted missiles and underwing fuel tanks, a generic fighter configuration, and a model of NASA-Ames' 12-ft Pressure Wind Tunnel.

Madson, Michael D.↗

Perturbation of a normal shock - A test case for unsteady Euler codes

Shock tracking capabilities of modern nonoscillatory schemes for solving the Euler equations are studied by comparing their solutions to solutions obtained using the method of characteristics on two test problems. Both test problems involve accelerating and decelerating shock motions. The schemes tested include first order upwind, second order SONIC-A, third order ENO, and third order TONIC-A. Results show that the second and higher order schemes give nearly identical shock trajectories that are closer to the characteristic solution than the first order upwind scheme. All schemes approach the characteristic solution upon mesh refinement.

Suresh, Ambady↗

Coronal shock waves

The formation of magnetohydrodynamic (MHD) shocks near the leading edge of Coronal Mass Ejections (CME's) is considered using results from both shock theory and numerical simulations. By using just the shock theory results and a simplified model corona, the following is shown: a slow shock front should be concave upward (away from the solar surface); a configuration containing intermediate shocks should also be concave upward near the center of the CME and concave downward at some distance from the CME centerline; and a fast shock configuration should be concave downward. These results are verified with numerical simulations of the MHD equations for the propagation of shocks through a quiescent coronal streamer. By examining the characteristic wave speeds in the ambient streamer, it is estimated that the various shock configurations should form for CME speeds within the following speed increments: slow shocks, 200 to 300 km/sec; intermediate shocks, 300 to 900 km/sec; and fast shocks, greater than 900 km/sec.

Steinolfson, Richard S.↗

SeaWiFS technical report series. Volume 23: SeaWiFS prelaunch radiometric calibration and spectral characterization

Based on the operating characteristics of the Sea-viewing Wide Field-of-view Sensor (SeaWiFS), calibration equations have been developed that allow conversion of the counts from the radiometer into Earth-existing radiances. These radiances are the geophysical properties the instrument has been designed to measure. SeaWiFS uses bilinear gains to allow high sensitivity measurements of ocean-leaving radiances and low sensitivity measurements of radiances from clouds, which are much brighter than the ocean. The calculation of these bilinear gains is central to the calibration equations. Several other factors within these equations are also included. Among these are the spectral responses of the eight SeaWiFS bands. A band's spectral response includes the ability of the band to isolate a portion of the electromagnetic spectrum and the amount of light that lies outside of that region. The latter is termed out-of-band response. In the calibration procedure, some of the counts from the instrument are produced by radiance in the out-of-band region. The number of those counts for each band is a function of the spectral shape of the source. For the SeaWiFS calibration equations, the out-of-band responses are converted from those for the laboratory source into those for a source with the spectral shape of solar flux. The solar flux, unlike the laboratory calibration, approximates the spectral shape of the Earth-existing radiance from the oceans. This conversion modifies the results from the laboratory radiometric calibration by 1-4 percent, depending on the band. These and other factors in the SeaWiFS calibration equations are presented here, both for users of the SeaWiFS data set and for researchers making ground-based radiance measurements in support of Sea WiFS.

Barnes, Robert A.↗

Long-term velocity enhancements in the solar wind

Throughout most of the last three solar cycles, the Pioneer 10, Pioneer 11, Voyager 2, IMP 8, and Pioneer Venus Orbiter spacecraft have observed long-term enhancements in solar wind velocity. These enhancements are typically on the order of 100-200 km/s, with durations on the order of several months to over a year. They are observed over a range of heliocentric distances that ranges from 0.72 to more than 60 AU, which suggests that they are a characteristic feature throughout the entire heliosphere, at least in the vicinity of the solar equator. They appear to be related to the 'long term velocity shifts' reported by Gazis [1987], but are much more widespread. Since the last solar minimum, they have recurred with the 13-year periodicity reported by Richardson et al [1994], but prior to the last solar minimum there were long intervals were his periodicity was different or absent. We examine and characterize these long-term velocity enhancements and compare them to shorter-term variations in the solar wind such as CMEs, interaction regions, merged interaction regions (MIRs) and global merged interaction regions (GMIRs).

Gazis, P. R.↗

Computational Role of Tunneling in a Programmable Quantum Annealer

Quantum tunneling is a phenomenon in which a quantum state tunnels through energy barriers above the energy of the state itself. Tunneling has been hypothesized as an advantageous physical resource for optimization. Here we present the first experimental evidence of a computational role of multiqubit quantum tunneling in the evolution of a programmable quantum annealer. We developed a theoretical model based on a NIBA Quantum Master Equation to describe the multi-qubit dissipative cotunneling effects under the complex noise characteristics of such quantum devices.We start by considering a computational primitive, the simplest non-convex optimization problem consisting of just one global and one local minimum. The quantum evolutions enable tunneling to the global minimum while the corresponding classical paths are trapped in a false minimum. In our study the non-convex potentials are realized by frustrated networks of qubit clusters with strong intra-cluster coupling. We show that the collective effect of the quantum environment is suppressed in the critical phase during the evolution where quantum tunneling decides the right path to solution. In a later stage dissipation facilitates the multiqubit cotunneling leading to the solution state. The predictions of the model accurately describe the experimental data from the D-WaveII quantum annealer at NASA Ames. In our computational primitive the temperature dependence of the probability of success in the quantum model is opposite to that of the classical paths with thermal hopping. Specially, we provide an analysis of an optimization problem with sixteen qubits,demonstrating eight qubit cotunneling that increases success probabilities. Furthermore, we report results for larger problems with up to 200 qubits that contain the primitive as subproblems.

hard problems↗

Isokinetic TWC Evaporator Probe: Calculations and Systemic Uncertainty Analysis

A new Isokinetic Total Water Content Evaporator (IKP2) was downsized from a prototype instrument, specifically to make airborne measurements of hydrometeor total water content (TWC) in deep tropical convective clouds to assess the new ice crystal Appendix D icing envelope. The probe underwent numerous laboratory and wind tunnel investigations to ensure reliable operation under the difficult high altitude/speed/TWC conditions under which other TWC instruments have been known to either fail, or have unknown performance characteristics and the results are presented in a companion paper. This paper presents the equations used to determine the total water content (TWC) of the sampled atmosphere from the values measured by the IKP2 or necessary ancillary data from other instruments. The uncertainty in the final TWC is determined by propagating the uncertainty in the measured values through the calculations to the final result. Two techniques were used and the results compared. The first is a typical analytical method of propagating uncertainty and the second performs a Monte Carlo simulation. The results are very similar with differences that are insignificant for practical purposes. The uncertainty is between 2 percent and 3 percent at most practical operating conditions. The capture efficiency of the IKP2 was also examined based on a computational fluid dynamic simulation of the original IKP and scaled down to the IKP2. Particles above 24 microns were found to have a capture efficiency greater than 99 percent at all operating conditions.

Aircraft Icing↗

Modeling the Effects of a Backward-Facing Step on Boundary-Layer Transition

We model transition to turbulence in a two-dimensional boundary layer downstream of a backward-facing step (BFS) along a flat plate. With the goal of evaluating the available engineering models for predicting the effects of step excrescences on the transition characteristics, two separate methodologies are used to monitor the streamwise shift in the transition onset location as the step height and the flow speed are varied across the range of a previously reported experiment involving step-height-to-local-displacement-thickness ratios of 0 < h/δ* < 1.6. Unlike the variable N -factor method from the previous literature, both of these methods are general in scope and do not involve any empirical correlations that are specific to step excrescences. The first of these techniques involves an N -factor method that directly accounts for the change in boundary-layer instability characteristics due to the step. Stability computations using the harmonic linearized Navier-Stokes equations (HLNSE), which fully account for the nonparallel-mean-flow effects close to the BFS, indicate that the measured transition locations at nearly all test conditions ( h/δ* < 1.3) correlate well with a computed N -factor of N tr = 7.6, demonstrating a successful stability-based transition criterion related to step excrescences. Linear stability theory, which does not account for nonparallel effects, demonstrates reasonable agreement with the HLNSE results, yielding good predictions for the overall trends, but predicts a somewhat earlier onset of transition than HLNSE. The other methodology used in this work involves transport-equation-based transition models. We first show that the Langtry-Menter y - Re θt transition model cannot accurately predict the location of transition onset for moderate BFS heights because it is unable to accurately account for the flow history effects. Along with the Langtry-Menter transition model, we also show the amplification factor transport model does not produce accurate transition locations for subsonic flow over steps even though it accounts for some flow history effects.

Transition↗

Modeling the Effects of a Backward-Facing Step on Boundary-Layer Transition

We model transition to turbulence in a two-dimensional boundary layer downstream of a backward-facing step (BFS) along a flat plate. With the goal of evaluating the available engineering models for predicting the effects of step excrescences on the transition characteristics, two separate methodologies are used to monitor the streamwise shift in the transition onset location as the step height and the flow speed are varied across the range of a previously reported experiment involving step-height-to-local-displacement-thickness ratios of 0 < h/δ* < 1.6. Unlike the variable N -factor method from the previous literature, both of these methods are general in scope and do not involve any empirical correlations that are specific to step excrescences. The first of these techniques involves an N -factor method that directly accounts for the change in boundary-layer instability characteristics due to the step. Stability computations using the harmonic linearized Navier-Stokes equations (HLNSE), which fully account for the nonparallel-mean-flow effects close to the BFS, indicate that the measured transition locations at nearly all test conditions ( h/δ* < 1.3) correlate well with a computed N -factor of N tr = 7.6, demonstrating a successful stability-based transition criterion related to step excrescences. Linear stability theory, which does not account for nonparallel effects, demonstrates reasonable agreement with the HLNSE results, yielding good predictions for the overall trends, but predicts a somewhat earlier onset of transition than HLNSE. The other methodology used in this work involves transport-equation-based transition models. We first show that the Langtry-Menter y - Re θt transition model cannot accurately predict the location of transition onset for moderate BFS heights because it is unable to accurately account for the flow history effects. Along with the Langtry-Menter transition model, we also show the amplification factor transport model does not produce accurate transition locations for subsonic flow over steps even though it accounts for some flow history effects.

Transition↗