Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “stiff problem”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
NASA NTRS1983

Application of the method of lines for solutions of the Navier-Stokes equations using a nonuniform grid distribution

The feasibility of the method of lines for solutions of physical problems requiring nonuniform grid distributions is investigated. To attain this, it is also necessary to investigate the stiffness characteristics of the pertinent equations. For specific applications, the governing equations considered are those for viscous, incompressible, two dimensional and axisymmetric flows. These equations are transformed from the physical domain having a variable mesh to a computational domain with a uniform mesh. The two governing partial differential equations are the vorticity and stream function equations. The method of lines is used to solve the vorticity equation and the successive over relaxation technique is used to solve the stream function equation. The method is applied to three laminar flow problems: the flow in ducts, curved-wall diffusers, and a driven cavity. Results obtained for different flow conditions are in good agreement with available analytical and numerical solutions. The viability and validity of the method of lines are demonstrated by its application to Navier-Stokes equations in the physical domain having a variable mesh.

Abolhassani, J. S.↗
NASA NTRS1991

Resolved rate and torque control schemes for large scale space based kinematically redundant manipulators

Resolved rate control of kinematically redundant ground based manipulators is a challenging problem. The structural, actuator, and control loop frequency characteristics of industrial grade robots generally allow operation with resolved rate control; a rate command is achievable with good accuracy. However, space based manipulators are different, typically have less structural stiffness, more motor and joint friction, and lower control loop cycle frequencies. These undesirable characteristics present a considerable Point of Resolution (POR) control problem for space based, kinematically redundant manipulators for the following reason: a kinematically redundant manipulator requires an arbitrary constraint to solve for the joint rate commands. A space manipulator will not respond to joint rate commands because of these characteristics. A space based manipulator simulation, including free end rigid body dynamics, motor dynamics, motor striction/friction, gearbox backlash, joint striction/friction, and Space Station Remote Manipulator System type configuration parameters, is used to evaluate the performance of a documented resolved rate control law. Alternate schemes which include torque control are also evaluated.

Robert W Bailey↗
DOE OSTI2026

Physicochemical and Performance Characterization of Six Commercial Organic Solvent Nanofiltration Membranes

This work introduces a novel, gradient-free metamaterial design method based on Gaussian process regression to represent the density field of a unit cell. The dimension of the design space is determined by the covariance matrix dimension in the Gaussian process regression. We propose compressing this matrix using an autoencoder, enabling the decoder to generate the density field and effectively reduce the originally large design space to a lower-dimensional subspace. In this compressed space, we employ an active learning method, Bayesian Adaptive Direct Search (BADS), for efficient exploration of the design space. We demonstrate that for simple 2D designs aimed at maximizing unit cell stiffness, our method yields results comparable to those of standard topology optimization. Furthermore, we extend our approach to various mechanical problems, from linear elasticity to hyperelastic large deformation and elasto-plasticity under finite deformation, to 3D metamaterial design. This illustrates the method’s versatility and effectiveness across a range of applications.

Wu, Haoran↗
NASA NTRS1962

SOME SHELL STABILITY PROBLEMS IN MISSILE AND SPACE VEHICLE ANALYSIS

A brief discussion is presented of three structural shell stability studies in progress at Lockheed Missiles & Space Company: the stabilizing influence of solid propellants, the effect of cushion stiffness on buckling under cushion loading, and the snap-through buckling of axially loaded cylinders. Recent analytical results are compared with test data. Also included is a brief description of a few shell stability problems of present interest to LMSC for which satisfactory methods of analysis are not available.

MISSILE↗
NASA NTRS1986

Physical and numerical sources of computational inefficiency in integration of chemical kinetic rate equations: Etiology, treatment and prognosis

The design of a very fast, automatic black-box code for homogeneous, gas-phase chemical kinetics problems requires an understanding of the physical and numerical sources of computational inefficiency. Some major sources reviewed in this report are stiffness of the governing ordinary differential equations (ODE's) and its detection, choice of appropriate method (i.e., integration algorithm plus step-size control strategy), nonphysical initial conditions, and too frequent evaluation of thermochemical and kinetic properties. Specific techniques are recommended (and some advised against) for improving or overcoming the identified problem areas. It is argued that, because reactive species increase exponentially with time during induction, and all species exhibit asymptotic, exponential decay with time during equilibration, exponential-fitted integration algorithms are inherently more accurate for kinetics modeling than classical, polynomial-interpolant methods for the same computational work. But current codes using the exponential-fitted method lack the sophisticated stepsize-control logic of existing black-box ODE solver codes, such as EPISODE and LSODE. The ultimate chemical kinetics code does not exist yet, but the general characteristics of such a code are becoming apparent.

Pratt, D. T.↗
NASA NTRS1989

Parallel linear equation solvers for finite element computations

The overall objective of this research is to develop efficient methods for the solution of linear and nonlinear systems of equations on parallel and supercomputers, and to apply these methods to the solution of problems in structural analysis. Attention has been given so far only to linear equations. The methods considered for the solution of the stiffness equation Kx=f have been Choleski factorization and the conjugate gradient iteration with SSOR and Incomplete Choleski preconditioning. More detail on these methods will be given on subsequent slides. These methods have been used to solve for the static displacements for the mast and panel focus problems in conjunction with the CSM testbed system based on NICE/SPAR.

Ortega, James M.↗
NASA NTRS1972

Nonapplicability of linear finite element programs to the stress analysis of tires

A static finite element stress analysis of an inflated radial car tire was carried out. The deformed shape of the sidewall presents outward bulging. The analysis of a homogeneous isotropic toroidal shell shows that the problem is common to all solids of this type. The study suggests that the geometric stiffness due to the inflation pressure has to be taken into account. Also, the resulting large displacements make it necessary for the geometry to be updated at each load step.

Durand, M.↗
NASA NTRSDate not supplied

Error Localization Examples: Looking for a Needle in a Hay-stack

Finite element models (FEM) are routinely developed and used during fabrication of high dollar-value hardware. NASA as part of the pre-flight certification of launch vehicles routinely conducts vibration and static tests to calibrate models used for flight-risk assessments. As part of the calibration process, certain areas in the model are modified, using engineering judgment and sensitivity analysis, to match the test results. Unfortunately, tools to identify problem areas in the FEM using test data directly are scarce. Over the years, Error Localization Algorithms (ELA) have been proposed with very limited success. Recently, the Analytical Dynamics Model Improvement (ADMI) algorithm, which computes closed-form mass and stiffness corrections to match the test data exactly, have been shown to be effective for error localization. The paper will present several FEM example problems where ELA is used with simulated test data to determine FEM problem areas. For each example, the correct answer is shown along with ELA results. It is shown that the ELA process is able to identify general problem areas in the FEM, which are consistent with known model perturbations. However, in most cases the ELA identified area of improvement is larger than the true answer. Nonetheless, with proper optimization tools, calibration results using the ELA identified areas provide excellent results.

error localization↗
NASA NTRS1993

The effect of nonsymmetric pressure stiffness on the dynamic characteristics of Solid Rocket Motors

This paper discusses the effect of pressure on the dynamics of pre-stiffened structures such as the Advanced Solid Rocket Motor (ASRM). Previous work in which the stiffness terms resulting from constant pressure were derived has been extended to enable modeling of nonconstant pressure applied over nonenclosed volumes. These conditions will result in nonsymmetric terms in the global stiffness matrix which will not cancel out. Three new pressure stiffness elements incorporating these nonsymmetric terms have been implemented as dummy elements in COSMIC NASTRAN and have been tested on various simple examples as well as an existing ASRM NASTRAN finite element model. The results indicate that for all load cases of practical interest to the ASRM program, the nonsymmetric terms have very little effect on the dynamic characteristics. In addition, the pressure stiffness elements developed in the previous work which assumed constant pressure gave virtually the same results as the new elements even for problems in which the pressures are not constant. The original elements appear to work well as long as the pressure gradient across any individual element is no larger than about 0.75 psi/inch. The new elements are therefore most useful for determining the conditions under which the original pressure stiffness elements can be used.

Christensen, Eric R.↗
NASA NTRS1983

On a finite-difference method for solving transient viscous flow problems

A method has been developed to solve the unsteady, compressible Navier-Stokes equation with the property of consistency and the ability of minimizing the equation stiffness. It relies on innovative extensions of the state-of-the-art finite-difference techniques and is composed of: (1) the upwind scheme for split-flux and the central scheme for conventional flux terms in the inviscid and viscous regions, respectively; (2) the characteristic treatment of both inviscid and viscous boundaries; (3) an ADI procedure compatible with interior and boundary points; and (4) a scalar matrix coefficient including viscous terms. The performance of this method is assessed with four sample problems; namely, a standing shock in the Laval duct, a shock reflected from the wall, the shock-induced boundary-layer separation, and a transient internal nozzle flow. The results from the present method, an existing hybrid block method, and a well-known two-step explicit method are compared and discussed. It is concluded that this method has an optimal trade-off between the solution accuracy and computational economy, and other desirable properties for analyzing transient viscous flow problems.

Li, C. P.↗
NASA NTRS1988

Extension and validation of a method for locating damaged members in large space trusses

The damage location approach employs the control system capabilities for the structure to test the structure and measure the dynamic response. The measurements are then used in a system identification algorithm to produce a model of the damaged structure. The model is compared to one for the undamaged structure to find regions of reduced stiffness which indicate the location of damage. Kabe's 3,4 stiffness matrix adjustment method was the central identification algorithm. The strength of his method is that, with minimal data, it preserves the representation of the physical connectivity of the structure in the resulting model of the damaged truss. However, extensive storage and computational effort were required as a result. Extension of the damage location method to overcome these problems is the first part of the current work. The central system identification algorithm is replaced with the MSMT method of stiffness matrix adjustment which was previously derived by generalizing an optimal-update secant method form quasi-Newton approaches for nonlinear optimization. Validation of the extended damage location method is the second goal.

Smith, Suzanne Weaver↗
NASA NTRS1957

Torsional stiffness of thin-walled shells having reinforcing cores and rectangular, triangular, or diamond cross section

A theoretical investigation has been made of the Saint-Venant torsion of certain composite bars. These bars are composed of two materials -- one material in the form of a thin-walled cylindrical shell and the other material in the form of a core which fills the interior of the shell and is bonded to it. An approximate boundary-value problem is formulated on assumptions similar to those of the theory of torsion of hollow thin-walled shells (Bredt theory). This boundary-value problem is solved exactly for a rectangular cross section and approximately for slender triangular and diamond cross sections. Results for the torsional stiffness constants are presented graphically.

Mccomb, Harvey G , Jr↗
NASA NTRS1985

Relative displacement measurements for two-body problems

While high-sensitivity moire interferometry is very effective for displacement measurements in continuous bodies, special difficulties arise with two-body problems. Fringes may become discontinuous at the boundary between the bodies and as a result the relative displacements cannot be extracted from the fringe patterns. In this work, flexible bridges are used between the bodies to provide a continuous path for fringe counting. The bridges are made of relatively low modulus material and have a low stiffness geometry, such that they deform under minimal load. The experimental methods used for obtaining relative displacement measurements in two-body problems are described for the two cases; dovetail joint in turbo-machinery, and thick adherend lap joint.

Post, D.↗
NASA NTRS1985

An Overview of Latest Model Reduction and Control Methods of Large Flexible Space Structures

The latest trends and theoretical developments involved with the modeling and control of Large Flexible Space Structures (LFSS) are described. The paper addresses first the basic problems, characteristics, and difficulties inherent in modeling and control of LFSS. Major sources of difficulties and errors are the stiffness and damping operators of the dynamic model. Extensions of Linear Quadratic Gaussian (LQG) theory as applied to LFSS are presented, including frequency-shaped cost functionals and perturbation methods. The minimum data/maximum entropy approach which uses a stochastic design model to overcome difficulties found in the LQG-based methods is described. Latest trends in system theory including balanced realization and singular-value analysis are used to determine reduced order controllers and models. Ad hoc methods such as component cost analysis and modal cost analysis are discussed in context with the closed-loop reduction problem of controller order versus performance. The minimum data/maximum entropy approach also addresses controller order versus performance. Those areas of control science and large scale systems that appear to have an important role in understanding and solving LFSS modeling and control are also identified.

Santiago, J. M.↗
NASA NTRS1973

NASTRAN buckling study of a linear induction motor reaction rail

NASTRAN was used to study problems associated with the installation of a linear induction motor reaction rail test track. Specific problems studied include determination of the critical axial compressive buckling stress and establishment of the lateral stiffness of the reaction rail under combined loads. NASTRAN results were compared with experimentally obtained values and satisfactory agreement was obtained. The reaction rail was found to buckle at an axial compressive stress of 11,400 pounds per square inch. The results of this investigation were used to select procedures for installation of the reaction rail.

Williams, J. G.↗
NASA NTRS1975

Mixed isoparametric elements for Saint-Venant torsion

Mixed isoparametric elements are presented for the Saint-Venant torsion problem of laminated and anisotropic bars. Both triangular and quadrilateral elements are considered. The 'generalized' element stiffness matrix is obtained by using a modified form of the Hellinger-Reissner mixed variational principle. Group-theoretic techniques are used in conjunction with computerized symbolic integration to obtain analytic expressions for the stiffness coefficients. The accuracy of the mixed isoparametric elements developed is demonstrated by means of numerical examples, and their advantages over commonly used stress and displacement elements are discussed.

Noor, A. K.↗
NASA NTRS1977

Recent advances in aerodynamic energy concept for flutter suppression and gust alleviation using active controls

Control laws are derived, by using realizable transfer functions, which permit relaxation of the stability requirements of the aerodynamic energy concept. The resulting aerodynamic eigenvalues indicate that both the trailing edge and the leading edge-trailing edge control systems can be made more effective. These control laws permit the introduction of aerodynamic damping and stiffness terms in accordance with the requirements of any specific system. Flutter suppression and gust alleviation problems can now be treated by either a trailing edge control system or by a leading edge-trailing edge control system by using the aerodynamic energy concept. Results are applicable to a wide class of aircraft operating at subsonic Mach numbers.

Nissim, E.↗
NASA NTRS1986

Implicit methods for computing chemically reacting flow

The backward Euler scheme was used to solve a large system of inviscid flow and chemical rate equations in three spatial coordinates. The flow equations were integrated simultaneously in time by a conventional ADI factorization technique, then the species equations were solved by either simultaneous or successive techniques. The methods were evaluated in their efficiency and robustness for a hypersonic flow problem involving an aerobrake configuration. It was found that both implicit methods can effectively reduce the stiffness associated with the chemical production term and that the successive solution for the species was as stable as the simultaneous solution. The latter method is more economical because the computation time varies linearly with the number of species.

Li, C. P.↗