DETERMINATION OF LIGAMENT QUALITY FACTORS IN ADDITIVELY MANUFACTURED LATTICE STRUCTURES USING IN-SITU COMPRESSION TESTING MICRO-CT
Explore the source record for details and available documents.
DINOVA, VINCENT A.↗
SEARCH · Search NASA
Results for “compressibility factor”
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.
Characterizing The Mechanical Properties of Polyurethane Foams
The Department of Energy has many radioactive facilities that are on the path of deactivation and decommissioning (D and D). These facilities can sit cold and dark for many years before final disposition, and must be maintained to ensure no radiological release occurs in the interim. Improvements and additions to D and D tool sets can greatly: Save time and money, Reduce worker risk. Conventional fixatives widely used often take the form of paints or films that are not readily applicable to 3-dimensional void spaces. Foams are one promising platform that may offer solutions to a number of contaminated problem sets such as: Gloveboxes, Pipes, Tanks. SRNL is working to characterize various commercial foams that would encapsulate the interior volume of a given space and are capable of immobilizing any remaining contamination. One key performance metric of these foams is how well they will perform in an accident scenario. To this end, SRNL is researching the mechanical properties of these foams to ensure that the material can withstand the environment of application while maintaining structural integrity. ASTM standard E3191-18 served as a guiding document for this project, outlining the requirements that foaming fixatives used for the mitigation of radioactive contamination need to meet before being implemented. Objective: Quantification of the mechanical properties of 6 commercial polyurethane foams was performed using multiple ASTM standards to record measurements for: Compression testing: Flexible Foams, Force required to produce 50% compression, Rigid Foams, Compressive and apparent modulus, Point of 10% core deformation, The 'Zero Deformation' point, Compressive strength, Yield point. Tensile testing: Tensile strength, Tensile stress, Percent elongation. Experiment 1: An electromechanical compressive tester (MTS Criterion Series 43) was used to evaluate 6 foams (4 flexible, 2 rigid). Per ASTM D1621, each rigid foam was compressed at 10% of the measured thickness per minute until the sample was 13% of it's original thickness. Per ASTM D35/4, each flexible foam was pre-flexed twice to 80% original thickness at a rate of 250 mm/min, then compressed to 80% original thickness at a rate of 50 mm/min. Experiment 2: An electromechanical tensile tester (MTS Criterion Series 43) was used to evaluate 6 foams (4 flexible, 2 rigid). Per ASTM D1623, rigid foams were pulled apart at a rate of 1.27 mm/min until the sample broke. Per ASTM D35/4, flexible foams were pulled apart at a rate of 500 mm/min until the sample broke. The strongest material in both compression and tensile testing scenarios was found to be the rigid intumescent polyurethane Hilti foam. The experiments revealed that the Hilti foam in a tensile scenario had a peak stress value that was larger than the closest competitor by a factor of 2.3 and a compressive yield point that was larger than the closest competitor by a factor of 1.4, indicating that the Hilti foam is the best choice for implementation in mechanically harsh environments. The performance metrics measured can serve as a basis for future mechanical tests that would help set relevant ASTM standards (E3191) for intumescent polyurethane foams in fixating applications. Tests like surface adhesion, impact, and flame tests would serve as a better indicator as to how this material would perform in environmentally harsh scenarios often found in decommissioned nuclear facilities. Further tests of the foams' intumescent properties would also be important should these foams be implemented in environmentally harsh scenarios.
12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
Homomorphic data compression for real time photon correlation analysis
The construction of highly coherent X-ray sources, combined with next-generation detectors that are larger and faster, has enabled new research opportunities across the scientific landscape. Among the techniques that benefit most from these advancements is X-ray photon correlation spectroscopy (XPCS), where faster acquisition unlocks the ability to study faster dynamics within samples. However, faster acquisition on larger detectors also introduces unprecedented challenges for online data processing and offline data storage. Such challenges are particularly prominent for XPCS, where real time analyses require simultaneous calculation of all the previously acquired data in the time series. We present a homomorphic compression scheme to effectively reduce the computational time and memory space required for XPCS analysis. Leveraging similarities in the mathematical expression between a matrix-based compression algorithm and the correlation calculation, our approach allows direct operation on the compressed data without their decompression. The offline compression scheme extends storage capacity by a factor of 40 while preserving key features in the lossy compressed data. Meanwhile, the online compression scheme reduces the computational time to below 1 ms, enabling real time calculation of the correlation functions at kHz framerate. Our demonstration of a homomorphic compression of scientific data provides an effective solution to the big data challenge at coherent light sources. Beyond the example shown in this work, the framework can be extended to facilitate real-time operations directly on a compressed data stream for other techniques.
36 MATERIALS SCIENCE↗
A graphics processing unit accelerated sparse direct solver and preconditioner with block low rank compression
We present the GPU implementation efforts and challenges of the sparse solver package STRUMPACK. The code is made publicly available on github with a permissive BSD license. STRUMPACK implements an approximate multifrontal solver, a sparse LU factorization which makes use of compression methods to accelerate time to solution and reduce memory usage. Multiple compression schemes based on rank-structured and hierarchical matrix approximations are supported, including hierarchically semi-separable, hierarchically off-diagonal butterfly, and block low rank. Here, in this paper, we present the GPU implementation of the block low rank (BLR) compression method within a multifrontal solver. Our GPU implementation relies on highly optimized vendor libraries such as cuBLAS and cuSOLVER for NVIDIA GPUs, rocBLAS and rocSOLVER for AMD GPUs and the Intel oneAPI Math Kernel Library (oneMKL) for Intel GPUs. Additionally, we rely on external open source libraries such as SLATE (Software for Linear Algebra Targeting Exascale), MAGMA (Matrix Algebra on GPU and Multi-core Architectures), and KBLAS (KAUST BLAS). SLATE is used as a GPU-capable ScaLAPACK replacement. From MAGMA we use variable sized batched dense linear algebra operations such as GEMM, TRSM and LU with partial pivoting. KBLAS provides efficient (batched) low rank matrix compression for NVIDIA GPUs using an adaptive randomized sampling scheme. The resulting sparse solver and preconditioner runs on NVIDIA, AMD and Intel GPUs. Interfaces are available from PETSc, Trilinos and MFEM, or the solver can be used directly in user code. We report results for a range of benchmark applications, using the Perlmutter system from NERSC, Frontier from ORNL, and Aurora from ALCF. For a high frequency wave equation on a regular mesh, using 32 Perlmutter compute nodes, the factorization phase of the exact GPU solver is about 6.5× faster compared to the CPU-only solver. The BLR-enabled GPU solver is about 13.8× faster than the CPU exact solver. For a collection of SuiteSparse matrices, the STRUMPACK exact factorization on a single GPU is on average 1.9× faster than NVIDIA’s cuDSS solver.
97 MATHEMATICS AND COMPUTING↗
Discrepancy in the Critical State Void Ratio of Poorly Graded Sand due to Shear Strain Localization
The critical state (CS) concept is a theoretical framework that models the constitutive behavior of soils, including sand and other granular materials. It supports the notion of a unique postfailure state, where the soil ultimately experiences continuous shearing with no change in the plastic volumetric strain. However, the published literature has frequently noted the nonconvergence of sand specimens with different initial densities to a unique CS in the compression plane due to many factors such as specimen fabric, particle morphology, breakage, and grain size distribution. This paper examines the CS for poorly graded (uniform) glass beads and 3 different types of silica sands using 50 conventional triaxial compression (CTC) experiments, 12 oedometer tests, and in situ synchrotron microcomputed tomography (SMT) scans for 10 CTC experiments. The results of the 50 CTC experiments revealed a diffused CS zone in the compression plane, which was further examined using the in situ SMT scans. A thorough three-dimensional image analysis of the SMT scans accurately quantified the evolution of the local void ratio (elocal) versus axial compression within zones of intensive shearing toward the center of the specimen. The evolution of the void ratio was also measured using the entire volume of the specimen (eglobal). At the CS, the elocal/eglobal ratio was assessed to be ~1.25 when a single shear band developed within the scanned specimens and ~1.1-1.15 for specimens that failed via external bulging that was internally manifested by the development of multiple shear bands. This finding suggests that the CS zone in the compression plane can be attributed to the common wrong consideration of eglobal evolution in lieu of elocal within the developing shear bands. Furthermore, the lack of shear band development in uniaxial compression has made the results of the oedometer test reliable in quantifying the CS parameters in the compression plane.
58 GEOSCIENCES↗
Extreme Compression for Large Scale Data Store
For the last 5 years Accelogic pioneered and perfected a radically new theory of numerical computing codenamed “Compressive Computing”, which has an extremely profound impact on real-world computer science [1]. At the core of this new theory is the discovery of one of its fundamental theorems which states that, under very general conditions, the vast majority (typically between 70% and 80%) of the bits used in modern large-scale numerical computations are absolutely irrelevant for the accuracy of the end result. This theory of Compressive Computing provides mechanisms able to identify (with high intelligence and surgical accuracy) the number of bits (i.e., the precision) that can be used to represent numbers without affecting the substance of the end results, as they are computed and vary in real time. The bottom line outcome would be to provide a state-of-the-art compression algorithm that surpasses those currently available in the ROOT framework, with the purpose of enabling substantial economic and operational gains (including speedup) for High Energy and Nuclear Physics data storage/analysis. In our initial studies, a factor of nearly x4 (3.9) compression was achieved with RHIC/STAR data where ROOT compression managed only x1.4.In this contribution, we will present our concepts of “functionally lossless compression”, have a glance at examples and achievements in other communities, present the results and outcome of our current, ongoing R&D, as well as present a high-level view of our plan to move forward with a ROOT implementation that would deliver a basic solution readily integrated into HENP applications. As a collaboration of experimental scientists, private industry, and the ROOT Team, our aim is to capitalize on the substantial success delivered by the initial effort and produce a robust technology properly packaged as an open-source tool that could be used by virtually every experiment around the world as means for improving data management and accessibility.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Progress toward Accelogic compression in ROOT
For the last 7 years, Accelogic pioneered and perfected a radically new theory of numerical computing codenamed “Compressive Computing”, which has an extremely profound impact on real-world computer science [1]. At the core of this new theory is the discovery of one of its fundamental theorems which states that, under very general conditions, the vast majority (typically between 70% and 80%) of the bits used in modern large-scale numerical computations are absolutely irrelevant for the accuracy of the end result. This theory of Compressive Computing provides mechanisms able to identify (with high intelligence and surgical accuracy) the number of bits (i.e., the precision) that can be used to represent numbers without affecting the substance of the end results, as they are computed and vary in real time. The bottom-line outcome will be to provide state-of-the-art compression algorithms --and accompanying software libraries-- able to surpass the performance of the compression engines currently available in the ROOT [7] framework. The resulting technology has the capability to enable substantial economic and operational gains (including speedup) for High Energy and Nuclear Physics data storage/analysis. In our initial studies, a factor of nearly x4 (3.9) compression was achieved with RHIC/STAR data where ROOT compression managed only x1.4 [6].As a collaboration of experimental scientists, private industry, and the ROOT Team, our aim is to capitalize on the substantial success delivered by the initial effort and produce a robust technology properly packaged as an open-source tool that could be used by virtually every experiment around the world as means for improving data management and accessibility.In this contribution, we will present our efforts integrating our concepts of “functionally lossless compression” within the ROOT framework implementation, with the purpose of producing a basic solution readily integrated into HENP applications. We will also present our progress applying this compression through realistic examples of analysis from both the STAR and CMS experiments.
Canal, Ph.↗
Effect of Moisture and Feedstock Variability on the Rheological Behavior of Corn Stover Particles
Continuous feeding, processing, and handling of biomass powders is pivotal to the economic viability of integrated biorefineries. However, current challenges associated with the operational reliability of bulk solids handling and transport greatly impact the process economics and ultimately the widespread commercialization of integrated biorefineries. In this work, we examine the effect of moisture and feedstock variability on the flow behavior of corn stover biomass particles. The total flow energy, compressibility, shear properties, and wall friction angles were measured for corn stover samples A and B containing 0%, 15%, 25%, 50%, and 75% (mass fraction) moisture contents using a FT4 powder rheometer. In general, the flowability of both A and B was reduced when moisture was present as indicated by the stability and variable flow rate, compressibility, and shear tests. The 15% moisture sample had the highest flow energy, revealing the interplay between the increased surface tension and looser packing both of which were caused by liquid bridging. The 75% moisture sample had the highest compressibility and the lowest flowability factor (ff c ). The trend in wall friction angle was found to be dependent upon the surface hydrophobicity of the wall material. The wall friction angle on the hydrophilic, stainless-steel surface increased with moisture and therefore poses additional handling challenges. On the other hand, the wall friction angle on the hydrophobic polymer surfaces were the lowest for samples with intermediate moisture contents. Sample B had greater bulk density, smaller compressibility, and greater flowability than the sample A , as suggested by the compressibility and shear tests. The wall friction angle of sample B was higher than or equal to that of sample A depending on both the type and the surface roughness of the wall material.
09 BIOMASS FUELS↗
Coordination changes in liquid tin under shock compression determined using in situ femtosecond x-ray diffraction
Little is known regarding the liquid structure of materials compressed to extreme conditions, and even less is known about liquid structures undergoing rapid compression on nanosecond timescales. Here, we report on liquid structure factor and radial distribution function measurements of tin shock compressed to 84(19) GPa. High-quality, femtosecond x-ray diffraction measurements at the Linac Coherent Light Source were used to extract the liquid diffuse scattering signal. From the radial distribution function, we find that the structural evolution of the liquid with increasing pressure mimics the evolution of the solid phase. With increasing pressure, we find that the liquid structure evolves from a complex structure, with a low coordination number, to a simple liquid structure with a coordination number of ~12. We provide a pathway for future experiments to study liquids at elevated pressures using high-energy lasers to shock compress materials beyond the reach of static diamond anvil cell techniques.
X-ray diffraction↗
Spaced Orthocyclic Winding Pattern for Improved Die Compressed Coils
Winding fill, coil fill, or packing factor, is critical to the performance of many types of electromagnetic devices from electric machines to actuators, transformers, inductors, etc. One technique to increase the packing of conductors in the winding area is to compress round magnet wire conductors in a die, so that they deform and air voids between coil turns are removed. The winding layout of the coil turns on a bobbin before compression has a direct impact on the achievable coil fill and the damage to the coil insulation or conductors after compression. It is thought that winding the turns in an orthocyclic pattern on the bobbin before compression would minimize damage to the turns. This however causes turns to displace during the winding compression process. A new pre-compression winding pattern is proposed, a spaced orthocyclic winding, to minimize turn displacement and damage. A design of experiments study has been carried out to optimize this winding method and to compare it to the traditional orthocyclic winding.
33 ADVANCED PROPULSION SYSTEMS↗
Ablating Pellets for Areal Density Symmetry Control in Indirectly Driven Inertial Confinement Fusion Target Designs
Asymmetry of the compressed fusion fuel configuration, as characterized by a variation in areal density, in inertial confinement fusion (ICF) implosions is a known performance limitation. While ignition has been achieved on the National Ignition Facility (NIF), implosion symmetry was not perfect, and subsequent experiments that have extended fusion performance still struggle with symmetry control and would, in principle, perform much better if the areal density uniformity of the fuel was improved by eliminating thin regions in the DT (Deuterium-Tritium) fuel that develop by the time of peak compression. In indirect drive (IDD), altering the time-dependent laser cone fraction and cross-beam energy transfer are the usual tactics for controlling implosion asymmetry in IDD target designs, but it may be that those tools have been pushed to their limits. In this Letter, we present a new tactic for controlling the implosion using satellite ablating pellets to manipulate the time-dependent x-ray field around an ICF capsule. In theory and in simulations the action of these pellets is quite effective at reducing the swing of x-ray drive from waist hot to pole hot, thus achieving a more one-dimensional fuel configuration at peak compression. Simulations show that a factor of ∼2 in yield performance can be obtained via this technique, when applied to past record implosions on the NIF. As a result, some speculative alternate uses of the proposed target concept are also mentioned in the conclusion of this Letter.
Implosion symmetry↗
Bit-GraphBLAS: Bit-Level Optimizations of Matrix-Centric Graph Processing on GPU
In the graph data structure like adjacency matrix, the connectivity of two nodes can be sufficiently represented using only 1 bit, but they are generally treated as 32-bit full-precision in state-of-the-art graph frameworks to adopt common sparse format such as CSR. Meanwhile, bit-level parallelism has recently be explored to have high-performance potential and low storage requirement on GPUs with dense bit-tiles. To fill the gap, our solution is a hierarchical storage format that contains the bit-indexing base and dense bit-tile units. Inherently, the granularity of the bit-tile is an essential factor in achieving both storage compression and GPU parallelism. How to find a sweet spot that trades off between avoiding sparsity and exploiting is comprehensively researched in this work. In the experiment, we evaluate the proposed storage format and algorithms on modern generation GPUs, including Pascal and Volta, to figure out critical software co-designs in conjunction with existing hardware-specific optimization.
Chen, Jou-An↗
IACMI Project 4.7: Pultruded Textile Carbon Fiber for Spar Caps (Final Report)
The primary objective of this project was to demonstrate the potential to significantly reduce the cost of wind turbine blades with carbon fiber reinforced polymer (CFRP) structure. Applicability of textile carbon fibers (TCF) were evaluated for use in pultruded spar cap (SC) elements as a path to cost reduction for utility scale wind turbine blades. In earlier work for the Department of Energy (DOE) Wind Energy Technologies Office (WETO), a collaboration of Sandia National Laboratory (SNL), Oak Ridge National Laboratory (ORNL), and Montana State University has demonstrated potential for pultruded TCF to compete with infused fiberglass and commercially available carbon fiber pultruded sections for spar cap construction. In the design cases evaluated, the TCF sections fared well when compared on cost per unit composite stiffness and cost per unit composite compressive strength for those designs [1]. Both stiffness and compressive strength tend to be key factors in the design of blade composite Spar Cap which carry the bulk of the blade structural loads in bending. Spar Cap design tends to distribute largely symmetric tensile and compressive stresses to opposite sides of the spar structure, but since carbon fiber composite compressive strength is typically 20-50% lower than tensile strength, the compressive loading reaches failure levels well before the tensile loading. Stiffness is critical in containing the large tip deflection in high wind loading situations. However, materials and process development were very limited in the earlier study and the work in this project was expanded to make the comparative information more representative of what will be required in order to make further inroads towards implementation. Similar to that study, this project team confirmed that the primary materials of interest for pultruded spar cap elements should be thermoset (TS) resins reinforced by carbon fibers, utilizing as high a percentage of TCF as practical to benchmark cost and performance against commercial carbon fibers. To make the closest comparison possible and eliminate specific test article size, resin selection, and equipment/operational nuance effects, the team planned to pultrude sections with 100% commercially available carbon fiber (Panex 35 carbon fiber from Zoltek) as well as samples utilizing high fractions of TCF. The resin system chosen was based on formulations recommended by large wind industry supplier Hexion and consisted of Hexion resin RSL-4597, curing agent CCA-138, and internal mold release additive 117, along with common kaolin filler ASP400P from BASF. As commonly deployed in spar cap configurations, the team had a mold built to pultrude a rectangular spar cap element of 100mm width and 3mm thickness. The extremely limited number of samples produced for the earlier study were produced with a “generic” epoxy utilized for a variety of applications by the pultruder contracted to produce test articles for demonstration purposes. More importantly, those samples were produced at a fiber fraction only slightly over 50%. Based on feedback from our industrial advisory team for that project and strongly recommended by this project team, the consensus is that it is highly desirable to obtain fiber fractions of 65-68% for significant penetration in wind blade spar cap. Although this requirement has yet to be exhaustively confirmed in readily available information, this was established as a project goal and informally decided we needed to exceed 60% fiber fraction to gain serious industry consideration. Previous TCF pultrusion trials have been challenged by the lack of robust TCF packages, resulting in non-uniform tension across and between tows, as well as excess labor and waste for removal of interleaved paper. The non-uniform tension and associated intermingling of tows in textile acrylic fiber tows and associated difficulties created from broken filaments in carbon fiber conversion inhibit the ordered packing necessary to enhance fiber fraction elevation. (These “cross-overs” are not considered undesirable for textile applications and there is some sense that they might be advantageous for those applications). In addition to work that is ongoing at the acrylic fiber manufacturers to improve their formats, The Institute of Advanced Composites Manufacturing Innovation (IACMI) Project 6.12 (report PA16-0349-6.12-01) [2] has developed and demonstrated a more robust packaging and creeling approach that at least partially addresses these issues, thus improving control of the TCF feed into the pultrusion unit. It was hoped that these and other improvements currently being implemented would allow us to achieve fiber fractions at least approaching these fiber fraction targets. During this project, sections utilizing 100% commercially available carbon fiber reinforcement were produced as a baseline, as well as sections reinforced with about 94% TCF and the balance being commercially available fiber for comparison. The most important finding was that similar to results reported in the earlier WETO-funded project and results from tests of TCF reported at IACMI meetings, this work demonstrated that sections pultruded with TCF in an epoxy resin frequently utilized in actual spar cap production had stiffness and compressive strengths largely comparable to similar sections pultruded with a commercially available carbon fiber also frequently utilized in the wind industry. Although the amount of that data is limited, some of the tensile strength results were actually closer than would have been expected based on fiber strength results provided by the TCF and commercial fiber producers. The actual test data are reported and discussed in detail in Section 5. The pultruded sections dominated by TCF reinforcement were approximately 8-10% lower in fiber fraction than for the sections produced using commercial fiber alone, making direct comparison difficult. The COVID-19 project has provided significant insight into the current state-of-the-art with various TCF product forms. The data obtained in this project will guide the planned improvements at the precursor level, especially in attaining uniform tensioning and payout to facilitate enhanced fiber fractions and overall processability of the TCF composites. The project team is providing guidance to stakeholders concerning the attributes, needs, and potential demand for TCF in wind blade spar caps. Results achieved in this project are consistent with findings in the related work cited [1] and support this guidance and the high potential for this product type. TCF precursor-producing partners continue to express interest in enhancing their product forms and the team looks forward to working with these improved materials as they become available.
42 ENGINEERING↗
Compressive strength improvements from noncircular carbon fibers: A numerical study
The benefits of high-performance unidirectional carbon fiber composites are limited in many cost-driven industries due to the high cost relative to alternative reinforcement fibers. Low-cost carbon fibers have been previously proposed, but the longitudinal compressive strength continues to be a limiting factor or studies are based on simplifications that warrant further analysis. A micromechanical model is used to (1) determine if the longitudinal compressive strength of composites can be improved with noncircular carbon fiber shapes and (2) characterize why some shapes are stronger than others in compression. In comparison to circular fibers, the results suggest that the strength can be increased by 10%–13% by using a specific six-lobe fiber shape and by 6%–9% for a three-lobe fiber shape. A slight increase is predicted in the compressive strength of the study two-lobe fiber but has the highest uncertainty and sensitivity to fiber orientation and misalignment direction. The underlying mechanism governing the compressive failure of the composites was linked to the unique stress fields created by the lobes, particularly the pressure stress in the matrix. Importantly, this work provides mechanics-based evidence of strength improvements from noncircular fiber shapes and insight on how matrix yielding is altered with alternative fiber shapes.
36 MATERIALS SCIENCE↗
Deep learning for in situ data compression of large turbulent flow simulations
As the size of turbulent flow simulations continues to grow, in situ data compression is becoming increasingly important for visualization, analysis, and restart checkpointing. For these applications, single-pass compression techniques with low computational and communication overhead are crucial. In this paper we present a deep-learning approach to in situ compression using an autoencoder architecture that is customized for three-dimensional turbulent flows and is well suited for contemporary heterogeneous computing resources. The autoencoder is compared against a recently introduced randomized single-pass singular value decomposition (SVD) for three different canonical turbulent flows: decaying homogeneous isotropic turbulence, a Taylor-Green vortex, and turbulent channel flow. Our proposed fully convolutional autoencoder architecture compresses turbulent flow snapshots by a factor of 64 with a single pass, allows for arbitrarily sized input fields, is cheaper to compute than the randomized single-pass SVD for typical simulation sizes, performs well on unseen flow configurations, and has been made publicly available. The results reported here show that the autoencoder dramatically outperforms a randomized single-pass SVD with similar compression ratio and yields comparable performance to a higher-rank decomposition with an order of magnitude less compression in regard to preserving a number of important statistical quantities such as turbulent kinetic energy, enstrophy, and Reynolds stresses.
97 MATHEMATICS AND COMPUTING↗
AEflow (Autoencoder fluid flow compression network) [SWR-22-29]
As the size of turbulent flow simulations continues to grow, in situ data compression is becoming increasingly important for visualization, analysis, and restart checkpointing. For these applications, single-pass compression techniques with low computational and communication overhead are crucial. In this paper we present a deep-learning approach to in situ compression using an autoencoder architecture that is customized for three-dimensional turbulent flows and is well suited for contemporary heterogeneous computing resources. The autoencoder is compared against a recently introduced randomized single-pass singular value decomposition (SVD) for three different canonical turbulent flows: decaying homogeneous isotropic turbulence, a Taylor-Green vortex, and turbulent channel flow. Our proposed fully convolutional autoencoder architecture compresses turbulent flow snapshots by a factor of 64 with a single pass, allows for arbitrarily sized input fields, is cheaper to compute than the randomized single-pass SVD for typical simulation sizes, performs well on unseen flow configurations, and has been made publicly available. The results reported here show that the autoencoder dramatically outperforms a randomized single-pass SVD with similar compression ratio and yields comparable performance to a higher-rank decomposition with an order of magnitude less compression in regard to preserving a number of important statistical quantities such as turbulent kinetic energy, enstrophy, and Reynolds stresses.
King, Ryan↗
Two-Phase Fluid Flow Properties of Rough Fractures With Heterogeneous Wettability: Analysis With Lattice Boltzmann Simulations
Fractures are conduits for fluid flow in low-permeability geological formations. Multiphase flow properties of fractures are important in natural processes and in engineering applications such as the evaluation of the sealing capacity of caprocks and productivity of hydrocarbon-bearing tight rocks. Investigations of flow and transport through fractures typically focus on the effects of fracture geometric and mechanical factors such as aperture, roughness, and compressibility. The wettability of the fracture surfaces and its influence on microscale interfacial phenomena and macroscale effective transport properties are seldom studied. In this study, we investigated the effect of heterogeneous wetting properties on the displacement of water by supercritical CO2 through a series of lattice Boltzmann method simulations. The results show the evolution of the CO2 plume within a fracture is controlled by both the roughness of the aperture field and the wetting distribution. We combined these factors into a capillary pressure map that can be related to the macroscopic flow behavior of the fracture. We observed that heterogeneous wetting distributions promote the residual trapping of water where lower capillary pressures allowed for isolated water pockets in higher capillary pressure zones. Analysis of fracture unsteady relative permeability shows the effect of wetting on permeability evolution and provides support for the viscous-coupling relative permeability model. Finally, analysis of the steady-state relative permeability and saturation demonstrates a strong correlation between permeability and the standard deviation of the capillary pressure field. Thus, characterizing the distribution of wetting properties of fractures is crucial to understanding multiphase fracture flow and transport properties.
58 GEOSCIENCES↗
Nonnegative canonical tensor decomposition with linear constraints: nnCANDELINC
Abstract There is an emerging interest for tensor factorization applications in big‐data analytics and machine learning. To speed up the factorization of extra‐large datasets, organized in multidimensional arrays (also known as tensors), easy to compute compression‐based tensor representations, such as, Tucker and tensor train formats, are used to approximate the initial large‐tensor. Further, tensor factorization is used to extract latent features that can facilitate discoveries of new mechanisms and signatures hidden in the data, where the explainability of the latent features is of principal importance. Nonnegative tensor factorization extracts latent features that are naturally sparse and parts of the data, which makes them easily interpretable. However, to take into account available domain knowledge and subject matter expertise, often additional constraints need to be imposed, which lead us to canonical decomposition with linear constraints (CANDELINC), a canonical polyadic decomposition with rank deficient factors. In CANDELINC, Tucker compression is used as a preprocessing step, which lead to a larger residual error but to more explainable latent features. Here, we propose a nonnegative CANDELINC (nnCANDELINC) accomplished via a specific nonnegative Tucker decomposition; we refer to as minimal or canonical nonnegative Tucker. We derive several results required to understand the specificity of nnCANDELINC, focusing on the difficulties of preserving the nonnegative rank of a tensor to its Tucker core and comparing the real valued to nonnegative case. Finally, we demonstrate nnCANDELINC performance on synthetic and real‐world examples.
97 MATHEMATICS AND COMPUTING↗