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 73 records · Page 4

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↗

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↗

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↗

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↗

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↗

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↗

A broadly-applicable unified closure relation for Taylor bubble rise velocity in pipes with stagnant liquid

Taylor bubble velocity in slug flow is a closure relation which significantly affects the prediction of liquid holdup (or void fraction) and pressure gradient in mechanistic models of slug flow for oil and gas pipe applications. In this work, we use a validated Computational Fluid Dynamics (CFD) approach to simulate the motion of Taylor bubbles in pipes; the interface is tracked with a Level-Set method implemented in a commercial code. A large numerical database is generated covering the most ample range of fluid properties and pipe inclination angles explored to date (Eo ∊ [10, 700],Mo ∊ [1 x 10 -6 , 5 x 10 3 , and θ ∊ [0°, 90°]). A unified Taylor bubble rise velocity correlation is extracted from that database. The new correlation predicts the numerical database with 8.6% absolute average relative error and a coefficient of determination R 2 = 0.97 and other available experimental data with 13.0% absolute average relative error and R 2 = 0.84 outperforming existing correlations and models.

Mechanics↗

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↗

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↗

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↗

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↗

TH-AB-BRB-05: Using a Research Real-Time Control Interface to Go Beyond Dynamic MLC Tracking

Current state-of-the art digital C-arm medical linear accelerators are capable of delivering radiation treatments with high level of automation, which affords coordinated motions of gantry, couch, and multileaf collimator (MLC) with dose rate modulations. The new machine capacity has shown the potential to bring substantially improved radiation dosimetry and/or delivery efficiency to many challenging diseases. Combining an integrated beam orientation optimization algorithm with automated machine navigation, markedly improved dose conformity has been achieved using 4ρ therapy. Trajectory modulated radiation therapy (TMAT) can be used to deliver highly conformal dose to partial breast or to carve complex dose distribution for therapy involving extended volumes such as total marrow and total lymph node treatment. Dynamic electron arc radiotherapy (DEAR) not only overcomes the deficiencies of conventional electron therapy in dose conformity and homogeneity but also achieves so without patient-specific shields. The combination of MLC and couch tracking provides improved motion management of thoracic and abdominal tumors. A substantial body of work has been done in these technological advances for clinical translation. The proposed symposium will provide a timely review of these exciting opportunities. Learning Objectives: Recognize the potential of using digitally controlled linacs for clinically significant improvements in delivered dose distributions for various treatment sites. Identify existing approaches to treatment planning, optimization and delivery for treatment techniques utilizing the advanced functions of digital linacs and venues for further development and improvement. Understand methods for testing and validating delivery system performance. Identify tools available on current delivery systems for implementation and control for such treatments. Obtain the update in clinical applications, trials and regulatory approval. K. Sheng, NIH U19AI067769, NIH R43CA183390, NIH R01CA188300, Varian Medical Systems V. Yu, Varian Medical Systems, AAPM Summer Undergraduate Fellowship, NSF graduate fellowship S. Nill, Elekta AB. Cancer Research UK under Programme C33589/A19727, NIHR Biomedical Research Centre at The Royal Marsden and The Institute of Cancer Research.

60 APPLIED LIFE SCIENCES↗

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↗

Multiscale modelling of evolving foams

We present a set of multi-scale interlinked algorithms to model the dynamics of evolving foams. These algorithms couple the key effects of macroscopic bubble rearrangement, thin film drainage, and membrane rupture. For each of the mechanisms, we construct consistent and accurate algorithms, and couple them together to work across the wide range of space and time scales that occur in foam dynamics. These algorithms include second order finite difference projection methods for computing incompressible fluid flow on the macroscale, second order finite element methods to solve thin film drainage equations in the lamellae and Plateau borders, multiphase Voronoi Implicit Interface Methods to track interconnected membrane boundaries and capture topological changes, and Lagrangian particle methods for conservative liquid redistribution during rearrangement and rupture. We derive a full set of numerical approximations that are coupled via interface jump conditions and flux boundary conditions, and show convergence for the individual mechanisms. We demonstrate our approach by computing a variety of foam dynamics, including coupled evolution of three-dimensional bubble clusters attached to an anchored membrane and collapse of a foam cluster.

97 MATHEMATICS AND COMPUTING↗

Comparison of Mutual Awareness in Analog vs. Digital Control Rooms

Control rooms in nuclear power plants are complex, collaborative working environments rife with potential for human error. As control rooms evolve from analog to digital interface environments, crew communication strategies must change as well. Operators in the analog control room operate the plant and communicate with each other with the help of control boards with large control panels, meters, and switches. Reactor operators could choose a particular set of controls on the particular control board which corresponds to the part of the system of interest. Visual display units (VDU), which became a fixture in the control rooms with the shift to digital human-machine interfaces (HMI), bring benefits for safety and flexibility combining information from various sources. With the increase in automation and the use of digital HMIs operators received their personal sitting workstations where they are able to navigate and operate the plant. The amount of information increased as well as the complexity of monitoring and analyzing the system states. Also computerized systems cause operators to get absorbed by their personal screens and reduce interaction between the crew members. Therefore, the collaboration and coordination of task demands requires vigilance and mutual awareness. The level of automation has an influence on the level of mutual awareness experienced by operators in the control rooms. Teams function to carry out highly interdependent tasks, requiring synthesis and sharing of information to achieve common goals. Improving the performance of the team necessitates the operators to focus on the system, the surrounding environment, and situational awareness about the system states. Collaboration of team members is critical when employing the digital operating systems as the acute awareness of the overall progress of the group and individual activities of the operators drives the team to the desired end-states of the plant operations. Out-of-the loop (OOTL) performance problems can be referred to as the lack of supervisory control, system monitoring, and information processing in the digital HMIs. Under some circumstances Senior Reactor Operators (SRO) struggle to follow the actions of other Reactor Operators (RO) because with the new digital interfaces SROs lose track of how the procedures are carried out by ROs, which results in the loss of three way communication. Mutual awareness plays a significant role in keeping the SROs in the loop of plant operations. This paper will investigate, through literature review and expert interviews, the impact and human errors attributed to the increasing complexity of new digital main control rooms. Comparison of the new digital interfaces with the analog interfaces in nuclear power plants has the potential to shed light on increasing mutual awareness and three way communications between the SROs and ROs. Further, identification of the human errors and human actions can provide beneficial information, suggesting significant and critical improvements that can be implemented in the main control rooms to increase safety and reliability of new technology. Additionally, insights on mutual awareness for analog vs. digital control rooms can serve to improve human reliability.

99 GENERAL AND MISCELLANEOUS↗

The Yale–Potsdam Stellar Isochrones

We introduce the Yale–Potsdam Stellar Isochrones (YaPSI), a new grid of stellar evolution tracks and isochrones of solar-scaled composition. In an effort to improve the Yonsei–Yale database, special emphasis is placed on the construction of accurate low-mass models (M{sub ∗}<0.6 M{sub ⊙}), and in particular on their mass–luminosity and mass–radius relations, both crucial for characterizing exoplanet-host stars, and, in turn, their planetary systems. The YaPSI models cover the mass range 0.15–5.0 M{sub ⊙} densely enough to permit detailed interpolation in mass, and the metallicity and helium abundance ranges [Fe/H] = −1.5 to +0.3 and Y {sub 0} = 0.25–0.37 are specified independently of each other (i.e., no fixed ΔY/ΔZ relation is assumed). The evolutionary tracks are calculated from the pre-main sequence up to the tip of the red giant branch. The isochrones, with ages between 1 Myr and 20 Gyr, provide UBVRI colors in the Johnson–Cousins system, and JHK colors in the homogenized Bessell and Brett system, derived from two different semi-empirical T{sub eff}–color calibrations from the literature. We also provide utility codes, such as an isochrone interpolator, in age, metallicity, and helium content, and an interface of the tracks with an open-source Monte Carlo Markov-Chain tool for the analysis of individual stars. Finally, we present comparisons of the YaPSI models with the best empirical mass–luminosity and mass–radius relations available to date, as well as isochrone fitting of well-studied stellar clusters.

79 ASTRONOMY AND ASTROPHYSICS↗