Search NASA⌕ Search

SEARCH · Search NASA

Results for “grid computing”

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 595 records · Page 33

SUPIN: A Tool for the Aerodynamic Design and Analysis of Supersonic Inlets

The Supersonic Inlet Design and Analysis Tool (SUPIN) is a computational program that performs geometric modeling and aerodynamic design and analysis of a set of supersonic inlets whose types include the axisymmetric pitot, three-dimensional pitot, axisymmetric, two-dimensional, two-dimensional bifurcated-duct, and streamline-traced inlets. The inlets are modeled by dividing the inlet flowpath into components that start at the freestream and ends at the engine face. The primary components of the inlet include the external supersonic diffuser, throat section, and subsonic diffuser. Each component is characterized by a set of inputs that include geometric and aerodynamic factors. Examples of the geometric factors are angles, lengths, area ratios, and coordinates that describe the geometry of the components. Examples of the aerodynamic factors are Mach numbers, flow ratios, and pressure ratios that set constraints or desired aerodynamic properties of the component. These input factors are specified within a text-based input data file that is read into SUPIN. SUPIN is a Fortran 95 program that performs the inlet design and analysis as a batch or non-interactive process. SUPIN writes output data files describing the inlet geometry and aerodynamic performance. SUPIN uses analytic, empirical, and numerical methods to perform the inlet design and analysis. The inlet geometry is constructed using planar and three-dimensional geometry constructs and is represented by coordinates, angles, areas, profiles, and surfaces. Surfaces are output in the form of a Plot3D surface grid file. SUPIN does have some capability to automatically generate planar and three-dimensional, multi-block, structured grids for computational fluid dynamics (CFD) analysis. The aerodynamic analysis performed by SUPIN computes quasi-one-dimensional properties, such as Mach numbers, pressures, and temperatures, through the inlet flowpath. Planar aerodynamic properties, such as those for shock and Mach waves, are computed for the supersonic compression of the external supersonic diffuser and throat section. The inlet performance is characterized by the inlet flow rates, total pressure recovery, and inlet drag. The primary intent of this document is to describe the usage of SUPIN. The methods used within SUPIN are documented and may serve as a useful reference for understanding the aerodynamic design and analysis of supersonic inlets.

supersonic inlets↗

WTK-LED: The WIND Toolkit Long-Term Ensemble Dataset

To satisfy a wide group of stakeholders across various wind energy disciplines, including but not limited to stakeholders in the distributed and utility scale wind industry, the new emerging airborne wind energy field, grid integration, power systems modeling, environmental modeling, and researchers in academia, and to close some of the gaps that current public datasets have, we aimed at developing an updated version of the meteorological WIND Toolkit, named WIND Toolkit Long-term Ensemble Dataset (WTK-LED), which is a meteorological dataset providing time series every 5 min and 2 km, including model uncertainty of wind speed at every modeling grid point so that users are provided with a range of possible wind speeds every 2 km. The data were produced using the Weather Research and Forecasting Model (WRF). The vertical grid used in WTK-LED includes many vertical layers in the atmospheric boundary layer to provide information of atmospheric quantities across the rotor layer of utility scale and distributed wind turbines. The WTK-LED includes: 1) Numerical simulations covering the continental United States, Alaska, and Hawaii, with high-resolution data being available for 3 years (2018-2020). 2) Climate simulations from Argonne National Laboratories covering the North American continent, including Alaska, Canada, and most of Mexico and the Caribbean Islands. These simulations complement the new WTK-LED to offer a 4-km dataset covering 20 years, from 2001-2020. 3) Specific long-term,high-resolution offshore simulations have been conducted separately for the US coasts, Hawaii, and the Great Lakes, leading to the 2023 National Offshore Wind data set. This report focuses on a description of the land-based WTK-LED for CONUS, Hawaii, and Alaska, for the 3-year 2-km/5-min dataset and the 20-year 4-km/hourly dataset, as well as the uncertainty quantification method. We also provide limited validation results. Based on our results to date, we suggest use cases and applications for each dataset of the WTK-LED.

17 WIND ENERGY↗

Adjoint-Based, Three-Dimensional Error Prediction and Grid Adaptation

Engineering computational fluid dynamics (CFD) analysis and design applications focus on output functions (e.g., lift, drag). Errors in these output functions are generally unknown and conservatively accurate solutions may be computed. Computable error estimates can offer the possibility to minimize computational work for a prescribed error tolerance. Such an estimate can be computed by solving the flow equations and the linear adjoint problem for the functional of interest. The computational mesh can be modified to minimize the uncertainty of a computed error estimate. This robust mesh-adaptation procedure automatically terminates when the simulation is within a user specified error tolerance. This procedure for estimating and adapting to error in a functional is demonstrated for three-dimensional Euler problems. An adaptive mesh procedure that links to a Computer Aided Design (CAD) surface representation is demonstrated for wing, wing-body, and extruded high lift airfoil configurations. The error estimation and adaptation procedure yielded corrected functions that are as accurate as functions calculated on uniformly refined grids with ten times as many grid points.

Park, Michael A.↗

Subsampling of Parametric Models with Bifidelity Boosting

Least squares regression is a ubiquitous tool for building emulators (a.k.a. surrogate models) of problems across science and engineering for purposes such as design space exploration and uncertainty quantification. When the regression data are generated using an experimental design process (e.g., a quadrature grid) involving computationally expensive models, or when the data size is large, sketching techniques have shown promise at reducing the cost of the construction of the regression model while ensuring accuracy comparable to that of the full data. However, random sketching strategies, such as those based on leverage scores, lead to regression errors that are random and may exhibit large variability. To mitigate this issue, we present a novel boosting approach that leverages cheaper, lower-fidelity data of the problem at hand to identify the best sketch among a set of candidate sketches. This in turn specifies the sketch of the intended high-fidelity model and the associated data. We provide theoretical analyses of this bifidelity boosting (BFB) approach and discuss the conditions the low- and high-fidelity data must satisfy for a successful boosting. In doing so, we derive a bound on the residual norm of the BFB sketched solution relating it to its ideal, but computationally expensive, high-fidelity boosted counterpart. Finally, empirical results on both manufactured and PDE data corroborate the theoretical analyses and illustrate the efficacy of the BFB solution in reducing the regression error, as compared to the nonboosted solution.

97 MATHEMATICS AND COMPUTING↗

DER Cybersecurity Standards: Assessment and Gap Analysis

The purpose of this report is to share the comprehensive gap analysis of existing cybersecurity standards applicable to Distributed Energy Resources (DERs) within the electric power sector. This analysis aims to identify critical deficiencies in current standards, assess their alignment with industry needs, and provide actionable recommendations for enhancing cybersecurity measures. The scope encompasses various DER technologies, including solar, wind, energy storage, and hydrogen fuel cells, and emphasizes the significance of establishing robust cybersecurity frameworks and standards to safeguard these increasingly integrated systems. The report provides valuable insights for stakeholders in the DER ecosystem, including manufacturers, utilities, and regulators. It underscores the importance of continued development and refinement of cybersecurity standards to keep up with the technical advances in DERs and associated cybersecurity challenges. The analysis evaluated IEC, IEEE, ISA, ISO, and UL standards relevant to DER cybersecurity. Standards were assessed on their coverage of key requirements including data availability, integrity, confidentiality, access control, authentication, encryption, and system hardening. For each standard, the analysis assessed its alignment with current industry practices, regulatory compliance, effectiveness in addressing known risks, coverage of emerging risks, and how it promotes interoperability. The evaluation also considered potential integration challenges and barriers to adoption.

97 MATHEMATICS AND COMPUTING↗

A transformation method for deriving from a photograph, position and heading of a vehicle in a plane

Equations have been derived that transform perspectively viewed planar surface coordinates, as seen in a photograph, into coordinates of the original plane surface. These transformation equations are developed in terms of nine geometric variables that define the photographic setup and are redefined in terms of eight parameters. The parameters are then treated as independent quantities that fully characterize the transformation and are expressed directly in terms of the four corner coordinates of a reference rectangle in the object plane and their coordinates as seen in a photograph. Vehicle position is determined by transforming the perspectively viewed coordinate position of a representative vehicle target into runway coordinates. Vehicle heading is determined from the runway coordinates of two vehicle target points. When the targets are elevated above the plane of the reference grid, the computation of the heading angle is unaffected; however, the computation of the target position may require adjustment of two parameters. Methods are given for adjusting the parameters for elevation and an example is included for both nonelevated and elevated target conditions.

Sleeper, R. K.↗

Beam efflux measurements

Measurements have been made of the high energy thrust ions, (Group I), high angle/high energy ions (Group II), and high angle/low energy ions (Group IV) of a mercury electron bombardment thruster in the angular divergence range from 0 deg to greater than 90 deg. The measurements have been made as a function of thrust ion current, propellant utilization efficiency, bombardment discharge voltage, screen and accelerator grid potential (accel-decel ratio) and neutralizer keeper potential. The shape of the Group IV (charge exchange) ion plume has remained essentially fixed within the range of variation of the engine operation parameters. The magnitude of the charge exchange ion flux scales with thrust ion current, for good propellant utilization conditions. For fixed thrust ion current, charge exchange ion flux increases for diminishing propellant utilization efficiency. Facility effects influence experimental accuracies within the range of propellant utilization efficiency used in the experiments. The flux of high angle/high energy Group II ions is significantly diminished by the use of minimum decel voltages on the accelerator grid. A computer model of charge exchange ion production and motion has been developed. The program allows computation of charge exchange ion volume production rate, total production rate, and charge exchange ion trajectories for "genuine" and "facilities effects" particles. In the computed flux deposition patterns, the Group I and Group IV ion plumes exhibit a counter motion.

Komatsu, G. K.↗

Some mesh generation requirements and methods

Discretized solution algorithms, which find solutions of field equations in a two or three dimensional field, generally use meshes which are fitted to the field boundary to allow convenient formulation of boundary conditions there. A mesh is defined to be the image of a rectangular grid in computational space under a mesh mapping which maps computational space into physical space. It is not necessary that all of computational space be mapped onto the region of interest in physical space. Parts of it can be excised to give a better fit to the boundary. Many different excisions can be made to fit a single boundary; the choice depends on the mesh arrangement desired in the field.

Dickson, L. J.↗

Numerical investigation of turbulent channel flow

Fully developed turbulent channel flow was simulated numerically at Reynolds number 13800, based on centerline velocity and channel halt width. The large-scale flow field was obtained by directly integrating the filtered, three dimensional, time dependent, Navier-Stokes equations. The small-scale field motions were simulated through an eddy viscosity model. The calculations were carried out on the ILLIAC IV computer with up to 516,096 grid points. The computed flow field was used to study the statistical properties of the flow as well as its time dependent features. The agreement of the computed mean velocity profile, turbulence statistics, and detailed flow structures with experimental data is good. The resolvable portion of the statistical correlations appearing in the Reynolds stress equations are calculated. Particular attention is given to the examination of the flow structure in the vicinity of the wall.

Moin, P.↗

Numerical simulation of wall-bounded turbulent shear flows

Developments in three dimensional, time dependent numerical simulation of turbulent flows bounded by a wall are reviewed. Both direct and large eddy simulation techniques are considered within the same computational framework. The computational spatial grid requirements as dictated by the known structure of turbulent boundary layers are presented. The numerical methods currently in use are reviewed and some of the features of these algorithms, including spatial differencing and accuracy, time advancement, and data management are discussed. A selection of the results of the recent calculations of turbulent channel flow, including the effects of system rotation and transpiration on the flow are included.

Moin, P.↗

Analysis of inviscid and viscous flows in cascades with an explicit multiple-grid algorithm

A rapid technique is used for calculating inviscid and viscous flows in turbomachinery cascades. The Euler and thin-layer Navier-Stokes equations are solved using the original explicit MacCormack algorithm. The Baldwin-Lomax eddy viscosity model is used for turbulent flows. Convergence to a steady state is accelerated by use of a variable time-step and a multiple-grid scheme. Computer time is reduced through vectorization. Details of the numerical method are presented along with computed results for two low-speed wind tunnel turning vanes, a space shuttle fuel pump turbine rotor, and a supersonic inflow compressor rotor. The method can predict subtle viscous flow phenomena in cascades and is fast enough to be used as a design tool.

Chima, R. V.↗

Numerical simulation of wall-bounded turbulent shear flows

Developments in three dimensional, time dependent numerical simulation of turbulent flows bounded by a wall are reviewed. Both direct and large eddy simulation techniques are considered within the same computational framework. The computational spatial grid requirements as dictated by the known structure of turbulent boundary layers are presented. The numerical methods currently in use are reviewed and some of the features of these algorithms, including spatial differencing and accuracy, time advancement, and data management are discussed. A selection of the results of the recent calculations of turbulent channel flow, including the effects of system rotation and transpiration on the flow are included. Previously announced in STAR as N82-28577

Moin, P.↗

Analysis of inviscid and viscous flows in cascades with an explicit multiple-grid algorithm

A rapid technique is used for calculating inviscid and viscous flows in turbomachinery cascades. The Euler and thin-layer Navier-Stokes equations are solved using the original explicit MacCormack algorithm. The Baldwin-Lomax eddy viscosity model is used for turbulent flows. Convergence to a steady state is accelerated by use of a variable time-step and a multiple-grid scheme. Computer time is reduced through vectorization. Details of the numerical method are presented along with computed results for two low-speed wind tunnel turning vanes, a space shuttle fuel pump turbine rotor, and a supersonic inflow compressor rotor. The method can predict subtle viscous flow phenomena in cascades and is fast enough to be used as a design tool. Previously announced in STAR as N84-22527

Chima, R. V.↗

An adaptive finite element scheme for the Euler and Navier-Stokes equations

The implementation of two explicit finite element schemes for the solution of the compressible Euler and Navier-Stokes equations is presented. The schemes can be employed with general unstructured triangular meshes in two dimensions. Either scheme can therefore be used as the basic solver in a solution adaptive mesh procedure in a direct manner. The particular adaptive approach which is advocated here is intended for the solution of steady state problems only and involves an adaptive regeneration of the grid at prescribed stages during the false transient. The grid regeneration is accomplished by a mesh generator which has the capability of generating triangular grids over computational domains of arbitrary shape. The procedure is illustrated by solving transonic flows over multi-airfoil configurations and high speed flows, involving shock interactions, past circular cylinders.

Morgan, K.↗

Climatology of the stratospheric polar vortex and planetary wave breaking

The distribution of Ertel's potential vorticity (PV) on the 850 K isentropic surface is used to establish a climatology for the transient evolution of the planetary scale circulation in the Northern Hemisphere winter midstratosphere. PV distributions are computed from gridded NMC daily temperature and height maps for the 10 and 30 mb levels, and show that a very good approximation for 850 K PV can be derived from 10 mb heights and temperatures alone. It is assumed that reversals of the latitudinal gradient of PV, localized in longitude and latitude may be regarded as signatures of planetary wave breaking. Wave breaking identified by such signatures tends to occur mainly in the vicinity of the Aleutian anticyclone, with a secondary maximum over Europe. The area of the polar vortex, defined as the area enclosed by PV contours greater than a certain critical value, is strongly influenced by wave breaking. Erosion of the polar vortex due to transport and mixing of PV leads to a preconditioned state, when defined in terms of vortex area, that always occurs prior to major stratospheric warmings. During winters with little PV transport or mixing, the vortex area evolves rather uniformly in response to radiative forcing. During winters with major sudden warmings, the wave breaking signature as defined here first appears at low values of PV, then rapidly moves toward higher values as the vortex area is reduced and the 'surf-zone' structure becomes well defined.

Baldwin, Mark P.↗