Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “POLAR COORDINATE”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Axisymmetric hydrodynamics in numerical relativity: treating coordinate singularity, artificial heating and modeling MHD instabilities

Two-dimensional axisymmetric simulations of binary neutron star (BNS) merger remnant are a cheap alternative to 3D simulations. To maintain realism for secular timescales, simulations must avoid accumulated errors from drifts in conserved quantities and artificial heating, and they must model turbulent transport in a way that remains plausible throughout the evolution. It is also crucial to avoid numerical artifacts due to the polar coordinate axis singularity. Methods that behave well near the axis often break flux-conservative form of the hydrodynamic equations, resulting in significant drifts in conserved quantities. We present a flux-conservative scheme that maintains smoothness near the axis without sacrificing conservative formulation of the equations or incurring drifts in conserved global quantities. We compare the numerical performance of different treatments of the hydrodynamic equations when evolving a hypermassive neutron star resembling the remnant of a BNS merger. These simulations demonstrate that the new scheme combines the axis smoothness of non-conservative methods with the mass and angular momentum conservation of other conservative methods on $\sim 10^2$ ms timescales of viscous and neutrino-driven evolution. Because fluid profiles remain smooth in the remnant interior, it is possible to remove artificial heating by evolving the entropy density. We show how physical heating and cooling terms can be easily calculated from source terms of the conservative evolution variables and demonstrate our implementation. Finally, we discuss and implement improvements to the effective viscosity scheme to better model the effect of magnetohydrodynamic instabilities as the remnant evolves.

axisymmetry↗
DOE OSTI2025

Dynamics of McMillan mappings II. axially symmetric map

Here, in this article, we investigate the transverse dynamics of a single particle in a model integrable accelerator lattice, based on a McMillan axially-symmetric electron lens. Although the McMillan e-lens has been considered as a device potentially capable of mitigating collective space charge forces, some of its fundamental properties have not been described yet. The main goal of our work is to close this gap and understand the limitations and potentials of this device. It is worth mentioning that the McMillan axially symmetric map provides the first-order approximations of dynamics for a general linear lattice plus an arbitrary thin lens with motion separable in polar coordinates. Therefore, advancements in its understanding should give us a better picture of more generic and not necessarily integrable round beams. In the first part of the article, we classify all possible regimes with stable trajectories and find the canonical action-angle variables. This provides an evaluation of the dynamical aperture, Poincaré rotation numbers as functions of amplitudes, and thus determines the spread in nonlinear tunes. Also, we provide a parameterization of invariant curves, allowing for the immediate determination of the map image forward and backward in time. The second part investigates the particle dynamics as a function of system parameters. We show that there are three fundamentally different configurations of the accelerator optics causing different regimes of nonlinear oscillations. Each regime is considered in great detail, including the limiting cases of large and small amplitudes. In addition, we analyze the dynamics in Cartesian coordinates and provide a description of observable variables and corresponding spectra.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

PCM1 coordinates centrosome asymmetry with polarized endosome dynamics to regulate daughter cell fate

Vertebrate radial glia progenitors (RGPs) balance self-renewal and differentiation through asymmetric cell division (ACD), which involves unequal centrosome inheritance. How centrosome asymmetry directs cell fate remains poorly understood. Here, we identify Pericentriolar material 1 (Pcm1) as a key player in this process. In zebrafish embryonic RGPs, Pcm1 is asymmetrically associated with Cep83, a mother centrosome marker. Using in vivo time-lapse imaging and nanoscale-resolution expansion microscopy, we detect Pcm1 on Notch ligand-containing endosomes, where it interacts–either directly or indirectly–with Par-3 and dynein. Loss of pcm1 disrupts endosome dynamics, increasing neuronal differentiation at the expense of RGP self-renewal. Mechanistically, Pcm1 facilitates the transition from Rab5b to Rab11a and promotes the assembly of Par-3 and dynein macromolecular complexes on recycling endosomes. Furthermore, we find conserved PARD3-PCM1-CEP83-RAB11 associations in human cortical brain organoids. Our findings uncover that Pcm1 links centrosome asymmetry to polarized endosome trafficking, thereby regulating RGP fate decisions.

Cell fate and cell lineage↗
DOE OSTI2025

Towards a Deeper Fundamental Understanding of (Al,Sc)N Ferroelectric Nitrides

Density functional theory (DFT) calculations, within the virtual crystal alloy approximation, are performed, along with the development of a Landau-type model employing a symmetry-allowed analytical expression of the internal energy and having parameters determined from first principles, to investigate properties and energetics of Al1-xScxN ferroelectric nitrides in their hexagonal forms. These DFT computations and this model predict the existence of two different types of minima, namely, the fourfold-coordinated wurtzite (WZ) polar structure and a five-fold coordinated paraelectric hexagonal phase (denoted as H5), for any Sc composition up to 40%. The H5 minimum progressively becomes the lowest-energy state within hexagonal symmetry as the Sc concentration increases from 0 to 0.4. Furthermore, the model points to several key findings. Examples include the crucial role of the coupling between polarization and strains to create the WZ minimum, in addition to polar and elastic energies, and that the origin of the H5 state overcoming the WZ phase as the global minimum within hexagonal symmetry when increasing the Sc composition mostly lies in the compositional dependency of only two parameters-one linked to the polarization and another one being purely elastic in nature. Other examples are that forcing Al1-xScxN systems to have no or a weak change in lattice parameters when heating them allows us to reproduce their finite-temperature polar properties well and that a value of the axial ratio close to that of the ideal WZ structure implies a large polarization at low temperatures but not necessarily at high temperatures because of the ordered-disordered character of the temperature-induced formation of the WZ state. Such findings should allow for a better fundamental understanding of (Al,Sc)N ferroelectric nitrides, which may be used to design efficient devices having, e.g., low operating voltages.

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

Room temperature crystal field splitting of curium resolved by circularly polarized luminescence spectroscopy

Coordination of Cm(III) with a chiral decadentate ligand N,N,N',N'-tetrakis[(6-carboxypyridin-2-yl)methyl]-1,2-diaminocyclohexane (tpadac) generated complexes with strong luminescence allowing for the unprecedented measurement of well-resolved Cm(III) circularly polarized luminescence spectra. Quantitative resolution of the electronic structure of the [Cm(tpadac)][K] complexes was achieved at room temperature, highlighting the strength of the combination of luminescence and circularly polarized luminescence spectroscopies to unravel the fundamental electronic structure of Cm(III). These results are a clear demonstration that these spectroscopies are powerful yet simple tools for the fundamental understanding of electronic structure, which opens the door to future investigations of other Cm(III) complexes in geometries relevant to nuclear applications, and even other 5f-elements.

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

Trans-Influence in Dinuclear Pt(III) Complexes: Electronic Structure, σ-Donation, and Pt–Pt Spin–Spin Coupling

This study investigates the trans influence in dinuclear platinum(III) complexes using a combined approach of ab initio molecular dynamics and natural localized molecular orbital (NLMO) analysis. Focusing on pivalamidate-bridged Pt III complexes with axial ligands of varying σ-donation strength, it is quantified how ligand−metal interactions propagate through the Pt−Pt bond, and how they affect bond polarization, axial water coordination, and 1 J PtPt spin−spin coupling constants. NLMO analysis reveals quantitatively that strong σ-donating ligands polarize the Pt−Pt bond, shifting the electron density toward the opposite platinum center. The polarization mechanism is identified as the primary reason for the observed reduction of 1 J PtPt , because the bond polarization diminishes the transmission of the nuclear magnetic spin-induced electron spin density through the Pt−Pt bond. Additionally, the destabilization of axial water coordination at the opposite Pt site can be rationalized through a polarizationinduced Pt IV − Pt II -like mixed-valence character.

Ab initio molecular dynamics↗
DOE OSTI2026

Cosolvent-tuned interactions in ionic liquids: A vibrational and quantum-chemical study of ethylene glycol ratio effects

Ionic liquids (ILs) are attractive media for CO 2 capture but remain limited by viscosity and cost. Blending ILs with ethylene glycol (EG) is a practical route to mitigate these constraints, yet the molecular origins of cosolvent effects and their dependence on composition are not well resolved. We combine Fourier-transform infrared (FT-IR) spectroscopy with quantum-chemical (DFT) analysis to elucidate how the IL:EG molar ratio modulates intermolecular interactions and electronic structure. Computed vibrational frequencies enable mode assignment and deconvolution of overlapping bands, revealing systematic, ratio-dependent shifts and broadenings in (i) EG O–H stretching, (ii) cation and EG C–H stretchings (imidazolium C2–H, C4–H, C5–H, methyl and ethyl groups, -CH2 of EG), (iii) anion signature modes (e.g., CN motifs), and (iv) EG C–O and C–C stretchings, consistent with the redistribution of hydrogen-bonding networks. Molecular electrostatic potential (MESP) maps quantify attenuation of extreme potential regions with increasing EG, indicating progressive screening of cation–anion electrostatic interactions. Quantum Theory of Atoms in Molecules (QTAIM) identifies emergent bond critical points between EG and the IL ions, while Reduced Density Gradient–Noncovalent Interaction (RDG–NCI) analysis differentiates strong directional hydrogen bonds from dispersive contacts across compositions. Together, these results show that EG fraction controls a switch from predominantly ion–ion to mixed ion–EG coordination, altering local polarity and polarizability that underlie the observed FT-IR trends. The framework provides composition–structure–spectrum relationships that can guide rational selection of IL:EG ratios to balance favorable molecular interactions with practical performance targets in scalable CO 2 capture systems.

DAC↗
DOE OSTI2025

LDRD Abbreviated report: High-Order General-Discrete-Ordinates Method Enabling Efficient Deterministic Transport in Hydrodynamic Simulations

Deterministic transport simulations for national-security and energy applications often operate in high-dimensional phase-space, where accuracy and cost both become major challenges. A common numerical artifact in such problems is the “ray-effect,” which appears as unphysical streaks. Beyond misinterpretation, these artifacts can contaminate tightly coupled physics, such as fluid dynamics, radiation-hydrodynamics, and laser-plasma interactions, eroding the predictive capability of entire multiphysics workflows. Our objective was to make high-dimension studies practical on modern hardware while mitigating the ray-effect without relying on prohibitively expensive sampling approaches such as Monte Carlo methods. We developed the Generic Discretization Library (GenDiL), a Graphics Processing Unit (GPU)-first framework that uses high-order Discontinuous Galerkin (DG) methods and matrix-free algorithms to reduce memory usage and improve computational efficiency, critical for phase-space simulations. GenDiL supports phase-space adaptivity in both mesh size and polynomial order (hp-adaptivity) to place resolution only where it is needed. A central capability is Local Dimensional Refinement (LDR), which couples lower-dimension continuum models to higher-dimension kinetic models through stable and conservative interfaces, so that high-fidelity physics is applied only in regions where it is essential. Building on the GenDiL framework, we developed the General SN (GSN) family of algorithms as a true generalization of the polar SN approach (discrete ordinates, often denoted SN). Rather than tying discrete ordinates to a specific polar change of coordinates, GSN formulates transport on an arbitrary change of coordinates chosen to reduce ray-effect. We studied two complementary variants: an analytic variant, where the coordinate map is prescribed in advance by a closed-form function; and a data-driven variant, where a quantity of interest, such as the net flux, guides the coordinate system. GenDiL provides the library infrastructure for efficient GPU execution, but the GSN concept is algorithmic and independent of any one library. Across representative high-dimension tests, including non-symmetric solutions, both variants delivered strong ray-effect mitigation at practical cost, moving four- to six-dimensional analysis toward repeatable, routine studies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Reducing uncertainty of polar to midlatitude linkages using DOE’s E3SM in a coordinated model-experiment setting

This project brought DOE’s climate modeling effort with the Energy Exascale Earth System Model into the Polar Amplification Model Intercomparison Project (PAMIP), which is part of the sixth and latest Coupled Model Intercomparison Project, CMIP6. PAMIP examines the causes and consequences of polar amplification, when external forcing results in a larger temperature increase in high latitudes than the global average, in a coordinated set of model experiments with a common modeling protocol. Our teams from UC Irvine and the University of Toronto have designed, carried out, analyzed, and disseminated PAMIP output from the Energy Exascale Earth System Model (E3SM) and the Community Earth System Model (CESM). PAMIP’s ongoing stream of significant new results have advanced progress in the community’s understanding and led to new outstanding research questions that have motivated further work. PAMIP has led to improved consensus on the atmospheric response to sea ice loss. The important finding is that for a similar sea ice anomaly forcing, the simulated atmospheric response in the troposphere is remarkably consistent among the 16 models’ runs analyzed. The zonal-mean tropospheric response consists of a very robust equatorward shift of the westerly flow in mid-latitudes. However, while the multi-model mean response is robust, it has a weak amplitude relative to internal variability. We identified a weakness in the models (including E3SM) in terms of their response to sea-ice forcing that is related to eddy forcing (or nonlinear dynamical effects) at mid- to high latitudes. In fact, E3SM is an outlier in terms of the models participating, and in that sense, it turned out to be a vital participant model. We found that reductions in energy transport due sea-ice loss and involving dry air only are compensated by increases in moist energy transport from warmer sea surface temperature in midlatitudes. This leads to a large spread in energy transport into the Arctic and is a potential source of spread in Arctic amplification. We identified an important role that climate modes, including tropical modes of variability (El Nino and the Southern Oscillation (ENSO); the Quasi-biennial Oscillation (QBO)) play in the response to sea-ice anomalies, including in ocean coupled experiments. Similarly, we identified and quantified the contribution of sea-ice thickness to the atmospheric response compared to the response to sea-ice extent only. We found that it is important to run large ensembles and even with an ensemble size of 100 simulations the response is largely influenced by internal variability. We demonstrated convincingly that Ural blocking, not sea-ice loss, provides the weakening of the stratospheric polar vortex in fall/early winter and a negative phase of the North Atlantic Oscillation that can last for up to two months. However, sea-ice anomalies can influence the background flow so that the response to Ural blocking is more persistent under low sea-ice conditions in the Barents/Kara Sea than high sea ice Atmospheric model hierarchies that progressively add individual processes have a long history in providing dynamical insight for modeling the atmosphere. Similarly, coupled model hierarchies that progressively add individual ocean processes can provide insights into the workings of the coupled climate system, however such hierarchies have not been available except for a non-dynamical slab ocean model. Because of the missing processes, surface flux corrections must be added to produce a target climate. In this project, we managed to overcome this problem and develop a globally coupled ocean model hierarchy in CESM that can turn on and off the processes of mixed-layer entrainment and Ekman flow. We used the hierarchy to study the impact of Arctic sea-ice loss on the climate system. We find that the effect of mixed-layer entrainment on ocean heat uptake influences the atmospheric circulation by shifting the latitudinal positions of the mid-latitude westerly jet and the Intertropical Convergence Zone (ITCZ). In quadrupled CO 2 experiments, we studied how air-sea coupling affects the response of tropical rainfall under global warming. In order to identify the importance of individual ocean processes, we used the hierarchy of ocean models to separate the effects of seasonal mixed-layer entrainment, wind-driven Ekman flows, and frictional flows. We showed that including Ekman and frictional flows allows our simulation to produce the Pacific Ocean's enhanced equatorial warming pattern and equatorward ITCZ contraction noted in previous climate simulations. We also showed that the frictional flow, which has yet to receive much attention, is as important as the Ekman flow in generating equatorial heat convergence.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Stiffness anisotropy coordinates supracellular contractility driving long-range myotube-ECM alignment

The ability of cells to organize into tissues with proper structure and function requires the effective coordination of proliferation, migration, polarization, and differentiation across length scales. Skeletal muscle is innately anisotropic; however, few biomaterials can emulate mechanical anisotropy to determine its influence on tissue patterning without introducing confounding topography. Here, we demonstrate that substrate stiffness anisotropy coordinates contractility-driven collective cellular dynamics resulting in C2C12 myotube alignment over millimeter-scale distances. When cultured on mechanically anisotropic liquid crystalline polymer networks (LCNs) lacking topography, C2C12 myoblasts collectively polarize in the stiffest direction. Cellular coordination is amplified through reciprocal cell-ECM dynamics that emerge during fusion, driving global myotube-ECM ordering. Conversely, myotube alignment was restricted to small local domains with no directional preference on mechanically isotropic LCNs of the same chemical formulation. These findings provide valuable insights for designing biomaterials that mimic anisotropic microenvironments and underscore the importance of stiffness anisotropy in orchestrating tissue morphogenesis.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Constructing nested coordinates inside strongly shaped toroids using an action principle

A new approach for constructing polar-like boundary-conforming coordinates inside a toroid with strongly shaped cross-sections is presented. A coordinate mapping is obtained through a variational approach, which involves identifying extremal points of a proposed action in the mapping space from [0,2π] 2 ×[0,1] to a toroidal domain in $\mathbb{R}$ 3 . This approach employs an action built on the squared Jacobian and radial length. Extensive testing is conducted on general toroidal boundaries using a global Fourier–Zernike basis via action minimisation. The results demonstrate successful coordinate construction capable of accurately describing strongly shaped toroidal domains. The coordinate construction is successfully applied to the computation of three-dimensional magnetohydrodynamic equilibria in the GVEC code where the use of traditional coordinate construction by interpolation from the boundary failed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Multi-Ion Complexes and Competition from Bulky BArF – Anions for Chloride Binding in Ion Pairing Conditions

Cationic charges have long been used to enhance anion binding. Embedding charge introduces strong ion pairing for target anions but also for off-target ions, ultimately generating a mixture of multi-ion species that are hard to identify and quantify. While many sidestep this problem using polar solvents and weakly coordinating ions, these approaches exclude a substantial cross-section of conditions found in applications spanning recognition, assembly, separations, templation, and catalysis. To confront this complexity, we study the binding of an anion to a cationic receptor featuring low shape complementarity in a low-polarity solvent to maximize ionic interactions. We prepared the receptor as a salt of the weakly coordinating tetrakis­[3,5-bis­(trifluoromethyl)­phenyl]­borate (BArF – ) anion and studied the binding of small chloride (Cl – ) and medium-sized iodide (I – ) anions. Surprisingly, the use of the bulky BArF – anion does not suppress ion pairing interactions, with 65% of the receptor being paired at 0.5 mM in dichloromethane. We observe multi-ion receptor-Cl – complexes (2:1, 1:1, 1:2), reinforcing the complexity that emerges when working in low-polarity media. Here, we reveal the dependence of affinity on anion charge density and size and that bulky BArF – counteranions compete for chloride binding. These studies reveal the noninnocence of BArF – anions and strategies to quantify multi-ion species.

Ion pairing↗
DOE OSTI2024

Cosmology Large Angular Scale Surveyor (CLASS): 90 GHz Telescope Pointing, Beam Profile, Window Function, and Polarization Performance

The Cosmology Large Angular Scale Surveyor (CLASS) is a telescope array that observes the cosmic microwave background (CMB) over ∼75% of the sky from the Atacama Desert, Chile, at frequency bands centered near 40, 90, 150, and 220 GHz. CLASS measures the large angular scale CMB polarization to constrain the tensor-to-scalar ratio and the optical depth to last scattering. This paper presents the optical characterization of the 90 GHz telescope. Observations of the Moon establish the pointing while dedicated observations of Jupiter are used for beam calibration. The standard deviations of the pointing error in azimuth, elevation, and boresight angle are 1$^{'}_.$3, 2$^{'}_.$1, and 2$^{'}_.$0, respectively, over the first 3 yr of observations. This corresponds to a pointing uncertainty ∼7% of the beam’s full width at half-maximum (FWHM). The effective azimuthally symmetrized instrument 1D beam estimated at 90 GHz has an FWHM of 0$^°_.$620 ± 0$^°_.$003 and a solid angle of 138.7 ± 0.6(stats.) ± 1.1(sys.) μsr integrated to a radius of 4°. The corresponding beam window function drops to $b^2_ℓ$ = 0.93, 0.71, 0.14 at ℓ = 30, 100, 300, respectively. Far-sidelobes are studied using detector-centered intensity maps of the Moon and measured to be at a level of 10 −3 or below relative to the peak. The polarization angle of Tau A estimated from preliminary survey maps is 149°.6 ± 0°.2(stats.) in equatorial coordinates. The instrumental temperature-to-polarization (T → P) leakage fraction, inferred from per-detector demodulated Jupiter scan data, has a monopole component at the level of 1.7 × 10 −3 , a dipole component with an amplitude of 4.3 × 10 −3 , with no evidence of quadrupolar leakage.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Circularly Polarized Attosecond Pulses Enabled by an Azimuthal Phase and Polarization Grating

High-harmonic generation (HHG) is an extreme nonlinear optical process that can map the properties of an infrared driving laser beam onto short wavelength attosecond pulse trains. However, current techniques for generating circularly polarized high harmonics for probing magnetic materials and chiral systems have limitations: two-color collinear counter-rotating driving lasers result in a low cutoff photon energy, while single-color non-collinear counter-rotating schemes suffer from low conversion efficiency. In this work, we generate circularly polarized attosecond pulse trains by using a structured laser driver which has a rotating polarization and phase grating along the azimuthal coordinate. Furthermore, our experimental and numerical results demonstrate the production of left and right circularly polarized harmonics, which naturally separate upon propagation. Our approach uses a single laser color in a collinear geometry, that can be scaled for high efficiency. Simulations show this scheme can extend into the soft x-ray region when driven by mid-infrared driving lasers, while preserving the same high phase-matching cutoff photon energy as for linearly-polarized high harmonics.

attosecond science↗
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 OSTI2026

Quantifying Outer- and Inner-Coordination Sphere Effects Using Uranium Redox Chemistry in Molten Salt Solutions

Defining the relative influence of intramolecular and intermolecular forces is a fundamental problem in chemistry that is difficult to quantify. To address this challenge, we developed a method to evaluate the relative impact of direct chemical bonding in the inner-coordination sphere vs effects from cations in the outer-coordination sphere by comparative analysis of uranium redox reactivity in various molten salts. We observed that outer-coordination sphere cations (M 1+ ) and inner-coordination sphere anions (X 1– ) both affected uranium redox reactivity, with more polarizing M 1+ and larger X 1– favoring uranium in low oxidation states. Changing M 1+ (Li, Na, K) shifted the U IV + e 1– ⇌ U III (U IV/III ) and U III + 3e 1– → U 0 metal potentials by +330 and +240 mV, respectively. Changing X 1– (Cl, Br, I) caused larger shifts of +440 mV for the U IV/III redox potential and +1060 mV for the U 0 metal deposition potential. Using Coulomb’s Law, we correlated these potentials with electrostatic interactions between UIII and the molten salt. This model provided a facile way of predicting redox chemistry within molten salts.

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

FOILPOLARS (Grassmannian Foil Shape Sweeps for Polar Generation) [SWR-26-095]

FOILPOLARS (Grassmannian Foil Shape Sweeps for Polar Generation): Multifidelity aerodynamic polar data generation for hydrofoil/tidal-turbine airfoil sections. Foilpolars ties together three pieces: *AeroSandbox supplies the baseline airfoil coordinates (UIUC database). *G2Aero parameterizes those shapes on a Grassmannian manifold (Karcher mean + PGA basis) and samples new perturbed shapes around that basis. *XFoil (panel method) and NeuralFoil (neural-network surrogate, shipped with AeroSandbox) each solve the resulting shapes for lift, drag, moment, and pressure at the swept angles of attack, Reynolds numbers, and n_crit values. Design optimization of foil shapes in a computationally efficient way requires polars data across many candidate shapes, not just a handful of baseline foils. However, high-fidelity CFD at that scale is too costly, and naive shape perturbation strays from realistic geometries. FOILPOLARS addresses this by loading baseline airfoils (via AeroSandbox) and mapping them onto a Grassmannian manifold (via G2Aero), computing a Karcher mean and principal geodesic analysis (PGA) basis. New shapes are sampled by perturbing PGA coefficients, keeping them close to the manifold of realistic foils. Each sampled shape is evaluated across a configurable sweep of angle of attack, Reynolds number, and critical amplification factor using two solvers: XFoil (panel method) and NeuralFoil (neural-network surrogate), producing a paired dataset of lift, drag, moment, pressure, convergence, and confidence, indexed alongside each shape's PGA coefficients and shared Grassmannian basis in a single xarray dataset. From this, FOILPOLARS produces convergence summaries and comparison plots per shape, Reynolds number, and n_crit. A command-line interface exposes each pipeline stage independently, supporting data-driven design, optimization, and machine-learning workflows for foils.

Sandhu, Rimple [National Laboratory of the Rockies↗