Search NASA⌕ Search

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.

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

At least 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

Experiments on Diffusion Flame Structure of a Laminar Vortex Ring

The study of flame-vortex interactions provides one of the means to better understand turbulent combustion, and allows for canonical configurations that contain the fundamental elements found in turbulent flames, These include concentrated vorticity, entrainment and mixing, strain and nonequilibrium phenomena, diffusion and differential diffusion, partial premixing and diluent effects, and heat release effects. In flame- vortex configurations, these fundamental elements can be studied under more controlled conditions than is possible in direct investigations of turbulent flames. Since the paper of Marble, the problem of the flame-vortex interaction has received considerable attention theoretically, numerically and experimentally. Several configurations exist for study of the premixed flame/vortex ring interaction but more limited results have been obtained to date for the diffusion flame/vortex ring case. The setup of Chen and Dahm, which is conceptually similar to that of Karagozian and Manda and Karagozian, Suganuma and Strom where the ring is composed of fuel and air and combustion begins during the ring formation process, is used in the current study. However, it is essential to conduct the experiments in microgravity to remove the asymmetries caused by buoyancy and thus obtain highly symmetric and repeatable interactions. In previous studies it was found that the flame structure of the vortex ring was similar to that obtained analytically by Karagozian and Manda. Dilution of propane with nitrogen led mainly to a reduction in flame luminosities, flame burnout times were affected by both fuel volumes and amount of dilution, and a simple model of the burnout times was developed. In this paper, a discussion on reacting ring displacement and flame burnout time will be given, and the flame structures of vortex rings containing ethane and air will be compared to those of propane reacting in air.

Chen, Shin-Juh↗
DOE OSTI2026

Hierarchical-embedding autoencoder with a predictor as efficient architecture for learning time-evolution in multi-scale turbulent flows

We introduce a scale-aware, data-driven deep learning modeling framework for accurately predicting the time evolution of multi-scale turbulent plasma and liquid flows. The approach is motivated by the idea of scale separation. Structures of vastly different length scales emerge in these systems, and interactions between these structures occur only locally. To exploit this structure, the flow state is transformed by a hierarchical, fully convolutional autoencoder, not into a single embedding layer as in conventional convolutional surrogate models, but into a series of embedding layers. A stepwise training strategy ensures that fine-scale features are encoded on a high-resolution grid, while larger structures are represented on progressively coarser layers. The time evolution predictor advances all embedding layers in sync, capturing local interactions between features at the same scale as well as between all scales. This approach enables efficient modeling of multi-scale systems since negligible interactions between distant, small-scale structures do not need to be directly modeled. Our hierarchical-embedding autoencoder with a predictor framework is evaluated on canonical examples of multi-scale turbulence: two-dimensional Kolmogorov flow and Hasegawa–Wakatani plasma turbulence. In both cases, the proposed framework significantly improves predictive accuracy relative to conventional convolutional network architectures. A significant improvement in prediction accuracy was observed for crucial statistical characteristics of the Hasegawa–Wakatani plasma as well as for individual trajectories of the Kolmogorov flow turbulence. Importantly, the model's rollout for the Hasegawa–Wakatani problem demonstrates a four-order-of-magnitude speedup compared to traditional numerical solvers.

Khrabry, Alexander I. [Princeton Univ., NJ (United↗
NASA NTRS1979

Linear state feedback, quadratic weights, and closed loop eigenstructures

Results are given on the relationships between closed loop eigenstructures, state feedback gain matrices of the linear state feedback problem, and quadratic weights of the linear quadratic regulator. Equations are derived for the angles of general multivariable root loci and linear quadratic optimal root loci, including angles of departure and approach. The generalized eigenvalue problem is used for the first time to compute angles of approach. Equations are also derived to find the sensitivity of closed loop eigenvalues and the directional derivatives of closed loop eigenvectors (with respect to a scalar multiplying the feedback gain matrix or the quadratic control weight). An equivalence class of quadratic weights that produce the same asymptotic eigenstructure is defined, sufficient conditions to be in it are given, a canonical element is defined, and an algorithm to find it is given. The behavior of the optimal root locus in the nonasymptotic region is shown to be different for quadratic weights with the same asymptotic properties.

Thompson, P. M.↗
NASA NTRS1984

Equivalence of the generalized Lie-Hori method and the method of averaging

In this investigation, a comparison is made of two methods for developing perturbation theories for non-canonical dynamical systems. The methods compared are the generalized Lie-Hori method and the method of averaging. In the comparison presented here, the equivalence of the methods up to the second order in the small parameter is shown. However, the approach used can be extended to demonstrate the equivalence for higher orders. To illustrate the equivalence Duffing's equation, the van der Pol equation and the oscillator with quadratic damping problem are solved using each method.

Ahmed, A. H.↗
NASA NTRS1972

Control system optimization studies. Volume 2: High frequency cutoff filter analysis

The problem of digital implementation of a cutoff filter is approached with consideration to word length, sampling rate, accuracy requirements, computing time and hardware restrictions. Computing time and hardware requirements for four possible programming forms for the linear portions of the filter are determined. Upper bounds for the steady state system output error due to quantization for digital control systems containing a digital network programmed both in the direct form and in the canonical form are derived. This is accomplished by defining a set of error equations in the z domain and then applying the final value theorem to the solution. Quantization error was found to depend upon the digital word length, sampling rate, and system time constants. The error bound developed may be used to estimate the digital word length and sampling rate required to achieve a given system specification. From the quantization error accumulation, computing time and hardware point of view, and the fact that complex poles and zeros must be realized, the canonical form of programming seems preferable.

Fong, M. H.↗
DOE OSTI2024

SEAS Communication Engine: An Extensible, Flexible Wrapper for Co-Simulation Agents

When modeling and analyzing the power grid and other large scale systems, researchers often express scenarios as optimization problems and feed them into advanced software solvers. In order to allow multiple solvers to communicate with each other and share data from different domains, the National Renewable Energy Laboratory (NREL) and associated Department of Energy (DOE) labs have developed a software framework called the Hierarchical Engine for Large-scale Infrastructure Co-Simulation (HELICS). HELICS allows cosimulation via a collection of client libraries for different languages that can be called from the appropriate optimization software. However, these client libraries do not provide a higher level of abstraction beyond reading and writing data off of the shared HELICS bus. In this paper, we describe a new software library called the SEAS Communication Engine that exposes a higher-level API for running cosimulation problems. The SEAS Engine provides a class-based abstraction on top of the Python HELICS client, in order to allow users to implement their domain-specific cosimulations without needing to interact with core HELICS primitives. This will make adoption of HELICS and cosimulation in general easier, by exposing a simpler API. In the second part of the paper, we validate our library on a collection of different simulation examples, including the canonical IEEE 13 Bus Feeder. Lastly, we demonstrate using the SEAS Engine to directly call domain-specific code written in the Julia programming language. Our hope is that this will serve as a template for easily calling software in different programming languages via the SEAS Engine, thereby avoiding code duplication and complexity.

co-simulation↗
NASA NTRS1974

Analytical satellite theories based on a new set of canonical elements

A new analytical satellite theory is presented. Instead of the 6 classical elements of Delaunay, a set of 8 canonical elements is used. Whereas the time is the independent variable in classical theory, the true anomaly is the independent variable in the new theory. The new approach has four features: (1) The amount of formulas in the solution is reduced considerably. (2) The first order results are almost as accurate as second order results in classical theory. (3) The theory is easier to understand from a didactical point of view. (4) The problems connected with the inaccuracy of the mean motion that are typical for classical satellite theory are no longer present. The new elements are applied to analytical solutions of the zonal oblateness problem and to the problem of the 24 hour satellite.

Scheifele, G.↗
NASA NTRS1976

The elimination of short and intermediate period terms from the problem of a high altitude earth satellite

The Delaunay-Similar elements of Scheifele (1976) are applied to the problem of an earth satellite that is perturbed by the sun, moon, and J2. All three effects are assumed to be the same order of magnitude. Since the external body terms depend explicitly on time, the time element appears as an additional angle variable. Also, the eccentric anomaly is used as a noncanonical auxiliary variable. A solution to the first Von Zeipel equation allows simultaneous elimination of short- and intermediate-period terms. The canonical transformation to mean elements is defined by a generating function that is a series involving Bessel coefficients.

Graf, O. F., Jr.↗
NASA NTRS1990

Mixed variational formulation of finite element analysis of acoustoelastic/slosh fluid-structure interaction

A general three-field variational principle is obtained for the motion of an acoustic fluid enclosed in a rigid or flexible container by the method of canonical decomposition applied to a modified form of the wave equation in the displacement potential. The general principle is specialized to a mixed two-field principle that contains the fluid displacement potential and pressure as independent fields. This principle contains a free parameter alpha. Semidiscrete finite-element equations of motion based on this principle are displayed and applied to the transient response and free-vibrations of the coupled fluid-structure problem. It is shown that a particular setting of alpha yields a rich set of formulations that can be customized to fit physical and computational requirements. The variational principle is then extended to handle slosh motions in a uniform gravity field, and used to derived semidiscrete equations of motion that account for such effects.

Felippa, C. A.↗
NASA NTRS1991

Mixed variational formulations of finite element analysis of elastoacoustic/slosh fluid-structure interaction

A general three-field variational principle is obtained for the motion of an acoustic fluid enclosed in a rigid or flexible container by the method of canonical decomposition applied to a modified form of the wave equation in the displacement potential. The general principle is specialized to a mixed two-field principle that contains the fluid displacement potential and pressure as independent fields. This principle contains a free parameter alpha. Semidiscrete finite-element equations of motion based on this principle are displayed and applied to the transient response and free-vibrations of the coupled fluid-structure problem. It is shown that a particular setting of alpha yields a rich set of formulations that can be customized to fit physical and computational requirements. The variational principle is then extended to handle slosh motions in a uniform gravity field, and used to derive semidiscrete equations of motion that account for such effects.

Felippa, Carlos A.↗
NASA NTRS1993

Tsien's method for generating non-Keplerian trajectories. Part 2: The question of thrust to orbit a sphere and the restricted three-body problem

Tsien's method is extended to treat the orbital motion of a body undergoing accelerations and decelerations. A generalized solution is discussed for the generalized case where a body undergoes azimuthal and radial thrust and the problem is further simplified for azimuthal thrust alone. Judicious selection of thrust could generate either an elliptic or hyperbolic trajectory. This is unexpected especially when the body has only enough energy for a lower state trajectory. The methodology is extended treating the problem of vehicle thrust for orbiting a sphere and vehicle thrust within the classical restricted three-body problem. Results for the latter situation can produce hyperbolic trajectories through eigen value decomposition. Since eigen values for no-thrust can be imaginary, thrust can generate real eigen values to describe hyperbolic trajectories. Keplerian dynamics appears to represent but a small subset of a much larger non-Keplerian domain especially when thrust effects are considered. The need for high thrust long duration space-based propulsion systems for changing a trajectory's canonical form is clearly demonstrated.

Murad, P. A.↗
DOE OSTI2025

A Comparison of GPU-Accelerated Multiphase CFD Solvers on the Polaris Supercomputer: Part 1

This report is in support of the Innovative and Novel Computational Impact on Theory and Experiment (INCITE) program sponsored by the U.S. Department of Energy (USDOE). With INCITE-level resources, one project, titled BubblyFlow, was granted computational resources for the 2025 calendar year on the Polaris supercomputer at the Argonne Leadership Computing Facility (ALCF). The project aims to conduct simulations to understand the fundamental characteristics of turbulent bubbly flow phenomena in nature. Staff at the ALCF and Argonne’s Computational Science division, along with collaborators at the City College of New York and University of Illinois at Chicago, helped a summer student to assess the accuracy and performance of two high performance computing (HPC) codes. Both codes, ImExLBM and FluTAS, are fundamentally different in their mathematical and numerical modeling. However, both may be used to solve the same physical problem. The collaboration sought to better understand the differences between both codes in terms of accuracy and efficiency. This would ultimately help the BubblyFlow project better utilize resources and establish a knowledge-base of code capabilities in future simulation campaigns. We compare ImExLBM and FluTAS, two high-performance multiphase computational fluid dynamics (CFD) solvers, in terms of physical fidelity, time-to-solution, and parallel efficiency. We validate ImExLBM (Implicit-Explicit Lattice Boltzmann Method) against a canonical benchmark and assess it’s performance relative to FluTAS (Fluid Transport Accelerated Solver), a well-established open-source CFD code.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

A Scaling Study for Incompressible Multispecies Solver in Vertex-CFD

Multispecies incompressible flows occur widely in engineering and environmental applications, such as chemical reactors, fuel cells, ocean mixing, and biomedical systems. However, accurately resolving the complex transport and mixing phenomena associated with multiple interacting species remains computationally challenging, especially for large-scale problems. In this study, we present a robust, high-performance computing--enabled multispecies incompressible Navier–Stokes solver integrated within the Vertex-CFD framework. Our solver employs a fully coupled, implicit, finite element--based formulation that accurately captures the advection, diffusion, and interaction of multiple species in incompressible flows by leveraging the Kokkos library for parallel computing to achieve high computational efficiency. For pressure coupling, the entropically damped artificial compressibility method is utilized. We validated the solver against canonical test cases, including multispecies advection, diffusion, and Bateman systems; the results demonstrate second- and third-order spatial accuracy and consistent convergence. Additionally, we demonstrated the strong and weak scaling study results obtained on the leadership-class high-performance computing system, Frontier at Oak Ridge National Laboratory.

Oz, Furkan [ORNL] (ORCID:0000000265831724)↗
NASA NTRS1984

X-rays from accretion of red giant winds

X-ray observations of the late-type red giants Mira and R Aqr obtained with the Einstein Observatory are presented, and the general problems of white dwarf accretion from late-type giant winds is considered. The extremely low measured luminosities obtained for the two systems leads to the conclusion that the companions of Mira and R Aqr are most likely low-mass main sequence objects rather than white dwarfs as is usually assumed. The expected X-ray luminosities of true red giant/white dwarf systems are considered, and it is concluded that far too few have been detected if the canonical accretion scenario is adopted. A possible explanation of this situation in terms of grain-dominated Eddington-limited accretion is proposed.

Jura, M.↗
NASA NTRS1980

A note on a Lagrangian formulation for motion about the collinear points

A Lagrangian formulation for the three-dimensional motion of a satellite in the vicinity of the collinear points of the circular-restricted problem is reconsidered. It is shown that the influence of the primaries can be expressed in the form of two third-body disturbing functions. By use of this approach, the equations for the Lagrangian and for the motion itself are readily developed into highly compact expressions. All orders of the non-linear developments are shown to be easily obtainable using well-known recursive relationships. The resulting forms for these equations are well suited for use in the initial phase of canonical or non-canonical investigations.

Richardson, D. L.↗
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.↗