SEARCH · Search NASA

Results for “Convex geometry”

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 19 records

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Quantum simulation of quantum phase transitions using the convex geometry of reduced density matrices

Transitions of many-particle quantum systems between distinct phases at absolute-zero temperature, known as quantum phase transitions, require an exacting treatment of particle correlations. Here in this work, we present a general quantum-computing approach to quantum phase transitions that exploits the geometric structure of reduced density matrices. While typical approaches to quantum phase transitions examine discontinuities in the order parameters, the origin of phase transitions—their order parameters and symmetry breaking—can be understood geometrically in terms of the set of two-particle reduced density matrices (2-RDMs). The convex set of 2-RDMs provides a comprehensive map of the quantum system including its distinct phases as well as the transitions connecting these phases. Because 2-RDMs can potentially be computed on quantum computers at nonexponential cost, even when the quantum system is strongly correlated, they are ideally suited for a quantum-computing approach to quantum phase transitions. We compute the convex set of 2-RDMs for a Lipkin-Meshkov-Glick spin model on IBM superconducting-qubit quantum processors. Even though computations are limited to few-particle models due to device noise, comparisons with a classically solvable 1000-particle model reveal that the finite-particle quantum solutions capture the key features of the phase transitions including the strong correlation and the symmetry breaking.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2024

The geometry of the modular bootstrap

Abstract We explore the geometry behind the modular bootstrap and its image in the space of Taylor coefficients of the torus partition function. In the first part, we identify the geometry as an intersection of planes with the convex hull of moment curves onR + ⊗ℤ, with boundaries characterized by the total positivity of generalized Hankel matrices. We phrase the Hankel constraints as a semi-definite program, which has several advantages, such as the validity of bounds irrespective of spin truncation. We derive bounds on the gap, twist-gap, and the space of Taylor coefficients themselves. We find that if the gap is above$$ {\Delta }_{\textrm{gap}}^{\ast } $$ ∆ gap ∗ , where$$ \frac{c-1}{12}<{\Delta}_{\textrm{gap}}^{\ast }<\frac{c}{12} $$ c − 1 12 < Δ gap ∗ < c 12 , all coefficients become bounded on both sides and kinks develop in the space. In the second part, we propose an analytic method of imposing the integrality condition for the degeneracy number in the spinless bootstrap, which leads to a non-convex geometry. We find that even at very low derivative order this condition rules out regions otherwise allowed by bootstraps at high derivative order.

Physics↗

DOE OSTI2021

Upgrades and redesign of the National Ignition Facility’s soft x-ray opacity spectrometer (OpSpec)

The soft x-ray Opacity Spectrometer (OpSpec) used on the National Ignition Facility (NIF) has recently incorporated an elliptically shaped crystal. The original OpSpec used two convex cylindrical crystals for time-integrated measurements of point-projection spectra from 540 to 2100 eV. However, with the convex geometry, the low-energy portion of the spectrum suffered from high backgrounds due to scattered x-rays as well as reflections from alternate crystal planes. An elliptically shaped crystal allows an acceptance aperture at the crossover focus between the crystal and the detector, which reduces background and eliminates nearly all reflections from alternate crystal planes. The current elliptical design is an improvement from the convex cylindrical design but has a usable energy range from 900 to 2100 eV. In addition, OpSpec is currently used on 18 NIF shots/year, in which both crystals are typically damaged beyond reuse, so efficient production of 36 crystals/year is required. Design efforts to improve the existing system focus on mounting reliability, reducing crystal strain to increase survivability between mounting and shot time, and extending the energy range of the instrument down to 520 eV. The elliptical design, results, and future options are presented.

Wallace, M. S. (ORCID:0000000310871216)↗

DOE OSTI2022

Voronoi cell analysis: The shapes of particle systems

Many physical systems can be studied as collections of particles embedded in space, often evolving in time. Natural questions arise concerning how to characterize these arrangements—are they ordered or disordered? If they are ordered, how are they ordered and what kinds of defects do they possess? Voronoi tessellations, originally introduced to study problems in pure mathematics, have become a powerful and versatile tool for analyzing countless problems in pure and applied physics. Here we explain the basics of Voronoi tessellations and the shapes that they produce and describe how they can be used to characterize many physical systems.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

A discretized representation for Monte Carlo simulation of deformed semiflexible chains

In this study, we present a novel orientation discretization approach based on the rhombic triacontahedron for Monte Carlo simulations of semiflexible polymer chains, aiming at enhancing structural analysis through rheo-small-angle scattering (rheo-SAS). Our approach provides a more accurate representation of the geometric features of semiflexible chains under deformation, surpassing the capabilities of traditional lattice structures. Validation against the Kratky–Porod chain system demonstrated superior consistency, underscoring its potential to significantly improve the precision of uncovering geometric details from rheo-SAS data. Further, this approach opens new avenues for investigating the conformations of semiflexible polymers and mechanically induced phase transitions in more complex polymer structures, offering deeper insights into their behavior under various conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

A three-dimensional laser ray-tracing methodology for radiation-hydrodynamics simulations

We report on a methodology for performing laser ray-tracing in three spatial dimensions for radiation-hydrodynamics simulation codes. Our method, which is an extension of that developed in Haines et al., Comput. Fluids 201, 104478 (2020), utilizes an automatically generated separate mesh for the laser ray-tracing from the radiation-hydrodynamics mesh. This enables the laser mesh to be tailored to minimize ray noise with significantly fewer rays than would be required when the ray-tracing is performed on the radiation-hydrodynamics mesh, primarily by allowing the use of high-aspect-ratio cells that are not suitable for hydrodynamics solvers. For a planar target, we show that our method provides a ≈ 100× reduction in computational expense to achieve a fixed level of ray noise relative to ray-tracing directly on the radiation-hydrodynamics mesh. The relatively low ray requirement also enables efficient computation of cross-beam energy transfer. Each cell in the logically cubic laser mesh is a non-convex dodecahedron with triangular sides, and numerical integration of the ray trajectories and inverse bremsstrahlung is performed by mapping each cell to the unit cube. We will describe our methodology in detail as well as its implementation in the xRAGE radiation-hydrodynamics code, discuss performance, and present the results from applying the methodology to test problems with analytic solutions for laser ray-tracing through a quadratic density gradient with an analytic solution as well as for a laser-driven heat front. In 3D radiation-hydrodynamics simulations of laser-driven experiments performed on the National Ignition Facility, laser ray-tracing with our methodology uses less than 1% of total computational time while introducing acceptably low levels of ray noise.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

A plano–convex thick-lens velocity map imaging apparatus for direct, high resolution 3D momentum measurements of photoelectrons with ion time-of-flight coincidence

Since their inception, velocity map imaging (VMI) techniques have received continued interest in their expansion from 2D to 3D momentum measurements through either reconstructive or direct methods. Recently, much work has been devoted to the latter of these by relating electron time-of-flight (TOF) to the third momentum component. The challenge is having a timing resolution sufficient to resolve the structure in the narrow (<10 ns) electron TOF spread. Here, we build upon the work in VMI lens design and 3D VMI measurement by using a plano–convex thick-lens (PCTL) VMI in conjunction with an event-driven camera (TPX3CAM) providing TOF information for high resolution 3D electron momentum measurements. We perform simulations to show that, with the addition of a mesh electrode to the thick-lens geometry, the resulting plano–convex electrostatic field extends the detectable electron cutoff energy range while retaining the high resolution. This design also extends the electron TOF range, allowing for a better momentum resolution along this axis. We experimentally demonstrate these capabilities by examining above-threshold ionization in xenon, where the apparatus is shown to collect electrons of energy up to ∼7 eV with a TOF spread of ∼30 ns, both of which are improved compared to a previous work by factors of ∼1.4 and ∼3.75, respectively. Finally, the PCTL-VMI is equipped with a coincident ion TOF spectrometer, which is shown to effectively extract unique 3D momentum distributions for different ionic species in a gas mixture. These techniques have the potential to lend themselves to more advanced measurements involving systems where the electron momentum distributions possess non-trivial symmetries.

Instruments & Instrumentation↗

DOE OSTI2021

Analysis of the ratio of ℓ 1 and ℓ 2 norms in compressed sensing

We study the ratio of ℓ 1 and ℓ 2 norms ( ℓ 1 / ℓ 2 ) as a sparsity-promoting objective in compressed sensing. We first propose a novel criterion that guarantees that an s-sparse signal is the local minimizer of the ℓ 1 / ℓ 2 objective; our criterion is interpretable and useful in practice. We also give the first uniform recovery condition using a geometric characterization of the null space of the measurement matrix, and show that this condition is satisfied for a class of random matrices. We also present analysis on the robustness of the procedure when noise pollutes data. Numerical experiments are provided that compare ℓ 1 / ℓ 2 with some other popular non-convex methods in compressed sensing. Finally, we propose a novel initialization approach to accelerate the numerical optimization procedure. We call this initialization approach support selection, and we demonstrate that it empirically improves the performance of existing ℓ 1 / ℓ 2 algorithms.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Data-driven, structure-preserving approximations to entropy-based moment closures for kinetic equations

In this study, we present a data-driven approach for approximating entropy-based closures of moment systems from kinetic equations. The proposed closure learns the entropy function by fitting the map between the moments and the entropy of the moment system, and thus does not depend on the spacetime discretization of the moment system or specific problem configurations such as initial and boundary conditions. With convex and C 2 approximations, this data-driven closure inherits several structural properties from entropy-based closures, such as entropy dissipation, hyperbolicity, and H-Theorem. We construct convex approximations to the Maxwell–Boltzmann entropy using convex splines and neural networks, test them on the plane source benchmark problem for linear transport in slab geometry, and compare the results to the standard, entropy-based systems which solve a convex optimization problem to find the closure. Numerical results indicate that these data-driven closures provide accurate solutions in much less computation time than that required by the optimization routine.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Polyelectrolyte Surface Diffusion in a Nanoslit Geometry

The surface diffusion of poly-l-lysine (PLL) in a planar nanoslit was studied using convex lens-induced confinement (CLiC) single-molecule tracking microscopy. Here, three surface chemistries were employed to understand the interplay of electrostatic and short-range interactions: an amine-functionalized silica surface, an oligo(ethylene oxide) (OEG)-modified surface, and a 1:1 mixture of the two ligands. Effective surface diffusion coefficients increased rapidly with slit height until saturating for slit heights <30 nm. While diffusion at a semi-infinite interface was significantly faster for OEG surfaces, the diffusion coefficient increased most rapidly with slit height for amine-functionalized surfaces, resulting in surface diffusion within very thin slits being nearly independent of surface chemistry. Intermittent random walks were simulated within a planar slit geometry, using experimentally measured parameters obtained from diffusion at a single interface to account for the characteristic short-range interactions between PLL and each surface chemistry, and were in good agreement with experimental measurements.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Chiral Kirigami for Bend‐Tolerant Reconfigurable Hologram with Continuously Variable Chirality Measures

Abstract Despite the commonality of static holograms, the holography with multiple information layers and reconfigurable grey‐scale images at communication frequencies remain a confluence of scientific challenges. One well‐known difficulty is the simultaneous modulation of phase and amplitude of electromagnetic wavefronts with a high modulation depth. A less appreciated challenge is scrambling of the information and images with hologram bending. Here, this work shows that chirality‐guided pixelation of plasmonic kirigami sheets enables tunable multiplexed holography at terahertz (THz) frequencies. The convex and concave structures with slanted Au strips exhibit gradual variations in geometries facilitating modulation of light ellipticity reaching 40 deg. Real‐time switching of 3D images of the letter “ M ” and the Mona Lisa demonstrates the possibility of complex grey‐scale information content and importance of continuously variable mirror asymmetry. Microscale chirality measures of each pixel experiences little change with bending while retaining controllable reconfigurability upon stretching, which translates to remarkable resilience of chiral holograms to bending. Simplicity of their design with local chirality measures opens the door to information technologies with fault‐tolerant THz encryption, wearable holographic devices, and new communication technologies.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Consequences of Convex Nanopore Chemistry on Confined Water Dynamics

A basic understanding of confined water is crucial for developing selective ion transport and water purification membranes, yet the roles of nanopore geometry and functionality on confined water dynamics remain unresolved. We discuss the synthesis of perdeuterated ionic alkylsulfonate amphiphiles and their water-induced self-assembly into lyotropic liquid crystal (LLC) mesophases with well-defined, convex, sulfonate-lined nanopores. Quasielastic neutron scattering (QENS) measurements demonstrate that the water self-diffusion coefficients within these sulfonate-lined convex nanopores depend on the hydration level and amphiphile counterion identity (H + , K + , NMe 4 + ). The consistency of the observed counterion-dependent water dynamics trends with those of carboxylate LLCs is rationalized on the basis of similarities in the counterion spatial distributions in the water-filled channels, which we deduce from electron density maps derived from small-angle X-ray scattering (SAXS) analyses. These findings indicate that water diffusion is systematically faster in sulfonate-lined nanopores as compared to carboxylate-lined pores due to weaker water interactions with the softer and more hydrophobic–SO 3 – functionalities. These molecular-level insights into the relationships between convex pore wall chemical functionalities, hydrated counterions, and confined water diffusion may inform future development of new nanoporous media.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Computational optimal transport for molecular spectra: The fully continuous case

Computational optimal transport is used to analyze the difference between pairs of continuous molecular spectra. It is demonstrated that transport distances which are derived from this approach may be a more appropriate measure of the difference between two continuous spectra than more familiar measures of distance under many common circumstances. Associated with the transport distances is the transport map which provides a detailed analysis of the difference between two molecular spectra and is a key component of our study of quantitative differences between two continuous spectra. The use of optimal transport for comparing molecular spectra is developed in detail here with a set of model spectra, so that the discussion is self-contained. The difference between the transport distance and more common definitions of distance is elucidated for some well-chosen examples and it is shown where transport distances may be very useful alternatives to standard definitions of distance. The transport distance between a theoretical and experimental electronic absorption spectrum for SO 2 is studied and it is shown how the theoretical spectrum can be modified to fit the experimental spectrum better adjusting the theoretical band origin and the resolution of the theoretical spectrum. In conclusion, this analysis includes the calculation of transport maps between the theoretical and experimental spectra suggesting future applications of the methodology.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2026

Bio-inspired alula-based winglet design for enhanced heat transfer in high temperature fin-and-tube heat exchangers

Fin-and-tube heat exchangers (FTHEs) are widely used for high-temperature flue-gas heat recovery, but their performance is often limited by wake regions and non-uniform fin-surface temperatures. This study proposes and numerically evaluates four bio-inspired longitudinal vortex generator (VG) configurations in a high-temperature FTHE with flue-gas inlet temperature ∼1230 K: double-delta, curved double-delta, alula, and a new curved-alula geometry. The reference fin is not hydraulically plain; it already incorporates leading-edge separation columns and convex protrusions, so the alula-type winglets are assessed as downstream add-ons acting on a strongly disturbed flow. In a second step, perforations (one, two and three circular holes) are introduced into the curved-alula VGs to further tailor the flow field. Three-dimensional simulations with the Shear Stress Transpor (SST) $k - ω$ model, temperature-dependent flue-gas properties and conjugate conduction are carried out for gas-side Reynolds numbers $Re_g ≈ 8.0$ x $10^2 - 3.6$ x $10^3$ (mass flow rates 0.5 – 2.5 g/s), and the designs are compared in terms of surface heat flux, Nusselt number, friction factor and hydrothermal performance factor (HTPF). For this already-promoted fin, the additional downstream winglets provide moderate, incremental hydrothermal gains. At the highest Reynolds number, the best non-perforated design (curved-alula) increases surface heat flux from 1630.9 to 1794.7 kW/m² (∼ 10 % gain) and the Nusselt number from 227.6 to 242.6 (∼ 7 % gain), while the friction factor rises from 0.26 to about 0.30, yielding HTPF values close to unity (∼ 0.9 – 1.0). Introducing circular perforations into the curved-alula winglets acts mainly as a wake-bleeding refinement: the three-hole configuration provides a heat flux of 1824.7 kW/m² and a pressure drop of 127.9 Pa, with HTPF in the range ∼ 1.03 – 1.14 and a small (∼ 1 – 3 %) improvement over the solid curved-alula design. Flow-field analysis shows that the perforated curved-alula VGs shrink tube-wake regions, thin the thermal boundary layer and homogenize the fin-surface temperature (outlet-gas temperature ∼ 510 – 520 K and fin-surface temperature ∼ 420 – 421 K for the three-hole case). An optimal flue-gas mass flow rate of ∼ 1 g/s ($Re_g ≈ 1.5$ x $10^3$) is identified, beyond which additional heat-transfer gains are offset by rapidly increasing pressure losses. Overall, the results highlight that initial fin geometry and VG placement are as important as VG shape: alula-based winglets are expected to yield larger relative gains on simpler flat-fin layouts or when positioned closer to the fin leading edge and tube

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2022

GEOMWRITER - A Workflow for Automatically Generating MCNP Geometry for Neutron Guide Simulations

This report describes the functionality of the geomwriter.py script, which is a python script for producing complex MCNP CSG geometries from the IGES CAD format. It is specifically written to handle neutron supermirror guides, but can handle more complex objects through utilizing a convex hull algorithm. This work flow reduces the amount of time spent on generating the surface definitions for a complex neutron guide system with hundreds of segments such that changes to the neutron guide system design in McStas can be rapidly incorporated into an MCNP model and integrated with the shielding geometry present in that model.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2020

Oxidation of 3D-printed SiC in air and steam environments

Here, the high-temperature oxidation of additively manufactured and chemically vapor infiltrated (3D-printed SiC) has been compared to chemical vapor deposited (CVD) SiC. 100-h isothermal exposures were conducted at 1425° and 1300°C at 1 atm under both dry air and steam environments. A SiC reaction tube was utilized to reduce silica volatility. After steam oxidation at 1425° and 1300°C, on the 3D-printed SiC surface, which was intrinsically rougher than the CVD surface, scales were 70%–90% thicker at the convex regions compared to concave/flat regions. In the convex regions, large cracks perpendicular to the oxidizing interface were observed. After dry air oxidation, scale thicknesses were comparable between 3D-printed SiC and CVD SiC, regardless of geometry. Finite element modeling, conducted to elucidate the relationship between SiC geometry and ß- to α-cristobalite transformation stress, determined cristobalite transformation tensile stresses to be on the order of 103 MPa during cool down, assuming a 6 vol% reduction. Compared to flat SiC substrates, tensile transformation stresses were elevated at concave regions and relaxed at convex regions. Combined with specimen mass gain (accounting for the rougher surface) of 3D-printed SiC being 15%–32% higher for 3D-printed SiC after 1300°C and 1425°C steam oxidation, the work presented concludes that the increased oxidation of 3D-printed SiC is primarily caused by tensile hoop stresses driven by oxidation volume expansion. Lastly, the efficacy of the 3D-printing method is demonstrated through the production of tristructural isotropic imbedded 3D-printed SiC fuel forms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Convex Programs for Minimal-Area Problems

The minimal-area problem that defines string diagrams in closed string field theory asks for the metric of least area on a Riemann surface with the condition that all non-contractible closed curves have length at least 2π. This is an extremal length problem in conformal geometry as well as a problem in systolic geometry. Here, we consider the analogous minimal-area problem for homology classes of curves and, with the aid of calibrations and the max flow-min cut theorem, formulate it as a local convex program. We derive an equivalent dual program involving maximization of a concave functional. These two programs give new insights into the form of the minimal-area metric and are amenable to numerical solution. We explain how the homology problem can be modified to provide the solution to the original homotopy problem.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Curvature effects in the reaction-diffusion system governing frontal polymerization

Frontal polymerization is a fast and energy-efficient curing process of thermoset polymers involving the self-sustained propagation of a sharp front that transforms a liquid monomer or gel into a polymer. The temperature, speed, and geometry of the polymerization front are closely intertwined quantities that can affect the curing process and the mechanical properties of the polymer. This work investigates the effect of the front curvature on the front propagation speed based on a reaction-diffusion model. The proposed analytical relations between front speed and curvature are verified by finite element simulations and validated through experiments. The models show that a concave front has a higher front temperature and speed than its convex counterpart. When the curvature of a convex front reaches a critical value, the reaction can no longer be sustained, and the front is quenched. The relation between the speed of a steadily propagating front and its geometry is derived, showing good agreement with experimental observations. This study also provides analytical estimates of the thickness of the polymerization front and insight into the relationship between the front curvature and the kinetics of the exothermic reaction.

42 ENGINEERING↗