Search NASASearch

SEARCH · Search NASA

Results for “monotone stability”

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 37 records · Page 2

Separated flows near the nose of a body of revolution

The solution of the Navier-Stokes equations for the problem of cross-flow separataion about a deforming cylinder was achieved by iteration. It was shown that the separation starts at the rear stagnation point and the point of primary separation moves upstram along the cylinder surface. A general method of linear stability analysis for nonparallel external flows was constructed, which consists of representing the eigenfunctions with complete orthogonal sets and forms characteristic equations with the Galerkin method. The method was applied to the Kovasznay flow which is an exact solution of the Navier-Stokes equation. The results show that when the critical parameter is exceeded, there are only a few isolated unstable eigen-frequencies. Another exact solution is shown to be absolutely and monotonically stable with respect to infinitesimal disturbances of all frequencies. The flow is also globally, asymptotically, and monotonically stable in the mean with respect o three-dimensional disturbances. This result forms the sound foundation of rigorous stability analysis for nonparallel flows, and provides an invaluable test ground for future studies of nonparallel flows in which the basic states do not posses exact solutions. The application of this method to the study of the formation of spiral vorticies near the nose of a rotating body of revolution is underway. The same method will be applied to the stability analysis of reversed flow over a plate with suction.

Lin, S. P.

A Dynamical Model for The Jupiter Thermocline

A new conceptual model for the observed banded wind structure of the giant outer planets is proposed, assuming a smoothly distributed potential vorticity (PV) atop a layer of deeply seated, latitudinally variable stratification. With a vanishing horizontal entropy contrast presumably enforced at the bottom of the flow layer by the underlying convection, the thermal-wind balance of observed cloud-top motions implies a mapping of constant potential temperature surfaces, nearly vertical at upper tropospheric levels, down to a deeper, flatter, but variably sloped 'thermocline.' For a generally colder-poleward distribution of isentropes consistent with the strong equatorial jets of both Jupiter and Saturn, this kind of mapping would also imply a poleward decreasing static stability. The proposed temperature-stability distribution is just the reverse of the Earth's troposphere, where the strongest latitudinal potential temperature gradients are at the bottom instead of the top, with static stability generally increasing toward the pole. The warmer-stabler correlation for Jupiter would be consistent, however, with a nearly monotonic distribution of PV from low to high latitudes, over planetary scales for which the planetary vorticity dominates the relative vorticity of the jets. In this way efficient PV mixing for the isentropically bounded thermocline can account for the dynamical maintenance of the cyclonic flanks of the equatorial jet. Over latitudinal intervals comparable to the internal deformation scale, however, the gradient of the absolute vorticity is dominated by the flow curvature, with the local latitudinal gradient of the static stability proportional to the local potential temperature gradient and therefore the local variation in geostrophic velocity. The resulting correlation of velocities and vorticity gradients implies an alternation of the flow over the internal deformation radius set by the vertical entropy contrast.

Allison, Michael

Correlation between Mechanical Behavior and Actuator-type Performance of Ni-Ti-Pd High-temperature Shape Memory Alloys

High-temperature shape memory alloys in the NiTiPd system are being investigated as lower cost alternatives to NiTiPt alloys for use in compact solid-state actuators for the aerospace, automotive, and power generation industries. A range of ternary NiTiPd alloys containing 15 to 46 at.% Pd has been processed and actuator mimicking tests (thermal cycling under load) were used to measure transformation temperatures, work behavior, and dimensional stability. With increasing Pd content, the work output of the material decreased, while the amount of permanent strain resulting from each load-biased thermal cycle increased. Monotonic isothermal tension testing of the high-temperature austenite and low temperature martensite phases was used to partially explain these behaviors, where a mismatch in yield strength between the austenite and martensite phases was observed at high Pd levels. Moreover, to further understand the source of the permanent strain at lower Pd levels, strain recovery tests were conducted to determine the onset of plastic deformation in the martensite phase. Consequently, the work behavior and dimensional stability during thermal cycling under load of the various NiTiPd alloys is discussed in relation to the deformation behavior of the materials as revealed by the strain recovery and monotonic tension tests.

Bigelow, Glen S.

The influence of thermal inertia on temperatures and frost stability on Triton

It is presently argued, in view of (1) a thermal inertia model for the surface of Triton which (like previous ones) predicts a monotonic recession of permanent N2 deposits toward the poles and very little seasonal N2 frost in the southern hemisphere, and (2) new spectroscopic evidence for nonvolatile CO2 on Triton's bright southern hemisphere, that much of that bright southern material is not N2. Such bright southern hemisphere volatiles may allow the formation of seasonal frosts, thereby helping to explain the observed spectroscopic changes of Triton during the last decade.

Spencer, John R.

An 8-b 1.3-MHz successive-approximation A/D converter

A new successive-approximation A/D converter is described. It combines a string of equal-valued polysilicon resistors and a set of ratioed capacitors in a unique circuit configuration so that a high sampling rate is achieved. The comparator is realized by a chopper-stabilized amplifier to reduce the effect of the offset voltages of MOS amplifiers. The converter performs an 8-b monotonic conversion with a differential nonlinearity less than 1 LSB in 770 nsec. The die area is 3750 sq mil. This new conversion technique can also be utilized in a pipelined A/D converter (Temes, 1985) and enables it to achieve high speed.

Hadidi, KH.

Positivity-preserving numerical schemes for multidimensional advection

This report describes the construction of an explicit, single time-step, conservative, finite-volume method for multidimensional advective flow, based on a uniformly third-order polynomial interpolation algorithm (UTOPIA). Particular attention is paid to the problem of flow-to-grid angle-dependent, anisotropic distortion typical of one-dimensional schemes used component-wise. The third-order multidimensional scheme automatically includes certain cross-difference terms that guarantee good isotropy (and stability). However, above first-order, polynomial-based advection schemes do not preserve positivity (the multidimensional analogue of monotonicity). For this reason, a multidimensional generalization of the first author's universal flux-limiter is sought. This is a very challenging problem. A simple flux-limiter can be found; but this introduces strong anisotropic distortion. A more sophisticated technique, limiting part of the flux and then restoring the isotropy-maintaining cross-terms afterwards, gives more satisfactory results. Test cases are confined to two dimensions; three-dimensional extensions are briefly discussed.

Leonard, B. P.

Ab Initio Molecular Dynamics Insights into Stress Corrosion Cracking and Dissolution of Metal Oxides

Oxide phases such as α-Fe2O3 (hematite) and α-Al2O3 (corundum) are highly insoluble in water; however, subcritical crack growth has been observed in humidity nonetheless. Chemically induced bond breaking at the crack tip appears unlikely due to sterically hindered molecular transport. The molecular mechanics of a crack in corundum with a reactive force field reveal minimal lattice trapping, leading to bond breaking before sufficient space opens for water transport. To address this, we model a pre-built blunt crack with space for H2O molecule adsorption at the tip and show that it reduces fracture toughness by lowering the critical J-integral. Then, we explore stress-enhanced dissolution to understand the mechanism of crack tip blunting in the oxide/water system. Density functional theory combined with metadynamics was employed to describe atomic dissolution from flat hematite and corundum surfaces in pure water. Strain accelerates dissolution, stabilizing intermediate states with broken bonds before full atom detachment, while the free energy profile of unstrained surfaces is almost monotonic. The atomistic calculations provided input for a kinetic model, predicting the shape evolution of a blunt crack tip, which displays three distinct regimes: (i) dissolution primarily away from the tip, (ii) enhanced blunting near but not at the apex, and (iii) sharpening near the apex. The transition between regimes occurs at a low strain, highlighting the critical role of water in the subcritical crack growth of oxide scales, with dissolution as the fundamental microscopic mechanism behind this process.

Chemistry

Fracture Analysis of Semi-Elliptical Surface Cracks in Ductile Materials

Accurate life assessment of structural components may require advanced life prediction criteria and methodologies. Structural components often exhibit several different types of defects, among the most prevalent being surface cracks. A semi-elliptical surface crack subjected to monotonic loading will exhibit stable crack growth until the crack has reached a critical size, at which the crack loses stability and fracture ensues (Newman, 2000). The shape and geometry of the flaw are among the most influential factors. When considering simpler crack configurations, such as a through-the-thickness crack, a three-dimensional (3D) geometry may be modeled under the approximation of two-dimensional (2D) plane stress or plane strain. The more complex surface crack is typically modeled numerically with the Finite Element Method (FEM). A semi-elliptical surface crack is illustrated in Figure 1-1.

Daniewicz, S. R.

Input-output-controlled nonlinear equation solvers

To upgrade the efficiency and stability of the successive substitution (SS) and Newton-Raphson (NR) schemes, the concept of input-output-controlled solvers (IOCS) is introduced. By employing the formal properties of the constrained version of the SS and NR schemes, the IOCS algorithm can handle indefiniteness of the system Jacobian, can maintain iterate monotonicity, and provide for separate control of load incrementation and iterate excursions, as well as having other features. To illustrate the algorithmic properties, the results for several benchmark examples are presented. These define the associated numerical efficiency and stability of the IOCS.

Padovan, Joseph

Spectral State Evolution of 4U 1820-30: the Stability of the Spectral Index of Comptonization Tail

We analyze the X-ray spectra and their timing properties of the compact Xray binary 4U 1820-30. We establish spectral transitions in this source seen with BeppoSAX and the Rossi X-ray Timing Explorer (RXTE). During the RXTE observations (1996 - 2009), the source were approximately approximately 75% of its time in the soft state making the lower banana and upper banana transitions combined with long-term low-high state transitions. We reveal that all of the X-ray spectra of 4U 1820-30 are fit by a composition of a thermal (blackbody) component, a Comptonization component (COMPTB) and a Gaussian-line component. Thus using this spectral analysis we find that the photon power-law index Gamma of the Comptonization component is almost unchangeable (Gamma approximately 2) while the electron temperature kTe changes from 2.9 to 21 keV during these spectral events. We also establish that for these spectral events the normalization of COMPTB component (which is proportional to mass accretion rate ˙M) increases by factor 8 when kTe decreases from 21 keV to 2.9 keV. Before this index stability effect was also found analyzing X-ray data for Z-source GX 340+0 and for atolls, 4U 1728-34, GX 3+1. Thus, we can suggest that this spectral stability property is a spectral signature of an accreting neutron star source. On the other hand in a black hole binary Γ monotonically increases with ˙Mand ultimately its value saturates at large ˙M.

evolution

Spectral stability of Taylor's vortex array

In a global sense it is shown that the two-dimensional Taylor vortex array, an exact solution of the Navier-Stokes equation, is absolutely and monotonically stable with respect to infinitesimal disturbances of all discrete frequencies as long as the viscosity is positive. It is suggested that the Taylor vortex array may also be stable with respect to finite amplitude disturbances.

Lin, S. P.

High-Order Semi-Discrete Central-Upwind Schemes for Multi-Dimensional Hamilton-Jacobi Equations

We present the first fifth order, semi-discrete central upwind method for approximating solutions of multi-dimensional Hamilton-Jacobi equations. Unlike most of the commonly used high order upwind schemes, our scheme is formulated as a Godunov-type scheme. The scheme is based on the fluxes of Kurganov-Tadmor and Kurganov-Tadmor-Petrova, and is derived for an arbitrary number of space dimensions. A theorem establishing the monotonicity of these fluxes is provided. The spacial discretization is based on a weighted essentially non-oscillatory reconstruction of the derivative. The accuracy and stability properties of our scheme are demonstrated in a variety of examples. A comparison between our method and other fifth-order schemes for Hamilton-Jacobi equations shows that our method exhibits smaller errors without any increase in the complexity of the computations.

Bryson, Steve

Theoretical analyses of Baroclinic flows

A stability analysis of a thin horizontal rotating fluid layer which is subjected to arbitrary horizontal and vertical temperature gradients is presented. The basic state is a nonlinear Hadley cell which contains both Ekman and thermal boundary layers; it is given in closed form. The stability analysis is based on the linearized Navier-Stokes equations, and zonally symmetric perturbations in the form of waves propagating in the meridional direction are considered. Numerical methods were used for the stability problem. It was found that the instability sets in when the Richardson number is close to unity and that the critical Richardson number is a non-monotonic function of the Prandtl number. Further, it was found that the critical Richardson number decreases with increasing Ekman number until a critical value of the Ekman number is reached beyond which the fluid is stable.

Antar, B.

Improvements in the accuracy and stability of algorithms for the small-disturbance and full-potential equations applied to transonic flows

Numerical techniques that improve the accuracy and stability of algorithms for the small disturbance and full potential equations used to calculate transonic flows are described. For the small disturbance equation, the algorithm improvements are: (1) the use of monotone switches in the type dependent finite differencing, and (2) the use of stable and simple second order accurate spatial differencing; these improvements are for steady and unsteady transonic flows. For the steady full potential equation, the improvement is in the use of a monotone switch in the type dependent finite differencing of an approximate factorization (AF2) algorithm. All these improvements are implemented in present computer codes by making minor coding modifications.

Goorjian, P. M.

Microscopic Origin of Twist-Dependent Electron Transfer Rate in Bilayer Graphene

Using molecular simulation and continuum dielectric theory, we consider how electrochemical kinetics are modulated by the twist angle in bilayer graphene electrodes. By establishing a connection between the twist angle and the screening length of charge carriers within the electrode, we investigate how tunable metallicity modifies the statistics of the electron transfer energy gap. Constant potential molecular simulations show that the activation free energy for electron transfer increases with screening length, leading to a non-monotonic dependence on the twist angle. Here, the twist angle alters the density of states, tuning the number of thermally accessible channels for electron transfer and the reorganization energy by affecting the stability of the vertically excited state through attenuated image charge interactions. Understanding these effects allows us to express the Marcus rate of interfacial electron transfer as a function of the twist angle in a manner consistent with experimental observations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

Tensile and Compressive Constitutive Response of 316 Stainless Steel at Elevated Temperatures

Creep rate in compression is lower by factors of 2 to 10 than in tension if the microstructure of the two specimens is the same and are tested at equal temperatures and equal but opposite stresses. Such behavior is characteristic for monotonic creep and conditions involving cyclic creep. In the latter case creep rate in both tension and compression progressively increases from cycle to cycle, rendering questionable the possibility of expressing a time stabilized constitutive relationship. The difference in creep rates in tension and compression is considerably reduced if the tension specimen is first subjected to cycles of tensile creep (reversed by compressive plasticity), while the compression specimen is first subjected to cycles of compressive creep (reversed by tensile plasticity). In both cases, the test temperature is the same and the stresses are equal and opposite. Such reduction is a reflection of differences in microstructure of the specimens resulting from different prior mechanical history.

Manson, S. S.

Characterization of Magnetoacoustic Emission Related to Structural Properties of Ferromagnets

An extensive study of magnetoacoustic emission (MAE) properties has been performed over the past several years. As a result, the dependence of the spectral characteristics of MAE on certain microstructural variations and uniaxially applied stress in a particular type of low carbon steel are now well known. The embrittlement-causing concentration of certain atomic species, e. g., tin, sulphur, phosphorous etc., at the grain boundaries of this steel creates strong potential barriers resisting the motion of non-180 domain walls which is the source of MAE. (Since the only type of non-180 domain walls in this material are 90 domain walls, the term 90 domain wall will be used throughout this paper in place of non-180 domain wall.) An MAE burst produced during one-half cycle of a hysteresis loop at a low AC magnetic field frequency (e. g., 0.7 Hz) shows two sub-peaks; the leading peak is usually sharp and short-lived, while the trailing peak is usually smooth and long-lasting. It has been shown that the enhanced domain wall-defect interaction, due to the strengthened potential barriers, causes an increase in the asymmetry of the MAE signal by suppressing the leading sub-peak and amplifying the trailing sub-peak. This phenomena is due to the delayed motion of the 90 domain walls. The effect of a tensile stress applied parallel to the external AC magnetic field is to diminish the MAE. On the other hand, the amplitude of the MAE burst has been shown to be a non-monotonic function of the stress amplitude. Recently, our study has concentrated on obtaining quantitative values for parameters computed from the MAE spectra averaged over a sufficient number of cycles to achieve statistical stability. Nevertheless, certain fundamental elements of the MAE characteristics remain unexplained.

Namkung, M.

Tensile and compressive constitutive response of 316 stainless steel at elevated temperatures

It is demonstrated that creep rate of 316 SS is lower by factors of 2 to 10 in compression than in tension if the microstructure is the same and tests are conducted at identical temperatures and equal but opposite stresses. Such behavior was observed for both monotonic creep and conditions involving cyclic creep. In the latter case creep rate in both tension and compression progressively increases from cycle to cycle, rendering questionable the possibility of expressing a time-stabilized constitutive relationship. The difference in creep rates in tension and compression is considerably reduced if the tension specimen is first subjected to cycles of tensile creep (reversed by compressive plasticity), while the compression specimen is first subjected to cycles of compressive creep (reversed by tensile plasticity). In both cases, the test temperature is the same and the stresses are equal and opposite. Such reduction is a reflection of differences in microstructure of the specimens resulting from different prior mechanical history.

Manson, S. S.