Search NASA⌕ Search

SEARCH · Search NASA

Results for “spatially explicit modeling”

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 235 records · Page 13

Assessing Several Modeling Approaches in Depletion Studies of a TRISO-Fueled Microreactor

The impact of MCNP6 depletion resolution on core lifetime is examined in the context of the Snowflake microreactor with explicit TRISO fuel. The change in core lifetime and isotope mass as a result of different tracked isotopes, timesteps, and spatial regions is discussed. Calculation speed of a prototype MCNP delta tracking module is compared to a reactivity equivalent physical transformation (RPT). Using a single depletion material underpredicts core lifetime by 15%, and the resolution necessary to converge isotope mass greatly depends on the specific isotope, in addition to the size (and location) of the depletion region. The prototype delta tracking module decreases the CPU time of explicit TRISO criticality calculations by 30%, but does not always result in a speedup when used with depletion. Significant depletion speedup is obtained using RPT (50% faster), and all isotope masses agreed within three percent.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

The Use of RGPS Kinematic Data to Estimate Nonlinear Sea Ice Motion

In current simulations of the interaction between sea ice and its environment, large significance is placed on the deformation of the sea ice. Sea ice deformation is an important process in determining the sea ice thickness distribution across a wide range of space and time scales. Changes in the sea ice thickness distribution affect energy and mass fluxes between the atmosphere and ocean and also the strength of the ice. While most current ice models assume linear variation in the ice motion field to calculate strain, deformation of sea ice occurs through the opening, closing and shearing of ice along discrete linear features. New numerical models are being developed which explicitly account for discontinuities in ice motion, and the need for requisite data sets for model validation has emerged. Multiple buoy data sets, as well as satellite data, have been used to examine the movement and deformation of sea ice. Generally it has been found that the ice motion field has been represented well by buoy data, as well as satellite data over a broad range of scales. However, the underlying deformation (spatial variation in displacement) as represented by different data sets may vary. For the work presented here, sea ice motion In current simulations of the interaction between sea ice and its environment, large significance is placed on the deformation of the sea ice. Sea ice deformation is an important process in determining the sea ice thickness distribution across a wide range of space and time scales. Changes in the sea ice thickness distribution affect energy and mass fluxes between the atmosphere and ocean and also the strength of the ice. While most current ice models assume linear variation in the ice motion field to calculate strain, deformation of sea ice occurs through the opening, closing and shearing of ice along discrete linear features. New numerical models are being developed which explicitly account for discontinuities in ice motion, and the need for requisite data sets for model validation has emerged. Multiple buoy data sets, as well as satellite data, have been used to examine the movement and deformation of sea ice. Generally it has been found that the ice motion field has been represented well by buoy data, as well as satellite data over a broad range of scales. However, the underlying deformation (spatial variation in displacement) as represented by different data sets may vary. For the work presented here, sea ice motio

Pruis, M.↗

Explicitly Resolving Lightning and Electrification Processes from the 10-12 April 2019 Thundersnow Outbreak

The 10-12 April 2019 thundersnow (i.e., lightning within snowfall) outbreak was examined via ground- and space-based lightning observations and was simulated using a numerical weather prediction model with an explicit electrification parameterization. When compared to observations, the simulation propagated the synoptic snowband two to six hours faster while also exaggerating the 3-D reflectivity structure. Throughout the event, the simulation produced 1,733 thundersnow flashes which was less than what was observed by ground- and space-based lightning sensors. In general, simulated thundersnow flashes were spatially offset from the largest reflectivities within the synoptic snowband and tended to occur within elevated convection that traversed isentropically along the top of mid-level frontogenesis. These simulated thundersnow flashes were associated with a tripole charge structure with ice/snow hydrometeors contributing most to the main negative charge region. Both simulated and observed thundersnow flashes initiated in conditionally unstable environments. Lastly, a conceptual model was developed to explain the spatial separation between the largest reflectivities in the snowband and the occurrence of thundersnow. It is hypothesized that the spatial offset of thundersnow initiation from the reflectivity cores within the synoptic snowband arose from a thermal circulation – induced by mid-level frontogenesis – that advects positively charged ice/snow hydrometeors towards the surface and creates a nearly homogeneous vertical charge structure.

lightning↗

Integrated Bosch Process System Models for In-Situ Oxygen and Carbon Production

In-Situ Resource Utilization (ISRU) technology is a vital component to NASA’s mission of a sustainable presence on the Moon and Mars. Local resources can be leveraged to reduce resupply frequency and mass. Elements of the Bosch process, combined with the carbothermal reduction process, can produce oxygen on the lunar surface with minimal consumables. The Bosch process can also produce oxygen on the Martian surface by using the CO 2 -rich environment. Between both systems, adsorption pump, solar thermal energy, carbon formation reactor, and water recovery subsystems are modeled and integrated to create a functional model in MATLAB software. The model is used to simulate performance of the system and reduce mass, power, and volume requirements. This integrated system model provides a tool to scale ISRU technologies for oxygen and carbon production. The MATLAB model is created by developing a system of independent subsystem models that are solved for their quasi-steady state values which can be integrated with respect to time to determine the change in current states. A flexible time stepping method is used to ensure a high level of accuracy during periods of rapid change while still making use of a simple explicit integration method. The flexible time step is calculated for each independent subsystem and the minimum value from those is used as the overall time step. A flexible time step is calculated by dividing a resolution value, or the maximum change per time step, by the variables current rate of change. The maximum value from all points in space is used for subsystem models that contain multiple values. The process is done for every variable that is being monitored in each subsystem and the global minimum is used as that iteration’s timestep. Several assumptions used in the MATLAB model for fluid flow dynamics, such as 1-D gas flow through the sorption pump, are supported by modeling in Ansys Fluent software. The Lunar oxygen production system is outlined in Fig. 1. The carbothermal reduction subsystem uses solar energy to heat a mixture of lunar regolith and carbon powder to produce carbon monoxide. To begin, the carbon monoxide feeds to the modified Bosch subsystem along with hydrogen gas. The reactants then enter the carbon formation reactor where water and carbon powder are produced. Solar thermal energy is used to add energy to the reactor, but waste heat from the carbothermal process is another potential heat source. The water is collected and electrolyzed to produce hydrogen which reenters the Bosch subsystem, and the oxygen is stored for downstream use. The carbon powder is collected and feeds back into the carbothermal subsystem. The Martian oxygen production system uses the full Bosch process and is outlined in Fig 2. A CO 2 adsorption pump thermally cycles to scrub and pressurize CO 2 from the environment. Along with an initial supply of hydrogen, the reactants enter the Reverse Water Gas Shift Reactor (RWGSR) which produces carbon monoxide and water. Carbon monoxide and unreacted hydrogen enter the carbon formation reactor to produce water and carbon powder. The water is collected from both reactors and electrolyzed to reintroduce hydrogen and store oxygen for propellant production or life support. Carbon is removed from the carbon formation reactor and stored. The adsorption pump utilizes rapid cycle temperature swings within a stack of zeolite coated surfaces. The subsystem model solves 1-D quasi-steady conservation laws of the quasi-steady form, shown in Eq. 1, for the gas stream and heat exchange liquid to predict performance parameters such as breakthrough capacity and optimum cycle time. The source term S is used to capture interactions between the fluid flows and the sorbent. A quasi-steady-state scheme is used where no time derivatives appear in the governing equations, except for those in the source terms. This results in an autonomous system, where ∂F/∂x = ƒ(F). The fluxes F are provided at the inlet, and an explicit method is used to solve for the spatial distribution of F. The heat and mass flows to the sorbent are then extracted from the source terms. These flows are numerically integrated to produce a 1-D solution for the system’s state as a function of both time and space. The body of the adsorption pump is separated into two semi-independent models: the heat exchanger fluid flow and gas flow through the zeolite coated surfaces. Both models are solved using the above-described method to find a 1-D solution as a function of space and interact only once a timestep is taken. The interaction point is the sorbent through which all heat transfer between the two models must occur. Sorbent mass adsorption is calculated using the Lagergren model, shown in Eq. 2, where the transfer coefficient, λ D , is found by solving a system of nondimensionalized equations derived by using the heat and mass transfer analogy for transport phenomena. Using Grade 544 Type 13X zeolite as the sorbent material, the equilibrium concentration, θ eq , is calculated using the k-site Langmuir isotherm and fit parameters. Additionally, the enthalpy of adsorption used in the model is computed by interpolation of available data [1]. The subsystem model was validated using the Rapid Cycle Temperature Swing Adsorption (RC-TSA) pump. The solar thermal energy subsystem focuses on a solar concentrator concept with a heat exchanger to heat the reactants before entering the carbon formation reactor. The subsystem model assumes a fixed solar flux and reflector efficiency to calculate the reactant temperature given the incoming temperature, pressure, and exchanger geometry. The receiver is a custom manufactured series of copper blocks with serpentine channels to increase its surface area and the residence time of the reactants to heat up to 550 °C. The subsystem model was validated using a heat exchanger developed at NASA Glenn Research Center (GRC). The solar thermal energy subsystem focuses on a solar concentrator concept with a heat exchanger to heat the reactants before entering the carbon formation reactor. The subsystem model assumes a fixed solar flux and reflector efficiency to calculate the reactant temperature given the incoming temperature, pressure, and exchanger geometry. The receiver is a custom manufactured series of copper blocks with serpentine channels to increase its surface area and the residence time of the reactants to heat up to 550 °C. The subsystem model was validated using a heat exchanger developed at NASA Glenn Research Center (GRC).

In situ Resource Utilization↗

Some factors influencing radiation of sound from flow interaction with edges of finite surfaces

Edges of surfaces which are exposed to unsteady flow cause both strictly acoustic effects and hydrodynamic effects, in the form of generation of new hydrodynamic sources in the immediate vicinity of the edge. An analytical model is presented which develops the explicit sound-generation role of the velocity and Mach number of the eddy convection past the edge, and the importance of relative scale lengths of the turbulence, as well as the relative intensity of pressure fluctuations. The Mach number (velocity) effects show that the important paramater is the convection Mach number of the eddies. The effects of turbulence scale lengths, isotropy, and spatial density (separation) are shown to be important in determining the level and spectrum of edge sound radiated for the edge dipole mechanism. Experimental data is presented which provides support for the dipole edge noise model in terms of Mach number (velocity) scaling, parametric dependence on flow field parameter, directivity, and edge diffraction effects.

Hayden, R. E.↗

Deriving spatially explicit direct and indirect interaction networks from animal movement data

Abstract Quantifying spatiotemporally explicit interactions within animal populations facilitates the understanding of social structure and its relationship with ecological processes. Data from animal tracking technologies (Global Positioning Systems [“GPS”]) can circumvent longstanding challenges in the estimation of spatiotemporally explicit interactions, but the discrete nature and coarse temporal resolution of data mean that ephemeral interactions that occur between consecutive GPS locations go undetected. Here, we developed a method to quantify individual and spatial patterns of interaction using continuous‐time movement models (CTMMs) fit to GPS tracking data. We first applied CTMMs to infer the full movement trajectories at an arbitrarily fine temporal scale before estimating interactions, thus allowing inference of interactions occurring between observed GPS locations. Our framework then infers indirect interactions—individuals occurring at the same location, but at different times—while allowing the identification of indirect interactions to vary with ecological context based on CTMM outputs. We assessed the performance of our new method using simulations and illustrated its implementation by deriving disease‐relevant interaction networks for two behaviorally differentiated species, wild pigs ( Sus scrofa ) that can host African Swine Fever and mule deer ( Odocoileus hemionus ) that can host chronic wasting disease. Simulations showed that interactions derived from observed GPS data can be substantially underestimated when temporal resolution of movement data exceeds 30‐min intervals. Empirical application suggested that underestimation occurred in both interaction rates and their spatial distributions. CTMM‐Interaction method, which can introduce uncertainties, recovered majority of true interactions. Our method leverages advances in movement ecology to quantify fine‐scale spatiotemporal interactions between individuals from lower temporal resolution GPS data. It can be leveraged to infer dynamic social networks, transmission potential in disease systems, consumer–resource interactions, information sharing, and beyond. The method also sets the stage for future predictive models linking observed spatiotemporal interaction patterns to environmental drivers.

59 BASIC BIOLOGICAL SCIENCES↗

A function space approach to state and model error estimation for elliptic systems

An approach is advanced for the concurrent estimation of the state and of the model errors of a system described by elliptic equations. The estimates are obtained by a deterministic least-squares approach that seeks to minimize a quadratic functional of the model errors, or equivalently, to find the vector of smallest norm subject to linear constraints in a suitably defined function space. The minimum norm solution can be obtained by solving either a Fredholm integral equation of the second kind for the case with continuously distributed data or a related matrix equation for the problem with discretely located measurements. Solution of either one of these equations is obtained in a batch-processing mode in which all of the data is processed simultaneously or, in certain restricted geometries, in a spatially scanning mode in which the data is processed recursively. After the methods for computation of the optimal estimates are developed, an analysis of the second-order statistics of the estimates and of the corresponding estimation error is conducted. Based on this analysis, explicit expressions for the mean-square estimation error associated with both the state and model error estimates are then developed.

Rodriguez, G.↗

A function space approach to state and model error estimation for elliptic systems

An approach is advanced for the concurrent estimation of the state and of the model errors of a system described by elliptic equations. The estimates are obtained by a deterministic least-squares approach that seeks to minimize a quadratic functional of the model errors, or equivalently, to find the vector of smallest norm subject to linear constraints in a suitably defined function space. The minimum norm solution can be obtained by solving either a Fredholm integral equation of the second kind for the case with continuously distributed data or a related matrix equation for the problem with discretely located measurements. Solution of either one of these equations is obtained in a batch-processing mode in which all of the data is processed simultaneously or, in certain restricted geometries, in a spatially scanning mode in which the data is processed recursively. After the methods for computation of the optimal esimates are developed, an analysis of the second-order statistics of the estimates and of the corresponding estimation error is conducted. Based on this analysis, explicit expressions for the mean-square estimation error associated with both the state and model error estimates are then developed. While this paper focuses on theoretical developments, applications arising in the area of large structure static shape determination are contained in a closely related paper (Rodriguez and Scheid, 1982).

Rodriguez, G.↗

Evolution of solitary density waves in stellar winds of early-type stars: A simple explanation of discrete absorption component behavior

We model the evolution of a density shell propagating through the stellar wind of an early-type star, in order to investigate the effects of such shells on UV P Cygni line profiles. Unlike previous treatments, we solve the mass, momentum, and energy conservation equations, using an explicit time-differencing scheme, and present a parametric study of the density, velocity, and temperature response. Under the assumed conditions, relatively large spatial scale, large-amplitude density shells propagate as stable waves through the supersonic portion of the wind. Their dynamical behavior appears to mimic propagating 'solitary waves,' and they are found to accelerate at the same rate as the underlying steady state stellar wind (i.e., the shell rides the wind). These hydrodynamically stable structures quantitatively reproduce the anomalous 'discrete absorption component' (DAC) behavior observed in the winds of luminous early-type stars, as illustrated by comparisons of model predictions to an extensive International Ultraviolet Explorer (IUE) time series of spectra of zeta Puppis (O4f). From these comparisons, we find no conclusive evidence indicative of DACs accelerating at a significantly slower rate than the underlying stellar wind, contrary to earlier reports. In addition, these density shells are found to be consistent within the constraints set by the IR observations. We conclude that the concept of propagating density shells should be seriously reconsidered as a possible explanation of the DAC phenomenon in early-type stars.

Waldron, Wayne L.↗

Finite element algorithms for compressible flow computation on a supercomputer

Methods of applying computational fluid dynamics techniques to engineering problems are discussed. The linking factors between theoretical and applied research are the shape and weight functions which produce the Galerkin finite element schemes. The flow is unsteady, viscous, heat conducting and compressible, and steady-state flows are the asymptotic limit of unsteady flows. All flows are described by the time-averaged Navier-Stokes equations (NSE) with mass and energy conservation. Local curvilinear intrinsic coordinates are applied to discretize the NSE in arbitrary geometric domains. A modified weighted residuals approach defines the discrete analogs of the physical systems modeled and a general interpolants method is used to derive families of numerical models, both implicit and explicit, finite difference and finite elemnt, from a single point of departure. Time-dependent element approximations are achieved with a progressive assembly of generalized elements method. Hyperbolic steady-state Euler algorithms and quasi-parabolic are spatial marching algorithms for solving the discretized equations. Application of the techniques is illustrated through calculation of the compressible flow around a hypersonic flight vehicle at 100,000 ft at zero angle of attack at speeds of Mach 4-10.

Spradley, L. W.↗

Modeling of outgassing and matrix decomposition in carbon-phenolic composites

Work done in the period Jan. - June 1994 is summarized. Two threads of research have been followed. First, the thermodynamics approach was used to model the chemical and mechanical responses of composites exposed to high temperatures. The thermodynamics approach lends itself easily to the usage of variational principles. This thermodynamic-variational approach has been applied to the transpiration cooling problem. The second thread is the development of a better algorithm to solve the governing equations resulting from the modeling. Explicit finite difference method is explored for solving the governing nonlinear, partial differential equations. The method allows detailed material models to be included and solution on massively parallel supercomputers. To demonstrate the feasibility of the explicit scheme in solving nonlinear partial differential equations, a transpiration cooling problem was solved. Some interesting transient behaviors were captured such as stress waves and small spatial oscillations of transient pressure distribution.

Mcmanus, Hugh L.↗

Computing aerodynamic sound using advanced statistical turbulence theories

It is noted that the calculation of turbulence-generated aerodynamic sound requires knowledge of the spatial and temporal variation of Q sub ij (xi sub k, tau), the two-point, two-time turbulent velocity correlations. A technique is presented to obtain an approximate form of these correlations based on closure of the Reynolds stress equations by modeling of higher order terms. The governing equations for Q sub ij are first developed for a general flow. The case of homogeneous, stationary turbulence in a unidirectional constant shear mean flow is then assumed. The required closure form for Q sub ij is selected which is capable of qualitatively reproducing experimentally observed behavior. This form contains separation time dependent scale factors as parameters and depends explicitly on spatial separation. The approximate forms of Q sub ij are used in the differential equations and integral moments are taken over the spatial domain. The velocity correlations are used in the Lighthill theory of aerodynamic sound by assuming normal joint probability.

Hecht, A. M.↗

Global Climate Forcing from Albedo Change Caused by Large-scale Deforestation and Reforestation: Quantification and Attribution of Geographic Variation

Large-scale deforestation and reforestation have contributed substantially to historical and contemporary global climate change in part through albedo-induced radiative forcing, with meaningful implications for forest management aiming to mitigate climate change. Associated warming or cooling varies widely across the globe due to a range of factors including forest type, snow cover, and insolation, but resulting geographic variation remain spoorly described and has been largely based on model assessments. This study provides an observation-based approach to quantify local and global radiative forcings from large-scale deforestation and reforestation and further examines mechanisms that result in the spatial heterogeneity of radiative forcing. We incorporate a new spatially and temporally explicit land cover-specific albedo product derived from Moderate Resolution Imaging Spectroradiometer with a historical land use data set (Land Use Harmonization product). Spatial variation in radiative forcing was attributed to four mechanisms, including the change in snow-covered albedo, change in snow-free albedo, snow cover fraction, and incoming solar radiation. We find an albedo-only radiative forcing (RF) of -0.819 W m(exp -2) if year 2000 forests were completely deforested and converted to croplands. Albedo RF from global reforestation of present-day croplands to recover year 1700 forests is estimated to be 0.161 W m)exp -2). Snow-cover fraction is identified as the primary factor in determining the spatial variation of radiative forcing in winter, while the magnitude of the change in snow-free albedo is the primary factor determining variations in summertime RF. Findings reinforce the notion that, for conifers at the snowier high latitudes, albedo RF diminishes the warming from forest loss and the cooling from forest gain more so than for other forest types, latitudes, and climate settings.

Jiao, Tong↗

Conjugate Compressible Fluid Flow and Heat Transfer in Ducts

A computational approach to modeling transient, compressible fluid flow with heat transfer in long, narrow ducts is presented. The primary application of the model is for analyzing fluid flow and heat transfer in solid propellant rocket motor nozzle joints during motor start-up, but the approach is relevant to a wide range of analyses involving rapid pressurization and filling of ducts. Fluid flow is modeled through solution of the spatially one-dimensional, transient Euler equations. Source terms are included in the governing equations to account for the effects of wall friction and heat transfer. The equation solver is fully-implicit, thus providing greater flexibility than an explicit solver. This approach allows for resolution of pressure wave effects on the flow as well as for fast calculation of the steady-state solution when a quasi-steady approach is sufficient. Solution of the one-dimensional Euler equations with source terms significantly reduces computational run times compared to general purpose computational fluid dynamics packages solving the Navier-Stokes equations with resolved boundary layers. In addition, conjugate heat transfer is more readily implemented using the approach described in this paper than with most general purpose computational fluid dynamics packages. The compressible flow code has been integrated with a transient heat transfer solver to analyze heat transfer between the fluid and surrounding structure. Conjugate fluid flow and heat transfer solutions are presented. The author is unaware of any previous work available in the open literature which uses the same approach described in this paper.

Cross, M. F.↗

Truck platooning in the U.S. national road network: A system-level modeling approach

Truck platooning enables a group of trucks to move close together, which helps reduce truck fuel use and increase effective road capacity. In this paper, a system-level equilibrium model is developed to characterize spontaneous truck platooning with coexistence of non-platooning vehicles in a network, by explicitly accounting for the interlocking relationship among platoon formation time, truck fuel saving, and increase in effective road capacity. To equilibrate the relationships, an algorithm is proposed which involves a diagonalization approach and a bush based algorithm to solve decomposed subproblems. The condition of proportionality is imposed to obtain unique traffic flows for each class of vehicles on road links. In addition, a spatially constrained multivariate clustering technique is employed to construct origin/destination zones that are smaller than the coarse Freight Analysis Framework (FAF) zones, while maintaining reasonable computational burden for network traffic assignment. Model implementation in the U.S. shows that platooning could lead to 7.9% fuel saving among platoonable trucks in 2025 and a comparable increase in effective capacity of platoonable road links, which would account for 60% of rural interstate roads. The fuel saving and road capacity improvement translate into an annual cost reduction of $\$$868 million for the U.S. intercity trucking sector and reduced road infrastructure investment needs worth $\$$4.8 billion. Extensive sensitivity analysis further reveals that fuel saving of platoonable trucks increases with platoon size but decreases with inter-truck distance in a platoon. Fuel saving potential suggests that priority should be given to rural rather than urban roads in deploying platooning technologies. As we expected, greater market penetration of platooning technologies means higher fuel saving and greater increase in effective road capacity.

33 ADVANCED PROPULSION SYSTEMS↗

Justice 40 Tool (J40 Tool) v1.0

The Justice 40 tool provides a quantitative framework to support decision-making around equitable energy interventions at the community level. The tool calculates the optimal portfolio of policy interventions that explicitly mitigates energy insecurity of an eligible population, by reducing its disproportionate energy burden. The place-based analysis assumes a spatial census tract-level resolution and distinguishes different sociodemographic groups within each tract. Instead of focusing on a specific technology, the underlying J40 model captures the combined effect of a set of policy interventions, currently including weatherization, rooftop solar, community solar and community wind. For each tract, the model chooses the optimal (least cost) combination of interventions to address the disproportionate burden, considering the specific population demographics and techno-economic potentials of technologies. Mathematically, this problem is formalized as an optimization model and formulated as a linear program.

Heleno, Miguel↗

The Poisson tensor completion non-parametric differential entropy estimator

We introduce the Poisson tensor completion (PTC) estimator, a non-parametric differential entropy estimator. The PTC estimator leverages inter-sample relationships to compute a low-rank Poisson tensor decomposition of the frequency histogram. Our crucial observation is that the histogram bins are an instance of a space partitioning of counts and thus can be identified with a spatial Poisson process. The Poisson tensor decomposition leads to a completion of the intensity measure over all bins—including those containing few to no samples—and leads to our proposed PTC differential entropy estimator. A Poisson tensor decomposition models the underlying distribution of the count data and guarantees non-negative estimated values and so can be safely used directly in entropy estimation. Our estimator is the first tensor-based estimator that exploits the underlying spatial Poisson process related to the histogram explicitly when estimating the probability density with low-rank tensor decompositions for the purpose of tensor completion. Furthermore, we demonstrate that our PTC estimator is a substantial improvement over standard histogram-based estimators for sub-Gaussian probability distributions because of the concentration of norm phenomenon.

42 ENGINEERING↗

Development of a grid-independent approximate Riemannsolver

A grid-independent approximate Riemann solver for use with the Euler and Navier-Stokes equations was introduced and explored. The two-dimensional Euler and Navier-Stokes equations are described in Cartesian and generalized coordinates, as well as the traveling wave form of the Euler equations. The spatial and temporal discretization are described for both explicit and implicit time-marching schemes. The grid-aligned flux function of Roe is outlined, while the 5-wave grid-independent flux function is derived. The stability and monotonicity analysis of the 5-wave model are presented. Two-dimensional results are provided and extended to three dimensions. The corresponding results are presented.

Rumsey, Christopher Lockwood↗