Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flow conservation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,387 records · Page 77

Isotherm Modeling and Techno-Economic Analysis of Contactor Technologies for New Tetraamine-Appended MOF for NGCC Applications

Recently, a family of tetraamine-functionalized metal-organic frameworks (MOF) has been reported as promising sorbent materials for capturing CO2 from flue gas conditions relevant to natural gas combined cycle (NGCC) applications. The main advantages of these materials are their two-step cooperative CO2 adsorption, which gives rise to unusual two step-shaped CO2 adsorption profiles and their high thermal stability. This work presents the modelling of the two-transition isotherm of the tetraamine-appended MOF, N,N'-bis(3-aminopropyl)-1,4-diaminobutane (3-4-3)-appended Mg2(dobpdc), and the techno-economic analysis (TEA) of carbon capture processes utilizing this sorbent. Due to the unusual isotherm shapes of the experimental CO2 adsorption data for tetraamine-appended Mg2(dobpdc) and the strong nonlinearity of CO2 loading with respect to temperature and pressure, we tested two different models which use logistic functions for representing the different isotherm behaviors in the different pressure ranges. The first model uses the quadratic isotherm model in the low-pressure region, the Langmuir isotherm model in the middle pressure range, and the dual site Langmuir isotherm model in the high-pressure range. To model the transition between regions we used the arctangent functions independent of temperature and the thermal effect was accounted by using the Clausius-Clapeyron relation]. The second model is an extension of the weighted dual-site Langmuir isotherm model presented by Hughes et al.. In this extended model, the dual-site Langmuir isotherm is employed in the three transition regions, using temperature-dependent logistic functions to activate or deactivate the isotherm model in the low, middle, and high-pressure ranges. Both models fit the experimental data quite well with root mean squared errors (RMSE’s) of 0.41 and 0.17 for model 1 and model 2, respectively. Since model 2 resulted in a lower RMSE, it was leveraged for the development of the gas/solid contactor models used by the TEA. Specifically, two different contactor models, an axial-flow fixed bed and moving bed contactor, were developed as part of this work. These models are dynamic, pressure-driven, and consist of mass, energy, and momentum conservation equations. A kinetic model was also developed by performing parameter estimation using experimental fixed bed breakthrough data. These models are then used to simulate CO2 capture processes from the flue gas generated from a ~600 gross MW NGCC power plant. A cost model was developed which considers the capital cost of the reactors and the significant operating costs such as steam and electricity. Using NETL’s Framework for the Optimization and Quantification of Uncertainty of Uncertainty and Surrogates tool (FOQUS), which has the capability of linking models built using numerous modelling platforms with derivative-free optimization solvers, a techno-economic optimization of the carbon capture processes was performed which minimizes the cost of capture.

Caballero, Daison↗

Transfer and Flushing Evaluation for Potential Precipitated Solids in the 241-AP-106 to EMF Transfer Pipeline for DFLAW

Waste transfer pipeline flushing is necessary to maintain pipeline performance, but added flush fluid increases the nuclear waste inventory. During the Direct Feed Low-Activity Waste (DFLAW) mission, waste transfer lines from the Interim Low-Activity Waste Storage Tank (i.e., 241-AP-106) to the Hanford Waste Treatment and Immobilization Plant (WTP) Low-Activity Waste (LAW) Facility must have the capability of being flushed to prevent accumulation of solids and to mitigate corrosion concerns. Operational experience at Hanford has shown that the solids precipitation occurs during processing operations from dilute Hanford liquids. Pacific Northwest National Laboratory was requested to evaluate the potential precipitated solid particle transport and flushing operational capabilities of the AP-106 through the Effluent Management Facility (EMF) low-point drain LAW feed pipeline. The evaluation results do not provide operational requirements, but rather provide a scoping basis for understanding the potential operational significance of solids precipitation in the pipeline. The salt and aluminum phase solids most likely to precipitate from the LAW feeds during cooling, evaporation, or mixing have particle densities that are estimated to range from 1.62 to 2.78 g/mL and a spherical particle size range of 8 to 2100 μm, depending on the solid phase considered. Application of the method to calculate the critical deposition velocity required by TFC-ENG-STD-26 at conservatively bounding estimates for the solids concentrations demonstrates that potentially precipitated solids may deposit on the bottom of the transfer pipe invert at the lower LAW feed pipeline flow rate, but the upper LAW feed pipeline flow rate will likely prevent deposition. The system pressure limit for the AP-106 through the EMF low-point drain LAW feed pipeline likely exceeds any pressure loss at a calculated critical deposition velocity even for the most adverse potential precipitated solid. However, should solids settle as a result of no-flow conditions, the flow capabilities of the AP-106 through the EMF low-point drain LAW feed pipeline are shown, based on the available literature, to be potentially inadequate for effective flushing operations to remove solids, even at solids concentrations below the maximum specified for LAW feed. Thus, stepwise solids accumulations over multiple transfers may be an issue, and eventual line plugging may occur. Further inadequacies may potentially be realized should the cohesive nature of certain precipitates be accounted for.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Gila Water Resources II: Using Earth Observations to Identify Wildfire Impacts on Hydrologic Functions and Recovery in the Gila National Forest

Wildfires have the potential to cause devastating and long-lasting impacts on ecological systems. In the Gila NF, wildfire events have occurred with increasing frequency and severity over recent years. These disturbances such as the historic Whitewater-Baldy Complex Fire (2012) and Silver Fire (2013) have raised concerns over post-fire soil erosion, flooding, debris flows, and vegetation recovery. Understanding the connection between burn events and ecological functions is crucial for developing effective land management practices within the Gila NF that ensure the conservation of the watershed. Our Gila Water Resources II team worked in partnership with the US Department of Agriculture (USDA) US Forest Service’s (USFS) Gila National Forest and Region 3. This project was designed to provide insight into the influence of wildfires on increased flooding events and to determine if restoration efforts in the Gila NF are having a beneficial and noticeable impact on recovery. The goals of this project included generating data-supported end products to inform land management decisions, including the prioritization of specific regions and time scales for post-fire restoration efforts. To understand recovery trends and hydrologic impact in the Gila NF between 2000 and 2020, this project used NASA EO and ancillary data, including, but not limited to, Landsat 5 Thematic Mapper (TM), Landsat 7 Enhanced Thematic Mapper Plus (ETM+), Landsat 8 Operational Land Imager (OLI), USGS stream gauge data, and data provided by USDA USFS’s Gila National Forest and Region 3.

Water Resources↗

Recent Advancements in the PATO Material Response Code

Introduction: Predicting the complicated multiphysics phenomena during atmospheric entry requires high-fidelity modeling tools to refine estimates of mission risks during entry. To this end, new capabilities are being added to the Porous-material Analysis Toolbox based on OpenFOAM (PATO). PATO is an open-source software for Computational Material Response (CMR) of reactive porous materials submitted to high-temperature environments. The objective of this work is to highlight current efforts to add to and improve upon the modeling capabilities of PATO. These include efforts to loosely couple PATO with other discipline specialized codes including hypersonic Computational Fluid Dynamics (CFD), to assess the interaction effects between pyrolysis gas blowing and the boundary layer, and Computational Solid Mechanics (CSM), to address modeling of mechanical erosion. Other refinements include surface phenomena modeling capabilities to address the effects of silicone-based coatings applied to the TPS during flight preparation, and a unified multiphase solver for a mixed porous-material and plain-fluid domain. Coupling CMR with CFD (CMR/CFD): A loose coupling between PATO and the Data Parallel Line Relaxation (DPLR) CFD code has been achieved by making use of a blowing boundary condition at the heatshield surface available in DPLR. Starting with heat flux estimates with no pyrolysis gas blowing at the surface, blowing gases are computed by the CMR and passed to the CFD such that aerothermal properties of the environment can be recomputed for a new CMR computation. This leads to an iterative process which is supplemented with an estimate of the radiative heat flux using the Nonequilibrium air radiation (NEQAIR) program. The entire iterative process is illustrated in Figure 1. This coupling strategy has been utilized in computing the MSL material response. The goal is to compare the coupled CMR/CFD results with material response results obtained using traditional blowing corrections. Coupling CMS with CMR: A mechanical erosion model is currently being implemented in PATO to account for the additional mass removal induced by high shear conditions. The modeling process at each timestep consists of updating the mechanical properties as a function of temperature and computing the stress tensor and displacement fields of the material. Then, a failure criteria model determines the regions in which the stress exceeds the ultimate strength values resulting in mesh movement to account for mass removal. This model allows the material response simulation to compute the recession due to both oxidation and shear-induced erosion. The model is demonstrated by computing material response of sphere-cone arc jet samples. Surface Modeling Capabilities: NuSil, a silicone-based coating, was sprayed onto the MSL and Mars 2020 heatshields to mitigate shedding of phenolic dust. To better understand the effects of the NuSil coating on the material response, a novel model has been implemented in PATO. In this model, the equilibrium of the charred NuSil surface is modeled as pure silica, and a constant offset, inspired by the classical spallation model, is added to the the char blowing rate and wall enthalpy to reproduce HyMETS experimental results. The model has also been used to estimate the 3D material response of the MSL heatshield. Unified Solver: In addition to the iterative loose coupling approach mentioned above, a multiphase unified solver is being developed to couple the environment (plain-fluid phase) and the porous-material phase. The solver is based on the volume averaged conservation of mass, momentum, and energy for the macroscale with closure models which include microscale effects through effective physicochemical properties. The unified solver has been used to compute flow through a porous plug and solve the Beavers and Joseph problem. Since the strong coupling between phases is inherent to this solver, modeling assumptions present in other coupling methods of material response are mitigated. This strategy also makes it feasible to capture the competition between surface and volume ablation in the same computational domain, which is usually not possible with other coupling approaches.

Material Response↗

International Space Station Bipropellant Plume Contamination Model Update for Short Thruster Pulse Widths

The International Space Station (ISS) Bipropellant Plume Contamination Model has been a vital tool for characterizing the thruster plume-induced contamination environment at the ISS but was not intended to be used for very short thruster pulse widths. This paper presents an updated model that ensures full start-up and shut-down phases are modeled for each thruster firing, when the majority of liquid phase contaminant mass is released. The updated ISS Bipropellant Plume Contamination Model prevents potential under-prediction of thruster plume-induced contamination due to visiting vehicle proximity operations and provides a way to take advantage of thruster start-up and shut-down data performance data gathered during thruster test programs, if available. The International Space Station (ISS) Bipropellant Plume Contamination Model developed by Boeing Space Environments is a semi-empirical model anchored in flight experiment data and has been a vital tool for characterizing the thruster plume-induced contamination environment at the ISS.[1] The current model utilizes flight experiment data from the Plume Impingement Contamination (PIC) and Shuttle Plume Impingement Flight Experiment (SPIFEX) studies, which include Orbiter 3870 N Primary Reaction Control System (PRCS) and Russian 130 N thrusters operating in pulse mode with 80-100 ms.[2] As the next generation of crew and cargo visiting vehicles are developed and arrive at ISS, minimum pulse widths of vehicle thrusters used for ISS proximity operations have decreased significantly below 80 ms. An update to the ISS Bipropellant Contamination Model is needed to prevent potential under-prediction of thruster plume-induced contamination for these very short pulse widths. Contamination due to thruster plumes occurs in the liquid phase (i.e., unburned or partially burned propellant in the plume). Liquid phase releases primarily occur during thruster start-up and shut-down phases, with the steady state phase contributing a small amount to the total contaminant mass released. The current ISS Bipropellant Plume Contamination Model includes functional dependencies on thruster parameters (thrust, mass flow rate, specific impulse (Isp)), distance to receiver surface, and angle off plume centerline. For modeling purposes, contaminant mass released scales linearly with pulse width. This is a conservative approach for pulse widths greater than 80 ms but could under-predict contamination for very short pulse widths (i.e., with little or no steady state phase). This paper presents a model update to add a functional dependency on commanded pulse width to ensure the portion of the pulse spent in start-up and shut-down phases is appropriately modeled, when the majority of liquid phase contaminant mass is released. The updated model prevents potential under-prediction of thruster plume-induced contamination due to visiting vehicle proximity operations while providing a way to take advantage of thruster start-up and shut-down data gathered during test programs, if available. Developing the plume model for a specific thruster using the ISS Bipropellant Plume Contamination Model must be done in consideration of all available thruster performance data. Example thruster data will be used to illustrate this point and options for implementing the model update with existing ISS Bipropellant Plume Contamination Model code will be discussed.

plume contamination↗

Recent Advancements in the PATO Material Response Code

Introduction: Predicting the complicated multiphysics phenomena during atmospheric entry requires high-fidelity modeling tools to refine estimates of mission risks during entry. To this end, new capabilities are being added to the Porous-material Analysis Toolbox based on OpenFOAM (PATO) [1,2,3]. PATO is an open-source software for Computational Material Response (CMR) of reactive porous materials submitted to high-temperature environments. The objective of this work is to highlight current efforts to add to and improve upon the modeling capabilities of PATO. These include efforts to loosely couple PATO with other discipline specialized codes including hypersonic Computational Fluid Dynamics (CFD), to assess the interaction effects between pyrolysis gas blowing and the boundary layer, and Computational Solid Mechanics (CSM), to address modeling of mechanical erosion. Other refinements include surface phenomena modeling capabilities to address the effects of silicone-based coatings applied to the TPS during flight preparation, and a unified multiphase solver for a mixed porous-material and plain-fluid domain. Coupling CMR with CFD (CMR/CFD): A loose coupling between PATO and the Data Parallel Line Relaxation (DPLR) [4] CFD code has been achieved by making use of a blowing boundary condition at the heatshield surface available in DPLR. Starting with heat flux estimates with no pyrolysis gas blowing at the surface, blowing gases are computed by the CMR and passed to the CFD such that aerothermal properties of the environment can be recomputed for a new CMR computation. This leads to an iterative process which is supplemented with an estimate of the radiative heat flux using the Nonequilibrium air radiation (NEQAIR) [5] program. The entire iterative process is illustrated in Figure 1. This coupling strategy has been utilized in computing the MSL material response. The goal is to compare the coupled CMR/CFD results with material response results obtained using traditional blowing corrections. Coupling CMS with CMR: A mechanical erosion model is currently being implemented in PATO to account for the additional mass removal induced by high shear conditions. The modeling process at each timestep consists of updating the mechanical properties as a function of temperature and computing the stress tensor and displacement fields of the material. Then, a failure criteria model determines the regions in which the stress exceeds the ultimate strength values resulting in mesh movement to account for mass removal. This model allows the material response simulation to compute the recession due to both oxidation and shear-induced erosion. The model is demonstrated by computing material response of sphere-cone arc jet samples. Surface Modeling Capabilities: NuSil, a silicone-based coating, was sprayed onto the MSL and Mars 2020 heatshields to mitigate shedding of phenolic dust. To better understand the effects of the NuSil coating on the material response, a novel model has been implemented in PATO. In this model, the equilibrium of the charred NuSil surface is modeled as pure silica, and a constant offset, inspired by the classical spallation model, is added to the the char blowing rate and wall enthalpy to reproduce HyMETS experimental results. The model has also been used to estimate the 3D material response of the MSL heatshield [6]. Unified Solver: In addition to the iterative loose coupling approach mentioned above, a multiphase unified solver is being developed to couple the environment (plain-fluid phase) and the porous-material phase. The solver is based on the volume averaged conservation of mass, momentum, and energy for the macroscale with closure models which include microscale effects through effective physicochemical properties. The unified solver has been used to compute flow through a porous plug and solve the Beavers and Joseph problem [7]. Since the strong coupling between phases is inherent to this solver, modeling assumptions present in other coupling methods of material response are mitigated. This strategy also makes it feasible to capture the competition between surface and volume ablation in the same computational domain, which is usually not possible with other coupling approaches.

Thermal Protection Systems↗

Experimental investigation of power transient flow boiling

To advance the mechanistic understanding of power transients in flow boiling, experimental investigations are performed to capture the thermal response of four different kinds of cladding materials, including FeCrAls, to Fuchs reactivity initiated accident transients and linear ramp power transients. The transition boiling regime is phenomenologically reflected in power-transient flow boiling, which is different from steady-state boiling. In light of this, the critical heat flux that is usually featured with temperature overshooting does not indicate thermal safety margin of cladding materials as conservatively as the maximum heat flux, which is greater than critical heat flux in the power transient boiling curve. This is physically attributed to the thermal energy deposition on the cladding wall. In addition, the thermal energy released from the cladding wall results in different boiling heat transfer coefficient for the decreasing and increasing power stages respectively during the Fuchs power transient. The maximum heat flux difference gap, which appears due to changes in cladding materials and transient time scales, is appreciable under an intermediate heat convection regime. However, this difference gap is gradually weakened by the progressive increasing of mass flux and/or inlet subcooling due to the enhanced dominance of heat convection over heat conduction. Finally, the small difference gap of maximum heat flux is found to be insignificant under the weak heat convection because of power transient induced annular flow instability.

42 ENGINEERING↗

A minimally-dissipative low-Mach number solver for complex reacting flows in Open$\mathrm{FOAM}$

Large eddy simulation (LES) has become the de-facto computational tool for modeling complex reacting flows, especially in gas turbine applications. However, readily usable general-purpose LES codes for complex geometries are typically academic or proprietary/commercial in nature. In this report the objective is to develop and disseminate an open source LES tool for low-Mach number turbulent combustion using the OpenFOAM framework. In particular, a collocated-mesh approach suited for unstructured grid formulation is provided. Unlike other fluid dynamics models, LES accuracy is intricately linked to so-called primary and secondary conservation properties of the numerical discretization schemes. This implies that although the solver only evolves equations for mass, momentum, and energy, the implied discrete equation for kinetic energy (square of velocity) should be minimally-dissipative. Here, a specific spatial and temporal discretization is imposed such that this kinetic energy dissipation is minimized. The method is demonstrated using manufactured solutions approach on regular and skewed meshes for a canonical flow and a lab-scale turbulent flow problem.

42 ENGINEERING↗

Field test and simulation evaluation of variable refrigerant flow systems performance

Our study aims to compare the performance of a Variable Refrigerant Flow (VRF) system with a Roof Top Unit, (RTU) and a variable-air-volume (VAV) system through field tests and energy simulations. The field test was conducted in as similar conditions as possible between the two systems, such as the installation and operating environments of heating, the ventilation and air conditioning (HVAC) system, including internal heat gain and outdoor conditions, including buildings to compare the performance of the two systems accurately. A VRF system and RTU were installed at the test building located in Oak Ridge, Tennessee, in the USA. The same internal heat gain was generated at the same operating time of the two systems using lighting, electric heaters, and humidifiers inside the building. The HVAC system was alternately operated between cooling and heating operations to acquire energy performance data and to compare energy usage. Furthermore, an hourly building energy simulation model was developed with regard to the VRF system and RTU, and then the model was calibrated using actual measured data. Then, annual energy consumption of the two systems were compared and analyzed using the calibrated model. Moreover, additional analysis was conducted when the controlled discharge air temperature in the RTU was changed. The field test result showed that when energy consumptions of two systems were compared at the same outdoor conditions, using the weather-normalized model, the VRF system exhibited an energy reduction of approximately 17% during cooling operation and of approximately 74% during heating operations. A comparison on the annual energy consumption using simulations showed that the VRF system reduced energy consumption more than that of the RTU by 60%. Furthermore, when the discharge air temperature in the RTU was controlled according to the outdoor air temperature, energy consumption of the RTU was reduced by 6% in cooling operations and by 18% in heating operation. As a result, energy consumption of the VRF system was reduced by more than that of the RTU by 55% approximately.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Mass transport limitations in polymer electrolyte water electrolyzers using spatially-resolved current measurement

Here this work utilizes spatially-resolved current measurements to provide insight into mass transport limitations in electrolyzers that are not observable from traditional polarization measurement. In this study, two types of flow-fields (parallel and triple-serpentine) and two types of diffusion media (patterned porous thin titanium foil LGDLs and Ir-coated titanium felt PTLs) were examined. A non-uniform current distribution dominated by mass transport limitations was observed to be instigated by the restriction of liquid water transport to catalyst sites. Additionally, conditions are revealed which yield similar polarization performance but dissimilar current distributions. In such cases, the transport limitations for different architectures and porous media affect polarization in different regions of the active area. Furthermore, the triple-serpentine flow-field used in this study performs better than the parallel flow-field under mass transport limited operating conditions. This indicates that the parallel flow-field used in this study is more susceptible to starvation than the triple-serpentine flow-field for the electrolyzer studied.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Passive Room-to-Room Air Transfer, Fresno, California (Fact Sheet)

Field testing was performed in a retrofit unoccupied test house in Fresno, California. Three air-based heating, ventilation, and air conditioning (HVAC) distribution systems - a typical airflow ducted system to the bedrooms, a low airflow ducted system to the bedrooms, and a system with no ductwork to the bedrooms - were evaluated during heating, cooling, and midseason conditions. The relative ability of each of the three systems was assessed with respect to relevant Air Conditioning Contractors of America (ACCA) and ASHRAE standards for house temperature uniformity and stability, respectively. Computational fluid dynamics (CFD) modeling also was performed and refined based on comparison to field test results to determine the air flow rate into the bedrooms of over-door and bottom-of-door air transfer grilles.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Separating Sensible and Latent Cooling with Electrically Charged Rotating Vortexes and Vapor Capturing Air Handler Technology

The dehumidification of air is a highly energy-intensive process, typically requiring vapor compression systems to cool the air below its dew point temperature. This project investigated air dehumidification achieved by electrostatic and dielectrophoresis forces in electrosprays and centrifugal principles in cyclone separators. As a practical implementation of these processes, Electrospray Vortical Flow eXchanger (EVFX) prototypes were designed, fabricated, and tested. Electrosprays were employed to energize water droplets, which then captured water vapor via dielectrophoresis. A cyclone separator was positioned downstream of the electrosprays. This combination increased the dehumidification efficiency and controlled the flow of the water vapor condensate removed from the air.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Electronically conductive polymer binder for lithium-ion battery electrode

A family of carboxylic acid groups containing fluorene/fluorenon copolymers is disclosed as binders of silicon particles in the fabrication of negative electrodes for use with lithium ion batteries. Triethyleneoxide side chains provide improved adhesion to materials such as, graphite, silicon, silicon alloy, tin, tin alloy. These binders enable the use of silicon as an electrode material as they significantly improve the cycle-ability of silicon by preventing electrode degradation over time. In particular, these polymers, which become conductive on first charge, bind to the silicon particles of the electrode, are flexible so as to better accommodate the expansion and contraction of the electrode during charge/discharge, and being conductive promote the flow battery current.

25 ENERGY STORAGE↗

Transonic wing design using potential-flow codes - Successes and failures

The state-of-the-art of transonic wing design by use of computer codes based on the potential flow-theory is presented. The capabilities and limitations of these codes are exemplified by several experiment-theory correlations, including an assessment of pressure distribution from isolated wing and wing-body codes. Computer codes using both conservative and nonconservative differencing schemes were used, and the effects of boundary-layer corrections were considered. Results showed that calculations from a full potential, isolated code correlate well with data from an isolated wing test, but may give poor predictions of the aerodynamic characteristics of some wing-body configurations. Boundary-layer correlations were found to have only moderate effects on experiment-theory correlation. Aeroelastic effects were considered important for high aspect ratio wings of low to moderate thickness, and viscous effects were minimal for typical cruise conditions, even for Reynold's numbers as low as two million. The effect of wind-tunnel walls on experiment-theory correlations remained inconclusive. A wing-body code was used to calculate the flow field about a wing-body configuration with body-mounted engines, typical of the transonic Biz-Jet aircraft.

Hicks, R. M.↗

HPC Institutional Computing Project: W15_lesreactiveflow KIVA-hpFE Development: A Robust and Accurate Engine Modeling Software

KIVA-hpFE is a high performance computer software for solving the physics of multi-species and multiphase turbulent reactive flow in complex geometries having immersed moving parts. The code is written in Fortran 90/95 and can be used on any computer platform with any popular complier. The code is in two versions, a serial version and a parallel version utilizing MPICH2 type Message Passing Interface (MPI or Intel MPI) for solving distributed domains. The parallel version is at least 30x faster than the serial version and much faster than our previous generation of parallel engine modeling software, by many factors. The 5th generation algorithm construction is a Galerkin type Finite Element Method (FEM) solving conservative momentum, species, and energy transport equations along with two-equation turbulent model k-ω Reynolds Averaged Navier-Stokes (RANS) model and a Vreman type dynamic Large Eddy Simulation (LES) method. The LES method is capable modeling transitional flow from laminar to fully turbulent; therefore, this LES method does not require special hybrid or blending to walls. The FEM projection method also uses a Petrov-Galerkin (P-G) stabilization along with pressure stabilization. We employ hierarchical basis sets, constructed on the fly with enrichment in areas associated with relatively larger error as determined by error estimation methods. In addition, when not using the hp-adaptive module, the code employs Lagrangian basis or shape functions. The shape functions are constructed for hexahedral, prismatic and tetrahedral elements. The software is designed to solve many types of reactive flow problems, from burners to internal combustion engines and turbines. In addition, the formulation allows for direct integration of solid bodies (conjugate heat transfer), as in heat transfer through housings, parts, cylinders. It can also easily be extended to stress modeling of solids, used in fluid structure interactions problems, solidification, porous media modeling and magneto hydrodynamics.

97 MATHEMATICS AND COMPUTING↗