SEARCH · Engineering Papers

Results for “Stiffness removal”

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 19 records

All sources NASA NTRS DOE OSTI

DOE OSTI2025

An iterative dynamic chemical stiffness removal method for reacting flow simulations

Abstract An iterative dynamic chemical stiffness removal method (IDCSR) based on quasi-steady-state approximation (QSSA) is proposed. The IDCSR method is built on a previously developed non-iterative method which has proved to work well for small timestep sizes. A novel iterative procedure is designed in IDCSR to enable explicit time integration of stiff chemistry at relatively large timestep sizes relevant to practical reacting flow simulations. The effectiveness of the iterative procedure is first demonstrated with a toy problem and homogeneous auto-ignition with fixed integration step sizes, showing that larger timestep sizes can be allowed for explicit time integration using IDCSR compared with the previous non-iterative method. IDCSR is then compared with existing explicit chemistry solvers for simulations of homogeneous auto-ignition and shows similar or lower computational cost but significantly higher accuracy across a wide range of timestep sizes. IDCSR is further combined with an automatic adaptive time-stepping scheme for simulations of 0-D homogeneous auto-ignition and a 2-D laminar lifted n -dodecane jet flame. For the 0-D auto-ignition simulations, IDCSR is shown to reduce both the error (by 43%–90%) and computational cost (by 6–15 times) compared with existing explicit solvers, while achieving speed-up factors of up to 400 compared with VODE for a wide range of timestep sizes and reaction mechanisms. For the 2-D jet flame simulations, speed-up factors of 15 and 31 for chemistry integration, and 5 and 9 for overall simulation, are achieved by IDCSR compared with CVODE with and without analytic Jacobian, respectively.

Xu, Chao (ORCID:0000000153074159)↗

NASA NTRS1989

Design of optimally smoothing multi-stage schemes for the Euler equations

In this paper, a method is developed for designing multi-stage schemes that give optimal damping of high-frequencies for a given spatial-differencing operator. The objective of the method is to design schemes that combine well with multi-grid acceleration. The schemes are tested on a nonlinear scalar equation, and compared to Runge-Kutta schemes with the maximum stable time-step. The optimally smoothing schemes perform better than the Runge-Kutta schemes, even on a single grid. The analysis is extended to the Euler equations in one space-dimension by use of 'characteristic time-stepping', which preconditions the equations, removing stiffness due to variations among characteristic speeds. Convergence rates independent of the number of cells in the finest grid are achieved for transonic flow with and without a shock. Characteristic time-stepping is shown to be preferable to local time-stepping, although use of the optimally damping schemes appears to enhance the performance of local time-stepping. The extension of the analysis to the two-dimensional Euler equations is hampered by the lack of a model for characteristic time-stepping in two dimensions. Some results for local time-stepping are presented.

Van Leer, Bram↗

NASA NTRS2012

Progress Toward Affordable High Fidelity Combustion Simulations Using Filtered Density Functions for Hypersonic Flows in Complex Geometries

Significant progress has been made in the development of subgrid scale (SGS) closures based on a filtered density function (FDF) for large eddy simulations (LES) of turbulent reacting flows. The FDF is the counterpart of the probability density function (PDF) method, which has proven effective in Reynolds averaged simulations (RAS). However, while systematic progress is being made advancing the FDF models for relatively simple flows and lab-scale flames, the application of these methods in complex geometries and high speed, wall-bounded flows with shocks remains a challenge. The key difficulties are the significant computational cost associated with solving the FDF transport equation and numerically stiff finite rate chemistry. For LES/FDF methods to make a more significant impact in practical applications a pragmatic approach must be taken that significantly reduces the computational cost while maintaining high modeling fidelity. An example of one such ongoing effort is at the NASA Langley Research Center, where the first generation FDF models, namely the scalar filtered mass density function (SFMDF) are being implemented into VULCAN, a production-quality RAS and LES solver widely used for design of high speed propulsion flowpaths. This effort leverages internal and external collaborations to reduce the overall computational cost of high fidelity simulations in VULCAN by: implementing high order methods that allow reduction in the total number of computational cells without loss in accuracy; implementing first generation of high fidelity scalar PDF/FDF models applicable to high-speed compressible flows; coupling RAS/PDF and LES/FDF into a hybrid framework to efficiently and accurately model the effects of combustion in the vicinity of the walls; developing efficient Lagrangian particle tracking algorithms to support robust solutions of the FDF equations for high speed flows; and utilizing finite rate chemistry parametrization, such as flamelet models, to reduce the number of transported reactive species and remove numerical stiffness. This paper briefly introduces the SFMDF model (highlighting key benefits and challenges), and discusses particle tracking for flows with shocks, the hybrid coupled RAS/PDF and LES/FDF model, flamelet generated manifolds (FGM) model, and the Irregularly Portioned Lagrangian Monte Carlo Finite Difference (IPLMCFD) methodology for scalable simulation of high-speed reacting compressible flows.

Drozda, Tomasz G.↗

DOE OSTI2026

Time-dependent-bases with local CUR decomposition method for accelerating turbulent combustion simulations

Here, this study presents a novel reduced-order modeling framework, Time-Dependent Bases with Local CUR decomposition (TDB-L-CUR), designed to efficiently and accurately approximate the species transport equations in reacting flow simulations. The method extends the existing TDB-CUR approach for chemically reacting flows (Jung et al. Comput. Methods Appl. Mech. Engrg. 437 (2025) 117758), which leverages matrix decomposition techniques to form a global-in-space, time-dependent low-dimensional manifold. While TDB-CUR performs well in homogeneous systems, it may be less well-suited to spatially heterogeneous systems such as turbulent flames, where higher-rank approximations are typically required. The proposed TDB-L-CUR framework introduces two methodological extensions to the baseline approach. First, it applies unsupervised clustering to partition the physical domain into distinct regions, enabling spatially localized manifold construction, thereby reducing the rank required for the reduced-order representation. Second, it incorporates a computational singular perturbation (CSP)-based scheme for identifying and penalizing fast species, allowing for spatio-temporally adaptive mitigation of chemical stiffness. The proposed framework is validated on a hierarchy of test cases, including a one-dimensional premixed flame, a two-dimensional nonpremixed ignition case with vortex interaction, and a three-dimensional turbulent premixed flame. TDB-L-CUR significantly improves accuracy over TDB-CUR while further reducing computational cost. The fully on-the-fly formulation of TDB-L-CUR (i.e., requiring no offline training or prior knowledge) makes it a robust and scalable tool for reduced-order modeling of reactive flows.

Local manifold↗

NASA NTRS1986

A new approximate LU factorization scheme for the Reynolds-averaged Navier-Stokes equations

A new approximate LU factorization scheme is developed to solve the steady state Reynolds-averaged Navier-Stokes (NS) equations. Central differencing is used for both implicit and explicit operator and special care is taken to obtain well-conditioned factors on the implicit side. The scheme is then analyzed and optimized according to a simple linear analysis. It is unconditionally stable for the model hyperbolic equation in both two- and three-dimensions. However, the requirement for well-conditioned factors has essentially limited the effective time step the scheme can achieve. Supersonic and transonic three-dimensinal flows past a hemisphere cylinder are computed to demonstrate the convergence characteristics of the scheme. A good convergence rate is achieved for the inviscid case. Finally, an explicit eigenvector annihilation procedure is adopted successfully to remove the stiffness caused by the fine grid spacing for viscous flows.

Pan, D.↗

NASA NTRS1987

Splitting methods for low Mach number Euler and Navier-Stokes equations

Examined are some splitting techniques for low Mach number Euler flows. Shortcomings of some of the proposed methods are pointed out and an explanation for their inadequacy suggested. A symmetric splitting for both the Euler and Navier-Stokes equations is then presented which removes the stiffness of these equations when the Mach number is small. The splitting is shown to be stable.

Abarbanel, Saul↗

NASA NTRS1989

Splitting methods for low Mach number Euler and Navier-Stokes equations

Abarbanel, Saul↗

NASA NTRS1985

Prestressed Thermal-Protection Panels

Panels held securely with minimum of mounting hardware. Each panel held in place by single screw that pulls it into flat shape from its original shallow-dish shape. Shape and prestressing make panel stiff: resists vibration and withstands large mechanical loads. Panel shape and mounting arrangement not limited to thermal-protection systems but also used on aircraft, building walls, or wherever large surfaces must be covered with stiff, flat sheets easily removed for maintenance.

Dunn, T. J.↗

NASA NTRS1992

Drilling and Digging Techniques for the Early Lunar Outpost

The theme of this workshop is lunar resource assessment. Topics include identification, quantification, and location of useful elements on and below the lunar surface. The objective of this paper is to look at another side of the issue--how to remove soil from the stiff lunar-soil matrix once useful deposits are located. The goal of this paper is to cause those who think that digging or excavating on the Moon is a trivial problem to rethink the reasons for their opinions. Another goal is to encourage them to view total reliance upon terrestrial heuristics with suspicion. This paper will focus primarily upon digging.

Boles, Walter W.↗

DOE OSTI2022

A machining digital twin for hybrid manufacturing

Hybrid manufacturing consisting of metal additively manufactured preforms and computer numerical control (CNC) machining has been established to be an effective method for high material use rates. However, hybrid manufacturing introduces unique challenges. Near-net shape designs are typically selected, which result in a smaller margin for part placement within the stock and stringent requirements for work coordinate system identification. Additionally, less stock material reduces the preform stiffness, which limits the material removal rates during machining. This paper demonstrates a digital twin for CNC machining of a wire arc additively manufactured preform that implements: 1) structured light scanning for stock model identification and tool path generation; 2) a fused filament fabrication apparatus to attach temporary fiducials and scan targets to the preform that enable coordinate system definition for both the CAM and CNC machine; 3) preform and tool tip frequency response function measurements to enable stable milling parameter selection; and 4) post-manufacturing measurements of geometry, surface finish, and structural dynamics to confirm designer intent. These efforts define key components of the machining digital twin for hybrid manufacturing.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Self-Propulsion by Directed Explosive Emulsification

An active droplet system, programmed to repeatedly move autonomously at a specific velocity in a well-defined direction, is demonstrated. Coulombic energy is stored in oversaturated interfacial assemblies of charged nanoparticle-surfactants by an applied DC electric field and can be released on demand. Spontaneous emulsification is suppressed by an increase in the stiffness of the oversaturated assemblies. Rapidly removing the field releases the stored energy in an explosive event that propels the droplet, where thousands of charged microdroplets are ballistically ejected from the surface of the parent droplet. The ejection is made directional by a symmetry breaking of the interfacial assembly, and the combined interaction force of the microdroplet plume on one side of the droplet propels the droplet distances tens of times its size, making the droplet active. The propulsion is autonomous, repeatable, and agnostic to the chemical composition of the nanoparticles. The symmetry-breaking in the nanoparticle assembly controls the microdroplet velocity and direction of propulsion. Finally, this mechanism of droplet propulsion will advance soft micro-robotics, establishes a new type of active matter, and introduces new vehicles for compartmentalized delivery.

36 MATERIALS SCIENCE↗

NASA NTRS1985

Characteristics of the Elevation Drive Suspension of the 64-meter Antennas

The elevation axis drive unit of the DSN 64-m antennas is arranged so that its output pinion can self-align to the bull gear. The design is described and conflicting desiderata are discussed. Static stability criterion for the elevation drive, derivation of strut stiffness ratios, and an outline for removing, installing, and adjusting the vertical struts are included.

Mcginness, H. D.↗

NASA NTRS1989

Model reduction in the physical coordinate system

In the dynamics modeling of a flexible structure, finite element analysis employs reduction techniques, such as Guyan's reduction, to remove some of the insignificant physical coordinates, thus producing a dynamics model that has smaller mass and stiffness matrices. But this reduction is limited in the sense that it removes certain degrees of freedom at a node points themselves in the model. From the standpoint of linear control design, the resultant model is still too large despite the reduction. Thus, some form of the model reduction is frequently used in control design by approximating a large dynamical system with a fewer number of state variables. However, a problem arises from the placement of sensors and actuators in the reduced model, because a model usually undergoes, before being reduced, some form of coordinate transformations that do not preserve the physical meanings of the states. To correct such a problem, a method is developed that expresses a reduced model in terms of a subset of the original states. The proposed method starts with a dynamic model that is originated and reduced in finite element analysis. Then the model is converted to the state space form, and reduced again by the internal balancing method. At this point, being in the balanced coordinate system, the states in the reduced model have no apparent resemblance to those of the original model. Through another coordinate transformation that is developed, however, this reduced model is expressed by a subset of the original states.

Yae, K. Harold↗

NASA NTRS1981

A fully redundant power hinge for LANDSAT-D appendages

The configuration and testing of a power driven hinge for deployment of the solar array and antenna boom for the LANDSAT-D spacecraft is discussed. The hinge is fully mechanically and electrically redundant and, thereby, can sustain a single point failure of any one motor (or its power supply), speed reducer, or bearing set without loss of its ability to function. This design utilizes the capability of the stepper motor drive to remove the flexibility of the drive train from the joint stiffness equation when the hinge is loaded against its stop. This feature precludes gapping of the joint under spacecraft maneuver loads even in the absence of a latching feature. Thus, retraction is easily accomplished by motor reversal without the need for a solenoid function to remove the latch.

Mamrol, F. E.↗

DOE OSTI2020

Eigenfrequency constrained topology optimization of finite strain hyperelastic structures

This paper incorporates hyperelastic materials, nonlinear kinematics, and preloads in eigenfrequency constrained density–based topology optimization. The formulation allows for initial finite deformations and subsequent small harmonic oscillations. The optimization problem is solved by the method of moving asymptotes, and the gradients are calculated using the adjoint method. Both simple and degenerate eigenfrequencies are considered in the sensitivity analysis. A well-posed topology optimization problem is formulated by filtering the volume fraction field. Numerical issues associated with excessive distortion and spurious eigenmodes in void regions are reduced by removing low volume fraction elements. The optimization objective is to maximize stiffness subject to a lower bound on the fundamental eigenfrequency. Numerical examples show that the eigenfrequencies drastically change with the load magnitude, and that the optimization is able to produce designs with the desired fundamental eigenfrequency.

42 ENGINEERING↗

NASA NTRS1987

Bendable Routing Tool

Tool extends routing bit into internal cavities or passages so burrs and similar defects removed. Bent so inserted through curving channels. Copper sheath gives stiffness to flexible shaft but can itself be bent. Several types of routing bits attached to shaft.

Mayer, Walter↗

NASA NTRS1993

A 69 cent look at thermoplastic softening

The objectives are to demonstrate the change in mechanical properties of a thermoplastic polymer when the polymer is heated above the glass transition temperature. The instructor or student can perform this demonstration. In preparation for the demonstration, fill the flask or beaker 3/4 full with water and bring the water to a boil. First observe the initial stiffness of the toothbrush by slightly flexing it. Note that it is fairly stiff and retains its initial shape when the forces are removed. If you were to apply enough force, the toothbrush would break in a somewhat brittle fashion. Now place the toothbrush in the boiling water for 2 minutes. After the 2 minutes are up, remove the toothbrush from the boiling water with the washcloth. At this time, you will be able to bend the toothbrush easily. When it cools it will retain this bent shape. You've just demonstrated heating the polymer above the glass transition temperature, demonstrating thermoplastic softening. Instructor notes are provided.

Vanasupa, Linda S.↗

NASA NTRS1986

Structural performance of space station trusses with missing members

Structural performance of orthogonal tetrahedral and Warren-type full truss beams and platforms are compared. In addition, degradation of truss structural performance is determined for beams, platforms and a space station when individual struts are removed from the trusses. The truss beam, space station, and truss platform analytical models used in the studies are described. Stiffness degradation of the trusses due to single strut failures is determined using flexible body vibration modes. Ease of strut replacement is assessed by removing a strut and examining the truss deflection at the resulting gap due to applied forces. Finally, the reduction in truss beam strength due to a missing longeron is determined for a space station transverse boom model.

Dorsey, J. T.↗