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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Micromechanical modeling of damage growth in titanium based metal-matrix composites

The thermomechanical behavior of continuous-fiber reinforced titanium based metal-matrix composites (MMC) is studied using the finite element method. A thermoviscoplastic unified state variable constitutive theory is employed to capture inelastic and strain-rate sensitive behavior in the Timetal-21s matrix. The SCS-6 fibers are modeled as thermoplastic. The effects of residual stresses generated during the consolidation process on the tensile response of the composites are investigated. Unidirectional and cross-ply geometries are considered. Differences between the tensile responses in composites with perfectly bonded and completely debonded fiber/matrix interfaces are discussed. Model simulations for the completely debonded-interface condition are shown to correlate well with experimental results.

Sherwood, James A.↗
NASA NTRS2002

System Model for MEMS based Laser Ultrasonic Receiver

A need has been identified for more advanced nondestructive Evaluation technologies for assuring the integrity of airframe structures, wiring, etc. Laser ultrasonic inspection instruments have been shown to detect flaws in structures. However, these instruments are generally too bulky to be used in the confined spaces that are typical of aerospace vehicles. Microsystems technology is one key to reducing the size of current instruments and enabling increased inspection coverage in areas that were previously inaccessible due to instrument size and weight. This paper investigates the system modeling of a Micro OptoElectroMechanical System (MOEMS) based laser ultrasonic receiver. The system model is constructed in software using MATLAB s dynamical simulator, Simulink. The optical components are modeled using geometrical matrix methods and include some image processing. The system model includes a test bench which simulates input stimuli and models the behavior of the material under test.

Wilson, William C.↗
DOE OSTI2023

Modeling Workloads of a Linear Electromagnetic Code for Load Balancing Matrix Assembly

This report presents our work to model the workloads of a linear electromagnetic application based on the method of moments in the frequency domain to effectively load balance the matrix assembly. This application is particularly challenging to load balance due to its lack of persistent iterative behavior, its operation under tight memory constraint (where the matrix may fill 80% of memory on each node), and the algorithmic complexity of the computational method. This report describes the first step in our work to apply an inspector-executor approach for load balancing workloads where key parameters are exposed during the inspector phase and a pre-trained model is applied to predict relative task weights for the load balancer.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Investigation of Coupled Processes in Fractures and the Bordering Matrix via a Micro‐Continuum Reactive Transport Model

Abstract In multi‐mineral fractured rocks, the altered porous layer on the fracture surface resulting from preferential dissolution of the fast‐reacting minerals can have profound impacts on subsequent chemical‐physical alteration of the fractures. This study adopts the micro‐continuum approach to provide further understanding of reactive transport processes in the altered layer (AL), and mass exchanges with the bordering matrix and fracture. The modeling framework couples the Darcy‐Brinkman‐Stokes (DBS) solver in COMSOL Multiphysics and the geochemical modeling capability of CrunchFlow. Three‐dimensional steady state simulations with systematically varied chemical‐physical parameters of the AL were performed to examine the impacts of individual factors and processes. Our simulation results confirm previous observations that dissolution of the fast‐reacting mineral (i.e., calcite) is largely controlled by diffusion across the AL. We also show that dissolution of the slow‐reacting mineral (i.e., dolomite), which controls AL development and fracture enlargement, increases with surface area and has a complex dependence on different local rate‐limiting processes. In particular, advection can result in evident spatial variations in the local dissolution rates of dolomite, although it does not affect the bulk chemistry significantly. The difference in the spatial patterns between simulations with and without advection in the AL is more noticeable in the locations with smaller apertures, with up to 20% difference in local reaction rates. Therefore, it is important to include a full depiction of advection, diffusion, and reactions for accurately capturing local dynamics that control long‐term fracture evolution.

58 GEOSCIENCES↗
DOE OSTI2016

New type of hill-top inflation

We suggest a new type of hill-top inflation originating from the initial conditions in the form of the microcanonical density matrix for the cosmological model with a large number of quantum fields conformally coupled to gravity. Initial conditions for inflation are set up by cosmological instantons describing underbarrier oscillations in the vicinity of the inflaton potential maximum. These periodic oscillations of the inflaton field and cosmological scale factor are obtained within the approximation of two coupled oscillators subject to the slow roll regime in the Euclidean time. This regime is characterized by rapid oscillations of the scale factor on the background of a slowly varying inflaton, which guarantees smallness of slow roll parameters ϵ and η of the following inflation stage. A hill-like shape of the inflaton potential is shown to be generated by logarithmic loop corrections to the tree-level asymptotically shift-invariant potential in the non-minimal Higgs inflation model and R{sup 2}-gravity. The solution to the problem of hierarchy between the Planckian scale and the inflation scale is discussed within the concept of conformal higher spin fields, which also suggests the mechanism bringing the model below the gravitational cutoff and, thus, protecting it from large graviton loop corrections.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2016

$B^0_{(s)}$-mixing matrix elements from lattice QCD for the Standard Model and beyond

We calculate—for the first time in three-flavor lattice QCD—the hadronic matrix elements of all five local operators that contribute to neutral B 0 - and B s -meson mixing in and beyond the Standard Model. We present a complete error budget for each matrix element and also provide the full set of correlations among the matrix elements. We also present the corresponding bag parameters and their correlations, as well as specific combinations of the mixing matrix elements that enter the expression for the neutral B-meson width difference. We obtain the most precise determination to date of the SU(3)-breaking ratio ξ=1.206(18)(6), where the second error stems from the omission of charm-sea quarks, while the first encompasses all other uncertainties. The threefold reduction in total uncertainty, relative to the 2013 Flavor Lattice Averaging Group results, tightens the constraint from B mixing on the Cabibbo-Kobayashi-Maskawa (CKM) unitarity triangle. Our calculation employs gauge-field ensembles generated by the MILC Collaboration with four lattice spacings and pion masses close to the physical value. We use the asqtad-improved staggered action for the light-valence quarks and the Fermilab method for the bottom quark. We use heavy-light meson chiral perturbation theory modified to include lattice-spacing effects to extrapolate the five matrix elements to the physical point. We combine our results with experimental measurements of the neutral B-meson oscillation frequencies to determine the CKM matrix elements |V td | = 8.00(34)(8)×10 -3 , |V ts | = 39.0(1.2)(0.4)×10 -3 , and |V td /V ts | = 0.2052(31)(10), which differ from CKM-unitarity expectations by about 2σ. In addition, these results and others from flavor-changing-neutral currents point towards an emerging tension between weak processes that are mediated at the loop and tree levels.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS2017

Modeling of Different Fiber Type and Content SiC/SiC Minicomposites Creep Behavior

Silicon Carbide based Ceramic Matrix Composites (CMCs) are attractive materials for use in high-temperature applications in the aerospace and nuclear industries. However, creep damage mechanism in CMCs is the most dominant mechanism at elevated temperatures. Consequently, the tensile creep behavior of Hi-Nicalon, Hi-Nicalon Type S SiC fibers and Chemical vapor infiltrated Silicon Carbide matrix (CVI-SiC) were characterized and creep parameters were extracted from creep experiments. Some fiber creep tests were performed in inert environment at 1200 C on individual fibers. Creep behavior of different fiber content pristine and precracked Hi-Nicalon and Hi-Nicalon Type S reinforced minicomposites with BN interphases and CVI-SiC matrix were then modelled using the creep data found in this study and the literature and compared with creep experiments results for the pristine and precracked Hi-Nicalon and Hi-Nicalon Type S minicomposites. Finally, the effects of load-sharing and matrix cracking on CMC creep behavior will be discussed.

ceramic matrix minicomposites↗
NASA NTRS2003

Modeling of 3-D Woven Ceramic Matrix Composites

Three different approaches are being pursued at the NASA Glenn Research Center to predict the nanostructural behavior of three-dimensional woven ceramic matrix composites. These are: a micromechanics-based approach using W-CEMCAN (Woven Ceramic Matrix Composite Analyzer), a laminate analogy method and a structural frame approach (based on the finite element method). All three techniques are applied to predict the thermomechanical properties of a three-dimensional woven angle interlock C/SiC composite. The properties are predicted for room temperature and 1100 C and the predicted properties are compared to measurements. General observations regarding the three approaches for three-dimensional composite modeling are discussed.

Murthy, Pappu L. N.↗
DOE OSTI2015

Shared Dosimetry Error in Epidemiological Dose-Response Analyses

Radiation dose reconstruction systems for large-scale epidemiological studies are sophisticated both in providing estimates of dose and in representing dosimetry uncertainty. For example, a computer program was used by the Hanford Thyroid Disease Study to provide 100 realizations of possible dose to study participants. The variation in realizations reflected the range of possible dose for each cohort member consistent with the data on dose determinates in the cohort. Another example is the Mayak Worker Dosimetry System 2013 which estimates both external and internal exposures and provides multiple realizations of "possible" dose history to workers given dose determinants. This paper takes up the problem of dealing with complex dosimetry systems that provide multiple realizations of dose in an epidemiologic analysis. In this paper we derive expected scores and the information matrix for a model used widely in radiation epidemiology, namely the linear excess relative risk (ERR) model that allows for a linear dose response (risk in relation to radiation) and distinguishes between modifiers of background rates and of the excess risk due to exposure. We show that treating the mean dose for each individual (calculated by averaging over the realizations) as if it was true dose (ignoring both shared and unshared dosimetry errors) gives asymptotically unbiased estimates (i.e. the score has expectation zero) and valid tests of the null hypothesis that the ERR slope β is zero. Although the score is unbiased the information matrix (and hence the standard errors of the estimate of β) is biased for β≠0 when ignoring errors in dose estimates, and we show how to adjust the information matrix to remove this bias, using the multiple realizations of dose. Use of these methods for several studies, including the Mayak Worker Cohort and the U.S. Atomic Veterans Study, is discussed.

epidemiological analysis↗
DOE OSTI2014

TCDD induces dermal accumulation of keratinocyte-derived matrix metalloproteinase-10 in an organotypic model of human skin

The epidermis of skin is the first line of defense against the environment. A three dimensional model of human skin was used to investigate tissue-specific phenotypes induced by the environmental contaminant, 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD). Continuous treatment of organotypic cultures of human keratinocytes with TCDD resulted in intracellular spaces between keratinocytes of the basal and immediately suprabasal layers as well as thinning of the basement membrane, in addition to the previously reported hyperkeratinization. These tissue remodeling events were preceded temporally by changes in expression of the extracellular matrix degrading enzyme, matrix metalloproteinase-10 (MMP-10). In organotypic cultures MMP-10 mRNA and protein were highly induced following TCDD treatment. Q-PCR and immunoblot results from TCDD-treated monolayer cultures, as well as indirect immunofluorescence and immunoblot analysis of TCDD-treated organotypic cultures, showed that MMP-10 was specifically contributed by the epidermal keratinocytes but not the dermal fibroblasts. Keratinocyte-derived MMP-10 protein accumulated over time in the dermal compartment of organotypic cultures. TCDD-induced epidermal phenotypes in organotypic cultures were attenuated by the keratinocyte-specific expression of tissue inhibitor of metalloproteinase-1, a known inhibitor of MMP-10. These studies suggest that MMP-10 and possibly other MMP-10-activated MMPs are responsible for the phenotypes exhibited in the basement membrane, the basal keratinocyte layer, and the cornified layer of TCDD-treated organotypic cultures. Our studies reveal a novel mechanism by which the epithelial–stromal microenvironment is altered in a tissue-specific manner thereby inducing structural and functional pathology in the interfollicular epidermis of human skin. - Highlights: • TCDD causes hyperkeratosis and basement membrane changes in a model of human skin. • TCDD induces MMP-10 expression in organotypic cultures of human keratinocytes. • Keratinocyte-expressed MMP-10 accumulates in the dermal compartment. • Keratinocyte K14 promoter-driven TIMP-1 expression ablates TCDD-induced phenotypes.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Predicting ductility of Mg/SiCp nanocomposite under multiaxial loading conditions based on unit cell modeling

This article presents an investigation on the ductility of Mg/SiCp under nine different multiaxial loading conditions based on a unit cell (or called representative volumetric element, RVE) finite element model. The constitutive relationship of magnesium matrix is characterized by a model combined Lode angle dependent plasticity model with the CPB06 plasticity model, which is able to describe the tension-compression asymmetry in magnesium matrix. The debonding of the matrix-particle interface and matrix fracture are considered by the modified Mohr-Coulomb (MMC) fracture model. Simulation results show that the interface debonding and subsequent matrix fracture are responsible for the ductility and fracture of nanocomposite. The computational results with parametric studies are compared and well correlated with experimental data, and it finally gives a stress-based full fracture locus of Mg/SiCp nanocomposite.

42 ENGINEERING↗
DOE OSTI2016

Model Predictive Control of A Matrix-Converter Based Solid State Transformer for Utility Grid Interaction

The matrix converter solid state transformer (MC-SST), formed from the back-to-back connection of two three-to-single-phase matrix converters, is studied for use in the interconnection of two ac grids. The matrix converter topology provides a light weight and low volume single-stage bidirectional ac-ac power conversion without the need for a dc link. Thus, the lifetime limitations of dc-bus storage capacitors are avoided. However, space vector modulation of this type of MC-SST requires to compute vectors for each of the two MCs, which must be carefully coordinated to avoid commutation failure. An additional controller is also required to control power exchange between the two ac grids. In this paper, model predictive control (MPC) is proposed for an MC-SST connecting two different ac power grids. The proposed MPC predicts the circuit variables based on the discrete model of MC-SST system and the cost function is formulated so that the optimal switch vector for the next sample period is selected, thereby generating the required grid currents for the SST. Simulation and experimental studies are carried out to demonstrate the effectiveness and simplicity of the proposed MPC for such MC-SST-based grid interfacing systems.

Xue, Yaosuo↗
DOE OSTI2016

Role of vertex corrections in the matrix formulation of the random phase approximation for the multiorbital Hubbard model

In the framework of a multiorbital Hubbard model description of superconductivity, a matrix formulation of the superconducting pairing interaction that has been widely used is designed to treat spin, charge, and orbital fluctuations within a random phase approximation (RPA). In terms of Feynman diagrams, this takes into account particle-hole ladder and bubble contributions as expected. It turns out, however, that this matrix formulation also generates additional terms which have the diagrammatic structure of vertex corrections. Furthermore we examine these terms and discuss the relationship between the matrix-RPA superconducting pairing interaction and the Feynman diagrams that it sums.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Numerical Modeling of Fiber-Reinforced Metal Matrix Composite Processing by the Liquid Route: Literature Contribution

One of the technologies used to produce metal matrix composites (MMCs) is liquid route processing. One solution is to inject a liquid metal under pressure or at constant rate through a fibrous preform. This foundry technique overcomes the problem of the wettability of ceramic fibers by liquid metal. The liquid route can also be used to produce semiproducts by coating a filament with a molten metal. These processes involve physical phenomena combined with mass and heat transfer and phase change. The phase change phenomena related to solidification and also to the melting of the metal during the process notably result in modifications to the permeability of porous media, in gaps in impregnation, in the appearance of defects (porosities), and in segregation in the final product. In this article, we provide a state-of-the-art review of numerical models and simulation developed to study these physical phenomena involved in MMC processing by the liquid route.

36 MATERIALS SCIENCE↗
NASA NTRS1980

Mechanisms of fatigue damage in boron/aluminum composites

Tensile fatigue tests were conducted on several laminates of boron/aluminum (6061-0). In laminates with 0 deg fibers on the outside, an analysis that identifies "shakedown" conditions predicted the stress amplitude below which no fatigue damage accumulated. A fatigue damage accumulation model which relates matrix fatigue cracking and the overall laminate properties is described. A model for the saturation damage stage development is presented, that identical laminates, tested in directions 90 deg apart (such that one layup has 90 deg outer plies and the other 0 deg), have different fatigue behaviors due to the stacking sequence. The 90 deg plies on the surface develop cracks earlier than predicted by shakedown. An attempt was made to explain this stacking sequence effect. Variable load history effects on the fatigue damage response were investigated. Tests reveal that for a given stress ratio the specimen seeks the saturation damage state for the largest stress range to which it is subjected. It was also found that little damage is generated by shifting a given stress range down, whereas significant damage may be created by shifting it upward. The laminate stresses were always tensile.

Johnson, W. S.↗
NASA NTRS1989

Modeling and control system design and analysis tools for flexible structures

Described here are Boeing software tools used for the development of control laws of flexible structures. The Boeing Company has developed a software tool called Modern Control Software Package (MPAC). MPAC provides the environment necessary for linear model development, analysis, and controller design for large models of flexible structures. There are two features of MPAC which are particularly appropriate for use with large models: (1) numerical accuracy and (2) label-driven nature. With the first feature MPAC uses double precision arithmetic for all numerical operations and relies on EISPAC and LINPACK for the numerical foundation. With the second feature, all MPAC model inputs, outputs, and states are referenced by user-defined labels. This feature allows model modification while maintaining the same state, input, and output names. In addition, there is no need for the user to keep track of a model variable's matrix row and colunm locations. There is a wide range of model manipulation, analysis, and design features within the numerically robust and flexible environment provided by MPAC. Models can be built or modified using either state space or transfer function representations. Existing models can be combined via parallel, series, and feedback connections; and loops of a closed-loop model may be broken for analysis.

Anissipour, Amir A.↗
DOE OSTI2018

Nonequilibrium Steady-State Transport in Quantum Impurity Models: A Thermofield and Quantum Quench Approach Using Matrix Product States

Here, the numerical renormalization group (NRG) is tailored to describe interacting impurity models in equilibrium, but it faces limitations for steady-state nonequilibrium, arising, e.g., due to an applied bias voltage. We show that these limitations can be overcome by describing the thermal leads using a thermofield approach, integrating out high energy modes using NRG, and then treating the nonequilibrium dynamics at low energies using a quench protocol, implemented using the time-dependent density matrix renormalization group. This yields quantitatively reliable results for the current (with errors ≲3%) down to the exponentially small energy scales characteristic of impurity models. We present results of benchmark quality for the temperature and magnetic field dependence of the zero-bias conductance peak for the single-impurity Anderson model.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗