Results for “Primitive equations” — Page 20

NASA NTRS1984

Use of shape-preserving interpolation methods in surface modeling

In many large-scale scientific computations, it is necessary to use surface models based on information provided at only a finite number of points (rather than determined everywhere via an analytic formula). As an example, an equation of state (EOS) table may provide values of pressure as a function of temperature and density for a particular material. These values, while known quite accurately, are typically known only on a rectangular (but generally quite nonuniform) mesh in (T,d)-space. Thus interpolation methods are necessary to completely determine the EOS surface. The most primitive EOS interpolation scheme is bilinear interpolation. This has the advantages of depending only on local information, so that changes in data remote from a mesh element have no effect on the surface over the element, and of preserving shape information, such as monotonicity. Most scientific calculations, however, require greater smoothness. Standard higher-order interpolation schemes, such as Coons patches or bicubic splines, while providing the requisite smoothness, tend to produce surfaces that are not physically reasonable. This means that the interpolant may have bumps or wiggles that are not supported by the data. The mathematical quantification of ideas such as physically reasonable and visually pleasing is examined.

Ftitsch, F. N.↗

NASA NTRS1982

A Computer Code for Swirling Turbulent Axisymmetric Recirculating Flows in Practical Isothermal Combustor Geometries

A primitive pressure-velocity variable finite difference computer code was developed to predict swirling recirculating inert turbulent flows in axisymmetric combustors in general, and for application to a specific idealized combustion chamber with sudden or gradual expansion. The technique involves a staggered grid system for axial and radial velocities, a line relaxation procedure for efficient solution of the equations, a two-equation k-epsilon turbulence model, a stairstep boundary representation of the expansion flow, and realistic accommodation of swirl effects. A user's manual, dealing with the computational problem, showing how the mathematical basis and computational scheme may be translated into a computer program is presented. A flow chart, FORTRAN IV listing, notes about various subroutines and a user's guide are supplied as an aid to prospective users of the code.

Lilley, D. G.↗

NASA NTRS1984

Space structure vibration modes: How many exist? Which ones are important?

This report attempts to shed some light on the two issues raised in the title, namely, how many vibration modes does a real structure have, and which of these modes are important? The surprise-free answers to these two questions are, respectively, an infinite number and the first several modes. The author argues that the absurd subspace (all but the first billion modes) is not a strength of continuum modeling, but, in fact, a weakness. Partial differential equations are not real structures, only mathematical models. This note also explains (1) that the PDE model and the finite element model are, in fact, the same model, the latter being a numerical method for dealing with the former, (2) that modes may be selected on dynamical grounds other than frequency alone, and (3) that long slender rods are useful as primitive cases but dangerous to extrapolate from.

Hughes, P. C.↗

NASA NTRS1985

A Numerical Method of Calculating Propeller Noise Including Acoustic Nonlinear Effects

Using the transonic flow fields(s) generated by the NASPROP-E computer code for an eight blade SR3-series propeller, a theoretical method is investigated to calculate the total noise values and frequency content in the acoustic near and far field without using the Ffowcs Williams - Hawkings equation. The flow field is numerically generated using an implicit three dimensional Euler equation solver in weak conservation law form. Numerical damping is required by the differencing method for stability in three dimensions, and the influence of the damping on the calculated acoustic values is investigated. The acoustic near field is solved by integrating with respect to time the pressure oscillations induced at a stationary observer location. The acoustic far field is calculated from the near field primitive variables as generated by NASPROP-E computer code using a method involving a perturbation velocity potential as suggested by Hawkings in the calculation of the acoustic pressure time-history at a specified far field observed location. the methodologies described are valid for calculating total noise levels and are applicable to any propeller geometry for which a flow field solution is available.

Korkan, K. D.↗

NASA NTRS1987

MPF: A portable message passing facility for shared memory multiprocessors

The design, implementation, and performance evaluation of a message passing facility (MPF) for shared memory multiprocessors are presented. The MPF is based on a message passing model conceptually similar to conversations. Participants (parallel processors) can enter or leave a conversation at any time. The message passing primitives for this model are implemented as a portable library of C function calls. The MPF is currently operational on a Sequent Balance 21000, and several parallel applications were developed and tested. Several simple benchmark programs are presented to establish interprocess communication performance for common patterns of interprocess communication. Finally, performance figures are presented for two parallel applications, linear systems solution, and iterative solution of partial differential equations.

Malony, Allen D.↗

NASA NTRS1989

Pb isotopes in anorthositic breccias 67075 and 62237 - A search for primitive lunar lead

The lunar Pb isotopic evolution is explored based on an investigation of the U-Th-Pb systematics of the primitive Apollo 16 anorthositic breccias 76075 and 62237. The isotope data for acid-leached residues from 67075 lie within or on a triangle formed by plagioclase, opaques, and a pyroxen-olivine mixture, and the isotope data for 62237 lie within or on a triangle formed by plagioclase, olivine, and pyroxene. The results suggest that the moon had high first-stage values for the Pb evolution equations from the onset.

Premo, W. R.↗

NASA NTRS1975

Numerical solution of the Navier-Stokes equations for arbitrary two-dimensional multi-element airfoils

Abstracts are presented on a method of numerical solution of the Navier-Stokes equation for the flow about arbitrary airfoils, using a numerically generated curvilinear coordinate system having a coordinate line coincident with the body contour. Results of continuing research are reported and include: application of the Navier-Stokes solution in the vorticity-stream function formulation to a number of single airfoils at Reynolds numbers up to 2000; programming of the Navier-Stokes solution for multiple airfoils in the primitive variable formulation; testing of the potential flow solution of multiple bodies; and development of a generalized coordinate system program.

Thompson, J. F.↗

NASA NTRS1987

Initial planetary base construction techniques and machine implementation

Conceptual designs of (1) initial planetary base structures, and (2) an unmanned machine to perform the construction of these structures using materials local to the planet are presented. Rock melting is suggested as a possible technique to be used by the machine in fabricating roads, platforms, and interlocking bricks. Identification of problem areas in machine design and materials processing is accomplished. The feasibility of the designs is contingent upon favorable results of an analysis of the engineering behavior of the product materials. The analysis requires knowledge of several parameters for solution of the constitutive equations of the theory of elasticity. An initial collection of these parameters is presented which helps to define research needed to perform a realistic feasibility study. A qualitative approach to estimating power and mass lift requirements for the proposed machine is used which employs specifications of currently available equipment. An initial, unmanned mission scenario is discussed with emphasis on identifying uncompleted tasks and suggesting design considerations for vehicles and primitive structures which use the products of the machine processing.

Crockford, William W.↗

NASA NTRS2000

Evaluation of the Lattice-Boltzmann Equation Solver PowerFLOW for Aerodynamic Applications

A careful comparison of the performance of a commercially available Lattice-Boltzmann Equation solver (Power-FLOW) was made with a conventional, block-structured computational fluid-dynamics code (CFL3D) for the flow over a two-dimensional NACA-0012 airfoil. The results suggest that the version of PowerFLOW used in the investigation produced solutions with large errors in the computed flow field; these errors are attributed to inadequate resolution of the boundary layer for reasons related to grid resolution and primitive turbulence modeling. The requirement of square grid cells in the PowerFLOW calculations limited the number of points that could be used to span the boundary layer on the wing and still keep the computation size small enough to fit on the available computers. Although not discussed in detail, disappointing results were also obtained with PowerFLOW for a cavity flow and for the flow around a generic helicopter configuration.

Lockard, David P.↗

NASA NTRS1982

A basic code for the prediction of transient three-dimensional turbulent flowfields

A primitive pressure-velocity finite difference code has been developed to predict transient three-dimensional turbulent flow. The code is a simplified yet effective prediction procedure for use by persons with little experience in computational fluid dynamics, and into which user-oriented complexities can easily be added. The method is based on the transient two-dimensional Los Alamos SOLA prediction technique for laminar flows. Turbulence is simulated by means of a two-equation turbulence model; species diffusion and buoyancy are also included. Two applications of the code are presented to local destratification near the release structure of a reservoir and to the deflection of a jet entering normally into a uniform cross-flow. Predicted results exhibit good agreement with experimental data, showing that a useful characterization of fully three-dimensional flows is now available.

Busnaina, A. A.↗

NASA NTRS2008

Computational Analyses in Support of Sub-scale Diffuser Testing for the A-3 Facility: Steady Predictions - Part 1

Simulation technology can play an important role in rocket engine test facility design and development by assessing risks, providing analysis of dynamic pressure and thermal loads, identifying failure modes and predicting anomalous behavior of critical systems. Advanced numerical tools assume greater significance in supporting testing and design of high altitude testing facilities and plume induced testing environments of high thrust engines because of the greater inter-dependence and synergy in the functioning of the different sub-systems. This is especially true for facilities such as the proposed A-3 facility at NASA SSC because of a challenging operating envelope linked to variable throttle conditions at relatively low chamber pressures. Facility designs in this case will require a complex network of diffuser ducts, steam ejector trains, fast operating valves, cooling water systems and flow diverters that need to be characterized for steady state performance. In this paper, we will demonstrate with the use of CFD analyses s advanced capability to evaluate supersonic diffuser and steam ejector performance in a sub-scale A-3 facility at NASA Stennis Space Center (SSC) where extensive testing was performed. Furthermore, the focus in this paper relates to modeling of critical sub-systems and components used in facilities such as the A-3 facility. The work here will address deficiencies in empirical models and current CFD analyses that are used for design of supersonic diffusers/turning vanes/ejectors as well as analyses for confined plumes and venting processes. The primary areas that will be addressed are: (1) supersonic diffuser performance including analyses of thermal loads (2) accurate shock capturing in the diffuser duct; (3) effect of turning duct on the performance of the facility (4) prediction of mass flow rates and performance classification for steam ejectors (5) comparisons with test data from sub-scale diffuser testing and assessment of confidence levels in CFD based flowpath modeling of the facility. The analyses tools used here expand on the multi-element unstructured CFD which has been tailored and validated for impingement dynamics of dry plumes, complex valve/feed systems, and high pressure propellant delivery systems used in engine and component test stands at NASA SSC. The analyses performed in the evaluation of the sub-scale diffuser facility explored several important factors that influence modeling and understanding of facility operation such as (a) importance of modeling the facility with Real Gas approximation, (b) approximating the cluster of steam ejector nozzles as a single annular nozzle, (c) existence of mixed subsonic/supersonic flow downstream of the turning duct, and (d) inadequacy of two-equation turbulence models in predicting the correct pressurization in the turning duct and expansion of the second stage steam ejectors. The procedure used for modeling the facility was as follows: (i) The engine, test cell and first stage ejectors were simulated with an axisymmetric approximation (ii) the turning duct, second stage ejectors and the piping downstream of the second stage ejectors were analyzed with a three-dimensional simulation utilizing a half-plane symmetry approximation. The solution i.e. primitive variables such as pressure, velocity components, temperature and turbulence quantities were passed from the first computational domain and specified as a supersonic boundary condition for the second simulation. (iii) The third domain comprised of the exit diffuser and the region in the vicinity of the facility (primary included to get the correct shock structure at the exit of the facility and entrainment characteristics). The first set of simulations comprising the engine, test cell and first stage ejectors was carried out both as a turbulent real gas calculation as well as a turbulent perfect gas calculation. A comparison for the two cases (Real Turbulent and Perfect gas turbulent) of the Ma Number distribution and temperature distributions are shown in Figures 1 and 2 respectively. The Mach Number distribution shows small yet distinct differences between the two cases such as locations of shocks/shock reflections and a slightly different impingement point on the wall of the diffuser from the expansion at the exit of the nozzle. Similarly the temperature distribution indicates different flow recirculation patterns in the test cell. Both cases capture all the essential flow phenomena such as the shock-boundary layer interaction, plume expansion, expansion of the first stage ejectors, mixing between the engine plume and the first stage ejector flow and pressurization due to the first stage ejectors. The final paper will discuss thermal loads on the walls of the diffuser and cooling mechanisms investigated.

Allgood, Daniel C.↗

NASA NTRS2010

Computational Analyses in Support of Sub-scale Diffuser Testing for the A-3 Facility: Steady Predictions - Part 1

Simulation technology can play an important role in rocket engine test facility design and development by assessing risks, providing analysis of dynamic pressure and thermal loads, identifying failure modes and predicting anomalous behavior of critical systems. Advanced numerical tools assume greater significance in supporting testing and design of high altitude testing facilities and plume induced testing environments of high thrust engines because of the greater inter-dependence and synergy in the functioning of the different sub-systems. This is especially true for facilities such as the proposed A-3 facility at NASA SSC because of a challenging operating envelope linked to variable throttle conditions at relatively low chamber pressures. Facility designs in this case will require a complex network of diffuser ducts, steam ejector trains, fast operating valves, cooling water systems and flow diverters that need to be characterized for steady state performance. In this paper, we will demonstrate with the use of CFD analyses s advanced capability to evaluate supersonic diffuser and steam ejector performance in a sub-scale A-3 facility at NASA Stennis Space Center (SSC) where extensive testing was performed. Furthermore, the focus in this paper relates to modeling of critical sub-systems and components used in facilities such as the A-3 facility. The work here will address deficiencies in empirical models and current CFD analyses that are used for design of supersonic diffusers/turning vanes/ejectors as well as analyses for confined plumes and venting processes. The primary areas that will be addressed are: (1) supersonic diffuser performance including analyses of thermal loads (2) accurate shock capturing in the diffuser duct; (3) effect of turning duct on the performance of the facility (4) prediction of mass flow rates and performance classification for steam ejectors (5) comparisons with test data from sub-scale diffuser testing and assessment of confidence levels in CFD based flowpath modeling of the facility. The analyses tools used here expand on the multi-element unstructured CFD which has been tailored and validated for impingement dynamics of dry plumes, complex valve/feed systems, and high pressure propellant delivery systems used in engine and component test stands at NASA SSC. The analyses performed in the evaluation of the sub-scale diffuser facility explored several important factors that influence modeling and understanding of facility operation such as (a) importance of modeling the facility with Real Gas approximation, (b) approximating the cluster of steam ejector nozzles as a single annular nozzle, (c) existence of mixed subsonic/supersonic flow downstream of the turning duct, and (d) inadequacy of two-equation turbulence models in predicting the correct pressurization in the turning duct and expansion of the second stage steam ejectors. The procedure used for modeling the facility was as follows: (i) The engine, test cell and first stage ejectors were simulated with an axisymmetric approximation (ii) the turning duct, second stage ejectors and the piping downstream of the second stage ejectors were analyzed with a three-dimensional simulation utilizing a half-plane symmetry approximation. The solution i.e. primitive variables such as pressure, velocity components, temperature and turbulence quantities were passed from the first computational domain and specified as a supersonic boundary condition for the second simulation. (iii) The third domain comprised of the exit diffuser and the region in the vicinity of the facility (primary included to get the correct shock structure at the exit of the facility and entrainment characteristics). The first set of simulations comprising the engine, test cell and first stage ejectors was carried out both as a turbulent real gas calculation as well as a turbulent perfect gas calculation. A comparison for the two cases (Real Turbulent and Perfect gas turbulent) of the Ma Number distribution and temperature distributions are shown in Figures 1 and 2 respectively.

Allgood, Daniel C.↗

NASA NTRS1972

Origin, evolution and present thermal state of the moon.

The relative absence of lunar volcanism in the last 3,000,000,000 years and the Apollo 15 heat flow measurement suggest that present-day temperatures in the moon are approximately steady-state to depths of about 100 km. An exponential distribution of heat sources with depth may then be scaled by equating the surface heat flow to the integrated heat production of this exterior shell. Presumed present-day interior temperatures, as well as the present-day surface heat flow of about 30 erg/sq cm-sec, may be obtained with an initial temperature roughly corresponding to the Apollo 11 basalt solidus and exponential scaling of heat sources. The concentrations of U for an originally homogeneous moon are estimated to be about 0.009 micrograms per gram, close to that measured for eucrites and inferred for primitive inclusions of the Allende meteorite. The estimated homogeneous concentrations of U, the chemistry of the lunar surface material and inferences to modest depth, and the short accretion time of the moon necessary to provide large-scale differentiation at 4.6 AE suggest that the moon had its origin in the rapid accretion of compounds first condensing from the protoplanetary nebula.

Hanks, T. C.↗

NASA NTRS1999

Molecular evolution inferred from small subunit rRNA sequences: what does it tell us about phylogenetic relationships and taxonomy of the parabasalids?

The Parabasala are a primitive group of protists divided into two classes: the trichomonads and the hypermastigids. Until recently, phylogeny and taxonomy of parabasalids were mainly based on the comparative analysis of morphological characters primarily linked to the development of their cytoskeleton. Recent use of molecular markers, such as small subunit (SSU) rRNA has led to now insights into the systematics of the Parabasala and other groups of prolists. An updated phylogeny based on SSU rRNA is provided and compared to that inferred from ultrastructural data. The SSU rRNA phylogeny contradicts the dogma equating simple characters with pumitive characters. Hypermastigids, possessing a hyperdeveloped cytoskeleton, exhibit the most basal emergence in the parabasalid lineage. Other observations emerge from the SSU rRNA analysis, such as the secondary loss of some cytoskeleton structures in all representatives of the Monocercomonadidae, the existence of secondarily free living taxa (reversibility of parasitism) and the evidence against the co-evolution of the endobiotic parabasalids and their animal hosts. According to phylogenies based on SSU rRNA, all the trichomonad families are not monophyletic groups, putting into question the validity of current taxonomic assignments. The precise branching order of some taxa remains unclear, but this issue can possibly be addressed by the molecular analysis of additional parabasalids. The goal of such additional analyses would be to propose, in a near future, a revision of the taxonomy of this group of protists that takes into account both molecular and morphological data.

NASA Discipline Exobiology↗

NASA NTRS1989

Automatic determination of fault effects on aircraft functionality

The problem of determining the behavior of physical systems subsequent to the occurrence of malfunctions is discussed. It is established that while it was reasonable to assume that the most important fault behavior modes of primitive components and simple subsystems could be known and predicted, interactions within composite systems reached levels of complexity that precluded the use of traditional rule-based expert system techniques. Reasoning from first principles, i.e., on the basis of causal models of the physical system, was required. The first question that arises is, of course, how the causal information required for such reasoning should be represented. The bond graphs presented here occupy a position intermediate between qualitative and quantitative models, allowing the automatic derivation of Kuipers-like qualitative constraint models as well as state equations. Their most salient feature, however, is that entities corresponding to components and interactions in the physical system are explicitly represented in the bond graph model, thus permitting systematic model updates to reflect malfunctions. Researchers show how this is done, as well as presenting a number of techniques for obtaining qualitative information from the state equations derivable from bond graph models. One insight is the fact that one of the most important advantages of the bond graph ontology is the highly systematic approach to model construction it imposes on the modeler, who is forced to classify the relevant physical entities into a small number of categories, and to look for two highly specific types of interactions among them. The systematic nature of bond graph model construction facilitates the process to the point where the guidelines are sufficiently specific to be followed by modelers who are not domain experts. As a result, models of a given system constructed by different modelers will have extensive similarities. Researchers conclude by pointing out that the ease of updating bond graph models to reflect malfunctions is a manifestation of the systematic nature of bond graph construction, and the regularity of the relationship between bond graph models and physical reality.

Feyock, Stefan↗

NASA NTRS2012

Simulation of Oxygen Disintegration and Mixing With Hydrogen or Helium at Supercritical Pressure

The simulation of high-pressure turbulent flows, where the pressure, p, is larger than the critical value, p(sub c), for the species under consideration, is relevant to a wide array of propulsion systems, e.g. gas turbine, diesel, and liquid rocket engines. Most turbulence models, however, have been developed for atmospheric-p turbulent flows. The difference between atmospheric-p and supercritical-p turbulence is that, in the former situation, the coupling between dynamics and thermodynamics is moderate to negligible, but for the latter it is very significant, and can dominate the flow characteristics. The reason for this stems from the mathematical form of the equation of state (EOS), which is the perfect-gas EOS in the former case, and the real-gas EOS in the latter case. For flows at supercritical pressure, p, the large eddy simulation (LES) equations consist of the differential conservation equations coupled with a real-gas EOS. The equations use transport properties that depend on the thermodynamic variables. Compared to previous LES models, the differential equations contain not only the subgrid scale (SGS) fluxes, but also new SGS terms, each denoted as a correction. These additional terms, typically assumed null for atmospheric pressure flows, stem from filtering the differential governing equations, and represent differences between a filtered term and the same term computed as a function of the filtered flow field. In particular, the energy equation contains a heat-flux correction (q-correction) that is the difference between the filtered divergence of the heat flux and the divergence of the heat flux computed as a function of the filtered flow field. In a previous study, there was only partial success in modeling the q-correction term, but in this innovation, success has been achieved by using a different modeling approach. This analysis, based on a temporal mixing layer Direct Numerical Simulation database, shows that the focus in modeling the q-correction should be on reconstructing the primitive variable gradients rather than their coefficients, and proposes the approximate deconvolution model (ADM) as an effective means of flow field reconstruction for LES heat flux calculation. Further, results for a study conducted for temporal mixing layers initially containing oxygen in the lower stream, and hydrogen or helium in the upper stream, show that, for any LES, including SGS-flux models (constant-coefficient Gradient or Scale-Similarity models, dynamic-coefficient Smagorinsky/Yoshizawa or mixed Smagorinsky/Yoshizawa/Gradient models), the inclusion of the q-correction in the LES leads to the theoretical maximum reduction of the SGS heat-flux difference. The remaining error in modeling this new subgrid term is thus irreducible.

Bellan, Josette↗

NASA NTRS2013

Circular-Orbit Maintenance Strategies for Primitive Body Orbiters

For missions to smaller primitive bodies, solar radiation pressure (SRP) is a significant perturbation to Keplerian dynamics. For most orbits, SRP drives large oscillations in orbit eccentricity, which leads to large perturbations from the irregular gravity field at periapsis. Ultimately, chaotic motion results that often escapes or impacts that body. This paper presents an orbit maintenance strategy to keep the orbit eccentricity small, thus avoiding the destabilizing secondary interaction with the gravity field. An estimate of the frequency and magnitude of the required maneuvers as a function of the orbit and body parameters is derived from the analytic perturbation equations.

comet↗

NASA NTRS1988

Developing fluid flow in a curved duct of square cross-section and its fully developed dual solutions

Developing fluid flow in a curved duct of square cross-section is studied numerically by a factored ADI finite-difference method on a staggered grid. A central-difference scheme with primitive variables is used inside the computational domain to reduce numerical diffusion. Two Reynolds numbers, 574 and 790, based upon a bulk velocity and hydraulic diameter are chosen for curvature ratios of 1/6.45 and 1/2.3, respectively. It is found that the secondary flow is far more complicated than expected, with the appearance of at least two pairs of vortices. Main-flow separation is also observed for the higher curvature ratio. Furthermore, it is observed that the flow develops into two quite different states downstream, depending upon the inlet conditions. Solutions of the fully developed Navier-Stokes equations is shown to be not unique beyond a certain critical Reynolds number. Developing flow seems to evolve into the fully developed state along a particular branch into which the fully developed solution bifurcates.

Soh, W. Y.↗