Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “matrix approximation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,009 records · Page 56

Braided Composite Threaded Fasteners

Fasteners of carbon/carbon or other fiber/matrix composite materials made by braiding-and-molding process. Preform of braided fibers impregnated with matrix material pressed and cured in internally threaded, split mold to make externally threaded rod. Process costs less, produces stronger fasteners by avoiding breakage of fibers, and adaptable to mass production. Intended for use with ceramic and composite-material structural parts at temperatures approximately greater than 2,500 degrees F (approximately greater than 1,400 degrees C), at which metal fasteners cannot be used.

Sawyer, James Wayne↗

Effect of grid aspect ratio on convergence

The effects of high-aspect-ratio grids on algorithm convergence are considered by means of vector stability theory and computational experiments. The results indicate that approximately factored implicit schemes experience convergence deterioration because of nonoptimum local time-stepping procedures and increased need for viscous preconditioning. Based on this insight, an enhanced algorithm is devised using improved selection of the local time step, appropriate definition of the viscous preconditioning matrix, and proper implementation of the boundary conditions. The new algorithm provides uniformly efficient convergence at all aspect ratios for both Euler and Navier-Stokes computations for a variety of test problems.

Buelow, P. E. O.↗

On substructuring algorithms and solution techniques for the numerical approximation of partial differential equations

Substructuring methods are in common use in mechanics problems where typically the associated linear systems of algebraic equations are positive definite. Here these methods are extended to problems which lead to nonpositive definite, nonsymmetric matrices. The extension is based on an algorithm which carries out the block Gauss elimination procedure without the need for interchanges even when a pivot matrix is singular. Examples are provided wherein the method is used in connection with finite element solutions of the stationary Stokes equations and the Helmholtz equation, and dual methods for second-order elliptic equations.

Gunzburger, M. D.↗

Binary Cation Matrix Electrolyte and Its Effect on Solid Electrolyte Interphase Suppression and Evolution of Si Anode

An unstable solid electrolyte interphase (SEI) has been recognized as one of the biggest challenges to commercializing silicon (Si) anodes for high-energy-density batteries. This work thoroughly investigates a binary cation matrix of Mg 2+ +Li + electrolyte and its role in SEI development, suppression, and evolution of a Si anode. Findings demonstrate that introducing Mg ions dramatically reduces the SEI growth before lithiation occurs, primarily due to the suppression of solvent reduction, particularly ethylene carbonate (EC) reduction. The Mg 2+ alters the Li + cation solvation environment as EC preferably participates in the oxophyllic Mg 2+ solvation sheath, thereby altering the solvent reduction process, resulting in a distinct SEI formation mechanism. The initial SEI formation before lithiation is reduced by 70% in the electrolyte with the presence of Mg 2+ cations. While the SEI continues to develop in the postlithiation, the inclusion of Mg ions results in an approximately 80% reduction in the postlithiation SEI growth. Continuous electrochemical cycling reveals that Mg 2+ plays a crucial role in stabilizing the deep-lithiated Si phases, which effectively mitigates side reactions, resulting in controlled SEI growth and stable interphase while eliminating complex Li x Si y formation. Mg ions promote the development of a notably more rigid and homogeneous SEI, characterized by a reduced dissipation (ΔD) in the Mg 2+ +Li + ion matrix compared to the solely Li + system. In conclusion, this report reveals how the Mg 2+ +Li + ion matrix affects the SEI evolution, viscoelastic properties, and electrochemical behavior at the Si interface in real time, laying the groundwork for devising strategies to enhance the performance and longevity of Si-based next-generation battery systems.

25 ENERGY STORAGE↗

Benchmark calculations of electron impact electronic excitation of the hydrogen molecule

We present benchmark integrated and differential cross-sections for electron collisions with H 2 using two different theoretical approaches, namely, the R-matrix and molecular convergent close-coupling. This is similar to comparative studies conducted on electron–atom collisions for H, He and Mg. Electron impact excitation to the ${\mathrm{b}}^{3}{{\Sigma}}_{\mathrm{u}}^{+}$, ${\mathrm{a}}^{3}{{\Sigma}}_{\mathrm{g}}^{+}$, ${\mathrm{B}}^{1}{{\Sigma}}_{\mathrm{u}}^{+}$, c3Πu, ${\mathrm{E}\mathrm{F}}^{1}{{\Sigma}}_{\mathrm{g}}^{+}$, C1Πu, ${\mathrm{e}}^{3}{{\Sigma}}_{\mathrm{u}}^{+}$, ${\mathrm{h}}^{3}{{\Sigma}}_{\mathrm{g}}^{+}$, ${\mathrm{B}}^{\prime 1}{{\Sigma}}_{\mathrm{u}}^{+}$ and d3Πu excited electronic states are considered. Calculations are presented in both the fixed nuclei and adiabatic nuclei approximations, where the latter is shown only for the ${\mathrm{b}}^{3}{{\Sigma}}_{\mathrm{u}}^{+}$ state. Good agreement is found for all transitions presented. Where available, we compare with existing experimental and recommended data.

74 ATOMIC AND MOLECULAR PHYSICS↗

Mean-field dynamo as a quantum-like modulational instability

Presented here is a novel formulation of the mean-field dynamo as a modulational instability of magnetohydrodynamic (MHD) turbulence. This formulation, termed mean-field wave kinetics (MFWK), is based on the Weyl symbol calculus and allows describing the interaction between the mean fields (magnetic field and fluid velocity) and turbulence without requiring scale separation that is commonly assumed in the literature. The turbulence is described by the Wigner–Moyal equation for the spectrum of the two-point correlation matrix (Wigner matrix) of magnetic-field and velocity fluctuations and depicts the turbulence as an effective plasma of quantum-like particles that interact via the mean fields. Eddy–eddy interactions, which serve as ‘collisions’ in this effective plasma, are modelled within the standard minimal tau approximation to aid comparison with existing theories. Using MFWK, the non-local electromotive force is calculated for generic turbulence from first principles, modulo the limitations of MFWK. This result is then used to study, both analytically and numerically, the modulational modes of MHD turbulence, which appear as linear instabilities of the said effective quantum-like plasma of fluctuations. The standard α 2 -dynamo and other known results are reproduced as special cases. A new dynamo effect is predicted that is driven by correlations between the turbulent flow velocity and the turbulent current.

astrophysical plasmas↗

Glycolate Analysis in Tank 22: Developing and Testing Analytical Methods for the Savannah River Site Liquid Waste System

Researchers developed and tested a range of techniques to support a low mg/L Limit of Quantitation (LOQ) of glycolate in radioactive Tank 22 waste solution by Ion Chromatography (IC) and by proton NMR. For the IC method, Dionex OnGuard II cartridges were tested as a means of removing alkali earth and transition metals that can interfere during ion chromatography applications especially when analyzing for low-concentration, chelating analytes. Additionally, concentrations of nitrate in the raw Tank 22 sample (5,000 to 10,000 mg/L) were managed by using reasonable levels of sample dilution. The resulting IC performance quality was enhanced by improving the baseline, increasing sensitivity [Limit of Quantitation (LOQ) of 12 mg/L and Limit of Detection (LOD) of 4.0 mg/L], and resolving analytes into well-defined Gaussian peaks when using the Dionex OnGuard II H⁺ cartridges to remove matrix interferences. High concentrations of nitrate limited the performance of the IC method. To achieve an acceptable baseline, samples with high nitrate content require more dilution, resulting in higher detection limits. Tests on samples with higher concentrations of nitrate (Tank 30 and 32 supernate with approximately 150000 mg/L nitrate) suggested that IC detection limits for glycolate in these samples would be > 500 mg/L. Thus, an LOQ of 12 mg/L is not feasible by this IC method on evaporator feed samples. An alternative method of glycolate analysis, using proton nuclear magnetic resonance (H NMR), was developed by the research team. In initial tests, the H NMR technique provided reasonable quantitation of glycolate in Tank 22 conditions by direct observation of the liquid. The H NMR method may provide improved detection limits for solutions with higher nitrate concentrations. We recommend further development of this analysis for high nitrate LWS samples such as evaporator feed and evaporator drop tank content. A method for pretreatment of samples using crystalline silicotitanate (CST) was developed and tested. The objective of the pretreatment was to facilitate analysis of samples with higher levels of radioactivity. The pretreatment did not influence the glycolate concentration in solution, and the pretreatment is projected to reduce Cs 137 activity in a sample by a factor of 16,200. The double strike CST pretreatment was demonstrated and would allow milliliters of higher activity samples to be transferred from the Shielded Cells and handled in a containment unit for glycolate analysis. The various studies validated IC and H NMR methods for glycolate analysis, defined the range of applicability, and demonstrated key supporting analytical protocols. Based on the results, high quality glycolate analysis of Tank 22 is feasible down to approximately 12 mg/L, with the potential for broader applicability of the methods to other conditions in the Savannah River Site Liquid Waste System (LWS).

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

A Simulation-based Method for Correcting Mode Coupling in CMB Angular Power Spectra

Modern cosmic microwave background (CMB) analysis pipelines regularly employ complex time-domain filters, beam models, masking, and other techniques during the production of sky maps and their corresponding angular power spectra. However, these processes can generate couplings between multipoles from the same spectrum and from different spectra, in addition to the typical power attenuation. Within the context of pseudo-C ℓ based, MASTER-style analyses, the net effect of the time-domain filtering is commonly approximated by a multiplicative transfer function, F ℓ , that can fail to capture mode mixing and is dependent on the spectrum of the signal. To address these shortcomings, we have developed a simulation-based spectral correction approach that constructs a two-dimensional transfer matrix, ${J}_{{\ell }{\ell }^{\prime} }$, which contains information about mode mixing in addition to mode attenuation. We demonstrate the application of this approach on data from the first flight of the Spider balloon-borne CMB experiment.

79 ASTRONOMY AND ASTROPHYSICS↗

Polynominal Interpolation Methods for Viscous Flow Calculations

Higher-order collocation procedures resulting in tridiagonal matrix systems are derived from polynomial spline interpolation and by Hermitian (Taylor series) finite-difference discretization. The similarities and special features of these different developments are discussed. The governing systems apply for both uniform and variable meshes. Hybrid schemes resulting from two different polynomial approximations for the first and second derivatives lead to a nonuniform mesh extension of the so-called compact or Pad? difference technique (Hermite 4). A variety of fourth-order methods are described and the Hermitian approach is extended to sixth-order (Hermite 6). The appropriate spline boundary conditions are derived for all procedures. For central finite differences, this leads to a two-point, second-order accurate generalization of the commonly used three-point end-difference formula. Solutions with several spline and Hermite procedures are presented for the boundary layer equations, with and without mass transfer, and for the incompressible viscous flow in a driven cavity. Divergence and nondivergence equations are considered for the cavity. Among the fourth-order techniques, it is shown that spline 4 has the smallest truncation error. The spline 4 procedure generally requires one-quarter the number of mesh points in a given coordinate direction as a central finite-difference calculation of equal accuracy. The Hermite 6 procedure leads to remarkably accurate boundary layer solutions.

Rubin, S. G.↗

Recursive inversion of externally defined linear systems

The approximate inversion of an internally unknown linear system, given by its impulse response sequence, by an inverse system having a finite impulse response, is considered. The recursive least squares procedure is shown to have an exact initialization, based on the triangular Toeplitz structure of the matrix involved. The proposed approach also suggests solutions to the problems of system identification and compensation.

Bach, Ralph E., Jr.↗

An improved error assessment for the GEM-T1 gravitational model

Several tests were designed to determine the correct error variances for the Goddard Earth Model (GEM)-T1 gravitational solution which was derived exclusively from satellite tracking data. The basic method employs both wholly independent and dependent subset data solutions and produces a full field coefficient estimate of the model uncertainties. The GEM-T1 errors were further analyzed using a method based upon eigenvalue-eigenvector analysis which calibrates the entire covariance matrix. Dependent satellite and independent altimetric and surface gravity data sets, as well as independent satellite deep resonance information, confirm essentially the same error assessment. These calibrations (utilizing each of the major data subsets within the solution) yield very stable calibration factors which vary by approximately 10 percent over the range of tests employed. Measurements of gravity anomalies obtained from altimetry were also used directly as observations to show that GEM-T1 is calibrated. The mathematical representation of the covariance error in the presence of unmodeled systematic error effects in the data is analyzed and an optimum weighting technique is developed for these conditions. This technique yields an internal self-calibration of the error model, a process which GEM-T1 is shown to approximate.

Lerch, F. J.↗

Application of a two-point exponential approximation method in optimizing rotorcraft airframe structures

The Fadel (1990) two-point exponential approximation technique for rotorcraft airframe structural optimization is presently examined by comparing its results with the well established linear and hybrid first-order methods. Attention is given to both the elastic-line and FEM models of an AH-1G airframe. The Direct Matrix Abstraction Programming Language is used to compute sensitivity derivatives and motivate the CONMIN optimizer to produce the design change.

Sareen, Ashish K.↗

The effect of grid aspect ratio on convergence

Vector stability theory and computational experiments are used to examine the effects of high aspect ratio grids on algorithm convergence. Emphasis is placed on identifying the degree of convergence deterioration that ADI schemes experience as grid aspect ratio increases and on deriving methods for circumventing this problem. Approximately factored implicit schemes are found to experience convergence deterioration because of nonoptimum local time stepping procedures and the increased need for viscous preconditioning. An enhanced algorithm is devised on the basis of improved selection of the local time step, appropriate definition of viscous preconditioning matrix, and proper implementation of boundary conditions. The novel algorithm provides uniformly efficient convergence at all aspect ratios for both Euler and Navier-Stokes computations for a variety of test problems.

Buelow, P. E. O.↗

High-Performance SiC/SiC Ceramic Composite Systems Developed for 1315 C (2400 F) Engine Components

As structural materials for hot-section components in advanced aerospace and land-based gas turbine engines, silicon carbide (SiC) ceramic matrix composites reinforced by high performance SiC fibers offer a variety of performance advantages over current bill-of-materials, such as nickel-based superalloys. These advantages are based on the SiC/SiC composites displaying higher temperature capability for a given structural load, lower density (approximately 30- to 50-percent metal density), and lower thermal expansion. These properties should, in turn, result in many important engine benefits, such as reduced component cooling air requirements, simpler component design, reduced support structure weight, improved fuel efficiency, reduced emissions, higher blade frequencies, reduced blade clearances, and higher thrust. Under the NASA Ultra-Efficient Engine Technology (UEET) Project, much progress has been made at the NASA Glenn Research Center in identifying and optimizing two highperformance SiC/SiC composite systems. The table compares typical properties of oxide/oxide panels and SiC/SiC panels formed by the random stacking of balanced 0 degrees/90 degrees fabric pieces reinforced by the indicated fiber types. The Glenn SiC/SiC systems A and B (shaded area of the table) were reinforced by the Sylramic-iBN SiC fiber, which was produced at Glenn by thermal treatment of the commercial Sylramic SiC fiber (Dow Corning, Midland, MI; ref. 2). The treatment process (1) removes boron from the Sylramic fiber, thereby improving fiber creep, rupture, and oxidation resistance and (2) allows the boron to react with nitrogen to form a thin in situ grown BN coating on the fiber surface, thereby providing an oxidation-resistant buffer layer between contacting fibers in the fabric and the final composite. The fabric stacks for all SiC/SiC panels were provided to GE Power Systems Composites for chemical vapor infiltration of Glenn designed BN fiber coatings and conventional SiC matrices. Composite panels with system B were heat treated at Glenn, and the pores that remained open were filled by silicon melt infiltration (MI). Panels with system A and the other SiC/SiC systems were not heat treated, and remaining open pores in these systems were filled with SiC slurry and silicon MI.

DiCarlo, James A.↗

Tunable Fiber Bragg Grating Ring Lasers using Macro Fiber Composite Actuators

The research reported herein includes the fabrication of a tunable optical fiber Bragg grating (FBG) fiber ring laser (FRL)1 from commercially available components as a high-speed alternative tunable laser source for NASA Langley s optical frequency domain reflectometer (OFDR) interrogator, which reads low reflectivity FBG sensors. A Macro-Fiber Composite (MFC) actuator invented at NASA Langley Research Center (LaRC) was selected to tune the laser. MFC actuators use a piezoelectric sheet cut into uniaxially aligned rectangular piezo-fibers surrounded by a polymer matrix and incorporate interdigitated electrodes to deliver electric fields along the length of the piezo-fibers. This configuration enables MFC actuators to produce displacements larger than the original uncut piezoelectric sheet. The FBG filter was sandwiched between two MFC actuators, and when strained, produced approximately 3.62 nm of wavelength shift in the FRL when biasing the MFC actuators from 500 V to 2000 V. This tunability range is comparable to that of other tunable lasers and is adequate for interrogating FBG sensors using OFDR technology. Three different FRL configurations were studied. Configuration A examined the importance of erbium-doped fiber length and output coupling. Configuration B demonstrated the importance of the FBG filter. Configuration C added an output coupler to increase the output power and to isolate the filter. Only configuration C was tuned because it offered the best optical power output of the three configurations. Use of Plastic Optical Fiber (POF) FBG s holds promise for enhanced tunability in future research.

Geddis, Demetris L.↗

Post-translational control of collagen fibrillogenesis in mineralizing cultures of chick osteoblasts

Cultured osteoblasts from chick embryo calvaria were used as a model system to investigate the post-translational extracellular mechanisms controlling the macroassembly of collagen fibrils. The results of these studies demonstrated that cultured osteoblasts secreted a collagenous extracellular matrix that assembled and mineralized in a defined temporal and spatial sequence. The assembly of collagen occurred in a polarized fashion, such that successive orthogonal arrays of fibrils formed between successive cell layers proceeding from the culture surface toward the media. Mineralization followed in the same manner, being observed first in the deepest and oldest fibril layers. Collagen fibrillogenesis, the kinetics of cross-link formation, and collagen stability in the extracellular matrix of the cultures were examined over a 30 day culture period. Between days 8 and 12 in culture, collagen fibril diameters increased from < 30 nm to an average of 30-45 nm. Thereafter, diameters ranged in size from 20 to 200 nm. Quantitation of the collagen cross-linking residues, hydroxylysyl pyridinoline (HP) and lysyl pyridinoline (LP), showed that these mature cross-links increased from undetectable levels to concentrations found in normal chick bone. Analysis of the kinetics of their formation by pulse-chase labeling the cultures with [3H]lysine showed a doubling time of approximately 5 days. The relationships between cross-link formation, fibrillogenesis, and collagen stability were examined in cultures treated with beta-aminopropionitrile (beta-APN), a potent inhibitor of lysyl oxidase and cross-link formation. In beta-APN-treated cultures, total collagen synthesis was increased twofold, with no change in mRNA levels for type I collagen, whereas the amount of collagen accumulated in the cell layer was decreased by 50% and mineral deposition was reduced. The rate of collagen retention in the matrix was assessed by pulse-chase analysis of [3H]proline over a 16 day period in control and beta-APN-treated cultures. In control cultures, about 20% of the labeled collagen was lost from the cell layers over a 16 day period compared with > 80% in the presence of beta-APN. The beta-APN-treated cultures also showed a wider diversity of fibril diameters with a median in the > 45-60 nm range. In summary, these data suggest that cross-linking and assembly of collagen fibrils secreted by osteoblasts in vitro occur in a fashion similar to that found in vivo. The rate of cross-link formation is relatively constant and may be correlated with increasing collagen mass.(ABSTRACT TRUNCATED AT 400 WORDS).

NASA Discipline Musculoskeletal↗

Exploration with Scalable Gaussian Process Reinforcement Learning

Exploration is a challenging problem in reinforcement learning (RL), especially in environments with sparse rewards. Quantifying and utilizing the parametric uncertainty has been shown to be paramount for successful exploration [Osband et al., 2018]. Bayesian, or approximately Bayesian, methods present a principled means of estimating the parametric uncertainty in RL problems. Gaussian processes, nonparametric Bayesian models, are often impractical due to poor scalability and computational bottlenecks. We introduce a scalable Gaussian process RL (GPRL) method which directly induces sparsity in the covariance matrix to facilitate faster computation. This is a departure from previous GPRL methods which instead rely on data reduction and subsampling. We compare various covariance-based exploration techniques (Thompson sampling, upper confidence bound, and probabilistic maximum variance) which leverage our scalable GP framework in sparse reward environments. Finally, we show favorable comparison against the bootstrapped deep Q-Network.

97 MATHEMATICS AND COMPUTING↗

On the stability of Galerkin methods for initial-boundary value problems for hyperbolic systems

The stability of approximating the solution of mixed initial-boundary value problems for hyperbolic systems by semidiscrete Galerkin methods is studied. It is shown that a particular straightforward Galerkin method yields an unstable approximation, and that this numerical instability is caused by an improper treatment of the boundary. Stable schemes are then presented, one of which differs from the unstable scheme only insofar as the treatment of the boundary is concerned. These stable schemes make use of a particular matrix which symmetrizes the differential system. It is therefore shown that the use of this matrix is crucial to the stability of the computations as well as for obtaining a priori bounds on the energy of the continuous system. This symmetrizing matrix is also related to the diagonalizing matrix for the system of hyperbolic equations and to the Liapunov matrix for the system of ordinary differential equations resulting from the application of Galerkin's method.

Gunzburger, M. D.↗