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 991 records · Page 55

Combining cluster number counts and galaxy clustering

The abundance of clusters and the clustering of galaxies are two of the important cosmological probes for current and future large scale surveys of galaxies, such as the Dark Energy Survey. In order to combine them one has to account for the fact that they are not independent quantities, since they probe the same density field. It is important to develop a good understanding of their correlation in order to extract parameter constraints. We present a detailed modelling of the joint covariance matrix between cluster number counts and the galaxy angular power spectrum. We employ the framework of the halo model complemented by a Halo Occupation Distribution model (HOD). We demonstrate the importance of accounting for non-Gaussianity to produce accurate covariance predictions. Indeed, we show that the non-Gaussian covariance becomes dominant at small scales, low redshifts or high cluster masses. We discuss in particular the case of the super-sample covariance (SSC), including the effects of galaxy shot-noise, halo second order bias and non-local bias. We demonstrate that the SSC obeys mathematical inequalities and positivity. Using the joint covariance matrix and a Fisher matrix methodology, we examine the prospects of combining these two probes to constrain cosmological and HOD parameters. We find that the combination indeed results in noticeably better constraints, with improvements of order 20% on cosmological parameters compared to the best single probe, and even greater improvement on HOD parameters, with reduction of error bars by a factor 1.4-4.8. This happens in particular because the cross-covariance introduces a synergy between the probes on small scales. We conclude that accounting for non-Gaussian effects is required for the joint analysis of these observables in galaxy surveys.

79 ASTRONOMY AND ASTROPHYSICS↗

Assessment of CO 2 Storage Potential in Naturally Fractured Reservoirs With Dual–Porosity Models

Naturally Fractured Reservoirs (NFR’s) have received little attention as potential CO 2 storage sites. Two main facts deter from storage projects in fractured reservoirs: (1) CO 2 tends to be nonwetting in target formations and capillary forces will keep CO 2 in the fractures, which typically have low pore volume; and (2) the high conductivity of the fractures may lead to increased spatial spreading of the CO 2 plume. Numerical simulations are a powerful tool to understand the physics behind brine-CO 2 flow in NFR’s. Dual-porosity models are typically used to simulate multiphase flow in fractured formations. However, existing dual-porosity models are based on crude approximations of the matrix-fracture fluid transfer processes and often fail to capture the dynamics of fluid exchange accurately. Therefore, more accurate transfer functions are needed in order to evaluate the CO 2 transfer to the matrix. This work presents an assessment of CO 2 storage potential in NFR’s using dual-porosity models. We investigate the impact of a system of fractures on storage in a saline aquifer, by analyzing the time scales of brine drainage by CO 2 in the matrix blocks and the maximum CO 2 that can be stored in the rock matrix. A new model to estimate drainage time scales is developed and used in a transfer function for dual-porosity simulations. Furthermore, we then analyze how injection rates should be limited in order to avoid early spill of CO 2 (lost control of the plume) on a conceptual anticline model. Numerical simulations on the anticline show that naturally fractured reservoirs may be used to store CO 2 .

54 ENVIRONMENTAL SCIENCES↗

Modeling of Time-Dependent Strength Degradation of SiC/SiC Ceramic Matrix Composites via MAC/GMC Computer Code

Silicon carbide fiber reinforced silicon carbide (SiC/SiC) ceramic matrix composites (CMCs) display time-dependent strength degradation at intermediate temperatures (600 to 900 degrees Centigrade). This is generally believed to be an oxidation induced phenomenon. The understanding of the effect of temperature with environment (oxidation) is key towards development of SiC/SiC CMCs with a reliable load carrying capacity. Various theories have been proposed to explain the strength degradation. One suggests that the boron nitride (BN) coating deposited on the fibers oxidizes causing fusion of fibers. Another theory proposes that the SiC fibers are oxidized forming a silica scale leading to premature fiber failure. A more recent theory suggests that SiC fiber strength is intrinsically time-dependent due to slow crack growth of flaws in the fibers. An empirical model, termed as a "fiber classic model," which is based on a standard slow crack growth type power-law, has been implemented within NASA's micromechanics-based MAC/GMC computer code as a user routine. Model parameters for this "classic model" were calibrated from stress-rupture data of Hi-Nicalon (TM) monofilaments using the maximum likelihood estimation (MLE) technique. This new capability in the MAC/GMC computer code was then used to predict the stress-rupture behavior of Hi-Nicalon (TM) tows as well as 2-D SiC/SiC composites reinforced with Hi-Nicalon (TM) fibers. Results demonstrate that the MAC/GMC with this new capability successfully predicts the time-to-failure vs. applied stress within the intermediate temperature range at various scales as well as laminated composites in an oxidizing environment.

Environmental degradation↗

Micromechanical Modeling of the Thermal Expansion of Graphite/copper Composites with Nonuniform Microstructure

Two micromechanical models were developed to investigate the thermal expansion of graphite/copper (Gr/Cu) composites. The models incorporate the effects of temperature-dependent material properties, matrix inelasticity, initial residual stresses due to processing history, and nonuniform fiber distribution. The first model is based on the multiple concentric cylinder geometry, with each cylinder treated as a two-phase composite with a characteristic fiber volume fractions. By altering the fiber volume fraction of the individual cylinders, unidirectional composites with radially nonuniform fiber distributions can be investigated using this model. The second model is based on the inelastic lamination theory. By varying the fiber content in the individual laminae, composites with nonuniform fiber distribution in the thickness direction can be investigated. In both models, the properties of the individual regions (cylinders or laminae) are calculated using the method of cells micromechanical model. Classical incremental plasticity theory is used to model the inelastic response of the copper matrix at the microlevel. The models were used to characterize the effects of nonuniform fiber distribution on the thermal expansion of Gr/Cu. These effects were compared to the effects of matrix plasticity, choice of stress-free temperature, and slight fiber misalignment. It was found that the radially nonuniform fiber distribution has little effect on the thermal expansion of Gr/Cu but could become significant for composites with large fiber-matrix transverse CTE and Young's modulus mismatch. The effect of nonuniform fiber distribution in the through-thickness direction of a laminate was more significant, but only approached that of the stress-free temperature for the most extreme cases that include large amounts of bending. Subsequent comparison with experimental thermal expansion data indicated the need for more accurate characterization of the graphite fiber thermomechanical properties. Correlation with cyclic data revealed the presence of a mechanism not considered in the developed models. The predicted response did, however, exhibit ratcheting behavior that has been observed experimentally in Gr/Cu. Finally, simulation of the actual fiber distribution of particular specimens had little effect on the predicted thermal expansion.

Bednarcyk, Brett A.↗

Manifestation of Nonuniversality of Lepton Interactions in Spontaneously Violated Mirror Symmetry

Data from the LHCb experiments are indicative of a substantial distinction between the B → K (or K*) + e{sup +}e{sup −} and B → K (or K*) + μ{sup +}μ{sup −} branching ratios (April 2017). The branching ratio for the e+e− channel is substantially greater than that for the μ{sup +}μ{sup −} channel, whereas Standard Model (SM) calculations require that they be equal to each other. The above distinction may suggest the existence of a new interaction changing generations and discriminating between leptons that has couplings that are much greater than and are inverse in strength to the SM fermion couplings to the Higgs boson. Under conditions of spontaneously violated mirror symmetry, the coupling of SM particles to the second Higgs scalar that is inevitably present there and which is in principle heavy possesses precisely these properties. An inverse character of this coupling and its strength are not an additional hypothesis but a necessary condition for qualitatively reproducing, in addition, the observed charged-lepton mass hierarchy and the structure of weak lepton mixing—the Pontecorvo–Maki–Nakagawa–Sakata (PMNS) matrix. Within the mirror model being considered, all properties of the new interaction, including its inverse character, are due to the hierarchical character of the quark and charged-lepton mass spectrum.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Probability Bounds Analysis Applied to Multi-Purpose Crew Vehicle Nonlinearity

The Multi-Purpose Crew Vehicle (MPCV) Program Orion vehicle finite element model (FEM) was updated based on a modal test performed by Lockheed Martin. Due to nonlinearity observed in the test results, linear low force level (LL) and high force level (HL) FEMs were developed for use during various Space Launch System (SLS) flight regimes depending on expected forcing levels. Uncertainty models were derived for the combined MPCV and MPCV Stage Adaptor LL and HL Hurty/Craig-Bampton (HCB) components based on the MPCV structural test article Configuration 4 modal test-analysis correlation results. Subsequently, system-level uncertainty quantification analyses were performed using both models for various SLS flight configurations to determine the impact of the nonlinearity on important system metrics. The system metrics included both transfer functions associated with attitude control and dynamic loads associated with aerodynamic buffeting during ascent. In each case, an independent Monte Carlo (MC) analysis was performed, and no attempt was made to combine the results. The Hybrid Parametric Variation (HPV) method was used to develop the LL and HL MPCV HCB uncertainty models. The HPV method provides both parametric and non-parametric components of uncertainty. The non-parametric uncertainty accounts for the difference in model-form between the linearized analytical model and the corresponding linearized component test results in the form of mode shapes and frequencies at that force level. This linear model-form uncertainty is implemented in the HPV method using random matrix theory. However, the HPV uncertainty models developed for the linear LL and HL MPCV components do not account for the nonlinearity in the MPCV. With respect to the linearized models, this nonlinearity is also an uncertainty in model form, but in this case, it must be treated independently as an epistemic uncertainty. It represents a lack of knowledge, in contrast to an aleatory uncertainty due to the randomness of a variable. In the case of an epistemic variable, the true value is unknown, only the interval within which it lies is known. Epistemic uncertainty can be reduced with increased knowledge, while in general, aleatory uncertainty cannot. This work combines the epistemic uncertainty due to the MPCV nonlinearity with the parametric and non-parametric uncertainty within the HPV method using a second order propagation approach. The LL and HL test data is augmented with surrogate test data derived from a nonlinear MPCV representation. The impact of the MPCV nonlinearity on system response statistics is determined using a series of cumulative distribution functions in the form of a horsetail plot, or p-box. This results in an interval of probabilities for a specific response value, or an interval of response values at a specific probability.

Daniel C Kammer↗

Statistical scaling relationships and size effects in the strength and creep rupture of fibrous composites

In this presentation we discuss a new theoretical model and supporting experimental results for the strength and lifetime in creep rupture of unidirectional, carbon fiber/epoxy matrix composites at ambient conditions. First we review the 'standard' Weibull/power-law methodology that has been standard practice. Then we discuss features of a recent model which build on the statistical aspects of fiber strength, micromechanical aspects of stress transfer around fiber breaks, and time-dependent creep of the matrix. The model is applied to 'microcomposites' consisting of seven fibers in a matrix for which strength and creep-rupture data are available. The model yields Weibull distributions in an envelope format for both strength and lifetime. The respective shape, scale and power-law parameters depend on such parameters as the Weibull shape parameter for fiber strength, the exponent for matrix creep, the effective load transfer length (which grows in time due to matrix creep) and the critical cluster size for failed fibers. The experimental results are consistent with the theory, though time-dependent debonding appears to be part of the failure process.

Phoenix, S. Leigh↗

Gradient-based constrained optimization using a database of linear reduced-order models

A methodology grounded in model reduction is presented for accelerating the gradient-based solution of a family of linear or nonlinear constrained optimization problems where the constraints include at least one linear Partial Differential Equation (PDE). A key component of this methodology is the construction, during an offline phase, of a database of pointwise, linear, Projection-based Reduced-Order Models (PROM)s associated with a design parameter space and the linear PDE(s). A parameter sampling procedure based on an appropriate saturation assumption is proposed to maximize the efficiency of such a database of PROMs. A real-time method is also presented for interpolating at any queried but unsampled parameter vector in the design parameter space the relevant sensitivities of a PROM. The practical feasibility, computational advantages, and performance of the proposed methodology are demonstrated for several realistic, nonlinear, aerodynamic shape optimization problems governed by linear aeroelastic constraints.

97 MATHEMATICS AND COMPUTING↗

Elastic and Piezoelectric Properties of Boron Nitride Nanotube Composites: Finite Element Model - Part II

This paper is the second part of a two-part series where the first part presents a molecular dynamics model of a single Boron Nitride Nanotube (BNNT) and this paper scales up to multiple BNNTs in a polymer matrix. This paper presents finite element (FE) models to investigate the effective elastic and piezoelectric properties of (BNNT) nanocomposites. The nanocomposites studied in this paper are thin films of polymer matrix with aligned co-planar BNNTs. The FE modelling approach provides a computationally efficient way to gain an understanding of the material properties. We examine several FE models to identify the most suitable models and investigate the effective properties with respect to the BNNT volume fraction and the number of nanotube walls. The FE models are constructed to represent aligned and randomly distributed BNNTs in a matrix of resin using 2D and 3D hollow and 3D filled cylinders. The homogenisation approach is employed to determine the overall elastic and piezoelectric constants for a range of volume fractions. These models are compared with an analytical model based on Mori-Tanaka formulation suitable for finite length cylindrical inclusions. The model applies to primarily single-wall BNNTs but is also extended to multi-wall BNNTs, for which preliminary results will be presented. Results from the Part 1 of this series can help to establish a constitutive relationship for input into the finite element model to enable the modeling of multiple BNNTs in a polymer matrix.

Kim, H. Alicia↗

Analysis of the phase control of the ITER ICRH antenna array. Influence on the load resilience and radiated power spectrum

The paper analyses how the phasing of the ITER ICRH 24 strap array evolves from the power sources up to the strap currents of the antenna. The study of the phasing control and coherence through the feeding circuits with prematching and automatic matching and decoupling network is made by modeling starting from the TOPICA matrix of the antenna array for a low coupling plasma profile and for current drive phasing (worst case for mutual coupling effects). The main results of the analysis are: (i) the strap current amplitude is well controlled by the antinode V{sub max} amplitude of the feeding lines, (ii) the best toroidal phasing control is done by the adjustment of the mean phase of V{sub max} of each poloidal straps column, (iii) with well adjusted system the largest strap current phasing error is ±20°, (iv) the effect on load resilience remains well below the maximum affordable VSWR of the generators, (v) the effect on the radiated power spectrum versus k{sub //} computed by means of the coupling code ANTITER II remains small for the considered cases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Tungsten isotopes and the origin of chondrules and chondrites

Chondrules and matrix from carbonaceous chondrites exhibit complementary nucleosynthetic W isotope anomalies that result from the depletion of a metallic s-process carrier in the chondrules, and the enrichment of this carrier in the matrix. The complementarity is difficult to reconcile with an origin of chondrules in protoplanetary impacts and also with models in which chondrules and matrix formed independently of each other in distinct regions of the disk. Instead, the complementarity indicates that chondrules formed by localized melting of dust aggregates in the solar nebula. The Hf–W ages for metal-silicate fractionation in CV and CR chondrites are 2.2 ± 0.8 Ma and 3.6 ± 0.6 Ma after formation of Ca-Al-rich inclusions, and are indistinguishable from Al–Mg ages for CV and CR chondrules. The good agreement between these ages strongly suggests that 26 Al was homogeneously distributed in the solar protoplanetary disk and that therefore Al–Mg ages are chronologically meaningful. The concordant Al–Mg and Hf–W ages reveal that chondrule formation (as dated by Al–Mg) was associated with metal-silicate fractionation (as dated by Hf–W), both within a given chondrite but also among the different subgroups of ordinary chondrites. These age data indicate that chondrules from a given chondrite group formed in a narrow time interval of <1 Ma, and that chondrule formation and chondrite accretion were closely linked in time and space. The rapid accretion of chondrules into a chondrite parent body is consistent with the isotopic complementarity, which requires that neither chondrules nor matrix were lost prior to chondrite accretion. Here combined, these observations suggest that chondrule formation was an important step in the accretion of planetesimals.

58 GEOSCIENCES↗

Forward problem solution as the operator of filtered and back projection matrix to reconstruct the various method of data collection and the object element model in electrical impedance tomography

Back projection reconstruction has been implemented to get the dynamical image in electrical impedance tomography. However the implementation is still limited in method of adjacent data collection and circular object element model. The study aims to develop the methods of back projection as reconstruction method that has the high speed, accuracy, and flexibility, which can be used for various methods of data collection and model of the object element. The proposed method uses the forward problem solution as the operator of filtered and back projection matrix. This is done through a simulation study on several methods of data collection and various models of the object element. The results indicate that the developed method is capable of producing images, fastly and accurately for reconstruction of the various methods of data collection and models of the object element.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Estimating the Pyrolysis Kinetic Parameters of Coal, Biomass, and Their Blends: A Comparative Study

In this work, the pyrolysis kinetic parameters of two coal ranks (DECS-25 lignite and DECS-38 sub-bituminous), two biomass materials (corn stover and switchgrass), and their respective blends were investigated at various heating rates ranging between 5 and 40 °C/min using thermogravimetric analysis. Complex models for devolatilization of the feedstocks were solved for obtaining and predicting the global kinetic parameters. Distributed activation energy model (method 1) and matrix inversion algorithm (method 2) were utilized and compared for this purpose. The results indicate that the matrix inversion algorithm predicts the kinetic parameters such that the weight loss characteristics can be best represented for both single fuels as well as that of blended materials. The algorithm can also be used for determining the number of reactions occurring in the devolatilization temperature interval. The number of reactions occurring during the devolatilization of blended materials falls between those that occur during the devolatilization of single fuels and the number of reactions gradually decrease with increase of biomass concentration in the blend. In addition, weight loss characteristics of fuel blends at unknown heating rates can be effectively predicted within 1% error through the use of this algorithm.

01 COAL, LIGNITE, AND PEAT↗

Grid Impact Analysis of Heavy-Duty Electric Vehicle Charging Stations

This paper presents a grid impact analysis of heavy-duty electric vehicle (EV) charging stations. In this work, heavy-duty EVs have battery capacities high enough to provide a range of 250-500 miles on a single charge, such as long-haul trucks. Heavy-duty EVs will require extremely fast charging rates to reduce charging time and will induce very high charging loads (at the multiple-megawatt scale) if several vehicles charge at the same time. Therefore, analysis is needed to understand the impact of charging station loads on the electric power grid and set the baseline for developing mitigation plans and necessary system upgrades. We develop a systematic procedure to analyze the potential impact of the placement of charging stations on the grid. Charging load is modeled using a DC fast-charging station model. A voltage load sensitivity matrix approach is leveraged to investigate the challenges of placing charging stations on the feeder. Given the charging load profiles and suggested charging station locations, time-series simulations are performed on various connection points on the feeder to understand the impact. The analysis is performed on both the IEEE 34-bus system and a realistic feeder from California. Initial mitigation solutions are developed based on insights from this analysis.

47 OTHER INSTRUMENTATION↗

Scale-Dependent Fracture-Matrix Interactions And Their Impact on Radionuclide Transport - Final Report

Matrix diffusion and adsorption within a rock matrix are widely regarded as important mechanisms for retarding the transport of radionuclides and other solutes in fractured rock (e.g., Neretnieks, 1980; Tang et al., 1981; Maloszewski and Zuber, 1985; Novakowski and Lapcevic, 1994; Jardine et al., 1999; Zhou and Xie, 2003; Reimus et al., 2003a,b). When remediation options are being evaluated for old sources of contamination, where a large fraction of contaminants reside within the rock matrix, slow diffusion out of the matrix greatly increases the difficulty and timeframe of remediation. Estimating the rates of solute exchange between fractures and the adjacent rock matrix is a critical factor in quantifying immobilization and/or remobilization of DOE-relevant contaminants within the subsurface. In principle, the most rigorous approach to modeling solute transport with fracture-matrix interaction would be based on local-scale coupled advection-diffusion/dispersion equations for the rock matrix and in discrete fractures that comprise the fracture network (Discrete Fracture Network and Matrix approach, hereinafter referred to as DFNM approach), fully resolving aperture variability in fractures and matrix property heterogeneity. However, such approaches are computationally demanding, and thus, many predictive models rely upon simplified models. These models typically idealize fracture rock masses as a single fracture or system of parallel fractures interacting with slabs of porous matrix or as a mobile-immobile or multi-rate mass transfer system. These idealizations provide tractable approaches for interpreting tracer tests and predicting contaminant mobility, but rely upon a fitted effective matrix diffusivity or mass-transfer coefficients. However, because these fitted parameters are based upon simplified conceptual models, their effectiveness at predicting long-term transport processes remains uncertain. Evidence of scale dependence of effective matrix diffusion coefficients obtained from tracer tests highlights this point and suggests that the underlying mechanisms and relationship between rock and fracture properties are not fully understood in large complex fracture networks. In this project, we developed a high-resolution DFN model of solute transport in fracture networks to explore and quantify the mechanisms that control transport in complex fracture networks and how these may give rise to observed scale-dependent matrix diffusion coefficients. Results demonstrate that small scale heterogeneity in the flow field caused by local aperture variability within individual fractures can lead to long-tailed breakthrough curves indicative of matrix diffusion, even in the absence of interactions with the fracture matrix. Furthermore, the temporal and spatial scale dependence of these processes highlights the inability of short-term tracer tests to estimate transport parameters that will control long-term fate and transport of contaminants in fractured aquifers.

58 GEOSCIENCES↗

Emergence of p+ip superconductivity in two-dimensional doped Dirac systems

Searching for the $p + ip$ superconducting (SC) state has become a fascinating subject in condensed matter physics recently, as a dream application awaiting in topological quantum computation. Unfortunately, so far there is no universal principle for realizing $p + ip$ in generic solid-state systems. Here we report a theoretical discovery of a $p + ip$ SC ground state (coexisting with ferromagnetic order) in the honeycomb lattice Hubbard model in the extremely strong-coupling limit (e.g., infinite U) at low doping $(δ < 0.2)$, by using both the state-of-art Grassmann tensor product state approach and a continuum quantum field theory approach. Our discovery suggests a mechanism for the $p + ip$ SC state in generic strongly correlated systems based on spin-charge separation and the charge current–current coupling scenario, which opens a door towards experimental realization. The $p + ip$ SC state has an instability towards a potential non-Fermi liquid with a large but finite U . Nevertheless, by applying an in-plane Zeeman field, such a $p + ip$ SC state can be stabilized with finite U in a very wide range of doping. Relevant realistic materials are also proposed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Phase transitions in the Rényi entropies of a 2+1D large- N interacting vector model

Here we consider the reduced density matrix for a disk in the ground state of the interacting fixed points of the large- N O(N) vector model in 2 + 1 dimensions. Using the map to the free energy on $\mathbb H_2 × S^1$, we show that there is an instability in the Rényi entropies at n = 1 in the N → ∞ limit, indicating a phase transition at or near this Rényi parameter. We close with a discussion of the finite- N case.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Ground-state properties of the periodic Anderson model

The ground-state energy, hybridization matrix element, local moment, and spin-density correlations of a one-dimensional, finite-chain, periodic, symmetric Anderson model are obtained by numerical simulations and compared with perturbation theory and strong-coupling results. It is found that the local f-electron spins are compensated by correlation with other f-electrons as well as band electrons leading to a nonmagnetic ground state.

Blankenbecler, R.↗