Search NASA⌕ Search

SEARCH · Search NASA

Results for “Three dimensional computational fluid dynamics model”

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 37 records · Page 2

Forest structural complexity and ignition pattern influence simulated prescribed fire effects

Background: Forest structural characteristics, the burning environment, and the choice of ignition pattern each influence prescribed fire behaviors and resulting fire effects; however, few studies examine the influences and interactions of these factors. Understanding how interactions among these drivers can influence prescribed fire behavior and effects is crucial for executing prescribed fires that can safely and effectively meet management objectives. To analyze the interactions between the fuels complex and ignition patterns, we used FIRETEC, a three-dimensional computational fluid dynamics fire behavior model, to simulate fire behavior and effects across a range of horizontal and vertical forest structural complexities. For each forest structure, we then simulated three different prescribed fires each with a unique ignition pattern: strip-head, dot, and alternating dot. Results: Forest structural complexity and ignition pattern affected the proportions of simulated crown scorch, consumption, and damage for prescribed fires in a dry, fire-prone ecosystem. Prescribed fires in forests with complex canopy structures resulted in increased crown consumption, scorch, and damage compared to less spatially complex forests. The choice of using a strip-head ignition pattern over either a dot or alternating-dot pattern increased the degree of crown foliage scorched and damaged, though did not affect the proportion of crown consumed. We found no evidence of an interaction between forest structural complexity and ignition pattern on canopy fuel consumption, scorch, or damage. Conclusions: We found that forest structure and ignition pattern, two powerful drivers of fire behavior that forest managers can readily account for or even manipulate, can be leveraged to influence fire behavior and the resultant fire effects of prescribed fire. These simulation findings have critical implications for how managers can plan and perform forest thinning and prescribed burn treatments to meet risk management or ecological objectives.

54 ENVIRONMENTAL SCIENCES↗

A Study of Propane Combustion in a Spark-Ignited Cooperative Fuel Research (CFR) Engine

Liquefied petroleum gas (LPG), whose primary composition is propane, is a promising candidate for heavy-duty vehicle applications as a diesel fuel alternative due to its CO 2 reduction potential and high knock resistance. To realize diesel-like efficiencies, spark-ignited LPG engines are proposed to operate near knock-limit over a wide range of operating conditions, which necessitates an investigation of fuel-engine interactions that leads to end-gas autoignition with propane combustion. This work presents both experimental and numerical studies of stoichiometric propane combustion in a sparkignited (SI) cooperative fuel research (CFR) engine. Engine experiments are initially conducted at different compression ratio (CR) values, and the effects of CR on engine combustion are characterized. A three-pressure analysis (TPA) model based on the two-zone combustion concept is developed in GT-Power and validated using test results to estimate in-cylinder wall temperatures, residual gas fraction, etc. This model is further utilized to examine end-gas chemistry by enabling the SI turbulent flame combustion and unburned gas chemical kinetics modules. Finally, a three-dimensional (3D) computational fluid dynamic (CFD) model of the CFR engine is developed in CONVERGE, where the G-equation and SAGE detailed chemical kinetics models are implemented for combustion modeling. Here, a 153 species reduced chemical kinetics mechanism derived from the detailed NUIGMech1.1 mechanism based on the ignition delay and laminar flame speed (LFS) studies is used to generate an LFS lookup table and to describe end-gas autoignition chemistry. Multi-cycle Reynolds-averaged Navier-Stokes (RANS) simulations are then performed for the tested CRs, and the numerical model is shown to be capable of predicting the propane combustion characteristics, particularly the end-gas autoignition behavior.

20 FOSSIL-FUELED POWER PLANTS↗

ROM-Based Surrogate Systems Modeling of EBR-II

We report the System Analysis Module (SAM), developed and maintained by Argonne National Laboratory, is designed to provide whole-plant transient safety analysis capabilities for a number of advanced non-light water reactors, including sodium-cooled fast reactor (SFR), lead-cooled fast reactor (LFR), and molten salt reactor (MSR)/fluoride-salt-cooled high-temperature reactor (FHR) designs. SAM is primarily constructed as a systems-level analysis tool, with the potential to incorporate reduced order models from three-dimensional computational fluid dynamics (CFD) simulations to improve characterization of complex, multidimensional physics. It is recognized that the computational expense associated with CFD can be intractable for various engineering analyses, such as uncertainty quantification, inference, and design optimization. This paper explores the reducibility of a SAM model using recent advances in randomized linear algebra techniques, which attempt to find recurring patterns in the various realizations generated by a model after randomly perturbing all its input parameters. The reduction is described in terms of fewer degrees of freedom (DOFs), referred to as the active DOFs, for the model variables such as input model parameters and model responses. The results indicate that there is significant room for additional reduction that may be leveraged for additional computational gains when employing SAM for engineering-intensive analyses that require repeated model executions. Different from physics-based reduction approaches, the proposed approach allows one to estimate upper bounds on the reduction errors, which are rigorously developed in this work. Finally, different methods for surrogate model construction, such as regression and neural network-based training, are employed to correlate the input and output active DOFs, which are related back to the original variables using matrix-based linear transformations.

42 ENGINEERING↗

Experimental evaluation of an advanced Space Shuttle main engine hot-gas manifold design concept

This study, using an extensively modified, full-scale space shuttle main engine (SSME) hot-gas manifold (HGM), established a detailed aerodynamic data base to support development of an advanced, three-dimensional, fluid-dynamic analysis computer model. In addition, the advanced SSME hot-gas manifold design used in this study demonstrated improved flow environment (uniformity) in the fuel side turbine exit and transfer duct exit regions. Major modifications were incorporated in the full-scale HGM flow test article model using two large transfer ducts on the fuel turbine side of the HGM in place of the three small transfer ducts in the present design. Other model features included an increases in the flow areas downstream of the 180-degree turn and in the fishbowl regions.

Pelaccio, D. G.↗

Reliability of Lumped Thermodynamic Hydrogen Fueling Model under Slow-Fill Conditions

This study verifies the reliability of lumped thermodynamic fuel cell electric vehicle (FCEV) tank model under considerably slow-fill conditions. Many countries put efforts into research and development of FCEVs. As part of the efforts, thermodynamic hydrogen fueling models have been developed to understand the hydrogen temperature in the onboard tanks of FCEVs during the fueling process. Most of the models treat the hydrogen inside the tanks to be a lumped system and assume the hydrogen temperature to be uniform throughout the tanks. In other words, the hydrogen temperature is treated as average bulk temperature. This study certifies whether treating the hydrogen temperature in a tank as the bulk gas temperature is suitable by evaluating the temperature distributions inside the tank by a three-dimensional computational fluid dynamics (3D CFD) model.

DIRECT ENERGY CONVERSION,HYDROGEN↗

A Numerical Study on the Energy Performance of a Novel Furnace With Acidic Gas Trap Absorbers

Natural gas furnaces are widely used in US residential and commercial building markets. An important issue for natural gas furnaces is serious corrosion and fouling problems caused by acidic gas, such as SOx. An advanced adsorption technology based on acidic gas trap (AGT) absorbers offers the possibility to remove SOx acidic gas from natural gas furnaces with high efficiency and low cost, thereby enabling the development of condensing furnaces without the use of expensive corrosion resistant materials in the heat exchanger. A three-dimensional (3D) computational fluid dynamics (CFD) model has been developed to evaluate the heat transfer performance of a furnace with AGT absorbers and to compare it with a baseline conventional furnace without the AGT. Moreover, an axisymmetric model has been built focusing on the absorbing process in the AGT. The baseline conventional furnace used for the study is a commercial condensing furnace (Rheem 92% AFUE 84,000 BTU Multi-Position Gas Furnace). This furnace was completely disassembled, and the dimensions of each part were carefully measured and used to build a detailed CFD model. A model representing the new furnace, incorporating the AGT absorbers, was developed by adding the AGT system to the conventional furnace model. For the CFD analysis, a mixture model was employed to characterize the heat and mass transfer during the condensing process in the furnace while considering three components—air, water vapor and liquid water. Condensation takes place in the condensing heat exchanger, where water vapor changes phase to liquid water, and the latent heat is thus used in the furnace for useful heating. The simulation results characterize the energy performance of both the conventional furnace and the novel furnace with AGT absorbers, as well as the reactive processing in the AGT. These results provide insightful guidance for the development of the AGT absorber-based furnace from the perspective of its energy performance and will be used to further optimize this novel furnace design.

Laclair, Tim↗

Application of the kappa-omega Turbulence Model to Quasi-Three-Dimensional Turbomachinery Flows

Many computational fluid dynamics codes for turbomachinery use the Baldwin-Lomax (B-L) turbulence model. It is easy to implement in two dimensions and works well for predicting overall turbomachinery performance. However, it is awkward to implement in three dimensions, often has difficulty finding the length scale, has a crude transition model, and neglects freestream turbulence, surface roughness, and mass injection. The kappa-omega model developed by Wilcox is an appealing alternative for several reasons. First, it is the only two-equation model that can be integrated to the wall without requiring damping functions or the distance to the wall, and hence, should behave well numerically. Second, the effects of freestream turbulence, surface roughness, and mass injection are easily included in the model. Finally, transition can be simulated using the low Reynolds number version of the model. Menter applied the kappa-model to external flows and showed very good results for flows with adverse pressure gradients. Liu and Zheng described their implementation of the kappa-model in a cascade code that included an area change term to account for endwall convergence. They validated the model for a flat plate, and compared the B-L and kappa-models to measured surface pressures for a low-pressure turbine cascade. Since they did not use the low Reynolds number version of the model, their results showed problems resulting from early transition. In this Note the low Reynolds number kappa-model was incorporated in the author's quasi-three-dimensional turbomachinery analysis code. The code includes the effects of rotation, radius change, and stream-surface thickness variation, and also includes the B-L turbulence model. The kappa-omega model was implemented using many of Menter's recommendations and an implicit approximate-factorization scheme described by Baldwin and Barth. The model was tested for a transonic compressor with rotation and variable stream-surface radius and height, and for a transonic turbine vane with transition and heat transfer. Results were compared to the B-L model and to experimental data.

Chima, Rodrick V.↗

Numerical Analysis of Fuel Effects on Advanced Compression Ignition Using a Cooperative Fuel Research Engine Computational Fluid Dynamics Model

Growing environmental concerns and demand for a better fuel economy are driving forces that motivate the research for more advanced engines. Multi-mode combustion strategies have gained attention for their potential to provide high thermal efficiency and low emissions for light-duty applications. These strategies target optimizing the engine performance by correlating different combustion modes to load operating conditions. The extension from boosted spark ignition (SI) mode at high loads to advanced compression ignition (ACI) mode at low loads can be achieved by increasing the compression ratio and utilizing intake air heating. Further, in order to enable an accurate control of intake charge condition for ACI mode and rapid mode-switches, it is essential to gain fundamental insights into the autoignition process. Within the scope of ACI, homogeneous charge compression ignition (HCCI) mode is of significant interest. It is known for its potential benefits, operation at low fuel consumption, low NOx, and particulate matter (PM) emissions. In the present work, a virtual Cooperative Fuel Research (CFR) engine model is used to analyze fuel effects on ACI combustion. In particular, the effect of fuel octane sensitivity (S) (at constant Research Octane Number (RON)) on autoignition propensity is assessed under beyond-RON (BRON) and beyond-MON (BMON) ACI conditions. The three-dimensional CFR engine computational fluid dynamics (CFD) model employs a finite-rate chemistry approach with a multi-zone binning strategy to capture autoignition. Two binary blends with Research Octane Number (RON) of 90 are chosen for this study: primary reference fuel (PRF) with S = 0 and toluene–heptane (TH) blend with S = 10.8, representing paraffinic and aromatic gasoline surrogates. Reduced mechanisms for these blends are generated from a detailed gasoline surrogate kinetic mechanism. Simulation results with the reduced mechanisms are validated against experimental data from an in-house CFR engine, with respect to in-cylinder pressure, heat release rate, and combustion phasing. Thereafter, the sensitivity of combustion behavior to ACI operating condition (BRON versus BMON), air-fuel ratio (λ = 2 and 3), and engine speed (600 and 900 rpm) is analyzed for both fuels. It is shown that the sensitivity of a fuel’s autoignition characteristics to λ and engine speed significantly differs at BRON and BMON conditions. Moreover, this sensitivity is found to vary among fuels, despite the same RON. It is also observed that the presence of low-temperature heat release (LTHR) under BRON condition leads to more sequential autoignition and longer combustion duration than BMON condition. Finally, the study indicates that the octane index (OI) fails to capture the trend in the variation of autoignition propensity with S under the BMON condition.

33 ADVANCED PROPULSION SYSTEMS↗

3-D CFD Simulation and Validation of Oxygen-Rich Hydrocarbon Combustion in a Gas-Centered Swirl Coaxial Injector using a Flamelet-Based Approach

Injector design is a critical part of the development of a rocket Thrust Chamber Assembly (TCA). Proper detailed injector design can maximize propulsion efficiency while minimizing the potential for failures in the combustion chamber. Traditional design and analysis methods for hydrocarbon-fuel injector elements are based heavily on empirical data and models developed from heritage hardware tests. Using this limited set of data produces challenges when trying to design a new propulsion system where the operating conditions may greatly differ from heritage applications. Time-accurate, Three-Dimensional (3-D) Computational Fluid Dynamics (CFD) modeling of combusting flows inside of injectors has long been a goal of the fluid analysis group at Marshall Space Flight Center (MSFC) and the larger CFD modeling community. CFD simulation can provide insight into the design and function of an injector that cannot be obtained easily through testing or empirical comparisons to existing hardware. However, the traditional finite-rate chemistry modeling approach utilized to simulate combusting flows for complex fuels, such as Rocket Propellant-2 (RP-2), is prohibitively expensive and time consuming even with a large amount of computational resources. MSFC has been working, in partnership with Streamline Numerics, Inc., to develop a computationally efficient, flamelet-based approach for modeling complex combusting flow applications. In this work, a flamelet modeling approach is used to simulate time-accurate, 3-D, combusting flow inside a single Gas Centered Swirl Coaxial (GCSC) injector using the flow solver, Loci-STREAM. CFD simulations were performed for several different injector geometries. Results of the CFD analysis helped guide the design of the injector from an initial concept to a tested prototype. The results of the CFD analysis are compared to data gathered from several hot-fire, single element injector tests performed in the Air Force Research Lab EC-1 test facility located at Edwards Air Force Base.

Richardson, Brian↗

Evaluating Chemical Kinetics Predictions for Propane Using 3-D and 0-D Models in a Boosted Spark-Ignited Engine

Propane has been shown to be a promising alternative fuel to reduce emissions while simultaneously achieving high efficiencies in medium- and heavy-duty engines. These high-power density applications require boosted engines which, combined with high compression ratio, can lead to auto-ignition and knock. While three-dimensional (3-D) computational fluid dynamics (CFD) models are often used for resolving the complex fluid flow in engines, these models can become computationally expensive when simulating detailed chemical kinetics. Likewise, zero-dimensional (0-D) models are computationally concise enough for kinetics development, but lack any flow-field information which governs the flame propagation processes in spark ignition (SI) engines. This work presents a comprehensive comparison between 3-D and 0-D closed cycle simulations at knocking conditions in a high compression ratio high stroke-to-bore ratio propane engine. In order to initialize the flow-field for the 3-D closed cycle (intake valve closing, (IVC) to exhaust valve opening, (EVO)) simulation, a motored multi-cycle 3-D model was run using Converge to create a map at IVC, reducing the computational time. The map allowed a non-homogeneous 3-D closed cycle simulation to be satisfactorily validated against experiments, while a homogeneous case using only the turbulence field mapping was also simulated, mimicking 0-D modeling. The 3-D simulations were used to prescribe the initial conditions (e.g., IVC thermodynamics, speciation, burn-rate profile) for a 2-zone 0-D SI engine model in Chemkin Pro for both cases. It was found that 2-zone 0-D modeling underpredicted the knock onset timing, likely due to the lack of thermal stratification in the unburned gas region. Future work will carry multi-zone 0-D modeling to capture the fuel auto-ignition in the unburned region.

Douvry-Rabjeau, Julien [Oakland University, Roches↗

Dissertation Defense Computational Fluid Dynamics Uncertainty Analysis for Payload Fairing Spacecraft Environmental Control Systems

Spacecraft thermal protection systems are at risk of being damaged due to airflow produced from Environmental Control Systems. There are inherent uncertainties and errors associated with using Computational Fluid Dynamics to predict the airflow field around a spacecraft from the Environmental Control System. This paper describes an approach to quantify the uncertainty in using Computational Fluid Dynamics to predict airflow speeds around an encapsulated spacecraft without the use of test data. Quantifying the uncertainty in analytical predictions is imperative to the success of any simulation-based product. The method could provide an alternative to traditional "validation by test only" mentality. This method could be extended to other disciplines and has potential to provide uncertainty for any numerical simulation, thus lowering the cost of performing these verifications while increasing the confidence in those predictions. Spacecraft requirements can include a maximum airflow speed to protect delicate instruments during ground processing. Computational Fluid Dynamics can be used to verify these requirements; however, the model must be validated by test data. This research includes the following three objectives and methods. Objective one is develop, model, and perform a Computational Fluid Dynamics analysis of three (3) generic, non-proprietary, environmental control systems and spacecraft configurations. Several commercially available and open source solvers have the capability to model the turbulent, highly three-dimensional, incompressible flow regime. The proposed method uses FLUENT, STARCCM+, and OPENFOAM. Objective two is to perform an uncertainty analysis of the Computational Fluid Dynamics model using the methodology found in "Comprehensive Approach to Verification and Validation of Computational Fluid Dynamics Simulations". This method requires three separate grids and solutions, which quantify the error bars around Computational Fluid Dynamics predictions. The method accounts for all uncertainty terms from both numerical and input variables. Objective three is to compile a table of uncertainty parameters that could be used to estimate the error in a Computational Fluid Dynamics model of the Environmental Control System /spacecraft system. Previous studies have looked at the uncertainty in a Computational Fluid Dynamics model for a single output variable at a single point, for example the re-attachment length of a backward facing step. For the flow regime being analyzed (turbulent, three-dimensional, incompressible), the error at a single point can propagate into the solution both via flow physics and numerical methods. Calculating the uncertainty in using Computational Fluid Dynamics to accurately predict airflow speeds around encapsulated spacecraft in is imperative to the success of future missions.

CFD↗

Dissertation Defense: Computational Fluid Dynamics Uncertainty Analysis for Payload Fairing Spacecraft Environmental Control Systems

Spacecraft thermal protection systems are at risk of being damaged due to airflow produced from Environmental Control Systems. There are inherent uncertainties and errors associated with using Computational Fluid Dynamics to predict the airflow field around a spacecraft from the Environmental Control System. This paper describes an approach to quantify the uncertainty in using Computational Fluid Dynamics to predict airflow speeds around an encapsulated spacecraft without the use of test data. Quantifying the uncertainty in analytical predictions is imperative to the success of any simulation-based product. The method could provide an alternative to traditional validation by test only mentality. This method could be extended to other disciplines and has potential to provide uncertainty for any numerical simulation, thus lowering the cost of performing these verifications while increasing the confidence in those predictions.Spacecraft requirements can include a maximum airflow speed to protect delicate instruments during ground processing. Computational Fluid Dynamics can be used to verify these requirements; however, the model must be validated by test data. This research includes the following three objectives and methods. Objective one is develop, model, and perform a Computational Fluid Dynamics analysis of three (3) generic, non-proprietary, environmental control systems and spacecraft configurations. Several commercially available and open source solvers have the capability to model the turbulent, highly three-dimensional, incompressible flow regime. The proposed method uses FLUENT, STARCCM+, and OPENFOAM. Objective two is to perform an uncertainty analysis of the Computational Fluid Dynamics model using the methodology found in Comprehensive Approach to Verification and Validation of Computational Fluid Dynamics Simulations. This method requires three separate grids and solutions, which quantify the error bars around Computational Fluid Dynamics predictions. The method accounts for all uncertainty terms from both numerical and input variables. Objective three is to compile a table of uncertainty parameters that could be used to estimate the error in a Computational Fluid Dynamics model of the Environmental Control System spacecraft system.Previous studies have looked at the uncertainty in a Computational Fluid Dynamics model for a single output variable at a single point, for example the re-attachment length of a backward facing step. For the flow regime being analyzed (turbulent, three-dimensional, incompressible), the error at a single point can propagate into the solution both via flow physics and numerical methods. Calculating the uncertainty in using Computational Fluid Dynamics to accurately predict airflow speeds around encapsulated spacecraft in is imperative to the success of future missions.

Fluid Dynamics↗

A CFD Analysis of Easterly Wind Flow Impacting the Vehicle Assembly Building

In an attempt to explain the high loss of panels from the south face of the Vehicle Assembly Building (VAB) during Hurricane Frances, a three-dimensional computational fluid dynamics (3-D CFD) model was developed to simulate local velocity and pressure distributions resulting from such a storm. A preconditioned compressible Navier-Stokes flow solver 1 was used to compute the flow field around the VAB complex, including the Launch Control Center, the Low and High Bays of the VAB, and several outbuildings in the immediate LC-39 area. The mapping of the forces and velocities on and along the affected faces of the VAB correlated surprisingly well with the extensive damage areas realized on both on the south face and on the southeast section of the roof. The model results were also consistent with the minimal damage seen on the east, north, and west faces of the structure.

Vu, B. T.↗

Analysis of thermo-chemical nonequilibrium models for carbon dioxide flows

The aerothermodynamics of thermochemical nonequilibrium carbon dioxide flows is studied. The chemical kinetics models of McKenzie and Park are implemented in separate three-dimensional computational fluid dynamics codes. The codes incorporate a five-species gas model characterized by a translational-rotational and a vibrational temperature. Solutions are obtained for flow over finite length elliptical and circular cylinders. The computed flowfields are then employed to calculate Mach-Zehnder interferograms for comparison with experimental data. The accuracy of the chemical kinetics models is determined through this comparison. Also, the methodology of the three-dimensional thermochemical nonequilibrium code is verified by the reproduction of the experiments.

Rock, Stacey G.↗

Surface Flow Visualization of the High-Lift Common Research Model

A 10% scale version of the High-Lift Common Research Model (CRM-HL) was tested in the NASA Langley 14- by 22-Foot Subsonic Tunnel (14x22) in support of the NASA Advanced Air Transport Technology (AATT) Project. The CRM-HL experiment included various configurations such as conventional and simple-hinged flaps, with and without engine nacelle/pylon, with and without nacelle chine, different Active Flow Control (AFC) methods (sweeping jets, alternating pulsed jets, and preconditioned boundary layer blowing), and their various parameters. This particular study is focused on the surface flow visualization of the conventional CRM-HL model at landing configuration. The conventional CRM-HL model with the single-slotted Fowler flap system serves as a baseline for the AFC-enabled simplified high-lift configuration as well as a high-lift technology development platform due to its publicly open geometry. Surface flow visualizations were performed using fluorescent minitufts, which were found to be nonintrusive to the aerodynamic performance. Tuft flow visualizations are supplemented with the relevant pressure and force measurements in order to understand the flow characteristics developed on the conventional CRM- HL model. In addition, three dimensional, unsteady, compressible Computational Fluid Dynamic (CFD) simulations were performed for selective cases. The surface streamlines and transverse velocity fluctuations obtained by the CFD simulations are qualitatively compared to the tuft direction and tuft unsteadiness, respectively. Force measurements of the CRM-HL model show performance degradation at higher angles of attack. Surface flow visualizations revealed the performance loss due to the nacelle/pylon wake that grows with angle of attack and eventually promotes flow separation over the inboard wing. This performance loss was successfully recovered by placing a chine on the engine nacelle.

Koklu, Mehti↗

Development of a Virtual CFR Engine Model for Knocking Combustion Analysis

Knock is a major bottleneck to achieving higher thermal efficiency in spark ignition (SI) engines. The overall tendency to knock is highly dependent on fuel anti-knock quality as well as engine operating conditions. It is, therefore, critical to gain a better understanding of fuel-engine interactions in order to develop robust knock mitigation strategies. In the present work, a numerical model based on three-dimensional (3-D) computational fluid dynamics (CFD) was developed to capture knock in a Cooperative Fuel Research (CFR) engine. For combustion modeling, a hybrid approach incorporating the G-equation model to track turbulent flame propagation, and a homogeneous reactor multi-zone model to predict end-gas auto-ignition ahead of the flame front and post-flame oxidation in the burned zone, was employed. In addition, a hybrid methodology was implemented wherein a laminar flame speed lookup table generated a priori from a chemical kinetic mechanism could be used to provide flame speed as an input to the G-equation model, instead of using conventional empirical correlations. Multi-cycle Reynolds-averaged Navier-Stokes (RANS) simulations were performed for two different spark timings (STs) corresponding to non-knocking and knocking conditions, with other operating conditions kept the same as those of a standard research octane number (RON) test. Isooctane was considered as the fuel for the numerical study. Two different reduced kinetic mechanisms were employed to describe end-gas auto-ignition chemistry and to generate the flame speed lookup table. Experimental data, including intake/exhaust boundary conditions, was provided by an isooctane ST sweep study conducted in an in-house CFR engine. Moreover, cylinder wall/valve/port surface temperatures and residual gas fraction (RGF) were estimated using a well-calibrated one-dimensional (1-D) model. On the other hand, a novel methodology was also developed to analyze experimental data for the knocking case and identify the most representative cycle. For the non-knocking case, a good agreement was found between experiment and CFD simulation, with respect to cycle-averaged values of 10% burn point (CA10), 50% burn point (CA50), and peak pressure magnitude/location. The virtual CFR engine model was also demonstrated to be capable of predicting average knock characteristics for the knocking case, such as knock point, knock intensity, and energy of resonance, with good accuracy.

42 ENGINEERING↗

Development of Transient Coupled Chemical-Thermal-Fluid Multiphysics Simulation for Unsealed, Vented Aluminum-clad Spent Nuclear Fuel Storage Canister

As the second step towards modeling three-dimensional (3D) multi-physics computational fluid dynamics (CFD) of the unsealed containers in the IFSF facility, the previous methodology used for developing the sealed DOE standard canister model was utilized here. This 3D CFD model considers unsealed canister with small mass exchange with surrounding air as time evolves in the simulation and contains a larger network of chemical reactions to consider due to several of the side reactions occurring in the presence of water vapor, nitrogen, oxygen and carbon dioxide rather than a helium backfilled canister. Similar to the sealed system, the set of full radiolytic chemical reactions of bulk moist air were established, and then reduced to a more manageable set of equations for the computational efficiency of the 3D CFD simulation. Due to the similar size of the unsealed canisters to the sealed DOE canisters, it is assumed that the chemical and thermal response times to changes in the conditions remains similar such that the same time-stepping algorithm can be utilized. This report documents the theory of the model, workflow to establish radiolysis reaction network, and initial simulations of the evolutions of thermal fields, hydrogen gas and nitric acid concentrations within unsealed canisters in conditions relevant to INL INTEC CPP-603 over a long period of time.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Time-Accurate Computational Fluid Dynamics Analysis of a Supersonic Turbine to Assess Structural Mode Excitation

Gas turbines can experience significant unsteady loading during operation which can cause damage such as blade cracking or bearing wear. While computational fluid dynamics (CFD) modeling of full three-dimensional gas turbine flow to verify meanline design performance is standard practice within the industry, the use of time-accurate CFD makes it possible to assess the unsteady loading which can be especially problematic at operating conditions where unsteady fluid loading aligns with structural modes. The Fluid Dynamics branch (ER42) at NASA Marshall Space Flight Center has developed a time-accurate CFD analysis methodology using the unstructured, density-based, Loci-CHEM solver which is capable of simulating the relative motion required for turbomachinery applications through sliding interface meshes. The highly-parallel unstructured solver allows for full three-dimensional, time-accurate simulation of turbines including upstream and downstream components such as manifolds, nozzles, and stators, any of which could excite a turbine structural mode. Recently, Loci-CHEM was used to analyze the FASTRAC supersonic turbine to assess the loading on the rotor. Overall performance and time-averaged loading on the rotor was computed, including the axial force over the entire rotor disk which is required for force balance calculations on the turbopump. The FASTRAC turbine has supersonic nozzles upstream of the rotor that result in significant unsteady loading on the rotor blades at multiples of the nozzle pass frequency. The unsteady blade pressure loading was decomposed into the frequency domain for direct use in structural forced response analysis. Utilizing this type of CFD analysis coupled with structural analysis has proven to be a valuable tool in identifying structural modes that are being excited by unsteady fluid loading.

M.V. Fischels↗