Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Coordinate transformations”

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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Faster-than-real-time Simulation with Demonstration for Resilient DER Integration

The US electric grid is facing operational, stability, and security challenges. Transmission system operators need some measure of visibility into distribution system renewable generation. Distribution system generation needs to support transmission system voltage. The grid is experiencing an expansion in measurement systems. How to take full advantage of this expansion and defend against attacks, both cyber and physical, poses additional challenges. The Faster-than-real-time Simulation with demonstration for Resilient DER Integration project set out to do the following: a. Flatten the voltage profile through the feeders and system for cost saving and voltage stabilization needs. b. Increase the amount of intermittent distributed energy resources (IDERs) that could be deployed on a utility feeder and provide 100% or more energy needed for the demands on that feeder, and c. based on an accurate model (Digital Twin) of the utilities system, be able to detect any abnormalities on the utilities distribution system. To manage the voltage and increase IDER penetration (a,b), Graph Trace Analysis is employed in a time-series, optimal power flow to coordinate the time-varying feedback control setpoints of a distribution feeder’s utility control devices. Under the coordinated control are a Load Tap Changing Transformer, a voltage regulator, and five switched capacitor banks. The feeder serves over 2000 customers, the feeder secondaries are modeled, and the feeder has 2.3 MW of PV generation, corresponding to a 17.4% penetration of PV generation. The feeder model has over 12,000 components, where every customer load bus and PV generator are modeled. The accuracy of the power flow solution is compared against historical meter voltage measurements, the improvement in conservation voltage reduction energy savings as a function of the coordinated control desired voltage profile is investigated, and the increase in PV penetration of the coordinated control over the existing control is presented. To achieve improved control performance while observing system operation constraints, bellwether Advanced Metering Infrastructure (AMI) voltage measurements are used to adjust the desired voltage profile used by the optimal power flow analysis. To detect and alleviate or negate attacks or failures on the distribution and transmission utility grids (c) the grid needs to be resilient and self-healing. In this project software was designed to do just that. At the center of the software is an Integrated System Model (ISM) that spans from transmission to secondary distribution. The ISM is employed in real-time abnormality detection, voltage stability forecasting, and multi-mode control. Testing results are presented for: 1—attacks on utility infrastructure; 2—energy savings from optimal control; 3—distribution system control response during a low voltage transmission system event; 4—cyber-attacks on PV inverters, where physical inverters are used in hard-ware-in-the-simulation-loop studies. Contributions of this work include real-time analysis that spans from three-phase transmission through secondary distribution; an approach for detecting abnormalities that employs measurements from three independent measurement systems; and a multi-mode distribution system control that responds to cyber-attacks, physical attacks, equipment failures, and transmission system needs.

Integrated System Model, Graph Trace Analysis, Adv↗
NASA NTRS1993

Use of Fourier transforms for asynoptic mapping: Applications to the Upper Atmosphere Research Satellite microwave limb sounder

Fourier analysis has been applied to data obtained from limb viewing instruments on the Upper Atmosphere Research Satellite. A coordinate system rotation facilitates the efficient computation of Fourier transforms in the temporal and longitudinal domains. Fields such as ozone (O3), chlorine monoxide (ClO), temperature, and water vapor have been transformed by this process. The transforms have been inverted to provide maps of these quantities at selected times, providing a method of accurate time interpolation. Maps obtained by this process show evidence of both horizontal and vertical transport of important trace species such as O3 and ClO. An examination of the polar regions indicates that large-scale planetary variations are likely to play a significant role in transporting midstratospheric O3 into the polar regions. There is also evidence that downward transport occurs, providing a means of moving O3 into the polar vortex at lower altitudes. The transforms themselves show the structure and propagation characteristics of wave variations.

Elson, Lee S.↗
NASA NTRS1980

A three-dimensional body-fitted coordinate system for flow field calculations on asymmetric nosetips

A three dimensional body-fitted coordinate system developed for use in the calculation of inviscid flows over ablated, asymmetric reentry vehicle nosetips is described. Because of the potential geometric asymmetries, no standard coordinate system (e.g., spherical, axisymmetric reference surface-normal) is capable of being closely aligned with the nosetip surface. To generate a 3-D, body-fitted coordinate system an analytic mapping procedure is applied that is conformal within each meridional plane of the nosetip; these transformations are then coupled circumferentially to yield a three dimensional coordinate system. The mappings used are defined in terms of hinge points, which are points selected to approximate the body contours in each meridional plane. The selection of appropriate hinge points was automated to facilitate the use of the resulting nosetip flow field code.

Hall, D. W.↗
NASA NTRS1974

Finite element solution theory for three-dimensional boundary flows

A finite element algorithm is derived for the numerical solution of a three-dimensional flow field described by a system of initial-valued, elliptic boundary value partial differential equations. The familiar three-dimensional boundary layer equations belong to this description when diffusional processes in only one coordinate direction are important. The finite element algorithm transforms the original description into large order systems of ordinary differential equations written for the dependent variables discretized at node points of an arbitrarily irregular computational lattice. The generalized elliptic boundary conditions is piecewise valid for each dependent variable on boundaries that need not explicitly coincide with coordinate surfaces. Solutions for sample problems in laminar and turbulent boundary flows illustrate favorable solution accuracy, convergence, and versatility.

Baker, A. J.↗
DOE OSTI2025

Extending the Mott–Gurney law to one-dimensional nonplanar diodes using point transformations

Recent studies have applied variational calculus, conformal mapping, and point transformations to generalize the one-dimensional (1D) space-charge limited current density (SCLCD) and electron emission mechanisms to nonplanar geometries; however, these assessments have focused on extending the Child–Langmuir law (CLL) for SCLCD in vacuum. Since the charge in the diode is independent of the coordinate system (i.e., covariant), we apply bijective point transformations to extend the Mott–Gurney law (MGL) for the SCLCD in a collisional or semiconductor gap to nonplanar 1D geometries. This yields a modified MGL that replaces the Cartesian gap distance with a canonical gap distance that may be written generally in terms of geometric scale factors that are known for multiple geometries. We tabulate results for common geometries. Such an approach may be applied to any current density, including non-space-charge limited gaps and SCLCD that may fall between the CLL and MGL.

Garner, Allen L. [Purdue Univ., West Lafayette, IN↗
DOE OSTI2024

Solution of the Schrödinger equation for quasi-one-dimensional materials using helical waves

We formulate and implement a spectral method for solving the Schrödinger equation, as it applies to quasi-one-dimensional materials and structures. This allows for computation of the electronic structure of important technological materials such as nanotubes (of arbitrary chirality), nanowires, nanoribbons, chiral nanoassemblies, nanosprings and nanocoils, in an accurate, efficient and systematic manner. Our work is motivated by the observation that one of the most successful methods for carrying out electronic structure calculations of bulk/crystalline systems — the plane-wave method — is a spectral method based on eigenfunction expansion. Our scheme avoids computationally onerous approximations involving periodic supercells often employed in conventional plane-wave calculations of quasi-one-dimensional materials, and also overcomes several limitations of other discretization strategies, e.g., those based on finite differences and atomic orbitals. The basis functions in our method — called helical waves (or twisted waves) — are eigenfunctions of the Laplacian with symmetry adapted boundary conditions, and are expressible in terms of plane waves and Bessel functions in helical coordinates. We describe the setup of fast transforms to carry out discretization of the governing equations using our basis set, and the use of matrix-free iterative diagonalization to obtain the electronic eigenstates. Miscellaneous computational details, including the choice of eigensolvers, use of a preconditioning scheme, evaluation of oscillatory radial integrals and the imposition of a kinetic energy cutoff are discussed. We have implemented these strategies into a computational package called HelicES (Helical Electronic Structure). We demonstrate the utility of our method in carrying out systematic electronic structure calculations of various quasi-one-dimensional materials through numerous examples involving nanotubes, nanoribbons and nanowires. We also explore the convergence properties of our method, and assess its accuracy and computational efficiency by comparison against reference finite difference, transfer matrix method and plane-wave results. We anticipate that our method will find applications in computational nanomechanics and multiscale modeling, for carrying out transport calculations of interest to the field of semiconductor devices, and for the discovery of novel chiral phases of matter that are of relevance to the burgeoning quantum hardware industry.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Coordination environment tuning of nickel sites by oxyanions to optimize methanol electro-oxidation activity

To achieve zero-carbon economy, advanced anode catalysts are desirable for hydrogen production and biomass upgrading powered by renewable energy. Ni-based non-precious electrocatalysts are considered as potential candidates because of intrinsic redox attributes, but in-depth understanding and rational design of Ni site coordination still remain challenging. Here, we perform anodic electrochemical oxidation of Ni-metalloids (NiP x , NiS x , and NiSe x ) to in-situ construct different oxyanion-coordinated amorphous nickel oxyhydroxides (NiOOH-TO x ), among which NiOOH-PO x shows optimal local coordination environment and boosts electrocatalytic activity of Ni sites towards selective oxidation of methanol to formate. Experiments and theoretical results demonstrate that NiOOH-PO x possesses improved adsorption of OH* and methanol, and favors the formation of CH 3 O* intermediates. The coordinated phosphate oxyanions effectively tailor the d band center of Ni sites and increases Ni-O covalency, promoting the catalytic activity. This study provides additional insights into modulation of active-center coordination environment via oxyanions for organic molecules transformation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Theory, observations, and simulations of kinetic entropy in a magnetotail electron diffusion region

In this work, we examine velocity-space kinetic entropy, a spatially local measure of entropy for systems out of thermal equilibrium, during an encounter of an electron diffusion region at a magnetic reconnection site in Earth’s magnetotail by the Magnetospheric Multiscale (MMS) mission. We start by generalizing the theory of kinetic entropy to the case of non-uniform velocity space grids and transforming the equations into spheri- cal energy coordinates useful to experimental plasma detectors. The theory is then applied to MMS data and compared to particle-in-cell simulations of reconnection. We demonstrate that the entropy-based non-Maxwellianity measure from the MMS data is of sufficiently high precision to reliably identify non-Maxwellian distributions and therefore the measurements when kinetic effects are most significant. By comparing two different non-Maxwellian measures, we show that total entropy density suffers from “information loss” because it lacks a dependence on the velocity space grid, and so has lost information about how well a distribution function is resolved. Local velocity-space kinetic entropy density recovers this information. We quantify information loss and argue that the considerations needed to minimize it are crucial for instruments designed to measure distribution functions in situ.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Atomic protein structure refinement using all-atom graph representations and SE(3)-equivariant graph transformer

Abstract Motivation The state-of-art protein structure prediction methods such as AlphaFold are being widely used to predict structures of uncharacterized proteins in biomedical research. There is a significant need to further improve the quality and nativeness of the predicted structures to enhance their usability. In this work, we develop ATOMRefine, a deep learning-based, end-to-end, all-atom protein structural model refinement method. It uses a SE(3)-equivariant graph transformer network to directly refine protein atomic coordinates in a predicted tertiary structure represented as a molecular graph. Results The method is first trained and tested on the structural models in AlphaFoldDB whose experimental structures are known, and then blindly tested on 69 CASP14 regular targets and 7 CASP14 refinement targets. ATOMRefine improves the quality of both backbone atoms and all-atom conformation of the initial structural models generated by AlphaFold. It also performs better than two state-of-the-art refinement methods in multiple evaluation metrics including an all-atom model quality score—the MolProbity score based on the analysis of all-atom contacts, bond length, atom clashes, torsion angles, and side-chain rotamers. As ATOMRefine can refine a protein structure quickly, it provides a viable, fast solution for improving protein geometry and fixing structural errors of predicted structures through direct coordinate refinement. Availability and implementation The source code of ATOMRefine is available in the GitHub repository (https://github.com/BioinfoMachineLearning/ATOMRefine). All the required data for training and testing are available at https://doi.org/10.5281/zenodo.6944368.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Doppler Lidar Motion-Correction (DLMC) Value-Added Product Report

The U.S. Department of Energy Atmospheric Radiation Measurement (ARM) second ARM Mobile Facility (AMF2) Doppler lidar (S/N 0319-160) was deployed on the German ice breaker Polarstern during the Multidisciplinary Drifting Observatory for the Study of Arctic Climate (MOSAiC) campaign during 2019-2020 (see Figure 1). This was the first deployment of the new AMF2 Doppler lidar, as well as the first ship-based deployment of a Doppler lidar (DL) by ARM. In contrast to land-based deployments, the lidar’s heading (a.k.a. home point) and tilt are constantly changing, and lidar’s radial velocity measurements are impacted by the ship’s motion. Since the beam directions are reported relative to the instrument’s frame of reference, derivation of higher-order data products such as wind speed and direction require that the lidar’s attitude and translational velocity be properly accounted for. The Doppler Lidar Motion-Correction (DLMC) Value-Added-Product (VAP) was developed specifically for ship-based deployments of the Doppler lidar. This VAP combines raw uncorrected data from the DL and simultaneous measurements from the ARM Navigation system (NAV) (Walton 2019) to transform the beam angles from the lidar coordinate system to an Earth-fixed coordinate system. The VAP also removes the contribution of the lidar’s platform velocity from the radial velocity measurements. This report documents the methods used by the DLMC to perform these corrections.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1972

Analysis of heat transfer in a porous cooled wall with variable pressure and temperature along the coolant exit boundary

Fluid from a reservior at constant pressure and temperature is forced through a porous wall of uniform thickness. The boundary through which the fluid exits has specified variations in pressure and temperature along it in one direction so that the flow and heat transfer are two-dimensional. The local fluid and matrix temperatures are assumed to be equal and therefore a single energy equation governs the temperature distribution within the wall. The solution is obtained by transforming this energy equation into potential plane coordinates, which results in a separable equation. A technique yielding an integral equation is used to adapt the general solution so that it satisfies the variable-pressure boundary condition. Analytical expressions are given for the normal exit velocity and heat flux along the exit boundary. Illustrative examples are carried out which indicate to what extent the solution is locally one-dimensional.

Siegel, R.↗
NASA NTRS1976

ISIS-2 satellite imagery and auroral morphology

Auroral morphology is emphasized over auroral dynamics in a paper describing conspicuous auroral features picked up by the ISIS-2 scanning photometer. Results of improved programs designed to transform the data into a corrected geomagnetic coordinate frame and generate latitude profiles of auroral intensities at different magnetic local times are reported. The diffuse aurora and its relation to the morphology of discrete aurorae is given special attention.

Anger, C. D.↗
NASA NTRS1978

ISEE-C attitude determination using fine sun sensor data only

Techniques developed to determine the spin axis attitude using Fine Sun Sensor (FSS) data only are described. At any given time, the Sun angle specifies the orientation of the spin axis relative to the sunline. The instantaneous time rate of change of the sun angle is directly proportional to the orientation of the spin axis relative to a reference plane that is normal to the ecliptic. Thus, the spin axis attitude can be determined when sufficient data has been collected to accurately measure the rate of change of the sun angle. The uncertainties can be computed directly from the uncertainties in the coefficients of the smoothed sun angle curve. The FSS-only technique is unique in that ephemeris vectors are required only to transform the attitude results to more conventional coordinate frames. The combination of the mission geometry and the FSS accuracy make ISEE-C an ideal mission for applying this method. However, the technique can be used on other missions, such as spin stabilized geosynchronous missions.

Gunshol, L. P.↗
NASA NTRS1982

Hybrid Position/Force Control of Robot Manipulators

In proposed method for task-oriented control of robot manipulator, position and force error signals for each task degree of freedom are used to calculate appropriate control parameters in task coordinates. Position and force error signals are transformed and summed to create drive signal for each actuator. New hybrid control technique does not require operator to supply complex transform matrices. Control trajectories are easily visualized in terms of task to be performed.

Raibert, M. H.↗
NASA NTRS1982

Bubble shapes in steady axisymmetric flows at intermediate Reynolds number

The shape of a gas bubble which rises through a quiescent incompressible, Newtonian fluid at intermediate Reynolds numbers is considered. Exact numerical solutions for the velocity and pressure fields, as well as the bubble shape, are obtained using finite difference techniques and a numerically generated transformation to an orthogonal, boundary-fitted coordinate system. No restriction is placed on the allowable magnitude of deformation.

Ryskin, G.↗
NASA NTRS1983

Applications of the contravariant form of the Navier-Stokes equations

The contravariant Navier-Stokes equations in weak conservation form are well suited to certain fluid flow analysis problems. Three dimensional contravariant momentum equations may be used to obtain Navier-Stokes equations in weak conservation form on a nonplanar two dimensional surface with varying streamsheet thickness. Thus a three dimensional flow can be simulated with two dimensional equations to obtain a quasi-three dimensional solution for viscous flow. When the Navier-Stokes equations on the two dimensional nonplanar surface are transformed to a generalized body fitted mesh coordinate system, the resulting equations are similar to the equations for a body fitted mesh coordinate system on the Euclidean plane. Contravariant momentum components are also useful for analyzing compressible, three dimensional viscous flow through an internal duct by parabolic marching. This type of flow is efficiently analyzed by parabolic marching methods, where the streamwise momentum equation is uncoupled from the two crossflow momentum equations. This can be done, even for ducts with a large amount of turning, if the Navier-Stokes equations are written with contravariant components.

Katsanis, T.↗
NASA NTRS1984

The Io decametric emission cone

The Voyager 1 and 2 Planetary Radio Astronomy observations of Jovian decametric radiation (DAM) were transformed into a fixed Io-Jupiter coordinate system for Io longitudes every 10 deg apart starting from 200 to 260 deg (seven plots per spacecaft). An analysis of the data in this coordinate system has shown that the Io-dependent DAM emissions are generated in hollow cones at all frequencies over the fixed Io longitudes studied. The edges of the emission cones are, in many cases, broad in longitude (30 deg or more), with distinct differences found between the characteristics of DAM emission at the edges of the emission cones. The shapes of the DAM emission cones are explained by using three-dimensional ray-tracing calculations within a model Jovian atmosphere.

Green, J. L.↗