Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “MATRICES”

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 1,027 records · Page 57

Restoring the Dunes: Using Satellite Imagery to Help Restore Turtle Nesting Grounds In Puerto Rico After Hurricane Maria

The North coast of Puerto Rico sustained heavy damage from Hurricane Maria and other subsequent tropical storms. This damage included erosion of beaches and foredunes. The foredunes are a delicate ecological environment, which provides important nesting grounds for endangered turtles. Local ecologists have been working to restore these dunes with mixed results. Presented here is an evaluation of dune restoration sites, and strategies to improve the success of subsequent restoration efforts. This involves both historical satellite image analysis and field work to evaluate sediment migration pathways. Local ecologists use a technique called Biomimicry to restore dunes. This involves inserting sticks or plants into the ground, simulating the effects of natural vegetation to slow wind speeds and deposit wind-blown sand. Manually planted vegetation can follow to stabilize the trapped sand. While often successful, this technique occasionally fails to accumulate significant amounts of sand. Satellite imagery shows that one of the controlling factors is the orientation of the dune and beach to the prevailing wind. Dunes that are oriented parallel to the prevailing wind are relatively slower-growing due to sediment bypass. Satellite imagery also indicates that the fastest growing dunes are in close proximity to major sediment sources. This is complicated by large dams, which limit the amount of sediment delivered to the coast. Field observation shows that biomimicry matrices are occasionally located too close to the swash zone, where storms over-top the dunes, poisoning stabilizing vegetation with salt water. Another common observation was that many dunes were starting to form behind fallen trees, suggesting that a different type of biomimicry might be successful. Satellite and field observation suggest that strategies to protect nearshore environments and improve dune restoration biomimicry efforts include: 1) orient biomimicry matrices at a higher angle to prevailing winds. 2) Prioritize sites near river mouths, particularly those without upstream dams. 3) Test other orientations of biomimicry, such as horizontal boards to mimic fallen trees. 4) Improve sediment supply by preventing further dam development, and releasing trapped sediment.

Paul Bremner↗

A Model-Based Framework for NASA Science Mission Formulation

This paper details an effort to implement NASA systems engineering standards and practices using a Digital System Model (DSM), which we also refer to as the system architecture model (SAM). Creating a SAM at the beginning of a design effort helps systems engineer identify errors, inconsistencies, and miscommunications as early as possible. Such issues can then be resolved before it becomes costly and requires significant rework to do so. A pre-formulation SAM also enables advanced trade studies at the earliest stage of concept development. This results in “win-wins” where architecture changes can reduce cost without decreasing effectiveness, or increase effectiveness without increasing cost. SAMs also streamline the creation of project documentation and facilitate superior dialogue between science, engineering, and management stakeholders. Common artifacts such as Master Equipment Lists (MELs), Science Traceability Matrices (STMs), and Mission Traceability Matrices (MTMs), can be automatically maintained through requirements and structural models in the SAM. Key performance parameters (and changes to them) can be tied to simulations in the SAM, allowing far more rapid and extensive trade space exploration. As will be shown in this work, these advantages have been realized for the first time in an actual NASA Goddard Space Flight Center (GSFC) pre-phase A and phase A concept study. Modelbased design tools have been developed in a general, modular manner to maximize reuse on future programs. Focus is placed on structural and requirements architecture modeling. Data is ingested into the SAM from existing discipline model outputs (e.g. MS Excel spreadsheets) and used to update a SysML structural model. A mission requirements model is created within the SAM, containing mission specific as well as standard GSFC mission requirements. With the linked requirements model and structural model, automated compliance checking will be performed as mission parameters evolve without the need for discipline engineers to work directly with the SAM in SysML.

MBSE↗

Effect of salinity on the preservation of biomarkers in hypersaline microbial mat kerogens

Kerogen is the residue that remains after the minerals and organic compounds have been extracted from a rock using acid and organic solvents, respectively [1]. It is not only the largest organic carbon sink in the geological record but the best-preserved fraction due to its recalcitrance and ability to retain syngenetic source information during deposition [2]. Due to the complex, macromolecular, cross-linked polymeric structure of kerogen, it is difficult to analyse intact, often requiring breakdown into smaller subunits through thermal (pyrolysis) or chemical treatment (chemolysis) [3]. Saline and hypersaline environments host diverse microbial communities, which can trap and bind organic matrices such as cellular components (e.g., extracellular polymeric substances, cell wall, peptidoglycan, sheath material, etc.) with sediments, forming layered, laminated structures called microbial mats. Here we present depth profiles of two hypersaline cores through a microbial mat with differing total salinities. We looked at the composition of organic compounds released from kerogen during high temperature, hydrogen-assisted, pyrolysis to understand the effect of salinity on the composition and timing of their incorporation [4]. Field Site: The Exportadora del Sal S. A. (ESSA) saltern of Guerrero Negro (Baja California Sur, Mexico) lies within a typical sabhka environment adjacent to Laguna Ojo de Liebre. Note the spatial extent of the 13 major ponds of the ESSA system. The GN ponds thrive under well-monitored and stable conditions with steady state accretion and degradation of mat layers (around 0.5–1.4 cm per year). Organic mats typically attain up to 10 cm in thickness, depending on the location, which represents about 60 years of growth since these ponds were established. Sedimentation rate is rapid in these non-lithifying and non-mineralising mats, 1–1.5 cm per year. This environment is ideal for investigating in situ lipid preservation and diagenesis without the influence of constantly changing physical parameters. Methods: P4n5 and P5AB cores were collected as ~9 cm and 6.5 cm fresh cores, respectively, in June 2001 and September 2010, respectively from the ESSA saltern. The P4n5 core (9-9.1% salinity) was taken from Pond 4 near 5 and sub-divided into 10 layers. The P5AB core (11.2% salinity) was taken from Pond 5A near 5B and sub-divided into 8 layers. Lyophilised microbial mat powders were solvent-extracted using a modified Bligh-Dyer method and analysed by gas chromatography-mass spectrometry (GC-MS). Additionally, a subset of layers—layer 4 and layer 7—which represented active biomass and sediment-processed, respectively, were subjected to mild-acid methanolysis. The respective residues from conventional solvent extraction and acid methanolysis were loaded with a molybdenum sulfide catalyst and placed into a stainless-steel reactor. Samples were run on the hydropyrolysis set up (heated to 500 °C with 13-15 MPa of H2 pressure), extracted, and analysed by GC-MS. Results: Molecular profiles from the two ponds were significantly different and reflect the contributions of different photosynthetic and respiratory microbial communities to preserved organic matter. For example, the higher salinity P5AB core showed evidence of greater archaeal inputs (similarly observed in lab culture experiments [5]) as well as potential kerogen-bound carotenoids/carotenoid rearranged products or fragments. A hydrophobic emulsion formed during the acid methanolysis processing of layer 4 from both P4n5 and P5AB and layer 7 of P5AB. Hydropyrolysis of this hydrophobic residues released a greater abundance of polycyclic lipids relative to the pre-extracted control. Implications and Future Work: Catalytic hydropyrolysis of kerogen can rapidly generate abundant saturated pyrolysate products from the bound biomarker pool without altering the structures or stereochemistries of the products [6]. Chemical processing of solvent-extractable residues indicated that cellular matrices such as extracellular polymeric substances (EPS) may play a role in the sequestration of polycyclic lipid biomarkers such as steranes and hopanes. Their higher relative abundance compared to the pre-extracted control in the higher salinity layers and cores preliminarily aligns with this hypothesis, although further experiments are required to confirm this. It has been well-documented that EPS plays an important role in mineralisation, specifically carbonate formation, in microbial ecosystems. EPS can enhance calcium carbonate precipitation by providing diffusion-limited sites that create alkalinity gradients in response to microbial processes [7]. It has been experimentally demonstrated that salinity influences the total amount of EPS (both loosely- and tightly-bound) which increases with increasing salinity [8]. The higher salinity at Guerrero Negro switches the microbial population from filamentous Microcoleus to being dominated by Phormidium, Oscillatoria, and unicellular cyanobacteria. We demonstrated that lipid binding into kerogen via strong covalent linkages occurs at the very earliest stages of sedimentary diagenesis. We are currently investigating the influence of salinity on organic preservation through experimental and modelling approaches. Understanding how key biomarkers transform into preserved organic matter (i.e., from precursor biolipids to bound geolipids) in brine ecosystems will aid the search for organic biosignatures on other planetary bodies, especially Icy Moons and modern Mars.

hypersaline↗

Evaluation of the Interference of Tenax®TA Adsorbent With Dimethylformamide Dimethyl Acetal Reagent for Gas Chromatography-Dragonfly Mass Spectrometry and Future Gas Chromatography-Mass Spectrometry in Situ Analysis

Among future space missions, national aeronautics and space administration (NASA) selected two of them to analyze the diversity in organic content within Martian and Titan soil samples using a gas chromatograph – mass spectrometer (GC–MS) instrument. The Dragonfly space mission is planned to be launched in 2027 to Titan's surface and explore the Shangri-La surface region for years. One of the main goals of this mission is to understand the past and actual abundant prebiotic chemistry on Titan, which is not well characterized yet. The ExoMars space mission is planned to be launched in 2028 to Mars’ surface and explore the Oxia Planum and Mawrth Vallis region for years. The main objectives focus on the exploration of the subsurface soil samples, potentially richer in organics, that might be relevant for the search of past life traces on Mars where irradiation does not impact the matrices and organics. One recently used sample pre-treatment for gas chromatography – mass spectrometry analysis is planned on both space missions to detect refractory organic molecules of interest for astrobiology. This pre-treatment is called derivatization and uses a chemical reagent – called dimethylformamide dimethyl acetal (DMF-DMA) – to sublimate organic compounds keeping them safe from thermal degradation and conserving the chirality of the molecules extracted from Titan or Mars’ matrices. Indeed, the detection of building blocks of life or enantiomeric excess of some organics (e.g. amino acids) after DMF-DMA pre-treatment and GC–MS analyses would be both bioindicators. The main results highlighted by our work on DMF-DMA and Tenax®TA interaction and efficiency to detect organic compounds at ppb levels in a fast and single preparation are first that Tenax®TA did not show the onset of degradation until after 150 experiments – a 120 h at 300 °C experiment – which greatly exceeds the experimental lifetimes for the DraMS and GC-space in situ investigations. Tenax®TA polymer and DMF-DMA produce many by-products (about 70 and 46, respectively, depending on the activation temperature). Further, the interaction between the two leads to the production of 22 additional by-products from DMF-DMA degradation, but these listed by-products do not prevent the detection of trace-level organic molecules after their efficient derivatization and volatilization by DMF-DMA in the oven ahead the GC–MS trap and column.

DraMS-Dragonfly mission↗

Approximate Inverse Chain Preconditioner: Iteration Count Case Study for Spectral Support Solvers

As the growing availability of computational power slows, there has been an increasing reliance on algorithmic advances. However, faster algorithms alone will not necessarily bridge the gap in allowing computational scientists to study problems at the edge of scientific discovery in the next several decades. Often, it is necessary to simplify or precondition solvers to accelerate the study of large systems of linear equations commonly seen in a number of scientific fields. Preconditioning a problem to increase efficiency is often seen as the best approach; yet, preconditioners which are fast, smart, and efficient do not always exist. Following the progress of [1], we present a new preconditioner for symmetric diagonally dominant (SDD) systems of linear equations. These systems are common in certain PDEs, network science, and supervised learning among others. Based on spectral support graph theory, this new preconditioner builds off of the work of [2], computing and applying a V-cycle chain of approximate inverse matrices. This preconditioner approach is both algebraic in nature as well as hierarchically-constrained depending on the condition number of the system to be solved. Due to its generation of an Approximate Inverse Chain of matrices, we refer to this as the AIC preconditioner. We further accelerate the AIC preconditioner by utilizing precomputations to simplify setup and multiplications in the con-text of an iterative Krylov-subspace solver. While these iterative solvers can greatly reduce solution time, the number of iterations can grow large quickly in the absence of good preconditioners. Initial results for the AIC preconditioner have shown a very large reduction in iteration counts for SDD systems as compared to standard preconditioners such as Incomplete Cholesky (ICC) and Multigrid (MG). We further show significant reduction in iteration counts against the more advanced Combinatorial Multigrid (CMG) preconditioner. We have further developed no-fill sparsification techniques to ensure that the computational cost of applying the AIC preconditioner does not grow prohibitively large as the depth of the V-cycle grows for systems with larger condition numbers. Our numerical results have shown that these sparsifiers maintain the sparsity structure of our system while also displaying significant reductions in iteration counts.1 2

97 MATHEMATICS AND COMPUTING↗

Study on Fracture of Fiber-Reinforced Polymeric Composites Using Spiral Notch Torsion Test

Fiber-reinforced polymer-matrix composites have been widely applied in various areas due to their excellent mechanical properties. In this study, carbon fiber-reinforced polymer-matrix composites and glass fiber-reinforced polymer-matrix composites failure tests under quasi static torsion have been demonstrated by using spiral notch torsion test methodology. Three types of cylindrical specimens with different fiber orientations were prepared and tested. For both composites, after increase stage, A type exhibited limited stress variations before the sudden failure, whereas B type showed gradual failure with torque drops. Displacement contours, stress contours and energy release rates were acquired using the anisotropic elastic finite element models. Microscopic characterization exhibited that the failure of A type occurred along interfaces between fibers and matrices, whereas circular cracks in B type propagated alternatively between polymeric matrices and interfaces. Type C failure is similar to type A whereas with a curvature observed on the fracture surface.

Du, Fen↗

Approximate Inverse Chain Preconditioner: Iteration Count Case Study for Spectral Support Solvers

As the growing availability of computational power slows, there has been an increasing reliance on algorithmic advances. However, faster algorithms alone will not necessarily bridge the gap in allowing computational scientists to study problems at the edge of scientific discovery in the next several decades. Often, it is necessary to simplify or precondition solvers to accelerate the study of large systems of linear equations commonly seen in a number of scientific fields. Preconditioning a problem to increase efficiency is often seen as the best approach; yet, preconditioners which are fast, smart, and efficient do not always exist. Following the progress of [1], we present a new preconditioner for symmetric diagonally dominant (SDD) systems of linear equations. These systems are common in certain PDEs, network science, and supervised learning among others. Based on spectral support graph theory, this new preconditioner builds off of the work of [2], computing and applying a V-cycle chain of approximate inverse matrices. This preconditioner approach is both algebraic in nature as well as hierarchically-constrained depending on the condition number of the system to be solved. Due to its generation of an Approximate Inverse Chain of matrices, we refer to this as the AIC preconditioner. We further accelerate the AIC preconditioner by utilizing precomputations to simplify setup and multiplications in the con-text of an iterative Krylov-subspace solver. While these iterative solvers can greatly reduce solution time, the number of iterations can grow large quickly in the absence of good preconditioners. Initial results for the AIC preconditioner have shown a very large reduction in iteration counts for SDD systems as compared to standard preconditioners such as Incomplete Cholesky (ICC) and Multigrid (MG). We further show significant reduction in iteration counts against the more advanced Combinatorial Multigrid (CMG) preconditioner. We have further developed no-fill sparsification techniques to ensure that the computational cost of applying the AIC preconditioner does not grow prohibitively large as the depth of the V-cycle grows for systems with larger condition numbers. Our numerical results have shown that these sparsifiers maintain the sparsity structure of our system while also displaying significant reductions in iteration counts.1 2

97 MATHEMATICS AND COMPUTING↗

Nonlinear Matrix Approximation with Radial Basis Function Components

We introduce and investigate matrix approximation by decomposition into a sum of radial basis function (RBF) components. An RBF component is a generalization of the outer product between a pair of vectors, where an RBF function replaces the scalar multiplication between individual vector elements. Even though the RBF functions are positive definite, the summation across components is not restricted to convex combinations and allows us to compute the decomposition for any real matrix that is not necessarily symmetric or positive definite. We formulate the problem of seeking such a decomposition as an optimization problem with a nonlinear and non-convex loss function. Several modern versions of the gradient descent method, including their scalable stochastic counterparts, are used to solve this problem. We provide extensive empirical evidence of the effectiveness of the RBF decomposition and that of the gradient-based fitting algorithm. While being conceptually motivated by singular value decomposition (SVD), our proposed nonlinear counterpart outperforms SVD by drastically reducing the memory required to approximate a data matrix with the same L2 error for a wide range of matrix types. For example, it leads to 2 to 6 times memory save for Gaussian noise, graph adjacency matrices, and kernel matrices. Moreover, this proximity-based decomposition can offer additional interpretability in applications that involve, e.g., capturing the inner low-dimensional structure of the data, retaining graph connectivity structure, and preserving the acutance of images.

Rebrova, Elizaveta↗

Structure Perception in 3D Point Clouds: Preprint

Understanding human perception is critical to the design of ef- fective visualizations. The relative benefits of using 2D versus 3D techniques for data visualization is a complex decision space, with varying levels of uncertainty and disagreement in both the liter- ature and in practice. This study aims to add easily reproducible, empirical evidence on the role of depth cues in perceiving structures or patterns in 3D point clouds. We describe a method to synthesize a 3D point cloud that contains a 3D structure, where 2D projec- tions of the data strongly resemble a Gaussian distribution. We performed a within-subjects structure identification study with 128 participants that compared scatterplot matrices (canonical 2D projections) and 3D scatterplots under three types of motion: rota- tion, xy-translation, and z-translation. We found that users could consistently identify three separate hidden structures under ro- tation, while those structures remained hidden in the scatterplot matrices and under translation. This work contributes a set of 3D point clouds that provide definitive examples of 3D patterns per- ceptible in 3D scatterplots under rotation but imperceptible in 2D scatterplots.

data analysis↗

Fast simultaneous feasibility testing and applications for electrical power grids

Fast simultaneous feasibility testing (SFT) for management of an electrical power grid is achieved through various innovations. The computation problem relates to evaluation of candidate solutions for external power flows into a power grid, with respect to predetermined constraints and contingencies. Storage and computations are reduced by formulating the problem in terms of transactional nodes (e.g. third party connections for generators and loads) instead of the larger number of bus nodes. Further advantages are achieved by precomputing matrices that can be reused across multiple SFT invocations, organizing matrices and operations to reduce storage and computation, and eliminating branches that have no contingency violations. The fast SFT further enables new applications for integrating SFT with a Security-Constrained Unit Commitment (SCUC) optimizer for efficient single pass solutions; operating SFT as a service to multiple SCUC solvers; or integration with Security-Constrained Economic Dispatch (SCED) optimization.

97 MATHEMATICS AND COMPUTING↗

Synthesis and Characterization of Polymeric Pore Former Microspheres using Modified Polymer Resin

Porosity in polymer parts is often desired to improve mechanical properties. Individual pores must be small enough that bulk mechanical properties are retained but significant enough to reduce the density of the material. Filler materials can be incorporated into polymer matrices and subsequently removed to create voids. The materials used to create the voids are known as pore-formers. The present research aims to create spherical pores in a silicone matrix by inserting and removing core materials or fillers post-cure. A method was developed to incorporate pore formers into an uncured polymer resin and remove the materials from the matrix. Polyethylene Glycol (PEG) and Polylactic Acid (PLA) microspheres are created and incorporated into silicone matrices as pore formers. Microspheres with controllable size distributions are formed via ultrasonic spray pyrolysis and emulsion formation. Microspheres are removed from the matrix by solvent extraction and calcination. PLA microspheres were modified by hydroxyl addition to the surface to study their solubility in water. PEG-PLA block copolymers are also synthesized and investigated, comparing their physical and chemical properties to the monomers to identify whether the block copolymer microspheres will be more easily removable than their monomeric microspheres. Moreover, other properties including thermal stability, rheology, and mechanical response are studied. The present dissertation provides a comprehensive understanding of the synthesis and characterization of PEG and PLA-based pore formers, toward the future development of feedstock for additive manufacturing porous structures.

Damptey, Ransford↗

A Laboratory Assessment and Optimization of Nanoparticle-Stabilized CO2 Foams for Hydraulic Fracturing Applications

In this research, we have conducted a comprehensive assessment and optimization of aqueous CO2 foams for fracturing applications. The performance of CO2 foams is evaluated using three major matrices: foam stability, foam apparent viscosity, and leakoff. The response of these performance matrices to different internal factors and external conditions is evaluated. These factors and conditions include foam quality, nanoparticle concentration, salinity, shear rate, surfactant type and concentration, polymer concentration, and temperature. Before foam generation, the colloidal stability of the nanosilica dispersion is evaluated at different types of electrolytes and under high temperatures. With a good salinity resistance and thermal stability, the nanosilica is then used for the entire set of foam generation experiments. The foam configuration with optimized parameters is used for the leakoff test, which serves as a supplementary section to further demonstrate the advantages of the proposed foam system. Then the final commercialization section provides a general view of what should be considered to commercialize the research results. In general, it is found that the aqueous CO2 foam system has demonstrated a good potential to be used as a fracturing fluid. The design of the foam system needs to consider formations properties and formation fluid properties etc. to realize the optimal performance of the system on a laboratory scale.

02 PETROLEUM↗

Potential Advantages to Employing an Onsite Laboratory at FUSRAP Sites Contaminated with Radionuclides and Metals - 20353

Various factors that support an evaluation of the relative costs, advantages and disadvantages of utilizing an onsite laboratory or an offsite laboratory for remediation of Formerly Utilized Sites Remedial Action Program (FUSRAP) sites contaminated with radionuclides and metals are presented. In terms of how much more quickly an onsite laboratory provides results versus an offsite laboratory, the results of this study provide tools for calculating the amount of cost savings per day for excavation of contaminated soils for various scenarios including: - State, municipal or local property remediation where the remediation schedule becomes time-critical due to the extra cost and inconvenience to the town or municipality caused by the need for increased security and safety personnel, possible road closures and detours, and the inconvenience to the public due to modified or decreased services of State or municipal facilities; - Residential property where the remediation schedule becomes time-critical due to the inconvenience to the homeowner as well as the importance of maintaining good community relations; - Commercial property where the remediation schedule becomes time-critical due to the hardship and potential loss of revenue endured by the property owner; - Open area where treatment of excavation water is required, and - Open area where treatment of excavation water is not required. For every scenario, the number of matrices, radionuclide and metal contaminants of concern and the project duration are important factors. The costs for the onsite laboratory included the capital costs associated with setting up the lab such as trailer(s) or building(s), power and water hookup, analytical testing equipment, balances, computers, desks and bookcases; support equipment such as hot plates, filtration assemblies, safety equipment, ovens, etc.; and supplies such as glassware, pipettes, filters, etc. and ongoing costs such as labor, chemicals and reagents, standards, cost of the laboratory obtaining certification and maintenance. Offsite laboratory sample testing costs are generally very competitive and typically substantially lower than onsite laboratory costs if only the sample testing costs are considered. The great disadvantage of using an offsite laboratory is that they often cannot guarantee rapid turnaround times for results which creates project delays associated with the remediation process including delays determining the completeness of the characterization and remediation processes, delays in backfilling a property and additional costs for treating wastewater (excavation water, stormwater or both). These costs will vary greatly depending upon the area of the country in which the project is located due to the wide range of labor costs for the onsite laboratory and field workers, as well as the contaminants of interest, matrices and most importantly the size of the impacted area and the number of properties located in that area. This paper provides guidelines for calculating the cost savings that may be realized if an onsite laboratory can generate results more quickly than an offsite laboratory by estimating the costs associated with the remediation process for the various scenarios noted above such as work crew labor costs, excavator and dump truck equipment rentals, wastewater treatment costs (if applicable) and the additional costs incurred by homeowners, commercial property owners and towns/municipalities that are reimbursed by the potentially responsible parties. Key assumptions involve the typical quickest turnaround time provided by the onsite and offsite laboratories. Costs for the operation of a wastewater treatment system and the remediation of a site area or property were provided by personnel who have been involved in the operation of a large FUSRAP site in New Jersey for many years. This study provides valuable information for determining the most cost-effective approach for analysis of project samples and therefore may be used on upcoming FUSRAP projects. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Nuclear data uncertainty quantification for the nuclide inventory of a Calvert Cliffs spent fuel sample

The impact of nuclear data cross section uncertainties and covariance matrices on the nuclide vector of spent nuclear fuel was investigated; This exercise was carried out for the Calvert Cliffs fuel assembly D047 benchmark available in the SFCOMPO database. Sample P irradiated in rod MKP109 for 4 cycles up to a burnup of approximately 44 GWd/MTU was selected for the analysis. Nuclear data uncertainties were taken from the most recent libraries released by evaluation projects JEFF, ENDF/B and JENDL, and were propagated using the SANDY code via a stochastic sampling approach. This paper provides a quantification of the uncertainty on the concentration of several actinides and fission products relevant for spent fuel management. Uncertainties generally below 5 % were predicted for the concentrations of most of the uranium, neptunium and plutonium isotopes relevant for SNF applications. Curium isotopes carry larger uncertainties that might exceed 10 %. The contribution of cross section uncertainties on the concentrations of fission products was found to be marginal with the exception of a few nuclides. These results can be significantly affected by the lack of evaluated covariance matrices for the capture cross section of several fission products. Burnup tracers such as {sup 148}Nd and {sup 137}Cs have negligible uncertainties because of the power normalisation imposed in every stochastic calculation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

autoGEMM: Pushing the Limits of Irregular Matrix Multiplication on Arm Architectures

This paper presents an open-source library that pushes the limits of performance portability for irregular General Matrix Multiplication (GEMM) on the widely-used Arm architectures. Our library, autoGEMM, is designed to support a wide range of Arm processors: from edge devices to HPC-grade CPUs. autoGEMM generates optimized kernels for various hardware configurations by auto-combining fragments of autogenerated micro-kernels that employ hand-written optimizations to maximize computational efficiency. We optimize the kernel pipeline by tuning the register reuse and the data load/store overlapping. In addition, we use a dynamic tiling scheme to generate balanced tile shapes. Finally, we position autoGEMM on top of the TVM framework where our dynamic tiling scheme prunes the search space for TVM to identify the optimal combination of parameters for code optimization. Evaluations on five different classes of Arm chips demonstrate the advantages of autoGEMM. For small matrices, autoGEMM achieves 98% of peak and up to 2.0x speedup over state-of-the-art libraries such as LIBXSMM and LibShalom. For irregular matrices (i.e. tall skinny and long rectangles), autoGEMM is 1.3-2.0x faster than widely-used libraries such as OpenBLAS and Eigen. autoGEMM is publicly available at: https://github.com/wudu98/autoGEMM.

Wu, Du↗

Online Alpha Monitoring of High Cs-137 Hanford and SRS High Level Waste with Tensioned Metastable Fluid Detectors

The Department of Energy’s Hanford and Savannah River Sites maintain millions of gallons of caustic supernate and salt high activity waste in their high-level waste (HLW) tank farm inventories. The Savannah River Site is currently treating this waste with a calixarene-based solvent extraction of Cs-137 to reduce these inventories. Hanford is employing an at-tank crystalline silico titanate (CST) based solid phase extraction methodology to reduce their liquid HLW inventories. Due to the high solubility of Cs-137 and the relative insolubility of the actinides in these caustic waste forms, the beta to alpha radioactivity ratio can often exceed six orders of magnitude in the feed solutions to these treatment processes. This unique characteristic leads to significant technical challenges in making rapid gross alpha measurements in the presence of the overwhelming beta, gamma, as well as dissolved sodium salt in these HLW matrices. Conventional radioanalytical techniques, such as liquid scintillation analysis or gas flow proportional counting require significant radiochemistry preparation prior to the radiometric measurements for gross alpha activity. These required pretreatments render these technologies untenable for rapid quantification of gross alpha activity that could be required to support on or at-line measurements ensuring a waste stream will meet regulatory requirements. The radiation measurement properties of Tensioned Metastable Fluid Detectors (TMFDs) have been studied by Purdue University’s Taleyarkhan research group for well over a decade. Fluids tensioned to the appropriate degree will rupture when struck by radiation, resulting in a measurable cavitation event. The negative pressure generating this tension can be adjusted by centrifugal rotation or by acoustic means in such a way that these cavitation events can be generated from alpha radiation but will not be generated by beta or gamma radiation. Purdue University and the Savannah River National Laboratory are currently collaborating to develop a gamma/beta blind, spectroscopic alpha measurement system based on the Tensioned Metastable Fluid Detector technology to provide a potential solution for performing rapid gross alpha measurements on these high gamma/beta sample matrices. Measurements using the Indirect Drive Acoustically Tensioned Metastable Fluid Detectors developed as part of this collaboration were performed with an alpha emitting radionuclide. Successful determination of gross alpha activity was observed, indicating a potential pathway for rapid gross alpha measurements in remote-handled shielded cells or in process situations requiring online alpha monitoring. Measurements using this system have been conducted on high beta activity solutions, demonstrating the beta blind capability of this system. Measurements are currently underway to test the system’s capability to measure gross alpha activity on Savannah River Site high level waste high Cs-137 samples that have been measured by the SRNL radiochemistry team. This work was supported by the DOE EM Technology Development program.

DiPrete, David [Savannah River National Laboratory↗

Method used in Tokunaga, T. et al. 2018 for simulating time- and depth- dependent temperature variations in soil and sediment in Rifle, Colorado, 2013-2016.

This data package contains temperature, moisture content, matric potential, and carbon dioxide (CO2) concentration data from the Rifle, Colorado field site and from a 2 meter tall laboratory soil/sediment column designed to simulate field conditions. In addition, redox potential data from the laboratory column are included. The purpose of the data and paper is to demonstrate how a laboratory column was constructed to reproduce the dynamics of nonisothermal conditions in the field, such that CO2 production above a variable-depth water table can be studied in the laboratory. The data package contains 3 data files corresponding to the figures in the related paper (Tokunaga, T. et al. 2018) and 1 README file with the paper information.

54 ENVIRONMENTAL SCIENCES↗

Accelerating structural dynamics simulations with localised phenomena through matrix compression and projection‐based model order reduction

In this work, a novel approach is introduced for accelerating the solution of structural dynamics problems in the presence of localised phenomena, such as cracks. For this category of problems, conventional projection-based Model Order Reduction (MOR) methods are either limited with respect to the range of system configurations that can be represented or require frequent solutions of the Full Order Model (FOM) to update the low-dimensional spaces, in which solutions are represented. In the proposed approach, low-dimensional spaces, constructed for the healthy structure, are enriched with appropriately selected columns of the flexibility matrix of the system. It can be shown that these spaces contain the solution to the original problem for the static case, while their dimension is much smaller. In order to allow their online construction for arbitrary localised features, the full flexibility matrix of the system should be available. To this end, a hierarchical representation is used for the matrices involved, allowing to compute the flexibility matrix efficiently and with reduced memory requirements. The resulting method offers significant speedups, without sacrificing the flexibility and accuracy of the full order model. The performance and limitations of the approach are studied through a series of examples in structural dynamics.

fracture mechanics↗