Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “surface codes”

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 91 records · Page 5

Recent Developments in Grid Generation and Force Integration Technology for Overset Grids

Recent developments in algorithms and software tools for generating overset grids for complex configurations are described. These include the overset surface grid generation code SURGRD and version 2.0 of the hyperbolic volume grid generation code HYPGEN. The SURGRD code is in beta test mode where the new features include the capability to march over a collection of panel networks, a variety of ways to control the side boundaries and the marching step sizes and distance, a more robust projection scheme and an interpolation option. New features in version 2.0 of HYPGEN include a wider range of boundary condition types. The code also allows the user to specify different marching step sizes and distance for each point on the surface grid. A scheme that takes into account of the overlapped zones on the body surface for the purpose of forces and moments computation is also briefly described, The process involves the following two software modules: MIXSUR - a composite grid generation module to produce a collection of quadrilaterals and triangles on which pressure and viscous stresses are to be integrated, and OVERINT - a forces and moments integration module.

Chan, William M.↗

The Project for Intercomparison of Land-surface Parameterization Schemes

The Project for Intercomparison of Land-surface Parameterization Schemes (PILPS) is described and the first stage science plan outlined. PILPS is a project designed to improve the parameterization of the continental surface, especially the hydrological, energy, momentum, and carbon exchanges with the atmosphere. The PILPS Science Plan incorporates enhanced documentation, comparison, and validation of continental surface parameterization schemes by community participation. Potential participants include code developers, code users, and those who can provide datasets for validation and who have expertise of value in this exercise. PILPS is an important activity because existing intercomparisons, although piecemeal, demonstrate that there are significant differences in the formulation of individual processes in the available land surface schemes. These differences are comparable to other recognized differences among current global climate models such as cloud and convection parameterizations. It is also clear that too few sensitivity studies have been undertaken with the result that there is not yet enough information to indicate which simplifications or omissions are important for the near-surface continental climate, hydrology, and biogeochemistry. PILPS emphasizes sensitivity studies with and intercomparisons of existing land surface codes and the development of areally extensive datasets for their testing and validation.

Henderson-Sellers, A.↗

Coupling Carbon Oxidation and Surface Recession in Direct-Simulation Monte Carlo Code, SPARTA

Ablative thermal protection system (TPS) materials for spacecraft are composites that are often made out of carbon-based reinforcement and a polymeric matrix. They endure high-temperature oxidation and surface recession when re-entering Earth’s atmosphere. Ablation is the result of many coupled and competing thermal, mechanical, and chemical phenomena, and it is difficult to isolate the role of each on the overall degradation of the TPS. Here we develop an ablation model for material recession coupled explicitly to finite rate carbon oxidation in complex microstructures. In this work, Stochastic PArallel Rarified-gas Time-accurate Analyzer (SPARTA), a direct-simulation Monte Carlo (DSMC) code, is modified to allow oxidation-driven ablation of implicitly defined carbon surfaces. In SPARTA, implicit surfaces are generated from the grid corner point values via a marching cubes algorithm, therefore creating a new set of surface elements every time ablation is performed. The finite-rate oxidation model developed by Gopalan et. al, was adapted to tally surface reactions and other surface data on a per-grid cell basis. The ablation functionality was also adjusted so once the reactions have occurred, the number of reactions leading to CO formation can be converted to corner point reduction values; therefore, carbon removal is directly proportional to surface recession. We also develop robust algorithms which handle the evolution of the flow cells and solid material regions, including split cells (flow cell divided in two by a solid surface). Finally, we demonstrate our implicit chemistry model for 2D and 3D geometries by producing reaction statistics and detailed visualization of oxidation-induced material recession at the microscale.

V Arias↗

Automatic Conversion of Conceptual Geometry to CFD Geometry for Aircraft Design

Conceptual aircraft design is usually based on simple analysis codes. Its objective is to provide an overall system performance of the developed concept, while preliminary aircraft design uses high-fidelity analysis tools such as computational fluid dynamics (CFD) analysis codes or finite element structural analysis codes. In some applications, such as low-boom supersonic concept development, it is important to be able to explore a variety of drastically different configurations while using CFD analysis to check whether a given configuration can be tailored to have a low-boom ground signature. It poses an extremely challenging problem of integrating CFD analysis in conceptual design. This presentation will discuss a computer code, called iPatch, for automatic conversion of conceptual geometry to CFD geometry. In general, conceptual aircraft geometry is not as well-defined as a CAD geometry model. In particular, a conceptual aircraft geometry model usually does not define the intersection curves for the connecting surfaces. The computer code iPatch eliminates the gap between conceptual geometry and CFD geometry by accomplishing the following three tasks automatically: (1) use bicubic B-splines to extrapolate (if necessary) each surface in a conceptual geometry so that all the independently defined geometry components (such as wing and fuselage) can be intersected to form a watertight CFD geometry, (2) compute the intersection curves of surface patches at any resolution (up to 10-7 accuracy) specified by users, and (3) write the B-spline surface patches and the corresponding boundary points for the watertight CFD geometry in the format that can be directly exported to the meshing tool VGRID in the CFD software TetrUSS. As a result, conceptual designers can get quick feedback on the aerodynamic characteristics of their concepts, which will allow them to understand some subtlety in their concepts and to be able to assess their concepts with a higher degree of confidence. This integration of CFD analysis in conceptual aircraft design will greatly eliminate some uncertainty due to simple analysis codes used to develop the concepts and improve the feasibility/credibility of the final concept. The presentation will highlight the mathematical challenges of accomplishing the aforementioned three tasks and the computational algorithms used by iPatch.

Li, Wu↗

Intelligent Patching of Conceptual Geometry for CFD Analysis

The iPatch computer code for intelligently patching surface grids was developed to convert conceptual geometry to computational fluid dynamics (CFD) geometry (see figure). It automatically uses bicubic B-splines to extrapolate (if necessary) each surface in a conceptual geometry so that all the independently defined geometric components (such as wing and fuselage) can be intersected to form a watertight CFD geometry. The software also computes the intersection curves of surface patches at any resolution (up to 10.4 accuracy) specified by the user, and it writes the B-spline surface patches, and the corresponding boundary points, for the watertight CFD geometry in the format that can be directly used by the grid generation tool VGRID. iPatch requires that input geometry be in PLOT3D format where each component surface is defined by a rectangular grid {(x(i,j), y(i,j), z(i,j)):1less than or equal to i less than or equal to m, 1 less than or equal to j less than or equal to n} that represents a smooth B-spline surface. All surfaces in the PLOT3D file conceptually represent a watertight geometry of components of an aircraft on the half-space y greater than or equal to 0. Overlapping surfaces are not allowed, but could be fixed by a utility code "fixp3d". The fixp3d utility code first finds the two grid lines on the two surface grids that are closest to each other in Hausdorff distance (a metric to measure the discrepancies of two sets); then uses one of the grid lines as the transition line, extending grid lines on one grid to the other grid to form a merged grid. Any two connecting surfaces shall have a "visually" common boundary curve, or can be described by an intersection relationship defined in a geometry specification file. The intersection of two surfaces can be at a conceptual level. However, the intersection is directional (along either i or j index direction), and each intersecting grid line (or its spine extrapolation) on the first surface should intersect the second surface. No two intersection relationships will result in a common intersection point of three surfaces. The output files of iPatch are IGES, d3m, and mapbc files that define the CFD geometry in VGRID format. The IGES file gives the NURBS definition of the outer mold line in the geometry. The d3m file defines how the outer mold line is broken into surface patches whose boundary curves are defined by points. The mapbc file specifies what the boundary condition is on each patch and the corresponding NURBS surface definition of each non-planar patch in the IGES file.

Li, Wu↗

Chemical nonequilibrium Navier-Stokes solutions for hypersonic flow over an ablating graphite nosetip

The general boundary conditions including mass and energy balances of chemically equilibrated or nonequilibrated gas adjacent to ablating surfaces have been derived. A computer procedure based on these conditions was developed and interfaced with the Navier-Stokes solver for predictions of the flow field, surface temperature, and surface ablation rates over re-entry space vehicles with ablating Thermal Protection Systems (TPS). The Navier-Stokes solver with general surface thermochemistry boundary conditions can predict more realistic solutions and provide useful information for the design of TPS. A test case with a proposed hypersonic test vehicle configuration and associated free stream conditions was developed. Solutions with various surface boundary conditions were obtained, and the effect of nonequilibrium gas as well as surface chemistry on surface heating and ablation rate were examined. The solutions of the GASP code with complete ablating surface conditions were compared with those of the ASC code. The direction of future work is also discussed.

Chen, Y. K.↗

Computer code for scattering from impedance bodies of revolution. Part 3: Surface impedance with s and phi variation. Analytical and numerical results

The third phase of the development of the computer codes for scattering by coated bodies that has been part of an ongoing effort in the Electromagnetics Laboratory of the Electrical Engineering and Computer Science Department at the University of Illinois at Chicago is described. The work reported discusses the analytical and numerical results for the scattering of an obliquely incident plane wave by impedance bodies of revolution with phi variation of the surface impedance. Integral equation formulation of the problem is considered. All three types of integral equations, electric field, magnetic field, and combined field, are considered. These equations are solved numerically via the method of moments with parametric elements. Both TE and TM polarization of the incident plane wave are considered. The surface impedance is allowed to vary along both the profile of the scatterer and in the phi direction. Computer code developed for this purpose determines the electric surface current as well as the bistatic radar cross section. The results obtained with this code were validated by comparing the results with available results for specific scatterers such as the perfectly conducting sphere. Results for the cone-sphere and cone-cylinder-sphere for the case of an axially incident plane were validated by comparing the results with the results with those obtained in the first phase of this project. Results for body of revolution scatterers with an abrupt change in the surface impedance along both the profile of the scatterer and the phi direction are presented.

Uslenghi, Piergiorgio L. E.↗

Code Calibration Applied to the TCA High-Lift Model in the 14 x 22 Wind Tunnel (Simulation With and Without Model Post-Mount)

The objective of this study is to calibrate a Navier-Stokes code for the TCA (30/10) baseline configuration (partial span leading edge flaps were deflected at 30 degs. and all the trailing edge flaps were deflected at 10 degs). The computational results for several angles of attack are compared with experimental force, moments, and surface pressures. The code used in this study is CFL3D; mesh sequencing and multi-grid were used to full advantage to accelerate convergence. A multi-grid approach was used similar to that used for the Reference H configuration allowing point-to-point matching across all the trailingedge block interfaces. From past experiences with the Reference H (ie, good force, moment, and pressure comparisons were obtained), it was assumed that the mounting system would produce small effects; hence, it was not initially modeled. However, comparisons of lower surface pressures indicated the post mount significantly influenced the lower surface pressures, so the post geometry was inserted into the existing grid using Chimera (overset grids).

Lessard, Wendy B.↗

Numerical Simulations of Cosmogenic Neutron Production and Transport in Planetary Surfaces

The numerical simulation code MCNPX was used to calculate the production and transport of cosmic-ray-produced neutrons in Mars and meteorites. These calculations help to understand the processes involved and the parameters that control the neutron fluxes. Results are presented here for neutrons in Mars and for the distribution of cosmic-ray neutrons and protons in a 50 cm-radius L-chondrite.

Kim, Kyeong J.↗

Space Debris Surfaces - Probability of no penetration versus impact velocity and obliquity

A collection of computer codes called Space Debris Surfaces (SD-SURF), have been developed to assist in the design and analysis of space debris protection systems. An SD-SURF analysis will show which obliquities and velocities are most likely to cause a penetration to help the analyst select a shield design best suited to the predominant penetration mechanism. Examples of the interaction between space vehicle geometry, the space debris environment, and the penetration and critical damage ballistic limit surfaces of the shield under consideration are presented.

Elfer, N.↗

Hypersonic code efficiency and validation studies

Renewed interest in hypersonic and supersonic flows spurred the development of the Compressible Navier-Stokes (CNS) code. Originally developed for external flows, CNS was modified to enable it to also be applied to internal high speed flows. In the initial phase of this study CNS was applied to both internal flow applications and fellow researchers were taught to run CNS. The second phase of this research was the development of surface grids over various aircraft configurations for the High Speed Research Program (HSRP). The complex nature of these configurations required the development of improved surface grid generation techniques. A significant portion of the grid generation effort was devoted to testing and recommending modifications to early versions of the S3D surface grid generation code.

Bennett, Bradford C.↗

Numerical solution of supersonic viscous flow over blunt delta wings

A general parabolized Navier-Stokes code has been developed to compute the steady supersonic viscous flow around arbitrary body shapes at high angles of attack. A nonorthogonal three-dimensional coordinate frame permits the code to march with solution surfaces which are the most appropriate. The code has been used to calculate the laminar flow over a slab delta wing with 70 deg sweep at angles of attack up to 41.5 deg and Mach numbers of 6.8 and 9.6. The computed shock shapes, surface pressures and heat transfer coefficients are compared with experiment and show agreement.

Tannehill, J. C.↗

Comparison of CFD Predictions with Shuttle Global Flight Thermal Imagery and Discrete Surface Measurements

Surface temperature measurements from the STS-119 boundary-layer transition experiment on the space shuttle orbiter Discovery provide a rare opportunity to assess turbulent CFD models at hypersonic flight conditions. This flight data was acquired by on-board thermocouples and by infrared images taken off-board by the Hypersonic Thermodynamic Infrared Measurements (HYTHIRM) team, and is suitable for hypersonic CFD turbulence assessment between Mach 6 and 14. The primary assessment is for the Baldwin-Lomax and Cebeci-Smith algebraic turbulence models in the DPLR and LAURA CFD codes, respectively. A secondary assessment is made of the Shear-Stress Transport (SST) two-equation turbulence model in the DPLR code. Based upon surface temperature comparisons at eleven thermocouple locations, the algebraic-model turbulent CFD results average 4% lower than the measurements for Mach numbers less than 11. For Mach numbers greater than 11, the algebraic-model turbulent CFD results average 5% higher than the three available thermocouple measurements. Surface temperature predictions from the two SST cases were consistently 3 4% higher than the algebraic-model results. The thermocouple temperatures exhibit a change in trend with Mach number at about Mach 11; this trend is not reflected in the CFD results. Because the temperature trends from the turbulent CFD simulations and the flight data diverge above Mach 11, extrapolation of the turbulent CFD accuracy to higher Mach numbers is not recommended.

Wood, William A.↗

Transonic Drag Prediction on a DLR-F6 Transport Configuration Using Unstructured Grid Solvers

A second international AIAA Drag Prediction Workshop (DPW-II) was organized and held in Orlando Florida on June 21-22, 2003. The primary purpose was to inves- tigate the code-to-code uncertainty. address the sensitivity of the drag prediction to grid size and quantify the uncertainty in predicting nacelle/pylon drag increments at a transonic cruise condition. This paper presents an in-depth analysis of the DPW-II computational results from three state-of-the-art unstructured grid Navier-Stokes flow solvers exercised on similar families of tetrahedral grids. The flow solvers are USM3D - a tetrahedral cell-centered upwind solver. FUN3D - a tetrahedral node-centered upwind solver, and NSU3D - a general element node-centered central-differenced solver. For the wingbody, the total drag predicted for a constant-lift transonic cruise condition showed a decrease in code-to-code variation with grid refinement as expected. For the same flight condition, the wing/body/nacelle/pylon total drag and the nacelle/pylon drag increment predicted showed an increase in code-to-code variation with grid refinement. Although the range in total drag for the wingbody fine grids was only 5 counts, a code-to-code comparison of surface pressures and surface restricted streamlines indicated that the three solvers were not all converging to the same flow solutions- different shock locations and separation patterns were evident. Similarly, the wing/body/nacelle/pylon solutions did not appear to be converging to the same flow solutions. Overall, grid refinement did not consistently improve the correlation with experimental data for either the wingbody or the wing/body/nacelle pylon configuration. Although the absolute values of total drag predicted by two of the solvers for the medium and fine grids did not compare well with the experiment, the incremental drag predictions were within plus or minus 3 counts of the experimental data. The correlation with experimental incremental drag was not significantly changed by specifying transition. Although the sources of code-to-code variation in force and moment predictions for the three unstructured grid codes have not yet been identified, the current study reinforces the necessity of applying multiple codes to the same application to assess uncertainty.

Lee-Rausch, E. M.↗

Parallelized Carbon Oxidation and Surface Recession Model in Direct-Simulation Monte Carlo Code, SPARTA

Ablative thermal protection system (TPS) materials for spacecraft are composites that often consist of a carbon-based reinforcement and a polymeric matrix. During Earth re-entry, they endure high-temperature oxidation and surface recession. Oxidation is an important mechanism for ablation, sometimes leading to the weakening, spallation, or failure of the oxidized fibers at the surface and in the char layer. However, more details are required including accurate material properties of the fiber microstructure, whether the fibers recede homogenously or localized at pits, and the role of pyrolysis outgassing in order to evaluate the role oxidation plays in the degradation and failure mechanisms of these materials. In this work, we demonstrate a parallelizable oxidation-driven ablation model developed for detailed, large-scale simulations in the DSMC code SPARTA. We also develop robust algorithms which handle the conservation of the surface state after an ablation step. Finally, we verify our model for both simple and more complex chemistry as well as microstructures with reaction statistics, oxidation depth calculations, and detailed visualization of oxidation-induced material recession at the microscale.

V Arias↗

Land Surface Temperature Measurements from EOD MODIS Data

We made more tests of the version 2.0 daily Level 2 and Level 3 Land-Surface Temperature (LST) code (PGE 16) jointly with the MODIS Science Data Support Team (SDST). After making minor changes a few times, the PGE16 code has been successfully integrated and tested by MODIS SDST, and recently has passed the inspection at the Goddard Distributed Active Archive Center (DAAC). We conducted a field campaign in the area of Mono Lake, California on March 10, 1998, in order to validate the MODIS LST algorithm in cold and dry conditions. Two MODIS Airborne Simulator (MAS) flights were completed during the field campaign, one before noon, and another around 10 pm PST. The weather condition for the daytime flight was perfect: clear sky, the column water vapor measured by radiosonde around 0.3 cm, and wind speed less than a half meter per second. The quality of MAS data is good for both day and night flights. We analyzed the noise equivalent temperature difference (NE(delta)T) and the calibration accuracy of the seven MAS thermal infrared (TIR) bands, that are used in the MODIS day/night LST algorithm, with daytime MAS data over four flat homogeneous study areas: two on Grant Lake (covered with ice and snow, respectively), one on Mono Lake, and another on the snow field site where we made field measurements. NE(delta)T ranges from 0.2 to 0.6 k for bands 42, 45, 46, and 48. It ranges from 0.8 to 1.1 K for bands 30-32. The day and night MAS data have been used to retrieve surface temperature and emissivities in these bands. A simple method to correct the effect of night thin cirrus has been incorporated into the day/night LST algorithm in dry atmospheric conditions. We compared the retrieved surface temperatures with those measured with TIR spectrometer, radiometers and thermistors in the snow test site, and the retrieved emissivity images with topographic map. The daytime LST values match well within 1 K. The night LST retrieved from MAS data is 3.3 K colder than those from field measurements most likely because of the effect of haze at night. The good agreement among the regional averaged surface temperatures obtained from LST values retrieved at different resolutions increased our confidence in the MODIS day/night LST algorithm.

Wan, Zheng-Ming↗

Three-dimensional canard-wing shape optimization in aircraft cruise and maneuver environments

This paper demonstrates a numerical technique for canard-wing shape optimization at two operating conditions. For purposes of simplicity, a mean surface wing paneling code is employed for the aerodynamic calculations. The optimization procedures are based on the method of feasible directions. The shape functions for describing the thickness, camber, and twist are based on polynomial representations. The primary design requirements imposed restrictions on the canard and wing volumes and on the lift coefficients at the operating conditions. Results indicate that significant improvements in minimum drag and lift-to-drag ratio are possible with reasonable aircraft geometries. Calculations were done for supersonic speeds with Mach numbers ranging from 1 to 6. Planforms were mainly of a delta shape with aspect ratio of 1.

De Silva, B. M. E.↗

Summary and recent results from the NASA advanced High Speed Propeller Research Program

Advanced high-speed propellers offer large performance improvements for aircraft that cruise in the Mach 0.7 to 0.8 speed regime. The current status of the NASA research program on high-speed propeller aerodynamics, acoustics, and aeroelastics is described. Recent wind tunnel results for five 8- to 10-blade advanced models are compared with analytical predictions. Test results show that blade sweep was important in achieving net efficiencies near 80 percent at Mach 0.8 and reducing near-field cruise noise by dB. Lifting line and lifting surface aerodynamic analysis codes are under development and some initial lifting line results are compared with propeller force and probe data. Some initial laser velocimeter measurements of the flow field velocities of an 8-bladed 45 deg swept propeller are shown. Experimental aeroelastic results indicate that cascade effects and blade sweep strongly affect propeller aeroelastic characteristics. Comparisons of propeller near-field noise data with linear acoustic theory indicate that the theory adequate predicts near-field noise for subsonic tip speeds but overpredicts the noise for supersonic tip speeds. Potential large gains in propeller efficiency of 7 to 11 percent at Mach 0.8 may be possible with advanced counter-rotation propellers.

Mitchell, G. A.↗