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.
Stress and Damage in Polymer Matrix Composite Materials Due to Material Degradation at High Temperatures
This report describes analytical methods for calculating stresses and damage caused by degradation of the matrix constituent in polymer matrix composite materials. Laminate geometry, material properties, and matrix degradation states are specified as functions of position and time. Matrix shrinkage and property changes are modeled as functions of the degradation states. The model is incorporated into an existing composite mechanics computer code. Stresses, strains, and deformations at the laminate, ply, and micro levels are calculated, and from these calculations it is determined if there is failure of any kind. The rationale for the model (based on published experimental work) is presented, its integration into the laminate analysis code is outlined, and example results are given, with comparisons to existing material and structural data. The mechanisms behind the changes in properties and in surface cracking during long-term aging of polyimide matrix composites are clarified. High-temperature-material test methods are also evaluated.
Effects of Adiabatic Heating on the High Strain Rate Deformation of Polymer Matrix Composites
Polymer matrix composites (PMCs) are increasingly being used in aerospace structures that are expected to experience complex dynamic loading conditions throughout their lifetime. As such, a detailed understanding of the high strain rate behavior of the constituents, particularly the strain rate, temperature, and pressure dependent polymer matrix, is paramount. In this paper, preliminary efforts in modeling experimentally observed temperature rises due to plastic deformation in PMCs subjected to dynamic loading are presented. To this end, an existing isothermal viscoplastic polymer constitutive formulation is extended to model adiabatic conditions by incorporating temperature dependent elastic properties and modifying the components of the inelastic strain rate tensor to explicitly depend on temperature. It is demonstrated that the modified polymer constitutive model is capable of capturing strain rate and temperature dependent yield as well as thermal softening associated with the conversion of plastic work to heat at high rates of strain. The modified constitutive model is then embedded within a strength of materials based micromechanics framework to investigate the manifestation of matrix thermal softening, due to the conversion of plastic work to heat, on the high strain rate response of a T700Epon 862 (T700E862) unidirectional composite. Adiabatic model predictions for high strain rate composite longitudinal tensile, transverse tensile, and in-plane shear loading are presented. Results show a substantial deviation from isothermal conditions; significant thermal softening is observed for matrix dominated deformation modes (transverse tension and in-plane shear), highlighting the importance of accounting for the conversion of plastic work to heat in the polymer matrix in the high strain rate analysis of PMC structures.
Thermoplastic Ribbon-Ply Bonding Model
The aim of the present work was to identify key variables in rapid weldbonding of thermoplastic tow (ribbon) and their relationship to matrix polymer properties and to ribbon microstructure. Theoretical models for viscosity, establishment of ply-ply contact, instantaneous (Velcro) bonding, molecular interdiffusion (healing), void growth suppression, and gap filling were reviewed and synthesized. Consideration of the theoretical bonding mechanisms and length scales and of the experimental weld/peel data allow the prediction of such quantities as the time and pressure required to achieve good contact between a ribbon and a flat substrate, the time dependence of bond strength, pressures needed to prevent void growth from dissolved moisture and conditions for filling gaps and smoothing overlaps.
Stripe order enhanced superconductivity in the Hubbard model
Unidirectional (“stripe”) charge density wave order has now been established as a ubiquitous feature in the phase diagram of the cuprate high-temperature superconductors, where it generally competes with superconductivity. Nonetheless, on theoretical grounds it has been conjectured that stripe order (or other forms of “optimal” inhomogeneity) may play an essential positive role in the mechanism of high-temperature superconductivity. Here, we report density matrix renormalization group studies of the Hubbard model on long four- and six-leg cylinders, where the hopping matrix elements transverse to the long direction are periodically modulated—mimicking the effect of putative period 2 stripe order. We find that even modest amplitude modulations can enhance the long-distance superconducting correlations by many orders of magnitude and drive the system into a phase with a substantial spin gap and superconducting quasi–long-range order with a Luttinger exponent, K sc ~ 1.
Electronic specific heat capacities and entropies from density matrix quantum Monte Carlo using Gaussian process regression to find gradients of noisy data
In this work, we present a machine learning approach to calculating electronic specific heat capacities for a variety of benchmark molecular systems. Our models are based on data from density matrix quantum Monte Carlo, which is a stochastic method that can calculate the electronic energy at finite temperature. As these energies typically have noise, numerical derivatives of the energy can be challenging to find reliably. In order to circumvent this problem, we use Gaussian process regression to model the energy and use analytical derivatives to produce the specific heat capacity. From there, we also calculate the entropy by numerical integration. We compare our results to cubic splines and finite differences in a variety of molecules in which Hamiltonians can be diagonalized exactly with full configuration interaction. We finally apply this method to look at larger molecules where exact diagonalization is not possible and make comparisons with more approximate ways to calculate the specific heat capacity and entropy.
Two-component quantum Hall effects in topological flat bands
Here in this paper, we study quantum Hall states for two-component particles (hardcore bosons and fermions) loading in topological lattice models. By tuning the interplay of interspecies and intraspecies interactions, we demonstrate that two-component fractional quantum Hall states emerge at certain fractional filling factors ν = 1/2 for fermions (ν = 2/3 for bosons) in the lowest Chern band, classified by features from ground states including the unique Chern number matrix (inverse of the K matrix), the fractional charge and spin pumpings, and two parallel propagating edge modes. Moreover, we also apply our strategy to two-component fermions at integer filling factor ν = 2 , where a possible topological Neel antiferromagnetic phase is under intense debate very recently. For the typical π -flux checkerboard lattice, by tuning the onsite Hubbard repulsion, we establish a first-order phase transition directly from a two-component fermionic ν = 2 quantum Hall state at weak interaction to a topologically trivial antiferromagnetic insulator at strong interaction, and therefore exclude the possibility of an intermediate topological phase for our system.
Renormalization of the Standard Model Effective Field Theory from geometry
S-matrix elements are invariant under field redefinitions of the Lagrangian. They are determined by geometric quantities such as the curvature of the field-space manifold of scalar and gauge fields. We present a formalism where scalar and gauge fields are treated together, with a metric on the combined space of both types of fields. Scalar and gauge scattering amplitudes are given by the Riemann curvature Rijkl of this combined space, with indices i, j, k, l chosen to be scalar or gauge indices depending on the type of external particle. One-loop divergences can also be computed in terms of geometric invariants of the combined space, which greatly simplifies the computation of renormalization group equations. We apply our formalism to the Standard Model Effective Field Theory (SMEFT), and compute the renormalization group equations for even-parity bosonic operators to mass dimension eight.
A Global Modeling Study on Carbonaceous Aerosol Microphysical Characteristics and Radiative Effects
Recently, attention has been drawn towards black carbon aerosols as a short-term climate warming mitigation candidate. However the global and regional impacts of the direct, indirect and semi-direct aerosol effects are highly uncertain, due to the complex nature of aerosol evolution and the way that mixed, aged aerosols interact with clouds and radiation. A detailed aerosol microphysical scheme, MATRIX, embedded within the GISS climate model is used in this study to present a quantitative assessment of the impact of microphysical processes involving black carbon, such as emission size distributions and optical properties on aerosol cloud activation and radiative effects. Our best estimate for net direct and indirect aerosol radiative flux change between 1750 and 2000 is -0.56 W/m2. However, the direct and indirect aerosol effects are quite sensitive to the black and organic carbon size distribution and consequential mixing state. The net radiative flux change can vary between -0.32 to -0.75 W/m2 depending on these carbonaceous particle properties at emission. Taking into account internally mixed black carbon particles let us simulate correct aerosol absorption. Absorption of black carbon aerosols is amplified by sulfate and nitrate coatings and, even more strongly, by organic coatings. Black carbon mitigation scenarios generally showed reduced radiative fluxeswhen sources with a large proportion of black carbon, such as diesel, are reduced; however reducing sources with a larger organic carbon component as well, such as bio-fuels, does not necessarily lead to a reduction in positive radiative flux.
Phase diagram of the one-dimensional t 1 - t 2 - J model: Ferromagnetism, triplet pairing, and charge and pair density waves
Here, we present a density matrix renormalization group (DMRG) study of an extended t-J model with hopping to the first and second neighbors -- the one dimensional t 1 -t 2 -J model. The full phase diagram as a function of the density n and exchange strength J, for both positive and negative values of t 2 , is obtained. For t 2 = -0.5 we observe that, in the strongly interacting region, Nagaoka Ferromagnetism (FM) is accompanied by a triplet pair density wave (PDW) upon doping. As the spin exchange J increases, a charge density wave (CDW) phase emerges and then gives way to singlet superconductivity (SC). This phase behaves as a singlet PDW with vanishing spacial average of the order parameter and a spin gap. When t 2 = 0.5, the physics is basically reminiscent of the conventional t-J model, undergoing a transition from a metallic to a SC phase as a function of J.
Magnetic properties and pairing tendencies of the iron-based superconducting ladder BaFe 2 S 3 : Combined ab initio and density matrix renormalization group study
The recent discovery of superconductivity under high pressure in the two-leg ladder compound BaFe 2 S 3 [H. Takahashi et al., Nat. Mater. 14, 1008 (2015)] opens a broad avenue of research, because it represents the first report of pairing tendencies in a quasi-one-dimensional iron-based high-critical-temperature superconductor. Similarly, as in the case of the cuprates, ladders and chains can be far more accurately studied using many-body techniques and model Hamiltonians than their layered counterparts, particularly if several orbitals are active. In this publication, we derive a two-orbital Hubbard model from first principles that describes individual ladders of BaFe 2 S 3 . The model is studied with the density matrix renormalization group. These first reported results are exciting for two reasons: (i) at half-filling, ferromagnetic order emerges as the dominant magnetic pattern along the rungs of the ladder, and antiferromagnetic order along the legs, in excellent agreement with neutron experiments; and (ii) with hole doping, pairs form in the strong coupling regime, as found by studying the binding energy of two holes doped on the half-filled system. In addition, orbital selective Mott phase characteristics develop with doping, with only oneWannier orbital receiving the hole carriers while the other remains half-filled. Lastly, these results suggest that the analysis of models for iron-based two-leg ladders could clarify the origin of pairing tendencies and other exotic properties of iron-based high-critical-temperature superconductors in general.
Full-dimensional diabatic potential energy surfaces including dissociation: The {sup 2}E{sup ″} state of NO{sub 3}
A scheme to produce accurate full-dimensional coupled diabatic potential energy surfaces including dissociative regions and suitable for dynamical calculations is proposed. The scheme is successfully applied to model the two-sheeted surface of the {sup 2}E{sup ″} state of the NO{sub 3} radical. An accurate potential energy surface for the NO{sub 3}{sup −} anion ground state is developed as well. Both surfaces are based on high-level ab initio calculations. The model consists of a diabatic potential matrix, which is expanded to higher order in terms of symmetry polynomials of symmetry coordinates. The choice of coordinates is key for the accuracy of the obtained potential energy surfaces and is discussed in detail. A second central aspect is the generation of reference data to fit the expansion coefficients of the model for which a stochastic approach is proposed. A third ingredient is a new and simple scheme to handle problematic regions of the potential energy surfaces, resulting from the massive undersampling by the reference data unavoidable for high-dimensional problems. The final analytical diabatic surfaces are used to compute the lowest vibrational levels of NO{sub 3}{sup −} and the photo-electron detachment spectrum of NO{sub 3}{sup −} leading to the neutral radical in the {sup 2}E{sup ″} state by full dimensional multi-surface wave-packet propagation for NO{sub 3} performed using the Multi-Configuration Time Dependent Hartree method. The achieved agreement of the simulations with available experimental data demonstrates the power of the proposed scheme and the high quality of the obtained potential energy surfaces.
Coreflood Tests to Evaluate Enhanced Oil Recovery Potential of Wettability-Altering Surfactants for Oil-Wet Heterogeneous Carbonate Reservoirs
Summary Oil-wetness and heterogeneity are two main factors that result in low oil recovery (OR) by waterflood in carbonate reservoirs. The injected water is likely to flow through high-permeability regions and bypass the oil in the low-permeability matrix. In this study, systematic coreflood tests were carried out in both “homogeneous” cores and “heterogeneous” cores with a wettability-altering surfactant. The homogeneous coreflood tests were conducted to evaluate surfactant retention, as well as to compare tertiary surfactant flooding with secondary surfactant flooding. The heterogeneous coreflood test was proposed to model bypassing in low-permeability matrix during waterfloods, and dynamic imbibition of surfactant into the low-permeability matrix. Surfactant retention results suggest that retention increases as initial oil saturation decreases. The retention of selected surfactant in the target reservoir cores was measured to be within a range of 0.07–0.12 mg/g-rock, which is economically acceptable. The results of homogeneous coreflood tests showed that both secondary waterflood and secondary surfactant flood can achieve high OR (>50%) from relatively homogeneous oil-wet cores. A shut-in phase after the surfactant injection resulted in a surge in oil production, which suggests that enough time should be given for wettability alteration by surfactants. The results of heterogeneous coreflood tests showed that more oil is bypassed in the tighter matrix by waterflood if the permeability is higher in the flooded layer and this bypassed oil is the target for the wettability-altering surfactant floods. Slow wettability-altering surfactant injection leads to imbibition into bypassed regions. When the oil-wet carbonate reservoirs have large unswept regions after waterflood, wettability-altering surfactants can significantly improve OR if enough time is given for imbibition.
Inhibition of TGF-β pathway reverts extracellular matrix remodeling in T. cruzi-infected cardiac spheroids
Highlights: • A model of cardiac spheroids to study fibrosis in pathological processes is proposed. • The study evaluates a TGF-β inhibitor in restoring extracellular matrix balance disrupted by T. cruzi infection. • Inhibition of TGF-β signaling reduces fibrosis and hypertrophy in infected-cardiac spheroids. • The extracellular matrix remodeling is achieved through regulation of MMP-2 activity and TIMP-1 expression. Chagasic cardiomyopathy (CC) is the main manifestation of Chagas Disease (CD). CC is a progressive dysfunctional illness, in which transforming growth factor beta (TGF-β) plays a central role in fibrogenesis and hypertrophy. In the present study, we tested in a three-dimensional (3D) model of cardiac cells culture (named cardiac spheroids), capable of mimicking the aspects of fibrosis and hypertrophy observed in CC, the role of TGF-β pathway inhibition in restoring extracellular matrix (ECM) balance disrupted by T. cruzi infection. Treatment of T. cruzi-infected cardiac spheroids with SB 431542, a selective inhibitor of TGF-β type I receptor, resulted in a reduction in the size of spheroids, which was accompanied by a decrease in parasite load and in fibronectin expression. The inhibition of TGF-β pathway also promoted an increase in the activity of matrix metalloproteinase (MMP)-2 and a decrease in tissue inhibitor of matrix metalloproteinase (TIMP)-1 expression, which may be one of the mechanisms regulating extracellular matrix remodeling. Therefore, our study provides new insights into the molecular mechanisms by which inhibition of TGF-β signaling reverts fibrosis and hypertrophy generated by T. cruzi during CC and also highlights the use of cardiac spheroids as a valuable tool for the study of fibrogenesis and anti-fibrotic compounds.
The Influence of Shallow Subsurface Properties on Particle Motion in Acoustic-Seismic Coupling
Atmospheric acoustic waves transmit energy into the solid Earth through air-to-ground coupling. These waves are recorded by seismic sensors and provide insight into both atmospheric phenomena and subsurface properties. Interpreting these signals is often challenging because they are modulated by subsurface structure and the incidence angle of the acoustic wave. This study examines acoustic--seismic coupling generated by the 2012 Camp Minden Explosion, which was recorded by hundreds of seismoacoustic stations. We apply a novel technique to quantify the seismic particle motion, model coupled waves with a propagator matrix approach, and apply a Bayesian inversion to infer properties of the shallow subsurface. Our analysis reveals that prograde motion is widespread and focused in low shear-wave velocity regions, such as the Mississippi Embayment, and retrograde motion is more common in higher shear-velocity areas. Inversion results at some stations produce plausible subsurface models with strong waveform fits, while inversions at other sites are less successful. These results indicate prograde particle motion in air-to-ground coupled waves is more prevalent than previously recognized and may serve as a diagnostic for shallow velocity structure. Our comprehensive modeling and inversion framework provides a potential method to extract layered near-surface properties from acoustic-seismic coupling observations.
Production of heavy sterile neutrinos from vector boson decay at electroweak temperatures
Here, in the standard model extended with a seesaw mass matrix, we study the production of sterile neutrinos from the decay of vector bosons at temperatures near the masses of the electroweak bosons. We derive a general quantum kinetic equation for the production of sterile neutrinos and their effective mixing angles, which is applicable over a wide range of temperature, to all orders in interactions of the standard model and to leading order in a small mixing angle for the neutrinos. We emphasize the relation between the production rate and Landau damping at one-loop order and show that production rates and effective mixing angles depend sensitively upon the neutrino’s helicity. Sterile neutrinos with positive helicity interact more weakly with the medium than those with negative helicity, and their effective mixing angle is not modified significantly. Negative helicity states couple more strongly to the vector bosons, but their mixing angle is strongly suppressed by the medium. Consequently, if the mass of the sterile neutrino is ≲ 8.35 MeV , there are fewer states with negative helicity produced than those with positive helicity. There is an Mikheyev-Smirnov-Wolfenstein-type resonance in the absence of lepton asymmetry, but due to screening by the damping rate, the production rate is not enhanced. Sterile neutrinos with negative helicity freeze out at T$-\atop{f}$ ≃ 5 GeV , whereas positive helicity neutrinos freeze out at T$+\atop{f}$≃ 8 GeV , with both distributions far from thermal. As the temperature decreases, due to competition between a decreasing production rate and an increasing mixing angle, the distribution function for states with negative helicity is broader in momentum and hotter than that for those with positive helicity. Sterile neutrinos produced via vector boson decay do not satisfy the abundance, lifetime, and cosmological constraints to be the sole dark matter component in the Universe. Massive sterile neutrinos produced via vector boson decay might solve the 7 Li problem, albeit at the very edge of the possible parameter space. A heavy sterile neutrino with a mass of a few MeV could decay into light sterile neutrinos, of a few keV in mass, that contribute to warm dark matter. In conclusion, we argue that heavy sterile neutrinos with lifetime ≤1/H 0 reach local thermodynamic equilibrium.
Retrieval of Microphysical Characteristics of Particles in Atmospheres of Distant Comets from Ground-Based Polarimetry
We summarize unique aperture data on the degree of linear polarization observed for distant comets C/2010 S1, C/2010 R1, C/2011 KP36, C/2012 J1, C/2013 V4, and C/2014 A4 with heliocentric distances exceeding 3 AU. Observations have been carried out at the 6-m telescope of the Special Astrophysical Observatory of the Russian Academy of Sciences (Nizhnij Arkhyz, Russia) during the period from 2011 to 2016. The measured negative polarization proves to be significantly larger in absolute value than what is typically observed for comets close to the Sun. We compare the new observational data with the results of numerical modeling performed with the T-matrix and superposition T-matrix methods. In our computer simulations, we assume the cometary coma to be an optically thin cloud containing particles in the form of spheroids, fractal aggregates composed of spherical monomers, and mixtures of spheroids and aggregate particles. We obtain a good semi-quantitative agreement between all polarimetric data for the observed distant comets and the results of numerical modeling for the following models of the cometary dust: (i) a mixture of submicrometer water-ice oblate spheroids with aggregates composed of submicrometer silicate monomers; and (ii) a mixture of submicrometer water-ice oblate spheroids and aggregates consisting of both silicate and organic monomers. The microphysical parameters of these models are presented and discussed.
Diffusion analysis for two-phase metal-matrix composite
Diffusion controlled filament matrix interaction in a metal matrix composite, where the filaments and matrix comprise a two phase binary alloy system, was mathematically modeled. The problem of a diffusion controlled, two phase moving interface by means of a one dimensional, variable grid, finite difference technique was analyzed. Concentration dependent diffusion coefficients and equilibrium solubility limits were used, and the change in filament diameter and compositional changes in the matrix were calculated as a function of exposure time at elevated temperatures. With the tungsten nickel (W-Ni) system as a model composite system, unidirectional composites containing from 0.06 to 0.44 initial filament volume fraction were modeled. Compositional changes in the matrix were calculated by superposition of the contributions from neighboring filaments. Alternate methods for determining compositional changes between first and second nearest neighbor filaments were also considered. The results show the relative importance of filament volume fraction, filament diameter, exposure temperature, and exposure time as they affect the rate and extent of filament matrix interaction.