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 55 records · Page 3

A multi-layer discrete-ordinate method for vector radiative transfer in a vertically-inhomogeneous, emitting and scattering atmosphere. I - Theory. II - Application

A theory is developed for discretizing the vector integro-differential radiative transfer equation including both solar and thermal radiation. A complete solution and boundary equations are obtained using the discrete-ordinate method. An efficient numerical procedure is presented for calculating the phase matrix and achieving computational stability. With natural light used as a beam source, the Stokes parameters from the model proposed here are compared with the analytical solutions of Chandrasekhar (1960) for a Rayleigh scattering atmosphere. The model is then applied to microwave frequencies with a thermal source, and the brightness temperatures are compared with those from Stamnes'(1988) radiative transfer model.

Weng, Fuzhong↗

Radiative transfer theory for active remote sensing of a layer of nonspherical particles

The radiative transfer theory is applied to calculate the scattering by a layer of randomly positioned and oriented nonspherical particles. The scattering amplitude functions of each individual particle are calculated with Waterman's T matrix method, which utilizes vector spherical wave functions for expansion of incident, scattered, and surface fields. The orientation of the particles is described by a probability density function of the Eulerian angles of rotation. A rotation matrix is used to relate the T matrix of the principal frame to that of the natural frame of the particle. The extinction matrix and phase matrix of the radiative transfer equations are expressed in terms of the T matrix elements. The extinction matrix for nonspherical particles is generally nondiagonal. There are only two attenuation rates in a specified direction of propagation. The radiative transfer equations are solved by an iterative method to first order in albedo. Numerical results are illustrated as functions of incidence angle and frequency with applications to active remote sensing.

Tsang, L.↗

Water vapor profile retrieval simulation studies for the HIRS/MSU and AMTS/MSU sounders

A combined physical statistical method for the retrieval of atmospheric water vapor column densities has been developed at the Goddard Laboratory for Atmospheric Sciences. In this method differences between observed and calculated clear column brightness temperatures are used in a regression formalism to obtain estimates of integrated column densities in given atmospheric layers. The physics of the radiative transfer problem is contained in the calculated clear column brightness temperatures and in accounting for the effects of clouds on the observed radiances. The regression matrix is equivalent to a set of empirical sensitivity factors.

Reuter, D.↗

Computational architecture for integrated controls and structures design

To facilitate the development of control structure interaction (CSI) design methodology, a computational architecture for interdisciplinary design of active structures is presented. The emphasis of the computational procedure is to exploit existing sparse matrix structural analysis techniques, in-core data transfer with control synthesis programs, and versatility in the optimization methodology to avoid unnecessary structural or control calculations. The architecture is designed such that all required structure, control and optimization analyses are performed within one program. Hence, the optimization strategy is not unduly constrained by cold starts of existing structural analysis and control synthesis packages.

Belvin, W. Keith↗

Control-structure-thermal interactions in analysis of lunar telescopes

The lunar telescope project was an excellent model for the CSTI study because a telescope is a very sensitive instrument, and thermal expansion or mechanical vibration of the mirror assemblies will rapidly degrade the resolution of the device. Consequently, the interactions are strongly coupled. The lunar surface experiences very large temperature variations that range from approximately -180 C to over 100 C. Although the optical assemblies of the telescopes will be well insulated, the temperature of the mirrors will inevitably fluctuate in a similar cycle, but of much smaller magnitude. In order to obtain images of high quality and clarity, allowable thermal deformations of any point on a mirror must be less than 1 micron. Initial estimates indicate that this corresponds to a temperature variation of much less than 1 deg through the thickness of the mirror. Therefore, a lunar telescope design will most probably include active thermal control, a means of controlling the shape of the mirrors, or a combination of both systems. Historically, the design of a complex vehicle was primarily a sequential process in which the basic structure was defined without concurrent detailed analyses or other subsystems. The basic configuration was then passed to the different teams responsible for each subsystem, and their task was to produce a workable solution without requiring major alterations to any principal components or subsystems. Consequently, the final design of the vehicle was not always the most efficient, owing to the fact that each subsystem design was partially constrained by the previous work. This procedure was necessary at the time because the analysis process was extremely time-consuming and had to be started over with each significant alteration of the vehicle. With recent advances in the power and capacity of small computers, and the parallel development of powerful software in structural, thermal, and control system analysis, it is now possible to produce very detailed analyses of intermediate designs in a much shorter period of time. The subsystems can thus be designed concurrently, and alterations in the overall design can be quickly adopted into each analysis; the design becomes an iterative process in which it is much easier to experiment with new ideas, configurations, and components. Concurrent engineering has the potential to produce efficient, highly capable designs because the effect of one subystem on another can be assessed in much more detail at a very early point in the program. The research program consisted of several tasks: scale a prototype telescope assembly to a 1 m aperture, develop a model of the telescope assembly by using finite element (FEM) codes that are available on site, determine structural deflections of the mirror surfaces due to the temperature variations, develop a prototype control system to maintain the proper shape of the optical elements, and most important of all, demonstrate the concurrent engineering approach with this example. In addition, the software used for the finite element models and thermal analysis was relatively new within the Program Development Office and had yet to be applied to systems this large or complex; understanding the software and modifying it for use with this project was also required. The I-DEAS software by Structural Dynamics Research Corporation (SDRC) was used to build the finite element models, and TMG developed by Maya Heat Transfer Technologies, Ltd. (which runs as an I-DEAS module) was used for the thermal model calculations. All control system development was accomplished with MATRIX(sub X) by Integrated Systems, Inc.

Thompson, Roger C.↗

Matrix operator theory of radiative transfer. I - Rayleigh scattering.

An entirely rigorous method for the solution of the equations for radiative transfer based on the matrix operator theory is reviewed. The advantages of the present method are: (1) all orders of the reflection and transmission matrices are calculated at once; (2) layers of any thickness may be combined, so that a realistic model of the atmosphere can be developed from any arbitrary number of layers, each with different properties and thicknesses; (3) calculations can readily be made for large optical depths and with highly anisotropic phase functions; (4) results are obtained for any desired value of the surface albedo including the value unity and for a large number of polar and azimuthal angles; (5) all fundamental equations can be interpreted immediately in terms of the physical interactions appropriate to the problem; and (6) both upward and downward radiance can be calculated at interior points from relatively simple expressions.

Plass, G. N.↗

Multiple scattering effects on wave propagation due to rain

The paper examines the multiple scattering of microwaves due to rain and the resulting incoherent intensity for a plane wave incident on a plane-parallel rain region. General formulations of the transfer equation using Stokes parameters are presented, and an extinction matrix which takes into account depolarization effects and nonspherical droplet shape is introduced. Scattering characteristics are calculated for spherical droplets using the Mie solution and the Laws-Parsons distributions. The ratios of incoherent to coherent received powers are calculated at 30 GHz for rain rates of 12.5, 50, 100, and 150 mm/h, rain thicknesses of 3 km and 10 km, and fields of view of 1.5, 5, and 15 deg.

Ishimaru, A.↗

A test of the applicability of independent scattering to high albedo planetary regoliths

We show that 'independent scattering' is a useful approximation for high albedo particles whose size and packing density are similar to typical particles in the lunar regolith. Laboratory measurements of the intensity and linear polarization of light scattered from a laboratory sample of glass spheres of know size and composition are compared to radiative transfer calculations of the same observable quantities. Mie scattering is integrated over the size distribution of the particles to determine the mean phase and polarization functions, or phase matrix, of the particle. Assuming that the particles scatter independently, the 'doubling method' is used to rigorously calculate multiple scattering in an optically thick layer of these anisotropically scattering particles. All of the major features 'predicted' in the calculated intensity (double peaks at small phase angles) and polarization (negative branch at small phase angles, large positive peak near 20 degrees phase, and small polarization elsewhere) are observed in the laboratory measurements, with good quantitative agreement indicated at phase angles less than 90 degrees. Even though the particles are supported by physical contact with each other in the sample, as are the particles in planetary regoliths, the independent scattering calculation yields a good approximation to both the intensity and polarization. The physical parameters input to the calculation are only the size distribution of the particles and their complex index of refraction (composition). Significant advantages of this approach are that the phase matrix is calculated from basic physical principles and both the intensity and polarization are determined simultaneously. This model may have broad applications to the interpretation of photometry, spectroscopy, and polarimetry of the icy regoliths of high albedo satellites. The intent of this effort is to perform a controlled experiment that tests the utility of the independent scattering assumption

Goguen, J. D.↗

A general and computationally fast formulation for radiative transfer with scattering

A general formulation of monocromatic radiative transfer with scattering has been developed for plane-parallel geometry. The inhomogeneous and nonisothermal medium absorbs, emits, and anisotropically scatters radiation. Surfaces can emit and scatter radiation in any specified manner. The solution procedure uses the fact that phase incoherent scattering is linear in radiative sources. Certain basic scattering functions are then defined and calculated by an adding computer code using matrix algebra. These scattering functions are weighted by the temperature field and summed (superimposed) to obtain the solution for any specific problem. Numerical results for exiting intensities and one-sided heat fluxes from general media bound by one arbitrary surface are presented. These parametric studies demonstrate the effects of scattering particles and surfaces on radiative transfer from inhomogeneous and nonisothermal media. Application of the formulation to radiative equilibrium is also discussed. The conclusion is that all problems in plane-parallel radiative transfer with scattering can be solved by a common and computationally fast algorithm based on this formulation.

Cogley, A. C.↗

Electron scattering by highly polar molecules. III - CsCl

Utilizing a crossed electron-beam-molecular-beam scattering geometry, relative values of differential electron scattering cross sections for cesium chloride at 5 and 20 eV electron impact energies and at scattering angles between 10 and 120 deg have been measured. These relative cross sections have been normalized to the cross section at 15 deg scattering angle calculated by the hybrid S-matrix technique. In the angular range between 0 and 10 deg and between 120 and 180 deg extrapolations have been made to obtain integral and momentum transfer cross sections. An energy-loss spectrum is also presented which gives various spectral features lying between the 4 and 10 eV regions in CsCl.

Vuskovic, L.↗

Effects of multiple scattering and surface albedo on the photochemistry of the troposphere

The effect of treatment of incoming solar radiation on the photochemistry of the troposphere is discussed. A one dimensional photochemical model of the troposphere containing the species of the nitrogen, oxygen, carbon, hydrogen, and sulfur families was developed. The vertical flux is simulated by use of the parameterized eddy diffusion coefficients. The photochemical model is coupled to a radiative transfer model that calculates the radiation field due to the incoming solar radiation which initiates much of the photochemistry of the troposphere. Vertical profiles of tropospheric species were compared with the Leighton approximation, radiative transfer, matrix inversion model. The radiative transfer code includes the effects of multiple scattering due to molecules and aerosols, pure absorption, and surface albedo on the transfer of incoming solar radiation. It is indicated that significant differences exist for several key photolysis frequencies and species number density profiles between the Leighton approximation and the profiles generated with, radiative transfer, matrix inversion technique. Most species show enhanced vertical profiles when the more realistic treatment of the incoming solar radiation field is included

Augustsson, T. R.↗

Boundary layer simulator (BLIMPJ) improvement

Improvements were made to advance the current Boundary Layer Integral Matrix Procedure - Version J (BLIMPJ) containing previously modeled simplified calculation methods accounting for surface roughness, relaminarization, and thick boundary layer effects. These improvements have potential applications in the design of the future Orbit Transfer Vehicles (OTV) engines. The coded wall roughness module was validated successfully against measured data from the LeRC 60-15 degree half-angle conical nozzle. The relaminarization model was modified to include wall cooling and roughness based on a set of test data from the same nozzle. The model was compared against typical data sets for validation in different flow regimes. The thrust loss calculation for thick boundary layer effects was computerized where the iterations between TDK and BLIMPJ proceeded automatically until convergence was achieved. An example is given for a typical OTV nozzle with an area ratio 1293.

Praharaj, S. C.↗

Response of discretely stiffened structures

This paper presents a theoretical study on dynamic response of complex discretely stiffened and interconnected structures to random loads. Modal and transfer matrix techniques are utilized to analyze these multidimensional structures. Application of this formulation for response and noise generation studies in space station and fuselage-wing design is discussed. The response spectral densities are calculated for several typical discretely stiffened and interconnected structural systems. It is found that this approach allows for easy parametric evaluation and that by proper selection of structural parameters and loading conditions, the vibration response can be reduced.

Vaicaitis, R.↗

Asymptotic unbounded root loci - Formulas and computation

A new geometric way of computing the asymptotic behavior of unbounded root loci of a strictly proper linear time-invariant control system as loop gain goes to infinity is presented. Properties of certain restricted linear maps and nested restrictions of linear maps are developed, and formulas are obtained for the leading coefficient of the asymptotic values of the unbounded multivariable root loci are obtained in terms of eigenvalues of those maps. Published results and a certain simple null structure assumption are used to relate these asymptotic values to the structure at infinity of the Smith-McMillan form of the open loop transfer function. Explicit matrix formulas for the more abstract derived formulas are given and additional geometric insights are developed with orthogonal projections and singular value decomposition. Formulas for the pivots of the unbounded root loci are calculated and shown to have the same form as the coefficients of the unbounded asymptotic root loci.

Sastry, S. S.↗

Thermal histories of the samples of two KOSI comet nucleus simulation experiments

Temperatures recorded during two KOSI comet nucleus simulation experiments strongly suggest that heat transport by vapor flow into the interior of the sample is very important. Two comet nucleus simulation experiments have been done by the KOSI team in a big space simulator. The thermal evolution of the sample during insolation and the results of simplified thermal evolution calculations are discussed. The observed thermal histories cannot be explained by a simple model with heat transferred by heat conduction at a constant conductivity, so a coupled heat and mass transfer problem was considered. The porous ice matrix was assumed to have a constant thermal conductivity and to be in thermal equilibrium with vapor in the pores, the internal pressure being the vapor pressure. The vapor was modelled as an ideal gas because, at the temperatures relevant to the problem, the mean free path length of the vapor molecules is large in comparison with the pore dimensions. The heat capacity at constant volume per unit mass of the two phase mixture was also assumed constant. The vapor was allowed to flow and transfer heat in response to an internal pressure gradient.

Spohn, T.↗

Modal analysis of turborotors using planar modes - Theory

The generalized dynamic equations of motion have been obtained by the direct stiffness method for multimass flexible rotor bearing systems including the effects of gyroscopic moments, disk skew, and rotor acceleration. A set of undamped critical speed mode shapes calculated from the average horizontal and vertical bearing stiffness is used to transform the equations of motion into a set of coupled modal equations of motion. The modal equations are coupled by the generalized bearing coefficients and the gyroscopic moments. An analysis using only undamped critical speeds or decoupled modal analysis assuming proportional damping may lead to erroneous results. This paper presents a rapid method of calculating rotor resonance speeds with their corresponding amplification factors, stability and unbalance response of turborotors. Examples of the application of this modal approach are presented and results are compared to those of other methods such as matrix transfer analysis.

Gunter, E. J.↗

Optical Control of Intersubband Absorption in a Multiple Quantum Well-Embedded Semiconductor Microcravity

Optical intersubband response of a multiple quantum well (MQW)-embedded microcavity driven by a coherent pump field is studied theoretically. The n-type doped MQW structure with three subbands in the conduction band is sandwiched between a semi-infinite medium and a distributed Bragg reflector (DBR). A strong pump field couples the two upper subbands and a weak field probes the two lower subbands. To describe the optical response of the MQW-embedded microcavity, we adopt a semi-classical nonlocal response theory. Taking into account the pump-probe interaction, we derive the probe-induced current density associated with intersubband transitions from the single-particle density-matrix formalism. By incorporating the current density into the Maxwell equation, we solve the probe local field exactly by means of Green's function technique and the transfer-matrix method. We obtain an exact expression for the probe absorption coefficient of the microcavity. For a GaAs/Al(sub x)Ga(sub 1-x)As MQW structure sandwiched between a GaAs/AlAs DBR and vacuum, we performed numerical calculations of the probe absorption spectra for different parameters such as pump intensity, pump detuning, and cavity length. We find that the probe spectrum is strongly dependent on these parameters. In particular, we find that the combination of the cavity effect and the Autler-Townes effect results in a triplet in the optical spectrum of the MQW system. The optical absorption peak value and its location can be feasibly controlled by varying the pump intensity and detuning.

Liu, Ansheng↗

Algorithms for Efficient Computation of Transfer Functions for Large Order Flexible Systems

An efficient and robust computational scheme is given for the calculation of the frequency response function of a large order, flexible system implemented with a linear, time invariant control system. Advantage is taken of the highly structured sparsity of the system matrix of the plant based on a model of the structure using normal mode coordinates. The computational time per frequency point of the new computational scheme is a linear function of system size, a significant improvement over traditional, still-matrix techniques whose computational times per frequency point range from quadratic to cubic functions of system size. This permits the practical frequency domain analysis of systems of much larger order than by traditional, full-matrix techniques. Formulations are given for both open- and closed-loop systems. Numerical examples are presented showing the advantages of the present formulation over traditional approaches, both in speed and in accuracy. Using a model with 703 structural modes, the present method was up to two orders of magnitude faster than a traditional method. The present method generally showed good to excellent accuracy throughout the range of test frequencies, while traditional methods gave adequate accuracy for lower frequencies, but generally deteriorated in performance at higher frequencies with worst case errors being many orders of magnitude times the correct values.

Maghami, Peiman G.↗