SEARCH · Engineering Papers

Results for “Matrix factorization”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 955 records · Page 53

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Time-Lapse Integration at FWU: Fluids, Rock Physics, Numerical Model Integration, and Field Data Comparison

We present the current status of time-lapse seismic integration at the Farnsworth (FWU) CO2 WAG (water-alternating-gas) EOR (Enhanced Oil Recovery) project at Ochiltree County, northwest Texas. As a potential carbon sequestration mechanism, CO2 WAG projects will be subject to some degree of monitoring and verification, either as a regulatory requirement or to qualify for economic incentives. In order to evaluate the viability of time-lapse seismic as a monitoring method the Southwest Partnership (SWP) has conducted time-lapse seismic monitoring at FWU using the 3D Vertical Seismic Profiling (VSP) method. The efficacy of seismic time-lapse depends on a number of key factors, which vary widely from one application to another. Most important among these are the thermophysical properties of the original fluid in place and the displacing fluid, followed by the petrophysical properties of the rock matrix, which together determine the effective elastic properties of the rock fluid system. We present systematic analysis of fluid thermodynamics and resulting thermophysical properties, petrophysics and rock frame elastic properties, and elastic property modeling through fluid substitution using data collected at FWU. These analyses will be framed in realistic scenarios presented by the FWU CO2 WAG development. The resulting fluid/rock physics models will be applied to output from the calibrated FWU compositional reservoir simulation model to forward model the time-lapse seismic response. Modeled results are compared with field time-lapse seismic measurements and strategies for numerical model feedback/update are discussed. While mechanical effects are neglected in the work presented here, complementary parallel studies are underway in which laboratory measurements are introduced to introduce stress dependence of matrix elastic moduli.

4D↗

NASA NTRS1989

Shear testing of fiber reinforced metal matrix composites

This paper outlines the elements of a combined experimental/analytical methodology for accurate shear characterization of unidirectional composites in the linear and nonlinear range, with particular attention devoted to metal matrix composites. It is illustrated that consistent results can be obtained for a large class of composites from two commonly employed shear test methods currently in use by composites researchers when the influence of various factors that affect the determination of the actual shear response is properly analyzed. These factors include the effects of material anisotropy, specimen geometry, manner of load introduction, and test fixture design on the stress and deformation fields in the test section of off-axis and Iosipescu specimens. Common errors associated with the measurement of deformation fields and calculation of stress fields are discussed and quantified. Particular problems in the determination of the shear response of unidirectional boron/aluminum using the Iosipescu test are illustrated and discussed.

Pindera, Marek-Jerzy↗

NASA NTRS1993

Finite-Difference Algorithms For Computing Sound Waves

Governing equations considered as matrix system. Method variant of method described in "Scheme for Finite-Difference Computations of Waves" (ARC-12970). Present method begins with matrix-vector formulation of fundamental equations, involving first-order partial derivatives of primitive variables with respect to space and time. Particular matrix formulation places time and spatial coordinates on equal footing, so governing equations considered as matrix system and treated as unit. Spatial and temporal discretizations not treated separately as in other finite-difference methods, instead treated together by linking spatial-grid interval and time step via common scale factor related to speed of sound.

Davis, Sanford↗

NASA NTRS1997

Microphysical Retrievals Over Stratiform Rain Using Measurements from an Airborne Dual-Wavelength Radar-Radiometer

The need to understand the complementarity of the radar and radiometer is important not only to the Tropical Rain Measuring Mission (TRMM) program but to a growing number of multi-instrumented airborne experiment that combine single or dual-frequency radars with multichannel radiometers. The method of analysis used in this study begins with the derivation of dual-wavelength radar equations for the estimation of a two-parameter drop size distribution (DSD). Defining a "storm model" as the set of parameters that characterize snow density, cloud water, water vapor, and features of the melting layer, then to each storm model there will usually correspond a set of range-profiled drop size distributions that are approximate solutions of the radar equations. To test these solutions, a radiative transfer model is used to compute the brightness temperatures for the radiometric frequencies of interest. A storm model or class of storm models is considered optimum if it provides the best reproduction of the radar and radiometer measurements. Tests of the method are made for stratiform rain using simulated storm models as well as measured airborne data. Preliminary results show that the best correspondence between the measured and estimated radar profiles usually can be obtained by using a moderate snow density (0.1-0.2 g/cu cm), the Maxwell-Garnett mixing formula for partially melted hydrometeors (water matrix with snow inclusions), and low to moderate values of the integrated cloud liquid water (less than 1 kg/sq m). The storm-model parameters that yield the best reproductions of the measured radar reflectivity factors also provide brightness temperatures at 10 GHz that agree well with the measurements. On the other hand, the correspondence between the measured and modeled values usually worsens in going to the higher frequency channels at 19 and 34 GHz. In searching for possible reasons for the discrepancies, It is found that changes in the DSD parameter Mu, the radar constants, or the path-integrated attenuation can affect the high frequency channels significantly. In particular, parameters that cause only modest increases in the median mass diameter of the snow, and which have a minor effect on the radar returns or the low frequency brightness temperature, can produce a strong cooling of the 34 GHz brightness temperature.

Meneghini, Robert↗

DOE OSTI2021

Compressed Expanded Natural Graphite (CENG) Processing for PCM Composites

The use of phase change materials (PCMs) in thermal energy storage applications has received considerable attention in recent decades. Organic PCMs are popular due to their high latent heat of fusion, noncorrosive properties, and relative stability over many charge and discharge cycles. A primary limitation of these materials is their low thermal conductivity. This has led researchers to develop various methods to increase thermal conductivity by seeding PCM with or impregnating them into more conductive materials. One method is to impregnate PCMs into compressed expanded natural graphite (CENG) matrices, which can improve thermal conductivity by a factor of 100. CENG matrices have received particular interest due to their low cost, high porosity, small (nano/micro) pore size, high pore density, high thermal conductivity, and ability to be compressed into many geometries. PCM/CENG matrix composites have been extensively studied; however, the effect that CENG processing has on PCM saturation and the overall matrix thermal conductivity has not been well investigated. This processing includes four major steps including graphite intercalation, thermal shock, compression, and PCM saturation. Intercalation involves soaking graphite flakes in sulfuric and/or nitric acid to intercalate the acids between the graphene layers. The graphite flakes are then subjected to a high-temperature thermal shock, during which the intercalated acid is gasified rapidly, pushing the graphene layers apart, resulting in accordion-shaped graphite "worms". The "worms" are then compacted to a desired bulk density and then soaked with molten PCM until fully saturated. The properties of the produced CENG matrix, and its ability to allow PCM permeation, are sensitive to the processing parameters, namely, the thermal shock temperature and exposure time, as well as the matrix apparent density or porosity. Here, we study the effect of the thermal shock conditions necessary to expand intercalated graphite flakes on PCM saturation and the expanded graphite's thermal conductivity and morphology. We found that the thermal shock temperature exhibits the greatest influence. At greater shock temperatures, SEM images showed that expanded graphite worms exhibited greater density of pores, thus increasing total surface area within the matrices. Increasing thermal shock temperature yielded greater overall PCM saturation, as well as an increased rate of saturation. Improvements in PCM saturation rate and overall saturation are obtained as the shock temperature is increased. Longer exposure to thermal shock also improves initial saturation rates and is beneficial if a shortened impregnation time is needed. Thermal shock conditions did not impact thermal conductivity; however, conductivity was largely affected by matrix porosity. A local maximum in axial thermal conductivity was observed at around 83% porosity, which is similar to that observed in previous studies.

compressed expanded natural graphite↗

NASA NTRS2009

Metal Decorated Multi-Walled Carbon Nanotube/Polyimide Composites with High Dielectric Constants and Low Loss Factors

The measurement of observable electromagnetic phenomena in materials and their derived intrinsic electrical material properties are of prime importance in the discovery and development of material systems for electronic and aerospace applications. Nanocomposite materials comprised of metal decorated multi-walled carbon nanotubes (MWCNTs) were prepared by a facile method and characterized. Metal particles such as silver(Ag), platinum(Pt) and palladium(Pd) with diameters ranging from less than 5 to over 50 nanometers were distributed randomly on the MWCNTs. The present study is focused on silver decorated MWCNTs dispersed in a polyimide matrix. The Ag-containing MWCNTs were melt mixed into Ultem(TradeMark) and the mixture extruded as ribbons. The extruded ribbons exhibited a moderate to high degree of MWCNT alignment as determined by HRSEM. These ribbons were then fabricated into test specimens while maintaining MWCNT alignment and subsequently characterized for electrical and electromagnetic properties at 8-12 GHz. The results of the electromagnetic characterization showed that certain sample configurations exhibited a decoupling of the permittivity (epsilon ) and loss factor (epsilon") indicating that these properties could be tailored within certain limits. The decoupling and independent control of these fundamental electrical material parameters offers a new class of materials with potential applications in electronics, microwave engineering and optics.

Elliott, Holly A.↗

DOE OSTI2020

Modelling dislocation-graphene interactions in a BCC Fe matrix by molecular dynamics simulations and gradient plasticity theory

We report graphene nanosheets (GNS) can enhance the strength and ductility of metal-based composites as they can obstruct the propagation of dislocations. The present article employs Molecular Dynamics (MD) simulations to investigate dislocation-GNS interaction mechanisms and possible influencing factors, including the number of GNS layers, the thickness of the metallic amorphous layer and the C - C bond strength. The results indicated that the shear strength of the metal/GNS interface and the bending stiffness of GNS determined the ability of GNS to block dislocation transmission. A physically based phenomenological parameter that can capture such dislocation-GNS interactions is the mechanical interface energy that has been put forth within gradient plasticity. By fitting the theoretical expressions to the simulation data, it was possible to obtain estimates for the mechanical interface energy for the GNS. It was found that increasing the GNS layers and adding an amorphous layer resulted in a strengthening in the stress–strain response and increased the value of this interfacial parameter. This indicates that the mechanical interfacial energy can be a unified measure for capturing and tuning the strength of various interfaces such as grain boundaries, GNS, amorphous-crystalline interface and bimetal interfaces.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Renormalization and matching for the Collins-Soper kernel from lattice QCD

The Collins-Soper kernel, which governs the energy evolution of transverse- momentum dependent parton distribution functions (TMDPDFs), is required to accurately predict Drell-Yan like processes at small transverse momentum, and is a key ingredient for extracting TMDPDFs from experiment. Earlier we proposed a method to calculate this kernel from ratios of the so-called quasi-TMDPDFs determined with lattice QCD, which are defined as hadronic matrix elements of staple-shaped Euclidean Wilson line operators. Here we provide the one-loop renormalization of these operators in a regularization-independent momentum subtraction (RI'/MOM) scheme, as well as the conversion factor from the RI'/MOM-renormalized quasi-TMDPDF to the $ \overline{\mathrm{MS}} $ scheme. We also propose a procedure for calculating the Collins-Soper kernel directly from position space correlators, which simplifies the lattice determination.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2020

Relations between b → cτν decay modes in scalar models

As a consequence of the Ward identity for hadronic matrix elements, we find relations between the differential decay rates of semileptonic decay modes with the underlying quark-level transition b → cτν, which are valid in scalar models. The decay-mode dependent scalar form factor is the only necessary theoretical ingredient for the relations. Otherwise, they combine measurable decay rates as a function of the invariant mass-squared of the lepton pair $q^2$ in such a way that a universal decay-mode independent function is found for decays to vector and pseudoscalar mesons, respectively. This can be applied to the decays $B → D^{\ast}τv, B_s → D^\ast_sτv, B_c → J/ψτv$ and $B → Dτv, B_s → D_sτv, B_c → η_cτv$, with implications for $R(D^{(\ast)}), R(D^{(\ast)}_s), R(J/ψ), R(η_c)$, and $\mathcal{B}(B_c → τv)$. The slope and curvature of the characteristic $q^2$-dependence is proportional to scalar new physics parameters, facilitating their straight forward extraction, complementary to global fits.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

QCD analysis of Δ S = 0 hadronic parity violation

We present a QCD analysis of the effective weak Hamiltonian at hadronic energy scales for strangeness-nonchanging (ΔS = 0) hadronic processes. Performing a leading-order renormalization group analysis in QCD from the W to the O(2 GeV) energy scale, we derive the pertinent effective Hamiltonian for hadronic parity violation, including the effects of both neutral and charged weak currents. We compute the complete renormalization group evolution of all isosectors and the evolution through heavy-flavor thresholds for the first time. We show that the additional four-quark operators that enter below the W mass scale from QCD operator mixing effects form a closed set, and they result in a 12 x 12 anomalous dimension matrix. Computing the resulting effective Hamiltonian and comparing to earlier results, we affirm the importance of operator mixing effects and find, as an example, that the parity-violating pion-nucleon coupling constant, using the factorization Ansatz and an assessment of the pertinent quark charge of the nucleon in lattice QCD at the 2 GeV scale, is in better agreement with recent experiments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2025

Complete definition of 𝑁 →Δ transition generalized parton distributions

We revisit the definition of the leading-twist chiral-even generalized parton distributions (GPDs) for 𝑁 →Δ baryon transitions. We identify and address deficiencies in previous definitions of the transition GPDs inspired by the transition form factors of the vector and axial vector currents. Through systematic analysis of all possible covariant structures, respecting discrete symmetries and the baryon spinor equations of motion, we derive complete sets of independent structures for the transition matrix elements of the vector and axial vector partonic operators. They contain additional structures proportional to the light-cone vector, corresponding to transition GPDs of vanishing first moment, which were not included in previous parametrizations. Their presence is confirmed independently by the light-front multipole expansion and the cross-channel SO(3) partial-wave analysis of the transition matrix elements. Our analysis provides a complete definition of the 𝑁 →Δ transition GPDs for use in theoretical and phenomenological studies.

baryons↗

DOE OSTI2022

Large-Scale Optimization with Linear Equality Constraints Using Reduced Compact Representation

For optimization problems with linear equality constraints, we prove that the (1,1) block of the inverse KKT matrix remains unchanged when projected onto the nullspace of the constraint matrix. In this work, we develop reduced compact representations of the limited-memory inverse BFGS Hessian to compute search directions efficiently when the constraint Jacobian is sparse. Orthogonal projections are implemented by a sparse QR factorization or a preconditioned LSQR iteration. In numerical experiments two proposed trust-region algorithms improve in computation times, often significantly, compared to previous implementations of related algorithms and compared to IPOPT.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1984

A solution of Mx(double dot) + Cx(dot) + Kx = 0 applicable to the design of active dampers

A solution is presented for the equations of motion for the damped linear oscillator, Mx(double dot) + Cx(dot) + Kx = 0. The algorithm solves a transformed set of equations in terms of the modal variables of the undamped system and, at the same time, solves the adjoint equation of the transformed problem. The adjoint solution is normalized to give the inverse of the solution matrix of the transformed problem. The normalized inverse is useful in design for direct computation of sensitivity derivatives of damping ratios with respect to damping rates. The algorithm is programmed to reduce storage requirements by a factor of three-fourths compared to standard complex eigenvalue subroutines. A numerical example is included.

Thurston, G. A.↗

NASA NTRS1984

Human factors of visual displays

Several human factors issues in visual displays are addressed in this report. They are as follows: (1) the importance of luminance range and contrast; (2) uniformity of visual displays; (3) image quality; (4) color contrast; and (5) dot matrix fonts.

Snyder, H. L.↗

NASA NTRS1986

Implicit solution methods in computational fluid dynamics

A class of implicit finite difference schemes of the Beam and Warming approximate factorization type will be addressed. The development and analysis of various aspects of this class of schemes will be given along with the motivations behind many of the choices. Various acceleration and efficiency modifications such as matrix reduction, diagonalization and flux split schemes will be presented. The methods are demonstrated in fully vectorized codes for a CRAY type architecture. The emphasis will be on the Euler equations in generalized coordinates.

Pulliam, Thomas H.↗

NASA NTRS1989

Non linear viscoelasticity applied for the study of durability of polymer matrix composites

A methodology is described for the durability analysis of polymer matrix composites, based on nonlinear viscoelasticity theory. The durability analysis is performed on the basis of a certain number of tests carried out on limited and, if possible, short time scale by the use of accelerating factors. The method was applied to thermomatrix composites under uniaxial and biaxial loadings, showing satisfactory agreement between the life-time predictions and the published data on real-time behavior.

Cardon, A.↗

NASA NTRS2010

Design of Orion Soil Impact Study using the Modern Design of Experiments

Two conventional One Factor At a Time (OFAT) test matrices under consideration for an Orion Landing System subscale soil impact study are reviewed. Certain weaknesses in the designs, systemic to OFAT experiment designs generally, are identified. An alternative test matrix is proposed that is based in the Modern Design of Experiments (MDOE), which achieves certain synergies by combining the original two test matrices into one. The attendant resource savings are quantified and the impact on uncertainty is discussed.

DeLoach, Richard↗

NASA NTRS2018

Polarization Ray Tracing and Polarization Aberration Compensation in Reflective Astronomical Telescopes

Assessing and optimizing polarization performance in the context of ray-based optical design can be challenging. We describe an approach to this problem that decouples polarization effects from optical system geometry for reflective systems. Each surface's polarization properties are parameterized in terms of their impact on retardance and diattenuation in the small angle-of-incidence limit, separating polarization assessment from the task of coating design. A low-resolution ray trace of the system is adequate to determine ray geometry at each interface, which can then be interpolated to rapidly evaluate net Jones Matrix pupil functions. Coating behavior can be easily varied using the ellipsometric parameters to investigate impacts and compensation. Desired values of these parameters can then be specified as constraints in coating design. Investigation with candidate telescope optical designs for LUVOIR show baseline root-mean-square wavefront errors in the nm range for the on-diagonal Jones matrix component, and throughputs of tens of parts per million. Promising possibilities for compensation using a purpose-designed coating on the secondary mirror are discussed, which reduce the on-diagonal wavefront error by a factor 20, with accompanying but more modest reductions in coupling into off-diagonal terms.

Sabatke, Derek↗