Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “large 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 1,117 records · Page 62

Cloud Regimes as a Tool for Systematic Study of Various Aerosol-Cloud-Precipitation Interactions

Systematic changes of clouds and precipitation are notoriously difficult to ascribe to aerosols. This presentation will showcase yet one more attempt to at least credibly detect the signal of aerosol-cloud-precipitation interactions. We surmise that the concept of cloud regimes (CRs) is appropriate to conduct such an investigation. Previous studies focused on what we call here dynamical CRs, and while we continue to adopt those too for our analysis, we have found that a different way of organizing cloud systems, namely via microphysical regimes is also promising. Our analysis relies on MODIS Collection 6 Level-3 data for clouds and aerosols, and TRMM-TMPA data for precipitation. The regimes are derived by applying clustering analysis on MODIS joint histograms, and once each grid cell is assigned a regime, aerosol and precipitation data can be spatiotemporally matched and composited by regime. The composites of various cloud and precipitation variables for high (upper quartile of distribution) and low (lower quartile) aerosol loadings can then be contrasted. We seek evidence of aerosol effects both in regimes with large fractions of deep ice-rich clouds, as well as regimes where low liquid phase clouds dominate. Signals can be seen, especially when the analysis is broken by land-ocean and when additional filters are applied, but there are of course caveats which will be discussed.

cloud regimes↗

Composite sizing and ply orientation for stiffness requirements using a large finite element structural model

A NASTRAN bulk dataset preprocessor was developed to facilitate the integration of filamentary composite laminate properties into composite structural resizing for stiffness requirements. The NASCOMP system generates delta stiffness and delta mass matrices for input to the flutter derivative program. The flutter baseline analysis, derivative calculations, and stiffness and mass matrix updates are controlled by engineer defined processes under an operating system called CBUS. A multi-layered design variable grid system permits high fidelity resizing without excessive computer cost. The NASCOMP system uses ply layup drawings for basic input. The aeroelastic resizing for stiffness capability was used during an actual design exercise.

Radovcich, N. A.↗

Sampling Versus Filtering in Large-Eddy Simulations

A LES formalism in which the filter operator is replaced by a sampling operator is proposed. The unknown quantities that appear in the LES equations originate only from inadequate resolution (Discretization errors). The resulting viewpoint seems to make a link between finite difference approaches and finite element methods. Sampling operators are shown to commute with nonlinearities and to be purely projective. Moreover, their use allows an unambiguous definition of the LES numerical grid. The price to pay is that sampling never commutes with spatial derivatives and the commutation errors must be modeled. It is shown that models for the discretization errors may be treated using the dynamic procedure. Preliminary results, using the Smagorinsky model, are very encouraging.

Debliquy, O.↗

Trends and pacing items in computational aerodynamics

A perspective is presented of trends in computational aerodynamics, and of important technology development items that pace future advanced applications. From a survey of AIAA Journal papers published during the past two decades, the growth trends and the progressively increasing emphasis on code development for viscous, compressible, turbulent flow are illustrated. These trends are reflected in the chronology of introduction by the aerospace industry of new computational methods in aircraft design. Key pacing items outlined are: automatic grid generation for nonlinear inviscid computations; advanced computers, improved efficiency of numerical methods, and improved turbulence models for Reynolds-averaged Navier-Stokes computations; advanced computers, time-dependent three-dimensional law-of-the-wall, code development, improved efficiency of numerical methods, and improved subgrid-scale turbulence modeling for large eddy simulations.

Chapman, D. R.↗

Computation of the tip vortex flowfield for advanced aircraft propellers

The tip vortex flowfield plays a significant role in the performance of advanced aircraft propellers. The flowfield in the tip region is complex, three-dimensional and viscous with large secondary velocities. An analysis is presented using an approximate set of equations which contains the physics required by the tip vortex flowfield, but which does not require the resources of the full Navier-Stokes equations. A computer code was developed to predict the tip vortex flowfield of advanced aircraft propellers. A grid generation package was developed to allow specification of a variety of advanced aircraft propeller shapes. Calculations of the tip vortex generation on an SR3 type blade at high Reynolds numbers were made using this code and a parametric study was performed to show the effect of tip thickness on tip vortex intensity. In addition, calculations of the tip vortex generation on a NACA 0012 type blade were made, including the flowfield downstream of the blade trailing edge. Comparison of flowfield calculations with experimental data from an F4 blade was made. A user's manual was also prepared for the computer code (NASA CR-182178).

Tsai, Tommy M.↗

The effect of severe storms on the ice cover of the northern Tatarskiy Strait

Passive microwave images from the Special Sensor Microwave Imager are used to study the volume of ice and sea-bottom water in the Japan Sea as affected by winds and severe storms. The data set comprises brightness temperatures gridded on a polar stereographic projection, and the processing is accomplished with a linear algorithm by Cavalieri et al. (1983) based on the vertically polarized 37-GHz channel. The expressions for calculating heat fluxes and downwelling radiation are given, and ice-cover fluctuations are correlated with severe storm events. The storms generate large transient polynya that occur simultaneously with the strongest heat fluxes, and severe storms are found to contribute about 25 percent of the annual introduction of 25 cu km of ice in the region. The ice production could lead to the renewal of enough sea-bottom water to account for the C-14 data provided, and the generation of Japan Sea bottom water is found to vary directly with storm activity.

Martin, Seelye↗

The Diurnal Cycle in TOGA-COARE: Regional Scale Model Simulations

The diurnal variation of precipitation processes over the tropics is a well-known phenomenon and has been studied using surface rainfall data, radar reflectivity data, and satellite-derived cloudiness and precipitation. Recently, analyzed observations from Tropical Oceans and Global Atmosphere Coupled Ocean-Atmosphere Response Experiment (TOGA COARE) in the tropical western Pacific ocean to study the relevant mechanisms producing diurnal variation of precipitation. They found that the diurnal Sea surface temperature (SST) cycle is important for afternoon showers in the undisturbed periods and diurnal radiative processes for nocturnal rainfall. Cloud resolving models (CRMS) have been used to determine the mechanisms associated with diurnal variation of precipitating processes. CRMs allow explicit cloud-radiation and air-sea interactive processes. However, CRMs can be only used for idealized simulations (i.e., no feedback between clouds and their embedded large-scale environments; cyclic lateral boundary conditions and idealized initial conditions). In this study, the Penn State/NCAR Mesoscale Model (MM5) with improved physics (i.e., cloud microphysics, radiation, land-soil-vegetation-surface processes, and TOGA COARE flux scheme) and a multiple level nesting technique (covers the TOGA COARE LSA/IFA with a 54 km grid and can nest down to 18, 6 and possibly even 2 km) will be adopted for studying the diurnal variations of rainfall. We will examine precipitation processes over open ocean and over land. We will also perform sensitivity tests to determine how the radiative forcing and diurnal SST cycle affects the development of convection.

Tao, W.-K.↗

Implementation and Characterization of Three-Dimensional Particle-in-Cell Codes on Multiple-Instruction-Multiple-Data Massively Parallel Supercomputers

A three-dimensional electrostatic particle-in-cell (PIC) plasma simulation code has been developed on coarse-grain distributed-memory massively parallel computers with message passing communications. Our implementation is the generalization to three-dimensions of the general concurrent particle-in-cell (GCPIC) algorithm. In the GCPIC algorithm, the particle computation is divided among the processors using a domain decomposition of the simulation domain. In a three-dimensional simulation, the domain can be partitioned into one-, two-, or three-dimensional subdomains ("slabs," "rods," or "cubes") and we investigate the efficiency of the parallel implementation of the push for all three choices. The present implementation runs on the Intel Touchstone Delta machine at Caltech; a multiple-instruction-multiple-data (MIMD) parallel computer with 512 nodes. We find that the parallel efficiency of the push is very high, with the ratio of communication to computation time in the range 0.3%-10.0%. The highest efficiency (> 99%) occurs for a large, scaled problem with 64(sup 3) particles per processing node (approximately 134 million particles of 512 nodes) which has a push time of about 250 ns per particle per time step. We have also developed expressions for the timing of the code which are a function of both code parameters (number of grid points, particles, etc.) and machine-dependent parameters (effective FLOP rate, and the effective interprocessor bandwidths for the communication of particles and grid points). These expressions can be used to estimate the performance of scaled problems--including those with inhomogeneous plasmas--to other parallel machines once the machine-dependent parameters are known.

Lyster, P. M.↗

Software and Algorithms for Biomedical Image Data Processing and Visualization

A new software equipped with novel image processing algorithms and graphical-user-interface (GUI) tools has been designed for automated analysis and processing of large amounts of biomedical image data. The software, called PlaqTrak, has been specifically used for analysis of plaque on teeth of patients. New algorithms have been developed and implemented to segment teeth of interest from surrounding gum, and a real-time image-based morphing procedure is used to automatically overlay a grid onto each segmented tooth. Pattern recognition methods are used to classify plaque from surrounding gum and enamel, while ignoring glare effects due to the reflection of camera light and ambient light from enamel regions. The PlaqTrak system integrates these components into a single software suite with an easy-to-use GUI (see Figure 1) that allows users to do an end-to-end run of a patient s record, including tooth segmentation of all teeth, grid morphing of each segmented tooth, and plaque classification of each tooth image. The automated and accurate processing of the captured images to segment each tooth [see Figure 2(a)] and then detect plaque on a tooth-by-tooth basis is a critical component of the PlaqTrak system to do clinical trials and analysis with minimal human intervention. These features offer distinct advantages over other competing systems that analyze groups of teeth or synthetic teeth. PlaqTrak divides each segmented tooth into eight regions using an advanced graphics morphing procedure [see results on a chipped tooth in Figure 2(b)], and a pattern recognition classifier is then used to locate plaque [red regions in Figure 2(d)] and enamel regions. The morphing allows analysis within regions of teeth, thereby facilitating detailed statistical analysis such as the amount of plaque present on the biting surfaces on teeth. This software system is applicable to a host of biomedical applications, such as cell analysis and life detection, or robotic applications, such as product inspection or assembly of parts in space and industry.

Talukder, Ashit↗

Sinter of uniform, predictable, blemish-free nickel plaque for large aerospace nickel cadmium cells

A series of nickel slurry compositions were tested. Important slurry parameters were found to be the nature of the binder, a pore former and the method of mixing. A slow roll mixing which is non-turbulent successfully eliminated entrapped air so that bubbles and pockets were avoided in the sinter. A slurry applicator was developed which enabled an equal quantity of slurry to be applied to both sides of the grid. Sintering in a furnace having a graded atmosphere characteristic, ranging from oxidizing to strongly reducing, improved adhesion of porous sinter to grid and resulted in a uniform welding of nickel particles to each other throughout the plaque. Sintering was carried out in a horizontal furnace having three heating zones and 16 heating control circuits. Tests used for plaque evaluation include (1) appearance, (2) grid location and adhesion, (3) mechanical strength, (4) thickness, (5) weight per unit area, (6) void volume per unit area, (7) surface area and (8) electrical resistance. Plaque material was impregnated using Heliotek proprietary processes and 100 AH cells were fabricated.

Seiger, H. N.↗

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↗

The Analysis and Design of Low Boom Configurations Using CFD and Numerical Optimization Techniques

The use of computational fluid dynamics (CFD) for the analysis of sonic booms generated by aircraft has been shown to increase the accuracy and reliability of predictions. CFD takes into account important three-dimensional and nonlinear effects that are generally neglected by modified linear theory (MLT) methods. Up to the present time, CFD methods have been primarily used for analysis or prediction. Some investigators have used CFD to impact the design of low boom configurations using trial and error methods. One investigator developed a hybrid design method using a combination of Modified Linear Theory (e.g. F-functions) and CFD to provide equivalent area due to lift driven by a numerical optimizer to redesign or modify an existing configuration to achieve a shaped sonic boom signature. A three-dimensional design methodology has not yet been developed that completely uses nonlinear methods or CFD. Constrained numerical optimization techniques have existed for some time. Many of these methods use gradients to search for the minimum of a specified objective function subject to a variety of design variable bounds, linear and nonlinear constraints. Gradient based design optimization methods require the determination of the objective function gradients with respect to each of the design variables. These optimization methods are efficient and work well if the gradients can be obtained analytically. If analytical gradients are not available, the objective gradients or derivatives with respect to the design variables must be obtained numerically. To obtain numerical gradients, say, for 10 design variables, might require anywhere from 10 to 20 objective function evaluations. Typically, 5-10 global iterations of the optimizer are required to minimize the objective function. In terms of using CFD as a design optimization tool, the numerical evaluation of gradients can require anywhere from 100 to 200 CFD computations per design for only 10 design variables. If one CFD computation requires an hour of computational time on a Cray computer, one can see that the use of constrained numerical optimization quickly becomes impractical.Hence, in order to practically couple a numerical design optimization technique with a CFD method, the CFD method must be extremely efficient with running times on the order of only minutes. The CFD Euler code developed under NASA sponsorship and referred to as MIM3D-SB for the most part fulfills these efficiency requirements. Analysis of wing- body configurations can be computed in a matter of a few minutes. The present study will concentrate on the feasibility of the use of this CFD code in conjunction with a numerical design optimization technique for the sonic boom reduction of candidate HSCT configurations. A preliminary supersonic aircraft design system has been established that utilizes the numerical design optimization code NPSOL developed at Stanford University coupled with the supersonic NUM3D-SB CFD code. Many questions still need to be answered in regard to using CFD and numerical optimizers as design tools. There are difficulties related to both the CFD codes and the numerical optimizers. Numerical optimizers can converge to a local minima rather than a global minima. This behavior is largely a function of the initial guess in the design space. The optimizer also is searching for a minimum of the function in terms of its derivative without any regard to the actual function value. Numerically (i.e. CFD) determined gradients can also generate spurious numerical local minima. In addition, for the sonic boom problem, grid fineness will also determine the accuracy of the final design solution. Design optimization methods work well on problems defined by continuous objective functions. The sonic boom signature design problem is not necessarily defined by a continuous objective function. The signature can have a variety of shapes; i.e. from N-wave to multiple shocks. The far-field or ground signature may not transition continuously from one shape to another and hence, may exhibit discontinuous behavior. This is also a source of difficulty in using design optimization methods.In the following sections, several low boom and one reference aircraft configuration will be analyzed to predict their sonic boom signature characteristics. Modifications to some of these designs will also be presented to demonstrate the feasibility of using CFD as a design tool and to demonstrate the feasibility of designing shaped sonic boom signatures. Design modifications to some configurations will be presented to demonstrate the feasibility of achieving shaped signatures with reduced levels and not necessarily to represent realistic or aerodynamically efficient design modifications. Fuselage volume or camber are used as design variables in order to have a minimal effect on the primary wing aerodynamics. The paper will also seek to demonstrate whether a hybrid or ramped signature is feasible to achieve. For the low-boom configurations, the CFD predicted signatures will be compared qualitatively to their MLT design signatures.

Siclari, Michael J.↗

A Priori Subgrid Analysis of Temporal Mixing Layers with Evaporating Droplets

Subgrid analysis of a transitional temporal mixing layer with evaporating droplets has been performed using three sets of results from a Direct Numerical Simulation (DNS) database, with Reynolds numbers (based on initial vorticity thickness) as large as 600 and with droplet mass loadings as large as 0.5. In the DNS, the gas phase is computed using a Eulerian formulation, with Lagrangian droplet tracking. The Large Eddy Simulation (LES) equations corresponding to the DNS are first derived, and key assumptions in deriving them are first confirmed by computing the terms using the DNS database. Since LES of this flow requires the computation of unfiltered gas-phase variables at droplet locations from filtered gas-phase variables at the grid points, it is proposed to model these by assuming the gas-phase variables to be the sum of the filtered variables and a correction based on the filtered standard deviation; this correction is then computed from the Subgrid Scale (SGS) standard deviation. This model predicts the unfiltered variables at droplet locations considerably better than simply interpolating the filtered variables. Three methods are investigated for modeling the SGS standard deviation: the Smagorinsky approach, the Gradient model and the Scale-Similarity formulation. When the proportionality constant inherent in the SGS models is properly calculated, the Gradient and Scale-Similarity methods give results in excellent agreement with the DNS.

Okongo, Nora↗

A Comparative Study of Three Methodologies for Modeling Dynamic Stall

During the past two decades, there has been an increased reliance on the use of computational fluid dynamics methods for modeling rotors in high speed forward flight. Computational methods are being developed for modeling the shock induced loads on the advancing side, first-principles based modeling of the trailing wake evolution, and for retreating blade stall. The retreating blade dynamic stall problem has received particular attention, because the large variations in lift and pitching moments encountered in dynamic stall can lead to blade vibrations and pitch link fatigue. Restricting to aerodynamics, the numerical prediction of dynamic stall is still a complex and challenging CFD problem, that, even in two dimensions at low speed, gathers the major difficulties of aerodynamics, such as the grid resolution requirements for the viscous phenomena at leading-edge bubbles or in mixing-layers, the bias of the numerical viscosity, and the major difficulties of the physical modeling, such as the turbulence models, the transition models, whose both determinant influences, already present in static maximal-lift or stall computations, are emphasized by the dynamic aspect of the phenomena.

Sankar, L.↗

Software for Viewing Landsat Mosaic Images

A Windows-based computer program has been written to enable novice users (especially educators and students) to view images of large areas of the Earth (e.g., the continental United States) generated from image data acquired in the Landsat observations performed circa the year 1990. The large-area images are constructed as mosaics from the original Landsat images, which were acquired in several wavelength bands and each of which spans an area (in effect, one tile of a mosaic) of .5 in latitude by .6 in longitude. Whereas the original Landsat data are registered on a universal transverse Mercator (UTM) grid, the program converts the UTM coordinates of a mouse pointer in the image to latitude and longitude, which are continuously updated and displayed as the pointer is moved. The mosaic image currently on display can be exported as a Windows bitmap file. Other images (e.g., of state boundaries or interstate highways) can be overlaid on Landsat mosaics. The program interacts with the user via standard toolbar, keyboard, and mouse user interfaces. The program is supplied on a compact disk along with tutorial and educational information.

Watts, Zack↗

Comparison of Aero-Propulsive Performance Predictions for Distributed Propulsion Configurations

NASA's X-57 "Maxwell" flight demonstrator incorporates distributed electric propulsion technologies in a design that will achieve a significant reduction in energy used in cruise flight. A substantial portion of these energy savings come from beneficial aerodynamic-propulsion interaction. Previous research has shown the benefits of particular instantiations of distributed propulsion, such as the use of wingtip-mounted cruise propellers and leading edge high-lift propellers. However, these benefits have not been reduced to a generalized design or analysis approach suitable for large-scale design exploration. This paper discusses the rapid, "design-order" toolchains developed to investigate the large, complex tradespace of candidate geometries for the X-57. Due to the lack of an appropriate, rigorous set of validation data, the results of these tools were compared to three different computational flow solvers for selected wing and propulsion geometries. The comparisons were conducted using a common input geometry, but otherwise different input grids and, when appropriate, different flow assumptions to bound the comparisons. The results of these studies showed that the X-57 distributed propulsion wing should be able to meet the as-designed performance in cruise flight, while also meeting or exceeding targets for high-lift generation in low-speed flight.

Borer, Nicholas K.↗

Comparisons of Mixing Efficiency for the Strut Fuel Injector Obtained from Large-Eddy and Reynolds-Averaged Simulations, and Experiments

Mixing efficiency is obtained for a strut fuel injector at hypervelocity flow conditions by using large-eddy simulations (LES), Reynolds-averaged simulations (RAS), and experiments. The injector and flow conditions have been previously investigated by using RAS and experiments as a part of the Enhanced Injection and Mixing Project (EIMP) at the NASA Langley Research Center (LaRC). Because the fidelity of LES is a strong function of the grid, the mixing efficiency is obtained on two grids, the coarser of which is a factor of two coarser in each of the three dimensions with respect to the fine grid. The RAS uses the two-equation linear eddy viscosity and diffusivity modeling of Menter. In RAS, the species diffusivity model exhibits a strong dependence on the turbulent Schmidt number, which is often adjusted until some metric of engineering interest, such as the mixing efficiency, matches the experimental data. In the absence of experimental data, scale-resolving simulations, such as LES, have been proposed as surrogates for experiments that could provide the data needed to “calibrate” the turbulent Schmidt number in the RAS models. This approach is followed because LES requires significantly more computational resources (CPU, data storage, and time) than RAS, making it prohibitive for use in many engineering applications and specifically for parameter exploration or optimization. Here we examine the mixing efficiency obtained from several RAS with different values of the turbulent Schmidt number, and compare the results with those obtained from the LES and experiments. In addition, the least squares fitting approach was used to demonstrate how to obtain an estimate for the turbulent Schmidt number from LES analytically. These estimates were then used together with prior knowledge about RAS model sensitivity to select a turbulence model that was expected to best match the LES data.

LES↗

Comparisons of Mixing Efficiency for the Strut Fuel Injector Obtained from Large-Eddy and Reynolds-Averaged Simulations, and Experiments

Mixing efficiency is obtained for a strut fuel injector at hypervelocity flow conditions by using large-eddy simulations (LES), Reynolds-averaged simulations (RAS), and experiments. The injector and flow conditions have been previously investigated by using RAS and experiments as a part of the Enhanced Injection and Mixing Project (EIMP) at the NASA Langley Research Center (LaRC). Because the fidelity of LES is a strong function of the grid, the mixing efficiency is obtained on two grids, the coarser of which is a factor of two coarser in each of the three dimensions with respect to the fine grid. The RAS uses the two-equation linear eddy viscosity and diffusivity modeling of Menter. In RAS, the species diffusivity model exhibits a strong dependence on the turbulent Schmidt number, which is often adjusted until some metric of engineering interest, such as the mixing efficiency, matches the experimental data. In the absence of experimental data, scale-resolving simulations, such as LES, have been proposed as surrogates for experiments that could provide the data needed to “calibrate” the turbulent Schmidt number in the RAS models. This approach is followed because LES requires significantly more computational resources (CPU, data storage, and time) than RAS, making it prohibitive for use in many engineering applications and specifically for parameter exploration or optimization. Here we examine the mixing efficiency obtained from several RAS with different values of the turbulent Schmidt number, and compare the results with those obtained from the LES and experiments. In addition, the least squares fitting approach was used to demonstrate how to obtain an estimate for the turbulent Schmidt number from LES analytically. These estimates were then used together with prior knowledge about RAS model sensitivity to select a turbulence model that was expected to best match the LES data.

LES↗