On the solution of error covariance difference equations by means of canonical decomposition and z transform.
Canonical decomposition of nonlinear error covariance difference equation derived for discrete estimation problems
NONLINEAR EQUATION↗
SEARCH · Engineering Papers
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.
Photometric theory for wide-angle phenomena
An examination is made of the problem posed by wide-angle photographic photometry, in order to extract a photometric-morphological history of Comet P/Halley. Photometric solutions are presently achieved over wide angles through a generalization of an assumption-free moment-sum method. Standard stars in the field allow a complete solution to be obtained for extinction, sky brightness, and the characteristic curve. After formulating Newton's method for the solution of the general nonlinear least-square problem, an implementation is undertaken for a canonical data set. Attention is given to the problem of random and systematic photometric errors.
Usher, Peter D.↗
A robust spectral element implementation of the $k - τ$ RANS model in Nek5000/NekRS
The $k - ω$ Reynolds Averaged Navier Stokes (RANS) model is one of the industry standard approaches for modeling of turbulent flows. It performs better than the $k - ϵ$ model for low Reynolds number flows and is also more suitable for boundary layers with adverse pressure gradients. Major drawback of the model, however, is that the asymptotic value of $ω$ at the walls is singular, necessitating the use of a contrived “sufficiently” large value for $ω$ as the boundary condition for its transport equation. Here, this invariably leads to the solution being sensitive to near wall grid spacing. While an acceptable solution for low order (finite volume) methods, the excessive near wall gradients lead to persistent numerical stability issues in high order codes. To alleviate the problem, specifically in the context of the high order spectral element code Nek5000, a regularized $k - ω$ approach was formulated in our prior work (Tomboulides et al., 2018). The formulation, however, relies on the use of wall distance and its gradients for modeling the closure terms and can pose problems for simulations in complex geometries. This work presents a novel implementation of the $k - τ$ RANS model in Nek5000, where $τ = 1/ω$, eliminating the need for regularization, owing to the asymptotically bounded behavior of the source terms in the $τ$ transport equation, and also eliminating dependence on wall distance. Robustness and stability of the $k - τ$ model is ensured through implicit treatment of the source terms and their careful numerical implementation and demonstrated through several cases aimed at verification and validation. Studies include both canonical and engineering relevant problems, viz., turbulent channel flow, pipe flow, backward facing step, flow over NACA 0012 airfoil and flow in a T-junction. Results from the $k - τ$ model are shown to be consistent with regularized $k - ω$ model and also with the $k - ω$ SST model in OpenFOAM (for select studies). Comparison with experimental data is also shown, where available, to bolster validation efforts for the $k - τ$ model implementation through prediction of key turbulent quantities of interest.
Nek5000↗
A conservative discontinuous Galerkin algorithm for particle kinetics on smooth manifolds
A novel, conservative discontinuous Galerkin algorithm is presented for particle kinetics on manifolds. The motion of particles on the manifold is represented using both canonical and non-canonical Hamiltonian formulations. Our schemes apply to both formulations, but the canonical formulation results in a particularly efficient scheme that also conserves particle density and energy exactly. The collisionless update is coupled to a Bhatnagar-Gross-Krook (BGK) collision operator that provides a simplified model for relaxation to local thermodynamic equilibrium. An iterative scheme is constructed to ensure collisional invariants (density, momentum and energy) are preserved numerically. Rotation of the manifold is incorporated by modifying the Hamiltonian while ensuring a canonical formulation. Several test problems, including a kinetic version of the classical Sod shock problem, Kelvin-Helmholtz instability on the surfaces of a sphere and a hyperboloid, with and without rotations, are presented. A prospectus for further development of this approach to simulation of kinetic theory in general relativity is presented.
Discontinuous Galerkin↗
Numerical studies of the margin of vortices with decaying cores
The merging of vortices to a single one is a canonical incompressible viscous flow problem. The merging process begins when the core sizes or the vortices are comparable to their distances and ends when the contour lines of constant vorticity lines are circularized around one center. Approximate solutions to this problem are constructed by adapting the asymptotic solutions for distinct vortices. For the early stage of merging, the next-order terms in the asymptotic solutions are added to the leading term. For the later stage of merging, the vorticity distribution is reinitialized by vortices with overlapping core structures guided by the 'rule of merging' and the velocity of the 'vortex centers' are then defined by a minimum principle. To show the accuracy of the approximate solution, it is compared with the finite-difference solution.
Liu, G. C.↗
A 3-D Nodal-Averaged Gradient Approach for Unstructured-Grid Cell-Centered Finite-Volume Methods for Application to Turbulent Hypersonic Flow
A 2-D nodal weighted least-squares gradient method and a related face-averaged nodal gradient approach that were developed for use with triangular grids are extended to 3-D for use with tetrahedral grids. In addition, a method, developed in 2-D, to stabilize the iterative convergence of these methods on quadrilateral cells is described and extended to 3-D and remedies are investigated to determine the nodal gradient averaging approach most suitable for use with grids made up of hexahedral, prismatic, pyramidal and tetrahedral cells. Moreover, due to an interest in hypersonic flow, a robust multidimensional gradient limiter procedure that is consistent with the stencil used to construct the nodal gradients is described. Finally, we demonstrate that the resulting 3-D methods are sufficiently robust for use in scramjet computations through the solution of three canonical turbulent hypersonic flow problems as well as a physically realistic 3-D scramjet inlet geometry.
Jeffery A White↗
Canonical solutions for unsteady flow fields
The initial value problem of one-dimensional gas-dynamics involving discontinuous, nonuniform initial data is discussed. Canonical solutions which are valid in a small x, t region aroung a discontinuity, and which include the first order effects of nonuniformities in the data, are derived explicitly. The theory is derived by considering a group of elementary piston problems. Solutions with a shock or with a centered expansion wave are worked out individually in order to relate initial flow properties and their gradients to the speed and acceleration of the discontinuity waves. They are then combined to represent the solution of a general initial value problem by regarding the piston path as a contact line. In addition, problems with chemical reaction are discussed in terms of elementary piston problems which involve strong detonation waves, Chapman-Jouguet detonation waves, and deflagration waves.
Liu, G. C.↗
Canonical forms of multidimensional steady inviscid flows
Canonical forms and canonical variables for inviscid flow problems are derived. In these forms the components of the system governed by different types of operators (elliptic and hyperbolic) are separated. Both the incompressible and compressible cases are analyzed, and their similarities and differences are discussed. The canonical forms obtained are block upper triangular operator form in which the elliptic and non-elliptic parts reside in different blocks. The full nonlinear equations are treated without using any linearization process. This form enables a better analysis of the equations as well as better numerical treatment. These forms are the analog of the decomposition of the one dimensional Euler equations into characteristic directions and Riemann invariants.
Taasan, Shlomo↗
Passive stabilization for large space systems
The optimal tuning of multiple tuned-mass dampers for the transient vibration damping of large space structures is investigated. A multidisciplinary approach is used. Structural dynamic techniques are applied to gain physical insight into absorber/structure interaction and to optimize specific cases. Modern control theory and parameter optimization techniques are applied to the general optimization problem. A design procedure for multi-absorber multi-DOF vibration damping problems is presented. Classical dynamic models are extended to investigate the effects of absorber placement, existing structural damping, and absorber cross-coupling on the optimal design synthesis. The control design process for the general optimization problem is formulated as a linear output feedback control problem via the development of a feedback control canonical form. The techniques are applied to sample micro-g and pointing problems on the NASA dual keel space station.
Sesak, J. R.↗
LuGo: An enhanced quantum phase estimation implementation
Quantum Phase Estimation (QPE) is a cardinal algorithm in quantum computing that plays a crucial role in various applications, including cryptography, molecular simulation, and solving systems of linear equations. However, the standard implementation of QPE faces challenges related to time complexity and circuit depth, which limit its practicality for large-scale computations. We introduce LuGo, a novel framework designed to enhance the performance of QPE by reducing circuit duplication, as well as using parallelization techniques to achieve faster generation of the QPE circuit and gate reduction. We validate the effectiveness of our framework by generating quantum linear solver circuits, which require both QPE and inverse QPE, to solve linear systems of equations. LuGo achieves significant improvements in both computational efficiency and hardware requirements without compromising on accuracy. Compared to a standard QPE implementation, LuGo reduces time consumption to generate a circuit that solves a 2 6 × 2 6 system matrix by a factor of 50.68 and over 31× reduction of quantum gates and circuit depth, with no fidelity loss on an ideal quantum simulator. Furthermore, we demonstrated the versatility and scalability of LuGo enabled HHL algorithm by simulating a canonical Hele-Shaw fluid problem using a quantum simulator. With these advantages, LuGo paves the way for more efficient implementations of QPE, enabling broader applications across several quantum computing domains.
Quantum algorithm↗
Linear and regular celestial mechanics- Perturbed two-body motion, numerical methods, canonical theory
Book on linear and ordinary celestial mechanics covering perturbed two body motion, numerical methods, canonical theory and initial value problems
Scheifele, G.↗
High Order Numerical Methods for LES of Turbulent Flows with Shocks
Simulation of turbulent flows with shocks employing explicit subgrid-scale (SGS) filtering may encounter a loss of accuracy in the vicinity of a shock. In this work we perform a comparative study of different approaches to reduce this loss of accuracy within the framework of the dynamic Germano SGS model. One of the possible approaches is to apply Harten's subcell resolution procedure to locate and sharpen the shock, and to use a one-sided test filter at the grid points adjacent to the exact shock location. The other considered approach is local disabling of the SGS terms in the vicinity of the shock location. In this study we use a canonical shock-turbulence interaction problem for comparison of the considered modifications of the SGS filtering procedure. For the considered test case both approaches show a similar improvement in the accuracy near the shock.
High Order↗
High Order Numerical Methods for LES of Turbulent Flows with Shocks
Simulation of turbulent flows with shocks employing explicit subgrid-scale (SGS) filtering may encounter a loss of accuracy in the vicinity of a shock. In this work we perform a comparative study of different approaches to reduce this loss of accuracy within the framework of the dynamic Germano SGS model. One of the possible approaches is to apply Harten’s subcell resolution procedure to locate and sharpen the shock, and to use a one-sided test filter at the grid points adjacent to the exact shock location. The other considered approach is local disabling of the SGS terms in the vicinity of the shock location. In this study we use a canonical shock-turbulence interaction problem for comparison of the considered modifications of the SGS filtering procedure. For the considered test case both approaches show a similar improvement in the accuracy near the shock.
LES↗
NASA Taxonomies for Searching Problem Reports and FMEAs
Many types of hazard and risk analyses are used during the life cycle of complex systems, including Failure Modes and Effects Analysis (FMEA), Hazard Analysis, Fault Tree and Event Tree Analysis, Probabilistic Risk Assessment, Reliability Analysis and analysis of Problem Reporting and Corrective Action (PRACA) databases. The success of these methods depends on the availability of input data and the analysts knowledge. Standard nomenclature can increase the reusability of hazard, risk and problem data. When nomenclature in the source texts is not standard, taxonomies with mapping words (sets of rough synonyms) can be combined with semantic search to identify items and tag them with metadata based on a rich standard nomenclature. Semantic search uses word meanings in the context of parsed phrases to find matches. The NASA taxonomies provide the word meanings. Spacecraft taxonomies and ontologies (generalization hierarchies with attributes and relationships, based on terms meanings) are being developed for types of subsystems, functions, entities, hazards and failures. The ontologies are broad and general, covering hardware, software and human systems. Semantic search of Space Station texts was used to validate and extend the taxonomies. The taxonomies have also been used to extract system connectivity (interaction) models and functions from requirements text. Now the Reconciler semantic search tool and the taxonomies are being applied to improve search in the Space Shuttle PRACA database, to discover recurring patterns of failure. Usual methods of string search and keyword search fall short because the entries are terse and have numerous shortcuts (irregular abbreviations, nonstandard acronyms, cryptic codes) and modifier words cannot be used in sentence context to refine the search. The limited and fixed FMEA categories associated with the entries do not make the fine distinctions needed in the search. The approach assigns PRACA report titles to problem classes in the taxonomy. Each ontology class includes mapping words - near-synonyms naming different manifestations of that problem class. The mapping words for Problems, Entities and Functions are converted to a canonical form plus any of a small set of modifier words (e.g. non-uniformity NOT + UNIFORM.) The report titles are parsed as sentences if possible, or treated as a flat sequence of word tokens if parsing fails. When canonical forms in the title match mapping words, the PRACA entry is associated with the corresponding Problem, Entity or Function in the ontology. The user can search for types of failures associated with types of equipment, clustering by type of problem (e.g., all bearings found with problems of being uneven: rough, irregular, gritty ). The results could also be used for tagging PRACA report entries with rich metadata. This approach could also be applied to searching and tagging failure modes, failure effects and mitigations in FMEAs. In the pilot work, parsing 52K+ truncated titles (the test cases that were available), has resulted in identification of both a type of equipment and type of problem in about 75% of the cases. The results are displayed in a manner analogous to Google search results. The effort has also led to the enrichment of the taxonomy, adding some new categories and many new mapping words. Further work would make enhancements that have been identified for improving the clustering and further reducing the false alarm rate. (In searching for recurring problems, good clustering is more important than reducing false alarms). Searching complete PRACA reports should lead to immediate improvement.
Malin, Jane T.↗
Riccati group invariants of linear hamiltonian systems
The action of the Riccati group on the Riccati differential equation is associated with the action of a subgroup of the symplectic group on a set of hamiltonian matrices. Within this framework various sets of canonical forms are developed for the matrix coefficients of the Riccati differential equation. The canonical forms presented are valid for arbitrary Kronecker indices, and it is shown that the Kronecker indices are invariants for this group action. These canonical forms are useful for studying problems arising in the areas of optimal decentralized control and the spectral theory of optimal control problems.
Garzia, M. R.↗
Canonical-variables multigrid method for steady-state Euler equations
A novel approach for the solution of inviscid flow problems for subsonic compressible flows is described. The approach is based on canonical forms of the equations, in which subsystems governed by hyperbolic operators are separated from those governed by elliptic ones. The discretizations used as well as the iterative techniques for the different subsystems are inherently different. Hyperbolic parts, which describe, in general, propagation phenomena, are discretized using upwind schemes and are solved by marching techniques. Elliptic parts, which are directionally unbiased, are discretized using h-elliptic central discretizations, and are solved by pointwise relaxations together with coarse grid acceleration. The resulting discretization schemes introduce artificial viscosity only for the hyperbolic parts of the system; thus a smaller total artificial viscosity is used, while the multigrid solvers used are much more efficient. Solutions of the subsonic compressible Euler equations are achieved at the same efficiency as the full potential equation.
Taasan, Shlomo↗
Intrastep, Stage-Value Predictors for Diagonally-Implicit Runge-Kutta Methods
To better identify the necessary attributes of good stage-value predictors (SVPs), numerous SVPs are designed for an existing: ESDIRK4(3)7L[2]SA [26] and a new: ESDIRK4(3)8L[2]SA scheme.1 Both are stifflyaccurate, stage-order two, explicit, singly-diagonally implicit Runge–Kutta (ESDIRK) schemes. Tradeoffs are studied in the parameter spaces enforcing the constraints on accuracy, linear stability, nonlinear stability and coefficient size to determine which objectives correlate with effective predictors. The SVPs are tested on three challenging external aerodynamics problems [107 − 108 degrees of freedom (DoFs)], each with a different level of stiffness. The problems include two 3D airfoils simulations and one canonical turbulence simulation. All simulations use the compressible Navier-Stokes equations (CNSE). An entropy stable spectral collocation formulation is used for discretizing the spatial terms in the equations. Simulations are performed at a wide variety of temporal error tolerances. Problems that are sufficiently stiff (e.g., lax temporal error tolerances) benefit from SVPs designed with second-order accuracy and stability properties: A-stability, and L-stability, rather than high accuracy constraints. Simulations with modest stiffness (e.g., strict error tolerances) are better suited for SVPs designed using high accuracy constraints. Designing SVPs with enhanced stability properties is tedious but worthwhile. Simulation times are reduced with optimal SVPs by as much as 100% on some stages, with combined stepwise improvements of between 50 − 100% for both methods. A comparative study is performed with the two aforementioned methods as well as four other ESDIRKs. The newly designed ESDIRK4(3)8L[2]SA with γ ≈ 1/10, proves to be the most efficient of the six tested ESDIRK schemes simulating the CNSE.
Diagonally-Implicit Runge-Kutta↗
The uniform asymptotic swallowtail approximation - Practical methods for oscillating integrals with four coalescing saddle points
Methods that can be used in the numerical implementation of the uniform swallowtail approximation are described. An explicit expression for that approximation is presented to the lowest order, showing that there are three problems which must be overcome in practice before the approximation can be applied to any given problem. It is shown that a recently developed quadrature method can be used for the accurate numerical evaluation of the swallowtail canonical integral and its partial derivatives. Isometric plots of these are presented to illustrate some of their properties. The problem of obtaining the arguments of the swallowtail integral from an analytical function of its argument is considered, describing two methods of solving this problem. The asymptotic evaluation of the butterfly canonical integral is addressed.
Connor, J. N. L.↗