Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “grid codes”

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 469 records · Page 26

An Icy Kuiper-Belt Around the Young Solar-Type Star HD 181327

HD 181327 is a young Main Sequence F5/F6 V star belonging to the Beta Pictoris moving group (age approx 12 Myr). It harbors an optically thin belt of circumstellar material at approx90 AU, presumed to result from collisions in a populat.ion of unseen planetesimals. Aims. We aim to study the dust properties in the belt in great details, and to constrain the gas-to-dust ratio. Methods. We obtained far-IR photometric observations of HD 181327 with the PACS instrument onboard the Herschel Space Observatory, complemented by new 3.2 nun observations carried with the ATCA array. The geometry of the belt is constrained with newly reduced HST /NICMOS scattered light images that break the degeneracy between the disk geometry and the dust properties. We then use the radiative transfer code GRaTer to compute a large grid of dust models, and we apply a Bayesian inference method to identify the grain models that best reproduce the SED. We attempt to detect the oxygen and ionized carbon fine-structure lines with Herschel/PACS spectroscopy, providing observables to our photochemical code ProDiMo. Results. The HST observations confirm that the dust is confined in a narrow belt. The continuum is detected with Herschel/PACS completing nicely the SED in the far-infrared. The disk is marginally resolved with both PACS and ATCA. A medium integration of the gas spectral lines only provides upper limits on the [OI] and [CII] line fluxes. We show that the HD 181327 dust disk consists of micron-sized grains of porous amorphous silicates and carbonaceous material surrounded by an import.ant layer of ice for a total dust mass of approx 0.05 stellar Mass. We discuss evidences that the grains consists of fluffy aggregates. The upper limits on the gas atomic lines do not provide unambiguous constraints: only if the PAH abundance is high, the gas mass must be lower than approx 17 Stellar Mass Conclusions. Despite the weak constraints on the gas disk, the age of HD 181327 and the properties of the dust disk suggest that it has passed the stage of gaseous planets formation. The dust reveals a population of icy planetesimals, similar to the primitive Edgeworth-Kuiper Belt, that may be a source for the future delivery of water and volatiles onto forming terrestrial planets.

Lebreton, J.↗

An Icy Kuiper Belt Around the Young Solar-type Star HD 181327

Context. HD 181327 is a young main sequence F5/F6 V star belonging to the Beta Pictoris moving group (age approx.. 12 Myr). It harbors an optically thin belt of circumstellar material at radius approx.. 90 AU, presumed to result from collisions in a population of unseen planetesimals. Aims. We aim to study the dust properties in the belt in details, and to constrain the gas-to-dust ratio. Methods. We obtained far-infrared photometric observations of HD 181327 with the PACS instrument onboard the Herschel Space Observatory, complemented by new 3.2 mm observations carried with the ATCA array. The geometry of the belt is constrained with newly reduced HST/NICMOS scattered light images that allow the degeneracy between the disk geometry and the dust properties to be broken. We then use the radiative transfer code GRaTeR to compute a large grid of models, and we identify the grain models that best reproduce the spectral energy distribution (SED) through a Bayesian analysis. We attempt to detect the oxygen and ionized carbon fine-structure lines with Herschel/PACS spectroscopy, providing observables to our photochemical code ProDiMo. Results. The HST observations confirm that the dust is confined in a narrow belt. The continuum is detected with Herschel/PACS completing nicely the SED in the far-infrared. The disk is marginally resolved with both PACS and ATCA. A medium integration of the gas spectral lines only provides upper limits on the [OI] and [CII] line fluxes.We show that the HD 181327 dust disk consists of micron-sized grains of porous amorphous silicates and carbonaceous material surrounded by an important layer of ice, for a total dust mass of approx.. 0.05 Solar Mass (in grains up to 1 mm). We discuss evidences that the grains consists of fluffy aggregates. The upper limits on the gas atomic lines do not provide unambiguous constraints: only if the PAH abundance is high, the gas mass must be lower than approx. 17 Solar Mass. Conclusions. Despite the weak constraints on the gas disk, the age of HD 181327 and the properties of the dust disk suggest that it has passed the stage of gaseous planets formation. The dust reveals a population of icy planetesimals, similar to the primitive Edgeworth-Kuiper belt, that may be a source for the future delivery of water and volatiles onto forming terrestrial planets.

HD 181327↗

Time-dependent Navier-Stokes computations of separated flows over airfoils

The application of both a central-difference and an upwind implicit approximate factorization Navier-Stokes algorithm to highly separated flow is described. Using the thin-layer approximation, both algorithms are employed to solve the low Reynolds number laminar flow around a circular cylinder with periodic shedding. Results agree with experimental Strouhal numbers to within about 10 percent. A grid refinement study is conducted using the upwind code for both the thin-layer and complete Navier-Stokes approximations. The circular cylinder solution is shown to be highly dependent upon both grid density outward from the body and grid extent. Finally, the thin-layer central-difference code is used to predict the time-accurate unsteady flow about two stalled airfoils. The NACA 0012 airfoil turbulent solutions are periodic with Strouhal numbers of about 0.1; a strong vortex is periodically shed from the leading edge. The LS(1)-0417 airfoil turbulent solution shows a high frequency oscillatory trailing edge separation. The 20 degree NACA 0012 airfoil solution is massively separated, but shows no evidence of periodicity.

Rumsey, C. L.↗

Flowfield Comparisons from Three Navier-Stokes Solvers for an Axisymmetric Separate Flow Jet

To meet new noise reduction goals, many concepts to enhance mixing in the exhaust jets of turbofan engines are being studied. Accurate steady state flowfield predictions from state-of-the-art computational fluid dynamics (CFD) solvers are needed as input to the latest noise prediction codes. The main intent of this paper was to ascertain that similar Navier-Stokes solvers run at different sites would yield comparable results for an axisymmetric two-stream nozzle case. Predictions from the WIND and the NPARC codes are compared to previously reported experimental data and results from the CRAFT Navier-Stokes solver. Similar k-epsilon turbulence models were employed in each solver, and identical computational grids were used. Agreement between experimental data and predictions from each code was generally good for mean values. All three codes underpredict the maximum value of turbulent kinetic energy. The predicted locations of the maximum turbulent kinetic energy were farther downstream than seen in the data. A grid study was conducted using the WIND code, and comments about convergence criteria and grid requirements for CFD solutions to be used as input for noise prediction computations are given. Additionally, noise predictions from the MGBK code, using the CFD results from the CRAFT code, NPARC, and WIND as input are compared to data.

Koch, L. Danielle↗

Addition of tabulated equation of state and neutrino leakage support to illinoisgrmhd

Here we have added support for realistic, microphysical, finite-temperature equations of state (EOS) and neutrino physics via a leakage scheme to illinoisgrmhd, an open-source GRMHD code for dynamical spacetimes in the einstein toolkit. These new features are provided by two new, nrpy+-based codes: nrpyeos, which performs highly efficient EOS table lookups and interpolations, and nrpyleakage, which implements a new, adaptive mesh refinement (AMR)-capable neutrino leakage scheme in the einstein toolkit. We have performed a series of strenuous validation tests that demonstrate the robustness of these new codes, particularly on the Cartesian AMR grids provided by carpet. Furthermore, we show results from fully dynamical GRMHD simulations of single unmagnetized neutron stars, and magnetized binary neutron star mergers. This new version of illinoisgrmhd, as well as nrpyeos and nrpyleakage, is pedagogically documented in jupyter notebooks and fully open source. The codes will be proposed for inclusion in an upcoming version of the einstein toolkit.

79 ASTRONOMY AND ASTROPHYSICS↗

Accurate modeling of parallel scientific computations

Scientific codes are usually parallelized by partitioning a grid among processors. To achieve top performance it is necessary to partition the grid so as to balance workload and minimize communication/synchronization costs. This problem is particularly acute when the grid is irregular, changes over the course of the computation, and is not known until load time. Critical mapping and remapping decisions rest on the ability to accurately predict performance, given a description of a grid and its partition. This paper discusses one approach to this problem, and illustrates its use on a one-dimensional fluids code. The models constructed are shown to be accurate, and are used to find optimal remapping schedules.

Nicol, David M.↗

The simulation of hypersonic flow using a zonal approach

The development of a Navier-Stokes code for the simulation of hypersonic flow will be described. The code uses a zonal philosophy which helps alleviate the grid generation problem for complex configurations and allows an efficient solution procedure. The code has been used to simulate equilibrium flow past a simple generic hypersonic configuration. Validation of the code on a blunt cone has also been conducted.

Flores, Jolen↗

The Impact of Back-Sputtered Carbon on the Accelerator Grid Wear Rates of the NEXT and NSTAR Ion Thrusters

A study was conducted to quantify the impact of back-sputtered carbon on the downstream accelerator grid erosion rates of the NEXT (NASA's Evolutionary Xenon Thruster) Long Duration Test (LDT1). A similar analysis that was conducted for the NSTAR (NASA's Solar Electric Propulsion Technology Applications Readiness Program) Life Demonstration Test (LDT2) was used as a foundation for the analysis developed herein. A new carbon surface coverage model was developed that accounted for multiple carbon adlayers before complete surface coverage is achieved. The resulting model requires knowledge of more model inputs, so they were conservatively estimated using the results of past thin film sputtering studies and particle reflection predictions. In addition, accelerator current densities across the grid were rigorously determined using an ion optics code to determine accelerator current distributions and an algorithm to determine beam current densities along a grid using downstream measurements. The improved analysis was applied to the NSTAR test results for evaluation. The improved analysis demonstrated that the impact of back-sputtered carbon on pit and groove wear rate for the NSTAR LDT2 was negligible throughout most of eroded grid radius. The improved analysis also predicted the accelerator current density for transition from net erosion to net deposition considerably more accurately than the original analysis. The improved analysis was used to estimate the impact of back-sputtered carbon on the accelerator grid pit and groove wear rate of the NEXT Long Duration Test (LDT1). Unlike the NSTAR analysis, the NEXT analysis was more challenging because the thruster was operated for extended durations at various operating conditions and was unavailable for measurements because the test is ongoing. As a result, the NEXT LDT1 estimates presented herein are considered preliminary until the results of future posttest analyses are incorporated. The worst-case impact of carbon back-sputtering was determined to be the full power operating condition, but the maximum impact of back-sputtered carbon was only a four percent reduction in wear rate. As a result, back-sputtered carbon is estimated to have an insignificant impact on the first failure mode of the NEXT LDT at all operating conditions.

Soulas, George C.↗

The Impact of Back-Sputtered Carbon on the Accelerator Grid Wear Rates of the NEXT and NSTAR Ion Thrusters

A study was conducted to quantify the impact of back-sputtered carbon on the downstream accelerator grid erosion rates of the NASA's Evolutionary Xenon Thruster (NEXT) Long Duration Test (LDT1). A similar analysis that was conducted for the NASA's Solar Electric Propulsion Technology Applications Readiness Program (NSTAR) Life Demonstration Test (LDT2) was used as a foundation for the analysis developed herein. A new carbon surface coverage model was developed that accounted for multiple carbon adlayers before complete surface coverage is achieved. The resulting model requires knowledge of more model inputs, so they were conservatively estimated using the results of past thin film sputtering studies and particle reflection predictions. In addition, accelerator current densities across the grid were rigorously determined using an ion optics code to determine accelerator current distributions and an algorithm to determine beam current densities along a grid using downstream measurements. The improved analysis was applied to the NSTAR test results for evaluation. The improved analysis demonstrated that the impact of back-sputtered carbon on pit and groove wear rate for the NSTAR LDT2 was negligible throughout most of eroded grid radius. The improved analysis also predicted the accelerator current density for transition from net erosion to net deposition considerably more accurately than the original analysis. The improved analysis was used to estimate the impact of back-sputtered carbon on the accelerator grid pit and groove wear rate of the NEXT Long Duration Test (LDT1). Unlike the NSTAR analysis, the NEXT analysis was more challenging because the thruster was operated for extended durations at various operating conditions and was unavailable for measurements because the test is ongoing. As a result, the NEXT LDT1 estimates presented herein are considered preliminary until the results of future post-test analyses are incorporated. The worst-case impact of carbon back-sputtering was determined to be the full power operating condition, but the maximum impact of back-sputtered carbon was only a 4 percent reduction in wear rate. As a result, back-sputtered carbon is estimated to have an insignificant impact on the first failure mode of the NEXT LDT1 at all operating conditions.

Soulas, George C.↗

Representing the Future Role of Hydropower and Pumped Storage Hydropower (PSH) in Electricity Planning Tools

Existing tools for long-term electric sector planning struggle to represent hydropower's nuanced site-specific technical and operating characteristics, which depend on technical specifications as well as water management practices and regulations. As a result, long-term planning models and tools insufficiently characterize hydropower value and incentives, and they cannot fully represent the role hydropower can play in a future electricity system that could include a high penetration of variable wind and solar generation, battery storage, and other low-carbon technologies. This presentation demonstrates the culmination of a multi-year effort to enhance hydropower representations in electricity planning models at the National Renewable Energy Laboratory (NREL), as part of the U.S. Department of Energy (USDOE) HydroWIRES Initiative. New modeling techniques are demonstrated using the NREL Regional Energy Deployment System (ReEDS), an open-access electric sector capacity expansion model used extensively in a wide range of technology deployment and integration analysis, including the 2016 USDOE Hydropower Vision. ReEDS uses a least-cost optimization approach to understand investment and operation of electricity generation, storage, and transmission technologies under future scenarios of electricity technology innovation, demand, policy, and other sectoral drivers. ReEDS was modified to better represent value and opportunities for both pumped storage hydropower (PSH) and hydropower systems without pumping. We incorporated a new national closed-loop PSH resource and cost assessment to explore new PSH deployment opportunities and added plant-level data to better represent the existing PSH fleet. New upgrade pathways enable opportunities for enhanced hydropower flexibility by adding pumps, upgrading dispatchability, increasing capacity, or increasing energy availability. The model was also modified to better represent the value of long-duration energy storage beyond diurnal time scales, allowing both hydropower and PSH to better balance energy supply and demand variations in high-renewable systems. These new features are demonstrated under reference and high-renewable futures and a range of sensitivity scenarios to understand which hydropower and PSH deployment and upgrade opportunities are the most attractive. These scenarios indicate potential for new closed-loop PSH deployment and for hydropower flexibility improvements to have important impacts on long-term electricity system emissions and economic outcomes. Increasing flexibility of the existing hydropower fleet can reduce the need to invest in new flexible grid technologies and help achieve decarbonization goals. Systems with sufficient energy storage could also be valuable for balancing seasonal differences in renewable energy availability, particularly from solar energy. The methods developed for ReEDS and subsequent scenario results reveal important considerations for future hydropower and grid system planning, and all data and code is freely available in a public code repository for use throughout the hydropower industry.

capacity expansion↗

Efficient Cache use for Stencil Operations on Structured Discretization Grids

We derive tight bounds on the cache misses for evaluation of explicit stencil operators on structured grids. Our lower bound is based on the isoperimetrical property of the discrete octahedron. Our upper bound is based on a good surface to volume ratio of a parallelepiped spanned by a reduced basis of the interference lattice of a grid. Measurements show that our algorithm typically reduces the number of cache misses by a factor of three, relative to a compiler optimized code. We show that stencil calculations on grids whose interference lattice have a short vector feature abnormally high numbers of cache misses. We call such grids unfavorable and suggest to avoid these in computations by appropriate padding. By direct measurements on a MIPS R10000 processor we show a good correlation between abnormally high numbers of cache misses and unfavorable three-dimensional grids.

Frumkin, Michael↗

Numerical Simulation of Ion Thruster Optics

A three-dimensional simulation code (ffx) designed to analyze ion thruster optics is described. It is an extension of an earlier code and includes special features like the ability to model a wide range of grid geometries, cusp details, and mis-aligned aperture pairs to name a few. However, the principle reason for advancing the code was in the study of ion optics erosion. Ground based testing of ion thruster optics, essential to the understanding of the processes of grid erosion, can be time consuming and costly. Simulation codes that can accurately predict grid lifetimes and the physical mechanisms of grid erosion can be of great utility in the development of future ion thruster optics designed for more ambitious applications. Results of simulations are presented that describe wear profiles for several standard and nonstandard aperture geometries, such as those grid sets with square- or slotted-hole layout patterns. The goal of this paper will be to introduce the methods employed in the ffx code and to briefly demonstrate their use.

Rawlin, Vincent K.↗

Multitasking the INS3D-LU code on the Cray Y-MP

This paper presents the results of multitasking the INS3D-LU code on eight processors. The code is a full Navier-Stokes solver for incompressible fluid in three dimensional generalized coordinates using a lower-upper symmetric-Gauss-Seidel implicit scheme. This code has been fully vectorized on oblique planes of sweep and parallelized using autotasking with some directives and minor modifications. The timing results for five grid sizes are presented and analyzed. The code has achieved a processing rate of over one Gflops.

Fatoohi, Rod↗

Body conforming grids for general unsteady airfoil motion

A method for computing a new body conforming grid at each time step of an unsteady flow calculation is presented. The method is based on a spring system analogy. It simply updates the grid from the previous time step, hence it requires that an initial grid be created using another method. The method is designed to be independent of the flow equations used, so it can be easily transferred between aerodynamic codes with the same type of initial grid (i.e., O-type). The initial development was done using an O-type grid. It was then extended to C-type grids which required only minor changes to the routine. Results are presented for both O-type and C-type grids. The quality of the updated grids were examined by observing their movement and the flow results they produced. The effect due to the shape of the prescribed grid lines for the O-type grids was examined for a NACA 0012 airfoil undergoing harmonic pitch. The time history of the pressure coefficients of this case was compared with the interpolation scheme and experimental data.

Levine, M.↗

Black box multigrid

The black box multigrid (BOXMG) code, which only needs specification of the matrix problem for application in the multigrid method was investigated. It is contended that a major problem with the multigrid method is that each new grid configuration requires a major programming effort to develop a code that specifically handles that grid configuration. The SOR and ICCG methods only specify the matrix problem, no matter what the grid configuration. It is concluded that the BOXMG does everything else necessary to set up the auxiliary coarser problems to achieve a multigrid solution.

Dendy, J. E., Jr.↗

PARCEQ2D heat transfer grid sensitivity analysis

The material presented in this paper is an extension of two-dimensional Aeroassist Flight Experiment (AFE) results shown previously. This study has focused on the heating rate calculations to the AFE obtained from an equilibrium real gas code, with attention placed on the sensitivity of grid dependence and wall temperature. Heat transfer results calculated by the PARCEQ2D code compare well with those computed by other researchers. Temperature convergence in the case of kinetic transport has been accomplished by increasing the wall temperature gradually from 300 K to the wall temperature of 1700 K.

Saladino, Anthony J.↗

User Manual for Beta Version of TURBO-GRD: A Software System for Interactive Two-Dimensional Boundary/ Field Grid Generation, Modification, and Refinement

TURBO-GRD is a software system for interactive two-dimensional boundary/field grid generation. modification, and refinement. Its features allow users to explicitly control grid quality locally and globally. The grid control can be achieved interactively by using control points that the user picks and moves on the workstation monitor or by direct stretching and refining. The techniques used in the code are the control point form of algebraic grid generation, a damped cubic spline for edge meshing and parametric mapping between physical and computational domains. It also performs elliptic grid smoothing and free-form boundary control for boundary geometry manipulation. Internal block boundaries are constructed and shaped by using Bezier curve. Because TURBO-GRD is a highly interactive code, users can read in an initial solution, display its solution contour in the background of the grid and control net, and exercise grid modification using the solution contour as a guide. This process can be called an interactive solution-adaptive grid generation.

Choo, Yung K.↗

Creating a Simple Single Computational Approach to Modeling Rarefied and Continuum Flow About Aerospace Vehicles

We proposed to create a single computational code incorporating methods that can model both rarefied and continuum flow to enable the efficient simulation of flow about space craft and high altitude hypersonic aerospace vehicles. The code was to use a single grid structure that permits a smooth transition between the continuum and rarefied portions of the flow. Developing an appropriate computational boundary between the two regions represented a major challenge. The primary approach chosen involves coupling a four-speed Lattice Boltzmann model for the continuum flow with the DSMC method in the rarefied regime. We also explored the possibility of using a standard finite difference Navier Stokes solver for the continuum flow. With the resulting code we will ultimately investigate three-dimensional plume impingement effects, a subject of critical importance to NASA and related to the work of Drs. Forrest Lumpkin, Steve Fitzgerald and Jay Le Beau at Johnson Space Center. Below is a brief background on the project and a summary of the results as of the end of the grant.

Goldstein, David B.↗