Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Matrix completion”

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 487 records · Page 27

Calculation of compressible flow about three-dimensional inlets with auxiliary inlets, slats and vanes by means of a panel method

An efficient and user-oriented method has been constructed for calculating flow in and about complex inlet configurations. Efficiency is attained by: the use of a panel method, a technique of superposition for obtaining solutions at any inlet operating condition, and employment of an advanced matrix-iteration technique for solving large full systems of equations, including the nonlinear equations for the Kutta condition. User concerns are addressed by the provision of several novel graphical output options that, taken together, yield a more complete comprehension of the flowfield than had been possible previously. Examples of these features are presented for some complicated configurations, and where possible, comparisons are made between calculation and experiment.

Hess, J. L.↗

Eigenstructure assignment for a thrust-vectored high angle-of-attack aircraft

Eigenstructure assignment is utilized to design flight control laws for a thrust-vectored aircraft at several different angles of attack. An interesting characteristic of the aircraft model is that the control distribution matrix is rank-deficient. Also, the effectiveness of the control inputs varies with the angle of attack. A pseudocontrol strategy is used to reduce the control space to two dimensions. After the eigenstructure assignment design is complete, the controller is mapped back to the original five-dimensional control space. The designs are shown to exhibit acceptable multivariable stability margins at the aircraft inputs.

Sobel, Kenneth M.↗

Improving the Performance of DGEMM with MoA and Cache-Blocking: Preprint

The goal of this paper is to demonstrate performance enhancements of the high performance dense linear algebra matrix-matrix multiply DGEMM kernel, widely implemented by vendors in the basic linear algebra subroutine BLAS library. The mathematics of arrays (MoA) paradigm due to Mullin (1988) results in contiguous memory accesses in combination with Church-Rosser complete language constructs optimized for target processor architectures [3]. Our performance studies demonstrate that the MoA implementation of DGEMM combined with optimal cache-blocking strategies results in at least a 25% performance gain on both Intel Xeon Skylake and IBM Power-9 processors over the vendor supplied Intel MKL and IBM ESSL basic linear algebra libraries. Results are presented for the NREL Eagle and ORNL Summit supercomputers.

cache-blocking↗

A summary of recent developments in transportation hazard classification activities for ammonium perchlorate

The transportation hazard classification of Ammonium Perchlorate is discussed. A test program was completed and data were forwarded to retain a Class 5.1 designation (oxidizer) for AP which is shipped internationally. As a follow-on to the initial team effort to conduct AP tests existing data were examined and a matrix which catalogs test parameters and findings was compiled. A collection of test protocols is developed to standardize test methods for energetic materials of all types. The actions to date are summarized; the participating organizations and their roles as presently understood; specific findings on AP (matrix); and issues, lessons learned, and potential actions of particular interest to the propulsion community which may evolve as a result of future U.N. propellant transportation classification activities.

Koller, A. M., Jr.↗

An improved model for the Earth's gravity field

An improved model for the Earth's gravity field, TEG-1, was determined using data sets from fourteen satellites, spanning the inclination ranges from 15 to 115 deg, and global surface gravity anomaly data. The satellite measurements include laser ranging data, Doppler range-rate data, and satellite-to-ocean radar altimeter data measurements, which include the direct height measurement and the differenced measurements at ground track crossings (crossover measurements). Also determined was another gravity field model, TEG-1S, which included all the data sets in TEG-1 with the exception of direct altimeter data. The effort has included an intense scrutiny of the gravity field solution methodology. The estimated parameters included geopotential coefficients complete to degree and order 50 with selected higher order coefficients, ocean and solid Earth tide parameters, Doppler tracking station coordinates and the quasi-stationary sea surface topography. Extensive error analysis and calibration of the formal covariance matrix indicate that the gravity field model is a significant improvement over previous models and can be used for general applications in geodesy.

Tapley, B. D.↗

Containerless low gravity processing of glass forming and immiscible alloys

Under normal one-g conditions immiscible alloys segregate extensively during solidification due to sedimentation of the more dense of the immiscible liquid phases. Immiscible (hypermonotectic) gold-rhodium alloys were processed in the 100 meter drop tube under low gravity, containerless conditions to determine the feasibility of producing dispersed structures. Three alloy compositions were utilized. Alloys containing 10 percent by volume of the gold-rich hypermonotectic phase exhibited a tendency for the gold-rich liquid to wet the outer surface of the samples. This wetting tendency led to extensive segregation in several cases. Alloys containing 80 and 90 percent by volume of the gold-rich phase possessed completely different microstructures from the 10 percent samples when processed under low-g, containerless conditions. Several samples exhibited microstructures consisting of well dispersed 2 to 3 microns diameter rhodium-rich spheres in a gold-rich matrix.

Andrews, J. Barry↗

Transition operators in electromagnetic-wave diffraction theory - General theory

A formal theory is developed for the scattering of time-harmonic electromagnetic waves from impenetrable immobile obstacles with given linear, homogeneous, and generally nonlocal boundary conditions of Leontovich (impedance) type for the wave of the obstacle's surface. The theory is modeled on the complete Green's function and the transition (T) operator in time-independent formal scattering theory of nonrelativistic quantum mechanics. An expression for the differential scattering cross section for plane electromagnetic waves is derived in terms of certain matrix elements of the T operator for the obstacle.

Hahne, G. E.↗

An algorithm for the basis of the finite Fourier transform

The Finite Fourier Transformation matrix (F.F.T.) plays a central role in the formulation of quantum mechanics in a finite dimensional space studied by the author over the past couple of decades. An outstanding problem which still remains open is to find a complete basis for F.F.T. In this paper we suggest a simple algorithm to find the eigenvectors of F.T.T.

Santhanam, Thalanayar S.↗

Utilization of Hemp Processing Waste for 3D Printing of Biocomposites

Unlike stem biomass, the residues after the extraction of cannabidiol (CBD) oil from hemp flower are challenging to utilize because of their high extractive content (~ 40%, mainly lipids) and are typically considered waste and landfilled. This study presented a novel approach to effectively valorize this underutilized hemp processing waste via chemical processing for three-dimensional (3D) printing applications. Hemp processing waste was processed with sodium hydroxide (NaOH) to control extractives for solving nozzle clogging and then applied as a biofiller in polylactic acid (PLA) composites to improve the mechanical strength. The novelty of this work lies in demonstrating that controlled extractive removal via NaOH treatment not only improves processability but also enhances mechanical performance in 3D-printed biocomposites. We systematically investigated the effect of the processed biofiller content (2.5–10 wt%) on the mechanical and thermal properties of the biocomposites. The decrease in the content of extractives reduced the non-structural components and improved the surface compatibility of the hemp waste with the PLA matrix, thereby enhancing the polymer-biofiller interactions. The best performance was achieved at 2.5 wt% loading, where Young’s modulus increased from 2.3 GPa to 2.6 GPa and tensile strength from 42.7 MPa to 48.8 MPa. Interestingly, the complete removal of extractives also reduced the mechanical strength of their biocomposites, indicating the interfacial adhesion effects of extractives. Furthermore, this study provides new insights into balancing extractive content for optimal mechanical properties, offering a sustainable solution for waste valorization in additive manufacturing.

Additive manufacturing↗

ZeroCAL: Eliminating Carbon Dioxide Emissions from Limestone’s Decomposition to Decarbonize Cement Production

Limestone (calcite, CaCO 3 ) is an abundant and cost-effective source of calcium oxide (CaO) for cement and lime production. However, the thermochemical decomposition of limestone (~800 °C, 1 bar) to produce lime (CaO) results in substantial carbon dioxide (CO 2(g) ) emissions and energy use, i.e., ~1 tonne [t] of CO 2 and ~1.4 MWh per t of CaO produced. Here, we describe a new pathway to use CaCO 3 as a Ca source to make hydrated lime (portlandite, Ca(OH) 2 ) at ambient conditions (p, T) while nearly eliminating process CO 2(g) emissions (as low as 1.5 mol. % of the CO 2 in the precursor CaCO 3 , equivalent to 9 kg of CO 2(g) per t of Ca(OH) 2 ) within an aqueous flowelectrolysis/ pH-swing process that coproduces hydrogen (H 2(g) ) and oxygen (O 2(g) ). Because Ca(OH) 2 is a zero-carbon precursor for cement and lime production, this approach represents a significant advancement in the production of zero-carbon cement. The Zero CArbon Lime (ZeroCAL) process includes dissolution, separation/recovery, and electrolysis stages according to the following steps: (Step 1) chelator (e.g., ethylenediaminetetraacetic acid, EDTA)-promoted dissolution of CaCO 3 and complexation of Ca 2+ under basic (>pH 9) conditions, (Step 2a) Ca enrichment and separation using nanofiltration (NF), which allows separation of the Ca-EDTA complex from the accompanying bicarbonate (HCO 3 – ) species, (Step 2b) acidity-promoted decomplexation of Ca from EDTA, which allows near-complete chelator recovery and the formation of a Ca-enriched stream, and (Step 3) rapid precipitation of Ca(OH) 2 from the Ca-enriched stream using electrolytically produced alkalinity. These reactions can be conducted in a seawater matrix yielding coproducts including hydrochloric acid (HCl) and sodium bicarbonate (NaHCO 3 ), resulting from electrolysis and limestone dissolution, respectively. Careful analysis of the reaction stoichiometries and energy balances indicates that approximately 1.35 t of CaCO 3 , 1.09 t of water, 0.79 t of sodium chloride (NaCl), and ~2 MWh of electrical energy are required to produce 1 t of Ca(OH) 2 , with significant opportunity for process intensification. This approach has major implications for decarbonizing cement production within a paradigm that emphasizes the use of existing cement plants and electrification of industrial operations, while also creating approaches for alkalinity production that enable cost-effective and scalable CO 2 mineralization via Ca(OH) 2 carbonation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Polarimetric signatures of a canopy of dielectric cylinders based on first and second order vector radiative transfer theory

Complete polarimetric signatures of a canopy of dielectric cylinders overlying a homogeneous half space are studied with the first and second order solutions of the vector radiative transfer theory. The vector radiative transfer equations contain a general nondiagonal extinction matrix and a phase matrix. The energy conservation issue is addressed by calculating the elements of the extinction matrix and the elements of the phase matrix in a manner that is consistent with energy conservation. Two methods are used. In the first method, the surface fields and the internal fields of the dielectric cylinder are calculated by using the fields of an infinite cylinder. The phase matrix is calculated and the extinction matrix is calculated by summing the absorption and scattering to ensure energy conservation. In the second method, the method of moments is used to calculate the elements of the extinction and phase matrices. The Mueller matrix based on the first order and second order multiple scattering solutions of the vector radiative transfer equation are calculated. Results from the two methods are compared. The vector radiative transfer equations, combined with the solution based on method of moments, obey both energy conservation and reciprocity. The polarimetric signatures, copolarized and depolarized return, degree of polarization, and phase differences are studied as a function of the orientation, sizes, and dielectric properties of the cylinders. It is shown that second order scattering is generally important for vegetation canopy at C band and can be important at L band for some cases.

Tsang, Leung↗

Spacecraft Formation Flying Maneuvers Using Linear Quadratic Regulation With No Radial Axis Inputs

Regarding multiple spacecraft formation flying, the observation has been made that control thrust need only be applied coplanar to the local horizon to achieve complete controllability of a two-satellite (leader-follower) formation. A formulation of orbital dynamics using the state of one satellite relative to another is used. Without the need for thrust along the radial (zenith-nadir) axis of the relative reference frame, propulsion system simplifications and weight reduction may be accomplished. This work focuses on the validation of this control system on its own merits, and in comparison to a related system which does provide thrust along the radial axis of the relative frame. Maneuver simulations are performed using commercial ODE solvers to propagate the Keplerian dynamics of a controlled satellite relative to an uncontrolled leader. These short maneuver simulations demonstrate the capacity of the controller to perform changes from one formation geometry to another. Control algorithm performance is evaluated based on measures such as the fuel required to complete a maneuver and the maximum acceleration required by the controller. Based on this evaluation, the exclusion of the radial axis of control still allows enough control authority to use Linear Quadratic Regulator (LQR) techniques to design a gain matrix of adequate performance over finite maneuvers. Additional simulations are conducted including perturbations and using no radial control inputs. A major conclusion presented is that control inputs along the three axes have significantly different relationships to the governing orbital dynamics that may be exploited using LQR.

Starin, Scott R.↗

Constructing a Simulation Surrogate with Partially Observed Output

Gaussian process surrogates are a popular alternative to directly using computationally expensive simulation models. When the simulation output consists of many responses, dimension-reduction techniques are often employed to construct these surrogates. However, surrogate methods with dimension reduction generally rely on complete output training data. This article proposes a new Gaussian process surrogate method that permits the use of partially observed output while remaining computationally efficient. The new method involves the imputation of missing values and the adjustment of the covariance matrix used for Gaussian process inference. The resulting surrogate represents the available responses, disregards the missing responses, and provides meaningful uncertainty quantification. In conclusion, the proposed approach is shown to offer sharper inference than alternatives in a simulation study and a case study where an energy density functional model that frequently returns incomplete output is calibrated.

42 ENGINEERING↗

Supersymmetric Virasoro minimal strings

A random matrix model definition of a family of N = 1 supersymmetric extensions of the Virasoro minimal string of Collier, Eberhardt, Mühlmann, and Rodriguez is presented. An analysis of the defining string equations shows that the models all naturally have unambiguous nonperturbative completions, which are explicitly supplied by the double-scaled orthogonal polynomial techniques employed. Perturbatively, the multiloop correlation functions of the model define a special supersymmetric class of “quantum volumes,” generalizing the prototype case, some of which are computed. Published by the American Physical Society 2024

Astronomy & Astrophysics↗

Space-based Measurements of Neutron Lifetime: A Novel Approach to a Long-standing Problem

Free neutrons decay via the weak interaction with a mean lifetime of around 15 minutes. Knowledge of this lifetime is important as it provides constraints on the unitarity of the CKM matrix and is a key parameter for studies of Big-Bang nucleosynthesis. However, current laboratory measurements differ by as much as 5σ. We are investigating a new technique to measure neutron lifetime: space-based neutron spectroscopy. We will complete a proof-of-principle demonstration using data from NASA’s MESSENGER mission to Mercury.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

UAE6 - Wind Tunnel Tests Data - UAE6 - Sequence O - Raw Data

Sequence O: Sin AOA, Parked (P) This sequence was designed to quantify the blade 3-D unsteady aerodynamic response in the absence of rotational influences by varying blade pitch angle. This test sequence used an upwind, rigid turbine with a 0° cone angle. The wind speed ranged from 18.9 m/s to 39.3 m/s, and the yaw angle was 0° throughout the sequence. The rotor was parked with the instrumented blade fixed at 0° azimuth, and the rotor lock was installed (see Appendix A). Blade and probe pressure measurements were collected. The teeter dampers were replaced with rigid links, and these two channels were flagged as not applicable by setting the measured values in the data file to -99999.99 Nm. The teeter link load cell was pre-tensioned to 40,000 N. Sinusoidal angle-of-attack histories were designed to emulate 2-D S809 dynamic data obtained at the OSU wind tunnel (Reuss Ramsey et al. 1995) as well as 1P yawed flow angle-of-attack variations the turbine could be expected to encounter in routine operation. The reduced frequency (K), mean angle of attack (αm), and oscillation amplitude (αω) were computed for each of the five primary span locations. The file lengths varied according to the time required to complete 40 oscillations at the specified frequency. The file name convention used the sequence designation O, followed by two digits for the span location (30, 47, 80, or 95), followed by two digits representing the condition specified in the test matrix, followed by the repetition digit. Due to time constraints, the test conditions corresponding to the 63% span location were eliminated. Test points with a static pitch angle of 3° were collected throughout the sequence to ascertain instrumentation fidelity.

17 WIND ENERGY↗

Fluid-structure finite-element vibrational analysis

A fluid finite element has been developed for a quasi-compressible fluid. Both kinetic and potential energy are expressed as functions of nodal displacements. Thus, the formulation is similar to that used for structural elements, with the only differences being that the fluid can possess gravitational potential, and the constitutive equations for fluid contain no shear coefficients. Using this approach, structural and fluid elements can be used interchangeably in existing efficient sparse-matrix structural computer programs such as SPAR. The theoretical development of the element formulations and the relationships of the local and global coordinates are shown. Solutions of fluid slosh, liquid compressibility, and coupled fluid-shell oscillation problems which were completed using a temporary digital computer program are shown. The frequency correlation of the solutions with classical theory is excellent.

Feng, G. C.↗

Space fabrication: Graphite composite truss welding and cap forming subsystems

An automated beam builder for the fabrication of space structures is described. The beam builder forms a triangular truss 1.3 meters on a side. Flat strips of preconsolidated graphite fiber fabric in a polysulfone matrix are coiled in a storage canister. Heaters raise the material to forming temperature then the structural cap section is formed by a series of rollers. After cooling, cross members and diagonal tension cords are ultrasonically welded in place to complete the truss. The stability of fabricated structures and composite materials is also examined.

Jenkins, L. M.↗