Multi-Material Interface Tracking for Enriched Finite Element Methods.
Abstract not provided.
Fisher, Koa↗
SEARCH · Search NASA
Results for “Interface Tracking Method”
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.
A Thermo‐Flow‐Mechanics‐Fracture Model Coupling a Phase‐Field Interface Approach and Thermo‐Fluid‐Structure Interaction
This work proposes a novel approach for coupling non-isothermal fluid dynamics with fracture mechanics to capture thermal effects within fluid-filled fractures accurately. This method addresses critical aspects of calculating fracture width in enhanced geothermal systems, where the temperature effects of fractures are crucial. The proposed algorithm features an iterative coupling between an interface-capturing phase-field fracture method and interface-tracking thermo-fluid-structure interaction using arbitrary Lagrangian–Eulerian coordinates. We use a phase-field approach to represent fractures and reconstruct the geometry to frame a thermo-fluid-structure interaction problem, resulting in pressure and temperature fields that drive fracture propagation. We developed a novel phase-field interface model accounting for thermal effects, enabling the coupling of quantities specific to the fluid-filled fracture with the phase-field model through the interface between the fracture and the intact solid domain. We provide several numerical examples to demonstrate the capabilities of the proposed algorithm. In particular, we analyze mesh convergence of our phase-field interface model, investigate the effects of temperature on crack width and volume in a static regime, and highlight the method's potential for modeling slowly propagating fractures.
fracture↗
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↗
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↗
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↗
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↗
A computational model for molten corium spreading and solidification
Not Available
Kucala, Alec (ORCID:0000000160502166)↗
Nested Dissection Interface Reconstruction in Pececillo
A nested dissection method for interface reconstruction in a volume tracking framework has been implemented in Pececillo. This method provides a significant improvement over the traditional onion-skin method, which does not appropriately handle T-shaped multimaterial intersections and dynamic contact lines present in additive manufacturing simulations. The resulting implementation lays the groundwork for further re- search in numerical contact angle estimates.
97 MATHEMATICS AND COMPUTING↗
Numerical comparison of bubbling in a waste glass melter
Radioactive tank waste is scheduled for vitrification at the Hanford Tank Waste Treatment and Immobilization Plant, also known as the Vit Plant, being constructed at the Hanford Site. Testing of the pilot-scale, joule-heated DuraMelter 1200 at the Vitreous State Laboratory has shown that bubbling increases the melt rate, and as a result, melter throughput. Computational fluid dynamics (CFD) models of this pilot-scale waste glass melter were developed to improve our understanding of the processes that occur within the melter to aid in process optimization and troubleshooting of the Vit Plant melters. Unfortunately, model validation is complicated by the difficulty of obtaining suitable experimental data for operational melters attributed to the large-scale, radioactivity and high temperatures at which these melters are operated. This study focuses on confirming the fidelity of the CFD models to accurately reproduce the bubbling behavior. Because of the paucity of experimental data at the resolution required for CFD validation, a code-to-code comparison was used to evaluate two common approaches for simulating flows of immiscible fluids on numerical grids and resolving multiphase interfaces. Here, the volume of fluid and level set methods are used to track the dynamically evolving interfaces between the molten glass and the air bubbles. To aid in the verification of the results of these codes, a comparison of the bubble behavior, growth, and frequency of bubble generation are presented and a grid convergence study is performed for the two approaches. In addition to comparing the numerical results to available experimental data and observations from pilot-scale testing, the results are evaluated in relation to laboratory experiments with bubbles injected into viscous fluids, empirical correlations obtained from the published literature, and the basic laws of fluid dynamics. Furthermore, the results of this work can assist in validating waste glass melter models and provide a better understanding of the flow patterns within the WTP melters.
12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
Deformation of contact surfaces in a vacuum interrupter after high-current interruptions
In a high-current interruption, the contact surface in a vacuum interrupter might be severely damaged by constricted vacuum arcs causing a molten area on it. As a result, a protrusion will be initiated by a transient recovery voltage after current zero, enhancing the local electric field and making breakdowns occur easier. The objective of this paper is to simulate the deformation process on the molten area under a high electric field by adopting the finite element method. A time-dependent Electrohydrodynamic model was established, and the liquid-gas interface was tracked by the level-set method. From the results, the liquid metal can be deformed to a Taylor cone if the applied electric field is above a critical value. This value is correlated to the initial geometry of the liquid metal, which increases as the size of the liquid metal decreases. Moreover, the buildup time of a Taylor cone obeys the power law t = k × E{sup −3}, where E is the initial electric field and k is a coefficient related to the material property, indicating a temporal self-similar characteristic. In addition, the influence of temperature has little impact on the deformation but has great impact on electron emission. Finally, the possible reason to initiate a delayed breakdown is associated with the deformation. The breakdown does not occur immediately when the voltage is just applied upon the gap but is postponed to several milliseconds later when the tip is formed on the liquid metal.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Nested Dissection Interface Reconstruction in Pececillo
A nested dissection method for interface reconstruction in a volume tracking framework has been implemented in Pececillo, a mini-app for Truchas, which is the ASC code for casting and additive manufacturing. This method provides a significant improvement over the traditional onion-skin method, which does not appropriately handle T-shaped multimaterial intersections and dynamic contact lines present in additive manufacturing simulations. The resulting implementation lays the groundwork for further research in contact angle estimates and surface tension calculations.
97 MATHEMATICS AND COMPUTING↗
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↗
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↗
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↗
A plasma–vacuum interface tracking algorithm for magnetohydrodynamic simulations of coaxial plasma accelerators
The resistive Magneto-hydrodynamic (MHD) model describes the behavior of a strongly ionized plasma in the presence of external electric and magnetic fields. For problems involving the plasma with an open boundary to very low pressure/vacuum regions, the continuum assumption is no longer valid in the entire domain of interest. For example, this is seen in a plasma accelerator where a dense plasma expands into a vacuum background. A common practice to deal with this issue is to assign small values of density and pressure to the vacuum regions and proceed to solve the continuum based MHD equations throughout the domain. We show that this approach fails to produce solutions consistent with the physics of the problem and can give rise to unacceptable artifacts such as spurious shocks. Here, we develop a plasma–vacuum interface tracking approach to mitigate this problem. The plasma–vacuum interface is tracked in a physically consistent manner and the MHD equations are solved only in the regions that contain the plasma. The interface tracking is achieved using a face-flux formulation derived from the theoretical solution to a 1D free expansion problem. Coupled with a threshold based approach, the interface tracking is implemented for both explicit and implicit time stepping frameworks on generalized unstructured grids. Simulations of magnetized thermal plasma jets expanding into a vacuum background indicate plume profiles devoid of the unphysical shock obtained using the small background density approach. In the context of resistive MHD simulations, the interface tracking approach overcomes the numerical stiffness induced by specifying a large background resistivity in the vacuum regions, resulting in significant wall-clock time gains.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Particle Tracking Methods for Battery Precipitation Reactions
Precipitation and deposition reactions at solid–liquid interfaces play a key role in a number of battery chemistries, including Li-ion, so-called “anode free” batteries, zinc-based battery chemistries, and lithium–sulfur, among others. Although models with heterogeneous nucleation and growth phenomena are present in the literature, papers have not to date provided much detail on the numerical algorithms used to track the temporal evolution of the particle size distribution of deposits on electrode surfaces. In this paper we examine several approaches to discretize and track the particle size distribution, demonstrating that common approaches lead to anomalous flattening of the particle size distribution. We conclude by presenting an algorithm that preserves the appropriate particle size distribution during particle growth.
Algorithms↗
Implementation of a High-Fidelity Interface Resolving Method in Nek5000
The development and utilization of computational fluid dynamics (CFD) models for large, high-temperature electric melters in the Waste Treatment and Immobilization Plant (WTP) in eastern Washington State have proven to be valuable for various purposes. These models allow for a better understanding of the physio-chemical processes occurring within the melter vessels and can contribute to improving operational efficiency, throughput, and addressing operational issues related to vitrification. The CFD models employed for these melter vessels incorporate multiphase fluid flow and heat transfer simulations in different regions, including the plenum, cold cap, and molten glass regions. As the tank waste and glass formers are introduced into the melter, a reacting batch layer known as the cold cap forms on top of the molten glass. To enhance the melt rate, forced convection bubblers located at the bottom of the melters generate convection currents that help homogenize the molten glass and provide heat to the cold cap. As the bubbles rise through the highly viscous glass, they adopt a spherical-cap shape [1]. Meanwhile, the conversion of the batch to glass generates significant amounts of gases (such as water vapor, carbon dioxide, sulfur dioxide, and NOx) due to thermal decomposition [2]. These gases become trapped between the cold cap and molten glass, forming a foam layer [3]. For modeling multiphase flow in CFD and heat transfer simulations of waste glass melters at different scales, efforts are underway to augment the capabilities of the Nek5000 [4] and NekRS [5] open-source codes [6]. Nek5000/NekRS is a scalable and efficient spectral element code that has been successfully applied to a wide range of fluid dynamics problems. By leveraging the Nek5000/NekRS software, it becomes possible to model the melter systems more affordably and with lower computational requirements compared to currently utilized commercial CFD software. The specific objectives of this ongoing effort include: 1. Implementation of a level set method in Nek5000/NekRS: The level set method is a numerical technique commonly employed in CFD simulations to track and represent the interface between different phases or materials accurately. By incorporating this method into Nek5000/NekRS, the ability to simulate multiphase flows in waste glass melters at a high level of fidelity can be achieved. 2. Demonstration of capability for air bubbling through molten glass: As part of the development process, a specific case of air bubbling through molten glass will be simulated using the augmented Nek5000/NekRS code. This demonstration aims to showcase the ability of the software to accurately capture and analyze the complex phenomena involved in the multiphase flow within waste glass melters. By achieving these objectives, the improved Nek5000/NekRS code will offer a powerful computational tool for simulating and analyzing waste glass melter systems, enabling better understanding, optimization, and troubleshooting of these vitrification processes. The ability to accurately model and simulate multiphase flows has broad relevance across many industries and scientific domains, and the improved functionality can contribute to advancements in various fields beyond waste glass melter simulations.
12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
A sharp interface method for compressible liquid–vapor flow with phase transition and surface tension
The numerical approximation of non-isothermal liquid–vapor flow within the compressible regime is a difficult task because complex physical effects at the phase interfaces can govern the global flow behavior. We present a sharp interface approach which treats the interface as a shock-wave like discontinuity. Any mixing of fluid phases is avoided by using the flow solver in the bulk regions only, and a ghost-fluid approach close to the interface. The coupling states for the numerical solution in the bulk regions are determined by the solution of local two-phase Riemann problems across the interface. The Riemann solution accounts for the relevant physics by enforcing appropriate jump conditions at the phase boundary. A wide variety of interface effects can be handled in a thermodynamically consistent way. This includes surface tension or mass/energy transfer by phase transition. Moreover, the local normal speed of the interface, which is needed to calculate the time evolution of the interface, is given by the Riemann solution. The interface tracking itself is based on a level-set method. The focus in this paper is the description of the two-phase Riemann solver and its usage within the sharp interface approach. One-dimensional problems are selected to validate the approach. Finally, the three-dimensional simulation of a wobbling droplet and a shock droplet interaction in two dimensions are shown. In both problems phase transition and surface tension determine the global bulk behavior.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗