Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Matrix elements”

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 235 records · Page 13

Locating and Quantifying Broadband Fan Sources Using In-Duct Microphones

In-duct beamforming techniques have been developed for locating broadband noise sources on a low-speed fan and quantifying the acoustic power in the inlet and aft fan ducts. The NASA Glenn Research Center's Advanced Noise Control Fan was used as a test bed. Several of the blades were modified to provide a broadband source to evaluate the efficacy of the in-duct beamforming technique. Phased arrays consisting of rings and line arrays of microphones were employed. For the imaging, the data were mathematically resampled in the frame of reference of the rotating fan. For both the imaging and power measurement steps, array steering vectors were computed using annular duct modal expansions, selected subsets of the cross spectral matrix elements were used, and the DAMAS and CLEAN-SC deconvolution algorithms were applied.

Dougherty, Robert P.↗

Graph Theory Roots of Spatial Operators for Kinematics and Dynamics

Spatial operators have been used to analyze the dynamics of robotic multibody systems and to develop novel computational dynamics algorithms. Mass matrix factorization, inversion, diagonalization, and linearization are among several new insights obtained using such operators. While initially developed for serial rigid body manipulators, the spatial operators and the related mathematical analysis have been shown to extend very broadly including to tree and closed topology systems, to systems with flexible joints, links, etc. This work uses concepts from graph theory to explore the mathematical foundations of spatial operators. The goal is to study and characterize the properties of the spatial operators at an abstract level so that they can be applied to a broader range of dynamics problems. The rich mathematical properties of the kinematics and dynamics of robotic multibody systems has been an area of strong research interest for several decades. These properties are important to understand the inherent physical behavior of systems, for stability and control analysis, for the development of computational algorithms, and for model development of faithful models. Recurring patterns in spatial operators leads one to ask the more abstract question about the properties and characteristics of spatial operators that make them so broadly applicable. The idea is to step back from the specific application systems, and understand more deeply the generic requirements and properties of spatial operators, so that the insights and techniques are readily available across different kinematics and dynamics problems. In this work, techniques from graph theory were used to explore the abstract basis for the spatial operators. The close relationship between the mathematical properties of adjacency matrices for graphs and those of spatial operators and their kernels were established. The connections hold across very basic requirements on the system topology, the nature of the component bodies, the indexing schemes, etc. The relationship of the underlying structure is intimately connected with efficient, recursive computational algorithms. The results provide the foundational groundwork for a much broader look at the key problems in kinematics and dynamics. The properties of general graphs and trees of nodes and edge were examined, as well as the properties of adjacency matrices that are used to describe graph connectivity. The nilpotency property of such matrices for directed trees was reviewed, and the adjacency matrices were generalized to the notion of block weighted adjacency matrices that support block matrix elements. This leads us to the development of the notion of Spatial Kernel Operator SKO kernels. These kernels provide the basis for the development of SKO resolvent operators.

Jain, Abhinandan↗

A Multiple Sphere T-Matrix Fortran Code for Use on Parallel Computer Clusters

A general-purpose Fortran-90 code for calculation of the electromagnetic scattering and absorption properties of multiple sphere clusters is described. The code can calculate the efficiency factors and scattering matrix elements of the cluster for either fixed or random orientation with respect to the incident beam and for plane wave or localized- approximation Gaussian incident fields. In addition, the code can calculate maps of the electric field both interior and exterior to the spheres.The code is written with message passing interface instructions to enable the use on distributed memory compute clusters, and for such platforms the code can make feasible the calculation of absorption, scattering, and general EM characteristics of systems containing several thousand spheres.

Mackowski, D. W.↗

Spectral Line Parameters Including Temperature Dependences of Self- and Air-Broadening in the 2 (left arrow) 0 Band of CO at 2.3 micrometers

Temperature dependences of pressure-broadened half-width and pressure-induced shift coefficients along with accurate positions and intensities have been determined for transitions in the 2<--0 band of C-12 O-16 from analyzing high-resolution and high signal-to-noise spectra recorded with two different Fourier transform spectrometers. A total of 28 spectra, 16 self-broadened and 12 air-broadened, recorded using high- purity (greater than or equal to 99.5% C-12-enriched) CO samples and CO diluted with dry air(research grade) at different temperatures and pressures, were analyzed simultaneously to maximize the accuracy of the retrieved parameters. The sample temperatures ranged from 150 to 298K and the total pressures varied between 5 and 700 Torr. A multispectrum nonlinear least squares spectrum fitting technique was used to adjust the rovibrational constants (G, B, D, etc.) and intensity parameters (including Herman-Wallis coefficients), rather than determining individual line positions and intensities. Self-and air-broadened Lorentz half-width coefficients, their temperature dependence exponents, self- and air-pressure-induced shift coefficients, their temperature dependences, self- and air-line mixing coefficients, their temperature dependences and speed dependence have been retrieved from the analysis. Speed-dependent line shapes with line mixing employing off-diagonal relaxation matrix element formalism were needed to minimize the fit residuals. This study presents a precise and complete set of spectral line parameters that consistently reproduce the spectrum of carbon monoxide over terrestrial atmospheric conditions.

Devi, V. Malathy↗

Bound-free Spectra for Diatomic Molecules

It is now recognized that prediction of radiative heating of entering space craft requires explicit treatment of the radiation field from the infrared (IR) to the vacuum ultra violet (VUV). While at low temperatures and longer wavelengths, molecular radiation is well described by bound-bound transitions, in the short wavelength, high temperature regime, bound-free transitions can play an important role. In this work we describe first principles calculations we have carried out for bound-bound and bound-free transitions in N2, O2, C2, CO, CN, NO, and N2+. Compared to bound ]bound transitions, bound-free transitions have several particularities that make them different to deal with. These include more complicated line shapes and a dependence of emission intensity on both bound state diatomic and atomic concentrations. These will be discussed in detail below. The general procedure we used was the same for all species. The first step is to generate potential energy curves, transition moments, and coupling matrix elements by carrying out ab initio electronic structure calculations. These calculations are expensive, and thus approximations need to be made in order to make the calculations tractable. The only practical method we have to carry out these calculations is the internally contracted multi-reference configuration interaction (icMRCI) method as implemented in the program suite Molpro. This is a widely used method for these kinds of calculations, and is capable of generating very accurate results. With this method, we must first of choose which electrons to correlate, the one-electron basis to use, and then how to generate the molecular orbitals.

Schwenke, David W.↗

T-Matrix Modeling of Linear Depolarization by Morphologically Complex Soot and Soot-Containing Aerosols

We use state-of-the-art public-domain Fortran codes based on the T-matrix method to calculate orientation and ensemble averaged scattering matrix elements for a variety of morphologically complex black carbon (BC) and BC-containing aerosol particles, with a special emphasis on the linear depolarization ratio (LDR). We explain theoretically the quasi-Rayleigh LDR peak at side-scattering angles typical of low-density soot fractals and conclude that the measurement of this feature enables one to evaluate the compactness state of BC clusters and trace the evolution of low-density fluffy fractals into densely packed aggregates. We show that small backscattering LDRs measured with groundbased, airborne, and spaceborne lidars for fresh smoke generally agree with the values predicted theoretically for fluffy BC fractals and densely packed near-spheroidal BC aggregates. To reproduce higher lidar LDRs observed for aged smoke, one needs alternative particle models such as shape mixtures of BC spheroids or cylinders.

atmospheric radiation↗

Two Dimensional Symmetric Correlation Functions of the S Operator and Two Dimensional Fourier Transforms: Considering the Line Coupling for P and R Lines of Linear Molecules

The refinement of the Robert-Bonamy (RB) formalism by considering the line coupling for isotropic Raman Q lines of linear molecules developed in our previous study [Q. Ma, C. Boulet, and R. H. Tipping, J. Chem. Phys. 139, 034305 (2013)] has been extended to infrared P and R lines. In these calculations, the main task is to derive diagonal and off-diagonal matrix elements of the Liouville operator iS1 − S2 introduced in the formalism. When one considers the line coupling for isotropic Raman Q lines where their initial and final rotational quantum numbers are identical, the derivations of off-diagonal elements do not require extra correlation functions of the ^S operator and their Fourier transforms except for those used in deriving diagonal elements. In contrast, the derivations for infrared P and R lines become more difficult because they require a lot of new correlation functions and their Fourier transforms. By introducing two dimensional correlation functions labeled by two tensor ranks and making variable changes to become even functions, the derivations only require the latters' two dimensional Fourier transforms evaluated at two modulation frequencies characterizing the averaged energy gap and the frequency detuning between the two coupled transitions. With the coordinate representation, it is easy to accurately derive these two dimensional correlation functions. Meanwhile, by using the sampling theory one is able to effectively evaluate their two dimensional Fourier transforms. Thus, the obstacles in considering the line coupling for P and R lines have been overcome. Numerical calculations have been carried out for the half-widths of both the isotropic Raman Q lines and the infrared P and R lines of C2H2 broadened by N2. In comparison with values derived from the RB formalism, new calculated values are significantly reduced and become closer to measurements.

infrared radiation↗

Effects on Calculated Half-Widths and Shifts from the Line Coupling for Asymmetric-Top Molecules

The refinement of the Robert-Bonamy formalism by considering the line coupling for linear molecules developed in our previous studies [Q. Ma, C. Boulet, and R. H. Tipping, J. Chem. Phys. 139, 034305 (2013); 140, 104304 (2014)] have been extended to asymmetric-top molecules. For H2O immersed in N2 bath, the line coupling selection rules applicable for the pure rotational band to determine whether two specified lines are coupled or not are established. Meanwhile, because the coupling strengths are determined by relative importance of off-diagonal matrix elements versus diagonal elements of the operator -iS1 -S2, quantitative tools are developed with which one is able to remove weakly coupled lines from consideration. By applying these tools, we have found that within reasonable tolerances, most of the H2O lines in the pure rotational band are not coupled. This reflects the fact that differences of energy levels of the H2O states are pretty large. But, there are several dozen strongly coupled lines and they can be categorized into different groups such that the line couplings occur only within the same groups. In practice, to identify those strongly coupled lines and to confine them into sub-linespaces are crucial steps in considering the line coupling. We have calculated half-widths and shifts for some groups, including the line coupling. Based on these calculations, one can conclude that for most of the H2O lines, it is unnecessary to consider the line coupling. However, for several dozens of lines, effects on the calculated half-widths from the line coupling are small, but remain noticeable and reductions of calculated half-widths due to including the line coupling could reach to 5%. Meanwhile, effects on the calculated shifts are very significant and variations of calculated shifts could be as large as 25%.

formalism↗

Air- and Self-Broadened Half Widths, Pressure-Induced Shifts, and Line Mixing in the Nu(sub 2) Band of (12)CH4

Lorentz self- and air-broadened half width and pressure-induced shift coefficients and their dependences on temperature have been measured from laboratory absorption spectra for nearly 130 transitions in the nu(sub 2) band of (12)CH4. In addition line mixing coefficients (using the relaxation matrix element formalism) for both self- and airbroadening were experimentally determined for the first time for a small number of transitions in this band. Accurate line positions and absolute line intensities were also determined. These parameters were obtained by analyzing high-resolution (approx. 0.003 to 0.01 per cm) laboratory spectra of high-purity natural CH4 and air-broadened CH4 recorded at temperatures between 226 and 297 K using the McMath-Pierce Fourier transform spectrometer (FTS) located at the National Solar Observatory on Kitt Peak, Arizona. A multispectrum nonlinear least squares technique was used to fit short (5-15 per cm) spectral intervals in 24-29 spectra simultaneously. Parameters were determined for nu(sub 2) transitions up to J" = 16. The variations of the measured broadening and shift parameters with the rotational quantum number index and tetrahedral symmetry species are examined. The present results are also compared with previous measurements available in the literature.

Smith, M. A. H.↗

Estimating Terra MODIS Polarization Effect Using Ocean Data

Terra MODIS has been known since pre-launch to have polarization sensitivity, particularly in shortest-wavelength bands 8 and 9. On-orbit reflectance trending of pseudo-invariant sites show a variation in reflectance as a function of band and scan mirror angle of incidence consistent with time-dependent polarization effects from the rotating double-sided scan mirror. The MODIS Characterization Support Team [MCST] estimates the Mueller matrix trending from this variation as observed from a single desert site, but this effect is not included in Collection 6 [C6] calibration. Here we extend the MCSTs current polarization sensitivity monitoring to two ocean sites distributed over latitude to helpestimate the uncertainties in the derived Mueller matrix. The Mueller matrix elements derived for polarization-sensitive Band 8 for a given site are found to be fairly insensitive to surface brdf modeling. The site-to-site variation is a measure of the uncertainty in the Mueller estimation.Results for band 8 show that the polarization correction reduces mirror-side striping by up to 50% and reduces the instrument polarization effect on reflectance time series of an ocean target.

Wald, Andrew E.↗

Self- and Air-Broadened Line Shapes in the 2v3 P and R Branches of 12CH4

In this paper we report line shape parameters of 12CH4 for several hundred 2V(sub 3) transitions in the spectral regions 5891-5996 cm( exp -1) (P branch) and 6015-6115 cm(exp -1) (R branch). Air- and self-broadening coefficients were measured as a function of temperature; line mixing via off-diagonal relaxation matrix element coefficients was also obtained for 47 transition pairs. In total, nearly 1517 positions and intensities were retrieved, but many transitions were too weak for the line shape study. For this analysis, we used 25 high-resolution (0.0056 and 0.0067 cm(ex[ -1) and high signal-to-noise (S/N) spectra of high-purity 12CH4 and the same high-purity 12CH4 broadened by dry air recorded at different sample temperatures between 130 K and 295 K with the Bruker IFS 125HR Fourier transform spectrometer at JPL. Three different absorption cells were used (1) a White cell set to a path length of 13.09 m for room temperature data, (2) a single-pass 0.2038 m long coolable cell (for self-broadening) and (3) a multipass cell with 20.941 m total path coolable Herriott cell (for air-broadening). In total there were 13 spectra with pure 12CH4 (0.27-599 Torr) and 12 air-broadened spectra with total sample pressures of 80-805 Torr and volume mixing ratios (VMR) of methane between 0.18 and 1.0. An interactive multispectrum nonlinear least-squares technique was employed to fit the individual P10-P1 and R0-R10 manifolds in all the spectra simultaneously. Results obtained from the present analysis are compared to other recent measurements.

Devi, V. Malathy↗

The Effect of the Spin-Forbidden Co((sup 1) Sigma plus) plus O((sup 3) P) Yields CO2 (1 Sigma (sub G) plus) Recombination Reaction on Afterbody Heating of Mars Entry Vehicles

Vibrationally excited CO2, formed by two-body recombination from CO((sup 1) sigma plus) and O((sup 3) P) in the wake behind spacecraft entering the Martian atmosphere reaction, is potentially responsible for the higher than anticipated radiative heating of the backshell, compared to pre-flight predictions. This process involves a spin-forbidden transition of the transient triplet CO2 molecule to the longer-lived singlet. To accurately predict the singlet-triplet transition probability and estimate the thermal rate coefficient of the recombination reaction, ab initio methods were used to compute the first singlet and three lowest triplet CO2 potential energy surfaces and the spin-orbit coupling matrix elements between these states. Analytical fits to these four potential energy surfaces were generated for surface hopping trajectory calculations, using Tully's fewest switches surface hopping algorithm. Preliminary results for the trajectory calculations are presented. The calculated probability of a CO((sup 1) sigma plus) and O((sup 3) P) collision leading to singlet CO2 formation is on the order of 10 (sup -4). The predicted flowfield conditions for various Mars entry scenarios predict temperatures in the range of 1000 degrees Kelvin - 4000 degrees Kelvin and pressures in the range of 300-2500 pascals at the shoulder and in the wake, which is consistent with a heavy-particle collision frequency of 10 (sup 6) to 10 (sup 7) per second. Owing to this low collision frequency, it is likely that CO((sup 1) sigma plus) molecules formed by this mechanism will mostly be frozen in a highly nonequilibrium rovibrational energy state until they relax by photoemission.

CO2 recombination↗

Line shape parameters of PH3 transitions: Theoretical studies of self-broadened widths and line mixing effects

Line mixing effects have been calculated in various parallel and perpendicular bands of self-broadened PH3 lines and compared with recent experimental data. The theoretical approach is an extension to symmetric tops with high inversion barrier of the formalism previously developed for NH3 [Q. Ma and C. Boulet, J. Chem. Phys. 144, 224303 (2016)]. The model takes into account the non-diagonality of the scattering operator within the line space as well as, in a correct way, the double degeneracy of the j, k levels when k ≠ 0. Transitions between such levels should be considered as doublets whose components may be coupled by the line mixing process. It has been shown that, at low pressure, the inversion of the experimental data will strongly depend on the splitting between the two components of a doublet. When it is significant, one can measure independently both the width of one component and the intra-doublet coupling matrix element. Otherwise, one can only measure the sum of these two elements. Comparisons with measurements show that the present formalism leads to accurate predictions of the experimental line shapes.

Remote sensing↗

Theoretical Study of CH3Cl-N2 Line Shapes in the ν1 Band. Line Mixing Effects in QR Doublets and QQk Sub-Branches

Line coupling and line mixing effects of CH3Cl lines in the ν1 band perturbed by N2 have been investigated. We have taken into account the non-diagonality of the exp(−S2) operator within specified line spaces as well as the k-degeneracy of the transitions (due to the double degeneracy of the j,k levels with k ≠ 0). These transitions should be considered as doublets coupled by the line mixing process. A new problem appears in the calculation when the atom-atom potential model is introduced. In order to overcome this difficulty, a pragmatic approach is proposed. Comparisons of theoretically calculated matrix elements of W with measurements of QR(j,k) doublets as well as some QQk sub-branches, which are strongly affected by line mixing, have been made. The results show that the formalism improved in this way leads to rather accurate predictions.

Methyl chloride↗

Simultaneous Stoquasticity

Stoquastic Hamiltonians play a role in the computational complexity of the local Hamiltonian problem as well as the study of classical simulability. In particular, stoquastic Hamiltonians can be straightforwardly simulated using Monte Carlo techniques. We address the question of whether two or more Hamiltonians may be made simultaneously stoquastic via a unitary transformation. This question has important implications for the complexity of simulating quantum annealing where quantum advantage is related to the stoquasticity of the Hamiltonians involved in the anneal. We find that for almost all problems no such unitary exists and show that the problem of determining the existence of such a unitary is equivalent to identifying if there is a solution to a system of polynomial (in)equalities in the matrix elements of the initial and transformed Hamiltonians. Solving such a system of equations is NP-hard. We highlight a geometric understanding of this problem in terms of a collection of generalized Bloch vectors.

Jacob Bringewatt↗

Simultaneous Stoquasticity

Stoquastic Hamiltonians play a role in the computational complexity of the local Hamiltonian problem as well as the study of classical simulability. In particular, stoquastic Hamiltonians can be straightforwardly simulated using Monte Carlo techniques. We address the question of whether two or more Hamiltonians may be made simultaneously stoquastic via a unitary transformation. This question has important implications for the complexity of simulating quantum annealing where quantum advantage is related to the stoquasticity of the Hamiltonians involved in the anneal. We find that for almost all problems no such unitary exists and show that the problem of determining the existence of such a unitary is equivalent to identifying if there is a solution to a system of polynomial (in)equalities in the matrix elements of the initial and transformed Hamiltonians. Solving such a system of equations is NP-hard. We highlight a geometric understanding of this problem in terms of a collection of generalized Bloch vectors.

Monte Carlo↗

Electron-H2 Collisions Studied Using the Finite Element Z-Matrix Method

We have applied the Z-matrix method, using a mixed basis of finite elements and Gaussians, to study e-H2 elastic and inelastic collisions. Special attention is paid to the quality of the basis set and the treatment of electron correlation. The calculated cross sections are invariant, to machine accuracy, with respect to the choice of parameters a, b, d, e as long as they satisfy Equation (3). However, the log derivative approach, i.e., the choice a = -e = 1, b = d = 0 appears to converge slightly faster than other choices. The cross sections agree well with previous theoretical results. Comparison will be made with available experimental data.

Huo, Winifred M.↗