Search NASA⌕ Search

SEARCH · Search NASA

Results for “transfer matrix calculations”

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 37 records · Page 2

Theoretical and material studies on thin-film electroluminescent devices

A theoretical study of resonant tunneling in multilayered heterostructures is presented based on an exact solution of the Schroedinger equation under the application of a constant electric field. By use of the transfer matrix approach, the transmissivity of the structure is determined as a function of the incident electron energy. The approach presented is easily extended to many layer structures where it is more accurate than other existing transfer matrix or WKB models. The transmission resonances are compared to the bound state energies calculated for a finite square well under bias using either an asymmetric square well model or the exact solution of an infinite square well under the application of an electric field. The results show good agreement with other existing models as well as with the bound state energies. The calculations were then applied to a new superlattice structure, the variablly spaced superlattice energy filter, (VSSEP) which is designed such that under bias the spatial quantization levels fully align. Based on these calculations, a new class of resonant tunneling superlattice devices can be designed.

Summers, C. J.↗

Development and testing of methodology for evaluating the performance of multi-input/multi-output digital control systems

A Controller Performance Evaluation (CPE) methodology for multi-input/multi-output digital control systems was developed and tested on an aeroelastic wind-tunnel model. Modern signal processing methods were used to implement control laws and to acquire time domain data of the whole system (controller and plant) from which appropriate transfer matrices of the control system could be generated. Matrix computational procedures were used to calculate singular values of return-difference matrices at the plant input and output points to evaluate the performance of the control system. The CPE procedures effectively identified potentially destabilizing controllers and confirmed the satisfactory performance of stabilizing ones.

Pototzky, Anthony S.↗

Thermohydrodynamic analysis of cryogenic liquid turbulent flow fluid film bearings, phase 2

The Phase 2 (1994) Annual Progress Report presents two major report sections describing the thermal analysis of tilting- and flexure-pad hybrid bearings, and the unsteady flow and transient response of a point mass rotor supported on fluid film bearings. A literature review on the subject of two-phase flow in fluid film bearings and part of the proposed work for 1995 are also included. The programs delivered at the end of 1994 are named hydroflext and hydrotran. Both codes are fully compatible with the hydrosealt (1993) program. The new programs retain the same calculating options of hydrosealt plus the added bearing geometries, and unsteady flow and transient forced response. Refer to the hydroflext & hydrotran User's Manual and Tutorial for basic information on the analysis and instructions to run the programs. The Examples Handbook contains the test bearing cases along with comparisons with experimental data or published analytical values. The following major tasks were completed in 1994 (Phase 2): (1) extension of the thermohydrodynamic analysis and development of computer program hydroflext to model various bearing geometries, namely, tilting-pad hydrodynamic journal bearings, flexure-pad cylindrical bearings (hydrostatic and hydrodynamic), and cylindrical pad bearings with a simple elastic matrix (ideal foil bearings); (2) improved thermal model including radial heat transfer through the bearing stator; (3) calculation of the unsteady bulk-flow field in fluid film bearings and the transient response of a point mass rotor supported on bearings; and (4) a literature review on the subject of two-phase flows and homogeneous-mixture flows in thin-film geometries.

Sanandres, Luis↗

Critical speeds and forced response solutions for active magnetic bearing turbomachinery, part 1

The prediction of critical speeds and forced response of active magnetic bearing turbomachinery is of great interest due to the increased use of this new and promising technology. Calculating the system undamped critical speeds and forced response is important to all those who are involved in the design of the active magnetic bearing system. An extended Jeffcott model which was used as an approximate solution to a more accurate transfer matrix procedure is presented. Theory behind a two-degree-of freedom extended Jeffcoat model is presented. Results of the natural frequency calculation are shown followed by the results of the forced response calculation. The system response was predicted for two types of forcing. A constant magnitude excitation with a wide frequency variation was applied at the bearings as one forcing function. The normal unbalance force at the midspan was the second source of excitation. The results of this extended Jeffcott solution gives useful design guidance for the influence of the first and third modes of a symmetric rotor system.

Keesee, J.↗

Higher harmonic control analysis for vibration reduction of helicopter rotor systems

An advanced higher harmonic control (HHC) analysis has been developed and applied to investigate its effect on vibration reduction levels, blade and control system fatigue loads, rotor performance, and power requirements of servo-actuators. The analysis is based on a finite element method in space and time. A nonlinear time domain unsteady aerodynamic model, based on the indicial response formulation, is used to calculate the airloads. The rotor induced inflow is computed using a free wake model. The vehicle trim controls and blade steady responses are solved as one coupled solution using a modified Newton method. A linear frequency-domain quasi-steady transfer matrix is used to relate the harmonics of the vibratory hub loads to the harmonics of the HHC inputs. Optimal HHC is calculated from the minimization of the vibratory hub loads expressed in term of a quadratic performance index. Predicted vibratory hub shears are correlated with wind tunnel data. The fixed-gain HHC controller suppresses completely the vibratory hub shears for most of steady or quasi-steady flight conditions. HHC actuator amplitudes and power increase significantly at high forward speeds (above 100 knots). Due to the applied HHC, the blade torsional stresses and control loads are increased substantially. For flight conditions where the blades are stalled considerably, the HHC input-output model is quite nonlinear. For such cases, the adaptive-gain controller is effective in suppressing vibratory hub loads, even though HHC may actually increase stall areas on the rotor disk. The fixed-gain controller performs poorly for such flight conditions. Comparison study of different rotor systems indicates that a soft-inplane hingeless rotor requires less actuator power at high speeds (above 130 knots) than an articulated rotor, and a stiff-inplane hingeless rotor generally requires more actuator power than an articulated or a soft-inplane hingeless rotor. Parametric studies for a hingeless rotor operating in a transition flight regime and for an articulated rotor operating at the level-flight boundary (high speed and high thrust conditions) indicate that blade parameters including flap, lag, torsion stiffness distributions, linear pretwist, chordwise offset of center-of-mass from elastic axis and chordwise offset of elastic axis from aerodynamic center can be selected to minimize the actuator power requirements for HHC.

Nguyen, Khanh Q.↗

A simplified analytical solution for thermal response of a one-dimensional, steady state transpiration cooling system in radiative and convective environment

A simplified analytical method for calculation of thermal response within a transpiration-cooled porous heat shield material in an intense radiative-convective heating environment is presented. The essential assumptions of the radiative and convective transfer processes in the heat shield matrix are the two-temperature approximation and the specified radiative-convective heatings of the front surface. Sample calculations for porous silica with CO2 injection are presented for some typical parameters of mass injection rate, porosity, and material thickness. The effect of these parameters on the cooling system is discussed.

Kubota, H.↗

Polarized Bidirectional Reflectance of Optically Thick Sparse Particulate Layers: an Efficient Numerically Exact Radiative-Transfer Solution

We describe a simple yet efficient numerical algorithm for computing polarized bidirectional reflectance of an optically thick (semi-infinite), macroscopically flat layer composed of statistically isotropic and mirror symmetric random particles. The spatial distribution of the particles is assumed to be sparse, random, and statistically uniform. The 44 Stokes reflection matrix is calculated by iterating the Ambartsumian's vector nonlinear integral equation. The result is a numerically exact solution of the vector radiative transfer equation and as such fully satisfies the energy conservation law and the fundamental reciprocity relation. Since this technique bypasses the computation of the internal radiation field, it is very fast and highly accurate. The FORTRAN implementation of the technique is publicly available on the World Wide Web at http://www.giss.nasa.gov/staff/ mmishchenko/brf. It can be combined with several existing computer programs providing the requisite single-scattering properties of spherical or morphologically complex particles and applied to a wide range of optical characterization problems. Benchmark results obtained with this program can be used for testing alternative solvers of the vector radiative transfer equation.

radiative transfer↗

Working Toward Nanotube Composites

One of the most attractive applications of single-wall carbon nanotubes (SWNT) is found in the area of structural materials. Nanotubes have a unique combination of high strength, modulus, and elongation to failure, and therefore have potential to significantly enhance the mechanical properties of today's composites. This is especially attractive for the aerospace industry looking for any chance to save weight. This is why NASA has chosen to tackle this difficult application of SWNT. Nanotube properties differ significantly from that of conventional carbon fibers, and a whole new set of problems, including adhesion and dispersion in the adhesive polymer matrix, must be resolved in order to engineer superior composite materials. From recent work on a variety of applications it is obvious that the wide range of research in nanotubes will lead to advances in physics, chemistry, and engineering. However, the possibility of ultralightweight structures is what causes dreamers to really get excited. One of the important issues in composite engineering is aspect ratio of the fibers, since it affects load transfer in composites. Nanotube length was a gray area for years, since they are formed in bundles, making it impossible to monitor individual nanotube length. Even though bundles are observed to be tens and hundreds of microns long, they can be built of relatively short tubes weakly bound by Van der Waals forces. Nanotube length can be affected by subsequent purification and ultrasound processing, which has been necessary in order to disperse nanotubes and introduce them into a polymer matrix. Some calculations show that nanotubes with 10(exp 5) aspect ratio may be necessary to achieve good load transfer. We show here that nanotubes produced in our laser system are as much as tens of microns long and get cut into lengths of hundreds of nanometers during ultrasound processing. Nanotube length was measured by AFM on pristine nanotube specimens as well, as after sonication. In each case great care was taken to measure individual nanotubes, rather than bundles. Pristine nanotubes were collected on quartz substrates placed directly in the laser oven and exposed for 0.5 s. This results in an equal mix of bundles and individual nanotubes. Nanotube length measurements were limited practically by the lateral span of the AFM scanner, but nanotube length is certainly in excess of 20 micrometers.

Arepalli, Sivaram↗

Vibrations and structureborne noise in space station

Analytical models were developed to predict vibrations and structureborne noise generation of cylindrical and rectangular acoustic enclosures. These models are then used to determine structural vibration levels and interior noise to random point input forces. The guidelines developed could provide preliminary information on acoustical and vibrational environments in space station habitability modules under orbital operations. The structural models include single wall monocoque shell, double wall shell, stiffened orthotropic shell, descretely stiffened flat panels, and a coupled system composed of a cantilever beam structure and a stiffened sidewall. Aluminum and fiber reinforced composite materials are considered for single and double wall shells. The end caps of the cylindrical enclosures are modeled either as single or double wall circular plates. Sound generation in the interior space is calculated by coupling the structural vibrations to the acoustic field in the enclosure. Modal methods and transfer matrix techniques are used to obtain structural vibrations. Parametric studies are performed to determine the sensitivity of interior noise environment to changes in input, geometric and structural conditions.

Vaicaitis, R.↗

A case study of view-factor rectification procedures for diffuse-gray radiation enclosure computations

The view factors which are used in diffuse-gray radiation enclosure calculations are often computed by approximate numerical integrations. These approximately calculated view factors will usually not satisfy the important physical constraints of reciprocity and closure. In this paper several view-factor rectification algorithms are reviewed and a rectification algorithm based on a least-squares numerical filtering scheme is proposed with both weighted and unweighted classes. A Monte-Carlo investigation is undertaken to study the propagation of view-factor and surface-area uncertainties into the heat transfer results of the diffuse-gray enclosure calculations. It is found that the weighted least-squares algorithm is vastly superior to the other rectification schemes for the reduction of the heat-flux sensitivities to view-factor uncertainties. In a sample problem, which has proven to be very sensitive to uncertainties in view factor, the heat transfer calculations with weighted least-squares rectified view factors are very good with an original view-factor matrix computed to only one-digit accuracy. All of the algorithms had roughly equivalent effects on the reduction in sensitivity to area uncertainty in this case study.

Taylor, Robert P.↗

Theoretical Investigation of OCN(-) Charge Transfer Complexes in Condensed Phase Media: Spectroscopic Properties in Amorphous Ice

Density functional theory (DFT) calculations of cyanate (OCN(-)) charge-transfer complexes were performed to model the "XCN" feature observed in interstellar icy grain mantles. OCN(-) charge-transfer complexes were formed from precursor combinations of HNCO or HOCN with either NH3 or H2O. Three different solvation strategies for realistically modeling the ice matrix environment were explored, including (1) continuum solvation, (2) pure DFT cluster calculations, and (3) an ONIOM DFT/PM3 cluster calculation. The model complexes were evaluated by their ability to reproduce seven spectroscopic measurements associated with XCN: the band origin of the OCN(-) asymmetric stretching mode, shifts in that frequency due to isotopic substitutions of C, N, O, and H, plus two weak features. The continuum solvent field method produced results consistent with some of the experimental data but failed to account for other behavior due to its limited capacity to describe molecular interactions with solvent. DFT cluster calculations successfully reproduced the available spectroscopic measurements very well. In particular, the deuterium shift showed excellent agreement in complexes where OCN(-) was fully solvated. Detailed studies of representative complexes including from two to twelve water molecules allowed the exploration of various possible solvation structures and provided insights into solvation trends. Moreover, complexes arising from cyanic or isocyanic acid in pure water suggested an alternative mechanism for the formation of OCN(-) charge-transfer complexes without the need for a strong base such as NH3 to be present. An extended ONIOM (B3LYP/PM3) cluster calculation was also performed to assess the impact of a more realistic environment on HNCO dissociation in pure water.

Park, Jin-Young↗

Multicyclic control of a helicopter rotor considering the influence of vibration, loads, and control motion

Weighted multiple linear regression is used to establish a transfer function matrix relationship between higher harmonic control inputs and transducer vibration outputs for a controllable twist rotor. Data used in the regression were taken from the test of a KAMAN controllable twist rotor conducted in the Ames Research Center's 40- by 80-Foot Wind Tunnel in June 1977. Optimal controls to minimize fixed system vibrational levels are calculated using linear quadratic regulatory theory with a control deflection penalty included in the performance criteria. Control sensitivity to changes in control travel, forward speed, and lift and propulsive forces is examined. It is found that the linear transfer matrix is a strong function of forward speed and a weak function of lift and propulsive force. An open-loop strategy is proposed for systems with limited control travel.

Brown, T. J.↗

Theory of resonant tunneling in a variably spaced multiquantum well structure - An Airy function approach

A theoretical study of resonant tunneling in multilayered heterostructures is presented based on an exact solution of the Schroedinger equation under the application of a constant electric field. By use of the transfer matrix approach, the transmissivity of the structure is determined as a function of the incident electron energy. The approach presented herein is easily extended to many layer structures where it is more accurate than other existing transfer matrix or Wentzel-Kramers-Brillouin (WKB) models. The transmission resonances are compared to the bound-state energies calculated for a finite square well under bias using either an asymmetric square-well model or the exact solution of an infinite square well under the application of an electric field. The results show good agreement with other existing models as well as with the bound-state energies. The calculations were then applied to a new superlattice structure, the variably spaced superlattice energy filter, which is designed such that under bias the spatial quantization levels fully align. Based on these calculations, a new class of resonant tunneling superlattice devices can be designed.

Brennan, K. F.↗

Momentum loss in proton-nucleus and nucleus-nucleus collisions

An optical model description, based on multiple scattering theory, of longitudinal momentum loss in proton-nucleus and nucleus-nucleus collisions is presented. The crucial role of the imaginary component of the nucleon-nucleon transition matrix in accounting for longitudinal momentum transfer is demonstrated. Results obtained with this model are compared with Intranuclear Cascade (INC) calculations, as well as with predictions from Vlasov-Uehling-Uhlenbeck (VUU) and quantum molecular dynamics (QMD) simulations. Comparisons are also made with experimental data where available. These indicate that the present model is adequate to account for longitudinal momentum transfer in both proton-nucleus and nucleus-nucleus collisions over a wide range of energies.

Khan, Ferdous↗

Electromagnetic scattering calculations on the Intel Touchstone Delta

During the first year's operation of the Intel Touchstone Delta system, software which solves the electric field integral equations for fields scattered from arbitrarily shaped objects has been transferred to the Delta. To fully realize the Delta's resources, an out-of-core dense matrix solution algorithm that utilizes some or all of the 90 Gbyte of concurrent file system (CFS) has been used. The largest calculation completed to date computes the fields scattered from a perfectly conducting sphere modeled by 48,672 unknown functions, resulting in a complex valued dense matrix needing 37.9 Gbyte of storage. The out-of-core LU matrix factorization algorithm was executed in 8.25 h at a rate of 10.35 Gflops. Total time to complete the calculation was 19.7 h-the additional time was used to compute the 48,672 x 48,672 matrix entries, solve the system for a given excitation, and compute observable quantities. The calculation was performed in 64-b precision.

Cwik, Tom↗

Structural analyses of engine wall cooling concepts and materials

The severe thermal environments under which hypersonic aircraft such as the National Aerospace Plane (NASP) will operate require cooling of the engine walls, especially in the combustor. A preliminary assessment is made of some candidate materials based on structural analyses for a number of convective cooling configurations. Three materials were studied: graphite/copper and tungsten/copper composite alloys with 50 percent fiber volume fractions and a wrought cobalt-base superalloy, Haynes 188. Anisotropic mechanical and thermal properties for the composites were obtained from a computer code, ICAN, which determines the composite material properties from the individual properties of the fiber and matrix materials. The structural analyses were performed by using the MARC nonlinear finite element code. Heat transfer analyses were conducted to calculate the metal temperature distributions.

Kaufman, Albert↗

Tunneling calculations for GaAs-Al(x)Ga(1-x) as graded band-gap sawtooth superlattices

Quantum mechanical tunneling calculations for sawtooth (linearly graded band-gap) and step-barrier AlGaAs superlattices were performed by means of a transfer matrix method, within the effective mass approximation. The transmission coefficient and tunneling current versus applied voltage were computed for several representative structures. Particular consideration was given to effective mass variations. The tunneling properties of step and sawtooth superlattices show some qualitative similarities. Both structures exhibit resonant tunneling, however, because they deform differently under applied fields, the J-V curves differ.

Forrest, Kathrine A.↗

Experimental investigation of two-phase heat transfer in a porous matrix.

One-dimensional two-phase flow transpiration cooling through porous metal is studied experimentally. The experimental data is compared with a previous one-dimensional analysis. Good agreement with calculated temperature distribution is obtained as long as the basic assumptions of the analytical model are satisfied. Deviations from the basic assumptions are caused by nonhomogeneous and oscillating flow conditions. Preliminary derivation of nondimensional parameters which characterize the stable and unstable flow conditions is given. Superheated liquid droplets observed sputtering from the heated surface indicated incomplete evaporation at heat fluxes well in access of the latent energy transport. A parameter is developed to account for the nonequilibrium thermodynamic effects. Measured and calculated pressure drops show contradicting trends which are attributed to capillary forces.

Von Reth, R.↗