Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multiphase 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.

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

At least 1,153 records · Page 64
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Probing the Outflowing Multiphase Gas ∼1 kpc below the Galactic Center

Comparison of interstellar medium (ISM) absorption in the UV spectrum of LS 4825, a B1 Ib−II star d = 21 ± 5 kpc from the Sun toward l = 1.°67 and b = −6.°63, with ISM absorption toward an aligned foreground star at d < 7.0 ± 1.7 kpc, allows us to isolate and study gas associated with the Milky Way nuclear wind. Spectra from the Space Telescope Imaging Spectrograph show low-ionization absorption out to d < 7 kpc (e.g., O i, C ii, Mg ii, Si ii, Fe ii, S ii) only between 0 and 40 km s{sup −1}, while absorption at d > 7 kpc, ∼1 kpc below the Galactic plane, is complex and spans −290 to +94 km s{sup −1}. The intermediate and high ions Si iii, C iv, Si iv, and N v show extremely strong absorption with multiple components from −283 to 107 km s{sup −1}, implying that the ISM ∼1 kpc below the Galactic center has a substantial reservoir of plasma and more gas containing C iv and N v than in the Carina OB1 association at z = 0 kpc. Abundances and physical conditions are presented for many absorption components. The high ion absorption traces cooling transition temperature plasma probably driven by the outflowing hot gas, while the extraordinarily large thermal pressure, p / k ∼ 10{sup 5} cm{sup −3} K{sup −1}, in an absorption component at −114 km s{sup −1} probably arises from the ram pressure of the outflowing hot gas. The observations are consistent with a flow whose ionization structure in the high ions can be understood through a combination of nonequilibrium radiative cooling and turbulent mixing.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRSDate not supplied

Comprehensive Digital Twin of a Macro-fluidic Electrochemical Reactor to Optimize the Electrochemical-based Recovery of O2 from Metabolic CO2

Future long-duration missions will require a sustainable and efficient system capable of yielding a minimum of 75% O 2 recovery from metabolic CO 2 to achieve self-sufficiency for long space missions beyond Earth's low orbit. A Macro-fluidic Electrochemical Reactor (MFECR) development effort to electrochemically recover O 2 from CO 2 is underway at NASA Marshall Space Flight Center (MSFC) to increase current O 2 recovery efficiency and reduce air revitalization (AR) system complexity at the International Space Station (ISS) habitat and future long missions. The authors have developed and deployed a digital twin of an actual single cell of the MFECR via a comprehensive 3D multiphysics model that thoroughly replicates the exact configuration and fluid/material domains of the MFECR. This model's electrochemical physics consists of multicomponent-multiphase electrochemical-driven reactions leading to CO 2 conversion to C 2 H 4 and CO along with the formation of H 2 on the cathode in parallel with the generation of O 2 and H 2 O on the anode. This electrochemical model is coupled with all the physics phenomena involved in the process, including but not limited to fluid and non-ideal mass transfer of reactant and product species in free and porous media, convective/conduction/radiative heat transfer, and conduction of DC electrical current with Joule heating generation. The digital twin has proved to be an essential tool for performing different qualitative studies, including the effect of reducing the height (increasing the respect ratio) of the serpentine walls, leading to a further redesign of the EDU (the newly redesigned MFECR EDU has been fabricated and is expected to be used in future tests), the evaluation of setting different MFECRs connected in series, and the assessment of feeding air directly to the MFECR skipping the preceding metabolic CO 2 separation from air. The MFECR's test stand is fully automated and equipped with several inline measurements (flow, pressure, temperature, pH, component concentration) systems on all six MFECR's IO streams, allowing reliable experimental validation of the model, parametric determination of all electrochemical reactions, and process optimization.

Jesus A Dominguez↗
DOE OSTI2016

High Fidelity BWR Fuel Simulations

This report describes the Consortium for Advanced Simulation of Light Water Reactors (CASL) work conducted for completion of the Thermal Hydraulics Methods (THM) Level 3 milestone THM.CFD.P13.03: High Fidelity BWR Fuel Simulation. High fidelity computational fluid dynamics (CFD) simulation for Boiling Water Reactor (BWR) was conducted to investigate the applicability and robustness performance of BWR closures. As a preliminary study, a CFD model with simplified Ferrule spacer grid geometry of NUPEC BWR Full-size Fine-mesh Bundle Test (BFBT) benchmark has been implemented. Performance of multiphase segregated solver with baseline boiling closures has been evaluated. Although the mean values of void fraction and exit quality of CFD result for BFBT case 4101-61 agreed with experimental data, the local void distribution was not predicted accurately. The mesh quality was one of the critical factors to obtain converged result. The stability and robustness of the simulation was mainly affected by the mesh quality, combination of BWR closure models. In addition, the CFD modeling of fully-detailed spacer grid geometry with mixing vane is necessary for improving the accuracy of CFD simulation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRSDate not supplied

A CFD Framework for Film Boiling and Quenching in Cryogenic Chilldown

Tests of cryogenic chilldown in a vertical stainless steel tube are simulated using computational fluid dynamics (CFD). Subcooled liquid nitrogen enters the room-temperature tube at one end, resulting in film boiling and eventual quenching of the tube walls. A CFD-based modeling framework employing the mixture multiphase model has been developed to simulate this chilldown process and applied to validate it against experimental data. Extensive sensitivity studies have been performed to create a consistent modeling framework, and a two test conditions are simulated to test the modeling capability. The importance of the rewetting temperature and its relation to the foam limit temperature is discussed and is considered in the CFD model. The presented results demonstrate a qualitatively consistent representation of the film boiling and quenching regime. However, the quantitative accuracy of the simulations varies with the test conditions. Simulated vapor volume fraction distributions are also presented, although the current lack of experimental volume fraction data does not permit a comparison.

Computational Fluid Dynamics↗
NASA NTRSDate not supplied

A CFD Framework for Film Boiling and Quenching in Cryogenic Chilldown

Tests of cryogenic chilldown in a vertical stainless-steel tube are simulated using computational fluid dynamics (CFD). Subcooled liquid nitrogen enters a room-temperature tube at one end, resulting in film boiling and eventual quenching of the tube walls. A CFDbased modeling framework employing the mixture multiphase model was developed to simulate this chilldown process and was applied to validate it against experimental data. Extensive sensitivity studies have been performed to create a consistent modeling framework, and a range of test conditions are simulated to assess the modeling capability. The importance of the rewetting temperature and its relation to the foam limit temperature is discussed and considered in the CFD framework. The presented results demonstrate a qualitatively consistent representation of the film boiling and quenching regime. However, the quantitative accuracy of the simulations varies with the test conditions. In addition, simulated distributions of the vapor volume fraction are presented, but the current lack of experimental volume fraction data prevents comparisons.

Computational Fluid Dynamics↗
NASA NTRSDate not supplied

A CFD Framework for Film Boiling and Quenching in Cryogenic Chilldown

Tests of cryogenic chilldown in a vertical stainless-steel tube are simulated using computational fluid dynamics (CFD). Subcooled liquid nitrogen enters a room-temperature tube at one end, resulting in film boiling and eventual quenching of the tube walls. A CFDbased modeling framework employing the mixture multiphase model was developed to simulate this chilldown process and was applied to validate it against experimental data. Extensive sensitivity studies have been performed to create a consistent modeling framework, and a range of test conditions are simulated to assess the modeling capability. The importance of the rewetting temperature and its relation to the foam limit temperature is discussed and considered in the CFD framework. The presented results demonstrate a qualitatively consistent representation of the film boiling and quenching regime. However, the quantitative accuracy of the simulations varies with the test conditions. In addition, simulated distributions of the vapor volume fraction are presented, but the current lack of experimental volume fraction data prevents comparisons.

Computational Fluid Dynamics↗
DOE OSTI2013

Mass transfer effects in a gasification riser

In the development of multiphase reacting computational fluid dynamics (CFD) codes, a number of simplifications were incorporated into the codes and models. One of these simplifications was the use of a simplistic mass transfer correlation for the faster reactions and omission of mass transfer effects completely on the moderate speed and slow speed reactions such as those in a fluidized bed gasifier. Another problem that has propagated is that the mass transfer correlation used in the codes is not universal and is being used far from its developed bubbling fluidized bed regime when applied to circulating fluidized bed (CFB) riser reactors. These problems are true for the major CFD codes. To alleviate this problem, a mechanistic based mass transfer coefficient algorithm has been developed based upon an earlier work by Breault et al. This fundamental approach uses the local hydrodynamics to predict a local, time varying mass transfer coefficient. The predicted mass transfer coefficients and the corresponding Sherwood numbers agree well with literature data and are typically about an order of magnitude lower than the correlation noted above. The incorporation of the new mass transfer model gives the expected behavior for all the gasification reactions evaluated in the paper. At the expected and typical design values for the solid flow rate in a CFB riser gasifier an ANOVA analysis has shown the predictions from the new code to be significantly different from the original code predictions. The new algorithm should be used such that the conversions are not over predicted. Additionally, its behaviors with changes in solid flow rate are consistent with the changes in the hydrodynamics.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2020

Phase distribution in Fischer-Tropsch mimicked slurry bubble column via computed tomography

This study aims to investigate the phases distribution and hydrodynamics in an air-Therminol LT-glass beads slurry bubble column reactor in bubbly and churn-turbulent flow regimes by implementing a gamma-ray computed tomography. It is found that gas holdup increases with increasing superficial gas velocity and operating pressure while increasing solids loading has an opposite effect on gas holdup measurements. It is noticed that the pressure effect is more significant at high solids loading and the solids loading is relatively more substantial at lower pressure. Solids loading has the most significant a stronger effect on the solids holdup profile than on the gas holdup profile. Current results were compared with the Sedimentation-Dispersion Model (SDM) predictions. Overall, the findings of this study advance the knowledge of the true behavior of phases distribution in FT mimicked slurry bubble columns which are pertinent in designing such reactors.

42 ENGINEERING↗
DOE OSTI2016

The effect of impeller type on silica sol formation in laboratory scale agitated tank

The multiphase polymerization reaction of the silica sol formation produced from silicic acid and potassium hydroxide solutions in laboratory scale agitated tank was studied. The reactor is equipped with four segmental baffle and top entering impeller. The inside diameter of reactor is 9 cm, the baffle width is 0.9 cm, and the impeller position is 3 cm from tank bottom. The diameter of standard six blades Rushton and three blades marine propeller impellers are 5 cm. The silicic acid solution was made from 0.2 volume fraction of water glass (sodium silicate) solution in which the sodium ion was exchanged by hydrogen ion from cation resin. The reactor initially filled with 286 ml silicic acid solution was operated in semi batch mode and the temperature was kept constant in 60 °C. The 3 ml/minute of 1 M potassium hydroxide solution was added into stirred tank and the solution was stirred. The impeller rotational speed was varied from 100 until 700 rpm. This titration was stopped if the solution in stirred tank had reached the pH of 10-The morphology of the silica particles in the silica sol product was analyzed by Scanning Electron Microscope (SEM). The size of silica particles in silica sol was measured based on the SEM image. The silica particle obtained in this research was amorphous particle and the shape was roughly cylinder. The flow field generated by different impeller gave significant effect on particle size and shape. The smallest geometric mean of length and diameter of particle (4.92 µm and 2.42 µm, respectively) was generated in reactor with marine propeller at 600 rpm. The reactor with Rushton impeller produced particle which the geometric mean of length and diameter of particle was 4.85 µm and 2.36 µm, respectively, at 150 rpm.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

Capillary pressure-saturation relations in quartz and carbonate sands: Limitations for correlating capillary and wettability influences on air, oil, and supercritical CO 2 trapping

Capillary pressure (P c )-saturation (S w ) relations are essential for predicting equilibrium and flow of immiscible fluid pairs in soils and deeper geologic formations. In systems that are difficult to measure, behavior is often estimated based on capillary scaling of easily measured P c –S w relations (e.g., air-water, and oil-water), yet the reliability of such approximations needs to be examined. In this study, 17 sets of brine drainage and imbibition curves were measured with air-brine, decane-brine, and supercritical (sc) CO 2 -brine in homogeneous quartz and carbonate sands, using porous plate systems under ambient (0.1 MPa, 23°C) and reservoir (12.0 MPa, 45°C) conditions. Comparisons between these measurements showed significant differences in residual nonwetting phase saturation, S nw,r . Through applying capillary scaling, changes in interfacial properties were indicated, particularly wettability. With respect to the residual trapping of the nonwetting phases, S nwr, CO2 > S nwr, decane > S nwr, air . Decane-brine and scCO 2 -brine Pc–Sw curves deviated significantly from predictions assuming hydrophilic interactions. Moreover, neither the scaled capillary behavior nor S nw,r for scCO 2 -brine were well represented by decane-brine, apparently because of differences in wettability and viscosities, indicating limitations for using decane (and other organic liquids) as a surrogate fluid in studies intended to apply to geological carbon sequestration. Thus, challenges remain in applying scaling for predicting capillary trapping and multiphase displacement processes across such diverse fields as vadose zone hydrology, enhanced oil recovery, and geologic carbon sequestration.

58 GEOSCIENCES↗
DOE OSTI2021

Saturation Carrying Capacity for Group A Particles in a Circulating Fluidized Bed

Empirical models continue to play a significant role in the design process of multiphase chemical reactors, particularly riser reactors in circulating fluidized bed (CFB) processes. It is imperative that accurate, industrial relevant correlations are developed to aid these design efforts. Using poor correlations could result in startup issues and significant redesign work. In this work, a new correlation is proposed to predict the saturation carrying capacity of Geldart Group A particles. This new correlation improves upon the currently available correlations for these materials and covers a broad range of Geldart Group A particles (particle diameters from 52 to 70 µm, and Archimedes numbers ranging from 5 to 20), superficial gas velocities (1 to 4 m/s), and riser diameters (0.066 to 0.3048 m). The new correlation has an Absolute Average Percent Deviation of only 17.6%, making it the most accurate correlation for Geldart Group A particles in the current literature.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRSDate not supplied

Modeling Autogenous Pressurization and Draining of a Cryogenic Storage Tank in Normal Gravity

A two-phase CFD model for autogenous pressurization and draining of a cryogenic storage tank is presented using both the Sharp Interface and Volume-Of-Fluid (VOF) approaches for capturing the front and the associated interfacial heat, mass and momentum transfer between the liquid and vapor regions. Both models are validated against data provided by the Cryogenic Propellant Storage and Transfer (CPST) Engineering Development Unit (EDU) experiment. The results of the autogenous pressurization are presented first, focusing on the phase change and turbulence effects on the tank pressure and temperature predictions. Both the Sharp Interface (SI-CFD) and VOF (VOF-CFD) multiphase models predict tank pressure during pressurization within 3% of the measured values. The sensitivity of key physical and numerical parameters of the problem are tested using the Sharp Interface model. Effects of the accommodation coefficient (AC), and the computational grid structure are studied. The second part of this paper is devoted to validating the SI-VOF model, with an enhanced capability of moving the liquid-vapor interface, against the draining segment of the EDU experiment. The VOF-CFD model was also used to simulate tank draining and its results were compared with the results of the Sharp Interface model with the moving interface. Both models predict tank pressure during draining within 3.5% of the measured values. Pressure decrease rate is underpredicted by both models during the first 100 seconds of draining but matches the experimental rate for the rest of the simulation.

Computational Fluid Dynamics↗
NASA NTRSDate not supplied

Modeling Autogenous Pressurization and Draining of a Cryogenic Storage Tank in Normal Gravity

A two-phase CFD model for autogenous pressurization and draining of a cryogenic storage tank is presented using both the Sharp Interface and Volume-Of-Fluid (VOF) approaches for capturing the front and the associated interfacial heat, mass and momentum transfer between the liquid and vapor regions. Both models are validated against data provided by the Cryogenic Propellant Storage and Transfer (CPST) Engineering Development Unit (EDU) experiment1. The results of the autogenous pressurization are presented first, focusing on the phase change and turbulence effects on the tank pressure and temperature predictions. Both the Sharp Interface (SI-CFD) and VOF (VOF-CFD) multiphase models predict tank pressure during pressurization within 3% of the measured values. The sensitivity of key physical and numerical parameters of the problem are tested using the Sharp Interface model. Effects of the accommodation coefficient (AC), and the computational grid structure are studied. The second part of this paper is devoted to validating the SI-VOF model, with an enhanced capability of moving the liquid-vapor interface, against the draining segment of the EDU experiment. The VOF-CFD model was also used to simulate tank draining and its results were compared with the results of the Sharp Interface model with the moving interface. Both models predict tank pressure during draining within 3.5% of the measured values. Pressure decrease rate is underpredicted by both models during the first 100 seconds of draining but matches the experimental rate for the rest of the simulation.

Computational Fluid Dynamics↗
DOE OSTI2013

Multiphasic analysis of the temporal development of the distal gut microbiota in patients following ileal pouch anal anastomosis

The indigenous gut microbiota are thought to play a crucial role in the development and maintenance of the abnormal inflammatory responses that are the hallmark of inflammatory bowel disease. Direct tests of the role of the gut microbiome in these disorders are typically limited by the fact that sampling of the microbiota generally occurs once disease has become manifest. This limitation could potentially be circumvented by studying patients who undergo total proctocolectomy with ileal pouch anal anastomosis (IPAA) for the definitive treatment of ulcerative colitis. A subset of patients who undergo IPAA develops an inflammatory condition known as pouchitis, which is thought to mirror the pathogenesis of ulcerative colitis. Following the development of the microbiome of the pouch would allow characterization of the microbial community that predates the development of overt disease. We monitored the development of the pouch microbiota in four patients who underwent IPAA. Mucosal and luminal samples were obtained prior to takedown of the diverting ileostomy and compared to samples obtained 2, 4 and 8 weeks after intestinal continuity had been restored. Through the combined analysis of 16S rRNA-encoding gene amplicons, targeted 16S amplification and microbial cultivation, we observed major changes in structure and function of the pouch microbiota following ileostomy. There is a relative increase in anaerobic microorganisms with the capacity for fermentation of complex carbohydrates, which corresponds to the physical stasis of intestinal contents in the ileal pouch. Compared to the microbiome structure encountered in the colonic mucosa of healthy individuals, the pouch microbial community in three of the four individuals was quite distinct. In the fourth patient, a community that was much like that seen in a healthy colon was established, and this patient also had the most benign clinical course of the four patients, without the development of pouchitis 2 years after IPAA. The microbiota that inhabit the ileal-anal pouch of patients who undergo IPAA for treatment of ulcerative colitis demonstrate significant structural and functional changes related to the restoration of fecal flow. Our preliminary results suggest once the pouch has assumed the physiologic role previously played by the intact colon, the precise structure and function of the pouch microbiome, relative to a normal colonic microbiota, will determine if there is establishment of a stable, healthy mucosal environment or the reinitiation of the pathogenic cascade that results in intestinal inflammation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2015

Non-intrusive, high-resolution, real-time, two-dimensional imaging of multiphase materials using acoustic array sensors

Two parallel multi-element ultrasonic acoustic arrays combined with sets of focal laws for acoustic signal generation and a classical tomographic inversion algorithm are used to generate real-time two-dimensional micro seismic acoustic images of multiphase materials. Proof of concept and calibration measurements were performed for single phase and two phase liquids, uniform polyvinyl chloride (PVC) plates, and aluminum cylinders imbedded in PVC plates. Measurement artefacts, arising from the limited range of viewing angles, and the compromise between data acquisition rate and image quality are discussed. The angle range of scanning and the image resolution were varied, and the effects on the quality of the reproduction of the speed of sound profiles of model solids and liquids with known geometries and compositions were analysed in detail. The best image quality results were obtained for a scanning angle range of [−35°, 35°] at a step size of 2.5° post processed to generate images on a 40 μm square grid. The data acquisition time for high quality images with a 30 mm × 40 mm view field is 10 min. Representation of two-phase solids with large differences in speed of sound between phases and where one phase is dispersed in the form of macroscopic objects (greater than 1 mm in diameter) proved to be the most difficult to image accurately. Liquid-liquid and liquid-vapor phase boundaries, in micro porous solids by contrast, were more readily defined. Displacement of air by water and water by heptane in natural porous limestone provides illustrative kinetic examples. Measurement results with these realistic cases demonstrate the feasibility of the technique to monitor in real time and on the micrometer length scale local composition and flow of organic liquids in inorganic porous media, one of many envisioned engineering applications. Improvement of data acquisition rate is an area for future collaborative study.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2018

Gen-0 BWR Modeling Application and Validation (L3 Milestone THM.P18.01)

This milestone report describes the Consortium for Advanced Simulation of Light Water Reactors (CASL) work conducted for completion of the thermal hydraulics Methods (THM) Level 3 milestone L3:THM.P18.01 Gen-0 BWR Modeling Application and Validation. The extension of existing boiling closure designed for Pressurized Water Reactors (PWRs) to the high void fraction condition of Boiling Water Reactors (BWRs) is one of the challenges for Multiphase Computational Fluid Dynamics (M-CFD) applications. In the efforts aimed at developing the advanced boiling closure for BWRs, the baseline closure model, so-called ‘Zero’ closure, has been developed by leveraging the experience gained through the Numerical Nuclear Reactor project. This baseline closure model is physics-based, simplest and robust, but it does not contain any explicit descriptions of the complex high void fraction mechanisms. This is a continuation of the CASL thermal-hydraulic study to develop and assess the CFD-BWR models. In FY-17, the applicability has been evaluated against the international OECD/NRC BFBT benchmark data. In the present work, we assessed the capability of CFD-BWR model for void drift phenomena in the modern BWR fuel designs which make use of both full-length and part-length rods.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2016

An adaptive lattice Boltzmann scheme for modeling two-fluid-phase flow in porous medium systems

In this paper, we formulate a multiple-relaxation-time (MRT) lattice-Boltzmann method (LBM) to simulate two-fluid-phase flow in porous medium systems. The MRT LBM is applied to simulate the displacement of a wetting fluid by a nonwetting fluid in a system corresponding to a microfluidic cell. Analysis of the simulation shows widely varying time scales for the dynamics of fluid pressures, fluid saturations, and interfacial curvatures that are typical characteristics of such systems. Displacement phenomena include Haines jumps, which are relatively short duration isolated events of rapid fluid displacement driven by capillary instability. An adaptive algorithm is advanced using a level-set method to locate interfaces and estimate their rate of advancement. Because the displacement dynamics are confined to the interfacial regions for a majority of the relaxation time, the computational effort is focused on these regions. The proposed algorithm is shown to reduce computational effort by an order of magnitude, while yielding essentially identical solutions to a conventional fully coupled approach. Finally, the challenges posed by Haines jumps are also resolved by the adaptive algorithm. Possible extensions to the advanced method are discussed.

42 ENGINEERING↗
DOE OSTI2018

MFIX-DEM Phi: Performance and Capability Improvements Towards Industrial Grade Open-source DEM Framework with Integrated Uncertainty Quantification

Two major challenges hinder the effective use and adoption of multiphase computational fluid dynamics tools by the industry. The first is the need for significant computational resources, which is inversely proportional to the accuracy of solutions due to computational intensity of the algorithms. The second barrier is assessing the prediction credibility and confidence in the simulation results. In this project, a multi-tiered approach has been proposed under four broad activities to overcome these challenges while addressing all of the objectives outlined in FOA-0001238 through Phases 1 and 2 of the project. The present report consists of the results for only Phase 1, which was the funded performance period. From the start the project, all of the objectives outlined in FOA were addressed through four major activity tasks in an integrated and balanced fashion to improve adoption of MFIX suite of solvers for industrial use. The first task aimed to improve the performance of MFIX-DEM specifically targeting to acquire the peak performance on Intel Xeon and Xeon Phi based systems, which are expected to be one of the primary high-performance computing platforms both affordable and available for the industrial users in the next two to five years. However, due to a number of changes in course of the project, the scope of the performance improvements related task was significantly reduced to avoid duplicate work. Hence, more emphasis was placed on the other three tasks as discussed below. The second task aimed at physical modeling enhancements through implementation of polydispersity capability and validation of heat transfer models in MFIX. An extended verification and validation (V&V) study was performed for the new polydispersity feature implemented in MFIX-DEM both for granular and coupled gas-solid flows. The features of the polydispersity capability and results for an industrially relevant problem were disseminated through journal papers (one published and one under review at the time of writing of the final technical report). As part of the validation efforts, another industrially relevant problem of interest based on rotary drums was studied for several modes of heat transfer and results were presented in conferences. Third task was aimed towards an important and unique contribution of the project, which was to develop a unified uncertainty quantification framework by integrating MFIX-DEM with a graphical user interface (GUI) driven uncertainty quantification (UQ) engine, i.e., MFIX-GUI and PSUADE. The goal was to enable a user with only modest knowledge of statistics to effectively utilize the UQ framework offered with MFIX-DEM Phi to perform UQ analysis routinely. For Phase 1, a proof-of-concept demonstration of the proposed framework was completed and shared. Direct industry involvement was one of the key virtues of this project, which was performed through forth task. For this purpose, even at the proposal stage, the project team received strong interest in the proposed capabilities from two major corporations, which were further expanded throughout Phase 1 and a new collaboration with another major corporation from chemical industry was also initiated. The level of interest received and continued collaboration for the project during Phase 1 clearly shows the relevance and potential impact of the project for the industrial users.

42 ENGINEERING↗