Search NASA⌕ Search

SEARCH · Search NASA

Results for “parallel data structures”

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 505 records · Page 28

Analysis Ready Data in Analytics Optimized Data Stores for Analysis of Big Earth Data in the Cloud

Cloud computing offers the possibility of making the analysis of Big Data approachable for a wider community due to affordable access to computing power, an ecosystem of usable tools for parallel processing, and migration of many large datasets to archives in the cloud, allowing data-proximal computing. Generally, data analysis acceleration in the cloud comes from running multiple nodes in a split-combine-apply strategy. Data systems such as the Earth Observing System Data and Information System are in a position to "pre-split" the data by storing them in a data store that is optimized for data parallel computing, i.e., an Analytics-Optimized Data Store (AODS). A variety of approaches to AODS are possible, from highly scalable databases to scalable filesystems to data formats optimized for cloud access (e.g., zarr and cloud-optimized datasets), with the optimal choice dependent on both the types of analysis and the geospatial structure of the data. A key question is how much preprocessing of the data to do, both before splitting and as the first part of the apply step. Again, the geospatial structure of the data and the analysis type influence the decision, with the added complexity of the user type. Trans-disciplinary users who are not well-versed in the nuances of quality-filtering and georeferencing of remote sensing orbit/swath/scene data tend to ask for more highly processed data, relying on the data provider to make sensible decisions on preprocessing parameters. (This accounts for the popularity of "Level 3" gridded data, despite the lower spatial resolution it provides.) In this case, data can be preprocessed before the split, resulting in higher performance in the rest of the "apply" step, which can be transformative for use cases such as interactive data exploration at scale. Discipline researchers who are experienced with remote sensing data often prefer more flexibility in customizing the preprocessing data into Analysis Ready Data, resulting in more need for on-the-fly preprocessing.

Lynnes, Christopher↗

Investigation of Benchmark $k$ eff Sensitivity and Uncertainty for 239 Pu fission in Specific Energy Ranges

Nuclear data at intermediate energies (from 1 to 100s of keV) are evaluated based on scarce differential data and theory unable to capture physics’ expected structure. There is also a lack of integral data. This is a known deficiency and is challenging to address. Calculated effective multiplication factor, k eff , values for intermediate energy experiments are ~25× further from experiment than for fast energies and are often well outside the experimental uncertainties. The goal of the PARADIGM (PARallel Approach of Differential and InteGral Measurements) project is to significantly re duce the uncertainties of intermediate energy nuclear data for 239 Pu. To this end, PARADIGM simultaneously optimizes experiments at both the Los Alamos Neutron Science Center (LANSCE) and National Criticality Experiments Research Center (NCERC). The combined set of data will inform new intermediate-energy nuclear data. By execution of differential and integral experiments, establishment of new theory, and undertaking nuclear data evaluation in parallel, the timeline to deliver improved nuclear data to users will be reduced significantly that is to three years. For the PARADIGM project, it was decided to optimize an integral experiment for two neutron energy ranges, within the full intermediate energy range. The low energy range goes from 1 to 30 keV, while the higher energy range goes from 30 to 600 keV. This work focuses on nuclear data sensitivities and uncertainties for 239 Pu fission for existing experiments in the International Criticality Safety Benchmark Evaluation Project (ICSBEP). When designing new experiments, it is important to understand what benchmarks currently exist. For a more traditional experiment design (in which a specific application model(s) exists), comparisons would be made between the application model(s) and existing benchmarks. For PARADIGM, there is no specific application model, but instead the specific nuclear data reaction and energy ranges of interest can be explored for existing benchmarks.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Chocolate tablet aspects of cytherean Meshkenet Tessera

Meshkenet Tessera structures were mapped from Magellan data and several resemblances to chocolate tablet boudinage were found. The complex fault sets display polyphase tectonic sequences of a few main deformation phases. Shear and tension have contributed to the areal deformation. Main faults cut the 1600-km long Meshkenet Tessera highland into bar-like blocks which have ridge and groove pattern oriented along or at high angles to the faults. The first approach to the surface block deformation is an assumption of initial parallel shear faulting followed by a chocolate tablet boudinage. Major faults which cut Meshkenet Tessera into rectangular blocks have been active repetitively while two progressive or superposed boudinage set formations have taken place at high angles during the relaxational or flattening type deformation of the area. Chocolate tablet boudinage is caused by a layer-parallel two-dimensional extension resulting in fracturing of the competent layer. Such structures, defined by two sets of boudin neck lines at right angles to each other, have been described by a number of authors. They develop in a flattening type of bulk deformation or during superposed deformation where the rock is elongated in two dimensions parallel to the surface. This is an attempt to describe and understand the formation and development of structures of Meshkenet Tessera which has complicated fault structures.

Raitala, J.↗

A physics-informed and hierarchically regularized data-driven model for predicting fluid flow through porous media

This paper presents a new deep learning data-driven model for predicting structure dependent pore-fluid velocity fields in rock. The model is based on a Convolutional Auto-Encoder (CAE) artificial neural network capable of learning from image data generated by direct numerical simulations of fluid flow through pore-structures, such as by Lattice Boltzmann or molecular dynamics methods. The main novelty of the model in comparison to previous CAE-based data-driven approaches consists of three parts. The first is a methodology for decomposing the full-domain of the porous media into sub-regions, or “sub-domains”, in order to reduce the overall size of the CAE, batch process the sub-domains in parallel, and enable the CAE to learn local and generalizable nonlinear mappings of pore-fluid velocities. The second consists of embedding the finite difference solutions of the incompressible Navier-Stokes and continuity equations into convolutional layers prior to the CAE in order to provide the CAE with knowledge of fluid dynamics physics (PhyFlow). The third main novelty is that the training of the CAE is regularized with a hierarchical loss function that encourages the learning of fluid flow patterns (in a way similar to ranked modes in principal component analysis), ranking from most to least important. This is shown to increase the stability in learning, reduce over-fitting, and promote interpretability of the CAE neural network layers (HierCAE). The comprehensive new data-driven model, which we call the PhyFlow-HierCAE model, is shown to exhibit improved accuracy and generalizability of flow field predictions over conventional CAE models, attributable to the embedded physical knowledge and the hierarchical regularization, as well as realize orders of magnitude speed-ups in computation times as a surrogate for the direct numerical simulations. Examples of training and forward predictions on unseen pore-structures are provided and evaluated for data from Lattice Boltzmann and molecular dynamics simulations of pore-fluid flow. The model is shown to be a fast and accurate emulator (or “surrogate”) for predicting effective permeability of unseen pore-structures based on learning from relatively small direct numerical simulation datasets.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

HIF-dependent and reversible nucleosome disassembly in hypoxia-inducible gene promoters

Highlights: • Nucleosome-free regions (NFRs) are established in some gene promoters in a hypoxia-dependent manner. • Hypoxia-inducible NFR (iNFR) formation requires hypoxia-inducible transcription factor (HIF) activity. • Reoxygenation induces the reassembly of nucleosome structures in iNFRs. • Nucleosome reassembly in iNFRs requires SIN3A. Hypoxia causes dramatic changes in gene expression profiles, and the mechanism of hypoxia-inducible transcription has been analyzed for use as a model system of stress-inducible gene regulation. In this study, changes in chromatin organization in promoters of hypoxia-inducible genes were investigated during hypoxia-reoxygenation conditions. Most of the hypoxia-inducible gene promoters were hypersensitive to DNase I under both normal and hypoxic conditions, and our data indicate an immediate recruitment of transcription factors under hypoxic conditions. In some of the hypoxia-inducible promoters, nucleosome-free DNA regions (NFRs) were established in parallel with hypoxia-induced transcription. We also show that the hypoxia-inducible formation of NFRs requires that hypoxia-inducible transcription factors (HIFs) bind to the promoters together with the transcriptional coactivator CBP. Within 1 h after the hypoxia exposure was ended (reoxygenation), HIF complexes were dissociated from the promoter regions. Within 24 h of reoxygenation, the hypoxia-induced transcription returned to basal levels and the nucleosome structure was reassembled in the hypoxia-inducible NFRs. Nucleosome reassembly required the function of the transcriptional coregulator SIN3A. Thus, reversible changes in nucleosome organization mediated by transcription factors are notable features of stress-inducible gene regulation.

60 APPLIED LIFE SCIENCES↗

Development of NDE/NDT Tools for High-Volume & High-Speed Inspection of CFRP Structures in Automotive Manufacturing

Main advantages of the air-coupled ultrasound testing (ACUT) and electromagnetic testing (EMT) techniques for NDE of CFRP composites were non-contact sensing, scalability for high-speed inspection, cost-effectiveness, and non-hazardous operation. Despite these advantages, no systems that would satisfy the project requirements were commercially available. Hence, one of the major efforts of the Michigan State University (MSU) team at the initial stage of the project was to close this technological gap by developing, optimizing, and validating array sensors that would provide sufficient sensitivity, spatial coverage, and resolution for robust defect detection. Optimization of the ACUT and EMT sensor designs was performed using experimentally validated finite element models. Initial experiments using array probes were conducted on relatively flat CFRP samples. In parallel, the MSU team designed and assembled a portable platform with two robotic arms. The robots were equipped with newly designed sensors that enabled high-speed NDE of curved CFRP parts. Presently, the developed robotic platform can be used as a demo/template NDE system, which is easily adaptable to manufacturing environments and in-line NDE. The ACUT NDE system developed by the MSU team used a high-power 4-channel pulser receiver for parallel data acquisition. The array probes were designed by stacking commercially available ACUT transducers, which operated in the frequency range between 100 kHz and 500 kHz. MSU optimized the excitation procedure and developed wave focusing cones so as to reduce the crosstalk between the transducers and to provide higher pulse repletion frequency (PRF). The through-transmission (TT) and single-side access (SSA) inspection modes were successfully implemented. In the TT-ACUT, structural defects in CFRP were detected by passing ultrasonic waves through the test part. Hence, the ACUT transmitters and receivers needed to be placed on the opposite sides of the test part. In the SSA-ACUT, guided waves (GW) were excited in the test part using the transmitters and were sensed by the receivers from the same side. Multi-channel TT-ACUT and SSA-ACUT provided high-speed NDE, and were successfully validated on CFRP test samples with interlaminar delaminations and other embedded defects The EM techniques developed by the MSU team included: 1) eddy current testing (ECT), 2) capacitive imaging (CI) and hybrid dual-mode imaging. In ECT, structural damage was detected in CFRP using coils sensor arrays. In ECT, the excitation magnetic field is generated by passing an alternating current through a coil, which is placed above the test sample. The excitation field penetrates the conductive sample and induces the eddy currents in its transect. In turn, the eddy currents generate the reaction field, which affects the total field sensed by a coil. Hence, the presence of structural flaws will alter the eddy current flow and the picked-up signal. ECT is mostly sensitive to local changes of the electric conductivity of the test sample, and CFRPs are mostly conductive in the direction of carbon fibers. Hence, ECT was well suited for the detection of fiber damage/fiber irregularities. The MSU team developed printed circuit boards (PCB) with coil sensor arrays optimized for NDE of CFRP. Unlike most commercial probes designed for ECT of metallic structures, the MSU array probes were designed for operation in [1-10] MHz frequency range, which was optimal for low-conductive CFRP. Multiple sensing topologies (coil groups excitation/sensing arrangements) were implemented and successfully validated. Capacitive Imaging (CI) technique developed by MSU was complementary to ECT. In contrast to ECT, which was sensitive to local changes of the electrical conductivity, the CI was sensitive to local changes of the dielectric constant. Therefore, CI could provide information about matrix damage/matrix irregularities in CFRP. The MSU CI sensor arrays were made of multiple circular or rectangular open-plate capacitors printed on PCB. Sensors of this type are not commercially available. In addition to ECT and CI, the MSU team developed a hybrid (dual-mode) inductive/capacitive measurement technique that synergistically combined the benefits of inductive and capacitive sensing for rapid NDE of fiber reinforced polymer (FRP) composite structures. Fiber damage and fiber irregularities in FRPs were detected by configuring hybrid sensors as coil sensors. Similarly, matrix damage, matrix irregularities and interlaminar delaminations were detected by configuring hybrid sensors as capacitive sensors. ECT and CI were performed sequentially by means of electronic switching. Hence, eliminating the need for mounting two separate sensor arrays on the probe. Portable robotic platform was developed by MSU for multi-technique high-speed NDE of CFRP test parts. The platform had two 6-axis robots, which enabled inspection of curved parts in approximately a 6×6×6 ft 3 active scan area. On the software side, the MSU team integrated scripts for NDE hardware control with scripts for robot motion control. MSU also implemented automated path planning for the robots, reconstruction of part’s surfaces via stereovision, 3D rendering of inspection data, and image processing algorithms for enhanced defect detection. Automotive composite parts manufactured by Plasan Composites from Phase I were used to validate the ACUT and EMT techniques on representative testbeds. Among those parts were three X-braces for a Dodge Viper, one composite calibration plaque with known defects at known locations, and four other test sections, including sections from a front splitter, a corner section from a composite hood, and a high-pressure RTM panel made using non crimp fabric. Other test samples included CFRP and GFRP calibration plates with fiber/matrix defects fabricated at MSU/CVRC.

36 MATERIALS SCIENCE↗

Large-Scale Welding Process Simulation by GPU Parallelized Computing

The computational design of industrially relevant welded structures is extremely time consuming due to coupled physics and high nonlinearity. Previously, most welding distortion and residual stress simulations have been limited to small coupons and reduced order (from three-dimensional [3D] to two-dimensional [2D]), or inherent strain approximations were used for large structures. In this current study, an explicit finite element code based on a graphics processing unit was utilized to perform 3D transient thermomechanical simulation of structural components during welding. Laser brazing of aluminum alloy panels as representative of automotive manufacturing scenarios was simulated to predict out-of-plane distortion under different clamping conditions. The predicted deformation pattern and magnitude were validated by laser scanning data of physical assemblies. In addition, the code was used to investigate residual stresses developed during multipass arc welding of a nuclear industry pressurizer surge nozzle and subsequent welding repair where a 3D simulation was necessary. Taking the experimental data as reference, the 3D model predicted better residual stress distribution than a typical 2D asymmetrical model. Stress evolution in welding repair was also presented and discussed in this study. Furthermore, the efficient numerical model made it feasible to use integrated computational welding engineering to simulate welding processes for large-scale structures.

97 MATHEMATICS AND COMPUTING↗

Pursuing Heteroleptic Ligand Design Principles for Photoactive Fe Complexes with Ultrafast X-ray Emission and Variable-Temperature Optical Spectroscopies

Understanding the key parameters that govern the photophysical and photochemical properties of transition metal complexes is essential for the development of efficient photosensitizers for photocatalytic applications. Achieving this objective necessitates clear and detailed investigations of their electronic excited states, for which time-resolved metal Kβ X-ray emission spectroscopy (XES) has proven highly effective. Here, we present a time-resolved Fe Kβ XES study of a heteroleptic Fe(II) polypyridyl carbene complex, [Fe(phen) 2 (C 4 H 10 N 4 )] 2+ (1; phen = 1,10-phenanthroline), utilizing both the valence-to-core and Kβ mainline spectral regions, complemented by variable-temperature transient optical absorption (VT-TA) spectroscopy. Detailed analysis of the time-resolved Kβ XES data, supported by density functional theory (DFT) calculations and an Eyring analysis of the VT-TA data, reveals parallel excited state relaxation dynamics that support an assignment of the long-lived excited state to a triplet metal-centered state. Placing these results in the context of prior studies of heteroleptic Fe(II) polypyridyl cyanide complexes motivated a series of DFT calculations to investigate the effects of ligand structural flexibility and arrangement. These calculations reinforce the experimentally derived conclusion that constraining structural flexibility with multidentate ligands significantly impacts the excited state relaxation dynamics. Furthermore, our study emphasizes that the arrangement of strong field ligands in heteroleptic complexes substantially affects the energy of Jahn–Teller active triplet metal-centered states in low-spin d 6 metal complexes. Together, these findings provide synthetic design principles for extending metal-to-ligand charge transfer excited state lifetimes of heteroleptic Fe complexes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Crystal Structure of Borrelia turicatae protein, BTA121, a differentially regulated gene in the tick-mammalian transmission cycle of relapsing fever spirochetes

Tick-borne relapsing fever (RF) borreliosis is a neglected disease that is often misdiagnosed. RF species circulating in the United States include Borrelia turicatae, which is transmitted by argasid ticks. Environmental adaptation by RF Borrelia is poorly understood, however our previous studies indicated differential regulation of B. turicatae genes localized on the 150 kb linear megaplasmid during the tick-mammalian transmission cycle, including bta121. This gene is up-regulated by B. turicatae in the tick versus the mammal, and the encoded protein (BTA121) is predicted to be surface localized. The structure of BTA121 was solved by single-wavelength anomalous dispersion (SAD) using selenomethionine-derivative protein. The topology of BTA121 is unique with four helical domains organized into two helical bundles. Due to the sequence similarity of several genes on the megaplasmid, BTA121 can serve as a model for their tertiary structures. BTA121 has large interconnected tunnels and cavities that can accommodate ligands, notably long parallel helices, which have a large hydrophobic central pocket. Preliminary in-vitro studies suggest that BTA121 binds lipids, notably palmitate with a similar order of binding affinity as tablysin-15, a known palmitate-binding protein. The reported data will guide mechanistic studies to determine the role of BTA121 in the tick-mammalian transmission cycle of B. turicatae.

59 BASIC BIOLOGICAL SCIENCES↗

mesoflow [SWR-22-56]

Mesoflow is a continuum scale simulation tool developed specifically for modeling transport and chemistry at the mesoscale. Our solver utilizes Cartesian block-structured adaptive mesh refinement to resolve complex surface morphologies (of catalysts/biomass particles among others) directly obtained from X-ray tomography data. An immersed boundary based formulation enables rapid representation of complex geometries prevalent in most mesoporous interfaces. The solver is developed on top of open-source performance portable library, AMReX, providing parallel execution capabilities on current and upcoming high-performance-computing (HPC) architectures. Our flexible software framework enables integration of complex chemical mechanisms at heterogenous interfaces and time-split algorithms for circumventing highly disparate reaction and flow time-scales. Our current studies indicate a ten-fold performance gain by using graphics-processing-units (GPU) compared to a single processor for representative problem sizes (2 million cell mesh).

Sitaraman, Hariswaran↗

ARM-IRL: Adaptive Resilience Metric Quantification Using Inverse Reinforcement Learning

The resilience of safety-critical systems is gaining importance due to the rise in cyber and physical threats, especially within critical infrastructure. Traditional static resilience metrics may not capture dynamic system states, leading to inaccurate assessments and ineffective responses to cyber threats. This work aims to develop a data-driven, adaptive method for resilience metric learning. We propose a data-driven approach using inverse reinforcement learning (IRL) to learn a single, adaptive resilience metric. The method infers a reward function from expert control actions. Unlike previous approaches using static weights or fuzzy logic, this work applies adversarial inverse reinforcement learning (AIRL), training a generator and discriminator in parallel to learn the reward structure and derive an optimal policy. The proposed approach is evaluated on multiple scenarios: optimal communication network rerouting, power distribution network reconfiguration, and cyber–physical restoration of critical loads using the IEEE 123-bus system. The adaptive, learned resilience metric enables faster critical load restoration in comparison to conventional RL approaches.

97 MATHEMATICS AND COMPUTING↗

Eu(III) and Cm(III) Complexation by the Aminocarboxylates NTA, EDTA, and EGTA Studied with NMR, TRLFS, and ITC—An Improved Approach to More Robust Thermodynamics

The complex formation of Eu(III) and Cm(III) was studied via tetradentate, hexadentate, and octadentate coordinating ligands of the aminopolycarboxylate family, viz., nitrilotriacetate (NTA 3- ), ethylenediaminetetraacetate (EDTA 4- ), and ethylene glycol-bis(2-aminoethyl ether)-N,N,N',N'-tetraacetate (EGTA 4- ), respectively. Based on the complexones’ pK a values obtained from 1 H nuclear magnetic resonance (NMR) spectroscopic pH titration, complex formation constants were determined by means of the parallel-factor-analysis-assisted evaluation of Eu(III) and Cm(III) time-resolved laser-induced fluorescence spectroscopy (TRLFS). This was complemented by isothermal titration calorimetry (ITC), providing the enthalpy and entropy of the complex formation. This allowed us to obtain genuine species along with their molecular structures and corresponding reliable thermodynamic data. The three investigated complexones formed 1:1 complexes with both Eu(III) and Cm(III). Besides the established Eu(III)–NTA 1:1 and 1:2 complexes, we observed, for the first time, the existence of a Eu(III)–NTA 2:2 complex of millimolar metal and ligand concentrations. Demonstrated for thermodynamic studies on Eu(III) and Cm(III) interaction with complexones, the utilized approach is commonly applicable to many other metal–ligand systems, even to high-affinity ligands.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The Flare Genesis Experiment

Using the Flare Genesis Experiment (FGE), a balloon-borne observatory with an 80-cm solar telescope we observed the active region NOAA 8844 on January 25, 2000 for several hours. FGE was equipped with a vector polarimeter and a tunable Fabry-Perot narrow-band filter. It recorded time series of filtergrams, vector magnetograms, and Dopplergrams at the Ca(I) 6122.2 angstrom line, and H-alpha filtergrams with a cadence between 2.5 and 7.5 minutes. At the time of the observations, NOAA 8844 was located at approximately 5 N 30 W. The region was rapidly growing during the observations; new magnetic flux was constantly emerging in three supergranules near its center. We describe in detail how the FGE data were analyzed and report on the structure and behavior of peculiar moving dipolar features (MDFs) observed in the active region. In longitudinal magnetograms, the MDFs appeared to be small dipoles in the emerging fields. The east-west orientation of their polarities was opposite that of the sunspots. The dipoles were oriented parallel to their direction of motion, which was in most cases towards the sunspots. Previously, dipolar moving magnetic features have only been observed flowing out from sunspots. Vector magnetograms show that the magnetic field of each MDF negative part was less inclined to the local horizontal than the ones of the positive part. We identify the MDFs as undulations, or stitches, where the emerging flux ropes are still tied to the photosphere. We present a U-loop model that can account for their unusual structure and behavior, and it shows how emerging flux can shed its entrained mass.

Rust, D. M.↗

Internal structure of the San Jacinto fault zone at the Ramona Reservation, north of Anza, California, from dense array seismic data

SUMMARY We image the internal structure of the San Jacinto fault zone (SJFZ) near Anza, California, with seismic data recorded by two dense arrays (RA and RR) from ∼42 000 local and ∼180 teleseismic events occurring between 2012 and 2017. The RA linear array has short aperture (∼470 m long with 12 strong motion sensors) and recorded for the entire analysed time window, whereas the RR is a large three-component nodal array (97 geophones across a ∼2.4 km × 1.4 km area) that operated for about a month in September–October 2016. The SJFZ at the site contains three near-parallel surface traces F1, F2 and F3 from SW to NE that have accommodated several Mw > 6 earthquakes in the past 15 000 yr. Waveform changes in the fault normal direction indicate structural discontinuities that are consistent with the three fault surface traces. Relative slowness from local events and delay time analysis of teleseismic arrivals in the fault normal direction suggest a slower SW side than the NE with a core damage zone between F1 and F2. This core damage zone causes ∼0.05 s delay at stations RR26–31 in the teleseismic P arrivals compared with the SW-most station, and generates both P- and S-type fault zone trapped waves. Inversion of S trapped waves indicates the core damaged structure is ∼100 m wide, ∼4 km deep with a Q value of ∼20 and 40 per cent S-wave velocity reduction compared with bounding rocks. Fault zone head waves observed at stations SW of F3 indicate a local bimaterial interface that separates the locally faster NE block from the broad damage zone in the SW at shallow depth and merges with a deep interface that separates the regionally faster NE block from rocks to the SW with slower velocities at greater depth. The multiscale structural components observed at the site are related to the geological and earthquake rupture history at the site, and provide important information on the preferred NW propagation of earthquake ruptures on the San Jacinto fault.

Geochemistry & Geophysics↗

Overview of NASA GRC Stirling Technology Development

The Stirling Radioisotope Generator (SRG) is currently being developed by Lockheed Martin Astronautics (LMA) under contract to the Depar1ment of Energy (DOE). The generator will be a high efficiency electric power source for NASA Space Science missions with the capability to operate in the vacuum of deep space or in an atmosphere such as on the surface of Mars. High system efficiency is obtained through the use of free-piston Stirling power conversion technology. Power output of the generator will be greater than 100 watts at the beginning of life with the decline in power being largely due to the decay of the plutonium heat source. In suppOl1 of the DOE SRG project, the NASA Glenn Research Center (GRC) has established a near-term technology effort to provide some of the critical data to ensure a successful transition to flight for what will be the first dynamic power system used in space. Initially, a limited number of technical areas were selected for the GRC effort, however this is now being expanded to more thoroughly cover a range of technical issues. The tasks include in-house testing of Stirling convertors and controllers, materials evaluation and heater head life assessment, structural dynamics, electromagnetic interference, organics evaluation, and reliability analysis. Most of these high-level tasks have several subtasks within. There is also an advanced technology effort that is complementary near-term technology effort. Many of the tests make use of the 55-We Technology Demonstration Convel10r (TDC). There have been multiple controller tests to support the LMA flight controller design effort. Preparation is continuing for a thermal/vacuum system demonstration. A pair of flight prototype TDC's have recently been placed on an extended test with unattended, continuous operation. Heater head life assessment efforts continue, with the material data being refined and the analysis moving toward the system perspective. Long-term magnet aging tests are continuing to characterize any possible aging in the strength or demagnetization resistance of the permanent magnets used in the linear alternator. In a parallel effort, higher performance magnets are also being evaluated. A reliability effort is being initiated that will help to guide the development activities with an increased focus on the necessary components and subsystems. Some other disciplines that are active in the GRC technology effort include structural dynamics, linear alternator analysis, EMI/EMC, controls, and mechanical design evaluation. This paper will provide an overview of some of the GRC technical efforts, including the current status, and a description of future efforts.

Schreiber, Jeffrey↗

A Tower-based Prototype VHF/UHF Radar for Subsurface Sensing: System Description and Data Inversion Results

Knowledge of subsurface characteristics such as permittivity variations and layering structure could provide a breakthrough in many terrestrial and planetary science disciplines. For Earth science, knowledge of subsurface and subcanopy soil moisture layers can enable the estimation of vertical flow in the soil column linking surface hydrologic processes with that in the subsurface. For planetary science, determining the existence of subsurface water and ice is regarded as one of the most critical information needs for the study of the origins of the solar system. The subsurface in general can be described as several near-parallel layers with rough interfaces. Each homogenous rough layer can be defined by its average thickness, permittivity, and rms interface roughness assuming a known surface spectral distribution. As the number and depth of layers increase, the number of measurements needed to invert for the layer unknowns also increases, and deeper penetration capability would be required. To nondestructively calculate the characteristics of the rough layers, a multifrequency polarimetric radar backscattering approach can be used. One such system is that we have developed for data prototyping of the Microwave Observatory of Subcanopy and Subsurface (MOSS) mission concept. A tower-mounted radar makes backscattering measurements at VHF, UHF, and L-band frequencies. The radar is a pulsed CW system, which uses the same wideband antenna to transmit and receive the signals at all three frequencies. To focus the beam at various incidence angles within the beamwidth of the antenna, the tower is moved vertically and measurements made at each position. The signals are coherently summed to achieve focusing and image formation in the subsurface. This requires an estimate of wave velocity profiles. To solve the inverse scattering problem for subsurface velocity profile simultaneously with radar focusing, we use an iterative technique based on a forward numerical solution of the layered rough surface problem. The layers are each defined in terms of a small number of unknown distributions as given above. An a priori estimate of the solution is first assumed, based on which the forward problem is solved for the backscattered measurements. This is compared with the measured data and using iterative techniques an update to the solution for the unknowns is calculated. The process continues until convergence is achieved. Numerical results will be shown using actual radar data acquired with the MOSS tower radar system in Arizona in Fall 2003, and compared with in-situ measurements.

Moghaddam, Mahta↗

Structure of current sheets in magnetic holes at 1 AU

Current density profiles in several types of interplanetary magnetic holes have been calculated using high-resolution Imp 6 magnetic field data (12.5 vector measurements/s), assuming that the currents flow in planar sheets and that the magnetic field varies only in the direction normal to the sheet. The planarity was verified in four holes which were observed by two suitably spaced spacecraft. Four types of simple magnetic holes are discussed, in which B varies nearly monotonically on each side of the hole. In two of the holes, B varies in intensity but not in direction as a result of currents normal to B. In the other two holes, B changes in both magnitude and direction as a result of currents both normal and parallel to B. The observed structures are found to be qualitatively consistent with the models of Burlaga and Lemaire (1978). Examples of complex irregular magnetic holes are also presented, and they are shown to contain multiple current sheets in which currents flow parallel to one another at various angles with respect to B. There is no model of such magnetic holes at present.

Fitzenreiter, R. J.↗

Multifrontal Non-negative Matrix Factorization

Non-negative matrix factorization (Nmf) is an important tool in high-performance large scale data analytics with applications ranging from community detection, recommender system, feature detection and linear and non-linear unmixing. While traditional Nmf works well when the data set is relatively dense, however, it may not extract sufficient structure when the data is extremely sparse. Specifically, traditional Nmf fails to exploit the structured sparsity of the large and sparse data sets resulting in dense factors. We propose a new algorithm for performing Nmf on sparse data that we call multifrontal Nmf (Mf-Nmf) since it borrows several ideas from the multifrontal method for unconstrained factorization (e.g. LU and QR). We also present an efficient shared memory parallel implementation of Mf-Nmf and discuss its performance and scalability. We conduct several experiments on synthetic and realworld datasets and demonstrate the usefulness of the algorithm by comparing it against standard baselines. We obtain a speedup of 1.2x to 19.5x on 24 cores with an average speed up of 10.3x across all the real world datasets.

Sao, Piyush↗