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 91 records · Page 5

Development of an Extended Reality (XR) Tool for Earth Science Visualization

In this presentation, we will discuss our work in adapting the NASA open source XR software, the Mixed Reality Exploration Toolkit (MRET), to an earth science domain. MRET is a NASA open source XR software for rapidly building extended reality (XR) environments for NASA domain problems, e.g., pulling in CAD models of thermal vac chamber and Roman Space Telescope to do fit checks. Primarily used for hardware integration & test, we have been adapting and extending MRET for science problems. Traditionally, scientists view and analyze the result of calculated or measured observables with static 1-D, 2-D or 3-D plots. It can be difficult to identify, track and understand the evolution of key features due to poor viewing angles and the nature of flat computer screens. Additionally, numerical models, such as the NASA GEOS climate model, are almost exclusively formulated and analyzed on Eulerian grids with points fixed in space and time. However, atmospheric phenomena such as convective clouds, hurricanes and wildfire smoke plumes move with the 3-D flow field, and it is often difficult and unnatural to understand these phenomena in an Eulerian reference frame as opposed to the Lagrangian reference frame in which nature operates. As part of an Earth Science Technology Office (ESTO) proposal, we have been adapting MRET to be a scientific exploration and analysis XR tool with integrated Lagrangian Dynamics (LD) for the Goddard Earth Observing System (GEOS) numerical weather prediction model. We believe this will help scientists identify, track, and understand the evolution of Earth Science phenomena. This presentation will discuss current results in our work in developing this XR tool for Earth Science.

Thomas G Grubb↗

Development of an Extended Reality (XR) Tool for Earth Science Visualization

In this presentation, we will discuss our work in adapting the NASA open source XR software, the Mixed Reality Exploration Toolkit (MRET), to an earth science domain. MRET is a NASA open source XR software for rapidly building extended reality (XR) environments for NASA domain problems, e.g., pulling in CAD models of thermal vac chamber and Roman Space Telescope to do fit checks. Primarily used for hardware integration & test, we have been adapting and extending MRET for science problems. Traditionally, scientists view and analyze the result of calculated or measured observables with static 1-D, 2-D or 3-D plots. It can be difficult to identify, track and understand the evolution of key features due to poor viewing angles and the nature of flat computer screens. Additionally, numerical models, such as the NASA GEOS climate model, are almost exclusively formulated and analyzed on Eulerian grids with points fixed in space and time. However, atmospheric phenomena such as convective clouds, hurricanes and wildfire smoke plumes move with the 3-D flow field, and it is often difficult and unnatural to understand these phenomena in an Eulerian reference frame as opposed to the Lagrangian reference frame in which nature operates. As part of an Earth Science Technology Office (ESTO) proposal, we have been adapting MRET to be a scientific exploration and analysis XR tool with integrated Lagrangian Dynamics (LD) for the Goddard Earth Observing System (GEOS) numerical weather prediction model. We believe this will help scientists identify, track, and understand the evolution of Earth Science phenomena. This presentation will discuss current results in our work in developing this XR tool for Earth Science.

Thomas Grubb↗

DPW-VI Results Using FUN3D with Focus on k-kL-MEAH2015 (k-kL) Turbulence Model

The Common Research Model wing-body configuration is investigated with the k-kL-MEAH2015 turbulence model implemented in FUN3D. This includes results presented at the Sixth Drag Prediction Workshop and additional results generated after the workshop with a nonlinear Quadratic Constitutive Relation (QCR) variant of the same turbulence model. The workshop provided grids are used, and a uniform grid refinement study is performed at the design condition. A large variation between results with and without a reconstruction limiter is exhibited on "medium" grid sizes, indicating that the medium grid size is too coarse for drawing conclusions in comparison with experiment. This variation is reduced with grid refinement. At a fixed angle of attack near design conditions, the QCR variant yielded decreased lift and drag compared with the linear eddy-viscosity model by an amount that was approximately constant with grid refinement. The k-kL-MEAH2015 turbulence model produced wing root junction flow behavior consistent with wind tunnel observations.

Abdol-Hamid, K. S.↗

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.↗

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.↗

Multigrid solution of the Navier-Stokes equations on highly stretched grids with defect correction

Relaxation-based multigrid solvers for the steady incompressible Navier-Stokes equations are examined to determine their computational speed and robustness. Four relaxation methods with a common discretization have been used as smoothers in a single tailored multigrid procedure. The equations are discretized on a staggered grid with first order upwind used for convection in the relaxation process on all grids and defect correction to second order central on the fine grid introduced once per multigrid cycle. A fixed W(1,1) cycle with full weighting of residuals is used in the FAS multigrid process. The resulting solvers have been applied to three 2D flow problems, over a range of Reynolds numbers, on both uniform and highly stretched grids. In all cases the L(sub 2) norm of the velocity changes is reduced to 10(exp -6) in a few 10's of fine grid sweeps. The results from this study are used to draw conclusions on the strengths and weaknesses of the individual relaxation schemes as well as those of the overall multigrid procedure when used as a solver on highly stretched grids.

Sockol, Peter M.↗

Shock effects from a large impact on the moon

The paper calculates the shock and shear deformation-induced internal energy distribution associated with a major basin-forming hypervelocity meteorite impact on the moon. The Hageman and Walsh formulation of the axisymmetric two-dimensional conservation equations in finite difference form is used, and the flow field induced upon impact of an iron meteorite traveling at 15 km/sec with a gabbroic anorthosite lunar crust is calculated for sequential time steps over a grid of hoop-shaped zones fixed in space. If an energy of 5 times 10 to the 32nd power ergs for the projectile energy is required to excavate a basin of Imbrium size, than a flow field of approximately 210 km in radius approximately 19 seconds after impact is indicated. The meteorite residue, the melt, and the ejecta are discussed, and the results are compared with the Gault and Heitowit (1963) formulation.

Okeefe, J. D.↗