SEARCH · Search NASA

Results for “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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 811 records · Page 45

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Compressing branch-and-bound trees

A branch-and-bound (BB) tree certifies a dual bound on the value of an integer program. In this work, we introduce the tree compression problem (TCP): Given a BB tree T that certifies a dual bound, can we obtain a smaller tree with the same (or stronger) bound by either (1) applying a different disjunction at some node in T or (2) removing leaves from T? Here we believe such post-hoc analysis of BB trees may assist in identifying helpful general disjunctions in BB algorithms. We initiate our study by considering computational complexity and limitations of TCP. We then conduct experiments to evaluate the compressibility of realistic branch-and-bound trees generated by commonly-used branching strategies, using both an exact and a heuristic compression algorithm.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2018

Effect of Compressive Mode I on the Mixed Mode I/II Fatigue Crack Growth Rate of 42CrMo4

Subsurface cracks in mechanical contact loading components are subjected to mixed mode I/II, so it is necessary to evaluate the fatigue behavior of materials under mixed mode loading. For this purpose, fatigue crack propagation tests are performed with compact tension shear specimens for several stress intensity factor (SIF) ratios of mode I and mode II. The effect of compressive mode I loading on mixed mode I/II crack growth rate and fracture surface is investigated. Tests are carried out for the pure mode I, pure mode II, and two different mixed mode loading angles. On the basis of the experimental results, mixed mode crack growth rate parameters are proposed according to Tanaka and Richard with Paris’ law. Results show neither Richard’s nor Tanaka’s equivalent SIFs are very useful because these SIFs depend strongly on the loading angle, but Richard’s equivalent SIF formula is more suitable than Tanaka’s formula. The compressive mode I causes the crack closure, and the friction force between the crack surfaces resists against the crack growth. In compressive loading with 45° angle, da/dN increases as K{sub eq} decreases.

36 MATERIALS SCIENCE↗

DOE OSTI2018

Investigation of Mild Steel Thin-Wall Tubes in Unfilled and Foam-Filled Triangle, Square, and Hexagonal Cross Sections Under Compression Load

This study has examined proposed structures with mild steel-reinforced LM30 aluminum (Al) alloy having diversely unfilled and 10 wt.% SiCp composite foam-filled tubes for improving axial compression performance. This class of material has novel physical, mechanical, and electrical properties along with low density. In the present experiment, Al alloy foams were prepared by the melt route technique using metal hydride powder as a foaming agent. Crash energy phenomena for diverse unfilled and foam-filled in mild steel thin-wall tubes (triangular, square and hexagonal) were studied as well. Compression deformation investigation was conducted at strain rates of 0.001-0.1/s for evaluating specific energy absorption (SEA) under axial loading conditions. The results were examined to measure plateau stress, maximum densification strain, and deformation mechanism of the materials. Specific energy absorption and total energy absorption capacities of the unfilled and filled sections were determined from the compressive stress–strain curves, which were then compared with each other.

36 MATERIALS SCIENCE↗

DOE OSTI2018

Effect of Plane Strain Compression and Subsequent Recrystallization Annealing on Microstructures and Phase Transformation of NiTiFe Shape Memory Alloy

The effect of plane strain compression and subsequent recrystallization annealing on microstructures and phase transformation of NiTiFe shape memory alloy (SMA) is investigated. Inhomogeneous plastic deformation at various deformation degrees occurs in NiTiFe SMA during plane strain compression. Nanocrystalline phase and amorphous phase increase as the deformation degree increases. B2 austenite, B19′ martensite, nanocrystalline and amorphous phases coexist in the NiTiFe samples subjected to large plastic strain. The static recrystallization mechanisms depend on the microstructures of the deformed NiTiFe samples. The static recrystallization mechanisms deal with nucleation and growth of the recrystallized grains, growth of nanocrystalline phase and crystallization of amorphous phase. Grain size, subgrain boundaries, geometrically necessary dislocation density and Schmid factor are captured on the basis of electron backscattered diffraction data. The process of recrystallization annealing cannot eliminate the deformation texture completely. The slip direction [110] is the most favorable slip direction in the recrystallized NiTiFe sample. Plane strain compression along with subsequent recrystallization annealing changes the phase transformation path of as-rolled NiTiFe SMA, and it results in the decreasing martensite transformation start temperature. The three annealed NiTiFe samples exhibit the similar phase transformation behavior since complete recrystallization annealing leads to the similar microstructures.

36 MATERIALS SCIENCE↗

DOE OSTI2015

Room temperature deformation mechanisms of alumina particles observed from in situ micro-compression and atomistic simulations.

Aerosol deposition (AD) is a solid-state deposition technology that has been developed to fabricate ceramic coatings nominally at room temperature. Sub-micron ceramic particles accelerated by pressurized gas impact, deform, and consolidate on substrates under vacuum. Ceramic particle consolidation in AD coatings is highly dependent on particle deformation and bonding; these behaviors are not well understood. In this work, atomistic simulations and in situ micro-compressions in the scanning electron microscope, and the transmission electron microscope (TEM) were utilized to investigate fundamental mechanisms responsible for plastic deformation/fracture of particles under applied compression. Results showed that highly defective micron-sized alumina particles, initially containing numerous dislocations or a grain boundary, exhibited no observable shape change before fracture/fragmentation. Simulations and experimental results indicated that particles containing a grain boundary only accommodate low strain energy per unit volume before crack nucleation and propagation. In contrast, nearly defect-free, sub-micron, single crystal alumina particles exhibited plastic deformation and fracture without fragmentation. Dislocation nucleation/motion, significant plastic deformation, and shape change were observed. Simulation and TEM in situ micro-compression results indicated that nearly defect-free particles accommodate high strain energy per unit volume associated with dislocation plasticity before fracture. As a result, the identified deformation mechanisms provide insight into feedstock design for AD.

36 MATERIALS SCIENCE↗

DOE OSTI2023

In-Situ X-Ray Imaging High Strain Rate Compression of Laminate Al-Graphene Composite and Mechanical Property Characterization

Owing to the opaque nature of the laminated structures, traditional high-speed optical camera cannot be used to detect the dynamic process of sub-surface deformation. In this article, we report a study of using high speed X-ray imaging to study the high strain rate deformation in laminated Al structures. We used a Kolsky bar apparatus to apply dynamic compression and a high-speed synchrotron X-ray phase contrast imaging (PCI) setup to conduct the in situ X-ray imaging study. The in situ X-ray imaging captures the shock wave propagation in the laminated structures. After shock compression, we characterized the microstructures by using transmission electron microscopy (TEM), which demonstrates an increase of dislocation density. In conclusion, the micro-pillar compression tests show that the yield strength at 0.2% offset of laminated Al-graphene composite has a significant increase of 67%, from 30 to 50 MPa, compared to laminate Al after shock loading.

36 MATERIALS SCIENCE↗

DOE OSTI2017

Cyclic compression response of micropillars extracted from textured nanocrystalline NiTi thin-walled tubes

Compression-compression cyclic deformation of nanocrystalline NiTi tubes intended for medical stents and with an outer diameter of 1 mm and wall thickness of 70 μm was studied using micropillars produced by FIB with the loading axis orthogonal to the tube axis. These micropillars were cycled in a displacement-controlled mode using a nanoindenter equipped with a flat punch to strain levels of 4, 6 and 8% in each cycle and specimens were subjected to several hundred cycles. Furthermore, the cyclic response of two NiTi tubes, one with Af of 17 °C and the other with an Af of -5 °C is compared. The texture of the tube with the Af of -5 °C was measured at the microscopic level using transmission electron microscopy and at the macroscopic level by X-ray diffraction and good agreement was noted. Characteristics such as i) a reduction in the forward transformation stress, ii) increase in maximum stress for a given displacement amplitude, and iii) a reduction in the hysteresis loop area, all with increasing number of cycles, observed typically during cyclic deformation of conventional macroscopic specimens, were captured in the micropillar cyclic tests. Our observations lead to the conclusion that micropillar compression testing in a cyclic mode can enable characterizing the orientation-dependent response in such small dimension components that see complex loading in service, and additionally provide an opportunity for calibrating constitutive equations in micromechanical models.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Investigating the compressive strength and strain localization of nanotwinned nickel alloys

Sputter deposited nickel-molybdenum-tungsten (Ni-Mo-W) thin films possess a beneficial suite of properties that stem from the extremely fine growth twins that form during the deposition process. Previously these materi- als were only characterized in tension, however, in this study in situ micropillar compression and post-mortem microstructural analysis of nanotwinned Ni 84 Mo 11 W 5 micropillars were employed to measure the compressive response and elucidate the attendant deformation mechanisms. The pillars exhibit Hookean behavior to com- pressive strengths of 3-3.5 GPa and the onset of non-linear plastic flow was manifest by discrete strain bursts and highly localized shear bands. Plastic deformation was concentrated at the top of the pillar, while the bulk of the micropillar was nominally unaffected. Post-mortem investigations indicate that at sufficiently high stresses shear banding is triggered, resulting in intense and highly localized plastic deformation that led to the formation of twin-free nanocrystalline grains within highly deformed shear bands. By contrast, the regions adjacent to the shear bands were unaffected. The absence of detwinning and dislocation glide mechanisms was unexpected and in direct contrast to what has been observed in nanotwinned Cu-Al. Post-mortem observations of the Ni-Mo-W micropillars suggest that the ultrafine twins create a unique form of dislocation starvation and source-limited plasticity. The ultrahigh compressive strength is governed by the triggering of shear bands rather than the activa- tion of dislocation glide. Finally, the specialized nature of plasticity in nanotwinned Ni-Mo-W is clear, even though the precise trigger for shear band formation in nanotwinned Ni-Mo-W remains to be identified.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Utilizing computer vision and artificial intelligence algorithms to predict and design the mechanical compression response of direct ink write 3D printed foam replacement structures

Additive Manufacturing (AM) of porous polymeric materials, such as foams, recently became a topic of intensive research due their unique combination of low density, impressive mechanical properties, and stress dissipation capabilities. Conventional methods for fabricating foams rely on complex and stochastic processes, making it challenging to achieve precise architectural control of structured porosity. In contrast, AM provides access to a wide range of printable materials, where precise spatial control over structured porosity can be modulated during the fabrication process enabling the production of foam replacement structures (FRS). Current approaches for designing FRS are based on intuitive understanding of their properties or an extensive number of finite element method (FEM) simulations. These approaches, however, are computationally expensive and time consuming. As such, in this work, we present a novel methodology for determining the mechanical compression response of direct ink write (DIW) 3D printed FRS using a simple cross-sectional image. By obtaining measurement data for a relatively small number of samples, an artificial neural network (ANN) was trained, and a computer vision algorithm was used to make inferences about foam compression characteristics from a single cross-sectional image. Finally, a genetic algorithm (GA) was used to solve the inverse design problem, generating the AM printing parameters that an engineer should use to achieve a desired compression response from a DIW printed FRS. The methods developed herein present an avenue for entirely autonomous design and analysis of additively manufactured structures using artificial intelligence.

36 MATERIALS SCIENCE↗

DOE OSTI2019

A New Computational Method for Membrane Compressibility: Bilayer Mechanical Thickness Revisited

Because lipid bilayers can bend and stretch in ways similar to thin elastic sheets, physical models of bilayer deformation have utilized mechanical constants such as the moduli for bending rigidity (κ C ) and area compressibility (K A ). However, the use of these models to quantify the energetics of membrane deformation associated with protein-membrane interactions, and the membrane response to stress is often hampered by the shortage of experimental data suitable for the estimation of the mechanical constants of various lipid mixtures. Although computational tools such as molecular dynamics simulations can provide alternative means to estimate K A values, current approaches suffer significant technical limitations. Here, in this study, we present a novel, to our knowledge, computational framework that allows for a direct estimation of K A values for individual bilayer leaflets. The theory is based on the concept of elasticity and derives K A from real-space analysis of local thickness fluctuations sampled in molecular dynamics simulations. We explore and validate the model on a large set of single and multicomponent bilayers of different lipid compositions and sizes, simulated at different temperatures. The calculated bilayer compressibility moduli agree with values estimated previously from experiments and those obtained from a standard computational method based on a series of constrained tension simulations. We further validate our framework in a comparison with an existing polymer brush model and confirm the polymer brush model’s predicted linear relationship with proportionality coefficient of 24, using elastic parameters calculated from the simulation trajectories. The robustness of the results that emerge from the method allows us to revisit the origins of the bilayer mechanical (compressible) thickness and in particular its dependence on acyl-chain unsaturation and the presence of cholesterol.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2016

Reducing Disk Storage of Full-3D Seismic Waveform Tomography (F3DT) Through Lossy Online Compression

Full-3D seismic waveform tomography (F3DT) is the latest seismic tomography technique that can assimilate broadband, multi-component seismic waveform observations into high-resolution 3D subsurface seismic structure models. The main drawback in the current F3DT implementation, in particular the scattering-integral implementation (F3DT-SI), is the high disk storage cost and the associated I/O overhead of archiving the 4D space-time wavefields of the receiver- or source-side strain tensors. The strain tensor fields are needed for computing the data sensitivity kernels, which are used for constructing the Jacobian matrix in the Gauss-Newton optimization algorithm. In this study, we have successfully integrated a lossy compression algorithm into our F3DT SI workflow to significantly reduce the disk space for storing the strain tensor fields. The compressor supports a user-specified tolerance for bounding the error, and can be integrated into our finite-difference wave-propagation simulation code used for computing the strain fields. The decompressor can be integrated into the kernel calculation code that reads the strain fields from the disk and compute the data sensitivity kernels. During the wave-propagation simulations, we compress the strain fields before writing them to the disk. To compute the data sensitivity kernels, we read the compressed strain fields from the disk and decompress them before using them in kernel calculations. Experiments using a realistic dataset in our California statewide F3DT project have shown that we can reduce the strain-field disk storage by at least an order of magnitude with acceptable loss, and also improve the overall I/O performance of the entire F3DT-SI workflow significantly. The integration of the lossy online compressor may potentially open up the possibilities of the wide adoption of F3DT-SI in routine seismic tomography practices in the near future.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2018

Autoignition of trans -decalin, a diesel surrogate compound: Rapid compression machine experiments and chemical kinetic modeling

Decahydronaphthalene (decalin), with both cis and trans isomers, is a bicyclic alkane that is found in aviation fuels, diesel fuels, and alternative fuels from tar sands and oil shales. Between the two decalin isomers, trans-decalin has a lower cetane number, is energetically more stable, and has a lower boiling point. Moreover, trans-decalin has often been chosen as a surrogate component to represent two-ring naphthenes in transportation fuels. Recognizing the importance of understanding the chemical kinetics of trans-decalin in the development of surrogate models, an experimental and modeling study has been conducted. Experimentally, the autoignition characteristics of trans-decalin were investigated using a rapid compression machine (RCM) by using trans-decalin/O2/N2 mixtures at compressed pressures of PC=10–25 bar, low-to-intermediate compressed temperatures of TC=620–895 K, and varying equivalence ratios of Φ=0.5, 1.0, and 2.0. These new experimental data demonstrate the effects of pressure, fuel loading, and oxygen concentration on autoignition of trans-decalin. The current RCM data of trans-decalin at lower temperatures were also found to complement well with the literature shock tube data of decalin (mixture of cis+trans) at higher temperatures. Furthermore, a chemical kinetic model for the oxidation of trans-decalin has been developed with new reaction rates and pathways, including, for the first time, a fully-detailed representation of low-temperature chemical kinetics for transdecalin. This model shows good agreement with the overall ignition delay results of the current RCM experiments and the literature shock tube studies. Chemical kinetic analyses of the developed model were further conducted to help identify the fuel decomposition pathways and the reactions controlling the autoignition at varying conditions.

10 SYNTHETIC FUELS↗

DOE OSTI2022

Capture of novel sp 3 hybridized Z-BN by compressing boron nitride nanotubes with small diameter

Experimental synthesis of new sp 3 hybridized carbon/boron nitride structures remains challenging despite that numerous sp 3 structures have been proposed in theory. Here, we showed that compressed multi-walled boron nitride nanotubes (MWBNNTs) and boron nitride peapods (C 60 @BNNTs) with small diameters could transform into a new sp 3 hybridized boron nitride allotrope (Z-BN). This strategy is considered from the topological transition point of view in boron nitride nanotubes upon compression. Due to the increased curvature in compressed small-diameter MWBNNTs, the uncommon 4- and 8-membered rings in Z-BN could be more favorably formed. And the irreversible tube collapse is proved to be a critical factor for the capture of the formed Z-BN, because of the competition between the resilience of tube before collapse and the stress limitation for the lattice stabilization of Z-BN upon decompression. In this case, Z-BN starts to form above 19.0 GPa, which is fully reversible below 45 GPa and finally becomes quenchable at 93.5 GPa. Consequently, this collapse-induced capture of the high-pressure phase could also be extended to other tubular materials for quenching novel sp 3 structures.

36 MATERIALS SCIENCE↗

DOE OSTI2018

Thermodynamic performance and cost optimization of a novel hybrid thermal-compressed air energy storage system design

Compressed Air Energy Storage (CAES) can potentially allow renewable energy sources to meet electricity demands as reliably as coal-fired power plants. However, conventional CAES systems rely on the combustion of natural gas, require large storage volumes, and operate at high pressures, which possess inherent problems such as high costs, strict geological locations, and the production of greenhouse gas emissions. A novel and patented hybrid thermal-compressed air energy storage (HT-CAES) design is presented which allows a portion of the available energy, from the grid or renewable sources, to operate a compressor and the remainder to be converted and stored in the form of heat, through joule heating in a sensible thermal storage medium. The HT-CAES design incudes a turbocharger unit that provides supplementary mass flow rate alongside the air storage. The hybrid design and the addition of a turbocharger have the beneficial effect of mitigating the shortcomings of conventional CAES systems and its derivatives by eliminating combustion emissions and reducing storage volumes, operating pressures, and costs. Storage efficiency and cost are the two key factors, which upon integration with renewable energies would allow the sources to operate as independent forms of sustainable energy. The potential of the HT-CAES design is illustrated through a thermodynamic optimization study, which outlines key variables that have a major impact on the performance and economics of the storage system. The optimization analysis quantifies the required distribution of energy between thermal and compressed air energy storage, for maximum efficiency, and for minimum cost. This study provides a roundtrip energy and exergy efficiency map of the storage system and illustrates a trade off that exists between its capital cost and performance.

25 ENERGY STORAGE↗

DOE OSTI2022

Combustion dynamics of crude and upgraded Thermal DeOxygenation oils in a compression ignition engine

Thermal DeOxygenation (TDO) is a robust thermochemical conversion scheme to produce hydrocarbons from biomass feedstock. The process targets the carbohydrate fraction of biomass to yield a broad mixture of primarily aromatic hydrocarbons within a boiling point range of 348–798 K with low oxygen content (<4 wt%). The resulting materials are amenable to traditional hydrotreating and distillation processes with approximately 70% by mass in the distillate fuel range. The simple conversion scheme combined with commercially viable upgrading routes makes TDO oils a feasible alternative fuel for transportation applications. Here this paper is the first systematic treatment of fit-for-purpose testing of TDO oils in a compression ignition engine. Blends of partially upgraded TDO oils (e.g. whole oil, distilled oil, hydrotreated oil and hydrotreated-distilled oil) are prepared with certified ultra-low sulfur diesel at 5%, 10%, 15% and 20% by volume. The resulting fuels are analyzed for fuel characteristics and combustion dynamics in an instrumented single-cylinder compression ignition engine. All fuel blends at 10% blend level by volume exhibited adequate engine performance. Hydrotreating of the oils is a necessary step to meet cetane number and EPA soot emissions requirements at 20% blend volume. Combined distillation and hydrotreatment meet or exceed all fuel specification and engine performance benchmarks. Partially upgraded TDO oils are suitable fuel options for compression ignition applications.

42 ENGINEERING↗

DOE OSTI2026

An analysis of alcohol compression ignition strategies on different engine platforms

Strategies that enable ignition of the low cetane alcohol fuels (methanol and ethanol) include intake heating (direct or via hot residual trapping), increased compression ratio (CR), multi-injection strategies, ignition enhancers, and active prechambers. The literature consistently shows alcohol compression ignition (CI) is more efficient than diesel at high loads irrespective of ignition strategy but shows contradictory results at low loads. To understand this, this work combines data from four different engine platforms ranging from a displaced volume of 0.4 to 2.1 L/cyl. using different alcohol CI ignition strategies. Using modeling to support experimental analysis, alcohol-specific efficiency penalties are highlighted to inform alcohol combustion system development. The alcohols’ lower energy density means more fuel must be injected near top dead center to achieve a given load, increasing the relative penalty of sensible/latent heating. Pilot injections help reduce this penalty. Intake heating incurs both a heat transfer penalty and a thermodynamic penalty due to a reduced ratio of specific heats during compression. A strategy that combines elevated CR, large pilot injections, and ignition enhancers can avoid this penalty. Alternatively, a prechamber strategy can achieve this without changing the CR, nor needing ignition enhancers in the fuel. At low loads, diesel burns in a partially premixed mode characterized by high burn rates and low heat losses, reducing alcohol-specific combustion benefits. Alcohol CI on the light-duty platform showed abnormally high heat transfer over diesel compared to other platforms, favoring a more premixed/partially premixed strategy despite a combustion efficiency tradeoff.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2025

Numerical simulation of compressible fluid-dynamics in the chamber of inertial fusion energy systems

Here, this paper aims to establish new and innovative modeling capabilities for analyzing chambers in Inertial Fusion Energy (IFE) systems. IFE is emerging as a promising method to achieve fusion power production, but several challenges must be overcome to develop an IFE pilot plant or deploy commercial IFE systems. These challenges are both theoretical and technical, encompassing a deeper understanding of the underlying physical phenomena and the development of new technologies and materials. One of the needs is to develop mathematical models to describe IFE systems and numerical tools to simulate them. This paper contributes to this endeavor by presenting a new OpenFOAM solver for IFE systems, focusing on gas dynamics in their chambers. The analysis and development of chamber designs will play a significant role in the transition from single-shot experiments to high-repetition rates, as there is a need to protect the chamber walls from the intense radiation fields produced by fusion reactions. A promising design option, normally referred to as thick wall chamber design, consists in using lithium or molten salt jet arrays within the chamber. A critical phenomenon is the venting of high-pressure gases from the center to the external part of the chamber, passing through the blanket array. This process involves the propagation and attenuation of strong pressure waves, requiring suitable modeling approaches for compressible fluid-dynamics. The solver proposed in this work implements a multi-material hydrodynamics model tailored to accurately describe the non-linear propagation of pressure waves while avoiding numerical oscillation issues typical of high-velocity compressible simulation. This solver is verified against numerical test cases, validated against experimental data, and applied to the analysis of the High-Yield Lithium-Injection Fusion-Energy (HYLIFE-I) concept. The relevance of this paper is threefold. Firstly, it contributes to developing and testing modeling approaches for compressible fluid-dynamics phenomena, with specific focus on the new and unexplored topic of IFE thick-liquid-wall blanket modeling. Secondly, it marks one of the first applications of the OpenFOAM library in the research field of IFE systems. Finally, the investigated problem is of practical interest for IFE developers, as it provides useful indications about relevant phenomena in pressure wave propagation in the chamber of these systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2014

A new limiting procedure for discontinuous Galerkin methods applied to compressible multiphase flows with shocks and interfaces

Although the Discontinuous Galerkin (dg) method has seen widespread use for compressible flow problems in a single fluid with constant material properties, it has yet to be implemented in a consistent fashion for compressible multiphase flows with shocks and interfaces. Specifically, it is challenging to design a scheme that meets the following requirements: conservation, high-order accuracy in smooth regions and non-oscillatory behavior at discontinuities (in particular, material interfaces). Following the interface-capturing approach of Abgrall [1], we model flows of multiple fluid components or phases using a single equation of state with variable material properties; discontinuities in these properties correspond to interfaces. Here, to represent compressible phenomena in solids, liquids, and gases, we present our analysis for equations of state belonging to the Mie–Grüneisen family. Within the dg framework, we propose a conservative, high-order accurate, and non-oscillatory limiting procedure, verified with simple multifluid and multiphase problems. We show analytically that two key elements are required to prevent spurious pressure oscillations at interfaces and maintain conservation: (i) the transport equation(s) describing the material properties must be solved in a non-conservative weak form, and (ii) the suitable variables must be limited (density, momentum, pressure, and appropriate properties entering the equation of state), coupled with a consistent reconstruction of the energy. Further, we introduce a physics-based discontinuity sensor to apply limiting in a solution-adaptive fashion. We verify this approach with one- and two-dimensional problems with shocks and interfaces, including high pressure and density ratios, for fluids obeying different equations of state to illustrate the robustness and versatility of the method. Lastly, the algorithm is implemented on parallel graphics processing units (gpu) to achieve high speedup.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗