Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Matrix Models”

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,747 records · Page 97
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Theory of optical transitions in π-conjugated macrocycles

We describe a theoretical and computational investigation of the optical properties of π-conjugated macrocycles. Since the low-energy excitations of these systems are Frenkel excitons that couple to high-frequency dispersionless phonons, we employ the quantized Frenkel-Holstein model and solve it via the density matrix renormalization group (DMRG) method. First we consider optical emission from perfectly circular systems. Owing to optical selection rules, such systems radiate via two mechanisms: (i) within the Condon approximation, by thermally induced emission from the optically allowed j = ± 1 states and (ii) beyond the Condon approximation, by emission from the j = 0 state via coupling with a totally non-symmetric phonon (namely, the Herzberg-Teller effect). Using perturbation theory, we derive an expression for the Herzberg-Teller correction and show via DMRG calculations that this expression soon fails as ħ ω/J and the size of the macrocycle increase. Next, we consider the role of broken symmetry caused by torsional disorder. In this case the quantum number j no longer labels eigenstates of angular momentum, but instead labels localized local exciton groundstates (LEGSs) or quasi-extended states (QEESs). As for linear polymers, LEGSs define chromophores, with the higher energy QEESs being extended over numerous LEGSs. Within the Condon approximation (i.e., neglecting the Herzberg-Teller correction) we show that increased disorder increases the emissive optical intensity, because all the LEGSs are optically active. We next consider the combined role of broken symmetry and curvature, by explicitly evaluating the Herzberg-Teller correction in disordered systems via the DMRG method. The Herzberg-Teller correction is most evident in the emission intensity ratio, I{sub 00}/I{sub 01}. In the Condon approximation I{sub 00}/I{sub 01} is a constant function of curvature, whereas in practice it vanishes for closed rings and only approaches a constant in the limit of vanishing curvature. We calculate the optical spectra of a model system, cyclo-poly(para-phenylene ethynylene), for different amounts of torsional disorder within and beyond the Condon approximation. We show how broken symmetry and the Herzberg-Teller effect explain the spectral features. The Herzberg-Teller correction to the 0–1 emission vibronic peak is always significant. Finally, we note the qualitative similarities between the optical properties of conformationally disordered linear polymers and macrocycles in the limit of sufficiently large disorder, because in both cases they are determined by the optical properties of curved chromophores.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2018

Fano Resonance by Means of the One-Dimensional Superconductor Photonic Crystals

In the present work, we introduce a study about the reflectance properties of a one-dimensional superconductor photonic crystal in a different manner. Our structure is designed from alternately layers of high-temperature superconducting material (BSCCO) and a conventional superconducting material (Nb) terminated by the dielectric cap layer. The investigated numerical results are essentially based on the basics of the transfer matrix method and the two-fluid model. The numerical results showed the appearance of two sharp peaks, which are referred to Fano and electromagnetic-induced reflectance resonances. Also, we demonstrate the effects of many parameters on the reflectance properties of Fano resonance and the electromagnetic-induced reflectance. The effects of the constituent layer thicknesses, the number of periods, the operating temperature, the angle of incidence, and the refractive index of the dielectric cap layer are all considered. Finally, the effects of hydrostatic pressure on the Fano and EIR resonances are studied and calculated. This structure can play an important role in optical switching devices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2013

Optimal system size for complex dynamics in random neural networks near criticality

In this article, we consider a model of dynamical agents coupled through a random connectivity matrix, as introduced by Sompolinsky et al. [Phys. Rev. Lett. 61(3), 259–262 (1988)] in the context of random neural networks. When system size is infinite, it is known that increasing the disorder parameter induces a phase transition leading to chaotic dynamics. We observe and investigate here a novel phenomenon in the sub-critical regime for finite size systems: the probability of observing complex dynamics is maximal for an intermediate system size when the disorder is close enough to criticality. We give a more general explanation of this type of system size resonance in the framework of extreme values theory for eigenvalues of random matrices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Ferromagnetism and doublon localization in a Wannier-Hubbard chain

We derive a “Wannier-Hubbard” model consisting of an array of overlapping atomic orbitals interacting via a local Coulomb interaction. Transforming to an orthogonal Wannier basis set, the resulting Hamiltonian displays long-range hopping and interactions, with new terms such as correlated hopping and ferromagnetic direct exchange, among others. We numerically study the one-dimensional version of the model at half-filling using the density matrix renormalization group method, unveiling a rich phase diagram as a function of the interaction U and the overlap s with metallic and ferromagnetic phases separated by a ferrimagnetic region. Our results indicate a path toward understanding emergent phases under pressure and beyond standard model Hamiltonians.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Signatures of pairing in the magnetic excitation spectrum of strongly correlated two-leg ladders [Signatures of pairing in the magnetic excitation spectrum of strongly correlated ladders]

Magnetic interactions are widely believed to play a crucial role in the microscopic mechanism leading to high critical temperature superconductivity. It is therefore important to study the signatures of pairing in the magnetic excitation spectrum of simple models known to show unconventional superconducting tendencies. Using the density matrix renormalization group technique, we calculate the dynamical spin structure factor S(k,ω) of a generalized t–U–J Hubbard model away from half filling in a two-leg ladder geometry. The addition of J enhances pairing tendencies. We analyze quantitatively the signatures of pairing in the magnetic excitation spectra. We found that the superconducting pair-correlation strength, that can be estimated independently from ground state properties, is closely correlated with the integrated low-energy magnetic spectral weight in the vicinity of (π,π). In this wave-vector region, robust spin incommensurate features develop with increasing doping. The branch of the spectrum with rung direction wave vector k rung =0 does not change substantially with doping where pairing dominates and thus plays a minor role. As a result, we discuss the implications of our results for neutron scattering experiments, where the spin excitation dynamics of hole-doped quasi-one-dimensional magnetic materials can be measured and also address implications for recent resonant inelastic x-ray scattering experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1987

Partitioned Matrices for Combined Linear Systems

Multiple-input, multiple-output mathematical models combined. Report presents mathematical background for combination of linear multiple-input, multiple-output equations of subsystems into overall matrix equations for whole system. Development of such system mathematical models essential for design of complicated control systems.

Duke, Eugene L.↗
DOE OSTI2017

Band Gaps for Elastic Wave Propagation in a Periodic Composite Beam Structure Incorporating Microstructure and Surface Energy Effects

Here, a new model for determining band gaps for elastic wave propagation in a periodic composite beam structure is developed using a non-classical Bernoulli–Euler beam model that incorporates the microstructure, surface energy and rotational inertia effects. The Bloch theorem and transfer matrix method for periodic structures are employed in the formulation. The new model reduces to the classical elasticity-based model when both the microstructure and surface energy effects are not considered. The band gaps predicted by the new model depend on the microstructure and surface elasticity of each constituent material, the unit cell size, the rotational inertia, and the volume fraction. To quantitatively illustrate the effects of these factors, a parametric study is conducted. The numerical results reveal that the band gap predicted by the current non-classical model is always larger than that predicted by the classical model when the beam thickness is very small, but the difference is diminishing as the thickness becomes large. Also, it is found that the first frequency for producing the band gap and the band gap size decrease with the increase of the unit cell length according to both the current and classical models. In addition, it is observed that the effect of the rotational inertia is larger when the exciting frequency is higher and the unit cell length is smaller. Furthermore, it is seen that the volume fraction has a significant effect on the band gap size, and large band gaps can be obtained by tailoring the volume fraction and material parameters.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Triaxiality in selenium-76

Background: Selenium and germanium nuclei are associated with both triaxiality and shape coexistence. The relative influence of these deformation effects on the low-lying nuclear structure remains the subject of much discussion with additional attention drawn to 76 Se and 76 Ge due to the potential for the observation of neutrinoless double- β decay. Purpose: Experimental observables related to the deformation of 76 Se are lacking in precision. The purpose of the present paper is to provide electric quadrupole matrix elements with improved precision in order to determine the deformation of low-lying states in a model-independent manner. Methods: Sub-barrier Coulomb excitation was employed at the reaccelerated beam facility of the National Superconducting Cyclotron Laboratory using the JANUS setup. Using this method nineteen E2 matrix elements were extracted. Results: Extracted matrix elements agree within uncertainties with those in the literature but with improved precision. Through both a comparison with geometric models and a model-independent evaluation of E2 matrix elements using rotational invariants, the ground state of 76 Se is best described as having a significant triaxial component, but is not maximally triaxially deformed. Conclusions: Selenium-76 exhibits a significant degree of triaxiality in its ground state. Finally, a detailed comparison with configuration-interaction calculations indicates that this can be well reproduced theoretically.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2015

Incorporation of Plasticity and Damage Into an Orthotropic Three-Dimensional Model with Tabulated Input Suitable for Use in Composite Impact Problems

The need for accurate material models to simulate the deformation, damage and failure of polymer matrix composites under impact conditions is becoming critical as these materials are gaining increased usage in the aerospace and automotive industries. While there are several composite material models currently available within commercial transient dynamic finite element codes, several features have been identified as being lacking in the currently available material models that could substantially enhance the predictive capability of the impact simulations. A specific desired feature pertains to the incorporation of both plasticity and damage within the material model. Another desired feature relates to using experimentally based tabulated stress-strain input to define the evolution of plasticity and damage as opposed to specifying discrete input properties (such as modulus and strength) and employing analytical functions to track the response of the material. To begin to address these needs, a combined plasticity and damage model suitable for use with both solid and shell elements is being developed for implementation within the commercial code LS-DYNA. The plasticity model is based on extending the Tsai-Wu composite failure model into a strain-hardening based orthotropic plasticity model with a non-associative flow rule. The evolution of the yield surface is determined based on tabulated stress-strain curves in the various normal and shear directions and is tracked using the effective plastic strain. The effective plastic strain is computed by using the non-associative flow rule in combination with appropriate numerical methods. To compute the evolution of damage, a strain equivalent semi-coupled formulation is used, in which a load in one direction results in a stiffness reduction in multiple coordinate directions. A specific laminated composite is examined to demonstrate the process of characterizing and analyzing the response of a composite using the developed model.

Impact↗
NASA NTRS1985

A Kalman filter for a two-dimensional shallow-water model

A two-dimensional Kalman filter is described for data assimilation for making weather forecasts. The filter is regarded as superior to the optimal interpolation method because the filter determines the forecast error covariance matrix exactly instead of using an approximation. A generalized time step is defined which includes expressions for one time step of the forecast model, the error covariance matrix, the gain matrix, and the evolution of the covariance matrix. Subsequent time steps are achieved by quantifying the forecast variables or employing a linear extrapolation from a current variable set, assuming the forecast dynamics are linear. Calculations for the evolution of the error covariance matrix are banded, i.e., are performed only with the elements significantly different from zero. Experimental results are provided from an application of the filter to a shallow-water simulation covering a 6000 x 6000 km grid.

Parrish, D. F.↗
DOE OSTI2019

Analytical local stress model for UMo/Al dispersion fuel

The stress evolution occurring in UMo/Al dispersion fuel is important since it affects the fuel performance during irradiation. In this study, a new analytical model was developed to predict the local stresses in UMo/Al dispersion fuel. In the model, a hypothetical unit sphere composed of a UMo fuel particle, interaction layer (IL), and Al matrix was considered, and the governing equations for the stress-strain relationship, strain-displacement, and mechanical equilibrium were established using a spherical coordinate system. The mathematical derivations were obtained for local stresses in the radial and circumferential direction using a thick-walled sphere model. Furthermore, this analytical model employed the stress distribution as the boundary condition, which is calculated using a finite element model with homogenized fuel meat. The developed model’s solution scheme was verified against Abaqus solutions obtained for two irradiated plates with heterogeneous meat. The model calculated consistent results for interfacial stresses in the IL and Al matrix, indicating that the newly developed model was reliable when simultaneously calculating fission gas pressurization on the UMo/IL/Al composite, and the use of the stress distribution in the homogenized fuel meat as the boundary condition for the analytical model was acceptable

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2019

Evidence of pair-density wave in spin-valley locked systems

Cooper pairs with a finite center-of-mass momentum form a remarkable state in which the superconducting order parameter is modulated periodically in space. Although intense interest in such a “pair-density wave” (PDW) state has emerged due to recent discoveries in high $T_c$ superconductors, there is little theoretical understanding of the mechanism driving this exotic state. The challenge is that many competing states lie close in energy in seemingly simple models, such as the Hubbard model, in the strongly correlated regime. Here, we show that inversion symmetry breaking and the resulting spin-valley locking can promote PDWs over more commonly found spin stripes through frustration against magnetic order. Specifically, we find the first robust evidence for a PDW within density matrix renormalization group simulation of a simple fermionic model. Our results point to a tantalizing possibility in hole-doped group VI transition metal dichalcogenides, with spin-valley locked band structure and moderate correlations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Geomechanical and Fluid Transport Properties

Complex changes to fluid and mechanical properties of a pyrolyzing formation could lead to significant subsidence. Capturing subsidence requires accounting for changes in porosity due to generation and drainage of fluid. Pores, permeability and pore pressure are generated during pyrolysis of kerogen into oil and gas. In a series reaction scheme, oil is an intermediate product of pyrolysis. A large portion of oil-range fluids may be lost in an in-situ operation if the oil produced cannot be drained quickly enough. Petroleum fluid quality and production are highly dependent on the relative rates of creation, destruction, and transport of oil and gas in a given volume. All of these are dependent on the heating rate of oil shale. If oil shale is heated quickly, fluid flow may increase, but at a cost to the amount of oil that is finally extracted. To account for the effect of pore fluid on mechanical deformation, an elastic model was combined with a Mohr-Coulomb failure criterion and Biot’s linear consolidation model. The overall model is poroelastic, because it accounts for rock deforming by considering the response of the bulk system and the pore fluids it contains. Biot’s linear model includes all ways that pore fluid volume can change in a system, including stress-strain interactions between pores and the bulk material. Additionally, the model captures the distinction between fluid and solid matrix compression and describes how bulk deformation is affected by fluid drainage. The fluid response to changing stress is affected by changes in pore pressure, while the bulk material responds to effective stress changes. To capture all of these effects, a number of theoretical and empirical models were employed in FLAC3D™, a rock-mechanics software package produced by Itasca Consulting Group. The basic capabilities of FLAC3D™ have been supplemented by customized algorithms for the simulation package. The most important contribution has been to develop a methodology for describing poro-thermo-mechanical properties of oil shale across the complete chronology of in situ pyrolysis. This report highlights those methods and results for inferring realistic oil shale properties.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2012

Superallowed nuclear beta decay: Precision measurements for basic physics

For 60 years, superallowed 0{sup +}{yields}0{sup +} nuclear beta decay has been used to probe the weak interaction, currently verifying the conservation of the vector current (CVC) to high precision ({+-}0.01%) and anchoring the most demanding available test of the unitarity of the Cabibbo-Kobayashi-Maskawa (CKM) matrix ({+-}0.06%), a fundamental pillar of the electroweak standard model. Each superallowed transition is characterized by its ft-value, a result obtained from three measured quantities: the total decay energy of the transition, its branching ratio, and the half-life of the parent state. Today's data set is composed of some 150 independent measurements of 13 separate superallowed transitions covering a wide range of parent nuclei from {sup 10}C to {sup 74}Rb. Excellent consistency among the average results for all 13 transitions - a prediction of CVC - also confirms the validity of the small transition-dependent theoretical corrections that have been applied to account for isospin symmetry breaking. With CVC consistency established, the value of the vector coupling constant, G{sub V}, has been extracted from the data and used to determine the top left element of the CKM matrix, V{sub ud}. With this result the top-row unitarity test of the CKM matrix yields the value 0.99995(61), a result that sets a tight limit on possible new physics beyond the standard model. To have any impact on these fundamental weak-interaction tests, any measurement must be made with a precision of 0.1% or better - a substantial experimental challenge well beyond the requirements of most nuclear physics measurements. I overview the current state of the field and outline some of the requirements that need to be met by experimentalists if they aim to make measurements with this high level of precision.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2014

Theoretical Development of an Orthotropic Elasto-Plastic Generalized Composite Material Model

The need for accurate material models to simulate the deformation, damage and failure of polymer matrix composites is becoming critical as these materials are gaining increased usage in the aerospace and automotive industries. While there are several composite material models currently available within LSDYNA (Livermore Software Technology Corporation), there are several features that have been identified that could improve the predictive capability of a composite model. To address these needs, a combined plasticity and damage model suitable for use with both solid and shell elements is being developed and is being implemented into LS-DYNA as MAT_213. A key feature of the improved material model is the use of tabulated stress-strain data in a variety of coordinate directions to fully define the stress-strain response of the material. To date, the model development efforts have focused on creating the plasticity portion of the model. The Tsai-Wu composite failure model has been generalized and extended to a strain-hardening based orthotropic yield function with a nonassociative flow rule. The coefficients of the yield function, and the stresses to be used in both the yield function and the flow rule, are computed based on the input stress-strain curves using the effective plastic strain as the tracking variable. The coefficients in the flow rule are computed based on the obtained stress-strain data. The developed material model is suitable for implementation within LS-DYNA for use in analyzing the nonlinear response of polymer composites.

Polymer Matrix Composites↗