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 973 records · Page 54

3D Representation of UAV-obstacle Collision Risk Under Off-nominal Conditions

Safe operations of autonomous unmanned aerial vehicles (UAVs) in low-altitude airspace with beyond visual line-of-sight (BVLOS) flights demand robust risk monitoring of airspace as well as of people and property on ground. One of the safety critical factors for UAV flights is the risk of collision with static and dynamic obstacles in proximity to its flight path. This paper presents a detailed formulation of risk of obstacle collision incorporating the effects of off-nominal conditions introduced by component failures, degraded controllability and environmental disturbances such as wind gusts. The risk is represented in terms of a matrix with rows corresponding to the likelihood of occurrence of collision and columns representing severity of collision to the vehicle and surrounding structures. Risk likelihood is generated using a Bayesian Belief Network (BBN) that compiles knowledge from related Failure Modes and Effects Analysis (FMEAs) and Subject Matter Experts (SMEs) to determine the probability of collision based on on-board sensor measurements indicative of vehicle health and controllability. Risk severity is computed utilizing a point-mass 3D kinematic model of the vehicle in presence of wind. The proposed risk factor is demonstrated on real flight data from experimental flights of an octocopter at NASA Langley Research Center in presence of simulated obstacles and wind conditions. Effect of varying wind conditions, level of controllability and obstacle measurement noise on the risk factor is demonstrated. The proposed approach enables risk-informed decision making for timely mitigation of current and future unsafe events in autonomous systems.

risk analysis↗

A Parallelized Oxidation-Driven Surface Recession Framework in DSMC Code, SPARTA

Spacecrafts rely on ablative thermal protection systems (TPS) made of composites consisting of a carbon-based reinforcement and a polymeric matrix. These materials are designed to withstand high-temperature oxidation and surface recession during re-entry into the Earth's atmosphere. However, ablation occurs due to a complex interplay of thermal, mechanical, and chemical factors, making it challenging to determine the individual impact of each on the TPS's overall degradation. In this study, we have developed an ablation model that can leverage a finite rate carbon oxidation model to predict material recession and surface states more accurately. Stochastic PArallel Rarified-gas Time-accurate Analyzer (SPARTA), a direct-simulation Monte Carlo (DSMC) code, is modified to allow oxidation-driven ablation of implicitly defined carbon surfaces. In SPARTA, implicit surfaces are generated from the grid corner point values via a marching cubes algorithm, therefore creating a new set of surface elements every time ablation is performed. The finite-rate oxidation model developed by Gopalan et. al can perform both gas-surface and pure-surface reactions and is now adapted to tally surface data on a per grid cell basis. The ablation functionality was also adjusted so once the reactions have occurred, the number of reactions leading to CO formation can be converted to corner point reduction values; therefore, carbon removal is directly proportional to surface recession. We also briefly discuss some unique challenges associated with parallelizing this dynamic surface state and geometry. Finally, we analyze the performance of this parallelized implicit chemistry model with simple 2D and 3D benchmark cases by producing surface state statistics, area changes over time, and visualization across a range of surface temperatures and processors with and without load-balancing.

DSMC↗

The GeneLab Buffet: A Bioinformatic MATRIX of MANGO and TOAST

The GeneLab data repository provides an unparalleled resource for exploring how spaceflight affects organisms with omics-level insights. However, two major interlinked challenges to capitalizing on the information within these data are their vast breadth and the often-specialized expertise that has been required in the past for their analysis. How do you compare responses within and between studies, especially if you are a non-bioinformatics specialist? This presentation will discuss how Space Biology data can be accessed using software to help provide these data resources to address research questions and generate new hypotheses. The presentation will cover a wide range of the available space life science tools but will focus on TOAST, MANGO, the MATRIX, RadBioApp and other interactive relational databases (https://genelab.nasa.gov/external-vis-apps). These exploration environments have been developed to search the GeneLab data repository for new insights that inform how model organisms respond to microgravity, radiation and other factors associated with spaceflight. The presentation will be interactive, and participants will have the opportunity to ask questions and learn more about the data viz and modeling tools that are available to them.

AstroBotany↗

Nonrenormalization and Operator Mixing via On-Shell Methods

Using on-shell methods, we present a new perturbative nonrenormalization theorem for operator mixing in massless four-dimensional quantum field theories. By examining how unitarity cuts of form factors encode anomalous dimensions, we show that longer operators are often restricted from renormalizing shorter operators at the first order where Feynman diagrams exist. The theorem applies quite generally and depends only on the field content of the operators involved. We apply our theorem to operators of dimension five through seven in the standard model effective field theory, including examples of nontrivial zeros in the anomalous-dimension matrix at one through four loops. The zeros at two and higher loops go beyond those previously explained using helicity selection rules. We also include explicit sample calculations at two loops.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Linear and Nonlinear Solvers for Simulating Multiphase Flow within Large-Scale Engineered Subsurface Systems

Simulation of multiphase flow in the subsurface is well-known to be computationally challenging. While there have been many studies that have explored approaches to overcoming these challenges, they often utilize relatively simple case studies. In this paper, we focus on the unique numerical challenges posed by modeling large-scale engineered subsurface systems, characterized by discrete features embedded in a heterogeneous natural subsurface setting. The man-made features such as shafts, tunnels, and barriers often cause multiple challenges in modeling the domain for multiphase porous media flow. This flow scenario can have a wide range of applications such as nuclear waste repositories, enhanced recovery of a petroleum reservoir, geothermal engineering, and carbon sequestration. An example of these severe numerical challenges is the case of performance assessment (PA) for Waste Isolation Pilot Plant (WIPP), the only operating deep geological repository in the US, which simulates extreme material properties of bedded salt rock formation and extreme contrast due to open excavation next to the formation. The models have extremes not only of permeability and porosity but also of the constitutive models needed for multiphase flow; additionally, they have process models like salt creep closure reducing porosity over time, fracturing in clay and anhydrite interbeds of the bedded salt, gas generation from the waste materials, and unintentional human borehole intrusions in some scenarios. Numerical simulations require the solution of coupled systems of nonlinear PDEs; in our work, we use the open-source simulator PFLOTRAN which is based on Finite Volume discretization. The solution of the nonlinear equations requires use of the Newton-Raphson iteration at each time step, which entails the solution of the linearized Jacobian system at each iteration. The effects of all the processes (i.e., large number of unknowns, highly nonlinear constitutive relations, large contrasts in material properties in short distances) lead to an ill-conditioned Jacobian matrix that severely challenges traditional linear solver, i.e., stabilized biconjugate gradient with block Jacobi incomplete LU preconditioner (BCGS-ILU) leading to non-convergence for traditional Newton-Raphson nonlinear solver causing unacceptably long computation time for each model. This paper presents linear solvers such as constrained pressure residual (CPR) two-stage preconditioner with alternate-block-factorization (ABF) and quasi- implicit pressure and explicit saturation (QIMPES) decouplers and flexible generalized residual solver (FGMRES). The new general-purpose nonlinear solver, Newton trust-region dogleg Cauchy (NTRDC), is also introduced to resolve extreme nonlinearities in the models. We demonstrate the effectiveness of each method relative to the default BCGS-Newton solver. The two best cases had nearly 50 times speed-up and achieved completion of a simulation in 14 hours that never completed due to non-convergence with the default solver. We also investigate the strong scalability of each method and discuss some of the deficiencies found for Block Jacobi preconditioner using parallel domain decomposition, and node packing effects of modern processor architecture.

Preconditioner, Nonlinear, Porous media, Multiphas↗

Study of extraterrestrial disposal of radioactive wastes. Part 2: Preliminary feasibility screening study of extraterrestrial disposal of radioactive wastes in concentrations, matrix materials, and containers designed for storage on earth

The results are reported of a preliminary feasibility screening study for providing long-term solutions to the problems of handling and managing radioactive wastes by extraterrestrial transportation of the wastes. Matrix materials and containers are discussed along with payloads, costs, and destinations for candidate space vehicles. The conclusions reached are: (1) Matrix material such as spray melt can be used without exceeding temperature limits of the matrix. (2) The cost in mills per kw hr electric, of space disposal of fission products is 4, 5, and 28 mills per kw hr for earth escape, solar orbit, and solar escape, respectively. (3) A major factor effecting cost is the earth storage time. Based on a normal operating condition design for solar escape, a storage time of more than sixty years is required to make the space disposal charge less than 10% of the bus-bar electric cost. (4) Based on a 10 year earth storage without further processing, the number of shuttle launches required would exceed one per day.

Hyland, R. E.↗

Partitioning sparse matrices with eigenvectors of graphs

The problem of computing a small vertex separator in a graph arises in the context of computing a good ordering for the parallel factorization of sparse, symmetric matrices. An algebraic approach for computing vertex separators is considered in this paper. It is shown that lower bounds on separator sizes can be obtained in terms of the eigenvalues of the Laplacian matrix associated with a graph. The Laplacian eigenvectors of grid graphs can be computed from Kronecker products involving the eigenvectors of path graphs, and these eigenvectors can be used to compute good separators in grid graphs. A heuristic algorithm is designed to compute a vertex separator in a general graph by first computing an edge separator in the graph from an eigenvector of the Laplacian matrix, and then using a maximum matching in a subgraph to compute the vertex separator. Results on the quality of the separators computed by the spectral algorithm are presented, and these are compared with separators obtained from other algorithms for computing separators. Finally, the time required to compute the Laplacian eigenvector is reported, and the accuracy with which the eigenvector must be computed to obtain good separators is considered. The spectral algorithm has the advantage that it can be implemented on a medium-size multiprocessor in a straightforward manner.

Pothen, Alex↗

Stability of Pseudobrookite-Type Titanium Oxides

Orthorhombic, (Bbmm), (Al, Fe, Cr, Ti)(sub 2) TiO5-(Mg, Fe)Ti2O5 solid solutions (pseudobrookites, s.l.) are found either as an oxidation product of ilmenite and/or spinel or a primary crystallizing phase in igneous and metamorphic rocks on Earth (e.g., basalt flows, crustal and mantle xenoliths, hornfels), and basaltic rocks on the Moon. Moreover, orthorhombic oxides are often part of the crystalline matrix in glass/ceramics with useful applications, and play a major role in the industrial production of TiO2. To fully exploit the potential of these compounds as petrogenetic indicators and/or useful materials we need to quantitatively understand the factors controlling their properties and stability, and thus, to extrapolate beyond the calibrating experiments. For that purpose, we need to combine thermochemistry, phase equilibrium, and in situ P-V-T-cation disorder experimental data that presently either are incomplete or lacking. Perhaps, the most complete data set is that for MgTi2O5 (karrooite) which allows the calibration of models for the Gibbs free energy of the MgTi2O5 as a function of pressure, temperature, and the Mg2+-Ti4+ distribution between the two nonequivalent octahedral sites. Consequently, the effect of cation disorder on MgTi2O5 stability, and the phase relations among MgTi2O5, other titanium oxides, and silicate minerals can be examined. Calculated phase relations in the Mg-Ti-Si-O system and phase equilibrium experiments in Fe-bearing compositions suggest that pseudobrookite-type oxides may be a more common in rocks than previously realized. However, homogeneous and heterogeneous equilibria, and crystallization paths likely affect their stability. For example, isobaric increases in temperature favor disordering and thus entropy-stabilization, in contrast, isothermal increases in pressure have the opposite effect. Although, currently, the potential effect of composition to cation disorder cannot be fully explored, it appears that enrichment in trivalent cations probably enhances entropy-stabilization and thus may increase the stability of (Al, Fe, Cr, Ti)-rich pseudobrookites relative to that of (Mg, Fe)-rich ones. In addition, high-temperature, nearly isothermal, decompression paths of olivine+orthopyroxene+oxide assemblages may favor pseudobrookites (s.l.) over rutile and/or ilmenite, in contrast, cooling at low pressures seems to favor ilmenite and/or rutile. In the case of crustal and mantle xenoliths, the presence or absence of orthorhombic oxides is probably controlled by reactions with olivine, orthopyroxene, ilmenite, and rutile. In oceanic mantle xenoliths such reactions may also involve a TiO2-enriched but not SiO2-enriched melt/fluid, because pseudobrookites (s.l.) would react with the SiO2-enriched melt/fluid to form orthopyroxene and rutile. Parenthetically, experiments and model calculations in the Mg-Ti-Si-O system suggest that low degree partial melting of low-TiO2 bulk compositions may produce Ti-enriched liquids in equilibrium with olivine, orthopyroxen ad=nd MgTi2O5, rutile or ilmenite.

Xirouchakis, Dimistrios↗

Experimental Characterization of an Electric Ducted Fan for the SUSAN 25% Scale Research Vehicle

An electric ducted fan (EDF) was tested in the NASA Langley 12-Foot Low-Speed Tunnel for the SUbsonic Single Aft eNgine (SUSAN) Electrofan 25% scale research vehicle, which incorporates distributed electric fans along the wings and serves as a scaled electrified transport aircraft testbed. The commercial-off-the-shelf EDF was selected as a propulsor that can meet the 25% scale vehicle geometry constraints, along with the performance requirements throughout the flight profile. The test included a series of one-factor-at-a-time (OFAT) and design of experiments (DOE) approaches that were implemented across a range of fan speed, voltage, freestream dynamic pressure, and incidence angle settings to characterize the performance of the EDF. Several test matrix designs and modeling approaches were considered. System identification methods were used to develop a static motor model, a dynamic motor model, an EDF thrust model, and a power consumption model. The model predictive capabilities were then evaluated with validation data not used in the modeling process. The results provide a comprehensive set of models developed to characterize the fan, and offer general insight into EDF modeling approaches.

electric ducted fan↗

Efficient Unitary Designs from Random Sums and Permutations

A unitary k-design is an ensemble of unitaries that matches the first k moments of the Haar measure. In this work, we provide two efficient constructions of k-designs on n-qubits using new random matrix theory techniques. Our first construction is based on exponentiating sums of random i.i.d. Hermitian matrices and uses O(k2n2)-many gates. In the spirit of central limit theorems, we show that this random sum approximates the Gaussian Unitary Ensemble (GUE). We then show that the product of just two exponentiated GUE matrices is already approximately Haar random. Our second construction is based on products of exponentiated sums of random permutations and uses Õ(k poly (n)) many gates. The k dependence is optimal (up to polylogarithmic factors) and is inherited from the efficiency of existing k-wise independent permutations. Furthermore, replacing random permutations with quantum-secure pseudorandom permutations (PRPs), we also obtain a pseudorandom unitary (PRU) ensemble that is secure under nonadaptive queries. A central feature of both proofs is a new connection between the polynomial method in quantum query complexity and the large-dimension (N) expansion in random matrix theory. In particular, the first construction uses the polynomial method to control high moments of certain random matrix ensembles without requiring delicate Weingarten calculations. In doing so, we define and solve a moment problem on the unit circle, asking whether a finite number of equally weighted points can reproduce a given set of moments. In our second construction, the key step is to exhibit an orthonormal basis for irreducible representations of the partition algebra that has a low-degree large-N expansion. This allows us to show that the distinguishing probability is a low-degree rational polynomial of the dimension N.

algebra↗

Microbe-Encapsulated Silica Gel Biosorbents for Selective Extraction of Scandium from Coal Byproducts

Establishing an inexpensive and environmentally friendly scandium (Sc) supply is critical for the development of clean energy technologies. However, refining Sc from Sc-bearing sources using current technologies poses economic and environmental challenges due to its low concentration relative to base metals and chemically similar lanthanides. In this work, we developed a biosorption-based flow-through process for extraction of Sc from low-grade feedstocks. A novel biosorbent material was synthesized by encapsulating Arthrobacter nicotianae, a bacterium that selectively adsorbs Sc, within a porous silica matrix. Through the use of several complementary microscopy techniques, we demonstrated a highly porous biosorbent structure with a high cell-loading density and a homogenous cell distribution. Batch adsorption assays revealed selective Sc adsorption over lanthanides and common based metals, with the exception of Fe(III), as evidenced by separation factors greater than 50. The MESG particles were packed into fixed-bed columns to enable adsorption under flow-through conditions, which demonstrated effective Sc extraction at flow rates up to 0.08 cm/s and high stability for reuse; greater than 95% of the adsorption capacity was maintained after 10 consecutive adsorption/desorption cycles. When applied to a lignite coal leachate, whose Fe content was depleted through pH-mediated precipitation, the MESG particles yielded Sc breakthrough at 30 bed volumes, whereas all other metals broke through after only a few bed volumes. Following desorption, an eluate with a 124-fold increase in Sc purity was achieved relative to the lignite leachate with Sc constituting 96.4% of the total REEs. This study established a rapid, facile, and scalable cell immobilization approach that enables the use of microbial biomass for Sc recovery from low-grade sources.

Dong, Ziye↗

Colossal piezoresistance in narrow-gap Eu 5 In 2 Sb 6

Piezoresistance, the change of the electrical resistance (R) of a material in response to an applied mechanical stress (σ), is the driving principle of electromechanical devices such as strain gauges, accelerometers, and cantilever force sensors. Enhanced piezoresistance has been traditionally observed in two classes of uncorrelated materials: nonmagnetic semiconductors and composite structures. In this study, we report the discovery of a remarkably large piezoresistance in Eu 5 In 2 Sb 6 single crystals, wherein anisotropic metallic clusters naturally form within a semiconducting matrix due to electronic interactions. Eu 5 In 2 Sb 6 shows a highly anisotropic piezoresistance, and uniaxial pressure along [001] of only 0.4 GPa leads to a resistivity drop of >99.95%, which results in a colossal piezoresistance factor of 5000×10 –11 Pa –1 . Our result not only reveals the role of interactions and phase separation in the realization of colossal piezoresistance, but it also highlights a route to multifunctional devices with large responses to both pressure and magnetic field.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Power counting to saturation

We present a description of saturation in small 𝑥 deep inelastic scattering from power counting in a top-down effective theory derived from quantum chromodynamics (QCD). A factorization formula isolates the universal physics of the nucleus at leading power in 𝑥. The onset of saturation is then understood as a breakdown in the expansion in an emergent power counting parameter, which is defined by the matrix element of a gauge invariant operator. We identify a new radiation mode, which enables us to extend previous literature by distinguishing the appearance of the saturation scale from the transition to nonlinear evolution.

Astronomy & Astrophysics↗

Optimizing Communication in 2D Grid-Based MPI Applications at Exascale

The new reality of exascale computing faces many challenges in achieving optimal performance on large numbers of nodes. A key challenge is the efficient utilization of the message-passing interface (MPI), a critical component for process communication. This paper explores communication optimization strategies to harness the GPU-accelerated architectures of these supercomputers. We focus on MPI applications where processors form a two-dimensional process grid, a common arrangement in applications involving dense matrix operations. This configuration offers a unique opportunity to implement innovative strategies to improve performance and maintain effective load distribution. We study two applications— Dist-FW (Apsp:all-pair-shortest-path) and HPL-MxP (LU factorization with Mixed precision)—on two accelerated systems: Summit (IBM Power and NVIDIA V100) and Frontier (AMD EPYC and MI250X). These supercomputers are operated by the Oak Ridge Leadership Computing Facility (OLCF) and are currently ranked #1 and #5 on the Top500 list. We show how to scale up both applications to exascale levels and tackle the MPI challenges related to implementation, synchronization, and performance. We also compare the performance of several communication strategies at an unprecedented scale. Accurately predicting application performance becomes crucial for cost reduction as the computational scale grows. To address this, we suggest a hyperbolic model as a better alternative to the traditional one-sided asymptotic model for predicting future application performance at such large scales.

Lu, Hao↗

A numerical study of bounds in the correlations of fractional quantum Hall states

We numerically compute the guiding center static structure factor \bar S(k) S ‾ ( k ) of various fractional quantum Hall (FQH) states to \mathcal{O}(k\ell)^6 𝒪 ( k ℓ ) 6 where k k is the wavenumber and \ell ℓ is the magnetic length. Employing density matrix renormalization group on an infinite cylinder of circumference L_y L y , we study the two-dimensional limit using L_y/\xi \gg 1 L y / ξ ≫ 1 , where \xi ξ is the correlation length. The main findings of our work are: 1) the ground states that deviate away from the ideal conformal block wavefunctions, do not saturate the Haldane bound, and 2) the coefficient of O(k\ell)^6 O ( k ℓ ) 6 term appears to be bounded above by a value predicted by field theories proposed in the literature. The first finding implies that the graviton mode is not maximally chiral for experimentally relevant FQH states.

Kumar, Prashant (ORCID:0000000258002768)↗

Learning an Algebriac Multrigrid Interpolation Operator Using a Modified GraphNet Architecture

This work, building on previous efforts, develops a suite of new graph neural network machine learning architectures that generate data-driven prolongators for use in Algebraic Multigrid (AMG). Algebraic Multigrid is a powerful and common technique for solving large, sparse linear systems. Its effectiveness is problem dependent and heavily depends on the choice of the prolongation operator, which interpolates the coarse mesh results onto a finer mesh. Previous work has used recent developments in graph neural networks to learn a prolongation operator from a given coefficient matrix. In this paper, we expand on previous work by exploring architectural enhancements of graph neural networks. A new method for generating a training set is developed which more closely aligns to the test set. Asymptotic error reduction factors are compared on a test suite of 3-dimensional Poisson problems with varying degrees of element stretching. Results show modest improvements in asymptotic error factor over both commonly chosen baselines and learning methods from previous work.

97 MATHEMATICS AND COMPUTING↗

Cratering and penetration experiments in teflon targets at velocities from 1 to 7 km/s

Approximately 20 sq m of protective thermal blankets, largely composed of Teflon, were retrieved from the Long Duration Exposure Facility after the spacecraft spent approximately 5.7 years in space. Examination of these blankets revealed that they contained thousands of hypervelocity impact features ranging from micron-sized craters to penetration holes several millimeters in diameter. We conducted impact experiments to reproduce such features and to understand the relationships between projectile size and the resulting crater or penetration hole diameter over a wide range of impact velocities. Such relationships are needed to derive the size and mass frequency distribution and flux of natural and man-made particles in low-earth orbit. Powder propellant and light-gas guns were used to launch soda-lime glass spheres into pure Teflon targets at velocities ranging from 1 to 7 km/s. Target thickness varied over more than three orders of magnitude from finite halfspace targets to very thin films. Cratering and penetration of massive Teflon targets is dominated by brittle failure and the development of extensive spall zones at the target's front and, if penetrated, the target's rear side. Mass removal by spallation at the back side of Teflon targets may be so severe that the absolute penetration hole diameter can become larger than that of a standard crater. The crater diameter in infinite halfspace Teflon targets increases, at otherwise constant impact conditions, with encounter velocity by a factor of V (exp 0.44). In contrast, the penetration hole size in very thin foils is essentially unaffected by impact velocity. Penetrations at target thicknesses intermediate to these extremes will scale with variable exponents of V. Our experimental matrix is sufficiently systematic and complete, up to 7 km/s, to make reasonable recommendations for velocity-scaling of Teflon craters and penetrations. We specifically suggest that cratering behavior and associated equations apply to all impacts in which the shock-pulse duration of the projectile is shorter than that assigned a unique projectile size, provided an impact velocity is known or assumed. This calibration seems superior to the traditional ballistic-limit approach.

Horz, Friedrich↗

Durability Challenges for Next Generation of Gas Turbine Engine Materials

Aggressive fuel burn and carbon dioxide emission reduction goals for future gas turbine engines will require higher overall pressure ratio, and a significant increase in turbine inlet temperature. These goals can be achieved by increasing temperature capability of turbine engine hot section materials and decreasing weight of fan section of the engine. NASA is currently developing several advanced hot section materials for increasing temperature capability of future gas turbine engines. The materials of interest include ceramic matrix composites with 1482 - 1648 C temperature capability, advanced disk alloys with 815 C capability, and low conductivity thermal barrier coatings with erosion resistance. The presentation will provide an overview of durability challenges with emphasis on the environmental factors affecting durability for the next generation of gas turbine engine materials. The environmental factors include gaseous atmosphere in gas turbine engines, molten salt and glass deposits from airborne contaminants, impact from foreign object damage, and erosion from ingestion of small particles.

Misra, Ajay K.↗