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 631 records · Page 35

Laboratory studies of the newly discovered infrared band at 4705.2 cm-1 (2.1253 micrometers) in the spectrum of Io: the tentative identification of CO2

We discuss over 120 laboratory experiments pertaining to the identification of the new absorption band discovered by Trafton et al. (1991) at 4705.2 cm-1 (2.1253 micrometers) in the spectrum of Io. It is shown that this band is not due to overtones or combinations of the fundamental bands associated with the molecules (or their chemical complexes) already identified on Io, namely, SO2, H2S, and H2O. Thus, this band is due to a new, previously unidentified, component of Io. Experiments also demonstrate that the band is not due to molecular H2 frozen in SO2 frosts. Since the frequency of this band is very close to the first overtone of the nu 3 asymmetric stretching mode of CO2, we have investigated the spectral behavior of CO2 under a variety of conditions appropriate for Io. The profile of the Io band is not consistent with the rotational envelope expected for single, freely rotating, gaseous CO2 under Io-like conditions. It was found that pure, solid CO2 and CO2 intimately mixed in a matrix of solid SO2 and H2S produce bands with similar widths (5-10 cm-1), but that these bands consistently fall at frequencies about 10-20 cm-1 (approximately 0.007 micrometer) lower than the Io band. CO2 in SO2 : H2S ices also produces several additional bands that are not in the Io spectra. The spectral fit improves, however, as the CO2 concentration in SO2 increases, suggesting that CO2-CO2 interactions might be involved. A series of Ar : CO2 and Kr : CO2 matrix isolation experiments, as well as laboratory work done elsewhere, show that CO2 clustering shifts the band position to higher frequencies and provides a better fit to the Io band. Various laboratory experiments have shown that gaseous CO2 molecules have a propensity to cluster between 80 and 100 K, temperatures similar to those found on the colder regions of Io. We thus tentatively identify the newly discovered Io band at 4705.2 cm-1 (2.1253 micrometers) with CO2 multimers or "clusters" on Io. Whether these clusters are buried within an SO2 frost, reside on the surface, or are in a residual, steady-state "atmospheric aerosol" population over local coldtraps is not entirely clear, although we presently favor the latter possibility. The size of these clusters is not well defined, but evidence suggests groups of more than four molecules are required. The absorption strength of the 2 nu 3 CO2 cluster overtone determined in the laboratory, in conjunction with the observed strength of the Io band, suggests that the disk-integrated abundance of CO2 is less than 1% that of the SO2. Studies of the sublimation behavior of CO2 indicate that it probably resides predominantly in the cooler areas (< 100 K) of Io. The relative constancy of the Io feature over a variety of orbital phases suggests that the polar regions may contain much of the material. Some consequences of the physical properties of CO2 under conditions pertinent to Io are discussed. The presence of CO2 clusters on Io could be verified by the detection of any one of several other infrared bands associated with the CO2 molecule, of which the strongest are the nu 3 12CO2 asymmetric stretch fundamental near 2350 cm-1 (4.25 micrometers) and the nu 2 bending mode fundamental near 660 cm-1 (15.1 micrometers). Weaker bands that may also be detectable include the nu 3 13CO2 asymmetric stretch fundamental near 2280 cm-1 (4.39 micrometers), the 2 nu 2 + nu 3 combination/overtone band near 3600 cm-1 (2.78 micrometers), and the nu 1 + nu 3 combination band near 3705 cm-1 (2.70 micrometers).

NASA Center ARC↗

Destructive PIE and Safety Testing of Six AGR-5/6/7 Capsule 2 Compacts

This study evaluates fission product retention and particle failure mechanisms in AGR-5/6/7 Capsule 2 uranium carbide and uranium oxide (UCO) tristructural isotropic (TRISO) fuel under high-temperature gas reactor accident-relevant conditions using high-temperature safety tests and destructive post-irradiation examination. Three Capsule 2 compacts were held isothermally at 1600°C for approximately 300 hours and one compact at 1800°C for approximately 300 hours; two additional compacts were examined in the as-irradiated state. Post-test deconsolidation–leach–burn–leach (DLBL) quantified nuclide inventories in matrix and particles. Individual particles were surveyed for radioisotope inventories, and microanalytical approaches resolved microstructural evolution and fission product distributions within the coating layers. At 1600°C, no krypton was detected above the minimum detectable limit, and cesium releases were far below a single particle equivalent, indicating the absence of full TRISO failure or SiC failures. Silver releases were limited and primarily reflected depleted post-irradiation inventories, consistent with prior compact-level exams indicating substantial in-pile 110m Ag loss. At 1800°C, cumulative 134 Cs release of approximately 2.5 particle equivalents and delayed 85 Kr totaling approximately 0.53 particle equivalents were consistent with one full TRISO failure and two SiC failures. Europium and strontium releases were roughly one order of magnitude higher than at 1600°C and comparable to AGR-1/AGR-2 high-temperature tests, with sustained late-hold rates indicating diffusion through intact coatings coupled with matrix depletion. Overall, AGR-5/6/7 Capsule 2 UCO fuel demonstrated fission product retention during safety testing consistent with prior AGR campaigns, while distinctive in-pile 110m Ag depletion and measurable 1600°C europium loss motivate targeted follow-on studies.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Enhancing photoionization rate calculations in low-temperature plasmas using spectral methods

Photoionization plays a central role in the development of streamer discharges and other non-equilibrium plasma phenomena. It creates seed electrons, which are essential for positive streamer propagation, allowing the ionization front to move forward. Because of this, accurate modeling of photoionization is very important for predicting streamer behavior and plasma evolution. The photoionization process in air (N 2 – O 2 mixture) is often described by the Zheleznyak model (1982). This model is usually solved through Helmholtz-type equations that approximate the Zheleznyak photoionization model (Zheleznyak et al. 1982) as Partial Differential Equations (PDEs). Conventional numerical methods, such as the Finite Difference Method (FDM) or Finite Volume Method (FVM), are widely used to solve these equations. Although they are prevalent, the computational cost of these methods due to their need for matrix operations and iterative solver is demanding. To address this challenge, this work develops a spectral solver based on the Fast Fourier Transform (FFT) combined with Discrete Cosine Transform (DCT) and Discrete Sine Transform (DST) to calculate the photoionization rate efficiently in an axisymmetric cylindrical domain. This method naturally satisfies the boundary conditions used in the model and converts the PDE into algebraic ones in spectral space. Thus, avoids the need for iterative matrix solvers. When compared with FDM results, it is demonstrated that the new solver not only maintains accuracy, but also reduces the computational cost, showing a performance increase of approximately 100 compared to FDM over a wide range of problem sizes. The method is parallelized using Message Passing Interface (MPI) and has been integrated into a fluid plasma model for streamer simulation. Here, this FFT-based approach provides a fast and reliable alternative for calculating photoionization in fluid models, helping large-scale plasma simulations run faster and efficiently, and allows higher-resolution simulation without extra computational cost.

Axisymmetric system↗

Enabling attractive-repulsive potentials in binary-collision-approximation monte-carlo codes for ion-surface interactions

Abstract Binary Collision Approximation (BCA) codes for ion-material interactions, such as SRIM, Tridyn, F-TRIDYN, and SDtrimSP, have historically been limited to screened Coulomb potentials even at low energies due to the difficulty in numerically solving the Distance of Closest Approach (DOCA) problem for attractive-repulsive potentials. Techniques such as direct n-body simulation or modifications to Newton’s method are either prohibitively costly or not guaranteed to work for all potentials. Advanced rootfinding techniques, such as companion matrix solvers, offer a solution. For many attractive-repulsive potentials, however, a companion matrix cannot be used directly, because there is no way to put the associated functions into a monomial basis form. A complementary technique is proxy rootfinding—by finding the best-fit polynomial approximant of a function, the zeros of the approximant can be guaranteed to be close to the zeros of the function. Using the Chebyshev basis and grid offers additional guarantees with regards to the quality of the approximation, the speed of convergence, and the avoidance of Runge’s phenomenon. By finding Chebyshev interpolants and using the Chebyshev-Frobenius companion matrix, the zeros of any real function on a bounded domain can be found. Here we show that using an Adaptive Chebyshev Proxy Rootfinder with Automatic Subdivision (ACPRAS) with appropriate scaling functions, numerical issues presented by attractive-repulsive potentials, including those of scale, can be handled. Using these techniques, we show that it is possible to include any physically reasonable interatomic potential in a BCA code, and to guarantee correctness of the resulting scattering angle calculations.

Materials Science↗

Matrix Theory for Data Association in PVS

Consider a collection of data generated by sensors from a set of aircraft. Data association is the process of connecting each sensor measurement with its corresponding aircraft. Furthermore once the data association has taken place, the state of the aircraft can be approximated using a Kalman filter. This talk aims to explore formal specification and verification of data association in the Prototype Verification System (PVS). Formal specification and verification of data association includes development of Kalman filters, Mahalanobis distance, and other topics of matrix analysis in PVS.

Linear Algebra↗

Improving the chi-squared approximation for bivariate normal tolerance regions

Let X be a two-dimensional random variable distributed according to N2(mu,Sigma) and let bar-X and S be the respective sample mean and covariance matrix calculated from N observations of X. Given a containment probability beta and a level of confidence gamma, we seek a number c, depending only on N, beta, and gamma such that the ellipsoid R = (x: (x - bar-X)'S(exp -1) (x - bar-X) less than or = c) is a tolerance region of content beta and level gamma; i.e., R has probability gamma of containing at least 100 beta percent of the distribution of X. Various approximations for c exist in the literature, but one of the simplest to compute -- a multiple of the ratio of certain chi-squared percentage points -- is badly biased for small N. For the bivariate normal case, most of the bias can be removed by simple adjustment using a factor A which depends on beta and gamma. This paper provides values of A for various beta and gamma so that the simple approximation for c can be made viable for any reasonable sample size. The methodology provides an illustrative example of how a combination of Monte-Carlo simulation and simple regression modelling can be used to improve an existing approximation.

Feiveson, Alan H.↗

Multimodal photon and neutron imaging devices with fast event-by-event multisampling readout and optimized particle type discrimination

A particle imaging method for distinguishing between types of incident particles, such as neutrons, photons, and alphas, and improving the position resolution of particle imaging devices with matrix readout. The method includes high frequency multisampling readout electronics that provides the sequences of multiple measurements for each detected event, resulting in recorded detailed waveform information describing the signals. Such detailed information is used to approximate each signal waveform with a parameterized function in which the extracted parameter sets determine the type of the incident particle in an optimized fashion. The detailed event-by-event multisampling information for each signal readout channel in the matrix readout of the radiation imaging devices improves and optimizes the position resolution for variable shapes of the signals. Such devices can be used in mixed radiation fields, creating a new class of multimodal photon and neutron imagers.

Degtiarenko, Pavel V.↗

Design of Silicide-Strengthened Nb–Si–Cr–(Mo) alloys for additive manufacturing

Three high–intermetallic volume Nb–Si–Cr–(Mo) alloys were designed using CALPHAD modeling with the goal of identifying high–specific strength, oxidation-resistant alloys that can be additively manufactured using powder bed fusion. The silicides Nb 5 Si 3 and Nb 9 Si 2 Cr3 were targeted as the primary strengthening phases, and the addition of Cr promoted the NbCr 2 phase. These alloys were cast and surface-processed with electron beam welding at different speeds to simulate additive manufacturing, and the phases and microstructures of both cast and welded regions were characterized. The weld processing was found to produce fine-grained microstructures in each alloy with fine-scale intermetallics uniformly distributed among a body-centered cubic Nb matrix. Microstructural refinement and hardness were found to increase with weld velocity; one alloy reached its highest hardness of approximately 16 GPa before the brittleness at higher velocities became detrimental. One alloy was found to be qualitatively the least brittle while also attaining a hardness of 13 GPa and was therefore identified as a good candidate for additive manufacturing.

36 MATERIALS SCIENCE↗

The polarized-signal density matrix: A practical way to recover molecular frame information from isotropic samples

We present a novel approach to model ultrafast time-dependent nonlinear optical polarization sensitive signals emitted from randomly oriented molecules. By projecting the laboratory-frame analyzer polarization axis into the molecular-frame and linking that axis with the density matrix through a tensor product, we demonstrate an approach to find a specific molecular orientation that yields a good approximation to simulated four-wave mixing signals produced by the same model but with averaging over molecular orientation.

Thurston, Richard L↗

An Application of Molecular Recognition for the Efficient Removal of Cesium from Hanford Nuclear Waste by Modular Solvent Extraction

In this work, experimental results leading to flowsheet design are presented showing how a calixarene-crown ether based solvent-extraction process can meet the challenge of cesium removal from nuclear tank wastes stored at the US Department of Energy Hanford site. Cleanup of legacy Cold War nuclear waste stored in underground tanks represents one of the greatest environmental challenges facing the US Department of Energy in terms of risk, cost, and effectiveness of applicable science and technology. Planning for the cleanup at the Hanford Site calls for the removal of the radioactive fission product 137Cs from its alkaline salt waste, including the use of modular processes that can be deployed near the tank farms. To meet the resulting need for extremely high selectivity, the Next-Generation Caustic-Side Solvent Extraction (NG-CSSX) process employing a calix[4]arene-crown ether in modified kerosene has been adapted to remove sub-millimolar cesium in competition with molar sodium and potassium in a high-nitrate alkaline matrix. Potassium loading in the solvent was determined in extraction, scrubbing, and stripping, leading to an empirical model closely approximating cesium distribution ratios for a variety of Hanford waste types. Process chemistry has been developed based on this molecular-recognition approach, focusing on the competitive effect of potassium loading and the mitigating process modifications needed, including extending the scrub section. The result is a modular flowsheet design that can achieve cesium decontamination factors well in excess of 15,000 even for the worst-case Hanford waste.

Williams, Neil [ORNL] (ORCID:000000023159226X)↗

A Polarization Tensor for Neutron Optics

A rank 2 dimension 3 tensor formulation for polarized neutron transport in a magnetic field is derived. Using this formulation, the rotating magnetic field problem is revisited to produce a single 3 × 3 vector-transformation matrix that describes the polarization transport. That result is then used to fashion a new strategy for a numerical approximation of polarization transport in an arbitrary magnetic field using time-dependent perturbation theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Thin film design of amorphous hafnium oxide nanocomposites enabling strong interfacial resistive switching uniformity

A design concept of phase-separated amorphous nanocomposite thin films is presented that realizes interfacial resistive switching (RS) in hafnium-oxide-based devices. The films are formed by incorporating an average of 7% Ba into hafnium oxide during pulsed laser deposition at temperatures ≤400°C. The added Ba prevents the films from crystallizing and leads to ~20-nm-thin films consisting of an amorphous HfO x host matrix interspersed with ~2-nm-wide, ~5-to-10-nm-pitch Ba-rich amorphous nanocolumns penetrating approximately two-thirds through the films. This restricts the RS to an interfacial Schottky-like energy barrier whose magnitude is tuned by ionic migration under an applied electric field. Resulting devices achieve stable cycle-to-cycle, device-to-device, and sample-to-sample reproducibility with a measured switching endurance of ≥10 4 cycles for a memory window ≥10 at switching voltages of ±2 V. Each device can be set to multiple intermediate resistance states, which enables synaptic spike-timing–dependent plasticity. The presented concept unlocks additional design variables for RS devices.

36 MATERIALS SCIENCE↗

The small motor plume/material impingement test program at MSFC utilizing nine different propellant formulations (April and May 1974), section 2

The relative effects of several candidate SSRM propellant formulations and their plume impingement effects on HRSI and RCC materials were evaluated. Nine solid propellant formulations were tested. The selected propellant matrix allowed an evaluation of propellants with and without metal additives, with and without burning rate catalyst, and low (approximately 1927 C) and high (approximately 2649 C) combustion temperatures. Motors were fired at a simulated SRB staging altitude of 3.96 km (130,000 ft) (nominal). The altitude pressure was predicted to drop approximately 0.6 km (20,000 ft) during a motor firing. All motors were loaded with 1.8 to 2.3 kg (4 to 5 lb) of propellant and burned for approximately 2 s.

Source record↗

Convergence of the SMI and the diagonally loaded SMI algorithms with weak interference

Approximations for the power levels at the output of an adaptive array that uses the diagonally loaded sample-matrix inversion (SMI) algorithm are derived. Diagonal loading is a technique where the diagonal of the covariance matrix is augmented with a positive or negative constant prior to inversion. This paper examines how the signal-to-interference-plus-noise ratio (SINR) and signal-to-interference ratio (SIR) at the array output vary with the number of samples taken when the input signals are CW. It is shown that positive loading produces more rapid convergence with a reduction in output SIR. Negative loading provides an improved SIR level, but positive loading produces more rapid convergence. Negative loading provides an improved SIR level, but the output power levels are erratic and slow to converge. Simulation results which verify the theoretical procedure are given.

Ganz, Matthew W.↗

Wave Number Selection for Incompressible Parallel Jet Flows Periodic in Space

The temporal instability of a spatially periodic parallel flow of an incompressible inviscid fluid for various jet velocity profiles is studied numerically using Floquet Analysis. The transition matrix at the end of a period is evaluated by direct numerical integration. For verification, a method based on approximating a continuous function by a series of step functions was used. Unstable solutions were found only over a limited range of wave numbers and have a band type structure. The results obtained are analogous to the behavior observed in systems exhibiting complexity at the edge of order and chaos.

Miles, Jeffrey Hilton↗

Characterization of SiC Fiber (SCS-6) Reinforced-Reaction-Formed Silicon Carbide Matrix Composites

Silicon carbide fiber (SCS-6) reinforced-reaction-formed silicon carbide matrix composites were fabricated using a reaction-forming process. Silicon-2 at.% niobium alloy was used as an infiltrant instead of pure silicon to reduce the amount of free silicon in the matrix after reaction forming. The matrix primarily consists of silicon carbide with a bimodal grain size distribution. Minority phases dispersed within the matrix are niobium disilicide (NbSi2), carbon, and silicon. Fiber pushout tests on these composites determined a debond stress of approximately 67 MPa and a frictional stress of approximately 60 MPa. A typical four-point flexural strength of the composite is 297 MPa (43.1 KSi). This composite shows tough behavior through fiber pullout.

Singh, M.↗

Optimization of an oligonucleotide microchip for microbial identification studies: a non-equilibrium dissociation approach

The utility of a high-density oligonucleotide microarray (microchip) for identifying strains of five closely related bacilli (Bacillus anthracis, Bacillus cereus, Bacillus mycoides, Bacillus medusa and Bacillus subtilis) was demonstrated using an approach that compares the non-equilibrium dissociation rates ('melting curves') of all probe-target duplexes simultaneously. For this study, a hierarchical set of 30 oligonucleotide probes targeting the 16S ribosomal RNA of these bacilli at multiple levels of specificity (approximate taxonomic ranks of domain, kingdom, order, genus and species) was designed and immobilized in a high-density matrix of gel pads on a glass slide. Reproducible melting curves for probes with different levels of specificity were obtained using an optimized salt concentration. Clear discrimination between perfect match (PM) and mismatch (MM) duplexes was achieved. By normalizing the signals to an internal standard (a universal probe), a more than twofold discrimination (> 2.4x) was achieved between PM and 1-MM duplexes at the dissociation temperature at which 50% of the probe-target duplexes remained intact. This provided excellent differentiation among representatives of different Bacillus species, both individually and in mixtures of two or three. The overall pattern of hybridization derived from this hierarchical probe set also provided a clear 'chip fingerprint' for each of these closely related Bacillus species.

Non-NASA Center↗

On Kalman filter solution of space-time interpolation

The approximate Kalman filtering algorithm presented in [1] for image sequence processing can introduce unacceptable negative eigenvalues in the information matrix and can have degraded performance in some applications. The improved algorithm presented in this note guarantees a positive definite information matrix, leading to more stable filter performance.

data↗