Search NASA⌕ Search

SEARCH · Search NASA

Results for “weak solutions”

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 145 records · Page 8

Bifurcation analysis of a simple analytic model of self-propagating star formation

We investigate the structure and stability of rotationally symmetric nonhomogeneous time-independent solutions derived from a simple analytic model of self-propagating star formation. For this purpose we employ two methodologies: We use bifurcation theoretical methods to prove the existence of nonhomogeneous axisymmetric stationary solutions of an appropriate nonlinear evolution equation for the stellar density. We show that the nonhomogeneous solution branch bifurcates from the homogeneous one at a critical parameter value of the star formation rate. Further, the analytical theory allows us to show that the new solution set is stable in the weakly nonlinear regime near the bifurcation point. To follow the solution branch further, we use numerical methods. The numerical calculation shows the structure and stability of these solutions. We conclude that no periodic time-dependent solutions of this special model exist, and no further bifurcations can be found. The same results have been found in simulations of stochastic self-propagating star formation based on similar models. Therefore, our findings provide a natural explanation, why long-lived large-scale structure have not been found in those simulations.

Neukirch, Thomas↗

Investigating Low-Altitude Constellations of Ad-Hoc Lunar PNT System for Distributed Spacecraft Autonomy

In this study, we examine a low-altitude Lunar Position, Navigation, and Timing (LPNT) constellations and the localization performance of Centralized Extended Kalman Filter (CEKF) and Decentralized Extended Kalman Filter (DEKF) algorithms. The primary investigation involves a 100-node swarm operating at a 100 km altitude, in contrast to previous studies that examined a 21-node asset in a frozen-orbit at 5,500 km. The autonomous operation of large-scale swarm is based on two-way Inter-Satellite Link (ISL) measurements, which involve pseudoranges and relative velocities among swarm nodes. We perform a numerical assessment of the two filtering approaches, utilizing ‘fully sampled’ measurements from all available assets as well as ‘two ISL’ measurements where each spacecraft is restricted to only two antennas. This research includes an analysis of CEKF under 2-ISL constraints and evaluates the performance of DEKF in a 100-node swarm, which has not been explored in previous studies. In addition, we examine the impact of increasing the sampling frequency for DEKF, showing that the update cycle can be shortened from a 10-minute interval. A novel approach for ‘2-ISL limited’ DEKF will also be introduced, using a matching formulation that exhaustively enumerates all potential matches. This study provides valuable insights into large-scale distributed swarm operations, considering various filter configurations, sampling frequencies, matching strategies, and scalability of CEKF and DEKF for low-altitude LPNT applications. The Lunar PNT technology plays a key role in providing reliable and robust navigation services on the Moon's surface and the South pole, where the primary Lunar missions are planned. To support upcoming Lunar missions, including small satellites from NASA's Commercial Lunar Payload Services program, the Lunar PNT system must be adaptable to smaller platforms like CubeSats. Driven by the growing involvement of public and private exploration partnerships, the traditional low Earth orbit missions are shifting to beyond geosynchronous orbit [1]. These upcoming missions aim to foster a sustainable and innovative exploration program, in collaboration with commercial and international partners, to facilitate human expansion throughout the solar system and return new knowledge and opportunities to Earth [2]. As part of this trend, there are increasing efforts to utilize science missions in Lunar orbit to develop a non-dedicated and ad-hoc PNT network system. Two traditional approaches, the Deep Space Network (DSN) and the weak signal Global Positioning System (GPS), are established deep-space navigation technologies for missions beyond the geosynchronous orbit. Beginning in 1958, the DSN was developed to communicate with the Explorer 1 spacecraft based on the use of radiometric tracking in spacecraft navigation [3]. The DSN is capable of providing nearly unfettered coverage to spacecraft beyond low-Earth orbit (LEO), however, increased space mission volume has created concerns about future expectations of DSN usage for spacecraft navigation [4]. For cislunar mission applications, the position accuracy using DSN achieves 100 m (3σ) with at least three geometrically diverse ground stations when using radiometric tracking alone [5]. The DSN's dependence on Earth-based ground stations restricts its operational capabilities to periods of Earth visibility. This limitation, coupled with its poor localization performance, renders the DSN unsuitable for future lunar missions that demand continuous tracking and precise positioning. To satisfy the increasing requirements of DSN in Lunar applications, spacecrafts are also required to improve their onboard antenna power and efficiency of the transmission. However, there is an important aggregate cost trade between adding capabilities to every spacecraft and adding to a capacity on the ground that serves multiple spacecraft [6]. A weak GPS system can provide PNT service while the user spacecraft is bound to the Moon, leveraging a single, steerable high gain antenna with the relatively narrow beam which includes all the sources in its field of view [7]. However, the higher the altitude the receiver is above the GPS constellations, the poorer and the weaker are the relative geometry and the received signal powers, respectively, leading to a significant navigation accuracy reduction [8]. The transmitted power becomes weaker with increasing distance from the Earth as well as signals tracked from one of the side lobes of the GPS antenna pattern. As a results, the number of visible satellites and relative geometric condition of the GPS satellites at very high altitude drops dramatically and reduces the navigation solution accuracy. Therefore, the weak GPS system is also not an ideal way to provide PNT service to upcoming Lunar missions when considering its limited geometric condition and the recued navigation accuracy. Another navigation approach on the Moon is being developed, similar to the Global Navigation Satellite System (GNSS) on Earth, aiming to offer navigation service with continuous 24/7 coverage across the entire Lunar surface. For example, lunar communications relay and navigation systems (LCRNS) by NASA and Lunar navigation satellite systems (LNSS) by JAXA are designed to serve as dedicated Position, Navigation, and Timing (PNT) systems for the Moon. However, designing a dedicated LNSS and PNT service involves additional challenges, which are unique to the lunar environment, including limited payload capacity for the CubeSat platform, i.e., the size, weight, and power (SWaP) of the onboard clock, limited lunar ground monitoring stations, and limited financial investment as compared to the legacy Earth-GPS [9]. NASA’s focus on utilizing CubeSat platforms on the Moon leads to an alternative Lunar navigation platform that leverages the existing Lunar science and exploration assets. The small satellites used in Lunar missions can be used to create a low-cost, autonomous, ad-hoc, and on-demand mission-centric Lunar PNT swarm capable of providing PNT services to these low-cost lunar missions [10]. As upcoming Lunar missions will often operate at low-altitude about 30 km to 100 km for scientific observations and mapping purposes, the low-altitude orbital constellations could be employed to create an ad-hoc Lunar PNT system. However, several issues must be addressed, such as the instability of these orbits, which often require maintenance or are only suitable for short-duration missions, operating for fewer than 90 days. Additionally, at an altitude of 100 km, the satellites have a limited period during which they are above the horizon and capable of providing PNT service to users. The implementation of a non-dedicated, ad-hoc Lunar navigation constellation facilitates on-demand PNT services. A preliminary study of ad-hoc Lunar PNT system was conducted using 21 spacecraft in 5,5000 km altitude frozen orbits to test its feasibility and a basic performance of orbital asset localization among ad-hoc Lunar constellations in small satellites format [10]. These swarm assets are designed for autonomous localization with minimal Earth interaction, reducing dependency on bandwidth and ground resources. The design in [10] demonstrated the feasibility of a decentralized PNT approach, specifically employing a DEKF approach for state estimation, which helps minimize onboard operating costs. The DEKF method distributes computation across individual satellites, which lightens the computational load while maintaining accuracy in orbit ephemeris and clock offsets, similar to centralized systems [11]. In a follow-on study [12], each spacecraft was limited to 2 communications antennae, forcing the selection of measurements and scheduling spacecraft activities to perform the measurements. A matching algorithm is implemented to select the best measurements and schedule position estimation updates. The decentralized localization performance is also investigated with increasing levels of network degradation for swarm assets considering the impact of intermittent and permanent communication failure, to demonstrate the robustness and fidelity of the decentralized Lunar PNT service [13]. This study confirmed that the ad-hoc PNT constellations in frozen orbit are highly robust and resilient to communication failures. However, unlike frozen orbit swarm assets, the low-altitude satellites have a limited ground view at an altitude of 100 km, where the ad-hoc Lunar constellation consists of 98 low-altitude satellites, evenly distributed across seven circular polar orbital planes, alongside two satellites in a frozen orbit at an altitude of 5,500 km (Figure 1). Therefore, the number of satellites visible to ground users is significantly limited in low-altitude orbit constellations. As each visibility of a spacecraft remains intact for only a few ticks before it moves out of the field of view, the ground user encounters challenges in maintaining continuous navigation service, resulting in sparse availability and provision of Lunar PNT system. Consequently, service availability is primarily restricted to the Lunar South Pole region (Figure 2). Given these limitations and concerns, the localization performance of low-altitude swarm assets will be assessed in this study. We focus on the investigation of the localization performance of low-altitude swarm assets and ground users near the Lunar South Pole. The overall flow of the Lunar PNT simulation incorporates the DEKF approach of asset localization and the weighted least-squares approach in user localization (Figure 3). The autonomous Lunar PNT simulation is primarily implemented in MATLAB, where the DEKF based on the matching scheduler is implemented with Google’s OR-tools as a model builder and Gurobi optimization tool as a backend solver. The General Mission Analysis Tool (GMAT) is utilized to generate ephemeris data for swarm assets, and accounts for satellite orbital details, mass, and perturbations like solar radiation pressure and drag coefficients. Each ephemeris dataset is produced in the Moon International Celestial Reference Frame (ICRF) inertial coordinate system. For state estimation, the distributed swarm assets rely on two-way Inter-Satellite Link (ISL) measurements, which involve tracking pseudoranges and relative velocities between visible satellites and anchor nodes during each observation. Numerical evaluations of the decentralized localization process are conducted to demonstrate the feasibility of the low-altitude PNT system in providing reliable navigation services. The main approach involves using DEKF and CEKF to localize 100 satellites in low-altitude constellations, where the CEKF is implemented to serve as a baseline for comparing the performance of distributed algorithms. In both cases, we evaluate ‘fully sampled’ measurements from all available assets, and ‘two ISL’ measurements when spacecraft are constrained to have only two antennas. We test four estimation techniques: CEKF fully sampled, CEKF two ISL, DEKF fully sampled, and DEKF two ISL filters. As the DEKF update cycle is comprised of network setup, communication, and computations, a global broadcast network and 2-way ISL network setup will take from 4 to 6 minutes as maximum [12]. In this simulation, the DEKF update cycle is set to 10 minutes, including a 4-minute latency for obtaining and computing the actual measurement updates. We experiment an increased update cycle to demonstrate the feasibility and evaluate the impact on localization performance using various tuning values for measurement noise covariances (Figures 4 and 5). By comparing centralized and decentralized approaches using a matching algorithm, we analyze the influence of cross-correlation factors in the covariance matrix, assuming 100% reliability of all assets and measurements. The increased frequency and the adjustments of tuning parameters reveal distinct error patterns between the two scenarios. The localization accuracy of the swarm assets and ground users is assessed by taking the median error across 100 assets and one ground user (84.9°S, 137.5°E) over 7-day simulation period (Table 1). Since the user localization accuracy is significantly affected by the performance of the swarm assets, it is crucial to maintain high localization accuracy within the swarm. This study will continue to explore decentralized filtering for autonomous LPNT operations, with further investigation of an 'iterative' matching approach which enumerates every valid matching pair, planned for the following month.

Yeji Kim↗

Hydrodynamization times of a holographic fluid far from equilibrium

Here we investigate several hydrodynamization times for an ensemble of different far-from-equilibrium solutions of the strongly coupled $\mathscr{N}$ = 4 supersymmetric Yang-Mills plasma undergoing Bjorken flow. For the ensemble of initial data analyzed in the present work, we find that, with typical tolerances between 3% to 5%, the average hydrodynamization time associated with the late time convergence of the pressure anisotropy to the corresponding Borel resummed hydrodynamic attractor is approximately equal to the average hydrodynamization time associated with the Navier-Stokes result, while both are shorter than the average hydrodynamization time associated with second-order hydrodynamics. On the other hand, we find that the entropy density of the different solutions coalesces to second-order hydrodynamics long before entering in the Navier-Stokes regime. A clear hierarchy between the different average hydrodynamization times of the Bjorken expanding fluid is established for the set of analyzed initial data, comprising also some solutions which, whilst satisfying the dominant and the weak energy conditions at the initial time, evolve such as to transiently violate one or both conditions when the fluid is still far from equilibrium. In particular, solutions violating the weak energy condition are generally found to take a longer time to enter in the hydrodynamic regime than the other solutions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Parity solution to the strong CP problem and a unified framework for inflation, baryogenesis, and dark matter

It has been known for some time that asymptotic parity invariance of weak interactions can provide a solution to the strong CP problem without the need for the axion. Left-right symmetric theories which employ a minimal Higgs sector consisting of a left-handed and a right-handed doublet is an example of such a theory wherein all fermion masses arise through a generalized seesaw mechanism. In this paper we present a way to understand the origin of matter-antimatter asymmetry as well as the dark matter content of the universe in these theories using the Affleck-Dine (AD) leptogenesis mechanism and inflaton decay, respectively. Three gauge singlet fermions are needed for this purpose, two of which help to implement the Dirac seesaw for neutrino masses while the third one becomes the non-thermal dark matter candidate. A soft lepton number breaking term involving the AD scalar field is used to generate lepton asymmetry which suffers no wash-out effects and maintains the Dirac nature of neutrinos. This framework thus provides a unified description of many of the unresolved puzzles of the standard model that require new physics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A species scale-driven breakdown of effective field theory in time-dependent string backgrounds

We present a novel way in which effective field theory (EFT) can break down in cosmological string backgrounds depending on the behavior of the quantum gravity cutoff in infinite distance limits, known as the species scale $Λ_s$. Namely, EFT can break down if the species scale $Λ_s$ falls off so rapidly as the Friedmann-Robertson-Walker (FRW) scale factor grows from some initial value $a_i$ to some final value $a_f$ that the physical momentum of an initial Hubble-sized perturbation $\sim H_i^{-1}$ grows to exceed the species scale. For EFT to remain valid, a new condition $H_i \frac{a_i}{a_f} \ll Λ_{s,f}$ must hold, which is distinct from Trans-Planckian conditions discussed in the literature. Using the universal relation $\frac{\nabla m}{m} \cdot \frac{\nabla Λ_s}{Λ_s} = \frac{1}{d-2}$ in the infinite distance limits of moduli space where $m$ is the mass scale of the lightest tower and $\nabla$ measures variations with respect to the canonical metric on moduli space, we show that spatially flat FRW solutions in the string landscape violate this condition or at best marginally satisfy it. However, we find that sufficiently large negative spatial curvature always avoids a breakdown. To avoid EFT breakdown, we derive an upper bound on the duration of quasi-de Sitter expansion that classically evolves to decelerated expansion. Our bound is proportional to the Trans-Planckian Censorship Conjecture (TCC) bound, with the advantage that it applies to any FRW solution in the string landscape. Finally, we distinguish EFT breakdown from TCC violation, the latter being a quantum gravity constraint rather than an EFT limitation. Perhaps our most surprising finding is that in any flat FRW solution that develops a weakly coupled string at future infinity the EFT inevitably breaks down.

Cosmology and Nongalactic Astrophysics (astro-ph.C↗

Electrochemical generation of oxygen. 1: The effects of anions and cations on hydrogen chemisorption and anodic oxide film formation on platinum electrode. 2: The effects of anions and cations on oxygen generation on platinum electrode

The effects were studied of anions and cations on hydrogen chemisorption and anodic oxide film formation on Pt by linear sweep voltammetry, and on oxygen generation on Pt by potentiostatic overpotential measurement. The hydrogen chemisorption and anodic oxide film formation regions are greatly influenced by anion adsorption. In acids, the strongly bound hydrogen occurs at more cathodic potential when chloride and sulfate are present. Sulfate affects the initial phase of oxide film formation by produced fine structure while chloride retards the oxide-film formation. In alkaline solutions, both strongly and weakly bound hydrogen are influenced by iodide, cyanide, and barium and calcium cations. These ions also influence the oxide film formation. Factors considered to explain these effects are discussed. The Tafel slope for oxygen generation was found to be independent on the oxide thickness and the presence of cations or anions. The catalytic activity indicated by the exchange current density was observed decreasing with increasing oxide layer thickness, only a minor dependence on the addition of certain cations and anions was found.

Huang, C. J.↗

Space Shuttle Main Engine development /May 1980/

The Space Shuttle Main Engine is the most advanced liquid propellant rocket engine ever developed. During the last year, significant progess has been made toward the identification and solution of several design weaknesses. The initial two test series to certify the engine's readiness for manned flight have been completed, and additional test series are in progress. Operation of the total propulsion system has been demonstrated by a successful series of flight-simulation tests using the three-engine cluster Main Propulsion Test Article. The first set of flight engines has been acceptance tested for use in Shuttle Orbiter 102 (Columbia) for the first manned orbital flight in late 1980 or early 1981. Engine development is continuing with initial demonstrations of full-power level (109%) operation capability.

Sanchini, D. J.↗

Collisional redistribution of radiation. III - The equation of motion for the correlation function and the scattered spectrum

Computations were made of the scattering of monochromatic radiation by a degenerate atom in the binary-collision approximation for field strengths whose products of the Rabi frequency for atomic transition and the duration of a strong collision are much less than 1. An expression of motion for the correlation function is derived which does not exclude the region where thermal correlations may be neglected; the equation is valid outside the quantum-regression regime, and has a straightforward solution for practical cases. Solutions for the weak-field linear response regime are presented in terms of generalized absorption and emission profiles which depend on the indices of the atomic multipoles.

Burnett, K.↗

Numerical model of long-lived Jovian vortices

The extension of the measured zonal velocity profile into the adiabatic interior of Jupiter, while eddies and large oval structures are confined to a shallow stably-stratified upper layer, are assumed in a nonlinear numerical model of long-lived Jovian vortices. In agreement of the observed flows of Jupiter, each vortex is stationary with respect to the shear flow at a critical latitude that is close to the latitude of the vortex center. The solutions obtained are strongly nonlinear, in contrast to the solitary wave solutions that are the weakly nonlinear extensions of ultralong linear waves. The merging of two stable vortices upon collision, rather than the non-interaction predicted by solitary wave theory, is in keeping with Jovian vortex observations. It is suggested that long-lived vortices maintain themselves against dissipation by absorbing smaller vortices produced by convection.

Ingersoll, A. P.↗

Inertial Range Dynamics in Boussinesq Turbulence

L'vov and Falkovich have shown that the dimensionally possible inertial range scaling laws for Boussinesq turbulence, Kolmogorov and Bolgiano scaling, describe steady states with constant flux of kinetic energy and of entropy respectively. These scaling laws are treated as similarity solutions of the direct interaction approximation for Boussinesq turbulence. The Kolmogorov scaling solution corresponds to a weak perturbation by gravity of a state in which the temperature is a passive scalar but in which a source of temperature fluctuations exists. Using standard inertial range balances, the renormalized viscosity and conductivity, turbulent Prandtl number, and spectral scaling law constants are computed for Bolgiano scaling.

Rubinstein, Robert↗

Weak-wave analysis of shock interaction with a slipstream

A weak wave analysis of shock interaction with a slipstream is presented. The theory is compared to that for the acoustic case and to the exact nonlinear analysis. Sample calculations indicate that the weak wave theory yields a good approximation to the exact solution when the shock waves are sufficiently weak that the associated entropy increase is negligible. A qualitative discussion of the case of counterflowing streams is also included.

Barger, Raymond L.↗

A simple model of a strong shock driven by a spherical or cylindrical piston

A simple model of piston-driven spherical and cylindrical shocks is suggested. The model is based on a consistent use of two factors: (a) an almost uniform pressure across the shocked layer and, (b) continuous geometrical stretching of the surface elements of the expanding piston. It turns out that for a uniform pre-shock medium the gas between the piston and the shock behaves essentially as an incompressible fluid. An algebraic equation for the shock vs piston position is obtained. Detailed evaluation of the accuracy of the proposed solution shows that its accuracy is a few percent for the adiabatic index γ = 5/3 (as in ideal plasma). A closed-form solution describing enhancement of a weak ambient magnetic field by the shock is presented. The proposed model of piston-driven shocks goes beyond the classical self-similar solutions in that it: (1) naturally covers an early, non-asymptotic dynamics and its transition to asymptotic regime; (2) allows for smooth radial density variation of an ambient gas of the form of bumps, dimples or ramps between two constant values; and (3) allows for smooth temporal variation of piston velocity of the form of bumps, dimples, or ramps. Furthermore, this simple and versatile model provides some new insights into a classical hydrodynamical problem.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Conservation of quasiparticles in weakly turbulent plasmas.

An expansion making use of the eikonal is shown to yield a solution to the equations of a general weak turbulent plasma which is weakly dependent on space and time. The method is used to derive a quasi-particle conservation equation for quasi-static perturbations of Vlasov plasmas with general equilibrium field configurations.

Bloomberg, H. W.↗

An intelligent decomposition approach for efficient design of non-hierarchic systems

The design process associated with large engineering systems requires an initial decomposition of the complex systems into subsystem modules which are coupled through transference of output data. The implementation of such a decomposition approach assumes the ability exists to determine what subsystems and interactions exist and what order of execution will be imposed during the analysis process. Unfortunately, this is quite often an extremely complex task which may be beyond human ability to efficiently achieve. Further, in optimizing such a coupled system, it is essential to be able to determine which interactions figure prominently enough to significantly affect the accuracy of the optimal solution. The ability to determine 'weak' versus 'strong' coupling strengths would aid the designer in deciding which couplings could be permanently removed from consideration or which could be temporarily suspended so as to achieve computational savings with minimal loss in solution accuracy. An approach that uses normalized sensitivities to quantify coupling strengths is presented. The approach is applied to a coupled system composed of analysis equations for verification purposes.

Bloebaum, Christina L.↗

A Simple Theory of Capillary-Gravity Wave Turbulence

Employing a recently proposed 'multi-wave interaction' theory [JFM, 243, 623-625], spectra of capillary-gravity waves are derived. This case is characterized by a rather high degree of nonlinearity and a complicated dispersion law. The resultant absence of scale invariance makes this and some other problems of wave turbulence (e.g., nonlinear Rossby waves) intractable by small-perturbation techniques, even in the weak turbulence limit. The analytical solution obtained in the present work is shown to be in good agreement with experimental data. Its low- and high-frequency limits yield power-laws characterizing spectra of purely gravity and capillary waves, respectively. In the limits of weak and strong linearity, these reduce of the Zakharov-Filonenko and Phillips spectra, respectively.

Glazman, Roman E.↗

A Perturbative Solution for Nonlinear Stratified Upwelling over a Frictional Slope

Abstract A perturbative solution of simplified primitive equations for nonlinear weakly stratified upwelling over a frictional slope is found that resolves the vertical structure of velocity fields and can satisfy Ertel’s potential vorticity conservation in the stratified inviscid interior. The solution uses assumptions consistent with the model proposed by Lentz and Chapman, including a steady-state, constant cross-shore density gradient, no alongshore gradients, laterally inviscid, and consideration of cross-shore advection of alongshore momentum. The solution resolves the vertical structure of velocity fields (including subsurface maxima of compensational flow, not resolved by Lentz and Chapman) and can satisfy Ertel’s potential vorticity conservation in the stratified inviscid interior. The dynamics are similar to Lentz and Chapman; bottom stress balances alongshore wind stress in a homogeneous density ocean and is replaced by nonlinear cross-shore transport of alongshore momentum as the Burger number (S=αN/f, whereα,N, andfare the bottom slope, buoyancy frequency, Coriolis frequency, respectively) increases. When the solution uses the empirical relation between cross-shore and vertical density gradients proposed by Lentz and Chapman, vorticity conservation is not satisfied and the nonlinear momentum transport estimated by the solution linearly increases withS, asymptotically matching Lentz and Chapman forS< 1. When the solution conserves interior potential vorticity, the momentum transport is proportional toS 2 forS< 1 and is in better agreement with numerical simulations.

Oceanography↗

A Taylor weak-statement algorithm for hyperbolic conservation laws

Finite element analysis, applied to computational fluid dynamics (CFD) problem classes, presents a formal procedure for establishing the ingredients of a discrete approximation numerical solution algorithm. A classical Galerkin weak-statement formulation, formed on a Taylor series extension of the conservation law system, is developed herein that embeds a set of parameters eligible for constraint according to specification of suitable norms. The derived family of Taylor weak statements is shown to contain, as special cases, over one dozen independently derived CFD algorithms published over the past several decades for the high speed flow problem class. A theoretical analysis is completed that facilitates direct qualitative comparisons. Numerical results for definitive linear and nonlinear test problems permit direct quantitative performance comparisons.

Baker, A. J.↗