SEARCH · Engineering Papers

Results for “fixed 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 217 records · Page 12

All sources NASA NTRS DOE OSTI

NASA NTRS1977

NASTRAN hydroelastic modal studies. Volume 3: NASTRAN 3-D hydroelastic analysis and MESHGEN user's manuals

The method used to formulate the fluid/structure equations is described in the updates to the NASTRAN Theoretical Manual. The basis for defining the fluid are three dimensional finite elements connected to fluid grid points defining the Eulerian pressure at a point fixed in space. The use of a pressure single degree of freedom at each point rather than three displacements allows a finer mesh of elements with a reasonable matrix order. In the formulation of the fluid/structure system the interior fluid degrees of freedom are transformed and removed from the solution matrices. The eigenvalues of the combination are extracted from small, fully dense, symmetric mass and stiffness matrices, efficiently processed with the Givens method. The solution matrices are defined by the free surface displacements and the reduced structure coordinates.

Source record↗

NASA NTRS2001

Mars Constellation Design

In this study, several different Mars constellation designs were evaluated. Various figures of merit were used to examine constellation configurations in order to optimize the performance of Mars Network. Plots were also obtained in order to make a clear comparison of the various figures of merit. In the following analysis, six figures of merit were used in the comparison (average time to achieve a 1 meter fix, maximum time to achieve a 1 meter fix, maximum wait time between satellite sightings, percent of time a satellite is in view, percent of cases to achieve a 1 meter fix in 10 minutes, and percent of cases to achieve a 1 meter fix in 1 hour). For the communication figures of merit (maximum wait time between satellite sightings and percent of time a satellite is in view), an initial run over a period of 10 days was used to produce the data. For the remaining figures of merit, the simulation was first restarted at least 100 times over a period of 1000 Earth days. Each time the simulation was continued until a 1 meter fix was achieved at all stations. The ground station grid was spaced 20 deg in longitude and 5 deg in latitude. In order to obtain the following graphs the values were first averaged over longitude and then an average of the positive and negative latitude values was made (i.e., -90 deg lat and 90 deg lat were averaged). In each case, the resulting figures of merit are shown. For some cases, the difference between the given case and the 'Eagles 2 Satellites' case is also given.

Source record↗

NASA NTRS2006

Calculation of Thermally-Induced Displacements in Spherically Domed Ion Engine Grids

An analytical method for predicting the thermally-induced normal and tangential displacements of spherically domed ion optics grids under an axisymmetric thermal loading is presented. A fixed edge support that could be thermally expanded is used for this analysis. Equations for the displacements both normal and tangential to the surface of the spherical shell are derived. A simplified equation for the displacement at the center of the spherical dome is also derived. The effects of plate perforation on displacements and stresses are determined by modeling the perforated plate as an equivalent solid plate with modified, or effective, material properties. Analytical model results are compared to the results from a finite element model. For the solid shell, comparisons showed that the analytical model produces results that closely match the finite element model results. The simplified equation for the normal displacement of the spherical dome center is also found to accurately predict this displacement. For the perforated shells, the analytical solution and simplified equation produce accurate results for materials with low thermal expansion coefficients.

Soulas, George C.↗

DOE OSTI2015

Effect of kinetic energy on the doping efficiency of cesium cations into superfluid helium droplets

We present an experimental investigation of the effect of kinetic energy on the ion doping efficiency of superfluid helium droplets using cesium cations from a thermionic emission source. The kinetic energy of Cs{sup +} is controlled by the bias voltage of a collection grid collinearly arranged with the droplet beam. Efficient doping from ions with kinetic energies from 20 eV up to 480 V has been observed in different sized helium droplets. The relative ion doping efficiency is determined by both the kinetic energy of the ions and the average size of the droplet beam. At a fixed source temperature, the number of doped droplets increases with increasing grid voltage, while the relative ion doping efficiency decreases. This result implies that not all ions are captured upon encountering with a sufficiently large droplet, a deviation from the near unity doping efficiency for closed shell neutral molecules. We propose that this drop in ion doping efficiency with kinetic energy is related to the limited deceleration rate inside a helium droplet. When the source temperature changes from 14 K to 17 K, the relative ion doping efficiency decreases rapidly, perhaps due to the lack of viable sized droplets. The size distribution of the Cs{sup +}-doped droplet beam can be measured by deflection and by energy filtering. The observed doped droplet size is about 5 × 10{sup 6} helium atoms when the source temperature is between 14 K and 17 K.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2013

Enhancing the Representation of Subgrid Land Surface Characteristics in Land Surface Models

Land surface heterogeneity has long been recognized and increasingly incorporated in the land surface modelling. In most existing land surface models, the spatial variability of surface cover is represented as subgrid composition of multiple surface cover types. In this study, we developed a new subgrid classification method (SGC) that accounts for the topographic variability of the vegetation cover. Each model grid cell was represented with a number of elevation classes and each elevation class was further described by a number of vegetation types. The numbers of elevation classes and vegetation types were variable and optimized for each model grid so that the spatial variability of both elevation and vegetation can be reasonably explained given a pre-determined total number of classes. The subgrid structure of the Community Land Model (CLM) was used as an example to illustrate the newly developed method in this study. With similar computational burden as the current subgrid vegetation representation in CLM, the new method is able to explain at least 80% of the total subgrid PFTs and greatly reduced the variations of elevation within each subgrid class compared to the baseline method where a single elevation class is assigned to each subgrid PFT. The new method was also evaluated against two other subgrid methods (SGC1 and SGC2) that assigned fixed numbers of elevation and vegetation classes for each model grid with different perspectives of surface cover classification. Implemented at five model resolutions (0.1°, 0.25°, 0.5°, 1.0° and 2.0°) with three maximum-allowed total number of classes N_class of 24, 18 and 12 representing different computational burdens over the North America (NA) continent, the new method showed variable performances compared to the SGC1 and SGC2 methods. However, the advantage of the SGC method over the other two methods clearly emerged at coarser model resolutions and with moderate computational intensity (N_class = 18) as it explained the most PFTs and elevation variability among the three subgrid methods. Spatially, the SGC method explained more elevation variability in topography-complex areas and more vegetation variability in flat areas. Furthermore, the variability of both elevation and vegetation explained by the new method was more spatially homogeneous regardless of the model resolutions and computational burdens. The SGC method will be implemented in CLM over the NA continent to assess its impacts on simulating land surface processes.

characteristics↗

NASA NTRS2014

A New Hybrid Spatio-temporal Model for Estimating Daily Multi-year PM2.5 Concentrations Across Northeastern USA Using High Resolution Aerosol Optical Depth Data

The use of satellite-based aerosol optical depth (AOD) to estimate fine particulate matter PM(sub 2.5) for epidemiology studies has increased substantially over the past few years. These recent studies often report moderate predictive power, which can generate downward bias in effect estimates. In addition, AOD measurements have only moderate spatial resolution, and have substantial missing data. We make use of recent advances in MODIS satellite data processing algorithms (Multi-Angle Implementation of Atmospheric Correction (MAIAC), which allow us to use 1 km (versus currently available 10 km) resolution AOD data.We developed and cross validated models to predict daily PM(sub 2.5) at a 1X 1 km resolution across the northeastern USA (New England, New York and New Jersey) for the years 2003-2011, allowing us to better differentiate daily and long term exposure between urban, suburban, and rural areas. Additionally, we developed an approach that allows us to generate daily high-resolution 200 m localized predictions representing deviations from the area 1 X 1 km grid predictions. We used mixed models regressing PM(sub 2.5) measurements against day-specific random intercepts, and fixed and random AOD and temperature slopes. We then use generalized additive mixed models with spatial smoothing to generate grid cell predictions when AOD was missing. Finally, to get 200 m localized predictions, we regressed the residuals from the final model for each monitor against the local spatial and temporal variables at each monitoring site. Our model performance was excellent (mean out-of-sample R(sup 2) = 0.88). The spatial and temporal components of the out-of-sample results also presented very good fits to the withheld data (R(sup 2) = 0.87, R(sup)2 = 0.87). In addition, our results revealed very little bias in the predicted concentrations (Slope of predictions versus withheld observations = 0.99). Our daily model results show high predictive accuracy at high spatial resolutions and will be useful in reconstructing exposure histories for epidemiological studies across this region.

Air pollution↗

NASA NTRS2004

Computational Aerothermodynamic Simulation Issues on Unstructured Grids

The synthesis of physical models for gas chemistry and turbulence from the structured grid codes LAURA and VULCAN into the unstructured grid code FUN3D is described. A directionally Symmetric, Total Variation Diminishing (STVD) algorithm and an entropy fix (eigenvalue limiter) keyed to local cell Reynolds number are introduced to improve solution quality for hypersonic aeroheating applications. A simple grid-adaptation procedure is incorporated within the flow solver. Simulations of flow over an ellipsoid (perfect gas, inviscid), Shuttle Orbiter (viscous, chemical nonequilibrium) and comparisons to the structured grid solvers LAURA (cylinder, Shuttle Orbiter) and VULCAN (flat plate) are presented to show current capabilities. The quality of heating in 3D stagnation regions is very sensitive to algorithm options in general, high aspect ratio tetrahedral elements complicate the simulation of high Reynolds number, viscous flow as compared to locally structured meshes aligned with the flow.

Gnoffo, Peter A.↗

DOE OSTI2020

Compare linear-system solver and preconditioner stacks with emphasis on GPU performance and propose phase-2 NGP solver development pathway

The goal of the ExaWind project is to enable predictive simulations of wind farms comprised of many megawatt-scale turbines situated in complex terrain. Predictive simulations will require computational fluid dynamics (CFD) simulations for which the mesh resolves the geometry of the turbines and captures the rotation and large deflections of blades. Whereas such simulations for a single turbine are arguably petascale class, multi-turbine wind farm simulations will require exascale-class resources. The primary physics codes in the ExaWind project are Nalu-Wind, which is an unstructured-grid solver for the acoustically incompressible Navier-Stokes equations, and OpenFAST, which is a whole-turbine simulation code. The Nalu-Wind model consists of the mass-continuity Poisson-type equation for pressure and a momentum equation for the velocity. For such modeling approaches, simulation times are dominated by linear-system setup and solution for the continuity and momentum systems. For the ExaWind challenge problem, the moving meshes greatly affect overall solver costs as reinitialization of matrices and recomputation of preconditioners is required at every time step. In this report we evaluated GPU-performance baselines for the linear solvers in the Trilinos and hypre solver stacks using two representative Nalu-Wind simulations: an atmospheric boundary layer precursor simulation on a structured mesh, and a fixed-wing simulation using unstructured overset meshes. Both strong-scaling and weak-scaling experiments were conducted on the OLCF supercomputer Summit and similar proxy clusters. We focused on the performance of multi-threaded Gauss-Seidel and two-stage Gauss-Seidel that are extensions of classical Gauss-Seidel; of one-reduce GMRES, a communication-reducing variant of the Krylov GMRES; and algebraic multigrid methods that incorporate the afore-mentioned methods. The team has established that AMG methods are capable of solving linear systems arising from the fixed-wing overset meshes on CPU, a critical intermediate result for ExaWind FY20 Q3 and Q4 milestones. For the fixed-wing strong-scaling study (model with 3M grid-points), the team identified that Nalu-Wind simulations with the new Trilinos and hypre solvers scale to modest GPU counts, maintaining above 70% efficiency up to 6 GPUs. However, there still remain significant bottlenecks to performance: matrix assembly (hypre), AMG setup (hypre and Trilinos) In the weak-scaling experiments (going from 0.4M to 211M gridpoints), it's shown that the solver apply phases are faster on GPUs, but that Nalu-Wind simulation times grow, primarily due to the multigrid-setup process. Finally, based on the report outcomes, we propose a linear solver path-forward for the remainder of the ExaWind project. Near term, the NREL team will continue their work on GPU-based linear-system assembly. They will also investigate how the use of alternatives to the NVIDIA UVM (unified virtual memory) paradigm affects performance. Longer term, the NREL team will evaluate algorithmic performance on other types of accelerators and merge their improvements back to the main hypre repository branch. Near term, the Trilinos team will address performance bottlenecks identified in this milestone, such as implementing a GPU-based segregated momentum solve and reusing matrix graphs across linear-system assembly phases. Longer term, the Trilinos team will do detailed analysis and optimization of multigrid setup.

17 WIND ENERGY↗

DOE OSTI2017

An Investigation into Solution Verification for CFD-DEM

This report presents the study of the convergence behavior of the computational fluid dynamicsdiscrete element method (CFD-DEM) method, specifically National Energy Technology Laboratory’s (NETL) open source MFiX code (MFiX-DEM) with a diffusion based particle-tocontinuum filtering scheme. In particular, this study focused on determining if the numerical method had a solution in the high-resolution limit where the grid size is smaller than the particle size. To address this uncertainty, fixed particle beds of two primary configurations were studied: i) fictitious beds where the particles are seeded with a random particle generator, and ii) instantaneous snapshots from a transient simulation of an experimentally relevant problem. Both problems considered a uniform inlet boundary and a pressure outflow. The CFD grid was refined from a few particle diameters down to 1/6 th of a particle diameter. The pressure drop between two vertical elevations, averaged across the bed cross-section was considered as the system response quantity of interest. A least-squares regression method was used to extrapolate the grid-dependent results to an approximate “grid-free” solution in the limit of infinite resolution. The results show that the diffusion based scheme does yield a converging solution. However, the convergence is more complicated than encountered in simpler, single-phase flow problems showing strong oscillations and, at times, oscillations superimposed on top of globally non-monotonic behavior. The challenging convergence behavior highlights the importance of using at least four grid resolutions in solution verification problems so that (over-determined) regression-based extrapolation methods may be applied to approximate the grid-free solution. The grid-free solution is very important in solution verification and VVUQ exercise in general as the difference between it and the reference solution largely determines the numerical uncertainty. By testing different randomized particle configurations of the same general problem (for the fictitious case) or different instances of freezing a transient simulation, the numerical uncertainties appeared to be on the same order of magnitude as ensemble or time averaging uncertainties. By testing different drag laws, almost all cases studied show that model form uncertainty in this one, very important closure relation was larger than the numerical uncertainty, at least with a reasonable CFD grid, roughly five particle diameters. In this study, the diffusion width (filtering length scale) was mostly set at a constant of six particle diameters. A few exploratory tests were performed to show that similar convergence behavior was observed for diffusion widths greater than approximately two particle diameters. However, this subject was not investigated in great detail because determining an appropriate filter size is really a validation question which must be determined by comparison to experimental or highly accurate numerical data. Future studies are being considered targeting solution verification of transient simulations as well as validation of the filter size with direct numerical simulation data.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Optimal electric-distribution-grid planning considering the demand-side flexibility of thermal building systems for a test case in Singapore

The planning of district-scale electric grids, i.e., distribution grids, has traditionally relied on finding the most cost-effective design such that they are able to supply the peak loads in a district. With the advent of electric demand side flexibility (DSF), there is the opportunity to reshape peak loads such that the investment cost of the electric grid decreases in exchange for a minor increase in the operation cost. This paper formulates an optimal planning approach for the electric grid at the district scale, which incorporates the DSF from thermal building systems, e.g., heating ventilation and air-conditioning (HVAC) systems. The problem is formulated as a mixed-integer linear program (MILP) and aims at minimizing the investment cost for the grid along with the operation cost of the flexible loads. This is subjected to the fixed electricity demand and thermal comfort constraints of building occupants. To this end, linear models for the thermal comfort in the buildings and the power flow in electric grid are considered. The approach is tested on a district planning test case based in Singapore, where the results show up to 30.9 % reductions in investment cost and up to 3.7 % reduction in total annualized cost. Urban planning authorities, developers and utility companies can all benefit from the presented approach to make optimized investment decisions. For building operators, the results point to the need of adopting their control systems for DSF.

Troitzsch, Sebastian↗

NASA NTRS1997

Correction of Single Frequency Altimeter Measurements for Ionosphere Delay

This study is a preliminary analysis of the accuracy of various ionosphere models to correct single frequency altimeter height measurements for Ionospheric path delay. In particular, research focused on adjusting empirical and parameterized ionosphere models in the parameterized real-time ionospheric specification model (PRISM) 1.2 using total electron content (TEC) data from the global positioning system (GPS). The types of GPS data used to adjust PRISM included GPS line-of-sight (LOS) TEC data mapped to the vertical, and a grid of GPS derived TEC data in a sun-fixed longitude frame. The adjusted PRISM TEC values, as well as predictions by IRI-90, a climatotogical model, were compared to TOPEX/Poseidon (T/P) TEC measurements from the dual-frequency altimeter for a number of T/P tracks. When adjusted with GPS LOS data, the PRISM empirical model predicted TEC over 24 1 h data sets for a given local time to with in a global error of 8.60 TECU rms during a midnight centered ionosphere and 9.74 TECU rms during a noon centered ionosphere. Using GPS derived sun-fixed TEC data, the PRISM parameterized model predicted TEC within an error of 8.47 TECU rms centered at midnight and 12.83 TECU rms centered at noon. From these best results, it is clear that the proposed requirement of 3-4 TECU global rms for TOPEX/Poseidon Follow-On will be very difficult to meet, even with a substantial increase in the number of GPS ground stations, with any realizable combination of the aforementioned models or data assimilation schemes.

Schreiner, William S.↗

DOE OSTI2018

An Assessment of High-Resolution Gridded Temperature Datasets over California

Abstract High-resolution gridded datasets are in high demand because they are spatially complete and include important finescale details. Previous assessments have been limited to two to three gridded datasets or analyzed the datasets only at the station locations. Here, eight high-resolution gridded temperature datasets are assessed two ways: at the stations, by comparing with Global Historical Climatology Network–Daily data; and away from the stations, using physical principles. This assessment includes six station-based datasets, one interpolated reanalysis, and one dynamically downscaled reanalysis. California is used as a test domain because of its complex terrain and coastlines, features known to differentiate gridded datasets. As expected, climatologies of station-based datasets agree closely with station data. However, away from stations, spread in climatologies can exceed 6°C. Some station-based datasets are very likely biased near the coast and in complex terrain, due to inaccurate lapse rates. Many station-based datasets have large unphysical trends (>1°C decade−1) due to unhomogenized or missing station data—an issue that has been fixed in some datasets by using homogenization algorithms. Meanwhile, reanalysis-based gridded datasets have systematic biases relative to station data. Dynamically downscaled reanalysis has smaller biases than interpolated reanalysis, and has more realistic variability and trends. Dynamical downscaling also captures snow–albedo feedback, which station-based datasets miss. Overall, these results indicate that 1) gridded dataset choice can be a substantial source of uncertainty, and 2) some datasets are better suited for certain applications.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Modeling simultaneous particle shrinkage, dissolution and breakage using the modified moving grid technique

Simultaneous shrinkage, dissolution and breakage are important particle size reduction phenomena that characterize processes like the reactive degradation of solid chemicals. The dynamics of the particle size distribution (PSD) for such processes are non-trivial to model due to the number expending processes brought about by the eventual dissolution of particles. To this end, Population Balance Model (PBM) resolved through the sectional techniques is the natural approach. Here, we introduce a modified Moving Grid technique (m-MGT) to accurately resolve the particle size reduction phenomena. Our technique mimics the perpetual particle shrinkage through a continuously left-moving size grid and incorporates a strategic grid removal routine to capture the disappearance of particles. Coupled with the Fixed Pivot (FP) discretization for breakage, our m-MGT not only preserves the moment-related properties, but also benchmarked very well against the analytical number densities and exhibited a minimum of first-order convergence in all assessed case studies.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Simulink Modeling and Dynamic Study of Fixed-Speed, Variable-Speed, and Ternary Pumped Storage Hydropower

Pumped Storage Hydropower (PSH) is one of the most popular energy storage technologies in the world. It uses an upper reservoir to store water which can be later used during high-demand. In the United States, most of the energy storage capability actually corresponds to PSH. Moreover, PSH also brings multiple benefits to grid operation. This report presents the Simulink models of three common PSH technologies: Fixed-Speed (FS), Variable-Speed (VS), and Ternary (T)-PSH. These models are available to the general public on this GitHub repository, which contains the MATLAB model initialization files, the Simulink model files, and supplementary MATLAB code used to obtain the figures in this work. For each PSH model, an introductory description of the model components and other relevant functionalities are provided. For further information regarding the models and the initialization parameters, the reader is referred to the shared files in the repository. This report also presents the dynamic behavior of each model. The response of such models to a load event is analyzed and matched with each model's features. A custom IEEE 39 bus case is employed for the FS and T-PSH simulations, while the VS-PSH is simulated on a simplified three-bus test system due to the computational complexity of the model. For the T-PSH, the steady-state and the switching between several operating modes are also studied in this work.

13 HYDRO ENERGY↗

DOE OSTI2013

Annual Fossil-Fuel CO2 Emissions: Isomass of Emissions Gridded by One Degree latitude by One Degree Longitude (V. 2013)

The monthly, isotopic (δ 13C) fossil-fuel CO2 emissions estimates from 1950-2013 provided in this database are derived from time series of global, regional, and national fossil-fuel CO2 emissions (Boden et al. 2016), the references therein, and the methodology described in Andres et al. (2011). The data accessible here take these tabular, national, mass-emissions data, multiply them by stable carbon isotopic signatures (δ 13C) as described in Andres et al. (2000), and distribute them spatially on a one degree latitude by one degree longitude grid. The within-country spatial distribution is achieved through a fixed population distribution as reported in Andres et al. (1996). Note that the mass-emissions data used here are based on fossil-fuel consumption estimates as these are more representative of within country emissions than fossil-fuel production estimates (see http://cdiac.ess-dive.lbl.gov/faq.html#Q10 for a description why emission totals based upon consumption differ from those based upon production).

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2015

Annual Fossil-Fuel CO2 Emissions: Isomass of Emissions Gridded by One Degree Latitude by One Degree Longitude (V. 2015)

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2013

Monthly Fossil-Fuel CO2 Emissions: Isomass of Emissions Gridded by One Degree Latitude by One Degree Longitude (1950 - 2010) (V. 2013)

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2015

Monthly Fossil-Fuel CO2 Emissions: Isomass of Emissions Gridded by One Degree Latitude by One Degree Longitude

The basic data provided in these data files are derived from time series of Global, Regional, and National Fossil-Fuel CO2 Emissions (http://cdiac.ess-dive.lbl.gov/trends/emis/overview_2011.html), the references therein, and the methodology described in Andres et al. (2011). The data accessible here take these tabular, national, mass-emissions data, multiply them by stable carbon isotopic signature (del 13C) as described in Andres et al. (2000), and distribute them spatially on a one degree latitude by one degree longitude grid. The within-country spatial distribution is achieved through a fixed population distribution as reported in Andres et al. (1996). Note that the mass-emissions data used here are based on fossil-fuel consumption estimates as these are more representative of within country emissions than fossil-fuel production estimates (see http://cdiac.ess-dive.lbl.gov/faq.html#Q10 for a description why emission totals based upon consumption differ from those based upon production).

54 ENVIRONMENTAL SCIENCES↗