Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computational fluid dynamics simulation”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Isolated Building Wake Experiment At Texas Tech University’s Wind Engineering Research Field Lab

Texas Tech University’s (TTU’s) Wind Engineering Research Field Lab (WERFL) building is an idealized rectangular building on the scale of a typical suburban house. It is surrounded almost entirely by unobstructed planes and is built on a turn table. WERFL is ideal for providing measurements of an idealized building under real-world conditions including atmospheric stability and variable relative wind angles. The intent of this data set is for verification and validation of computational fluid dynamics model simulations of flow around buildings. The primary focus of this dataset features 24 3D sonic anemometers (20 downwind and 4 upwind) and 10 2D sonic anemometers (9 downwind and 1 upwind) located in the immediate vicinity of the WERFL building. It also includes a variety of instruments intended to measure the ambient meteorological parameters and the undisturbed wind profile.

17 WIND ENERGY↗
DOE OSTI2022

Completion of initial reduced-order model for flammable refrigerant dispersal in residential spaces

Environmental regulations aimed at reducing global warming impacts of HVAC&R refrigerants have resulted in the phase-out of chlorofluorocarbons (CFCs) in 2010 and hydrochlorofluorocarbons (HCFCs) by 2030 in developed countries with additional restrictions on the use of hydrofluorocarbons (HFCs) set to take effect in 2036. Many of the remaining alternative refrigerants that have lower global warming potential (GWP) and that are suitable for use in HVAC&R systems (e.g., propane, difluoromethane) are flammable to some degree. Flammable refrigerants introduce new challenges and hazards to property and personal safety related to potential deflagration during system maintenance or due to leakage of the refrigerant accumulating in the conditioned space. Standards have been developed to set maximum charge limits for flammable refrigerants in HVAC&R systems; however, existing safety standards and codes still restrict their use. The bodies that maintain and update these codes need publicly available, science-based information to enable credible guidelines for setting safe charge limits for different flammable refrigerants in different HVAC&R applications. In 2016, the Alliance for Responsible Atmospheric Policy, the Air-Conditioning Heating and Refrigeration Institute (AHRI), ASHRAE, the U.S. Department of Energy (DOE), and the State of California began efforts to develop such information. As part of this effort, Oak Ridge National Laboratory (ORNL) began the current, ongoing project to examine imposed charge limits for flammable refrigerants and identify reasonable adjustments to these limits when found appropriate. Past tasks under this project have included development of experimentally vetted, computational fluid dynamics (CFD) simulation approaches to study the results of leakage of flammable refrigerants from various HVAC&R systems into different types of commercial and residential spaces. In this report, we discuss recent efforts to develop a predictive model of the flammable volume fraction and accumulated refrigerant mass in a single-room residential space resulting from the leak of a flammable refrigerant from a small room air conditioning (RAC) unit. The eventual goal is development of a model suitable for public release which could simulate a range of scenarios. In discussions with the AHRTI Flammable Refrigerant Subcommittee (FRS) at the beginning of this effort, a total of 9 input parameters were chosen for consideration including room area, room/door opening area, ventilation fan flow rate, unit/leak height, leak area, leak rate, total refrigerant charge, refrigerant molecular weight, and state of the unit fan.

42 ENGINEERING↗
DOE OSTI2023

Nanosecond Repetitively Pulsed (NRP) Plasmas: Relationship Between Induced Flow and Plasma Characteristics at Atmospheric Pressure (Final Technical Report)

The local flow fields induced by plasma discharges at atmospheric pressure can play a critical role in aerodynamic flow and combustion control and have significant implications for plasma-based technologies in medicine and environmental engineering. However, the coupling between the plasma characteristics, the geometry and the induced flow field is not well understood. The primary goal of this work is to determine the relationship between the plasma characteristics and the induced flow field for nanosecond pulsed discharges at atmospheric pressure. A comprehensive experimental investigation of plasmas produced by nanosecond high-voltage pulses in a pin-to-pin electrode configuration was conducted, employing several measurement techniques including microwave and laser Rayleigh scattering, optical emission spectroscopy, and high-speed imaging. The flow field induced by the nanosecond plasma was measured and characterized using measurements of the velocity (particle image velocimetry) and density (background oriented schlieren) with unprecedented spatiotemporal resolution. The experimental measurements were complemented by development of a one-dimensional model of the plasma discharge and a high-fidelity computational fluid dynamics (CFD) simulation of the subsequent plasma-induced flow field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Numerical Analysis and Flow Induced Vibration Studies in Pressurized Helium Gas Cooling Channels

Los Alamos National Laboratory (LANL) is working with NorthStar Medical Radioisotopes (NMR) on their efforts to produce Mo-99 from the irradiation of Mo-100 targets. The NMR target consists of an Inconel window that allows the electron beam to penetrate and irradiate a stack of Mo-100 discs. The irradiation process generates large amounts of heat and is cooled using pressurized helium gas flowing at 400 g/s. LANL provides both numerical analysis support and experimental support. Computational fluid dynamics (CFD) simulations for both 300 g/s and 400 g/s are performed on the NMR target to analyze the helium flow through the NMR system. The pressure drop, the mass flow rate through each of the disc and window channels and the force applied on the first disc are some of the quantities estimated using CFD. These results will be combined with Monte Carlo N-Particle (MCNP) estimates of heat flux into the system from the electron beam to perform conjugate heat transfer analysis (CHT). The experiments presented in this report focus on the effect of the flow induced vibrations from helium flowing through the coolant gaps between the discs. There is concern that the flow induced vibrations may lead to excessive mass loss of the enriched Mo-100 discs, and the goal would be to mitigate the mass loss. This work presents benchtop tests of flow through a subset of the target stack which involves high speed imaging and displacement measurements used to estimate disc vibration, as well as mass loss.

42 ENGINEERING↗
DOE OSTI2024

Experimental and Theoretical Evaluation of Feed Flow Collar Design for Shell Fed Hollow Fiber Membrane Modules

An experimental and theoretical study of module collar design is presented here. Hollow fiber membranes are prepared by dip coating a poly(vinylidene) (PVDF) support with a polydimethylsiloxane (PDMS) gutter layer and a Pebax 2533 selective layer. Fiber bundles with a well-defined fiber packing are prepared using a 3D printed module. A parallel fiber bundle consisting of 4-9 uniformly spaced fibers is created with printed tabs that align the fibers and create a tubesheet. The tabs are sealed within a printed case that possesses a series of external ports for gas introduction and removal. Uniquely, both port location and the use of a collar to assist fluid distribution in the shell can be varied for the same fiber bundle. Experimental measurements are compared to computational fluid dynamics (CFD) simulations. The experimental module design allows high-fidelity representation of the fiber bundle and module case in the simulations. Comparisons between experiment and simulation are in good agreement over a broad range of experimental conditions. The detrimental effect of having ports located too close, leading to stagnation regions, is captured as well as the beneficial effects of using a collar for shell-side fluid distribution around the fiber bundle. Such results help validate the use of CFD to develop high-performance module designs.

Tran, Thien↗
DOE OSTI2024

Performance Testing of a Moving-Bed Gasifier Using Coal, Biomass, and Waste Plastic Blends with Washed and Unwashed Legacy Coals and Other Waste Fuels to Generate White Hydrogen

The objective of this effort, primarily funded by the United States Department of Energy (DOE), and led by the Electric Power Research Institute, Inc. (EPRI), with support by Hamilton Maurer International (HMI) and Sotacarbo S.p.A. (Sotacarbo), has been to qualify coal, biomass, and plastic waste blends based on performance testing of selected fuel pellet compositions in a pilot-scale updraft moving-bed (UDMB) gasifier. The testing provided relevant data to advance the commercial-scale design of the moving-bed gasifier to be able to successfully use these feedstocks to produce hydrogen. In particular, the effects of waste plastics on feedstock development (i.e., blending and pelletizing) and the resulting products (i.e., syngas compositions, organic condensate production, and ash characteristics) are the focus. The gasifier used for testing is HMI’s moving-bed gasifier, which has been proven capable of gasifying nearly all coal ranks. It has also shown the ability in prior testing work to gasify wood chips (biomass). However, mixtures of these fuels with plastic wastes have not been prepared and gasified together. The three feedstocks were densified and pelletized by California Pellet Mill (CPM) to meet the feedstock size required by Sotacarbo’s 30mm ID UDMB gasifier, under contract to HMI. The technical tasks and results from this two-year research project included: (1) Feed Procurement and Preparation: Nine different tri-fuel pellets were prepared from varying compositions of fresh mined PRB coal, corn stover biomass, and car fluff waste plastics. Tri-fuel pellets were produced by CPM and shipped to Sotacarbo’s test facility in Carbonia, Sardinia, Italy. (2) Test Plan Development: A test plan was created to define the test runs to be performed. The test plan detailed the different UDMB gasification tests to be performed in Sotacarbo’s 12-inch ID pilot scale gasifier, the process monitoring instrumentation used, and the extractive samples recovered for analysis of the total gasification process mass and energy balance. (3) Gasifier Testing: Nine different gasification runs were performed in the pilot-scale gasifier at Sotacarbo using nine different fuel feedstock compositions generated from varying mixtures of PRB coal, biomass, and plastic wastes. The testing generated performance data on gasification reaction efficiency and performance, yielding relevant data for models used to scale up the gasifier design. This task also included work to refurbish and reassemble the pilot gasifier at Sotacarbo and perform a baseline 100% PRB coal run. (4) Data Analysis and Reporting: Review of the data, determination of figures of merit, and interpretation of the results are reported in the project’s final report, published in March 2024. The results show that all tri-fuel pellets gasified well and maintained structural integrity throughout the gasification process. The syngas generated can be shifted to hydrogen by using commercial syngas shifting technologies. (5) High Fidelity computational fluid dynamics (CFD) Simulation: The National Energy Technology Laboratory (NETL) team performed CFD simulations of the UDMB gasifier for two of the tri-fuel pellets gasified in Sotacarbo’s pilot scale gasifier. The kinetic mechanisms for the pyrolysis of each constituent, PRB coal, corn stover biomass, and waste plastics are based on thermogravimetric analysis performed by Sotacarbo. The gasification model was validated by comparing the predicted syngas composition at the exit of the gasifier with the measured syngas composition. In addition, the reactor’s measured internal temperature profile agreed well with the predicted internal reactor temperature profile. These results validate that the model can be used to predict the performance of the updraft moving bed gasifier for different feedstocks and operating conditions. This paper summarizes the results of the completed work in which the pelletizing procedure was validated to ensure the viability of the tri-fuel pellets for the gasification runs performed at Sotacarbo’s 30 mm UDMB gasifier. The gasification performance data from this series of nine runs will enable modeling of a full-scale HMI industrial scale gasifier supporting both combined heat and power, and Hydrogen production from coal (both fresh mined and legacy) combined with various biomass and waste plastics. Additionally, plans and progress on a follow-up project, being executed by the same project team, will be presented. In this project, a total of twenty (20) different feedstocks are being prepared from varying compositions of biomass (both woody biomass and corn stover) with a mixture of legacy coal waste, plastic waste, and refuse-derived fuel (RDF). The testing will provide information on gasification reaction efficiency/performance, yielding relevant data for models used to scale up the gasifier design to 50 megawatt electric (MWe) (equivalent hydrogen production). Tests will also be performed on a bench-scale fluidized-bed gasifier for comparison purposes. The results of this testing will be used to specify the range of feedstock blends that can be successfully gasified as well as quantify gasifier outputs based on specific blends.

08 HYDROGEN↗
DOE OSTI2021

Investigating the Effects of Chemical Mechanism on Soot Formation Under High-Pressure Fuel Pyrolysis

We performed Computational Fluid Dynamics (CFD) simulations using a Reynolds-Averaged Navier-Stokes (RANS) turbulence model of high-pressure spray pyrolysis with a detailed chemical kinetic mechanism encompassing pyrolysis of n-dodecane and formation of polycyclic aromatic hydrocarbons. We compare the results using the detailed mechanism and those found using several different reduced chemical mechanisms to experiments carried out in an optically accessible, high-pressure, constant-volume combustion chamber. Three different soot models implemented in the CONVERGE CFD software are used: an empirical soot model, a method of moments, and a discrete sectional method. There is a large variation in the prediction of the soot between different combinations of chemical mechanisms and soot model. Furthermore, the amount of soot produced from all models is substantially less than experimental measurements. All of this indicates that there is still substantial work that needs to be done to arrive at simulations that can be relied on to accurately predict soot formation.

02 PETROLEUM↗
DOE OSTI2025

Outflow Boundary Conditions for Turbine-Integrated Rotating Detonation Combustors

This study examines outflow boundary conditions (BCs) in computational fluid dynamics (CFD) simulations of a transition duct with and without guide vanes that converts supersonic flow exiting a rotating detonation combustor (RDC) to subsonic flow to drive a turbine. Since the flow exiting the transition duct has swirling shock waves with significant spatial and temporal variations in pressure, temperature, and Mach number, imposing proper BCs poses a challenge. To ensure all swirling shock waves exit the transition duct without creating non-physical reflected waves at its outlet, this study examined three outflow BCs: (1) the average pressure imposed at the duct’s outlet, (2) a nonreflecting BC (NRBC) with a specified average pressure imposed at the duct’s outlet, (3) the average pressure imposed at the outlet of an extension duct made up of a buffer layer and a sponge layer. This study is based on the three-dimensional, unsteady density-weighted-ensemble-averaged continuity, Navier–Stokes, and energy equations for a thermally perfect gas closed by the realizable k–ε model and “enhanced” wall functions. The results obtained show that imposing an average pressure at the transition duct’s outlet produces spurious waves that degrade the physical meaningfulness of the solution. When the NRBC was applied, swirling shock waves exited the duct’s outlet without creating spurious waves. However, its usage requires the gas to be thermally, as well as calorically, perfect, which this study shows could be a concern. By imposing the average pressure at the outlet of an extension duct, the gas does not need to be calorically perfect. The results obtained show the effects of the sponge layer’s length and coarsening ratio on damping nonuniformities in non-physical reflected waves to ensure the flow exiting the transition duct’s outlet can do so as if there are no boundaries present and has the desired average pressure—even though the BC is applied at the extension duct’s outlet.

gas turbines↗
DOE OSTI2020

Optimization of Window Positions for Wind-Driven Natural Ventilation Performance

This paper optimizes opening positions on building facades to maximize the natural ventilation’s potential for ventilation and cooling purposes. The paper demonstrates how to apply computational fluid dynamics (CFD) simulation results to architectural design processes, and how the CFD-driven decisions impact ventilation and cooling: (1) background: A CFD helps predict the natural ventilation’s potential, the integration of CFD results into design decision-making has not been actively practiced; (2) methods: Pressure data on building facades were obtained from CFD simulations and mapped into the 3D modeling environment, which were then used to identify optimal positions of two openings of a zone. The effect of the selected opening positions was validated with building energy simulations; (3) results: The cross-comparison study of different window positions based on different geographical locations quantified the impact on natural ventilation effectiveness; and (4) conclusions: The optimized window position was shown to be effective, and some optimal solutions contradicted the typical cross-ventilation strategy.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

Measurements and Prediction of Ash Deposition in a Cyclone-Fired Boiler Operating under Variable Load Conditions

Measurements of ash deposition rates were made between the secondary superheater and reheater sections of a 450 MW cyclone-fired lignite boiler as the operational load varied from 33 to 100%. Significant reductions in deposition rates with a decrease in operational load were observed. To uncover the causative mechanisms behind these observations, operational data from the power plant were used to carry out computational fluid dynamic (CFD) simulations of the boiler. After ascertaining that the gas temperatures and velocities at various sections within the boiler were being represented adequately, decoupled simulations of the ash deposition process on the deposit probe were carried out using a finely resolved boundary layer mesh. Fly ash particle size distribution (PSD) and its concentration for the decoupled calculations were determined from stand-alone cyclone barrel simulations. The ash partitioning (mass %) between the fly ash and slag was found to be ~50:50, which was in line with previous field observations, and it did not vary significantly across different cyclone loads. The predicted PSD of the deposit ash was concentrated in the size range 10–30 microns, which was in agreement with cross-sectional images of the deposit obtained from the measurements. At lower loads, sharp variations in the deposition rates were predicted in the gas temperature range 950–1150 K. The particle kinetic energy—particle viscosity-based capture methodology utilized in this study in conjunction with appropriate ash compositions, ash viscosity models and gas temperature estimates can help estimate slagging propensities at different loads reasonably well in these systems.

Krishnamoorthy, Gautham (ORCID:0000000255205092)↗
DOE OSTI2025

Performance Evaluation of Low-Cost Condensers Immersed in the Water Tank of a Heat Pump Water Heater

This paper reports the effect of low-cost immersed condensers on the performance of a heat pump water heater (HPWH) that does not require a water-circulating pump. The proposed design is based on the immersion of an L-style condenser coil through an opening in the top of the water tank. The novel design aims to eliminate the need for the water-circulating pumps, thereby substantially improving efficiency and reducing costs and maintenance of the HPWH system. Comprehensive computational fluid dynamics (CFD) simulations and experimental tests were performed, and the results showed that the CFD simulations were close to the experimental data. Further results confirm that the L-style condenser can introduce buoyancy-driven flow and eliminate temperature stratification in the studied HPWH water tank, thereby improving the heat transfer coefficient and coefficient of performance of the HPWH. The testing results further showed that the L-style condenser enabled more than 53% higher coefficient of performance and 400% higher heat transfer coefficient compared with the straight vertical condenser, and achieved much faster water heating. Furthermore, the effects of immersed condensers with different circuits on the HPWH performance were evaluated and compared. All results indicate that the current energy-saving and inexpensive HPWH technology is technically valid, and the benefits of performance improvement and low cost are attractive in future HPWH design.

Gao, Zhiming [ORNL] (ORCID:0000000271397995)↗
DOE OSTI2026

Investigation of Flux Spreading in a Light-Trapping, Planar-Cavity Receiver for Enclosed Solar Particle Heating

Concentrating solar thermal power (CSP) technology development has recently focused on increasing the operating temperatures to accommodate high efficiency power cycles and thermochemical processes. Inert solid particles as heat transfer media enable solar receivers to operate above 700 degrees Celsius resulting in increased system thermal efficiency compared to the conventional molten salt based CSP system. An open-cavity falling-particle solar receiver that can efficiently heat particles by direct heating from concentrated solar radiation faces challenges with large particle losses from wind and unable to support thermochemical reactions. A light-trapping, planar cavity reiver (LTPCR) where particles are indirectly heated can significantly minimize the particle losses during the operation, support thermochemical reactions, and offer scalability potential. The LTPCR features an array of vertical planar receiver/absorber panels arranged within a cavity configuration. Concentrated solar radiation from heliostats is focused onto the receiver walls, where heat is indirectly transferred to solid particles flowing inside the receiver channels. Heat transfer occurs through direct contact between the receiver panel walls and particles, and can be enhanced by fluidizing particles with air. This fluidization increases particle-wall contact and extends particle residence time, maximizing heat transfer efficiency. The unique vertical planar receiver structure originated from a near-blackbody tubular light absorber, effectively distributing the incoming solar beam spread across the panel walls and trapping light. This flux spreading effect, driven by cosine projection, converts high incident solar flux into a lower, more uniform heat flux on the panel walls. This redistribution enhances heat transfer efficiency between particle-wall or reaction gases-wall, while preventing localized overheating of the receiver panel. Indirect planar cavity solar receivers completely separate solid particles from the ambient environment that can greatly reduce the thermal losses in heated particles resulting in high efficiency at high temperatures above 700 degrees Celsius. This design ensures no particle losses to the environment during the operation while open-cavity designs can experience significant particle losses from wind. An experimental investigation was conducted to observe flux spreading on the receiver panel wall. A lab-scale prototype planar receiver, fabricated using Haynes 230 alloy, was tested under direct concentrated solar radiation using the high-flux solar furnace (HFSF) facility at NREL. The experiment was performed under normal peak radiative heat fluxes ranging from 800 to 1900 kW/m2. A temperature distribution on the panel wall was measured using a thermal imaging camera (FLIR A 6600). To prevent overheating at the receiver front tip, prism-shaped heat shields (Zircar UNIFROM C1) were placed in front of the receiver, and their influence on flux spreading was also studied. Absorbed flux distribution on the panel wall was modeled using SolTrace. The total solar power and flux distributions delivered from HFSF were determined based on the heliostat mirror optical properties, direct normal irradiance (DNI) on the on-sun testing days, peak flux measurement during the on-sun testing, and shutter/attenuator settings Due to the large incident angles of the solar beam on the panel wall, the angular optical properties of Haynes 230 alloy and Zircar heat shields were incorporated into the model. This flux distribution model was then integrated into a computational fluid dynamics (CFD) simulation to predict the receiver panel wall temperature, which was compared with the experimental measurements. Both prediction and measurements identified a temperature hotspot at the backside of the panel, indicating that the incident solar beam can fully reach to the rear of the receiver. The heat shields positioned at the front of the receiver effectively reduced the excessive temperature rise at the receiver front tip. Overall, the temperature was well distributed over the panel wall, with a minor hotspot at the back of the receiver. The model slightly overpredicted the temperature, possibly due to discrepancies in optical properties of the panel and an underprediction of thermal loss in the receiver. The advancement of the particle LTPCR offers a viable alternative to open-cavity receivers by addressing particle loss issues. Additionally, it presents a pathway for enabling solar thermochemical processes, extending CSP technology beyond power generation to fuel and chemical production.

14 SOLAR ENERGY↗
NASA NTRS1990

Supersonic flameholding by attached oblique shock waves

The object of this investigation was to explore various regimes of shock-induced combustion that may occur as a result of supersonic flow of a reactive gas mixture past a simple compression ramp or wedge. Attention was confined to combinations of wedge angles, fuel-air mixture ratios and freestream velocities that permit the resulting oblique shock or oblique detonation wave to remain attached. Computational fluid dynamic (CFD) simulations, both inviscid and viscous, were performed and analyzed in order to investigate two modes of compression ramp flameholding; namely, complete oblique detonation wave (ODW) and shock-induced combustion. In addition, the possibility of flame persistence or hysteresis was explored. It was preliminarily concluded that inclusion of viscous effects may be essential to accurately predict the occurrence of the various flameholding modes.

Fort, James A.↗
NASA NTRS1991

Integrated flow and structural modeling for rocket engine component test facility propellant systems

A set of PC-based computational Dynamic Fluid Flow Simulation models is presented for modeling facility gas and cryogenic systems. Data obtained provide important information regarding performance envelope parameters for the facility using different engine components; time-dependent valve setting for controlling steady-state, quasi-steady state, and transient profiles; optimum facility pipe and pipe component sizes and parameters; momentum transfer loads; and fluid conditions at critical points. A set of COSMIC NASTRAN-based finite element models is also presented to evaluate the loads and stresses on test facility piping systems from fluid and gaseous effects, thermal chill down, and occasional wind loads. The models are based on Apple Macintosh software which makes it possible to change numerous parameters.

Dequay, L.↗
NASA NTRS1991

Aerothermodynamic environment predictions in support of the Aeroassist Flight Experiment

A computational fluid dynamic (CFD) simulation has been made of the flow past an Aeroassist Flight Experiment (AFE) aerobrake model with a payload and sting. The simulation involves solving the complete Navier-Stokes equations in three dimensions. Comparisons with data taken in a Mach 10 wind tunnel test are made and agreement is shown to be very good. Points of comparison include the aerobrake surface pressure and heat transfer rate, the sting surface streamlines, and the sting heat transfer rate. The dependence of the solution on the grid is also explored. Finally, predictions are made for conditions which correspond to a hypersonic shock tunnel force and moment test to be carried out later this year.

Fay, John F.↗
NASA NTRS1997

Efficient Parallel Kernel Solvers for Computational Fluid Dynamics Applications

Distributed-memory parallel computers dominate today's parallel computing arena. These machines, such as Intel Paragon, IBM SP2, and Cray Origin2OO, have successfully delivered high performance computing power for solving some of the so-called "grand-challenge" problems. Despite initial success, parallel machines have not been widely accepted in production engineering environments due to the complexity of parallel programming. On a parallel computing system, a task has to be partitioned and distributed appropriately among processors to reduce communication cost and to attain load balance. More importantly, even with careful partitioning and mapping, the performance of an algorithm may still be unsatisfactory, since conventional sequential algorithms may be serial in nature and may not be implemented efficiently on parallel machines. In many cases, new algorithms have to be introduced to increase parallel performance. In order to achieve optimal performance, in addition to partitioning and mapping, a careful performance study should be conducted for a given application to find a good algorithm-machine combination. This process, however, is usually painful and elusive. The goal of this project is to design and develop efficient parallel algorithms for highly accurate Computational Fluid Dynamics (CFD) simulations and other engineering applications. The work plan is 1) developing highly accurate parallel numerical algorithms, 2) conduct preliminary testing to verify the effectiveness and potential of these algorithms, 3) incorporate newly developed algorithms into actual simulation packages. The work plan has well achieved. Two highly accurate, efficient Poisson solvers have been developed and tested based on two different approaches: (1) Adopting a mathematical geometry which has a better capacity to describe the fluid, (2) Using compact scheme to gain high order accuracy in numerical discretization. The previously developed Parallel Diagonal Dominant (PDD) algorithm and Reduced Parallel Diagonal Dominant (RPDD) algorithm have been carefully studied on different parallel platforms for different applications, and a NASA simulation code developed by Man M. Rai and his colleagues has been parallelized and implemented based on data dependency analysis. These achievements are addressed in detail in the paper.

Sun, Xian-He↗
NASA NTRS1999

Plane Smoothers for Multiblock Grids: Computational Aspects

Standard multigrid methods are not well suited for problems with anisotropic discrete operators, which can occur, for example, on grids that are stretched in order to resolve a boundary layer. One of the most efficient approaches to yield robust methods is the combination of standard coarsening with alternating-direction plane relaxation in the three dimensions. However, this approach may be difficult to implement in codes with multiblock structured grids because there may be no natural definition of global lines or planes. This inherent obstacle limits the range of an implicit smoother to only the portion of the computational domain in the current block. This report studies in detail, both numerically and analytically, the behavior of blockwise plane smoothers in order to provide guidance to engineers who use block-structured grids. The results obtained so far show alternating-direction plane smoothers to be very robust, even on multiblock grids. In common computational fluid dynamics multiblock simulations, where the number of subdomains crossed by the line of a strong anisotropy is low (up to four), textbook multigrid convergence rates can be obtained with a small overlap of cells between neighboring blocks.

Llorente, Ignacio M.↗
NASA NTRS2000

ICEG2D: An Integrated Software Package for Automated Prediction of Flow Fields for Single-Element Airfoils with Ice Accretion

An integrated software package, ICEG2D, was developed to automate computational fluid dynamics (CFD) simulations for single-element airfoils with ice accretion. ICEG2D is designed to automatically perform three primary functions: (1) generating a grid-ready, surface definition based on the geometrical characteristics of the iced airfoil surface, (2) generating a high-quality grid using the generated surface point distribution, and (3) generating the input and restart files needed to run the general purpose CFD solver NPARC. ICEG2D can be executed in batch mode using a script file or in an interactive mode by entering directives from a command line. This report summarizes activities completed in the first year of a three-year research and development program to address issues related to CFD simulations for aircraft components with ice accretion. Specifically, this document describes the technology employed in the software, the installation procedure, and a description of the operation of the software package. Validation of the geometry and grid generation modules of ICEG2D is also discussed.

Thompson, David S.↗