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

Conceptual design and structural analysis of the spectroscopy of the atmosphere using far infrared emission (SAFIRE) instrument

The conceptual design and structural analysis for the Spectroscopy of the Atmosphere using Far Infrared Emission (SAFIRE) Instrument are provided. SAFIRE, which is an international effort, is proposed for the Earth Observing Systems (EOS) program for atmospheric ozone studies. A concept was developed which meets mission requirements and is the product of numerous parametric studies and design/analysis iterations. Stiffness, thermal stability, and weight constraints led to a graphite/epoxy composite design for the optical bench and supporting struts. The structural configuration was determined by considering various mounting arrangements of the optical, cryo, and electronic components. Quasi-static, thermal, modal, and dynamic response analyses were performed, and the results are presented for the selected configuration.

Moses, Robert W.↗

Nonlinearly stable compact schemes for shock calculations

The applications of high-order, compact finite difference methods in shock calculations are discussed. The main concern is to define a local mean which will serve as a reference for introducing a local nonlinear limiting to control spurious numerical oscillations while maintaining the formal accuracy of the scheme. For scalar conservation laws, the resulting schemes can be proven total-variation stable in one space dimension and maximum-norm stable in multiple space dimensions. Numerical examples are shown to verify accuracy and stability of such schemes for problems containing shocks. These ideas can also be applied to other implicit schemes such as the continuous Galerkin finite element methods.

Cockburn, Bernardo↗

Time domain numerical calculations of unsteady vortical flows about a flat plate airfoil

A time domain numerical scheme is developed to solve for the unsteady flow about a flat plate airfoil due to imposed upstream, small amplitude, transverse velocity perturbations. The governing equation for the resulting unsteady potential is a homogeneous, constant coefficient, convective wave equation. Accurate solution of the problem requires the development of approximate boundary conditions which correctly model the physics of the unsteady flow in the far field. A uniformly valid far field boundary condition is developed, and numerical results are presented using this condition. The stability of the scheme is discussed, and the stability restriction for the scheme is established as a function of the Mach number. Finally, comparisons are made with the frequency domain calculation by Scott and Atassi, and the relative strengths and weaknesses of each approach are assessed.

Hariharan, S. I.↗

Designing reduced-order linear multivariable controllers using experimentally derived plant data

An iterative numerical algorithm for simultaneously improving multiple performance and stability robustness criteria for multivariable feedback systems is developed. The unsatisfied design criteria are improved by updating the free parameters of an initial, stabilizing controller's state-space matrices. Analytical expressions for the gradients of the design criteria are employed to determine a parameter correction that improves all of the feasible, unsatisfied design criteria at each iteration. A controller design is performed using the algorithm with experimentally derived data from a large space structure test facility. Experimental results of the controller's performance at the facility are presented.

Frazier, W. G.↗

Multigrid methods and high order finite difference for flow in transition - Effects of isolated and distributed roughness elements

The high order finite difference and multigrid methods have been successfully applied to direct numerical simulation (DNS) for flow transition in 3D channels and 3D boundary layers with 2D and 3D isolated and distributed roughness in a curvilinear coordinate system. A fourth-order finite difference technique on stretched and staggered grids, a fully-implicit time marching scheme, a semicoarsening multigrid method associated with line distributive relaxation scheme, and a new treatment of the outflow boundary condition, which needs only a very short buffer domain to damp all wave reflection, are developed. These approaches make the multigrid DNS code very accurate and efficient. This makes us not only able to do spatial DNS for the 3D channel and flat plate at low computational costs, but also able to do spatial DNS for transition in the 3D boundary layer with 3D single and multiple roughness elements. Numerical results show good agreement with the linear stability theory, the secondary instability theory, and a number of laboratory experiments.

Liu, C.↗

Multigrid methods for flow transition in three-dimensional boundary layers with surface roughness

The efficient multilevel adaptive method has been successfully applied to perform direct numerical simulations (DNS) of flow transition in 3-D channels and 3-D boundary layers with 2-D and 3-D isolated and distributed roughness in a curvilinear coordinate system. A fourth-order finite difference technique on stretched and staggered grids, a fully-implicit time marching scheme, a semi-coarsening multigrid method associated with line distributive relaxation scheme, and an improved outflow boundary-condition treatment, which needs only a very short buffer domain to damp all order-one wave reflections, are developed. These approaches make the multigrid DNS code very accurate and efficient. This allows us not only to be able to do spatial DNS for the 3-D channel and flat plate at low computational costs, but also to do spatial DNS for transition in the 3-D boundary layer with 3-D single and multiple roughness elements, which would have extremely high computational costs with conventional methods. Numerical results show good agreement with the linear stability theory, the secondary instability theory, and a number of laboratory experiments. The contribution of isolated and distributed roughness to transition is analyzed.

Liu, Chaoqun↗

Temporal model of an optically pumped co-doped solid state laser

Currently, research is being conducted on the optical properties of materials associated with the development of solid state lasers in the two micron region. In support of this effort, a mathematical model describing the energy transfer in a holmium laser sensitized with thulium is developed. In this paper, we establish some qualitative properties of the solution of the model, such as non-negativity, boundedness, and integrability. A local stability analysis is then performed from which conditions for asymptotic stability are attained. Finally, we report on our numerical analysis of the system and how it compares with experimental results.

Wangler, T. G.↗

Transient nucleation in glasses

Nucleation rates in condensed systems are frequently not at their steady state values. Such time dependent (or transient) nucleation is most clearly observed in devitrification studies of metallic and silicate glasses. The origin of transient nucleation and its role in the formation and stability of desired phases and microstructures are discussed. Numerical models of nucleation in isothermal and nonisothermal situations, based on the coupled differential equations describing cluster evolution within the classical theory, are presented. The importance of transient nucleation in glass formation and crystallization is discussed.

Kelton, K. F.↗

Astrophysical Adaptation of Points, the Precision Optical Interferometer in Space

POINTS (Precision Optical INTerferometer in Space) would perform microarcsecond optical astrometric measurements from space, yielding submicroarcsecond astrometric results from the mission. It comprises a pair of independent Michelson stellar interferometers and a laser metrology system that measures both the critical starlight paths and the angle between the baselines. The instrument has two baselines of 2 m, each with two subapertures of 35 cm; by articulating the angle between the baselines, it observes targets separated by 87 to 93 deg. POINTS does global astrometry, i.e., it measures widely separated targets, which yields closure calibration, numerous bright reference stars, and absolute parallax. Simplicity, stability, and the mitigation of systematic error are the central design themes. The instrument has only three moving-part mechanisms, and only one of these must move with sub-milliradian precision; the other two can tolerate a precision of several tenths of a degree. Optical surfaces preceding the beamsplitter or its fold flat are interferometrically critical; on each side of the interferometer, there are only three such. Thus, light loss and wavefront distortion are minimized. POINTS represents a minimalistic design developed ab initio for space. Since it is intended for astrometry, and therefore does not require the u-v-plane coverage of an imaging, instrument, each interferometer need have only two subapertures. The design relies on articulation of the angle between the interferometers and body pointing to select targets; the observations are restricted to the 'instrument plane.' That plane, which is fixed in the pointed instrument, is defined by the sensitive direction for the two interferometers. Thus, there is no need for siderostats and moving delay lines, which would have added many precision mechanisms with rolling and sliding parts that would be required to function throughout the mission. Further, there is no need for a third interferometer, as is required when out-of-plane observations are made. An instrument for astrometry, unlike those for imaging, can be compact and yet scientifically productive. The POINTS instrument is compact and therefore requires no deployment of precision structures, has no low-frequency (i.e., under 100 Hz) vibration modes, and is relatively easy to control thermally. Because of its small size and mass, it is easily and quickly repointed between observations. Further, because of the low mass, it can be economically launched into high Earth orbit which, in conjunction with a solar shield, yields nearly unrestricted sky coverage and a stable thermal environment.

Reasenberg, Robert D.↗

Perfectly Matched Layer for Linearized Euler Equations in Open and Ducted Domains

Recently, perfectly matched layer (PML) as an absorbing boundary condition has widespread applications. The idea was first introduced by Berenger for electromagnetic waves computations. In this paper, it is shown that the PML equations for the linearized Euler equations support unstable solutions when the mean flow has a component normal to the layer. To suppress such unstable solutions so as to render the PML concept useful for this class of problems, it is proposed that artificial selective damping terms be added to the discretized PML equations. It is demonstrated that with a proper choice of artificial mesh Reynolds number, the PML equations can be made stable. Numerical examples are provided to illustrate that the stabilized PML performs well as an absorbing boundary condition. In a ducted environment, the wave mode are dispersive. It will be shown that the group velocity and phase velocity of these modes can have opposite signs. This results in a confined environment, PML may not be suitable as an absorbing boundary condition.

Tam, Christopher K. W.↗

Redox Interactions between Iron and Carbon in Planetary Mantles: Implications for Degassing and Melting Processes

Carbon stability in planetary mantles has been studied by numerous authors because it is thought to be the source of C-bearing atmospheres and of C-rich lavas observed at the planetary surface. In the Earth, carbonaceous peridotites and eclogites compositions have been experimentally studied at mantle conditions [1] [2] [3]. [4] showed that the fO2 variations observed in martian meteorites can be explained by polybaric graphite-CO-CO2 equilibria in the Martian mantle. Based on thermodynamic calculations [4] and [5] inferred that the stable form of carbon in the source regions of the Martian basalts should be graphite (and/or diamond), and equilibrium with melts would be a source of CO2 for the martian atmosphere. Considering the high content of iron in the Martian mantle (approx.18.0 wt% FeO; [6]), compared to Earth s mantle (8.0 wt% FeO; [7]) Fe/C redox interactions should be studied in more detail.

Martin, A.↗

Development and Breakdown of Goertler Vortices in High Speed Boundary Layers

The nonlinear development of G rtler instability over a concave surface gives rise to a highly distorted stationary flow in the boundary layer that has strong velocity gradients in both spanwise and wall-normal directions. This distorted flow is susceptible to strong, high frequency secondary instability that leads to the onset of transition. For high Mach number flows, the boundary layer is also subject to the second mode instability. The nonlinear development of G rtler vortices and the ensuing growth and breakdown of secondary instability, the G rtler vortex interactions with second mode instabilities as well as oblique second mode interactions are examined in the context of both internal and external hypersonic configurations using nonlinear parabolized stability equations, 2-D eigenvalue analysis and direct numerical simulation. For G rtler vortex development inside the Purdue Mach 6 Ludwieg tube wind tunnel, multiple families of unstable secondary eigenmodes are identified and their linear and nonlinear evolution is examined. The computation of secondary instability is continued past the onset of transition to elucidate the physical mechanisms underlying the laminar breakdown process. Nonlinear breakdown scenarios associated with transition over a Mach 6 compression cone configuration are also explored.

Li, Fei↗

Instability Wave–Streak Interactions in a High Mach Number Boundary Layer at Flight Conditions

The interaction of stationary streaks undergoing nonmodal growth with modally unstable instability waves in a hypersonic boundary-layer flow is studied using numerical computations. The geometry and flow conditions are selected to match a relevant trajectory location from the ascent phase of the HIFiRE-1 flight experiment; namely, a 7 degree half-angle, circular cone with 2.5 mm nose radius, freestream Mach number equal to 5.30, unit Reynolds number equal to 13.42 m -1 , and wall-to-adiabatic temperature ratio of approximately 0.35 over most of the vehicle. This paper investigates the nonlinear evolution of initially linear optimal disturbances that evolve into finite-amplitude streaks, followed by an analysis of the modal instability characteristics of the perturbed, streaky boundary-layer flow. The investigation is performed with stationary direct numerical simulations (DNS) and plane-marching parabolized stability equations (PSE), in conjunction with partial-differential-equation-based planar eigenvalue analysis. The overall effect of streaks is to reduce the peak amplification factors of instability waves, indicating a possible downstream shift in the onset of laminar-turbulent transition. The present study confirms previous findings that the mean flow distortion of the nonlinear streak perturbation reduces the amplification rates of the Mack-mode instability. More importantly, however, the present results demonstrate that the spanwise varying component of the streak can produce a larger effect on the Mack-mode amplification. The study with selected azimuthal wavenumbers for the stationary streaks reveals that a wavenumber of approximately 1.4 times larger than the optimal wavenumber is more effective in stabilizing the planar Mack-mode instabilities. In the absence of unstable first-mode waves for the present cold-wall condition, transition onset is expected to be delayed until the peak streak amplitude increases to nearly 35 percent of the freestream velocity, when intrinsic instabilities of the boundary-layer streaks begin to dominate the transition process. For streak amplitudes below that limit a significant net stabilization is achieved, yielding a potential transition delay that can exceed 100 percent of the length of the laminar region in the uncontrolled case.

Boundary layer stability↗

On the stability of Kuzmin/Toomre disks

A simple two-component system consiting of a Kuzmin/Toomre disk embedded in a spherically symmetrical 'halo' is studied with Plummer's density distribution. The bounds on the mass fraction of the disk component which result from applying the semiempirical criterion against bar instabilities due to Ostriker and Peebles are studied analytically and the results are compared with recent numerical simulations. It is found that the locus t = 0.14, calculated using analytic approximations for a variety of rotation laws, agrees reasonably well with the neutral stability line derived from numerical results.

Frank, Juhan↗

Oscillatory Thermocapillary Convection

Stability analysis of thermocapillary convection in rectangular cavities is performed using direct numerical simulation. The influence of the Reynolds number (Re), the fluid Prandtl number (Pr) and cavity aspect ratio (Ar) on the motion is investigated. Neutral stability curves for transition to time-dependent convection are delineated in the Re-Ar plane for fluids with Pr=1.0, 4.4, 6.78, and 10. Several interesting features of these diagrams are discussed. One important conclusion is that Ar(epsilon tau) increases as Pr decreases. Thus large values of both Ar and Re are necessary to induce thermocapillary oscillations for small Pr fluids such as liquid metals and semiconductor melts. Energy analysis is also performed for the oscillatory flow in the neighborhood of critical points in order to gain insight into the mechanisms leading to instability.

Xu, Jieyong↗

Long-term risk analysis associated with nuclear waste disposal in space

An assessment and verification of previous analytic results on the long term risk of earth reentry for hazardous payloads is presented. The two areas were studied: (1) stability of nominal, near-circular storage orbits in the regions between Venus and earth and between earth and Mars, and (2) probability of earth reentry for off-nominal planet-crossing orbits resulting from deployment system failures. In the first area, numerical integrations of the equations of motion are compared with stability predications based on secular perturbation theory. The agreement is good in terms of the heliocentric distances covered and the general behavior of the orbital history, although certain near-resonance situations can lead to difficulty. In the second area, a Monte Carlo simulation of orbital evolution is used and the results compared with Opik's analytic theory of planetary encounters and collision statistics, with data verified to within a close order-of-magnitude.

Friedlander, A. L.↗

The properties of close multiple stars

A reexamination of the properties of close multiple systems is conducted, taking into account recent numerical, analytical, and observational results. The orbital elements of 43 spectroscopic multiple systems are presented in a table. Tables with photometry data of the systems and data describing their properties are also provided. Attention is given to period ratios, mass ratios, orbital coplanarity, stability, lithium abundances and age estimates, and prospects of future observations. Numerical studies indicate that fragmentation produces binary components with mass ratios usually close to 1.0, while fission usually produces mass ratios in the range 0.1-0.5. Using analytical stability criteria, all 27 close multiple systems considered are stable if they are corotating. If they are counter-rotating, four systems are possibly unstable. However, other factors suggest that even these four are probably stable. The solar-type close multiple systems ranged in age from 100 million years to 2,000 million years.

Fekel, F. C., Jr.↗

The influence of dynamic inflow and torsional flexibility on rotor damping in forward flight from symbolically generated equations

The combined effects of blade torsion and dynamic inflow on the aeroelastic stability of an elastic rotor blade in forward flight are studied. The governing sets of equations of motion (fully nonlinear, linearized, and multiblade equations) used in this study are derived symbolically using a program written in FORTRAN. Stability results are presented for different structural models with and without dynamic inflow. A combination of symbolic and numerical programs at the proper stage in the derivation process makes the obtainment of final stability results an efficient and straightforward procedure.

Reddy, T. S. R.↗