Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical 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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Properties of Nb x Ti (1–x) N thin films deposited on 300 mm silicon wafers for upscaling superconducting digital circuits

Scaling superconducting digital circuits requires fundamental changes in the current material set and fabrication process. The transition to 300 mm wafers and the implementation of advanced lithography are instrumental in facilitating mature CMOS processes, ensuring uniformity, and optimizing the yield. Here, this study explores the properties of Nb x Ti (1–x) N films fabricated by magnetron DC sputtering on 300 mm Si wafers. As a promising alternative to traditional Nb in device manufacturing, Nb x Ti (1–x) N offers numerous advantages, including enhanced stability and scalability to smaller dimensions, in both processing and design. As a ternary material, Nb x Ti (1–x) N allows engineering material parameters by changing deposition conditions. The engineered properties can be used to modulate device parameters through the stack and mitigate failure modes. We report characterization of Nb x Ti (1–x) N films at less than 2% thickness variability, 2.4% T c variability and 3% composition variability. Film resistivity (140–375 Ωcm) shows a strong correlation with the film oxygen content, while the critical temperature T c (4.6 K–14.1 K) is strongly affected by film stoichiometry and its microstructure has only a moderate effect on modifying T c . Our results offer insights about the interplay between film stoichiometry, film microstructure and critical temperature.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

The virtual element method for linear elastodynamics models: Design, analysis, and implementation

We design the conforming virtual element method for the numerical simulation of two dimensional time-dependent elastodynamics problems. We investigate the performance of the method both theoretically and numerically. We prove the stability and the convergence of the semi-discrete approximation in the energy norm and derive optimal error estimates. We also show the convergence in the L 2 norm. The performance of the virtual element method is assessed on a set of different computational meshes, including non-convex cells up to order four in the h-refinement setting. Exponential convergence is also experimentally seen in the p-refinement setting.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1962

GENERAL INSTABILITY OF ORTHOGONALLY STIFFENED CYLINDRICAL SHELLS

Earlier research at the National Aeronautics Research Institute (N.L.R.), Amsterdam, which forms the basis of recent work is reviewed. This early work refers to 2 schemes: the orthotropic shell and, in view of buckling modes where the half wave length is of the order of the ring distance, the shell with continuously distributed stringers and discrete rings. Linear theory is considered to be adequate for these structures, where the imperfections are small in comparison to the height of the ring sections. Recent developments account for pressure difference in addition to axial compression, for the correct stiffness matrix of skin panels in the post-buckling stage and for stringer bending due to hoop stresses in the skin, which are of importance as has been shown by the investigation of the post-buckling behaviour. Numerical data for the stiffness matrix of skin panels have been established. Numerical evaluation of the stability equation has not been performed as yet.

SHELL STABILITY↗
NASA NTRS1973

Numerical studies of unsteady two dimensional subsonic flows using the ICE method

A numerical program was developed to compute transient compressible and incompressible laminar flows in two dimensions with multicomponent mixing and chemical reaction. The algorithm used the Los Alamos Scientific Laboratory ICE (Implicit Continuous-Fluid Eulerian) method as its base. The program can compute both high and low speed compressible flows. The numerical program incorporating the stabilization techniques was quite successful in treating both old and new problems. Detailed calculations of coaxial flow very close to the entry plane were possible. The program treated complex flows such as the formation and downstream growth of a recirculation cell. An implicit solution of the species equation predicted mixing and reaction rates which compared favorably with the literature.

Wieber, P. R.↗
NASA NTRS1975

Upwind second-order difference schemes and applications in unsteady aerodynamic flows

Explicit second-order upwind difference schemes in combination with spatially symmetric schemes can produce larger stability bounds and better numerical resolution than symmetric schemes alone. However, if conservation form is essential, a special operator is required for transition between schemes. An operational approach has been devised for deriving transition operators so that strict conservation and local consistency are maintained. Various aspects of hybrid schemes are studied numerically for model linear and nonlinear equations. To demonstrate the utility of combining two different algorithms, MacCormack's explicit, noncentered, second-order method is combined with a completely upwind version, and numerical solutions of the Euler equations are obtained for two-dimensional, transonic flows with embedded supersonic regions and shock waves. The general utility of the operational approach for combining schemes is emphasized by deriving a second-order conservative scheme for the steady transonic small-disturbance potential equation.

Warming, R. F.↗
NASA NTRS1980

Heat- and chemical-resistant oxdiazole elastomers

Heat and chemical resistant polymers with triazine crosslinks are prepared by thermal condensation reactions to form 1,2,4-oxdiazole linkages. They are compounded with variety of fillers, extenders, and modifiers for numerous applications in which stability, impermiability to liquids and gases, good plasticity, and elasticity or rigidity are important.

Rosser, R. W.↗
NASA NTRS1982

Generation of the auroral kilometric radiation

Data collected from the S3-3 spacecraft in the auroral kilometric radiation (AKR) source region are employed to form a stability theory from which numerical results are obtained and discussed. The distribution function was found to be isotropic outside the auroral atmospheric loss-cone region, which was partially filled with upcoming electrons. A parallel electric field was observed to be modifying the loss-cone distribution. A model distribution function is formulated, along with a magnetic field model expressed in terms of the electron cyclotron frequency. Models for the parallel electric field are also introduced, the first with the field distributed over a broad altitude range, while the second assumes a potential drop only above 4000 km. The presence of the field is found to enhance the growth rate of the AKR. At higher altitudes, the cutoff frequency of the X mode is affected by the rising energetic electrons. Finally, the folding distances of spatial amplification are calculated.

Wu, C. S.↗
NASA NTRS1982

Mixed time integration methods for transient thermal analysis of structures

The computational methods used to predict and optimize the thermal structural behavior of aerospace vehicle structures are reviewed. In general, two classes of algorithms, implicit and explicit, are used in transient thermal analysis of structures. Each of these two methods has its own merits. Due to the different time scales of the mechanical and thermal responses, the selection of a time integration method can be a different yet critical factor in the efficient solution of such problems. Therefore mixed time integration methods for transient thermal analysis of structures are being developed. The computer implementation aspects and numerical evaluation of these mixed time implicit-explicit algorithms in thermal analysis of structures are presented. A computationally useful method of estimating the critical time step for linear quadrilateral element is also given. Numerical tests confirm the stability criterion and accuracy characteristics of the methods. The superiority of these mixed time methods to the fully implicit method or the fully explicit method is also demonstrated.

Liu, W. K.↗
NASA NTRS1983

Mixed time integration methods for transient thermal analysis of structures

The computational methods used to predict and optimize the thermal-structural behavior of aerospace vehicle structures are reviewed. In general, two classes of algorithms, implicit and explicit, are used in transient thermal analysis of structures. Each of these two methods has its own merits. Due to the different time scales of the mechanical and thermal responses, the selection of a time integration method can be a difficult yet critical factor in the efficient solution of such problems. Therefore mixed time integration methods for transient thermal analysis of structures are being developed. The computer implementation aspects and numerical evaluation of these mixed time implicit-explicit algorithms in thermal analysis of structures are presented. A computationally-useful method of estimating the critical time step for linear quadrilateral element is also given. Numerical tests confirm the stability criterion and accuracy characteristics of the methods. The superiority of these mixed time methods to the fully implicit method or the fully explicit method is also demonstrated.

Liu, W. K.↗
NASA NTRS1982

Flight testing the Rotor Systems Research Aircraft (RSRA)

The Rotor Systems Research Aircraft (RSRA) is a dedicated rotor test vehicle whose function is to fill the gap between theory, wind tunnel tests and flight verification data. Its flight test envelope has been designed to encompass the expected envelopes of future rotor systems under all flight conditions. The test configurations of the RSRA include pure helicopter and compound (winged helicopter) modes. In addition, should it become necessary to jettison an unstable rotor system in flight, the RSRA may be flown as a fixed wing aircraft. The heart of the RSRA's electronic flight control system is the TDY-43 computer, which can be programmed in numerous ways to change stability and control or force feel system gains. Computer programming changes allow the RSRA to be used as a five-degree-of-freedom inflight simulator for studying the handling qualities of research rotors.

Merrill, R. K.↗
NASA NTRS1983

Improvement of mixed time implicit-explicit algorithms for thermal analysis of structures

Computer implementation aspects and numerical evaluation of the recently introduced mixed time implicit-explicit algorithms in thermal analysis of structures are presented. A computationally useful method of estimating the critical time step for a linear quadrilateral element is given herein for the methods introduced by Liu and co-workers. Numerical tests confirm the stability criterion and accuracy characteristics of the methods. The superiority of these mixed time methods to the fully implicit method or the fully explicit method is also demonstrated.

Liu, W. K.↗
NASA NTRS1984

Effect of a forced Couette flow on coupled convective and morphological instabilities during unidirectional solidification

The effect of a forced Couette flow, parallel to a horizontal crystal-melt interface during directional solidification of an alloy of lead containing tin, on the onset of convective and morphological instabilities, is calculated numerically via a linear stability analysis. Such a flow does not affect perturbations with wave vectors perpendicular to the flow. For perturbations with wave vectors parallel to the flow, the onset of morphological instability is somewhat suppressed and thermosolutal convection is greatly suppressed. When instabilities occur, they are oscillatory and correspond to travelling waves. For values of the crystal growth velocity for which mixed morphological and convective modes occur, the presence of a forced flow produces sufficient decoupling to allow otherwise degenerate branches to be identified.

Coriell, S. R.↗
NASA NTRS1991

Waves produced from a harmonic point source in a supersonic boundary layer

The disturbance wave pattern produced by a harmonic point source in a compressible flat-plate boundary layer is computed using linear stability theory and direct numerical integration approach. Receptivity coefficients are computed for the spectrum of spanwise modes generated at the source. The effect of boundary layer growth on the development of linear waves is determined by using the method of multiple scales. Results are presented for Mach numbers of 0, 2, 4.5, and 7. It is found that disturbances spread in wedge-shaped regions behind the source and the wedge angle decreases with Mach number. The lateral spreading angle for the instability waves turns out to be quite close to the angle found experimentally for turbulence lateral contamination.

Balakumar, P.↗
NASA NTRS1991

Linear stability of free planetary waves in the presence of radiative-photochemical feedbacks

The diabatic effects of Newtonian cooling and ozone-dynamics interaction on the linear stability of free planetary waves in the atmosphere have been studied using a simple beta-plane model. The model couples radiative transfer, ozone advection, and ozone photochemistry with the quasi-geostrophic dynamical circulation. An analytical expression is derived which demonstrates the following: (1) the influence of meridional ozone advection on wave growth or decay depends on the wave and basic state vertical structures; and (2) photochemically accelerated cooling, which predominates in the upper stratosphere, augments the Newtonian cooling rate and is stabilizing. Attention is also given to the 1D linear stability problem which is numerically solved for a Charney basic state and for zonal mean basic states. It is shown that ozone heating generated by ozone-dynamics interaction in the stratosphere can reduce (enhance) the damping rates due to Newtonian cooling by as much as 50 percent for planetary waves of large vertical scale and maximum amplitude in the stratosphere.

Nathan, Terrence R.↗
NASA NTRS1993

An explicit solution to the optimal LQG problem for flexible structures with collocated rate sensors

We present a class of compensators in explicit form (not requiring numerical computer calculations) for stabilizing flexible structures with collocated rate sensors. They are based on the explicit solution, valid for both Continuum and FEM Models, of the LQG problem for minimizing mean square rate. They are robust with respect to system stability (will not destabilize modes even with mismatch of parameters), can be instrumented in state space form suitable for digital controllers, and can be specified directly from the structure modes and mode 'signature' (displacement vectors at sensor locations). Some simulation results are presented for the NASA LaRC Phase-Zero Evolutionary Model - a modal Trust model with 86 modes - showing damping ratios attainable as a function of compensator design parameters and complexity.

Balakrishnan, A. V.↗
NASA NTRS1992

The transition prediction toolkit: LST, SIT, PSE, DNS, and LES

The e(sup N) method for predicting transition onset is an amplitude ratio criterion that is on the verge of full maturation for three-dimensional, compressible, real gas flows. Many of the components for a more sophisticated, absolute amplitude criterion are now emerging: receptivity theory, secondary instability theory, parabolized stability equations approaches, direct numerical simulation and large-eddy simulation. This paper will provide a description of each of these new theoretical tools and provide indications of their current status.

Zang, Thomas A.↗
NASA NTRS1993

Slewing maneuvers and vibration control of space structures by feedforward/feedback moment-gyro controls

This paper presents a moment-gyro control approach to the maneuver and vibration suppression of a flexible truss arm undergoing a constant slewing motion. The overall slewing motion is triggered by a feedforward input, and a companion feedback controller is employed to augment the feedforward input and subsequently to control vibrations. The feedforward input for the given motion requirement is determined from the combined CMG (Control Momentum Gyro) devices and the desired rigid-body motion. The rigid-body dynamic model has enabled us to identify the attendant CMG momentum saturation constraints. The task for vibration control is carried out in two stages; first in the search of a suitable CMG placement along the beam span for various slewing maneuvers, and subsequently in the development of Liapunov-based control algorithms for CMG spin-stabilization. Both analytical and numerical results are presented to show the effectiveness of the present approach.

Yang, Li-Farn↗