Search NASA⌕ Search

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.

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

At least 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Implicit flux-split Euler schemes for unsteady aerodynamic analysis involving unstructured dynamic meshes

Improved algorithm for the solution of the time-dependent Euler equations are presented for unsteady 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 an implicit time-integration scheme using a Gauss-Seidel relaxation procedure which is computationally efficient for either steady or unsteady flow problems. For example, 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 and unsteady flow results are presented for the NACA 0012 airfoil to demonstrate applications of the new Euler solvers. The unsteady results were obtained for the airfoil pitching harmonically about the quarter chord. The resulting instantaneous pressure distributions and lift and moment coefficients during a cycle of motion compare well with experimental data. A description of the Euler solvers is presented along with results and comparisons which assess the capability.

Batina, John T.↗
NASA NTRS1989

Algorithm For Hypersonic Flow In Chemical Equilibrium

Implicit, finite-difference, shock-capturing algorithm calculates inviscid, hypersonic flows in chemical equilibrium. Implicit formulation chosen because overcomes limitation on mathematical stability encountered in explicit formulations. For dynamical portion of problem, Euler equations written in conservation-law form in Cartesian coordinate system for two-dimensional or axisymmetric flow. For chemical portion of problem, equilibrium state of gas at each point in computational grid determined by minimizing local Gibbs free energy, subject to local conservation of molecules, atoms, ions, and total enthalpy. Major advantage: resulting algorithm naturally stable and captures strong shocks without help of artificial-dissipation terms to damp out spurious numerical oscillations.

Palmer, Grant↗
NASA NTRS1993

A detailed study of mean-flow solutions for stability analysis of transitional flows

A finite-volume upwind-difference parabolized Navier-Stokes code is utilized to obtain laminar mean-flow solutions at Mach 3.5 on a half-angle cone of 5 deg at an angle-of-attack of 2 deg. A detailed study is conducted on this configuration; the main focus is the velocity profiles in the leeward and windward symmetry planes at various axial locations. Comparisons of the solution profiles are made with both a central-difference code that incorporates scalar and matrix dissipation models and another state-of-the-art upwind-difference finitevolume code. The results obtained emphasize the importance of using matrix dissipation models for schemes that require explicit artificial dissipation. These results also illustrate the accuracy and efficiency of the planeby-plane marching procedure for computing mean-flow solutions for predicting the onset of transition with linear instability.

Ramakrishnan, R.↗
NASA NTRS1988

Gravitational instability in two-phase disks and the origin of the moon

Two-phase disks may be gravitationally unstable at temperatures or surface densities at which a disk composed of either single phase would be highly stable. It is argued that two-phase disks can achieve a marginally unstable state (in addition to a highly unstable state that leads to fragmentation), limited by the ability of the photosphere to radiate the energy dissipated in the disk. A self-consistent prescription for the viscosity induced by the slow instabilities is provided. Two-phase disks are more centrally condensed than single-phase disks, and their secular cooling time may be comparable to their spreading time. A circumterrestrial disk of sufficient mass to form the moon provides a detailed example of all the preceding points. Its stability, structure, and dynamical evolution are investigated, and it is concluded that its spreading time is short (about 100 yr); the moon is formed molten, or partially molten; the moon's initial orbit lies in the earth's equatorial plane; and only a small fraction of the disk mass is lost in a wind, although this may represent a substantial fraction of volatiles. Most of these conclusions are independent of how the disk was formed, e.g., from a giant impact.

Thompson, Christopher↗
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 [1, 2, 3]. 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 [4, 5]. 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. 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, no slip boundary conditions, and the hybrid 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, thrust, ect.). These results are then used by the coefficient form of the mathematics module to determine the complex eigenvalues of the Acoustic Velocity Potential Equation (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 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 [2, 3]. 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, Sean R.↗
NASA NTRS1989

Modeling and control of flexible space platforms with articulated payloads

The first steps in developing a methodology for spacecraft control-structure interaction (CSI) optimization are identification and classification of anticipated missions, and the development of tractable mathematical models in each mission class. A mathematical model of a generic large flexible space platform (LFSP) with multiple independently pointed rigid payloads is considered. The objective is not to develop a general purpose numerical simulation, but rather to develop an analytically tractable mathematical model of such composite systems. The equations of motion for a single payload case are derived, and are linearized about zero steady-state. The resulting model is then extended to include multiple rigid payloads, yielding the desired analytical form. The mathematical models developed clearly show the internal inertial/elastic couplings, and are therefore suitable for analytical and numerical studies. A simple decentralized control law is proposed for fine pointing the payloads and LFSP attitude control, and simulation results are presented for an example problem. The decentralized controller is shown to be adequate for the example problem chosen, but does not, in general, guarantee stability. A centralized dissipative controller is then proposed, requiring a symmetric form of the composite system equations. Such a controller guarantees robust closed loop stability despite unmodeled elastic dynamics and parameter uncertainties.

Graves, Philip C.↗
NASA NTRS2022

Multi-Agent Control Approach to the Stability of Linear Deployable Systems

Collisions and undesired dynamic responses may occur during the expansion phase of deployable elements on a spacecraft. This problem is more pronounced in large and complex deployable elements where characterization of transient dynamics using testing and simulation can be challenging. This paper introduces a multi-agent control theory approach to address the problem of collisions and transient behavior in deployable lattice systems. We introduce the graph theoretic and relative error vector-based state-space structures to represent the geometry of a lattice-based deploy system. We show that under the assumption of almost-strict dissipativity, the closed-loop feedback gains of the relative error dynamics are the spring and damper constants in the structure. Moreover, we introduce the control perspective of shaping transient dynamics using eigenstructure assignment and balanced coordinate transformation.

Pablos, Juan L. de↗
DOE OSTI2022

CeO 2 Doping of Hf 0.5 Zr 0.5 O 2 Thin Films for High Endurance Ferroelectric Memories

Ferroelectric switching is demonstrated in CeO 2 -doped Hf 0.5 Zr 0.5 O 2 (HZCO) thin films with application in back-end-of-line compatible embedded memories. At low cerium oxide doping concentrations (2.0–5.6 mol%), the ferroelectric orthorhombic phase is stabilized after annealing at temperatures below 400 °C. HZCO ferroelectrics show reliable switching characteristics beyond 10 11 cycles in TiN/HZCO/TiN capacitors, several orders of magnitude greater than identically processed Hf 0.5 Zr 0.5 O 2 (HZO) capacitors, without sacrificing polarization and retention. Internal photoemission and photoconductivity experiments show that CeO 2 -doping introduces in-gap states in HZCO that are nearly aligned with TiN Fermi level, facilitating electron injection through these states. Furthermore, the enhanced average bulk conduction, which may lead to more uniform thermal dissipation in the HZCO films, delays irreversible degradation via breakdown that leads to device failure after repeated programming cycles.

36 MATERIALS SCIENCE↗
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.↗
DOE OSTI2020

Spectroscopic Studies of Charge-Transfer Character and Photoresponses of F 4 TCNQ-Based Donor–Acceptor Complexes

F 4 TCNQ (2,3,5,6-tetrafluoro-7,7,8,8-tetracyanoquinodimethane) is used widely as a hole-doping agent in photoresponsive organic semiconducting materials, yet relatively little is known about the photoresponses of the F 4 TCNQ· - anion generated via doping. Furthermore, there is still relatively little systematic exploration of how the properties of the local material or chemical environment impacts the driving force for generating these charge-transfer complexes. Here we present spectroscopic and photophysical studies of F 4 TCNQ in charge-transfer complexes (CTCs) with the electron donor N,N’-Diphenyl-N-N’-di-p-tolylbenzene-1,4-diamine (MPDA) both in dichloroethane solution and polystyrene matrices. Integer charge transfer (ICT) between donor and acceptor occurs readily in dichloroethane solvent to form F 4 TCNQ· - :MPDA + CTCs, due to a ~150 mV difference in MPDA + /MPDA and F 4 TCNQ/F 4 TCNQ· - reduction potentials. Ultrafast spectroscopic studies of the CTC as well as electrochemically generated F 4 TCNQ· - and MDPA + reveal that the photoresponses of these CTCs are dominated by that of the dopant anion, including rapid deactivation (800 fs) after excitation to the anion D 1 excited state, followed by slower (~10 ps) vibrational cooling in the anion D 0 state. Excitation to the higher-lying D 2 state results in a rapid relaxation to the D 1 state, in contrast to direct D 2 -D 0 relaxation previously observed for F 4 TCNQ· - in gas phase. CTCs embedded in polystyrene (PS) matrices are observed to lose their integer charge-transfer character upon evaporation of solvent, as evidenced by changes to electronic and vibrational absorption features associated with F 4 TCNQ· - . This change is attributed to the loss of solvent stabilization of the ion pair formed through the charge-transfer reaction. Ultrafast spectral measurements reveal that the photoresponses of the partial charge-transfer (PCT) species embedded in PS are still highly similar to those of the ICT species and unlike that of neutral F 4 TCNQ, implying the electronic properties of the PCT state are likewise dominated by properties of the reduced acceptor molecule. We conclude that excitation of ICT or PCT states introduces optical losses for photoresponses of doped organic semiconductor materials due to the large anion absorption cross section and its rapid, dissipative deactivation dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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.↗