Engineering Papers⌕ Search

SEARCH · Engineering Papers

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

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

At least 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Phase-field modeling of diffusion bonding in 316H stainless steel: Impact of processing conditions on grain morphology and bonding quality

A novel multi-phase, multi-component phase‐field model is presented to study the diffusion bonding of 316H stainless steel. Combined with targeted experimental investigations, this model simulates the bond-growth process and predicts the bonding quality. Unlike previous models, our approach captures the simultaneous evolution of voids and grain structures, while quantifying bonding quality using defined bonding ratio. A comprehensive analysis of bond process control is performed by changing temperature, pressure and surface roughness observing the resulting bond structure, which is consistent with experimental observations and analytical predictions. Temperature is determined to be the dominant factor, with the transition from a flat to a robust bond occurring between 1000 °C and 1050 °C. At the ideal bonding temperature of 1050 °C, a surface roughness exceeding 0.6 μm or an applied stress below 4 MPa results in poor bonding quality. Beyond this, higher pressures and smoother surfaces reduce void size, accelerate void shrinkage, and lead to improved bond integrity. This diffuse-interface model can be extended to other material systems if supplied with appropriate thermodynamic and kinetic data. In conclusion, this makes it an effective modeling platform for optimizing high-temperature diffusion bonding and developing reliable bonded components such as compact heat exchangers.

Diffusion bonding↗
DOE OSTI2022

FSISPH: An SPH formulation for impacts between dissimilar materials

Here, we present an SPH formulation with several new features designed to better model the fully-compressible interaction of dissimilar materials. We developed the new method to simulate the atmospheric entry and break-up of small celestial bodies in planetary atmospheres. The formulation uses a unity-based, density-energy discretization of the hydrodynamic conservation laws with linear-corrected kernel gradients. To account for variations in compressibility, we use an HLLC approximate Riemann solver to adjust the velocity gradient at material interfaces. To handle large transverse velocity discontinuities, we introduce a simple slip interface model that limits the artificial viscosity at material interfaces. Diffusion is optionally applied through the velocity gradient and this allows the density and specific thermal energy to evolve in a manner more consistent with the first law of thermodynamics in comparison to other more direct diffusion schemes. We also introduce a material-local second-order artificial conduction scheme used to smooth the specific thermal energy field. Material damage fits neatly under this framework by treating the damage front as a material interface. The method has been implemented as a solver, FSISPH, within the code, Spheral++, and is publicly available on github. We test our new solver on a number of classic shock, mixing, and multi-material problem. The components we outline can significantly improve accuracy of SPH for problems with sharp contact discontinuities.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

A regularized-interface method as a unified formulation for simulations of high-pressure multiphase flows

The injection of multi-species fluids into high-pressure and high-temperature environments beyond the species' critical points is commonly found in engineering applications. At these conditions, for immiscible species, both subcritical interfacial dynamics and supercritical mixing can coexist due to variations in temperature around the mixture critical point. The modeling of these complex transcritical phenomena for large-scale configurations is so far not possible. To address this issue, we propose the Regularized-Interface Method (RIM) as a unified formulation that can describe both sub- and supercritical processes as well as the transition between them. The proposed method is derived via filtering of the nanoscale interface-resolving formulation based on van der Waals' linear gradient theory. Thus, this approach allows for the consistent modeling of interfacial dynamics that vanishes at supercritical conditions, while significantly reducing the temporal and spatial resolution constraints of the original nanoscale formulation. The resulting RIM formulation is examined in interface-capturing simulations of sub-, trans-, and supercritical fuel injection processes, involving droplets and jets. Furthermore, these results highlight the importance of resolving spatio-temporal transitions from subcritical interfacial dynamics to supercritical mixing in high-pressure multiphase simulations, in contrast to commonly employed diffused-interface methods, where interfacial dynamics are often neglected.

Interface capturing↗
DOE OSTI2021

On a nonlocal Cahn–Hilliard model permitting sharp interfaces

A nonlocal Cahn–Hilliard model with a non-smooth potential of double-well obstacle type that promotes sharp interfaces in the solution is presented. To capture long-range interactions between particles, a nonlocal Ginzburg–Landau energy functional is defined which recovers the classical (local) model as the extent of nonlocal interactions vanish. In contrast to the local Cahn–Hilliard problem that always leads to diffuse interfaces, the proposed nonlocal model can lead to a strict separation into pure phases of the substance. In this work, the lack of smoothness of the potential is essential to guarantee the aforementioned sharp-interface property. Mathematically, this introduces additional inequality constraints that, in a weak formulation, lead to a coupled system of variational inequalities which at each time instance can be restated as a constrained optimization problem. We prove the well-posedness and regularity of the semi-discrete and continuous in time weak solutions, and derive the conditions under which pure phases are admitted. Moreover, we develop discretizations of the problem based on finite element methods and implicit–explicit time-stepping methods that can be realized efficiently. Finally, we illustrate our theoretical findings through several numerical experiments in one and two spatial dimensions that highlight the differences in features of local and nonlocal solutions and also the sharp interface properties of the nonlocal model.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1998

X-Ray Photoelectron Spectroscopy Study of the Heating Effects on Pd/6H-SiC Schottky Structure

X-ray photoelectron spectroscopy is used to study the effects of heat treatment on the Pd/6H-SiC Schottky diode structure. After heating the structure at 425 C for 140 h, a very thin surface layer of PdO mixed with SiO(x) formed on the palladium surface of the Schottky structure. Heat treatment promoted interfacial diffusion and reaction which significantly broadened the interfacial region. In the interfacial region, the palladium concentration decreases with depth, and the interfacial products are Pd(x)Si (x = 1,2,3,4). In the high Pd concentration regions, Pd4Si is the major silicide component while gr and Pd2Si are major components in the low Pd concentration region. At the center of the interface, where the total palladium concentration equals that of silicon, the concentrations of palladium associated with various palladium silicides (Pd(x)Si, x= 1,2,3,4) are approximately equal. The surface passivation layer composed of PdO and SiO, may significantly affect the electronic and catalytic properties of the surface of the Schottky diode which plays a major role in gas detection. The electronic properties of the Schottky structure may be dominated by a (Pd+Pd(x)Si)/SiC interface. In order to stabilize the properties of the Schottky structure the surface and interface diffusion and reactions must be controlled.

Chen, Liang-Yu↗
DOE OSTI2020

Deep learning of interface structures from simulated 4D STEM data: cation intermixing vs. roughening ∗

Abstract Interface structures in complex oxides remain an active area of condensed matter physics research, largely enabled by recent advances in scanning transmission electron microscopy (STEM). Yet the nature of the STEM contrast in which the structure is projected along the given direction precludes separation of possible structural models. Here, we utilize deep convolutional neural networks (DCNN) trained on simulated 4D STEM datasets to predict structural descriptors of interfaces. We focus on the widely studied interface between LaAlO 3 and SrTiO 3 , using dynamical diffraction theory and leveraging high performance computing to simulate thousands of possible 4D STEM datasets to train the DCNN to learn properties of the underlying structures on which the simulations are based. We test the DCNN on simulated data and show that it is possible (with >95% accuracy) to identify a physically rough from a chemically diffuse interface and create a DCNN regression model to predict step positions. We quantify the applicability of the model to different thicknesses and the transferability of the approach. The method shown here is general and can be applied for any inverse imaging problem where forward models are present.

42 ENGINEERING↗
NASA NTRS2000

Diffusing Wave Spectroscopy Used to Study Foams

The white appearance of familiar objects such as clouds, snow, milk, or foam is due to the random scattering of light by the sample. As we all know, pure water is clear and easily passes a beam of light. However, tiny water droplets, such as those in a cloud, scatter light because the air and water droplet have different indexes of refraction. When many droplets, or scattering sites, are present, the incident light is scattered in random directions and the sample takes on a milky white appearance. In a glass of milk, the scattering is due to small colloidal particles. The white appearance of shaving cream, or foam, is due to the scattering of light at the water-bubble interface. Diffusing wave spectroscopy (DWS) is a laser light-scattering technique used to noninvasively probe the particle dynamics in systems that strongly scatter light. The technique takes advantage of the diffuse nature of light, which is reflected or transmitted from samples such as foams, dense colloidal suspensions (such as paint and milk), emulsions, liquid crystals, sandpiles, and even biological tissues.

Zimmerli, Gregory A.↗
DOE OSTI2020

Laser powder bed fusion of Inconel 718 on 316 stainless steel

The joining of Inconel 718 and 316 stainless steel has notable industrial importance. However, traditional arc-welding has many issues. In particular, the formation of the detrimental Nb-rich Laves phase in the heat affected zone can promote crack initiation and propagation, which leads to poor mechanical properties. The small heat input and the rapid solidification of laser powder bed fusion (PBF) has the potential to suppress this problem. Despite the infancy of the multi-material PBF, it is of interest to investigate the relationship among the process parameter, the microstructure at the materials joining point, and the mechanical performance. In this study, we used high-speed x-ray imaging to in-situ monitor the printing of Inconel 718 tracks on a 316H SS substrate via laser powder bed fusion, followed by post-build characterization and testing with SEM, XRD and nanoindentation. We observed that a higher linear energy density resulted in a more diffused interface and a smaller compositional gradient, but as a trade-off, generated more keyhole pores and thermal cracks. Furthermore, we also showed that the built layer exhibited chemical inhomogeneity at two length scales: a tens-of-mu m-scale compositional vortex due to the melt flow, and a sub-mu m-scale Nb-rich dendritic structure due to the non-equilibrium solidification that led to the formation of the Laves phase. It was inferred that the chemical inhomogeneity may benefit the mechanical properties at the interface by providing more interlocking between the two materials and also by limiting the effect of the brittle Laves phase to potentially short range.

36 MATERIALS SCIENCE↗
DOE OSTI2025

An end-to-end deep learning method for solving nonlocal Allen–Cahn and Cahn–Hilliard phase-field models

Here, we propose an efficient end-to-end deep learning method for solving nonlocal Allen–Cahn (AC) and Cahn–Hilliard (CH) phase-field models. One motivation for this effort emanates from the fact that discretized partial differential equation-based AC or CH phase-field models result in diffuse interfaces between phases, with the only recourse for remediation is to severely refine the spatial grids in the vicinity of the true moving sharp interface whose width is determined by a grid-independent parameter that is substantially larger than the local grid size. In this work, we introduce non-mass conserving nonlocal AC or CH phase-field models with regular, logarithmic, or obstacle double-well potentials. Because of non-locality, some of these models feature totally sharp interfaces separating phases. The discretization of such models can lead to a transition between phases whose width is only a single grid cell wide. Another motivation is to use deep learning approaches to ameliorate the otherwise high cost of solving discretized nonlocal phase-field models. To this end, loss functions of the customized neural networks are defined using the residual of the fully discrete approximations of the AC or CH models, which results from applying a Fourier collocation method and a temporal semi-implicit approximation. To address the long-range interactions in the models, we tailor the architecture of the neural network by incorporating a nonlocal kernel as an input channel to the neural network model. We then provide the results of extensive computational experiments to illustrate the accuracy, predictive capabilities, and cost reductions of the proposed method.

42 ENGINEERING↗
DOE OSTI2020

Investigation of the discontinuous precipitation of U-Nb alloys via thermodynamic analysis and phase-field modeling

U-Nb’s discontinuous precipitation, $\gamma^{bcc}_{matrix}$, $\alpha^{orth}_{cellular}$, $\gamma\prime^{bcc}_{cellular}$, is intriguing in the sense that it allows formation and growth of the metastable $\gamma\prime$ phase during the course of its occurrence. Previous attempts to explain the thermodynamic origin of U-Nb’s discontinuous precipitation hypothesized that the energy of $\alpha$ forms an intermediate common tangent with the first potential of the double-well energy of $\gamma$ at the $\gamma\prime$ composition. In this work, we examine different possible mechanisms by which the discontinuous precipitation product in the U-Nb system can be stabilized. We put forward a mechanism by which the bulk free energy of the $\gamma$ can develop a non-monotonic dependency with composition. Additionally we posit that local contributions due to lattice mismatch between the $\alpha$ and $\gamma$ phases may be responsible for the generation of metastable states that may stabilized by thermodynamics as well as by kinetics. Our work suggests that local misfit strain tends to play a crucial role in the growth of the discontinuous precipitation product. Depending on the magnitude of strain developed at the $\alpha$/$\gamma\prime$ interfaces, either an increasing $\gamma\prime$ composition or a random distribution of $\gamma\prime$ compositions around the equiatomic value with respect to increasing temperature could be expected. Moreover, we show how it is possible to stabilize the discontinuous precipitation front through highly anisotropic and fast interface diffusion.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Modeling discontinuous dynamic recrystallization containing second phase particles in magnesium alloys utilizing phase field method

Pre-existing second phase particles can induce particle stimulated nucleation (PSN) and pinning effect on grain boundaries during DRX process of magnesium alloys. The interaction among pinning effect, PSN mechanism and strain-induced grain boundary migration (SIBM) nucleation mechanism imposes challenge in quantitative study of microstructure evolution. A phase field model that describes discontinuous dynamic recrystallization process considering second phase particles (PF-DDRX-SP) which mainly distribute along grain boundaries is developed. The second phase order parameter with diffusion interface is introduced to the free energy function. With considering the pinning effect of second phase particles on grain boundaries, the grain boundary energy couples the interaction between second phase and grains. The effect of second phase on dislocation density evolution is considered in the model, and the PSN and SIBM mechanism are both executed by limiting the nucleation position to the grain boundary and the second phase interface. For the validation, the initial topology consistent with initial microstructure is applied to the simulation of DRX kinetics and flow behavior of AZ80 magnesium alloy, and great reliability and accuracy of developed model is indicated. Furthermore, through the model, the microstructure evolution with different particle size and with different volume fraction of round particle during DRX process were simulated. Finally, the results show that particles with higher dispersion reduce the number of DRX nuclei, which indicate that second phase particles suppress the nucleation at grain boundaries because of the pinning effect. More importantly, the proposed model could also quantitatively predict the relationships between the parameters of second phase particles and DRX behaviors, and enable to optimize the initial second phase structure in a uniform grain structure during thermomechanical process.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Mesoscale modeling and semi-analytical approach for the microstructure-aware effective thermal conductivity of porous polygranular materials

Here we established a comprehensive modeling approach for investigating the microstructure-aware effective thermal conductivity ($κ_{eff}$) for porous microstructures containing solid particles and gaseous pores. Our approach combines the mesoscale computational modeling framework and the semi-analytical method, allowing for efficient prediction of $κ_{eff}$ for realistic porous microstructures, while considering complicated microstructural thermal conduction pathways effectively in the prediction. We used the diffuse-interface mesoscale computational model to generate extensive simulated $κ_{eff}$ data for realistic digital representations of microstructures with wide ranges of porosity ($f_p$), thermal conductivity of the gas phase ($κ_g$), and thermal conductivity of the solid phase ($κ_s$). From the simulated data, we identified two property variation regimes for $κ_{eff}$: (1) a slow $κ_{eff}$ increase for $κ_s ~ κ_g$; and (2) a faster $κ_{eff}$ increase for $κ_s \gg κ_g$. To capture the key features of the relationship between the microstructure and $κ_{eff}$, we derived a semi-analytical model by introducing structure and intensification factors. The two new factors incorporate the calibrated effective contribution of the solid volume with $κ_s$ and additional interfacial effects into the prediction of $κ_{eff}$, respectively, allowing for consideration of parallel, serial, and interfacial conduction mechanisms effectively. Using the selected simulation data, we quantified key model parameters within the semi-analytical model and verified that the parameterized model exhibits excellent agreement with simulated $κ_{eff}$ for the entire range of the parameter space.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Analysis of a new multispecies tumor growth model coupling 3D phase-fields with a 1D vascular network

In this study, we present and analyze a mathematical model for tumor growth incorporating ECM erosion, interstitial flow, and the effect of vascular flow and nutrient transport. The model is of phase-field or diffused-interface type in which multiple phases of cell species and other constituents are separated by smooth evolving interfaces. The model involves a mesoscale version of Darcy’s law to capture the flow mechanism in the tissue matrix. Modeling flow and transport processes in the vasculature supplying the healthy and cancerous tissue, one-dimensional (1D) equations are considered. Since the models governing the transport and flow processes are defined together with cell species models on a three-dimensional (3D) domain, we obtain a 3D–1D coupled model.

3D–1D coupled blood flow models↗
DOE OSTI2020

Leveraging Sequential Doping of Semiconducting Polymers to Enable Functionally Graded Materials for Organic Thermoelectrics

With the ability to modulate electronic properties through molecular doping coupled with ease in processability, semiconducting polymers are at the forefront in enabling organic thermoelectric devices for thermal energy management. In contrast to uniform thermoelectric material properties, an alternative route focuses on functionally graded materials (FGMs) where one spatially controls and optimizes transport properties across the length of a thermoelectric material. While primarily studied in the context of inorganic materials, the concept of FGMs for organic thermoelectrics has not been explored. Herein, we introduce how molecular doping of semiconducting polymers enables spatial compositional control of thin-film FGMs. Specifically, we utilize sequential vapor doping of poly[2,5-bis(3-tetradecylthiophen-2-yl) thieno [3,2-b]thiophene] (PBTTT) with the small molecule acceptor 2,3,5,6-tetrafluoro-7,7,8,8-tetracyanoquinodimethane (F4TCNQ) to fabricate the simplest form of FGMs – double- segmented thin films. The two thin-film segments are of equal length (7.5 mm), but each set to different doping levels. Our study focuses on understanding the thermoelectric properties (Seebeck coefficient, α, and electronic conductivity, σ) and structural properties (through X-ray scattering, UV-vis-NIR spectroscopy, and Raman spectroscopy) within and across the two segments. We observe the presence of a small diffuse interfacial region of 0.5 – 1 mm between the two segments where the doping level and transport properties vary continuously. Despite the diffuse interface, the measured effective Seebeck coefficient (αeff) across the two segments is simply the average of α within each segment. Importantly, this experimental result is consistent with reported mathematical models describing αeff in graded thermoelectric materials. Finally, our results demonstrate the facile fabrication and characterization of functionally graded organic thermoelectric materials, providing guidelines for further development on more complex FGMs.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Fluctuating hydrodynamics and the Rayleigh–Plateau instability

The Rayleigh–Plateau instability occurs when surface tension makes a fluid column become unstable to small perturbations. At nanometer scales, thermal fluctuations are comparable to interfacial energy densities. Consequently, at these scales, thermal fluctuations play a significant role in the dynamics of the instability. These microscopic effects have previously been investigated numerically using particle-based simulations, such as molecular dynamics (MD), and stochastic partial differential equation–based hydrodynamic models, such as stochastic lubrication theory. In this paper, we present an incompressible fluctuating hydrodynamics model with a diffuse-interface formulation for binary fluid mixtures designed for the study of stochastic interfacial phenomena. An efficient numerical algorithm is outlined and validated in numerical simulations of stable equilibrium interfaces. We present results from simulations of the Rayleigh–Plateau instability for long cylinders pinching into droplets for Ohnesorge numbers of Oh = 0.5 and 5.0. Both stochastic and perturbed deterministic simulations are analyzed and ensemble results show significant differences in the temporal evolution of the minimum radius near pinching. Short cylinders, with lengths less than their circumference, were also investigated. As previously observed in MD simulations, we find that thermal fluctuations cause these to pinch in cases where a perturbed cylinder would be stable deterministically. Finally, we show that the fluctuating hydrodynamics model can be applied to study a broader range of surface tension–driven phenomena.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Anisotropic Thermal Transport in Quasi-2D Ruddlesden-Popper Hybrid Perovskite Superlattices

Two-dimensional hybrid metal-halide perovskites (2D-MHPs) have emerged as important solution-processed semiconductors with favorable optical and electronic properties for diverse applications in photovoltaics, optoelectronics, and spintronics. The quasi-2D layered structures, featuring large acoustic impedance mismatches between the organic and inorganic sublattices, are expected to result in distinct and anisotropic thermal transport properties along the cross-plane and in-plane directions. Here, we introduce transducer-free vibrational-pump-visible-probe (VPVP) approaches that enable accurate quantification of anisotropic thermal transport properties in various archetypical single-crystalline 2D-MHPs. Specifically, using VPVP spectroscopy and VPVP microscopy, we measure the anisotropic thermal diffusivities of 2D-MHPs with systematically varied Pb-I octahedral layer thicknesses, as well as organic spacer types and lengths, revealing how these structural parameters alter the cross-plane and in-plane thermal transport properties in distinct ways. While diffuse interface scattering plays an important role in dictating cross-plane thermal transport, in-plane thermal transport is primarily determined by phonon transport within interconnected inorganic layers. Density functional theory incorporating four-phonon scatterings provides further insight into the low thermal conductivity and modest thermal conduction anisotropy in 2D-MHPs. Our work demonstrates a new all-optical and noncontact method, which requires minimal sample preparation and allows direct visualization of cross-plane and in-plane thermal transport, potentially compatible with sample environments. The demonstrated VPVP approaches can advance understanding of thermal transport in 2D-MHPs as well as wide-ranging hybrid and polymeric semiconductors beyond 2D-MHPs.

36 MATERIALS SCIENCE↗
NASA NTRS1981

Fabrication of titanium multiwall Thermal Protection System (TPS) curved panel

A method is described for fabricating curved T1-6Al-4V multiwall panels having a single radius of curvature of 30.5 cm (12 inches). A panel 30.5 x 30.5 x 1.7 cm (12 x 12 x .688 inch) having 4 dimpled sheets, 5 undimpled sheets, and stepped side closures was designed, fabricated and delivered for evaluation. Tool design for skin and dimpled sheet forming are described as well as layup and liquid interface diffusion bonding.

Blair, W.↗
NASA NTRS1982

Single-mode theory of diffusive layers in thermohaline convection

A two-layer configuration of thermohaline convection is studied, with the principal aim of explaining the observed independence of the buoyancy-flux ratio on the stability parameter when the latter is large. Temperature is destabilizing and salinity is stabilizing, so diffusive interfaces separate the convecting layers. The convection is treated in the single-mode approximation, with a prescribed horizontal planform and wavenumber. Surveys of numerical solutions are presented for a selection of Rayleigh numbers R, stability parameters lambda and horizontal wavenumbers. The solutions yield a buoyancy flux ratio chi that is insensitive to lambda, in accord with laboratory experiments. However chi increases with increasing R, in contradiction to laboratory observations.

Gough, D. O.↗