Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computational geometry”

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 721 records · Page 40

Effect of geometry on airfoil icing characteristics

A droplet trajectory computer code is used to predict the water droplet impingement characteristics of several low- and medium-speed airfoils. The maximum impingement efficiency, total collection efficiency, and limits of impingement are analyzed as functions of the airfoil geometry and freestream conditions. The airfoil geometry is represented by leading edge radius, maximum thickness, maximum camber, and angle of attack. The analysis shows that the primary effects are an increase in maximum impingement efficiency with a decrease in leading edge radius, a reduction in total collection efficiency for thicker airfoils, and a change in the limits of impingement for airfoils of different maximum camber.

Bragg, M. B.↗

Multiphase lattice Boltzmann simulations for porous media applications

Over the last two decades, lattice Boltzmann methods have become an increasingly popular tool to compute the flow in complex geometries such as porous media. In addition to single phase simulations allowing, for example, a precise quantification of the permeability of a porous sample, a number of extensions to the lattice Boltzmann method are available which allow to study multiphase and multicomponent flows on a pore scale level. In this work, we give an extensive overview on a number of these diffuse interface models and discuss their advantages and disadvantages. Furthermore, we shortly report on multiphase flows containing solid particles, as well as implementation details and optimization issues.

42 ENGINEERING↗

Thermally Effective and Efficient Cooling Technologies for Advanced Gas Turbine Systems

This report presents the research which has been undertaken under DOE grant DE-FE-0011875 as part of NETL’s UTSR (University Turbine System Research) Program. The experimental portion of this work was undertaken at the University of North Dakota (UND) under the direction of Dr. Forrest Ames. The computational work was undertaken at LSU (Louisiana State University), UM (University of Memphis), and at IIT (Illinois Institute of Technology) under the direction of Dr. Sumanta Acharya. The research program was organized into a four phase program which included: 1) development of a FDA/FEA (finite difference analysis/finite element analysis) cooling model, 2) experimental/computational investigation of appropriate internal cooling methods, 3) experimental/computational investigation of external heat transfer and film cooling, 4) development of an enhanced FDA/FEA model along with the development of a warm cascade test for comparison purposes. The experimental work at UND included research on both internal cooling methods and external heat transfer. The internal cooling work included the design and testing of multiple variable hole-size “incremental impingement” configurations along with testing of two converging rounded diamond pedestal arrays for the trailing edge region. The external heat transfer work included the building of a large scale cascade facility to test a vane design with a large leading edge and an aft loaded suction surface. UND has acquired both midspan and full suction surface heat transfer measurements as well as prepared for film cooling measurements under the grant. LSU/UM/IIT has developed and conducted wall resolved LES computations for the incremental impingement geometries and provided rational understanding of the effect of hole sizes and their distributions. Additionally they have provided design correlations for the incremental impingement heat transfer coefficients. They have also developed LES methods to predict flow through film cooling plenums, which have pedestals along with the development of slot film cooling downstream from the plenums. Additionally, they have developed wall resolved LES to predict vane surface heat transfer including bypass transition without the use of a transition model. Through these simulations they provide understanding of the aerodynamic and surface heat transfer behavior under a range of influences such as freestream turbulence. The present study has produced some very useful research outcomes.

01 COAL, LIGNITE, AND PEAT↗

Analytical investigation of rotor wake formation and geometry

A number of refinements in the computer code were worked out and tested. Three codes have been written to date. One program is for an isolated wing and is being used to compare with data for the vortex wake (Weston). The second code is for an isolated wing with a streamwise vortex passing above it. This program is being used to validate the computational procedure for incorporating the vortex into the Euler equation calculations. The third program is the hovering rotor code which is the overall objective of the research. The optimization calculations for a hovering helicopter rotor have been completed.

Miller, R.↗

Compilation of methods in orbital mechanics and solar geometry

This paper contains a collection of computational algorithms for determining geocentric ephemerides of Earth satellites, useful for both mission planning and data reduction applications. Special emphasis is placed on the computation of sidereal time, and on the determination of the geocentric coordinate of the center of the Sun, all to the accuracy found in the Astronomical Almanac. The report is completely self-contained in that no requirement is placed on any external source of information, and hence, these methods are ideal for computer application.

Buglia, James J.↗

Unsteady rotor aerodynamics using a vortex panel method

The problem of the potential incompressible flow about a helicopter rotor blade is solved using an unsteady vortex-panel method where the mutual interaction between the blade and the distorting free wake is taken into account. The present method alleviates the need to rely upon measured-wake geometries or p5escribed-wake models in order to calculate the airloads. A computer program has been developed which is capable of predicting the geometry of the time-dependent three-dimensional (3-D) wake and the instantaneous loadings for a single blade in hover, climb, and forward flight. The solution is obtained by using a time-accurate step-by-step procedure. The complex-wake geometry at any time is presented graphically with a computer graphics system. Calculated results are compared with published data for a rotor blade in both hover and forward flight. The code has also been applied to the study of the effect of changing blade tip geometry.

Crispin, Y.↗

Lobed Mixer Optimization for Advanced Ejector Geometries

The overall objectives are: 1) to pursue analytical, computational, and experimental studies that enhance basic understanding of forced mixing phenomena relevant to supersonic jet noise reduction, and 2) to integrate this enhanced understanding (analytical, computational, and empirical) into a design-oriented model of a mixer-ejector noise suppression system. The work is focused on ejector geometries and flow conditions typical of those being investigated in the NASA High Speed Research Program (HSRP). The research will be carried out in collaboration with the NASA HSRP Nozzle Integrated Technology Development (ITD) Team, and will both contribute to, and benefit from, the results of other HSRP research. The noise suppressor system model that is being developed under this grant is distinct from analytical tools developed by industry because it directly links details of lobe geometry to mixer-ejector performance. In addition, the model provides a 'technology road map to define gaps in the current understanding of various phenomena related to mixer-ejector design and to help prioritize research areas. This report describes research completed in the past year, as well as work proposed for the following year.

Waitz, Ian A.↗

Code Coverage Status of the ARC Code PERSENT

The Argonne Reactor Code (ARC) software system supports users in their fast reactor design goals by providing neutronic, thermal-hydraulic, and structural analysis capabilities. PERSENT fulfills the role of generating reactivity coefficients for a given time point of a REBUS calculation usable in a point kinetics based safety analysis capability. PERSENT also provides a sensitivity coefficient capability on eigenvalue, reactivity worth, and several other key coefficients that are used in the follow-on safety analysis. Given a co-variance matrix, PERSENT can carry out the uncertainty quantification to indicate the amount of error in the reactivity coefficients derived from the errors in the cross section measurements. With continued improvement of computational resources, many of the geometry modeling capabilities in DIF3D that were primarily used in low order schemes are not really needed anymore. Today, the diffusion and transport capabilities of DIF3D-VARIANT are primarily used in the reactor design process with some scattered usage of DIF3D-FD and DIF3D-Nodal. PERSENT is part of the ARC code system and is built around DIF3D-VARIANT and the flux solution it provides. The purpose of the present work is to identify a set of test problems for PERSENT and assess the code coverage of PERSENT for those test problems. PERSENT treats the DIF3D executable as an external executable and thus the code coverage considerations only need to focus on the PERSENT source code and only a fraction of the connected modules in the existing ARC software library. The goal is to document what parts of the existing PERSENT code are touched by the set of test problems and which are not. Because the verification work done on PERSENT was focused on the most common uses of PERSENT for fast reactor analysis, the code coverage assessment of those capabilities is the highest priority. This will ensure that nothing is being missed by the existing verification test problems that users of PERSENT rely upon. The code coverage analysis of PERSENT was performed with the Code Coverage Tool of the Intel Fortran compiler which requires modifications to the compilation of PERSENT. The detailed coverage tables are given for each submodule of PERSENT. Most of the uncovered parts/files could be easily ignored because they are either for error message and debugging output or not needed by PERSENT today. Only a few uncovered parts of PERSENT deserve extending the verification test suite.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Code Coverage Status of ARC Code-DIF3D

The Argonne Reactor Code (ARC) software system supports users in their fast reactor design goals by providing neutronic, thermal-hydraulic, and structural analysis capabilities. DIF3D plays a pivotal role in the ARC system as the primary homogenized assembly neutronic calculation methodology for fast reactor problems. Over its 40 years history, ARC software usage with DIF3D has been applied to numerous fast and thermal spectrum reactor analysis projects with good to excellent comparison against experiments. With continued improvement of computation resources, many of the geometry modeling capabilities in DIF3D that were primarily used in low order schemes are not really needed anymore. Today, the diffusion and transport capabilities of DIF3D-VARIANT are primarily used in the reactor design process with some scattered usage of DIF3D-FD and DIF3D-Nodal. In recent work, the DIF3D software verification was completed for DIF3D-FD and DIF3D-VARIANT on the geometry options used in the Versatile Test Reactor project. While we can be confident that these capabilities of DIF3D are well used and thus trusted, it does not demonstrate that all possible input options of DIF3D are actually working, but just those that were tested as part of VTR are and that they are correct. Thus, the purpose of the present work is to identify a set of test problems for DIF3D and assess the code coverage of DIF3D for those test problems. The goal is to document what parts of the existing DIF3D code are touched by the set of test problems and which are not. Because the verification work done on DIF3D-VARIANT and DIF3D-FD was focused on the most common uses of DIF3D for fast reactor analysis, the code coverage assessment of those capabilities is the highest priority. This will ensure that nothing is being missed by the existing verification test problems that DIF3D relies upon. The DIF3D-Nodal capability will also be inspected for code coverage as part of this work to further ensure that regular regression testing of DIF3D will trap any likely errors the end user might experience with the DIF3D software. The code coverage analysis of DIF3D was performed with the Code Coverage Tool of the Intel Fortran compiler which requires modifications to the compilation of DIF3D. The detailed coverage tables are given for each submodule of DIF3D separately, and for the submodules which are primarily developed for DIF3D, most of the source files could be at least partially touched. Most of the uncovered parts/files could be easily ignored, because they are either for error message and debugging output or obviously not needed by DIF3D. Out of the entire source codes of DIF3D, only a few uncovered modules deserve further investigation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A Computer Program for Calculating Three-Dimensional Compressible Laminar and Turbulent Boundary Layers on Arbitrary Wings

A computer program for calculating three dimensional compressible laminar and turbulent boundary layers on arbitrary wings is described and presented. The computer program consists of three separate programs, namely, a geometry program to represent the wing analytically, a velocity program to compute the external velocity components from a given experimental pressure distribution and a finite difference boundary layer method to solve the governing equations for compressible flows. To illustrate the usage of the computer program, three different test cases are presented and the preparation of the input data as well as the computed output data is discussed in some detail.

Cebeci, T.↗

Code Coverage Status of the ARC Code DIF3D

The Argonne Reactor Code (ARC) software system supports users in their fast reactor design goals by providing neutronic, thermal-hydraulic, and structural analysis capabilities. DIF3D plays a pivotal role in the ARC system as the primary homogenized assembly neutronic calculation methodology for fast reactor problems. Over its 40 years history, ARC software usage with DIF3D has been applied to numerous fast and thermal spectrum reactor analysis projects with good to excellent comparison against experiments. With continued improvement of computation resources, many of the geometry modeling capabilities in DIF3D that were primarily used in low order schemes are not really needed anymore. Today, the diffusion and transport capabilities of DIF3D-VARIANT are primarily used in the reactor design process with some scattered usage of DIF3D-FD and DIF3D-Nodal. In recent work, the DIF3D software verification was completed for DIF3D-FD and DIF3D-VARIANT on the geometry options used in the Versatile Test Reactor project. While we can be confident that these capabilities of DIF3D are well used and thus trusted, it does not demonstrate that all possible input options of DIF3D are actually working, but just those that were tested as part of VTR are and that they are correct. Thus, the purpose of the present work is to identify a set of test problems for DIF3D and assess the code coverage of DIF3D for those test problems. The goal is to document what parts of the existing DIF3D code are touched by the set of test problems and which are not. Because the verification work done on DIF3D-VARIANT and DIF3D-FD was focused on the most common uses of DIF3D for fast reactor analysis, the code coverage assessment of those capabilities is the highest priority. This will ensure that nothing is being missed by the existing verification test problems that DIF3D relies upon. The DIF3D-Nodal capability will also be inspected for code coverage as part of this work to further ensure that regular regression testing of DIF3D will trap any likely errors the end user might experience with the DIF3D software. The code coverage analysis of DIF3D was performed with the Code Coverage Tool of the Intel Fortran compiler which requires modifications to the compilation of DIF3D. The detailed coverage tables are given for each submodule of DIF3D separately, and for the submodules which are primarily developed for DIF3D, most of the source files could be at least partially touched. Most of the uncovered parts/files could be easily ignored, because they are either for error message and debugging output or obviously not needed by DIF3D. Out of the entire source codes of DIF3D, only a few uncovered modules deserve further investigation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Unstructured and adaptive mesh generation for high Reynolds number viscous flows

A method for generating and adaptively refining a highly stretched unstructured mesh suitable for the computation of high-Reynolds-number viscous flows about arbitrary two-dimensional geometries was developed. The method is based on the Delaunay triangulation of a predetermined set of points and employs a local mapping in order to achieve the high stretching rates required in the boundary-layer and wake regions. The initial mesh-point distribution is determined in a geometry-adaptive manner which clusters points in regions of high curvature and sharp corners. Adaptive mesh refinement is achieved by adding new points in regions of large flow gradients, and locally retriangulating; thus, obviating the need for global mesh regeneration. Initial and adapted meshes about complex multi-element airfoil geometries are shown and compressible flow solutions are computed on these meshes.

Mavriplis, Dimitri J.↗

A study of fracture mechanisms in ATD roller bearing

The purpose was to investigate how microstructures, especially anisotropy, affects internal stresses and the overall mechanical response of bearings. Samples with the stress axis along the aligned carbide direction possessed high modulus values compared to those with their axis perpendicular to the carbide stringers. The difference in the modulus was found to be more than a factor of two. A series of experiments was conducted on rolled samples to further investigate this effect; the two sets of results were consistent with each other. The degree of anisotropy of the microstructure in terms of the carbide and matrix orientations was determined using x-ray diffraction. The stress state determination was conducted using neutron diffraction. It was found that there was little variation in the distribution of the internal stresses amongst different samples, indicating that small changes in the processing and geometrical conditions did not result in significant variations in the internal stress. A nominal tensile hoop stress of 39 ksi was obtained for the inner raceway. Furthermore, during the course of the investigation, it was apparent that there was a need to determine the stress state induced by the shrink fitting process. Therefore, a series of experiments was conducted using strain gages to identify the stress distribution in a shrink fitting process in three different types of geometries. Correlations were obtained to estimate the highest stress values in the outer and inner groove geometry. A finite element program based on the ANSYS system was developed to compute the stress distribution in the inner raceway geometry. This analysis indicates that the highest tensile stress in the system occurs at the ID of the ring with a stress value of over 5 times that of the applied radial stress. Results from all these facets were correlated with one another. It appears that the material does not fail as a result of any one single factor, but results from a combination of the various factors investigated.

Zee, Ralph H.↗

Unstructured and adaptive mesh generation for high Reynolds number viscous flows

A method for generating and adaptively refining a highly stretched unstructured mesh suitable for the computation of high-Reynolds-number viscous flows about arbitrary two-dimensional geometries was developed. The method is based on the Delaunay triangulation of a predetermined set of points and employs a local mapping in order to achieve the high stretching rates required in the boundary-layer and wake regions. The initial mesh-point distribution is determined in a geometry-adaptive manner which clusters points in regions of high curvature and sharp corners. Adaptive mesh refinement is achieved by adding new points in regions of large flow gradients, and locally retriangulating; thus, obviating the need for global mesh regeneration. Initial and adapted meshes about complex multi-element airfoil geometries are shown and compressible flow solutions are computed on these meshes.

Mavriplis, D. J.↗

Transient Two-Dimensional Analysis of Side Load in Liquid Rocket Engine Nozzles

Two-dimensional planar and axisymmetric numerical investigations on the nozzle start-up side load physics were performed. The objective of this study is to develop a computational methodology to identify nozzle side load physics using simplified two-dimensional geometries, in order to come up with a computational strategy to eventually predict the three-dimensional side loads. The computational methodology is based on a multidimensional, finite-volume, viscous, chemically reacting, unstructured-grid, and pressure-based computational fluid dynamics formulation, and a transient inlet condition based on an engine system modeling. The side load physics captured in the low aspect-ratio, two-dimensional planar nozzle include the Coanda effect, afterburning wave, and the associated lip free-shock oscillation. Results of parametric studies indicate that equivalence ratio, combustion and ramp rate affect the side load physics. The side load physics inferred in the high aspect-ratio, axisymmetric nozzle study include the afterburning wave; transition from free-shock to restricted-shock separation, reverting back to free-shock separation, and transforming to restricted-shock separation again; and lip restricted-shock oscillation. The Mach disk loci and wall pressure history studies reconfirm that combustion and the associated thermodynamic properties affect the formation and duration of the asymmetric flow.

Wang, Ten-See↗

Turbulent Supersonic/Hypersonic Heating Correlations for Open and Closed Cavities

Supersonic/hypersonic laminar heating correlations that were developed for damage assessment analysis of atmospheric re-entry vehicles have been modified and extended to cover fully-turbulent conditions over rectangular cavity geometries that are aligned with the local velocity. Turbulent boundary layer properties were computationally determined and used to develop the cavity geometry parametrics and to correlate experimental closed cavity heating data to yield new relationships for the floor-averaged and centerline endwall peak-heating augmentation. With the form of the closed-cavity correlations established, historical data were used to develop new correlations for turbulent open-cavity heating.

Everhart, Joel L.↗

Turbulent Supersonic/Hypersonic Heating Correlations for Open and Closed Cavities

Supersonic/hypersonic laminar heating correlations that were developed for damage assessment analysis of atmospheric re-entry vehicles have been modified and extended to cover fully-turbulent conditions over rectangular cavity geometries that are aligned with the local velocity. Turbulent boundary layer properties were computationally determined and used to develop the cavity geometry parametrics and to correlate experimental closed cavity heating data to yield new relationships for the floor-averaged and centerline endwall peakheating augmentation. With the form of the closed-cavity correlations established, historical data were used to develop new correlations for turbulent open-cavity heating.

Everhart, Joel L.↗

Numerical Predictions of Mode Reflections in an Open Circular Duct: Comparison with Theory

The NASA Broadband Aeroacoustic Stator Simulation code was used to compute the acoustic field for higher-order modes in a circular duct geometry. To test the accuracy of the results computed by the code, the duct was terminated by an open end with an infinite flange or no flange. Both open end conditions have a theoretical solution that was used to compare with the computed results. Excellent comparison for reflection matrix values was achieved after suitable refinement of the grid at the open end. The study also revealed issues with the level of the mode amplitude introduced into the acoustic held from the source boundary and the amount of reflection that occurred at the source boundary when a general nonreflecting boundary condition was applied.

Duct Modes↗