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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Mapdfn.py

SAND2021-4824 O mapdfn.py takes the output of dfnWorks describing a discrete fracture network (DFN) and converts the DFN to an equivalent porous medium representation. This is done by looping through cells in a 3-D model domain, tagging those intersected by fractures, and calculating effective properties. Coordinate transformations are imported from an open source script (transformations.py) developed at the University of California. An associated driver script, mapdfn2pflotran.py, calls the functions in mapdfn.py and writes the mapdfn.py output to hdf5 files for use with PFLOTRAN. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525. SAND2021-4824 O

Stein, Emily↗
DOE OSTI2025

Water Mass Transformation Budgets in Finite‐Volume Generalized Vertical Coordinate Ocean Models

Water Mass Transformation (WMT) theory provides conceptual tools that in principle enable innovative analyses of numerical ocean models; in practice, however, these methods can be challenging to implement and interpret, and therefore remain under-utilized. Our aim is to demonstrate the feasibility of diagnosing all terms in the water mass budget and to exemplify their usefulness for scientific inquiry and model development by quantitatively relating water mass changes, overturning circulations, boundary fluxes, and interior mixing. We begin with a pedagogical derivation of key results of classical WMT theory. We then describe best practices for diagnosing each of the water mass budget terms from the output of Finite-Volume Generalized Vertical Coordinate (FV-GVC) ocean models, including the identification of a non-negligible remainder term as the spurious numerical mixing due to advection scheme discretization errors. We illustrate key aspects of the methodology through the analysis of a polygonal region of the Greater Baltic Sea in a regional demonstration simulation using the Modular Ocean Model v6 (MOM6). We verify the convergence of our WMT diagnostics by brute-force, comparing time-averaged (“offline”) diagnostics on various vertical grids to timestep-averaged (“online”) diagnostics on the native model grid. Finally, we briefly describe a stack of xarray-enabled Python packages for evaluating WMT budgets in FV-GVC models (culminating in the new xwmb package), which is intended to be model-agnostic and available for community use and development.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Diagonalizing the Born–Oppenheimer Hamiltonian via Moyal perturbation theory, nonadiabatic corrections, and translational degrees of freedom

This article describes a method for calculating higher order or nonadiabatic corrections in Born–Oppenheimer theory and its interaction with the translational degrees of freedom. The method uses the Wigner–Weyl correspondence to map nuclear operators into functions on the classical phase space and the Moyal star product to represent operator multiplication on those functions. These are explained in the body of the paper. The result is a power series in κ 2 , where κ = (m/M)1/4 is the usual Born–Oppenheimer parameter. The lowest order term is the usual Born–Oppenheimer approximation, while higher order terms are nonadiabatic corrections. These are needed in calculations of electronic currents, momenta, and densities. The separation of nuclear and electronic degrees of freedom takes place in the context of the exact symmetries (for an isolated molecule) of translations and rotations, and these, especially translations, are explicitly incorporated into our discussion. This article presents an independent derivation of the Moyal expansion in molecular Born–Oppenheimer theory. Here, we show how electronic currents and momenta can be calculated within the framework of Moyal perturbation theory; we derive the transformation laws of the electronic Hamiltonian, the electronic eigenstates, and the derivative couplings under translations; we discuss in detail the rectilinear motion of the molecular center of mass in the Born–Oppenheimer representation; and we show how the elimination of the translational components of the derivative couplings leads to a unitary transformation that has the effect of exactly separating the translational degrees of freedom.

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

Asymptotic vacuum solution at tokamak X-point tip

In the H-mode regime of diverted tokamaks, the presence of strong pressure gradients in the pedestal gives rise to a sizable bootstrap current, together with the Ohmic and Pfirsch–Schlueter currents, close to the separatrix. For such equilibria, the presence of finite current density close to the separatrix requires the reexamination of equilibrium properties. It is almost universally assumed that the two branches of the separatrix (the stable and unstable manifolds) are straight as they cross at the X-point. However, the opposite angles of the plasma-filled segment and vacuum one cannot be equal if the current density does not vanish at the separatrix on the plasma side. We solve this difficulty by chipping off a thin layer of plasma edge so that the sharp corner of the plasma-filled segment becomes a hyperbola. Using the conformal transformation, we found that in the assumption of a hyperbolic boundary, the X point moves beyond the plasma boundary to fall in the vacuum region. An acute angle of the plasma-filled segment leads to an obtuse opposite angle of vacuum segment and vice versa. In the case of an acute angle of the plasma-filled segment, the new X point shifts inside the X point formed by the asymptotes of a hyperbolic boundary; in the case of an obtuse angle of the plasma-filled segment, the new X point shifts outside the X point formed by the asymptotes of a hyperbolic plasma boundary. Furthermore, the results are important for understanding the X point features, which affect the tokamak edge stability and transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Geometry optimization speedup through a geodesic approach to internal coordinates

We present a new geodesic-based method for geometry optimization in a basis set of redundant internal coordinates. Overall, our method updates the molecular geometry by following the geodesic generated by a displacement vector on the internal coordinate manifold, which dramatically reduces the number of steps required to converge to a minimum. Our method can be implemented in any existing optimization code, requiring only implementation of derivatives of the Wilson B-matrix and the ability to numerically solve an ordinary differential equation.

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

Revisiting the pseudo-supercritical path method: An improved formulation for the alchemical calculation of solid–liquid coexistence

Alchemical free energy calculations via molecular dynamics have been applied to obtain thermodynamic properties related to solid–liquid equilibrium conditions, such as melting points. In recent years, the pseudo-supercritical path (PSCP) method has proved to be an important approach to melting point prediction due to its flexibility and applicability. In the present work, we propose improvements to the PSCP alchemical cycle to make it more compact and efficient through a concerted evaluation of different potential energies. The multistate Bennett acceptance ratio (MBAR) estimator was applied at all stages of the new cycle to provide greater accuracy and uniformity, which is essential concerning uncertainty calculations. In particular, for the multistate expansion stage from solid to liquid, we employed the MBAR estimator with a reduced energy function that allows affine transformations of coordinates. Free energy and mean derivative profiles were calculated at different cycle stages for argon, triazole, propenal, and the ionic liquid 1-ethyl-3-methyl-imidazolium hexafluorophosphate. Comparisons showed a better performance of the proposed method than the original PSCP cycle for systems with higher complexity, especially the ionic liquid. Furthermore, a detailed study of the expansion stage revealed that remapping the centers of mass of the molecules or ions is preferable to remapping the coordinates of each atom, yielding better overlap between adjacent states and improving the accuracy of the methodology.

25 ENERGY STORAGE↗
DOE OSTI2025

Modeling cross-beam energy transfer with sector ray tracing

Ray-based cross-beam energy transfer (CBET) models are an essential feature of the radiation-hydrodynamic codes used to simulate inertial confinement fusion implosions, but full 3D ray-based CBET calculations can have a prohibitively high computational cost. Sector ray tracing can be used to reduce the cost by orders of magnitude in cases where the coronal plasma and laser drive can be approximated as spherically symmetric. An extension of sector ray tracing (section ray tracing) can be used to relax the assumption of a spherically symmetric drive while still retaining most of the computational savings of a sector ray trace. We discuss the foundations of sector and section ray tracing and compare them to full ray tracing.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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 OSTI2024

Role of bath-induced many-body interactions in the dissipative phases of the Su-Schrieffer-Heeger model

The Su-Schrieffer-Heeger chain is a prototype example of a symmetry-protected topological insulator. Coupling it nonperturbatively to local thermal environments, either through the intercell or the intracell fermion tunneling elements, modifies the topological window. To understand this effect, we employ the recently developed reaction-coordinate polaron transform (RCPT) method, which allows treating system-bath interactions at arbitrary strengths. The effective system Hamiltonian, which is obtained via the RCPT, exposes the impact of the baths on the SSH chain through renormalization of tunneling elements and the generation of many-body interaction terms. By performing exact diagonalization and computing the ensemble geometric phase, a topological invariant, which is applicable even to systems at finite temperature, we distinguish the trivial band insulator (BI) from the topological insulator (TI) phases. Furthermore, through the RCPT mapping, we are able to pinpoint the main mechanism behind the extension of the parameter space for the TI or the BI phases (depending on the coupling scheme, intracell or intercell), which is the bath-induced, dimerized, many-body interaction. In conclusion, we also study the effect of on-site staggered potentials on the SSH phase diagram and discuss extensions of our method to higher dimensions.

1-dimensional systems↗
DOE OSTI2022

The Importance of a Moving Boundary Approach for Modeling the SEI Layer Growth to Predict Capacity Fade

One of the contributing factors to the aging of lithium-ion batteries is the growth of the solid-electrolyte interphase (SEI) layer. The growth of the SEI layer leads to the irreversible loss of lithium available for cycling and increases the resistance of the battery. Physics-based models in literature model the kinetically limited or solvent diffusion-limited growth. In such models, the interface resistance is a constant, and the contribution to the overpotential of the intercalation reaction from the SEI layer is considered to be ohmic. In this study, we propose a model that describes the growth of the SEI layer on the electrode surface as a moving interface. The transport of lithium ions and the solvent in the electrolyte are affected by this moving interface. The equations that govern the species transport and the potential drop across the SEI layer are derived from dilute solution theory and solved by transforming the coordinates of the moving boundary. The ion transport induces changes in the conductivity across the SEI layer, which affects the potential drop that arises due to its growth. The effects of this potential on capacity fade are studied over cycling the battery.

25 ENERGY STORAGE↗
DOE OSTI2021

Commercial Building Energy Code Field Study: Data Collection Methodology and Protocol

In support of the U.S. Department of Energy’s Commercial Buildings Energy Code Field Study, this data collection methodology and protocol provides guidance on all aspects of undertaking a compliance study, from development of a sampling plan to recruitment to code requirements and compliance checks for each energy code measure specified to be collected. The protocol also includes a data collection form that captures all key information needed for analysis of commercial energy code compliance. This methodology was developed by the Institute for Market Transformation in coordination with Pacific Northwest National Laboratory (PNNL) and the U.S. Department of Energy Building Energy Codes Program with the objective of assisting states, jurisdictions, utilities and others as they seek to measure and demonstrate compliance rates with energy codes in commercial buildings, as well as to target areas for improvement through increased energy code compliance and broader energy-efficiency programs. It is also intended to facilitate a consistent and replicable approach to research studies of this type and establish a transparent data set representing baseline construction practices across the U.S.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

Commercial Building Energy Code Field Study: Data Collection Methodology and Protocol

In support of the U.S. Department of Energy’s Commercial Buildings Energy Code Field Study, this data collection methodology and protocol provides guidance on all aspects of undertaking a compliance study, from development of a sampling plan to recruitment to code requirements and compliance checks for each energy code measure specified to be collected. The protocol also includes a data collection form that captures all key information needed for analysis of commercial energy code compliance. This methodology was developed by the Institute for Market Transformation in coordination with Pacific Northwest National Laboratory (PNNL) and the U.S. Department of Energy Building Energy Codes Program with the objective of assisting states, jurisdictions, utilities and others as they seek to measure and demonstrate compliance rates with energy codes in commercial buildings, as well as to target areas for improvement through increased energy code compliance and broader energy-efficiency programs. It is also intended to facilitate a consistent and replicable approach to research studies of this type and establish a transparent data set representing baseline construction practices across the U.S.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2026

Legacy Devices and Interoperability: Developing a Scientific Strategy for Inclusion and Representation of Older Devices

Legacy devices across grid, transportation, and building systems present complex modernization challenges that extend far beyond simple equipment obsolescence. Pulling from presentations and break-out group discussions from this workshop on legacy devices across all three sectors, this report examines the interconnected factors driving replacement decisions, documents successful modernization initiatives from recent years, and provides strategic guidance for coordinating infrastructure transformation across these critical sectors.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2026

ACDC (Automated Campbell Diagram Code) [SWR-26-042]

This application provides a web-based graphical user interface to generating Campbell Diagrams and visualizing mode shapes for OpenFAST turbine models. Determining the aeroelastic stability and dynamic characteristics of wind turbines is a critical step in turbine design and analysis. Historically, extracting natural frequencies and mode shapes from OpenFAST—the industry-standard whole-turbine simulation code—has been a fragmented and tedious process. It required manual model configuration, command-line linearization execution, and complex post-processing via proprietary scripts to handle rotating-frame dynamics. To address these workflow bottlenecks, we present the Automated Campbell Diagram Code (ACDC), an open-source graphical software tool developed by the National Laboratory of the Rockies (NLR) under the DOE-funded Distributed Wind Aeroelastic Modeling (dWAM) project. ACDC streamlines the end-to-end linearization and stability analysis workflow into a single, intuitive cross-platform application. The software guides users through OpenFAST model configuration, definition of operating points, and the automated execution of steady-state trim and linearization simulations. Under the hood, ACDC automates the complex mathematical post-processing steps required for rotating systems, including Multi-Blade Coordinate (MBC) transformations, eigenanalysis, and advanced modal tracking utilizing the Modal Assurance Criterion (MAC) and spectral clustering. Finally, ACDC processes these results to automatically generate Campbell diagrams and features a robust 3D visualization engine to animate full-system mode shapes. By eliminating the reliance on external post-processing environments and manual data manipulation, ACDC significantly accelerates dynamic analysis and lowers the barrier to entry for wind energy researchers and engineers.

Summerville, Brent [National Laboratory of the Roc↗
DOE OSTI2021

The controlling role of atmosphere in dawsonite versus gibbsite precipitation from tetrahedral aluminate species

In highly alkaline solution, aluminum speciates as the tetrahedrally coordinated aluminate monomer, Al(OH)4- and/or dimer Al2O(OH)62-, yet precipitates as octahedrally coordinated gibbsite (Al(OH)3). This tetrahedral to octahedral transformation governs Al precipitation, which is crucial to worldwide Al production, and to the processing of caustic high-level radioactive wastes. Despite its significance, the transformation pathway remains unknown. Here we explore the roles of atmospheric water and carbon dioxide in mediating the transformation of the tetrahedrally coordinated potassium aluminate dimer salt (K2Al2O(OH)6) to gibbsite versus potassium dawsonite (KAl(CO3)(OH)2). A combination of in-situ attenuated total reflection infrared spectroscopy, ex-situ micro X-ray diffraction, and multivariate curve resolution-alternating least squares chemometrics analysis reveals that humidity plays a key role in the transformation by limiting the amount of alkalinity neutralization by dissolved CO2. Lower humidity favors higher alkalinity and incorporation of carbonate species in the final Al product to form KAl(CO3)(OH)2. Higher humidity enables more acid generation that destabilizes dawsonite and favors gibbsite as the solubility limiting phase. Because the transformation was restricted to occur in thin water films, the results suggest that transition from tetra- to octahedrally coordinated Al does not have to occur in bulk solution, as has often been hypothesized, but may instead appear on the source mineral surface.

Dembowski, Mateusz↗
DOE OSTI2020

Influence of soluble oligomeric aluminum on precipitation in the Al-KOH-H2O system

The role of oligomeric aluminate (Al(OH)4-) species in the precipitation of aluminum phases such as gibbsite (a-Al(OH)3) from aqueous hydroxide solutions remains unclear and difficult to probe directly, despite its importance to developing accurate predictions of Al solubility in highly alkaline systems. Precipitation in this system entails a transition from predominantly tetrahedrally coordinated Al species in solution to octahedrally coordinated Al in gibbsite. Here we report a quantitative study of dissolved Al in the Al-KOH-H2O system using a combination of molecular spectroscopies. We established a relationship between changes in 27Al NMR chemical shifts and the relative intensity of Raman vibrational bands, indicative of variations in the ensemble speciation of Al in solution, and the formation of unique contact ion pair interactions with the aluminate dimer, Al2O(OH)62 . A strong correlation between the extent of Al oligomerization and the amount of solvated Al was demonstrated by systematically varying the KOH:Al molar ratio. The concentration of dissolved oligomeric Al in solution also directly impacted the particle size and morphology of gibbsite; high concentrations of dimeric Al2O(OH)62-, yielded smaller and more numerous anhedral to subhedral gibbsite particles, while low concentrations yielded fewer and larger euhedral gibbsite platelets. The collective observations suggest a key role of the Al2O(OH)62- dimer in promoting gibbsite precipitation from solution, where the potassium ion-paired dimer possibly catalyzes a more rapid transformation of Al coordination from tetrahedral in solution to octahedral in gibbsite.

Dembowski, Mateusz↗
DOE OSTI2021

Insulation Coordination Design for Grid-Connected Solid-State Transformers

The deployment of solid-state transformers (SSTs) in medium-voltage (MV) distribution systems is facing various challenges, especially the challenge of insulation coordination design against grid-originated lightning impulses. In this article, two challenges in existing insulation coordination designs for grid-connected SSTs are identified. One challenge is the mismatch between metal-oxide varistor (MOV) protective levels and SST insulation strength, the other challenge is the incompatibility of standard impulse test on SST protective structures. To address the MOV selection challenge, a novel lightning protection scheme is designed to protect a single-stage SST where the semiconductor modules are directly exposed to external lightning impulses. Here, the in-lab lightning impulse tests are performed to verify the overvoltage attenuation performance of the prototyped lightning protection scheme. To address the impulse test challenge, the surge withstand capability of the protected SST is comprehensively evaluated with a complete set of insulation coordination design procedures beyond the basic insulation level (BIL) test. After these two challenges are addressed, a discussion is presented on substituting conventional transformers with the protected SSTs into insulation-coordinated distribution systems to facilitate the field deployment of SSTs.

solid-state transformer (SST)↗