Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “interface capturing”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Voronoi Meshing to Accurately Capture Geological Structure in Subsurface Simulations

Mesh generation lies at the interface of geological modeling and reservoir simulation. Highly skewed or very small grid cells may be necessary to accurately capture the geometry of geological features, but the resulting poorly scaled or small grid cells can have a substantial negative impact on simulator accuracy and speed. One way to minimize numerical errors caused by gridding complex structures is to simulate on high-quality Voronoi meshes, which reduce grid orientation effects in fluid flow. This work presents a complete methodology to create Voronoi simulation grids, model fluid flow in complex geological systems, and visualize the results. A recently developed Voronoi meshing method that can automatically generate provably good unstructured meshes that conform to input surfaces creating closed volumes is used. Initially an analytical benchmark simulation is presented to validate the quality of the meshes and simulation results and demonstrate the superiority of simulation results using Voronoi meshes over flexed-hexahedral meshes on a domain with internal features. Next, meshes are created for test structures representing four of the most common geological features in the subsurface: layering, pinch-out, an interior lens that tapers to zero thickness on all sides and a fault with offset. Two benchmark flow simulations are run for each test structure. Finally, a realistic geological example for CO 2 injection into an anticline is simulated. Three realizations of the Voronoi mesh at the same resolution are generated for the simulations. Each mesh is highly refined near the injection wells and coarse in areas of less interest. These three meshes are used to model the CO 2 plume in the subsurface as it migrates to the top of the structure and then fills downward. Simulations on the meshes with randomly generated elements inside the input volumes each give slightly different fingering patterns for the viscous-unstable buoyant gas flow. The results presented in this work show a promising step towards utilizing fully automated Voronoi meshing for subsurface flow simulations in complex geology.

58 GEOSCIENCES↗
DOE OSTI2023

The Dynamics of Oxygen Ion Exchange in Epitaxial Strontium Cobaltite Bilayers

Abstract The exchange of ions across interfaces is key to the field of iontronics, where the properties of the device can be altered by the local ion concentration. This study investigates a complex oxide system where structural and electronic phase transitions can be driven by changes in the concentration of oxygen ions. In situ coherent X‐ray studies are conducted on epitaxial bilayers of insulating SrCoO 2.5 and metallic SrCoO 3 − δ . The diffusion of oxygen ions across the bilayer is studied with X‐ray photon correlation spectroscopy to capture the dynamical behavior of the interface in reducing and oxidizing environments. The behavior is strongly asymmetric, with much slower dynamics appearing in reducing versus oxidizing environments. According to the correlation functions determined from different points in reciprocal space, this study finds that the dynamics near the center of the SrCoO 2.5 crystal are generally similar to those near the heterointerfaces. The results suggest that the interface is stable and reversible, making SrCoO x a model system for the study of iontronic behavior.

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

Efficient numerical method for evaluating normal and anomalous time-domain equilibrium Green's functions in inhomogeneous systems

In this work we develop EPOCH (equilibrium propagator by orthogonal polynomial chain), a computationally efficient method to calculate the time-dependent equilibrium Green's functions, including the anomalous Green's functions of superconductors, to capture the time evolution in large inhomogeneous systems. The EPOCH method generalizes the Chebyshev wave-packet propagation method from quantum chemistry and efficiently incorporates the Fermi-Dirac statistics that is needed for equilibrium quantum condensed matter systems. The computational cost of EPOCH scales only linearly in the system degrees of freedom, generating an extremely efficient algorithm also for very large systems. We demonstrate the power of the EPOCH method by calculating the time evolution of an excitation near a superconductor–normal metal interface in two and three dimensions, capturing transmission as well as normal and Andreev reflections.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

A level-set immersed boundary method for reactive transport in complex topologies with moving interfaces

A simulation framework based on the level-set and the immersed boundary methods (LS-IBM) has been developed for reactive transport problems in porous media involving a moving solid-fluid interface. The interface movement due to surface reactions is tracked by the level-set method, while the immersed boundary method captures the momentum and mass transport at the interface. The proposed method is capable of accurately modeling transport near evolving boundaries in Cartesian grids. The framework formulation guarantees second order accuracy in space. Since the interface velocity is only defined at the moving boundary, an interface velocity propagation method is also proposed. The method can be applied to other moving interface problems of the “Stefan” type. Here, we validate the proposed LS-IBM both for flow and transport close to an immersed object with reactive boundaries as well as for crystal growth. Lastly, the proposed method provides a powerful tool to model more realistic problems involving moving reactive interfaces in complex domains.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Bonds over Electrons: Proton Coupled Electron Transfer at Solid–Solution Interfaces

This Perspective argues that most redox reactions of materials at an interface with a protic solution involve net proton-coupled electron transfer (PCET) (or other cation-coupled ET). This view contrasts with the traditional electron-transferfocused view of redox reactions at semiconductors, but redox processes at metal surfaces are often described as PCET. Taking a thermodynamic perspective, transfer of an electron is typically accompanied by a stoichiometric proton, much as the chemistry of lithium-ion batteries involves coupled transfers of e – and Li + . The PCET viewpoint implicates the surface–H bond dissociation free energy (BDFE) as the preeminent energetic parameter and its conceptual equivalents, the electrochemical ne – /nH + potential versus the reversible hydrogen electrode (RHE) and the free energy of hydrogenation, ΔG° H . These parameters capture the thermochemistry of PCET at interfaces better than electronic parameters such as Fermi energies, electron chemical potentials, flatband potentials, or band-edge energies. A unified picture of PCET at metal and semiconductor surfaces is presented. Exceptions, limitations, implications, and future directions motivated by this approach are described.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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↗
DOE OSTI2021

A consistent and conservative volume distribution algorithm and its applications to multiphase flows using Phase-Field models

In the present study, the multiphase volume distribution problem, where there can be an arbitrary number of phases, is addressed using a consistent and conservative volume distribution algorithm. The proposed algorithm satisfies the summation constraint, the conservation constraint, and the consistency of reduction. The first application of the volume distribution algorithm is to determine the Lagrange multipliers in multiphase Phase-Field models that enforce the mass conservation, and a multiphase conservative Allen-Cahn model that satisfies the consistency of reduction is developed. A corresponding consistent and conservative numerical scheme is developed for the model. The multiphase conservative Allen-Cahn model has a better ability than the multiphase Cahn-Hilliard model to preserve under-resolved structures. The second application is to develop a numerical procedure, called the boundedness mapping, to map the order parameters, obtained numerically from a multiphase model, into their physical interval, and at the same time to preserve the physical properties of the order parameters. Along with the consistent and conservative schemes for the multiphase Phase-Field models, the numerical solutions of the order parameters are reduction consistent, conservative, and bounded, which are theoretically analyzed and numerically validated. Then, the multiphase Phase-Field models are coupled with the momentum equation by satisfying the consistency of mass conservation and the consistency of mass and momentum transport, thanks to the consistent formulation. Finally, it is demonstrated that the proposed model and scheme converge to the sharp-interface solution and are capable of capturing the complicated multiphase dynamics even when there is a large density and/or viscosity ratio.

42 ENGINEERING↗
DOE OSTI2020

Air-Stable Calcium Cyanamide-Supported Ruthenium Catalyst for Ammonia Synthesis and Decomposition

Efficient ammonia synthesis and decomposition processes under mild conditions are important to meet the expanding demand in major applications of ammonia as the energy carrier and to provide feedstock for chemical industry. In this study, we report that air-stable calcium cyanamide-supported ruthenium (Ru/CaCN 2 ) works as an efficient and stable catalyst for ammonia synthesis and decomposition. Ru/CaCN 2 exhibits greater catalytic performances for both reactions than Ru/Ca 2 N electride and Ru–Cs/MgO. The kinetic analysis for ammonia synthesis suggests that Ru/CaCN 2 exhibits low apparent activation energy and high resistance to hydrogen poisoning, which has characteristics similar to the kinetic parameters of Ru-supported electride catalyst. H 2 -temperature programmed reaction and temperature programmed desorption revealed that Ru promoted the formation of CN 2 vacancies on the CaCN 2 surface which in turn capture hydrogen as H – ions during the reaction. Density functional theory calculations provide insights into how the formation of CN 2 vacancies is promoted by Ru, which leads to the decrease in the work function of the CaCN 2 surface and the hydrogen capture at the Ru-support interface. These results suggest that the high catalytic performance of Ru/CaCN 2 can be attributed to the formation of a quasi-electride structure at the Ru–CaCN 2 interface.

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

Understanding the Effect of Lead Iodide Excess on the Performance of Methylammonium Lead Iodide Perovskite Solar Cells

The presence of unreacted lead iodide in organic-inorganic lead halide perovskite solar cells is widely correlated with an increase in power conversion efficiency. We investigate the mechanism for this increase by identifying the role of surfaces and interfaces present between methylammonium lead iodide perovskite films and excess lead iodide. We show how type I and II band alignments arising under different conditions result in either passivation of surface defects or hole injection. Through first-principles simulations of solid-solid interfaces, we find that lead iodide captures holes from methylammonium lead iodide and modulates the formation of defects in the perovskite, affecting recombination. Using surface-sensitive optical spectroscopy techniques, such as transient reflectance and time-resolved photoluminescence, we show how excess lead iodide affects the diffusion and surface recombination velocity of charge carriers in methylammonium lead iodide films. Our coupled experimental and theoretical results elucidate the role of excess lead iodide in perovskite solar cells.

14 SOLAR ENERGY↗
DOE OSTI2021

In Situ X-ray Diffraction of High Density Polyethylene During Dynamic Compression: Polymer Chain Compression and Phase Transformations

Polymers are widely implemented materials often used for safety and security in our daily lives. However, we do not have a thorough understanding of the phase transformations of these materials when exposed to extreme condition, high strain rate, events. To study polymers in these severe environments, we performed plate impact experiments that impart a one-dimensional shock wave, at high pressures and temperatures, to observe the materials response. Generally, velocimetry is the experimental diagnostic used to capture the movement of high velocity interfaces to then calculate the mechanical variables associated with shock compression. However, when a material undergoes a phase transition, such as between two crystalline phases or melting to decomposition during shock compression, velocimetry does not explicitly indicate what structural change occurred. At the Dynamic Compression Sector at the Advanced Photon Source gas gun driven compression has been paired with x-ray diffraction, providing the ability to study both the shock parameters and structural changes concurrently. We explored high density polyethylene (HDPE) both on its reactants Hugoniot (<10.12 GPa) and within its transition region (24.3 GPa) where HDPE is hypothesized to decompose from reactants to products. With x-ray diffraction we observed both polymer chain compression and crystal structure reorientation from orthorhombic to monoclinic above 6 GPa, as well as either a possible phase change or chemical decomposition at 24.3 GPa. With these experiments we were able to determine the crystalline and amorphous Hugoniot loci at pressures below the polymer transition region.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Customs: Inspect and Report Python Imports (customs) v1.0.0

We have developed a simple, minimally invasive monitoring framework that enables us to capture key metrics from Python user processes. This is part of a larger effort to collect data across our entire data-intensive science workload at NERSC. The approach for Python is inspired by a similar framework developed for Blue Waters at NCSA. Both frameworks leverage standard Python features, sitecustomize and atexit, to capture Python imports of interest and other job data. The customs package decomposes the information capture into separate inspection and reporting interfaces.

Thomas, Rollin↗
DOE OSTI2025

Dissolved CO 2 Modulates the Electrochemical Capacitance on Gold Electrodes

The presence of CO 2 at an electrified interface between an aqueous electrolyte and a metal electrode is the prerequisite for many electrochemical CO 2 capture technologies. To understand the behavior of dissolved CO 2 at an aqueous electrified interface, we characterized the electrochemical interface of planar gold electrodes with cyclic voltammetry, electrochemical impedance spectroscopy (EIS), electrochemical surface plasmon resonance (EC–SPR), and attenuated total reflectance surface-enhanced infrared spectroscopy (ATR–SEIRAS). Under all investigated conditions, we observed a decrease in the electrochemical capacitance upon saturation of the electrolyte with CO 2 , as compared to an electrolyte saturated with Ar. EIS and EC–SPR showed that this capacitance reduction was also potential dependent: it reached a minimum near the point of zero charge and became more significant as the applied potential moved further away from the point of zero charge. Hybrid quantum–classical simulations of the gold/aqueous electrolyte interface indicate that bicarbonate decreases the capacitance and modifies the composition of the electric double layer. In addition to the binding of bicarbonate under positive bias, we propose that molecular CO 2 can be induced by applied potential to concentrate in the diffuse layer of the electric double layer, leading to a reduction in the electrochemical capacitance under both negative and positive bias. Furthermore, this work advances the understanding of non-Faradaic effects of dissolved CO 2 at aqueous electrified interfaces of relevance for electrochemical CO 2 capture.

25 ENERGY STORAGE↗
DOE OSTI2021

X-ray phase contrast and absorption imaging for the quantification of transient cavitation in high-speed nozzle flows

In this work, high-flux synchrotron radiation has been employed in a time-resolved manner to characterise the distinct topology features and dynamics of different cavitation regimes arising in a throttle orifice with an abrupt flow-entry contraction. Radiographs obtained though both X-ray phase- contrast and absorption imaging have been captured at 67,890 frames per second. The flow lied in the turbulent regime (Re=35,500), while moderate (CN=2.0) to well-established (CN=6.0) cavitation conditions were examined encompassing the cloud and vortical cavitation regimes with pertinent transient features, such as cloud-cavity shedding. X-ray Phase-Contrast Imaging (XPCI), exploiting the shift in the X-ray wave phase during interactions with matter, offers sharp-refractive index gradients in the interface region. Hence, it is suitable for capturing fine morphological fluctuations of transient cavitation structures. Nevertheless, the technique cannot provide information on the quantity of vapour within the orifice. Such data have been obtained utilising absorption imaging, where beam attenuation is not associated with scattering and refraction events, and hence can be explicitly correlated with the projected vapour thickness in line-of-sight measurements. A combination of the two methods is proposed, as it has been found is capable of quantifying the vapour content arising in the complex nozzle flow, while also faithfully illustrating the dynamics of the highly-transient cavitation features.

42 ENGINEERING↗
DOE OSTI2025

Understanding the Effects of Inhomogeneities at the Back Interface of CdTe‐Based Solar Cells Using 2D Modeling

One-dimensional modeling cannot capture lateral inhomogeneities in CdTe-based devices. Here, we use 2D modeling to investigate the role of varying energetics at the back interface. We consider improvements in the back interface layer (BIL) through either reducing back surface recombination velocity (BSRV) or decreasing the downward band bending near the back interface. We show that when the BSRV is reduced, but strong downward band bending remains, there is no change in the device performance until the BSRV of 90% of the back interface is improved by the BIL. On the other hand, any coverage with a BIL that improves band bending results in device improvements. We use band bending, back interface recombination current densities, and voltage dependent current flow through the device to understand these improvements. The modeling shows that lateral flow of carriers greatly affects device performance, which is not captured in parallel diode modeling, and demonstrates improved understanding with 2D modeling.

2D modelling↗
DOE OSTI2020

In Vivo Molecular Insights into Syntrophic Geobacter Aggregates

Direct interspecies electron transfer (DIET) has been considered as a novel and highly efficient strategy in both natural anaerobic environments and artificial microbial fuel cells. A syntrophic model consisting of Geobacter metallireducens and Geobacter sulfurreducens was studied in this work. We conducted in vivo molecular mapping of the outer surface of the syntrophic community as the interface of nutrients and energy exchange. System for Analysis at the Liquid Vacuum Interface combined with time-of-flight secondary ion mass spectrometry was employed to capture the molecular distribution of syntrophic Geobacter communities in the living and hydrated state. Principal component analysis with selected peaks revealed that syntrophic Geobacter aggregates were well differentiated from other control samples, including syntrophic planktonic cells, pure cultured planktonic cells, and single population biofilms. Our in vivo imaging indicated that a unique molecular surface was formed. Specifically, aromatic amino acids, phosphatidylethanolamine components, and large water clusters were identified as key components that favored the DIET of syntrophic Geobacter aggregates. Moreover, the molecular changes in depths of the Geobacter aggregates were captured using dynamic depth profiling. Our findings shed new light on the interface components supporting electron transfer in syntrophic communities based on in vivo molecular imaging.

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

Time-resolved electrical potential pump – X-ray photoelectron spectroscopy probe developments for investigating dynamic processes occurring at electrochemical interfaces

Electrode–electrolyte interfaces are of critical importance in several fields, including renewable energy, corrosion, and environmental chemistry. However, investigating these interfaces under operational conditions poses considerable challenges due to the limitations of the instrumentation employed. While recent advancements in in situ and operando techniques have enhanced our comprehension of the steady-state properties of solid-liquid interfaces, the dynamic behaviors of these systems remain inadequately explored. This study introduces a time-resolved X-ray photoelectron spectroscopy (XPS) technique designed to capture transient reaction intermediates and charging dynamics at electrified interfaces. The presented proof-of-principle study demonstrates that electrochemical processes, represented by an equivalent electrical circuit (EEC) model, can be probed and understood using square wave voltage pulses of a potentiostat synchronized to the modified data acquisition of an XPS setup. This method offers a valuable alternative to traditional pump–probe techniques, facilitating the investigation of a broader range of electrochemical systems. A dedicated software package for analyzing time- and energy-resolved XPS with a focus on extracting parameters of the EEC is geared towards benchmarking different EECs in future real-world electrochemical experiments.

Electrochemistry↗
DOE OSTI2019

A space-time adaptive finite element method with exponential time integrator for the phase field model of pitting corrosion

In this paper we propose a space-time adaptive finite element method for the phase field model of pitting corrosion, which is a parabolic partial differential equation system consisting of a phase variable and a concentration variable. A major challenge in solving this phase field model is that the problem is very stiff, which makes the time step size extremely small for standard temporal discretizations. Another difficulty is that a high spatial resolution is required to capture the steep gradients within the diffused interface, which results in very large number of degrees of freedom for uniform meshes. To overcome the stiffness of this model, we combine the Rosenbrock–Euler exponential integrator with Crank–Nicolson scheme for the temporal discretization. Moreover, by exploiting the fact that the speed of the corroding interface decreases with time, we derive an adaptive time stepping formula. For the spatial approximation, we propose a simple and efficient strategy to generate adaptive meshes that reduces the computational cost significantly. Thus, the proposed method utilizes local adaptivity and mesh refinement for efficient simulation of the corrosive dissolution over long times in heterogeneous media with complex microstructures. We also present an extensive set of numerical experiments in both two and three dimensional spaces to demonstrate efficiency and robustness of the proposed method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Development and testing of coarse-grained models for ultrasonic simulations of cast austenitic stainless steel

Ultrasonic inspection of cast austenitic stainless steel (CASS) in the nuclear industry is particularly challenging because of sound field scatter and attenuation caused by the coarse-grained microstructure. Modeling and simulation are important tools in ultrasonic testing, as they can be used to help address key aspects of inspections, such as developing new probe designs, predicting inspection reliability, and testing phased-array focal laws. However, developing a useful and reliable CASS model is challenging due to the many grain interfaces and crystalline orientations that must be captured. We demonstrate a method of creating a realistic CASS model that is usable in CIVA, a commercially available modeling and simulation software platform. Using polished and chemically etched sections, we generate models of a coarse-grained equiaxed specimen and a columnar specimen. We also test an alternative method of generating a coarse-grained model using Voronoi regions. We qualitatively compare sound field scatter and quantitatively compare sound field attenuation and beam partitioning in simulated sound fields to those of laboratory-measured sound fields. Results show that the Voronoi models perform as well as or better than the models based on actual grain morphology. Here we also show that model-to-model randomness in Voronoi grain structure can impact the magnitude of a simulated echo response by a factor of two or more. Although CASS models are potentially a good depiction of reality for a given scenario, they should not be considered representative since CASS morphology can change significantly from specimen to specimen or within the same specimen.

36 MATERIALS SCIENCE↗