Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computational fluid dynamics simulation”

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 1,531 records · Page 85
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Airfoil Table Generation and Validation for the VR-12 and SSC-A09 Airfoils and Quadrotor Performance Prediction

As part of NASA’s Urban Air Mobility (UAM) mission, the Revolutionary Vertical Lift Technology (RVLT) project is creating rotorcraft reference designs that can be used by the rotorcraft community to develop quiet, efficient, and safer air vehicles. At NASA, this effort is made possible by the use of the RVLT toolchain software. Airfoil tables are often a critical part of the conceptual design process as they inform both comprehensive analysis and Computational Fluid Dynamics (CFD) codes of 2D aerodynamic coefficients, such as lift, drag, and moment coefficients. These airfoil tables, typically generated during experimental testing, are often proprietary and thus not widely distributable; further, there is not always clear or available reference documentation that provides information on the test conditions. The airfoil tables from these experiments are referred to in this work as Legacy airfoil tables. The scope of this study is to generate airfoil tables with the NASA OVERset grid CFD FLOW solver (OVERFLOW) coupled with the AirFoil Table Generator (AFTGen) software and validate these airfoil tables based on performance predictions calculated in Comprehensive Hierarchical Aeromechanics Rotorcraft Model (CHARM). The overall goal is to provide CFD generated airfoil tables, validate them through comprehensive analysis, and publish the airfoil tables for public distribution. The RVLT quadrotor reference model was selected as the basis for all analyses. The quadrotor uses the VR-12 and SSC-A09 airfoils for its rotor blades. Simulations were run in OVERFLOW for the VR-12 and SSC-A09 airfoils for ±20° angles of attack, Mach numbers from 0.3 to 1.0, and a Reynolds-Mach proportionality constant of 1.023x107. The C81 format airfoil tables from AFTGen were then blended with Legacy NACA 0012 airfoil table data for ±20° to ±180° angle of attack range. CHARM was used to compute thrust coefficient, power coefficient, and figure of merit for an isolated rotor in hover using proprietary legacy airfoil tables for the VR-12 and SSC-A09, as well as with the OVERFLOW generated airfoil tables for the VR-12 and SSC-A09. Overall, the OVERFLOW and Legacy performance predictions agreed well, with the best correlation between simulated and Legacy data observed with the VR-12 airfoil tables.

Airfoil↗
NASA NTRS1987

Development and application of unified algorithms for problems in computational science

A framework is presented for developing computationally unified numerical algorithms for solving nonlinear equations that arise in modeling various problems in mathematical physics. The concept of computational unification is an attempt to encompass efficient solution procedures for computing various nonlinear phenomena that may occur in a given problem. For example, in Computational Fluid Dynamics (CFD), a unified algorithm will be one that allows for solutions to subsonic (elliptic), transonic (mixed elliptic-hyperbolic), and supersonic (hyperbolic) flows for both steady and unsteady problems. The objectives are: development of superior unified algorithms emphasizing accuracy and efficiency aspects; development of codes based on selected algorithms leading to validation; application of mature codes to realistic problems; and extension/application of CFD-based algorithms to problems in other areas of mathematical physics. The ultimate objective is to achieve integration of multidisciplinary technologies to enhance synergism in the design process through computational simulation. Specific unified algorithms for a hierarchy of gas dynamics equations and their applications to two other areas: electromagnetic scattering, and laser-materials interaction accounting for melting.

Shankar, Vijaya↗
NASA NTRS2005

Convergence Acceleration and Documentation of CFD Codes for Turbomachinery Applications

The development and analysis of turbomachinery components for industrial and aerospace applications has been greatly enhanced in recent years through the advent of computational fluid dynamics (CFD) codes and techniques. Although the use of this technology has greatly reduced the time required to perform analysis and design, there still remains much room for improvement in the process. In particular, there is a steep learning curve associated with most turbomachinery CFD codes, and the computation times need to be reduced in order to facilitate their integration into standard work processes. Two turbomachinery codes have recently been developed by Dr. Daniel Dorney (MSFC) and Dr. Douglas Sondak (Boston University). These codes are entitled Aardvark (for 2-D and quasi 3-D simulations) and Phantom (for 3-D simulations). The codes utilize the General Equation Set (GES), structured grid methodology, and overset O- and H-grids. The codes have been used with success by Drs. Dorney and Sondak, as well as others within the turbomachinery community, to analyze engine components and other geometries. One of the primary objectives of this study was to establish a set of parametric input values which will enhance convergence rates for steady state simulations, as well as reduce the runtime required for unsteady cases. The goal is to reduce the turnaround time for CFD simulations, thus permitting more design parametrics to be run within a given time period. In addition, other code enhancements to reduce runtimes were investigated and implemented. The other primary goal of the study was to develop enhanced users manuals for Aardvark and Phantom. These manuals are intended to answer most questions for new users, as well as provide valuable detailed information for the experienced user. The existence of detailed user s manuals will enable new users to become proficient with the codes, as well as reducing the dependency of new users on the code authors. In order to achieve the objectives listed, the following tasks were accomplished: 1) Parametric Study Of Preconditioning Parameters And Other Code Inputs; 2) Code Modifications To Reduce Runtimes; 3) Investigation Of Compiler Options To Reduce Code Runtime; and 4) Development/Enhancement of Users Manuals for Aardvark and Phantom

Marquart, Jed E.↗
NASA NTRSDate not supplied

NASA Langley FUN3D Analyses in Support of the 1st AIAA Stability and Control Prediction Workshop

This work summarizes the results of FUN3D analyses conducted for the 1st AIAA Stability and Control Workshop on behalf of participants from the NASA Langley Research Center. The workshop was created to establish best practices for the prediction of stability and control derivatives using computational fluid dynamics and assess the limitations of these methods when those best practices are applied. The inaugural workshop considered the ONERA version of the NASA/Boeing Common Research Model, which includes the wing, body, horizontal tail, and a vertical tail designed by ONERA. Wind tunnel data at small sideslip angles remain unpublished and served as ‘blind’ data for computational comparisons. The present research generated workshop test case data using the NASA FUN3D code, which is a parallelized, unstructured, node-based, finite-volume discretization, Reynolds-averaged Navier-Stokes flow solver. Steady- state numerical simulations were conducted for workshop test cases investigating the following: grid convergence, Mach number effect on static stability, wind tunnel sting increments, static stability-derivative calculations, and a sideslip angle sweep. Results were generated for two series of unstructured, mixed-element grids, one set provided by the workshop and another set created using the HeldenMesh grid generation software. The results provided include total- and component-level breakdowns of the force and moment coefficients, in addition to sectional pressure distributions for the wing and tail components for comparisons to wind tunnel data.

CFD↗
NASA NTRS2000

Mixing, Noise and Thrust Benefits Using Corrugated Designs

These projects are directed toward the analysis of several concepts for nozzle and inlet performance improvement and noise reduction from jet exhausts. Currently. The FM&AL also initiates new joint research between the HU/FM&AL, the Hyper-X Program Team at the LaRC, and the Central Institute of Aviation Motors (CIAM), Moscow, Russia in the field of optimization of fuel injection and mixing in air-breathing propulsion systems. The main results of theoretical, numerical simulation and experimental tests obtained in the previous research are in the papers and patents. The goals of the 14U/FM&AL programs are twofold: 1) to improve the working efficiency of the HU/FM&AL team in generating new innovative ideas and in conducting research in the field of fluid dynamics and acoustics, basically for improvement of supersonic and subsonic aircraft engines, and 2) to attract promising minority students to this research and training and, in cooperation with other HU departments, to teach them basic knowledge in Aerodynamics, Gas Dynamics, and Theoretical and Experimental Methods in Aeroacoustics and Computational Fluid Dynamics (CFD). The research at the HU/FM&AL supports reduction schemes associated with the emission of en 'ne pollutants for commercial aircraft and concepts for reduction of 91 observables for military aircraft. These research endeavors relate to the goals of the NASA Strategic Enterprise in Aeronautics concerning the development of environmentally acceptable aircraft. It is in this precise area, where the US aircraft industry, academia, and Government are in great need of trained professionals and which is a high priority goal of the Minority University Research and Education (MLTREP) Program, that the HU/FM&AL can make its most important contribution.

Morgan, Morris H., III↗
NASA NTRS2005

Numerical Modeling of Flow Distribution in Micro-Fluidics Systems

This paper describes an application of a general purpose computer program, GFSSP (Generalized Fluid System Simulation Program) for calculating flow distribution in a network of micro-channels. GFSSP employs a finite volume formulation of mass and momentum conservation equations in a network consisting of nodes and branches. Mass conservation equation is solved for pressures at the nodes while the momentum conservation equation is solved at the branches to calculate flowrate. The system of equations describing the fluid network is solved by a numerical method that is a combination of the Newton-Raphson and successive substitution methods. The numerical results have been compared with test data and detailed CFD (computational Fluid Dynamics) calculations. The agreement between test data and predictions is satisfactory. The discrepancies between the predictions and test data can be attributed to the frictional correlation which does not include the effect of surface tension or electro-kinetic effect.

Majumdar, Alok↗
NASA NTRS2009

Automated CFD for Generation of Airfoil Performance Tables

A method of automated computational fluid dynamics (CFD) has been invented for the generation of performance tables for an object subject to fluid flow. The method is applicable to the generation of tables that summarize the effects of two-dimensional flows about airfoils and that are in a format known in the art as C81. (A C81 airfoil performance table is a text file that lists coefficients of lift, drag, and pitching moment of an airfoil as functions of angle of attack for a range of Mach numbers.) The method makes it possible to efficiently generate and tabulate data from simulations of flows for parameter values spanning all operational ranges of actual or potential interest. In so doing, the method also enables filling of gaps and resolution of inconsistencies in C81 tables generated previously from incomplete experimental data or from theoretical calculations that involved questionable assumptions.

Strawn, Roger↗
NASA NTRSDate not supplied

Improved Models for Simulation of Cryogenic Spray Cooling for Propellant Tank Chilldown

This paper presents advances in spray boiling models for cryogenic fluids that exhibit strong variation in thermodynamic properties; it accounts for the interaction of the spray droplets with walls and the complex heat transfer and boiling processes that play a vital role in chilling solid surfaces. The model has been developed by corroborating and validating it with test data in liquid nitrogen obtained at the University of Connecticut. The test data includes detailed visualization and plate temperature measurements for liquid nitrogen sprays impinging on a warm test plate and subsequent transient chilldown of the plate. Simulations of transient plate chilldown in liquid nitrogen were carried out and temperature histories at various locations on the plate were compared with test measurements. The final paper will present sensitivity studies and additional comparisons with test data.

Computational Fluid Dynamics, Cryogenic, Spray Chi↗
NASA NTRSDate not supplied

Prediction of Broadband Blade-Wake Back-Scatter Noise from a Hovering Ideally Twisted Rotor using OVERFLOW2-ANOPP2

This work illustrates the use of a hybrid RANS/LES framework for the prediction of broadband noise via integral formulations from an ideally twisted hovering rotor. Four simulations were conducted, which investigated the effect of different turbulence modeling and spatiotemporal resolutions. The aerodynamic performance and aeroacoustics predicted by these four computational cases were compared against data acquired in the Small Hover Anechoic Chamber at the NASA Langley Research Center and against previous simulations conducted using the lattice-Boltzmann solver, PowerFLOW. The simulation case with the finest spatial resolution had the most favorable aerodynamic performance and broadband noise comparison with measured data; however, the simulation case utilizing the SST-DDES turbulence model predicted the most accurate tonal noise signature. Aerodynamic flowfield visualizations identified near-wake trailing edge vortex shedding, perpendicular blade-vortex interactions, and tip vortex formation, which are known to cause broadband noise. The fluid dynamics associated with a new broadband noise source, blade-wake back-scatter, were also identified. Lastly, a broadband noise source visualization technique was utilized to confirm the presence of blade-wake back-scatter as a dominant broadband noise source.

blade-wake back-scatter↗
NASA NTRS2001

Fluid Aspects of Solar Wind Disturbances Driven by Coronal Mass Ejections

Transient disturbances in the solar wind initiated by coronal eruptions have been modeled for many years, beginning with the self-similar analytical models of Parker and Simon and Axford. The first numerical computer code (one-dimensional, gas dynamic) to study disturbance propagation in the solar wind was developed in the late 1960s, and a variety of other codes ranging from simple one-dimensional gas dynamic codes through three-dimensional gas dynamic and magnetohydrodynamic codes have been developed in subsequent years. For the most part, these codes have been applied to the problem of disturbances driven by fast CMEs propagating into a structureless solar wind. Pizzo provided an excellent summary of the level of understanding achieved from such simulation studies through about 1984, and other reviews have subsequently become available. More recently, some attention has been focused on disturbances generated by slow CMEs, on disturbances driven by CMEs having high internal pressures, and disturbance propagation effects associated with a structured ambient solar wind. Our purpose here is to provide a brief tutorial on fluid aspects of solar wind disturbances derived from numerical gas dynamic simulations. For the most part we illustrate disturbance evolution by propagating idealized perturbations, mimicking different types of CMEs, into a structureless solar wind using a simple one-dimensional, adiabatic (except at shocks), gas dynamic code. The simulations begin outside the critical point where the solar wind becomes supersonic and thus do not address questions of how the CMEs themselves are initiated. Limited to one dimension (the radial direction), the simulation code predicts too strong an interaction between newly ejected solar material and the ambient wind because it neglects azimuthal and meridional motions of the plasma that help relieve pressure stresses. Moreover, the code ignores magnetic forces and thus also underestimates the speed with which pressure disturbances propagate in the wind.

Gosling, J. T.↗
NASA NTRS1998

Numerical Experiment with Time and Spatial Accuracy of Navier-Stokes Computation For Helicopter Problems

Helicopter flowfields are highly unsteady, nonlinear and three-dimensional. In forward flight and in hover, the rotor blades interact with the tip vortex and wake sheet developed by either itself or the other blades. This interaction, known as blade-vortex interactions (BVI), results in unsteady loading of the blades and can cause a distinctive acoustic signature. Accurate and cost-effective computational fluid dynamic solutions that capture blade-vortex interactions can help rotor designers and engineers to predict rotor performance and to develop designs for low acoustic signature. Such a predictive method must preserve a blade's shed vortex for several blade revolutions before being dissipated. A number of researchers have explored the requirements for this task. This paper will outline some new capabilities that have been added to the NASA Ames' OVERFLOW code to improve its overall accuracy for both vortex capturing and unsteady flows. To highlight these improvements, a number of case studies will be presented. These case studies consist of free convection of a 2-dimensional vortex, dynamically pitching 2-D airfoil including light-stall, and a full 3-D unsteady viscous solution of a helicopter rotor in forward flight In this study both central and upwind difference schemes are modified to be more accurate. Central difference scheme is chosen for this simulation because the flowfield is not dominated by strong shocks. The feature of shock-vortex interaction in such a flow is less important than the dominant blade-vortex interaction. The scheme is second-order accurate in time and solves the thin-layer Navier-Stokes equations in fully-implicit manner at each time-step. The spatial accuracy is either second and fourth-order central difference or third-order upwind difference using Roe-flux and MUSCLE scheme. This paper will highlight and demonstrate the methods for several sample cases and for a helicopter rotor. Preliminary computations on a rotor were performed by using this method and are in the process of documentation.

Ahmad, Jasim↗
NASA NTRS2004

Exhaust Simulation Testing of a Hypersonic Airbreathing Model at Transonic Speeds

An experimental study was performed to examine jet-effects for an airframe-integrated, scramjet-rocket combined-cycle vehicle configuration at transonic test conditions. This investigation was performed by testing an existing exhaust simulation wind tunnel model, known as Model 5B, in the NASA Langley 16-Ft. Transonic Tunnel. Tests were conducted at freestream Mach numbers from 0.7 to 1.2, at angles of attack from 2 to +14 degrees, and at up to seven nozzle static pressure ratio values for a set of horizontal-tail and body-flap deflections. The model aftbody, horizontal tails, and body flaps were extensively pressure instrumented to provide an understanding of jet-effects and control-surface/plume interactions, as well as for the development of analytical methodologies and calibration of computational fluid dynamic codes to predict this type of flow phenomenon. At all transonic test conditions examined, the exhaust flow at the exit of the internal nozzle was over-expanded, generating an exhaust plume that turned toward the aftbody. Pressure contour plots for the aftbody of Model 5B are presented for freestream transonic Mach numbers of 0.70, 0.95, and 1.20. These pressure data, along with shadowgraph images, indicated the impingement of an internal plume shock and at least one reflected shock onto the aftbody for all transonic conditions tested. These results also provided evidence of the highly three-dimensional nature of the aftbody exhaust flowfield. Parametric testing showed that angle-of-attack, static nozzle pressure ratio, and freestream Mach number all affected the exhaust-plume size, exhaust-flowfield shock structure, and the aftbody-pressure distribution, with Mach number having the largest effect. Integration of the aftbody pressure data showed large variations in the pitching moment throughout the transonic regime.

Huebner, Lawrence D.↗
NASA NTRS2019

Computational Fluid Dynamics Methods Used in the Development of the Space Launch System Liftoff and Transition Lineloads Databases

The objective of this paper is to document the reasoning and trade studies that supported the selection of appropriate tools for constructing aerodynamic lineload databases for the Liftoff and Transition phases of flight for launch vehicles. These decisions were made amid the maturation of an evolving workflow for generating databases on variants of the Space Launch System launch vehicle, with most being based on results from brief developmental studies performed in response to specific, unforeseen challenges that were encountered in analyzing a given configuration. This report is intended to provide a summary of the results and the decision-making processes chronologically over the design cycles of various configurations, starting with isolated free-air bodies for the Block 1 Crew, then the Block 1B Crew and Cargo configurations, and most recently the Block 1B Crew configuration in proximity to the launch tower. The results from these analyses led to the selection of the CREATE-AV Kestrel flowsolver for simulating these problems. The need to accurately capture the expected leeward-wake flow field characteristics required the use of Delayed Detached Eddy Simulation (DDES) method, for which the vorticity magnitude was employed as the solution Adaptive Mesh Refinement (AMR) function over the off-body Cartesian grid region. In addition, the Spalart-Allmaras (SA) model is used to account for the flow turbulence effects.

Ratnayake, Nalin A.↗
NASA NTRS1989

Rarefied-continuum gas dynamics transition for SUMS project

This program is to develop an analytic method for reducing SUMS data for the determination of the undisturbed atmosphere conditions ahead of the shuttle along its descending trajectory. It is divided into an internal flow problem, an external flow problem and their matching conditions. Since the existing method of Direct Simulation Monte Carlo (DSMC) failed completely for the internal flow problem, the emphasis is on the internal flow of a highly non-equilibrium, rarefied air through a short tube of a diameter much less than the gaseous mean free path. A two fluid model analysis of this internal flow problem has been developed and studied with typical results illustrated. A computer program for such an analysis and a technical paper published in Lecture Notes in Physics No. 323 (1989) are included as Appendices 3 and 4. A proposal for in situ determination of the surface accommodation coefficients sigma sub t and sigma e is included in Appendix 5 because of their importance in quantitative data reduction. A two fluid formulation for the external flow problem is included as Appendix 6 and a review article for AIAA on Hypersonic propulsion, much dependent on ambient atmospheric density, is also included as Appendix 7.

Cheng, Sin-I↗
NASA NTRS2003

Using CFD as Rocket Injector Design Tool: Recent Progress at Marshall Space Flight Center

The choice of tools used for injector design is in a transitional phase between exclusive reliance on the empirically based correlations and extensive use of computational fluid dynamics (CFD). The Next Generation Launch Technology (NGLT) Program goals emphasizing lower costs and increased reliability have produced a need to enable CFD as an injector design tool in a shorter time frame. This is the primary objective of the Staged Combustor Injector Technology Task currently under way at Marshall Space Flight Center (MSFC). The documentation of this effort begins with a very brief status of current injector design tools. MSFC's vision for use of CFD as a tool for combustion devices design is stated and discussed with emphasis on the injector. The concept of the Simulation Readiness Level (SRL), comprised of solution fidelity, robustness and accuracy, is introduced and discussed. This quantitative measurement is used to establish the gap between the current state of demonstrated capability and that necessary for regular use in the design process. MSFC's view of the validation process is presented and issues associated with obtaining the necessary data are noted and discussed. Three current experimental efforts aimed at generating validation data are presented. The importance of uncertainty analysis to understand the data quality is also demonstrated. First, a brief status of current injector design tools is provided as context for the current effort. Next, the MSFC vision for using CFD as an injector design tool is stated. A generic CFD-based injector design methodology is also outlined and briefly discussed. Three areas where MSFC is using injector CFD analyses for program support will be discussed. These include the Integrated Powerhead Development (IPD) engine which uses hydrogen and oxygen propellants in a full flow staged combustion (FFSC) cycle and the TR-107 and the RS84 engine both of which use RP-1 and oxygen in an ORSC cycle. Finally, an attempt is made to objectively summarize what progress has been made at MSFC in enabling CFD as an injector design tool.

Tucker, Kevin↗
NASA NTRS2011

CFD Study of the J2-X LOX Two-Bladed Inducer in the Presence of a Circumferential Groove

Results are presented of a computational fluid dynamics (CFD) study done in support of water flow experiments of the J-2X Oxidizer Turbopump (OTP) 2-bladed inducer with a circumferential groove that were conducted at Marshall Space Flight Center (MSFC). Sub-scale water flow testing results indicate that the circumferential groove greatly reduces synchronous cavitation and subsequent bearing loads at a minimal performance cost, but the energy reappears as high order cavitation (HOC) that spans a relatively large frequency range. Thus, HOC may have implications for the full-scale OTP inducer in terms of reduced structural margin at higher mode frequencies. Simulations using the LOCI-Stream CFD program were conducted in order to explore the root physical cause of the HOC. It was found that the axial recirculation pattern in the circumferential groove causes high-swirl fluid to interact with the nearly-axial incoming fluid just above the inducer blades. The high-shear interface between the fluids is Kelvin-Helmholtz unstable, resulting in trains of low pressure regions or 'pearls' forming near the upstream edge of the groove. When the pressure in these regions becomes low enough and they get cut by the blade leading edge, HOC is thought to occur. Although further work is required, the numerical models indicate that the root cause of HOC is hydrodynamic. That is, the pearls are always present, even when cavitation is not. Comparisons to ongoing water flow experiments will be discussed, as well as predictions for the full-scale OTP inducer.

Coker, Robert↗
NASA NTRS1995

The 3D Navier-Stokes analysis of a Mach 2.68 bifurcated rectangular mixed-compression inlet

The supersonic diffuser of a Mach 2.68 bifurcated, rectangular, mixed-compression inlet was analyzed using a three-dimensional (3D) Navier-Stokes flow solver. A two-equation turbulence model, and a porous bleed model based on unchoked bleed hole discharge coefficients were used. Comparisons were made with experimental data, inviscid theory, and two-dimensional Navier-Stokes analyses. The main objective was to gain insight into the inlet fluid dynamics. Examination of the computational results along with the experimental data suggest that the cowl shock-sidewall boundary layer interaction near the leading edge caused a substantial separation in the wind tunnel inlet model. As a result, the inlet performance may have been compromised by increased spillage and higher bleed mass flow requirements. The internal flow contained substantial waves that were not in the original inviscid design. 3D effects were fairly minor for this inlet at on-design conditions. Navier-Stokes analysis appears to be an useful tool for gaining insight into the inlet fluid dynamics. It provides a higher fidelity simulation of the flowfield than the original inviscid design, by taking into account boundary layers, porous bleed, and their interactions with shock waves.

Mizukami, M.↗
NASA NTRS2000

Computations of Boiling in Microgravity

The absence (or reduction) of gravity, can lead to major changes in boiling heat transfer. On Earth, convection has a major effect on the heat distribution ahead of an evaporation front, and buoyancy determines the motion of the growing bubbles. In microgravity, convection and buoyancy are absent or greatly reduced and the dynamics of the growing vapor bubbles can change in a fundamental way. In particular, the lack of redistribution of heat can lead to a large superheat and explosive growth of bubbles once they form. While considerable efforts have been devoted to examining boiling experimentally, including the effect of microgravity, theoretical and computational work have been limited. Here, the growth of boiling bubbles is studied by direct numerical simulations where the flow field is fully resolved and the effects of inertia, viscosity, surface deformation, heat conduction and convection, as well as the phase change, are fully accounted for. Boiling involves both fluid flow and heat transfer and thus requires the solution of the Navier-Stokes and the energy equations. The numerical method is based on writing one set of governing transport equations which is valid in both the liquid and vapor phases. This local, single-field formulation incorporates the effect of the interface in the governing equations as source terms acting only at the interface. These sources account for surface tension and latent heat in the equations for conservation of momentum and energy as well as mass transfer across the interface due to phase change. The single-field formulation naturally incorporates the correct mass, momentum and energy balances across the interface. Integration of the conservation equations across the interface directly yields the jump conditions derived in the local instant formulation for two-phase systems. In the numerical implementation, the conservation equations for the whole computational domain (both vapor and liquid) are solved using a stationary grid and the phase boundary is followed by a moving unstructured two-dimensional grid. While two-dimensional simulations have been used for preliminary studies and to examine the resolution requirement, the focus is on fully three-dimensional simulations. The numerical methodology, including the parallelization and grid refinement strategy is discussed, and preliminary results shown. For buoyancy driven flow, the heat transfer is in good agreement with experimental correlations. The changes when gravity is turned off and/or fluid shear is added are discussed, as well as the difference between simulations of a layer freely releasing bubbles versus simulations using only one wavelength initial perturbation. Figure 1 shows the early stages of the formation of a three-dimensional bubble from a thin vapor layer. The boundary conditions are periodic in the x and y direction, the bottom is a hot and the top allows a free outflow. The jagged edge of the surface close to the bottom of the computational domain is due to some of the surface elements being on the other side of the domain and some elements not plotted by our plotting routine. In the second figure, we show the temperature distribution through two perpendicular planes.

Tryggvason, G.↗