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,189 records · Page 66

Carbon Nanotube/Polymer Nanocomposites Flexible Stress and Strain Sensors

Conformable stress and strain sensors are required for monitoring the integrity of airframe structures as well as for sensing the mechanical stimuli in prosthetic arms. For this purpose, we have developed a series of piezoresistive single-wall carbon nanotube (SWCNT)/polymer nanocomposites. The electromechanical coupling of pressure with resistance changes in these nanocomposites is exceptionally greater than that of metallic piezoresistive materials. In fact, the piezoresistive stress coefficient (pi) of a SWCNT/polymer nanocomposite is approximately two orders of magnitude higher than that of a typical metallic piezoresistive. The piezoresistive stress coefficient is a function of the nanotube concentration wherein the maximum value occurs at a concentration just above the percolation threshold concentration (phi approx. 0.05 %). This response appears to originate from a change in intrinsic resistivity under compression/tension. A systematic study of the effect of the modulus of the polymer matrix on piezoresistivity allowed us to make flexible and conformable sensors for biomedical applications. The prototype haptic sensors using these nanocomposites are demonstrated. The piezocapacitive properties of SWCNT/polymer are also characterized by monitoring the capacitance change under pressure.

Kang, Jin Ho↗

Reduction of Raman Spectroscopy Data for H2-CO2-Air Tubular Flame Measurements

An approach of reducing spontaneous Raman spectroscopy data, based on the matrix inversion method which logically removes undesired sources of the Raman scattering signals including cross-talk due to the signals from other species, is discussed. The approach contains extensive calibration of temperature-dependent system response. The Raman spectra of several cold and hot calibration flames and their resulting calibration matrix are elaborated. The structure of H2-CO2-air tubular flames, in terms of major species (CO2, O2, N2, H2O, and H2) concentrations and temperature, was successfully measured with the optical technique. The Raman spectra data used in this study are taken from a uniquely designed experimental apparatus – a tubular burner and Raman spectroscopy system at Vanderbilt University. The characteristic uncertainties for chemical species are approximately ±2% by mole fraction for hot products and ±0.5% for room temperature reactants. The approach described here for hydrogen/air tubular flames may be adapted to hydrocarbon/air tubular flames when other necessary major species, for example, CH4 and CO for methane/air tubular flame, are considered in the whole process.

Raman Spectroscopy↗

Simulation and Experimental Validation of an Integrated Heat Pump – Thermal Energy Storage Using a Room-Temperature Phase Change Material

As the dependence on electrical heat pumps (HPs) and intermittent renewables increases, grid strains are expected to grow. This necessitates an energy storage system to reduce the mismatch between energy supply and demand. Thus, a proposed dual-mode commercially available 14.1 kW HP was integrated with a single 22°C phase change material (PCM) thermal storage system (TES) to load-shift both cooling and heating loads. The HP-TES system was manufactured and experimentally tested using a novel test matrix based on AHRI 210/240 psychrometric conditions. Furthermore, transient dual-mode system-level HP-TES models were developed in Modelica and validated using the experimental test conditions. Base HP cooling and heating experimental tests at ambient temperatures of 35°C and −8.3°C show that the modified HP-TES maintained the rated system capacity and performance. The HP-TES discharge provided approximately 30% and 50% reductions in cooling and heating demand, respectively. The transient HP-TES models predicted system capacity and total power input for discharge and recharge operating modes within ±4% mean percentage error, and recharge power input within ±2%, with maximum errors occurring at the equipment startup. During system operation, the sources of model deviations are first-order polynomial fits of the PCM digital scanning calorimetry (DSC) data and unaccounted supercooling in the PCM during solidification. Nonetheless, the model predictions agree with the experimental tests, demonstrating the availability of robust, accurate, and validated transient models that can be used for further validation and the development of system controls.

25 ENERGY STORAGE↗

Development of Novel Ferritic-Martensitic Steels with Superior Creep Properties for Power Plant Applications

A creep resistant martensitic steel, CPJ-7, was developed with an operating temperature approaching 650°C. Subsequently, another creep resistant martensitic steel, JMP, was designed with the potential to operate at, or slightly above, 650°C. This report describes the development of these alloys from early iterations of CPJ steels to the final CPJ-7 formulation as well as the more heat resistant JMP steel formulation. The design originated from computational modeling for phase stability and precipitate strengthening using fifteen constituent elements. Approximately forty heats of CPJ and ten of JMP, each weighing ~7 kg, were vacuum induction melted. A computationally optimized heat treatment schedule was developed to homogenize the ingots prior to hot forging and rolling prior to final normalization and tempering. Overall, wrought and cast versions of CPJ-7 present superior creep properties when compared to wrought and cast versions of COST alloys for steam turbine and wrought and cast versions of P91/92 for boiler applications. For instance, the Larson Miller Parameter curve for CPJ-7 at 650°C almost coincides with that of COST E at 620°C. The prolonged creep life was attributed to slowing down the process of the destabilization of the MX and M 23 C 6 precipitates at 650°C. On the other hand, the cast version of CPJ-7 also revealed superior mechanical performance, well above commercially available cast 9% Cr martensitic steel or derivatives, especially for creep life at 650°C. The casting process employed slow cooling to simulate the conditions of a thick wall full-size steam turbine casing but utilized a separate homogenization step prior to final normalization and tempering. To advance the development of CPJ-7 for commercial applications, a process was used to scale up the production of the alloy using vacuum induction melting (VIM) and electroslag remelting (ESR), which underlined the importance of melt processing control of intentionally designed minor and trace elements in these advanced alloys. Following the work on CPJ-7, the JMP steels were designed with higher Co for increased solid solution strengthening, Si for oxidation resistance and increased W (with low Mo content) for matrix strength and stability as well as solid solution strengthening. The JMP steels showed increases in creep life compared to CPJ-7 between 118 to 150% at 650°C for testing at various stresses between 138 MPa and 207 MPa. On a Larson-Miller plot, the performance of the JMP steels surpasses that of state-of-the-art MARBN and other MARBN-type steels. The influence of various elements within the composition of the alloys on the microstructure and mechanical properties are discussed. This report presents approximately 420,000 h of in-house creep testing, the equivalent of almost 50 years of cumulative creep tests.

36 MATERIALS SCIENCE↗

Eigenmode analysis of the sheared-flow Z-pinch

Experiments have demonstrated that a Z-pinch can persist for thousands of times longer than the growth time of global magnetohydrodynamic (MHD) instabilities such as the m=0 sausage and m=1 kink modes. These modes have growth times on the order of ta=a/vi, where vi is the ion thermal speed and a is the pinch radius. Axial flows with duz/dr ≲ vi/a have been measured during the stable period, and the commonly accepted theory is that this amount of shear is sufficient to stabilize these modes as predicted by numerical studies using the ideal MHD equations. However, these studies only consider specific equilibrium profiles that typically have a modest magnitude for the logarithmic pressure gradient, qP≡d ln P/d ln r, and may not represent experimental conditions. Linear stability of the sheared-flow Z-pinch is studied here via a direct eigen-decomposition of the matrix operator obtained from the linear ideal MHD equations. Several equilibrium profiles with a large variation of qP are examined. Considering a practical range of k, 1/3 ≲ ka ≲ 10, it is shown that the shear required to stabilize m=0 modes can be expressed as duz/dr≥Cγ0/(ka)α. Here, γ0=γ0(ka) is the profile-specific growth rate in the absence of shear, which scales approximately with |qP|. Both C and α are profile-specific constants, but C is order unity and α≈1. It is further demonstrated that even a large value of shear, duz/dr=3vi/a, is not sufficient to provide linear stabilization of the m=1 kink mode for all profiles considered. This result is in contrast to the currently accepted theory predicting stabilization at much lower shear, duz/dr=0.1vi/a, and suggests that the experimentally observed stability cannot be explained within the linear ideal-MHD model.

Angus, J. R. (ORCID:0000000314740002)↗

Comparing different density-matrix expansions for long-range pion exchange

Empirical energy density functionals (EDFs) are generally successful in describing nuclear properties across the table of nuclides. But their limitations motivate using the density-matrix expansion (DME) to embed long-range pion interactions into a Skyrme functional. Recent results on the impact of the pion were both encouraging and puzzling, necessitating a careful re-examination of the DME implementation. In this study, we take the first steps, focusing on two-body scalar terms in the DME. Exchange energies with long-range one-pion contributions are well approximated by all DME implementations considered, with preference for variants that do not truncate at two derivatives in every EDF term. The use of the DME for chiral pion contributions is therefore supported by this investigation. For scalar-isovector energies it is important to treat neutrons and protons separately. The results are found to apply under broad conditions, although self-consistency is not yet tested.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Theoretical research program to study transition metal trimers and embedded clusters

Small transition metal clusters were studied at a high level of approximation, including all the valence electrons in the calculation and extensive electron correlation, in order to understand the electronic structure of these small metal clusters. By comparison of dimers, trimers, and possibly higher clusters, the information obtained was used to provide insights into the electronic structure of bulk transition metals. Small metal clusters are currently of considerable experimental interest and some information is becomming available both from matrix electron spin resonance studies and from gas phase spectroscopy. Collaboration between theorists and experimentalists is thus expected to be especially profitable at this time since there is some experimental information which can serve to guide the theoretical work.

Walch, S. P.↗

Errors in reduction methods

A mathematical basis is given for comparing the relative merits of various techniques used to reduce the order of large linear and nonlinear dynamics problems during their numerical integration. In such techniques as Guyan-Irons, path derivatives, selected eigenvectors, Ritz vectors, etc., the nth order initial value problem of /y(dot) = f(y) for t greater than 0, y(0) given/ is typically reduced to the mth order (m is much less than n) problem of /z(dot) = g(z) for t greater than 0, z(0) given/ by the transformation y = Pz where P changes from technique to technique. This paper gives an explicit approximate expression for the reduction error e-i in terms of P and the Jacobian of f. It is shown that: (a) reduction techniques are more accurate when the time rate of change of the response y is relatively small; (b) the change in response between two successive stations contributes to the errors at future stations after the change in response is transformed by a filtering matrix H, defined in terms of P; (c) the error committed at a station propagates to future stations by a mixing and scaling matrix G, defined in terms of P, Jacobian and of f, and time increment h. The paper discusses the conditions under which the reduction errors may be minimized and gives guidelines for selecting the reduction basis vector, i.e., the columns of P.

Utku, S.↗

Processable, high temperature polymers from 1,4,5,8-tetrahydro-1,4;5,8-diepoxyanthracene and bis-dienes

1,4,5,8-tetrahydro-1,4;5,8-diepoxyanthracene reacts with various anthracene end-capped polyimide oligomers to form Diels-Alder cycloaddition copolymers. The polymers are soluble in common organic solvents, and have molecular weights of approximately 21,000 to 32,000. These resins exhibit lower weight loss in air than in nitrogen. This is suggested to be due to dehydration (loss of water ranges from 2 to 5 percent) at temperatures of 390 to 400 C to give thermooxidatively stable pentiptycene units along the polymer backbone. Because of their high softening points and good thermooxidative stability, the polymers have potential as processable, matrix resins for high temperature composite applications.

Meador, Mary Ann B.↗

Development of a steady potential solver for use with linearized, unsteady aerodynamic analyses

A full potential steady flow solver (SFLOW) developed explicitly for use with an inviscid unsteady aerodynamic analysis (LINFLO) is described. The steady solver uses the nonconservative form of the nonlinear potential flow equations together with an implicit, least squares, finite difference approximation to solve for the steady flow field. The difference equations were developed on a composite mesh which consists of a C grid embedded in a rectilinear (H grid) cascade mesh. The composite mesh is capable of resolving blade to blade and far field phenomena on the H grid, while accurately resolving local phenomena on the C grid. The resulting system of algebraic equations is arranged in matrix form using a sparse matrix package and solved by Newton's method. Steady and unsteady results are presented for two cascade configurations: a high speed compressor and a turbine with high exit Mach number.

Hoyniak, Daniel↗

The Influence of Correlated Crustal Signals in Modelling the Main Geomagnetic Field

Algorithms used in geomagnetic main-field modelling have for the most part treated the noise in the field measurements as if it were white. A major component of the noise consists of the field due to magnetization in the crust and it has been realized for some time that such signals are highly correlated at satellite altitude. Hence approximation by white noise, while of undoubted utility, is of unknown validity. In this paper we study two plausible statistical models for the crustal magnetization, in which the magnetization is a realization of a stationary, isotropic, random process. At a typical satellite altitude the associated fields exhibit significant correlation over ranges as great as 15 deg. or more, which introduces off-diagonal elements into the covariance matrix, elements that have usually been neglected in modelling procedures. Dealing with a full covariance matrix for a large data set would present a formidable computational challenge, but fortunately most of the entries in the covariance matrix are so small that they can be replaced by zeros. The resultant matrix comprises only about 3 per cent non-zero entries and thus we can take advantage of efficient sparse matrix techniques to solve the numerical system. We construct several main-field models based on vertical-component data from a selected 5 deg. by 5 deg. data set derived from the Magsat mission. Models with and without off-diagonal terms are compared.

Rygaard-Hjalsted, C.↗

Development of microchannel plates in advanced wind-tunnel instrumentation

Microchannel plate (MCP) electron multiplier dynamic range has been increased 3 to 4 orders of magnitude at ambient temperatures, through enhanced input count rate capability and reduced background or 'dark' noise. The previous upper limit of roughly 10(exp 7) - 10(exp 8) cm(exp -2)s(exp -1) at ambient has been extended to levels approach 10(exp 10) cm(exp -2)s(exp -1) under continuous dc operation. The lower limit, previously set by an irreducible background component (approximately 0.6 cm(exp -2)s(exp -1)), has been lowered to the cosmic ray limit of .01 cm(exp -2)s(exp -1). The high end improvement was achieved by conductively cooling a very low resistance MCP by bonding it to a heat sink, while maintaining pulse-counting operation with multianode readouts. The low-end improvement was achieved by removing all radioisotopes from the MCP matrix glass. The detectors will benefit optical and mass spectrometry, flow visualization, plasma diagnostics, magnetometry, and other high signal flux applications. Very low MCP background noise will benefit X-ray and UV astronomy, medical imaging, trace isotope mass spectrometry, and other applications where the signal flux is often extremely low.

Feller, W. Bruce↗

Changes in gravity inhibit lymphocyte locomotion through type I collagen

Immunity relies on the circulation of lymphocytes through many different tissues including blood vessels, lymphatic channels, and lymphoid organs. The ability of lymphocytes to traverse the interstitium in both nonlymphoid and lymphoid tissues can be determined in vitro by assaying their capacity to locomote through Type I collagen. In an attempt to characterize potential causes of microgravity-induced immunosuppression, we investigated the effects of simulated microgravity on human lymphocyte function in vitro using a specialized rotating-wall vessel culture system developed at the Johnson Space Center. This very low shear culture system randomizes gravitational vectors and provides an in vitro approximation of microgravity. In the randomized gravity of the rotating-wall vessel culture system, peripheral blood lymphocytes did not locomote through Type I collagen, whereas static cultures supported normal movement. Although cells remained viable during the entire culture period, peripheral blood lymphocytes transferred to unit gravity (static culture) after 6 h in the rotating-wall vessel culture system were slow to recover and locomote into collagen matrix. After 72 h in the rotating-wall vessel culture system and an additional 72 h in static culture, peripheral blood lymphocytes did not recover their ability to locomote. Loss of locomotory activity in rotating-wall vessel cultures appears to be related to changes in the activation state of the lymphocytes and the expression of adhesion molecules. Culture in the rotating-wall vessel system blunted the ability of peripheral blood lymphocytes to respond to polyclonal activation with phytohemagglutinin. Locomotory response remained intact when peripheral blood lymphocytes were activated by anti-CD3 antibody and interleukin-2 prior to introduction into the rotating-wall vessel culture system. Thus, in addition to the systemic stress factors that may affect immunity, isolated lymphocytes respond to gravitational changes by ceasing locomotion through model interstitium. These in vitro investigations suggest that microgravity induces non-stress-related changes in cell function that may be critical to immunity. Preliminary analysis of locomotion in true microgravity revealed a substantial inhibition of cellular movement in Type I collagen. Thus, the rotating-wall vessel culture system provides a model for analyzing the microgravity-induced inhibition of lymphocyte locomotion and the investigation of the mechanisms related to lymphocyte movement.

Weightlessness↗

Dynamics of Polymerization and Gelation in Epoxy Nanocomposites via X-ray Photon Correlation Spectroscopy

The details of the curing process of epoxy resins are notoriously difficult to ascertain, hampering polymer matrix composite manufacturing. Here, the curing of a series of six epoxy resins containing dilute spherical nanoparticles was investigated via X-ray photon correlation spectroscopy (XPCS). The resin formulation was varied to achieve a range of crosslink densities and topologies. Prior to gelation, the viscosities obtained from XPCS agreed with bulk rheology. A clear dynamic transition was detected in all samples, where the q-scaling changes from approximately q2 .2 to q 1 , indicating a transition from sub-diffusive to ballistic-like nanoparticle motion. For the samples with a high crosslink density, this change was abrupt and coincided with gelation, and there was no evidence of crosslink heterogeneity. Chain entanglement played a role in the transition for samples with lower crosslink density. Finally, this work reveals the potential of XPCS to deliver new quantitative insights into the curing process of thermosetting resins.

36 MATERIALS SCIENCE↗

Data-driven linear time advance operators for the acceleration of plasma physics simulation

In this study, we demonstrate the application of data-driven linear operator construction for time advance with a goal of accelerating plasma physics simulation. We apply dynamic mode decomposition (DMD) to data produced by the nonlinear SOLPS-ITER (Scrape-off Layer Plasma Simulator - International Thermonuclear Experimental Reactor) plasma boundary code suite in order to estimate a series of linear operators and monitor their predictive accuracy via online error analysis. We find that this approach defines when these dynamics can be represented by a sequence of approximate linear operators and is essential for providing consistent projections when compared to an unconstrained application. For linear diffusion and advection–diffusion fluid test problems, we construct and apply operators within explicit and implicit time advance schemes, demonstrating that stability can be robustly guaranteed in each case. We further investigate the use of the linear time advance operators within several integration methods including forward Euler, backward Euler, and the matrix exponential. The application of this method to simulation data from SOLPS-ITER, with varying levels of Markov chain Monte Carlo numerical noise, shows that constrained DMD operators yield a capability to identify, extract, and integrate a (slow) subset of the present timescales. Example applications show that for projected speedup factors of [Formula: see text], and [Formula: see text], a mean relative error of 3%, 5%, and 8% and maximum relative error less than 20% are achievable, which appears acceptable for typical SOLPS-ITER steady-state simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Classical shadows for quantum process tomography on near-term quantum computers

Quantum process tomography is a powerful tool for understanding quantum channels and characterizing the properties of quantum devices. Inspired by recent advances using classical shadows in quantum state tomography [H.-Y. Huang, R. Kueng, and J. Preskill, .], we have developed ShadowQPT, a classical shadow method for quantum process tomography. We introduce two related formulations with and without ancilla qubits. ShadowQPT stochastically reconstructs the Choi matrix of the device allowing for an classical evaluation of the device on arbitrary inputs with respect to arbitrary outputs. Using shadows, we then show how to compute overlaps, generate all k -weight reduced processes, and perform reconstruction via Hamiltonian learning. These latter two tasks are efficient for large systems as the number of quantum measurements needed scales only logarithmically with the number of qubits. A number of additional approximations and improvements are developed, including the use of a pair-factorized Clifford shadow and a series of postprocessing techniques that significantly enhance the accuracy for recovering the quantum channel. We have implemented ShadowQPT using both Pauli and Clifford measurements on the IonQ trapped ion quantum computer for quantum processes up to n = 4 qubits, and we achieved good performance. Published by the American Physical Society 2024

Levy, Ryan (ORCID:0000000349527156)↗

A NASTRAN implementation of the doubly asymptotic approximation for underwater shock response

A detailed description is given of how the decoupling approximation known as the doubly asymptotic approximation is implemented with NASTRAN to solve shock problems for submerged structures. The general approach involves locating the nonsymmetric terms (which couple structural and fluid variables) on the right hand side of the equations. This approach results in coefficient matrices of acceptable bandwidth but degrades numerical stability, requiring a smaller time step size than would otherwise be used. It is also shown how the structure's added (virtual) mass matrix, is calculated with NASTRAN.

Everstine, G. C.↗

Formation chronology for C2 meteorites

The times at which phyllosilicate matrix and euhedral olivines became associated have been determined for five C2 meteorites. The ages, calculated from fission track densities on crystal surfaces, are based on an initial Pu-244/U-238 ratio in the matrix material of 0.0154 at 4.6 b.y., and range from 4.22 b.y. for Nogoya to 4.42 b.y. for Murray. Unless the initial Pu-244/U-238 ratio was less than 0.004, the meteorites cannot have existed in their present form for 4.6 billion years. The measured ages place limits on the time when pre-compaction effects such as micrometeorite craters and solar flare tracks were produced, and they may approximately date the formation of the olivine crystals themselves.

Macdougall, J. D.↗