Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Derived Matrices”

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 325 records · Page 18

[PAH Cations as Viable Carriers of DIBs]

This report is intended to fill in the blanks in NASA's file system for our lab astro study of molecular ions of astrophysical interest. In order to give NASA what it needs for its files, I attach below the text of the section from our recent proposal to continue this work, in which we describe progress to date, including a large number of publications. Our initial studies were focused on PAH cations, which appear to be viable candidates as the carriers of the DIBs, an idea that has been supported by laboratory spectroscopy of PAH cations in inert matrices. Beginning with the simplest aromatic (benzene; C6H6) and moving progressively to larger species (naphthalene, C10OH8; pyrene, C16H10; and most recently chrysene, C18H12), we have been able to derive rate coefficients for reactions with neutral spices that are abundant in the diffuse interstellar medium.

Snow, Ted↗

Improved solution for ill-conditioned algebraic equations by epsilon decomposition

Matrix eigenvalue theory is presently used to examine the source of ill-conditioning in linear algebraic equations; the approach highlights the critical role played by the zero and near-zero eigenvalues and corresponding eigenvectors of poorly-conditioned systems. Insights derived from this approach are used to improve the recently developed epsilon-decomposition (E-D) solution procedure. The efficiency of E-D is significant for large matrices possessing small rank deficiency.

Ojalvo, Irving U.↗

Highlights of TOMS Version 9 Total Ozone Algorithm

The fundamental basis of TOMS total ozone algorithm was developed some 45 years ago by Dave and Mateer. It was designed to estimate total ozone from satellite measurements of the backscattered UV radiances at few discrete wavelengths in the Huggins ozone absorption band (310-340 nm). Over the years, as the need for higher accuracy in measuring total ozone from space has increased, several improvements to the basic algorithms have been made. They include: better correction for the effects of aerosols and clouds, an improved method to account for the variation in shape of ozone profiles with season, latitude, and total ozone, and a multi-wavelength correction for remaining profile shape errors. These improvements have made it possible to retrieve total ozone with just 3 spectral channels of moderate spectral resolution (approx. 1 nm) with accuracy comparable to state-of-the-art spectral fitting algorithms like DOAS that require high spectral resolution measurements at large number of wavelengths. One of the deficiencies of the TOMS algorithm has been that it doesn't provide an error estimate. This is a particular problem in high latitudes when the profile shape errors become significant and vary with latitude, season, total ozone, and instrument viewing geometry. The primary objective of the TOMS V9 algorithm is to account for these effects in estimating the error bars. This is done by a straightforward implementation of the Rodgers optimum estimation method using a priori ozone profiles and their error covariances matrices constructed using Aura MLS and ozonesonde data. The algorithm produces a vertical ozone profile that contains 1-2.5 pieces of information (degrees of freedom of signal) depending upon solar zenith angle (SZA). The profile is integrated to obtain the total column. We provide information that shows the altitude range in which the profile is best determined by the measurements. One can use this information in data assimilation and analysis. A side benefit of this algorithm is that it is considerably simpler than the present algorithm that uses a database of 1512 profiles to retrieve total ozone. These profiles are tedious to construct and modify. Though conceptually similar to the SBUV V8 algorithm that was developed about a decade ago, the SBUV and TOMS V9 algorithms differ in detail. The TOMS algorithm uses 3 wavelengths to retrieve the profile while the SBUV algorithm uses 6-9 wavelengths, so TOMS provides less profile information. However both algorithms have comparable total ozone information and TOMS V9 can be easily adapted to use additional wavelengths from instruments like GOME, OMI and OMPS to provide better profile information at smaller SZAs. The other significant difference between the two algorithms is that while the SBUV algorithm has been optimized for deriving monthly zonal means by making an appropriate choice of the a priori error covariance matrix, the TOMS algorithm has been optimized for tracking short-term variability using month and latitude dependent covariance matrices.

Bhartia, Pawan↗

Estimation of full modal damping matrices from complex test modes

This paper describes the refinements of a previously published method for estimating a full modal damping matrix from complex test modes. It also documents application of the refined method to a structure where complex test modes were derived by the ERA method from multi-input random vibration test data. A numerical example based on simulated test data is presented to demonstrate the validity of the method. The application using real data was not successful, presumably because of noise in the small phase angles of the measured complex modes. Alternative test and data reduction procedures are suggested as possible remedies to the problem. A careful analysis of measurement and data processing errors should be made to examine basic feasibility before implementing the alternative procedures. The ability to estimate a full modal damping matrix is considered important for the preflight estimation of on-orbit damping, and for the synthesis of structural damping from substructure tests.

Hasselman, T. K.↗

Classifying topology in photonic crystal slabs with radiative environments

Abstract In the recent years, photonic Chern materials have attracted substantial interest as they feature topological edge states that are robust against disorder, promising to realize defect-agnostic integrated photonic crystal slab devices. However, the out-of-plane radiative losses in those photonic Chern slabs has been previously neglected, yielding limited accuracy for predictions of these systems’ topological protection. Here, we develop a general framework for measuring the topological protection in photonic systems, such as in photonic crystal slabs, while accounting for in-plane and out-of-plane radiative losses. Our approach relies on the spectral localizer that combines the position and Hamiltonian matrices of the system to draw a real-picture of the system’s topology. This operator-based approach to topology allows us to use an effective Hamiltonian directly derived from the full-wave Maxwell equations after discretization via finite-elements method (FEM), resulting in the full account of all the system’s physical processes. As the spectral FEM-localizer is constructed solely from FEM discretization of the system’s master equation, the proposed framework is applicable to any physical system and is compatible with commonly used FEM software. Moving forward, we anticipate the generality of the method to aid in the topological classification of a broad range of complex physical systems.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

N -representability violations in truncated equation-of-motion coupled-cluster methods

One-electron reduced density matrices (1RDMs) from equation-of-motion (EOM) coupled-cluster with single and double excitations (CCSD) calculations are analyzed to assess their N-representability (i.e., whether they are derivable from a physical N-electron state). We identify EOM-CCSD stationary states whose 1RDMs violate either ensemble-state N-representability conditions or pure-state conditions known as generalized Pauli constraints. As such, these 1RDMs do not correspond to any physical N-electron state. Unphysical states are also encountered in the course of time-dependent EOM-CC simulations; when an external field drives transitions between a pair of stationary states with pure-state N-representable 1RDMs, the 1RDM of the time-dependent state can violate ensemble-state conditions. Furthermore, these observations point to potential challenges in interpreting the results of time-dependent EOM-CCSD simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Windowed multipole representation of R -matrix cross sections

Nuclear cross sections are basic inputs to any nuclear computation. Campaigns of experiments are fitted with the parametric R-matrix model of quantum nuclear interactions, and the resulting cross sections are documented—both pointwise and as resonance parameters (with uncertainties)—in standard evaluated nuclear data libraries (ENDF, JEFF, BROND, JENDL, CENDL, TENDL): these constitute our common knowledge of fundamental low-energy nuclear cross sections. In the past decade, a collaborative effort has been deployed to establish a new nuclear cross-section library format—the Windowed Multipole Library—with the goal of considerably reducing the computational cost of cross-section calculations in nuclear transport simulations. This work lays the theoretical foundations underpinning these efforts. From general R-matrix scattering theory, we derive the windowed multipole representation of nuclear cross sections. Though physically and mathematically equivalent to R-matrix cross sections, the windowed multipole representation is particularly well suited for subsequent temperature treatment of angle-integrated cross sections, in particular Doppler broadening, which is the averaging of cross sections over the thermal motion of the target atoms. Doppler broadening is of critical importance in neutron transport applications, as it ensures the stability of many nuclear reactors (negative thermal reactivity). Yet, Doppler broadening of nuclear cross sections has been a considerable bottleneck for nuclear transport computations, often requiring memory-costly pretabulations. We show that the windowed multipole representation can perform accurate Doppler broadening analytically (up to the first reaction threshold), from which we derive cross-section temperature derivatives to any order—all computable on the fly (without precalculations stored in memory). Furthermore, we here establish a way of converting the R-matrix resonance parameters uncertainty (covariance matrices) into windowed multipole parameters uncertainty. We show that generating stochastic nuclear cross sections by sampling from the resulting windowed multipole covariance matrix can reproduce the cross-section uncertainty in the original nuclear data file. The windowed multipole representation is therefore a novel nuclear physics formalism able to generate Doppler broadened stochastic nuclear cross sections on the fly, unlocking breakthrough computational gains for nuclear computations. Through this foundational paper, we hope to make the windowed multipole representation accessible, reproducible, and usable for the nuclear physics community, as well as provide the theoretical basis for future research on expanding its capabilities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Synthesis of improved phenolic resins

Twenty seven addition cured phenolic resin compositions were prepared and tested for their ability to give char residues comparable to state-of-the-art phenolic resins. Cyanate, epoxy, allyl, acrylate, methacrylate and ethynyl derivatized phenolic oligomers were investigated. The novolac-cyanate and propargyl-novolac resins provided anaerobic char yields at 800 C of 58 percent. A 59 percent char yield was obtained from modified epoxy novolacs. A phosphonitrilic derivative was found to be effective as an additive for increasing char yields. The novolac-cyanate, epoxy-novolac and methacrylate-epoxy-novolac systems were investigated as composite matrices with Thornel 300 graphite fiber. All three resins showed good potential as composite matrices. The free radical cured methacrylate-epoxy-novolac graphite composite provided short beam shear strengths at room temperature of 93.3 MPa (13.5 ksi). The novolac-cyanate graphite composite produced a short beam shear strength of 74 MPa (10.7 ksi) and flexural strength of 1302 MPa (189 ksi) at 177 C. Air heat aging of the novolac-cyanate and epoxy novolac based composites for 12 weeks at 204 C showed good property retention.

Delano, C. B.↗

Ada Namelist Package

Ada Namelist Package, developed for Ada programming language, enables calling program to read and write FORTRAN-style namelist files. Features are: handling of any combination of types defined by user; ability to read vectors, matrices, and slices of vectors and matrices; handling of mismatches between variables in namelist file and those in programmed list of namelist variables; and ability to avoid searching entire input file for each variable. Principle benefits derived by user: ability to read and write namelist-readable files, ability to detect most file errors in initialization phase, and organization keeping number of instantiated units to few packages rather than to many subprograms.

Klumpp, Allan R.↗

The direct effects of gravity on the control and output matrices of controlled structure models

The effects of gravity on the dynamic performance of structural control actuators and sensors are dual forms of an additive perturbation that can attenuate or amplify the device response (input or output). The modal modeling of these perturbations is derived for the general case of arbitrarily oriented devices and arbitrarily oriented planes of deformation. A nondimensional sensitivity analysis to identify the circumstances under which the effects of gravity are important is presented. Results show that gravity effects become important when the product of the ratio of the normalized modal slope and the modal displacement is comparable to the ratio of the gravitational acceleration and the product of the beam length and the squared eigenfrequency for a given mode.

Rey, Daniel A.↗

Hybrid Semiloof elements for plates and shells based upon a modified Hu-Washizu principle

Hybrid SemiLoof elements for plates and shells are developed based upon modified Hu-Washizu principle. In the new version of the assumed stress hybrid formulation the equilibrium equations are satisfied through the introduction of internal displacement parameters as Lagrange multipliers. The inversion of the resulting H-matrices is simplified particularly when the stresses are expressed in terms of natural coordinates. A 24-DOF triangular element and a 32-DOF quadrilateral element based on shallow shell theory are derived and evaluated.

Pian, T. H. H.↗

Derivation of phase statistics from the Mueller matrix

Experimental observations show a strong dependence of phase differences of scattering matrix elements on the physical parameters of random media. Here, the statistical behavior of the phase differences is studied for distributed targets. The pdfs of the phase differences are derived from the Mueller matrix of the target. In deriving the density functions, it is assumed that the real and imaginary parts of the copolarized and cross-polarized terms of the scattering matrix are jointly Gaussian and their covariance matrices are found in terms of the Mueller matrix elements. The functional forms of the copolarized and cross-polarized density functions are similar and are obtained independently. It is shown that the density function of the phase difference is completely determined in terms of only two parameters.

Sarabandi, K.↗

The Effects of Noise on Polarimetric SAR Data

Polarimetric SAR data can provide a great deal of information about the scattering behavior of the surface under observation. Polarimetric SAR systems often measure the scattering matrices of the areas under observation in linear polarizations (H and V). From the scattering matrix commonly used forms such as the covariance matrix and the Stokes matrix can be easily derived. Other measures derived from polarimetric SAR data include the standard deviation of texture, correlation coefficients between scattering matrix terms, and the mode and variance of phase differences between scattering matrix terms. The effects of additive system noise on these measurements is not often considered in the literature on this subject.

Polarimetric SAR Data analysis↗

The Effects of Noise on Polarimetric SAR Data

Polarimetric SAR data can provide a great deal of information about the scattering behavior of the surface under observation. Polarimetric SAR systems often measure the scattering matrices of the areas under observation in linear polarizations (H and V). From the scattering matrix commonly used forms such as the covariance matrix and the Stokes matrix can be easily derived. Other measures derived from polarimetric SAR data include correlation coefficients between scattering matrix terms and the mode and variance of phase differences between scattering matrix terms. The effects of additive system noise on these measurements is not often considered in the literature on this subject. In this paper, the effects of additive system noise on measurements derived from polarimetric SAR data will be examined. It will be shown how first-order noise effects can be removed and how second-order noise effects can be reduced for some measurements...

Freeman, A.↗

Retrieving Storm Electric Fields from Aircrfaft Field Mill Data: Part II: Applications

The Lagrange multiplier theory developed in Part I of this study is applied to complete a relative calibration of a Citation aircraft that is instrumented with six field mill sensors. When side constraints related to average fields are used, the Lagrange multiplier method performs well in computer simulations. For mill measurement errors of 1 V m(sup -1) and a 5 V m(sup -1) error in the mean fair-weather field function, the 3D storm electric field is retrieved to within an error of about 12%. A side constraint that involves estimating the detailed structure of the fair-weather field was also tested using computer simulations. For mill measurement errors of 1 V m(sup -l), the method retrieves the 3D storm field to within an error of about 8% if the fair-weather field estimate is typically within 1 V m(sup -1) of the true fair-weather field. Using this type of side constraint and data from fair-weather field maneuvers taken on 29 June 2001, the Citation aircraft was calibrated. Absolute calibration was completed using the pitch down method developed in Part I, and conventional analyses. The resulting calibration matrices were then used to retrieve storm electric fields during a Citation flight on 2 June 2001. The storm field results are encouraging and agree favorably in many respects with results derived from earlier (iterative) techniques of calibration.

Koshak, William↗

Generating accurate density matrices on the tangent space of a Grassmann manifold

Interpolating a density matrix from a set of known density matrices is not a trivial task. This is because a linear combination of density matrices does not necessarily correspond to another density matrix. In this Communication, density matrices are examined as objects of a Grassmann manifold. Although this manifold is not a vector space, its tangent space is a vector space. As a result, one can map the density matrices on this manifold to their corresponding vectors in the tangent space and then perform interpolations on that tangent space. The resulting interpolated vector can be mapped back to the Grassmann manifold, which can then be utilized (1) as an optimal initial guess for a self-consistent field (SCF) calculation or (2) to derive energy directly without time-consuming SCF iterations. Such a promising approach is denoted as Grassmann interpolation (G-Int). The hydrogen molecule has been used to illustrate that the described interpolated method in this work preserves the essential attributes of a density matrix. For phosphorus mononitride and ferrocene, it was demonstrated numerically that reference points for the definition of the corresponding tangent spaces can be chosen arbitrarily. In addition, the interpolated density matrices provide a superior and essentially converged initial guess for an SCF calculation to make the SCF procedure itself unnecessary. Finally, this accurate, efficient, robust, and systematically improved G-Int strategy has been used for the first time to generate highly accurate potential energy surfaces with fine details for the difficult case, ferrocene.

Chemistry↗

Carbamoylation as an Effective Tool in the Analysis of the Soman Nerve Agent Marker Pinacolyl Alcohol in Soil Matrices by EI-GC-MS and LC-HRMS

Pinacolyl alcohol (PA) is a Schedule 2 chemical commonly featured in most proficiency tests (PTs) administered by the Organisation for the Prohibition of Chemical Weapons (OPCW) due to its direct as well as forensic link to the nerve agent Soman. Therefore, its detection by Chemical Weapons Convention (CWC) inspection teams during on-site investigations is a strong indicator of the past or latent presence of Soman in the environment. Its small molecular weight (102), early elution time, and poor ionization profile make PA a challenging analyte to detect particularly at low concentrations (∼1–10 μg/g). In this work, 1,1′-carbonyldimidazole (CDI) has been used to effectively modify PA for the first time in two different soil matrices (Virginia type A soil and silt sediment) at two separate concentrations (1 and 10 μg/g) for its subsequent detection by EI-GC-MS and LC-HRMS methods. For the EI-GC-MS analysis, the PA carbamate derivative (PIC) exhibits improved chromatography relative to PA such as improved peak shape, increased molecular weight (196 for PIC and 102 for PA) and increased retention time (16.9 min for PIC and ∼4.1 min for PA). In addition, the derivatization also improves the detection of PA by LC-HRMS as the PIC product possesses protonation sites (i.e., imidazole ring) relative to none exhibited by PA. More importantly, the carbamoylation proceeds under mild conditions (55 °C, no base) and rapidly (3 h), characteristics that make it an appealing protocol for the analysis of PA during OPCW PTs or real case scenarios particularly in instances where it is present at low concentrations. It is anticipated that the protocol can be applied to the forensic analysis of this important Soman marker in various environmental matrices.

Chemistry↗