SEARCH · Search NASA

Results for “Stiff problems”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 685 records · Page 38

All sources NASA NTRS DOE OSTI

NASA NTRS1960

Some Divergence Characteristics of Low-Aspect-Ratio Wings at Transonic and Supersonic Speeds

The problem of chordwise, or camber, divergence at transonic and supersonic speeds is treated with primary emphasis on slender delta wings having a cantilever support at the trailing edge. Experimental and analytical results are presented for four wing models having apex half-angles of 5 deg, 10 deg, 15 deg, and 20 deg. A Mach number range from 0.8 to 7.3 is covered. The analytical results include calculations based on small-aspect-ratio theory, lifting-surface theory, and strip theory. A closed-form solution of the equilibrium equation is given, which is based on low-aspect-ratio theory but which applies only to certain stiffness distributions. Also presented is an iterative procedure for use with other aerodynamic theories and with arbitrary stiffness distribution.

Woolston, Donald S.↗

DOE OSTI2022

Stress-constrained topology optimization of lattice-like structures using component-wise reduced order models

We report lattice-like structures can provide a combination of high stiffness with light weight that is useful in many applications, but a resolved finite element mesh of such structures results in a computationally expensive discretization. This computational expense may be particularly burdensome in many-query applications, such as optimization. We develop a stress-constrained topology optimization method for lattice-like structures that uses component-wise reduced order models as a cheap surrogate, providing accurate computation of stress fields while greatly reducing run time relative to a full order model. We demonstrate the ability of our method to produce large reductions in mass while respecting a constraint on the maximum stress in a pair of test problems. The ROM methodology provides a speedup of about 150x in forward solves compared to full order static condensation and provides a relative error of less than 5% in the relaxed stress.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1993

Mass-additive modal test method for verification of constrained structural models

A method for deriving constrained or fixed-base modes and frequencies from free-free modes of a structure with mass-loaded boundaries is developed. Problems associated with design and development of test fixtures can be avoided with such an approach. The analytical methodology presented is used to assess applicability of the mass-additive method for three types of structures and to determine the accuracy of derived constrained modes and frequencies. Results show that mass loading of the boundaries enables local interface modes to be measured within a desired frequency bandwidth, thus allowing constrained modes to be derived with considerably fewer free-free modes than for unloaded boundaries. Good convergence was obtained for a simple beam and a truss-like Shuttle payload, both of which had well-spaced modes and stiff interface support structures. Slow convergence was obtained for a space station module prototype, a shell-like structure having high modal density.

Admire, John R.↗

NASA NTRS1965

The Dynamics of Solid Propellant Rocket Motors

From the point of view of dynamics, a solid propellant rocket motor is a unique structure in that it is composed of a substantial mass of propellant material case-bonded to a relatively massless, thin-walled cylinder. The mechanical properties of the propellant are such that it contributes little to the stiffness of the composite structure but it does contribute to the dynamic characteristics of the structure because of its mass. Furthermore, due to the viscoelastic character of the propellant, it can be expected to provide considerable damping to the system. At this point in the development of the art of design of solid propellant rocket motors there are no clear-cut or well-founded methods to evaluate quantitatively the contributions of the propellant to the dynamic responses of the composite structure. It is clear that unless such methods are devised, it will be difficult to arrive at accurate missile designs. In this paper a survey of the dynamic problems of solid propellant rocket motors is presented starting from the simplest model thereof and proceeding, step-by-step, to the consideration of more sophisticated and realistic models. The consideration is restricted to infinitesimal deformations of propellant grains with linear mechanical properties. Substantial progress has been achieved towards the solution of many important dynamical problems. We shall attempt to summarize the pertinent developments and indicate along which lines we feel future study should proceed. A few illustrative solutions are included in areas wherein progress has been substantial.

DYNAMIC RESPONSE↗

NASA NTRS1981

Graph theory approach to the eigenvalue problem of large space structures

Graph theory is used to obtain numerical solutions to eigenvalue problems of large space structures (LSS) characterized by a state vector of large dimensions. The LSS are considered as large, flexible systems requiring both orientation and surface shape control. Graphic interpretation of the determinant of a matrix is employed to reduce a higher dimensional matrix into combinations of smaller dimensional sub-matrices. The reduction is implemented by means of a Boolean equivalent of the original matrices formulated to obtain smaller dimensional equivalents of the original numerical matrix. Computation time becomes less and more accurate solutions are possible. An example is provided in the form of a free-free square plate. Linearized system equations and numerical values of a stiffness matrix are presented, featuring a state vector with 16 components.

Reddy, A. S. S. R.↗

DOE OSTI2021

Learning constitutive relations using symmetric positive definite neural networks

In this work, we present a new neural-network architecture, called the Cholesky-factored symmetric positive definite neural network (SPD-NN), for modeling constitutive relations in computational mechanics. Instead of directly predicting the stress of the material, the SPD-NN trains a neural network to predict the Cholesky factor of the tangent stiffness matrix, based on which the stress is calculated in incremental form. As a result of this special structure, SPD-NN weakly imposes convexity on the strain energy function, satisfies the second order work criterion (Hill's criterion) and time consistency for path-dependent materials, and therefore improves numerical stability, especially when the SPD-NN is used in finite element simulations. Depending on the types of available data, we propose two training methods, namely direct training for strain and stress pairs and indirect training for loads and displacement pairs. We demonstrate the effectiveness of SPD-NN on hyperelastic, elasto-plastic, and multiscale fiber-reinforced plate problems from solid mechanics. The generality and robustness of SPD-NN make it a promising tool for a wide range of constitutive modeling applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2016

Bifurcation, chaos, and scan instability in dynamic atomic force microscopy

The dynamical motion at any point on the cantilever of an atomic force microscope can be expressed quite generally as a superposition of simple harmonic oscillators corresponding to the vibrational modes allowed by the cantilever shape. Central to the dynamical equations is the representation of the cantilever-sample interaction force as a polynomial expansion with coefficients that account for the interaction force “stiffness,” the cantilever-to-sample energy transfer, and the displacement amplitude of cantilever oscillation. Renormalization of the cantilever beam model shows that for a given cantilever drive frequency cantilever dynamics can be accurately represented by a single nonlinear mass-spring model with frequency-dependent stiffness and damping coefficients [S. A. Cantrell and J. H. Cantrell, J. Appl. Phys. 110, 094314 (2011)]. Application of the Melnikov method to the renormalized dynamical equation is shown to predict a cascade of period doubling bifurcations with increasing cantilever drive force that terminates in chaos. The threshold value of the drive force necessary to initiate bifurcation is shown to depend strongly on the cantilever setpoint and drive frequency, effective damping coefficient, nonlinearity of the cantilever-sample interaction force, and the displacement amplitude of cantilever oscillation. The model predicts the experimentally observed interruptions of the bifurcation cascade for cantilevers of sufficiently large stiffness. Operational factors leading to the loss of image quality in dynamic atomic force microscopy are addressed, and guidelines for optimizing scan stability are proposed using a quantitative analysis based on system dynamical parameters and choice of feedback loop parameter.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS2019

Toward Ultralight High Strength Structural Materials via Collapsed Carbon Nanotube Bonding

The growing commercial availability of carbon nanotube (CNT) macro-assemblies such as sheet and yarn is making their use in structural composite components increasingly feasible. However, the mechanical properties of these materials continue to trail those of state-of-the-art carbon fiber composites due to relatively weak inter-tube load transfer. Forming covalent links between adjacent CNTs promises to mitigate this problem, but it has proven difficult in practice to introduce them chemically within densified and aligned CNT materials due to their low permeability. To avoid this limitation, this work explores the combination of pulsed electrical current, temperature, and pressure to introduce inter-CNT bonds. Reactive molecular dynamics simulations identify the most probable locations, configurations, and conditions for inter-nanotube bonds to form. This process is shown to introduce covalent linkages within the CNT material that manifest as improved macroscale mechanical properties. The magnitude of this effect increases with increasing levels of prealignment of the CNT material, promising a new synthesis pathway to ultralight structural materials with specific strengths and stiffnesses exceeding 1 and 100 GPa/(g/cu.cm), respectively.

Jensen, Benjamin D.↗

NASA NTRS1995

Parametric studies of stitching effectiveness for preventing substructure disbond

A methodology is desired that will allow a designer to select appropriate amounts of through-thickness reinforcement needed to meet design requirements. The goal is to use a relatively simple analysis to minimize the amount of testing that needs to be performed, and to make test results from simple configurations applicable to more general structures. Using this methodology one should be able to optimize the selection of stitching materials, the weight of the yarn, and the stitching density. The analysis approach is to treat substructure disbond as a crack propagation problem. In this approach, the stitches have little influence until a delamination begins to grow. Once the delamination reaches, or extends beyond a stitch, the stitch serves to reduce the strain-energy-release-rate (G) at the crack tip for a given applied load. The reduced G can then be compared to the unstitched materials toughness to predict the load required to further extend the crack. The current model treats the stitch as a simple spring which responds to displacements in the vertical (through-thickness) direction. In concept, this approach is similar to that proposed by other authors. Test results indicate that the model should be refined to include the shearing stiffness of the stitch. The strain-energy-release-rate calculations are performed using a code which uses interconnected higher-order plates to model built-up composite cross-sections. When plates are stacked vertically, the interfacial tractions between the plates can be computed. The plate differential equations are solved in closed-form. The code, called SUBLAM, was developed as part of this section in one dimension. Because of this limitation, rows of stitches are treated as a two-dimensional sheet. The spring stiffness of a row of stitches can be estimated from the stitch material, weight, and density. As a practical and conservative approach, we can assume that the stitch is bonded until a crack passes the stitch location. After the crack passes, it is fully bonded. A series of tests were performed to exercise this methodology and incorporated an attached flange such that the sudden change in thickness initiated a delamination. The analysis was used to estimate the materials' critical G from that of the unstitched specimens. With this data, a prediction was made for the load required to delaminate the stitched specimens. Using the methodology, design charts have been created for simplified geometries. These charts give stitch force, along with G(sub 1) and G(sub 2) as as function of the stitch spring stiffness. Using the charts, it should be possible to determine the stitch spring stiffness and strength required to reduce the G to a desired level. From these parameters, the actual stitching material, weight, and density can be computed.

Flanagan, Gerry↗

DOE OSTI2024

Robust 3D multi-material hydrodynamics using discontinuous Galerkin methods

A high-order discontinuous Galerkin (DG) method is presented for nonequilibrium multi-material (m ≥ 2) flow with sharp interfaces. Material interfaces are reconstructed using the algebraic THINC approach, resulting in a sharp interface resolution. The system assumes stiff velocity relaxation and pressure nonequilibrium. The presented DG method uses Dubiner's orthogonal basis functions on tetrahedral elements. This results in a unique combination of sharp multimaterial interfaces and high-order accurate solutions in smooth single-material regions. A novel shock indicator based on the interface conservation condition is introduced to mark regions with discontinuities. Slope limiting techniques are applied only in these regions so that nonphysical oscillations are eliminated while maintaining high-order accuracy in smooth regions. A local projection is applied on the limited solution to ensure discrete closure law preservation. The effectiveness of this novel limiting strategy is demonstrated for complex three-dimensional multi-material problems, where robustness of the method is critical. The presented numerical problems demonstrate that more accurate and efficient multi-material solutions can be obtained by the DG method, as compared to second-order finite volume methods.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Symbolic construction of the chemical Jacobian of quasi-steady state (QSS) chemistries for Exascale computing platforms

The Quasi-Steady State Approximation (QSSA) can be an effective tool for reducing the size and stiffness of chemical mechanisms for implementation in computational reacting flow solvers. However, for many applications, the resulting model still requires implicit methods for efficient time integration. Here, in this paper, we outline an approach to formulating the QSSA reduction that is coupled with a strategy to generate C++ source code to evaluate the net species production rates, and the chemical Jacobian. The code-generation component employs a symbolic approach enabling a simple and effective strategy to analytically compute the chemical Jacobian. For computational tractability, the symbolic approach needs to be paired with common subexpression elimination which can negatively affect memory usage. Several solutions are outlined and successfully tested on a 3D multipulse ignition problem, thus allowing portable application across chemical model sizes and GPU capabilities. The implementation of the proposed method is available at https://github.com/AMReX-Combustion/PelePhysics under an open-source license.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1978

Solution of the subsynchronous whirl problem in the high-pressure hydrogen turbomachinery of the Space Shuttle Main Engine

Subsynchronous whirl of the high-pressure fuel turbopump limited operation of the Space Shuttle Main Engine for some months in early 1976. The means by which this problem was successfully attacked is of particular interest to the rotor-dynamics community, not only because this machinery was the highest power-to-weight ratio known (77,000 hp...760 pounds), but because of the multiple forcing functions involved and the means, both analytical and experimental, which were utilized in separating variables, pointing toward successful solutions, and evaluating results. The general means of identifying fundamental characteristics, analyzing data, and conducting computer investigations are delineated. The results of analysis and testing are discussed. Since whirl inception occurs at a shaft speed greater than twice the first system critical, this was increased by stiffening the shaft and bearing supports; the system damping and system stiffness was additionally increased by proper interstage seal design to the extent that the instability threshold is now beyond the operating range.

Ek, M. C.↗

NASA NTRS1989

A magnetic damper for first mode vibration reduction in multimass flexible rotors

Many rotating machines such as compressors, turbines and pumps have long thin shafts with resulting vibration problems, and would benefit from additional damping near the center of the shaft. Magnetic dampers have the potential to be employed in these machines because they can operate in the working fluid environment unlike conventional bearings. An experimental test rig is described which was set up with a long thin shaft and several masses to represent a flexible shaft machine. An active magnetic damper was placed in three locations: near the midspan, near one end disk, and close to the bearing. With typical control parameter settings, the midspan location reduced the first mode vibration 82 percent, the disk location reduced it 75 percent and the bearing location attained a 74 percent reduction. Magnetic damper stiffness and damping values used to obtain these reductions were only a few percent of the bearing stiffness and damping values. A theoretical model of both the rotor and the damper was developed and compared to the measured results. The agreement was good.

Kasarda, M. E. F.↗

DOE OSTI2021

Component-wise reduced order model lattice-type structure design

Lattice-type structures can provide a combination of stiffness with light weight that is desirable in a variety of applications. Design optimization of these structures must rely on approximations of the governing physics to render solution of a mathematical model feasible. In this paper, we propose a topology optimization (TO) formulation that approximates the governing physics using component-wise reduced order modeling as introduced in Huynh et al. (2013); Eftang and Patera (2013), which can reduce solution time by multiple orders of magnitude over a full-order finite element model while providing a relative error in the solution of 1%. In addition, the offline training data set from such component-wise models is reusable, allowing its application to many design problems for only the cost of a single offline training phase, and the component-wise method is nearly embarrassingly parallel. We also show how the parameterization chosen in our optimization allows a simplification of the component-wise reduced order model (CWROM) not noted in previous literature, for further speedup of the optimization process. Furthermore, the sensitivity of the compliance with respect to the particular parameterization is derived solely at the component level. In numerical examples, we demonstrate a 1000x speedup over a full-order FEM model with relative error of 1% and show minimum compliance designs for two different cantilever beam examples, one smaller and one larger. Finally, error bounds for the displacement field, compliance, and compliance sensitivity of the CWROM are derived.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1991

On the application of ENO scheme with subcell resolution to conservation laws with stiff source terms

Two approaches are used to extend the essentially non-oscillatory (ENO) schemes to treat conservation laws with stiff source terms. One approach is the application of the Strang time-splitting method. Here the basic ENO scheme and the Harten modification using subcell resolution (SR), ENO/SR scheme, are extended this way. The other approach is a direct method and a modification of the ENO/SR. Here the technique of ENO reconstruction with subcell resolution is used to locate the discontinuity within a cell and the time evolution is then accomplished by solving the differential equation along characteristics locally and advancing in the characteristic direction. This scheme is denoted ENO/SRCD (subcell resolution - characteristic direction). All the schemes are tested on the equation of LeVeque and Yee (NASA-TM-100075, 1988) modeling reacting flow problems. Numerical results show that these schemes handle this intriguing model problem very well, especially with ENO/SRCD which produces perfect resolution at the discontinuity.

Chang, Shih-Hung↗

NASA NTRS2005

Tools for Designing and Analyzing Structures

Structural Design and Analysis Toolset is a collection of approximately 26 Microsoft Excel spreadsheet programs, each of which performs calculations within a different subdiscipline of structural design and analysis. These programs present input and output data in user-friendly, menu-driven formats. Although these programs cannot solve complex cases like those treated by larger finite element codes, these programs do yield quick solutions to numerous common problems more rapidly than the finite element codes, thereby making it possible to quickly perform multiple preliminary analyses - e.g., to establish approximate limits prior to detailed analyses by the larger finite element codes. These programs perform different types of calculations, as follows: 1. determination of geometric properties for a variety of standard structural components; 2. analysis of static, vibrational, and thermal- gradient loads and deflections in certain structures (mostly beams and, in the case of thermal-gradients, mirrors); 3. kinetic energies of fans; 4. detailed analysis of stress and buckling in beams, plates, columns, and a variety of shell structures; and 5. temperature dependent properties of materials, including figures of merit that characterize strength, stiffness, and deformation response to thermal gradients

Luz, Paul L.↗

NASA NTRS1985

Rotor-to-stator rub vibration in centrifugal compressor

One example of excessive vibration encountered during loading of a centrifugal compressor train (H type compressor with HP casing) is discussed. An investigation was made of the effects of the dynamic load on the bearing stiffness and the rotor-bearing system critical speed. The high vibration occurred at a "threshold load," but the machine didn't run smoothly due to rubs even when it had passed through the threshold load. The acquisition and discussion of the data taken in the field as well as a description of the case history which utilizes background information to identify the malfunction conditions is presented. The analysis shows that the failures, including full reverse precession rub and exact one half subharmonic vibration, were caused by the oversize bearings and displacement of the rotor center due to foundation deformation and misalignment between gear shafts, etc. The corrective actions taken to alleviate excessive vibration and the problems which remain to be solved are also presented.

Gao, J. J.↗

NASA NTRS2022

Mars 2020 Perseverance Rover SHERLOC Instrument Isolation System

The NASA Jet Propulsion Laboratory (JPL) successfully landed the Mars 2020 Perseverance Rover at Jezero Crater on the Martian surface. Perseverance’s main mission objective is to cache Martian rock and regolith samples in hermetically sealed tubes to be brought back to Earth by future missions. To achieve this goal the rover is equipped with a 2-meter-long Robotic Arm (RA) which manipulates the Turret to interact with the surface. The Turret is comprised of science instruments and tools that enable surface sampling and science. The backbone of the Turret hardware is the Rotary Percussive Coring Drill. The science instruments, which are directly mounted to the structural housing of this drill, must be able to withstand not only the launch and entry, descent and landing (EDL) loads of the vehicle but also the dynamic percussive environment, induced by the drill throughout surface operations. This paper discusses the technical hardware design, development and testing of a vibration isolation system to protect the Scanning Habitable Environments with Raman and Luminescence for Organics and Chemicals (SHERLOC) Instrument from these dynamic environments. The hardware design aspect of this problem was limited by both size and mass constraints. The Turret is tightly configured and the physical space in which hardware can reside is restricted due to multiple axes of surface interaction points. Another driving requirement is the extreme range of thermal non-operational environment from -135C to +90C. This, in addition to SHERLOC being a heavier instrument than previously integrated on a turret, made the use of isolation systems employed on previous rover missions problematic. The final design incorporates custom wire mesh springs preloaded within titanium hexapod struts equipped with flexure ends. Custom versions of commercially available wire mesh springs were designed and tested. These custom springs differed in dimension, density and stiffness from the off the shelf options. The newly designed springs went through various stages of flight acceptance testing. Detailed physics-based modeling to understand what loads the isolation system would experience was used to define the test conditions for the springs and entire isolation subsystem. The testing program developed around this was extensive to fully characterize the isolation performance over cleanliness levels, temperature and life. Another major aspect of this work was developing a method of replicating the drill’s percussive environment for performance and life testing of Turret mounted hardware. Attaching the instrument to an operating drill for testing was not a viable option. A method of using a Highly Accelerated Life Test (HALT) table was developed and applied across the project as the most representative and tunable option to physically simulate this new environment. This method, along with launch environment random vibration testing, was used to build an entire dynamic qualification and flight hardware test program. The results of the design, modeling, testing and analysis of the SHERLOC Isolation System is described throughout the following paper.

Krafchak, Todd↗