SEARCH · Search NASA

Results for “dissipatively stabilized states”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3

All sources NASA NTRS DOE OSTI

NASA NTRS1999

Entropy Splitting and Numerical Dissipation

A rigorous stability estimate for arbitrary order of accuracy of spatial central difference schemes for initial-boundary value problems of nonlinear symmetrizable systems of hyperbolic conservation laws was established recently by Olsson and Oliger (1994) and Olsson (1995) and was applied to the two-dimensional compressible Euler equations for a perfect gas by Gerritsen and Olsson (1996) and Gerritsen (1996). The basic building block in developing the stability estimate is a generalized energy approach based on a special splitting of the flux derivative via a convex entropy function and certain homogeneous properties. Due to some of the unique properties of the compressible Euler equations for a perfect gas, the splitting resulted in the sum of a conservative portion and a non-conservative portion of the flux derivative. hereafter referred to as the "Entropy Splitting." There are several potential desirable attributes and side benefits of the entropy splitting for the compressible Euler equations that were not fully explored in Gerritsen and Olsson. The paper has several objectives. The first is to investigate the choice of the arbitrary parameter that determines the amount of splitting and its dependence on the type of physics of current interest to computational fluid dynamics. The second is to investigate in what manner the splitting affects the nonlinear stability of the central schemes for long time integrations of unsteady flows such as in nonlinear aeroacoustics and turbulence dynamics. If numerical dissipation indeed is needed to stabilize the central scheme, can the splitting help minimize the numerical dissipation compared to its un-split cousin? Extensive numerical study on the vortex preservation capability of the splitting in conjunction with central schemes for long time integrations will be presented. The third is to study the effect of the non-conservative proportion of splitting in obtaining the correct shock location for high speed complex shock-turbulence interactions. The fourth is to determine if this method can be extended to other physical equations of state and other evolutionary equation sets. If numerical dissipation is needed, the Yee, Sandham, and Djomehri (1999) numerical dissipation is employed. The Yee et al. schemes fit in the Olsson and Oliger framework.

Yee, H. C.↗

NASA NTRS1985

Effect of Scattering on Instabilities in Line-driven Stellar Winds

Line driven O-B stellar winds are unstable to perturbations of short spatial wavelength and the growth rates for such instabilities are very rapid. The nature of the nonlinear development of this instability is unknown, but might possibly be one of blobs of gas driven through ambient gas or a quasiregular train of outward moving shocks. In either case the resulting dissipation of mechanical energy might explain the observed anomalous heating n O-B stars as evidenced by their X-ray emission and high ionization state. It also might explain the observed fine structure of the absorption lines and their time variability. The driving due to the absorption of the stellar continuum flux was considered and the effects of the diffuse, scattered radiation field were neglected. It was shown that under certain special conditions the effect of scattering could reduce the instability growth rate to zero. A stability analysis that includes scattering, but that uses the more physically realistic assumption of complete redistribution instead of coherent scattering, and that includes the effects of transverse velocity gradients, which become important as the flow moves away from the stellar surface is presented. It is found that the instability is eliminated right at the base of the wind, but that as the flow moves outward the instability rate rapidly becomes equal to a substantial fraction of the calculated value, the fraction asymptotically reaching 80% at large radii. Since this still implies many e folds in a characteristic outflow time, the primary conclusion that these winds are highly unstable is unchanged.

Owocki, S. P.↗

NASA NTRS2015

Solid Rocket Motor Combustion Instability Modeling in COMSOL Multiphysics

Combustion instability modeling of Solid Rocket Motors (SRM) remains a topic of active research. Many rockets display violent fluctuations in pressure, velocity, and temperature originating from the complex interactions between the combustion process, acoustics, and steady-state gas dynamics. Recent advances in defining the energy transport of disturbances within steady flow-fields have been applied by combustion stability modelers to improve the analysis framework. Employing this more accurate global energy balance requires a higher fidelity model of the SRM flow-field and acoustic mode shapes. The current industry standard analysis tool utilizes a one dimensional analysis of the time dependent fluid dynamics along with a quasi-three dimensional propellant grain regression model to determine the SRM ballistics. The code then couples with another application that calculates the eigenvalues of the one dimensional homogenous wave equation. The mean flow parameters and acoustic normal modes are coupled to evaluate the stability theory developed and popularized by Culick. The assumption of a linear, non-dissipative wave in a quiescent fluid remains valid while acoustic amplitudes are small and local gas velocities stay below Mach 0.2. The current study employs the COMSOL Multiphysics finite element framework to model the steady flow-field parameters and acoustic normal modes of a generic SRM. This work builds upon previous efforts to verify the use of the acoustic velocity potential equation (AVPE) laid out by Campos. The acoustic velocity potential (psi) describing the acoustic wave motion in the presence of an inhomogeneous steady high-speed flow is defined by, del squared psi - (lambda/c) squared psi - M x [M x del((del)(psi))] - 2((lambda)(M)/c + M x del(M) x (del)(psi) - 2(lambda)(psi)[M x del(1/c)] = 0. with M as the Mach vector, c as the speed of sound, and λ as the complex eigenvalue. The study requires one way coupling of the CFD High Mach Number Flow (HMNF) and mathematics module. The HMNF module evaluates the gas flow inside of a SRM using St. Robert's law to model the solid propellant burn rate, slip boundary conditions, and the supersonic outflow condition. Results from the HMNF model are verified by comparing the pertinent ballistics parameters with the industry standard code outputs (i.e. pressure drop, axial velocity, exit velocity). These results are then used by the coefficient form of the mathematics module to determine the complex eigenvalues of the AVPE. The mathematics model is truncated at the nozzle sonic line, where a zero flux boundary condition is self-satisfying. The remaining boundaries are modeled with a zero flux boundary condition, assuming zero acoustic absorption on all surfaces. The one way coupled analysis is perform four times utilizing geometries determined through traditional SRM modeling procedures. The results of the steady-state CFD and AVPE analyses are used to calculate the linear acoustic growth rate as is defined by Flandro and Jacob. In order to verify the process implemented within COMSOL we first employ the Culick theory and compare the results with the industry standard. After the process is verified, the Flandro/Jacob energy balance theory is employed and results displayed.

Fischbach, S. R.↗

NASA NTRS1989

Thermal state of an ice shell on Europa

The thickness of the ice shell presently hypothesized for Europa, which is decoupled from a silicate core by a liquid water layer, is calculated as a function of colatitude and longitude under suitable assumptions for two plausible ice rheology behaviors: that of Maxwell type, and that of generalized flow-law rheology. Due to the dissipation rate's pronounced temperature dependence, virtually all tidal dissipation is found to concentrate in the lowest few kilometers of the shell. While for some parameter choices an insulating regolith that raises the near-surface temperature by more than a few tens of degrees C may stabilize the shell against polar wander, a modest regolith may enhance the probability of such wander's occurrence through the reduction of the shell's retarding friction.

Ojakangas, Gregory W.↗

NASA NTRS1976

Thickness of magnetic structures associated with the earth's bow shock

The thickness of magnetic structures associated with quasi-perpendicular bow shocks has been observed to vary over more than 2 orders of magnitude. In this paper, a criterion is employed which states that the magnetic structure of perpendicular shocks adjusts itself to keep electrostatic waves excited by the relative drift of electrons and ions at marginal stability. The marginal stability boundary and measured plasma parameters are used to calculate the expected thickness of bow-shock magnetic structures. The calculated and measured thicknesses are compared and found to be in good agreement over a range from less than 10 km to greater than 500 km. The conclusion is that the thickness of magnetic structures associated with quasi-perpendicular bow shocks is determined primarily by dissipation due to electron-ion streaming instabilities rather than by dispersive effects.

Morse, D. L.↗

NASA NTRS1990

Three-dimensional flux-split Euler schemes involving unstructured dynamic meshes

Improved algorithms for the solution of the 3-D time dependent Euler equations are presented for aerodynamic analysis involving unstructured dynamic meshes. The improvements were developed recently to the spatial and temporal discretizations used by unstructured grid flow solvers. The spatial discretization involves a flux split approach which is naturally dissipative and captures shock waves sharply with at most one grid point within the shock structure. The temporal discretization involves either an explicit time integration scheme using a multistage Runge-Kutta procedure or an implicit time integration scheme using a Gauss-Seidel relaxation procedure which is computationally efficient for either steady or unsteady flow problems. With the implicit Gauss-Seidel procedure, very large time steps may be used for rapid convergence to steady state, and the step size for unsteady cases may be selected for temporal accuracy rather than for numerical stability. Steady flow results are presented for both the NACA 0012 airfoil and the ONERA M6 wing to demonstrate applications of the new Euler solvers. A description of the Euler solvers is presented along with results and comparisons which assess the capability.

Batina, John T.↗

NASA NTRS1981

Numerical method for gas dynamics combining characteristic and conservation concepts

An efficient implicit numerical method that solves the compressible Navier-Stokes equations in arbitrary curvilinear coordinates by the finite-volume technique is presented. An intrinsically dissipative difference scheme and a fully implicit treatment of boundary conditions, based on characteristic and conservation concepts, are used to improve stability and accuracy. Efficiency is achieved by using a diagonal form of the implicit algorithm and spatially varying time-steps. Comparisons of various schemes and methods are presented for one- and two-dimensional flows, including transonic separated flow past a thick circular-arc airfoil in a channel. The new method is equal to or better than a version of MacCormack's hybrid method in accuracy and it converges to a steady state up to an order of magnitude faster.

Coakley, T. J.↗

NASA NTRS1991

Accuracy of an unstructured-grid upwind-Euler algorithm for the ONERA M6 wing

Improved algorithms for the solution of the three-dimensional, time-dependent Euler equations are presented for aerodynamic analysis involving unstructured dynamic meshes. The improvements have been developed recently to the spatial and temporal discretizations used by unstructured-grid flow solvers. The spatial discretization involves a flux-split approach that is naturally dissipative and captures shock waves sharply with at most one grid point within the shock structure. The temporal discretization involves either an explicit time-integration scheme using a multistage Runge-Kutta procedure or an implicit time-integration scheme using a Gauss-Seidel relaxation procedure, which is computationally efficient for either steady or unsteady flow problems. With the implicit Gauss-Seidel procedure, very large time steps may be used for rapid convergence to steady state, and the step size for unsteady cases may be selected for temporal accuracy rather than for numerical stability. Steady flow results are presented for both the NACA 0012 airfoil and the Office National d'Etudes et de Recherches Aerospatiales M6 wing to demonstrate applications of the new Euler solvers. The paper presents a description of the Euler solvers along with results and comparisons that assess the capability.

Batina, John T.↗

NASA NTRS1989

A critical comparison of turbulence models for homogeneous shear flows in a rotating frame

A variety of turbulence models, including five second-order closure models and four two equation models, are tested for the problem of homogeneous turbulent shear flow in a rotating frame. The model predictions for the time evolution of the turbulent kinetic energy and dissipation rate, as well as those for the equilibrium states, are compared with the results of physical and numerical experiments. Most of the two-equation models predict the same results for all rotation rates (omega/S) in which there is an exponential time growth of the turbulent kinetic energy and dissipation rate. The second-order closures are qualitatively superior since, consistent with physical and numerical experiments, they only predict this type of unstable flow for intermediate rotation rates in the range -0.1 less than or equal to omega/S less than or equal to 1.6. For rotation rates outside this range, there is an exchange of stabilities with a solution whose kinetic energy and dissipation rate decay with time. Although the second-order closures are superior to the two-equation models, there are still problems with the quantitative accuracy of their predictions.

Speziale, Charles G.↗

NASA NTRS1989

A critical comparison of turbulence models for homogeneous shear flows in a rotating frame

Speziale, Charles G.↗

NASA NTRS1990

A critical comparison of turbulence models for homogeneous shear flows in a rotating frame

Speziale, Charles G.↗

NASA NTRS1990

Three-dimensional flux-split Euler schemes involving unstructured dynamic meshes

Improved algorithms for the solution of the three-dimensional time-dependent Euler equations are presented for aerodynamic analysis involving unstructured dynamic meshes. The improvements have been developed recently to the spatial and temporal discretizations used by unstructured grid flow solvers. The spatial discretization involves a flux-split approach which is naturally dissipative and captures shock waves sharply with at most one grid point within the shock structure. The temporal discretization involves either an explicit time-integration scheme using a multi-stage Runge-Kutta procedure or an implicit time-integration scheme using a Gauss-Seidel relaxation procedure which is computationally efficient for either steady or unsteady flow problems. With the implicit Gauss-Seidel procedure, very large time steps may be used for rapid convergence to steady state, and the step size for unsteady cases may be selected for temporal accuracy rather than for numerical stability. Steady flow results are presented for both the NACA 0012 airfoil and the ONERA M6 wing to demonstrate applications of the new Euler solvers. The paper presents a description of the Euler solvers along with results and comparisons which assess the capability.

Batina, John T.↗

NASA NTRS1994

Rapidly rotating neutron stars in general relativity: Realistic equations of state

We construct equilibrium sequences of rotating neutron stars in general relativity. We compare results for 14 nuclear matter equations of state. We determine a number of important physical parameters for such stars, including the maximum mass and maximum spin rate. The stability of the configurations to quasi-radial perturbations is assessed. We employ a numerical scheme particularly well suited to handle rapid rotation and large departures from spherical symmetry. We provide an extensive tabulation of models for future reference. Two classes of evolutionary sequences of fixed baryon rest mass and entropy are explored: normal sequences, which behave very much like Newtonian sequences, and supramassive sequences, which exist for neutron stars solely because of general relativistic effects. Adiabatic dissipation of energy and angular momentum causes a star to evolve in quasi-stationary fashion along an evolutionary sequence. Supramassive sequences have masses exceeding the maximum mass of a nonrotating neutron star. A supramassive star evolves toward eventual catastrophic collapse to a black hole. Prior to collapse, the star actually spins up as it loses angular momentum, an effect that may provide an observable precursor to gravitational collapse to a black hole.

Cook, Gregory B.↗

NASA NTRS1994

Nonlinear Asymptotic Integration Algorithms for One-dimensional Autonomous Dissipative First-order Odes

Nonlinear asymptotic integrators are applied to one-dimensional, nonlinear, autonomous, dissipative, ordinary differential equations. These integrators, including a one-step explicit, a one-step implicit, and a one- and two-step midpoint algorithm, are designed to follow the asymptotic behavior of a system approaching a steady state. The methods require that the differential equation be written in a particular asymptotic form. This is always possible for a one-dimensional equation with a globally asymptotic steady state. In this case, conditions are obtained to guarantee that the implicit algorithms are well defined. Further conditions are determined for the implicit methods to be contractive. These methods are all first order accurate, while under certain conditions the midpoint algorithms may also become second order accurate. The stability of each method is investigated and an estimate of the local error is provided.

Haslach, Henry W., Jr.↗

NASA NTRS2003

Characterization of High-Frequency Excitation of a Wake by Simulation

Insights into the effects of high-frequency forcing on free shear layer evolution are gained through analysis of several direct numerical simulations. High-frequency forcing of a fully turbulent plane wake results in only a weak transient effect. On the other hand, significant changes in the developed turbulent state may result when high-frequency forcing is applied to a transitional wake. The impacts of varying the characteristics of the high-frequency forcing are examined, particularly, the streamwise wavenumber band in which forcing is applied and the initial amplitude of the forcing. The high-frequency excitation is found to increase the dissipation rate of turbulent kinetic energy, to reduce the turbulent kinetic energy production rate, and to reduce the turbulent kinetic energy suppression increases with forcing amplitude once a threshold level has been reached. For a given initial forcing energy, the largest reduction in turbulent kinetic energy density was achieved by forcing wavenumbers that are about two to three times the neutral wavenumber determined from linear stability theory.

Cain, Alan B.↗

NASA NTRSDate not supplied

Tracking Critical Thermal Metrics throughout the Life Cycle of a Large Observatory Thermal Model

Observatory thermal models for large, complex missions, such as the Wide Field InfraRed Survey Telescope (WFIRST) mission, produce an immense amount of data to be processed. Configuration management of the model throughout the project life cycle has mainly focused on which versions of the subsystem models form the current observatory level configuration. However, the results produced by the model are not nearly as well tracked and traceable. Given the various states of design maturity for each of the components in the WFIRST design, an updated component model is nearly ready to be integrated at the next higher level of assembly about every month or two. With each subsystem model delivery, the observatory model needs to remove the old component, integrate the new one, execute the model, and inspect the results. Usually, this inspection focuses primarily on the newly integrated component. Recently, a Metric Tracking Spreadsheet was developed to help provide a “big picture” view of the entire observatory highlighting key parameters critical to mission performance. This spreadsheet helps track impacts on subsystems by updates of other subsystems that were not intuitively obvious. Metrics tracked include: absorbed environmental loading (to determine effectiveness of sunshield), temperatures of critical avionics, electrical dissipations, heater power predictions, stability of critical optics, parasitic heat leaks in cryogenic region, high level heat flows between elements, and model run time. Producing this data for the same operational configuration with each model update has helped produce a trail of data to evaluate the impact of model updates. While the metrics selected are specific for WFIRST, other large, complex observatories could be well served to establish their own metrics early in the project life cycle to track to quickly assess the impact of any subsystem on other subsystems or the overall system itself.

Thermal Desktop↗

NASA NTRS2009

Liquid-Propellant Rocket Engine Throttling: A Comprehensive Review

Liquid-Propellant Rocket Engines (LREs) are capable of on-command variable thrust or thrust modulation, an operability advantage that has been studied intermittently since the late 1930s. Throttleable LREs can be used for planetary entry and descent, space rendezvous, orbital maneuvering including orientation and stabilization in space, and hovering and hazard avoidance during planetary landing. Other applications have included control of aircraft rocket engines, limiting of vehicle acceleration or velocity using retrograde rockets, and ballistic missile defense trajectory control. Throttleable LREs can also continuously follow the most economical thrust curve in a given situation, compared to discrete throttling changes over a few select operating points. The effects of variable thrust on the mechanics and dynamics of an LRE as well as difficulties and issues surrounding the throttling process are important aspects of throttling behavior. This review provides a detailed survey of LRE throttling centered around engines from the United States. Several LRE throttling methods are discussed, including high-pressure-drop systems, dual-injector manifolds, gas injection, multiple chambers, pulse modulation, throat throttling, movable injector components, and hydrodynamically dissipative injectors. Concerns and issues surrounding each method are examined, and the advantages and shortcomings compared.

Casiano, Matthew↗

NASA NTRS2009

A CFD Approach to Modeling Spacecraft Fuel Slosh

Energy dissipation and resonant coupling from sloshing fuel in spacecraft fuel tanks is a problem that occurs in the design of many spacecraft. In the case of a spin stabilized spacecraft, this energy dissipation can cause a growth in the spacecrafts' nutation (wobble) that may lead to disastrous consequences for the mission. Even in non-spinning spacecraft, coupling between the spacecraft or upper stage flight control system and an unanticipated slosh resonance can result in catastrophe. By using a Computational Fluid Dynamics (CFD) solver such as Fluent, a model for this fuel slosh can be created. The accuracy of the model must be tested by comparing its results to an experimental test case. Such a model will allow for the variation of many different parameters such as fluid viscosity and gravitational field, yielding a deeper understanding of spacecraft slosh dynamics. In order to gain a better understanding of the dynamics behind sloshing fluids, the Launch Services Program (LSP) at the NASA Kennedy Space Center (KSC) is interested in finding ways to better model this behavior. Thanks to past research, a state-of-the-art fuel slosh research facility was designed and fabricated at Embry Riddle Aeronautical University (ERAU). This test facility has produced interesting results and a fairly reliable parameter estimation process to predict the necessary values that accurately characterize a mechanical pendulum analog model. The current study at ERAU uses a different approach to model the free surface sloshing of liquid in a spherical tank using Computational Fluid Dynamics (CFD) methods. Using a software package called Fluent, a model was created to simulate the sloshing motion of the propellant. This finite volume program uses a technique called the Volume of Fluid (VOF) method to model the interaction between two fluids [4]. For the case of free surface slosh, the two fluids are the propellant and air. As the fuel sloshes around in the tank, it naturally displaces the air. Using the conservation of mass, momentum, and energy equations, as well as the VOF equations, one can predict the behavior of the sloshing fluid and calculate the forces, pressure gradients, and velocity field for the entire liquid as a function of time.

Marsell, Brandon↗