Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “CFD modeling”

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 253 records · Page 14

High-Fidelity and High-Performance Computational Simulations for Rapid Design Optimization of Sulfur Thermal Energy Storage

Industrial process heating (IPH) accounts for approximately 70% of US manufacturing energy use and is primarily produced by fossil fuel combustion. Approximately 1500 TWht (approximately 60%) of IPH demand is in the temperature range of 100-300. Industrial applications in this temperature range include drying, hydrothermal processing, thermal enhanced oil recovery, food and beverage, bioethanol production, etc. Cost-effective thermal energy storage (TES) that increases the utilization of waste and renewable heat (solar, geothermal, etc.) could provide significant energy savings and reliable heat sources, decrease emissions, and increase US manufacturing competitiveness through reductions in fuel consumption. TES development has historically been dominated by technologies suitable for deployment with concentrating solar power (CSP). State-of-the-art thermal storage deployed commercially with power tower CSP plants uses a 60%/40% NaNO3/KNO3 molten salt and operates between temperatures of approximately 280 degrees Celsius and 570 degrees Celsius using a two-tank configuration. However, these nitrate salts are unsuitable for operation outside of this temperature range due to a high freezing point of approximately 220 degrees Celsius, and limits on high-temperature salt stability and corrosion resistance of containment alloys. Other materials being investigated for TES include those based on: (1) sensible energy storage (various molten salt compositions, inert solid particles, rocks or pebble beds, sulfur, water, concrete, graphite, etc.), (2) latent energy storage in materials that undergo solid-liquid phase change at relevant temperatures (organic materials for low-temperature applications, inorganic salts and/or metals for high-temperature applications), or (3) thermochemical energy storage (hydrides, hydroxides, carbonates, metal oxides, etc.). The application temperature and challenges pertaining to storage material and/or containment cost, energy density, long-term thermal and cyclic stability, and charge/discharge heat transfer effectiveness drive material selection for a given IPH or electricity generation application. Sulfur is a cheap commodity at $80/ton compared to $1100 - 1300/ton for conventional salts. When using a metric of storage cost per kWh, sulfur costs around 2-3 $/kWh. Previous sulfur TES development focused on high temperature (>600 degrees) concentrated solar power applications with sulfur encapsulated in pipes and flow of gaseous HTF (air) in the shell side. However, for lower-temperature IPH applications in the range of approximately 100-300 degrees Celsius Element 16 adopted a compact and scalable TES design with molten sulfur in the shell and HTF pipes submerged in the molten sulfur bath. The low-cost molten sulfur TES for dispatchable IPH has deployment potential for broad applications. The spatial and temporal evolution of the HTF and sulfur temperature is critical to the TES system performance, and thus detailed modeling can improve understanding of the performance and facilitate design improvements. Using high performance computing and computational fluid dynamics (CFD) a low-cost molten sulfur thermal energy storage (TES) system for industrial process heating (IPH) applications was developed. The unique challenges in CFD modeling of sulfur TES are the sharp property changes of sulfur relevant to the working temperatures. Above 159, liquid sulfur undergoes polymerization, and the viscosity of sulfur rapidly increases by several orders of magnitude between 159 degrees Celsius and 188 degrees Celsius, followed by a decrease in viscosity beyond 188 degrees Celsius due to thermal bound dissociation. In addition, various concentrations of H2S impurities can also modify sulfur viscosity. This numerical challenge is especially relevant to transient simulation of the sulfur TES charging and discharging processes as the extreme property variations limit the applicability of traditional heat transfer correlations. Transient CFD simulations including the temperature-dependent sulfur properties and geometric complexity of the TES design were used to predict the effect of natural convection during charging and discharging on the heat transfer process, sulfur temperature uniformity, charge/discharge rates, and performance of the storage devices. The CFD model was validated with experimental results for a full charge and discharge cycle. The work will show 3D and 2D simulation comparisons aimed to facilitate rapid design iterations and a machine learning based design optimization approach.

CFD↗

A multiphase flow model of water droplets dielectrophoretic-induced air dehumidification phenomena

Air humidity in indoor spaces plays a critical role in human comfort and health. Dehumidification systems are used for building humidity controls, but they can take significant energy consumption, especially in geographic locations with high outdoor humidity and warm climates. Consequently, there is a growing demand for innovative dehumidification processes that consume minimal energy. Dielectrophoretic air dehumidification represents one such promising approach. However, it has not garnered significant attention due to the absence of engineering models and simulation tools capable of evaluating its performance and limitations at large-scale airflows. A new numerical multiphase CFD model, which is also experimentally validated, is developed in a customized Reacting Foam solver based on OpenFOAM® version 9. The newly developed model seeks to decrease substantial energy consumption and lower costs by leveraging the dielectrophoretic phenomenon to regulate moisture levels in the air. The solver integrates a hybrid Eulerian-Lagrangian framework to track the droplet's trajectory and growth rate while solving the continuum equations for the moist air. An electrospray produces electrically charged droplets, which grow during their in-flight trajectories as water vapor condenses onto their surfaces. The role of electrostatic forces in promoting vapor condensation within a high-gradient electrical field is investigated, and the dielectrophoretic vapor nucleation process on charged water droplets is discussed. The CFD model was validated against results from the literature and from proof-of-concept experiments conducted by the authors, which showed a 2 % air dehumidification with a single electrospray and airflow rate of 5 cubic feet per minute. The simulation results indicated that augmenting the number of electrically charged spray droplets increased the dehumidification of the air to 25 %. The initial mean droplet diameter, the orientation of the injector and relative humidity significantly influence the assessment of dehumidification. As a result, scaling up this approach to larger airflow volumes is identified as a potential future research direction.

42 ENGINEERING↗

Multidimensional Modeling of Mixture Formation in a Hydrogen-Fueled Heavy-Duty Optical Engine With Direct Injection

Hydrogen (H 2 ), as a carbon-free fuel, is considered as one of the most promising solutions to reduce the carbon footprint of hard-to-decarbonize energy and transportation sectors. As such, hydrogen-fueled internal combustion engines (H 2 ICEs) have recently been receiving increasing attention, particularly in applications such as on-road/off-road heavy-duty transport and combined heat and power. The direct injection (DI) of gaseous hydrogen into the combustion chamber offers great potential for achieving high power density and high engine efficiency, while mitigating the risk of backfire and reducing pre-ignition. However, the numerical simulation of H 2 DI system remains a formidable challenge associated with the high computational cost of reproducing compressible supersonic flow and shocks in narrow injector passages and in near-nozzle regions. In general, there is a lack of well-established and validated practices for the modeling of high-pressure H 2 DI in large-bore engines. Here, to this end, this study focuses on computational fluid dynamics (CFD) modeling of the mixture formation process in a heavy-duty optical engine employing a medium-pressure H 2 DI system. Both large eddy simulations (LES) and Reynolds Averaged Navier–Stokes (RANS) simulations are performed and evaluated against optical data. Gaseous hydrogen is injected into the combustion chamber via a centrally located outward opening hollow-cone injector at a pressure of 40 bar. Simulations are carried out for two injection timings, namely, −120 and −60 °CA. The numerical predictions for H 2 distribution in different horizontal and vertical planes during the compression stroke are systematically compared against optical data obtained through planar laser-induced fluorescence (PLIF) measurements. Overall, the LES approach using the Dynamic Structure model is found to have good predictive capabilities for the early jet penetration in terms of length and shape, as well as the later H 2 distributions. However, the unsteady RANS approach with the renormalization group $k - ϵ$ model, which is widely used by industry to model heavy-duty ICEs, significantly underpredicts the H 2 mixing, even at similar mesh resolution to that used in LES. These results indicate that there is a need for the improvement of mixing submodels within the RANS approach when applied to H 2 DI simulations.

LES↗

Recent Improvements in Pronghorn for Advanced Reactor Modeling

Pronghorn is a thermal-hydraulics computational tool developed using the Idaho National Laboratory's Multiphysics Object-Oriented Simulation Environment (MOOSE). It is designed to support Computational Fluid Dynamics (CFD) modeling, ranging from subchannel and porous media analysis to Reynolds Averaged Navier-Stokes (RANS) turbulence modeling. As an integral part of the MOOSE-based suite of tools, Pronghorn seamlessly couples with other MOOSE-based applications to simulate a variety of physical phenomena. This article highlights recent significant enhancements to Pronghorn's CFD modeling capabilities and demonstrates their application to advanced nuclear reactor designs. The recent improvements in Pronghorn primarily focus on modifications to its turbulence modeling capabilities, near-wall corrections and numerical schemes. In terms of turbulence modeling, the two-equation $k-\epsilon$ and $k-\omega$ SST models have been implemented and validated with both equilibrium and non-equilibrium wall treatments. Additionally, corrections for wall roughness, and curvature, and wall-channeling in pebble beds have been introduced in the near-wall modeling. These developments enable more accurate simulations of advanced nuclear reactors. Two case studies are presented in this work: a pool-type Molten Chloride Reactor and a salt-cooled Pebble-Bed High Temperature Reactor. In both cases, the previous models in Pronghorn are compared with the new implementations, demonstrating the improved accuracy achieved with the updated models.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Prototype Modeling for a Light-Trapping Planar-Cavity Enclosed Particle Solar Receiver

Concentrating solar thermal (CST) systems present a promising avenue for affordable and reliable energy production. Solar receivers are key components that determine the efficiency and longevity of these systems. Particle-based solar receivers have emerged as a compelling alternative to traditional technologies, offering several advantages that address limitations in current CST systems. This is especially true as next-generation CST technologies target applications including electricity generation, thermochemical processes, and industrial process heat, many of which necessitate higher operating temperatures than current commercial molten salt systems. Molten-salt thermal energy storage (TES) systems, commonly used in CSP, face challenges related to freezing and corrosion. Particle-based TES systems, in contrast, do not experience these issues, as particles are stable at high temperatures, exceeding 1000 degrees Celsius. This capability allows for a wider range of applications, including those requiring higher temperatures for industrial processes and efficient electricity generation. A novel innovation in particle-based solar receiver technology is the light-trapping planar cavity receiver (LTPCR) configuration developed by NREL. The LTPCR design consists of small cavity-like structures using opaque planar surfaces, enabling efficient capture and absorption of solar energy. A high incident flux concentration at the cavity aperture is absorbed on the receiver walls, and subsequently transferred to particles on the inside of cavities. The particles flow through the system, forming a fluidized bed inside of the receiver panels, effectively capturing the absorbed solar heat. Air is used as a fluidizing medium in this process to enhance particle heat transfer and mixing. The effectiveness of this design lies in its ability to manage solar flux conditions and ensure high solar-to-thermal receiver efficiency. A 100-kW prototype is currently being tested at the King Saud University in Saudi Arabia to assess the receiver performance. A range of modeling analyses for the optical, thermal, and mechanical effects were conducted to assess the performance of the receiver under on-sun conditions. The solar flux resulting from the KSU heliostat field was modeled using NREL SolTrace software and produced up to 600 kW/m2 at the receiver aperture. The solar flux absorbed on the receiver walls was then used within a computational fluid dynamics (CFD) model to predict wall temperature distributions along with radiation and convection loss. A two-phase CFD model was developed for the fluidized bed of silica sand inside the receiver panels to predict local wall-to-particle heat transfer coefficients, particle temperature distributions, and outlet temperature of the particles. We have also conducted analyses to understand the thermomechanical behavior of these innovative enclosed light-trapping solar receivers optimized for particle heating. We used finite element analysis (FEA) to predict the receiver's performance using temperature distributions obtained from CFD and based on the resulting stress profiles, evaluated creep-fatigue damage with a goal of achieving a 30-year service life. Analysis showed a significant impact of the particle-to-wall heat transfer coefficients (HTCs) on receiver performance, with higher HTCs resulting in reduced stress and increased lifespan. For instance, when using Inconel 740H, increasing the HTC from 800 W/m2 K to 1400 W/m2 K increased the creep life from 4,000 hours to over 100,000 hours. This highlights the importance of understanding and optimizing heat transfer in the design of high-efficiency receivers.

14 SOLAR ENERGY↗

Experiments and CFD simulation of mass transfer and hydrodynamics in a cylindrical bubble column

Flow visualization, Particle Image Velocimetry (PIV) and dissolved oxygen charging experiments are carried out in a small-scale cylindrical air-water bubble column at different superficial gas velocities. Void fraction profile, total gas holdup, liquid phase velocity field and volumetric mass transfer coefficient are obtained. Details of each experiment are given, with a new procedure for obtaining a 2D void fraction profile from flow visualization by image processing and time averaging which could provide useful information for CFD model validation. The experimental results are analyzed and compared with simulations where good agreement is found for both local and global flow patterns. The averaged governing equations and constitutive relations used in the simulation are discussed briefly, with emphasis on the proper comparison between simulations and experiments based on the derivation of the two-fluid model. Finally, by conducting three sets of experiments related to phase distribution, velocity field and mass transfer in the same apparatus, a validated CFD model for multiphase mass transfer is established and can be used in the future design of a multi-phase reactor.

42 ENGINEERING↗

Modifications to Solar Titan-130 Combustion Systems for Efficient, High Turndown Operation

The project team of Southwest Research Institute® (SwRI®), Solar Turbines Incorporated (Solar), the Electric Power Research Institute (EPRI), the University of California, Irvine (UCI), and the Georgia Institute of Technology (Georgia Tech) investigated methods to allow higher efficiency part-load operation of a Solar Titan 130 gas turbine. The objective was to develop a low-emission combustion system capable of sustaining combustion and avoiding lean blowout during high turndown operation, which would allow the gas turbine to operate as efficiently as possible at part load. Currently, electric utility markets are beginning to experience substantial increases in renewable energy generation. Some of these renewable energy sources have highly variable output in an uncontrolled manner. In order to maintain grid stability, there is a need for power plants to ramp up power to the grid rapidly to make up for drops in renewable generation. This is often termed spinning reserve, but the size of this reserve may need to increase as renewable penetration into the electric utility market increases. Small combined heat and power (CHP) power plants provide a promising option for meeting this spinning reserve requirement. In order to operate in spinning reserve while still meeting the heat requirements for the CHP, the gas turbine needs to operate efficiently at very low loads. Efficient, high turndown operations in this engine are limited by the lean flammability limit of the premixed combustion system. This project sought enhance the lean operability range of the Titan 130 combustor. First, the project team participated in a brainstorming activity and ultimately selected two concepts to explore: fuel augmentation with hydrogen (H2) to improve the stability at lean operating conditions and modifications to the fuel nozzle to improve the emissions performance at lean operating conditions. Analytical and laboratory investigations were accomplished by UCI to investigate the efficacy of H2 addition at improving lean blow out (LBO) limits and the resulting emissions. These investigations used a variety of chemical reactor network (CRN) and CFD models, validated against laboratory data, to model the impact of H 2 and inform the experimental efforts accomplished by SwRI and Solar. Ultimately, both the CRN and CFD models yielded generally good agreement with the experimental data below a particular temperature threshold. Atmospheric tests of a full-scale T130 annular combustor were performed at SwRI facilities in San Antonio, Texas, to investigate the use of H 2 addition. For these tests, the T130 combustion system remained largely unchanged; minor modifications were performed to the fuel ducting to allow for the safe use of H 2 . The test ultimately demonstrated that the addition of H 2 to the fuel mixture significantly increased the AFR ratio at which the combustor could operate. This improvement to the LBO limit should allow for less use of compressor bleed and less throttling needed by the inlet guide vanes (IGV). This in turn could result in more efficient operation of the gas turbine at lower load points. The second modification explored in this work was a direct modification to the T130 injector. The project team hypothesized that modifications to the pilot of the T130 injector could provide lower emissions at high turn-down operations. These modifications were manufactured and explored by the team at Solar. High pressure rig tests, originally slated to occur at SwRI, were ultimately accomplished by Solar to maintain overall project budget and mitigate cost growth attributable to supply chain issues and inflation. The pressurized rig tests ultimately showed that the SwRI Project No. 18.24153 - DE-EE0008415 Page 2 Final Technical Report January 24, 2024 modifications did not significantly alter the performance of the combustion system at the high turn-down conditions; both the modified injectors and the baseline configuration exhibited elevated emissions comparted to the full-load operating condition. A final set of studies performed by EPRI investigated the benefit-cost of flexible CHP as well as a grid interconnection study for the California Independent System Operator (CAISO) grid. These studies considered: traditional CHP with no spinning reserve available for on-demand grid support, 50% flexible CHP where 50% of the machine’s capacity is consumed by on-site baseload operations while providing an additional 50% capacity for on-demand grid support, and 70% flexible CHP where 70% of capacity is consumed on-site by baseload operations and 30% is available for on-demand grid support. In all cases, the analyses showed a benefit-to-cost ratio greater than unity implying a positive net present value for all configurations. However, the traditional CHP showed the most economic benefit. These results are sensitive to several factors, many of which are not fully known and may vary over time. Thus site owners must be convinced that taking up the increased costs and risks from flexible CHP would be worth implementing. As the grid in California and across the country transition to incorporate larger renewable energy generation, flexible CHP can provide much needed operating reserves and dispatchability. Alternative fuel options, such as hydrogen blending and biofuels, may also lower carbon intensities of CHP. Flexible CHP should be examined in the evolving market to understand innovative business models, changes market rules and services, and new technologies.

20 FOSSIL-FUELED POWER PLANTS↗

Space Shuttle and Launch Pad Lift-Off Debris Transport Analysis: SRB Plume-Driven

This paper discusses the Space Shuttle Lift-Off model developed for potential Lift-Off Debris transport. A critical Lift-Off portion of the flight is defined from approximately 1.5 sec after SRB Ignition up to 'Tower Clear', where exhaust plume interactions with the Launch Pad occur. A CFD model containing the Space Shuttle and Launch Pad geometry has been constructed and executed. The CFD model works in conjunction with a debris particle transport model and a debris particle impact damage tolerance model. These models have been used to assess the effects of the Space Shuttle plumes, the wind environment, their interactions with the Launch Pad, and their ultimate effect on potential debris during Lift-Off. Emphasis in this paper is on potential debris that might be caught by the SRB plumes.

West, Jeff↗

Evaluation of Howard A. Hanson Dam Juvenile Fish Passage and Survival Study Live Fish Injury Assessment, Sensor Fish, and BioPA Modeling Tasks

The live fish injury assessment, Sensor Fish, and BioPA modeling study tasks were conducted by researchers from Pacific Northwest National Laboratory (PNNL). The four tasks were part of the larger Evaluation of Howard A. Hanson Dam (HAHD) Juvenile Fish Passage and Survival study, which had six total tasks. To achieve study objectives for each of the four tasks, field work occurred at Green Peter Dam (GPR) to evaluate the highest elevation steep slope bypass pipe, at HAHD to evaluate baseline conditions of the horseshoe tunnel, and at PNNL’s Aquatic Research Laboratory (ARL) to evaluate simulated dam passage conditions (i.e., shear forces and collision). Each of these evaluations utilized live fish injury assessment, Sensor Fish, and BioPA modeling. Live fish injury assessment and survival (tagged with and without balloon or passive integrated transponder [PIT] tags) was correlated with Sensor Fish to determine thresholds. The CFD analyses were then performed, and the computed values were compared to the corresponding measured values of Sensor Fish data. The results of the overall injury and survival of fish was also used in the validation of the CFD modeling method. Collectively, the results will aid in future modeling of fish passage at HAHD. Results from these tasks can be used by biologists, engineers, resource managers, and regional decision-makers to inform baseline conditions under current operations and the engineering design of the new FPF at HAHD. This draft report contains initial data and results from the four tasks. Table 8 1, Table 8 2, and Table 8 3, and Figure 8 1, Figure 8 2, and Figure 8 3 depict the CFD modeling findings for the GPR steep slope bypass, HAHD horseshoe tunnel, and laboratory testing. Table 8 4, Table 8 5, and Table 8 6 depict the Sensor Fish findings for the GPR steep slope bypass and HAHD horseshoe tunnel testing. The Mv values observed in the HAHD were significantly lower compared to the laboratory experiments conducted at PNNL. Currently, investigations are underway to understand the reasons for this disparity and to establish an appropriate threshold value for Mv. Survival predictions presented in the tables below should be considered preliminary and should not be used until further analyses and adjustments are completed. The next steps for modeling will include the flow regime, (i.e., density of flow regimes due to water and air mixing ) to continue to improve on the threshold value for Mv.

13 HYDRO ENERGY↗

HPC-Enabled Optimization of High Temperature Heat Exchangers (CRADA Final Report)

This project was a collaborative effort between Lawrence Livermore National Security, LLC (LLNS) as manager and operator of Lawrence Livermore National Laboratory (LLNL) and Materials Sciences, LLC, to develop a technology for design and optimization of heat exchangers using powerful desktop and laptop computers. The project was originally designated as a 12-month project, and consisted of 3# major tasks and the following 8# major deliverables: 1) CFD models of 3D heat exchangers based on existing and new geometry. 2) Validation against experimental data provided by MSC and published in the literature. 3) CFD models of 3D unit cells based on TPMS. 4) Surrogate models capable of delivering the gradients of the homogenized properties with respect to the parametrization. 5) 3D design methodology using TO algorithms. 6) Conventional reference and topology optimized designs. 7) 3D optimized designs stored in a 3D printer build format. 8) Verification of the improved performance. All of the deliverables for this project were successfully completed with two no-cost time extensions.

36 MATERIALS SCIENCE↗

HPC-Enabled Optimization of High Temperature Heat Exchangers (CRADA Final Report)

This project was a collaborative effort between Lawrence Livermore National Security, LLC (LLNS) as manager and operator of Lawrence Livermore National Laboratory (LLNL) and Materials Sciences, LLC, to develop a technology for design and optimization of heat exchangers using powerful desktop and laptop computers. The project was originally designated as a 12-month project, and consisted of 3# major tasks and the following 8# major deliverables: 1) CFD models of 3D heat exchangers based on existing and new geometry. 2) Validation against experimental data provided by MSC and published in the literature. 3) CFD models of 3D unit cells based on TPMS. 4) Surrogate models capable of delivering the gradients of the homogenized properties with respect to the parametrization. 5) 3D design methodology using TO algorithms. 6) Conventional reference and topology optimized designs. 7) 3D optimized designs stored in a 3D printer build format. 8) Verification of the improved performance. All of the deliverables for this project were successfully completed with two no-cost time extensions.

13 HYDRO ENERGY↗

ADVANCED THERMAL-HYDRAULIC MODEL OF HEAT RECOVERY STEAM GENERATORS

Steam-side oxidation at elevated temperatures in heat recovery steam generators (HRSGs) can lead to overheating of local tubes, potential oxide scale exfoliation, and damage to critical equipment. Predicting the oxide growth distribution requires detailed knowledge of the distribution of tube metal temperatures. To address these challenges, a computational fluid dynamics (CFD) model of a high-pressure superheater (HPSH) was developed within the Siemens’ STAR-CCM+ CFD commercial software. The HPSH section of the steam generator was selected for modeling because of its high oxidation rates.The model development included several steps. In the first step, an oxidation model based on a parabolic oxidation relationship was implemented. The oxidation growth predicted by the model was then verified for a full-load boiler operation case. In the second step, a three-dimensional CFD model was developed to account for the thermal and flow physics involved in HRSGs. The CFD model was then verified and validated on a simplified geometry that contained six row-tubes which mimicked the tube arrangements of an HPSH section. Numerical results compared well with published experimental results. In the third step, the oxidation model was coupled to the HRSG’s CFD model, and a comparison study was performed for two cases: (a) regular finless tubes and (b) finned tubes. Preliminary results showed that the oxide thickness predicted by the case with tube fins was higher compared to that of the finless tube case, thus demonstrating the need to model the effects of the tube fins on oxide thickness.

Panicker, Nithin↗

Numerical study of the SSME nozzle flow fields during transient operations: A comparison of the animated results with test

A computational fluid dynamics (CFD) model has been applied to study the transient flow phenomena of the nozzle and exhaust plume of the Space Shuttle Main Engine (SSME), fired at sea level. The CFD model is a time accurate, pressure based, reactive flow solver. A six-species hydrogen/oxygen equilibrium chemistry is used to describe the chemical-thermodynamics. An adaptive upwinding scheme is employed for the spatial discretization, and a predictor, multiple corrector method is used for the temporal solution. Both engine start-up and shut-down processes were simulated. The elapse time is approximately five seconds for both cases. The computed results were animated and compared with the test. The images for the animation were created with PLOT3D and FAST and then animated with ABEKAS. The hysteresis effects, and the issues of free-shock separation, restricted-shock separation and the end-effects were addressed.

Wang, Ten-See↗

Computational Fluid Dynamics Simulations to Assess Spatial Variability and Optimal Ventilation Scenarios for Biological Laboratory Exposures

A significant amount of uncertainty exists regarding potential human exposure to laboratory biomaterials and organisms in Biosafety Level 2 (BSL-2) research laboratories. Computational fluid dynamics (CFD) modeling is proposed as a way to better understand potential impacts of different combinations of biomaterials, laboratory manipulations, and exposure routes on risks to laboratory workers. Here, in this study, we use CFD models to simulate airborne concentrations of contaminants in an actual BSL-2 laboratory under different configurations. Results show that ventilation configuration, sampling location, and contaminant source location can significantly impact airborne concentrations and exposures. Depending on the source location and airflow patterns, the transient and time-integrated concentrations varied by several orders of magnitude. Contaminant plumes from sources located near a return vent (or exhaust like a fume hood or ventilated biosafety cabinet) are likely to be more contained than sources that are further from the exhaust. Having a direct flow between the source and the exhaust (through-flow condition) may reduce potential exposures to individuals outside the air flow path. Designing a BSL-2 room with ventilation and airflow patterns that maximize through-flow conditions to the return/exhaust vents and minimize dispersion and mixing throughout the room is, therefore, recommended. CFD simulations can also be used to assist in characterizing the impacts of supply and return vent locations, room layout, and source locations on spatial and temporal contaminant concentrations. In addition, proper placement of particle sensors can also be informed by CFD simulations to provide additional characterization and monitoring of potential exposures in BSL-2 facilities.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗

Effect of cold cap coverage and emissivity on the plenum temperature in a pilot-scale waste vitrification melter

Integrated computational fluid dynamics (CFD) models are being developed to model the complex physics occurring within the high-level waste melter for vitrification of legacy tank waste at the Hanford site. Here, we present a validation of the integrated CFD model by using data from two experimental runs in a pilot-scale melter. While the model uses several simplifying assumptions (such as constant heat sinks from a cooling jacket and inleakage of ambient air, steady state feed-to-batch conversion heat, and a cold cap model with a simplified shape), it closely predicts the molten glass (1150°C and 1175°C) and plenum temperatures (550°C) obtained from thermocouples during two pilot-scale tests, with an average cold cap coverage of 80%. Additional simulations were performed to explore the sensitivity of the predicted plenum temperatures to variations in cold cap coverage (fraction of melt surface covered by the glass batch) and batch emissivity. The plenum temperature was found to be in the range of 606°C when cold cap coverage decreased from 95% to 70%. Cold cap emissivity had a smaller effect, increasing the plenum temperature by as much as 179°C when cold cap emissivity increased from 0.2 to 0.8. Maintaining a high cold cap coverage without overfeeding is important for a sustained melter operation with high glass throughput. This work provides a tool for achieving that goal in terms of correlating the plenum temperature with the cold cap coverage.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Identification and Computation of Individual Propeller Acoustics of the Joby Aviation Aircraft

The individual propeller sources of the Joby Aviation aircraft are separated using an order tracking filter on measurement data and computed using high- fidelity computational fluid dynamics (CFD) to quantify the relative contributions to the total noise. Computational results are verified against experimental data for hover and steady level flight at 60 knots and demonstrate very good agreement of both overall sound pressure level magnitude and directivity. The propeller source separation for hover shows good agreement between the methods except for the inboard propellers, and it is expected that the assumptions in the CFD model do not accurately capture the flight test condition. Both methods are able to separate the higher tonal levels generated by the tail propellers due to the aerodynamic interaction. Overall, considering the discrepancies between the CFD model and flight test, the agreement in Vold-Kalman-filtered results and prediction yields confidence that the separation method is successful.

eVTOL↗

Assessing Critical Conditions for Scour Near Obstructions using Bed Shear, Particle Onset of Motion Balances, and CFD-DEM Modeling of Granular Beds

Computational Fluid Dynamics combined with a Discrete Element Method is one of the computational methods that can be used to model multiphase flows. In this method various phases, gas and liquid or solid, are present in the same computational domain. The local averaged Navier–Stokes equations determine the flow of the continuous phase fluid and are solved using the traditional CFD finite volume approach. DEM is based on a Lagrangian formulation, which solves the equations of motion, expressed in ordinary differential equations, for representative particles as they move in space and time. The interactions between the continuous fluid phase and discrete solid phase are modeled with the use of Newton’s laws of motion via drag force. The particles interact with each other and with the boundaries of the fluid continuum, and the resulting contact forces are included in the equations of motion. The properties of solid particles and boundaries are treated as elastic bodies, with specified density, elastic modulus, and Poisson’s ratio. Particle shapes may vary from single spherical particles to more complex-shaped composite particles. The particles may be introduced into the domain by random or structured injection at a point, surface, or volume, depending on the application. More details on the formulation can be found in the Simcenter STAR-CCM+ User’s Manual and OpenFOAM website.

97 MATHEMATICS AND COMPUTING↗

HYTHIRM Radiance Modeling and Image Analyses in Support of STS-119, STS-125 and STS-128 Space Shuttle Hypersonic Re-entries

We provide the first geometrically accurate (i.e., 3-D) temperature maps of the entire windward surface of the Space Shuttle during hypersonic reentry. To accomplish this task we began with estimated surface temperatures derived from CFD models at integral high Mach numbers and used them, the Shuttle's surface properties and reasonable estimates of the sensor-to-target geometry to predict the emitted spectral radiance from the surface (in units of W sr-1 m-2 nm-1). These data were converted to sensor counts using properties of the sensor (e.g. aperture, spectral band, and various efficiencies), the expected background, and the atmosphere transmission to inform the optimal settings for the near-infrared and midwave IR cameras on the Cast Glance aircraft. Once these data were collected, calibrated, edited, registered and co-added we formed both 2-D maps of the scene in the above units and 3-D maps of the bottom surface in temperature that could be compared with not only the initial inputs but also thermocouple data from the Shuttle itself. The 3-D temperature mapping process was based on the initial radiance modeling process. Here temperatures were guessed for each node in a well-resolved 3-D framework, a radiance model was produced and compared to the processed imagery, and corrections to the temperature were estimated until the iterative process converged. This process did very well in characterizing the temperature structure of the large asymmetric boundary layer transition the covered much of the starboard bottom surface of STS-119 Discovery. Both internally estimated accuracies and differences with CFD models and thermocouple measurements are at most a few percent. The technique did less well characterizing the temperature structure of the turbulent wedge behind the trip due to limitations in understanding the true sensor resolution. (Note: Those less inclined to read the entire paper are encouraged to read an Executive Summary provided at the end.)

Gibson, David M.↗