SEARCH · Search NASA

Results for “CANONICAL PROBLEM”

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 181 records · Page 10

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Reproducibility of fixed-node diffusion Monte Carlo across diverse community codes: The case of water–methane dimer

Fixed-node diffusion quantum Monte Carlo (FN-DMC) is a widely trusted many-body method for solving the Schrödinger equation, known for its reliable predictions of material and molecular properties. Furthermore, its excellent scalability with system complexity and near-perfect utilization of computational power make FN-DMC ideally positioned to leverage new advances in computing to address increasingly complex scientific problems. Even though the method is widely used as a computational gold standard, reproducibility across the numerous FN-DMC code implementations has yet to be demonstrated. This difficulty stems from the diverse array of DMC algorithms and trial wave functions, compounded by the method’s inherent stochastic nature. Here, this study represents a community-wide effort to assess the reproducibility of the method, affirming that yes, FN-DMC is reproducible (when handled with care). Using the water–methane dimer as the canonical test case, we compare results from eleven different FN-DMC codes and show that the approximations to treat the non-locality of pseudopotentials are the primary source of the discrepancies between them. In particular, we demonstrate that, for the same choice of determinantal component in the trial wave function, reliable and reproducible predictions can be achieved by employing the T-move, the determinant locality approximation, or the determinant T-move schemes, while the older locality approximation leads to considerable variability in results. These findings demonstrate that, with appropriate choices of algorithmic details, fixed-node DMC is reproducible across diverse community codes—highlighting the maturity and robustness of the method as a tool for open and reliable computational science.

Della Pia, Flaviano [Univ. of Cambridge (United Ki↗

NASA NTRS1985

Research on output feedback control

In designing fixed order compensators, an output feedback formulation has been adopted by suitably augmenting the system description to include the compensator states. However, the minimization of the performance index over the range of possible compensator descriptions was impeded due to the nonuniqueness of the compensator transfer function. A controller canonical form of the compensator was chosen to reduce the number of free parameters to its minimal number in the optimization. In the MIMO case, the controller form requires a prespecified set of ascending controllability indices. This constraint on the compensator structure is rather innocuous in relation to the increase in convergence rate of the optimization. Moreover, the controller form is easily relatable to a unique controller transfer function description. This structure of the compensator does not require penalizing the compensator states for a nonzero or coupled solution, a problem that occurs when following a standard output feedback synthesis formulation.

Calise, A. J.↗

NASA NTRS1975

On the problems associated with the rapidly varying compact nonthermal sources

Difficulties and possible interpretations are discussed for observations of extragalactic compact nonthermal radio sources. The sources considered are those which have been identified optically with normal and peculiar galaxies, quasars, N galaxies, and BL Lacertae objects. Optical, infrared, and radio flux variations in different sources are noted along with difficulties associated with determining angular sizes from interferometric measurements. It is suggested that the ultimate origin of the infrared flux is probably a nonthermal synchrotron mechanism in some sources and absorption and reradiation of ultraviolet synchrotron radiation by dust in others. Model calculations of compact sources are reviewed, emphasizing equipartition of energy, total energies, the importance of the Compton effect, and sources which present severe difficulties for the theory of a canonical incoherent electron-synchrotron source.

Burbidge, G. R.↗

DOE OSTI2026

From Quantum Time to Manifestly Covariant QFT: On the Need for a Quantum-Action-Based Quantization

In quantum time (QT) schemes, time is promoted to a degree of freedom, allowing Lorentz covariance to be made explicit for single particles. We ask whether this can be lifted to QFT so that Lorentz covariance becomes manifest at the Hilbert-space level, rather than being hidden as in the standard canonical formulation. We address this question by proposing a second-quantized approach in which the elementary particle is the QT particle itself, leading naturally to the notion of spacetime field algebras and of quantum action. We show, however, that a naive many-body construction runs into inconsistencies. To pinpoint their origin we introduce a classical counterpart of the second-quantized formalism, spacetime classical mechanics (SCM), and prove a no-go theorem: Dirac quantization of SCM collapses back to standard QFT and therefore hides covariance. We circumvent this problem by presenting a quantum-action-based quantization that yields a spacetime version of quantum mechanics (SQM), making covariance manifest for (interacting) QFTs. Finally, we show that this resolution is tied to a genuine spacetime generalization of the notion of a quantum state, required by causality and closely connected to recent “states over time” proposals and, in dS/CFT–motivated settings, to microscopic notions of timelike entanglement and emergent time.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

Testing the scalar weak gravity conjecture in no-scale supergravity

We explore possible extensions of the Weak Gravity Conjecture (WGC) to scalar field theories. To avoid charged black hole remnants, the WGC requires the existence of a particle with a mass m ≤ gqM P , with charge q and U(1) gauge coupling g, allowing the decay to shed the black hole charge. Although there is no obvious problem that arises in the absence of a U(1) charge, it has been postulated that gravity must remain the weakest force even when extended to scalar interactions. Quantifying this conjecture may be done by comparing scalar and gravitational amplitudes, or as we advocate here by comparing scattering cross sections. In theories with non-trivial field space geometries, by working out examples with perturbation theory around arbitrary field values and performing tadpole resummations, we argue that the conjecture must be applied only at extrema of the scalar potential (when expressed in locally canonical coordinates). We consider several toy models in the context of no-scale supergravity and also consider examples of inflationary models.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2025

Integrated Molten Salt Reactor Modeling Capabilities in NEAMS Thermal Hydraulics Tools

The DOE neams program supports a full range of computational thermal fluids analysis capabilities and code developments for a broad range of advanced reactor concepts. The research and development approach under the thermal fluids technical area synergistically combines three length and time scales in a hierarchical multi-scale approach. To enable multi-scale thermal fluids capability using these codes, a key joint effort has been underway to develop an integrated system- and engineering-scale thermal fluids analysis capability, through integration of SAM and Pronghorn codes, both based on the MOOSE framework. This report summarizes recent advances in developing an integrated system- and engineering-scale modeling capability for the msr concept, which has gained significant interest in recent years. A consistent framework was established by coupling Pronghorn and SAM through the Saline interface, with thermophysical properties provided by the Molten Salt Thermal Property Database (MSTDB-TP). Further improvements were made to the coupling schemes and domain-overlapping strategies, enhancing the stability and robustness of multi-code simulations. Verification and validation efforts demonstrate the accuracy of this integration across a range of benchmark problems, including one-dimensional heated pipe flows, three-dimensional natural convection loops with evolving isotopic compositions, and \gls{msre} demonstration cases. Within Pronghorn, new capabilities were introduced to model corrosion and noble-metal plating phenomena, supported by an extended thermal-hydraulics framework and refined turbulence treatments. To capture two-phase flow behavior, a multiphase Euler–Euler model was implemented in Pronghorn, including advanced closure relations, high-resolution advection techniques, and capillary force reconstruction. Preliminary verification cases confirm the fidelity of the approach, while planned validation efforts target canonical multiphase benchmarks and application to msr components such as the msre pump bowl. Finally, updates to SAM’s msr mass transfer modeling were extended to consider noble gas migration into porous structures like graphite. The point kinetics model was updated to include reactivity feedback contributions from any defined species, such as xenon. The gas transport model was expanded for applicability to gas mixtures, bubble efflux phenomena, and species transport between liquid and gas phases. A selection of multi-scale Sherwood number correlations from MOSCATO/NekRS and multi-phase correlations from literature have been added for improved accuracy in calculating mass transfer coefficients. A companion effort on developing system-level redox corrosion has also been incorporated into SAM. Collectively, these enhancements strengthen the predictive capability of SAM and Pronghorn for simulating MSR thermal-hydraulics, corrosion, multiphase behavior, and fission-product transport, providing a more complete toolset for design, safety analysis, and licensing support of next-generation \gls{msr}s.

42 - ENGINEERING↗

DOE OSTI2026

Efficient CP Rounding Using Alternating Least Squares with QR Decomposition

The CANDECOMP/PARAFAC (CP) decomposition is widely used for analyzing multidimensional data, and the alternating least squares (CP-ALS) algorithm is a common method for its computation. CP rounding is the problem of computing a lower-rank CP decomposition of an input already in a higher-rank CP format. While the normal equations (NE) approach in CP-ALS is efficient for the CP rounding problem and frequently used, it becomes unstable in the presence of ill-conditioned subproblems. This paper presents a new QR-based CP-ALS method for CP rounding that preserves both numerical stability and computational efficiency. Here, our experiments show that the proposed method offers significant speedup over a previous QR-based approach and the Tensor Toolbox's NE-based implementation, particularly for higher-order tensors. Furthermore, our approach demonstrates a marked reduction in error for ill-conditioned problems, with error reductions several orders of magnitude smaller compared to the NE-based method, while achieving faster convergence and more accurate solutions. By using a more numerically stable approach, we can solve more problems in reduced working precision, which enables further reduction in time to solution.

CANDECOMP/PARAFAC↗

DOE OSTI2025

Nonenzymatic RNA copying with a potentially primordial genetic alphabet

Nonenzymatic RNA copying is thought to have been responsible for the replication of genetic information during the origin of life. However, chemical copying with the canonical nucleotides (A, U, G, and C) strongly favors the incorporation of G and C and disfavors the incorporation of A and especially U because of the stronger G:C vs. A:U base pair and the weaker stacking interactions of U. Recent advances in prebiotic chemistry suggest that the 2-thiopyrimidines were precursors to the canonical pyrimidines, raising the possibility that they may have played an important early role in RNA copying chemistry. Furthermore, 2-thiouridine (s 2 U) and inosine (I) form by deamination of 2-thiocytidine (s 2 C) and A, respectively. We used thermodynamic and crystallographic analyses to compare the I:s 2 C and A:s 2 U base pairs. We find that the I:s 2 C base pair is isomorphic and isoenergetic with the A:s 2 U base pair. The I:s 2 C base pair is weaker than a canonical G:C base pair, while the A:s 2 U base pair is stronger than the canonical A:U base pair, so that a genetic alphabet consisting of s 2 U, s 2 C, I, and A generates RNA duplexes with uniform base pairing energies. Consistent with these results, kinetic analysis of nonenzymatic template-directed primer extension reactions reveals that s 2 C and s 2 U substrates bind similarly to I and A in the template, and vice versa. Our work supports the plausibility of a potentially primordial genetic alphabet consisting of s 2 U, s 2 C, I, and A and offers a potential solution to the long-standing problem of biased nucleotide incorporation during nonenzymatic template copying.

Science & Technology - Other Topics↗

DOE OSTI2025

Scalar flux transport models for self-similar turbulent mixing

A common approach to closing turbulent species flux in multicomponent Reynolds-averaged Navier-Stokes models is to use the standard gradient diffusion approximation. While such an approach has been shown to work well when applied to many canonical turbulent mixing configurations, a gradient diffusion approach is fundamentally limited in its ability to capture complex phenomena such as countergradient transport. For this reason, complicated mixing applications may benefit by treating the turbulent diffusivity with a model transport equation in a manner analogous to second-moment momentum closure in Reynolds-stress transport models. Here, the present work explores the development and application of two different scalar flux transport (SFT) models. Self-similarity constraints are derived for these models, and they are evaluated against gradient-diffusion-based models in several one- and two-dimensional problems of turbulent mixing. It is found that the new SFT models out-perform gradient diffusion models in problems involving rapid acceleration reversal and in problems involving anisotropic transport of materials. In addition, it is found that even a hybrid-SFT approach, in which an SFT equation is utilized along with a gradient diffusion closure, provides some measure of improvement over models that transport the mass flux rather than the scalar flux.

Reynolds-averaged Navier Stokes↗

DOE OSTI2024

The Average Spectrum Norm and Near-Optimal Tensor Completion

We propose the average spectrum norm to study the minimum number of measurements required to approximate a multidimensional array (i.e., sample complexity) via low-rank tensor recovery. Our focus is on the tensor completion problem, where the aim is to estimate a multiway array using a subset of tensor entries corrupted by noise. Our average spectrum norm-based analysis provides near-optimal sample complexities, exhibiting dependence on the ambient dimensions and rank that do not suffer from exponential scaling as the order increases.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1994

A Modular, Portable Model of Image Fidelity

There is a persistent need for a trustworthy model of perceptual image fidelity, especially in applications such as image compression and display design. A fidelity model provides a measure of the visual discriminability of two images. Ahumada has previously shown that the existing fidelity models may be categorized according to their inclusion of various canonical properties, such as a contrast sensitivity function, spatial frequency channels, etc. This suggests that research would be aided by the availability of a modular model, in which these components could be easily inserted or removed. A further impediment to research in this area has been that most models are written in low-level languages and are consequently large, non-portable, and difficult to understand, modify, and maintain. We therefore believe research would also be aided by models written in high-level languages. To serve both of these purposes, and to honor our conference host for his lifetime dedication to the problem of image quality. Global brightness and its effect on perceptual image quality. We offer a modular model written in the high-level language Mathematica. We will demonstrate this model and show how it may be modified.

Watson, Andrew B.↗

NASA NTRS1995

Safety analysis

We are engaged in a research program in safety-critical computing that is based on two case studies. We use these case studies to provide application-specific details of the various research issues, and as targets for evaluation of research ideas. The first case study is the Magnetic Stereotaxis System (MSS), an investigational device for performing human neurosurgery being developed in a joint effort between the Department of Physics at the University of Virginia and the Department of Neurosurgery at the University of Iowa. The system operates by manipulating a small permanent magnet (known as a 'seed') within the brain using an externally applied magnetic field. By varying the magnitude and gradient of the external magnetic field, the seed can be moved along a non-linear path and positioned at a site requiring therapy, e.g., a tumor. The magnetic field required for movement through brain tissue is extremely high, and is generated by a set of six superconducting magnets located in a housing surrounding the patient's head. The system uses two X-ray cameras positioned at right angles to detect in real time the locations of the seed and of X-ray opaque markers affixed to the patient's skull. the X-ray images are used to locate the objects of interest in a canonical frame of reference. the second case study is the University of Virginia Research Nuclear Reactor (UVAR). It is a 2 MW thermal, concrete-walled pool reactor. The system operates using 20 to 25 plate-type fuel assemblies placed on a rectangular grid plate. There are three scramable safety rods, and one non-scramable regulating rod that can be put in automatic mode. It was originally constructed in 1959 as a 1 MW system, and it was upgraded to 2 MW in 1973. Though only a research reactor rather than a power reactor, the issues raised are significant and can be related to the problems faced by full-scale reactor systems.

Knight, John C.↗

NASA NTRS1974

Optimal compensator structure for linear time-invariant plant with inaccessible states

The problem is considered of designing an optimal linear time-invariant dynamic compensator for the regulation of an n-th order linear time-invariant plant with m independent outputs. The initial plant state is characterized by its first and second moments, and the cost is usual quadratic infinite-time penalty on the state and control, averaged over the initial plant and compensator states. The compensator is based on a minimal-order Luenberger observer and consequently has fixed dimension n-m. Necessary and sufficient conditions are derived for optimality of the compensator gains. The optimal compensator is shown to be unique if the plant has a particular canonical form, and, in general, for any arbitrary plant, the class of all optimal compensators is precisely determined.

Blanvillain, P. J. P.↗

NASA NTRS1994

A uniform geometrical optics and an extended uniform geometrical theory of diffraction for evaluating high frequency EM fields near smooth caustics and composite shadow boundaries

A uniform geometrical optics (UGO) and an extended uniform geometrical theory of diffraction (EUTD) are developed for evaluating high frequency electromagnetic (EM) fields within transition regions associated with a two and three dimensional smooth caustic of reflected rays and a composite shadow boundary formed by the caustic termination or the confluence of the caustic with the reflection shadow boundary (RSB). The UGO is a uniform version of the classic geometrical optics (GO). It retains the simple ray optical expressions of classic GO and employs a new set of uniform reflection coefficients. The UGO also includes a uniform version of the complex GO ray field that exists on the dark side of the smooth caustic. The EUTD is an extension of the classic uniform geometrical theory of diffraction (UTD) and accounts for the non-ray optical behavior of the UGO reflected field near caustics by using a two-variable transition function in the expressions for the edge diffraction coefficients. It also uniformly recovers the classic UTD behavior of the edge diffracted field outside the composite shadow boundary transition region. The approach employed for constructing the UGO/EUTD solution is based on a spatial domain physical optics (PO) radiation integral representation for the fields which is then reduced using uniform asymptotic procedures. The UGO/EUTD analysis is also employed to investigate the far-zone RCS problem of plane wave scattering from two and three dimensional polynomial defined surfaces, and uniform reflection, zero-curvature, and edge diffraction coefficients are derived. Numerical results for the scattering and diffraction from cubic and fourth order polynomial strips are also shown and the UGO/EUTD solution is validated by comparison to an independent moment method (MM) solution. The UGO/EUTD solution is also compared with the classic GO/UTD solution. The failure of the classic techniques near caustics and composite shadow boundaries is clearly demonstrated and it is shown that the UGO/EUTD results remain valid and uniformly reduce to the classic results away from the transition regions. Mathematical details on the asymptotic properties and efficient numerical evaluation of the canonical functions involved in the UGO/EUTD expressions are also provided.

Constantinides, E. D.↗

DOE OSTI2025

Shock tube simulations for the three-layer Richtmyer–Meshkov instability with single-mode and multimode perturbations

While the canonical two-component, single-mode Richtmyer–Meshkov instability (RMI) has been extensively studied, relatively less work has focused on the effects of an additional intermediate-density middle layer. This work investigates such three-material RMI configurations at two Atwood number scenarios using the ARES hydrodynamics code. After validation against previous experimental and computational studies, setups corresponding to recent three-layer shock tube experiments are simulated. Cases with both single-mode and multimode perturbations are studied to quantify mixing across the interface between the materials with highest and intermediate density. In particular, this work is able to comprehensibly examine differences between two- and three-dimensional setups for the single-mode and multimode problems. Observations from previous two-layer investigations still apply in the three-layer setup, but over the time horizons considered, there appears to be insufficient nonlinear mode coupling to create significant differences between two- and three-dimensional simulations following the first passage of a shock. Finally, additional reshock simulations have additional nonlinear growth that does result in expected differences between two- and three-dimensional cases in this three-layer setup, but significant differences do not manifest during the time horizon studied.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

Cosmohedra

It has been a long-standing challenge to find a geometric object underlying the cosmological wavefunction for Tr(ϕ 3 ) theory, generalizing associahedra and surfacehedra for scattering amplitudes. In this note, we describe a new class of polytopes — “cosmohedra” — that provide a natural solution to this problem. The faces of associahedra capture the combinatorics of non-overlapping chords of the momentum polygon, reflecting all partial factorizations of amplitudes. Cosmohedra are far richer — instead of non-overlapping chords, their faces capture the “russian doll” structure of non-overlapping subpolygons that determine the wavefunction. We show that cosmohedra are intimately related to associahedra, obtained by “blowing up” faces of the associahedron in a simple way. We give a full combinatorial description of cosmohedron faces and their factorization properties, and provide an explicit realization in terms of facet inequalities that further “shave” the facet inequalities of the associahedron. We also discuss a novel way for computing the wavefunction from cosmohedron geometry that extends the usual connection with polytope canonical forms. We illustrate cosmohedra with examples at tree-level and one loop; the close connection to surfacehedra suggests the generalization to all loop orders. Moving beyond the wavefunction, we briefly describe “cosmological correlahedra” for full correlators, which are one higher-dimensional polytopes, interpolating between associahedra and cosmohedra on opposite facets in an extra direction associated with the total energy. We speculate on how the existence of cosmohedra might suggest a “stringy” formulation for the cosmological wavefunction/correlators, generalizing the way in which the Minkowski sum decomposition of associahedra naturally extend particle to string amplitudes.

Scattering Amplitudes↗

NASA NTRS2008

Learning from Space Entrepreneurs

The early days of rocketry and space exploration in the United States were marked by incredibly rapid progress: a seemingly endless parade of firsts. Not coincidentally, this period also saw more than its fair share of failure, especially in the infamous "kaputnik" days prior to the successful launch of Explorer. Without a standard canon of known quantities to turn to, the early pioneers of rocketry and space flight were forced to dream up new ideas that ranged from the elegant to the bizarre and to accept the fact that the price of radical progress is occasional failure. Nowadays, rapid prototyping and testing have slowed, as we rely more and more on the extensive knowledge pined by our predecessors and on the embarrassment of riches modern engineers get from computational modeling and computer assisted design. In many cases, this leads to much improved or phenomenally more efficient designs. It also, however, fosters a culture so terrified of failure that we over-engineer and overanalyze everything, often tweaking designs for decades before a new system takes flight. (This is not a problem unique to rockets; the same phenomenon seems to have occurred in high-performance jets.) This is one reason why it was possible for President Kennedy to dream of the completion of the Mercury and Gemini missions and a successful landing on the moon in under a decade, while returning to the moon may take nearly twice as long. Lacking access to the tremendous computational resources of the national space program-and, just as importantly, removed from the harsh judgment of public shareholders or congressional appropriations committees-the hungry entrepreneurs who compete for our prizes tend not to display such fear of failure. Instead, most of them follow a rapid "build, test, fly" program. They are willing to throw a handful of concepts against the wall and see what sticks. They often go from drawing on the back of a napkin to firing engines or even flying vehicles in a matter of weeks or months, learning valuable lessons along the way. Indeed, our teams have repeatedly learned many of the most valuable lessons after only a few moments of working with real hardware-lessons that could never have been learned from a CAD drawing, like finding the failure modes of different welding practices or tracking down the interference between an onboard camera and a GPS unit. As Paul Breed, the leader of a Northrop Grumman Lunar Lander Challenge team (playfully called Unreasonable Rocket), is fond of saying, "In computer simulations the plumbing never leaks. In real life, it always does."

Pomerantz, William↗

NASA NTRS1975

On the tesseral-harmonics resonance problem in artificial-satellite theory

The longitude-dependent part of the geopotential usually gives rise only to short-period effects in the motion of an artificial satellite. However, when the motion of the satellite is commensurable with that of the earth, the path of the satellite repeats itself relative to the earth and perturbations build up at each passage of the satellite in the same spot, so that there can be important long-period effects. In order to take these effects into account in deriving a theoretical solution to the equations of motion of an artificial satellite, it is necessary to select terms in the longitude-dependent part of the geopotential that will contribute significantly to the perturbations. Attempts made to obtain a selection that is valid in a general case, regardless of the initial eccentricity of the orbit and of the order of the resonance, are reported. The solution to the equations of motion of an artificial satellite, in a geopotential thus determined, is then derived by using Hori's method by Lie series, which, by its properties regarding canonical invariance, has proved advantageous in the classical theory.

Romanowicz, B. A.↗