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,477 records · Page 82
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Online monitoring of solutions within microfluidic chips: simultaneous Raman and UV-Vis absorption spectroscopies

Microfluidics is an appealing analytical tool in the global effort to close the nuclear fuel cycle. Using a microfluidic chip permits analysis of greatly reduced sample volumes compared to what is necessary for traditional analytical methods. There is a commensurate reduction in disposal volume and cost. Development of spectroscopic methods practicable on the microscale are necessary to take advantage of the microchip configuration. This study uses simultaneously applied UV-vis and micro-Raman spectroscopies adapted to function on the microscale to analyze in situ both the Nd3+ (UV-vis active) and HNO3 (Raman active) concentrations in the same sample. An adjustable translation platform was designed to hold the micro-Raman probe above and perpendicular to the chip face, and the UV-vis probe in the plane of the chip. These complimentary spectral techniques when processed through multivariate PLS (Partial Least Squares) models gave an accurate picture of the widely varying solution concentrations as a function of time for each solution component. Solution matrix effects can drastically alter analyte signatures as measured by both UV-vis absorbance and Raman spectroscopy. PLS methods successfully modeled these spectral changes and accurately measured concentrations of components of interest within the microfluidic chip.

on-line monitoring, microfluidics, Raman, UV-vis↗
DOE OSTI2014

Universal tight binding model for chemical reactions in solution and at surfaces. I. Organic molecules

As is now well established, a first order expansion of the Hohenberg–Kohn total energy density functional about a trial input density, namely, the Harris–Foulkes functional, can be used to rationalize a non self consistent tight binding model. If the expansion is taken to second order then the energy and electron density matrix need to be calculated self consistently and from this functional one can derive a charge self consistent tight binding theory. In this paper we have used this to describe a polarizable ion tight binding model which has the benefit of treating charge transfer in point multipoles. This admits a ready description of ionic polarizability and crystal field splitting. It is necessary in constructing such a model to find a number of parameters that mimic their more exact counterparts in the density functional theory. We describe in detail how this is done using a combination of intuition, exact analytical fitting, and a genetic optimization algorithm. Having obtained model parameters we show that this constitutes a transferable scheme that can be applied rather universally to small and medium sized organic molecules. We have shown that the model gives a good account of static structural and dynamic vibrational properties of a library of molecules, and finally we demonstrate the model's capability by showing a real time simulation of an enolization reaction in aqueous solution. In two subsequent papers, we show that the model is a great deal more general in that it will describe solvents and solid substrates and that therefore we have created a self consistent quantum mechanical scheme that may be applied to simulations in heterogeneous catalysis.

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

Effect of double heterogeneity treatment on neutronics modeling of HTGR unit cell - 18524

The High Temperature Gas-cooled Reactor (HTGR) utilizes unique coated fuel particle design, which is embedded within a graphite matrix and shaped into spherical or hexagonal fuel element. The feature of the inherent double heterogeneity in the HTGR fuel element therefore raises difficulties in modeling and requires special simulation strategies to address its impact on the self-shielding calculation. In this study, neutronics calculations with various levels of heterogeneity processing were performed for both spherical and hexagonal HTGR fuel elements. The results indicate that ignoring either level of double heterogeneity can lead to significant underestimation of self-shielding effect and consequently yield as much as 8% deviation in the prediction of the multiplication factor. Fission and capture rates were compared for the spherical HTGR fuel model and over 30% deviation was found in the resonance energy range if homogenized model used. Further investigation on the neutron spectrum demonstrated similar tendency for both types of fuel elements when applying the simplified self-shielding models, that is, {approx}20% underestimation of neutron flux in 6.5-6.75 eV energy interval was observed due to the incomplete evaluation of the cross section resonance self-shielding effect. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

CSPlib: A performance portable parallel software toolkit for analyzing complex kinetic mechanisms

Computational singular perturbation (CSP) is a method to analyze dynamical systems. It targets the decoupling of fast and slow dynamics using an alternate linear expansion of the right-hand side of the governing equations based on eigenanalysis of the associated Jacobian matrix. This representation facilitates diagnostic analysis, detection and control of stiffness, and the development of simplified models. For this work, we have implemented CSP in a C++ open-source library CSPlib using the Kokkos parallel programming model to address portability across diverse heterogeneous computing platforms, i.e., multi/many-core CPUs and GPUs. We describe the CSPlib implementation and present its computational performance across different computing platforms using several test problems. Specifically, we test the CSPlib performance for a constant pressure ignition reactor model on different architectures, including IBM Power 9, Intel Xeon Skylake, and NVIDIA V100 GPU. The size of the chemical kinetic mechanism is varied in these tests. As expected, the Jacobian matrix evaluation, the eigensolution of the Jacobian matrix, and matrix inversion are the most expensive computational tasks. When considering the higher throughput characteristic of GPUs, GPUs performs better for small matrices with higher occupancy rate. CPUs gain more advantages from the higher performance of well-tuned and optimized linear algebra libraries such as OpenBLAS.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

On the reinforcement homogenization in CNT/metal matrix composites during severe plastic deformation

Highlights: • Carbon nanotube-reinforced composites were processed by severe plastic deformation. • The evolution of the second phase distribution homogeneity is thoroughly analyzed. • The analysis is carried out with increasing equivalent strains. • A combination of two methods is applied and discussed. • A model that predicts the equivalent strain necessary for homogeneity is proposed. - Abstract: Carbon nanotube (CNT)-reinforced nickel matrix composites with different concentrations were processed by high pressure torsion (HPT). We thoroughly characterized the CNT agglomerates' spatial arrangement at different stages of deformation in order to extract information valuable for the optimization of the processing parameters and to elucidate the mechanisms involved during the processing of particle reinforced metal matrix composites by HPT. From the electron micrographs taken on the radial direction with increasing equivalent strains, we observed that CNT agglomerates debond by relative sliding between CNT during HPT, becoming spherical at higher stages of deformation. Furthermore, we introduced a model for the prediction of the minimum strain required for a homogeneous distribution of a second phase during HPT, which can be correlated to the material's three-dimensional structure and agrees well with the experimental data.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Modeling of Magnetic Properties of Magnetorheological Elastomers Using JA Hysteresis Model

Magnetorheological elastomers (MREs) are composite materials that consist of magnetically permeable particles in a nonmagnetic polymeric matrix. Under the influence of an external magnetic field, a reversible deformation change occurs in the mechanical properties of these materials. Due to their coupled magnetomechanical response, these materials have been found suitable for a range of applications including tunable vibration absorbers, sensors, and actuators. Notably, improvement of such devices are prerequisites to efficient energy conversion systems, hence the need to understand further the MRE technology. The Jiles-Atherton (JA) theory takes into consideration the magneto-coupling experienced by effective domains in a magnetic material. Algorithm based on the theory yields five model parameters; saturation magnetization (M s ), domain density (a), domain coupling (α), loss coefficient (k), and reversibility (c). Using JA theory, model parameters were calculated and linked to the physical attributes of Fe powder and isotropic MRE. The results show that the calculated Ms for the MRE is reasonably related to that of the Fe powder by a factor of the particle's volume fraction used in the MRE. Further, the calculated k, a, and α provided support for the reduced pinning factor, domain density, and increased domain coupling in the MRE due to the changes in the domain structure between the two materials. From the calculated JA parameters, finite-element modeling (FEM) of the MRE hysteresis loop was performed. The analysis showed that the modeled magnetic properties including coercivity, remanence, and coordinates of the hysteresis loop tip vary with geometric position.

42 ENGINEERING↗
DOE OSTI2023

Assessing Multi-Dimensional Impacts of Achieving Sustainability Goals by Projecting the Sustainable Agriculture Matrix Into the Future

The concept of sustainability inherently spans multiple spatial scales, sectors, variables, and time horizons. This study links a recently developed method of assessing present-day agricultural sustainability across environmental, economic, and social dimensions with a process-based integrated assessment model, in order to allow forward-looking analysis of sustainability by region and scenario. The sustainable agriculture matrix estimates present-day agricultural sustainability at the national level using 18 indicator variables, of which this study estimates nine to the year 2100, using an enhanced version of the Global Change Analysis Model. Scenarios include a reference scenario, and scenarios that apply the following measures, both individually and in combination, that are thought to improve sustainability: yield intensification, transition toward more plant-based (“flexitarian”) diets, and economy-wide greenhouse gas emissions mitigation. The scenarios illustrate considerable complexity and tradeoffs inherent to efforts to improve agricultural sustainability in all regions globally. For example, yield intensification typically increases nitrogen pollution, flexitarian diets can reduce agricultural output, and greenhouse gas mitigation efforts may either increase deforestation or crowd out crop and livestock production due to consequent bioenergy demands. However, there is considerable inter-regional heterogeneity in the responses, and the importance of such secondary responses also differs by region. The analysis and post-processing methods developed in this study allow quantification and visualization of the absolute and relative magnitude of the tradeoffs between agricultural sustainability indicator variables across regions, time periods, and scenarios.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Analysis and modeling of tungsten emission and net erosion in the DIII-D divertor using updated atomic data

Tungsten (W) is one of the leading candidate materials for plasma-facing components. However, its main drawback is its high radiative efficiency; if W penetrates the plasma, it can lead to core degradation or even collapse. Since eroded tungsten tends to ionize in the sheath and redeposit promptly, the net erosion flux that escapes prompt redeposition can differ significantly from the gross erosion. This work presents a modeling framework to estimate net erosion and photon emission from W coatings exposed to the lower divertor of DIII-D using the DiMES material exposure probe. The approach couples RustBCA for sputtering yields with a Monte Carlo transport code (LPTMC) that models redeposition and W emission. Computation is carried out with new atomic data, based on R-matrix and Mons calculations, leading to lower ionization probabilities and a twofold increase in net erosion estimates compared to calculations done with OPEN-ADAS atomic data. The model results are benchmarked against experimental measurements, showing quantitative agreement for erosion, although the trends in W emission are reproduced only qualitatively. The model is also used to assess whether W II emission can serve as a direct measurement of the net erosion of W in the lower divertor of DIII-D. Simulations show that this is not valid if the electron pressure is above ~120 Pa or if the toroidal length of the eroded material is smaller than the parallel-to-B distance traveled by impurity ions before steady-state conditions are reached. Finally, simulations suggest that when W is sputtered by carbon ions with high impact energies (≳300 eV) in DIII-D, W net erosion scales with W gross erosion and can be numerically approximated using W I flux alone as input.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Magnetic excitations, nonclassicality, and quantum wake spin dynamics in the Hubbard chain

Recent work has demonstrated that quantum Fisher information (QFI), a witness of multipartite entanglement, and magnetic Van Hove correlations G(r,t), a probe of local real-space real-time spin dynamics, can be successfully extracted from inelastic neutron scattering on spin systems through accurate measurements of the dynamical spin structure factor S(k,ω). In this work, we apply theoretically these ideas to the half-filled Hubbard chain with nearest-neighbor hopping, away from the strong-coupling limit. This model has nontrivial redistribution of spectral weight in S(k,ω) going from the noninteracting limit (U = 0) to strong coupling (U→∞), where it reduces to the Heisenberg quantum spin chain. We use the density matrix renormalization group to find S(k,ω), from which QFI is then calculated. We find that QFI grows with U. With realistic energy resolution it becomes capable of witnessing bipartite entanglement above U=2.5 (in units of the hopping), where it also changes slope. This point is also proximate to slope changes of the bandwidth W(U) and the half-chain von Neumann entanglement entropy. We compute G(r,t) by Fourier transforming S(k,ω). The results indicate a crossover in the short-time short-distance dynamics at low U characterized by ferromagnetic light-cone wavefronts, to a Heisenberg-type behavior at large U featuring antiferromagnetic light cones and spatially period-doubled antiferromagnetism. We find this crossover has largely been completed by U=3. Our results thus provide evidence that, in several aspects, the strong-coupling limit of the Hubbard chain is reached qualitatively already at a relatively modest interaction strength. We discuss experimental candidates for observing the G(r,t) dynamics found at low U.

1-dimensional systems↗
DOE OSTI2015

Multireference linearized coupled cluster theory for strongly correlated systems using matrix product states

We propose a multireference linearized coupled cluster theory using matrix product states (MPSs-LCC) which provides remarkably accurate ground-state energies, at a computational cost that has the same scaling as multireference configuration interaction singles and doubles, for a wide variety of electronic Hamiltonians. These range from first-row dimers at equilibrium and stretched geometries to highly multireference systems such as the chromium dimer and lattice models such as periodic two-dimensional 1-band and 3-band Hubbard models. The MPS-LCC theory shows a speed up of several orders of magnitude over the usual Density Matrix Renormalization Group (DMRG) algorithm while delivering energies in excellent agreement with converged DMRG calculations. Also, in all the benchmark calculations presented here, MPS-LCC outperformed the commonly used multi-reference quantum chemistry methods in some cases giving energies in excess of an order of magnitude more accurate. As a size-extensive method that can treat large active spaces, MPS-LCC opens up the use of multireference quantum chemical techniques in strongly correlated ab initio Hamiltonians, including two- and three-dimensional solids.

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

Hydrologic Model Analyses for Decision Support at the Los Alamos National Laboratory Site - 18363

Modeling hydrologic flow and transport in the subsurface presents a number of challenges including large datasets that contain limited information about the properties of the subsurface and models with a large number of parameters characterizing the properties of the subsurface and the ambient hydrologic conditions (e.g., groundwater recharge). This, in turn, leads to substantial uncertainty in predictions about quantities of interest (e.g., the concentration of a contaminant at points of compliance). We present a number of methods and tools that have recently been developed to address these problems. The model calibration tools are capable of assimilating large datasets (for example, on the order of 1,000,000 model calibration targets) into numerical models, calibrating models with a large number of model parameters (on the order of 1,000,000), or both. These tools leverage recently developed methods in numerical linear algebra such as matrix-free methods and randomized matrix sketching. The uncertainty quantification tools are highly-parallel and enable us to robustly quantify the uncertainties in model parameters (and therefore uncertainty in quantities of interest). Having these capabilities in model analysis tools enables more informed decision making related to remediation of groundwater contaminants. The model analyses that we present are related to the Los Alamos National Laboratory (LANL) chromium site. The analyses are performed by MADS (Model Analysis and Decision Support; http://mads.lanl.gov). (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

A Modified Network Alteration Model to Predict Quasi-Static Behavior of the Cross-Linked Polymers During Hydrolytic Aging

Abstract The network alteration theory is utilized to capture the effects of hydrolytic aging on the mechanical integrity of the polymer matrix. Hydrolytic degradation generally considered as a result of two simultaneous chemical phenomena; (i) reduction of the cross-links, and (ii) active chain loss. Assuming the alteration of the chain density is identical to the change of the number average molecular weight along the aging trajectory, the strain energy of polymer matrix can be rewritten as a function of deformation, deformation history, storage time, and aging temperature. The proposed model includes only five physically inspired material parameters. Thus, while it is computationally efficient, it shows good agreement with own experimental data, which performed on various range of accelerated aging temperatures and times. With respect to its computational efficiency, simplicity, accuracy, and interpret-ability, the model is a good choice for advanced implementations in FE programs.

Bahrololoumi, Amir↗
DOE OSTI2014

A 3D Orthotropic Strain-Rate Dependent Elastic Damage Material Model.

A three dimensional orthotropic elastic constitutive model with continuum damage and cohesive based fracture is implemented for a general polymer matrix composite lamina. The formulation assumes the possibility of distributed (continuum) damage followed b y localized damage. The current damage activation functions are simply partially interactive quadratic strain criteria . However, the code structure allows for changes in the functions without extraordinary effort. The material model formulation, implementation, characterization and use cases are presented.

English, Shawn Allen↗
DOE OSTI2014

Comparison of Fission Product Release Predictions Using PARFUME with Results from the AGR-1 Irradiation Experiment

This report documents comparisons between post-irradiation examination measurements and model predictions of silver (Ag), cesium (Cs), and strontium (Sr) release from selected tristructural isotropic (TRISO) fuel particles and compacts during the first irradiation test of the Advanced Gas Reactor program that occurred from December 2006 to November 2009 in the Advanced Test Reactor (ATR) at Idaho National Laboratory (INL). The modeling was performed using the particle fuel model computer code PARFUME (PARticle FUel ModEl) developed at INL. PARFUME is an advanced gas-cooled reactor fuel performance modeling and analysis code (Miller 2009). It has been developed as an integrated mechanistic code that evaluates the thermal, mechanical, and physico-chemical behavior of fuel particles during irradiation to determine the failure probability of a population of fuel particles given the particle-to-particle statistical variations in physical dimensions and material properties that arise from the fuel fabrication process, accounting for all viable mechanisms that can lead to particle failure. The code also determines the diffusion of fission products from the fuel through the particle coating layers, and through the fuel matrix to the coolant boundary. The subsequent release of fission products is calculated at the compact level (release of fission products from the compact) but it can be assessed at the particle level by adjusting the diffusivity in the fuel matrix to very high values. Furthermore, the diffusivity of each layer can be individually set to a high value (typically 10-6 m2/s) to simulate a failed layer with no capability of fission product retention. In this study, the comparison to PIE focused on fission product release and because of the lack of failure in the irradiation, the probability of particle failure was not calculated. During the AGR-1 irradiation campaign, the fuel kernel produced and released fission products, which migrated through the successive layers of the TRISO-coated particle and potentially through the compact matrix. The release of these fission products was measured in PIE and modeled with PARFUME.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2016

Uncertainty analysis in geospatial merit matrix–based hydropower resource assessment

Hydraulic head and mean annual streamflow, two main input parameters in hydropower resource assessment, are not measured at every point along the stream. Translation and interpolation are used to derive these parameters, resulting in uncertainties. This study estimates the uncertainties and their effects on model output parameters: the total potential power and the number of potential locations (stream-reach). These parameters are quantified through Monte Carlo Simulation (MCS) linking with a geospatial merit matrix based hydropower resource assessment (GMM-HRA) Model. The methodology is applied to flat, mild, and steep terrains. Results show that the uncertainty associated with the hydraulic head is within 20% for mild and steep terrains, and the uncertainty associated with streamflow is around 16% for all three terrains. Output uncertainty increases as input uncertainty increases. However, output uncertainty is around 10% to 20% of the input uncertainty, demonstrating the robustness of the GMM-HRA model. Hydraulic head is more sensitive to output parameters in steep terrain than in flat and mild terrains. Furthermore, mean annual streamflow is more sensitive to output parameters in flat terrain.

13 HYDRO ENERGY↗
DOE OSTI2013

Scaling behavior and phase diagram of a PT-symmetric non-Hermitian Bose-Hubbard system

We study scaling behavior and phase diagram of a PT-symmetric non-Hermitian Bose-Hubbard model. In the free interaction case, using both analytical and numerical approaches, the metric operator for many-particle is constructed. The derived properties of the metric operator, similarity matrix and equivalent Hamiltonian reflect the fact that all the matrix elements change dramatically with diverging derivatives near the exceptional point. In the nonzero interaction case, it is found that even small on-site interaction can break the PT symmetry drastically. It is demonstrated that the scaling law can be established for the exceptional point in both small and large interaction limit. Based on perturbation and numerical methods, we also find that the phase diagram shows rich structure: there exist multiple regions of unbroken PT symmetry. - Highlights: Black-Right-Pointing-Pointer Generalization is done to many-particle sector of PT-symmetric non-Hermitian system. Black-Right-Pointing-Pointer The metric operator, similarity matrix and equivalent Hamiltonian possess diverging derivatives near the exceptional point. Black-Right-Pointing-Pointer The scaling law is established for the exceptional point of a PT-symmetric non-Hermitian Bose-Hubbard model. Black-Right-Pointing-Pointer The phase diagram shows that multiple regions of unbroken PT symmetry exist. Black-Right-Pointing-Pointer It is shown that even small on-site interaction can break the PT symmetry drastically.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Towards optimal extraction of cosmological information from nonlinear data

One of the main unsolved problems of cosmology is how to maximize the extraction of information from nonlinear data. If the data are nonlinear the usual approach is to employ a sequence of statistics (N-point statistics, counting statistics of clusters, density peaks or voids etc.), along with the corresponding covariance matrices. However, this approach is computationally prohibitive and has not been shown to be exhaustive in terms of information content. Here we instead develop a hierarchical Bayesian approach, expanding the likelihood around the maximum posterior of linear modes, which we solve for using optimization methods. By integrating out the modes using perturbative expansion of the likelihood we construct an initial power spectrum estimator, which for a fixed forward model contains all the cosmological information if the initial modes are gaussian distributed. We develop a method to construct the window and covariance matrix such that the estimator is explicitly unbiased and nearly optimal. We then generalize the method to include the forward model parameters, including cosmological and nuisance parameters, and primordial non-gaussianity. We apply the method in the simplified context of nonlinear structure formation, using either simplified 2-LPT dynamics or N-body simulations as the nonlinear mapping between linear and nonlinear density, and 2-LPT dynamics in the optimization steps used to reconstruct the initial density modes. We demonstrate that the method gives an unbiased estimator of the initial power spectrum, providing among other a near optimal reconstruction of linear baryonic acoustic oscillations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Wave Energy Converter Effects on Wave Fields: Evaluation of SNL-SWAN and Sensitivity Studies in Monterey Bay CA.

A modified version of an indust ry standard wave modeling tool was evaluated, optimized, and utilized to investigate model sensitivity to input parameters a nd wave energy converter ( WEC ) array deployment scenarios. Wave propagation was investigated d ownstream of the WECs to evaluate overall near - and far - field effects of WEC arrays. The sensitivity study illustrate d that wave direction and WEC device type we r e most sensitive to the variation in the model parameters examined in this study . Generally, the changes in wave height we re the primary alteration caused by the presence of a WEC array. Specifically, W EC device type and subsequently their size directly re sult ed in wave height variations; however, it is important to utilize ongoing laboratory studies and future field tests to determine the most appropriate power matrix values for a particular WEC device and configuration in order to improve modeling results .

Roberts, Jesse D.↗