Engineering Papers⌕ Search

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.

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

At least 289 records · Page 16

A real time spectrum to dose conversion system

A system has been developed which permits the determination of dose in real time or near real time directly from the pulse-height output of a radiation spectrometer. The technique involves the use of the resolution matrix of a spectrometer, the radiation energy-to-dose conversion function, and the geometrical factors, although the order of matrix operations is reversed. The new technique yields a result which is mathematically identical to the standard method while requiring no matrix manipulations or resolution matrix storage in the remote computer. It utilizes only a single function for each type dose required and each geometric factor involved.

Farmer, B. J.↗

UDU/T/ covariance factorization for Kalman filtering

There has been strong motivation to produce numerically stable formulations of the Kalman filter algorithms because it has long been known that the original discrete-time Kalman formulas are numerically unreliable. Numerical instability can be avoided by propagating certain factors of the estimate error covariance matrix rather than the covariance matrix itself. This paper documents filter algorithms that correspond to the covariance factorization P = UDU(T), where U is a unit upper triangular matrix and D is diagonal. Emphasis is on computational efficiency and numerical stability, since these properties are of key importance in real-time filter applications. The history of square-root and U-D covariance filters is reviewed. Simple examples are given to illustrate the numerical inadequacy of the Kalman covariance filter algorithms; these examples show how factorization techniques can give improved computational reliability.

Thornton, C. L.↗

Multifunctional and biologically active matrices from multicomponent polymeric solutions

A functionalized electrospun matrix for the controlled-release of biologically active agents, such as growth factors, is presented. The functionalized matrix comprises a matrix polymer, a compatibilizing polymer and a biomolecule or other small functioning molecule. In certain aspects the electrospun polymer fibers comprise at least one biologically active molecule functionalized with low molecular weight heparin.

Kiick, Kristi L.↗

An empirical model for polarized and cross-polarized scattering from a vegetation layer

An empirical model for scattering from a vegetation layer above an irregular ground surface is developed in terms of the first-order solution for like-polarized scattering and the second-order solution for cross-polarized scattering. The effects of multiple scattering within the layer and at the surface-volume boundary are compensated by using a correction factor based on the matrix doubling method. The major feature of this model is that all parameters in the model are physical parameters of the vegetation medium. There are no regression parameters. Comparisons of this empirical model with theoretical matrix-doubling method and radar measurements indicate good agreements in polarization, angular trends, and k sub a up to 4, where k is the wave number and a is the disk radius. The computational time is shortened by a factor of 8, relative to the theoretical model calculation.

Liu, H. L.↗

Parathyroid hormone-dependent signaling pathways regulating genes in bone cells

Parathyroid hormone (PTH) is an 84-amino-acid polypeptide hormone functioning as a major mediator of bone remodeling and as an essential regulator of calcium homeostasis. PTH and PTH-related protein (PTHrP) indirectly activate osteoclasts resulting in increased bone resorption. During this process, PTH changes the phenotype of the osteoblast from a cell involved in bone formation to one directing bone resorption. In addition to these catabolic effects, PTH has been demonstrated to be an anabolic factor in skeletal tissue and in vitro. As a result, PTH has potential medical application to the treatment of osteoporosis, since intermittent administration of PTH stimulates bone formation. Activation of osteoblasts by PTH results in expression of genes important for the degradation of the extracellular matrix, production of growth factors, and stimulation and recruitment of osteoclasts. The ability of PTH to drive changes in gene expression is dependent upon activation of transcription factors such as the activator protein-1 family, RUNX2, and cAMP response element binding protein (CREB). Much of the regulation of these processes by PTH is protein kinase A (PKA)-dependent. However, while PKA is linked to many of the changes in gene expression directed by PTH, PKA activation has been shown to inhibit mitogen-activated protein kinase (MAPK) and proliferation of osteoblasts. It is now known that stimulation of MAPK and proliferation by PTH at low concentrations is protein kinase C (PKC)-dependent in both osteoblastic and kidney cells. Furthermore, PTH has been demonstrated to regulate components of the cell cycle. However, whether this regulation requires PKC and/or extracellular signal-regulated kinases or whether PTH is able to stimulate other components of the cell cycle is unknown. It is possible that stimulation of this signaling pathway by PTH mediates a unique pattern of gene expression resulting in proliferation in osteoblastic and kidney cells; however, specific examples of this are still unknown. This review will focus on what is known about PTH-mediated cell signaling, and discuss the established or putative PTH-regulated pattern of gene expression in osteoblastic cells following treatment with catabolic (high) or anabolic (low) concentrations of the hormone.

NASA Discipline Musculoskeletal↗

Matrix toughness, long-term behavior, and damage tolerance of notched graphite fiber-reinforced composite materials

The long-term behavior of notched graphite-fiber-reinforced composite laminates with brittle or tough matrix materials and different fiber architectures was investigated using damage measurements and stiffness change, residual strength, and life data. The fiber/matrix materials included T300/5208, AS4/3501-6, AS4/1808, AS4/PEEK, and C3000/PMR-15 matrices and unidirectional tape and woven cloth fiber architectures. Results of damage evaluation and of residual strength measurements during the fatigue damage development showed that the long-term behavior and damage tolerance are controlled by a number of interacting factors such as the matrix toughness, fiber architecture, loading levels, and damage types and distributions.

Bakis, C. E.↗

An algorithm for the solution of dynamic linear programs

The algorithm's objective is to efficiently solve Dynamic Linear Programs (DLP) by taking advantage of their special staircase structure. This algorithm constitutes a stepping stone to an improved algorithm for solving Dynamic Quadratic Programs, which, in turn, would make the nonlinear programming method of Successive Quadratic Programs more practical for solving trajectory optimization problems. The ultimate goal is to being trajectory optimization solution speeds into the realm of real-time control. The algorithm exploits the staircase nature of the large constraint matrix of the equality-constrained DLPs encountered when solving inequality-constrained DLPs by an active set approach. A numerically-stable, staircase QL factorization of the staircase constraint matrix is carried out starting from its last rows and columns. The resulting recursion is like the time-varying Riccati equation from multi-stage LQR theory. The resulting factorization increases the efficiency of all of the typical LP solution operations over that of a dense matrix LP code. At the same time numerical stability is ensured. The algorithm also takes advantage of dynamic programming ideas about the cost-to-go by relaxing active pseudo constraints in a backwards sweeping process. This further decreases the cost per update of the LP rank-1 updating procedure, although it may result in more changes of the active set that if pseudo constraints were relaxed in a non-stagewise fashion. The usual stability of closed-loop Linear/Quadratic optimally-controlled systems, if it carries over to strictly linear cost functions, implies that the saving due to reduced factor update effort may outweigh the cost of an increased number of updates. An aerospace example is presented in which a ground-to-ground rocket's distance is maximized. This example demonstrates the applicability of this class of algorithms to aerospace guidance. It also sheds light on the efficacy of the proposed pseudo constraint relaxation scheme.

Psiaki, Mark L.↗

The optimal fiber volume fraction and fiber-matrix property compatibility in fiber reinforced composites

Although the question of minimum or critical fiber volume fraction beyond which a composite can then be strengthened due to addition of fibers has been dealt with by several investigators for both continuous and short fiber composites, a study of maximum or optimal fiber volume fraction at which the composite reaches its highest strength has not been reported yet. The present analysis has investigated this issue for short fiber case based on the well-known shear lag (the elastic stress transfer) theory as the first step. Using the relationships obtained, the minimum spacing between fibers is determined upon which the maximum fiber volume fraction can be calculated, depending on the fiber packing forms within the composites. The effects on the value of this maximum fiber volume fraction due to such factors as fiber and matrix properties, fiber aspect ratio and fiber packing forms are discussed. Furthermore, combined with the previous analysis on the minimum fiber volume fraction, this maximum fiber volume fraction can be used to examine the property compatibility of fiber and matrix in forming a composite. This is deemed to be useful for composite design. Finally some examples are provided to illustrate the results.

Pan, Ning↗

Regulation of expression of collagenase-3 in normal, differentiating rat osteoblasts

We investigated the regulation of collagenase-3 expression in normal, differentiating rat osteoblasts. Fetal rat calvarial cell cultures showed an increase in alkaline phosphatase activity reaching maximal levels between 7-14 days post-confluence, then declining with the onset of mineralization. Collagenase-3 mRNA was just detectable after proliferation ceased at day 7, increased up to day 21, and declined at later ages. Postconfluent cells maintained in non-mineralizing medium expressed collagenase-3 but did not show the developmental increase exhibited by cells switched to mineralization medium. Cells maintained in non-mineralizing medium continued to proliferate; cells in mineralization medium ceased proliferation. In addition, collagenase-3 mRNA was not detected in subcultured cells allowed to remineralize. These results suggest that enhanced accumulation of collagenase-3 mRNA is triggered by cessation of proliferation or acquisition of a mineralized extracellular matrix and that other factors may also be required. After initiation of basal expression, parathyroid hormone (PTH) caused a dose-dependent increase in collagenase-3 mRNA. Both the cyclic adenosine monophosphate (cAMP) analogue, 8-bromo-cAMP (8-Br-cAMP), and the protein kinase C (PKC) activator, phorbol myristate acetate, increased collagenase-3 expression, while the calcium ionophore, ionomycin, did not, suggesting that PTH was acting through the protein kinase A (PKA) and PKC pathways. Inhibition of protein synthesis with cycloheximide caused an increase in basal collagenase-3 expression but blocked the effect of PTH, suggesting that an inhibitory factor prevents basal expression while an inductive factor is involved with PTH action. In summary, collagenase-3 is expressed in mineralized osteoblasts and cessation of proliferation and initiation of mineralization are triggers for collagenase-3 expression. PTH also stimulates expression of the enzyme through both PKA and PKC pathways in the mineralizing osteoblast. Copyright 1999 Wiley-Liss, Inc.

Non-NASA Center↗

Detailed analysis of excited-state systematics in a lattice QCD calculation of 𝑔 𝐴

Excited state contamination remains one of the most challenging sources of systematic uncertainty to control in lattice QCD calculations of nucleon matrix elements and form factors: early time separations are contaminated by excited states and late times suffer from an exponentially bad signal-to-noise problem. High-statistics calculations at large time separations ≳ 1 fm are commonly used to combat these issues. In this work, focusing on g A , we explore the alternative strategy of utilizing a large number of relatively low-statistics calculations at short to medium time separations (0.2–1 fm), combined with a multistate analysis. On an ensemble with a pion mass of approximately 310 MeV and a lattice spacing of approximately 0.09 fm, we find this provides a more robust and economical method of quantifying and controlling the excited state systematic uncertainty. A quantitative separation of various types of excited states enables the identification of the transition matrix elements as the dominant contamination. The excited state contamination of the Feynman-Hellmann correlation function is found to reduce to the 1% level at approximately 1 fm while, for the more standard three-point functions, this does not occur until after 2 fm. Critical to our findings is the use of a global minimization, rather than fixing the spectrum from the two-point functions and using them as input to the three-point analysis. We find that the ground state parameters determined in such a global analysis are stable against variations in the excited state model, the number of excited states, and the truncation of early-time or late-time numerical data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

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

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

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

General Bayesian Inference over the Stiefel Manifold via the Givens Transform

We introduce here an approach based on the Givens representation for posterior inference in statistical models with orthogonal matrix parameters, such as factor models and probabilistic principal component analysis (PPCA). We show how the Givens representation can be used to develop practical methods for transforming densities over the Stiefel manifold into densities over subsets of Euclidean space. We demonstrate how to deal with issues arising from the topology of the Stiefel manifold and how to inexpensively compute the change-of-measure terms. We introduce an auxiliary parameter approach that limits the impact of topological issues. We provide both analysis of our methods and numerical examples demonstrating the effectiveness of the approach. We also discuss how our Givens representation can be used to define general classes of distributions over the space of orthogonal matrices. We then give demonstrations on several examples showing how the Givens approach performs in practice in comparison with other methods.

97 MATHEMATICS AND COMPUTING↗

The continuum and leading twist limits of parton distribution functions in lattice QCD

In this study, we present continuum limit results for the unpolarized parton distribution function of the nucleon computed in lattice QCD. This study is the first continuum limit using the pseudo-PDF approach with Short Distance Factorization for factorizing lattice QCD calculable matrix elements. Our findings are also compared with the pertinent phenomenological determinations. Inter alia, we are employing the summation Generalized Eigenvalue Problem (sGEVP) technique in order to optimize our control over the excited state contamination which can be one of the most serious systematic errors in this type of calculations. A crucial novel ingredient of our analysis is the parameterization of systematic errors using Jacobi polynomials to characterize and remove both lattice spacing and higher twist contaminations, as well as the leading twist distribution. This method can be expanded in further studies to remove all other systematic errors.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Comprehensive review of 2 β decay half-lives

Here, the double-beta (2 β )-decay is the rarest nuclear physics process, and its experimental half-lives (T 1/2 ) exceed the age of the Universe from nine to fourteen orders of magnitude. Double-beta decay was observed, and its half-life was measured in 14 parent nuclei using direct, radiochemical, and geochemical methods. The decay observables are analyzed using the Evaluated Nuclear Structure Data File (ENSDF) procedures, and the recommended T 1/2 were deduced. Using the calculated values of phase factors, the effective nuclear matrix elements were extracted and compared with available data. Thousands of theoretical and experimental works have been dedicated to these topics in the last 85 years, and we present two data sets of recommended values to encapsulate the results.

2β-decay↗

Skeletal model reduction with forced optimally time dependent modes

Skeletal model reduction based on local sensitivity analysis of time dependent systems is presented in which sensitivities are modeled by forced optimally time dependent (f-OTD) modes. The f-OTD factorizes the sensitivity coefficient matrix into a compressed format as the product of two skinny matrices, i.e. f-OTD modes and f-OTD coefficients. The modes create a low-dimensional, time dependent, orthonormal basis which capture the directions of the phase space associated with most dominant sensitivities. These directions highlight the instantaneous active species, and reaction paths. Evolution equations for the f-OTD modes and coefficients are derived, and the implementation of f-OTD for skeletal reduction is described. For demonstration, skeletal reduction is conducted of the constant pressure ethylene-air burning in a zero-dimensional reactor, and new reduced models are generated. The laminar flame speed, the ignition delay, and the extinction curve as predicted by the models are compared against some existing skeletal models in literature for the same detailed model. Finally, the results demonstrate the capability of f-OTD to eliminate unimportant reactions and species in a systematic, efficient and accurate manner.

42 ENGINEERING↗

Tailoring grain growth and solid solution strengthening of single-phase CrCoNi medium-entropy alloys by solute selection

In the present study, we selected solutes to be added to the CrCoNi medium-entropy alloy (MEA) based on the mismatch of self-diffusion activation energy (SDQ) between the alloying elements and constituent elements of the matrix, and then investigated their grain growth behavior and mechanical properties. Mo and Al were selected as the solutes for investigation primarily because they have higher and lower SDQ, respectively, than those of the matrix elements; a secondary factor was their higher and lower shear modulus. Their concentrations were fixed at 3 at.% each because previous work had shown these compositions to be single-phase solid solutions with the face-centered cubic structure. Three alloys were produced by arc melting, casting, homogenizing, cold rolling and annealing at various temperatures and times to produce samples with different grain sizes. They were (a) the base alloy CrCoNi, (b) the base alloy plus 3 at.% Mo, and (c) the base alloy plus 3 at.% Al. The activation energies for grain growth of the CrCoNi, CrCoNi-3Mo and CrCoNi-3Al MEAs were found to be ~251, ~368 and ~219 kJ/mol, respectively, consistent with the notion that elements with higher SDQ (in this study Mo) retard grain growth (likely by a solute-drag effect), whereas those with lower values (Al) accelerate grain growth. The room-temperature tensile properties show that Mo increases the yield strength by ~40% but Al addition has a smaller strengthening effect consistent with their relative shear moduli. The yield strength as a function of grain size for the three singlephase MEAs follows the classical Hall-Petch relationship with much higher slopes (> 600 MPa μm -0.5 ) than traditional solid solutions. Furthermore, this work shows that the grain growth kinetics and solid solution strengthening of the CrCoNi MEA can be tuned by selecting solute elements that have appropriate diffusion and physical properties.

36 MATERIALS SCIENCE↗

Operator-level quantum acceleration of non-logconcave sampling

Sampling from probability distributions of the form 𝝈 ∝ e −𝜷V , where V is a continuous potential, is a fundamental task across physics, chemistry, biology, computer science, and statistics. However, when V is nonconvex, the resulting distribution becomes non-logconcave, and classical methods such as Langevin dynamics often exhibit poor performance. We introduce a quantum algorithm that provably accelerates a broad class of continuous-time sampling dynamics. For Langevin dynamics, our method encodes the target Gibbs measure into the amplitudes of aquantum state, identified as the kernel of a block matrix derived from a factorization of the Witten Laplacian operator. This connection enables Gibbs sampling via singular value thresholding and yields up to a quartic quantum speedup over best-knownclassical Langevin-based methods in the non-logconcave setting. Building on this framework, we further develop the first quantum algorithm that accelerates replica exchange Langevin diffusion, a widely used method for sampling from complex, rugged energy landscapes.

97 MATHEMATICS AND COMPUTING↗

Overview of Potential Hazards in Electric Aircraft Charging Infrastructure

With increasing efforts in electrification of Advanced Air Mobility and electric aircraft, there is a growing need to install new infrastructure to support them. This report highlights potential hazards (non exhaustive) that the installation must be aware of and be prepared to mitigate with increased electrical equipment on site. The hazards can be due to natural causes as well as non natural causes. The report aims to provide starting guidelines for awareness of these hazards and to help plan for mitigation of the same.

24 POWER TRANSMISSION AND DISTRIBUTION↗