Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computational Fluid Dynamic (CFD)”

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,243 records · Page 69
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Development of Segregated Thermal-Hydraulics Solvers in MOOSE

The simulation of fluid flows is an essential part of the design and analysis of nuclear systems. Algorithms able to simulate flows at different fidelity levels are available in the Multiphysics Object-Oriented Simulation Environment (MOOSE) and MOOSE-based applications such as Pronghorn \cite{novak2018pronghorn}, Pronghorn-Subchannel, RELAP-7, and SAM. Currently, significant effort is being invested in the development of coarse-mesh Computational Fluid Dynamics (CFD) capabilities within MOOSE and Pronghorn for the simulation of Generation IV nuclear reactors. Traditionally, the solution algorithms in MOOSE have relied on Newton or quasi-Newton methods (such as the preconditioned Jacobian-free Newton-Krylov method) where residuals and Jacobians (or approximations thereof) are constructed. Both Newton and quasi-Newton methods require the solution of a linear system at each nonlinear Newton iteration with the Jacobian as the system matrix. The Jacobian contains blocks originating from all variables in the problem (i.e., for thermal-hydraulics at least pressure, velocities, and temperature). Due to the formulation of the problem in a general multiphysics setting on unstructured mesh, creating a good preconditioner for the linear system can be challenging, thus many fluid applications have utilized direct solver-based methods such as LU factorization. However, with increasing system size and complexity in multi-dimensional problems, the direct solution of linear systems becomes computationally expensive both in execution time and and memory. For this reason, recent effort has focused on adapting segregated solution algorithms for CFD problems in MOOSE. These algorithms use fixed-point iteration between segregated systems whose assembly and preconditioning are easier those of the monolithic system. Initial results show that the segregated solution algorithm outperforms the monolithic approach in terms of memory usage and for large 3D problems in terms of CPU time as well.

42 ENGINEERING↗
DOE OSTI2024

Development of Segregated Thermal-Hydraulics Solvers in MOOSE

The simulation of fluid flows is an essential part of the design and analysis of nuclear systems. Algorithms able to simulate flows at different fidelity levels are available in the Multiphysics Object-Oriented Simulation Environment (MOOSE) and MOOSE-based applications such as Pronghorn \cite{novak2018pronghorn}, Pronghorn-Subchannel, RELAP-7, and SAM. Currently, significant effort is being invested in the development of coarse-mesh Computational Fluid Dynamics (CFD) capabilities within MOOSE and Pronghorn for the simulation of Generation IV nuclear reactors. Traditionally, the solution algorithms in MOOSE have relied on Newton or quasi-Newton methods (such as the preconditioned Jacobian-free Newton-Krylov method) where residuals and Jacobians (or approximations thereof) are constructed. Both Newton and quasi-Newton methods require the solution of a linear system at each nonlinear Newton iteration with the Jacobian as the system matrix. The Jacobian contains blocks originating from all variables in the problem (i.e., for thermal-hydraulics at least pressure, velocities, and temperature). Due to the formulation of the problem in a general multiphysics setting on unstructured mesh, creating a good preconditioner for the linear system can be challenging, thus many fluid applications have utilized direct solver-based methods such as LU factorization. However, with increasing system size and complexity in multi-dimensional problems, the direct solution of linear systems becomes computationally expensive both in execution time and and memory. For this reason, recent effort has focused on adapting segregated solution algorithms for CFD problems in MOOSE. These algorithms use fixed-point iteration between segregated systems whose assembly and preconditioning are easier those of the monolithic system. Initial results show that the segregated solution algorithm outperforms the monolithic approach in terms of memory usage and for large 3D problems in terms of CPU time as well.

42 ENGINEERING↗
DOE OSTI2020

Multi-Scale Thermal Analysis for Design of SiC-Based Medium Voltage Motor Drive: Preprint

This paper describes a multi-scale thermal analysis approach for the design of an air-cooled 1.7-kV SiC MOSFET-based medium-voltage variable-speed motor drive. The National Renewable Energy Laboratory supported the project objectives by providing thermal management design support to the project partner The Ohio State University (OSU). The scope of the models and limited resources required efficient and flexible thermal models to be developed. Two modeling techniques are described that significantly reduced model run time and enabled more complex models to be run faster while retaining needed accuracy. The first technique uses the effectiveness-NTU method to extract convection boundary conditions from a computational fluid dynamics (CFD) model that can be applied to a fast-running finite element analysis model. The second is a porous media technique that enables system-level CFD simulations that incorporate effects from heat exchangers (e.g., pin fin heat sinks) that run in a fraction of the time required for fully resolved CFD simulations. The multi-scale approach to the thermal analysis enabled fast and accurate simulation for the converter design ranging from the die level up to the full system with 36 submodules. The modeling results were validated against experimental data from system tests performed by OSU.

DIRECT ENERGY CONVERSION↗
NASA NTRS2014

Long Hole Film Cooling Dataset for CFD Development - Flow and Film Effectiveness

An experiment investigating flow and heat transfer of long (length to diameter ratio of 18) cylindrical film cooling holes has been completed. In this paper, the thermal field in the flow and on the surface of the film cooled flat plate is presented for nominal freestream turbulence intensities of 1.5 and 8 percent. The holes are inclined at 30 deg above the downstream direction, injecting chilled air of density ratio 1.0 onto the surface of a flat plate. The diameter of the hole is 0.75 in. (approx. 0.02 m) with center to center spacing (pitch) of 3 hole diameters. Coolant was injected into the mainstream flow at nominal blowing ratios of 0.5, 1.0, 1.5, and 2.0. The Reynolds number of the freestream was approximately 11,000 based on hole diameter. Thermocouple surveys were used to characterize the thermal field. Infrared thermography was used to determine the adiabatic film effectiveness on the plate. Hotwire anemometry was used to provide flowfield physics and turbulence measurements. The results are compared to existing data in the literature. The aim of this work is to produce a benchmark dataset for Computational Fluid Dynamics (CFD) development to eliminate the effects of hole length to diameter ratio and to improve resolution in the near-hole region. In this report, a Time Filtered Navier Stokes (TFNS), also known as Partially Resolved Navier Stokes (PRNS), method that was implemented in the Glenn-HT code is used to model coolant-mainstream interaction. This method is a high fidelity unsteady method that aims to represent large scale flow features and mixing more accurately.

Shyam, Vikram↗
NASA NTRS2013

Long Hole Film Cooling Dataset for CFD Development : Infrared Thermography and Thermocouple Surveys - Part 1

An experiment investigating flow and heat transfer of long (length to diameter ratio of 18) cylindrical film cooling holes has been completed. In this paper, the thermal field in the flow and on the surface of the film cooled flat plate is presented for nominal freestream turbulence intensities of 1.5 and 8 percent. The holes are inclined at 30deg above the downstream direction, injecting chilled air of density ratio 1.0 onto the surface of a flat plate. The diameter of the hole is 0.75 in. (0.01905 m) with center to center spacing (pitch) of 3 hole diameters. Coolant was injected into the mainstream flow at nominal blowing ratios of 0.5, 1.0, 1.5, and 2.0. The Reynolds number of the freestream was approximately 11,000 based on hole diameter. Thermocouple surveys were used to characterize the thermal field. Infrared thermography was used to determine the adiabatic film effectiveness on the plate. Hotwire anemometry was used to provide flowfield physics and turbulence measurements. The results are compared to existing data in the literature. The aim of this work is to produce a benchmark dataset for Computational Fluid Dynamics (CFD) development to eliminate the effects of hole length to diameter ratio and to improve resolution in the near-hole region. In this report, a Time-Filtered Navier Stokes (TFNS), also known as Partially Resolved Navier Stokes (PRNS), method that was implemented in the Glenn-HT code is used to model coolant-mainstream interaction. This method is a high fidelity unsteady method that aims to represent large scale flow features and mixing more accurately.

Shyam, Vikram↗
NASA NTRS1996

Application of Two-Phase CFD to the Design and Analysis of a Subscale Motor Experiment to Evaluate Propellant Slag Production

The redesigned solid rocket motor (RSRM) Pressure Perturbation Investigation Team concluded that the cause of recent pressure spikes during both static and flight motor burns was the expulsion of molten aluminum oxide slag from a pool which collects in the aft end of the motor around the submerged nozzle nose during the last half of motor operation. It is suspected that some motors produce more slag than others due to differences in aluminum oxide agglomerate particle sizes which may relate to subtle differences in propellant ingredient characteristics such as particle size distribution, contaminants, or processing variations. In order to determine the effect of suspect propellant ingredient characteristics on the propensity for slag production in a real motor environment, a subscale motor experiment was designed. An existing 5 inch ballistic test motor was selected as the basic test vehicle due to low cost and quick turn around times. The standard converging/diverging nozzle was replaced with a submerged nozzle nose design to provide a positive trap for the slag which would increase both the quantity and repeatability of measured slag weights. Computational fluid dynamics (CFD) was used to assess a variety of submerged nose configurations to identify the design which possessed the best capability to reliably collect slag. CFD also was used to assure that the final selected nozzle design would result in flow field characteristics such as dividing streamline location, nose attach point, and separated flow structure which would have similtude with the RSRM submerged nozzle nose flow field. It also was decided to spin the 5 inch motor about its longitudinal axis to further enhance slag collection quantities. Again, CFD was used to select an appropriate spin rate along with other considerations, including the avoidance of burn rate enhancement from radial acceleration effects.

Whitesides, R. Harold↗
NASA NTRSDate not supplied

CFD Modeling & Validation of Single Bubble Growth During Cavitation in Zero Gravity

Cryogenic Fluid Management (CFM) in zero gravity represents a vital area of research as NASA pursues future long-duration space missions. Success of these missions relies heavily on the ability to store and transfer the cryogenic liquid fuels efficiently and safely. In the Zero Boil-Off Tank (ZBOT-1) Experiment in 2018, unexpected and intense phase change occurred at the screened Liquid Acquisition Device (LAD) during tank pressure control that was later attributed to cavitation. The cavitation ultimately resulted in pump failure. Since cryogen cavitation during storage and transfer of propellants in microgravity has been a source of serious concern for the aerospace community, a series of drop tower experiments with liquid methane were also performed at the University of Bremen to carefully study cavitation from a single bubble nucleus in an engineered micron-sized wall cavity. As the test cell is depressurized during the 10-second drop, cavitation causes bubble growth from the nucleus. Here, in a collaboration between the Case Western Reserve University and University of Bremen, we present a two-phase Computational Fluid Dynamics (CFD) model of bubble growth from a single wall cavity during depressurization in microgravity. Both 2D axisymmetric and 3D parametric simulation case studies are performed to show the evolution of the bubble growth. The model uses a Volume of Fluid (VOF) approach with and without adaptive mesh refinement and adaptive time stepping to capture the growth dynamics using an explicit formulation with algebraic reconstruction of the interface. CFD results show good agreement with corresponding experiment. Based on the velocity, temperature, and volume fraction fields generated by the validation case studies, detailed discussion and analysis of fluid flow and heat transfer during bubble growth process are presented to delineate the intricacies of the cavitation process. Strengths and shortcomings of the present numerical approach in capturing this phenomenon are also highlighted.

Two Phase Flow↗
NASA NTRS2008

Software Performs Complex Design Analysis

Designers use computational fluid dynamics (CFD) to gain greater understanding of the fluid flow phenomena involved in components being designed. They also use finite element analysis (FEA) as a tool to help gain greater understanding of the structural response of components to loads, stresses and strains, and the prediction of failure modes. Automated CFD and FEA engineering design has centered on shape optimization, which has been hindered by two major problems: 1) inadequate shape parameterization algorithms, and 2) inadequate algorithms for CFD and FEA grid modification. Working with software engineers at Stennis Space Center, a NASA commercial partner, Optimal Solutions Software LLC, was able to utilize its revolutionary, one-of-a-kind arbitrary shape deformation (ASD) capability-a major advancement in solving these two aforementioned problems-to optimize the shapes of complex pipe components that transport highly sensitive fluids. The ASD technology solves the problem of inadequate shape parameterization algorithms by allowing the CFD designers to freely create their own shape parameters, therefore eliminating the restriction of only being able to use the computer-aided design (CAD) parameters. The problem of inadequate algorithms for CFD grid modification is solved by the fact that the new software performs a smooth volumetric deformation. This eliminates the extremely costly process of having to remesh the grid for every shape change desired. The program can perform a design change in a markedly reduced amount of time, a process that would traditionally involve the designer returning to the CAD model to reshape and then remesh the shapes, something that has been known to take hours, days-even weeks or months-depending upon the size of the model.

Source record↗
NASA NTRS2014

Film Cooled Recession of SiC/SiC Ceramic Matrix Composites: Test Development, CFD Modeling and Experimental Observations

SiCSiC ceramic matrix composites (CMCs) systems will play a crucial role in next generation turbine engines for hot-section component applications because of their ability to significantly increase engine operating temperatures, reduce engine weight and cooling requirements. However, the environmental stability of Si-based ceramics in high pressure, high velocity turbine engine combustion environment is of major concern. The water vapor containing combustion gas leads to accelerated oxidation and corrosion of the SiC based ceramics due to the water vapor reactions with silica (SiO2) scales forming non-protective volatile hydroxide species, resulting in recession of the ceramic components. Although environmental barrier coatings are being developed to help protect the CMC components, there is a need to better understand the fundamental recession behavior of in more realistic cooled engine component environments.In this paper, we describe a comprehensive film cooled high pressure burner rig based testing approach, by using standardized film cooled SiCSiC disc test specimen configurations. The SiCSiC specimens were designed for implementing the burner rig testing in turbine engine relevant combustion environments, obtaining generic film cooled recession rate data under the combustion water vapor conditions, and helping developing the Computational Fluid Dynamics (CFD) film cooled models and performing model validation. Factors affecting the film cooled recession such as temperature, water vapor concentration, combustion gas velocity, and pressure are particularly investigated and modeled, and compared with impingement cooling only recession data in similar combustion flow environments. The experimental and modeling work will help predict the SiCSiC CMC recession behavior, and developing durable CMC systems in complex turbine engine operating conditions.

Recession↗
NASA NTRS2000

Mixing and NOx Emission Calculations of Confined Reacting Jet Flows in Cylindrical and Annular Ducts

Rapid mixing of cold lateral jets with hot cross-stream flows in confined configurations is of practical interest in gas turbine combustors as it strongly affects combustor exit temperature quality, and gaseous emissions in for example rich-lean combustion. It is therefore important to further improve our fundamental understanding of the important processes of dilution jet mixing especially when the injected jet mass flow rate exceeds that of the cross-stream. The results reported in this report describe some of the main flow characteristics which develop in the mixing process in a cylindrical duct. A three-dimensional computational fluid dynamics (CFD) code has been used to predict the mixing flow field characteristics and NOx emission in a quench section of a rich-burn/quick-mix/lean-burn (RQL) combustor. Sixty configurations have been analyzed in both circular and annular geometries in a fully reacting environment simulating the operating condition of an actual RQL gas turbine combustion liner. The evaluation matrix was constructed by varying the number of orifices per row and orifice shape. Other parameters such as J (momentum-flux ratio), MR (mass flowrate ratio), DR (density ratio), and mixer sector orifice ACd (effective orifice area) were maintained constant throughout the entire study. The results indicate that the mixing flow field can be correlated with the NOx production if they are referenced with the stoichiometric equivalence ratio value and not the equilibrium value. The mixing flowfields in both circular and annular mixers are different. The penetration of equal jets in both annular and circular geometries is vastly different which significantly affects the performance of the mixing section. In the computational results with the circular mixer, most of the NOx formation occurred behind the orifice starting at the orifice wake region. General trends have been observed in the NOx production as the number of orifices is changed and this appears to be common for all hole configurations and mixer types (circular or annular). The performance of any orifice shape (in producing minimum NOx) appears to be acceptable if the number of orifices can be freely varied in order to attain the optimum jet penetration.

Oechsle, Victor L.↗
DOE OSTI2024

Modeling Flow and Particle Deposition in a Spent Nuclear Fuel Assembly

CFD (Computational Fluid Dynamic) simulation of aerosol-laden natural convective flow and particle deposition in a spent fuel storage canister with 37 assemblies is currently computationally prohibitive. PWR (Pressurized Water Reactor) assemblies have up to 289 pins or tubes with several spacer grids to align the pins. Spacer grids with mixing vanes induce swirling during operation to increase heat transfer. Each spacer grid contains hundreds of small structures such as retaining clips, channel walls, and openings. The largest canisters store 37 PWR assemblies thus, there are numerous pins, tubes, and spacer grids for which the flow region between and around these structures needs to be determined along with the movement and deposition of aerosol particles. Because of the complicated geometry, modeling the intricate flow even for just one assembly is currently impractical. Nonetheless, we are developing techniques for a practical model to assess the natural aerosol particle deposition process in a canister in the event that a release occurs from one or more fuel pins. In the previous work it was demonstrated that CFD can model the flow through a PWR spacer grid with mixing vanes, including particle deposition, in a reasonable amount of time on a personal computer. In this work, the analysis is extended to include the bypass region between an assembly and the canister basket walls. It is shown that the flow velocity in the bypass region is about three times that of the interstitial region between the pins. The lengths before and after the spacer grid are also extended to determine when the flow becomes fully developed. In addition, the approach of computationally “stitching together” segments of an assembly is demonstrated with the plan to ultimately model a full assembly. The fraction of particles that are deposited in a segment with a spacer grid is determined as a function of particle size and flow velocity.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

ENGINE COMBUSTION SYSTEM OPTIMIZATION USING CFD AND MACHINE LEARNING: A METHODOLOGICAL APPROACH

Gasoline compression ignition (GCI) engines are considered an attractive alternative to traditional spark-ignition and diesel engines. In this work, a Machine Learning-Grid Gradient Algorithm (ML-GGA) approach was developed to optimize the performance of internal combustion engines. Machine learning (ML) offers a pathway to transform complex physical processes that occur in a combustion engine into compact informational processes. The developed ML-GGA model was compared with a recently developed Machine learning Genetic Algorithm (ML-GA). Detailed investigations of optimization solver parameters and variables limits extension were performed in the present ML-GGA model to improve the accuracy and robustness of the optimization process. Detailed descriptions of the different procedures, optimization tools and criteria that must be followed for a successful output are provided here. The developed ML-GGA approach was used to optimize the operating conditions (case 1) and the piston bowl design (case 2) of a heavy-duty diesel engine running on a gasoline fuel with a Research Octane Number (RON) of 80. The ML-GGA approach yielded > 2% improvements in the merit function, compared to the optimum obtained from a thorough computational fluid dynamics (CFD) guided system optimization. The predictions from the ML-GGA approach were validated with engine CFD simulations. This study demonstrates the potential of ML-GGA to significantly reduce the time needed for optimization problems, without loss in accuracy compared to traditional approaches.

engine optimization↗
NASA NTRS2003

Numerical Models of Human Circulatory System under Altered Gravity: Brain Circulation

A computational fluid dynamics (CFD) approach is presented to model the blood flow through the human circulatory system under altered gravity conditions. Models required for CFD simulation relevant to major hemodynamic issues are introduced such as non-Newtonian flow models governed by red blood cells, a model for arterial wall motion due to fluid-wall interactions, a vascular bed model for outflow boundary conditions, and a model for auto-regulation mechanism. The three-dimensional unsteady incompressible Navier-Stokes equations coupled with these models are solved iteratively using the pseudocompressibility method and dual time stepping. Moving wall boundary conditions from the first-order fluid-wall interaction model are used to study the influence of arterial wall distensibility on flow patterns and wall shear stresses during the heart pulse. A vascular bed modeling utilizing the analogy with electric circuits is coupled with an auto-regulation algorithm for multiple outflow boundaries. For the treatment of complex geometry, a chimera overset grid technique is adopted to obtain connectivity between arterial branches. For code validation, computed results are compared with experimental data for steady and unsteady non-Newtonian flows. Good agreement is obtained for both cases. In sin-type Gravity Benchmark Problems, gravity source terms are added to the Navier-Stokes equations to study the effect of gravitational variation on the human circulatory system. This computational approach is then applied to localized blood flows through a realistic carotid bifurcation and two Circle of Willis models, one using an idealized geometry and the other model using an anatomical data set. A three- dimensional anatomical Circle of Willis configuration is reconstructed from human-specific magnetic resonance images using an image segmentation method. The blood flow through these Circle of Willis models is simulated to provide means for studying gravitational effects on the brain circulation under auto-regulation.

Kim, Chang Sung↗
NASA NTRS2014

CFD Analysis of Thermal Control System Using NX Thermal and Flow

The Thermal Control Subsystem (TCS) is a key part of the Advanced Plant Habitat (APH) for the International Space Station (ISS). The purpose of this subsystem is to provide thermal control, mainly cooling, to the other APH subsystems. One of these subsystems, the Environmental Control Subsystem (ECS), controls the temperature and humidity of the growth chamber (GC) air to optimize the growth of plants in the habitat. The TCS provides thermal control to the ECS with three cold plates, which use Thermoelectric Coolers (TECs) to heat or cool water as needed to control the air temperature in the ECS system. In order to optimize the TCS design, pressure drop and heat transfer analyses were needed. The analysis for this system was performed in Siemens NX Thermal/Flow software (Version 8.5). NX Thermal/Flow has the ability to perform 1D or 3D flow solutions. The 1D flow solver can be used to represent simple geometries, such as pipes and tubes. The 1D flow method also has the ability to simulate either fluid only or fluid and wall regions. The 3D flow solver is similar to other Computational Fluid Dynamic (CFD) software. TCS performance was analyzed using both the 1D and 3D solvers. Each method produced different results, which will be evaluated and discussed.

CFD↗
NASA NTRS1992

Application of computational fluid dynamics to the design of the film cooled STME subscale nozzle for the National Launch System

The status of computational fluid dynamics (CFD) calculations for the Space Transportation Main Engine (STME) film/dump cooled nozzle design is presented, with an emphasis on the timely impact of CFD on the design of the sub-scale nozzle coolant system. The following aspects of the sub-scale coolant delivery system were analyzed with CFD: 1) a design trade study of a mechanical flow splitting device for uniform distribution of the subsonic cavity flow, 2) a design trade study of the subsonic cavity lip to achieve film integrity, and 3) an analysis of the primary flow interaction with the core/secondary coolant streams. All design calculations were performed with the Generalized Aerodynamic Simulation Program (GASP), a 3-D, multi-block, generalized Navier-Stokes code capable of solving with frozen, finite-rate or equilibrium chemical kinetics. The initial design of the subsonic cavity flow used square posts to distribute the sonic orifice jets into a uniform flow. Calculations for this design indicated that an unacceptable mal-distribution of film occurred. Design modifications involving curved and slotted posts were computed in an effort to uniformly distribute the secondary coolant flow. Analysis of these configurations showed that although the flowfield improved in uniformity, it was still unacceptable, especially at higher feed pressures. Results from these studies were then incorporated into a design that resulted in the insertion of a porous metal ring into the subsonic cavity. Subsequent water flow model studies showed that this concept was successful in uniformly distributing flow exiting the cavity. In addition to the design of the subsonic cavity, CFD was also used to analyze the secondary coolant lip and the primary flow interaction with the core/secondary coolant streams. A series of calculations were first performed to modify the subsonic cavity lip contour. The flow over the modified lip was then computed simultaneously with the primary injectors to determine the impact of the subsonic coolant stream on the primary slot jets. Pressure, temperature, velocity, and coolant mass fraction contours are presented for these configurations. Information is given in viewgraph form.

Garrett, Joseph L.↗
DOE OSTI2020

Application of GOTHIC to Groundwater Transport Analysis - 20152

Migration of pollutants and hazardous wastes, potentially containing radioactive isotopes, via groundwater transport is a concern at most waste cleanup sites. Predictive analysis can be used to evaluate mitigating actions intended to minimize impact on the environment and public exposure. GOTHIC is a multipurpose thermal hydraulics code that is used extensively in the nuclear industry for design, licensing and operation evaluations. It combines the capabilities of typical one-dimensional system codes and the essential features of Computational Fluid Dynamics (CFD) codes for three-dimensional analysis. There are other codes that are specifically developed for groundwater transport analysis and the results presented here are consistent with prior analyses. However, GOTHIC has some unique features that offer advantages for applications related to nuclear waste. Most importantly, it has been developed and maintained under a Quality Assurance program in compliance with the requirements of 10CFR50 Appendix B [2] and applicable portions of ASME NQA-1 [3] since 1995. Available GOTHIC capabilities that make the code especially useful for groundwater transport of nuclear materials include: - Tracking of any number of tracer elements for contaminants and other species of interest; - Radioactive decay and progeny of tracer elements; - Adsorption/desorption of tracer elements; - Tracking of any number of dissolved gases; - Release and absorption for dissolved gases; - Vapor phase tracking; - Non-Newtonian fluid modeling. The general porous body modeling approach makes GOTHIC well suited to groundwater transport analysis. The multi-region modeling approach used by GOTHIC simplifies model construction for regions of varying hydrologic characteristics and focuses the computational effort on regions of particular interest while simultaneously capturing the macroscopic response and any feedback effects across the larger domain. The applicability of GOTHIC to groundwater transport applications is demonstrated by comparing code results with available analytic or semi-analytic solutions for groundwater behavior. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS1998

Virtual and Experimental Visualization of Flows in Packed Beds of Spheres Simulating Porous Media Flows

A videotape presentation of flow in a packed bed of spheres is provided. The flow experiment consisted of three principal elements: (1) an oil tunnel 76.2 mm by 76.2 mm in cross section, (2) a packed bed of spheres in regular and irregular arrays, and (3) a flow characterization methodology, either (a) full flow field tracking (FFFT) or (b) computational fluid dynamic (CFD) simulation. The refraction indices of the oil and the test array of spheres were closely matched, and the flow was seeded with aluminum oxide particles. Planar laser light provided a two-dimensional projection of the flow field, and a traverse simulated a three-dimensional image of the entire flow field. Light focusing and reflection rendered the spheres black, permitting visualization of the planar circular interfaces in both the axial and transverse directions. Flows were observed near the wall-sphere interface and within the set of spheres. The CFD model required that a representative section of a packed bed be formed and gridded, enclosing and cutting six spheres so that symmetry conditions could be imposed at all cross-boundaries. Simulations had to be made with the flow direction at right angles to that used in the experiments, however, to take advantage of flow symmetry. Careful attention to detail was required for proper gridding. The flow field was three-dimensional and complex to describe, yet the most prominent finding was flow threads, as computed in the representative 'cube' of spheres with face symmetry and conclusively demonstrated experimentally herein. Random packing and bed voids tended to disrupt the laminar flow, creating vortices.

Hendricks, R. C.↗
NASA NTRS1996

Aerodynamic shape optimization using control theory

Aerodynamic shape design has long persisted as a difficult scientific challenge due its highly nonlinear flow physics and daunting geometric complexity. However, with the emergence of Computational Fluid Dynamics (CFD) it has become possible to make accurate predictions of flows which are not dominated by viscous effects. It is thus worthwhile to explore the extension of CFD methods for flow analysis to the treatment of aerodynamic shape design. Two new aerodynamic shape design methods are developed which combine existing CFD technology, optimal control theory, and numerical optimization techniques. Flow analysis methods for the potential flow equation and the Euler equations form the basis of the two respective design methods. In each case, optimal control theory is used to derive the adjoint differential equations, the solution of which provides the necessary gradient information to a numerical optimization method much more efficiently then by conventional finite differencing. Each technique uses a quasi-Newton numerical optimization algorithm to drive an aerodynamic objective function toward a minimum. An analytic grid perturbation method is developed to modify body fitted meshes to accommodate shape changes during the design process. Both Hicks-Henne perturbation functions and B-spline control points are explored as suitable design variables. The new methods prove to be computationally efficient and robust, and can be used for practical airfoil design including geometric and aerodynamic constraints. Objective functions are chosen to allow both inverse design to a target pressure distribution and wave drag minimization. Several design cases are presented for each method illustrating its practicality and efficiency. These include non-lifting and lifting airfoils operating at both subsonic and transonic conditions.

Reuther, James↗