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 685 records · Page 38

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Optical analogues to the equatorial Kerr–Newman black hole

Optical analogues to black holes allow the investigation of general relativity in a laboratory setting. Previous works have considered analogues to Schwarzschild black holes in an isotropic coordinate system; the major drawback is that required material properties diverge at the horizon. We present the dielectric permittivity and permeability tensors that exactly reproduce the equatorial Kerr–Newman metric, as well as the gradient-index material that reproduces equatorial Kerr–Newman null geodesics. Importantly, the radial profile of the scalar refractive index is finite along all trajectories except at the point of rotation reversal for counter-rotating geodesics. Construction of these analogues is feasible with available ordinary materials. A finite-difference frequency-domain solver of Maxwell’s equations is used to simulate light trajectories around a variety of Kerr–Newman black holes. For reasonably sized experimental systems, ray tracing confirms that null geodesics can be well-approximated in the lab, even when allowing for imperfect construction and experimental error.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

X-point effects on the ideal MHD modes in tokamaks in the description of dual-poloidal-region safety factor

The flux coordinates with dual-region safety factor (q) in the poloidal direction are developed in this work. The X-point effects on the ideal MHD modes in tokamaks are then analyzed using this coordinate system. Since the X-point effects mainly affect the edge region, the modes localized at the tokamak edge are particularly examined. Two types of modes are studied. The first is related to the conventional peeling or peeling-ballooning modes. The mode existence aligned with the local magnetic field in the poloidally core region as observed experimentally is confirmed. The X points are shown to contribute to a stabilizing effect for the conventionally treated modes with the surface-averaged q and with the tokamak edge portion truncated. The other is the axisymmetric modes localized in the vicinity of X points, which can affect the cross-field-line transport near the X points. The existence of axisymmetric modes points to the possibility of applying a toroidally axisymmetric resonant magnetic perturbation (RMP) in the X-point area for mitigating the edge localized modes, which can be an alternative to the current RMP design. The dual q description also has important implications for the existing non-axisymmetric RMP concept. It helps to understand why the RMP suppression of edge localized modes is difficult to achieve in the double-null tokamak configurations and points to the possibility of further improving the current RMP concept by considering the alignment to the local q.

Fourier analysis↗

DOE OSTI2025

A Grad-Shafranov model for compact quasisymmetric stellarators

A Grad-Shafranov equation (GSE) valid for compact quasisymmetric stellarators is derived by an asymptotic expansion around a vacuum field carried to first order. We obtain an equation for the existence of flux surfaces leading up to the GSE. The flux surface label must simultaneously satisfy the existence equation and the GSE, which generally leads to an overdetermined problem. We show how the overdetermined problem can be resolved within our model for a class of hybrid devices similar to that studied by Henneberg and Plunk [Phys. Rev. Res. 6, L022052 (2024)]. We are also able to solve the existence equation for flux surfaces analytically in the most general case by introducing a special coordinate system. This will enable us to carry out an optimization seeking to minimize the error in our GSE while obeying the flux surface existence equation, which will allow us to find solutions outside the class of hybrid devices. This will allow for a coarse-grained approximate search in the space of quasisymmetric equilibria that should be faster than a conventional stellarator optimization. Nevertheless, it would still be necessary to fine-tune the approximate solutions using conventional tools to obtain a more precise optimized equilibrium.

Nikulsin, N. (ORCID:0000000318611777)↗

DOE OSTI2025

Reference Correlations for the Density and Viscosity of Molten Alkali and Alkaline Earth Fluoride Salts

While there is a significant body of literature pertaining to thermophysical property measurements of molten salts, there is often a wide degree of variability among independent measurements of the same compounds. As such, the scientific community benefits greatly from an unbiased, independent assessment of duplicate datasets, so that reference correlations which describe these thermophysical properties as functions of temperature can be determined and then commonly used by researchers, scientists, and engineers. With regard to molten fluoride compounds, a significant time has elapsed since density and viscosity reference correlations have been determined; Janz conducted the most recent effort, in 1988, to provide reference correlations for the densities and viscosities of molten fluoride compounds via the National Standard Reference Data System coordinated by the National Bureau of Standards. Since then, new data have been published for molten fluoride compounds, and a new precedent has surfaced for putting forth reference correlations that involve fitting to multiple primary datasets. In this work, reference correlations are put forth for molten alkali and alkaline earth fluoride compounds in an effort to provide updated, improved correlations for general use. For molten alkali fluoride densities, estimated uncertainties with a 95% confidence interval are summarized as follows: LiF (0.63%), NaF (0.48%), KF (0.76%), RbF (0.93%), and CsF (0.75%). For molten alkaline earth fluoride densities, an estimated uncertainty was not able to be quantified for BeF 2 because of limited data; however, estimated uncertainties with a 95% confidence interval are summarized as follows for the remaining alkaline earth fluorides: MgF 2 (1.5%), CaF 2 (0.92%), SrF 2 (1.6%), and BaF 2 (0.23%). For molten alkali fluoride viscosities, uncertainty was not able to be quantified for RbF and CsF because of limited data; however, estimated uncertainties with a 95% confidence interval are summarized as follows for the remaining alkali fluorides: LiF (4.4%), NaF (3.0%), and KF (4.0%). For molten alkaline earth fluoride viscosities, limited consistent data resulted in the recommendation of single datasets (from literature) that are deemed to be the most trustworthy based on the quality of the underlying experimental studies.

Birri, A. [Oak Ridge National Laboratory (ORNL), O↗

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 OSTI2026

Novel approach to general curvilinear coordinates for plasma fluid applications

In general geometry, plasma fluid equations include nonlinear geometric sources associated with fictitious forces, which pose significant challenges to computer simulations. We reformulate the plasma fluid hierarchy to rigorously preserve geometry and conservation properties critical to numerical simulations, while concealing the geometric sources. In their discrete form, the reformulated models conserve mass, angular momentum, and energy naturally, by simple analogy with the continuum equations. These conservation properties have minimal requirements in discrete space, namely, the anti-symmetry of the first derivative and the orthogonality of the scalar and cross products. By decoupling magnetic geometry, coordinate systems, and numerical discretization, this enables maximum flexibility while preserving physics fidelity. As a testbed, we apply the novel representation to the resistive magnetohydrodynamic system, which involves a complete set of curvilinear operations. We verify the correctness of the approach using steady state liquid metal flows and the classic Orszag–Tang vortex.

Halpern, Federico D. [General Atomics, San Diego, ↗

DOE OSTI2022

Reduction of transport due to magnetic shear in gyrokinetic simulations of the scrape-off layer

We report the effect of varying magnetic shear on scrape-off layer turbulence and profiles is studied via electromagnetic gyrokinetic simulations of a helical scrape-off layer model. We develop a model helical geometry with magnetic shear and a corresponding field-aligned coordinate system, which is used for simulations with the Gkeyll code. We find that perpendicular transport is reduced in cases with stronger shear, resulting in higher peak particle and heat fluxes to the endplates. Electromagnetic effects slightly increase transport in strong shear cases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Disentangling Alzheimer’s disease neurodegeneration from typical brain ageing using machine learning

Abstract Neuroimaging biomarkers that distinguish between changes due to typical brain ageing and Alzheimer’s disease are valuable for determining how much each contributes to cognitive decline. Supervised machine learning models can derive multivariate patterns of brain change related to the two processes, including the Spatial Patterns of Atrophy for Recognition of Alzheimer’s Disease (SPARE-AD) and of Brain Aging (SPARE-BA) scores investigated herein. However, the substantial overlap between brain regions affected in the two processes confounds measuring them independently. We present a methodology, and associated results, towards disentangling the two. T1-weighted MRI scans of 4054 participants (48–95 years) with Alzheimer’s disease, mild cognitive impairment (MCI), or cognitively normal (CN) diagnoses from the Imaging-based coordinate SysTem for AGIng and NeurodeGenerative diseases (iSTAGING) consortium were analysed. Multiple sets of SPARE scores were investigated, in order to probe imaging signatures of certain clinically or molecularly defined sub-cohorts. First, a subset of clinical Alzheimer’s disease patients (n = 718) and age- and sex-matched CN adults (n = 718) were selected based purely on clinical diagnoses to train SPARE-BA1 (regression of age using CN individuals) and SPARE-AD1 (classification of CN versus Alzheimer’s disease) models. Second, analogous groups were selected based on clinical and molecular markers to train SPARE-BA2 and SPARE-AD2 models: amyloid-positive Alzheimer’s disease continuum group (n = 718; consisting of amyloid-positive Alzheimer’s disease, amyloid-positive MCI, amyloid- and tau-positive CN individuals) and amyloid-negative CN group (n = 718). Finally, the combined group of the Alzheimer’s disease continuum and amyloid-negative CN individuals was used to train SPARE-BA3 model, with the intention to estimate brain age regardless of Alzheimer’s disease-related brain changes. The disentangled SPARE models, SPARE-AD2 and SPARE-BA3, derived brain patterns that were more specific to the two types of brain changes. The correlation between the SPARE-BA Gap (SPARE-BA minus chronological age) and SPARE-AD was significantly reduced after the decoupling (r = 0.56–0.06). The correlation of disentangled SPARE-AD was non-inferior to amyloid- and tau-related measurements and to the number of APOE ε4 alleles but was lower to Alzheimer’s disease-related psychometric test scores, suggesting the contribution of advanced brain ageing to the latter. The disentangled SPARE-BA was consistently less correlated with Alzheimer’s disease-related clinical, molecular and genetic variables. By employing conservative molecular diagnoses and introducing Alzheimer’s disease continuum cases to the SPARE-BA model training, we achieved more dissociable neuroanatomical biomarkers of typical brain ageing and Alzheimer’s disease.

Hwang, Gyujoon↗

DOE OSTI2025

Rotational excitation and de-excitation of magnesium mono-sulphide (MgS) by He collisions

ABSTRACT Magnesium mono-sulphide (MgS) plays a crucial role in astrochemical processes, particularly in the interstellar medium (ISM), where metal-sulphur chemistry influences molecular formation and evolution. This work presents a detailed study of the rotational excitation and de-excitation of MgS through collisions with helium (He) atoms, the second most abundant species in the ISM after hydrogen, which includes both atomic (H) and molecular forms (H2). The focus on MgS–He collisions arises from He's high abundance, chemical inertness, and simpler electronic structure, which make it well suited for quantum scattering calculations. These characteristics establish He as an ideal candidate for initial studies, providing fundamental data for future investigations involving H2. The study uses quantum scattering methods to calculate the collisional rate coefficients over a broad temperature range. These rates are critical for interpreting observational data on MgS and predicting its abundance in space. The interaction potential between MgS and He is calculated using the rigid rotor approximation and the Jacobi coordinate system, employing the CCSD(T)-F12a/aug-cc-pVTZ method for accurate two-dimensional potential energy surface. The study explores the anisotropic nature of the MgS–He interaction, which favours odd Δj rotational transitions at low collision energies. The inelastic cross-sections for rotational transitions involving up to 16 rotational levels of MgS were computed up to 1000 cm−1, enabling the calculation of rate coefficients up to 150 K for Δj = ±1, ±2, and ± 3 rotational transitions. The results show that Δj = 1 transitions dominate at low temperatures, while Δj = 2 transitions become more significant at higher temperatures. This study provides valuable data for interpreting future astrophysical observations of MgS. The findings also propose new rotational transitions for MgS detection in space, enhancing our ability to track and study this molecule in various cosmic environments.

Hendaoui, Hamza (ORCID:0000000218641872)↗

DOE OSTI2020

Effective field theory for deformed odd-mass nuclei

We develop an effective field theory (EFT) for deformed odd-mass nuclei. These are described as an axially symmetric core to which a nucleon is coupled. In the coordinate system fixed to the core the nucleon is subject to an axially symmetric potential. Power counting is based on the separation of scales between low-lying rotations and higher-lying states of the core. In leading order, core and nucleon are coupled by universal derivative terms. These comprise a covariant derivative and gauge potentials which account for Coriolis forces and relate to Berry-phase phenomena. At leading order, the EFT combines the particle-rotor and Nilsson models. Here, we work out the EFT up to next-to-leading order and illustrate the results in 239 Pu and 187 Os . At leading order, odd-mass nuclei with rotational bandheads that are close in energy and differ by one unit of angular momentum are triaxially deformed. For bandheads that are well separated in energy, triaxiality becomes a subleading effect. The EFT developed in this paper presents a model-independent approach to the particle-rotor system that is capable of systematic improvement.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

Handing off the outcome of binary neutron star mergers for accurate and long-term postmerger simulations

In this work, we perform binary neutron star (BNS) merger simulations in full dynamical general relativity with IllinoisGRMHD, on a Cartesian grid with adaptive-mesh refinement. After the remnant black hole has become nearly stationary, the evolution of the surrounding accretion disk on Cartesian grids over long timescales (~1s) is suboptimal, as Cartesian coordinates over-resolve the angular coordinates at large distances, and the accreting plasma flows obliquely across coordinate lines dissipating angular momentum artificially from the disk. To address this, we present the HandOff, a set of computational tools that enables the transfer of general relativistic magnetohydrodynamic (GRMHD) and spacetime data from IllinoisGRMHD to harm3d, a GRMHD code that specializes in modeling black hole accretion disks in static spacetimes over long timescales, making use of general coordinate systems with spherical topology. We demonstrate that the handoff allows for a smooth and reliable transition of GRMHD fields and spacetime data, enabling us to efficiently and reliably evolve BNS dynamics well beyond merger. We also discuss future plans, which involve incorporating advanced equations of state and neutrino physics into BNS simulations using the HandOff approach.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

DXRD : a user-friendly suite of two- and multiple-beam dynamical X-ray diffraction programs

The DXRD program suite consisting of a series of dynamical theory programs is introduced for computing dynamical X-ray diffraction from single crystals. Its interactive graphical user interfaces (GUIs) allow general users to make complicated calculations with minimal effort. It can calculate plane-wave Darwin curves of single crystals (or multiple crystals) for both the Bragg and Laue cases, including grazing-incidence diffraction and backward diffraction (with Bragg angles approaching 90°). It is also capable of simulating rocking curves for divergent incident X-ray beams with finite bandwidths. A unique feature of DXRD is that it provides a convenient GUI-based multiple-beam diffraction program that can accurately compute arbitrary N-beam diffraction of any geometry using a universal 4N × 4N matrix method. DXRD also provides a mapping program for plotting all the multiple-beam diffraction lines (monochromator glitches) in the azimuth–energy coordinate system. All these functions make DXRD a convenient and powerful software tool for designing crystal-based synchrotron/X-ray optics (monochromators, analyzers, polarizers, phase plates etc.) and for crystal characterization, X-ray spectroscopy and X-ray diffraction teaching.

Bragg reflection↗

DOE OSTI2026

Framework for X-ray mirror surface shape fitting

For accurate characterization of grazing-incidence X-ray mirrors, we present a comprehensive framework to fit measured surface shapes (either slope or height) of X-ray mirrors used in synchrotron radiation and free-electron laser facilities. We summarize the closed-form expressions of some typical surface shapes of X-ray mirrors including elliptic cylinders, hyperbolic cylinders, ellipsoids, hyperboloids, and diaboloids. This framework is composed of four layers: definition of standard shapes with closed-form expressions, generation of theoretical surface with pose parameters (six degrees of freedom defining an object's position and orientation relative to a coordinate system), parameter optimization with the ability to select which parameters are fit and which are held constant, and the development of user-friendly fitting function wrappers for particular fitting tasks. A few practical fitting examples are demonstrated to verify the effectiveness of the proposed fitting framework. We discuss the physical meanings of the fitting parameters, and provide several examples using the elliptic cylinder and ellipsoid shapes to highlight some features of the framework. Moreover, we provide the presented framework as open-source codes (MATLAB and Python codes available at https://github.com/nsls2omf/xmf) to the community to encourage academic collaboration and further improvements.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Fast and Invertible Simplicial Approximation of Magnetic‐Following Interpolation for Visualizing Fusion Plasma Simulation Data

We introduce a fast and invertible approximation for fusion plasma simulation data represented as 2D planar meshes with connectivities approximating magnetic field lines along the toroidal dimension in deformed 3D toroidal spaces. Scientific variables (e.g., density and temperature) in these fusion data are interpolated following a complex magnetic-field-line-following scheme in the toroidal space represented by a cylindrical coordinate system. This deformation in the 3D space poses challenges for root-finding and interpolation. To this end, we propose a novel paradigm for visualizing and analyzing such data based on a newly developed algorithm for constructing a 3D simplicial mesh within the deformed 3D space. Our algorithm generates a tetrahedral mesh that connects the 2D meshes using tetrahedra while adhering to the constraints on node connectivities imposed by the magnetic field-line scheme. Specifically, we first divide the space into smaller partitions to reduce complexity based on the input geometries and constraints on connectivities. Then, we independently search for a feasible tetrahedralization of each partition, considering nonconvexity. We demonstrate our method with two X-Point Gyrokinetic Code (XGC) simulation datasets on the International Thermonuclear Experimental Reactor (ITER) and Wendelstein 7-X (W7-X), and use an ocean simulation dataset to substantiate broader applicability of our method. An open source implementation of our algorithm is available at https://github.com/rcrcarissa/DeformedSpaceTet.

Ren, Congrong [The Ohio State Univ., Columbus, OH ↗

DOE OSTI2020

Validation of finite-element models using full-field experimental data: Levelling finite-element analysis data through a digital image correlation engine

Full-field data from digital image correlation (DIC) provide rich information for finite-element analysis (FEA) validation. However, there are several inherent inconsistencies between FEA and DIC data that must be rectified before meaningful, quantitative comparisons can be made, including strain formulations, coordinate systems, data locations, strain calculation algorithms, spatial resolutions and data filtering. As such, in this paper, we investigate two full-field validation approaches: (1) the direct interpolation approach, which addresses the first three inconsistencies by interpolating the quantity of interest from one mesh to the other, and (2) the proposed DIC-levelling approach, which addresses all six inconsistencies simultaneously by processing the FEA data through a stereo-DIC simulator to ‘level’ the FEA data to the DIC data in a regularisation sense. Synthetic ‘experimental’ DIC data were generated based on a reference FEA of an exemplar test specimen. The direct interpolation approach was applied, and significant strain errorswere computed, even though therewas no model form error, because the filtering effect of theDIC enginewas neglected. In contrast, the levelling approach provided accurate validation results, with no strain error when no model form error was present. Next, model form error was purposefully introduced via a mismatch of boundary conditions. With the direct interpolation approach, the mismatch in boundary conditions was completely obfuscated, while with the levelling approach, it was clearly observed. Finally, the ‘experimental’ DIC datawere purposefully misaligned slightly fromthe FEA data. Both validation techniques suffered from the misalignment, thus motivating continued efforts to develop a robust alignment process. In conclusion, direct interpolation is insufficient, and the proposed levelling approach is required to ensure that the FEA and the DIC data have the same spatial resolution and data filtering. Only after the FEA data have been ‘levelled’ to the DIC data can meaningful, quantitative error maps be computed.

42 ENGINEERING↗

DOE OSTI2024

Development of Modified Perturbation Solutions to the One-Phase Stefan Problems With a Convective Boundary

The classical Stefan problem is used to track the moving solid-liquid interface during the freezing process. Perturbation theory has often been applied to find an approximate analytical solution due to the nonlinearity of the moving interface. However, the Stefan number (i.e., the sensible over latent heat) must be small and usually less than 0.01 to assume the perturbation expansion, which in turn limits the thermal engineering applications. In this study, a modified perturbation solution is developed by adding a correction term after the leading-order solution to be valid for a much wider range of Stefan numbers (i.e., 0.01 less than or equal to Ste less than or equal to 1). Specifically, a one-phase Stefan problem is first formulated subjected to a convective boundary in the Cartesian, cylindrical, and spherical coordinate systems for inward solidification. The leading-order solution is calculated based on the regular perturbation theory, while the correction term is obtained using the Monte-Carlo method and a multi-variant regression. Results show that the correction term has a linear relationship with the Stefan number and is not significantly influenced by the Biot number. The proposed modified perturbation solution can accurately and rapidly predict the nonlinear moving interface motion for the freezing process.

DIRECT ENERGY CONVERSION,MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Analytical expressions of the surface shape of ‘diaboloid’ mirrors

Modern deterministic polishing processes allow fabrication of x-ray optics with almost any arbitrary aspherical surface shape. Among these optics, the so called "diaboloid"mirror is of special interest. The diaboloid mirror that converts a cylindrical wave to a spherical wave would improve focusing in x-ray beamlines implementing a diffraction element between a parabolic cylinder and a toroidal mirror. The replacement of the toroidal mirror in existing beamlines by the diaboloid mirror would mitigate aberrations. The shape of the diaboloid mirror is usually calculated numerically based on a truncated polynomial solution of the optical path problem. Here, we present an exact analytical solution for the shape of a diaboloid mirror as a function of the conjugate parameters of the mirror placed in a beamline. The derived analytical expressions for the diaboloid mirror in both the canonical and mirror-based coordinate systems are implemented in ray-tracing simulations to verify the beamline performances.

Yashchuk, Valeriy V.↗

DOE OSTI2024

Observation and differential cross section measurement of neutral current DIS events with an empty hemisphere in the Breit frame

The Breit frame provides a natural frame to analyze lepton–proton scattering events. In this reference frame, the parton model hard interactions between a quark and an exchanged boson defines the coordinate system such that the struck quark is back-scattered along the virtual photon momentum direction. In Quantum Chromodynamics (QCD), higher order perturbative or non-perturbative effects can change this picture drastically. As Bjorken-x decreases below one half, a rather peculiar event signature is predicted with increasing probability, where no radiation is present in one of the two Breit-frame hemispheres and all emissions are to be found in the other hemisphere. At higher orders in $\alpha$ s or in the presence of soft QCD effects, predictions of the rate of these events are far from trivial, and that motivates measurements with real data. We report on the first observation of the empty current hemisphere events in electron–proton collisions at the HERA collider using data recorded with the H1 detector at a center-of-mass energy of 319 GeV. The fraction of inclusive neutral-current DIS events with an empty hemisphere is found to be 0.0112 ± 3.9% stat ± 4.5% syst ± 1.6% mod in the selected kinematic region of 150 < Q 2 < 1500 GeV 2 and inelasticity 0.14 < y < 0.7. The data sample corresponds to an integrated luminosity of 351.1 pb -1 , sufficient to enable differential cross section measurements of these events. The results show an enhanced discriminating power at lower Bjorken-x among different Monte Carlo event generator predictions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗