Search NASA⌕ Search

SEARCH · Search NASA

Results for “Sharp 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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
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 OSTI2022

3rd harmonic magnetometry assessment of NbTiN-based SIS structures

In the quest for alternative superconducting materials to bring accelerator cavity performance beyond the bulk niobium (Nb) intrinsic limits, a promising concept proposes that superconductor-insulator-superconductor (SIS) thin film structures can delay magnetic flux penetration in accelerator cavities to higher fields [1]. NbTiN is a candidate superconductor for such structures. We have demonstrated high quality NbTiN and AlN deposited by reactive direct current magnetron sputtering (DCMS), both for individual layers and multilayers. Interface quality has been assessed for bi-layer stacks with various NbTiN and AlN thicknesses from 500 and 30 nm down to 3 and 1 nm. These SIS structures show continued sharp interfaces. The Hfp enhancement of the films was examined with 3rd harmonic magnetometry. The system was designed and built in an ongoing collaboration with CEA Saclay. It can measure 1? to 2? samples on a temperature controlled stage. This contribution presents the assessment of the first penetration field enhancement with 3rd harmonic magnetometry for standalone films and multilayer nanostructures.

Valente, Anne-Marie↗
DOE OSTI2025

Multi-material ALE remap with interface sharpening using high-order matrix-free finite element methods

The arbitrary Lagrangian-Eulerian (ALE) technique involves remapping field quantities from a Lagrangian mesh to an optimized mesh in a conservative, accurate and bounds-preserving manner. For methods based on arbitrary order finite elements, as described in a reference, material volume fractions are advected in pseudo-time using flux-corrected transport (FCT) without any form of interface reconstruction. In practice, this can lead to excessive propagation of small volume fractions throughout the domain. In addition, this method requires assembly of a global advection matrix to compute the bounds-preserving low-order FCT solution. In this work, we introduce a new approach for ALE remap using a high-order matrix-free technique which incorporates a flux modification to sharpen material interfaces in a conservative manner. Our approach begins with computing a bounds-preserving low-order solution to the ALE remap equations at the element level. We then compute a sharp interface solution (not guaranteed to be bounds-preserving) which comes from solving an augmented version of the ALE remap equations with a conservative flux modification which acts to sharpen material volume fractions based on their gradients and transport directions. Using the sharp interface solution, we make global corrections to the bounds-preserving solution while maintaining preservation of bounds. By blending with the sharpened solution at the global level we are able to globally conserve mass without hindering the remap pseudo-time step. This new interface-aware ALE remap method is based entirely on partial assembly techniques where globally assembled matrix operators are no longer needed, resulting in a globally matrix-free FCT method for multi-material, multi-field ALE remap with high performance on GPU architectures. We present results of our new remap method on 1D, 2D and 3D benchmarks and describe the algorithmic tailoring for GPU architectures that was developed.

Vargas, Arturo [Lawrence Livermore National Labora↗
NASA NTRSDate not supplied

CFD Model Development of a Cryogenic Storage Tank Self-Pressurization in Normal Gravity and Validation against SHIIVER Experiment

Two-phase flow and heat transfer simulations with interfacial phase change of the Structural Heat Intercept, Insulation, and Vibration Evaluation Rig (SHIIVER) self-pressurization experiment were conducted using storage tank CFD model in the framework of the ANSYS Fluent CFD code. The simulations were performed for the 70% fill level case with MLI on domes and no vapor cooling. All the phase change calculations in these simulations were generated by in-house Schrage-based evaporation-condensation model. The calculations were performed using both the Volume of Fluid (VOF) and Sharp Interface multiphase 2D axisymmetric models. A number of parametric and sensitivity studies were performed to check the various aspects of the CFD model. These studies helped to understand the effects of varying several parameters on the tank pressure and temperature during self-pressurization. Turbulence modeling; turbulence damping at the interface; and using constant vs. temperature dependent fluid properties were shown to have the most profound influence on predicted tank pressures and temperatures. Current study indicates that including phase change at the interface into the computational model is crucial for accurate prediction of the tank self-pressurization process. The effect of accommodation coefficient was also studied. Tank pressure values predicted by the VOF, and Sharp Interface models are within 4% of the experimental ones.

Computational Fluid Dynamics↗
DOE OSTI2022

Smoothed boundary method for simulating incompressible flow in complex geometries

Simulating flow through porous media with explicit considerations of complex microstructures is very challenging using conventional sharp-interface methods because of the difficulties in generating meshes conformal to complex geometries. In this work, a diffuse interface embedded boundary method known as the Smoothed Boundary Method (SBM) is utilized to facilitate simulations of fluid dynamics involving complex geometries. In diffuse-interface methods, the geometry is described by a domain parameter. The SBM allows the straightforward reformulation of the time-dependent Navier-Stokes equations in terms of this domain parameter, using only algebraic identities. Thus, enforcing the appropriate boundary conditions at the irregular embedded boundary is greatly simplified. Adaptive mesh refinement is used to increase the accuracy of the diffuse interface method by allowing thinner interfacial thickness to be used in the domain parameter. Here, the SBM-formulated Navier-Stokes equations are solved with the Finite Difference Method on refined mesh systems. Sharp-interface Finite Element Method simulations using the commercial software COMSOL on body-conforming meshes are also provided for comparison. Favorable agreement between the two methods is observed. Since it is no longer necessary for the mesh to conform to the complex geometry, the grid system for the SBM simulations can be generated rapidly and without additional manual interventions, making the entire simulation process more expedient.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Twin nucleation and growth in hexagonal close-packed metals: The role of slip-mediated plasticity on twin embryo formation and evolution

Twinning is a key deformation mechanism in hexagonal close-packed (hcp) metals, which are typified by a lack of easily activated slip systems that can accommodate a general state of loading. Pragmatically, the nucleation and evolution of twin domains occurs concomitantly with slip, such that the eventual twin network is conditioned by both external and internal stresses resulting, among others, from the evolution of dislocations. However our understanding of the interplay between dislocations, and twin nucleation and stability remains limited. This work focuses on elucidating the influence of dislocation-mediated plasticity on the formation, growth, and stability of twin embryos. First, a new mesoscale spectral crystal plasticity-twinning framework is developed and used to quantify the change in the free energy landscape following the nucleation {$10\bar{1}2$} of twins in Mg. Representative twin morphologies are modeled under conditions of limited and profuse slip activity, which are emulative of small- and bulk-scale samples, respectively. Then, the driving traction profiles around twin embryos are investigated via a sharp interface approach to obtain insights into how concurrent slip-mediated plasticity can influence the growth/stabilization of nanometric twins. The driving traction profiles are further utilized to determine the stability of twin embryos post loading. The initial dislocation density in the samples, and within the different domains (i.e., twin vs. parent), is seen to have a significant effect on the twin nucleation stress. Namely, the activation of high levels of concomitant plasticity in the parent grain is seen to significantly drive the formation of nanometric twin nuclei at stresses as low as 250 MPa. Further, the sharp interface analysis reveals that profuse plasticity in the parent grain simultaneously alters the forward and back stresses, such that the magnitude/polarity of the driving tractions become increasingly favorable for nanometric twin growth when slip is active. Finally, prior plasticity in the parent grain is seen to result in favorable driving tractions for nanometric twin growth, even at applied stresses as low as ~ 100 MPa. In conclusion, these results are in stark contrast to a case without any dislocations, wherein applied stresses as high as ~ 500 MPa are necessary to grow the twin domains.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Multilayer interface tracking model of pure tungsten oxidation

Here, we present a numerical model to predict oxide scale growth on tungsten surfaces under exposure to oxygen at high temperatures. The model captures the formation of four thermodynamically-compatible oxide sublayers, WO 2 , WO 2.72 , WO 2.9 , and WO 3 , on top of the metal substrate. Oxide layer growth is simulated by tracking the oxide/oxide and oxide/metal interfaces using a sharp-interface Stefan model coupled to diffusion kinetics. The model is parameterized using selected experimental measurements and electronic structure calculations of the diffusivities of all the oxide subphases involved. We simulate oxide growth at temperatures of 600°C and above, extracting the power law growth exponents in each case, which we find to deviate from classical parabolic growth in several cases. We conduct a comparison of the model predictions with an extensive experimental data set, with reasonable agreement at most temperatures. While many gaps in our understanding still exist, this work is a first attempt at embedding the thermodynamic and kinetic complexity of tungsten oxide growth into a comprehensive mesoscale kinetic model that attempts to capture the essential features of tungsten oxidation to fill existing knowledge gaps and guide and enhance future tungsten oxidation models.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Nonequilibrium statistical thermodynamics of multicomponent interfaces

Nonequilibrium interfacial thermodynamics has important implications for crucial biological, physical, and industrial-scale transport processes. Here, we discuss a theory of local equilibrium for multiphase multicomponent interfaces that builds upon the “sharp” interface concept first introduced by Gibbs, allowing for a description of nonequilibrium interfacial processes such as those arising in evaporation, condensation, adsorption, etc. By requiring that the thermodynamics be insensitive to the precise location of the dividing surface, one can identify conditions for local equilibrium and develop methods for measuring the values of intensive variables at the interface. We then use extensive, high-precision nonequilibrium molecular dynamics (NEMD) simulations to verify the theory and establish the validity of the local equilibrium hypothesis. In particular, we demonstrate that equilibrium equations of state are also valid out of equilibrium, and can be used to determine interfacial temperature and chemical potential(s) that are consistent with nonequilibrium generalizations of the Clapeyron and Gibbs adsorption equations. We also show, for example, that, far from equilibrium, temperature or chemical potential differences need not be uniform across an interface and may instead exhibit pronounced discontinuities. However, even in these circumstances, we demonstrate that the local equilibrium hypothesis and its implications remain valid. These results provide a thermodynamic foundation and computational tools for studying or revisiting a wide variety of interfacial transport phenomena.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Advanced Computational Modeling of High-Level Waste Vitrification at the Hanford Site

The U.S. Department of Energy (DOE) has selected vitrification for stabilizing legacy tank waste at the Hanford site, where radioactive waste from plutonium production was historically stored in underground tanks. This waste will be separated into low-activity waste (LAW) and high-level waste (HLW) fractions and processed at the Waste Treatment and Immobilization Plant (WTP). At WTP, glass melters are used for the vitrification of radioactive tank waste, transforming it into a stable borosilicate glass form for safe long-term storage. The melter vessel is constructed from highly durable and heat-resistant materials, where the vitrification process occurs. The main regions that are modeled are the melt pool, plenum, cold cap, riser/discharge chamber, and surrounding structure with insulation layers. Forced convection induced by air bubblers at the base of the melter ensure uniform temperature distribution and provide heat to the cold cap layer. The cold cap is a region of reacting batch feed that floats on top of the molten glass and is where the batch-to-glass reactions occur. Joule heating provided by electrodes mounted along the vertical walls of the melter and immersed directly in the glass, generates the necessary heat for the net endothermic conversion processes that occur in the cold cap. The high temperatures, radioactivity, and opaque nature of the glass prevent direct observation inside the melters. Therefore, computational models are essential for providing insight into factors that affect melter throughput. Thermocouples in the plenum provide operators with plenum temperature measurements. Operational adjustments include bubbling rate, voltage supplied to the electrodes, feed adjustments, and glass removal rate. Different computational fluid dynamics (CFD) models have been developed, each serving a specific purpose. There are CFD models of different scale melters, as well as models that capture the two-phase flow interfaces of rising bubbles in the molten glass or models with a simplified molten glass region so that the surrounding structure and plenum can be feasibly incorporated. Pilot-scale melter models have been developed to serve as validation of the methods employed in the simulation of the full-scale WTP melters. Models incorporating resolved bubbling are used to develop momentum source terms to implement into a single phase, multi-region, steady-state flow model that is being validated by measured process parameters such as glass production rate, voltage, input power, plenum temperatures, etc. The resolved bubbling model uses the multiphase volume of fluid approach to model the system with a high-resolution interface capturing scheme to maintain sharp interfaces between the molten glass and the air phase. The suite of CFD models is continually being improved to incorporate more realistic physics and achieve faster turnaround time. For example, an incremental controller is implemented to automatically adjust electrode voltage within the simulation to a molten glass set point temperature of 1150°C. Newer models feature improved meshes to ensure conformal meshes between regions and eliminate unnecessary mesh refinement in areas that are not of interest (such as boundary layers in offgas ports). Instead of explicitly modeling the structural, refractory, and insulation layers of the melter, a thermal resistance approach is used with published correlations used for boundary conditions. The development of robust and efficient CFD models will be instrumental in enabling the WTP to successfully fulfill its mission of safely stabilizing legacy nuclear waste.

12 - MGMT OF RADIOACTIVE AND NON-RADIOACTIVE WASTE↗
DOE OSTI2023

A general non-Fourier Stefan problem formulation that accounts for memory effects

The Stefan problem is the classical model of a melting phase change. In heterogeneous systems, such phase changes can exhibit non-Fourier (anomalous) behaviors, where the advance of the melt interface does not follow the expected time scaling. These situations can be modeled by replacing the derivatives, in the governing partial differential equations, with fractional order derivatives. In particular, replacing the time derivatives leads to non-Fourier models that account for memory effects in the system. In this work, by using appropriate time convolution integrals, a general thermodynamic balance statement for melting phase problems, explicitly accounting for memory effects, is developed. From this balance, a general model formulation applicable to problems involving melting over a temperature range (i.e., a mushy region) is derived. A key component in this model is the representation of memory effects through the use of fractional derivative based constitutive models of the enthalpy and heat flux. Further, on shrinking the mushy region to a single isotherm, a general sharp interface melting model is obtained. Here, in contrast to the classic Stefan problem, the fractional derivatives induce a natural regularization, such that the constitutive models for enthalpy and heat flux are continuous at the melt interface; a result confirmed through numerical simulation. To further support the theoretical findings, a physical example of a non-Fourier Stefan problem is presented. Overall the development and results in this paper underscore the importance of explicitly relating the development of fractional calculus models to the appropriate thermodynamic balance statements.

42 ENGINEERING↗
NASA NTRS1993

A numerical treatment of geodynamic viscous flow problems involving the advection of material interfaces

Effective numerical treatment of multicomponent viscous flow problems involving the advection of sharp interfaces between materials of differing physical properties requires correction techniques to prevent spurious diffusion and dispersion. We develop a particular algorithm, based on modern shock-capture techniques, employing a two-step nonlinear method. The first step involves the global application of a high-order upwind scheme to a hyperbolic advection equation used to model the distribution of distinct material components in a flow field. The second step is corrective and involves the application of a global filter designed to remove dispersion errors that result from the advection of discontinuities (e.g., material interfaces) by high-order, minimally dissipative schemes. The filter introduces no additional diffusion error. Nonuniform viscosity across a material interface is allowed for by the implementation of a compositionally weighted-inverse interface viscosity scheme. The combined method approaches the optimal accuracy of modern shock-capture techniques with a minimal increase in computational time and memory. A key advantage of this method is its simplicity to incorporate into preexisting codes be they finite difference, element, or volume of two or three dimensions.

Lenardic, A.↗
DOE OSTI2026

Quadratic magnetoresistance across the insulator–metal transition in V4O7 thin films

We present a comparative magnetotransport study of V4O7 thin films grown simultaneously on single-crystal Al2O3 (r-cut) and amorphous SiO2 substrates, to examine substrate-induced effects on electronic transport across the insulator–metal transition (IMT). Temperature-dependent resistivity reveals a broad electronic ordering regime below the IMT, with both films exhibiting similar transition temperatures but differing by nearly an order of magnitude in absolute resistivity. Despite this strong substrate dependence of the resistivity scale, the magnetoresistance remains predominantly quadratic in magnetic field over wide temperature intervals for both substrates. Analysis of the quadratic magnetoresistance defines an effective mobility scale that evolves systematically across the IMT and exhibits a similar temperature dependence for the two films. Structural characterization, including atomic-resolution scanning transmission electron microscopy of the V4O7/Al2O3 interface and azimuthal X-ray diffraction measurements, reveals a sharp interface and well-defined in-plane crystallographic registry for the r-cut film, in contrast to the rotationally averaged in-plane structure expected for films grown on amorphous SiO2. These results indicate that while substrate-controlled structural registry modifies the resistivity scale, the quadratic magnetotransport response across the IMT remains largely robust, providing a consistent diagnostic of electronic transport in V4O7 thin films.

25 ENERGY STORAGE↗
DOE OSTI2022

A Phase-Field Model for In-Space Manufacturing of Binary Alloys

The integrity of the final printed components is mostly dictated by the adhesion between the particles and phases that form upon solidification, which is a major problem in printing metallic parts using available In-Space Manufacturing (ISM) technologies based on the Fused Deposition Modeling (FDM) methodology. Understanding the melting/solidification process helps increase particle adherence and allows to produce components with greater mechanical integrity. We developed a phase-field model of solidification for binary alloys. The phase-field approach is unique in capturing the microstructure with computationally tractable costs. The developed phase-field model of solidification of binary alloys satisfies the stability conditions at all temperatures. The suggested model is tuned for Ni-Cu alloy feedstocks. We derived the Ginzburg-Landau equations governing the phase transformation kinetics and solved them analytically for the dilute solution. We calculated the concentration profile as a function of interface velocity for a one-dimensional steady-state diffuse interface neglecting elasticity and obtained the partition coefficient, k, as a function of interface velocity. Numerical simulations for the diluted solution are used to study the interface velocity as a function of undercooling for the classic sharp interface model, partitionless solidification, and thin interface.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Extending Interfaces in 3D to Achieve Superior Nanoscale Strength in Ti/Nb Nanolaminates

Tuning the atomic-level structure of nanolaminates enables high strength, increased deformability, and the ability to absorb and mitigate damage due to varied crystalline defects, including dislocations. Here, we present the enhanced strength of Ti/Nb nanolaminates containing thick 3D interfaces (3DIs), relative to their chemically abrupt 2D counterparts. We examine the effects of crystallographic alignment and compositional gradients on mechanical behavior via experiments and phase-field-dislocation dynamics (PFDD) modeling. Mechanical testing reveals that nanolaminates containing thicker 3DIs demonstrate a 28% hardness enhancement compared with sharp-interface nanolaminates. PFDD modeling shows that the critical resolved shear stress (CRSS) increases with the 3DI thickness. Gradual compositional transitions in 3DIs were confirmed via scanning transmission electron microscopy and high-resolution transmission electron microscopy, showing sharp crystallographic transitions and a heightened interface topography. The findings establish a positive function between the 3DI thickness and mechanical robustness for hexagonal-closest-packed-containing composites, emphasizing the role of defect–interface interactions in tailoring the mechanical performance and providing a foundation for future interfacial engineering.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Simulating hindered grain boundary diffusion using the smoothed boundary method

Abstract Grain boundaries can greatly affect the transport properties of polycrystalline materials, particularly when the grain size approaches the nanoscale. While grain boundaries often enhance diffusion by providing a fast pathway for chemical transport, some material systems, such as those of solid oxide fuel cells and battery cathode particles, exhibit the opposite behavior, where grain boundaries act to hinder diffusion. To facilitate the study of systems with hindered grain boundary diffusion, we propose a model that utilizes the smoothed boundary method to simulate the dynamic concentration evolution in polycrystalline systems. The model employs domain parameters with diffuse interfaces to describe the grains, thereby enabling solutions with explicit consideration of their complex geometries. The intrinsic error arising from the diffuse interface approach employed in our proposed model is explored by comparing the results against a sharp interface model for a variety of parameter sets. Finally, two case studies are considered to demonstrate potential applications of the model. First, a nanocrystalline yttria-stabilized zirconia solid oxide fuel cell system is investigated, and the effective diffusivities are extracted from the simulation results and are compared to the values obtained through mean-field approximations. Second, the concentration evolution during lithiation of a polycrystalline battery cathode particle is simulated to demonstrate the method’s capability.

Materials Science↗
DOE OSTI2022

Ephemeral antibubbles: Spatiotemporal evolution from direct numerical simulations

Antibubbles, which consist of a shell of a low-density fluid inside a high-density fluid, have several promising applications. We show, via extensive direct numerical simulations (DNSs), in both two and three dimensions, that the spatiotemporal evolution of antibubbles can be described naturally by the coupled Cahn-Hilliard-Navier-Stokes (CHNS) equations for a binary fluid. Our DNSs capture elegantly the gravity-induced thinning and breakup of an antibubble via the time evolution of the Cahn-Hilliard scalar-order-parameter field $\phi$, which varies continuously across interfaces, so we do not have to enforce complicated boundary conditions at the moving antibubble interfaces. To ensure that our results are robust, we supplement our CHNS simulations with sharp-interface volume-of-fluid DNSs. We track the thickness of the antibubble and calculate the dependence of the lifetime of an antibubble on several parameters; we show that our DNS results agree with various experimental results; in particular, the velocity with which the arms of the antibubble retract after breakup scales as σ 1/2 , where σ is the surface tension.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Localized Excitons in NbSe 2 -MoSe 2 Heterostructures

Neutral and charged excitons (trions) in atomically thin materials offer important capabilities for photonics, from ultrafast photodetectors to highly efficient light-emitting diodes and lasers. Recent studies of van der Waals (vdW) heterostructures comprised of dissimilar monolayer materials have uncovered a wealth of optical phenomena that are predominantly governed by interlayer interactions. Here, we examine the optical properties in NbSe 2 -MoSe 2 vdW hetero- structures, which provide an important model system to study metal–semiconductor interfaces, a common element in optoelectronics. Through low-temperature photoluminescence (PL) microscopy, we discover a sharp emission feature, L1, that is localized at the NbSe 2 -capped regions of MoSe 2 . L1 is observed at energies below the commonly studied MoSe 2 excitons and trions and exhibits temperature- and power-dependent PL consistent with exciton localization in a confining potential. This PL feature is robust, observed in a variety of samples fabricated with different stacking geometries and cleaning procedures. Using first-principles calculations, we reveal that the confinement potential required for exciton localization naturally arises from the in-plane band bending due to the changes in the electron affinity between pristine MoSe 2 and NbSe 2 -MoSe 2 heterostructure. Here, we discuss the implications of our studies for atomically thin optoelectronics devices with atomically sharp interfaces and tunable electronic structures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Dislocation dynamics in heterogeneous nanostructured materials

Crystalline materials can be strengthened by introducing dissimilar phases that impede dislocation glide. At the same time, the changes in microstructure and chemistry usually make the materials less ductile. One way to circumvent the strength–ductility dilemma is to take advantage of heterogeneous nanophases which simultaneously serve as dislocation barriers and sources. Owing to their superior mechanical properties, heterogeneous nanostructured materials (HNMs) have attracted a lot of attention worldwide. Nevertheless, it has been difficult to characterize dislocation dynamics in HNMs using classical continuum models, mainly due to the challenges in describing the elastic and plastic heterogeneity among the phases. Here, in this work, we advance a phase-field dislocation dynamics (PFDD) model to treat multi-phase materials, consisting of phases differing in composition, structural order, and size in the same system. We then apply the advanced PFDD model to exploring two important but divergent materials design problems in HNMs: dislocation/obstacle interactions and dislocation/interface interactions. Results show that the interactions between a dislocation and distribution of obstacles varying in structure and composition cannot be understood by simply interpolating from their individual interactions with a dislocation. It is also found that materials containing interfaces with nanoscale thicknesses and compositional gradients have a much higher dislocation bypass stress than those with sharp interfaces, providing an explanation for the simultaneous high strength and toughness of thick interface-containing nanolaminates as observed in recent experiments.

36 MATERIALS SCIENCE↗