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 235 records · Page 13

Piezo-Antenna: An Electromagnetic-Ultrasonic Data Transfer System

Electronic systems contained within sealed metallic enclosures, such as nuclear waste containers where wire penetration is not an option, require a through-wall data transfer technology to communicate. One such technology utilizes piezoelectric ultrasonic transducers to communicate through the barrier using elastic waves. Currently demonstrated ultrasonic through-wall data communication techniques have relied on having a direct interface with the ultrasonic transducers on either side. Here, this work studies the coupling of an antenna to such an ultrasonic communication system. A computationally efficient model of a piezobarrier ultrasonic communication system and antenna using a simplified antenna model and a transfer matrix piezobarrier model is utilized to perform a representative parameter space search. The system efficiency and bandwidth trends are shown for the system parameters as well as operating frequency. The tradeoff between efficiency and bandwidth is also presented. A system configuration within the parameter space is fabricated and demonstrated to achieve a data rate of up to 100 kb/s.

Antenna↗

Finite Element Modeling of the Fiber-Matrix Interface in Polymer Composites

Polymer composites are used in numerous industries due to their high specific strength and high specific stiffness. Composites have markedly different properties than both the reinforcement and the matrix. Of the several factors that govern the final properties of the composite, the interface is an important factor that influences the stress transfer between the fiber and matrix. The present study is an effort to characterize and model the fiber-matrix interface in polymer matrix composites. Finite element models were developed to study the interfacial behavior during pull-out of a single fiber in continuous fiber-reinforced polymer composites. A three-dimensional (3D) unit-cell cohesive damage model (CDM) for the fiber/matrix interface debonding was employed to investigate the effect of interface/sizing coverage on the fiber. Furthermore, a two-dimensional (2D) axisymmetric model was used to (a) analyze the sensitivity of interface stiffness, interface strength, friction coefficient, and fiber length via a parametric study; and (b) study the shear stress distribution across the fiber-interface-matrix zone. It was determined that the force required to debond a single fiber from the matrix is three times higher if there is adequate distribution of the sizing on the fiber. The parametric study indicated that cohesive strength was the most influential factor in debonding. Moreover, the stress distribution model showed the debonding mechanism of the interface. It was observed that the interface debonded first from the matrix and remained in contact with the fiber even when the fiber was completely pulled out.

36 MATERIALS SCIENCE↗

Mobile Sensing for Wind Field Estimation in Wind Farms: Preprint

This paper introduces a novel approach for estimating the wind field over an entire wind farm using a mobile sensor to collect limited amounts of data. The proposed method estimates the boundary conditions of a simplified turbine wake model by computing the model sensitivity matrix and using a recursive least-squares algorithm to recover the model parameters from the wind field measurements. To address the fact that it is not practical to take measurements across the entire wind farm, the proposed method classifies each area on the map based on its sensitivity to parameter variations. This classification is then used to generate a suitable path for a mobile sensor, which is charged with collecting data for the recursive least-squares algorithm. The proposed framework can successfully estimate the model boundary conditions using just the measurements collected along the path of the mobile sensor. This preliminary result paves the way for using real-time wind field estimates for the coordinated control of all the turbines within a wind farm.

estimation↗

Electron-impact excitation of Fe{sup 2+}: A comparison of intermediate coupling frame transformation, Breit-Pauli and Dirac R-matrix calculations

Modeling the spectral emission of low-charge iron group ions enables the diagnostic determination of the local physical conditions of many cool plasma environments such as those found in H II regions, planetary nebulae, active galactic nuclei, etc. Electron-impact excitation drives the population of the emitting levels and, hence, their emissivities. By carrying-out Breit-Pauli and intermediate coupling frame transformation (ICFT) R-matrix calculations for the electron-impact excitation of Fe{sup 2+}, which both use the exact same atomic structure and the same close-coupling expansion, we demonstrate the validity of the application of the powerful ICFT method to low-charge iron group ions. This is in contradiction to the finding of Bautista et al., who carried-out ICFT and Dirac R-matrix calculations for the same ion. We discuss possible reasons.

79 ASTRONOMY AND ASTROPHYSICS↗

SALSA3D: A Tomographic Model of Compressional Wave Slowness in the Earth’s Mantle for Improved Travel-Time Prediction and Travel-Time Prediction Uncertainty

The task of monitoring the Earth for nuclear explosions relies heavily on seismic data to detect, locate, and characterize suspected nuclear tests. In this study, motivated by the need to locate suspected explosions as accurately and precisely as possible, we developed a tomographic model of the compressional wave slowness in the Earth’s mantle with primary focus on the accuracy and precision of travel-time predictions for P and Pn ray paths through the model. Path-dependent travel-time prediction uncertainties are obtained by computing the full 3D model covariance matrix and then integrating slowness variance and covariance along ray paths from source to receiver. Path-dependent travel-time prediction uncertainties reflect the amount of seismic data that was used in tomography with very low values for paths represented by abundant data in the tomographic data set and very high values for paths through portions of the model that were poorly sampled by the tomography data set. The pattern of travel-time prediction uncertainty is a direct result of the off-diagonal terms of the model covariance matrix and underscores the importance of incorporating the full model covariance matrix in the determination of travel-time prediction uncertainty. In addition, the computed pattern of uncertainty differs significantly from that of 1D distance-dependent travel-time uncertainties computed using traditional methods, which are only appropriate for use with travel times computed through 1D velocity models.

58 GEOSCIENCES↗

Microstructure in matrix elements

We investigate the simple model of Pennington, Shenker, Stanford and Yang for modeling the density matrix of Hawking radiation, but further include dynamics for EOW branes behind the horizon. This allows interactions that scatter one interior state to another, and also allows EOW loops. At strong coupling, we find that EOW states are no longer random; the ensemble has collapsed, and coupling constants encode the microscopic matrix elements of Hawking radiation. This suggests strong interior dynamics are important for understanding evaporating black holes, without any ensemble average. In this concrete model the density matrix of the radiation deviates from the thermal state, small off-diagonal fluctuations encode equivalences between naively orthogonal states, and bound the entropy from above. For almost evaporated black holes the off-diagonal terms become as large as the diagonal ones, eventually giving a pure state. We also find the unique analytic formula for all Renyi entropies.

2D gravity↗

An ICME Modeling Framework for Titanium/Tungsten-Carbide Metal Matrix Composites

This report describes a collaborative project to develop a validated, predictive model for the high temperature mechanical properties of a titanium-matrix, tungsten-carbide/cobalt-reinforced metal matrix composite. The modeling approach was to first develop a detailed, microstructural model linking the material structure and the interfacial debonding properties to the effective properties of the material. The project then completed a throughput simulation campaign to generate a large number of simulations for discrete microstructures and different debonding parameters. Finally, the project trained a fast, Gaussian process surrogate model against this simulation database to provide a quick model linking the material compositions, structure, and processing parameters to the resulting material properties. This model was validated against high temperature tensile test data on a few particular composite compositions. The tests validate the model predictions for ultimate tensile strength and uniform elongation/ductility, meaning the final surrogate model can now be used to tune the material composition and processing parameters to identify optimal composite compositions for particular applications.

36 MATERIALS SCIENCE↗

A semi-analytical method for simulating matrix diffusion in numerical transport models

We report a semi-analytical approximation for transient matrix diffusion is developed for use in numerical contaminant transport simulators. This method is an adaptation and extension of the heat conduction method of Vinsome and Westerveld (1980) used to simulate heat losses during thermally enhanced oil recovery. The semi-analytical method is used in place of discretization of the low permeability materials, and it represents the concentration profile in the low permeability materials with a fitting function that is adjusted in each element at each time-step. The resulting matrix diffusion fluxes are added to the numerical model as linear concentration-dependent source/sink terms. Since only the high permeability zones need to be discretized, the numerical formulation is extremely efficient compared to traditional approaches that require discretization of both the high and low permeability zones. The semi-analytical method compares favorably with the analytical solution for transient one-dimensional diffusion with first order decay, with a two-layer aquifer/aquitard solution, with the solution for transport in a fracture with matrix diffusion and decay, and with a fully numerical solution for transport in a thin sand zone bounded by clay with variable decay rates.

54 ENVIRONMENTAL SCIENCES↗

Direct Feed High-Level Waste APPS Model Glass Testing (DFHLW APPS) Matrix

This report summarizes the data collected during the batching and melting of the Direct Feed High-Level Waste APPS Model Glass Matrix (DFHLW APPS) to serve as a quality-assured validation of the Aspen Process Performance Simulation (APPS) formulation method. Of 15 glasses tested, 12 satisfied all target property constraints. Two glasses, APPS-05 and -06, formed nepheline on canister centerline cooling heat-treatment and failed the Product Consistency Test response limits. Glass APPS-07-2 formed unacceptably high concentrations of crystals (primarily Na3Nd(PO4)2) when heat treated at 950 °C. All other glasses were found to be satisfactory. The measured property values were compared to predicted values from a set of current models. In many cases the current models were found to be inadequate for design of DFHLW glasses. These models are being adjusted to correct for mispredictions. Other models, e.g., density, toxicity characteristic leaching procedure, and sulfur solubility, are adequate for formulation of DFHLW glasses.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Fast matrix algebra for Bayesian model calibration

In Bayesian model calibration, evaluation of the likelihood function usually involves finding the inverse and determinant of a covariance matrix. When Markov Chain Monte Carlo (MCMC) methods are used to sample from the posterior, hundreds of thousands of likelihood evaluations may be required. In this paper, we demonstrate that the structure of the covariance matrix can be exploited, leading to substantial time savings in practice. Here, we also derive two simple equations for approximating the inverse of the covariance matrix in this setting, which can be computed in near-quadratic time. The practical implications of these strategies are demonstrated using a simple numerical case study and the "quack" R package. For a covariance matrix with 1000 rows, application of these strategies for a million likelihood evaluations leads to a speedup of roughly 4000 compared to the naive implementation

97 MATHEMATICS AND COMPUTING↗

Model predictive control for matrix converter operating in current control mode with load current estimation

A matrix converter system operating in current control mode is provided. The matrix converter includes a switching matrix coupled between a low voltage side and a high voltage side, wherein the matrix converter is coupled at its low voltage side to a generator for receiving an input power including an input current and transforming the input power into an output power at its high voltage side, wherein a load is coupled to the high voltage side. The control system includes a load observer and a model predictive controller (MPC). The load observer is configured to estimate a load current that flows to the load from the high voltage side of the matrix converter as a function of a switching state of the switching matrix, an output voltage output at the high voltage side of the matrix converter, and the input current. The MPC is configured to select the switching state of the switching matrix to meet one or more control objectives defined by minimization of a multi-objective function that tracks the output voltage and the input current using the estimated load current.

24 POWER TRANSMISSION AND DISTRIBUTION↗

An ICME Modeling Framework for Metal Matrix Composites Focusing on Ultrahigh Temperature Matrix Material and Tungsten Carbide Reinforcement Particulate (Final CRADA Report)

Ultra High Temperature Metal Matrix Composites (UHT-MMCs) are of interest to the due to their potential defense and aerospace applications. These materials consist of a metal matrix reinforced with a stiff ceramic or cermet. Based on a preliminary literature survey, this project will focused on a pure titanium (Ti) matrix reinforced with tungsten carbide-cobalt (WC-Co). UHT-MMC technology is still at a nascent stage of development and requires extensive experimental work to optimize the structure and resulting material properties of a particular UHT-MMC material system. An accurate and complete high-performance computing (HPC) model of the Ti/WC-Co material system could greatly accelerate the development and deployment Ti/WC-Co materials by connecting processing parameters to material properties. This type of model could then be used in an integrated computational materials engineering approach to optimize the material structure to achieve a targeted level of performance.

36 MATERIALS SCIENCE↗

Phase-Field Modeling of Mechanical Damages in Ceramic Matrix Composites

Developed a phase-field model for mechanical damages in CMCs, which incorporates the CMC microstructures, matrix cracking, fiber breakage, and interfacial sliding. Two types of mechanisms, fiber bridging and fiber pull-out, are considered. The obtained simulation results agree with experimental observations and an analytical solution. Simulation results suggest that the performance of CMCs would be enhanced with thicker fibers, longer fibers, and higher fiber density. Opposite trends of interfacial sliding resistances are suggested for the two types of situations. In reality, a mixture of the two situations may exist, and then an intermediate interfacial sliding resistance may be optimal.

Xue, Fei↗

Scale-Dependent Fracture-Matrix Interactions and Their Impact on Radionuclide Transport: Development of efficient particle-tracking methods

Matrix Diffusion and Adsorption within a rock matrix are important mechanisms for retarding transport of radionuclides in fractured rock. Due to computational limitations and difficulties in characterizing complex subsurface systems, diffusive exchange between a fracture network and surrounding rock matrix is often modeled using simplified conceptual representations. There is significant uncertainty in “effective” parameters used in these models, such as the “effective matrix diffusivity”. Often, these parameters are estimated by fitting sparse breakthrough data, and estimated values fall outside meaningful ranges, because simplified interpretive models do not consider complex three-dimensional flow. There is limited understanding of the relationship between the effective parameters and rock mass characteristics including network structure and matrix properties. There is also evidence for an apparent scale-dependence in “effective matrix diffusion” coefficients. These observations raise questions on whether fracture-matrix interaction parameters estimated from small-scale tracer tests can be used for predicting radionuclide fate and transport at the scale of DOE field sites. High-resolution three-dimensional Discrete-Fracture-Network-Matrix (DFNM) models based on well-defined local scale transport equations can help to address some of these questions. Due to tremendous advances in computational technology over the last 10 years, DFNM modeling in relatively large domains is now feasible. The overarching objective of our research is to use DFNM modeling to improve fundamental understanding of how effective parameters in conceptual models are related to fracture network structure and matrix properties. An advanced three-dimensional DFNM model is being developed, which combines upscaled particle-tracking algorithms for fracture-matrix interaction and a parallel fracture-network flow simulator. The particle-tracking algorithms allow complexity in flow fields at different scales, and track transport across fracture-matrix interfaces based on rigorous local approximations to the transport equations. This modeling approach can incorporate aperture variability, multi-scale preferential flow and matrix heterogeneity. We developed efficient particle-tracking methods for handling matrix diffusion and adsorption on fracture walls and demonstrated their efficiency for use within the context of large-scale complex fracture network models with variability in apertures across a network of fractures and within individual fractures.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

The Effect of the Prior and the Experimental Design on the Inference of the Precision Matrix in Gaussian Chain Graph Models

Here, we investigate whether (and how) experimental design could aid in the estimation of the precision matrix in a Gaussian chain graph model, especially the interplay between the design, the effect of the experiment and prior knowledge about the effect. Estimation of the precision matrix is a fundamental task to infer biological graphical structures like microbial networks. We compare the marginal posterior precision of the precision matrix under four priors: flat, conjugate Normal-Wishart, Normal-MGIG and a general independent. Under the flat and conjugate priors, the Laplace-approximated posterior precision is not a function of the design matrix rendering useless any efforts to find an optimal experimental design to infer the precision matrix. In contrast, the Normal-MGIG and general independent priors do allow for the search of optimal experimental designs, yet there is a sharp upper bound on the information that can be extracted from a given experiment. We confirm our theoretical findings via a simulation study comparing (i) the KL divergence between prior and posterior and (ii) the Stein’s loss difference of MAPs between random and no experiment. Our findings provide practical advice for domain scientists conducting experiments to better infer the precision matrix as a representation of a biological network.

54 ENVIRONMENTAL SCIENCES↗

Osteoporosis Treatments Affect Bone Matrix Maturation in a Rat Model of Induced Cortical Remodeling

ABSTRACT To test how osteoporosis drugs affect bone matrix maturation during cortical bone remodeling, 72 pregnant rats were switched from a 0.4% to a 0.01% calcium diet at parturition for a 23‐day lactation period. At weaning, eight dams were sacrificed to establish baseline values, while the remaining dams were returned to 0.4% calcium and treated with vehicle (saline), sodium fluoride (NaF), zoledronic acid (ZA), or sclerostin antibody (Scl‐Ab) for either 7 or 28 days (eight animals per group per time point). Femora were examined by μCT, dynamic histomorphometry, Fourier transform infrared imaging, and three‐point bending of notched specimens. Cortical porosity decreased in all groups from baseline to day 28. Intracortical mineralizing surface (MS/BS) and mineral apposition rate (MAR), as well as the mineral‐to‐matrix ratio were unaffected by treatment, but intracortical crystallinity was increased in the ZA group at day 10 compared with vehicle. Cortical area increased in all groups over 28 days mainly because of an addition of bone at the endocortical surface. Endocortical MS/BS did not vary among the groups, but endocortical MAR was suppressed in the NaF group at day 2 and elevated in the Scl‐Ab group at day 4 compared with vehicle. Endocortical mineral‐to‐matrix ratio was increased at days 5 and 10 following NaF treatment and endocortical crystallinity was increased at day 5 following ZA treatment compared with vehicle. Fracture toughness did not differ among the groups. Thus, the treatments affected matrix maturation more strongly at the endocortical then intracortical envelope. In this model of induced remodeling, the bone formation phase is synchronized at multiple sites, facilitating study of the effects of drugs or other bone‐targeting agents on matrix maturation independent of their effects on the initiation of remodeling. © 2020 The Authors. JBMR Plus published by Wiley Periodicals, Inc. on behalf of American Society for Bone and Mineral Research.

59 BASIC BIOLOGICAL SCIENCES↗

On time-dependent backgrounds in 1 + 1 dimensional string theory

In perturbative string theory, one is generally interested in asymptotic observables, such as the S-matrix in flat spacetime, and boundary correlation functions in anti-de Sitter spacetime. However, there are backgrounds in which such observables do not exist. We study examples of such backgrounds in 1 + 1 dimensional string theory. In these examples, the Liouville wall accelerates and can become spacelike in the past and/or future. When that happens, the corresponding null infinity, at which the standard scattering states are defined, is shielded by the Liouville wall. We compute scattering and particle production amplitudes in these backgrounds in the region in parameter space where the wall remains timelike, and discuss the continuation of this picture to the spacelike regime. We also discuss the physics from the point of view of the dynamics of free fermions in backgrounds with a time-dependent Fermi surface.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗