Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ANSYS Fluent”

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 55 records · Page 3

Integrated Heat Exchanger-Phase Change Material Thermal Energy Storage System

The purpose of this study is to experimentally investigate the thermal performance of an innovative thermal energy storage (TES) system that combines the advantages of the phase-change material (PCM)/graphite foam latent heat TES medium developed at Argonne National Laboratory (Argonne) and the internally supported plate-fin (ISPE) cell architecture heat transfer fluid (HTF) flow channels developed at Brayton Energy (Brayton). Several essential tasks were accomplished: (1) Thermal property characterization. Thermal properties of the graphite foam were characterized, providing necessary data for experimental result analysis and numerical simulation. (2) Design and optimization of lab-scale test module. Based on Brayton’s full-scale heat exchanger (HX)-TES system, the experimental test module was designed, optimized, and fabricated. (3) Thermal performance testing and data analysis. Five cycle tests were successfully conducted—including one with approximately 3.5 psig of pressure applied to the diaphragms—to investigate the thermal performance of the experimental test module for charging and discharging. Temperature profiles were generated for each charging test and discharging test as a function of time. The temperature profiles clearly show three TES stages: sensible heat (temperature increase), latent heat (melting), and sensible heat (temperature increase) for the charging process. Similarly, the temperature profiles clearly show three thermal energy release stages: sensible heat (temperature decrease), latent heat (solidification), and sensible heat (temperature decrease). Melting and solidification of the PCM generally occurred in relatively narrow temperature ranges, indicated by the flattened temperature regions in the temperature profiles. These phase changes ranged approximately 3°C for melting and 3.5°C for solidification. The charging and discharging temperature profiles were similar for similar experimental parameter tests whether or not pressure was applied to the diaphragm to eliminate the gap between the HX surface and the TES subsystem. This indicates that the effect of a small gap between the HX surface and the TES subsystem is insignificant for charging and discharging. (4) Comparison of experimental data and simulation results. We compared the experimental data to the numerical simulation results. Numerical simulations were conducted by using the ANSYS FLUENT 2019 R3 commercial computational fluid dynamics software. The predicted phase-change times agreed reasonably well with those from the experimental data. In most cases, the estimated time differences between the relative phase changes were within 16%. The predicted start and end times for the charging process agreed well with those from the experimental data. However, the simulation results showed earlier start and end times than the experimental data for the discharging process. Overall, the experimental data and its comparison with the simulation predictions verified the technical viability of the integrated ISPF HX-PCM/graphite foam latent-heat TES system.

25 ENERGY STORAGE↗

CFD Analysis of RLUOB Zone 1 HEPA Filter Plenum and Testing Manifolds

The purpose of this computational fluid dynamics (CFD) analysis is to ensure that Camfil Farr’s upstream and downstream injection and sampling manifolds can meet or exceed the requirements outlined in the ASME AG-1 1997 a(2000) Code on Nuclear Air and Gas Treatment for the testing of HEPA and adsorbent filters. This paper will present a numerical simulation of airflow in the Radiological Laboratory Utility Office Building (RLUOB) zone 1 HEPA filter plenum and testing manifolds using the commercial CFD software ANSYS FLUENT 2020R1. The CFD analysis focuses on the investigation of the air flow distribution and air-aerosol mixing uniformity. The evaluation was done for all steps of the modeling process: grid generation, physics setup, simulation, and post-processing. The mass flow rate in each section of the zone 1 injection and sampling manifolds is also reported.

42 ENGINEERING↗

Perovskite Sorbent Oxygen Separation Modeling with MFiX

This document chronicles the development and implementation of computational kinetic rate models that capture absorption and desorption characteristics of the National Energy Technology Laboratory (NETL) developed perovskite, Sr 1-x Ca x FeO 3-δ . Two paths to create accurate kinetic rates were followed: (1) an isothermal rate approach where thermogravimetric (TGA) data are recast as oxygen capacities through a pseudo-second order Lagergren equation (He et al., 2009); and (2) a more traditional Arrhenius approach where experimental data are fit with a power law model to derive associate activation energies (Bulfin et al., 2020a). For reference, the mathematics and associate experimental strategies that support these derivations are included in this report. In addition, computational fluid dynamics (CFD) models were developed to utilize both kinetic rate formulations and applied to simulate oxygen uptake and release in small scale scenarios. The program Multiphase Flow with interphase eXchanges (MFiX) was used to create: (1) discrete element method (DEM) simulations of a single tube of granular perovskite experiencing isothermal O 2 -absorption; and desorption and (2) two-fluid-model (TFM) non-isothermal simulations of perovskite O 2 -absorption and desorption tubes that share a wall. Conjugate heat transfer between steel walled tubes and the perovskite bed were managed with user-defined functions. As the project moves to simulating larger scale devices that will require more robust conjugate heat transfer methods, developed kinetic rates and associate methodologies have been recast for use in the ANSYS Fluent CFD program.

36 MATERIALS SCIENCE↗

A numerical evaluation of the ambient air temperature in the Electron-Ion Collider tunnel

The Electron-Ion Collider (EIC) is a next-generation collider-accelerator that may require consistent operating temperature conditions for the beams within the accelerator tunnels to maintain stable operation. Variations in ambient temperature within the tunnel can cause thermal expansion of beampipe and component supports and can negatively affect the tunnel equipment, impacting the stability of the beamline. Modifications will be made to the Relativistic Heavy Ion Collider (RHIC) at Brookhaven National Laboratory (BNL) to create the EIC, which necessitates a temperature model that addresses these modifications. To approach this problem, the consistency of temperature changes in different tunnel sections was first evaluated by plotting RHIC tunnel temperature data at various times of the day and year. From this data, a tunnel section was selected and a 2D temperature model was created for RHIC, EIC, and EIC with added cooling configurations. Soil temperature data was analyzed to determine the maximum, average, and mode soil temperatures, which were used as boundary conditions in different temperature scenarios. Computational fluid dynamics modeling was used to create 2D temperature profiles for the configurations. From this model, the predicted temperatures indicate that further analysis is required to validate the boundary conditions and benchmark the current conditions to allow the prediction of the tunnel ambient conditions at EIC. This research can be used as a preliminary model to create an EIC tunnel cooling system that will increase the operational stability of the EIC. As a result of my work this summer, I have become familiar with computational fluid dynamics, including creating fluid dynamic simulations using ANSYS Fluent and related software. I have also learned about the project process required for planning large-scale engineering projects.

43 PARTICLE ACCELERATORS↗

Numerical investigation of the AP1000 response following loss-of-coolant accident using PCTRAN and CFD to support R&D of SMRs

Small Modular Reactors (SMRs) present a promising solution for the future of sustainable energy, offering advantages such as reduced waste generation, advanced passive safety features, and potential cost efficiencies. Ongoing research focuses on the design and development of SMRs, addressing challenges through numerical simulations and experimental test data. Given that many next-generation reactors are first-of-a-kind (FOK), current operational reactors serve as valuable benchmarks for understanding various thermal-hydraulic phenomena during postulated design basis accidents (DBAs), including loss-of-coolant accidents (LOCA), main steam line breaks (MSLB), and steam generator tube ruptures (SGTR). This study employs the AP1000 model Personal Computer Transient Analyzer (PCTRAN) to simulate reactor responses to small-break LOCA scenarios, specifically analyzing break sizes of 2 and 10 inches. Additionally, ANSYS FLUENT software is utilized to assess containment responses to large-break LOCAs, focusing on the quantification of decay heat removal via natural convection. Verification of the PCTRAN results is achieved using Westinghouse data, with CFD results corroborating the findings. The outcomes demonstrate strong agreement with Westinghouse data, confirming the accuracy of the simulations.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Assessment and Improvement Proposals for High-Pressure Rinsing Processes during Cleanroom Preparation of Superconducting Radio-Frequency Cavities for Particle Accelerators

The goal of this internship is to develop a fluid dynamics characterization of the high-pressure flow used for cleaning the internal surface of the SRF cavities. Initially, the High Pressure Rinsing (HPR) process is introduced, focusing on the problems of the SRF cavities and how this method can be used to increase their acceleration efficiency and performance. The characterization of the problem is then presented, including an introduction to the multiphase flow modeling and the analysis of the Fermilab facility for the HPR. To accomplish the goal, a preliminary estimation of the variables was performed using MATLAB. Then, multiple 2D simulations were developed with Ansys Fluent, and the influence of the main parameter was investigated. The simulation results are then presented and analysed, with a focus on the possible future developments of the HPR process fluid dynamic characterization.

43 PARTICLE ACCELERATORS↗

Effects of turbulence-radiation coupling in pressurized oxy-combustion

Concerns over climate change have led to numerous efforts in developing low-carbon energy technologies. Pressurized oxy-combustion (POC) is a promising candidate to reduce carbon emission in power generation. Due to the strong impact of pressure on thermal radiation, heat transfer in POC differs significantly from the situation of conventional atmospheric pressure combustion. Thus thermal radiation in POC needs to be investigated to aid new combustor development. The present computational work is a step in this direction, initiating a systematic analysis of thermal radiation and heat transfer in a pilot-scale POC combustor, which has been developed at Washington University in St. Louis (WUSTL). In a POC process, pulverized coal is burned under elevated pressure and O2-CO2 environment. While most previous computational works on flame radiation focused on atmospheric pressure condition, this work considered a pressurized flame. Specifically, a 15-bar POC combustor of power 50 kW is modeled employing the Ansys FLUENT commercial platform, using both Reynolds-averaged Navier-Stokes (RANS) modeling and large-eddy simulation (LES). A recently published and validated global radiation model is used to predict the radiative property of the flue gas. The discrete ordinates (DO) radiation model is chosen to solve the radiative transport equation. Incident radiation on the walls of the combustor is identified and investigated. It is revealed that for this pilot-scale, pressurized combustor, thermal radiation exhibits similar patterns in the RANS and LES models. However, the magnitudes of radiation are different in both models, presumably because interactions between radiation and turbulence are embedded in the LES model but not in the RANS model. This difference thereby underlines the major goal of this research: investigating the impact of such turbulence-to-radiation coupling on thermal radiation in a pilot-scale POC setting. Most previous studies on turbulence-to-radiation interactions focused on bench-scale, atmospheric pressure flames. This work is the first effort to extend this study to a pilot-scale, pressurized flame.

Li, Lei↗

Positron Injector for Ce+BAF

A baseline concept for a continuous wave (CW) polarized positron injector was developed for the CEBAF at Jefferson Lab. This concept is based on the generation of CW longitudinally polarized positrons by a high-current, polarized electron beam 1 mA, 130-370 MeV, and 90% longitudinal polarization that passes through a rotating, water-cooled, tungsten target. The simulations for the Ce+BAF injector at the Low Energy Recirculator Facility (LERF) are performed, including positron beam generation, capture, energy selection, and acceleration to 123 MeV. The positron yield (or positron current) and longitudinal polarization are calculated considering the longitudinal and transverse CEBAF acceptances (<1% energy spread, 1.2 mm bunch length and normalized emittance of 100 mm·mrad). The impact of target thickness, electron beam energy, and transverse size on positron yield within the required emittance limit is evaluated. Target thermal and structural FEM analyses (ANSYS Fluent) are performed to determine the maximum electron beam current and minimum transverse beam size on the target.

Benesch, Jay [Thomas Jefferson National Accelerato↗

Seawater Acidification and Bubble Plume Dispersion from Accidental Subsea CO 2 Pipeline Rupture: A Multiphase CFD Study

If a CO 2 reservoir or transmission pipeline were to leak, both the surrounding ecology and maritime traffic safety could be put at risk. To better understand and prepare for this risk, multiphase Computational Fluid Dynamics (CFD) models were built in ANSYS Fluent to capture the behavior of a leak once it enters the water. A 3D Eulerian–Eulerian model was used for validation, while a simplified 2D model was applied to simulate conditions at a 50-m depth. The models integrate bubble dynamics, gas holdup, CO 2 dissolution, dissolved species transport, and seawater acidification into a unified CFD framework. Mass transfer was calculated using the Hughmark correlation, and local seawater temperature and salinity were factored in to determine dissociation behavior and the relevant Henry’s Law constant. To confirm the 3D model’s accuracy, results were checked against two experimental datasets: the QICS field study and the Hauser Tank experiments. The team also modeled a hypothetical release scenario at the High Island 10L site and compared the results with earlier published work. The results show that at a depth of 50 m, the surrounding water column can completely absorb a CO 2 release at a rate of 35 kg/s, since the gas dissolves into the seawater as it rises toward the surface. Beyond confirming this mitigation capacity, the simulations shed light on how a leak would actually unfold in the environment, including the shape and movement of the rising bubble plume, how much CO 2 dissolves along the way, and the resulting shifts in seawater pH and pCO 2 . Together, this provides a practical framework for assessing how CO 2 leaks could affect marine environments in the Gulf of Mexico.

54 ENVIRONMENTAL SCIENCES↗

A Computational Study for the Utilization of Jet Pulsations in Gas Turbine Film Cooling and Flow Control

This report is the second part of a three-part final report of research performed under an NRA cooperative Agreement contract. The first part is NASA/CR-2012-217415. The third part is NASA/CR-2012-217417. Jets have been utilized in various turbomachinery applications in order to improve gas turbines performance. Jet pulsation is a promising technique because of the reduction in the amount of air removed from compressor. In this work two areas of pulsed jets applications were computationally investigated using the commercial code Fluent (ANSYS, Inc.); the first one is film cooling of High Pressure Turbine (HPT) blades and second one is flow separation control over Low Pressure Turbine (LPT) airfoil using Vortex Generator Jets (VGJ). Using pulsed jets for film cooling purposes can help to improve the effectiveness and thus allow higher turbine inlet temperature. Effects of the film hole geometry, blowing ratio and density ratio of the jet, pulsation frequency and duty cycle of blowing on the film cooling effectiveness were investigated. As for the low-pressure turbine (LPT) stages, the boundary layer separation on the suction side of airfoils can occur due to strong adverse pressure gradients. The problem is exacerbated as airfoil loading is increased. Active flow control could provide a means for minimizing separation under conditions where it is most severe (low Reynolds number), without causing additional losses under other conditions (high Reynolds number). The effects of the jet geometry, blowing ratio, density ratio, pulsation frequency and duty cycle on the size of the separated region were examined in this work. The results from Reynolds Averaged Navier-Stokes and Large Eddy Simulation computational approaches were compared with the experimental data.

Kartuzova, Olga V.↗

Flow Control Under Low-Pressure Turbine Conditions Using Pulsed Jets

This publication is the final report of research performed under an NRA/Cooperative Interagency Agreement, and includes a supplemental CD-ROM with detailed data. It is complemented by NASA/CR-2012-217416 and NASA/CR-2012-217417 which include a Ph.D. Dissertation and an M.S. thesis respectively, performed under this contract. In this study the effects of unsteady wakes and flow control using vortex generator jets (VGJs) were studied experimentally and computationally on the flow over the L1A low pressure turbine (LPT) airfoil. The experimental facility was a six passage linear cascade in a low speed wind tunnel at the U.S. Naval Academy. In parallel, computational work using the commercial code FLUENT (ANSYS, Inc.) was performed at Cleveland State University, using Unsteady Reynolds Averaged Navier Stokes (URANS) and Large Eddy Simulations (LES) methods. In the first phase of the work, the baseline flow was documented under steady inflow conditions without flow control. URANS calculations were done using a variety of turbulence models. In the second phase of the work, flow control was added using steady and pulsed vortex generator jets. The VGJs successfully suppressed separation and reduced aerodynamic losses. Pulsed operation was more effective and mass flow requirements are very low. Numerical simulations of the VGJs cases showed that URANS failed to capture the effect of the jets. LES results were generally better. In the third phase, effects of unsteady wakes were studied. Computations with URANS and LES captured the wake effect and generally predicted separation and reattachment to match the experiments. Quantitatively the results were mixed. In the final phase of the study, wakes and VGJs were combined and synchronized using various timing schemes. The timing of the jets with respect to the wakes had some effect, but in general once the disturbance frequency was high enough to control separation, the timing was not very important.

Volino, Ralph J.↗

Flow Control Under Low-Pressure Turbine Conditions Using Pulsed Jets: Experimental Data Archive

This publication is the final report of research performed under an NRA/Cooperative Interagency Agreement, and includes a supplemental CD-ROM with detailed data. It is complemented by NASA/CR-2012-217416 and NASA/CR-2012-217417 which include a Ph.D. Dissertation and an M.S. thesis respectively, performed under this contract. In this study the effects of unsteady wakes and flow control using vortex generator jets (VGJs) were studied experimentally and computationally on the flow over the L1A low pressure turbine (LPT) airfoil. The experimental facility was a six passage linear cascade in a low speed wind tunnel at the U.S. Naval Academy. In parallel, computational work using the commercial code FLUENT (ANSYS, Inc.) was performed at Cleveland State University, using Unsteady Reynolds Averaged Navier Stokes (URANS) and Large Eddy Simulations (LES) methods. In the first phase of the work, the baseline flow was documented under steady inflow conditions without flow control. URANS calculations were done using a variety of turbulence models. In the second phase of the work, flow control was added using steady and pulsed vortex generator jets. The VGJs successfully suppressed separation and reduced aerodynamic losses. Pulsed operation was more effective and mass flow requirements are very low. Numerical simulations of the VGJs cases showed that URANS failed to capture the effect of the jets. LES results were generally better. In the third phase, effects of unsteady wakes were studied. Computations with URANS and LES captured the wake effect and generally predicted separation and reattachment to match the experiments. Quantitatively the results were mixed. In the final phase of the study, wakes and VGJs were combined and synchronized using various timing schemes. The timing of the jets with respect to the wakes had some effect, but in general once the disturbance frequency was high enough to control separation, the timing was not very important. This is the supplemental CD-ROM

Volino, Ralph J.↗

Computational Investigation of Novel Tip Leakage Mitigation Methods for High Pressure Turbine Blades

This paper presents preliminary findings on a possible approach to reducing tip leakage losses. In this paper a computational study was conducted on the Energy Efficient Engine (EEE) High Pressure Turbine (HPT) rotor tip geometry using the commercial numerical solver ANSYS FLUENT. The flow solver was validated against aerodynamic data acquired in the NASA Transonic Turbine Blade Cascade facility. The scope of the ongoing study is to computationally investigate how the tip leakage and overall blade losses are affected by (1) injection from the tip near the pressure side, (2) injection from the tip surface at the camber line, and (3) injection from the tip surface into the tip separation bubble. The objective is to identify the locations on the tip surface at which to place appropriately configured blowing keeping in mind the film cooling application of tip blowing holes. The validation was conducted at Reynolds numbers of 85,000, 343,000, and 685,000 and at engine realistic flow conditions. The coolant injection simulations were conducted at a Reynolds number of 343,000 based on blade chord and inlet velocity and utilized the SST turbulence model in FLUENT. The key parameters examined are the number of jets, jet angle and jet location. A coolant to inlet pressure ratio of 1.0 was studied for angles of +30 deg, -30 deg, and 90 deg to the local free stream on the tip. For the 3 hole configuration, 3 holes spaced 3 hole diameters apart with length to diameter ratio of 1.5 were used. A simulation including 11 holes along the entire mean camber line is also presented (30 deg toward suction side). In addition, the effect of a single hole is also compared to a flat tip with no injection. The results provide insight into tip flow control methods and can be used to guide further investigation into tip flow control. As noted in past research it is concluded that reducing leakage flow is not necessarily synonymous with reducing losses due to leakage.

Ibrahim, Mounir↗

Computational Investigation of Novel Tip Leakage Mitigation Methods for High Pressure Turbine Blades

This paper presents preliminary findings on a possible approach to reducing tip leakage losses. In this paper a computational study was conducted on the EEE (Energy Efficient Engine) HPT (High Pressure Turbine) rotor tip geometry using the commercial numerical solver ANSYS FLUENT. The flow solver was validated against aerodynamic data acquired in the NASA Transonic Turbine Blade Cascade facility. The scope of the ongoing study is to computationally investigate how the tip leakage and overall blade losses are affected by 1. injection from the tip near the pressure side, 2. injection from the tip surface at the camber line, and 3. injection from the tip surface into the tip separation bubble. The objective is to identify the locations on the tip surface at which to place appropriately configured blowing keeping in mind the film cooling application of tip blowing holes. The validation was conducted at Reynolds numbers of 85,000, 343,000 and 685,000 and at engine realistic flow conditions. The coolant injection simulations were conducted at a Reynolds number of 343,000 based on blade chord and inlet velocity and utilized the SST turbulence model in FLUENT. The key parameters examined are the number of jets, jet angle and jet location. A coolant to inlet pressure ratio of 1.0 was studied for angles of +30 deg., -30 deg. and 90 deg. to the local free stream on the tip. For the 3 hole configuration, 3 holes spaced 3 hole diameters apart with length to diameter ratio of 1.5 were used. A simulation including 11 holes along the entire mean camber line is also presented (30 degrees toward suction side). In addition, the effect of a single hole is also compared to a flat tip with no injection. The results provide insight into tip flow control methods and can be used to guide further investigation into tip flow control. As noted in past research it is concluded that reducing leakage flow is not necessarily synonymous with reducing losses due to leakage.

Ibrahim, Mounir↗

An Internal Thermal Environment Model of an Aluminized Solid Rocket Motor with Experimental Validation

Due to the severity of the internal solid rocket motor (SRM) environment, very few direct measurements of that environment exist; therefore, the appearance of such data provides a unique opportunity to assess current thermal/fluid modeling capabilities. As part of a previous study of SRM internal insulation performance, the internal thermal environment of a laboratory-scale SRM featuring aluminized propellant was characterized with two types of custom heat-flux calorimeters: one that measured the total heat flux to a graphite slab within the SRM chamber and another that measured the thermal radiation flux. Therefore, in the current study, a thermal/fluid model of this lab-scale SRM was constructed using ANSYS Fluent to predict not only the flow field structure within the SRM and the convective heat transfer to the interior walls, but also the resulting dispersion of alumina droplets and the radiative heat transfer to the interior walls. The dispersion of alumina droplets within the SRM chamber was determined by employing the Lagrangian discrete phase model that was fully coupled to the Eulerian gas-phase flow. The P1-approximation was engaged to model the radiative heat transfer through the SRM chamber where the radiative contributions of the gas phase were ignored and the aggregate radiative properties of the alumina dispersion were computed from the radiative properties of its individual constituent droplets, which were sourced from literature. The convective and radiative heat fluxes computed from the thermal/fluid model were then compared with those measured in the lab-scale SRM test firings and the modeling approach evaluated.

Martin, Heath T.↗

Modeling of Non-Isothermal Cryogenic Fluid Sloshing

A computational fluid dynamic model was used to simulate the thermal destratification in an upright self-pressurized cryostat approximately half-filled with liquid nitrogen and subjected to forced sinusoidal lateral shaking. A full three-dimensional computational grid was used to model the tank dynamics, fluid flow and thermodynamics using the ANSYS Fluent code. A non-inertial grid was used which required the addition of momentum and energy source terms to account for the inertial forces, energy transfer and wall reaction forces produced by the shaken tank. The kinetics-based Schrage mass transfer model provided the interfacial mass transfer due to evaporation and condensation at the sloshing interface. The dynamic behavior of the sloshing interface, its amplitude and transition to different wave modes, provided insight into the fluid process at the interface. The tank pressure evolution and temperature profiles compared relatively well with the shaken cryostat experimental test data provided by the Centre National D'Etudes Spatiales.

Cryogenics↗

Investigation of Gas Seeding for Planar Laser-Induced Fluorescence in Hypersonic Boundary Layers

Numerical simulations of the gas-seeding strategies required for planar laser-induced fluorescence in a Mach 10 (approximately Mach 8.2 postshock) airflow were performed. The work was performed to understand and quantify the adverse effects associated with gas seeding and to assess various types of seed gas that could potentially be used in future experiments. In prior experiments, NO and NO2 were injected through a slot near the leading edge of a flatplate wedge model used in NASA Langley Research Center's 31 in. Mach 10 air tunnel facility. In this paper, nitric oxide, krypton, and iodine gases were simulated at various injection rates. Simulations showing the deflection of the velocity boundary layer for each of the cases are presented. Streamwise distributions of velocity and concentration boundary-layer thicknesses, as well as vertical distributions of velocity, temperature, and mass distributions, are presented for each of the cases. A comparison between simulated streamwise velocity profiles and experimentally obtained molecular tagging velocimetry profiles using a nitric oxide seeding strategy is performed to verify the influence of such a strategy on the boundary layer. The relative merits of the different seeding strategies are discussed. The results from a custom solver based on OpenFOAM version 2.2.1 are compared against results obtained from ANSYS® Fluent version 6.3.

Arisman, C. J.↗