SEARCH · Engineering Papers

Results for “Models based on symmetries”

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 577 records · Page 32

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Machine Learning-Driven Conservative-to-Primitive Conversion in Hybrid Piecewise Polytropic and Tabulated Equations of State

We present a novel machine learning (ML)-based method to accelerate conservative-to-primitive inversion, focusing on hybrid piecewise polytropic and tabulated equations of state. Traditional root-finding techniques are computationally expensive, particularly for large-scale relativistic hydrodynamics simulations. To address this, we employ feedforward neural networks (NNC2PS and NNC2PL), trained in PyTorch (2.0+) and optimized for GPU inference using NVIDIA TensorRT (8.4.1), achieving significant speedups with minimal accuracy loss. The NNC2PS model achieves 𝐿 1 and 𝐿 ∞ errors of 4.54 × 10 −7 and 3.44 × 10−6, respectively, while the NNC2PL model exhibits even lower error values. TensorRT optimization with mixed-precision deployment substantially accelerates performance compared to traditional root-finding methods. Specifically, the mixed-precision TensorRT engine for NNC2PS achieves inference speeds approximately 400 times faster than a traditional single-threaded CPU implementation for a dataset size of 1,000,000 points. Ideal parallelization across an entire compute node in the Delta supercomputer (dual AMD 64-core 2.45 GHz Milan processors and 8 NVIDIA A100 GPUs with 40 GB HBM2 RAM and NVLink) predicts a 25-fold speedup for TensorRT over an optimally parallelized numerical method when processing 8 million data points. Moreover, the ML method exhibits sub-linear scaling with increasing dataset sizes. We release the scientific software developed, enabling further validation and extension of our findings. By exploiting the underlying symmetries within the equation of state, these findings highlight the potential of ML, combined with GPU optimization and model quantization, to accelerate conservative-to-primitive inversion in relativistic hydrodynamics simulations.

conservative-to-primitive conversion↗

DOE OSTI2026

Testing Classical Properties from Quantum Data

Many properties of Boolean functions can be tested far more efficiently than the function itself can be learned. However, this dramatic advantage often disappears when testers are limited to random samples of ƒ instead of adaptively chosen queries to f. In this work we investigate the quantum version of this restriction: quantum algorithms that test properties of a Boolean function f solely from copies of either the function state |ƒ⟩ ∝ ∑ x |x, ƒ(x)⟩ or the phase state |(-1) ƒ ⟩ ∝ ∑ x (-1) ƒ(x) |x⟩. For monotonicity, symmetry, and triangle-freeness, we show passive quantum testers are unboundedly or super-polynomially better than their classical passive testing counterparts. They are competitive with classic query -based testers in each case. Our new testers use techniques beyond quantum Fourier sampling, and it turns out this is necessary: we show a certain class of bent functions can be tested from 𝒪(1) function states but has a sample complexity lower bound of 2 Ω(n) for any tester relying exclusively on Fourier and classical samples. Our passive quantum testers are competitive with classical query -based testers, but this isn't universal: we exhibit a testing problem that can be solved from 𝒪(1) classical queries but requires Ω(2 n/2 ) function state copies. The Forrelation problem provides a separation of the same magnitude in the opposite direction, so we conclude that quantum data and classical queries are "maximally incomparable" resources for testing. We also begin the study of lower bounds for testing from quantum data. For quantum monotonicity testing, we prove that the ensembles of [Goldreich et al., 2000; Black, 2024], which give exponential lower bounds for classical sample-based testing, do not yield any nontrivial lower bounds for testing from quantum data. New insights specific to quantum data will be required for proving copy complexity lower bounds for testing in this model.

Boolean Functions↗

DOE OSTI2025

Flux-Closure Domain Structures in Ferroelectric K 0.5 Na 0.5 NbO 3 Thin Films

Topological domain structures in ferroelectric materials have garnered increasing attention due to their intriguing physical properties and promising applications. While most existing topological structures in ferroelectric perovskite oxides originate from tetragonal or rhombohedral bulk phases, much less is understood about their counterparts in orthorhombic ferroelectrics. Here, in this work, we employ ferroelectric K 0.5 Na 0.5 NbO 3 (KNN) thin films as a model system and leverage phase-field simulations to theoretically predict the static structures and dynamic behaviors of three types of flux-closure domain configurations: in-plane (Type-I), out-of-plane (Type-II), and superdomain (Type-III) flux-closure structures. We systematically investigate the effects of finite size, misfit strains, and electrical boundary conditions on the formation and switching of these topological structures. For the Type-I structure, size reduction or small misfit strain facilitates a transition of the flux-closure pattern to polar vortices. Type-II structures emerge under open-circuit electrical boundary conditions of the film, forming at the junctions of specific domain walls with the film surface or the film–substrate interface. The formation mechanisms of these two flux-closure structures are rationalized from an energy perspective. We demonstrated switching capabilities of Type-II and Type-III structures by obtaining polarization–electric field hysteresis loops. Our simulations also reveal a reversible electric-field-induced transition between the orthorhombic and rhombohedral ferroelectric phases with a checkerboard domain pattern, during which the integrity of the flux-closure structure is preserved. These findings provide theoretical insights and practical guidance for identifying and manipulating topological structures in low-symmetry ferroelectrics, paving the way for developing energy-efficient microelectronic devices based on topological structures.

P-E loop↗

DOE OSTI2024

Search for Baryon Number Violation in Top Quark Production and Decay Using Proton-Proton Collisions at $\sqrt{𝑠}$ = 13 TeV

A search is presented for baryon number violating interactions in top quark production and decay. The analysis uses data from proton-proton collisions at a center-of-mass energy of 13 TeV, collected with the CMS detector at the LHC with an integrated luminosity of 138 fb −1 . Candidate events are selected by requiring two oppositely charged leptons (electrons or muons) and exactly one jet identified as originating from a bottom quark. Multivariate discriminants are used to separate the signal from the background. No significant deviation from the standard model prediction is observed. Upper limits are placed on the strength of baryon number violating couplings. For the first time the production of single top quarks via baryon number violating interactions is studied. This allows the search to set the most stringent constraints to date on the branching fraction of the top quark decay to a lepton, an up-type quark (𝑢 or 𝑐), and a down-type quark (𝑑, 𝑠, or 𝑏). The results improve the previous bounds by 3 to 6 orders of magnitude based on the fermion flavor combination of the baryon number violating interactions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

Symmetry breaking in an extended O(2) model

Motivated by attempts to quantum simulate lattice models with continuous Abelian symmetries using discrete approximations, we study an extended-O(2) model in two dimensions that differs from the ordinary O(2) model by the addition of an explicit symmetry breaking term − h q cos ( q φ ) . Its coupling h q allows to smoothly interpolate between the O(2) model ( h q = 0 ) and a q -state clock model ( h q → ∞ ). In the latter case, a q -state clock model can also be defined for noninteger values of q . Thus, such a limit can also be considered as an analytic continuation of an ordinary q -state clock model to noninteger q . In previous work, we established the phase diagram for noninteger q in the infinite coupling limit ( h q → ∞ ). We showed that there is a second-order phase transition at low temperature and a crossover at high temperature. In this work, we seek to establish the phase diagram at finite values of the coupling using Monte Carlo and tensor methods. We show that for noninteger q , the second-order phase transition at low temperature and crossover at high temperature persist to finite coupling. For integer q = 2 , 3, 4, we know there is a second-order phase transition at infinite coupling (i.e. the well-known clock models). At finite coupling, we find that the critical exponents for q = 3 , 4 vary with the coupling, and for q = 4 the transition may turn into a Berezinskii-Kosterlitz-Thouless transition at small coupling. We comment on the similarities and differences of the phase diagrams with those of quantum simulators of the Abelian-Higgs model based on ladder-shaped arrays of Rydberg atoms. Published by the American Physical Society 2024

Astronomy & Astrophysics↗

DOE OSTI2024

Capturing CO 2 in Quadrupolar Binding Pockets: Broadband Microwave Spectroscopy of Pyrimidine-(CO 2 ) n , n = 1,2

Pyrimidine has two in-plane CH(δ+)/N̈(δ–)/CH(δ+) binding sites that are complementary to the (δ–/2δ+/δ–) quadrupole moment of CO 2 . For this study, we recorded broadband microwave spectra over the 7.5–17.5 GHz range for pyrimidine-(CO 2 ) n with n = 1 and 2 formed in a supersonic expansion. Based on fits of the rotational transitions, including nuclear hyperfine splitting due to the two 14 N nuclei, we have assigned 313 hyperfine components across 105 rotational transitions for the n = 1 complex and 208 hyperfine components across 105 rotational transitions for the n = 2 complex. The pyrimidine-CO 2 complex is planar, with CO 2 occupying one of the quadrupolar binding sites, forming a structure in which the CO 2 is stabilized in the plane by interactions with the C–H hydrogens adjacent to the nitrogen atom. This structure is closely analogous to that of the pyridine-CO 2 complex studied previously by (Doran, J. L. J. Mol. Struct. 2012, 1019, 191–195). The fit to the n = 2 cluster gives rotational constants consistent with a planar cluster of C 2v symmetry in which the second CO 2 molecule binds in the second quadrupolar binding pocket on the opposite side of the ring. The calculated total binding energy in pyrimidine-CO 2 is –13.7 kJ mol –1 , including corrections for basis set superposition error and zero-point energy, at the CCSD(T)/ 6-311++G(3df,2p) level, while that in pyrimidine-(CO 2 ) 2 is almost exactly double that size, indicating little interaction between the two CO 2 molecules in the two binding sites. The enthalpy, entropy, and free energy of binding are also calculated at 300 K within the harmonic oscillator/rigid-rotor model. This model is shown to lack quantitative accuracy when it is applied to the formation of weakly bound complexes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS2002

Design, Fabrication, Optical Testing, and Performance of Diamond Machined Aspheric Mirrors for Ground-Based Near-IR Astronomy

Challenges in fabrication and testing have historically limited the choice of surfaces available for the design of reflective optical instruments. Spherical and conic mirrors are common, but, for future science instruments, more degrees of freedom are necessary to meet challenging performance and packaging requirements. These instruments will be composed of unusual aspheres located far off-axis with large spherical departure, and some designs will require asymmetric surface profiles. In particular, single-surface astigmatism correction in spectrographs necessitates a toroidal surface, which lacks an axis of rotational symmetry. We describe the design, fabrication, optical testing, and performance of three rotationally symmetric, off-axis, aspheric mirrors and one toroidal, off-axis, biconic camera mirror on aluminum substrates for the Infrared Multi-Object Spectrograph (IRMOS) instrument. IRMOS is a facility instrument for the Kitt Peak National Observatory's Mayall Telescope (3.8 m) and an engineering prototype for a possible design of the Next Generation Space Telescope/Multi-Object Spectrograph. The symmetric mirrors range in aperture from 94x86 mm to 286x269 mm and in f-number from 0.9 to 2.4. They are various off-axis, convex and concave, prolate and oblate ellipsoids. The concave biconic mirror has a 94x76 mm aperture, Rx=377 mm, kx=0.0778, Ry=407 mm, and ky=0.1265 and is decentered. by -2 mm in x and 227 mm in y. The mirrors have an aspect ratio of approximately 4:1. The surface error fabrication tolerances are less than 63.3 nm RMS figure error and less than 10 nm RMS microroughness. The mirrors are attached to the instrument bench via a semi-kinematic, integral flexure mount. We describe mirror design, diamond machining, the results of figure testing using computer-generated holograms, and imaging and scattered light modeling and performance.

Ohl, Raymond G.↗

DOE OSTI2024

Lie–Poisson Neural Networks (LPNets): Data-based computing of Hamiltonian systems with symmetries

An accurate data-based prediction of the long-term evolution of Hamiltonian systems requires a network that preserves the appropriate structure under each time step. Every Hamiltonian system contains two essential ingredients: the Poisson bracket and the Hamiltonian. Hamiltonian systems with symmetries, whose paradigm examples are the Lie–Poisson systems, have been shown to describe a broad category of physical phenomena, from satellite motion to underwater vehicles, fluids, geophysical applications, complex fluids, and plasma physics. The Poisson bracket in these systems comes from the symmetries, while the Hamiltonian comes from the underlying physics. We view the symmetry of the system as primary, hence the Lie–Poisson bracket is known exactly, whereas the Hamiltonian is regarded as coming from physics and is considered not known, or known approximately. Using this approach, we develop a network based on transformations that exactly preserve the Poisson bracket and the special functions of the Lie–Poisson systems (Casimirs) to machine precision. We present two flavors of such systems: one, where the parameters of transformations are computed from data using a dense neural network (LPNets), and another, where the composition of transformations is used as building blocks (G-LPNets). We also show how to adapt these methods to a larger class of Poisson brackets. We apply the resulting methods to several examples, such as rigid body (satellite) motion, underwater vehicles, a particle in a magnetic field, and others. The methods developed in this paper are important for the construction of accurate data-based methods for simulating the long-term dynamics of physical systems.

Data-based modeling↗

DOE OSTI2023

New Theories for Neutrino Masses and Dark Matter

The main goal of the project is to investigate new theories for neutrino masses and dark matter, and understand their testability in different experiments. The PI proposed several mechanisms for neutrino masses in the context of B-L theories, left-right symmetric theories, supersymmetric theories and grand unified theories. The PI would like to investigate new ways to understand the origin of neutrino masses where the seesaw scale is in the multi-TeV region. One can have a multi-TeV seesaw scale in scenarios where the same symmetry relevant for neutrino masses defines the dark matter relic density, or in the minimal supersymmetric theory based on local B-L. In these theories one can have new ways to understand the testability of the origin of neutrino masses at colliders and low energy experiments. In these theories one predicts the possibility to observe lepton number violating signatures at the LHC, one can predict large contributions to lepton number violating processes such as mu to e conversion, neutrinoless double beta decay experiments and others. The nature of the dark matter in the Universe is one of the most important problems in cosmology. Two of the most popular candidates are the Axions and the Weakly Interacting Massive Particles (WIMPs). The PI would like to investigate the implications of a new electroweak theta term similar to the QCD vacuum angle but in the SU(2) gauge sector of the SM which is physical if the baryon and lepton numbers are broken symmetries. We will investigate the implications of having a new electroweak theta term dark matter portal, the implications for baryogenesis and understand the constraints coming from different experiments such as the searches for electric dipole moments and axion experiments. The implications of having the electroweak theta term in theories with different sources of baryon and lepton number violation will be investigated. WIMPs are perhaps the most appealing candidates for many reasons. The PI would like to investigate the properties of WIMPs in different gauge theories, study the implications of the cosmological bounds on the WIMPs relic density for the symmetry breaking scale in models for new physics. The PI would like to investigate in great detail the predictions for gamma lines, neutrino lines and other striking signatures which can be present in these theories. These studies can have a profound impact in the testability of different theories for physics beyond the Standard Model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2020

The fermionic universal one-loop effective action

Recent development of path integral matching techniques based on the covariant derivative expansion has made manifest a universal structure of one-loop effective Lagrangians. The universal terms can be computed once and for all to serve as a reference for one-loop matching calculations and to ease their automation. Here we present the fermionic universal one-loop effective action (UOLEA), resulting from integrating out heavy fermions (Dirac or Majorana) with scalar, pseudo-scalar, vector and axial-vector couplings. We also clarify the relation of the new terms computed here to terms previously computed in the literature and those that remain to complete the UOLEA. Our results can be readily used to efficiently obtain analytical expressions for effective operators arising from heavy fermion loops.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2021

SPiRaL: a multiresolution global tomography model of seismic wave speeds and radial anisotropy variations in the crust and mantle

SUMMARY SPiRaL is a joint global-scale model of wave speeds (P and S) and anisotropy (vertical transverse isotropy, VTI) variations in the crust and mantle. The model is comprised of >2.1 million nodes with five parameters at each node that capture velocity variations for P- and S-waves travelling at arbitrary directions in transversely isotropic media with a vertical symmetry axis (VTI). The crust (including ice, water, sediments and crystalline layers) is directly incorporated into the model. The default node spacing is approximately 2° in the lower mantle and 1° in the crust and upper mantle. The grid is refined with ∼0.25° minimum node spacing in highly sampled regions of the crust and upper mantle throughout North America and Eurasia. The data considered in the construction of SPiRaL includes millions of body wave traveltimes (crustal, regional and teleseismic phases with multiples) and surface wave (Rayleigh and Love) dispersion. A multiresolution inversion approach is employed to capture long-wavelength heterogeneities commonly depicted in global-scale tomography images as well as more localized details that are typically resolved in more focused regional-scale studies. Our previous work has demonstrated that such global-scale models with regional-scale detail can accurately predict both teleseismic and regional body wave traveltimes, which is necessary for more accurate location of small seismic events that may have limited signal at teleseismic distances. SPiRaL was constructed to predict traveltimes for event location and long-period waveform dispersion for seismic source inversion applications in regions without sufficiently tuned models. SPiRaL may also serve as a starting model for full-waveform inversion (FWI) with the goal of fitting waves with periods 10–50 s over multiple broad regions (thousands of kilometres) and potentially the globe. To gain insight to this possibility, we simulated waveforms for a small set of events using SPiRaL and independent waveform-based models for comparison. For the events tested, the performance of the traveltime-based SPiRaL model is shown to be generally on par with regional 3-D waveform-based models in three regions (western United States, Middle East, Korean Peninsula) suggesting SPiRaL may serve as a starting model for FWI over broad regions.

58 GEOSCIENCES↗

DOE OSTI2020

One Structure, Two Elements—LuGe 2 Superconductor vs Ordinary Metallic Conductor LuSn 2 . A Case Study on How Site-Selective Germanium for Tin Atom Substitution Leads to Modulating of the Charge Distribution

The substitution of chemically similar elements in a given crystal structure is an effective way to enhance physical properties, but the understanding on such improvements is usually impeded because the substitutions are random, and the roles of the different atoms cannot be distinguished by crystallographic symmetry. Herein, we provide a detailed crystallographic analysis and property measurements for the continuous solid solutions LuGe x Sn 2–x (0 < x < 2). The results show that there is no apparent change of the global symmetry, with the end-members LuGe 2 and LuSn 2 , as well as the intermediate LuGe x Sn 2–x compositions adopting the ZrSi 2 type structure (space group Cmcm, Pearson index oC12). Yet, the refinements of the crystal structures from single-crystal X-ray diffraction data show that Ge–Sn atom substitutions are not random, but occur preferentially at the zig-zag chain. The patterned distribution of two group 14 elements leads to a significant variation in chemical bonding and charge ordering within the other structural fragment, the 2D square nets, thereby resulting in tuned electron transport. The enhancement is greater than the typical Bloch-Gruneisen model, and more akin to the parallel-resistor model. Magnetization measurements on single crystals show bulk superconductivity in all LuGe x Sn 2–x samples with shielding fractions as high as 90%. Specific heat data confirm the effect to originate from residual metallic tin in the material, indicating that Sn atom substitutions in the 2D square nets cause disruptions of the hypervalent bonding and local anisotropy, which ultimately leads to vanishing of the superconducting state in the end-member LuGe 2 . This work sheds light on how the complexity in chemical interactions by two different carbon congeners leads to changes in the physical properties and how they can be correlated with the induced charge distribution. Furthermore, these studies also provide a general approach to modulation of charge density, and thus, of emerging physical properties in other classes of intermetallic systems based on the main-group elements of groups 13 to 15.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS2017

Partitioning of F and Cl Between Apatite and a Synthetic Shergottite Liquid (QUE 94201) at 4 Gpa from 1300 TO 1500 C

Apatite [Ca5(PO4)3(F,Cl,OH)] is present in a wide range of planetary materials. Due to the presence of volatiles within its crystal structure (Xsite), many recent studies have attempted to use apatite to constrain the volatile contents of planetary magmas and mantle sources. In order to use the volatile contents of apatite to accurately determine the abundances of volatiles in coexisting silicate melt or fluids, thermodynamic models for the apatite solid solution and for the apatite components in multicomponent silicate melts and fluids are required. Although some thermodynamic models for apatite have been developed, they are incomplete. Furthermore, no mixing model is available for all of the apatite components in silicate melts or fluids, especially for F and Cl components. Several experimental studies have investigated the apatite-melt and apatite-fluid partitioning behavior of F, Cl, and OH in terrestrial and planetary systems, which have determined that apatite-melt partitioning of volatiles are best described as exchange equilibria similar to Fe-Mg partitioning between olivine and silicate melt. However, McCubbin et al., recently reported that the exchange coefficients vary in portions of apatite compositional space where F, Cl, and OH do not mix ideally in apatite. In particular, solution calorimetry data of apatite compositions along the F-Cl join exhibit substantial excess enthalpies of mixing, and McCubbin et al. reported substantial deviations in the Cl-F exchange Kd along the F-Cl apatite join that could be explained by the preferential incorporation of F into apatite. In the present study, we assess the effect of apatite crystal chemistry on F-Cl exchange equilibria between apatite and melt at 4 GPa over the temperature range of 1300-1500 C. The goal of these experiments is to assess the variation in the Ap-melt Cl-F exchange Kd over a broad range of F:Cl ratios in apatite. The results of these experiments could be used to understand at what composition apatite shifts from a hexagonal unit cell with space group P63/m to a unit cell with monoclinic symmetry within space group P21/b. We anticipate that this transition occurs at >70% chlorapatite based on solution calorimetry data.

McCubbin, F. M.↗

DOE OSTI2020

Residual lattice strain in quartzites as a potential palaeo-piezometer

SUMMARY If a crystal lattice is subjected to a stress, it becomes distorted and no longer represents the ideal crystal symmetry, and if the stress introduces defects such as dislocations, some of this distortion is preserved after the applied stress is removed. In this study, we investigate lattice distortion in quartz at the micron scale with synchrotron X-ray Laue diffraction. From Laue images the local deviatoric strain tensor is derived and corresponding stresses are calculated based on elastic properties. The method is applied to metasedimentary quartzites from the Bergell Alps that were deformed at conditions of greenschist facies metamorphism. The residual palaeostrain is represented in maps of the deviatoric strain tensor components and with deviatoric strain axis pole figures. Data suggest overall shortening perpendicular to the schistosity plane but with considerable asymmetry relative to foliation and lineation, probably attributed to simple shear. Crystallographic pole figures from Laue diffraction agree with neutron diffraction and EBSD measurements and display quartz c-axes girdle distributions with maxima also perpendicular to schistosity. The method shows promise to be used as a palaeo-piezometer to unravel the stress field during tectonic deformation.

58 GEOSCIENCES↗

DOE OSTI2025

Composite quarks and leptons with embedded QCD

We construct a model of quark and lepton compositeness based on an 𝑆⁢𝑈⁡(15) gauge interaction that confines chiral preons, which are also charged under the weakly coupled 𝑆⁢𝑈⁢(4) PS ×𝑆⁢𝑈⁢(2) 𝐿 ×𝑆⁢𝑈⁢(2) 𝑅 gauge group. The breaking of the latter, down to the Standard Model group, is achieved by scalar 𝑆⁢𝑈⁡(15) bound states at a scale in the 30–100 TeV range. The embedding of the QCD gauge group in 𝑆⁢𝑈⁢(4) PS slows down the running of 𝛼 𝑠 in the UV. We estimate the effects of the strongly coupled 𝑆⁢𝑈⁡(15) dynamics on the running of the 𝑆⁢𝑈⁢(4) PS ×𝑆⁢𝑈⁢(2) 𝐿 ×𝑆⁢𝑈⁢(2) 𝑅 gauge couplings, which likely remain perturbative beyond the compositeness scale of about 10 3 –10 4 TeV, and even above a unification scale. A composite vectorlike lepton doublet acquires a mass in the TeV range probed at future colliders, and an extended Higgs sector arises from six-preon bound states.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS2006

Numerical Modeling of Flow Control in a Boundary-Layer-Ingesting Offset Inlet Diffuser at Transonic Mach Numbers

This paper will investigate the validation of the NASA developed, Reynolds-averaged Navier-Stokes (RANS) flow solver, OVERFLOW, for a boundary-layer-ingesting (BLI) offset (S-shaped) inlet in transonic flow with passive and active flow control devices as well as a baseline case. Numerical simulations are compared to wind tunnel results of a BLI inlet experiment conducted at the NASA Langley 0.3-Meter Transonic Cryogenic Tunnel. Comparisons of inlet flow distortion, pressure recovery, and inlet wall pressures are performed. The numerical simulations are compared to the BLI inlet data at a free-stream Mach number of 0.85 and a Reynolds number of approximately 2 million based on the fanface diameter. The numerical simulations with and without tunnel walls are performed, quantifying tunnel wall effects on the BLI inlet flow. A comparison is made between the numerical simulations and the BLI inlet experiment for the baseline and VG vane cases at various inlet mass flow rates. A comparison is also made to a BLI inlet jet configuration for varying actuator mass flow rates at a fixed inlet mass flow rate. Overall, the numerical simulations were able to predict the baseline circumferential flow distortion, DPCP avg, very well within the designed operating range of the BLI inlet. A comparison of the average total pressure recovery showed that the simulations were able to predict the trends but had a negative 0.01 offset when compared to the experimental levels. Numerical simulations of the baseline inlet flow also showed good agreement with the experimental inlet centerline surface pressures. The vane case showed that the CFD predicted the correct trends in the circumferential distortion levels for varying inlet mass flow but had a distortion level that was nearly twice as large as the experiment. Comparison to circumferential distortion measurements for a 15 deg clocked 40 probe rake indicated that the circumferential distortion levels are very sensitive to the symmetry of the flow and that a misalignment of the vanes in the experiment could have resulted in this difference. The numerical simulations of the BLI inlet with jets showed good agreement with the circumferential inlet distortion levels for a range of jet actuator mass flow ratios at a fixed inlet mass flow rate. The CFD simulations for the jet case also predicted an average total pressure recovery offset that was 0.01 lower than the experiment as was seen in the baseline. Comparisons of the flow features for the jet cases revealed that the CFD predicted a much larger vortex at the engine fan-face when compare to the experiment.

Allan, Brian G.↗

DOE OSTI2025

Advection algorithms for quantum neutrino moment transport

Neutrino transport in compact objects is an inherently challenging multidimensional problem. Here, this difficulty is compounded if one includes flavor transformation—an intrinsically quantum phenomenon requiring one to follow the coherence between flavors and thus necessitating the introduction of complex numbers. To reduce the computational burden, simulations of compact objects that include neutrino transport often make use of momentum-angle-integrated moments (the lowest order ones being commonly referred to as the energy density and flux) and these quantities can be generalized to include neutrino flavor, i.e., they become quantum moments. Numerous finite-volume approaches to solving the moment evolution equations for classical neutrino transport have been developed based on solving a Riemann problem at cell interfaces. In this paper we describe our generalization of a Riemann solver for quantum moments, specifically decomposing complex numbers in terms of a (signed) magnitude and phase instead of real and imaginary parts. We then test our new algorithm in numerous cases showing a neutrino fast flavor instability, varying from toy models with analytic solutions to snapshots from neutron star merger simulations. Compared to previous algorithms for neutrino transport with flavor mixing, we find uniformly smaller growth rates of the flavor transformation along with concomitantly larger length-scales, and that the results are a better match with the growth rates seen from multiangle codes.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2006

Numerical Modeling of Flow Control in a Boundary-Layer-Ingesting Offset Inlet Diffuser at Transonic Mach Numbers

This paper will investigate the validation of a NASA developed, Reynolds-averaged Navier-Stokes (RANS) flow solver, OVERFLOW, for a boundary-layer-ingesting (BLI) offset (S-shaped) inlet in transonic flow with passive and active flow control devices as well as the baseline case. Numerical simulations are compared to wind tunnel results of a BLI inlet conducted at the NASA Langley 0.3-Meter Transonic Cryogenic Tunnel. Comparisons of inlet flow distortion, pressure recovery, and inlet wall pressures are performed. The numerical simulations are compared to the BLI inlet data at a freestream Mach number of 0.85 and a Reynolds number of approximately 2 million based on the length of the fan-face diameter. The numerical simulations with and without wind tunnel walls are performed, quantifying effects of the tunnel walls on the BLI inlet flow measurements. The wind tunnel test evaluated several different combinations of jet locations and mass flow rates as well as a vortex generator (VG) vane case. The numerical simulations will be performed on a single jet configuration for varying actuator mass flow rates at a fix inlet mass flow condition. Validation of the numerical simulations for the VG vane case will also be performed for varying inlet mass flow rates. Overall, the numerical simulations were able to predict the baseline circumferential flow distortion, DPCPavg, very well for comparisons made within the designed operating range of the BLI inlet. However the CFD simulations did predict a total pressure recovery that was 0.01 lower than the experiment. Numerical simulations of the baseline inlet flow also showed good agreement with the experimental inlet centerline surface pressures. The vane case showed that the CFD predicted the correct trends in the circumferential distortion for varying inlet mass flow but had a distortion level that was nearly twice as large as the experiment. Comparison to circumferential distortion measurements for a 15 deg clocked 40 probe rake indicated that the circumferential distortion levels are very sensitive to the symmetry of the flow and that a miss alignment of the vanes in the experiment could have resulted in this difference. The numerical simulations of the BLI inlet with jets showed good agreement with the circumferential inlet distortion levels for a range of jet actuator mass flow ratios at a fixed inlet mass flow rate. The CFD simulations for the jet case also predicted an average total pressure recovery that was 0.01 lower than the experiment as was seen in the baseline. Comparison of the flow features the jet case revealed that the CFD predicted a much larger vortex at the engine fan-face when compare to the experiment.

Allan Brian G.↗