Search NASA⌕ Search

SEARCH · Search NASA

Results for “Conservative variables”

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 19 records

Conservative-variable average states for equilibrium gas multi-dimensional fluxes

Modern split component evaluations of the flux vector Jacobians are thoroughly analyzed for equilibrium-gas average-state determinations. It is shown that all such derivations satisfy a fundamental eigenvalue consistency theorem. A conservative-variable average state is then developed for arbitrary equilibrium-gas equations of state and curvilinear-coordinate fluxes. Original expressions for eigenvalues, sound speed, Mach number, and eigenvectors are then determined for a general average Jacobian, and it is shown that the average eigenvalues, Mach number, and eigenvectors may not coincide with their classical pointwise counterparts. A general equilibrium-gas equation of state is then discussed for conservative-variable computational fluid dynamics (CFD) Euler formulations. The associated derivations lead to unique compatibility relations that constrain the pressure Jacobian derivatives. Thereafter, alternative forms for the pressure variation and average sound speed are developed in terms of two average pressure Jacobian derivatives. Significantly, no additional degree of freedom exists in the determination of these two average partial derivatives of pressure. Therefore, they are simultaneously computed exactly without any auxiliary relation, hence without any geometric solution projection or arbitrary scale factors. Several alternative formulations are then compared and key differences highlighted with emphasis on the determination of the pressure variation and average sound speed. The relevant underlying assumptions are identified, including some subtle approximations that are inherently employed in published average-state procedures. Finally, a representative test case is discussed for which an intrinsically exact average state is determined. This exact state is then compared with the predictions of recent methods, and their inherent approximations are appropriately quantified.

Iannelli, G. S.↗

Conserved variable analysis of the convective boundary layer thermodynamic structure over the tropical oceans

An analysis of FGGE dropwindsonde data using conserved thermodynamic variables shows mixing line structures for the convective boundary layer over the equatorial Pacific. Deeper boundary layers show a double structure. Reversals of the gradients of mixing ratio and equivalent potential temperature above the boundary-layer top are present in all the averages and suggest that the origin of the air sinking into the boundary layer needs further study.

Betts, Alan K.↗

Choice of Variables and Preconditioning for Time Dependent Problems

We consider the use of low speed preconditioning for time dependent problems. These are solved using a dual time step approach. We consider the effect of this dual time step on the parameter of the low speed preconditioning. In addition, we compare the use of two sets of variables, conservation and primitive variables, to solve the system. We show the effect of these choices on both the convergence to a steady state and the accuracy of the numerical solutions for low Mach number steady state and time dependent flows.

Turkel, Eli↗

Thermodynamic structure of the stratocumulus-capped boundary layer on 7 July, 1987

As part of project First ISCCP Regional Experiment (FIRE), a mission was carried out on 7 July 1987 to study the thermodynamic structure of a boundary layer which is in transition from a clear to a cloudy state. The National Center for Aeronautical Research (NCAR) Electra flew a pattern in tight coordination with the NASA ER-2 aircraft near 122 West, 31.6 North off the coast of California. A description is given here of the thermodynamic structure. The purpose is to derive the entrainment rate and the fluxes of the thermodynamic variables. To this end researchers represent the data in conserved variable diagrams.

Boers, Reinout↗

A generalized hyperbolic marching technique for three-dimensional supersonic flow with shocks

The current paper introduces a generalized, nonorthogonal coordinate system for application to supersonic flows with such large angles between streamlines of the initial Cauchy data and the marching direction that the velocity component in the marching direction is actually subsonic. Finite-difference operators are developed to solve the three-dimensional, steady, inviscid conservation equations of fluid flow referenced to this general frame. Moreover, a new method of aligning the difference mesh to the bow shock wave is presented which eliminates both the necessity for differencing the free-stream flow properties and the spurious fluctuations that often arise in conservative variables when they are differenced across the discontinuity. Results are given for the case of flow past a conical surface.

Rizzi, A. W.↗

Numerical methods for the calculation of three-dimensional nozzle exhaust flow fields

Numerical codes developed for the calculation of three-dimensional nozzle exhaust flow fields associated with hypersonic airbreathing aircraft are described. Both codes employ reference plane grid networks with respect to three coordinate systems. Program CHAR3D is a characteristic code utilizing a new wave preserving network within the reference planes, while program BIGMAC is a finite difference code utilizing conservation variables and a one-sided difference algorithm. Secondary waves are numerically captured by both codes, while the underexpansion shock and plume boundary are treated discretely. The exhaust gas properties consist of hydrogen-air combustion product mixtures in local chemical equilibrium. Nozzle contours are treated by a newly developed geometry package based on dual cubic splines. Results are presented for simple configurations demonstrating two- and three-dimensional multiple wave interactions.

Dash, S. M.↗

An implicit factored scheme for the compressible Navier-Stokes equations

An implicit finite difference scheme is developed for the numerical solution of the compressible Navier-Stokes equations in conservation-law form. The algorithm is second-order-time accurate, noniterative, and spatially factored. In order to obtain an efficient factored algorithm, the spatial cross-derivatives are evaluated explicitly. However, the algorithm is unconditionally stable and, although a three-time-level scheme, requires only two-time-levels of data storage. The algorithm is constructed in a 'delta' form (i.e., increments of the conserved variables and fluxes) that provides a direct derivation of the scheme and leads to an efficient computational algorithm. In addition, the delta form has the advantageous property of a steady-state (if one exists) independent of the size of the time step. Numerical results are presented for a two-dimensional shock boundary-layer interaction problem.

Beam, R. M.↗

Algorithm for Unsteady Potential Flow About Airfoils

Implicit finite-difference scheme efficiently computes unsteady potential flow about airfoils. Formulation uses density and velocity potential as dependent variables. Conservation form is retained to assure that shock wave location and speed are computed correctly. Scheme fills need for method to calculate efficiently unsteady potential flow about airfoils and to predict flutter and other unsteady aeroelastic phenomena in transonic flow regimes.

Chipman, R.↗

SSME Turbopump Turbine Computations

A two-dimensional viscous code was developed to be used in the prediction of the flow in the SSME high-pressure turbopump blade passages. The rotor viscous code (RVC) employs a four-step Runge-Kutta scheme to solve the two-dimensional, thin-layer Navier-Stokes equations. The Baldwin-Lomax eddy-viscosity model is used for these turbulent flow calculations. A viable method was developed to use the relative exit conditions from an upstream blade row as the inlet conditions to the next blade row. The blade loading diagrams are compared with the meridional values obtained from an in-house quasithree-dimensional inviscid code. Periodic boundary conditions are imposed on a body-fitted C-grid computed by using the GRAPE GRids about Airfoils using Poisson's Equation (GRAPE) code. Total pressure, total temperature, and flow angle are specified at the inlet. The upstream-running Riemann invariant is extrapolated from the interior. Static pressure is specified at the exit such that mass flow is conserved from blade row to blade row, and the conservative variables are extrapolated from the interior. For viscous flows the noslip condition is imposed at the wall. The normal momentum equation gives the pressure at the wall. The density at the wall is obtained from the wall total temperature.

Jorgenson, P. G. E.↗

Flow structure capturing on overset patched meshes

The present paper describes the application of a simple, robust stable implicit approach to solutions of the conservative equations of gasdynamics on either composite or overset meshes to capture flow structures such as shocks with high resolution in complex geometric domains. Without requiring special flux conservative operators, but rather, interpolating conservative variable data at mesh boundaries, the well posed upwind method provides accurate and oscillation free solutions, even when shocks cross the interior patch boundaries. In three problems with flow complexities that require grid refinement, the paper demonstrates the capability to conveniently carry out for gasdynamics the adaptive refined meshing strategy in overset patches proposed by Berger and Oliger (1984), and it extends this technique to rapidly convergent implicit methods for the Euler and Navier-Stokes equations. The numerical experiments show concretely, in a realistic aerodynamic problem, the savings in mesh points (about an order of magnitude here in two dimensions) for similar accuracy that flow structure aligned adaptive patched meshing affords compared to uniform grid refinement.

Venkatapathy, E.↗

Saturation point structure of marine stratocumulus clouds

An investigation of the microstructure of a Pacific stratocumulus capped boundary layer is presented. A complex structure of three branches, identified using conserved variable diagrams, is found to correspond well to a conceptual model for the unstable, radiatively cooled cloud topped boundary layer. A simple conditional sampling method was used to identify saturation point pairs for ascending and descending branches of the internal boundary layer circulation. Results indicate a primary circulation scale of 5 km and provide a reasonable cloud top entrainment rate of 1 cm/s.

Boers, Reinout↗

A stratocumulus thermodynamic analysis: July 5 case study

On July 5 (NCAR Electra flight 4, Mission 186-G) the Electra flew a single aircraft mission which consisted of cross and along-wind legs at 6 different altitudes between 10:43 to 16:00 PDT (17:43 to 23:00 GMT). The leg length was kept short (8 to 10 minutes) to permit maximum vertical resolution, and there were 8 soundings. Observer notes report a thin, solid stratocumulus cloud deck which gradually became more broken in the afternoon. Winds were from the north at 10 to 13/ms throughout the night. Sea surface temperature measurements and conservative variable analyses for several of the July 5 legs are presented. These results are preliminary to a study of the thermodynamic budget on July 5; they indicate that: the sea surface temperature dropped more than 1 K (from 17.3 to 15.9 C) over the course of the flight (18:01 and at 21:51 GMT); mixing lines for each of the horizontal subcloud legs show the effect of a strong north-south gradient in SST; and there is a clear demarcation over a transition of 5 to 10 km between air to the south and cooler, moister air to the north. The FSSP measurements indicate there are small clouds/scud 250 m below cloud base on the cold northern side of this transition. The transition is seen in the saturation point diagrams at 984, 959, and 946 mb. There is no corresponding change in the horizontal wind across the transition regions.

Austin, Philip↗

A unified view of convective transports by stratocumulus clouds, shallow cumulus clouds, and deep convection

A bulk planetary boundary layer (PBL) model was developed with a simple internal vertical structure and a simple second-order closure, designed for use as a PBL parameterization in a large-scale model. The model allows the mean fields to vary with height within the PBL, and so must address the vertical profiles of the turbulent fluxes, going beyond the usual mixed-layer assumption that the fluxes of conservative variables are linear with height. This is accomplished using the same convective mass flux approach that has also been used in cumulus parameterizations. The purpose is to show that such a mass flux model can include, in a single framework, the compensating subsidence concept, downgradient mixing, and well-mixed layers.

Randall, David A.↗

Cold-air outbreak during GALE - Lidar observations and modeling of boundary layer dynamics

Two cold-air outbreaks were studied during the Genesis of Atlantic Lows Experiment. A lidar system was operated to observe the boundary layer evolution and the development of clouds. On the first day (January 30, 1986) boundary layer rise was less than 50 percent of the value for the second day (March 2, 1986). On the first day only a thin broken cloud cover formed, while on the second day a thick solid cloud deck formed - although the average moisture content was 60 percent of that on the first day. A trajectory slab model was employed to simulate the evolution of the layer over the ocean near the east Atlantic shore. The model allows for vertical gradients in conservative variables under neutrally buoyant conditions. The primary effect of these assumptions, which are based on observed thermodynamic profiles, is to reduce cloudiness to be more in line with observations. Boundary-layer depth was reasonably well predicted as was sensible and latent heat flux.

Boers, Reinout↗

FIFE atmospheric boundary layer budget methods

The budget methods and the mixed layer model employed to analyze the aircraft data from the First ISLSCP Field Experiment (FIFE) are described. Vector budgets for the mixed layer are discussed on conserved variable diagrams. Theoretical solutions are presented for the critical surface Bowen ratio that produces no boundary layer moistening or equivalent potential temperature rise as a function of the Bowen ratio at the inversion.

Betts, A. K.↗

A multi-dimensional kinetic-based upwind solver for the Euler equations

A multidimensional kinetic fluctuation-splitting scheme has been developed for the Euler equations. The scheme is based on an N-scheme discretization of the Boltzmann equation at the kinetic level for triangulated Cartesian meshes with a diagonal-adaptive strategy. The resulting Euler scheme is a cell-vertex fluctuation-splitting scheme where fluctuations in the conserved-variable vector Q are obtained as moments of the fluctuation in the Maxwellian velocity distribution function at the kinetic level. Encouraging preliminary results have been obtained for perfect gases on Cartesian meshes with first-order spatial accuracy. The present approach represents an improvement to the well-established dimensionally-split upwind schemes.

Eppard, W. M.↗

Comparison of Mixing Calculations for Reacting and Non-Reacting Flows in a Cylindrical Duct

A production 3-D elliptic flow code has been used to calculate non-reacting and reacting flow fields in an experimental mixing section relevant to a rich burn/quick mix/lean burn (RQL) combustion system. A number of test cases have been run to assess the effects of the variation in the number of orifices, mass flow ratio, and rich-zone equivalence ratio on the flow field and mixing rates. The calculated normalized temperature profiles for the non-reacting flow field agree qualitatively well with the normalized conserved variable isopleths for the reacting flow field indicating that non-reacting mixing experiments are appropriate for screening and ranking potential rapid mixing concepts. For a given set of jet momentum-flux ratio, mass flow ratio, and density ratio (J, MR, and DR), the reacting flow calculations show a reduced level of mixing compared to the non-reacting cases. In addition, the rich-zone equivalence ratio has noticeable effect on the mixing flow characteristics for reacting flows.

Oechsle, V. L.↗

A Second Law Based Unstructured Finite Volume Procedure for Generalized Flow Simulation

An unstructured finite volume procedure has been developed for steady and transient thermo-fluid dynamic analysis of fluid systems and components. The procedure is applicable for a flow network consisting of pipes and various fittings where flow is assumed to be one dimensional. It can also be used to simulate flow in a component by modeling a multi-dimensional flow using the same numerical scheme. The flow domain is discretized into a number of interconnected control volumes located arbitrarily in space. The conservation equations for each control volume account for the transport of mass, momentum and entropy from the neighboring control volumes. In addition, they also include the sources of each conserved variable and time dependent terms. The source term of entropy equation contains entropy generation due to heat transfer and fluid friction. Thermodynamic properties are computed from the equation of state of a real fluid. The system of equations is solved by a hybrid numerical method which is a combination of simultaneous Newton-Raphson and successive substitution schemes. The paper also describes the application and verification of the procedure by comparing its predictions with the analytical and numerical solution of several benchmark problems.

Majumdar, Alok↗