SEARCH · Engineering Papers
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.
A pair of flat grid blanks are clamped together at their edges with an impervious metal sheet on top. All of the blanks and sheets are dished simultaneously by forcing fluid to inflate an elastic sheet which contacts the bottom grid blank. A second impervious metal sheet is inserted between the two grid blanks if the grids have high percentage open areas. The dished grids are stress relieved simultaneously.
Two endurance tests of a 5-cm mercury bombardment thruster are reported. Both tests used a translational screen-grid system with the beam vectored 10 degrees. The first test lasted 141 hours and the second test operated for 2026 hours. In each test two fused silica samples (solar cell covers), 2.0 cm by 2.1 cm, were placed in shielded holders to detect materials sputtered from the thruster. Spectral optical properties between 0.398 and 2.16 microns were measured on each sample, both before and after the endurance tests. The deposition on each sample was spectrographically analyzed to determine the type of materials sputtered from the thruster. It was found that sputtering from the neutralizer is highly dependent on its position with respect to the beam edge. The sputtering from the accelerator grid of the translational screen-grid system of the 2026 hour test was sufficient to form an opaque film on the sample located in the direction opposite to the vectored beam.
NASA Lewis is currently engaged in a research effort as a team member of the High Alpha Technology Program (HATP) within NASA. This program utilizes a specially equipped F/A-18, the High Alpha Research Vehicle (HARV), in an ambitious effort to improve the maneuverability of high-performance military aircraft at low subsonic speed, high angle of attack conditions. The overall objective of the Lewis effort is to develop inlet technology that will ensure efficient airflow delivery to the engine during these maneuvers. One part of the Lewis approach utilizes computational fluid dynamics codes to predict the installed performance of inlets for these highly maneuverable aircraft. Full Navier-Stokes (FNS) calculations on the installed F/A-18 inlet at 30 degrees angle of attack, 0 degrees yaw, and a freestream Mach number of 0.2 have been obtained in this study using an algebraic turbulence model with two grids (original and revised). Results obtained with the original grid were used to determine where further grid refinements and additional geometry were needed. In order to account properly for the external effects, the forebody, leading edge extension (LEX), ramp, and wing were included with inlet geometry. In the original grid, the diverter, LEX slot, and leading edge flap were not included due to insufficient geometry definition, but were included in a revised grid. In addition, a thin-layer Navier-Stokes (TLNS) code is used with the revised grid and the numerical results are compared to those obtained with the FNS code. The TLNS code was used to evaluate the effects on the solution using a code with more recent CFD developments such as upwinding with TVD schemes versus central differencing with artificial dissipation. The calculations are compared to a limited amount of available experimental data. The predicted forebody/fuselage surface static pressures compared well with data of all solutions. The predicted trajectory of the vortex generated under the LEX was different for each solution. These discrepancies are attributed to differences in the grid resolution and turbulence modeling. All solutions predict that this vortex is ingested by the inlet. The predicted inlet total pressure recoveries are lower than data and the distortions are higher than data. The results obtained with the revised grid were significantly improved from the original grid results. The original grid results indicated the ingested vortex migrated to the engine face and caused additional distortions to those already present due to secondary flow development. The revised grid results indicate that the ingested vortex is dissipated along the inlet duct inboard wall. The TLNS results indicate the flow at the engine face was much more distorted than the FNS results and is attributed to the pole boundary condition introducing numerical distortions into the flow field.
An elliptic grid generation scheme is described which produces a curvilinear, boundary-fitted grid for highly swept wings with sharp leading edges. The three-dimensional integration domain is segmented in successive crossflow sections wherein the physical coordinates are transformed to computational coordinates by solving an elliptic set of two quasi-linear differential equations. Grid points are clustered in the vicinity of the wing surface, especially in the neighborhood of the leading edges. The formulation of the source terms in the equations governing the transformation contains adjustable parameters which are determined from limiting forms of the elliptic equations along the boundaries. The source terms are assigned to each grid point by interpolation of these parameters into the interior domain. The formulation does not depend on the boundary shape or on the distribution of the grid points along the boundaries. For the first time, the discretized governing equations have been solved using a fast AF1 iteration scheme.
This study investigates the effects of grid topology and grid adaptation on numerical solutions of the Navier-Stokes equations. In the first part of this study, a general procedure is presented for computation of high-speed flow over complex three-dimensional configurations. The flow field is simulated on the surface of a Butler wing in a uniform stream. Results are presented for Mach number 3.5 and a Reynolds number of 2,000,000. The O-type and H-type grids have been used for this study, and the results are compared together and with other theoretical and experimental results. The results demonstrate that while the H-type grid is suitable for the leading and trailing edges, a more accurate solution can be obtained for the middle part of the wing with an O-type grid. In the second part of this study, methods of grid adaption are reviewed and a method is developed with the capability of adapting to several variables. This method is based on a variational approach and is an algebraic method. Also, the method has been formulated in such a way that there is no need for any matrix inversion. This method is used in conjunction with the calculation of hypersonic flow over a blunt-nose body. A movie has been produced which shows simultaneously the transient behavior of the solution and the grid adaption.
This paper describes the spectral analysis of integrated pressures on patches of the scale model of the Space Launch System (SLS) Block 1B crew vehicle with the Unsteady Pressure-Sensitive Paint (uPSP) measurements, which were collected in the ascent transient aerodynamics tests with the Unitary Plan Wind Tunnel 11-by-11-foot Transonic Wind Tunnel in September 2019 at NASA Ames Research Center. Recent research has demonstrated that uPSP can be an essential tool for the assessment of the unsteady, aerodynamic phenomena. The work described in this paper is a part of NASA’s development of a new state-of-the-art uPSP capability in production wind tunnels. In this paper, 108 patches are defined with x station values and azimuth angles of the scale model. For each patch, the polygons are determined from the surface cells of the grid of the model, clipped with the edges of the patch, and each of the polygons is divided into triangles. The inputs of the pressure integration algorithm are the time series of pressure coefficients on the vertices of the grid of the model, generated by the uPSP data processing program from the videos taken with 4 Phantom high-speed cameras during the wind tunnel tests. The integrated pressure of the patch is determined as the ratio of the sum of the forces on the triangles over the sum of the areas of the triangles. For each of the test cases, the Cross Power Spectral Density (CPSD) and Magnitude-Squared Coherence (msCohere) are computed for frequencies over 1/3 octave bands from the time series of the integrated pressures on the patches. The spectral analysis outputs of different cases of the wind tunnel tests are analyzed and the coherence of patches is investigated over frequencies, x station values, azimuth angles, Mach numbers and SLS vehicle attitudes. The pressure integration and spectral analysis described in this paper were executed on the NASA Pleiades supercomputer. Funding for this research was provided by the NASA Aeroscience Evaluation and Test Capabilities (AETC) Project.
Three-dimensional thin-layer Navier-Stokes computations are presented for the F/A-18 configuration. The modeled configuration includes an accurate surface representation of the fuselage, leading-edge extension (LEX), and wing, both with and without leading-edge flap deflection. A multiblock structured grid strategy is employed to decompose the computational flowfield domain around the subject configuration. Steady-state solutions are obtained from an algorithm that solves the compressible Navier-Stokes equations with an upwind-biased, flux-difference splitting approach. The results presented are based on a fully turbulent flow assumption, simulating the high Reynolds number flow conditions that correspond to a recent FI A-18 flight experiment. Good agreements between the computations and the flight test results are obtained for both surface flow patterns as well as surface pressure distributions. Furthermore, a correlation between the computed LEX vortex-core and the flight test results, observed by way of smoke visualization, is also presented.
Turbulent thin-layer, Reynolds-Averaged Navier-Stokes solutions, based on a multi-block structured grid, are presented for a 65 deg delta wing having either a sharp leading edge (SLE) or blunt leading edge (BLE) geometry. The primary objective of the study is to assess the prediction capability of the method for simulating the leading-edge flow separation and the ensuing vortex flow characteristics. Computational results are obtained for two angles of attack of approximately 13 and 20 deg, at free-stream Mach number of 0.40 and Reynolds number of 6 million based on the wing mean aerodynamic chord. The effects of two turbulence models of Baldwin-Lomax with Degani-Schiff (BL/DS) and the Spalart-Allmaras (SA) on the numerical results are also discussed. The computations also explore the effects of two numerical flux-splitting schemes, i.e., flux difference splitting (fds) and flux vector splitting (fvs), on the solution development and convergence characteristics. The resulting trends in solution sensitivity to grid resolution for the selected leading-edge geometries, angles of attack, turbulence models and flux splitting schemes are also presented. The validity of the numerical results is evaluated against a unique set of experimental wind-tunnel data that was obtained in the National Transonic Facility at the NASA Langley Research Center.
A composite grid was generated for flows in turbines. It consisted of the C-grid (or O-grid) in the immediate vicinity of the blade and the H-grid in the middle of the blade passage between the C-grids and in the upstream region. This new composite grid provides better smoothness, resolution, and orthogonality than any single grid for a typical turbine blade with a large camber and rounded leading and trailing edges. The C-H (or O-H) composite grid has an unusual grid point that is connected to more than four neighboring nodes in two dimensions (more than six neighboring nodes in three dimensions). A finite-volume lower-upper (LU) implicit scheme to be used on this grid poses no problem and requires no special treatment because each interior cell of this composite grid has only four neighboring cells in two dimensions (six cells in three dimensions). The LU implicit scheme was demonstrated to be efficient and robust for external flows in a broad flow regime and can be easily applied to internal flows and extended from two to three dimensions.
The NASA Short-term Prediction Research and Transition (SPoRT) Center has been running a real-time version of the Land Information System (LIS) since summer 2010 (hereafter, SPoRTLIS). The real-time SPoRT-LIS runs the Noah land surface model (LSM) in an offline capacity apart from a numerical weather prediction model, using input atmospheric and precipitation analyses (i.e., "forcings") to drive the Noah LSM integration at 3-km resolution. Its objectives are to (1) produce local-scale information about the soil state for NOAA/National Weather Service (NWS) situational awareness applications such as drought monitoring and assessing flood potential, and (2) provide land surface initialization fields for local modeling initiatives. The current domain extent has been limited by the input atmospheric analyses that drive the Noah LSM integration within SPoRT-LIS, specifically the National Centers for Environmental Prediction (NCEP) Stage IV precipitation analyses. Due to the nature of the geographical edges of the Stage IV precipitation grid and its limitations in the western U.S., the SPoRT-LIS was originally confined to a domain fully nested within the Stage IV grid, over the southeastern half of the Conterminous United States (CONUS). In order to expand the real-time SPoRT-LIS to a full CONUS domain, alternative precipitation forcing datasets were explored in year-long, offline comparison runs of the Noah LSM. Based on results of these comparison simulations, we chose to implement the radar/gauge-based precipitation analyses from the National Severe Storms Laboratory as a replacement to the Stage IV product. The Multi-Radar Multi-Sensor (MRMS; formerly known as the National Mosaic and multi-sensor Quantitative precipitation estimate) product has full CONUS coverage at higher-resolution, thereby providing better coverage and greater detail than that of the Stage IV product. This paper will describe the expanded/upgraded SPoRT-LIS, present comparisons between the original and upgraded SPoRT-LIS, and discuss the path forward for future collaboration opportunities with SPoRT partners in the NWS.
An Observing System Simulation Experiment (OSSE) is a pure modeling study used when actual observations are too expensive or difficult to obtain. OSSEs are valuable tools for determining the potential impact of new observing systems on numerical weather forecasts and for evaluation of data assimilation systems (DAS). An OSSE has been developed at the NASA Global Modeling and Assimilation Office (GMAO, Errico et al 2013). The GMAO OSSE uses a 13-month integration of the European Centre for Medium- Range Weather Forecasts 2005 operational model at T511/L91 resolution for the Nature Run (NR). Synthetic observations have been updated so that they are based on real observations during the summer of 2013. The emulated observation types include AMSU-A, MHS, IASI, AIRS, and HIRS4 radiance data, GPS-RO, and conventional types including aircraft, rawinsonde, profiler, surface, and satellite winds. The synthetic satellite wind observations are colocated with the NR cloud fields, and the rawinsondes are advected during ascent using the NR wind fields. Data counts for the synthetic observations are matched as closely as possible to real data counts, as shown in Figure 2. Errors are added to the synthetic observations to emulate representativeness and instrument errors. The synthetic errors are calibrated so that the statistics of observation innovation and analysis increments in the OSSE are similar to the same statistics for assimilation of real observations, in an iterative method described by Errico et al (2013). The standard deviations of observation minus forecast (xo-H(xb)) are compared for the OSSE and real data in Figure 3. The synthetic errors include both random, uncorrelated errors, and an additional correlated error component for some observational types. Vertically correlated errors are included for conventional sounding data and GPS-RO, and channel correlated errors are introduced to AIRS and IASI (Figure 4). HIRS, AMSU-A, and MHS have a component of horizontally correlated error. The forecast model used by the GMAO OSSE is the Goddard Earth Observing System Model, Version 5 (GEOS-5) with Gridpoint Statistical Interpolation (GSI) DAS. The model version has been updated to v. 5.13.3, corresponding to the current operational model. Forecasts are run on a cube-sphere grid with 180 points along each edge of the cube (approximately 0.5 degree horizontal resolution) with 72 vertical levels. The DAS is cycled at 6-hour intervals, with 240 hour forecasts launched daily at 0000 UTC. Evaluation of the forecasting skill for July and August is currently underway. Prior versions of the GMAO OSSE have been found to have greater forecasting skill than real world forecasts. It is anticipated that similar forecast skill will be found in the updated OSSE.
A computational fluid dynamics (CFD) method has been employed to compute vortical flows around slender wing/body configurations. The emphasis of the paper is on the effectiveness of an adaptive grid procedure in "capturing" concentrated vortices generated at sharp edges or flow separation lines of lifting surfaces flying at high angles of attack. The method is based on a tetrahedral unstructured grid technology developed at the NASA Langley Research Center. Two steady-state, subsonic, inviscid and Navier-Stokes flow test cases are presented to demonstrate the applicability of the method for solving vortical flow problems. The first test case concerns vortex flow over a simple 65 delta wing with different values of leading-edge radius. Although the geometry is quite simple, it poses a challenging problem for computing vortices originating from blunt leading edges. The second case is that of a more complex fighter configuration. The superiority of the adapted solutions in capturing the vortex flow structure over the conventional unadapted results is demonstrated by comparisons with the wind-tunnel experimental data. The study shows that numerical prediction of vortical flows is highly sensitive to the local grid resolution and that the implementation of grid adaptation is essential when applying CFD methods to such complicated flow problems.
A computational fluid dynamics (CFD) method has been employed to compute vortical flows around slender wing/body configurations. The emphasis of the paper is on the effectiveness of an adaptive grid procedure in "capturing" concentrated vortices generated at sharp edges or flow separation lines of lifting surfaces flying at high angles of attack. The method is based on a tetrahedral unstructured grid technology developed at the NASA Langley Research Center. Two steady-state, subsonic, inviscid and Navier-Stokes flow test cases are presented to demonstrate the applicability of the method for solving practical vortical flow problems. The first test case concerns vortex flow over a simple 65' delta wing with different values of leading-edge bluntness, and the second case is that of a more complex fighter configuration. The superiority of the adapted solutions in capturing the vortex flow structure over the conventional unadapted results is demonstrated by comparisons with the windtunnel experimental data. The study shows that numerical prediction of vortical flows is highly sensitive to the local grid resolution and that the implementation of grid adaptation is essential when applying CFD methods to such complicated flow problems.
A computational fluid dynamics (CFD) method has been employed to compute vortical flows around slender wing/body configurations. The emphasis of the paper is on the effectiveness of an adaptive grid procedure in "capturing" concentrated vortices generated at sharp edges or flow separation lines of lifting surfaces flying at high angles of attack. The method is based on a tetrahedral unstructured grid technology developed at the NASA Langley Research Center. Two steady-state, subsonic, inviscid and Navier-Stokes flow test cases are presented to demonstrate the applicability of the method for solving practical vortical flow problems. The first test case concerns vortex flow over a simple 65deg delta wing with different values of leading-edge bluntness, and the second case is that of a more complex fighter configuration. The superiority of the adapted solutions in capturing the vortex flow structure over the conventional unadapted results is demonstrated by comparisons with the windtunnel experimental data. The study shows that numerical prediction of vortical flows is highly sensitive to the local grid resolution and that the implementation of grid adaptation is essential when applying CFD methods to such complicated flow problems.
A computational fluid dynamics (CFD) method has been employed to compute vortical flows around slender wing/body configurations. The emphasis of the paper is on the effectiveness of an adaptive grid procedure in "capturing" concentrated vortices generated at sharp edges or flow separation lines of lifting surfaces flying at high angles of attack. The method is based on a tetrahedral unstructured grid technology developed at the NASA Langley Research Center. Two steady-state, subsonic, inviscid and Navier-Stokes flow test cases are presented to demonstrate the applicability of the method for solving practical vortical flow problems. The first test case concerns vortex flow over a simple 65 delta wing with different values of leading-edge radius. Although the geometry is quite simple, it poses a challenging problem for computing vortices originating from blunt leading edges. The second case is that of a more complex fighter configuration. The superiority of the adapted solutions in capturing the vortex flow structure over the conventional unadapted results is demonstrated by comparisons with the wind-tunnel experimental data. The study shows that numerical prediction of vortical flows is highly sensitive to the local grid resolution and that the implementation of grid adaptation is essential when applying CFD methods to such complicated flow problems.
In supersonic and hypersonic flow computations, well-resolved boundary layers are essential for accurate quantification of surface heating and transition prediction, particularly via linear stability analysis. Grid design for hypersonic flows with shocks, boundary-layer separation, and/or complex mean flow features incorporating spanwise/azimuthal inhomogeneities is a difficult issue. In comparison to a manual grid adaptation procedure, an autonomous grid adaptation technique offers significant improvements in computing time and solution quality. The VULCAN-CFD solver already includes a validated procedure for automatic adaption of structured grids to the bow shock. The present focus is on implementing an automatic boundary-layer adaptation capability in VULCAN-CFD that adapts structured, multiblock grids to both the bow shock and the boundary layer at the same time. The boundary-layer adaptation algorithm allows the user to specify the number of cells within the boundary layer, along with the input parameters used for detecting the edge of the boundary layer, namely, the variable used in the edge detection criterion,the edge detection method, the detection direction, and the relaxation factor used during the morphing of the grid. The algorithm automatically distributes grid points along the wall-normal direction to achieve a smooth variation in grid spacing from the edge of the boundary layer to a"junction" location within the outer part of the grid. Illustrative results are presented for three different high-speed configurations: the two-dimensional flow over a cylinder at Mach 17.6 and unit Reynolds number of Re=0.38x10^6 m^-1, the axisymmetric flow over a cone-cylinder-flare model at Mach 6.0 and Re = 10.5×10^6 m^-1, and the three-dimensional flow over a blunt, 7-degree half-angle cone at 5-degree angle of attack in a Mach 9.79 flow with Re = 17.1×10^6 m^-1. The automated boundary-layer adaptation is shown to provide an adequate grid topology that is aligned with the bow shock in the outer part of the grid and also resolves the viscous boundary-layer region close to the surface.
A scheme is presented for the automatic generation of structured overset meshes ongeometries that are defined by Boundary Representation (BRep) solids. The surface meshsystem consists of face, edge and node meshes corresponding to the three respective basicBRep entity types. A cut-cell method is introduced to improve robustness of the on-geometrydetermination test for a face mesh grid point. A geometric component tagging scheme is utilizedto enhance local grid point distribution on a configuration with a large range of geometricscales. A cap grid topology is automatically utilized around the trailing edge of wing and tailtips to enhance mesh quality and to enable more effective surface coverage. Robustness of thehyperbolic surface marching method is improved by replacing the point projection scheme witha walking scheme. Relaxation of surface grid spacing at concave corners enables automatedhigh quality hyperbolic volume mesh creation. Domain connectivity is automatically performedon the surface mesh system. A variety of test cases are presented including a re-entry capsule,two models of the Juncture Flow Experiment wing-body, five rotorcraft concept vehicles, andvarious components of the High-Lift Common Research Model from the High-Lift PredictionWorkshop 4.
Geometry modeling and grid generation (GMGG) have played and will continue to play an important role in computational aerosciences. During the past two decades, tremendous progress has occurred in GMGG; however, GMGG is still the biggest bottleneck to routine applications for complicated Computational Fluid Dynamics (CFD) and Computational Structures Mechanics (CSM) models for analysis, design, and optimization. We are still far from incorporating GMGG tools in a design and optimization environment for complicated configurations. It is still a challenging task to parameterize an existing model in today's Computer-Aided Design (CAD) systems, and the models created are not always good enough for automatic grid generation tools. Designers may believe their models are complete and accurate, but unseen imperfections (e.g., gaps, unwanted wiggles, free edges, slivers, and transition cracks) often cause problems in gridding for CSM and CFD. Despite many advances in grid generation, the process is still the most labor-intensive and time-consuming part of the computational aerosciences for analysis, design, and optimization. In an ideal design environment, a design engineer would use a parametric model to evaluate alternative designs effortlessly and optimize an existing design for a new set of design objectives and constraints. For this ideal environment to be realized, the GMGG tools must have the following characteristics: (1) be automated, (2) provide consistent geometry across all disciplines, (3) be parametric, and (4) provide sensitivity derivatives. This paper will review the status of GMGG for analysis, design, and optimization processes, and it will focus on some emerging ideas that will advance the GMGG toward the ideal design environment.