SEARCH · Engineering Papers

Results for “INTERFACE”

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 721 records · Page 40

All sources NASA NTRS DOE OSTI

DOE OSTI2013

An interface tracking model for droplet electrocoalescence.

This report describes an Early Career Laboratory Directed Research and Development (LDRD) project to develop an interface tracking model for droplet electrocoalescence. Many fluid-based technologies rely on electrical fields to control the motion of droplets, e.g. microfluidic devices for high-speed droplet sorting, solution separation for chemical detectors, and purification of biodiesel fuel. Precise control over droplets is crucial to these applications. However, electric fields can induce complex and unpredictable fluid dynamics. Recent experiments (Ristenpart et al. 2009) have demonstrated that oppositely charged droplets bounce rather than coalesce in the presence of strong electric fields. A transient aqueous bridge forms between approaching drops prior to pinch-off. This observation applies to many types of fluids, but neither theory nor experiments have been able to offer a satisfactory explanation. Analytic hydrodynamic approximations for interfaces become invalid near coalescence, and therefore detailed numerical simulations are necessary. This is a computationally challenging problem that involves tracking a moving interface and solving complex multi-physics and multi-scale dynamics, which are beyond the capabilities of most state-of-the-art simulations. An interface-tracking model for electro-coalescence can provide a new perspective to a variety of applications in which interfacial physics are coupled with electrodynamics, including electro-osmosis, fabrication of microelectronics, fuel atomization, oil dehydration, nuclear waste reprocessing and solution separation for chemical detectors. We present a conformal decomposition finite element (CDFEM) interface-tracking method for the electrohydrodynamics of two-phase flow to demonstrate electro-coalescence. CDFEM is a sharp interface method that decomposes elements along fluid-fluid boundaries and uses a level set function to represent the interface.

Erickson, Lindsay Crowl↗

DOE OSTI2014

Area of Interest 1, CO 2 at the Interface: Nature and Dynamics of the Reservoir/Caprock Contact and Implications for Carbon Storage Performance

We examined the influence of geologic features present at the reservoir/caprock interface on the transmission of supercritical CO 2 into and through caprock. We focused on the case of deformation-band faults in reservoir lithologies that intersect the interface and transition to opening-mode fractures in caprock lithologies. Deformation-band faults are exceeding common in potential CO 2 injection units and our fieldwork in Utah indicates that this sort of transition is common. To quantify the impact of these interface features on flow and transport we first described the sedimentology and permeability characteristics of selected sites along the Navajo Sandstone (reservoir lithology) and Carmel Formation (caprock lithology) interface, and along the Slickrock Member (reservoir lithology) and Earthy Member (caprock lithology) of the Entrada Sandstone interface, and used this information to construct conceptual permeability models for numerical analysis. We then examined the impact of these structures on flow using single-phase and multiphase numerical flow models for these study sites. Key findings include: (1) Deformation-band faults strongly compartmentalize the reservoir and largely block cross-fault flow of supercritical CO 2 . (2) Significant flow of CO 2 through the fractures is possible, however, the magnitude is dependent on the small-scale geometry of the contact between the opening-mode fracture and the deformation band fault. (3) Due to the presence of permeable units in the caprock, caprock units are capable of storing significant volumes of CO 2 , particularly when the fracture network does not extend all the way through the caprock. The large-scale distribution of these deformation-bandfault-to-opening-mode-fractures is related to the curvature of the beds, with greater densities of fractures in high curvature regions. We also examined core and outcrops from the Mount Simon Sandstone and Eau Claire Formation reservoir/caprock interface in order to extend our work to a reservoir/caprock pair this is currently being assessed for long-term carbon storage. These analyses indicate that interface features similar to those observed at the Utah sites 3 were not observed. Although not directly related to our main study topic, one byproduct of our investigation is documentation of exceptionally high degrees of heterogeneity in the pore-size distribution of the Mount Simon Sandstone. This suggests that the unit has a greater-than-normal potential for residual trapping of supercritical CO 2 .

54 ENVIRONMENTAL SCIENCES↗

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↗

DOE OSTI2018

Diffuse-Interface Capturing Methods for Compressible Two-Phase Flows

Simulation of compressible flows became a routine activity with the appearance of shock-/contact-capturing methods. These methods can determine all waves, particularly discontinuous ones. However, additional difficulties may appear in two-phase and multimaterial flows due to the abrupt variation of thermodynamic properties across the interfacial region, with discontinuous thermodynamical representations at the interfaces. To overcome this difficulty, researchers have developed augmented systems of governing equations to extend the capturing strategy. These extended systems, reviewed here, are termed diffuse-interface models, because they are designed to compute flow variables correctly in numerically diffused zones surrounding interfaces. In particular, they facilitate coupling the dynamics on both sides of the (diffuse) interfaces and tend to the proper pure fluid–governing equations far from the interfaces. This strategy has become efficient for contact interfaces separating fluids that are governed by different equations of state, in the presence or absence of capillary effects, and with phase change. As a result, more sophisticated materials than fluids (e.g., elastic–plastic materials) have been considered as well.

42 ENGINEERING↗

DOE OSTI2020

Ion redistributions at interfaces facilitate nucleation and growth of branched Ag 3 PO 4 polypods

Branched structures are of spectacular interest due to their improved light-harvesting, direct carrier transportation pathway, high surface area, and controllable electronic structures. Understanding their growth mechanisms and controlling factors enables the design of materials with improved functions. Interfaces of liquid-solid and air-liquid are known to change the chemical and physical properties of the liquid solution. However, their role in crystal growth, especially of branched structures, is seldom investigated. In this work we synthesized branched silver phosphate polypods in a thin layer of aqueous solution at room temperature taking advantage of the unique properties near the interfaces. The branched crystals nucleate and grow with a liquid thickness of 810 µm or below requiring the presence of both the gas-liquid and liquid-solid interfaces. The hydrophilicity of the solid substrate also facilitates nucleation and growth. Our results indicate that the synergism of long-range electric fields (ion distributions) at the interfaces promotes not only the nucleation and growth of silver phosphate but also the formation of branched structures. This work can be a reference for facilitating crystal nucleation and growth and controlling structures with specific morphologies, such as branched ones, by understanding and mimicking the properties at the interfaces.

36 MATERIALS SCIENCE↗

DOE OSTI2026

All‐Solid‐State Batteries With Mechanically Stable Interfaces Consisting of a Zero‐Strain Cation‐Disordered Rocksalt Cathode

Interface stabilization is critical to the development of working all‐solid‐state batteries. Rigid cathode/solid electrolyte interfaces often disintegrate due to anisotropic volume change of cathode‐active materials, resulting in irreversible capacity loss. Herein, we demonstrate that Li 1.211 Mo 0.467 Cr 0.3 O 2 (LMCO), a pioneering cation‐disordered rocksalt oxide (DRX) cathode that has intrinsically small volume change upon lithium intercalation, can be integrated with a thiophosphate‐based solid electrolyte for all‐solid‐state batteries. Interface stability of the all‐solid LMCO cell was investigated by electrochemical impedance spectroscopy, X‐ray micro‐computed tomography, and electron microscopy. Since LMCO was initially synthesized as a layered phase exhibiting a large volume change, interface disintegration can be observable in the early cycles. As layered LMCO phase‐transformed into DRX LMCO in subsequent cycles, reintegration of the interfaces occurs within a pressurized cell as a result of its zero‐stain behavior. Consequently, the DRX LMCO cathode maintains interface integrity, and thus electrical wiring, over an extended number of cycles, leading to improved capacity retention with small internal cell resistance.

DRX↗

DOE OSTI2022

A Many-Body Perturbation Theory Approach to Energy Band Alignment at the Crystalline Tetracene–Silicon Interface

Hybrid inorganic–organic semiconductor (HIOS) interfaces are of interest for new photovoltaic devices operating above the Shockley–Queisser limit. Predicting energy band alignment at the interfaces is crucial for their design, but represents a challenging problem due to the large scales of the system, the energy precision required and a wide range of physical phenomena that occur at the interface. To tackle this problem, many-body perturbation theory in the non-self-consistent GW approximation, orbital relaxation corrections for organic semiconductors, and line-up potential method for inorganic semiconductors which allows for tractable and accurate computing of energy band alignment in crystalline van-der-Waals hybrid inorganic–organic semiconductor interfaces are used. In this work, crystalline tetracene physisorbed on the clean hydrogen-passivated 1 × 2 reconstructed (100) silicon surface is studied. Using this computational approach, it is found that the energy band alignment is determined by an interplay of the mutual dynamic dielectric screening of two materials and the formation of a dipole layer due to a weak hybridization of atomic/molecular orbitals at the interface. The significant role of the exchange-correlation effects in predicting band offsets for the hybrid inorganic–organic semiconductor interfaces is also emphasized.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2016

An implementation and evaluation of the MPI 3.0 one-sided communication interface

The Q1 Message Passing Interface (MPI) 3.0 standard includes a significant revision to MPI’s remote memory access (RMA) interface, which provides support for one-sided communication. MPI-3 RMA is expected to greatly enhance the usability and performance ofMPI RMA.We present the first complete implementation of MPI-3 RMA and document implementation techniques and performance optimization opportunities enabled by the new interface. Our implementation targets messaging-based networks and is publicly available in the latest release of the MPICH MPI implementation. Here using this implementation, we explore the performance impact of new MPI-3 functionality and semantics. Results indicate that the MPI-3 RMA interface provides significant advantages over the MPI-2 interface by enabling increased communication concurrency through relaxed semantics in the interface and additional routines that provide new window types, synchronization modes, and atomic operations.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2019

Propagation of Spin Waves Through an Interface Between Ferromagnetic and Antiferromagnetic Materials

Boundary conditions for order parameters at an interface between ferromagnetic (FM) and two-sublattice antiferromagnetic (AFM) materials are obtained in the continuous medium approximation similarly to the approach which allows one to take into account the finite thickness of the FM/FM interface, which is much less than spin wavelength. Three order parameters are considered inside an interface of finite thickness with the magnetization M of FM, magnetizations of both sublattices M{sub 1} and M{sub 2} of AFM. The uniform and non-uniform exchange between all order parameters are taken into account to the interface energy. Using these boundary conditions, the excitation of a surface evanescent spin wave is considered in AFM when the spin wave in FM falls onto this interface. The coefficients and the phases of transmission and reflection of spin wave through the FM/AFM interface are derived.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Geophysical methods reveal the soil architecture and subsurface stratigraphic heterogeneities across land-lake interfaces along Lake Erie

Abstract Purpose The land-lake interface is a unique zone where terrestrial and aquatic ecosystems meet, forming part of the Earth’s most geochemically and biologically active zones. The unique characteristics of this interface are yet to be properly understood due to the inherently high spatiotemporal variability of subsurface properties, which are difficult to capture with the traditional soil sampling methods. Geophysical methods offer non-invasive techniques to capture variabilities in soil properties at a high resolution across various spatiotemporal scales. Methods We combined electromagnetic induction (EMI), electrical resistivity tomography (ERT), and ground penetrating radar (GPR) with data from soil cores and in situ sensors to investigate hydrostratigraphic heterogeneities across land-lake interfaces along the western basin of Lake Erie. Results EMI revealed high spatial heterogeneities in ECa distribution across the land-lake interfaces, with higher values in the wetland and transition zones compared to the upland zone. Soil ECa maps matched soil maps from a public database with the hydric soil units delineated as high conductivity zones (ECa > 40 mS/m). ERT and GPR showed vertical variation in soil properties with clear stratigraphic boundaries, and correlation of ERT profiles with lithologs from piezometers revealed the stratigraphic units of silt–clay and till sequence down to 3.5 m depth which are consistent with the surficial geology of the study area. Conclusions These results validate the use of multiple geophysical methods for extrapolating soil properties and mapping stratigraphic structures at land-lake interfaces, thereby providing the missing information required to improve the earth system model (ESM) of coastal interfaces.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2015

Irradiation-induced formation of a spinel phase at the FeCr/MgO interface

Oxide dispersion strengthened ferritic/martensitic alloys have attracted significant attention for their potential uses in future nuclear reactors and storage vessels, as the metal/oxide interfaces act as stable high-strength sinks for point defects while also dispersing helium. Here, in order to unravel the evolution and interplay of interface structure and chemistry upon irradiation in these types of materials, an atomically sharp FeCr/MgO interface was synthesized at 500 °C and separately annealed and irradiated with Ni 3+ ions at 500 °C. After annealing, a slight enrichment of Cr atoms was observed at the interface, but no other structural changes were found. However, under irradiation, sufficient Cr diffuses across the interface into the MgO to form a Cr-enriched transition layer that contains spinel precipitates. First-principles calculations indicate that it is energetically favorable to incorporate Cr, but not Fe, substitutionally into MgO. Furthermore, our results indicate that irradiation can be used to form new phases and complexions at interfaces, which may have different radiation tolerance than the pristine structures.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Effects of three-dimensional Cu/Nb interfaces on strengthening and shear banding in nanoscale metallic multilayers

Manipulation of the atomic-scale structure of two-dimensional (2D) interfaces have been shown to provide nanocomposites with enhanced strength, deformability, and radiation damage resistance. In comparison with 2D interfaces, here we investigate the mechanical response of nanocomposites containing three-dimensional (3D) Cu/Nb interfaces consisting of a chemical/structural gradient separating pure Cu and Nb layers, through which the lattice mismatch between face-centered cubic Cu and body-centered cubic Nb is accommodated over a distance of several nanometers. Furthermore, it is demonstrated that 3D interfaces increase the yield and flow strength by 50% and 22%, respectively, over composites containing 2D interfaces at similar layer thicknesses. After 14% compressive strain, the onset of shear banding results in co-deformation of both Cu and Nb phases within and outside of the shear band. We conclude with a discussion of the role of interface structure in shear band formation and growth in 3D Cu/Nb.

3-dimensional interface↗

DOE OSTI2022

Characterization of the atomic-level structure of γ-alumina and (111) Pt/γ-alumina interfaces

The atomic-level structure of platinum/γ-alumina interfaces is characterized in a model system of dense γ-alumina embedded with faceted Pt NPs produced by implantation of platinum ions into sapphire followed by thermal annealing in air at 800 °C. Aberration-corrected scanning transmission electron microscopy (STEM) was used to collect atomic-resolution images, which are compared to STEM image simulations of two experimentally-based bulk models of γ-alumina by Smrčok et al. and Zhou and Snyder. A density functional theory (DFT) based model of (111) interfaces with different chemical terminations (O, Al 1 , Al 2 ) of the γ-alumina developed by Oware Sarfo et al. is also compared to experimental STEM data from Pt/γ-alumina interfaces. Further, the Smrčok γ-alumina model provides a better fit than the Zhou structure to the bulk of the γ-alumina. The oxygen-terminated Oware Sarfo model best fits the experimental data and is a very good model close to the interface. However, the fit of the interface model to the experimental data is poorer beyond the third atomic layer in the γ-alumina. This is attributed to compromises required in the design of the model to limit the cell size and computational time for DFT calculations. Understanding the accuracy and limits of the structural models of γ-alumina and Pt/γ-alumina interfaces is important to further the understanding of the structure/property relationships in this system.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Multi-material hydrodynamics with algebraic sharp interface capturing

A finite volume method for Eulerian multi-material hydrodynamics with sharp interface capturing is presented here. The pressure-temperature non-equilibrium multi-material system with finite-rate pressure relaxation in mixed-cells is considered here. This pressure closure facilitates material-property-dependent pressure relaxation, rather than instantaneous pressure equilibration, which in turn allows the use of unsplit high-order time-integrators. A modified tangent of hyperbola for interface capturing (THINC) method is used to reconstruct multi-material (> 2) interfaces, on three-dimensional unstructured meshes. A simple modification which extends the THINC reconstruction to interfaces between more than two materials is proposed. It is demonstrated that the modified THINC can capture multi-material interfaces within 2–4 tetrahedral cells. Since no geometric reconstructions are required by the THINC method, the presented multi-material method is algorithmically simple, and computationally efficient. Consistent reconstructions of conserved quantities at material interfaces ensure that conservation and closure laws are satisfied at the discrete level. Through a suite of test problems solved on unstructured meshes, it is demonstrated that the presented method is a promising candidate for accurate and efficient multi-material hydrodynamics computations.

97 MATHEMATICS AND COMPUTING↗