Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 OSTI2024

A mineral precipitation model based on the volume of fluid method

A novel volume of fluid method is presented for mineral precipitation coupled with fluid flow and reactive transport. The approach describes the fluid-solid interface as a smooth transitional region, which is designed to provide the same precipitation rate and viscous drag force as a sharp interface. Specifically, the governing equation of mineral precipitation is discretized by an upwind scheme, and a rigorous effective viscosity model is derived around the interface. The model is validated against analytical solutions for mineral precipitation in channel and ring-shaped structures. It also compares well with interface tracking simulations of advection-diffusion-reaction problems. Here, the methodology is finally employed to model mineral precipitation in fracture networks, which is challenging due to the low porosity and complex geometry. Compared to other approaches, the proposed model has a concise algorithm and contains no free parameters. In the modeling, only the pore space requires meshing, which improves the computational effciency especially for low-porosity media.

58 GEOSCIENCES↗
DOE OSTI2018

Effects of fuel viscosity on the primary breakup dynamics of a high-speed liquid jet with comparison to X-ray radiography

One of the major concerns in combustion engines is the sensitivity of engine performance to fuel prop- erties. Recent works have shown that even slight differences in fuel properties can cause significant changes in performance and emission of an engine. In order to design the combustion engines with multi-fuel flexi- bilities, the precise assessment of fuel sensitivity on liquid jet atomization process is a prerequisite since the resulting fuel/air mixture is critical to the subsequent combustion process. The present study is focusing on the effect of physical fuel properties, mostly viscosity difference, on the breakup process of the liquid jet injected into still air. Two different jet fuels, CAT-A2 and CAT-C3, are considered here as surrogates for a fossil-based fuel and a bio-derived high-viscosity alternative fuel. The simulations are performed using the volume-of-fluid (VoF) interface tracking method coupled to Lagrangian particle method in order to capture the breakup instabilities of jets and the resulting droplets. The investigations take the actual geometry of the injector into account to resolve the unsteady flow phenomena inside the nozzle that impact the turbu- lence transition and atomization. The simulation results are compared to the experimental measurement using X-ray radiography. Both simulation and X-ray measurements consistently describe the effects of different fuels on the fundamental properties of atomization including the breakup length, transverse liquid volume fraction and the droplet sauter-mean-diameter. The application of a Detailed Numerical Simulation approach complemented by unique X-ray diagnostics is novel and providing new understanding and research directions in engine spray dynamics

Breakup instability↗
DOE OSTI2024

Pore‐Scale Modeling of Reactive Transport with Coupled Mineral Dissolution and Precipitation

Abstract We present a new pore‐scale model for multicomponent advective‐diffusive transport with coupled mineral dissolution and precipitation. Both dissolution and precipitation are captured simultaneously by introducing a phase transformation vector field representing the direction and magnitude of the overall phase change. An effective viscosity model is adopted in simulating fluid flow during mineral dissolution‐precipitation that can accurately capture the velocity field without introducing any empirical parameters. The proposed approach is validated against analytical solutions and interface tracking simulations in simplified structures. After validation, the proposed approach is employed in modeling realistic rocks where mineral dissolution and precipitation are dominant at different locations. We have identified three regimes for mineral dissolution‐precipitation coupling: (a) compact dissolution‐precipitation where dissolution is dominant near the inlet and precipitation is dominant near the outlet, (b) wormhole dissolution with clustered precipitation where dissolution generates wormholes in the main flow paths and precipitation clogs the secondary flow paths, and (c) dissolution dominant where all solid grains are gradually dissolved. In the three regimes, the proposed approach provides reliable porosity‐permeability relationships that cannot be described well by traditional macroscale models. We find that the permeability can increase while the overall porosity decreases when the main flow paths are expanded by dissolution and adjacent pore spaces are clogged by precipitation.

58 GEOSCIENCES↗
DOE OSTI2016

A high-fidelity approach towards simulation of pool boiling

A novel numerical approach is developed to simulate the multiscale problem of pool-boiling phase change. The particular focus is to develop a simulation technique that is capable of predicting the heat transfer and hydrodynamic characteristics of nucleate boiling and the transition to critical heat flux on surfaces of arbitrary shape and roughness distribution addressing a critical need to design enhanced boiling heat transfer surfaces. The macro-scale of the phase change and bubble dynamics is addressed through employing off-the-shelf Computational Fluid Dynamics (CFD) methods for interface tracking and interphase mass and energy transfer. The micro-scale of the microlayer, which forms at early stage of bubble nucleation near the wall, is resolved through asymptotic approximation of the thin-film theory which provides a closed-form solution for the distribution of the micro-layer and its influence on the evaporation process. In addition, the sub-grid surface roughness is represented stochastically through probabilistic density functions and its role in bubble nucleation and growth is then represented based on the thermodynamics of nucleation process. This combination of deterministic CFD, local approximation, and stochastic representation allows the simulation of pool boiling on any surface with known roughness and enhancement characteristics. The numerical model is validated for dynamics and hydrothermal characteristics of a single nucleated bubble on a flat surface against available literature data. In addition, the prediction of pool-boiling heat transfer coefficient is verified against experimental measurements as well as reputable correlations for various roughness distributions and different surface orientations. Finally, the model is employed to demonstrate pool-boiling phenomenon on enhanced structures with reentrance cavities and to explore the effect of enhancement feature design on thermal and hydrodynamic characteristics of these surfaces.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Spatially Resolved Recombination Analysis of CuIn $_\text{x}$ Ga $_\text{1-x}$ Se $_{2}$ Absorbers With Alkali Postdeposition Treatments

In this contribution, we probe spatial variations in charge-carrier recombination in CuIn x Ga 1-x Se 2 (CIGS) absorbers grown on soda-lime glass (SLG) and alkali-free sapphire substrates with NaF and KF postdeposition treatments (PDTs). Temperature-and illumination-dependent device measurements are used to track interface recombination and recombination in the quasi-neutral region. The analysis of these data reveals that the benefit of alkali PDTs depends on the substrate: interface recombination is reduced in devices grown on sapphire substrates, whereas recombination in the quasi-neutral regions is reduced in devices grown on SLG substrates. Cathodoluminescence (CL) spectrum imaging is used to study the spatial distribution of recombination with respect to the grain structure. The grain-boundary CL contrast is similar in films with no PDT, NaF PDT, or KF PDT. A reduced grain-boundary contrast is observed with a NaF + KF PDT; however, suggesting a reduced recombination strength at the grain boundaries (GBs) for combined NaF + KF treatment. CL spectra indicate band tailing, consistent with the fluctuating potential model. Fluctuating potentials are believed to reduce open-circuit voltage, but their spatial distribution has not been studied. Here, CL spectrum imaging data are used to generate maps of the root-mean-square value of the potential energy fluctuations-γ. These maps reveal a bimodal ..gamma.. distribution for all samples: ..gamma.. is generally in the range ~15-50 meV or ~100-180 meV. The higher ..gamma.. range is more significantly affected by the PDTs; after the PDTs, it is strongly correlated with GBs. The lower ..gamma.. range is correlated with higher emission intensity regions, typically grain interiors, and increases in area fraction after the PDTs. These results demonstrate how spatially resolved luminescence and device characterization measurements can be used to monitor changes in recombination in CIGS films and photovoltaic devices. Such measurements can complement empirical device optimization and help improve device performance.

14 SOLAR ENERGY↗
DOE OSTI2021

picassodev/picasso

Picasso is a C++ template library for the development of exascale-class particle-in-cell methods. It contains data structures for particle and grid management using the ECP Kokkos and Cabana libraries, abstractions for particle and grid discretizations, and tools for interface tracking and reconstruction with particles.

Slattery, Stuart↗
DOE OSTI2016

Two-phase Flow Regime Transition Study using Level-Set Method

In Boiling Water Reactors (BWRs) the coolant exiting the core often has void fractions as high as 50%, resulting in high gas volume fraction two-phase flow regimes, such as annular or churn-turbulent flow. As the coolant boils along the subchannels it will experience flow regime transitions from bubbly to slug, slug to churn, and, finally, churn to annular. While these transition points are known to occur as a higher void develops along the flow path, the driving physics of the transition phenomenon is not well understood. Currently, experiment based correlations are the most reliable method for predicting where in the core the flow regime transitions will occur. However, experimentally, it's ever difficult to observe the details of the two-phase flow at prototypical pressures/temperatures, or, as with high fidelity experiments, the conditions differ from those in a reactor. Being able to predict these flow regime transitions is critical to safe and efficient reactor operations as well as maintaining fuel rod integrity by avoiding dryout. The characterization of flow regimes has been extensively studied, producing a wide set of flow regime maps based on characteristic parameters such as vapor and liquid superficial velocity, time dependent pressure drop, interfacial area density, etc. While these studies are useful in predicting the flow regime for certain geometry and flow conditions, they do not provide a universal predictive capability for any geometry under any conditions. Additionally, it is difficult to achieve reactor realistic pressure and temperature conditions in experiments. Direct Numerical Simulation (DNS) coupled with Interface Tracking Methods (ITMs) becomes advantageous in this case because the simulations are not limited by the material properties or to the real world conditions. DNS allows for analysis of flows that would be impossible to achieve, or observe in detail, experimentally. The finite element based DNS code, PHASTA (parallel, hierarchic, higher-order accurate, adaptive, stabilized transient analysis flow solver), used in this work utilizes the level set ITM and has been previously shown to properly resolve bubbly, slug, and annular flows, but not transitions between the two-phase flow regimes. The present work uses PHASTA to analyze flow in a pipe for various void fractions and superficial velocities. Flow regime transitions were induced to demonstrate the transient behavior of this phenomenon. With the unique advanced analysis capabilities developed in PHASTA that allow for individual bubble tracking, it becomes possible to extract new localized information about the flow behavior that is currently not attainable from experiments. The potential of this capability is to uncover the underlying physics of the flow regime transition phenomenon. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2026

Reduced-Order CFD Modeling to Support Waste Loading Optimization in Hanford WTP Vitrification

The U.S. DOE Hanford Site stores over 56 million gallons of radioactive liquid tank waste that must be treated and immobilized for long-term disposal The Waste Treatment and Immobilization Plant (WTP) will vitrify this waste by feeding it into Joule-heated melters, where it is incorporated into a stable borosilicate glass Computational fluid dynamics (CFD) simulations of glass melters can provide insight into the maximum achievable waste loading under varying melter operating conditions Fully resolved VOF multiphase simulations were used as the reference model to capture bubble-driven convection in the melter, including bubble formation, rise behavior, and induced glass melt circulation Effective bubble column diameter and rise velocity were extracted from the resolved simulations, compared with empirical correlations, and refit across relevant viscosity and gas flow rate conditions Explicit gas–liquid interface tracking was replaced with a single-phase momentum source term model, enabling faster steady-state CFD simulations while preserving the dominant hydrodynamic effects of bubbling New empirical correlations were developed for effective bubble column diameter and bubble rise velocity by fitting resolved simulation data across expected melter viscosity and gas flow rate ranges, providing improved inputs for the momentum source term model compared with existing literature correlations The momentum source term model reduced fluid-domain mesh size by 89% and achieved an 8.4× computational speedup relative to resolved bubbling simulations The validated momentum source term approach enables prediction of process-relevant heat transfer behavior in the integrated melter model, including heat transfer from the molten glass to the cold cap, plenum, refractory walls, and surrounding structural regions under varying melter operating conditions

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

Sharp front tracking with geometric interface reconstruction

Here, this paper presents a novel sharp front-tracking method designed to address limitations in classical front-tracking approaches, specifically their reliance on smooth interpolation kernels and extended stencils for coupling the front and fluid mesh. In contrast, the proposed method employs exclusively sharp, localized interpolation and spreading kernels, restricting the coupling to the interfacial fluid cells–those containing the interface/front. This localized coupling is achieved by integrating a divergence-preserving velocity interpolation method with a piecewise parabolic interface calculation (PPIC) and a polyhedron intersection algorithm to compute the indicator function and local interface curvature. Surface tension is computed using the Continuum Surface Force (CSF) method, maintaining consistency with the sharp representation. Additionally, we propose an efficient local roughness smoothing implementation to account for surface mesh undulations, which is easily applicable to any triangulated surface mesh. Building on our previous work, the primary innovation of this study lies in the localization of the coupling for both the indicator function and surface tension calculations. By reducing the interface thickness on the fluid mesh to a single cell, as opposed to the 4–5 cell spans typical in classical methods, the proposed sharp front-tracking method achieves a highly localized and accurate representation of the interface. This sharper representation mitigates parasitic currents and improves force balancing, making it particularly suitable for scenarios where the interface plays a critical role, such as microfluidics, fluid-fluid interactions, and fluid-structure interactions. The proposed method is comprehensively validated and tested on canonical interfacial flow problems, including stationary and translating Laplace equilibria, oscillating droplets, and rising bubbles. The presented results demonstrate that the sharp front-tracking method significantly outperforms the classical approach in terms of accuracy, stability, and computational efficiency. Notably, parasitic currents are reduced by approximately two orders of magnitude and stable results are obtained for parameter ranges where classical front tracking fails to converge.

42 ENGINEERING↗
DOE OSTI2024

Role of interfaces on the mechanical response of accumulative roll bonded nanometallic laminates investigated via dislocation dynamics simulations

Unraveling the effects of continuous dislocation interactions with interfaces, particularly at the nanometer length scales, is key to a broader understanding of plasticity, to material design and to material certification. To this end, this work proposes a novel discrete dislocation dynamics-based model for dislocation interface interactions tracking the fate of residual dislocation on interfaces. This new approach is used to predict the impact of dislocation/interface reactions on the overall mechanical behavior of accumulative roll bonded nanometallic laminates. The framework considers the dynamic evolution of the interface concurrent with a large network of dislocations, thus, accounting for the local short and long range effects of the dislocations under the external boundary conditions. Specifically, this study focuses on two-phase Fe/Cu nanometallic laminates, and investigates the role of the underlying elastic and plastic contrast of the Fe and the Cu layers on the composite response of the material. Moreover, the role of initial microstructures, resulting from processing is also investigated. Subsequently, the model is used to examine the effect of layer thickness and interface orientation relationship on the residual stresses of the relaxed microstructure. The associated mechanical response of these laminates are compared when loaded under normal direction compression, as well as shear compression. Finally, this work predicts a dominant effect of the layer thickness, as compared to the interface orientation relationship, on the macroscopic response and on the residual stresses of these nanolaminates, while the local dislocation transmission propensity through the interface is significantly influenced by the corresponding orientation relationship.

36 MATERIALS SCIENCE↗
DOE OSTI2018

A Digitally Interfaced Analog Correlation Filter System for Object Tracking Applications

This study reports a 24×57 correlation filter system for object tracking applications. While digital interfacing of the input and output data enabled a standard and flexible way of communication with pre- and post-processing digital blocks, the multiply-accumulate (MAC) operations were performed in the analog domain to save power and area. The proposed system utilizes non-volatile floating-gate memories to store filter coefficients. The chip was fabricated in a 0.13- μm CMOS process and occupies 3.23 mm 2 of silicon area. The system dissipates 388.4 μW of power at a throughput of 11.3 kVec/s, achieving an energy efficiency of 25.2 pJ/MAC. Finally, experimental results for a custom filter designed to detect vehicles are presented.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

SU-F-T-518: Development and Characterization of a Gated Treatment System Implemented with An In-House Optical Tracking System and the Elekta Response Interface

Purpose: To report the development and characterization of the first in-house gating system implemented with an optical tracking system (OTS) and the Elekta Response™ interface. Methods: The Response™ connects a patient tracking system with a linac, enabling the tracking system to control radiation delivery. The developed system uses an in-house OTS to monitor patient breathing. The OTS consists of two infrared-based cameras, tracking markers affixed on patient. It achieves gated or breath-held (BH) treatment by calling beam ON/OFF functions in the Response™ dynamic-link library (DLL). A 4D motion phantom was used to evaluate its dosimetric and time delay characteristics. Two FF- and two FFF-IMRT beams were delivered in non-gated, BH and gated mode. The sinusoidal gating signal had a 6 sec period and 15 mm amplitude. The duty cycle included 10%, 20%, 30% and 50%. The BH signal was adapted from the sinusoidal wave by inserting 15 sec BHs. Each delivery was measured with a 2D diode array (MapCHECK™) and compared with the non-gated delivery using gamma analysis (3%). The beam ON/OFF time was captured using the service graphing utility of the linac. Results: The gated treatments were successfully delivered except the 10% duty cycle. The BH delivery had perfect agreement (100%) with non-gated delivery; the agreement of gated delivery decreased from 99% to 88% as duty cycle reduced from 50% to 20%. The beam on/off delay was on average 0.25/0.06 sec. The delivery time for the 50%, 30% and 20% duty cycle increased by 29%, 71% and 139%, respectively. No dosimetric or time delay difference was noticed between FF- and FFF-IMRT beams. Conclusion: The in-house gating system was successfully developed with dosimetric and time delay characteristics in line with published results for commercial systems. It will be an important platform for further research and clinical development of gated treatment.

43 PARTICLE ACCELERATORS↗
DOE OSTI2018

Numerical Simulation of Multiphase Magnetohydrodynamic Flow and Deformation of Electrolyte–Metal Interface in Aluminum Electrolysis Cells

A computational fluid dynamics based multiphase magnetohydrodynamic (MHD) flow model for simulating the melt flow and bath–metal interface deformation in realistic aluminum reduction cells is presented. The model accounts for the complex physics of the MHD problem in aluminum reduction cells by coupling two immiscible fluids, electromagnetic field, Lorentz force, flow turbulence, and complex cell geometry with large length scale. Especially, the deformation of bath–metal interface is tracked directly in the simulation, and the condition of constant anode–cathode distance (ACD) is maintained by moving anode bottom dynamically with the deforming bath–metal interface. The metal pad deformation and melt flow predicted by the current model are compared to the predictions using a simplified model where the bath–metal interface is assumed flat. The effects of the induced electric current due to fluid flow and the magnetic field due to the interior cell current on the metal pad deformation and melt flow are investigated. The presented model extends the conventional simplified box model by including detailed cell geometry such as the ledge profile and all channels (side, central, and cross-channels). The simulations show the model sensitivity to different side ledge profiles and the cross-channel width by comparing the predicted melt flow and metal pad heaving. In addition, the model dependencies upon the reduction cell operation conditions such as ACD, current distribution on cathode surface and open/closed channel top, are discussed.

42 ENGINEERING↗
DOE OSTI2023

Tracking Molecular Transport Across Oil/Aqueous Interfaces: Insight into “Antagonistic” Binding in Solvent Extraction

Liquid/liquid (L/L) interfaces play a key, yet poorly understood, role in a range of complex chemical phenomena where time-evolving interfacial structures and transient supramolecular assemblies act as gatekeepers to function. Here, for this study, we employ surface-specific vibrational sum frequency generation combined with neutron and X-ray scattering methods to track the transport of dioctyl phosphoric acid (DOP) and di-(2-ethylhexyl) phosphoric acid (DEHPA) ligands used in solvent extraction at buried oil/aqueous interfaces away from equilibrium. Our results show evidence for a dynamic interfacial restructuring at low ligand concentrations in contrast to expectation. These time-varying interfaces arise from the transport of sparingly soluble interfacial ligands into the neighboring aqueous phase. These results support a proposed “antagonistic” role of ligand complexation in the aqueous phase that could serve as a holdback mechanism in kinetic liquid extractions. These findings provide new insights into interfacially controlled chemical transport at L/L interfaces and how these interfaces vary chemically, structurally, and temporally in a concentration-dependent manner and present potential avenues to design selective kinetic separations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗