Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Stability boundary”

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 1,207 records · Page 67

Calculation of two-dimensional inlet flow fields by an implicit method including viscous effects: User's manual

Inlet flow fields for airbreathing missiles are calculated by the adaptation of a two dimensional computational method developed for the flow around airfoils. A supersonic free stream is assumed to allow the forebody calculation to be uncoupled from the inlet calculation. The inlet calculation employs an implicit, time marching finite difference procedure to solve the thin layer Navier-Stokes equations formulated in body fitted coordinates. The mathematical formulation of the problem and the solution algorithm are given. Numerical stability and accuracy as well as the initial and boundary conditions used are discussed. Instructions for program use and operation along with the overall program logic are also given.

Biringen, S.↗

Nonparallel stability of three-dimensional flows

The linear stability of three-dimensional incompressible, isothermal, nonparallel boundary-layer flows has been investigated. The method of multiple scales is used to derive the partial-differential equations that describe the spatial modulations of the amplitude, phase and wavenumber of a disturbance. Group velocities are used to determine the disturbance growth direction. The envelope method is used to calculate the logarithmic amplitude growth rate N. The theory is applied to the flows over a swept-back tapered wing with boundary-layer suction. Results of such analysis for the X-21 wing are discussed. It is found that the nonparallel effects for this wing is substantial.

Padhye, A. R.↗

On the thermal stability of coronal loop plasma

The stability to thermal perturbation of static models of coronal loops is considered including the effects of cool, radiatively stable material at the loop base. The linear stability turns out to be sensitive only to the boundary conditions assumed on the velocity at the loop base. The question of the appropriate boundary conditions is discussed, and it is concluded that the free surface condition (the pressure perturbation vanishes), rather than the rigid wall (the velocity vanishes), is relevant to the solar case. The static models are found to be thermally unstable, with a growth time of the order of the coronal cooking time. The physical implications of these results for the solar corona and transition region are examined.

Antiochos, S. K.↗

Conservative supra-characteristics method for splitting the hyperbolic systems of gasdynamics with computed boundaries for real and perfect gases

Implicit methods developed by Beam and Warming (1978 and Briley and McDonald (1977) make it possible to overcome the hyperbolic stiffness of the conservative compressible Navier-Stokes equations in the fine wall region computational mesh for high Reynold's number flow. Certain difficulties related to the use of these methods could be overcome by employing an approach reported by Roe (1981). In the present investigation Roe's conceptual framework has been adopted for constructing globally conservative finite difference methods. A globally conservative upwind finite difference method (CSCM) consisting of both implicit interior point and boundary point equations is constructed from a new characteristics based flux difference splitting. It is found that the employed upwind eigenvector split scheme which combines fully coupled implicit interior point and boundary point approximations has the desired properties of robust stability and accuracy.

Lombard, C. K.↗

Application of the implicit MacCormack scheme to the parabolized Navier-Stokes equations

MacCormack's implicit finite-difference scheme was used to solve the two-dimensional parabolized Navier-Stokes (PNS) equations. This method for solving the PNS equations does not require the inversion of block tridiagonal systems of algebraic equations and permits the original explicit MacCormack scheme to be employed in those regions where implicit treatment is not needed. The advantages and disadvantages of the present adaptation are discussed in relation to those of the conventional Beam-Warming scheme for a flat plate boundary layer test case. Comparisons are made for accuracy, stability, computer time, computer storage, and ease of implementation. The present method was also applied to a second test case of hypersonic laminar flow over a 15% compression corner. The computed results compare favorably with experiment and a numerical solution of the complete Navier-Stokes equations.

Lawrence, J. L.↗

Dynamics of dendritic pattern formation

This paper is a brief review of recent and ongoing developments in the theory of dendritic pattern formation. A marginal stability hypothesis has achieved some success in predicting dendritic growth rates, but a fully nonlinear dynamic theory of the dendritic growth mechanism has yet to be developed. A possible clue for the construction of such a theory is the existence of a class of one-dimensional models in which a pattern is found to spread into an unstable region by propagating at the slowest stable speed. This mathematical realization of the marginal stability mechanism has led to the development of a boundary layer model which describes localized dynamics of fully two- or three-dimensional solidification fronts and which retains some of the mathematical features of the propagating systems. Preliminary results of numerical calculations using the boundary layer model indicate that it is capable of probing fully nonlinear features of dendritic pattern formation.

Langer, J. S.↗

The effect of gravity modulation on thermosolutal convection in an infinite layer of fluid

The effect of time-periodic vertical gravity modulation on the onset of thermosolutal convection in an infinite horizontal layer with stress-free boundaries is investigated using Floquet theory for the linear stability analysis. Situations for which the fluid layer is stably stratified in either the fingering or diffusive regimes of double-diffusive convection are considered. Results are presented both with and without steady background acceleration. Modulation may stabilize an unstable base solution or destabilize a stable base solution. In addition to synchronous and subharmonic response to the modulation frequency, instability in the double diffusive system can occur via a complex conjugate mode. In the diffusive regime, where oscillatory onset occurs in the unmodulated system, regions of resonant instability occur and exhibit strong coupling with the unmodulated oscillatory frequency. The response to modulation of the fundamental instability of the unmodulated system is described both analytically and numerically; in the double-diffusive system this mode persists under subcritical conditions as a high-frequency lobe.

Saunders, B. V.↗

MHD modeling of magnetotail instability for anisotropic pressure

Results of three-dimensional MHD simulations of magnetotail dynamics with anisotropic pressure are presented. The pressure tensor is assumed to be gyrotropic, satisfying a modified double adiabatic approximation including Ohmic heating. It is found that these constraints tend to stabilize the tail with respect to the resistive tearing instability. Including an increasing level of additional isotropization results in a gradual transition to the fast instability found in an isotropic model, which is discussed for comparison. Possible isotropization mechanisms yielding magnetotail instability are discussed also. The stabilizing anisotropies are found primarily in the plasma sheet boundary region. This fact indicates a possibly important role of isotropization mechanisms operating in this region in destabilizing the magnetotail and initiating a tearing instability without the necessity of changes in fluctuation levels in the neutral sheet.

Hesse, Michael↗

Convection of tin in a Bridgman system. I - Flow characterization by effective diffusivity measurements

An electrochemical titration method was used to investigate the dynamic states in a cylindrical layer of convecting tin. The liquid tin was contained in a cell, with curved boundaries made of quartz and flat boundaries made of a solid state electrolyte - yttria-stabilized zirconia (YSZ). The electrolyte acted as a window through which a trace amount of oxygen could be pumped in or out by the application of a constant voltage. The concentration at the YSZ interface was monitored by operating the electrochemical cell in the galvanic mode. Experimentally determined effective diffusivities of oxygen were compared with the molecular diffusivity. Dynamic states in the convective flow were thus inferred. Temperature measurements were simultaneously made in order to identify the onset of oscillations from a steady convective regime. The experiments were conducted for two different aspect ratios for various imposed temperature gradients and two different orientations with respect to gravity. Transcritical states were identified and comparison to two-dimensional numerical models were made.

Sears, B.↗

Effect of Surface Waviness on Transition in Three-Dimensional Boundary-Layer Flow

The effect of a surface wave on transition in three-dimensional boundary-layer flow over an infinite swept wing was studied. The mean flow computed using interacting boundary-layer theory, and transition was predicted using linear stability theory coupled with the empirical eN method. It was found that decreasing the wave height, sweep angle, or freestream unit Reynolds number, and increasing the freestream Mach number or suction level all stabilized the flow and moved transition onset to downstream locations.

Masad, Jamal A.↗

Flow Visualization of Low Prandtl Number Fluids using Electrochemical Measurements

It is well established that residual flows exist in contained liquid metal processes. In 1-g processing, buoyancy forces often drive these flows and their magnitudes can be substantial. It is also known that residual flows can exist during microgravity processing, and although greatly reduced in magnitude, they can influence the properties of the processed materials. Unfortunately, there are very few techniques to visualize flows in opaque, high temperature liquid metals, and those available are not easily adapted to flight investigation. In this study, a novel technique is developed that uses liquid tin as the model fluid and solid-state electrochemical cells constructed from Yttria-Stabilized Zirconia (YSZ) to establish and measure dissolved oxygen boundary conditions. The melt serves as a common electrode for each of the electrochemical cells in this design, while independent reference electrodes are maintained at the outside surfaces of the electrolyte. By constructing isolated electrochemical cells at various locations along the container walls, oxygen is introduced or extracted by imposing a known electrical potential or passing a given current between the melt and the reference electrode. This programmed titration then establishes a known oxygen concentration boundary condition at the selected electrolyte-melt interface. Using the other cells, the concentration of oxygen at the electrolyte-melt interface is also monitored by measuring the open-circuit potentials developed between the melt and reference electrodes. Thus the electrochemical cells serve to both establish boundary conditions for the passive tracer and sense its path. Rayleigh-Benard convection was used to validate the electrochemical approach to flow visualization. Thus, a numerical characterization of the second critical Rayleigh numbers in liquid tin was conducted for a variety of Cartesian aspect ratios. The extremely low Prandtl number of tin represents the lowest value studied numerically. Additionally, flow field oscillations are visualized and the effect of tilt on convecting systems is quantified. Experimental studies of the effect of convection in liquid tin are presented. Three geometries are studied: (1) double electrochemical cell with vertical concentration gradients; (2) double cell with horizontal concentration gradients; and (3) multiple cells with vertical temperature gradients. The first critical Rayleigh number transition is detected with geometry (1) and it is concluded that current measurements are not as affected by convection as EMF measurements. The system is compared with numerical simulations in geometry (2), and oscillating convection is detected with geometry (3).

Crunkleton, D.↗

The Road to Mach 10: A History of the X-43A Hypersonic Flight Test Program at NASA Dryden...Toward the Future

In terms of technology, the X-43A/Hyper-X represented a singular milestone. After nearly a half century of high hopes, studies, wind tunnel tests, proposals, and canceled projects, a scramjet-powered vehicle had flown. The performance of the engine qualified the scramjet design tools and scaling laws. In turn, the theoretical calculations and ground testing could be used to design more advanced engine concepts. Just as important, both the scramjet and vehicle systems had successfully operated in the variable temperatures and densities of the atmosphere. The X-43A systems were able to maintain the exact flight conditions necessary for the scramjet to operate properly. Control deflections to correct the engine-induced moments were close to pre-flight predictions. When the unexpected occurred, such as when the vehicle pitched up during the cowl opening on the second flight, the control system was sufficiently designed to correct the situation. The airframe and wing structure, the thermal protection material, and the internal conditions of the X-43A performed largely as predicted. The HXLV thermal anomaly during the ascent on the third flight and "the Mach 8 unpleasantness" during the descent indicated that the HXLV and X-43A were not as resilient to aerodynamic heating as expected. The X-43A 's airframe drag and lift both were slightly higher than predicted, but still within preflight uncertainty predictions. The stability and control were as predicted, as was the boundary layer transition. The biggest aerodynamic worry before the flight was the separation of the HXLV and the X- 43A. After all was said and done, this went exactly as predicted, proving that non-symmetrical/high-dynamic pressure stage separations could be performed. This in turn meant that two-stage-to-orbit vehicles employing this technology were feasible. The Hyper-X program also served as a training ground for a new generation of scramjet and hypersonic researchers. This included both NASA and contractor personnel, providing them with experience in ground testing and component development; vehicle design, construction, integration, system checkout, and, ultimately, flight testing and data analysis. Additionally, researchers learned the practical details of running a project within finite budget and time limits, about the ambiguousness of risk assessment, and about the need to spend a significant amount of time and effort dealing with engineering problems, such as those with the FAS, that have nothing to do with the project's research goals. Finally, all those who worked on the X-43A project now know what it is like to spend years transforming an idea into a functional vehicle, only for it to be lost in a matter of seconds. And then to go through years of work to correct the problems, to face the possibility that still more might exist, and finally to savor the triumph of two successful flights. For those who will work on the hypersonic projects that emerge in coming years, these experiences may prove to be the most valuable of all.

Peebles, Curtis↗

Fomalhaut's Debris Disk and Planet: Constraining the Mass of Formalhaut B from Disk Morphology

Following the optical imaging of exoplanet candidate Fomalhaut b (Fom b), we present a numerical model of how Fomalhaut's debris disk is gravitationally shaped by a single interior planet. The model is simple, adaptable to other debris disks, and can be extended to accommodate multiple planets. If Fom b is the dominant perturber of the belt, then to produce the observed disk morphology it must have a mass M(sub pl) < 3M(sub J), an orbital semimajor axis a(sub pl) > 101.5AU, and an orbital eccentricity e(sub pl) = 0.11 - 0.13. These conclusions are independent of Fom b's photometry. To not disrupt the disk, a greater mass for Fom b demands a smaller orbit farther removed from the disk; thus, future astrometric measurement of Fom b's orbit, combined with our model of planet-disk interaction, can be used to determine the mass more precisely. The inner edge of the debris disk at a approximately equals 133AU lies at the periphery of Fom b's chaotic zone, and the mean disk eccentricity of e approximately equals 0.11 is secularly forced by the planet, supporting predictions made prior to the discovery of Fom b. However, previous mass constraints based on disk morphology rely on several oversimplifications. We explain why our constraint is more reliable. It is based on a global model of the disk that is not restricted to the planet's chaotic zone boundary. Moreover, we screen disk parent bodies for dynamical stability over the system age of approximately 100 Myr, and model them separately from their dust grain progeny; the latter's orbits are strongly affected by radiation pressure and their lifetimes are limited to approximately 0.1 Myr by destructive grain-grain collisions. The single planet model predicts that planet and disk orbits be apsidally aligned. Fomalhaut b's nominal space velocity does not bear this out, but the astrometric uncertainties are difficult to quantify. Even if the apsidal misalignment proves real, our calculated upper mass limit of 3 M(sub J) still holds. Parent bodies are evacuated from mean-motion resonances with Fom b; these empty resonances are akin to the Kirkwood gaps opened by Jupiter. The belt contains at least 3M(sub Earth) of solids that are grinding down to dust, their velocity dispersions stirred so strongly by Fom b that collisions are destructive. Such a large mass in solids is consistent with Fom b having formed in situ.

Chiang, E.↗

Structure and Dynamics of Domains in Ferroelectric Nanostructures – Phase-Field Modeling (Final Technical Report)

Domain pattern formation is one of the most common phenomena in nature and is a topic of immense interest in many fields ranging from materials science and physics to chemistry and biology. This DOE sponsored research project explored the basic science concerning the thermodynamic stability of mesoscale polarization domain patterns and their temporal evolution mechanisms during formation and subsequent switching in ferroelectric nanostructures and heterostructures. The project employed the computational phase-field method in combination of microelasticity and electrostatic theories. The research was carried out in close collaboration with a number of experimental groups who used High Resolution Transmission Electron Microscopy (HRTEM), in situ TEM with Scanning Probe Microscopy (SPM), or Piezoresponse Force Microscopy (PFM) to characterize the domain structures and dynamics of ferroelectric thin films and heterostructures and who grow high-quality ferroelectric and multiferroic thin films using advanced growth techniques such as Molecular Beam Epitaxy (MBE), Pulsed Laser Deposition (PLLD), and sputtering. The research efforts of the project help establish the phase-field method as the most powerful method for understanding and predicting domain structures in ferroelectric thin films and nanostructures. The findings of the project led to the basic understanding of stability and switching mechanisms of ferroelectric domains under different mechanical boundary conditions and under either homogeneous capacitor configurations or local fields using metallic probes as electrodes and the emergence of charged domain walls during domain switching. The project predicted the spatial length scales, temperature ranges, and electromechanical conditions for different polar states in heterostructures and guided the discovery of both transient and stable novel polarization states containing vortex lattices in oxide superlattices. The project resulted in 134 journal publications and 6 PhD theses with all the PhD graduates currently working in either academia or industry within the United States. The basic understanding on the stability of mesoscale polar states and pattern evolution achieved by the project improved our ability to control and engineer properties of ferroelectric thin films and heterostructures for potential applications in nanoscale electronic devices.

36 MATERIALS SCIENCE↗

On the control, stability, and waiting time in a slotted ALOHA random-access system

This paper explores some of the boundaries in performance of slotted ALOHA systems by analyzing a simple and almost optimal centrally supervised control. The control results in a very simple Markov chain model and allows an examination of stability, conditional waiting time distribution of transmitting terminals, and many other system measures. The key to the simplicity is to have a probability of successful packet transmission that is independent of the number of transmitting terminals. In considering waiting time, we calculate the mean and other moments of the waiting time of a terminal when it enters the system to find (n - 1) other terminals already there competing for the channel. Under this control, the average time is proportional to n. The control requires exact knowledge of the number of terminals contending for the channel, and hence is not implementable, except as an approximation.

Ferguson, M. J.↗

Quantum spin helices more stable than the ground state: Onset of helical protection

Topological magnetic structures are promising candidates for resilient information storage. An elementary example is spin helices in one-dimensional easy-plane quantum magnets. To quantify their stability, we numerically implement the stochastic Schrödinger equation and time-dependent perturbation theory for spin chains with fluctuating local magnetic fields. Here, we find two classes of quantum spin helices that can reach and even exceed ground-state stability: spin-current-maximizing helices and, for fine-tuned boundary conditions, the recently discovered “phantom helices.” Beyond that, we show that the helicity itself (left or right rotating) is even more stable. We explain these findings by separated helical sectors and connect them to topological sectors in continuous spin systems. The resulting helical protection mechanism is a promising phenomenon for stabilizing helical quantum structures, e.g., in ultracold atoms and solid-state systems. We also identify a third type of phantom helix in the system.

1-dimensional spin chains↗

Stability investigations of airfoil flow by global analysis

As the result of global, non-parallel flow stability analysis the single value of the disturbance growth-rate and respective frequency is obtained. This complex value characterizes the stability of the whole flow configuration and is not referred to any particular flow pattern. The global analysis assures that all the flow elements (wake, boundary and shear layer) are taken into account. The physical phenomena connected with the wake instability are properly reproduced by the global analysis. This enhances the investigations of instability of any 2-D flows, including ones in which the boundary layer instability effects are known to be of dominating importance. Assuming fully 2-D disturbance form, the global linear stability problem is formulated. The system of partial differential equations is solved for the eigenvalues and eigenvectors. The equations, written in the pure stream function formulation, are discretized via FDM using a curvilinear coordinate system. The complex eigenvalues and corresponding eigenvectors are evaluated by an iterative method. The investigations performed for various Reynolds numbers emphasize that the wake instability develops into the Karman vortex street. This phenomenon is shown to be connected with the first mode obtained from the non-parallel flow stability analysis. The higher modes are reflecting different physical phenomena as for example Tollmien-Schlichting waves, originating in the boundary layer and having the tendency to emerge as instabilities for the growing Reynolds number. The investigations are carried out for a circular cylinder, oblong ellipsis and airfoil. It is shown that the onset of the wake instability, the waves in the boundary layer, the shear layer instability are different solutions of the same eigenvalue problem, formulated using the non-parallel theory. The analysis offers large potential possibilities as the generalization of methods used till now for the stability analysis.

Morzynski, Marek↗

Goertler instability in compressible boundary layers along curved surfaces with suction and cooling

The Goertler instability of the laminar compressible boundary layer flows along concave surfaces is investigated. The linearized disturbance equations for the three-dimensional, counter-rotating streamwise vortices in two-dimensional boundary layers are presented in an orthogonal curvilinear coordinate. The basic approximation of the disturbance equations, that includes the effect of the growth of the boundary layer, is considered and solved numerically. The effect of compressibility on critical stability limits, growth rates, and amplitude ratios of the vortices is evaluated for a range of Mach numbers for 0 to 5. The effect of wall cooling and suction of the boundary layer on the development of Goertler vortices is investigated for different Mach numbers.

El-Hady, N.↗