Engineering Papers⌕ Search

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.

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

At least 127 records · Page 7

Interpolated Sounding and Gridded Sounding Value-Added Products

Standard Atmospheric Radiation Measurement (ARM) Climate Research Facility sounding files provide atmospheric state data in one dimension of increasing time and height per sonde launch. Many applications require a quick estimate of the atmospheric state at higher time resolution. The INTERPOLATEDSONDE (i.e., Interpolated Sounding) Value-Added Product (VAP) transforms sounding data into continuous daily files on a fixed time-height grid, at 1-minute time resolution, on 332 levels, from the surface up to a limit of approximately 40 km. The grid extends that high so the full height of soundings can be captured; however, most soundings terminate at an altitude between 25 and 30 km, above which no data is provided. Between soundings, the VAP linearly interpolates atmospheric state variables in time for each height level. In addition, INTERPOLATEDSONDE provides relative humidity scaled to microwave radiometer (MWR) observations.The INTERPOLATEDSONDE VAP, a continuous time-height grid of relative humidity-corrected sounding data, is intended to provide input to higher-order products, such as the Merged Soundings (MERGESONDE; Troyan 2012) VAP, which extends INTERPOLATEDSONDE by incorporating model data. The INTERPOLATEDSONDE VAP also is used to correct gaseous attenuation of radar reflectivity in products such as the KAZRCOR VAP.

54 ENVIRONMENTAL SCIENCES↗

Unleashing the Frequency: Multi-Megawatt Demonstration of 100% Renewable Power Systems with Decentralized Communication-Less Control Scheme

Power systems, which range in size from small microgrids to island systems to large regional grids, are typically managed by a central controller that requires complex communication methods and can be unreliable and pose cyber security risks in certain applications, especially when controlling a larger number of nodes. We propose an inherently robust, scalable method of integration using multiple energy storage systems and distributed energy resources, which does not require any means of dedicated communication. This method moves beyond the paradigm of controlling grid frequency at a fixed value (e.g., 60 Hz), instead allowing the frequency to fluctuate within certain limits (e.g., 59.6-60.4 Hz). With a greater operating range, the frequency can carry necessary information from energy storage systems to highly variable distributed energy resources like photovoltaics, wind, hydro, etc.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Coastal Zone Color Scanner (CZCS): Imagery of near-surface phytoplankton pigment concentrations from the first coastal ocean dynamics experiment (CODE-1), March - July 1981

As part of the first Coastal Ocean Dynamics Experiment, images of ocean color were collected from late March until late July, 1981, by the Coastal Zone Color Scanner aboard Nimbus-7. Images that had sufficient cloud-free area to be of interest were processed to yield near-surface phytoplankton pigment concentrations. These images were then remapped to a fixed equal-area grid. This report contains photographs of the digital images and a brief description of the processing methods.

Abbott, M. R.↗

Evaluation of Boundary Conditions for the Gust-Cascade Problem

Using a high-order accuracy finite-difference time-domain algorithm, the acoustic scattering from a flat-plate cascade is computed. Keeping the grid and time step fixed, the effect of four different boundary conditions on the accuracy and stability of the computed solution is compared.

Hixon, R.↗

Numerical Simulations of Crystal Growth of an Alloy Under Microgravity Conditions

The directional solidification of a dilute binary alloy (Bi-1.0 at.%Sn)is investigated. Results are obtained at a gravity level of I pg. Computations are performed in two dimensions with a fixed, non-uniform grid. The simulation involves a solution of the species concentration equation (modified to account for solute rejection at the interface) and energy equation (modified to account for phase-change) for both the solid and liquid phases, in addition to the constitutive equations for describing convective flow in the melt. The effects of conductive heat transfer in the ampoule and in a capillary tube in the sample are included. To gauge the effects of including this growth capillary tube in the apparatus, simulations both with and without the capillary tube are presented and compared. Fully transient simulations have been performed; no simplifying steady-state approximations are used, however, the influence of solute on the melting temperature at the interface is not included. Both thermal and solutal convective cells are seen to form. Convective velocities are significantly damped inside the capillary, causing less segregation due to convection. As solidification proceeds beyond the capillary tube, longitudinal segregation arises as a result of the change in cross-sectional area of solidifying material. The magnitudes of the velocities in this cell increase significantly once the solid/liquid front passes beyond the end of the capillary tube; this causes a corresponding increase in the level of radial solute segregation in the solidified material.

Simpson, James E.↗

Towards Verification of Unstructured-Grid Solvers

New methodology for verification of finite-volume computational methods using unstructured grids is presented. The discretization order properties are studied in computational windows, easily constructed within a collection of grids or a single grid. Tests are performed within each window and address a combination of problem-, solution-, and discretization/grid-related features affecting discretization error convergence. The windows can be adjusted to isolate particular elements of the computational scheme, such as the interior discretization, the boundary discretization, or singularities. Studies can use traditional grid-refinement computations within a fixed window or downscaling, a recently-introduced technique in which computations are made within windows contracting toward a focal point of interest. Grids within the windows are constrained to be consistently refined, allowing a meaningful assessment of asymptotic error convergence on unstructured grids. Demonstrations of the method are shown, including a comparative accuracy assessment of commonly-used schemes on general mixed grids and the identification of local accuracy deterioration at boundary intersections. Recommendations to enable attainment of design-order discretization errors for large-scale computational simulations are given.

Thomas, James L.↗

High-Frequency Density Oscillations from a Plasma Source Used for Simulating Low-Earth Orbit Plasma Environment

We present data from ground-based, vacuum-chamber tests demonstrating the ability to modulate the output of a plasma source capable of producing a low-Earth orbit (LEO) type plasma. We obtained plasma oscillations up to 2.5 kHz impingent on stationary test equipment, which corresponds to meter-level ionospheric structures in LEO. This plasma source is, therefore, suitable for developing scientific instruments that measure the LEO plasma environment, in situ, with meter-level spatial resolution. Measurements were made using a fixed-bias collector and an electrometer sampling at 40 kHz. A mechanical aperture was established at the output of the plasma source via two concentric grids. The outer grid was free to rotate in the azimuthal direction with respect to the fixed inner grid. An identical, alternating hole pattern in the two grids resulted in a variable aperture that cycles through 90 open/close cycles per revolution. The frequency of the plasma oscillations is limited by the mechanism used to spin the grids and the bearing assembly on which the grids rotate. Higher frequencies are obtainable by upgrading the drive mechanism, allowing the possibility of centimeter-level spatial resolution.

McTernan, Jesse K.↗

Transformation of two and three-dimensional regions by elliptic systems

A reliable linear system is presented for grid generation in 2-D and 3-D. The method is robust in the sense that convergence is guaranteed but is not as reliable as other nonlinear elliptic methods in generating nonfolding grids. The construction of nonfolding grids depends on having reasonable approximations of cell aspect ratios and an appropriate distribution of grid points on the boundary of the region. Some guidelines are included on approximating the aspect ratios, but little help is offered on setting up the boundary grid other than to say that in 2-D the boundary correspondence should be close to that generated by a conformal mapping. It is assumed that the functions which control the grid distribution depend only on the computational variables and not on the physical variables. Whether this is actually the case depends on how the grid is constructed. In a dynamic adaptive procedure where the grid is constructed in the process of solving a fluid flow problem, the grid is usually updated at fixed iteration counts using the current value of the control function. Since the control function is not being updated during the iteration of the grid equations, the grid construction is a linear procedure. However, in the case of a static adaptive procedure where a trial solution is computed and used to construct an adaptive grid, the control functions may be recomputed at every step of the grid iteration.

Mastin, C. Wayne↗

The use of solution adaptive grids in solving partial differential equations

The grid point distribution used in solving a partial differential equation using a numerical method has a substantial influence on the quality of the solution. An adaptive grid which adjusts as the solution changes provides the best results when the number of grid points available for use during the calculation is fixed. Basic concepts used in generating and applying adaptive grids are reviewed in this paper, and examples illustrating applications of these concepts are presented.

Anderson, D. A.↗

Improved Understanding of Multicentury Greenland Ice Sheet Response to Strong Warming in the Coupled CESM2‐CISM2 With Regional Grid Refinement

The simulation of ice sheet‐climate interactions, such as surface mass balance fluxes, is sensitive to model grid resolution. Here we simulate the multi‐century evolution of the Greenland Ice Sheet (GrIS) and its interaction with the climate using the Community Earth System Model version 2.2 (CESM2.2) including an interactive GrIS component (the Community Ice Sheet Model v2.1 [CISM2.1]) under an idealized warming scenario (atmospheric CO 2 increases by 1% yr -1 until quadrupling the pre‐industrial level and then is held fixed). A variable‐resolution (VR) grid with 1/4° regional refinement over the broader Arctic and 1° resolution elsewhere is applied to the atmosphere and land components, and the results are compared with conventional 1° lat‐lon grid simulations to investigate the impact of grid refinement. Compared with the 1° runs, the VR run features a slower rate of surface melt, especially over the western and northern GrIS, where the ice surface slopes gently toward the periphery. This difference pattern originates primarily from higher snow albedo and, thus, weaker albedo feedback in the VR run. The VR grid better captures the CISM ice sheet topography by reducing elevation discrepancies between CAM and CISM and is, therefore, less reliant on the downscaling algorithm, which is known to underestimate albedo gradients. The sea level rise contribution from the GrIS in the VR run is 53 mm by year 150 and 831 mm by year 350, approximately 40% and 20% less than that of the 1° runs, respectively.

Earth System Model↗

Adaptive critic design-based reinforcement learning approach in controlling virtual inertia-based grid-connected inverters

In this report, an adaptive critic design (ACD) approach is proposed to control the phase and voltage of a grid-connected virtual synchronous generator (VSG). The penetration of fast responding inertia-less power converters significantly affect the stability of the power system, especially weak systems such as micro grids. The concept of virtual inertia addresses this concern by virtually emulating the behavior of a synchronous generator. However, the conventional VSG is designed based on two conditions: (i) fixed operating point and (ii) inductive grid connections. The performance of VSGs in low-voltage semi-resistive microgrids is far from optimal. To overcome the aforementioned concerns, a heuristic dynamic programing (HDP) approach is proposed to optimally control grid-connected VSGs. The neural-network-based inherence of the HDP enables the proposed technique to adapt to any impedance angle. The HDP controller includes two subnetworks: (i) the action network that controls the system optimally and (ii) the critic network, which evaluates the effectiveness of the action network. The simulation and experimental results are provided to evaluate the effectiveness of the proposed technique. As shown, the HDP-based approach illustrates a better performance in comparison with the conventional PI-based VSG in various operating conditions.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Tri-Level Linear Programming Model for Automatic Load Shedding Using Spectral Clustering

Traditional load shedding schemes can be inadequate in grids with high renewable penetration, leading to unstable events and unnecessary grid islanding. Although for both manual and automatic operating modes load shedding areas have been predefined by grid operators, they have remained fixed, and may be sub-optimal due to dynamic operating conditions. In this work, a distributed tri-level linear programming model for automatic load shedding to avoid system islanding is presented. Preventing islanding is preferred because it reduces the need for additional load shedding besides the disconnection of transmission lines between islands. This is crucial as maintaining the local generation-demand balance is necessary to preserve frequency stability. Furthermore, uneven distribution of generation resources among islands can lead to increased load shedding, causing economic and reliability challenges. This issue is further compounded in modern power systems heavily dependent on non-dispatchable resources like wind and solar. The upper-level model uses complex power flow measurements to determine the system areas to shed load depending on actual operating conditions using a spectral clustering approach. The mid-level model estimates the area system state, while the lower-level model determines the locations and load values to be shed. The solution is practical and promising for real-world applications.

Baquedano-Aguilar, Mario D.↗

Solution of steady and unsteady transonic-vortex flows using Euler and full-potential equations

Two methods are presented for inviscid transonic flows: unsteady Euler equations in a rotating frame of reference for transonic-vortex flows and integral solution of full-potential equation with and without embedded Euler domains for transonic airfoil flows. The computational results covered: steady and unsteady conical vortex flows; 3-D steady transonic vortex flow; and transonic airfoil flows. The results are in good agreement with other computational results and experimental data. The rotating frame of reference solution is potentially efficient as compared with the space fixed reference formulation with dynamic gridding. The integral equation solution with embedded Euler domain is computationally efficient and as accurate as the Euler equations.

Kandil, Osama A.↗

Adaptive computations of multispecies mixing between scramjet nozzle flows and hypersonic freestream

The objective of this paper is to compute the expansion of a supersonic flow through an internal-external nozzle and its viscous mixing with the hypersonic flow of air. The supersonic jet may be that of a multispecies gas other than air. Calculations are performed for one case where both flows are those of air, and another case where a mixture of freon-12 and argon is discharged supersonically to mix with the hypersonic airflow. Comparisons are made between these two cases with respect to gas compositions, and fixed versus flow-adaptive grids. All the computational results are compared successfully with the wind-tunnel tests results.

Baysa, Oktay↗

A numerical solution method for acoustic radiation from axisymmetric bodies

A new and very efficient numerical method for solving equations of the Helmholtz type is specialized for problems having axisymmetric geometry. It is then demonstrated by application to the classical problem of acoustic radiation from a vibrating piston set in a stationary infinite plane. The method utilizes 'Green's Function Discretization', to obtain an accurate resolution of the waves using only 2-3 points per wave. Locally valid free space Green's functions, used in the discretization step, are obtained by quadrature. Results are computed for a range of grid spacing/piston radius ratios at a frequency parameter, omega R/c(sub 0), of 2 pi. In this case, the minimum required grid resolution appears to be fixed by the need to resolve a step boundary condition at the piston edge rather than by the length scale imposed by the wave length of the acoustic radiation. It is also demonstrated that a local near-field radiation boundary procedure allows the domain to be truncated very near the radiating source with little effect on the solution.

Caruthers, John E.↗

Combined Global Navigation Satellite Systems in the Space Service Volume

Besides providing position, velocity, and timing (PVT) for terrestrial users, the Global Positioning System (GPS) is also being used to provide PVT information for earth orbiting satellites. In 2006, F. H. Bauer, et. al., defined the Space Service Volume in the paper GPS in the Space Service Volume , presented at ION s 19th international Technical Meeting of the Satellite Division, and looked at GPS coverage for orbiting satellites. With GLONASS already operational, and the first satellites of the Galileo and Beidou/COMPASS constellations already in orbit, it is time to look at the use of the new Global Navigation Satellite Systems (GNSS) coming into service to provide PVT information for earth orbiting satellites. This presentation extends GPS in the Space Service Volume by examining the coverage capability of combinations of the new constellations with GPS GPS was first explored as a system for refining the position, velocity, and timing of other spacecraft equipped with GPS receivers in the early eighties. Because of this, a new GPS utility developed beyond the original purpose of providing position, velocity, and timing services for land, maritime, and aerial applications. GPS signals are now received and processed by spacecraft both above and below the GPS constellation, including signals that spill over the limb of the earth. Support of GPS space applications is now part of the system plan for GPS, and support of the Space Service Volume by other GNSS providers has been proposed to the UN International Committee on GNSS (ICG). GPS has been demonstrated to provide decimeter level position accuracy in real-time for satellites in low Earth orbit (centimeter level in non-real-time applications). GPS has been proven useful for satellites in geosynchronous orbit, and also for satellites in highly elliptical orbits. Depending on how many satellites are in view, one can keep time locked to the GNSS standard, and through that to Universal Time as long as at least one satellite is in view (the longest duration with no satellites in view is important in determining the maximum clock drift from GNSS time). Instantaneous position requires four satellites in view, but because orbital motion is predictable, it is possible to build up knowledge of the orbital position gradually through time without a need for constant four satellite coverage. However, it is desirable to have four satellite coverage when performing satellite maneuvers, since there can be significant changes in velocity, leading to large changes in orbit parameter, causing substantial divergence in position over time. The Space Service Volume has been defined as the volume between three thousand km altitude and geosynchronous altitude, and can be divided into medium orbit services between three thousand km altitude and eight thousand km altitude, and high orbit services above eight thousand km. The Terrestrial Service Volume includes the Earth s surface, the atmosphere, and space below the altitude of three thousand km. The Terrestrial Service Volume is the volume within which the GNSS systems will have very similar performance to the Earth surface, and satellites need only use the signals specified to provide terrestrial performance. Above three thousand km the use of signals passing by the Earth s limb becomes important, so it is desirable to have additional information on signal strength, phase delay, and group delay covering wider beam angles than are needed for terrestrial service (and which can be obtained by monitoring GNSS signals from the Earth s surface). This presentation will look at each of the new GNSS constellations in combination with GPS (GLONASS with GPS, Galileo with GPS, Beidou/COMPASS with GPS), and also at the combination of all four GNSS systems. The presentation will largely follow the format of GPS in the Space Service Volume , presenting data on the availability of one, two, three, or four of the various combinations of GNSS constellation satelles at approximately two thousand grid points evenly spaced and fixed in longitude and latitude, the duration of the longest single-fold outages (intervals when no satellites are available), and the duration of the longest four-fold outages (intervals when fewer than four satellites are available) at several altitudes. Following the original paper, we will use the altitudes of three hundred km (typical LEO satellite, and within the Terrestrial Service Volume), at three thousand km (border between Terrestrial Service Volume and Space Service Volume), at eight thousand km (the border between medium and high orbit service within the Space Service Volume), at fifteen thousand km (just below the GNSS constellations), at twenty five thousand km (just above the GNSS constellations), at thirty six thousand, five hundred km (limit of Space Service Volume definition, geosynchronous altitude), and at seventy thousand km (to show the potential usefulness of GNSS beyond geosynchronous altitude).

Force, Dale A.↗

Ionospheric calibration for single frequency altimeter measurements

This study is a preliminary analysis of the effectiveness (in terms of altimeter calibration accuracy) of various ionosphere models and the Global Positioning System (GPS) to calibrate single frequency altimeter height measurements for ionospheric path delay. In particular, the research focused on ingesting GPS Total Electron Content (TEC) data into the physical Parameterized Real-Time Ionospheric Specification Model (PRISM), which estimates the composition of the ionosphere using independent empirical and physical models and has the capability of adjusting to additional ionospheric measurements. Two types of GPS data were used to adjust the PRISM model: GPS receiver station data mapped from line-of-sight observations to the vertical at the point of interest and a grid map (generated at the Jet Propulsion Laboratory) of GPS derived TEC in a sun-fixed longitude frame. The adjusted PRISM TEC values, as well as predictions by the International Reference Ionosphere (IRI-90), a climatological (monthly mean) model of the ionosphere, were compared to TOPEX dual-frequency TEC measurements (considered as truth) for a number of TOPEX sub-satellite tracks. For a 13.6 GHz altimeter, a Total Electron Content (TEC) of 1 TECU 10(exp 16) electrons/sq m corresponds to approximately 0.218 centimeters of range delay. A maximum expected TEC (at solar maximum or during solar storms) of 10(exp 18) electrons/sq m will create 22 centimeters of range delay. Compared with the TOPEX data, the PRISM predictions were generally accurate within the TECU when the sub-satellite track of interest passed within 300 to 400 km of the GPS TEC data or when the track passed through a night-time ionosphere. If neither was the case, in particular if the track passed through a local noon ionosphere, the PRISM values differed by more than 10 TECU and by as much as 40 TECU. The IRI-90 model, with no current ability to unseat GPS data, predicted TEC to a slightly higher error of 12 TECU. The performance of PRISM is very promising for predicting TEC and will prove useful for calibrating single frequency altimeter height measurements for ionospheric path delay. When adjusted to the GPS line-of-sight data the PRISM URSI empirical model predicted TEC over a day's period to within a global error of 8.60 TECU rms during a nighttime ionosphere and 9.74 TECU rms during the day. When adjusted to the GPS derived TEC grid, the PRISM parametrized model predicted TEC to within an error of 8.47 TECU rms for a nighttime ionosphere and 12.83 TECU rms during the day. However, the grid cannot be considered globally due to the lack of sufficient numbers of GPS stations and large latitude gaps in GPS data. It is the opinion of the authors that using the PRISM model and adjusting to the global sun-fixed TEC grid regenerated with a localized weighted interpolation offers the best possibility of meeting the 10 TECU global rms (or 2 cm at 13.6 GHz) ionosphere range correction accuracy requirement of TOPEX/Poseidon and should be the subject of further study. However, it is clear that the anticipated requirement of 34 TECU global rms for TOPEX/Poseidon Follow-On (corresponding to the TOPEX/Poseidon performance) can not be met with any realizable combination of existing models and data assimilation schemes.

Schreiner, William S.↗

Grid generation for general 2-D regions using hyperbolic equations

A method for applying a hyperbolic grid generation scheme to the construction of meshes in general 2-D regions has been developed. This approach, which follows the theory developed by Steger and Chaussee (1980) and the algorithm outlined by Kinsey and Barth (1984), is based on improving local grid control. This is accomplished by adding an angle control source term to the equations and using a new algorithm for computing the volume source term. These modifications lead to superior methods for fixing the 'local' problems of hyperbolic grid generation, namely, propagation of initial discontinuities and formation of grid shocks (crossing grid lines). More importantly, a method for solving the global problem of constraining the grid with more than one boundary (internal grid generation) has been developed. These algorithms have been implemented in an interactive grid generation program and the results for several geometries are presented and discussed.

Cordova, Jeffrey Q.↗