Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “data compression”

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,279 records · Page 71

Compressive strength of the plastic-bonded explosive surrogate idoxuridine under quasistatic and dynamic compression

Inert simulant materials, or “mocks”, are often used as surrogates for plastic-bonded explosives (PBX) in non-detonative tests in order to mitigate hazards. Mocks should reproduce as many non-detonative properties of the explosive as possible, including structural behavior in a variety of thermal and mechanical environments. Recently, the molecular crystal idoxuridine (IDOX) has been identified as an ideal mock for the main component in the explosive polymer-matrix composite PBX 9501, and has performed favorably under quasistatic loading conditions. Here, in order to assess robustness over a range of mechanical environments, plastic-bonded IDOX was compression tested from 0.001/s to 2000/s strain rates and compared to PBX 9501 historical data. Plastic-bonded IDOX showed good agreement to PBX 9501 across these strain rates, justifying continued development and production as a mock.

Burch, Alexandra C.↗

A Multiwavelength Investigation of PSR J2229+6114 and its Pulsar Wind Nebula in the Radio, X-Ray, and Gamma-Ray Bands

Abstract G106.3+2.7, commonly considered to be a composite supernova remnant (SNR), is characterized by a boomerang-shaped pulsar wind nebula (PWN) and two distinct (“head” and “tail”) regions in the radio band. A discovery of very-high-energy gamma-ray emission (E γ > 100 GeV) followed by the recent detection of ultrahigh-energy gamma-ray emission (E γ > 100 TeV) from the tail region suggests that G106.3+2.7 is a PeVatron candidate. We present a comprehensive multiwavelength study of the Boomerang PWN (100″ around PSR J2229+6114) using archival radio and Chandra data obtained two decades ago, a new NuSTAR X-ray observation from 2020, and upper limits on gamma-ray fluxes obtained by Fermi-LAT and VERITAS observatories. The NuSTAR observation allowed us to detect a 51.67 ms spin period from the pulsar PSR J2229+6114 and the PWN emission characterized by a power-law model with Γ = 1.52 ± 0.06 up to 20 keV. Contrary to the previous radio study by Kothes et al., we prefer a much lower PWNB-field (B∼ 3μG) and larger distance (d∼ 8 kpc) based on (1) the nonvarying X-ray flux over the last two decades, (2) the energy-dependent X-ray size of the PWN resulting from synchrotron burn-off, and (3) the multiwavelength spectral energy distribution (SED) data. Our SED model suggests that the PWN is currently re-expanding after being compressed by the SNR reverse shock ∼1000 yr ago. In this case, the head region should be formed by GeV–TeV electrons injected earlier by the pulsar propagating into the low-density environment.

Astronomy & Astrophysics↗

Low-Carbon District Heating: Performance Modeling of Hybrid Solar, Heat Pump, and Thermal Storage Systems for District Thermal Energy in the United States

District heating requires thermal energy in the temperature range of 40 degrees C - 120 degrees C. Typically, the thermal energy input for these systems has largely been met through fossil energy. However, the temperature range is low enough that it presents an opportunity for low-carbon technologies such as solar thermal and electrified thermal generators like heat pumps to decarbonize the heat generation. In this paper, a heat pump model was applied to estimate the performance and economics of a real-world low-carbon district heating substation. This system is comprised of a flat plate solar collector field paired with a mechanical vapor compression heat pump and hot water thermal storage, augmented by gas-fired boilers. Plant data was used to tune the model and estimate the system's benefits in terms of both standard financial metrics (IRR and payback), and environmental metrics, including avoided CO2 emissions. The model is subsequently employed to estimate the technical and economic potential of solar + heat pump + ther-mal storage hybrid systems as retrofits for district heating systems in eight US Markets.

district heating↗

Crush Strengths of Stainless Steel Honeycomb for Packaging

Honeycomb structures are used as energy absorbers to dissipate kinetic energy. Their low density to high axial strength and inexpensive cost make them ideal for the energy absorption application. Mechanical properties of lightweight materials are well known compared to more dense materials like stainless steel. Determines the viability of using 304 Stainless Steel honeycomb in future SRNL radioactive material packages. Using a combination of the material properties and testing procedures both static and dynamic crush strength of stainless steel honeycomb will be determined. Both static and dynamic crush strengths of stainless steel honeycomb need to be quantified because strain rate plays a role in the material strength. The stainless steel honeycomb being tested will be cut into sections of about 2.0 inches in width, 2.3 inches in length and 2.5 inches in height. The cell size of the material is 0.25 in and the ribbon thickness 0.003 in. The material has a yield strength of 48 ksi and an ultimate strength of 103.2 ksi, with an elastic modulus of 29,300 ksi. Compression testing to determine the static crush strength will use an Instron compression test machine. This Instron machine uses a 20 ksi capacity load cell, while recording data at varying low strain rates. Dynamic impact testing will involve dropping a 43.8 lb impactor from a height around 5 feet onto the stainless steel honeycomb. The honeycomb structure will be on a rigid surface and the impactor will be dropped via a quick release hook. Dynamic crush strength is determined from the weight of the impactor, height of the drop and difference in height of the honeycomb after the drop.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Paul‐Mohr‐Coulomb failure surface of rock in the brittle regime

Abstract The Paul‐Mohr‐Coulomb failure criterion includes the intermediate principal stress σ II and friction angles at the limiting stress states of σ II = σ III and σ II = σ I , where σ I and σ III are major and minor principal stresses. Conventional triaxial compression ( σ II = σ III ), extension ( σ II = σ I ), and plane strain ( σ I ≠ σ II ≠ σ III ) experiments were performed on dry rock. The failure data were plotted in principal stress space, and material parameters were determined in the context of two internal friction angles and the theoretical uniform triaxial (all‐around equal) tensile strength. Assuming isotropy, the triaxial compression and extension results were used to construct a six‐sided pyramidal failure surface, and the extension friction angle was larger than the compression friction angle, a sufficient but not necessary condition of the intermediate stress effect. To capture the behavior of the rock in multiaxial loading, the Paul‐Mohr‐Coulomb criterion was extended to form a 12‐sided pyramid with best fit planes.

Makhnenko, Roman Y.↗

High pressure structure studies of 6H-SrIrO{sub 3} and the octahedral tilting in 3C-SrIrO{sub 3} towards a post-perovskite

The high pressure behaviors of the two perovskite structures (hexagonal 6H-SrIrO{sub 3} and orthorhombic 3C-SrIrO{sub 3}) have been studied in diamond anvil cells to 43 and 60 GPa, respectively, using synchrotron powder X-ray diffraction. 6H-SrIrO{sub 3} was first synthesized at ambient pressure and subsequently transformed into 3C-SrIrO{sub 3} in a large volume press at 8.8 GPa and 1000 °C. Both polymorphs were found to retain the initial symmetry up to the highest pressures measured, but in the case of 6H-SrIrO{sub 3}, two anomalies were identified: a change in the axial compressibilities at 24 GPa and a change in both the axial and volume compressibilities at 32 GPa. Fitting a 3rd order Birch-Murnaghan equation of state to the obtained P-V data yielded bulk moduli of K{sub 0}=151.5(12) GPa (fitted range 0<P<30 GPa) for 6H-SrIrO{sub 3} and K{sub 0}=187.1(9) GPa for 3C-SrIrO{sub 3}. Analysis of the structural parameters for 6H-SrIrO{sub 3} aided by F-f plots suggests the anomalies are caused by changes in the compression mechanism. Comparison of the two polymorphs reveals that 6H-SrIrO{sub 3} becomes less compressible than 3C-SrIrO{sub 3} above 32 GPa as a result of the mechanistic change, and a crossing of their P-V curves is avoided. For 3C-SrIrO{sub 3}, analysis of the octahedral tilt angles shows that these increase monotonically from the ambient value of 7.23(6) to 23.0(2)° at 60 GPa suggesting that a transition to a post-perovskite is approached. - Graphical abstract: The structure of 3C-SrIrO{sub 3} synthesized at high pressure has been studied in a diamond anvil cell and is found to maintain its orthorhombic structure up to 60 GPa but distorts drastically due to octahedral tilting, strongly suggesting that a post-perovskite transition is possible. The compression is compared to 6H-SrIrO{sub 3} synthesized at ambient pressure, which is initially found to be more compressible but undergoes changes in its compression mechanism and 3C remains denser throughout compression. Display Omitted - Highlights: • Structures of 6H- and 3C-SrIrO{sub 3} were studied to pressures of 43 and 60 GPa, respectively. • Both structures maintain the initial symmetry up to the highest pressures measured. • 6H-SrIrO{sub 3} is initially more compressible than 3C-SrIrO{sub 3}. • Due to changes in the compression mechanism in 6H-SrIrO{sub 3} the 3C-SrIrO{sub 3} remains denser throughout. • 3C-SrIrO{sub 3} show large octahedral tilting, strongly suggesting a post-perovskite transition is imminent.

36 MATERIALS SCIENCE↗

Public Reference Data for Megawatt-Scale Hydrogen Electrolysis - Simulated Marine Hydrokinetic Tidal Turbine

The U.S. Department of Energy and National Laboratory of the Rockies (NLR) demonstrate hydrogen electrolysis, hydrogen compression and storage, and variable hydrogen fuel cell power production using megawatt-scale equipment at NLR’s Flatirons Campus as part of the Advanced Research on Integrated Energy Systems (ARIES) initiative. This dataset is part of that effort and is intended for academic, national laboratory, industrial, and other stakeholders to plan, design, and validate models of megawatt-scale hydrogen technologies and diverse energy infrastructure. These data provide a baseline for how existing hydrogen electrolysis technologies perform when coupled with other energy technologies. This dataset contains inputs and outputs from simulations of a floating marine hydrokinetic turbine over approximately half a tidal cycle (~6.6 hours). Inflow conditions were derived from field measurements in Alaska’s Cook Inlet and represent a tidal environment in which the current speed ramps from near 0 m/s to a peak of 3 m/s and back. The original acoustic doppler current profiler dataset is publicly available on the Marine and Hydrokinetic Data Repository. In a full tidal cycle, the flow reverses and the rotor would reorient; this reversal was not modeled. In the Cook Inlet campaign , turbulence intensity was similar in both directions. Two inflow cases are included. In the first case, labeled “raw” in the files, the measured current time series was used directly in the InflowWind module of OpenFAST. Speed and direction were applied as a function of time and elevation, uniformly in the horizontal direction. With full spatial coherence, this approach captures high turbulent variability and results in pronounced power fluctuations, so it is considered a conservative, near-worst-case representation of loading. In the second case, labeled “average” in the files, a 30-minute moving average was applied to extract the slowly varying mean speed. The residual fluctuations about this mean were used to generate spatially varying, full-field turbulence inputs with TurbSim, giving a more physically realistic representation of the inflow across the rotor disk. Two random realizations were used to produce distinct inflow conditions for two OpenFAST simulations representing a two-turbine array. The same turbulence intensity is applied across the full time series, producing larger fluctuations at the start and end, where the mean speed is low. The second case is the more appropriate framework for performance and power assessment but overpredicts turbulence at lower flow speeds and underpredicts it at higher speeds. As the floating platform moves and the rotor changes its x-position, Taylor’s frozen turbulence hypothesis used by InflowWind assumes a constant rather than a time-varying mean velocity, introducing some inaccuracy in the velocity plane sampling. The turbine modeled is the 500-kW Reference Model 1, a horizontal-axis two-bladed hydrokinetic turbine on a four-column floating semisubmersible substructure . Simulations were performed using OpenFAST v4.1 with the Reference Open Source Controller (ROSCO) v2.10. All input files required to reproduce the simulations are included. The electrolyzer is a 1.25-MW proton exchange membrane type MC250 system manufactured by Nel . This unit supports up to 2.5 MW, but NLR has only a single 1.25-MW stack. The datasets report hydrogen balance-of-plant and system data, all captured at 1 Hz, including hydrogen mass production measured with an Emerson Coriolis flow meter. The system controls hydrogen production by varying direct current applied to the stack, from a maximum of 3,000 A to a minimum safe operating current of 300 A, or 10%. Because the current–voltage characteristic changes as the stack ages and efficiency degrades, the actual minimum safe operating power changes over time. The simulated tidal turbine time series data was translated from power (kilowatts) to current (amperes) using a curve fit with calibration data and sent to the electrolyzer power supply at 1-Hz. Each zip file represents a single tidal electrolysis experiment and is named: {technology}_{inflow method}_{number of 500 kW tidal turbines connected} For instance, “tidal-500kW-RM1_average_2.zip” is a 6-hour experiment using the 500-kW tidal reference model, scaled by 2x (1-MW) to better match the electrolyzer maximum of 1.25MW, fed with the 30-minute moving average current case. Each zip folder contains the following files: A .csv file of raw data. An .xlsx file explaining all the fields in the raw data. A .png plot showing the time series of hydrogen production in kilograms per hour, electrolysis power consumption, and input wave power. A .csv file combines all tidal profiles as "combined_tidal_experiments.csv." A separate experiment, “characterization_200.zip,” shows the MC250 electrolyzer steady-state response with 30-minute load steps over 5 hours and is accessible with this entry.

08 HYDROGEN↗

Multifidelity validation of digital surrogates using variable-density turbulent mixing models

In this study, ensembles of experimental data are presented and utilized to compare and validate two models used in the simulation of variable-density (Atwood = 0.22), compressible turbulent mixing. Though models of this kind (Reynolds averaged NavierStokes and large-eddy simulations) have been validated extensively with more canonical flows in previous studies, here the present approach offers novelty in the complexity of the geometry, the ensemble-based validation, and the uniformity of the computational framework on which the models are tested. Moreover, all experimental and computational tasks were completed by the authors which has led to a tightly coupled experimental configuration with its “digital twin.” The experimental divergent-shock-tube facility and its data acquisition methods are described and replicated in simulation space. A 2D Euler model which neglects the turbulent mixing at the interface is optimized to experimental data using a Gaussian process. This model then serves as the basis for both the 2D RANS and 3D LES studies that make comparisons to the mixing-layer data from the experiment. A relatively simple RANS model is shown to produce good agreement with experimental data only at late flow development times. The LES ensembles generally show good agreement with experimental data but display sensitivity to the characterization of initial conditions. Resolution-dependent behavior is also observed for certain higher-order statistics of interest. Overall, the LES model successfully captures the effects of divergent geometry, compressibility, and combined nonlinear instabilities inherent to the problem. The successful prediction of mixing width and its growth rate highlight the existence of three distinct regimes in the development of the instability, each with similarities to previously studied instabilities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Uniaxial Compression of Cellular Materials at a 10 -1 s -1 Strain Rate Simultaneously with Synchrotron X-ray Computed Tomographic Imaging

The topic is presented as a series of slides. Motivation for the work included the following: X-ray tomography is a fantastic technique for characterizing a material’s starting structure as well as for non-destructive, in situ experiments to investigate material response; 3D X-ray tomography is needed to fully characterize the morphology of cellular materials; and synchrotron micro-CT can capture 3D images without pausing experiment. Among the conclusions reached are these: High-rate radiographic and tomographic imaging (0.25 s 3D frame rate) using synchrotron CT can capture full 3D images of hyper-elastic materials at a 10 -2 strain rate; dynamic true in situ uniaxial loading can be accurately captured; the three stages of compression can be imaged: bending, buckling, and breaking; implementation of linear modeling is completed; meshes have been imported into LANL modeling codes--testing and validation is underway and direct comparison and validation between in situ data and modeled mechanical response is possible.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Ultrabright x-ray laser scattering for dynamic warm dense matter physics

In megabar shock waves, materials compress and undergo a phase transition to a dense charged-particle system that is dominated by strong correlations and quantum effects. This complex state, known as warm dense matter, exists in planetary interiors and many laboratory experiments (for example, during high-power laser interactions with solids or the compression phase of inertial confinement fusion implosions). Here, we apply record peak brightness X-rays at the Linac Coherent Light Source to resolve ionic interactions at atomic (ångström) scale lengths and to determine their physical properties. Our in situ measurements characterize the compressed lattice and resolve the transition to warm dense matter, demonstrating that short-range repulsion between ions must be accounted for to obtain accurate structure factor and equation of state data. Additionally, the unique properties of the X-ray laser provide plasmon spectra that yield the temperature and density with unprecedented precision at micrometre-scale resolution in dynamic compression experiments.

laser-produced plasmas↗

Performance assessment of 700-bar compressed hydrogen storage for light duty fuel cell vehicles

In this study, type 4 700-bar compressed hydrogen storage tanks were modeled using ABAQUS. The finite element model was first calibrated against data for 35-L subscale test tanks to obtain the composite translation efficiency, and then applied to full sized tanks. Two variations of the baseline T700/epoxy composite were considered in which the epoxy was replaced with a low cost vinyl ester resin and low cost resin with an alternate sizing. The results showed that the reduction in composite weight was attributed primarily to the lower density of the resin and higher fiber volume fraction in the composite due to increased squeeze-out with the lower viscosity vinyl ester resin. The system gravimetric and volumetric capacities for the onboard storage system that holds 5.6 kg H 2 are 4.2 wt% (1.40 kWh/kg) and 24.4 g-H 2 /L (0.81 kWh/L), respectively. The system capacities increase and carbon fiber requirement decreases if the in-tank amount of unrecoverable hydrogen is reduced by lowering the tank "empty" pressure. Models of an alternate tank design showed potential 4-7% saving in composite usage for tanks with a length-to-diameter (L/D) ratio of 2.8-3.0 but no saving for L/D of 1.7. Lastly, a boss with smaller opening and longer flange does not appear to reduce the amount of helical windings.

08 HYDROGEN↗

A Parametric, Data-Driven, Non-Intrusive Reduced-Order Model Framework for Crystal Plasticity Simulations of Voids

The influence of the internal structure at micrometer length scales on the deformation of polycrystalline materials can be effectively captured using crystal plasticity finite element methods (CPFEM). However, the complexity and nonlinearity of the deformation equations CPFEM solves demand significant computational power and resources to achieve accurate predictions, limiting its broader application. To address this challenge, we have identified a reduced-order representation of the complex data in order to establish a computationally efficient reduced-order models (ROM) and drastically reduce the computational expense of CPFEM. Specifically, in this work, we developed a parametric, data-driven, and non-intrusive ROM framework for CPFEM using proper orthogonal decomposition (POD) and sparse variational Gaussian process (SVGP) regression for single-crystal microstructures under tensile loading conditions. The developed protocol enables one to compress field into a latent/low-dimensional space described by principal component analysis (PCA) via the singular value decomposition (SVD) algorithm. As a result, the high-dimensional data are reduced to a significantly smaller amount of dimensions with POD bases and POD coefficients. Furthermore, we deployed an ensemble of SVGPs—extended from the classical Gaussian process (GP) regression for scalability and handling big data—in a massively parallel manner to train and predict latent POD coefficients using known POD bases from a set of previously obtained simulations results. Lastly, using the predicted POD coefficients, we reconstructed the full-field results and showed reasonable agreement compared with the true values obtained from running CPFEM. The developed framework is validated with a set of CPFEM simulations of a single embedded void in single-crystal aluminum alloy. While the framework is broadly applicable, this work specifically focuses on single-crystal microstructures, a single load case (e.g., tensile), and a specific void geometry (spherical).

Anisotropy↗

Development, Verification, and Validation of an OpenFOAM-Based Solver for Modeling Inertial Fusion Energy Chambers

Our work seeks to introduce a computational tool tailored to the physics of inertial fusion energy chambers, in particular, those concepts based on thick liquid walls. In this approach, the structural materials are protected by several neutron mean-free-paths of renewable liquid and thus will be able to survive much longer than un-shielded walls, with virtually all structures lasting for the life of the plant and enabling the use of commercially available and qualified materials. The OpenFOAM-based solver named rhoCentralFoam has been used as a starting point. rhoCentralFoam belongs to the standard OpenFOAM solver toolset. It is a high-speed, explicit compressible flow solver with shock-capturing capability. While the main features have been retained, the solver had to be restructured to make use of tabular data for equations of states, a necessary addition to model the complex thermo-physical properties of ionized gasses. This entailed the need to change the independent state variables used by the solver, resulting in a new thermodynamic library and slightly different solution algorithm. Moreover, a radiation heat transfer model based on the P-1 approximation was added to the solver. The solver is verified against an analytical solution from the Sedov-Taylor-Neumann test problem to showcase the ability of the hydrodynamic solvers to handle strong shocks, whereas the P-1 model was verified using a simple one-dimensional problem with an analytical solution. Additionally, a validation case involving shock-wave propagation through jet array is presented, and the results are compared with experimental data from the open literature. Lastly, in order to showcase the utility of the solver for practical cases, we applied the refined solver to two representative scenarios: gas venting within the HYLIFE-II chamber and the compression of the gas following the partial ablation of the liquid wall.

Chamber dynamics↗

Evaluation of the robustness of the preprocessing technique improving reversible compressibility of CT images: Tested on various CT examinations

Purpose: To modify the preprocessing technique, which was previously proposed, improving compressibility of computed tomography (CT) images to cover the diversity of three dimensional configurations of different body parts and to evaluate the robustness of the technique in terms of segmentation correctness and increase in reversible compression ratio (CR) for various CT examinations.Methods: This study had institutional review board approval with waiver of informed patient consent. A preprocessing technique was previously proposed to improve the compressibility of CT images by replacing pixel values outside the body region with a constant value resulting in maximizing data redundancy. Since the technique was developed aiming at only chest CT images, the authors modified the segmentation method to cover the diversity of three dimensional configurations of different body parts. The modified version was evaluated as follows. In randomly selected 368 CT examinations (352 787 images), each image was preprocessed by using the modified preprocessing technique. Radiologists visually confirmed whether the segmented region covers the body region or not. The images with and without the preprocessing were reversibly compressed using Joint Photographic Experts Group (JPEG), JPEG2000 two-dimensional (2D), and JPEG2000 three-dimensional (3D) compressions. The percentage increase in CR per examination (CR{sub I}) was measured.Results: The rate of correct segmentation was 100.0% (95% CI: 99.9%, 100.0%) for all the examinations. The median of CR{sub I} were 26.1% (95% CI: 24.9%, 27.1%), 40.2% (38.5%, 41.1%), and 34.5% (32.7%, 36.2%) in JPEG, JPEG2000 2D, and JPEG2000 3D, respectively.Conclusions: In various CT examinations, the modified preprocessing technique can increase in the CR by 25% or more without concerning about degradation of diagnostic information.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

A discontinuous Galerkin spectral element method for compressible reacting flows

High-order methods have recently been shown to be an effective tool for high-fidelity flow computations like direct numerical simulations and large-eddy simulations because of their strong balance between accuracy and computational cost. In this work, a high-order discontinuous Galerkin spectral element method (DGSEM) is developed to solve the chemically reacting Navier-Stokes equations. To handle the disparate length and time scales associated with these equations, we develop a novel method which combines the spectral accuracy of the SEM with the flexibility of the DG approach. The framework, implemented in the spectral element code Nek5000, is well suited to capture turbulence in smooth regions of the flow, while maintaining numerical stability in the presence of shocks. An entropy-residual based artificial viscosity is added to smooth shocked regions of flow, and a positivity-preserving limiter is implemented to suppress non-physical oscillations. These enhancements support the numerical stability of the hydrodynamic sub-step, which is decoupled from the chemistry integration through a second-order operator splitting method. Here, a series of smooth and discontinuous validation cases are presented in increasing physical and computational complexity for both inviscid and viscous flows. In particular, simulations of canonical one-dimensional and two-dimensional detonations are performed, and the high-order numerical results are validated against available literature data. Additional validation studies are carried out for classical three-dimensional numerical simulations of incompressible and compressible turbulent flows.

Compressible reacting flows↗

Metadata compression

A computer-executable method, system, and computer program product for managing a Parallel-Logged File System (PLFS) on a data storage system, wherein the data storage system is in communication with an application that updates the data storage system with checkpoint data, the computer-executable method, system, and computer program product comprising receiving a first checkpoint including data from the application; analyzing the data in the first checkpoint to determine whether a pattern exists in the checkpoint data; and upon a positive determination, adding the pattern into PLFS.

97 MATHEMATICS AND COMPUTING↗

Alternative Fuels Data Center Fuel Properties Comparison

This chart compares the physical fuel properties and considerations associated with gasoline/E10, low sulfur diesel, biodiesel, propane, compressed natural gas, liquefied natural gas, ethanol, methanol, hydrogen, and electricity for use as vehicle fuels.

47 OTHER INSTRUMENTATION↗

Lossy checkpoint compression in full waveform inversion: a case study with ZFPv0.5.5 and the overthrust model

This paper proposes a new method that combines checkpointing methods with error-controlled lossy compression for large-scale high-performance full-waveform inversion (FWI), an inverse problem commonly used in geophysical exploration. This combination can significantly reduce data movement, allowing a reduction in run time as well as peak memory. In the exascale computing era, frequent data transfer (e.g., memory bandwidth, PCIe bandwidth for GPUs, or network) is the performance bottleneck rather than the peak FLOPS of the processing unit. Like many other adjoint-based optimization problems, FWI is costly in terms of the number of floating-point operations, large memory footprint during backpropagation, and data transfer overheads. Past work for adjoint methods has developed checkpointing methods that reduce the peak memory requirements during backpropagation at the cost of additional floating-point computations. Combining this traditional checkpointing with error-controlled lossy compression, we explore the three-way tradeoff between memory, precision, and time to solution. We investigate how approximation errors introduced by lossy compression of the forward solution impact the objective function gradient and final inverted solution. Empirical results from these numerical experiments indicate that high lossy-compression rates (compression factors ranging up to 100) have a relatively minor impact on convergence rates and the quality of the final solution.

58 GEOSCIENCES↗