Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Matrix Theory”

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 1,477 records · Page 82

Ab Initio Uncertainty Quantification of Neutrinoless Double-Beta Decay in 76 Ge

The observation of neutrinoless double-beta (0⁢𝜈⁢𝛽⁢𝛽) decay would offer proof of lepton number violation, demonstrating that neutrinos are Majorana particles, while also helping us understand why there is more matter than antimatter in the Universe. If the decay is driven by the exchange of the three known light neutrinos, a discovery would, in addition, link the observed decay rate to the neutrino mass scale through a theoretical quantity known as the nuclear matrix element (NME). Accurate values of the NMEs for all nuclei considered for use in 0⁢𝜈⁢𝛽⁢𝛽 experiments are therefore crucial for designing and interpreting those experiments. Here, we report the first comprehensive ab initio uncertainty quantification of the 0⁢𝜈⁢𝛽⁢𝛽-decay NME, in the key nucleus 76 Ge. Here, our method employs nuclear strong and weak interactions derived within chiral effective field theory and recently developed many-body emulators. Our result, with a conservative treatment of uncertainty, is an NME of 2.60$^{+1.28}_{−1.36}$, which, together with the best-existing half-life sensitivity and phase-space factor, sets an upper limit for effective neutrino mass of 187$^{+205}_{−62}$ meV. The result is important for designing next generation germanium detectors aiming to cover the entire inverted hierarchy region of neutrino masses.

Ab initio calculations↗

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

An analytic solution for the orbital perturbations of the Venus Radar Mapper due to gravitational harmonics

Hill's variational equations are solved analytically for the orbital perturbations of a spacecraft nominally in an elliptic orbit around a non-spherical body. The rotation of the central planet about its spin-axis is not considered in the analysis. The perturbations are restricted to the planetary gravitational harmonics only. An extremely simple algorithm is derived to transform the spherical harmonic potentials to the orbital coordinate system, and the resulting accelerations are shown to be simply trigonometric functions of the true anomaly. With the principal matrix solution for the differential equations of the adjoint system given in closed form, the orthogonality of the trigonometric functions makes it possible to obtain an analytic solution for the non-homogeneous problem, at intervals of 2 pi in true anomaly. The solution for orbital perturbations can be extended over several revolutions by applying well-known results from Floquet's theory. The technique is demonstrated with results presented on the spacecraft periapsis altitude for the forthcoming Venus Radar Mapper Mission.

Vijayaraghavan, A.↗

Nonlinear Elastic Effects on the Energy Flux Deviation of Ultrasonic Waves in GR/EP Composites

In isotropic materials, the direction of the energy flux (energy per unit time per unit area) of an ultrasonic plane wave is always along the same direction as the normal to the wave front. In anisotropic materials, however, this is true only along symmetry directions. Along other directions, the energy flux of the wave deviates from the intended direction of propagation. This phenomenon is known as energy flux deviation and is illustrated. The direction of the energy flux is dependent on the elastic coefficients of the material. This effect has been demonstrated in many anisotropic crystalline materials. In transparent quartz crystals, Schlieren photographs have been obtained which allow visualization of the ultrasonic waves and the energy flux deviation. The energy flux deviation in graphite/epoxy (gr/ep) composite materials can be quite large because of their high anisotropy. The flux deviation angle has been calculated for unidirectional gr/ep composites as a function of both fiber orientation and fiber volume content. Experimental measurements have also been made in unidirectional composites. It has been further demonstrated that changes in composite materials which alter the elastic properties such as moisture absorption by the matrix or fiber degradation, can be detected nondestructively by measurements of the energy flux shift. In this research, the effects of nonlinear elasticity on energy flux deviation in unidirectional gr/ep composites were studied. Because of elastic nonlinearity, the angle of the energy flux deviation was shown to be a function of applied stress. This shift in flux deviation was modeled using acoustoelastic theory and the previously measured second and third order elastic stiffness coefficients for T300/5208 gr/ep. Two conditions of applied uniaxial stress were considered. In the first case, the direction of applied uniaxial stress was along the fiber axis (x3) while in the second case it was perpendicular to the fiber axis along the laminate stacking direction (x1).

Prosser, William H.↗

Application of the NASA Multiscale Analysis Tool: Multiscale Integration and Interoperability

The NASA Multiscale Analysis Tool (NASMAT) was developed recently to allow a wide variety of multiscale analysis problems to be effectively and efficiently solved. The architecture of NASMAT was established specifically to enable parallelized, “plug-and-play” functionality to reduce the complexity associated with adding new features to the code in the future and to allow end users to rapidly implement and evaluate user-defined capabilities. Additionally, the tool utilizes recursive data structures and subroutines to allow for an arbitrary number of length scales when performing multiscale analyses of heterogeneous materials. These features permit the rapid integration of user-defined capabilities (e.g., a material model, micromechanics approach, or failure theory) at all stages within a NASMAT calculation while leveraging built-in techniques where needed. Additionally, these features allow NASMAT to both be called from an external program as well as call an external program. This paper specifically focuses on the multiscale integration and interoperability of NASMAT with other analysis techniques through an illustrative, multiscale analysis of a 3D woven polymer matrix composite (PMC).

NASMAT↗

The adiabatic semiclassical perturbation theory for vibrationally inelastic scattering. I - Collinear calculations. II - Three-dimensional treatment

A semiclassical approximation to treat vibrationally inelastic scattering is developed. The vibrational basis set used is adiabatic with respect to a reference potential which is chosen to be as close as possible to the true potential and also gives easily obtainable solutions to the vibrational wave equation. The radial wave functions are obtained using the WKB approximation, and the coupled Schroedinger equations are solved by a first-order perturbation method to yield a phase shift matrix which is exponentiated to give the full scattering matrix. Results were obtained for all the cases computed by Secrest and Johnson and by Clark and Dickinson, and the agreement is better than 10% for half of the cross-sections and rarely off by more than a factor of 2.

Cross, R. J., Jr.↗

Assessing the theory-data tension in neutrino-induced charged pion production: the effect of final-state nucleon distortion

Pion production on nuclei constitutes a significant part of the total cross section in experiments involving few-GeV neutrinos. Combined analyses of data on deuterium and heavier nuclei points to tensions between the bubble chamber data and the data of the MINER$\nu$A experiment, which are often ascribed to unspecified nuclear effects. To understand the origin of these tensions, a microscopic quantum mechanical framework is needed to compute nuclear matrix elements. We use the local approximation to the relativistic distorted wave impulse approximation (RDWIA) to assess the role of final-state nucleon distortion. To perform this comparison under conditions relevant to neutrino experiments, we compute cross sections for the MINER$\nu$A and T2K charged pion production datasets. The inclusion of nucleon distortion leads to a reduction of the cross section up to 10%, but to no significant change in shape of the flux-averaged cross sections. Results with and without distortion compare favorably to experimental data, with the exception of the low-$Q^2$ MINER$\nu$A $\pi^+$ data. We point out that hydrogen target data from BEBC is also overpredicted at low-$Q^2$, and that the discrepancy is similar in shape and magnitude to what is found in comparison to MINER$\nu$A data. Including nucleon distortion alone cannot explain the overprediction of low-$Q^2$ cross sections measured by MINER$\nu$A. The similar overprediction of BEBC data on hydrogen means that it is impossible to ascribe this discrepancy solely to a nuclear effect. Axial couplings and their $Q^2$ dependence should ideally be derived from more precise data on hydrogen and deuterium.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

CrysFieldExplorer : rapid optimization of the crystal field Hamiltonian

A new approach to the fast optimization of crystal electric field (CEF) parameters to fit experimental data is presented. This approach is implemented in a lightweight Python-based program, CrysFieldExplorer. The main novelty of the method is the development of a unique loss function, referred to as the spectrum characteristic loss (L Spectrum ), which is based on the characteristic polynomial of the Hamiltonian matrix. Particle swarm optimization and a covariance matrix adaptation evolution strategy are used to find the minimum of the total loss function. It is demonstrated that CrysFieldExplorer can perform direct fitting of CEF parameters to any experimental data such as a neutron spectrum, susceptibility or magnetization measurements etc. CrysFieldExplorer can handle a large number of non-zero CEF parameters and reveal multiple local and global minimum solutions. Crystal field theory, the loss function, and the implementation and limitations of the program are discussed within the context of two examples.

42 ENGINEERING↗

Statistical analysis of the modal properties of large structural systems.

A theory is developed to predict eigenvalue and eigenvector uncertainty in large dynamic models. The uncertainty is based on physical property uncertainty and should not be confused with numerical roundoff, although the method can be extended to include the latter. The theory, when implemented on a computer, is used to analyze the uncertainties in frequencies and mode shapes based on uncertainties in mass, stiffness, modulus of elasticity, etc. The method incorporates a linear statistical model which is quite adequate for handling property uncertainties of 10% or more. The model is not limited to small systems but uses certain statistical assumptions as well as selective matrix manipulations to keep the size of all matrix operations to within the number of degrees of freedom of the system. Examples are given for two longitudinal vibration problems, and the results are supported by a Monte Carlo analysis.

Collins, J. D.↗

Closed-form solutions for a class of optimal quadratic regulator problems with terminal constraints

Closed-form solutions are derived for coupled Riccati-like matrix differential equations describing the solution of a class of optimal finite time quadratic regulator problems with terminal constraints. Analytical solutions are obtained for the feedback gains and the closed-loop response trajectory. A computational procedure is presented which introduces new variables for efficient computation of the terminal control law. Two examples are given to illustrate the validity and usefulness of the theory.

Juang, J.-N.↗

Radar polarimetry for geoscience applications

The present volume on radar polarimetry for geoscience applications discusses wave properties and polarization, scattering matrix representation for simple targets, scattering models for point and distributed targets, polarimetric scatterometer systems and measurements, polarimetric radar system design, and polarimetric SAR applications. Attention is given to plane waves in a lossless homogeneous medium-wave polarization, polarization synthesis and response, and coordinate system transformations. Topics addressed include high- and low-frequency scattering, rough-surface scattering models, radiative transfer theory and deficiencies thereof, solutions for the radiative transfer equation, and a radiative transfer model for a forest canopy. Also discussed are network analyzer-based polarimetric scatterometers, calibration of polarimetric scatterometers, synthesized polarization response of distributed targets, and measurement of the propagation parameters of a forest canopy.

Ulaby, Fawwaz T.↗

Effect of Thin-Film Adhesives on Mode I Interlaminar Fracture Toughness in Carbon Fiber Composites with Shape Memory Alloy Inserts

A single sheet of NiTi shape memory alloy (SMA) was introduced within a unidirectional HexPly® 8552/IM7 (Hexcel®) polymer matrix composite (PMC) panel in conjunction with multiple thin-film adhesives to promote the interfacial bond strength between the SMA and PMC. A double cantilever beam (DCB) test was performed in accordance with the ASTM D5528 method for evaluation of Mode I interlaminar fracture toughness of unidirectional fiber-reinforced PMCs. The modal acoustic emissions (MAEs) were monitored during testing with two acoustic sensors attached to the specimens. The composite panels were subjected to a C-scan before testing and examined using optical and scanning electron microscopy (SEM) techniques after part failure. The data were used in conjunction with modified beam theory (MBT), the compliance calibration (CC) method, and the modified compliance calibration (MCC) method. The Mode I interlaminar toughness (G(sub IC)) values for control specimens were higher than previously reported and are attributed to extensive fiber bridging during testing. The presence of adhesives with SMA inserts stabilized crack propagation during DCB testing. The results reveal a new phenomenon of SMA bridging, whereby crack propagation would switch from one side of the SMA insert to the other, thus increasing the load and G(sub IC) values of specimens containing SMA.

shape memory alloys↗

Improved Kelbg Potentials for Z > 1 and Application to Carbon Plasmas

In this work, we present a general form for the electron‐ion diffractive potential derived from the quantum pair density matrix and fit to the improved Kelbg potential for atomic numbers up to $Z = 54$. We apply classical molecular dynamics using the improved Kelbg potential for carbon with various forms of the Pauli potential to compute internal energies and pressures for hot, dense plasma conditions. Our results are compared to an equation of state model based on path integral Monte Carlo and density functional theory simulations to examine the extent to which the improved Kelbg potential reproduces the internal energy and pressure of carbon plasmas. The regions of validity for carbon agree generally with those derived previously for hydrogen once pressure ionization effects are incorporated. Based on our carbon results and previously published hydrogen studies, we discuss the general applicability and limitations of these potentials for equation of state studies in warm dense matter and high energy density plasmas.

general physics↗

Full-dimensional quantum stereodynamics of the non-adiabatic quenching of OH(A2Σ+) by H2

Abstract The Born–Oppenheimer approximation, assuming separable nuclear and electronic motion, is widely adopted for characterizing chemical reactions in a single electronic state. However, the breakdown of the Born–Oppenheimer approximation is omnipresent in chemistry, and a detailed understanding of the non-adiabatic dynamics is still incomplete. Here we investigate the non-adiabatic quenching of electronically excited OH(A 2 Σ + ) molecules by H 2 molecules using full-dimensional quantum dynamics calculations for zero total nuclear angular momentum using a high-quality diabatic-potential-energy matrix. Good agreement with experimental observations is found for the OH(X 2 Π) ro-vibrational distribution, and the non-adiabatic dynamics are shown to be controlled by stereodynamics, namely the relative orientation of the two reactants. The uncovering of a major (in)elastic channel, neglected in a previous analysis but confirmed by a recent experiment, resolves a long-standing experiment–theory disagreement concerning the branching ratio of the two electronic quenching channels.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A time-dependent momentum-resolved scattering approach to core-level spectroscopies

While new light sources allow for unprecedented resolution in experiments with X-rays, a theoretical understanding of the scattering cross-section is lacking. In the particular case of strongly correlated electron systems, numerical techniques are quite limited, since conventional approaches rely on calculating a response function (Kramers-Heisenberg formula) that is obtained from a perturbative analysis of scattering processes in the frequency domain. This requires a knowledge of a full set of eigenstates in order to account for all intermediate processes away from equilibrium, limiting the applicability to small tractable systems. In this work, we present an alternative paradigm, recasting the problem in the time domain and explicitly solving the time-dependent Schrödinger equation without the limitations of perturbation theory: a faithful simulation of the scattering processes taking place in actual experiments, including photons and core electrons. We show how this approach can yield the full time and momentum resolved Resonant Inelastic X-Ray Scattering (RIXS) spectrum of strongly interacting many-body systems. We demonstrate the formalism with an application to Mott insulating Hubbard chains using the time-dependent density matrix renormalization group method, which does not require a priory knowledge of the eigenstates and can solve very large systems with dozens of orbitals. This approach can readily be applied to systems out of equilibrium without modification and generalized to other spectroscopies.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Modeling Nucleation and Grain Growth in the Solar Nebula: Initial Progress Report

The primitive solar nebula was a violent and chaotic environment where high energy collisions, lightning, shocks and magnetic re-connection events rapidly vaporized some fraction of nebular dust, melted larger particles while leaving the largest grains virtually undisturbed. At the same time, some tiny grains containing very easily disturbed noble gas signatures (e.g., small, pre-solar graphite or SiC particles) never experienced this violence, yet can be found directly adjacent to much larger meteoritic components (chondrules or CAIs) that did. Additional components in the matrix of the most primitive carbonaceous chondrites and in some chondritic porous interplanetary dust particles include tiny nebular condensates, aggregates of condensates and partially annealed aggregates. Grains formed in violent transient events in the solar nebula did not come to equilibrium with their surroundings. To understand the formation and textures of these materials as well as their nebular abundances we must rely on Nucleation Theory and kinetic models of grain growth, coagulation and annealing. Such models have been very uncertain in the past: we will discuss the steps we are taking to increase their reliability.

Nuth, Joseph A.↗

Comparison of Multiscale Method of Cells-Based Models for Predicting Elastic Properties of Filament Wound C/C-SiC

Three different multiscale models, based on the method of cells (generalized and high fidelity) micromechanics models were developed and used to predict the elastic properties of C/C-SiC composites. In particular, the following multiscale modeling strategies were employed: Concurrent multiscale modeling of all phases using the generalized method of cells, synergistic (two-way coupling in space) multiscale modeling with the generalized method of cells, and hierarchical (one-way coupling in space) multiscale modeling with the high fidelity generalized method of cells. The three models are validated against data from a hierarchical multiscale finite element model in the literature for a repeating unit cell of C/C-SiC. Furthermore, the multiscale models are used in conjunction with classical lamination theory to predict the stiffness of C/C-SiC plates manufactured via a wet filament winding and liquid silicon infiltration process recently developed by the German Aerospace Institute.

multiscale modeling↗