Search NASA⌕ Search

SEARCH · Search NASA

Results for “nodal formulation”

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 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS2018

VCCT with Progressive Nodal Release for Simulating Mixed-Mode Delamination: Formulation, Algorithmic Improvements and Implications

To simulate the propagation of cracks using the Virtual Crack Closure Technique (VCCT), intermediate crack positions between existing node pairs can be modeled by progressively releasing the nodes using fracture mechanics principles. This approach has been implemented in commercial software and can be applied to simulate delamination growth. Recent research has provided insight into the synchronized mixed mode nodal release, extending and formalizing the methodology for delamination growth. This paper documents the development of the Jacobians associated with delamination propagation under quasi-static mixed-mode conditions, leading to improvements in convergence of the implicit nonlinear solvers used.

Mabson, G. E.↗
NASA NTRS2018

VCCT with Progressive Nodal Release for Simulating Mixed-Mode Delamination: Formulation, Algorithmic Improvements and Implications

To simulate the propagation of cracks using the Virtual Crack Closure Technique (VCCT), intermediate crack positions between existing node pairs can be modeled by progressively releasing the nodes using fracture mechanics principles. This approach has been implemented in commercial software and can be applied to simulate delamination growth. Recent research has provided insight into the synchronized mixed mode nodal release, extending and formalizing the methodology for delamination growth. This paper documents the development of the Jacobians associated with delamination propagation under quasi-static mixed-mode conditions, leading to improvements in convergence of the implicit nonlinear solvers used.

G E Mabson↗
NASA NTRS1992

Single-phase power distribution system power flow and fault analysis

Alternative methods for power flow and fault analysis of single-phase distribution systems are presented. The algorithms for both power flow and fault analysis utilize a generalized approach to network modeling. The generalized admittance matrix, formed using elements of linear graph theory, is an accurate network model for all possible single-phase network configurations. Unlike the standard nodal admittance matrix formulation algorithms, the generalized approach uses generalized component models for the transmission line and transformer. The standard assumption of a common node voltage reference point is not required to construct the generalized admittance matrix. Therefore, truly accurate simulation results can be obtained for networks that cannot be modeled using traditional techniques.

Halpin, S. M.↗
NASA NTRS1976

Performance of a recoverable tug for planetary missions including use of perigee propulsion and corrections for nodal regression

Mathematical equations are derived by using the Maximum Principle to obtain the maximum payload capability of a reusable tug for planetary missions. The mathematical formulation includes correction for nodal precession of the space shuttle orbit. The tug performs this nodal correction in returning to this precessed orbit. The sample case analyzed represents an inner planet mission as defined by the declination (fixed) and right ascension of the outgoing asymptote and the mission energy. Payload capability is derived for a typical cryogenic tug and the sample case with and without perigee propulsion. Optimal trajectory profiles and some important orbital elements are also discussed.

Borsody, J.↗
NASA NTRS1976

Performance of a recoverable Tug for planetary missions including use of perigee propulsion and corrections for nodal regression

Mathematical equations are derived by using the Maximum Principle to obtain the maximum payload capability of a reusable Tug for planetary missions. The mathematical formulation includes correction for nodal precession of the Space Shuttle orbit. The Tug performs this nodal correction in returning to this precessed orbit. The sample case analyzed represents an inner planet mission as defined by the declination (fixed) and right ascension of the outgoing asymptote and the mission energy. Payload capability is derived for a typical cryogenic Tug and the sample case with and without perigee propulsion. Optimal trajectory profiles and some important orbital elements are also discussed.

Borsody, J.↗
NASA NTRS1974

Fluid-structure finite-element vibrational analysis

A fluid finite element has been developed for a quasi-compressible fluid. Both kinetic and potential energy are expressed as functions of nodal displacements. Thus, the formulation is similar to that used for structural elements, with the only differences being that the fluid can possess gravitational potential, and the constitutive equations for fluid contain no shear coefficients. Using this approach, structural and fluid elements can be used interchangeably in existing efficient sparse-matrix structural computer programs such as SPAR. The theoretical development of the element formulations and the relationships of the local and global coordinates are shown. Solutions of fluid slosh, liquid compressibility, and coupled fluid-shell oscillation problems which were completed using a temporary digital computer program are shown. The frequency correlation of the solutions with classical theory is excellent.

Feng, G. C.↗
NASA NTRS1991

The MHOST finite element program: 3-D inelastic analysis methods for hot section components. Volume 1: Theoretical manual

Formulations and algorithms implemented in the MHOST finite element program are discussed. The code uses a novel concept of the mixed iterative solution technique for the efficient 3-D computations of turbine engine hot section components. The general framework of variational formulation and solution algorithms are discussed which were derived from the mixed three field Hu-Washizu principle. This formulation enables the use of nodal interpolation for coordinates, displacements, strains, and stresses. Algorithmic description of the mixed iterative method includes variations for the quasi static, transient dynamic and buckling analyses. The global-local analysis procedure referred to as the subelement refinement is developed in the framework of the mixed iterative solution, of which the detail is presented. The numerically integrated isoparametric elements implemented in the framework is discussed. Methods to filter certain parts of strain and project the element discontinuous quantities to the nodes are developed for a family of linear elements. Integration algorithms are described for linear and nonlinear equations included in MHOST program.

Nakazawa, Shohei↗
NASA NTRS1991

Computational methods for shock waves in three-dimensional supersonic flow

A central artificial-viscosity and an upwind-biased difference method are contrived to solve the Euler equations for flowfields over typical spacecrafts. The spatial discretization is based on either nodal or cell-vertex formulation in the domain extending from free stream to the end of the vehicle. The outer boundary is treated as a bow shock in the first method but is placed in the free stream in the second, which captures both bow and internal shocks using an approximate Riemann solver based on high-order extrapolation to the cell face. These methods were tested for the Shuttle and Hermes orbiters at wind-tunnel conditions and angles of attack ranging from 0 to 60 deg. The artificial-viscosity method incorporated with a shock-fitting procedure shows smeared crossflow and wing-shock positions and required 15 percent more CPU per node than the upwind method. Greater flexibility and robustness is demonstrated by the latter on a fixed grid for all cases considered.

Li, Chien-Peng↗
NASA NTRS1996

A viscoelastic higher-order beam finite element

A viscoelastic internal variable constitutive theory is applied to a higher-order elastic beam theory and finite element formulation. The behavior of the viscous material in the beam is approximately modeled as a Maxwell solid. The finite element formulation requires additional sets of nodal variables for each relaxation time constant needed by the Maxwell solid. Recent developments in modeling viscoelastic material behavior with strain variables that are conjugate to the elastic strain measures are combined with advances in modeling through-the-thickness stresses and strains in thick beams. The result is a viscous thick-beam finite element that possesses superior characteristics for transient analysis since its nodal viscous forces are not linearly dependent an the nodal velocities, which is the case when damping matrices are used. Instead, the nodal viscous forces are directly dependent on the material's relaxation spectrum and the history of the nodal variables through a differential form of the constitutive law for a Maxwell solid. The thick beam quasistatic analysis is explored herein as a first step towards developing more complex viscoelastic models for thick plates and shells, and for dynamic analyses. The internal variable constitutive theory is derived directly from the Boltzmann superposition theorem. The mechanical strains and the conjugate internal strains are shown to be related through a system of first-order, ordinary differential equations. The total time-dependent stress is the superposition of its elastic and viscous components. Equations of motion for the solid are derived from the virtual work principle using the total time-dependent stress. Numerical examples for the problems of relaxation, creep, and cyclic creep are carried out for a beam made from an orthotropic Maxwell solid.

Johnson, Arthur R.↗
NASA NTRS1975

Long period nodal motion of sun synchronous orbits

An approximative model is formulated for assessing these perturbations that significantly affect long term modal motion of sun synchronous orbits. Computer simulations with several independent computer programs consider zonal and tesseral gravitational harmonics, third body gravitational disturbances induced by the sun and the moon, and atmospheric drag. A pendulum model consisting of evenzonal harmonics through order 4 and solar gravity dominated nodal motion approximation. This pendulum motion results from solar gravity inducing an inclination oscillation which couples into the nodal precession induced by the earth's oblateness. The pendulum model correlated well with simulations observed flight data.

Duck, K. I.↗
NASA NTRS1986

Senstitivty analysis and optimization of nodal point placement for vibration reduction

A method is developed for sensitivity analysis and optimization of nodal point locations in connection with vibration reduction. A straightforward derivation of the expression for the derivative of nodal locations is given, and the role of the derivative in assessing design trends is demonstrated. An optimization process is developed which uses added lumped masses on the structure as design variables to move the node to a preselected location - for example, where low response amplitude is required or to a point which makes the mode shape nearly orthogonal to the force distribution, thereby minimizing the generalized force. The optimization formulation leads to values for added masses that adjust a nodal location while minimizing the total amount of added mass required to do so. As an example, the node of the second mode of a cantilever box beam is relocated to coincide with the centroid of a prescribed force distribution, thereby reducing the generalized force substantially without adding excessive mass. A comparison with an optimization formulation that directly minimizes the generalized force indicates that nodal placement gives essentially a minimum generalized force when the node is appropriately placed.

Pritchard, J. I.↗
NASA NTRS1987

Sensitivity derivatives and optimization of nodal point locations for vibration reduction

A method is developed for sensitivity analysis and optimization of nodal point locations in connection with vibration reduction. A straightforward derivation of the expression for the derivative of nodal locations is given, and the role of the derivative in assessing design trends is demonstrated. An optimization process is developed which uses added lumped masses on the structure as design variables to move the node to a preselected location; for example, where low response amplitude is required or to a point which makes the mode shape nearly orthogonal to the force distribution, thereby minimizing the generalized force. The optimization formulation leads to values for added masses that adjust a nodal location while minimizing the total amount of added mass required to do so. As an example, the node of the second mode of a cantilever box beam is relocated to coincide with the centroid of a prescribed force distribution, thereby reducing the generalized force substantially without adding excessive mass. A comparison with an optimization formulation that directly minimizes the generalized force indicates that nodal placement gives essentially a minimum generalized force when the node is appropriately placed.

Pritchard, Jocelyn I.↗
NASA NTRS1987

Sensitivity analysis and optimization of nodal point placement for vibration reduction

A method is developed for sensitivity analysis and optimization of nodal point locations in connection with vibration reduction. A straightforward derivation of the expression for the derivative of nodal locations is given, and the role of the derivative in assessing design trends is demonstrated. An optimization process is developed which uses added lumped masses on the structure as design variables to move the node to a preselected location - for example, where low response amplitude is required or to a point which makes the mode shape nearly orthogonal to the force distribution, thereby minimizing the generalized force. The optimization formulation leads to values for added masses that adjust a nodal location while minimizing the total amount of added mass required to do so. As an example, the node of the second mode of a cantilever box beam is relocated to coincide with the centroid of a prescribed force distribution, thereby reducing the generalized force substantially without adding excessive mass. A comparison with an optimization formulation that directly minimizes the generalized force indicates that nodal placement gives essentially a minimum generalized force when the node is appropriately placed.

Pritchard, J. I.↗
NASA NTRS1998

WARP3D-Release 10.8: Dynamic Nonlinear Analysis of Solids using a Preconditioned Conjugate Gradient Software Architecture

This report describes theoretical background material and commands necessary to use the WARP3D finite element code. WARP3D is under continuing development as a research code for the solution of very large-scale, 3-D solid models subjected to static and dynamic loads. Specific features in the code oriented toward the investigation of ductile fracture in metals include a robust finite strain formulation, a general J-integral computation facility (with inertia, face loading), an element extinction facility to model crack growth, nonlinear material models including viscoplastic effects, and the Gurson-Tver-gaard dilatant plasticity model for void growth. The nonlinear, dynamic equilibrium equations are solved using an incremental-iterative, implicit formulation with full Newton iterations to eliminate residual nodal forces. The history integration of the nonlinear equations of motion is accomplished with Newmarks Beta method. A central feature of WARP3D involves the use of a linear-preconditioned conjugate gradient (LPCG) solver implemented in an element-by-element format to replace a conventional direct linear equation solver. This software architecture dramatically reduces both the memory requirements and CPU time for very large, nonlinear solid models since formation of the assembled (dynamic) stiffness matrix is avoided. Analyses thus exhibit the numerical stability for large time (load) steps provided by the implicit formulation coupled with the low memory requirements characteristic of an explicit code. In addition to the much lower memory requirements of the LPCG solver, the CPU time required for solution of the linear equations during each Newton iteration is generally one-half or less of the CPU time required for a traditional direct solver. All other computational aspects of the code (element stiffnesses, element strains, stress updating, element internal forces) are implemented in the element-by- element, blocked architecture. This greatly improves vectorization of the code on uni-processor hardware and enables straightforward parallel-vector processing of element blocks on multi-processor hardware.

Koppenhoefer, Kyle C.↗
NASA NTRS1999

Unified CFD Methods Via Flowfield-Dependent Variation Theory

This paper addresses the flowfield-dependent variation (FDV) methods in which complex physical phenomena are taken into account in the final form of partial differential equations to be solved so that finite difference methods (FDM) or finite element methods (FEM) themselves will not dictate the physics, but rather are no more than simply the options how to discretize between adjacent nodal points or within an element. The variation parameters introduced in the formulation are calculated from the current flowfield based on changes of Mach numbers, Reynolds numbers, Peclet numbers, and Damkohler numbers between adjacent nodal points, which play many significant roles such as adjusting the governing equations (hyperbolic, parabolic, and/or e!liptic), resolving various physical phenomena, and controlling the accuracy and stability of the numerical solution. The theory is verified by a number of example problems addressing the physical implications of the variation parameters which resemble the flowfield itself, shock capturing mechanism, transitions and interactions between inviscid/viscous, compressibility/incompressibility, and laminar/turbulent flows.

Chung, T. J.↗
NASA NTRS1989

Specialty functions singularity mechanics problems

The focus is in the development of more accurate and efficient advanced methods for solution of singular problems encountered in mechanics. At present, finite element methods in conjunction with special functions, boolean sum and blending interpolations are being considered. In dealing with systems which contain a singularity, special finite elements are being formulated to be used in singular regions. Further, special transition elements are being formulated to couple the special element to the mesh that models the rest of the system, and to be used in conjunction with 1-D, 2-D and 3-D elements within the same mesh. Computational simulation with a least squares fit is being utilized to construct special elements, if there is an unknown singularity in the system. A novel approach is taken in formulation of the elements in that: (1) the material properties are modified to include time, temperature, coordinate and stress dependant behavior within the element; (2) material properties vary at nodal points of the elements; (3) a hidden-symbolic computation scheme is developed and utilized in formulating the elements; and (4) special functions and boolean sum are utilized in order to interpolate the field variables and their derivatives along the boundary of the elements. It may be noted that the proposed methods are also applicable to fluids and coupled problems.

Sarigul, Nesrin↗
NASA NTRS2007

Flow Separation Side Loads Excitation of Rocket Nozzle FEM

Modern rocket nozzles are designed to operate over a wide range of altitudes, and are also built with large aspect ratios to enable high efficiencies. Nozzles designed to operate over specific regions of a trajectory are being replaced in modern launch vehicles by those that are designed to operate from earth to orbit. This is happening in parallel with modern manufacturing and wall cooling techniques allowing for larger aspect ratio nozzles to be produced. Such nozzles, though operating over a large range of altitudes and ambient pressures, are typically designed for one specific altitude. Above that altitude the nozzle flow is 'underexpanded' and below that altitude, the nozzle flow is 'overexpanded'. In both conditions the nozzle produces less than the maximum possible thrust at that altitude. Usually the nozzle design altitude is well above sea level, leaving the nozzle flow in an overexpanded state for its start up as well as for its ground testing where, if it is a reusable nozzle such as the Space Shuttle Main Engine (SSME), the nozzle will operate for the majority of its life. Overexpansion in a rocket nozzle presents the critical, and sometimes design driving, problem of flow separation induced side loads. To increase their understanding of nozzle side loads, engineers at MSFC began an investigation in 2000 into the phenomenon through a task entitled "Characterization and Accurate Modeling of Rocket Engine Nozzle Side Loads", led by A. Brown. The stated objective of this study was to develop a methodology to accurately predict the character and magnitude of nozzle side loads. The study included further hot-fire testing of the MC-l engine, cold flow testing of subscale nozzles, CFD analyses of both hot-fire and cold flow nozzle testing, and finite element (fe.) analysis of the MC-1 engine and cold flow tested nozzles. A follow on task included an effort to formulate a simplified methodology for modeling a side load during a two nodal diameter fluid/structure interaction for a single moment in time.

Smalley, Kurt B.↗
NASA NTRS1984

On the development of hierarchical solution strategies for nonlinear finite element formulations

This paper develops a hierarchical type solution scheme which can handle the field equations associated with nonlinear finite element simulations. The overall procedure possesses various levels of application namely degree of freedom, nodal, elemental, substructural as well as global. In particular iteration, updating, assembly and solution control occurs at the various hierarchical levels. Due to the manner of formulation, the degree of matrix inversion depends on the size of the various hierarchical partitioned groups. In this context, degree of freedom partitioning requires no inversion. To benchmark the overall scheme, the results of several numerical examples are presented.

Padovan, J.↗