SEARCH · Search NASA

Results for “Two Phase Flow”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,063 records · Page 59

All sources NASA NTRS DOE OSTI

DOE OSTI2014

Thickness-based adaptive mesh refinement methods for multi-phase flow simulations with thin regions

In numerical simulations of multi-scale, multi-phase flows, grid refinement is required to resolve regions with small scales. A notable example is liquid-jet atomization and subsequent droplet dynamics. It is essential to characterize the detailed flow physics with variable length scales with high fidelity, in order to elucidate the underlying mechanisms. In this paper, two thickness-based mesh refinement schemes are developed based on distance- and topology-oriented criteria for thin regions with confining wall/plane of symmetry and in any situation, respectively. Both techniques are implemented in a general framework with a volume-of-fluid formulation and an adaptive-mesh-refinement capability. The distance-oriented technique compares against a critical value, the ratio of an interfacial cell size to the distance between the mass center of the cell and a reference plane. The topology-oriented technique is developed from digital topology theories to handle more general conditions. The requirement for interfacial mesh refinement can be detected swiftly, without the need of thickness information, equation solving, variable averaging or mesh repairing. The mesh refinement level increases smoothly on demand in thin regions. The schemes have been verified and validated against several benchmark cases to demonstrate their effectiveness and robustness. These include the dynamics of colliding droplets, droplet motions in a microchannel, and atomization of liquid impinging jets. Overall, the thickness-based refinement technique provides highly adaptive meshes for problems with thin regions in an efficient and fully automatic manner.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

Atomic resolution scanning transmission electron microscopy at liquid helium temperatures for quantum materials

Fundamental quantum phenomena in condensed matter, ranging from correlated electron systems to quantum information processors, manifest their emergent characteristics and behaviors predominantly at low temperatures. This necessitates the use of liquid helium (LHe) cooling for experimental observation. Atomic resolution scanning transmission electron microscopy combined with LHe cooling (cryo-STEM) provides a powerful characterization technique to probe local atomic structural modulations and their coupling with charge, spin and orbital degrees-of-freedom in quantum materials. However, achieving atomic resolution in cryo-STEM is exceptionally challenging, primarily due to sample drifts arising from temperature changes and noises associated with LHe bubbling, turbulent gas flow, etc. In this work, we demonstrate atomic resolution cryo-STEM imaging at LHe temperatures using a commercial side-entry LHe cooling holder. Firstly, we examine STEM imaging performance as a function of He gas flow rate, identifying two primary noise sources: He-gas pulsing and He-gas bubbling. Secondly, we propose two strategies to achieve low noise conditions for atomic resolution STEM imaging: either by temporarily suppressing He gas flow rate using the needle valve or by acquiring images during the natural warming process. Lastly, we show the applications of image acquisition methods and image processing techniques in investigating structural phase transitions in Cr 2 Ge 2 Te 6 , CuIr 2 S 4 , and CrCl 3 . In conclusion, our findings represent an advance in the field of atomic resolution electron microscopy imaging for quantum materials and devices at LHe temperatures, which can be applied to other commercial side-entry LHe cooling TEM holders.

36 MATERIALS SCIENCE↗

DOE OSTI2017

Optimal Power Flow in Multiphase Radial Networks with Delta Connections: Preprint

This paper focuses on multiphase radial distribution networks with mixed wye and delta connections, and proposes a semidefinite relaxation of the AC optimal power flow (OPF) problem. Two multiphase power-flow models are developed to facilitate the integration of delta-connected generation units/loads in the OPF problem. The first model extends traditional branch flow models - and it is referred to as extended branch flow model (EBFM). The second model leverages a linear relationship between per-phase power injections and delta connections, which holds under a balanced voltage approximation (BVA). Based on these models, pertinent OPF problems are formulated and relaxed to semidefinite programs (SDPs). Numerical studies on IEEE test feeders show that SDP relaxations can be solved efficiently by a generic optimization solver. Numerical evidences indicate that solving the resultant SDP under BVA is faster than under EBFM. Moreover, both SDP solutions are numerically exact with respect to voltages and branch flows. It is also shown that the SDP solution under BVA has a small optimality gap, while the BVA model is accurate in the sense that it reflects actual system voltages.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2016

The NETL MFiX Suite of multiphase flow models: A brief review and recent applications of MFiX-TFM to fossil energy technologies

Here, the MFiX suite of multiphase computational fluid dynamics (CFD) codes is being developed at U.S. Department of Energy's National Energy Technology Laboratory (NETL). It includes several different approaches to multiphase simulation: MFiX-TFM, a two-fluid (Eulerian–Eulerian) model; MFiX-DEM, an Eulerian fluid model with a Lagrangian Discrete Element Model for the solids phase; and MFiX-PIC, Eulerian fluid model with Lagrangian particle ‘parcels’ representing particle groups. These models are undergoing continuous development and application, with verification, validation, and uncertainty quantification (VV&UQ) as integrated activities. After a brief summary of recent progress in the verification, validation and uncertainty quantification (VV&UQ), this article highlights two recent accomplishments in the application of MFiX-TFM to fossil energy technology development. First, recent application of MFiX to the pilot-scale KBR TRIG™ Transport Gasifier located at DOE's National Carbon Capture Center (NCCC) is described. Gasifier performance over a range of operating conditions was modeled and compared to NCCC operational data to validate the ability of the model to predict parametric behavior. Second, comparison of code predictions at a detailed fundamental scale is presented studying solid sorbents for the post-combustion capture of CO 2 from flue gas. Specifically designed NETL experiments are being used to validate hydrodynamics and chemical kinetics for the sorbent-based carbon capture process.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2015

The cementitious barriers partnership experimental programs and software advancing DOE's waste disposal/tank closure efforts - 15436

The U.S. Department of Energy Environmental Management (DOE-EM) Office of Tank Waste Management sponsored Cementitious Barriers Partnership (CBP) is chartered with providing the technical basis for implementing cement-based waste forms and radioactive waste containment structures for long-term disposal. DOE needs in this area include the following to support progress in final treatment and disposal of legacy waste and closure of High-Level Waste (HLW) tanks in the DOE complex: - Long-term performance predictions and - Flow sheet development and flow sheet enhancements - Conceptual designs for new disposal facilities The DOE-EM Cementitious Barriers Partnership is producing software and experimental programs resulting in new methods and data needed for end-users involved with environmental cleanup and waste disposal. Both the modeling tools and the experimental data have already benefited the DOE sites in the areas of performance assessments by increasing confidence backed up with modeling support, leaching methods, and transport properties developed for actual DOE materials. In 2014, the CBP Partnership released the CBP Software Toolbox - 'Version 2.0' which provides concrete degradation models for: 1) sulfate attack, 2) carbonation, 3) chloride initiated rebar corrosion, and includes constituent leaching. These models are applicable and can be used by both DOE and the Nuclear Regulatory Commission (NRC) for service life and long-term performance evaluations and predictions of nuclear and radioactive waste containment structures across the DOE complex (Figure 1) including Future SRS salt-stone and HLW tank performance assessments and special analyses - Hanford site HLW tank closure projects and other projects in which cementitious barriers are required - The Advanced Simulation Capability for Environmental Management (ASCEM) project which requires source terms from cementitious containment structures as input to their flow simulations - Regulatory reviews of DOE performance assessments - Nuclear Regulatory Commission reviews of commercial nuclear power plant (NPP) structures which are part of the overall US Energy Security program to extend the service life of NPPs In addition, the CBP experimental programs have had a significant impact on the DOE complex by providing specific data unique to DOE sodium salt wastes at Hanford and SRS which are not readily available in the literature. Two recent experimental programs on cementitious phase characterization and on technetium (Tc) mobility have provided significant conclusions as summarized below: Recent mineralogy characterization discussed in this paper illustrates that sodium salt waste form matrices are somewhat similar to but not the same as those found in blended cement matrices which to date have been used in long-term thermodynamic modeling and contaminant sequestration as a first approximation. Utilizing the CBP generated data in long-term performance predictions provides for a more defensible technical basis in performance evaluations. - In addition, recent experimental studies related to technetium mobility indicate that conventional leaching protocols may not be conservative for direct disposal of Tc- containing waste forms in vadose zone environments. These results have the potential to influence the current Hanford supplemental waste treatment flow sheet and disposal conceptual design. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2015

Predictions for the reverberating spectral line from a newly formed black hole accretion disk: Case of tidal disruption flares

Tidal disruption events (TDEs) can be perfect probes of dormant supermassive black holes in normal galaxies. During the rising phase, the accretion luminosity can increase by orders of magnitude in several weeks, and the emergent ionizing radiation illuminates the fresh accretion flow. In this paper, we simulated the evolution of the expected spectral line profile of iron due to such a flare by using a ray-tracing code with effects of general relativity taken into account. We found that the time-dependent profile changes significantly with black hole spin, inclination angle with respect to the black hole equatorial plane, and the expansion velocity of the ionization front. At low values of spin, a “loop” feature appears in the line profile versus time plot when the inclination is no less than 30° and the expansion velocity v{sub exp} is no less than half the speed of light, owing to a shadow in the emission of the truncated disk. In the light curve two peaks occur depending on the inclination angle. At large v{sub exp}, a shallow “nose” feature may develop ahead of the loop; its duration depends on the expansion velocity and the inclination angle. We explore the entire interval of black hole spin parameter ranging from extreme prograde to extreme retrograde rotation, −1<a<1. In the prograde case, a low-energy tail appears to be more pronounced in the evolving centroid energy of the line. Our results demonstrate the importance of searching for X-ray spectral lines in the early phase of TDE flares in order to constrain black hole mass and spin, as well as properties of the innermost accretion flow.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2013

Statistics of particle time-temperature histories :

Progress toward predictions of the statistics of particle time-temperature histories is presented. These predictions are to be made using Lagrangian particle models within the one-dimensional turbulence (ODT) model. In the present reporting period we have further characterized the performance, behavior and capabilities of the particle dispersion models that were added to the ODT model in the first period. We have also extended the capabilities in two manners. First we provide alternate implementations of the particle transport process within ODT; within this context the original implementation is referred to as the type-I and the new implementations are referred to as the type-C and type-IC interactions. Second we have developed and implemented models for two-way coupling between the particle and fluid phase. This allows us to predict the reduced rate of turbulent mixing associated with particle dissipation of energy and similar phenomena. Work in characterizing these capabilities has taken place in homogeneous decaying turbulence, in free shear layers, in jets and in channel flow with walls, and selected results are presented.

Hewson, John C.↗

DOE OSTI2021

White paper on NIF Discovery Science and frontier regimes of HEDP

Over the past two decades, experimental capabilities at high energy density (HED) facilities such as the National Ignition Facility (NIF), Omega, and Omega EP lasers, and the Z pulsed power facility are able to make precision measurements of fundamental quantities such as equations of state (EOS), material phase, opacities, and ionization levels Z. Flow based quantities are also being studied, such as hydrodynamic instabilities, evolution into turbulence, turbulent dynamo magnetic field amplification, magnetic reconnection; plasma instabilities, plasma wakefield particle acceleration; and relativistic plasma phenomena, such as (e + ,e - ) pair creation. Experimental regimes can be created that reproduce conditions relevant to planetary and exoplanet interiors, stellar interiors, the interiors of brown dwarfs, and the envelope conditions of white dwarf stars. Diagnostic techniques have been developed hand in hand with these new experimental regimes that can be accessed. In this white paper, a brief summary of key experimental achievements from the NIF Discovery Science (basic science) program will be presented; followed by a discussion of current challenges, needs, and potential paths forward.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Reaction Temperature Manipulation as a Process Intensification Approach for CO 2 Absorption

Reactor temperature manipulation to increase product yields of chemical reactions is a known technique used in many industrial processes. In the case of exothermic chemical reactions, the well-known Le Chatelier’s principle predicts that a decrease in temperature will displace the chemical reaction toward the formation of products by increasing the value of the equilibrium constant. The reverse is true for endothermic reactions. Reactor temperature manipulation in an industrial system, however, affects the values of many variables, including physical properties, transport parameters, reaction kinetic parameters, etc. In the case of reactive absorption, some variables change with increasing temperatures due to solute absorption, while others change in such a way that the solute absorption rate decreases. For example, temperature drop increases product formation for exothermic reactions but reduces the value of transport parameters, leading to decreasing interfacial concentrations and absorption rates. Therefore, temperature manipulation strategies must be designed carefully to achieve the process goals. In this work, we theoretically study the use of temperature as a tool to increase CO 2 absorption by solvents in a semi-batch reactor. A computer code has been developed and validated using reported experimental data. Calculated results demonstrate an increase in absorbed CO 2 of more than 28% with respect to the highest temperature used. Despite high agitation and high gas flow rate, the system is mass transfer controlled at short times, becoming kinetically controlled as time increases. An operating strategy to decrease cooling energy costs is also proposed. This study reveals that reactor temperature manipulation can be an effective process to improve CO 2 absorption by solvents in two-phase semi-batch reactors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2019

Deposition Mitigation in Flowing Systems Using Coatings

Abstract In pipelines, solid compounds including gas hydrates and asphaltenes may form/precipitate and accumulate on the pipe surface, leading to a gradual stenosis of the flowline. As a result, production may become increasingly difficult or possibly interrupted if mitigation efforts are not enacted. Typically, injected chemicals will either inhibit nucleation or dissolve already-formed deposits to restore original flow conditions back to the system; however, this can be a costly option. More recently, management strategies have been proposed where solids are handled in a controlled fashion rather than completely avoided. One such proposed management strategy as suggested for wall deposit formation is the use of coatings. Here, coatings can provide a low surface energy layer on the pipe wall, which restricts liquid and solid accumulation, allowing for a stable slurry flow through a system. This study utilized two material formulations within several experimental setups to probe their interactions with water, gas hydrate, asphaltene, and crude oil. The results serve as part of an ongoing investigation into a surface treatment formulation that can be tested on larger-scale, fully flowing systems, which could be ultimately implemented into real-world production scenarios. The first surface treatment is a water-based polymeric surface that displays repellency to both oil and water phases (omniphobic). Testing of this material consisted of water contact angle measurements and static asphaltene/crude oil deposition quantification at atmospheric conditions, as well as visual confirmation of hydrate deposition prevention at high pressures. Additionally, an experimental superomniphobic surface treatment, which displays elevated resiliency to both water and hydrocarbons, was also examined within the asphaltene/crude oil test as a comparison to the omniphobic surface treatment. Static contact angle results showed that the omniphobic surface treatment had reduced surface interaction with water droplets in air, increasing the low contact angles of corroded surfaces (0-31°) to slightly hydrophobic conditions of 91.5°. Additionally, rocking cells tests indicated that these omniphobic surface treatments may prevent gas hydrate deposition under high-pressure, semi-flowing conditions. Multiple tests found that formed hydrate agglomerants did not deposit for at least 48 and 72 hours. Finally, static deposition tests conducted in crude oil with forced asphaltene precipitation suggested that the omniphobic surface treatment displayed a resistance to both asphaltenes and crude oil when compared to untreated and superomniphobic surfaces.

03 NATURAL GAS↗

DOE OSTI2020

Thin hybrid capillary two-phase cooling system

A novel hybrid two-phase cooling system was developed that integrated a mechanically pumped two-phase loop with a capillary-driven two-phase cooling device. The latter cooling mechanism was based on evaporation/boiling from wick structures made by sintering copper particles on the interior surfaces of a copper cold plate. The cold plate provided cooling to two surfaces and each of them included four heaters in series. The novelty of the developed technology was preventing flooding of the evaporator wicks by isolating the evaporation surface from the pumped liquid flow that fed them. This arrangement allowed for a high liquid feed flow rate much greater than would be allowed by a capillary pumped system while maintaining a low thermal resistance at the evaporation surface. Using this approach, the cooling system removed over 850 W with a low pumping power below 1.0 W while using R245fa as the working fluid. The equivalent heat fluxes exceeded 970 W/cm 2 over areas less than 0.12 cm 2 . The measured thermal resistance was as low as 0.09 K-cm 2 /W. The presented thermal management solution enables an increase in the power of high heat flux electronic devices beyond the state-of-the-art.

42 ENGINEERING↗

DOE OSTI2020

Computational Fluid Dynamics Analysis of a 12 MW Circulating Fluidized Bed Rise

A three-dimensional (3D) two fluid model (TFM) model of the 12 MW circulating fluidized bed (CFB) riser at Chalmers University was developed using the open-source Multiphase Flow with Interphase eXchanges (MFiX) Software Suite. The hydrodynamic behavior of the riser was validated by comparing axial pressure profiles from experimental measurements with the results from the simulations. Under some of the operating conditions where there was no external circulation of solids, the pressure profiles from the simulation matched well with the experimental measurements. For the cases where there was external circulation of solids, the pressure profiles with the monodisperse solid phase modeling underpredicted the solids concentration in the freeboard region. With the inclusion of polydispersity in the model, the predicted pressure profiles agreed well with the experimental measurements. A two-dimensional (2D) model of the riser was used to perform simulations of coal combustion. The predicted flue gas outlet compositions of oxygen and carbon dioxide compared reasonably to experiments. The experiments also reported some concentration of carbon monoxide which was overpredicted in the simulations. Additional experimental data such as axial oxygen concentration profile and pressure profiles were also compared with the experiment. It was observed that the oxygen concentration in the bottom bed from the simulations was higher than in the experiments, which suggested that most of the oxidation reactions occurred in the bottom bed. The concentration of volatiles and char in the simulations was lower in the bottom bed than in the experiments, which suggest that there is not sufficient penetration of the solid phase species into the bed.

20 FOSSIL-FUELED POWER PLANTS↗

DOE OSTI2015

The Cementitious Barriers Partnership Experimental Programs and Software Advancing DOE’s Waste Disposal/Tank Closure Efforts – 15436

The U.S. Department of Energy Environmental Management (DOE-EM) Office of Tank Waste Management-sponsored Cementitious Barriers Partnership (CBP) is chartered with providing the technical basis for implementing cement-based waste forms and radioactive waste containment structures for long-term disposal. DOE needs in this area include the following to support progress in final treatment and disposal of legacy waste and closure of High-Level Waste (HLW) tanks in the DOE complex: long-term performance predictions, flow sheet development and flow sheet enhancements, and conceptual designs for new disposal facilities. The DOE-EM Cementitious Barriers Partnership is producing software and experimental programs resulting in new methods and data needed for end-users involved with environmental cleanup and waste disposal. Both the modeling tools and the experimental data have already benefited the DOE sites in the areas of performance assessments by increasing confidence backed up with modeling support, leaching methods, and transport properties developed for actual DOE materials. In 2014, the CBP Partnership released the CBP Software Toolbox –“Version 2.0” which provides concrete degradation models for 1) sulfate attack, 2) carbonation, and 3) chloride initiated rebar corrosion, and includes constituent leaching. These models are applicable and can be used by both DOE and the Nuclear Regulatory Commission (NRC) for service life and long-term performance evaluations and predictions of nuclear and radioactive waste containment structures across the DOE complex, including future SRS Saltstone and HLW tank performance assessments and special analyses, Hanford site HLW tank closure projects and other projects in which cementitious barriers are required, the Advanced Simulation Capability for Environmental Management (ASCEM) project which requires source terms from cementitious containment structures as input to their flow simulations, regulatory reviews of DOE performance assessments, and Nuclear Regulatory Commission reviews of commercial nuclear power plant (NPP) structures which are part of the overall US Energy Security program to extend the service life of NPPs. In addition, the CBP experimental programs have had a significant impact on the DOE complex by providing specific data unique to DOE sodium salt wastes at Hanford and SRS which are not readily available in the literature. Two recent experimental programs on cementitious phase characterization and on technetium (Tc) mobility have provided significant conclusions as follows: recent mineralogy characterization discussed in this paper illustrates that sodium salt waste form matrices are somewhat similar to but not the same as those found in blended cement matrices which to date have been used in long-term thermodynamic modeling and contaminant sequestration as a first approximation. Utilizing the CBP generated data in long-term performance predictions provides for a more defensible technical basis in performance evaluations. In addition, recent experimental studies related to technetium mobility indicate that conventional leaching protocols may not be conservative for direct disposal of Tc-containing waste forms in vadose zone environments. These results have the potential to influence the current Hanford supplemental waste treatment flow sheet and disposal conceptual design.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2017

Atmospheric Pressure Ionization Using a High Voltage Target Compared to Electrospray Ionization

A new atmospheric pressure ionization (API) source, viz. UniSpray, was evaluated for mass spectrometry (MS) analysis of pharmaceutical compounds by head-to-head comparison with electrospray ionization (ESI) on the same high-resolution MS system. The atmospheric pressure ionization source is composed of a grounded nebulizer spraying onto a high voltage, cylindrical stainless steel target. Molecules are ionized in a similar fashion to electrospray ionization, predominantly producing protonated or deprotonated species. Adduct formation (e.g., proton and sodium adducts) and in-source fragmentation is shown to be almost identical between the two sources. The performance of the new API source was compared with electrospray by infusion of a mix of 22 pharmaceutical compounds with a wide variety of functional groups and physico-chemical properties (molecular weight, logP, and pKa) in more than 100 different conditions (mobile phase strength, solvents, pH, and flow rate). The new API source shows an intensity gain of a factor 2.2 compared with ESI considering all conditions on all compounds tested. Finally, some hypotheses on the ionization mechanism, similarities, and differences with ESI, are discussed. .

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2015

Observation and numerical modeling of chromospheric evaporation during the impulsive phase of a solar flare

We have studied the chromospheric evaporation flow during the impulsive phase of the flare by using the Hinode/EUV Imaging Spectrometer observation and 1D hydrodynamic numerical simulation coupled to the time-dependent ionization. The observation clearly shows that the strong redshift can be observed at the base of the flaring loop only during the impulsive phase. We performed two different numerical simulations to reproduce the strong downflows in FeXII and FeXV during the impulsive phase. By changing the thermal conduction coefficient, we carried out the numerical calculation of chromospheric evaporation in the thermal conduction dominant regime (conductivity coefficient κ{sub 0} = classical value) and the enthalpy flux dominant regime (κ{sub 0} = 0.1 × classical value). The chromospheric evaporation calculation in the enthalpy flux dominant regime could reproduce the strong redshift at the base of the flare during the impulsive phase. This result might indicate that the thermal conduction can be strongly suppressed in some cases of flare. We also find that time-dependent ionization effect is important to reproduce the strong downflows in Fe XII and Fe XV.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2014

Noiseless Vlasov–Poisson simulations with linearly transformed particles

We introduce a deterministic discrete-particle simulation approach, the Linearly-Transformed Particle-In-Cell (LTPIC) method, that employs linear deformations of the particles to reduce the noise traditionally associated with particle schemes. Formally, transforming the particles is justified by local first order expansions of the characteristic flow in phase space. In practice the method amounts of using deformation matrices within the particle shape functions; these matrices are updated via local evaluations of the forward numerical flow. Because it is necessary to periodically remap the particles on a regular grid to avoid excessively deforming their shapes, the method can be seen as a development of Denavit's Forward Semi-Lagrangian (FSL) scheme (Denavit, 1972 [8]). However, it has recently been established (Campos Pinto, 2012 [20]) that the underlying Linearly-Transformed Particle scheme converges for abstract transport problems, with no need to remap the particles; deforming the particles can thus be seen as a way to significantly lower the remapping frequency needed in the FSL schemes, and hence the associated numerical diffusion. To couple the method with electrostatic field solvers, two specific charge deposition schemes are examined, and their performance compared with that of the standard deposition method. Finally, numerical 1d1v simulations involving benchmark test cases and halo formation in an initially mismatched thermal sheet beam demonstrate some advantages of our LTPIC scheme over the classical PIC and FSL methods. Benchmarked test cases also indicate that, for numerical choices involving similar computational effort, the LTPIC method is capable of accuracy comparable to or exceeding that of state-of-the-art, high-resolution Vlasov schemes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2013

Heterogeneous scalable framework for multiphase flows

Two categories of challenges confront the developer of computational spray models: those related to the computation and those related to the physics. Regarding the computation, the trend towards heterogeneous, multi- and many-core platforms will require considerable re-engineering of codes written for the current supercomputing platforms. Regarding the physics, accurate methods for transferring mass, momentum and energy from the dispersed phase onto the carrier fluid grid have so far eluded modelers. Significant challenges also lie at the intersection between these two categories. To be competitive, any physics model must be expressible in a parallel algorithm that performs well on evolving computer platforms. This work created an application based on a software architecture where the physics and software concerns are separated in a way that adds flexibility to both. The develop spray-tracking package includes an application programming interface (API) that abstracts away the platform-dependent parallelization concerns, enabling the scientific programmer to write serial code that the API resolves into parallel processes and threads of execution. The project also developed the infrastructure required to provide similar API s to other application. The API allow object-oriented Fortran applications direct interaction with Trilinos to support memory management of distributed objects in central processing units (CPU) and graphic processing units (GPU) nodes for applications using C++.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2015

Incompressible multiphase flow and encapsulation simulations using the moment-of-fluid method

A moment-of-fluid method is presented for computing solutions to incompressible multiphase flows in which the number of materials can be greater than two. In this work, the multimaterial moment-of-fluid interface representation technique is applied to simulating surface tension effects at points where three materials meet. The advection terms are solved using a directionally split cell integrated semi-Lagrangian algorithm, and the projection method is used to evaluate the pressure gradient force term. The underlying computational grid is a dynamic block-structured adaptive grid. The new method is applied to multiphase problems illustrating contact-line dynamics, triple junctions, and encapsulation in order to demonstrate its capabilities. Examples are given in two-dimensional, three-dimensional axisymmetric (R–Z), and three-dimensional (X–Y–Z) coordinate systems.

97 MATHEMATICS AND COMPUTING↗