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 73 records · Page 4

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.↗

Targeting Ballistic Lunar Capture Trajectories Using Periodic Orbits in the Sun-Earth CRTBP

A particular periodic orbit in the Earth-Sun circular restricted three body problem is shown to have the characteristics needed for a ballistic lunar capture transfer. An injection from a circular parking orbit into the periodic orbit serves as an initial guess for a targeting algorithm. By targeting appropriate parameters incrementally in increasingly complicated force models and using precise derivatives calculated from the state transition matrix, a reliable algorithm is produced. Ballistic lunar capture trajectories in restricted four body systems are shown to be able to be produced in a systematic way.

Cooley, D.S.↗

Radiative transfer in realistic planetary atmospheres

The research accomplished during this period is briefly summarized. The interior radiances within an optically deep absorbing medium scattering according to the Haze L phase function is discussed along with a method for calculating the radiance and color of the twilight sky. The application of the matrix operator method to calculations of radiance, polarization, and ellipticity of the radiation scattered from homogeneous layers scattering is reported. Reports, and publications are listed.

Plass, G. N.↗

Evaluation and development of new VAS remote sensing algorithms

An algorithm, developed to permit the simultaneous retrieval of temperature/moisture profiles and surface skin temperature by direct analytical solution of the radiative transfer equation, alleviates the problem associated with the interdependencies of water vapor retrieval on temperature retrieval and their dependencies on surface emissions. Simultaneous solution of all these quantities is achieved in one calculation using the available radiance observations. Since only a single matrix inversion is required for the specification of all parameters, the solution is computationally efficient. Ancillary observations of temperature and/or moisture from surface sensors on aircraft can be readily incorporated into the solution. Simulation tests of the method indicate improved performance over the previous iterative technique, particularly for the lower troposphere and for water vapor.

Smith, W. L.↗

Expansion of Tabulated Scattering Matrices in Generalized Spherical Functions

An efficient way to solve the vector radiative transfer equation for plane-parallel turbid media is to Fourier-decompose it in azimuth. This methodology is typically based on the analytical computation of the Fourier components of the phase matrix and is predicated on the knowledge of the coefficients appearing in the expansion of the normalized scattering matrix in generalized spherical functions. Quite often the expansion coefficients have to be determined from tabulated values of the scattering matrix obtained from measurements or calculated by solving the Maxwell equations. In such cases one needs an efficient and accurate computer procedure converting a tabulated scattering matrix into the corresponding set of expansion coefficients. This short communication summarizes the theoretical basis of this procedure and serves as the user guide to a simple public-domain FORTRAN program.

scattering matrix↗

Effect of ohmic, mass-transfer, and kinetic resistances on linear-sweep voltammetry in a cylindrical-pore electrode

A means is presented to account for the effect of ohmic, mass-transfer, and kinetic resistances on linear-sweep voltammograms by modeling a pore in a porous matrix as a cylindrical-pore electrode, and solving the mass and charge conservation equations in the context of this geometry for the simply redox reaction O + ne(-) yield R where both O and R are soluble species. Both analytical and numerical techniques are used to solve the governing equations. The calculated peak currents and potentials are correlated by empirical formulas to the measurable parameters: sweep rate, concentration of the redox species, diffusion coefficient, conductivity of the electrolyte, and pore dimensions. Using the correlations, a methodology is established for determining if the redox reaction kinetics are irreversible or reversible (Nernstian). If the reaction is irreversible, it is shown how the standard rate constant and the transfer coefficient may be extracted from linear-sweep voltammetry data, or, if the reaction is reversible, how the number of electrons transferred may be deduced.

Weidner, John W.↗

Modifying real convolutional codes for protecting digital filtering systems

A novel method is proposed for protecting digital filters from temporary and permanent failures that are not easily detected by conventional fault-tolerant computer design principles, on the basis of the error-detecting properties of real convolutional codes. Erroneous behavior is detected by externally comparing the calculated and regenerated parity samples. Great simplifications are obtainable by modifying the code structure to yield simplified parity channels with finite impulse response structures. A matrix equation involving the original parity values of the code and the polynomial of the digital filter's transfer function is formed, and row manipulations separate this equation into a set of homogeneous equations constraining the modifying scaling coefficients and another set which defines the code parity values' implementation.

Redinbo, G. R.↗

The updated statistical inversion technique to the evaluation of Umkehr observations

In the present study the standard retrieval Umkehr method to estimate the vertical distribution of ozone was updated using a statistical approach to the mathematical inversion scheme. The vertical ozone profile covariance matrix was used as a priori information for the inverse problem. A new method of the ozonesonde data organization according to air mass types helped to improve the covariance matrix quality. A retrieval method was developed using eigenvector technique. An optimal vertical ozone profile resolution was determined from the mathematical inversion scheme analysis based on the same technique. The sun radiation transfer was accounted for multiple scattering and atmospheric sphericity in this calculation. The retrievals using actual Umkehr Dobson spectrophotometer observations were also performed to provide the comparison of the standard and updated methods with concurrent ozone sound data at Boulder U.S. The comparison has revealed that the present method has some advantages in both resolution and accuracy, as compared to the standard one, especially for the atmospheric layers below ozone maximum.

Frolov, Alexander D.↗

Multi-Frequency NMM3D Simulations of Wave Propagation in Vegetation for Remote Sensing of Soil Moisture

To investigate the feasibility of using multi-frequencies to fur- ther improve the soil moisture retrieval, the recently devel- oped hybrid method is used to perform full-wave scattering simulations of a wheat field at three different frequencies of L-, S- and C- bands. In the hybrid method, the multiple scattering within a single wheat plant is first captured using the T-matrix of vec- tor cylindrical waves calculated from the full-wave solutions of HFSS. In the second step, the scattering among different plants is considered using the Foldy-Lax equations of multiple scattering theory (MST) and the extracted T-matrix. The transmission of microwaves through wheat field at the L-, S- and C-bands are calculated using the hybrid method. Results show that: (1) the transmission obtained from the full-wave simulations is much larger than those computed from the radiative transfer equations (RTE) model, (2) the hybrid method transmission has a weaker frequency dependence than those of the RTE, and (3) the attenuation caused by the vegetation layer would saturate with an increase in frequency.

Yueh, Simon↗

Minimum impulse three-body trajectories.

A rapid and accurate method of calculating optimal impulsive transfers in the restricted problem of three bodies has been developed. The technique combines a multi-conic method of trajectory integration with primer vector theory and an accelerated gradient method of trajectory optimization. A unique feature is that the state transition matrix and the primer vector are found analytical without additional integrations or differentiations. The method has been applied to the determination of optimal two and three impulse transfers between the L2 libration point and circular orbits about both the earth and the moon.

D'Amario, L.↗

Multiple scattered radiation emerging from Rayleigh and continental haze layers. I - Radiance, polarization, and neutral points

The matrix operator method was used to calculate the polarization of radiation scattered on layers of various optical thicknesses, with results compared for Rayleigh scattering and for scattering from a continental haze. In both cases, there are neutral points arising from the zeros of the polarization of single scattered photons at scattering angles of zero and 180 degrees. The angular position of these Rayleigh-like neutral points (RNP) in the sky shows appreciable variation with the optical thickness of the scattering layer for a Rayleigh phase matrix, but only a small variation for haze L phase matrix. Another type of neutral point exists for non-Rayleigh phase functions that is associated with the zeros of the polarization for single scattering which occurs between the end points of the curve. A comparison of radiances calculated from the complete theory of radiative transfer using Stokes vectors with those obtained from the scalar theory shows that differences of the order of 23% may be obtained for Rayleigh scattering, while the largest difference found for a haze L phase function was of the order of 0.1%.

Kattawar, G. W.↗

Calculation of the Frequency Distribution of the Energy Deposition in DNA Volumes by Heavy Ions

Radiation quality effects are largely determined by energy deposition in small volumes of characteristic sizes less than 10 nm representative of short-segments of DNA, the DNA nucleosome, or molecules initiating oxidative stress in the nucleus, mitochondria, or extra-cellular matrix. On this scale, qualitatively distinct types of molecular damage are possible for high linear energy transfer (LET) radiation such as heavy ions compared to low LET radiation. Unique types of DNA lesions or oxidative damages are the likely outcome of the energy deposition. The frequency distribution for energy imparted to 1-20 nm targets per unit dose or particle fluence is a useful descriptor and can be evaluated as a function of impact parameter from an ions track. In this work, the simulation of 1-Gy irradiation of a cubic volume of 5 micron by: 1) 450 (1)H(+) ions, 300 MeV; 2) 10 (12)C(6+) ions, 290 MeV/amu and 3) (56)Fe(26+) ions, 1000 MeV/amu was done with the Monte-Carlo simulation code RITRACKS. Cylindrical targets are generated in the irradiated volume, with random orientation. The frequency distribution curves of the energy deposited in the targets is obtained. For small targets (i.e. <25 nm size), the probability of an ion to hit a target is very small; therefore a large number of tracks and targets as well as a large number of histories are necessary to obtain statistically significant results. This simulation is very time-consuming and is difficult to perform by using the original version of RITRACKS. Consequently, the code RITRACKS was adapted to use multiple CPU on a workstation or on a computer cluster. To validate the simulation results, similar calculations were performed using targets with fixed position and orientation, for which experimental data are available [5]. Since the probability of single- and double-strand breaks in DNA as function of energy deposited is well know, the results that were obtained can be used to estimate the yield of DSB, and can be extended to include other targeted or non-target effects.

Plante, Ianik↗

Computer programs for calculation of matrix stability and frequency response from a state-space system description

FORTRAN computer subroutines stemming from requirements to process state variable system equations for systems of high order are presented. They find the characteristic equation of a matrix using the method of Danilevsky, the number of roots with positive real parts using the Routh-Horwitz alternate formulation, convert a state variable system description to a Laplace transfer function using the method of Bollinger, and evaluate that transfer function and obtain its frequency response. A sample problem is presented to demonstrate use of the subroutines.

Seidel, R. C.↗

Single and multiple scattering contributions to circumsolar radiation

The contributions to the angular distribution of the almucantar radiance in the forward direction due to multiple scattering are compared to those due to single scattering. The contributions have been calculated by a computer code employing the Gauss-Seidel iterative approach to the solution of the radiative transfer equation for a plane parallel atmosphere composed of air molecules, aerosol particles, and ozone. The code is similar to that of Dave (1972) except in the construction of the source matrix. In the near-forward direction the multiple scattering contributions are significant for optical depths of the order of 0.4. The shape of the angular distribution of almucantar radiance to 10 degrees is less sensitive to multiple scattering.

Box, M. A.↗

Pressure transfer function of a JT15D nozzle due to acoustic and convected entropy fluctuations

An acoustic transmission matrix analysis of sound propagation in a variable area duct with and without flow is extended to include convected entropy fluctuations. The boundary conditions used in the analysis are a transfer function relating entropy and pressure at the nozzle inlet and the nozzle exit impedance. The nozzle pressure transfer function calculated is compared with JT15D turbofan engine nozzle data. The one dimensional theory for sound propagation in a variable area nozzle with flow but without convected entropy is good at the low engine speeds where the nozzle exit Mach number is low (M=0.2) and the duct exit impedance model is good. The effect of convected entropy appears to be so negligible that it is obscured by the inaccuracy of the nozzle exit impedance model, the lack of information on the magnitude of the convected entropy and its phase relationship with the pressure, and the scatter in the data. An improved duct exit impedance model is required at the higher engine speeds where the nozzle exit Mach number is high (M=0.56) and at low frequencies (below 120 Hz).

Miles, J. H.↗

Recent advances at NASA in calculating the electronic spectra of diatomic molecules

Advanced entry vehicles, such as the proposed Aero-assisted Orbital Transfer Vehicle, provide new and challenging problems for spectroscopy. Large portions of the flow field about such vehicles will be characterized by chemical and thermal nonequilibrium. Only by considering the actual overlap of the atomic and rotational lines emitted by the species present can the impact of radiative transport within the flow field be assessed correctly. To help make such an assessment, a new computer program is described that can generate high-resolution, line-by-line spectra for any spin-allowed transitions in diatomic molecules. The program includes the matrix elements for the rotational energy and distortion to the fourth order; the spin-orbit, spin-spin, and spin-rotation interactions to first order; and the lambda splitting by a perturbation calculation. An overview of the Computational Chemistry Branch at Ames Research Center is also presented.

Whiting, Ellis E.↗

Analytical and experimental study of vibrations in a gear transmission

An analytical simulation of the dynamics of a gear transmission system is presented and compared to experimental results from a gear noise test rig at the NASA Lewis Research Center. The analytical procedure developed couples the dynamic behaviors of the rotor-bearing-gear system with the response of the gearbox structure. The modal synthesis method is used in solving the overall dynamics of the system. Locally each rotor-gear stage is modeled as an individual rotor-bearing system using the matrix transfer technique. The dynamics of each individual rotor are coupled with other rotor stages through the nonlinear gear mesh forces and with the gearbox structure through bearing support systems. The modal characteristics of the gearbox structure are evaluated using the finite element procedure. A variable time steping integration routine is used to calculate the overall time transient behavior of the system in modal coordinates. The global dynamic behavior of the system is expressed in a generalized coordinate system. Transient and steady state vibrations of the gearbox system are presented in the time and frequency domains. The vibration characteristics of a simple single mesh gear noise test rig is modeled. The numerical simulations are compared to experimental data measured under typical operating conditions. The comparison of system natural frequencies, peak vibration amplitudes, and gear mesh frequencies are generally in good agreement.

Choy, F. K.↗