Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “weak grids”

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

Flow Quality and Boundary Layer Transition

The widely held view is that transition to turbulence in the Blasius boundary layer occurs via amplification and eventual nonlinear breakdown of initially small amplitude instabilities i.e. Tollmien-Schlichting (TS) waves. However this scenario is only observed for low amplitude free-stream turbulence levels, i.e. u/U < 0.1%. Bypass of linear TS instability mechanism occurs for higher EST levels, yet considerable differences exist between the few experiments carefully designed to assess the effect of EST on transition. The consensus is that EST leads to longitudinal streaks that form near the leading edge in the boundary layer . These streaks appeal to be regions of concentrated streamwise vorticity and they are often referred to as Klebanoff modes. The importance of mean flow free-stream nonuniformity (FSN) is not as widely appreciated as EST for characterizing wind tunnel flow quality. Here it is shown that, although the v like generated by a d=50micron wire located upstream of the contraction (Re(sub d)=6.6, x/d=45,000) is immeasurably small by the time it interacts with the leading edge in the test section, it is responsible for generation of a pair of weak streamwise vortices in the boundary layer downstream. The characteristics of these wake-induced vortices and their effect on TS waves are demonstrated. Small remnant FSN variations are also shown to exist downstream of a turbulence grid. The question arises Are the adverse effects introduced by the turbulence grid caused by FST or by small remnant FSN variations?

Watmuff, Jonathan H↗

The Determination of the Large-Scale Circulation of the Pacific Ocean from Satellite Altimetry using Model Green's Functions

A Green's function method for obtaining an estimate of the ocean circulation using both a general circulation model and altimetric data is demonstrated. The fundamental assumption is that the model is so accurate that the differences between the observations and the model-estimated fields obey a linear dynamics. In the present case, the calculations are demonstrated for model/data differences occurring on very a large scale, where the linearization hypothesis appears to be a good one. A semi-automatic linearization of the Bryan/Cox general circulation model is effected by calculating the model response to a series of isolated (in both space and time) geostrophically balanced vortices. These resulting impulse responses or 'Green's functions' then provide the kernels for a linear inverse problem. The method is first demonstrated with a set of 'twin experiments' and then with real data spanning the entire model domain and a year of TOPEX/POSEIDON observations. Our present focus is on the estimate of the time-mean and annual cycle of the model. Residuals of the inversion/assimilation are largest in the western tropical Pacific, and are believed to reflect primarily geoid error. Vertical resolution diminishes with depth with 1 year of data. The model mean is modified such that the subtropical gyre is weakened by about 1 cm/s and the center of the gyre shifted southward by about 10 deg. Corrections to the flow field at the annual cycle suggest that the dynamical response is weak except in the tropics, where the estimated seasonal cycle of the low-latitude current system is of the order of 2 cm/s. The underestimation of observed fluctuations can be related to the inversion on the coarse spatial grid, which does not permit full resolution of the tropical physics. The methodology is easily extended to higher resolution, to use of spatially correlated errors, and to other data types.

Stammer, Detlef↗

Cloud-Precipitation Hybrid Regimes and their Projection onto IMERG Precipitation Data

We extend and enhance the concept of the Cloud Regimes (CRs) developed from two-dimensional joint histograms of cloud optical thickness and cloud top pressure from the Moderate Resolution Imaging Spectroradiometer (MODIS), by adding precipitation information in order to better understand cloud-precipitation relationships. Taking advantage of the high-resolution Integrated Multi-satellitE Retrievals for GPM (IMERG) precipitation dataset, cloud-precipitation “hybrid” regimes are derived by implementing the k-means clustering algorithm with advanced initialization and objective measures to determine the most optimal clusters. By expressing precipitation rates within 1-degree grid cell as histograms and making choices on the relative weight of cloud and precipitation, we could obtain several editions of hybrid cloud-precipitation regimes (CPRs), and examine their characteristics. In the deep tropics, when precipitation is weighted weakly, the cloud part of the hybrid entroids resembles the centroid of cloud-only regimes, but still tightens the cloud-precipitation relationship by decreasing the precipitation variability of each regime. As precipitation weight progressively increases, the shape of the cloudy part of the hybrid centroids becomes blunter, while the precipitation part of the centroids sharpens. In the case where cloud and precipitation are weighted equally, the CPRs representing high clouds with intermediate to heavy precipitation exhibit distinct features in the precipitation parts of the centroids, which allows us to project them onto the 30-minly IMERG domain. Such a projection can be used to overcome the temporal sparseness of MODIS cloud observations, which leads to great application potential for various convection-focused studies, including diurnal cycle analysis.

cloud-precipitation↗

Moire deflectometry with deferred analysis

A simplified version of moire deflectometry for flow visualization is presented. The modification comprises replacement of the second grating, the ground glass screen, and the camera with a plate film holder to obtain a shadow picture. Postanalysis proceeds by placing a grid transparency of the same periodicity of the remaining glass grating on the negative for projection viewing. The rotational angle between the grating and the plate alters the moire fringe periodicity. Use of the method in studying flows is predicted to reveal weak density gradient areas and shocks.

Rhodes, D. B.↗

Relationship between clouds and sea surface temperatures in the western tropical Pacific

Analysis of four years of earth radiation budget, cloud, and sea surface temperature data confirms that cloud parameters change dramatically when and where sea surface temperatures increase above approximately 300K. These results are based upon monthly mean values within 2.5 deg x 2.5 deg grid points over the 'warm pool' region of the western tropical Pacific. The question of whether sea surface temperatures are influenced, in turn, by the radiative effects of these clouds (Ramanathan and Collins) is less clear. Such a feedback, if it exists, is weak. The reason why clouds might have so little influence, despite large changes in their longwave and shortwave radiative effects, might be that the sea surface responds to both the longwave heating and the shortwave cooling effects of clouds, and the two effects nearly cancel. There are strong correlations between the rate of change of sea surface temperature and any of the radiation budget parameters that are highly correlated with the incident solar flux-implying that season and latitude are the critical factors determining sea surface temperatures. With the seasonal or both seasonal and latitudinal variations removed, the rate of change of sea surface temperature shows no correlation with cloud-related parameters in the western tropical Pacific.

Arking, Albert↗

Validation of a Node-Centered Wall Function Model for the Unstructured Flow Code FUN3D

In this paper, the implementation of two wall function models in the Reynolds averaged Navier-Stokes (RANS) computational uid dynamics (CFD) code FUN3D is described. FUN3D is a node centered method for solving the three-dimensional Navier-Stokes equations on unstructured computational grids. The first wall function model, based on the work of Knopp et al., is used in conjunction with the one-equation turbulence model of Spalart-Allmaras. The second wall function model, also based on the work of Knopp, is used in conjunction with the two-equation k-! turbulence model of Menter. The wall function models compute the wall momentum and energy flux, which are used to weakly enforce the wall velocity and pressure flux boundary conditions in the mean flow momentum and energy equations. These wall conditions are implemented in an implicit form where the contribution of the wall function model to the Jacobian are also included. The boundary conditions of the turbulence transport equations are enforced explicitly (strongly) on all solid boundaries. The use of the wall function models is demonstrated on four test cases: a at plate boundary layer, a subsonic di user, a 2D airfoil, and a 3D semi-span wing. Where possible, different near-wall viscous spacing tactics are examined. Iterative residual convergence was obtained in most cases. Solution results are compared with theoretical and experimental data for several variations of grid spacing. In general, very good comparisons with data were achieved.

Carlson, Jan-Renee↗

Micro-Ramps for External Compression Low-Boom Inlets

The application of vortex generators for flow control in an external compression, axisymmetric, low-boom concept inlet was investigated using RANS simulations with three-dimensional (3-D), structured, chimera (overset) grids and the WIND-US code. The low-boom inlet design is based on previous scale model 1- by 1-ft wind tunnel tests and features a zero-angle cowl and relaxed isentropic compression centerbody spike, resulting in defocused oblique shocks and a weak terminating normal shock. Validation of the methodology was first performed for micro-ramps in supersonic flow on a flat plate with and without oblique shocks. For the inlet configuration, simulations with several types of vortex generators were conducted for positions both upstream and downstream of the terminating normal shock. The performance parameters included incompressible axisymmetric shape factor, separation area, inlet pressure recovery, and massflow ratio. The design of experiments (DOE) methodology was used to select device size and location, analyze the resulting data, and determine the optimal choice of device geometry. The optimum upstream configuration was found to substantially reduce the post-shock separation area but did not significantly impact recovery at the aerodynamic interface plane (AIP). Downstream device placement allowed for fuller boundary layer velocity profiles and reduced distortion. This resulted in an improved pressure recovery and massflow ratio at the AIP compared to the baseline solid-wall configuration.

Rybalko, Michael↗

Extension of transonic flow computational concepts in the analysis of cavitated bearings

An analogy between the mathematical modeling of transonic potential flow and the flow in a cavitating bearing is described. Based on the similarities, characteristics of the cavitated region and jump conditions across the film reformation and rupture fronts are developed using the method of weak solutions. The mathematical analogy is extended by utilizing a few computational concepts of transonic flow to numerically model the cavitating bearing. Methods of shock fitting and shock capturing are discussed. Various procedures used in transonic flow computations are adapted to bearing cavitation applications, for example, type differencing, grid transformation, an approximate factorization technique, and Newton's iteration method. These concepts have proved to be successful and have vastly improved the efficiency of numerical modeling of cavitated bearings.

Vijayaraghavan, D.↗

Extension of transonic flow computational concepts in the analysis of cavitated bearings

An analogy between the mathematical modeling of transonic potential flow and the flow in a cavitating bearing is described. Based on the similarities, characteristics of the cavitated region and jump conditions across the film reformation and rupture fronts are developed using the method of weak solutions. The mathematical analogy is extended by utilizing a few computational concepts of transonic flow to numerically model the cavitating bearing. Methods of shock fitting and shock capturing are discussed. Various procedures used in transonic flow computations are adapted to bearing cavitation applications, for example, type differencing, grid transformation, an approximate factorization technique, and Newton's iteration method. These concepts have proved to be successful and have vastly improved the efficiency of numerical modeling of cavitated bearings.

Vijayaraghavan, D.↗

HLPW-4/GMGW-3: Wall-Modeled LES and Lattice-Boltzmann Technology Focus Group Workshop Summary

A summary of the nine submissions to the Wall-Modeled LES and Lattice-Boltzmann(WMLESLB) Technical Focus Group (TFG) at the 4th High lift Prediction Workshop is provided. The focus of this TFG was to assess the current capabilities of WMLES and Lattice-Boltzmann methods on a complex high-lift configuration across a wide range of angles of attack. Preliminary analysis of the submitted data suggests that>250𝑀spatial degrees of freedom are needed to accurately predict pitching moments at high angles of attack due to large pressure gradients present on the outboard slat and main element for𝛼 >17◦(corrected for free-air).While some scatter is reported in pitching moments at the low-angles of attack (𝛼 <11◦) - likely caused by differences in flap separation possibly due to low Reynolds number effects - excellent agreement is observed between the submissions near the𝐶𝐿,max state. Objective superiority over steady state RANS can be seen in terms of lack of excess outboard separation; a majority of the good quality WMLES and LB submissions predict wedge-shaped separation patterns consistent with the experimental oil flow. Differences in the onset of stall mechanism in the free-air configuration for𝛼 >20◦is reported with two distinct topologies observed. Topology A is characterized by the onset of corner-flow separation which progressively grows to produce a pitch break in free-air with an angle of attack increase from𝛼=20.55◦to𝛼=21.47◦+𝜖where𝜖is a small perturbation (varying between submissions and likely to be within±0.3◦). Topology Bis characterized by boundary layer weakness emanating from the inboard side of the wing-pylon juncture substantially larger than any weakness in the wing-body juncture (due to corner-flow);submissions in this category do not show any tendency for occurrence of a pitch break for the free-air configuration within the vicinity of𝛼=21.47◦. The in-tunnel simulations submitted by3 participants using different discretizations, grids and closure models show excellent agreement with the experiment in terms of a) integrated loads, b) surface flow-topology, and c) mechanism for the onset of inboard stall. Further evidence is provided to demonstrate both qualitative and quantitative superiority of all 3 WMLES submissions over a single steady state RANS submission to the workshop.

TTT↗

Large eddy simulation of shock turbulence interaction

A nonconservative formulation of the energy equation (solving for internal energy) was used to perform large eddy simulations of compressible turbulence in Moin et al. due to its simplicity in implementing SGS models compared to the conservative formulation (solving for total energy). In problems with shocks in the domain, however, the total energy formulation is preferred due to its conservative nature. A conservative set of equations for the LES were derived from the nonconservative equations derived by Moin et al. Performance of the conservative formulation was compared with the experiment on decaying grid-generated turbulence as well as with the filtered DNS field. Various shock-capturing schemes were tested, and an ENO shock-capturing scheme of Shu and Osher was chosen for the simulation of shock/turbulence interaction. The scheme was tested and validated against the data base generated by DNS of weak shock waves. The results obtained with the essentially nonoscillatory (ENO) scheme were within 5 percent from the DNS results, and used less than 25 percent of the CPU time used in the DNS.

Lee, Sangsan↗

A Finite Element Theory for Predicting the Attenuation of Extended-Reacting Liners

A non-modal finite element theory for predicting the attenuation of an extended-reacting liner containing a porous facesheet and located in a no-flow duct is presented. The mathematical approach is to solve separate wave equations in the liner and duct airway and to couple these two solutions by invoking kinematic constraints at the facesheet that are consistent with a continuum theory of fluid motion. Given the liner intrinsic properties, a weak Galerkin finite element formulation with cubic polynomial basis functions is used as the basis for generating a discrete system of acoustic equations that are solved to obtain the coupled acoustic field. A state-of-the-art, asymmetric, parallel, sparse equation solver is implemented that allows tens of thousands of grid points to be analyzed. A grid refinement study is presented to show that the predicted attenuation converges. Excellent comparison of the numerically predicted attenuation to that of a mode theory (using a Haynes 25 metal foam liner) is used to validate the computational approach. Simulations are also presented for fifteen porous plate, extended-reacting liners. The construction of some of the porous plate liners suggest that they should behave as resonant liners while the construction of others suggest that they should behave as broadband attenuators. In each case the finite element theory is observed to predict the proper attenuation trend.

Watson, W. R.↗

Ultraviolet and optical studies of binaries with luminous cool primaries and hot companions. IV - Further IUE detections

We have obtained IUE spectra for 31 middle and late-type giant and supergiant stars whose TD-1 fluxes or ground-based spectra indicate the presence of a hot component, or whose radial velocities indicate an unseen component. Stellar components earlier than type F1 were detected in 22 cases. While 20 of the hot secondaries are seen weakly in optical spectra, two are UV discoveries: HD 58134 and HD 183864. The hot companions are classified accurately by temperature class from their far-UV spectra. The interstellar extinction of each system and the relative luminosities of the components are derived from the UV and optical fluxes, using a new grid of UV intrinsic colors for hot dwarfs. We find that many giant stars apparently have companions which are too hot and hence too luminous for consistency with the primary's spectral classification.

Parsons, Sidney B.↗

Comparison of TRMM and Global Precipitation Climatology Project (GPCP) Precipitation Analyses

This paper describes recent results of using Tropical Rainfall Measuring Mission (TRMM) (launched in November 1997) information as the key calibration tool in a merged analysis on a 1 x 1' latitude/longitude monthly scale based on multiple satellite sources and raingauge analyses. The TRMM-based product is compared with the community-based Global Precipitation Climatology Project (GPCP) results. The long-term GPCP analysis is compared to the new TRMM-based analysis which uses the most accurate TRMM information to calibrate the estimates from the Special Sensor Microwave/Imager (SSM/I) and geosynchronous IR observations and merges those estimates together with the TRMM and gauge information to produce accurate rainfall estimates with the increased sampling provided by the combined satellite information. The comparison with TRMM results on a month-to-month basis should clarify the strengths and weaknesses of the long-term GPCP product in the tropics and point to how to improve the monitoring analysis. Preliminary results from the TRMM merged satellite analysis indicates fairly close agreement with the GPCP estimates. The GPCP analysis is done at 2.5 degree latitude/longitude resolution and interpolated to a 1 degree grid for comparison with the TRMM analysis. As expected the same features are evident in both panels, but there are subtle differences in the magnitudes. Focusing on the Pacific Ocean Inter-Tropical Convergence Zone (ITCZ) one can see the TRMM-based estimates having higher peak values and lower values in the ITCZ periphery. These attributes also show up in the statistics, where GPCP>TRMM at low values (below 10 mm/d) and TRMM>GPCP at high values (greater than 15 mm/d). The area in the Indian Ocean which shows consistently higher values of TRMM over GPCP needs to be examined carefully to determine if the lack of geosynchronous data has led to a difference in the two analyses. By the time of the meeting over a year of TRMM products will be available for comparison. Global tropical and regional values will be compared. Both products will be compared to TRMM validation site data over land and water. The results should begin to determine the use of the TRMM estimates in the evaluation of the GPCP analysis.

Adler, Robert F.↗

Evaluation of Isolated Fuselage and Rotor-Fuselage Interaction Using CFD

The US Army Aeroflightdynamics Directorate (AFDD), the French Office National d'Etudes et de Recherches Aerospatiales (ONERA) and the Georgia Institute of Technology (GIT) are working under the United States/France Memorandum of Agreement on Helicopter Aeromechanics to study rotorcraft aeromechanics issues of interest to both nations. As a task under this agreement, a comparative study of the Dauphin 365N helicopter has been undertaken to analyze the capabilities and weaknesses of state-of-the-art computational fluid dynamics (CFD) codes, with the aim of fuselage performance prediction and investigation of rotor-fuselage interaction. Three CFD flow solvers applied on three meshes provide similar results in terms of pressure coefficient. Force predictions vary somewhat. This paper presents details on the grid sensitivity and the low Mach number preconditioning influence. The importance of taking into account the wind tunnel strut and the rotor hub is shown. The pressure coefficients along top and bottom centerlines of the fuselage are in good agreement with the experiment except in the area aft of the hub. There remains a discrepancy between the computed forces and the experimental data due in part to modeling inaccuracies. Rotor- fuselage interactions are performed using uniform and non-uniform actuator disk models in order to simulate the rotor downwash.

Renaud, Thomas↗

GMS-based"Future Time" Rainfall Data Assimilation for Mesoscale Weather Prediction over Korean Peninsula and Future Prospects with GPM Satellite Measurements

This study examines the use of satellite-derived nowcasted (short-term forecasted) rainfall over 3-hour time periods to gain an equivalent time increment in initializing a nonhydrostatic mesoscale model used for predicting convective rainfall events over the Korean peninsula. Infrared (IR) window measurements from the Japanese Geostationary Meteorological Satellite (GMS) are used to specify latent heating for a spinup period of the model - but in future time -- thus initializing in advance of actual time in the framework of a prediction scenario. The main scientific objective of the study is to investigate the strengths and weaknesses of this approach insofar as data assimilation, in which the nowcasted assimilation data are derived independently of the prognostic model itself. Although there have been various recent improvements in formulating the dynamics, thermodynamics, and microphysics of mesoscale models, as well as computer advances which allow the use of high resolution cloud-resolving grids and explicit latent heating over regional domains, spinup remains at the forefront of unresolved mesoscale modeling problems. In general, non-realistic spinup limits the skill in predicting the spatial-temporal distribution of convection and precipitation, primarily in the early hours of a. forecast, stemming from standard prognostic variables not representing the initial diabatic heating field produced by the ambient convection and cloud fields. The long-term goal of this research is to improve short-range (12-hour) quantitative precipitation forecasting (QPF) over the Korean peninsula through the use of innovative data assimilation methods based on geosynchronous satellite measurements. As a step in th~s direction, a non-standard data assimilation experiment in conjunction with GMS-retrieved nowcasted rainfall information introduced to the mesoscale model is conducted. The 3-hourly precipitation forecast information is assimilated through nudging the associated diabatic heating during the early stages of a forecast period. This procedure is expected to enhance details in the moisture field during model integration, and thus improve spinup performance, assuming the errors in the future time latent heating data ate less than intrinsic model background errors.

Smith, Eric A.↗

Electron lithography STAR design guidelines. Part 3: The mosaic transistor array applied to custom microprocessors. Part 4: Stores logic arrays, SLAs implemented with clocked CMOS

The Mosaic Transistor Array is an extension of the STAR system developed by NASA which has dedicated field cells designed to be specifically used in semicustom microprocessor applications. The Sandia radiation hard bulk CMOS process is utilized in order to satisfy the requirements of space flights. A design philosophy is developed which utilizes the strengths and recognizes the weaknesses of the Sandia process. A style of circuitry is developed which incorporates the low power and high drive capability of CMOS. In addition the density achieved is better than that for classic CMOS, although not as good as for NMOS. The basic logic functions for a data path are designed with compatible interface to the STAR grid system. In this manner either random logic or PLA type structures can be utilized for the control logic.

Trotter, J. D.↗

Navier-Stokes turbine heat transfer predictions using two-equation turbulence closures

Navier-Stokes calculations were carried out in order to predict the heat-transfer rates on turbine blades. The calculations were performed using TRAF2D which is a k-epsilon, explicit, finite volume mass-averaged Navier-Stokes solver. Turbulence was modeled using Coakley's q-omega and Chien's k-epsilon two-equation models and the Baldwin-Lomax algebraic model. The model equations along with the flow equations were solved explicitly on a nonperiodic C grid. Implicit residual smoothing (IRS) or a combination of multigrid technique and IRS was applied to enhance convergence rates. Calculations were performed to predict the Stanton number distributions on the first stage vane and blade row as well as the second stage vane row of the SSME high-pressure fuel turbine. The comparison serves to highlight the weaknesses of the turbulence models for use in turbomachinery heat-transfer calculations.

Ameri, Ali A.↗