Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Distribution grid”

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 865 records · Page 48

Relative efficiency and accuracy of two Navier-Stokes codes for simulating attached transonic flow over wings

Two codes which solve the 3-D Thin Layer Navier-Stokes (TLNS) equations are used to compute the steady state flow for two test cases representing typical finite wings at transonic conditions. Several grids of C-O topology and varying point densities are used to determine the effects of grid refinement. After a description of each code and test case, standards for determining code efficiency and accuracy are defined and applied to determine the relative performance of the two codes in predicting turbulent transonic wing flows. Comparisons of computed surface pressure distributions with experimental data are made.

Bonhaus, Daryl L.↗

A numerical code for a three-dimensional magnetospheric MHD equilibrium model

Two dimensional and three dimensional MHD equilibrium models were begun for Earth's magnetosphere. The original proposal was motivated by realizing that global, purely data based models of Earth's magnetosphere are inadequate for studying the underlying plasma physical principles according to which the magnetosphere evolves on the quasi-static convection time scale. Complex numerical grid generation schemes were established for a 3-D Poisson solver, and a robust Grad-Shafranov solver was coded for high beta MHD equilibria. Thus, the effects were calculated of both the magnetopause geometry and boundary conditions on the magnetotail current distribution.

Voigt, G.-H.↗

A Grid-resolved Analysis of Base Flowfield for a Four-Engine Clustered Nozzle Configuration

The objective of this study is to propose a computational methodology that can effectively anchor the base flowfield of a four-engine clustered nozzle configuration. This computational methodology is based on a three-dimensional, viscous flow, pressure-based computational fluid dynamics (CFD) formulation. For efficient CFD calculation, a Prandtl-Meyer solution treatment is applied to the algebraic grid lines for initial plume expansion resolution. As the solution evolves, the computational grid is adapted to the pertinent flow gradients. The CFD model employs an upwind scheme in which second- and fourth-order central differencing schemes with artificial dissipation are used. The computed quantitative base flow properties such as the radial base pressure distributions, model centerline static pressure, Mach number and impact pressure variations, and base pressure characteristic curve agreed reasonably well with those of the measurement.

Wang, Ten-See↗

Numerical Modeling of Liquid-Vapor Phase Change

We implemented a two- and three-dimensional finite difference/front tracking technique to solve liquid-vapor phase change problems. The mathematical and the numerical features of the method were explained in great detail in our previous reports, Briefly, we used a single formula representation which incorporated jump conditions into the governing equations. The interfacial terms were distributed as singular terms using delta functions so that the governing equations would be the same as conventional conservation equations away from the interface and in the vicinity of the interface they would provide correct jump conditions. We used a fixed staggered grid to discretize these equations and an unstructured grid to explicitly track the front. While in two dimensions the front was simply a connection of small line segments, in three dimensions it was represented by a connection of small triangular elements. The equations were written in conservative forms and during the course of computations we used regriding to control the size of the elements of the unstructured grid. Moreover, we implemented a coalescence in two dimensions which allowed the merging of different fronts or two segments of the same front when they were sufficiently close. We used our code to study thermocapillary migration of bubbles, burst of bubbles at a free surface, buoyancy-driven interactions of bubbles, evaporation of drops, rapid evaporation of an interface, planar solidification of an undercooled melt, dendritic solidification, and a host of other problems cited in the reference.

Esmaeeli, Asghar↗

Computation of Coupled Thermal-Fluid Problems in Distributed Memory Environment

The thermal-fluid coupling problems are very important to aerospace and engineering applications. Instead of analyzing heat transfer and fluid flow separately, this study merged two well-accepted engineering solution methods, SINDA for thermal analysis and FDNS for fluid flow simulation, into a unified multi-disciplinary thermal fluid prediction method. A fully conservative patched grid interface algorithm for arbitrary two-dimensional and three-dimensional geometry has been developed. The state-of-the-art parallel computing concept was used to couple SINDA and FDNS for the communication of boundary conditions through PVM (Parallel Virtual Machine) libraries. Therefore, the thermal analysis performed by SINDA and the fluid flow calculated by FDNS are fully coupled to obtain steady state or transient solutions. The natural convection between two thick-walled eccentric tubes was calculated and the predicted results match the experiment data perfectly. A 3-D rocket engine model and a real 3-D SSME geometry were used to test the current model, and the reasonable temperature field was obtained.

Wei, H.↗

Distributed Relaxation Multigrid and Defect Correction Applied to the Compressible Navier-Stokes Equations

The distributed-relaxation multigrid and defect- correction methods are applied to the two- dimensional compressible Navier-Stokes equations. The formulation is intended for high Reynolds number applications and several applications are made at a laminar Reynolds number of 10,000. A staggered- grid arrangement of variables is used; the coupled pressure and internal energy equations are solved together with multigrid, requiring a block 2x2 matrix solution. Textbook multigrid efficiencies are attained for incompressible and slightly compressible simulations of the boundary layer on a flat plate. Textbook efficiencies are obtained for compressible simulations up to Mach numbers of 0.7 for a viscous wake simulation.

Thomas, J. L.↗

Transformation of OODT CAS to Perform Larger Tasks

A computer program denoted OODT CAS has been transformed to enable performance of larger tasks that involve greatly increased data volumes and increasingly intensive processing of data on heterogeneous, geographically dispersed computers. Prior to the transformation, OODT CAS (also alternatively denoted, simply, 'CAS') [wherein 'OODT' signifies 'Object-Oriented Data Technology' and 'CAS' signifies 'Catalog and Archive Service'] was a proven software component used to manage scientific data from spaceflight missions. In the transformation, CAS was split into two separate components representing its canonical capabilities: file management and workflow management. In addition, CAS was augmented by addition of a resource-management component. This third component enables CAS to manage heterogeneous computing by use of diverse resources, including high-performance clusters of computers, commodity computing hardware, and grid computing infrastructures. CAS is now more easily maintainable, evolvable, and reusable. These components can be used separately or, taking advantage of synergies, can be used together. Other elements of the transformation included addition of a separate Web presentation layer that supports distribution of data products via Really Simple Syndication (RSS) feeds, and provision for full Resource Description Framework (RDF) exports of metadata.

Mattmann, Chris↗

A Coupled Fluid-Structure Interaction Analysis of Solid Rocket Motor with Flexible Inhibitors

A capability to couple NASA production CFD code, Loci/CHEM, with CFDRC's structural finite element code, CoBi, has been developed. This paper summarizes the efforts in applying the installed coupling software to demonstrate/investigate fluid-structure interaction (FSI) between pressure wave and flexible inhibitor inside reusable solid rocket motor (RSRM). First a unified governing equation for both fluid and structure is presented, then an Eulerian-Lagrangian framework is described to satisfy the interfacial continuity requirements. The features of fluid solver, Loci/CHEM and structural solver, CoBi, are discussed before the coupling methodology of the solvers is described. The simulation uses production level CFD LES turbulence model with a grid resolution of 80 million cells. The flexible inhibitor is modeled with full 3D shell elements. Verifications against analytical solutions of structural model under steady uniform pressure condition and under dynamic condition of modal analysis show excellent agreements in terms of displacement distribution and eigen modal frequencies. The preliminary coupled result shows that due to acoustic coupling, the dynamics of one of the more flexible inhibitors shift from its first modal frequency to the first acoustic frequency of the solid rocket motor.

Yang, H. Q.↗

Effects of Surface BRDF on the OMI Cloud and NO2 Retrievals: A New Approach Based on Geometry-Dependent Lambertian Equivalent Reflectivity (GLER) Derived from MODIS

The Ozone Monitoring Instrument (OMI) cloud and NO2 algorithms use a monthly gridded surface reflectivity climatology that does not depend upon the observation geometry. In reality, reflection of incoming direct and diffuse solar light from land or ocean surfaces is sensitive to the sun sensor geometry. This dependence is described by the bidirectional reflectance distribution function (BRDF). To account for the BRDF, we propose to use a new concept of geometry-dependent Lambertian equivalent reflectivity (GLER). Implementation within the existing OMI cloud and NO2 retrieval infrastructure requires changes only to the input surface reflectivity database. GLER is calculated using a vector radiative transfer model with high spatial resolution BRDF information from MODIS over land and the Cox Munk slope distribution over ocean with a contribution from water-leaving radiance. We compare GLER and climatological LER at 466 nm, which is used in the OMI O2-O2cloud algorithm to derive effective cloud fractions. A detailed comparison of the cloud fractions and pressures derived with climatological and GLERs is carried out. GLER and corresponding retrieved cloud products are then used as input to the OMI NO2 algorithm. We find that replacing the climatological OMI-based LERs with GLERs can increase NO2 vertical columns by up to 50 % in highly polluted areas; the differences include both BRDF effects and biases between the MODIS and OMI-based surface reflectance data sets. Only minor changes to NO2 columns (within 5 %) are found over unpolluted and overcast areas.

MODIS BRDF OMI clouds NO2 GLER↗

Radar Reflectivity Simulated by a 2-D Spectra Bin Model: Sensitivity of Cloud-aerosol Interaction

The Goddard Cumulus Ensemble (GCE) model with bin spectra microphysics is used to simulate mesoscale convective systems.The model uses explicit bins to represent size spectra of cloud nuclei, water drops, ice crystals, snow and graupel. Each hydrometeorite category is described by 33 mass bins. The simulations provide a unique data set of simulated raindrop size distribution in a realistic dynamic frame. Calculations of radar parameters using simulated drop size distribution serve as an evaluation of numerical model performance. In addition, the GCE bin spectra modes is a very useful tool to study uncertainties related to radar observations; all the environmental parameters are precisely known. In this presentation, we concentrate on the discussion of Z-R (ZDR-R) relation in the simulated systems. Due to computational limitations, the spectra bin model has been run in two dimensions with 31 stretched vertical layers and 1026 horizontal grid points (1 km resolution). Two different cases, one in midlatitude continent, the other in tropical ocean, have been simulated. The continental case is a strong convection which lasted for two hours. The oceanic case is a persistent system with more than 10 hours' life span. It is shown that the simulated Z-R (ZDR-R) relations generally agree with observations using radar and rain gauge data. The spatial and temporal variations of Z-R relation in different locations are also analyzed. Impact of aerosols on cloud formation and raindrop size distribution was studied. Both clean (low CCN) and dirty (high CCN) cases are simulated. The Z-R relation is shown to vary considerable in the initial CCN concentrations.

Li, Kiaowen↗

Linear trends in Northern Hemisphere tropospheric geopotential height and temperature patterns

Gridded National Meteorological Center data for 500 mb geopotential height and 300-500 mb and 500-700 mb thickness for the period 1951-1978 are subjected to linear trend analyses. The analyses are carried out for each calendar month. Significant geographical and seasonal distributions of cooling and warming patterns are found. An atmospheric cooling trend over the North Pacific during the winter months is seen in a region where oceanic cooling has also been observed, but planetary-wave adjustments rather than ocean-atmosphere feedback mechanisms appear to be dominant in the atmospheric cooling on climatic time scales. Consistently large temperature trends are also seen over the continent of Asia. Comparisons between thickness trends in the layer 300-500 mb and those in the layer 500-700 mb reveal pronounced patterns of stabilization and destabilization.

Reiter, E. R.↗

A viscous-inviscid interactive procedure for rotational flow in cascades of two-dimensional airfoils of arbitrary shape

A viscous-inviscid interactive calculation procedure is developed for application to flow in cascades of two-dimensional airfoils. This procedure has essentially three components. First, a numerical solution of the Euler equations which can accommodate an arbitrarily specified cascade geometry is carried out on a nonorthogonal curvileanear grid mesh that is fitted to the geometry of the cascade. A method of grid generation has been used which relies in part on a succession of conformal mappings. Second, a viscous solution for use in boundary layer and wake regions has been programmed. Finally, an interactive scheme which takes the form of a source-sink distribution along the blade surface and wake centerline is employed. Results have been obtained with this procedure for several cascade flow situations, and some comparisons with experiment are presented.

Johnston, W.↗

Tectonics of the Valhalla Structure on Callisto

The surficial morphology of Callisto is dominated by large concentric patterns of ridges, scarps, and furrows. The largest one, Valhalla, shows a smooth central area 350 km in radius, surrounded by concentric sinuous ridges which extend as far as 600km from the center. Beyond the ridged area, scarps and furrows are observed as far as 2000 km. The global azimuthal distribution of the scarps and furrows of the outer ring (beyond 700 km from the center) was analyzed. It is proposed that the local geometry and the global distribution of scarps and furrows around Valhalla are the results of the reactivation of an old pattern by the Valhalla event. The coincidence of these local directions of a preexisting pattern (NW-SE and NE-SW) with the directions of the grids observed on most bodies observed so far in the solar system would indicate that this local pattern is a part of a global calistean grid.

Thomas, P. G.↗

Accounting for the Effects of Surface BRDF on Satellite Cloud and Trace-Gas Retrievals: A New Approach Based on Geometry-Dependent Lambertian-Equivalent Reflectivity Applied to OMI Algorithms

Most satellite nadir ultraviolet and visible cloud, aerosol, and trace-gas algorithms make use of climatological surface reflectivity databases. For example, cloud and NO2 retrievals for the Ozone Monitoring Instrument (OMI) use monthly gridded surface reflectivity climatologies that do not depend upon the observation geometry. In reality, reflection of incoming direct and diffuse solar light from land or ocean surfaces is sensitive to the sun-sensor geometry. This dependence is described by the bidirectional reflectance distribution function (BRDF). To account for the BRDF, we propose to use a new concept of geometry-dependent Lambertian equivalent reflectivity (LER). Implementation within the existing OMI cloud and NO2 retrieval infrastructure requires changes only to the input surface reflectivity database. The geometry-dependent LER is calculated using a vector radiative transfer model with high spatial resolution BRDF information from the Moderate Resolution Imaging Spectroradiometer (MODIS) over land and the Cox-Munk slope distribution over ocean with a contribution from water-leaving radiance. We compare the geometry-dependent and climatological LERs for two wavelengths, 354 and 466 nm, that are used in OMI cloud algorithms to derive cloud fractions. A detailed comparison of the cloud fractions and pressures derived with climatological and geometry-dependent LERs is carried out. Geometry-dependent LER and corresponding retrieved cloud products are then used as inputs to our OMI NO2 algorithm. We find that replacing the climatological OMI-based LERs with geometry-dependent LERs can increase NO2 vertical columns by up to 50% in highly polluted areas; the differences include both BRDF effects and biases between the MODIS and OMI-based surface reflectance data sets. Only minor changes to NO2 columns (within 5 %) are found over unpolluted and overcast areas.

Lambert-equivalent reflectivity↗

Radar Altimetry as a Proxy for Determining Terrestrial Water Storage Variability in Tropical Basins

The Gravity Recovery and Climate Experiment (GRACE) mission has provided us with unforeseen information on terrestrial water-storage (TWS) variability, contributing to our understanding of global hydrological processes, including hydrological extreme events and anthropogenic impacts on water storage. Attempts to decompose GRACE-based TWS signals into its different water storage layers, i.e., surface water storage (SWS), soil moisture, groundwater and snow, have shown that SWS is a principal component, particularly in the tropics, where major rivers flow over arid regions at high latitudes. Here, we demonstrate that water levels, measured with radar altimeters at a limited number of locations, can be used to reconstruct gridded GRACE-based TWS signals in the Amazon basin, at spatial resolutions ranging from 0.5 to 3°, with mean absolute errors (MAE) as low as 2.5 cm and correlations as high as 0.98. We show that, at 3° spatial resolution, spatially-distributed TWS time series can be precisely reconstructed with as few as 41 water-level time series located within the basin. The proposed approach is competitive when compared to existing TWS estimates derived from physically based and computationally expensive methods. Also, a validation experiment indicates that TWS estimates can be extrapolated to periods beyond that of the model regression with low errors. The approach is robust, based on regression models and interpolation techniques, and offers a new possibility to reproduce spatially and temporally distributed TWS that could be used to fill inter-mission gaps and to extend GRACE-based TWS time series beyond its timespan.

Terrestrial water storage↗

Visual SLAM Using Variance Grid Maps

An algorithm denoted Gamma-SLAM performs further processing, in real time, of preprocessed digitized images acquired by a stereoscopic pair of electronic cameras aboard an off-road robotic ground vehicle to build accurate maps of the terrain and determine the location of the vehicle with respect to the maps. Part of the name of the algorithm reflects the fact that the process of building the maps and determining the location with respect to them is denoted simultaneous localization and mapping (SLAM). Most prior real-time SLAM algorithms have been limited in applicability to (1) systems equipped with scanning laser range finders as the primary sensors in (2) indoor environments (or relatively simply structured outdoor environments). The few prior vision-based SLAM algorithms have been feature-based and not suitable for real-time applications and, hence, not suitable for autonomous navigation on irregularly structured terrain. The Gamma-SLAM algorithm incorporates two key innovations: Visual odometry (in contradistinction to wheel odometry) is used to estimate the motion of the vehicle. An elevation variance map (in contradistinction to an occupancy or an elevation map) is used to represent the terrain. The Gamma-SLAM algorithm makes use of a Rao-Blackwellized particle filter (RBPF) from Bayesian estimation theory for maintaining a distribution over poses and maps. The core idea of the RBPF approach is that the SLAM problem can be factored into two parts: (1) finding the distribution over robot trajectories, and (2) finding the map conditioned on any given trajectory. The factorization involves the use of a particle filter in which each particle encodes both a possible trajectory and a map conditioned on that trajectory. The base estimate of the trajectory is derived from visual odometry, and the map conditioned on that trajectory is a Cartesian grid of elevation variances. In comparison with traditional occupancy or elevation grid maps, the grid elevation variance maps are much better for representing the structure of vegetated or rocky terrain.

Howard, Andrew B.↗

Characterization of Errors in Satellite-Based HCHO/NO 2 Tropospheric Column Ratios With Respect to Chemistry, Column-to-PBL Translation, Spatial Representation, and Retrieval Uncertainties

The availability of formaldehyde (HCHO) (a proxy for volatile organic compound reactivity) and nitrogen dioxide (NO 2 ) (a proxy for nitrogen oxides) tropospheric columns from ultraviolet–visible (UV–Vis) satellites has motivated many to use their ratios to gain some insights into the near-surface ozone sensitivity. Strong emphasis has been placed on the challenges that come with transforming what is being observed in the tropospheric column to what is actually in the planetary boundary layer (PBL) and near the surface; however, little attention has been paid to other sources of error such as chemistry, spatial representation, and retrieval uncertainties. Here we leverage a wide spectrum of tools and data to quantify those errors carefully. Concerning the chemistry error, a well-characterized box model constrained by more than 500 h of aircraft data from NASA's air quality campaigns is used to simulate the ratio of the chemical loss of HO 2 + RO 2 (LRO x ) to the chemical loss of NO x (LNO x ). Subsequently, we challenge the predictive power of HCHO/NO 2 ratios (FNRs), which are commonly applied in current research, in detecting the underlying ozone regimes by comparing them to LRO x /LNO x . FNRs show a strongly linear (R 2 =0.94) relationship with LRO x /LNO x , but only on the logarithmic scale. Following the baseline (i.e., ln(LRO x /LNO x ) = −1.0 ± 0.2) with the model and mechanism (CB06, r2) used for segregating NO x -sensitive from VOC-sensitive regimes, we observe a broad range of FNR thresholds ranging from 1 to 4. The transitioning ratios strictly follow a Gaussian distribution with a mean and standard deviation of 1.8 and 0.4, respectively. This implies that the FNR has an inherent 20 % standard error (1σ) resulting from not accurately describing the RO x –HO x cycle. We calculate high ozone production rates (PO 3 ) dominated by large HCHO × NO 2 concentration levels, a new proxy for the abundance of ozone precursors. The relationship between PO 3 and HCHO × NO 2 becomes more pronounced when moving towards NO x -sensitive regions due to nonlinear chemistry; our results indicate that there is fruitful information in the HCHO × NO 2 metric that has not been utilized in ozone studies. The vast amount of vertical information on HCHO and NO 2 concentrations from the air quality campaigns enables us to parameterize the vertical shapes of FNRs using a second-order rational function permitting an analytical solution for an altitude adjustment factor to partition the tropospheric columns into the PBL region. We propose a mathematical solution to the spatial representation error based on modeling isotropic semivariograms. Based on summertime-averaged data, the Ozone Monitoring Instrument (OMI) loses 12 % of its spatial information at its native resolution with respect to a high-resolution sensor like the TROPOspheric Monitoring Instrument (TROPOMI) (> 5.5 × 3.5 km 2 ). A pixel with a grid size of 216 km 2 fails at capturing ∼ 65 % of the spatial information in FNRs at a 50 km length scale comparable to the size of a large urban center (e.g., Los Angeles). We ultimately leverage a large suite of in situ and ground-based remote sensing measurements to draw the error distributions of daily TROPOMI and OMI tropospheric NO 2 and HCHO columns. At a 68 % confidence interval (1σ), errors pertaining to daily TROPOMI observations, either HCHO or tropospheric NO 2 columns, should be above 1.2–1.5 × 10 16 molec. cm −2 to attain a 20 %–30 % standard error in the ratio. This level of error is almost non-achievable with the OMI given its large error in HCHO. The satellite column retrieval error is the largest contributor to the total error (40 %–90 %) in the FNRs. Due to a stronger signal in cities, the total relative error (< 50 %) tends to be mild, whereas areas with low vegetation and anthropogenic sources (e.g., the Rocky Mountains) are markedly uncertain (> 100 %). Our study suggests that continuing development in the retrieval algorithm and sensor design and calibration is essential to be able to advance the application of FNRs beyond a qualitative metric.

Amir H. Souri↗

Investigation of the effects of external current systems on the MAGSAT data utilizing grid cell modeling techniques

The feasibility of modeling magnetic fields due to certain electrical currents flowing in the Earth's ionosphere and magnetosphere was studied. Initial efforts were devoted to reading MAGSAT data tapes in preparation for further analysis of the MAGSAT data. Further efforts concern a modeling procedure developed to compute the magnetic field at satellite orbit due to hypothesized current distributions in the ionosphere and magnetosphere. This technique utilizes a linear current element representation of the large scale space current system.

Klumpar, D. M.↗