Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “INTERFACES”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Screw dislocation impingement and slip transfer at fcc-bcc semicoherent interfaces

Understanding of dislocation-interface interactions at the atomistic level could provide useful insights into the mechanisms behind slip transfer, a key information for higher length-scale formulations. Here, a systematic study of screw dislocation impingement on both Fe/Al and Fe/Cu interfaces with Kurdjumov-Sachs orientation relationship has been carried out using atomistic simulations. For the “strong” Fe/Al interface, there is appreciable dislocation delocalization, in which the change in dislocation content provided by Nye tensor analysis. Interestingly, the delocalization is significant enough that large simulation cells must be employed in order to avoid boundary effects. When the impinging dislocations are far away from the boundaries, the large accumulation of shear strains nearby the impinging dislocations at the interface causes dislocation nucleation/emission at these places. While for the “weak” Fe/Cu interface, the dislocation delocalization at the interface is notably long-range, extending more than 20 nm away from where the initial dislocations entered the interface.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Regenerative Solid Interfaces Enhance High-Performance All- Solid-State Lithium Batteries

All-solid-state Lithium batteries (ASSLBs) have gained widespread attention in recent years. However, their performance is still largely limited by the poor transport properties and stability of the solid-solid interfaces (SSIs). In this study, we report a new strategy for creating a reversible interface with good conductivity and self-adaptive mechanical properties for high-performance all-solid-state lithium batteries (ASSLBs). The interface is formed in situ from nanosized lithium iodide (LiI), a component of the solid-state electrolyte (SSE), which flows plastically along the SSE interfaces under high pressure due to its high ductility. Moreover, LiI segregates and enriches along the Li/SSE interfaces, reversibly enhancing Li plating/stripping due to its lithophilicity and high ionic conductivity. This dynamic LiI interface enables stable operation of metallic Li anode (>1000 h) at high current densities and elevated temperatures, and long-term cycling of all-solid-state Li-sulfur batteries (>250 cycles) with high sulfur utilization rate (>1400 mAh g-1) and areal capacity (>2 mAh cm-2). This study reveals a unique role of LiI in building robust SSIs and provides new insights into the design of novel SSEs and interfaces for high-performance ASSLBs.

Yu, Zhaoxin↗
DOE OSTI2025

Interface Evolutionand Long-Term Performance of NegativeCarbon Fiber Structural Electrodes

Abstract Laminated structural batteries present a transformative solution to reducing weight constraints in electric vehicles. These structural batteries are based on a multifunctional material that incorporates an energy storage function within a carbon fiber-reinforced polymer. Despite the potential of this technology, the intricate morphology of fiber–matrix or electrode–electrolyte interfaces and the impact of long-term cycling at low current rates (C-rates) on these interfaces remain insufficiently understood. This study addresses these critical knowledge gaps by examining the influence of matrix composition on the long-term electrochemical performance of structural battery electrodes and exploring advanced techniques to investigate carbon fiber–matrix interfaces. Localized imaging and X-ray scattering techniques were used to characterize morphological changes at the electrode–electrolyte interfaces by analyzing negative structural electrodes. The findings revealed that the matrix composition influences long-term electrochemical behavior and fiber–matrix interface formation. While the intrinsic properties of carbon fibers largely remain unaffected by long-term cycling, cycling promotes debonding at fiber–matrix interfaces. Nonetheless, residual regions of adhesion persist, underscoring the potential for preserving multifunctionality even under prolonged cycling conditions. These insights advance the understanding of interface dynamics, which is critical for optimizing structural battery technologies.

Chemistry↗
DOE OSTI2021

Reversible assembly of silica nanoparticles at water–hydrocarbon interfaces controlled by SDS surfactant

Achieving reversible and tunable assembly of silica nanoparticles at liquid–liquid interfaces is vital for a wide range of scientific and technological applications including sustainable subsurface energy applications, catalysis, drug delivery and material synthesis. In this study, we report the mechanisms controlling the assembly of silica nanoparticles (dia. 50 nm and 100 nm) at water–heptane and water–toluene interfaces using sodium dodecyl sulfate (SDS) surfactant with concentrations ranging from 0.001–0.1 wt% using operando ultrasmall/small-angle X-ray scattering, cryogenic scanning electron microscopy imaging and classical molecular dynamics simulations. The results show that the assembly of silica nanoparticles at water–hydrocarbon interfaces can be tuned by controlling the concentrations of SDS. Silica nanoparticles are found to: (a) dominate the interfaces in the absence of interfacial SDS molecules, (b) coexist with SDS at the interfaces at low surfactant concentration of 0.001 wt% and (c) migrate toward the aqueous phase at a high SDS concentration of 0.1 wt%. Energetic analyses suggest that the van der Waals and electrostatic interactions between silica nanoparticles and SDS surfactants increase with SDS concentration. However, the favorable van der Waals and electrostatic interactions between the silica nanoparticles and toluene or heptane decrease with increasing SDS concentration. As a result, the silica nanoparticles migrate away from the water–hydrocarbon interface and towards bulk water at higher SDS concentrations. These calibrated investigations reveal the mechanistic basis for tuning silica nanoparticle assembly at complex interfaces.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2022

Plasma parameters and the reduction potential at a plasma-liquid interface

We report nonthermal plasmas in contact with liquids have been shown to generate a variety of reactive species capable of initiating reduction–oxidation (redox) reactions at the electrochemically active plasma–liquid interface. In conventional electrochemical cells, selective redox chemistry is achieved by controlling the reduction potential at the solid electrode–electrolyte interface by applying a bias via an external circuit. In the case of plasma–liquid systems, an analogous means of tuning the reduction potential near the interface has not clearly been identified. When treated as a floating surface, the liquid is expected to adopt a net negative charge to balance the flux of hot electrons and relatively cold positive ions. The reduction potential near the plasma–liquid interface is hypothesized to be proportional to the floating potential, which can be approximated using an analytical model provided the plasma parameters are known. Herein, we present a framework for correlating the electron density and electron temperature of a noble gas plasma jet to the reduction potential near the plasma–liquid interface. The plasma parameters were acquired for an argon atmospheric plasma jet in contact with an aqueous solution by means of laser Thomson scattering. The reduction potential was determined using identical reference electrodes to measure the potential difference between the plasma–liquid interface and bulk solution. Interestingly, the measured reduction potentials near the plasma–liquid interface were found to be in good agreement with the model-predicted values determined using the plasma parameters obtained from the Thomson scattering experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Assembly of polyelectrolyte star block copolymers at the oil–water interface

To understand and resolve adsorption, reconfiguration, and equilibrium conformations of charged star copolymers, we carried out an integrated experimental and coarse-grained molecular dynamics simulation study of the assembly process at the oil–water interface. This is important to guide development of novel surfactants or amphiphiles for chemical transformations and separations. The star block copolymer consisted of arms that are comprised of hydrophilic–hydrophobic block copolymers that are covalently tethered via the hydrophobic blocks to one point. The hydrophobic core represents polystyrene (PS) chains, while the hydrophilic corona represents quaternized poly(2-vinylpyridine) (P2VP) chains. The P2VP is modeled to become protonated when in contact with an acidic aqueous phase, thereby massively increasing the hydrophilicity of this block, and changing the nature of the star at the oil–water interface. This results in a configurational change whereby the chains comprising the hydrophilic corona are significantly stretched into the aqueous phase, while the hydrophobic core remains solubilized in the oil phase. In the simulations, we followed the kinetics of the anchoring and assembly of the star block copolymer at the interface, monitoring the lateral assembly, and the subsequent reconfiguration of the star via changes in the interfacial tension that varies as the degree-of-protonation increases. At low fractions of protonation, the arm cannot fully partition into the aqueous side of the interface and instead interacts with other arms in the oil phase forming a network near the interface. These insights were used to interpret the non-monotonic dependence of pH with the asymptotic interfacial tension from pendant drop tensiometry experiments and spectral signatures of aromatic stretches seen in vibrational sum frequency generation (SFG) spectroscopy. We describe the relationship of interfacial tension to the star assembly via the Frumkin isotherm, which phenomenologically describes anti-cooperativity in adsorbing stars to the interface due to crowding. Although our model explicitly considers long-range electrostatics, the contribution of electrostatics to interfacial tension is small and brought about by strong counterion condensation at the interface. Finally, these results provide key insights into resolving the adsorption, reconfiguration, and equilibrium conformations of charged star block copolymers as surfactants.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Design optimization for Richtmyer–Meshkov instability suppression at shock-compressed material interfaces

We report the Richtmyer–Meshkov instability (RMI) is a phenomenon that occurs at the interface of two substances of different densities due to an impulsive acceleration, such as a shock wave passing through this interface. Under these conditions, the instability can be seen as interface perturbations begin to grow into narrow jets or spikes of one substance that propagate into the other. In some cases, this interface may involve an elastic–plastic material, which can play a significant role in the development and behavior of the RMI. The ability to effectively control RMI jetting and spike growth is one major limiting factor in technological challenges, such as inertial confinement fusion, that involve using high-pressure shock waves to implode a fuel target. The propagation of RMI growth can lead to increased asymmetry in this implosion process and significantly reduce the obtained energy yield. We use hydrodynamics simulations of impactor shock-compression experiments and methods based in design optimization to suppress RMI spike growth by altering the geometry and other properties of a shock-compressed elastic–plastic material target that shares an interface with atmospheric air. These hydrodynamics simulations use an arbitrary Lagrangian–Eulerian method with a high-order finite element approach. Our results demonstrate that RMI suppression can be achieved by intentionally creating a separate upstream interface instability to counteract the growth of long narrow RMI spikes at an interface with initial perturbations.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Influence of nanoscale interfaces on the dynamic deformation and spall failure of Cu–Fe alloy microstructures

Additively manufactured immiscible Cu–Fe alloys can exhibit a distribution of nanoscale interfaces due to the distribution of nanoscale clusters in the equilibrium and metastable phases. Molecular dynamics simulations investigate the role of such interfaces on the phase stability and transformation behavior during shock compression, as well as the mechanisms of damage nucleation during spall failure. The model multiphase Cu–Fe systems studied here comprise a distribution of Fe clusters in an FCC Cu matrix, as well as Cu clusters in a BCC Fe matrix. The length scales of the nanoscale clusters determine the energetics of the interfaces that can result in FCC/BCC (equilibrium) or FCC/FCC (metastable) interfaces in the Cu system matrix, and BCC/FCC (equilibrium) or BCC/BCC (metastable) interfaces in the Fe matrix. The MD simulations demonstrate that nanoscale metastable interface microstructures can induce plastic deformation and also stabilize phases or suppress phase transformations in metastable phase clusters. In contrast, equilibrium interfaces can influence phase transformation thresholds and serve as additional void nucleation sites during failure.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Temporal evolution of scalar modes in Richtmyer–Meshkov instability of inclined interface using high-speed PIV and PLIF measurements at 60 kHz

Abstract The current work presents simultaneous, high-speed measurements at 60,000 fields per second of velocity and mole fraction using particle image velocimetry (PIV) and planar laser induced acetone-fluorescence in a Richtmyer–Meshkov instability of an inclined interface (Atwood number, At = 0.22). Specifically, around 2 ms of temporal evolution of the vortex structures and their associated scalar modes immediately following the interface-reshock interaction is presented. Two initial interface conditions are discussed—(a) a sharp, inclined ‘single mode’ interface and (b) a ‘multi-mode’ interface where small perturbations are imposed on the single mode case. A 2D wavelet decomposition of the scalar flow field shows a highly intermittent distribution of small-scale variance throughout the interface even at late times. These are correlated strongly with the vortex structures and local turbulence intensity, where each small-scale scalar mode is sandwiched between two co-rotating vortex structures. This indicates that the interstitial regions between the vortices are significant hotspots of entrainment, which is then dispersed by the induced, counter-flow velocity fields. The multimode case demonstrates similar organization at large scales, while the scalar field is much more homogeneous at smaller scales. These observations highlight the importance of capturing the early time vortex evolution to accurately estimate any late time intermittency, especially where deposition of intense vorticity on sharp interfaces is present.

Engineering↗
DOE OSTI2021

Enhanced Interface-Driven Perpendicular Magnetic Anisotropy by Symmetry Control in Oxide Superlattices

Perpendicular magnetic anisotropy (PMA) has recently been shown to emerge at interfaces of 3d and 5d transition-metal oxides (TMOs). However, strategies to systematically stabilize such interface-driven PMA still remains elusive, hindering further applications of this design approach. Here, tuning crystal symmetry is shown to be an effective means to engineer this interfacial phenomenon. The evolution of PMA strength as a function of ferromagnetic oxide thickness quantitatively reveals the competition between volume- and interface-specific contributions that determine the magnetic anisotropy. By applying different degrees of epitaxial strain, the relative contributions to PMA are modulated, clearly revealing their correlations with crystal symmetries. To be more specific, the volume anisotropy energy is found to be correlated with the tetragonal distortion of the ferromagnetic layer, while the interface anisotropy energy is mainly modulated by the octahedral tilting at the interface. With these insights, superlattices with enhanced interface-driven PMA and higher Curie temperature are realized. These findings reveal a route to engineering interface-driven PMA and associated magnetic phenomena in TMO heterostructures for future spintronic applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

MILK : a Python scripting interface to MAUD for automation of Rietveld analysis

Modern diffraction experiments ( e.g. in situ parametric studies) present scientists with many diffraction patterns to analyze. Interactive analyses via graphical user interfaces tend to slow down obtaining quantitative results such as lattice parameters and phase fractions. Furthermore, Rietveld refinement strategies ( i.e. the parameter turn-on-off sequences) tend to be instrument specific or even specific to a given dataset, such that selection of strategies can become a bottleneck for efficient data analysis. Managing multi-histogram datasets such as from multi-bank neutron diffractometers or caked 2D synchrotron data presents additional challenges due to the large number of histogram-specific parameters. To overcome these challenges in the Rietveld software Material Analysis Using Diffraction ( MAUD ), the MAUD Interface Language Kit ( MILK ) is developed along with an updated text batch interface for MAUD . The open-source software MILK is computer-platform independent and is packaged as a Python library that interfaces with MAUD . Using MILK , model selection ( e.g. various texture or peak-broadening models), Rietveld parameter manipulation and distributed parallel batch computing can be performed through a high-level Python interface. A high-level interface enables analysis workflows to be easily programmed, shared and applied to large datasets, and external tools to be integrated with MAUD . Through modification to the MAUD batch interface, plot and data exports have been improved. The resulting hierarchical folders from Rietveld refinements with MILK are compatible with Cinema: Debye–Scherrer , a tool for visualizing and inspecting the results of multi-parameter analyses of large quantities of diffraction data. In this manuscript, the combined Python scripting and visualization capability of MILK is demonstrated with a quantitative texture and phase analysis of data collected at the HIPPO neutron diffractometer.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Signatures of a spin-active interface and a locally enhanced Zeeman field in a superconductor-chiral material heterostructure

A localized Zeeman field, intensified at heterostructure interfaces, could play a crucial role in a broad area including spintronics and unconventional superconductors. Conventionally, the generation of a local Zeeman field is achieved through magnetic exchange coupling with a magnetic material. However, magnetic elements often introduce defects, which could weaken or destroy superconductivity. Alternatively, the coupling between a superconductor with strong spin-orbit coupling and a nonmagnetic chiral material could serve as a promising approach to generate a spin-active interface. Here, we leverage an interface superconductor, namely, induced superconductivity in noble metal surface states, to probe the spin-active interface. Our results unveil an enhanced interface Zeeman field, which selectively closes the surface superconducting gap while preserving the bulk superconducting pairing. The chiral material, i.e., trigonal tellurium, also induces Andreev bound states (ABS) exhibiting spin polarization. The field dependence of ABS manifests a substantially enhanced interface Landég-factor (g eff ~ 12), thereby corroborating the enhanced interface Zeeman energy.

Science & Technology - Other Topics↗
DOE OSTI2023

Degradation Effects at the Porous Transport Layer/Catalyst Layer Interface in Polymer Electrolyte Membrane Water Electrolyzer

The porous transport layer (PTL)/catalyst layer (CL) interface plays a crucial role in the achievement of high performance and efficiency in polymer electrolyte membrane water electrolyzers (PEMWEs). This study investigated the effects of the PTL/CL interface on the degradation of membrane electrode assemblies (MEAs) during a 4000 h test, comparing the MEAs assembled with uncoated and Ir-coated Ti PTLs. Our results show that compared to an uncoated PTL/CL interface, an optimized interface formed when using a platinum group metal (PGM) coating, i.e., an iridium layer at the PTL/CL interface, and reduced the degradation of the MEA. The agglomeration and formation of voids and cracks could be found for both MEAs after the long-term test, but the incorporation of an Ir coating on the PTL did not affect the morphology change or oxidation of IrO x in the catalyst layer. In addition, our studies suggest that the ionomer loss and restructuring of the anodic MEA can also be reduced by Ir coating of the PTL/CL interface. Optimization of the PTL/CL interface improves the performance and durability of a PEMWE.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2020

An Integrated Energy Systems Prototype Human-System Interface for a Steam Extraction Loop System to Support Joint Electricity-Hydrogen Flexible Operations

Due to increasing economic competition from renewables and combined-cycle natural gas plants, nuclear power plants are looking toward flexible operations to enhance their cost competitiveness. The Integrated Energy Systems project under the Light Water Reactor Sustainability program of the U.S. Department of Energy focuses on joint electricity-hydrogen flexible operations. Joint electricity-hydrogen flexible operations entail the nuclear power plant diverting thermal energy via main steam to a hydrogen production plant located nearby. The steam serves to enhance the efficiency of the hydrogen production. Furthermore, high temperature electrolysis requires a large amount of electricity, which the plant can also provide. The plant provides steam and electricity to the hydrogen plant throughout the day, but during peak demand hours the nuclear power plant returns to solely providing electricity to meet the high demand. Through this flexible concept of operations, the plant can optimize the thermal energy it produces without having to maneuver the power of the plant. This report documents the human factors process to design and develop a prototype human-system interface for the steam extraction loop that serves as the conduit between the nuclear power plant and the adjacent hydrogen plant. The design process followed the human factors guidelines set by NUREG-0711, Human Factors Engineering Program Review Model (O’Hara, Higgins, & Fleger, 2012), and expanded upon by the Guideline for Operational Nuclear Usability and Knowledge Elicitation (GONUKE; Boring, Ulrich, Joe, & Lew, 2015; Boring, Lew, & Ulrich, 2016). The design process entailed operator interviews to determine the concept of operations for the steam extraction loop, a review and adaptation of digital interface design concepts developed by the team in prior projects, an iterative design process, and reviews conducted by both operators and human factors experts. Several versions of the prototype human-system interface were developed. Operators were interviewed to determine what design features they found useful and would like to see in the interface. The design underwent a review by human factors experts against NUREG-0700, Human Interface Design Review Guidelines (U.S. Nuclear Regulatory Commission, 2019), to ensure compliance with the latest human factors standards for digital interfaces in nuclear applications. The design was then prototyped as a functional windows-based application integrated with the Generic Pressurized Water Reactor simulator modified to include the steam extraction loop. The simulation is supported by the Human Systems Simulation Laboratory at Idaho National Laboratory, which supports operator-in-the-loop testing. This is an ongoing project and the next phase of the project entails performing an operator-in-the-loop usability study to evaluate the interface and examine the proposed concept of operations to extraction steam from the nuclear power plant for delivery to the coupled hydrogen production plant.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2026

Evaluation of RMI Ejecta Sourcing from Reconstructed Interfaces in 2D FLAG Flyer-Plate Simulations

This study investigates the performance of RMI ejecta sourcing from a reconstructed interface in two-dimensional FLAG ejecta flyer-plate simulations. Enabling the RMI ejecta model to operate on reconstructed interfaces has been a longstanding challenge in FLAG, limiting its applicability to experiments in which interface reconstruction is not required. Addressing this limitation is therefore an important step toward expanding the model’s predictive capability. The first phase of this effort focused on identifying a flyer-plate/target configuration that produces large ejecta clouds from two distinct shocks using the standard RMI ejecta model. A series of simulations examined the effects of flyer composition, layer thickness, and mesh resolution on shock structure and ejecta formation. VISAR traces were used to evaluate shock timing and identify configurations that produced a clear double-shock response. A copper middle-layer flyer with a doubled tin layer yielded the most favorable shock characteristics and was selected as the baseline configuration. After confirming this configuration performed as expected when running a reference FLAG version, it became the basis for three simulations: (1) a standard two-dimensional ejecta flyer-plate calculation using a reference FLAG version, (2) the same calculation with the reconstructed interface option enabled, and (3) a calculation employing ALE strategies to force the ejecta model to operate on the reconstructed interface. While the baseline simulation produced the expected ejecta response, the reconstructed interface approach did not result in measurable ejecta, indicating that this capability remains under development. These results highlight both the progress made and the remaining challenges in extending the RMI ejecta model to reconstructed interfaces.

36 MATERIALS SCIENCE↗
DOE OSTI2023

In situ mechanical characterization of a single cell-cell adhesion interface under large strain

A method of measuring a stress-strain curve in a cell-cell adhesion interface, the method including: providing a structure including a first movable island supported by a first beam, a second movable island supported by a second beam, and a gap therebetween connected by a pair of cells forming a junction, the pair of cells comprising a cell-cell adhesion interface having an initial length defined by a distance between nuclei of the pair of cells; moving the second movable island with a defined displacement; determining a displacement of the first movable island based on moving the second movable island; calculating a difference between the displacement of the first movable island and the defined displacement of the second movable island based on moving the second movable island; determining an applied strain in the cell-cell adhesion interface between the pair of cells based on the difference divided by the initial length of the cell-cell adhesion interface; calculating a force between the cell-cell adhesion interface of the pair of cells based on the displacement of the first movable island; calculating a stress in the cell-cell adhesion interface between the pair of cells based on the force; and determining the stress-strain curve of the cell-cell adhesion interface between the pair of cells by plotting the calculated stress against the applied strain.

Yang, Ruiguo↗
DOE OSTI2022

An investigation of the multi-mode Richtmyer-Meshkov instability at a gas/HE interface using Pagosa

Here, in this work, we present a hydrocode Pagosa and explore the Richtmyer-Meshkov Instability (RMI) at an air/high explosive (HE) interface for the first time that is important but has not received much attention yet in the high explosive safety field. Thus, the presented Pagosa can be expected to predict the whole deflagration-to-detonation transition (DDT) process in future. In Pagosa, spatial discretization is implemented on cubic staggered grids by computing different variables at the vertex and the cell center, respectively, a special operator-splitting technique is employed to reduce the computational cost, and an artificial viscosity is added to handle the discontinuous shock waves in our simulations. To quantitatively evaluate the capability of Pagosa to solve these kinds of instabilities, the single mode Rayleigh-Taylor instability (RTI) and the multimode RMI at an air/SF 6 interface are first performed, respectively. The Pagosa results are compared with the related numerical solutions in the existing references and the experimental result. Moreover, a theoretical derivation of growth of RTI is also provided based on our numerical method. Subsequently, we explore the multi-mode RMI at an air/HE interface as well as the effects of several factors using Pagosa. Numerical results show that Pagosa is a powerful toolset to generate the right structures and the amplitude of RTI and RMI at an air/SF 6 interface. The solid HE can be penetrated by a strong shock wave and forms RMI deformations. The RMI at an air/HE interface behaves very different than at an air/SF 6 interface, periodic, decreased oscillation is observed due to material character, and is very sensitive to the initial simulation settings, that is, a tiny change in physical quantities will lead to a remarkable RMI structure, which is also observed in a shock bubble interaction. The findings in this work are significant, and will present a new insight for the high explosive field.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Assessment of diffuse-interface methods for compressible multiphase fluid flows and elastic-plastic deformation in solids

This work describes three diffuse-interface methods for the simulation of immiscible, compressible multiphase fluid flows and elastic-plastic deformation in solids. The first method is the localized-artificial-diffusivity approach of Cook, Subramaniam et al., and Adler and Lele, in which artificial diffusion terms are added to the individual phase mass fraction transport equations and are coupled with the other conservation equations. The second method is the gradient-form approach that is based on the quasi-conservative method of Shukla et al., in which the diffusion and sharpening terms (together called regularization terms) are added to the individual phase volume fraction transport equations and are coupled with the other conservation equations. The third approach is the divergence-form approach that is based on the fully conservative method of Jain et al., in which the regularization terms are added to the individual phase volume fraction transport equations and are coupled with the other conservation equations. In the present study, all three diffuse-interface methods are used in conjunction with a four-equation, multicomponent mixture model, in which pressure and temperature equilibria are assumed among the various phases. The primary objective of this work is to compare these three methods in terms of their ability to: maintain constant interface thickness throughout the simulation; conserve mass, momentum, and energy; and maintain accurate interface shape for long-time integration. The second objective of this work is to consistently extend these methods to model interfaces between solid materials with strength. To assess and compare the methods, they are used to simulate a wide variety of problems, including (1) advection of an air bubble in water, (2) shock interaction with a helium bubble in air, (3) shock interaction and the collapse of an air bubble in water, and (4) Richtmyer–Meshkov instability of a copper–aluminum interface. The current work focuses on comparing these methods in the limit of relatively coarse grid resolution, which illustrates the true performance of these methods. In conclusion, this is because it is rarely practical to use hundreds of grid points to resolve a single bubble or drop in large-scale simulations of engineering interest.

97 MATHEMATICS AND COMPUTING↗