Search NASA⌕ Search

SEARCH · Search NASA

Results for “grid alignment”

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 163 records · Page 9

Barriers to Achieving Textbook Multigrid Efficiency (TME) in CFD

As a guide to attaining this optimal performance for general CFD problems, the table below lists every foreseen kind of computational difficulty for achieving that goal, together with the possible ways for resolving that difficulty, their current state of development, and references. Included in the table are staggered and nonstaggered, conservative and nonconservative discretizations of viscous and inviscid, incompressible and compressible flows at various Mach numbers, as well as a simple (algebraic) turbulence model and comments on chemically reacting flows. The listing of associated computational barriers involves: non-alignment of streamlines or sonic characteristics with the grids; recirculating flows; stagnation points; discretization and relaxation on and near shocks and boundaries; far-field artificial boundary conditions; small-scale singularities (meaning important features, such as the complete airplane, which are not visible on some of the coarse grids); large grid aspect ratios; boundary layer resolution; and grid adaption.

Brandt, Achi↗

Unsteady Pressure-Sensitive-Paint Shot Noise Reduction

In the study of unsteady Pressure-Sensitive Paint (uPSP) image data sets, pixel intensity values collected by Phantom cameras from multiple perspectives are used to calculate the surface pressure of an experimental model. This paper focuses on mitigation of shot noise which is introduced into the image data set by the cameras. Shot noise impacts the quality of data collected, increasing the system error in its pressure value calculations. This paper characterises statistical methods to reduce shot noise within experimental data by taking an average or weighted average across an area of pixels of their respective counts. Each combined-pixel has a higher total effective intensity value and therefore exhibits reduced noise throughout its time history. This is demonstrated by calculating the standard deviations of the combined pixel values through time, which decrease as the combined pixel area increases. Finally, these methods are compared across experimental wind tunnel conditions to examine their effectiveness at different Mach numbers. Ultimately, the results are used to recommend a blur filter type and size which will be applied to the raw image data prior to processing, as well as a minimum camera resolution for future experiments. Downstream of the filter application, the software to convert the 2D image data sets into surface pressure readings on a 3D grid use a ratio of decimation that assigns multiple pixels to each 3D grid node. This decimation ratio will also be aligned with the size of the blur filter, resulting in a more manageable data set size and compromised spatial resolution. Combined with the effect of the blur filter, the resulting pixel intensity readings will exhibit significantly reduced shot noise, reducing the magnitude of uncertainty and error in the final calculated pressure. For future tests, the camera spatial resolution must be sufficient to capture small flow and model features even after the recommended blur filter is applied. The image data set is taken from a uPSP experiment of the Space Launch System Block 1 cargo vehicle model in September 2019.

Lucy Tang↗

Unsteady Pressure-Sensitive Paint (UPSP) Shot Noise Reduction

In the study of unsteady Pressure-Sensitive Paint (uPSP) image data sets, pixel intensity values collected by Phantom cameras from multiple perspectives are used to calculate the surface pressure of an experimental model. This paper focuses on mitigation of shot noise which is introduced into the image data set by the cameras. Shot noise impacts the quality of data collected, increasing the system error in its pressure value calculations. This paper characterises statistical methods to reduce shot noise within experimental data by taking an average or weighted average across an area of pixels of their respective counts. Each combined-pixel has a higher total effective intensity value and therefore exhibits reduced noise throughout its time history. This is demonstrated by calculating the standard deviations of the combined pixel values through time, which decrease as the combined pixel area increases. Finally, these methods are compared across experimental wind tunnel conditions to examine their effectiveness at different Mach numbers. Ultimately, the results are used to recommend a blur filter type and size which will be applied to the raw image data prior to processing, as well as a minimum camera resolution for future experiments. Downstream of the filter application, the software to convert the 2D image data sets into surface pressure readings on a 3D grid use a ratio of decimation that assigns multiple pixels to each 3D grid node. This decimation ratio will also be aligned with the size of the blur filter, resulting in a more manageable data set size and compromised spatial resolution. Combined with the effect of the blur filter, the resulting pixel intensity readings will exhibit significantly reduced shot noise, reducing the magnitude of uncertainty and error in the final calculated pressure. For future tests, the camera spatial resolution must be sufficient to capture small flow and model features even after the recommended blur filter is applied. The image data set is taken from a uPSP experiment of the Space Launch System Block 1 cargo vehicle model in September 2019.

Lucy Tang↗

Electricity Markets and Long-Duration Energy Storage: A Survey of Grid Services and Revenue Streams

Purpose of Review Long Duration Energy Storage (LDES) is increasingly viewed as a potential resource for providing grid services that enhance the stability and flexibility of electricity systems. While some LDES services are integrated into existing market frameworks, traditional mechanisms may not fully account for their operational characteristics, potentially leading to undervaluation. Within this context, this paper reviews the literature and industry practices to assess potential grid services for LDES, evaluates existing compensation mechanisms, and identifies challenges to full market integration. Recent Findings We first review existing literature and identify key grid services unique to LDES, including enhancing grid resilience during extreme weather events, enabling long-term energy shifting, and providing flexible and firm energy in systems with limited dispatchable resources. Here, we also review how LDES services are compensated in current market frameworks and the challenges associated with the full realization of LDES values. Additionally, we summarize market mechanisms for storage technologies across U.S. wholesale markets. We find that some markets are adjusting incentive structures, such as incorporating storage duration in capacity accreditation, to better align with system needs and LDES contributions to the grid. However, further refinements in capacity remuneration and dispatch timeframes may be needed for more effective realization of LDES value. Summary This review evaluates potential grid services for LDES, examines existing compensation mechanisms for LDES technologies, and identifies gaps between these mechanisms and LDES operational characteristics. The review concludes by outlining potential market enhancements for more effective LDES integration and articulating additional research needs to support its efficient participation in future power systems.

Flexible resources↗

USM3D-ME Contributions to the 5th AIAA High Lift Prediction Workshop

This paper presents the results of Reynolds-averaged Navier-Stokes (RANS) simulations conducted by NASA’s flow solver, mixed-element USM3D (USM3D-ME), for the 5th AIAA High-Lift Prediction Workshop. As part of the Fixed-Grid RANS Technology Focus Group (TFG), these simulations were performed to assess the accuracy and efficiency of the USM3D-ME solutions in predicting high-lift flows. The High-Lift Common Research Model (CRM-HL) served as the primary geometry. Several CRM-HL configurations were used for three case studies: a verification study (Case 1), a configuration buildup study (Case 2), and a Reynolds-number variation study (Case 3). Overall, USM3D-ME RANS results aligned with the solutions selected by the Fixed-Grid RANS TFG and available wind tunnel data Simulations for Cases 1 and Configuration 2.1 achieved machine-zero residual convergence, with aerodynamic coefficients converging to steady-state values. However, Configurations 2.2-2.4 and Case 3 encountered iterative- and grid-convergence challenges, particularly at high angles of attack. Compared with the experimental data available for Configurations 2.2-2.4, close agreement was demonstrated at low angles of attack. However, for angles of attack approaching the maximum lift conditions, the predicted lift coefficient and pitching moment deviated from experimental values. The drag-coefficient predictions were in a relatively good agreement, however, slight overpredictions were observed at the highest angle of attack corresponding to the maximum-lift condition. Although iterative convergence for Configurations 2.2-2.4 at high angles of attack remains a persistent challenge, averaging aerodynamic coefficients over the last 5000 iterations yielded satisfactory agreement with the available wind tunnel experimental data. During the workshop, the lack of iterative convergence was attributed to the vortex structures emanating from the slat brackets. To investigate this issue further, post-workshop simulations were conducted on Configuration 2.2. In one study, RANS simulations were performed on a simplified geometry with the slat brackets removed. The second study focused on performing URANS simulations on the original Configuration 2.2 geometry. Preliminary results from both studies are presented and compared with wind tunnel data for Configuration 2.2. Consistent with the findings of other participants in the Fixed-Grid RANS TFG, this study emphasizes the necessity for further exploration and advancement in RANS technology for predicting high-lift flows.

CFD↗

USM3D-ME Contributions to the 5th AIAA High Lift Prediction Workshop

This paper presents the results of Reynolds-averaged Navier-Stokes (RANS) simulations conducted by the NASA flow solver, mixed-element USM3D (USM3D-ME), for the 5th AIAA High-Lift Prediction Workshop. As part of the Fixed-Grid RANS Technology Focus Group (TFG), these simulations were performed to assess the accuracy and efficiency of the USM3D-ME solutions in predicting high-lift flows. The High-Lift Common Research Model (CRM-HL) served as the primary geometry. Several CRM-HL configurations were used for three case studies: a verification study (Case 1), a configuration buildup study (Case 2), and a Reynolds-number variation study (Case 3). Overall, USM3D-ME RANS results aligned with the solutions selected by the Fixed-Grid RANS TFG and available wind tunnel data. Simulations for Case 1 and Configuration 2.1 achieved machine-zero residual convergence, with aerodynamic coefficients converging to steady-state values. However, Configurations 2.2-2.4 and Case 3 encountered iterative- and grid-convergence challenges, particularly at high angles of attack. Compared with the experimental data available for Configurations 2.2-2.4, close agreement was demonstrated at low angles of attack. However, for angles of attack approaching the maximum lift conditions, the predicted lift coefficient and pitching moment deviated from experimental values. The drag-coefficient predictions were in relatively good agreement, however, slight overpredictions were observed at the highest angle of attack corresponding to the maximum-lift condition. Although iterative convergence for Configurations 2.2-2.4 at high angles of attack remains a persistent challenge, averaging aerodynamic coefficients over the last 5000 iterations yielded satisfactory agreement with the available wind tunnel experimental data. During the workshop, the lack of iterative convergence was attributed to the vortex structures emanating from the slat brackets. To investigate this issue further, post-workshop simulations were conducted on Configuration 2.2. In one study, RANS simulations were performed on a simplified geometry with the slat brackets removed. The second study focused on performing unsteady RANS (URANS) simulations on the original Configuration 2.2 geometry. Preliminary results from both studies are presented and compared with wind tunnel data for Configuration 2.2. Consistent with the findings of other participants in the Fixed-Grid RANS TFG, this study emphasizes the necessity for further exploration and advancement in RANS technology for predicting high-lift flows.

Aerodynamics↗

CHESS 2025: Orthorectified airborne RGB imagery from NEON AOP surveys

This dataset provides Level 1 (L1) and Level 3 (L3) orthorectified Red-Green-Blue (RGB) imagery collected for the 2025 Colorado Headwaters Ecological Spectroscopy Study (CHESS). This high-resolution imagery is a photographic record of red, green, and blue visible light from sunlight reflected off of the Earth’s surface. The data comprise full-color images of the ground surface and are primarily intended to provide context to imaging spectroscopy and light detection and ranging (LiDAR) data. The data were acquired over three study domains in the Upper Gunnison river basin: the upper East River watershed (CRBU); Almont Triangle and Taylor Canyon (ALMO); and Upper Taylor River watershed (UPTA) between 2025-06-13 and 2025-07-15. RGB images were acquired using the PhaseOne IXM-RS150F high-resolution digital camera onboard the National Ecological Observatory Network (NEON) Airborne Observation Platform (AOP). The package data include both an L1 product comprising one camera frame per file and an L3 mosaic aligned to the Universal Transverse Mercator (UTM) Zone 13N grid and the World Geodetic System (WGS) 84 projection. Both products are provided in geotif (.tif) format at 0.1 m ground resolution. The bulk of the imagery was collected during the main CHESS field campaign from June 13 to July 15, 2025. Additional images of a portion of the Upper Taylor (UPTA) domain were collected on September 18, 2025, to fill gaps in imagery identified after the main campaign was complete. RGB camera imagery is not radiometrically calibrated, and therefore pixel values should not be exploited for scientific analysis. Pixel values have undergone a manual adjustment to enhance feature identification. The imagery is rigorously geolocated which does allow for reliable geometric information to be retrieved. To generate the orthorectified imagery, the NEON AOP camera captured visible spectrum in red, green, and blue bands. The raw images were then processed using NEON’s camera orthorectification workflow. A boresight calibration flight was made to build a complete camera, distortion, and alignment model. Color balance/white balance and exposure correction were applied to the raw RGB images. The corrected images were orthorectified by ray-tracing image pixels to a lidar-derived digital surface model (DSM) mesh using the refined camera model, outputting orthorectified raster pixels on a regular grid. Flightline-level data were mosaicked by selecting per-pixel contributions from overlapping orthorectified images using line-of-sight (LOS) zenith angle minimization to reduce edge distortions. CHESS Project Description: The Colorado Headwaters Ecological Spectroscopy Study (CHESS) comprised a multi-week airborne remote sensing and field observation campaign in the Upper Gunnison Basin, Colorado, conducted in June and July of 2025. Airborne remote sensing was conducted by the National Ecological Observatory Network Airborne Observation Platform (NEON AOP), concurrent with a field campaign run by the Rocky Mountain Biological Laboratory (RMBL), the Lawrence Berkeley National Laboratory (LBNL) and SLAC National Accelerator Laboratory Watershed Function Science Focus Area (SFA), and NASA-JPL (Jet Propulsion Laboratory) Earth Surface Mineral Dust Source Investigation (EMIT) program. Between June 10 and July 18, 2025, the NEON AOP flight team collected high-resolution aerial imaging spectroscopy and Light Detection and Ranging (LiDAR) data over three domains: the Upper East River (CRBU), Almont Triangle (ALMO), and the Upper Taylor Basin (UPTA). In coordination with the flights, a field campaign acquired ground-truth observations, including observations of vegetation composition, foliar traits, forest demography, and subsurface properties in 18 core sampling areas within the domains. Additional surface water observations were taken at over 380 point locations. All CHESS campaign datasets can be found within the CHESS ESS-DIVE data portal: https://data.ess-dive.lbl.gov/portals/chess. Funding Acknowledgement: Field and remote-sensing data acquisition was performed under a grant from the National Aeronautics and Space Administration (80NSSC24K1005). This work was also supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

2018 NEON and 2025 CHESS Campaigns↗

Numerical Simulations of a Quiet SuperSonic Technology (QueSST) Aircraft Preliminary Design

Reynolds Averaged Navier-Stokes (RANS) simulations were performed on a Lockheed Martin Quiet SuperSonic Technology (QueSST) aircraft preliminary design to assess inlet performance. The FUN3D flow solver and its adjoint-based grid refinement capability were used for the simulations in hopes of determining internal "best practices" for predicting inlet performance on top-aft-mounted inlets. Several parameters were explored including tetrahedral vs. pentahedral cells in/around the boundary-layer regions, an engine axis-aligned linear pressure sensor vs. a pressure box objective as the grid adaptation metric, and the number of grid adaptation cycles performed. Additional simulations were performed on manually refined grids for comparison with the adjoint-based adapted grids. Results showed poor agreement in predicted inlet performance on the refined grids compared to experimental data. This was true regardless of whether the refinement was adjoint-based or manual, the cell type in/near the boundary-layer regions, or the grid adaptation metric used. In addition, the 40-probe total pressure recovery was shown to decrease asymptotically as the number of adaptation cycles is increased. Solutions on the unadapted grids generally had better agreement with experimental data than their refined grid counterparts.

Inlet Performance↗

Numerical Simulations of a Quiet SuperSonic Technology (QueSST) Aircraft Preliminary Design

Reynolds Averaged Navier-Stokes (RANS) simulations were performed on the Lockheed Martin Quiet SuperSonic Technology (QueSST) aircraft preliminary design to assess inlet performance. The FUN3D flow solver and its adjoint-based grid refinement capability was used for the simulations in hopes of determining internal "best practices" for predicting inlet performance on top-aft-mounted inlets. Several parameters were explored including tetrahedral vs. pentahedral cells in/around the boundary-layer regions, an engine axis-aligned linear pressure sensor vs. a pressure box objective as the grid adaptation metric, and the number of grid adaptation cycles performed. Additional simulations were performed on manually refined grids for comparison with the adjoint-based adapted grids

Adjoint-based Grid Adaptation↗

Numerical Simulations of a Quiet SuperSonic Technology (QueSST) Aircraft Preliminary Design

Reynolds Averaged Navier-Stokes (RANS) simulations were performed on the Lockheed Martin Quiet SuperSonic Technology (QueSST) aircraft preliminary design to assess inlet performance. The FUN3D flow solver and its adjoint-based grid refinement capability was used for the simulations in hopes of determining internal "best practices" for predicting inlet performance on top-aft-mounted inlets. Several parameters were explored including tetrahedral vs. pentahedral cells in/around the boundary-layer regions, an engine axis-aligned linear pressure sensor vs. a pressure box objective as the grid adaptation metric, and the number of grid adaptation cycles performed. Additional simulations were performed on manually refined grids for comparison with the adjoint-based adapted grids

Adjoint-based Grid Adaptation↗

High Resolution Image Reconstruction from Projection of Low Resolution Images DIffering in Subpixel Shifts

In this paper, we demonstrate a simple algorithm that projects low resolution (LR) images differing in subpixel shifts on a high resolution (HR) also called super resolution (SR) grid. The algorithm is very effective in accuracy as well as time efficiency. A number of spatial interpolation techniques using nearest neighbor, inverse-distance weighted averages, Radial Basis Functions (RBF) etc. used in projection yield comparable results. For best accuracy of reconstructing SR image by a factor of two requires four LR images differing in four independent subpixel shifts. The algorithm has two steps: i) registration of low resolution images and (ii) shifting the low resolution images to align with reference image and projecting them on high resolution grid based on the shifts of each low resolution image using different interpolation techniques. Experiments are conducted by simulating low resolution images by subpixel shifts and subsampling of original high resolution image and the reconstructing the high resolution images from the simulated low resolution images. The results of accuracy of reconstruction are compared by using mean squared error measure between original high resolution image and reconstructed image. The algorithm was tested on remote sensing images and found to outperform previously proposed techniques such as Iterative Back Projection algorithm (IBP), Maximum Likelihood (ML), and Maximum a posterior (MAP) algorithms. The algorithm is robust and is not overly sensitive to the registration inaccuracies.

A 'digital' technique for manual extraction of data from aerial photography

The interpretation procedure described uses a grid cell approach. In addition, a random point is located in each cell. The procedure required that the cell/point grid be established on a base map, and identical grids be made to precisely match the scale of the photographic frames. The grid is then positioned on the photography by visual alignment to obvious features. Several alignments on one frame are sometimes required to make a precise match of all points to be interpreted. This system inherently corrects for distortions in the photography. Interpretation is then done cell by cell. In order to meet the time constraints, first order interpretation should be maintained. The data is put onto coding forms, along with other appropriate data, if desired. This 'digital' manual interpretation technique has proven to be efficient, and time and cost effective, while meeting strict requirements for data format and accuracy.

Istvan, L. B.↗

High Voltage Modeling and Testing of Transformer, Line Interface Devices, and Bulk System Components Under Electromagnetic Pulse, Geomagnetic Disturbance, and other Abnormal Transients

This project aims at quantitatively analyzing possible knowledge gaps related to the risk of Electromagnetic Pulse (EMP) impact on high-voltage, power transformers (transmission-class), that represent the most critical power grid assets. In relation to the electric infrastructure, some past studies indicate that the grid is overall fairly robust and resilient against most EMP threats. This assessment stems, in part, from estimates of the peak voltage that can be induced on a transmission power line for a nominal EMP fast-transient, compared with the typically very high basic insulation level that is featured by the high-voltage infrastructure. The project provides additional insight toward a more comprehensive risk evaluation of the grid EMP and Geomagnetic Disturbances (GMD) vulnerabilities. This work is aligned with the U.S. Dept. of Energy current research thrust in support of the development of a highly secure and resilient power infrastructure. Project milestones have been devised to address the fundamental technical issues in support of protection of large power transformer against EMP threats. More specifically, this project included parallel research tasks focused on the following topics: analysis of EMP-induced, fast-transient, current and voltages on a variety of power line configurations and EMP impact scenarios, detailed characterization of EMP damage potential due to the voltage surge, from the impact on the power line, to the transformer windings, GMD and EMP slow transients and impact from harmonics generation and reactive power flow on a complex grid layout, and identification of risk mitigation paths, with emphasis on GIC damage prevention for large power transformers, and on impact protection options from a grid stability perspective.

24 POWER TRANSMISSION AND DISTRIBUTION↗

The Effects of Neutral Inertia on Ionospheric Currents in the High-Latitude Thermosphere Following a Geomagnetic Storm

Results of an experimental and theoretical investigation into the effects of the time dependent neutral wind flywheel on high-latitude ionospheric electrodynamics are presented. The results extend our previous work which used the National Center for Atmospheric Research Thermosphere/Ionosphere General Circulation Model (NCAR TIGCM) to theoretically simulate flywheel effects in the aftermath of a geomagnetic storm. The previous results indicated that the neutral circulation, set up by ion-neutral momentum coupling in the main phase of a geomagnetic storm, is maintained for several hours after the main phase has ended and may dominate height-integrated Hall currents and field-aligned currents for up to 4-5 hours. We extend the work of Deng et al. to include comparisons between the calculated time-dependent ionospheric Hall current system in the storm-time recovery period and that measured by instruments on board the Dynamics Explorer 2 (DE 2) satellite. Also, comparisons are made between calculated field-aligned currents and those derived from DE 2 magnetometer measurements. These calculations also allow us to calculate the power transfer rate (sometimes called the Poynting flux) between the magnetosphere and ionosphere. The following conclusions have been drawn: (1) Neutral winds can contribute significantly to the horizontal ionospheric current system in the period immediately following the main phase of a geomagnetic storm, especially over the magnetic polar cap and in regions of ion drift shear. (2) Neutral winds drive Hall currents that flow in the opposite direction to those driven by ion drifts. (3) The overall morphology of the calculated field-aligned current system agrees with previously published observations for the interplanetary magnetic field (IMF) B(sub Z) southward conditions, although the region I and region 2 currents are smeared by the TI(ICM model grid resolution. (4) Neutral winds can make significant contributions to the field-aligned current system when B(sub Z) northward conditions prevail following the main phase of a storm, but can account for only a fraction of the observed currents. (5) DE 2 measurements provide a demonstration of "local" (satellite-altitude) flywheel effects. (6) On the assumption that the magnetosphere acts as an insulator, we calculate neutral-wind-induced polarization electric fields of approx. 20-30 kV in the period immediately following the geomagnetic storm.

Deng, W.↗

BGAs for High Reliability Applications

Ball Grid Array (BGA) is an important technology for utilizing higher pin counts, without the attendant handling and processing problems of the peripheral leaded packages. They are also robust in processing because of their higher pitch (0.050 inches typical), better lead rigidity, and self-alignment characteristics during reflow processing.

Ball Grid Array (BGA) lead regidity reflow process↗

Hexagonal Pixels and Indexing Scheme for Binary Images

A scheme for resampling binaryimage data from a rectangular grid to a regular hexagonal grid and an associated tree-structured pixel-indexing scheme keyed to the level of resolution have been devised. This scheme could be utilized in conjunction with appropriate image-data-processing algorithms to enable automated retrieval and/or recognition of images. For some purposes, this scheme is superior to a prior scheme that relies on rectangular pixels: one example of such a purpose is recognition of fingerprints, which can be approximated more closely by use of line segments along hexagonal axes than by line segments along rectangular axes. This scheme could also be combined with algorithms for query-image-based retrieval of images via the Internet. A binary image on a rectangular grid is generated by raster scanning or by sampling on a stationary grid of rectangular pixels. In either case, each pixel (each cell in the rectangular grid) is denoted as either bright or dark, depending on whether the light level in the pixel is above or below a prescribed threshold. The binary data on such an image are stored in a matrix form that lends itself readily to searches of line segments aligned with either or both of the perpendicular coordinate axes. The first step in resampling onto a regular hexagonal grid is to make the resolution of the hexagonal grid fine enough to capture all the binaryimage detail from the rectangular grid. In practice, this amounts to choosing a hexagonal-cell width equal to or less than a third of the rectangular- cell width. Once the data have been resampled onto the hexagonal grid, the image can readily be checked for line segments aligned with the hexagonal coordinate axes, which typically lie at angles of 30deg, 90deg, and 150deg with respect to say, the horizontal rectangular coordinate axis. Optionally, one can then rotate the rectangular image by 90deg, then again sample onto the hexagonal grid and check for line segments at angles of 0deg, 60deg, and 120deg to the original horizontal coordinate axis. The net result is that one has checked for line segments at angular intervals of 30deg. For even finer angular resolution, one could, for example, then rotate the rectangular-grid image +/-45deg before sampling to perform checking for line segments at angular intervals of 15deg.

Johnson, Gordon G.↗

Development of an implicit electromagnetic capability for a hybrid gyrokinetic ion-fluid electron model

Here we report on the development and implementation of a hybrid kinetic ion–fluid electron model for electromagnetic COGENT simulations of edge plasmas. COGENT is a finite-volume gyrokinetic code that employs a locally field-aligned coordinate system combined with a mapped multi-block grid technology to handle strongly anisotropic edge plasma turbulence. The simulation model involves the long-wavelength limit of the ion gyrokinetic equation coupled to the vorticity and Ohm's law equations for the electromagnetic field perturbations. In order to handle the fast Alfvén wave time scales, an implicit-explicit time integration approach with a physics-based preconditioner is used. The model is successfully applied to the simulations of ion-scale resistive-drift ballooning turbulence in a toroidal annulus geometry. Substantial speed-up over a fully explicit time integration approach is observed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗