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,135 records · Page 63
All sourcesNASA NTRSDOE OSTI
NASA NTRS2004

Associative Flow Rule Used to Include Hydrostatic Stress Effects in Analysis of Strain-Rate-Dependent Deformation of Polymer Matrix Composites

designing reliable composite engine cases that are lighter than the metal cases in current use. The types of polymer matrix composites that are likely to be used in such an application have a deformation response that is nonlinear and that varies with strain rate. The nonlinearity and the strain-rate dependence of the composite response are due primarily to the matrix constituent. Therefore, in developing material models to be used in the design of impact-resistant composite engine cases, the deformation of the polymer matrix must be correctly analyzed. However, unlike in metals, the nonlinear response of polymers depends on the hydrostatic stresses, which must be accounted for within an analytical model. By applying micromechanics techniques along with given fiber properties, one can also determine the effects of the hydrostatic stresses in the polymer on the overall composite deformation response. First efforts to account for the hydrostatic stress effects in the composite deformation applied purely empirical methods that relied on composite-level data. In later efforts, to allow polymer properties to be characterized solely on the basis of polymer data, researchers at the NASA Glenn Research Center developed equations to model the polymers that were based on a non-associative flow rule, and efforts to use these equations to simulate the deformation of representative polymer materials were reasonably successful. However, these equations were found to have difficulty in correctly analyzing the multiaxial stress states found in the polymer matrix constituent of a composite material. To correct these difficulties, and to allow for the accurate simulation of the nonlinear strain-rate-dependent deformation analysis of polymer matrix composites, in the efforts reported here Glenn researchers reformulated the polymer constitutive equations from basic principles using the concept of an associative flow rule. These revised equations were characterized and validated in an experimental program carried out through a university grant with the Ohio State University, wherein tensile and shear deformation data were obtained for a representative polymer for strain rates ranging from quasi-static to high rates of several hundred per second. Tensile deformation data also were obtained over a variety of strain rates and fiber orientation angles for a representative polymer matrix composite composed using the polymer.

Goldberg, Robert K.↗
NASA NTRS2000

A Study of Influencing Factors on the Tensile Response of a Titanium Matrix Composite With Weak Interfacial Bonding

The generalized method of cells micromechanics model is utilized to analyze the tensile stress-strain response of a representative titanium matrix composite with weak interfacial bonding. The fiber/matrix interface is modeled through application of a displacement discontinuity between the fiber and matrix once a critical debonding stress has been exceeded. Unidirectional composites with loading parallel and perpendicular to the fibers are examined, as well as a cross-ply laminate. For each of the laminates studied, analytically obtained results are compared to experimental data. The application of residual stresses through a cool-down process was found to have a significant effect on the tensile response. For the unidirectional laminate with loading applied perpendicular to the fibers, fiber packing and fiber shape were shown to have a significant effect on the predicted tensile response. Furthermore, the interface was characterized through the use of semi-emperical parameters including an interfacial compliance and a "debond stress;" defined as the stress level across the interface which activates fiber/matrix debonding. The results in this paper demonstrate that if architectural factors are correctly accounted for and the interface is appropriately characterized, the macro-level composite behavior can be correctly predicted without modifying any of the fiber or matrix constituent properties.

Goldberg, Robert K.↗
NASA NTRS1992

Spatial orientation of the vestibular system

1. A simplified three-dimensional state space model of visual vestibular interaction was formulated. Matrix and dynamical system operators representing coupling from the semicircular canals and the visual system to the velocity storage integrator were incorporated into the model. 2. It was postulated that the system matrix for a tilted position was a composition of two linear transformations of the system matrix for the upright position. One transformation modifies the eigenvalues of the system matrix while another rotates the pitch and roll eigenvectors with the head, while maintaining the yaw axis eigenvector approximately spatially invariant. Using this representation, the response characteristics of the pitch, roll, and yaw eye velocity were obtained in terms of the eigenvalues and associated eigenvectors. 3. Using OKAN data obtained from monkeys and comparing to the model predictions, the eigenvalues and eigenvectors of the system matrix were identified as a function of tilt to the side or of tilt to the prone positions, using a modification of the Marquardt algorithm. The yaw eigenvector for right-side-down tilt and for downward pitch cross-coupling was approximately 30 degrees from the spatial vertical. For the prone position, the eigenvector was computed to be approximately 20 degrees relative to the spatial vertical. For both side-down and prone positions, oblique OKN induced along eigenvector directions generated OKAN which decayed to zero along a straight line with approximately a single time constant. This was verified by a spectral analysis of the residual sequence about the straight line fit to the decaying data. The residual sequence was associated with a narrow autocorrelation function and a wide power spectrum. 4. Parameters found using the Marquardt algorithm were incorporated into the model. Diagonal matrices in a head coordinate frame were introduced to represent the direct pathway and the coupling of the visual system to the integrator. Model simulations predicted the behavior of yaw and pitch OKN and OKAN when the animal was upright, as well as the cross-coupling in the tilted position. The trajectories in velocity space were also accurately simulated. 5. There were similarities between the monkey eigenvectors and human perception of the spatial vertical. For side-down tilts and downward eye velocity cross-coupling, there was only an Aubert (A) effect. For upward eye velocity cross-coupling there were both Muller (E) and Aubert (A) effects. The mean of the eigenvectors for upward and downward eye velocities overlay human 1 x g perceptual data.(ABSTRACT TRUNCATED AT 400 WORDS).

NASA Discipline Neuroscience↗
DOE OSTI2017

Architecture Independent Integrated Early Performance and Energy Estimation

The end of Dennard’s scaling, coupled to the necessity to keep increasing computational performance in constrained power envelopes, is leading developers towards the need to opti- mize applications not only for performance, but also for energy consumption. In such a scenario, solutions that enable estimating performance and energy consumption of a program as early as possible, even during development and without the necessity to access the final target hardware, can become fundamental tools to reduce the design space, the optimization effort, and provide potential opportunities to automate the optimization process itself. In this position paper, we discuss the design of a performance and energy estimation framework based on a retargetable com- piler. Our proposed approach targets the Intermediate Repre- sentation (IR) of the LLVM compiler. The IR based approach enables, after a model is built, estimation from a host different from the target architecture, even considering dynamic infor- mation. We present the rationale behind such a framework, identify opportunities and components readily available (such as the provided static analysis, quick ways to add instrumentation for dynamic profiling, and the multiple backends), and aspects that still need further research efforts (such as more effective non-linear, or machine learning based models). We also discuss a hand-developed case study, based on a simple sparse matrix vector multiplication with a linear model, to motivate the needs for such a framework.

Estimation, Compilers, Performance Analysis↗
NASA NTRS1984

Wave normal directions of chorus near the equatorial source region

OGO 5 triaxial search coil magnetometer data are used to determine the wave propagation directions of postmidnight chorus in the near-equatorial region at L shells of 6 to 7. The methods used to estimate the wave normal directions involved minimum variance, the imaginary part of the cross-spectral matrix, the eigenvector of the Hermitian cross-spectral matrix, and a fitting of dispersion relations for one-wave and two-wave models to the cross-spectral matrix. Wave propagation at all frequencies within chorus tones is found to occur most frequently along the magnetic field with median and average cone angles of 9.1 deg and 12.2 deg, respectively. It is concluded that for waves propagating parallel to B, wave growth is maximum.

Goldstein, B. E.↗
DOE OSTI2021

Fragility of the nematic spin liquid induced by diagonal couplings in the square-lattice SU(3) model

Here, we present a large-scale density matrix renormalization group (DMRG) study of the spin-1 SU(3) bilinear-biquadratic model on the square lattice, which was suggested to host a nematic spin liquid state in recent DMRG calculations. We report that this spin liquid appears to strongly compete with a three-sublattice magnetic order. To further study the competition between the two states, and the reason of the emergent nematic spin liquid, we included an additional next-nearest-neighbor SU(3) symmetric interactions along one of the two plaquette diagonal directions. This allows to tune the square-lattice model to the triangular-lattice model. By computing spin correlation functions and various order parameters, we find that the three-sublattice order may develop at infinitesimal additional new coupling, at least within the precision of our study. Compared with the previous findings that the nematic spin liquid is stable in extended parameter regions with additional couplings that respect the lattice symmetries of the square lattice, we argue that here the diagonal couplings, which frustrate the bipartite-lattice structure, rapidly suppress the two-sublattice fluctuations and the three-sublattice order thus wins. This numerical result is consistent with the conjecture that the nematic spin liquid emerges from the competition between two- and three-sublattice fluctuations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Dark Energy Survey Year 1 results: measurement of the galaxy angular power spectrum

We use data from the first-year observations of the DES collaboration to measure the galaxy angular power spectrum (APS), and search for its BAO feature. We test our methodology in a sample of 1800 DES Y1-like mock catalogues. We use the pseudo-C ℓ method to estimate the APS and the mock catalogues to estimate its covariance matrix. We use templates to model the measured spectra and estimate template parameters firstly from the C ℓ ’s of the mocks using two different methods, a maximum likelihood estimator and a Markov Chain Monte Carlo, finding consistent results with a good reduced χ 2 . Robustness tests are performed to estimate the impact of different choices of settings used in our analysis. Finally, we apply our method to a galaxy sample constructed from DES Y1 data specifically for LSS studies. This catalogue comprises galaxies within an effective area of 1318 deg 2 and 0.6 < z < 1.0. We find that the DES Y1 data favour a model with BAO at the 2.6σ C.L. However, the goodness of fit is somewhat poor, with χ 2 /(d.o.f.) = 1.49. We identify a possible cause showing that using a theoretical covariance matrix obtained from C ℓ ’s that are better adjusted to data results in an improved value of χ 2 /(dof) = 1.36 which is similar to the value obtained with the real-space analysis. Lastly, our results correspond to a distance measurement of D A (z eff = 0.81)/r d = 10.65 ± 0.49, consistent with the main DES BAO findings. This is a companion paper to the main DES BAO article showing the details of the harmonic space analysis.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Modeling and Analysis of FCM UN TRISO Fuel Using the PARFUME Code

The PARFUME (PARticle FUel ModEl) modeling code was used to assess the overall fuel performance of uranium nitride (UN) tri-structural isotropic (TRISO) ceramic fuel in the frame of the design and development of Fully Ceramic Matrix (FCM) fuel. A specific modeling of a TRISO particle with UN kernel was developed with PARFUME, and its behavior was assessed in irradiation conditions typical of a Light Water Reactor (LWR). The calculations were used to access the dimensional changes of the fuel particle layers and kernel, including the formation of an internal gap. The survivability of the UN TRISO particle was estimated depending on the strain behavior of the constituent materials at high fast fluence and burn-up. For nominal cases, internal gas pressure and representative thermal profiles across the kernel and layers were determined along with stress levels in the pyrolytic carbon (PyC) and silicon carbide (SiC) layers. These parameters were then used to evaluate fuel particle failure probabilities. Results of the study show that the survivability of UN TRISO fuel under LWR irradiation conditions might only be guaranteed if the kernel and PyC swelling rates are limited at high fast fluence and burn-up. These material properties are unknown at the irradiation levels expected to be reached by UN TRISO fuel in LWRs. Therefore, more effort is needed to determine them and positively conclude on the applicability of FCM fuel to LWRs.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2017

Data Synthesis and Data Assimilation at Global Change Experiments and Fluxnet Toward Improving Land Process Models

The project was conducted during the period from 7/1/2012 to 6/30/2017 with three major tasks: (1) data synthesis and development of data assimilation (DA) techniques to constrain modeled ecosystem feedback to climate change; (2) applications of DA techniques to improve process models at different scales from ecosystem to regions and the globe; and 3) improvements of modeling soil carbon (C) dynamics by land surface models. During this period, we have synthesized published data from soil incubation experiments (e.g., Chen et al., 2016; Xu et al., 2016; Feng et al., 2016), global change experiments (e.g., Li et al., 2013; Shi et al., 2015, 2016; Liang et al., 2016) and fluxnet (e.g., Niu et al., 2012., Xia et al., 2015; Li et al., 2016). These data have been organized into multiple data products and have been used to identify general mechanisms and estimate parameters for model improvement. We used the data sets that we collected and the DA techniques to improve model performance of both ecosystem models and global land models. The objectives are: 1) to improve model simulations of litter and soil carbon storage (e.g., Schädel et al., 2013; Hararuk and Luo, 2014; Hararuk et al., 2014; Liang et al., 2015); 2) to explore the effects of CO 2 , warming and precipitation on ecosystem processes (e.g., van Groenigen et al., 2014; Shi et al., 2015, 2016; Feng et al., 2017); and 3) to estimate parameters variability in different ecosystems (e.g., Li et al., 2016). We developed a traceability framework, which was based on matrix approaches and decomposed the modeled steady-state terrestrial ecosystem carbon storage capacity into four can trace the difference in ecosystem carbon storage capacity among different biomes to four traceable components: net primary productivity (NPP), baseline C residence times, environmental scalars and climate forcing (Xia et al., 2013). With this framework, we can diagnose the differences in modeled carbon storage across ecosystems, biomes, and models. This framework has been successfully implemented by several global land models, such as CABLE (Xia et al., 2013), LPJ-GUESS (Ahlström et al., 2015), CLM (Hararuk et al., 2014; Huang et al., 2017, submitted; Shi et al., 2017, submitted), and ORCHIDEE (Huang et al., 2017, unpublished). Moreover, we have identified the theoretical foundation of the determinants of transient C storage dynamics by adding another term, C storage potential, to the steady-state traceability framework (Luo et al., 2017). The theoretical foundation of transient C storage dynamics has been applied to develop a transient traceability framework to explore the traceable components of transient C storage dynamics responded to the rising CO 2 and climate change in the two contrasting ecosystem types Duke needleleaved forest and Harvard deciduous broadleaved forest (Jiang et al., 2017, in revision). Overall, with the data synthesis, data assimilation techniques, and the steady-state and transient traceability frameworks, we have greatly improved land process models for predicting responses and feedback of terrestrial C dynamics to global change. The matrix approaches has the potential to be applied in theoretical research on nitrogen and phosphorus cycle, and therefore, the coupling of carbon-nitrogen-phosphorus.

58 GEOSCIENCES↗
NASA NTRS2017

Interaction of Delaminations and Matrix Cracks in a CFRP Plate, Part II: Simulation Using an Enriched Shell Finite Element Model

Numerical simulations are presented of a recently developed test which creates multiple delaminations in a CFRP laminate specimen that grow and interact via transverse matrix cracks. A novel shell element enriched with the Floating Node Method, and a damage algorithm based on the Virtual Crack Closure Technique, were used to successfully simulate the tests. Additionally, a 3D high mesh fidelity model based on cohesive zones and continuum damage mechanics was used to simulate the tests and act as a representative of other similar state-of-the-art high mesh fidelity modeling techniques to compare to the enriched shell element. The enriched shell and high mesh fidelity models had similar levels of accuracy and generally matched the experimental data. With runtimes of 36 minutes for the shell model and 55 hours for the high mesh fidelity model, the shell model is 92 times faster than the high- fidelity simulation.

Mark W McElroy↗
DOE OSTI2015

High-fidelity micro-scale modeling of the thermo-visco-plastic behavior of carbon fiber polymer matrix composites

An experimentally validated micro-scale analysis of the visco-thermo-mechanical behavior of polymer matrix composites under different loads is proposed. A new constitutive law for the matrix material is developed taking into account the pressure dependence of the material as well as strain-rate and temperature dependence. Capturing the matrix behavior under multi-axial stress states is concluded to be essential to accurately predict the composite material behavior, even when considering simple load cases such as transverse compression and/or shear. Without any calibration procedure at the composite level, good agreement with the experimental data is observed for different loading conditions, including strain-rate dependency. Using this validated micro-scale model, a three-dimensional simulation of the formation of a kink band under longitudinal compression of the composite is conducted. A new evidence at micro-scale is found supporting the hypothesis that shear stresses transferred between fibers and matrix are particularly important in the formation of the kink band.

36 MATERIALS SCIENCE↗
NASA NTRS1991

Thermal histories of CO3 chondrites - Application of olivine diffusion modelling to parent body metamorphism

The petrologic sequence observed in the CO3 chondrite group has been suggested to be the result of thermal metamorphism on a parent body. A model developed to examine the possibility that chondrule and matrix olivines equilibrated in situ, during parent body metamorphism is presented. The model considers Fe-Mg interdiffusion between chondrule and matrix olivines. Zoning profiles comparable to those observed in chondrule olivines from partially equilibrated members of the series are reproduced successfully. Metamorphism of CO3 chondrites on a parent body is therefore a viable model for the observed equilibration. Results indicate that peak metamorphic temperatures experienced by the CO3 chondrites were around 500 C, and that the range of peak temperatures between unequilibrated and equilibrated subtypes was relatively narrow, around 100 C.

Jones, Rhian H.↗
DOE OSTI2025

Superconductivity in the lightly doped Hubbard model on the cylindrical honeycomb lattice

Here, we have performed large-scale density-matrix renormalization group studies of the lightly doped Hubbard model on the honeycomb lattice on long three- and four-leg cylinders. We find that the ground state of the system upon lightly doping is consistent with that of a superconducting state with coexisting quasi-long-range superconducting and charge density wave orders. Both the superconducting and charge density wave correlations decay as a power law at long distances with corresponding exponents 𝐾 𝑠⁢𝑐 < 2 and 𝐾 𝑐 < 2. On the contrary, the spin-spin and single-particle correlations decay exponentially, although with relatively long correlation lengths.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

INITIAL ANALYSIS OF TRANSIENT POWER TIME LAG DUE TO HETEROGENEITY WITHIN THE TREAT FUEL MATRIX.

The topic Nuclear Safety encompasses a broad spectrum of focal areas within the nuclear industry; one specific aspect centers on the performance and integrity of nuclear fuel during a reactivity insertion accident (RIA). This specific accident has proven to be fundamentally difficult to theoretically characterize due to the numerous empirically driven characteristics that quantify the fuel and reactor performance. The Transient Reactor Test (TREAT) facility was designed and operated to better understand fuel behavior under extreme (i.e. accident) conditions; it was shutdown in 1994. Recently, efforts have been underway to commission the TREAT facility to continue testing of advanced accident tolerant fuels (i.e. recently developed fuel concepts). To aid in the restart effort, new simulation tools are being used to investigate the behavior of nuclear fuels during facility’s transient events. This study focuses specifically on the characterizing modeled effects of fuel particles within the fuel matrix of the TREAT. The objective of this study was to (1) identify the impact of modeled heterogeneity within the fuel matrix during a transient event, and (2) demonstrate acceptable modeling processes for the purpose of TREAT safety analyses, specific to fuel matrix and particle size. Hypothetically, a fuel that is dominantly heterogeneous will demonstrate a clearly different temporal heating response to that of a modeled homogeneous fuel. This time difference is a result of the uniqueness of the thermal diffusivity within the fuel particle and fuel matrix. Using MOOSE/BISON to simulate the temperature time-lag effect of fuel particle diameter during a transient event, a comparison of the average graphite moderator temperature surrounding a spherical particle of fuel was made for both types of fuel simulations. This comparison showed that at a given time and with a specific fuel particle diameter, the fuel particle (heterogeneous) simulation and the homogeneous simulation were related by a multiplier relative to the average moderator temperature. As time increases the multiplier is comparable to the factor found in a previous analytical study from literature. The implementation of this multiplier and the method of analysis may be employed to remove assumptions and increase fidelity for future research on the effect of fuel particles during transient events.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2018

Mice overexpressing latent matrix metalloproteinase-2 develop lung emphysema after short-term exposure to cigarette smoke extract

Highlights: • Protease/anti-protease imbalance is an important event in the pathogenesis of lung emphysema. • Long time exposure to cigarette smoke is generally required to induce mouse model of lung emphysema. • Mice overexpressing matrix metalloproteinase-2 develop lung emphysema after short exposure to cigarette smoke extract. Chronic obstructive pulmonary disease is the major growing cause of mortality and morbidity worldwide, and it is going to become the third most common cause of death by 2020. Chronic obstructive pulmonary disease is pathologically characterized by lung emphysema and small airway inflammation. Animal models are very important to get insights into the disease pathogenesis but current models of chronic obstructive pulmonary disease take a long time to develop. The need of a new model is compelling. In the present study we focus on the role of matrix metalloproteinases in the pathogenesis of chronic obstructive pulmonary disease and hypothesized that lung overexpression of latent matrix metalloproteinases-2 would allow the development of emphysema after short-term exposure to cigarette smoke extract inhalation. Human latent matrix metalloproteinases-2 transgenic mouse expressing high level of the protein in the lungs and wild type mouse were exposed to aerosolized cigarette smoke extract for two weeks. Transgenic mice showed significant lung emphysematous changes, increased infiltration of inflammatory cells and enhanced lung concentrations of inflammatory cytokines in the lungs compared to their wild type counterparts after inhalation of cigarette smoke extract. This novel mouse model will be a very useful tool for evaluating the mechanistic pathways and for development of novel therapies in cigarette smoke-associated lung emphysema.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1992

Towards the development of micromechanics equations for ceramic matrix composites via fiber substructuring

A generic unit cell model which includes a unique fiber substructuring concept is proposed for the development of micromechanics equations for continuous fiber reinforcement ceramic composites. The unit cell consists of three constituents: fiber, matrix, and an interphase. In the present approach, the unit cell is further subdivided into several slices and the equations of micromechanics are derived for each slice. These are subsequently integrated to obtain ply level properties. A stand alone computer code containing the micromechanics model as a module is currently being developed specifically for the analysis of ceramic matrix composites. Towards this development, equivalent ply property results for a SiC/Ti-15-3 composite with 0.5 fiber volume ratio are presented and compared with those obtained from customary micromechanics models to illustrate the concept. Also, comparisons with limited experimental data for the ceramic matrix composite, SiC/RBSN (Reaction Bonded Silicon Nitride) with a 0.3 fiber volume ratio are given to validate the concepts.

Murthy, P. L. N.↗