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 613 records · Page 34

Local stress within a granular molecular solvent matrix, a mechanism for individual ion hydration

The activity of water is important in many natural systems and directed processes. To better understand water activity, the hydration of individual ions is approximated for a recently developed speciation-based solution model. Here a matrix of values for hydration of complete salts (anion and cation) was used as source data to calculate individual ion hydration values. Some electrolytes were excluded from calculations of individual ion hydration due to the prevalence of ion pairing and resulting liberation of water upon ion pairing. Evaluation of the periodic trends emerging from individual ion hydration values suggests that ionic volume, relative to the volume of molecular water, dictates both modeled hydration requirements (related to vapor-liquid-equilibrium-derived activity) as well as properties historically described by hydrodynamic radius (ion diffusion and solution viscosity). When water is conceptualized as a granular three-dimensional matrix, solutes may function as substitutional defects (dilation or constriction) inducing local compressive or tensile stress. These results have implications on the mechanism of the Hofmeister series and series reversals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Superelement methods in high temperature metal matrix composites

A study into fiber fracture and debonding in metal matrix composites is conducted using the finite element method. The superelement finite element technique was used to model a metal matrix composite under various loading condition and with varying degrees of fiber debonding. The use of superelement saved many man hours by allowing for alteration of only the primary superelement to manipulate partial bonding for the entire model. The composite's material properties were calculated and the effects of fiber debonding on these properties were noted. The internal stress state of the composite while under various loads was also studied. Special interest was devoted to the change in stress state as a result of increasing fiber debonding.

Trowbridge, Daniel↗

Internal state variable approach for predicting stiffness reductions in fibrous laminated composites with matrix cracks

A mathematical model utilizing the internal state variable concept is proposed for predicting the upper bound of the reduced axial stiffnesses in cross-ply laminates with matrix cracks. The axial crack opening displacement is explicitly expressed in terms of the observable axial strain and the undamaged material properties. A crack parameter representing the effect of matrix cracks on the observable axial Young's modulus is calculated for glass/epoxy and graphite/epoxy material systems. The results show that the matrix crack opening displacement and the effective Young's modulus depend not on the crack length, but on its ratio to the crack spacing.

Lee, Jong-Won↗

General structure of democratic mass matrix of quark sector in E{sub 6} model

An extension of the Standard Model (SM) fermion sector, which is inspired by the E{sub 6} Grand Unified Theory (GUT) model, might be a good candidate to explain a number of unanswered questions in SM. Existence of the isosinglet quarks might explain great mass difference of bottom and top quarks. Also, democracy on mass matrix elements is a natural approach in SM. In this study, we have given general structure of Democratic Mass Matrix (DMM) of quark sector in E6 model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Pairing tendencies in the doped Kitaev-Heisenberg model

Here, we study the impact of hole doping on the Kitaev-Heisenberg model on the honeycomb lattice. We investigate the pairing tendencies and correlation functions in the framework of a t - J - K model using density matrix renormalization group calculations on three-leg cylinders. In the case of the pure Kitaev model, which realizes a quantum spin-liquid phase at half-filling, we find that binding of two holes only occurs at low values of the hopping, where the holes are slow. We have theoretically verified that pair formation occurs in the limit of immobile holes, where the pure Kitaev model remains exactly solvable. When we instead fix the hopping at an intermediate, more realistic, value, and vary the Heisenberg and Kitaev interaction strengths, we find pairing tendencies only in the Néel phase. This is in contrast to prior mean-field calculations, highlighting the importance of accounting for the kinetic energy of dopants in generalized Kitaev models. Interestingly, we also find signatures of pair-density wave formation over the studied range of model parameters, namely, a periodic modulation of the charge density as well as the spin-spin and pair-pair correlations in real space. Moreover, we present a comparative study of the different correlations as a function of doping. We finally discuss the potential for experimentally observing the studied physics in quantum materials and heterostructures.

36 MATERIALS SCIENCE↗

Modeling Host-Pathogen Interactions in the Context of the Microenvironment: Three-Dimensional Cell Culture Comes of Age

Tissues and organs provide the structural and biochemical landscapes upon which microbial pathogens and commensals function to regulate health and disease. While flat two-dimensional (2-D) monolayers composed of a single cell type have provided important insight into understanding host-pathogen interactions and infectious disease mechanisms, these reductionist models lack many essential features present in the native host microenvironment that are known to regulate infection, including three-dimensional (3-D) architecture, multicellular complexity, commensal microbiota, gas exchange and nutrient gradients, and physiologically relevant biomechanical forces (e.g., fluid shear, stretch, compression). Here, a major challenge in tissue engineering for infectious disease research is recreating this dynamic 3-D microenvironment (biological, chemical, and physical/mechanical) to more accurately model the initiation and progression of host-pathogen interactions in the laboratory. Here we review selected 3-D models of human intestinal mucosa, which represent a major portal of entry for infectious pathogens and an important niche for commensal microbiota. We highlight seminal studies that have used these models to interrogate host-pathogen interactions and infectious disease mechanisms, and we present this literature in the appropriate historical context. Models discussed include 3-D organotypic cultures engineered in the rotating wall vessel (RWV) bioreactor, extracellular matrix (ECM)-embedded/organoid models, and organ-on-a-chip (OAC) models. Collectively, these technologies provide a more physiologically relevant and predictive framework for investigating infectious disease mechanisms and antimicrobial therapies at the intersection of the host, microbe, and their local microenvironments.

3-D↗

Accurate Holdup Calculations with Predictive Modeling & Data Integration

In facilities that process special nuclear material (SNM) it is important to account accurately for the fissile material that enters and leaves the plant. Although there are many stages and processes through which materials must be traced and measured, the focus of this project is material that is “held-up” in equipment, pipes, and ducts during normal operation and that can accumulate over time into significant quantities. Accurately estimating the holdup is essential for proper SNM accounting (vis-à-vis nuclear non-proliferation), criticality and radiation safety, waste management, and efficient plant operation. Usually it is not possible to directly measure the holdup quantity and location, so these must be inferred from measured radiation fields, primarily gamma and less frequently neutrons. Current methods to quantify holdup, i.e. Generalized Geometry Holdup (GGH), primarily rely on simple source configurations and crude radiation transport models aided by ad hoc correction factors. This project seeks an alternate method of performing measurement-based holdup calculations using a predictive model that employs state-of-the-art radiation transport codes capable of accurately simulating such situations. Inverse and data assimilation methods use the forward transport model to search for a source configuration that best matches the measured data and simultaneously provide an estimate of the level of confidence in the correctness of such configuration. In this work the holdup problem is re-interpreted as an inverse problem that is under-determined, hence may permit multiple solutions. A probabilistic approach is applied to solving the resulting inverse problem. This approach rates possible solutions according to their plausibility given the measurements and initial information. This is accomplished through the use of Bayes’ Theorem that resolves the issue of multiple solutions by giving an estimate of the probability of observing each possible solution. To use Bayes’ Theorem, one must have a model y(x) that maps the state variables x (the solution in this case) to the measurements y. In this case, the unknown state variables are the configuration and composition of the heldup SNM. The measurements are the detector readings. Thus, the natural model is neutral-particle radiation transport where a wealth of computational tools exists for performing these simulations accurately and efficiently. The combination of predictive model and Bayesian inference forms the Data Integration with Modeled Predictions (DIMP) method that serves as foundation for this project. The cost functional describing the model-to-data misfit is computed via a norm created by the inverse of the covariance matrix of the model parameters and responses. Since the model y(x) for the holdup problem is nonlinear, a nonlinear optimization on Q is conducted via Newton-type iterative methods to find the optimal values of the model parameters x. This project comprised a collaboration between NC State University (NCSU), the University of South Carolina (USC), and Oak Ridge National Laboratory (ORNL). The project was originally proposed in seven main tasks with an eighth contingency task to be performed if time and funding permitted; in fact time did not permit commencement of the contingency task and it was not performed. The remaining tasks involved holdup analysis with gamma detection strategies and separately with neutrons based on coincidence counting. Early in the project, and upon consultation with experts in coincidence counting it became evident that this approach is not viable for holdup applications and this task was replaced with an alternative, but valuable investigation that was carried out by the USC partner. Nevertheless, the experimental 4 measurements at ORNL of both gamma and neutron sources for the purpose of constructing Detector Response Functions (DRFs) with the associated uncertainties were indeed completed.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

A model of the thermal-spike mechanism in graphite/epoxy laminates

The influence of a thermal spike on a moisture-saturated graphite/epoxy composite was studied in detail. A single thermal spike from 25 C to 132 C was found to produce damage as evidenced by a significant increase in the level of moisture saturation in the composite. Approximately half of this increase remained after a vacuum anneal at 150 C for 7 days, suggesting the presence of an irreversible damage component. Subsequent thermal spikes created less and less additional moisture absorption, with the cumulative effect being a maximum or limiting moisture capacity of the composite. These observations are explained in terms of a model previously developed to explain the reverse thermal effect of moisture absorption in epoxy and epoxy matrix composites. This model, based on the inverse temperature dependence of free volume, contributes an improved understanding of thermal-spike effects in graphite/epoxy composites.

Adamson, M. J.↗

Modeling 2-D jets impinging on Stirling regenerators

The extent to which flow leaving Stirling coolers or heaters in the form of high-velocity jets penetrate the regenerator matrix is visually modeled using a computer program. Two-dimensional laminar jets are shown impinging on regenerator samples of variable permeability ranging from no matrix at all to matrices dense enough to stop the jet dead on. The results lend credibility to a simple tension for flow uniformity as a function of penetration depth.

Gedeon, David↗

Slushy weightings for the optimal pilot model

A pilot model is described which accounts for the effect of motion cues in a well defined visual tracking task. The effect of visual and motion cues are accounted for in the model in two ways. First, the observation matrix in the pilot model is structured to account for the visual and motion inputs presented to the pilot. Secondly, the weightings in the quadratic cost function associated with the pilot model are modified to account for the pilot's perception of the variables he considers important in the task. Analytic results obtained using the pilot model are compared to experimental results and in general good agreement is demonstrated. The analytic model yields small improvements in tracking performance with the addition of motion cues for easily controlled task dynamics and large improvements in tracking performance with the addition of motion cues for difficult task dynamics.

Dillow, J. D.↗

A NASTRAN DMAP alter for the coupling of modal and physical coordinate substructures

A method is described to derive a generalized coordinate model consisting of flexible cantilever modes and rigid body modes from a physical coordinate model using a direct matrix abstraction procedure alter to rigid format 3. This model can readily be coupled to other substructures using modal synthesis techniques. It allows the use of a reduced size model for structural analyses while maintaining the capability of recovering the accelerations, forces, stresses, etc., from the original, large, complex model. This output recovery is accomplished with the use of a loads transformation matrix which relates the output parameters to the modal coordinate accelerations. In addition, a method is described to synthesize structural models consisting of hybrid coordinates for use in dynamic response analyses where one structure is described using physical coordinates and the other using generalized modal coordinates.

Wilson, T. L.↗

Finite Element Model for Failure Study of Two-Dimensional Triaxially Braided Composite

A new three-dimensional finite element model of two-dimensional triaxially braided composites is presented in this paper. This meso-scale modeling technique is used to examine and predict the deformation and damage observed in tests of straight sided specimens. A unit cell based approach is used to take into account the braiding architecture as well as the mechanical properties of the fiber tows, the matrix and the fiber tow-matrix interface. A 0 deg / plus or minus 60 deg. braiding configuration has been investigated by conducting static finite element analyses. Failure initiation and progressive degradation has been simulated in the fiber tows by use of the Hashin failure criteria and a damage evolution law. The fiber tow-matrix interface was modeled by using a cohesive zone approach to capture any fiber-matrix debonding. By comparing the analytical results to those obtained experimentally, the applicability of the developed model was assessed and the failure process was investigated.

Li, Xuetao↗

Scalable ab initio fragmentation methods based on a truncated expansion of the non-orthogonal molecular orbital model

An alternative formulation of the non-orthogonal molecular orbital model of electronic structure theory is developed based on the expansion of the inverse molecular orbital overlap matrix. From this model, a hierarchy of ab initio fragment-based quantum chemistry methods, referred to as the nth-order expanded non-orthogonal molecular orbital methods, are developed using a minimal number of approximations, each of which is frequently employed in intermolecular interaction theory. These novel methods are compared to existing fragment-based quantum chemistry methods, and the implications of those significant differences, where they exist, between the methods developed herein and those already existing methods are examined in detail. Computational complexities and theoretical scaling are also analyzed and discussed. Future extensions for the hierarchy of methods, to account for additional intrafragment and interfragment interactions, are outlined.

Westheimer, Bryce M. (ORCID:0000000267645714)↗

Automated RF Phase Adjustment for Beam Stabilization in the Fermilab Linac

The Fermilab Linac experiences longitudinal beam phase drift, leading to increased particle loss, conventionally corrected through labor-intensive manual RF adjustments. This project explores machine learning-based automation for drift correction, employing a prototype-based classification approach. Our model utilizes a 34-dimensional feature set (RF settings and BPM readings) and leverages a 7x27 response matrix for system modeling. To overcome limited real-world data, we generate synthetic data, enhancing model training and generalizability. Custom loss functions, including a surrogate energy-consistent loss and a temporal smoothness constraint, ensure physically plausible drift predictions. The goal is a robust system for autonomous phase adjustments, ensuring stable beam acceleration and reduced manual intervention.

Chichili, R. R. [Illinois U., Chicago]↗

Automated RF Phase Adjustment for Beam Stabilization in the Fermilab Linac

The Fermilab Linac experiences longitudinal beam phase drift, leading to increased particle loss, conventionally cor- rected through labor-intensive manual RF adjustments. This project explores machine learning-based automation for drift correction, employing a prototype-based classification approach. Our model utilizes a 34-dimensional feature set (RF settings and BPM readings) and leverages a 7x27 response matrix for system modeling. To overcome limited real-world data, we generate synthetic data, enhancing model training and generalizability. Custom loss functions, including a sur- rogate energy-consistent loss and a temporal smoothness constraint, ensure physically plausible drift predictions. The goal is a robust system for autonomous phase adjustments, ensuring stable beam acceleration and reduced manual intervention.

Chichili, R. R. [U. Illinois, Chicago]↗

A conceptual model of the thermal-spike mechanism in graphte/epoxy laminates

The influence of a thermal spike on a moisture-saturated graphite/epoxy composite was studied in detail. A single thermal spike from 25 to 132 C was found to produce damage, as evidenced by a significant increase in the level of moisture saturation in the composite. Approximately half of this increase remained after a vacuum-anneal at 150 C for seven days, suggesting the presence of an irreversible damage component. Subsequent thermal spikes created less and less additional moisture absorption, with the cumulative effect being a maximum or limiting moisture capacity of the composite. These observations are explained in terms of a model previously developed to explain the reverse thermal effect of moisture absorption in epoxy and epoxy-matrix composites. This model, based on the inverse temperature dependence of free volume, improves current understanding of thermal spike effects in graphite/epoxy composites.

Adamson, M. J.↗

A model of filament-wound thin cylinders

A model was developed for simulating he manufacturing process of filament-wound cylinders made of a thermoset matrix composite. The model relates the process variables (winding speed, fiber tension, applied temperature) to the parameters characterizing the composite cylinder and the mandrel. The model is applicable to cylinders for which the diameter is large compared to the wall thickness. The model was implemented by a user-friendly computer code suitable for generating numerical results.

Calius, Emilio P.↗

Single-phase power distribution system power flow and fault analysis

Alternative methods for power flow and fault analysis of single-phase distribution systems are presented. The algorithms for both power flow and fault analysis utilize a generalized approach to network modeling. The generalized admittance matrix, formed using elements of linear graph theory, is an accurate network model for all possible single-phase network configurations. Unlike the standard nodal admittance matrix formulation algorithms, the generalized approach uses generalized component models for the transmission line and transformer. The standard assumption of a common node voltage reference point is not required to construct the generalized admittance matrix. Therefore, truly accurate simulation results can be obtained for networks that cannot be modeled using traditional techniques.

Halpin, S. M.↗