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,621 records · Page 90
All sourcesNASA NTRSDOE OSTI
NASA NTRS2001

Numerical Model of Flame Spread Over Solids in Microgravity: A Supplementary Tool for Designing a Space Experiment

The recently developed numerical model of concurrent-flow flame spread over thin solids has been used as a simulation tool to help the designs of a space experiment. The two-dimensional and three-dimensional, steady form of the compressible Navier-Stokes equations with chemical reactions are solved. With the coupled multi-dimensional solver of the radiative heat transfer, the model is capable of answering a number of questions regarding the experiment concept and the hardware designs. In this paper, the capabilities of the numerical model are demonstrated by providing the guidance for several experimental designing issues. The test matrix and operating conditions of the experiment are estimated through the modeling results. The three-dimensional calculations are made to simulate the flame-spreading experiment with realistic hardware configuration. The computed detailed flame structures provide the insight to the data collection. In addition, the heating load and the requirements of the product exhaust cleanup for the flow tunnel are estimated with the model. We anticipate that using this simulation tool will enable a more efficient and successful space experiment to be conducted.

Shih, Hsin-Yi↗
DOE OSTI2019

A ROBUST MECHANISTIC APPROACH TO PREDICTION OF DEPARTURE FROM NUCLEATE BOILING

departure from nucleate boiling (DNB) sets a regulatory limit for licensing pressurized water reactors (PWRs). Under such conditions, the heated surface is permanently blanketed by a vapor film, leading to a sharp deterioration of the heat transfer coefficient at the heater/coolant interface and an abrupt temperature rise. Unfortunately, the path for an accurate, robust prediction of DNB has been elusive due to lack of consensus on its triggering mechanism. This work reviews existing physics-driven modeling tools in the first place. An evolutionary channel-scale mechanistic model that leverages key assumptions in the relatively well-accepted mechanisms of liquid sublayer dryout and near-wall bubble crowding is then proposed. Detailed validation of the proposed model has demonstrated its improved predictive capabilities over previous data/physics-driven models with regard to an extensive DNB-specific CHF test matrix covering a wide range of flow conditions. The unique feature of the proposed model lies in its ability to predict DNB for different heater geometries (including round tube, annulus, and rod bundle), which is essential in deciphering fuel performance metrics from different facilities and reactor types. The proposed model will be further implemented in fuel performance codes to help improve modeling of transient DNB scenarios such as during a reactivity-initiated accident (RIA).

Zhao, Xingang↗
NASA NTRS2006

Adapting Covariance Propagation to Account for the Presence of Modeled and Unmodeled Maneuvers

This paper explores techniques that can be used to adapt the standard linearized propagation of an orbital covariance matrix to the case where there is a maneuver and an associated execution uncertainty. A Monte Carlo technique is used to construct a final orbital covariance matrix for a 'propagate-burn-propagate' process that takes into account initial state uncertainty and execution uncertainties in the maneuver magnitude. This final orbital covariance matrix is regarded as 'truth' and comparisons between it and three methods using modified linearized covariance propagation are made. The first method accounts for the maneuver by modeling its nominal effect within the state transition matrix but excludes the execution uncertainty by omitting a process noise matrix from the computation. In the second method, the maneuver is not modeled but the uncertainty in its magnitude is accounted for by the inclusion of a process noise matrix. In the third method, which is essentially a hybrid of the first two, the nominal portion of the maneuver is included via the state transition matrix while a process noise matrix is used to account for the magnitude uncertainty. Since this method also correctly accounts for the presence of the maneuver in the nominal orbit, it is the best method for applications involving the computation of times of closest approach and the corresponding probability of collision, Pc. However, applications for the two other methods exist and are briefly discussed. Despite the fact that the process model ('propagate-burn-propagate') that was studied was very simple - point-mass gravitational effects due to the Earth combined with an impulsive delta-V in the velocity direction for the maneuver - generalizations to more complex scenarios, including high fidelity force models, finite duration maneuvers, and maneuver pointing errors, are straightforward and are discussed in the conclusion.

Schiff, Conrad↗
NASA NTRS1992

Transient thermal stress recovery for structural models

A method for computing transient thermal stress vectors from temperature vectors is described. The three step procedure involves the use of NASTRAN to generate an influence coefficient matrix which relates temperatures to stresses in the structural model. The transient thermal stresses are then recovered and sorted for maximum and minimum values. Verification data for the procedure is also provided.

Walls, William↗
NASA NTRS2008

Effect of Adding a Regenerator to Kornhauser's MIT "Two-Space" (Gas-Spring+Heat Exchanger) Test Rig

This study employed entropy-based second law post-processing analysis to characterize the various thermodynamic losses inside a 3-space solution domain (gas spring+heat exchanger+regenerator) operating under conditions of oscillating pressure and oscillating flow. The 3- space solution domain is adapted from the 2-space solution domain (gas spring+heat exchanger) in Kornhauser's MIT test rig by modifying the heat exchanger space to include a porous regenerator system. A thermal nonequilibrium model which assumes that the regenerator porous matrix and gas average temperatures can differ by several degrees at a given axial location and time during the cycle is employed. An important and primary objective of this study is the development and application of a thermodynamic loss post-processor to characterize the major thermodynamic losses inside the 3-space model. It is anticipated that the experience gained from thermodynamic loss analysis of the simple 3-space model can be extrapolated to more complex systems like the Stirling engine. It is hoped that successful development of loss post-processors will facilitate the improvement of the optimization capability of Stirling engine analysis codes through better understanding of the heat transfer and power losses. It is also anticipated that the incorporation of a successful thermal nonequilibrium model of the regenerator in Stirling engine CFD analysis codes, will improve our ability to accurately model Stirling regenerators relative to current multidimensional thermal-equilibrium porous media models.

Ebiana, Asuquo B.↗
DOE OSTI2022

Machine Learning Moment Closure Models for the Radiative Transfer Equation III: Enforcing Hyperbolicity and Physical Characteristic Speeds

This is the third paper in a series in which we develop machine learning (ML) moment closure models for the radiative transfer equation. In our previous work (Huang et al. in J Comput Phys 453:110941, 2022), we proposed an approach to learn the gradient of the unclosed high order moment, which performs much better than learning the moment itself and the conventional $P_N$ closure. However, while the ML moment closure has better accuracy, it is not able to guarantee hyperbolicity and has issues with long time stability. In our second paper (Huang et al., in: Machine learning moment closure models for the radiative transfer equation II: enforcing global hyperbolicity in gradient based closures, 2021. arXiv:2105.14410), we identified a symmetrizer which leads to conditions that enforce that the gradient based ML closure is symmetrizable hyperbolic and stable over long time. The limitation of this approach is that in practice the highest moment can only be related to four, or fewer, lower moments. In this paper, we propose a new method to enforce the hyperbolicity of the ML closure model. Motivated by the observation that the coefficient matrix of the closure system is a lower Hessenberg matrix, we relate its eigenvalues to the roots of an associated polynomial. Here, we design two new neural network architectures based on this relation. The ML closure model resulting from the first neural network is weakly hyperbolic and guarantees the physical characteristic speeds, i.e., the eigenvalues are bounded by the speed of light. The second model is strictly hyperbolic and does not guarantee the boundedness of the eigenvalues. Several benchmark tests including the Gaussian source problem and the two-material problem show the good accuracy, stability and generalizability of our hyperbolic ML closure model.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Identification and Quantification of Carbon Phases in Conversion Fuel for the Transient Reactor Test Facility

As part of an overall effort to convert US research reactors to low-enriched uranium (LEU) fuel use, a LEU conversion fuel is being designed for the Transient Reactor Test Facility (TREAT) at the Idaho National Laboratory. TREAT fuel compacts are comprised of UO2 fuel particles in a graphitic matrix material. In order to refine heat transfer modeling, as well as determine other physical and nuclear characteristics of the fuel, the amount and type of graphite and non-graphite phases within the fuel matrix must be known. In this study, we performed a series of complementary analyses, designed to allow detailed characterization of the graphite and phenolic resin based fuel matrix. Methods included Scanning Electron and Transmission Electron Microscopies, Raman spectroscopy, X-ray Diffraction, and Dual-Beam Focused Ion Beam Tomography. Our results indicate that no single characterization technique will yield all of the desired information; however, through the use of statistical and empirical data analysis, such as curve fitting, partial least squares regression, volume extrapolation and spectra peak ratios, a degree of certainty for the quantity of each phase can be obtained.

10 SYNTHETIC FUELS↗
DOE OSTI2016

Identification and Quantification of Carbon Phases in Conversion Fuel for the Transient Reactor Test Facility

As part of an overall effort to convert US research reactors to low-enriched uranium (LEU) fuel, a LEU conversion fuel is being designed for the Transient Reactor Test Facility (TREAT) at the Idaho National Laboratory. TREAT fuel compacts are comprised of UO{sub 2} fuel particles in a graphitic matrix material. In order to refine heat transfer modeling, as well as determine other physical and nuclear characteristics of the fuel, the amount and type of graphite and non-graphite phases within the fuel matrix must be known. In this study, we performed a series of complementary analyses, designed to allow detailed characterization of the graphite and phenolic resin based fuel matrix. Methods included Scanning Electron Microscopy, Raman spectroscopy, and X-ray Diffraction. Our results indicate that no single characterization technique will yield all of the desired information; however, through the use of statistical and empirical data analysis, such as curve fitting, partial least squares regression, volume extrapolation and spectra peak ratios, a degree of certainty for the quantity of each phase can potentially be obtained. Although no one type of analysis has provided a complete determination of the quantities and identities of all carbon phases present in the fuel matrix, additional information has been gathered and is being used to refine future analysis methods. From the three types of quantification performed, results have been within 1% (92.8%, 93.02% and 93.6% for BSE image analysis, XRD and Raman spectroscopy, respectively). More work remains to identify any trace crystalline carbonaceous phases that may be present, quantify any traces that are identified, and develop and use further complementary bulk phase quantification techniques. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1997

State-Space System Realization with Input- and Output-Data Correlation

This paper introduces a general version of the information matrix consisting of the autocorrelation and cross-correlation matrices of the shifted input and output data. Based on the concept of data correlation, a new system realization algorithm is developed to create a model directly from input and output data. The algorithm starts by computing a special type of correlation matrix derived from the information matrix. The special correlation matrix provides information on the system-observability matrix and the state-vector correlation. A system model is then developed from the observability matrix in conjunction with other algebraic manipulations. This approach leads to several different algorithms for computing system matrices for use in representing the system model. The relationship of the new algorithms with other realization algorithms in the time and frequency domains is established with matrix factorization of the information matrix. Several examples are given to illustrate the validity and usefulness of these new algorithms.

Juang, Jer-Nan↗
NASA NTRS2005

Mechanical Properties of Nanostructured Materials Determined Through Molecular Modeling Techniques

The potential for gains in material properties over conventional materials has motivated an effort to develop novel nanostructured materials for aerospace applications. These novel materials typically consist of a polymer matrix reinforced with particles on the nanometer length scale. In this study, molecular modeling is used to construct fully atomistic models of a carbon nanotube embedded in an epoxy polymer matrix. Functionalization of the nanotube which consists of the introduction of direct chemical bonding between the polymer matrix and the nanotube, hence providing a load transfer mechanism, is systematically varied. The relative effectiveness of functionalization in a nanostructured material may depend on a variety of factors related to the details of the chemical bonding and the polymer structure at the nanotube-polymer interface. The objective of this modeling is to determine what influence the details of functionalization of the carbon nanotube with the polymer matrix has on the resulting mechanical properties. By considering a range of degree of functionalization, the structure-property relationships of these materials is examined and mechanical properties of these models are calculated using standard techniques.

Clancy, Thomas C.↗
DOE OSTI2018

Universality for 1d Random Band Matrices: Sigma-Model Approximation

The paper continues the development of the rigorous supersymmetric transfer matrix approach to the random band matrices started in (J Stat Phys 164:1233–1260, 2016; Commun Math Phys 351:1009–1044, 2017). We consider random Hermitian block band matrices consisting of W×W random Gaussian blocks (parametrized by j,k∈Λ=[1,n]{sup d}∩ℤ{sup d}) with a fixed entry’s variance J{sub jk}=δ{sub j,k}W{sup −1}+βΔ{sub j,k}W{sup −2}, β>0 in each block. Taking the limit W→∞ with fixed n and β, we derive the sigma-model approximation of the second correlation function similar to Efetov’s one. Then, considering the limit β,n→∞, we prove that in the dimension d=1 the behaviour of the sigma-model approximation in the bulk of the spectrum, as β≫n, is determined by the classical Wigner–Dyson statistics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Investigation of the discontinuous precipitation of U-Nb alloys via thermodynamic analysis and phase-field modeling

U-Nb’s discontinuous precipitation, $\gamma^{bcc}_{matrix}$, $\alpha^{orth}_{cellular}$, $\gamma\prime^{bcc}_{cellular}$, is intriguing in the sense that it allows formation and growth of the metastable $\gamma\prime$ phase during the course of its occurrence. Previous attempts to explain the thermodynamic origin of U-Nb’s discontinuous precipitation hypothesized that the energy of $\alpha$ forms an intermediate common tangent with the first potential of the double-well energy of $\gamma$ at the $\gamma\prime$ composition. In this work, we examine different possible mechanisms by which the discontinuous precipitation product in the U-Nb system can be stabilized. We put forward a mechanism by which the bulk free energy of the $\gamma$ can develop a non-monotonic dependency with composition. Additionally we posit that local contributions due to lattice mismatch between the $\alpha$ and $\gamma$ phases may be responsible for the generation of metastable states that may stabilized by thermodynamics as well as by kinetics. Our work suggests that local misfit strain tends to play a crucial role in the growth of the discontinuous precipitation product. Depending on the magnitude of strain developed at the $\alpha$/$\gamma\prime$ interfaces, either an increasing $\gamma\prime$ composition or a random distribution of $\gamma\prime$ compositions around the equiatomic value with respect to increasing temperature could be expected. Moreover, we show how it is possible to stabilize the discontinuous precipitation front through highly anisotropic and fast interface diffusion.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Water Remediation: PVA-Based Magnetic Gels as Efficient Devices to Heavy Metal Removal

Scientific and technological researches are devoted to obtain materials capable of retaining different kinds of pollutants, contributing to contamination solutions. In this context, hydrogels have emerged as great candidates because of their excellent absorption properties as well as good mechanical, thermal and chemical properties. More specifically, ferrogels (magnetic gels) present the extra advantage of being easily manipulated by a permanent magnet. Here, we present the results derived from the application of ferrogels as efficient tools to extract heavy metal pollutants from wastewater samples. The gels were prepared following the method of freezing and thawing of a polyvinyl alcohol aqueous solution with magnetic nanoparticles coated with polyacrylic acid. Ferrogels were fully characterized and their ability to retain Cu{sup 2+} and Cd{sup 2+}, as model heavy metals, was studied. Thus kinetics and mechanisms of adsorption were evaluated and modeled. The concentration of MNPs on the PVA matrix was key to improve the adsorption capability (approximately the double of retention is improved by the MNPs addition). The adsorption kinetics was determined as pseudo-second order model, whereas the Langmuir model was the most appropriate to explain the behavior of the gels. Finally reuse ability was evaluated to determine the real potential of these materials, the ferrogels demonstrated high efficiency up to about five cycles, retaining about 80–90% of their initial adsorption capability. All the results indicated that the materials are promising candidates able to compete with the commercial technology regarding to water remediation.

36 MATERIALS SCIENCE↗
NASA NTRS2003

Comparison of Flow-Dependent and Static Error Correlation Models in the DAO Ozone Data Assimilation System

In a data assimilation system the forecast error covariance matrix governs the way in which the data information is spread throughout the model grid. Implementation of a correct method of assigning covariances is expected to have an impact on the analysis results. The simplest models assume that correlations are constant in time and isotropic or nearly isotropic. In such models the analysis depends on the dynamics only through assumed error standard deviations. In applications to atmospheric tracer data assimilation this may lead to inaccuracies, especially in regions with strong wind shears or high gradient of potential vorticity, as well as in areas where no data are available. In order to overcome this problem we have developed a flow-dependent covariance model that is based on short term evolution of error correlations. The presentation compares performance of a static and a flow-dependent model applied to a global three- dimensional ozone data assimilation system developed at NASA s Data Assimilation Office. We will present some results of validation against WMO balloon-borne sondes and the Polar Ozone and Aerosol Measurement (POAM) III instrument. Experiments show that allowing forecast error correlations to evolve with the flow results in positive impact on assimilated ozone within the regions where data were not assimilated, particularly at high latitudes in both hemispheres and in the troposphere. We will also discuss statistical characteristics of both models; in particular we will argue that including evolution of error correlations leads to stronger internal consistency of a data assimilation ,

Wargan, K.↗
NASA NTRS1989

The evolution and discharge of electric fields within a thunderstorm

An analysis of the present three-dimensional thunderstorm electrical model and its finite-difference approximations indicates unconditional stability for the discretization that results from the approximation of the spatial derivatives by a box-schemelike method and of the temporal derivative by either a backward-difference or Crank-Nicholson scheme. Lightning propagation is treated through numerical techniques based on the inverse-matrix modification formula and Cholesky updates. The model is applied to a storm observed at the Kennedy Space Center in 1978, and numerical comparisons are conducted between the model and the theoretical results obtained by Wilson (1920) and Holzer and Saxon (1952).

Hager, William W.↗
NASA NTRS1992

Random-anisotropy Blume-Emery-Griffiths model

The results are described of studies of a random-anisotropy Blume-Emery-Griffiths spin-1 Ising model using mean-field theory, transfer-matrix calculations, and position-space renormalization-group calculations. The interplay between the quenched randomness of the anisotropy and the annealed disorder introduced by the spin-1 model leads to a rich phase diagram with a variety of phase transitions and reentrant behavior. The results may be relevant to the study of the phase separation of He-3 - He-4 mixtures in porous media in the vicinity of the superfluid transition.

Maritan, Amos↗
DOE OSTI2019

Axial Tensile Failure Analysis of SiC{sub f}/Ti Composite Based on Continuum Cohesive Zone Model

In this paper, the failure behavior of SiC fiber-reinforced Ti-6Al-4V matrix composite (SiC{sub f}/Ti) under longitudinal loads was studied experimentally and theoretically. Experimental results show that the interface of SiC{sub f}/Ti composite is continuous and integral. Through SEM observations, it was found that transverse matrix cracks initiated from the debonding region of the interface. A three-dimensional representative volume element was developed to simulate the failure process of SiC{sub f}/Ti-6Al-4V composite. A continuum cohesive zone model is employed to describe the debonding behavior at the interface, compared with a discrete one. More accurate simulated results were obtained by a continuum cohesive zone model, corresponding to experimental results. The transverse cracking behavior of the matrix was simulated. The stress–strain curves predicted by the continuum cohesive zone model show a good agreement with the experimental results. The fiber volume fraction, interface strength and interface toughness are investigated to exhibit their influences upon the longitudinal mechanical characteristics of SiC{sub f}/Ti-6Al-4V composite.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Transport Upscaling under Flow Heterogeneity and Matrix-Diffusion in Three-Dimensional Discrete Fracture Networks

For this work, we investigate the combined effects of network scale flow variability and retention due to matrix-diffusion on the scaling behavior of transport through fractured media. Two of the principal mechanisms controlling the transport of solutes through fractured low-permeability media are broad distributions of flow velocities and retention times in the solid matrix. We study the relative impact of these two processes under different initial conditions using a set of three-dimensional discrete fracture network simulations. We use these simulations to develop and calibrate an upscaled continuous time random walk (CTRW) approach for advective transport based on an Ornstein-Uhlenbeck model for the particle velocities that accounts for the fracture-matrix coupling using a compound Poisson process. This CTRW model can be conditioned on the initial solute distribution and allows to observe late-time scaling behavior at distances beyond what is feasible to observe using high-fidelity direct numerical simulations. We determine that the initial distribution of particles leads to marked differences in the persistent long-term scale behavior in the solute travel time distributions, even those undergoing retention due to matrix diffusion through implementation and analysis of the model.

54 ENVIRONMENTAL SCIENCES↗