Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “matrix multiplication”

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 883 records · Page 49

Identification and Evaluation of Integration and Cross Cutting Issues Across HRP Risks

The HRP Integrated Research Plan contains the research plans for the 32 risks requiring research to characterize and mitigate. These risks to human health and performance in spaceflight are identified by evidence and each one focuses on a single aspect of human physiology or performance. They are further categorized by aspects of the spaceflight environment, such as altered gravity or space radiation, that that play a major role in their likelihood and consequence. From its inception the "integrate" in the Research Plan has denoted the integrated nature of risks to human health and performance, the connectedness of physiological systems within the human body regardless of the spaceflight environment, and the integrated response of the human body to the spaceflight environment. Common characteristics of the spaceflight environment include altered gravity, atmospheres and light/dark cycles, space radiation, isolation, noise, and periods of high or low workload. Long term exposure to this unique environment produces a suite of physiological effects such as stress; vision, neurocognitive and anthropometric changes; circadian misalignment; fluid shifts, deconditioning; immune dysregulation; and altered nutritional requirements. Matrix diagraming was used to systematically identify, analyze and rate the many-to-many relationships between environmental characteristics and the suite of physiological effects. It was also to identify patterns in the relationships of common physiological effects to each other. Analyses of patterns or relationships in these diagrams help to identify issues that cut across multiple risks. Cross cutting issues benefit from a multidisciplinary approach that synthesizes concepts or data from two or more disciplines to identify and characterize risk factors or develop countermeasures relevant to multiple risks. They also help to illuminate possible problem areas that may arise when a countermeasure impacts risks other than those which it was developed to mitigate, or identify groupings of physiological changes that are likely to occur that may impact the overall risk posture.

Steinberg, S. L.↗

Parallel Computation of the Jacobian Matrix for Nonlinear Equation Solvers Using MATLAB

Demonstrating speedup for parallel code on a multicore shared memory PC can be challenging in MATLAB due to underlying parallel operations that are often opaque to the user. This can limit potential for improvement of serial code even for the so-called embarrassingly parallel applications. One such application is the computation of the Jacobian matrix inherent to most nonlinear equation solvers. Computation of this matrix represents the primary bottleneck in nonlinear solver speed such that commercial finite element (FE) and multi-body-dynamic (MBD) codes attempt to minimize computations. A timing study using MATLAB's Parallel Computing Toolbox was performed for numerical computation of the Jacobian. Several approaches for implementing parallel code were investigated while only the single program multiple data (spmd) method using composite objects provided positive results. Parallel code speedup is demonstrated but the goal of linear speedup through the addition of processors was not achieved due to PC architecture.

Rose, Geoffrey K.↗

CMPLE: Correlation Modeling to Decode Photosynthesis Using the Minorize–Maximize Algorithm

In plant genomic experiments, correlations among various biological traits (phenotypes) give new insights into how genetic diversity may have tuned biological processes to enhance fitness under diverse conditions. Consequently, knowing how the correlations are affected by genetic (G) and environmental (E) factors helps develop climate-resilient plants. However, the current literature lacks any method for assessing the effect of predictors on pairwise correlations among multiple phenotypes together with easily interpretable model parameters. To address this need, we propose to model pairwise correlations directly in terms of G and E and develop a computationally efficient inference procedure. Two major novelties in our methodology are (1) the use of a composite pairwise likelihood method to avoid the positive definiteness restriction on the correlation matrix and (2) the use of a novel Minorize–Maximize (MM) algorithm for the efficient estimation of a large number of parameters. The proposed method shows excellent numerical performance on synthetic datasets. Here, the analysis of the motivating data on cowpea reveals that the rates of solar energy storage by photosynthesis (the aggregate trait) are differentially affected by different genetic loci through two distinct processes: “photoinhibition” which results from photodamage caused by excess light, and “photoprotection” which protects plants from photodamage but also results in energy loss.

Correlation modeling↗

Specific Hardening Function Definition and Characterization of a Multimechanism Generalized Potential-based Viscoelastoplasticity Model

Given the previous complete-potential structure framework together with the notion of strain- and stress-partitioning in terms of separate contributions of several submechanisms (viscoelastic and viscoplastic) to the thermodynamic functions (stored energy and dissipation) a detailed viscoelastoplastic multimechanism characterization of a specific hardening functional form of the model is presented and discussed. TIMETAL 21S is the material of choice as a comprehensive test matrix, including creep, relaxation, constant strain-rate tension tests, etc. are available at various temperatures. Discussion of these correlations tests, together with comparisons to several other experimental results, are given to assess the performance and predictive capabilities of the present model particularly with regard to the notion of hardening saturation as well as the interaction of multiplicity of dissipative (reversible/irreversible) mechanisms.

Arnold, S. M.↗

Experimental Setup and Parameters for Testing Uncoated and EBC-Coated CMCs under Thermal Gradients Induced by Laser Heating and Backside Air Cooling

An experimental setup to induce through thickness thermal gradients in ceramic matrix composite (CMC) specimens without and with an applied environmental barrier coating (EBC) is described in detail. Components of the experimental set up include a laser for the heat source on the front of the specimen, a cooling air nozzle on the back of the specimen with pyrometers to measure temperatures on both sides, and a full-field thermal image camera to ensure uniformity of the temperature distribution in the gage section of the specimen. Test specimen is loaded in a servohydraulic system, which can superimpose either static or dynamic loads on the specimen to investigate creep and fatigue behaviors of CMCs under through thickness thermal gradient conditions. Mechanical strains on the front and backsides of the test specimen were measured during the thermal gradient tests with superimposed mechanical loading using two contact extensometers. Selected through thickness thermal gradient test results generated with the experimental setup under different nominal front side temperatures and multiple constant laser power levels on uncoated and EBC coated CMCs and procedures used to generate those results are described in this report. Two examples of tests conducted using the test system on uncoated CMCs under creep and sustained-peak, low-cycle fatigue with through thickness thermal gradient conditions are provided.

ceramic matrix composite, environmental barrier co↗

ACTS Ka-Band Earth Stations: Technology, Performance, and Lessons Learned

The Advanced Communications Technology Satellite (ACTS) Project invested heavily in prototype Ka-band satellite ground terminals to conduct an experiments program with the ACTS satellite. The ACTS experiment's program proposed to validate Ka-band satellite and ground station technology. demonstrate future telecommunication services. demonstrate commercial viability and market acceptability of these new services, evaluate system networking and processing technology, and characterize Ka-band propagation effects, including development of techniques to mitigate signal fading. This paper will present a summary of the fixed ground terminals developed by the NASA Glenn Research Center and its industry partners, emphasizing the technology and performance of the terminals (Part 1) and the lessons learned throughout their six year operation including the inclined orbit phase of operations (Full Report). An overview of the Ka-band technology and components developed for the ACTS ground stations is presented. Next. the performance of the ground station technology and its evolution during the ACTS campaign are discussed to illustrate the technical tradeoffs made during the program and highlight technical advances by industry to support the ACTS experiments program and terminal operations. Finally. lessons learned during development and operation of the user terminals are discussed for consideration of commercial adoption into future Ka-band systems. The fixed ground stations used for experiments by government, academic, and commercial entities used reflector based offset-fed antenna systems ranging in size from 0.35m to 3.4m antenna diameter. Gateway earth stations included two systems, referred to as the NASA Ground Station (NGS) and the Link Evaluation Terminal (LET). The NGS provides tracking, telemetry, and control (TT&C) and Time Division Multiple Access (TDMA) network control functions. The LET supports technology verification and high data rate experiments. The ground stations successfully demonstrated many services and applications at Ka-band in three different modes of operation: circuit switched TDMA using the satellite on-board processor, satellite switched SS-TDMA applications using the on-board Microwave Switch Matrix (MSM), and conventional transponder (bent-pipe) operation. Data rates ranged from 4.8 kbps up to 622 Mbps. Experiments included: 1) low rate (4.8- 1 00's kbps) remote data acquisition and control using small earth stations, 2) moderate rate (1-45 Mbps) experiments included full duplex voice and video conferencing and both full duplex and asymmetric data rate protocol and network evaluation using mid-size ground stations, and 3) link characterization experiments and high data rate (155-622 Mbps) terrestrial and satellite interoperability application experiments conducted by a consortium of experimenters using the large transportable ground stations.

Reinhart, Richard C.↗

A Generic Guidance and Control Structure for Six-Degree-of-Freedom Conceptual Aircraft Design

A control system framework is presented for both real-time and batch six-degree-of-freedom simulation. This framework allows stabilization and control with multiple command options, from body rate control to waypoint guidance. Also, pilot commands can be used to operate the simulation in a pilot-in-the-loop environment. This control system framework is created by using direct vehicle state feedback with nonlinear dynamic inversion. A direct control allocation scheme is used to command aircraft effectors. Online B-matrix estimation is used in the control allocation algorithm for maximum algorithm flexibility. Primary uses for this framework include conceptual design and early preliminary design of aircraft, where vehicle models change rapidly and a knowledge of vehicle six-degree-of-freedom performance is required. A simulated airbreathing hypersonic vehicle and a simulated high performance fighter are controlled to demonstrate the flexibility and utility of the control system.

Cotting, M. Christopher↗

A Generic Guidance and Control Structure for Six-Degree-of-Freedom Conceptual Aircraft Design

A control system framework is presented for both real-time and batch six-degree-of-freedom simulation. This framework allows stabilization and control with multiple command options, from body rate control to waypoint guidance. Also, pilot commands can be used to operate the simulation in a pilot-in-the-loop environment. This control system framework is created by using direct vehicle state feedback with nonlinear dynamic inversion. A direct control allocation scheme is used to command aircraft effectors. Online B-matrix estimation is used in the control allocation algorithm for maximum algorithm flexibility. Primary uses for this framework include conceptual design and early preliminary design of aircraft, where vehicle models change rapidly and a knowledge of vehicle six-degree-of-freedom performance is required. A simulated airbreathing hypersonic vehicle and a simulated high performance fighter are controlled to demonstrate the flexibility and utility of the control system.

Cotting, M. Christopher↗

Space Weathering of Genesis Mission Solar-wind Collectors with Inferences for Weathering on Airless Bodies

Samples from the Genesis Solar Wind Sample Return (NASA Discovery 5) are a unique opportunity to study the initial stages of space weathering, i.e., the physical and chemical effects of solar-wind irradiation. Arrays of collectors containing multiple materials were each exposed to a different solar-wind regime (fast, slow, bulk, or coronal mass ejection) at the L1 point for long durations (years). Materials exposed to the solar wind included metals, semiconductors, and insulators. Although the time of exposure was obviously short relative to samples having extraterrestrial origins, optical properties, surface chemical properties, and matrix structure have changed in many collectors due to exposure to solar wind. The thickness of amorphous zones, where present, appears to correspond with the depth of the peak of the solar-wind H distribution in each regime. Damage from high-energy particles was negligible because the collectors were 700 μm or less in thickness and shielded from the back by the spacecraft. Micrometeorite impacts and sputtering were also negligible because of the short exposure times. Our current results are preliminary: We hope future workers will extend this study both to support Genesis characterization efforts and to further understand space-weathering processes on a geologic timescale.

space weathering↗

Three-dimensional imaging of pion using lattice QCD: generalized parton distributions

In this work, we report a lattice calculation of x-dependent valence pion generalized parton distributions (GPDs) at zero skewness with multiple values of the momentum transfer −t. The calculations are based on an N f = 2 + 1 gauge ensemble of highly improved staggered quarks with Wilson-Clover valence fermion. The lattice spacing is 0.04 fm, and the pion valence mass is tuned to be 300 MeV. We determine the Lorentz-invariant amplitudes of the quasi-GPD matrix elements for both symmetric and asymmetric momenta transfers with similar values and show the equivalence of both frames. Then, focusing on the asymmetric frame, we utilize a hybrid scheme to renormalize the quasi-GPD matrix elements obtained from the lattice calculations. After the Fourier transforms, the quasi-GPDs are then matched to the light-cone GPDs within the framework of large momentum effective theory with improved matching, including the next-to-next-to-leading order perturbative corrections, and leading renormalon and renormalization group resummations. We also present the 3-dimensional image of the pion in impact-parameter space through the Fourier transform of the momentum transfer −t.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Quantitative and three-dimensional assessment of holdup material

Nuclear material deposited in equipment, transfer lines, and ventilation systems of a processing facility is usually referred to as holdup. Operators need to know the location and amount of holdup in their facility for radiation safety and material accountability. In this work, we propose to use an array of detectors co-axial to the inspected pipe to measure the holdup material. This method is implementable into an automated system capable of crawling on surfaces and pipes of various curvatures, which would enable faster, easier, and more accurate holdup safeguards measurements. We first demonstrated that the current holdup assay procedure could lead to a non-negligible bias in the estimate of special nuclear material mass, due to the simplified assumption of deposited geometry introduced by the Generalized Geometry Holdup (GGH) model. The new approach consists of imaging the inner holdup material by characterizing the detector array’s response and unfolding it from the measured light output. Our experimental proof of principle consists of three NaI(Tl) detectors surrounding an aluminum pipe containing two cesium-137 ( 137 Cs) sources. We derived the source distribution inside the pipe by first calculating the detector response matrix using a method adaptive to the surface geometry of the object containing the measured holdup material. Creating a matrix of the detector array’s measured counts, we then proceed to solve an inverse problem, resulting in an accurately located source position and activity distribution within the response matrix’s spatial resolution. Furthermore, we then developed a simulated model of the envisioned experimental setup, which accurately described both the activity and position of the source in 2D. Finally, we extended our model onto a discretized three-dimensional model of the system, encompassing 36 detectors. The model was simulated to measure the source distribution and activity of multiple sources in pipes of varied geometries, which could then be conveniently translated to a real-time holdup measurement protocol. For the 3D simulation of four different source geometries, the model accurately localized the source position in 3D, while the activity retained a maximum relative error of ±5.32%.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

A call to standardize metrics for monitoring baleen whales near marine construction activities

Effective monitoring is necessary to protect marine mammal species during the construction of offshore infrastructure. The tools for detecting or monitoring marine mammals span traditional (e.g., visual observers, optical cameras), to newer (e.g., passive acoustic monitoring, infrared cameras, tags), and emerging (e.g., satellite imagery, environmental DNA, dimethyl sulfide concentration) technologies. Some are better suited for use during offshore development; however, peer-reviewed literature does not typically evaluate and report on the performance of these various technologies. We define a minimum set of metrics related to efficacy (i.e., confusion matrix, precision and recall, probability of missed mitigation), detection range (i.e., maximum and reliable detection range, spatial resolution), and data delivery (i.e., detection latency, system reliability, temporal resolution) that we recommend are needed to assess the utility of monitoring technologies for this purpose. Following a literature review of relevant studies, we highlight which publications reported these metrics and used multiple technologies to compare relative performance. We also emphasize the benefits of multi-modal approaches and recommend performance assessments through modeling or large-scale collaborative field testing. These metrics will standardize data collection, reporting, and analysis; promote consistent and comparable results; and foster collaboration among developers, regulatory agencies, and scientists. This may lead to the co-development of technology that achieves multiple goals, has greater application, and can answer research questions while collecting data to fulfill permitting requirements. These metrics may also inform decisions on what systems regulatory agencies might consider using and reduce monitoring costs, which is critical to support the marine sector's rapid growth alongside marine mammal conservation.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Data assimilation in operator algebras

We develop an algebraic framework for sequential data assimilation of partially observed dynamical systems. In this framework, Bayesian data assimilation is embedded in a nonabelian operator algebra, which provides a representation of observables by multiplication operators and probability densities by density operators (quantum states). In the algebraic approach, the forecast step of data assimilation is represented by a quantum operation induced by the Koopman operator of the dynamical system. Moreover, the analysis step is described by a quantum effect, which generalizes the Bayesian observational update rule. Projecting this formulation to finite-dimensional matrix algebras leads to computational schemes that are i) automatically positivity-preserving and ii) amenable to consistent data-driven approximation using kernel methods for machine learning. Moreover, these methods are natural candidates for implementation on quantum computers. Applications to the Lorenz 96 multiscale system and the El Niño Southern Oscillation in a climate model show promising results in terms of forecast skill and uncertainty quantification.

97 MATHEMATICS AND COMPUTING↗

A Block-Based Triangle Counting Algorithm on Heterogeneous Environments

Triangle counting is a fundamental building block in graph algorithms. In this article, we propose a block-based triangle counting algorithm to reduce data movement during both sequential and parallel execution. Our block-based formulation makes the algorithm naturally suitable for heterogeneous architectures. The problem of partitioning the adjacency matrix of a graph is well-studied. Our task decomposition goes one step further: it partitions the set of triangles in the graph. By streaming these small tasks to compute resources, we can solve problems that do not fit on a device. We demonstrate the effectiveness of our approach by providing an implementation on a compute node with multiple sockets, cores and GPUs. The current state-of-the-art in triangle enumeration processes the Friendster graph in 2.1 seconds, not including data copy time between CPU and GPU. Using that metric, our approach is 20 percent faster. When copy times are included, our algorithm takes 3.2 seconds. This is 5.6 times faster than the fastest published CPU-only time.

97 MATHEMATICS AND COMPUTING↗

Core-flood Effluent and Shale Surface Chemistries in Predicting Interaction between Shale, Brine, and Reactive Fluid

We report field and laboratory observations to date indicate that the efficiency of hydraulic fracturing, as it relates to hydrocarbon recovery, depends significantly on geochemical alterations to rock surfaces that diminish accessibility by partial or total plugging of the pore and fracture networks. This is caused by mineral scale deposition such as coating of fracture surfaces with precipitates, particle migration, and/or crack closure due to dissolution under stress. In reactive flow-through experiments, mineral reactions in response to acidic fluid injection significantly reduced system porosity and core permeability. The present study focuses on changes to fluid chemistry and shale surfaces (inlet and fracture walls) resulting from shale-fluid interactions and integrating these findings for an improved estimate of transport-related consequences. The reacted shale surfaces were examined by spatially-resolved scanning electron microscopy - energy dispersive spectroscopy (SEM-EDS) analysis. Importantly, inductively coupled plasma - mass spectrometry/optical emission spectroscopy (ICP-MS/OES) was utilized to probe the chemical evolution of the core-flood effluents. The three study cores selected from the Marcellus formation represent different mineralogies and structural features. In flow-through experiments, lab-generated brine and HCl-based fracture fluid (pH=2) were injected sequentially under effective stress (up to 500 psi) at reservoir temperature (80°C). SEM-EDS results confirmed by the ICP concentration trends showed significant Fe hydroxide precipitates in clay- and pyrite-rich outcrop samples due to partial oxidation of Fe-bearing phases in the case of intrusion of low salinity water-based fluids. Porosity reduction in the MSEEL (Marcellus Shale Energy and Environmental Laboratory) carbonate-rich sample is related to compaction of cores under stress due to matrix softening with substantial dissolution, and pore-filling by hydroxides, as well as barite and salts. Despite the same fluid compositions and experimental conditions used for both MSEEL samples, barite precipitation was much more intense in the MSEEL clay-rich sample due to its greater sorption capacity and additional sulfate source as well as fissile nature with multiple lengthwise cracks. ICP tests revealed time-resolved concentration trends in produced brine and reactive fluids that in turn complemented the pre-/post-reaction SEM-EDS observations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Mechanistic Approach to Analyzing and Improving Unconventional Hydrocarbon Production [Slides]

DOE research is developing the physical basis and tools needed to manage pressure effectively to increase recovery efficiency. By coupling fast, accurate physics with machine learning, DOE is producing science-based platforms any operator can use. DOE’s research portfolio is targeting hydrocarbon transport at multiple scales, with the goal of increasing recovery efficiency. DOE’s research has led to new, fast & accurate platforms for predicting gas production from fractured shales. Using data from the MSEEL-I site to calibrate our physics-based model, we have early results on pressure management. We have shown that both mechanical and chemical processes in the matrix can negatively impact production.

08 HYDROGEN↗

A Block-Based Triangle Counting Algorithm on Heterogeneous Environments

Triangle counting is a fundamental building block in graph algorithms. In this paper, we propose a block-based triangle counting algorithm to reduce data movement during both sequential and parallel execution. Our block-based formulation makes the algorithm naturally suitable for heterogeneous architectures. The problem of partitioning the adjacency matrix of a graph is well-studied. Our task decomposition goes one step further: it partitions the set of triangles in the graph. By streaming these small tasks to compute resources, we can solve problems that do not fit on a device. We demonstrate the effectiveness of our approach by providing an implementation on a compute node with multiple sockets, cores and GPUs. The current state-of-the-art in triangle enumeration processes the Friendster graph in 2.1 seconds, not including data copy time between CPU and GPU. Using that metric, our approach is 20 percent faster. When copy times are included, our algorithm takes 3.2 seconds. This is 5.6 times faster than the fastest published CPU-only time.

97 MATHEMATICS AND COMPUTING↗

Versatile soil gas concentration and isotope monitoring: optimization and integration of novel soil gas probes with online trace gas detection

Abstract. Gas concentrations and isotopic signatures can unveil microbial metabolisms and their responses to environmental changes in soil. Currently, few methods measure in situ soil trace gases such as the products of nitrogen and carbon cycling or volatile organic compounds (VOCs) that constrain microbial biochemical processes like nitrification, methanogenesis, respiration, and microbial communication. Versatile trace gas sampling systems that integrate soil probes with sensitive trace gas analyzers could fill this gap with in situ soil gas measurements that resolve spatial (centimeters) and temporal (minutes) patterns. We developed a system that integrates new porous and hydrophobic sintered polytetrafluoroethylene (sPTFE) diffusive soil gas probes that non-disruptively collect soil gas samples with a transfer system to direct gas from multiple probes to one or more central gas analyzer(s) such as laser and mass spectrometers. Here, we demonstrate the feasibility and versatility of this automated multiprobe system for soil gas measurements of isotopic ratios of nitrous oxide (δ18O, δ15N, and the 15N site preference of N2O), methane, carbon dioxide (δ13C), and VOCs. First, we used an inert silica matrix to challenge probe measurements under controlled gas conditions. By changing and controlling system flow parameters, including the probe flow rate, we optimized recovery of representative soil gas samples while reducing sampling artifacts on subsurface concentrations. Second, we used this system to provide a real-time window into the impact of environmental manipulation of irrigation and soil redox conditions on in situ N2O and VOC concentrations. Moreover, to reveal the dynamics in the stable isotope ratios of N2O (i.e., 14N14N16O, 14N15N16O, 15N14N16O, and 14N14N18O), we developed a new high-precision laser spectrometer with a reduced sample volume demand. Our integrated system – a tunable infrared laser direct absorption spectrometry (TILDAS) in parallel with Vocus proton transfer reaction mass spectrometry (PTR-MS), in line with sPTFE soil gas probes – successfully quantified isotopic signatures for N2O, CO2, and VOCs in real time as responses to changes in the dry–wetting cycle and redox conditions. Broadening the collection of trace gases that can be monitored in the subsurface is critical for monitoring biogeochemical cycles, ecosystem health, and management practices at scales relevant to the soil system.

54 ENVIRONMENTAL SCIENCES↗