SEARCH · Engineering Papers

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4

All sources NASA NTRS DOE OSTI

DOE OSTI2026

The Physics Imposed on a Streaming Operator by Spherical Transport Problems

The streaming operator, which generates a displacement of a particle on a straight line at a constant speed in transport theory, is derived algebraically from a spherical coordinate formulation of Newton’s second law. This derivation leads to an operator that has more partial derivatives than a Cartesian coordinate formulation of the operator. The additional partial derivatives, which are with respect to the normalized velocity variables of a particle, take into account the intrinsic curvature of a ball. Moreover, these partial derivatives mitigate ray effects, which arise when a finite number of normalized velocities (also called directions or discrete ordinates) are used to simulate a continuous S 2 sphere of directions, by rotating the polar axis of the S 2 sphere into the radial direction of the coordinate system. As a consequence of this rotation, the number of actual discrete ordinates is greatly amplified to an enormous number of effective discrete ordinates by a multiplier that is equal to the number of patches that partitions a spherical surface. In addition to the derivation of the streaming operator, we provide in closed form a solution to the system of characteristic equations that is equivalent to the streaming operator. Furthermore, the solution to the system of characteristic equations enables the construction of an integral operator that is the inverse to the streaming operator. Examples in which ray effects are immensely mitigated by spherical coordinates are presented.

integral operator↗

DOE OSTI2020

Axisymmetric hydrodynamics in numerical relativity using a multipatch method

In this paper, we describe a method of implementing the axisymmetric evolution of general-relativistic hydrodynamics and magnetohydrodynamics through modification of a multipatch grid scheme. In order to ease the computational requirements required to evolve the post-merger phase of systems involving binary compact massive objects in numerical relativity, it is often beneficial to take advantage of these system's tendency to rapidly settle into states that are nearly axisymmetric, allowing for 2D evolution of secular timescales. We implement this scheme in the spectral Einstein code and show the results of application of this method to four test systems including viscosity, magnetic fields, and neutrino radiation transport. Our results show that this method can be used to quickly allow already existing 3D infrastructure that makes use of local coordinate system transformations to be made to run in axisymmetric 2D with the flexible grid creation capabilities of multipatch methods. Our code tests include a simple model of a binary neutron star postmerger remnant, for which we confirm the formation of a massive torus which is a promising source of post-merger ejecta.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Impact of FERC Order 2222 on DER Participation Rules in US Electricity Markets

Electricity markets in the bulk grid are beginning to implement market mechanisms that support the procurement of flexible capabilities from wide range of technologies, including distributed energy resources (DERs). The flexibility of these resources will help counterbalance supply uncertainties from large-scale integration of variable renewable generation. To encourage development of distributed and aggregated market participants, FERC Order 2222 was issued in September 2020 to require each Independent System Operator (ISO) in the US to implement rules that enable broader participation from aggregations of DERs in the bulk market. The following paper first describes the generic design of ISO markets before introducing the new market participation rules that ISOs have proposed for compliance with Order 2222. The paper then describes how software performance issues may continue to affect the eligibility requirements and offer structures for DER aggregations participating in ISOs, noting that continued research on computational methods may help reduce burdens for DER integration. The prospects for transmission and distribution system coordination is second major issue discussed, which will require minor changes to existing processes in the short term. In the longer term, there is more opportunity for more wide-ranging reforms, such as the development of a Distribution System Operator (DSO) framework. Newly proposed market rules may affect how Transactive Energy Systems (TES) will help facilitate efficient formation of DER aggregations and operation of the individual DERs within an aggregation. Within the TES context, the challenge is to fully understand how resource eligibility and operational and planning coordination methods will affect the design and implementation of TES.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2021

Diaboloidal mirrors: algebraic solution and surface shape approximations

A new type of optical element that can focus a cylindrical wave to a point focus (or vice versa) is analytically described. Such waves are, for example, produced in a beamline where light is collimated in one direction and then doubly focused by a single optic. A classical example in X-ray optics is the collimated two-crystal monochromator, with toroidal mirror refocusing. The element here replaces the toroid, and in such a system provides completely aberration free, point-to-point imaging of rays from the on-axis source point. We present an analytic solution for the mirror shape in its laboratory coordinate system with zero slope at the centre, and approximate solutions, based on bending an oblique circular cone and a bent right circular cylinder, that may facilitate fabrication and metrology.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2022

Curved lattices of crystals formed in glass

Atomic arrangements systemically deviating from the classic definition of single crystal lattice have been known in spherulites for over a century. However, their unusual curved lattice could not be exploited for useful applications due to their uncontrolled orientations and lattice curvatures. Heating of glass with focused laser, electron or x-ray beam has yielded well-defined conditions for forming single crystal architectures. Furthermore, by controlling the conditions of irradiation and confinement by glass matrix for a given system, desired lattice curvature (rotations and twists) can be reproducibly introduced into these crystals. Here, we review several examples of laser-fabricated crystals with curved lattices in glass using a common coordinate system, which helps to establish correlations between processing parameters and lattice curvature. Here, the mechanisms of lattice deformations are discussed along with potential methods to control the orientation of the initially formed seeds that determine the lattice orientations of the end product. The recent developments of laser-fabrication of single crystal architecture in glass with controlled heating profile appear promising for engineering the curvature of their lattices.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Endogenous Interface Pricing for Consistent Transmission–Distribution Co-Optimization With Discrete Distribution Controls

This paper proposes an endogenous interface pricing model for day-ahead transmission–distribution co-optimization that co-determines the interface locational marginal price (LMP) and the transmission–distribution exchange, ensuring price–dispatch consistency while optimally scheduling discrete distribution controls. The formulation couples a DC optimal power flow (OPF) with a branch-flow AC OPF that schedules distributed energy resources (DERs), tap-changer settings, capacitor banks (CBs), and multi-period energy storage systems (ESSs) under feeder voltage and current limits, and is solved as a mixed-integer second-order cone program (MISOCP). In a T14–D33 system, coordinated device scheduling recovers about 90% of the distribution-to-transmission export achievable in a reference case that ignores distribution network (DN) limits, while satisfying a 1.05 p.u. voltage upper bound. In a T39–D34/D37/D123 system, a sequential decoupled benchmark produces interface LMP distortions up to 12.5% and a 7.28% mismatch in net export energy, whereas the proposed model removes these distortions and the associated settlement mismatches. Second-order cone (SOC) relaxation gaps remain below $10^{-3}$ in all cases.

Noh, Seung-Gil↗

DOE OSTI2024

Systems Level Fuel Cycle Modeling in TMAP8 - A Demonstration

The tritium migration analysis program (TMAP) has been used for tritium inventory tracking and analysis for several years, and the Multiphysics Object Oriented Simulation Environment (MOOSE)-based TMAP8 has several improvements over TMAP4 and TMAP7, such as support for multiple dimensions and non-cartesian coordinate systems, as well as interoperability with sub-apps at higher and lower length scales. We demonstrate that TMAP8 has the additional capacity to solve systems-level problems using zero-dimensional ordinary differential equations by reproducing a literature model which describes a systems-level fuel-cycle of tritium inventory in a hypothetical fusion power plant. The capacity to run several coupled multi-scale physics calculations as part of a single package will be necessary for accurate blanket and fuel-cycle design.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

A Sheaf Theoretical Approach to Uncertainty Quantification of Heterogeneous Geolocation Information

Integration of multiple, heterogeneous sensors is a challenging problem across a range of applications. Prominent among these are multi-target tracking, where one must combine observations from different sensor types in a meaningful and efficient way to track multiple targets. Because different sensors have differing error models, we seek a theoretically justified quantification of the agreement among ensembles of sensors, both overall for a sensor collection, and also at a fine-grained level specifying pairwise and multi-way interactions among sensors. We demonstrate that the theory of mathematical sheaves provides a unified answer to this need, supporting both quantitative and qualitative data. Furthermore, the theory provides algorithms to globalize data across the network of deployed sensors, and to diagnose issues when the data do not globalize cleanly. We demonstrate and illustrate the utility of sheaf-based tracking models based on experimental data of a wild population of black bears in Asheville, North Carolina. A measurement model involving four sensors deployed among the bears and the team of scientists charged with tracking their location is deployed. This provides a sheaf-based integration model which is small enough to fully interpret, but of sufficient complexity to demonstrate the sheaf’s ability to recover a holistic picture of the locations and behaviors of both individual bears and the bear-human tracking system. A statistical approach was developed in parallel for comparison, a dynamic linear model which was estimated using a Kalman filter. This approach also recovered bear and human locations and sensor accuracies. When the observations are normalized into a common coordinate system, the structure of the dynamic linear observation model recapitulates the structure of the sheaf model, demonstrating the canonicity of the sheaf-based approach. However, when the observations are not so normalized, the sheaf model still remains valid.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Rigid-flexible coupling multi-body dynamics modeling of a semi-submersible floating offshore wind turbine

A 14-degree-of-freedom (14-DOF) flexible multibody dynamics model is developed and verified for a semi-submersible floating offshore wind turbine (FOWT). The model considers the coupled dynamics of the platform, tower, nacelle, blades, and mooring subjected to external wind and wave loads. The platform is simplified as a rigid body connected to the seabed by the mooring system. The tower and blade are modeled as flexible cantilever beams. The energy method is used to derive the governing equations of motion, where the kinetic energy, potential energy, and work done by external wind-wave forces are all deduced in a global coordinate system. The 5-MW baseline semi-submersible FOWT is used to verify the derived model against the results simulated from FAST developed by the US National Renewable Energy Laboratory (NREL) at two scenarios: free decay state and different wind-wave load cases. Results show that the established FOWT model can well reflect the vibration characteristics of FAST model. The application of the model to the control of platform pitch with a tuned mass damper is studied. Furthermore, the simplified model could provide a low-order method for the structural dynamics analysis and advanced vibration control design for the multi-body components of the semi-submersible FOWT in the future.

17 WIND ENERGY↗

DOE OSTI2020

Solving gyrokinetic systems with higher-order time dependence

We discuss theoretical and numerical aspects of gyrokinetics as a Lagrangian field theory when the field perturbation is introduced into the symplectic part. A consequence is that the field equations and particle equations of motion in general depend on the time derivatives of the field. The most well-known example is when the parallel vector potential is introduced as a perturbation, where a time derivative of the field arises only in the equations of motion, so an explicit equation for the fields may still be written. We will consider the conceptually more problematic case where the time-dependent fields appear in both the field equations and equations of motion, but where the additional term in the field equations is formally small. The conceptual issues were described by Burby (J. Plasma Phys., vol. 82 (3), 2016, 905820304): these terms lead to apparent additional degrees of freedom to the problem, so that the electric field now requires an initial condition, which is not required in low-frequency (Darwin) Vlasov–Maxwell equations. Also, the small terms in the Euler–Lagrange equations are a singular perturbation, and these two issues are interlinked. For well-behaved problems the apparent additional degrees of freedom are spurious, and the physically relevant solution may be directly identified. Because we needed to assume that the system is well behaved for small perturbations when deriving gyrokinetic theory, we must continue to assume that when solving it, and the physical solutions are thus the regular ones. The spurious nature of the singular degrees of freedom may also be seen by changing coordinate systems so the varying field appears only in the Hamiltonian. We then describe how methods appropriate for singular perturbation theory may be used to solve these asymptotic equations numerically. We then describe a proof-of-principle implementation of these methods for an electrostatic strong-flow gyrokinetic system; two basic test cases are presented to illustrate code functionality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

The Integrated Reference Region Analysis for Parallel DFIGs’ Interfacing Inductors

Although the traditional design of doubly-fed induction generators (DFIGs)’ interfacing inductors consider the peak ripple of the grid-side converter (GSC)’s output current, it does not consider the inductors’ impact on DFIGs’ smallsignal stability. Located in series between the GSC and the stator, the inappropriate selection of the interfacing inductor can easily result to system instability. Therefore, this paper first proposes the integrated reference region analysis for small wind farm's interfacing inductors to improve the traditional design method. Firstly, a new detailed parallel DFIGs’ smallsignal model that considers output currents’ coupling is built in d-q coordinate system with state-space approach. The model focuses on representing the operating states of parallel DFIGs in wind farm. Secondly, the linearized state-space matrix of parallel DFIGs is decomposed into nominal-value matrix and location matrix. Considering the traditional design requirements, the integrated reference region for interfacing inductor is proposed through spectral radius and bialternate matrix sum (BMS). Furthermore, it can provide better guidance for parameter selecting and stabilization method researches. Finally, the simulation and experimental results show that the proposed integrated reference region for parallel DFIGs’ interfacing inductors is accurate and instructive.

47 OTHER INSTRUMENTATION↗

DOE OSTI2026

Lethality is Local, but Survival is Systemic: Temporal and Multi-Organ Responses to Chlorine Gas Exposure in a Murine Model

Chlorine gas (Cl2) is a highly toxic chemical associated with both localized lung injury and systemic health effects. While pulmonary damage has been well characterized, the systemic inflammatory and metabolic responses remain poorly understood. We aimed to define the temporal and multi-organ responses to Cl2 exposure in a murine model, with a focus on identifying spatiotemporal inflammation and its impact on survival and lethality. SKH1 mice were exposed for 10 min to varying concentrations of Cl2 (94.4–810 ppm, representative of non-lethal, LD10, and LD50 doses) and monitored for respiratory function, perfusion, and acidosis using organ-specific imaging. At multiple time points (40 min, 6 h, 24 h, and 7 d), we measured phosphoproteins, cytokines, chemokines, growth factors, and metabolic hormones in the lungs, heart, cortex, and plasma. Statistical modeling and logistic regression were used to identify biomarkers associated with lethality and survival. We found that lung injury was the primary cause of potential lethality, particularly via early phosphoprotein signaling disruptions. However, survival correlated with early systemic coordination of inflammatory and metabolic signals across organs. Perfusion and acidosis imaging were strongly associated with chemokine and hormone responses. Key survival-associated plasma biomarkers included decreased insulin, increased ghrelin, and decreased eotaxin. While potential lethality from Cl2 exposure is locally driven by pulmonary injury, survival depends on systemic, multi-organ responses that occur rapidly post-exposure. Within this model, our findings identify a potential therapeutic window to enhance survival and suggest candidate biomarkers that may be explored translationally for both triage and treatment of chlorine-related incidents.

chlorine gas↗

DOE OSTI2020

Gaussian null coordinates, near-horizon geometry and conserved charges on the horizon of extremal nondilatonic black p-branes

In this paper, we examine the emergence of conserved charges on the horizon of a particular class of extremal nondilatonic black [Formula: see text]-branes (which reduce to extremal dilatonic black holes in [Formula: see text] dimensions upon toroidal compactification) in the presence of a probe massless scalar field in the bulk. This result is achieved by writing the black [Formula: see text]-brane geometry in a Gaussian null coordinate system which allows us to get a nonsingular near-horizon geometry description. We find that the near-horizon geometry is [Formula: see text] and that the [Formula: see text] section has an internal structure which can be seen as a warped product of [Formula: see text] in Gaussian null coordinates. We show that the bulk scalar field satisfying the field equations is expanded in terms of nonnormalizable and normalizable modes, which for certain suitable quantization conditions are well-behaved at the boundary of [Formula: see text] space. Furthermore, we show that picking the normalizable modes results in the existence of conserved quantities on the horizon. We discuss the impact of these conserved quantities in the late-time regime.

Astronomy & Astrophysics↗

DOE OSTI2020

Near-axis expansion of stellarator equilibrium at arbitrary order in the distance to the axis

A direct construction of equilibrium magnetic fields with toroidal topology at arbitrary order in the distance from the magnetic axis is carried out, yielding an analytical framework able to explore the landscape of possible magnetic flux surfaces in the vicinity of the axis. This framework can provide meaningful analytical insight into the character of high-aspect-ratio stellarator shapes, such as the dependence of the rotational transform and the plasma beta limit on geometrical properties of the resulting flux surfaces. The approach developed here is based on an asymptotic expansion on the inverse aspect ratio of the ideal magnetohydrodynamics equation. The analysis is simplified by using an orthogonal coordinate system relative to the Frenet–Serret frame at the magnetic axis. The magnetic field vector, the toroidal magnetic flux, the current density, the field line label and the rotational transform are derived at arbitrary order in the expansion parameter. Moreover, a comparison with a near-axis expansion formalism employing an inverse coordinate method based on Boozer coordinates (the so-called Garren–Boozer construction) is made, where both methods are shown to agree at lowest order. Finally, as a practical example, a numerical solution using a W7-X equilibrium is presented, and a comparison between the lowest-order solution and the W7-X magnetic field is performed.

Physics↗

DOE OSTI2021

Improved Prediction of Cold-Air Pools in the Weather Research and Forecasting Model Using a Truly Horizontal Diffusion Scheme for Potential Temperature

The terrain-following vertical coordinate system used by many atmospheric models, including the Weather Research and Forecasting (WRF) Model, is prone to errors in regions of complex terrain. These errors stem, in part, from the calculation of horizontal gradients within the diffusion term of the momentum or scalar evolution equations. In WRF, such gradients can be calculated along coordinate surfaces, or using metric terms that help account for grid skewness. However, neither of these options ensures a truly horizontal gradient calculation, especially if a grid cell is skewed enough that the heights of the neighboring grid points used in the calculation fall outside the vertical range of the cell. In this work, an improved scheme that uses Taylor series approximations to vertically interpolate variables to the level necessary for a truly horizontal gradient calculation is implemented in WRF for the diffusion of potential temperature. The scheme is validated using an atmosphere-at-rest configuration, in which spurious flows develop only as a result of numerical errors and can thus be used as a proxy for model performance. Following validation, the method is applied to the simulation of cold-air pools (CAPs), which occur in regions of complex terrain and are characterized by strong near-surface temperature gradients. Using the truly horizontal scheme, idealized simulations demonstrate reduced numerical mixing in a quiescent CAP, and a realistic case study in the Columbia River basin shows a reduction in positive wind speed bias by up to roughly 20% compared to observations from the Second Wind Forecast Improvement Project.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Continuum gyrokinetic simulations of edge plasmas in single-null geometries

The first continuum gyrokinetic calculations of electrostatic ion scale turbulence are presented for the case of a diverted tokamak geometry. The simulation model solves the long-wavelength limit of the full-F gyrokinetic equation for ion species coupled to the quasi-neutrality equation for electrostatic potential variations, where a fluid model is used for an electron response. Furthermore, in order to facilitate simulations of highly-anisotropic microturbulence in the presence of strong magnetic shear and a magnetic X-point, a numerical algorithm utilizing a locally field-aligned multiblock coordinate system has been developed and implemented in the 5D finite-volume code COGENT. In this approach, the toroidal direction is divided into blocks, such that within each block, the cells are field-aligned and a non-matching grid interface is allowed at block boundaries. The toroidal angle corresponds to the “coarse” field-aligned coordinate, whereas the poloidal cross section, comprised of the radial and poloidal directions, is finely gridded to resolve short-scale perpendicular turbulence structures and to support accurate re-mapping (interpolation) at block boundaries. The 5D simulations explore cross-separatrix ion scale turbulence in the presence of a self-consistent radial electric field and address the effects of magnetic-shear stabilization in the X-point region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

A conservative implicit-PIC scheme for the hybrid kinetic-ion fluid-electron plasma model on curvilinear meshes

We report that the hybrid kinetic-ion fluid-electron plasma model is widely used to study challenging multi-scale problems in space and laboratory plasma physics. Here, a novel conservative scheme for this model employing implicit particle-in-cell techniques is extended to arbitrary coordinate systems via curvilinear maps from logical to physical space. The scheme features a fully non-linear electromagnetic formulation with a multi-rate time advance - including sub-cycling and orbit-averaging for the kinetic ions. By careful choice of compatible particle-based kinetic-ion and mesh-based fluid-electron discretizations in curvilinear coordinates, as well as particle-mesh interpolations and implicit midpoint time advance, the scheme is proven to conserve total energy for arbitrary curvilinear meshes. In the electrostatic limit, the method is also proven to conserve total momentum for arbitrary curvilinear meshes. Although momentum is not conserved for arbitrary curvilinear meshes in the electromagnetic case, it is for an important subset of Cartesian tensor-packed meshes. The scheme and its novel conservation properties are demonstrated for several challenging numerical problems using different curvilinear meshes, including a merging flux-rope simulation for a space weather application, and a helical m = 1 mode simulation for magnetic fusion energy application.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

The absence of ray-effects in the discrete ordinate solution to the transport equation in spherical coordinates in multi-dimensions

The streaming operator of the transport equation is derived for spherical coordinates by starting from Newton’s second law for a free particle expressed in spherical coordinates. We shall show that the partial derivatives with respect to the velocity variables of the particle, which are absent in the Cartesian coordinate formulation of the transport equation, arise in the spherical coordinate formulation of the transport equation in response to the centrifugal force which prevents a free particle from ‘falling into the origin’ of the coordinate system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗