Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Singularity”

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,639 records · Page 91
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

New results on the realizability of Reynolds stress turbulence closures

The realizability of Reynolds stress models in homogeneous turbulence is critically assessed from a theoretical standpoint. It is proven that a well known second-order closure formulated by Shih and Lumley using the strong realizability constraints of Schumann is, in fact, not a realizable model. The problem arises from the failure to properly satisfy the necessary positive second time derivative constraint when a principal Reynolds stress vanishes - a fatal flaw that becomes apparent when the non-analytic terms in their model are made single-valued as required on physical grounds. It is furthermore shown that the centrifugal acceleration generated by rotations of the principal axes of the Reynolds stress tensor can make the second derivative singular at the most extreme limits of realizable turbulence. This previously overlooked effect appears to make it impossible to identically satisfy the strong form of realizability in any version of the present generation of second-order closures. On the other hand, models properly formulated to satisfy the weak form of realizability - wherein states of one or two component turbulence are not accessible in finite time are found to be realizable. However, unlike the simpler and more commonly used second order closures, these models can be ill-behaved near the extreme limits of realizable turbulence due to the way that higher-degree nonlinearities are often unnecessarily introduced to satisfy realizability. Illustrative computations of homogeneous shear flows are presented to demonstrate these points which can have important implications for turbulence modeling.

Speziale, Charles G.↗
NASA NTRS1993

Nested Krylov methods and preserving the orthogonality

Recently the GMRESR inner-outer iteraction scheme for the solution of linear systems of equations was proposed by Van der Vorst and Vuik. Similar methods have been proposed by Axelsson and Vassilevski and Saad (FGMRES). The outer iteration is GCR, which minimizes the residual over a given set of direction vectors. The inner iteration is GMRES, which at each step computes a new direction vector by approximately solving the residual equation. However, the optimality of the approximation over the space of outer search directions is ignored in the inner GMRES iteration. This leads to suboptimal corrections to the solution in the outer iteration, as components of the outer iteration directions may reenter in the inner iteration process. Therefore we propose to preserve the orthogonality relations of GCR in the inner GMRES iteration. This gives optimal corrections; however, it involves working with a singular, non-symmetric operator. We will discuss some important properties, and we will show by experiments that, in terms of matrix vector products, this modification (almost) always leads to better convergence. However, because we do more orthogonalizations, it does not always give an improved performance in CPU-time. Furthermore, we will discuss efficient implementations as well as the truncation possibilities of the outer GCR process. The experimental results indicate that for such methods it is advantageous to preserve the orthogonality in the inner iteration. Of course we can also use iteration schemes other than GMRES as the inner method; methods with short recurrences like GICGSTAB are of interest.

Desturler, Eric↗
NASA NTRS1993

Energy transfer and constrained simulations in isotropic turbulence

The defining characteristic of turbulent flows is their ability to dissipate energy, even in the limit of zero viscosity. The Euler equations, if constrained in such a way that the velocity derivatives remain bounded, conserve energy. But when they arise as the limit of the Navier-Stokes (NS) equations, when the Reynolds number goes to infinity, there is persuasive empirical evidence that the gradients become singular as just the right function of Re for the dissipation to remain non-zero and to approach a well defined limit. It is generally believed that this limiting value of the dissipation is a property of the Euler equations themselves, independent of the particular dissipative mechanism involved, and that it can be normalized with the large scale properties of the turbulent flow (e.g. the kinetic energy per unit volume u'(exp 2)/2, and the integral scale L) without reference to the Reynolds number or to other dissipative quantities. This is usually taken to imply that the low wave number end of the energy spectrum, far from the dissipative range, is also independent of the particular mechanism chosen to dispose of the energy transfer. In the following sections, we present some numerical experiments on the effect of substituting different dissipation models into the truncated Euler equations. We will see that the effect is mainly felt in the 'near dissipation' range of the energy spectrum, but that this range can be quite wide in some cases, contaminating a substantial range of wave numbers. In the process, we will develop a 'practical' approximation to the subgrid energy transfer in isotropic turbulence, and we will gain insight into the structure of the nonlinear interactions among turbulent scales of comparable size, and into the nature of energy backscatter. Some considerations on future research directions are offered at the end.

Jimenez, Javier↗
NASA NTRS1994

The effects of anisotropy on the nonlinear behavior of bridged cracks in long strips

A model which can be used to predict the two-dimensional nonlinear behavior of bridged cracks in orthotropic strips is presented. The results obtained using a singular integral equation formulation which incorporates the anisotropy rigorously show that, although the effects of anisotropy are significant, the nondimensional quantities employed by Cox and Marshall can generate nearly universal results (R-curves, for example) for different levels of relative anisotropy. The role of composite constituent properties in the behavior of bridged cracks is clarified.

Ballarini, R.↗
NASA NTRS1993

A global approach to kinematic path planning to robots with holonomic and nonholonomic constraints

Robots in applications may be subject to holonomic or nonholonomic constraints. Examples of holonomic constraints include a manipulator constrained through the contact with the environment, e.g., inserting a part, turning a crank, etc., and multiple manipulators constrained through a common payload. Examples of nonholonomic constraints include no-slip constraints on mobile robot wheels, local normal rotation constraints for soft finger and rolling contacts in grasping, and conservation of angular momentum of in-orbit space robots. The above examples all involve equality constraints; in applications, there are usually additional inequality constraints such as robot joint limits, self collision and environment collision avoidance constraints, steering angle constraints in mobile robots, etc. The problem of finding a kinematically feasible path that satisfies a given set of holonomic and nonholonomic constraints, of both equality and inequality types is addressed. The path planning problem is first posed as a finite time nonlinear control problem. This problem is subsequently transformed to a static root finding problem in an augmented space which can then be iteratively solved. The algorithm has shown promising results in planning feasible paths for redundant arms satisfying Cartesian path following and goal endpoint specifications, and mobile vehicles with multiple trailers. In contrast to local approaches, this algorithm is less prone to problems such as singularities and local minima.

Divelbiss, Adam↗
NASA NTRS1993

Inertial-space disturbance rejection for space-based manipulators

The implementation of a disturbance rejection controller for a 6-DOF PUMA manipulator mounted on a 3-DOF platform was described. A control algorithm is designed to track the desired position and attitude of the end-effector in inertial space, subject to unknown disturbances in the platform axes. Experimental results are presented for step, sinusoidal, and random disturbances in the platform rotational axis and in the neighborhood of kinematic singularities. Robotic manipulators were proposed as a means of reducing the amount of extra vehicular activity time required for space station assembly and maintenance. The proposed scenario involves a robotic manipulator attached to some mobile platform, such as a spacecraft, satellite, or the space station itself. Disturbances in the platform position and attitude may prevent the manipulator from successfully completing the task. The possibility of using the manipulator to compensate for platform disturbances was explored. The problem of controlling a robotic manipulator on a mobile platform has received considerable attention in the past few years. Joshi and Desrochers designed a nonlinear feedback control law to carry out tasks (with respect to the robot base frame) in the presence of roll, pitch and yaw disturbances in the platform axes. Dubowsky, Vance, and Torres proposed a time-optimal planning algorithm for a robotic manipulator mounted on a spacecraft, subject to saturation limits in the attitude control reaction jets. Papadopoulos and Dubowsky developed a general framework for analyzing the control of free-floating space manipulator systems. Most recently, Torres and Dubowsky have presented a technique called the enhanced disturbance map to find manipulator trajectories that reduce the effect of disturbances in the spacecraft position and attitude. One common assumption in the literature is that the disturbance signal is exactly known. If this is the case, then the end-effector location can be calculated without relying on direct end-point sensing. However, this assumption is invalid if there is a significant delay in the platform position and attitude measurements, or if the kinematics of the platform are not well known, or if the platform is a non-rigid structure. In the more likely case that only the nominal platform location and upper bound on the disturbance signal are known, direct end-point sensing is needed to measure the end-effector location.

Holt, Kevin↗
NASA NTRS1993

A method for the modelling of porous and solid wind tunnel walls in computational fluid dynamics codes

Porous wall wind tunnels have been used for several decades and have proven effective in reducing wall interference effects in both low speed and transonic testing. They allow for testing through Mach 1, reduce blockage effects and reduce shock wave reflections in the test section. Their usefulness in developing computational fluid dynamics (CFD) codes has been limited, however, by the difficulties associated with modelling the effect of a porous wall in CFD codes. Previous approaches to modelling porous wall effects have depended either upon a simplified linear boundary condition, which has proven inadequate, or upon detailed measurements of the normal velocity near the wall, which require extensive wind tunnel time. The current work was initiated in an effort to find a simple, accurate method of modelling a porous wall boundary condition in CFD codes. The development of such a method would allow data from porous wall wind tunnels to be used more readily in validating CFD codes. This would be beneficial when transonic validations are desired, or when large models are used to achieve high Reynolds numbers in testing. A computational and experimental study was undertaken to investigate a new method of modelling solid and porous wall boundary conditions in CFD codes. The method utilized experimental measurements at the walls to develop a flow field solution based on the method of singularities. This flow field solution was then imposed as a pressure boundary condition in a CFD simulation of the internal flow field. The effectiveness of this method in describing the effect of porosity changes on the wall was investigated. Also, the effectiveness of this method when only sparse experimental measurements were available has been investigated. The current work demonstrated this approach for low speed flows and compared the results with experimental data obtained from a heavily instrumented variable porosity test section. The approach developed was simple, computationally inexpensive, and did not require extensive or intrusive measurements of the boundary conditions during the wind tunnel test. It may be applied to both solid and porous wall wind tunnel tests.

Beutner, Thomas John↗
NASA NTRS1993

Waves in space plasma dipole antenna subsystem

The Waves In Space Plasma (WISP) flight experiment requires a 50-meter-long deployable dipole antenna subsystem (DASS) to radiate radio frequencies from the STS Orbiter cargo bay. The transmissions are to excite outer ionospheric plasma between the dipole and a free-flying receiver (Spartan) for scientific purposes. This report describes the singular DASS design requirements and how the resulting design satisfies them. A jettison latch is described in some detail. The latch releases the antenna in case of any problems which might prevent the bay doors from closing for re-entry and landing of the Orbiter.

Thomson, Mark↗
NASA NTRS1992

Analysis of crustal structure of Venus utilizing residual Line-of-Sight (LOS) gravity acceleration and surface topography data. A trial of global modeling of Venus gravity field using harmonic spline method

To construct Venus' gravity disturbance field (or gravity anomaly) with the spacecraft-observer line of site (LOS) acceleration perturbation data, both a global and a local approach can be used. The global approach, e.g., spherical harmonic coefficients, and the local approach, e.g., the integral operator method, based on geodetic techniques are generally not the same, so that they must be used separately for mapping long wavelength features and short wavelength features. Harmonic spline, as an interpolation and extrapolation technique, is intrinsically flexible to both global and local mapping of a potential field. Theoretically, it preserves the information of the potential field up to the bound by sampling theorem regardless of whether it is global or local mapping, and is never bothered with truncation errors. The improvement of harmonic spline methodology for global mapping is reported. New basis functions, a singular value decomposition (SVD) based modification to Parker & Shure's numerical procedure, and preliminary results are presented.

Fang, Ming↗
NASA NTRS1992

Digital image profilers for detecting faint sources which have bright companions

For this program, an image profiling system was developed which offers the potential for detecting extremely faint optical sources that are located in close proximity to bright companions. The approach employed is novel in three respects. First, it does not require an optical system wherein extraordinary measures must be taken to minimize diffraction and scatter. Second, it does not require detectors possessing either extreme uniformity in sensitivity or extreme temporal stability. Finally, the system can readily be calibrated, or nulled, in space by testing against an unresolved singular stellar source.

Morris, Elena↗
NASA NTRS1992

Tension fracture of laminates for transport fuselage. Part 1: Material screening

Transport fuselage structures are designed to contain pressure following a large penetrating damage event. Application of composites to fuselage structures requires a data base and supporting analysis on tension damage tolerance. Tests with 430 fracture specimens were used to accomplish the following: (1) identify critical material and laminate variables affecting notch sensitivity, (2) evaluate composite failure criteria, and (3) recommend a screening test method. Variables studied included fiber type, matrix toughness, lamination manufacturing process, and intraply hybridization. The laminates found to have the lowest notch sensitivity were manufactured using automated tow placement. This suggests a possible relationship between the stress distribution and repeatable levels of material inhomogeneity that are larger than found in traditional tape laminates. Laminates with the highest notch sensitivity consisted of toughened matrix materials that were resistant to a splitting phenomena that reduces stress concentrations in load bearing plies. Parameters for conventional fracture criteria were found to increase with the crack length of the smallest notch sizes studied. Most materials and laminate combinations followed less than a square root singularity for the largest crack sizes studied. Specimen geometry, notch type, and notch size were evaluated in developing a screening test procedure. Results indicate that a range of notch sizes must be tested to determine notch sensitivity.

Walker, T. H.↗
NASA NTRS1994

The effect of inlet swirl on the dynamics of long annular seals in centrifugal pumps

This paper describes additional results from a continuing research program which aims to identify the dynamics of long annular seals in centrifugal pumps. A seal test rig designed at Heriot-Watt University and commissioned at Weir Pumps Research Laboratory in Alloa permits the identification of mass, stiffness, and damping coefficients using a least-squares technique based on the singular value decomposition method. The analysis is carried out in the time domain using a multi-fiequency forcing function. The experimental method relies on the forced excitation of a flexibly supported stator by two hydraulic shakers. Running through the stator embodying two symmetrical balance drum seals is a rigid rotor supported in rolling element bearings. The only physical connection between shaft and stator is the pair of annular gaps filled with pressurized water discharged axially. The experimental coefficients obtained from the tests are compared with theoretical values.

Ismail, M.↗
NASA NTRS1992

Characterization of dynamic stall phenomenon using two-dimensional unsteady Navier-Stokes equations

Among the new significant aspects of the present work are: (1) the treatment of the far-field boundary; (2) the use of C-grid topology, with the branch-cut singularity treated analytically; (3) evaluation of the effect of the envelope of prevailing initial states, and finally; (4) the ability to employ streakline/pathline 'visualization' to probe the unsteady features prevailing in vortex-dominated flows. The far-field boundary is placed at infinity, using appropriate grid stretching. This contributes to the accuracy of the solutions, but raised a number of important issues which needed to be resolved; this includes determining the equivalent time-dependent circulation for the pitching airfoil. A secondary counter-clockwise vortex erupts from within the boundary layer and immediately pinches off the energetic leading-edge shear layer which then, through hydrodynamic instability, rolls up into the dynamic stall vortex. The streakline/pathline visualization serves to provide information for insight into the physics of the unsteady separated flow.

Ghia, K. N.↗
NASA NTRS1994

Tensor-GMRES method for large sparse systems of nonlinear equations

This paper introduces a tensor-Krylov method, the tensor-GMRES method, for large sparse systems of nonlinear equations. This method is a coupling of tensor model formation and solution techniques for nonlinear equations with Krylov subspace projection techniques for unsymmetric systems of linear equations. Traditional tensor methods for nonlinear equations are based on a quadratic model of the nonlinear function, a standard linear model augmented by a simple second order term. These methods are shown to be significantly more efficient than standard methods both on nonsingular problems and on problems where the Jacobian matrix at the solution is singular. A major disadvantage of the traditional tensor methods is that the solution of the tensor model requires the factorization of the Jacobian matrix, which may not be suitable for problems where the Jacobian matrix is large and has a 'bad' sparsity structure for an efficient factorization. We overcome this difficulty by forming and solving the tensor model using an extension of a Newton-GMRES scheme. Like traditional tensor methods, we show that the new tensor method has significant computational advantages over the analogous Newton counterpart. Consistent with Krylov subspace based methods, the new tensor method does not depend on the factorization of the Jacobian matrix. As a matter of fact, the Jacobian matrix is never needed explicitly.

Feng, Dan↗
NASA NTRS1993

On the cause of thin current sheets in the near-Earth magnetotail and their possible significance for magnetospheric substorms

The formation of thin current sheets in the near-Earth magnetotail during substorm growth phases is addressed in terms of a simple model. An appropriate part of the unperturbed magnetotail is represented by a plane sheet model. Perturbations are applied to the upper and to the left boundaries, representing the magnetopause and the near-Earth tail boundary, where the perturbation at the latter models the interaction between the tail and the inner magnetosphere. Treating the perturbation as ideal (dissipation-free), we found that singular current sheets develop in the midplane of the tail. The analytical results are explored numerically. Using realistic dimensions of the domain considered, the influence of the earthward boundary on current sheet formation dominates. It is argued that current sheet formation of this type plays an important role in the processes associated with the onset of magnetospheric substorms.

Schindler, K.↗
NASA NTRS1993

The piecewise-linear predictor-corrector code - A Lagrangian-remap method for astrophysical flows

We describe a time-explicit finite-difference algorithm for solving the nonlinear fluid equations. The method is similar to existing Eulerian schemes in its use of operator-splitting and artificial viscosity, except that we solve the Lagrangian equations of motion with a predictor-corrector and then remap onto a fixed Eulerian grid. The remap is formulated to eliminate errors associated with coordinate singularities, with a general prescription for remaps of arbitrary order. We perform a comprehensive series of tests on standard problems. Self-convergence tests show that the code has a second-order rate of convergence in smooth, two-dimensional flow, with pressure forces, gravity, and curvilinear geometry included. While not as accurate on idealized problems as high-order Riemann-solving schemes, the predictor-corrector Lagrangian-remap code has great flexibility for application to a variety of astrophysical problems.

Lufkin, Eric A.↗
NASA NTRS1993

The piecewise-parabolic method in curvilinear coordinates

We derive interpolation formulae for a third-order finite difference method in curvilinear, orthogonal coordinate systems. These formulae serve as a supplement to Colella and Woodward's PPM scheme for problems where the coordinate origin is included in the computational domain. Numerical examples of the improved accuracy of the advection scheme near coordinate singularities are shown.

Blondin, John M.↗
NASA NTRS1992

Balanced systems and structures - Reduction, assignment, and perturbations

New results on properties of balanced linear systems and structures are presented. Balanced representation is defined for systems with poles at an imaginary axis or at the origin. Grammians do not exist in this case, but antigrammians are introduced which make balanced reduction possible. System grammians of specified properties are obtained by assigning the sensor and actuator configuration. The latter is determined using input-output assignment procedures introduced for general systems and specified to structures. A system is said to be uniformly balanced if all its Hankel singular values are equal.

Gawronski, Wodek↗