SEARCH · Search NASA
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.
Improved numerical procedure more than twice as fast as previous methods. Unified numerical algorithm efficiently solves free-vibration problems of stationary or spinning structures with or without viscous or structural damping. Algorithm used to solve static problems involving multiple loads and to solve quadratic matrix eigenvalue problems associated with finite-dynamic-element structural discretization.
The paper first presents the details of the development of a new six-noded plane triangular finite dynamic element. A block Lanczos algorithm is developed next for the accurate and efficient solution of the quadratic matrix eigenvalue problem associated with the finite dynamic element formulation. The resulting computer program fully exploits matrix sparsity inherent in such a discretization and proves to be most efficient for the extraction of the usually required first few roots and vectors, including repeated ones. Most importantly, the present eigenproblem solution is shown to be comparable to that of the corresponding finite element analysis, thereby rendering the associated dynamic element method rather attractive owing to superior convergence characteristics of such elements, presented herein.
The Data System Dynamic Simulator (DSDS) is a discrete event simulation tool. It was developed for NASA for the specific purpose of evaluating candidate architectures for data systems of the Space Station era. DSDS provides three methods for meeting this requirement. First, the user has access to a library of standard pre-programmed elements. These elements represent tailorable components of NASA data systems and can be connected in any logical manner. Secondly, DSDS supports the development of additional elements. This allows the more sophisticated DSDS user the option of extending the standard element set. Thirdly, DSDS supports the use of data streams simulation. Data streams is the name given to a technique that ignores packet boundaries, but is sensitive to rate changes. Because rate changes are rare compared to packet arrivals in a typical NASA data system, data stream simulations require a fraction of the CPU run time. Additionally, the data stream technique is considerably more accurate than another commonly-used optimization technique.
A surrogate model methodology is described for predicting, during flight, the residual strength of aircraft structures that sustain discrete-source damage. Starting with design of experiment, an artificial neural network is developed that takes as input discrete-source damage parameters and outputs a prediction of the structural residual strength. Target residual strength values used to train the artificial neural network are derived from 3D finite element-based fracture simulations. Two ductile fracture simulations are presented to show that crack growth and residual strength are determined more accurately in discrete-source damage cases by using an elastic-plastic fracture framework rather than a linear-elastic fracture mechanics-based method. Improving accuracy of the residual strength training data does, in turn, improve accuracy of the surrogate model. When combined, the surrogate model methodology and high fidelity fracture simulation framework provide useful tools for adaptive flight technology.
An overview of the capabilities of the CHarring Ablator Response (CHAR) code is presented. CHAR is a one-, two-, and three-dimensional unstructured continuous Galerkin finite-element heat conduction and ablation solver with both direct and inverse modes. Additionally, CHAR includes a coupled linear thermoelastic solver for determination of internal stresses induced from the temperature field and surface loading. Background on the development process, governing equations, material models, discretization techniques, and numerical methods is provided. Special focus is put on the available boundary conditions including thermochemical ablation, surface-to-surface radiation exchange, and flowfield coupling. Finally, a discussion of ongoing development efforts is presented.
An overview of the capabilities of the CHarring Ablator Response (CHAR) code is presented. CHAR is a one-, two-, and three-dimensional unstructured continuous Galerkin finite-element heat conduction and ablation solver with both direct and inverse modes. Additionally, CHAR includes a coupled linear thermoelastic solver for determination of internal stresses induced from the temperature field and surface loading. Background on the development process, governing equations, material models, discretization techniques, and numerical methods is provided. Special focus is put on the available boundary conditions including thermochemical ablation and contact interfaces, and example simulations are included. Finally, a discussion of ongoing development efforts is presented.
A mesoscale progressive damage analysis was applied to the analysis of recently tested adhesively bonded sandwich panels, in which ten nominally identical tension loaded panels failed in three different failure modes, with peakloads ranging from 13.0 –16.1 kips. A parametric finite element model for progressive damage analysis of an adhesively bonded joint (ABJ) considering matrix and fiber damage, and delamination cracks in unidirectional and fabric composite sections of the joint, as well as core crushing, and adhesive fracture was developed. The model was developed using Abaqus/Explicit and utilized continuum damage mechanics material models for intralaminar damage in the facesheet and doubler plies, cohesive elements for delamination and adhesive debonding, and discretized cells considering plasticity to model crushing of the aluminum honeycomb core. Analysis predictions obtained with a pristine model agreed well with the average experimental peak load, strain in the joint, and post-test damage states. Then, various methods to incorporate stochastic material variability into ABJ model were evaluated. These methods included simulation of pre-existing matrix damage, pre-existing adhesive porosity, and uniform modification of adhesive properties. These material and manufacturing defects changed the damage mode predictions in some cases. The analysis of stochastic strength due to the considered defects revealed unintuitive damage mode interaction in the ABJ, explaining some structure-property relationships observed experimentally.
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.
This viewgraph presentation provides information on optimizing the travel distance between two points on a curved surface. The presentation addresses the single source shortest path problem, fast algorithms for estimating the eikonal equation, fast schemes and barrier theorems, and the discontinuous Galerkin method, including hyperbolic causality, finite element method, scalars, and marching the discontinuous Galerkin Eikonal approximation.
Numerical simulations are essential to evaluate the performance and safety of engineered subsurface systems such as geological carbon storage sites, enhanced geothermal fields, and oil and gas reservoirs. A key challenge lies in accurately modeling the frictional contact behavior along fault surfaces. This problem involves inequality constraints that arise from the physics of frictional slip, requiring specialized numerical methods to handle the resulting highly nonlinear and path-dependent behavior. Here, in this work, we address this challenge using an Augmented Lagrangian Method (ALM) implemented via the Uzawa algorithm. The formulation employs mixed finite element spaces, combining low-order piecewise linear displacements within the 3D domain cells with piecewise constant tractions defined on the fault surfaces. Furthermore, to ensure stability and satisfy the inf-sup condition, the discrete displacement space is enriched with face bubble functions on both sides of the contact interfaces. This approach offers several advantages over other stabilization techniques that rely on additional terms, and it integrates naturally in the Uzawa framework.
Finite element methods are described for modeling high speed compressible flows with strong advection, problems important to aerodynamics. The situations are characterized by high pressure and temperature gradients, transients and the appearance of discontinuities, factors which require mesh refinement during computations. Techniques are developed for temporal and spatial discretization of a model problem. Several observations are made regarding the explicit and implicit features of the calculations, the use of the Lax-Wendroff scheme to produce a mass-matrix for obtaining accurate results for transients, methods of performing stability analyses, and simplification techniques. Examples are provided of solving the nonlinear shallow-water equations and describing compressible flows, particularly transonic flows. Domain splitting is defined for improving the calculations at each time step and in different parts of the flow regime while simultaneously advancing the calculations towards a solution.
It is well-known that stresses (and strains) calculated by a displacement-based finite element analysis are generally not as accurate as the displacements. In addition, the calculated stress field is typically discontinuous at element interfaces. Because the stresses are typically of more interest than the displacements, several procedures have been proposed to obtain a smooth stress field, given the finite element stresses, and to improve the accuracy. Hinton and Irons introduced global least squares smoothing of discrete data defined on a plane using a finite element formulation. Tessler and co-workers recently developed a conceptually similar formulation for smoothing of two-dimensional data based on a discrete least square approximation with a penalty constraint. The penalty constraint results in a stress field which is C(exp 1)-continuous, a result not previously obtained. The approach requires additional, 'smoothing' finite element analysis and for their two-dimensional application, they used a conforming C(exp 0)-continuous triangular finite element based on a conforming plate element. This paper presents the results of a detailed investigation into the application of Tessler's smoothing procedure to the smoothing of finite element stresses from one-dimensional problems. Although the one-dimensional formulation has some practical applicability, such as in truss, beam, axisymmetric mechanics, and one-dimensional heat conduction, the primary motivation for developing the one-dimensional smoothing case is to explore the characteristics of the general smoothing strategy. In particular, it is used to describe the behavior of the method and to explore the suitability of criteria proposed for the smoothing analysis. Prior to presenting numerical results, the variational formulation of the smoothing strategy is presented and a criterion for the smoothing analysis is described.
The numerical solution of partial differential equations about irregular geometries and with varying characteristic scales has created the need for coordinate systems and associated transformations which reflect both geometric and physical requirements. The process of finding coordinate transformations in discrete representations is called 'grid generation'. The present investigation is concerned with three algebraic grid generation methods. The methods include transfinite interpolation, the multisurface method, and the two-boundary technique. Interpolation formulas, in terms of homotopic mappings and constraints in terms of point positions and/or derivatives are the essential elements of the techniques. The methods are relatively simple to understand, they are explicit and do not require extensive computational effort, and they have a high degree of generality.
New shock-capturing finite difference approximations for solving two scalar conservation law nonlinear partial differential equations describing inviscid, isentropic, compressible flows of aerodynamics at transonic speeds are presented. A global linear stability theorem is applied to these schemes in order to derive a necessary and sufficient condition for the finite element method. A technique is proposed to render the described approximations total variation-stable by applying the flux limiters to the nonlinear terms of the difference equation dimension by dimension. An entropy theorem applying to the approximations is proved, and an implicit, forward Euler-type time discretization of the approximation is presented. Results of some numerical experiments using the approximations are reported.
Of the two finite-element models presently used to investigate the effect of the conservation of integral invariants, by means of finite-element discretization schemes of the shallow-water equations, in order to serve as a paradigm of long-term atmospheric-model integrations, the first employs rectangular elements and conserves total energy, while the second uses triangular elements and a high-accuracy two-stage Numerov-Galerkin method. Attention is given to critical times for numerical nonlinear instability, as well as to the determination of the critical degree of dissipation entailed by the achievement of stable long-term integrations. Relative computational efficiency and accuracy comparisons are presented for the two finite-element schemes.
A new method has been developed to calculate the steady flow and structural deformations for fluid/structure interaction problems. The discretized fluid dynamic and structural equations are regarded as a single set of coupled, nonlinear, algebraic equations. The equilibrium solution is directly obtained using Newton's method. The governing equations used for the fluid flow are the two-dimensional Navier-Stokes equations, and a finite-element model is used to represent the structure. This paper describes the analytical method and presents sample calculations demonstrating the technique. The results show rapid convergence and good agreement with experimental data.
We consider the discontinuous Galerkin (DG) finite element discretization of first order systems of conservation laws derivable as moments of the kinetic Boltzmann equation. This includes well known conservation law systems such as the Euler For the class of first order nonlinear conservation laws equipped with an entropy extension, an energy analysis of the DG method for the Cauchy initial value problem is developed. Using this DG energy analysis, several new variants of existing numerical flux functions are derived and shown to be energy stable.
The silicon substrate carrier was created so that a large-area array (in this case 62,000+ elements of a microshutter array) and a variety of discrete passive and active devices could be mounted on a single board, similar to a printed circuit board. However, the density and number of interconnects far exceeds the capabilities of printed circuit board technology. To overcome this hurdle, a method was developed to fabricate this carrier out of silicon and implement silicon integrated circuit (IC) technology. This method achieves a large number of high-density metal interconnects; a 100-percent yield over a 6-in. (approximately equal to 15-cm) diameter wafer (one unit per wafer); a rigid, thermally compatible structure (all components and operating conditions) to cryogenic temperatures; re-workability and component replaceability, if required; and the ability to precisely cut large-area holes through the substrate. A method that would employ indium bump technology along with wafer-scale integration onto a silicon carrier was also developed. By establishing a silicon-based version of a printed circuit board, the objectives could be met with one solution. The silicon substrate would be 2 mm thick to survive the environmental loads of a launch. More than 2,300 metal traces and over 1,500 individual wire bonds are required. To mate the microshutter array to the silicon substrate, more than 10,000 indium bumps are required. A window was cut in the substrate to allow the light signal to pass through the substrate and reach the microshutter array. The substrate was also the receptacle for multiple unpackaged IC die wire-bonded directly to the substrate (thus conserving space over conventionally packaged die). Unique features of this technology include the implementation of a 2-mmthick silicon wafer to withstand extreme mechanical loads (from a rocket launch); integrated polysilicon resistor heaters directly on the substrate; the precise formation of an open aperture (approximately equal to 3x3cm) without any crack propagation; implementation of IR transmission blocking techniques; and compatibility with indium bump bonding. Although designed for the microshutter arrays for the NIRSpec instrument on the James Webb Space Telescope, these substrates can be linked to microshutter applications in the photomask generation and stepper equipment used to make ICs and microelectromechanical system (MEMS) devices.