Modeling and Simulation of Bubbling Fluidized Beds Containing Particle Mixtures
A comprehensive model is presented for mathematically describing the isothermal, non-reactive, fluid dynamics of a mixture of particles in a gas.
SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
A comprehensive model is presented for mathematically describing the isothermal, non-reactive, fluid dynamics of a mixture of particles in a gas.
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.
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.
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.
A model for high-pressure binary diffusion coefficient calculation is proposed based on considerations originating from recasting both the low pressure kinetic theory and the Stokes-Einstein infinite dilution expressions into forms consistent with corresponding states theory.
In the search for efficient energy storage battery technologies, designing stable electrolytes has been a long-standing challenge. Electrolytes based on molten salt eutectics are known for their stability with minimum parasitic reactions when compared to their widely used organic counterparts. However, the operating temperatures of these molten salt electrolyte-based batteries are dictated by the melting point of the eutectic mixtures. Design and high throughput screening of low melting temperature eutectic molten salt mixtures have been hindered by the lack of computational models. In this work, we develop thermodynamic models to predict the eutectic points of several molten salt mixtures. The framework of the COSMO-SAC model is used for the predictions and is compared with experimental data and other thermodynamic approaches. Rapid thermodynamics-based approaches, as shown in this study, can accelerate the discovery of new materials, complementing experimental techniques.
Abstract. Deep learning (DL) and machine learning (ML) are widely used in hydrological modelling, which plays a critical role in improving the accuracy of hydrological predictions. However, the trade-off between model performance and computational cost has always been a challenge for hydrologists when selecting a suitable model, particularly for probabilistic post-processing with large ensemble members. This study aims to systematically compare the quantile regression forest (QRF) model and countable mixtures of asymmetric Laplacians long short-term memory (CMAL-LSTM) model as hydrological probabilistic post-processors. Specifically, we evaluate their ability in dealing with biased streamflow simulations driven by three satellite precipitation products across 522 nested sub-basins of the Yalong River basin in China. Model performance is comprehensively assessed using a series of scoring metrics from both probabilistic and deterministic perspectives. Our results show that the QRF model and the CMAL-LSTM model are comparable in terms of probabilistic prediction, and their performances are closely related to the flow accumulation area (FAA) of the sub-basin. The QRF model outperforms the CMAL-LSTM model in most sub-basins with smaller FAA, while the CMAL-LSTM model has an undebatable advantage in sub-basins with FAA larger than 60 000 km2 in the Yalong River basin. In terms of deterministic predictions, the CMAL-LSTM model is preferred, especially when the raw streamflow is poorly simulated and used as input. However, setting aside the differences in model performance, the QRF model with 100-member quantiles demonstrates a noteworthy advantage by exhibiting a 50 % reduction in computation time compared to the CMAL-LSTM model with the same ensemble members in all experiments. As a result, this study provides insights into model selection in hydrological post-processing and the trade-offs between model performance and computational efficiency. The findings highlight the importance of considering the specific application scenario, such as the catchment size and the required accuracy level, when selecting a suitable model for hydrological post-processing.
The advanced thermal-hydraulic system code, System Analysis Module (SAM), was originally developed for the modeling of single-phase flow in advanced reactors. It has since been expanded to include a four-equation drift flux model for the modeling of two-phase flows containing a noncondensable gas. The model was expanded to support the modeling of molten salt reactor (MSR) designs in which the fuel is directly dissolved in the circulating coolant. These designs have shown that circulating gas bubbles can play an important role in the management of fission products and the operational behavior of the reactor. A drift flux model was implemented to more accurately capture the localized behavior of the void in the core and its impact on the mass transfer of fission products. A thorough assessment of the new model was performed by developing a verification and validation test suite. Verification problems were designed to test all major terms in the new governing equations. The new model converged to the correct solution at the expected order of accuracy for all verification cases. The validation cases included a wide range of flow and void conditions in different pipe geometries. Although higher void experiments show a slight underprediction of void by the drift flux model, experiments that aim to reproduce Molten Salt Reactor Experiment (MSRE) experimental conditions show good agreement with the model. The gas transport model was activated for a SAM model of the MSRE to demonstrate that it can be used in a more complex model. Finally, this gas transport model will be used along with an interfacial area transport equation being implemented in SAM for the prediction of mass transport behavior in MSR conditions.
An accurate prediction of the probability density function (PDF) of the mixture fraction is crucial to the prediction of combustion since mixing plays an important role in turbulent non-premixed and partially premixed flames. This work provides an assessment of the large-eddy simulation (LES)/filtered density function (FDF) method for the prediction of the PDF of the mixture fraction. The advantage of the LES/FDF method is that it provides the full predictions of the statistical distribution of scalars including the mixture fraction. The predictive accuracy of the method for the PDF is yet to be fully validated. Assessing the prediction of the PDF of the mixture fraction, a conserved scalar, is an important starting point. The Sydney/Sandia inhomogeneous inlet jet flame is used as a test case. A quick comparison shows that the LES/FDF predicted PDF shapes of the mixture fraction deviate significantly from the commonly presumed Beta-PDF as well as from the experimental data in the flame. Here, to examine the source of the discrepancy, we clarify the different PDF definitions used in the comparison among the predictions, measurements, and the presumed shape PDFs. The discrepancy observed from the comparison is largely reconciled by clarifying the difference between the PDFs that are examined. The PDF of the resolved mixture fraction is shown to be close to the Beta-PDF in both the measurements and predictions, while the PDF directly deduced from the LES/FDF particles deviates significantly from the Beta-PDF. A multimodal PDF analysis and a pseudo convergence analysis are conducted to provide plausible evidence to support the predicted multimodal PDF shapes. The sub-filter scale FDF is shown to be close to the Beta-PDF too through the construction of a synthesized PDF, which supports the common presumed Beta-PDF assumption used in the presumed PDF methods when combined with LES.
Here, we report a mechano-chemical model for anodic degradation during fast-charging of nickel-manganese-cobalt (NMC)/graphite (C) cell due to SEI growth, lithium plating/stripping, dead lithium storage, and film fracture of composite SEI and plated lithium film. Degradation of the battery is analyzed for a range of charging rates from 1 to 6 C-rates, and the influence of plating mechanisms – lithium plating and dead lithium deposition and recovery during stripping – on the film resistance of the anode are accounted for in the model. Dynamic evolution of the interfacial properties is modeled using rule-of-mixture approach. Model predictions of plating associated stress fields are used to compute critical energy release rate for film cracking. The results indicate an increased tendency of fracture for thinner SEI film with lithium plating at higher charging rates. The process of reforming the cracked film absorbs a significant portion of the electrode current thereby reducing the cell capacity and plating efficiency. The mechano-chemical model provides an extensive analytical framework for understanding the synergistic coupling of anodic degradation mechanisms, prognosticating conditions of SEI failure, and evaluating the capacity fade and efficiency of lithium-ion battery.
Abstract In recent years, hydrogen-carrying compounds have accrued interest as an alternative to traditional fossil fuels due to their function as zero-emission fuels. As such, there is interest in investigating hydrogen-carrying compounds to improve understanding of the fuels’ characteristics for use in high-pressure systems. In the current study, the oxidation of ammonia/natural gas/hydrogen mixtures was carried out to study CO formation profiles as well as the ignition delay times behind reflected shock waves in order to refine chemical kinetic models. Experiments were carried out in the University of Central Florida’s shock tube facility by utilizing chemiluminescence to obtain OH* emission and laser absorption spectroscopy to obtain CO profiles. Experimental results were then compared with the GRI 3.0 mechanism, as well as the proprietary UCF 2022 mechanism utilizing CHEMKIN-Pro software. In general, both models were able to capture the trend in autoignition delay times and CO time histories for natural gas and ammonia mixtures. However, for ammonia-hydrogen mixtures, GRI 3.0 failed to predict ignition delay times whereas the UCF 2022 mechanism was able to capture the IDTs within the uncertainty limits of the experiments. A sensitivity analysis was conducted for different mixtures to understand the important reactions at the experimental conditions. Finally, a reaction pathway analysis was carried out to understand important ammonia decomposition pathways in the presence of hydrogen and natural gas.
Current plans for the River Protection Project envision starting to vitrifying low-activity waste (LAW) by 2023 using a Direct Feed Low-Activity Waste (DFLAW) approach and subsequently using a full-pretreatment approach. The Hanford Tank Waste Treatment and Immobilization Plant (WTP) LAW Facility will be operated and controlled using a LAW glass formulation algorithm (GFA), which requires several inputs based on research and development results. LAW glass property-composition models for several product quality and processing properties are key inputs for the LAW GFA. It is envisioned that the preliminary LAW GFA discussed by Kim and Vienna (2012) will be used for commissioning and initial radioactive operations of the WTP LAW Facility under Bechtel National, Inc. using the DFLAW approach. Then, an updated LAW GFA will be developed for implementation by the WTP operating contractor that takes over after WTP LAW Facility commissioning. This report documents the enhanced LAW glass property-composition models developed for use in the updated LAW GFA. The properties for which models were developed include Product Consistency Test (PCT) response, Vapor Hydration Test (VHT) response, viscosity at 1150 °C, electrical conductivity at 1150 °C, melter SO 3 tolerance at 1150 °C, and K-3 refractory corrosion at 1208 °C. Table S.1 lists the tables in this report that contain the recommended models for each of these properties. The model types recommended include partial quadratic mixture (PQM) models for viscosity, electrical conductivity, melter SO 3 tolerance and K-3 corrosion, bias corrected PQM model (bcPQM) for PCT, and logistic PQM model for VHT. The fits of model and validation subsets were found to be well predicted by the recommended models.
Current plans for the River Protection Project envision starting to vitrifying low-activity waste (LAW) by 2023 using a Direct Feed Low-Activity Waste (DFLAW) approach and subsequently using a full-pretreatment approach. The Hanford Tank Waste Treatment and Immobilization Plant (WTP) LAW Facility will be operated and controlled using a LAW glass formulation algorithm (GFA), which requires several inputs based on research and development results. LAW glass property-composition models for several product quality and processing properties are key inputs for the LAW GFA. It is envisioned that the preliminary LAW GFA discussed by Kim and Vienna (2012) will be used for commissioning and initial radioactive operations of the WTP LAW Facility under Bechtel National, Inc. using the DFLAW approach. Then, an updated LAW GFA will be developed for implementation by the WTP operating contractor that takes over after WTP LAW Facility commissioning. This report documents the enhanced LAW glass property-composition models developed for use in the updated LAW GFA. The properties for which models were developed include Product Consistency Test (PCT) response, Vapor Hydration Test (VHT) response, viscosity at 1150 °C, electrical conductivity at 1150 °C, melter SO3 tolerance at 1150 °C, and K-3 refractory corrosion at 1208 °C. Table S.1 lists the tables in this report that contain the recommended models for each of these properties. The model types recommended include partial quadratic mixture (PQM) models for viscosity, electrical conductivity, melter SO 3 tolerance and K-3 corrosion, bias corrected PQM model (bcPQM) for PCT, and logistic PQM model for VHT. The fits of model and validation subsets were found to be well predicted by the recommended models.
The dielectric properties of vegetation material--primarily agricultural plants--as a function of moisture content and microwave frequency was measured and used to develop dielectric mixing models for the vegetation-water mixture. A model for the loss factor of a vegetation canopy was also developed. During the first phase of this investigation, three waveguide transmission systems covering from 1 to 2-GHz, from 3.5 to 6.5-GHz, and from 7.5 to 8.5-GHz bands were constructed and calibrated. By measuring the magnitude and phase of the field transmission coefficient of a given sample, it was possible to calculate the real and imaginary parts of the complex dielectric constant of the sample. Measurements were made for numerous samples of leaves and stalks of wheat and corn, and for wheat heads. Also, dielectric measurements were made of the liquid included in the vegetation material after extraction by mechanical means. The propagation loss is more than an order of magnitude greater than had previously been assumed. Various types of dielectric mixing models were investigated in terms of the available data, and a propagation model was developed and evaluated against direct canopy attenuation measurements. The canopy measurements were made by transmitting a signal from a radar antenna mounted atop a truck-mounted boom, and using a small antenna mounted on a rail beneath the canopy to receive it.
Oxidation of ATJ/n-heptane blends was studied over a wide range of test conditions, using single pulse shock tubes. Test conditions were designed to study the effect of pressure (4 and 50 bar), and fuel loading (~100-1400 ppm) on the oxidation of the blends across a wide range of temperatures (800-1300 K). These effects were observed by measuring concentrations of intermediate species formed using in-line Gas Chromatography (GC) and GCxGC TOF-MS. Results showed that increasing the initial fuel load does not have a significant impact on the results. However, increasing the pressure shifts the mixture reactivity to a lower temperature by about 150 K and causes the fuel to be oxidized instead of decomposing pyrolytically as at lower pressures. Additional experiments for pure ATJ and pure n-heptane were performed at conditions matching the ATJ/n-heptane 50 bar experiments to analyze the differences between the pure and mixed fuels. Speciation measurements were compared against predictions from a detailed kinetic model. The ability of the kinetic model to capture the effects of varying experimental test conditions on the evolution of intermediate species is discussed and kinetic analyses have been conducted to identify the important reaction pathways.
The effects of grain size and shape on the reflectance spectra of mineral mixtures are investigated to improve a reflectance model called the isograin model, whose prototype was proposed by M. Kinoshita in 1985. The sample powder was assumed to consist of an infinite number of layers, each of which has the same thickness with the grain size d.
We present direct numerical simulations of a three-layer Rayleigh–Taylor instability (RTI) problem with a configuration based on the experiments of Suchandra & Ranjan ( J. Fluid Mech. , vol. 974, 2023, A35) and Jacobs & Dalziel ( J. Fluid Mech. , vol. 542, 2005, pp. 251–279). The problem consists of a layer of light fluid between two layers of heavy fluid with an Atwood number of 0.3. These simulations are first validated through comparison with available experimental data. The validated simulations are then utilized to analyse statistics in this three-component flow. First, length scales are examined utilizing spectra and two-point spatial correlations of velocity and species concentration fluctuations. Next, joint probability density functions (p.d.f.s) of species concentration are compared against several model p.d.f.s representing generalizations of the bivariate beta distribution. Notably, the joint p.d.f.s do not appear to be accurately described by a Dirichlet distribution, indicating the marginal distributions do not conform to a beta distribution. Finally, similarity of the present configuration to three-component mixing found in inertial confinement fusion (ICF) applications is exploited to develop and validate an improved model for the impact of multicomponent mixing on thermonuclear (TN) reaction rates. A single time instant from the present simulations is chosen for a TN burn calculation under the hypothetical assumption of ICF materials and temperatures. Total TN output from this second calculation is then compared against the prediction of the improved model. The new model is found to accurately predict TN reaction rates in both premixed and non-premixed configurations.