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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Multi-channel 1 to 2 Matrix Elements in Finite Volume

We derive a model-independent expression for finite-volume matrix elements. Specifically, we present a relativistic, non-perturbative analysis of the matrix element of an external current between a one-scalar in-state and a two-scalar out-state. Our result, which is valid for energies below higher-particle inelastic thresholds, generalizes the Lellouch-Luscher formula in two ways: we allow the external current to inject arbitrary momentum into the system and we allow for the final state to be composed an arbitrary number of strongly coupled two-particle states with arbitrary partial waves (including partial-wave mixing induced by the volume). We also illustrate how our general result can be applied to some key examples, such as heavy meson decays and meson photo production.

Walker-Loud, Andre↗
DOE OSTI2018

Analytical Methods of Approximating the Fisher Information Matrix for the Lognormal Distribution

Modelling the matter density in the universe through time has become one of the most important open problems in cosmology. One distribution that works particularly well is the lognormal distribution, where if Y follows a lognormal distribution, ln Y follows a normal distribution. In many cases, we want to know how much information a variable carries about some unknown parameter; this can be accomplished by calculating the Fisher information matrix. However, for the lognormal distribution, the Fisher matrix is much more complicated to evaluate in closed form. In this paper, we will look at some analytical methods to evaluate the Fisher matrix.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Singlet fission in linear chains of molecules

We develop a model configuration interaction Hamiltonian to study the electronic structure of a chain of molecules undergoing singlet fission. We first consider models for dimer and trimer and then we use a matrix partitioning technique to build models of arbitrary size able to describe the relevant electronic structure for singlet fission in linear aggregates. We find that the multi-excitonic state (ME) is stabilized at short inter-monomer distance and the extent of this stabilization depends upon the size of orbital coupling between neighboring monomers. We also find that the coupling between ME states located on different molecules is extremely small leading to bandwidths in the order of ∼10 meV. This observation suggests that multi-exciton states are extremely localized by electron-phonon coupling and that singlet fission involves the transition between a relatively delocalized Frenkel exciton and a strongly localized multi-exciton state. We adopt the methodology commonly used to study non-radiative transitions to describe the singlet fission dynamics in these aggregates and we discuss the limit of validity of the approach. The results indicate that the phenomenology of singlet fission in molecular crystals is different in many important ways from what is observed in isolated dimers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Magnetic states of the quasi-one-dimensional iron chalcogenide Ba 2 FeS 3

Quasi-one-dimensional iron-based ladders and chains, with the 3d iron electronic density n=6, are attracting considerable attention. Recently, a new iron chain system Ba 2 FeS 3 , also with n=6, was prepared under high-pressure and high-temperature conditions. Here the magnetic and electronic phase diagrams are theoretically studied for this quasi-one-dimensional compound. Based on first-principles calculations, a strongly anisotropic one-dimensional electronic band behavior near the Fermi level was observed. In addition, a three-orbital electronic Hubbard model for this chain was constructed. Introducing the Hubbard and Hund couplings and studying the model via the density matrix renormalization group (DMRG) method, we studied the ground-state phase diagram. A robust staggered ↑-↓-↑-↓ AFM region was unveiled in the chain direction, consistent with our density functional theory (DFT) calculations. Furthermore, at intermediate Hubbard U coupling strengths, this system was found to display an orbital selective Mott phase (OSMP) with one localized orbital and two itinerant metallic orbitals. At very large U/W (W=bandwidth), the system displays Mott insulator characteristics, with two orbitals half-filled and one doubly occupied. Our results for high pressure Ba 2 FeS 3 provide guidance to experimentalists and theorists working on this one-dimensional iron chalcogenide chain material.

1-dimensional systems↗
DOE OSTI2016

Large covariance matrices: smooth models from the two-point correlation function

We introduce a new method for estimating the covariance matrix for the galaxy correlation function in surveys of large-scale structure. Our method combines simple theoretical results with a realistic characterization of the survey to dramatically reduce noise in the covariance matrix. For example, with an investment of only ≈1000 CPU hours we can produce a model covariance matrix with noise levels that would otherwise require ~35 000 mocks. Non-Gaussian contributions to the model are calibrated against mock catalogues, after which the model covariance is found to be in impressive agreement with the mock covariance matrix. Since calibration of this method requires fewer mocks than brute force approaches, we believe that it could dramatically reduce the number of mocks required to analyse future surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Rolling Element Bearing Stiffness Matrix Determination (Presentation)

Current theoretical bearing models differ in their stiffness estimates because of different model assumptions. In this study, a finite element/contact mechanics model is developed for rolling element bearings with the focus of obtaining accurate bearing stiffness for a wide range of bearing types and parameters. A combined surface integral and finite element method is used to solve for the contact mechanics between the rolling elements and races. This model captures the time-dependent characteristics of the bearing contact due to the orbital motion of the rolling elements. A numerical method is developed to determine the full bearing stiffness matrix corresponding to two radial, one axial, and two angular coordinates; the rotation about the shaft axis is free by design. This proposed stiffness determination method is validated against experiments in the literature and compared to existing analytical models and widely used advanced computational methods. The fully-populated stiffness matrix demonstrates the coupling between bearing radial, axial, and tilting bearing deflections.

17 WIND ENERGY↗
DOE OSTI2021

Tracing the Mott-Hubbard transition in one-dimensional Hubbard models without Umklapp scattering

Here we apply the density-matrix renormalization group (DMRG) method to a one-dimensional Hubbard model that lacks Umklapp scattering and thus provides an ideal case to study the Mott-Hubbard transition analytically and numerically. The model has a linear dispersion and displays a metal-to-insulator transition when the Hubbard interaction U equals the bandwidth U c = W, where the single-particle gap opens linearly, Δ(U ≥ W) = U - W. The simple nature of the elementary excitations permits us to determine numerically with high accuracy the critical interaction strength and the gap function in the thermodynamic limit. The jump discontinuity of the momentum distribution $n_k$ at the Fermi wave number $k_F$ = 0 cannot be used to locate accurately U c from finite-size systems. However, the slope of $n_k$ at the band edges $k_B$ = ±π, reveals the formation of a single-particle bound state which can be used to determine U c reliably from $n_k$ using accurate finite-size data.

1-dimensional systems↗
DOE OSTI2022

A New Model for Simulating the Imbibition of a Wetting-Phase Fluid in a Matrix-Fracture Dual Connectivity System

The imbibition experiment is an effective approach for measuring petrophysical properties of porous media, with many such experiments performed over the past decade. Quite some empirical, analytical, and numerical models have been developed to simulate spontaneous imbibition of the wetting phase fluid into porous media, but limitations still exist. In previous studies, the imbibition process has been considered to give a piston-like displacement or the porous medium modeled as multiply-sized pores linked with bonds; both approaches fail to yield comprehensive results due to their neglect of the presence of irregular fractures or nonuniform flow paths through the matrix. By building a numerical model for simulating laboratory-scale experimental data, we performed imbibition tests on several fractured Barnett Shale samples having fractures either parallel ( P ) or transverse ( T ) to the bedding plane and used MATLAB to build a new numerical model by combining the imbibition process in fractures and the matrix using concepts from percolation theory. The experimental data show that the rocks with P -direction fractures have a more steady increase of imbibition rates than the case of T -direction one. As the shale matrix with low pore connectivity hampers the upward water movement, the imbibition rate of shales with T -direction fractures will decrease suddenly after the bottom layer in contact with water is saturated during the initial period. This wetting phase movement (WPM) model can simulate 3D porous media with 2D fractures. The rate of imbibition by fractured porous media is associated with physical parameters such as porosity and fracture distribution (e.g., the number and angle of fractures). Using Monte Carlo methods, we examined fracture parameters and predicted elapsed time and cumulative water imbibition, for the Barnett Shale samples. The results show that the rate of imbibed water mass is sensitive to the number of fractures directly connected to water source, and the connectivity between two neighboring grid cells is a key parameter for the wetting-front progression. The findings of this study can help to better understand the imbibition process with multiple influencing processes and factors in fractured-matrix rocks. Although the experiments, data simulation, and prediction results are based only on Barnett Shale samples, the model is readily applicable to imbibition tests of other fractured rocks to show the spatial and temporal behavior during a dynamic imbibition process that are not easily captured experimentally.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2013

Measurement of baryon acoustic oscillations in the Lyman-α forest fluctuations in BOSS data release 9

We use the Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 9 (DR9) to detect and measure the position of the Baryonic Acoustic Oscillation (BAO) feature in the three-dimensional correlation function in the Lyman-α flux fluctuations at a redshift z{sub eff} = 2.4. The feature is clearly detected at significance between 3 and 5 sigma (depending on the broadband model and method of error covariance matrix estimation) and is consistent with predictions of the standard ΛCDM model. We assess the biases in our method, stability of the error covariance matrix and possible systematic effects. We fit the resulting correlation function with several models that decouple the broadband and acoustic scale information. For an isotropic dilation factor, we measure 100 × (α{sub iso} − 1) = −1.6{sup +2.0+4.3+7.4}{sub −2.0−4.1−6.8} (stat.) ±1.0 (syst.) (multiple statistical errors denote 1,2 and 3 sigma confidence limits) with respect to the acoustic scale in the fiducial cosmological model (flat ΛCDM with Ω{sub m} = 0.27, h = 0.7). When fitting separately for the radial and transversal dilation factors we find marginalised constraints 100 × (α{sub ||} − 1) = −1.3{sup +3.5+7.6+12.3}{sub −3.3−6.7−10.2} (stat.) ±2.0 (syst.) and 100 × (α{sub p}erpendicular − 1) = −2.2{sup +7.4+17}{sub −7.1−15} (stat.) ±3.0 (syst.). The dilation factor measurements are significantly correlated with cross-correlation coefficient of ∼ −0.55. Errors become significantly non-Gaussian for deviations over 3 standard deviations from best fit value. Because of the data cuts and analysis method, these measurements give tighter constraints than a previous BAO analysis of the BOSS DR9 Lyman-α sample, providing an important consistency test of the standard cosmological model in a new redshift regime.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Evaluation of Lower East Fork Poplar Creek Mercury Sources

This report summarizes a 3-year research project undertaken to better understand the nature and magnitude of mercury (Hg) fluxes in East Fork Poplar Creek (EFPC). This project addresses the requirements of Action Plan 1 in the 2011 Oak Ridge Reservation-wide Comprehensive Environmental Response, Compensation, and Liability Act Five Year Review (FYR). The Action Plan is designed to address a twofold 2011 FYR issue: (1) new information suggests mobilization of mercury from the upper and lower EFPC streambeds and stream banks is the primary source of mercury export during high-flow conditions, and (2) the current Record of Decision did not address the entire hydrologic system and creek bank or creek bed sediments. To obtain a more robust watershed-scale understanding of mercury sources and processes in lower EFPC (LEFPC), new field and laboratory studies were coupled with existing data from multiple US Department of Energy programs to develop a dynamic watershed and bioaccumulation model. LEFPC field studies for the project focused primarily on quantification of streambank erosion and an evaluation of mercury dynamics in shallow groundwater adjacent to LEFPC and potential connection to the surface water. The approach to the stream bank study was innovative in using imagery from kayak floats’ surveys from the headwaters to the mouth of EFPC to estimate erosion, coupled with detailed bank soil mercury analyses. The goal of new field assessments and modeling was to generate a more holistic and quantitative understanding of the watershed and the sources, flux, concentration, transformation, and bioaccumulation of inorganic mercury (IHg) and methylmercury (MeHg). Model development used a hybrid approach that dynamically linked a spreadsheet-based physical and chemical watershed model to a systems dynamics, mercury bioaccumulation model for key fish species. The watershed model tracks total Hg and MeHg fluxes and concentrations by examining upstream inputs, floodplain runoff, floodplain leaching, bank soil erosion, and periphyton matrix dynamics. The bioaccumulation model tracks the feeding, growth, and mercury assimilation of representative individual fish through their typical life span using key inputs of fish size, water temperature, and diet. The LEFPC watershed was divided into five modeling reaches, and fluxes and concentrations are assessed at this spatial scale. Following are the key findings of the field and laboratory studies and the watershed and bioaccumulation modeling: • The greatest flux of total mercury (HgT) in LEFPC is related to stormflow transport of Hg-contaminated solids entering the creek because of bank erosion in the upper reaches of the creek. • The second greatest flux originates from upper EFPC (Station 17 representing the exit stream sampling point near the boundary of the Y-12 Complex), and appears to control base flow fluxes. • The observed increase in MeHg concentration and flux from upstream to downstream is related primarily to instream methylation by periphyton and other biological activity. • A meaningful substantial reduction of the HgT flux in LEFPC would require addressing the flux of HgT originating from bank erosion and from Station 17. • Actions to reduce LEFPC floodplain leaching and runoff would not produce much of an impact on HgT or MeHg concentrations or fluxes unless other major sources are eliminated first. This project addresses the Action Plan goal to evaluate the role of LEFPC bank soil sources and to consider the entire EFPC hydrologic system. Model conclusions are dependent on the data available at the time of this assessment. However, a robust understanding and quantification for some mercury-related parameters and relationships is still lacking; there is a continued need for field data collection and modeling improvements. Model predictions should be viewed cautiously, with comparisons of the magnitude of predictions between scenarios being more valid than absolute predictions of concentrations or fluxes. With continued updates and refinement, the watershed-scale model can be a useful, valuable tool for future EFPC research prioritization, technology development, and remedial decision-making.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Optimal laser wavelength for efficient laser power converter operation over temperature

A temperature dependent modeling study is conducted on a GaAs laser power converter to identify the optimal incident laser wavelength for optical power transmission. Furthermore, the respective temperature dependent maximal conversion efficiencies in the radiative limit as well as in a practically achievable limit are presented. The model is based on the transfer matrix method coupled to a two-diode model, and is calibrated to experimental data of a GaAs photovoltaic device over laser irradiance and temperature. Since the laser wavelength does not strongly influence the open circuit voltage of the laser power converter, the optimal laser wavelength is determined to be in the range where the external quantum efficiency is maximal, but weighted by the photon flux of the laser.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Mass spectrometry imaging of natural carbonyl products directly from agar-based microbial interactions using 4-APEBA derivatization

Aliphatic carboxylic acids, aldehydes, and ketones play diverse roles in microbial adaptation to their microenvironment, from excretion as toxins to adaptive metabolites for membrane fluidity. However, the spatial distribution of these molecules throughout biofilms and how microbes in these environments exchange these molecules remain elusive for many of these bioactive species due to inefficient molecular imaging strategies. Herein, we apply on-tissue chemical derivatization (OTCD) using 4-(2-((4-bromophenethyl)dimethylammonio)ethoxy)benzenaminium dibromide (4-APEBA) on a co-culture of a soil bacterium (Bacillus subtilis NCIB 3610) and fungus (Fusarium sp. DS 682) grown on agar as our model system. Using matrix-assisted laser desorption/ionization mass spectrometry imaging (MALDI-MSI), we spatially resolved more than 300 different metabolites containing carbonyl groups within this model system. Various spatial patterns are observable in these species, which indicate possible extracellular or intercellular processes of the metabolites and their up- or downregulation during microbial interaction. The unique chemistry of our approach allowed us to bring additional confidence in accurate carbonyl identification, especially when multiple isomeric candidates were possible, and this provided the ability to generate hypotheses about the potential role of some aliphatic carbonyls in this B. subtilis/Fusarium sp. interaction. The results shown here demonstrate the utility of 4-ABEBA-based OTCD MALDI-MSI in probing interkingdom interactions directly from microbial co-cultures, and these methods will enable future microbial interaction studies with expanded metabolic coverage.

4-APEBA↗
DOE OSTI2023

Quantification of the Impact of Acidified Brine on Fracture-Matrix Transport in a Naturally Fractured Shale Using in Situ Imaging and Modeling

Understanding flow, transport, chemical reactions, and hydromechanical processes in fractured geologic materials is key for optimizing a range of subsurface processes including carbon dioxide and hydrogen storage, unconventional energy resource extraction, and geothermal energy recovery. Flow and transport processes in naturally fractured shale rocks have been challenging to characterize due to experimental complexity and the multiscale nature of quantifying continuum scale descriptions of mass exchange between micrometer-scale fractures and nanometer-scale pores. In this study, we use positron emission tomography (PET) to image the transport of a conservative tracer in a naturally fractured Wolfcamp shale core before and after the core was exposed to low pH brine conditions. Image-based experimental observations are interpreted by fitting an analytical transport model to fracture-containing voxels in the core. Results of this analysis indicate subtle increases in matrix diffusivity and a slightly more uniform fracture velocity distribution following exposure to low pH conditions. These observations are compared with a multicomponent one-dimensional reactive transport model that indicates the capacity for a 10% increase in porosity at the fracture-matrix interface as a result of the low pH brine exposure. This porosity change is the result of the dissolution of carbonate minerals in the shale matrix to low pH conditions. Further, this image-based workflow represents a new approach for quantifying spatially resolved fracture-matrix transport processes and provides a foundation for future work to better understand the role of coupled transport, reaction, and mechanical processes in naturally fractured rocks.

58 GEOSCIENCES↗
DOE OSTI2017

Intrinsic properties and strengthening mechanism of monocrystalline Ni-containing ternary concentrated solid solutions

Ternary single-phase concentrated solid solution alloys (SP-CSAs), so-called "medium entropy alloys", not only possess notable mechanical and physical properties but also form a model system linking the relatively simple binary alloys to the complex high entropy alloys. Our knowledge of their intrinsic properties is vital to understand the material behavior and to prompt future applications. To this end, three model alloys NiCoFe, NiCoCr, and NiFe-20Cr have been selected and grown as single crystals. We measured their elastic constants using an ultrasonic method, and several key materials properties, such as shear modulus, bulk modulus, elastic anisotropy, and Debye temperatures have been derived. Furthermore, nanoindentation tests have been performed on these three alloys together with Ni, NiCo and NiFe on their (100) surface, to investigate the strengthening mechanisms. NiCoCr has the highest hardness, NiFe, NiCoFe and NiFe-20Cr share a similar hardness that is apparently lower than NiCoCr; NiCo has the lowest hardness in the alloys, which is similar to elemental Ni. The Labusch-type solid solution model has been applied to interpret the nanoindentation data, with two approaches used to calculate the lattice mismatch. Finally, by adopting an interatomic spacing matrix method, the Labusch model can reasonably predict the hardening effects for the whole set of materials.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Robust Statistical Approach for Determination of Graphite Nitridation Using Bayesian Model Comparison

A better estimation of surface reaction efficiency of semiconductor-grade graphite with atomic nitrogen, as well as the calibration error are calculated using Bayesian updating based on experimental data. Compared with a conventional deterministic model, the stochastic model approach is a powerful tool in the sense that the model is capable of taking into account underlying error correlations among the data quantities. In this paper, we investigate four different stochastic models (called “stochastic system model classes” herein) corresponding to different descriptions of modeling and measurement error structures, given one deterministic physical model. These stochastic system model classes differ in the covariance matrix structure that is used in the uncertainty model to represent uncertainties associated with the physical model and experimental measurements. For each model class, Bayesian inference is used to estimate the posterior probabilities of the physical model parameters as well as of the stochastic model parameters. Model comparison and selection are then applied based on two measures including Bayesian evidence and Bayesian information criterion, as well as the deviance information criterion. Both measures suggest the stochastic model class, which considers that a correlation between errors in two data quantities among different data points is the most plausible. With the stochastic model class, the range of uncertainty in surface reaction efficiency is estimated to be about two orders of magnitude at [Formula: see text].

Engineering↗
DOE OSTI2020

Modeling Aerial Gamma-Ray Backgrounds Using Non-negative Matrix Factorization

Airborne gamma-ray surveys are useful for many applications, ranging from geology and mining to public health and nuclear security. In all these contexts, the ability to decompose a measured spectrum into a linear combination of background source terms can provide useful insights into the data and lead to improvements in the techniques that use spectral energy windows. Multiple methods for the linear decomposition of spectra exist but are subject to various drawbacks, such as allowing negative photon fluxes or requiring detailed Monte Carlo modeling. In this work, we propose using non-negative matrix factorization (NMF) as a data-driven approach to spectral decomposition. Using aerial surveys that include flights over water, we demonstrate that the mathematical approach of NMF finds physically relevant structure in the aerial gamma-ray background, namely, that measured spectra can be expressed as the sum of nearby terrestrial emission, distant terrestrial emission, and radon and cosmic emission. Furthermore, these NMF background components are compared with the background components obtained by noise-adjusted singular value decomposition (NASVD), which contain negative photon fluxes and, thus, do not represent the emission spectra in as straightforward a way. Finally, we comment on the potential areas of research that are enabled by NMF decompositions, such as new approaches to spectral anomaly detection and data fusion.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Phase-field modeling of rate-dependent fluid-driven fracture initiation and propagation

The rate-dependent behavior associated with deformation and fracturing of materials, such as natural rocks, poses significant challenges for modeling. In addition to the complications of the viscoelastic response, the speed of fracture propagation reflects micromechanical mechanisms in the fracture process zone(FPZ). In order to represent these complicated behaviors, a thermodynamically consistent, rate-dependent fracture model is required. Based on rigorous thermodynamic principles, we derive a rate-dependent phase-field mechanical model coupled with single-phase fluid flow in both the matrix and the fracture. The model is guaranteed to satisfy energy conservation during fracture propagation. Here, the system of equations is solved using the introduced solution procedure and a novel preconditioner that accounts for the complex fluid-structure interaction. The proposed phase-field model is tested against several benchmark problems on solid-fluid coupling, fluid-driven fracture propagation and rate-dependent viscoelastic deformation. The model serves as a strong basis for investigating rate-dependent fracturing experiments and for making predictions of material behaviors under new conditions.

42 ENGINEERING↗
DOE OSTI2015

Quantitative validation of carbon-fiber laminate low velocity impact simulations

Simulations of low velocity impact with a flat cylindrical indenter upon a carbon fiber fabric reinforced polymer laminate are rigorously validated. Comparison of the impact energy absorption between the model and experiment is used as the validation metric. Additionally, non-destructive evaluation, including ultrasonic scans and three-dimensional computed tomography, provide qualitative validation of the models. The simulations include delamination, matrix cracks and fiber breaks. An orthotropic damage and failure constitutive model, capable of predicting progressive damage and failure, is developed in conjunction and described. An ensemble of simulations incorporating model parameter uncertainties is used to predict a response distribution which is then compared to experimental output using appropriate statistical methods. Lastly, the model form errors are exposed and corrected for use in an additional blind validation analysis. The result is a quantifiable confidence in material characterization and model physics when simulating low velocity impact in structures of interest.

36 MATERIALS SCIENCE↗