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,045 records · Page 58

Differential continuum damage mechanics models for creep and fatigue of unidirectional metal matrix composites

Three multiaxial isothermal continuum damage mechanics models for creep, fatigue, and creep/fatigue interaction of a unidirectional metal matrix composite volume element are presented, only one of which will be discussed in depth. Each model is phenomenological and stress based, with varying degrees of complexity to accurately predict the initiation and propagation of intergranular and transgranular defects over a wide range of loading conditions. The development of these models is founded on the definition of an initially transversely isotropic fatigue limit surface, static fracture surface, normalized stress amplitude function and isochronous creep damage failure surface, from which both fatigue and creep damage evolutionary laws can be obtained. The anisotropy of each model is defined through physically meaningful invariants reflecting the local stress and material orientation. All three transversely isotropic models have been shown, when taken to their isotropic limit, to directly simplify to previously developed and validated creep and fatigue continuum damage theories. Results of a nondimensional parametric study illustrate (1) the flexibility of the present formulation when attempting to characterize a large class of composite materials, and (2) its ability to predict anticipated qualitative trends in the fatigue behavior of unidirectional metal matrix composites. Additionally, the potential for the inclusion of various micromechanical effects (e.g., fiber/matrix bond strength, fiber volume fraction, etc.), into the phenomenological anisotropic parameters is noted, as well as a detailed discussion regarding the necessary exploratory and characterization experiments needed to utilize the featured damage theories.

Arnold, S. M.↗

Continuous phase transition between bosonic integer quantum Hall liquid and a trivial insulator: Evidence for deconfined quantum criticality

The deconfined quantum critical point, a prototype Landau-forbidden transition, could exist, in principle, in the phase transitions involving a symmetry-protected topological phase; however, examples of such kinds of transitions in physical systems are rare beyond one-dimensional systems. Here, using a density-matrix renormalization-group calculation, we unveil a bosonic integer quantum Hall phase in a two-dimensional correlated honeycomb lattice, by full identification of its internal structure from the topological K matrix. Furthermore, we demonstrate that imbalanced periodic chemical potentials can destroy the bosonic integer quantum Hall state and drive it into a featureless trivial (Mott) insulator, where all physical observables evolve smoothly across the critical point. At the critical point, the entanglement entropy reveals a characteristic scaling behavior, which is consistent with the critical field theory as an emergent QED3 with two flavors of Dirac fermions.

36 MATERIALS SCIENCE↗

Towards a nonperturbative construction of the S-matrix

We present a nonperturbative recipe for directly computing the S-matrix in strongly-coupled QFTs. The method makes use of spectral data obtained in a Hamiltonian framework and can be applied to a wide range of theories, including potentially QCD. We demonstrate the utility of this prescription in the specific example of the 2+1d O(N) model at large N, using energy eigenstates computed with Hamiltonian truncation to reproduce the full 2 → 2 scattering amplitude for arbitrary (complex) center-of-mass energy.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Electronic energy gap closure and metal-insulator transition in dense liquid hydrogen

Here, using quantum Monte Carlo (QMC) calculations, we investigate the insulator-metal transition observed in liquid hydrogen at high pressure. Below the critical temperature of the transition from the molecular to the atomic liquid, the fundamental electronic gap closure occurs abruptly, with a small discontinuity reflecting the weak first-order transition in the thermodynamic equation of state. Above the critical temperature, molecular dissociation sets in while the gap is still open. When the gap closes, the decay of the off-diagonal reduced density matrix shows that the liquid enters a gapless, but localized, phase: there is a crossover between the insulating and the metallic liquids. Compared to different density functional theory (DFT) functionals, our QMC calculations provide larger values for the fundamental gap and the electronic density of states close to the band edges, indicating that optical properties from DFT potentially benefit from error cancellations.

36 MATERIALS SCIENCE↗

Geometry program for aerodynamic lifting surface theory

A computer program that provides the geometry and boundary conditions appropriate for an analysis of a lifting, thin wing with control surfaces in linearized, subsonic, steady flow is presented. The kernel function method lifting surface theory is applied. The data which is generated by the program is stored on disk files or tapes for later use by programs which calculate an influence matrix, plot the wing planform, and evaluate the loads on the wing. In addition to processing data for subsequent use in a lifting surface analysis, the program is useful for computing area and mean geometric chords of the wing and control surfaces.

Medan, R. T.↗

Thermo-viscoelastic analysis of composite materials, volume 1

Advanced composite materials, especially graphite/epoxy, are being applied to aircraft structures in order to improve performance and save weight. An important consideration in composite design is the residual strength of a structure containing holes, delaminations, or interlaminar damage when subjected to compressive loads. Recent studies have revealed the importance of viscoelastic effects in polymer-based composites. The viscoelastic effect is particularly significant at elevated temperature/moisture conditions since the matrix material is strongly affected by the environment. The solution of viscoelastic problems in composites was limited to special cases which can be solved by classical lamination theory. A finite element procedure is presented for calculating time-dependent stresses and strains in composite structures with general configurations and complicated boundary conditions. Using this procedure the in-plane and interlaminar stress distributions and histories in notched and unnotched composites were obtained for mechanical and thermal loads. Both two-dimensional and three-dimensional viscoelastic problems are analyzed. The effects of layup orientation and load spectrum on creep response and stress relaxation were also studied.

Lin, K. Y.↗

Numerically stable algorithm for discrete-ordinate-method radiative transfer in multiple scattering and emitting layered media

The transfer of monochromatic radiation in a scattering, absorbing, and emitting plane-parallel medium with a specified bidirectional reflectivity at the lower boundary is considered. The equations and boundary conditions are summarized. The numerical implementation of the theory is discussed with attention given to the reliable and efficient computation of eigenvalues and eigenvectors. Ways of avoiding fatal overflows and ill-conditioning in the matrix inversion needed to determine the integration constants are also presented.

Stamnes, Knut↗

A line-of-sight performance criterion for controller design of a proposed laboratory model

A line-of-sight performance criterion is derived for a proposed Controls Structures Interaction model, and its many uses in the control design process for fine pointing control are illustrated. A linearized line-of-sight (LOS) criterion is used for direct controller design and as a performance measure to judge different control methodologies. Numerical simulation results are shown where the three approaches: linear quadratic Gaussian theory, robust eigensystem assignment, and local velocity feedback are used for vibration control. Results indicate that the linear quadratic Gaussian controller, which incorporates a linearized LOS weighting matrix directly, yields good performance without wasting energy to control motions that have no influence on the LOS.

Lim, Kyong B.↗

Scattering from dielectric structures above impedance surfaces and resistive sheets

Interest in understanding of electromagnetic interaction with rough surfaces has prompted the study of scattering from typical dielectric humps over impedance surfaces. It is shown that the Green's function of the problem for a resistive sheet resembles that of the impedance surface. Hence both problems are considered here. A numerical solution for the scattered field of a two-dimensional dielectric object, possibly inhomogeneous, with arbitrary cross section above the impedance surface or resistive sheet is sought. First the Green's function of the problem is derived based on the exact image theory. This form of the Green's function is amenable to numerical computation. Then the induced polarization currents are calculated by casting the integral equations into a matrix equation via the method of moments. Numerical problems in calculation of the Green's function when both source and observation points are close to the surface are discussed. Comparison of numerical results with a perturbation solution shows excellent agreement between the two methods.

Sarabandi, Kamal↗

Quarkonium semiclassical transport in quark-gluon plasma: factorization and quantum correction

We study quarkonium transport in the quark-gluon plasma by using the potential nonrelativistic QCD (pNRQCD) effective field theory and the framework of open quantum systems. We argue that the coupling between quarkonium and the thermal bath is weak using separation of scales, so the initial density matrix of the total system factorizes and the time evolution of the subsystem is Markovian. We derive the semiclassical Boltzmann equation for quarkonium by applying a Wigner transform to the Lindblad equation and carrying out a semiclassical expansion. We resum relevant interactions to all orders in the coupling constant at leading power of the nonrelativistic and multipole expansions. The derivation is valid for both weakly coupled and strongly coupled quark-gluon plasmas. We find reaction rates in the transport equation factorize into a quarkonium dipole transition function and a chromoelectric gluon distribution function. For the differential reaction rate, the definition of the momentum dependent chromoelectric gluon distribution function involves staple-shaped Wilson lines. For the inclusive reaction rate, the Wilson lines collapse into a straight line along the real time axis and the distribution becomes momentum independent. The relation between the two Wilson lines is analogous to the relation between the Wilson lines appearing in the gluon parton distribution function (PDF) and the gluon transverse momentum dependent parton distribution function (TMDPDF). The centrality dependence of the quarkonium nuclear modification factor measured by experiments probes the momentum independent distribution while the transverse momentum dependence and measurements of the azimuthal angular anisotropy may be able to probe the momentum dependent one. We discuss one way to indirectly constrain the quarkonium in-medium real potential by using the factorization formula and lattice calculations. The leading quantum correction to the semiclassical transport equation of quarkonium is also worked out. The study can be easily generalized to quarkonium transport in cold nuclear matter, which is relevant for quarkonium production in eA collisions in the future Electron-Ion Collider.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Stochastic microstructural analysis of ceramic matrix composites using a high-fidelity multiscale framework

This paper presents a high-fidelity multiscale framework to simulate the mechanical behavior of ceramic matrix composites (CMCs), accounting for the complex micromorphology captured using detailed material characterization. First, high-resolution micrographs are obtained for the specific carbon fiber silicon-carbide-nitride matrix (C/SiNC) CMC to characterize the variability of the architectural features and manufacturing-induced defects at the microscale. An image processing algorithm is then used to precisely estimate the size and distribution of all subscale features and defects from the micrographs. The information is then used to generate a three-dimensional stochastic representative volume element (SRVE) to reconstruct microscale constituents accounting for the variability. As a result, the generated SRVEs are simulated using the high-fidelity generalized method of cells (HFGMC) micromechanics theory to investigate the effects of defects on the elastic properties of C/SiNC CMCs.

Ceramic Matrix Composites↗

Quantum electron transport in degenerate donor–acceptor systems

In this work, we develop a mathematically rigorous theory for the quantum transfer processes in degenerate donor–acceptor dimers in contact with a thermal environment. We explicitly calculate the transfer rates and the acceptor population efficiency. The latter depends critically on the initial donor state. Additionally, we show that quantum coherence in the initial state enhances the transfer process. If the electron is initially shared coherently by the donor levels, then the efficiency can reach values close to 100%, while an incoherent initial donor state will significantly suppress the efficiency. The results are useful for a better understanding of the quantum electron transport in many chemical, solid state, and biological systems with complex degenerate and quasi-degenerate energy landscapes.

36 MATERIALS SCIENCE↗

2 νββ spectrum in chiral effective field theory

We investigate two-neutrino double beta decay (2νββ) in chiral effective field theory. We find contributions from weak magnetism and double-weak pion-exchange at next-to-leading-order in the chiral power counting. We discuss the impact of the chiral corrections on the electron spectra and find that they should be included in analyses of 2νββ decay that aim to uncover new physics signatures in the electron spectrum. We illustrate this point by revisiting the effect of sterile neutrinos and non-standard charged interactions. We also find that the pion-exchange contributions involve nuclear matrix elements that are related to those appearing in neutrinoless double beta decay (0νββ). We investigate whether the 0νββ nuclear matrix elements can be obtained from detailed measurements of the energy spectrum of the outgoing electrons in 2νββ transitions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Performance of an atomic mean-field spin–orbit approach within exact two-component theory for perturbative treatment of spin–orbit coupling

The paper reports development of an atomic mean-field (AMF) spin–orbit approach within the spin-free exact two-component theory in its one-electron variant (SFX2C-1e), a new relativistic quantum-chemical method for treating spin-orbit coupling in molecules. The effective one-electron spin–orbit integral matrix in the four-component representation is assembled as a direct sum of one-centre spin–orbit integral matrices with the mean-field two-electron contributions evaluated using atomic SFX2C-1e Hartree–Fock density matrices. It is then transformed into two-component representation using analytic SFX2C-1e energy derivative formulation. The resulting two-component spin–orbit integral matrix is by design suitable for use in perturbative calculations of spin–orbit coupling, treating SFX2C-1e wavefunctions as unperturbed states. The accuracy of the present AMF approach has been demonstrated using benchmark calculations of spin–orbit splittings for representative diatomic radicals at the equation-of-motion coupled-cluster singles and doubles level. In conclusion, to demonstrate the applicability and accuracy of the present perturbative spin–orbit scheme in calculations of challenging heavy-element containing systems, a thorough computational investigation of six low-lying electronic states of ThO+ is reported.

74 ATOMIC AND MOLECULAR PHYSICS↗

Calculation of ion–ion mutual neutralization rate constants using Landau–Zener theory coupled with trajectory simulations for Ar + –Cl − , Br − , I −

In this computational study, we self-consistently calculate the rate constants of mutual neutralization reactions by incorporating the electron transfer probability, using Landau–Zener state transition theory with inputs derived from ab initio quantum chemistry calculations, into classical trajectory simulations. Electronic structure calculations are done using correlation consistent basis sets with multi-reference configuration interaction to map all the molecular electronic states below the ion-dissociation limit as a function of the distance between the reacting species. Our electronic structure calculations have been significantly improved from our previous work through improved selection of molecular electronic configurations maintaining a fine grid of 1a 0 over a wide range of bond lengths and accurate treatment of spin–orbit couplings. Non-adiabatic coupling matrix elements are calculated with the three-point central difference method near each avoided crossing to estimate the exact crossing point R x and coupling parameter H if , which are inputs to the multi-channel Landau–Zener theory to calculate the electron transition probability. Our approach is applied to estimate the mutual neutralization rate constants for the following ion pairs: Ar + –Cl − , Ar + –Br − , Ar + –I − at ∼133 Pa. Furthermore, our predictions are compared against the experimental data reported. It is seen that the improvement in the electronic structure calculation results in excellent agreement between the simulation results and the available experimental data to within a factor of ∼2 or ∼±50%.

Complete-active space self-consistent field↗

Validation of three-dimensional incompressible spatial direct numerical simulation code: A comparison with linear stability and parabolic stability equation theories for boundary-layer transition on a flat plate

Spatially evolving instabilities in a boundary layer on a flat plate are computed by direct numerical simulation (DNS) of the incompressible Navier-Stokes equations. In a truncated physical domain, a nonstaggered mesh is used for the grid. A Chebyshev-collocation method is used normal to the wall; finite difference and compact difference methods are used in the streamwise direction; and a Fourier series is used in the spanwise direction. For time stepping, implicit Crank-Nicolson and explicit Runge-Kutta schemes are used to the time-splitting method. The influence-matrix technique is used to solve the pressure equation. At the outflow boundary, the buffer-domain technique is used to prevent convective wave reflection or upstream propagation of information from the boundary. Results of the DNS are compared with those from both linear stability theory (LST) and parabolized stability equation (PSE) theory. Computed disturbance amplitudes and phases are in very good agreement with those of LST (for small inflow disturbance amplitudes). A measure of the sensitivity of the inflow condition is demonstrated with both LST and PSE theory used to approximate inflows. Although the DNS numerics are very different than those of PSE theory, the results are in good agreement. A small discrepancy in the results that does occur is likely a result of the variation in PSE boundary condition treatment in the far field. Finally, a small-amplitude wave triad is forced at the inflow, and simulation results are compared with those of LST. Again, very good agreement is found between DNS and LST results for the 3-D simulations, the implication being that the disturbance amplitudes are sufficiently small that nonlinear interactions are negligible.

Joslin, Ronald D.↗

Coherent spin states and emergent de Sitter quasinormal modes

As a toy model for the microscopic description of matter in de Sitter space, we consider a Hamiltonian acting on the spin-j representation of SU(2). This is a model with a finite-dimensional Hilbert space, from which quasinormal modes emerge in the large-spin limit. The path integral over coherent spin states can be evaluated at the semiclassical level and from it we find the single-particle de Sitter density of states, including 1/j corrections. Along the way, we discuss the use of quasinormal modes in quantum mechanics, starting from the paradigmatic upside-down harmonic oscillator.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A General Materials Data Science Framework for Quantitative 2D Analysis of Particle Growth from Image Sequences

Abstract Phase transformations are a challenging problem in materials science, which lead to changes in properties and may impact performance of material systems in various applications. We introduce a general framework for the analysis of particle growth kinetics by utilizing concepts from machine learning and graph theory. As a model system, we use image sequences of atomic force microscopy showing the crystallization of an amorphous fluoroelastomer film. To identify crystalline particles in an amorphous matrix and track the temporal evolution of the particle dispersion, we have developed quantitative methods of 2D analysis. 700 image sequences were analyzed using a neural network architecture, achieving 0.97 pixel-wise classification accuracy as a measure of the correctly classified pixels. The growth kinetics of isolated and impinged particles were tracked throughout time using these image sequences. The relationship between image sequences and spatiotemporal graph representations was explored to identify the proximity of crystallites from each other. The framework enables the analysis of all image sequences without the requirement of sampling for specific particles or timesteps for various materials systems.

36 MATERIALS SCIENCE↗