Search NASA⌕ Search

SEARCH · Search NASA

Results for “Irregular 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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS2008

Mapping Snow Grain Size over Greenland from MODIS

This paper presents a new automatic algorithm to derive optical snow grain size (SGS) at 1 km resolution using Moderate Resolution Imaging Spectroradiometer (MODIS) measurements. Differently from previous approaches, snow grains are not assumed to be spherical but a fractal approach is used to account for their irregular shape. The retrieval is conceptually based on an analytical asymptotic radiative transfer model which predicts spectral bidirectional snow reflectance as a function of the grain size and ice absorption. The analytical form of solution leads to an explicit and fast retrieval algorithm. The time series analysis of derived SGS shows a good sensitivity to snow metamorphism, including melting and snow precipitation events. Preprocessing is performed by a Multi-Angle Implementation of Atmospheric Correction (MAIAC) algorithm, which includes gridding MODIS data to 1 km resolution, water vapor retrieval, cloud masking and an atmospheric correction. MAIAC cloud mask (CM) is a new algorithm based on a time series of gridded MODIS measurements and an image-based rather than pixel-based processing. Extensive processing of MODIS TERRA data over Greenland shows a robust performance of CM algorithm in discrimination of clouds over bright snow and ice. As part of the validation analysis, SGS derived from MODIS over selected sites in 2004 was compared to the microwave brightness temperature measurements of SSM\I radiometer, which is sensitive to the amount of liquid water in the snowpack. The comparison showed a good qualitative agreement, with both datasets detecting two main periods of snowmelt. Additionally, MODIS SGS was compared with predictions of the snow model CROCUS driven by measurements of the automatic whether stations of the Greenland Climate Network. We found that CROCUS grain size is on average a factor of two larger than MODIS-derived SGS. Overall, the agreement between CROCUS and MODIS results was satisfactory, in particular before and during the first melting period in mid-June. Following detailed time series analysis of SGS for four permanent sites, the paper presents SGS maps over the Greenland ice sheet for the March-September period of 2004.

Lyapustin, Alexei↗
NASA NTRSDate not supplied

Determining True Sensor Spatial Resolution of Very High Resolution Optical Imagery

Some satellite data is delivered in images with gridded pixels. This gridded pixel size is often assumed to be the spatial resolution of the satellite sensor; however, this is not always the case. An image can be grided to any arbitrary pixel size, but the sensor resolution will remain constant. For example, an image with a pixel grid size much smaller than the sensor resolution will appear blurry along what should be sharp transitions. This discrepancy between an image’s pixel size and true sensor spatial resolution can be the source of much confusion and even misinformation among data users, which may lead them to waste time and resources on using images that do not suit their spatial resolution needs. This presentation will highlight our evaluation of the true spatial resolution of various government and commercial images in the pixel size range of 0.3 m to 60 m. Images evaluated include ESA’s Sentinel-2 (60 m, 20 m, 10 m pixels), USGS’s Landsat 8/9 (30 m & 15 m pixels), Planet’s SuperDoves (3 m pixels), BlackSky’s Globals (~1 m pixels), and the optical bands of Maxar’s WorldView-2 (2.4 m – 0.41 m pixels) and WorldView-3 (1.38 m – 0.31 m pixels). Our evaluation of true sensor spatial resolution, or ‘footprint size’ is based on the sensor’s line spread function (LSF). We calculate the width at half the height of the LSF to find the full width at half maximum (FWHM). The FWHM is how we report sensor spatial resolution. Different objects are examined for constructing the LSF depending on the sensor spatial resolution. Coarser resolution sensors in this evaluation such as Sentinel-2 and Landsat 8/9 are examined at bridges over a dark water background. The bright bridge acts as a line impulse, giving a sensor’s line spread function (LSF) in one direction. Additionally, we simulate the impacts of bridge width on the apparent LSF to obtain a true LSF without the effects of bridge width for these sensors. Finer resolution sensors will image the irregularities in bridges such as trusses, sidewalks, and in some cases painted lines, interfering with the LSF construction. Instead, these sensors are evaluated at large (60 m – 140 m) black and white checkerboards known as Cal/Val sites. At these locations, the image’s transition from black to white is extracted as an edge spread function (ESF). We calculate the derivative of this ESF to obtain the sensor’s LSF. From there, we find the FWHM as we do for the coarser resolution images. With the FWHM and pixel size, we determine how over- or under-sampled the images are. When the ratio of a sensor’s spatial resolution and the gridded image’s pixel size is less than 1, the image is considered under-sampled. In this case, each pixel’s information is unique but only a portion of that pixel’s ground area has been measured. On the other side, if the ratio is greater than 1, the image is considered over-sampled. That is, each pixel’s information is sourced from within the ground extent of the pixel and some extent outside additionally. We will show the true spatial resolution and the extent of over-/under-sampling in the imagery from ESA’s Sentinel-2 (60 m, 20 m, 10 m pixels), USGS’s Landsat 8/9 (30 m & 15 m pixels), Planet’s SuperDoves (3 m pixels), BlackSky’s Globals (~1 m pixels), and the optical bands of Maxar’s WorldView-2 (2.4 m – 0.41 m pixels) and WorldView-3 (1.38 m – 0.31 m pixels).

Alana Semple↗
DOE OSTI2024

Refactoring the elastic–viscous–plastic solver from the sea ice model CICE v6.5.1 for improved performance

This study focuses on the performance of the elastic–viscous–plastic (EVP) dynamical solver within the sea ice model, CICE v6.5.1. The study has been conducted in two steps. First, the standard EVP solver was extracted from CICE for experiments with refactored versions, which are used for performance testing. Second, one refactored version was integrated and tested in the full CICE model to demonstrate that the new algorithms do not significantly impact the physical results. The study reveals two dominant bottlenecks, namely (1) the number of Message Parsing Interface (MPI) and Open Multi-Processing (OpenMP) synchronization points required for halo exchanges during each time step combined with the irregular domain of active sea ice points and (2) the lack of single-instruction, multiple-data (SIMD) code generation. The standard EVP solver has been refactored based on two generic patterns. The first pattern exposes how general finite differences on masked multi-dimensional arrays can be expressed in order to produce significantly better code generation by changing the memory access pattern from random access to direct access. The second pattern takes an alternative approach to handle static grid properties. The measured single-core performance improvement is more than a factor of 5 compared to the standard implementation. The refactored implementation of strong scales on the Intel® Xeon® Scalable Processors series node until the available bandwidth of the node is used. For the Intel® Xeon® CPU Max series, there is sufficient bandwidth to allow the strong scaling to continue for all the cores on the node, resulting in a single-node improvement factor of 35 over the standard implementation. This study also demonstrates improved performance on GPU processors.

58 GEOSCIENCES↗
DOE OSTI2023

A Real-Time ANPC Inverter Digital Twin with Integrated Design-For-Trust

The demand for renewable energy has increased over the last few years, and so has the demand for greater expectations within the energy market. This increasing trend has been accompanied by more significant usage of internet-connected devices (IoT), leading to critical electrical infrastructure being connected to the internet. Implementing internet connectivity with such devices and systems provides benefits such as improving the system's performance, facilitating irregularity and anomaly mitigation, and providing additional situational awareness for enhanced decision-making. However, enhancing the connected system with IoT introduces a drawback – a greater vulnerability to cyber-attacks. Cyber-attacks targeting critical infrastructure in the electrical sector have occurred in the United States and Ukraine. These cyber-attacks highlight and expose vulnerabilities that a system inherits when connecting to the internet. These attacks left thousands of customers without electricity for hours until operators could regain control of the electric utility grid. Therefore, to address the vulnerabilities of an internet-connected power electronic device, this work focused on the hardware layer of the system. Implementing a cyber-control system inside the hardware layer can significantly reduce the possibility of an attacker patching malicious controller firmware into a photovoltaic grid-connected inverter, thus mitigating the likelihood that the inverter becomes inactive a cyber-attack scenario. With this mitigation technique, if a cyberattack is successful and an attacker gains control of the network, a cyber-defense technique is in place to mitigate the impact of the cyber-attack. This additional protection layer was developed based on an innovative concept known as Digital Twin (DT). A DT, in this case, replicates an Active-Neutral Point Clamped (ANPC) inverter and was designed using a hardware language known as VHDL (Very High-Speed Integrated Circuit Hardware Description Language) and applied to Field-Programmable-GateArray (FPGA). The DT is embedded within the FPGA and contained in a controller board, the UCB (Unified Controller Board), developed by the University of Arkansas electrical engineering team. This UCB also contains two Digital Signal Processors (DSPs) responsible for generating associated signals to control an authentic physical inverter. These DSP signals are received and processed by the FPGA that implements the DT of an ANPC; in other words, it simulates in realtime the expected output of an actual ANPC inverter using the signals from the DSP. When a new firmware is ready to be patched, the DT provides output signals simulating behavior that a real ANPC inverter would generate with the new firmware. The new firmware is tested to check if it meets all the operational requirements established using a Design-For-Trust technique (DFTr). If the new firmware fails in at least one of the DFT tests, it is considered malicious and must be rejected. This work is divided into sections, such as Background, which explains the pieces that were used and the strategy behind this work; Process and Procedure, which explains the methodology that was adopted to prove the reliability and effectiveness of this work; Results and Discussion, where the simulations and results are described and explained; followed by Conclusion and Future work section, which concludes this work and adds possible future projects to continue this work further.

do Amaral Custodio, Paulo Vitor↗
NASA NTRS1989

Methods for local gravity field approximation

The most widely known modern method for estimating gravity field values from observed data is least-squares collocation. Its advantages are that it can make estimates at arbitrary locations based on irregularly spaced observations, and that it makes use of statistical information about errors in the input data while providing corresponding information about the quality of the output estimates. Disadvantages of collocation include the necessity of inverting square matrices of dimension equal to the number of data values and the need to assume covariance models for the gravity field and the data errors. Fourier methods are an important alternative to collocation; having the advantage of greater computational efficiency, but requiring data estimates to be on a regular grid and not using or providing statistical accuracy information. The GEOFAST algorithm is an implementation of collocation that achieves high computational efficiency by transforming the estimation equations into the frequency domain where an accurate approximation may be made to reduce the workload. The forward and inverse Fast Fourier Transforms (FFTs) are utilized. The accuracy and computational efficiency of the GEOFAST algorithm is demonstrated using two sets of synthetic gravity data: marine gravity for an ocean trench region including wavelengths longer than 200 km; and local land gravity containing wavelengths as short as 5 km. These results are discussed along with issues such as the advantages of first removing reference field models before carrying out the estimation algorithm.

Sailor, R. V.↗
NASA NTRS1993

Three-dimensional numerical simulation of near-surface flows over the Martian north polar cap

Measurements made by Viking Lander VL-2 (48 N) have shown that the near-surface wind and temperature regime on Mars displays striking similarities to terrestrial counterparts. The diurnal radiative cycle is responsible for establishment of a well-defined thermal circulation in which downslope (Katabatic) flows prevail during the nighttime hours and weak upslope (anabatic) conditions prevail during the daytime. Previous work has indicated that the slope flows are much like those found on Earth, particularly the Katabatic winds, which show striking similarities to drainage flows observed over Antarctica. The low-level wind regime appears to be an important factor in the scouring of the martian landscape. The north polar cap shows evidence of eolian features such as dunes, frost streaks, and grooves from Viking imagery. The direction of the prevailing wind can in cases be inferred from the eolian features. We examine the thermally induced flows that result from the radiative heating and cooling of the martian north polar region using a comprehensive three-dimensional atmospheric mesoscale numerical model. The same model has been used previously for simulation of Antarctic Katabatic winds. The model equations are written in terrain-following coordinates to allow for irregular terrain; prognostic equations include the flux forms of the horizontal momentum equations, temperature, continuity. A surface energy budget equation is also incorporated in which the surface temperature is determined. Explicit parameterization of both terrestrial (longwave) and solar (shortwave) radiation is included. Turbulent transfer of heat and momentum in the martian atmosphere is assumed to follow the similarity expressions in the surface boundary layer on Earth. The terrain heights for the martian north polar region have been obtained from the U.S. Geological Survey map and digitized onto a 57x57 grid with a spacing of 75 km. The resulting terrain map is shown in Fig. 1. The vertical grid consists of 15 levels (delta equals 0.998, 0.99, 0.98, 0.97, 0.96, 0.94, 0.92, 0.90, 0.85, 0.775, 0.70, 0.60, 0.50, 0.30, 0.10). The high resolution in the lower atmosphere is necessary to capture details of the boundary layer flows. The lowest level corresponds to a height of approximately 20 m above the ground, the seconds level 100 m.

Parish, Thomas R.↗
DOE OSTI2019

Long-term Spot Stability in the Post-common-envelope Binary QS Vir

We observed the post-common-envelope eclipsing binary with a white dwarf component, QS Vir, using the 1.88 m telescope of Kotammia Observatory in Egypt. The new observations were analyzed together with all multicolor light curves available online (sampling a period of 25 yr), using a full-feature binary system modeling software based on Roche geometry. This is the first time complete photometric modeling was done with most of these data. QS Vir is a detached system, with the red dwarf component underfilling its Roche lobe by a small margin. All light curves feature out-of-eclipse variability that is associated with ellipsoidal variation, mutual irradiation and irregularities in surface brightness of the tidally distorted and magnetically active red dwarf. We tested models with one, two, and three dark spots and found that one spot is sufficient to account for the light curve asymmetry in all data sets, although this does not rule out the presence of multiple spots. We also found that a single spotted model cannot fit light curves observed simultaneously in different filters. Instead, each filter requires a different spot configuration. To thoroughly explore the parameter space of spot locations, we devised a grid-search procedure and used it to find consistent solutions. Based on this, we conclude that the dark spot responsible for light curve distortions has been stable for the past 15 yr, after a major migration that happened between 1993 and 2002, possibly due to a flip-flop event.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1984

Heat transfer distributions around nominal ice accretion shapes formed on a cylinder in the NASA Lewis Icing Research Tunnel

Local heat transfer coefficients were obtained on irregular cylindrical shapes which typify the accretion of ice on circular cylinders in cross flow. The shapes were 2, 5, and 15 min accumulations of glaze ice and 15 min accumulation of rime ice. These icing shapes were averaged axially to obtain a nominal shape of constant cross section for the heat transfer tests. Heat transfer coefficients were also measured around the cylinder with no ice accretion. The models were run in a 15.2 x 68.6 cm (6 x 27 in.) wind tunnel at several velocities. The models were also run with a turbulence producing grid which gave about 3.5 percent turbulence. The effect of roughness was also simulated with sand grains glued to the surface. Results are presented as Nusselt number versus angle from the stagnation line for the smooth and rough models for both high and low levels of free stream turbulence. Roughness of the surface in the region prior to flow separation plays a major role in determining the heat transfer distribution. Free stream turbulence does not affect the distribution of heat transfer in this region but raises the level by a nearly uniform amount. For the rime shape, roughness had a larger effect in the near wedge shaped region past the initial separation point.

Van Fossen, G. J.↗
DOE OSTI2022

Smoothed boundary method for simulating incompressible flow in complex geometries

Simulating flow through porous media with explicit considerations of complex microstructures is very challenging using conventional sharp-interface methods because of the difficulties in generating meshes conformal to complex geometries. In this work, a diffuse interface embedded boundary method known as the Smoothed Boundary Method (SBM) is utilized to facilitate simulations of fluid dynamics involving complex geometries. In diffuse-interface methods, the geometry is described by a domain parameter. The SBM allows the straightforward reformulation of the time-dependent Navier-Stokes equations in terms of this domain parameter, using only algebraic identities. Thus, enforcing the appropriate boundary conditions at the irregular embedded boundary is greatly simplified. Adaptive mesh refinement is used to increase the accuracy of the diffuse interface method by allowing thinner interfacial thickness to be used in the domain parameter. Here, the SBM-formulated Navier-Stokes equations are solved with the Finite Difference Method on refined mesh systems. Sharp-interface Finite Element Method simulations using the commercial software COMSOL on body-conforming meshes are also provided for comparison. Favorable agreement between the two methods is observed. Since it is no longer necessary for the mesh to conform to the complex geometry, the grid system for the SBM simulations can be generated rapidly and without additional manual interventions, making the entire simulation process more expedient.

36 MATERIALS SCIENCE↗
DOE OSTI2018

De novo design of a fluorescence-activating β-barrel

The regular arrangements of β-strands around a central axis in β-barrels and of α-helices in coiled coils contrast with the irregular tertiary structures of most globular proteins, and have fascinated structural biologists since they were first discovered. Simple parametric models have been used to design a wide range of α-helical coiled-coil structures, but to date there has been no success with β-barrels. Here we show that accurate de novo design of β-barrels requires considerable symmetry-breaking to achieve continuous hydrogen-bond connectivity and eliminate backbone strain. We then build ensembles of β-barrel backbone models with cavity shapes that match the fluorogenic compound DFHBI, and use a hierarchical grid-based search method to simultaneously optimize the rigid-body placement of DFHBI in these cavities and the identities of the surrounding amino acids to achieve high shape and chemical complementarity. The designs have high structural accuracy and bind and fluorescently activate DFHBI in vitro and in Escherichia coli, yeast and mammalian cells. Finally, this de novo design of small-molecule binding activity, using backbones custom-built to bind the ligand, should enable the design of increasingly sophisticated ligand-binding proteins, sensors and catalysts that are not limited by the backbone geometries available in known protein structures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2001

Predictability of the California Current System

The physical and biological oceanography of the Southern California Bight (SCB), a highly productive subregion of the California Current System (CCS) that extends from Point Conception, California, south to Ensenada, Mexico, continues to be extensively studied. For example, the California Cooperative Oceanic Fisheries Investigations (CalCOFI) program has sampled this region for over 50 years, providing an unparalleled time series of physical and biological data. However, our understanding of what physical processes control the large-scale and mesoscale variations in these properties is incomplete. In particular, the non-synoptic and relatively coarse spatial sampling (70km) of the hydrographic grid does not completely resolve the mesoscale eddy field (Figure 1a). Moreover, these unresolved physical variations exert a dominant influence on the evolution of the ecosystem. In recent years, additional datasets that partially sample the SCB have become available. Acoustic Doppler Current Profiler (ADCP) measurements, which now sample upper-ocean velocity between stations, and sea level observations along TOPEX tracks give a more complete picture of the mesoscale variability. However, both TOPEX and ADCP are well-sampled only along the cruise or orbit tracks and coarsely sampled in time and between tracks. Surface Lagrangian drifters also sample the region, although irregularly in time and space. SeaWiFS provides estimates of upper-ocean chlorophyll-a (chl-alpha), usually giving nearly complete coverage for week-long intervals, depending on cloud coverage. Historical ocean color data from the Coastal Zone Color Scanner (CZCS) has been used extensively to determine phytoplankton patterns and variability, characterize the primary production across the SCB coastal fronts, and describe the seasonal and interannual variability in pigment concentrations. As in CalCOFI, these studies described much of the observed structures and their variability over relatively large space and time scales.

Miller, Arthur J.↗
DOE OSTI2023

Power System Event Detection Using the Energy Detector: A Performance Analysis

As the grid becomes smarter, the need to accurately detect, predict, and classify waveform phenomena is growing. When it comes to detecting high-frequency behaviors, i.e. transients, it is especially important to employ an event detection system that is able to accurately uncover these types of disturbances that would otherwise be lost with traditional hardware. In this paper, we first present the energy detector; a waveform event detection system that adaptively monitors a signal's energy and picks out high-frequency events that deviate from the nominal state. Secondly, we evaluate the performance of this detector against waveform data that have been corrupted by sensor irregularities. Using Oak Ridge National Laboratory's sensor testbed, we are able to show the results of the detector's performance against events that have been corrupted by three distinct sensor types, and examine how these results change with multiple trials. The results show excellent performance when detecting the beginning of an anomalous event with an average of less than 1% error.

Ekti, Ali Riza↗
DOE OSTI2021

Asynchronous distributed-memory task-parallel algorithm for compressible flows on unstructured 3D Eulerian grids

Here, we discuss the implementation of a finite element method, used to numerically solve the Euler equations of compressible flows, using an asynchronous runtime system (RTS). The algorithm is implemented for distributed-memory machines, using stationary unstructured 3D meshes, combining data-, and task-parallelism on top of the Charm++ RTS. Charm++’s execution model is asynchronous by default, allowing arbitrary overlap of computation and communication. Task-parallelism allows scheduling parts of an algorithm independently of, or dependent on, each other. Built-in automatic load balancing enables continuous redistribution of computational load by migration of work units based on real-time CPU load measurement. The RTS also features automatic checkpointing, fault tolerance, resilience against hardware failure, and supports power-, and energy-aware computation. We demonstrate scalability up to 25 x 10 9 cells at $\mathscr{O}$10 4 compute cores and the benefits of automatic load balancing for irregular workloads. The full source code with documentation is available at https://quinoacomputing.org.

42 ENGINEERING↗
DOE OSTI2014

Reliable and Efficient Distributed Checkpointing System for Grid Environments

In Fine-Grained Cycle Sharing (FGCS) systems, machine owners voluntarily share their unused CPU cycles with guest jobs, as long as their performance degradation is tolerable. However, unpredictable evictions of guest jobs lead to fluctuating completion times. Checkpoint-recovery is an attractive mechanism for recovering from such “failures”. Today’s FGCS systems often use expensive, high-performance dedicated checkpoint servers. However, in geographically distributed clusters, this may incur high checkpoint transfer latencies. Here we present a distributed checkpointing system called FALCON that uses available disk resources of the FGCS machines as shared checkpoint repositories. However, an unavailable storage host may lead to loss of checkpoint data. Therefore, we model the failures of a storage host and develop a prediction algorithm for choosing reliable checkpoint repositories. We experiment with FALCON in the university-wide Condor testbed at Purdue and show improved and consistent performance for guest jobs in the presence of irregular resource availability.

checkpoint↗
DOE OSTI2021

Ensemble Learning, Prediction and Li-Ion Cell Charging Cycle Divergence

In recent years, the pervasive use of lithium ion (Li-ion) batteries in applications such as cell phones, laptop computers, electric vehicles, and grid energy storage systems has prompted the development of specialized battery management systems (BMS). The primary goal of a BMS is to maintain a reliable and safe battery power source while maximizing the calendar life and performance of the cells. To maintain safe operation, a BMS should be programmed to minimize degradation and prevent damage to a Li-ion cell, which can lead to thermal runaway. Cell damage can occur over time if a BMS is not properly configured to avoid overcharging and discharging. To prevent cell damage, efficient and accurate cell charging cycle characteristics algorithms must be employed. In this paper, computationally efficient and accurate ensemble learning algorithms capable of detecting Li-ion cell charging irregularities are described. Additionally, it is shown using machine and deep learning that it is possible to accurately and efficiently detect when a cell has experienced thermal and electrical stress due to cell overcharging by measuring charging cycle divergence.

25 ENERGY STORAGE↗
DOE OSTI2012

NONPARAMETRIC BAYESIAN ESTIMATION OF PERIODIC LIGHT CURVES

Many astronomical phenomena exhibit patterns that have periodic behavior. An important step when analyzing data from such processes is the problem of identifying the period: estimating the period of a periodic function based on noisy observations made at irregularly spaced time points. This problem is still a difficult challenge despite extensive study in different disciplines. This paper makes several contributions toward solving this problem. First, we present a nonparametric Bayesian model for period finding, based on Gaussian Processes (GPs), that does not make assumptions on the shape of the periodic function. As our experiments demonstrate, the new model leads to significantly better results in period estimation especially when the light curve does not exhibit sinusoidal shape. Second, we develop a new algorithm for parameter optimization for GP which is useful when the likelihood function is very sensitive to the parameters with numerous local minima, as in the case of period estimation. The algorithm combines gradient optimization with grid search and incorporates several mechanisms to overcome the high computational complexity of GP. Third, we develop a novel approach for using domain knowledge, in the form of a probabilistic generative model, and incorporate it into the period estimation algorithm. Experimental results validate our approach showing significant improvement over existing methods.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2023

Is the Global MHD Modeling of the Magnetosphere Adequate for GIC Prediction: the May 27–28, 2017 Storm

Practical steps taken by the international community to reduce the damage to technological systems from space weather include the development of numerical models capable of real-time predictions of electromagnetic disturbances at the Earth’s surface. Here we examine the feasibility of a version of the Space Weather Modeling Framework (SWMF) global MHD simulation code similar to that used by the NOAA Space Weather Prediction Center to predict the level of geomagnetic field variability, and consequently geomagnetically induced currents (GICs). We consider the contribution of geomagnetic disturbances to the bursts of GIC in the electric power line of the Kola Peninsula during the May 27–28, 2017 storm and compare the observations with results of the global MHD model. During the maximal disturbance magnetic field variations at East Scandinavian stations become more irregular, as intense Pi3 pulsations are superposed on the magnetic bay. These pulsations are not quasi-sinusoidal waves like typical Pc5 pulsations, but they are rather a quasi-periodic sequence of magnetic impulses with time scales ~5–15 min. During this period with elevated Pi3 activity very high values of GIC were recorded (variations >100 A) in the electric power transmission line. The SWMF modeling reasonably well reproduces the global magnetospheric parameters, such as SYM-H index or cross-polar potential. However, the magnetic field variability dB/dt in the East Scandinavia predicted by the modeling has turned out to be more than order of magnitude less than that observed. Thus, the version of SWMF with the grid used by NOAA SWPC still cannot adequately predict for the May 27–28 event the fine structure of the storm/substorm—Pi3 geomagnetic disturbances, and consequently the magnitude of the GIC that they drive.

geomagnetically induced currents↗
NASA NTRS1998

Virtual and Experimental Visualization of Flows in Packed Beds of Spheres Simulating Porous Media Flows

A videotape presentation of flow in a packed bed of spheres is provided. The flow experiment consisted of three principal elements: (1) an oil tunnel 76.2 mm by 76.2 mm in cross section, (2) a packed bed of spheres in regular and irregular arrays, and (3) a flow characterization methodology, either (a) full flow field tracking (FFFT) or (b) computational fluid dynamic (CFD) simulation. The refraction indices of the oil and the test array of spheres were closely matched, and the flow was seeded with aluminum oxide particles. Planar laser light provided a two-dimensional projection of the flow field, and a traverse simulated a three-dimensional image of the entire flow field. Light focusing and reflection rendered the spheres black, permitting visualization of the planar circular interfaces in both the axial and transverse directions. Flows were observed near the wall-sphere interface and within the set of spheres. The CFD model required that a representative section of a packed bed be formed and gridded, enclosing and cutting six spheres so that symmetry conditions could be imposed at all cross-boundaries. Simulations had to be made with the flow direction at right angles to that used in the experiments, however, to take advantage of flow symmetry. Careful attention to detail was required for proper gridding. The flow field was three-dimensional and complex to describe, yet the most prominent finding was flow threads, as computed in the representative 'cube' of spheres with face symmetry and conclusively demonstrated experimentally herein. Random packing and bed voids tended to disrupt the laminar flow, creating vortices.

Hendricks, R. C.↗