Search NASA⌕ Search

SEARCH · Search NASA

Results for “coordinate system”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Generalized Boozer coordinates: A natural coordinate system for quasisymmetry

We prove the existence of a straight-field-line coordinate system we call generalized Boozer coordinates. This coordinate system exists for magnetic fields with nested toroidal flux surfaces provided ∮ d l / B   ( j · ∇ ψ ) = 0 , where symbols have their usual meaning, and the integral is taken along closed magnetic field lines. All quasisymmetric fields, regardless of their associated form of equilibria, must satisfy this condition. This coordinate system presents itself as a convenient form to describe general quasisymmetric configurations and their properties. Insight can be gained analytically into the difference between strong and weak forms of quasisymmetry, as well as axisymmetry, and the interaction of quasisymmetry with different forms of equilibria.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Constructing field-aligned coordinate systems for gyrokinetic simulations of tokamaks in X-point geometries

Structures in tokamak plasmas are elongated along the direction of the magnetic field and short in the directions perpendicular to the magnetic field. Many tokamak simulation codes take advantage of this by using a field-aligned coordinate system. However, field-aligned coordinate systems have a coordinate singularity at magnetic X-points where the poloidal magnetic field vanishes, which makes it difficult to use field-aligned coordinate systems when simulating the core and scrape-off layer simultaneously. Here, we present an algorithm for grid generation and computing geometric quantities in a standard field-aligned coordinate system that avoids the singularity and allows one to conduct two-dimensional gyrokinetic axisymmetric simulations in X-point geometries. Convergence tests of advection, boundary value problems and geometric quantities all show greater than first-order convergence even in the vicinity of the X-point. We also demonstrate the geometric consistency of our algorithm with an example simulation of the spherical tokamak for energy production, which shows machine-precision particle conservation.

fusion plasma↗
DOE OSTI2022

Pitch-Angle Diffusion in the Earth’s Magnetosphere Organized by the Mozer-Transformed Coordinate System

In the Earth’s dipole magnetosphere finite-gyroradius effects produce a shift of the atmospheric loss cone away from the direction of the magnetic field. This loss-cone shift is theoretically described by the “Mozer transform”, which is based upon the curvature drift of particles crossing the equatorial plane. For positive ions the northern and southern loss cones both shift westward and for electrons the northern and southern loss cones both shift eastward. This loss-cone shift is part of a coordinate-system transform, with the transformed coordinates better organizing the behavior of particle orbits in the dipole magnetic field (e.g. first adiabatic invariants, mirror heights, and bounce times). In this report it is demonstrated that the transformed coordinate system also properly organizes pitch-angle diffusion. This improved organization of the diffusion is true whether the angular scattering is produced by plasma-wave scattering or by field-line-curvature (FLC) scattering. It is shown that FLC scattering and the loss cone shift are linked, so that if FLC scattering is occurring, there is a loss cone shifted away from the magnetic-field direction and the Mozer-transformed coordinates are needed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Explicit structure-preserving geometric particle-in-cell algorithm in curvilinear orthogonal coordinate systems and its applications to whole-device 6D kinetic simulations of tokamak physics

Explicit structure-preserving geometric particle-in-cell (PIC) algorithm in curvilinear orthogonal coordinate systems is developed. The work reported represents a further development of the structure-preserving geometric PIC algorithm achieving the goal of practical applications in magnetic fusion research. The algorithm is constructed by discretizing the field theory for the system of charged particles and electromagnetic field using Whitney forms, discrete exterior calculus, and explicit non-canonical symplectic integration. In addition to the truncated infinitely dimensional symplectic structure, the algorithm preserves exactly many important physical symmetries and conservation laws, such as local energy conservation, gauge symmetry and the corresponding local charge conservation. As a result, the algorithm possesses the long-term accuracy and fidelity required for first-principles-based simulations of the multiscale tokamak physics. The algorithm has been implemented in the SymPIC code, which is designed for high-efficiency massively-parallel PIC simulations in modern clusters. The code has been applied to carry out whole-device 6D kinetic simulation studies of tokamak physics. A self-consistent kinetic steady state for fusion plasma in the tokamak geometry is numerically found with a predominately diagonal and anisotropic pressure tensor. The state also admits a steady-state sub-sonic ion flow in the range of 10 km s -1 , agreeing with experimental observations and analytical calculations Kinetic ballooning instability in the self-consistent kinetic steady state is simulated. It is shown that high-n ballooning modes have larger growth rates than low-n global modes, and in the nonlinear phase the modes saturate approximately in 5 ion transit times at the 2% level by the E × B flow generated by the instability. These results are consistent with early and recent electromagnetic gyrokinetic simulations.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

QUIC-URB and QUIC-fire extension to complex terrain: Development of a terrain-following coordinate system

Ensemble-based approaches to prescribed fire planning cannot be supported by CFD-based models like FIRETEC and WFDS because they are too computationally expensive and cannot leverage LES approaches like CAWFE and WRF-SFIRE because too coarse of resolution. QUIC-Fire was developed to fill this gap but it cannot currently address complex terrain, typical for instance of the Western United States. In this paper, we describe the extension of the diagnostic wind model QUIC-URB, the wind engine of QUIC-Fire, to a terrain-following coordinate system. In particular, the paper presents the mathematical derivation of the wind solver leading to a linear system of equations that are solved through the successive over-relaxation method. The model is validated against a standard test used in previous works (the Askervein Hill) and against a new dataset from measurements in the Socorro Mountains, New Mexico. The terrain-following implementation captured the correct phenomenology for the isolated Askervein Hill, with a wind speed up at the top of the hill. We report the model agreed well with measurements on the upwind side of the peak, but overestimated speed-up on the downwind side of the hill. This is due to the inability of the model to generate flow separation and wake-eddy dynamics. On a common laptop, the divergence-free wind field was obtained in 6 s, making the solver appealing for coupled fire–atmosphere simulations. The Socorro Mountain was highly complex, with many cliff faces, peaks, and valleys. Although the model captures the magnitude and direction of inlet and outlet areas of the domain, it performs rather poorly in the valley region and in the regions near the steep cliffs. Hence, the model shows good agreement with data in areas of open sloped terrain but lacks in areas where flow separation and thermally driven effects may be present (neither effect was addressed in this work). Results highlight that future work should focus on the implementation of parameterizations of wake-eddies, similar to QUIC-URB’s building parameterizations, and on thermodynamic-driven flow.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Texture measurements on quartz single crystals to validate coordinate systems for neutron time-of-flight texture analysis

In crystallographic texture analysis, ensuring that sample directions are preserved from experiment to the resulting orientation distribution is crucial to obtain physical meaning from diffraction data. This work details a procedure to ensure instrument and sample coordinates are consistent when analyzing diffraction data with a Rietveld refinement using the texture analysis software MAUD . A quartz crystal is measured on the HIPPO diffractometer at Los Alamos National Laboratory for this purpose. The methods described here can be applied to any diffraction instrument measuring orientation distributions in polycrystalline materials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Data-Driven Distribution System Coordinated PV Inverter Control Using Deep Reinforcement Learning

The deployment of distributed solar photovoltaic (PV) systems has increased consistently over the past decades. High penetrations of PVs could cause a series of adverse grid impacts, such as voltage violations. The recent development of smart inverter technologies rises the incentives of developing PV control solutions that regulate the inverter output power and seeking the optimization on system operational objectives. This paper proposes a data-driven control solution based on deep reinforcement learning (DRL) to optimize PV inverters for voltage regulation. The proposed solution can minimize PV real power curtailment while maintaining network voltage at an acceptable range. Comparison results between the proposed DRL control algorithms with deep deterministic policy gradient (DDPG) and volt-var control on a real feeder in west Colorado highlight the advantage of the proposed framework in controlling the system voltage while minimizing the PV real power curtailment.

deep reinforcement learning↗
DOE OSTI2022

Control system for multi-system coordination via a single reference governor

This report describes the improvements to the Feasible Actuator Range Modifier (FARM) component of the RAVEN-based HYBRID framework for analysis of Integrated Energy Systems (IES). FARM supports the HERON plug-in that solves the power dispatch problem. The solution involves economically optimal dispatches that satisfy the limits on production variables and corresponding rates of variation (explicit constraints) as well as the limits on the process variables tied to the service life of equipment (implicit constraints). The problem can be addressed as a two-stage process, i.e., HERON power dispatcher estimates a solution that meets explicit constraints (low-resolution physics), whereas FARM uses the simulation outcomes of HYBRID high-fidelity model to capture the system dynamic response and enforce implicit constraints (high-resolution physics). The initial version of the code (FARM-Alpha) was released by Argonne National Laboratory in January 2021 followed by FARM-Beta (January 2022) and FARM-Gamma (April 2022).

42 ENGINEERING↗
DOE OSTI2025

Secure State Estimation with Asynchronous Measurements for Coordinated Cyber Attack Detection in Active Distribution Systems

Coordinated cyber attacks tamper with measurement data to disrupt the situational awareness of active distribution systems. Various sensors report measurements asynchronously at different rates, which introduces challenges during state estimation. In addition, this forces cyber intruders to exert greater effort to compromise multiple communication channels and launch coordinated attacks. Therefore, multi-channel and asynchronous measurements could be harnessed to develop more secure cyber defense strategies. In this paper, a prediction-correction-based multi-rate observer is designed to exploit the value of asynchronous measurements for the detection of coordinated false data injection (FDI) attacks. First, a time-function-dependent prediction-correction strategy is proposed to adjust the sampling interval for each sensor’s measurement. Then, an observer is designed based on the trade-off between estimation error and the optimal period of the most recent sampling instant, with the convergence of estimation error with the maximum permitted sampling interval. Moreover, the conditions for exponential stability are developed using the Lyapunov–Krasovskii functional technique. Next, a coordinated FDI attack detection strategy is developed based on the dual nonlinear minimization problem. The proposed attack detection and secure state estimation strategies are tested on the IEEE 13-node system. Simulation results show that these schemes are effective in enhancing attack detection based on asynchronous measurements or compromised data.

asynchronous measurements↗
DOE OSTI2024

Model-Based Detection of Coordinated Attacks (DCA) in Distribution Systems

The fast-paced growth in digitization of smart grid components enhances system observability and remote-control capabilities through efficient communication. However, enhanced connectivity results in heightened system vulnerability towards cybersecurity risks in the cyber-physical power system. Coordinated cyber-attacks (CCA), when undetected, lead to system-wide impact in terms of large disturbances or widespread outages. Detecting CCA in the cyber layer is critical to thwart cyber-attacks in real-time before the attack impacts the physical system. The challenge of locating CCA stems from the complex grid dynamics, making it difficult to distinguish between normal operational variations and cyber-attack impact. CCA often employs multiple attack vectors targeting geographically distributed components, further complicating CCA identification. Existing research in intrusion detection is primarily focused on the transmission network and limited to detecting individual attacks. In this paper, a novel proactive DCA strategy is proposed for early detection of CCA by establishing correlations among distinct attack events through model-based reinforcement learning that utilizes abductive reasoning to conclude the attacker goal. The solution includes understanding the system model, learning the system dynamics, and correlating individual cyber-attacks to extract the attacker’s objective. The developed learning algorithm identifies the most probable attack path to reach the attacker’s objective by predicting the next attack steps. A DNP3-based cyber-physical co-simulation testbed is developed to test the proposed algorithm using the IEEE 13-node test feeder.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Tidal Energy Resource Characterization, Velocity and Turbulence Measurements, Processed Data, Cook Inlet, AK, 2021

This submission contains processed datasets from a long-term deployment of 3 moorings and a transect survey of the proposed tidal energy site off the East Forelands in Cook Inlet, AK. The long-term mooring datasets were created from 8 instruments mounted on a Terrasond High Energy Oceanographic Mooring (THEOM) bottom lander and two Mid-Water Mooring (MWM) Stablemoor buoys from 1 July 2021 to 31 August 2021 (60 days). The west-most mooring (MWM1) was deployed at 60.720225 N, 151.436196 W in ~50 m of water. The middle mooring (THEOM) was deployed at 60.720703 N, 151.429500 W in ~52 m of water. The east-most buoy (MWM2) was deployed at 60.720081 N, 151.420896 W in ~50 m of water. Each Stablemoor carried three instruments: 1. A Nortek Vector acoustic Doppler velocimeter (ADV) mounted at the Stablemoor's nose. Data were recorded at 8 Hz on a 5 minute duty cycle every 20 minutes. Data was motion-corrected using the internal IMU and external ADCP bottom-track data and then bin-averaged into 4 minute bins and converted to the Principal (streamwise, cross-stream, vertical) coordinate system. (Note: 30 seconds were trimmed from the beginning and end of each 5 minute duty cycle to account for the filter end-effects from turning on and turning off the IMU.) 2. A down-looking Nortek Signature 1000 kHz acoustic Doppler current profiler (ADCP) mounted in the first Stablemoor instrument well. Data were recorded in 2 Hz with 5-beam burst and bottom-track enabled. Processed data has been averaged into 10 minute bins and converted into the Principal coordinate system. 3. An up-looking Nortek Signature 1000 kHz acoustic Doppler current profiler (ADCP) mounted in the second Stablemoor instrument well. Data were recorded at 4 Hz with 5 beam burst enabled. Processed data has been averaged into 10 minute bins and converted into the Principal coordinate system. Note: the down-facing ADCP on MWM1 failed on July 10th, 2021, only recording 9 days of data. Because ADV motion-correction required bottom track, the ADV from MWM1 also only has 9 days processed. Additionally, only 25 days of data were processed from the MWM2 ADV because it appeared to have been impacted by debris on 7/25. Two instruments were mounted on the THEOM (see MHKDR link further below for THEOM raw data): 4. A Nortek Vector acoustic Doppler velocimeter (ADV). Data were recorded at 8 Hz on a 5 minute duty cycle every 20 minutes. Data was bin-averaged into 5 minute bins, and converted to the Principal coordinate system. 5. A Nortek Signature 500 kHz acoustic Doppler current profiler (ADCP). Data were recorded in 4 Hz in the beam coordinate system from all 5 beams. Processed data has been averaged into 10 minutes bins and converted to the Principal coordinate system.

16 TIDAL AND WAVE POWER↗
DOE OSTI2023

How to Flatten the Earth?

Since the passage of the Geospatial Data Act (GDA) in 2018, federal agencies, including the Department of Energy (DOE), have been working to ensure compliance and outline guidance for the geospatial datasets and affiliated products publicly released from their agencies. At DOE this sparked work to develop a geospatial data management strategy and implementation plan; both of which spotlighted the need for the agency to set preferred guidance for the agency to use to inform the release of any public facing geospatial data and products. One of the immediate needs was set around establishing preferred guidance for coordinate systems, given the foundational role they play in all geospatial data, as they help reference data to the earth. DOE’s Geospatial Information Officer (GIO) and the Geospatial Sciences Steering Committee (GSSC) worked together to summarize best practices and current DOE uses to outline recommendations for coordinate systems to use for public geospatial data and products. This poster summarizes the role of coordinate systems, trends in coordinate system use from across the Department that were used to set preferred guidance for the Department to consider when releasing public-facing geospatial data and products.

Bauer, Jennifer↗
DOE OSTI2021

Tidal Energy Resource Characterization, Bottom Lander Measurements, Cook Inlet, AK, 2021

These datasets are from tidal resource characterization measurements collected on the Terrasond High Energy Oceanographic Mooring (THEOM) from 1 July 2021 to 30 August 2021 (60 days) in Cook Inlet, Alaska. The lander was deployed at 60.7207031 N, 151.4294998 W in ~50 m of water. The dataset contains raw and processed data from the following two instruments: 1. A Nortek Signature 500 kHz acoustic Doppler current profiler (ADCP). Data were recorded in 4 Hz in the beam coordinate system from all 5 beams. Processed data has been averaged into 5 minutes bins and converted to the East-North-Up (ENU) coordinate system. 2. A Nortek Vector acoustic Doppler velocimeter (ADV). Data were recorded at 8 Hz in the beam coordinate system. Processed data has been averaged into 5 minutes bins and converted to the Streamwise - Cross-stream - Vertical (Principal) coordinate system. Turbulence statistics were calculated from 5-minute bins, with an FFT length equal to the bin length, and saved in the processed dataset. Data was read and analyzed using the DOLfYN (version 1.0.2) python package and saved in MATLAB (.mat) and netCDF (.nc) file formats. Files containing analyzed data (".b1") were standardized using the TSDAT (version 0.4.2) python package. NetCDF files can be opened using DOLfYN (e.g., `dat = dolfyn.load(''*.nc")`) or the xarray python package (e.g. `dat = xarray.open_dataset("*.nc"). All distances are in meters (e.g., depth, range, etc), and all velocities in m/s. See the DOLfYN documentation linked in the submission, and/or the Nortek documentation for additional details.

16 TIDAL AND WAVE POWER↗
DOE OSTI2022

Sharing tuples across independent coordination namespace systems

A system and method for federating a tuple storage database across multiple coordinated namespace (CNS) extended memory storage systems allowing the sharing of tuples and tuple data across independent systems. The method provides a federation service for multiple coordination namespace systems. The method retrieves a tuple from connected independent CNS systems wherein a local CNS Controller sends a read request to the local gatekeeper to retrieve a first tuple and creates a local pending remote record. The local gatekeeper at a requesting node sends a broadcast query to a plurality of remote gatekeepers for the tuple and Remote gatekeepers at remote nodes query in its local CNS for the tuple. The Local gatekeeper process at the requesting node receives results from a plurality of remote gatekeepers for the said tuple and selects one remote gatekeeper to receive the requested tuple and broadcasts a read for tuple data with selected gatekeeper.

Jacob, Philip↗
DOE OSTI2023

Building-to-Grid Modeling Framework with a Case Study of Battery Systems

Coordinating electricity load shifting operations of multiple buildings has the potential to achieve significant benefits for both building owners and grid operation. In support of this research, this paper develops a novel building-to-grid modeling framework, which uses EnergyPlus to develop physics-based building energy models and GridLAB-D to develop distribution-level grid models. The data exchanges between buildings and distribution-level grid are conducted using Python script and HELICS co-simulation platform. Then, by using this modeling framework, this paper conducts a case study of battery systems in a connected community, which consists of 10 multi-family buildings with appropriately sized batteries. The three scenarios, baseline operation, tolerant to price fluctuations, and opportune usage, are studied. The result displays that, by adopting onsite batteries, both electricity costs for building owners and aggregated peak power demands are reduced.

Building-to-grid, connected community, co-simulati↗
DOE OSTI2026

PAVplotR

PAVplotR is an R package that accepts as input a haploid VCF and, for a user-specified bin size, calculates proportional presence-absence variation (PAV). First, all variants in the file are considered to calculate an expanded coordinate system that accounts for insertions relative to the reference sequence. Next, for each sample in the VCF, proportional PAV is summarized per bin in the expanded coordinate system and finally, PAV is plotted such that rows are samples and columns are bins with the x-axis (chromosome position) labeled according to the reference haplotype's original coordinate system.

Harder, Avril↗