Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Programming, Linear”

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,639 records · Page 91

Design of a celestial Thomson-scattering X-ray polarimeter.

A general discussion of the scientific importance and status of stellar X-ray polarimetry is presented. A stellar X-ray polarimeter designed to fit into the bottom half of the NASA OSO-1 wheel compartment or other similar spacecraft is described. In this design, the linear polarization is obtained as a function of energy. The sensitivity of the polarimeter in the 4-24 keV energy range was optimized with the aid of a Monte Carlo simulation computer program and is given for several important celestial X-ray sources. Estimates of sensitivity thresholds for a much larger polarimeter, suitable for flight in the NASA High-Energy Astronomy Observatory (HEAO), are also given. The minimum detectable polarization for several X-ray sources is given.

Landecker, P. B.↗

Performance Analysis and Design Synthesis (PADS) computer program. Volume 3: User manual

The two-fold purpose of the Performance Analysis and Design Synthesis (PADS) computer program is discussed. The program can size launch vehicles in conjunction with calculus-of-variations optimal trajectories and can also be used as a general purpose branched trajectory optimization program. For trajectory optimization alone or with sizing, PADS has two trajectory modules. The first trajectory module uses the method of steepest descent. The second module uses the method of quasi-linearization, which requires a starting solution from the first trajectory module.

Source record↗

Analysis of high aspect ratio jet flap wings of arbitrary geometry.

Paper presents a design technique for rapidly computing lift, induced drag, and spanwise loading of unswept jet flap wings of arbitrary thickness, chord, twist, blowing, and jet angle, including discontinuities. Linear theory is used, extending Spence's method for elliptically loaded jet flap wings. Curves for uniformly blown rectangular wings are presented for direct performance estimation. Arbitrary planforms require a simple computer program. Method of reducing wing to equivalent stretched, twisted, unblown planform for hand calculation is also given. Results correlate with limited existing data, and show lifting line theory is reasonable down to aspect ratios of 5.

Lissaman, P. B. S.↗

Transmission of sound through nonuniform circular ducts with compressible mean flows

An acoustic theory is developed to determine the sound transmission and attenuation through an infinite, hard-walled or lined, circular duct carrying compressible, sheared, mean flows and having a variable cross section. The theory is applicable to large as well as small axial variations, as long as the mean flow does not separate. Although the theory is described for circular ducts, it is applicable to other duct configurations - annular, two dimensional, and rectangular. The theory is described for the linear problem, but the technique is general and has the advantage of being applicable to the nonlinear case as well as the linear case. The technique is based on solving for the envelopes of the quasi-parallel acoustic modes that exist in the duct instead of solving for the actual wave. A computer program was developed. The mean flow model consists of a one dimensional flow in the core and a quarter-sine profile in the boundary layer. Results are presented for the reflection and transmission coefficients in ducts with varying slopes and carrying different mean flows.

Nayfeh, A. H.↗

Model verification of large structural systems

A computer program for the application of parameter identification on the structural dynamic models of space shuttle and other large models with hundreds of degrees of freedom is described. Finite element, dynamic, analytic, and modal models are used to represent the structural system. The interface with math models is such that output from any structural analysis program applied to any structural configuration can be used directly. Processed data from either sine-sweep tests or resonant dwell tests are directly usable. The program uses measured modal data to condition the prior analystic model so as to improve the frequency match between model and test. A Bayesian estimator generates an improved analytical model and a linear estimator is used in an iterative fashion on highly nonlinear equations. Mass and stiffness scaling parameters are generated for an improved finite element model, and the optimum set of parameters is obtained in one step.

Lee, L. T.↗

Viscosity and thermal conductivity of model Jupiter atmospheres

The viscosity and thermal conductivity coefficient are estimated for three models of the atmosphere of Jupiter: a heavy model consisting of 22% helium and 78% hydrogen, a nominal model consisting of 11% helium and 89% hydrogen, and a light model consisting of pure hydrogen. The effect of trace elements is neglected. Linearized approximations are used for the transport coefficients of the mixtures; these are found to be in almost constant ratio to the values for pure hydrogen, independent of temperature. Short Basic language programs for computing the coefficients are listed.

Hansen, C. F.↗

Recent advances in the solution of three-dimensional flow over wings with leading edge vortex separation

Recent advances in a panel method for the solution of three-dimensional flow about wing and wing-body combinations with leading-edge vortex separation are presented. These advances were achieved as part of an ultimately successful assault on two shortcomings of the method, namely convergence failures in seemingly random cases, and overprediction of lift coefficient for high aspect-ratio wings. Advances include the implementation of improved panel numerics for the purpose of eliminating the highly non-linear effects of ring vortices around doublet panel edges, and the development of a least squares procedure for damping vortex sheet geometry update instabilities. A variety of cases generated by the computer program implementing the method are presented. These cases are of two types. The first type consists of numerical studies, which verify the underlying mathematical assumptions of the method and moreover show that the results are strongly invariant with respect to such user dependent input as wing panel layout, initial sheet shape, sheet rollup, etc. The second type consists of cases run for the purpose of comparing computed results with experimental data, and these comparisons verify the underlying physical assumptions made by the method.

Johnson, F. T.↗

Comparison of inlet suppressor data with approximate theory based on cutoff ratio

Inlet suppressor far-field directivity suppression was quantitatively compared with that predicted using an approximate linear design and evaluation method based upon mode cutoff ratio. The experimental data was obtained using a series of cylindrical point-reacting inlet liners on a YF102 engine. The theoretical prediction program is based upon simplified sound propagation concepts derived from exact calculations. These indicate that all of the controlling phenomenon can be approximately correlated with mode cutoff ratio which itself is intimately related to the angles of propagation within the duct. The theory-data comparisons are intended to point out possible deficiencies in the approximate theory which may be corrected. After all theoretical refinements are made, then empirical corrections can be applied.

Rice, E. J.↗

Unsteady Subsonic Loadings Due to Control-Surface Motion

RHOIV computer program predicts unsteady lifting-surface loadings caused by motions of leading edge and trailing-edge control surfaces having sealed gaps at hinge lines and side edges. Analysis is based on linear, subsonic, potential-flow theory using kernel function method. Linear combinations of pressure distributions are used that are continuous except at planform edges and hinge lines. Loading solution is obtained by downwash collocation for deflection or motion of main wing and control surfaces.

Rowe, W. S.↗

The influence of leading-edge load alleviation on supersonic wing design

A theoretical and experimental program to assess the effect of leading-edge load constraints on wing design and performance was conducted. For a planform characterized by a highly swept leading edge on the inboard region, linear theory was used to design camber surfaces which produced minimum drag-due-to-lift at the design lift coefficient of 0.08 and a design Mach number of 2.4. In an effort to delay the formation of leading edge vortices which often occur on highly swept wings, two approaches were used in the design criteria to limit the loadings on the leading edge. One wing was constrained to have the normal Mach number less than one everywhere along the leading edge and the second wing was constrained to have a pressure coefficient of zero on the leading edge. Force tests were run on the two constrained wings, on a flat reference wing and on an optimized wing with no leading edge constraints. All wings had identical planforms and thicknesses and were tested over a range of Mach numbers from 1.8 to 2.8 and a range in angles of attack from -5 deg to 8 deg. A comparison of the experimental performance of these four models is shown. Correlations of these results with theoretical predictions and flow visualization photographs are also included.

Darden, C. M.↗

Cost estimates for flat plate and concentrator collector arrays

The current module and installation costs for the U.S. National Photovoltaic Program's grid-connected systems are significantly higher than required for economic viability of this alternative. Attention is accordingly given to the prospects for installed module cost reductions in flat plate, linear focus Fresnel concentrator, and point focus Fresnel concentrator candidate systems. Cost projections indicate that all three systems would meet near-term and midterm goals, provided that module costs of $2.80/W(p) and $0.70/W(p), respectively, are met. The point focus Fresnel system emerges as the most viable for the near term.

Shimada, K.↗

Modeling of transitional flows

An effort directed at developing improved transitional models was initiated. The focus of this work was concentrated on the critical assessment of a popular existing transitional model developed by McDonald and Fish in 1972. The objective of this effort was to identify the shortcomings of the McDonald-Fish model and to use the insights gained to suggest modifications or alterations of the basic model. In order to evaluate the transitional model, a compressible boundary layer code was required. Accordingly, a two-dimensional compressible boundary layer code was developed. The program was based on a three-point fully implicit finite difference algorithm where the equations were solved in an uncoupled manner with second order extrapolation used to evaluate the non-linear coefficients. Iteration was offered as an option if the extrapolation error could not be tolerated. The differencing scheme was arranged to be second order in both spatial directions on an arbitrarily stretched mesh. A variety of boundary condition options were implemented including specification of an external pressure gradient, specification of a wall temperature distribution, and specification of an external temperature distribution. Overall the results of the initial phase of this work indicate that the McDonald-Fish model does a poor job at predicting the details of the turbulent flow structure during the transition region.

Lund, Thomas S.↗

PROTEUS two-dimensional Navier-Stokes computer code, version 1.0. Volume 2: User's guide

A new computer code was developed to solve the two-dimensional or axisymmetric, Reynolds averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The thin-layer or Euler equations may also be solved. Turbulence is modeled using an algebraic eddy viscosity model. The objective was to develop a code for aerospace applications that is easy to use and easy to modify. Code readability, modularity, and documentation were emphasized. The equations are written in nonorthogonal body-fitted coordinates, and solved by marching in time using a fully-coupled alternating direction-implicit procedure with generalized first- or second-order time differencing. All terms are linearized using second-order Taylor series. The boundary conditions are treated implicitly, and may be steady, unsteady, or spatially periodic. Simple Cartesian or polar grids may be generated internally by the program. More complex geometries require an externally generated computational coordinate system. The documentation is divided into three volumes. Volume 2 is the User's Guide, and describes the program's general features, the input and output, the procedure for setting up initial conditions, the computer resource requirements, the diagnostic messages that may be generated, the job control language used to run the program, and several test cases.

Towne, Charles E.↗

PROTEUS two-dimensional Navier-Stokes computer code, version 1.0. Volume 3: Programmer's reference

A new computer code was developed to solve the 2-D or axisymmetric, Reynolds-averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The thin-layer or Euler equations may also be solved. Turbulence is modeled using an algebraic eddy viscosity model. The objective was to develop a code for aerospace applications that is easy to use and easy to modify. Code readability, modularity, and documentation were emphasized. The equations are written in nonorthogonal body-fitted coordinates, and solved by marching in time using a fully-coupled alternating-direction-implicit procedure with generalized first- or second-order time differencing. All terms are linearized using second-order Taylor series. The boundary conditions are treated implicitly, and may be steady, unsteady, or spatially periodic. Simple Cartesian or polar grids may be generated internally by the program. More complex geometries require an externally generated computational coordinate system. The documentation is divided into three volumes. Volume 3 is the Programmer's Reference, and describes the program structure, the FORTRAN variables stored in common blocks, and the details of each subprogram.

Towne, Charles E.↗

Stratus: An interactive steady state mixed layer model for personal computers

A steady-state, horizontally homogeneous, cloud-topped marine boundary layer model based primarily on the work of Lilly (1968) and Schubert et al., (1979) is presented. The conservative thermodynamic variables are equivalent potential temperature, theta(sub e), and total water mixing ratio, q + l. Some of the differences between this and Lilly's (1968) model are: radiation is allowed to penetrate into the boundary layer; cloud top values of longwave radiation, equivalent potential temperature, and water vapor mixing ratio are linear functions of height derived from climatological data at California coastal stations; and the closure assumption assumes a weighted average of Lilly's (1968) maximum and minimum entrainment theories. This model was programmed in FORTRAN and will run interactively on an IBM-compatible personal computer. The program allows the user to specify the geographical location, the wind speed, the sea-surface temperature, the large scale horizontal divergence, and the initial guess for cloud top height. Output includes the steady state values of cloud top and cloud base height, mixed layer equivalent potential temperature and total water mixing ratio, and the associated convective and radiative fluxes.

Guinn, Thomas A.↗

A Worst-Case Approach for On-Line Flutter Prediction

Worst-case flutter margins may be computed for a linear model with respect to a set of uncertainty operators using the structured singular value. This paper considers an on-line implementation to compute these robust margins in a flight test program. Uncertainty descriptions are updated at test points to account for unmodeled time-varying dynamics of the airplane by ensuring the robust model is not invalidated by measured flight data. Robust margins computed with respect to this uncertainty remain conservative to the changing dynamics throughout the flight. A simulation clearly demonstrates this method can improve the efficiency of flight testing by accurately predicting the flutter margin to improve safety while reducing the necessary flight time.

Lind, Rick C.↗

A Very High Order, Adaptable MESA Implementation for Aeroacoustic Computations

Since computational efficiency and wave resolution scale with accuracy, the ideal would be infinitely high accuracy for problems with widely varying wavelength scales. Currently, many of the computational aeroacoustics methods are limited to 4th order accurate Runge-Kutta methods in time which limits their resolution and efficiency. However, a new procedure for implementing the Modified Expansion Solution Approximation (MESA) schemes, based upon Hermitian divided differences, is presented which extends the effective accuracy of the MESA schemes to 57th order in space and time when using 128 bit floating point precision. This new approach has the advantages of reducing round-off error, being easy to program. and is more computationally efficient when compared to previous approaches. Its accuracy is limited only by the floating point hardware. The advantages of this new approach are demonstrated by solving the linearized Euler equations in an open bi-periodic domain. A 500th order MESA scheme can now be created in seconds, making these schemes ideally suited for the next generation of high performance 256-bit (double quadruple) or higher precision computers. This ease of creation makes it possible to adapt the algorithm to the mesh in time instead of its converse: this is ideal for resolving varying wavelength scales which occur in noise generation simulations. And finally, the sources of round-off error which effect the very high order methods are examined and remedies provided that effectively increase the accuracy of the MESA schemes while using current computer technology.

Dydson, Roger W.↗

Deriving the Coronal Magnetic Field Using Parametric Transformation Analysis

When plasma-beta greater than 1 then the gas pressure dominates over the magnetic pressure. This ratio as a function along the coronal magnetic field lines varies from beta greater than 1 in the photosphere at the base of the field lines, to beta much less than 1 in the mid-corona, to beta greater than 1 in the upper corona. Almost all magnetic field extrapolations do not or cannot take into account the full range of beta. They essentially assume beta much less than 1, since the full boundary conditions do not exist in the beta greater than 1 regions. We use a basic parametric representation of the magnetic field lines such that the field lines can be manipulated to match linear features in the EUV and SXR coronal images in a least squares sense. This research employs free-form deformation mathematics to generate the associated coronal magnetic field. In our research program, the complex magnetic field topology uses Parametric Transformation Analysis (PTA) which is a new and innovative method to describe the coronal fields that we are developing. In this technique the field lines can be viewed as being embedded in a plastic medium, the frozen-in-field-line concept. As the medium is deformed the field lines are similarly deformed. However the advantage of the PTA method is that the field line movement represents a transformation of one magnetic field solution into another magnetic field solution. When fully implemented, this method will allow the resulting magnetic field solution to fully match the magnetic field lines with EUV/SXR coronal loops by minimizing the differences in direction and dispersion of a collection of PTA magnetic field lines and observed field lines. The derived magnetic field will then allow beta greater than 1 regions to be included, the electric currents to be calculated, and the Lorentz force to be determined. The advantage of this technique is that the solution is: (1) independent of the upper and side boundary conditions, (2) allows non-vanishing magnetic forces, and (3) provides a global magnetic field solution, which contains high- and low-beta regimes and maximizes the similarity between the field lines structure and all the coronal images of the region. The coronal image analysis is crucial to the investigation and for the first time these images can be exploited to derive the coronal magnetic field in a well-posed mathematical formulation. This program is an outgrowth of an investigation in which an extrapolated potential field was required to be "inflated" in order to have the field lines match the Yohkoh/SXT images. The field lines were radially stretched resulting in a better match to the coronal loops of an active region. The PTA method of radial and non-radial deformations of field lines to provide a match to the EUV/SXR images will be presented.

Gary, G. Allen↗