Search NASA⌕ Search

SEARCH · Search NASA

Results for “interface tracking”

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 OSTI2020

Annular Flow Simulation Supported by Iterative In-Memory Mesh Adaptation

Various flow regimes exist in a boiling water reactor (BWR) as the steam quality increases in the uprising coolant flow, from bubbly flow, slug/churn flow, to annular flow. The annular flow is characterized by the presence of a fast-moving gas core and the surrounding liquid film flowing on the conduit wall. Additionally, entrained droplets can be observed in the gas core with ingested bubbles in the liquid film. The dynamics occurring on the wavy interface between the liquid film and gas core plays a crucial role in affecting the heat transfer rate and pressure drop within the BWR core. However, a fundamental understanding of annular flow is still lacking, partly due to the difficulty in obtaining detailed local data in annular flow experiments. In the current study, a novel simulation framework is developed for the annular flow by coupling a computational fluid dynamics flow solver with state-of-the-art meshing software. The gas-liquid interface is tracked with the level set method. Based on the computed flow solutions, the computational mesh is dynamically adapted in memory to meet the local mesh resolution requirement. This iterative simulation-adaptation framework can ensure the fine mesh resolution across the interface, which not only helps mitigate the mass conservation degradation known to level set methods but also improves the representation of dramatic interface topological changes such as wave breaking and droplet entrainment. The present investigation will shed light onto the complex interfacial processes involved in annular flow and generate much needed simulation data for annular flow modeling.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Determining valence band offsets in heterojunctions using a single core-level x-ray photoelectron spectrum

The behavior and functionality of semiconductor heterojunctions depend critically on the alignments of the valence and conduction bands at the various interfaces. Traditionally band alignment has been measured by x-ray photoelectron spectroscopy using pairs of distinct core levels, one from each side of the interface, to track the valence band maxima in the two materials. Here, we demonstrate that it is possible to determine band alignment using a single core-level photoelectron spectrum for an element that is common to both materials. The energy splitting between the photoemission features originating in the two materials is shown to be dominated by the difference in electrostatic potential across the interface, leading to a reliable determination of the band alignment.

Wang, Le↗
DOE OSTI2018

Numerical Simulation of Slag Eye Formation and Slag Entrapment in a Bottom-Blown Argon-Stirred Ladle

A transient mathematical model is developed for simulating the bubble-steel-slag-top gas four-phase flow in a bottom-blown argon-stirred ladle with a 70-ton capacity. The Lagrangian discrete phase model (DPM) is used for describing the moving behavior of bubbles in the steel and slag. To observe the formation process of slag eye, the volume of fluid (VOF) model is used to track the interfaces between three incompressible phases: metal/slag, metal/gas, and slag/gas. The complex multiphase turbulent flow induced by bubble-liquid interactions is solved by a large eddy simulation (LES) model. Slag eye area and slag droplet dispersion are investigated under different gas flow rates. The results show that the movement of bubbles, formation and collapse of slag eye, volatility of steel/slag interface and behavior of slag entrapment can be properly predicted in the current model. When the gas flow rate is 300 L/min, the circulation driven by the bubble plume will stir the entire ladle adequately and form a slag eye of the right size. At the same time, it will not cause strong erosion to the ladle wall, and the fluctuation of the interface is of adequate intensity, which will be helpful for improving the desulfurization efficiency; the slag entrapment behavior can also be decreased. Interestingly, with the motion of liquid steel circulation, the collision and coalescence of dispersed slag droplets occur during the floating process in the vicinity of the wall.

36 MATERIALS SCIENCE↗
DOE OSTI2019

High Order Anchoring and Reinitialization of Level Set Function for Simulating Interface Motion

A second order interface anchoring method has been developed and used with fast sweeping algorithm for reinitialization of a level set function. The algebraic anchoring formulation ensures that the location of the actual interface is preserved, leading to better mass conservation property. It also provides high order accurate algebraic constraint for solving the Eikonal equation on a finite difference grid. Geometric properties of the interface such as surface normal and curvature are subsequently computed from the reinitialized distance function. Various analytical functions for modeling distortion in level set field are considered and accuracy of reinitialization is evaluated using first and second order anchoring schemes. It is also shown that accurate computation of interface curvature requires a high order anchor in addition to a high order fast sweeping method. Mass conservation property of reinitialization is also analyzed by considering test problems from literature including the classic Rider–Kothe single vortex problem. The formulation is suitable for efficient parallelization for both distributed memory and on-node shared memory parallel systems. In conclusion, scalability and performance of the reinitialization scheme on multiple architectures are demonstrated.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Automatic signal decoding and sensor stability of a 3-electrode mixed-potential sensor for NO x /NH 3 quantification

Sensors to detect mixtures of NO x /NH 3 are needed to monitor emissions of diesel automobiles where a selective catalytic reduction system uses an NH 3 mediated reaction to reduce NO x . We report on the application of a three electrode La 0.8 Sr 0.2 CrO 3 , Au 0.5 Pd 0.5 , Pt mixed potential sensor using yttria-stabilized-zirconia (YSZ) as a solid electrolyte to NO x /NH 3 sensing. Artificial neural networks were used to automatically decode the concentrations of NO x /NH 3 and errors of less than 15% are achieved. The optimal architecture for ANN decoding and the maximum density of training data points are also determined. The stability of the sensor was monitored by electrochemical impedance spectroscopy. The impedance associated with YSZ oxygen ion conduction and the electrochemical reactions at the three-phase interface are tracked for a period of over 100 days.

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

Phase transformation kinetics model for metals

We develop a new model for phase transformation kinetics in metals by generalizing the Levitas–Preston (LP) phase field model of martensite phase transformations (see Levitas and Preston (2002a,b) and Levitas et al. (2003)) to arbitrary pressure. Furthermore, we account for and track: the interface speed of the pressure-driven phase transformation, properties of critical nuclei, as well as nucleation at grain sites and on dislocations and homogeneous nucleation. The volume fraction evolution of each phase is described by employing KJMA (Kolmogorov, 1937; Johnson and Mehl, 1939; Avrami, 1939, 1940, 1941) kinetic theory. We then test our new model for iron under ramp loading conditions and compare our predictions for the α → ϵ iron phase transition to experimental data of Smith et al. (2013). In conclusion, more than one combination of material and model parameters (such as dislocation density and interface speed) led to good agreement of our simulations to the experimental data, thus highlighting the importance of having accurate microstructure data for the sample under consideration.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Exact signed distance fields using parallel Fast Sweeping Method

Signed distance fields are often used in multiphysics simulations to track material interfaces. We present a simple methodology based on the fast sweeping method to generate the exact signed distance from triangular meshes and linear paths on Cartesian grids. The methodology propagates the closest primitive to the boundary to the rest of the domain following the characteristics. A local upwind criterion is used to decide between the new and existing closest primitive at each grid point while capturing the correct sign of the global function. The methodology has optimal computational complexity and runs efficiently in distributed-memory architectures. We include 2D and 3D test cases along with a resolution study up to 0.512 trillion zones and 1,000 computer cores. The solution strategy can also be applied to other types of meshes or collections of primitives.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Extended Electronic Structure Inhomogeneity Created by Double Chain Layer Defects Surrounding Columnar Tracks in Heavy-Ion Irradiated YBa{sub 2}Cu{sub 3}O{sub 7-{delta}}.

In YBa2Cu3O7-delta (YBCO), heavy-ion irradiation creates continuous amorphous tracks that are highly effective for vortex pinning. However, the electronic structure landscape of defects and consequently their vortex pinning roles are unclear. Here, we show double chain layer (DCL) defects which intersect the columnar tracks are additionally produced by high energy Pb ion irradiation in YBCO. The DCL defects are similar to 29 nm long, about three times the columnar defect diameter. The electronic structures of the DCL and columnar track defects were determined using atomic-resolution scanning transmission electron microscopy (STEM) and high resolution electron energy loss spectroscopy (EELS). Results show a decrease in the oxygen and copper content at the interface between the track and the YBCO matrix, resulting in interfacial strain. For the DCL defects, the STEM/EELS study revealed a localized electron doped CuO2 plane next to the DCL defect. The DCL defects thus further extend the electronic inhomogeneity into YBCO.

critical current↗
DOE OSTI2022

Development of machine learning framework for interface force closures based on bubble tracking data

Interfacial force closures in the two-fluid model play a critical role for the predictive capabilities of void fraction distribution. However, the practices of interfacial force modeling have long been challenged by the inherent physical complexity of the two-phase flows. The rapidly expanding computational capabilities in the recent years have made high-fidelity data from the interface-captured direct numerical simulation become more available, and hence potential for data-driven interfacial force modeling has prevailed. In this work, we established a data-driven modeling framework integrated to the HZDR multiphase Eulerian-Eulerian framework for computational fluid dynamics simulations. The data-driven framework is verified in a benchmark problem, where a feedforward neural network managed to capture the non-linear mapping between bubble Reynolds number and drag coefficient and reproduce the void distribution resulting from the baseline model in the test case. The second focus is on utilizing the bubble tracking data set to form a closure for the bubble drag in the turbulent bubbly flow, in which the drag coefficient is set to be correlated with the bubble Reynolds number and the Eötvös number. Pseudo-steady state filtering in the Frenet Frame was carried out to obtain the drag coefficient from the turbulent bubbly flow data. The performance of the data-driven drag model is also examined through a case study, where improvement of model’s prediction near-wall is regarded necessary. In conclusion, discussion and further plans of investigation are provided.

42 ENGINEERING↗
DOE OSTI2021

Effect of topology changes on the breakup of a periodic liquid jet

Here, the breakup of a periodic jet is examined computationally, using a front-tracking/finite-volume method, where the interface is represented by connected marker points moving with the fluid, while the governing equations are solved on a fixed grid. Tracking the interface allows control of whether topology changes take place or not. The Reynolds and Capillary numbers are kept relatively low ($Re = 150$ and $Ca = 2$) so most of the flow is well resolved. The effect of topology changes is examined by following the jet until it has mostly disintegrated, for different “coalescence criterion,” based on the thickness of thin films and threads. The evolution of both two-dimensional and fully three-dimensional flows is examined. It is found that although there is a significant difference between the evolution when no breakup takes place and when it does, once breakup takes place the evolution is relatively insensitive to exactly how it is triggered for a range of coalescence criterion, and any differences are mostly confined to the smallest scales.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Modeling of High-Temperature Corrosion of Zirconium Alloys Using the eXtended Finite Element Method (X-FEM)

Oxidation modeling in modern nuclear fuel performance codes is currently limited by the lack of coupling with mechanics, thus preventing proper description of how high-temperature oxidation impacts mechanical properties. This is mostly due to the fact that the finite difference formalism adopted in corrosion models is incompatible with the direct coupling with mechanics in the finite element modeling employed in modern nuclear fuel performance codes. In this study, a physically based zirconium alloy corrosion model called the Coupled-Current Charge Compensation (C4) model, which was initially developed for operating temperature conditions, has been updated to include high-temperature corrosion in order to provide additional critical information (e.g., oxygen concentration profile) under loss-of-coolant accident (LOCA) conditions—information lacking in existing empirical models. The C4 model was implemented in the MOOSE finite-element framework developed at Idaho National Laboratory, enabling it to be used in the BISON nuclear fuel performance code based on the MOOSE framework. To precisely track the different interfaces at a relatively low computational cost, the eXtended Finite Element Method (X-FEM) was applied in MOOSE. The model’s results were compared to those of existing empirical models as well as metallographic analysis of high-temperature oxidized Zircaloy-4 coupons. Oxygen diffusivities in the a and ß phases resulting from this comparison closely agree with those found in the literature. The C4 model implemented with X-FEM in MOOSE now has the capability to accurately predict oxide, oxygen-stabilized a, and prior ß phase layer growth kinetics under isothermal exposure at high temperature (1000–1500°C). Furthermore, in contrast with the empirical models, the C4 model accounts for the finite thickness of the fuel cladding. It can predict the oxygen concentration profile evolution through the whole cladding, enabling evaluation of the remaining ductile thickness—a crucial variable for modeling the mechanical behavior of the fuel cladding under LOCA. Furthermore, this implementation allows direct coupling with mechanics, at a low computing cost, using finite-element-based nuclear fuel performance codes such as BISON.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Tip-Enhanced Raman Chemical and Chemical Reaction Imaging in H 2 O with Sub-3-nm Spatial Resolution

Reproducible chemical and chemical reaction nano-imaging at solid-liquid interfaces remains challenging, particularly when resolutions on the order of a few nm are sought after. In this work, we demonstrate the latter through liquid tip-enhanced Raman (TER) measurements that target gold nanoplates functionalized with 4-mercaptobenzonitrile (MBN). In addition to chemical imaging and local optical field nano-visualization with high spatial resolution, we observe the signatures of 4-mercaptobenzoic acid, which forms as a result of plasmon-induced hydrolysis of MBN. Evidently, the solvent leads to distinct plasmon-induced/enhanced chemical reaction pathways that have not been documented. Here this work shows that such reactions that take place at solid-liquid interfaces can be tracked with a record sub-3 nm spatial resolution via TER spectral nano-imaging in liquids.

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

Operando Raman spectroscopy uncovers hydroxide and CO species enhance ethanol selectivity during pulsed CO 2 electroreduction

Pulsed CO 2 electroreduction (CO 2 RR) has recently emerged as a facile way to in situ tune the product selectivity, in particular toward ethanol, without re-designing the catalytic system. However, in-depth mechanistic understanding requires comprehensive operando time-resolved studies to identify the kinetics and dynamics of the electrocatalytic interface. Here, we track the adsorbates and the catalyst state of pre-reduced Cu 2 O nanocubes (~30 nm) during pulsed CO 2 RR using sub-second time-resolved operando Raman spectroscopy. By screening a variety of product-steering pulse length conditions, we unravel the critical role of co-adsorbed OH and CO on the Cu surface next to the oxidative formation of Cu-O ad or CuO x /(OH) y species, impacting the kinetics of CO adsorption and boosting the ethanol selectivity. However, a too low OH ad coverage following the formation of bulk-like Cu 2 O induces a significant increase in the C 1 selectivity, while a too high OH ad coverage poisons the surface for C-C coupling. Thus, we unveil the importance of co-adsorbed OH on the alcohol formation under CO 2 RR conditions and thereby, pave the way for improved catalyst design and operating conditions.

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

Near-field models and simulations of pellet ablation in tokamaks

Here, detailed numerical studies of the ablation of a single neon pellet in the plasma disruption mitigation parameter space have been performed. Simulations were carried out using FronTier, a hydrodynamic/low magnetic Reynolds number MHD code with explicit tracking of material interfaces. FronTier's physics models resolve the pellet surface ablation and the formation of a dense, cold cloud of ablated material, the deposition of energy from hot plasma electrons, the expansion of the ablation cloud along magnetic field lines, and the radiation losses. A local thermodynamic equilibrium model based on Saha equations has been used to resolve atomic processes in the cloud and Redlich–Kwong corrections to the ideal gas equation of state for cold and dense gases, which have been used near the pellet surface. The FronTier pellet code is a next generation of the code described. It has been validated against the semi-analytic improved neutral gas shielding model in the 1D spherically symmetric approximation. Main results include quantification of the influence of atomic processes and Redlich–Kwong corrections on the pellet ablation in spherically symmetric approximation and verification of analytic scaling laws in a broad range of pellet and plasma parameters. Using axially symmetric MHD simulations, the properties of ablation channels and the reduction of pellet ablation rates in magnetic fields of increasing strength have been studied. While the main emphasis has been given to neon pellets for the plasma disruption mitigation, selected results on deuterium fueling pellets have also been presented.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Rancor Microworld Simulation Environment for Nuclear Process Control

The Rancor Microworld is a simplified nuclear power plant simulator. It represents the process of an integrated pressurized light water reactor core with a secondary Rankine cycle turbine/generator. The reactor core is governed by simplified core physics. Operators have the ability to raise and lower control rods to change the reactivity and consequently generate more or less electricity. The microworld provides a digital human machine interface for operating and monitoring the plant efficiently to maximize profits within the simulation. The simulation supports multi-unit operations analogous to small modular reactor plants and supports faults on various components within each unit to simulate some common training scenarios such as reactor and turbine trip events. The microworld includes experimental modules to log process parameters and both automated and user control inputs and events, administer situation awareness probes, and interface with eye tracking software.

Ulrich, ThomasA.↗