Search NASA⌕ Search

SEARCH · Search NASA

Results for “Interface problems”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Numerical simulation of phase transition problems with explicit interface tracking

Phase change is ubiquitous in nature and industrial processes. Started from the Stefan problem, it is a topic with a long history in applied mathematics and sciences and continues to generate outstanding mathematical problems. For instance, the explicit tracking of the Gibbs dividing surface between phases is still a grand challenge. Our work has been motivated by such challenge and here we report on progress made in solving the governing equations of continuum transport in the presence of a moving interface by the front tracking method. The most pressing issue is the accounting of topological changes suffered by the interface between phases wherein break up and/or merge takes place. The underlying physics of topological changes require the incorporation of space-time subscales not at reach at the moment. Therefore we use heuristic geometrical arguments to reconnect phases in space. This heuristic approach provides new insight in various applications and it is extensible to include subscale physics and chemistry in the future. We demonstrate the method on applications such as simulating freezing, melting, dissolution, and precipitation. The later examples also include the coupling of the phase transition solution with the Navier-Stokes equations for the effect of flow convection.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Scattering of surface modes propagating along metamaterial interface

The problem of surface mode scattering from an iris diaphragm is studied by the integral equation technique. The theory is illustrated by the example of scattering of mode, propagating along metamaterial interface. The characteristics of radiation modes excited by the diaphragm are considered. The surface mode reflection and transmission coefficients, the near-field distribution and far-field radiation pattern are calculated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Nonoverlapping Localized Exponential Time Differencing Methods for Diffusion Problems

In this paper, we propose nonoverlapping localized exponential time differencing (ETD) methods for diffusion problems. The model time-dependent diffusion equation is first reformulated on subdomains based on the nonoverlapping domain decomposition, in which Neumann boundary conditions are imposed on the interfaces for the subdomain problems and Dirichlet type conditions are enforced to form a space-time interface problem. After spatial discretization by standard central finite differences and temporal integration with the first or second order ETD methods, the fully discrete interface problem is obtained. Such an interface problem is then solved iteratively either at each time step or over the whole time interval: the former involves the solution of stationary problems in each subdomain at each iteration while the latter involves the solution of time-dependent subdomain problems at each iteration. For both approaches, we prove that localized ETD solutions conserve mass exactly and converge in time to the exact space semidiscrete solution. Numerical experiments in two dimensions are also presented to illustrate the performance of the proposed methods.

Mass conservation↗
DOE OSTI2018

A gradient-based optimization approach for the detection of partially connected surfaces using vibration tests

The integrity of engineering structures is often compromised by embedded surfaces that result from incomplete bonding during the manufacturing process, or initiation of damage from fatigue or impact processes. Examples include delaminations in composite materials, incomplete weld bonds when joining two components, and internal crack planes that may form when a structure is damaged. In many cases the areas of the structure in question may not be easily accessible, thus precluding the direct assessment of structural integrity. In this paper, we present a gradient-based, partial differential equation (PDE)-constrained optimization approach for solving the inverse problem of interface detection in the context of steady-state dynamics. An objective function is defined that represents the difference between the model predictions of structural response at a set of spatial locations, and the experimentally measured responses. One of the contributions of our work is a novel representation of the design variables using a density field that takes values in the range [0, 1] and raised to an integer exponent that promotes solutions to be near the extrema of the range. The density field is combined with the penalty method for enforcing a zero gap condition and realizing partially bonded surfaces. Furthermore the use of the penalty method with a density field representation leads to objective functions that are continuously differentiable with respect to the unknown parameters, enabling the use of efficient gradient-based optimization algorithms. Numerical examples of delaminated plates are presented to demonstrate the feasibility of the approach.

42 ENGINEERING↗
DOE OSTI2022

A coupling approach for linear elasticity problems with spatially non-coincident discretized interfaces

Here we present a new method for coupled linear elasticity problems whose finite element discretization may lead to spatially non-coincident discretized interfaces. Our approach combines the classical Dirichlet–Neumann coupling formulation with a new set of discretized interface conditions obtained through Taylor series expansions. We show that these conditions ensure linear consistency of the coupled finite element solution. We then formulate an iterative solution method for the coupled discrete system and apply the new coupling approach to two representative settings for which we also provide several numerical illustrations. The first setting is a mesh-tying problem in which both coupled structures have the same Lamé parameters whereas the second setting is an interface problem for which the Lamé parameters in the two coupled structures are different.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Non-conformal interface-cohesive modeling with the shifted boundary method

The accurate simulation of boundary- and interface-dominated problems on complex geometries remains challenging when boundary- or interface-fitted meshes are difficult to generate, particularly for curved boundaries, polycrystalline microstructures, and dense interface networks. The Shifted Boundary Method (SBM) alleviates this meshing burden by shifting the enforcement of boundary conditions from the true boundary to a nearby surrogate boundary and recovering the effect of the true boundary through geometric correction terms, thereby enabling standard finite element spaces on non-boundary-fitted meshes. In this report, we develop a general shiftedboundary and shifted-interface framework within the open-source MOOSE framework. We first present a general SBM implementation for complex geometries on non-boundary-fitted meshes. We then adopt the Shifted Interface Method (SIM) for internal interfaces and develop a unified shifted-interface treatment in which the interface law is enforced on a surrogate interface and the effect of the true interface is recovered through shifted jumps, fluxes, and tractions. This perspective brings scalar thermal-contact and vector-valued cohesive-zone mechanics into a single framework, the latter realized as the Shifted Cohesive Zone Method (SCZM) and coupled with history-dependent constitutive models from NEML2. We further extend the MOOSE mesh infrastructure to support cohesive-zone calculations on distributed meshes. The framework is verified and demonstrated through three progressive studies: Poisson’s equation on a smoothed starshaped domain, a manufactured thermal-contact problem on a non-interface-fitted mesh, and a two-dimensional polycrystalline representative volume element combining crystal plasticity with cohesive grain-boundary interfaces. Across these studies, the shifted formulations reproduce boundary- and interface-fitted reference solutions with high fidelity, indicating that the proposed framework provides an accurate and efficient route to boundary- and interface-dominated simulations on arbitrary geometries without requiring fitted meshes.

Yang, Cheng-Hau↗
DOE OSTI2019

Domain decomposition and multiscale mortar mixed finite element methods for linear elasticity with weak stress symmetry

Two non-overlapping domain decomposition methods are presented for the mixed finite element formulation of linear elasticity with weakly enforced stress symmetry. The methods utilize either displacement or normal stress Lagrange multiplier to impose interface continuity of normal stress or displacement, respectively. By eliminating the interior subdomain variables, the global problem is reduced to an interface problem, which is then solved by an iterative procedure. The condition number of the resulting algebraic interface problem is analyzed for both methods. A multiscale mortar mixed finite element method for the problem of interest on non-matching multiblock grids is also studied. It uses a coarse scale mortar finite element space on the non-matching interfaces to approximate the trace of the displacement and impose weakly the continuity of normal stress. A priori error analysis is performed. It is shown that, with appropriate choice of the mortar space, optimal convergence on the fine scale is obtained for the stress, displacement, and rotation, as well as some superconvergence for the displacement. Computational results are presented in confirmation of the theory of all proposed methods.

Mathematics↗
DOE OSTI2025

Adaptive Interface-PINNs (AdaI-PINNs) for transient diffusion: Applications to forward and inverse problems in heterogeneous media

We model transient diffusion in heterogeneous materials using a novel physics-informed neural networks framework (PINNs) termed Adaptive interface physics-informed neural networks or AdaI-PINNs (Roy et al. arXiv preprint arXiv:2406.04626, 2024). AdaI-PINNs utilize different activation functions with trainable slopes tailored to each material region within the computational domain, allowing for a fully automated and adaptive PINNs approach to model interface problems with strongly and weakly discontinuous solutions. To enhance its performance in highly heterogeneous transient diffusion systems, we prescribe a suite of robust practices, including appropriate non-dimensionalization of equations, a biased sampling method, Glorot initialization, and the hard enforcement of boundary and initial conditions. Here we evaluate the efficacy of the proposed method on several benchmark forward and inverse problems. Comparative studies on one-dimensional and two-dimensional benchmark problems reveal that the modified AdaI-PINNs outperform its unmodified counterpart, achieving root-mean-square errors that are at least two orders of magnitude better in forward problems. For inverse problems, the maximum errors in the approximated diffusion coefficients by modified AdaI-PINNs are four orders of magnitude better than those of the unmodified version. Additionally, modified AdaI-PINNs demonstrate improved stability in problems with large material mismatches.

42 ENGINEERING↗
DOE OSTI2022

Mechanisms of Interface Cleaning in Heterostructures Made from Polymer‐Contaminated Graphene

Abstract Heterostructures obtained from layered assembly of 2D materials such as graphene and hexagonal boron nitride have potential in the development of new electronic devices. Whereas various materials techniques can now produce macroscopic scale graphene, the construction of similar size heterostructures with atomically clean interfaces is still unrealized. A primary barrier has been the inability to remove polymeric residues from the interfaces that arise between layers when fabricating heterostructures. Here, the interface cleaning problem of polymer‐contaminated heterostructures is experimentally studied from an energy viewpoint. With this approach, it is established that the interface cleaning mechanism involves a combination of thermally activated polymer residue mobilization and their mechanical actuation. This framework allows a systematic approach for fabricating record large‐area clean heterostructures from polymer‐contaminated graphene. These heterostructures provide state‐of‐the‐art electronic performance. This study opens new strategies for the scalable production of layered materials heterostructures.

Huang, Zhujun↗
DOE OSTI2019

Density estimation techniques for multiscale coupling of kinetic models of the plasma material interface

In this work we analyze two classes of Density-Estimation techniques which can be used to consistently couple different kinetic models of the plasma-material interface, intended as the region of plasma immediately interacting with the first surface layers of a material wall. In particular, we handle the general problem of interfacing a continuum multi-species Vlasov-Poisson-BGK plasma model to discrete surface erosion models. The continuum model solves for the energy-angle distributions of the particles striking the surface, which are then driving the surface response. A modification to the classical Binary-Collision Approximation (BCA) method is here utilized as a prototype discrete model of the surface, to provide boundary conditions and impurity distributions representative of the material behavior during plasma irradiation. The numerical tests revealed the superior convergence properties of Gaussian Mixture Models over Kernel Density Estimation methods, with Gaussian Mixtures and Epanechnikov-KDEs both being up to two orders of magnitude faster than Gaussian-KDEs. The methodology here presented allows a self-consistent treatment of the plasma-material interface in magnetic fusion devices, including both the near-surface plasma (plasma sheath and presheath) in magnetized conditions, and surface effects such as sputtering, back-scattering, and ion implantation. The same coupling techniques can also be utilized for other discrete material models such as Molecular Dynamics.

42 ENGINEERING↗
DOE OSTI2016

Mechanics of finite cracks in dissimilar anisotropic elastic media considering interfacial elasticity

Interfacial crack fields and singularities in bimaterial interfaces (i.e., grain boundaries or dissimilar materials interfaces) are considered through a general formulation for two-dimensional (2-D) anisotropic elasticity while accounting for the interfacial structure by means of an interfacial elasticity paradigm. The interfacial elasticity formulation introduces boundary conditions that are effectively equivalent to those for a weakly bounded interface. This formalism considers the 2-D crack-tip elastic fields using complex variable techniques. While the consideration of the interfacial elasticity does not affect the order of the singularity, it modifies the oscillatory effects associated with problems involving interface cracks. Constructive or destructive “interferences” are directly affected by the interface structure and its elastic response. Furthermore, this general formulation provides an insight on the physical significance and the obvious coupling between the interface structure and the associated mechanical fields in the vicinity of the crack tip.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Instability of surface electron cyclotron TM-modes influenced by non-monochromatic alternating electric field

The influence of non-monochromaticity of an external alternating electric field on excitation of TM eigenmodes at harmonics of the electron cyclotron frequency is considered here. These TM-modes propagate along the plasma interface in a metal waveguide. An external static constant magnetic field is oriented perpendicularly to the plasma interface. The problem is solved theoretically using the kinetic Vlasov-Boltzmann equation for description of plasma particles motion and the Maxwell equations for description of the electromagnetic mode fields. The external alternating electric field is supposed to be a superposition of two waves, whose amplitudes are different and their frequencies correlate as 2:1. An infinite set of equations for electric field harmonics of these modes is derived with the aid of nonlinear boundary conditions. This set is solved using the wave packet approach consisting of the main harmonic frequency and two nearest satellite temporal harmonics. Analytical studies of the obtained set of equations allow one to find two different regimes of parametric instability, namely, enhancement and suppression of the instability. Numerical analysis of the instability is carried out for the three first electron cyclotron harmonics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Approximate static condensation algorithm for solving multi-material diffusion problems on meshes non-aligned with material interfaces

In this study, we describe a new algorithm for solving multi-material diffusion problem when material interfaces are not aligned with the mesh. In this case interface reconstruction methods are used to construct approximate representation of interfaces between materials. They produce so-called multi-material cells, in which materials are represented by material polygons that contain only one material. The reconstructed interface is not continuous between cells. Finally, we suggest the new method for solving multi-material diffusion problems on such meshes and compare its performance with known homogenization methods.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Groundwater flow and heat transport for systems undergoing freeze-thaw: Intercomparison of numerical simulators for 2D test cases

In high-elevation, boreal and arctic regions, hydrological processes and associated water bodies can be strongly influenced by the distribution of permafrost. Recent field and modeling studies indicate that a fully-coupled multidimensional thermo-hydraulic approach is required to accurately model the evolution of these permafrost-impacted landscapes and groundwater systems. However, the relatively new and complex numerical codes being developed for coupled non-linear freeze-thaw systems require verification. Here in this paper, this issue is addressed by means of an intercomparison of thirteen numerical codes for two-dimensional test cases with several performance metrics (PMs). These codes comprise a wide range of numerical approaches, spatial and temporal discretization strategies, and computational efficiencies. Results suggest that the codes provide robust results for the test cases considered and that minor discrepancies are explained by computational precision. However, larger discrepancies are observed for some PMs resulting from differences in the governing equations, discretization issues, or in the freezing curve used by some codes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Continuous plating/stripping behavior of solid-state lithium metal anode in a 3D ion-conductive framework

Significance Solid-state lithium metal anode possesses great promise owing to its high energy density and improved safety. This work introduces a strategy for constructing a 3D Li metal anode, which is hosted in a solid-state ion-conducting host and shows a safe and dendrite-free plating/stripping behavior. The host is based on a garnet-type ion-conductive framework and a bottom-deposited Cu current collector. The Li anode is plated within the solid garnet framework from the bottom Cu layer and shows a dendrite-free deposition behavior, effectively averting the dendrite penetration issue. Owing to the 3D ion-conductive host, the volume change and interface contact problem of the Li anode have been significantly mitigated, realizing a high-capacity and safe Li metal anode for solid-state high-energy-density batteries.

25 ENERGY STORAGE↗
DOE OSTI2014

Quantification of Wellbore Leakage Risk Using Non-destructive Borehole Logging Techniques

Well integrity is important at all potential CCS locations and may play a crucial role establishing leakage risk in areas where there is a high density of existing wells that could be impacted by the storage operations including depleted petroleum fields where EOR or CCS will occur. To address a need for risk quantification methods that can be directly applied to individual wells using borehole logging tools a study was conducted using data from five wells in Wyoming. The objectives of the study were: Objective 1: Develop methods to establish the baseline flow parameters (porosity and permeability or mobility) from individual measurements of the material properties and defects in a well. Objective 2: Develop a correlation between field flow-property data and cement logs that can be used to establish the flow-properties of well materials and well features using cement mapping tools. Objective 3: Establish a method that uses the flow-property model (Objective 2) to analyze the statistical uncertainties associated with individual well leakage that can provide basis for uncertainty in risk calculations. The project objectives were met through the logging of five wells in Carbon and Natrona County Wyoming to collect data that was used to estimate individual and average well flow properties and model the results using ultrasonic data collected during the logging. Three of the five wells provided data on point and average flow properties for well annuli. Data from the other two wells were used to create models of cement permeability and test whether information collected in one well could be used to characterize another well. The results of the in-situ point measurements were confirmed by the lab measurements sidewall cores collected near the same depths Objective 1 was met using the data collected through logging, testing, and sampling. The methods were developed that can establish baseline flow parameters of wells by both point and average test methods. The methods to estimate the flow properties modeling of point pressure tests, modeling of vertical interference tests, and laboratory measurement of cased-hole sidewall cores The wells were in sufficiently good shape to allow the development of the characterization methods while still having enough defects to study differences in results as they relate to well integrity. Samples and tests analyzed from three of five wells studied in showed the cements were largely intact and had not degraded from exposure native brines. Log results taken in conjunction with the core measurements indicate that interfaces and/or problems with cement placement due to eccentering provide preferential flow paths for fluids, which can increase the effective permeability of the barrier several orders of magnitude above the permeability of intact cement. The results of the maps created using logging tools indicating that the cement condition and bond are generally good identify a need for more research to understand how logs can be used to predicteffective well permeabilities such as those measured by the VITs in this study.

42 ENGINEERING↗